Correlations and predictions of carboxylic acid pKa values using intermolecular structure and properties of hydrogen-bonded complexes

Density functional theory calculations have been preformed on a series of hydrogen-bonded complexes of substituted aliphatic and aromatic carboxylic acids with ammonia. Molecular properties, particularly those related to hydrogen bonding, have been carefully examined for their interdependence as well as dependence on the acidity of the acid. The bond length and stretching frequency of the hydroxyl group and the hydrogen-bond length and energy of the complex are shown to be highly correlated with each other and are linearly correlated with available literature pKa values of the carboxylic acids. The linear correlations resulting from the fit to the available pKa values can be used to predict the pKa values of similar carboxylic acids. The pKa values so predicted using the different molecular properties are highly consistent and in good agreement with the literature values. This study suggests that calculated molecular properties of hydrogen-bonded complexes allow effective and systematic prediction of pKa values for a large range of organic acids using the established linear correlations. This approach is unique in its capability to determine the acidity of a particular functional group or the local acidity within a large molecular system such as a protein.     

Theoretical study of hydrogen-bonded complexes of chlorophenols with water or ammonia: Correlations and predictions of pKa values

A simple practical method for predicting the acidity constants (as pKa values) of chlorophenols is proposed based on density functional theory calculations of a series of hydrogen-bonded complexes of phenol and 19 different congeners of chlorophenol, with a single probe molecule, either water or ammonia. Relevant structural parameters and molecular properties of these complexes, primarily involving the acidic hydroxyl group, are examined and plotted against the known pKa values of 14 chlorophenols and phenol. Strong linear correlations are found for these compounds. Such correlations are used to determine the pKa values of five chlorophenols whose experimental acidities have large uncertainties. Similar predicted pKa values are obtained by using different structural parameters and molecular properties for the complexes with either probe molecule. The study may be extended to determine the acidity of other compounds with a single acidic functional group.     

Molecular structures and properties of the complete series of bromophenols: Density functional theory calculations

The complete series of 19 bromophenols have been studied by density functional theory (DFT) calculations at the B3LYP/6-311G++(d,p) level. The molecular structures and properties of bromophenols are strongly influenced by intramolecular hydrogen bonding of ortho-bromine, steric and inductive effects of substituted bromine, and other intramolecular electrostatic interactions. Systematic trends in several structural parameters and molecular properties of bromophenols have been found with the increasing number of bromine substitutions, including increase in O-H bond length, decrease in C-O bond length, red shift in O-H stretching frequency, and blue shift in O-H torsional frequency. Correlations among several key molecular parameters as well as those with available aqueous pKa values are examined. Comparisons with chlorophenols have indicated that the inductive effect of substituted bromine appears larger and bromophenols are slightly stronger acids than chlorophenols.     

NMR parameters and geometries of OHN and ODN hydrogen bonds of pyridine-acid complexes

In this paper, equations are proposed which relate various NMR parameters of OHN hydrogen-bonded pyridine-acid complexes to their bond valences which are in turn correlated with their hydrogen-bond geometries. As the valence bond model is strictly valid only for weak hydrogen bonds appropriate empirical correction factors are proposed which take into account anharmonic zero-point energy vibrations. The correction factors are different for OHN and ODN hydrogen bonds and depend on whether a double or a single well potential is realized in the strong hydrogen-bond regime. One correction factor was determined from the known experimental structure of a very strong OHN hydrogen bond between pentachlorophenol and 4-methylpyridine, determined by the neutron diffraction method. The remaining correction factors which allow one also to describe H/D isotope effects on the NMR parameters and geometries of OHN hydrogen bond were determined by analysing the NMR parameters of the series of protonated and deuterated pyridine- and collidine-acid complexes. The method may be used in the future to establish hydrogen-bond geometries in biologically relevant functional OHN hydrogen bonds.     

Varying acidity of aqua ligands in dependence on the microenvironment in mononucleobase (nb) complexes of type cis- and trans-[Pt(NH3)2(nb)(H2O)]n+

