| 查看: 764 | 回复: 2 | |||
[交流]
【求助】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 |
回复此楼» 猜你喜欢
26年申博自荐-计算机视觉
已经有4人回复
-
考博
已经有4人回复
-
药化及相关博士的申请
已经有3人回复
-
一篇MDPI论文改变了学习工作和生活
已经有4人回复
-
一个化合物的合成路线：CAS:367929-02-0 名称：8β-乙烯基雌二醇
已经有4人回复
-
太白金星有点烦
已经有3人回复
-
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业
已经有3人回复
-
收到国自然专家邀请后几年才会有本子送过来评
已经有3人回复
-
有没有快的中文核心比较快录用的，纳米材料光催化
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何在lammps里加入新的势函数形式？
已经有6人回复
- 各位高手，怎么在windows下编译LAMMPS？
已经有4人回复
- lammps的时间步长问题
已经有16人回复
- 薄膜在lammps条件下的谐振
已经有5人回复
- LAMMPS Error:lost atoms
已经有7人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- 【求助】LAMMPS、NAMD和CHARMM谁最适合生物大分子模拟？
已经有8人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
2026年江西师范大学人文地理与低碳经济方向拟招收1名“申请入学制”学术博士研究生
+2/162
-
华侨大学全日制硕士博士（可不需要统考过了四级就行）-厦门校区-明天截止-急招
+1/129
-
西澳大学 环境工程 全奖博士 2027 CSC博士招生（接收联培培养博士生）
+2/118
-
环境工程专业实习生|纯科研
+1/80
-
精细结构常数α=1/137的素数起源：论文投稿审核历程（A202604-1381）
+1/75
-
安农大李亚辉教授招收化学、生物、农药学、植保类方向申请考核制博士生
+3/59
-
2026年江西师范大学药学院陈芬儿院士课题组招收智能药学博士生
+1/32
-
西北工业大学物理科学与技术学院-诚聘【AI与光学交叉方向】博士后
+1/31
-
山东大学齐鲁交通学院博士招生--压缩空气储能研究方向
+2/30
-
智合健物课题组2026年博士生招生(5月6日—2026年5月13日)
+1/28
-
树脂真空干燥除水时的泡沫问题
+2/23
-
深圳理工大学-湖南大学项目招收钙钛矿方向博士生
+1/11
-
杨亚研究员课题组（广西大学-纳米能源研究中心）招收2026年物理学专业学术型博士
+1/5
-
招收2026年秋季入学博士生1名（河北工业大学/北京科技大学联合 增材制造/生物材料）
+1/5
-
双一流高校南林化学工程学院-国家级青年人才团队招2026级博士（最后一波）
+1/5
-
中山大学材料科学与工程学院周剑课题组2026年9月入学博士生招生需求
+1/3
-
上海交通大学陈向洋计算化学课题组诚聘2027级博士研究生
+1/3
-
中国药科大学陈达威教授2026招收药理学学术型博士生一名(截止日期5月12日)
+1/2
-
华南理工大学材料学院招收2026年博士生（有机、高分子合成、二氧化碳转化）
+1/2
-
综述投稿，we两个月
+1/1
2楼2010-11-26 12:54:25
★ ★ ★ ★
cxz314(金币+5):非常感谢！我在核实一下 2010-11-27 19:57:37
ghcacj(金币+4):谢谢 2010-11-29 18:23:05
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, 都是约化单位么？建议你再核实一下。不过我不是很了解这个势，也可能没有问题。 |
3楼2010-11-27 17:39:21