1
|
Zhang X, Li C, Ye HZ, Berkelbach TC, Chan GKL. Performant automatic differentiation of local coupled cluster theories: Response properties and ab initio molecular dynamics. J Chem Phys 2024; 161:014109. [PMID: 38949583 DOI: 10.1063/5.0212274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2024] [Accepted: 06/12/2024] [Indexed: 07/02/2024] Open
Abstract
In this work, we introduce a differentiable implementation of the local natural orbital coupled cluster (LNO-CC) method within the automatic differentiation framework of the PySCFAD package. The implementation is comprehensively tuned for enhanced performance, which enables the calculation of first-order static response properties on medium-sized molecular systems using coupled cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. We evaluate the accuracy of our method by benchmarking it against the canonical CCSD(T) reference for nuclear gradients, dipole moments, and geometry optimizations. In addition, we demonstrate the possibility of property calculations for chemically interesting systems through the computation of bond orders and Mössbauer spectroscopy parameters for a [NiFe]-hydrogenase active site model, along with the simulation of infrared spectra via ab initio LNO-CC molecular dynamics for a protonated water hexamer.
Collapse
Affiliation(s)
- Xing Zhang
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Chenghan Li
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
2
|
Liao K, Ding L, Schilling C. Quantum Information Orbitals (QIO): Unveiling Intrinsic Many-Body Complexity by Compressing Single-Body Triviality. J Phys Chem Lett 2024; 15:6782-6790. [PMID: 38913404 DOI: 10.1021/acs.jpclett.4c01105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/25/2024]
Abstract
The simultaneous treatment of static and dynamic correlations in strongly correlated electron systems is a critical challenge. In particular, finding a universal scheme for identifying a single-particle orbital basis that minimizes the representational complexity of the many-body wave function is a formidable and longstanding problem. As a contribution toward its solution, we show that the total orbital correlation actually reveals and quantifies the intrinsic complexity of the wave function, once it is minimized via orbital rotations. To demonstrate the power of this concept in practice, an iterative scheme is proposed to optimize the orbitals by minimizing the total orbital correlation calculated by the tailored coupled cluster singles and doubles (TCCSD) ansatz. The optimized orbitals enable the limited TCCSD ansatz to capture more nontrivial information on the many-body wave function, indicated by the improved wave function and energy. An initial application of this scheme shows great improvement of TCCSD in predicting the singlet ground state potential energy curves of the strongly correlated C2 and Cr2 molecule.
Collapse
Affiliation(s)
- Ke Liao
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Lexin Ding
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| | - Christian Schilling
- Faculty of Physics, Arnold Sommerfeld Centre for Theoretical Physics (ASC), Ludwig-Maximilians-Universität München, Theresienstr. 37, 80333 München, Germany
- Munich Center for Quantum Science and Technology (MCQST), Schellingstrasse 4, 80799 München, Germany
| |
Collapse
|
3
|
Rahrt R, Hein-Janke B, Amarasinghe KN, Shafique M, Feldt M, Guo L, Harvey JN, Pollice R, Koszinowski K, Mata RA. The Fe-MAN Challenge: Ferrates-Microkinetic Assessment of Numerical Quantum Chemistry. J Phys Chem A 2024; 128:4663-4673. [PMID: 38832568 PMCID: PMC11182345 DOI: 10.1021/acs.jpca.4c01361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
Organometallic species, such as organoferrate ions, are prototypical nucleophiles prone to reacting with a wide range of electrophiles, including proton donors. In solution, the operation of dynamic equilibria and the simultaneous presence of several organometallic species severely complicate the analysis of these fundamentally important reactions. This can be overcome by gas-phase experiments on mass-selected ions, which allow for the determination of the microscopic reactivity of the target species. In this contribution, we focus on the reactivity of a series of trisarylferrate complexes toward 2,2,2-trifluoroethanol and 2,2-difluoroethanol. By means of mass-spectrometric measurements, we determined the experimental bimolecular rate constants kexp of the gas-phase protolysis reactions of the trisarylferrate anions FePh3- and FeMes3- with the aforementioned acids. Based on these experiments, we carried out a dual blind challenge, inviting theoretical groups to submit their best predictions for the activation barriers and/or theoretical rate constants ktheo. This provides a unique opportunity to evaluate different computational protocols under minimal bias and sets the stage for further benchmarking of quantum chemical methods and data-driven approaches in the future.
Collapse
Affiliation(s)
- Rene Rahrt
- Institut
für Organische und Biomolekulare Chemie, Universität Göttingen, Tammannstr. 2, Göttingen 37077, Germany
| | - Björn Hein-Janke
- Institut
für Physikalische Chemie, Universität
Göttingen, Tammannstr.
6, Göttingen 37077, Germany
| | - Kosala N. Amarasinghe
- Leibniz
Institute for Catalysis (LIKAT), Albert-Einstein-Str. 29A, Rostock 18059, Germany
| | - Muhammad Shafique
- Leibniz
Institute for Catalysis (LIKAT), Albert-Einstein-Str. 29A, Rostock 18059, Germany
| | - Milica Feldt
- Leibniz
Institute for Catalysis (LIKAT), Albert-Einstein-Str. 29A, Rostock 18059, Germany
| | - Luxuan Guo
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven B-3001, Belgium
| | - Jeremy N. Harvey
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, Leuven B-3001, Belgium
| | - Robert Pollice
- Stratingh
Institute for Chemistry, University of Groningen, Nijenborgh 4, Groningen 9747 AG, The
Netherlands
| | - Konrad Koszinowski
- Institut
für Organische und Biomolekulare Chemie, Universität Göttingen, Tammannstr. 2, Göttingen 37077, Germany
| | - Ricardo A. Mata
- Institut
für Physikalische Chemie, Universität
Göttingen, Tammannstr.
6, Göttingen 37077, Germany
| |
Collapse
|
4
|
Panadés-Barrueta RL, Duflot D, Soto J, Martínez-Núñez E, Peláez D. Automatic Determination of the Non-Covalent Stable Conformations of the NO 2-Pyrene Cluster in Full Dimensionality (81D) Using the vdW-TSSCDS Approach. Chemphyschem 2024; 25:e202301001. [PMID: 38662437 DOI: 10.1002/cphc.202301001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2023] [Revised: 02/27/2024] [Indexed: 05/24/2024]
Abstract
We present the detailed topographical characterisation (stationary points and minimum energy paths connecting them) of the full dimensional (81D) intermolecular potential energy surface associated with the non-covalent interactions between the NO2 radical and the pyrene (C16H10) molecule. The whole procedure is (quasi) fully automated. We have used our recent algorithm vdW-TSSCDS as implemented on the freely-available AutoMekin software package. To this end, a series of inexpensive classical trajectories using forces from a low-level (semi-empirical) theory are used to sample the configuration space of the system in the search for candidates to first order saddle points. These guess structures are determined by means of a graph-theory based algorithm using the concept of adjacency matrix. Low-level optimizations are followed by re-optimizations at a final high-level of theory (DFT and CCSD(T)-F12 in our case.). The resulting set of stationary points and paths connecting them constitutes the so-called reaction network. In the case of NO2-pyrene, this network exhibits four major basins which can be characterized by their point-group symmetry. A central one, of global C2 symmetry, comprises the global minimum (as well as all other permutationally related conformers) together with the corresponding C2v saddle points connecting them. This central basin is connected to three others of lower C1 symmetry. The latter can be distinguished by the projection of the position of the NO2 nitrogen atom on the pyrene plane in combination with the relative orientation of the oxygen pair pointing either inwards, outwards, upwards or downwards.
Collapse
Affiliation(s)
- Ramón L Panadés-Barrueta
- Faculty of Chemistry and Food Chemistry, Technische Universität Dresden, 01069, Dresden, Germany
| | - Denis Duflot
- Univ. Lille, CNRS, UMR 8523, PhLAM - Physique des Lasers Atomes et Molécules, F-59000, Lille, France
| | - Juan Soto
- Departamento de Química Física, Facultad de Ciencias, Universidad de Málaga, Málaga, Spain
| | - Emilio Martínez-Núñez
- Departamento de Química Física, Facultade de Química, Universidade de Santiago de Compostela, 15782, Santiago de Compostela, Spain
| | - Daniel Peláez
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, 91405, Orsay, France
| |
Collapse
|
5
|
Leng JG, Sharples TR, Fournier M, McKendrick KG, Craciunescu L, Paterson MJ, Costen ML. Inelastic scattering of NO(A 2Σ +) + CO 2: rotation-rotation pair-correlated differential cross sections. Faraday Discuss 2024. [PMID: 38757419 DOI: 10.1039/d3fd00162h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2024]
Abstract
A crossed beam velocity-map ion-imaging apparatus has been used to determine differential cross sections (DCSs) for the rotationally inelastic scattering of NO(A2Σ+, v = 0, j = 0.5) with CO2, as a function of both NO(A, v = 0, N') final state and the coincident final rotational energy of the CO2. The DCSs are dominated by forward-peaked scattering for all N', with significant rotational excitation of CO2, and a small backward scattered peak is also observed for all final N'. However, no rotational rainbow scattering is observed and there is no evidence for significant product rotational angular momentum polarization. New ab initio potential energy surface calculations at the PNO-CCSD(T)-F12b level of theory report strong attractive forces at long ranges with significant anisotropy relative to both NO and CO2. The absence of rotational rainbow scattering is consistent with removal of low-impact-parameter collisions via electronic quenching, in agreement with the literature quenching rates of NO(A) by CO2 and recent electronic structure calculations. We propose that high-impact-parameter collisions, that do not lead to quenching, experience strong anisotropic attractive forces that lead to significant rotational excitation in both NO and CO2, depolarizing product angular momentum while leading to forward and backward glory scattering.
Collapse
Affiliation(s)
- Joseph G Leng
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, EH14 4AS, UK.
| | - Thomas R Sharples
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, EH14 4AS, UK.
| | - Martin Fournier
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, EH14 4AS, UK.
| | - Kenneth G McKendrick
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, EH14 4AS, UK.
| | - Luca Craciunescu
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, EH14 4AS, UK.
| | - Martin J Paterson
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, EH14 4AS, UK.
| | - Matthew L Costen
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, EH14 4AS, UK.
| |
Collapse
|
6
|
Gasevic T, Bursch M, Ma Q, Grimme S, Werner HJ, Hansen A. The p-block challenge: assessing quantum chemistry methods for inorganic heterocycle dimerizations. Phys Chem Chem Phys 2024; 26:13884-13908. [PMID: 38661329 DOI: 10.1039/d3cp06217a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
The elements of the p-block of the periodic table are of high interest in various chemical and technical applications like frustrated Lewis-pairs (FLP) or opto-electronics. However, high-quality benchmark data to assess approximate density functional theory (DFT) for their theoretical description are sparse. In this work, we present a benchmark set of 604 dimerization energies of 302 "inorganic benzenes" composed of all non-carbon p-block elements of main groups III to VI up to polonium. This so-called IHD302 test set comprises two classes of structures formed by covalent bonding and by weaker donor-acceptor (WDA) interactions, respectively. Generating reliable reference data with ab initio methods is challenging due to large electron correlation contributions, core-valence correlation effects, and especially the slow basis set convergence. To compute reference values for these dimerization reactions, after thorough testing, we applied a computational protocol using state-of-the-art explicitly correlated local coupled cluster theory termed PNO-LCCSD(T)-F12/cc-VTZ-PP-F12(corr.). It includes a basis set correction at the PNO-LMP2-F12/aug-cc-pwCVTZ level. Based on these reference data, we assess 26 DFT methods in combination with three different dispersion corrections and the def2-QZVPP basis set, five composite DFT approaches, and five semi-empirical quantum mechanical methods. For the covalent dimerizations, the r2SCAN-D4 meta-GGA, the r2SCAN0-D4 and ωB97M-V hybrids, and the revDSD-PBEP86-D4 double-hybrid functional are found to be the best-performing methods among the evaluated functionals of the respective class. However, since def2 basis sets for the 4th period are not associated to relativistic pseudo-potentials, we obtained significant errors in the covalent dimerization energies (up to 6 kcal mol-1) for molecules containing p-block elements of the 4th period. Significant improvements were achieved for systems containing 4th row elements by using ECP10MDF pseudopotentials along with re-contracted aug-cc-pVQZ-PP-KS basis sets introduced in this work with the contraction coefficients taken from atomic DFT (PBE0) calculations. Overall, the IHD302 set represents a challenge to contemporary quantum chemical methods. This is due to a large number of spatially close p-element bonds which are underrepresented in other benchmark sets, and the partial covalent bonding character for the WDA interactions. The IHD302 set may be helpful to develop more robust and transferable approximate quantum chemical methods in the future.
Collapse
Affiliation(s)
- Thomas Gasevic
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Markus Bursch
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
- FACCTs GmbH, 50677, Koeln, Germany
| | - Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany.
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany.
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| |
Collapse
|
7
|
Türkmen I, Dolg M. Linear Scaling Incremental Scheme for Correlation Energies with Embedding Generated Virtuals. J Chem Theory Comput 2024; 20:3154-3168. [PMID: 38588492 DOI: 10.1021/acs.jctc.3c01377] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/10/2024]
Abstract
A novel incremental scheme is presented including an incremental expansion of the virtual space for the calculation of electron correlation energies, which is compatible with any size-extensive correlation method and scales asymptotically linear for large molecules. The performance is studied for organic molecules, water clusters, and a La(III)-water complex, where the compatibility with pseudopotentials is also examined. The computational requirements are already reduced tremendously for medium-sized water clusters and hydrocarbons with respect to the canonical CCSD as well as the ordinary incremental scheme references. Correlation energies within chemical accuracy have been observed for all studied systems. The novelty of the method is that relatively small virtual spaces are used in combination with tuples of localized occupied spaces. The corresponding orthonormal occupied and virtual orbitals are obtained from QM/QM embedding calculations and can thus be used with standard quantum chemistry codes for correlation calculations. It is presented how relevant virtual spaces are selected and the correlation energies are linked in the new virtual space expansion.
Collapse
Affiliation(s)
- Ilyas Türkmen
- Institute for Theoretical Chemistry, University of Cologne, Greinstr. 4, D-50939 Cologne, Germany
| | - Michael Dolg
- Institute for Theoretical Chemistry, University of Cologne, Greinstr. 4, D-50939 Cologne, Germany
| |
Collapse
|
8
|
Mata RA, Zhanabekova T, Obenchain DA, Suhm MA. Dispersion Control over Molecule Cohesion: Exploiting and Dissecting the Tipping Power of Aromatic Rings. Acc Chem Res 2024; 57:1077-1086. [PMID: 38537179 PMCID: PMC11025128 DOI: 10.1021/acs.accounts.3c00664] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Revised: 02/26/2024] [Accepted: 02/27/2024] [Indexed: 04/17/2024]
Abstract
ConspectusWe have learned over the past years how London dispersion forces can be effectively used to influence or even qualitatively tip the structure of aggregates and the conformation of single molecules. This happens despite the fact that single dispersion contacts are much weaker than competing polar forces. It is a classical case of strength by numbers, with the importance of London dispersion forces scaling with the system size. Knowledge about the tipping points, however difficult to attain, is necessary for a rational design of intermolecular forces. One requires a careful assessment of the competing interactions, either by sensitive spectroscopic techniques for the study of the isolated molecules and aggregates or by theoretical approaches. Of particular interest are the systems close to the tipping point, when dispersion interactions barely outweigh or approach the strength of the other interactions. Such subtle cases are important milestones for a scale-up to realistic multi-interaction situations encountered in the fields of life and materials science. In searching for examples that provide ideal competing interactions in complexes and small clusters, aromatic systems can offer a diverse set of molecules with a variation of dispersion and electrostatic forces that control the dominant and peripheral interactions. Our combined spectroscopic and theoretical investigations provide valuable insights into the balance of intermolecular forces because they typically allow us to switch the aromatic substituent on and off. High-resolution rotational spectroscopy serves as a benchmark for molecular structures, as correct calculations should be based on correct geometries. When discussing the competition with other noncovalent interactions, obvious competitors are directional hydrogen bonds. As a second counterweight to aryl interactions, we will discuss aurophilic/metallophilic interactions, which also have a strong stabilization with a small number of atoms involved. Vibrational spectroscopy is most sensitive to interactions of light atoms, and the competition of OH hydrogen bonds with dispersion forces in a molecular aggregate can be judged well by the OH stretching frequency. Experiments in the gas phase are ideal for gauging the accuracy of quantum chemical predictions free of solvent forces. A tight collaboration utilizing these three methods allows experiment vs experiment vs theory benchmarking of the overall influence of dispersion in molecular structures and energetics.
Collapse
Affiliation(s)
- Ricardo A. Mata
- Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Tlektes Zhanabekova
- Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Daniel A. Obenchain
- Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Martin A. Suhm
- Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| |
Collapse
|
9
|
Sorathia K, Frantzov D, Tew DP. Improved CPS and CBS Extrapolation of PNO-CCSD(T) Energies: The MOBH35 and ISOL24 Data Sets. J Chem Theory Comput 2024; 20:2740-2750. [PMID: 38513261 PMCID: PMC11008106 DOI: 10.1021/acs.jctc.3c00974] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Revised: 03/06/2024] [Accepted: 03/08/2024] [Indexed: 03/23/2024]
Abstract
Computation of heats of reaction of large molecules is now feasible using the domain-based pair natural orbital (PNO)-coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] theory. However, to obtain agreement within 1 kcal/mol of experiment, it is necessary to eliminate basis set incompleteness error, which comprises both the AO basis set error and the PNO truncation error. Our investigation into the convergence to the canonical limit of PNO-CCSD(T) energies with the PNO truncation threshold T shows that errors follow the model E ( T ) = E + A T 1 / 2 . Therefore, PNO truncation errors can be eliminated using a simple two-point CPS extrapolation to the canonical limit so that subsequent CBS extrapolation is not limited by the residual PNO truncation error. Using the ISOL24 and MOBH35 data sets, we find that PNO truncation errors are larger for molecules with significant static correlation and that it is necessary to use very tight thresholds of T = 10 - 8 to ensure that errors do not exceed 1 kcal/mol. We present a lower-cost extrapolation scheme that uses information from small basis sets to estimate the PNO truncation errors for larger basis sets. In this way, the canonical limit of CCSD(T) calculations on sizable molecules with large basis sets can be reliably estimated in a practical way. Using this approach, we report near complete basis set (CBS)-CCSD(T) reaction energies for the full ISOL24 and MOBH35 data sets.
Collapse
Affiliation(s)
- Kesha Sorathia
- University of Oxford, South Parks Road, Oxford OX1 3QZ, U.K.
| | - Damyan Frantzov
- University of Oxford, South Parks Road, Oxford OX1 3QZ, U.K.
| | - David P. Tew
- University of Oxford, South Parks Road, Oxford OX1 3QZ, U.K.
| |
Collapse
|
10
|
Klostermann SV, Kappler J, Waigum A, Buchmeiser MR, Köhn A, Kästner J. The reduction behavior of sulfurized polyacrylonitrile (SPAN) in lithium-sulfur batteries using a carbonate electrolyte: a computational study. Phys Chem Chem Phys 2024; 26:9998-10007. [PMID: 38477497 DOI: 10.1039/d3cp06248a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/14/2024]
Abstract
Lithium-sulfur batteries (LSBs) have attracted attention due to their high theoretical energy density. This and various other advantages, such as the availability and non-toxicity of sulfur, raise interest in LSBs against the background of the energy revolution. However, a polysulfide shuttle mechanism can adversely affect the electrochemical performance of the cell. The sulfur redox properties are influenced, for example, by the electrolyte and the cathode material. Here, a computational study of the discharge process of an LSB with sulfurized poly(acrylonitrile) (SPAN) as the cathode material in combination with a carbonate electrolyte is presented. The nucleation of produced solid Li2S is compared to soluble Li2S. Dominating species are determined by comparing the Gibbs free energy of several species. We found that multiple lithiation steps occur before each Li2S detachment, preventing longer-chain polysulfide cleavage and a polysulfide shuttle. Through nucleating on the nitrogen-rich backbone of SPAN, Li2S units are stabilized by interactions with each other and with the nitrogen atoms. Experimental data show a potential drop and plateau during discharge, which is consistent with the calculated discharge profiles of SPAN with both soluble and nucleated Li2S, and hints at a direct solid-solid transition in the Li-SPAN cell during discharge when using carbonate-based electrolytes.
Collapse
Affiliation(s)
- S V Klostermann
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany.
| | - J Kappler
- Institute of Polymer Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - A Waigum
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany.
| | - M R Buchmeiser
- Institute of Polymer Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - A Köhn
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany.
| | - J Kästner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany.
| |
Collapse
|
11
|
Plett C, Grimme S, Hansen A. Conformational energies of biomolecules in solution: Extending the MPCONF196 benchmark with explicit water molecules. J Comput Chem 2024; 45:419-429. [PMID: 37982322 DOI: 10.1002/jcc.27248] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2023] [Revised: 10/16/2023] [Accepted: 10/17/2023] [Indexed: 11/21/2023]
Abstract
A prerequisite for the computational prediction of molecular properties like conformational energies of biomolecules is a reliable, robust, and computationally affordable method usually selected according to its performance for relevant benchmark sets. However, most of these sets comprise molecules in the gas phase and do not cover interactions with a solvent, even though biomolecules typically occur in aqueous solution. To address this issue, we introduce a with explicit water molecules solvated version of a gas-phase benchmark set containing 196 conformers of 13 peptides and other relevant macrocycles, namely MPCONF196 [J. Řezáč et al., JCTC 2018, 14, 1254-1266], and provide very accurate PNO-LCCSD(T)-F12b/AVQZ' reference values. The novel solvMPCONF196 benchmark set features two additional challenges beyond the description of conformers in the gas phase: conformer-water and water-water interactions. The overall best performing method for this set is the double hybrid revDSDPBEP86-D4/def2-QZVPP yielding conformational energies of almost coupled cluster quality. Furthermore, some (meta-)GGAs and hybrid functionals like B97M-V and ω B97M-D with a large basis set reproduce the coupled cluster reference with an MAD below 1 kcal mol- 1 . If more efficient methods are required, the composite DFT-method r2 SCAN-3c (MAD of 1.2 kcal mol- 1 ) is a good alternative, and when conformational energies of polypeptides or macrocycles with more than 500-1000 atoms are in the focus, the semi-empirical GFN2-xTB or the MMFF94 force field (for very large systems) are recommended.
Collapse
Affiliation(s)
- Christoph Plett
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Universität Bonn, Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Universität Bonn, Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Universität Bonn, Bonn, Germany
| |
Collapse
|
12
|
Hillers-Bendtsen AE, Jensen F, Mikkelsen KV, Olsen J, Jørgensen P. Cluster perturbation theory IX: Perturbation series for the coupled cluster singles and doubles ground state energy. J Chem Phys 2024; 160:104108. [PMID: 38477336 DOI: 10.1063/5.0192388] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2023] [Accepted: 02/17/2024] [Indexed: 03/14/2024] Open
Abstract
In this paper, we develop and analyze a number of perturbation series that target the coupled cluster singles and doubles (CCSD) ground state energy. We show how classical Møller-Plesset perturbation theory series can be restructured to target the CCSD energy based on a reference CCS calculation and how the corresponding cluster perturbation series differs from the classical Møller-Plesset perturbation series. Subsequently, we reformulate these series using the coupled cluster Lagrangian framework to obtain series, where fourth and fifth order energies are determined only using parameters through second order. To test the methods, we perform a series of test calculations on molecular photoswitches of both total energies and reaction energies. We find that the fifth order reaction energies are of CCSD quality and that they are of comparable accuracy to state-of-the-art approximations to the CCSD energy based on local pair natural orbitals. The advantage of the present approach over local correlation methods is the absence of user defined threshold parameters for neglecting or approximating contributions to the correlation energy. Fixed threshold parameters lead to discontinuous energy surfaces, although this effect is often small enough to be ignored, but the present approach has a differentiable energy that will facilitate derivation and implementation of gradients and higher derivatives. A further advantage is that the calculation of the perturbation correction is non-iterative and can, therefore, be calculated in parallel, leading to a short time-to-solution.
Collapse
Affiliation(s)
| | - Frank Jensen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK 8000 Aarhus C, Denmark
| | - Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK 2100 Copenhagen Ø, Denmark
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK 8000 Aarhus C, Denmark
| | - Poul Jørgensen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK 8000 Aarhus C, Denmark
| |
Collapse
|
13
|
Gasevic T, Bamberg M, Wicke J, Bolte M, Virovets A, Lerner HW, Grimme S, Hansen A, Wagner M, Bursch M. Confined Lewis Pairs: Investigation of the X - →Si 20 Interaction in Halogen-Encapsulating Silafulleranes. Angew Chem Int Ed Engl 2024; 63:e202314238. [PMID: 38059423 DOI: 10.1002/anie.202314238] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2023] [Revised: 12/04/2023] [Accepted: 12/05/2023] [Indexed: 12/08/2023]
Abstract
A joint theoretical and experimental study on 32 endohedral silafullerane derivatives [X@Si20 Y20 ]- (X=F-I; Y=F-I, H, Me, Et) andT h ${T_h }$ -[Cl@Si20 H12 Y8 ]- (Y=F-I) is presented. First, we evaluated the structure-determining template effect of Cl- in a systematic series of concave silapolyquinane model systems. Second, we investigated the X- →Si20 interaction energy (E int ${E_{{\rm{int}}} }$ ) as a function of X- and Y and found the largestE int ${E_{{\rm{int}}} }$ values for electron-withdrawing exohedral substituents Y. Given that X- ions can be considered as Lewis bases and empty Si20 Y20 clusters as Lewis acids, we classify our inseparable host-guest complexes [X@Si20 Y20 ]- as "confined Lewis pairs". Third, 35 Cl NMR spectroscopy proved to be highly diagnostic for an experimental assessment of the Cl- →Si20 interaction as the paramagnetic shielding and, in turn, δ ${\delta }$ (35 Cl) of the endohedral Cl- ion correlate inversely withE int ${E_{{\rm{int}}} }$ . Finally, we disclose the synthesis of [PPN][Cl@Si20 Y20 ] (Y=Me, Et, Br) and provide a thorough characterization of these new silafulleranes.
Collapse
Affiliation(s)
- Thomas Gasevic
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstraße 4, 53115, Bonn, Germany
| | - Marcel Bamberg
- Institut für Anorganische und Analytische Chemie, Goethe-Universität Frankfurt am Main, Max-von-Laue-Straße 7, 60438, Frankfurt am Main, Germany
| | - Julius Wicke
- Institut für Anorganische und Analytische Chemie, Goethe-Universität Frankfurt am Main, Max-von-Laue-Straße 7, 60438, Frankfurt am Main, Germany
| | - Michael Bolte
- Institut für Anorganische und Analytische Chemie, Goethe-Universität Frankfurt am Main, Max-von-Laue-Straße 7, 60438, Frankfurt am Main, Germany
| | - Alexander Virovets
- Institut für Anorganische und Analytische Chemie, Goethe-Universität Frankfurt am Main, Max-von-Laue-Straße 7, 60438, Frankfurt am Main, Germany
| | - Hans-Wolfram Lerner
- Institut für Anorganische und Analytische Chemie, Goethe-Universität Frankfurt am Main, Max-von-Laue-Straße 7, 60438, Frankfurt am Main, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstraße 4, 53115, Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstraße 4, 53115, Bonn, Germany
| | - Matthias Wagner
- Institut für Anorganische und Analytische Chemie, Goethe-Universität Frankfurt am Main, Max-von-Laue-Straße 7, 60438, Frankfurt am Main, Germany
| | - Markus Bursch
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470, Mülheim an der Ruhr, Germany
| |
Collapse
|
14
|
Paenurk E, Chen P. Robustness of Threshold Collision-Induced Dissociation Simulations for Bond Dissociation Energies. J Phys Chem A 2024; 128:333-342. [PMID: 38155581 DOI: 10.1021/acs.jpca.3c06862] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2023]
Abstract
The threshold collision-induced dissociation (T-CID) method is the workhorse for gas-phase bond dissociation energy (BDE) measurements. However, T-CID does not measure BDEs directly; instead, BDEs are obtained by fitting simulated data to the experimental data. We previously observed several large discrepancies between the computed and experimental BDEs. To analyze the reliability of the experimental values, we previously reported a study of the dissociation rate models in the simulation. Here, we report a study of the collision simulation part, specifically in the L-CID (ligand CID) program. We show that the BDE values are robust even to intentionally introduced mistakes in the simulations, varying in most cases by less than 3 kcal mol-1. The most significant exception is the collisional energy transfer (CET) simulation, which led to deviations larger than 10 kcal mol-1. However, we found that the BDEs obtained with explicitly simulated CET distributions deviated by only 3 kcal mol-1 from those simulated with the original model. Collectively, our results suggest that the T-CID-derived BDE values are robust and are likely to be accurate.
Collapse
Affiliation(s)
- Eno Paenurk
- Department of Chemistry and Applied Biosciences, ETH Zürich, Zürich 8093, Switzerland
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Peter Chen
- Department of Chemistry and Applied Biosciences, ETH Zürich, Zürich 8093, Switzerland
| |
Collapse
|
15
|
Bursch M, Grimme S, Hansen A. Influence of Steric and Dispersion Interactions on the Thermochemistry of Crowded (Fluoro)alkyl Compounds. Acc Chem Res 2024; 57:153-163. [PMID: 38102118 DOI: 10.1021/acs.accounts.3c00634] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2023]
Abstract
ConspectusAlkanes play a pivotal role in industrial, environmental, and biological processes. They are characterized by their carbon-carbon single-bond structure, remarkable stability, and conformational diversity. Fluorination of such compounds imparts unique physicochemical properties that often enhance pharmacokinetic profiles, metabolic stability, and receptor interactions while keeping beneficial properties. However, such per- and polyfluoroalkyl substances (PFAS) show a persistent presence in the environment and potential adverse health effects, which propelled them to the forefront of global environmental and health discussions. Alkyl compounds are also prototypical for stereoelectronic (SE) effects that are widely applied in chemistry. Substituents are typically described as electron-density-donating/withdrawing and/or responsible for sterically interacting with reagents or strategic groups in the molecule. That alkane branching can result in higher stability compared to less-branched isomers has been investigated in detail also by testing quantum chemical methods, in particular density functional theory (DFT). Alkane branching results in spatially compact structures with close intramolecular contacts so that at a specific size the detailed balance of attractive London dispersion and covalent versus repulsive Pauli exchange interactions shifts to new, chemically unfragile situations. This may lead to dissociation at room temperature and opens the central question: what is the smallest crowed alkane that cannot be made synthetically? In this Account, we try to shed light on the interplay among the various (free) energy components for crowded (fluoro)alkane dissociation. In this context, homolytic cleavage of the central C-C bond in a series of model alkanes of increasing size with tert-butyl (tBu), adamantyl (Ad), and [1.1.1]propellanyl (Prop) substituents is investigated. Reference energies are calculated at the PNO-LCCSD(T)-F12b level and used to benchmark the performance of contemporary DFT functionals. In line with previous conclusions, the application of dispersion corrections to density functionals is mandatory. For crowed structures, the accurate description of the midrange correlation effects, specifically repulsive van der Waals interactions, is crucial, and we observed that the density-dependent VV10 correction is superior to D4 in this context, although the asymptotic region is better described by the latter. The best available dispersion-inclusive functionals show systematic and reasonably small residual errors and can be safely applied to large systems (>100 atoms), for which coupled cluster methods with large basis sets are not computationally feasible anymore. For qualitatively correct predictions of synthetic accessibility under equilibrium conditions (free energy), the inclusion of thermostatistical (entropy) contributions is also essential. According to our results, tetra-tert-butylmethane (C17tBu) is the largest and most crowded system with a positive dissociation free energy and should be synthesizable. The difference between hydrogenated and perfluorinated systems originates from the increase in the steric repulsion of spatially close substituents, which is not compensated to the same extent by attractive orbital and dispersion interactions. A sometimes-assumed similar steric demand for fluorine and hydrogen atoms is not corroborated by our investigations on crowded systems. Perfluorination is found to substantially decrease thermal stability, rendering perfluorinated hexamethylethane (C8tBuF) the last potentially stable representative.
Collapse
Affiliation(s)
- Markus Bursch
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Beringstraße 4, D-53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Beringstraße 4, D-53115 Bonn, Germany
| |
Collapse
|
16
|
Wang Y, Guo Y, Neese F, Valeev EF, Li W, Li S. Cluster-in-Molecule Approach with Explicitly Correlated Methods for Large Molecules. J Chem Theory Comput 2023; 19:8076-8089. [PMID: 37920973 DOI: 10.1021/acs.jctc.3c00627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2023]
Abstract
In this article, we present a series of explicitly correlated local correlation methods developed under the cluster-in-molecule (CIM) framework, including explicitly correlated second-order Møller-Plesset perturbation (MP2), coupled-cluster singles and doubles (CCSD), domain-based local pair natural orbital CCSD (DLPNO-CCSD), and DLPNO-CCSD with perturbative triples (DLPNO-CCSD(T)). In these methods, F12 correction is decomposed into contributions from each occupied local molecular orbital and then evaluated independently in a given cluster, which consists of a subset of localized orbitals. These newly developed methods allow F12 calculations of large molecules (up to 145 atoms for quasi-one-dimensional systems) on a single node. We use these methods to investigate the relative stability between extended and folded alkane C30H62, the relative stability of four secondary structures of a polyglycine Ace(Gly)10NH2, and the binding energies of two host-guest complexes. The results demonstrate that the combination of CIM with F12 methods is a promising way to investigate large molecules with small basis set errors.
Collapse
Affiliation(s)
- Yuqi Wang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| |
Collapse
|
17
|
Shi B, Zen A, Kapil V, Nagy PR, Grüneis A, Michaelides A. Many-Body Methods for Surface Chemistry Come of Age: Achieving Consensus with Experiments. J Am Chem Soc 2023; 145:25372-25381. [PMID: 37948071 PMCID: PMC10683001 DOI: 10.1021/jacs.3c09616] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2023] [Revised: 10/15/2023] [Accepted: 10/17/2023] [Indexed: 11/12/2023]
Abstract
The adsorption energy of a molecule onto the surface of a material underpins a wide array of applications, spanning heterogeneous catalysis, gas storage, and many more. It is the key quantity where experimental measurements and theoretical calculations meet, with agreement being necessary for reliable predictions of chemical reaction rates and mechanisms. The prototypical molecule-surface system is CO adsorbed on MgO, but despite intense scrutiny from theory and experiment, there is still no consensus on its adsorption energy. In particular, the large cost of accurate many-body methods makes reaching converged theoretical estimates difficult, generating a wide range of values. In this work, we address this challenge, leveraging the latest advances in diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] to obtain accurate predictions for CO on MgO. These reliable theoretical estimates allow us to evaluate the inconsistencies in published temperature-programed desorption experiments, revealing that they arise from variations in employed pre-exponential factors. Utilizing this insight, we derive new experimental estimates of the (electronic) adsorption energy with a (more) precise pre-exponential factor. As a culmination of all of this effort, we are able to reach a consensus between multiple theoretical calculations and multiple experiments for the first time. In addition, we show that our recently developed cluster-based CCSD(T) approach provides a low-cost route toward achieving accurate adsorption energies. This sets the stage for affordable and reliable theoretical predictions of chemical reactions on surfaces to guide the realization of new catalysts and gas storage materials.
Collapse
Affiliation(s)
- Benjamin
X. Shi
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| | - Andrea Zen
- Dipartimento
di Fisica Ettore Pancini, Università
di Napoli Federico II, Monte S. Angelo, I-80126 Napoli, Italy
- Department
of Earth Sciences, University College London, Gower Street, WC1E 6BT London, U.K.
| | - Venkat Kapil
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| | - Péter R. Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Müegyetem rkp. 3, H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Andreas Grüneis
- Institute
for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| | - Angelos Michaelides
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, CB2 1EW Cambridge, U.K.
| |
Collapse
|
18
|
Werner HJ, Hansen A. Accurate Calculation of Isomerization and Conformational Energies of Larger Molecules Using Explicitly Correlated Local Coupled Cluster Methods in Molpro and ORCA. J Chem Theory Comput 2023; 19:7007-7030. [PMID: 37486154 DOI: 10.1021/acs.jctc.3c00270] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2023]
Abstract
An overview of the approximations in the explicitly correlated local coupled cluster methods PNO-LCCSD(T)-F12 in Molpro and DLPNO-CCSD(T)F12 in ORCA is given. Options to select the domains of projected atomic orbitals (PAOs), pair natural orbitals (PNOs), and triples natural orbitals (TNOs) in both programs are described and compared in detail. The two programs are applied to compute isomerization and conformational energies of the ISOL24 and ACONFL test sets, where the former is part of the GMTKN55 benchmark suite. Thorough studies of basis set effects are presented for selected systems. These revealed large intramolecular basis set superposition effects that make it practically impossible to reliably determine the complete basis set (CBS) limits without including explicitly correlated terms. The latter strongly reduce the basis set dependence and at the same time also errors caused by the local domain approximations. On the basis of these studies, the PNO-LCCSD(T)-F12 method is applied to determine new reference energies for the above-mentioned benchmark sets. We are confident that our results should agree within a few tenths of a kcal mol-1 with the (unknown) CCSD(T)/CBS values, which therefore allowed us to define computational settings for accurate explicitly correlated local coupled cluster methods with moderate computational effort. With these protocols, especially PNO-LCCSD(T)-F12b/AVTZ', reliable reference values for comprehensive benchmark sets can be generated efficiently. This can significantly advance the development and evaluation of the performance of approximate electronic structure methods, especially improved density functional approximations or machine learning approaches.
Collapse
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstrasse 4, D-53115 Bonn, Germany
| |
Collapse
|
19
|
Peyton BG, Wang Z, Crawford TD. Reduced Scaling Real-Time Coupled Cluster Theory. J Phys Chem A 2023; 127:8486-8499. [PMID: 37782945 DOI: 10.1021/acs.jpca.3c05151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/04/2023]
Abstract
Real-time coupled cluster (CC) methods have several advantages over their frequency-domain counterparts, namely, response and equation of motion CC theories. Broadband spectra, strong fields, and pulse manipulation allow for the simulation of complex spectroscopies that are unreachable using frequency-domain approaches. Due to the high-order polynomial scaling, the required numerical time propagation of the CC residual expressions is a computationally demanding process. This scaling may be reduced by local correlation schemes, which aim to reduce the size of the (virtual) orbital space by truncation according to user-defined parameters. We present the first application of local correlation to real-time CC. As in previous studies of locally correlated frequency-domain CC, traditional local correlation schemes are of limited utility for field-dependent properties; however, a perturbation-aware scheme proves promising. A detailed analysis of the amplitude dynamics suggests that the main challenge is a strong time dependence of the wave function sparsity.
Collapse
Affiliation(s)
- Benjamin G Peyton
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Zhe Wang
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - T Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
20
|
Craciunescu L, Liane EM, Kirrander A, Paterson MJ. Excited-state van der Waals potential energy surfaces for the NO A2Σ+ + CO2X1Σg+ collision complex. J Chem Phys 2023; 159:124303. [PMID: 38127380 DOI: 10.1063/5.0165769] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Accepted: 09/06/2023] [Indexed: 12/23/2023] Open
Abstract
Excited state van der Waals (vdW) potential energy surfaces (PESs) of the NO A2Σ+ + CO2X1Σg+ system are thoroughly investigated using coupled cluster theory and complete active space perturbation theory to second order (CASPT2). First, it is shown that pair natural orbital coupled cluster singles and doubles with perturbative triples yields comparable accuracy compared to CCSD(T) for molecular properties and vdW-minima at a fraction of computational cost of the latter. Using this method in conjunction with highly diffuse basis sets and counterpoise correction for basis set superposition error, the PESs for different intermolecular orientations are investigated. These show numerous vdW-wells, interconnected for all geometries except one, with a maximum depth of up to 830 cm-1; considerably deeper than those on the ground state surface. Multi-reference effects are investigated with CASPT2 calculations. The long-range vdW-surfaces support recent experimental observations relating to rotational energy transfer due the anisotropy in the potentials.
Collapse
Affiliation(s)
- Luca Craciunescu
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, Scotland EH14 4AS, United Kingdom
| | - Eirik M Liane
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, OX1 3QZ Oxford, United Kingdom
| | - Adam Kirrander
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, OX1 3QZ Oxford, United Kingdom
| | - Martin J Paterson
- Institute of Chemical Sciences, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, Scotland EH14 4AS, United Kingdom
| |
Collapse
|
21
|
Semidalas E, Martin JML. Correlation Consistent Basis Sets for Explicitly Correlated Theory: The Transition Metals. J Chem Theory Comput 2023; 19:5806-5820. [PMID: 37540641 PMCID: PMC10500978 DOI: 10.1021/acs.jctc.3c00506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2023] [Indexed: 08/06/2023]
Abstract
We present correlation consistent basis sets for explicitly correlated (F12) calculations, denoted VnZ(-PP)-F12-wis (n = D,T), for the d-block elements. The cc-pVDZ-F12-wis basis set is contracted to [8s7p5d2f] for the 3d-block, while its ECP counterpart for the 4d and 5d-blocks, cc-pVDZ-PP-F12-wis, is contracted to [6s6p5d2f]. The corresponding contracted sizes for cc-pVTZ(-PP)-F12-wis are [9s8p6d3f2g] for the 3d-block elements and [7s7p6d3f2g] for the 4d and 5d-block elements. Our VnZ(-PP)-F12-wis basis sets are evaluated on challenging test sets for metal-organic barrier heights (MOBH35) and group-11 metal clusters (CUAGAU-2). In F12 calculations, they are found to be about as close to the complete basis set limit as the combination of standard cc-pVnZ-F12 on main-group elements with the standard aug-cc-pV(n+1)Z(-PP) basis sets on the transition metal(s). While our basis sets are somewhat more compact than aug-cc-pV(n+1)Z(-PP), the CPU time benefit is negligible for catalytic complexes that contain only one or two transition metals among dozens of main-group elements; however, it is somewhat more significant for metal clusters.
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Molecular Chemistry
and Materials Science, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
| | - Jan M. L. Martin
- Department of Molecular Chemistry
and Materials Science, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
22
|
Mohammadi MD, Patsalidis N, Bhowmick S, Harmandaris VA, Biskos G. Adsorption of air pollutants onto silver and gold atomic clusters: DFT and PNO-LCCSD-F12 calculations. RSC Adv 2023; 13:18014-18024. [PMID: 37323454 PMCID: PMC10265589 DOI: 10.1039/d3ra01717f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Accepted: 05/05/2023] [Indexed: 06/17/2023] Open
Abstract
We provide a comprehensive investigation of intermolecular interactions between atmospheric gaseous pollutants, including CH4, CO, CO2, NO, NO2, SO2, as well as H2O and Agn (n = 1-22) or Aun (n = 1-20) atomic clusters. The optimized geometries of all the systems investigated in our study were determined using density functional theory (DFT) with M06-2X functional and SDD basis set. The PNO-LCCSD-F12/SDD method was used for more accurate single-point energy calculations. Compared to their isolated states, the structures of the Agn and Aun clusters undergo severe deformations upon adsorption of the gaseous species, which become more significant as the size of the clusters decreases. Considering that, in addition to adsorption energy, we have determined the interaction and deformation energy of all the systems. All our calculations consistently show that among the gaseous species examined, SO2 and NO2 exhibit a higher preference for adsorption on both types of clusters, with a slightly higher preference for the Ag clusters compared to the Au clusters, with the SO2/Ag16 system exhibiting the lowest adsorption energy. The type of intermolecular interactions was investigated through wave function analyses, including natural bond orbital (NBO) and quantum theory of atoms in molecules (QTAIM), showing that NO2 and SO2 are chemisorbed on the Agn and Aun atomic clusters, whereas the other gas molecules exhibit a much weaker interaction with them. The reported data can be used as input parameters for molecular dynamics simulations to study the selectivity of atomic clusters towards specific gases under ambient conditions, as well as to design materials that take advantage of the studied intermolecular interactions.
Collapse
Affiliation(s)
| | - Nikolaos Patsalidis
- Computation-based Science and Technology Research Centre, The Cyprus Institute Nicosia 2121 Cyprus
| | - Somnath Bhowmick
- Climate and Atmosphere Research Centre, The Cyprus Institute Nicosia 2121 Cyprus
| | - Vagelis A Harmandaris
- Computation-based Science and Technology Research Centre, The Cyprus Institute Nicosia 2121 Cyprus
- Institute of Applied and Computational Mathematics (IACM), Foundation for Research and Technology Hellas, (FORTH), IACM/FORTH GR-71110 Heraklion Crete Greece
- Department of Mathematics and Applied Mathematics, University of Crete GR-71409 Heraklion Crete Greece
| | - George Biskos
- Climate and Atmosphere Research Centre, The Cyprus Institute Nicosia 2121 Cyprus
- Faculty of Civil Engineering and Geosciences, Delft University of Technology Delft 2628 CN The Netherlands
| |
Collapse
|
23
|
Nash HW, Shaw RA, Hill JG. Correlation consistent auxiliary basis sets in density fitting Hartree-Fock: The atoms sodium through argon revisited. J Comput Chem 2023; 44:1119-1128. [PMID: 36636897 PMCID: PMC10947126 DOI: 10.1002/jcc.27069] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2022] [Revised: 11/29/2022] [Accepted: 12/14/2022] [Indexed: 01/14/2023]
Abstract
We present a series of auxiliary basis sets, for the elements Na to Ar, for use in density-fitted Hartree-Fock calculations with the correlation consistent cc-pV(n + d)Z orbital basis sets. Benchmarking on total molecular energies, reaction energies and the spectroscopic constants of the SO molecule demonstrate that the new sets address the deficiencies of using existing auxiliary sets in combination with these orbital basis sets. We also report auxiliary basis sets for Na and Mg matched to cc-pVnZ, along with recommendations for pairing auxiliary basis sets to the cc-pVnZ-F12 basis sets for Hartree-Fock calculations.
Collapse
Affiliation(s)
- Harry W. Nash
- Department of ChemistryUniversity of SheffieldSheffieldUK
| | - Robert A. Shaw
- Department of ChemistryUniversity of SheffieldSheffieldUK
| | | |
Collapse
|
24
|
Altun A, Riplinger C, Neese F, Bistoni G. Exploring the Accuracy Limits of PNO-Based Local Coupled-Cluster Calculations for Transition-Metal Complexes. J Chem Theory Comput 2023; 19:2039-2047. [PMID: 36917767 PMCID: PMC10100528 DOI: 10.1021/acs.jctc.3c00087] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/15/2023]
Abstract
While the domain-based local pair natural orbital coupled-cluster method with singles, doubles, and perturbative triples (DLPNO-CCSD(T)) has proven instrumental for computing energies and properties of large and complex systems accurately, calculations on first-row transition metals with a complex electronic structure remain challenging. In this work, we identify and address the two main error sources that influence the DLPNO-CCSD(T) accuracy in this context, namely, (i) correlation effects from the 3s and 3p semicore orbitals and (ii) dynamic correlation-induced orbital relaxation (DCIOR) effects that are not described by the local MP2 guess. We present a computational strategy that allows us to completely eliminate the DLPNO error associated with semicore correlation effects, while increasing, at the same time, the efficiency of the method. As regards the DCIOR effects, we introduce a diagnostic for estimating the deviation between DLPNO-CCSD(T) and canonical CCSD(T) for systems with significant orbital relaxation.
Collapse
Affiliation(s)
- Ahmet Altun
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | | | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Department of Chemistry, Biology, and Biotechnology, University of Perugia, 06123 Perugia, Italy
| |
Collapse
|
25
|
Kállay M, Horváth RA, Gyevi-Nagy L, Nagy PR. Basis Set Limit CCSD(T) Energies for Extended Molecules via a Reduced-Cost Explicitly Correlated Approach. J Chem Theory Comput 2022; 19:174-189. [PMID: 36576419 PMCID: PMC9835832 DOI: 10.1021/acs.jctc.2c01031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Several approximations are introduced and tested to reduce the computational expenses of the explicitly correlated coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method for both closed and open-shell species. First, the well-established frozen natural orbital (FNO) technique is adapted to explicitly correlated CC approaches. Second, our natural auxiliary function (NAF) scheme is employed to reduce the size of the auxiliary basis required for the density fitting approximation regularly used in explicitly correlated calculations. Third, a new approach, termed the natural auxiliary basis (NAB) approximation, is proposed to decrease the size of the auxiliary basis needed for the expansion of the explicitly correlated geminals. The performance of the above approximations and that of the combined FNO-NAF-NAB approach are tested for atomization and reaction energies. Our results show that overall speedups of 7-, 5-, and 3-times can be achieved with double-, triple-, and quadruple-ζ basis sets, respectively, without any loss in accuracy. The new method can provide, e.g., reaction energies and barrier heights well within chemical accuracy for molecules with more than 40 atoms within a few days using a few dozen processor cores, and calculations with 50+ atoms are still feasible. These routinely affordable computations considerably extend the reach of explicitly correlated CCSD(T).
Collapse
Affiliation(s)
- Mihály Kállay
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,
| | - Réka A. Horváth
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - László Gyevi-Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Péter R. Nagy
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
26
|
Wang Y, Ni Z, Neese F, Li W, Guo Y, Li S. Cluster-in-Molecule Method Combined with the Domain-Based Local Pair Natural Orbital Approach for Electron Correlation Calculations of Periodic Systems. J Chem Theory Comput 2022; 18:6510-6521. [PMID: 36240189 DOI: 10.1021/acs.jctc.2c00412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The cluster-in-molecule (CIM) method was extended to systems with periodic boundary conditions (PBCs) in a previous work (PBC-CIM) [J. Chem. Theory Comput.2019, 15, 2933], which is able to compute the electronic structures of periodic systems at second-order Møller-Plesset perturbation theory (MP2) and coupled cluster singles and doubles (CCSD) levels. However, the high computational costs of CCSD with respect to the size of clusters limit the usage of PBC-CIM to crystals with small or medium unit cells. In this work, we further develop the PBC-CIM method by employing the domain-based local pair natural orbital (DLPNO) methods for the electron correlation calculations of clusters to reduce the computational costs. The combined approach allows CCSD with perturbative triples, denoted as CCSD(T), to be computationally available for accurate descriptions of periodic systems. The distant-pair correction is also implemented to improve the accuracy of PBC-CIM. As in the molecular cases, the distant pair correction significantly improves the accuracy of various PBC-CIM methods with few additional costs. The PBC-CIM-DLPNO-CCSD(T) approach has been applied to investigate the optimized lattice parameter of the cubic LiCl crystal and two adsorption problems (CO on the NaCl(100) surface and H2O on the h-BN surface). The results show that the CIM-DLPNO-CCSD(T) method offers accurate and efficient descriptions for the studied systems. Another application to the cohesive energy of the acetic acid crystal reveals that large basis sets are necessary for reliable calculations on the cohesive energies of molecular crystals.
Collapse
Affiliation(s)
- Yuqi Wang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| | - Zhigang Ni
- College of Material, Chemistry and Chemical Engineering, Key Laboratory of Organosilicon Chemistry and Material Technology of Ministry of Education, Hangzhou Normal University, Hangzhou311121, P. R. China
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, Mülheim an der RuhrD-45470, Germany
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong266237, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| |
Collapse
|
27
|
Saitow M, Uemura K, Yanai T. A local pair-natural orbital-based complete-active space perturbation theory using orthogonal localized virtual molecular orbitals. J Chem Phys 2022; 157:084101. [DOI: 10.1063/5.0094777] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The multireference second-order perturbation theory (CASPT2) is known to deliver a quantitative description of various complex electronic states. Despite its near-size-consistent nature, the applicability of the CASPT2 method to large, real-life systems is mostly hindered by large computational and storage costs for the two-external tensors, such as two-electron integrals, amplitudes, and residuum. To this end, Menezes and co-workers developed a reduced-scaling CASPT2 scheme by incorporating the local pair-natural orbital (PNO) representation of the many-body wave functions using non-orthonormal projected atomic orbitals (PAOs) into the CASPT theory [F. Menezes et al., J. Chem. Phys. 145, 124115 (2016)]. Alternatively, in this paper, we develop a new PNO-based CASPT2 scheme using the orthonormal localized virtual molecular orbitals (LVMOs) and assess its performance and accuracy in comparison with the conventional PAO-based counterpart. Albeit the compactness, the LVMOs were considered to perform somewhat poorly compared to PAOs in the local correlation framework because they caused enormously large orbital domains. In this work, we show that the size of LVMO domains can be rendered comparable to or even smaller than that of PAOs by the use of the differential overlap integrals for domain construction. Optimality of the MOs from the CASSCF treatment is a key to reducing the LVMO domain size for the multireference case. Due to the augmented Hessian-based localization algorithm, an additional computational cost for obtaining the LVMOs is relatively minor. We demonstrate that the LVMO-based PNO-CASPT2 method is routinely applicable to large, real-life molecules such as Menshutkin SN2 reaction in a single-walled carbon nanotube reaction field.
Collapse
Affiliation(s)
- Masaaki Saitow
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| | - Kazuma Uemura
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| | - Takeshi Yanai
- Department of Chemistry, Graduate School of Science, Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
- Institute of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furocho, Chikusa Ward, Nagoya, Aichi 464-8601, Japan
| |
Collapse
|
28
|
Menezes F, Popowicz GM. ULYSSES: An Efficient and Easy to Use Semiempirical Library for C+. J Chem Inf Model 2022; 62:3685-3694. [PMID: 35930308 DOI: 10.1021/acs.jcim.2c00757] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We introduce ULYSSES, a user-friendly and robust C++ library for semiempirical quantum chemical calculations. In its current version, ULYSSES is equipped with a large set of different semiempirical models, most of which are based on the Neglect of Diatomic Differential Overlap (NDDO) approximation. Empirical corrections for dispersion and hydrogen bonding are available for most methods, so that higher quality is achieved in the calculation of energies of nonbonded complexes. The library is furthermore equipped with geometry optimization, as well as modules for calculating molecular properties of general interest. Ideal gas thermodynamics is available and allows single structure as well as conformer (multistructure) averaged properties to be calculated. We offer the possibility to use several vibrational partition functions according to the nature of interactions being studied: for covalent systems, the traditional harmonic oscillator approximation is available; for nonbonded complexes, we systematically extended the partition function proposed by Grimme for all thermodynamic functions. The library is also capable of running Born-Oppenheimer molecular dynamics.
Collapse
Affiliation(s)
- Filipe Menezes
- Institute of Structural Biology, Helmholtz Zentrum Muenchen, Ingolstaedter Landstrasse 1, 85764 Neuherberg, Germany
| | - Grzegorz M Popowicz
- Institute of Structural Biology, Helmholtz Zentrum Muenchen, Ingolstaedter Landstrasse 1, 85764 Neuherberg, Germany
| |
Collapse
|
29
|
Altun A, Neese F, Bistoni G. Open-Shell Variant of the London Dispersion-Corrected Hartree-Fock Method (HFLD) for the Quantification and Analysis of Noncovalent Interaction Energies. J Chem Theory Comput 2022; 18:2292-2307. [PMID: 35167304 PMCID: PMC9009084 DOI: 10.1021/acs.jctc.1c01295] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The London dispersion
(LD)-corrected Hartree–Fock (HF) method
(HFLD) is an ab initio approach for the quantification
and analysis of noncovalent interactions (NCIs) in large systems that
is based on the domain-based local pair natural orbital coupled-cluster
(DLPNO-CC) theory. In the original HFLD paper, we discussed the implementation,
accuracy, and efficiency of its closed-shell variant. Herein, an extension
of this method to open-shell molecular systems is presented. Its accuracy
is tested on challenging benchmark sets for NCIs, using CCSD(T) energies
at the estimated complete basis set limit as reference. The HFLD scheme
was found to be as accurate as the best-performing dispersion-corrected
exchange-correlation functionals, while being nonempirical and equally
efficient. In addition, it can be combined with the well-established
local energy decomposition (LED) for the analysis of NCIs, thus yielding
additional physical insights.
Collapse
Affiliation(s)
- Ahmet Altun
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany.,Department of Chemistry, Biology and Biotechnology, University of Perugia, 06123 Perugia, Italy
| |
Collapse
|
30
|
Rudshteyn B, Weber JL, Coskun D, Devlaminck PA, Zhang S, Reichman DR, Shee J, Friesner RA. Calculation of Metallocene Ionization Potentials via Auxiliary Field Quantum Monte Carlo: Toward Benchmark Quantum Chemistry for Transition Metals. J Chem Theory Comput 2022; 18:2845-2862. [PMID: 35377642 PMCID: PMC9123894 DOI: 10.1021/acs.jctc.1c01071] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
The accurate ab initio prediction of ionization energies is essential to understanding the electrochemistry of transition metal complexes in both materials science and biological applications. However, such predictions have been complicated by the scarcity of gas phase experimental data, the relatively large size of the relevant molecules, and the presence of strong electron correlation effects. In this work, we apply all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multideterminant trial wave functions to six metallocene complexes to compare the computed adiabatic and vertical ionization energies with experimental results. We find that ph-AFQMC yields mean absolute errors (MAEs) of 1.69 ± 1.02 kcal/mol for the adiabatic energies and 2.85 ± 1.13 kcal/mol for the vertical energies. We also carry out density functional theory (DFT) calculations using a variety of functionals, which yields MAEs of 3.62-6.98 kcal/mol and 3.31-9.88 kcal/mol, as well as one variant of localized coupled cluster calculations (DLPNO-CCSD(T0) with moderate PNO cutoffs), which has MAEs of 4.96 and 6.08 kcal/mol, respectively. We also test the reliability of DLPNO-CCSD(T0) and DFT on acetylacetonate (acac) complexes for adiabatic energies measured in the same manner experimentally, and we find higher MAEs, ranging from 4.56 to 10.99 kcal/mol (with a different ordering) for DFT and 6.97 kcal/mol for DLPNO-CCSD(T0). Finally, by utilizing experimental solvation energies, we show that accurate reduction potentials in solution for the metallocene series can be obtained from the AFQMC gas phase results.
Collapse
Affiliation(s)
- Benjamin Rudshteyn
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Dilek Coskun
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Pierre A Devlaminck
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 Fifth Avenue, New York, New York 10010, United States.,Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
31
|
Meuwly M. Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─ Quo Vadis?. J Phys Chem B 2022; 126:2155-2167. [PMID: 35286087 DOI: 10.1021/acs.jpcb.2c00212] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.
Collapse
Affiliation(s)
- Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, 4056 Basel, Switzerland
| |
Collapse
|
32
|
Rask AE, Zimmerman PM. The many-body electronic interactions of Fe(II)–porphyrin. J Chem Phys 2022; 156:094110. [DOI: 10.1063/5.0079310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Fe(II)–porphyrin complexes exhibit a diverse range of electronic interactions between the metal and macrocycle. Herein, the incremental full configuration interaction method is applied to the entire space of valence orbitals of a Fe(II)–porphyrin model using a modest basis set. A novel visualization framework is proposed to analyze individual many-body contributions to the correlation energy, providing detailed maps of this complex’s highly correlated electronic structure. This technique is used to parse the numerous interactions of two low-lying triplet states (3A2g and 3Eg) and to show that strong metal d–d and macrocycle π–π orbital interactions preferentially stabilize the 3A2g state. d–π interactions, on the other hand, preferentially stabilize the 3Eg state and primarily appear when correlating six electrons at a time. Ultimately, the Fe(II)–porphyrin model’s full set of 88 valence electrons are correlated in 275 orbitals, showing the interactions up to the 4-body level, which covers the great majority of correlations in this system.
Collapse
Affiliation(s)
- A. E. Rask
- Department of Chemistry, University of Michigan, 930 N. University Ave., Ann Arbor, Michigan 48109, USA
| | - P. M. Zimmerman
- Department of Chemistry, University of Michigan, 930 N. University Ave., Ann Arbor, Michigan 48109, USA
| |
Collapse
|
33
|
Hohenstein EG, Fales BS, Parrish RM, Martínez TJ. Rank-reduced coupled-cluster. III. Tensor hypercontraction of the doubles amplitudes. J Chem Phys 2022; 156:054102. [DOI: 10.1063/5.0077770] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Edward G. Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | - B. Scott Fales
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| | | | - Todd J. Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| |
Collapse
|
34
|
Semidalas E, Martin JM. The MOBH35 Metal–Organic Barrier Heights Reconsidered: Performance of Local-Orbital Coupled Cluster Approaches in Different Static Correlation Regimes. J Chem Theory Comput 2022; 18:883-898. [PMID: 35045709 PMCID: PMC8830049 DOI: 10.1021/acs.jctc.1c01126] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
We have revisited
the MOBH35 (Metal–Organic Barrier Heights,
35 reactions) benchmark [Iron, Janes, J.
Phys. Chem. A, 2019, 123 ( (17), ), 3761−378130973722; ibid. 2019, 123, 6379–6380] for realistic organometallic catalytic reactions, using both canonical
CCSD(T) and localized orbital approximations to it. For low levels
of static correlation, all of DLPNO-CCSD(T), PNO-LCCSD(T), and LNO-CCSD(T)
perform well; for moderately strong levels of static correlation,
DLPNO-CCSD(T) and (T1) may break down catastrophically,
and PNO-LCCSD(T) is vulnerable as well. In contrast, LNO-CCSD(T) converges
smoothly to the canonical CCSD(T) answer with increasingly tight convergence
settings. The only two reactions for which our revised MOBH35 reference
values differ substantially from the original ones are reaction 9
and to a lesser extent 8, both involving iron. For the purpose of
evaluating density functional theory (DFT) methods for MOBH35, it
would be best to remove reaction 9 entirely as its severe level of
static correlation makes it just too demanding for a test. The magnitude
of the difference between DLPNO-CCSD(T) and DLPNO-CCSD(T1) is a reasonably good predictor for errors in DLPNO-CCSD(T1) compared to canonical CCSD(T); otherwise, monitoring all of T1, D1, max|tiA|, and 1/(εLUMO – εHOMO) should provide adequate warning
for potential problems. Our conclusions are not specific to the def2-SVP
basis set but are largely conserved for the larger def2-TZVPP, as
they are for the smaller def2-SV(P): the latter may be an economical
choice for calibrating against canonical CCSD(T). Finally, diagnostics
for static correlation are statistically clustered into groups corresponding
to (1) importance of single excitations in the wavefunction; (2a)
the small band gap, weakly separated from (2b) correlation entropy;
and (3) thermochemical importance of correlation energy, as well as
the slope of the DFT reaction energy with respect to the percentage
of HF exchange. Finally, a variable reduction analysis reveals that
much information on the multireference character is provided by T1, IND/Itot, and the exchange-based diagnostic A100[TPSS].
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot 7610001, Israel
| | - Jan M.L. Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot 7610001, Israel
| |
Collapse
|
35
|
Stöckl Y, Fellmeth T, Bauer F, Wank B, Frey W, Claasen B, Zens A, Köhn A, Laschat S. Chasing polycyclic natural products: 5/6/5‐ or 5/6/6‐carbotricyclic scaffold construction via stereodivergent Diels‐Alder reaction of chiral hydrindanes and their boron complexes. European J Org Chem 2022. [DOI: 10.1002/ejoc.202101416] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Yannick Stöckl
- Universität Stuttgart: Universitat Stuttgart Institut für Organische Chemie GERMANY
| | - Thomas Fellmeth
- University of Stuttgart: Universitat Stuttgart Institut für Theoretische Chemie GERMANY
| | - Florian Bauer
- University of Stuttgart: Universitat Stuttgart Institut für Theoretische Chemie GERMANY
| | - Bianca Wank
- Universität Stuttgart: Universitat Stuttgart Institut für Organische Chemie GERMANY
| | - Wolfgang Frey
- Universität Stuttgart: Universitat Stuttgart Institut für Organische Chemie GERMANY
| | - Birgit Claasen
- Universität Stuttgart: Universitat Stuttgart Institut für Organische Chemie GERMANY
| | - Anna Zens
- Universität Stuttgart: Universitat Stuttgart Institut für Organische Chemie GERMANY
| | - Andreas Köhn
- Universität Stuttgart: Universitat Stuttgart Institut für Theoretische Chemie GERMANY
| | - Sabine Laschat
- Universität Stuttgart Institut für Organische Chemie Pfaffenwaldring 55 70569 Stuttgart GERMANY
| |
Collapse
|
36
|
Altun A, Ghosh S, Riplinger C, Neese F, Bistoni G. Addressing the System-Size Dependence of the Local Approximation Error in Coupled-Cluster Calculations. J Phys Chem A 2021; 125:9932-9939. [PMID: 34730360 PMCID: PMC8607505 DOI: 10.1021/acs.jpca.1c09106] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Over the last two decades, the local approximation has been successfully used to extend the range of applicability of the "gold standard" singles and doubles coupled-cluster method with perturbative triples CCSD(T) to systems with hundreds of atoms. The local approximation error grows in absolute value with the increasing system size, i.e., by increasing the number of electron pairs in the system. In this study, we demonstrate that the recently introduced two-point extrapolation scheme for approaching the complete pair natural orbital (PNOs) space limit in domain-based pair natural orbital CCSD(T) calculations drastically reduces the dependence of the error on the system size, thus opening up unprecedented opportunities for the calculation of benchmark quality relative energies for large systems.
Collapse
Affiliation(s)
- Ahmet Altun
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Soumen Ghosh
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | | | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
37
|
Haldar S, Dutta AK. An efficient Fock space multi-reference coupled cluster method based on natural orbitals: Theory, implementation, and benchmark. J Chem Phys 2021; 155:014105. [PMID: 34241374 DOI: 10.1063/5.0054171] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a natural orbital-based implementation of the intermediate Hamiltonian Fock space coupled-cluster method for the (1, 1) sector of Fock space. The use of natural orbitals significantly reduces the computational cost and can automatically choose an appropriate set of active orbitals. The new method retains the charge transfer separability of the original intermediate Hamiltonian Fock space coupled-cluster method and gives excellent performance for valence, Rydberg, and charge-transfer excited states. It offers significant computational advantages over the popular equation of motion coupled cluster method for excited states dominated by single excitations.
Collapse
Affiliation(s)
- Soumi Haldar
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
38
|
Kowalski K, Bair R, Bauman NP, Boschen JS, Bylaska EJ, Daily J, de Jong WA, Dunning T, Govind N, Harrison RJ, Keçeli M, Keipert K, Krishnamoorthy S, Kumar S, Mutlu E, Palmer B, Panyala A, Peng B, Richard RM, Straatsma TP, Sushko P, Valeev EF, Valiev M, van Dam HJJ, Waldrop JM, Williams-Young DB, Yang C, Zalewski M, Windus TL. From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape. Chem Rev 2021; 121:4962-4998. [PMID: 33788546 DOI: 10.1021/acs.chemrev.0c00998] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems. As the hardware and software have evolved, so have the theoretical and computational chemistry methods and algorithms. Parallel computers clearly changed the common computing paradigm in the late 1970s and 80s, and the field has again seen a paradigm shift with the advent of graphical processing units. This review explores the challenges and some of the solutions in transforming software from the terascale to the petascale and now to the upcoming exascale computers. While discussing the field in general, NWChem and its redesign, NWChemEx, will be highlighted as one of the early codesign projects to take advantage of massively parallel computers and emerging software standards to enable large scientific challenges to be tackled.
Collapse
Affiliation(s)
- Karol Kowalski
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Raymond Bair
- Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Nicholas P Bauman
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | | - Eric J Bylaska
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jeff Daily
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Wibe A de Jong
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Thom Dunning
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Niranjan Govind
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Robert J Harrison
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, United States
| | - Murat Keçeli
- Argonne National Laboratory, Lemont, Illinois 60439, United States
| | | | | | - Suraj Kumar
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Erdal Mutlu
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Bruce Palmer
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Ajay Panyala
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Bo Peng
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | | - T P Straatsma
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6373, United States
| | - Peter Sushko
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Marat Valiev
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | | | | | | - Chao Yang
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Marcin Zalewski
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Theresa L Windus
- Department of Chemistry, Iowa State University and Ames Laboratory, Ames, Iowa 50011, United States
| |
Collapse
|
39
|
Paenurk E, Chen P. Modeling Gas-Phase Unimolecular Dissociation for Bond Dissociation Energies: Comparison of Statistical Rate Models within RRKM Theory. J Phys Chem A 2021; 125:1927-1940. [PMID: 33635061 DOI: 10.1021/acs.jpca.1c00183] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Rice-Ramsperger-Kassel-Marcus (RRKM) theory provides a simple yet powerful rate theory for calculating microcanonical rate constants. In particular, it has found widespread use in combination with gas-phase kinetic experiments of unimolecular dissociations to extract experimental bond dissociation energies (BDEs). We have previously found several discrepancies between the computed BDE values and the respective experimental ones, obtained with our empirical rate model, named L-CID. To investigate the reliability of our rate model, we conducted a theoretical analysis and comparison of the performance of conventional rate models and L-CID within the RRKM framework. Using the previously published microcanonical rate data as well as reaction cross-section data, we show that the BDE values obtained with the L-CID model agree with the ones from the other rate models within the expected uncertainty bounds. Based on this agreement, we discuss the possible rationalization of the good performance of the L-CID model.
Collapse
Affiliation(s)
- Eno Paenurk
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| | - Peter Chen
- Department of Chemistry and Applied Biosciences, ETH Zürich, 8093 Zürich, Switzerland
| |
Collapse
|
40
|
Ma Q, Werner HJ. Scalable Electron Correlation Methods. 8. Explicitly Correlated Open-Shell Coupled-Cluster with Pair Natural Orbitals PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. J Chem Theory Comput 2021; 17:902-926. [PMID: 33405921 DOI: 10.1021/acs.jctc.0c01129] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
We present explicitly correlated open-shell pair natural orbital local coupled-cluster methods, PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. The methods are extensions of our previously reported PNO-R/UCCSD methods (J. Chem. Theory Comput., 2020, 16, 3135-3151, https://pubs.acs.org/doi/10.1021/acs.jctc.0c00192) with additions of explicit correlation and perturbative triples corrections. The explicit correlation treatment follows the spin-orbital CCSD-F12b theory using Ansatz 3*A, which is found to yield comparable or better basis set convergence than the more rigorous Ansatz 3C in computed ionization potentials and reaction energies using double- to quaduple-ζ basis sets. The perturbative triples correction is adapted from the spin-orbital (T) theory to use triples natural orbitals (TNOs). To address the coupling due to off-diagonal Fock matrix elements, the local triples amplitudes are iteratively solved using small domains of TNOs, and a semicanonical (T0) domain correction with larger domains is applied to reduce the domain errors. The performance of the methods is demonstrated through benchmark calculations on ionization potentials, radical stabilization energies, reaction energies of fragmentations and rearrangements in radical cations, and spin-state energy differences of iron complexes. For a few test sets where canonical calculations are feasible, PNO-RCCSD(T)-F12 results agree with the canonical ones to within 0.4 kcal mol-1, and this maximum error is reduced to below 0.2 kcal mol-1 when large local domains are used. For larger systems, results using different thresholds for the local approximations are compared to demonstrate that 1 kcal mol-1 level of accuracy can be achieved using our default settings. For a couple of difficult cases, it is demonstrated that the errors from individual approximations are only a fraction of 1 kcal mol-1, and the overall accuracy of the method does not rely on error compensations. In contrast to canonical calculations, the use of spin-orbitals does not lead to a significant increase of computational time and memory usage in the most expensive steps of PNO-R/UCCSD(T)-F12 calculations. The only exception is the iterative solution of the (T) amplitudes, which can be avoided without significant errors by using a perturbative treatment of the off-diagonal coupling, known as (T1) approximation. For most systems, even the semicanonical approximation (T0) leads only to small errors in relative energies. Our program is well parallelized and capable of computing accurate correlation energies for molecules with 100-200 atoms using augmented triple-ζ basis sets in less than a day of elapsed time on a small computer cluster.
Collapse
Affiliation(s)
- Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
41
|
Tew DP. Principal domains in F12 explicitly correlated theory. ADVANCES IN QUANTUM CHEMISTRY 2021. [DOI: 10.1016/bs.aiq.2021.06.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
42
|
Mukherjee M, Tripathi D, Brehm M, Riplinger C, Dutta AK. Efficient EOM-CC-based Protocol for the Calculation of Electron Affinity of Solvated Nucleobases: Uracil as a Case Study. J Chem Theory Comput 2020; 17:105-116. [DOI: 10.1021/acs.jctc.0c00655] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
| | - Divya Tripathi
- Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Martin Brehm
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | | | | |
Collapse
|
43
|
Calvin JA, Peng C, Rishi V, Kumar A, Valeev EF. Many-Body Quantum Chemistry on Massively Parallel Computers. Chem Rev 2020; 121:1203-1231. [DOI: 10.1021/acs.chemrev.0c00006] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Justus A. Calvin
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Chong Peng
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Varun Rishi
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
44
|
Sorathia K, Tew DP. Basis set extrapolation in pair natural orbital theories. J Chem Phys 2020; 153:174112. [PMID: 33167642 DOI: 10.1063/5.0022077] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present the results of a benchmark study of the effect of Pair Natural Orbital (PNO) truncation errors on the performance of basis set extrapolation. We find that reliable conclusions from the application of Helgaker's extrapolation method are only obtained when using tight PNO thresholds of at least 10-7. The use of looser thresholds introduces a significant risk of observing a false basis set convergence and underestimating the residual basis set errors. We propose an alternative extrapolation approach based on the PNO truncation level that only requires a single basis set and show that it is a viable alternative to hierarchical basis set extrapolation methods.
Collapse
Affiliation(s)
- Kesha Sorathia
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, GermanyUniversity of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - David P Tew
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| |
Collapse
|
45
|
Scott CJC, Di Remigio R, Crawford TD, Thom AJW. Theory and implementation of a novel stochastic approach to coupled cluster. J Chem Phys 2020; 153:144117. [DOI: 10.1063/5.0026513] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Charles J. C. Scott
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| | - Roberto Di Remigio
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - T. Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
46
|
Altun A, Neese F, Bistoni G. Extrapolation to the Limit of a Complete Pair Natural Orbital Space in Local Coupled-Cluster Calculations. J Chem Theory Comput 2020; 16:6142-6149. [PMID: 32897712 PMCID: PMC7586325 DOI: 10.1021/acs.jctc.0c00344] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
![]()
The domain-based
local pair natural orbital (PNO) coupled-cluster
DLPNO-CCSD(T) method allows one to perform single point energy calculations
for systems with hundreds of atoms while retaining essentially the
accuracy of its canonical counterpart, with errors that are typically
smaller than 1 kcal/mol for relative energies. Crucial to the accuracy
and efficiency of the method is a proper definition of the virtual
space in which the coupled-cluster equations are solved, which is
spanned by a highly compact set of pair natural orbitals (PNOs) that
are specific for each electron pair. The dimension of the PNO space
is controlled by the TCutPNO threshold:
only PNOs with an occupation number greater than TCutPNO are included in the correlation space of a given
electron pair, whilst the remaining PNOs are discarded. To keep the
error of the method small, a conservative TCutPNO value is used in standard DLPNO-CCSD(T) calculations. This often
leads to unnecessarily large PNO spaces, which limits the efficiency
of the method. Herein, we introduce a new computational strategy to
approach the complete PNO space limit (for a given basis set) that
consists in extrapolating the results obtained with different TCutPNO values. The method is validated on the
GMTKN55 set using canonical CCSD(T) data as the reference. Our results
demonstrate that a simple two-point extrapolation scheme can be used
to significantly increase the efficiency and accuracy of DLPNO-CCSD(T)
calculations, thus extending the range of applicability of the technique.
Collapse
Affiliation(s)
- Ahmet Altun
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
47
|
Förster A, Visscher L. Double hybrid DFT calculations with Slater type orbitals. J Comput Chem 2020; 41:1660-1684. [PMID: 32297682 PMCID: PMC7317772 DOI: 10.1002/jcc.26209] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Revised: 03/31/2020] [Accepted: 04/01/2020] [Indexed: 12/14/2022]
Abstract
On a comprehensive database with 1,644 datapoints, covering several aspects of main-group as well as of transition metal chemistry, we assess the performance of 60 density functional approximations (DFA), among them 36 double hybrids (DH). All calculations are performed using a Slater type orbital (STO) basis set of triple-ζ (TZ) quality and the highly efficient pair atomic resolution of the identity approach for the exchange- and Coulomb-term of the KS matrix (PARI-K and PARI-J, respectively) and for the evaluation of the MP2 energy correction (PARI-MP2). Employing the quadratic scaling SOS-AO-PARI-MP2 algorithm, DHs based on the spin-opposite-scaled (SOS) MP2 approximation are benchmarked against a database of large molecules. We evaluate the accuracy of STO/PARI calculations for B3LYP as well as for the DH B2GP-PLYP and show that the combined basis set and PARI-error is comparable to the one obtained using the well-known def2-TZVPP Gaussian-type basis set in conjunction with global density fitting. While quadruple-ζ (QZ) calculations are currently not feasible for PARI-MP2 due to numerical issues, we show that, on the TZ level, Jacob's ladder for classifying DFAs is reproduced. However, while the best DHs are more accurate than the best hybrids, the improvements are less pronounced than the ones commonly found on the QZ level. For conformers of organic molecules and noncovalent interactions where very high accuracy is required for qualitatively correct results, DHs provide only small improvements over hybrids, while they still excel in thermochemistry, kinetics, transition metal chemistry and the description of strained organic systems.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical ChemistryVrije UniversiteitAmsterdamThe Netherlands
| | - Lucas Visscher
- Theoretical ChemistryVrije UniversiteitAmsterdamThe Netherlands
| |
Collapse
|
48
|
Sylvetsky N, Banerjee A, Alonso M, Martin JML. Performance of Localized Coupled Cluster Methods in a Moderately Strong Correlation Regime: Hückel-Möbius Interconversions in Expanded Porphyrins. J Chem Theory Comput 2020; 16:3641-3653. [PMID: 32338891 PMCID: PMC7304861 DOI: 10.1021/acs.jctc.0c00297] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
![]()
Localized
orbital coupled cluster theory has recently emerged as
a nonempirical alternative to DFT for large systems. Intuitively,
one might expect such methods to perform less well for highly delocalized
systems. In the present work, we apply both canonical CCSD(T) approximations
and a variety of localized approximations to a set of flexible expanded
porphyrins—macrocycles that can switch between Hückel,
figure-eight, and Möbius topologies under external stimuli.
Both minima and isomerization transition states are considered. We
find that Möbius(-like) structures have much stronger static
correlation character than the remaining structures, and that this
causes significant errors in DLPNO-CCSD(T) and even DLPNO-CCSD(T1) approaches, unless TightPNO cutoffs are employed. If sub-kcal
mol–1 accuracy with respect to canonical relative
energies is required even for Möbius-type systems (or other
systems plagued by strong static correlation), then Nagy and Kallay’s
LNO-CCSD(T) method with “tight” settings is the suitable
localized approach. We propose the present POLYPYR21 data set as a
benchmark for localized orbital methods or, more broadly, for the
ability of lower-level methods to handle energetics with strongly
varying degrees of static correlation.
Collapse
Affiliation(s)
- Nitai Sylvetsky
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| | - Ambar Banerjee
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| | - Mercedes Alonso
- Department of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, 1050 Brussels, Belgium
| | - Jan M L Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| |
Collapse
|
49
|
Ma Q, Werner HJ. Scalable Electron Correlation Methods. 7. Local Open-Shell Coupled-Cluster Methods Using Pair Natural Orbitals: PNO-RCCSD and PNO-UCCSD. J Chem Theory Comput 2020; 16:3135-3151. [PMID: 32275428 DOI: 10.1021/acs.jctc.0c00192] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present well-parallelized local implementations of high-spin open-shell coupled cluster methods with single and double excitations (CCSD) using pair natural orbitals (PNOs). The methods are based on the spin-orbital coupled cluster theory using restricted open-shell Hartree-Fock (ROHF) reference functions. Two variants, namely, PNO-UCCSD and PNO-RCCSD are implemented and compared. In PNO-UCCSD, the coupled cluster amplitudes are spin-unrestricted, while in PNO-RCCSD the linear terms are spin-adapted by a spin-projection approach as described in J. Chem. Phys. 1993, 99, 5219-5227. Near linear scaling of the computational cost with the number of correlated electrons is achieved by applying domain and pair approximations. The PNOs are spin-independent and obtained using a semicanonical spin-restricted MP2 approximation with large domains of projected atomic orbitals (PAOs). The pair approximations of our previously described closed-shell PNO-LCCSD method are carefully revised so that they are compatible to the UCCSD theory, and PNO-UCCSD or PNO-RCCSD calculations for closed-shell molecules yield exactly the same results as corresponding spin-free closed-shell PNO-LCCSD calculations. The convergence of the results with respect to the thresholds and options that control the domain and pair approximations is demonstrated. It is found that large domains are required for the single excitations in open-shell calculations in order to obtain converged results. In general, the errors of relative energies caused by the local approximations can be reduced to below 1 kcal mol-1, even for difficult cases. Presently, PNO-RCCSD and PNO-UCCSD calculations for molecules with 100-200 atoms and augmented triple-ζ basis sets can be carried out in a few hours of elapsed time using ∼100 CPU cores. In addition, the program is also capable of performing distinguishable cluster (PNO-RDCSD and PNO-UDCSC) calculations. The present work is a critical step in developing fully local open-shell PNO-RCCSD(T)-F12 methods.
Collapse
Affiliation(s)
- Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
50
|
Haldar S, Dutta AK. A Multilayer Approach to the Equation of Motion Coupled-Cluster Method for the Electron Affinity. J Phys Chem A 2020; 124:3947-3962. [DOI: 10.1021/acs.jpca.0c01793] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Soumi Haldar
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|