Resolving an apparent discrepancy between theory and experiment: spin-spin coupling constants for FCCF

Ab initio equation of motion coupled cluster singles and doubles (EOM-CCSD) and second-order polarization propagator approximation (SOPPA) calculations have been performed to evaluate spin-spin coupling constants for FCCF (difluoroethyne). The computed EOM-CCSD value of (3)J(F-F) obtained at the experimental geometry of this molecule supports the previously reported experimental value of 2.1 Hz, thereby resolving an apparent discrepancy between theory and experiment. This coupling constant exhibits a strong dependence on the C-C and C-F distances, and its small positive value results from a sensitive balance of paramagnetic spin-orbit (PSO) and spin-dipole (SD) terms. The three other unique FCCF coupling constants (1)J(C-C), (1)J(C-F), and (2)J(C-F) have also been reported and compared with experimental data. While (1)J(C-F) is in agreement with experiment, the computed value of (2)J(C-F) is larger than our estimate of the experimental coupling constant.     

Do traditional, chlorine-shared, and ion-pair halogen bonds exist? An ab initio investigation of FCl:CNX complexes

Ab initio MP2/aug'-cc-pVTZ calculations have been carried out to determine the structures, binding energies, and bonding of complexes FCl:CNX, with X = CN, NC, NO(2), F, CF(3), Cl, Br, H, CCF, CCH, CH(3), SiH(3), Li, and Na. Equation-of-motion coupled cluster calculations have also been carried out to determine the coupling constants (1)J(F-Cl), (1X)J(Cl-C), and (2X)J(F-C) across these halogen bonds. As the strength of the base is systematically increased, the nature of the halogen bond changes from traditional, to chlorine-shared, to ion-pair. The type of halogen bond present in a complex can be readily determined from its structure, binding energy, AIM bonding analyses, and spin-spin coupling constants. Coupling constants across halogen bonds are compared with corresponding coupling constants across traditional, proton-shared, and ion-pair hydrogen bonds.     

Probing the proton-transfer coordinate of complexes with F--H...P hydrogen bonds using one- and two-bond spin-spin coupling constants

Scalar coupling constants have been computed using the EOM-CCSD method for equilibrium structures of complexes stabilized by F--H...P hydrogen bonds, as well as structures along the proton-transfer coordinates of these complexes. Variations in the signs and absolute values of (1)J(F--H), (1h)J(H--P) and (2h)J(F--P) have been analyzed and interpreted in terms of changing hydrogen bond type. Of the three phosphorus bases (phosphine, trimethylphosphine and phosphinine) investigated in this study, trimethylphosphine forms the strongest complex with FH, and has the largest two-bond F--P coupling constant. Among the relatively simple phosphorus bases, it would appear to be a leading candidate for experimental NMR study. Similarities and differences are noted between the corresponding coupling constants (J) and the reduced coupling constants (K) across F--H...P and F--H...N hydrogen bonds.     

About the calculation of exchange coupling constants using density-functional theory: the role of the self-interaction error

The effect of the correction of the self-interaction error on the calculation of exchange coupling constants with methods based on density-functional theory has been tested in simple model systems. The inclusion of the self-interaction correction cancels the nondynamical correlation energy contributions simulated by the commonly used functionals. Hence, such correction should be important in the accurate determination of exchange coupling constants. We have also tested several recent functionals to calculate exchange coupling constants in transition-metal complexes, such as meta-GGA functionals or new formulations of hybrid functionals. The influence of the basis set and of the use of pseudopotentials on the calculated J values has also been evaluated for a Fe(III) dinuclear complex in which the paramagnetic centers bear several unpaired electrons.     

Substitution and protonation effects on spin-spin coupling constants in prototypical aromatic rings: C6H6, C5H5N and C5H5P

Ab initio equation-of-motion coupled cluster calculations have been carried out to evaluate one-, two-, and three-bond 13C-13C, 15N-13C, 31P-13C coupling constants in benzene, pyridine, pyridinium, phosphinine, and phosphininium. The introduction of N or P heteroatoms into the aromatic ring not only changes the magnitudes of the corresponding X-C coupling constants (J, for X = C, N, or P) but also the signs and magnitudes of corresponding reduced coupling constants (K). Protonation of the heteroatoms also produces dramatic changes in coupling constants and, by removing the lone pair of electrons from the sigma-electron framework, leads to the same signs for corresponding reduced coupling constants for benzene, pyridinium, and phosphininium. C-C coupling constants are rather insensitive to the presence of the heteroatoms and protonation. All terms that contribute to the total coupling constant (except for the diamagnetic spin-orbit (DSO) term) must be computed if good agreement with experimental data is to be obtained.     

