New look at the Badger-Bauer rule: Correlations of spectroscopic IR and NMR parameters with hydrogen bond energy and geometry. FHF complexes

Phosphine selenides: versatile NMR probes for analyzing hydrogen OH⋯Se and halogen I⋯Se bonds

Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool for studying the structure and dynamics of various non-covalent interactions. However, often spectral parameters that are applicable for estimation of parameters of one type of non-covalent interaction will be inapplicable for another. Therefore, researchers are compelled to use spectral parameters that are specifically tailored to the type of non-covalent interaction being studied. This complexity makes it difficult to compare different types of non-covalent interactions with each other and, consequently, to establish a strict unified classification for them. This pioneering study proposes to use phosphine selenides as universal probes for investigating hydrogen and halogen bonding in solution. The study was carried out using the example of triethylphosphine selenide Et3PSe complexes with hydrogen bonds of Se⋯HO type and R3PSe (where R: Me, Et, n-Bu, t-Bu and Ph) with halogen bonds of Se⋯X type (where X: I and Br) in solution. The presence of non-covalent interactions was confirmed experimentally by means of 1H, 31P and 77Se NMR, as well as by quantum chemical calculation methods (optimization: PW6B95-D3/def2-QZVP; NMR: B97-2/pcsSeg-2).

Unveiling the electronic structure peculiarities of phosphine selenides as NMR probes for non-covalent interactions: an experimental and theoretical study

In this work, R3PSe (R = Me, Et, n-Bu, t-Bu and Ph) were studied experimentally using NMR spectroscopy in solution and the solid-state in combinaton with quantum chemical methods. The study shows that the NMR parameters of these phosphine selenides, such as δP, δSe, and 1JPSe, are sensitive to subtle changes in the electronic environment of the P and Se atoms. Consequently, phosphine selenides R3PSe can serve as promising spectral probes for the detection and quantitative investigation of various non-covalent interactions. Additionally, the variations of R in phosphine selenides influence the observed NMR spectral parameters, primarily through effects such as π-backdonation and hyperconjugation, which have been observed experimentally and confirmed theoretically.

Conformational preference of 2-(4-methoxyphenyl)ethanol studied by supersonic jet spectroscopy: Intramolecular OH/π interaction

A π-type hydrogen bonding between the OH group and the π electron is a crucial factor for the conformational preference of the molecular structure with a flexible group. However, the information on the effect of the substituent on the OH/π interaction is insufficient. The laser-induced fluorescence (LIF) excitation, the dispersed fluorescence (DF), the IR-UV hole-burning, and the IR dip spectra of jet-cooled 2-(4-methoxyphenyl)ethanol were measured for the first time. Almost all bands observed in the spectral region of 35 550-36 500 cm-1 in the LIF excitation spectrum were successfully assigned with the DF and the IR-UV hole-burning spectra coupled with the quantum chemical calculation at M06-2x/6-311G and MP2/6-311G levels. Five conformers were found in the LIF excitation spectrum. The most stable conformer was Ggπ, and the second most stable conformer was Ggπ' (the trans rotamer of the methoxy group for Ggπ). Ggπ and Ggπ' had the OH group directed toward the π electron system of the benzene ring. The OH stretching frequency of Ggπ/Ggπ' of MPE in the IR dip spectra was red-shifted against that of Ggπ of phenylethanol, indicating that the introduction of the methoxy group would enhance the intramolecular OH/π interaction. In addition, the torsional vibration between the benzene ring and the side chain (-CH2CH2OH) (mode 63) was observed in the DF spectra of the Ggπ-00 and Ggπ'-00 band excitation, but their intensities were rather different, resulting from the different orientation of the OH group for each conformer toward the π electron system. The methoxy group would increase the negative charge on the benzene ring and would enhance the intramolecular OH/π interaction through the electrostatic interaction.

Experimental confirmation of the Badger-Bauer rule in the protonated methanol clusters: weak hydrogen bond formation as a measure of terminal OH acidity in hydrogen bond networks

We report a linear correlation between the OH stretch frequency shift of the protonated methanol cluster, H+(MeOH)n, upon the π-hydrogen bond formation with benzene and the enthalpy change in clustering of H+(MeOH)n to H+(MeOH)n+1. This result suggests a new method to explore hydrogen bond strength in hydrogen bond networks.

Evolution of vibrational bands upon gradual protonation/deprotonation of arsinic acid H2As(O)OH in media of different polarity

