A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and openshell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr 2 dimer, exploring zeolitecatalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.Keywords quantum chemistry, software, electronic structure theory, density functional theory, electron correlation, computational modelling, Q-Chem Disciplines Chemistry CommentsThis article is from Molecular Physics: An International Journal at the Interface Between Chemistry and Physics 113 (2015): 184, doi:10.1080/00268976.2014. RightsWorks produced by employees of the U.S. Government as part of their official duties are not copyrighted within the U.S. The content of this document is not copyrighted. Authors 185A summary of the technical advances that are incorporated in the fourth major release of the Q-CHEM quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly corre...

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design.

The high symmetry and stability of phenalenyl systems, both as the planar pi-radical (P*) and as the pi-cation (P+), are desirable characteristics of prototypical aromatic donor/acceptor pairs that encourage their use as (binary) models for the study of intermolecular interactions extant in stacked molecular arrays. Thus, quantitative ESR spectroscopy of the paramagnetic P* identifies its spontaneous self-association to the diamagnetic P2, previously characterized as the stacked pi-dimer by X-ray crystallography. Likewise, the rapid cross-association of P* with the closed-shell P+ leads to the stacked pi-dimer cation P2*+ with the "doubled" ESR spectrum diagnostic of complete (odd) electron delocalization. These pi-associations are confirmed by UV-vis studies that reveal diagnostic near-IR bands of both P2 and P2*+-strongly reminiscent of intermolecular charge-transfer absorptions in related aromatic (donor/acceptor) pi-associations. Ab initio molecular-orbital calculations for the pi-dimer P2 predict a binding energy of DeltaED = -11 kcal mol(-1), which is in accord with the experimental enthalpy change of DeltaHD = -9.5 kcal mol(-1) in dichloromethane solution. Most importantly, the calculations reproduce the intermonomer spacings and reveal the delicate interplay of attractive covalent and dispersion forces, balanced against the repulsions between filled orbitals. For comparison, the binding energy in the structurally related cationic pi-pimer P2*+ is calculated to be significantly larger with DeltaEP approximately -20 kcal mol(-1) (gas phase), owing to favorable electrostatic interactions not present in the neutral pi-dimer (which outweigh the partial loss of covalent interactions). As a result, our theoretical formulation can correctly account for the experimental enthalpy change in solution of DeltaHP = -6.5 kcal mol(-1) by the inclusion of differential ionic solvation in the formation of the pi-pimer.

Molecular dynamics simulations are presented for condensed-phase electron transfer (ET) systems where the electronic polarizability of both the solvent and the solute is incorporated. The solute polarizability is allowed to change with electronic transition. The results display notable deviation from the standard free energy parabolas of traditional ET theories. A new three-parameter ET model is applied, and the theory is shown to accurately model the free energy surfaces. This paper presents conclusive evidence that the traditional theory for the free energy barrier of ET reactions requires modification.

Electronic-structure calculations for the self-association of phenalenyl radical (P*) predict the formation of dimeric species (sigma-P2) in which both moieties are connected by a sigma-bond with rP-P approximately 1.59 A and bond dissociation enthalpy of DeltaH(D) approximately 16 kcal mol(-1). Such an unusually weak sigma-bond is related to the loss of aromatic stabilization energy of approximately 34 kcal mol(-1) per phenalenyl moiety, largely owing to rehybridization. Ab initio calculations also reveal that the corresponding (one-electron) bond between phenalenyl radical and its closed-shell cation in sigma-P2+* is unstable relative to dissociation. Time-dependent DFT computations indicate the absence of any (strongly allowed) electronic transition in the visible region of the absorption spectrum of phenalenyl sigma-dimer. Such theoretical predictions are supported by experimental (ESR and UV-NIR) spectroscopic studies, in which the availability of a series of sterically hindered phenalenyl radicals allows definitive separations of the sigma-dimerization process from interference by pi-dimerization. As such, the thermodynamic parameters (determined from the temperature dependence of the ESR signals) with DeltaH(D) = 14 kcal mol(-1) and DeltaS(D) = 52 e.u. can be assigned to the formation of the colorless sigma-dimer. Similar results are obtained for all phenalenyl derivatives (provided their substitution patterns allow sigma-bond formation) to confirm the energetic preference of sigma-dimerization over pi-dimerization.

No abstract

We demonstrate a block-tensor based implementation of coupled-cluster valence-bond singles and doubles (CCVB-SD) [Small, D. W.; Head-Gordon M. J. Chem. Phys. 2012, 137, 114103] which is a simple modification to restricted CCSD (RCCSD) that provides a qualitatively correct description of valence correlations even in strongly correlated systems. We derive the Λ-equation of CCVB-SD and the corresponding unrelaxed density matrices. The resulting production-level implementation is applied to oligoacenes, correlating up to 318 electrons in 318 orbitals. CCVB-SD shows a qualitative agreement with exact methods for short acenes and reaches the bulk limit of oligoacenes in terms of natural orbital occupation numbers, whereas RCCSD shows nonvariational behavior even for relatively short acenes. A significant reduction in polyradicaloid character is found when correlating all valence electrons instead of only the π-electrons.

In this work, we revisited the idea of using the coupled-cluster ground state formalism to target excited states. Our main focus was targeting doubly excited states and double core hole states. Typical equation-of-motion (EOM) approaches for obtaining these states struggle without higher-order excitations than doubles. We showed that by using a non-aufbau determinant optimized via the maximum overlap method the CC ground state solver can target higher energy states. Furthermore, just with singles and doubles (i.e., CCSD), we demonstrated that the accuracy of ∆CCSD and ∆CCSD(T) far surpasses that of EOM-CCSD for doubly excited states. The accuracy of ∆CCSD(T) is nearly exact for doubly excited states considered in this work.For double core hole states, we used an improved ansatz for greater numerical stability by freezing core hole orbitals. The improved methods, core valence separation (CVS)-∆CCSD and CVS-∆CCSD(T), were applied to the calculation of the double ionization potential of small molecules. Even without relativistic corrections, we observed qualitatively accurate results with CVS-∆CCSD and CVS-∆CCSD(T). Remaining challenges in ∆CC include the description of open-shell singlet excited states with the single-reference CC ground state formalism as well as excited states with genuine multi-reference character. The tools and intuition developed in this work may serve as a stepping stone towards directly targeting arbitrary excited states using ground state CC methods.