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.

The Cambridge Structural Database and structural dynamics

With the availability of the computer readable information in the Cambridge Structural Database (CSD), wide ranging, largely automated comparisons of fragment, molecular, and crystal structures have become possible. They show that the distributions of interatomic distances, angles, and torsion angles for a given structural fragment occurring in different environments are highly correlated among themselves and with other observables such as spectroscopic signals, reaction and activation energies. The correlations often extend continuously over large ranges of parameter values. They are reminiscent of bond breaking and forming reactions, polyhedral rearrangements, and conformational changes. They map-qualitatively-the regions of the structural parameter space in which molecular dynamics take place, namely, the low energy regions of the respective (free) energy surfaces. The extension and continuous nature of the correlations provides an organizing principle of large groups of structural data and suggests a reconsideration of traditional definitions and descriptions of bonds, "nonbonded" and "noncovalent" interactions in terms of Lewis acids interacting with Lewis bases. These aspects are illustrated with selected examples of historic importance and with some later developments. It seems that the amount of information in the CSD (and other structural databases) and the knowledge on the nature of, and the correlations within, this body of information should allow one-in the near future-to make credible interpolations and possibly predictions of structures and their properties with machine learning methods.

Effects of charge-assisted hydrogen bond on sorption and co-sorption of pharmaceutical contaminants on carbonaceous materials: Spectroscopic and theoretical studies

Charge-assisted hydrogen bond (CAHB) is a key mechanism that affects the environmental behavior of pharmaceutical pollutants (PCs). However, the strength and stability of various CAHBs, and their effects on the co-sorption behavior of PCs are still unclear. Herein, DFT calculation with different solvent models including two implicit solvent model (PCM and SMD), and one explicit solvent model (ESM) were applied in this study, to investigate the effects of different CAHBs on the sorption and co-sorption behavior of four PCs (e.g., clofibric acid, p-aminobenzoic acid, acetaminophen, and sulfamerazine) on three model carbonaceous materials. First, the appearance of new peaks in the very low field of 1H NMR, and the blue shift of OH and NH2 peaks in FTIR indicated that CAHBs were indeed formed between PCs and carbonaceous materials. Next, according to the principal component analysis and correlation analysis of parameters (e.g., ΔEads, bond length, bond angle, Egap, and ΔG) of these CAHBs calculated by the DFT with different solvent models, the results showed that SMD is the optimal model for calculating the strength and stability of CAHBs by DFT, and the strength and stability of CAHBs formed between PCs and carbonaceous materials in this study were in the order of homonuclear [O⋯H⋯O]- CAHB > heteronuclear [O⋯HN]-/[N⋯HO]+ type of CAHB > homonuclear [N⋯H⋯N]+. Also, the co-sorption behavior of different PCs co-existing in binary systems further confirmed that, all above types of CAHBs formed between PCs and carbonaceous materials can produce obvious competition effect on the co-existing PCs that only OHB formed between them. This study not only reveals the environmental behavior of co-existing PCs, but also provides a theoretical basis for the design of obligate sorption materials for PCs in the natural environment.

Evidence of extremely short hydrogen bond in the homoconjugated ferrocene-1,1'-diyl-bisphosphinic acid anion: sign change of the H/D isotope effect on the 31P NMR chemical shift

The structure of the two intramolecular hydrogen-bonded motifs within ferrocene-1,1'-diyl-bisphosphinic acid, a cyclic dimer and a homoconjugated anion, has been experimentally revealed by NMR spectroscopy for a solution in a low-freezing polar aprotic CDF3/CDF2Cl medium at 100 K. Structure elucidation was achieved with the help of the H/D isotope effects on the 1H, 2H and, for the first time, 31P NMR chemical shifts. The questions of bridging hydron localization and origins of opposite signs of H/D isotope effects on the 31P NMR chemical shifts in the cyclic dimer and homoconjugated anion have also been addressed by DFT calculations, including numerical solution of the Schrödinger equation for proton and deuteron vibrations in the anharmonic double well potentials.

Coexistence of Cationic and Anionic Phosphate Moieties in Solids: Unusual but Not Impossible

Phosphoric acid is commonly known either as a neutral molecule or as an anion (phosphate). We theoretically confirm by ab initio molecular dynamics simulations (AIMD) that a cationic form H4PO4+ coexists with the anionic form H2PO4- in the same salt. This paradoxical situation is achieved by partial substitution of Cs+ by H4PO4+ in CsH2PO4. Thus, HnPO4 acts simultaneously as both the positive and the negative ion of the salt. We analyze the dynamical protonation pattern within the unusual hydrogen bond network that is established between the ions. Our AIMD simulations show that a conventional assignment of protonation states of the phosphate groups is not meaningful. Instead, a better description of the protonation situation is achieved by an efficiently fractional assignment of the strongly hydrogen-bonded protons to both its nearest and next-nearest oxygen neighbors.

Weak, Broken, but Working-Intramolecular Hydrogen Bond in 2,2'-bipyridine

From an academic and practical point of view, it is desirable to be able to assess the possibility of the proton exchange of a given molecular system just by knowing the positions of the proton acceptor and the proton donor. This study addresses the difference between intramolecular hydrogen bonds in 2,2'-bipyridinium and 1,10-phenanthrolinium. Solid-state ¹⁵N NMR and model calculations show that these hydrogen bonds are weak; their energies are 25 kJ/mol and 15 kJ/mol, respectively. Neither these hydrogen bonds nor N-H stretches can be responsible for the fast reversible proton transfer observed for 2,2'-bipyridinium in a polar solvent down to 115 K. This process must have been caused by an external force, which was a fluctuating electric field present in the solution. However, these hydrogen bonds are the grain that tips the scales precisely because they are an integral part of a large system of interactions, including both intramolecular interactions and environmental influence.

