Complete basis set B3LYP NMR calculations of CDCl3 solvent's water fine spectral details

Extreme NMR shielding in fluoro-nitrogen cations

The structure and NMR shielding of a set of N-F containing cations is reported to near-quantitative accuracy from extensive ab initio calculations. Currently, the shortest experimentally confirmed N-F bond is 1.2461(10) Å in NNF+, however CCSD(T)-F12b/cc-pVQZ-F12 optimised geometries suggest that even shorter N-F bonds are possible for both monocations (1.236 Å, HNF+) and dications (1.098 Å, NF2+). NMR shielding constants have been calculated in a composite manner with individual components from coupled-cluster expansions up to CCSDTQP and basis sets up to aug-cc-pCV8Z, together with vibrational and relativistic corrections. 15N and 19F NMR chemical shifts correlate well with available experimental data. Extreme 19F chemical shifts are predicted for HNF+ (1628.9 ppm) and NH2F2+ (1298.0 ppm), which are by far the largest 19F chemical shifts ever reported and well outside the known range of +865 ppm (F2O2) to -448 ppm (ClF). The 15N chemical shift of -1283.07 ppm in HNF+ is similarly extreme, being well outside the known range of 15N chemical shifts of -730 to 260 ppm (CH3NO2 reference). This work highlights the application of state-of-the-art theoretical techniques, and provides accurate NMR properties of both isolated and yet unknown N-F cations, which can serve to guide and supplement NMR experimentation.

NMR parameters of FNNF as a test for coupled-cluster methods: CCSDT shielding and CC3 spin-spin coupling

NMR shielding and spin-spin coupling constants of cis and trans isomers of FNNF have been determined to near-quantitative accuracy from ab initio calculations. The FNNF system, containing multiple N-F bonds and fluorine atoms, provides a severe test of computational methods. Coupled-cluster methods were used with large basis sets and complete basis set (CBS) extrapolations of the equilibrium geometry results, with vibrational and relativistic corrections. Shielding constants were calculated with basis sets as large as aug-cc-pCV7Z, together with coupled-cluster expansions up to CCSDT, at the all-electron CCSD(T)/aug-cc-pCVQZ optimized geometries. Spin-spin coupling constants have been determined with specialized versions of the correlation consistent basis sets ccJ-pVXZ, further augmented with diffuse functions. All-electron coupled-cluster methods up to CC3 were applied in these calculations. The results of this work highlight the application of state-of-the-art theoretical techniques, and provide the most accurate NMR properties of FNNF to date, which can serve to guide and supplement NMR experimentation.

Local aromaticity mapping in the vicinity of planar and nonplanar molecules

We report on nucleus-independent magnetic shielding (NICS) scans over the centers of six- and five-membered rings in selected metal phthalocyanines (MPc) and fullerene C60 for more accurate characterization of local aromaticity in these compounds. Detailed tests were conducted on model aromatic molecules including benzene, pyrrole, indole, isoindole, and carbazole and subsequently applied to H₂ Pc, ZnPc, Al(OH)Pc, and CuPc. Similar behavior of three selected magnetic probes, Bq, ³ He, and ⁷ Li⁺ , approaching perpendicularly the ring centers, was observed. For better visualization of shielding zone over the centers of aromatic rings, we introduced a simple mathematical procedure: the first and second derivatives of scan curves with respect to magnetic probe position enabled their additional examination. It allowed an easier localization of curve minimum and discrimination between areas in space varying by the magnetic field magnitude and to illustrate local aromaticity of two different kinds of rings in MPc with better resolution. Our results supported earlier reports on very low aromaticity indexes of pyrrole ring incorporated into MPc and significant aromaticity of the central macrocycle. No direct dependence between harmonic oscillator model of aromaticity and NICS was observed. Instead, a correlation between position of scan curve minimum and its magnitude were observed. In addition, the NICS values and ³ He chemical shifts in the middle of neutral C60 and C60^6- anion agreed well with the reported experimental NMR values for He@C60 and He@C60^6- .

Theoretical calculations of carbon-hydrogen spin-spin coupling constants

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.

Convergence of nuclear magnetic shieldings and one-bond 1 J(11 B1 H) indirect spin-spin coupling constants in small boron molecules

Self-consistent field Hartree-Fock, density functional theory, and coupled-cluster calculations of the nuclear magnetic shielding constants of BH and BH₃ molecules have been conducted to characterize the convergence of individual results obtained with correlation-consistent and polarization-consistent basis sets. The individual ¹¹ B and ¹ H NMR parameters were estimated in the complete basis set limit and compared with benchmark results. Only the KT3 density functional accurately reproduced ¹¹ B shielding in BH molecule.

