Ab initio and relativistic DFT study of spin–rotation and NMR shielding constants in XF6 molecules, X = S, Se, Te, Mo, and W

Reducing Exact Two-Component Theory for NMR Couplings to a One-Component Approach: Efficiency and Accuracy

The self-consistent and complex spin-orbit exact two-component (X2C) formalism for NMR spin-spin coupling constants [ J. Chem. Theory Comput. 17, 2021, 3874-3994] is reduced to a scalar one-component ansatz. This way, the first-order response term can be partitioned into the Fermi-contact (FC) and spin-dipole (SD) interactions as well as the paramagnetic spin-orbit (PSO) contribution. The FC+SD terms are real and symmetric, while the PSO term is purely imaginary and antisymmetric. The relativistic one-component approach is combined with a modern density functional treatment up to local hybrid functionals including the response of the current density. Computational demands are reduced by factors of 8-24 as shown for a large tin compound consisting of 137 atoms. Limitations of the current ansatz are critically assessed for Sn, Pb, Pd, and Pt compounds, i.e. the one-component treatment is not sufficient for tin compounds featuring a few heavy halogen atoms.

A Computational Study on Closed-Shell Molecular Hexafluorides MF6 (M=S, Se, Te, Po, Xe, Rn, Cr, Mo, W, U) - Molecular Structure, Anharmonic Frequency Calculations, and Prediction of the NdF6 Molecule

Quantum chemical methods were used to study the molecular structure and anharmonic IR spectra of the experimentally known closed-shell molecular hexafluorides MF₆ (M=S, Se, Te, Xe, Mo, W, U). First, the molecular structures and harmonic frequencies were investigated using Density Functional Theory (DFT) with all-electron basis sets and explicitly considering the influence of spin-orbit coupling. Second, anharmonic frequencies and IR intensities were calculated with the CCSD(T) coupled cluster method and compared, where available, with IR spectra recorded by us. These comparisons showed satisfactory results. The anharmonic IR spectra provide means for identifying experimentally too little studied or unknown MF₆ molecules with M=Cr, Po, Rn. To the best of our knowledge, we predict the NdF₆ molecule for the first time and show it to be a true local minimum on the potential energy surface. We used intrinsic bond orbital (IBO) analyses to characterize the bonding situation in comparison with the UF₆ molecule.

Quantum chemical calculations of 77 Se and 125 Te nuclear magnetic resonance spectral parameters and their structural applications

An accurate quantum chemical (QC) modeling of ⁷⁷ Se and ¹²⁵ Te nuclear magnetic resonance (NMR) spectra is deeply involved in the NMR structural assignment for selenium and tellurium compounds that are of utmost importance both in organic and inorganic chemistry nowadays due to their huge application potential in many fields, like biology, medicine, and metallurgy. The main interest of this review is focused on the progress in QC computations of ⁷⁷ Se and ¹²⁵ Te NMR chemical shifts and indirect spin-spin coupling constants involving these nuclei. Different computational methodologies that have been used to simulate the NMR spectra of selenium and tellurium compounds since the middle of the 1990s are discussed with a strong emphasis on their accuracy. A special accent is placed on the calculations resorting to the relativistic methodologies, because taking into account the relativistic effects appreciably influences the precision of NMR calculations of selenium and, especially, tellurium compounds. Stereochemical applications of quantum chemical calculations of ⁷⁷ Se and ¹²⁵ Te NMR parameters are discussed so as to exemplify the importance of integrated approach of experimental and computational NMR techniques.

Correlated ab initio calculations of one-bond 31 P77 Se and 31 P125 Te spin-spin coupling constants in a series of PSe and PTe systems accounting for relativistic effects (part 2)

Synthetic chalcogen-phosphorus chemistry permanently makes new challenges to computational Nuclear Magnetic Resonance (NMR) spectroscopy, which has proven to be a powerful tool of structural analysis of chalcogen-phosphorus compounds. This paper reports on the calculations of one-bond ³¹ P⁷⁷ Se and ³¹ P¹²⁵ Te NMR spin-spin coupling constants (SSCCs) in the series of phosphine selenides and tellurides. The applicability of the combined computational approach to the one-bond ³¹ P⁷⁷ Se and ³¹ P¹²⁵ Te SSCCs, incorporating the composite nonrelativistic scheme, built of high-accuracy correlated SOPPA (CC2) and Coupled Cluster Single and Double (CCSD) methods and the Density Functional Theory (DFT) relativistic corrections (four-component level), was examined against the experiment and another scheme based on the four-component relativistic DFT method. A special J-oriented basis set (acv3z-J) for selenium and tellurium atoms, developed previously by the authors, was used throughout the NMR calculations in this work at the first time. The proposed computational methodologies (combined and 'pure') provided a reasonable accuracy for ³¹ P⁷⁷ Se and ³¹ P¹²⁵ Te SSCCs against experimental data, characterizing by the mean absolute percentage errors of about 4% and 1%, and 12% and 8% for selenium and tellurium species, respectively. The present study reports typical relativistic corrections to ⁷⁷ Se³¹ P and ¹²⁵ Te³¹ P SSCCs, calculated within the four-component DFT formalism for a broad series of tertiary phosphine selenides and tellurides with different substituents at phosphorus.

