Search

The computation of magnetizability tensors using gauge-including atomic orbitals is discussed in the context of Cholesky decomposition for the two-electron repulsion integrals with a focus on the involved doubly differentiated integrals. Three schemes for their handling are suggested: the first exploits the DF aspect of Cholesky decomposition, the second uses expressions obtained by differentiating the CD expression for the unperturbed two electron integrals, while the third addresses the issue that the first two schemes are not able to represent the doubly differentiated integrals with arbitrary accuracy. This scheme uses a separate Cholesky decomposition for the cross terms in the doubly differentiated two-electron integrals. Test calculations reveal that all three schemes are able to represent the integrals with similar accuracy and yield indistinguishable results for the values of the computed magnetizability tensor elements. Thus, we recommend our first scheme which has the lowest computational cost for routine computations. The applicability of our CD schemes is further shown in large-scale Hartree-Fock calculations of the magnetizability tensor of coronene (C24H12) with a doubly polarized triple-zeta basis consisting of 684 basis functions.

Local electronic-structure methods in quantum chemistry operate on the ability to compress electron correlations more efficiently in a basis of spatially localized molecular orbitals than in a parent set of canonical orbitals. However, many typical choices of localized orbitals tend to be related by selected, near-exact symmetry operations whenever a molecule belongs to a point group, a feature which remains largely unexploited in most local correlation methods. The present Letter demonstrates how to leverage a recent unitary protocol for enforcing symmetry properties among localized orbitals to yield a high-accuracy estimate of the exact ground-state correlation energy of benzene ($D_{6h}$) in correlation-consistent polarized basis sets of both double- and triple-$\zeta$ quality. Through an initial application to many-body expanded full configuration interaction (MBE-FCI) theory, we show how molecular point-group symmetry can lead to computational savings that are inversely proportional to the order of the point group in a manner generally applicable to the acceleration of modern local correlation methods., Comment: 13+5 pages, 5 figures. SI as an ancillary file

The unfolding of molecular complexes or biomolecules under the influence of external mechanical forces can routinely be simulated with atomistic resolution. To obtain a match of the characteristic time scales with those of experimental force spectroscopy, often coarse graining procedures are employed. Here, building on a previous study, we apply the adaptice resolution scheme (AdResS) to force probe molecular dynamics (FPMD) simulations using two model systems as examples. One system is the previously investigated calix[4]arene dimer that shows reversible one-step unfolding and the other example is provided by a small peptide, a $\beta$-alanine octamer in methanol solvent. The mechanical unfolding of this peptide proceeds via a metastable intermediate and therefore represents a first step towards a complex unfolding pathway. In addition to increasing the complexity of the relevant conformational changes we study the impact of the methodology used for coarse graining. Apart from a standard technique, the iterative Boltzmann inversion method, we apply an ideal gas approximation and therefore we replace the solvent by a non-interacting system of spherical particles. In all cases we find excellent agreement between the results of FPMD simulations performed fully atomistically and the AdResS simulations also in the case of fast pulling. This holds for all details of the unfolding pathways like the distributions of the characteristic forces and also the sequence of hydrogen-bond opening in case of the $\beta$-alanine octamer. Therefore, the methodology is very well suited to simulate the mechanical unfolding of systems of experimental relevance also in the presence of protic solvents., Comment: 22 pages, 11 figures

We present an implementation for the use of Cholesky decomposition (CD) of two-electron integrals within the spin-free Dirac-Coulomb (SFDC) scheme that enables to perform high-accuracy coupled-cluster (CC) calculations at costs almost comparable to those of their non-relativistic counterparts. While for non-relativistic CC calculations atomic-orbital (AO) based algorithms, due to their significantly reduced disk-space requirements, are the key to efficient large-scale computations, such algorithms are less advantageous in the SFDC case due to their increased computational cost on that case. Here, molecular-orbital (MO) based algorithms exploiting the CD of the two-electron integrals allow to reduce disk-space requirements, and lead to computational cost in the CC step that are more or less the same as in the non-relativistic case. The only remaining overhead in a CD-SFDC-CC calculation are due to the need to compute additional two-electron integrals, the somewhat higher cost of the Hartree-Fock calculation in the SFDC case, and additional cost in the transformation of the Cholesky vectors from the AO to the MO representation. However, these additional costs typically amount to less than 5-15 % of the total wall time and are thus acceptable. We illustrate the efficiency of our CD scheme for SFDC-CC calculations on a series of illustrative calculations for the X(CO)$_4$ molecules with X = Ni, Pd, Pt.