NMR shielding constants in group 15 trifluorides

By combining large basis and complete basis set (CBS) extrapolations of the coupled-cluster equilibrium geometry results with rovibrational and relativistic corrections, we demonstrate that it is possible to achieve near-quantitative accuracy for the NMR shielding constants in three group 15 trifluorides – NF₃, PF₃ and AsF₃.

Calculation of 15N NMR chemical shifts: Recent advances and perspectives

Recent advances in computation of ¹⁵N 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 ¹⁵N NMR are illustrated for the series of neutral and protonated open-chain nitrogen-containing compounds and nitrogen heterocycles, coordination and intermolecular complexes.

Normal halogen dependence of 13 C NMR chemical shifts of halogenomethanes revisited at the four-component relativistic level

The 'Normal Halogen Dependence' of ¹³ C NMR chemical shifts in the series of halogenomethanes is revisited at the four-component relativistic level. Calculations of ¹³ C NMR chemical shifts of 70 halogenomethanes have been carried out at the density functional theory (DFT) and MP2 levels with taking into account relativistic effects using the four-component relativistic theory of Dirac-Coulomb within the different computational methods (4RPA, 4OPW91) and hybrid computational schemes (MP2 + 4RPA, MP2 + 4OPW91). The most efficient computational protocols are derived for practical purposes. Relativistic shielding effect reaches as much as several hundreds of ppm for heavy halogenomethanes, and to account for this effect in comparison with experiment at the qualitative level, relativistic Dyall's basis sets of triple-zeta quality or higher are to be used within the framework of the four-component relativistic theory taking into account solvent effects. Relativistic geometrical optimization (as compared with the non-relativistic level) is essential for the molecules containing at least two iodines at one carbon atom. Copyright © 2016 John Wiley & Sons, Ltd.

Predicting the structure and vibrational frequencies of ethylene using harmonic and anharmonic approaches at the Kohn-Sham complete basis set limit

In this work, regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn–Sham complete basis set (KS CBS) limit are demonstrated for the first time. The performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with two Pople basis sets (6-311++G** and 6-311++G(3df,2pd)), the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, 6) was tested.

The BLYP-calculated harmonic frequencies were found to be markedly closer than the B3LYP-calculated harmonic frequencies to the experimentally derived values, while the calculated anharmonic frequencies consistently underestimated the observed wavenumbers. The different basis set families gave very similar estimated values for the CBS parameters. The anharmonic frequencies calculated with B3LYP/aug-pc-3 were consistently significantly higher than those obtained with the pc-3 basis set; applying the aug-pcseg-n basis set family alleviated this problem. Utilization of B3LYP/aug-pcseg-n basis sets instead of B3LYP/aug-cc-pVXZ, which is computationally less expensive, is suggested for medium-sized molecules. Harmonic BLYP/pc-2 calculations produced fairly accurate ethylene frequencies.

Graphical Abstract

MP2 calculation of (77) Se NMR chemical shifts taking into account relativistic corrections

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%.

First example of a high-level correlated calculation of the indirect spin-spin coupling constants involving tellurium: tellurophene and divinyl telluride

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.

3He NMR: from free gas to its encapsulation in fullerene

The (3)He nuclear magnetic shieldings were calculated for single helium atom, its dimer, simple models of fullerene cages (He@Cn), and single wall carbon nanotubes. The performances of several levels of theory (HF, MP2, DFT-VSXC, CCSD, CCSD(T), and CCSDT) were tested. Two sets of polarization-consistent basis sets were used (pcS-n and aug-pcS-n), and an estimate of (3)He nuclear magnetic shieldings in the complete basis set limit using a two-parameter fit was established. Theoretical (3)He results reproduced accurately previously reported theoretical values for helium gas, dimer, and helium probe inside several fullerene cages. Excellent agreement with experimental values was achieved. (3)He nuclear magnetic shieldings of single helium atom approaching various points of benzene ring were tested, and an impact of (3)He confinement within fullerene cages of different size on the (3)He chemical shift was determined.

From CCSD(T)/aug-cc-pVTZ-J to CCSD(T) complete basis set limit isotropic nuclear magnetic shieldings via affordable DFT/CBS calculations