Acid Dissociation in (HX)n (H2 O)n Clusters (X=F, Cl, Br, I; n=2, 3)

In this work, we analyze the interactions between two or three hydrogen halide molecules and the same number of water moieties through a systematic exploration of their potential energy surfaces. Our results indicate that the most stable HF and HCl aggregates do not experience dissociation of any of the acid fragments, even with three water molecules. In contrast, in the HBr and HI clusters, one of the acid fragments does dissociate. While the global minimum of (HBr)₃ (H₂ O)₃ is a hydrogen-bridged bihalide anion (BrHBr^- ), which is persistent at temperatures up to 203 K, the lowest energy structure of (HI)₃ (H₂ O)₃ has a separated ion pair, but the motif with a bihalide anion (IHI^- ) is only 0.2 kcal mol^-1 above the global minimum. Among the more stable structures is a broad spectrum of contacts, including water⋯water, HX⋯water, and HX⋯HX hydrogen bonds, halogen bonds, ionic and long-range X⋯H contacts.

From Local Covalent Bonding to Extended Electric Field Interactions in Proton Hydration

Seemingly simple yet surprisingly difficult to probe, excess protons in water constitute complex quantum objects with strong interactions with the extended and dynamically changing hydrogen-bonding network of the liquid. Proton hydration plays pivotal roles in energy transport in hydrogen fuel cells and signal transduction in transmembrane proteins. While geometries and stoichiometry have been widely addressed in both experiment and theory, the electronic structure of these specific hydrated proton complexes has remained elusive. Here we show, layer by layer, how utilizing novel flatjet technology for accurate x-ray spectroscopic measurements and combining infrared spectral analysis and calculations, we find orbital-specific markers that distinguish two main electronic-structure effects: Local orbital interactions determine covalent bonding between the proton and neigbouring water molecules, while orbital-energy shifts measure the strength of the extended electric field of the proton.

Proton-mediated photoprotection mechanism in photosystem II

Photo-induced charge separation, which is terminated by electron transfer from the primary quinone Q_A to the secondary quinone Q_B, provides the driving force for O₂ evolution in photosystem II (PSII). However, the backward charge recombination using the same electron-transfer pathway leads to the triplet chlorophyll formation, generating harmful singlet-oxygen species. Here, we investigated the molecular mechanism of proton-mediated Q_A ^⋅- stabilization. Quantum mechanical/molecular mechanical (QM/MM) calculations show that in response to the loss of the bicarbonate ligand, a low-barrier H-bond forms between D2-His214 and Q_A ^⋅-. The migration of the proton from D2-His214 toward Q_A ^⋅- stabilizes Q_A ^⋅-. The release of the bicarbonate ligand from the binding Fe²⁺ site is an energetically uphill process, whereas the bidentate-to-monodentate reorientation is almost isoenergetic. These suggest that the bicarbonate protonation and decomposition may be a basis of the mechanism of photoprotection via Q_A ^⋅-/Q_AH^⋅ stabilization, increasing the Q_A redox potential and activating a charge-recombination pathway that does not generate the harmful singlet oxygen.

Direct Spectroscopic Evidence for Charge-Assisted Hydrogen-Bond Formation between Ionizable Organic Chemicals and Carbonaceous Materials

The direct evidence for the formation of charge-assisted hydrogen bond (CAHB) between the charged groups of ionizable organic chemicals (IOCs) and carbonaceous materials with similar proton affinity remains elusive. We therefore selected three pharmaceutical contaminants (PCs) as representative IOCs to provide the direct evidence of CAHB formation between IOCs and functionalized carbon nanotubes (CNTs) and its intensity/contribution to PC sorption on CNTs by NMR, FTIR, and DFT analyses. Sorption of PCs on functionalized CNTs resulted in the FTIR characteristic peak that appeared at a higher frequency (3787 cm^-1) and the ¹H NMR characteristic peak that emerged at an extremely low-field region (<18.0 ppm), which can be used as the direct spectroscopic evidence for CAHB formation. Both homonuclear CAHB (HM-CAHB, e.g., [O-H···O]^-) and heteronuclear CAHB (HT-CAHB, e.g., [N⁺-H···O^-]/[O-H···N]⁺) exhibited a much higher sorption energy (|E_ads| ≥ 56.24 kJ/mol) than ordinary hydrogen bond (OHB, |E_ads| ≤ 6.136 kJ/mol), leading to a greater sorption contribution (HM-/HT-CAHB ≥ 42.3%, OHB ≤ 36.5%) and irreversibility (hysteresis index: HM-/HT-CAHB ≥ 1.69, OHB ≤ 0.43) of PCs on CNTs. This work presents the direct evidence for CAHB formation between IOCs and CNTs, which is significant for understanding and predicting the environmental fate and risk of IOCs, thus providing new insights for controlling their pollution using specifically designed carbonaceous materials.

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

A phosphonic acid anion and acid dimer dianion stabilized by proton transfer in OHN hydrogen bonds - models of structural motifs in blend polymer membranes

The structure of intermolecular hydrogen-bonded complexes formed between tert-butylphosphonic acid and trimethylpyridine molecules has been experimentally studied as the simplest model system of the structural motifs in blend proton-conducting polymer membranes based on phosphonic acid residues. The stoichiometry of the formed complexes and proton positions in OHO and OHN hydrogen bonds were established by the H/D isotope effects and temperature dependences of the signals in ¹H and ³¹P NMR spectra. Two structural motifs, namely, 1 : 2 and 2 : 2 acid-base complexes, were identified at the low temperature in a polar aprotic environment. In the 1 : 2 complex, one proton has passed through the hydrogen bond center creating a chain of two cooperatively coupled OHN bonds, while in the 2 : 2 complex both OHN bonds are zwitterionic and anti-cooperatively coupled to each other via a dianionic cyclic dimer of phosphonic acid in the middle. The dianionic cyclic dimer is metastable by itself, but under the used experimental conditions it is stabilized by complexation with two trimethylpyridinium cations. Additionally, quantum chemical calculations using the DFT method were carried out to support the experimental data.

