| 3 | 1/1 | 返回列表 |
| 查看: 3171 | 回復(fù): 2 | ||
小卒無(wú)名金蟲(chóng) (小有名氣)
|
[求助]
Lammps 中NEB計(jì)算
|
|
根據(jù)http://simulation.haotui.com/viewthread.php?tid=40495中的例子,編寫(xiě)的用于尋找過(guò)渡態(tài)的程序,程序跑不通,顯示:cann't use NEB with a single replica,對(duì)于官網(wǎng)手冊(cè)也是一頭霧水。問(wèn)題:1,不知道怎么指示初態(tài)和終態(tài);2,如何保存中間過(guò)程中的原子坐標(biāo)及能量。 atom_modify map array sort 0 0.0 # ---------- Create Atoms --------------------- read_data readdata.1 variable u uloop 2 # ---------- Define Interatomic Potential --------------------- pair_style eam/alloy pair_coeff * * Cu01.eam.alloy Cu neighbor 2.0 bin neigh_modify delay 10 check yes # ---------- Define Settings --------------------- compute perpot all pe/atom compute perkin all ke/atom compute pot all pe compute kin all ke #---------------Fix nvt----------------# #velocity all set 0.0 0.0 0.0 mom yes rot no fix 1 all setforce 0 0 0 fix 2 all nvt temp 0.0001 0.0001 10.0 drag 0.2 #fix 3 all neb 1.0 minimize 1e-25 1e-25 5000 10000 min_style quickmin #---------------Output Thermo----------------# reset_timestep 0 timestep 0.01 thermo 10 thermo_style custom step lx ly lz press pxx pyy pzz vol temp c_pot c_kin etotal dump 1 all custom 100 dump.neb.$u id type x y z neb 0.1 0.01 1000 500 1 coords.final next u jump in print "All done" fianlfile #LAMMPS Description 3980 913 1.74218 2.37651 1.47087 909 1.67493 4.3525 3.19397 3120 1.91223 0.409578 3.2502 915 3.45571 4.08904 1.32395 3125 3.73239 0.639876 1.40948 3131 3.75153 2.3175 3.21221 908 -0.050082 2.2276 3.3366 3130 5.55566 2.38966 1.39588 3127 5.16336 0.426954 3.35545 917 5.35775 4.12671 3.19239 3070 7.09337 0.47682 1.38471 3072 7.20937 2.33175 3.10631 845 7.06057 4.28264 1.32245 3068 8.79948 2.44301 1.33783 3064 8.8503 0.5081 3.04721 3063 10.771 0.613592 1.34068 848 8.89819 4.17629 3.06949 846 10.807 4.28234 1.37612 854 10.6779 2.2451 3.18697 3079 12.5435 2.4197 1.3716 |
金蟲(chóng) (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-29 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿:西北大學(xué),英一數(shù)一408-284分求調(diào)劑 +3 | 12.27 2026-03-27 | 3/150 |
|
|
[考研] 343求調(diào)劑 +6 | 愛(ài)羈絆 2026-03-29 | 6/300 |
|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +3 | Nadiums 2026-03-25 | 4/200 |
|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 085600,材料與化工321分求調(diào)劑 +9 | 大饞小子 2026-03-28 | 9/450 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +5 | niko- 2026-03-27 | 5/250 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見(jiàn)α 2026-03-25 | 3/150 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專(zhuān)碩324分,求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|