| 查看: 3159 | 回复: 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 |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有229人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有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