24小时热门版块排行榜    

查看: 3692  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
信息提示
请填处理意见