1
|
Jessen LM, Sauer SPA. On the performance of HRPA(D) for NMR spin-spin coupling constants: Smaller molecules, aromatic and fluoroaromatic compounds. J Chem Phys 2024; 160:064102. [PMID: 38341775 DOI: 10.1063/5.0189932] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2023] [Accepted: 01/15/2024] [Indexed: 02/13/2024] Open
Abstract
In this study, the performance of the doubles-corrected higher random-phase approximation [HRPA(D)] has been investigated in calculations of nuclear magnetic resonance spin-spin coupling constants (SSCCs) for 58 molecules with the experimental values used as the reference values. HRPA(D) is an approximation to the second-order polarization propagator approximation (SOPPA) and is, therefore, computationally less expensive than SOPPA. HRPA(D) performs comparable and sometimes even better than SOPPA, and therefore, when calculating SSCCs, it should be considered as an alternative to SOPPA. Furthermore, it was investigated whether a coupled-cluster singles, doubles and perturbative triples [CCSD(T)] or Møller-Plesset second order (MP2) geometry optimization was optimal for a SOPPA and a HRPA(D) SSCC calculation for eight smaller molecules. CCSD(T) is the optimal geometry optimization for the SOPPA calculation, and MP2 was optimal for HRPA(D) SSCC calculations.
Collapse
Affiliation(s)
- Louise Møller Jessen
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, DK-2100 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
|
Karpov VV, Antonov AS, Tupikina EY. Choice of computational protocol for carbon-lithium spin-spin coupling constants 1 J CLi. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2022; 60:985-995. [PMID: 35881390 DOI: 10.1002/mrc.5299] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/29/2022] [Revised: 07/19/2022] [Accepted: 07/21/2022] [Indexed: 06/15/2023]
Abstract
In this work, we tested various computational schemes for calculations of 1 JCLi constants with a high accuracy. On the example of six organolithium reagents (phenyllithium monomer and dimer, monomer s-butyllithium, monomers of 1- and 2-lithionaphthalenes, and a methyllithium tetramer), the following aspects are discussed: (i) the role of a model system geometry, (ii) influence of solvent effects, and (iii) the choice of a functional and basis set. Practical recommendations for calculation of 1 JCLi with an accuracy ±2 Hz are formulated.
Collapse
Affiliation(s)
- Valerii V Karpov
- Institute of Chemistry, St. Petersburg State University, St. Petersburg, Russia
| | - Alexander S Antonov
- Institute of Chemistry, St. Petersburg State University, St. Petersburg, Russia
| | - Elena Yu Tupikina
- Institute of Chemistry, St. Petersburg State University, St. Petersburg, Russia
| |
Collapse
|
4
|
Jakubowska K, Pecul M, Ruud K. Vibrational Corrections to NMR Spin-Spin Coupling Constants from Relativistic Four-Component DFT Calculations. J Phys Chem A 2022; 126:7013-7020. [PMID: 36135807 PMCID: PMC9549459 DOI: 10.1021/acs.jpca.2c05019] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Zero-point vibrational
(ZPV) corrections to the nuclear spin–spin
coupling constants have been calculated using four-component Dirac–Kohn–Sham
DFT for H2X (where X = O, S, Se, Te, Po), XH3 (where X = N, P, As, Sb, Bi), and XH4 (where X = C, Si,
Ge, Sn, and Pb) molecules and for HC≡CPbH3. The
main goal was to study the influence of relativistic effects on the
ZPV corrections and thus results calculated at relativistic and nonrelativistic
approaches have been compared. The effects of relativity become notable
for the ZPV corrections to the spin–spin coupling constants
for compounds with lighter elements (selenium and germanium) than
for the spin–spin coupling constants themselves. In the case
of molecules containing heavier atoms, for instance BiH3 and PbH4, relativistic effects play a crucial role on
the results and approximating ZPV corrections by the nonrelativistic
results may lead to larger errors than omitting ZPV corrections altogether.
Collapse
Affiliation(s)
| | - Magdalena Pecul
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Poland
| | - Kenneth Ruud
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT─The Arctic University of Norway, N-9037 Tromsø, Norway.,Norwegian Defence Research Establishment, P.O. Box 25, 2027 Kjeller, Norway
| |
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
|
Reactions of a Bis(pentalene)dititanium complex with alkenes; the molecular structure of the butadiene complex [Ti2(µ: η5,η5-Pn††)2(μ: η2,η2-s-trans-C4H6)] (Pn†† = 1,4-(Si Pr3)2-C8H4). Polyhedron 2022. [DOI: 10.1016/j.poly.2021.115574] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
8
|
Giovanetti MD, Bitencourt LFF, Cormanich R, Sauer SPA. On the Unexpected Accuracy of the M06L Functional in the Calculation of 1JFC Spin-Spin Coupling Constants. J Chem Theory Comput 2021; 17:7712-7723. [PMID: 34751577 DOI: 10.1021/acs.jctc.1c00287] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
One-bond spin-spin coupling constants (SSCCs) between F and C are computed with density functional theory (DFT). Surprisingly, M06L stands out for its striking accuracy, outperforming any other investigated functional, including PBE0, otherwise considered one of the most reliable for couplings involving F. Although the computation of nuclear magnetic resonance (NMR) parameters involving F is known to be a challenging task, even with a rather small basis set as pcJ-1, M06L provides results with a MAD = 11.7 Hz, whereas the average deviation gets as much as 5 times larger for PBE0 (MAD = 60.0 Hz). In the context of SSCCs on the order of 300 Hz, this is particularly remarkable. We find that the accuracy of M06L/pcJ-1 in predicting 1JFC constants does not stem from a well-suited exchange or correlation part of the functional. Instead, it is believed to arise from a fortuitous cancellation of errors, as revealed by investigating the convergence of the basis set. Our findings also indicate that 1JFC constants are highly dependent on the amount of exact exchange included in the expression of the functional, with large fractions being critically important to achieving satisfactory results. Studying the effects of the geometry on 1JFC, we find that optimizing the geometry at the level of theory used to calculate SSCCs generally improves the quality of the results, although the combination of a M06-2X/aug-cc-pVTZ geometry with M06L/pcJ-1 1JFC constants best reproduces the experimental data for organofluorine compounds (with the exception of fluoroalkenes).
Collapse
Affiliation(s)
- Marinella de Giovanetti
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| | | | - Rodrigo Cormanich
- Chemistry Institute, State University of Campinas, P.O. Box 6154, 13083-970 Campinas, SP, Brazil
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| |
Collapse
|
9
|
Møller CHS, Schnack-Petersen AK, Sauer SPA. RPA(D) and HRPA(D): calculation of carbon–carbon spin–spin coupling constants for saturated cycloalkanes. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1757773] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
10
|
Adamson J, Nazarski RB, Jarvet J, Pehk T, Aav R. Shortfall of B3LYP in Reproducing NMR J
CH
Couplings in Some Isomeric Epoxy Structures with Strong Stereoelectronic Effects: A Benchmark Study on DFT Functionals. Chemphyschem 2018; 19:631-642. [DOI: 10.1002/cphc.201701125] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2017] [Revised: 11/17/2017] [Indexed: 01/13/2023]
Affiliation(s)
- Jasper Adamson
- National Institute of Chemical Physics and Biophysics; Akadeemia tee 23 12618 Tallinn Estonia
| | - Ryszard B. Nazarski
- Theoretical and Structural Chemistry Group; Department of Physical Chemistry; Faculty of Chemistry; University of Lodz; Pomorska 163/165 90-236 Łódź Poland
| | - Jüri Jarvet
- National Institute of Chemical Physics and Biophysics; Akadeemia tee 23 12618 Tallinn Estonia
- Current address: Department of Biochemistry and Biophysics; Arrhenius Laboratories; Stockholm University; Svante Arrhenius väg 16 10691 Stockholm Sweden
| | - Tõnis Pehk
- National Institute of Chemical Physics and Biophysics; Akadeemia tee 23 12618 Tallinn Estonia
| | - Riina Aav
- Department of Chemistry and Biotechnology; Tallinn University of Technology; Akadeemia tee 15 12618 Tallinn Estonia
| |
Collapse
|
11
|
Faber R, Kaminsky J, Sauer SPA. Rovibrational and Temperature Effects in Theoretical Studies of NMR Parameters. GAS PHASE NMR 2016. [DOI: 10.1039/9781782623816-00218] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The demand for high precision calculations of NMR shieldings (or their related values, chemical shifts δ) and spin–spin coupling constants facilitating and supporting detailed interpretations of NMR spectra increases hand in hand with the development of computational techniques and hardware resources. Highly sophisticated calculations including even relativistic effects are nowadays possible for these properties. However, NMR parameters depend not only on molecular structure and environment but also on molecular flexibility and temperature and the apparent success of theoretical predictions for molecular equilibrium geometries creates a demand for zero-point vibrational and temperature corrections. In this chapter we describe briefly the theory behind rovibrational corrections and review then some important contributions to this field.
Collapse
Affiliation(s)
- Rasmus Faber
- Department of Chemistry, University of Copenhagen Universitetsparken 5 DK-2100 Copenhagen Ø Denmark
| | - Jakub Kaminsky
- Department of Molecular Spectroscopy, Institute of Organic Chemistry and Biochemistry 166 10 Prague Czech Republic
| | - Stephan P. A. Sauer
- Department of Chemistry, University of Copenhagen Universitetsparken 5 DK-2100 Copenhagen Ø Denmark
| |
Collapse
|
12
|
Suardíaz R, Pérez C, Crespo-Otero R, García de la Vega JM, Fabián JS. Influence of Density Functionals and Basis Sets on One-Bond Carbon-Carbon NMR Spin-Spin Coupling Constants. J Chem Theory Comput 2015; 4:448-56. [PMID: 26620785 DOI: 10.1021/ct7003287] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The basis set and the functional dependence of one-bond carbon-carbon NMR spin-spin coupling constants (SSCC) have been analyzed using density functional theory. Four basis sets (6-311G**, TZVP, EPR-III, and aug-cc-pVTZ-J) and four functionals (PBE, PW91, B3LYP, and B3P86) are tested by comparison with 70 experimental values corresponding to 49 molecules that represent multiple types of hybridization of the carbon atoms. The two hybrid functionals B3P86 and B3LYP combined either EPR-III or aug-cc-pVTZ-J basis sets lead to the best accuracy of calculated SSCC. However, a simple linear regression allows for the obtaining of scaled coupling constants that fit much better with the experimental data and where the differences between the different basis sets and/or functional results are significantly reduced. For large molecules the TZVP basis set can be an appropriate election presenting a good compromise between quality of results and computational cost.
Collapse
Affiliation(s)
- R Suardíaz
- Departmento de Química Física, Facultad de Química, Universidad de la Habana, La Habana 10400, Cuba, and Departmento de Química Física Aplicada, Facultad de Ciencias, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - C Pérez
- Departmento de Química Física, Facultad de Química, Universidad de la Habana, La Habana 10400, Cuba, and Departmento de Química Física Aplicada, Facultad de Ciencias, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - R Crespo-Otero
- Departmento de Química Física, Facultad de Química, Universidad de la Habana, La Habana 10400, Cuba, and Departmento de Química Física Aplicada, Facultad de Ciencias, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - José M García de la Vega
- Departmento de Química Física, Facultad de Química, Universidad de la Habana, La Habana 10400, Cuba, and Departmento de Química Física Aplicada, Facultad de Ciencias, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Jesús San Fabián
- Departmento de Química Física, Facultad de Química, Universidad de la Habana, La Habana 10400, Cuba, and Departmento de Química Física Aplicada, Facultad de Ciencias, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| |
Collapse
|
13
|
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
|
14
|
Xu H, Saebo S, Pittman CU. The potential energy surface of singlet cyclobutadiene and substituted analogs: a coupled-cluster study. Struct Chem 2013. [DOI: 10.1007/s11224-013-0352-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
15
|
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]
|
16
|
Zarycz N, Aucar GA. The analysis of NMR J-couplings of saturated and unsaturated compounds by the localized second order polarization propagator approach method. J Chem Phys 2012; 136:174115. [PMID: 22583218 DOI: 10.1063/1.3697844] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Calculations of NMR J-coupling with polarization propagators are not invariant under unitary transformations at second order level of approach, second order polarization propagator approach (SOPPA). They are only invariant at first order or random phase level of approach (RPA). We performed "localized" SOPPA (Loc-SOPPA), calculations of J-couplings applying two different schemes for the localization of molecular orbitals(LMO): Foster-Boys and Pipek-Mezey. We show here that results of such Loc-SOPPA calculations are different though not much: they are less than 6% different in the worst case. Therefore it is possible to apply them with confidence in the analysis of the transmission of different coupling mechanisms within the molecule. We are able now to get reliable information on what LMOs are the most important (and so which are not important) for a given J-coupling in a molecule. This information can then be used for selecting which are the paths that should be described with the highest possible accuracy for that J-coupling calculation. A few unsaturated compounds are analyzed: ethene, trans-difluoroethene or DiF-ethene, and imine. It is shown that different lone pairs (of p(z) or p(x/y) type) are responsible for the vicinal F-F J-coupling in DiF-ethene; and also the fact that the main LP contributor is not the same for the fermi contact and the spin-dipolar mechanisms. We also studied phosphorous containing compounds such as phosphine and cis-propylene phosphine. In both cases the analysis of the main LMO contributing to one-bond P-H coupling and through-space P-C coupling were performed. The above mentioned unsaturated molecular systems have quasiinstability problems that arise at RPA level of approach. We show here that they are mostly originated in the antibonding π∗ LMO, corresponding to the C=C or C=N double bonds. We performed the analysis of the origin of quasiinstabilities for the SD mechanism. The contribution of each kind of excitation terms to SOPPA calculations were considered, meaning the main contributions by single and double excitations. It is shown that one can get more than 97% of the total electron correlation contribution when including terms that mainly contain single excitations (though double-excitation matrix elements should still be calculated).
Collapse
Affiliation(s)
- Natalia Zarycz
- Department of Physics, Northeastern University of Argentina and IMIT Institute, Avenida Libertad 5460 - Corrientes, W 3404 AAS Argentina
| | | |
Collapse
|
17
|
Sneskov K, Stanton JF. Effects of vibrational averaging on coupled cluster calculations of spin–spin coupling constants for hydrocarbons. Mol Phys 2012. [DOI: 10.1080/00268976.2012.678904] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Kristian Sneskov
- a The Lundbeck Foundation Center for Theoretical Chemistry and Department of Chemistry , Aarhus University , Langelandsgade 140, DK-8000 Aarhus C , Denmark
| | - John F. Stanton
- b Institute for Theoretical Chemistry and Department of Chemistry and Biochemistry , The University of Texas at Austin , Austin , Texas 78712-0165 , USA
| |
Collapse
|
18
|
Helgaker T, Coriani S, Jørgensen P, Kristensen K, Olsen J, Ruud K. Recent Advances in Wave Function-Based Methods of Molecular-Property Calculations. Chem Rev 2012; 112:543-631. [DOI: 10.1021/cr2002239] [Citation(s) in RCA: 463] [Impact Index Per Article: 38.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Trygve Helgaker
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Sonia Coriani
- Dipartimento di Scienze Chimiche e Farmaceutiche, Università degli Studi di Trieste, Via Giorgieri 1, I-34127 Trieste, Italy
| | - Poul Jørgensen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Kasper Kristensen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Jeppe Olsen
- Lundbeck Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | - Kenneth Ruud
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø, N-9037 Tromsø, Norway
| |
Collapse
|
19
|
Bally T, Rablen PR. Quantum-Chemical Simulation of 1H NMR Spectra. 2. Comparison of DFT-Based Procedures for Computing Proton–Proton Coupling Constants in Organic Molecules. J Org Chem 2011; 76:4818-30. [DOI: 10.1021/jo200513q] [Citation(s) in RCA: 144] [Impact Index Per Article: 11.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Thomas Bally
- Department of Chemistry, University of Fribourg, CH-1700 Fribourg, Switzerland
| | - Paul R. Rablen
- Department of Chemistry & Biochemistry, Swarthmore College, Swarthmore, Pennsylvania 19081-1397, United States
| |
Collapse
|
20
|
Benchmarking SOPPA(CC2) for the calculation of indirect nuclear spin–spin coupling constants: Carbocycles. Chem Phys 2011. [DOI: 10.1016/j.chemphys.2011.01.006] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
21
|
Kjaer H, Sauer SPA, Kongsted J. Benchmarking NMR indirect nuclear spin-spin coupling constants: SOPPA, SOPPA(CC2), and SOPPA(CCSD) versus CCSD. J Chem Phys 2011; 133:144106. [PMID: 20949986 DOI: 10.1063/1.3483197] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurate calculations of NMR indirect nuclear spin-spin coupling constants require especially optimized basis sets and correlated wave function methods such as CCSD or SOPPA(CCSD). Both methods scale as N(6), where N is the number of orbitals, which prevents routine applications to molecules with more than 10-15 nonhydrogen atoms. We have therefore developed a modification of the SOPPA(CCSD) method in which the CCSD singles and doubles amplitudes are replaced by CC2 singles and doubles amplitudes. This new method, called SOPPA(CC2), scales only as N(5), like the original SOPPA-method. The performance of the SOPPA(CC2) method for the calculation of indirect nuclear spin-spin coupling constants is compared to SOPPA and SOPPA(CCSD) employing a set of benchmark molecules. We also investigate the basis set dependence by employing three different basis sets optimized for spin-spin coupling constants, namely the HuzIV-su4, ccJ-pVTZ, and ccJ-pVQZ basis sets. The results of the corresponding CCSD calculations are used as a theoretical reference.
Collapse
Affiliation(s)
- Hanna Kjaer
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark.
| | | | | |
Collapse
|
22
|
Colherinhas G, Fonseca T, Castro M. 13C chemical shifts of polyacetylene chains with charged conformational defects: A GIAO–DFT study. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.01.014] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
23
|
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
|
24
|
Yachmenev A, Yurchenko SN, Paidarová I, Jensen P, Thiel W, Sauer SPA. Thermal averaging of the indirect nuclear spin-spin coupling constants of ammonia: The importance of the large amplitude inversion mode. J Chem Phys 2010; 132:114305. [DOI: 10.1063/1.3359850] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
25
|
Esteban AL, Díez E, Galache MP, San Fabián J, Casanueva J, Contreras RH. Vibrational contributions to vicinal proton–proton coupling constants3JHH. Mol Phys 2010. [DOI: 10.1080/00268971003630687] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
26
|
Alkorta I, Blanco F, Elguero J. A SOPPA theoretical study of the spin–spin coupling constants of all fluorobenzenes C6HnF6−n (n=0–5). J Mol Struct 2010. [DOI: 10.1016/j.molstruc.2009.11.019] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
27
|
|
28
|
Jackowski K, Makulski W, Szyprowska A, Antusek A, Jaszuński M. Temperature dependence of the (1)J((11)B(19)F) spin-spin coupling in BF(3) molecule. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2009; 47:857-861. [PMID: 19637211 DOI: 10.1002/mrc.2484] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
The (1)J((11)B(19)F) spin-spin coupling of gaseous BF(3) was observed in (11)B NMR spectra as a function of density in a wide range of temperatures. Following the extrapolation of the measured values to the zero-density limit, the coupling constant free from intermolecular effects (1)J(0)((11)B(19)F) was obtained for each temperature. In contrast to previous investigations, the final results indicate a nonlinear dependence of (1)J(0)((11)B(19)F) on temperature. In the corresponding ab initio calculations of spin-spin coupling constants performed at the coupled cluster singles and doubles (CCSD) level to obtain a reliable result for this coupling constant we had to take into account large vibrational corrections.
Collapse
Affiliation(s)
- Karol Jackowski
- Laboratory of NMR Spectroscopy, Department of Chemistry, Warsaw University, 02-093 Warszawa, Poland.
| | | | | | | | | |
Collapse
|
29
|
Møgelhøj A, Aidas K, Mikkelsen KV, Sauer SPA, Kongsted J. Prediction of spin-spin coupling constants in solution based on combined density functional theory/molecular mechanics. J Chem Phys 2009; 130:134508. [PMID: 19355752 DOI: 10.1063/1.3098255] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present theory and implementation of calculation of spin-spin coupling constants within combined quantum mechanics/molecular mechanics methods. Special attention is given to the role of explicit solvent polarization as well as the molecular consequences due to hydrogen bonding. The model is generally applicable but is here implemented for the case of density functional theory. First applications to liquid water and acetylene in aqueous solution are presented. Good agreement between theory and experiment is obtained in both cases, thereby showing the strength of our approach. Finally, spin-spin coupling constants across hydrogen bonds are discussed considering for the first time the role of an explicit solvent on this class of spin-spin couplings.
Collapse
Affiliation(s)
- Andreas Møgelhøj
- Department of Chemistry, H. C. Ørsted Institute, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| | | | | | | | | |
Collapse
|
30
|
Del Bene JE, Alkorta I, Elguero J. A Systematic Comparison of Second-Order Polarization Propagator Approximation and Equation-of-Motion Coupled Cluster Singles and Doubles C−C, C−N, N−N, C−H, and N−H Spin−Spin Coupling Constants. J Phys Chem A 2009; 113:12411-20. [DOI: 10.1021/jp902089g] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
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
|
31
|
Dracínský M, Kaminský J, Bour P. Relative importance of first and second derivatives of nuclear magnetic resonance chemical shifts and spin-spin coupling constants for vibrational averaging. J Chem Phys 2009; 130:094106. [PMID: 19275395 DOI: 10.1063/1.3081317] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Relative importance of anharmonic corrections to molecular vibrational energies, nuclear magnetic resonance (NMR) chemical shifts, and J-coupling constants was assessed for a model set of methane derivatives, differently charged alanine forms, and sugar models. Molecular quartic force fields and NMR parameter derivatives were obtained quantum mechanically by a numerical differentiation. In most cases the harmonic vibrational function combined with the property second derivatives provided the largest correction of the equilibrium values, while anharmonic corrections (third and fourth energy derivatives) were found less important. The most computationally expensive off-diagonal quartic energy derivatives involving four different coordinates provided a negligible contribution. The vibrational corrections of NMR shifts were small and yielded a convincing improvement only for very accurate wave function calculations. For the indirect spin-spin coupling constants the averaging significantly improved already the equilibrium values obtained at the density functional theory level. Both first and complete second shielding derivatives were found important for the shift corrections, while for the J-coupling constants the vibrational parts were dominated by the diagonal second derivatives. The vibrational corrections were also applied to some isotopic effects, where the corrected values reasonably well reproduced the experiment, but only if a full second-order expansion of the NMR parameters was included. Contributions of individual vibrational modes for the averaging are discussed. Similar behavior was found for the methane derivatives, and for the larger and polar molecules. The vibrational averaging thus facilitates interpretation of previous experimental results and suggests that it can make future molecular structural studies more reliable. Because of the lengthy numerical differentiation required to compute the NMR parameter derivatives their analytical implementation in future quantum chemistry packages is desirable.
Collapse
Affiliation(s)
- Martin Dracínský
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo namesti 2, 166 10 Prague, Czech Republic.
| | | | | |
Collapse
|
32
|
Del Bene JE, Alkorta I, Elguero J. 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 Molecules with C, N, and O Double and Triple Bonds and Selected F-Substituted Derivatives. J Chem Theory Comput 2008; 5:208-16. [DOI: 10.1021/ct800321b] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
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
|
Hansen MB, Kongsted J, Toffoli D, Christiansen O. Vibrational Contributions to Indirect Spin−Spin Coupling Constants Calculated via Variational Anharmonic Approaches. J Phys Chem A 2008; 112:8436-45. [DOI: 10.1021/jp804306s] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Mikkel B. Hansen
- The Lundbeck Foundation Center for Theoretical Chemistry and Center for Oxygen Microscopy and Imaging, Department of Chemistry, University of Århus, Langelandsgade 140, DK-8000 Århus C, Denmark, and Department of Theoretical Chemistry, Chemical Center, University of Lund, P.O. Box 124, S-221 00 Lund, Sweden
| | - Jacob Kongsted
- The Lundbeck Foundation Center for Theoretical Chemistry and Center for Oxygen Microscopy and Imaging, Department of Chemistry, University of Århus, Langelandsgade 140, DK-8000 Århus C, Denmark, and Department of Theoretical Chemistry, Chemical Center, University of Lund, P.O. Box 124, S-221 00 Lund, Sweden
| | - Daniele Toffoli
- The Lundbeck Foundation Center for Theoretical Chemistry and Center for Oxygen Microscopy and Imaging, Department of Chemistry, University of Århus, Langelandsgade 140, DK-8000 Århus C, Denmark, and Department of Theoretical Chemistry, Chemical Center, University of Lund, P.O. Box 124, S-221 00 Lund, Sweden
| | - Ove Christiansen
- The Lundbeck Foundation Center for Theoretical Chemistry and Center for Oxygen Microscopy and Imaging, Department of Chemistry, University of Århus, Langelandsgade 140, DK-8000 Århus C, Denmark, and Department of Theoretical Chemistry, Chemical Center, University of Lund, P.O. Box 124, S-221 00 Lund, Sweden
| |
Collapse
|
34
|
Benda L, Prokop A, Blechta V, Schraml J, Špirko V. Probing the Flexibility of Internal Rotation in Silylated Phenols with the NMR Scalar Spin−Spin Coupling Constants. J Phys Chem A 2008; 112:5167-74. [DOI: 10.1021/jp711809e] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ladislav Benda
- Institute of Organic Chemistry and Biochemistry of the ASCR, v.v.i., Flemingovo sq. 2., 166 10 Prague 6, Czech Republic, and Institute of Chemical Process Fundamentals of the ASCR, v.v.i., Rozvojová 135, 165 02, Prague 6, Czech Republic
| | - Alexandr Prokop
- Institute of Organic Chemistry and Biochemistry of the ASCR, v.v.i., Flemingovo sq. 2., 166 10 Prague 6, Czech Republic, and Institute of Chemical Process Fundamentals of the ASCR, v.v.i., Rozvojová 135, 165 02, Prague 6, Czech Republic
| | - Vratislav Blechta
- Institute of Organic Chemistry and Biochemistry of the ASCR, v.v.i., Flemingovo sq. 2., 166 10 Prague 6, Czech Republic, and Institute of Chemical Process Fundamentals of the ASCR, v.v.i., Rozvojová 135, 165 02, Prague 6, Czech Republic
| | - Jan Schraml
- Institute of Organic Chemistry and Biochemistry of the ASCR, v.v.i., Flemingovo sq. 2., 166 10 Prague 6, Czech Republic, and Institute of Chemical Process Fundamentals of the ASCR, v.v.i., Rozvojová 135, 165 02, Prague 6, Czech Republic
| | - Vladimír Špirko
- Institute of Organic Chemistry and Biochemistry of the ASCR, v.v.i., Flemingovo sq. 2., 166 10 Prague 6, Czech Republic, and Institute of Chemical Process Fundamentals of the ASCR, v.v.i., Rozvojová 135, 165 02, Prague 6, Czech Republic
| |
Collapse
|
35
|
San Fabián J, Díez E, García de la Vega JM, Suardíaz R. Approximating correlation effects in multiconfigurational self-consistent field calculations of spin-spin coupling constants. J Chem Phys 2008; 128:084108. [PMID: 18315034 DOI: 10.1063/1.2834210] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- J San Fabián
- Departamento de Química Física Aplicada, Facultad de Ciencias, Universidad Autónoma de Madrid, E-28049 Madrid, Spain.
| | | | | | | |
Collapse
|
36
|
Mort BC, Autschbach J. A Pragmatic Recipe for the Treatment of Hindered Rotations in the Vibrational Averaging of Molecular Properties. Chemphyschem 2008; 9:159-70. [DOI: 10.1002/cphc.200700628] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
37
|
Mort BC, Autschbach J. Magnitude of zero-point vibrational corrections to optical rotation in rigid organic molecules: a time-dependent density functional study. J Phys Chem A 2007; 109:8617-23. [PMID: 16834261 DOI: 10.1021/jp051685y] [Citation(s) in RCA: 94] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The zero-point vibrational corrections (ZPVCs) to the optical rotation of 22 rigid organic molecules have been calculated using time-dependent density functional theory with the B3LYP hybrid functional. We outline an implementation for calculating ZPVCs that can be applied with a variety of quantum chemistry programs and methods. It is shown that the ZPVCs to optical rotation have a wide range of values and can be quite significant depending on the molecule. On average, it has been determined that vibrational corrections can account for about 20% of the optical rotation for the equilibrium value. It is also concluded that vibrational effects alone cannot be the only factor in improving the calculated values of optical rotation with respect to experimental data measured in solution.
Collapse
Affiliation(s)
- Brendan C Mort
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, USA
| | | |
Collapse
|
38
|
Contreras RH, Esteban AL, Díez E, Lochert IJ, Della EW, Tormena CF. Experimental and DFT studies on the transmission mechanisms of analogous NMR JCH and JCC couplings in 1-X- and 1-X-3-methylbicyclo[1.1.1]-pentanes. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2007; 45:572-7. [PMID: 17534878 DOI: 10.1002/mrc.2009] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
The main aim of this work is to compare the transmission mechanisms for the Fermi contact term of spin-spin couplings, SSCCs, in series 1-X-bicyclo[1.1.1]-pentane, (1), and 1-X-3-methylbicyclo[1.1.1]pentane, (2), and from that comparison to gain insight into some subtle aspects of the FC transmission. To this end, 18 members of the latter series were isotopically enriched in (13)C at the methyl position and the following couplings were measured; 1JC3CMe, 3JC1CMe and 4JCXCMe. These three types of SSCCs in (2) are compared, respectively, with 1JC3H3, 3JC1H3 and 4JCXH in (1); these latter values were taken from previous works. Since electron delocalization plays an important role in the transmission of the FC interaction, the natural bond orbital (NBO) method is employed to quantify electron delocalization interactions within selected members of series (1) and (2). It is found that 1JC3H3 SSCCs in (1) is more efficiently transmitted than 1JC3CMe SSCCs in (2). On the other hand, 3JC1H3 and 4JCXH SSCCs in (1) are notably less efficiently transmitted than 3JC1CMe and 4JCXCMe SSCCs in (2), although substituent effects on these two SSCCs show the opposite trends. These different efficiencies are rationalized in terms of different sigma-hyperconjugative interactions in both series of compounds.
Collapse
Affiliation(s)
- Rubén H Contreras
- Department of Physics, FCEyN, University of Buenos Aires, Ciudad Universitaria, Pabellón 1 (C1428EHA) Buenos Aires, Argentina
| | | | | | | | | | | |
Collapse
|
39
|
Mort BC, Autschbach J. Vibrational Corrections to Magneto-Optical Rotation: A Computational Study. J Phys Chem A 2007; 111:5563-71. [PMID: 17539612 DOI: 10.1021/jp070448n] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Vibrational corrections to the Verdet constants of nine molecules (H2, N2, CO, H2O, CH4, benzene, toluene, p-xylene, and o-xylene) were calculated with pure density functional theory (DFT), hybrid DFT, and an approximate coupled-cluster theory. Comparisons are made for the accuracy of the vibrational averages among different methods and with respect to experimental data where available. It is found that vibrational corrections to magneto-optical rotation can be as large as 10% of the equilibrium value. Hybrid DFT with the B3LYP hybrid functional offers reasonable accuracy at a relatively inexpensive computational cost for accurate calculations of vibrationally averaged Verdet constants.
Collapse
Affiliation(s)
- Brendan C Mort
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, NY 14260-3000, USA
| | | |
Collapse
|
40
|
Mort BC, Autschbach J. Temperature Dependence of the Optical Rotation in Six Bicyclic Organic Molecules Calculated by Vibrational Averaging. Chemphyschem 2007; 8:605-16. [PMID: 17304606 DOI: 10.1002/cphc.200600757] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The vibrational corrections and the temperature dependence of the specific rotation of six rigid organic molecules (alpha-pinene, beta-pinene, cis-pinane, camphene, camphor, and fenchone) were calculated at three wavelengths using hybrid time-dependent density functional theory (TDDFT). A technique for calculating the temperature dependence of the vibrational average of a molecular property has been applied to obtain the specific rotation of the molecules as a function of temperature. For cases in which accurate equilibrium optical rotations can be obtained as a "base value," and for which there is little effect from solvation, accurate predictions of the trends in the temperature-dependence of the specific rotations can be calculated. For other cases, the method can be used to extract purely vibrational contributions to the overall temperature dependence of optical rotation.
Collapse
Affiliation(s)
- Brendan C Mort
- Department of Chemistry, State University of New York at Buffalo, Buffalo, New York 14260-3000, USA
| | | |
Collapse
|
41
|
Wrackmeyer B, Hernández ZG, Herberhold M. 1-Cyclohepta-2,4,6-trienyl-selanes--a 77Se NMR study: indirect nuclear 77Se--13C spin-spin coupling constants and application of density functional theory (DFT) calculations. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2007; 45:198-204. [PMID: 17183533 DOI: 10.1002/mrc.1946] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
1-Cyclohepta-2,4,6-trienyl-selanes Se(C(7)H(7))(2) (2c), R--Se--C(7)H(7) with R = Bu, (t)Bu, Ph, 4-F--C(6)H(4) (12a,b,c,d) were prepared by the reaction of the corresponding silanes, Si(SeMe(3))(2) and R--Se--SiMe(3), respectively, with tropylium bromide C(7)H(7)Br. In spite of the low stability of the selanes even in dilute solutions and at low temperature, they could be characterised by their (1)H, (13)C and (77)Se NMR parameters. Coupling constants (1)J((77)Se,(13)C) were measured and calculated by DFT methods at the B3LYP/6-311+G(d,p) level of theory. The comparison of experimental and calculated coupling constants (1)J((77)Se,(13)C) included numerous selenium carbon compounds with largely different Se--C bonds, revealing a satisfactory agreement. Both the spin-dipole (SD) and the paramagnetic spin-orbital (PSO) terms contributed significantly to the spin-spin coupling interaction, in addition to the Fermi contact (FC) term.
Collapse
Affiliation(s)
- Bernd Wrackmeyer
- Anorganische Chemie II, Universität Bayreuth, D-95440 Bayreuth, Germany.
| | | | | |
Collapse
|
42
|
van Mourik T, Dingley AJ. Geometry Dependence of Spin–Spin Couplings in Cyanamide by DFT Analysis. Chemphyschem 2007; 8:288-96. [PMID: 17221902 DOI: 10.1002/cphc.200600489] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
There have been numerous theoretical and experimental investigations examining NMR parameters related to non-amino N-H...N H-bonded moieties in both biological and chemical contexts. In contrast, little information on the geometry dependence of NMR parameters related to the biologically important H-bond donor amino group is available. Herein, the geometric dependencies of the one-bond amino N-H spin-spin coupling constants [(1)J(NH)] in the cyanamide monomer and dimer have been computed with B3LYP and the aug-cc-pVTZ-su0 basis set. In an isolated planar cyanamide molecule, the |(1)J(NH)| couplings were found to increase as the N-H bond lengthened. In contrast, in the planar cyanamide dimer the size of the H-bonded amino N-H coupling (|(1)J(N(d)H(d))|) decreased with increasing N(d)H(d) bond length. The |(1)J(N(d)H(d))| coupling was larger than the |(1)J(N(d)H(free))| coupling for N(d)H(d) distances up to 1.18 A (for a fixed N(d)H(free) distance of 1.006 A). Hence, the decrease of |(1)J(NH)| with increasing N-H distance, as well as the larger value of |(1)J(N(d)H(d))| compared to |(1)J(N(d)H(free))|, were only observed for situations where the amino group is involved in an H-bonding interaction. This is attributed to electron redistribution induced by the presence of the second cyanamide molecule. Similar electron-redistribution effects are thought to be responsible for the observed distance dependence of computed (1)J(NH) couplings of H-bonded amino groups in near-planar G-quartet structures. Here, the |(1)J(NH)| couplings of the amino N-H bonds decreased with increasing N-H bond length whereas the |(1)J(N(d)H(d))| couplings are approximately 7 Hz larger than the |(1)J(N(d)H(free))| couplings, despite the longer N(d)-H(d) bond length.
Collapse
Affiliation(s)
- Tanja van Mourik
- Chemistry Department, University College London, 20 Gordon Street, London WC1H 0AJ, UK.
| | | |
Collapse
|
43
|
Helgaker T, Ola B. Lutnæs, Jaszuński M. Density-Functional and Coupled-Cluster Singles-and-Doubles Calculations of the Nuclear Shielding and Indirect Nuclear Spin−Spin Coupling Constants of o-Benzyne. J Chem Theory Comput 2006; 3:86-94. [DOI: 10.1021/ct600234n] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Trygve Helgaker
- Department of Chemistry, University of Durham, South Road, Durham, DH1 3LE, United Kingdom
| | - Ola B. Lutnæs
- Department of Chemistry, University of Oslo, P.O.B. 1033 Blindern, N-0315 Oslo, Norway
| | - Michał Jaszuński
- Institute of Organic Chemistry, Polish Academy of Sciences, 01-224 Warszawa, Kasprzaka 44, Poland
| |
Collapse
|
44
|
Kongsted J, Christiansen O. Automatic generation of force fields and property surfaces for use in variational vibrational calculations of anharmonic vibrational energies and zero-point vibrational averaged properties. J Chem Phys 2006; 125:124108. [PMID: 17014167 DOI: 10.1063/1.2352734] [Citation(s) in RCA: 65] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An automatic and general procedure for the calculation of geometrical derivatives of the energy and general property surfaces for molecular systems is developed and implemented. General expressions for an n-mode representation are derived, where the n-mode representation includes only the couplings between n or less degrees of freedom. The general expressions are specialized to derivative force fields and property surfaces, and a scheme for calculation of the numerical derivatives is implemented. The implementation is interfaced to electronic structure programs and may be used for both ground and excited electronic states. The implementation is done in the context of a vibrational structure program and can be used in combination with vibrational self-consistent field (VSCF), vibrational configuration interaction (VCI), vibrational Moller-Plesset, and vibrational coupled cluster calculations of anharmonic wave functions and calculation of vibrational averaged properties at the VSCF and VCI levels. Sample calculations are presented for fundamental vibrational energies and vibrationally averaged dipole moments and frequency dependent polarizabilities and hyperpolarizabilities of water and formaldehyde.
Collapse
Affiliation(s)
- Jacob Kongsted
- Department of Chemistry, University of Aarhus, DK-8000 Aarhus C, Denmark.
| | | |
Collapse
|
45
|
|
46
|
Watson MA, Sałek P, Macak P, Jaszuński M, Helgaker T. The calculation of indirect nuclear spin-spin coupling constants in large molecules. Chemistry 2006; 10:4627-39. [PMID: 15378642 DOI: 10.1002/chem.200306065] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
We present calculations of indirect nuclear spin-spin coupling constants in large molecular systems, performed using density functional theory. Such calculations, which have become possible because of the use of linear-scaling techniques in the evaluation of the Coulomb and exchange-correlation contributions to the electronic energy, allow us to study indirect spin-spin couplings in molecules of biological interest, without having to construct artificial model systems. In addition to presenting a statistical analysis of the large number of short-range coupling constants in large molecular systems, we analyse the asymptotic dependence of the indirect nuclear spin-spin coupling constants on the internuclear separation. In particular, we demonstrate that, in a sufficiently large one-electron basis set, the indirect spin-spin coupling constants become proportional to the inverse cube of the internuclear separation, even though the diamagnetic and paramagnetic spin-orbit contributions to the spin-spin coupling constants separately decay as the inverse square of this separation. By contrast, the triplet Fermi contact and spin-dipole contributions to the indirect spin-spin coupling constants decay exponentially and as the inverse cube of the internuclear separation, respectively. Thus, whereas short-range indirect spin-spin coupling constants are usually dominated by the Fermi contact contribution, long-range coupling constants are always dominated by the negative diamagnetic spin-orbit contribution and by the positive paramagnetic spin-orbit contribution, with small spin-dipole and negligible Fermi contact contributions.
Collapse
Affiliation(s)
- Mark A Watson
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK
| | | | | | | | | |
Collapse
|
47
|
Dodziuk H, Jaszuński M, Schilf W. 1H and 13C NMR chemical shifts and spin-spin coupling constants in trans- and cis-decalins. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2005; 43:639-46. [PMID: 15915544 DOI: 10.1002/mrc.1598] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
The NMR parameters characterizing the spectra of trans- and cis-decalins were determined from theoretical calculations and experimental spectra. The calculated values of the shielding constants are in good agreement with the measured chemical shifts, with a small but noticeable difference in accuracy for the bridgehead atoms. Of all the spin-spin coupling constants, only most of (1)J(C,C) and (1)J(C,H) values could be extracted from the spectra, and the corresponding computed values are in good agreement with experiment. It appears that the applied density functional theory (DFT) approach overestimates slightly the J(C,C) coupling and underestimates the differences between one-bond (1)J(C,H) coupling constants. For all these constants [J(C,C), J(C,H) and J(H,H)] through one to three bonds, which could not be obtained experimentally, the predicted values are in good agreement with the general rules relating spin-spin coupling to the number and spatial arrangement of the intervening bonds.
Collapse
Affiliation(s)
- Helena Dodziuk
- Institute of Physical Chemistry, Polish Academy of Sciences, Kasprzaka 44/52, 01-224 Warsaw, Poland.
| | | | | |
Collapse
|
48
|
Ramalho TC, Taft CA. Thermal and solvent effects on the NMR and UV parameters of some bioreductive drugs. J Chem Phys 2005; 123:054319. [PMID: 16108651 DOI: 10.1063/1.1996577] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
15N NMR chemical shifts and n-->pi* electronic transition energy for metronidazole (1) has been calculated and compared with experimental data. A detailed computational study of 1 is presented, with special attention to the performance of various theoretical methods for reproducing spectroscopic parameters in solution. The most sophisticated approach involves density functional based on the Car-Parrinello molecular dynamics simulations of 1 in aqueous solution (BP86 level) and averaging chemical shifts and deltaE(n-->pi*) over snapshots from the trajectory. In the NMR and UV calculations for these snapshots (performed at the B3LYP level), a small number of discrete water molecules are retained, and the remaining bulk solution effects are included via a polarizable continuum model (PCM). A good agreement with experiment is also obtained using static geometry optimization and NMR computation of pristine 1 employing a PCM approach. Further theoretical predictions are also reported for 17O NMR and deltaE(n-->pi*) of three hydroxycinnamic acid derivatives, which suggest that it is essential to incorporate the dynamics and solvent effects for NMR and UV calculations in the condensed phase.
Collapse
Affiliation(s)
- Teodorico C Ramalho
- Departamento de Físico-Química and Laboratório de Impacto de Fótons e Elétrons, Instituto de Química, Universidade Federal do Rio de Janeiro, Ilha do Fundão, Centro de Technologia, Bloco A, Lab. 406, Rio de Janeiro 21949-900, Rio de Janeiro, Brazil
| | | |
Collapse
|
49
|
Krivdin LB, Scherbina NA, Istomina NV. Non-empirical calculations of NMR indirect carbon-carbon coupling constants. Part 12--aliphatic and alicyclic oximes. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2005; 43:435-443. [PMID: 15809975 DOI: 10.1002/mrc.1572] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
One-bond carbon-carbon coupling constants were calculated in a series of nine aliphatic and alicyclic oximes at the SOPPA (second-order polarization propagator approach) level in good agreement with the available experimental data, and several unknown couplings were predicted with high reliability. The experimental difference between J(C,C) of the corresponding carbon-carbon bonds in cis and trans orientations to the nitrogen lone pair is very well reproduced at the SOPPA level, and this provides an additional tool in the configurational assignment at the C=N bond in oximes and related systems.
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
|
50
|
Díez * E, Casanueva J, San Fabián J, Esteban AL, Galache MP, Barone V, Peralta JE, Contreras RH. Prediction of vicinal proton–proton coupling constants3JHHfrom density functional theory calculations. Mol Phys 2005. [DOI: 10.1080/00268970412331333131] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|