Nuclear spin–spin coupling constants from regular approximate relativistic density functional calculations. II. Spin–orbit coupling effects and anisotropies

Giant spin–orbit effects on 1H and 13C NMR shifts for uranium(vi) complexes revisited: role of the exchange–correlation response kernel, bonding analyses, and new predictions

Visiting the previously predicted giant spin–orbit-induced ¹H and ¹³C shifts in U(vi) complexes with improved methodology retains the reported unusual shift ranges, provides better understanding of the observations and gives improved confidence in the predictions.

[P3Se4](+): A Binary Phosphorus-Selenium Cation

Although a fairly large number of binary group 15/16 element cations have been reported, no example involving phosphorus in combination with a group 16 element has been synthesized and characterized to date. In this contribution is reported the synthesis and structural characterization of the first example of such a cation, namely a nortricyclane-type [P3Se4](+). This cation has been independently discovered by three groups through three different synthetic routes, as described herein. The molecular and electronic structure of the [P3Se4](+) cage and its crystal properties in the solid state have been characterized comprehensively by using X-ray diffraction, Raman, and nuclear magnetic resonance spectroscopies, as well as quantum chemical calculations.

The angular overlap model extended for two-open-shell f and d electrons

We discuss the applicability of the Angular Overlap Model (AOM) to evaluate the electronic structure of lanthanide compounds, which are currently the subject of incredible interest in the field of luminescent materials. The functioning of phosphors is well established by the f-d transitions, which requires the investigation of both the ground 4f(n) and excited 4f(n-1)5d(1) electron configurations of the lanthanides. The computational approach to the problem is based on the effective Hamiltonian adjusted from ligand field theory, but not restricted to it. The AOM parameterization implies the chemical bonding concept. Focusing our interest on this interaction, we take the advantages offered by modern computational tools to extract AOM parameters, which ensure the transparency of the theoretical determination and convey chemical intuitiveness of the non-empirical results. The given model contributes to the understanding of lanthanides in modern phosphors with high or low site symmetry and presents a non-empirical approach using a less sophisticated computational procedure for the rather complex problem of the ligand field of both 4f and 5d open shells.

Interpretation of the longitudinal (13)C nuclear spin relaxation and chemical shift data for five bromoazaheterocycles supported by nonrelativistic and relativistic DFT calculations

The longitudinal relaxation times of (13)C nuclei and NOE enhancement factors for 2-bromopyridine (1), 6-bromo-9-methylpurine (2), 3,5-dibromopyridine (3), 2,4-dibromopyrimidine (4), and 2,4,6-tribromopyrimidine (5) have been measured at 25 °C and B0 = 11.7 T. The most important contributions to the overall relaxation rates of nonbrominated carbons, i.e., the relaxation rates due to the (13)C-(1)H dipolar interactions and the shielding anisotropy mechanism, have been separated out. For 3 and 5, additionally, the T2,Q((14)N) values have been established from (14)N NMR line widths. All of these data have been used to determine rotational diffusion tensors for the investigated molecules. The measured saturation recovery curves of brominated carbons have been decomposed into two components to yield relaxation times, which after proper corrections provided parameters characterizing the scalar relaxation of the second kind for (13)C nuclei of (79)Br- and (81)Br-bonded carbons. These parameters and theoretically calculated quadrupole coupling constants for bromine nuclei have allowed the values of one-bond (13)C-(79)Br spin-spin coupling constants to be calculated. Independently, the coupling constants and magnetic shielding constants of the carbon nuclei have been calculated theoretically using the nonrelativistic and relativistic DFT methods F/6-311++G(2d,p)/PCM and so-ZORA/F/TZ2P/COSMO (F = BHandH or B3LYP), respectively. The agreement between the experimental and theoretical values of these parameters is remarkably dependent on the theoretical method used.

What factors influence the metal-proton spin-spin coupling constants in mercury- and cadmium-substutited rubredoxin?

