diff --git a/Manuscript/seniority.tex b/Manuscript/seniority.tex index 6af9aa3..b010b5b 100644 --- a/Manuscript/seniority.tex +++ b/Manuscript/seniority.tex @@ -94,21 +94,20 @@ that best recover the correlation energy, understood as the energy difference be %that lead as fast as possible to the FCI limit. The most well-known and popular class of CI methods is excitation-based, -where one accounts for all determinants generated by exciting up to $d$ electrons from a given close-shell reference, which is usually the restricted HF solution, but does not have to. -In this way, the excitation degree $d$ parameter defines the hierarchy +where one accounts for all determinants generated by exciting up to $e$ electrons from a given close-shell reference, which is usually the restricted HF solution, but does not have to. +In this way, the excitation degree $e$ parameter defines the sequence CI with single excitations (CIS), CI with single and double excitations (CISD), CI with single, double, and triple excitations (CISDT), and so on. -% scaling is based on the excitation degree $d$. -The excitation-based CI hierarchy manages to quickly recover weak (dynamic) correlation effects, but struggles at strong (static) correlation regimes. +Excitation-based CI manages to quickly recover weak (dynamic) correlation effects, but struggles at strong (static) correlation regimes. Importantly, the number of determinants $N_{det}$ (which control the computational cost) scale polynomially with the number of electrons $N$ as $N^{2d}$. %This means that the contribution of higher excitations become progressively smaller. %In turn, seniority-based CI is specially targeted to describe static correlation. %\fk{Still have to work in this paragraph.} Alternatively, CI methods based on the seniority number \cite{Ring_1980} have been proposed \cite{Bytautas_2011}. -In short, the seniority number $\Omega$ is the number of unpaired electrons in a given determinant. -The seniority zero ($\Omega = 0$) sector has been shown to be the most important for static correlation, while higher sectors tend to contribute progressively less ~\cite{Bytautas_2011,Bytautas_2015,Alcoba_2014b,Alcoba_2014}. +In short, the seniority number $s$ is the number of unpaired electrons in a given determinant. +The seniority zero ($s = 0$) sector has been shown to be the most important for static correlation, while higher sectors tend to contribute progressively less ~\cite{Bytautas_2011,Bytautas_2015,Alcoba_2014b,Alcoba_2014}. % scaling -However, already at the CI$\Omega$0 level the number of determinants scale exponentially with $N$, since excitations of all excitation degrees $d$ are included. +However, already at the CIs0 level the number of determinants scale exponentially with $N$, since excitations of all excitation degrees $e$ are included. Therefore, despite the encouraging successes of seniority-based CI methods, their unfourable computational scaling restricts applications to very small systems. Besides CI, other methods that exploit the concpet of seniority number have been pursued. \cite{Henderson_2014,Chen_2015,Bytautas_2018} % Seniority Number in Valence Bond Theory @@ -122,15 +121,15 @@ When targeting static correlation, seniority-based CI methods tend to have a bet The latter class of methods, in contrast, are well-suited for recovering dynamic correlation, and only at polynomial cost with system size. % tackling Ideally, we aim for a method that captures most of both static and dynamic correlation, with as few determinants as possible. -With this goal in mind, we propose a new partitioning of the Hilbert space, named configuration interaction order (CIo). -It combines both the excitation degree $d$ and the seniority number $\Omega$ into one single parameter $o$ (order), +With this goal in mind, we propose a new partitioning of the Hilbert space, named hierarchy configuration interaction (hCI). +It combines both the excitation degree $e$ and the seniority number $s$ into one single hierarchy parameter $h$, \begin{equation} - \label{eq:o} - o = \frac{d+\Omega/2}{2}, + \label{eq:h} + h = \frac{e+s/2}{2}, \end{equation} which assumes half-integer values. % open-shell -Fig.~\ref{fig:allCI} shows how the Hilbert space is populated in excitation-based CI, seniority-based CI, and our hybrid CIo methods. +Fig.~\ref{fig:allCI} shows how the Hilbert space is populated in excitation-based CI, seniority-based CI, and our hybrid hCI methods. %%% FIG 1 %%% \begin{figure}[h!] @@ -148,12 +147,12 @@ Fig.~\ref{fig:allCI} shows how the Hilbert space is populated in excitation-base \end{subfigure} % \hfill \begin{subfigure}[b]{0.48\linewidth} -\includegraphics[width=\linewidth]{table_CIo} +\includegraphics[width=\linewidth]{table_hCI} \caption{Hybrid excitation-seniority CI.} - \label{fig:CIo} + \label{fig:hCI} \end{subfigure} - \caption{Partionining of the full Hilbert space into blocks of specific excitation degree $d$ (with respect to a closed-shell determinant) and seniority number $\Omega$. - Each of three classes of CI methods truncate this $d$-$\Omega$ map differently, and each color tone represents the added determinants at a given CI level.} + \caption{Partionining of the full Hilbert space into blocks of specific excitation degree $e$ (with respect to a closed-shell determinant) and seniority number $s$. + Each of three classes of CI methods truncate this $e$-$s$ map differently, and each color tone represents the added determinants at a given CI level.} \label{fig:allCI} \end{figure} %%% %%% %%% @@ -162,50 +161,50 @@ We have three key justifications for this new CI hierarchy. The first one is physical. We know that low degree excitations and low seniority sectors, when looked at individually, often have the most important contribution to the FCI expansion. %carry the most important weights. -By combining $d$ and $\Omega$ as is eq.~\ref{eq:o}, we ensure that both directions in the excitation-seniority map (see Fig.~\ref{fig:allCI}) will be contemplated. -Rather than filling the map top-bottom (as in excitation-based CI) or left-right (as in seniority-based CI), the CIo hierarchy fills it diagonally. +By combining $e$ and $s$ as is eq.~\ref{eq:h}, we ensure that both directions in the excitation-seniority map (see Fig.~\ref{fig:allCI}) will be contemplated. +Rather than filling the map top-bottom (as in excitation-based CI) or left-right (as in seniority-based CI), the hCI methods fills it diagonally. In this sense, we hope to recover dynamic correlation by moving right in the map (increasing the excitation degree while keeping a low seniority number), at the same time as static correlation, by moving down (increasing the seniority number while keeping a low excitation degree). %dynamic correlation is recovered with traditional CI. The second justification is computational. %Fig.~\ref{fig:scaling} also illustrates how the number of determinants within each block scales with the number of occupied orbitals $O$ and the number of virtual orbitals $V$. -In the CIo class of methods, each next level of theory accomodates additional determinants from different excitation-seniority sectors (each block of Fig.~\ref{fig:allCI}). -The key realization of the CIo hierarchy is that the number of additional determinants presents the same scaling with respect to $N$, for all excitation-seniority sectors entering at a given order $o$. +In the hCI class of methods, each next level of theory accomodates additional determinants from different excitation-seniority sectors (each block of Fig.~\ref{fig:allCI}). +The key realization behind hCI is that the number of additional determinants presents the same scaling with respect to $N$, for all excitation-seniority sectors entering at a given order $o$. %to $O$ and $V$, for all excitation-seniority sectors of a given order $o$. -%This computational realization represents the second justification for the introduction of the CIo method. -This further justifies the parameter $o$ as being the simple average between $d$ and $\Omega/2$. +%This computational realization represents the second justification for the introduction of the hCI method. +This further justifies the parameter $o$ as being the simple average between $e$ and $s/2$. -Each level of excitation-based CI has a CIo counterpart with the same scaling of $N_{det}$ with respect to $N$. -%However, CIo counts with additional half-integer levels of theory, with no parallel in excitation-based CI. -For example, in both CIo2 and CISD we have $N_{det} \sim N^4$, whereas in CIo3 and CISDT, $N_{det} \sim N^6$, and so on. -%the number of determinants of CIo2 and CISD scale as $O^2V^2$, those of CIo3 and CISDT scale as $O^3V^3$, and so on. -From this computational perspective, the CIo hierarchy can be seen as a more natural choice than the traditional excitation-based CI, -because if one can afford for, say, the $N_{det} \sim N^6$ cost of a CISDT calculation, than one can probably afford a CIo3 calculation, which has the same computational scaling. -Of course, in practice an integer-$o$ CIo method will have more determinants than its excitation-based counterpart (though the same scaling), -and thus one should first ensure whether including the lower-triangular blocks (going from CISDT to CIo3 in our example) +Each level of excitation-based CI has a hCI counterpart with the same scaling of $N_{det}$ with respect to $N$. +%However, hCI counts with additional half-integer levels of theory, with no parallel in excitation-based CI. +For example, in both hCI2 and CISD we have $N_{det} \sim N^4$, whereas in hCI3 and CISDT, $N_{det} \sim N^6$, and so on. +%the number of determinants of hCI2 and CISD scale as $O^2V^2$, those of hCI3 and CISDT scale as $O^3V^3$, and so on. +From this computational perspective, hCI can be seen as a more natural choice than the traditional excitation-based CI, +because if one can afford for, say, the $N_{det} \sim N^6$ cost of a CISDT calculation, than one can probably afford a hCI3 calculation, which has the same computational scaling. +Of course, in practice an integer-$o$ hCI method will have more determinants than its excitation-based counterpart (though the same scaling), +and thus one should first ensure whether including the lower-triangular blocks (going from CISDT to hCI3 in our example) is a better strategy than adding the next column (going from CISDT to CISDTQ). Therefore, here we decided to discuss the results in terms of the number of determinants, rather than the computational scaling, -which could make the comparison somewhat biased toward the CIo hierarchy. +which could make the comparison somewhat biased toward hCI. % -The lowest level in the CIo hierarchy (CIo1) parallels CIS of excitation-based CI. +The lowest level in hCI (hCI1) parallels CIS of excitation-based CI. However, the single excitations do not connect with the reference, at least for HF orbitals, and therefore CIS provides the same energy as HF. -In contrast, the paired doubles excitations of CIo1 do connect with the reference (as well as the singles, indireclty via the doubles). +In contrast, the paired doubles excitations of hCI1 do connect with the reference (as well as the singles, indireclty via the doubles). Therefore, while the HF-based lowest level of excitation-based CI (CIS) does not improve with respect to the mean-field HF wave function, -the CIo1 counterpart already represents a minimally correlated model, with the very favourable $N_{det} \sim N^2$ scaling. +the hCI1 counterpart already represents a minimally correlated model, with the very favourable $N_{det} \sim N^2$ scaling. %number of determinants scaling only as $OV$. % -In addition, CIo allows for half-integer orders $o$, with no parallel in excitation-based CI. +In addition, hCI allows for half-integer orders $o$, with no parallel in excitation-based CI. This gives extra flexibility in terms of choice of method. %when evaluating the computational cost and desired accuracy of a calculation. For a particular application with excitation-based CI, CISD might be too inaccurate, for example, while the price for the improved accuracy of CISDT might be too high. -With the CIo hierarchy, CIo2.5 represents an alternative, being more accurate than CIo2 and less expensive than CIo3. +hCI2.5 could represent an alternative, being more accurate than hCI2 and less expensive than hCI3. -Finally, the third justification for our CIo method is empirical and closely related to the computational motivation. +Finally, the third justification for our hCI method is empirical and closely related to the computational motivation. There are many possible ways to populate the Hilbert space starting from the a given reference determinant, and one can in principle formulate any systematic recipe that includes progressively more determinants. Besides a physical or computational perspective, the question of what makes for a good recipe can be framed empirically. -Does our CIo class of methods perform better than excitation-based or seniority-based CI, +Does our hCI class of methods perform better than excitation-based or seniority-based CI, in the sense of recovering most of the correlation energy with the least computational effort? A hybrid approach based on both excitation degree and seniority number has been proposed. \cite{Alcoba_2014,Raemdonck_2015,Alcoba_2018} @@ -215,8 +214,8 @@ For the union case, the number of determinants grows exponentially with $N$, while in the intersection approach the Hilbert space is filled rectangle-wise in our excitation-seniority map. In the latter case, the scaling of $N_{det}$ would be dominated by the rightmost bottom block. Bytautas et al.\cite{Bytautas_2015} explored a different hybrid scheme combining determinants from a complete active space and with a maximum seniority number. -In comparison to previous approaches, our hybrid CIo scheme has two key advantages. -First, it is defined by a single parameter that unifies excitation degree and seniority number (eq.\ref{eq:o}). +In comparison to previous approaches, our hybrid hCI scheme has two key advantages. +First, it is defined by a single parameter that unifies excitation degree and seniority number (eq.\ref{eq:h}). And second, each next level includes all classes of determinants sharing the same scaling with system size, as discussed before, thus keeping the method at a polynomial scaling. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @@ -224,7 +223,7 @@ And second, each next level includes all classes of determinants sharing the sam \label{sec:compdet} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -The CIo method was implemented in {\QP} via a straightforward adaptation of the \textit{configuration interaction using a perturbative selection made iteratively} (CIPSI) algorithm \cite{Huron_1973,Giner_2013,Giner_2015}, +The hCI method was implemented in {\QP} via a straightforward adaptation of the \textit{configuration interaction using a perturbative selection made iteratively} (CIPSI) algorithm \cite{Huron_1973,Giner_2013,Giner_2015}, by allowing only for determinants at a given order $o$. In practice, the CI energy is converged (within a chosen threshold of) with considerably fewer determinants than the formal number of determinants at a given $o$. The traditional excitation-based CI and the FCI calculations presented here were also performed with the CIPSI algorithm implemented in {\QP}. \cite{Huron_1973,Giner_2013,Giner_2015,Garniron_2019} @@ -238,7 +237,7 @@ Namely, whenever the eigenvalue of the orbital rotation Hessian is negative and then this gradient component is replaced by $g_0 |g_i|/g_i$. While we can never ensure that the obtained solutions are global minima in the orbital parameter space, we verified that in all cases surveyed here, the stationary solutions are real minima (rather than maxima or stationary points). -Here we assess the performance of the hierarchy of CIo methods against its excitation-based and seniority-based counterparts. +Here we assess the performance of the hCI methods against its excitation-based and seniority-based counterparts. To do so, we calculated the potential energy curves (PECs) for a total of 8 systems: \ce{HF}, \ce{F2}, \ce{N2}, %\ce{Be2}, \ce{H2O}, @@ -296,7 +295,7 @@ Thus, while the NPE probes the similarity regarding the shape of the PECs, the d %%% FIG 2 %%% \begin{figure}[h!] \includegraphics[width=\linewidth]{plot_stat} - \caption{Nonparallelity errors as function of the number of determinants, for the three classes of CI methods: seniority-based CI (blue), excitation-based CI (red), and our proposed hybrid CIo (green). + \caption{Nonparallelity errors as function of the number of determinants, for the three classes of CI methods: seniority-based CI (blue), excitation-based CI (red), and our proposed hybrid hCI (green). } \label{fig:plot_stat} \end{figure} @@ -307,31 +306,31 @@ Thus, while the NPE probes the similarity regarding the shape of the PECs, the d %For different CI approaches, Fig.~\ref{fig:N2_pes} shows PECs and their differences with respect to FCI, as well as the NPE and distance errors. %The associated differences with respect to the FCI result can be seen in the Supporting Information. -The main result contained in Fig.~\ref{fig:plot_stat} concerns the overall faster convergence of the CIo methods when compared to excitation-based and seniority-based CI methods. +The main result contained in Fig.~\ref{fig:plot_stat} concerns the overall faster convergence of the hCI methods when compared to excitation-based and seniority-based CI methods. This is observed both for single bond breaking (\ce{HF} and \ce{F2}) as well as the more challenging double (ethylene) and triple (\ce{N2}) bond breaking. The convergence with respect to the number of determinants is slower in the latter cases, irrespective of the class of CI methods, as would be expected. -But more importantly, the superiority of the CIo methods appear to be highlighted in the multiple bond break systems. +But more importantly, the superiority of the hCI methods appear to be highlighted in the multiple bond break systems. %Unless stated otherwise, from here on the performance of each method is probed by their NPE. Later we discuss other metrics. -For the four systems (more so for ethylene and \ce{N2}), CIo2 is better than CISD, two methods where the number of determinants scales as $N^4$. -CIo2.5 is better than CISDT, despite its lower computational cost, whereas CIo3 is much better than CISDT, and comparable in accuracy with CISDTQ. -Inspection of the PECs (see SI) reveal that the lower NPE in the CIo results stem mostly from the contribution of the dissociation region. -This result demonstrates the importance of higher-order excitations with low seniority number in this strong correlation regime, which are accounted for in CIo but not in excitation-based CI (for a given scaling with the number of determinants). +For the four systems (more so for ethylene and \ce{N2}), hCI2 is better than CISD, two methods where the number of determinants scales as $N^4$. +hCI2.5 is better than CISDT, despite its lower computational cost, whereas hCI3 is much better than CISDT, and comparable in accuracy with CISDTQ. +Inspection of the PECs (see SI) reveal that the lower NPE in the hCI results stem mostly from the contribution of the dissociation region. +This result demonstrates the importance of higher-order excitations with low seniority number in this strong correlation regime, which are accounted for in hCI but not in excitation-based CI (for a given scaling with the number of determinants). %The situation at the Franck-Condon region will be discussed later. % -Meanwhile, the first level of seniority-based CI (CI$\Omega$0, which is the same as doubly-occupied CI) tends to offer a rather low NPE when compare to the other CI methods with a similar number of determinants (CIo2.5 and CISDT). -However, convergence is clearly slower for the next levels in this hierarchy (CI$\Omega$2 and CI$\Omega$4), while excitation-based CI and specially CIo methods converge faster. +Meanwhile, the first level of seniority-based CI (CIs0, which is the same as doubly-occupied CI) tends to offer a rather low NPE when compare to the other CI methods with a similar number of determinants (hCI2.5 and CISDT). +However, convergence is clearly slower for the next levels (CIs2 and CIs4), while excitation-based CI and specially hCI methods converge faster. % -For the symmetric dissociation of linear \ce{H4} and \ce{H8} the performance of CIo and excitation-based CI are similar, both being superior to seniority-based CI. +For the symmetric dissociation of linear \ce{H4} and \ce{H8} the performance of hCI and excitation-based CI are similar, both being superior to seniority-based CI. -It is worth mentioning the surprisingly good performance of the CIo1 and CIo1.5 methods. +It is worth mentioning the surprisingly good performance of the hCI1 and hCI1.5 methods. For \ce{HF}, \ce{F2}, and ethylene, they presented lower NPEs than the much more expensive CISDT method, being slightly higher in the case of \ce{N2}. -For the same systems we also see the NPEs increase from CIo1.5 to CIo2, and decreasing to lower values only at the CIo3 level. +For the same systems we also see the NPEs increase from hCI1.5 to hCI2, and decreasing to lower values only at the hCI3 level. Both findings are not observed for \ce{H4} and \ce{H8}. -It seems that both the relative success of CIo1 and CIo1.5 methods as well as the relative worsening of the CIo2 method decrease as progressively more bonds are being broken (compare for instance \ce{F2}, \ce{N2}, and \ce{H8} in Fig.~\ref{fig:plot_stat}). +It seems that both the relative success of hCI1 and hCI1.5 methods as well as the relative worsening of the hCI2 method decrease as progressively more bonds are being broken (compare for instance \ce{F2}, \ce{N2}, and \ce{H8} in Fig.~\ref{fig:plot_stat}). This is because \fk{in progress...} -Even than, it is important to remember that even the CIo2 method remains superior to its excitation-based counterpart. +Even than, it is important to remember that even the hCI2 method remains superior to its excitation-based counterpart. %Whereas in excitation-based CI, the NPE always decrease as one moves to higher orders, @@ -343,29 +342,29 @@ In Fig.~\ref{fig:xe} and \ref{fig:freq}, we present the convergence of the equil with respect to the number of determinants, for the three types of CI approaches. %, vibrational frequencies, and dissociation energies, % -For \ce{F2}, the CIo method has an overall better convergence than the excitation-based CI counterpart, and much better than seniority-based CI. -The values oscillate around the FCI limit in CIo, whereas the convergence is monotonic in the two CI alternatives. -Interstingly, CIo1 and specially CIo1.5, two methods with a modest computational cost, provide very accurate equilibrium geometries and vibrational frequencies, +For \ce{F2}, the hCI method has an overall better convergence than the excitation-based CI counterpart, and much better than seniority-based CI. +The values oscillate around the FCI limit in hCI, whereas the convergence is monotonic in the two CI alternatives. +Interstingly, hCI1 and specially hCI1.5, two methods with a modest computational cost, provide very accurate equilibrium geometries and vibrational frequencies, Orbital optimization does not change the overall picture. -It does, however, lead to a more monotonic convergence in the case of CIo, which is not necessarily an advantage. -In particular, ooCIo1 and ooCIo1.5 are less accurate than their non-optimized counterparts. +It does, however, lead to a more monotonic convergence in the case of hCI, which is not necessarily an advantage. +In particular, oohCI1 and oohCI1.5 are less accurate than their non-optimized counterparts. % For \ce{HF} (results in the Supporting Information), -CIo and excitation-based CI are comparable to each other and superior to seniority-based CI, at least for HF orbitals. -Orbital optimization significantly improves the case for seniority-based CI, and leads to slightly better convergence for CIo with respect to excitation-based CI. +hCI and excitation-based CI are comparable to each other and superior to seniority-based CI, at least for HF orbitals. +Orbital optimization significantly improves the case for seniority-based CI, and leads to slightly better convergence for hCI with respect to excitation-based CI. -In the case of \ce{N2}, CIo and excitation-based CI present similar convergence behaviours, both being superior to seniority-based CI. -Also, CIo is slightly better than excitation-based CI for HF orbitals, whereas both are equally good with orbital optimization. -%the advantages of CIo are less evident, though stil present. +In the case of \ce{N2}, hCI and excitation-based CI present similar convergence behaviours, both being superior to seniority-based CI. +Also, hCI is slightly better than excitation-based CI for HF orbitals, whereas both are equally good with orbital optimization. +%the advantages of hCI are less evident, though stil present. % %A somewhat better convergence is also observed in the case of ethylene (see SI). The same conclusions hold for ethylene, \ce{H4}, and \ce{H8} (see SI). -Most of the times, the convergence of CIo either exceeds or is comparable to that of excitation-based CI. +Most of the times, the convergence of hCI either exceeds or is comparable to that of excitation-based CI. %%% FIG 3 %%% \begin{figure}[h!] \includegraphics[width=\linewidth]{xe} - \caption{Equilibrium geometries as function of the number of determinants, for the three classes of CI methods: seniority-based CI (blue), excitation-based CI (red), and our proposed hybrid CIo (green). + \caption{Equilibrium geometries as function of the number of determinants, for the three classes of CI methods: seniority-based CI (blue), excitation-based CI (red), and our proposed hybrid hCI (green). } \label{fig:xe} \end{figure} @@ -374,7 +373,7 @@ Most of the times, the convergence of CIo either exceeds or is comparable to tha %%% FIG 4 %%% \begin{figure}[h!] \includegraphics[width=\linewidth]{freq} - \caption{Vibrational frequencies (or force constants) as function of the number of determinants, for the three classes of CI methods: seniority-based CI (blue), excitation-based CI (red), and our proposed hybrid CIo (green). + \caption{Vibrational frequencies (or force constants) as function of the number of determinants, for the three classes of CI methods: seniority-based CI (blue), excitation-based CI (red), and our proposed hybrid hCI (green). } \label{fig:freq} \end{figure} @@ -395,9 +394,9 @@ Now we discuss the role of further optimizing the orbitals for each given CI met \fk{in progress...} -In CIo, each next step in the hierarchy brings different blocks of determinants which share the same computational scaling with respect to the number of electrons. +In hCI, each next step of the hierarchy brings different blocks of determinants which share the same computational scaling with respect to the number of electrons. -One of our key findings is that the NPE decreases faster with our hybrid CIo method than with either excitation-based or seniority-based CI. +One of our key findings is that the NPE decreases faster with our hybrid hCI method than with either excitation-based or seniority-based CI. One important conclusion is that orbital optimization is not necessarily a recommended strategy, depending on the properties one is interested in. While orbital optimization will certainly improve the energy at a particular geometry, such improvement may vary largely on the geometry, which may or may not decrease the NPE. @@ -406,10 +405,10 @@ and may imply in a significant computational burden (associated with the calcula In this sense, stepping up in the CI hierarchy might be a more straightforward and possibly cheaper alternative than optimizing the orbitals. One interesting possibility to explore is to first employ a low order CI method to optimize the orbitals, and then to employ this set of orbitals at a higher level of CI. -The CIo hierarchy of methods presented here offers several interesting possibilities to explore. -One is to investigate the performance of CIo or some adaptation of it for excited states. +The hCI pathways presented here offers several interesting possibilities to explore. +One is to investigate the performance of hCI or some adaptation of it for excited states. Another is to develop coupled cluster methods based on an analogous hybrid excitation-seniority truncation of the excitation operator. -One could also test the performance of CIo wave functions for Quantum Monte Carlo simulations. +One could also test the performance of hCI wave functions for Quantum Monte Carlo simulations. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{acknowledgements}