Search

An efficient implementation of the quadratic unitary coupled-cluster singles and doubles (qUCCSD) scheme for calculations of electronic ground and excited states using an unrestricted molecular spin-orbital formulation and an efficient tensor contraction library is reported. The accuracy of the qUCCSD scheme and the efficiency of the present implementation are demonstrated using extensive benchmark calculations of excitation energies and an application to S

Dual fluorescence in 4-(dimethylamino)benzonitrile (DMABN) and its derivatives in polar solvents has been studied extensively for the past several decades. An intramolecular charge transfer (ICT) minimum on the excited state potential energy surface, in addition to the localized low-energy (LE) minimum, has been proposed as a mechanism for this dual fluorescence, with large geometric relaxation and molecular orbital reorganization a key feature of the ICT pathway. Herein, we have used both equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) and time-dependent density functional (TDDFT) methods to investigate the landscape of excited state potential energy surfaces across a number of geometric conformations proposed as ICT structures. In order to correlate these geometries and valence excited states in terms of potential experimental observables, we have calculated the nitrogen K-edge ground and excited state absorption spectra for each of the predicted "signpost" structures, and identified several key spectral features which could be used to interpret a future time-resolved x-ray absorption experiment.

The phenomenon of orbital relaxation upon excitation of core electrons is a major problem in the linear-response treatment of core-hole spectroscopies. Rather than addressing relaxation through direct dynamical correlation of the excited state via equation-of-motion coupled cluster theory (EOMEE-CC), we extend the alternative similarity-transformed equation-of-motion coupled cluster theory (STEOMEE-CC) by including the core-valence separation (CVS) and correlation of triple excitations only within the calculation of core ionization energies. This new method, CVS-STEOMEE-CCSD+cT, significantly improves on CVS-EOMEE-CCSD and unmodified CVS-STEOMEE-CCSD when compared to full CVS-EOM-CCSDT for K-edge core-excitation energies of a set of small molecules. The improvement in both absolute and relative (shifted) peak positions is nearly as good as that for transition-potential coupled cluster (TP-CC), which includes an explicit treatment of orbital relaxation, and CVS-EOMEE-CCSD*, which includes a perturbative treatment of triple excitations.

In the framework of the computational determination of highly-accurate vertical excitation energies in small organic compounds, we explore the possibilities offered by the equation-of-motion formalism relying on the approximate fourth-order coupled-cluster (CC) method, CC4. We demonstrate, using an extended set of more than 200 reference values based on CC including up to quadruples excitations (CCSDTQ), that CC4 is an excellent approximation to CCSDTQ for excited states with a dominant contribution from single excitations with an average deviation as small as $0.003$ eV. We next assess the accuracy of several additive basis set correction schemes, in which vertical excitation energies obtained with a compact basis set and a high-order CC method are corrected with lower-order CC calculations performed in larger basis sets. Such strategies are found to be overall very beneficial, though their accuracy depend significantly on the actual scheme. Finally, CC4 is employed to improve several theoretical best estimates of the QUEST database for molecules containing between four and six (non-hydrogen) atoms, for which previous estimates were computed at the CCSDT level., 12 pages, 1 figure (supporting information available)

Community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.

The overriding theme of John Stanton's prestigious career has been the close collaboration between experiment and theory. This theme began early in John's academic training as a Ph.D. student in th...

In addition to tensor contractions, one of the most pronounced computational bottlenecks in the nonorthogonally spin-adapted forms of the quantum chemistry methods CCSDT and CCSDTQ, and their approximate forms—including CCSD(T) and CCSDT(Q)—are spin summations. At a first sight, spin summations are operations similar to tensor transpositions, but a closer look reveals additional challenges to high-performance calculations, including temporal locality and scattered memory accesses. This article explores a sequence of algorithmic solutions for spin summations, each exploiting individual properties of either the underlying hardware (e.g., caches, vectorization) or the problem itself (e.g., factorizability). The final algorithm combines the advantages of all the solutions while avoiding their drawbacks; this algorithm achieves high performance through parallelization and vectorization, and by exploiting the temporal locality inherent to spin summations. Combined, these optimizations result in speedups between 2.4× and 5.5× over the NCC quantum chemistry software package. In addition to such a performance boost, our algorithm can perform the spin summations in-place , thus reducing the memory footprint by 2× over an out-of-place variant.

