Intermediate state representation#

Derivation via intermediate states#

As the name suggest, the intermediate state representation (ISR) approach consists of constructing the ADC matrix with the help of intermediate states \(\ket{\tilde{\Psi}_I}\), obtained by applying excitation operators to the ground state \(\ket{0}\). In second quantization, the excitation operator is written as \(\hat{C}_I=\{ \hat{a}_a^\dagger\hat{a}_i;\hat{a}_a^\dagger\hat{a}_b^\dagger\hat{a}_j\hat{a}_i, a<b, i<j;... \}\), where the indices \(a,b...\) refer to unoccupied orbitals, while \(i,j...\) represent occupied orbitals [MS96, Sch91, ST04]. Schematic representations of single and double excitations, which are the only two excitation classes that are needed for ADC up to third order in perturbation theory, are depicted in below figure:

../../_images/isr_adc_matrix.svg

Fig. 30 Illustration of (a) single excitations, (b) double excitations, (c) the structure of the ADC(2) matrix. The numbers in parenthesis indicate the highest order of perturbation theory used to describe each particular block.#

The intermediate states \(\ket{\tilde{\Psi}_I}\) are obtained by first applying \(\hat{C}_I\) to the many-body ground state:

\[ \ket{\Psi_I^{0}}=\hat{C}_I\ket{\Psi_0} \, , \]

