24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cxz314

木虫 (正式写手)

应助: 1 (幼儿园)
贵宾: 0.2
金币: 3815.1
帖子: 997
在线: 189小时
虫号: 475745

[交流] 【求助】Gaussian core model 的LAMMPS模拟

大家好
我用LAMMPS模拟Gaussian core model ，发现模拟后粒子都集中在一个非常小的区域，都重叠在一起了。不知道是为什么。请各位指点一下
下面是我的LAMMPS模拟的in 文件
#GCM AT 0.04
units  lj
boundary  p p p
atom_style  atomic

lattice fcc 3.419952
region box block 0 8.0 0 8.0 0 8.0
create_box  1 box
create_atoms  1 box

pair_style  gauss  3.2
pair_coeff  * * 1.0 0.9
mass       1 1.0

timestep  0.05
thermo  100

neighbor  0.5 bin
neigh_modify  every 4 delay 0 check yes

velocity  all create 0.04 4928459 rot yes dist gaussian
fix       1 all nvt temp 0.04 0.04 1.5 drag 0.3
thermo_style custom step temp press vol etotal
log       log_1.gcm
run    500000

dump  1 all custom 100 dump_1.gcm id x y z
#fix    3 all msd 100 msd_1.gcm
run    10000
unfix    1
clear
next    x
jump    in.eam

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

如何在lammps里加入新的势函数形式？已经有6人回复
各位高手，怎么在windows下编译LAMMPS？已经有4人回复
lammps的时间步长问题已经有16人回复
薄膜在lammps条件下的谐振已经有5人回复
LAMMPS Error:lost atoms 已经有7人回复
寻求lammps的应用欢迎大家畅所欲言已经有8人回复
【求助】LAMMPS、NAMD和CHARMM谁最适合生物大分子模拟？已经有8人回复
【求助】Lammps 多机并行的问题已经有11人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-11-26 09:01:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★ ★ ★
cxz314(金币+5):非常感谢！我在核实一下 2010-11-27 19:57:37
ghcacj(金币+4):谢谢 2010-11-29 18:23:05

cxz314(金币+1):谢谢！是这样的Gaussian core model 的约化单位和LJ势的约化单位一样的

这个。。。我有点不明白，pair_style gauss的势函数形式是 E=Aexp(-Br^2) ，这样的势，要用约化单位么？那你的sigma，epsilon之类是如何设置的？都是1么？那你的 lattice fcc 3.419952，还有region, 都是约化单位么？建议你再核实一下。不过我不是很了解这个势，也可能没有问题。