计算化学公社

标题: molpro 计算遇到的不收敛问题 [打印本页]

作者
Author: willrmrex 时间: 2017-2-20 08:57
标题: molpro 计算遇到的不收敛问题
请教一下大家，在molpro 计算的时候经常遇到两个不收敛问题
一个是在优化结构的时候遇到NO CONVERGENCE OF CP-MCSCF AFTER900 ITEATIONS AND 9 RESTARTS; FINAL ACTEST= 0.40D-03
这样的错误提示，这个问题在求助的时候有人给出的建议是
1. This error occurs despite smooth convergence of the MCSCF wave function.2. This error is exacerbated when the number of states in thestate-averaged wave function are increased.
有的时候发现减少态平均的态数时的确能够进一步的计算，但是在这个时候又会遇到另一个问题 ? Error ? No convergence ? The problem occurs in Multi对于这样的错误提示，一般可以通过增加态平均的态数来解决，但是这样又会导致第一个问题的出现，这个真是不知道该怎么弄？还请大家指点指点。

作者
Author: woshiwhm 时间: 2021-7-28 15:38
您好，请问您说的Error ? No convergence ? The problem occurs in Multi对于这样的错误提示这个是如何解决的

作者
Author: zjxitcc 时间: 2021-7-28 16:33

woshiwhm 发表于 2021-7-28 15:38
您好，请问您说的Error ? No convergence ? The problem occurs in Multi对于这样的错误提示这个是如何解决 ...

最好截图展示，信息量更多。

作者
Author: woshiwhm 时间: 2021-7-29 09:41

zjxitcc 发表于 2021-7-28 16:33
最好截图展示，信息量更多。

老师，您好，这是我的报错信息。问题出现在opt这一模块。
想用Molpro的CAS优化构型，优化一个态。麻烦您指点我一下!
输入文件：
basis,def2-tzvp;
{multi;maxiter,40;occ,26;closed,18;wf,44,1,0;state,1;
cpmcscf,grad,1.1,spin=0;
}
optg;
pop;
put,xyz;
输出文件：
Geometry optimization using default procedure for command MULTI

Geometry written to block  1 of record 701

Making model hessian for Z-matrix coordinates

Entering Rational Function Geometry Optimization

ITER. ENERGY(OLD) ENERGY(NEW)    DE       GRADMAX    GRADNORM GRADRMS    STEPMAX    STEPLEN    STEPRMS CPU-time
1  -433.73892615  -433.73910241 -0.00017625  0.00451311  0.00883943  0.00266519  0.03399086  0.03639405  0.00330855    32.20
2  -433.73910241  -433.73919861 -0.00009620  0.00241764  0.00478614  0.00144307  0.05292750  0.05464132  0.00496739    67.12
3  -433.73919861  -433.73920701 -0.00000840  0.00204689  0.00320668  0.00096685  0.01269162  0.01275413  0.00115947 103.94
4  -433.73920701  -433.73921242 -0.00000541  0.00199273  0.00297908  0.00089823  0.01160023  0.01200452  0.00109132 151.87
5  -433.73921242  -433.73921463 -0.00000221  0.00200715  0.00269531  0.00081267  0.00248745  0.00301397  0.00027400 185.04
6  -433.73921463  -433.73922324 -0.00000861  0.00187604  0.00255853  0.00077143  0.00567059  0.00958687  0.00087153 215.77
? Error
? No convergence
? The problem occurs in Multi

GLOBAL ERROR fehler on processor 0

作者
Author: wzkchem5 时间: 2021-7-29 16:27

woshiwhm 发表于 2021-7-29 02:41
老师，您好，这是我的报错信息。问题出现在opt这一模块。
想用Molpro的CAS优化构型，优化一个态。麻烦您 ...

虽然报错在opt，但是报错信息提示是CASSCF不收敛。
尝试一下以结构优化最新输出的结构为初始结构，继续优化，有一定概率能解决。如果还不能解决，贴最后一步结构优化的Multi输出，不要贴opt的输出

作者
Author: zjxitcc 时间: 2021-7-29 17:28

wzkchem5 发表于 2021-7-29 16:27
虽然报错在opt，但是报错信息提示是CASSCF不收敛。
尝试一下以结构优化最新输出的结构为初始结构，继续 ...

他这个只能提取出结构再算个单点，方能看出CASSCF情况如何。结构优化过程中是没有CASSCF输出的。

作者
Author: wzkchem5 时间: 2021-7-29 17:52

zjxitcc 发表于 2021-7-29 10:28
他这个只能提取出结构再算个单点，方能看出CASSCF情况如何。结构优化过程中是没有CASSCF输出的。

有输出吧，输出到log文件，和主要的输出不在一起

作者
Author: 小崔 时间: 2024-11-14 10:42

