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

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.

Shared Hydrogen Atom Location and Chemical Composition in Picolinic Acid and Pyridoxal Hydrochloride Derivatives Determined by X-ray Crystallography

Three new single-crystal structures were isolated for picolinic acid (2), the trifluoroacetate salt of picolinic acid (1), and pyridoxal hydrochloride (3). These compounds displayed unconventional crystallographic features that must be considered when structural refinements are carried out. Thus, the generated Fourier differences map obtained with the diffraction data collected at 100 K was crucial to visualize electron densities, which were balanced by either one hydrogen atom or a hydrogen atom with an occupancy factor of 1/2 located between either two carboxylate moieties, two phenolic oxygen atoms, or two pyridinic nitrogen atoms. Moreover, NMR studies were conducted to analyze the bulk chemical composition of single crystals of 2-pyridinecarboxylic acid obtained from the gem-diol/hemiacetal forms and the polymerization products after the treatment of 2-pyridinecarboxaldehyde with TFA:H₂O (1) or a diluted Cu(NO₃)₂ solution (2). The quantitative yield of the pyridoxal hydrochloride crystalline material (3) obtained from a diluted CuCl₂ solution was exhaustively characterized by solid-state NMR methods. These methods allowed the resolution of the signals corresponding to the protons of the hydroxyl moiety of the intramolecular hemiacetal group and the phenolic hydrogen. Theoretical calculations using DFT methods were done to complement the atomic location of the hydrogen atoms obtained from the X-ray analysis.

Step-by-step from amorphous phosphate to nano-structured calcium hydroxyapatite: monitoring by solid-state 1H and 31P NMR and spin dynamics

The solid-state ¹H, ³¹P NMR spectra and cross-polarization (CP MAS) kinetics in the series of samples containing amorphous phosphate phase (AMP), composite of AMP + nano-structured calcium hydroxyapatite (nano-CaHA) and high-crystalline nano-CaHA were studied under moderate spinning rates (5-30 kHz). The combined analysis of the solid-state ¹H and ³¹P NMR spectra provides the possibility to determine the hydration numbers of the components and the phase composition index. A broad set of spin dynamics models (isotropic/anisotropic, relaxing/non-relaxing, secular/semi-non-secular) was applied and fitted to the experimental CP MAS data. The anisotropic model with the angular averaging of dipolar coupling was applied for AMP and nano-CaHA for the first time. It was deduced that the spin diffusion in AMP is close to isotropic, whereas it is highly anisotropic in nano-CaHA being close to the Ising-type. This can be caused by the different number of internuclear interactions that must be explicitly considered in the spin system for AMP (I-S spin pair) and nano-CaHA (I_N-S spin system with N ≥ 2). The P-H distance in nano-CaHA was found to be significantly shorter than its crystallographic value. An underestimation can be caused by several factors, among those - proton conductivity via a large-amplitude motion of protons (O-H tumbling and the short-range diffusion) that occurs along OH^- chains. The P-H distance deduced for AMP, i.e. the compound with HPO₄^2- as the dominant structure, is fairly well matched to the crystallographic data. This means that the CP MAS kinetics is a capable technique to obtain complementary information on the proton localization in H-bonds and the proton transfer in the cases where traditional structure determination methods fail.

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.

NMR Properties of the Cyanide Anion, a Quasisymmetric Two-Faced Hydrogen Bonding Acceptor

The isotopically enriched cyanide anion, (13C≡15N)−, has a great potential as the NMR probe of non-covalent interactions. However, hydrogen cyanide is highly toxic and can decompose explosively. It is therefore desirable to be able to theoretically estimate any valuable results of certain experiments in advance in order to carry out experimental studies only for the most suitable molecular systems. We report the effect of hydrogen bonding on NMR properties of 15N≡13CH···X and 13C≡15NH···X hydrogen bonding complexes in solution, where X = 19F, 15N, and O=31P, calculated at the ωB97XD/def2tzvp and the polarizable continuum model (PCM) approximations. In many cases, the isotropic 13C and 15N chemical shieldings of the cyanide anion are not the most informative NMR properties of such complexes. Instead, the anisotropy of these chemical shieldings and the values of scalar coupling constants, including those across hydrogen bonds, can be used to characterize the geometry of such complexes in solids and solutions. 1J(15N13C) strongly correlates with the length of the N≡C bond.

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

Towards a new pseudo-quantitative approach to evaluate the ionization response of nitrogen compounds in complex matrices

