High-Resolution Rotational Spectroscopy of a Cyclic Ether - American

Jan 27, 2012 - resolution rotational study of a crown ether, 1,4,7,10,13-pentaoxacyclopentadecane. (15-crown-5), by means of molecular beam Fourier ...
0 downloads 0 Views 816KB Size
Letter pubs.acs.org/JPCL

High-Resolution Rotational Spectroscopy of a Cyclic Ether F. Gámez,† B. Martínez-Haya,*,† S. Blanco,‡ J. C. López,*,‡ and J. L. Alonso‡ †

Department of Physical, Chemical and Natural Systems, Universidad Pablo de Olavide, 41013 Seville, Spain Grupo de Espectroscopía Molecular, Edificio Quifima, Á rea de Química Física, Campus Miguel Delibes, Universidad de Valladolid, 47011 Valladolid, Spain



S Supporting Information *

ABSTRACT: The conformational landscape of crown ethers has constituted a central topic in the development of host−guest supramolecular chemistry. We report a highresolution rotational study of a crown ether, 1,4,7,10,13-pentaoxacyclopentadecane (15-crown-5), by means of molecular beam Fourier transform microwave spectroscopy. The considerable size and the broad range of conformations allowed by the flexibility of the cyclic backbone of this ether pose important challenges to spectroscopy approaches. In this investigation, three stable rotamers of the 15-crown-5 ether have been identified and characterized through their rotational constants and centrifugal distortion coefficients. Ab initio quantum calculations at the MP2 level predict these conformers as the most stable ones for the title system and reproduce accurately their distinct structural features. The results pave the ground for an extensive survey of the conformational landscape of the 15-crown-5 and related cyclic ethers in the near term. SECTION: Kinetics, Spectroscopy

C

conformers are internally cooled, but the conformer distribution before the expansion is roughly preserved, provided that interconversion barriers between conformers are sufficiently high.10 Under these conditions, the spectroscopic properties of the different conformers can be probed in a local environment of virtual isolation. Previous studies on neutral crown ether systems in supersonic expansions have explored the vibronic spectra of benzo-15-crown-511 and dibenzo-18-crown-612 and of their complexes with solvent molecules,13,14 by means of laserinduced fluorescence and double-resonance laser techniques. We report here what to our knowledge constitutes the first high-resolution rotational study of a crown ether. Microwave spectroscopy is applied in combination with quantum calculations to characterize the most stable rotamers of the 15-crown-5 molecule (henceforth, 15c5; systematic name: 1,4,7,10,13-pentaoxacyclopentadecane). The different conformational species have unique moments of inertia and give rise to separate rotational spectra. Even in conformationally challenging systems as the title one they can be distinguished by rotational spectroscopy. This has been recently illustrated for amino acids15 and with the assignment of nine low-energy conformers in γ-aminobutyric acid.16 The molecular-beam Fourier transform microwave spectrometer (MB-FTMW) employed in this study has been described elsewhere.17 A supersonic jet of 15c5 in a Ne seeding bath is generated by vaporizing the crown ether at 120 °C. The supersonic expansion relaxes the internal energy of the 15c5

rown ethers stand as one cornerstone molecular class in host−guest supramolecular chemistry and are the building blocks of a broad range of modern materials in separation and catalytic science.1,2 The most remarkable feature of crown ethers is the ability of their flexible backbone of adopting a broad range of open (ring-like) and folded (cagelike) coordination structures with metal and molecular cations.3−7 The stability of the cationic complexes of crown ethers is sustained by the net negative charge conferred to the crown cavity by the lone-pair orbitals of the ether oxygens. This is in marked contrast with the structures displayed by uncomplexed crown ethers, in which part of the oxygen sites are oriented outward from the ring.8 In this way, the mutual electrostatic repulsions between the oxygens are minimized, and weak intramolecular C−H···O bonds are allowed to form across the crown cavity. Hence, the stabilization of the cationic complex involves a “conformational switch” of the crown ether backbone between two qualitatively different structures wellseparated in conformational space: concerted turns of several dihedral angles of the initially free crown ether are required to build an optimum coordination arrangement of the oxygen sites around the cation. The barriers between the free and complexed conformations can be expected to drive the kinetics of the cation binding process. Given the ability of crown ethers to form molecular complexes easily, it is advantageous to characterize their intrinsic preferred structures under isolation conditions. A low internal temperature of the molecules is also essential because the flexibility of crown ethers leads to a distribution of their equilibrium population among a wealth of low-energy conformers.8 The preparation of molecules in supersonic jets fulfills both of these requirements.9 In the jet, the individual © 2012 American Chemical Society

