| 查看: 729 | 回复: 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 |
回复此楼» 猜你喜欢
交叉科学部支持青年基金，对三无青椒是个机会吗？
已经有7人回复
-
青椒八年已不青，大家都被折磨成啥样了？
已经有15人回复
-
多组分精馏求助
已经有4人回复
-
免疫学博士有名额，速联系
已经有4人回复
-
国家基金申请书模板内插入图片不可调整大小?
已经有6人回复
-
国家级人才课题组招收2026年入学博士
已经有5人回复
-
Fe3O4@SiO2合成
已经有6人回复
-
青年基金C终止
已经有4人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何在lammps里加入新的势函数形式？
已经有6人回复
- 各位高手，怎么在windows下编译LAMMPS？
已经有4人回复
- lammps的时间步长问题
已经有16人回复
- 薄膜在lammps条件下的谐振
已经有5人回复
- LAMMPS Error:lost atoms
已经有7人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- 【求助】LAMMPS、NAMD和CHARMM谁最适合生物大分子模拟？
已经有8人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
» 抢金币啦！回帖就可以得到:
-
丙烯液相
+1/81
-
诚聘 有机光电材料计算方向 博士后、科研助理、访问学生
+1/78
-
中国石油大学（北京）2026年博士招生2名：化学、材料、石油工程：油田化学
+1/73
-
西北工业大学民航学院招博士与硕士复合材料方向
+1/72
-
[求助] 半导体物理PPT课件 东华理工机电学院 彭新村老师
+1/71
-
征婚
+1/60
-
以色列理工-生物质塑料等催化转化及流体力学方向---全奖博士研究生和科研助理
+2/40
-
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生
+1/28
-
中国科学院上海光学精密机械研究所 特种强激光薄膜课题组
+1/25
-
青岛大学招收少数民族【少干计划】生物与医药博士研究生
+1/24
-
复旦大学化学系凡勇教授/张凡教授团队招聘博士后
+1/15
-
美国密苏里大学堪萨斯城分校（UMKC）生物材料诚聘全奖博士
+1/11
-
诚征女友，非诚勿扰
+3/11
-
浙江大学 “分子智造”课题组 诚聘 博士后及科研助理
+1/8
-
法国斯特拉斯堡大学有机光伏全奖博士招聘
+1/7
-
电子科技大学崔春华课题组招收物理化学背景博士生1名-申请考核制
+2/6
-
上海交通大学化学化工学院张智涛课题组诚聘博士后
+1/3
-
福建师范大学柔性电子学院 院士团队招2026级博士 光电器件、发光传感忆阻器
+1/2
-
澳科大诚招2026年秋季生物材料全奖博士研究生（今日16:30线上宣讲会）
+1/2
-
广东以色列理工学院博士/硕士招生-通过稀疏观测用数据驱动方法预测湍流
+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