应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 【求助】dlpoly算CO2
北京石油化工学院2026年研究生招生接收调剂公告
81/1返回列表
查看: 1029  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

paladinlee

铁杆木虫 (正式写手)

[交流] 【求助】dlpoly算CO2 已有4人参与

现在用dlpoly算CO2，我用类似水的模型，不过角度为180.类似水的sample里Constrain C-O，C-O，O-O。出现error105：shake algorithm failed to converage。
请问怎么解决？或者谁有更好的模型。
field文件如下：
UNITS K
MOLECULES        1
CO2
NUMMOLS      4000
ATOMS       3
C         12.01115     0.5888
O1        15.9994     -0.2944
O1        15.9994     -0.2944
CONSTRAINTS  3
1       2       1.163
1       3       1.163
2       3       2.326
FINISH
VDW       3
O1      O1      lj  82.656   3.000  
C       C       lj  28.845   2.7918
C       O1      lj  48.8284  2.8959
close
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
我是胖人，不是粗人
1楼 2010-12-17 14:49:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zh1987hs(金币+3):谢谢 2010-12-17 19:53:34
由于你键长键角都固定了，可以使用rigid body，不需要使用constraint。具体设置可以查看说明书。
出现error 105，一般是时间步长设置太大了，可以减小时间步试试。
赞一下(2人) 回复此楼
2楼2010-12-17 16:50:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paladinlee

铁杆木虫 (正式写手)
改小了试算，还是有问题，我感觉是CO2力场模型的问题。O-O键该怎么考虑？
赞一下 回复此楼
我是胖人，不是粗人
3楼2010-12-18 09:39:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+3):谢谢 2010-12-18 20:34:52
限制O-O 键的目的是为了固定O-C-O键

把
CONSTRAINTS  3
1       2       1.163
1       3       1.163
2       3       2.326
换成
rigid 1
3 1 2 3
试试
不解释，含义自己看说明书

[ Last edited by 043114076 on 2010-12-18 at 15:18 ]
赞一下(2人) 回复此楼
4楼2010-12-18 15:16:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paladinlee

铁杆木虫 (正式写手)
谢谢！
我现在用angle固定O-O180度。CONSTRAINT C-O键，已经算了一半了...
我先看看结果吧~
赞一下 回复此楼
我是胖人，不是粗人
5楼2010-12-19 10:35:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meixinqing

金虫 (著名写手)

梅心晴｜一只蘑菇
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ghcacj(金币+2): 谢谢 2011-02-05 20:46:34
据developer说，CO2没法做rigid的constraint 因为是linear。。。设置constraint的时候总会出问题。。。我明天再去确认一下。。。。
赞一下(2人) 回复此楼
由于诡异的原因msn偶尔上不去，由于众所周知的原因qq基本不登陆，私聊的pmgmail地址加gtalk/*emuch是个很好很和谐的论坛*/
6楼2011-02-04 06:01:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paladinlee

铁杆木虫 (正式写手)
这确实是个问题，我用angle算的误差太大，rigid算的时间长的令人无法接受。
赞一下 回复此楼
我是胖人，不是粗人
7楼2011-02-04 07:38:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saphy

银虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
御剑江湖(金币+1): 谢谢 2011-04-19 22:31:44
如果使用CONSTRAINTS的话，最好先跑一个仅含本体CO2的MD模拟来优化一下CO2的结构
将优化后的CO2分子插入到材料跑MD时，就不会再出shake的错误
赞一下(2人) 回复此楼
8楼2011-04-19 21:03:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 paladinlee 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学求调剂 一志愿沪9，326分 +4 刘墨墨 2026-04-06 4/200 2026-04-06 18:43 by macy2011
[考研] 070300化学279求调剂 +16 哈哈哈^_^ 2026-03-31 19/950 2026-04-06 14:14 by 无际的草原
[考研] 0703调剂 +13 拾玖壹 2026-04-04 15/750 2026-04-06 10:35 by 蓝云思雨
[考研] 372分，材料与化工，一志愿湖南大学，求调剂 +3 蓝笺片 2026-04-01 3/150 2026-04-06 09:04 by 无际的草原
[考研] 求调剂 +5 chenxrlkx 2026-04-05 7/350 2026-04-06 07:54 by houyaoxu
[考研] 找调剂 +10 楚乔乔 2026-04-01 10/500 2026-04-05 22:19 by syh9288
[考研] 0860 求调剂 一志愿国科大 348 分 +3 WiiiP 2026-04-03 3/150 2026-04-05 17:43 by Ecowxq666！
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 材料专硕322分 +13 哈哈哈吼吼吼哈 2026-04-01 13/650 2026-04-03 16:08 by 哦哦123
[考研] 求调剂 +4 15064154688 2026-04-03 5/250 2026-04-03 15:07 by zrongyan
[考研] 319求调剂 +18 太容易1018 2026-04-01 18/900 2026-04-03 11:18 by linyelide
[考博] 申博求助 +3 Reee1Llll 2026-04-01 3/150 2026-04-02 22:29 by 这是一个无聊的�
[考研] 初试301，代码085701环境工程，本硕一致，四六级已过，有二区一作，共发表5篇论文 +6 axibli 2026-04-01 6/300 2026-04-02 13:42 by Ecowxq666！
[考研] 314求调剂 +11 1xiaojun23 2026-03-31 12/600 2026-04-02 12:31 by 1xiaojun23
[考研] 08开头看过来！！！ +4 wwwwffffff 2026-03-31 6/300 2026-04-02 11:42 by 均值回归
[考研] 284求调剂 +12 小熊～～ 2026-03-31 12/600 2026-04-01 20:23 by 花??
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 一志愿浙江大学工科动力工程370,数一121,专业课135，现在能去哪里 +3 080700调剂 2026-03-30 4/200 2026-03-31 12:00 by KLMY666
[考研] 一志愿食品科学与工程083200求调剂 +4 XQTJZ 2026-03-30 4/200 2026-03-31 04:10 by fmesaito
[考研] 哈尔滨工业大学材料与化工专硕378求调剂 +3 塔比乌斯 2026-03-30 3/150 2026-03-30 22:55 by 无际的草原
信息提示
关闭
请填处理意见
关闭