Search

\textsc{Quantum Package} is an open-source programming environment for quantum chemistry specially designed for wave function methods. Its main goal is the development of determinant-driven selected configuration interaction (sCI) methods and multi-reference second-order perturbation theory (PT2). The determinant-driven framework allows the programmer to include any arbitrary set of determinants in the reference space, hence providing greater methodological freedoms. The sCI method implemented in \textsc{Quantum Package} is based on the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm which complements the variational sCI energy with a PT2 correction. Additional external plugins have been recently added to perform calculations with multireference coupled cluster theory and range-separated density-functional theory. All the programs are developed with the IRPF90 code generator, which simplifies collaborative work and the development of new features. \textsc{Quantum Package} strives to allow easy implementation and experimentation of new methods, while making parallel computation as simple and efficient as possible on modern supercomputer architectures. Currently, the code enables, routinely, to realize runs on roughly 2\,000 CPU cores, with tens of millions of determinants in the reference space. Moreover, we have been able to push up to 12\,288 cores in order to test its parallel efficiency. In the present manuscript, we also introduce some key new developments: i) a renormalized second-order perturbative correction for efficient extrapolation to the full CI limit, and ii) a stochastic version of the CIPSI selection performed simultaneously to the PT2 calculation at no extra cost., Comment: 21 pages, 8 figures

Striving to define very accurate vertical transition energies, we perform both high-level coupled cluster (CC) calculations (up to CCSDTQP) and selected configuration interaction (sCI) calculations (up to several millions of determinants) for 18 small compounds. By systematically increasing the order of the CC expansion, the number of determinants in the CI expansion as well as the size of the one-electron basis set, we have been able to reach near full CI (FCI) quality transition energies. These calculations are carried out on CC3/aug-cc-pVTZ geometries, using a series of increasingly large atomic basis sets systematically including diffuse functions. In this way, we define a list of 110 transition energies for states of various characters to be used as references for further calculations. Benchmark transition energies are provided at the aug-pVTZ level as well as with additional basis set corrections, in order to obtain results close to the complete basis set limit. These reference data are used to benchmark a series of twelve excited-state wave function methods accounting for double and triple contributions, namely ADC(2), ADC(3), CIS(D), CIS(D$_\infty$), CC2, STEOM-CCSD, CCSD, CCSDR(3), CCSDT-3, CC3, CCSDT and CCSDTQ. It turns out that CCSDTQ yields a negligible difference with the extrapolated CI values with a mean absolute error as small as $0.01$ eV, whereas the coupled cluster approaches including iterative triples are also very accurate (mean absolute error of $0.03$ eV). Consequently, CCSDT-3 and CC3 can be used to define reliable benchmarks. This observation does not hold for ADC(3) that delivers quite large errors for this set of small compounds, with a clear tendency to overcorrect its second-order version, ADC(2). Finally, we discuss the possibility to use basis set extrapolation approaches so as to tackle more easily larger compounds., Comment: 57 pages, 2 figures

We present a comprehensive theoretical study of the physical phenomena that determine the relative energies of the three of the lowest electronic states of each of the square-planar copper complexes $\cucl$, $\cunh$ and $\cuwater$, and present a detailed analysis of the extent to which truncated configuration interaction (CI) and coupled cluster (CC) theories succeed in predicing the excitation energies. We find that ligand-metal charge transfer (CT) single excitations play a crucial role in the correct determination of the properties of these systems, even though the CT processes first occur at fourth order in perturbation theory, and propose a suitable choice of minimal active space for describing these systems with multi-reference theories. CCSD energy differences agree very well with near full CI values even though the $T_1$ diagnostics are large, which casts doubt on the usefulness of singles-amplitude based multi-reference diagnostics. CISD severely underestimates the excitation energies and the failure is a direct consequence of the size-inconsisency errors in CISD. Finally, we present reference values for the energy differences computed using explicitly correlated CCSD(T) and BCCD(T) theory., Comment: 33 pages, 14 figures

Selected configuration interaction (sCI) methods including second-order perturbative corrections provide near full CI (FCI) quality energies with only a small fraction of the determinants of the FCI space. Here, we introduce both a state-specific and a multi-state sCI method based on the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm. The present method revises the reference (internal) space under the effect of its interaction with the outer space via the construction of an effective Hamiltonian, following the shifted-Bk philosophy of Davidson and coworkers. In particular, the multi-state algorithm removes the storage bottleneck of the effective Hamiltonian via a low-rank factorization of the dressing matrix., Comment: 6 pages, 1 figure. Supplementary material available

