Andrei
Filippov
*ab,
Oleg I.
Gnezdilov
b,
Nicklas
Hjalmarsson
c,
Oleg N.
Antzutkin
a,
Sergei
Glavatskih
de,
István
Furó
f and
Mark W.
Rutland
cg
aChemistry of Interfaces, Luleå University of Technology, SE-97187 Luleå, Sweden. E-mail: Andrei.Filippov@ltu.se; Tel: +46 73-6782225
bInstitute of Physics, Kazan Federal University, 420008 Kazan, Russia
cSurface and Corrosion Science, KTH Royal Institute of Technology, SE-10044 Stockholm, Sweden
dSystem and Component Design, KTH Royal Institute of Technology, SE-10044 Stockholm, Sweden
eDepartment of Mechanical Construction and Production, Ghent University, B-9000 Ghent, Belgium
fApplied Physical Chemistry, KTH Royal Institute of Technology, SE-10044 Stockholm, Sweden
gRISE Technical Research Institute of Sweden, Chemistry Materials and Surfaces, Box 5607, SE-11486 Stockholm, Sweden
First published on 7th September 2017
Diffusion of EAN confined between polar glass plates separated by a few micrometers is higher by a factor of ca. 2 as compared to bulk values. Formation of a new phase, different to the bulk, was suggested.
Recently, small-angle neutron scattering revealed that EAN itself exhibits an inherent amphiphilic nanostructure in the pure liquid state.9 For the lamellar structures deemed as most probable, the calculated Bragg spacing is approximately twice the ion-pair dimension of the IL, which suggests that the IL is structured on the length scale of the ions, with the (hydrogenous) alkyl groups associated together and segregated from the H-bonded ionic moieties –NH3+ (or –ND3+) and NO3−. Based on this finding, it has been suggested that surfaces may induce alignment of such locally ordered domains, creating a persistent molecular layering.10 Vibrational sum frequency spectroscopy and X-ray reflectivity studies11 confirmed the existence in EAN of significant interfacial structures within a ca. 3 nm thick layer. Surface structures have also been found in other ILs,12–14 sometimes persisting up to tens of nanometers into the bulk13,14 with a clear dependence on the molecular nature of the surface. Exceptionally extensive surface-induced structures with a thickness of up to 2 μm were recently demonstrated for a number of imidazolium bis(trifluoromethylsulfonyl)imide ILs by Anaredy and Shaw.15 Much less attention has been paid to the structure and dynamics of ILs, which can be formed in a micrometer-scale confinement.
The purpose of this work was to study the dynamic features exhibited by EAN in the presence of polar surfaces and micrometer-scale confinement. For this purpose, we exploited NMR diffusion measurements, previously applied to a wide range of ILs.16–22 Of particular relevance, diffusion NMR has been used to study the molecular dynamics in ILs under confinement (coincidentally, all in silica pores).18,23–27 In comparison to NMR relaxation studies, diffusion NMR data are typically simpler to interpret and are more robust regarding experimental artefacts.28,29 Both simplicity and robustness are of great advantage for ILs with their frequently complex behaviour. In contrast to other transport methods, NMR diffusion data can be obtained selectively for different constituting moieties, in the present case the ethylammonium cation and the nitrate anion. This latter feature, as also shown here, can be aided by suitable isotope enrichment (deuterium, 2H, and 15N).
EAN (see Fig. 1) consists of an ethylammonium (EA) cation and a nitrate (NI) anion. All samples were synthesized as described in the ESI,† following previously established protocols.30,31 While Sample 1 of EAN has a natural abundance (n.a.) of all isotopes: 1H (99.98%), 13C (1.108%) and 15N (0.37%) in both EA and NI, the nitrate anion was 15N enriched (∼98 atom%) in Sample 2 of EAN. Sample 3 of EAN was prepared by admixing EAN-D5, i.e. selectively deuterated at CD3- and CD2-positions (∼98 atom% of 2H) of the cation, with Sample 1 at a composition of 20/80 mol% of EAN-D5/EAN (n.a.). Prior to experiments, each sample was degassed under vacuum (at a pressure of less than 10−3 mbar and temperature of 313 K) for 60 hours to remove absorbed water and gases. The chemical composition of each sample (including impurities) was established by liquid 1H and 13C NMR and is presented in the ESI.† In summary, the hydrogen and carbon contents associated with the ethylammonium cation were approximately 99% in Sample 1 and more than 99.6% in Sample 2.
![]() | ||
Fig. 1 The chemical structure of EAN consisting of an ethylammonium cation and a nitrate anion, with nitrogen (blue), oxygen (red), carbon (dark grey) and hydrogen (light grey) atoms shown. |
In a series of experiments, samples of EAN were confined between thin glass plates arranged in stacks (Fig. S1, ESI†). The plates used (5 × 14 × 0.1 mm, Thermo Scientific Menzel Gläser, Menzel GmbH, Germany) were cleaned carefully before being used (see the ESI†). Contact angle measurements with Milli-Q water provided a contact angle near 0° (see the ESI†), indicating hydrophilic glass surfaces. A stack of glass plates filled with EAN was prepared in a glove box in a dry nitrogen atmosphere. A drop of EAN was placed on the first plate, then it was covered by the second plate by placing another drop of EAN on the top, and so on until the thickness of the stack reached approximately 5 mm. Excess of EAN from the sides of the stack was removed by wiping. The plates were thereafter placed in a sample cell of square cross section. The mean distance (spacing) between the glass plates was estimated by weighing the introduced EAN (see the ESI†), which yielded d ∼ 3.8–4.5 μm. Measuring the thickness directly and subtracting it from the total glass thickness (see the ESI†) indicates that this is a consistent value.
The 1H (Samples 1 and 2) and 15N (Sample 2) NMR self-diffusion measurements of EAN in the bulk were performed using a Pulsed-Field-Gradient (PFG) probe Diff50 (Bruker), and a Bruker Ascend Aeon WB 400 (Bruker BioSpin AG, Fällanden, Switzerland) NMR spectrometer was used with working frequencies of 400.27 MHz for 1H and 40.56 MHz for 15N. Data were processed using the Bruker TopSpin 3.5 software. Bulk samples of EAN (approximately 300 μl each) were placed in a standard 5 mm glass NMR sample tube and sealed with a plastic cap to avoid contact with air. Before each measurement, the sample was equilibrated at the specified temperature for 20 min. The 2H NMR spectra of Sample 3 were obtained using the same spectrometer with the quadrupole-echo pulse sequence.
The 1H NMR diffusion measurements of EAN confined in the glass stack were performed using a NMR goniometer probe, which was used to macroscopically align the films in the sample stack at angles 0° and 90° with respect to the applied magnetic field gradient, the latter being co-parallel to the main magnetic field.
The diffusional decays (DDs) were recorded using the stimulated echo pulse sequence. For single-component isotropic diffusion, the signal intensity A changes as:32,33
![]() | (1) |
The anion lacks hydrogen and carbon. Hence the 1H and 13C NMR spectra of EAN arise exclusively from the EA cation (Fig. S2–S5 in the ESI†). These spectra demonstrated three resonance lines, which were assigned to the NH3+, CH2 and CH3 moieties, in accordance with previously published data.34 As expected, the 15N NMR spectrum of the nitrate anion in Sample 2 revealed a single broad resonance line (Fig. S6 in the ESI†) at ca. 380 ppm, in accordance with literature data.35
The DDs of the 1H signals of the EA cation in Sample 1 (Fig. S7, ESI†) and Sample 2, as well as the DD of the 15N signal of the nitrate anion (Fig. S8, ESI†), are all linear in the semi-logarithmic scale. Therefore, the mobility of each ion can be described by a single diffusion coefficient, the values of which (D) are tabulated in Table S1 (ESI†). The D values of the EA cation in both these samples coincide; therefore, any slight difference in the composition as that for Samples 1 and 2 does not influence the cation translation dynamics. An important observation is that the D values of the nitrate anion are higher by a factor of ∼1.7 than those of the EA cation. This means that the diffusive motion of cations and anions in bulk EAN is not tightly coupled (that is, ion pairing is not dominant). In dense ionic systems, self-diffusion is influenced strongly by electrostatic interactions and is no longer defined solely, as in simple aprotic liquids built up by neutral molecules, by molecular size.21 In EAN, hydrogen bonding, too, is a significant intermolecular force with concomitant solvophobicity.9 To exemplify this point, one should recall that the viscosity of EAN is ca. 30 times higher than that of nitropropane or similar liquids.36,37 In addition, EAN exhibits a nano-scale order9 instead of a random local structure. Hence, it is difficult to pinpoint a single dominant reason for the factor-1.7 difference. We note that in another IL with nano-scale ordering, one could find in magnitude a rather similar (factor-2) ratio between anion and cation diffusion coefficients.22 Yet, in that system, this ratio seemed to depend strongly on the nano-scale order that clearly shows the importance of local structure for the translational dynamics of the ions.
Fig. 2 shows the DDs of the 1H signals of the EA cation with the field gradient directed along and normal to glass plates. Under current conditions, the DDs recorded are caused by diffusive molecular displacements in the IL films along and normal to the glass plates, respectively.
![]() | ||
Fig. 2 The DDs of the 1H NMR signals of the EA cation recorded at 293 K by pulsed-field-gradient stimulated echo experiments in bulk EAN (solid squares, D0 = 3.0 × 10−11 m2 s−1) and in EAN films confined between parallel glass plates, with gradient directions along (A) and perpendicular to (B) the plates. The DDs were obtained by increasing the gradient amplitude g, and the different colored symbols indicate data recorded at different diffusion times td. Best fits (see the ESI†) of suitable expressions for restricted diffusion to the data in (B) are presented by the lines of corresponding colors. |
In the direction along the plates (Fig. 2A), the DDs maintain a linear form typical of the bulk state, but the slope (dotted line) of the decays and, consequently, the value of the derived self-diffusion coefficient D‖ = D* is a factor of ca. 1.9 higher than D0 in the bulk (the decay for which is also shown for comparison). The DD is invariant to the changing diffusion time in the range of 50–1000 ms.
In the direction normal to the plates (Fig. 2B, colored symbols) and short diffusion times (3 and 10 ms), the signal decays faster than that in bulk (black squares). The initial slopes of the DDs for these short diffusion times correspond to D*. This means that the diffusion is accelerated to a similar extent along as well as normal to the plates, D⊥ = D‖ = D*.
The forms of DDs (Fig. 2B) deviate from the single-exponential behavior typical for the bulk (recall eqn (1), shown here as a solid line and black squares), and the extent of this deviation increases as the diffusion time increases from 3 ms to 1 s. The slopes of these DDs decrease with an increase in diffusion time, revealing that molecular displacements across the IL films confined between glass plates are hindered.28,32 The hindrance is constituted by the glass plates, whose separation of ca. 4 μm is comparable to the diffusional path-length of ions in EAN, L ≈ √(2D*·td),38 which is, in the current case, between 1.8 and 5.8 μm for td ranging from 30 to 300 ms, respectively. As is well known, in such situations the mean displacement levels off, while the apparent mean diffusion coefficient decreases with increasing td revealing the same qualitative behavior as shown in Fig. 2B. Planar restriction (diffusion along the plane normal between parallel plates) is one of the simplest regular geometries, for which the DD of the confined liquids has been analytically solved,39 under the assumption of elastic collisions of diffusing molecules with the wall. Fits of that theoretical expression under assumptions concerning the distribution of plate–plate separation in the prepared sample stack (see Fig. S9 in the ESI†) provided a mean distance between planes of 4.1 μm that agrees rather well with the experimentally estimated value of ∼3.8 μm. In summary, the data in Fig. 2B are well explained by simple mechanical hindrance exerted by the glass plates on the diffusive displacement of cations (with D*) in the IL film.
The DDs in Fig. 2A that arise from diffusive displacements in the IL films along the plates exhibit no dependence on the diffusion time and are thereby not influenced by any restrictions. This feature is as expected for free diffusion of ions along the IL films. Yet, in remarkable contrast, diffusion along the plates seems to be characterized by a diffusion coefficient D‖ = D* that is higher by a factor of ∼1.9 than the value of D0 observed in the bulk. This increase is so large that it is comparable to that obtained in the bulk by increasing the temperature by ca. 15 K.
15N NMR diffusion measurements of EAN confined between plates were performed in order to investigate the effect of confinement on the diffusion of the nitrate anion. However, weak signals were obtained due to the lower filling factor of EAN between the plates as compared to the bulk sample (the ratio of the filling factors was estimated to be 0.001). The gyromagnetic ratio of 15N is lower than that of 1H by a factor of 10, while the number of 15N spins in EAN is lower than the number of protons by a factor of 8. This is the main reason that we could measure diffusion of the cation (1H) but not of the anion (15N).
In addition to PFG NMR experiments, the T1 and T2 NMR relaxation times of protons were also measured for the bulk and confined EAN (see the ESI† for details). Remarkably, the transverse relaxation time T2 of –NH3 protons in EA cations in the confinement between glass plates is shorter than that of EA cations in the bulk by a factor of ca. 22 (see Fig. S10 of the ESI†), which can be explained, for example, by a significant acceleration in the exchange of NH3 protons in the whole layer of EAN confined between polar glass plates. Note that the presence of Si–OH and H2O protons at glass surfaces can only partly be invoked in the model explaining such a dramatic change in T2, since the vast amount of EAN ions is located tens to hundreds of ionic radii away from the surfaces.
After having made these remarkable observations, it became imperative to consider and rule out possible experimental artifacts. Hence, the experimental set-ups and samples were carefully tested. Particular attention was paid to the possible effect of water; since the stack contained many glass plates, one could assume that water adsorbed at the glass surface may have dissolved in the IL films and contributed to faster molecular dynamics. To exclude this as a significant contribution, the diffusion coefficient in bulk EAN with increasing water content has been measured (see Fig. S11–S13 in the ESI†). From these data and from the 1H NMR spectra, from which the approximate water content in EAN could be estimated, we can conclude that the effect of water in the observed phenomenon is negligible (approximately 8 wt% water would be required to modify the diffusivity to the required extent, whereas the water content was shown to be <0.1 wt%). Another possible source of artefact could be a capillary force, which might force the liquid to flow in narrow spaces. However, the samples were prepared some time before performing the experiments; therefore it seems unlikely that the experiments were not performed under equilibrium conditions. The only source of molecular motion was the thermal energy. Gravity might also lead to the flow in the thin layer of a liquid. In this case, the effect would depend on the orientation of the stack of glass plates relative to the gravity field of the Earth. No effects of orientation on the measured diffusion coefficient were observed. In summary, we conclude that the observed acceleration of diffusion between the confining glass plates is not an artifact, but an utterly surprising molecular (or structural) effect. This may affect not only the diffusion of EAN but also other physical parameters such as density. In the following we provide, with the support of some additional data, a tentative explanation to this observation.
Previous studies have established the notion that some surfaces can exert an effect on ILs that permeate the bulk for tens of nanometers. This effect was typically described as an imposed molecular order.9,11–15 Other studies have shown that changes in short-range molecular order can have an influence on the macroscopic (displacements over micrometers) self-diffusion of the cations and anions in ILs.22 On this basis, we provide a tentative interpretation of our findings.
The diffusion path length over 3 ms with D0 = 3 × 10−11 m2 s−1, the bulk diffusion coefficient, is (2·D0·td)0.5 ∼ 0.42 μm, or if we take D* = 5.6 × 10−11 m2 s−1, the diffusion path is ∼0.58 μm. However, the DDs obtained for diffusion along the plates at 3 ms are single-exponential and coincide with the DDs obtained at longer diffusion times (Fig. 2A). This indicates that EAN diffusion is the same as in the center of a layer and near the plates (at a proximity of ∼0.5 μm to the plates). Therefore, diffusion is accelerated not only in the thin layers near the surfaces, but in the whole thickness of the layer, even at the distance as far as ∼2 μm from the surface.
Regarding the molecular order within the surface layer, if it is orientational akin to that exhibited by liquid crystals it may leave an imprint in NMR line splitting if quadrupolar nuclei, such as 2H, are used as an internal (in EA-D5) or external (D2O as a solute in EAN) spin-probe. This is the case if the structure within the layer is similar to that in liquid crystal phases with non-zero second-rank order parameter S, such as lamellar/smectic phases. In contrast, phases with S = 0, such as cubic LCs, do not exhibit line splitting. To investigate this, we recorded the 2H NMR spectra of Sample 3 prepared with EAN-D5. As shown by the NMR spectra (see Fig. S14 in the ESI†), there is no detectable quadrupolar splitting in the 2H NMR spectra. This finding suggests that either the surface layer is thin, in which case the 2H NMR signal from it remains hidden in the background (at λ = 60 nm, only ∼3% of the NMR signal arises from the surface layer), or that the molecular order is such that the order parameter S remains zero.
The diffusivity of EA cations is clearly shown to be isotropic when confined between glass plates. The bulk diffusion properties are highly consistent with the literature consensus of a sponge-like self-assembly structure. Consideration of spin echo data, where no significant splitting is observed, further allows the conclusion to be drawn that the majority of cations are unaffected by the glass surfaces. Water is often present in ionic liquids and is fully compatible with EAN. In this work, it is categorically shown by exhaustive measurements that water associated with the polar glass plates cannot be solely responsible for the dramatic change in diffusivity compared to the bulk. At the same time, the confinement leads to acceleration in the exchange of NH3 protons in the whole layer. The most plausible idea for the accelerated diffusion of ions is that the nano-/microscale structural features of EAN confined between polar glass plates are considerably different from the bulk phase structure of EAN. We suggest that the mechanism of the increased rate of relaxation (decreased T2) of –NH3 protons is an exchange process. Indeed, the maximum effect occurs only with “exchangeable” –NH3 protons; however, it is not the same with “non-exchangeable” –CH2– and –CH3 protons. –NH3 protons are really generalized between –NH3 of the EA cation and the NO3− anion. In bulk EAN, they demonstrate a narrow line (as shown in Fig. S2 in the ESI†). Therefore, they are in the condition of “fast exchange”, where the difference in the chemical shift frequencies of the two sites is much less than the exchange rate.40 This is a case of the so-called “motional averaging” exchange regime. The acceleration of the rate of T2 relaxation of –NH3 protons for EAN placed between polar glass plates (see Table S2 in the ESI†) and broadening of the corresponding line in the 1H NMR spectrum (Fig. S12 in the ESI†) indicate that the exchange rate decreases and now –NH3 protons are in the “intermediate exchange” regime where the difference in the chemical shift frequencies of the two sites is of the order of the exchange rate. This transition from the “motional averaging” to the “intermediate” regime can be related to the suggested structural change of EAN occurring between polar glass plates. Indeed, in the sponge-like phase of bulk EAN, the ions –NH3+ and NO3−, as suggested, are located in a separate micro-phase, which is different from another micro-phase containing –CH2– and –CH3 groups. This hinders diffusion but facilitates the proton exchange process. If the structure changes to a more homogeneous one, it might facilitate the diffusion of ions but will hinder proton exchange.
More work and other experimental techniques, such as the effect of surface roughness, surface polarity and variable hydrophobicity of the ionic liquid constituents, are needed to further clarify this remarkable phenomenon. Further experiments with different ionic liquids (firstly, ammonium nitrates), different substrates, different inter-plate gapes, without and with added solvents and salts would help to address the physics of the system.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c7cp01772c |
This journal is © the Owner Societies 2017 |