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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 用modredundant限定原子距离错误
北京石油化工学院2026年研究生招生接收调剂公告
81/1返回列表
查看: 1185  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

mcv

木虫 (小有名气)

[求助] 用modredundant限定原子距离错误

我把配合物的金属变成假原子Bq，然后限定Bq和N的键长键角做opt，这没问题
但我要在两个金属之间的空挡加点东西以保证结构不会向中间塌陷
于是在中间加了个Ar原子，然后用同样的方法固定Ar的位置，输入文件如下
但是这样之后，算一段时间就提示 Error in internal coordinate system的错误了

输入文件
%chk=p7013.chk
%mem=900MW
%nprocshared=8
# opt=modredundant b3lyp/6-31g(d) nosymm geom=connectivity

HUMFIF

0 1
Bq                 2.18600000   10.68130000    4.37690000
Bq                 4.35130000   10.70660000   -2.29340000
N                  1.64300000   10.75860000    6.29980000
N                  4.23450000   10.81010000    5.17930000
N                  1.29640000   12.57560000    4.23130000
N                  1.54000000    8.68010000    4.42090000
N                  4.99040000   10.63860000   -4.17870000
N                  3.88440000   12.71850000   -2.69540000
N                  3.54220000    8.80450000   -2.66750000
N                  6.49000000   10.52280000   -1.78600000
C                  2.19140000   13.65490000    2.12600000
C                  3.63720000    7.39070000    0.70950000
H                  4.54250000    7.10620000    0.67040000
C                  2.99530000    7.75230000   -0.46730000
C                  3.00610000    7.42780000    1.93030000
H                  3.46810000    7.15480000    2.71320000
C                  3.42660000   13.96510000    1.59510000
H                  4.14220000   14.20240000    2.17500000
C                  1.96290000   13.75350000    3.60620000
H                  2.83470000   13.89070000    4.05150000
H                  1.41100000   14.55550000    3.78480000
C                  2.63240000   13.54770000   -0.66060000
C                  1.17920000   13.30320000    1.23060000
H                  0.31370000   13.10740000    1.57070000
C                  2.78880000   13.54770000   -2.14320000
H                  1.93670000   13.23600000   -2.53950000
H                  2.92870000   14.48110000   -2.43680000
C                  1.03870000    8.23410000    0.86360000
H                  0.14220000    8.54860000    0.90280000
C                  1.68150000    7.87100000    2.02330000
C                  5.80520000   10.33260000    3.28080000
C                  3.63490000   13.93220000    0.21530000
H                  4.48480000   14.17950000   -0.13210000
C                  6.38290000   11.93800000    1.57310000
H                  6.49160000   12.84580000    1.31620000
C                  0.48000000   11.67930000    6.35120000
H                  0.28180000   11.91510000    7.29070000
H                 -0.31590000   11.23470000    5.96710000
C                  0.91970000    7.88100000    3.33220000
H                  4.49480000   11.69780000    5.08150000
H                  0.00890000    8.23410000    3.16340000
H                  0.81950000    6.94610000    3.64530000
C                  1.67710000    8.14690000   -0.35470000
H                  1.19350000    8.36560000   -1.14500000
C                  6.09580000    9.32630000    2.39520000
H                  5.99580000    8.42280000    2.66670000
C                  1.38830000   13.22890000   -0.10520000
H                  0.67830000   12.95730000   -0.67280000
C                  1.23080000    9.37910000    6.73290000
H                  0.53090000   12.43830000    3.71870000
H                  2.02280000    8.88450000    7.06320000
H                  0.57930000    9.44630000    7.47660000
C                  3.69110000    7.62370000   -1.79330000
H                  4.65380000    7.46500000   -1.63180000
H                  3.32980000    6.83170000   -2.26550000
C                  6.53690000    9.61070000    1.08630000
H                  6.75280000    8.90310000    0.48930000
C                  3.79140000   12.74850000   -4.17140000
H                  2.93230000   12.35400000   -4.45760000
H                  2.28640000    8.19120000    4.68270000
H                  3.82090000   13.68630000   -4.48450000
C                  2.79950000   11.24180000    7.08760000
H                  2.64810000   11.04740000    8.04670000
H                  2.88190000   12.22250000    6.98490000
C                  0.79300000   12.94150000    5.56100000
H                 -0.02620000   13.49050000    5.47050000
H                  1.47360000   13.47760000    6.04050000
C                  4.89860000   11.99800000   -4.77080000
H                  5.74770000   12.48270000   -4.61910000
H                  4.75510000   11.92510000   -5.74690000
C                  4.09380000    9.68650000   -4.89550000
H                  4.49350000    9.44630000   -5.76890000
H                  3.21920000   10.11680000   -5.06680000
C                  3.89630000    8.44570000   -4.06860000
H                  3.17420000    7.89530000   -4.46250000
H                  4.72840000    7.90960000   -4.07350000
C                  5.35990000    9.99670000    4.67530000
H                  5.09690000    9.04320000    4.70230000
H                  6.12900000   10.11110000    5.28700000
H                  4.66160000   13.18880000   -2.49550000
C                  6.66190000   10.94730000    0.67280000
C                  0.61320000    8.63150000    5.59280000
H                  0.45170000    7.69090000    5.85460000
H                 -0.25370000    9.04180000    5.35300000
C                  5.94760000   11.64350000    2.85020000
H                  5.73970000   12.35400000    3.44470000
C                  6.40430000   10.17250000   -4.18600000
H                  6.82290000   10.40550000   -5.05210000
H                  6.42600000    9.18760000   -4.08820000
H                  2.62300000    8.95170000   -2.69610000
C                  4.05180000   10.58850000    6.64240000
H                  4.00920000    9.61790000    6.83070000
H                  4.82310000   10.96740000    7.13660000
C                  7.18290000   10.80720000   -3.06310000
H                  8.10100000   10.43840000   -3.03860000
H                  7.24600000   11.78500000   -3.20500000
C                  7.12230000   11.29470000   -0.70950000
H                  6.94980000   12.25540000   -0.86610000
H                  8.10280000   11.15610000   -0.75840000
H                  6.64240000    9.62650000   -1.59510000
Ar                 3.26865000   10.69395000    1.04175000