The indirect metal-proton spin-spin coupling constants between protons in cysteine groups and the mercury or cadmium nucleus have been calculated for a small model of Me-rubredoxin complex (Me = Cd, Hg) by means of density functional theory with zeroth-order regular approximation Hamiltonian (DFT-ZORA). The calculated spin-spin coupling constants, in spite of the moderate size of the model system, are in good agreement with the values measured in NMR experiment, which are in the 0.29-0.56 Hz range for the Cd complex and in the 0.57-2.20 Hz range for the Hg complex. The robustness of the chosen method has been verified by calculations with a number of different exchange-correlation functionals and basis sets. Additionally, it has been shown that the short- and long-distance metal-proton coupling constants are affected mainly by the values of the metal-proton distance and the H-N-C-C dihedral angle.

Relativistic calculations of magnetic resonance parameters: background and some recent developments

This article outlines some basic concepts of relativistic quantum chemistry and recent developments of relativistic methods for the calculation of the molecular properties that define the basic parameters of magnetic resonance spectroscopic techniques, i.e. nuclear magnetic resonance shielding, indirect nuclear spin-spin coupling and electric field gradients (nuclear quadrupole coupling), as well as with electron paramagnetic resonance g-factors and electron-nucleus hyperfine coupling. Density functional theory (DFT) has been very successful in molecular property calculations, despite a number of problems related to approximations in the functionals. In particular, for heavy-element systems, the large electron count and the need for a relativistic treatment often render the application of correlated wave function ab initio methods impracticable. Selected applications of DFT in relativistic calculation of magnetic resonance parameters are reviewed.

Quantum-chemical analyses of aromaticity, UV spectra, and NMR chemical shifts in plumbacyclopentadienylidenes stabilized by Lewis bases

We carried out a series of zeroth-order regular approximation (ZORA)-density functional theory (DFT) and ZORA-time-dependent (TD)-DFT calculations for molecular geometries, NMR chemical shifts, nucleus-independent chemical shifts (NICS), and electronic transition energies of plumbacyclopentadienylidenes stabilized by several Lewis bases, (Ph)2 ((t) BuMe2 Si)2 C4 PbL1 L2 (L1, L2 = tetrahydrofuran, Pyridine, N-heterocyclic carbene), and their model molecules. We mainly discussed the Lewis-base effect on the aromaticity of these complexes. The NICS was used to examine the aromaticity. The NICS values showed that the aromaticity of these complexes increases when the donation from the Lewis bases to Pb becomes large. This trend seems to be reasonable when the 4n-Huckel rule is applied to the fractional π-electron number. The calculated (13)C- and (207)Pb-NMR chemical shifts and the calculated UV transition energies reasonably reproduced the experimental trends. We found a specific relationship between the (13)C-NMR chemical shifts and the transition energies. As we expected, the relativistic effect was essential to reproduce a trend not only in the (207)Pb-NMR chemical shifts and J[Pb-C] but also in the (13)C-NMR chemical shifts of carbons adjacent to the lead atom.

Formation, structural characterization, and reactions of a unique cyclotrimeric vicinal Lewis pair containing (C6F5)2P-Lewis base and (C6F5)BH-Lewis acid components

A unique vicinal Lewis pair cyclotrimer has been synthesized and characterized with respect to its structure and reactivity.

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

Donor-free phosphenium-metal(0)-halides with unsymmetrically bridging phosphenium ligands

