Isotropic Nuclear Magnetic Shielding Research Articles

Effects of the Hubbard potential on the NMR shielding and optoelectronic properties of BiMnVO5 compound

This study explores the nuclear magnetic shielding, chemical shifts, and the optoelectronic properties of the BiMnVO5 compound using the full-potential linearized augmented plane wave method within the generalized gradient approximation by employing the Hubbard model (GGA + U). The 209Bi and 51V chemical shifts and bandgap values of the BiMnVO5 compound in a triclinic crystal structure are found to be directly related to Hubbard potential. The relationship between the isotropic nuclear magnetic shielding σiso and chemical shift δiso is obtained with a slope of 1.0231 and − 0.00188 for 209Bi and 51V atoms, respectively. It is also observed that the bandgap, isotropic nuclear magnetic shielding, and chemical shifts increase with the change in Hubbard potentials (U) of 3, 4, 5, 6, and 7.

Solvent Dependent Nuclear Magnetic Resonance Molecular Parameters Based on a Polarization Consistent Screened Range Separated Hybrid Density Functional Theory Framework.

Nuclear magnetic resonance (NMR) properties of solvated molecules are significantly affected by the solvent. We, therefore, employ a polarization consistent framework that efficiently addresses the solvent polarizing environment effects. Toward this goal a dielectric screened range separated hybrid (SRSH) functional is invoked with a polarizable continuum model (PCM) to properly represent the orbital gap in the condensed phase. We build on the success of range separated hybrid (RSH) functionals to address the erroneous tendency of traditional density functional theory (DFT) to collapse the orbital gap. Recently, the impact of RSH that properly opens up the orbital gap in gas-phase calculations on NMR properties has been assessed. Here, we report the use of SRSH-PCM that produces properly solute orbital gaps in calculating isotropic nuclear magnetic shielding and chemical shift parameters of molecular systems in the condensed phase. We show that in contrast to simpler DFT-PCM approaches, SRSH-PCM successfully follows expected dielectric constant trends.

Efficiently Computing NMR 1H and 13C Chemical Shifts of Saccharides in Aqueous Environment.

Determining the structure of saccharides in their native environment is crucial to understanding their function and more accurately targeting their utilization. Nuclear magnetic resonance observables such as the nuclear Overhauser effect or spin-spin coupling constants are routinely utilized to study saccharides in their native water environment. However, while highly sensitive to the local environment, chemical shifts are mostly overlooked, despite being commonly measured for compounds identification. Although chemical shifts carry considerable structural information, their direct association with structure is notoriously difficult. This is mostly due to the similarity in the chemical nature of most saccharides causing similar physicochemical environments close to sugar C and H atoms, resulting in comparable chemical shifts. The rise of computational power allows one to compute reliable chemical shifts and use them to determine atomistic details of these sugars in solution. However, any prediction is severely limited by the computational protocol used and its accuracy. In this work, we studied a set of 31 saccharides on which we evaluated various computational protocols to calculate the total number of 375 1H and 327 13C chemical shifts of sugars in an aqueous environment. Our study proposes two cost-effective protocols for simulating 1H and 13C chemical shifts that we recommend for further use. These protocols can help with the interpretation of experimental spectra, but we also show that they are also capable of structure prediction independently. This is possible because of the low mean absolute deviations of calculated shifts from the experiment (0.06 ppm for 1H and 1.09 ppm for 13C). We explore different solvation methods, basis sets, and optimization schemes to reach such accuracy. A correct sampling of the conformation phase space of flexible sugar molecules is also key to obtaining accurately converged theoretical chemical shifts. The linear regression method was applied to convert the calculated isotropic nuclear magnetic shielding constants to simulated chemical shifts comparable with the experiment. The achieved level of accuracy can help in utilizing chemical shifts for elucidating the 3D atomistic structure of saccharides in aqueous solutions. All linear regression parameters obtained on our extensive set of sugars for all the tested protocols can be reutilized in future works.

Assessment of the Performance of Optimally Tuned Range‐Separated Hybrid Functionals for Nuclear Magnetic Shielding Calculations

