loos/CBD
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

typos and others for Enzo

Browse Source
This commit is contained in:
Pierre-Francois Loos 2022-03-14 15:29:31 +01:00
parent 2be233eea6
commit 363127773e
1 changed files with 33 additions and 35 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

68
Manuscript/CBD.tex
View File

@ -199,6 +199,7 @@
\title{Reference energies for cyclobutadiene}
\author{Enzo \surname{Monino}}
\email{emonino@irsamc.ups-tlse.fr}
\affiliation{\LCPQ}
\author{Martial \surname{Boggio-Pasqua}}
\affiliation{\LCPQ}
@ -211,7 +212,7 @@
\affiliation{\LCPQ}
\begin{abstract}
The cyclobutadiene (CBD) molecule represents a playground for ground state and excited states methods. Indeed due to the high symmetry of the molecule, especially at the square geometry ($D_{4h}$) but also at the rectangular structure ($D_{2h}$), the ground state and the excited states of the CBD exhibit multiconfigurational character where single reference methods such as adiabatic time-dependent density functional theory (TD-DFT) or equation-of-motion coupled cluster (EOM-CC) show difficulty to describe them. In this work we provide an extensive study of the autoisomerization barrier (AB), where the rectangular ($D_{2h}$) and the square geometry ($D_{4h}$) are needed, and of the vertical excitations energies of the CBD molecule using a large range of methods and basis set. In order to tackle the problem of multiconfigurational character presents in the AB and in the vertical excitation energies selected configuration interaction (SCI) and multi reference (CASSCF,CASPT2,NEVPT2) calculations are performed. Moreover coupled cluster calculations such as CCSD, CC3, CCSDT, CC4 and CCSDTQ are added to the set of methods. To complete the study we provide spin-flip (SF) results, which are known to give correct description of multiconfigurational character states, in the TD-DFT framework where numerous exchange-correlation functionals are considered, we also add algebraic diagrammatic construction (ADC) calculations in the SF formalism where we use the ADC(2)-s, ADC(2)-x and ADC(3) schemes. A theoretical best estimate (TBE) is given for the AB and for each vertical energies.
The cyclobutadiene molecule represents a playground for ground state and excited states methods. Indeed, due to the high symmetry of the molecule, especially at the square geometry ($D_{4h}$) but also at the rectangular structure ($D_{2h}$), the ground state and the excited states of cyclobutadiene exhibit multiconfigurational character where single reference methods such as adiabatic time-dependent density functional theory (TD-DFT) or equation-of-motion coupled cluster (EOM-CC) show difficulty to describe them. In this work we provide an extensive study of the autoisomerization barrier, where the rectangular ($D_{2h}$) and the square geometry ($D_{4h}$) are needed, and of the vertical excitations energies of cyclobutadiene using a large range of methods and basis sets. In order to tackle the problem of multiconfigurational character presents in the autoisomerization barrier and in the vertical excitation energies selected configuration interaction and multireference (CASSCF, CASPT2, and NEVPT2) calculations are performed. Moreover coupled cluster calculations such as CCSD, CC3, CCSDT, CC4 and CCSDTQ are added to the set of methods. To complete the study we provide spin-flip results, which are known to give correct description of multiconfigurational character states, in the TD-DFT framework where numerous exchange-correlation functionals are considered, we also add algebraic diagrammatic construction (ADC) calculations in the spin-flip formalism where we use the ADC(2)-s, ADC(2)-x and ADC(3) schemes. A theoretical best estimate is defined for the autoisomerization barrier and for each vertical energies.
\end{abstract}
\maketitle
@ -221,7 +222,7 @@ The cyclobutadiene (CBD) molecule represents a playground for ground state and e
\label{sec:intro}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Despite the fact that excited states are involved in ubiquitious processes such as photochemistry, catalysis or in solar cell technology, none of the many methods existing is the reference in providing accurate excitation energies. Indeed, each method has its own flaws and there are so many chemical scenario that can occur, so it is still one of the biggest challenge in theoretical chemistry. Speaking of difficult task, cyclobutadiene (CBD) molecule has been a real challenge for experimental and theoretical chemists for many decades \cite{bally_1980}. Due to his antiaromaticity \cite{minkin_1994} and his large angular strain \cite{baeyer_1885} the CBD molecule presents a high reactivity which made the synthesis of this molecule a particularly difficult exercise. Hückel molecular orbital theory gives a triplet state with square ($D_{4h}$) geometry for the ground state of the CBD,with the two singly occupied frontier orbitals that are degenerated by symmetry. This degeneracy is lifted by the Jahn-Teller effect, meaning by distortion of the molecule (lowering symmetry), and gives a singlet state with rectangular ($D_{2h}$) geometry for the ground state.
Despite the fact that excited states are involved in ubiquitous processes such as photochemistry, catalysis or in solar cell technology, none of the many methods existing is the reference in providing accurate excitation energies. Indeed, each method has its own flaws and there are so many chemical scenario that can occur, so it is still one of the biggest challenge in theoretical chemistry. Speaking of difficult task, cyclobutadiene (CBD) molecule has been a real challenge for experimental and theoretical chemists for many decades \cite{bally_1980}. Due to his antiaromaticity \cite{minkin_1994} and his large angular strain \cite{baeyer_1885} the CBD molecule presents a high reactivity which made the synthesis of this molecule a particularly difficult exercise. H\"uckel molecular orbital theory gives a triplet state with square ($D_{4h}$) geometry for the ground state of the CBD, with the two singly occupied frontier orbitals that are degenerated by symmetry. This degeneracy is lifted by the Jahn-Teller effect, meaning by distortion of the molecule (lowering symmetry), and gives a singlet state with rectangular ($D_{2h}$) geometry for the ground state.
Indeed, synthetic work from Pettis and co-workers \cite{reeves_1969} gives a rectangular geometry to the singlet ground state of CBD and then was confirmed by experimental works \cite{irngartinger_1983,ermer_1983,kreile_1986}.
At the ground state structrure ($D_{2h}$), the ${}^1A_g$ state has a weak multi-configurational character because of the well separated frontier orbitals and can be described by single-reference methods. But at the square ($D_{4h}$) geometry, the singlet state ${}^1B_{1g}$ has two singly occupied frontier orbitals that are degenerated so has a two-configurational character and single-reference methods are unreliable to describe it. The singlet ($D_{4h}$) is a transition state in the automerization reaction between the two rectangular structures (see Fig.\ref{fig:CBD}). The autoisomerization barrier (AB) for the CBD molecule is defined as the energy difference between the singlet ground state of the square ($D_{4h}$) structure and the singlet ground state of the rectangular ($D_{2h}$) geometry. The energy of this barrier was predicted, experimentally, in the range of 1.6-10 kcal.mol$^{-1}$ \cite{whitman_1982} and multi-reference calculations gave an energy barrier in the range of 6-7 kcal.mol$^{-1}$ \cite{eckert-maksic_2006}. All the specificities of the CBD molecule make it a real playground for excited-states methods.
@ -234,7 +235,7 @@ Another way to deal with double excitations is to use high level truncation of t
An alternative to multiconfigurational and CC methods is the use of selected CI (SCI) methods for the computation of transition energies for singly and doubly excited states that are known to reach near full CI energies for small molecules. These methods allow to avoid an exponential increase of the size of the CI expansion by retaining the most energetically relevant determinants only, using a second-order energetic criterion to select perturbatively determinants in the FCI space.
Finally, to describe doubly excited states, one can think of spin-flip formalism established by Krylov in 2001 \cite{casanova_2020}. To briefly introduce the spin-flip idea we can present it like: instead of taking the singlet ground state as reference, the reference configuration is taken as the lowest triplet state. So one can access the singlet ground state and the singlet doubly-excited state via a spin-flip deexcitation and excitation (respectively), the difference of these two excitation energies providing an estimate of the double excitation. Obviously spin-flip methods have their own flaws, especially the spin contamination \cite{casanova_2020} (i.e., an artificial mixing of electronic states of differents spin multiplicities) due to spin incompleteness of the spin-flip expansion and by spin contamination of the reference configuration. One can adress part of this problem by expansion of the excitation order but with an increase of the computational cost or by complementing the spin-incomplete configuration set with the missing configurations.
Finally, to describe doubly excited states, one can think of spin-flip formalism established by Krylov in 2001 \cite{casanova_2020}. To briefly introduce the spin-flip idea we can present it like: instead of taking the singlet ground state as reference, the reference configuration is taken as the lowest triplet state. So one can access the singlet ground state and the singlet doubly-excited state via a spin-flip deexcitation and excitation (respectively), the difference of these two excitation energies providing an estimate of the double excitation. Obviously spin-flip methods have their own flaws, especially the spin contamination \cite{casanova_2020} (i.e., an artificial mixing of electronic states of different spin multiplicities) due to spin incompleteness of the spin-flip expansion and by spin contamination of the reference configuration. One can address part of this problem by expansion of the excitation order but with an increase of the computational cost or by complementing the spin-incomplete configuration set with the missing configurations.
In the present work we investigate ${}^1A_g$, $1{}^3B_{1g}$, $1{}^1B_{1g}$, $2{}^1A_{g}$ and ${}^1B_{1g}$, $1{}^3A_{2g}$, $2{}^1A_{1g}$,$1{}^1B_{2g}$ excited states for the $D_{2h}$ and $D_{4h}$ geometries, respectively. Computational details are reported in Section \ref{sec:compmet} for SCI (Subsection \ref{sec:SCI}), EOM-CC (Subsection \ref{sec:CC}), multiconfigurational (Subsection \ref{sec:Multi}) and spin-flip (Subsection \ref{sec:sf}) methods. Section \ref{sec:res} is devoted to the discussion of our results, first we consider the ground state property studied which is the AB (Subsection \ref{sec:auto}) and then we study the excited states (Subsection \ref{sec:states}) of the CBD molecule for both geometries.
@ -245,40 +246,40 @@ In the present work we investigate ${}^1A_g$, $1{}^3B_{1g}$, $1{}^1B_{1g}$, $2{}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Computational Details}
\section{Computational details}
\label{sec:compmet}
%The system under investigation in this work is the cyclobutadiene (CBD) molecule, rectangular ($D_{2h}$) and square ($D_{4h}$) geometries are considered. The ($D_{2h}$) geometry is obtained at the CC3 level without frozen core using the aug-cc-pVTZ and the ($D_{4h}$) geometry is obtained at the RO-CCSD(T) level using aug-cc-pVTZ again without frozen core. All the calculations are performed using four basis set, the 6-31+G(d) basis and the aug-cc-pVXZ with X$=$D, T, Q. In the following we will use the notation AVXZ for the aug-cc-pVXZ basis, again with X$=$D, T, Q. The $\%T_1$ metric that gives the percentage of single excitation calculated at the CC3 level in \textcolor{red}{DALTON} allows to characterize the various states.Throughout all this work, spin-flip and spin-conserved calculations are performed with a UHF reference.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Selected Configuration Interaction}
\subsection{Selected configuration interaction calculations}
\label{sec:SCI}
States energies and excitations energies calculations in the SCI framework are performed using QUANTUM PACKAGE \cite{garniron_2019} where the CIPSI algorithm is implemented. The CIPSI algorithm allows to avoid the exponential increase of the CI expansion. To treat electronic states in the same way we use a state-averaged formalism meaning that the ground and excited states are represented with the same number and same set of determinants but using different CI coefficients. Then the SCI energy is the sum of two terms, the variational energy obtained by diagonalization of the CI matrix in the reference space and a second-order perturbative correction which estimates the contribution of the determinants not included in the CI space (estimate error in the truncation). It is possible to estimate the FCI limit for the total energies and compute the corresponding transition energies by extrapolating this second-order correction to zero. Extrapolation brings error and we can estime this one by energy difference between excitation energies obtained with the largest SCI wave function and the FCI extrapolated value. These errors provide a rough idea of the quality of the FCI extrapolation and cannot be seen as true bar error, they are reported in the following tables.
States energies and excitations energies calculations in the SCI framework are performed using QUANTUM PACKAGE \cite{garniron_2019} where the CIPSI algorithm is implemented. The CIPSI algorithm allows to avoid the exponential increase of the CI expansion. To treat electronic states in the same way we use a state-averaged formalism meaning that the ground and excited states are represented with the same number and same set of determinants but using different CI coefficients. Then the SCI energy is the sum of two terms, the variational energy obtained by diagonalization of the CI matrix in the reference space and a second-order perturbative correction which estimates the contribution of the determinants not included in the CI space (estimate error in the truncation). It is possible to estimate the FCI limit for the total energies and compute the corresponding transition energies by extrapolating this second-order correction to zero. Extrapolation brings error and we can estimate this one by energy difference between excitation energies obtained with the largest SCI wave function and the FCI extrapolated value. These errors provide a rough idea of the quality of the FCI extrapolation and cannot be seen as true bar error, they are reported in the following tables.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Coupled-Cluster}
\subsection{Coupled-Cluster calculations}
\label{sec:CC}
Different flavours of coupled-cluster (CC) calculations are performed using different codes. Indeed, CC theory provides a hierarchy of methods that provide increasingly accurate energies via the increase of the maximum excitation degree of the cluster operator. Without any truncation of the cluster operator one has the full CC (FCC) that is equivalent to the full configuration interaction (FCI) giving the exact energy and wave function of the system for a fixed atomic basis set. However, due to the computational exponential scaling with system size we have to use truncated CC methods. The CC with singles and doubles (CCSD), CC with singles, doubles and triples (CCSDT) calculations are achieved with \textcolor{red}{CFOUR}. The calculations in the context of CC response theory or ``approximate'' series (CC3,CC4) are performed with \textcolor{red}{DALTON} \cite{aidas_2014}. The CC with singles, doubles, triples and quadruples (CCSDTQ) are done with the \textcolor{red}{CFOUR} code. The CC2 \cite{christiansen_1995a,hattig_2000}, CC3 \cite{christiansen_1995b,koch_1995} and CC4 \cite{kallay_2005,matthews_2020} methods can be seen as cheaper approximations of CCSD \cite{purvis_1982}, CCSDT \cite{noga_1987} and CCSDTQ \cite{kucharski_1991a} by skipping the most expensive terms and avoiding the storage of higher-excitations amplitudes.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Multiconfigurational methods}
\subsection{Multiconfigurational calculations}
\label{sec:Multi}
State-averaged complete-active-space self-consistent field (SA-CASSCF) calculations are performed with \textcolor{red}{MOLPRO} \cite{werner_2012}. On top of those, NEVPT2 calculations, both partially contracted (PC) and strongly contracted (SC) scheme are considered. The PC-NEVPT2 is theoretically more accurate to the strongly contracted version due to the larger number of pertubers and greater flexibility. CASPT2 is performed and extended multistate (XMS) CASPT2 for strong mixing between states with same spin and spatial symmetries is also performed.
State-averaged complete-active-space self-consistent field (SA-CASSCF) calculations are performed with \textcolor{red}{MOLPRO} \cite{werner_2012}. On top of those, NEVPT2 calculations, both partially contracted (PC) and strongly contracted (SC) scheme are considered. The PC-NEVPT2 is theoretically more accurate to the strongly contracted version due to the larger number of perturbers and greater flexibility. CASPT2 is performed and extended multistate (XMS) CASPT2 for strong mixing between states with same spin and spatial symmetries is also performed.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Spin-Flip}
\subsection{Spin-flip calculations}
\label{sec:sf}
In both structures the CBD has a singlet ground state, for the spin-flip calculations we consider the lowest triplet state as reference. Spin-flip techniques are broadly accessible and here, among them, we explore algebraic-diagrammatic construction \cite{schirmer_1982} (ADC) using standard ADC(2)-s \cite{trofimov_1997,dreuw_2015} and extended ADC(2)-x \cite{dreuw_2015} schemes as well as the ADC(3) \cite{dreuw_2015,trofimov_2002,harbach_2014} scheme. We also use spin-flip within the TD-DFT \cite{casida_1995} framework. The standard and extended spin-flip ADC(2) (SF-ADC(2)-s and SF-ADC(2)-x respectively) and SF-ADC(3) are performed using Q-CHEM 5.2.1 \cite{shao_2015}. Spin-flip TD-DFT calculations are also performed using Q-CHEM 5.2.1. The B3LYP \cite{becke_1988b,lee_1988a,becke_1993b}, PBE0 \cite{adamo_1999a,ernzerhof_1999} and BH\&HLYP hybrid GGA functionals are considered, they contain 20\%, 25\%, 50\% of exact exchange and are labeled, respectively, as SF-BLYP, SF-B3LYP, SF-PBE0, SF-BH\&HLYP. We also have done spin-flip TD-DFT calculations using range-separated hybrid (RSH) functionals as: CAM-B3LYP \cite{yanai_2004a}, LC-$\omega$PBE08 \cite{weintraub_2009a} and $\omega$B97X-V \cite{mardirossian_2014}. The main difference here between these RSH functionals is the amount of exact-exchange at long-range: 75$\%$ for CAM-B3LYP and 100$\%$ for LC-$\omega$PBE08 and $\omega$B97X-V. To complete the use of spin-flip TD-DFT we also considered the hybrid meta-GGA functional M06-2X \cite{zhao_2008} and the RSH meta-GGA functional M11 \cite{peverati_2011}. Note that all SF-TD-DFT calculations are done within the TDA approximation.
%EOM-SF-CCSD and EOM-SF-CC(2,3) are also performed with Q-CHEM 5.2.1.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Theoretical Best Estimate (TBE)}
All the calculations are performed using four basis set, the 6-31+G(d) basis and the aug-cc-pVXZ with X$=$D, T, Q \cite{dunning_1989}. In the following we will use the notation AVXZ for the aug-cc-pVXZ basis, again with X$=$D, T, Q. For each studied quantity, i.e., the autoisomerisation barrier and the vertical excitations, we provide a theoretical best estimate (TBE). These TBEs are provided using extrapolated CCSDTQ/AVTZ values when possible and using NEVPT2(12,12) otherwise. The extrapolation of the CCSDTQ/AVTZ values is done using two schemes. The first one uses CC4 values for the extrapolation and proceed as follows
\subsection{Theoretical best estimates}
All the calculations are performed using four basis set, the 6-31+G(d) basis and the aug-cc-pVXZ with X $=$ D, T, Q \cite{dunning_1989}. In the following we will use the notation AVXZ for the aug-cc-pVXZ basis, again with X $=$ D, T, Q. For each studied quantity, i.e., the autoisomerisation barrier and the vertical excitations, we provide a theoretical best estimate (TBE). These TBEs are provided using extrapolated CCSDTQ/AVTZ values when possible and using NEVPT2(12,12) otherwise. The extrapolation of the CCSDTQ/AVTZ values is done using two schemes. The first one uses CC4 values for the extrapolation and proceed as follows
\begin{equation}
\label{eq:AVTZ}
@ -315,12 +316,12 @@ Two different sets of geometries obtained with different level of theory are con
%%% TABLE I %%%
\begin{squeezetable}
\begin{table}
\caption{Optimized geometries of the CBD molecule. Bond lenghts are in angstr\"om and angles are in degree.}
\caption{Optimized geometries of the CBD molecule. Bond lengths are in \si{\angstrom} and angles are in degree.}
\label{tab:geometries}
\begin{ruledtabular}
\begin{tabular}{llrrr}
Method & C=C & C-C & C-H & H-C=C\fnm[1] \\
Method & \ce{C=C} & \ce{C-C} & \ce{C-H} & \ce{H-C=C}\fnm[1] \\
\hline
$D_{2h}$\\
\hline
@ -338,8 +339,8 @@ Two different sets of geometries obtained with different level of theory are con
RO-CCSD(T)/AVTZ & 1.439 & 1.439 & 1.075 & 135.00
\end{tabular}
\end{ruledtabular}
\fnt[1]{Angle between the C-H bond and the C=C bond.}
\fnt[2]{From reference \cite{manohar_2008}}
\fnt[1]{Angle between the \ce{C-H} bond and the \ce{C=C} bond.}
\fnt[2]{From Ref.~\onlinecite{manohar_2008}.}
\end{table}
\end{squeezetable}
%%% %%% %%% %%%
@ -350,7 +351,7 @@ Two different sets of geometries obtained with different level of theory are con
%================================================
\subsection{Autoisomerization barrier}
\label{sec:auto}
The results for the calculation of the automerization barrier energy with and without spin-flip methods are shown in Table \ref{tab:auto_standard}. Two types of methods are used with spin-flip, SF-TD-DFT with several functionals and SF-ADC with the variants SF-ADC(2)-s, SF-ADC(2)-x and SF-ADC(3). First, one can see that there is large variations of the energy between the different functionals. Indeed if we look at the hybrid functionals B3LYP and BH\&HLYP we can see that there is a difference of around 7 \kcalmol throught all the bases, the difference in energy between the B3LYP and PBE0 functionals is much smaller with around 1.5 \kcalmol throught all the bases. We find a similar behaviour regarding the RSH functionals, we find a difference of about 8-9 \kcalmol between the M06-2X and the CAM-B3LYP functionals for all bases. The results between the CAM-B3LYP and the $\omega$B97X-V are very close in energy with a difference of around 0.1-0.2 \kcalmol . The energy difference between the M11 and the M06-2X functionals is larger with 0.6-0.9 \kcalmol for the AVXZ bases and with 1.7 \kcalmol for the 6-31+G(d) basis. For the SF-ADC methods the energy differences are smaller with 1.7-2.0 \kcalmol between the ADC(2)-s and the ADC(2)-x schemes, 0.9-1.6 \kcalmol between the ADC(2)-s and the ADC(3) schemes and 0.4-0.8 \kcalmol between the ADC(2)-x and the ADC(3) schemes.
The results for the calculation of the automerization barrier energy with and without spin-flip methods are shown in Table \ref{tab:auto_standard}. Two types of methods are used with spin-flip, SF-TD-DFT with several functionals and SF-ADC with the variants SF-ADC(2)-s, SF-ADC(2)-x and SF-ADC(3). First, one can see that there is large variations of the energy between the different functionals. Indeed if we look at the hybrid functionals B3LYP and BH\&HLYP we can see that there is a difference of around 7 \kcalmol through all the bases, the difference in energy between the B3LYP and PBE0 functionals is much smaller with around 1.5 \kcalmol through all the bases. We find a similar behavior regarding the RSH functionals, we find a difference of about 8-9 \kcalmol between the M06-2X and the CAM-B3LYP functionals for all bases. The results between the CAM-B3LYP and the $\omega$B97X-V are very close in energy with a difference of around 0.1-0.2 \kcalmol . The energy difference between the M11 and the M06-2X functionals is larger with 0.6-0.9 \kcalmol for the AVXZ bases and with 1.7 \kcalmol for the 6-31+G(d) basis. For the SF-ADC methods the energy differences are smaller with 1.7-2.0 \kcalmol between the ADC(2)-s and the ADC(2)-x schemes, 0.9-1.6 \kcalmol between the ADC(2)-s and the ADC(3) schemes and 0.4-0.8 \kcalmol between the ADC(2)-x and the ADC(3) schemes.
Then we compare results for multireference methods, we can see a difference of about 2.91-3.22 \kcalmol through all the bases between the CASSCF(12,12) and the CASPT2(12,12) methods. These differences can be explained by the well known lack of dynamical correlation for the CASSCF method that is typically important for close-shell molecules which is the case for the ground states of the rectangular and square geometries of the CBD molecule. The results between the CASPT2(12,12) and the NEVPT2(12,12) are much closer with an energy difference of around 0.12-0.23 \kcalmol for all the bases. Finally the last results shown in Table \ref{tab:auto_standard} are the CC ones, for the autoisomerization barrier energy we consider the CCSD, CCSDT, CCSDTQ methods and the approximations of CCDT and of CCSDTQ, the CC3 and the CC4 methods, respectively. We can see that the CCSD values are higher than the other CC methods with an energy difference of around 1.05-1.24 \kcalmol between the CCSD and the CCSDT methods. The CCSDT and CCSDTQ autoisomerization barrier energies are closer with 0.25 \kcalmol of energy difference. The energy difference between the CCSDT and its approximation CC3 is about 0.67-0.8 \kcalmol for all the bases whereas the energy difference between the CCSDTQ and its approximate version CC4 is 0.11 \kcalmol.
%%% TABLE I %%%
@ -385,8 +386,6 @@ CCSDT & $7.26$ & $7.64$ & $8.68$ &$\left[ 8.86\right]$\fnm[1] \\
CC4 & $7.40$ & $7.78$ & $\left[ 8.82\right]$\fnm[2] & $\left[ 9.00\right]$\fnm[3]\\
CCSDTQ & $7.51$ & $\left[ 7.89\right]$\fnm[4]& $\left[ 8.93\right]$\fnm[5]& $\left[ 9.11\right]$\fnm[6]\\
CIPSI & $7.91\pm 0.21$ & $8.58\pm 0.14$ & & \\
\end{tabular}
\end{ruledtabular}
\fnt[1]{Value obtained using CCSDT/AVTZ corrected by the difference between CC3/AVQZ and CC3/AVTZ.}
@ -419,13 +418,13 @@ Figure \ref{fig:AB} shows the autoisomerization barrier (AB) energies for the CB
\label{sec:states}
%All the calculations are performed using four basis set, the 6-31+G(d) basis and the aug-cc-pVXZ with X$=$D, T, Q. In the following we will use the notation AVXZ for the aug-cc-pVXZ basis, again with X$=$D, T, Q.
\subsubsection{D2h geometry}
Table \ref{tab:sf_tddft_D2h} shows the results obtained for the vertical transition energies using spin-flip methods and Table \ref{tab:D2h} shows the results obtained with the standard methods. We discuss first the SF-TD-DFT values with hybrid functionals. For the B3LYP functional we can see that the energy differences for each state and throughout all bases are small with the largest one for the $1\,{}^3B_{1g}$ state with 0.012 eV. The same observation can be done for the PBE0 and BH\&HLYP functionals, we can also oberve that adding exact exchange to the functional (20\% of exact exchange for B3LYP, 25\% for PBE0 and 50\% for BH\&HLYP increase the energy difference between the $1\,{}^1B_{1g}$ and the $1\,{}^3B_{1g}$ states. Put another way, increasing the amount of exact exchange in the functional reduces the energy difference between the $1\,{}^3B_{1g}$ and the $\text{X}\,{}^1A_{g}$ states. We can also notice that the vertical energies of the different states do not vary in the same way when adding exact exchange, for instance the energy variation for the $1\,{}^3B_{1g}$ state from PBE0 to BH\&HLYP is around 0.1 eV whereas for the $1\,{}^1B_{1g}$ and the $2\,{}^1A_{1g}$ states this energy variation is about 0.4-0.5 eV and 0.34 eV respectively. For the RSH functionals we can not make the same observation, indeed we can see that for the CAM-B3LYP functional we have that the energy difference between the $1\,{}^1B_{1g}$ and the $1\,{}^3B_{1g}$ states is larger than for the $\omega$B97X-V and LC-$\omega $PBE08 functionals despite the fact that the latter ones have a bigger amount of exact exchange. However we can observe that we have a small energy difference for each state throughout the bases and for each RSH functional. The M06-2X functional is an Hybrid meta-GGA functional and contains 54\% of Hartree-Fock exchange, we can compare the energies with the BH\&HLYP functional and we can see that the energy differences are small with around 0.03-0.08 eV. We can notice that the upper bound of 0.08 eV in the energy differences is due to the $1\,{}^3B_{1g}$ state. The M11 vertical energies are close to the BH\&HLYP ones for the triplet and the first singlet excited states with 0.01-0.02 eV and 0.08-0.10 eV of energy difference, respectively. For the $2\,{}^1A_{1g}$ state the M11 energies are closer to the $\omega$B97X-V ones with 0.05-0.09 eV of energy difference.
\subsubsection{$D_{2h}$ geometry}
Table \ref{tab:sf_tddft_D2h} shows the results obtained for the vertical transition energies using spin-flip methods and Table \ref{tab:D2h} shows the results obtained with the standard methods. We discuss first the SF-TD-DFT values with hybrid functionals. For the B3LYP functional we can see that the energy differences for each state and throughout all bases are small with the largest one for the $1\,{}^3B_{1g}$ state with 0.012 eV. The same observation can be done for the PBE0 and BH\&HLYP functionals, we can also observe that adding exact exchange to the functional (20\% of exact exchange for B3LYP, 25\% for PBE0 and 50\% for BH\&HLYP increase the energy difference between the $1\,{}^1B_{1g}$ and the $1\,{}^3B_{1g}$ states. Put another way, increasing the amount of exact exchange in the functional reduces the energy difference between the $1\,{}^3B_{1g}$ and the $\text{X}\,{}^1A_{g}$ states. We can also notice that the vertical energies of the different states do not vary in the same way when adding exact exchange, for instance the energy variation for the $1\,{}^3B_{1g}$ state from PBE0 to BH\&HLYP is around 0.1 eV whereas for the $1\,{}^1B_{1g}$ and the $2\,{}^1A_{1g}$ states this energy variation is about 0.4-0.5 eV and 0.34 eV respectively. For the RSH functionals we can not make the same observation, indeed we can see that for the CAM-B3LYP functional we have that the energy difference between the $1\,{}^1B_{1g}$ and the $1\,{}^3B_{1g}$ states is larger than for the $\omega$B97X-V and LC-$\omega $PBE08 functionals despite the fact that the latter ones have a bigger amount of exact exchange. However we can observe that we have a small energy difference for each state throughout the bases and for each RSH functional. The M06-2X functional is an Hybrid meta-GGA functional and contains 54\% of Hartree-Fock exchange, we can compare the energies with the BH\&HLYP functional and we can see that the energy differences are small with around 0.03-0.08 eV. We can notice that the upper bound of 0.08 eV in the energy differences is due to the $1\,{}^3B_{1g}$ state. The M11 vertical energies are close to the BH\&HLYP ones for the triplet and the first singlet excited states with 0.01-0.02 eV and 0.08-0.10 eV of energy difference, respectively. For the $2\,{}^1A_{1g}$ state the M11 energies are closer to the $\omega$B97X-V ones with 0.05-0.09 eV of energy difference.
Then we discuss the various ADC scheme (ADC(2)-s, ADC(2)-x and ADC(3)) results. For ADC(2) we have vertical energy differencies of about 0.03 eV for the $1\,{}^3B_{1g}$ state and around 0.06 eV for $2\,{}^1A_{1g}$ state throughout all bases. However for the $1\,{}^1B_{1g}$ state we have an energy difference of about 0.2 eV. For the ADC(2)-x and ADC(3) schemes the calculations with the AVQZ basis are not feasible with our resources. With ADC(2)-x we have similar vertical energies for the triplet and the $1\,{}^1B_{1g}$ states but for the $2\,{}^1A_{1g}$ state the energy difference between the ADC(2) and ADC(2)-x schemes is about 0.4-0.5 eV. The ADC(3) values are closer from the ADC(2) than the ADC(2)-x except for the triplet state for which we have a energy difference of about 0.09-0.14 eV. Now, we look at Table \ref{tab:D2h} to discuss the results of standard methods. First we discuss the CC values, we have computed the vertical energies at the CCSDT and CCSDTQ level as well as their approximation the CC3 and CC4 methods, respectively. We can notice that for the $1\,{}^3B_{1g}$ and the $1\,{}^1B_{1g}$ states the CCSDT and the CC3 values are close with an energy difference of 0.009-0.02 eV for all bases. The energy difference is larger for the $2\,{}^1A_{1g}$ state with around 0.35-0.38 eV. The same observation can be done for CCSDTQ and CC4 with similar vertical energies for all bases. Note that the $1\,{}^3B_{1g}$ state can not be described with these methods.
Then we discuss the various ADC scheme (ADC(2)-s, ADC(2)-x and ADC(3)) results. For ADC(2) we have vertical energy differences of about 0.03 eV for the $1\,{}^3B_{1g}$ state and around 0.06 eV for $2\,{}^1A_{1g}$ state throughout all bases. However for the $1\,{}^1B_{1g}$ state we have an energy difference of about 0.2 eV. For the ADC(2)-x and ADC(3) schemes the calculations with the AVQZ basis are not feasible with our resources. With ADC(2)-x we have similar vertical energies for the triplet and the $1\,{}^1B_{1g}$ states but for the $2\,{}^1A_{1g}$ state the energy difference between the ADC(2) and ADC(2)-x schemes is about 0.4-0.5 eV. The ADC(3) values are closer from the ADC(2) than the ADC(2)-x except for the triplet state for which we have a energy difference of about 0.09-0.14 eV. Now, we look at Table \ref{tab:D2h} to discuss the results of standard methods. First we discuss the CC values, we have computed the vertical energies at the CCSDT and CCSDTQ level as well as their approximation the CC3 and CC4 methods, respectively. We can notice that for the $1\,{}^3B_{1g}$ and the $1\,{}^1B_{1g}$ states the CCSDT and the CC3 values are close with an energy difference of 0.009-0.02 eV for all bases. The energy difference is larger for the $2\,{}^1A_{1g}$ state with around 0.35-0.38 eV. The same observation can be done for CCSDTQ and CC4 with similar vertical energies for all bases. Note that the $1\,{}^3B_{1g}$ state can not be described with these methods.
Then we review the vertical energies obtained with multireference methods. The smallest active space considered is four electrons in four orbitals, for CASSCF(4,4) we have small energy variations throughout bases for the $1\,{}^3B_{1g}$ and the $2\,{}^1A_{1g}$ states but a larger variation for the $1\,{}^1B_{1g}$ state with around 0.1 eV. We can observe that we have the inversion of the states compared to all methods discussed so far between the $2\,{}^1A_{1g}$ and $1\,{}^1B_{1g}$ states with $1\,{}^1B_{1g}$ higher than $2\,{}^1A_{1g}$ due to the lack of dynamical correlation in the CASSCF methods. The $1\,{}^1B_{1g}$ state values in CASSCF(4,4) is much higher than for any of the other methods discussed so far. With CASPT2(4,4) we retrieve the right ordering between the states and we see large energy differences with the CASSCF values. Indeed, we have approximatively 0.22-0.25 eV of energy difference for the triplet state for all bases and 0.32-0.36 eV for the $2\,{}^1A_{1g}$ state, the largest energy difference is for the $1\,{}^1B_{1g}$ state with 1.5-1.6 eV.
For the XMS-CASPT2(4,4) only the $2\,{}^1A_{1g}$ state is described with values similar than for the CAPST2(4,4) method. For the NEVPT2(4,4) methods (SC-NEVPT2 and PC-NEVPT2) the vertical energies are similar for the $1\,{}^1B_{1g}$ and the $2\,{}^1A_{1g}$ states with approximatively 0.002-0.003 eV and 0.02-0.03 eV of energy difference for all bases, respectively. The energy difference for the $1\,{}^1B_{1g}$ state is slightly larger with 0.05 eV for all bases. Note that for this state the vertical energy varies of 0.23 eV from the 6-31+G(d) basis to the AVDZ one. Then we use a larger active space with twelve electrons in twelve orbitals, the CASSF(12,12) values are close to the the CASSCF(4,4) for the triplet state with 0.01-0.02 eV of energy differences. For the $2\,{}^1A_{1g}$ state we have an energy difference of about 0.2 eV between the CASSCF(4,4) and the CASSCF(12,12) values. We can notice that increasing the size of the active space gives the right ordering for the states and we have an energy difference of around 0.7 eV for the $1\,{}^1B_{1g}$ state between CASSCF(4,4) and the CASSCF(12,12) values. The CASPT2(12,12) method decreases the energy of all states compared to the CASSCF(12,12) method, again the decrease is not the same for all states. We have a diminution from CASSCF to CASPT2 of about 0.17 to 0.2 eV for the $1\,{}^3B_{1g}$ and the $2\,{}^1A_{1g}$ states and for the different bases. Again, the energy difference for the $2\,{}^1A_{1g}$ state is larger with 0.5-0.7 eV depending on the basis. In a similar way than with XMS-CASPT2(4,4), the XMS-CASPT(12,12) only describes the $2\,{}^1A_{1g}$ state and the vertical energies for this state are close to the CASPT(12,12) values. For the NEVPT2(12,12) schemes we see that for the triplet and the $2\,{}^1A_{1g}$ states the energies are similar with an energy difference between the SC-NEVPT2(12,12) and the PC-NEVPT2(12,12) values of about 0.03-0.04 eV and 0.02-0.03 eV respectively.
For the XMS-CASPT2(4,4) only the $2\,{}^1A_{1g}$ state is described with values similar than for the CAPST2(4,4) method. For the NEVPT2(4,4) methods (SC-NEVPT2 and PC-NEVPT2) the vertical energies are similar for the $1\,{}^1B_{1g}$ and the $2\,{}^1A_{1g}$ states with approximatively 0.002-0.003 eV and 0.02-0.03 eV of energy difference for all bases, respectively. The energy difference for the $1\,{}^1B_{1g}$ state is slightly larger with \SI{0.05}{\eV} for all bases. Note that for this state the vertical energy varies of 0.23 eV from the 6-31+G(d) basis to the AVDZ one. Then we use a larger active space with twelve electrons in twelve orbitals, the CASSF(12,12) values are close to the CASSCF(4,4) value for the triplet state with 0.01-0.02 eV of energy differences. For the $2\,{}^1A_{1g}$ state we have an energy difference of about 0.2 eV between the CASSCF(4,4) and the CASSCF(12,12) values. We can notice that increasing the size of the active space gives the right ordering for the states and we have an energy difference of around 0.7 eV for the $1\,{}^1B_{1g}$ state between CASSCF(4,4) and the CASSCF(12,12) values. The CASPT2(12,12) method decreases the energy of all states compared to the CASSCF(12,12) method, again the decrease is not the same for all states. We have a diminution from CASSCF to CASPT2 of about 0.17 to 0.2 eV for the $1\,{}^3B_{1g}$ and the $2\,{}^1A_{1g}$ states and for the different bases. Again, the energy difference for the $2\,{}^1A_{1g}$ state is larger with 0.5-0.7 eV depending on the basis. In a similar way than with XMS-CASPT2(4,4), the XMS-CASPT(12,12) only describes the $2\,{}^1A_{1g}$ state and the vertical energies for this state are close to the CASPT(12,12) values. For the NEVPT2(12,12) schemes we see that for the triplet and the $2\,{}^1A_{1g}$ states the energies are similar with an energy difference between the SC-NEVPT2(12,12) and the PC-NEVPT2(12,12) values of about 0.03-0.04 eV and 0.02-0.03 eV respectively.
%%% TABLE II %%%
\begin{squeezetable}
@ -488,13 +487,12 @@ SF-ADC(2)-x & 6-31+G(d) & $1.557$ & $3.232$ & $3.728$ \\
& AVTZ & $1.539$ & $3.031$ & $3.703$ \\[0.1cm]
SF-ADC(3) & 6-31+G(d) & $1.435$ & $3.352$ & $4.242$ \\
& AVDZ & $1.422$ & $3.180$ & $4.208$ \\
& AVTZ & $1.419$ & $3.162$ & $4.224$ \\[0.1cm]
& AVTZ & $1.419$ & $3.162$ & $4.224$ \\
% SF-EOM-CCSD & 6-31+G(d) & $1.663$ & $3.515$ & $4.275$ \\
% & AVDZ & $1.611$ & $3.315$ & $3.856$ \\
% & AVTZ & $1.609$ & $3.293$ & $4.245$ \\[0.1cm]
%SF-EOM-CC(2,3) & 6-31+G(d) & $1.490$ & $3.333$ & $4.061$ \\
%& AVDZ & $1.464$ & $3.156$ & $4.027$ \\
\end{tabular}
\end{ruledtabular}
\end{table}
@ -597,7 +595,7 @@ Figure \ref{fig:D2h} shows the vertical energies of the studied excited states d
%%% %%% %%% %%%
\subsubsection{D4h geometry}
\subsubsection{$D_{4h}$ geometry}
\label{sec:D4h}
Table \ref{tab:sf_D4h} shows vertical energies obtained with spin-flip methods and Table \ref{tab:D4h} vertical energies obtained with standard methods. As for the previous geometry we start with the SF-TD-DFT results with hybrid functionals. We can first notice that for the B3LYP and the PBE0 functionals we have a wrong ordering of the first triplet state $1\,{}^3A_{2g}$ and the ground state $\text{X}\,{}^1B_{1g}$. We retrieve the good ordering with the BH\&HLYP functional, so adding exact exchange to the functional allows us to have the right ordering between these two states. For the B3LYP and the PBE0 functionals we have that the energy differences for each states and for all bases are small with 0.004-0.007 eV for the triplet state $1\,{}^3A_{2g}$. We have 0.015-0.021 eV of energy difference for the $2\,{}^1A_{1g}$ state through all bases, we can notice that this state is around 0.13 eV (considering all bases) higher with the PBE0 functional. We can make the same observation for the $1\,{}^1B_{2g}$ state for the B3LYP and the PBE0 functionals, indeed we have small energy differences for all bases and the state is around 0.14-0.15 eV for the PBE0 functional. For the BH\&HLYP functional the $2\,{}^1A_{1g}$ and $1\,{}^1B_{2g}$ states are higher in energy than for the two other hybrid functionals with about 0.65-0.69 eV higher for the $2\,{}^1A_{1g}$ state and 0.75-0.77 eV for the $1\,{}^1B_{2g}$ state compared to the PBE0 functional. Then, we have the RSH functionals CAM-B3LYP, $\omega$B97X-V and LC-$\omega$PBE08. For these functionals the vertical energies are similar for the $2\,{}^1A_{1g}$ and $1\,{}^1B_{2g}$ states with a maximum energy difference of 0.01-0.02 eV for the $2\,{}^1A_{1g}$ state and 0.005-0.009 eV for the $1\,{}^1B_{2g}$ state considering all bases. The maximum energy difference for the triplet state is larger with 0.047-0.057 eV for all bases. Note that the vertical energies obtained with the RSH functionals are close to the PBE0 ones except that we have the right ordering between the triplet state $1\,{}^3A_{2g}$ and the ground state $\text{X}\,{}^1B_{1g}$. The M06-2X results are closer to the BH\&HLYP ones but the M06-2X vertical energies are always higher than the BH\&HLYP ones. We can notice that the M06-2X energies for the $2\,{}^1A_{1g}$ state are close to the BH\&HLYP energies for the $1\,{}^1B_{2g}$ state. For these two states, when we compared the results obtained with the M06-2X and the BH\&HLYP functionals, we have an energy difference of 0.16-0.17 eV for the $2\,{}^1A_{1g}$ state and 0.17-0.18 eV for the $1\,{}^1B_{2g}$ state considering all bases. For the triplet state $1\,{}^3A_{2g}$ the energy differences are smaller with 0.03-0.04 eV for all bases. The M11 vertical energies are very close to the M06-2X ones for the triplet state with a maximum energy difference of 0.003 eV considering all bases, and are closer to the BH\&HLYP results for the two other states with 0.06-0.07 eV and 0.07-0.08 eV of energy difference for the $2\,{}^1A_{1g}$ and $1\,{}^1B_{2g}$ states, respectively. Then we discuss the various ADC scheme results, note that we were not able to obtain the vertical energies with the AVQZ basis due to computational resources. For the ADC(2)-s scheme we can see that the energy difference for the triplet state are smaller than for the two other states, indeed we have an energy difference of 0.09 eV for the triplet state whereas we have 0.15 eV and 0.25 eV for the $2\,{}^1A_{1g}$ and $1\,{}^1B_{2g}$ states, again when considering all bases. The energy difference for each state and through the bases are similar for the two other ADC schemes. We can notice a large variation of the vertical energies for the $2\,{}^1A_{1g}$ state between ADC(2)-s and ADC(2)-x with around 0.52-0.58 eV through all bases. The ADC(3) vertical energies are very similar to the ADC(2) ones for the $1\,{}^1B_{2g}$ state with an energy difference of 0.01-0.02 eV for all bases, whereas we have an energy difference of 0.04-0.11 eV and 0.17-0.22 eV for the $2\,{}^1A_{1g}$ and $1\,{}^1B_{2g}$ states, respectively.
@ -760,7 +758,7 @@ Figure \ref{fig:D4h} shows the vertical energies of the studied excited states d
\end{figure*}
%%% %%% %%% %%%
\subsubsection{TBE}
\subsubsection{Theoretical Best Estimates}
\label{sec:TBE}
Table \ref{tab:TBE} shows the energy differences, for the autoisomerization barrier (AB) and the different states considered, between the various methods treated and the Theoretical Best Estimate (TBE) at the AVTZ level for the AB and the states. The percentage \% T1 shown in parentheses for the excited states of the $D_{2h}$ geometry is a metric that gives the percentage of single excitation calculated at the CC3/AVTZ level and it allows us to characterize the transition. First, we look at the AB energy difference. SF-TD-DFT shows large variations of the energy with errors of 1.42-10.81 \kcalmol compared to the TBE value. SF-ADC schemes provide smaller errors with 0.30-1.44 \kcalmol where we have that the SF-ADC(2)-x gives a worse error than the SF-ADC(2)-s method. CC methods also give small energy differences with 0.11-1.05 \kcalmol and where the CC4 provides an energy very close to the TBE one.
@ -770,7 +768,7 @@ For the $1\,{}^1B_{1g} $ state of the $(D_{2h})$ structure we see that all the x
Then, for the $2\,{}^1A_{g} $ state we obtain closer results of the SF-TD-DFT methods to the TBE than in the case of the $1\,{}^1B_{1g} $ state. Indeed, we have an energy difference of about 0.01-0.34 eV for the $2\,{}^1A_{g} $ state whereas we have 0.35-0.93 eV for the $1\,{}^1B_{1g} $ state. The ADC schemes give the same error to the TBE value than for the other singlet state with 0.02 eV for the ADC(2) scheme and 0.07 eV for the ADC(3) one. The ADC(2)-x scheme provides a larger error with 0.45 eV of energy difference. Here, the CC methods manifest more variations with 0.63 eV for the CC3 value and 0.28 eV for the CCSDT compared to the TBE values. The CC4 method provides a small error with less than 0.01 eV of energy difference. The multiconfigurational methods globally give smaller error than for the other singlet state with, for the two active spaces, 0.03-0.12 eV compared to the TBE value. We can notice that CC3 and CCSDT provide larger energy errors for the $2\,{}^1A_{g} $ state than for the $1\,{}^1B_{1g} $ state, this is due to the strong multiconfigurational character of the $2\,{}^1A_{g} $ state whereas the $1\,{}^1B_{1g} $ state has a very weak multiconfigurational character. It is interesting to see that SF-TD-DFT and SF-ADC methods give better results compared to the TBE than CC3 and even CCSDT meaning that spin-flip methods are able to describe double excited states. Note that multireference methods obviously give better results too for the $2\,{}^1A_{g} $ state.
Finally we look at the vertical energy errors for the $D_{4h}$ structure. First, we consider the $1\,{}^3A_{2g} $ state, the SF-TD-DFT methods give errors of about 0.07-1.6 eV where the largest energy differences are provided by the hybrid functionals. The ADC schemes give similar error with around 0.06-1.1 eV of energy difference. For the CC methods we have an energy error of 0.06 eV for CCSD and less than 0.01 eV for CCSDT. Then for the multireference methods with the four by four active space we have for CASSCF(4,4) 0.29 eV of error and 0.02 eV for CASPT2(4,4), again CASPT2 demonstrates its improvement compared to CASSCF. The other methods provide energy differences of about 0.12-0.13 eV. A larger active space shows again an improvement with 0.23 eV of error for CASSCF(12,12) and around 0.01-0.04 eV for the other multireference methods. CIPSI provides similar error with 0.02 eV. Then, we look at the $2\,{}^1A_{1g}$ state where the SF-TD-DFT shows large variations of error depending on the functionals, the energy error is about 0.10-1.03 eV. The ADC schemes give better errors with around 0.07-0.41 eV and where again the ADC(2)-x does not improve the ADC(2)-s energy but also gives worse results. For the CC methods we have energy errors of about 0.10-0.16 eV and CC4 provides really close energy to the TBE one with 0.01 eV of error. For the multireference methods we globally have an improvement of the energies from the four by four to the twelve by twelve active space with errors of 0.01-0.73 eV and 0.02-0.44 eV respectively with the largest errors coming from the CASSCF method. Lastly, we look at the $1\,{}^1B_{2g}$ state where we have globally larger errors. The SF-TD-DFT exhibits errors of 0.43-1.50 eV whereas ADC schemes give errors of 0.18-0.30 eV. CC3 and CCSDT provide energy differences of 0.50-0.69 eV and the CC4 shows again close energy to the CCSDTQ TBE energy with 0.01 eV of error. The multireference methods give energy differences of 0.38-1.39 eV and 0.11-0.60 eV for the four by four and twelve by twelve active spaces respectively. We can notice that we have larger variations for the vertical energies of the square CBD than for the rectangular one. This can be explained by the fact that because of the degeneracy in the $D_{4h}$ structure it leads to strong multiconfigurational character states where single reference methods are unreliable. We can also see that for the CC methods we have a better description of the $2\,{}^1A_{1g}$ state than the $1\,{}^1B_{2g}$ state, this can be related, as previously said in Subsubsection \ref{sec:D4h}, to the fact that $1\,{}^1B_{2g}$ corresponds to a double excitation from the reference state. To obtain an improved description of the $1\,{}^1B_{2g}$ state we have to include quadruples. At the end of Table \ref{tab:TBE} we show some literature results obtain from Ref ~\cite{lefrancois_2015,manohar_2008} where the cc-pVTZ basis is used. The SF-ADC(2)-s, SF-ADC(2)-x and SF-ADC(3)results are presented and are consistent with our results with the exact same schemes but with the AVTZ basis.
Finally we look at the vertical energy errors for the $D_{4h}$ structure. First, we consider the $1\,{}^3A_{2g} $ state, the SF-TD-DFT methods give errors of about 0.07-1.6 eV where the largest energy differences are provided by the hybrid functionals. The ADC schemes give similar error with around 0.06-1.1 eV of energy difference. For the CC methods we have an energy error of 0.06 eV for CCSD and less than 0.01 eV for CCSDT. Then for the multireference methods with the four by four active space we have for CASSCF(4,4) 0.29 eV of error and 0.02 eV for CASPT2(4,4), again CASPT2 demonstrates its improvement compared to CASSCF. The other methods provide energy differences of about 0.12-0.13 eV. A larger active space shows again an improvement with 0.23 eV of error for CASSCF(12,12) and around 0.01-0.04 eV for the other multireference methods. CIPSI provides similar error with 0.02 eV. Then, we look at the $2\,{}^1A_{1g}$ state where the SF-TD-DFT shows large variations of error depending on the functionals, the energy error is about 0.10-1.03 eV. The ADC schemes give better errors with around 0.07-0.41 eV and where again the ADC(2)-x does not improve the ADC(2)-s energy but also gives worse results. For the CC methods we have energy errors of about 0.10-0.16 eV and CC4 provides really close energy to the TBE one with 0.01 eV of error. For the multireference methods we globally have an improvement of the energies from the four by four to the twelve by twelve active space with errors of 0.01-0.73 eV and 0.02-0.44 eV respectively with the largest errors coming from the CASSCF method. Lastly, we look at the $1\,{}^1B_{2g}$ state where we have globally larger errors. The SF-TD-DFT exhibits errors of 0.43-1.50 eV whereas ADC schemes give errors of 0.18-0.30 eV. CC3 and CCSDT provide energy differences of 0.50-0.69 eV and the CC4 shows again close energy to the CCSDTQ TBE energy with 0.01 eV of error. The multireference methods give energy differences of 0.38-1.39 eV and 0.11-0.60 eV for the four by four and twelve by twelve active spaces respectively. We can notice that we have larger variations for the vertical energies of the square CBD than for the rectangular one. This can be explained by the fact that because of the degeneracy in the $D_{4h}$ structure it leads to strong multiconfigurational character states where single reference methods are unreliable. We can also see that for the CC methods we have a better description of the $2\,{}^1A_{1g}$ state than the $1\,{}^1B_{2g}$ state, this can be related, as previously said in Subsubsection \ref{sec:D4h}, to the fact that $1\,{}^1B_{2g}$ corresponds to a double excitation from the reference state. To obtain an improved description of the $1\,{}^1B_{2g}$ state we have to include quadruples. At the end of Table \ref{tab:TBE} we show some literature results obtain from Ref.~\onlinecite{lefrancois_2015,manohar_2008} where the cc-pVTZ basis is used. The SF-ADC(2)-s, SF-ADC(2)-x and SF-ADC(3)results are presented and are consistent with our results with the exact same schemes but with the AVTZ basis.
%Here again we can make the same comment for the $2\,{}^1A_{1g}$ and $1\,{}^1B_{2g}$ states of the square CBD than the $1\,{}^1B_{1g}$ and $2\,{}^1A_{g}$ states of the rectangular CBD. The first state ($2\,{}^1A_{1g}$) has a strong multiconfigurational character
@ -840,9 +838,9 @@ Finally we look at the vertical energy errors for the $D_{4h}$ structure. First,
\begin{ruledtabular}
\begin{tabular}{lrrrrrrr}
%\begin{tabular}{*{1}{*{8}{l}}}
&\mc{3}{r}{$D_{2h}$ excitation energies (eV)} & \mc{3}{r}{$D_{4h}$ excitation energies (eV)} \\
&&\mc{3}{c}{$D_{2h}$ excitation energies (eV)} & \mc{3}{c}{$D_{4h}$ excitation energies (eV)} \\
\cline{3-5} \cline{6-8}
Method & AB & $1\,{}^3B_{1g} $~(98.7 \%) & $1\,{}^1B_{1g} $ (95.0 \%)& $2\,{}^1A_{g} $(0.84 \%) & $1\,{}^3A_{2g} $ & $2\,{}^1A_{1g} $ & $1\,{}^1B_{2g} $ \\
Method & AB & $1\,{}^3B_{1g} $~(99\%) & $1\,{}^1B_{1g} $ (95\%)& $2\,{}^1A_{g} $(1\%) & $1\,{}^3A_{2g} $ & $2\,{}^1A_{1g} $ & $1\,{}^1B_{2g} $ \\
\hline
SF-TD-B3LYP & $10.41$ & $0.241$ & $-0.926$ & $-0.161$ & $-0.164$ & $-1.028$ & $-1.501$ \\
SF-TD-PBE0 & $8.95$ & $0.220$ & $-0.829$ & $-0.068$ & $-0.163$ & $-0.903$ & $-1.357$ \\
@ -873,7 +871,7 @@ SC-NEVPT2(12,12) & $-0.65$ & $0.039$ & $0.063$ & $-0.063$ & $0.021$ & $0.046$ &
PC-NEVPT2(12,12) & & $0.000$ & $-0.062$ & $-0.093$ & $-0.013$ & $-0.024$ & $-0.278$ \\[0.1cm]
%CIPSI & & $-0.001\pm 0.030$ & $0.017\pm 0.035$ & $-0.120\pm 0.090$ & $0.025\pm 0.029$ & $0.130\pm 0.050$ & \\
\bf{TBE} & $\left[\bf{8.93}\right]$\fnm[1] & $\left[\bf{1.462}\right]$\fnm[2] & $\left[\bf{3.125}\right]$\fnm[3] & $\left[\bf{4.149}\right]$\fnm[3] & $\left[\bf{0.144}\right]$\fnm[4] & $\left[\bf{1.500}\right]$\fnm[3] & $\left[\bf{2.034}\right]$\fnm[3] \\[0.1cm]
lit. & $8.53$\fnm[5] & $1.573$\fnm[5] & $3.208$\fnm[5] & $4.247$\fnm[5] & $0.266$\fnm[5] & $1.664$\fnm[5] & $1.910$\fnm[5] \\
Literature & $8.53$\fnm[5] & $1.573$\fnm[5] & $3.208$\fnm[5] & $4.247$\fnm[5] & $0.266$\fnm[5] & $1.664$\fnm[5] & $1.910$\fnm[5] \\
& $10.35$\fnm[6] & $1.576$\fnm[6] & $3.141$\fnm[6] & $3.796$\fnm[6] & $0.217$\fnm[6] & $1.123$\fnm[6] & $1.799$\fnm[6]\\
& $9.58$ \fnm[7]& $1.456$\fnm[7] & $3.285$\fnm[7] & $4.334$\fnm[7] & $0.083$\fnm[7] & $1.621$\fnm[7] & $1.930$\fnm[7] \\
& $7.48$\fnm[8]& $1.654$\fnm[8] & $3.416$\fnm[8] & $4.360$\fnm[8] & $0.369$\fnm[8] & $1.824$\fnm[8] & $2.143$\fnm[8] \\
@ -884,10 +882,10 @@ lit. & $8.53$\fnm[5] & $1.573$\fnm[5] & $3.208$\fnm[5] & $4.247$\fnm[5] & $0.266
\fnt[2]{Value obtained using the NEVPT2(12,12) one.}
\fnt[3]{Value obtained using CCSDTQ/AVDZ corrected by the difference between CC4/AVTZ and CC4/AVDZ.}
\fnt[4]{Value obtained using CCSDTQ/AVDZ corrected by the difference between CCSDT/AVTZ and CCSDT/AVDZ.}
\fnt[5]{Value obtained from Ref ~\cite{lefrancois_2015} at the SF-ADC(2)-s/cc-pVTZ level with the geometry obtained at the CCSD(T)/cc-pVTZ level.}
\fnt[6]{Value obtained from Ref ~\cite{lefrancois_2015} at the SF-ADC(2)-x/cc-pVTZ level with the geometry obtained at the CCSD(T)/cc-pVTZ level.}
\fnt[7]{Value obtained from Ref ~\cite{lefrancois_2015} at the SF-ADC(3)/cc-pVTZ level with the geometry obtained at the CCSD(T)/cc-pVTZ level.}
\fnt[8]{Value obtained from Ref ~\cite{manohar_2008} at the EOM-SF-CCSD/cc-pVTZ level with the geometry obtained at the CCSD(T)/cc-pVTZ level.}
\fnt[5]{Value obtained from Ref.~\onlinecite{lefrancois_2015} at the SF-ADC(2)-s/cc-pVTZ level with the geometry obtained at the CCSD(T)/cc-pVTZ level.}
\fnt[6]{Value obtained from Ref.~\onlinecite{lefrancois_2015} at the SF-ADC(2)-x/cc-pVTZ level with the geometry obtained at the CCSD(T)/cc-pVTZ level.}
\fnt[7]{Value obtained from Ref.~\onlinecite{lefrancois_2015} at the SF-ADC(3)/cc-pVTZ level with the geometry obtained at the CCSD(T)/cc-pVTZ level.}
\fnt[8]{Value obtained from Ref.~\onlinecite{manohar_2008} at the EOM-SF-CCSD/cc-pVTZ level with the geometry obtained at the CCSD(T)/cc-pVTZ level.}
\end{table*}
\end{squeezetable}