Yoshiya
Inokuchi
*a,
Takayuki
Ebata
a,
Toshiaki
Ikeda
a,
Takeharu
Haino
a,
Tetsunari
Kimura
b,
Hao
Guo
b and
Yuji
Furutani
b
aDepartment of Chemistry, Graduate School of Science, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan. E-mail: y-inokuchi@hiroshima-u.ac.jp; Tel: +81-82-424-7101
bInstitute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan
First published on 28th August 2015
We demonstrate a powerful spectroscopic technique, surface-enhanced infrared absorption (SEIRA) spectroscopy, not only for detecting host–guest complexes in solution but also for examining the relationship between the guest selectivity, complex structure, and solvent effect. We synthesize thiol derivatives of 15-crown-5 and 18-crown-6 [2-(6-mercaptohexyloxy)methyl-15-crown-5 (15C5-C1OC6-SH) and 2-(6-mercaptohexyloxy)methyl-18-crown-6 (18C6-C1OC6-SH)], which are adsorbed on gold surfaces through S–Au bonds. The IR difference spectra of the M+·15C5-C1OC6 (M = Li, Na, K, Rb, and Cs) complexes on gold are observed using aqueous solutions of MCl by SEIRA spectroscopy. The spectra show a noticeable change in the C–O stretching vibration at around 1100 cm−1. The spectral patterns of M+·15C5-C1OC6 are similar for Li+ and Na+, and for K+, Rb+, and Cs+; the interaction between the metal ions and 15C5-C1OC6 changes drastically between Na+ and K+ in the series of alkali metal ions. On the other hand, the equilibrium constant for complex formation determined by the IR intensity shows clear preference for Na+ ions. We also observe the IR difference spectra of M+·18C6-C1OC6 in methanol and compare them with those in water. The spectral patterns in methanol are almost the same as those in water, but the equilibrium constant in methanol does not show preference for any ion, different from the K+ preference in water. From these findings we attribute the origin of the ion selectivity of 15C5 and 18C6 in solution to the interaction between the metal ions and the crown ethers in the complexes or the solvation energy of free ions. In the case of 15C5-C1OC6 in water, the preference of Na+ over K+, Rb+, and Cs+ can be attributed to the strength of the interaction or the size matching between metal ions and 15C5-C1OC6; the Na+ selectivity over Li+ ions is dominated by the solvation energy of free ions. For 18C6-C1OC6 in methanol, the equilibrium constant for complex formation becomes much bigger in methanol than that in water and loses the selectivity in methanol, because the solvation energy in methanol is fairly smaller than that in water, predominating the contribution from the strength of the interaction between metal ions and 18C6-C1OC6. The IR spectra measured by SEIRA spectroscopy are quite sensitive to properties of host–guest complexes such as the intermolecular interaction, the structure, and the orientation against the gold surface. However, the evidence for guest selectivity emerges primarily in the intensity of the spectra, rather than band positions or spectral patterns in the IR spectra.
More recently, we have obtained the IR spectra of metal ion–CE complexes in solutions by surface-enhanced IR absorption (SEIRA) spectroscopy.10 In these measurements, CEs are bound to gold surfaces through S–Au chemical bonds, and solutions containing metal ions are placed on the surfaces. We found that the equilibrium constant of 18C6 tagged on the gold surfaces is the largest for K+ among alkali metal ions. This result is in good agreement with previously reported ones in solutions,1b,c showing that this method is well suitable for studying the encapsulation process in solutions. Thanks to the enhancement of the IR absorption by gold surfaces and its confinement to the immediate vicinity of the surfaces,11 the formation of metal ion–CE complexes can be detected with very high sensitivity, avoiding interference from strong absorption of solvent. The high sensitivity nature of SEIRA spectroscopy also makes it possible to measure the IR spectra with a very wide range of guest concentrations. This will provide the equilibrium constant for the complex formation from the IR spectra and reveal the relationship between the complex structure and the ion selectivity as in the case of a biological sample.11b,d Since the encapsulated metal ions can be removed and changed to other ions very easily with the CEs kept on gold surfaces, systematic and prompt measurements of the IR spectra are possible, which is essential for quantitative measurements with a series of metal ions. In this experiment, it is not necessary to dissolve CEs and metal ions in the same solvent. This enables us to study the solvent effect on encapsulation by using a variety of solvents. In addition, the tagging of CEs on gold surfaces is of great advantage for applications such as ion sensing.12
In this paper, we examine the crown size and solvent dependence on the complex formation of CEs by SEIRA spectroscopy. We synthesize thiol derivatives of 15C5 and 18C6, which are adsorbed on gold surfaces through S–Au chemical bonds. Two types of solutions, aqueous and methanol solutions of alkali metal chloride salts, are placed on the surface, and IR spectra are measured by SEIRA spectroscopy. We examine the spectral change (frequency shift and intensity change) of 15C5 and 18C6 occurring upon encapsulation of the metal ions. The change in the IR intensity induced by complex formation is dependent on the concentration of metal ions in solutions, which provides the equilibrium constant for complex formation. The results of M+·15C5 in water and M+·18C6 in methanol are compared with those of M+·18C6 in water.10
We perform SEIRA spectroscopy of CEs tagged on gold surfaces by using a home-made vacuum deposition system and an FTIR spectrometer (Vertex 70, Bruker) coupled with a single reflection Si attenuated total reflection (ATR) system (VeeMAX II, PIKE Technologies). The details of the preparation of gold surfaces on a Si prism have been described in previous studies of Guo et al.13 Briefly, the gold surfaces are formed on an unheated Si crystal by thermal deposition of a gold wire from a tungsten basket in a vacuum chamber with a deposition rate of ∼0.005 nm s−1. The thickness of the gold surfaces is ∼7 nm, with which the enhancement of IR absorption will be optimum and the distortion of IR spectral features will be minimized.13 A home-made Teflon chamber is mounted on top of the Si crystal to introduce and retain liquid samples on the gold surfaces. Modification of the gold surfaces with CEs (1) or (2) is performed by monitoring IR absorption of CEs. A solution of (1) or (2) dissolved in dimethylsulfoxide (DMSO) with a concentration of ∼5 × 10−3 M is introduced on gold surfaces. The CEs (1) and (2) are tagged on the gold surfaces by forming S–Au chemical bonds. We leave the solution of CEs on gold surfaces for 30 minutes for achieving complete coverage of the gold surfaces with CEs. Fig. 1 shows the schematic structure of CEs on gold surfaces. Hereafter the gold surfaces modified with the CEs (1) and (2) are called 15C5-C1OC6 and 18C6-C1OC6, respectively. After removing the DMSO solution and washing the surfaces several times with pure water or methanol, FTIR difference spectra are measured repeatedly using solutions of alkali metal chloride salts and pure solvent with a 4 cm−1 spectral resolution at room temperature. In every exchange process of metal ions on the same gold surface, we confirm the removal of the metal ions from the gold surface by observing the complete disappearance of the IR bands due to the metal ion–crown ether complexes.
We calculate the structure and IR spectra of the metal ion complexes of 15C5-C1OC6 by terminating the end of the hydrocarbon chain with a CH3 group.10 The initial conformational search is performed for the Li+·15C5-C1OC6-CH3, Na+·15C5-C1OC6-CH3, K+·15C5-C1OC6-CH3, and bare 18C6-C1OC6-CH3 under the conditions of water solutions by using the CONFLEX program package and the MMFF94s force field with a search limit of 1 kcal mol−1, which produces more than 1000 conformers.14 The most stable ∼20 conformers found using the force field calculations are then optimized at the M05-2X/6-31+G(d) level with the polarizable continuum model (PCM) of water using the GAUSSIAN09 program package.15 The CONFLEX program does not include the force field of Rb+ and Cs+. For the Rb+, and Cs+ complexes, initial forms for the geometry optimization are produced by replacing the K+ ion in the stable conformers of the K+·15C5-C1OC6-CH3 complex with Rb+ or Cs+ ions. Vibrational analysis is carried out for the optimized structures at the same computational levels. Calculated frequencies at the M05-2X/6-31+G(d) level are scaled with a factor of 0.9389 for comparison with the observed IR spectra.10 The scaling factor was determined so as to reproduce the C–O stretching vibration of the K+·18C6-C1OC6 complex in water.10
Fig. 2 (a–e) The IR difference spectra of the M+·15C5-C1OC6 (M = Li, Na, K, Rb, and Cs) complexes on the gold surface using aqueous MCl solutions. The concentration of the MCl salts in water ranges from 10−6 to 1 M for Li+ and Na+ and from 10−3 to 1 M for K+, Rb+, and Cs+. The blue curves are the IR difference spectra recorded at 1 M. (f) The IR difference spectrum of the K+·18C6-C1OC6 complex on the gold surface with the aqueous solution of KCl at 100 mM (ref. 10). |
Fig. 3 The model simulation of band contours in IR difference spectra between CE complexes and bare CEs (see the text). |
We examine the dependence of the IR signal intensity on the concentration of MCl in aqueous solutions to obtain the equilibrium constant for the complex formation between M+ and 15C5-C1OC6. Fig. 4 shows the plot of the IR difference intensity at around 1100 cm−1 in Fig. 2 as a function of the concentration of MCl. The data are fitted by Hill equations (solid curves), and the apparent dissociation constant (KD) of the M+·15C5-C1OC6 complexes is obtained.11b The curves in Fig. 4 show a noticeable difference between the metal ions. The titration curve of the Na+ complex (Fig. 4b) shows the smallest KD value (∼2.8 × 10−5 M) among the alkali metal ion complexes in Fig. 4. The data points of the Na+ complex are highly scattered at the concentration more than 10−2 M; this is due to distorted band shapes in the spectra of the Na+ complex at higher concentration (see Fig. 2b). One possibility for the distorted band contour is that the Na+ ion may form 2:1 or higher complexes with 15C5 at higher concentration, changing the band shapes. The reciprocal of KD (KD−1) gives the equilibrium constant for the complex formation. The values of KD−1 for M+·15C5-C1OC6 in aqueous solutions are plotted in Fig. 5a (blue circles), which are compared with those for M+·18C6-C1OC6 (red circles) obtained in our previous study.10 The KD−1 value of Na+·15C5-C1OC6 is more than two orders of magnitude larger than those of the Li+, K+, Rb+, and Cs+ complexes, indicating the preference of Na+ for the encapsulation by 15C5, which is in good contrast to the K+ preference by 18C6.10
Fig. 4 The peak-to-peak amplitude at ∼1100 cm−1 in the IR difference spectra of the M+·15C5-C1OC6 (M = Li, Na, K, Rb, and Cs) complexes (Fig. 2) as a function of the concentration of the M+ ions in water (closed circles). The data are reproduced by Hill equations (solid curves). |
Fig. 5 The reciprocal of the apparent KD values determined by fitting the intensity at around 1100 cm−1 in the IR difference spectra of (a) M+·15C5-C1OC6 in water and (b) M+·18C6-C1OC6 in methanol (blue circles). For comparison, the data of M+·18C6-C1OC6 in water are shown in both panels (red circles) (ref. 10). |
The optimized (the most stable) structures of bare 15C5-C1OC6-CH3 and the M+·15C5-C1OC6-CH3 complexes are shown in Fig. 6. In the case of the Li+ and Na+ complexes, the metal ions are located in the 15C5 cavity. For the K+, Rb+, and Cs+ complexes, the ions are too large to enter the 15C5 cavity. The C–O–C parts of the crown rings tend to orient the oxygen atoms to the metal ions. The C–O–C frames of the Li+ and Na+ complexes form a planar-like configuration with the metal ions. For the K+, Rb+, and Cs+ complexes, since the metal ions are not in the center of the cavity but on it, the C–O–C frames show more bumpy features. Fig. 7 displays the calculated IR spectra of bare 15C5-C1OC6-CH3 and the M+·15C5-C1OC6-CH3 complexes in the 1000–1500 cm−1 region. A scaling factor of 0.9389 is employed for the calculated vibrational frequencies, and the solid curves are the IR spectra reproduced by providing a Lorentzian component with a FWHM of 10 cm−1 to each vibration.10 The IR spectrum of bare 15C5-C1OC6-CH3 (Fig. 7a) has a maximum intensity of the C–O stretching vibrations at 1122 cm−1, and it shifts to lower frequency with Li+ ions (1108 cm−1). The maximum position shifts further to lower frequency for Na+ ions (1103 cm−1), and then it is located closer to the position of bare species for K+ at 1112 cm−1, Rb+ at 1115 cm−1, and Cs+ at 1118 cm−1. The red shift of the C–O stretching vibration upon the complex formation was observed also for other CEs.9a,l,10 By using the calculated IR spectra, we obtain the IR difference spectra (Fig. 8). For the production of the IR difference spectra of the Li+ and Na+ complexes, the calculated IR intensity of the complexes is multiplied by 2 before the subtraction. For the K+, Rb+, and Cs+ spectra, on the other hand, the IR intensity of bare 15C5-C1OC6-CH3 is multiplied by 2. These simulated spectra in Fig. 8 well reproduce the spectral features of the observed ones (Fig. 2). In the simulated difference spectra of Li+ and Na+ (Fig. 8a and b), a positive band emerge at around 1106 cm−1, whereas the difference spectra of K+, Rb+, and Cs+ (Fig. 8c–e) show a negative signal at ∼1123 cm−1, rather than a positive one like the Li+ and Na+ complexes. It should be noted again that the main origin of the difference in the spectral features of the M+·15C5-C1OC6 complexes (positive signals for Li+ and Na+ and negative ones for K+, Rb+, and Cs+, Fig. 2a–e) is the difference of the relative IR intensity between the metal ion complexes and bare ones. As mentioned above, the calculated IR intensity for the M+·15C5-C1OC6 complexes in Fig. 7 cannot be used as it is for the reproduction of the IR difference spectra in Fig. 2 or 8. The discrepancy in the IR intensity calculated and observed can be explained by the orientation of the 15C5 frame against the gold surfaces after the inclusion of the metal ions. As seen in Fig. 6, the K+, Rb+, and Cs+ ions are displaced largely from the 15C5 cavity and likely to be exposed to the water phase more than the Li+ and Na+ ions. As a result, the M+·15C5 part in the K+, Rb+, and Cs+ complexes can take an orientation almost parallel to the water phase, or the gold surfaces, as drawn in Fig. S2a of the ESI.† For the M+·15C5-C1OC6 complexes, the dipole derivatives of the C–O stretching vibrations are almost parallel to the 15C5 cavity (Fig. S2b and c of the ESI†). Osawa and coworkers demonstrated that in SEIRA spectroscopy a vibrational mode with its dipole derivative perpendicular to the gold surface is preferentially enhanced.11c Hence, the C–O stretching vibration of the M+·15C5-C1OC6 (M = K, Rb, and Cs) complexes is less enhanced due to the parallel orientation, resulting in weaker IR intensity. The origin of the more enhanced nature of the Li+ and Na+ complexes is not clear at the present stage, but probably the M+·15C5 part can take a perpendicular orientation to the gold surface with Li+ and Na+, due to the repulsive force between the neighboring complexes.
Fig. 6 The optimized structure of bare 15C5-C1OC6-CH3 and the M+·15C5-C1OC6-CH3 complexes. The numbers in the figure show the M+⋯O distance (Å). |
Fig. 8 The IR difference spectra made by subtracting the simulated IR spectrum of bare 15C5-C1OC6-CH3 in Fig. 7a from those of the M+·15C5-C1OC6-CH3 complexes in Fig. 7b–f. In the simulation of the Li+ and Na+ complexes, the calculated IR intensity of the M+·15C5-C1OC6-CH3 complexes is multiplied by 2 before the subtraction, while for the K+, Rb+, and Cs+ spectra, the IR intensity of bare 15C5-C1OC6-CH3 is multiplied by 2 (see the text). |
We then return to the results in Fig. 2 and 5a. The IR spectral features reflect the strength of the interaction between the metal ions and the crown ether in the complexes. From the similarity of the IR spectra in Fig. 2, the Li+ and Na+ ions undergo comparable interaction with 15C5-C1OC6, and the interaction in the K+, Rb+, and Cs+ complexes is similar to each other and should be weaker than that in the Li+ and Na+ complexes. On the other hand, the equilibrium constant determined from the IR signal intensity shows clear preference for Na+ over the other metal ions. From these findings, we can attribute the ion selectivity of Na+ by 15C5-C1OC6 to the energy before and after complex formation. Fig. 9 shows the schematic energy diagram proposed for the M+·15C5-C1OC6 complexes in water. The left side of Fig. 9 shows the energy levels of free ions and 15C5-C1OC6 in water, which are drawn on the basis of reported hydration molar free energy of alkali metal ions: −124.0 (Li+), −100.1 (Na+), −82.5 (K+), −77.4 (Rb+), −69.7 (Cs+) kcal mol−1.16 Our IR spectra probe the complex structure or the strength of the interaction in M+·15C5-C1OC6 on the right side of Fig. 9. The spectral resemblance in Fig. 2 indicates that the interaction energy between M+ and 15C5-C1OC6 is similar for Li+ and Na+, and for K+, Rb+, and Cs+, which is illustrated on the right side of Fig. 9. The equilibrium constant for the complex formation is governed by the difference in free energy (ΔG). The observed preference of Na+ over K+, Rb+, and Cs+ can be attributed to the strength of the interaction or the size matching in size between the metal ions and 15C5-C1OC6 (the right side of Fig. 9). In contrast, higher selectivity for Na+ over Li+ ions is attributed to the difference of solvation energy in aqueous solution between Na+ and Li+ (the left side of Fig. 9).9h,m
Footnote |
† Electronic supplementary information (ESI) available: The IR absorption spectra of the thiol derivatives used in this study (Fig. S1). The proposed orientation of the Rb+·15C5 component at the interface, and the dipole derivative of the most intense C–O stretching band for Na+·15C5-C1OC6-CH3 and Rb+·15C5-C1OC6-CH3 (Fig. S2). The concentration dependence of the peak-to-peak intensity at around 1100 cm−1 in the IR difference spectra of the M+·18C6-C1OC6 complexes using the methanol solutions of MCl (Fig. S3). The optimized structures of bare 15C5-C1OC6-CH3 and the M+·15C5-C1OC6-CH3 (M = Li, Na, K, Rb, and Cs) complexes in water (Fig. S4–S9). See DOI: 10.1039/c5nj01787d |
This journal is © The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2015 |