A comparison of two-component and four-component approaches for calculations of spin-spin coupling constants and NMR shielding constants of transition metal cyanides

Relativistic density functional theory (DFT) calculations of nuclear spin-spin coupling constants and shielding constants have been performed for selected transition metal (11th and 12th group of periodic table) and thallium cyanides. The calculations have been carried out using zeroth-order regular approximation (ZORA) Hamiltonian and four-component Dirac-Kohn-Sham (DKS) theory with different nonrelativistic exchange-correlation functionals. Two recent approaches for representing the magnetic balance (MB) between the large and small components of four-component spinors, namely, mDKS-RMB and sMB, have been employed for shielding tensor calculations and their results have been compared. Relativistic effects have also been analysed in terms of scalar and spin-orbit contributions at the two-component level of theory, including discussion of heavy-atom-on-light-atom effects for (1)J(CN), σ(C), and σ(N). The results for molecules containing metals from 4th row of periodic table show that relativistic effects for them are small (especially for spin-spin coupling constants). The biggest effects are observed for the 6th row where nonrelativistic theory reproduces only about 50%-70% of the two-component ZORA results for (1)J(MeC) and about 75% for heavy metal shielding constants. It is important to employ a full Dirac picture for calculations of heavy metal shielding constants, since ZORA reproduces only 75%-90% of the DKS results. Smaller discrepancies between ZORA-DFT and DKS are observed for nuclear spin-spin coupling constants. No significant differences are observed between the results obtained using mDKS-RMB and sMB approaches for magnetic balance in four-component calculations of the shielding constants.     

A theoretical study of the NMR spin-spin coupling constants of the complexes [(NC)(5)Pt-Tl(CN)(n)](n-) (n = 0-3) and [(NC)(5)Pt-Tl-Pt(CN)(5)](3-): a lesson on environmental effects

The molecular geometries and the nuclear spin-spin coupling constants of the complexes [(NC)(5)Pt-Tl(CN)(n)](n-), n = 0-3, and the related system [(NC)(5)Pt-Tl-Pt(CN)(5)](3-) are studied. These complexes have received considerable interest since the first characterization of the n = 1 system by Glaser and co-workers in 1995 [J. Am. Chem. Soc. 1995, 117, 7550-7551]. For instance, these systems exhibit outstanding NMR properties, such as extremely large Pt-Tl spin-spin coupling constants. For the present work, all nuclear spin-spin coupling constants J(Pt-Tl), J(Pt-C), and J(Tl-C) have been computed by means of a two-component relativistic density functional approach. It is demonstrated by the application of increasingly accurate computational models that both the huge J(Pt-Tl) for the complex (NC)(5)Pt-Tl and the whole experimental trend among the series are entirely due to solvent effects. An approximate inclusion of the bulk solvent effects by means of a continuum model, in addition to the direct coordination, proves to be crucial. Similarly drastic effects are reported for the coupling constants between the heavy atoms and the carbon nuclei. A computational model employing the statistical average of orbital-dependent model potentials (SAOP) in addition to the solvent effects allows to accurately reproduce the experimental coupling constants within reasonable limits.     

Broken-symmetry density functional theory investigation on bis-nitronyl nitroxide diradicals: influence of length and aromaticity of couplers

A series of nitronyl nitroxide (NN) diradicals with linear conjugated couplers and another series with aromatic couplers have been investigated by the broken-symmetry (BS) DFT approach. The overlap integral between the magnetically active orbitals in the BS state has been explicitly computed and used for the evaluation of the magnetic exchange coupling constant (J). The calculated J values are in very good agreement with the observed values in the literature. The magnitude of J depends on the length of the coupler as well as the conformation of the radical units. The aromaticity of the spacer decreases the strength of the exchange coupling constant. The SOMO-SOMO energy splitting analysis, where SOMO stands for the singly occupied molecular orbital, and the calculation of electron paramagnetic resonance (EPR) parameters have also been carried out. The computed hyperfine coupling constants support the intramolecular magnetic interactions. The nature of magnetic exchange coupling constant can also be predicted from the shape of the SOMOs as well as the spin alternation rule in the unrestricted Hartree-Fock (UHF) treatment. It is found that pi-conjugation along with the spin-polarization plays the major role in controlling the magnitude and sign of the coupling constant.     

