计算化学公社

标题: 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/) Powered by Discuz! X3.3