Self-association of diphenylpnictoginic acids in solution and solid state: covalent vs. hydrogen bonding

Triphenylpnictogens were oxidized to access diphenylpnictioginic acids Ph₂XOOH (X = P, As, Sb, Bi). It was shown that oxidation with chloramine-T does not lead to the cleavage of a C-pnictogen bond. The preliminary reductive cleavage with sodium in liquid ammonia followed by the oxidation with hydrogen peroxide was successfully utilised for the synthesis of diphenylphosphinic and diphenylarsinic acids. It was shown that in solid state (by means of XRD), all diphenylpnictoginic acids form polymeric chains. Diphenylbismuthinic and diphenylantimonic acids form polymeric covalent adducts, while diphenylphosphinic and diphenylarsinic chains are associated through hydrogen bonding. Unlike diphenylphosphinic acid, diphenilarsinic acid forms two polymorphs of hydrogen-bonded infinite chains. In solution in a polar aprotic solvent diphenylarsinic acid, similarly to dimethylarsinic, forms hydrogen-bonded cyclic dimers together with a small amount of cyclic trimers.

An N⋯H⋯N low-barrier hydrogen bond preorganizes the catalytic site of aspartate aminotransferase to facilitate the second half-reaction

Pyridoxal 5′-phosphate (PLP)-dependent enzymes have been extensively studied for their ability to fine-tune PLP cofactor electronics to promote a wide array of chemistries. Neutron crystallography offers a straightforward approach to studying the electronic states of PLP and the electrostatics of enzyme active sites, responsible for the reaction specificities, by enabling direct visualization of hydrogen atom positions. Here we report a room-temperature joint X-ray/neutron structure of aspartate aminotransferase (AAT) with pyridoxamine 5′-phosphate (PMP), the cofactor product of the first half reaction catalyzed by the enzyme. Between PMP N_SB and catalytic Lys258 Nζ amino groups an equally shared deuterium is observed in an apparent low-barrier hydrogen bond (LBHB). Density functional theory calculations were performed to provide further evidence of this LBHB interaction. The structural arrangement and the juxtaposition of PMP and Lys258, facilitated by the LBHB, suggests active site preorganization for the incoming ketoacid substrate that initiates the second half-reaction.

The neutron structure of pyridoxal 5′-phosphate-dependent enzyme aspartate aminotransferase with pyridoxamine 5′-phosphate (PMP) reveals a low-barrier hydrogen bond between the amino groups of PMP and catalytic Lys258, preorganizing the active site for catalysis

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

Dissociation of HCl in water nanoclusters: an energy decomposition analysis perspective

As known, small HCl-water nanoclusters display a particular dissociation behaviour, whereby at least four water molecules are required for the ionic dissociation of HCl. In this work, we examine how intermolecular interactions promote the ionic dissociation of such nanoclusters. To this end, a set of 45 HCl-water nanoclusters with up to four water molecules is introduced. Energy decomposition analysis based on absolutely localized molecular orbitals (ALMO-EDA) is employed in order to study the importance of frozen interaction, dispersion, polarization, and charge-transfer for the dissociation. The vertical ALMO-EDA scheme is applied to HCl-water clusters along a proton-transfer coordinate varying the amount of spectator water molecules. The corresponding ALMO-EDA results show a clear preference for the dissociated cluster only in the presence of four water molecules. Our analysis of adiabatic ALMO-EDA results reveals a push-pull mechanism for the destabilization of the HCl bond based on the synergy between forward and backward charge-transfer.

Modeling of the Response of Hydrogen Bond Properties on an External Electric Field: Geometry, NMR Chemical Shift, Spin-Spin Scalar Coupling

The response of the geometric and NMR properties of molecular systems to an external electric field has been studied theoretically in a wide field range. It has been shown that this adduct under field approach can be used to model the geometric and spectral changes experienced by molecular systems in polar media if the system in question has one and only one bond, the polarizability of which significantly exceeds the polarizability of other bonds. If this requirement is met, then it becomes possible to model even extreme cases, for example, proton dissociation in hydrogen halides. This requirement is fulfilled for many complexes with one hydrogen bond. For such complexes, this approach can be used to facilitate a detailed analysis of spectral changes associated with geometric changes in the hydrogen bond. For example, in hydrogen-bonded complexes of isocyanide C≡¹⁵N-¹H⋯X, ¹J(¹⁵N¹H) depends exclusively on the N-H distance, while δ(¹⁵N) is also slightly influenced by the nature of X.

Classification of So-Called Non-Covalent Interactions Based on VSEPR Model

The variety of interactions have been analyzed in numerous studies. They are often compared with the hydrogen bond that is crucial in numerous chemical and biological processes. One can mention such interactions as the halogen bond, pnicogen bond, and others that may be classified as σ-hole bonds. However, not only σ-holes may act as Lewis acid centers. Numerous species are characterized by the occurrence of π-holes, which also may play a role of the electron acceptor. The situation is complicated since numerous interactions, such as the pnicogen bond or the chalcogen bond, for example, may be classified as a σ-hole bond or π-hole bond; it ultimately depends on the configuration at the Lewis acid centre. The disadvantage of classifications of interactions is also connected with their names, derived from the names of groups such as halogen and tetrel bonds or from single elements such as hydrogen and carbon bonds. The chaos is aggravated by the properties of elements. For example, a hydrogen atom can act as the Lewis acid or as the Lewis base site if it is positively or negatively charged, respectively. Hence names of the corresponding interactions occur in literature, namely hydrogen bonds and hydride bonds. There are other numerous disadvantages connected with classifications and names of interactions; these are discussed in this study. Several studies show that the majority of interactions are ruled by the same mechanisms related to the electron charge shifts, and that the occurrence of numerous interactions leads to specific changes in geometries of interacting species. These changes follow the rules of the valence-shell electron-pair repulsion model (VSEPR). That is why the simple classification of interactions based on VSEPR is proposed here. This classification is still open since numerous processes and interactions not discussed in this study may be included within it.