Reactions of (cod)MCl2 (cod = 1,5 cyclooctadiene, M = Pd, Pt) with N-heterocyclic secondary phosphines or diphosphines produced complexes [(NHP)MCl]2 (NHP = N-heterocyclic phosphenium). The Pd complex was also accessible from a chlorophosphine precursor and Pd2(dba)3. Single-crystal X-ray diffraction studies established the presence of dinuclear complexes that contain μ-bridging NHP ligands in an unsymmetrical binding mode and display a surprising change in metal coordination geometry from distorted trigonal (M = Pd) to T-shaped (M = Pt). DFT calculations on model compounds reproduced these structural features for the Pt complex but predicted an unusual C2v-symmetric molecular structure with two different metal coordination environments for the Pd species. The deviation between this structure and the actual centrosymmetric geometry is accounted for by the prediction of a flat energy hypersurface, which permits large distortions in the orientation of the NHP ligands at very low energetic cost. The DFT results and spectroscopic studies suggest that the title compounds should be described as phosphenium-metal(0)-halides rather than conventional phosphido complexes of divalent metal cations and indicate that the NHP ligands receive net charge donation from the metals but retain a distinct cationic character. The unsymmetric NHP binding mode is associated with an unequal distribution of σ-donor/π-acceptor contributions in the two M-P bonds. Preliminary studies indicate that reactions of the Pd complex with phosphine donors provide a viable source of ligand-stabilized, zerovalent metal atoms and metal(0)-halide fragments.

Computational study and molecular orbital analysis of NMR shielding, spin-spin coupling, and electric field gradients of azido platinum complexes

(195)Pt, (14)N, and (15)N NMR data for five azido (N(3)(-)) complexes are studied using relativistic density functional theory (DFT). Good agreement with experiment is obtained for Pt and N chemical shifts as well as Pt-N J-coupling constants. Calculated (14)N electric field gradients (EFGs) reflect experimentally observed trends for the line broadening of azido (14)N NMR signals. A localized molecular orbital analysis of the nitrogen EFGs and chemical shifts is performed to explain some interesting trends seen experimentally and in the first-principles calculations: (i) (14)N NMR signals for the Pt-coordinating (N(α)) nuclei in the azido ligands are much broader than for the central (N(β)) or terminal (N(γ)) atoms. The N(β) signals are particularly narrow; (ii) compared to N(γ), the N(α) nuclei are particularly strongly shielded; (iii) N(β) nuclei have much larger chemical shifts than N(α) and N(γ) ; and (iv) The Pt-N(α) J-coupling constants are small in magnitude when considering the formal sp hybridization of N(α). It is found that for N(α) a significant shielding reduction due to formation of the dative N(α)-Pt bond is counterbalanced by an increased shielding from spin-orbit (SO) coupling originating at Pt. Upon coordination, the strongly delocalized π system of free azide localizes somewhat on N(β) and N(γ). This effect, along with rehybridization at N(α) upon bond formation with Pt, is shown to cause a deshielding of N(γ) relative to N(α) and a strong increase of the EFG at N(α). The large 2p character of the azide σ bonds is responsible for the particularly high N(β) chemical shifts. The nitrogen s-character of the Pt-N(α) bond is low, which is the reason for the small J-coupling. Similar bonding situations are likely to be found in azide complexes with other transition metals.

Observation of scalar nuclear spin–spin coupling in van der Waals complexes

Scalar couplings between covalently bound nuclear spins are a ubiquitous feature in nuclear magnetic resonance (NMR) experiments, imparting valuable information to NMR spectra regarding molecular structure and conformation. Such couplings arise due to a second-order hyperfine interaction, and, in principle, the same mechanism should lead to scalar couplings between nuclear spins in unbound van der Waals complexes. Here, we report the first observation of scalar couplings between nuclei in van der Waals complexes. Our measurements are performed in a solution of hyperpolarized 129Xe and pentane, using superconducting quantum interference devices to detect NMR in 10 mG fields, and are in good agreement with calculations based on density functional theory. van der Waals forces play an important role in many physical phenomena. The techniques presented here may provide a new method for probing such interactions.

1H and 13C{1H} NMR spectral parameters of sulfur mustards, nitrogen mustards, and lewisites: computing and predicting of reference spectra for chemical identification