A DFT study on the magnetostructural property of ferromagnetic heteroverdazyl diradicals with phenylene coupler

We have theoretically designed five different m-phenylene coupled high-spin bis-heteroverdazyl diradicals and their analogous p-phenylene coupled low-spin positional isomers. The geometry-based aromaticity index, harmonic oscillator model of aromaticity (HOMA) values for both the couplers (local HOMA), and the whole diradicals (global HOMA) have been calculated for all the diradicals. We also qualitatively relate these HOMA values with the intramolecular magnetic exchange coupling constants (J), calculated using a broken symmetry approach within unrestricted density functional theory framework. Structural aromaticity index HOMA of linkage specific benzene rings in our designed diradical systems shows that the aromatic character depends on the planarity of the molecule and it controls the sign and magnitude of J. The predicted J values are explained on the basis of spin polarization maps, average dihedral angles, and magnetic orbitals. The effect of the spin leakage phenomenon on magnetic exchange coupling constant and that on HOMA values of certain phosphaverdazyl systems has been explicitly discussed. In addition, a similar comparison is made between the calculated exchange coupling constants and corresponding HOMA values. The main novelty of this work stands on the consideration of the aromatic behavior by means of the geometrical index HOMA. We also estimate another aromaticity index, nucleus independent chemical shift (NICS) values for the phenylene coupler in each diradical to measure aromaticity and compare its value with that of HOMA. The ground state stabilities of these diradicals have also been compared.     

Theoretical prediction and assignment of vicinal 1H-1H coupling constants of diastereomeric 3-alkoxy-6,7-epoxy-2-oxabicyclo[3.3.0]octanes

Spin-spin coupling constants between nuclei in NMR spectroscopy reflect their spatial arrangement. A number of calculation methods, applying different levels of theory, have been developed to support the stereochemical assignment of novel compounds. Nevertheless, revisions of the assignment of structures in the literature are not rare. In the present work, the reliability of the calculation methods amenable for a theoretical prediction of spin-spin coupling constants of vicinal protons to support correct stereochemical assignment of substitution at five-membered rings of 3-alkoxy-6,7-epoxy-2-oxabicyclo[3.3.0]octanes was studied. Experimental (3)J(H,H) coupling constants were compared with the coupling constants calculated for all possible diastereomers. The fully quantum chemical approach provided theoretical (3)J(H,H) coupling constants with an absolute deviation of no more than 1.1 Hz for 91% of the experimentally studied coupled spins, whereas the methods without quantum chemical geometry optimization resulted in completely unreliable predictions. Consequently, for a reliable stereochemical assignment of small and medium size molecules, the protocol for calculating the coupling constants based on the results of the quantum chemical geometry optimization is recommended.     

Ab initio EOM-CCSD spin-spin coupling constants for hydrogen-bonded formamide complexes: bridging complexes with NH3, (NH3)2, H2O, (H2O)2, FH, and (FH)2

EOM-CCSD spin-spin coupling constants across hydrogen bonds have been computed for complexes in which NH3, H2O, and FH molecules and their hydrogen-bonded dimers form bridging complexes in the amide region of formamide. The formamide one-bond N-H coupling constant [(1)J(N-H)] across N-H...X hydrogen bonds increases in absolute value upon complexation. The signs of the one-bond coupling constants (1h)J(H-X) indicate that these complexes are stabilized by traditional hydrogen bonds. The two-bond coupling constants for hydrogen bonds with N-H as the donor [(2h)J(N-X)] and the carbonyl oxygen as the acceptor [(2h)J(X-O)] increase in absolute value in the formamide/dimer relative to the corresponding formamide/monomer complex as the hydrogen bonds acquire increased proton-shared character. The largest changes in coupling constants are found for complexes of formamide with FH and (FH)2, suggesting that bridging FH monomers and dimers in particular could be useful NMR spectroscopic probes of amide hydrogen bonding.     

Complexes with N-H(+)-P hydrogen bonds: structures, binding energies, and spin-spin coupling constants

