Modern quantum chemistry with [Open]Molcas

Artículo escrito por un elevado número de autores, sólo se referencian el que aparece en primer lugar, los autores pertenecientes a la UAM y el nombre del grupo de colaboración, si lo hubiereThe following article appeared in The Journal of Chemical Physics 152.21 (2020): 214117 and may be found at https://doi.org/10.1063/5.0004835MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree–Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functionsF.A. acknowledges financial support from the EU-H2020 research and innovation programme under Grant Agreement No. 654360 within the framework of the NFFA-Europe Transnational Access Activity. Part of this work was performed, thanks to computer resources provided by CINECA, under Project No. HPC-EUROPA3 (Grant No. INFRAIA-2016-1-730897), with the support of the EC Research Innovation Action of the H2020 Programme. D.-C.S. and J.A. acknowledge support from the U.S. Department of Energy, Office of Basic Energy Sciences, Heavy Element Chemistry program, under Grant No. DE-SC0001136. S.B. acknowledges support from the Swiss National Science Foundation (Grant No. P2SKP2_184034). A.B. is grateful for support from ETH Zurich (ETH Fellowship No. FEL-49 18-1). M.R. acknowledges support from the Swiss National Science Foundation (Project No. 200021_182400). L.D.V., L.P.-G., and M.Ol. acknowledge a MIUR (Ministero dell’Istruzione, dell’Università e della Ricerca) grant “Dipartimento di Eccellenza 2018-2022.” M.Ol. acknowledges NSF Grant No. CHE-CLP-1710191. M.D. and M.L. acknowledges support from the Olle Engkvist Foundation. E.D.L. and V.V. acknowledge computational resources provided by SNIC through LUNARC and NSC. T.B.P. acknowledges support from the Research Council of Norway through its Centres of Excellence scheme, Project No. 262695, and through Research Grant No. 240698. K.P. acknowledges financial support from KU Leuven through Grant No. C14/15/052. L.S. acknowledges financial support from Ministerio de Economía y Competitividad, Spain (Dirección General de Investigación y Gestión del Plan Nacional de I+D+i, Grant No. MAT2017-83553- P). J.S.-M. acknowledges support from the EU-H2020 Marie Curie Actions (AttoDNA, FP8-MSCA-IF, Grant No. 747662). I.S. gratefully acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (Grant No. 678169 PhotoMutant). L.U. and X.G. gratefully acknowledge scientific Grant Nos. R-143-000-A80- 114 and R-143-000-A65-133 from the National University of Singapore. Computational resources of the NSCC (ASPIRE-1, Grant No. 11001278) were used for this study

Accurate electronic properties and intercalation voltages of olivine-type Li-ion cathode materials from extended Hubbard functionals

The design of novel cathode materials for Li-ion batteries would greatly benefit from accurate first-principles predictions of structural, electronic, and magnetic properties as well as intercalation voltages in compounds containing transition-metal elements. For such systems, density-functional theory (DFT) with standard (semi-)local exchange-correlation functionals is of limited use as it often fails due to strong self-interaction errors that are especially relevant in the partially filled $d$ shells. Here, we perform a detailed comparative study of the phospho-olivine cathode materials Li$_x$MnPO$_4$, Li$_x$FePO$_4$, and the mixed transition metal Li$_x$Mn$_{1/2}$Fe$_{1/2}$PO$_4$ ($x=0, 1/4, 1/2, 3/4, 1$) using four electronic-structure methods: DFT, DFT+$U$, DFT+$U$+$V$, and HSE06. We show that DFT+$U$+$V$, with onsite $U$ and intersite $V$ Hubbard parameters determined from first principles and self-consistently with respect to the structural parameters by means of density-functional perturbation theory (linear response), provides the most accurate description of the electronic structure of these challenging compounds. In particular, we demonstrate that DFT+$U$+$V$ displays very clearly "digital" changes in oxidation states of the transition-metal ions in all compounds, including the mixed-valence phases occurring at intermediate Li concentrations, leading to voltages in remarkable agreement with experiments. We show that the inclusion of intersite Hubbard interactions is essential for the accurate prediction of thermodynamic quantities, balancing the drive for localization induced by the onsite $U$ with intersite $V$ orbital hybridizations

Pulay forces in density-functional theory with extended Hubbard functionals: From nonorthogonalized to orthogonalized manifolds

