Multiwfn forum

anupthms000

Member

Registered: 2020-05-23

Posts: 3

Regarding Dynamic and nondynamic correlation index

Can we get the nondynamic and dynamic correlation index from Unrestricted DFT studies using multiwfn. 

Look forward to hear from you.

Regards
Anup

sobereva

Tian Lu (Multiwfn developer)

From: Beijing

Registered: 2017-09-11

Posts: 1,468

Website

Re: Regarding Dynamic and nondynamic correlation index

Usually, the two indices should be calculated using multi-configurational methods, as illustrated in Section 4.A.6 of Multiwfn manual, since non-integer natural orbital occupations can be obtained in this case.

In my opinion, if you use UKS to calculate open-shell singlet system such as biradical via symmetry-broken formalism, and then obtain the so-called UNO (unrestricted natural orbital, which is generated by diagonalizing total density matrix of UKS wavefunction), the occupation numbers may also be used in nondynamic correlation index to judge nondynamic correlation, but I didn't have a try. More specifically, after loading an .fch/fchk file containing UKS wavefunction, you should input commands like follows

...ignored
 14 Adaptive natural density partitioning (AdNDP) analysis
 15 Fuzzy atomic space analysis
 16 Charge decomposition analysis (CDA) and plot orbital interaction diagram
 17 Basin analysis                    18 Electron excitation analysis
 19 Orbital localization analysis     20 Visual study of weak interaction
 21 Energy decomposition analysis
 100 Other functions (Part 1)         200 Other functions (Part 2)
 300 Other functions (Part 3)
200

               ============ Other functions (Part 2) ============
 0 Return
 1 Calculate core-valence bifurcation (CVB) index and related quantities
 2 Calculate atomic and bond dipole moments in Hilbert space
 3 Generate cube file for multiple orbital wavefunctions
 4 Generate iso-chemical shielding surfaces (ICSS) and related quantities
 5 Plot radial distribution function for a real space function
 6 Analyze correspondence between orbitals in two wavefunctions
 7 Parse output of (hyper)polarizability task of Gaussian
 8 Study (hyper)polarizability by sum-over-states (SOS) method
 9 Calculate average bond length and average coordinate number
 10 Output various kinds of integral between orbitals
 11 Calculate center, the first and second moments of a real space function
 12 Calculate energy index (EI) or bond polarity index (BPI)
 13 Evaluate orbital contributions to density difference or other grid data
 14 Domain analysis (Obtaining properties within isosurfaces of a function)
 15 Calculate electron correlation index (PCCP, 18, 24015)
 16 Generate natural orbitals based on the density matrix in .fch/.fchk file
 17 Calculate Coulomb and exchange integrals between two orbitals
 18 Calculate bond length/order alternation (BLA/BOA)
 20 Bond order density (BOD) and natural adaptive orbital (NAdO) analyses
16
 Input type of density, e.g. SCF, MP2, CI, CC, MP4...
 e.g. If the .fch was produced under MP2, you may input "SCF" or "MP2"
SCF
 Select natural orbitals you want to obtain
 1 Spatial natural orbitals (diagonalization of total density matrix)
 2 Alpha and beta natural orbitals (diagonalization of respective DM)
 3 Spin natural orbitals (diagonalization of spin density matrix)
1
 Loading density matrix...
 Density matrix was loaded from .fch/.fchk file

 Generating NOs, please wait...
 Occupation numbers:
    2.000000    2.000000    1.999999    1.999999    1.999977    1.999974
    1.999965    1.999963    1.999845    1.999798    1.999779    1.999773
    1.999188    1.999005    1.998891    1.309524    0.690476    0.001109
    0.000995    0.000812    0.000227    0.000221    0.000202    0.000155
    0.000037    0.000035    0.000026    0.000023    0.000001    0.000001
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
    0.000000    0.000000    0.000000    0.000000    0.000000    0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000   -0.000000
   -0.000000   -0.000000
 Done! Basis function information now correspond to natural orbitals

 If next you intend to analyze real space functions based on the NOs, you should
 export new.molden in current folder and then reload it, so that GTF information
 will also correspond to NOs
 Would you like to do this immediately? (y/n)

After exporting and load the new.molden containing UNO orbitals, you can perform dynamic/non-dynamic correlation index calculation as exemplified in Section 4.A.6 of Multiwfn manual.