Substituent effects on the carbon���carbon bonds of hydrocarbons have been a topic of interest within the past seven decades as resultant information would enable one to tune the activity of CC bonds. However, current assessments of the C���C, C=C, and C���C bond strength of acetylene, ethylene, and ethane as well as their derivatives rely on indirect measures such as bond length and bond dissociation enthalpy. In this work, we introduce a quantitative measure of the intrinsic strength of C���C, C=C, and C���C bonds for a set of 40 hydrocarbon systems consisting of 3 parent structures, 36 hydrocarbon derivatives involving CF3, CH3, CHO, F, NH2, or OH groups, and a conjugated system, based on vibrational spectroscopy. Local mode force constants ka(CC) were computed at the CCSD(T)/cc-pVTZ level of theory for systems 1���32 and 34���40 and CCSD(T)/cc-pVDZ for 33. From ka(CC), we derived related bond strength orders BSO n(CC) in order to provide quantitative measures of intrinsic bond strength. Topological electron density and natural population analyses were carried out as to analyze the nature of these bonds and complement bond strength measures. For substituted hydrocarbon systems, we found the strengthening/weakening of the CC bonds occurs as the covalent nature of the bond increases/decreases by means of varying charge delocalizations. Our findings provide new guidelines for desirably modulating C���C, C=C, and C���C bond strength and for the design of prospective pathways for bond cleavage reactions.

The analytic gradient theory for both iterative and non-iterative coupled-cluster approximations that include connected quadruple excitations is presented. These methods include, in particular, CCSDT(Q), which is an analog of the well-known CCSD(T) method which starts from the full CCSDT method rather than CCSD. The resulting methods are implemented in the CFOUR program suite, and pilot applications are presented for the equilibrium geometries and harmonic vibrational frequencies of the simplest Criegee intermediate, CH$_2$OO, as well as to the isomerization pathway between dimethylcarbene and propene. While all methods are seen to approximate the full CCSDTQ results well for "well-behaved" systems, the more difficult case of the Criegee intermediate shows that CCSDT(Q), as well as certain iterative approximations, display problematic behavior., 37 pages, 1 figure, 2 tables

The concept of the direct product decomposition (DPD) is extended to arbitrary tensors while maintaining the same theoretical reduction in storage and computation. Additionally, the structure of th...

An up-to-date overview of the CFOUR program system is given. After providing a brief outline of the evolution of the program since its inception in 1989, a comprehensive presentation is given of its well-known capabilities for high-level coupled-cluster theory and its application to molecular properties. Subsequent to this generally well-known background information, much of the remaining content focuses on lesser-known capabilities of CFOUR, most of which have become available to the public only recently or will become available in the near future. Each of these new features is illustrated by a representative example, with additional discussion targeted to educating users as to classes of applications that are now enabled by these capabilities. Finally, some speculation about future directions is given, and the mode of distribution and support for CFOUR are outlined. ispartof: JOURNAL OF CHEMICAL PHYSICS vol:152 issue:21 ispartof: location:United States status: published

A number of iterative and perturbative approximations to the full equation-of-motion coupled cluster method with single, double, and triple excitations (EOM-CCSDT) are evaluated in the context of calculating the K-edge core-excitation and core-ionisation energies of several small molecules. Several of these methods are found to accurately reproduce the full EOM-CCSDT energies well, in particular the EOM-CCSD* method which, when using the core-valence separation (CVS) ansatz, scales rigorously with the sixth power of molecular size. These results highlight that very accurate core excitation and ionisation energies for molecules with first-row atoms can be computed at a cost not much larger than that for EOM-CCSD.

