24小时热门版块排行榜    

查看: 1994  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cj4566

木虫 (正式写手)

[求助] ONIOM的Error imposing constraints错误

最近用ONIOM对氧化锌cluster优化时总出错,高手出来指点一下吧,谢谢了。
我按照小木虫上的经验,先将体系都设成low层,用分子动力学进行了测试,和分层时候的错误是一样的。我的输入文件:
%chk=5.chk
%mem=1500MB
%nprocs=6
#p opt pm3 iop(5/13=1) test

Zn60O60(100)

0 1
O               -1   -6.56590000    1.62460000    4.59430000 L
O               -1   -6.56590000    4.87390000    4.59430000 L
O                0   -3.75190000    3.24930000    4.59430000 L
O                0   -0.93800000    4.87390000    4.59430000 L
O               -1   -6.56590000    8.12320000    4.59430000 L
O                0   -3.75190000    6.49850000    4.59430000 L
O                0   -0.93800000    8.12320000    4.59430000 L
O               -1   -6.56590000   11.37240000    4.59430000 L
O                0   -3.75190000    9.74780000    4.59430000 L
O                0   -0.93800000   11.37240000    4.59430000 L
O                0   -3.75190000   12.99710000    4.59430000 L
O                0   -0.93800000   14.62170000    4.59430000 L
Zn              -1   -7.50390000    0.00000000    5.20540000 L
O               -1   -7.50390000    0.00000000    7.19700000 L
Zn               0   -4.68990000    1.62460000    5.20540000 L
O                0   -4.68990000    1.62460000    7.19700000 L
Zn              -1   -6.56590000    1.62460000    7.80820000 L
O               -1   -6.56590000    1.62460000    9.79980000 L
Zn              -1   -7.50390000    3.24930000    5.20540000 L
O               -1   -7.50390000    3.24930000    7.19700000 L
Zn               0   -4.68990000    4.87390000    5.20540000 L
O                0   -4.68990000    4.87390000    7.19700000 L
Zn               0   -1.87600000    3.24930000    5.20540000 L
O                0   -1.87600000    3.24930000    7.19700000 L
Zn              -1   -6.56590000    4.87390000    7.80820000 L
O               -1   -6.56590000    4.87390000    9.79980000 L
Zn               0   -3.75190000    3.24930000    7.80820000 L
O                0   -3.75190000    3.24930000    9.79980000 L
Zn               0   -0.93800000    4.87390000    7.80820000 L
O                0   -0.93800000    4.87390000    9.79980000 L
Zn               0   -7.50390000    6.49850000    5.20540000 L
O               -1   -7.50390000    6.49850000    7.19700000 L
Zn               0   -4.68990000    8.12320000    5.20540000 L
O                0   -4.68990000    8.12320000    7.19700000 L
Zn               0   -1.87600000    6.49850000    5.20540000 L
O                0   -1.87600000    6.49850000    7.19700000 L
Zn              -1   -6.56590000    8.12320000    7.80820000 L
O               -1   -6.56590000    8.12320000    9.79980000 L
Zn               0   -3.75190000    6.49850000    7.80820000 L
O                0   -3.75190000    6.49850000    9.79980000 L
Zn               0   -0.93800000    8.12320000    7.80820000 L
O                0   -0.93800000    8.12320000    9.79980000 L
Zn               0   -7.50390000    9.74780000    5.20540000 L
O               -1   -7.50390000    9.74780000    7.19700000 L
Zn               0   -4.68990000   11.37240000    5.20540000 L
O                0   -4.68990000   11.37240000    7.19700000 L
Zn               0   -1.87600000    9.74780000    5.20540000 L
O                0   -1.87600000    9.74780000    7.19700000 L
Zn              -1   -6.56590000   11.37240000    7.80820000 L
O               -1   -6.56590000   11.37240000    9.79980000 L
Zn               0   -3.75190000    9.74780000    7.80820000 L
O                0   -3.75190000    9.74780000    9.79980000 L
Zn               0   -0.93800000   11.37240000    7.80820000 L
O                0   -0.93800000   11.37240000    9.79980000 L
Zn               0   -1.87600000   12.99710000    5.20540000 L
O                0   -1.87600000   12.99710000    7.19700000 L
Zn               0   -3.75190000   12.99710000    7.80820000 L
O                0   -3.75190000   12.99710000    9.79980000 L
Zn               0   -0.93800000   14.62170000    7.80820000 L
O                0   -0.93800000   14.62170000    9.79980000 L
Zn              -1   -7.50390000    0.00000000   10.41090000 L
O               -1   -7.50390000    0.00000000   12.40250000 L
Zn               0   -4.68990000    1.62460000   10.41090000 L
O                0   -4.68990000    1.62460000   12.40250000 L
Zn              -1   -6.56590000    1.62460000   13.01360000 L
O               -1   -6.56590000    1.62460000   15.00520000 L
Zn              -1   -7.50390000    3.24930000   10.41090000 L
O               -1   -7.50390000    3.24930000   12.40250000 L
Zn               0   -4.68990000    4.87390000   10.41090000 L
O                0   -4.68990000    4.87390000   12.40250000 L
Zn               0   -1.87600000    3.24930000   10.41090000 L
O                0   -1.87600000    3.24930000   12.40250000 L
Zn              -1   -6.56590000    4.87390000   13.01360000 L
O               -1   -6.56590000    4.87390000   15.00520000 L
Zn               0   -3.75190000    3.24930000   13.01360000 L
O                0   -3.75190000    3.24930000   15.00520000 L
Zn               0   -0.93800000    4.87390000   13.01360000 L
O                0   -0.93800000    4.87390000   15.00520000 L
Zn              -1   -7.50390000    6.49850000   10.41090000 L
O               -1   -7.50390000    6.49850000   12.40250000 L
Zn               0   -4.68990000    8.12320000   10.41090000 L
O                0   -4.68990000    8.12320000   12.40250000 L
Zn               0   -1.87600000    6.49850000   10.41090000 L
O                0   -1.87600000    6.49850000   12.40250000 L
Zn              -1   -6.56590000    8.12320000   13.01360000 L
O               -1   -6.56590000    8.12320000   15.00520000 L
Zn               0   -3.75190000    6.49850000   13.01360000 L
O                0   -3.75190000    6.49850000   15.00520000 L
Zn               0   -0.93800000    8.12320000   13.01360000 L
O                0   -0.93800000    8.12320000   15.00520000 L
Zn              -1   -7.50390000    9.74780000   10.41090000 L
O               -1   -7.50390000    9.74780000   12.40250000 L
Zn               0   -4.68990000   11.37240000   10.41090000 L
O                0   -4.68990000   11.37240000   12.40250000 L
Zn               0   -1.87600000    9.74780000   10.41090000 L
O                0   -1.87600000    9.74780000   12.40250000 L
Zn              -1   -6.56590000   11.37240000   13.01360000 L
O               -1   -6.56590000   11.37240000   15.00520000 L
Zn               0   -3.75190000    9.74780000   13.01360000 L
O                0   -3.75190000    9.74780000   15.00520000 L
Zn               0   -0.93800000   11.37240000   13.01360000 L
O                0   -0.93800000   11.37240000   15.00520000 L
Zn               0   -1.87600000   12.99710000   10.41090000 L
O                0   -1.87600000   12.99710000   12.40250000 L
Zn               0   -3.75190000   12.99710000   13.01360000 L
O                0   -3.75190000   12.99710000   15.00520000 L
Zn               0   -0.93800000   14.62170000   13.01360000 L
O                0   -0.93800000   14.62170000   15.00520000 L
Zn              -1   -7.50390000    0.00000000   15.61630000 L
Zn               0   -4.68990000    1.62460000   15.61630000 L
Zn              -1   -7.50390000    3.24930000   15.61630000 L
Zn               0   -4.68990000    4.87390000   15.61630000 L
Zn               0   -1.87600000    3.24930000   15.61630000 L
Zn              -1   -7.50390000    6.49850000   15.61630000 L
Zn               0   -4.68990000    8.12320000   15.61630000 L
Zn               0   -1.87600000    6.49850000   15.61630000 L
Zn              -1   -7.50390000    9.74780000   15.61630000 L
Zn               0   -4.68990000   11.37240000   15.61630000 L
Zn               0   -1.87600000    9.74780000   15.61630000 L
Zn               0   -1.87600000   12.99710000   15.61630000 L


