| 查看: 649 | 回复: 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 |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
找到一些相关的精华帖子，希望有用哦~
薄膜在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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
