计算化学公社

标题: opt=readopt冻结原子优化失败 [打印本页]

作者
Author: BinWang 时间: 2025-8-30 17:02
标题: opt=readopt冻结原子优化失败
各位老师：
   我按照帖子http://bbs.keinsci.com/thread-9022-1-1.html，中讲的，冻结了部分原子优化，输入文件如下：
%chk=Au11.chk
%nproc=16
%mem=60GB
# opt=readopt blyp genecp  SCRF=(PCM,Solvent=Acetonitrile)

Au11-cluster

3 1
Au                9.77986084 11.12905299 7.76181044
Au             10.00547431 10.59981187 10.21203617
Au                9.15261962 12.90117902 9.70753372
Au             11.74171098 12.20232157 9.43037118
Au             11.66227505 9.47253661 8.77788352
Au                9.07147758 8.66101460 9.04767320
Au             12.27761408 10.35859962 11.32287991
Au             10.37578264 12.14394804 12.05256930
Au                7.57253316 10.78703775 9.60047497
Au                8.50281300 10.12329785 12.07379419
Au             10.47528574 8.42242246 11.35948015
P                9.40914629 12.10243256 5.89483251
P                8.35521369 14.76606492 9.04306263
P                13.50364093 13.23996483 8.86032546
。。。。。。。。。。。此处太长省略
。。。。。。。。。。
H                12.90567990 7.62444442 10.97016636

notatoms=C,H

Au 0
Lanl2DZ
****
P C H 0
STO-3G
****

Au 0
Lanl2DZ

我的目的是冻结前面所有的Au原子，结果发现优化过程中Au-Au的键长发生了轻微的变化，比如：Au（1-4）=2.79Anstong，优化了六步以后成了2.84Anstong，是这个方法有Bug吗？还是我输入上有错。我用这个基组是因为团簇太大做最初步的优化，各位老师勿喷

作者
Author: BinWang 时间: 2025-8-30 22:38
用第一种方法冻结没有原子移动的问题

，直接在原子符号后面添加-1

作者
Author: sobereva 时间: 2025-8-31 07:58
BLYP描述当前体系非常不适合。哪怕PBE都是更好的选择
STO-3G预优化毫无意义，对P连d极化都没有，结果一塌糊涂。预优化也至少得选个结果还能看得过去的，要不然还不如不预优化
notatoms=C,H冻结的是C和H，和你的目的完全相反

作者
Author: BinWang 时间: 2025-9-1 10:18

sobereva 发表于 2025-8-31 07:58
BLYP描述当前体系非常不适合。哪怕PBE都是更好的选择
STO-3G预优化毫无意义，对P连d极化都没有，结果一塌 ...

，非常感谢Sob老师，是我太矬了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)