Received: December 22, 2011 Accepted: January 27, 2012 Published: January 27, 2012 482

dx.doi.org/10.1021/jz201678w | J. Phys. Chem. Lett. 2012, 3, 482−485

The Journal of Physical Chemistry Letters

Letter

used for the prolate top rotamers 1 and 3, whereas the similar reduction in the IIIl representation was employed for the oblate top rotamer 2. For this latter rotamer, the centrifugal distortion constants δJ and δK were found to be nondeterminable and were kept fixed to zero. The fitted transitions match the experimental lines within a standard deviation of σFIT < 2 kHz. The best-fit rotational constants and centrifugal distortion coefficients of the three rotamers identified in this study are given in Table 1. It must be remarked that the small values of the rotational constants that characterize the 15c5 rotamers constitute an important challenge for their spectroscopic study in the gigahertz range. The observed rotamers of 15c5 have been assigned to the three lowest energy conformers predicted by the ab initio computations. These conformers are represented in Figure 2 and are described through their rotational parameters and dipole moments in Table 1. (Further structural parameters, interatomic distances, and bond angles are provided in Table S4 as Supporting Information.) The principal rotational constants of the computational conformers 1, 2, and 3 (labeled in order of increasing free energy) resemble closely those of the three rotamers characterized in the experiments. None of the other MP2 conformers found in our study matched the rotational constants of the experimental rotamers with a comparable level of accuracy. Quite good agreement exists as well between the experimental and MP2 values for the quartic centrifugal distortion coefficients. Also importantly, the components of the dipole moment of each conformer are consistent with the type and intensities of the lines observed experimentally. For rotamers 1 and 3, strong μb-type rotational transitions are detected, in agreement with a dipole moment aligned along the b principal axis. In addition, weaker transitions of μa type for rotamer 1 and of μc type for rotamer 3 could be assigned, again in consonance with the relative magnitude of the dipole moment components. Finally, rotamer 2 displayed strong μatype transitions, in accord with a dominant component of the dipole moment along the a axis. The MP2 calculation predicts conformer 1 to have a Gibbs free energy lower by only 1.1 and 1.6 kJ mol−1 (87 and 137 cm−1) than the next higher conformers 2 and 3, respectively. These conformers are illustrative of the remarkable diversity of energetically similar structures that the 15c5 ether can adopt. Conformers 1 and 2 are coincident with the conformers 1 and 11, respectively, reported by Hill and Feller.8 Conformer 3 had not been described in any previous work. The particularly small value for the rotational constant C for the three conformers is related to open ring conformations of the ether backbone, providing a relatively large moment of inertia with respect to the principal axis perpendicular to the ring. The two prolate conformers 1 and 3 have similar features, as reflected in the values of the experimental rotational constants. Each of these conformers displays significant differences between its three planar inertial moments. This is indicative of different mass distributions out of each of the inertial planes. The small values of the planar inertial moment Pcc = (Ia + Ib − Ic)/2 = ∑imici2 (mi and ci are the mass and c coordinate of atom i), which gives the mass extension out of the ab inertial plane, indicate that the heavy atoms (O and C) lay at positions not far from that plane. Similarly, the comparably large values found for Paa indicate that the mass distribution is extended mainly out of the bc plane. Conformers 1 and 3 also have in common a roughly coplanar arrangement of the O4−C−C−O5 segment of the ring (see atom numbering in Figure 2). Such arrangement

