This commit is contained in:
Pierre-Francois Loos 2022-03-20 18:53:38 +01:00
parent 2058b3cd37
commit 587c20d629
2 changed files with 65 additions and 50 deletions

View File

@ -193,60 +193,15 @@ deviation of the errors (SDE), as well as largest positive and negative deviatio
\end{figure} \end{figure}
%%% %%% %%% %%% %%% %%% %%% %%%
\begin{table*}
\caption{Statistical quantities (in eV), considering the 265 ``safe'' TBEs (out of 284) as reference, for various multi-reference methods.
All {error} values are in eV.}
\label{tab:stat}
\begin{ruledtabular}
\begin{tabular}{lccccccc}
&CASSCF\fnm[1] &CASPT2\fnm[1] &CASPT2\fnm[1] &CASPT3\fnm[1] &CASPT3\fnm[1] &SC-NEVPT2\fnm[2] &PC-NEVPT2\fnm[2]\\
& &(IPEA) &(NOIPEA) &(IPEA) &(NOIPEA)\\
\hline
MSE &$0.11$ &$0.06$ &$-0.26$ &$0,10$ &$0,05$ &$0,13$ &$0,09$\\
SDE &$0.58$ &$0.14$ &$0.21$ &$0.13$ &$0.13$ &$0.14$ &$0.14$\\
RMSE &$0.61$ &$0.16$ &$0.33$ &$0.17$ &$0.14$ &$0.19$ &$0.17$\\
MAE &$0.48$ &$0.11$ &$0.27$ &$0.11$ &$0.09$ &$0.15$ &$0.13$\\
Max($+$) &$2.15$ &$0.71$ &$0.30$ &$0.93$ &$0.79$ &$0.65$ &$0.46$\\
Max($-$) &$-1.18$ &$-0.32$ &$-1.02$ &$-0.28$ &$-0.36$ &$-0.38$ &$-0.57$\\
\end{tabular}
\end{ruledtabular}
\fnt[1]{Values from the present study.}
\fnt[2]{Values taken from Ref.~\onlinecite{Sarka_2022}.}
\end{table*}
\begin{table*}
\caption{MAEs determined for several subsets of transitions computed with various multi-reference methods.
Count is the number of excited states considered in each subset.}
\label{tab:stat_class}
\begin{ruledtabular}
\begin{tabular}{lcccccccc}
Transitions & Count &CASSCF\fnm[1] &CASPT2\fnm[1] &CASPT2\fnm[1] &CASPT3\fnm[1] &CASPT3\fnm[1] &SC-NEVPT2\fnm[2] &PC-NEVPT2\fnm[2]\\
& & &(IPEA) &(NOIPEA) &(IPEA) &(NOIPEA)\\
\hline
Singlet &174 &0.57 &0.14 &0.27 &0.14 &0.12 &0.16 &0.14\\
Triplet &110 &0.34 &0.07 &0.29 &0.07 &0.06 &0.13 &0.11\\
Valence &206 &0.45 &0.11 &0.33 &0.13 &0.10 &0.15 &0.12\\
Rydberg &78 &0.54 &0.13 &0.13 &0.08 &0.07 &0.14 &0.15\\
$n \to \pis$ &78 &0.44 &0.08 &0.44 &0.13 &0.10 &0.12 &0.10\\
$\pi \to \pis$ &119 &0.46 &0.12 &0.27 &0.13 &0.10 &0.18 &0.14\\
Double &9 &0.46 &0.11 &0.22 &0.12 &0.09 &0.14 &0.13\\
1-3 non-H atoms &39 &0.38 &0.07 &0.21 &0.06 &0.05 &0.10 &0.08\\
4 non-H atoms &94 &0.46 &0.11 &0.22 &0.12 &0.09 &0.14 &0.13\\
5-6 non-H atoms &151 &0.51 &0.12 &0.33 &0.12 &0.11 &0.17 &0.15\\
\end{tabular}
\end{ruledtabular}
\fnt[1]{Values from the present study.}
\fnt[2]{Values taken from Ref.~\onlinecite{Sarka_2022}.}
\end{table*}
A detailed discussion of each individual molecule can be found in Ref.~\onlinecite{Sarka_2022} where we also report relevant values from the literature. A detailed discussion of each individual molecule can be found in Ref.~\onlinecite{Sarka_2022} where we also report relevant values from the literature.
Here, we focus on global trends. Here, we focus on global trends.
The exhaustive list of CASPT2 and CASPT3 transitions can be found in Table \ref{tab:BigTab} and are represented in Fig.~\ref{fig:PT2_vs_PT3}. The exhaustive list of CASPT2 and CASPT3 transitions can be found in Table \ref{tab:BigTab} and are represented in Fig.~\ref{fig:PT2_vs_PT3}.
\begin{longtable*}{cllccccccccc} \begin{longtable*}{cllccccccccc}
\caption{Vertical excitation energies (in \si{\eV}) computed at various levels of theory. \caption{Vertical excitation energies (in \si{\eV}) computed with various multi-reference methods.
The reference TBEs of the QUEST database, their percentage of single excitations $\%T_1$ involved in the transition (computed at the CC3 level), their nature The reference TBEs of the QUEST database, their percentage of single excitations $\%T_1$ involved in the transition (computed at the CC3 level), their nature
(V and R stand for valence and Rydberg, respectively) are also reported. (V and R stand for valence and Rydberg, respectively) are also reported.
TBEs listed as ``safe'' are assumed to be chemically accurate (\ie, absolute error below \SI{0.05}{\eV}).
[F] indicates a fluorescence transition, \ie, a vertical transition energy computed from an excited-state equilibrium geometry. [F] indicates a fluorescence transition, \ie, a vertical transition energy computed from an excited-state equilibrium geometry.
\label{tab:BigTab}} \label{tab:BigTab}}
\\ \\
@ -551,6 +506,66 @@ The reference TBEs of the QUEST database, their percentage of single excitations
284 & &$^3A_2'(\pi,\pis)$ &V &97.6 &6.62 &Y &7.70 &6.62 &6.12 &6.76 &6.68\\ 284 & &$^3A_2'(\pi,\pis)$ &V &97.6 &6.62 &Y &7.70 &6.62 &6.12 &6.76 &6.68\\
\end{longtable*} \end{longtable*}
\begin{table*}
\caption{Statistical quantities (in eV), considering the 265 ``safe'' TBEs (out of 284) as reference, for various multi-reference methods.
All {error} values are in eV.}
\label{tab:stat}
\begin{ruledtabular}
\begin{tabular}{lccccccc}
&CASSCF\fnm[1] &CASPT2\fnm[1] &CASPT2\fnm[1] &CASPT3\fnm[1] &CASPT3\fnm[1] &SC-NEVPT2\fnm[2] &PC-NEVPT2\fnm[2]\\
& &(IPEA) &(NOIPEA) &(IPEA) &(NOIPEA)\\
\hline
MSE &$0.11$ &$0.06$ &$-0.26$ &$0,10$ &$0,05$ &$0,13$ &$0,09$\\
SDE &$0.58$ &$0.14$ &$0.21$ &$0.13$ &$0.13$ &$0.14$ &$0.14$\\
RMSE &$0.61$ &$0.16$ &$0.33$ &$0.17$ &$0.14$ &$0.19$ &$0.17$\\
MAE &$0.48$ &$0.11$ &$0.27$ &$0.11$ &$0.09$ &$0.15$ &$0.13$\\
Max($+$) &$2.15$ &$0.71$ &$0.30$ &$0.93$ &$0.79$ &$0.65$ &$0.46$\\
Max($-$) &$-1.18$ &$-0.32$ &$-1.02$ &$-0.28$ &$-0.36$ &$-0.38$ &$-0.57$\\
\end{tabular}
\end{ruledtabular}
\fnt[1]{Values from the present study.}
\fnt[2]{Values taken from Ref.~\onlinecite{Sarka_2022}.}
\end{table*}
\begin{table*}
\caption{MAEs determined for several subsets of transitions computed with various multi-reference methods.
Count is the number of excited states considered in each subset.}
\label{tab:stat_class}
\begin{ruledtabular}
\begin{tabular}{lcccccccc}
Transitions & Count &CASSCF\fnm[1] &CASPT2\fnm[1] &CASPT2\fnm[1] &CASPT3\fnm[1] &CASPT3\fnm[1] &SC-NEVPT2\fnm[2] &PC-NEVPT2\fnm[2]\\
& & &(IPEA) &(NOIPEA) &(IPEA) &(NOIPEA)\\
\hline
Singlet &174 &0.57 &0.14 &0.27 &0.14 &0.12 &0.16 &0.14\\
Triplet &110 &0.34 &0.07 &0.29 &0.07 &0.06 &0.13 &0.11\\
Valence &206 &0.45 &0.11 &0.33 &0.13 &0.10 &0.15 &0.12\\
Rydberg &78 &0.54 &0.13 &0.13 &0.08 &0.07 &0.14 &0.15\\
$n \to \pis$ &78 &0.44 &0.08 &0.44 &0.13 &0.10 &0.12 &0.10\\
$\pi \to \pis$ &119 &0.46 &0.12 &0.27 &0.13 &0.10 &0.18 &0.14\\
Double &9 &0.46 &0.11 &0.22 &0.12 &0.09 &0.14 &0.13\\
1-3 non-H atoms &39 &0.38 &0.07 &0.21 &0.06 &0.05 &0.10 &0.08\\
4 non-H atoms &94 &0.46 &0.11 &0.22 &0.12 &0.09 &0.14 &0.13\\
5-6 non-H atoms &151 &0.51 &0.12 &0.33 &0.12 &0.11 &0.17 &0.15\\
\end{tabular}
\end{ruledtabular}
\fnt[1]{Values from the present study.}
\fnt[2]{Values taken from Ref.~\onlinecite{Sarka_2022}.}
\end{table*}
\begin{table*}
\caption{CASPT2 and CASPT3 timings (in seconds) for a selection of systems and states.}
\begin{ruledtabular}
\begin{tabular}{llcccccccc}
System & State &Active &\# electrons &\# basis &\# CAS &\# contracted &\# uncontracted &CPU &CPU\\
& &Space & &functions &det. &config. &config. &CASPT2 &CASPT3 \\
\hline
Acetone &$^1A_2(n,\pis)$ &(6e,6o) &32 &322 &104 &$3.86 \times 10^6$ &$1.49 \times 10^8$ &12.50 &33.25\\
Pyrrole &$^1A_2(\pi,3s)$ &(6e,6o) &36 &345 &96 &$4.79 \times 10^6$ &$2.04 \times 10^8$ &13.24 &49.36\\
Imidazole &$^1A''(\pi,3s)$ &(8e,7o) &36 &322 &600 &$1.46 \times 10^7$ &$1.82 \times 10^9$ &193.93 &282.62\\
Pyrazine &$^1B_{3u}(n,\pis)$ &(10e,8o) &42 &368 &392 &$6.95 \times 10^6$ &$6.38 \times 10^8$ &29.58 &174.96\\
\end{tabular}
\end{ruledtabular}
\end{table*}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Conclusion} \section{Conclusion}

View File

@ -516,9 +516,9 @@ $^1A''[F](\pi,\pis)$ &($a'$:6,$a''$:3) &($A'$:1,$A''$:1) &0.72$^a$ &0.69$^a$ &0.
\end{tabular} \end{tabular}
\end{ruledtabular} \end{ruledtabular}
\flushleft \flushleft
$^a$Using reference $(8e,7o)$ active space including valence $\pi$, $\nN$ and $3s$ orbitals. $^a$Using reference (8e,7o) active space including valence $\pi$, $\nN$ and $3s$ orbitals.
$^b$Using reference $(6e,9o)$ active space including valence $\pi$ and four $3p_z$ orbitals. $^b$Using reference (6e,9o) active space including valence $\pi$ and four $3p_z$ orbitals.
$^c$Using reference $(8e,10o)$ active space including valence $\pi$, $\nN$, $3s$ and three $3p$ orbitals. $^c$Using reference (8e,10o) active space including valence $\pi$, $\nN$, $3s$ and three $3p$ orbitals.
\end{table*} \end{table*}
\begin{table*} \begin{table*}