We report the first investigation of the performance of EOM-CC4 -- an approximate equation-of-motion coupled-cluster model which includes iterative quadruple excitations -- for vertical excitation energies in molecular systems. By considering a set of 28 excited states in 10 small molecules for which we have computed CCSDTQP and FCI reference energies, we show that, in the case of excited states with a dominant contribution from the single excitations, CC4 yields excitation energies with sub-kJ~mol$^{-1}$ accuracy (i.e., error below $0.01$ eV), in very close agreement with its more expensive CCSDTQ parent. Therefore, if one aims at high accuracy, CC4 stands as a highly competitive approximate method to model molecular excited states, with a significant improvement over both CC3 and CCSDT. Our results also evidence that, although the same qualitative conclusions hold, one cannot reach the same level of accuracy for transitions with a dominant contribution from the double excitations., 6 pages, 2 figures

The problem of orbital relaxation in computational core-hole spectroscopies, including x-ray absorption and x-ray photoionization, has long plagued linear response approaches, including equation-of-motion coupled cluster with singles and doubles (EOM-CCSD). Instead of addressing this problem by including additional electron correlation, we propose an explicit treatment of orbital relaxation via the use of "transition potential" reference orbitals, leading to a transition-potential coupled cluster (TP-CC) family of methods. One member of this family in particular, TP-CCSD(1/2), is found to essentially eliminate the orbital relaxation error and achieve the same level of accuracy for core-hole spectra as is typically expected of EOM-CCSD in the valence region. These results show that very accurate x-ray absorption spectra for molecules with first-row atoms can be computed at a cost essentially the same as that for EOM-CCSD., 25 pages, 5 figures, 1 table

The method of sub-iteration, which was previously applied to the higher-order coupled cluster amplitude equations, is extended to the case of the coupled cluster $\Lambda$ equations. The sub-iteration procedure for the $\Lambda$ equations is found to be highly similar to that for the amplitude equations, and to exhibit a similar improvement in rate of convergence relative to extrapolation of all $\hat{T}$ or $\hat{\Lambda}$ amplitudes using DIIS. A method of dynamic damping is also presented which is found to effectively recover rapid convergence in the case of oscillatory behavior in the amplitude or $\Lambda$ equations. Together, these techniques allow for the convergence of both the amplitude and $\Lambda$ equations necessary for the calculation of analytic gradients and properties of higher-order coupled cluster methods without the high memory or disk I/O cost of full DIIS extrapolation., Comment: 18 pages, 8 figures

Benchmark scalar-relativistic core-valence-separated (CVS) equation-of-motion coupled-cluster ionization potential (EOMIP-CC) calculations of 21 K-edge ionization energies of C, O, N, and F in 14 molecules are reported. The CVS-EOMIP-CC methods are shown to be numerically more stable and more accurate than the parent EOMIP-CC methods, even when the calculations using the latter can be tightly converged. The superior performance of the CVS scheme is attributed to the exclusion of spurious couplings between core-ionized states and valence continuum states. Systematic improvement of computed K-edge ionization energies within the CVS-EOMIP-CC hierarchy, including the CC singles and doubles (CCSD) method, the CC singles, doubles, and triples (CCSDT) method, and the CC singles, doubles, triples, and quadruples (CCSDTQ) method, is demonstrated, with CCSDTQ yielding essentially quantitative results. Maximum absolute deviations between computed and experimental results amount to 2.54 eV for CCSD/cc-pCVQZ, 0.54 eV for CCSDT/cc-pCVQZ, and 0.23 eV for CCSDT/cc-pCVQZ augmented with quadruples contributions using the cc-pVTZ basis sets. The corresponding standard deviations are 1.91 eV for CCSD/cc-pCVQZ, 0.18 eV for CCSDT/cc-pCVQZ, and 0.10 eV for CCSDT/cc-pCVQZ augmented with quadruples contributions using the cc-pVTZ basis sets. Finally, CVS-EOMIP-CCSDT/cc-pCVTZ calculations of core ionization energies in CH3CN and CH3NC are reported, and experimental re-investigation of carbon 1s ionization energies in CH3CN is suggested.