NMR absolute shielding scales and nuclear magnetic dipole moments of transition metal nuclei

This work reports new, accurate nuclear magnetic dipole moments for transition metal nuclei where the long-standing systematic error due to obsolete diamagnetic correction has been eliminated by ab initio calculations of NMR shielding constants.

Electronic Structure Predictions of the Energetic Properties of Tellurium Fluorides

The heats of formation, bond dissociation energies (BDEs), fluoride affinities (FA), fluorocation affinities (FCA), electron affinities (EA), and ionization energies (IP) of TeF _n ( n = 1-6) have been predicted using the Feller-Peterson-Dixon (FPD) approach. To benchmark the approach, the bond dissociation energies of Te₂ and TeO, the heats of formation of Te₂, TeH₂, TeO, and TeO₂, and the electron affinity for TeO and TeO₂ were calculated as there are experimental thermodynamic data available for these tellurium compounds, which allow confirmation of the heat of formation of Te gas as Δ H_f,0K(Te) = 50.7 ± 0.6 kcal/mol. Spin-orbit corrections are required for good results and cannot be ignored. A comparison among fluoride affinities, fluorocation affinities, electron affinities, and ionization energies of TeF _n and SeF _n is reported.

Relativistic heavy atom effect on 13 C NMR chemical shifts initiated by adjacent multiple chalcogens

In this paper, we have investigated the cumulative peculiarity of the "heavy atom on light atom" effect on the ¹³ C NMR chemical shifts, initiated by the adjacent chalcogens. For this purpose, the most accurate hybrid computational scheme for the calculation of chemical shifts of carbon nuclei, directly bonded with several heavy chalcogens, is introduced and attested on the representative series of molecules. The best hybrid scheme combines the nonrelativistic coupled cluster-based approach with the different types of corrections, including vibrational, solvent, and relativistic. The dependences of the total relativistic corrections to carbon shielding constants in 2 series of model compounds, namely, X═¹³ C═Y (X, Y = O, S, Se, Te) and C(XH)_m (YH)_n (ZH)_p (QH)_s H_1-m H_1-n H_1-p H_1-s (X, Y, Z, Q = S, Se, Te and m, n, p, s = 0, 1), on the total atomic number of the adjacent chalcogens have been obtained.

Prediction of low-field nuclear singlet lifetimes with molecular dynamics and quantum-chemical property surface

Molecular dynamics and quantum chemistry methods are implemented to quantify nuclear spin-1/2 pair singlet-state relaxation rates. Illustrated is the relevant spin-internal-motion mechanism (SIM).

Relativistic Calculations of Nuclear Magnetic Resonance Parameters

Relativistic effects are important for the accurate evaluation of the observables of nuclear magnetic resonance (NMR) spectroscopy, the nuclear magnetic shielding and the indirect spin–spin coupling tensors. Some of the most notable relativistic effects, in particular for light elements in the vicinity of heavy nuclei, are due to spin–orbit effects, an effect difficult to evaluate when starting from a non-relativistic wavefunction. Two- and four-component relativistic methods include spin–orbit effects variationally, and the recent improvements in the computational efficiency of these methods open new opportunities for accurate calculations of NMR parameters also for molecules with heavy elements. We here present an overview of the different approximations that have been introduced for calculating relativistic effects with two- and four-component methods and how these methods can be used to calculate the NMR parameters. We will also give some examples of systems that have been studied computationally with two- and four-component relativistic methods and discuss the importance of relativistic effects on the shielding and indirect spin–spin coupling constants.

Theoretical analysis of NMR shieldings in XSe and XTe (X = Si, Ge, Sn and Pb): the spin-rotation constant saga

How the electronic contribution to the spin-rotation constant is close to the paramagnetic contribution of the NMR absolute shielding constant?

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

Spin-rotation and NMR shielding constants in XF molecules (X = B, Al, Ga, In, and Tl)

Accurate spin-rotation and absolute shielding constants in a series of XF molecules (X = (11)B, (27)Al, (69)Ga, (115)In, and (205)Tl) determined using high-level ab initio coupled-cluster and four-component relativistic density-functional theory (DFT) calculations are presented. The accuracy of the results is established by comparing the relativistically and vibrationally corrected calculated values with available experimental data; for spin-rotation and shielding constants for which no experimental data exist, we provide new and reliable values. For both properties, our results can be considered as reference values against which more approximate quantum-chemical methods can be benchmarked.