Eva
Gougoula
a,
Charlotte N.
Cummings
b,
Chris
Medcraft
c,
Juliane
Heitkämper
d and
Nicholas R.
Walker
*b
aDeutsches Elektronen-Synchrotron DESY, Notkestr. 85, 22607 Hamburg, Germany
bChemistry-School of Natural and Environmental Sciences, Newcastle University, Bedson Building, Newcastle-upon-Tyne, NE1 7RU, UK. E-mail: nick.walker@newcastle.ac.uk
cSchool of Chemistry, UNSW Sydney, Sydney, New South Wales 2052, Australia
dInstitute of Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, D-70569, Stuttgart, Germany
First published on 6th May 2022
Broadband microwave spectra have been recorded between 7.0 and 18.5 GHz for N-methylimidazole⋯H2O and 2-methylimidazole⋯H2O complexes. Each complex was generated by co-expansion of low concentrations of methylimidazole and H2O in argon buffer gas. The rotational spectra of five isotopologues of each complex have been assigned and analysed to determine rotational constants (A0, B0, C0), centrifugal distortion constants (DJ, DJK) and parameters that describe the internal rotation of the CH3 group. The results allow the determination of parameters in the (r0) molecular geometry of each complex. H2O is the hydrogen bond donor and the pyridinic nitrogen of imidazole is the hydrogen bond acceptor in each case. The ∠(O–Hb⋯N3) angles are 177(5)° and 166.3(28)° for N-methylimidazole⋯H2O and 2-methylimidazole⋯H2O respectively. These results are consistent with the presence of a weak electrostatic interaction between the oxygen atom of H2O and the hydrogen atom (or CH3 group) attached to the C2 carbon atom of imidazole, and with the results of density functional theory calculations. The (V3) barrier to internal rotation of the CH3 group within N-methylimidazole⋯H2O is essentially unchanged from the value of this parameter for the N-methylimidazole monomer. The same parameter is significantly higher for the 2-methylimidazole⋯H2O complex than for the 2-methylimidazole monomer as a consequence of the weak electrostatic interaction between the O atom and the CH3 group of 2-methylimidazole.
A recent work by our group described the geometries of two isomers of a complex formed between imidazole and H2O.4 One of these isomers (denoted as imid⋯H2O) contains a hydrogen bond between H2O and the pyridinic nitrogen, N3, of the imidazole ring (atom numbering of imidazole is shown in Fig. 1). This work and recent results from theoretical calculations5 confirm similarities between the structure of imid⋯H2O and those of other complexes formed of water and a heteroaromatic molecule such as thiazole,6 isoxazole7 or pyrimidine.8 Specifically, structures of these complexes usually feature H2O acting as a proton donor, a heteroatom of the ring acting as a proton acceptor and a secondary interaction where the oxygen atom of water binds to a hydrogen atom of a C–H group on the ring. The present work will show that the geometry of a complex formed of 2-methylimidazole (2-MI) and water is determined by interactions which are very similar to those described above. However, in this case, the secondary interaction is between the oxygen atom of H2O and the CH3 group attached to C2 of the imidazole ring.
Fig. 1 Equilibrium (re) geometries of N-MI⋯H2O (left) and 2-MI⋯H2O (right) calculated at the ωB97XD/aug-cc-pVQZ level. |
An aim of the present work is to perform the first experimental determination of the change in a (V3) barrier to internal rotation of CH3 which is caused by formation of an electrostatic interaction between H2O and a CH3 group. It thus begins experimental investigation of how electrostatic interactions of CH3 with one or more H2O molecules influence barriers to internal rotation. Alongside works that have explored the geometrical preferences of water-containing molecular complexes, the present work contributes to knowledge of the effects of solvation on conformational flexibility. It will be shown that V3 of the CH3 group of 2-MI increases significantly on formation of a complex between H2O and 2-MI. In contrast, V3 of CH3 in N-MI⋯H2O (where N-MI is N-methylimidazole) is very similar to the value of this parameter in the isolated N-MI monomer.9 It will be shown that these observations are consistent with the geometries of the complexes. An electrostatic (through-space) interaction between the H2O molecule and CH3 group is present in the geometry of 2-MI⋯H2O but not in N-MI⋯H2O. Experimentally-determined structural parameters that define the position and orientation of the water molecule will also be reported for both the complexes examined by this work.
The rotational spectra of N-MI⋯H2O and 2-MI⋯H2O were recorded between 7.0 and 18.5 GHz. The CP-FTMW spectrometer mixes a microwave pulse that linearly sweeps from 12.0 to 0.5 GHz over a duration of 1 μs against the 19 GHz reference signal of a Phase-locked Dielectric Resonant Oscillator (PDRO). Chirped pulses are generated from a 20 GS s−1 arbitrary waveform generator (AWG) (Tektronix AWG 7102). A low pass filter selects the 7.0–18.5 GHz sideband which is amplified by a 300 W Traveling-Wave Tube Amplifier (TWT) prior to its introduction from a horn antenna into the vacuum chamber perpendicular to the direction of the expanding gas jet. Subsequent to the molecular polarisation, the free induction decay (which has duration of 20 μs) of the molecular emission is detected by a second horn antenna and digitally recorded by a 100 GS s−1 oscilloscope (Tektronix DPO72304XS). Successive free induction decays are co-added in the time domain. The polarisation pulse and free induction decay are each sufficiently short that eight distinct measurements of the broadband microwave spectrum are routinely performed following each gas introduction pulse while taking advantage of the “fast frame” mode of the oscilloscope. A high resolution window function was used by the Fourier transform employed herein such that linewidths of transitions in the frequency domain spectrum are approximately 100 kHz for a well-isolated line at full-width half-maximum. This linewidth correlates with an estimated accuracy of 10 kHz in the measurement of line centre frequencies. Phase coherence in the time domain and accuracy in transition frequencies were provided by an Rb-clock (SRS FS725) to which the AWG, the PDRO and the oscilloscope were phase-locked.
Optimisations of the geometries of N-MI⋯H2O and 2-MI⋯H2O were performed using the Gaussian09 package.13 The harmonic hybrid functional14–16 of Becke, Lee, Yang, Parr, B3LYP, in conjunction with Grimme's dispersion correction effects17 and damping function,18 D3BJ, was initially used alongside Dunning's19,20 augmented triple-ζ aug-cc-pVTZ basis set. Geometry optimisations were subsequently performed with the long range corrected hybrid functional,21 ωB97X-D, of Chai and Gordon and Dunning's augmented quadrupole-ζ basis set with tight convergence criteria. The same level of theory was used during a previous study of imid⋯H2O.4 Structural parameters that describe the interaction between H2O and imidazole were assumed equal to those calculated for imid⋯H2O in the starting point of the geometry optimisation for each of N-MI⋯H2O and 2-MI⋯H2O. In each case, it was assumed that the molecule is oriented such that ∠(O–Hb⋯N3–C2) = 0°. The V3 barrier to internal rotation was calculated by scanning the ∠(H–C6–N1–C2) or ∠(H–C6–C2–N3) dihedral angle as appropriate to each complex. The results of the geometry optimisation for each complex are shown in Table S1 (ESI†) and Fig. 1 while the calculated rotational constants (Ae, Be, Ce), nuclear quadrupole coupling constants, dipole moment components (|μa|, |μb| and |μc|) and (V3) barriers to internal rotation are summarized in Table S2 (ESI†) which also shows the differences between calculated and experimentally-determined parameters. Calculations of anharmonic force fields employed the B3LYP functional with the D3BJ correction and the jun-cc-pVDZ basis set.
A previous experimental study observed two isomers of a 1:1 complex formed between water and imidazole.4 Water was the hydrogen bond donor in the first of these isomers (denoted as imid⋯H2O) and the hydrogen bond acceptor in the second (denoted as H2O⋯imid). The bifunctional nature of the imidazole ring implies that the method employed herein might allow the generation of H2O⋯2-MI as well as 2-MI⋯H2O. Ultimately, conclusive proof that the observed spectrum should be assigned to 2-MI⋯H2O was obtained by measuring the shifts in rotational constants on isotopic substitution. Spectra were measured for 2-MI⋯H216O, 2-MI⋯H218O, 2-MI⋯DOH, 2-MI⋯HOD and 2-MI⋯D2O isotopologues of 2-MI⋯H2O. All measured isotopic shifts were highly consistent with assignment of 2-MI⋯H2O as the carrier of the observed spectrum. The experiments did not reveal evidence for the formation of H2O⋯2-MI. Beyond the spectra identified and discussed above, a number of transitions observed during the present work assign to products generated by the fragmentation of imidazole. The same range of fragmentation products was observed during previous studies4,22 of (imid)2 and imid⋯H2O which were performed by the same method.
Partially resolved hyperfine structure (Fig. 2, bottom panels) was observed as consistent with the presence of two nitrogen nuclei (I = 1 for 14N) in each complex. Initially, A-species transition frequencies were fit with Watson's S-reduced Hamiltonian,23 implemented in PGOPHER;24
(1) |
The nuclear quadrupole moment dyadic is represented by Q and the dyadic of the electric field gradient is ΔE. Matrix elements were constructed in the coupled asymmetric rotor basis, J + IN1 = F1, F + IN3 = F2, and diagonalised in blocks of the quantum number, F. Components of the nuclear quadrupole coupling tensor of each quadrupolar nitrogen nucleus are denoted by χaa and (χbb − χcc) with nitrogen atoms numbered according to the convention when labelling heteroaromatic molecules with N1 and N3 denoting the pyrrolic and pyridinic nitrogen atoms respectively. Hyperfine structure was observed to be extensively blended but values of χaa(N1) and χaa(N3) were fitted for every isotopologue of each complex. Values of [χbb(N1) − χcc(N1)] for N-MI⋯DOH, N-MI⋯HOD and N-MI⋯D2O were fixed equal to the experimental result for N-MI⋯H216O. Values of [χbb(N3) − χcc(N3)] were held fixed for most isotopologues of each complex. The procedure yielded satisfactory agreement between simulated and experimentally-observed line shapes (Fig. 2, lower panel), and between the experimentally-determined and DFT-calculated results. The electron distribution at each nitrogen atom is confirmed to be very similar to that for the corresponding atom of the imidazole monomer (Table S5, ESI†). The final cycle of each spectroscopic fit is provided within the supplementary data. The standard deviations of all fits of A-species transitions (Tables S3 and S4, ESI†) are between 15 and 35 kHz. Spectroscopic parameters are thus determined with high precision in spite of the line blending and assumed values of hyperfine parameters. This is consistent with previous work which found that good agreement was achieved between DFT-calculated and experimental results for χaa and (χbb − χcc) for nitrogen atoms within N-heterocyclic rings.
Subsequent to the identification and assignment of E-species transitions, combined fits of spectroscopic parameters to the frequencies of both A- and E- species transitions were performed by two different methods. Previous papers have noted that the determined value of V3 can vary depending on the Hamiltonian model employed. Fits performed in XIAM34 employ a Hamiltonian that is constructed in the principal axis system. XIAM treats each torsional state separately without accounting for interactions between different torsional states. The BELGI-Cs program35 employs the alternative reference framework of the rho axis method and allows the inclusion of higher-order terms. Parameters were fitted using each of XIAM and BELGI-Cs to explore dependencies of the results on the models employed. The experiments probed complexes in their vibrational ground states preventing the independent determination of F0 and V3. The former was therefore fixed at the value for the corresponding methylimidazole monomer in the fits performed for each of N-MI⋯H2O and 2-MI⋯H2O. It was also assumed that the axis of the CH3 rotor lies in the ab plane. Nuclear quadrupole coupling of multiple quadrupolar nuclei is not currently implemented by either XIAM or BELGI-Cs so a bespoke procedure was used to establish line centre frequencies for A- and E-species transitions. First, the rotational and centrifugal distortion constants determined by fitting A species transitions (those in Tables S3 and S4, ESI†) were used to determine the unsplit line centres of A-species rotational transitions. Next, the shift of each E-species transition was measured relative to its A-species counterpart while assisted by PGOPHER24 visualisation tools that aid spectral assignment. Specifically, these tools allow (i) hyperfine structure to be toggled on or off in the displayed simulation, (ii) display of a simulated spectrum which correctly mirrors the instrument resolution and (iii) simulated and experimental spectra to be overlaid and/or offset in order to match simulated features with observations. The precisely-measured interval between each A-species transition and its E-species counterpart was thus used to determine the unsplit line centre frequency for each E-species transition. Fits performed using XIAM determined rotational constants and centrifugal distortion constants in the principal axis system, those performed using BELGI-Cs determined constants in the rho axis system. The linelist used by each BELGI-Cs fit was the same as that used by the XIAM fit for the corresponding isotopologue. A summary of results for one isotopologue of each complex is provided in Table 1 with full details for all isotopologues provided in Tables S6–S9 (ESI†).
N-MI⋯H216O | 2-MI⋯H216O | |||
---|---|---|---|---|
XIAM | BELGI-Cs | XIAM | BELGI-Cs | |
a Rotational constants in the principal axis system after transformation of the rho axis system inertia tensor determined by BELGI-Cs. b Values in the rho axis system. Centrifugal distortion constants determined by BELGI-Cs cannot be directly compared with those determined by XIAM because of the different models employed. c Could not be determined by fitting and therefore fixed to the result determined by XIAM. d Fixed to the result obtained for the N-MI⋯D2O isotopologue (Table S7, ESI). | ||||
A 0 (MHz) | 5010.78(21) | 5011(10)a | 4233.93(21) | 4233.2(50)a |
B 0 (MHz) | 1409.5628(67) | 1409.5(15)a | 1732.5621(61) | 1731.00(95)a |
C 0 (MHz) | 1107.2884(60) | 1107.3099(43)a | 1240.3104(66) | 1238.69(18)a |
D J (kHz) | 2.392(33) | 3.111(49)b | 0.233(81) | [0.233]bc |
D JK (kHz) | −18.57(40) | −20.6(5)b | 7.28(78) | 4.33(71)b |
d 1 (kHz) | 0.820(42) | [0.959]bd | — | — |
F 0 (GHz) | [157.929] | [157.929] | [157.690] | [157.690] |
V 3 (cm−1) | 182.23(10) | 173.6(16) | 154.99(8) | 150.68(90) |
D π2J(kHz) | — | — | −114(7) | — |
F v (MHz) | — | — | — | 5.49(60) |
∠(i, b) (°) | 47.16(10) | 46.3(3) | 36.55(6) | 36.56(5) |
σ RMS (kHz) | 50 | 55 | 64 | 79 |
N A/NE | 26/19 | 26/19 | 19/14 | 19/14 |
The internal rotation parameters included in XIAM fits were the barrier to internal rotation of the CH3 group, V3, and the angle between the axis of the CH3 rotor and the b inertial axis which is denoted by ∠(i, b). A term which couples the methyl rotation with the overall molecular rotation and centrifugal distortion, Dπ2J, was included when fitting the data for 2-MI⋯H2O. The BELGI-Cs fits determined the off-diagonal component of the inertia tensor in the rho axis system, DAB, and the internal rotation parameters, V3 and ρ for each isotopologue of each complex. A term that describes the J(J + 1) dependence of the V3 barrier, FV, was also fitted for isotopologues of 2-MI⋯H2O. The σRMS of the XIAM and BELGI-CS fits were higher than those obtained when fitting the A-species transitions in PGOPHER and higher than would typically be expected given the experimental linewidth of 100 kHz because of blending of different hyperfine components in the observed spectra and the need to estimate rotational line centre frequencies for the XIAM and BELGI-Cs fits.
The rotational constants and distortion constants output by BELGI-CS fits are in the rho axis system so the inertia tensor is not diagonalised and rotational constants cannot be directly compared with those determined in the principal axis system. The inertia tensor determined by the BELGI-Cs fits was transformed into the principal axis system for the purposes of the comparison shown in rows 1–3 of Table 1. The value of V3 for the CH3 group of N-MI⋯H216O is determined to be 182.23(10) cm−1 by XIAM with other isotopologues having V3 ranging from 182.08(8) cm−1 to 182.41(7) cm−1. These results are very similar to the value of this parameter for the N-MI monomer9 which was determined to be 185.104(11) cm−1 through fitting in XIAM. The same parameter for N-MI⋯H216O is 173.6(16) cm−1 when fitted using BELGI-Cs with other isotopologues having V3 ranging from 173.2(2) cm−1 to 177.2(7) cm−1. The XIAM fits yielded a result for V3 of 154.99(8) cm−1 for 2-MI⋯H216O with results spanning from 152.65(4) cm−1 to 156.02(12) cm−1 for other isotopologues. The fits performed using BELGI-CS determined the parameter to be 150.68(90) cm−1 for 2-MI⋯H216O with values determined for other isotopologues ranging from 152.58(11) cm−1 to 155.5(24) cm−1. The results determined for V3 of the CH3 group in 2-MI⋯H216O are all significantly higher than the value of the parameter for the 2-MI monomer which was reported to be 122.7529(38) cm−1.9 The reasons for this difference will be explored further in Section 4.
Values of ∠(i, b) fitted using XIAM and BELGI-Cs are in excellent agreement for all isotopologues. This parameter ranges from 46.06(3) to 47.54(7)° in the fits of internal rotor parameters for N-MI⋯H2O whereas the value implied by the r0 fit of the molecular geometry (to be described in Section 3.3) is 46.1°. The values of ∠(i, b) determined in the fits of internal rotor parameters for isotopologues of 2-MI⋯H2O range from 34.86(3) to 36.56(5)° while the r0 fit of the molecular geometry implies a result of 35.3° for this parameter. This high level of agreement provides confidence in the results of the internal rotor fits while also justifying assumptions that will be made to determine the r0 geometries of N-MI⋯H2O and 2-MI⋯H2O.
Δ0 = Icc − Iaa − Ibb | (2) |
The small changes in Δ0 when H2O attaches to either N-MI or 2-MI thus confirm that each complex forms in a geometry where all heavy atoms lie within the ab plane. In general, contributions of in-plane vibrations to Δ0 are negative whereas out-of-plane vibrations make positive contributions. The magnitudes of changes in Δ0 on addition of H2O to each of imidazole,4 thiazole6 and isoxazole7 are −0.02, −0.05 and +0.2 u Å2 respectively which are broadly consistent with the observations of the present work. In striking contrast, the change in Δ0 when attaching H2O to pyridine36 has been observed to be greater than 8 u Å2. There are significant changes in Δ0 upon isotopic substitution at Hnb which are presented in Tables S6 and S8 (ESI†) and discussed further below.
Determination of the rotational constants for five isotopologues of each complex allows the calculation of substitution (rs) coordinates for the hydrogen and oxygen atoms of the water molecule using Kisiel's program “KRA” available from the PROSPE website.37 The rotational constants obtained from the XIAM fits of A- and E-species transitions were used to calculate rs coordinates for Hb, Hnb and O which are shown alongside Costain errors38 in Table 2. The Kraitchman method determines the magnitudes of atomic coordinates while the signs are those inferred from the results of the DFT calculations. The results unambiguously confirm that the oxygen atom lies in the ab plane in each complex and that the molecular geometry of each of N-MI⋯H2O and 2-MI⋯H2O contains a hydrogen bond between water (acting as a proton donor) and the pyridinic nitrogen.
Method | a/Å | b/Å | c/Å | |
---|---|---|---|---|
a r e geometries are calculated at the ωB97XD/aug-cc-pVQZ level. b Numbers in parentheses are Costain errors.38 c Imaginary values were obtained by the rs method for coordinates indicated in square brackets which are assumed equal to zero. | ||||
N-MI⋯H2O | ||||
Hb | r e (calc.)a | 2.6691 | –0.1954 | –0.0341 |
r s (exp.) | 2.7785(8)b | [0]c | 0.05(4) | |
O | r e (calc.) | 3.5881 | –0.5061 | –0.0938 |
r s (exp.) | 3.4036(5) | –0.650(3) | [0] | |
Hnb | re (calc.) | 3.9972 | –0.2262 | 0.7223 |
rs (exp.) | 3.9163(7) | –0.728(4) | 0.462(6) | |
2-MI⋯H2O | ||||
Hb | r e (calc.) | 2.1105 | –0.7120 | –0.0300 |
r s (exp.) | 2.064(1) | –0.784(3) | 0.03(7) | |
O | r e (calc.) | 3.0707 | –0.5535 | –0.0549 |
r s (exp.) | 3.0709(6) | –0.538(3) | [0] | |
Hnb | r e (calc.) | 3.4265 | –1.0444 | 0.6827 |
r s (exp.) | 3.5032(6) | –1.181(2) | 0.359(6) |
The r0 method39 (as implemented within Kisiel's STRFIT program37), which allows a determination of structural parameters from zero-point rotational constants, was used to determine the orientation of the water molecule relative to the imidazole ring. The geometrical parameters that are internal to the N-MI and 2-MI monomers9 must first be assumed equal to their values in the DFT-calculated (re) geometries of N-MI⋯H2O and 2-MI⋯H2O respectively. The experimentally-determined rotational constants (A0, B0 and C0) determined by XIAM now allow fitting of the length of the intermolecular hydrogen bond, r(Hb⋯N3), and two angles which define the position and orientation of the water molecule, ∠(Hb⋯N3–C2) and ∠(O–Hb⋯N3), for each of N-MI⋯H2O and 2-MI⋯H2O. Results were determined while assuming ∠(Hnb–O–Hb⋯N3) = 180° and ∠(O–Hb⋯N3–C2) = 0° for each complex. These dihedral angles lead to sets of structural parameters (shown in Table 3) that are highly consistent with the results of the DFT calculations. An alternative assumption that ∠(O–Hb⋯N3–C2) = 180° for 2-MI⋯H2O would yield r(Hb⋯N3) = 1.957(23) Å, ∠(Hb⋯N3-C2) = 107.3(18)° and ∠(O-Hb⋯N3) = 154.4(60)° which are somewhat different from the DFT-calculated results. The assumption that ∠(O–Hb⋯N3–C2) = 180° for N-MI⋯H2O would yield ∠(O–Hb⋯N3)= 127(6)° which is certainly not consistent with the DFT calculated results. Fits performed while assuming that ∠(Hnb–O–Hb⋯N3) = 0° do not converge for either complex.
Molecule | Parameter | Method | Value |
---|---|---|---|
a r e geometries are calculated at the ωB97XD/aug-cc-pVQZ level. b Numbers in parentheses are one standard deviation in units of the final significant figure. | |||
N-MI⋯H2O | r(Hb⋯N3)/Å | r e (calc.)a | 1.897 |
r 0 (exp.) | 1.922(4)b | ||
∠(Hb⋯N3–C2)/° | r e (calc.) | 115.7 | |
r 0 (exp.) | 101.0(16) | ||
∠(O–Hb⋯N3)/° | r e (calc.) | 170.1 | |
r 0 (exp.) | 177(5) | ||
2-MI⋯H2O | r(Hb⋯N3)/Å | r e (calc.) | 1.8807 |
r 0 (exp.) | 1.923(5) | ||
∠(Hb⋯N3–C2)/° | r e (calc.) | 114.3 | |
r 0 (exp.) | 116.9(9) | ||
∠(O–Hb⋯N3)/° | r e (calc.) | 165.4 | |
r 0 (exp.) | 166.3(28) |
The DFT calculations of the present work identify that Hnb lies outside the plane of the imidazole ring in each of N-MI⋯H2O and 2-MI⋯H2O in the equilibrium (re) geometry. Previous works4,6 have found that similar “free” hydrogen atoms in related complexes undergo rapid zero-point vibrational motions such that they cannot be precisely located by r0 or rs methods. Comparing the results obtained by the rs and re methods, the a and b coordinates of Hb and O are highly consistent. The Kraitchman method is less precise when used to calculate the position of an atom that is close to an inertial axis. N-MI⋯H2O represents a particularly challenging case because both hydrogen atoms are very near to the a inertial axis. Changes in zero-point vibrational motions on isotopic substitution also introduce inaccuracies and these are particularly severe when the rs coordinates of a hydrogen atom are determined. These sources of inaccuracy explain small variations in the results for the b-coordinate of Hb obtained by the re and rs methods for each complex. Significant variation in the |b| and |c| coordinates of Hnb as determined by the r0, rs and re methods results from the rapid, zero-point vibrational motions mentioned earlier. Variation in the inertial defects measured for different isotopologues is informative here (Tables S6 and S8, ESI†). For each of N-MI⋯H2O and 2-MI⋯H2O, the Δ0 of isotopologues that contain H216O, DOH or H218O are very similar (within 0.02 u Å2 of each other) supporting the conclusion that Hb and O are localised within the plane of the imidazole ring in the zero-point state. However, isotopologues that are deuterated in the Hnb position have |Δ0| which are between 0.2 and 0.4 u Å2 higher than those of isotopologues which are not deuterated at Hnb. Evidently, out-of-plane contributions to zero-point vibrational motions mean that Hnb is not localised to the ab plane on the timescale of the molecular rotation. To facilitate broader comparisons, Tables S10 and S11 (ESI†) provide atomic coordinates calculated by the r0 method.
The analysis thus far has used the experimentally-determined rotational constants for the zero-point vibrational state of each complex to determine r0 and rs geometrical parameters. These r0 and rs parameters cannot necessarily be expected to agree with re (equlibrium) parameters with very high precision because the former are influenced by zero-point vibrational motions. However, it is possible to calculate semi-experimental equilibrium rotational constants5,40 (Bie) from experimentally-determined rotational constants (Bi0), where i indicates permutation over the a, b, c inertial axes, using
(3) |
During the present work, calculations of anharmonic force fields to yield values of αir for each of N-MI⋯H2O and 2-MI⋯H2O were performed using the B3LYP functional, the partially-augmented double-ζ jun-cc-pVDZ basis set and Grimme's DFT-D3 scheme with the Becke–Johnson (BJ) damping function for the treatment of dispersion effects. Prior to calculation of the anharmonic force field, the geometry of the complex was re-optimised at this level of theory. This scheme for calculation of the anharmonic force field was previously used by Melli et al during their recent study of imid⋯H2O. Values of αir and hence the semi-experimental equilibrium rotational constants, Bie, were calculated for isotopologues of each complex as shown in Tables S12 and S13 (ESI†). A least squares fit of the values of r(Hb⋯N3), ∠(Hb⋯N3–C2) and ∠(O–Hb⋯N3) to the semi-experimental equilibrium rotational constants (Bie) determined for the various isotopologues of 2-MI⋯H2O yielded values of 1.912(10) Å, 116.6(12)° and 160.8(38)° respectively for these parameters. A determination of rSEe parameters for N-MI⋯H2O by the same method yielded r(Hb⋯N3) = 1.904(20) Å, ∠(Hb⋯N3–C2) = 104.0(42)° and ∠(O–Hb⋯N3) = 168(13)°. For each complex, the rSEe value of r(Hb⋯N3) is in better agreement with the DFT-calculated re result than is the r0 parameter. However, particularly for N-MI⋯H2O, the determined rSEe angles are associated with high uncertainties. The discussion of Section 4 will therefore compare values of r0 parameters determined with high precision for imid⋯H2O, N-MI⋯H2O and 2-MI⋯H2O. Differing zero-point vibrational motions within each complex mean that precise agreement should not necessarily be expected in comparisons of r0 parameters determined for different complexes.
Melandri et al., McGlone et al., Caminati et al., McKenzie et al. and Li et al. identified non-linear hydrogen bonds in pyrimidine⋯H2O,8 isoxazole⋯H2O,7 pyridazine⋯H2O,41 pyrazine⋯H2O,42 pyridine⋯H2O36 and thiazole⋯H2O6 respectively. Our group recently reported that imid⋯H2O4 contains a non-linear hydrogen bond where ∠(O–Hb⋯N3) = 172.1(26)° in the r0 geometry. In each of the studies listed above, the non-linearity of the hydrogen bond results from a weak electrostatic interaction between the oxygen of H2O and a hydrogen on the C2 position of the aromatic ring. The local environments at the pyridinic nitrogen atoms of imidazole and N-methylimidazole are very similar. Unsurprisingly, therefore, the r0 results of ∠(Hb⋯N3–C2) = 101.1(16)° and ∠(O–Hb⋯N3) = 177(5)° reported herein for N-MI⋯H2O are similar to the results reported for imid⋯H2O. The oxygen atom is a distance of 3.081(6) Å from the hydrogen attached to C2 in the r0 geometry of N-MI⋯H2O while the separation between these two atoms in imid⋯H2O is 3.13(8) Å. The separation between oxygen and the hydrogen atoms of the CH3 in 2-MI⋯H2O depends on the value of the dihedral angle, ∠(N3–C2–C6–H), which is scanned by the CH3 internal rotation. When ∠(N3–C2–C6–H) = 0°, the separation between oxygen and the nearest hydrogen is 2.566(3) Å while the assumption that ∠(N3–C2–C6–H) = 60° places two hydrogen atoms only 3.175(3) Å from the oxygen atom. Evidently, oxygen interacts with its neighbouring hydrogen atoms over shorter distances (on average) within 2-MI⋯H2O than within imid⋯H2O or N-MI⋯H2O. These shorter interaction distances lead to ∠(O–Hb⋯N3) = 166.3(28)° and hence enhanced non-linearity of the primary hydrogen bond in 2-MI⋯H2O relative to either imid⋯H2O or N-MI⋯H2O. The hydrogen bond7 in isoxazole⋯H2O was reported to be significantly longer, and less linear, than those identified in imid⋯H2O, N-MI⋯H2O and 2-MI⋯H2O (Table 4).
Barriers to internal rotation have been extensively explored for CH3 substituents of aromatic rings. Where there is two-fold symmetry about the axis of a CH3 rotor, electronic overlap between a π-like orbital on the CH3 and π-orbitals on the ring leads to a very low barrier to internal rotation.43 The magnitude of the barrier increases with the asymmetry of participating molecular orbitals such that it is significantly higher for 4- or 5-methylimidazole than for 2-methylimidazole.9 The same trend is apparent in published results for the 2-, 4- and 5- isomers of methylthiazole30–32 and methyloxazole.33 The value for V3 of the N-methylimidazole monomer9 (determined using XIAM) is 185.104(11) cm−1 which is very similar to the results of 182.21(12) cm−1 and 173.6(16) cm−1 for N-MI⋯H216O determined through fits performed using XIAM and BELGI-Cs respectively. The water molecule is on the opposite side of the molecule so does not interact through space with the CH3 group within N-MI⋯H2O. The evidence of the experiments reported herein is that the attachment of H2O to N-MI does not significantly change the electronic environment at the CH3 group either. The results for N-MI⋯H2O are thus consistent with previous studies where (V3) barriers to internal rotation of CH3 groups were not observed to change in response to the attachment of H2O to a remote site of the molecule.
The V3 barrier determined for 2-MI⋯H216O (determined to be 154.99(8) cm−1 by XIAM and 150.68(90) cm−1 by BELGI-Cs) is significantly greater than the value determined for the 2-MI monomer9 by XIAM (122.7529(38) cm−1). An important question is whether this difference arises primarily because of electronic or steric effects. Does the attachment of H2O induce changes in the symmetry of molecular orbitals on the imidazole ring? Alternatively, is the presence of a weak electrostatic interaction between the oxygen atom and the CH3 group of 2-MI sufficient to explain the change in V3 that follows the attachment of H2O? It is impossible to exclude a contribution for electronic effects but these did not lead to a significant change in V3 on attachment of H2O to N-MI. It is therefore more likely that the weak electrostatic interaction between the oxygen atom and the CH3 group is the predominant cause of the difference between the V3 determined for 2-MI and 2-MI⋯H2O. It is noted that an earlier analysis44 of the microwave spectrum of the toluene⋯SO2 complex also observed a large increase (from 4.9 cm−1 in free toluene to 87.3 cm−1 in toluene⋯SO2) in the barrier to internal rotation of the CH3 group on formation of that complex. The authors of that study modelled electrostatic interactions between SO2 and the CH3 group and found these to be sufficient to explain the observed change in the barrier. The authors did not exclude the possibility that electronic effects also contribute.
Footnote |
† Electronic supplementary information (ESI) available: Results of spectroscopic fits, tables of parameters. See DOI: https://doi.org/10.1039/d1cp05526g |
This journal is © the Owner Societies 2022 |