Attacking boron nucleophiles: NMR properties of five-membered diazaborole rings

This paper reports computed NMR spectral data for the diazaborole anion (C2H4B1N2-) and the corresponding neutral five-membered rings with B-H (C2H5B1N2, diazaborole) and B-Li (C2H4B1Li1N2, Li-diazaborole) bonds, which are the central moieties of newly synthesized nucleophilic organoboryl five-membered rings, recently reported by Segawa et al. (Science 2006, 314, 113). Both spin-spin coupling constants and chemical shifts were obtained using high-level ab initio calculations. These data are a necessary complement to the very scarce experimental information available.     

Computed coupling constants in X(CH3)nH(4-n) moieties where X = 13C and 15N+, and n = 0-4: comparisons with experimental data

Seventy-three unique spin-spin coupling constants have been analyzed for the ten species in the two series X(CH3)nH(4-n), where the central atom X is 13C or 15N+. Thirty-seven experimental values have been obtained from the literature, and several new coupling constants have been measured for the methyl-substituted ammonium ions. Both DFT with the B3LYP functional and ab initio EOM-CCSD calculations have been carried out on these same systems. Coupling constants computed by these two methods are in agreement with experimental values. Some problems related to coupling constants for the cationic ammonium systems have been resolved when these were recomputed at EOM-CCSD for complexes in which NH4+ is hydrogen-bonded to H2O molecules.     

Can changes in one-bond spin-spin coupling constants in acids be related to gas-phase proton affinities of bases?

Ab initio EOM-CCSD calculations have been carried out in search of acids with one-bond spin-spin coupling constants which could serve as probes of molecular basicity upon complex formation. Only 1J(F-H) and 1J(B-Li) in complexes F-H...Y and H2B-Li...Y, respectively, are able to differentiate among the protonation energies of bases containing Y. Changes in 1J(F-H) upon complex formation are related to the protonation energies of bases, but only when these energies are about 200 kcal/mol or greater. On the other hand, changes in 1J(B-Li) upon complexation are related to base protonation energies in the range between 100 and about 190 kcal/mol. Thus, these two one-bond coupling constants 1J(F-H) and 1J(B-Li) are complementary probes. For the first time, the computed NMR property of a gas-phase one-bond spin-spin coupling constant has been related to the experimental gas-phase thermodynamic property of proton affinity.     

DFT calculations of 31P spin-spin coupling constants and chemical shift in dioxaphosphorinanes

One-bond heteronuclear spin-spin coupling constants (1)J(PX) (X=H, O, S, Se, C and N) between the phosphorus atom and axial and equatorial substituents in dioxaphosphorinanes are computed using density functional theory (DFT). The experimental values of these coupling constants for a variety of substituents can be applied to identify different diastereoisomers. The DFT calculations confirm the systematic trend observed in experiment, and indicate that the computed (1)J(PX) coupling constants are related to the length of the axial and equatorial bonds. A similar relation between the phosphorus chemical shift and the R(PX) bond length appears to be valid, with the exception of selenium substituents.     

A study on the intermolecular hydrogen bonds of alpha-glycylglycine in its actual crystalline phase using ab initio calculated 14N and 2H nuclear quadrupole coupling constants

alpha-Glycylglycine in its actual crystalline phase is studied by ab initio calculated nuclear quadrupole coupling constants. These physical quantities are computed for 2H and 14N in the hydrogen bonds. The type of hydrogen bond is the N-H...O type. The computations are performed with the RHF and B3LYP methods and 6-31++G** and 6-311++G** basis sets using the Gaussian 98 program. Values of the calculated nuclear quadrupole coupling constants are shown in Tables 1-3. The aim of this work is the study of 2H and 14N quadrupole coupling constants which contribute in the CON2H...O=CN2H type of hydrogen bond. The computed nuclear quadrupole coupling constants of 2H nuclei meet the related experimental values. In addition, the computed chi value of 14N belonging to the -CO-14NH- group agrees well with values obtained experimentally. However, there are some discrepancies between calculated 14N chi values of the N+H3 residue and experiments. Also, the values of these physical parameters are calculated for >C2H2 of alpha-glycylglycine in its crystalline phase. Calculations for these parameters are carried out in a single molecule using X-ray diffraction coordinates, too.     

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.     

Dielectric permittivity and temperature effects on spin-spin couplings studied on acetonitrile