Ab-initio MP2/aug'-cc-pVTZ calculations have been performed to determine the structures and binding energies of proton-bound complexes stabilized by N-H+-P hydrogen bonds and to investigate the nature of the proton-transfer coordinate in these systems. Double minima are found only when the difference between the protonation energies of the N and P bases is less than about 4 kcal/mol. The isomer in which the protonated nitrogen base is the donor lies lower on the potential surface and also has a greater binding energy relative to the corresponding isolated monomers. Equation-of-motion coupled cluster singles and doubles (EOM-CCSD) calculations have been employed to obtain one- and two-bond spin-spin coupling constants across these hydrogen bonds. Two-bond coupling constants (2h)J(N-P) correlate with N-P distances, irrespective of whether the donor ion is N-H+ or P-H+. One-bond coupling constants (1)J(N-H) and (1h)J(H-P) for complexes stabilized by N-H+...P hydrogen bonds correlate with corresponding distances, but similar correlations are not found for (1)J(P-H) and (1h)J(H-N) for complexes with P-H+...N hydrogen bonds. Negative values of (1h)K(H-N) and (1h)K(H-P) indicate that the hydrogen bonds in these complexes are traditional. Comparisons are made with complexes stabilized by N-H+-N and P-H+-P hydrogen bonds.     

One-bond spin-spin coupling constants of X-1H proton donors in complexes with X-H-Y hydrogen bonds, for X = 13C, 15N, 17O, and 19F: predictions, comparisons, and relationships among 1J X-H, 1K X-H, and X-H distances

Ab initio calculations at the equation-of-motion coupled cluster (EOM-CCSD) level of theory have been carried out to investigate one-bond (13)C-(1)H, (15)N-(1)H, (17)O-(1)H, and (19)F-(1)H coupling constants in a systematic study of monomers and hydrogen-bonded complexes. Computed coupling constants ((1)J(X-H)) for monomers are in good agreement with available experimental data. All reduced Fermi-contact terms and reduced coupling constants ((1)K(X-H)) for monomers and complexes are positive. Plots of (1)K(X-H) versus the X-H distance for the 16 monomers and the 64 complexes in which these monomers are proton donors exhibit significant scatter. However, a linear relationship has been demonstrated for the first time between coupling constants and X-H distances for different X atoms by plotting the ratios of the coupling constants for complexes and corresponding monomers versus the ratios of distances for complexes and corresponding monomers times the square of the Pauling electronegativity. Since the ratio removes the dependence of coupling constants on the magnetogyric ratios of X, this relationship holds for both (1)K(X-H) and (1)J(X-H). The decrease in reduced coupling constants ((1)K(X-H)) as the X-H distance increases is due primarily to the increased proton-shared character of the hydrogen bond.     

Characterization of monofluorinated polycyclic aromatic compounds by 1H, 13C and 19F NMR spectroscopy

Monofluorinated polycyclic aromatic hydrocarbons (F-PAHs) have attracted much attention in analytical, environmental, toxicological and mechanistic studies because of their physico-chemical properties, which are closely similar to those of the parent PAHs. Because of this, full NMR characterization has become of interest. Complete 1H, 13C and 19F NMR chemical shifts, and also 1J(H,C), (n)J(C,F), (n)J(H,F) and (n)J(H,H) coupling constants, have been assigned for the F-PAHs 1-fluoronaphthalene, 2-fluorofluorene, 5-fluoroacenaphthylene, 2-fluorophenanthrene, 3-fluorophenanthrene, 3-fluorofluoranthene, 1-fluoropyrene, 1-fluorochrysene, 2-fluorochrysene, 3-fluorochrysene and 9-fluorobenzo[k]fluoranthene. To allow comparison with the corresponding parent PAHs, the 1H and 13C chemical shifts of acenaphthylene, phenanthrene, fluoranthene, pyrene and benzo[k]fluoranthene were determined. Chemical shift increments and the effects on the coupling constants from the fluorine substitution are discussed.     

A theoretical study of the large Hg-Hg spin-spin coupling constants in Hg(2)(2+), Hg(3)(2+), and Hg(2)(2+)-crown ether complexes