The (1)H and (13)C{(1)H} chemical shifts and (1)H spin-spin couplings of sulfur mustards, nitrogen mustards, and lewisites scheduled in the Chemical Weapons Convention, and those of bis(2-chloromethyl)disulfide, were determined in CDCl(3), CD(2)Cl(2), and (CD(3))(2)CO. Accurate parameters of this kind of series can be used for evaluating the current molecular modeling programs and the chemical shift and coupling constant prediction possibilities of the programs. Several prediction tests were made with commercial programs, and the results are reported here.

Residual dipolar coupling between quadrupolar nuclei under magic-angle spinning and double-rotation conditions

Residual dipolar couplings between spin-1/2 and quadrupolar nuclei are often observed and exploited in the magic-angle spinning (MAS) NMR spectra of spin-1/2 nuclei. These orientation-dependent splittings contain information on the dipolar interaction, which can be translated into structural information. The same type of splittings may also be observed for pairs of quadrupolar nuclei, although information is often difficult to extract from the quadrupolar-broadened lineshapes. Here, the complete theory for describing the dipolar coupling between two quadrupolar nuclei in the frequency domain by Hamiltonian diagonalization is given. The theory is developed under MAS and double-rotation (DOR) conditions, and is valid for any spin quantum numbers, quadrupolar coupling constants, asymmetry parameters, and tensor orientations at both nuclei. All terms in the dipolar Hamiltonian become partially secular and contribute to the NMR spectrum. The theory is validated using experimental 11B and 35/37Cl NMR experiments carried out on powdered B-chlorocatecholborane, where both MAS and DOR are used to help separate effects of the quadrupolar interaction from those of the dipolar interaction. It is shown that the lineshapes are sensitive to the quadrupolar coupling constant of both nuclei and to the J coupling (including its sign). From these experiments, the dipolar coupling constant for a heteronuclear spin pair of quadrupolar nuclei may be obtained as well as the sign of the quadrupolar coupling constant of the perturbing nucleus; these are two parameters that are difficult to obtain experimentally otherwise.

[Pt@Pb12]2– — A challenging system for relativistic density functional theory calculations of 195Pt and 207Pb NMR parameters

We report computations of NMR chemical shifts and indirect spin-spin coupling constants (J couplings) for the [Pt@Pb12]2– “superatom”. The system is strongly influenced by relativistic effects. The Pt–Pb coupling constant is predicted to be negative, with its magnitude being in reasonable agreement with experiment. Pt and Pb chemical shifts also agree reasonably well with experiment. The Pb shielding tensor is strongly anisotropic, with a large deshielding principal component dominated by magnetic coupling between frontier orbitals of the cluster that resemble atomic g orbitals. The NMR parameters are sensitive to approximations made in the computations and require the inclusion of spin-orbit coupling in the theoretical model to achieve reliable results. Computing the NMR parameters of the compact [Pt@Pb12]2– system with its many electrons proves to be a challenging test case for relativistic density functional methods.

Conformational analysis and diastereotopic assignments in the series of selenium-containing heterocycles by means of 77Se-1H spin-spin coupling constants: a combined theoretical and experimental study

A combined theoretical and experimental study on the stereochemical behavior of (77)Se-(1)H spin-spin coupling constants has been performed at the second-order polarization propagator approach level together with heteronuclear multiple-bond correlation technique in the series of selenium-containing four-, five- and six-membered heterocycles including the derivatives of thiaselenetane, selenasilole, thiaselenole, thiaselenolane and dihydrothiaselenine. Geminal and vicinal (77)Se-(1)H spin-spin couplings were shown to have the pronounced stereochemical dependences in respect with the topology of the coupling pathway, internal rotation of the side-chain substituents and ring inversion providing a straightforward tool for the conformational analysis and diastereotopic assignments in the chiral organoselenium compounds.

A ZORA-DFT and NLMO study of the one-bond fluorine–X indirect nuclear spin-spin coupling tensors for various VSEPR geometries

