版块导航
正在加载中...
客户端APP下载
论文辅导
申博辅导
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求,自2017年10月1日起,未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用,请尽快对帐号进行手机号验证,感谢您的理解与支持!
24小时热门版块排行榜
>
论坛更新日志
(3521)
>
文献求助
(402)
>
导师招生
(214)
>
虫友互识
(176)
>
休闲灌水
(128)
>
论文道贺祈福
(101)
>
考博
(89)
>
招聘信息布告栏
(78)
>
论文投稿
(73)
>
硕博家园
(72)
>
基金申请
(49)
>
绿色求助(高悬赏)
(35)
>
博后之家
(34)
>
外文书籍求助
(22)
>
教师之家
(22)
>
SciFinder/Reaxys
(20)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
分子模拟
»
LAMMPS/DL_POLY
»
lammps中LJ势能下单位问题
4
1/1
返回列表
查看: 4748 | 回复: 5
查看全部回帖
@他人
存档
新回复提醒
(忽略)
收藏
在APP中查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖
jidanaoyun
新虫
(初入文坛)
应助: 0
(幼儿园)
金币: 29.1
帖子: 16
在线: 32.1小时
虫号: 2160782
注册: 2012-12-02
专业: 工程热力学
[
求助
]
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的单位制?
读便手册也没找到答案 哎。。。。
回复此楼
» 猜你喜欢
北京211副教授,35岁,想重新出发,去国外做博后,怎么样?
已经有5人回复
论文终于录用啦!满足毕业条件了
已经有25人回复
2026年机械制造与材料应用国际会议 (ICMMMA 2026)
已经有3人回复
自荐读博
已经有3人回复
求助:我三月中下旬出站,青基依托单位怎么办?
已经有5人回复
不自信的我
已经有5人回复
磺酰氟产物,毕不了业了!
已经有4人回复
投稿Elsevier的杂志(返修),总是在选择OA和subscription界面被踢皮球
已经有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人回复
1楼
2014-10-30 15:59:30
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
jidanaoyun
新虫
(初入文坛)
应助: 0
(幼儿园)
金币: 29.1
帖子: 16
在线: 32.1小时
虫号: 2160782
注册: 2012-12-02
专业: 工程热力学
引用回帖:
2楼
:
Originally posted by
ktwyw
at 2014-10-30 18:41:34
A和B是由你决定的,你愿意让哪个做单位都行,关键是单位选定之后,做单位的那个就是1,另外一个相应改变。比如sigma_A=2sigma_B,如果A作单位,B就是0.5;如果B作单位,A就是2。另外设置的时候注意各自的截断半径以 ...
感谢前辈赐教,我不知道怎么设置A为1,是不是如下命令就能将sigmaA设为1呢?
pair_coeff 1 1 1.0 1.0 2.5
pair_coeff 2 2 1.4564 0.0252 2.5
pair_modify mix arithmetic
如果我的第一行改为:pair_coeff 1 1 2.0 2.0 2.5,其他参数不变,那盒子边长又应该是多少呢(仍以上述盒子边长[0 0 0]-[1 1 1]为例)
赞
一下
回复此楼
3楼
2014-10-31 10:28:10
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
jidanaoyun
新虫
(初入文坛)
应助: 0
(幼儿园)
金币: 29.1
帖子: 16
在线: 32.1小时
虫号: 2160782
注册: 2012-12-02
专业: 工程热力学
引用回帖:
4楼
:
Originally posted by
ktwyw
at 2014-10-31 19:24:05
还是1sigma,不过这时的sigma既不对应A也不对应B。另外你设置cutoff的时候要具体看各种例子的直径,2.5不是万能的
大体有点懂,还是略糊涂
我已经把程序全改成了metal单位制的了,最起码可以运行了。。
谢谢解惑!
现在出了个问题,我fix npt,run之后温度远远偏离要维持的温度(fix npt是同时对气液两相做的约束),请问前辈知不知道应该怎么解决啊谢谢啦
赞
一下
回复此楼
5楼
2014-11-01 00:40:48
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
相关版块跳转
第一性原理
量子化学
计算模拟
分子模拟
仿真模拟
程序语言
我要订阅楼主
jidanaoyun
的主题更新
4
1/1
返回列表
如果回帖内容含有宣传信息,请如实选中。否则帐号将被全论坛禁言
普通表情
龙
兔
虎
猫
百度网盘
|
360云盘
|
千易网盘
|
华为网盘
在新窗口页面中打开自己喜欢的网盘网站,将文件上传后,然后将下载链接复制到帖子内容中就可以了。
信息提示
关闭
请填处理意见
关闭
确定