An extensive analysis has been carried out of the performance of standard families of basis sets with the hierarchy of coupled cluster methods CC2, CCSD, CC3, and CCSDT in computing selected Oxygen, Carbon, and Nitrogen K-edge (vertical) core excitation and ionization energies within a core-valence separated scheme in the molecules water, ammonia, and carbon monoxide. Complete basis set limits for the excitation energies have been estimated via different basis set extrapolation schemes. The importance of scalar relativistic effects has been established within the spin-free exact two-component theory in its one-electron variant (SFX2C-1e).

Goldstone diagrams, interpreted in terms of orbital (rather than spin-orbital) integrals and amplitudes, provide an efficient representation of many-body theories for closed-shell reference functions. Such diagrams are well-known in the context of many-body perturbation theory, but their usefulness extends to theories such as coupled cluster as well. However, the Goldstone diagrams are more numerous than their Brandow or Hugenholtz counterparts due to a lack of antisymmetry in the orbital quantities, and the resultant equations can always be factorized into a more efficient form. Here we present a new diagrammatic interpretation for closed-shell orbital theories that starts directly from the original Brandow diagrams. The resulting equations are in terms of orbital quantities, but are based on “spin-summed” amplitudes and integrals, allowing them to be cast into an efficient and factorized form from the start. The derivation of the new diagrammatic formulation proceeds directly from the usual spin-orbital diagrammatic rules. A discussion of the differences between the spin-orbital and traditional Goldstone forms is given, along with several illustrative examples of the new diagrammatics.