Double Proton Tautomerism via Intra- or Intermolecular Pathways? The Case of Tetramethyl Reductic Acid Studied by Dynamic NMR: Hydrogen Bond Association, Solvent and Kinetic H/D Isotope Effects

Using dynamic liquid-state NMR spectroscopy a degenerate double proton tautomerism was detected in tetramethyl reductic acid (TMRA) dissolved in toluene-d₈ and in CD₂Cl₂. Similar to vitamin C, TMRA belongs to the class of reductones of biologically important compounds. The tautomerism involves an intramolecular HH transfer that interconverts the peripheric and the central positions of the two OH groups. It is slow in the NMR time scale around 200 K and fast at room temperature. Pseudo-first-order rate constants of the HH transfer and of the HD transfer after suitable deuteration were obtained by line shape analyses. Interestingly, the chemical shifts were found to be temperature dependent carrying information about an equilibrium between a hydrogen bonded dimer and a monomer forming two weak intramolecular hydrogen bonds. The structures of the monomer and the dimer are discussed. The latter may consist of several rapidly interconverting hydrogen-bonded associates. A way was found to obtain the enthalpies and entropies of dissociation, which allowed us to convert the pseudo-first-order rate constants of the reaction mixture into first-order rate constants of the tautomerization of the monomer. Surprisingly, these intrinsic rate constants were the same for toluene-d₈ and CD₂Cl₂, but in the latter solvent more monomer is formed. This finding is attributed to the dipole moment of the TMRA monomer, compensated in the dimer, and to the larger dielectric constant of CD₂Cl₂. Within the margin of error, the kinetic HH/HD isotope effects were found to be of the order of 3 but independent of temperature. That finding indicates a stepwise HH transfer involving a tunnel mechanism along a double barrier pathway. The Arrhenius curves were described in terms of the Bell–Limbach tunneling model.

Combined NMR and UV-Vis Spectroscopic Studies of Models for the Hydrogen Bond System in the Active Site of Photoactive Yellow Protein: H-Bond Cooperativity and Medium Effects

Intramolecular hydrogen bonds in aprotic media were studied by combined (simultaneous) NMR and UV-vis spectroscopy. The species under investigation were anionic and featured single or coupled H-bonds between, for example, carboxylic groups and phenolic oxygen atoms (COO···H···OC)^-, among phenolic oxygen atoms (CO···H···OC)^-, and hydrogen bond chains between a carboxylic group and two phenolic oxygen atoms (COO···H···(OC)···H···OC)^-. The last anion may be regarded as a small molecule model for the hydrogen bond system in the active site of wild-type photoactive yellow protein (PYP) while the others mimic the corresponding H-bonds in site-selective mutants. Proton positions in isolated hydrogen bonds and hydrogen bond chains were assessed by calculations for vacuum conditions and spectroscopically for the two media, CD₂Cl₂ and the liquefied gas mixture CDClF₂/CDF₃ at low temperatures. NMR parameters allow for the estimation of time-averaged H-bond geometries, and optical spectra give additional information about geometry distributions. Comparison of the results from the various systems revealed the effects of the formation of hydrogen bond chains and changes of medium conditions on the geometry of individual H-bonds. In particular, the proton in a hydrogen bond to a carboxylic group shifts from the phenolic oxygen atom in the system COO^-···H-OC to the carboxylic group in COO-H···(OC)^-···H-OC as a result of hydrogen bond formation to the additional phenolic donor. Increase in medium polarity may, however, induce the conversion of a structure of a type COO-H···(OC)^-···H-OC to the type COO^-···H-(OC)···H-OC. Application of these results obtained from the model systems to PYP suggests that both cooperative effects within the hydrogen bond chain and a low-polarity protein environment are prerequisites for the stabilization of negative charge on the cofactor and hence for the spectral tuning of the photoreceptor.

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

Infrared and NMR Spectroscopic Fingerprints of the Asymmetric H7 + O3 Complex in Solution

Infrared (IR) absorption in the 1000–3700 cm⁻¹ range and ¹H NMR spectroscopy reveal the existence of an asymmetric protonated water trimer, H₇⁺O_3, in acetonitrile. The core H₇⁺O₃ motif persists in larger protonated water clusters in acetonitrile up to at least 8 water molecules. Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations reveal irreversible proton transport promoted by propagating the asymmetric H₇⁺O₃ structure in solution. The QM/MM calculations allow for the successful simulation of the measured IR absorption spectra of H₇⁺O₃ in the OH stretch region, which reaffirms the assignment of the H₇⁺O₃ spectra to a hybrid‐complex structure: a protonated water dimer strongly hydrogen‐bonded to a third water molecule with the proton exchanging between the two possible shared‐proton Zundel‐like centers. The H₇⁺O₃ structure lends itself to promoting irreversible proton transport in presence of even one additional water molecule. We demonstrate how continuously evolving H₇⁺O₃ structures may support proton transport within larger water solvates.

Self-Assembly of Hydrogen-Bonded Cage Tetramers of Phosphonic Acid