Dielectric permittivity (epsilon) and temperature effects on indirect spin-spin coupling constants were studied using acetonitrile as a probe molecule. Experiments were accompanied by hybrid DFT (density functional theory) studies, where the solvent was modeled using the polarization continuum model. Owing to its numerous types of J-couplings, acetonitrile is a very convenient molecule against which various basis sets can be tested or the best basis set can be selected for a given study. The results show reasonable agreement between calculated and experimental values. According to our data, scalar spin-spin coupling constants undergo substantial shifts at lower values of the dielectric constant. Thus J-coupling values are not transferable between measurements made at differing epsilon-conditions, and the assumption of the epsilon-independence of the J-coupling can lead to crucial mistakes in experiments using low-epsilon media. Dielectric permittivity also causes small geometric fluctuations within the molecule, which themselves can affect J-coupling values. Examinations of the results computed with frozen and relaxed geometries show that geometry mediation mostly affects the spin-dipole term of the J-coupling; hence, for accurate evaluation of the latter, frozen geometries are not acceptable. Another interesting fact revealed is the connection between the solvent dielectric properties and the temperature-dependence slopes of J-couplings in corresponding media.     

Dynamics of dilute water in carbon tetrachloride

A dilute solution of water in a hydrophobic solvent, such as carbon tetrachloride (CCl4), presents an opportunity to study the rotational properties of water without the complicating effects of hydrogen bonds. We report here the results of theoretical, experimental, and semiempirical studies of a 0.03 mole percent solution of water in CCl4. It is shown that for this solution there are negligible water-water interactions or water-CCl4 interactions; theoretical and experimental values for proton NMR chemical shifts (deltaH) are used to confirm the minimal interactions between water and the CCl4. Calculated ab initio values and semiempirical values for oxygen-17 and deuterium quadrupole coupling constants (chi) of water/CCl4 clusters are reported. Experimental values for the 17O, 2H, and 1H NMR spin-lattice relaxation times, T1, of 0.03 mole percent water in dilute CCl4 solution at 291 K are 94+/-3 ms, 7.0+/-0.2 s, and 12.6+/-0.4 s, respectively. These T1 values for bulk water are also referenced. "Experimental" values for the quadrupole coupling constants and relaxation times are used to obtain accurate, experimental values for the rotational correlation times for two orthogonal vectors in the water molecule. The average correlation time, tauc, for the position vector of 17O (orthogonal to the plane of the molecule) in monomer water, H2(17)O, is 91 fs. The average value for the deuterium correlation time for the deuterium vector in 2H2O is 104 fs; this vector is along the OD bond. These values indicate that the motion of monomer water in CCl4 is anisotropic. At 291 K, the oxygen rotational correlation time in bulk 2H2(17)O is 2.4 ps, the deuterium rotational correlation time in the same molecule is 3.25 ps. (Ropp, J.; Lawrence, C.; Farrar, T. C.; Skinner, J. L. J. Am. Chem. Soc. 2001, 123, 8047.) These values are a factor of about 20 longer than the tauc value for dilute monomer water in CCl4.     

Structure and dynamics of water confined in dimethyl sulfoxide.

We study the structure and dynamics of hydrogen-bonded complexes of H2O/D2O and dimethyl sulfoxide (DMSO) by infrared spectroscopy, NMR spectroscopy and ab initio calculations. We find that single water molecules occur in two configurations. For one half of the water monomers both OH/OD groups form strong hydrogen bonds to DMSO molecules, whereas for the other half only one of the two OH/OD groups is hydrogen-bonded to a solvent molecule. The H-bond strength between water and DMSO is in the order of that in bulk water. NMR deuteron relaxation rates and calculated deuteron quadrupole coupling constants yield rotational correlation times of water. The molecular reorientation of water monomers in DMSO is two-and-a-half times slower than in bulk water. This result can be explained by local structure behavior.     

Ozonic acid and its ionic salts: ab initio probing of the O(4)(2)(-) dianion

The pyramidal O(4)(2)(-) dianion is valence isoelectronic to the well-known ClO(3)(-) and SO(3)(2)(-) anions, and yet it has not been observed. The synthesis of any molecule containing such a dianion would represent a major breakthrough in making molecules containing more than three covalently bound oxygen atoms. We found that the parent H(2)O(4) ozonic acid is unstable in the form of the valence isoelectronic sulfurous acid H(2)SO(3). Our quantum chemical probing of the Li(2)O(4) ionic salt molecule is inconclusive. However, we found that the specially designed FLi(3)O(4) gas phase molecule is a true metastable species and could be considered as the first molecule containing the O(4)(2)(-) dianion. Our theoretical prediction of the first compound containing the tetraatomic covalently bound O(4)(2)(-) dianion opens the possibility to even more oxygen rich compounds, which will have a great potential as high density oxygen storage.     