Aqua ligands in mixed aqua/nucleobase metal complexes are potential sites of acid-base catalysis and/or, when present as hydroxo ligands, can directly be involved in hydrolysis reactions. pKa values of close to 7 are consequently of particular interest and potential significance. Here we report on the differential acidity of aqua complexes in model nucleobase (nb) complexes of cis- and trans-[Pt(NH3)2 (nb)(H2O)]n+ and discuss reasons as to why the nb in cis complexes influences the pKa (pKa 4.8-7.0), whereas in trans complexes the pKa values are rather constant (pKa approximately 5.2-5.3). The results of DFT calculations of a series of mono(nucleobase) complexes derived from cis-Pt(NH3)2 are critically examined with regard to the role of exocyclic groups of nucleobases in stabilizing aqua/hydroxo ligands through intracomplex hydrogen bond formation. This applies in particular to the exocyclic amino groups of nucleobases, for which gas-phase calculations suggest that they may act as H bond acceptors in certain cases, yet in the condensed phase this appears not to be the case.     

Quantitative relationships between bond lengths,stretching vibrational frequencies,bond force constants,and bond orders in the hydrogen-bonded complexes involving hydrogen halides

To uncover the correlation between the bond length change and the corresponding stretching frequency shift of the proton donor D–H upon hydrogen bond formation, a series of hydrogen-bonded complexes involving HF and HCl which exhibit the characteristics of red-shifted hydrogen bond were investigated at the MP2/aug-cc-pVTZ, M062X/aug-cc-pVTZ, and B3LYP/aug-cc-pVTZ(GD3) levels of theory with CP optimizations. A statistical analysis of these complexes leads to the quantitative illustrations of the relations between bond length and stretching vibrational frequency, between bond length and bond force constant, between stretching vibrational frequency and bond force constant, between bond length and bond order for hydrohalides in a mathematical way, which would provide valuable insights into the explanation of the geometrical and spectroscopic behaviors during hydrogen bond formation.     

Gas sorption and diffusion in hydrogen-bonded polymers. I. Vinyl alcohol/vinyl butyral copolymers

A series of vinyl alcohol/vinyl butyral copolymers was examined to assess the effect of internal hydrogen bonding on gas sorption and diffusion. Sorption and permeation measurements for carbon dioxide and methane were performed on four vinyl alcohol/vinyl butyral copolymers. Upon comparing the various data, it was found that hydrogen-bonded copolymers exhibit a much wider variation in diffusion coefficient than non-hydrogen-bonded copolymers. The fractional free volumes of the studied copolymers were considerably lower than expected based on values of the diffusion coefficient. This may be due to the fact that predicted occupied volumes are too large and that the effect of internal hydrogen bonding is not accounted for properly. Using a relationship between infrared spectral shift and hydrogen-bond length, fractional free volumes in the hydrogen bonded copolymers were correlated with the interatomic spacing associated with the hydrogen bond. This implies that the average length of a hydrogen bond can be used as a measure of chain packing in hydrogen-bonded polymers. © 1993 John Wiley & Sons, Inc.     

The effect of hydrogen bonding on structural parameters ii. an ab initio study of the monohydrated formamide molecule

Four hydrogen-bonded formamide-water complexes have been studied by ab initio calculations, two where the amino group acts as a donor and two where the carbonyl oxygen is an acceptor. The results indicate that the effect on the conjugated NCO fragment depends on both the type and the energy of the hydrogen bond formed. Although, in all cases the formation of a hydrogen bond leads to increased conjugation, expressed as a shortening of the CN bond and a corresponding lengthening of the CO bond, there is a significant difference in the effect of the two types of hydrogen bonds. This difference may be explained by changes in the electron populations. In two of the complexes the effect of varying the hydrogen bond length has been studied in some detail. It is found that the effect on the conjugated system depends on the length of the hydrogen bond, and analytical expressions have been found for the variations of the CO and CN bond lengths with changes in the hydrogen bond length. Potential functions for the N-H β O and O-H β O hydrogen bonds have also been derived.     