The self-association of phosphonic acids with general formula RP(O)(OH)2 in solution state remains largely unexplored. The general understanding is that such molecules form multiple intermolecular hydrogen bonds, but the stoichiometry of self-associates and the bonding motifs are unclear. In this work, we report the results of the study of self-association of tert-butylphosphonic acid using low temperature liquid-state 1H and 31P NMR spectroscopy (100 K; CDF3/CDF2Cl) and density functional theory (DFT) calculations. For the first time, we demonstrate conclusively that polar aprotic medium tert-butylphosphonic acid forms highly symmetric cage-like tetramers held by eight OHO hydrogen bonds, which makes the complex quite stable. In these associates. each phosphonic acid molecule is bonded to three other molecules by forming two hydrogen bonds as proton donor and two hydrogen bonds as proton acceptor. Though the structure of such cage-like tetramers is close to tetrahedral, the formal symmetry of the self-associate is C2.

Hydrogen Bond and Other Lewis Acid-Lewis Base Interactions as Preliminary Stages of Chemical Reactions

Various Lewis acid-Lewis base interactions are discussed as initiating chemical reactions and processes. For example, the hydrogen bond is often a preliminary stage of the proton transfer process or the tetrel and pnicogen bonds lead sometimes to the SN2 reactions. There are numerous characteristics of interactions being first stages of reactions; one can observe a meaningful electron charge transfer from the Lewis base unit to the Lewis acid; such interactions possess at least partly covalent character, one can mention other features. The results of different methods and approaches that are applied in numerous studies to describe the character of interactions are presented here. These are, for example, the results of the Quantum Theory of Atoms in Molecules, of the decomposition of the energy of interaction or of the structure-correlation method.

H/D Isotope Effects on 1H-NMR Chemical Shifts in Cyclic Heterodimers and Heterotrimers of Phosphinic and Phosphoric Acids

Hydrogen-bonded heterocomplexes formed by POOH-containing acids (diphenylphosphoric 1, dimethylphosphoric 2, diphenylphosphinic 3, and dimethylphosphinic 4) are studied by the low-temperature (100 K) ¹H-NMR and ³¹P-NMR using liquefied gases CDF₃/CDF₂Cl as a solvent. Formation of cyclic dimers and cyclic trimers consisting of molecules of two different acids is confirmed by the analysis of vicinal H/D isotope effects (changes in the bridging proton chemical shift, δH, after the deuteration of a neighboring H-bond). Acids 1 and 4 (or 1 and 3) form heterotrimers with very strong (short) H-bonds (δH ca. 17 ppm). While in the case of all heterotrimers the H-bonds are cyclically arranged head-to-tail, ···O=P-O-H···O=P-O-H···, and thus their cooperative coupling is expected, the signs of vicinal H/D isotope effects indicate an effective anticooperativity, presumably due to steric factors: when one of the H-bonds is elongated upon deuteration, the structure of the heterotrimer adjusts by shortening the neighboring hydrogen bonds. We also demonstrate the formation of cyclic tetramers: in the case of acids 1 and 4 the structure has alternating molecules of 1 and 4 in the cycle, while in case of acids 1 and 3 the cycle has two molecules of 1 followed by two molecules of 3.

Microscopic Difference of Hydrogen Double-minimum Potential Well Detected by Hydroxyl Group in Hydrogen-bonded System

We investigate the microscopic structure of hydrogen double-well potentials in a hydrogen-bonded ferroelectric system exposed to radioactive particles of hydrogen-ion beams. The hydrogen-bonded system is ubiquitous, forming the base of organic-inorganic materials and the double-helix structure of DNA inside biological materials. In order to determine the difference of microscopic environments, an atomic-scale level analysis of solid-state ¹H high-resolution nuclear magnetic resonance (NMR) spectra was performed. The hydrogen environments of inorganic systems represent the Morse potentials and wave function of the eigen state and eigen-state energy derived from the Schrödinger equation. The wave functions for the real space of the localized hydrogen derived from the approximated solutions in view of the atomic scale by using quantum mechanics are manifested by a difference in the charge-density distribution.

Adduct under Field-A Qualitative Approach to Account for Solvent Effect on Hydrogen Bonding

The location of a mobile proton in acid-base complexes in aprotic solvents can be predicted using a simplified Adduct under Field (AuF) approach, where solute–solvent effects on the geometry of hydrogen bond are simulated using a fictitious external electric field. The parameters of the field have been estimated using experimental data on acid-base complexes in CDF₃/CDClF₂. With some limitations, they can be applied to the chemically similar CHCl₃ and CH₂Cl₂. The obtained data indicate that the solute–solvent effects are critically important regardless of the type of complexes. The temperature dependences of the strength and fluctuation rate of the field explain the behavior of experimentally measured parameters.

Symmetry and 1H NMR chemical shifts of short hydrogen bonds: impact of electronic and nuclear quantum effects

Electronic and nuclear quantum effects determine the symmetry and highly downfield ¹H NMR chemical shifts of short hydrogen bonds.

QM/MM simulations of organic phosphorus adsorption at the diaspore-water interface

Phosphorus (P) immobilization and thus its availability for plants are mainly affected by the strong interaction of phosphates with soil components especially soil mineral surfaces. The related reactions have been studied extensively via sorption experiments especially by carrying out adsorption of ortho-phosphates onto Fe-oxide surfaces. But a molecular-level understanding of the P-binding mechanisms at the mineral-water interface is still lacking, especially for forest eco-systems. Therefore, the current contribution provides an investigation of the molecular binding mechanisms for two abundant phosphates in forest soils, inositol hexaphosphate (IHP) and glycerolphosphate (GP), at the diaspore mineral surface. Here a hybrid electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) based molecular dynamics simulation has been applied to explore the diaspore-IHP/GP-water interactions. The results provide evidence for the formation of different P-diaspore binding motifs involving monodentate (M) and bidentate (B) for GP and two (2M) as well as three (3M) monodentates for IHP. The interaction energy results indicated the abundance of the GP B motif compared to the M one. The IHP 3M motif has a higher total interaction energy compared to its 2M motif, but exhibits a lower interaction energy per bond. Compared to GP, IHP exhibited stronger interaction with the surface as well as with water. Water was found to play an important role in controlling these diaspore-IHP/GP-water interactions. The interfacial water molecules form moderately strong H-bonds (HBs) with GP and IHP as well as with the diaspore surface. For all the diaspore-IHP/GP-water complexes, the interaction of water with the diaspore exceeds that with the studied phosphates. Furthermore, some water molecules form covalent bonds with diaspore Al atoms while others dissociate at the surface to protons and hydroxyl groups leading to proton transfer processes. Finally, the current results confirm the previous experimental conclusions indicating the importance of the number of phosphate groups, HBs, and proton transfers in controlling the P-binding at soil mineral surfaces.

