24小时热门版块排行榜    

查看: 3281  |  回复: 11

herarysara

木虫 (小有名气)

[求助] lammps计算中出现原子距离太近现象

在用lammps计算气体反应时,发现体系中有些分子中键长太短。如氢气分子键长只有0.5,O-H只有0.7等等,体系温度很高;求解决办法。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meihai3166

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
温度太高只有缩小时间步长了。。这应该是惟一可能的方法,要么就不用
ToughGuy
2楼2013-03-15 00:11:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meihai3166

木虫 (小有名气)

引用回帖:
2楼: Originally posted by meihai3166 at 2013-03-15 00:11:28
温度太高只有缩小时间步长了。。这应该是惟一可能的方法,要么就不用

要么就只能用其他方法,这分子动力学的硬伤。
ToughGuy
3楼2013-03-15 00:15:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)

引用回帖:
2楼: Originally posted by meihai3166 at 2013-03-15 00:11:28
温度太高只有缩小时间步长了。。这应该是惟一可能的方法,要么就不用

原来是0.1已经够小了,改成0.01都出问题,看来没得救了。
4楼2013-03-15 08:53:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)

引用回帖:
3楼: Originally posted by meihai3166 at 2013-03-15 00:15:33
要么就只能用其他方法,这分子动力学的硬伤。...

其他方法是嘛意思?
5楼2013-03-15 08:54:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

倔强的坚果墙

铁杆木虫 (正式写手)

看你的势参数是否能准确表达原子间的作用势,这个不准备跑MD也会出现温度上升很大的情况
像小强一样活着
6楼2013-03-15 09:26:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meihai3166

木虫 (小有名气)

引用回帖:
5楼: Originally posted by herarysara at 2013-03-15 08:54:09
其他方法是嘛意思?...

准连续介质力学之类的,多尺度方面的方法,能跨越原子尺度到达微观尺度的就行
ToughGuy
7楼2013-03-15 09:51:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meihai3166

木虫 (小有名气)

引用回帖:
4楼: Originally posted by herarysara at 2013-03-15 08:53:53
原来是0.1已经够小了,改成0.01都出问题,看来没得救了。...

你做了弛豫没,如果弛豫没问题,势输入错误的可能性不大。
另一方面检查作用势是否合理吧,因为大部分的作用势都是在平衡态下拟合的,而且是常温,高温下的情况是否能描述,这个本来就具有争议,但大多人就这么用的。

» 本帖已获得的红花(最新10朵)

ToughGuy
8楼2013-03-15 10:02:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)

引用回帖:
6楼: Originally posted by 倔强的坚果墙 at 2013-03-15 09:26:47
看你的势参数是否能准确表达原子间的作用势,这个不准备跑MD也会出现温度上升很大的情况

感谢你的回复;我体系是很普通的CHO元素,作用势应该没问题,估计高温下可能真不适用。
9楼2013-03-15 10:19:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

herarysara

木虫 (小有名气)

送鲜花一朵
引用回帖:
8楼: Originally posted by meihai3166 at 2013-03-15 10:02:27
你做了弛豫没,如果弛豫没问题,势输入错误的可能性不大。
另一方面检查作用势是否合理吧,因为大部分的作用势都是在平衡态下拟合的,而且是常温,高温下的情况是否能描述,这个本来就具有争议,但大多人就这么用 ...

感谢回复,我做过弛豫的,而且是缓慢升温,到目标温度后用的nve,当体系开始反应后,急剧升温达到上万K,感觉是温度失控了,不知道怎样控制;如果是高温下,怎么判定势场的适用性?
10楼2013-03-15 10:21:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 herarysara 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 326求调剂 +3 上岸的小葡 2026-03-15 4/200 2026-03-15 18:50 by 无际的草原
[考研] 311求调剂 +5 冬十三 2026-03-15 5/250 2026-03-15 18:38 by 无际的草原
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 337一志愿华南理工材料求调剂(有希望2吗?) +3 mysdl 2026-03-09 3/150 2026-03-14 02:53 by JourneyLucky
[考研] 求调剂! +4 朔朔话 2026-03-09 4/200 2026-03-14 01:38 by JourneyLucky
[考研] 考研材料与化工,求调剂 +8 戏精丹丹丹 2026-03-09 8/400 2026-03-14 01:14 by JourneyLucky
[考研] 材料工程专硕,一志愿中国矿业大学,总分314,求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 327分求调剂086 +4 西红柿?小帅 2026-03-09 7/350 2026-03-10 14:47 by ruiyingmiao
[考研] 一志愿山东大学,总分327,英语二79,有论文,有竞赛,已过四六级 +3 木木目目1 2026-03-09 3/150 2026-03-09 19:52 by yuningshan
信息提示
请填处理意见