wzkchem5 发表于 2021-7-29 16:27
虽然报错在opt，但是报错信息提示是CASSCF不收敛。
尝试一下以结构优化最新输出的结构为初始结构，继续 ...

老师，您好，我用molpro的CASSCF优化激发态也出现了这样的报错，我用您说的'尝试一下以结构优化最新输出的结构为初始结构，继续优化'试了一下，出现了新的报错，我不知道在输出文件的哪里看Multi输出，麻烦您指导一下，我是量化新人，第一次使用公社，不规范的地方望您多多海涵。
输入文件：
wf,charge=0,spin=2
hf
{matrop;
read,mo,ORB,file=t.a;
save,mo,2140.2,ORBITALS}
{CASSCF;occ,14;closed,6}
{optg ! 进行scf结构优化
inactive,B1}
.log输出文件：
** WVFN ****  MAXIMUM NUMBER OF ITERATIONS REACHED

Gradient information for state 1.1 saved on record    5000.1

First order charge density matrix for state 1.1 saved on record  2140.2 (density set 1)

Results for state 1.1
=====================

!MCSCF STATE 1.1 Energy          -152.061557616437
Nuclear energy                      32.01678675
Kinetic energy                      152.02845618
One electron energy                -272.97907346
Two electron energy                88.90072910
Virial ratio                         2.00021773

!MCSCF STATE 1.1 Dipole moment       0.00194330 -0.12416778    1.60471908
Dipole moment /Debye                0.00493904 -0.31558235    4.07852192

Natural orbital dump at molpro section  2140.2 (Orbital set 2)

out输出文件：
PROGRAM * OPT (Geometry optimization)    Authors: F. Eckert and H.-J. Werner

BFGS update of hessian

*** Long output written to logfile /home/wangzg/workplace/czh/T1/t1/2/t.log ***

Geometry optimization using default procedure for command CASSCF

Inactive variables: B1

Active variables:    B2 A1 B3 A2 D1 B4
                     A3 D2 B5 A4 D3

Geometry written to block  1 of record 700

Making model hessian for Z-matrix coordinates

Quadratic Steepest Descent - Minimum Search

ITER. ENERGY(OLD) ENERGY(NEW)    DE       GRADMAX    GRADNORM GRADRMS    STEPMAX    STEPLEN    STEPRMS CPU-time
1  -151.97831269  -151.98961657 -0.01130387  0.08346969  0.10308766  0.03108210  0.19259522  0.36127505  0.03284319    92.61
2  -151.98961657  -151.99473043 -0.00511386  0.04600310  0.05059245  0.01525420  0.21678605  0.39770846  0.03615531 822.24
3  -151.99473043  -152.00793972 -0.01320930  0.03644359  0.04393714  0.01324755  0.35054316  0.49737552  0.04521596 1488.91
4  -152.00793972  -151.99356913    0.01437059  0.04885319  0.06687444  0.02016340  0.41373066  0.59877291  0.05443390 1745.75
5  -151.99356913  -152.04013073 -0.04656160  0.25525699  0.26601291  0.08020591  0.33117923  0.51210959  0.04655542 1963.77
6  -152.04013073  -152.05207650 -0.01194576  0.03032119  0.03551300  0.01070757  0.57413387  0.58938523  0.05358048 2178.45
7  -152.05207650  -152.05868160 -0.00660510  0.01475822  0.02101393  0.00633594  0.43188135  0.59275722  0.05388702 2580.70
8  -152.05868160  -152.06147814 -0.00279654  0.01151213  0.01537327  0.00463522  0.56051326  0.59408819  0.05400802 2906.83
no convergence in dis1, ss= 34479.1189927373
9  -152.06147814  -152.06262614 -0.00114800  0.00672327  0.01226785  0.00369890  0.47141013  0.57115496  0.05192318 3352.32
no convergence in dis1, ss= 1.71306811792390
no convergence in dis1, ss= 1.65496819392303
  10  -152.06262614  -152.06347639 -0.00085025  0.01615456  0.01860446  0.00560946  0.23195504  0.28524206  0.02593110 3593.95
no convergence in dis1, ss= 1.59809133780385
no convergence in dis1, ss= 1.53675511482986
  11  -152.06347639  -152.06369089 -0.00021450  0.00317651  0.00485690  0.00146441  0.27763894  0.31272701  0.02842973 3909.90
no convergence in dis1, ss= 1.20580559613949
no convergence in dis1, ss= 1.15020678815215
  12  -152.06369089  -152.06379176 -0.00010088  0.00536931  0.00626921  0.00189024  0.18766654  0.21345682  0.01940517 4072.17
no convergence in dis1, ss= 2.51051892963889
  13  -152.06379176  -152.06290684    0.00088493  0.00345425  0.00539767  0.00162746  0.37024491  0.52029513  0.04729956 4313.15