AbstractThe performance of optimally tuned range‐separated hybrid (OT‐RSH) functional calculations in predicting accurate isotropic nuclear magnetic shielding (σ) and chemical shift values is examined. To that end, the results of OT‐RSH and other approximate density functional theory calculations are assessed against recently published benchmark CCSD(T) calculations for a test set consisting of several molecules and bond types. It is found that for atoms in single bonds with a large paramagnetic contribution to σ, OT‐RSH offers a significant improvement in prediction of shielding constants over popular semi‐local and hybrid density functionals, yielding non‐empirical results that are as accurate as those of semi‐empirical density functionals often used for prediction of shielding constants. This success is attributed to the improved fundamental gap prediction of the OT‐RSH approach. For atoms in multiple bonds, however, larger errors often persist. By comparing OT‐RSH and recently reported double‐hybrid functional results, the remaining difficulties are traced to significant non‐local correlation.

Influence of Solvent Representation on Nuclear Shielding Calculations of Protonation States of Small Biological Molecules.

In this study, we assess the influence of solvation on the accuracy and reliability of isotropic nuclear magnetic shielding calculations for amino acids in comparison to experimental data. We focus particularly on the performance of solvation methods for different protonation states, as biological molecules occur almost exclusively in aqueous solution and are subject to protonation with pH. We identify significant shortcomings of current implicit solvent models and present a hybrid solvation approach that improves agreement with experimental data by taking into account the presence of direct interactions between amino acid protonation state and water molecules.

Sensitivity of 3He NMR Parameters to the Proximity of Heterocyclic Rings. The Helium–Furan Dimer

The potential of the NMR technique in probing the interaction between a noble gas atom (used as a spectroscopic probe) and five-membered heterocyclic systems was theoretically studied. Selected 3He NMR parameters of the free noble gas atom and its dimer were calculated using density functional theory with three functionals (M06-2X, CAM-B3LYP, and APFD), combined with dedicated polarization-consistent basis sets, and critically compared with theoretical high level calculations from the literature. The change of isotropic nuclear magnetic shielding of the free noble gas atom in the presence of furan (3He placed above the middle of the ring and used as NMR probe) increased with its proximity to the heterocyclic ring.

Estimation of isotropic nuclear magnetic shieldings in the CCSD(T) and MP2 complete basis set limit using affordable correlation calculations

A linear correlation between isotropic nuclear magnetic shielding constants for seven model molecules (CH2 O, H2 O, HF, F2 , HCN, SiH4 and H2 S) calculated with 37 methods (34 density functionals, RHF, MP2 and CCSD(T)), with affordable pcS-2 basis set and corresponding complete basis set results, estimated from calculations with the family of polarization-consistent pcS-n basis sets is reported. This dependence was also supported by inspection of profiles of deviation between CBS estimated nuclear shieldings and shieldings obtained with the significantly smaller basis sets pcS-2 and aug-cc-pVTZ-J for the selected set of 37 calculation methods. It was possible to formulate a practical approach of estimating the values of isotropic nuclear magnetic shielding constants at the CCSD(T)/CBS and MP2/CBS levels from affordable CCSD(T)/pcS-2, MP2/pcS-2 and DFT/CBS calculations with pcS-n basis sets. The proposed method leads to a fairly accurate estimation of nuclear magnetic shieldings and considerable saving of computational efforts.

NMR characterization of substituent effects in cation–π interactions

Substituent effects in C6H5X…M+ (M+=Li+, Na+, K+, and NH4+) systems have been characterized using isotropic nuclear magnetic shielding constants (σ). Good linear relationship between σ and the Hammett substituent constant σp is established which suggests that the substituent effects can be quantified by NMR experiments. Irrespective of the electronic nature of X, meta carbons, and all the hydrogens are always deshielded and this property is useful for detecting cation–π interaction. The σ-scan plots revealed that the deshielding effect of the cation on the aromatic ring is significant even at large cation–π distance (∼4.5Å).

NMR nuclear magnetic shielding anisotropy of linear molecules within the linear response within the elimination of the small component approach

The influence of the spin-Zeeman (SZ) operator in the evaluation of the spin-orbit effect on the nuclear magnetic shielding tensor in the context of the linear response within the elimination of the small component approach is critically discussed. It is shown that such term yields no contribution to the isotropic nuclear magnetic shielding constant, but it may be of great importance in the determination of individual tensor components, and particularly of the tensor anisotropy. In particular, an interesting relation between the SZ and orbital Zeeman contributions to the spin-orbit effect for the case of linear molecules is shown to hold. Numerical examples for the BrH, IH, and XeF(2) molecules are presented which show that, provided the SZ term is taken into account, results of the individual shielding tensor components and the tensor anisotropy are in good agreement with those obtained by other theoretical methods, and particularly by the Dirac-Hartree-Fock approach.

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: A sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7+/-0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm.