It is shown that a linear correlation exists between nuclear shielding constants for nine small inorganic and organic molecules (N(2), CO, CO(2), NH(3), CH(4), C(2)H(2), C(2)H(4), C(2)H(6) and C(6)H(6)) calculated with 47 methods (42 DFT methods, RHF, MP2, SOPPA, SOPPA(CCSD), CCSD(T)) and the aug-cc-pVTZ-J basis set and corresponding complete basis set results, estimated from calculations with the family of polarization-consistent pcS-n basis sets. This implies that the remaining basis set error of the aug-cc-pVTZ-J basis set is very similar in DFT and CCSD(T) calculations. As the aug-cc-pVTZ-J basis set is significantly smaller, CCSD(T)/aug-cc-pVTZ-J calculations allow in combination with affordable DFT/pcS-n complete basis set calculations the prediction of nuclear shieldings at the CCSD(T) level of nearly similar accuracy as those, obtained by fitting results obtained from computationally demanding pcS-n calculations at the CCSD(T) limit. A significant saving of computational efforts can thus be achieved by scaling inexpensive CCSD(T)/aug-cc-pVTZ-J calculations of nuclear isotropic shieldings with affordable DFT complete basis set limit corrections.

H2O, H2, HF, F2 and F2O nuclear magnetic shielding constants and indirect nuclear spin-spin coupling constants (SSCCs) in the BHandH/pcJ-n and BHandH/XZP Kohn-Sham limits

Good performance of segmented contracted basis sets XZP, where X = D, T, Q and 5, for obtaining H(2)O, H(2), HF, F(2) and F(2)O nuclear isotropic shielding constants in the BHandH Kohn-Sham basis set limit was shown. The results of two- and three-parameter complete basis set limit extrapolation schemes were compared with experimental results, earlier literature data and benchmark ab initio results. Similar convergence patterns of shieldings obtained from calculations using general purpose XZP basis sets and from polarization-consistent basis sets pcS-n and pcJ-n, where n = 0, 1, 2, 3 and 4, designed to accurately predict magnetic properties were observed. On the contrary, the SSCCs were more sensitive to the XZP basis set size and generally less accurate than those estimated using pcJ-n basis set family. The BHandH density functional markedly outperforms B3LYP method in predicting heavy atom shieldings and SSCCs values in the studied systems.

Complete basis set prediction of methanol isotropic nuclear magnetic shieldings and indirect nuclear spin-spin coupling constants (SSCC) using polarization-consistent and XZP basis sets and B3LYP and BHandH density functionals

Efficient B3LYP and BHandH density functionals were used to estimate methanol's nuclear magnetic isotropic shieldings and spin-spin coupling constants in the basis set limit. Polarization-consistent pcS-n and pcJ-n (n = 0, 1, 2, 3 and 4), and segmented contracted XZP, where X = D, T, Q and 5, basis sets were used and the results fitted with simple mathematical formulas. The performance of the methods was assessed from comparison with experiment and higher level calculations. 1J(CH) and 3J(HH) values were determined from very diluted solutions in deuterochloroform and compared with theoretical predictions. The agreement between complete basis set (CBS) density functional theory (DFT) predicted isotropic shieldings and spin-spin values and experiment was good. The BHandH/pcS-n methanol shieldings obtained using structures optimized at the same level of theory are approaching the accuracy of the advanced coupled-cluster-singles-doubles-approximate triples (CCSD(T)) calculations.

Prediction of water's isotropic nuclear shieldings and indirect nuclear spin-spin coupling constants (SSCCs) using correlation-consistent and polarization-consistent basis sets in the Kohn-Sham basis set limit

Density functional theory (DFT) was used to estimate water's isotropic nuclear shieldings and indirect nuclear spin-spin coupling constants (SSCCs) in the Kohn-Sham (KS) complete basis set (CBS) limit. Correlation-consistent cc-pVxZ and cc-pCVxZ (x = D, T, Q, 5, and 6), and their modified versions (ccJ-pVxZ, unc-ccJ-pVxZ, and aug-cc-pVTZ-J) and polarization-consistent pc-n and pcJ-n (n = 0, 1, 2, 3, and 4) basis sets were used, and the results fitted with a simple mathematical formula. The performance of over 20 studied density functionals was assessed from comparison with the experiment. The agreement between the CBS DFT-predicted isotropic shieldings, spin-spin values, and the experimental values was good and similar for the modified correlation-consistent and polarization-consistent basis sets. The BHandH method predicted the most accurate (1)H, (17)O isotropic shieldings and (1)J(OH) coupling constant (deviations from experiment of about -0.2 and -1 ppm and 0.6 Hz, respectively). The performance of BHandH for predicting water isotropic shieldings and (1)J(OH) is similar to the more advanced methods, second-order polarization propagator approximation (SOPPA) and SOPPA(CCSD), in the basis set limit.