molpro 计算遇到的不收敛问题

willrmrex

请教一下大家，在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对于这样的错误提示，一般可以通过增加态平均的态数来解决，但是这样又会导致第一个问题的出现，这个真是不知道该怎么弄？还请大家指点指点。

woshiwhm

您好，请问您说的Error ? No convergence ? The problem occurs in Multi对于这样的错误提示这个是如何解决的

zjxitcc

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

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

woshiwhm

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

wzkchem5

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

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

zjxitcc

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

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

wzkchem5

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

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

小崔

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

wzkchem5

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

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

小崔

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里读取上一步优化的波函数，请老师指教一下。

wzkchem5

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

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