Magnetic descriptors of hydrogen bonds in malonaldehyde and its derivatives

The nature of the hydrogen bond, HB, as such is still unknown, though a few of its most fundamental features has been uncovered during the last few decades. At the moment, it is possible to obtain reliable results for only a few of its broadest properties, like magnetic properties. They could give new insights into the physics underlying the strength and features of HBs. In this article we analyze the electronic origin of the NMR spectroscopic parameters of malonaldehyde, MA, and some substituted MAs. These substituted MAs are such that the H-bonds are assisted by one of two phenomena: resonance, RAHB, or charge, CAHB. We have studied the dependences of these parameters on two of the main factors which contribute the most to both phenomena, the geometrical and electronic factors, and found out how they can be used to characterize RAHB or CAHB by means of reliable theoretical calculations. We show that in the set of compounds analyzed here (i) the shielding of the proton of the H-bond can be used as a measure of the strength of the HB and (ii) the relation between the contact and non-contact mechanisms of J-couplings between donor and acceptor atoms is a reliable descriptor of whether the H-bond is resonance assisted or charge assisted.

Conformational Mobility and Proton Transfer in Hydrogen-Bonded Dimers and Trimers of Phosphinic and Phosphoric Acids

The monomers, H-bonded cyclic dimers, and trimers of five acids were studied by density functional theory calculations, such as hypophosphorous acid (H₂POOH, 1), dimethylphosphinic acid (Me₂POOH, 2), phenylphosphinic acid (PhHPOOH, 3), dimethylphosphoric acid ((MeO)₂POOH, 4), and diphenylphosphoric acid ((PhO)₂POOH, 5). Particular attention was paid to the conformational manifold existing due to the internal degrees of freedom: proton transfer (PT), puckering ("twist") within the ring of H-bonds, and mobility of the substituents (namely, -Ph, -OMe, and -OPh rotations). For acid 3, the number of conformers is additionally increased because of the varying relative orientation of nonequivalent substituents in cyclic complexes. We show that ³¹P NMR chemical shifts (δP) are very sensitive to the details of the conformation, spanning ranges from ca. 1 ppm (for trimers of acids 1 and 2) to ca. 12 ppm (for trimers of 4). The energy barriers for the transitions between conformers are rather low (<6 kcal/mol for PTs, <2.5 kcal/mol for puckerings, and ca. <3 kcal/mol for rotations of substituents), such that the fast exchange regime in the NMR timescale and subsequent δP averaging are expected. Correlations are proposed linking the change of average δP with the H-bond energy, showing the slope of ca. 4 ppm per kcal/mol. The sensitivity of δP to the OPO angle and the OPOH dihedral angle and the geometries of both H-bonds formed by the POOH moiety are analyzed.

Solvent effects on acid-base complexes. What is more important: A macroscopic reaction field or solute-solvent interactions?

Can the geometry of an acid-base complex in solution be reproduced in calculations using an implicit accounting for the solvent effect in the form of a macroscopic reaction field? The answer is, "Yes, it can." Is this field equal to the real electric field experienced by the complex in solution? The answer is, "No, it is not." How can the geometry be correct under wrong conditions? This question is answered using density functional theory modeling of geometric and NMR parameters of pyridine⋯HF⋯(HCF₃)_n adducts in the absence and presence of an external electric field. This adduct under field approach shows that the N⋯H distance is a function of the H-F distance whatever method is used to change the geometry of the latter. An explicit account for solute-solvent interactions is required to get a realistic value of the solvent reaction field. Besides that, this approach reveals how certain NMR parameters depend on the solvent reaction field, the solute-solvent interactions, and the geometry of the N⋯H-F hydrogen bond. For some of them, the obtained dependences are far from self-evident.

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.

Influence of Hydrogen Bonds in 1:1 Complexes of Phosphinic Acids with Substituted Pyridines on 1H and 31P NMR Chemical Shifts

Two series of 1:1 complexes with strong OHN hydrogen bonds formed by dimethylphosphinic and phenylphosphinic acids with 10 substituted pyridines were studied experimentally by liquid state NMR spectroscopy at 100 K in solution in a low-freezing polar aprotic solvent mixture CDF₃/CDClF₂. The hydrogen bond geometries were estimated using previously established correlations linking ¹H NMR chemical shifts of bridging protons with the O···H and H···N interatomic distances. A new correlation is proposed allowing one to estimate the interatomic distance within the OHN bridge from the displacement of ³¹P NMR signal upon complexation. We show that the values of ³¹P NMR chemical shifts are affected by an additional CH···O hydrogen bond formed between the P═O group of the acid and ortho-CH proton of the substituted pyridines. Breaking of this bond in the case of 2,6-disubstituted bases shifts the ³¹P NMR signal by ca. 1.5 ppm to the high field.

Hydration and Hydroxylation of MgO in Solution: NMR Identification of Proton-Containing Intermediate Phases