and then performing a Gram–Schmidt orthogonalization procedure with respect to lower excitation classes (including the ground state) to obtain precursor states \(\ket{\Psi_I^{\#}}\), which can then be orthonormalized symmetrically according to [Wen16]:

\[ %:label: eq:ISdefinition \ket{\tilde{\Psi}_I}=\sum_J\ket{\Psi_J^{\#}}S_{IJ}^{-1/2}\, , \]

where \(S_{IJ}=\braket{\Psi_I^{\#}|\Psi_J^{\#}}\) are overlap integrals of the precursor states.

The elements of the ADC matrix \(\mathbf{M}\) are obtained as matrix elements of the shifted Hamiltonian in the basis of the intermediate states:

\[ %:label: eq:Mdef M_{IJ}=\bra{\tilde{\Psi}_I}\hat{H}-E_0\ket{\tilde{\Psi}_J}=\sum_{K,L}S_{IK}^{-1/2}\bra{ \Psi_K^{\#}}\hat{H}-E_0\ket{\Psi_L^{\#}}S_{LJ}^{-1/2} \, , \]

where \(E_0\) is the ground-state energy. This representation of the (shifted) Hamiltonian leads to a Hermitian eigenvalue equation,

\[ %:label: eq:adc_eigenvalue_eq \mathbf{MY} = \mathbf{Y\Omega} \, , \quad \mathbf{Y}^\dagger \mathbf{Y} = \mathbf{1} \, , \]

the solution of which yields vertical excitation energies (\(\omega_n = E_n - E_0\)) as eigenvalues, collected in the diagonal matrix \(\mathbf{\Omega}\), and the corresponding excitation vectors as eigenvectors \(\mathbf{Y}_n\), collected in the columns of \(\mathbf{Y}\).

Having obtained an expression for the ADC matrix, we return to the series expansion of the polarization propagator. In the same way as the propagator is expanded in series, also the matrix elements can be written in terms of orders of perturbation [Wen16]:

\[ %:label: eq:Mseries M_{IJ}^{(k+l+m)} = \sum_{K,L}\left(S_{IK}^{-1/2}\right)^{(k)}\left(\bra{\Psi_K^{\#}}\hat{H}-E_0\ket{\Psi_L^{\#}}\right)^{(l)}\left(S_{LJ}^{-1/2}\right)^{(m)}\, , \]

where \(k\) and \(m\) are the orders of perturbation theory used for the overlap matrices \(S_{IK}\) and \(S_{LJ}\), \(l\) is the order used for the matrix elements of the shifted Hamiltonian in the basis of precursor states, and the sum \(k + l + m = n\) represents the order of the contribution to the ADC matrix \(\mathbf{M}\). In order to get all contributions of a given order \(n\), one needs to sum \(k,l,m\) over all terms for which \(k+l+m = n\).

The effective transition amplitudes \(\mathbf{f}\) can analogously be obtained in the ISR via

\[ f_{I,pq} = \bra{\tilde{\Psi}_{I}} \hat{a}_p^\dagger \hat{a}_q \ket{\Psi_0} \]

and the transition density matrices or transition amplitudes \(\mathbf{x}\) by contracting with the eigenvectors, \(\mathbf{x} = \mathbf{Y}^\dagger \mathbf{f}\). Note that since \(\mathbf{f}\) is a transition amplitude and not an expectation value, it is not symmetric: \(f_{I,pq} \neq f_{I,qp}\).

Structure of the ADC matrix#

As is shown here, the ADC eigenvalue equation in the space of single (S) and double (D) excitations can be written as

\[ \mathbf{M}_{\text{eff}} \mathbf{X}_{\text{S}} = \big[ \mathbf{M}_{\text{SS}} + \mathbf{M}_{\text{SD}} (\omega \mathbf{1} - \mathbf{M}_{\text{DD}})^{-1} \mathbf{M}_{\text{DS}} \big] \mathbf{X}_{\text{S}} = \omega \mathbf{X}_{\text{S}} \, , \]

where the subscript denote the corresponding subblocks of the ADC matrix and \(\mathbf{X}_{\text{S}}\) is the singles part of the eigenvector. From this, the order structure of the ADC matrix can be explained as follows. The leading contributions of \(\mathbf{M}_{\text{SS}}\) and \(\mathbf{M}_{\text{DD}}\), and thus also of \((\omega \mathbf{1} - \mathbf{M}_{\text{DD}})^{-1}\), are of zeroth order, whereas the leading contributions of \(\mathbf{M}_{\text{SD}}\) and \(\mathbf{M}_{\text{DS}}\) are of first order. Hence, if excited states dominated by single excitations are desired through, say, second order, then \(\mathbf{M}_{\text{eff}}\) needs to be correct through this order, meaning that \(\mathbf{M}_{\text{SS}}\) needs to be correct through second order, \(\mathbf{M}_{\text{SD}}\) and \(\mathbf{M}_{\text{DS}}\) need to be correct through first order, whereas \(\mathbf{M}_{\text{DD}}\) only needs to be correct through zeroth order. Expanding \(\mathbf{M}_{\text{DD}}\) through first order would lead to a third-order term since the coupling blocks are at least of first order. Equivalently, expanding the coupling blocks through second order would lead to third-order terms, which are neglected in a second-order scheme. A consistent third-order method requires each block one order higher in perturbation theory, i.e. the singles-singles block through third oder, the coupling blocks through second order and the doubles-doubles block through first order. These findings are depicted schematically in Fig. 31 for ADC(\(n\)) schemes up to third order. ADC(2)-x represents an ad hoc extension of ADC(2), where the first-order terms in the doubles-doubles block from ADC(3) are taken into account in order to improve the description of doubly excited states. However, this generally leads to an imbalanced description of singly and doubly excited states [DW15].

../../_images/adc_matrices_0_to_3.svg

Fig. 31 Structure of the ADC(\(n\)) matrix at different orders \(n\). The numbers indicate the orders of perturbation theory used to expand the corresponding block.#

Explicit expressions for ADC(2)#

Using above expansion of \(M_{IJ}\) in combination with specific classes of excitation operators and truncating the series at the desired order, various levels of ADC theory are obtained. One aspect to note is that the excitation classes needed to construct a specific ADC level are directly connected to the order of perturbation theory. This can be easily seen in Fig. 29, where the zeroth and first order terms are related to single excitations (only one particle-hole pair is involved), while second order terms involve double excitations (two particle-hole pairs are involved). To illustrate this further, we list the explicit expressions for the ADC matrix \(\mathbf{M}\) in spin-orbitals up to second order [Wen16, Wor09]:

\[\begin{split} %:label: eq:adcmat_ph_ph_0 M_{ia,jb}^{(0)} &=(\varepsilon_a - \varepsilon_i) \delta_{ab} \delta_{ij} \\ %:label: eq:adcmat_ph_ph_1 M_{ia,jb}^{(1)} &= -\braket{ja||ib} \\ %:label: eq:adcmat_ph_ph_2 %M_{ia,jb}^{(2)} &= \frac{1}{4} \delta_{ij} \sum_{c,k,l} \left[ \frac{\braket{ac||kl} \braket{kl||bc}}{\epsilon_a + %\epsilon_c - \epsilon_k - \epsilon_l} + \frac{\braket{ac||kl} \braket{kl||bc}}{\epsilon_b + \epsilon_c - %\epsilon_k - \epsilon_l} \right] \nonumber \\ %&+\frac{1}{4} \delta_{ab} \sum_{c,d,k} \left[ \frac{\braket{cd||ik} \braket{jk||cd}}{\epsilon_c + \epsilon_d - %\epsilon_i - \epsilon_k} + \frac{{\braket{cd||ik} \braket{jk||cd}}}{\epsilon_c + \epsilon_d - \epsilon_j - %\epsilon_k} \right] \nonumber \\ %&-\frac{1}{2}\sum_{c,k} \left[ \frac{{\braket{ac||ik} \braket{jk||bc}}}{\epsilon_a + \epsilon_c - \epsilon_i - %\epsilon_k} + \frac{{\braket{ac||ik} \braket{jk||bc}}}{\epsilon_b + \epsilon_c - \epsilon_j - \epsilon_k} \right] %\\ M_{ia,jb}^{(2)} &= \frac{1}{4} \delta_{ij} \sum_{c,k,l} \left[ t_{ackl} \braket{kl||bc} + t_{bckl}^* \braket{ac||kl} \right] \nonumber \\ &+\frac{1}{4} \delta_{ab} \sum_{c,d,k} \left[ t_{cdik} \braket{jk||cd} + t_{jkcd}^* \braket{cd||ik} \right] \nonumber \\ &-\frac{1}{2}\sum_{c,k} \left[ t_{acik} \braket{jk||bc} + t_{jkbc}^* \braket{ac||ik} \right] \\ %:label: eq:adcmat_ph_2p2h_1 M_{ia,klcd}^{(1)} &= \braket{kl||id}\delta_{ac}-\braket{kl||ic}\delta_{ad}-\braket{al||cd}\delta_{ik}+{\braket{ak||cd}\delta_{il}}\label{eq:adcmat_ph_2p2h_1} \\ %:label: eq:c2p2hph_1 M_{ijab,kc}^{(1)} &= \braket{kb||ij}\delta_{ac}-\braket{ka||ij}\delta_{bc}-\braket{ab||cj}\delta_{ik}+{\braket{ab||ci}\delta_{jk}}\label{c2p2hph_1} \\ %:label: eq:adcmat_2p2h_2p2h_0 M_{ijab,klcd}^{(0)} &= (\varepsilon_a + \varepsilon_b - \varepsilon_i - \varepsilon_j) \delta_{ac} \delta_{bd}\delta_{ik}\delta_{jl}\,, \end{split}\]

where \(\varepsilon_p\) are Hartree–Fock (HF) orbital energies, \(\braket{pq||rs}\) are anti-symmetrized two-electron integrals in physicists’ notation [SO12], \(\delta_{pq}\) is the Kronecker delta, and the \(t\)-amplitudes are \(t_{ijab} = \frac{\langle ab || ij \rangle}{\varepsilon_{a} + \varepsilon_{b} - \varepsilon_{i} - \varepsilon_{j}}\).

Explicit expressions for the effective transition amplitudes \(\mathbf{f}\) can be derived in an analogous manner. The necessary blocks through second order in perturbation theory are given as

\[\begin{split} f_{jb,ai}^{(0)} &= \delta_{ij} \delta_{ab} \\ f_{jb,ia}^{(1)} &= - t_{ijab} \\ f_{ka,ij}^{(2)} &= \frac{\delta_{kj}}{\varepsilon_a - \varepsilon_i} \Big[ \frac12 \sum_{jbc} t_{ijbc} \braket{ja||bc} + \frac12 \sum_{jkb} t_{jkab} \braket{jk||ib} \Big] \\ f_{klab,ij}^{(1)} &= - \delta_{ki} t_{jlab} \\ f_{ic,ab}^{(2)} &= \frac{- \delta_{ac}}{\varepsilon_b - \varepsilon_i} \Big[ \frac12 \sum_{jcd} t_{ijcd} \braket{jb||cd} + \frac12 \sum_{jkc} t_{jkbc} \braket{jk||ic} \Big] \\ f_{ijcd,ab}^{(1)} &= \delta_{ac} t_{ijbd} \\ f_{jb,ai}^{(2)} &= \frac12 \sum_{kc} t_{ikac}^* t_{jkbc} - \frac12 \delta_{ij} \gamma_{ab}^{(2)} + \frac12 \delta_{ab} \gamma_{ij}^{(2)} \\ %- \delta_{ij} \frac14 \sum_{klc} t_{klac}^* t_{klbc} - \delta_{ab} \frac14 \sum_{kcd} t_{ikcd}^* t_{jkcd} \\ f_{jb,ia}^{(2)} &= - \frac{1}{\varepsilon_a + \varepsilon_b - \varepsilon_i - \varepsilon_j} \Big[ (1 - \hat{P}_{ij})(1 - \hat{P}_{ab}) \big[ \sum_{kc} t_{ikac} \braket{kb||jc} \big] \\ &- \frac12 \sum_{cd} t_{ijcd} \braket{ab||cd} - \frac12 \sum_{kl} t_{klab} \braket{kl||ij} \Big] \, , \end{split}\]

where the operator \(\hat{P}_{pq}\) permutes the indices \(p\) and \(q\) in the following expression, and the second-order corrections to the ground-state density matrix \(\gamma_{ab}^{(2)}\) and \(\gamma_{ij}^{(2)}\) were defined here.

The structure of the ADC(2) matrix is depicted in Fig. 30(c). In principle, the ADC(2) matrix contains all the possible single and double excitations which can be constructed for the system of interest (using a particular basis set). However, to calculate all these excitations would be computationally very expensive or practically impossible for all but the smallest of systems. In practice, therefore, only the lowest \(n\) excited states are ever calculated by means of iterative diagonalization algorithms, where \(n\) is the number of states requested by the user. This means that the space of valence excitations is easily accessible, but makes the space of core excitations impossible to reach, except for molecules with very few electrons. An approach to overcome this problem will be discussed in more detail in the next section.

Excited-state properties#

A distinct advantage of the ISR over the classical propagator approach is that it gives direct access to excited-state wave functions by expanding it in the intermediate-state basis as

\[ \ket{\Psi_n} = \sum_{J} Y_{Jn} \ket{\tilde{\Psi}_J} \, , \]

where the elements of the eigenvectors are the expansion coefficients, \(Y_{Jn} = \braket{\tilde{\Psi}_J | \Psi_n}\). This immediately offers the opportunity to calculate physical properties \(D_n\) of electronically excited state \(n\) via [ST04]

\[ D_n = \bra{\Psi_n} \hat{D} \ket{\Psi_n} = \mathbf{Y}_n^\dagger \, \tilde{\mathbf{D}} \, \mathbf{Y}_n = \sum_{IJ} Y_{In}^* \, \tilde{D}_{IJ} \, Y_{Jn} \, , \]

where \(\tilde{\mathbf{D}}\) is the representation of the operator \(\hat{D}\) corresponding to the observable in the intermediate state basis,

\[ \tilde{D}_{IJ} = \bra{\tilde{\Psi}_I} \hat{D} \ket{\tilde{\Psi}_J} \, . \]

The matrix \(\tilde{\mathbf{D}}\) has a perturbation expansion analogous to that of \(\mathbf{M}\) [ST04]. It should be noted that properties calculated in this manner generally differ from those calculated as derivatives of the energy \(E_n\) [HRDH19]. Transition moments between two different excited states (\(m \neq n\)) can be obtained in a completely analogous manner as

\[ T_{mn} = \bra{\Psi_m} \hat{D} \ket{\Psi_n} = \mathbf{Y}_m^\dagger \, \tilde{\mathbf{D}} \, \mathbf{Y}_n \, . \]

While explicit expressions for the ADC matrix \(\mathbf{M}\) are available through third order, the effective transition moments \(\mathbf{f}\) and property matrix \(\tilde{\mathbf{D}}\) are only available through second order. Combining the eigenvectors \(\mathbf{Y}\) (and eigenvalues) of the ADC(3) matrix together with the second-order descriptions of \(\mathbf{f}\) and \(\tilde{\mathbf{D}}\) yields transition moments, oscillator strengths and other properties at a level referred to as “ADC(3/2)” [ST04].

Size consistency#

As was discussed in previous sections, truncated CI schemes suffer from the size-consistency problem, while MP2 does not. For excited-state methods, we need to define size consistency somewhat differently. Imagine again two systems that are very far apart, such that they do not interact. For the ground-state energy, a size-consistent method has to give the same result for the composite system as the sum of the individual fragments. Concerning excitation energies, the result for one of the fragments should be the same when we apply the method to the individual fragment or to the composite system and look at the local excitations on the respective fragment (to be more specific, this is referred to as size intensivity). We will check this in the following for a system consisting of two small molecules for ADC(1) and ADC(2).

import veloxchem as vlx
import gator
from gator.adconedriver import AdcOneDriver
from gator.adctwodriver import AdcTwoDriver
import numpy as np

np.set_printoptions(precision=5, suppress=True)
* Warning * Environment variable OMP_NUM_THREADS not set.
* Warning * Setting OMP_NUM_THREADS to 8.
* Warning * Setting MKL_THREADING_LAYER to "GNU".
# LiH molecule
lih_xyz="""2

Li  0.000000   0.000000   0.000000
H   0.000000   0.000000   1.000000
"""
lih = vlx.Molecule.from_xyz_string(lih_xyz)

# Water molecule
h2o_xyz = """3
                                                                                                                        
O    0.000000000000        0.000000000000        0.000000000000                         
H    0.000000000000        0.740848095288        0.582094932012                         
H    0.000000000000       -0.740848095288        0.582094932012
"""
h2o = vlx.Molecule.from_xyz_string(h2o_xyz)

# Basis set
basis_set_label = "6-31g"
basis_lih = vlx.MolecularBasis.read(lih, basis_set_label)
basis_h2o = vlx.MolecularBasis.read(h2o, basis_set_label)

Calculate HF reference states#

# SCF will be run by VeloxChem through Gator
scf_lih = gator.run_scf(lih, basis_lih, conv_thresh=1e-10, verbose=False)
scf_h2o = gator.run_scf(h2o, basis_h2o, conv_thresh=1e-10, verbose=False)
e_scf_lih = scf_lih.get_scf_energy()
e_scf_h2o = scf_h2o.get_scf_energy()
                                                                                                                          
                                            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-10                                                              
                   ERI Screening Scheme            : Cauchy Schwarz + Density                                             
                   ERI Screening Mode              : Dynamic                                                              
                   ERI Screening Threshold         : 1.0e-12                                                              
                   Linear Dependence Threshold     : 1.0e-06                                                              
                                                                                                                          
* Info * Nuclear repulsion energy: 1.5875316327 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.03 sec.                                                                   
                                                                                                                          
* Info * SAD initial guess computed in 0.00 sec.                                                                          
                                                                                                                          
* Info * Starting Reduced Basis SCF calculation...                                                                        
* Info * ...done. SCF energy in reduced basis set: -7.846922313983 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        -7.867950663866    0.0000000000      0.06853478      0.01397032      0.00000000                
                  2        -7.871025271945   -0.0030746081      0.01619204      0.00435197      0.13927999                
                  3        -7.871354381385   -0.0003291094      0.00163331      0.00029700      0.03163457                
                  4        -7.871355542872   -0.0000011615      0.00024747      0.00004421      0.00227037                
                  5        -7.871355574879   -0.0000000320      0.00005485      0.00001197      0.00043636                
                  6        -7.871355578105   -0.0000000032      0.00000137      0.00000026      0.00013962                
                  7        -7.871355578106   -0.0000000000      0.00000012      0.00000002      0.00000122                
                  8        -7.871355578106    0.0000000000      0.00000001      0.00000000      0.00000017                
                  9        -7.871355578106   -0.0000000000      0.00000000      0.00000000      0.00000001                
                 10        -7.871355578106   -0.0000000000      0.00000000      0.00000000      0.00000000                
                                                                                                                          
               *** SCF converged in 10 iterations. Time: 0.04 sec.                                                        
                                                                                                                          
               Spin-Restricted Hartree-Fock:                                                                              
               -----------------------------                                                                              
               Total Energy                       :       -7.8713555781 a.u.                                              
               Electronic Energy                  :       -9.4588872108 a.u.                                              
               Nuclear Repulsion Energy           :        1.5875316327 a.u.                                              
               ------------------------------------                                                                       
               Gradient Norm                      :        0.0000000000 a.u.                                              
                                                                                                                          
                                                                                                                          
               Ground State Information                                                                                   
               ------------------------                                                                                   
               Charge of Molecule            :  0.0                                                                       
               Multiplicity (2S+1)           :  1.0                                                                       
               Magnetic Quantum Number (M_S) :  0.0                                                                       
                                                                                                                          
                                                                                                                          
                                                 Spin Restricted Orbitals                                                 
                                                 ------------------------                                                 
                                                                                                                          
               Molecular Orbital No.   1:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:   -2.44320 a.u.                                                                  
               (   1 Li  1s  :     1.00)                                                                                  
                                                                                                                          
               Molecular Orbital No.   2:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:   -0.32129 a.u.                                                                  
               (   1 Li  1s  :     0.23) (   1 Li  2s  :    -0.29) (   1 Li  1p0 :    -0.38)                              
               (   2 H   1s  :    -0.36) (   2 H   2s  :    -0.31)                                                        
                                                                                                                          
               Molecular Orbital No.   3:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.01856 a.u.                                                                  
               (   1 Li  2s  :    -0.16) (   1 Li  3s  :     0.96) (   1 Li  2p0 :    -0.50)                              
               (   2 H   2s  :    -0.15)                                                                                  
                                                                                                                          
               Molecular Orbital No.   4:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.06554 a.u.                                                                  
               (   1 Li  2p+1:    -1.01)                                                                                  
                                                                                                                          
               Molecular Orbital No.   5:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.06554 a.u.                                                                  
               (   1 Li  2p-1:     1.01)                                                                                  
                                                                                                                          
               Molecular Orbital No.   6:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.13091 a.u.                                                                  
               (   1 Li  2s  :     0.27) (   1 Li  3s  :    -0.83) (   1 Li  1p0 :     0.71)                              
               (   1 Li  2p0 :    -1.33) (   2 H   2s  :     0.31)                                                        
                                                                                                                          
               Molecular Orbital No.   7:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.20610 a.u.                                                                  
               (   1 Li  2s  :    -1.69) (   1 Li  3s  :     1.75) (   1 Li  1p0 :     1.09)                              
               (   1 Li  2p0 :    -0.67) (   2 H   2s  :    -0.32)                                                        
                                                                                                                          
SCF converged in 10 iterations.
Total Energy: -7.8713555781 au
                                                                                                                          
                                            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-10                                                              
                   ERI Screening Scheme            : Cauchy Schwarz + Density                                             
                   ERI Screening Mode              : Dynamic                                                              
                   ERI Screening Threshold         : 1.0e-12                                                              
                   Linear Dependence Threshold     : 1.0e-06                                                              
                                                                                                                          
* Info * Nuclear repulsion energy: 9.3436381577 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 * SAD initial guess computed in 0.00 sec.                                                                          
                                                                                                                          
* Info * Starting Reduced Basis SCF calculation...                                                                        
* Info * ...done. SCF energy in reduced basis set: -75.983338616414 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       -75.983338639967    0.0000000000      0.00025895      0.00005820      0.00000000                
                  2       -75.983338647052   -0.0000000071      0.00006338      0.00001527      0.00006095                
                  3       -75.983338648327   -0.0000000013      0.00000770      0.00000104      0.00005132                
                  4       -75.983338648340   -0.0000000000      0.00000056      0.00000014      0.00000520                
                  5       -75.983338648340   -0.0000000000      0.00000004      0.00000001      0.00000050                
                  6       -75.983338648340    0.0000000000      0.00000000      0.00000000      0.00000004                
                  7       -75.983338648340    0.0000000000      0.00000000      0.00000000      0.00000000                
                  8       -75.983338648340   -0.0000000000      0.00000000      0.00000000      0.00000000                
                                                                                                                          
               *** SCF converged in 8 iterations. Time: 0.04 sec.                                                         
                                                                                                                          
               Spin-Restricted Hartree-Fock:                                                                              
               -----------------------------                                                                              
               Total Energy                       :      -75.9833386483 a.u.                                              
               Electronic Energy                  :      -85.3269768060 a.u.                                              
               Nuclear Repulsion Energy           :        9.3436381577 a.u.                                              
               ------------------------------------                                                                       
               Gradient Norm                      :        0.0000000000 a.u.                                              
                                                                                                                          
                                                                                                                          
               Ground State Information                                                                                   
               ------------------------                                                                                   
               Charge of Molecule            :  0.0                                                                       
               Multiplicity (2S+1)           :  1.0                                                                       
               Magnetic Quantum Number (M_S) :  0.0                                                                       
                                                                                                                          
                                                                                                                          
                                                 Spin Restricted Orbitals                                                 
                                                 ------------------------                                                 
                                                                                                                          
               Molecular Orbital No.   1:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:  -20.55797 a.u.                                                                  
               (   1 O   1s  :    -1.00)                                                                                  
                                                                                                                          
               Molecular Orbital No.   2:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:   -1.36562 a.u.                                                                  
               (   1 O   1s  :     0.21) (   1 O   2s  :    -0.47) (   1 O   3s  :    -0.47)                              
                                                                                                                          
               Molecular Orbital No.   3:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:   -0.71725 a.u.                                                                  
               (   1 O   1p-1:     0.51) (   1 O   2p-1:     0.27) (   2 H   1s  :     0.27)                              
               (   3 H   1s  :    -0.27)                                                                                  
                                                                                                                          
               Molecular Orbital No.   4:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:   -0.56447 a.u.                                                                  
               (   1 O   2s  :     0.18) (   1 O   3s  :     0.31) (   1 O   1p0 :    -0.55)                              
               (   1 O   2p0 :    -0.40)                                                                                  
                                                                                                                          
               Molecular Orbital No.   5:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:   -0.50264 a.u.                                                                  
               (   1 O   1p+1:     0.64) (   1 O   2p+1:     0.51)                                                        
                                                                                                                          
               Molecular Orbital No.   6:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.20696 a.u.                                                                  
               (   1 O   3s  :    -1.21) (   1 O   1p0 :    -0.23) (   1 O   2p0 :    -0.47)                              
               (   2 H   2s  :     1.00) (   3 H   2s  :     1.00)                                                        
                                                                                                                          
               Molecular Orbital No.   7:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.30383 a.u.                                                                  
               (   1 O   1p-1:    -0.33) (   1 O   2p-1:    -0.83) (   2 H   2s  :     1.42)                              
               (   3 H   2s  :    -1.42)                                                                                  
                                                                                                                          
               Molecular Orbital No.   8:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    1.05847 a.u.                                                                  
               (   1 O   2p-1:    -0.83) (   2 H   1s  :     0.97) (   2 H   2s  :    -0.43)                              
               (   3 H   1s  :    -0.97) (   3 H   2s  :     0.43)                                                        
                                                                                                                          
               Molecular Orbital No.   9:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    1.16412 a.u.                                                                  
               (   1 O   1p+1:     0.96) (   1 O   2p+1:    -1.04)                                                        
                                                                                                                          
               Molecular Orbital No.  10:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    1.20359 a.u.                                                                  
               (   1 O   1p0 :    -0.98) (   1 O   2p0 :     0.92) (   2 H   1s  :     0.32)                              
               (   2 H   2s  :    -0.38) (   3 H   1s  :     0.32) (   3 H   2s  :    -0.38)                              
                                                                                                                          
SCF converged in 8 iterations.
Total Energy: -75.9833386483 au

Excitation energies of individual systems#

# Number of excited states for H2O, twice as many for LiH, thrice for composite system
nexc = 3

# Calculate ADC(1) and ADC(2) excitation energies with Gator
adc1_drv = AdcOneDriver(scf_lih.comm)
adc1_drv.conv_thresh = 1e-5
adc1_results_h2o = adc1_drv.compute(h2o, basis_h2o, scf_h2o.scf_tensors)
adc1_drv.nstates = 2 * nexc
adc1_results_lih = adc1_drv.compute(lih, basis_lih, scf_lih.scf_tensors)

adc2_drv = AdcTwoDriver(scf_lih.comm, scf_lih.ostream)
adc2_drv.conv_thresh = 1e-5
adc2_results_h2o = adc2_drv.compute(h2o, basis_h2o, scf_h2o.scf_tensors)
adc2_drv.nstates = 2 * nexc
adc2_results_lih = adc2_drv.compute(lih, basis_lih, scf_lih.scf_tensors)
                                                                                                                          
                                                   ADC(1) Driver Setup                                                    
                                                  =====================                                                   
                                                                                                                          
                               Number Of Excited States  : 3                                                              
                               Max. Number Of Iterations : 50                                                             
                               Convergence Threshold     : 1.0e-05                                                        
                               ERI screening scheme      : Cauchy Schwarz + Density                                       
                               ERI Screening Threshold   : 1.0e-15                                                        
                                                                                                                          
* Info * Number of occupied orbitals: 5                                                                                   
* Info * Number of virtual orbitals: 8                                                                                    
                                                                                                                          
                                                MO Integrals Driver Setup                                                 
                                               ===========================                                                
                                                                                                                          
                               ERI screening scheme        : QQ_DEN                                                       
                               ERI Screening Threshold     : 1.0e-15                                                      
                               Batch Size of Fock Matrices : 100                                                          
                                                                                                                          
* Info * Processing Fock builds for the OO block...                                                                       
* Info *   batch 1/1                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the OO block done in 0.01 sec. Load imb.: 0.0 %                                     
                                                                                                                          
                   *** Iteration:   1 * Reduced Space:    3 * Residues (Max,Min): 1.49e-01 and 1.15e-01                   
                                                                                                                          
                   State  1: 0.36669368 au Residual Norm: 0.11478601                                                      
                   State  2: 0.44840498 au Residual Norm: 0.13775605                                                      
                   State  3: 0.46360595 au Residual Norm: 0.14905941                                                      
                                                                                                                          
                   *** Iteration:   2 * Reduced Space:    6 * Residues (Max,Min): 4.52e-02 and 1.08e-03                   
                                                                                                                          
                   State  1: 0.35291511 au Residual Norm: 0.01136344                                                      
                   State  2: 0.42544431 au Residual Norm: 0.00107731                                                      
                   State  3: 0.44554970 au Residual Norm: 0.04521708                                                      
                                                                                                                          
                   *** Iteration:   3 * Reduced Space:    9 * Residues (Max,Min): 7.67e-03 and 1.36e-05                   
                                                                                                                          
                   State  1: 0.35279762 au Residual Norm: 0.00151889                                                      
                   State  2: 0.42544293 au Residual Norm: 0.00001364                                                      
                   State  3: 0.44369379 au Residual Norm: 0.00767191                                                      
                                                                                                                          
                   *** Iteration:   4 * Reduced Space:   12 * Residues (Max,Min): 2.29e-03 and 9.39e-16                   
                                                                                                                          
                   State  1: 0.35279702 au Residual Norm: 0.00033748                                                      
                   State  2: 0.42544293 au Residual Norm: 0.00000000                                                      
                   State  3: 0.44360836 au Residual Norm: 0.00228847                                                      
                                                                                                                          
                   *** Iteration:   5 * Reduced Space:   15 * Residues (Max,Min): 9.00e-04 and 4.47e-15                   
                                                                                                                          
                   State  1: 0.35279701 au Residual Norm: 0.00000656                                                      
                   State  2: 0.42544293 au Residual Norm: 0.00000000                                                      
                   State  3: 0.44360608 au Residual Norm: 0.00089980                                                      
                                                                                                                          
                   *** Iteration:   6 * Reduced Space:   18 * Residues (Max,Min): 1.10e-04 and 8.68e-15                   
                                                                                                                          
                   State  1: 0.35279701 au Residual Norm: 0.00000010                                                      
                   State  2: 0.42544293 au Residual Norm: 0.00000000                                                      
                   State  3: 0.44360599 au Residual Norm: 0.00010969                                                      
                                                                                                                          
                   *** Iteration:   7 * Reduced Space:   21 * Residues (Max,Min): 4.35e-06 and 8.62e-15                   
                                                                                                                          
                   State  1: 0.35279701 au Residual Norm: 0.00000000                                                      
                   State  2: 0.42544293 au Residual Norm: 0.00000000                                                      
                   State  3: 0.44360599 au Residual Norm: 0.00000435                                                      
                                                                                                                          
               *** 3 excited states converged in 7 iterations. Time: 0.03 sec.                                            
                                                                                                                          
                                                                                                                          
               ADC(1) excited states                                                                                      
               ---------------------                                                                                      
               Excited State    S1:      0.35279701 a.u.      9.60010 eV    Osc.Str. 0.01663                              
               Excited State    S2:      0.42544293 a.u.     11.57689 eV    Osc.Str. 0.00000                              
               Excited State    S3:      0.44360599 a.u.     12.07113 eV    Osc.Str. 0.11240                              
                                                                                                                          
                                                                                                                          
                                                                                                                          
                                                   ADC(1) Driver Setup                                                    
                                                  =====================                                                   
                                                                                                                          
                               Number Of Excited States  : 6                                                              
                               Max. Number Of Iterations : 50                                                             
                               Convergence Threshold     : 1.0e-05                                                        
                               ERI screening scheme      : Cauchy Schwarz + Density                                       
                               ERI Screening Threshold   : 1.0e-15                                                        
                                                                                                                          
* Info * Number of occupied orbitals: 2                                                                                   
* Info * Number of virtual orbitals: 9                                                                                    
                                                                                                                          
                                                MO Integrals Driver Setup                                                 
                                               ===========================                                                
                                                                                                                          
                               ERI screening scheme        : QQ_DEN                                                       
                               ERI Screening Threshold     : 1.0e-15                                                      
                               Batch Size of Fock Matrices : 100                                                          
                                                                                                                          
* Info * Processing Fock builds for the OO block...                                                                       
* Info *   batch 1/1                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the OO block done in 0.00 sec. Load imb.: 0.0 %                                     
                                                                                                                          
                   *** Iteration:   1 * Reduced Space:    6 * Residues (Max,Min): 3.60e-02 and 1.47e-02                   
                                                                                                                          
                   State  1: 0.15889294 au Residual Norm: 0.01952183                                                      
                   State  2: 0.20608326 au Residual Norm: 0.01468372                                                      
                   State  3: 0.21659287 au Residual Norm: 0.03602681                                                      
                   State  4: 0.28273182 au Residual Norm: 0.01773908                                                      
                   State  5: 0.33269430 au Residual Norm: 0.01820955                                                      
                   State  6: 0.33539479 au Residual Norm: 0.02560730                                                      
                                                                                                                          
                   *** Iteration:   2 * Reduced Space:   12 * Residues (Max,Min): 4.60e-02 and 2.22e-14                   
                                                                                                                          
                   State  1: 0.15870155 au Residual Norm: 0.00048111                                                      
                   State  2: 0.20597994 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20599656 au Residual Norm: 0.00587184                                                      
                   State  4: 0.28256950 au Residual Norm: 0.00084005                                                      
                   State  5: 0.33252586 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33360191 au Residual Norm: 0.04595089                                                      
                                                                                                                          
                   *** Iteration:   3 * Reduced Space:   16 * Residues (Max,Min): 4.26e-05 and 7.92e-12                   
                                                                                                                          
                   State  1: 0.15870142 au Residual Norm: 0.00001296                                                      
                   State  2: 0.20597994 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20597994 au Residual Norm: 0.00000075                                                      
                   State  4: 0.28256876 au Residual Norm: 0.00004260                                                      
                   State  5: 0.33252586 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33252586 au Residual Norm: 0.00001253                                                      
                                                                                                                          
                   *** Iteration:   4 * Reduced Space:   18 * Residues (Max,Min): 1.24e-15 and 1.33e-16                   
                                                                                                                          
                   State  1: 0.15870142 au Residual Norm: 0.00000000                                                      
                   State  2: 0.20597994 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20597994 au Residual Norm: 0.00000000                                                      
                   State  4: 0.28256876 au Residual Norm: 0.00000000                                                      
                   State  5: 0.33252586 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33252586 au Residual Norm: 0.00000000                                                      
                                                                                                                          
               *** 6 excited states converged in 4 iterations. Time: 0.01 sec.                                            
                                                                                                                          
                                                                                                                          
               ADC(1) excited states                                                                                      
               ---------------------                                                                                      
               Excited State    S1:      0.15870142 a.u.      4.31849 eV    Osc.Str. 0.02843                              
               Excited State    S2:      0.20597994 a.u.      5.60500 eV    Osc.Str. 0.25169                              
               Excited State    S3:      0.20597994 a.u.      5.60500 eV    Osc.Str. 0.25169                              
               Excited State    S4:      0.28256876 a.u.      7.68909 eV    Osc.Str. 0.04341                              
               Excited State    S5:      0.33252586 a.u.      9.04849 eV    Osc.Str. 0.02652                              
               Excited State    S6:      0.33252586 a.u.      9.04849 eV    Osc.Str. 0.02652                              
                                                                                                                          
                                                                                                                          
                                                MO Integrals Driver Setup                                                 
                                               ===========================                                                
                                                                                                                          
                               ERI screening scheme        : QQ_DEN                                                       
                               ERI Screening Threshold     : 1.0e-15                                                      
                               Batch Size of Fock Matrices : 100                                                          
                                                                                                                          
* Info * Processing Fock builds for the OO block...                                                                       
* Info *   batch 1/1                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the OO block done in 0.01 sec. Load imb.: 0.0 %                                     
                                                                                                                          
* Info * Processing Fock builds for the VV block...                                                                       
* Info *   batch 1/1                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the VV block done in 0.01 sec. Load imb.: 0.0 %                                     
                                                                                                                          
* Info * Processing Fock builds for the OV block...                                                                       
* Info *   batch 1/1                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the OV block done in 0.01 sec. Load imb.: 0.0 %                                     
                                                                                                                          
* Info * 146.78 kB of memory used for MO integrals on the master node.                                                    
* Info * 4.06 GB of memory available on the master node.                                                                  
* Info * Time spent in calculating memory usage: 0.00 sec.                                                                
                                                                                                                          
                                                                                                                          
                                                   ADC(2) Driver Setup                                                    
                                                  =====================                                                   
                                                                                                                          
                               Number Of Excited States  : 3                                                              
                               Max. Number Of Iterations : 150                                                            
                               Convergence Threshold     : 1.0e-05                                                        
                               ERI screening scheme      : Cauchy Schwarz + Density                                       
                               ERI Screening Threshold   : 1.0e-15                                                        
                                                                                                                          
                                                                                                                          
                                                   ADC(1) Driver Setup                                                    
                                                  =====================                                                   
                                                                                                                          
                               Number Of Excited States  : 3                                                              
                               Max. Number Of Iterations : 50                                                             
                               Convergence Threshold     : 1.0e-04                                                        
                               ERI screening scheme      : Cauchy Schwarz + Density                                       
                               ERI Screening Threshold   : 1.0e-15                                                        
                                                                                                                          
* Info * Number of occupied orbitals: 5                                                                                   
* Info * Number of virtual orbitals: 8                                                                                    
                                                                                                                          
                   *** Iteration:   1 * Reduced Space:    3 * Residues (Max,Min): 1.49e-01 and 1.15e-01                   
                                                                                                                          
                   State  1: 0.36669368 au Residual Norm: 0.11478601                                                      
                   State  2: 0.44840498 au Residual Norm: 0.13775605                                                      
                   State  3: 0.46360595 au Residual Norm: 0.14905941                                                      
                                                                                                                          
                   *** Iteration:   2 * Reduced Space:    6 * Residues (Max,Min): 4.52e-02 and 1.08e-03                   
                                                                                                                          
                   State  1: 0.35291511 au Residual Norm: 0.01136344                                                      
                   State  2: 0.42544431 au Residual Norm: 0.00107731                                                      
                   State  3: 0.44554970 au Residual Norm: 0.04521708                                                      
                                                                                                                          
                   *** Iteration:   3 * Reduced Space:    9 * Residues (Max,Min): 7.67e-03 and 1.36e-05                   
                                                                                                                          
                   State  1: 0.35279762 au Residual Norm: 0.00151889                                                      
                   State  2: 0.42544293 au Residual Norm: 0.00001364                                                      
                   State  3: 0.44369379 au Residual Norm: 0.00767191                                                      
                                                                                                                          
                   *** Iteration:   4 * Reduced Space:   12 * Residues (Max,Min): 2.29e-03 and 6.27e-16                   
                                                                                                                          
                   State  1: 0.35279702 au Residual Norm: 0.00033748                                                      
                   State  2: 0.42544293 au Residual Norm: 0.00000000                                                      
                   State  3: 0.44360836 au Residual Norm: 0.00228847                                                      
                                                                                                                          
                   *** Iteration:   5 * Reduced Space:   15 * Residues (Max,Min): 8.30e-04 and 2.34e-15                   
                                                                                                                          
                   State  1: 0.35279701 au Residual Norm: 0.00000656                                                      
                   State  2: 0.42544293 au Residual Norm: 0.00000000                                                      
                   State  3: 0.44360607 au Residual Norm: 0.00082971                                                      
                                                                                                                          
                   *** Iteration:   6 * Reduced Space:   18 * Residues (Max,Min): 1.04e-04 and 6.80e-15                   
                                                                                                                          
                   State  1: 0.35279701 au Residual Norm: 0.00000010                                                      
                   State  2: 0.42544293 au Residual Norm: 0.00000000                                                      
                   State  3: 0.44360599 au Residual Norm: 0.00010364                                                      
                                                                                                                          
                   *** Iteration:   7 * Reduced Space:   20 * Residues (Max,Min): 4.05e-06 and 5.71e-15                   
                                                                                                                          
                   State  1: 0.35279701 au Residual Norm: 0.00000000                                                      
                   State  2: 0.42544293 au Residual Norm: 0.00000000                                                      
                   State  3: 0.44360599 au Residual Norm: 0.00000405                                                      
                                                                                                                          
               *** 3 excited states converged in 7 iterations. Time: 0.01 sec.                                            
                                                                                                                          
                                                                                                                          
               ADC(1) excited states                                                                                      
               ---------------------                                                                                      
               Excited State    S1:      0.35279701 a.u.      9.60010 eV    Osc.Str. 0.01663                              
               Excited State    S2:      0.42544293 a.u.     11.57689 eV    Osc.Str. 0.00000                              
               Excited State    S3:      0.44360599 a.u.     12.07113 eV    Osc.Str. 0.11240                              
                                                                                                                          
                                                                                                                          
* Info * Number of occupied orbitals: 5                                                                                   
* Info * Number of virtual orbitals: 8                                                                                    
                                                                                                                          
* Info * MP2 correlation energy: -0.1274706653 a.u. Time: 0.00 sec.                                                       
                                                                                                                          
* Info * Initial guess generated in 0.02 sec.                                                                             
                                                                                                                          
                   *** Iteration:   1 * Reduced Space:    3 * Residual Norm: 6.71e-02                                     
                                                                                                                          
                   State  1:     0.35279701 au  Eigenvalue Error: 3.98e-02                                                
                                                                                                                          
                   *** Iteration:   2 * Reduced Space:    6 * Residual Norm: 5.50e-03                                     
                                                                                                                          
                   State  1:     0.31462536 au  Eigenvalue Error: 2.87e-03                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:   3 * Reduced Space:    6 * Residual Norm: 3.42e-04                                     
                                                                                                                          
                   State  1:     0.31187300 au  Eigenvalue Error: 1.74e-05                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:   4 * Reduced Space:    6 * Residual Norm: 6.09e-05                                     
                                                                                                                          
                   State  1:     0.31185626 au  Eigenvalue Error: 1.01e-07                                                
                                                                                                                          
                   *** Iteration:   5 * Reduced Space:    9 * Residual Norm: 1.21e-05                                     
                                                                                                                          
                   State  1:     0.31185616 au  Eigenvalue Error: 9.50e-10                                                
                                                                                                                          
                   *** Iteration:   6 * Reduced Space:   12 * Residual Norm: 3.49e-07                                     
                                                                                                                          
                   State  1:     0.31185616 au  converged                                                                 
                                                                                                                          
* Info * State 1 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:   7 * Reduced Space:   15 * Residual Norm: 3.61e-03                                     
                                                                                                                          
                   State  1:     0.31185616 au  converged                                                                 
                   State  2:     0.40125793 au  Eigenvalue Error: 3.21e-03                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:   8 * Reduced Space:    6 * Residual Norm: 5.23e-05                                     
                                                                                                                          
                   State  1:     0.31185616 au  converged                                                                 
                   State  2:     0.39816818 au  Eigenvalue Error: 3.08e-06                                                
                                                                                                                          
                   *** Iteration:   9 * Reduced Space:    9 * Residual Norm: 9.73e-07                                     
                                                                                                                          
                   State  1:     0.31185616 au  converged                                                                 
                   State  2:     0.39816521 au  Eigenvalue Error: 3.57e-09                                                
                                                                                                                          
                   *** Iteration:  10 * Reduced Space:   12 * Residual Norm: 5.74e-08                                     
                                                                                                                          
                   State  1:     0.31185616 au  converged                                                                 
                   State  2:     0.39816520 au  converged                                                                 
                                                                                                                          
* Info * State 2 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:  11 * Reduced Space:   15 * Residual Norm: 2.40e-04                                     
                                                                                                                          
                   State  1:     0.31185616 au  converged                                                                 
                   State  2:     0.39816520 au  converged                                                                 
                   State  3:     0.40309707 au  Eigenvalue Error: 2.19e-04                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:  12 * Reduced Space:    6 * Residual Norm: 1.01e-05                                     
                                                                                                                          
                   State  1:     0.31185616 au  converged                                                                 
                   State  2:     0.39816520 au  converged                                                                 
                   State  3:     0.40288697 au  Eigenvalue Error: 1.01e-08                                                
                                                                                                                          
                   *** Iteration:  13 * Reduced Space:    9 * Residual Norm: 2.89e-06                                     
                                                                                                                          
                   State  1:     0.31185616 au  converged                                                                 
                   State  2:     0.39816520 au  converged                                                                 
                   State  3:     0.40288696 au  converged                                                                 
                                                                                                                          
* Info * State 3 is converged.                                                                                            
                                                                                                                          
* Info * Total number of sigma builds: 64                                                                                 
* Info * Average time per sigma build: 0.00 sec                                                                           
                                                                                                                          
* Info * Total number of d_sigma builds: 13                                                                               
* Info * Average time per d_sigma build: 0.00 sec                                                                         
                                                                                                                          
               *** 3 excited states converged in 13 iterations.  Time: 0.39 sec.                                          
                                                                                                                          
                                                                                                                          
               ADC(2) excited states                                                                                      
               ---------------------                                                                                      
               Excited State    S1:      0.31185616 a.u.      8.48604 eV    |v1|^2=0.9595                                 
               Excited State    S2:      0.39816520 a.u.     10.83463 eV    |v1|^2=0.9638                                 
               Excited State    S3:      0.40288696 a.u.     10.96311 eV    |v1|^2=0.9573                                 
                                                                                                                          
               State S1                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO         0.9772                                                                            
               HOMO     -> LUMO+4       0.0541                                                                            
                                                                                                                          
               State S2                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+1       0.9765                                                                            
               HOMO     -> LUMO+6      -0.0824                                                                            
               HOMO     -> LUMO+2       0.0588                                                                            
                                                                                                                          
               State S3                                                                                                   
               -------------------------------                                                                            
               HOMO-1   -> LUMO         0.9722                                                                            
               HOMO-2   -> LUMO+1       0.0749                                                                            
                                                                                                                          
                                                                                                                          
                                                MO Integrals Driver Setup                                                 
                                               ===========================                                                
                                                                                                                          
                               ERI screening scheme        : QQ_DEN                                                       
                               ERI Screening Threshold     : 1.0e-15                                                      
                               Batch Size of Fock Matrices : 100                                                          
                                                                                                                          
* Info * Processing Fock builds for the OO block...                                                                       
* Info *   batch 1/1                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the OO block done in 0.01 sec. Load imb.: 0.0 %                                     
                                                                                                                          
* Info * Processing Fock builds for the VV block...                                                                       
* Info *   batch 1/1                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the VV block done in 0.01 sec. Load imb.: 0.0 %                                     
                                                                                                                          
* Info * Processing Fock builds for the OV block...                                                                       
* Info *   batch 1/1                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the OV block done in 0.01 sec. Load imb.: 0.0 %                                     
                                                                                                                          
* Info * 71.95 kB of memory used for MO integrals on the master node.                                                     
* Info * 4.06 GB of memory available on the master node.                                                                  
* Info * Time spent in calculating memory usage: 0.00 sec.                                                                
                                                                                                                          
                                                                                                                          
                                                   ADC(2) Driver Setup                                                    
                                                  =====================                                                   
                                                                                                                          
                               Number Of Excited States  : 6                                                              
                               Max. Number Of Iterations : 150                                                            
                               Convergence Threshold     : 1.0e-05                                                        
                               ERI screening scheme      : Cauchy Schwarz + Density                                       
                               ERI Screening Threshold   : 1.0e-15                                                        
                                                                                                                          
                                                                                                                          
                                                   ADC(1) Driver Setup                                                    
                                                  =====================                                                   
                                                                                                                          
                               Number Of Excited States  : 6                                                              
                               Max. Number Of Iterations : 50                                                             
                               Convergence Threshold     : 1.0e-04                                                        
                               ERI screening scheme      : Cauchy Schwarz + Density                                       
                               ERI Screening Threshold   : 1.0e-15                                                        
                                                                                                                          
* Info * Number of occupied orbitals: 2                                                                                   
* Info * Number of virtual orbitals: 9                                                                                    
                                                                                                                          
                   *** Iteration:   1 * Reduced Space:    6 * Residues (Max,Min): 3.60e-02 and 1.47e-02                   
                                                                                                                          
                   State  1: 0.15889294 au Residual Norm: 0.01952183                                                      
                   State  2: 0.20608326 au Residual Norm: 0.01468372                                                      
                   State  3: 0.21659287 au Residual Norm: 0.03602681                                                      
                   State  4: 0.28273182 au Residual Norm: 0.01773908                                                      
                   State  5: 0.33269430 au Residual Norm: 0.01820955                                                      
                   State  6: 0.33539479 au Residual Norm: 0.02560730                                                      
                                                                                                                          
                   *** Iteration:   2 * Reduced Space:   12 * Residues (Max,Min): 4.60e-02 and 7.39e-15                   
                                                                                                                          
                   State  1: 0.15870155 au Residual Norm: 0.00048111                                                      
                   State  2: 0.20597994 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20599656 au Residual Norm: 0.00587184                                                      
                   State  4: 0.28256950 au Residual Norm: 0.00084005                                                      
                   State  5: 0.33252586 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33360191 au Residual Norm: 0.04595089                                                      
                                                                                                                          
                   *** Iteration:   3 * Reduced Space:   16 * Residues (Max,Min): 4.62e-05 and 4.74e-12                   
                                                                                                                          
                   State  1: 0.15870142 au Residual Norm: 0.00001296                                                      
                   State  2: 0.20597994 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20597994 au Residual Norm: 0.00000322                                                      
                   State  4: 0.28256876 au Residual Norm: 0.00004261                                                      
                   State  5: 0.33252586 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33252586 au Residual Norm: 0.00004621                                                      
                                                                                                                          
               *** 6 excited states converged in 3 iterations. Time: 0.01 sec.                                            
                                                                                                                          
                                                                                                                          
               ADC(1) excited states                                                                                      
               ---------------------                                                                                      
               Excited State    S1:      0.15870142 a.u.      4.31849 eV    Osc.Str. 0.02843                              
               Excited State    S2:      0.20597994 a.u.      5.60500 eV    Osc.Str. 0.25169                              
               Excited State    S3:      0.20597994 a.u.      5.60500 eV    Osc.Str. 0.25169                              
               Excited State    S4:      0.28256876 a.u.      7.68909 eV    Osc.Str. 0.04341                              
               Excited State    S5:      0.33252586 a.u.      9.04849 eV    Osc.Str. 0.02652                              
               Excited State    S6:      0.33252586 a.u.      9.04849 eV    Osc.Str. 0.02652                              
                                                                                                                          
                                                                                                                          
* Info * Number of occupied orbitals: 2                                                                                   
* Info * Number of virtual orbitals: 9                                                                                    
                                                                                                                          
* Info * MP2 correlation energy: -0.0136828971 a.u. Time: 0.00 sec.                                                       
                                                                                                                          
* Info * Initial guess generated in 0.01 sec.                                                                             
                                                                                                                          
                   *** Iteration:   1 * Reduced Space:    6 * Residual Norm: 1.82e-02                                     
                                                                                                                          
                   State  1:     0.15870142 au  Eigenvalue Error: 1.76e-02                                                
                                                                                                                          
                   *** Iteration:   2 * Reduced Space:   12 * Residual Norm: 1.26e-03                                     
                                                                                                                          
                   State  1:     0.14174508 au  Eigenvalue Error: 6.10e-05                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:   3 * Reduced Space:   12 * Residual Norm: 6.40e-05                                     
                                                                                                                          
                   State  1:     0.14168623 au  Eigenvalue Error: 1.61e-06                                                
                                                                                                                          
                   *** Iteration:   4 * Reduced Space:   16 * Residual Norm: 1.44e-06                                     
                                                                                                                          
                   State  1:     0.14168468 au  Eigenvalue Error: 2.14e-09                                                
                                                                                                                          
                   *** Iteration:   5 * Reduced Space:   18 * Residual Norm: 1.22e-08                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                                                                                                                          
* Info * State 1 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:   6 * Reduced Space:   18 * Residual Norm: 1.68e-03                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18797246 au  Eigenvalue Error: 1.66e-03                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:   7 * Reduced Space:   12 * Residual Norm: 3.55e-07                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637063 au  Eigenvalue Error: 3.55e-07                                                
                                                                                                                          
                   *** Iteration:   8 * Reduced Space:   17 * Residual Norm: 1.65e-09                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637029 au  converged                                                                 
                                                                                                                          
* Info * State 2 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:   9 * Reduced Space:   18 * Residual Norm: 1.36e-08                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637029 au  converged                                                                 
                   State  3:     0.18637028 au  Eigenvalue Error: 1.35e-08                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:  10 * Reduced Space:   12 * Residual Norm: 3.00e-10                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637029 au  converged                                                                 
                   State  3:     0.18637029 au  converged                                                                 
                                                                                                                          
* Info * State 3 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:  11 * Reduced Space:   17 * Residual Norm: 2.81e-03                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637029 au  converged                                                                 
                   State  3:     0.18637029 au  converged                                                                 
                   State  4:     0.26955011 au  Eigenvalue Error: 2.69e-03                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:  12 * Reduced Space:   12 * Residual Norm: 3.82e-05                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637029 au  converged                                                                 
                   State  3:     0.18637029 au  converged                                                                 
                   State  4:     0.26694944 au  Eigenvalue Error: 1.60e-06                                                
                                                                                                                          
                   *** Iteration:  13 * Reduced Space:   17 * Residual Norm: 1.94e-07                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637029 au  converged                                                                 
                   State  3:     0.18637029 au  converged                                                                 
                   State  4:     0.26694790 au  converged                                                                 
                                                                                                                          
* Info * State 4 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:  14 * Reduced Space:   18 * Residual Norm: 1.85e-03                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637029 au  converged                                                                 
                   State  3:     0.18637029 au  converged                                                                 
                   State  4:     0.26694790 au  converged                                                                 
                   State  5:     0.32063129 au  Eigenvalue Error: 1.75e-03                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:  15 * Reduced Space:   12 * Residual Norm: 3.53e-05                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637029 au  converged                                                                 
                   State  3:     0.18637029 au  converged                                                                 
                   State  4:     0.26694790 au  converged                                                                 
                   State  5:     0.31893708 au  Eigenvalue Error: 1.71e-06                                                
                                                                                                                          
                   *** Iteration:  16 * Reduced Space:   17 * Residual Norm: 2.07e-08                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637029 au  converged                                                                 
                   State  3:     0.18637029 au  converged                                                                 
                   State  4:     0.26694790 au  converged                                                                 
                   State  5:     0.31893873 au  converged                                                                 
                                                                                                                          
* Info * State 5 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:  17 * Reduced Space:   18 * Residual Norm: 3.06e-05                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637029 au  converged                                                                 
                   State  3:     0.18637029 au  converged                                                                 
                   State  4:     0.26694790 au  converged                                                                 
                   State  5:     0.31893873 au  converged                                                                 
                   State  6:     0.31996133 au  Eigenvalue Error: 2.96e-05                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:  18 * Reduced Space:   12 * Residual Norm: 3.43e-10                                     
                                                                                                                          
                   State  1:     0.14168467 au  converged                                                                 
                   State  2:     0.18637029 au  converged                                                                 
                   State  3:     0.18637029 au  converged                                                                 
                   State  4:     0.26694790 au  converged                                                                 
                   State  5:     0.31893873 au  converged                                                                 
                   State  6:     0.31993255 au  converged                                                                 
                                                                                                                          
* Info * State 6 is converged.                                                                                            
                                                                                                                          
* Info * Total number of sigma builds: 131                                                                                
* Info * Average time per sigma build: 0.00 sec                                                                           
                                                                                                                          
* Info * Total number of d_sigma builds: 18                                                                               
* Info * Average time per d_sigma build: 0.00 sec                                                                         
                                                                                                                          
               *** 6 excited states converged in 18 iterations.  Time: 0.44 sec.                                          
                                                                                                                          
                                                                                                                          
               ADC(2) excited states                                                                                      
               ---------------------                                                                                      
               Excited State    S1:      0.14168467 a.u.      3.85544 eV    |v1|^2=0.9645                                 
               Excited State    S2:      0.18637029 a.u.      5.07139 eV    |v1|^2=0.9634                                 
               Excited State    S3:      0.18637029 a.u.      5.07139 eV    |v1|^2=0.9634                                 
               Excited State    S4:      0.26694790 a.u.      7.26402 eV    |v1|^2=0.9652                                 
               Excited State    S5:      0.31893873 a.u.      8.67876 eV    |v1|^2=0.9662                                 
               Excited State    S6:      0.31993255 a.u.      8.70581 eV    |v1|^2=0.9719                                 
                                                                                                                          
               State S1                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO         0.9487                                                                            
               HOMO     -> LUMO+4      -0.2366                                                                            
               HOMO     -> LUMO+3      -0.0915                                                                            
                                                                                                                          
               State S2                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+2       0.9006                                                                            
               HOMO     -> LUMO+1      -0.2807                                                                            
               HOMO     -> LUMO+6       0.2637                                                                            
               HOMO     -> LUMO+5      -0.0627                                                                            
                                                                                                                          
               State S3                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+1      -0.9006                                                                            
               HOMO     -> LUMO+2      -0.2807                                                                            
               HOMO     -> LUMO+5      -0.2637                                                                            
               HOMO     -> LUMO+6      -0.0627                                                                            
                                                                                                                          
               State S4                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+3       0.9756                                                                            
               HOMO     -> LUMO+4      -0.0876                                                                            
               HOMO     -> LUMO         0.0741                                                                            
                                                                                                                          
               State S5                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+4      -0.9512                                                                            
               HOMO     -> LUMO        -0.2365                                                                            
               HOMO     -> LUMO+3      -0.0707                                                                            
                                                                                                                          
               State S6                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+6       0.9232                                                                            
               HOMO     -> LUMO+2      -0.2551                                                                            
               HOMO     -> LUMO+5      -0.2191                                                                            
               HOMO     -> LUMO+1       0.0794                                                                            
                                                                                                                          
                                                                                                                          
print("LiH ADC(1) excitation energies:\n", adc1_results_lih['eigenvalues'])
print("\nH2O ADC(1) excitation energies:\n", adc1_results_h2o['eigenvalues'])
print("\nLiH ADC(2) excitation energies:\n", adc2_results_lih['eigenvalues'])
print("\nH2O ADC(2) excitation energies:\n", adc2_results_h2o['eigenvalues'])
LiH ADC(1) excitation energies:
 [0.1587  0.20598 0.20598 0.28257 0.33253 0.33253]

H2O ADC(1) excitation energies:
 [0.3528  0.42544 0.44361]

LiH ADC(2) excitation energies:
 [0.14168 0.18637 0.18637 0.26695 0.31894 0.31993]

H2O ADC(2) excitation energies:
 [0.31186 0.39817 0.40289]

Composite system#

# LiH and H2O molecules 100 Å apart
lih_h2o_xyz="""5

Li  0.000000             0.000000          0.000000
H   0.000000             0.000000          1.000000
O   100.000000000000     0.000000000000    0.000000000000                         
H   100.000000000000     0.740848095288    0.582094932012                         
H   100.000000000000    -0.740848095288    0.582094932012
"""
lih_h2o = vlx.Molecule.from_xyz_string(lih_h2o_xyz)
basis_lih_h2o = vlx.MolecularBasis.read(lih_h2o, basis_set_label)
# Run SCF of composite system and compare to sum of individual SCF energies
scf_lih_h2o = gator.run_scf(lih_h2o, basis_lih_h2o, conv_thresh=1e-10, verbose=False)
print("Sum of SCF:  ", e_scf_lih + e_scf_h2o, "au")
                                                                                                                          
                                            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-10                                                              
                   ERI Screening Scheme            : Cauchy Schwarz + Density                                             
                   ERI Screening Mode              : Dynamic                                                              
                   ERI Screening Threshold         : 1.0e-12                                                              
                   Linear Dependence Threshold     : 1.0e-06                                                              
                                                                                                                          
* Info * Nuclear repulsion energy: 11.1428367662 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 * SAD initial guess computed in 0.00 sec.                                                                          
                                                                                                                          
* Info * Starting Reduced Basis SCF calculation...                                                                        
* Info * ...done. SCF energy in reduced basis set: -83.830261078187 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       -83.851291164519    0.0000000000      0.06853064      0.01397404      0.00000000                
                  2       -83.854364302432   -0.0030731379      0.01618988      0.00435148      0.13926994                
                  3       -83.854693323837   -0.0003290214      0.00163475      0.00029701      0.03162949                
                  4       -83.854694487659   -0.0000011638      0.00025024      0.00004407      0.00227174                
                  5       -83.854694519691   -0.0000000320      0.00005777      0.00001199      0.00043591                
                  6       -83.854694522926   -0.0000000032      0.00001852      0.00000337      0.00013749                
                  7       -83.854694522959   -0.0000000000      0.00000075      0.00000011      0.00000514                
                  8       -83.854694522959   -0.0000000000      0.00000010      0.00000002      0.00000053                
                  9       -83.854694522959   -0.0000000000      0.00000003      0.00000001      0.00000013                
                 10       -83.854694522959   -0.0000000000      0.00000000      0.00000000      0.00000004                
                 11       -83.854694522959    0.0000000000      0.00000000      0.00000000      0.00000000                
                 12       -83.854694522959   -0.0000000000      0.00000000      0.00000000      0.00000000                
                                                                                                                          
               *** SCF converged in 12 iterations. Time: 0.08 sec.                                                        
                                                                                                                          
               Spin-Restricted Hartree-Fock:                                                                              
               -----------------------------                                                                              
               Total Energy                       :      -83.8546945230 a.u.                                              
               Electronic Energy                  :      -94.9975312892 a.u.                                              
               Nuclear Repulsion Energy           :       11.1428367662 a.u.                                              
               ------------------------------------                                                                       
               Gradient Norm                      :        0.0000000000 a.u.                                              
                                                                                                                          
                                                                                                                          
               Ground State Information                                                                                   
               ------------------------                                                                                   
               Charge of Molecule            :  0.0                                                                       
               Multiplicity (2S+1)           :  1.0                                                                       
               Magnetic Quantum Number (M_S) :  0.0                                                                       
                                                                                                                          
                                                                                                                          
                                                 Spin Restricted Orbitals                                                 
                                                 ------------------------                                                 
                                                                                                                          
               Molecular Orbital No.   3:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:   -1.36562 a.u.                                                                  
               (   3 O   1s  :     0.21) (   3 O   2s  :    -0.47) (   3 O   3s  :    -0.47)                              
                                                                                                                          
               Molecular Orbital No.   4:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:   -0.71725 a.u.                                                                  
               (   3 O   1p-1:    -0.51) (   3 O   2p-1:    -0.27) (   4 H   1s  :    -0.27)                              
               (   5 H   1s  :     0.27)                                                                                  
                                                                                                                          
               Molecular Orbital No.   5:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:   -0.56447 a.u.                                                                  
               (   3 O   2s  :    -0.18) (   3 O   3s  :    -0.31) (   3 O   1p0 :     0.55)                              
               (   3 O   2p0 :     0.40)                                                                                  
                                                                                                                          
               Molecular Orbital No.   6:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:   -0.50264 a.u.                                                                  
               (   3 O   1p+1:    -0.64) (   3 O   2p+1:    -0.51)                                                        
                                                                                                                          
               Molecular Orbital No.   7:                                                                                 
               --------------------------                                                                                 
               Occupation: 2.000 Energy:   -0.32129 a.u.                                                                  
               (   1 Li  1s  :    -0.23) (   1 Li  2s  :     0.29) (   1 Li  1p0 :     0.38)                              
               (   2 H   1s  :     0.36) (   2 H   2s  :     0.31)                                                        
                                                                                                                          
               Molecular Orbital No.   8:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.01856 a.u.                                                                  
               (   1 Li  2s  :     0.16) (   1 Li  3s  :    -0.96) (   1 Li  2p0 :     0.50)                              
               (   2 H   2s  :     0.15)                                                                                  
                                                                                                                          
               Molecular Orbital No.   9:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.06554 a.u.                                                                  
               (   1 Li  2p-1:    -1.01)                                                                                  
                                                                                                                          
               Molecular Orbital No.  10:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.06554 a.u.                                                                  
               (   1 Li  2p+1:     1.01)                                                                                  
                                                                                                                          
               Molecular Orbital No.  11:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.13091 a.u.                                                                  
               (   1 Li  2s  :    -0.27) (   1 Li  3s  :     0.83) (   1 Li  1p0 :    -0.71)                              
               (   1 Li  2p0 :     1.33) (   2 H   2s  :    -0.31)                                                        
                                                                                                                          
               Molecular Orbital No.  12:                                                                                 
               --------------------------                                                                                 
               Occupation: 0.000 Energy:    0.20610 a.u.                                                                  
               (   1 Li  2s  :     1.69) (   1 Li  3s  :    -1.75) (   1 Li  1p0 :    -1.09)                              
               (   1 Li  2p0 :     0.67) (   2 H   2s  :     0.32)                                                        
                                                                                                                          
SCF converged in 12 iterations.
Total Energy: -83.8546945230 au
Sum of SCF:  
 -83.85469422644599 au
print("Sum of SCF:  ", e_scf_lih + e_scf_h2o, "au")
Sum of SCF:   -83.85469422644599 au
adc1_drv.nstates = 3 * nexc
adc1_results_lih_h2o = adc1_drv.compute(lih_h2o, basis_lih_h2o, scf_lih_h2o.scf_tensors)
adc2_drv.nstates = 3 * nexc
adc2_results_lih_h2o = adc2_drv.compute(lih_h2o, basis_lih_h2o, scf_lih_h2o.scf_tensors)
                                                                                                                          
                                                   ADC(1) Driver Setup                                                    
                                                  =====================                                                   
                                                                                                                          
                               Number Of Excited States  : 9                                                              
                               Max. Number Of Iterations : 50                                                             
                               Convergence Threshold     : 1.0e-05                                                        
                               ERI screening scheme      : Cauchy Schwarz + Density                                       
                               ERI Screening Threshold   : 1.0e-15                                                        
                                                                                                                          
* Info * Number of occupied orbitals: 7                                                                                   
* Info * Number of virtual orbitals: 17                                                                                   
                                                                                                                          
                                                MO Integrals Driver Setup                                                 
                                               ===========================                                                
                                                                                                                          
                               ERI screening scheme        : QQ_DEN                                                       
                               ERI Screening Threshold     : 1.0e-15                                                      
                               Batch Size of Fock Matrices : 100                                                          
                                                                                                                          
* Info * Processing Fock builds for the OO block...                                                                       
* Info *   batch 1/1                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the OO block done in 0.02 sec. Load imb.: 0.0 %                                     
                                                                                                                          
                   *** Iteration:   1 * Reduced Space:    9 * Residues (Max,Min): 2.56e-02 and 1.73e-05                   
                                                                                                                          
                   State  1: 0.15889350 au Residual Norm: 0.01952182                                                      
                   State  2: 0.20608352 au Residual Norm: 0.01468372                                                      
                   State  3: 0.20608352 au Residual Norm: 0.01468372                                                      
                   State  4: 0.28273172 au Residual Norm: 0.01773904                                                      
                   State  5: 0.33269457 au Residual Norm: 0.01820955                                                      
                   State  6: 0.33269457 au Residual Norm: 0.01820955                                                      
                   State  7: 0.33539541 au Residual Norm: 0.02560734                                                      
                   State  8: 0.51590645 au Residual Norm: 0.00006457                                                      
                   State  9: 0.52295746 au Residual Norm: 0.00001727                                                      
                                                                                                                          
                   *** Iteration:   2 * Reduced Space:   18 * Residues (Max,Min): 9.59e-04 and 3.22e-10                   
                                                                                                                          
                   State  1: 0.15870211 au Residual Norm: 0.00048111                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000013                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000016                                                      
                   State  4: 0.28256940 au Residual Norm: 0.00084006                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000019                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000027                                                      
                   State  7: 0.33476883 au Residual Norm: 0.00095895                                                      
                   State  8: 0.51590637 au Residual Norm: 0.00000200                                                      
                   State  9: 0.52295746 au Residual Norm: 0.00000000                                                      
                                                                                                                          
                   *** Iteration:   3 * Reduced Space:   27 * Residues (Max,Min): 1.28e-01 and 1.13e-08                   
                                                                                                                          
                   State  1: 0.15870199 au Residual Norm: 0.00001119                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000001                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000004                                                      
                   State  4: 0.28256867 au Residual Norm: 0.00002931                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000001                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000009                                                      
                   State  7: 0.33476780 au Residual Norm: 0.00004318                                                      
                   State  8: 0.36507782 au Residual Norm: 0.12766600                                                      
                   State  9: 0.51590637 au Residual Norm: 0.00000002                                                      
                                                                                                                          
                   *** Iteration:   4 * Reduced Space:   35 * Residues (Max,Min): 1.87e-01 and 2.18e-09                   
                                                                                                                          
                   State  1: 0.15870199 au Residual Norm: 0.00000001                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000001                                                      
                   State  4: 0.28256867 au Residual Norm: 0.00000001                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000003                                                      
                   State  7: 0.33476780 au Residual Norm: 0.00000001                                                      
                   State  8: 0.35283297 au Residual Norm: 0.02332168                                                      
                   State  9: 0.47917004 au Residual Norm: 0.18714435                                                      
                                                                                                                          
                   *** Iteration:   5 * Reduced Space:   44 * Residues (Max,Min): 2.33e-02 and 3.36e-11                   
                                                                                                                          
                   State  1: 0.15870199 au Residual Norm: 0.00000000                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  4: 0.28256867 au Residual Norm: 0.00000000                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  7: 0.33476780 au Residual Norm: 0.00000000                                                      
                   State  8: 0.35279733 au Residual Norm: 0.00006478                                                      
                   State  9: 0.44392307 au Residual Norm: 0.02334940                                                      
                                                                                                                          
                   *** Iteration:   6 * Reduced Space:   53 * Residues (Max,Min): 2.63e-04 and 1.00e-12                   
                                                                                                                          
                   State  1: 0.15870199 au Residual Norm: 0.00000000                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  4: 0.28256867 au Residual Norm: 0.00000000                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  7: 0.33476780 au Residual Norm: 0.00000000                                                      
                   State  8: 0.35279733 au Residual Norm: 0.00000209                                                      
                   State  9: 0.44360634 au Residual Norm: 0.00026329                                                      
                                                                                                                          
                   *** Iteration:   7 * Reduced Space:   58 * Residues (Max,Min): 6.63e-10 and 2.58e-16                   
                                                                                                                          
                   State  1: 0.15870199 au Residual Norm: 0.00000000                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  4: 0.28256867 au Residual Norm: 0.00000000                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  7: 0.33476780 au Residual Norm: 0.00000000                                                      
                   State  8: 0.35279733 au Residual Norm: 0.00000000                                                      
                   State  9: 0.44360633 au Residual Norm: 0.00000000                                                      
                                                                                                                          
               *** 9 excited states converged in 7 iterations. Time: 0.07 sec.                                            
                                                                                                                          
                                                                                                                          
               ADC(1) excited states                                                                                      
               ---------------------                                                                                      
               Excited State    S1:      0.15870199 a.u.      4.31850 eV    Osc.Str. 0.02843                              
               Excited State    S2:      0.20598021 a.u.      5.60501 eV    Osc.Str. 0.25169                              
               Excited State    S3:      0.20598021 a.u.      5.60501 eV    Osc.Str. 0.25169                              
               Excited State    S4:      0.28256867 a.u.      7.68909 eV    Osc.Str. 0.04341                              
               Excited State    S5:      0.33252612 a.u.      9.04850 eV    Osc.Str. 0.02652                              
               Excited State    S6:      0.33252612 a.u.      9.04850 eV    Osc.Str. 0.02652                              
               Excited State    S7:      0.33476780 a.u.      9.10950 eV    Osc.Str. 0.00002                              
               Excited State    S8:      0.35279733 a.u.      9.60010 eV    Osc.Str. 0.01663                              
               Excited State    S9:      0.44360633 a.u.     12.07114 eV    Osc.Str. 0.11240                              
                                                                                                                          
                                                                                                                          
                                                MO Integrals Driver Setup                                                 
                                               ===========================                                                
                                                                                                                          
                               ERI screening scheme        : QQ_DEN                                                       
                               ERI Screening Threshold     : 1.0e-15                                                      
                               Batch Size of Fock Matrices : 100                                                          
                                                                                                                          
* Info * Processing Fock builds for the OO block...                                                                       
* Info *   batch 1/1                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the OO block done in 0.02 sec. Load imb.: 0.0 %                                     
                                                                                                                          
* Info * Processing Fock builds for the VV block...                                                                       
* Info *   batch 1/2                                                                                                      
* Info *   batch 2/2                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the VV block done in 0.05 sec. Load imb.: 0.0 %                                     
                                                                                                                          
* Info * Processing Fock builds for the OV block...                                                                       
* Info *   batch 1/2                                                                                                      
* Info *   batch 2/2                                                                                                      
                                                                                                                          
* Info * Integrals transformation for the OV block done in 0.07 sec. Load imb.: 0.0 %                                     
                                                                                                                          
* Info * 1.17 MB of memory used for MO integrals on the master node.                                                      
* Info * 4.05 GB of memory available on the master node.                                                                  
* Info * Time spent in calculating memory usage: 0.01 sec.                                                                
                                                                                                                          
                                                                                                                          
                                                   ADC(2) Driver Setup                                                    
                                                  =====================                                                   
                                                                                                                          
                               Number Of Excited States  : 9                                                              
                               Max. Number Of Iterations : 150                                                            
                               Convergence Threshold     : 1.0e-05                                                        
                               ERI screening scheme      : Cauchy Schwarz + Density                                       
                               ERI Screening Threshold   : 1.0e-15                                                        
                                                                                                                          
                                                                                                                          
                                                   ADC(1) Driver Setup                                                    
                                                  =====================                                                   
                                                                                                                          
                               Number Of Excited States  : 9                                                              
                               Max. Number Of Iterations : 50                                                             
                               Convergence Threshold     : 1.0e-04                                                        
                               ERI screening scheme      : Cauchy Schwarz + Density                                       
                               ERI Screening Threshold   : 1.0e-15                                                        
                                                                                                                          
* Info * Number of occupied orbitals: 7                                                                                   
* Info * Number of virtual orbitals: 17                                                                                   
                                                                                                                          
                   *** Iteration:   1 * Reduced Space:    9 * Residues (Max,Min): 2.56e-02 and 1.73e-05                   
                                                                                                                          
                   State  1: 0.15889350 au Residual Norm: 0.01952182                                                      
                   State  2: 0.20608352 au Residual Norm: 0.01468372                                                      
                   State  3: 0.20608352 au Residual Norm: 0.01468372                                                      
                   State  4: 0.28273172 au Residual Norm: 0.01773904                                                      
                   State  5: 0.33269457 au Residual Norm: 0.01820955                                                      
                   State  6: 0.33269457 au Residual Norm: 0.01820955                                                      
                   State  7: 0.33539541 au Residual Norm: 0.02560734                                                      
                   State  8: 0.51590645 au Residual Norm: 0.00006457                                                      
                   State  9: 0.52295746 au Residual Norm: 0.00001727                                                      
                                                                                                                          
                   *** Iteration:   2 * Reduced Space:   18 * Residues (Max,Min): 9.59e-04 and 3.22e-10                   
                                                                                                                          
                   State  1: 0.15870211 au Residual Norm: 0.00048111                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000013                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000016                                                      
                   State  4: 0.28256940 au Residual Norm: 0.00084006                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000019                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000027                                                      
                   State  7: 0.33476883 au Residual Norm: 0.00095895                                                      
                   State  8: 0.51590637 au Residual Norm: 0.00000200                                                      
                   State  9: 0.52295746 au Residual Norm: 0.00000000                                                      
                                                                                                                          
                   *** Iteration:   3 * Reduced Space:   27 * Residues (Max,Min): 1.28e-01 and 1.13e-08                   
                                                                                                                          
                   State  1: 0.15870199 au Residual Norm: 0.00001119                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000001                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000004                                                      
                   State  4: 0.28256867 au Residual Norm: 0.00002931                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000001                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000009                                                      
                   State  7: 0.33476780 au Residual Norm: 0.00004318                                                      
                   State  8: 0.36507782 au Residual Norm: 0.12766600                                                      
                   State  9: 0.51590637 au Residual Norm: 0.00000002                                                      
                                                                                                                          
                   *** Iteration:   4 * Reduced Space:   35 * Residues (Max,Min): 1.87e-01 and 2.18e-09                   
                                                                                                                          
                   State  1: 0.15870199 au Residual Norm: 0.00000001                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000001                                                      
                   State  4: 0.28256867 au Residual Norm: 0.00000001                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000003                                                      
                   State  7: 0.33476780 au Residual Norm: 0.00000001                                                      
                   State  8: 0.35283297 au Residual Norm: 0.02332168                                                      
                   State  9: 0.47917004 au Residual Norm: 0.18714435                                                      
                                                                                                                          
                   *** Iteration:   5 * Reduced Space:   44 * Residues (Max,Min): 2.33e-02 and 3.36e-11                   
                                                                                                                          
                   State  1: 0.15870199 au Residual Norm: 0.00000000                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  4: 0.28256867 au Residual Norm: 0.00000000                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  7: 0.33476780 au Residual Norm: 0.00000000                                                      
                   State  8: 0.35279733 au Residual Norm: 0.00006478                                                      
                   State  9: 0.44392307 au Residual Norm: 0.02334941                                                      
                                                                                                                          
                   *** Iteration:   6 * Reduced Space:   53 * Residues (Max,Min): 2.63e-04 and 1.00e-12                   
                                                                                                                          
                   State  1: 0.15870199 au Residual Norm: 0.00000000                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  4: 0.28256867 au Residual Norm: 0.00000000                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  7: 0.33476780 au Residual Norm: 0.00000000                                                      
                   State  8: 0.35279733 au Residual Norm: 0.00000209                                                      
                   State  9: 0.44360634 au Residual Norm: 0.00026329                                                      
                                                                                                                          
                   *** Iteration:   7 * Reduced Space:   58 * Residues (Max,Min): 6.49e-10 and 1.88e-16                   
                                                                                                                          
                   State  1: 0.15870199 au Residual Norm: 0.00000000                                                      
                   State  2: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  3: 0.20598021 au Residual Norm: 0.00000000                                                      
                   State  4: 0.28256867 au Residual Norm: 0.00000000                                                      
                   State  5: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  6: 0.33252612 au Residual Norm: 0.00000000                                                      
                   State  7: 0.33476780 au Residual Norm: 0.00000000                                                      
                   State  8: 0.35279733 au Residual Norm: 0.00000000                                                      
                   State  9: 0.44360633 au Residual Norm: 0.00000000                                                      
                                                                                                                          
               *** 9 excited states converged in 7 iterations. Time: 0.05 sec.                                            
                                                                                                                          
                                                                                                                          
               ADC(1) excited states                                                                                      
               ---------------------                                                                                      
               Excited State    S1:      0.15870199 a.u.      4.31850 eV    Osc.Str. 0.02843                              
               Excited State    S2:      0.20598021 a.u.      5.60501 eV    Osc.Str. 0.25169                              
               Excited State    S3:      0.20598021 a.u.      5.60501 eV    Osc.Str. 0.25169                              
               Excited State    S4:      0.28256867 a.u.      7.68909 eV    Osc.Str. 0.04341                              
               Excited State    S5:      0.33252612 a.u.      9.04850 eV    Osc.Str. 0.02652                              
               Excited State    S6:      0.33252612 a.u.      9.04850 eV    Osc.Str. 0.02652                              
               Excited State    S7:      0.33476780 a.u.      9.10950 eV    Osc.Str. 0.00002                              
               Excited State    S8:      0.35279733 a.u.      9.60010 eV    Osc.Str. 0.01663                              
               Excited State    S9:      0.44360633 a.u.     12.07114 eV    Osc.Str. 0.11240                              
                                                                                                                          
                                                                                                                          
* Info * Number of occupied orbitals: 7                                                                                   
* Info * Number of virtual orbitals: 17                                                                                   
                                                                                                                          
* Info * MP2 correlation energy: -0.1411535484 a.u. Time: 0.00 sec.                                                       
                                                                                                                          
* Info * Initial guess generated in 0.05 sec.                                                                             
                                                                                                                          
                   *** Iteration:   1 * Reduced Space:    9 * Residual Norm: 1.80e-02                                     
                                                                                                                          
                   State  1:     0.15870199 au  Eigenvalue Error: 1.76e-02                                                
                                                                                                                          
                   *** Iteration:   2 * Reduced Space:   18 * Residual Norm: 1.57e-04                                     
                                                                                                                          
                   State  1:     0.14171087 au  Eigenvalue Error: 2.65e-05                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:   3 * Reduced Space:   18 * Residual Norm: 5.33e-06                                     
                                                                                                                          
                   State  1:     0.14168532 au  Eigenvalue Error: 5.76e-08                                                
                                                                                                                          
                   *** Iteration:   4 * Reduced Space:   27 * Residual Norm: 8.89e-08                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                                                                                                                          
* Info * State 1 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:   5 * Reduced Space:   34 * Residual Norm: 1.68e-03                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18797277 au  Eigenvalue Error: 1.66e-03                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:   6 * Reduced Space:   18 * Residual Norm: 3.55e-07                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637091 au  Eigenvalue Error: 3.55e-07                                                
                                                                                                                          
                   *** Iteration:   7 * Reduced Space:   27 * Residual Norm: 1.65e-09                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                                                                                                                          
* Info * State 2 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:   8 * Reduced Space:   34 * Residual Norm: 1.36e-08                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637055 au  Eigenvalue Error: 1.35e-08                                                
                                                                                                                          
                   *** Iteration:   9 * Reduced Space:   43 * Residual Norm: 1.65e-09                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                                                                                                                          
* Info * State 3 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:  10 * Reduced Space:   50 * Residual Norm: 2.81e-03                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26954997 au  Eigenvalue Error: 2.69e-03                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:  11 * Reduced Space:   18 * Residual Norm: 3.82e-05                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694932 au  Eigenvalue Error: 1.60e-06                                                
                                                                                                                          
                   *** Iteration:  12 * Reduced Space:   27 * Residual Norm: 3.64e-07                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                                                                                                                          
* Info * State 4 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:  13 * Reduced Space:   34 * Residual Norm: 2.10e-03                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31370842 au  Eigenvalue Error: 1.93e-03                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:  14 * Reduced Space:   18 * Residual Norm: 6.76e-05                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185716 au  Eigenvalue Error: 7.04e-07                                                
                                                                                                                          
                   *** Iteration:  15 * Reduced Space:   27 * Residual Norm: 1.86e-06                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  Eigenvalue Error: 3.91e-09                                                
                                                                                                                          
                   *** Iteration:  16 * Reduced Space:   34 * Residual Norm: 4.27e-07                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  converged                                                                 
                                                                                                                          
* Info * State 5 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:  17 * Reduced Space:   43 * Residual Norm: 2.69e-04                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  converged                                                                 
                   State  6:     0.31918451 au  Eigenvalue Error: 2.54e-04                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:  18 * Reduced Space:   18 * Residual Norm: 5.19e-06                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  converged                                                                 
                   State  6:     0.31893934 au  Eigenvalue Error: 3.87e-08                                                
                                                                                                                          
                   *** Iteration:  19 * Reduced Space:   27 * Residual Norm: 3.10e-08                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  converged                                                                 
                   State  6:     0.31893937 au  converged                                                                 
                                                                                                                          
* Info * State 6 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:  20 * Reduced Space:   34 * Residual Norm: 3.06e-05                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  converged                                                                 
                   State  6:     0.31893937 au  converged                                                                 
                   State  7:     0.31996155 au  Eigenvalue Error: 2.96e-05                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:  21 * Reduced Space:   18 * Residual Norm: 3.42e-10                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  converged                                                                 
                   State  6:     0.31893937 au  converged                                                                 
                   State  7:     0.31993282 au  converged                                                                 
                                                                                                                          
* Info * State 7 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:  22 * Reduced Space:   27 * Residual Norm: 3.98e-11                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  converged                                                                 
                   State  6:     0.31893937 au  converged                                                                 
                   State  7:     0.31993282 au  converged                                                                 
                   State  8:     0.31993282 au  converged                                                                 
                                                                                                                          
* Info * State 8 is converged.                                                                                            
                                                                                                                          
                   *** Iteration:  23 * Reduced Space:   34 * Residual Norm: 3.28e-03                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  converged                                                                 
                   State  6:     0.31893937 au  converged                                                                 
                   State  7:     0.31993282 au  converged                                                                 
                   State  8:     0.31993282 au  converged                                                                 
                   State  9:     0.40098055 au  Eigenvalue Error: 2.92e-03                                                
                                                                                                                          
* Info * Collapsing subspace...                                                                                           
                                                                                                                          
                   *** Iteration:  24 * Reduced Space:   18 * Residual Norm: 4.76e-05                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  converged                                                                 
                   State  6:     0.31893937 au  converged                                                                 
                   State  7:     0.31993282 au  converged                                                                 
                   State  8:     0.31993282 au  converged                                                                 
                   State  9:     0.39816796 au  Eigenvalue Error: 2.55e-06                                                
                                                                                                                          
                   *** Iteration:  25 * Reduced Space:   27 * Residual Norm: 8.86e-07                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  converged                                                                 
                   State  6:     0.31893937 au  converged                                                                 
                   State  7:     0.31993282 au  converged                                                                 
                   State  8:     0.31993282 au  converged                                                                 
                   State  9:     0.39816550 au  Eigenvalue Error: 2.96e-09                                                
                                                                                                                          
                   *** Iteration:  26 * Reduced Space:   34 * Residual Norm: 4.76e-08                                     
                                                                                                                          
                   State  1:     0.14168526 au  converged                                                                 
                   State  2:     0.18637056 au  converged                                                                 
                   State  3:     0.18637057 au  converged                                                                 
                   State  4:     0.26694777 au  converged                                                                 
                   State  5:     0.31185648 au  converged                                                                 
                   State  6:     0.31893937 au  converged                                                                 
                   State  7:     0.31993282 au  converged                                                                 
                   State  8:     0.31993282 au  converged                                                                 
                   State  9:     0.39816549 au  converged                                                                 
                                                                                                                          
* Info * State 9 is converged.                                                                                            
                                                                                                                          
* Info * Total number of sigma builds: 307                                                                                
* Info * Average time per sigma build: 0.01 sec                                                                           
                                                                                                                          
* Info * Total number of d_sigma builds: 26                                                                               
* Info * Average time per d_sigma build: 0.01 sec                                                                         
                                                                                                                          
               *** 9 excited states converged in 26 iterations.  Time: 4.47 sec.                                          
                                                                                                                          
                                                                                                                          
               ADC(2) excited states                                                                                      
               ---------------------                                                                                      
               Excited State    S1:      0.14168526 a.u.      3.85545 eV    |v1|^2=0.9645                                 
               Excited State    S2:      0.18637056 a.u.      5.07140 eV    |v1|^2=0.9634                                 
               Excited State    S3:      0.18637057 a.u.      5.07140 eV    |v1|^2=0.9634                                 
               Excited State    S4:      0.26694777 a.u.      7.26402 eV    |v1|^2=0.9652                                 
               Excited State    S5:      0.31185648 a.u.      8.48605 eV    |v1|^2=0.9595                                 
               Excited State    S6:      0.31893937 a.u.      8.67878 eV    |v1|^2=0.9662                                 
               Excited State    S7:      0.31993282 a.u.      8.70582 eV    |v1|^2=0.9719                                 
               Excited State    S8:      0.31993282 a.u.      8.70582 eV    |v1|^2=0.9719                                 
               Excited State    S9:      0.39816549 a.u.     10.83464 eV    |v1|^2=0.9638                                 
                                                                                                                          
               State S1                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO         0.9487                                                                            
               HOMO     -> LUMO+4      -0.2366                                                                            
               HOMO     -> LUMO+3      -0.0915                                                                            
                                                                                                                          
               State S2                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+1      -0.9433                                                                            
               HOMO     -> LUMO+6      -0.2710                                                                            
                                                                                                                          
               State S3                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+2      -0.9433                                                                            
               HOMO     -> LUMO+7       0.2710                                                                            
                                                                                                                          
               State S4                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+3       0.9756                                                                            
               HOMO     -> LUMO+4      -0.0876                                                                            
               HOMO     -> LUMO         0.0741                                                                            
                                                                                                                          
               State S5                                                                                                   
               -------------------------------                                                                            
               HOMO-1   -> LUMO+5      -0.9772                                                                            
               HOMO-1   -> LUMO+12      0.0541                                                                            
                                                                                                                          
               State S6                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+4       0.9512                                                                            
               HOMO     -> LUMO         0.2365                                                                            
               HOMO     -> LUMO+3       0.0707                                                                            
                                                                                                                          
               State S7                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+6      -0.9489                                                                            
               HOMO     -> LUMO+1       0.2672                                                                            
                                                                                                                          
               State S8                                                                                                   
               -------------------------------                                                                            
               HOMO     -> LUMO+7       0.9489                                                                            
               HOMO     -> LUMO+2       0.2672                                                                            
                                                                                                                          
               State S9                                                                                                   
               -------------------------------                                                                            
               HOMO-1   -> LUMO+8      -0.9765                                                                            
               HOMO-1   -> LUMO+14     -0.0824                                                                            
               HOMO-1   -> LUMO+10     -0.0588                                                                            
                                                                                                                          
                                                                                                                          
print("LiH ADC(1) excitation energies:\n", adc1_results_lih['eigenvalues'])
print("\nH2O ADC(1) excitation energies:\n", adc1_results_h2o['eigenvalues'])
print("\nLiH and H2O ADC(1) excitation energies:\n", adc1_results_lih_h2o['eigenvalues'])
print("\nLiH ADC(2) excitation energies:\n", adc2_results_lih['eigenvalues'])
print("\nH2O ADC(2) excitation energies:\n", adc2_results_h2o['eigenvalues'])
print("\nLiH and H2O ADC(2) excitation energies:\n", adc2_results_lih_h2o['eigenvalues'])
LiH ADC(1) excitation energies:
 [0.1587  0.20598 0.20598 0.28257 0.33253 0.33253]

H2O ADC(1) excitation energies:
 [0.3528  0.42544 0.44361]

LiH and H2O ADC(1) excitation energies:
 [0.1587  0.20598 0.20598 0.28257 0.33253 0.33253 0.33477 0.3528  0.44361]

LiH ADC(2) excitation energies:
 [0.14168 0.18637 0.18637 0.26695 0.31894 0.31993]

H2O ADC(2) excitation energies:
 [0.31186 0.39817 0.40289]

LiH and H2O ADC(2) excitation energies:
 [0.14169 0.18637 0.18637 0.26695 0.31186 0.31894 0.31993 0.31993 0.39817]

As we can see, the excitation energies of the individual LiH molecule occur also in the composite system, both at the ADC(1) and ADC(2) levels. At higher excitation energies, also the ones from water occur again. This means that the ADC(\(n\)) approaches are indeed size consistent (or size intensive). The interested reader can confirm the same finding also for the transition moments or oscillator strengths.

Note

Since ADC(1) is in fact equivalent to configuration interaction singles (CIS) for excitation energies and excited-state properties, this means that CIS is the only truncated CI approach that is size consistent.