Hydrogen-bonded complexes of serotonin with methanol and ethanol: a DFT study

Density functional theoretical studies on hydrogen-bonded complexes of serotonin with methanol/ethanol have been carried out in a systematic way. The conformational analysis led to ten stable conformers that can be either gauche or anti depending on the dihedral angle values taken by ethylamine side chain and the 5-hydroxyl group. Serotonin-molecules strongly bind with ethanol than methanol. Ethylamine side chain is the most reactive site in both methanol/ethanol complexes and it is responsible for the stability order. The topological parameters, electron density, and Laplacian of electron density show excellent correlation with the hydrogen bond length. Natural bond orbital analysis confirms C–H···O hydrogen bond formed between the serotonin–alcohol complexes to be red shift in nature except for Gph(out)anti complex both with methanol and ethanol to be blue shifted. The energy decomposition analysis reveals that strong interactions between serotonin and ethanol/methanol are due to the attractive contributions from the electrostatic component.     

The effect of hydrogen bonding on structural parameters : IV. An ab initio study of formamide—ammonia complexes

Four hydrogen-bonded formamide—ammonia complexes have been studied by ab initio calculations, two where the amino group acts as a donor and two where the carbonyl oxygen is acceptor. In all cases the formation of a hydrogen bond leads to increased conjugation, expressed as a shortening of the CN bond and a corresponding lengthening of the CO bond. In two of the complexes the effect of varying the hydrogen bond length has been studied in some detail. It is found that the effect on the conjugated system depends on the length of the hydrogen bond. Potential functions for the N—HN and N—HO hydrogen bonds have been derived.     

Study of hydrogen-bonded clusters of 2-methoxyphenol–water

Various hydrogen-bonded clusters of 2-methoxyphenol (2MP) with water have been analyzed using ab initio methods and Atoms in Molecules (AIM) theory. The intramolecular hydrogen bond energy (and enthalpy) for 2MP was evaluated from two different methods. The results of rotational barriers method are in better agreement with experimental data. Binding energies, vibrational frequencies and geometrical parameters were examined and compared for these complexes. It was shown that in the most stable complex, water acts both as a donor and an acceptor. The “bifurcated” complex was shown to be relatively stable based on energy values. Atoms in Molecules and Natural Bond Orbital (NBO) analysis were used to confirm the existence of hydrogen bonds and to compare the strengths of them. The results obtained from quantum mechanical, AIM and NBO calculations are in agreement with each other.     

Quantification of hyperconjugative effect on the proton donor X-H bond length changes in the red- and blueshifted hydrogen-bonded complexes

A whole dataset containing 55 hydrogen bonds were studied at the MP2∕aug-cc-pVTZ level of theory. The changes of geometries and stretching vibrational frequencies show that there are 31 redshifted and 24 blueshifted hydrogen-bonded complexes. Natural bond orbital analysis was carried out at the B3LYP∕aug-cc-pVTZ level of theory to obtain the electron densities in the bonding and antibonding orbitals of the proton donor X-H bond, which are closely associated with its bond length. Based on their relationship, a generally applicable method considering both the electron densities in the bonding and antibonding orbitals of the proton donor X-H bond has been developed to quantitatively describe the hyperconjugative effect on the X-H bond length changes in these hydrogen-bonded complexes.     

Ultrafast intermolecular hydrogen bond dynamics in the excited state of fluorenone

