GIC设置不生效

luzujia

我想使两个分子在优化过程中保持平行，关键词为："# B3LYP/6-31+G* OPT(cartesian,GIC) nosymm EM=GD3BJ"，末尾表达式为:
Cx1=X(20)
Cy1=Y(20)
Cz1=Z(20)
Cx2=X(27)
Cy2=Y(27)
Cz2=Z(27)
Cx3=X(33)
Cy3=Y(33)
Cz3=Z(33)
Aplane(freeze,Min=-0.1,Max=0.1)=(Cy2-Cy1)*(Cz3-Cz1)-(Cz2-Cz1)*(Cy3-Cy1)
Bplane(freeze,Min=-0.1,Max=0.1)=(Cz2-Cz1)*(Cx3-Cx1)-(Cx2-Cx1)*(Cz3-Cz1)

同时通过在原子后加“-1”的方式固定第一个分子。需要补充说明的是，这个体系的初始结构，第一个分子几何中心处于坐标原点，两个分子平行同时与XY平面平行。
遇到的问题是优化结果根本不满足我设置的GIC限定条件。于是我逐步简化限制，最后只剩下GIC条件：
Cz1=Z(20) freeze
Cz2=Z(27) freeze
Cz3=Z(33) freeze

就是只限制某几个原子的Z坐标。就这样还是不生效，实在无法找到原因出在哪。请各位老师指教！

下面附上输入文件。

luzujia

lanthanum 发表于 2025-10-11 00:17
opt=addgic

是这样的！谢谢老师！

lanthanum

opt=addgic