| 查看: 3226 | 回复: 2 | ||
[求助]
Lammps 中NEB计算
|
|
根据http://simulation.haotui.com/viewthread.php?tid=40495中的例子，编写的用于寻找过渡态的程序，程序跑不通，显示：cann't use NEB with a single replica,对于官网手册也是一头雾水。问题：1，不知道怎么指示初态和终态；2，如何保存中间过程中的原子坐标及能量。 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 |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有269人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- “慢”扩散的模拟计算？？？
已经有8人回复
- Lammps计算溶液介电常数（dielectric constant）
已经有17人回复
- Lammps计算油水界面，油相密度、平衡等问题求指导
已经有6人回复
- 用Lammps模拟计算水的密度
已经有10人回复
- 求助如何计算速度自相关函数。
已经有12人回复
- 怎样生成movie文件
已经有6人回复
- 为什么同样的设置，计算出的结果有时会相差很大呢？
已经有8人回复
- 由照片中的距离计算出实际距离
已经有13人回复
- 求助LAMMPS关于MSD计算的问题
已经有17人回复
- lammps计算中出现原子距离太近现象
已经有11人回复
- LAMMPS计算无反应，nan
已经有4人回复
- LAMMPS并行计算的问题（cpu——time关系）
已经有17人回复
- linux下的计算精度的设置
已经有5人回复
- 用MS建模导入LAMMPS计算，是不是不能用compass力场？
已经有9人回复
- lammps求助：计算含有不同原子体系的热导率的in文件该怎么修改
已经有7人回复
- lammps的多尺度计算
已经有15人回复
- 【讨论】LAMMPS 晶界计算 输入脚本
已经有4人回复
- 【讨论】lammps 用sw势计算硅的熔点怎么这么高（2400k),附in文件
已经有13人回复
2楼2015-07-08 13:11:55
3楼2015-09-10 17:10:17
50