1
|
Bruder F, Franzke YJ, Weigend F. Paramagnetic NMR Shielding Tensors Based on Scalar Exact Two-Component and Spin-Orbit Perturbation Theory. J Phys Chem A 2022; 126:5050-5069. [PMID: 35857421 DOI: 10.1021/acs.jpca.2c03579] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The temperature-dependent Fermi-contact and pseudocontact terms are important contributions to the paramagnetic NMR shielding tensor. Herein, we augment the scalar-relativistic (local) exact two-component (X2C) framework with spin-orbit perturbation theory including the screened nuclear spin-orbit correction for the EPR hyperfine coupling and g tensor to compute these temperature-dependent terms. The accuracy of this perturbative ansatz is assessed with the self-consistent spin-orbit two-component and four-component treatments serving as reference. This shows that the Fermi-contact and pseudocontact interaction is sufficiently described for paramagnetic NMR shifts; however, larger deviations are found for the EPR spectra and the principle components of the EPR properties of heavy elements. The impact of the perturbative treatment is further compared to that of the density functional approximation and the basis set. Large-scale calculations are routinely possible with the multipole-accelerated resolution of the identity approximation and the seminumerical exchange approximation, as shown for [CeTi6O3(OiPr)9(salicylate)6].
Collapse
Affiliation(s)
- Florian Bruder
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Florian Weigend
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| |
Collapse
|
2
|
Franzke YJ, Yu JM. Quasi-Relativistic Calculation of EPR g Tensors with Derivatives of the Decoupling Transformation, Gauge-Including Atomic Orbitals, and Magnetic Balance. J Chem Theory Comput 2022; 18:2246-2266. [PMID: 35354319 DOI: 10.1021/acs.jctc.1c01175] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
We present an exact two-component (X2C) ansatz for the EPR g tensor using gauge-including atomic orbitals (GIAOs) and a magnetically balanced basis set expansion. In contrast to previous X2C and four-component relativistic ansätze for the g tensor, this implementation results in a gauge-origin-invariant formalism. Furthermore, the derivatives of the relativistic decoupling matrix are incorporated to form the complete analytical derivative of the X2C Hamiltonian. To reduce the associated computational costs, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The quasi-relativistic X2C and DLU-X2C Hamiltonians accurately reproduce the results of the parent four-component relativistic theory when accounting for two-electron picture-change effects with the modified screened nuclear spin-orbit approximation in the respective one-electron integrals and integral derivatives. According to our benchmark studies, the uncontracted Dyall and segmented-contracted Karlsruhe x2c-type basis sets perform well when compared to large even-tempered basis sets. Moreover, (range-separated) hybrid density functional approximations such as LC-ωPBE and ωB97X-D are needed to match the experimental findings. The impact of the GIAOs depends on the distribution of the spin density, and their use may change the Δg shifts by 10-50% as shown for [(C5Me5)2Y(μ-S)2Mo(μ-S)2Y(C5Me5)2]-. Routine calculations of large molecules are possible with widely available and comparably low-cost hardware as demonstrated for [Pt(C6Cl5)4]- with 3003 basis functions and three spin-(1/2) La(II) and Lu(II) compounds, for which we observe good agreement with the experimental findings.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Jason M Yu
- Department of Chemistry, University of California─Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| |
Collapse
|
3
|
Glasbrenner M, Vogler S, Ochsenfeld C. Gauge-origin dependence in electronic g-tensor calculations. J Chem Phys 2018; 148:214101. [DOI: 10.1063/1.5028454] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michael Glasbrenner
- Chair of Theoretical Chemistry and Center for Integrated Protein Science Munich (CIPSM), Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, 81377 Munich, Germany
| | - Sigurd Vogler
- Chair of Theoretical Chemistry and Center for Integrated Protein Science Munich (CIPSM), Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, 81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry and Center for Integrated Protein Science Munich (CIPSM), Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, 81377 Munich, Germany
| |
Collapse
|
4
|
Vad MS, Pedersen MN, Nørager A, Jensen HJA. Correlated four-component EPR g-tensors for doublet molecules. J Chem Phys 2013; 138:214106. [DOI: 10.1063/1.4804352] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
|
5
|
Verma P, Autschbach J. Variational versus Perturbational Treatment of Spin–Orbit Coupling in Relativistic Density Functional Calculations of Electronic g Factors: Effects from Spin-Polarization and Exact Exchange. J Chem Theory Comput 2013; 9:1052-67. [DOI: 10.1021/ct3009864] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Prakash Verma
- Department of Chemistry, State University of New York
at Buffalo, Buffalo, New York 14260-3000
| | - Jochen Autschbach
- Department of Chemistry, State University of New York
at Buffalo, Buffalo, New York 14260-3000
| |
Collapse
|
6
|
Bonhomme C, Gervais C, Babonneau F, Coelho C, Pourpoint F, Azaïs T, Ashbrook SE, Griffin JM, Yates JR, Mauri F, Pickard CJ. First-principles calculation of NMR parameters using the gauge including projector augmented wave method: a chemist's point of view. Chem Rev 2012; 112:5733-79. [PMID: 23113537 DOI: 10.1021/cr300108a] [Citation(s) in RCA: 326] [Impact Index Per Article: 27.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Affiliation(s)
- Christian Bonhomme
- Laboratoire de Chimie de la Matière Condensée de Paris, Université Pierre et Marie Curie, CNRS UMR, Collège de France, France.
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
7
|
Rokhsana D, Howells AE, Dooley DM, Szilagyi RK. Role of the Tyr-Cys Cross-link to the Active Site Properties of Galactose Oxidase. Inorg Chem 2012; 51:3513-24. [DOI: 10.1021/ic2022769] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dalia Rokhsana
- Department of Chemistry
and Biochemistry, Montana State University, Bozeman, Montana 59717, United
States
| | - Alta E. Howells
- Department of Chemistry
and Biochemistry, Montana State University, Bozeman, Montana 59717, United
States
| | - David M. Dooley
- Department of Chemistry
and Biochemistry, Montana State University, Bozeman, Montana 59717, United
States
| | - Robert K. Szilagyi
- Department of Chemistry
and Biochemistry, Montana State University, Bozeman, Montana 59717, United
States
| |
Collapse
|
8
|
Aquino F, Pritchard B, Autschbach J. Scalar Relativistic Computations and Localized Orbital Analyses of Nuclear Hyperfine Coupling and Paramagnetic NMR Chemical Shifts. J Chem Theory Comput 2012; 8:598-609. [DOI: 10.1021/ct2008507] [Citation(s) in RCA: 73] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Fredy Aquino
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260−3000
| | - Ben Pritchard
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260−3000
| | - Jochen Autschbach
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260−3000
| |
Collapse
|
9
|
Aquino F, Govind N, Autschbach J. Scalar Relativistic Computations of Nuclear Magnetic Shielding and g-Shifts with the Zeroth-Order Regular Approximation and Range-Separated Hybrid Density Functionals. J Chem Theory Comput 2011; 7:3278-92. [DOI: 10.1021/ct200408j] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Fredy Aquino
- Department of Chemistry, State University of New York at Buffalo, Buffalo, New York 14260-3000
| | - Niranjan Govind
- William R. Wiley Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, 902 Battelle Blvd, P.O. Box 999, Mail Stop K8-91 Richland, Washington 99352, United States
| | - Jochen Autschbach
- Department of Chemistry, State University of New York at Buffalo, Buffalo, New York 14260-3000
| |
Collapse
|
10
|
Autschbach J, Patchkovskii S, Pritchard B. Calculation of Hyperfine Tensors and Paramagnetic NMR Shifts Using the Relativistic Zeroth-Order Regular Approximation and Density Functional Theory. J Chem Theory Comput 2011; 7:2175-88. [DOI: 10.1021/ct200143w] [Citation(s) in RCA: 97] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jochen Autschbach
- Department of Chemistry, State University of New York at Buffalo, Buffalo, New York 14260-3000, United States
| | - Serguei Patchkovskii
- National Research Council of Canada, 100 Sussex Drive Ottawa, Ontario K1A 0R6, Canada
| | - Ben Pritchard
- Department of Chemistry, State University of New York at Buffalo, Buffalo, New York 14260-3000, United States
| |
Collapse
|
11
|
Autschbach J, Pritchard B. Calculation of molecular g-tensors using the zeroth-order regular approximation and density functional theory: expectation value versus linear response approaches. Theor Chem Acc 2011. [DOI: 10.1007/s00214-010-0880-x] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
|
12
|
Tatchen J, Kleinschmidt M, Marian CM. Calculating electron paramagnetic resonance g-matrices for triplet state molecules from multireference spin-orbit configuration interaction wave functions. J Chem Phys 2009; 130:154106. [DOI: 10.1063/1.3115965] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
|
13
|
Steimle TC, Wang H, Gengler JJ, Stoll M, Meijer G. The Zeeman effect in the (0,0) band of the A 7Pi-X 7Sigma(+) transition of manganese monohydride, MnH. J Chem Phys 2008; 129:164318. [PMID: 19045277 DOI: 10.1063/1.3000007] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Zeeman tuning of the P(1)(0) line (nu=17 568.35 cm(-1)) of the A (7)Pi-X (7)Sigma(+) (0,0) band of manganese monohydride, MnH, has been investigated. The laser induced fluorescence spectrum of a supersonic molecular beam sample was recorded at a resolution of approximately 40 MHz and with field strengths of up to 362.0 mT. The observed spectrum was successfully fitted using a traditional effective Zeeman Hamiltonian to determine an effective magnetic g-factor for the J=2 level of the F(1)-spin component of the A (7)Pi(v=0) state. Spectral predictions of the P(1)(0) line at field strengths used in magnetic trapping experiments are presented.
Collapse
Affiliation(s)
- Timothy C Steimle
- Department of Chemistry and Biochemistry, Arizona State University, Tempe, Arizona 85287-1604, USA.
| | | | | | | | | |
Collapse
|
14
|
Rinkevicius Z, de Almeida KJ, Oprea CI, Vahtras O, Ågren H, Ruud K. Degenerate Perturbation Theory for Electronic g Tensors: Leading-Order Relativistic Effects. J Chem Theory Comput 2008; 4:1810-28. [DOI: 10.1021/ct800053f] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Zilvinas Rinkevicius
- Department of Theoretical Chemistry, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - Katia Julia de Almeida
- Department of Theoretical Chemistry, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - Cornel I. Oprea
- Department of Theoretical Chemistry, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - Olav Vahtras
- Department of Theoretical Chemistry, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - Hans Ågren
- Department of Theoretical Chemistry, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | - Kenneth Ruud
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø, N-9037 Tromsø, Norway
| |
Collapse
|
15
|
|
16
|
Soncini A. Charge and Spin Currents in Open-Shell Molecules: A Unified Description of NMR and EPR Observables. J Chem Theory Comput 2007; 3:2243-57. [DOI: 10.1021/ct700169h] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Alessandro Soncini
- Department of Chemistry, Laboratory of Quantum Chemistry, Katholieke Universiteit Leuven, Celestijnenlaan 200F, B-3001 Heverlee, Belgium
| |
Collapse
|
17
|
Kacprzak S, Reviakine R, Kaupp M. Understanding the electon paramagnetic resonance parameters of protein-bound glycyl radicals. J Phys Chem B 2007; 111:820-31. [PMID: 17249826 DOI: 10.1021/jp0674571] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The number of enzymes that require a glycyl-based radical for their function is growing. Here, we provide systematic quantum-chemical studies of spin-density distributions, electronic g-tensors, and hyperfine couplings of various models of protein-bound glycyl radicals. Similarly to what is found in a companion paper on N-acetylglycyl, the small g-anisotropy for this delocalized, unsymmetrical system presents appreciable challenges to state-of-the-art computational methodology. This pertains to the quality of structure optimization, as well as to the choice of the spin-orbit Hamiltonian and the gauge origin of the magnetic vector potential. Environmental effects due to hydrogen bonding are complicated and depend in a subtle fashion on the different intramolecular hydrogen bonding for different conformations of the radical. Indeed, the conformation has the largest overall effect on the computed g-tensors (less so on the hyperfine tensors). This is discussed in the context of different g-tensors obtained by recent high-field electron paramagnetic resonance (EPR) measurements for three different enzymes. On the basis of results of a parallel calibration study for N-acetylglycyl, it is suggested that the glycyl radical observed for E. coli anaerobic RNR may have a fully extended conformation, which differs from those of the corresponding radicals in pyruvate formate-lyase or benzylsuccinate synthase.
Collapse
Affiliation(s)
- Sylwia Kacprzak
- Institut für Anorganische Chemie, Universität Würzburg, Am Hubland, D 97074 Würzburg, Germany
| | | | | |
Collapse
|
18
|
Kacprzak S, Reviakine R, Kaupp M. Understanding the EPR Parameters of Glycine-Derived Radicals: The Case of N-Acetylglycyl in the N-Acetylglycine Single-Crystal Environment. J Phys Chem B 2006; 111:811-9. [PMID: 17249825 DOI: 10.1021/jp0660379] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
As a step toward an in-depth understanding of the electron paramagnetic resonance parameters of glycyl radicals in proteins, the hyperfine tensors and, particularly, the g-tensor of N-acetylglcyl in the environment of a single crystal of N-acetylglycine have been studied by systematic state-of-the-art quantum chemical calculations on various suitable model systems. The quantitative computation of the g-tensors for such glycyl-derived radicals is a veritable challenge, mainly because of the very small g-anisotropy combined with a nonsymmetrical, delocalized spin-density distribution and several atoms with comparable spin-orbit contributions to the g-tensors. The choice of gauge origin of the magnetic vector potential, and of approximate spin-orbit operators, both turn out to be more critical than found in previous studies of g-tensors for organic radicals. Environmental effects, included by supermolecular hydrogen-bonded models, were found to be moderate, because of a partial compensation between the influences from intramolecular and intermolecular hydrogen bonds. The largest effects on the g-tensor are caused by the conformation of the radical. The density functional theory methods employed systematically overestimate both the Delta gx and Delta gy components of the g-tensor. This is important for parallel investigations on the protein-glycyl radicals. The 1H alpha and 13C alpha hyperfine couplings depend only slightly on the supermolecular model chosen and appear less sensitive probes of detailed structure and environment.
Collapse
Affiliation(s)
- Sylwia Kacprzak
- Institut für Anorganische Chemie, Universität Würzburg, Am Hubland, D 97074 Würzburg, Germany
| | | | | |
Collapse
|