molecules and carries them into the 55 cm confocal spherical mirror Fabry−Perot resonator. A short microwave pulse (0.3 μs, 10−300 mW) polarizes the species in the jet, and the subsequent coherent transient emission is recorded in the time domain (typically over ca. 400 μs at 40−100 ns intervals) and Fourier-transformed to the frequency domain. The frequency accuracy of the experiment is better than 5 kHz. The computational study involved an extensive survey of the conformational space of the isolated 15c5 by means of simulated annealing. The resulting low-energy molecular structures were optimized with the ab initio second-order Møller−Plesset (MP2) method and the basis set 6-311+ +G(d,p). As many as 15 stable conformers are predicted with MP2 free-energy differences within 8 kJ mol−1, corresponding to asymmetric rotors with a broad range of positive (oblate) and negative (prolate) values of the asymmetry parameter κ = (2B − A − C)/(A − C) and a variety of electric dipole moment values. Wide frequency scans were conducted to search for K−1 and K+1 doublet patterns that provide identifiable signatures for individual conformers. (See Figure 1 for examples of rotational

Figure 1. Typical J′K−1 ′ ,K+1 ′ → J″K−1 ″ ,K+1 ″ rotational lines measured for 15crown-5 in this work. The J′ = 12 ← J″ = 11 and J′ = 5←J″ = 4 progressions are assigned to rotamers 1 and 3, respectively. Each line displays a Doppler splitting of ca. 40 kHz. The absolute centerfrequency of each transition is determined from the half-sum of the two Doppler components with an accuracy better than 5 kHz.

lines illustrative of the overall spectral resolution.) In this way, the spectra of three different rotameric structures of the 15c5 crown ether, labeled 1, 2, and 3, were identified. Ensembles of several tens of lines were measured for each rotamer within the frequency range 4−10 GHz, probing rotational states J = 4−16. A list of the complete set of rotational transitions observed is provided as Supporting Information. The assignment of the experimental lines to rotational transitions of a specific conformer was performed through an iterative measurement/ fit procedure based on Watson’s semirigid rotor Hamiltonian18−20 up to quartic terms. In the final fits, the asymmetric top reduction of the Hamiltonian in the Ir representation was 483

dx.doi.org/10.1021/jz201678w | J. Phys. Chem. Lett. 2012, 3, 482−485

The Journal of Physical Chemistry Letters

Letter

Table 1. Effective Rotational Constants (A, B, C), Anisotropy Parameter (κ = (2B − A − C)/(A − C)), Planar Inertial Moments (Paa = (Ib + Ic − Ia)/2, etc.) and Quartic Centrifugal Distortion Coefficients (ΔJ, ΔJK, ΔK, δJ, δK) in the A-Reduced Ir Representation for Conformers 1 and 3 (Prolate Asymmetric Tops) and in the A-Reduced IIIl Representation for Conformer 2 (Oblate Asymmetric Top)a conformer 1 exp A/MHz B/MHz C/MHz κ Paa/(amu Å2) Pbb/(amu Å2) Pcc/(amu Å2) ΔJ/kHz ΔJK/kHz ΔK/kHz δJ/kHz δK/kHz |μa|/D |μb|/D |μc|/D number of lines σFIT/kHz

b

806.82516(19) 417.16283(22) 303.779894(42) −0.549210(1) 1124.36153(51) 539.27424(51) 87.10570(51) 0.0103(10) 0.0608(48) 0.1212(43) 0.00147(47) 0.0534(58)

conformer 3

conformer 2

MP2

exp

MP2

exp

MP2

828.5 416.9 307.7 −0.582 1122.3 520.1 89.9 0.0098 0.0675 0.112 0.0064 0.0522 0.9 2.3 0.6

710.44967(14) 434.54680(15) 287.053113(37) −0.303280(2) 1106.11432(38) 654.46265(38) 56.88837(38) 0.01017(571) 0.0435(25) 0.0731(318) 0.00352(28) 0.0426(21)

725.3 433.2 289.2 −0.339 1108.7 638.8 57.9 0.0100 0.0486 0.066 0.0034 0.046 0.1 4.2 1.6

595.89572(45) 560.83329(32) 342.00682(13) 0.723797(5) 765.35204(86) 712.33479(86) 135.76814(86) 0.0820(37) −0.1345(90) 0.0663(51) [0.0]c [0.0]

602.9 565.6 347.9 0.710 754.0 698.7 139.6 0.0883 −0.146 0.0715 −0.010 0.129 3.7 0.1 0.8

55 1.3

61 1.3

46 1.7

a Number of experimental lines and the standard deviation, σFIT, of the fits are indicated in the bottom entries. The theoretical rotational constants, centrifugal distortion coefficients, and components of the dipole moment (μa, μb, μc) of the MP2/6-311++G(d,p) conformers are indicated next to the values of the corresponding experimental rotamers. bStandard error in parentheses in units of the last digit. cValues in square brackets were kept fixed in the fit.