This computational work is devoted to the investigation (MP2/def2-TZVP) of the geometry and IR parameters of arsinic acid H₂AsOOH and its hydrogen-bonded complexes under vacuum and in media with different polarity. The medium effects were accounted for in two ways: (1) implicitly, using the IEFPCM model, varying the dielectric permittivity (ε) and (2) explicitly, by considering hydrogen-bonded complexes of H₂As(O)OH with various hydrogen bond donors (41 complexes) or acceptors (38 complexes), imitating a gradual transition to the As(OH)₂⁺ or AsO₂^- moiety, respectively. It was shown that the transition from vacuum to a medium with ε > 1 causes the As(O)OH fragment to lose its flatness. The solvent polar medium introduces significant changes in the geometry and IR spectral parameters of hydrogen-bonded complexes too: as the polarity of a medium increases, weak hydrogen bonds become weaker, and strong and medium hydrogen bonds become stronger; in the case of a complex with two hydrogen bonds cooperativity effects are observed. In almost all cases the driving force of these changes appears to be preferential solvation of charge-separated structures. In the limiting case of complete deprotonation (or conversely complete protonation) the vibrational frequencies of ν_AsO and ν_As-O turn into ν_As-O(asym) and ν_As-O(sym), respectively. In the intermediate cases the distance between ν_AsO and ν_As-O is sensitive to both implicit solvation and explicit solvation and the systematic changes of this distance can be used for estimation of the degree of proton transfer within the hydrogen bond.

Estimations of OH·N hydrogen bond length from positions and intensities of IR bands

In this computational work applicability of IR spectral parameters for evaluations of OH···N hydrogen bond length is discussed. For a set of 124 complexes with OH···N hydrogen bond formed by combinations of methanol/acetic acid and pyridine (and their fluorine substituted versions) geometries, energies and IR parameters were calculated at MP2/def2-TZVP level of theory. For a number of IR parameters (the shift of proton donor group stretching vibration Δν_s, increase of its intensity I, the low-frequency hydrogen bond stretching vibration ν_σ, bending in-plane δ and out-of-plane γ vibrations) equations linking them with interatomic distances are proposed, the robustness and accuracy of such equations are discussed. The enthalpy of OH···N hydrogen bond formation ΔH was also linked with electron density parameters in (3; -1) critical point.

Simultaneous Estimation of Two Coupled Hydrogen Bond Geometries from Pairs of Entangled NMR Parameters: The Test Case of 4-Hydroxypyridine Anion

The computational method for estimating the geometry of two coupled hydrogen bonds with geometries close to linear using a pair of spectral NMR parameters was proposed. The method was developed based on the quantum-chemical investigation of 61 complexes with two hydrogen bonds formed by oxygen and nitrogen atoms of the 4-hydroxypyridine anion with OH groups of substituted methanols. The main idea of the method is as follows: from the NMR chemical shifts of nuclei of atoms forming the 4-hydroxylpyridine anion, we select such pairs, whose values can be used for simultaneous determination of the geometry of two hydrogen bonds, despite the fact that every NMR parameter is sensitive to the geometry of each of the hydrogen bonds. For these parameters, two-dimensional maps of dependencies of NMR chemical shifts on interatomic distances in two hydrogen bonds were constructed. It is shown that, in addition to chemical shifts of the nitrogen atom and quaternary carbon, which are experimentally difficult to obtain, chemical shifts of the carbons and protons of the CH groups can be used. The performance of the proposed method was evaluated computationally as well on three additional complexes with substituted alcohols. It was found that, for all considered cases, hydrogen bond geometries estimated using two-dimensional correlations differed from those directly calculated by quantum-chemical methods by not more than 0.04 Å.

Phosphine oxides as NMR and IR spectroscopic probes for geometry and energy of PO···H–A hydrogen bonds

In this work we evaluate the possibility to use the NMR and IR spectral properties of P=O group to estimate the geometry and strength of hydrogen bonds which it forms...

Hydration of selenolate moiety: Ab initio investigation of properties of O-H⋯Se(-) hydrogen bonds in CH3 Se(-)⋯(H2 O)n clusters

This work is devoted to investigations of the influence of O-H···Se(-) hydrogen bonds on the electronic shells of selenolate R-Se(-) fragment (R═CH₃ ). The geometric, energetic and nuclear magnetic resonance (NMR) spectral parameters for various conformers of CH₃ Se(-)⋯(H₂ O)_n clusters with n = 0-6 were calculated at CCSD/aug-cc-pVDZ level of theory. For selenolate anion CH₃ Se(-) solvation free energy was calculated, and for water media it is equal to -71.41 kcal/mol. For O-H···Se(-) hydrogen bonds the proportionality coefficients between QTAIM parameters at (3; -1) bond critical point and the strength of an individual hydrogen bond ∆E were proposed. It was shown, that despite a relative weakness of O-H···Se(-) hydrogen bonds, the outer electronic shell of the selenium atom changes significantly upon formation of each hydrogen bond. This, in turn, cause the dramatic change of NMR parameters of selenium nuclei.

