diff --git a/Manuscript/CBD.tex b/Manuscript/CBD.tex index 7a72f5e..69a15eb 100644 --- a/Manuscript/CBD.tex +++ b/Manuscript/CBD.tex @@ -198,15 +198,17 @@ \sisetup{range-units=single} %states -\newcommand{\Ag}{$2{}^1A_g$} -\newcommand{\Aoneg}{$2{}^1A_{1g}$} -\newcommand{\Btwog}{$1{}^1B_{2g}$} -\newcommand{\Atwog}{$1{}^3A_{2g}$} +%D2h states +\newcommand{\oneAg}{$1{}^1A_g$} \newcommand{\tBoneg}{$1{}^3B_{1g}$} \newcommand{\sBoneg}{$1{}^1B_{1g}$} +\newcommand{\twoAg}{$2{}^1A_g$} - - +%D4h states +%\newcommand{\oneBoneg}{$1{}^1B_{1g}$} same label as the D2h state +\newcommand{\Atwog}{$1{}^3A_{2g}$} +\newcommand{\Aoneg}{$2{}^1A_{1g}$} +\newcommand{\Btwog}{$1{}^1B_{2g}$} \begin{document} @@ -254,7 +256,7 @@ This degeneracy is lifted by the so-called Jahn-Teller effect, \ie, by a descent In such as case, H\"uckel molecular orbital theory (correctly) predicts a closed-shell singlet ground state at the {\Dtwo} rectangular geometry. \titou{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 {\Dtwo} ground-state structure, the \titou{{\Ag}} state has a weak multi-configurational character with well-separated frontier orbitals that can be described by single-reference methods. +At the {\Dtwo} ground-state structure, the \titou{{\oneAg}} state has a weak multi-configurational character with well-separated frontier orbitals that can be described by single-reference methods. However, at the {\Dfour} square geometry, the singlet state {\sBoneg} 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 (\Dfour) 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 (\Dfour) structure and the singlet ground state of the rectangular (\Dtwo) geometry. @@ -262,9 +264,8 @@ The energy of this barrier was predicted, experimentally, in the range of 1.6-10 All the specificities of CBD make it a real playground for excited-states methods. The lowest-energy excited states of CBD in both geometries are represented in Fig.~\ref{fig:CBD}. -Are represented {\Ag} and {\tBoneg} states for the rectangular geometry and {\sBoneg} and {\Atwog} for the square one. -Due to energy scaling doubly excited state {\sBoneg} and {\Aoneg} for the {\Dtwo} and {\Dfour} structures, respectively, are not drawn. -Doubly excited states are known to be challenging to represent for adiabatic time-dependent density functional theory \cite{casida_1995} (TD-DFT) and even for state-of-the-art methods like the approximate third-order coupled-cluster (CC3) \cite{christiansen_1995,koch_1997} or equation-of-motion coupled-cluster with singles, doubles and triples (EOM-CCSDT). \cite{kucharski_1991,kallay_2004,hirata_2000,hirata_2004} +Are represented {\oneAg} and {\tBoneg} states for the rectangular geometry and {\sBoneg} and {\Atwog} for the square one. +Due to energy scaling, the {\sBoneg}, {\twoAg} and {\Aoneg}, {\Btwog}, states for the {\Dtwo} and {\Dfour} structures, respectively, are not drawn. The {\twoAg} and {\Aoneg} states are doubly excited states and these doubly excited states are known to be challenging to represent for adiabatic time-dependent density functional theory \cite{casida_1995} (TD-DFT) and even for state-of-the-art methods like the approximate third-order coupled-cluster (CC3) \cite{christiansen_1995,koch_1997} or equation-of-motion coupled-cluster with singles, doubles and triples (EOM-CCSDT). \cite{kucharski_1991,kallay_2004,hirata_2000,hirata_2004} In order to tackle the problems of multi-configurational character and double excitations several ways are explored. The most evident way that one can think about to describe multiconfigurational and double excitations are multiconfigurational methods. Among these methods, one can find complete active space self-consistent field (CASSCF) \cite{roos_1996}, the second perturbation-corrected variant (CASPT2) \cite{andersson_1990} and the second-order $n$-electron valence state perturbation theory (NEVPT2). \cite{angeli_2001,angeli_2001a,angeli_2002} @@ -283,12 +284,12 @@ So one can access the singlet ground state and the singlet doubly-excited state Obviously spin-flip methods have their own flaws, especially the spin contamination \cite{casanova_2020} (\ie, 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 $1{}^1A_{g}$, {\tBoneg}, {\sBoneg}, {\Ag} and {\sBoneg}, {\Atwog}, {\Aoneg},{\Btwog} excited states for the {\Dtwo} and {\Dfour} geometries, respectively. +In the present work we investigate {\oneAg}, {\tBoneg}, {\sBoneg}, {\twoAg} and {\sBoneg}, {\Atwog}, {\Aoneg},{\Btwog} excited states for the {\Dtwo} and {\Dfour} 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. \begin{figure} - \includegraphics[width=0.6\linewidth]{figure2.png} + \includegraphics[width=0.6\linewidth]{figure1.pdf} \caption{Pictorial representation of the ground and excited states of CBD and its properties under investigation. The singlet ground-state (S) and triplet (T) properties are represented in black and red, respectively.} \label{fig:CBD} @@ -372,13 +373,13 @@ Two different sets of geometries obtained with different level of theory are con \begin{tabular}{lllrrr} State & Method & \ce{C=C} & \ce{C-C} & \ce{C-H} & \ce{H-C=C}\fnm[1] \\ \hline - {\Dtwo} $({}^1 A_{g})$& + {\Dtwo} ({\oneAg}) & CASPT2(12,12)/aug-cc-pVTZ & 1.355 & 1.566 & 1.077 & 134.99 \\ &CC3/aug-cc-pVTZ & 1.344 & 1.565 & 1.076 & 135.08 \\ &CCSD(T)/cc-pVTZ & 1.343 & 1.566 & 1.074 & 135.09 \fnm[2]\\ - {\Dfour} $({}^1 B_{1g})$& + {\Dfour} ({\sBoneg}) & CASPT2(12,12)/aug-cc-pVTZ & 1.449 & 1.449 & 1.076 & 135.00 \\ - {\Dfour} $({}^3 A_{2g})$& + {\Dfour} ({\Atwog}) & CASPT2(12,12)/aug-cc-pVTZ & 1.445 & 1.445 & 1.076 & 135.00 \\ &RO-CCSD(T)/aug-cc-pVTZ & 1.439 & 1.439 & 1.075 & 135.00 \end{tabular} @@ -463,24 +464,24 @@ Figure \ref{fig:AB} shows the autoisomerization barrier (AB) energies for the CB %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{{\Dtwo} 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 {\tBoneg} state with \SI{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 {\sBoneg} and the {\tBoneg} states. Put another way, increasing the amount of exact exchange in the functional reduces the energy difference between the {\tBoneg} 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 {\tBoneg} state from PBE0 to BH\&HLYP is around \SI{0.1}{\eV} whereas for the {\sBoneg} and the {\Ag} states this energy variation is about \SIrange{0.4}{0.5}{\eV} and \SI{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 {\sBoneg} and the {\tBoneg} 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 \SIrange{0.03}{0.08}{\eV}. We can notice that the upper bound of \SI{0.08}{\eV} in the energy differences is due to the {\tBoneg} state. The M11 vertical energies are close to the BH\&HLYP ones for the triplet and the first singlet excited states with \SIrange{0.01}{0.02}{\eV} and \SIrange{0.08}{0.10}{\eV} of energy difference, respectively. For the {\Ag} state the M11 energies are closer to the $\omega$B97X-V ones with \SIrange{0.05}{0.09}{\eV} of energy difference. +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 {\tBoneg} state with \SI{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 {\sBoneg} and the {\tBoneg} states. Put another way, increasing the amount of exact exchange in the functional reduces the energy difference between the {\tBoneg} and the {\oneAg} 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 {\tBoneg} state from PBE0 to BH\&HLYP is around \SI{0.1}{\eV} whereas for the {\sBoneg} and the {\twoAg} states this energy variation is about \SIrange{0.4}{0.5}{\eV} and \SI{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 {\sBoneg} and the {\tBoneg} 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 \SIrange{0.03}{0.08}{\eV}. We can notice that the upper bound of \SI{0.08}{\eV} in the energy differences is due to the {\tBoneg} state. The M11 vertical energies are close to the BH\&HLYP ones for the triplet and the first singlet excited states with \SIrange{0.01}{0.02}{\eV} and \SIrange{0.08}{0.10}{\eV} of energy difference, respectively. For the {\twoAg} state the M11 energies are closer to the $\omega$B97X-V ones with \SIrange{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 differences of about 0.03 eV for the {\tBoneg} state and around 0.06 eV for {\Ag} state throughout all bases. However for the {\sBoneg} state we have an energy difference of about 0.2 eV. For the ADC(2)-x and ADC(3) schemes the calculations with the aug-cc-pVQZ basis are not feasible with our resources. With ADC(2)-x we have similar vertical energies for the triplet and the {\sBoneg} states but for the {\Ag} state the energy difference between the ADC(2) and ADC(2)-x schemes is about \SIrange{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 \SIrange{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 {\tBoneg} and the {\sBoneg} states the CCSDT and the CC3 values are close with an energy difference of \SIrange{0.01}{0.02}{\eV} for all bases. The energy difference is larger for the {\Aoneg} state with around \SIrange{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 {\tBoneg} 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 {\tBoneg} state and around 0.06 eV for {\twoAg} state throughout all bases. However for the {\sBoneg} state we have an energy difference of about 0.2 eV. For the ADC(2)-x and ADC(3) schemes the calculations with the aug-cc-pVQZ basis are not feasible with our resources. With ADC(2)-x we have similar vertical energies for the triplet and the {\sBoneg} states but for the {\twoAg} state the energy difference between the ADC(2) and ADC(2)-x schemes is about \SIrange{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 \SIrange{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 {\tBoneg} and the {\sBoneg} states the CCSDT and the CC3 values are close with an energy difference of \SIrange{0.01}{0.02}{\eV} for all bases. The energy difference is larger for the {\twoAg} state with around \SIrange{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 {\tBoneg} 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 {\tBoneg} and the {\Ag} states but a larger variation for the {\sBoneg} state with around \SI{0.1}{\eV}. We can observe that we have the inversion of the states compared to all methods discussed so far between the {\Ag} and {\sBoneg} states with {\sBoneg} higher than {\Ag} due to the lack of dynamical correlation in the CASSCF methods. The {\sBoneg} 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 \SIrange{0.22}{0.25}{\eV} of energy difference for the triplet state for all bases and \SIrange{0.32}{0.36}{\eV} for the {\Ag} state, the largest energy difference is for the {\sBoneg} state with \SIrange{1.5}{1.6}{\eV}. -For the XMS-CASPT2(4,4) only the {\Ag} 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 {\sBoneg} and the {\Ag} states with approximatively \SIrange{0.002}{0.003}{\eV} and \SIrange{0.02}{0.03}{\eV} of energy difference for all bases, respectively. The energy difference for the {\sBoneg} state is slightly larger with \SI{0.05}{\eV} for all bases. Note that for this state the vertical energy varies of \SI{0.23}{eV} from the 6-31+G(d) basis to the aug-cc-pVDZ 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 {\Ag} state we have an energy difference of about \SI{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 \SI{0.7}{\eV} for the {\sBoneg} 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 \SIrange{0.17}{0.2}{\eV} for the {\tBoneg} and the {\Ag} states and for the different bases. Again, the energy difference for the {\Ag} state is larger with \SIrange{0.5}{0.7}{\eV} depending on the basis. In a similar way than with XMS-CASPT2(4,4), XMS-CASPT(12,12) only describes the {\Ag} 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 {\Aoneg} states the energies are similar with an energy difference between the SC-NEVPT2(12,12) and the PC-NEVPT2(12,12) values of about \SIrange{0.03}{0.04}{\eV} and \SIrange{0.02}{0.03}{\eV} respectively. +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 {\tBoneg} and the {\twoAg} states but a larger variation for the {\sBoneg} state with around \SI{0.1}{\eV}. We can observe that we have the inversion of the states compared to all methods discussed so far between the {\twoAg} and {\sBoneg} states with {\sBoneg} higher than {\twoAg} due to the lack of dynamical correlation in the CASSCF methods. The {\sBoneg} 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 \SIrange{0.22}{0.25}{\eV} of energy difference for the triplet state for all bases and \SIrange{0.32}{0.36}{\eV} for the {\twoAg} state, the largest energy difference is for the {\sBoneg} state with \SIrange{1.5}{1.6}{\eV}. +For the XMS-CASPT2(4,4) only the {\twoAg} 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 {\sBoneg} and the {\twoAg} states with approximatively \SIrange{0.002}{0.003}{\eV} and \SIrange{0.02}{0.03}{\eV} of energy difference for all bases, respectively. The energy difference for the {\sBoneg} state is slightly larger with \SI{0.05}{\eV} for all bases. Note that for this state the vertical energy varies of \SI{0.23}{eV} from the 6-31+G(d) basis to the aug-cc-pVDZ 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 {\twoAg} state we have an energy difference of about \SI{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 \SI{0.7}{\eV} for the {\sBoneg} 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 \SIrange{0.17}{0.2}{\eV} for the {\tBoneg} and the {\twoAg} states and for the different bases. Again, the energy difference for the {\twoAg} state is larger with \SIrange{0.5}{0.7}{\eV} depending on the basis. In a similar way than with XMS-CASPT2(4,4), XMS-CASPT(12,12) only describes the {\twoAg} 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 {\twoAg} states the energies are similar with an energy difference between the SC-NEVPT2(12,12) and the PC-NEVPT2(12,12) values of about \SIrange{0.03}{0.04}{\eV} and \SIrange{0.02}{0.03}{\eV} respectively. %%% TABLE II %%% \begin{squeezetable} \begin{table} \caption{ - Spin-flip TD-DFT vertical excitation energies (with respect to the singlet $\text{X}\,{}^1A_{g}$ ground state) of the {\tBoneg}, {\sBoneg}, and {\Ag} states of CBD at the {\Dtwo} rectangular equilibrium geometry of the $\text{X}\,{}^1 A_{g}$ ground state. + Spin-flip TD-DFT vertical excitation energies (with respect to the singlet {\oneAg} ground state) of the {\tBoneg}, {\sBoneg}, and {\twoAg} states of CBD at the {\Dtwo} rectangular equilibrium geometry of the {\oneAg} ground state. \label{tab:sf_tddft_D2h}} \begin{ruledtabular} \begin{tabular}{llrrr} & \mc{4}{r}{Excitation energies (eV)} \hspace{0.5cm}\\ \cline{3-5} - Method & Basis & {\tBoneg} & {\sBoneg} & {\Ag} \\ + Method & Basis & {\tBoneg} & {\sBoneg} & {\twoAg} \\ \hline % SF-TD-BLYP & 6-31+G(d) & $1.829$ & $1.926$ & $3.755$ \\ % & aug-cc-pVDZ & $1.828$ & $1.927$ & $3.586$ \\ @@ -549,13 +550,13 @@ SF-ADC(3) & 6-31+G(d) & $1.435$ & $3.352$ & $4.242$ \\ \begin{squeezetable} \begin{table} \caption{ - Vertical excitation energies (with respect to the singlet $\text{X}\,{}^1A_{g}$ ground state) of the {\tBoneg}, {\sBoneg}, and{\Ag} states of CBD at the {\Dtwo} rectangular equilibrium geometry of the $\text{X}\,{}^1 A_{g}$ ground state. + Vertical excitation energies (with respect to the singlet {\oneAg} ground state) of the {\tBoneg}, {\sBoneg}, and {\twoAg} states of CBD at the {\Dtwo} rectangular equilibrium geometry of the {\oneAg} ground state. \label{tab:D2h}} \begin{ruledtabular} \begin{tabular}{llrrr} & & \mc{3}{c}{Excitation energies (eV)} \\ \cline{3-5} - Method & Basis & {\tBoneg} & {\sBoneg} & {\Ag} \\ + Method & Basis & {\tBoneg} & {\sBoneg} & {\twoAg} \\ \hline CC3 &6-31+G(d)& $1.420$ & $3.341$ & $4.658$ \\ & aug-cc-pVDZ & $1.396$ & $3.158$ & $4.711$ \\ @@ -627,13 +628,13 @@ CIPSI &6-31+G(d)& $1.486\pm 0.005$ & $3.348\pm 0.024$ & $4.084\pm 0.012$ \\ \end{squeezetable} %%% %%% %%% %%% -Figure \ref{fig:D2h} shows the vertical energies of the studied excited states described in Tables \ref{tab:sf_tddft_D2h} and \ref{tab:D2h}. We see that the SF-TD-DFT vertical energies are rather close to the TBE one for the doubly excited state {\Ag}. Indeed, as presented later in subsection \ref{sec:TBE}, SF methods are able to describe effectively states with strong multiconfigurational character. The same observation can be done for the SF-ADC values but with much better results for the two other states. Again we can notice that the SF-ADC(2)-x and SF-ADC(3) schemes do not improve the values. For the multireference methods we see that using the smallest active space do not provide a good description of the {\sBoneg} state and in the case of CASSCF(4,4), as previously said, we even have {\sBoneg} state above the {\Ag} state. For the CC methods taking into account the quadruple excitations is necessary to have a good description of the doubly excited state {\Ag}. +Figure \ref{fig:D2h} shows the vertical energies of the studied excited states described in Tables \ref{tab:sf_tddft_D2h} and \ref{tab:D2h}. We see that the SF-TD-DFT vertical energies are rather close to the TBE one for the doubly excited state {\twoAg}. Indeed, as presented later in subsection \ref{sec:TBE}, SF methods are able to describe effectively states with strong multiconfigurational character. The same observation can be done for the SF-ADC values but with much better results for the two other states. Again we can notice that the SF-ADC(2)-x and SF-ADC(3) schemes do not improve the values. For the multireference methods we see that using the smallest active space do not provide a good description of the {\sBoneg} state and in the case of CASSCF(4,4), as previously said, we even have {\sBoneg} state above the {\twoAg} state. For the CC methods taking into account the quadruple excitations is necessary to have a good description of the doubly excited state {\twoAg}. %%% FIGURE III %%% \begin{figure*} %width=0.8\linewidth \includegraphics[scale=0.5]{D2h.pdf} - \caption{Vertical energies of the {\tBoneg}, {\sBoneg} and {\Ag} states for the {\Dtwo} geometry using the aug-cc-pVTZ basis. Purple, orange, green, blue and black lines correspond to the SF-TD-DFT, SF-ADC, multireference, CC, and TBE values, respectively.} + \caption{Vertical energies of the {\tBoneg}, {\sBoneg} and {\twoAg} states for the {\Dtwo} geometry using the aug-cc-pVTZ basis. Purple, orange, green, blue and black lines correspond to the SF-TD-DFT, SF-ADC, multireference, CC, and TBE values, respectively.} \label{fig:D2h} \end{figure*} %%% %%% %%% %%% @@ -641,17 +642,17 @@ Figure \ref{fig:D2h} shows the vertical energies of the studied excited states d \subsubsection{{\Dfour} 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 {\Atwog} 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 \SIrange{0.004}{0.007}{\eV} for the triplet state {\Atwog}. We have \SIrange{0.015}{0.021}{\eV} of energy difference for the {\Aoneg} state through all bases, we can notice that this state is around \SI{0.13}{\eV} (considering all bases) higher with the PBE0 functional. We can make the same observation for the {\Btwog} state for the B3LYP and the PBE0 functionals, indeed we have small energy differences for all bases and the state is around \SIrange{0.14}{0.15}{\eV} for the PBE0 functional. For the BH\&HLYP functional the {\Aoneg} and {\Btwog} states are higher in energy than for the two other hybrid functionals with about \SIrange{0.65}{0.69}{\eV} higher for the {\Aoneg} state and \SIrange{0.75}{0.77}{\eV} for the {\Btwog} 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 {\Aoneg} and {\Btwog} states with a maximum energy difference of \SIrange{0.01}{0.02}{\eV} for the{\Aoneg} state and \SIrange{0.005}{0.009}{\eV} for the {\Btwog} state considering all bases. The maximum energy difference for the triplet state is larger with \SIrange{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 {\Atwog} 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 {\Aoneg} state are close to the BH\&HLYP energies for the {\Btwog} 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 \SIrange{0.16}{0.17}{\eV} for the {\Aoneg} state and \SIrange{0.17}{0.18}{\eV} for the {\Btwog} state considering all bases. For the triplet state {\Atwog} the energy differences are smaller with \SIrange{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 \SI{0.003}{\eV} considering all bases, and are closer to the BH\&HLYP results for the two other states with \SIrange{0.06}{0.07}{\eV} and \SIrange{0.07}{0.08}{\eV} of energy difference for the {\Aoneg} and {\Btwog} states, respectively. Then we discuss the various ADC scheme results, note that we were not able to obtain the vertical energies with the aug-cc-pVQZ 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 \SI{0.09}{\eV} for the triplet state whereas we have \SI{0.15}{\eV} and \SI{0.25}{eV} for the {\Aoneg} and {\Btwog} 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 {\Aoneg} state between ADC(2)-s and ADC(2)-x with around \SIrange{0.52}{0.58}{\eV} through all bases. The ADC(3) vertical energies are very similar to the ADC(2) ones for the {\Btwog} state with an energy difference of \SIrange{0.01}{0.02}{\eV} for all bases, whereas we have an energy difference of \SIrange{0.04}{0.11}{\eV} and \SIrange{0.17}{0.22}{\eV} for the {\Aoneg} and {\Btwog} states, respectively. +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 {\Atwog} and the ground state {\sBoneg}. 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 \SIrange{0.004}{0.007}{\eV} for the triplet state {\Atwog}. We have \SIrange{0.015}{0.021}{\eV} of energy difference for the {\Aoneg} state through all bases, we can notice that this state is around \SI{0.13}{\eV} (considering all bases) higher with the PBE0 functional. We can make the same observation for the {\Btwog} state for the B3LYP and the PBE0 functionals, indeed we have small energy differences for all bases and the state is around \SIrange{0.14}{0.15}{\eV} for the PBE0 functional. For the BH\&HLYP functional the {\Aoneg} and {\Btwog} states are higher in energy than for the two other hybrid functionals with about \SIrange{0.65}{0.69}{\eV} higher for the {\Aoneg} state and \SIrange{0.75}{0.77}{\eV} for the {\Btwog} 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 {\Aoneg} and {\Btwog} states with a maximum energy difference of \SIrange{0.01}{0.02}{\eV} for the{\Aoneg} state and \SIrange{0.005}{0.009}{\eV} for the {\Btwog} state considering all bases. The maximum energy difference for the triplet state is larger with \SIrange{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 {\Atwog} and the ground state {\sBoneg}. 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 {\Aoneg} state are close to the BH\&HLYP energies for the {\Btwog} 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 \SIrange{0.16}{0.17}{\eV} for the {\Aoneg} state and \SIrange{0.17}{0.18}{\eV} for the {\Btwog} state considering all bases. For the triplet state {\Atwog} the energy differences are smaller with \SIrange{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 \SI{0.003}{\eV} considering all bases, and are closer to the BH\&HLYP results for the two other states with \SIrange{0.06}{0.07}{\eV} and \SIrange{0.07}{0.08}{\eV} of energy difference for the {\Aoneg} and {\Btwog} states, respectively. Then we discuss the various ADC scheme results, note that we were not able to obtain the vertical energies with the aug-cc-pVQZ 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 \SI{0.09}{\eV} for the triplet state whereas we have \SI{0.15}{\eV} and \SI{0.25}{eV} for the {\Aoneg} and {\Btwog} 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 {\Aoneg} state between ADC(2)-s and ADC(2)-x with around \SIrange{0.52}{0.58}{\eV} through all bases. The ADC(3) vertical energies are very similar to the ADC(2) ones for the {\Btwog} state with an energy difference of \SIrange{0.01}{0.02}{\eV} for all bases, whereas we have an energy difference of \SIrange{0.04}{0.11}{\eV} and \SIrange{0.17}{0.22}{\eV} for the {\Aoneg} and {\Btwog} states, respectively. Now we look at Table \ref{tab:D4h} for vertical energies obtained with standard methods, we have first the various CC methods, for all the CC methods we were only able to reach the aug-cc-pVTZ basis. %For the CCSD method we cannot obtain the vertical energy for the $1\,{}^1B_{2g}$ state, we can notice a small energy variation through the bases for the triplet state with around 0.06 eV. The energy variation is larger for the $2\,{}^1A_{1g}$ state with 0.20 eV. -Again, we have computed the vertical energies at the CCSDT and CCSDTQ level as well as their approximation the CC3 and CC4 methods, respectively. Note that for CC we started from a restricted Hartree-Fock (RHF) reference and in that case the ground state is the {\Aoneg}state, then the $X\,{}^1B_{1g}$ state is the single deexcitation and the{\Btwog} state is the double excitation from our ground state. For the CC3 method we do not have the vertical energies for the triplet state {\Atwog}. Considering all bases for the {\Aoneg}and {\Btwog} states we have an energy difference of about \SI{0.15}{\eV} and \SI{0.12}{\eV}, respectively. The CCSDT energies are close to the CC3 ones for the {\Aoneg} state with an energy difference of around \SIrange{0.03}{0.06}{\eV} considering all bases. For the{\Btwog} state the energy difference between the CC3 and the CCSDT values is larger with \SIrange{0.18}{0.27}{\eV}. We can make a similar observation between the CC4 and the CCSDTQ values, for the {\Aoneg} state we have an energy difference of about \SI{0.01}{\eV} and this time we have smaller energy difference for the {\Btwog} with \SI{0.01}{\eV}. Then we discuss the multireference results and this time we were able to reach the aug-cc-pVQZ basis. Again we consider two different active spaces, four electrons in four orbitals and twelve electrons in twelve orbitals. If we compare the CASSCF(4,4) and CASPT2(4,4) values we can see large energy differences for the {\Aoneg} and {\Btwog} states, for the {\Aoneg} state we have an energy difference of about \SIrange{0.67}{0.74}{\eV} and \SIrange{1.65}{1.81}{\eV} for the {\Btwog}state. The energy difference is smaller for the triplet state with \SIrange{0.27}{0.31}{\eV}, we can notice that all the CASSCF(4,4) vertical energies are higher than the CASPT2(4,4) ones. Then we have the NEVPT2(4,4) methods (SC-NEVPT2(4,4) and PC-NEVPT2(4,4)) for which the vertical energies are quite similar, however we can notice that we have a different ordering of the {\Aoneg} and {\Btwog} states with {\Btwog} higher in energy than {\Aoneg} for the two NEVPT2(4,4) methods. Then we have the results for the same methods but with a larger active space. For the CASSCF(12,12) we have similar values for the triplet state with energy difference of about \SI{0.06}{\eV} for all bases but larger energy difference for the {\Aoneg} state with around \SIrange{0.28}{0.29}{\eV} and \SIrange{0.79}{0.81}{\eV} for the {\Btwog} state. Note that from CASSCF(4,4) to CASSCF(12,12) every vertical energy is smaller. We can make the contrary observation for the CASPT2 method where from CASPT2(4,4) to CASPT2(12,12) every vertical energy is larger. For the CASPT2(12,12) method we have similar values for the triplet state and for the {\Aoneg} state, considering all bases, with an energy difference of around \SIrange{0.05}{0.06}{\eV} and \SIrange{0.02}{0.05}{\eV} respectively. The energy difference is larger for the {\Btwog} state with about \SIrange{0.27}{0.29}{\eV}. Next we have the NEVPT2 methods, the first observation that we can make is that by increasing the size of the active space we retrieve the right ordering of the {\Aoneg} and {\Btwog} states. Again, every vertical energy is higher in the NEVPT2(12,12) case than for the NEVPT2(4,4) one. +Again, we have computed the vertical energies at the CCSDT and CCSDTQ level as well as their approximation the CC3 and CC4 methods, respectively. Note that for CC we started from a restricted Hartree-Fock (RHF) reference and in that case the ground state is the {\Aoneg}state, then the {\sBoneg} state is the single deexcitation and the {\Btwog} state is the double excitation from our ground state. For the CC3 method we do not have the vertical energies for the triplet state {\Atwog}. Considering all bases for the {\Aoneg}and {\Btwog} states we have an energy difference of about \SI{0.15}{\eV} and \SI{0.12}{\eV}, respectively. The CCSDT energies are close to the CC3 ones for the {\Aoneg} state with an energy difference of around \SIrange{0.03}{0.06}{\eV} considering all bases. For the{\Btwog} state the energy difference between the CC3 and the CCSDT values is larger with \SIrange{0.18}{0.27}{\eV}. We can make a similar observation between the CC4 and the CCSDTQ values, for the {\Aoneg} state we have an energy difference of about \SI{0.01}{\eV} and this time we have smaller energy difference for the {\Btwog} with \SI{0.01}{\eV}. Then we discuss the multireference results and this time we were able to reach the aug-cc-pVQZ basis. Again we consider two different active spaces, four electrons in four orbitals and twelve electrons in twelve orbitals. If we compare the CASSCF(4,4) and CASPT2(4,4) values we can see large energy differences for the {\Aoneg} and {\Btwog} states, for the {\Aoneg} state we have an energy difference of about \SIrange{0.67}{0.74}{\eV} and \SIrange{1.65}{1.81}{\eV} for the {\Btwog} state. The energy difference is smaller for the triplet state with \SIrange{0.27}{0.31}{\eV}, we can notice that all the CASSCF(4,4) vertical energies are higher than the CASPT2(4,4) ones. Then we have the NEVPT2(4,4) methods (SC-NEVPT2(4,4) and PC-NEVPT2(4,4)) for which the vertical energies are quite similar, however we can notice that we have a different ordering of the {\Aoneg} and {\Btwog} states with {\Btwog} higher in energy than {\Aoneg} for the two NEVPT2(4,4) methods. Then we have the results for the same methods but with a larger active space. For the CASSCF(12,12) we have similar values for the triplet state with energy difference of about \SI{0.06}{\eV} for all bases but larger energy difference for the {\Aoneg} state with around \SIrange{0.28}{0.29}{\eV} and \SIrange{0.79}{0.81}{\eV} for the {\Btwog} state. Note that from CASSCF(4,4) to CASSCF(12,12) every vertical energy is smaller. We can make the contrary observation for the CASPT2 method where from CASPT2(4,4) to CASPT2(12,12) every vertical energy is larger. For the CASPT2(12,12) method we have similar values for the triplet state and for the {\Aoneg} state, considering all bases, with an energy difference of around \SIrange{0.05}{0.06}{\eV} and \SIrange{0.02}{0.05}{\eV} respectively. The energy difference is larger for the {\Btwog} state with about \SIrange{0.27}{0.29}{\eV}. Next we have the NEVPT2 methods, the first observation that we can make is that by increasing the size of the active space we retrieve the right ordering of the {\Aoneg} and {\Btwog} states. Again, every vertical energy is higher in the NEVPT2(12,12) case than for the NEVPT2(4,4) one. %%% TABLE VI %%% \begin{squeezetable} \begin{table} \caption{ - Standard vertical excitation energies (with respect to the singlet $\text{X}\,{}^1B_{1g}$ ground state) of the {\Atwog}, {\Aoneg}, and {\Btwog} states of CBD at the \Dfour square-planar equilibrium geometry of the {\Atwog} state. + Standard vertical excitation energies (with respect to the singlet {\sBoneg} ground state) of the {\Atwog}, {\Aoneg}, and {\Btwog} states of CBD at the {\Dfour} square-planar equilibrium geometry of the {\Atwog} state. \label{tab:sf_D4h}} \begin{ruledtabular} \begin{tabular}{llrrr} @@ -716,7 +717,7 @@ SF-ADC(3) & 6-31+G(d) & $0.123$ & $1.650$ & $2.078$ \\ \begin{squeezetable} \begin{table} \caption{ - Standard vertical excitation energies (with respect to the singlet $\text{X}\,{}^1B_{1g}$ ground state) of the {\Atwog}, {\Aoneg}, and {\Btwog} states of CBD at the {\Dfour} square-planar equilibrium geometry of the {\Atwog} state. + Standard vertical excitation energies (with respect to the singlet {\sBoneg} ground state) of the {\Atwog}, {\Aoneg}, and {\Btwog} states of CBD at the {\Dfour} square-planar equilibrium geometry of the {\Atwog} state. \label{tab:D4h}} \begin{ruledtabular} \begin{tabular}{llrrr} @@ -791,7 +792,7 @@ CIPSI & 6-31+G(d) & $0.2010\pm 0.0030$ & $1.602\pm 0.007$ & $2.13\pm 0.04$ \\ \end{squeezetable} %%% %%% %%% %%% -Figure \ref{fig:D4h} shows the vertical energies of the studied excited states described in Tables \ref{tab:sf_D4h} and \ref{tab:D4h}. We see that all methods exhibit larger variations of the vertical energies due to the high symmetry of the \Dfour structure. For the triplet state most of the methods used are able to give a vertical energy close to the TBE one, we can although see that CASSCF, with the two different active spaces, show larger energy error than most of the DFT functionals used. Then for the, strongly multiconfigurational character, {\Aoneg} state we have a good description by the CC and multireference methods with the largest active space, except for CASSCF. The SF-ADC(2) and SF-ADC(3) schemes are also able to provide a good description of the {\Aoneg} state and even for the {\Btwog} we see that SF-ADC(2)-x give a worst result than SF-ADC(2). Multiconfigurational methods with the smallest active space do not demonstrate a good description of the two singlet excited states especially for the CASSCF method and for the the NEVPT2 methods that give the vertical energy of the {\Btwog} state below the {\Aoneg} one. Note that CASPT2 improve a lot the description of all the states compared to CASSCF. The various TD-DFT functionals are not able to describe correctly the two singlet excited states. +Figure \ref{fig:D4h} shows the vertical energies of the studied excited states described in Tables \ref{tab:sf_D4h} and \ref{tab:D4h}. We see that all methods exhibit larger variations of the vertical energies due to the high symmetry of the {\Dfour} structure. For the triplet state most of the methods used are able to give a vertical energy close to the TBE one, we can although see that CASSCF, with the two different active spaces, show larger energy error than most of the DFT functionals used. Then for the, strongly multiconfigurational character, {\Aoneg} state we have a good description by the CC and multireference methods with the largest active space, except for CASSCF. The SF-ADC(2) and SF-ADC(3) schemes are also able to provide a good description of the {\Aoneg} state and even for the {\Btwog} we see that SF-ADC(2)-x give a worst result than SF-ADC(2). Multiconfigurational methods with the smallest active space do not demonstrate a good description of the two singlet excited states especially for the CASSCF method and for the the NEVPT2 methods that give the vertical energy of the {\Btwog} state below the {\Aoneg} one. Note that CASPT2 improve a lot the description of all the states compared to CASSCF. The various TD-DFT functionals are not able to describe correctly the two singlet excited states. %%% FIGURE IV %%% \begin{figure*} @@ -810,7 +811,7 @@ Then we look at the vertical energy errors for the {\Dtwo} structure. First we c For the {\sBoneg} state of the {\Dtwo} structure we see that all the xc-functional underestimate the vertical excitation energy with energy differences of about \SIrange{0.35}{0.93}{\eV}. The ADC values are much closer to the TBE with energy differences around \SIrange{0.03}{0.09}{\eV}. Obviously, the CC vertical energies are close to the TBE one with around or less than \SI{0.01}{\eV} of energy difference. For the CASSCF(4,4) vertical energy we have a large difference of around \SI{1.42}{\eV} compared to the TBE value due to the lack of dynamical correlation in the CASSCF method. As previously seen the CAPT2(4,4) method correct this and we obtain a value of \SI{0.20}{\eV}. The others multireference methods in this active space give energy differences of around \SIrange{0.55}{0.76}{\eV} compared the the TBE reference. For the largest active space with twelve electrons in twelve orbitals we have an improvement of the vertical energies with \SI{0.72}{eV} of energy difference for the CASSCF(12,12) method and around 0.06 eV for the others multiconfigurational methods. -Then, for the {\Ag} state we obtain closer results of the SF-TD-DFT methods to the TBE than in the case of the {\sBoneg} state. Indeed, we have an energy difference of about \SIrange{0.01}{0.34}{\eV} for the {\Ag} state whereas we have \SIrange{0.35}{0.93}{\eV} for the {\sBoneg} state. The ADC schemes give the same error to the TBE value than for the other singlet state with \SI{0.02}{\eV} for the ADC(2) scheme and \SI{0.07}{\eV} for the ADC(3) one. The ADC(2)-x scheme provides a larger error with \SI{0.45}{\eV} of energy difference. Here, the CC methods manifest more variations with \SI{0.63}{\eV} for the CC3 value and \SI{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, \SIrange{0.03}{0.12}{\eV} compared to the TBE value. We can notice that CC3 and CCSDT provide larger energy errors for the {\Ag} state than for the {\sBoneg} state, this is due to the strong multiconfigurational character of the {\Ag} state whereas the {\sBoneg} 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 {\Ag} state. +Then, for the {\twoAg} state we obtain closer results of the SF-TD-DFT methods to the TBE than in the case of the {\sBoneg} state. Indeed, we have an energy difference of about \SIrange{0.01}{0.34}{\eV} for the {\twoAg} state whereas we have \SIrange{0.35}{0.93}{\eV} for the {\sBoneg} state. The ADC schemes give the same error to the TBE value than for the other singlet state with \SI{0.02}{\eV} for the ADC(2) scheme and \SI{0.07}{\eV} for the ADC(3) one. The ADC(2)-x scheme provides a larger error with \SI{0.45}{\eV} of energy difference. Here, the CC methods manifest more variations with \SI{0.63}{\eV} for the CC3 value and \SI{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, \SIrange{0.03}{0.12}{\eV} compared to the TBE value. We can notice that CC3 and CCSDT provide larger energy errors for the {\twoAg} state than for the {\sBoneg} state, this is due to the strong multiconfigurational character of the {\twoAg} state whereas the {\sBoneg} 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 {\twoAg} state. Finally we look at the vertical energy errors for the \Dfour structure. First, we consider the {\Atwog} state, the SF-TD-DFT methods give errors of about \SIrange{0.07}{1.6}{\eV} where the largest energy differences are provided by the hybrid functionals. The ADC schemes give similar error with around \SIrange{0.06}{1.1}{\eV} of energy difference. For the CC methods we have an energy error of \SI{0.06}{\eV} for CCSD and less than \SI{0.01}{\eV} for CCSDT. Then for the multireference methods with the four by four active space we have for CASSCF(4,4) \SI{0.29}{\eV} of error and \SI{0.02}{\eV} for CASPT2(4,4), again CASPT2 demonstrates its improvement compared to CASSCF. The other methods provide energy differences of about \SIrange{0.12}{0.13}{\eV}. A larger active space shows again an improvement with \SI{0.23}{\eV} of error for CASSCF(12,12) and around \SIrange{0.01}{0.04}{\eV} for the other multireference methods. CIPSI provides similar error with \SI{0.02}{\eV}. Then, we look at the {\Aoneg} state where the SF-TD-DFT shows large variations of error depending on the functionals, the energy error is about \SIrange{0.10}{1.03}{\eV}. The ADC schemes give better errors with around \SIrange{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 \SIrange{0.10}{0.16}{\eV} and CC4 provides really close energy to the TBE one with \SI{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 \SIrange{0.01}{0.73}{\eV} and \SIrange{0.02}{0.44}{\eV} respectively with the largest errors coming from the CASSCF method. Lastly, we look at the {\Btwog} state where we have globally larger errors. The SF-TD-DFT exhibits errors of \SIrange{0.43}{1.50}{\eV} whereas ADC schemes give errors of \SIrange{0.18}{0.30}{\eV}. CC3 and CCSDT provide energy differences of \SIrange{0.50}{0.69}{\eV} and the CC4 shows again close energy to the CCSDTQ TBE energy with \SI{0.01}{\eV} of error. The multireference methods give energy differences of \SIrange{0.38}{1.39}{\eV} and \SIrange{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 {\Dfour} 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 {\Aoneg} state than the {\Btwog} state, this can be related, as previously said in Subsubsection \ref{sec:D4h}, to the fact that {\Btwog} corresponds to a double excitation from the reference state. To obtain an improved description of the {\Btwog} 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 aug-cc-pVTZ basis. @@ -885,7 +886,7 @@ Finally we look at the vertical energy errors for the \Dfour structure. First, w %\begin{tabular}{*{1}{*{8}{l}}} & &\mc{3}{c}{{\Dtwo} excitation energies (eV)} & \mc{3}{c}{{\Dfour} excitation energies (eV)} \\ \cline{3-5} \cline{6-8} -Method & AB (\si{\kcalmol}) & {\tBoneg}(99\%) & {\sBoneg}(95\%)& {\Ag}(1\%) & $1\,{}^3A_{2g} $ & {\Aoneg} & {\Btwog} \\ +Method & AB (\si{\kcalmol}) & {\tBoneg}(99\%) & {\sBoneg}(95\%)& {\twoAg}(1\%) & {\Atwog} & {\Aoneg} & {\Btwog} \\ \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$ \\ @@ -950,7 +951,7 @@ In order to provide a benchmark of the AB and vertical energies we used coupled- With SF-TD-DFT the quality of the results are, of course, dependent on the functional but for the doubly excited states we have solid results. In SF-ADC we have very good results compared to the TBEs even for the doubly excited states, nevertheless the ADC(2)-x scheme give almost systematically worse results than the ADC(2)-s ones and using the ADC(3) scheme does not always provide better values. -The description of the excited states of the {\Dtwo} structure give rise to good agreement between the single reference and multiconfigurational methods due to the large T1 percentage of the first two excited states. When this percentage is much smaller as in the case of the doubly excited state {\Ag} the spin-flip methods show very good results within the ADC framework but even at the TD-DFT level spin-flip display good results compared to the TBE value. As said in the discussion, for the {\Dfour} geometry, the description of excited states is harder because of the strong multiconfigurational character where SF-TD-DFT can present more than 1 eV of error compared to the TBE. However, SF-ADC can show error of around \SIrange{0.1}{0.2}{\eV} which can be better than the multiconfigurational methods results. +The description of the excited states of the {\Dtwo} structure give rise to good agreement between the single reference and multiconfigurational methods due to the large T1 percentage of the first two excited states. When this percentage is much smaller as in the case of the doubly excited state {\twoAg} the spin-flip methods show very good results within the ADC framework but even at the TD-DFT level spin-flip display good results compared to the TBE value. As said in the discussion, for the {\Dfour} geometry, the description of excited states is harder because of the strong multiconfigurational character where SF-TD-DFT can present more than 1 eV of error compared to the TBE. However, SF-ADC can show error of around \SIrange{0.1}{0.2}{\eV} which can be better than the multiconfigurational methods results. %%%%%%%%%%%%%%%%%%%%%%%% \acknowledgements{ diff --git a/Manuscript/figure1.aux b/Manuscript/figure1.aux index b640121..f23e546 100644 --- a/Manuscript/figure1.aux +++ b/Manuscript/figure1.aux @@ -1,2 +1 @@ \relax -\gdef \@abspage@last{1} diff --git a/Manuscript/figure1.log b/Manuscript/figure1.log index 012515f..166748e 100644 --- a/Manuscript/figure1.log +++ b/Manuscript/figure1.log @@ -1,177 +1,180 @@ -This is pdfTeX, Version 3.141592653-2.6-1.40.22 (TeX Live 2021) (preloaded format=pdflatex 2021.11.19) 18 JAN 2022 11:56 +This is pdfTeX, Version 3.14159265-2.6-1.40.21 (TeX Live 2020) (preloaded format=pdflatex 2020.10.4) 22 MAR 2022 09:53 entering extended mode restricted \write18 enabled. - file:line:error style messages enabled. %&-line parsing enabled. **figure1.tex (./figure1.tex -LaTeX2e <2020-10-01> patch level 4 -L3 programming layer <2021-02-18> -(/usr/local/texlive/2021/texmf-dist/tex/latex/standalone/standalone.cls +LaTeX2e <2020-02-02> patch level 5 +L3 programming layer <2020-03-06> +(/usr/local/texlive/2020/texmf-dist/tex/latex/standalone/standalone.cls Document Class: standalone 2018/03/26 v1.3a Class to compile TeX sub-files stan dalone -(/usr/local/texlive/2021/texmf-dist/tex/latex/tools/shellesc.sty +(/usr/local/texlive/2020/texmf-dist/tex/latex/tools/shellesc.sty Package: shellesc 2019/11/08 v1.0c unified shell escape interface for LaTeX Package shellesc Info: Restricted shell escape enabled on input line 77. ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/iftex/ifluatex.sty +(/usr/local/texlive/2020/texmf-dist/tex/generic/iftex/ifluatex.sty Package: ifluatex 2019/10/25 v1.5 ifluatex legacy package. Use iftex instead. -(/usr/local/texlive/2021/texmf-dist/tex/generic/iftex/iftex.sty +(/usr/local/texlive/2020/texmf-dist/tex/generic/iftex/iftex.sty Package: iftex 2020/03/06 v1.0d TeX engine tests )) -(/usr/local/texlive/2021/texmf-dist/tex/latex/xkeyval/xkeyval.sty -Package: xkeyval 2020/11/20 v2.8 package option processing (HA) +(/usr/local/texlive/2020/texmf-dist/tex/latex/xkeyval/xkeyval.sty +Package: xkeyval 2014/12/03 v2.7a package option processing (HA) -(/usr/local/texlive/2021/texmf-dist/tex/generic/xkeyval/xkeyval.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/xkeyval/xkvutils.tex +(/usr/local/texlive/2020/texmf-dist/tex/generic/xkeyval/xkeyval.tex +(/usr/local/texlive/2020/texmf-dist/tex/generic/xkeyval/xkvutils.tex \XKV@toks=\toks15 \XKV@tempa@toks=\toks16 -(/usr/local/texlive/2021/texmf-dist/tex/generic/xkeyval/keyval.tex)) -\XKV@depth=\count179 +(/usr/local/texlive/2020/texmf-dist/tex/generic/xkeyval/keyval.tex)) +\XKV@depth=\count167 File: xkeyval.tex 2014/12/03 v2.7a key=value parser (HA) )) -\sa@internal=\count180 -\c@sapage=\count181 +\sa@internal=\count168 +\c@sapage=\count169 -(/usr/local/texlive/2021/texmf-dist/tex/latex/standalone/standalone.cfg +(/usr/local/texlive/2020/texmf-dist/tex/latex/standalone/standalone.cfg File: standalone.cfg 2018/03/26 v1.3a Default configuration file for 'standalon e' class ) -(/usr/local/texlive/2021/texmf-dist/tex/latex/base/article.cls -Document Class: article 2020/04/10 v1.4m Standard LaTeX document class -(/usr/local/texlive/2021/texmf-dist/tex/latex/base/size10.clo -File: size10.clo 2020/04/10 v1.4m Standard LaTeX file (size option) +(/usr/local/texlive/2020/texmf-dist/tex/latex/base/article.cls +Document Class: article 2019/12/20 v1.4l Standard LaTeX document class +(/usr/local/texlive/2020/texmf-dist/tex/latex/base/size10.clo +File: size10.clo 2019/12/20 v1.4l Standard LaTeX file (size option) ) -\c@part=\count182 -\c@section=\count183 -\c@subsection=\count184 -\c@subsubsection=\count185 -\c@paragraph=\count186 -\c@subparagraph=\count187 -\c@figure=\count188 -\c@table=\count189 +\c@part=\count170 +\c@section=\count171 +\c@subsection=\count172 +\c@subsubsection=\count173 +\c@paragraph=\count174 +\c@subparagraph=\count175 +\c@figure=\count176 +\c@table=\count177 \abovecaptionskip=\skip47 \belowcaptionskip=\skip48 -\bibindent=\dimen138 +\bibindent=\dimen134 ) -\sa@box=\box47 +\sa@box=\box45 ) -(/usr/local/texlive/2021/texmf-dist/tex/latex/endiagram/endiagram.sty -(/usr/local/texlive/2021/texmf-dist/tex/latex/l3kernel/expl3.sty -Package: expl3 2021-02-18 L3 programming layer (loader) +(/usr/local/texlive/2020/texmf-dist/tex/latex/endiagram/endiagram.sty +(/usr/local/texlive/2020/texmf-dist/tex/latex/l3kernel/expl3.sty +Package: expl3 2020-03-06 L3 programming layer (loader) -(/usr/local/texlive/2021/texmf-dist/tex/latex/l3backend/l3backend-pdftex.def -File: l3backend-pdftex.def 2021-03-18 L3 backend support: PDF output (pdfTeX) -\l__color_backend_stack_int=\count190 -\l__pdf_internal_box=\box48 -)) -(/usr/local/texlive/2021/texmf-dist/tex/latex/l3packages/xparse/xparse.sty -(/usr/local/texlive/2021/texmf-dist/tex/latex/l3packages/xparse/xparse-2020-10- -01.sty -(/usr/local/texlive/2021/texmf-dist/tex/latex/l3packages/xparse/xparse-generic. -tex))) -(/usr/local/texlive/2021/texmf-dist/tex/latex/l3packages/l3keys2e/l3keys2e.sty -Package: l3keys2e 2021-03-12 LaTeX2e option processing using LaTeX3 keys -) (/usr/local/texlive/2021/texmf-dist/tex/latex/pgf/frontendlayer/tikz.sty -(/usr/local/texlive/2021/texmf-dist/tex/latex/pgf/basiclayer/pgf.sty -(/usr/local/texlive/2021/texmf-dist/tex/latex/pgf/utilities/pgfrcs.sty -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/utilities/pgfutil-common.te +(/usr/local/texlive/2020/texmf-dist/tex/latex/l3backend/l3backend-pdfmode.def +File: l3backend-pdfmode.def 2020-03-12 L3 backend support: PDF mode +\l__kernel_color_stack_int=\count178 +\l__pdf_internal_box=\box46 +)) (/usr/local/texlive/2020/texmf-dist/tex/latex/l3packages/xparse/xparse.sty +Package: xparse 2020-03-06 L3 Experimental document command parser +\l__xparse_current_arg_int=\count179 +\g__xparse_grabber_int=\count180 +\l__xparse_m_args_int=\count181 +\l__xparse_v_nesting_int=\count182 +) +(/usr/local/texlive/2020/texmf-dist/tex/latex/l3packages/l3keys2e/l3keys2e.sty +Package: l3keys2e 2020-03-06 LaTeX2e option processing using LaTeX3 keys +) (/usr/local/texlive/2020/texmf-dist/tex/latex/pgf/frontendlayer/tikz.sty +(/usr/local/texlive/2020/texmf-dist/tex/latex/pgf/basiclayer/pgf.sty +(/usr/local/texlive/2020/texmf-dist/tex/latex/pgf/utilities/pgfrcs.sty +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/utilities/pgfutil-common.te x \pgfutil@everybye=\toks17 -\pgfutil@tempdima=\dimen139 -\pgfutil@tempdimb=\dimen140 +\pgfutil@tempdima=\dimen135 +\pgfutil@tempdimb=\dimen136 -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/utilities/pgfutil-common-li +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/utilities/pgfutil-common-li sts.tex)) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/utilities/pgfutil-latex.def -\pgfutil@abb=\box49 -) (/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/utilities/pgfrcs.code.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/pgf.revision.tex) -Package: pgfrcs 2020/12/27 v3.1.8b (3.1.8b) +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/utilities/pgfutil-latex.def +\pgfutil@abb=\box47 +(/usr/local/texlive/2020/texmf-dist/tex/latex/ms/everyshi.sty +Package: everyshi 2001/05/15 v3.00 EveryShipout Package (MS) )) -Package: pgf 2020/12/27 v3.1.8b (3.1.8b) +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/utilities/pgfrcs.code.tex +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/pgf.revision.tex) +Package: pgfrcs 2020/01/08 v3.1.5b (3.1.5b) +)) +Package: pgf 2020/01/08 v3.1.5b (3.1.5b) -(/usr/local/texlive/2021/texmf-dist/tex/latex/pgf/basiclayer/pgfcore.sty -(/usr/local/texlive/2021/texmf-dist/tex/latex/graphics/graphicx.sty -Package: graphicx 2020/09/09 v1.2b Enhanced LaTeX Graphics (DPC,SPQR) +(/usr/local/texlive/2020/texmf-dist/tex/latex/pgf/basiclayer/pgfcore.sty +(/usr/local/texlive/2020/texmf-dist/tex/latex/graphics/graphicx.sty +Package: graphicx 2019/11/30 v1.2a Enhanced LaTeX Graphics (DPC,SPQR) -(/usr/local/texlive/2021/texmf-dist/tex/latex/graphics/graphics.sty -Package: graphics 2020/08/30 v1.4c Standard LaTeX Graphics (DPC,SPQR) +(/usr/local/texlive/2020/texmf-dist/tex/latex/graphics/graphics.sty +Package: graphics 2019/11/30 v1.4a Standard LaTeX Graphics (DPC,SPQR) -(/usr/local/texlive/2021/texmf-dist/tex/latex/graphics/trig.sty +(/usr/local/texlive/2020/texmf-dist/tex/latex/graphics/trig.sty Package: trig 2016/01/03 v1.10 sin cos tan (DPC) ) -(/usr/local/texlive/2021/texmf-dist/tex/latex/graphics-cfg/graphics.cfg +(/usr/local/texlive/2020/texmf-dist/tex/latex/graphics-cfg/graphics.cfg File: graphics.cfg 2016/06/04 v1.11 sample graphics configuration ) Package graphics Info: Driver file: pdftex.def on input line 105. -(/usr/local/texlive/2021/texmf-dist/tex/latex/graphics-def/pdftex.def -File: pdftex.def 2020/10/05 v1.2a Graphics/color driver for pdftex +(/usr/local/texlive/2020/texmf-dist/tex/latex/graphics-def/pdftex.def +File: pdftex.def 2018/01/08 v1.0l Graphics/color driver for pdftex )) -\Gin@req@height=\dimen141 -\Gin@req@width=\dimen142 +\Gin@req@height=\dimen137 +\Gin@req@width=\dimen138 ) -(/usr/local/texlive/2021/texmf-dist/tex/latex/pgf/systemlayer/pgfsys.sty -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/systemlayer/pgfsys.code.tex -Package: pgfsys 2020/12/27 v3.1.8b (3.1.8b) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/utilities/pgfkeys.code.tex +(/usr/local/texlive/2020/texmf-dist/tex/latex/pgf/systemlayer/pgfsys.sty +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/systemlayer/pgfsys.code.tex +Package: pgfsys 2020/01/08 v3.1.5b (3.1.5b) +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/utilities/pgfkeys.code.tex \pgfkeys@pathtoks=\toks18 \pgfkeys@temptoks=\toks19 -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/utilities/pgfkeysfiltered.c +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/utilities/pgfkeysfiltered.c ode.tex \pgfkeys@tmptoks=\toks20 )) -\pgf@x=\dimen143 -\pgf@y=\dimen144 -\pgf@xa=\dimen145 -\pgf@ya=\dimen146 -\pgf@xb=\dimen147 -\pgf@yb=\dimen148 -\pgf@xc=\dimen149 -\pgf@yc=\dimen150 -\pgf@xd=\dimen151 -\pgf@yd=\dimen152 +\pgf@x=\dimen139 +\pgf@y=\dimen140 +\pgf@xa=\dimen141 +\pgf@ya=\dimen142 +\pgf@xb=\dimen143 +\pgf@yb=\dimen144 +\pgf@xc=\dimen145 +\pgf@yc=\dimen146 +\pgf@xd=\dimen147 +\pgf@yd=\dimen148 \w@pgf@writea=\write3 \r@pgf@reada=\read2 -\c@pgf@counta=\count191 -\c@pgf@countb=\count192 -\c@pgf@countc=\count193 -\c@pgf@countd=\count194 +\c@pgf@counta=\count183 +\c@pgf@countb=\count184 +\c@pgf@countc=\count185 +\c@pgf@countd=\count186 \t@pgf@toka=\toks21 \t@pgf@tokb=\toks22 \t@pgf@tokc=\toks23 -\pgf@sys@id@count=\count195 +\pgf@sys@id@count=\count187 -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/systemlayer/pgf.cfg -File: pgf.cfg 2020/12/27 v3.1.8b (3.1.8b) +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/systemlayer/pgf.cfg +File: pgf.cfg 2020/01/08 v3.1.5b (3.1.5b) ) Driver file for pgf: pgfsys-pdftex.def -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/systemlayer/pgfsys-pdftex.d +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/systemlayer/pgfsys-pdftex.d ef -File: pgfsys-pdftex.def 2020/12/27 v3.1.8b (3.1.8b) +File: pgfsys-pdftex.def 2020/01/08 v3.1.5b (3.1.5b) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/systemlayer/pgfsys-common-p +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/systemlayer/pgfsys-common-p df.def -File: pgfsys-common-pdf.def 2020/12/27 v3.1.8b (3.1.8b) +File: pgfsys-common-pdf.def 2020/01/08 v3.1.5b (3.1.5b) ))) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/systemlayer/pgfsyssoftpath. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/systemlayer/pgfsyssoftpath. code.tex -File: pgfsyssoftpath.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgfsyssoftpath@smallbuffer@items=\count196 -\pgfsyssoftpath@bigbuffer@items=\count197 +File: pgfsyssoftpath.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgfsyssoftpath@smallbuffer@items=\count188 +\pgfsyssoftpath@bigbuffer@items=\count189 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/systemlayer/pgfsysprotocol. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/systemlayer/pgfsysprotocol. code.tex -File: pgfsysprotocol.code.tex 2020/12/27 v3.1.8b (3.1.8b) -)) (/usr/local/texlive/2021/texmf-dist/tex/latex/xcolor/xcolor.sty +File: pgfsysprotocol.code.tex 2020/01/08 v3.1.5b (3.1.5b) +)) (/usr/local/texlive/2020/texmf-dist/tex/latex/xcolor/xcolor.sty Package: xcolor 2016/05/11 v2.12 LaTeX color extensions (UK) -(/usr/local/texlive/2021/texmf-dist/tex/latex/graphics-cfg/color.cfg +(/usr/local/texlive/2020/texmf-dist/tex/latex/graphics-cfg/color.cfg File: color.cfg 2016/01/02 v1.6 sample color configuration ) Package xcolor Info: Driver file: pdftex.def on input line 225. @@ -185,406 +188,406 @@ Package xcolor Info: Model `HSB' substituted by `hsb' on input line 1369. Package xcolor Info: Model `Gray' substituted by `gray' on input line 1370. Package xcolor Info: Model `wave' substituted by `hsb' on input line 1371. ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcore.code.tex -Package: pgfcore 2020/12/27 v3.1.8b (3.1.8b) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmath.code.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathcalc.code.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathutil.code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathparser.code.tex -\pgfmath@dimen=\dimen153 -\pgfmath@count=\count198 -\pgfmath@box=\box50 +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcore.code.tex +Package: pgfcore 2020/01/08 v3.1.5b (3.1.5b) +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmath.code.tex +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathcalc.code.tex +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathutil.code.tex) +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathparser.code.tex +\pgfmath@dimen=\dimen149 +\pgfmath@count=\count190 +\pgfmath@box=\box48 \pgfmath@toks=\toks24 \pgfmath@stack@operand=\toks25 \pgfmath@stack@operation=\toks26 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.code. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.code. tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.basic +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.basic .code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.trigo +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.trigo nometric.code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.rando +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.rando m.code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.compa +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.compa rison.code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.base. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.base. code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.round +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.round .code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.misc. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.misc. code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.integ +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathfunctions.integ erarithmetics.code.tex))) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmathfloat.code.tex -\c@pgfmathroundto@lastzeros=\count199 -)) (/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfint.code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcorepoints.co +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmathfloat.code.tex +\c@pgfmathroundto@lastzeros=\count191 +)) (/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfint.code.tex) +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcorepoints.co de.tex -File: pgfcorepoints.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgf@picminx=\dimen154 -\pgf@picmaxx=\dimen155 -\pgf@picminy=\dimen156 -\pgf@picmaxy=\dimen157 -\pgf@pathminx=\dimen158 -\pgf@pathmaxx=\dimen159 -\pgf@pathminy=\dimen160 -\pgf@pathmaxy=\dimen161 -\pgf@xx=\dimen162 -\pgf@xy=\dimen163 -\pgf@yx=\dimen164 -\pgf@yy=\dimen165 -\pgf@zx=\dimen166 -\pgf@zy=\dimen167 +File: pgfcorepoints.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgf@picminx=\dimen150 +\pgf@picmaxx=\dimen151 +\pgf@picminy=\dimen152 +\pgf@picmaxy=\dimen153 +\pgf@pathminx=\dimen154 +\pgf@pathmaxx=\dimen155 +\pgf@pathminy=\dimen156 +\pgf@pathmaxy=\dimen157 +\pgf@xx=\dimen158 +\pgf@xy=\dimen159 +\pgf@yx=\dimen160 +\pgf@yy=\dimen161 +\pgf@zx=\dimen162 +\pgf@zy=\dimen163 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcorepathconst +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcorepathconst ruct.code.tex -File: pgfcorepathconstruct.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgf@path@lastx=\dimen168 -\pgf@path@lasty=\dimen169 +File: pgfcorepathconstruct.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgf@path@lastx=\dimen164 +\pgf@path@lasty=\dimen165 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcorepathusage +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcorepathusage .code.tex -File: pgfcorepathusage.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgf@shorten@end@additional=\dimen170 -\pgf@shorten@start@additional=\dimen171 +File: pgfcorepathusage.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgf@shorten@end@additional=\dimen166 +\pgf@shorten@start@additional=\dimen167 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcorescopes.co +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcorescopes.co de.tex -File: pgfcorescopes.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgfpic=\box51 -\pgf@hbox=\box52 -\pgf@layerbox@main=\box53 -\pgf@picture@serial@count=\count266 +File: pgfcorescopes.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgfpic=\box49 +\pgf@hbox=\box50 +\pgf@layerbox@main=\box51 +\pgf@picture@serial@count=\count192 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcoregraphicst +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcoregraphicst ate.code.tex -File: pgfcoregraphicstate.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgflinewidth=\dimen172 +File: pgfcoregraphicstate.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgflinewidth=\dimen168 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcoretransform +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcoretransform ations.code.tex -File: pgfcoretransformations.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgf@pt@x=\dimen173 -\pgf@pt@y=\dimen174 -\pgf@pt@temp=\dimen175 +File: pgfcoretransformations.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgf@pt@x=\dimen169 +\pgf@pt@y=\dimen170 +\pgf@pt@temp=\dimen171 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcorequick.cod +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcorequick.cod e.tex -File: pgfcorequick.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: pgfcorequick.code.tex 2020/01/08 v3.1.5b (3.1.5b) ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcoreobjects.c +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcoreobjects.c ode.tex -File: pgfcoreobjects.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: pgfcoreobjects.code.tex 2020/01/08 v3.1.5b (3.1.5b) ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcorepathproce +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcorepathproce ssing.code.tex -File: pgfcorepathprocessing.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: pgfcorepathprocessing.code.tex 2020/01/08 v3.1.5b (3.1.5b) ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcorearrows.co +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcorearrows.co de.tex -File: pgfcorearrows.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgfarrowsep=\dimen176 +File: pgfcorearrows.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgfarrowsep=\dimen172 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcoreshade.cod +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcoreshade.cod e.tex -File: pgfcoreshade.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgf@max=\dimen177 -\pgf@sys@shading@range@num=\count267 -\pgf@shadingcount=\count268 +File: pgfcoreshade.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgf@max=\dimen173 +\pgf@sys@shading@range@num=\count193 +\pgf@shadingcount=\count194 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcoreimage.cod +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcoreimage.cod e.tex -File: pgfcoreimage.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: pgfcoreimage.code.tex 2020/01/08 v3.1.5b (3.1.5b) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcoreexternal. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcoreexternal. code.tex -File: pgfcoreexternal.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgfexternal@startupbox=\box54 +File: pgfcoreexternal.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgfexternal@startupbox=\box52 )) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcorelayers.co +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcorelayers.co de.tex -File: pgfcorelayers.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: pgfcorelayers.code.tex 2020/01/08 v3.1.5b (3.1.5b) ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcoretranspare +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcoretranspare ncy.code.tex -File: pgfcoretransparency.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: pgfcoretransparency.code.tex 2020/01/08 v3.1.5b (3.1.5b) ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcorepatterns. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcorepatterns. code.tex -File: pgfcorepatterns.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: pgfcorepatterns.code.tex 2020/01/08 v3.1.5b (3.1.5b) ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/basiclayer/pgfcorerdf.code. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/basiclayer/pgfcorerdf.code. tex -File: pgfcorerdf.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: pgfcorerdf.code.tex 2020/01/08 v3.1.5b (3.1.5b) ))) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/modules/pgfmoduleshapes.cod +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/modules/pgfmoduleshapes.cod e.tex -File: pgfmoduleshapes.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgfnodeparttextbox=\box55 +File: pgfmoduleshapes.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgfnodeparttextbox=\box53 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/modules/pgfmoduleplot.code. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/modules/pgfmoduleplot.code. tex -File: pgfmoduleplot.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: pgfmoduleplot.code.tex 2020/01/08 v3.1.5b (3.1.5b) ) -(/usr/local/texlive/2021/texmf-dist/tex/latex/pgf/compatibility/pgfcomp-version +(/usr/local/texlive/2020/texmf-dist/tex/latex/pgf/compatibility/pgfcomp-version -0-65.sty -Package: pgfcomp-version-0-65 2020/12/27 v3.1.8b (3.1.8b) -\pgf@nodesepstart=\dimen178 -\pgf@nodesepend=\dimen179 +Package: pgfcomp-version-0-65 2020/01/08 v3.1.5b (3.1.5b) +\pgf@nodesepstart=\dimen174 +\pgf@nodesepend=\dimen175 ) -(/usr/local/texlive/2021/texmf-dist/tex/latex/pgf/compatibility/pgfcomp-version +(/usr/local/texlive/2020/texmf-dist/tex/latex/pgf/compatibility/pgfcomp-version -1-18.sty -Package: pgfcomp-version-1-18 2020/12/27 v3.1.8b (3.1.8b) +Package: pgfcomp-version-1-18 2020/01/08 v3.1.5b (3.1.5b) )) -(/usr/local/texlive/2021/texmf-dist/tex/latex/pgf/utilities/pgffor.sty -(/usr/local/texlive/2021/texmf-dist/tex/latex/pgf/utilities/pgfkeys.sty -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/utilities/pgfkeys.code.tex) -) (/usr/local/texlive/2021/texmf-dist/tex/latex/pgf/math/pgfmath.sty -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmath.code.tex)) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/utilities/pgffor.code.tex -Package: pgffor 2020/12/27 v3.1.8b (3.1.8b) +(/usr/local/texlive/2020/texmf-dist/tex/latex/pgf/utilities/pgffor.sty +(/usr/local/texlive/2020/texmf-dist/tex/latex/pgf/utilities/pgfkeys.sty +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/utilities/pgfkeys.code.tex) +) (/usr/local/texlive/2020/texmf-dist/tex/latex/pgf/math/pgfmath.sty +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmath.code.tex)) +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/utilities/pgffor.code.tex +Package: pgffor 2020/01/08 v3.1.5b (3.1.5b) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/math/pgfmath.code.tex) -\pgffor@iter=\dimen180 -\pgffor@skip=\dimen181 +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/math/pgfmath.code.tex) +\pgffor@iter=\dimen176 +\pgffor@skip=\dimen177 \pgffor@stack=\toks27 \pgffor@toks=\toks28 )) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/frontendlayer/tikz/tikz.cod +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/frontendlayer/tikz/tikz.cod e.tex -Package: tikz 2020/12/27 v3.1.8b (3.1.8b) +Package: tikz 2020/01/08 v3.1.5b (3.1.5b) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/libraries/pgflibraryplothan +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/libraries/pgflibraryplothan dlers.code.tex -File: pgflibraryplothandlers.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgf@plot@mark@count=\count269 -\pgfplotmarksize=\dimen182 +File: pgflibraryplothandlers.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgf@plot@mark@count=\count195 +\pgfplotmarksize=\dimen178 ) -\tikz@lastx=\dimen183 -\tikz@lasty=\dimen184 -\tikz@lastxsaved=\dimen185 -\tikz@lastysaved=\dimen186 -\tikz@lastmovetox=\dimen187 -\tikz@lastmovetoy=\dimen188 -\tikzleveldistance=\dimen189 -\tikzsiblingdistance=\dimen190 -\tikz@figbox=\box56 -\tikz@figbox@bg=\box57 -\tikz@tempbox=\box58 -\tikz@tempbox@bg=\box59 -\tikztreelevel=\count270 -\tikznumberofchildren=\count271 -\tikznumberofcurrentchild=\count272 -\tikz@fig@count=\count273 +\tikz@lastx=\dimen179 +\tikz@lasty=\dimen180 +\tikz@lastxsaved=\dimen181 +\tikz@lastysaved=\dimen182 +\tikz@lastmovetox=\dimen183 +\tikz@lastmovetoy=\dimen184 +\tikzleveldistance=\dimen185 +\tikzsiblingdistance=\dimen186 +\tikz@figbox=\box54 +\tikz@figbox@bg=\box55 +\tikz@tempbox=\box56 +\tikz@tempbox@bg=\box57 +\tikztreelevel=\count196 +\tikznumberofchildren=\count197 +\tikznumberofcurrentchild=\count198 +\tikz@fig@count=\count199 -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/modules/pgfmodulematrix.cod +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/modules/pgfmodulematrix.cod e.tex -File: pgfmodulematrix.code.tex 2020/12/27 v3.1.8b (3.1.8b) -\pgfmatrixcurrentrow=\count274 -\pgfmatrixcurrentcolumn=\count275 -\pgf@matrix@numberofcolumns=\count276 +File: pgfmodulematrix.code.tex 2020/01/08 v3.1.5b (3.1.5b) +\pgfmatrixcurrentrow=\count266 +\pgfmatrixcurrentcolumn=\count267 +\pgf@matrix@numberofcolumns=\count268 ) -\tikz@expandcount=\count277 +\tikz@expandcount=\count269 -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie s/tikzlibrarytopaths.code.tex -File: tikzlibrarytopaths.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: tikzlibrarytopaths.code.tex 2020/01/08 v3.1.5b (3.1.5b) ))) -(/usr/local/texlive/2021/texmf-dist/tex/latex/siunitx/siunitx.sty +(/usr/local/texlive/2020/texmf-dist/tex/latex/siunitx/siunitx.sty Package: siunitx 2020/02/25 v2.8b A comprehensive (SI) units package -(/usr/local/texlive/2021/texmf-dist/tex/latex/amsmath/amstext.sty +(/usr/local/texlive/2020/texmf-dist/tex/latex/amsmath/amstext.sty Package: amstext 2000/06/29 v2.01 AMS text -(/usr/local/texlive/2021/texmf-dist/tex/latex/amsmath/amsgen.sty +(/usr/local/texlive/2020/texmf-dist/tex/latex/amsmath/amsgen.sty File: amsgen.sty 1999/11/30 v2.0 generic functions \@emptytoks=\toks29 -\ex@=\dimen191 +\ex@=\dimen187 )) -(/usr/local/texlive/2021/texmf-dist/tex/latex/tools/array.sty -Package: array 2020/10/01 v2.5c Tabular extension package (FMi) -\col@sep=\dimen192 -\ar@mcellbox=\box60 -\extrarowheight=\dimen193 +(/usr/local/texlive/2020/texmf-dist/tex/latex/tools/array.sty +Package: array 2019/08/31 v2.4l Tabular extension package (FMi) +\col@sep=\dimen188 +\ar@mcellbox=\box58 +\extrarowheight=\dimen189 \NC@list=\toks30 \extratabsurround=\skip49 \backup@length=\skip50 -\ar@cellbox=\box61 +\ar@cellbox=\box59 ) -\l__siunitx_tmp_box=\box62 -\l__siunitx_tmp_dim=\dimen194 -\l__siunitx_tmp_int=\count278 -\l__siunitx_number_mantissa_length_int=\count279 -\l__siunitx_number_uncert_length_int=\count280 -\l__siunitx_round_int=\count281 -\l__siunitx_process_decimal_int=\count282 -\l__siunitx_process_uncertainty_int=\count283 -\l__siunitx_process_fixed_int=\count284 -\l__siunitx_process_integer_min_int=\count285 -\l__siunitx_process_precision_int=\count286 -\l__siunitx_group_min_int=\count287 -\l__siunitx_angle_marker_box=\box63 -\l__siunitx_angle_unit_box=\box64 -\l__siunitx_angle_marker_dim=\dimen195 -\l__siunitx_angle_unit_dim=\dimen196 -\l__siunitx_unit_int=\count288 -\l__siunitx_unit_denominator_int=\count289 -\l__siunitx_unit_numerator_int=\count290 -\l__siunitx_unit_prefix_int=\count291 -\l__siunitx_unit_prefix_base_int=\count292 -\l__siunitx_unit_prefix_gram_int=\count293 -\l__siunitx_number_product_int=\count294 +\l__siunitx_tmp_box=\box60 +\l__siunitx_tmp_dim=\dimen190 +\l__siunitx_tmp_int=\count270 +\l__siunitx_number_mantissa_length_int=\count271 +\l__siunitx_number_uncert_length_int=\count272 +\l__siunitx_round_int=\count273 +\l__siunitx_process_decimal_int=\count274 +\l__siunitx_process_uncertainty_int=\count275 +\l__siunitx_process_fixed_int=\count276 +\l__siunitx_process_integer_min_int=\count277 +\l__siunitx_process_precision_int=\count278 +\l__siunitx_group_min_int=\count279 +\l__siunitx_angle_marker_box=\box61 +\l__siunitx_angle_unit_box=\box62 +\l__siunitx_angle_marker_dim=\dimen191 +\l__siunitx_angle_unit_dim=\dimen192 +\l__siunitx_unit_int=\count280 +\l__siunitx_unit_denominator_int=\count281 +\l__siunitx_unit_numerator_int=\count282 +\l__siunitx_unit_prefix_int=\count283 +\l__siunitx_unit_prefix_base_int=\count284 +\l__siunitx_unit_prefix_gram_int=\count285 +\l__siunitx_number_product_int=\count286 \c__siunitx_one_fill_skip=\skip51 \l__siunitx_table_unit_align_skip=\skip52 -\l__siunitx_table_exponent_dim=\dimen197 -\l__siunitx_table_integer_dim=\dimen198 -\l__siunitx_table_mantissa_dim=\dimen199 -\l__siunitx_table_marker_dim=\dimen256 -\l__siunitx_table_result_dim=\dimen257 -\l__siunitx_table_uncert_dim=\dimen258 -\l__siunitx_table_fill_pre_dim=\dimen259 -\l__siunitx_table_fill_post_dim=\dimen260 -\l__siunitx_table_fill_mid_dim=\dimen261 -\l__siunitx_table_pre_box=\box65 -\l__siunitx_table_post_box=\box66 -\l__siunitx_table_mantissa_box=\box67 -\l__siunitx_table_result_box=\box68 +\l__siunitx_table_exponent_dim=\dimen193 +\l__siunitx_table_integer_dim=\dimen194 +\l__siunitx_table_mantissa_dim=\dimen195 +\l__siunitx_table_marker_dim=\dimen196 +\l__siunitx_table_result_dim=\dimen197 +\l__siunitx_table_uncert_dim=\dimen198 +\l__siunitx_table_fill_pre_dim=\dimen199 +\l__siunitx_table_fill_post_dim=\dimen256 +\l__siunitx_table_fill_mid_dim=\dimen257 +\l__siunitx_table_pre_box=\box63 +\l__siunitx_table_post_box=\box64 +\l__siunitx_table_mantissa_box=\box65 +\l__siunitx_table_result_box=\box66 \l__siunitx_table_number_align_skip=\skip53 \l__siunitx_table_text_align_skip=\skip54 -(/usr/local/texlive/2021/texmf-dist/tex/latex/translator/translator.sty -Package: translator 2020-08-03 v1.12c Easy translation of strings in LaTeX +(/usr/local/texlive/2020/texmf-dist/tex/latex/translator/translator.sty +Package: translator 2019-05-31 v1.12a Easy translation of strings in LaTeX )) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie s/tikzlibrarycalc.code.tex -File: tikzlibrarycalc.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: tikzlibrarycalc.code.tex 2020/01/08 v3.1.5b (3.1.5b) ) Package: endiagram 2014/06/28 v0.1d Easy creation of potential energy curve dia grams -\l__endiagram_tmpa_dim=\dimen262 -\g__Endiagram_tmpa_dim=\dimen263 -\l__endiagram_points_int=\count295 -\g__endiagram_curve_int=\count296 -\g__endiagram_points_max_int=\count297 -\l__endiagram_points_counter_int=\count298 -\l__endiagram_niveau_counter_int=\count299 -\g__endiagram_ea_int=\count300 -\l__endiagram_ea_count_int=\count301 +\l__endiagram_tmpa_dim=\dimen258 +\g__Endiagram_tmpa_dim=\dimen259 +\l__endiagram_points_int=\count287 +\g__endiagram_curve_int=\count288 +\g__endiagram_points_max_int=\count289 +\l__endiagram_points_counter_int=\count290 +\l__endiagram_niveau_counter_int=\count291 +\g__endiagram_ea_int=\count292 +\l__endiagram_ea_count_int=\count293 ) -(/usr/local/texlive/2021/texmf-dist/tex/latex/pgfplots/pgfplots.sty -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplots.revision.tex) +(/usr/local/texlive/2020/texmf-dist/tex/latex/pgfplots/pgfplots.sty +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplots.revision.tex) Package: pgfplots 2020/02/29 v1.17 Data Visualization (1.17) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplots.code.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplotscore.code.tex +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplots.code.tex +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplotscore.code.tex \t@pgfplots@toka=\toks31 \t@pgfplots@tokb=\toks32 \t@pgfplots@tokc=\toks33 -\pgfplots@tmpa=\dimen264 -\c@pgfplots@coordindex=\count302 -\c@pgfplots@scanlineindex=\count303 +\pgfplots@tmpa=\dimen260 +\c@pgfplots@coordindex=\count294 +\c@pgfplots@scanlineindex=\count295 -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/sys/pgfplotssysgeneric +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/sys/pgfplotssysgeneric .code.tex)) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/libs/pgfplotslibrary.c +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/libs/pgfplotslibrary.c ode.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/oldpgfcompatib/pgfplot +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/oldpgfcompatib/pgfplot soldpgfsupp_loader.code.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/libraries/pgflibraryfpu.cod +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/libraries/pgflibraryfpu.cod e.tex)) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/util/pgfplotsutil.code +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/util/pgfplotsutil.code .tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/liststructure/pgfplots +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/liststructure/pgfplots liststructure.code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/liststructure/pgfplots +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/liststructure/pgfplots liststructureext.code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/liststructure/pgfplots +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/liststructure/pgfplots array.code.tex -\c@pgfplotsarray@tmp=\count304 +\c@pgfplotsarray@tmp=\count296 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/liststructure/pgfplots +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/liststructure/pgfplots matrix.code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/numtable/pgfplotstable +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/numtable/pgfplotstable shared.code.tex -\c@pgfplotstable@counta=\count305 +\c@pgfplotstable@counta=\count297 \t@pgfplotstable@a=\toks34 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/liststructure/pgfplots +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/liststructure/pgfplots deque.code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/util/pgfplotsbinary.co +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/util/pgfplotsbinary.co de.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/util/pgfplotsbinary.da +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/util/pgfplotsbinary.da ta.code.tex)) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/util/pgfplotsutil.verb +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/util/pgfplotsutil.verb .code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/libs/pgflibrarypgfplot +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/libs/pgflibrarypgfplot s.surfshading.code.tex -\c@pgfplotslibrarysurf@no=\count306 +\c@pgfplotslibrarysurf@no=\count298 -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/sys/pgflibrarypgfplots +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/sys/pgflibrarypgfplots .surfshading.pgfsys-pdftex.def))) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/util/pgfplotscolormap. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/util/pgfplotscolormap. code.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/util/pgfplotscolor.cod +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/util/pgfplotscolor.cod e.tex)) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplotsstackedplots.c +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplotsstackedplots.c ode.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplotsplothandlers.c +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplotsplothandlers.c ode.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplotsmeshplothandle +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplotsmeshplothandle r.code.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplotsmeshplotimage. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplotsmeshplotimage. code.tex))) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplots.scaling.code. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplots.scaling.code. tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplotscoordprocessin +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplotscoordprocessin g.code.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplots.errorbars.cod +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplots.errorbars.cod e.tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplots.markers.code. +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplots.markers.code. tex) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplotsticks.code.tex +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplotsticks.code.tex ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgfplots/pgfplots.paths.code.te +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgfplots/pgfplots.paths.code.te x) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie s/tikzlibrarydecorations.code.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/modules/pgfmoduledecoration +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/modules/pgfmoduledecoration s.code.tex -\pgfdecoratedcompleteddistance=\dimen265 -\pgfdecoratedremainingdistance=\dimen266 -\pgfdecoratedinputsegmentcompleteddistance=\dimen267 -\pgfdecoratedinputsegmentremainingdistance=\dimen268 -\pgf@decorate@distancetomove=\dimen269 -\pgf@decorate@repeatstate=\count307 -\pgfdecorationsegmentamplitude=\dimen270 -\pgfdecorationsegmentlength=\dimen271 +\pgfdecoratedcompleteddistance=\dimen261 +\pgfdecoratedremainingdistance=\dimen262 +\pgfdecoratedinputsegmentcompleteddistance=\dimen263 +\pgfdecoratedinputsegmentremainingdistance=\dimen264 +\pgf@decorate@distancetomove=\dimen265 +\pgf@decorate@repeatstate=\count299 +\pgfdecorationsegmentamplitude=\dimen266 +\pgfdecorationsegmentlength=\dimen267 ) -\tikz@lib@dec@box=\box69 +\tikz@lib@dec@box=\box67 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie s/tikzlibrarydecorations.pathmorphing.code.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/libraries/decorations/pgfli +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/libraries/decorations/pgfli brarydecorations.pathmorphing.code.tex)) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie s/tikzlibrarydecorations.pathreplacing.code.tex -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/libraries/decorations/pgfli +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/libraries/decorations/pgfli brarydecorations.pathreplacing.code.tex)) -\pgfplots@numplots=\count308 -\pgfplots@xmin@reg=\dimen272 -\pgfplots@xmax@reg=\dimen273 -\pgfplots@ymin@reg=\dimen274 -\pgfplots@ymax@reg=\dimen275 -\pgfplots@zmin@reg=\dimen276 -\pgfplots@zmax@reg=\dimen277 +\pgfplots@numplots=\count300 +\pgfplots@xmin@reg=\dimen268 +\pgfplots@xmax@reg=\dimen269 +\pgfplots@ymin@reg=\dimen270 +\pgfplots@ymax@reg=\dimen271 +\pgfplots@zmin@reg=\dimen272 +\pgfplots@zmax@reg=\dimen273 ) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/frontendlayer/tikz/librarie s/tikzlibraryplotmarks.code.tex -File: tikzlibraryplotmarks.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: tikzlibraryplotmarks.code.tex 2020/01/08 v3.1.5b (3.1.5b) -(/usr/local/texlive/2021/texmf-dist/tex/generic/pgf/libraries/pgflibraryplotmar +(/usr/local/texlive/2020/texmf-dist/tex/generic/pgf/libraries/pgflibraryplotmar ks.code.tex -File: pgflibraryplotmarks.code.tex 2020/12/27 v3.1.8b (3.1.8b) +File: pgflibraryplotmarks.code.tex 2020/01/08 v3.1.5b (3.1.5b) ))) (./figure1.aux) \openout1 = `figure1.aux'. @@ -602,36 +605,36 @@ LaTeX Font Info: Checking defaults for OMX/cmex/m/n on input line 6. LaTeX Font Info: ... okay on input line 6. LaTeX Font Info: Checking defaults for U/cmr/m/n on input line 6. LaTeX Font Info: ... okay on input line 6. - -(/usr/local/texlive/2021/texmf-dist/tex/context/base/mkii/supp-pdf.mkii + ABD: EveryShipout initializing macros +(/usr/local/texlive/2020/texmf-dist/tex/context/base/mkii/supp-pdf.mkii [Loading MPS to PDF converter (version 2006.09.02).] -\scratchcounter=\count309 -\scratchdimen=\dimen278 -\scratchbox=\box70 -\nofMPsegments=\count310 -\nofMParguments=\count311 +\scratchcounter=\count301 +\scratchdimen=\dimen274 +\scratchbox=\box68 +\nofMPsegments=\count302 +\nofMParguments=\count303 \everyMPshowfont=\toks35 -\MPscratchCnt=\count312 -\MPscratchDim=\dimen279 -\MPnumerator=\count313 -\makeMPintoPDFobject=\count314 +\MPscratchCnt=\count304 +\MPscratchDim=\dimen275 +\MPnumerator=\count305 +\makeMPintoPDFobject=\count306 \everyMPtoPDFconversion=\toks36 -) (/usr/local/texlive/2021/texmf-dist/tex/latex/epstopdf-pkg/epstopdf-base.sty +) (/usr/local/texlive/2020/texmf-dist/tex/latex/epstopdf-pkg/epstopdf-base.sty Package: epstopdf-base 2020-01-24 v2.11 Base part for package epstopdf Package epstopdf-base Info: Redefining graphics rule for `.eps' on input line 4 85. -(/usr/local/texlive/2021/texmf-dist/tex/latex/latexconfig/epstopdf-sys.cfg +(/usr/local/texlive/2020/texmf-dist/tex/latex/latexconfig/epstopdf-sys.cfg File: epstopdf-sys.cfg 2010/07/13 v1.3 Configuration of (r)epstopdf for TeX Liv e )) -(/usr/local/texlive/2021/texmf-dist/tex/latex/siunitx/siunitx-abbreviations.cfg -File: siunitx-abbreviations.cfg 2017/11/26 v2.7k siunitx: Abbreviated units -) -(/usr/local/texlive/2021/texmf-dist/tex/latex/translator/translator-basic-dicti +(/usr/local/texlive/2020/texmf-dist/tex/latex/translator/translator-basic-dicti onary-English.dict Dictionary: translator-basic-dictionary, Language: English ) +(/usr/local/texlive/2020/texmf-dist/tex/latex/siunitx/siunitx-abbreviations.cfg +File: siunitx-abbreviations.cfg 2017/11/26 v2.7k siunitx: Abbreviated units +) Package pgfplots Warning: running in backwards compatibility mode (unsuitable t ick labels; missing features). Consider writing \pgfplotsset{compat=1.17} into @@ -644,22 +647,22 @@ LaTeX Font Info: External font `cmex10' loaded for size (Font) <5> on input line 6. [1 -{/usr/local/texlive/2021/texmf-var/fonts/map/pdftex/updmap/pdftex.map}] +{/usr/local/texlive/2020/texmf-var/fonts/map/pdftex/updmap/pdftex.map}] (./figure1.aux) ) Here is how much of TeX's memory you used: - 25054 strings out of 478994 - 671157 string characters out of 5858183 - 1073899 words of memory out of 5000000 - 42168 multiletter control sequences out of 15000+600000 - 404839 words of font info for 33 fonts, out of 8000000 for 9000 + 25269 strings out of 480608 + 677395 string characters out of 5903630 + 1032025 words of memory out of 5000000 + 40789 multiletter control sequences out of 15000+600000 + 533747 words of font info for 30 fonts, out of 8000000 for 9000 1141 hyphenation exceptions out of 8191 - 128i,7n,130p,703b,649s stack positions out of 5000i,500n,10000p,200000b,80000s - -Output written on figure1.pdf (1 page, 35291 bytes). +Output written on figure1.pdf (1 page, 35294 bytes). PDF statistics: 27 PDF objects out of 1000 (max. 8388607) 19 compressed objects within 1 object stream diff --git a/Manuscript/figure1.pdf b/Manuscript/figure1.pdf index d49804e..c3a3e9d 100644 Binary files a/Manuscript/figure1.pdf and b/Manuscript/figure1.pdf differ diff --git a/Manuscript/figure1.synctex.gz b/Manuscript/figure1.synctex.gz index 5e0b0b9..4f8d778 100644 Binary files a/Manuscript/figure1.synctex.gz and b/Manuscript/figure1.synctex.gz differ diff --git a/Manuscript/figure1.tex b/Manuscript/figure1.tex index 4881c96..10ee1e7 100644 --- a/Manuscript/figure1.tex +++ b/Manuscript/figure1.tex @@ -62,7 +62,7 @@ %left% \draw (-2.5,-1.5) -- (-3.5,-1.5); - \node[] at (-4.1,-1.5) {$^1$A$_{g}$}; + \node[] at (-4.1,-1.5) {$1{}^1$A$_{g}$}; \draw[red] (-2.5,9.5) -- (-3.5,9.5); \node[red] at (-4.1,9.5) {$1~^3$B$_{1g}$}; @@ -70,14 +70,14 @@ %right% \draw (2.5,-1.5) -- (3.5,-1.5); - \node[] at (4.1,-1.5) {$^1$A$_{g}$}; + \node[] at (4.1,-1.5) {$1{}^1$A$_{g}$}; \draw[red] (2.5,9.5) -- (3.5,9.5); \node[red] at (4.1,9.5) {$1~^3$B$_{1g}$}; %above% \draw (-0.5,1.5) -- (0.5,1.5); - \node[] at (1,1.4) {$^1$B$_{1g}$}; + \node[] at (1,1.4) {$1{}^1$B$_{1g}$}; \draw[red] (-0.5,3.5) -- (0.5,3.5); \node[red] at (1.1,3.5) {$1~^3$A$_{2g}$};