Modeling solvent effects on electron-spin-resonance hyperfine couplings by frozen-density embedding

In this study, we investigate the performance of the frozen-density embedding scheme within density-functional theory [J. Phys. Chem. 97, 8050 (1993)] to model the solvent effects on the electron-spin-resonance hyperfine coupling constants (hfcc's) of the H2NO molecule. The hfcc's for this molecule depend critically on the out-of-plane bending angle of the NO bond from the molecular plane. Therefore, solvent effects can have an influence on both the electronic structure for a given configuration of solute and solvent molecules and on the probability for different solute (plus solvent) structures compared to the gas phase. For an accurate modeling of dynamic effects in solution, we employ the Car-Parrinello molecular-dynamics (CPMD) approach. A first-principles-based Monte Carlo scheme is used for the gas-phase simulation, in order to avoid problems in the thermal equilibration for this small molecule. Calculations of small H2NO-water clusters show that microsolvation effects of water molecules due to hydrogen bonding can be reproduced by frozen-density embedding calculations. Even simple sum-of-molecular-densities approaches for the frozen density lead to good results. This allows us to include also bulk solvent effects by performing frozen-density calculations with many explicit water molecules for snapshots from the CPMD simulation. The electronic effect of the solvent at a given structure is reproduced by the frozen-density embedding. Dynamic structural effects in solution are found to be similar to the gas phase. But the small differences in the average structures still induce significant changes in the computed shifts due to the strong dependence of the hyperfine coupling constants on the out-of-plane bending angle.     

DFT study of the NMR properties of xenon in covalent compounds and van der waals complexes-implications for the use of 129Xe as a molecular probe

The NMR properties (chemical shift and spin-spin coupling constants) of (129)Xe in covalent compounds and weakly bound complexes have been investigated by DFT methods including relativistic effects. For covalent species, a good agreement between experimental and calculated results is achieved without scalar relativistic effects, but their inclusion (with a triple-zeta, double-polarization basis set) leads to some improvement in the quality of the correlation. The spin-orbit coupling term has a significant effect on the shielding constant, but makes a small contribution to the chemical shift. Coupling constants contain substantial contributions from the Fermi contact and paramagnetic spin-orbit terms; unlike light nuclei the spin-dipole term is also large, whereas the diamagnetic spin-orbit term is negligible. For van der Waals dimers, the dependence of the xenon chemical shift and anisotropy is calculated as a function of the distance. Small (<1 Hz) but non-negligible through-space coupling constants between (129)Xe and (13)C or (1)H are predicted. Much larger couplings, of the order of few Hz, are calculated between xenon and (17)O in a model silicate residue.     

Structural and magnetic properties of vanadyl dichloride solvates: from molecular units to extended hydrogen-bonded solids

The preparation, structural characterization and magnetic properties of three solvent adducts of VOCl(2), trans-VOCl(2)(THF)(2)(H(2)O) (1), trans-VOCl(2)(H(2)O)(2).2Et(2)O (2) and cis-VOCl(2)(MeOH)(3) (3) are described. In these solids, hydrogen bonding among the inorganic complexes is the critical determinant of the formation of extended magnetic networks. Compound forms one-dimensional double chains where alternating monomers from the two branches of the chain are hydrogen bonded via the V-Cl ... H-O-V network (with an axial water molecule and equatorial chloride ions). Magnetic studies indicate no interaction among the vanadyl centers. The paramagnetism of 1 is consistent with the extension of the network from the hydrogen donor site of the axial water, which is orthogonal to the d(xy) magnetic orbital. Compound 2 forms one-dimensional chains with water molecules of adjacent monomers held together by hydrogen bonds to ether molecules (V-O-H ... O(ether) ... H -O-V). The chain network radiates only through the equatorial plane of the complex where the water molecules are located. The presence of the intervening solvent molecule between hydrogen bonds of the primary coordination sphere magnetically insulates metal centers and compound is also a simple paramagnet. Removal of the solvent turns on the magnetic interaction and neighboring spin centers couple antiferromagnetically. Compound 3 forms a layered structure via V-Cl ... H-O-V hydrogen bonding, where all the hydrogen donor sites participate in the formation of the network. The vanadyl spin centers, at distances of 5.5 and 6.5 A from each other, couple antiferromagnetically (J/k=-0.7 K). Thus, magnetic coupling among metal centers is achieved when the hydrogen bond network directly radiates from the coordination plane containing the magnetic orbital. These results further support the utility of hydrogen bond as a viable design element in the construction of low dimensional, magnetic solids.     

Intermolecular interactions in uracil–nitrous acid complexes: structures,binding energy,topological properties,and nuclear magnetic resonance study

Molecular interactions between uracil and nitrous acid (U–NA) [C₄N₂O₂H₄? NO₂H] have been studied using B3LYP, B3PW91, and MP2 methods with different basis sets. The optimized geometries, harmonic vibrational frequencies, charge transfer, topological properties of electron density, nucleus‐independent chemical shift (NICS), and nuclear magnetic resonance one‐ and two‐bonds spin–spin coupling constants were calculated for U–NA complexes. In interaction between U and NA, eight cyclic complexes were obtained with two intermolecular hydrogen bonds N(C)HU…N(O) and OH_NA…O_U. In these complexes, uracil (U) simultaneously acts as proton acceptor and proton donor. The most stable complexes labeled, UNA1 and UNA2, are formed via NH bond of U with highest acidity and CO group of U with lowest proton affinity. There is a relationship between hydrogen bond distances and the corresponding frequency shifts. The solvent effect on complexes stability was examined using B3LYP method with the aug‐cc‐pVDZ basis set by applying the polarizable continuum model (PCM). The binding energies in the gas phase have also been compared with solvation energies computed using the PCM. Natural bond orbital analysis shows that in all complexes, the charge transfer takes place from U to NA. The results predict that the Lone Pair (LP)(O)_U → σ*(O? H) and LP(N(O)_NA → σ*(N(C)? H)_U donor–acceptor interactions are most important interactions in these complexes. Atom in molecule analysis confirms that hydrogen bond contacts are electrostatic in nature and covalent nature of proton donor groups decreases upon complexation. The relationship between spin–spin coupling constant (^1hJ_H…_Y and ^2hJ_H…_Y) with interaction energy and electronic density at corresponding hydrogen bond critical points and H‐bonds distances are investigated. NICS used for indicating of aromaticity of U ring upon complexation. © 2013 Wiley Periodicals, Inc.     

(3h)J((15)N-(31)P) spin-spin coupling constants across N[bond]H....O[bond]P hydrogen bonds

Equation-of-motion coupled cluster singles and doubles (EOM-CCSD) calculations have been performed to evaluate three-bond (15)N-(31)P coupling constants ((3h)J(N[bond]P)) across N[bond]H....O[bond]P hydrogen bonds in model cationic and anionic complexes including NH(4)(+):OPH, NH(4)(+):OPH(3), NH(3):(-)O(2)PH(2), NFH(2):(-)O(2)PH(2), and NF(2)H:(-)O(2)PH(2). Three-bond coupling constants can be appreciable when the phosphorus is P(V), but are negligible with P(III). (3h)J(N[bond]P) values in complexes with cyclic or open structures are less than 1 Hz, a consequence of the nonlinear arrangement of N, H, O, and P atoms. For complexes with these structures, (3h)J(N[bond]P) may not be experimentally measurable. In contrast, complexes in which the N, H, O, and P atoms are collinear or nearly collinear have larger values of (3h)J(N[bond]P), even though the N[bond]P distances are longer than N[bond]P distances in cyclic and open structures. In linear complexes, (3h)J(N[bond]P) is dominated by the Fermi-contact term, which is distance dependent. Therefore, N[bond]P (and hydrogen-bonding N[bond]O) distances in these complexes can be determined from experimentally measured (15)N-(31)P coupling constants.     

Hydrogen bond geometries from electron paramagnetic resonance and electron-nuclear double resonance parameters: density functional study of quinone radical anion-solvent interactions

Density functional theory was used to study the impact of hydrogen bonding on the p-benzosemiquinone radical anion BQ(*-) in coordination with water or alcohol molecules. After complete geometry optimizations, (1)H, (13)C, and (17)O hyperfine as well as (2)H nuclear quadrupole coupling constants and the g-tensor were computed. The suitability of different model systems with one, two, four, and 20 water molecules was tested; best agreement between theory and experiment could be obtained for the largest model system. Q-band pulse (2)H electron-nuclear double resonance (ENDOR) experiments were performed on BQ(*-) in D(2)O. They compare very well with the spectra simulated by use of the theoretical values from density functional theory. For BQ(*-) in coordination with four water or alcohol molecules, rather similar hydrogen-bond lengths between 1.75 and 1.78 A were calculated. Thus, the computed electron paramagnetic resonance (EPR) parameters are hardly distinguishable for the different solvents, in agreement with experimental findings. Furthermore, the distance dependence of the EPR parameters on the hydrogen-bond length was studied. The nuclear quadrupole and the dipolar hyperfine coupling constants of the bridging hydrogens show the expected dependencies on the H-bond length R(O.H). A correlation was obtained for the g-tensor. It is shown that the point-dipole model is suitable for the estimation of hydrogen-bond lengths from anisotropic hyperfine coupling constants of the bridging (1)H nuclei for H-bond lengths larger than approximately 1.7 A. Furthermore, the estimation of H-bond lengths from (2)H nuclear quadrupole coupling constants of bridging deuterium nuclei by empirical relations is discussed.     

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.     

One-bond 1J(15N,H) coupling constants at sp2-hybridized nitrogen of Schiff bases,enaminones and similar compounds: A theoretical study

¹J(¹⁵N,H) coupling constants for enaminones and NH-forms of intramolecularly hydrogen-bonded Schiff bases as model compounds for sp²-hybridized nitrogen atoms are evaluated using density functional theory (DFT) to find the optimal functionals and basis sets. Ammonia is used as a test molecule and its one-bond coupling constant is compared with experiment. A methylamine Schiff base of a truncated molecule of gossypol is used for checking the performance of selected B3LYP, O3LYP, PBE, BHandH, and APFD density functionals and standard, modified, and dedicated basis sets for coupling constants. Both in vacuum and in chloroform, modeled by the simple continuum model of solvent, the modified basis sets predict significantly better the ¹J(¹⁵N,H) value in ammonia and in the methylamine Schiff base of a truncated molecule of gossypol than the standard basis sets. This procure is then used on a broad set of intramolecularly hydrogen-bonded molecules, and a good correlation between calculated and experimental one-bond NH coupling constants is obtained. The ¹J(¹⁵N,H) couplings are slightly overestimated. The calculated data show for hydrogen-bonded NH interatomic distances that the calculated values depend on the NH bond lengths. The shorter the bond lengths, the larger the ¹J(¹⁵N,H). A useful correlation between ¹J(¹⁵N,H) and NH bond length is derived that enables realistic predictions of one-bond NH coupling constants. The calculations reproduce experimentally observed trends for the studied molecules.     

Spectroscopic investigation of the structures of dialkyl tartrates and their cyclodextrin complexes

Structures of three dialkyl tartrates, namely, dimethyl tartrate, diethyl tartrate, and diisopropyl tartrate, in CCl4, dimethyl sulfoxide (DMSO)/DMSO-d6, and H2O/D2O solvents have been investigated using vibrational absorption (VA), vibrational circular dichroism (VCD), and optical rotatory dispersion (ORD). VA, VCD, and ORD spectra are found to be dependent on the solvent used. Density functional theory (DFT) calculations are used to interpret the experimental data in CCl4 and DMSO. The trans-COOR conformer with hydrogen bonding between the OH group and the C=O group attached to the same chiral carbon is dominant for dialkyl tartrates both in vacuum and in CCl4. The experimental VA, VCD, and ORD data of dialkyl-D-tartrates in CCl4 correlated well with those predicted for dimethyl-(S,S)-tartrate molecule as both isolated and solvated in CCl4. In DMSO solvent, dialkyl tartrate molecules favor formation of intermolecular hydrogen bonding with DMSO molecules. Clusters of dimethyl-(S,S)-tartrate, with one molecule of dimethyl-(S,S)-tartrate hydrogen bonded to two DMSO molecules, are used for the DFT calculations. A trans-COOR cluster and a trans-H cluster are needed to obtain a reasonable agreement between the predicted and experimental data of dimethyl tartrate in DMSO solvent. VA, VCD, and optical rotations are also measured for dialkyl tartrate-cyclodextrin complexes. It is noted that these properties are barely affected by complexation of dialkyl tartrates with cyclodextrins, indicating weak interaction between tartrates and cyclodextrin. Binding constants of alpha-CD and beta-CD with diethyl L-tartrate in both H2O and DMSO have been determined using isothermal titration calorimetry technique. The smaller binding constants (less than 100) confirmed the weak interaction between tartrates and cyclodextrin in the solution state.     

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.