Steady-state fluorescence and time-resolved absorption measurements in pico- and femtosecond time domain have been used to investigate the dynamics of hydrogen bond in the excited singlet (S(1)) state of fluorenone in alcoholic solvents. A comparison of the features of the steady-state fluorescence spectra of fluorenone in various kinds of media demonstrates that two spectroscopically distinct forms of fluorenone in the S(1) state, namely the non-hydrogen-bonded (or free) molecule as well as the hydrogen-bonded complex, are responsible for the dual-fluorescence behavior of fluorenone in solutions of normal alcoholic solvents at room temperature (298 K). However, in 2,2,2-trifluoroethanol (TFE), a strong hydrogen bond donating solvent, emission from only the hydrogen-bonded complex is observed. Significant differences have also been observed in the temporal evolution of the absorption spectroscopic properties of the S(1) state of fluorenone in protic and aprotic solvents following photoexcitation using 400 nm laser pulses. An ultrafast component representing the solvent-induced vibrational energy relaxation (VER) process has been associated with the dynamics of the S(1) state of fluorenone in all kinds of solvents. However, in protic solvents, in addition to the VER process, further evolution of the spectroscopic and dynamical properties of the S(1) state have been observed because of repositioning of the hydrogen bonds around the carbonyl group. In normal alcohols, two different kinds of hydrogen-bonded complex of the fluorenone-alcohol system with different orientations of the hydrogen bond with respect to the carbonyl group and the molecular plane of fluorenone have been predicted. On the other hand, in TFE, formation of only one kind of hydrogen-bonded complex has been observed. These observations have been supported by theoretical calculations of the geometries of the hydrogen-bonded complexes in the ground and the excited states of fluorenone. Linear correlation between the lifetimes of the equilibration process occurring because of repositioning of the hydrogen bonds and Debye or longitudinal relaxation times of the normal alcoholic solvents establish the fact that, in weakly hydrogen bond donating solvents, the hydrogen bond dynamics can be described as merely a solvation process. Whereas, in TFE, hydrogen bond dynamics is better described by a process of conversion between two distinct excited states, namely, the non-hydrogen-bonded form and the hydrogen-bonded complex.     

Vibrational spectroscopic properties of hydrogen bonded acetonitrile studied by DFT

Vibrational properties (band position, Infrared and Raman intensities) of the acetonitrile C[triple bond]N stretching mode were studied in 27 gas-phase medium intensity (length range: = 1.71-2.05 angstroms; -deltaE range = 13-48 kJ/mol) hydrogen-bonded 1:1 complexes of CH3CN with organic and inorganic acids using density functional theory (DFT) calculations [B3LYP-6-31++G(2d,2p)]. Furthermore, general characteristics of the hydrogen bonds and vibrational changes in the OH stretching band of the acids were also considered. Experimentally observed blue-shifts of the C[triple bond]N stretching band promoted by the hydrogen bonding, which shortens the triple bond length, are very well reproduced and quantitatively depend on the hydrogen bond length. Both predicted enhancement of the infrared and Raman nu(C[triple bond]N) band intensities are in good agreement with the experimental results. Infrared band intensity increase is a direct function of the hydrogen bond energy. However, the predicted increase in the Raman band intensity increase is a more complex function, depending simultaneously on the characteristics of both the hydrogen bond (C[triple bond]N bond length) and the H-donating acid polarizability. Accounting for these two parameters, the calculated nu(C[triple bond]N) Raman intensities of the complexes are explained with a mean error of +/- 2.4%.     

T-型氢键液晶复合物的制备与液晶行为

以不具有液晶行为的2,6-二[N,N′-二-(4-烷基苯甲酰基)]氨基吡啶(A系列)和4-正烷氧基苯甲酸(D系列)作为氢键液晶复合物的单体,组装成T-型氢键液晶系列复合物(AmDn)。用红外光谱对其结构进行了表征,用DSC及偏光显微方法对其液晶行为进行了研究。结果表明:所合成的21种复合物分子间存在氢键且都具有向列相。通过调整2,6-二[N,N′-二-(4-烷基苯甲酰基)]氨基吡啶分子上柔性烷基的长度和极性,可以有效地调节它与4-烷氧基苯甲酸形成的氢键复合物的液晶相变温度以及液晶态的稳定性;增加2,6-二[N,N′-二-(4-烷基苯甲酰基)]氨基吡啶分子上柔性烷基的长度,其复合物AmDn的液晶相温度范围趋于变窄,清亮点逐渐下降,其液晶态稳定性也逐渐下降;以2,6-二[N,N′-二-(4-烷基苯甲酰基)]氨基吡啶分子替代2,6-二[N,N′-二-(4-烷氧基苯甲酰基)]氨基吡啶分子,可以降低分子的极性,使其单体的熔点及其氢键复合物AmDn的相变温度下降。     