Nuclear spin-spin coupling constants (1)J(Hg-Hg) in the systems Hg(2)(2+) and Hg(3)(2+) represent the largest coupling constants so far observed in NMR experiments. We have performed a computational study on these ions, on Hg(2)(2+) complexes with 18-crown-6 and 15-crown-5, and on Hg(3)(2+) with solvent molecules and counterions. The results obtained with our recently developed program for the density functional computation of heavy nucleus spin-spin coupling constants are in good agreement with experiments. The data reveal that the bare ions Hg(2)(2+) and Hg(3)(2+) would afford much larger coupling constants than those experimentally observed, with an upper limit of approximately 0.9 MHz for Hg(2)(2+). This limit is much larger than that previously estimated by Hückel theory. It is demonstrated that in solution or due to complexation the experimentally determined values are much smaller than the free ion's coupling constants. With the help of intuitive MO arguments, it is illustrated how the environment strongly reduces the coupling constants in Hg(2)(2+) and Hg(3)(2+). The two-bond coupling constant (2)J(Hg-Hg) in Hg(3)(2+) is also examined.     

Influence of the correlation,aggregation, and solvation on ab initio computed Li-C,Li-N,and Li-Li NMR coupling constants

The (1)J and (3)J(C-Li), (1)J(N-Li), and (2)J(Li-Li) NMR coupling constants have been calculated for various homogeneous and heterogeneous aggregates of methyllithium and lithium dimethylamide at the HF and MP2 levels of calculation. Ethereal solvation has also been taken into account either through a continuum model or through the explicit introduction of Me(2)O molecules. The results obtained are in good general agreement with the experimental data available for methyllithium itself or model alkyllithiums and supports the empirical rule proposed by Bauer, Winchester, and Schleyer to evaluate (1)J(C-Li) provided that calculations include solvent and/or aggregation effects.     

Characterizing hydrogen bonding and proton transfer in 2:1 FH:NH3 and FH:collidine complexes through one- and two-bond spin-spin coupling constants across hydrogen bonds

Ab initio equation-of-motion coupled cluster singles and doubles calculations have been carried out on a variety of 2:1 FH:NH(3) complexes (F(b)H(b):F(a)H(a):NH(3)) to investigate the effects of structural changes on one- and two-bond spin-spin coupling constants across F(a)-H(a)-N and F(b)-H(b)-F(a) hydrogen bonds and to provide insight into experimentally measured coupling constants for 2:1 FH:collidine (2:1 FH:2,4,6-trimethylpyridine) complexes. Coupling constants have been computed for 2:1 FH:NH(3) equilibrium structures and proton-transferred perpendicular and open structures at 2:1 FH:NH(3), FH:pyridine, and FH:collidine geometries. (2h)J(Fa)(-)(N), (1)J(Fa)(-)(Ha), and (1h)J(Ha)(-)(N) exhibit expected dependencies on distances, angles, and the nature of the nitrogen base. In contrast, one- and two-bond coupling constants associated with the F(b)-H(b)-F(a) hydrogen bond, particularly (2h)J(F)()b(-)(F)()a, vary significantly depending on the F-F distance, the orientation of the hydrogen-bonded pair, and the nature of the complex (HF dimer versus the anion FHF(-)). The structure of the 2:1 FH:collidine complex proposed on the basis of experimentally measured coupling constants is supported by the computed coupling constants. This study of the structures of open proton-transferred 2:1 FH:NH(3), FH:pyridine, and FH:collidine complexes and the coupling constants computed for 2:1 FH:NH(3) complexes at these geometries provides insight into the role of the solvent in enhancing proton transfer across both N-H(a)-F(a) and F(b)-H(b)-F(a) hydrogen bonds.     

On achieving experimental accuracy from molecular dynamics simulations of flexible molecules: aqueous glycerol

The rotational isomeric states (RIS) of glycerol at infinite dilution have been characterized in the aqueous phase via a 1 micros conventional molecular dynamics (MD) simulation, a 40 ns enhanced sampling replica exchange molecular dynamics (REMD) simulation, and a reevaluation of the experimental NMR data. The MD and REMD simulations employed the GLYCAM06/AMBER force field with explicit treatment of solvation. The shorter time scale of the REMD sampling method gave rise to RIS and theoretical scalar 3J(HH) coupling constants that were comparable to those from the much longer traditional MD simulation. The 3J(HH) coupling constants computed from the MD methods were in excellent agreement with those observed experimentally. Despite the agreement between the computed and the experimental J-values, there were variations between the rotamer populations computed directly from the MD data and those derived from the experimental NMR data. The experimentally derived populations were determined utilizing limiting J-values from an analysis of NMR data from substituted ethane molecules and may not be completely appropriate for application in more complex molecules, such as glycerol. Here, new limiting J-values have been derived via a combined MD and quantum mechanical approach and were used to decompose the experimental 3J(HH) coupling constants into population distributions for the glycerol RIS.