1
|
Wu H, Hemmingsen L, Sauer SPA. On the geometry dependence of the nuclear magnetic resonance chemical shift of mercury in thiolate complexes: A relativistic density functional theory study. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2024; 62:648-669. [PMID: 38773942 DOI: 10.1002/mrc.5452] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Revised: 04/19/2024] [Accepted: 04/22/2024] [Indexed: 05/24/2024]
Abstract
Thiolate containing mercury(II) complexes of the general formula [Hg(SR)n ]2 - n have been of great interest since the toxicity of mercury was recognized. 199Hg nuclear magnetic resonance spectroscopy (NMR) is a powerful tool for characterization of mercury complexes. In this work, the Hg shielding constants in a series of [Hg(SR)n ]2 - n complexes are therefore investigated computationally with particular emphasis on their geometry dependence. Geometry optimizations and NMR chemical shift calculations are performed at the density functional theory (DFT) level with both the zeroth-order regular approximation (ZORA) and four-component relativistic methods. The four exchange-correlation (XC) functionals PBE0, PBE, B3LYP, and BLYP are used in combination with either Dyall's Gaussian-type (GTO) or Slater-type orbitals (STOs) basis sets. Comparing ZORA and four-component calculations, one observes that the calculated shielding constants for a given molecular geometry have a constant difference of ∼ 1070 ppm. This confirms that ZORA is an acceptable relativistic method to compute NMR chemical shifts. The combinations of four-component/PBE0/v3z and ZORA/PBE0/QZ4P are applied to explore the geometry dependence of the isotropic shielding. For a given coordination number, the distance between mercury and sulfur is the key factor affecting the shielding constant, while changes in bond and dihedral angles and even different side groups have relatively little impact.
Collapse
Affiliation(s)
- Haide Wu
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Lars Hemmingsen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
2
|
Gleeson R, Aggelund PA, Østergaard FC, Schaltz KF, Sauer SPA. Exploring Alternate Methods for the Calculation of High-Level Vibrational Corrections of NMR Spin-Spin Coupling Constants. J Chem Theory Comput 2024. [PMID: 38299500 DOI: 10.1021/acs.jctc.3c01223] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2024]
Abstract
Traditional nuclear magnetic resonance (NMR) calculations typically treat systems with a Born-Oppenheimer-derived electronic wave function that is solved for a fixed nuclear geometry. One can numerically account for this neglected nuclear motion by averaging over property values for all nuclear geometries with a vibrational wave function and adding this expectation value as a correction to an equilibrium geometry property value. Presented are benchmark coupled-cluster singles and doubles (CCSD) vibrational corrections to spin-spin coupling constants (SSCCs) computed at the level of vibrational second-order perturbation theory (VPT2) using the vibrational averaging driver of the CFOUR program. As CCSD calculations of vibrational corrections are very costly, cheaper electronic structure methods are explored via a newly developed Python vibrational averaging program within the Dalton Project. Namely, results obtained with the second-order polarization propagator approximation (SOPPA) and density functional theory (DFT) with the B3LYP and PBE0 exchange-correlation functionals are compared to the benchmark CCSD//CCSD(T) and experimental values. CCSD//CCSD(T) corrections are also combined with literature CC3 equilibrium geometry values to form the highest-order vibrationally corrected values available (i.e., CC3//CCSD(T) + CCSD//CCSD(T)). CCSD//CCSD(T) statistics showed favorable statistics in comparison to experimental values, albeit at an unfavorably high computational cost. A cheaper CCSD//CCSD(T) + B3LYP method showed quite similar mean absolute deviation (MAD) values as CCSD//CCSD(T), concluding that CCSD//CCSD(T) + B3LYP is optimal in terms of cost and accuracy. With reference to experimental values, a vibrational correction was not worth the cost for all of the other methods tested. Finally, deviation statistics showed that CC3//CCSD(T) + CCSD//CCSD(T) vibrational-corrected equilibrium values deteriorated in comparison to CCSD//CCSD(T) attributed to the use of a smaller basis set or lack of solvation effects for the CC3 equilibrium calculations.
Collapse
Affiliation(s)
- Ronan Gleeson
- Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark
| | - Patrick A Aggelund
- Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark
| | | | - Kasper F Schaltz
- Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark
| |
Collapse
|
3
|
Ukhanev SA, Fedorov SV, Krivdin LB. Computational 19 F NMR of trifluoromethyl derivatives of alkenes, pyrimidines, and indenes. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2023; 61:306-317. [PMID: 36740363 DOI: 10.1002/mrc.5335] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Revised: 02/01/2023] [Accepted: 02/02/2023] [Indexed: 06/18/2023]
Abstract
The 19 F NMR chemical shifts of 13 trifluoromethyl derivatives of alkenes, pyrimidines, and indenes were calculated at the DFT level using the BhandHLYP, BHandH, PBE, PBE0, O3LYP, B3LYP, KT2, and KT3 functionals in combination with the pcS-2 basis set. Best result was documented for the BHandHLYP functional: The mean absolute error (MAE) of 0.66 ppm for the scaled values was achieved for the range of about 20 ppm. Solvent, vibrational, and relativistic corrections were found to be rather small, especially when taken in combination, generally demonstrating a slight decrease in the difference between calculated and experimental fluorine chemical shifts. As a measure of the practical importance of these compounds, one should recall that the growing number of life science products that contain trifluoromethyl groups provides a continuing driving force for the development of an effective methodology that enables both regio- and stereoselective introduction of trifluoromethyl groups into both aliphatic and aromatic systems.
Collapse
Affiliation(s)
- Stepan A Ukhanev
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, Irkutsk, 664033, Russia
| | - Sergei V Fedorov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, Irkutsk, 664033, Russia
| | - Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, Irkutsk, 664033, Russia
| |
Collapse
|
4
|
Rzepiela K, Kaminský J, Buczek A, Broda MA, Kupka T. Electron Correlation or Basis Set Quality: How to Obtain Converged and Accurate NMR Shieldings for the Third-Row Elements? MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27238230. [PMID: 36500321 PMCID: PMC9737175 DOI: 10.3390/molecules27238230] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Revised: 11/17/2022] [Accepted: 11/22/2022] [Indexed: 11/29/2022]
Abstract
The quality of theoretical NMR shieldings calculated at the quantum-chemical level depends on various theoretical aspects, of which the basis set type and size are among the most important factors. Nevertheless, not much information is available on the basis set effect on theoretical shieldings of the NMR-active nuclei of the third row. Here, we report on the importance of proper basis set selection to obtain accurate and reliable NMR shielding parameters for nuclei from the third row of the periodic table. All calculations were performed on a set of eleven compounds containing the elements Na, Mg, Al, Si, P, S, or Cl. NMR shielding tensors were calculated using the SCF-HF, DFT-B3LYP, and CCSD(T) methods, combined with the Dunning valence aug-cc-pVXZ, core-valence aug-cc-pCVXZ, Jensen polarized-convergent aug-pcSseg-n and Karlsruhe x2c-Def2 basis set families. We also estimated the complete basis set limit (CBS) values of the NMR parameters. Widely scattered nuclear shieldings were observed for the Dunning polarized-valence basis set, which provides irregular convergence. We show that the use of Dunning core-valence or Jensen basis sets effectively reduces the scatter of theoretical NMR results and leads to their exponential-like convergence to CBS. We also assessed the effect of vibrational, temperature, and relativistic corrections on the predicted shieldings. For systems with single bonds, all corrections are relatively small, amounting to less than 4% of the CCSD(T)/CBS value. Vibrational and temperature corrections were less reliable for H3PO and HSiCH due to the high anharmonicity of the molecules. An abnormally high relativistic correction was observed for phosphorus in PN, reaching ~20% of the CCSD(T)/CBS value, while the correction was less than 7% for other tested molecules.
Collapse
Affiliation(s)
- Kacper Rzepiela
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
| | - Jakub Kaminský
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám. 2, 166 10 Prague, Czech Republic
- Correspondence: (J.K.); (T.K.)
| | - Aneta Buczek
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
| | - Małgorzata A. Broda
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
| | - Teobald Kupka
- Faculty of Chemistry, University of Opole, 48 Oleska Street, 46-052 Opole, Poland
- Correspondence: (J.K.); (T.K.)
| |
Collapse
|
5
|
Krivdin LB. Computational 1 H and 13 C NMR in structural and stereochemical studies. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2022; 60:733-828. [PMID: 35182410 DOI: 10.1002/mrc.5260] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Revised: 02/14/2022] [Accepted: 02/16/2022] [Indexed: 06/14/2023]
Abstract
Present review outlines the advances and perspectives of computational 1 H and 13 C NMR applied to the stereochemical studies of inorganic, organic, and bioorganic compounds, involving in particular natural products, carbohydrates, and carbonium ions. The first part of the review briefly outlines theoretical background of the modern computational methods applied to the calculation of chemical shifts and spin-spin coupling constants at the DFT and the non-empirical levels. The second part of the review deals with the achievements of the computational 1 H and 13 C NMR in the stereochemical investigation of a variety of inorganic, organic, and bioorganic compounds, providing in an abridged form the material partly discussed by the author in a series of parent reviews. Major attention is focused herewith on the publications of the recent years, which were not reviewed elsewhere.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russia
| |
Collapse
|
6
|
Quantum Chemical Approaches to the Calculation of NMR Parameters: From Fundamentals to Recent Advances. MAGNETOCHEMISTRY 2022. [DOI: 10.3390/magnetochemistry8050050] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Quantum chemical methods for the calculation of indirect NMR spin–spin coupling constants and chemical shifts are always in progress. They never stay the same due to permanently developing computational facilities, which open new perspectives and create new challenges every now and then. This review starts from the fundamentals of the nonrelativistic and relativistic theory of nuclear magnetic resonance parameters, and gradually moves towards the discussion of the most popular common and newly developed methodologies for quantum chemical modeling of NMR spectra.
Collapse
|
7
|
Rusakov YY, Rusakova IL. An efficient method for generating property-energy consistent basis sets. New pecJ- n ( n = 1, 2) basis sets for high-quality calculations of indirect nuclear spin-spin coupling constants involving 1H, 13C, 15N, and 19F nuclei. Phys Chem Chem Phys 2021; 23:14925-14939. [PMID: 34223856 DOI: 10.1039/d1cp01984h] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
This paper presents a new method of generating property-energy consistent (PEC) basis sets that can be applied to any arbitrary molecular property. The PEC method generates a basis set that is optimized for the molecular property under interest, providing the least possible total molecular energy. The main algorithm of the PEC approach involves Monte Carlo simulations to generate random exponents in the predetermined range. In this work, the PEC method is introduced in the example of generation of new pecJ-n (n = 1, 2) basis sets suited for high-quality correlated calculations of indirect nuclear spin-spin coupling constants involving the most popular NMR-active nuclei: 1H, 13C, 15N, and 19F.
Collapse
Affiliation(s)
- Yuriy Yu Rusakov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation.
| | - Irina L Rusakova
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation.
| |
Collapse
|
8
|
Glent-Madsen I, Reinholdt A, Bendix J, Sauer SPA. Importance of Relativistic Effects for Carbon as an NMR Reporter Nucleus in Carbide-Bridged [RuCPt] Complexes. Organometallics 2021. [DOI: 10.1021/acs.organomet.1c00079] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Iben Glent-Madsen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Anders Reinholdt
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Jesper Bendix
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | | |
Collapse
|
9
|
Krivdin LB. Computational NMR of Carbohydrates: Theoretical Background, Applications, and Perspectives. Molecules 2021; 26:molecules26092450. [PMID: 33922318 PMCID: PMC8122784 DOI: 10.3390/molecules26092450] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2021] [Revised: 04/19/2021] [Accepted: 04/19/2021] [Indexed: 12/11/2022] Open
Abstract
This review is written amid a marked progress in the calculation of NMR parameters of carbohydrates substantiated by a vast amount of experimental data coming from several laboratories worldwide. By no means are we trying to cover in the present compilation a huge amount of all available data. The main idea of the present review was only to outline general trends and perspectives in this dynamically developing area on the background of a marked progress in theoretical and computational NMR. Presented material is arranged in three basic sections: (1)-a brief theoretical introduction; (2)-applications and perspectives in computational NMR of monosaccharides; and (3)-calculation of NMR chemical shifts and spin-spin coupling constants of di- and polysaccharides.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russia
| |
Collapse
|
10
|
Theibich YA, Sauer SP, Leggio LL, Hedegård ED. Estimating the accuracy of calculated electron paramagnetic resonance hyperfine couplings for a lytic polysaccharide monooxygenase. Comput Struct Biotechnol J 2020; 19:555-567. [PMID: 33510861 PMCID: PMC7807142 DOI: 10.1016/j.csbj.2020.12.014] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2020] [Revised: 12/11/2020] [Accepted: 12/13/2020] [Indexed: 11/07/2022] Open
Abstract
Lytic polysaccharide monooxygenases (LPMOs) are enzymes that bind polysaccharides followed by an (oxidative) disruption of the polysaccharide surface, thereby boosting depolymerization. The binding process between the LPMO catalytic domain and polysaccharide is key to the mechanism and establishing structure-function relationships for this binding is therefore crucial. The hyperfine coupling constants (HFCs) from EPR spectroscopy have proven useful for this purpose. Unfortunately, EPR does not provide direct structural data and therefore the experimental EPR parameters have to be supported with parameters calculated with density functional theory. Yet, calculated HFCs are extremely sensitive to the employed computational setup. Using the LPMO Ls(AA9)A catalytic domain, we here quantify the importance of several choices in the computational setup, ranging from the use of specialized basis, the underlying structures, and the employed exchange-correlation functional. We show that specialized basis sets are an absolute necessity, and also that care has to be taken in the optimization of the underlying structure: only by allowing large parts of the protein around the active site to structurally relax could we obtain results that uniformly reproduced experimental trends. We compare our results to previously published X-ray structures and experimental HFCs for Ls(AA9)A as well as to recent experimental/theoretical results for another (AA10) family of LPMOs.
Collapse
Affiliation(s)
- Yusuf A. Theibich
- Department of Chemistry, University of University, Copenhagen, Denmark
| | | | - Leila Lo Leggio
- Department of Chemistry, University of University, Copenhagen, Denmark
| | - Erik D. Hedegård
- Division of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
11
|
Mack F, Schattenberg CJ, Kaupp M, Weigend F. Nuclear Spin–Spin Couplings: Efficient Evaluation of Exact Exchange and Extension to Local Hybrid Functionals. J Phys Chem A 2020; 124:8529-8539. [DOI: 10.1021/acs.jpca.0c06897] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Fabian Mack
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Caspar J. Schattenberg
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Florian Weigend
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| |
Collapse
|
12
|
Payard PA, Perego LA, Grimaud L, Ciofini I. A DFT Protocol for the Prediction of 31P NMR Chemical Shifts of Phosphine Ligands in First-Row Transition-Metal Complexes. Organometallics 2020. [DOI: 10.1021/acs.organomet.0c00309] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pierre-Adrien Payard
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Luca Alessandro Perego
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
- PSL University, Institute of Chemistry for Health and Life Sciences, I-CLeHS, CNRS-Chimie ParisTech, 11 rue P. et M. Curie, F-75005 Paris 05, France
| | - Laurence Grimaud
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Ilaria Ciofini
- PSL University, Institute of Chemistry for Health and Life Sciences, I-CLeHS, CNRS-Chimie ParisTech, 11 rue P. et M. Curie, F-75005 Paris 05, France
| |
Collapse
|
13
|
Rusakov YY, Rusakova IL, Semenov VA, Krivdin LB. A New Basis Set for the Calculation of 13C NMR Chemical Shifts within a Non-empirical Correlated Framework. J Phys Chem A 2020; 124:7322-7330. [DOI: 10.1021/acs.jpca.0c06038] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Yuriy Yu. Rusakov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| | - Irina L. Rusakova
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| | - Valentin A. Semenov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| | - Leonid B. Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| |
Collapse
|
14
|
Semenov VA, Samultsev DO, Krivdin LB. The 1 H and 13 C NMR chemical shifts of Strychnos alkaloids revisited at the DFT level. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2020; 58:532-539. [PMID: 31663170 DOI: 10.1002/mrc.4948] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/02/2019] [Revised: 09/25/2019] [Accepted: 09/27/2019] [Indexed: 06/10/2023]
Abstract
The density functional theory calculation of 1 H and 13 C NMR chemical shifts in a series of ten 10 classically known Strychnos alkaloids with a strychnine skeleton was performed at the PBE0/pcSseg-2//pcseg-2 level. It was found that calculated 1 H and 13 C NMR chemical shifts provided a markedly good correlation with experiment characterized by a mean absolute error of 0.08 ppm in the range of 7 ppm for protons and 1.67 ppm in the range of 150 ppm for carbons, so that a mean absolute percentage error was as small as ~1% in both cases.
Collapse
Affiliation(s)
- Valentin A Semenov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russia
| | - Dmitry O Samultsev
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russia
| | - Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russia
| |
Collapse
|
15
|
Krivdin LB. Computational 1 H NMR: Part 1. Theoretical background. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2019; 57:897-914. [PMID: 30963636 DOI: 10.1002/mrc.4873] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2019] [Revised: 03/31/2019] [Accepted: 04/01/2019] [Indexed: 06/09/2023]
Abstract
This is the first one of the three closely interrelated reviews to be published in Magnetic Resonance in Chemistry dealing with accordingly theoretical background, chemical applications, and biochemical studies of and by means of computational 1 H NMR. Presented in the first part of the review is a general outline of the modern theoretical methods and accuracy factors of computational 1 H NMR involving locally dense basis set schemes, solvent effects, vibrational corrections, and relativistic effects performed at the density functional theory and/or nonempirical levels. This review is dedicated to Prof. Stephan Sauer in view of his invaluable contribution to the field of computational nuclear magnetic resonance.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russia
- Angarsk State Technical University, Angarsk, Russia
| |
Collapse
|
16
|
Rusakov YY, Rusakova IL, Fedorov SV, Gray GA, Krivdin LB. Stereochemical Dependences of 31P–13C Spin–Spin Coupling Constants of Heterocyclic Phosphines. J Phys Chem A 2019; 123:6298-6303. [DOI: 10.1021/acs.jpca.9b05385] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Yuriy Y. Rusakov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| | - Irina L. Rusakova
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| | - Sergei V. Fedorov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| | - George A. Gray
- Agilent Technologies, 5301 Stevens Creek Blvd, Santa Clara, California 95051, United States
| | - Leonid B. Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| |
Collapse
|
17
|
Krivdin LB. Computational protocols for calculating 13C NMR chemical shifts. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2019; 112-113:103-156. [PMID: 31481156 DOI: 10.1016/j.pnmrs.2019.05.004] [Citation(s) in RCA: 58] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2019] [Revised: 05/19/2019] [Accepted: 05/20/2019] [Indexed: 06/10/2023]
Abstract
The most recent results dealing with the computation of 13C NMR chemical shifts in chemistry (small molecules, saturated, unsaturated and aromatic compounds, heterocycles, functional derivatives, coordination complexes, carbocations, and natural products) are reviewed, paying special attention to theoretical background and accuracy, the latter involving solvent effects, vibrational corrections, and relativistic effects.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russia.
| |
Collapse
|
18
|
Lino JBDR, Ramalho TC. Exploring Through-Space Spin-Spin Couplings for Quantum Information Processing: Facing the Challenge of Coherence Time and Control Quantum States. J Phys Chem A 2019; 123:1372-1379. [PMID: 30673241 DOI: 10.1021/acs.jpca.8b09425] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Nuclear magnetic resonance (NMR) is a powerful tool for studying quantum information processing (QIP). Recently quantum technologies have been proposed to overcome the challenges in large-scale NMR QIP. Furthermore, computational chemistry can promote its improvement. Nuclear spins-1/2 are natural qubits and have been used in most NMR quantum computation experiments. However, molecules that enable many qubits NMR QIP implementations should meet some requirements regarding their spectroscopic properties. Exceptionally large through-space (TS) P-P spin-spin coupling constants (SSCC or J) observed in 1,8-diphosphanaphthalenes (PPN) and in naphtho[1,8- cd]-1,2-dithiole phenylphosphines (NTP) were proposed and investigated to provide more accurate control within large-scale NMR QIP. Spectroscopic properties of PPN and NTP derivatives were explored by theoretical strategies using locally dense basis sets (LDBS). 31P chemical shifts (δ) calculated at the B3LYP/aug-cc-pVTZ-J level and TS P-P SSCCs at the PBE1PBE/pcJ-2 (LDBS-1) level are very close to the experimental data for the PPN molecule. Differently, for the NTP dimer, PBE1PBE/pcJ-2 (LDBS-2) predicts more accurate 31P δ, whereas PBE1PBE/Def2-TZVP (LDBS-1) forecasts more accurate TS P-P SSCCs. From our results, PPNo-F, PPNo-ethyl, and PPNo-NH2 were the best candidates for NMR QIP, in which the large TS SSCCS could face the need of long-time quantum gates implementations. Therefore, it could overcome natural limitations concerning the development of large-scale NMR.
Collapse
Affiliation(s)
| | - Teodorico Castro Ramalho
- Chemistry Department , Federal University of Lavras , 37200-000 Lavras , MG Brazil.,Center for Basic and Applied Research, Faculty of Informatics and Management , University Hradec Kralove , 50003 Hradec Kralove , Czech Republic
| |
Collapse
|
19
|
Rusakov YY, Rusakova IL, Krivdin LB. Relativistic heavy atom effect on the 31 P NMR parameters of phosphine chalcogenides. Part 1. Chemical shifts. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2018; 56:1061-1073. [PMID: 29775489 DOI: 10.1002/mrc.4752] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2018] [Revised: 04/27/2018] [Accepted: 05/08/2018] [Indexed: 06/08/2023]
Abstract
Four-component density functional theory calculations of 31 P NMR chemical shifts have been performed for the representative series of 56 phosphine chalcogenides in order to investigate an influence of different functional groups on the heavy atom relativistic effect on the NMR chemical shifts of light phosphorous atoms (Heavy Atom on Light Atom effect). The validity of the 4-component density functional theory approach used for the wide-scale calculations of the phosphorous chemical shifts in a wide series of phosphine chalcogenides has been confirmed on a small series of 5 representative compounds with the aid of high-quality coupled cluster singles and doubles calculations taking into account solvent, vibrational, and the relativistic corrections in comparison with the experiment.
Collapse
Affiliation(s)
- Yury Yu Rusakov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| | - Irina L Rusakova
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| | - Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| |
Collapse
|
20
|
Krivdin LB. Theoretical calculations of carbon-hydrogen spin-spin coupling constants. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2018; 108:17-73. [PMID: 30538048 DOI: 10.1016/j.pnmrs.2018.10.002] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/08/2018] [Revised: 10/09/2018] [Accepted: 10/11/2018] [Indexed: 06/09/2023]
Abstract
Structural applications of theoretical calculations of carbon-hydrogen spin-spin coupling constants are reviewed covering papers published mainly during the last 10-15 years with a special emphasis on the most notable studies of hybridization, substitution and stereoelectronic effects together with the investigation of hydrogen bonding and intermolecular interactions. The wide scope of different applications of calculated carbon-hydrogen couplings in the structural elucidation of particular classes of organic and bioorganic molecules is reviewed, concentrating mainly on saturated, unsaturated, aromatic and heteroaromatic compounds and their functional derivatives, as well as on natural compounds and carbohydrates. The review is dedicated to Professor Emeritus Michael Barfield in view of his invaluable pioneering contribution to this field.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A.E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russia.
| |
Collapse
|
21
|
Semenov VA, Samultsev DO, Krivdin LB. GIAO-DFT calculation of 15 N NMR chemical shifts of Schiff bases: Accuracy factors and protonation effects. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2018; 56:727-739. [PMID: 29427330 DOI: 10.1002/mrc.4721] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2018] [Revised: 01/23/2018] [Accepted: 01/25/2018] [Indexed: 06/08/2023]
Abstract
15 N NMR chemical shifts in the representative series of Schiff bases together with their protonated forms have been calculated at the density functional theory level in comparison with available experiment. A number of functionals and basis sets have been tested in terms of a better agreement with experiment. Complimentary to gas phase results, 2 solvation models, namely, a classical Tomasi's polarizable continuum model (PCM) and that in combination with an explicit inclusion of one molecule of solvent into calculation space to form supermolecule 1:1 (SM + PCM), were examined. Best results are achieved with PCM and SM + PCM models resulting in mean absolute errors of calculated 15 N NMR chemical shifts in the whole series of neutral and protonated Schiff bases of accordingly 5.2 and 5.8 ppm as compared with 15.2 ppm in gas phase for the range of about 200 ppm. Noticeable protonation effects (exceeding 100 ppm) in protonated Schiff bases are rationalized in terms of a general natural bond orbital approach.
Collapse
Affiliation(s)
- Valentin A Semenov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, Irkutsk, 664033, Russia
| | - Dmitry O Samultsev
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, Irkutsk, 664033, Russia
| | - Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, Irkutsk, 664033, Russia
| |
Collapse
|
22
|
Semenov VA, Samultsev DO, Krivdin LB. Substitution effects in the 15 N NMR chemical shifts of heterocyclic azines evaluated at the GIAO-DFT level. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2018; 56:767-774. [PMID: 29504638 DOI: 10.1002/mrc.4731] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2018] [Revised: 02/18/2018] [Accepted: 02/20/2018] [Indexed: 06/08/2023]
Abstract
A systematic study of the accuracy factors for the computation of 15 N NMR chemical shifts in comparison with available experiment in the series of 72 diverse heterocyclic azines substituted with a classical series of substituents (CH3 , F, Cl, Br, NH2 , OCH3 , SCH3 , COCH3 , CONH2 , COOH, and CN) providing marked electronic σ- and π-electronic effects and strongly affecting 15 N NMR chemical shifts is performed. The best computational scheme for heterocyclic azines at the DFT level was found to be KT3/pcS-3//pc-2 (IEF-PCM). A vast amount of unknown 15 N NMR chemical shifts was predicted using the best computational protocol for substituted heterocyclic azines, especially for trizine, tetrazine, and pentazine where experimental 15 N NMR chemical shifts are almost totally unknown throughout the series. It was found that substitution effects in the classical series of substituents providing typical σ- and π-electronic effects followed the expected trends, as derived from the correlations of experimental and calculated 15 N NMR chemical shifts with Swain-Lupton's F and R constants.
Collapse
Affiliation(s)
- Valentin A Semenov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| | - Dmitry O Samultsev
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| | - Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| |
Collapse
|
23
|
Rusakov YY, Rusakova IL, Semenov VA, Samultsev DO, Fedorov SV, Krivdin LB. Calculation of 15N and 31P NMR Chemical Shifts of Azoles, Phospholes, and Phosphazoles: A Gateway to Higher Accuracy at Less Computational Cost. J Phys Chem A 2018; 122:6746-6759. [DOI: 10.1021/acs.jpca.8b05161] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yury Yu. Rusakov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky Street 1, 664033 Irkutsk, Russia
| | - Irina L. Rusakova
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky Street 1, 664033 Irkutsk, Russia
| | - Valentin A. Semenov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky Street 1, 664033 Irkutsk, Russia
| | - Dmitry O. Samultsev
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky Street 1, 664033 Irkutsk, Russia
| | - Sergei V. Fedorov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky Street 1, 664033 Irkutsk, Russia
| | - Leonid B. Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky Street 1, 664033 Irkutsk, Russia
| |
Collapse
|
24
|
Krivdin LB. Carbon-carbon spin-spin coupling constants: Practical applications of theoretical calculations. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2018; 105:54-99. [PMID: 29548367 DOI: 10.1016/j.pnmrs.2018.03.001] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/14/2017] [Revised: 02/27/2018] [Accepted: 03/01/2018] [Indexed: 06/08/2023]
Abstract
Practical applications of theoretical calculations of carbon-carbon spin-spin coupling constants in particular classes of organic and bioorganic molecules are reviewed, concentrating mainly on saturated, unsaturated, aromatic and heteroaromatic compounds and their functional derivatives as well as on carbohydrates and natural compounds.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A.E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russia.
| |
Collapse
|
25
|
Relativistic effects in the NMR spectra of compounds containing heavy chalcogens. MENDELEEV COMMUNICATIONS 2018. [DOI: 10.1016/j.mencom.2018.01.001] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
26
|
Samultsev DO, Semenov VA, Krivdin LB. On the accuracy factors and computational cost of the GIAO-DFT calculation of 15 N NMR chemical shifts of amides. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2017; 55:1015-1021. [PMID: 28600816 DOI: 10.1002/mrc.4625] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/17/2017] [Revised: 06/04/2017] [Accepted: 06/06/2017] [Indexed: 06/07/2023]
Abstract
The main factors affecting the accuracy and computational cost of Gauge-independent Atomic Orbital-density functional theory (GIAO-DFT) calculation of 15 N NMR chemical shifts in the benchmark series of 16 amides are considered. Among those are the choice of the DFT functional and basis set, solvent effects, internal reference conversion factor and applicability of the locally dense basis set (LDBS) scheme. Solvent effects are treated within the polarizable continuum model (PCM) scheme as well as at supermolecular level with solvent molecules considered in explicit way. The best result is found for Keal and Tozer's KT3 functional used in combination with Jensen's pcS-3 basis set with taking into account solvent effects within the polarizable continuum model. The proposed LDBS scheme implies pcS-3 on nitrogen and pc-2 elsewhere in the molecule. The resulting mean average error for the calculated 15 N NMR chemical shifts is about 6 ppm. The application of the LDBS approach tested in a series of 16 amides results in a dramatic decrease in computational cost (more than an order of magnitude in time scale) with insignificant loss of accuracy.
Collapse
Affiliation(s)
- Dmitry O Samultsev
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, 664033, Irkutsk, Russia
| | - Valentin A Semenov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, 664033, Irkutsk, Russia
| | - Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, 664033, Irkutsk, Russia
| |
Collapse
|
27
|
Krivdin LB. Calculation of 15N NMR chemical shifts: Recent advances and perspectives. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2017; 102-103:98-119. [PMID: 29157495 DOI: 10.1016/j.pnmrs.2017.08.001] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/20/2017] [Revised: 08/21/2017] [Accepted: 08/21/2017] [Indexed: 06/07/2023]
Abstract
Recent advances in computation of 15N NMR chemical shifts are reviewed, concentrating mainly on practical aspects of computational protocols and accuracy factors. The review includes the discussion of the level of theory, the choice of density functionals and basis sets together with taking into account solvent effects, rovibrational corrections and relativistic effects. Computational aspects of 15N NMR are illustrated for the series of neutral and protonated open-chain nitrogen-containing compounds and nitrogen heterocycles, coordination and intermolecular complexes.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russia.
| |
Collapse
|
28
|
Rusakov YY, Rusakova IL, Krivdin LB. On the significant relativistic heavy atom effect on 13C NMR chemical shifts of β- and γ-carbons in seleno- and telluroketones. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1350292] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Yuriy Yu. Rusakov
- A.E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russian Federation
| | - Irina L. Rusakova
- A.E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russian Federation
| | - Leonid B. Krivdin
- A.E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russian Federation
| |
Collapse
|
29
|
Fedorov SV, Rusakov YY, Krivdin LB. Quantum-chemical calculations of NMR chemical shifts of organic molecules: XV. Relativistic calculations of 29Si NMR chemical shifts of silanes. RUSSIAN JOURNAL OF ORGANIC CHEMISTRY 2017. [DOI: 10.1134/s1070428017050013] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
30
|
Rusakova IL, Rusakov YY, Krivdin LB. Calculation of 125Te NMR Chemical Shifts at the Full Four-Component Relativistic Level with Taking into Account Solvent and Vibrational Corrections: A Gateway to Better Agreement with Experiment. J Phys Chem A 2017; 121:4793-4803. [DOI: 10.1021/acs.jpca.7b03198] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Irina L. Rusakova
- A.E. Favorsky Irkutsk
Institute
of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| | - Yuriy Yu. Rusakov
- A.E. Favorsky Irkutsk
Institute
of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| | - Leonid B. Krivdin
- A.E. Favorsky Irkutsk
Institute
of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| |
Collapse
|
31
|
Rusakova IL, Rusakov YY, Krivdin LB. First example of the correlated calculation of the one-bond tellurium-carbon spin-spin coupling constants: Relativistic effects, vibrational corrections, and solvent effects. J Comput Chem 2016; 37:1367-72. [DOI: 10.1002/jcc.24345] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2015] [Revised: 02/01/2016] [Accepted: 02/02/2016] [Indexed: 11/11/2022]
Affiliation(s)
- Irina L. Rusakova
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences; Irkutsk 664033 Russian Federation
| | - Yury Yu. Rusakov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences; Irkutsk 664033 Russian Federation
| | - Leonid B. Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences; Irkutsk 664033 Russian Federation
| |
Collapse
|
32
|
Del Bene JE, Alkorta I, Elguero J. A Systematic Comparison of Second-Order Polarization Propagator Approximation (SOPPA) and Equation-of-Motion Coupled Cluster Singles and Doubles (EOM-CCSD) Spin-Spin Coupling Constants for Selected Singly Bonded Molecules, and the Hydrides NH3, H2O, and HF and Their Protonated and Deprotonated Ions and Hydrogen-Bonded Complexes. J Chem Theory Comput 2015; 4:967-73. [PMID: 26621237 DOI: 10.1021/ct800111j] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Second-order polarization propagator approximation (SOPPA) and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) methods have been employed for the calculation of one-bond spin-spin coupling constants in series of small molecules and ions, and of one- and two-bond coupling constants across X-H···Y hydrogen bonds. For isolated molecules, one-bond SOPPA coupling constants (1)J(X-Y) involving (13)C, (15)N, (17)O, and (19)F have larger absolute values than corresponding EOM-CCSD coupling constants, with the EOM-CCSD values being in significantly better agreement with available experimental data. The difference between SOPPA and EOM-CCSD tends to increase as the number of nonbonding electrons on the coupled atoms increases, and the SOPPA values for O-F coupling are significantly in error. Similarly, the absolute values of SOPPA one-bond coupling constants (1)J(X-H) for the hydrides NH3, H2O, and FH and their protonated and deprotonated ions are greater than EOM-CCSD values, with the largest differences occurring for F-H coupling. One- and two-bond coupling constants (1)J(X-H), (1h)J(H-Y), and (2h)J(X-Y) across X-H···Y hydrogen bonds in neutral, protonated, and deprotonated complexes formed from the hydrides are similar at SOPPA and EOM-CCSD, with the largest differences again found for (1)J(F-H) in complexes with F-H as the proton donor, and (2h)J(F-F) for (FHF)(-). The signs of (1)J(X-H), (1h)J(H-Y), and (2h)J(X-Y) are the same at both levels of theory, as is their variation across the proton-transfer coordinate in F-H···NH3. SOPPA would appear to provide a reliable and more cost-effective alternative approach for computing coupling constants across hydrogen bonds, although couplings involving F may be problematic.
Collapse
Affiliation(s)
- Janet E Del Bene
- Department of Chemistry, Youngstown State University, Youngstown, Ohio 44555, and Instituto de Química Médica, CSIC, Juan de la Cierva, 3, E-28006 Madrid, Spain
| | - Ibon Alkorta
- Department of Chemistry, Youngstown State University, Youngstown, Ohio 44555, and Instituto de Química Médica, CSIC, Juan de la Cierva, 3, E-28006 Madrid, Spain
| | - José Elguero
- Department of Chemistry, Youngstown State University, Youngstown, Ohio 44555, and Instituto de Química Médica, CSIC, Juan de la Cierva, 3, E-28006 Madrid, Spain
| |
Collapse
|
33
|
Semenov VA, Samultsev DO, Rulev AY, Krivdin LB. Theoretical and experimental (15)N NMR study of enamine-imine tautomerism of 4-trifluoromethyl[b]benzo-1,4-diazepine system. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2015; 53:1031-1034. [PMID: 26290420 DOI: 10.1002/mrc.4296] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/09/2015] [Revised: 06/16/2015] [Accepted: 06/17/2015] [Indexed: 06/04/2023]
Abstract
The tautomeric structure of 4-trifluoromethyl[b]benzo-1,4-diazepine system in solution has been evaluated by means of the calculation of (15)N NMR chemical shifts of individual tautomers in comparison with the averaged experimental shifts to show that the enamine-imine equilibrium is entirely shifted toward the imine form. The adequacy of the theoretical level used for the computation of (15)N NMR chemical shifts in this case has been verified based on the benchmark calculations in the series of the push-pull and captodative enamines together with related azomethynes, which demonstrated a good to excellent agreement with experiment.
Collapse
Affiliation(s)
- Valentin A Semenov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| | - Dmitry O Samultsev
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| | - Alexander Yu Rulev
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| | - Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| |
Collapse
|
34
|
On novel magnetic probe for fullerene characterization: Theoretical studies on NMR parameters of free and confined in fullerenes HD and H 2 molecules. J Mol Graph Model 2015; 62:26-37. [DOI: 10.1016/j.jmgm.2015.08.009] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2015] [Revised: 08/05/2015] [Accepted: 08/24/2015] [Indexed: 11/22/2022]
|
35
|
Rusakov YY, Rusakova IL, Krivdin LB. MP2 calculation of (77) Se NMR chemical shifts taking into account relativistic corrections. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2015; 53:485-492. [PMID: 25998325 DOI: 10.1002/mrc.4226] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/17/2014] [Revised: 01/26/2015] [Accepted: 01/29/2015] [Indexed: 06/04/2023]
Abstract
The main factors affecting the accuracy and computational cost of the Second-order Möller-Plesset perturbation theory (MP2) calculation of (77) Se NMR chemical shifts (methods and basis sets, relativistic corrections, and solvent effects) are addressed with a special emphasis on relativistic effects. For the latter, paramagnetic contribution (390-466 ppm) dominates over diamagnetic term (192-198 ppm) resulting in a total shielding relativistic correction of about 230-260 ppm (some 15% of the total values of selenium absolute shielding constants). Diamagnetic term is practically constant, while paramagnetic contribution spans over 70-80 ppm. In the (77) Se NMR chemical shifts scale, relativistic corrections are about 20-30 ppm (some 5% of the total values of selenium chemical shifts). Solvent effects evaluated within the polarizable continuum solvation model are of the same order of magnitude as relativistic corrections (about 5%). For the practical calculations of (77) Se NMR chemical shifts of the medium-sized organoselenium compounds, the most efficient computational protocols employing relativistic Dyall's basis sets and taking into account relativistic and solvent corrections are suggested. The best result is characterized by a mean absolute error of 17 ppm for the span of (77) Se NMR chemical shifts reaching 2500 ppm resulting in a mean absolute percentage error of 0.7%.
Collapse
Affiliation(s)
- Yury Yu Rusakov
- Siberian Branch of the Russian Academy of Sciences, A. E. Favorsky Irkutsk Institute of Chemistry, Favorsky St 1, Irkutsk, 664033, Russia
| | - Irina L Rusakova
- Siberian Branch of the Russian Academy of Sciences, A. E. Favorsky Irkutsk Institute of Chemistry, Favorsky St 1, Irkutsk, 664033, Russia
| | - Leonid B Krivdin
- Siberian Branch of the Russian Academy of Sciences, A. E. Favorsky Irkutsk Institute of Chemistry, Favorsky St 1, Irkutsk, 664033, Russia
| |
Collapse
|
36
|
Fedorov SV, Rusakov YY, Krivdin LB. Towards the versatile DFT and MP2 computational schemes for 31P NMR chemical shifts taking into account relativistic corrections. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2014; 52:699-710. [PMID: 25155415 DOI: 10.1002/mrc.4122] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/07/2014] [Revised: 07/09/2014] [Accepted: 07/11/2014] [Indexed: 06/03/2023]
Abstract
The main factors affecting the accuracy and computational cost of the calculation of (31)P NMR chemical shifts in the representative series of organophosphorous compounds are examined at the density functional theory (DFT) and second-order Møller-Plesset perturbation theory (MP2) levels. At the DFT level, the best functionals for the calculation of (31)P NMR chemical shifts are those of Keal and Tozer, KT2 and KT3. Both at the DFT and MP2 levels, the most reliable basis sets are those of Jensen, pcS-2 or larger, and those of Pople, 6-311G(d,p) or larger. The reliable basis sets of Dunning's family are those of at least penta-zeta quality that precludes their practical consideration. An encouraging finding is that basically, the locally dense basis set approach resulting in a dramatic decrease in computational cost is justified in the calculation of (31)P NMR chemical shifts within the 1-2-ppm error. Relativistic corrections to (31)P NMR absolute shielding constants are of major importance reaching about 20-30 ppm (ca 7%) improving (not worsening!) the agreement of calculation with experiment. Further better agreement with the experiment by 1-2 ppm can be obtained by taking into account solvent effects within the integral equation formalism polarizable continuum model solvation scheme. We recommend the GIAO-DFT-KT2/pcS-3//pcS-2 scheme with relativistic corrections and solvent effects taken into account as the most versatile computational scheme for the calculation of (31)P NMR chemical shifts characterized by a mean absolute error of ca 9 ppm in the range of 550 ppm.
Collapse
Affiliation(s)
- Sergey V Fedorov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| | | | | |
Collapse
|
37
|
Rusakova IL, Rusakov YY, Krivdin LB. Relativistic effects in the one-bond spin-spin coupling constants involving selenium. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2014; 52:500-510. [PMID: 25043341 DOI: 10.1002/mrc.4112] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/15/2014] [Revised: 06/29/2014] [Accepted: 07/01/2014] [Indexed: 06/03/2023]
Abstract
One-bond spin-spin coupling constants involving selenium of seven different types, (1) J(Se,X), X = (1) H, (13) C, (15) N, (19) F, (29) Si, (31) P, and (77) Se, were calculated in the series of 14 representative compounds at the SOPPA(CCSD) level taking into account relativistic corrections evaluated both at the RPA and DFT levels of theory in comparison with experiment. Relativistic corrections were found to play a major role in the calculation of (1) J(Se,X) reaching as much as almost 170% of the total value of (1) J(Se,Se) and up to 60-70% for the rest of (1) J(Se,X). Scalar relativistic effects (Darwin and mass-velocity corrections) by far dominate over spin-orbit coupling in the total relativistic effects for all (1) J(Se,X). Taking into account relativistic corrections at both random phase approximation and density functional theory levels essentially improves the agreement of theoretical results with experiment. The most 'relativistic' (1) J(Se,Se) demonstrates a marked Karplus-type dihedral angle dependence with respect to the mutual orientation of the selenium lone pairs providing a powerful tool for stereochemical analysis of selenoorganic compounds.
Collapse
Affiliation(s)
- Irina L Rusakova
- Siberian Branch of the Russian Academy of Sciences, A. E. Favorsky Irkutsk Institute of Chemistry, Favorsky St. 1, 664033, Irkutsk, Russia
| | | | | |
Collapse
|
38
|
Rusakova IL, Rusakov YY, Krivdin LB. Nonempirical calculations of the one-bond (29)Si-(13)C spin-spin coupling constants taking into account relativistic and solvent corrections. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2014; 52:413-421. [PMID: 24796525 DOI: 10.1002/mrc.4080] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/11/2014] [Revised: 04/08/2014] [Accepted: 04/09/2014] [Indexed: 06/03/2023]
Abstract
The computational study of the one-bond (29)Si-(13)C spin-spin coupling constants has been performed at the second-order polarization propagator approximation (SOPPA) level in the series of 60 diverse silanes with a special focus on the main factors affecting the accuracy of the calculation including the level of theory, the quality of the basis set, and the contribution of solvent and relativistic effects. Among three SOPPA-based methods, SOPPA(MP2), SOPPA(CC2), and SOPPA(CCSD), the best result was achieved with SOPPA(CCSD) when used in combination with Sauer's basis set aug-cc-pVTZ-J characterized by the mean absolute error of calculated coupling constants against the experiment of ca 2 Hz in the range of ca 200 Hz. The SOPPA(CCSD)/aug-cc-pVTZ-J method is recommended as the most accurate and effective computational scheme for the calculation of (1)J(Si,C). The slightly less accurate but essentially more economical SOPPA(MP2)/aug-cc-pVTZ-J and/or SOPPA(CC2)/aug-cc-pVTZ-J methods are recommended for larger molecular systems. It was shown that solvent and relativistic corrections do not play a major role in the computation of the total values of (1)J(Si,C); however, taking them into account noticeably improves agreement with the experiment. The rovibrational corrections are estimated to be of about 1 Hz or 1-1.5% of the total value of (1)J(Si,C).
Collapse
Affiliation(s)
- Irina L Rusakova
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St 1, 664033, Irkutsk, Russia
| | | | | |
Collapse
|
39
|
Rusakov YY, Rusakova IL, Krivdin LB. Full four-component relativistic calculations of the one-bond 77Se-13C spin-spin coupling constants in the series of selenium heterocycles and their parent open-chain selenides. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2014; 52:214-221. [PMID: 24549877 DOI: 10.1002/mrc.4053] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/25/2013] [Revised: 01/19/2014] [Accepted: 01/23/2014] [Indexed: 06/03/2023]
Abstract
Four-component relativistic calculations of (77)Se-(13)C spin-spin coupling constants have been performed in the series of selenium heterocycles and their parent open-chain selenides. It has been found that relativistic effects play an essential role in the selenium-carbon coupling mechanism and could result in a contribution of as much as 15-25% of the total values of the one-bond selenium-carbon spin-spin coupling constants. In the overall contribution of the relativistic effects to the total values of (1)J(Se,C), the scalar relativistic corrections (negative in sign) by far dominate over the spin-orbit ones (positive in sign), the latter being of less than 5%, as compared to the former (ca 20%). A combination of nonrelativistic second-order polarization propagator approach (CC2) with the four-component relativistic density functional theory scheme is recommended as a versatile tool for the calculation of (1)J(Se,C). Solvent effects in the values of (1)J(Se,C) calculated within the polarizable continuum model for the solvents with different dielectric constants (ε 2.2-78.4) are next to negligible decreasing negative (1)J(Se,C) in absolute value by only about 1 Hz. The use of the locally dense basis set approach applied herewith for the calculation of (77)Se-(13)C spin-spin coupling constants is fully justified resulting in a dramatic decrease in computational cost with only 0.1-0.2-Hz loss of accuracy.
Collapse
Affiliation(s)
- Yury Yu Rusakov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033, Irkutsk, Russia
| | | | | |
Collapse
|
40
|
Samultsev DO, Semenov VA, Krivdin LB. On the accuracy of the GIAO-DFT calculation of 15N NMR chemical shifts of the nitrogen-containing heterocycles--a gateway to better agreement with experiment at lower computational cost. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2014; 52:222-230. [PMID: 24573615 DOI: 10.1002/mrc.4055] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2013] [Revised: 01/21/2014] [Accepted: 01/31/2014] [Indexed: 06/03/2023]
Abstract
The main factors affecting the accuracy and computational cost of the gauge-independent atomic orbital density functional theory (GIAO-DFT) calculation of (15)N NMR chemical shifts in the representative series of key nitrogen-containing heterocycles--azoles and azines--have been systematically analyzed. In the calculation of (15)N NMR chemical shifts, the best result has been achieved with the KT3 functional used in combination with Jensen's pcS-3 basis set (GIAO-DFT-KT3/pcS-3) resulting in the value of mean absolute error as small as 5 ppm for a range exceeding 270 ppm in a benchmark series of 23 compounds with an overall number of 41 different (15)N NMR chemical shifts. Another essential finding is that basically, the application of the locally dense basis set approach is justified in the calculation of (15)N NMR chemical shifts within the 3-4 ppm error that results in a dramatic decrease in computational cost. Based on the present data, we recommend GIAO-DFT-KT3/pcS-3//pc-2 as one of the most effective locally dense basis set schemes for the calculation of (15)N NMR chemical shifts.
Collapse
Affiliation(s)
- Dmitry O Samultsev
- Siberian Branch of the Russian Academy of Sciences, A. E. Favorsky Irkutsk Institute of Chemistry, Favorsky St. 1, 664033, Irkutsk, Russia
| | | | | |
Collapse
|
41
|
Rusakov YY, Krivdin LB, Østerstrøm FF, Sauer SPA, Potapov VA, Amosova SV. First example of a high-level correlated calculation of the indirect spin-spin coupling constants involving tellurium: tellurophene and divinyl telluride. Phys Chem Chem Phys 2014; 15:13101-13107. [PMID: 23824065 DOI: 10.1039/c3cp51462e] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This paper documents the very first example of a high-level correlated calculation of spin-spin coupling constants involving tellurium taking into account relativistic effects, vibrational corrections and solvent effects for medium sized organotellurium molecules. The (125)Te-(1)H spin-spin coupling constants of tellurophene and divinyl telluride were calculated at the SOPPA and DFT levels, in good agreement with experimental data. A new full-electron basis set, av3z-J, for tellurium derived from the "relativistic" Dyall's basis set, dyall.av3z, and specifically optimized for the correlated calculations of spin-spin coupling constants involving tellurium was developed. The SOPPA method shows a much better performance compared to DFT, if relativistic effects calculated within the ZORA scheme are taken into account. Vibrational and solvent corrections are next to negligible, while conformational averaging is of prime importance in the calculation of (125)Te-(1)H spin-spin couplings. Based on the performed calculations at the SOPPA(CCSD) level, a marked stereospecificity of geminal and vicinal (125)Te-(1)H spin-spin coupling constants originating in the orientational lone pair effect of tellurium has been established, which opens a new guideline in organotellurium stereochemistry.
Collapse
Affiliation(s)
- Yury Yu Rusakov
- A.E. Favorsky Irkutsk Institute of Chemistry, Russian Academy of Sciences, Siberian Branch, 1 Favorsky St., 664033 Irkutsk, Russia.
| | - Leonid B Krivdin
- A.E. Favorsky Irkutsk Institute of Chemistry, Russian Academy of Sciences, Siberian Branch, 1 Favorsky St., 664033 Irkutsk, Russia.
| | - Freja F Østerstrøm
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| | - Vladimir A Potapov
- A.E. Favorsky Irkutsk Institute of Chemistry, Russian Academy of Sciences, Siberian Branch, 1 Favorsky St., 664033 Irkutsk, Russia.
| | - Svetlana V Amosova
- A.E. Favorsky Irkutsk Institute of Chemistry, Russian Academy of Sciences, Siberian Branch, 1 Favorsky St., 664033 Irkutsk, Russia.
| |
Collapse
|
42
|
Aucar IA, Gómez SS, Melo JI, Giribet CC, Ruiz de Azúa MC. Theoretical study of the nuclear spin-molecular rotation coupling for relativistic electrons and non-relativistic nuclei. II. Quantitative results in HX (X=H,F,Cl,Br,I) compounds. J Chem Phys 2013; 138:134107. [DOI: 10.1063/1.4796461] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
43
|
Rusakov YY, Krivdin LB. Modern quantum chemical methods for calculating spin–spin coupling constants: theoretical basis and structural applications in chemistry. RUSSIAN CHEMICAL REVIEWS 2013. [DOI: 10.1070/rc2013v082n02abeh004350] [Citation(s) in RCA: 99] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
44
|
Frank A, Möller HM, Exner TE. Toward the Quantum Chemical Calculation of NMR Chemical Shifts of Proteins. 2. Level of Theory, Basis Set, and Solvents Model Dependence. J Chem Theory Comput 2012; 8:1480-92. [PMID: 26596758 DOI: 10.1021/ct200913r] [Citation(s) in RCA: 63] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
It has been demonstrated that the fragmentation scheme of our adjustable density matrix assembler (ADMA) approach for the quantum chemical calculations of very large systems is well-suited to calculate NMR chemical shifts of proteins [ Frank et al. Proteins2011, 79, 2189-2202 ]. The systematic investigation performed here on the influences of the level of theory, basis set size, inclusion or exclusion of an implicit solvent model, and the use of partial charges to describe additional parts of the macromolecule on the accuracy of NMR chemical shifts demonstrates that using a valence triple-ζ basis set leads to large improvement compared to the results given in the previous publication. Additionally, moving from the B3LYP to the mPW1PW91 density functional and including partial charges and implicit solvents gave the best results with mean absolute errors of 0.44 ppm for hydrogen atoms excluding H(N) atoms and between 1.53 and 3.44 ppm for carbon atoms depending on the size and also on the accuracy of the protein structure. Polar hydrogen and nitrogen atoms are more difficult to predict. For the first, explicit hydrogen bonds to the solvents need to be included and, for the latter, going beyond DFT to post-Hartree-Fock methods like MP2 is probably required. Even if empirical methods like SHIFTX+ show similar performance, our calculations give for the first time very reliable chemical shifts that can also be used for complexes of proteins with small-molecule ligands or DNA/RNA. Therefore, taking advantage of its ab initio nature, our approach opens new fields of application that would otherwise be largely inaccessible due to insufficient availability of data for empirical parametrization.
Collapse
Affiliation(s)
- Andrea Frank
- Department of Chemistry and Zukunftskolleg, University of Konstanz , D-78457 Konstanz, Germany
| | - Heiko M Möller
- Department of Chemistry and Zukunftskolleg, University of Konstanz , D-78457 Konstanz, Germany
| | - Thomas E Exner
- Department of Chemistry and Zukunftskolleg, University of Konstanz , D-78457 Konstanz, Germany.,Department of Pharmaceutical and Medicinal Chemistry, Institute of Pharmacy , Eberhard Karls University Tübingen, D-72076 Tübingen, Germany
| |
Collapse
|
45
|
Kupka T, Nieradka M, Stachów M, Pluta T, Nowak P, Kjær H, Kongsted J, Kaminsky J. Basis Set Convergence of Indirect Spin–Spin Coupling Constants in the Kohn–Sham Limit for Several Small Molecules. J Phys Chem A 2012; 116:3728-38. [DOI: 10.1021/jp212588h] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Teobald Kupka
- Faculty of Chemistry, University of Opole, 45-052 Opole, Poland
| | | | - Michał Stachów
- Faculty of Chemistry, University of Opole, 45-052 Opole, Poland
| | - Tadeusz Pluta
- Institute of Chemistry, University of Silesia, 40-007 Katowice, Poland
| | - Piotr Nowak
- Technical University of Opole, Opole, Poland
| | - Hanna Kjær
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry
and Pharmacy, University of Southern Denmark, DK-5230 Odense, Denmark
| | - Jakub Kaminsky
- Institute of Organic Chemistry and Biochemistry, 166 10 Prague, Czech Republic
| |
Collapse
|
46
|
NMR espectroscopic parameters of HX and Si(Sn)X4 (X=H, F, Cl, Br and I) and SnBr4−nIn model compounds. Chem Phys 2012. [DOI: 10.1016/j.chemphys.2011.08.020] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
47
|
Melo JI, Maldonado A, Aucar GA. Relativistic effects on the shielding of SnH2XY and PbH2XY (X, Y = F, Cl, Br and I) heavy atom–containing molecules. Theor Chem Acc 2011. [DOI: 10.1007/s00214-010-0886-4] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
48
|
Provasi PF, Sauer SPA. Optimized basis sets for the calculation of indirect nuclear spin-spin coupling constants involving the atoms B, Al, Si, P, and Cl. J Chem Phys 2010; 133:054308. [PMID: 20707533 DOI: 10.1063/1.3465553] [Citation(s) in RCA: 73] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The aug-cc-pVTZ-J series of basis sets for indirect nuclear spin-spin coupling constants has been extended to the atoms B, Al, Si, P, and Cl. The basis sets were obtained according to the scheme previously described by Provasi et al. [J. Chem. Phys. 115, 1324 (2001)]. First, the completely uncontracted correlation consistent aug-cc-pVTZ basis sets were extended with four tight s and three tight d functions. Second, the s and p basis functions were contracted with the molecular orbital coefficients of self-consistent-field calculations performed with the uncontracted basis sets on the simplest hydrides of each atom. As a first illustration, we have calculated the one-bond indirect spin-spin coupling constants in BH(4)(-), BF, AlH, AlF, SiH(4), SiF(4), PH(3), PF(3), H(2)S, SF(6), HCl, and ClF at the level of density functional theory using the Becke three parameter Lee-Yang-Parr and the second order polarization propagator approximation with coupled cluster singles and doubles amplitudes.
Collapse
Affiliation(s)
- Patricio F Provasi
- Department of Physics, University of Northeastern, Av. Libertad 5500, W 3404 AAS Corrientes, Argentina.
| | | |
Collapse
|
49
|
Aucar GA, Romero RH, Maldonado AF. Polarization propagators: A powerful theoretical tool for a deeper understanding of NMR spectroscopic parameters. INT REV PHYS CHEM 2010. [DOI: 10.1080/01442350903432865] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
50
|
Rusakov YY, Krivdin LB, Sauer SPA, Levanova EP, Levkovskaya GG. Structural trends of 77Se-1H spin-spin coupling constants and conformational behavior of 2-substituted selenophenes. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2010; 48:44-52. [PMID: 19877126 DOI: 10.1002/mrc.2537] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Experimental measurements and second-order polarization propagator approach (SOPPA) calculations of (77)Se-(1)H spin-spin coupling constants together with theoretical energy-based conformational analysis in the series of 2-substituted selenophenes have been carried out. A new basis set optimized for the calculation of (77)Se-(1)H spin-spin coupling constants has been introduced by extending the aug-cc-pVTZ-J basis for selenium. Most of the spin-spin coupling constants under study, especially vicinal (77)Se-(1)H couplings, demonstrated a remarkable stereochemical behavior with respect to the internal rotation of the substituent in the 2-position of the selenophene ring, which is of major importance in the stereochemical studies of the related organoselenium compounds.
Collapse
Affiliation(s)
- Yury Yu Rusakov
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russia
| | | | | | | | | |
Collapse
|