Ultra high-resolution mass spectrometry (FT-ICR MS) coupled to electrospray ionization (ESI) provides unprecedented molecular characterization of complex matrices such as petroleum products. However, ESI faces major ionization competition phenomena that prevent the absolute quantification of the compounds of interest. On the other hand, comprehensive two-dimensional gas chromatography (GC × GC) coupled to specific detectors (HRMS or NCD) is able to quantify the main families identified in these complex matrices. In this paper, this innovative dual approach has been used to evaluate the ionization response of nitrogen compounds in gas oils as a case study. To this extent, a large gas oil dataset has been analyzed by GC × GC/HRMS, GC × GC-NCD and ESI(+/−)-FT-ICR MS. Then, the concentrations obtained from GC × GC-NCD have been compared to those obtained from FT-ICR MS hence proving that strong ionization competitions are taking place and also depending on the origin of the sample. Finally, multilinear regressions (MLR) have been used to quantitatively predict nitrogen families from FT-ICR MS measurements as well as start rationalizing the ionization competition phenomena taking place between them in different types of gas oils.

1,3,5-Triaza-7-Phosphaadamantane (PTA) as a 31P NMR Probe for Organometallic Transition Metal Complexes in Solution

Due to the rigid structure of 1,3,5-triaza-7-phosphaadamantane (PTA), its ³¹P chemical shift solely depends on non-covalent interactions in which the molecule is involved. The maximum range of change caused by the most common of these, hydrogen bonding, is only 6 ppm, because the active site is one of the PTA nitrogen atoms. In contrast, when the PTA phosphorus atom is coordinated to a metal, the range of change exceeds 100 ppm. This feature can be used to support or reject specific structural models of organometallic transition metal complexes in solution by comparing the experimental and Density Functional Theory (DFT) calculated values of this ³¹P chemical shift. This approach has been tested on a variety of the metals of groups 8-12 and molecular structures. General recommendations for appropriate basis sets are reported.

Modeling of Solute-Solvent Interactions Using an External Electric Field-From Tautomeric Equilibrium in Nonpolar Solvents to the Dissociation of Alkali Metal Halides

An implicit account of the solvent effect can be carried out using traditional static quantum chemistry calculations by applying an external electric field to the studied molecular system. This approach allows one to distinguish between the effects of the macroscopic reaction field of the solvent and specific solute-solvent interactions. In this study, we report on the dependence of the simulation results on the use of the polarizable continuum approximation and on the importance of the solvent effect in nonpolar solvents. The latter was demonstrated using experimental data on tautomeric equilibria between the pyridone and hydroxypyridine forms of 2,6-di-tert-butyl-4-hydroxy-pyridine in cyclohexane and chloroform.

Electric field effect on 31P NMR magnetic shielding

Magnetic shielding depends on molecular structure and noncovalent interactions. This study shows that it is also measurably dependent on the electric field generated by surrounding molecules. This effect has been observed explicitly for ³¹P nucleus using the adduct under field approach. The results obtained indicate that the field strength experienced by molecules in crystals consisting of molecules with large dipole moments is similar to that in polar solvents. Therefore, magnetic shielding should explicitly depend on solvent polarity. It is important to note that this effect cannot be reproduced correctly within the polarizable continuum model approach.

Inter- vs. Intramolecular Hydrogen Bond Patterns and Proton Dynamics in Nitrophthalic Acid Associates

Noncovalent interactions are among the main tools of molecular engineering. Rational molecular design requires knowledge about a result of interplay between given structural moieties within a given phase state. We herein report a study of intra- and intermolecular interactions of 3-nitrophthalic and 4-nitrophthalic acids in the gas, liquid, and solid phases. A combination of the Infrared, Raman, Nuclear Magnetic Resonance, and Incoherent Inelastic Neutron Scattering spectroscopies and the Car-Parrinello Molecular Dynamics and Density Functional Theory calculations was used. This integrated approach made it possible to assess the balance of repulsive and attractive intramolecular interactions between adjacent carboxyl groups as well as to study the dependence of this balance on steric confinement and the effect of this balance on intermolecular interactions of the carboxyl groups.

For Whom a Puddle Is the Sea? Adsorption of Organic Guests on Hydrated MCM-41 Silica