Magnesium oxide (MgO)-engineered barriers used in subsurface applications will be exposed to high concentration brine environments and may form stable intermediate phases that can alter the effectiveness of the barrier. To explore the formation of these secondary intermediate phases, MgO was aged in water and three different brine solutions and characterized with X-ray diffraction (XRD) and ¹H magic angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy. After aging, there is ∼4% molar equivalent of a hydrogen-containing species formed. The ¹H MAS NMR spectra resolved multiple minor phases not visible in XRD, indicating that diverse disordered proton-containing environments are present in addition to crystalline Mg(OH)₂ brucite. Density functional theory (DFT) simulations for the proposed Mg-O-H-, Mg-Cl-O-H-, and Na-O-H-containing phases were performed to index resonances observed in the experimental ¹H MAS NMR spectra. Although the intermediate crystal structures exhibited overlapping ¹H NMR resonances in the spectra, Mg-O-H intermediates were attributed to the growth of resonances in the δ +1.0 to 0.0 ppm region, and Mg-Cl-O-H structures produced the increasing contributions of the δ = +2.5 to 5.0 ppm resonances in the chloride-containing brines. Overall, ¹H NMR analysis of aged MgO indicates the formation of a wide range of possible intermediate structures that cannot be observed or resolved in the XRD analysis.

Geometric H/D isotope effect in a series of organic salts involving short O–H⋯O hydrogen bonds between carboxyl and carboxylate groups

Noticeable elongations of donor–acceptor distances upon deuteration are confirmed in short O–H⋯O hydrogen bonds between carboxyl and carboxylate groups.

Refinement of labile hydrogen positions based on DFT calculations of 1H NMR chemical shifts: comparison with X-ray and neutron diffraction methods

Numerous gas phase electron diffraction, ultra-fast electron diffraction, X-ray and neutron diffraction experiments on β-dicarbonyl compounds exhibiting enol-enol tautomeric equilibrium, with emphasis on acetylacetone and dibenzoylmethane, have so far been reported with conflicting results on the structural details of the O-HO intramolecular hydrogen bond and resulted in alternative hypotheses on the intramolecular hydrogen bond potential function either a double minimum potential corresponding to two tautomeric forms in equilibrium or a single symmetrical one. We demonstrate herein, firstly, that the DFT calculated OH ¹H NMR chemical shifts of acetylacetone and dibenzoylmethane exhibit a strong linear dependence on the computed OO hydrogen bond length of ∼-50 ppm Å^-1 and as a function of the O-HO bond angle of ∼1 ppm per degree, upon the transfer of the hydrogen atom from the ground state toward the transition state. Secondly, the refinement of labile hydrogen atomic positions in intramolecular hydrogen bonds based on the root-mean-square deviation between experimentally determined and DFT calculated ¹H NMR chemical shifts in solution can provide high resolution structures of O-H and O(H)O bond lengths and O-HO bond angles with an accuracy of ∼10^-2 Å and ∼0.5°, respectively. Thirdly, the calculated ¹H NMR chemical shifts in solution of the two ground state tautomers in equilibrium of acetylacetone and dibenzoylmethane are in excellent agreement with the experimental value, even for moderate basis sets for energy minimization. In contrast, the single symmetrical structure in a strongly delocalized system is a transition state with calculated ¹H NMR chemical shifts which strongly deviate from the experimental value. Fourth, the DFT calculated ground state O-H bond lengths of acetylacetone and dibenzoylmethane are in quantitative agreement with the literature data which take into account the effect of quantum nuclear motion. The DFT structural results are critically discussed with respect to the state-of-the-art variable temperature X-ray and neutron diffraction methods.

Cyclic trimers of phosphinic acids in polar aprotic solvent: symmetry, chirality and H/D isotope effects on NMR chemical shifts

The hydrogen-bonded self-associates of dimethylphosphinic (1), diphenylphosphoric (2), phenylphosphinic (3), and bis(2,4,4-trimethylpentyl)phosphinic (4) acids have been studied by using liquid-state NMR down to 100 K in a low-freezing polar solvent, CDF₃/CDClF₂. The H/D isotope effects on ¹H NMR chemical shifts caused by partial deuteration of hydroxyl groups unambiguously reveal the stoichiometry of the self-associates and the cooperativity of their hydrogen bonds. In all cases, cyclic trimers are the dominant form, while cyclic dimers are present as a minor form for 1 and 2. Due to the asymmetry of substituents, cyclic trimers of 3 exist in two isomeric forms, depending on the orientation of the phenyl groups with respect to the plane of the hydrogen bonds. The racemic mixture of 4 leads to the coexistence of up to 64 isomers of cyclic trimers, many of which are chemically equivalent or effectively isochronous. The mole fractions of such isomers deviate from the statistically expected values. This feature could provide information about the relative stabilization energies of hydrogen-bonded chiral self-associates. The complexation of 4 with SbCl₅ (complex 5) suppresses the self-association and 5 exists exclusively in the monomeric form with chemically non-equivalent ³¹P nuclei in RS, SR and RR/SS forms.

Geometric isotope effect of deuteration in a hydrogen-bonded host-guest crystal

Deuteration of a hydrogen bond by replacing protium (H) with deuterium (D) can cause geometric changes in the hydrogen bond, known as the geometric H/D isotope effect (GIE). Understanding the GIEs on global structures and bulk properties is of great importance to study structure–property relationships of hydrogen-bonded systems. Here, we report a hydrogen-bonded host–guest crystal, imidazolium hydrogen terephthalate, that exemplifies striking GIEs on its hydrogen bonds, phases, and bulk dielectric transition property. Upon deuteration, the donor–acceptor distance in the O–H···O hydrogen bonds in the host structure is found to increase, which results in a change in the global hydrogen-bonded supramolecular structure and the emergence of a new phase (i.e., isotopic polymorphism). Consequently, the dynamics of the confined guest, which depend on the internal pressure exerted by the host framework, are substantially altered, showing a downward shift of the dielectric switching temperature.