1
2
3 36 1.0 49 1.0 62 1.0
4 40 1.0 77 1.0 91 1.0
5 19 1.0 50 1.0 65 1.0
6 39 1.0 60 1.0 82 1.0
7 68 1.0 71 1.0 87 1.0
8 25 1.0 58 1.0 80 1.0
9 53 1.0 74 1.0 90 1.0
10 94 1.0 97 1.0 100 1.0
11 17 2.0 19 1.0 23 1.5
12 13 1.0 14 1.5 15 2.0
13
14 43 2.0 53 1.0
15 16 1.0 30 1.5
16
17 18 1.0 32 1.5
18
19 20 1.0 21 1.0
20
21
22 25 1.0 32 2.0 47 1.5
23 24 1.0 47 2.0
24
25 26 1.0 27 1.0
26
27
28 29 1.0 30 2.0 43 2.0
29
30 39 1.0
31 45 2.0 77 1.0 85 1.5
32 33 1.0
33
34 35 1.0 81 2.0 85 2.0
35
36 37 1.0 38 1.0 65 1.0
37
38
39 41 1.0 42 1.0
40
41
42
43 44 1.0
44
45 46 1.0 56 1.5
46
47 48 1.0
48
49 51 1.0 52 1.0 82 1.0
50
51
52
53 54 1.0 55 1.0
54
55
56 57 1.0 81 1.5
57
58 59 1.0 61 1.0 68 1.0
59
60
61
62 63 1.0 64 1.0 91 1.0
63
64
65 66 1.0 67 1.0
66
67
68 69 1.0 70 1.0
69
70
71 72 1.0 73 1.0 74 1.0
72
73
74 75 1.0 76 1.0
75
76
77 78 1.0 79 1.0
78
79
80
81 97 1.0
82 83 1.0 84 1.0
83
84
85 86 1.0
86
87 88 1.0 89 1.0 94 1.0
88
89
90
91 92 1.0 93 1.0
92
93
94 95 1.0 96 1.0
95
96
97 98 1.0 99 1.0
98
99
100
101

