Molcas Forum

tararia

Member

Registered: 2015-12-03

Posts: 13

Active Space expansion calculation issue

Dear team,
                  I am attempting SA-CASSCF on mononuclear uranium (III) complex. I tried to expand active space comprising occupied 5f and unoccupied 6d orbitals. Initially, I tried to use following input :
Input looks like :

&RASSCF
Spin
4
Inactive
78
Nactel
3 1 1 
Ras2
7
Ras3
5
Alter
5
1 97 86
1 98 87
1 102 88
1 103 89
1 104 90
CiRoot
35 35 1

But, the RASSCF output could not converge with the following out put error :

###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    Large orbital rotation.                                              ###
 ###                                                                         ###
 ###                                                                         ###
 ###############################################################################
 ###############################################################################
      Molecular orbital  88 of symmetry 1 MO space 2  weight is    0.260486
 ###############################################################################
 ###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    Large orbital rotation.                                              ###
 ###                                                                         ###
 ###                                                                         ###
 ###############################################################################
 ###############################################################################
      Molecular orbital  89 of symmetry 1 MO space 2  weight is    0.351700
 ###############################################################################
 ###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    Large orbital rotation.                                              ###
 ###                                                                         ###
 ###                                                                         ###
 ###############################################################################
 ###############################################################################
      Molecular orbital  90 of symmetry 1 MO space 2  weight is    0.341123
 ###############################################################################
 ###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    Large orbital rotation.                                              ###
 ###                                                                         ###
 ###                                                                         ###
 ###############################################################################
 ###############################################################################
      Molecular orbital  94 of symmetry 1 MO space 3  weight is    0.292376
       84 500  100    0  -337.49905012    0.32E+00*  0.56E+00*  73  85 1  0.61E-01*  2.81   0.99    QN   YES    0:00:46
        No convergence in the CI section: MAXJT will be increased to                   500
       85 500  100    0  -334.31005635    0.32E+01*  0.64E-01   69 143 1 -0.34E+00*  0.38   0.00    SX    NO    0:01:00

      ========================================================================================================================
 ###############################################################################
 ###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    Rasscf energy diverges.                                              ###
 ###                                                                         ###
 ###                                                                         ###
                                                                                                   

Then, when I meticulously studied MOLCAS forum issues, I noticed with similar RASSCF energy divergence one problem was reported, and in compliance with that suggestions, I tried to incorporate the extended empty 6d orbitals into RAS2 space itself (considering RAS1 and RAS3 space as zero):
Input looks like :

&RASSCF
Spin
4
Inactive
78
Nactel
3 0 0
Ras2
12
Alter
5
1 97 86
1 98 87
1 102 88
1 103 89
1 104 90
CiRoot
35 35 1

Howeevr, even though, I ended up with the similar previous error:

 ###############################################################################
 ###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    Large orbital rotation.                                              ###
 ###                                                                         ###
 ###                                                                         ###
 ###############################################################################
 ###############################################################################
      Molecular orbital  88 of symmetry 1 MO space 2  weight is    0.260486
 ###############################################################################
 ###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    Large orbital rotation.                                              ###
 ###                                                                         ###
 ###                                                                         ###
 ###############################################################################
 ###############################################################################
      Molecular orbital  89 of symmetry 1 MO space 2  weight is    0.351700
 ###############################################################################
 ###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    Large orbital rotation.                                              ###
 ###                                                                         ###
 ###                                                                         ###
 ###############################################################################
 ###############################################################################
      Molecular orbital  90 of symmetry 1 MO space 2  weight is    0.341123
 ###############################################################################
 ###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    Large orbital rotation.                                              ###
 ###                                                                         ###
 ###                                                                         ###
 ###############################################################################
 ###############################################################################
      Molecular orbital  94 of symmetry 1 MO space 3  weight is    0.292376
       84 500  100    0  -337.49905012    0.32E+00*  0.56E+00*  73  85 1  0.61E-01*  2.81   0.99    QN   YES    0:00:46
        No convergence in the CI section: MAXJT will be increased to                   500
       85 500  100    0  -334.31005635    0.32E+01*  0.64E-01   69 143 1 -0.34E+00*  0.38   0.00    SX    NO    0:01:00

      ========================================================================================================================
 ###############################################################################
 ###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    Rasscf energy diverges.                                              ###
 ###                                                                         ###
                                                                                                                       

Please help me out how can I resolve this issue? Please help me out with your invaluable suggestions.

thanks

Last edited by tararia (2017-07-12 09:55:43)

Ignacio

Administrator

From: Uppsala

Registered: 2015-11-03

Posts: 1,085

Re: Active Space expansion calculation issue

Those are not necessarily errors, it just means the initial active space is not "stable" and is changing during the calculation. Does the calculation eventually converge? How does the active space look at the end (it could give you a hint about why the active space wants to change)? Have you tried with a minimal basis set?

tararia

Member

Registered: 2015-12-03

Posts: 13

Re: Active Space expansion calculation issue

Hi,
     No,  the calculation did not converge successfully and ended up with error.  I tried with ANO-RCC-VDZ for all the elements.

 ###############################################################################
 ###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    Something is wrong: Did CI fail?                                     ###
 ###                                                                         ###
 ###                                                                         ###
 ###############################################################################
 ###############################################################################
 ###############################################################################
 ###############################################################################
 ###                                                                         ###
 ###                                                                         ###
 ###    There were warnings during the execution                             ###
 ###    Please, check the output with care!                                  ###
 ###                                                                         ###
 ###                                                                         ###
 ###############################################################################
 ###############################################################################
--- Stop Module:  rasscf at Wed Jul 12 09:52:28 2017 /rc= _GENERAL_ERROR_ ---
--- Module rasscf spent 2 hours and 4 minutes
    ...................................................................................................

Please suggest how should I proceed?

thanks

quan.phung

Member

Registered: 2015-12-18

Posts: 2

Re: Active Space expansion calculation issue

Have you tried to ask for just one root? Then using the converged RASORB as a starting orbital for the calculation with 35 roots?