We present a wide-reaching revamp of the generalized many-body expanded full configuration interaction (MBE-FCI) method. First, we outline how to automatize the selection of reference active spaces whereby the inherent bias introduced through a manual identification is reduced, also within the context of traditional complete active space methods. Second, we allow for the use of compact orbital clusters as expansion objects, which works to circumvent the unfavorable scaling with the number of orbitals included in the space complementary to the reference orbitals. Finally, we present a new algorithm for ensuring that many-body expansions can be efficiently terminated while conservatively accounting for resulting errors. These developments are all tested on a variety of molecular systems and different orbital representations to illustrate the abilities of our algorithm to produce correlation energies within predetermined error bounds, significantly broadening the overall applicability of the MBE-FCI method., Comment: 29+7 pages, 8 figures, 1 table. SI as an ancillary file

We present a novel implementation of the complete active space self-consistent field (CASSCF) method that makes use of the many-body expanded full configuration interaction (MBE-FCI) method to incrementally approximate electronic structures within large active spaces. On the basis of a hybrid first-order algorithm employing both Super-CI and quasi-Newton strategies for the optimization of molecular orbitals, we demonstrate both computational efficacy and high accuracy of the resulting MBE-CASSCF method. We assess the performance of our implementation on a set of established numerical tests before applying MBE-CASSCF in the investigation of the triplet-quintet spin gap of iron(II) porphyrin with active spaces as large as 50 electrons in 50 orbitals., Comment: 27+8 pages, 6 figures. SI as an ancillary file

A novel implementation of the coupled-cluster singles and doubles (CCSD) approach is presented that is specifically tailored for the treatment of large, symmetric systems. It fully exploits Abelian point-group symmetry and the use of the Cholesky decomposition of the two-electron repulsion integrals. In accordance with modern CCSD algorithms, we propose two alternative strategies for the computation of the so-called particle-particle ladder term. The code is driven towards the optimal choice depending on the available hardware resources. As a large-scale application, we computed the frozen-core correlation energy of buckminsterfullerene (C$_{60}$) with a polarized valence triple-zeta basis set (240 correlated electrons in 1740 orbitals).

A Cholesky decomposition (CD)-based implementation of relativistic two-component coupled-cluster (CC) and equation-of-motion CC (EOM-CC) methods using an exact two-component Hamiltonian augmented with atomic-mean-field spin-orbit integrals (the X2CAMF scheme) is reported. The present CD-based implementation of X2CAMF-CC and EOM-CC methods employs atomic-orbital-based algorithms to avoid the construction of two-electron integrals and intermediates involving three and four virtual indices. Our CD-based implementation extends the applicability of X2CAMF-CC and EOM-CC methods to medium-sized molecules with the possibility to correlate around 1000 spinors. Benchmark calculations for uranium-containing small molecules have been performed to assess the dependence of the CC results on the Cholesky threshold. A Cholesky threshold of $10^{-4}$ is shown to be sufficient to maintain chemical accuracy. Example calculations to illustrate the capability of the CD-based relativistic CC methods are reported for the bond-dissociation energy of the uranium hexafluoride molecule, UF$_6$, with up to quadruple-zeta basis sets, and the lowest excitation energy in solvated uranyl ion [UO$_2^{2+}$(H$_2$O)$_{12}$].

We present an implementation of coupled-perturbed complete active space self-consistent field (CP-CASSCF) theory for the computation of nuclear magnetic resonance chemical shifts using gauge-including atomic orbitals and Cholesky decomposed two-electron integrals. The CP-CASSCF equations are solved using a direct algorithm where the magnetic Hessian matrix-vector product is expressed in terms of one-index transformed quantities. Numerical tests on systems with up to about 1300 basis functions provide information regarding both the computational efficiency and limitations of our implementation.