A new method for generating fitting grids for least-squares tensor hypercontraction (LS-THC) is presented. This method draws inspiration from the related interpolative separable density fitting (ISDF) technique, but uses only a pivoted Cholesky decomposition of the metric matrix, S, already computed as a matter of course in LS-THC. The size and quality of the resulting grid is controlled by a user-defined cutoff parameter and the size of the starting grid. Additionally, the Cholesky-based method provides an alternative and possible more numerically stable method for performing the least-squares fit. The quality of the grids produced is evaluated for LS-DF-THC-MP2 calculations on retinal and benzene, the former with a large starting grid and small cc-pVDZ basis set, and the latter with a wide range of grids and basis sets. The error and grid size is found to be well-controlled by either the cutoff parameter (with a large starting grid) or the starting grid size (with a tight cutoff) and highly predictable. The Cholesky-based method is also able to generate unique grids tailored to different charge distributions, for example the (ab|, (ai|, and (ij| distributions that arise in the molecular orbital integrals. While only the (ai| grid directly affects the MP2 energy, the relative sizes of the other grids are examined., Comment: 11 pages, 1 figure, supporting information

The validation of the quality of the description of excited electronic states is of special importance in quantum chemistry as the general reliability of ab initio methods shows a much larger variation for these states than for the ground state. In this study, we investigate the quality of excited state energy gradients and potential energy surfaces on selected systems, as provided by the single reference coupled cluster variants CC2, CCSD, CCSD(T)(a)*, and CC3. Gradients and surface plots that follow the Franck–Condon forces are compared to the respective CCSDT reference values, thereby establishing a useful strategy for judging each variant’s accuracy. The results reveal serious flaws of lower order methods - in particular, CC2 - in several situations where they otherwise give accurate vertical excitation energies, as well as excellent accuracy and consistency of the recently proposed CCSD(T)(a)* method.

A derivation of fourth-order vibrational perturbation theory (VPT4) based on the Watson Hamiltonian in dimensionless rectilinear normal coordinates is presented. Terms that are linear and cubic in the (nk + 12), with nk being the zeroth-order harmonic oscillator quantum numbers, appear in fourth order and extend the much simpler second-order vibrational perturbation theory model. The rather involved expressions for the fourth-order terms are derived with Rayleigh-Schrodinger perturbation theory, the process of verifying their correctness is described, and a computer code to generate the VPT4 constants from the potential energy surface derivatives is provided. The paper concludes with numerical examples featuring the H2O, Si2C, and cyclic-C3H2 molecules.A derivation of fourth-order vibrational perturbation theory (VPT4) based on the Watson Hamiltonian in dimensionless rectilinear normal coordinates is presented. Terms that are linear and cubic in the (nk + 12), with nk being the zeroth-order harmonic oscillator quantum numbers, appear in fourth order and extend the much simpler second-order vibrational perturbation theory model. The rather involved expressions for the fourth-order terms are derived with Rayleigh-Schrodinger perturbation theory, the process of verifying their correctness is described, and a computer code to generate the VPT4 constants from the potential energy surface derivatives is provided. The paper concludes with numerical examples featuring the H2O, Si2C, and cyclic-C3H2 molecules.

The fraction of the collisionally stabilized Criegee species CH2OO produced from the ozonolysis of ethylene is calculated using a two-dimensional (E, J)-grained master equation technique and semiclassical transition-state theory based on the potential energy surface obtained from high-accuracy quantum chemical calculations. Our calculated yield of 42 ± 6% for the stabilized CH2OO agrees well, within experimental error, with available (indirect) experimental results. Inclusion of angular momentum in the master equation is found to play an essential role in bringing the theoretical results into agreement with the experiment. Additionally, yields of HO and HO2 radical products are predicted to be 13 ± 6% and 17 ± 6%, respectively. In the kinetic simulation, the HO radical product is produced mostly from the stepwise decomposition mechanism of primary ozonide rather than from dissociation of hot CH2OO.

This document is the Accepted Manuscript version of a Published Work that appeared in final form in Journal of Chemical Theory and Computation, copyright © American Chemical Society after peer review and technical editing by the publisher. To access the final edited and published work see https://doi.org/10.1021/acs.jctc.7b01138. Geometric energy derivatives which rely on core-corrected focal-point energies extrapolated to the complete basis set (CBS) limit of coupled cluster theory with iterative and noniterative quadruple excitations, CCSDTQ and CCSDT(Q), are used as elements of molecular gradients and, in the case of CCSDT(Q), expansion coefficients of an anharmonic force field. These gradients are used to determine the CCSDTQ/CBS and CCSDT(Q)/CBS equilibrium structure of the S0 ground state of H2CO where excellent agreement is observed with previous work and experimentally derived results. A fourth-order expansion about this CCSDT(Q)/CBS reference geometry using the same level of theory produces an exceptional level of agreement to spectroscopically observed vibrational band origins with a MAE of 0.57 cm–1. Second-order vibrational perturbation theory (VPT2) and variational discrete variable representation (DVR) results are contrasted and discussed. Vibration–rotation, anharmonicity, and centrifugal distortion constants from the VPT2 analysis are reported and compared to previous work. Additionally, an initial application of a sum-over-states fourth-order vibrational perturbation theory (VPT4) formalism is employed herein, utilizing quintic and sextic derivatives obtained with a recursive algorithmic approach for response theory.

Matrix multiplication (GEMM) is a core operation to numerous scientific applications. Traditional implementations of Strassen-like fast matrix multiplication (FMM) algorithms often do not perform well except for very large matrix sizes, due to the increased cost of memory movement, which is particularly noticeable for non-square matrices. Such implementations also require considerable workspace and modifications to the standard BLAS interface. We propose a code generator framework to automatically implement a large family of FMM algorithms suitable for multiplications of arbitrary matrix sizes and shapes. By representing FMM with a triple of matrices [U, V, W] that capture the linear combinations of submatrices that are formed, we can use the Kronecker product to define a multi-level representation of Strassen-like algorithms. Incorporating the matrix additions that must be performed for Strassen-like algorithms into the inherent packing and micro-kernel operations inside GEMM avoids extra workspace and reduces the cost of memory movement. Adopting the same loop structures as high-performance GEMM implementations allows parallelization of all FMM algorithms with simple but efficient data parallelism without the overhead of task parallelism. We present a simple performance model for general FMM algorithms and compare actual performance of 20+ FMM algorithms to modeled predictions. Our implementations demonstrate a performance benefit over conventional GEMM on single core and multi-core systems. This study shows that Strassen-like fast matrix multiplication can be incorporated into libraries for practical use.

Precise calculation of molecular electronic wavefunctions by methods such as coupled-cluster requires the computation of tensor contractions, the cost of which has polynomial computational scaling with respect to the system and basis set sizes. Each contraction may be executed via matrix multiplication on a properly ordered and structured tensor. However, data transpositions are often needed to reorder the tensors for each contraction. Writing and optimizing distributed-memory kernels for each transposition and contraction is tedious since the number of contractions scales combinatorially with the number of tensor indices. We present a distributed-memory numerical library (Cyclops Tensor Framework (CTF)) that automatically manages tensor blocking and redistribution to perform any user-specified contractions. CTF serves as the distributed-memory contraction engine in Aquarius, a new program designed for high-accuracy and massively-parallel quantum chemical computations. Aquarius implements a range of coupled-cluster and related methods such as CCSD and CCSDT by writing the equations on top of a C++ templated domain-specific language. This DSL calls CTF directly to manage the data and perform the contractions. Our CCSD and CCSDT implementations achieve high parallel scalability on the BlueGene/Q and Cray XC30 supercomputer architectures showing that accurate electronic structure calculations can be effectively carried out on top of general distributed-memory tensor primitives. We introduce Cyclops Tensor Framework (CTF), a distributed-memory library for tensor contractions.CTF is able to perform tensor decomposition, redistribution, and contraction at runtime.CTF enables the expression of massively-parallel coupled-cluster methods via a concise tensor contraction interface.The quantum chemistry software suite Aquarius employs CTF to execute two coupled-cluster methods: CCSD and CCSDT.The Aquarius CCSD and CCSDT codes scale well on BlueGene/Q and Cray XC30, comparing favorably to NWChem.

Working out how to pack benzene in silico Many organic compounds crystallize in several different energetically similar packing arrangements, or polymorphs. This complicates processes such as drug formulation that rely on reproducible crystallization. Yang et al. have now achieved the long-standing goal of calculating a crystal packing arrangement from first principles to an accuracy that can distinguish polymorphs (see the Perspective by Price). They used benzene as a prototypical test case and applied quantum chemical methods that improve estimates of multibody interactions. The results bode well for future applications of theory to optimization of crystallization protocols. Science , this issue p. 640 ; see also p. 619

Tensor contraction (TC) is an important computational kernel widely used in numerous applications. It is a multi-dimensional generalization of matrix multiplication (GEMM). While Strassen's algorithm for GEMM is well studied in theory and practice, extending it to accelerate TC has not been previously pursued. Thus, we believe this to be the first paper to demonstrate how one can in practice speed up tensor contraction with Strassen's algorithm. By adopting a Block-Scatter-Matrix format, a novel matrix-centric tensor layout, we can conceptually view TC as GEMM for a general stride storage, with an implicit tensor-to-matrix transformation. This insight enables us to tailor a recent state-of-the-art implementation of Strassen's algorithm to TC, avoiding explicit transpositions (permutations) and extra workspace, and reducing the overhead of memory movement that is incurred. Performance benefits are demonstrated with a performance model as well as in practice on modern single core, multicore, and distributed memory parallel architectures, achieving up to 1.3x speedup. The resulting implementations can serve as a drop-in replacement for various applications with significant speedup.

Accurate methods for excited, ionized, and electron-attached states are critical to the study of many chemical species such as reactive intermediates, radicals, and ionized systems. The equation-of-motion coupled cluster singles, doubles, and triples (EOM-CCSDT) family of methods is very accurate (roughly similar in accuracy as for CCSDT calculations of the ground state), but the computational cost scales iteratively as the eighth power of the system size. Many approximations already exist, although most either correct only the excited state or require an iterative 𝒪(n7) procedure which can also be prohibitively expensive. In this paper, new methods, termed EOM-CCSD(T)(a) and EOM-CCSD(T)(a)*, are proposed which correct both the ground and excited states based on a shared effective Hamiltonian, and the latter of which includes only non-iterative corrections to both the CCSD and EOM-CCSD energies. These methods are found to significantly improve the description of excited and ionized potential energy surfac...

We use, for the first time, ab initio coupled-cluster theory to compute the spectral function of the uniform electron gas at a Wigner-Seitz radius of $r_\mathrm{s}=4$. The coupled-cluster approximations we employ go significantly beyond the diagrammatic content of state-of-the-art $GW$ theory. We compare our calculations extensively to $GW$ and $GW$-plus-cumulant theory, illustrating the strengths and weaknesses of these methods in capturing the quasiparticle and satellite features of the electron gas. Our accurate calculations further allow us to address the long-standing debate over the occupied bandwidth of metallic sodium. Our findings indicate that the future application of coupled-cluster theory to condensed phase material spectra is highly promising., Comment: 6 pages, 2 figures

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., 22 pages, 5 figures, 1 table, 1 supporting information (attached as an ancillary file)

The infrared spectrum of methyl nitrate has been reinvestigated. Analysis of the 0.125 cm −1 room temperature spectrum has been carried out with the assistance of high-level ab initio calculations using coupled-cluster theory. Use of second-order vibrational perturbation theory (VPT2), with allowances—when necessary—for resonances has permitted the complete assignment of all fundamental vibrational modes in this molecule. A few errors in previous assignments are corrected (notably the CH stretching modes that are now assigned as ν 1 = 3034.1 cm −1 ; ν 2 = 2961.5 cm −1 , and ν 13 = 3014.5 cm −1 and the methyl deformation ν 15 , assigned to 1155.7 cm −1 ), and the work also presents a highly accurate empirically adjusted force field that reproduces all fundamental frequencies to within 0.1 cm −1 . This Hamiltonian should be accurate enough to predict two- and three-quantum transitions to within a few cm −1 . The remaining fundamental frequencies (in cm −1 ) are ν 4 = 1667.0; ν 5 = 1431.8; ν 6 = 1291.2; ν 7 = 1177.9; ν 8 = 1018.5; ν 9 = 855.4; ν 10 = 658.3; ν 11 = 571.0; ν 12 = 344.2; ν 14 = 1444.7; ν 16 = 760.6; ν 17 = 204.5; and ν 18 = 133.7.

We consider two distinct coupled cluster (CC) perturbation series that both expand the difference between the energies of the CCSD (CC with single and double excitations) and CCSDT (CC with single, double, and triple excitations) models in orders of the M{\o}ller-Plesset fluctuation potential. We initially introduce the E-CCSD(T-$n$) series, in which the CCSD amplitude equations are satisfied at the expansion point, and compare it to the recently developed CCSD(T-$n$) series [J. Chem. Phys. 140, 064108 (2014)], in which not only the CCSD amplitude, but also the CCSD multiplier equations are satisfied at the expansion point. The computational scaling is similar for the two series, and both are term-wise size extensive with a formal convergence towards the CCSDT target energy. However, the two series are different, and the CCSD(T-$n$) series is found to exhibit a more rapid convergence up through the series, which we trace back to the fact that more information at the expansion point is utilized than for the E-CCSD(T-$n$) series. The present analysis can be generalized to any perturbation expansion representing the difference between a parent CC model and a higher-level target CC model. In general, we demonstrate that, whenever the parent parameters depend upon the perturbation operator, a perturbation expansion of the CC energy (where only parent amplitudes are used) differs from a perturbation expansion of the CC Lagrangian (where both parent amplitudes and parent multipliers are used). For the latter case, the bivariational Lagrangian formulation becomes more than a convenient mathematical tool, since it facilitates a different and faster convergent perturbation series than the simpler energy-based expansion., Comment: 26 pages, 1 figure, 1 supporting information (attached as an ancillary file)

Tensor computations--in particular tensor contraction (TC)--are important kernels in many scientific computing applications. Due to the fundamental similarity of TC to matrix multiplication (MM) and to the availability of optimized implementations such as the BLAS, tensor operations have traditionally been implemented in terms of BLAS operations, incurring both a performance and a storage overhead. Instead, we implement TC using the flexible BLIS framework, which allows for transposition (reshaping) of the tensor to be fused with internal partitioning and packing operations, requiring no explicit transposition operations or additional workspace. This implementation, TBLIS, achieves performance approaching that of MM, and in some cases considerably higher than that of traditional TC. Our implementation supports multithreading using an approach identical to that used for MM in BLIS, with similar performance characteristics. The complexity of managing tensor-to-matrix transformations is also handled automatically in our approach, greatly simplifying its use in scientific applications., Comment: 24 pages, 8 figures, uses pgfplots

The coupled-cluster singles and doubles treatment with a perturbative treatment of triple excitations known as CCSD(T) has been used in conjunction with a hierarchy of atomic natural orbital basis sets to study stretching levels in water up to 2 eV above the zero-point level. Agreement with experiment obtained with perturbation theory augmented by a simple treatment of resonances suggested by Lehmann is quite remarkable in spite of the well-known and strong Darling–Dennison resonances in this system. With the largest basis sets, deviation from experiment is ca. 20 cm−1 at energies in the 15,000–16,000 cm−1 range.

The benefits of what is alternatively called a nonorthogonally spin-adapted, spin-free, or orbital representation of the coupled cluster equations is discussed relative to orthogonally spin-adapted, spin-orbital, and spin-integrated theories. In particular, specific linear combinations of the orbital cluster amplitudes, denoted spin-summed amplitudes, are shown to reduce the number of contractions that must be explicitly performed and to simplify the expressions and their derivation. The computational efficiency of the spin-summed approach is discussed and compared to orthogonally spin-adapted and spin-integrated approaches. The spin-summed approach is shown to have significant computational advantages relative to the other methods for CCSDT and CCSDTQ.

We compare the numerical performance of various non-iterative coupled cluster (CC) quadruples models. The results collectively show how approaches that attempt to correct the CC singles and doubles energy for the combined effect of triple and quadruple excitations all fail at recovering the correlation energy of the full CC singles, doubles, triples, and quadruples (CCSDTQ) model to within sufficient accuracy. Such a level of accuracy is only achieved by models that make corrections to the full CC singles, doubles, and triples (CCSDT) energy for the isolated effect of quadruple excitations of which the CCSDT(Q–3) and CCSDT(Q–4) models of the Lagrangian-based CCSDT(Q–n) perturbation series are found to outperform alternative models that add either of the established [Q] and (Q) corrections to the CCSDT energy.

We report infrared spectra of nitromethane anion, CH3NO2 (-), in the region 700-2150 cm(-1), obtained by Ar predissociation spectroscopy and electron detachment spectroscopy. The data are interpreted in the framework of second-order vibrational perturbation theory based on coupled-cluster electronic structure calculations. The modes in the spectroscopic region studied here are mainly based on vibrations involving the heavier atoms; this work complements earlier studies on nitromethane anion that focused on the CH stretching region of the spectrum. Electron detachment begins at photon energies far below the adiabatic electron affinity due to thermal population of excited vibrational states.