Zeroth-order regular approximation (ZORA) density functional theory (DFT) calculations of one-bond X–19F indirect nuclear spin-spin coupling (J) tensors were performed on a series of fluorine-containing compounds covering several valence shell electron pair repulsion (VSEPR) theory geometries for which J, by symmetry, is not required to be axially symmetric. The calculations show that the antisymmetric components of J are only of the same order of magnitude as the principal components of the symmetric J-coupling tensor for a few geometries, and that in cases of approximate axial symmetry along the bond, J remains nearly axially symmetric with its unique component along the bond. In general, different species having the same nominal geometry tend to have similar tensor orientations, magnitudes of anisotropy of J relative to the isotropic coupling constant, as well as the same dominant contributions from the different coupling mechanisms. Structures are also systematically modified to determine how the tensor components depend on geometrical parameters. The isotropic coupling constants are subsequently interpreted using a natural localized molecular orbital (NLMO) approach. Our results could prove to be useful for future experimental characterizations of J tensors in systems having symmetry properties that do not force J to be axially symmetric or coincident with the dipolar coupling tensor.

NMR J-coupling constants in cisplatin derivatives studied by molecular dynamics and relativistic DFT

Solvent effects on J((195)Pt-(15)N) one-bond nuclear spin-spin coupling constants (J(PtN)) of cisplatin [cis-diamminedichloroplatinum(II)] and three cisplatin derivatives are investigated using a combination of density functional theory (DFT) based ab initio molecular dynamics (aiMD) and all-electron relativistic DFT NMR calculations employing the two-component relativistic zeroth-order regular approximation (ZORA). Good agreement with experiment is obtained when explicit solvent molecules are considered and when the computations are performed with a hybrid functional. Spin-orbit coupling causes only small effects on J(PtN) . Key factors contributing to the magnitude of coupling constants are elucidated, with the most significant being the presence of solvent as well as the quality of the density functional and basis set combination. The solvent effects are of the same magnitude as J(PtN) calculated for gas-phase geometries. However, the trends of J(PtN) among the complexes are already present in the gas phase. Results obtained with a continuum solvent model agree quite well with the aiMD results, provided that the Pt solvent-accessible radius is carefully chosen. The aiMD results support the existence of a partial hydrogen-bond-like inverse-hydration-type interaction affording a weak (1)J(Pt⋅⋅⋅H(w)) coupling between the complexes and the coordinating water molecule.

Solid-state (55)Mn NMR spectroscopy of bis(μ-oxo)dimanganese(IV) [Mn(2)O(2)(salpn)(2)], a model for the oxygen evolving complex in photosystem II

We have examined the antiferromagneticly coupled bis(μ-oxo)dimanganese(IV) complex [Mn(2)O(2)(salpn)(2)] (1) with (55)Mn solid-state NMR at cryogenic temperatures and first-principle theory. The extracted values of the (55)Mn quadrupole coupling constant, C(Q), and its asymmetry parameter, η(Q), for 1 are 24.7 MHz and 0.43, respectively. Further, there was a large anisotropic contribution to the shielding of each Mn(4+), i.e. a Δσ of 3375 ppm. Utilizing broken symmetry density functional theory, the predicted values of the electric field gradient (EFG) or equivalently the C(Q) and η(Q) at ZORA, PBE QZ4P all electron level of theory are 23.4 MHz and 0.68, respectively, in good agreement with experimental observations.

Structural trends of 77Se-1H spin-spin coupling constants and conformational behavior of 2-substituted selenophenes

Experimental measurements and second-order polarization propagator approach (SOPPA) calculations of (77)Se-(1)H spin-spin coupling constants together with theoretical energy-based conformational analysis in the series of 2-substituted selenophenes have been carried out. A new basis set optimized for the calculation of (77)Se-(1)H spin-spin coupling constants has been introduced by extending the aug-cc-pVTZ-J basis for selenium. Most of the spin-spin coupling constants under study, especially vicinal (77)Se-(1)H couplings, demonstrated a remarkable stereochemical behavior with respect to the internal rotation of the substituent in the 2-position of the selenophene ring, which is of major importance in the stereochemical studies of the related organoselenium compounds.