A rigorous analysis is carried out concerning the use of Cholesky decomposition (CD) of two-electron integrals in the case of quantum-chemical calculations with finite or perturbative magnetic fields and gauge-including atomic orbitals. We investigate in particular how permutational symmetry can be accounted for in such calculations and how this symmetry can be exploited to reduce the computational requirements. A modified CD procedure is suggested for the finite-field case that roughly halves the memory demands for the storage of the Cholesky vectors. The resulting symmetry of the Cholesky vectors also enables savings in the computational costs. For the derivative two-electron integrals in case of a perturbative magnetic field we derive CD expressions by means of a first-order Taylor expansion of the corresponding finite magnetic-field formulas with the field-free case as reference point. The perturbed Cholesky vectors are shown to be antisymmetric (as already proposed by Burger et al. (J. Chem. Phys., 155, 074105 (2021))) and the corresponding expressions enable significant savings in the required integral evaluations (by a factor of about four) as well as in the actual construction of the Cholesky vectors (by means of a two-step procedure similar to the one presented by Folkestad et al. (J. Chem. Phys., 150, 194112 (2019)) and Zhang et al. (J. Phys. Chem. A, 125, 4258-4265 (2021))). Numerical examples with cases involving several hundred basis functions verify our suggestions concerning CD in case of finite and perturbative magnetic fields.

We present an interface of the wavefunction-based quantum-chemical software CFOUR to the multiscale modeling framework MiMiC. Electrostatic embedding of the quantummechanical (QM) part is achieved by analytic evaluation of one-electron integrals in CFOUR, while the rest of the QM/MM operations are treated according to the previous MiMiC-based QM/MM implementation. Long-range electrostatic interactions are treated by a multipole expansion of the potential from the QM electron density to reduce the computational cost without loss of accuracy. Testing on model water/water systems, we verified that the CFOUR interface to MiMiC is robust, guaranteeing fast convergence of the SCF cycles and optimal conservation of the energy during the integration of the equations of motion. Finally, we verified that the CFOUR interface to MiMiC is compatible with the use of a QM/QM multiple time-step algorithm, which effectively reduces the cost of AIMD or QM/MM-MD simulations using higher level wavefunction-based approaches compared to cheaper density-functional theory based ones. The new wavefunction-based AIMD and QM/MM-MD implementation was tested and validated for a large number of wavefunction approaches, including Hartree-Fock and post-Hartree-Fock methods like Moller-Plesset, coupled cluster, and complete active space self-consistent field.