the value of Pcc is larger than those of the two prolate conformers. This renders a unique compact structure that features by far the largest B and C rotational constants and κ parameter among the lowest 15 conformers predicted by the MP2 computations. Such structure provides a homogeneous ensemble of distances between the oxygens (2.9 to 3.1 Å) and allows for efficient C−H···O interactions. It is interesting to note that ΔJ for conformer 2 is significantly higher than it is for the other two conformers. This result is found irrespective of the representation and type of Watson’s reduction employed to analyze the complex; for example, the analysis in the Ir representation gives a value of ΔJ of 0.0362(55) kHz (ab initio 0,032 kHz), which is still more than three times higher than found for conformers 1 and 3 in the same representation. It is interesting to evaluate the magnitude of the conformational “switch” that the free 15c5 ether structures described above must undergo to form their typical cationic complexes. The open structures of conformers 1 and 3 of the free 15c5 ether resemble to some extent the features of the 15c5 backbone in the complexes with cations of relative large size, such as the heavier alkali metals4 or with the molecular cations HO3+ and NH4+.5 Conformer 2 would be a priori closer to the folded structures adopted by the native crown ethers in their complexes with cations smaller than the cavity.3,21 Nevertheless, significant rearrangement of the ether backbone is still required to bring all oxygens to the interior of the cavity when forming the complex. In particular, cation binding leads to a significant shortening of the distances between consecutive oxygens in the ring. For instance, the 15c5−K+ complex displays an average O−O distance of 2.85 Å,4 in comparison with 3.0 to 3.1 Å for the conformers of the isolated 15c5 presently described. Crown ethers can display a variety of qualitatively different conformations within a narrow energetic range. Each of such

Figure 2. Representation of the MP2/6-311++G(d,p) conformers of the isolated 15-crown-5 ether, which are unequivocally identified as the rotamers 1, 2, and 3 characterized in the present microwave rotational spectroscopy experiments. Three orientations according to the principal axis (a, b, c) of the conformers are shown. Some relevant structure parameters are provided in Table 1 and in Table S4 of the Supporting Information.

orients one of the oxygens toward the outside of the ring cavity and places it at a particularly large distance from the other oxygen sites. A detailed inspection of these conformers also reveals differences in the location of the polar oxygen atoms along the ether backbone in each case, which results in noticeable differences in the values of the electric dipole moments. Conformer 2 is qualitatively different from conformers 1 and 3 and corresponds to an oblate asymmetric top. This conformer is characterized by a remarkable degree of folding of the cyclic backbone, as can be appreciated in Figure 2. Correspondingly, 484

dx.doi.org/10.1021/jz201678w | J. Phys. Chem. Lett. 2012, 3, 482−485

The Journal of Physical Chemistry Letters

Letter

conformations finds an optimum balance between the flexibility of the ether backbone and the intramolecular O···O repulsions and CH···O attractions. This study has shown that MB-FTMW spectroscopy and quantum MP2 methods constitute powerful tools to unveil the fine aspects of the conformational landscape of these systems. In the present experiments, there was an ensemble of weaker unassigned lines remaining in the spectra that are likely to belong to conformers of the 15-crown-5 ether lying higher in energy than the three ones characterized here. Current work is being devoted to the identification of such conformers. This task will involve a considerable extension of the spectroscopic survey and will be the topic of a comprehensive report in the near future.