We present a derivation of the exact expression for Pulay forces in density-functional theory calculations augmented with extended Hubbard functionals, and arising from the use of orthogonalized atomic orbitals as projectors for the Hubbard manifold. The derivative of the inverse square root of the orbital overlap matrix is obtained as a closed-form solution of the associated Lyapunov (Sylvester) equation. The expression for the resulting contribution to the forces is presented in the framework of ultrasoft pseudopotentials and the projector-augmented-wave method, and using a plane wave basis set. We have benchmarked the present implementation with respect to finite differences of total energies for the case of NiO, finding excellent agreement. Owing to the accuracy of Hubbard-corrected density-functional theory calculations - provided the Hubbard parameters are computed for the manifold under consideration - the present work paves the way for systematic studies of solid-state and molecular transition-metal and rare-earth compounds.Comment: 16 pages, 1 figur

Fast noniterative orbital localization for large molecules

We use Cholesky decomposition of the density matrix in atomic orbital basis to define a new set of occupied molecular orbital coefficients. Analysis of the resulting orbitals (“Cholesky molecular orbitals”) demonstrates their localized character inherited from the sparsity of the density matrix. Comparison with the results of traditional iterative localization schemes shows minor differences with respect to a number of suitable measures of locality, particularly the scaling with system size of orbital pair domains used in local correlation methods. The Cholesky procedure for generating orthonormal localized orbitals is noniterative and may be made linear scaling. Although our present implementation scales cubically, the algorithm is significantly faster than any of the conventional localization schemes. In addition, since this approach does not require starting orbitals, it will be useful in local correlation treatments on top of diagonalization-free Hartree-Fock optimization [email protected]

Accurate ab initio density fitting for multiconfigurational self-consistent field methods

Using Cholesky decomposition and density fitting to approximate the electron repulsion integrals, an implementation of the complete active space self-consistent field (CASSCF) method suitable for large-scale applications is presented. Sample calculations on benzene, diaquo-tetra-μ-acetato-dicopper(II), and diuraniumendofullerene demonstrate that the Cholesky and density fitting approximations allow larger basis sets and larger systems to be treated at the CASSCF level of theory with controllable accuracy. While strict error control is an inherent property of the Cholesky approximation, errors arising from the density fitting approach are managed by using a recently proposed class of auxiliary basis sets constructed from Cholesky decomposition of the atomic electron repulsion [email protected]

The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations

New Approaches to Large-Scale Electronic Structure Calculations

Recent years have witnessed a growing interest of the scientific community for the use of ab initio and density functional theory methods in theoretical studies of molecules containing many atoms. However, the `scaling wall' of some of the most accurate of such methods is often an obstacle for their applicability to systems of real-life interest, e.g. in biochemistry and nanotechnologies, one bottleneck being the evaluation and storage of the two-electron repulsion integrals. We have explored the possibility to avoid the expensive evaluation of the full integral matrix by designing and implementing approximate approaches based on Cholesky decomposition techniques. The results show the general applicability of this approximation to any quantum chemical model with substantial computational savings compared to conventional implementations -- not uncommonly of 1-2 orders of magnitude. At the same time, the loss of accuracy is minimal and can be systematically reduced at the price of sustainable additional costs. A major theoretical achievement of this study has been the reformulation of the Cholesky approximation in terms of solution to a "density fitting" variational problem. This has lead to an elegant formulation of the analytic derivatives of the Cholesky vectors, a long-standing lack of this technique. Moreover, this observation has paved the way to a new generation of accurate density fitting approximations free from biases and derived fully ab initio

p-Benzoquinone in aqueous solution: Stark shifts in spectra, asymmetry in solvent structure

Results from a simulation of p-benzoquinone (PBQ) in water is presented. An explicit solvent representation is used together with a multiconfigurational ab initio quantum chemical method. The electronic n -> pi* transitions are studied in aqueous solution and the two such transitions are both blue-shifted but to different degree. Both non-equilibrium and many-body effects are found to have decisive influence on the solvation: despite stronger hydrogen bonding between solute and solvent in an excited state than in the ground state, there is a blue-shift, and the solvent structure around the non-polar PBQ is asymmetric, which is argued to come from special many-body effects. The unusual result of strengthened hydrogen bonds in the excited state that follows from an excitation of a non-bonding electron on a proton acceptor, is explained by the near-linear Stark shift that is present in the transition

Quartic scaling evaluation of canonical scaled opposite spin second-order Moller-Plesset correlation energy using Cholesky decompositions

The scaled opposite spin second-order Moller-Plesset (SOS-MP2) energy expression is reformulated using Cholesky decomposition of the amplitude matrix. The resulting algorithm requires an auxiliary basis or Cholesky representation of the two-electron integrals and shows fourth-order scaling with system size. Based on an analysis of operation counts, we estimate that the present approach is computationally advantageous compared to the analogous fourth-order algorithms that employ Laplace transforms. (C) 2007 Elsevier B.V. All rights reserved