Nuclear magnetic shielding constants of liquid water: Insights from hybrid quantum mechanics/molecular mechanics models

We present a gauge-origin independent method for the calculation of nuclear magnetic shielding tensors of molecules in a structured and polarizable environment. The method is based on a combination of density functional theory (DFT) or Hartree-Fock wave functions with molecular mechanics. The method is unique in the sense that it includes three important properties that need to be fulfilled in accurate calculations of nuclear magnetic shielding constants: (i) the model includes electron correlation effects, (ii) the model uses gauge-including atomic orbitals to give gauge-origin independent results, and (iii) the effect of the environment is treated self-consistently using a discrete reaction-field methodology. The authors present sample calculations of the isotropic nuclear magnetic shielding constants of liquid water based on a large number of solute-solvent configurations derived from molecular dynamics simulations employing potentials which treat solvent polarization either explicitly or implicitly. For both the (17)O and (1)H isotropic shielding constants the best predicted results compare fairly well with the experimental data, i.e., they reproduce the experimental solvent shifts to within 4 ppm for the (17)O shielding and 1 ppm for the (1)H shielding.

Author index

The GIAO-SCF method for calculating isotropic nuclear magnetic shielding values has been utilized to explain certain features in the 1H-NMR spectrum of 2-methylene-8,8-dimethyl-1,4,6,10-tetraoxaspiro[4.5] decane. Population distributions of the low-energy conformers based on their ab initio energies were used to produce weighting factors for the individual calculated shielding values to calculate the weighted average of the shielding values for a complete set of conformers. The differences in 1H chemical shifts between the hydrogens of the two methyl groups and between the axial and equatorial hydrogens in 2-methylene-8,8-dimethyl-1,4,6,10-tetraoxaspiro[4.5] decane were shown to be due to energy differences between the chair and boat orientations of the six-membered ring and contribution from a twist-boat conformation. Results suggest a hypothesis that intramolecular differences in chemical shift might be calculated to a greater degree of accuracy than chemical shifts calculated relative to a standard.

Table of contents

The GIAO-SCF method for calculating isotropic nuclear magnetic shielding values has been utilized to explain certain features in the 1H-NMR spectrum of 2-methylene-8,8-dimethyl-1,4,6,10-tetraoxaspiro[4.5] decane. Population distributions of the low-energy conformers based on their ab initio energies were used to produce weighting factors for the individual calculated shielding values to calculate the weighted average of the shielding values for a complete set of conformers. The differences in 1H chemical shifts between the hydrogens of the two methyl groups and between the axial and equatorial hydrogens in 2-methylene-8,8-dimethyl-1,4,6,10-tetraoxaspiro[4.5] decane were shown to be due to energy differences between the chair and boat orientations of the six-membered ring and contribution from a twist-boat conformation. Results suggest a hypothesis that intramolecular differences in chemical shift might be calculated to a greater degree of accuracy than chemical shifts calculated relative to a standard.

Semiempirical and ab initio conformational analysis of 2-methylene-8,8-dimethyl-1,4,6,10-tetraoxaspiro[4.5] decane with application of GIAO-SCF methods to NMR spectrum interpretation

The GIAO-SCF method for calculating isotropic nuclear magnetic shielding values has been utilized to explain certain features in the 1H-NMR spectrum of 2-methylene-8,8-dimethyl-1,4,6,10-tetraoxaspiro[4.5] decane. Population distributions of the low-energy conformers based on their ab initio energies were used to produce weighting factors for the individual calculated shielding values to calculate the weighted average of the shielding values for a complete set of conformers. The differences in 1H chemical shifts between the hydrogens of the two methyl groups and between the axial and equatorial hydrogens in 2-methylene-8,8-dimethyl-1,4,6,10-tetraoxaspiro[4.5] decane were shown to be due to energy differences between the chair and boat orientations of the six-membered ring and contribution from a twist-boat conformation. Results suggest a hypothesis that intramolecular differences in chemical shift might be calculated to a greater degree of accuracy than chemical shifts calculated relative to a standard.