| 查看: 4895 | 回复: 5 | |||
[求助]
lammps中LJ势能下单位问题 已有1人参与
|
|
求助一个问题，抱着一点小希望能有前辈给予指点。金币不多，海涵。。。 在lammps中，系统中有两种原子A和B，原子间势能都是L-J势，系统单位为units lj 由于在lj单位制下，物理量都是无量纲的，其中 长度单位就是LJ势中的参数 sigma 此时比如我定义了盒子大小为[0 0 0 ]-[1 1 1]，那盒子的边长是 1*sigma A 还是 sigma B 呢 ？ 还是说，我只能用real的单位制？ 读便手册也没找到答案 哎。。。。 |
回复此楼» 猜你喜欢
306求调剂
已经有5人回复
-
284求调剂
已经有4人回复
-
寻找调剂
已经有4人回复
-
一志愿中南化学（0703）总分337求调剂
已经有10人回复
-
324求调剂
已经有6人回复
-
一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。
已经有3人回复
-
289求调剂
已经有7人回复
-
319求调剂
已经有4人回复
-
298求调剂一志愿211
已经有3人回复
-
求调剂一志愿海大，0703化学学硕304分，有大创项目，四级已过
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求H2S的LJ势能参数
已经有3人回复
- LAMMPS中描述C-Si-O势能
已经有5人回复
- MS粗粒化模型如何导入lammps
已经有10人回复
- 多原子分子LJ势参数sigma 和epsilon怎么算
已经有9人回复
- lammps中界面的建立
已经有9人回复
- 有没有用 lj势函数 模拟lammps的朋友？
已经有6人回复
- 请帮忙看一下lammps中是否有这个势函数？
已经有3人回复
- 请问gromacs导入EAM势能函数的问题和lammps导入martini粗粒化力场的问题
已经有4人回复
- ms建模指定力场导入lammps后是否需要引入势函数
已经有3人回复
- 分子模拟中体系势能和总能量为很大的正值
已经有8人回复
- 新手。。。求教。。如何定义改动的LJ势函数
已经有11人回复
- 关于towhee 回答zhang_jaj的问题
已经有44人回复
- 如何在lammps里加入新的势函数形式？
已经有6人回复
- lammps里面的pair_coeff设置问题
已经有3人回复
- lammps的多尺度计算
已经有15人回复
- 【分享】LAMMPS的语法中文解释
已经有137人回复
- 【求助】分子的LJ势能参数的拟合
已经有8人回复
- 【讨论】lammps运行eam/fs出错
已经有38人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
ktwyw
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 11 (小学生)
- 金币: 2727.6
- 散金: 500
- 红花: 41
- 帖子: 162
- 在线: 114.6小时
- 虫号: 143500
- 注册: 2005-12-22
- 专业: 高分子物理与高分子物理化
2楼2014-10-30 18:41:34
3楼2014-10-31 10:28:10
ktwyw
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 11 (小学生)
- 金币: 2727.6
- 散金: 500
- 红花: 41
- 帖子: 162
- 在线: 114.6小时
- 虫号: 143500
- 注册: 2005-12-22
- 专业: 高分子物理与高分子物理化
4楼2014-10-31 19:24:05
5楼2014-11-01 00:40:48
|
请问用lj势可以不涉及集体的原子么？如ar,H,C等。不涉及两种的具体原子该怎么写？我不知道怎么写，就用了具体的晶格常数，我这样写的对么？ # LJ 12-6 system with npt ensemble units lj atom_style atomic boundary p p p lattice fcc 1.073 lattice fcc 0.8 region box block 0 10 0 5 0 5 create_box 2 box create_atoms 1 box create_atoms 2 box mass 1 1.0 mass 2 1.0 velocity all create 0.01 872877 timestep 0.01 dump 1 all xyz 1000 lj.xyz pair_style lj/cut 2.5 pair_coeff 1 1 1.0 1.0 2.5 pair_coeff 2 2 1.6 1.1765 2.5 pair_coeff 1 2 0.2 1.0882 2.5 #pair_modify shift yes mix geometric neighbor 0.3 bin neigh_modify every 10 delay 0 check yes thermo 1000 fix 1 all npt temp 0.01 0.01 1.00 iso 0.0024 0.0024 100.0 run 50000 unfix 1 thermo 1000 运行结果总出现nan怎么回事？我想不用到集体的原子怎么写？谢谢 |
6楼2017-04-15 21:27:30
