| 查看: 3167 | 回复: 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 |
回复此楼» 猜你喜欢
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有243人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有1人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有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