Huaqi Zhan,
Yongjun Hu*,
Pengchao Wang and
Jiaxin Chen
MOE Key Laboratory of Laser Life Science & Institute of Laser Life Science, College of Biophotonics, South China Normal University, Guangzhou 510631, P. R. China. E-mail: yjhu@scnu.edu.cn
First published on 18th January 2017
Morpholine (NH(CH2CH2)2O) is a typical six-membered aliphatic heterocyclic compound. Herein, infrared plus vacuum ultraviolet (IR/VUV) single photon “soft” ionization spectroscopy was employed to study the structures of neutral morpholine and its monohydrated clusters. Theoretical calculations revealed that the structures containing equatorial-chair and axial-chair conformations were the most stable conformers in the gas phase, and this was confirmed by IR spectral analysis. Analysis of the observed and calculated spectra of the monohydrated clusters suggested that multiple conformers may co-exist in the molecular beam, and that the water molecule acts as a hydrogen donor. In the most stable structure, the hydrogen atom of the water molecule is bound to the NH group of the equatorial-chair conformer of morpholine. Moreover, the water molecules simultaneously serve as hydrogen bond donors for the NH group and interact with the CH group weakly. It is suggested that the weak intermolecular CH⋯O interaction is also responsible for the molecular stability.
Morpholine is a typical six-membered aliphatic heterocyclic compound with the molecular formula NH(CH2CH2)2O. It contains a heterocyclic ring that is present in many compounds of biological and pharmaceutical relevance.9 Most interestingly, the morpholine heterocycle contains both nitrogen and oxygen atoms. The carbon atoms adjacent to the N–H bond are saturated, ensuring sp3 hybridization at the nitrogen center, with the result that the N atom lies out of the plane defined by the four carbon atoms. This gives rise to two possible conformers where the N–H bond lies either axial or equatorial.10 We were also curious to see what would happen if a monomer provided two different acceptor sites. Morpholine is an amine as well as an ether, and each group can provide an electron pair to accept a water proton and form a hydrogen bond. Because the alkalinity of the amine group is considerably greater than that of the ether group, we expected that the N⋯H–O bonded complex would be preferred.
In general, the molecular structure of morpholine can be explored using several experimental techniques, such as NMR, IR, microwave, and Raman spectroscopy, and dipole moment and Kerr effect measurements, etc.11–14 However, different conclusions have been reached by different groups. For instance, the axial-chair conformation of morpholine was deduced from dipole moment and Kerr effect measurements. Allinger et al. and Baldock and Katritzky,15,16 suggested an equatorial-chair conformation. Smith and Shoulders17 could not distinguish the conformers that were present in morpholine using NMR measurements because the axial and equatorial conformations easily interconvert on the timescale of NMR spectroscopy. Baldock and Katritzky16 considered the equatorial positions to be more stable from IR measurements between the equatorial and axial conformation. Sloan and Kewley13 found many characteristic features for the equatorial-chair conformation from microwave spectroscopy but none for the axial-chair form. In 1989, Grabow et al.,14 repeated the microwave spectroscopy measurements for the normal and N-deuterated species at very high resolution and derived deuterium quadrupole coupling constants. The presence of two conformers was supported by Raman spectroscopy in the liquid state.18,19 However, the details of the structure of the morpholine molecule and its hydrated clusters in the gas phase need further studies.
It is well know that Fourier transform IR (FTIR) absorption spectroscopy is one of the most effective approaches to study molecule and cluster structures, which had made massive outstanding achievements,20–31 these include the groups of Suhm, Georges, and Asselin. In this paper, we introduce an alternative infrared absorption spectroscopy technique. Supersonic jet-laser spectroscopy is one of the most powerful tools for investigating the variety of conformers because it enables researchers to freeze conformational conversion using the ultracold, collision-free conditions in a jet.32 The frozen conformers are easily distinguished by IR dip spectroscopy, and their structures can be precisely determined. Recently, mass-selected optical spectroscopy has been widely applied in the investigation of gas-phase structures and dynamics. Combining mass-selected optical spectroscopy with the supersonic jet technique enables us to determine the geometrical structures and dynamics of molecules and clusters with high reliability. Vibrational spectroscopy in supersonic jets is an effective approach that can provide information on the geometric structures of molecules and clusters, and their intermolecular interactions.33–38 Furthermore, the vibrational frequencies and intensities of IR absorption modes, such as modes involving OH, NH, and CH groups, are generally very sensitive probes of hydrogen bonding and other intra- and intermolecular interactions. This makes them excellent probes for the investigation of detailed aspects of the structures of biologically relevant ions, radicals, and molecules.39 In addition, the extent of hydrogen bonding can be directly investigated by using IR spectroscopy to detect the OH and NH stretching modes of the relevant clusters.40 The OH and NH vibrations are very susceptible to local perturbations and show dramatic red shifts when the OH or NH group is strongly hydrogen bonded as an acceptor group, such as the OH group of an adjacent water or alcohol molecule.41–44
Modern quantum chemical calculations, for instance, density functional theory (DFT) and Møller Plesset second-order perturbation theoretical functional (MP2) calculations, can be used to optimize the structures and predict the energies and frequencies of monomer molecules and clusters.19,45,46 Recently, vibrational spectroscopy combined with quantum chemical calculations was successfully applied to distinguish molecular conformers, exploring conformational equilibrium in the condensed state and assign vibrational bands.32,47,48
In the present study, we combine IR spectroscopy with various theoretical methods to investigate the structure of morpholine and its monohydrated complexes. These structures were probed using IR spectroscopy with vacuum ultraviolet (VUV) photoionization detection.
Generation of the VUV 118 nm laser light and the IR light was reported in previous publications.50 The VUV radiation is generated as the ninth harmonic of the fundamental of the 1064 nm output of a Nd3+:YAG (yttrium aluminum garnet) laser. The fundamental 1064 nm is tripled to 355 nm, which is focused into a cell with Xe/Ar at a ratio 1:10 at approximately 250 Torr by a 25 cm lens. The 118 nm light generated is focused by a MgF2 lens into the ionization region of a TOFMS. This same lens defocuses the residual 355 nm light. The 118 nm light is about 5 μJ per pulse or 1012 photons per pulse53 and the conversion efficiency is roughly 1.2 × 10−5. This light is employed for single photon ionization of the molecules and clusters in the molecular beam. Tunable IR radiation 3–5 mJ per pulse in the 2800–3800 cm−1 range with a bandwidth of 2–3 cm−1 (broad band mode) is generated from an optical parametric oscillator/optical parametric amplifier (OPO/OPA) system (LaserVision) pumped by the Nd:YAG fundamental and counterpropagates with respect to the VUV beam. Both the VUV and IR laser beams are perpendicular to the molecular beam and to the ion flight path of the TOFMS.
Conformer | Erel (B3LYP) | Erel (MP2) | Erel (ωB97X-D) |
---|---|---|---|
Eq-chair | 0 | 0 | 0 |
Ax-chair | 4.3 | 4.8 | 3.1 |
Ax-boat | 26 | 26 | 25 |
Eq-boat | 30 | 35 | 32 |
Ma–H2O-1 | 0 | 0 | 0 |
Ma–H2O-2 | 6.4 | 4.8 | 6.1 |
Ma–H2O-3 | 10 | 7.7 | 12 |
Ma–H2O-4 | 25 | 19 | 22 |
Me–H2O-1 | 8.5 | 5.8 | 6.2 |
Me–H2O-2 | 12 | 9.6 | 10 |
Me–H2O-3 | 15 | 13 | 14 |
Me–H2O-4 | 30 | 23 | 24 |
Fig. 2 The four lowest-energy structures of morpholine calculated by ωB97X-D/cc-pVTZ. The relative energy in kJ mol−1 is listed above each structure. |
The observed and calculated IR spectra of the neutral morpholine monomer in the range of 2600–3600 cm−1 are shown in Fig. 3; this frequency range corresponds primarily to C–H and N–H stretching. The calculated spectra are based on the four conformations of neutral morpholine at the ωB97X-D/cc-pVTZ level, and the vibrational frequencies are scaled by 0.9415. Fig. 3(a) shows the observed IR spectrum of gas-phase neutral morpholine in the range of 2600–3150 cm−1; the C–H stretching bands are clearly observed as slightly broadened bands at about 2820, 2848, 2907, 2939, and 2954 cm−1. Theoretical calculations show that both equatorial-chair and axial-chair are the stablest structures among all of the conformers, and the spectrum calculation for the two chair conformers show good agreement with all of the observed spectral features in the studied range. In contrast, the calculated spectra for the other conformers do not reproduce the observed intense C–H stretching bands very well. The IR spectrum of neutral morpholine in the range of 3200–3500 cm−1 is also show in Fig. 3(a). The frequencies of this region correspond primarily to N–H stretching vibrations. The bands attributed to the N–H stretching of the axial and equatorial conformers are observed at 3320 and 3336 cm−1, respectively, and these two bands overlap at approximately 3330 cm−1. Theoretical calculation results in this region show good agreement with the observed spectra (Fig. 3). Generally, the results above reveal that both of the two chair conformers may coexist in the conformational equilibrium of the gas. Hence, only the two chair conformers are considered further in this paper. Furthermore, the theoretical calculations reveal that the equatorial-chair conformer is more stable, and thus it may be more favorable in the gas phase.
Fig. 3 Observed IR spectrum of neutral morpholine and the calculated spectrum based on the optimized structure shown (ωB97X-D/cc-pVTZ). The calculated frequencies are scaled by 0.9415. |
Fig. 4 The eight lowest-energy structures of morpholine monohydrated complexes calculated by ωB97X-D/cc-pVTZ. The relative energy in kJ mol−1 is listed above each structure. |
In the most stable structure of monohydrated morpholine, Me–H2O-1 (shown in Fig. 4), the morpholine molecule has an equatorial-chair conformation and the water molecule serves as a hydrogen bond donor for the NH group. The morpholine molecule also exists as the equatorial-chair conformer in Me–H2O-2 and Me–H2O-3, but the water molecule serves as a hydrogen bond donor for the ether oxygen. The calculated results reveal that the energies of Me–H2O-2 and Me–H2O-3 are about 6.1 and 12 kJ mol−1 higher than that of Me–H2O-1 at the ωB97X-D/cc-pVTZ level, respectively. The theoretical calculations presented in Fig. 4 thus show that the monohydrated morpholine complex with the equatorial-chair conformation is more stable. Furthermore, the most stable monohydrated clusters containing equatorial and axial morpholine conformers are Me–H2O-1 and Ma–H2O-1, respectively. This indicates that the water molecule preferentially serves as a hydrogen donor for the NH group in morpholine. This conclusion agrees well with the Fourier transform microwave study of Indris et al.,57 and it reflects the large intrinsic basicity of this functional group, which is much larger than that of the ether oxygen. However, the water molecule serves as a hydrogen acceptor in the least stable complexes (Me–H2O-4 and Ma–H2O-4). Therefore, there is only a low probability that the morpholine molecule serves as hydrogen donor in the monohydrated morpholine.
Fig. 5 presents the IR plus VUV spectrum recorded at 118 nm for monohydrated morpholine, together with the corresponding spectra calculated for the eight lowest-energy conformers of monohydrated morpholine. These spectra are based on the conformations optimized at the ωB97X-D/cc-pVTZ level of theory, and the vibrational frequencies were scaled by 0.9415. In the observed spectrum of monohydrated morpholine (Fig. 5(a)), a sharp intense dip peak is observed at 3714 cm−1, which corresponds to the free O–H stretch of water. The two broad bands at about 3309 and 3432 cm−1 must therefore result from hydrogen bonded. The feature at 3432 cm−1 can be assigned to O–H stretching hydrogen bond, and the feature at 3309 cm−1 to the N–H stretching hydrogen bond. These spectral features seem to indicate that the monohydrated morpholine has a ring-type structure, in which the water acts as a bridge and serves as hydrogen bond acceptor to the N–H group and a hydrogen bond donor to the ether oxygen. However, DFT and MP2 energy calculations do not support this assignment, and the geometry optimization shows that this is an unstable structure. As we can see in Fig. 5, the observed spectrum does not match very well with the simulated spectra of the monohydrated morpholine complexes that contain O–H⋯N (b, e) or O–H⋯O hydrogen bonds (c, d, f, g). However, the signals corresponding to O–H⋯N hydrogen bond stretching in simulated spectra (b) and (e) match with the signal in the observed spectrum at 3309 cm−1, and those corresponding to O–H⋯O hydrogen bond stretching in (c), (d), (f), and (g) match with the signal in the observed spectrum at 3432 cm−1. Therefore, we conclude that structures of monohydrated morpholine in which the water molecule behaves as a hydrogen bond donor towards the NH group and the ether oxygen co-exist in the molecular beam. The discrepancy in the calculated spectra may be explained by the anharmonic resonance with the combination and the overtone bands of lower frequency modes. In addition, the absence of an antisymmetric free OH stretch in the experimental spectrum indicates that the morpholine monohydrated clusters in which the morpholine molecule behaves as the hydrogen donor do not contribute to the IR spectrum. This is also supported by the theoretical calculations. In other words, these results imply that more than one monohydrated conformer is present in the beam, where it would be more appropriate for the morpholine molecule to behave as the hydrogen donor in the monohydrated clusters. Considering the cooling effect in the molecular beam, the low energy structures will be dominant at low temperature. Indris and coworkers employed FT microwave spectroscopy to investigate the morpholine–H2O complex.57 Their results demonstrate that the N⋯H–O hydrogen-bonded structure was found existing in a molecular beam, while the O⋯H–O hydrogen-bonded structure was not detected. It is worth noting that no sharp band related to free N–H stretching is seen in Fig. 5; the absence of this feature may be due to the fact that the N–H stretching signal is too weak and its position may overlap with the broad hydrogen bond signals.
Fig. 5 Observed IR spectrum of neutral monohydrated morpholine and the calculated spectrum based on the optimized structure shown (ωB97X-D/cc-pVTZ). The calculated frequencies are scaled by 0.9415. |
We can find that the water molecules form two different kinds of hydrogen bond (O–H⋯O and O–H⋯N) in the complexes with morpholine molecule in the spectra. And the hydrogen bond length of O–H⋯N calculated in Me–H2O-1 is 189 pm, which is longer than the O–H⋯O hydrogen bond in Me–H2O-2 (185 pm) (Fig. 6). These two hydrogen bonds attributable to strong hydrogen bonds.60 To substantiated the spectroscopy results and estimated the interaction strength, we carried out the natural bond orbitals (NBO) calculations with MP2/6-311++G(d,p) (Pop = nbo) level. The second-order perturbation stabilization energy E(2) obtained from natural bond orbital (NBO) analysis for the n(O) → σ*(O–H) hydrogen bond interaction in Me–H2O-2 (48 kJ mol−1) is more than the n(N) → σ*(O–H) hydrogen bond interaction in Me–H2O-1 (29 kJ mol−1). The stabilization energy gained from the hydrogen bond interaction, denoted E(2), is a measure of the interaction strength.61 This result means that the σ* orbital of the water O–H bonds overlap better with the oxygen lone pair orbital of morpholine. As we can see in Fig. 4, the oxygen atom of morpholine acts as hydrogen bond acceptor, while the water molecule serves as a hydrogen bond donor and seems to have a weakly intermolecular interacts with the CH groups at the same time (Fig. 6) in the Me–H2O-2 and Ma–H2O-2 structure. However, this phenomenon was not found in Me–H2O-3 and Ma–H2O-3 structure. Furthermore, the Me–H2O-2 and Ma–H2O-2 structure is more stable than Me–H2O-3 and Ma–H2O-3 structure in energy respectively. This suggests that the weak intermolecular interaction C–H⋯O may be contributes to the molecular stability of the complex. The occurrence of the weak intermolecular reactions improves the relative stability of these morpholine–water complexes. To prove the existence of this interaction, the NBO diagrams shows that the oxygen lone pair orbital of water and the σ* orbital of the C–H bonds overlap in Me–H2O-1 and Me–H2O-2 structure would lead to intermolecular interaction occurring (see Fig. 6). Moreover, the second-order perturbation stabilization energies, E(2), between the nonbonding orbital and the σ* orbital (Lewis type NBO) is 0.92, 0.42 kJ mol−1 in Me–H2O-1 and Me–H2O-2 respectively. And we can found that the σ* orbitals of the C–H bonds overlap with the oxygen lone pair orbitals in Me–H2O-1 and Me–H2O-2. Thus, the interactions among the oxygen lone pair orbitals and the C–H bonds are confirmed. These results imply that there are weak interactions C–H⋯O between the molecules (morpholine and water) besides the strong hydrogen bonding interactions. And this weak force may also play an important role in the intermolecular reactions and contribute to improve the relative stability of these conformations.
The results show that the two chair conformers (axial and equatorial) are the most stable structures among all of the conformers, which are the dominant in a supersonic expansion. The predicted IR spectra for these two conformers show good agreement with the observed spectral features. DFT and MP2 calculations show that monohydrated morpholine has several stable conformers. The obtained IR spectra suggest that multiple isomers may be formed in the expansion, and that morpholine prefers to act as the hydrogen acceptor in hydrogen bonding involving the NH group and the ether oxygen. Moreover, NBO analysis result shows that the weak intermolecular reactions CH⋯O exists in the complexes, and this weak interaction contribute to improve the relative stability of these conformations.
This journal is © The Royal Society of Chemistry 2017 |