| 查看: 602 | 回复: 0 | |||
wusongz金虫 (小有名气)
|
[求助]
lammps得到如图g(r)，为什么峰会劈裂？
|
本人初学者，实在搞不懂。 units metal atom_style charge boundary p p p variable x index 1250 print "--------Temperature=$x K---------" read_data data_amorphous_vasp.CU pair_style comb pair_coeff * * ffield.comb Cu Si neighbor 0.5 bin neigh_modify every 4 delay 0 check yes group CU id <> 1 3200 group SI id <> 3201 4000 compute A CU property/atom q compute Q CU reduce ave c_A compute B SI property/atom q compute W SI reduce ave c_B set group CU charge 0.0 set group SI charge 0.0 timestep 0.0001 #dump 1 all cfg 10 *.cfg id type xs ys zs id q vx vy vz fx fy fz #dump_modify 1 element Cu Si dump 2 all atom 100 CuSi_$x.lammpstrj thermo 1000 velocity all create $x 825577 thermo_style custom step temp epair etotal vol press pxx pyy pzz c_Q c_W min_modify dmax 0.05 line quadratic minimize 1.0e-14 1.0e-20 20 100 fix 1 all nvt temp $x $x 0.01 drag 0.2 fix 2 all qeq/comb 1 3e-3 run 10000 unfix 1 min_modify dmax 0.05 line quadratic minimize 1.0e-14 1.0e-20 100000 100000 fix 3 all box/relax x 0.0 y 0.0 z 0.0 vmax 0.001 fix 4 all npt temp $x $x 0.01 x 0.0 0.0 0.01 y 0.0 0.0 0.01 z 0.0 0.0 0.01 drag 0.2 fix 5 all qeq/comb 1 3e-3 run 2200000 clear next x jump in.comb.2 |
回复此楼» 猜你喜欢
存款400万可以在学校里躺平吗
已经有36人回复
-
英文综述是否需要润色及查重
已经有5人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
最失望的一年
已经有18人回复
找到一些相关的精华帖子，希望有用哦~
薄膜在lammps条件下的谐振
已经有5人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- lammps中固定键长
已经有4人回复
- 【求助】lammps如何测量模拟以后材料的晶格参数？
已经有7人回复
- 【讨论】lammps中关于fix ave/spatial命令的输出问题
已经有4人回复
- 【求助】lammps中多元化合物的晶格创建问题
已经有14人回复
- 【讨论】lammps运行eam/fs出错
已经有38人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
- 【求助】LAMMPS中非周期性边界条件和Born-Mayer-Huggins势函数不会设置
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人