输出文件:
New curvilinear step failed, DQL= 7.02D+00 SP=-2.06D-01.
OIter  1 Iteration     1 Try  1 RMS(Cart)=  0.07083751 RMS(Int)=  0.16326873
OIter  1 Iteration     2 Try  1 RMS(Cart)=  0.01446703 RMS(Int)=  0.10854532
OIter  1 Iteration     3 Try  1 RMS(Cart)=  0.01023580 RMS(Int)=  0.15224388
OIter  1 Iteration     4 Try  1 RMS(Cart)=  0.01403792 RMS(Int)=  0.13226878
Old curvilinear step not converged, using linear step
OIter  1 Iteration     1 Try  1 RMS(Cart)=  0.00036482 RMS(Int)=  0.00046129
OIter  1 Iteration     2 Try  1 RMS(Cart)=  0.00001781 RMS(Int)=  0.00048441
OIter  1 Iteration     3 Try  1 RMS(Cart)=  0.00001657 RMS(Int)=  0.00060400
OIter  1 Iteration     4 Try  1 RMS(Cart)=  0.00003574 RMS(Int)=  0.00105301
Old curvilinear step not converged, using linear step
Error imposing constraints
Error termination via Lnk1e in /pkg/suse11/gaussian/g09/l103.exe at Thu Jun 16 12:05:02 2011.
我冻结了部分原子,怀疑是冻结的问题,但是文献也是冻结的,文献优化完的结构几乎没变。

