| 查看: 3170 | 回复: 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 |
回复此楼» 猜你喜欢
光学工程学硕调剂信息
已经有32人回复
-
欢迎加入课题组
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有178人回复
-
散金币，求好运，祝面上顺利！
已经有79人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有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