Thermal and hydration effects on the mobility of compact and branched organic molecules and a bulky pharmaceutical substance loaded in submonolayer amounts onto mesoporous silica have been elucidated using ¹H and ³¹P solid-state NMR. In all cases, the ambient hydration has a stronger effect than an increase in temperature to 370 K for water-free silica. The effect of hydration depends on the guest and ranges from complete solvation to a silica-water-guest sandwich structure to a silica-guest/silica-water pattern. The mobility of the guests under different conditions has been described. The specific structure of the MCM-41 surface allows one to study very slow surface diffusion, a diffusivity of about 10^-15-10^-16 m²/s. The data reported are relevant to any nonfunctionalized silica, while the method used is applicable to any phosphor-containing guest on any host.

Gas Phase Computational Study of Diclofenac Adsorption on Chitosan Materials

Environmental pollution with non-steroidal anti-inflammatory drugs and their metabolites exposes living organisms on their long-lasting, damaging influence. Hence, the ways of non-steroidal anti-inflammatory drugs (NSAIDs) removal from soils and wastewater is sought for. Among the potential adsorbents, biopolymers are employed for their good availability, biodegradability and low costs. The first available theoretical modeling study of the interactions of diclofenac with models of pristine chitosan and its modified chains is presented here. Supermolecular interaction energy in chitosan:drug complexes is compared with the the mutual attraction of the chitosan dimers. Supermolecular interaction energy for the chitosan-diclofenac complexes is significantly lower than the mutual interaction between two chitosan chains, suggesting that the diclofenac molecule will encounter problems when penetrating into the chitosan material. However, its surface adsorption is feasible due to a large number of hydrogen bond donors and acceptors both in biopolymer and in diclofenac. Modification of chitosan material introducing long-distanced amino groups significantly influences the intramolecular interactions within a single polymer chain, thus blocking the access of diclofenac to the biopolymer backbone. The strongest attraction between two chitosan chains with two long-distanced amino groups can exceed 120 kcal/mol, while the modified chitosan:diclofenac interaction remains of the order of 20 to 40 kcal/mol.

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.

Hydrogen peroxide adducts of triarylphosphine oxides

Five new hydrogen peroxide adducts of phosphine oxides (p-Tol₃PO·H₂O₂)₂ (1), (o-Tol₃PO·H₂O₂)₂ (2), (o-Tol₂PhPO·H₂O₂)₂ (3), (p-Tol₃PO)₂·H₂O₂ (4), and (o-TolPh₂PO)₂·H₂O₂ (5), and the water adduct (o-Tol₂PhPO·H₂O)₂ (6) have been synthesized and fully characterized. Their single crystal X-ray structures have been determined and analyzed. The IR and ³¹P NMR data are in accordance with strong hydrogen bonding of the hydrogen peroxide. The mono- versus dimeric nature of the adduct assemblies has been investigated by DOSY NMR experiments. Raman spectroscopy of the symmetric adducts and the ν(O-O) stretching bands confirm the presence of hydrogen-bonded hydrogen peroxide in the solid materials. The solubilities in organic solvents have been quantified. Due to the high solubilities of 1-6 in organic solvents their ¹⁷O NMR spectra could be recorded in natural abundance, providing well-resolved signals for the P[double bond, length as m-dash]O and O-O groups. The adducts 1-5 have been probed regarding their stability in solution at 105 °C. The decomposition of the adduct 1 takes place by loss of the active oxygen atoms in two steps.

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.

Simplified calculation approaches designed to reproduce the geometry of hydrogen bonds in molecular complexes in aprotic solvents

The impact of the environment onto the geometry of hydrogen bonds can be critically important for the properties of the questioned molecular system. The paper reports on the design of calculation approaches capable to simulate the effect of aprotic polar solvents on the geometric and NMR parameters of intermolecular hydrogen bonds. A hydrogen fluoride and pyridine complex has been used as the main model system because the experimental estimates of these parameters are available for it. Specifically, F-H, F⋯N, and H-N distances, the values of ¹⁵N NMR shift, and spin-spin coupling constants ¹J(¹⁹F¹H), ^1hJ(¹H¹⁵N), and ^2hJ(¹⁹F¹⁵N) have been analyzed. Calculation approaches based on the gas-phase and the Polarizable Continuum Model (PCM) approximations and their combinations with geometric constraints and additional noncovalent interactions have been probed. The main result of this work is that the effect of an aprotic polar solvent on the geometry of a proton-donor⋯H⋯proton-acceptor complex cannot be reproduced under the PCM approximation if no correction for solvent-solute interactions is made. These interactions can be implicitly accounted for using a simple computational protocol.