Electron density and energy decomposition analysis in hydrogen-bonded complexes of azabenzenes with water, acetamide, and thioacetamide

Ab initio and density functional theoretical studies on hydrogen-bonded complexes of azabenzenes with water, acetamide, and thioacetamide have been carried out to explore the controversy involved in the relative order of their stability in a systematic way. The interaction energies of these complexes have been analyzed using the Morokuma energy decomposition method, and the nature of the various hydrogen bonds formed has been investigated through topological aspects using Bader's atom in a molecule (AIM) theory. Morokuma energy decomposition analysis reveals that the major contributions to the energetics are from the polarization (PL) and charge transfer (CT) energies. From the calculated topological results, excellent linear correlation is shown to exist between the hydrogen-bond length, electron density [rho(r)], and its Laplacian [nabla(2)rho(r)] at the bond critical points for all the complexes considered.     

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.     

Vibrational spectral diffusion and hydrogen bond dynamics in heavy water from first principles

We present a first-principles theoretical study of vibrational spectral diffusion and hydrogen bond dynamics in heavy water without using any empirical model potentials. The calculations are based on ab initio molecular dynamics simulations for trajectory generation and a time series analysis using the wavelet method for frequency calculations. It is found that, in deuterated water, although a one-to-one relation does not exist between the instantaneous frequency of an OD bond and the distance of its associated hydrogen bond, such a relation does hold on average. The dynamics of spectral diffusion is investigated by means of frequency-time correlation and spectral hole dynamics calculations. Both of these functions are found to have a short-time decay with a time scale of approximately 100 fs corresponding to dynamics of intact hydrogen bonds and a slower long-time decay with a time constant of approximately 2 ps corresponding to lifetimes of hydrogen bonds. The connection of the slower time scale to the dynamics of local structural relaxation is also discussed. The dynamics of hydrogen bond making is shown to have a rather fast time scale of approximately 100 fs; hence, it can also contribute to the short-time dynamics of spectral diffusion. A damped oscillation is also found at around 150-200 fs, which is shown to have come from underdamped intermolecular vibrations of a hydrogen-bonded water pair. Such assignments are confirmed by independent calculations of power spectra of intermolecular motion and hydrogen bond kinetics using the population correlation function formalism. The details of the time constants of frequency correlations and spectral shifts are found to depend on the frequencies of chosen OD bonds and are analyzed in terms of the dynamics of hydrogen bonds of varying strengths and also of free non-hydrogen-bonded OD groups.     

Theoretical investigation of in-plane hydrogen-bonded complexes of ammonia with partially substituted fluorobenzenes

Systematic investigation of in-plane hydrogen-bonded complexes of ammonia with partially substituted fluorobenzenes has revealed that fluorobenzene, difluorobenzene, and trifluorobenzene favor formation of cyclic complexes with a C-H...N-H...F-C binding motif. On the other hand, tetrafluorobenzene and pentafluorobenzene favor formation of linear C-H...N hydrogen-bonded complexes. The complete absence of exclusively linear N-H...F hydrogen-bonded complexes for the entire series indicates that C-F bond in fluorobenzenes is a reluctant hydrogen-bond acceptor. However, fluorine does hydrogen bond when cooperatively stabilized with C-H...N hydrogen bonds for the lower fluoro analogues. The propensity of fluorobenzenes to adapt to the C-H...N-H...F-C binding motif decreases with the progressive fluorination of the benzene ring and disappears completely when benzene ring is substituted with five or more fluorine atoms.