no convergence in dis1, ss=  0.578158015716827
  14  -152.06290684  -152.06377919 -0.00087236  0.01233575  0.02421531  0.00730119  0.31242666  0.46825988  0.04256908 4565.40
no convergence in dis1, ss= 1.27771661130760
? Error
? No convergence
? The problem occurs in Multi
采用您的办法尝试后的.out输出文件：
Variable memory set to  800000000 words,  buffer space 230000 words

ZSYMEL=NOSYM

SETTING B1          =       3.01636000  ANG
SETTING B2          =       0.98916628  ANG
SETTING B3          =       3.33133843  ANG
SETTING B4          =       0.94211489  ANG
SETTING B5          = 0.94212566ANG
SETTING A1          =       3.85601642  DEGREE
SETTING A2          =    162.13804551  DEGREE
SETTING A3          =    125.65374693  DEGREE
SETTING A4          =    125.68164401  DEGREE
SETTING D1          =       -0.30245021  DEGREE
SETTING D2          =    -101.49545360  DEGREE
SETTING D3          =    101.86484726  DEGREE
SETTING BASIS       = USERDEF

Recomputing integrals since basis changed

Using spherical harmonics

? Possible error in Z matrix
Line                   6  connection                   1  expr=
0.94212566ANG
? Apparently depends on a variable, but not one that is yet defined

Z matrix    Current Z-matrix
========

Records: 700.0 700.0
Maximum, total, and current numbers of geometry points 100 0 0
Unique atoms/degeneracy:

#    Tag    N1 Distance  N2 Angle    N3 Angle       J
Grp Chg    Distance    Angle       Angle                x          y          z

1  1 O1
   1  8.0 0.0000000 0.0000000 0.0000000             0.0000000 0.0000000 0.0000000
2  2 O2    1 B1
   1  8.0 3.0163600 0.0000000 0.0000000             0.0000000 0.0000000 0.0000000
3  3 H1    1 B2       2 A1
   2  1.0 0.9891663 3.8560164 0.0000000             0.0000000 0.0000000 0.0000000
4  4 H2    1 B3       2 A2       3 D1          0
   2  1.0 3.3313384  162.1380455 -0.3024502             0.0000000 0.0000000 0.0000000
5  5 H3    2 B4       1 A3       3 D2          0
   2  1.0 0.9421149  125.6537469 -101.4954536             0.0000000 0.0000000 0.0000000
6  6 H4    2 0.942125 1 A4       3 D3          0
   2  1.0 0.0000000  125.6816440  101.8648473             0.0000000 0.0000000 0.0000000

Variables on which Z-matrix depends:
Name:                B1             B2             A1             B3             A2             D1             B4             A3
Value:       0.0000000 A    0.0000000 A    0.0000000 D    0.0000000 A    0.0000000 D    0.0000000 D    0.0000000 A    0.0000000 D
Name:                D2             B3             A2             D1             B4             A3             D2
Value:       0.0000000 D    0.0000000 A    0.0000000 D    0.0000000 D    0.0000000 A    0.0000000 D    0.0000000 D    0.0000000
Name:                A4             D3
Value:       0.0000000 D    0.0000000 D
? Error
? Error in Z matrix
? The problem occurs in zmatrix.f:zmat_evaluate

作者
Author: wzkchem5 时间: 2024-11-14 11:16

小崔发表于 2024-11-14 03:42
老师，您好，我用molpro的CASSCF优化激发态也出现了这样的报错，我用您说的'尝试一下以结构优化最新输出 ...

尝试读取最后一步的波函数继续计算
你后来重新优化报错应该是因为复制错了结构，有的键长都等于0了

作者
Author: 小崔 时间: 2024-11-16 12:36

wzkchem5 发表于 2024-11-14 11:16
尝试读取最后一步的波函数继续计算
你后来重新优化报错应该是因为复制错了结构，有的键长都等于0了

老师，我用.log输出文件中优化的最后结构作为mokit选活性空间的初始结构，.log里的最后结构如下：
Current geometry (xyz format, in Angstrom)

6
CASSCF000/USERDEF  ENERGY=-152.06377919
O       0.0000000000       0.0000000000       0.0000000000
O       0.0000000000       0.0000000000       3.0163600000
H       0.0000000000       0.0665208288       0.9869269986
H       -0.0053938233       1.0217894549    -3.1707622071
H       0.7501636169    -0.1525605246       3.5655050636
H       -0.7489112655    -0.1573405340       3.5658840104

老师，我不会在molpro里读取上一步优化的波函数，请老师指教一下。

作者
Author: wzkchem5 时间: 2024-11-16 15:03

小崔发表于 2024-11-16 05:36
老师，我用.log输出文件中优化的最后结构作为mokit选活性空间的初始结构，.log里的最后结构如下：
Curre ...

在说明书里搜read、guess等关键词。这种问题必须会自己解决

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