On the influence of water molecules on the outer electronic shells of R-SeH, R-Se(-) and R-SeOH fragments in the selenocysteine amino acid residue

In this computational work (MP2/aug-cc-pVTZ) we investigated the features of the outer electronic shells of R-SeH, R-Se(-) and R-SeOH fragments (R = CH3), which can be considered as simplified models for the forms of the active centres of glutathione peroxidases GPx along their catalytic pathway (reduction of peroxides). It is shown that the preferential direction of a nucleophilic attack on the R-Se(-) fragment by a peroxide molecule is determined by the presence of the electron-depleted region of the selenium atom in front of the C-Se bond and nucleophilic attack can be facilitated by the solvation of R-Se(-) by water molecules. Such solvation does not block the direction of potential nucleophilic attack and also leads to the increase of the maximal value of the molecular electrostatic potential on the selenium atom. It was shown that the 77Se NMR chemical shift is sensitive both to the oxidation state and the hydration state of the selenium-containing fragment.

Actual Symmetry of Symmetric Molecular Adducts in the Gas Phase, Solution and in the Solid State

This review discusses molecular adducts, whose composition allows a symmetric structure. Such adducts are popular model systems, as they are useful for analyzing the effect of structure on the property selected for study since they allow one to reduce the number of parameters. The main objectives of this discussion are to evaluate the influence of the surroundings on the symmetry of these adducts, steric hindrances within the adducts, competition between different noncovalent interactions responsible for stabilizing the adducts, and experimental methods that can be used to study the symmetry at different time scales. This review considers the following central binding units: hydrogen (proton), halogen (anion), metal (cation), water (hydrogen peroxide).

Estimations of FH···X hydrogen bond energies from IR intensities: Iogansen's rule revisited

In this work the possibility of using the IR intensity of the stretching vibration ν_s of proton donor group for estimation of hydrogen bond strength was investigated. For a set of complexes with FH···X (X = F, N, O) hydrogen bonds in the wide range of energies (0.1-49.2 kcal/mol) vibrational frequencies ν_s and their intensities A were calculated (CCSD at complete basis set limit). The validity of the previously proposed linear proportionality between the intensification of the stretching vibration ν_s in IR spectra and hydrogen bond enthalpy -ΔH = 12.2 ∆ A (A. V. Iogansen, Spectrochimica Acta A 1999) was examined. It is shown that for a range of similar hydrogen bond types with complexation energies ∆E <15 kcal/mol the ∆E( ∆ A ) function remains similar to that proposed in the Iogansen's work, while upon strengthening this dependency becomes significantly nonlinear. We examined two other parameters ( ∆ A ν s and ∆ A ∙ m R ) related to IR intensity as descriptors of hydrogen bond strength which are proportional to transition dipole moment matrix element and mass-independent dipole moment derivative. It was found that the dependency ∆E( ∆ A ν s ) stays linear in the whole studied range of complexation energies and it can be used for evaluation of ∆E from infrared spectral data with the accuracy about 2 kcal/mol. The mass-independent product ∆ A ∙ m R is an appropriate descriptor for sets of complexes with various hydrogen bond types. Simple equations proposed in this work can be used for estimations of hydrogen bond strength in various systems, where experimental thermodynamic methods or direct calculations are difficult or even impossible.

O-H stretching frequency red shifts do not correlate with the dissociation energies in the dimethylether and dimethylsulfide complexes of phenol derivatives

In this perspective, we present a comprehensive report on the spectroscopic and computational investigations of the hydrogen bonded (H-bonded) complexes of Me2O and Me2S with seven para-substituted H-bond donor phenols. The salient finding was that although the dissociation energies, D0, of the Me2O complexes were consistently higher than those of the analogous Me2S complexes, the red-shifts in phenolic O-H frequencies, Δν(O-H), showed the exactly opposite trend. This is in contravention of the general perception that the red shift in the X-H stretching frequency in the X-HY hydrogen bonded complexes is a reliable indicator of H-bond strength (D0), a concept popularly known as the Badger-Bauer rule. This is also in contrast to the trend reported for the H-bonded complexes of H2S/H2O with several para substituted phenols of different pKa values wherein the oxygen centered hydrogen bonded (OCHB) complexes consistently showed higher Δν(O-H) and D0 compared to those of the analogous sulfur centered hydrogen bonded (SCHB) complexes. Our effort was to understand these intriguing observations based on the spectroscopic investigations of 1 : 1 complexes in combination with a variety of high level quantum chemical calculations. Ab initio calculations at the MP2 level and the DFT calculations using various dispersion corrected density functionals (including DFT-D3) were performed on counterpoise corrected surfaces to compute the dissociation energy, D0, of the H-bonded complexes. The importance of anharmonic frequency computations is underscored as they were able to correctly reproduce the observed trend in the relative OH frequency shifts unlike the harmonic frequency computations. We have attempted to find a unified correlation that would globally fit the observed red shifts in the O-H frequency with the H-bonding strength for the four bases, namely, H2S, H2O, Me2O, and Me2S, in this set of H-bond donors. It was found that the proton affinity normalized Δν(O-H) values scale very well with the H-bond strength.