(10) Ruoff, R. S.; Klots, T. D.; Emilsson, T.; Gutowsky, H. S. Relaxation of Conformers and Isomers in Seeded Supersonic Jets of Inert Gases. J. Chem. Phys. 1990, 93, 3142−3150. (11) Shubert, V. A.; James, W. H., III; Zwier, T. S. Jet-Cooled Electronic and Vibrational Spectroscopy of Crown Ethers: Benzo-15crown-5 Ether and 4′-Amino-benzo-15-crown-5 Ether. J. Phys. Chem. A 2009, 113, 8055−8066. (12) Kusaka, R.; Inokuchi, Y.; Ebata, T. Structure of Hydrated Clusters of Dibenzo-18-crown-6-ether in a Supersonic Jet-encapsulation of Water Molecules in the Crown Cavity. Phys. Chem. Chem. Phys. 2008, 10, 6238−6244. (13) Shubert, V. A.; Müller, C. W.; Zwier, T. S. Water’s Role in Reshaping a Macrocycle’s Binding Pocket: Infrared and Ultraviolet Spectroscopy of Benzo-15-crown-5-(H2O)n and 4′-Aminobenzo-15crown-5-(H2O)n, n = 1,2. J. Phys. Chem. A 2009, 113, 8067−8079. (14) Kusaka, R.; Kokubu, S.; Inokuchi, Y.; Haino, T.; Ebata, T. Structure of Host−Guest Complexes between Dibenzo-18-crown-6 and Water, Ammonia, Methanol, and Acetylene: Evidence of Molecular Recognition on the Complexation. Phys. Chem. Chem. Phys. 2011, 13, 6827−6836. (15) Alonso, J. L.; Perez, C.; Sanz, M. E.; López, J. C.; Blanco, S. Seven Structures of Threonine in the Gas Phase: A LA-MB-FTMW Study. Phys. Chem. Chem. Phys. 2009, 11, 617−627. (16) Blanco, S.; López, J. C.; Mata, S.; Alonso, J. L. Conformations of γ-Aminobutyric Acid (GABA): The Role of the n → π* Interaction. Angew. Chem., Int. Ed. 2010, 49, 9187−9192. (17) Alonso, J. L.; Lorenzo, F. J.; López, J. C.; Lesarri, A.; Mata, S.; Dreizler, H. Construction of a Molecular Beam Fourier Transform Microwave Spectrometer Used to Study the 2,5-Dihydrofuran-argon van der Waals Complex. Chem. Phys. 1997, 218, 267−275. (18) Watson, J. K. G. Aspects of Quartic and Sextic Centrifugal Effects on Rotational Energy Levels. In Vibrational Spectra and Structure: A Series of Advances; Durig, J. R., Ed.; Elsevier: New York, 1977; pp 1−89, Vol. 6. (19) Watson, J. K. G. Determination of Centrifugal Distortion Coefficients of Asymmetric-Top Molecules. J. Chem. Phys. 1967, 46, 1935−1949. (20) Kivelson, D.; Wilson, E. B. Approximate Treatment of the Effect of Centrifugal Distortion on the Rotational Energy Levels of Asymmetric-Rotor Molecules. J. Chem. Phys. 1952, 20, 1575−1579. (21) Gámez, F.; Hurtado, P.; Martínez-Haya, B.; Berden, G.; Oomens, J. Vibrational Study of Isolated 18-crown-6 Ether Complexes with Alkaline-Earth Metal Cations. Int. J. Mass Spectrom. 2011, 308, 217224.

ASSOCIATED CONTENT

S Supporting Information *

Complete list of the experimental lines measured for the three 15c5 rotamers characterized in this work and selected geometrical parameters of the MP2 conformers. This material is available free of charge via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected], [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The research leading to these results is supported by the Programme Consolider−Ingenio of the Ministry of Science and Innovation of Spain (MICINN) through project CSD200900038. B.M.H. acknowledges funding through projects P07FQM-02600 (Junta de Andalucia-FEDER) and CTQ200910477 (MICINN).



REFERENCES

(1) Pedersen, C. J. The Discovery of Crown Ethers. Science 1988, 241, 536−540. (2) Gokel, G. W. Crown Ethers and Cryptands; Royal Society of Chemistry: Cambridge, U.K., 1994. (3) Martínez-Haya, B.; Hurtado, P.; Hortal, A. R.; Hamad, S.; Steill, J. D.; Oomens, J. Emergence of Symmetry and Chirality in Crown Ether Complexes with Alkali Metal Cations. J. Phys. Chem. A 2010, 114, 7048−7054. (4) Hurtado, P.; Hortal, A. R.; Gámez, F.; Hamad, S.; Martínez-Haya, B. Gas-Phase Complexes of Cyclic and Linear Polyethers with Alkali Cations. Phys. Chem. Chem. Phys. 2010, 12, 13752−13758. (5) Hurtado, P.; Gámez, F.; Hamad, S.; Martínez-Haya, B.; Steill, J. D.; Oomens, J. Crown Ether Complexes with H3O+ and NH4+ Proton Localization and Proton Bridge Formation. J. Phys. Chem. A 2011, 115, 7275−7282. (6) Rodriguez, J. D.; Lisy, J. M. Probing Ionophore Selectivity in Argon-Tagged Hydrated Alkali Metal Ion−Crown Ether Systems. J. Am. Chem. Soc. 2011, 133, 11136−11146. (7) Inokuchi, Y.; Boyarkin, O. V.; Kusaka, R.; Haino, T.; Ebata, T.; Rizzo, T. R. UV and IR Spectroscopic Studies of Cold Alkali Metal Ion−Crown Ether Complexes in the Gas Phase. J. Am. Chem. Soc. 2011, 133, 12256−12263. (8) Hill, S. E.; Feller, D. Theoretical Study of Cation/Ether Complexes: 15-crown-5 and its Alkali Metal Complexes. Int. J. Mass Spectrom. 2000, 201, 41−58. (9) Atomic, Molecular Beam Methods; Scoles, G., Ed.; Oxford University: Oxford, U.K., 1992; Vol 2 485

dx.doi.org/10.1021/jz201678w | J. Phys. Chem. Lett. 2012, 3, 482−485