Deuterating a hydrogen bond can change the bond’s geometry, a phenomenon known as the geometric isotope effect (GIE). Here, the authors find that a hydrogen-bonded host–guest crystal, imidazolium hydrogen terephthalate, exhibits significant GIE on its hydrogen bonds, changing its crystal phases and bulk dielectric properties.

Polar solvent fluctuations drive proton transfer in hydrogen bonded complexes of carboxylic acid with pyridines: NMR, IR and ab initio MD study

We study a series of intermolecular hydrogen-bonded 1 : 1 complexes formed by chloroacetic acid with 19 substituted pyridines and one aliphatic amine dissolved in CD₂Cl₂ at low temperature by ¹H and ¹³C NMR and FTIR spectroscopy. The hydrogen bond geometries in these complexes vary from molecular (O-HN) to zwitterionic (O^-H-N⁺) ones, while NMR spectra show the formation of short strong hydrogen bonds in intermediate cases. Analysis of C[double bond, length as m-dash]O stretching and asymmetric CO₂^- stretching bands in FTIR spectra reveal the presence of proton tautomerism. On the basis of these data, we construct the overall proton transfer pathway. In addition to that, we also study by use of ab initio molecular dynamics the complex formed by chloroacetic acid with 2-methylpyridine, surrounded by 71 CD₂Cl₂ molecules, revealing a dual-maximum distribution of hydrogen bond geometries in solution. The analysis of the calculated trajectory shows that the proton jumps between molecular and zwitterionic forms are indeed driven by dipole-dipole solvent-solute interactions, but the primary cause of the jumps is the formation/breaking of weak CHO bonds from solvent molecules to oxygen atoms of the carboxylate group.

NMR Study of Solvation Effect on the Geometry of Proton-Bound Homodimers of Increasing Size

Hydrogen bond geometries in the proton-bound homodimers of quinoline and acridine derivatives in an aprotic polar solution have been experimentally studied using ¹H NMR at 120 K. The reported results show that an increase of the dielectric permittivity of the medium results in contraction of the N···N distance. The degree of contraction depends on the homodimer's size and its substituent-specific solvation features. Neither of these effects can be reproduced using conventional implicit solvent models employed in computational studies. In general, the N···N distance in the homodimers of pyridine, quinoline, and acridine derivatives decreases in the sequence gas phase > solid state > polar solvent.

Quantum mechanics models of the methanol dimer: OH⋯O hydrogen bonds of β-d-glucose moieties from crystallographic data

The interaction of two methanol molecules, simplified models of carbohydrates and cellulose, was examined using a variety of quantum mechanics (QM) levels of theory. Energy plots for hydrogen bonding distance (H⋯O) and angle (OH⋯O) were constructed. All but two experimental structures were located in stabilized areas on the vacuum phase energy plots. Each of the 399 models was analyzed with Bader's atoms-in-molecules (AIM) theory, which showed a widespread ability by the dimer models to form OH⋯O hydrogen bonds that have bond paths and Bond Critical Points. Continuum solvation calculations suggest that a portion of the energy-stabilized structures could occur in the presence of water. A survey of the Cambridge Structural Database (CSD) for all donor-acceptor interactions in β-D-glucose moieties examined the similarities and differences among the hydroxyl groups and acetal oxygen atoms that participate in hydrogen bonds. Comparable behavior was observed for the O2H, O3H, O4H, and O6H hydroxyls, acting either as acceptors or donors. Ring O atoms showed distinct hydrogen bonding behavior that favored mid-length hydrogen bonds.

Hydrogen Atomic Positions of O-H···O Hydrogen Bonds in Solution and in the Solid State: The Synergy of Quantum Chemical Calculations with ¹H-NMR Chemical Shifts and X-ray Diffraction Methods

The exact knowledge of hydrogen atomic positions of O-H···O hydrogen bonds in solution and in the solid state has been a major challenge in structural and physical organic chemistry. The objective of this review article is to summarize recent developments in the refinement of labile hydrogen positions with the use of: (i) density functional theory (DFT) calculations after a structure has been determined by X-ray from single crystals or from powders; (ii) ¹H-NMR chemical shifts as constraints in DFT calculations, and (iii) use of root-mean-square deviation between experimentally determined and DFT calculated ¹H-NMR chemical shifts considering the great sensitivity of ¹H-NMR shielding to hydrogen bonding properties.

How High Concentrations of Proteins Stabilize the Amorphous State of Calcium Orthophosphate: A Solid-State Nuclear Magnetic Resonance (NMR) Study of the Casein Case

Understanding how proteins stabilize amorphous calcium ortho-phosphate (ACP) phases is of great importance in biology and for pharmaceutical or food applications. Until now, most of the former investigations about ACP-protein stability and equilibrium were performed under conditions where ACP colloidal nanoclusters are surrounded by low to moderate concentrations of peptides or proteins (15-30 g L^-1). As a result, the question of ACP-protein interactions in highly concentrated protein systems has clearly been overlooked, whereas it corresponds to actual industrial conditions such as drying or membrane filtration in the dairy industry for instance. In this study, the structure of an ACP phase is monitored in association with one model phosphorylated protein (casein) using solid-state nuclear magnetic resonance (ssNMR) under two conditions of high protein concentration (300 and 400 g L^-1). At both concentrations and at 25 °C, it is found that the caseins maintain the mineral phase in an amorphous form with no detectable influence on its structure or size. Interestingly, and in both cases, a significant amount of the nonphosphorylated side chains interacts with ACP through hydrogen bonds. The number of these interacting side chains is found to be higher at the highest casein concentration. At 45 °C, which is a destabilizing temperature of ACP under protein-free conditions, the amorphous structure of the mineral phase is partially transformed at a casein concentration of 300 g L^-1, while it remains almost intact at a casein concentration of 400 g L^-1. Therefore, these results clearly indicate that increasing the concentration of proteins favors ACP-protein interactions and stabilizes the ACP clusters more efficiently.