Static versus dynamic correlation correlation#
Bond-breaking#
The most commonly used example to illustrate static vs dynamical correlation is H
import multipsi as mtp
import numpy as np
import matplotlib.pyplot as plt
import veloxchem as vlx
bohr_to_Å = 0.529177
# HF calculation of H, equivalent to CI
mol_str = "H 0.0000 0.0000 0.0000"
molecule = vlx.Molecule.read_molecule_string(mol_str, units='angstrom')
molecule.set_multiplicity(2)
basis = vlx.MolecularBasis.read(molecule, "cc-pvdz")
scf_drv = vlx.ScfUnrestrictedDriver()
scf_results = scf_drv.compute(molecule, basis)
E_H_hf = scf_drv.get_scf_energy()
# H2 scan
mol_template = """
H 0.0000 0.0000 0.0000
H 0.0000 0.0000 dist
"""
scf_drv = vlx.ScfRestrictedDriver()
ci_drv = mtp.CIDriver()
distlist = [0.5, 0.6, 0.7, 0.75, 0.8,0.9,1.0,1.2,1.5,2,2.5,3,3.5,4,4.5,5]
E_hf = []
E_FCI = []
# Scan over O-H distances
for dist in distlist:
mol_str = mol_template.replace("dist", str(dist))
molecule = vlx.Molecule.read_molecule_string(mol_str, units='angstrom')
basis = vlx.MolecularBasis.read(molecule, "cc-pvdz")
scf_results = scf_drv.compute(molecule, basis)
E_hf.append(scf_drv.get_scf_energy() - 2*E_H_hf)
space=mtp.OrbSpace(molecule,scf_drv.mol_orbs)
space.fci()
ci_results = ci_drv.compute(molecule,basis,space)
E_FCI.append(ci_results["energies"][0] - 2*E_H_hf)
# Save orbitals and densities for equilibrium value
if dist == 0.75:
eq_molecule = molecule
eq_orbitals = scf_drv.mol_orbs
eq_2density = ci_drv.get_active_2body_density(0)
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Unrestricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.0000000000 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.01 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -0.499278403420 a.u. Time: 0.02 sec.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -0.499278403420 0.0000000000 0.00000005 0.00000003 0.00000000
*** SCF converged in 1 iterations. Time: 0.01 sec.
Spin-Unrestricted Hartree-Fock:
-------------------------------
Total Energy : -0.4992784034 a.u.
Electronic Energy : -0.4992784034 a.u.
Nuclear Repulsion Energy : 0.0000000000 a.u.
------------------------------------
Gradient Norm : 0.0000000462 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 2
Magnetic Quantum Number (M_S) : 0.5
Expectation value of S**2 : 0.7500
Spin Unrestricted Alpha Orbitals
--------------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 1.000 Energy: -0.49928 a.u.
( 1 H 1s : 0.59) ( 1 H 2s : 0.50)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: 0.57244 a.u.
( 1 H 1s : 1.24) ( 1 H 2s : -1.28)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 1.44303 a.u.
( 1 H 1p0 : 1.00)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 1.44303 a.u.
( 1 H 1p+1: 1.00)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.44303 a.u.
( 1 H 1p-1: -1.00)
Spin Unrestricted Beta Orbitals
-------------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 0.000 Energy: 0.07335 a.u.
( 1 H 1s : -0.22) ( 1 H 2s : -0.84)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: 0.75329 a.u.
( 1 H 1s : 1.35) ( 1 H 2s : -1.09)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 1.57364 a.u.
( 1 H 1p0 : 1.00)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 1.57364 a.u.
( 1 H 1p+1: 1.00)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.57364 a.u.
( 1 H 1p-1: 1.00)
Ground State Dipole Moment
----------------------------
X : 0.000000 a.u. 0.000000 Debye
Y : 0.000000 a.u. 0.000000 Debye
Z : 0.000000 a.u. 0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 1.0583544218 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -1.044195881437 a.u. Time: 0.03 sec.
* Info * Overlap matrix computed in 0.01 sec.
* Info * Kinetic energy matrix computed in 0.01 sec.
* Info * Nuclear potential matrix computed in 0.01 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -1.048687184260 0.0000000000 0.03228141 0.00976213 0.00000000
2 -1.048800398179 -0.0001132139 0.00075218 0.00019687 0.00838820
3 -1.048800556215 -0.0000001580 0.00000986 0.00000323 0.00041723
4 -1.048800556227 -0.0000000000 0.00000001 0.00000000 0.00000326
*** SCF converged in 4 iterations. Time: 0.04 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -1.0488005562 a.u.
Electronic Energy : -2.1071549780 a.u.
Nuclear Repulsion Energy : 1.0583544218 a.u.
------------------------------------
Gradient Norm : 0.0000000051 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.67746 a.u.
( 1 H 1s : 0.43) ( 2 H 1s : 0.43)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: 0.23819 a.u.
( 1 H 2s : 2.99) ( 2 H 2s : -2.99)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.42154 a.u.
( 1 H 1s : -0.62) ( 1 H 2s : 0.72) ( 2 H 1s : -0.62)
( 2 H 2s : 0.72)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 1.16876 a.u.
( 1 H 1s : 0.50) ( 1 H 2s : -1.87) ( 1 H 1p0 : -0.52)
( 2 H 1s : -0.50) ( 2 H 2s : 1.87) ( 2 H 1p0 : -0.52)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.33275 a.u.
( 1 H 1p+1: -0.54) ( 2 H 1p+1: -0.54)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.33275 a.u.
( 1 H 1p-1: -0.54) ( 2 H 1p-1: -0.54)
Ground State Dipole Moment
----------------------------
X : -0.000000 a.u. -0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : 0.000000 a.u. 0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : -0.00 (multiplicity = 1.0 )
- Energy : -1.0793700508928015
- Natural orbitals
1.97719 0.00008 0.00501 0.00970 0.00380 0.00380 0.00013 0.00013 0.00013 0.00001
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.8819620182 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -1.102601746873 a.u. Time: 0.01 sec.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -1.106788495584 0.0000000000 0.03045659 0.00841139 0.00000000
2 -1.106892338747 -0.0001038432 0.00069812 0.00020665 0.00725834
3 -1.106892475138 -0.0000001364 0.00000920 0.00000317 0.00039011
4 -1.106892475152 -0.0000000000 0.00000001 0.00000000 0.00000386
*** SCF converged in 4 iterations. Time: 0.04 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -1.1068924752 a.u.
Electronic Energy : -1.9888544933 a.u.
Nuclear Repulsion Energy : 0.8819620182 a.u.
------------------------------------
Gradient Norm : 0.0000000057 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.63879 a.u.
( 1 H 1s : 0.42) ( 2 H 1s : 0.42)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: 0.22260 a.u.
( 1 H 2s : -2.47) ( 2 H 2s : 2.47)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.44547 a.u.
( 1 H 1s : -0.66) ( 1 H 2s : 0.72) ( 2 H 1s : -0.66)
( 2 H 2s : 0.72)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 1.05754 a.u.
( 1 H 1s : -0.58) ( 1 H 2s : 1.73) ( 1 H 1p0 : 0.50)
( 2 H 1s : 0.58) ( 2 H 2s : -1.73) ( 2 H 1p0 : 0.50)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.31355 a.u.
( 1 H 1p+1: 0.55) ( 2 H 1p+1: 0.55)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.31355 a.u.
( 1 H 1p-1: 0.55) ( 2 H 1p-1: 0.55)
Ground State Dipole Moment
----------------------------
X : 0.000000 a.u. 0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : 0.000000 a.u. 0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : -0.00 (multiplicity = 1.0 )
- Energy : -1.1391732932951417
- Natural orbitals
1.97340 0.01319 0.00564 0.00001 0.00360 0.00360 0.00016 0.00015 0.00015 0.00011
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.7559674441 a.u.
* Info * Overlap matrix computed in 0.01 sec.
* Info * Kinetic energy matrix computed in 0.01 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -1.122770998276 a.u. Time: 0.04 sec.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -1.126826026762 0.0000000000 0.02918562 0.00724027 0.00000000
2 -1.126924576362 -0.0000985496 0.00064391 0.00020404 0.00640409
3 -1.126924692283 -0.0000001159 0.00000970 0.00000342 0.00036324
4 -1.126924692303 -0.0000000000 0.00000002 0.00000001 0.00000477
*** SCF converged in 4 iterations. Time: 0.02 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -1.1269246923 a.u.
Electronic Energy : -1.8828921364 a.u.
Nuclear Repulsion Energy : 0.7559674441 a.u.
------------------------------------
Gradient Norm : 0.0000000192 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.60482 a.u.
( 1 H 1s : 0.41) ( 1 H 2s : 0.16) ( 2 H 1s : 0.41)
( 2 H 2s : 0.16)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: 0.20498 a.u.
( 1 H 2s : -2.09) ( 2 H 2s : 2.09)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.46968 a.u.
( 1 H 1s : -0.69) ( 1 H 2s : 0.73) ( 2 H 1s : -0.69)
( 2 H 2s : 0.73)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.96759 a.u.
( 1 H 1s : -0.71) ( 1 H 2s : 1.64) ( 1 H 1p0 : 0.44)
( 2 H 1s : 0.71) ( 2 H 2s : -1.64) ( 2 H 1p0 : 0.44)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.29782 a.u.
( 1 H 1p+1: -0.53) ( 1 H 1p-1: 0.22) ( 2 H 1p+1: -0.53)
( 2 H 1p-1: 0.22)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.29782 a.u.
( 1 H 1p+1: 0.22) ( 1 H 1p-1: 0.53) ( 2 H 1p+1: 0.22)
( 2 H 1p-1: 0.53)
Ground State Dipole Moment
----------------------------
X : 0.000000 a.u. 0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : 0.000000 a.u. 0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : 0.00 (multiplicity = 1.0 )
- Energy : -1.1609046824806644
- Natural orbitals
1.96871 0.01801 0.00602 0.00001 0.00331 0.00331 0.00019 0.00015 0.00015 0.00014
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.7055696145 a.u.
* Info * Overlap matrix computed in 0.01 sec.
* Info * Kinetic energy matrix computed in 0.01 sec.
* Info * Nuclear potential matrix computed in 0.01 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -1.124608367719 a.u. Time: 0.05 sec.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -1.128645660151 0.0000000000 0.02872496 0.00752087 0.00000000
2 -1.128743027026 -0.0000973669 0.00062088 0.00020101 0.00607432
3 -1.128743134729 -0.0000001077 0.00001024 0.00000362 0.00035226
4 -1.128743134753 -0.0000000000 0.00000002 0.00000001 0.00000526
*** SCF converged in 4 iterations. Time: 0.02 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -1.1287431348 a.u.
Electronic Energy : -1.8343127493 a.u.
Nuclear Repulsion Energy : 0.7055696145 a.u.
------------------------------------
Gradient Norm : 0.0000000210 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.58941 a.u.
( 1 H 1s : 0.40) ( 1 H 2s : 0.18) ( 2 H 1s : 0.40)
( 2 H 2s : 0.18)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: 0.19551 a.u.
( 1 H 1s : 0.15) ( 1 H 2s : 1.93) ( 2 H 1s : -0.15)
( 2 H 2s : -1.93)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.48172 a.u.
( 1 H 1s : -0.70) ( 1 H 2s : 0.73) ( 2 H 1s : -0.70)
( 2 H 2s : 0.73)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.93020 a.u.
( 1 H 1s : -0.77) ( 1 H 2s : 1.61) ( 1 H 1p0 : 0.40)
( 2 H 1s : 0.77) ( 2 H 2s : -1.61) ( 2 H 1p0 : 0.40)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.29182 a.u.
( 1 H 1p-1: 0.58) ( 2 H 1p-1: 0.58)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.29182 a.u.
( 1 H 1p+1: 0.58) ( 2 H 1p+1: 0.58)
Ground State Dipole Moment
----------------------------
X : -0.000000 a.u. -0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : 0.000000 a.u. 0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : 0.00 (multiplicity = 1.0 )
- Energy : -1.1635935606538248
- Natural orbitals
1.96588 0.02104 0.00611 0.00001 0.00314 0.00314 0.00020 0.00016 0.00016 0.00016
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.6614715136 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -1.122863755794 a.u. Time: 0.01 sec.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -1.126903786231 0.0000000000 0.02833405 0.00776828 0.00000000
2 -1.127000620590 -0.0000968344 0.00059929 0.00019720 0.00578585
3 -1.127000720853 -0.0000001003 0.00001092 0.00000385 0.00034224
4 -1.127000720883 -0.0000000000 0.00000001 0.00000000 0.00000582
*** SCF converged in 4 iterations. Time: 0.03 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -1.1270007209 a.u.
Electronic Energy : -1.7884722345 a.u.
Nuclear Repulsion Energy : 0.6614715136 a.u.
------------------------------------
Gradient Norm : 0.0000000052 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.57494 a.u.
( 1 H 1s : 0.40) ( 1 H 2s : 0.19) ( 2 H 1s : 0.40)
( 2 H 2s : 0.19)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: 0.18567 a.u.
( 1 H 1s : -0.16) ( 1 H 2s : -1.79) ( 2 H 1s : 0.16)
( 2 H 2s : 1.79)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.49365 a.u.
( 1 H 1s : 0.72) ( 1 H 2s : -0.73) ( 2 H 1s : 0.72)
( 2 H 2s : -0.73)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.89695 a.u.
( 1 H 1s : -0.83) ( 1 H 2s : 1.58) ( 1 H 1p0 : 0.36)
( 2 H 1s : 0.83) ( 2 H 2s : -1.58) ( 2 H 1p0 : 0.36)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.28732 a.u.
( 1 H 1p+1: -0.15) ( 1 H 1p-1: -0.57) ( 2 H 1p+1: -0.15)
( 2 H 1p-1: -0.57)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.28732 a.u.
( 1 H 1p+1: -0.57) ( 1 H 1p-1: 0.15) ( 2 H 1p+1: -0.57)
( 2 H 1p-1: 0.15)
Ground State Dipole Moment
----------------------------
X : -0.000000 a.u. -0.000000 Debye
Y : 0.000000 a.u. 0.000000 Debye
Z : -0.000000 a.u. -0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : -0.00 (multiplicity = 1.0 )
- Energy : -1.1627507630215623
- Natural orbitals
1.96266 0.02457 0.00613 0.00017 0.00296 0.00296 0.00022 0.00016 0.00016 0.00001
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.5879746788 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -1.112220769289 a.u. Time: 0.01 sec.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -1.116294268816 0.0000000000 0.02764470 0.00815308 0.00000000
2 -1.116391213357 -0.0000969445 0.00055877 0.00018803 0.00530657
3 -1.116391300399 -0.0000000870 0.00001254 0.00000437 0.00032420
4 -1.116391300440 -0.0000000000 0.00000000 0.00000000 0.00000699
*** SCF converged in 4 iterations. Time: 0.03 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -1.1163913004 a.u.
Electronic Energy : -1.7043659792 a.u.
Nuclear Repulsion Energy : 0.5879746788 a.u.
------------------------------------
Gradient Norm : 0.0000000013 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.54852 a.u.
( 1 H 1s : -0.39) ( 1 H 2s : -0.21) ( 2 H 1s : -0.39)
( 2 H 2s : -0.21)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: 0.16505 a.u.
( 1 H 1s : 0.18) ( 1 H 2s : 1.55) ( 2 H 1s : -0.18)
( 2 H 2s : -1.55)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.51701 a.u.
( 1 H 1s : -0.75) ( 1 H 2s : 0.73) ( 2 H 1s : -0.75)
( 2 H 2s : 0.73)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.84044 a.u.
( 1 H 1s : 0.91) ( 1 H 2s : -1.52) ( 1 H 1p0 : -0.27)
( 2 H 1s : -0.91) ( 2 H 2s : 1.52) ( 2 H 1p0 : -0.27)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.28349 a.u.
( 1 H 1p+1: 0.60) ( 2 H 1p+1: 0.60)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.28349 a.u.
( 1 H 1p-1: -0.60) ( 2 H 1p-1: -0.60)
Ground State Dipole Moment
----------------------------
X : -0.000000 a.u. -0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : -0.000000 a.u. -0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : 0.00 (multiplicity = 1.0 )
- Energy : -1.1540817061165265
- Natural orbitals
1.95470 0.03339 0.00599 0.00020 0.00257 0.00257 0.00025 0.00016 0.00016 0.00002
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.5291772109 a.u.
* Info * Overlap matrix computed in 0.01 sec.
* Info * Kinetic energy matrix computed in 0.01 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -1.095959327021 a.u. Time: 0.02 sec.
* Info * Overlap matrix computed in 0.01 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -1.100056242572 0.0000000000 0.02692675 0.00838888 0.00000000
2 -1.100153689745 -0.0000974472 0.00051916 0.00017723 0.00493758
3 -1.100153764817 -0.0000000751 0.00001420 0.00000490 0.00030702
4 -1.100153764873 -0.0000000001 0.00000000 0.00000000 0.00000815
*** SCF converged in 4 iterations. Time: 0.04 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -1.1001537649 a.u.
Electronic Energy : -1.6293309758 a.u.
Nuclear Repulsion Energy : 0.5291772109 a.u.
------------------------------------
Gradient Norm : 0.0000000003 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.52504 a.u.
( 1 H 1s : 0.38) ( 1 H 2s : 0.23) ( 2 H 1s : 0.38)
( 2 H 2s : 0.23)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: 0.14359 a.u.
( 1 H 1s : 0.20) ( 1 H 2s : 1.36) ( 2 H 1s : -0.20)
( 2 H 2s : -1.36)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.53961 a.u.
( 1 H 1s : -0.78) ( 1 H 2s : 0.74) ( 2 H 1s : -0.78)
( 2 H 2s : 0.74)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.79480 a.u.
( 1 H 1s : 0.95) ( 1 H 2s : -1.47) ( 1 H 1p0 : -0.20)
( 2 H 1s : -0.95) ( 2 H 2s : 1.47) ( 2 H 1p0 : -0.20)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.28713 a.u.
( 1 H 1p+1: 0.31) ( 1 H 1p-1: 0.55) ( 2 H 1p+1: 0.31)
( 2 H 1p-1: 0.55)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.28713 a.u.
( 1 H 1p+1: 0.55) ( 1 H 1p-1: -0.31) ( 2 H 1p+1: 0.55)
( 2 H 1p-1: -0.31)
Ground State Dipole Moment
----------------------------
X : -0.000000 a.u. -0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : 0.000000 a.u. 0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : -0.00 (multiplicity = 1.0 )
- Energy : -1.140073480871046
- Natural orbitals
1.94416 0.04503 0.00565 0.00022 0.00217 0.00217 0.00027 0.00016 0.00016 0.00002
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.4409810091 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -1.057057778249 a.u. Time: 0.03 sec.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -1.061017008512 0.0000000000 0.02498883 0.00837395 0.00000000
2 -1.061111945190 -0.0000949367 0.00043469 0.00015110 0.00450475
3 -1.061111997775 -0.0000000526 0.00001623 0.00000553 0.00026831
4 -1.061111997850 -0.0000000001 0.00000000 0.00000000 0.00000968
*** SCF converged in 4 iterations. Time: 0.01 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -1.0611119978 a.u.
Electronic Energy : -1.5020930069 a.u.
Nuclear Repulsion Energy : 0.4409810091 a.u.
------------------------------------
Gradient Norm : 0.0000000004 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.48526 a.u.
( 1 H 1s : 0.36) ( 1 H 2s : 0.27) ( 2 H 1s : 0.36)
( 2 H 2s : 0.27)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: 0.10025 a.u.
( 1 H 1s : 0.24) ( 1 H 2s : 1.08) ( 2 H 1s : -0.24)
( 2 H 2s : -1.08)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.58296 a.u.
( 1 H 1s : 0.83) ( 1 H 2s : -0.75) ( 2 H 1s : 0.83)
( 2 H 2s : -0.75)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.72928 a.u.
( 1 H 1s : 0.96) ( 1 H 2s : -1.37) ( 2 H 1s : -0.96)
( 2 H 2s : 1.37)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.31531 a.u.
( 1 H 1p+1: -0.61) ( 1 H 1p-1: 0.26) ( 2 H 1p+1: -0.61)
( 2 H 1p-1: 0.26)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.31531 a.u.
( 1 H 1p+1: 0.26) ( 1 H 1p-1: 0.61) ( 2 H 1p+1: 0.26)
( 2 H 1p-1: 0.61)
Ground State Dipole Moment
----------------------------
X : 0.000000 a.u. 0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : -0.000000 a.u. -0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : -0.00 (multiplicity = 1.0 )
- Energy : -1.1068551925482588
- Natural orbitals
1.91211 0.07966 0.00454 0.00024 0.00142 0.00142 0.00030 0.00015 0.00015 0.00002
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.3527848073 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -0.999019668155 a.u. Time: 0.02 sec.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -1.002120347046 0.0000000000 0.02044695 0.00718048 0.00000000
2 -1.002192720639 -0.0000723736 0.00029841 0.00010488 0.00418857
3 -1.002192745431 -0.0000000248 0.00001282 0.00000447 0.00019427
4 -1.002192745478 -0.0000000000 0.00000000 0.00000000 0.00000804
*** SCF converged in 4 iterations. Time: 0.03 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -1.0021927455 a.u.
Electronic Energy : -1.3549775527 a.u.
Nuclear Repulsion Energy : 0.3527848073 a.u.
------------------------------------
Gradient Norm : 0.0000000003 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.43916 a.u.
( 1 H 1s : -0.33) ( 1 H 2s : -0.32) ( 2 H 1s : -0.33)
( 2 H 2s : -0.32)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: 0.04030 a.u.
( 1 H 1s : 0.28) ( 1 H 2s : 0.81) ( 2 H 1s : -0.28)
( 2 H 2s : -0.81)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.64366 a.u.
( 1 H 1s : 0.91) ( 1 H 2s : -0.76) ( 2 H 1s : 0.91)
( 2 H 2s : -0.76)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.67571 a.u.
( 1 H 1s : -0.94) ( 1 H 2s : 1.25) ( 2 H 1s : 0.94)
( 2 H 2s : -1.25)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.27713 a.u.
( 1 H 1s : 0.24) ( 1 H 1p0 : -0.63) ( 2 H 1s : 0.24)
( 2 H 1p0 : 0.63)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.38302 a.u.
( 1 H 1p+1: -0.67) ( 1 H 1p-1: 0.15) ( 2 H 1p+1: -0.67)
( 2 H 1p-1: 0.15)
Ground State Dipole Moment
----------------------------
X : -0.000000 a.u. -0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : 0.000000 a.u. 0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : 0.00 (multiplicity = 1.0 )
- Energy : -1.0615349496253987
- Natural orbitals
1.82407 0.17146 0.00246 0.00020 0.00026 0.00065 0.00065 0.00011 0.00011 0.00003
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.2645886055 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -0.920706243669 a.u. Time: 0.01 sec.
* Info * Overlap matrix computed in 0.01 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -0.921886356635 0.0000000000 0.01150627 0.00406041 0.00000000
2 -0.921908589697 -0.0000222331 0.00012300 0.00004281 0.00260878
3 -0.921908594102 -0.0000000044 0.00000367 0.00000128 0.00008607
4 -0.921908594106 -0.0000000000 0.00000000 0.00000000 0.00000249
*** SCF converged in 4 iterations. Time: 0.05 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -0.9219085941 a.u.
Electronic Energy : -1.1864971996 a.u.
Nuclear Repulsion Energy : 0.2645886055 a.u.
------------------------------------
Gradient Norm : 0.0000000000 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.38571 a.u.
( 1 H 1s : 0.31) ( 1 H 2s : 0.38) ( 2 H 1s : 0.31)
( 2 H 2s : 0.38)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: -0.03585 a.u.
( 1 H 1s : -0.31) ( 1 H 2s : -0.60) ( 2 H 1s : 0.31)
( 2 H 2s : 0.60)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.65320 a.u.
( 1 H 1s : 0.93) ( 1 H 2s : -1.11) ( 2 H 1s : -0.93)
( 2 H 2s : 1.11)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.67021 a.u.
( 1 H 1s : 0.94) ( 1 H 2s : -0.75) ( 2 H 1s : 0.94)
( 2 H 2s : -0.75)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.34213 a.u.
( 1 H 1p0 : 0.69) ( 2 H 1p0 : -0.69)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.45332 a.u.
( 1 H 1p+1: 0.63) ( 1 H 1p-1: -0.31) ( 2 H 1p+1: 0.63)
( 2 H 1p-1: -0.31)
Ground State Dipole Moment
----------------------------
X : -0.000000 a.u. -0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : -0.000000 a.u. -0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : -0.00 (multiplicity = 1.0 )
- Energy : -1.01759411404388
- Natural orbitals
1.56605 0.43300 0.00006 0.00010 0.00041 0.00014 0.00014 0.00004 0.00004 0.00001
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.2116708844 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -0.864890569654 a.u. Time: 0.02 sec.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -0.865323130286 0.0000000000 0.00662539 0.00233203 0.00000000
2 -0.865330119957 -0.0000069897 0.00002402 0.00000842 0.00153364
3 -0.865330120137 -0.0000000002 0.00000130 0.00000045 0.00001924
4 -0.865330120138 -0.0000000000 0.00000000 0.00000000 0.00000098
*** SCF converged in 4 iterations. Time: 0.05 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -0.8653301201 a.u.
Electronic Energy : -1.0770010045 a.u.
Nuclear Repulsion Energy : 0.2116708844 a.u.
------------------------------------
Gradient Norm : 0.0000000000 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.34934 a.u.
( 1 H 1s : 0.30) ( 1 H 2s : 0.42) ( 2 H 1s : 0.30)
( 2 H 2s : 0.42)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: -0.08466 a.u.
( 1 H 1s : -0.32) ( 1 H 2s : -0.51) ( 2 H 1s : 0.32)
( 2 H 2s : 0.51)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.62463 a.u.
( 1 H 1s : -0.91) ( 1 H 2s : 0.75) ( 2 H 1s : -0.91)
( 2 H 2s : 0.75)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.66777 a.u.
( 1 H 1s : -0.94) ( 1 H 2s : 1.02) ( 2 H 1s : 0.94)
( 2 H 2s : -1.02)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.45319 a.u.
( 1 H 1p0 : 0.71) ( 2 H 1p0 : -0.71)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.46317 a.u.
( 1 H 1p+1: 0.71) ( 2 H 1p+1: 0.71)
Ground State Dipole Moment
----------------------------
X : 0.000000 a.u. 0.000000 Debye
Y : 0.000000 a.u. 0.000000 Debye
Z : -0.000000 a.u. -0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : 0.00 (multiplicity = 1.0 )
- Energy : -1.0031292512270575
- Natural orbitals
1.30747 0.69238 0.00002 0.00001 0.00006 0.00003 0.00003 0.00001 0.00001 0.00000
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.1763924036 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -0.826215752307 a.u. Time: 0.02 sec.
* Info * Overlap matrix computed in 0.01 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -0.826444448804 0.0000000000 0.00461717 0.00162670 0.00000000
2 -0.826447843878 -0.0000033951 0.00000214 0.00000076 0.00110952
3 -0.826447843880 -0.0000000000 0.00000021 0.00000007 0.00000161
*** SCF converged in 3 iterations. Time: 0.04 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -0.8264478439 a.u.
Electronic Energy : -1.0028402475 a.u.
Nuclear Repulsion Energy : 0.1763924036 a.u.
------------------------------------
Gradient Norm : 0.0000002057 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.32328 a.u.
( 1 H 1s : -0.30) ( 1 H 2s : -0.44) ( 2 H 1s : -0.30)
( 2 H 2s : -0.44)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: -0.11576 a.u.
( 1 H 1s : 0.32) ( 1 H 2s : 0.48) ( 2 H 1s : -0.32)
( 2 H 2s : -0.48)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.59484 a.u.
( 1 H 1s : -0.90) ( 1 H 2s : 0.77) ( 2 H 1s : -0.90)
( 2 H 2s : 0.77)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.67845 a.u.
( 1 H 1s : -0.94) ( 1 H 2s : 0.96) ( 2 H 1s : 0.94)
( 2 H 2s : -0.96)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.46177 a.u.
( 1 H 1p+1: 0.55) ( 1 H 1p-1: -0.44) ( 2 H 1p+1: 0.55)
( 2 H 1p-1: -0.44)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.46177 a.u.
( 1 H 1p+1: 0.44) ( 1 H 1p-1: 0.55) ( 2 H 1p+1: 0.44)
( 2 H 1p-1: 0.55)
Ground State Dipole Moment
----------------------------
X : 0.000000 a.u. 0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : 0.000000 a.u. 0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : -0.00 (multiplicity = 1.0 )
- Energy : -0.9995506185969365
- Natural orbitals
1.15006 0.84991 0.00000 0.00000 0.00001 0.00001 0.00001 0.00000 0.00000 0.00000
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.1511934888 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -0.799871798054 a.u. Time: 0.02 sec.
* Info * Overlap matrix computed in 0.01 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -0.799991332224 0.0000000000 0.00320033 0.00112755 0.00000000
2 -0.799992972433 -0.0000016402 0.00000479 0.00000169 0.00078939
3 -0.799992972441 -0.0000000000 0.00000017 0.00000006 0.00000407
*** SCF converged in 3 iterations. Time: 0.05 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -0.7999929724 a.u.
Electronic Energy : -0.9511864613 a.u.
Nuclear Repulsion Energy : 0.1511934888 a.u.
------------------------------------
Gradient Norm : 0.0000001692 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.30398 a.u.
( 1 H 1s : 0.30) ( 1 H 2s : 0.45) ( 2 H 1s : 0.30)
( 2 H 2s : 0.45)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: -0.13570 a.u.
( 1 H 1s : -0.32) ( 1 H 2s : -0.46) ( 2 H 1s : 0.32)
( 2 H 2s : 0.46)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.59106 a.u.
( 1 H 1s : -0.91) ( 1 H 2s : 0.80) ( 2 H 1s : -0.91)
( 2 H 2s : 0.80)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.67000 a.u.
( 1 H 1s : 0.94) ( 1 H 2s : -0.92) ( 2 H 1s : -0.94)
( 2 H 2s : 0.92)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.46240 a.u.
( 1 H 1p+1: 0.71) ( 2 H 1p+1: 0.71)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.46240 a.u.
( 1 H 1p-1: 0.71) ( 2 H 1p-1: 0.71)
Ground State Dipole Moment
----------------------------
X : -0.000000 a.u. -0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : -0.000000 a.u. -0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : -0.00 (multiplicity = 1.0 )
- Energy : -0.9987710557416053
- Natural orbitals
1.07018 0.92981 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.1322943027 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -0.782139853255 a.u. Time: 0.02 sec.
* Info * Overlap matrix computed in 0.01 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.01 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -0.782197472959 0.0000000000 0.00214017 0.00075393 0.00000000
2 -0.782198208427 -0.0000007355 0.00000243 0.00000086 0.00053449
3 -0.782198208429 -0.0000000000 0.00000007 0.00000003 0.00000212
*** SCF converged in 3 iterations. Time: 0.05 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -0.7821982084 a.u.
Electronic Energy : -0.9144925112 a.u.
Nuclear Repulsion Energy : 0.1322943027 a.u.
------------------------------------
Gradient Norm : 0.0000000737 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.28967 a.u.
( 1 H 1s : 0.30) ( 1 H 2s : 0.46) ( 2 H 1s : 0.30)
( 2 H 2s : 0.46)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: -0.14868 a.u.
( 1 H 1s : 0.32) ( 1 H 2s : 0.46) ( 2 H 1s : -0.32)
( 2 H 2s : -0.46)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.60146 a.u.
( 1 H 1s : 0.91) ( 1 H 2s : -0.83) ( 2 H 1s : 0.91)
( 2 H 2s : -0.83)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.65436 a.u.
( 1 H 1s : -0.93) ( 1 H 2s : 0.88) ( 2 H 1s : 0.93)
( 2 H 2s : -0.88)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.46422 a.u.
( 1 H 1p0 : 0.71) ( 2 H 1p0 : -0.71)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.46443 a.u.
( 1 H 1p+1: 0.66) ( 1 H 1p-1: -0.26) ( 2 H 1p+1: 0.66)
( 2 H 1p-1: -0.26)
Ground State Dipole Moment
----------------------------
X : 0.000000 a.u. 0.000000 Debye
Y : -0.000000 a.u. -0.000000 Debye
Z : -0.000000 a.u. -0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : 0.00 (multiplicity = 1.0 )
- Energy : -0.9986061861413567
- Natural orbitals
1.03178 0.96822 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.1175949358 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -0.770215954354 a.u. Time: 0.02 sec.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -0.770245883432 0.0000000000 0.00149442 0.00052645 0.00000000
2 -0.770246242419 -0.0000003590 0.00000095 0.00000034 0.00037504
*** SCF converged in 2 iterations. Time: 0.02 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -0.7702462424 a.u.
Electronic Energy : -0.8878411782 a.u.
Nuclear Repulsion Energy : 0.1175949358 a.u.
------------------------------------
Gradient Norm : 0.0000009489 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.27925 a.u.
( 1 H 1s : 0.31) ( 1 H 2s : 0.46) ( 2 H 1s : 0.31)
( 2 H 2s : 0.46)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: -0.15756 a.u.
( 1 H 1s : -0.31) ( 1 H 2s : -0.46) ( 2 H 1s : 0.31)
( 2 H 2s : 0.46)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.61352 a.u.
( 1 H 1s : 0.92) ( 1 H 2s : -0.84) ( 2 H 1s : 0.92)
( 2 H 2s : -0.84)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.64164 a.u.
( 1 H 1s : 0.92) ( 1 H 2s : -0.87) ( 2 H 1s : -0.92)
( 2 H 2s : 0.87)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.46607 a.u.
( 1 H 1p0 : 0.71) ( 2 H 1p0 : -0.71)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.46635 a.u.
( 1 H 1p+1: 0.70) ( 2 H 1p+1: 0.70)
Ground State Dipole Moment
----------------------------
X : 0.000000 a.u. 0.000000 Debye
Y : 0.000000 a.u. 0.000000 Debye
Z : 0.000000 a.u. 0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : 0.00 (multiplicity = 1.0 )
- Energy : -0.9985690947312191
- Natural orbitals
1.01383 0.98617 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
Self Consistent Field Driver Setup
====================================
Wave Function Model : Spin-Restricted Hartree-Fock
Initial Guess Model : Superposition of Atomic Densities
Convergence Accelerator : Two Level Direct Inversion of Iterative Subspace
Max. Number of Iterations : 50
Max. Number of Error Vectors : 10
Convergence Threshold : 1.0e-06
ERI Screening Threshold : 1.0e-12
Linear Dependence Threshold : 1.0e-06
* Info * Nuclear repulsion energy: 0.1058354422 a.u.
* Info * Overlap matrix computed in 0.00 sec.
* Info * Kinetic energy matrix computed in 0.00 sec.
* Info * Nuclear potential matrix computed in 0.00 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
* Info * Starting Reduced Basis SCF calculation...
* Info * ...done. SCF energy in reduced basis set: -0.762026297205 a.u. Time: 0.04 sec.
* Info * Overlap matrix computed in 0.01 sec.
* Info * Kinetic energy matrix computed in 0.01 sec.
* Info * Nuclear potential matrix computed in 0.01 sec.
* Info * Orthogonalization matrix computed in 0.00 sec.
Iter. | Hartree-Fock Energy | Energy Change | Gradient Norm | Max. Gradient | Density Change
--------------------------------------------------------------------------------------------
1 -0.762044195778 0.0000000000 0.00112546 0.00039650 0.00000000
2 -0.762044399480 -0.0000002037 0.00000038 0.00000013 0.00028293
*** SCF converged in 2 iterations. Time: 0.05 sec.
Spin-Restricted Hartree-Fock:
-----------------------------
Total Energy : -0.7620443995 a.u.
Electronic Energy : -0.8678798417 a.u.
Nuclear Repulsion Energy : 0.1058354422 a.u.
------------------------------------
Gradient Norm : 0.0000003765 a.u.
Ground State Information
------------------------
Charge of Molecule : 0.0
Multiplicity (2S+1) : 1
Magnetic Quantum Number (M_S) : 0.0
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -0.27168 a.u.
( 1 H 1s : 0.31) ( 1 H 2s : 0.46) ( 2 H 1s : 0.31)
( 2 H 2s : 0.46)
Molecular Orbital No. 2:
--------------------------
Occupation: 0.000 Energy: -0.16408 a.u.
( 1 H 1s : 0.31) ( 1 H 2s : 0.46) ( 2 H 1s : -0.31)
( 2 H 2s : -0.46)
Molecular Orbital No. 3:
--------------------------
Occupation: 0.000 Energy: 0.62166 a.u.
( 1 H 1s : -0.92) ( 1 H 2s : 0.85) ( 2 H 1s : -0.92)
( 2 H 2s : 0.85)
Molecular Orbital No. 4:
--------------------------
Occupation: 0.000 Energy: 0.63412 a.u.
( 1 H 1s : 0.92) ( 1 H 2s : -0.86) ( 2 H 1s : -0.92)
( 2 H 2s : 0.86)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.000 Energy: 1.46736 a.u.
( 1 H 1p0 : 0.71) ( 2 H 1p0 : -0.71)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.000 Energy: 1.46760 a.u.
( 1 H 1p+1: -0.27) ( 1 H 1p-1: 0.65) ( 2 H 1p+1: -0.27)
( 2 H 1p-1: 0.65)
Ground State Dipole Moment
----------------------------
X : -0.000000 a.u. -0.000000 Debye
Y : 0.000000 a.u. 0.000000 Debye
Z : 0.000000 a.u. 0.000000 Debye
Total : 0.000000 a.u. 0.000000 Debye
Configuration Interaction Driver
==================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 0
Number of active orbitals: 10
Number of virtual orbitals: 0
This is a CASSCF wavefunction: CAS(2,10)
CI expansion:
-------------
Number of determinants: 100
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Solved by explicit diagonalization
** Convergence reached in 0 iterations
Final results
-------------
* State 1
- S^2 : -0.00 (multiplicity = 1.0 )
- Energy : -0.9985598772014671
- Natural orbitals
1.00568 0.99432 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000
plt.figure(figsize=(6,4))
plt.title('Energy during H$_2$ dissociation')
x = np.array(distlist)
y = np.array(E_FCI)
z = np.array(E_hf)
plt.axhline(y = 0, color = "k", linestyle = 'dashed')
plt.plot(x,y, label='FullCI')
plt.plot(x,z, label='RHF')
plt.arrow(x[3], z[3], 0, y[3]-z[3], length_includes_head = True, head_width = 0.05, head_length = 0.01, facecolor='r', edgecolor='r')
plt.annotate('Dynamical correlation', xy = (x[3]+0.1, 0.5*(z[3]+y[3])), color = 'r')
plt.arrow(x[-1], z[-1], 0, y[-1]-z[-1], length_includes_head = True, head_width = 0.05, head_length = 0.01, facecolor='r', edgecolor='r')
plt.annotate('Static correlation', xy = (x[-1]-1.25, 0.5*(z[-1]+y[-1])), color = 'r')
plt.xlabel("H-H distance (Å)")
plt.ylabel("Energy (Ha)")
plt.legend()
plt.tight_layout(); plt.show()
Following Löwdin, we usually call correlation the difference between HF and FullCI. At short distance, this corresponds to our intuitive idea of electrons dynamically avoiding each other. However, at long distance, we expect one electron on each hydrogen atom, and there is no dynamical correlation. The error of Hartree-Fock is that even if one electron is on the first hydrogen, the second one is still equally shared between both hydrogens. We can visualize this using the functions written in “Electron correlation” chapter:
n=100
norb = space.n_active
origin = np.zeros((n, 3))
coords = np.zeros((n, 3))
def get_orbital_values(coords, molecule, basis, mol_orbs):
orbital_values = []
for i in range(norb):
values = np.array(vis_drv.get_mo(coords, molecule, basis, mol_orbs.alpha_to_numpy(), i))
orbital_values.append(values)
return orbital_values
def get_n12(ci_twopart_denmat, orbital_values_1, orbital_values_2):
hf_twopart_denmat = np.zeros((norb,norb,norb,norb))
hf_twopart_denmat[0,0,0,0] = 2.0 # first MO is doubly occupied
# Two-particle densities for FCI
n12_hf = np.zeros(n)
n12_fci = np.zeros(n)
for i in range(norb):
for j in range(norb):
for k in range(norb):
for l in range(norb):
n12_hf += (
hf_twopart_denmat[i, j, k, l]
* orbital_values_1[i]
* orbital_values_1[j]
* orbital_values_2[k]
* orbital_values_2[l]
)
n12_fci += (
ci_twopart_denmat[i, j, k, l]
* orbital_values_1[i]
* orbital_values_1[j]
* orbital_values_2[k]
* orbital_values_2[l]
)
return n12_hf, n12_fci
z = np.linspace(-1, 2, n)
coords[:, 2] = z
vis_drv = vlx.VisualizationDriver()
orbital_values = get_orbital_values(coords/bohr_to_Å, eq_molecule, basis, eq_orbitals)
orbital_values_at_H1 = get_orbital_values([[0, 0, 0]], eq_molecule, basis, eq_orbitals)
n12_hf, n12_fci = get_n12(eq_2density, orbital_values_at_H1, orbital_values)
plt.figure(figsize=(6,4))
plt.title('Pair density at equilibrium distance')
plt.plot(z,n12_hf, label='HF')
plt.plot(z,n12_fci, label='FullCI')
plt.axvline(x = 0, color = 'k', linestyle = ':')
plt.xlabel("Internuclear axis coordinate (Å)")
plt.ylabel("Two-particle density (a.u.)")
plt.legend()
plt.tight_layout(); plt.show()
z = np.linspace(-2, 8, n)
coords[:, 2] = z
orbital_values = get_orbital_values(coords/bohr_to_Å, molecule, basis, scf_drv.mol_orbs)
orbital_values_at_H1 = get_orbital_values([[0, 0, 0]], molecule, basis, scf_drv.mol_orbs)
ci_twopart_denmat = ci_drv.get_active_2body_density(0)
n12_hf, n12_fci = get_n12(ci_twopart_denmat, orbital_values_at_H1, orbital_values)
plt.figure(figsize=(6,4))
plt.title('Pair density at dissociation')
plt.plot(z,n12_hf, label='HF')
plt.plot(z,n12_fci, label='FullCI')
plt.axvline(x = 0, color = 'k', linestyle = ':')
plt.xlabel("Internuclear axis coordinate (Å)")
plt.ylabel("Two-particle density (a.u.)")
plt.legend()
plt.tight_layout(); plt.show()
This correlation is called by contrast “static” or “non-dynamical” or sometimes simply “strong” correlation. In practice both type of correlation have the same physical origin. If one electron is on one atom, the second electron is more likely to be on the second one. However, at long distances, the presence of one electron on one nucleus guarantees that the second electron is on the other one. Both types of correlation are also treated with the same excited determinants. But the distinction is useful as they may guide the way we calculate such systems.
While H
At equilibrium distance#
While true bond breakings are rare in chemistry, it turns out that static/strong correlation can occur even at equilibrium distance. One way to easily detect such thing is by looking at occupation numbers in a CI (or MCSCF) calculations. As explained in the MCSCF chapter, correlation will show as a difference in occupation number from the integers we expect. For standard dynamical correlation, the occupation numbers will typically be around 1.98 for occupied orbitals and 0.02 for virtual orbitals, showing that the CI expansion contains some configurations exciting from the occupied to the virtual, but that the dominant configuration is still the Hartree-Fock reference.
However, in some cases, we get higher deviations from 2 and 0, for example let’s look at ethene, including the carbon-carbon
mol_str = """
H 0.000000 0.923274 1.238289
H 0.000000 -0.923274 1.238289
H 0.000000 0.923274 -1.238289
H 0.000000 -0.923274 -1.238289
C 0.000000 0.000000 0.668188
C 0.000000 0.000000 -0.668188
"""
molecule = vlx.Molecule.read_str(mol_str)
molecule.set_charge(0)
molecule.set_multiplicity(1)
basis = vlx.MolecularBasis.read(molecule,"def2-sv(p)")
eht = mtp.OrbitalGuess()
guessorb = eht.compute(molecule,basis)
space = mtp.OrbSpace(molecule,guessorb)
space.cas_orbitals([5,6,8,13])
mcscfdrv = mtp.McscfDriver()
mcresults = mcscfdrv.compute(molecule,basis,space)
Multi-Configurational Self-Consistent Field Driver
====================================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 6
Number of active orbitals: 4
Number of virtual orbitals: 26
This is a CASSCF wavefunction: CAS(4,4)
CI expansion:
-------------
Number of determinants: 36
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Number of states : 1
State-specific calculation
- State of interest : 1
Max. iterations : 50
BFGS window : 5
Convergence thresholds:
- Energy change : 1e-08
- Gradient sq. norm : 1e-08
MCSCF Iterations
-------------------
Iter. | Average Energy | E. Change | Grad. Norm | CI Iter. | Trust rad. | Time
---------------------------------------------------------------------------------
1 -77.869303461 0.0e+00 4.3e-01 0 0.40 0:00:00
2 -77.989696943 -1.2e-01 6.1e-02 0 0.40 0:00:00
3 -78.017442330 -2.8e-02 3.4e-03 0 0.48 0:00:00
4 -78.019216051 -1.8e-03 1.1e-04 0 0.58 0:00:00
5 -78.019343633 -1.3e-04 1.5e-05 0 0.58 0:00:00
6 -78.019364156 -2.1e-05 2.2e-06 0 0.69 0:00:00
7 -78.019366204 -2.0e-06 4.2e-07 0 0.80 0:00:00
8 -78.019366517 -3.1e-07 3.2e-08 0 0.80 0:00:00
9 -78.019366537 -2.0e-08 1.7e-09 0 0.80 0:00:00
10 -78.019366538 -1.2e-09 6.5e-11 0 0.80 0:00:00
** Convergence reached in 10 iterations
11 -78.019366538 -5.8e-11 6.2e-12 0 0.80 0:00:00
Final results
-------------
* State 1
- S^2 : 0.00 (multiplicity = 1.0 )
- Energy : -78.01936653836086
- Natural orbitals
1.98365 1.92629 0.07350 0.01656
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 4:
--------------------------
Occupation: 2.000 Energy: -0.77861 a.u.
( 1 H 1s : -0.19) ( 2 H 1s : -0.19) ( 3 H 1s : 0.19)
( 4 H 1s : 0.19) ( 5 C 1s : 0.15) ( 5 C 2s : -0.32)
( 5 C 1p0 : -0.17) ( 6 C 1s : -0.15) ( 6 C 2s : 0.32)
( 6 C 1p0 : -0.17)
Molecular Orbital No. 5:
--------------------------
Occupation: 2.000 Energy: -0.63704 a.u.
( 1 H 1s : -0.20) ( 2 H 1s : 0.20) ( 3 H 1s : -0.20)
( 4 H 1s : 0.20) ( 5 C 1p-1: -0.32) ( 6 C 1p-1: -0.32)
Molecular Orbital No. 6:
--------------------------
Occupation: 2.000 Energy: -0.50339 a.u.
( 1 H 1s : -0.25) ( 2 H 1s : 0.25) ( 3 H 1s : 0.25)
( 4 H 1s : -0.25) ( 5 C 1p-1: -0.30) ( 5 C 2p-1: -0.16)
( 6 C 1p-1: 0.30) ( 6 C 2p-1: 0.16)
Molecular Orbital No. 7:
--------------------------
Occupation: 1.984 Energy: -0.74951 a.u.
( 5 C 2s : -0.28) ( 5 C 1p0 : 0.39) ( 6 C 2s : -0.28)
( 6 C 1p0 : -0.39)
Molecular Orbital No. 8:
--------------------------
Occupation: 1.926 Energy: -0.36398 a.u.
( 5 C 1p+1: 0.40) ( 5 C 2p+1: 0.29) ( 6 C 1p+1: 0.40)
( 6 C 2p+1: 0.29)
Molecular Orbital No. 9:
--------------------------
Occupation: 0.074 Energy: 0.26140 a.u.
( 5 C 1p+1: 0.63) ( 5 C 2p+1: 0.31) ( 6 C 1p+1: -0.63)
( 6 C 2p+1: -0.31)
Molecular Orbital No. 10:
--------------------------
Occupation: 0.017 Energy: 0.96850 a.u.
( 5 C 1s : -0.22) ( 5 C 2s : 0.91) ( 5 C 1p0 : -0.84)
( 5 C 2p0 : -0.24) ( 6 C 1s : 0.22) ( 6 C 2s : -0.91)
( 6 C 1p0 : -0.84) ( 6 C 2p0 : -0.24)
Molecular Orbital No. 11:
--------------------------
Occupation: 0.000 Energy: 0.20754 a.u.
( 1 H 2s : -1.12) ( 2 H 2s : -1.12) ( 3 H 2s : -1.12)
( 4 H 2s : -1.12) ( 5 C 3s : 1.60) ( 5 C 2p0 : 0.57)
( 6 C 3s : 1.60) ( 6 C 2p0 : -0.57)
Molecular Orbital No. 12:
--------------------------
Occupation: 0.000 Energy: 0.23435 a.u.
( 1 H 2s : -1.16) ( 2 H 2s : 1.16) ( 3 H 2s : -1.16)
( 4 H 2s : 1.16) ( 5 C 1p-1: 0.22) ( 5 C 2p-1: 0.76)
( 6 C 1p-1: 0.22) ( 6 C 2p-1: 0.76)
Molecular Orbital No. 13:
--------------------------
Occupation: 0.000 Energy: 0.25745 a.u.
( 1 H 2s : 1.29) ( 2 H 2s : 1.29) ( 3 H 2s : -1.29)
( 4 H 2s : -1.29) ( 5 C 3s : -1.52) ( 5 C 1p0 : -0.18)
( 5 C 2p0 : -0.70) ( 6 C 3s : 1.52) ( 6 C 1p0 : -0.18)
( 6 C 2p0 : -0.70)
viewer = mtp.OrbitalViewer()
viewer.plot(molecule, basis, space)
print(mcresults["natural_occupations"][0])
[1.98365467 1.92628683 0.07349913 0.01655936]
The first and last orbitals are the
from IPython.display import Image
Image(filename="../../img/tutorials/ethene_bond_breaking.png", width=400)
The reason for this is that at the equilibrium distance, the
coords = np.zeros((n, 3))
z = np.linspace(-1.5, 1.5, n)
coords[:, 2] = z
coords_shifted = np.zeros((n, 3))
coords_shifted[:, 0] = 0.35 # Pi orbitals are 0 in the xy plane, so we compute them at shifted x value
coords_shifted[:, 2] = z
sigma = np.array(vis_drv.get_mo(coords/bohr_to_Å, molecule, basis, space.molecular_orbitals.alpha_to_numpy(), 6))
pi = np.array(vis_drv.get_mo(coords_shifted/bohr_to_Å, molecule, basis, space.molecular_orbitals.alpha_to_numpy(), 7))
# Make sure the sign is the same in the middle of the bond
sigma*= np.sign(sigma[50])
pi*= np.sign(pi[50])
sigma_max = abs(z[np.argmax(sigma)])
pi_max = abs(z[np.argmax(pi)])
plt.figure(figsize=(6,4))
plt.title('Absolute value of the orbital profiles')
plt.plot(z,sigma, label='$\sigma$')
plt.plot(z,pi, label='$\pi$')
plt.xlabel("Internuclear axis coordinate (Å)")
plt.ylabel("Two-particle density (a.u.)")
plt.axhline(y = 0, color = "k")
plt.axvline(x = 0.668188, color = 'k', linestyle = ':')
plt.axvline(x = -0.668188, color = 'k', linestyle = ':')
plt.legend()
plt.tight_layout(); plt.show()
<>:21: SyntaxWarning: invalid escape sequence '\s'
<>:22: SyntaxWarning: invalid escape sequence '\p'
<>:21: SyntaxWarning: invalid escape sequence '\s'
<>:22: SyntaxWarning: invalid escape sequence '\p'
/tmp/ipykernel_194837/1833864943.py:21: SyntaxWarning: invalid escape sequence '\s'
plt.plot(z,sigma, label='$\sigma$')
/tmp/ipykernel_194837/1833864943.py:22: SyntaxWarning: invalid escape sequence '\p'
plt.plot(z,pi, label='$\pi$')
It is clear in this orbital profile that the distance between the maxima of the
This situation is very general, and whenever bonds get “frustrated” due to the electronic repulsion of other bonds, they may be stretched and display static correlation even in the equilibrium structure. While the
mol_str = """
B 0.0 0.0 -0.80494
B 0.0 0.0 0.80494
"""
molecule = vlx.Molecule.read_str(mol_str)
molecule.set_charge(0)
molecule.set_multiplicity(3)
basis = vlx.MolecularBasis.read(molecule,"def2-sv(p)")
eht = mtp.OrbitalGuess()
guessorb = eht.compute(molecule,basis)
space = mtp.OrbSpace(molecule,guessorb)
space.cas(6,8)
mcscfdrv = mtp.McscfDriver()
mcresults = mcscfdrv.compute(molecule,basis,space)
Multi-Configurational Self-Consistent Field Driver
====================================================
Active space definition:
------------------------
Number of inactive (occupied) orbitals: 2
Number of active orbitals: 8
Number of virtual orbitals: 18
This is a CASSCF wavefunction: CAS(6,8)
CI expansion:
-------------
Number of determinants: 1960
╭────────────────────────────────────╮
│ Driver settings │
╰────────────────────────────────────╯
Number of states : 1
State-specific calculation
- State of interest : 1
Max. iterations : 50
BFGS window : 5
Convergence thresholds:
- Energy change : 1e-08
- Gradient sq. norm : 1e-08
MCSCF Iterations
-------------------
Iter. | Average Energy | E. Change | Grad. Norm | CI Iter. | Trust rad. | Time
---------------------------------------------------------------------------------
1 -49.135597057 0.0e+00 6.3e-02 10 0.40 0:00:00
2 -49.157696418 -2.2e-02 5.1e-03 10 0.40 0:00:00
3 -49.160489940 -2.8e-03 1.9e-04 10 0.48 0:00:00
4 -49.160601386 -1.1e-04 6.1e-06 10 0.58 0:00:00
5 -49.160607928 -6.5e-06 4.5e-07 9 0.58 0:00:00
6 -49.160608619 -6.9e-07 4.7e-08 9 0.69 0:00:00
7 -49.160608651 -3.3e-08 2.1e-09 9 0.80 0:00:00
8 -49.160608654 -2.4e-09 1.3e-10 9 0.80 0:00:00
** Convergence reached in 8 iterations
9 -49.160608654 -2.0e-10 1.7e-11 9 0.80 0:00:00
Final results
-------------
* State 1
- S^2 : 2.00 (multiplicity = 3.0 )
- Energy : -49.16060865387955
- Natural orbitals
1.96302 1.64311 0.96618 0.96618 0.32575 0.05890 0.05890 0.01795
Spin Restricted Orbitals
------------------------
Molecular Orbital No. 1:
--------------------------
Occupation: 2.000 Energy: -7.65972 a.u.
( 1 B 1s : 0.70) ( 2 B 1s : 0.71)
Molecular Orbital No. 2:
--------------------------
Occupation: 2.000 Energy: -7.65969 a.u.
( 1 B 1s : -0.71) ( 2 B 1s : 0.70)
Molecular Orbital No. 3:
--------------------------
Occupation: 1.963 Energy: -0.66298 a.u.
( 1 B 1s : 0.19) ( 1 B 2s : 0.43) ( 1 B 1p0 : -0.19)
( 2 B 1s : 0.19) ( 2 B 2s : 0.43) ( 2 B 1p0 : 0.19)
Molecular Orbital No. 4:
--------------------------
Occupation: 1.643 Energy: -0.29492 a.u.
( 1 B 1s : -0.18) ( 1 B 2s : -0.46) ( 1 B 3s : -0.27)
( 1 B 1p0 : -0.21) ( 2 B 1s : 0.18) ( 2 B 2s : 0.46)
( 2 B 3s : 0.27) ( 2 B 1p0 : -0.21)
Molecular Orbital No. 5:
--------------------------
Occupation: 0.966 Energy: -0.13452 a.u.
( 1 B 1p+1: -0.43) ( 1 B 2p+1: 0.25) ( 2 B 1p+1: -0.43)
( 2 B 2p+1: 0.25)
Molecular Orbital No. 6:
--------------------------
Occupation: 0.966 Energy: -0.13452 a.u.
( 1 B 1p-1: -0.43) ( 1 B 2p-1: 0.25) ( 2 B 1p-1: -0.43)
( 2 B 2p-1: 0.25)
Molecular Orbital No. 7:
--------------------------
Occupation: 0.326 Energy: -0.04079 a.u.
( 1 B 2s : -0.28) ( 1 B 1p0 : -0.49) ( 1 B 2p0 : 0.21)
( 2 B 2s : -0.28) ( 2 B 1p0 : 0.49) ( 2 B 2p0 : -0.21)
Molecular Orbital No. 8:
--------------------------
Occupation: 0.059 Energy: 0.21444 a.u.
( 1 B 1p-1: 0.63) ( 1 B 2p-1: -0.32) ( 2 B 1p-1: -0.63)
( 2 B 2p-1: 0.32)
Molecular Orbital No. 9:
--------------------------
Occupation: 0.059 Energy: 0.21444 a.u.
( 1 B 1p+1: 0.63) ( 1 B 2p+1: -0.32) ( 2 B 1p+1: -0.63)
( 2 B 2p+1: 0.32)
Molecular Orbital No. 10:
--------------------------
Occupation: 0.018 Energy: 0.63202 a.u.
( 1 B 1s : 0.19) ( 1 B 2s : 0.94) ( 1 B 3s : 0.26)
( 1 B 1p0 : -0.78) ( 1 B 2p0 : 0.42) ( 2 B 1s : -0.19)
( 2 B 2s : -0.94) ( 2 B 3s : -0.26) ( 2 B 1p0 : -0.78)
( 2 B 2p0 : 0.42)
Molecular Orbital No. 11:
--------------------------
Occupation: 0.000 Energy: 0.32558 a.u.
( 1 B 1s : -0.17) ( 1 B 2s : -1.17) ( 1 B 3s : 1.00)
( 1 B 2p0 : 0.20) ( 2 B 1s : -0.17) ( 2 B 2s : -1.17)
( 2 B 3s : 1.00) ( 2 B 2p0 : -0.20)
Molecular Orbital No. 12:
--------------------------
Occupation: 0.000 Energy: 0.40442 a.u.
( 1 B 3s : 0.26) ( 1 B 1p0 : -0.53) ( 1 B 2p0 : -1.00)
( 2 B 3s : 0.26) ( 2 B 1p0 : 0.53) ( 2 B 2p0 : 1.00)
print(mcresults["natural_occupations"][0])
[1.96302192 1.64311455 0.96617833 0.96617832 0.32574991 0.05890328
0.05890325 0.01795044]