Enhancing Structural Rigidity via a Strategy Involving Protons for Creating Water-Resistant Mn4+-Doped Fluoride Phosphors

The poor water resistance property of a commercial Mn⁴⁺-activated narrow-band red-emitting fluoride phosphor restricts its promising applications in high-performance white LEDs and wide-gamut displays. Herein, we develop a structural rigidity-enhancing strategy using a novel KHF₂:Mn⁴⁺ precursor as a Mn source to construct a proton-containing water-resistant phosphor K₂(H)TiF₆:Mn⁴⁺ (KHTFM). The parasitic [HMnF₆]^- complexes in the interstitial site from the fall off the KHF₂:Mn⁴⁺ are also transferred to the K₂TiF₆ host by ion exchange to form KHTFM with rigid bonding networks, improving the water resistance and thermostability of the sample. The KHTFM sample retains at least 92% of the original emission value after 180 min of water immersion, while the non-water-resistant K₂TiF₆:Mn⁴⁺(KTFM) phosphor maintains only 23%. Therefore, these findings not only illustrate the effect of protons on fluoride but also provide a novel insight into commercial water-resistant fluoride phosphors.

Relationships between NMR shifts and interaction energies in biphenyls, alkanes, aza-alkanes, and oxa-alkanes with X─H… H─Y and X─H… Z (X, Y = C or N; Z = N or O) hydrogen bonding

Hydrogen-hydrogen C─H^… H─C bonding between the bay-area hydrogens in biphenyls, and more generally in congested alkanes, very strained polycyclic alkanes, and cis-2-butene, has been investigated by calculation of proton nuclear magnetic resonance (NMR) shifts and atom-atom interaction energies. Computed NMR shifts for all protons in the biphenyl derivatives correlate very well with experimental data, with zero intercept, unit slope, and a root mean square deviation of 0.06 ppm. For some congested alkanes, there is generally good agreement between computed values for a selected conformer and the experimental data, when it is available. In both cases, the shift of a given proton or pair of protons tends to increase with the corresponding interaction energy. Computed NMR shift differences for methylene protons in polycyclic alkanes, where one is involved in a very short contact ("in") and the other is not ("out"), show a rough correlation with the corresponding C─H^… H─C exchange energies. The "in" and "in,in" isomers of selected aza- and diaza-cycloalkanes, respectively, are X─H^… H─N hydrogen bonded, whereas the "out" and "in,out" isomers display X─H^… N hydrogen bonds (X = C or N). Oxa-alkanes and the "in" isomers of aza-oxa-alkanes are X─H^… O hydrogen bonded. There is a very good general correlation, including both N─H^… H─Y (Y = C or N) and N─H^… Z (Z = N or O) interactions, for NH proton shifts against the exchange energy. For "in" CH protons, the data for the different C─H^… H─Y and C─H^… Z interactions are much more dispersed and the overall shift/exchange energy correlation is less satisfactory.

Correlations of NHN hydrogen bond energy with geometry and 1H NMR chemical shift difference of NH protons for aniline complexes

In this computational work, we propose to use the NMR chemical shift difference of NH₂ protons for 1:1 complexes formed by aniline and nitrogen-containing proton acceptors for the estimation of the hydrogen bond energy and geometry (N⋯H and N⋯N distances). The proposed correlations could be applied to other aromatic amines as well, in a gas phase, a solution, or a solid state, for both inter- and intramolecular hydrogen bonds. We considered a set of 21 complexes with the NHN hydrogen bond without proton transfer, including hydrogen bonds from weak to medium strong ones (2-21 kcal/mol), with neutral or anionic bases and with sp³ and sp² hybridized nitrogen proton acceptors. For each complex apart from direct hydrogen bond energy calculation, we have tested several other ways to estimate the energy: (a) using a correlation between NH stretching band intensity and hydrogen bond energy and (b) using correlations between electron density properties at (3, -1) bond critical point (quantum theory of atoms in molecules analysis) and hydrogen bond energy. Besides for the studied type of complexes, we obtained refined linear correlations linking the local electron kinetic (G) and potential (V) energy densities with the hydrogen bond energy.