We present the implementation of a quadratically convergent Self-consistent field (QCSCF) algorithm based on an adaptive trust-radius optimization scheme for restricted open-shell Hartree-Fock (ROHF), restricted Hartree-Fock (RHF), and unrestricted Hartree-Fock (UHF) references. The algorithm can exploit Cholesky decomposition (CD) of the two-electron integrals to allow calculations on larger systems. The most important feature of the QCSCF code lies in its black-box nature -- probably the most important quality desired by a generic user. As shown for pilot applications, it does not require one to tune the self-consistent field (SCF) parameters (damping, Pulay's DIIS, and other similar techniques) in difficult-to-converge molecules. Also, it can be used to obtain a very thigh convergence with extended basis set - a situation often needed when computing high-order molecular properties - where the standard SCF algorithm starts to oscillate. Nevertheless, trouble may appear even with a QCSCF solver. In this respect, we discuss what can go wrong, focusing on the multiple UHF solutions of ortho-benzyne

We report on a formulation and implementation of a scheme to compute NMR shieldings at second-order Moller-Plesset (MP2) perturbation theory using gauge-including atomic orbitals (GIAOs) to ensure gauge-origin independence and Cholesky decomposition (CD) to handle unperturbed as well as perturbed two-electron integrals. We investigate the accuracy of the CD for the derivatives of the two-electron integrals with respect to an external magnetic field as well as for the computed NMR shieldings, before we illustrate the applicability of our CD based GIAO-MP2 scheme in calculations involving up to about one hundred atoms and more than one thousand basis functions., Comment: 12 pages, 2 figures

Since in coupled-cluster (CC) theory ground-state and excitation energies are eigenvalues of a non-Hermitian matrix, these energies can in principle take on complex values. In this paper we discuss the appearance of complex energy values in CC calculations from a mathematical perspective. We analyze the behaviour of the eigenvalues of Hermitian matrices that are perturbed (in a non-Hermitian manner) by a real parameter. Based on these results we show that for CC calculations with real-valued Hamiltonian matrices the ground-state energy generally takes a real value. Furthermore, we show that in the case of real-valued Hamiltonian matrices complex excitation energies only occur in the context of conical intersections. In such a case, unphysical consequences are encountered such as a wrong dimension of the intersection seam, large numerical deviations from full configuration-interaction (FCI) results, and the square-root-like behaviour of the potential surfaces near the conical intersection. In the case of CC calculations with complex-valued Hamiltonian matrix elements, it turns out that complex energy values are to be expected for ground and excited states when no symmetry is present. We confirm the occurrence of complex energies by sample calculations using a six-state model and by CC calculations for the H2O molecule in a strong magnetic field. We furthermore show that symmetry can prevent the occurrence of complex energy values. Lastly, we demonstrate that in most cases the real part of the complex energy values provides a very good approximation to the FCI energy., Comment: 19 pages, 8 figures

More than 130 pure rotational transitions of $^{46}$TiO, $^{47}$TiO, $^{48}$TiO, $^{49}$TiO, $^{50}$TiO, and $^{48}$Ti$^{18}$O are recorded using a high-resolution mm-wave supersonic jet spectrometer in combination with a laser ablation source. For the first time a mass-independent Dunham-like analysis is performed encompassing rare titanium monoxide isotopologues, and are compared to results from high-accuracy quantum-chemical calculations. The obtained parametrization reveals for titanium monoxide effects due to deviations from the Born-Oppenheimer approximation. Additionally, the dominant titanium properties enable an insight into the electronic structure of TiO by analyzing its hyperfine interactions. Further, based on the mass-independent analysis, the frequency positions of the pure rotational transitions of the short lived rare isotopologue $^{44}$TiO are predicted with high accuracy, i.e., on a sub-MHz uncertainty level. This allows for dedicated radio-astronomical searches of this species in core-collapse environments of supernovae., Comment: 29 pages, 3 figures

In this contribution, we present the implementation of a second-order CASSCF algorithm in conjunction with the Cholesky decomposition of the two-electron repulsion integrals. The algorithm, called Norm-Extended Optimization, guarantees convergence of the optimization, but it involves the full Hessian of the wavefunction and is therefore computationally expensive. Coupling the second-order procedure with the Cholesky decomposition leads to a significant reduction in the computational cost, reduced memory requirements, and an improved parallel performance. As a result, CASSCF calculations of larger molecular systems become possible as a routine task. The performance of the new implementation is illustrated by means of benchmark calculations on molecules of increasing size, with up to about 3000 basis functions and 14 active orbitals.

Molecular simulations of the forced unfolding and refolding of biomolecules or molecular complexes allow to gain important kinetic, structural and thermodynamic information about the folding process and the underlying energy landscape. In force probe molecular dynamics (FPMD) simulations, one pulls one end of the molecule with a constant velocity in order to induce the relevant conformational transitions. Since the extended configuration of the system has to fit into the simulation box together with the solvent such simulations are very time consuming. Here, we apply a hybrid scheme in which the solute is treated with atomistic resolution and the solvent molecules far away from the solute are described in a coarse-grained manner. We use the adaptive resolution scheme (AdResS) that has very successfully been applied to various examples of equilibrium simulations. We perform FPMD simulations using AdResS on a well studied system, a dimer formed from mechanically interlocked calixarene capsules. The results of the multiscale simulations are compared to all-atom simulations of the identical system and we observe that the size of the region in which atomistic resolution is required depends on the pulling velocity, i.e. the particular non-equilibrium situation. For large pulling velocities a larger all atom region is required. Our results show that multiscale simulations can be applied also in the strong non-equilibrium situations that the system experiences in FPMD simulations., Comment: 22 pages, 12 figures

The recent many-body expanded full configuration interaction (MBE-FCI) method is reviewed by critically assessing its advantages and drawbacks in the context of contemporary near-exact electronic structure theory. Besides providing a succinct summary of the history of MBE-FCI to date within a generalized and unified theoretical setting, its finer algorithmic details are discussed alongside our optimized computational implementation of the theory. A selected few of the most recent applications of MBE-FCI are revisited, before we close by outlining its future research directions as well as its place among modern near-exact wave function-based methods., Comment: 18+8 pages, 5 figures. Invited review to WIREs Comput Mol Sci

The recently proposed many-body expanded full configuration interaction (MBE-FCI) method is extended to excited states and static first-order properties different from total, ground state correlation energies. Results are presented for excitation energies and (transition) dipole moments of two prototypical, heteronuclear diatomics---LiH and MgO---in augmented correlation consistent basis sets of up to quadruple-$\zeta$ quality. Given that MBE-FCI properties are evaluated without recourse to a sampled wave function and the storage of corresponding reduced density matrices, the memory overhead associated with the calculation of general first-order properties only scales with the dimension of the desired property. In combination with the demonstrated performance, the present developments are bound to admit a wide range of future applications by means of many-body expanded treatments of electron correlation., Comment: 22 pages, 2 figures, 2 tables. SI added as ancillary file

We report on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground state energy of the benzene molecule in a standard correlation-consistent basis set of double-$\zeta$ quality. As a broad international endeavour, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire offered by modern electronic structure theory. In our assessment, the evaluated high-level methods are all found to qualitatively agree on a final correlation energy, with most methods yielding an estimate of the FCI value around $-863$ m$E_{\text{H}}$. However, we find the root-mean-square deviation of the energies from the studied methods to be considerable (1.3 m$E_{\text{H}}$), which in light of the acclaimed performance of each of the methods for smaller molecular systems clearly displays the challenges faced in extending reliable, near-exact correlation methods to larger systems. While the discrepancies exposed by our study thus emphasize the fact that the current state-of-the-art approaches leave room for improvement, we still expect the present assessment to provide a valuable community resource for benchmark and calibration purposes going forward., Comment: 29 pages, 1 figure, 2 tables. SI as an ancillary file

Very accurate transition frequencies of HC$_5$N were determined between 5.3 and 21.4 GHz with a Fourier transform microwave spectrometer. The molecules were generated by passing a mixture of HC$_3$N and C$_2$H$_2$ highly diluted in neon through a discharge valve followed by supersonic expansion into the Fabry-Perot cavity of the spectrometer. The accuracies of the data permitted us to improve the experimental $^{14}$N nuclear quadrupole coupling parameter considerably and the first experimental determination of the $^{14}$N nuclear spin-rotation parameter. The transition frequencies are also well suited to determine in astronomical observations the local speed of rest velocities in molecular clouds with high fidelity. The same setup was used to study HC$_7$N, albeit with modest improvement of the experimental $^{14}$N nuclear quadrupole coupling parameter. Quantum chemical calculations were carried out to determine $^{14}$N nuclear quadrupole and spin-rotation coupling parameters of HC$_5$N, HC$_7$N, and related molecules. These calculations included evaluation of vibrational and relativistic corrections to the non-relativistic equilibrium quadrupole coupling parameters; their considerations improved the agreement between calculated and experimental values substantially., Comment: 9 pages; J. Mol. Spectrosc., accepted

Facilitated by a rigorous partitioning of a molecular system's orbital basis into two fundamental subspaces - a reference and an expansion space, both with orbitals of unspecified occupancy - we generalize our recently introduced many-body expanded full configuration interaction (MBE-FCI) method to allow for electron-rich model and molecular systems dominated by both weak and strong correlation to be addressed. By employing minimal or even empty reference spaces, we show through calculations on the one-dimensional Hubbard model with up to 46 lattice sites, the chromium dimer, and the benzene molecule how near-exact results may be obtained in a entirely unbiased manner for chemical and physical problems of not only academic, but also applied chemical interest. Given the massive parallelism and overall accuracy of the resulting method, we argue that generalized MBE-FCI theory possesses an immense potential to yield near-exact correlation energies for molecular systems of unprecedented size, composition, and complexity in the years to come., Comment: 19 pages, 1 TOC figure, 4 main figures, 1 SI attached as an ancillary file. UPDATE-v2 (Dec. 19): Discussion of Hubbard results now focussed on comparison again DMRG reference numbers

In this second part of our series on the recently proposed many-body expanded full configuration interaction (MBE-FCI) method, we introduce the concept of multideterminantal expansion references. Through theoretical arguments and numerical validations, the use of this class of starting points is shown to result in a focussed compression of the MBE decomposition of the FCI energy, thus allowing chemical problems dominated by strong correlation to be addressed by the method. The general applicability and performance enhancements of MBE-FCI are verified for standard stress tests such as the bond dissociations in H$_2$O, N$_2$, C$_2$, and a linear H$_{10}$ chain. Furthermore, the benefits of employing a multideterminantal expansion reference in accelerating calculations of high accuracy are discussed, with an emphasis on calculations in extended basis sets. As an illustration of this latter quality of the MBE-FCI method, results for H$_2$O and C$_2$ in basis sets ranging from double- to pentuple-$\zeta$ quality are presented, demonstrating near-ideal parallel scaling on up to almost $25000$ processing units., Comment: 41 pages, 4 tables, 10 figures, 1 SI attached as an ancillary file

Over the course of the past few decades, the field of computational chemistry has managed to manifest itself as a key complement to more traditional lab-oriented chemistry. This is particularly true in the wake of the recent renaissance of full configuration interaction (FCI)-level methodologies, albeit only if these can prove themselves sufficiently robust and versatile to be routinely applied to a variety of chemical problems of interest. In the present series of works, performance and feature enhancements of one such avenue towards FCI-level results for medium to large one-electron basis sets, the recently introduced many-body expanded full configuration interaction (MBE-FCI) formalism [J. Phys. Chem. Lett., 8, 4633 (2017)], will be presented. Specifically, in this opening part of the series, the capabilities of the MBE-FCI method in producing near-exact ground state energies for weakly correlated molecules of any spin multiplicity will be demonstrated., Comment: 38 pages, 7 tables, 3 figures, 1 SI attached as an ancillary file

In the present letter, it is demonstrated how full configuration interaction (FCI) results in extended basis sets may be obtained to within sub-kJ/mol accuracy by decomposing the energy in terms of many-body expansions in the virtual orbitals of the molecular system at hand. This extension of the FCI application range lends itself to two unique features of the current approach, namely that the total energy calculation can be performed entirely within considerably reduced orbital subspaces and may be so by means of embarrassingly parallel programming. Facilitated by a rigorous and methodical screening protocol and further aided by expansion points different from the Hartree-Fock solution, all-electron numerical results are reported for H$_2$O in polarized core-valence basis sets ranging from double-$\zeta$ (10 $e$, 28 $o$) to quadruple-$\zeta$ (10 $e$, 144 $o$) quality., Comment: 20 pages, 3 figures, 1 table. * With respect to the original arXiv version (v1), the present version of the letter contains updated results. The original TZ and QZ values were unfortunately in error due to a subtle PySCF bug, which has since then been fixed

We extend our assessment of the potential of perturbative coupled cluster (CC) expansions for a test set of open-shell atoms and organic radicals to the description of quadruple excitations. Namely, the second- through sixth-order models of the recently proposed CCSDT(Q-n) quadruples series [J. Chem. Phys. 140, 064108 (2014)] are compared to the prominent CCSDT(Q) and lambda-CCSDT(Q) models. From a comparison of the models in terms of their recovery of total CC singles, doubles, triples, and quadruples (CCSDTQ) energies, we find that the performance of the CCSDT(Q-n) models is independent of the reference used (unrestricted or restricted (open-shell) Hartree-Fock), in contrast to the CCSDT(Q) and lambda-CCSDT(Q) models, for which the accuracy is strongly dependent on the spin of the molecular ground state. By further comparing the ability of the models to recover relative CCSDTQ atomization energies, the discrepancy between them is found to be even more pronounced, stressing how a balanced description of both closed- and open-shell species - as found in the CCSDT(Q-n) models - is indeed of paramount importance if any perturbative CC model is to be of chemical relevance for high-accuracy applications. In particular, the third-order CCSDT(Q-3) model is found to offer an encouraging alternative to the existing choices of quadruples models used in modern computational thermochemistry, since the model is still only of moderate cost, albeit markedly more costly than, e.g., the CCSDT(Q) and lambda-CCSDT(Q) models., Comment: 18 pages, 3 figures, 1 supporting information (attached as an ancillary file)

The accuracy at which total energies of open-shell atoms and organic radicals may be calculated is assessed for selected coupled cluster perturbative triples expansions, all of which augment the coupled cluster singles and doubles (CCSD) energy by a non-iterative correction for the effect of triple excitations. Namely, the second- through sixth-order models of the recently proposed CCSD(T-n) triples series [J. Chem. Phys. 140, 064108 (2014)] are compared to the acclaimed CCSD(T) model for both unrestricted as well as restricted open-shell Hartree-Fock (UHF/ROHF) reference determinants. By comparing UHF- and ROHF-based statistical results for a test set of 18 modest-sized open-shell species with comparable RHF-based results, no behavioral differences are observed for the higher-order models of the CCSD(T-n) series in their correlated descriptions of closed- and open-shell species. In particular, we find that the convergence rate throughout the series towards the coupled cluster singles, doubles, and triples (CCSDT) solution is identical for the two cases. For the CCSD(T) model, on the other hand, not only its numerical consistency, but also its established, yet fortuitous cancellation of errors breaks down in the transition from closed- to open-shell systems. The higher-order CCSD(T-n) models (orders n>3) thus offer a consistent and significant improvement in accuracy relative to CCSDT over the CCSD(T) model, equally for RHF, UHF, and ROHF reference determinants, albeit at an increased computational cost., Comment: 22 pages, 5 figures, 1 table, 1 supporting information (attached as an ancillary file)

An implementation of coupled-cluster (CC) theory to treat atoms and molecules in finite magnetic fields is presented. The main challenges stem from the magnetic-field dependence in the Hamiltonian, or, more precisely, the appearance of the angular momentum operator, due to which the wave function becomes complex and which introduces a gauge-origin dependence. For this reason, an implementation of a complex CC code is required together with the use of gauge-including atomic orbitals to ensure gauge-origin independence. Results of coupled-cluster singles--doubles--perturbative-triples (CCSD(T)) calculations are presented for atoms and molecules with a focus on the dependence of correlation and binding energies on the magnetic field.

We have observed the rotational ground-state (J = 1-0) transitions of DCO+, HN13C and DNC with the IRAM 30m telescope toward the dark cloud LDN 1512 which has exceptionally narrow lines permitting hyperfine splitting to be resolved in part. The measured splittings of 50-300 kHz are used to derive nuclear quadrupole and spin-rotation parameters for these species. The measurements are supplemented by high-level quantum-chemical calculations using coupled-cluster techniques and large atomic-orbital basis sets. We find eQq = +151.12 (400) kHz and C_I = -1.12 (43) kHz for DCO+, eQq = 272.5 (51) kHz for HN13C, and eQq(D) = 265.9 (83) kHz and eQq(N) = 288.2 (71) kHz for DNC. The numbers for DNC are consistent with previous laboratory data, while our constants for DCO+ are somewhat smaller than previous results based on astronomical data. For both DCO+ and DNC, our results are more accurate than previous determinations. Our results are in good agreement with the corresponding best theoretical estimates. We also derive updated rotational constants for HN13C: B = 43545.6000 (47) MHz and D = 93.7 (20) kHz. The hyperfine splittings of the DCO+, DNC and HN13C J = 1-0 lines range over 0.47-1.28 km/s, which is comparable to typical line widths in pre-stellar cores and to systematic gas motions on ~1000 AU scales in protostellar cores. We present tabular information to allow inclusion of the hyperfine splitting in astronomical data interpretation. The large differences in the 14N quadrupole parameters of DNC and HN13C have been traced to differences in the vibrational corrections caused by significant non-rigidity of these molecules, particularly along the bending coordinate., Comment: Accepted by A&A; 11 pages, 3 figures

We are proposing a new computational thermochemistry protocol denoted W3 theory, as a successor to W1 and W2 theory proposed earlier [Martin and De Oliveira, J. Chem. Phys. 111, 1843 (1999)]. The new method is both more accurate overall (error statistics for total atomization energies approximately cut in half) and more robust (particularly towards systems exhibiting significant nondynamical correlation) than W2 theory. The cardinal improvement rests in an approximate account for post-CCSD(T) correlation effects. Iterative T_3 (connected triple excitations) effects exhibit a basis set convergence behavior similar to the T_3 contribution overall. They almost universally decrease molecular binding energies. Their inclusion in isolation yields less accurate results than CCSD(T) nearly across the board: it is only when T_4 (connected quadruple excitations) effects are included that superior performance is achieved. $T_4$ effects systematically increase molecular binding energies. Their basis set convergence is quite rapid, and even CCSDTQ/cc-pVDZ scaled by an empirical factor of 1.2532 will yield a quite passable quadruples contribution. The effect of still higher-order excitations was gauged for a subset of molecules (notably the eight-valence electron systems): T_5 (connected quintuple excitations) contributions reach 0.3 kcal/mol for the pathologically multireference X ^1\Sigma^+_g state of C_2 but are quite small for other systems. A variety of avenues for achieving accuracy beyond that of W3 theory were explored, to no significant avail. W3 thus appears to represent a good compromise between accuracy and computational cost for those seeking a robust method for computational thermochemistry in the kJ/mol accuracy range on small systems., Comment: J. Chem. Phys., in press (306406JCP)