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
|
Wilhelmer R, Diez M, Krondorfer JK, Hauser AW. Molecular Pseudorotation in Phthalocyanines as a Tool for Magnetic Field Control at the Nanoscale. J Am Chem Soc 2024; 146:14620-14632. [PMID: 38743819 PMCID: PMC11140746 DOI: 10.1021/jacs.4c01915] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2024] [Revised: 04/29/2024] [Accepted: 04/30/2024] [Indexed: 05/16/2024]
Abstract
Metal phthalocyanines, a highly versatile class of aromatic, planar, macrocyclic molecules with a chelated central metal ion, are topical objects of ongoing research and particularly interesting due to their magnetic properties. However, while the current focus lies almost exclusively on spin-Zeeman-related effects, the high symmetry of the molecule and its circular shape suggests the exploitation of light-induced excitation of 2-fold degenerate vibrational states in order to generate, switch, and manipulate magnetic fields at the nanoscale. The underlying mechanism is a molecular pseudorotation that can be triggered by infrared pulses and gives rise to a quantized, small, but controllable magnetic dipole moment. We investigate the optical stimulation of vibrationally induced molecular magnetism and estimate changes in the magnetic shielding constants for confirmation by future experiments.
Collapse
Affiliation(s)
- Raphael Wilhelmer
- Institute of Experimental Physics, Graz University of Technology, Petersgasse 16, A-8010 Graz, Austria
| | - Matthias Diez
- Institute of Experimental Physics, Graz University of Technology, Petersgasse 16, A-8010 Graz, Austria
| | - Johannes K Krondorfer
- Institute of Experimental Physics, Graz University of Technology, Petersgasse 16, A-8010 Graz, Austria
| | - Andreas W Hauser
- Institute of Experimental Physics, Graz University of Technology, Petersgasse 16, A-8010 Graz, Austria
| |
Collapse
|
3
|
Schattenberg C, Wodyński A, Åström H, Sundholm D, Kaupp M, Lehtola S. Revisiting Gauge-Independent Kinetic Energy Densities in Meta-GGAs and Local Hybrid Calculations of Magnetizabilities. J Phys Chem A 2023; 127:10896-10907. [PMID: 38100678 PMCID: PMC10758120 DOI: 10.1021/acs.jpca.3c06244] [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/17/2023] [Revised: 10/27/2023] [Accepted: 11/08/2023] [Indexed: 12/17/2023]
Abstract
In a recent study [J. Chem. Theory Comput. 2021, 17, 1457-1468], some of us examined the accuracy of magnetizabilities calculated with density functionals representing the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA (mGGA), as well as global hybrid (GH) and range-separated (RS) hybrid functionals by assessment against accurate reference values obtained with coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)]. Our study was later extended to local hybrid (LH) functionals by Holzer et al. [J. Chem. Theory Comput. 2021, 17, 2928-2947]; in this work, we examine a larger selection of LH functionals, also including range-separated LH (RSLH) functionals and strong-correlation LH (scLH) functionals. Holzer et al. also studied the importance of the physically correct handling of the magnetic gauge dependence of the kinetic energy density (τ) in mGGA calculations by comparing the Maximoff-Scuseria formulation of τ used in our aforementioned study to the more physical current-density extension derived by Dobson. In this work, we also revisit this comparison with a larger selection of mGGA functionals. We find that the newly tested LH, RSLH, and scLH functionals outperform all of the functionals considered in the previous studies. The various LH functionals afford the seven lowest mean absolute errors while also showing remarkably small standard deviations and mean errors. Most strikingly, the best two functionals are scLHs that also perform remarkably well in cases with significant multiconfigurational character, such as the ozone molecule, which is traditionally excluded from statistical error evaluations due to its large errors with common density functionals.
Collapse
Affiliation(s)
- Caspar
J. Schattenberg
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Artur Wodyński
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Hugo Åström
- Department
of Chemistry, Faculty of Science, University
of Helsinki, P.O. Box 55
(A.I. Virtanens plats 1), University of Helsinki FI-00014, Finland
| | - Dage Sundholm
- Department
of Chemistry, Faculty of Science, University
of Helsinki, P.O. Box 55
(A.I. Virtanens plats 1), University of Helsinki FI-00014, Finland
| | - Martin Kaupp
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Susi Lehtola
- Department
of Chemistry, Faculty of Science, University
of Helsinki, P.O. Box 55
(A.I. Virtanens plats 1), University of Helsinki FI-00014, Finland
- Molecular
Sciences Software Institute, Blacksburg, Virginia 24061, United States
| |
Collapse
|
4
|
Sundholm D, Dimitrova M, Berger RJF. Current density and molecular magnetic properties. Chem Commun (Camb) 2021; 57:12362-12378. [PMID: 34726205 DOI: 10.1039/d1cc03350f] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We give an overview of the molecular response to an external magnetic field perturbing quantum mechanical systems. We present state-of-the-art methods for calculating magnetically-induced current-density susceptibilities. We discuss the essence and properties of current-density susceptibilities and how molecular magnetic properties can be calculated from them. We also review the theory of spin-current densities, how relativity affects current densities and magnetic properties. An overview of the magnetic ring-current criterion for aromaticity is given, which has implications on theoretical and experimental research. The recently reported theory of antiaromaticity and how molecular symmetry affects the magnetic response are discussed and applied to closed-shell paramagnetic molecules. The topology of magnetically induced current densities and its consequences for molecular magnetic properties are also presented with twisted and toroidal molecules as examples.
Collapse
Affiliation(s)
- Dage Sundholm
- Department of Chemistry, Faculty of Science, FI-00014 University of Helsinki, P.O. Box 55, A. I. Virtasen aukio 1, Finland.
| | - Maria Dimitrova
- Department of Chemistry, Faculty of Science, FI-00014 University of Helsinki, P.O. Box 55, A. I. Virtasen aukio 1, Finland. .,Chemistry of Materials, Paris-Lodron University of Salzburg, Jakob-Haringerstr. 2A, A-5020 Salzburg, Austria
| | - Raphael J F Berger
- Chemistry of Materials, Paris-Lodron University of Salzburg, Jakob-Haringerstr. 2A, A-5020 Salzburg, Austria
| |
Collapse
|
5
|
Stauch T, Ganoe B, Wong J, Lee J, Rettig A, Liang J, Li J, Epifanovsky E, Head-Gordon T, Head-Gordon M. Molecular magnetisabilities computed via finite fields: assessing alternatives to MP2 and revisiting magnetic exaltations in aromatic and antiaromatic species. Mol Phys 2021; 119. [DOI: 10.1080/00268976.2021.1990426] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Tim Stauch
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Brad Ganoe
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Jonathan Wong
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Joonho Lee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Adam Rettig
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Jiashu Liang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Jie Li
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | | | - Teresa Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, CA, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| |
Collapse
|
6
|
Stoychev GL, Auer AA, Gauss J, Neese F. DLPNO-MP2 second derivatives for the computation of polarizabilities and NMR shieldings. J Chem Phys 2021; 154:164110. [PMID: 33940835 DOI: 10.1063/5.0047125] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a derivation and efficient implementation of the formally complete analytic second derivatives for the domain-based local pair natural orbital second order Møller-Plesset perturbation theory (MP2) method, applicable to electric or magnetic field-response properties but not yet to harmonic frequencies. We also discuss the occurrence and avoidance of numerical instability issues related to singular linear equation systems and near linear dependences in the projected atomic orbital domains. A series of benchmark calculations on medium-sized systems is performed to assess the effect of the local approximation on calculated nuclear magnetic resonance shieldings and the static dipole polarizabilities. Relative deviations from the resolution of the identity-based MP2 (RI-MP2) reference for both properties are below 0.5% with the default truncation thresholds. For large systems, our implementation achieves quadratic effective scaling, is more efficient than RI-MP2 starting at 280 correlated electrons, and is never more than 5-20 times slower than the equivalent Hartree-Fock property calculation. The largest calculation performed here was on the vancomycin molecule with 176 atoms, 542 correlated electrons, and 4700 basis functions and took 3.3 days on 12 central processing unit cores.
Collapse
Affiliation(s)
- Georgi L Stoychev
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Alexander A Auer
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, 55128 Mainz, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
7
|
Lehtola S, Dimitrova M, Fliegl H, Sundholm D. Benchmarking Magnetizabilities with Recent Density Functionals. J Chem Theory Comput 2021; 17:1457-1468. [PMID: 33599491 PMCID: PMC8023670 DOI: 10.1021/acs.jctc.0c01190] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2020] [Indexed: 12/21/2022]
Abstract
We have assessed the accuracy of the magnetic properties of a set of 51 density functional approximations, including both recently published and already established functionals. The accuracy assessment considers a series of 27 small molecules and is based on comparing the predicted magnetizabilities to literature reference values calculated using coupled-cluster theory with full singles and doubles and perturbative triples [CCSD(T)] employing large basis sets. The most accurate magnetizabilities, defined as the smallest mean absolute error, are obtained with the BHandHLYP functional. Three of the six studied Berkeley functionals and the three range-separated Florida functionals also yield accurate magnetizabilities. Also, some older functionals like CAM-B3LYP, KT1, BHLYP (BHandH), B3LYP, and PBE0 perform rather well. In contrast, unsatisfactory performance is generally obtained with Minnesota functionals, which are therefore not recommended for calculations of magnetically induced current density susceptibilities and related magnetic properties such as magnetizabilities and nuclear magnetic shieldings. We also demonstrate that magnetizabilities can be calculated by numerical integration of magnetizability density; we have implemented this approach as a new feature in the gauge-including magnetically induced current (GIMIC) method. Magnetizabilities can be calculated from magnetically induced current density susceptibilities within this approach even when analytical approaches for magnetizabilities as the second derivative of the energy have not been implemented. The magnetizability density can also be visualized, providing additional information that is not otherwise easily accessible on the spatial origin of magnetizabilities.
Collapse
Affiliation(s)
- Susi Lehtola
- Department
of Chemistry, University of Helsinki, P.O. Box 55, A.I. Virtanens plats
1, FI-00014 University
of Helsinki, Finland
- Molecular
Sciences Software Institute, Blacksburg, Virginia 24061, United States
| | - Maria Dimitrova
- Department
of Chemistry, University of Helsinki, P.O. Box 55, A.I. Virtanens plats
1, FI-00014 University
of Helsinki, Finland
| | - Heike Fliegl
- Institute
of Nanotechnology, KIT, Hermann-von-Helmholtz Platz 1, D-76344 Eggenstein-Leopoldshafen, Germany
| | - Dage Sundholm
- Department
of Chemistry, University of Helsinki, P.O. Box 55, A.I. Virtanens plats
1, FI-00014 University
of Helsinki, Finland
| |
Collapse
|
8
|
Pinski P, Neese F. Analytical gradient for the domain-based local pair natural orbital second order Møller-Plesset perturbation theory method (DLPNO-MP2). J Chem Phys 2019; 150:164102. [DOI: 10.1063/1.5086544] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Peter Pinski
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
9
|
Reimann S, Borgoo A, Austad J, Tellgren EI, Teale AM, Helgaker T, Stopkowicz S. Kohn–Sham energy decomposition for molecules in a magnetic field. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1495849] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Sarah Reimann
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - Alex Borgoo
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - Jon Austad
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - Erik I. Tellgren
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
| | - Andrew M. Teale
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Oslo, Norway
- Centre for Advanced Study at the Norwegian Academy of Science and Letters, Oslo, Norway
- School of Chemistry, University of Nottingham, Nottingham, UK
| | - Trygve Helgaker
- Centre for Advanced Study at the Norwegian Academy of Science and Letters, Oslo, Norway
- School of Chemistry, University of Nottingham, Nottingham, UK
| | - Stella Stopkowicz
- Centre for Advanced Study at the Norwegian Academy of Science and Letters, Oslo, Norway
- Institut für Physikalische Chemie, Johannes Gutenberg-Universität Mainz, Mainz, Germany
| |
Collapse
|
10
|
Field-Theodore TE, Wilson DJD. A coupled cluster study of the magnetisability, rotational g-tensor and quadrupole moment of NF3, PF3 and AsF3. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1349346] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Terri E. Field-Theodore
- Department of Chemistry and Physics, La Trobe Institute for Molecular Science, La Trobe University, Melbourne, Victoria, Australia
| | - David J. D. Wilson
- Department of Chemistry and Physics, La Trobe Institute for Molecular Science, La Trobe University, Melbourne, Victoria, Australia
| |
Collapse
|
11
|
Valiev RR, Fliegl H, Sundholm D. Closed-shell paramagnetic porphyrinoids. Chem Commun (Camb) 2017; 53:9866-9869. [DOI: 10.1039/c7cc05232d] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Magnetizabilities and magnetically induced ring-current strength susceptibilities have been calculated at the Hartree–Fock, density functional theory and second order Møller–Plesset levels for a number of antiaromatic closed-shell carbaporphyrins, carbathiaporphyrins and isophlorins.
Collapse
Affiliation(s)
- Rashid R. Valiev
- Tomsk State University
- Tomsk
- Russian Federation
- Department of Chemistry
- University of Helsinki
| | - Heike Fliegl
- Centre for Theoretical and Computational Chemistry (CTCC)
- Department of Chemistry
- University of Oslo
- 0315 Oslo
- Norway
| | - Dage Sundholm
- Department of Chemistry
- University of Helsinki
- Finland
| |
Collapse
|
12
|
Mandal A, Hunt KLC. Gauge-invariant expectation values of the energy of a molecule in an electromagnetic field. J Chem Phys 2016; 144:044109. [PMID: 26827204 DOI: 10.1063/1.4938564] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we show that the full Hamiltonian for a molecule in an electromagnetic field can be separated into a molecular Hamiltonian and a field Hamiltonian, both with gauge-invariant expectation values. The expectation value of the molecular Hamiltonian gives physically meaningful results for the energy of a molecule in a time-dependent applied field. In contrast, the usual partitioning of the full Hamiltonian into molecular and field terms introduces an arbitrary gauge-dependent potential into the molecular Hamiltonian and leaves a gauge-dependent form of the Hamiltonian for the field. With the usual partitioning of the Hamiltonian, this same problem of gauge dependence arises even in the absence of an applied field, as we show explicitly by considering a gauge transformation from zero applied field and zero external potentials to zero applied field, but non-zero external vector and scalar potentials. We resolve this problem and also remove the gauge dependence from the Hamiltonian for a molecule in a non-zero applied field and from the field Hamiltonian, by repartitioning the full Hamiltonian. It is possible to remove the gauge dependence because the interaction of the molecular charges with the gauge potential cancels identically with a gauge-dependent term in the usual form of the field Hamiltonian. We treat the electromagnetic field classically and treat the molecule quantum mechanically, but nonrelativistically. Our derivation starts from the Lagrangian for a set of charged particles and an electromagnetic field, with the particle coordinates, the vector potential, the scalar potential, and their time derivatives treated as the variables in the Lagrangian. We construct the full Hamiltonian using a Lagrange multiplier method originally suggested by Dirac, partition this Hamiltonian into a molecular term Hm and a field term Hf, and show that both Hm and Hf have gauge-independent expectation values. Any gauge may be chosen for the calculations; but following our partitioning, the expectation values of the molecular Hamiltonian are identical to those obtained directly in the Coulomb gauge. As a corollary of this result, the power absorbed by a molecule from a time-dependent, applied electromagnetic field is equal to the time derivative of the non-adiabatic term in the molecular energy, in any gauge.
Collapse
Affiliation(s)
- Anirban Mandal
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Katharine L C Hunt
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
13
|
Reynolds RD, Shiozaki T. Fully relativistic self-consistent field under a magnetic field. Phys Chem Chem Phys 2015; 17:14280-3. [DOI: 10.1039/c4cp04027a] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Four-component Dirac–Hartree–Fock method with gauge-including atomic orbitals.
Collapse
Affiliation(s)
| | - Toru Shiozaki
- Department of Chemistry
- Northwestern University
- Evanston
- USA
| |
Collapse
|