In diffusion Monte Carlo (DMC) methods, the nodes (or zeroes) of the trial wave function dictate the magnitude of the fixed-node (FN) error. Within standard DMC implementations, they emanate from short multideterminant expansions, \textit{stochastically} optimized in the presence of a Jastrow factor. Here, following a recent proposal, we follow an alternative route by considering the nodes of selected configuration interaction (sCI) expansions built with the CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively) algorithm. In contrast to standard implementations, these nodes can be \textit{systematically} and \textit{deterministically} improved by increasing the size of the sCI expansion. The present methodology is used to investigate the properties of the transition metal sulfide molecule FeS. This apparently simple molecule has been shown to be particularly challenging for electronic structure theory methods due to the proximity of two low-energy quintet electronic states of different spatial symmetry. In particular, we show that, at the triple-zeta basis set level, all sCI results --- including those extrapolated at the full CI (FCI) limit --- disagree with experiment, yielding an electronic ground state of $^{5}\Sigma^+$ symmetry. Performing FN-DMC simulation with sCI nodes, we show that the correct $^{5}\Delta$ ground state is obtained if sufficiently large expansions are used. Moreover, we show that one can systematically get accurate potential energy surfaces and reproduce the experimental dissociation energy as well as other spectroscopic constants., Comment: 8 pages, 2 figure and 4 tables

A hybrid stochastic-deterministic approach for computing the second-order perturbative contribution $E^{(2)}$ within multireference perturbation theory (MRPT) is presented. The idea at the heart of our hybrid scheme --- based on a reformulation of $E^{(2)}$ as a sum of elementary contributions associated with each determinant of the MR wave function --- is to split $E^{(2)}$ into a stochastic and a deterministic part. During the simulation, the stochastic part is gradually reduced by dynamically increasing the deterministic part until one reaches the desired accuracy. In sharp contrast with a purely stochastic MC scheme where the error decreases indefinitely as $t^{-1/2}$ (where $t$ is the computational time), the statistical error in our hybrid algorithm displays a polynomial decay $\sim t^{-n}$ with $n=3-4$ in the examples considered here. If desired, the calculation can be carried on until the stochastic part entirely vanishes. In that case, the exact result is obtained with no error bar and no noticeable computational overhead compared to the fully-deterministic calculation. The method is illustrated on the F$_2$ and Cr$_2$ molecules. Even for the largest case corresponding to the Cr$_2$ molecule treated with the cc-pVQZ basis set, very accurate results are obtained for $E^{(2)}$ for an active space of (28e,176o) and a MR wave function including up to $2 \times 10^7$ determinants., Comment: 8 pages, 5 figures

The present paper introduces a new multi-reference perturbation approach developed at second order, based on a Jeziorsky-Mokhorst expansion using individual Slater determinants as perturbers. Thanks to this choice of perturbers, an effective Hamiltonian may be built, allowing for the dressing of the Hamiltonian matrix within the reference space, assumed here to be a CAS-CI. Such a formulation accounts then for the coupling between the static and dynamic correlation effects. With our new definition of zeroth-order energies, these two approaches are strictly size-extensive provided that local orbitals are used, as numerically illustrated here and formally demonstrated in the appendix. Also, the present formalism allows for the factorization of all double excitation operators, just as in internally contracted approaches, strongly reducing the computational cost of these two approaches with respect to other determinant-based perturbation theories. The accuracy of these methods has been investigated on ground-state potential curves up to full dissociation limits for a set of six molecules involving single, double and triple bond breaking. The spectroscopic constants obtained with the present methods are found to be in very good agreement with the full configuration interaction (FCI) results. As the present formalism does not use any parameter or numerically unstable operation, the curves obtained with the two methods are smooth all along the dissociation path., Comment: 4 figures, 18 pages

A central difficulty of state-specific Multi-Reference Coupled Cluster (MR-CC) formalisms concerns the definition of the amplitudes of the single and double excitation operators appearing in the exponential wave operator. If the reference space is a complete active space (CAS) the number of these amplitudes is larger than the number of singly and doubly excited determinants on which one may project the eigenequation, and one must impose additional conditions. The present work first defines a state-specific reference-independent operator $\hat{\tilde{T}}^m$ which acting on the CAS component of the wave function $|\Psi_0^m \rangle$ maximizes the overlap between $(1+\hat{\tilde{T}}^m)|\Psi_0^m \rangle$ and the eigenvector of the CAS-SD CI matrix $|\Psi_{\rm CAS-SD}^m \rangle$. This operator may be used to generate approximate coefficients of the Triples and Quadruples, and a dressing of the CAS-SD CI matrix, according to the intermediate Hamiltonian formalism. The process may be iterated to convergence. As a refinement towards a strict Coupled Cluster formalism, one may exploit reference-independent amplitudes provided by $(1+\hat{\tilde{T}}^m)|\Psi_0^m \rangle$ to define a reference-dependent operator $\hat{T}^m$ by fitting the eigenvector of the (dressed) CAS-SD CI matrix. The two variants, which are internally uncontracted, give rather similar results. The new MR-CC version has been tested on the ground state potential energy curves of 6 molecules (up to triple-bond breaking) and a two excited states. The non-parallelism error with respect to the Full-CI curves is of the order of 1 m$E_{\rm h}$., Comment: 11 pages