B 1 3 F
B 1 101 F
B 101 2 F
A 1 101 2 F
B 1 4 F
B 1 5 F
B 1 6 F
B 2 7 F
B 2 8 F
B 2 9 F
B 2 10 F
A 1 4 40 F
A 1 5 50 F
A 1 6 60 F
A 2 8 80 F
A 2 9 90 F
A 2 10 100 F
A 1 3 36 F
A 1 3 62 F
A 1 3 49 F
A 2 7 68 F
A 2 7 71 F
A 2 7 87 F


出错部分：
SCF Done:  E(RB3LYP) =  -2369.10591723     A.U. after   14 cycles
...
Population analysis using the SCF
...
Dipole moment (field-independent basis, Debye)
...
Center     Atomic                   Forces (Hartrees/Bohr)
...
Cartesian Forces:  Max     0.136805184 RMS     0.060172568

GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
NTrRot=    -1 NTRed=  2936 NAtoms=   101 NSkip=  2639 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e
回复此楼

» 猜你喜欢
1楼 2013-06-27 00:07:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
gkf高: 金币+2 2013-06-27 10:54:58
Ar原子虽然是惰性原子，但它不是假原子，它是一个电子已经完全配对的原子，具有极强的电子排斥力，它在体系中波函数是要参加计算的，即使是形成范德华作用，也需要具有范德华半径，现在加在结构中固定，自然说坐标有误，不能计算。
赞一下(1人) 回复此楼
2楼2013-06-27 08:25:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mcv

木虫 (小有名气)
引用回帖:
2楼: Originally posted by zhou2009 at 2013-06-27 08:25:25
Ar原子虽然是惰性原子，但它不是假原子，它是一个电子已经完全配对的原子，具有极强的电子排斥力，它在体系中波函数是要参加计算的，即使是形成范德华作用，也需要具有范德华半径，现在加在结构中固定，自然说坐标有 ...

配合物本身是两个金属之间配位一个阴离子，也就是M-X-M
我现在想只看不同M...M距离下，配体的构象不同，于是删掉阴离子，然后添加惰性原子占位，如果要达成这个目标，有什么建议么
另外，原来没添加Ar的时候，我是用B 1 2 S 5 0.300000一次扫一段M...M距离，现在虽然没成功，但是貌似高斯没有d1=d2的功能，只能手动一次算一个，有更好的方案么
赞一下 回复此楼
3楼2013-06-27 08:49:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
你上述还没有把问题说清楚。
赞一下 回复此楼
4楼2013-06-27 09:54:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mcv

木虫 (小有名气)
gkf高: 金币+3 2013-06-27 10:55:11
gkf高: 金币-3, 粗心了~ 2013-06-27 10:57:14
就是L是个笼子，里面可以包双核金属，然后金属之间又可以结合一个阴离子
我现在想看配合物在包结不同尺寸阴离子时，L的构象以及不同构象的能量（ccdc有一部分不同M...M的结构）
所以把金属原子换成不会因为靠的太近而产生相互作用的假原子，并且通过固定M-N的键长键角保证配位环境不变
之前用B 1 2 S 5 0.300000一次扫一段M...M距离（1 2是M，即Bq的序号）
但是由于中间阴离子被删了，所以苯环在优化的过程中会扭向L的中间