---------------
如果不冻结原子,结构变化很大,不是我想要的了,而且出现下面错误:
Error in internal coordinate system.
Error termination via Lnk1e in /pkg/suse11/gaussian/g09/l103.exe at Thu Jun 16 22:37:39 2011.

[ Last edited by cj4566 on 2011-6-16 at 20:41 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冬天里的骄阳

铁杆木虫 (正式写手)

【答案】应助回帖

cj4566(金币+10): 太强了,谢谢! 2011-06-17 16:03:27
一般L103这样的错误,就是坐标的问题,g03是没有这样的问题的,g09才有,关键词写上opt=cartesian就没问题
7楼2011-06-17 15:13:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

cj4566

木虫 (正式写手)

怎么没人回答呢?
2楼2011-06-17 07:09:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冬天里的骄阳

铁杆木虫 (正式写手)

【答案】应助回帖

cj4566(金币+10): 谢谢! 2011-06-17 11:10:36
参考文献J. Phys. Chem. C, 2009, 113 (4), pp 1474–1485
话说PM3并不是分子动力学吧
你这个结构我优化了一下,没有你说的那个错误
3楼2011-06-17 10:31:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冬天里的骄阳

铁杆木虫 (正式写手)

【答案】应助回帖

补充:
lz的冻结选择并不好,起码我按照你给定的格式优化,结构变化很大,按参考文献再认真选择合适的冻结原子吧
4楼2011-06-17 10:34:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601求调剂总分293英一数二 +3 钢铁大炮 2026-03-24 3/150 2026-03-24 22:03 by bingxueer79
[考研] 材料与化工考研调剂 +7 孅華 2026-03-22 7/350 2026-03-24 21:04 by greychen00
[考研] 086003食品工程求调剂 +5 淼淼111 2026-03-24 5/250 2026-03-24 20:53 by lailaisimei
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +4 曼111 2026-03-24 5/250 2026-03-24 15:44 by 星空星月
[考研] 277分求调剂,跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 【双一流院校新能源、环境材料,材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 7/350 2026-03-24 11:23 by 种大树
[考研] 一志愿哈工大,085400,320,求调剂 +3 gdlf9999 2026-03-24 3/150 2026-03-24 10:08 by 搏击518
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +5 jiajunX 2026-03-22 5/250 2026-03-23 18:07 by YMU施老师
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 306求调剂 +9 chuanzhu川烛 2026-03-18 9/450 2026-03-23 13:17 by luoyongfeng
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
信息提示
请填处理意见