Irena Majerz*a and
Matthias J. Gutmannb
aFaculty of Pharmacy, Wroclaw Medical University, Borowska 211a, 50-556 Wroclaw, Poland. E-mail: majerz@yahoo.com; irena.majerz@umed.wroc.pl
bRutherford Appleton Laboratory, ISIS Facility, Chilton Didcot, Oxfordshire OX11 0QX, UK
First published on 2nd November 2015
The neutron structure of 3-methylpyridinium 2,6-dichloro-4-nitrophenolate measured at different temperatures is used to discuss the structural parameters of the OHN intermolecular hydrogen bridges with the proton moving. The experimental and optimized structures, with the precise location of the proton, are used to explain the mechanism of proton transfer in the OHN hydrogen bridge.
Among the intermolecular hydrogen bonds, the most typical are the complexes of phenols with tertiary amines. They form single hydrogen bonds without chains or bifurcated hydrogen bonds with strengths covering the whole range of hydrogen bonds, from molecular to very strong bonds with the proton at the center of the hydrogen bridge, up to ionic complexes with the proton shifted to the acceptor nitrogen atom. The complexes of phenols with tertiary amines are used as model complexes to investigate the interatomic OHN hydrogen bond. Throughout the long period of their investigation, many correlations linking the degree of proton transfer with structural,6 physicochemical,7 and spectroscopic parameters8 have been published. The most essential step in investigating the hydrogen bond is to determine its structure, with the precise location of the proton engaged in the hydrogen bond.
The structural parameters which can characterize the degree of proton transfer in the OHN hydrogen bond in phenol–tertiary amine complexes are the bond lengths connected with the hydrogen bridge, i.e. O⋯N, OH, and NH. The most characteristic is the O⋯N distance. For the strongest hydrogen bonds it is very short, below 2.6 Å. The proton in strong hydrogen bonds is located in the middle of the bridge, so the OH and NH distances are similar. Elongation of O⋯N is connected with a shifting of the proton to the donor or to the acceptor molecule, according to the properties of the donor and acceptor in the hydrogen bond complex. The OHN angle exhibits a general tendency to be linear when the O⋯N bridge is shortest. Some distances of the phenol molecule are also sensitive to the proton transfer, so the C–O distance can be used as a measure of the degree of proton transfer.
The mutual dependencies of O⋯N with OH and NH shown in9 have parabolic shape, and the correlation between OH and NH is sigmoidal. All these distances can be correlated with the value ΔpKa (ΔpKa = pKb base − pKa acid).
The dependencies of the structural parameters of the hydrogen bridge obtained experimentally for phenol–amine complexes have been compared with these parameters predicted using theoretical methods.10 Experimental structural parameters follow those calculated at the B3LYP/6-311++G** level. In this situation, one may expect that on the basis of the ΔpKa value of the phenol–amine complex the degree of proton transfer in the OHN bond can be predicted as well as all the structural parameters sensitive to the location of the proton in the hydrogen bridge.
Among numerous complexes with intermolecular hydrogen bond only few are characterized by very short distance between the proton donor and proton acceptor. Neutron structures of these complexes with the O⋯N distances shorter than 2.6 Å measured at different temperatures confirmed possibility of the proton moving along the hydrogen bridge upon the temperature.11,12
In this paper, the neutron structure of the of 3-methylpyridinium 2,6-dichloro-4-nitrophenolate that is a typical phenol–pyridine complex is shown. The ΔpKa value of this complex equals 1.54, suggesting a short O⋯N distance and with possible symmetric location of the proton in the hydrogen bridge. Previously investigated X-ray structure showed that the parameters of the OHN hydrogen bond (O⋯ − 2.544(4), OH – 1.60(3), NH – 0.97(3), CO – 1.285(4) Å, OHN – 165(1)° (ref. 13)) are typical for very short hydrogen bond in which the proton can move along the hydrogen bridge upon the temperature change. The location of the hydrogen bonded proton by X-ray diffraction has very low precision, and therefore to investigate the details of the proton moving the neutron structure of 3-methylpyridinium 2,6-dichloro-4-nitrophenolate has been measured. The aim of this study is a precise determination of the structural parameters of the hydrogen bond in the investigated complex at different temperature. The short hydrogen bond has been characterized by calculation of the potential energy curves for the proton moving in the hydrogen bridge. Combination of the experimental neutron structures and theoretical structures optimized for different location of the proton in the hydrogen bridge gives possibility to follow the proton transfer mechanism investigated by NBO analysis.
Neutron diffraction data in the temperature range from 10 K to 290 K were collected using the SXD instrument at the ISIS facility at Rutherford Appleton Laboratory (UK).14 To facilitate efficient data collection, two crystals were comounted and exposed simultaneously in the neutron beam, similar as was done in.15 Data were treated using the locally available SXD2001 software16 and refinements carried out using Jana2006.17 Details of the data collection and refinement parameters are summarized in Table 1.
Diffractometer | SXD neutron time-of-flight Laue diffractometer | ||||||||||
Wavelength (Å) | 0.37–8.8 | ||||||||||
Compound | 3-Methylpyridinium 2,6-dichloro-4-nitrophenolate, C12H10N2O3Cl2 | ||||||||||
Molecular weight | 301.1 | ||||||||||
Unit cell | Monoclinic, space group P21/c, Z = 4 | ||||||||||
T (K) | 10 | 30 | 50 | 80 | 110 | 140 | 170 | 200 | 230 | 260 | 290 |
a (Å) | 7.7629(16) | 7.762(2) | 7.766(2) | 7.780(2) | 7.792(2) | 7.808(2) | 7.827(2) | 7.846(2) | 7.865(2) | 7.889(2) | 7.911(2) |
b (Å) | 23.525(5) | 23.519(7) | 23.537(7) | 23.541(7) | 23.565(7) | 23.582(7) | 23.619(7) | 23.647(7) | 23.681(8) | 23.720(8) | 23.760(8) |
c (Å) | 7.3900(18) | 7.387(2) | 7.395(2) | 7.405(2) | 7.419(2) | 7.432(2) | 7.447(2) | 7.460(2) | 7.472(3) | 7.487(3) | 7.502(3) |
β (°) | 111.49(2) | 111.47(2) | 111.51(2) | 111.59(2) | 111.70(2) | 111.81(2) | 111.88(2) | 112.01(2) | 112.10(2) | 112.22(3) | 112.33(2) |
V (Å3) | 1255.8(5) | 1255.0(7) | 1257.6(7) | 1261.1(7) | 1265.8(7) | 1270.4(6) | 1277.6(7) | 1283.2(7) | 1289.6(7) | 1297.0(8) | 1304.4(7) |
ρcalc (g cm−3) | 1.59 | 1.59 | 1.59 | 1.59 | 1.59 | 1.57 | 1.57 | 1.56 | 1.55 | 1.54 | 1.53 |
Number of refls | 10![]() |
9693 | 9008 | 7985 | 7478 | 6697 | 5943 | 5593 | 4943 | 4651 | 4056 |
Unique refls I > 2σ | 9069 | 8591 | 7930 | 6946 | 6546 | 5808 | 5219 | 4973 | 4370 | 4120 | 3638 |
N parameters | 268 | ||||||||||
R(F) | 0.0937 | 0.0991 | 0.0901 | 0.0794 | 0.0891 | 0.0962 | 0.0884 | 0.0861 | 0.0850 | 0.0883 | 0.0817 |
wR(F2) | 0.1775 | 0.1871 | 0.1601 | 0.1603 | 0.1638 | 0.1849 | 0.1649 | 0.1591 | 0.1606 | 0.1708 | 0.1525 |
Extinction coeff. | 33.8(7) | 35.1(8) | 36.3(7) | 37.2(8) | 35.5(9) | 40.9(9) | 45.0(10) | 48.3(12) | 46.7(13) | 50.6(14) | 53.4(15) |
Absorption coeff. | 2.770 + 0.0155λ |
Ab initio calculations based on the Born–Oppenheimer approximation were carried out at the B3LYP/6-311++G** level using the Gaussian 09 program.18 For the experimental structures the proton in the hydrogen bond was moved with the step of 0.05 Å, within the range 0.8–2 Å to calculate the potential energy surfaces. The vibrational time-independent Schrödinger equation for the OH coordinate was solved numerically using the two-boundary-conditions approach basing on the Numerov algorithm19 with the related variational-theory-based program.20 The structure optimized when the OH distance was kept constant and the experimental structures were used to perform the Natural Bond Orbital (NBO) analysis using NBO5.0 program.21 The wave function evaluated for each molecule was used as the input to the AIMALL program.22
![]() | ||
Fig. 1 Crystal structure of 3-methylpyridinium 2,6-dichloro-4-nitrophenolate from single crystal neutron diffraction. Thermal displacement ellipsoids are drawn at the 50% probability level. |
T [K] | N⋯O [Å] | OH [Å] | NH [Å] | CO [Å] | OHN [°] |
---|---|---|---|---|---|
290 | 2.533(5) | 1.386(8) | 1.167(8) | 1.271(6) | 167.0(7) |
270 | 2.532(4) | 1.400(8) | 1.152(8) | 1.270(5) | 165.8(7) |
230 | 2.532(4) | 1.398(8) | 1.151(8) | 1.272(5) | 166.8(7) |
200 | 2.531(4) | 1.406(8) | 1.141(7) | 1.269(4) | 167.3(7) |
170 | 2.531(3) | 1.405(6) | 1.143(5) | 1.275(3) | 166.8(5) |
140 | 2.531(3) | 1.412(5) | 1.137(5) | 1.271(3) | 166.3(5) |
110 | 2.523(3) | 1.436(5) | 1.116(5) | 1.281(4) | 162.7(5) |
80 | 2.531(3) | 1.418(5) | 1.131(5) | 1.277(3) | 166.3(5) |
50 | 2.537(2) | 1.429(4) | 1.126(4) | 1.275(2) | 166.3(4) |
30 | 2.5416(6) | 1.4210(3) | 1.1363(3) | 1.2759(4) | 167.24(2) |
10 | 2.5340(4) | 1.4264(2) | 1.1294(2) | 1.2726(3) | 164.93(2) |
According to the ΔpKa value of this complex, that equals 1.54, the OHN hydrogen bond should be very short, below 2.6 Å. Indeed, the O⋯N bond length in 3-methylpyridinium 2,6-dichloro-4-nitrophenol measured at different temperatures changes from 2.5416(6) to 2.523(3) Å. Very short distance between the proton donor and proton acceptor is usually connected with location of the proton in the middle of the distance between the proton donor and proton acceptor. Despite this expectation, the proton in the OHN hydrogen bond in 3-methylpyridinium 2,6-dichloro-4-nitrophenolate is shifted to the acceptor. Temperature increasing moves the proton to the central location in the hydrogen bond. Probably at higher temperatures the OH and NH distances can be equalized. Previously the X-ray structure of the complex of 2,6-dichloro-4-nitrophenol with 3-methylpyridine was investigated,13 and the N⋯O, OH, NH, and CO bond lengths were 2.544(4), 1.60(3), 0.93(3), and 1.285(4) Å respectively. Very short N⋯O bond length should be connected with location of the proton in the middle of the O⋯N distance, but both: neutron and X-ray structural measurements locate the proton at the acceptor. Neutron diffraction allows a precise determination of the proton location in the hydrogen bridge. X-ray diffraction measurement of the structure of 3-methylpyridinium 2,6-dichloro-4-nitrophenolate showed that the proton is shifted to the acceptor. Precisely determined proton position in the neutron measurement shows that the OH and NH distances are more equalized and the difference between OH and NH is less significant than the determined in X-ray diffraction. Nonsymmetrical location of the proton can be explained by special properties of the strong hydrogen bond. Unlike to molecular and ionic hydrogen bonds, the strong hydrogen bond is sensitive to temperature, environment, and isotopic substitution. Neutron diffraction confirms sensitivity of the proton location to temperature. Probably electric properties of the crystal of the investigated compound are responsible for shifting of the proton to the acceptor atom.
The cell parameters in Table 1 undergo systematic changes with temperature (a = 1E-06T2 + 0.0002T + 7.7559, R2 = 0.9972, b = 2E-06T2 + 0.0002T + 23.518, R2 = 0.9958, c = 4E-07T2 + 0.0003T + 7.3809, R2 = 0.9942, β = 3E-06T2 + 0.0021T + 111.42, R2 = 0.9938) so a phase transition is excluded. The distances connected with the OHN hydrogen bond listed in Table 2 change not continuously. Limited temperature changes of the N⋯O distance are close to 3σ value. Changes of NH and OH are more significant and linear correlation between NH and OH (NH = −0.9007OH + 2.411, R2 = 0.9488) reflects the proton moving in the hydrogen bond caused by temperature.
Correlations linking the geometrical hydrogen bond parameters have been investigated for very long time.6,23 Fig. 2 collects correlations for the neutron structures of O⋯N hydrogen bond complexes with the proton moving upon temperature compared with theoretical curves.24 The experimental data from the neutron structures with precisely determined proton position reproduces and confirms theoretical correlations linking the parameters of the hydrogen bridge: the Bond Order Reaction Curve (BORC) curve25 in Fig. 2a and the theoretical curve for the hydrogen bond length and deviation of the proton from central location between the proton donor and acceptor – Fig. 2b. When the experimental data confirm the theoretical correlation linking OH and NH, the theoretical curve in Fig. 2b calculated without taking into account anharmonicity of the hydrogen bridge vibrations, can be compared with the correlation obtained with the experimental neutron data. Anharmonicity correction can be neglected for molecular and ionic OHN hydrogen bonds but is significant for the strongest hydrogen bonds. Except anharmonicity correction also the quantum nature of the proton in short hydrogen bond influences the O⋯N distance and it can be expected that for the investigated complex this effect should be about −0.02 Å.26 Experimental neutron data for the complexes with proton moving in the hydrogen bridge allow reproducing the shape of the experimental curve and deliver the anharmonicity corrections.
![]() | ||
Fig. 2 Correlations of the geometric parameter of the strong OHN hydrogen bonds with proton moving in the hydrogen bridge: (a) – relation between NH and OH bond lengths, (b) – dependence of N⋯O bridge length as a function of the difference between OH and NH. The experimental neutron bridge parameters are compared with the theoretical curve calculated according to BORC curve (a) and p1 = exp(−(r1 − r1o)/b1), p2 = exp(−(r2 − r2o)/b2), p1 + p2 = 0, r1 = rOH, r2 = rNH, where p1 and p2 are the corresponding valence bond orders of the diatomic units, b1 = 0.371, r1o = 0.942, b2 = 0.385, r2o = 0.992 Å (b). Experimental data: • – pyridinium 2,4-dinitrobenzoate,27 o – 4-methylpyridine – pentachlorophenol complex,28 □ – 1![]() ![]() |
The CO bond length of the phenolic part of the complex, used as a measure of the degree of proton transfer in the intermolecular hydrogen bond, changes from 1.34 Å in molecular complexes to 1.27 Å in ionic complexes. For very short hydrogen bonds the CO bond length should be between these two values, although it has been shown that CO bond length as a function of OH depends on the substituents of the phenol ring.10 In 3-methylpyridinium 2,6-dichloro-4-nitrophenolate the CO bond length found in neutron measurement is typical of ionic hydrogen-bonded complexes, which is in agreement with the location of the proton at the acceptor when this same length for the X-ray structure equals 1.285(4) Å. The CO and OH bond lengths of the 2,6-dichloro-4-nitrophenol – 3-methylpyridine complex do not fulfill the general sigmoidal correlation10,31 linking the OH and CO bond lengths.
Both curves are characteristic of a strong hydrogen bond. A very broad minimum spreads over a broad range of OH. The lowest energy value is located at an OH distance which shifts the proton to the acceptor as in the X-ray structure or to proton donor as in the neutron structure. The change in the O⋯N bond length from 2.544(4) to 2.533(5) Å causes a significant difference in the OH distance, characteristic for the energy minimum. The shape of the curve does not make a determination of the height of the energy barrier between the two minimums possible because the second minimum is not well shaped. It is characteristic that the small difference in the geometry of the molecule with regard to the neutron and X-ray structures causes a shift of the minimum of the potential energy curve for proton movement in the hydrogen bridge. For the X-ray structure (O⋯N = 2.544(4) Å) the minimum on the potential energy curve is located at the OH bond length of 1.56 Å for and for the neutron structure (O⋯N = 2.533(5) Å) the minimum is shifted to the OH bond length of 1.05 Å. In both curves, the first energy level is located at higher energy compared with the top of the barrier on the potential energy curve, so the proton can move easily along the hydrogen bridge and it is very easy to shift it from the donor to the acceptor atom. For the curves in Fig. 3 the highest frequencies of vibrational transitions are equal 4523 cm−1 for the curves calculated for O⋯N distance characteristic for the neutron structure and 3975 cm−1 for the distance characteristic for the X-ray structure. Analysis of the vibrational properties of strong hydrogen bond should involve tunneling and quantum nuclear effects32 but even rough analysis exhibits difference between both vibrational transitions and emphasizes sensitivity of the experimental properties to the O⋯N distance.
To investigate the energy as a function of two parameters of the hydrogen bridge, the potential energy surfaces presented in Fig. 4 have been calculated. The potential energy surfaces have been calculated for relaxed geometry except the fixed distances.
Comparison of the energy plots in Fig. 4: O⋯N versus OH and O⋯N versus NH shows that it is very easy to change the location of the proton in the hydrogen bridge, but more energy is needed to elongate the O⋯N distance. For this reason, temperature can change the OH and NH bond lengths when the O⋯N distance is the same. The lowest energy is characteristic of the location of the proton at the oxygen atom. The second, very broad minimum appears at the acceptor, and the energy difference between the two minimums is only 2.37 kcal mol−1. The flat shape of the two-dimensional energy surface does not allow a determination of the barrier between the two minimums and illustrates that for a very strong hydrogen bond, each location of the proton along the OHN hydrogen bridge is possible, taking into account the shape of the energy surface.
NBO analysis performed previously12 for the experimental neutron structures of complex of 3,5-dinitrobenzoic acid with 3,5-dimethylpyridine and pyridinium 2,4-dinitrobenzoate with the proton moving upon temperature change showed significant rearrangement of molecular orbitals caused by the proton shifting in the hydrogen bond. In 3-methylpyridinium 2,6-dichloro-4-nitrophenolate the proton location in the hydrogen bond is sensitive to temperature but moving of the proton is rather limited and, independent of temperature, the hydrogen bond has ionic character. To follow through the proton moving from molecular to proton transfer complexes with particular noticing of the proton close to the hydrogen bond middle, NBO analysis has been performed for theoretical structures with continuous proton moving from the donor to acceptor. Also for the experimental neutron structures the NBO analysis has been done. Investigation of molecular orbitals of the proton donor and proton acceptor explains the mechanism of the proton transfer. The investigated NLMOs are presented in Fig. 5.
Percentage of the parent NBOs in NLMOs is a measure of delocalization of the orbitals that changes with moving of the proton in the hydrogen bond. Among the orbitals in Fig. 5 only the orbital “d” is insensitive to proton transfer. Other orbitals are strongly delocalized except the OH and NH bonds that exist only for typical molecular or typical ionic hydrogen bonds respectively. For strong hydrogen bonds orbitals representing OH and NH bonds are replaced by delocalized orbital “c” (Fig. 5). It is characteristic that slight elongation of OH bond is sufficient to replace the OH bond orbital with orbital “c” typical for strong hydrogen bond. For the strongest hydrogen bonds with location of the proton in between the donor – acceptor distance, both oxygen lone pairs participating in the proton transfer as well as the lone pair of the acceptor are characterized by similar delocalization. The points representing the experimental structures confirm the theoretical correlations in Fig. 6.
![]() | ||
Fig. 6 % of the parent NBO orbitals in delocalized NLMO. a, b, c, e, f refers to the orbitals in Fig. 5. The dots represent the experimental neutron structures of 3-methylpyridinium 2,6-dichloro-4-nitrophenolate. |
Proton moving in the hydrogen bond is reflected in the participation of the proton in the NLMOs of the donor and acceptor. Percentage of the proton in orbitals presented in Fig. 7 illustrates the mechanism of proton transfer. For the strong hydrogen bond the proton is not engaged in OH or NH bond but participates in donor and acceptor lone pairs. Percentage of proton in NLMO is related to delocalization of the lone pairs and if the molecular orbital is more delocalized, the contribution of the proton is more significant.
Shifting of the proton to the middle of the O⋯N distance is connected with similar delocalization of the molecular orbitals participating in proton transfer so the proton can be similarly engaged in the electron cloud of the proton donor and acceptor.
Limited amount of the experimental structures with the proton moving is not sufficient to follow the proton transfer from donor to acceptor. Combination of the data for experimental structures with correlations obtained theoretically allows investigation of the proton transfer mechanism for complete range of the proton transfer.
The eigen values of the Hessian matrix of the electron density at BCP (λ1, λ2, λ3) are a source of information on stability of the bond and energetic properties of the electrons at BCP. The ellipticity (ε = (λ1/λ2) − 1) of electron density at BCP describes the deviation from spherical concentration of the electron density at the BCP. The kinetic and potential energy density can be extrapolated as: G(r)BCP = 15.3λ3 and V(r)BCP = 35.1(λ1 + λ2).38 The potential energy of the electrons at BCP expresses the pressure exerted on the electrons by other electrons; the kinetic energy expresses the mobility of the electrons at the BCP and their pressure exerted on other electrons. The Laplacian of the electron density at BCP can be used to classify the interatomic interactions as the shared-shell (∇2ρ(r) < 0) and the closed-shell interaction (∇2ρ(r) > 0). For the closed shell interaction and the values of |V(r)BCP|/G(r)BCP between 1 and 2 the hydrogen bond has partially covalent character.39
The AIM parameters characterizing the OH and NH bond critical points for the neutron structures of 3-methylpyridinium 2,6-dichloro-4-nitrophenolate are collected in Table 3. In all the analyzed structures, the OHN hydrogen bond is strong so the values of the electron density are significantly higher than the value of 0.002 [a.u.]. Other criteria of hydrogen bond existence are also fulfilled. The values of the electron density, Laplacian, potential, kinetic energy of the electrons at the BCPs of the O–H and N–H bonds indicate existence of strong hydrogen bond.40 The ratio |V(r)BCP|/G(r)BCP for OH bond indicates a noncovalent interaction with participation of covalent character.
Temp [K] | Bond | ρ(r) | ∇2ρ(r) | G(r)BCP [kJ mol−1] | V(r)BCP [kJ mol−1] | ε |
---|---|---|---|---|---|---|
10 | O–H | 0.0931 | 0.126 | 8.496 | −15.079 | 0.012 |
N–H | 0.239 | −1.027 | 11.184 | −61.695 | 0.021 | |
30 | O–H | 0.094 | 0.121 | 8.580 | −15.437 | 0.014 |
N–H | 0.235 | −0.981 | 11.142 | −59.979 | 0.021 | |
50 | O–H | 0.093 | 0.128 | 8.466 | −14.946 | 0.014 |
N–H | 0.241 | −1.048 | 11.197 | −62.671 | 0.021 | |
80 | O–H | 0.095 | 0.124 | 8.684 | −15.563 | 0.013 |
N–H | 0.238 | −1.014 | 11.216 | −61.316 | 0.021 | |
110 | O–H | 0.091 | 0.134 | 8.395 | −14.563 | 0.012 |
N–H | 0.248 | −1.119 | 11.310 | −65.230 | 0.020 | |
140 | O–H | 0.097 | 0.120 | 8.805 | −15.974 | 0.013 |
N–H | 0.234 | −0.972 | 11.250 | −59.747 | 0.021 | |
170 | O–H | 0.098 | 0.116 | 8.938 | −16.434 | 0.013 |
N–H | 0.230 | −0.933 | 11.154 | −58.340 | 0.021 | |
200 | O–H | 0.098 | 0.119 | 8.908 | −16.272 | 0.012 |
N–H | 0.232 | −0.950 | 11.146 | −58.915 | 0.021 | |
230 | O–H | 0.100 | 0.112 | 9.041 | −16.827 | 0.013 |
N–H | 0.225 | −0.883 | 11.073 | −56.378 | 0.021 | |
260 | O–H | 0.099 | 0.111 | 8.986 | −16.725 | 0.012 |
N–H | 0.225 | −0.879 | 11.068 | −56.237 | 0.021 | |
290 | O–H | 0.103 | 0.098 | 9.244 | −17.750 | 0.012 |
N–H | 0.216 | −0.792 | 10.978 | −52.976 | 0.021 |
In Fig. 8 is presented dependency of ρ(r)OH and ρ(r)NH on the OH bond length calculated for 3-methylpyridinium 2,6-dichloro-4-nitrophenolate with gradual shift of the proton realized by change of the OH bond length and optimization of other structural parameters. The obtained curves are identical with previously analyzed dependencies of AIM parameters for the complex of trimethylamine with 2,6-dichloro-4-nitrophenol.41 Analogous dependencies for other phenols are sensitive to the proton donor. Replacing of trimethylamine with 3-methylpyridine does not change the correlation of ρ(r)OH and ρ(r)NH on the OH bond length that is in agreement with the higher sensitivity of the proton donor than the proton acceptor to the shifting of the proton.
![]() | ||
Fig. 8 Theoretical electron density (ρ(r)) at the OH and NH (3,−1) critical point as a function of OH distance. The values for experimental neutron structure are marked as Δ. |
Electron density values calculated for the experimental neutron structures of 3-methylpyridinium 2,6-dichloro-4-nitrophenolate are located on the curve obtained for the optimized structures. This fact confirms the general character of the dependence of the electron density on the distance common for calculated and experimental results.
In Fig. 9 electron densities at BCP of the OH and NH bonds calculated for all intermolecular OHN neutron structures measured at different temperatures have been compiled. The electron densities at OH and NH BCP have been calculated and correlated with the difference of OH and NH bond length used as a measure of the hydrogen bond strength. For the strongest hydrogen bonds the difference between OH and NH bond length should be close to zero although the crossing point at OH–NH equal −0.0468 Å and the electron density of 0.1579 [a.u.] reflects nonequivalence of OH and NH bonds. The correlations in Fig. 9 include the data for phenols and benzoic acids and are common for OHN complexes independent of the proton donor. Electron densities at bond critical points are common for all OH proton donors and N proton acceptors and insensitive to particular molecules in the intermolecular hydrogen bond.
![]() | ||
Fig. 9 Relationship of ρ(r) at the OH and NH BCPs with OH – NH difference. • – pyridinium 2,4-dinitrobenzoate,27 o – 4-methylpyridine – pentachlorophenol complex,28 □ – 1![]() ![]() |
Footnote |
† Electronic supplementary information (ESI) available. CCDC 1060082–1060092. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c5ra06733b |
This journal is © The Royal Society of Chemistry 2015 |