所以现在在中间加了个稀有气体原子，但是增加固定Bq-Ar的键长键角之后就在log里面说Berny optimization后Error in internal coordinate system.
赞一下 回复此楼
5楼2013-06-27 10:05:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

L是个笼子，我在GV中看见了。
里面包双核金属，然后金属之间又可以结合一个阴离子，要分析它们与笼子的化学作用，才形成笼子的形状，把金属和阴离子都去掉了就纯粹只是笼子了，自然有它自己的构象。加Ar是不行的。
笼子中间装什么，之间都有作用的，去掉，整个体系就全变了，没有可比性了。反正这样的考察设计不合理。
赞一下 回复此楼
6楼2013-06-27 11:19:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mcv

木虫 (小有名气)
是这样，配合物内部的M-X的键级分析以及X和苯环的弱作用我都算过了
然后我想分析下笼子的间隔基团不同及其在不同M...M下L的能量，主要是从纯配体的构象上分析，也就是之前文献中一直提到但没细说的配体的刚柔性
之前有人做过了间位苯环的那个L随着配位阴离子的尺寸不同log K的大小趋势，他说不同M...M距离对识别比较关键，只有合适大小才有最大的log K

我就是想算出比较合适的Bq...Bq距离和他对比，并且找出未报道的对位苯环的合适的M...M距离
赞一下 回复此楼
7楼2013-06-27 11:32:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mcv

木虫 (小有名气)
继续google发现这个问题的产生是由于程序对一条线上的原子会产生问题，然后说g09 c01之后的版本可以解决，然后就去找c01了
赞一下 回复此楼
8楼2013-06-28 15:05:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mcv 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 300求调剂，材料科学英一数二 +6 leaflight 2026-03-24 6/300 2026-03-28 21:45 by sanrepian
[考研] 一志愿北京化工大学材料与化工（085600）296求调剂 +11 稻妻小编 2026-03-26 11/550 2026-03-28 21:04 by 唐沐儿
[考研] 求调剂一志愿海大，0703化学学硕304分，有大创项目，四级已过 +7 幸运哩哩 2026-03-22 11/550 2026-03-28 20:27 by 唐沐儿
[考研] 0856材料化工调剂 总分330 +11 zhubinhao 2026-03-27 11/550 2026-03-28 20:22 by 418490947
[考研] 食品工程专硕一志愿中海洋309求调剂 +4 小张zxy张 2026-03-26 8/400 2026-03-28 19:25 by lbsjt
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-27 8/400 2026-03-28 12:42 by 唐沐儿
[考研] 调剂 +3 好好读书。 2026-03-28 3/150 2026-03-28 12:04 by 王保杰33
[考研] 材料277求调剂 +7 min3 2026-03-24 7/350 2026-03-28 11:39 by xuxiang
[考研] 291求调剂 +7 孅華 2026-03-22 7/350 2026-03-28 04:02 by fmesaito
[考研] 275求调剂 +10 jjjjjjjjjjl 2026-03-27 10/500 2026-03-27 23:47 by barnett0632
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +7 18419759900 2026-03-25 8/400 2026-03-27 16:38 by 18419759900
[考博] 26申博 +3 加油冲啊！ 2026-03-26 3/150 2026-03-27 15:38 by cls512
[考研] 考研化学308分求调剂 +10 你好明天你好 2026-03-23 12/600 2026-03-27 14:43 by shangxh
[考研] 调剂推荐 +5 清酒714 2026-03-26 6/300 2026-03-27 11:12 by 不吃魚的貓
[考研] 材料学硕333求调剂 +8 北道巷 2026-03-24 8/400 2026-03-27 10:18 by 我是小康
[考研] 材料学硕，求调剂 6+5 糖葫芦888ll 2026-03-22 10/500 2026-03-27 08:18 by hypershenger
[考研] 318求调剂 +5 plum李子 2026-03-21 8/400 2026-03-25 09:26 by aa331100
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
信息提示
关闭
请填处理意见
关闭