 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1430 | 回復(fù): 0 | ||
| 【懸賞金幣】回答本帖問題，作者懶羊羊要當(dāng)大廚將贈(zèng)送您 5 個(gè)金幣 | ||
懶羊羊要當(dāng)大廚新蟲 (初入文壇)
|
[求助]
LAMMPS計(jì)算外加電場(chǎng)作用下reaxff
|
|
|
各位老師好，我在計(jì)算微波作用下焦炭CO2氣化反應(yīng)時(shí)，根據(jù)參考文獻(xiàn)通過外加電場(chǎng)來模擬微波的作用。令人沮喪的是，在外加電場(chǎng)作用下，反應(yīng)不僅沒有得到促進(jìn)，反而呈現(xiàn)抑制作用，這與文獻(xiàn)計(jì)算結(jié)果和實(shí)驗(yàn)結(jié)果相反。文獻(xiàn)使用的計(jì)算軟件是ADS，在偶然的機(jī)會(huì)下，我發(fā)現(xiàn)在在ADS官方網(wǎng)站上有這樣一句話“LAMMPS計(jì)算外加電場(chǎng)作用下反應(yīng)力場(chǎng)計(jì)算，僅考慮了電場(chǎng)的力的作用，而沒有考慮電場(chǎng)能量的作用（https://v.fermitech.com.cn/wiki/doku.php?id=adf%3Aapplywallpotential&sid_for_share=99125_4）”。這確實(shí)有可能造成我的計(jì)算結(jié)果與實(shí)驗(yàn)和文獻(xiàn)相反，希望各位老師能為我指點(diǎn)。下面是我計(jì)算所使用的in文件 units real dimension 3 boundary p p f neighbor 2 bin neigh_modify every 10 delay 0 check yes timestep 0.1 atom_style full pair_style reaxff NULL safezone 5.0 mincap 600 minhbonds 1000 variable T equal 3000.0 variable r_p string "result/0.04" variable Ez equal 0.04 log ${r_p}/logout.lammps read_data H2O_c.data fix zwalls all wall/reflect zlo EDGE zhi EDGE pair_coeff * * ffield.reax.CHON C N O H fix 1 all qeq/reaxff 1 0.0 10.0 1.0e-6 reaxff min_style sd minimize 1e-8 1e-8 10000 10000 velocity all create 300.0 4928459 rot yes dist gaussian reset_timestep 0 thermo 200 thermo_style custom step temp pe ke v_Ez thermo_modify flush yes fix 2 all nvt temp 300.0 300.0 10 fix 3 all reaxff/species 1 1 200 ${r_p}/NVT.out element C N O H cutoff 3 3 0.5 fix 4 all reaxff/bonds 200 ${r_p}/bonds_nvt.txt dump 1 all xyz 200 ${r_p}/NVT.lammpstrj dump_modify 1 element C N O H run 200000 unfix 2 unfix 3 unfix 4 undump 1 reset_timestep 0 fix 5 all nvt temp 300 ${T} 10.0 fix 6 all reaxff/species 1 1 200 ${r_p}/heat.out element C N O H cutoff 3 3 0.5 fix 7 all reaxff/bonds 200 ${r_p}/bonds_heat.txt run 200000 unfix 5 unfix 6 unfix 7 reset_timestep 0 fix e all efield 0.0 0.0 v_Ez fix t all nve fix 8 all temp/berendsen ${T} ${T} 1000 dump 2 all xyz 200 ${r_p}/gas.lammpstrj dump_modify 2 element C N O H fix 9 all reaxff/species 1 1 200 ${r_p}/gas.out element C N O H cutoff 3 3 0.5 fix 10 all reaxff/bonds 200 ${r_p}/bonds_gas.txt run 5000000 |
回復(fù)此樓» 猜你喜歡
新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生
已經(jīng)有5人回復(fù)
-
材料專碩290求調(diào)劑
已經(jīng)有4人回復(fù)
-
不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來
已經(jīng)有5人回復(fù)
-
招調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿蘇大材料工程專碩293求調(diào)劑
已經(jīng)有3人回復(fù)
-
274求調(diào)劑
已經(jīng)有9人回復(fù)
-
085602高分子方向求調(diào)劑
已經(jīng)有3人回復(fù)
-
歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組
已經(jīng)有22人回復(fù)
-
289求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 懶羊羊要當(dāng)大廚 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招調(diào)劑 +6 | 帆船哥 2026-03-04 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 材料調(diào)劑 +5 | ws 上岸鴨 2026-03-05 | 5/250 |
|
|
|
[考研] 085602 293分求調(diào)劑 +3 | SivanNano. 2026-03-05 | 3/150 |
|
|
|
[考研] 材料調(diào)劑 +4 | L9370 2026-03-05 | 4/200 |
|
|
[考研] 材料與化工304求調(diào)劑 +5 | 邱gl 2026-03-05 | 8/400 |
|
|
|
[考研] 080500材料科學(xué)與工程 +11 | 202114020319 2026-03-03 | 11/550 |
|
|
|
[考研] 0856材料與化工求調(diào)劑！ +5 | 化工考生111 2026-03-04 | 9/450 |
|
|
|
[考研] 一志愿東北大學(xué)080500，英一數(shù)二343分求調(diào)劑 +5 | Seele_v5 2026-03-04 | 5/250 |
|
|
|
[考研] 0856材料與化工，270求調(diào)劑 +17 | YXCT 2026-03-01 | 20/1000 |
|
|
|
[考研] 085600材料與化工 298 調(diào)劑 +9 | 小西笑嘻嘻 2026-03-03 | 9/450 |
|
|
|
[考研] 293求調(diào)劑 +3 | 是樂渝哇 2026-03-04 | 3/150 |
|
|
|
[考研] 085600 材料與化工 298 +14 | 小西笑嘻嘻 2026-03-03 | 14/700 |
|
|
|
[考研] 武理材料工程302調(diào)劑 +3 | Doleres 2026-03-03 | 6/300 |
|
|
|
[考研] 0703 總分319求調(diào)劑 +4 | Xinyuu 2026-03-02 | 4/200 |
|
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
|
[考研] 324求調(diào)劑 +4 | wxz2 2026-03-03 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|
|
|
[考研] 材料調(diào)劑 +6 | 愛擦汗的可樂冰 2026-02-28 | 7/350 |
|
|
|
[考研] 295復(fù)試調(diào)劑 +3 | 簡木ChuFront 2026-03-01 | 3/150 |
|
|
|
[論文投稿]
求助coordination chemistry reviews 的寫作模板
 10+3
|
ljplijiapeng 2026-02-27 | 4/200 |
|