Romain Lhermerout

Susan Perkin
email: susan.perkin@chem.ox.ac.uk

Nanoconfined ionic liquids: disentangling electrostatic & viscous forces

Recent reports of surface forces across nanoconfined ionic liquids have revealed the existence of an anomalously long-ranged interaction apparently of electrostatic origin. Ionic liquids are viscous and therefore it is important to inspect rigorously whether the observed repulsive forces are indeed equilibrium forces or, rather, arise from the viscous force during drainage of the fluid between two confining surfaces. In this paper we present our direct measurements of surface forces between mica sheets approaching in the ionic liquid [C2C1Im][NTf2], exploring three orders of magnitude in approach velocity. Trajectories are systematically fitted by solving the equation of motion, allowing us to disentangle the viscous and equilibrium contributions. First, we find that the drainage obeys classical hydrodynamics with a negative slip boundary condition in the range of the structural force, implying that a nm-thick portion of the liquid in the vicinity of the solid surface is composed of ordered molecules that do not contribute to the flow. Secondly, we show that a long-range static force must indeed be invoked, in addition to the viscous force, in order to describe the data quantitatively. This equilibrium interaction decays exponentially and with decay length is in agreement with the screening length reported for the same system in previous studies. In those studies the decay was simply checked to be independent of velocity and measured at low approach rate, rather than explicitly taking account of viscous effects: we explain why this gives indistinguishable outcomes for the screening length by noting that the viscous force is linear to very good approximation over a wide range of distances.

I. INTRODUCTION

The Surface Force Apparatus (SFA) or Balance (SFB) was originally designed by Tabor and Winterton to measure Van der Waals interactions between macroscopic bodies surrounded by air or vacuum [START_REF] Israelachvili | Intermolecular and Surface Forces[END_REF][START_REF] Tabor | The direct measurement of normal and retarded van der waals forces[END_REF]. Subsequent extensions by Israelachvili led to measurement of surface interactions across liquids [START_REF] Israelachvili | Measurement of forces between two mica surfaces in aqueous electrolyte solutions in the range 0-100 nm[END_REF] and many varieties of soft matter. Aqueous solutions of ions have been used extensively to study electrostatic interactions [START_REF] Israelachvili | Measurement of forces between two mica surfaces in aqueous electrolyte solutions in the range 0-100 nm[END_REF][START_REF] Horn | Double-layer and solvation forces measured in a molten salt and its mixtures with water[END_REF][START_REF] Pashley | DLVO and hydration forces between mica surfaces in Li+, Na+, K+, and Cs+ electrolyte solutions: A correlation of double-layer and hydration forces with surface cation exchange properties[END_REF][START_REF] Pashley | Molecular Layering of Water in Thin Films between Mica Surfaces and Its Relation to Hydration Forces[END_REF]; in particular this allowed validation of the DLVO theory of colloidal stability in dilute electrolytes. In these works, the viscosity and velocity were small enough to safely neglect viscous effects and to consider only equilibrium forces.

Separately, Chan & Horn took advantage of the capabilities of the SFA to investigate the viscous force in nanoconfined liquids [START_REF] Chan | The drainage of thin liquid films between solid surfaces[END_REF]. They focused on purely dissipative forces by using apolar liquids, and showed that drainage is quantitatively well-described by classical hydrodynamics down to the nanoscale. For the liquids used, that completely wet mica surfaces, the best agreement was obtained by considering a non-slip boundary condition a few nanometers away from the solid surface in the liquid phase, thereafter reported as negative slip. From this observation emerged the picture of a stopped layer, namely a thin film in the vicinity of the solid surface that is part of the liquid phase but does not contribute to the fluid flow. Since the pioneering work of Chan & Horn, many studies looked at the hydrodynamic boundary condition with the SFA, generalizing the analysis to many systems [START_REF] Israelachvili | Measurement of the Viscosity of Liquids in Very Thin Films[END_REF][START_REF] Horn | Molecular organization and viscosity of a thin film of molten polymer between two surfaces as probed by force measurements[END_REF][START_REF] Israelachvili | Shear Properties and Structure of Simple Liquids in Molecularly Thin Films: The Transition from Bulk (Continuum) to Molecular Behavior with Decreasing Film Thickness[END_REF][START_REF] Horn | Surface forces and viscosity of water measured between silica sheets[END_REF], some of them presenting partial slip [START_REF] Vinogradova | Drainage of a Thin Liquid Film Confined between Hydrophobic Surfaces[END_REF][START_REF] Bocquet | Nanofluidics, from bulk to interfaces[END_REF] or elastohydrodynamic deformation [START_REF] Leroy | Hydrodynamic interactions for the measurement of thin film elastic properties[END_REF][START_REF] Leroy | Hydrodynamic Interaction between a Spherical Particle and an Elastic Surface: A Gentle Probe for Soft Thin Films[END_REF][START_REF] Villey | Effect of Surface Elasticity on the Rheology of Nanometric Liquids[END_REF][START_REF] Wang | Out-of-Contact Elastohydrodynamic Deformation due to Lubrication Forces[END_REF][START_REF] Wang | Elastic deformation during dynamic force measurements in viscous fluids[END_REF][START_REF] Wang | Elastic deformation of soft coatings due to lubrication forces[END_REF].

Very recently ionic liquids (ILs) have been the focus of intense study because of their exceptional properties (negligible volatility, nonflammability, high ionic conductivity, wide electrochemical and thermal windows) which make them extremely promising for a range of engineering and electrochemical applications (energy storage, lubrication, catalysis, electrowetting, etc.) [START_REF] Fedorov | Ionic Liquids at Electrified Interfaces[END_REF][START_REF] Hayes | Structure and Nanostructure in Ionic Liquids[END_REF]. Forces of apparently electrostatic origin have been reported with very substantially longer range than anticipated based on existing theories, leading to an active debate in the community as to their interpretation [START_REF] Gebbie | Ionic liquids behave as dilute electrolyte solutions[END_REF][START_REF] Perkin | Is a stern and diffuse layer model appropriate to ionic liquids at surfaces?[END_REF][START_REF] Gebbie | Reply to Perkin et al.: Experimental observations demonstrate that ionic liquids form both bound (Stern) and diffuse electric double layers[END_REF][START_REF] Espinosa-Marzal | Ionic Liquids Confined in Hydrophilic Nanocontacts: Structure and Lubricity in the Presence of Water[END_REF][START_REF] Cheng | Characterizing the Influence of Water on Charging and Layering at Electrified Ionic-Liquid/Solid Interfaces[END_REF][START_REF] Lee | Are Room-Temperature Ionic Liquids Dilute Electrolytes?[END_REF][START_REF] Gebbie | Long-range electrostatic screening in ionic liquids[END_REF][START_REF] Smith | The Electrostatic Screening Length in Concentrated Electrolytes Increases with Concentration[END_REF][START_REF] Hjalmarsson | Switchable long-range double layer force observed in a protic ionic liquid[END_REF][START_REF] Comtet | Nanoscale capillary freezing of ionic liquids confined between metallic interfaces and the role of electronic screening[END_REF][START_REF] Gebbie | Long range electrostatic forces in ionic liquids[END_REF][START_REF] Lee | Underscreening in concentrated electrolytes[END_REF][START_REF] Lee | Scaling Analysis of the Screening Length in Concentrated Electrolytes[END_REF]. To measure such long-range forces, these studies used the same protocol as for aqueous solutions: using low enough velocities for the dynamical effects to be negligible, i.e. checking that the measured force was independent of velocity over a range of slow velocities. However ionic liquids have dynamic viscosity typically 100 times greater than that of water, consequently one can estimate the viscous force to be of the same order of magnitude as the reported electrostatic forces at the approach velocities employed. It is therefore necessary to formally disentangle the equilibrium and viscous contribution to the measured forces [START_REF] Horn | A reevaluation of forces measured across thin polymer films: Nonequilibrium and pinning effects[END_REF]. More generally, it is of interest to ask whether classical hydrodynamics applies to nanoconfined ionic liquids, and whether a slip or stick boundary condition applies. Indeed, ionic liquids are composed solely of charged molecules with strong electrostatic interactions that could potentially modify the drainage rate because of coupled electro-diffusio-osmotic mechanisms. In particular, when the liquid is confined at a scale within the range of a long-range electrostatic interaction, electroneutrality is broken and one can expect substantial electrokinetic effects [START_REF] Bocquet | Nanofluidics, from bulk to interfaces[END_REF][START_REF] Levine | Theory of electrokinetic flow in a narrow parallel-plate channel[END_REF].

Thus a systematic study is needed in order to tackle these dynamical aspects. This paper is dedicated to answering the following questions. Does classical hydrodynamics apply for nanoconfined ILs and what is the slip boundary condition? When considering quantitatively the viscous force during drainage, does ILs still exhibit a long-range electrostatic interaction? We report SFB measurements exploring 3 orders of magnitude in velocity, and the equations of motion combining viscous and electrostatic forces have been systematically solved to analyze the data. Such approach allowed us to disentangle static and dynamic contributions to the total force when the ionic liquid, here

[C 2 C 1 Im][NTf 2 ]
, is confined between approaching mica surfaces. We found first no significant deviation from Reynolds theory of fluid drainage, and that classical hydrodynamics holds down to the nanoscale with a negative slip length of 4.2 ± 0.7 nm, i.e. a layer of a few molecules that are stopped in the vicinity of the solid surface. We then show that an additional equilibrium, long-range exponentially decaying contribution has to be invoked to describe the total force quantitatively. The decay length of this interaction is found to be 8.0 ± 0.5 nm, in agreement with the previously reported values for this system. We finally come back to the interpretation of previous studies, and explain why ignoring a non-negligible viscous force turns out nonetheless to provide a correct evaluation of the range of the electrostatic force.

The manuscript is organized as follows. The experimental set-up and the protocol of analysis is presented in section II; the measurement at the highest velocity is shown in section III A in order to characterize the purely viscous response; data at very low velocities are then used in section III B to reveal the equilibrium electrostatic repulsion; and an intermediate velocity is finally explored in section III C to comment on previous analyses. The principle of the SFB is illustrated in Figure 1(a), and explained in detail elsewhere [START_REF] Israelachvili | Intermolecular and Surface Forces[END_REF][START_REF] Klein | Simple liquids confined to molecularly thin layers. I. Confinement-induced liquid-to-solid phase transitions[END_REF][START_REF] Perkin | Forces between Mica Surfaces, Prepared in Different Ways, Across Aqueous and Nonaqueous Liquids Confined to Molecularly Thin Films[END_REF]. Briefly, muscovite mica is cleaved to produce micrometer-thick facets that are molecularly smooth, which are then backsilvered and glued onto cylindrical glass lenses with an epoxy resin (EPON 1004, Shell Chemicals). Two lenses made with mica of the same thickness are arranged in a cross-cylinder geometry, and separated by a liquid bridge. The liquid used is 1-ethyl-

II. MATERIALS AND METHODS

3-methylimidazolium bis[(trifluoromethane)sulfonyl]imid [C 2 C 1 Im][NTf 2 ]
(Iolitec, 99%), which chemical structure is given in Figure 1(c), notably because it is one of the less viscous ILs (dynamic viscosity η = 37 mPa.s at 25 • C [START_REF] Yu | Viscosity of ionic liquids: Database, observation, and quantitative structure-property relationship analysis[END_REF]).

The liquid is dried in a Schlenck line at 60 • C and 2 × 10 -2 mbar for 10 hours, and injected in the chamber just after. Typically one hour is waited before measurements for the apparatus to be in thermal equilibrium with the room regulated to 25 • C.

The top lens is translated with a piezo-electric tube (Ferroperm PZ 29) that allows a very smooth and linear motion over a typical range of 200 nm, at velocity v controlled between 0.07 nm/s (limited by drifts) and 100 nm/s (limited by the camera frame rate). White light interferences between the silver mirrors are observed though a spectrometer, and the so-called FECO (Fringes of Equal Chromatic Order) pattern is recorded. Knowing the refractive index of

[C 2 C 1 Im][NTf 2 ] (1.42251 at 25 • C [40]
), it is possible to extract the radius of curvature R = 8.85 mm and the liquid thickness D [START_REF] Israelachvili | Thin Film Studies Using Multiple-Beam Interferometry[END_REF]. The bottom lens is mounted on a spring (spring constant k = 125 N/m), deflection of which allows us to deduce the normal force between the two surfaces.

Since the last publication of our group reporting force profiles [START_REF] Smith | Switching the Structural Force in Ionic Liquid-Solvent Mixtures by Varying Composition[END_REF], our setup has been improved to significantly reduce the noise level on the signal D(t). First, a QImaging Retiga R6 camera provides a high resolution (2688 × 2200 px 2 ), a better signal to noise ratio and a FPS of 14 Hz. Secondly the apex position of the FECO fringes is extracted from a local parabolic fitting. Thirdly, a box made of expanded polystyrene sheets surrounds the SFB in order to damp the vibrations coming from air. These modifications allows us to achieve a RMS noise of 0.02 nm and so very clean signals (see Figure 1(b)) that are ideal for fitting weak long-range forces.

In this version of the SFB, only the distance D between the lenses is measured directly, and the total force is deduced from it with the relation F (t) = k (D(t) -D 0 + vt), where D 0 is the starting distance, which is supposed to be large enough to have F (0) = 0. Here the velocity v (resp. force F ) has been taken positive for approaching motion (resp for repulsive interaction). A difficulty is that the velocity v of the top lens is not precisely known in advance, as there is a 10% fluctuation from run to run, and a calibration of such uncertainty prevents from reliably identifying weak long-range forces. The method commonly used to solve this problem consists in fitting linearly the D(t) trajectory at distances large enough to have negligible long-range forces, so that the fitted slope corresponds to v. This method is excellent when dealing with aqueous solutions and velocities such that the viscous force is below the limit of detection of the apparatus, but does not necessarily apply for viscous liquids and velocities such that viscous force is measurable. In this paper, we follow the approach of Chan & Horn by solving the equation of motion of the bottom lens and fitting the D(t) trajectories [START_REF] Chan | The drainage of thin liquid films between solid surfaces[END_REF].

Neglecting the inertia of the bottom lens, the equation of motion in the frame of the lab is simply the balance of the total force F exerted by the top lens across the liquid and of the spring restoring force. For our system, we model the total force by the superposition of two forces: an equilibrium, electrostatic force, and a dynamic, viscous force. A strong assumption is that these two contributions are supposed to be totally independent, i.e. no electrokinetic effect is included.

To describe the viscous force, we suppose that the nanoconfined IL behaves like a charge-less, Newtonian liquid in an incompressible flow between undeformable solids. In the lubrication approximation (D R), the viscous force is simply given by Reynolds formula. As suggested by Chan & Horn, we use a modified version of the Reynolds formula, with hydrodynamic zero shifted by the distance b in the liquid (see the viscous term in equation 1). This "negative slip length" can be seen as the presence of a "stopped layer" on each solid surface, i.e. liquid films of thickness b/2 that do not contribute to the drainage flow, on top of which a no-slip boundary condition applies.

The electrostatic interaction is modelled by an exponential function of decay length λ s and amplitude 2πRA, where 2πR is a geometrical factor coming from the Derjaguin approximation. Such an exponential law is expected to hold in the context of the linear Debye-Hückel theory for diluted electrolytes while imposing a constant potential boundary condition between two symmetric planar surfaces [START_REF] Israelachvili | Intermolecular and Surface Forces[END_REF], which is clearly not the current situation. Although the origin of the exponentially decaying force in ionic liquids and concentrated electrolytes is still under discussion, it appears clear that experiments to date agree on the empirical functional form [START_REF] Gebbie | Long range electrostatic forces in ionic liquids[END_REF][START_REF] Lee | Underscreening in concentrated electrolytes[END_REF][START_REF] Lee | Scaling Analysis of the Screening Length in Concentrated Electrolytes[END_REF]. Thus the equation of motion finally reads:

k (D(t) -D 0 + vt) = F (t) = - 6πηR 2 Ḋ D -b + 2πRAe -D/λs . (1)
Note that the structural force, observed experimentally for this system (Figure 1(b)), as well as the Van der Waals interaction are not included in the model for the sake of simplicity. Each of these is short ranged compared to the viscous and electrostatic interactions in pure ionic liquid and neglecting them does not alter the fitting of data in this

work.

An analytic solution only exists in the purely viscous case (A = 0) [START_REF] Chan | The drainage of thin liquid films between solid surfaces[END_REF], therefore the equation of motion is integrated numerically, using the Runge-Kutta method implemented in C++. Four parameters are known: the radius R and the initial distances D 0 are measured, the spring constant k has been calibrated prior to the experiment and the viscosity η is known from literature. The four other parameters are manually tuned to fit the data: first the velocity v is fitted at large distances, secondly the hydrodynamic zero b, the screening length λ s and the amplitude A are fitted at small scales, and finally all the parameters are refined iteratively. For a single trajectory, the parameters b, λ s and A are interdependent, but as we will see later the exploration of 3 decades of velocities allows to disentangle the role of each parameter. In short, high velocities provide first the determination of the hydrodynamic zero b only, from there the screening length λ s and the amplitude A can be fitted unequivocally at very low velocities.

After this fitting procedure, when the imposed motion D 0 -vt and so the spring deflection are precisely known, the total force can be calculated and separated into viscous and electrostatic contributions.

For the sake of clarity, only data obtained at four different approach velocities are shown in the following section. In practice, very many approaches have been carried out experimentally to ensure the reproducibility of the observations and the robustness of the interpretations.

III. RESULTS AND DISCUSSION

The force profile obtained upon approach at the lowest velocity v = 0.075 ± 0.010 nm/s is shown in Figure 1(b).

Here the rescaled force F/R is used, because according to the Derjaguin approximation this quantity is proportional to the interaction energy acting between the same system in the plane-plane geometry, and is convenient to compare measurements performed in different geometries.

At short distances (D < 10 nm typically), one observes a so-called structural force, with characteristic steps. This phenomenon is known to be the signature of the ordering of the molecules in layers in the confined film. A jump-in corresponds to the sudden squeeze-out of a layer, which -for ionic liquids -must consist of equal numbers of positive and negative ions to preserve electroneutrality. Here we show only the repulsive steps which contribute to the overall oscillatory structural force law, and we only probe the range down to 4nm: the present experiment was optimised for weak and long range forces (using a piezo with short range). Other studies have probed the full oscillations in the past [START_REF] Smith | Switching the Structural Force in Ionic Liquid-Solvent Mixtures by Varying Composition[END_REF][START_REF] Smith | Influence of Lithium Solutes on Double-Layer Structure of Ionic Liquids[END_REF]. Our present measurements are consistent with the previous observations in terms of the positions and amplitude at the maxima. The closest (small D) maximum explored here, indicated by the green dashed line in Figure 1(b), corresponds to the sixth layer and -with our now-improved resolution as described above -we are able to probe a further seven layers. That is to say, steps in the force profile corresponding to molecular layers (structural forces) are measured out to 13 layers in total. The remaining imperfections in the signal are horizontal "spikes", that are due to mechanical vibrations that are not completely damped. In Figure 1(d) the location of the maxima D i are plotted against the layer index i. The linearity of the curve obtained is very good, illustrating the regularity of the layers. The slope is the average layer thickness, and is found to be 0.71 ± 0.01 nm. This is close to a previous AFM determination and to the predicted ion pair diameter of 0.75 ± 0.01 nm [START_REF] Hayes | At the interface: solvation and designing ionic liquids[END_REF]. Our value is also slightly smaller than the 0.78 ± 0.01 nm previously reported [START_REF] Smith | Influence of Lithium Solutes on Double-Layer Structure of Ionic Liquids[END_REF], although in that case the layer thickness was extracted from the locations of minima and from deepest layers.

At larger distances (D > 10 nm typically) we also observe a long-range force, with no distinguishable layers, decaying much slower than the structural force. This monotonic force will be our focus for the remainder, in which we consider its equilibrium and dynamic contributions.

A. High velocities: characterization of the dynamical response while the top lens is moving at v = 95 ± 5 nm/s, the relative velocity Ḋ between the lenses does not exceed 12 nm/s.

Because of this huge viscous force, the piezo-electric tube reaches its maximum range and stops before the contact between the two lenses, at a time indicated by the vertical gray line. Nevertheless, the approach towards contact continues after, because of the elastic energy stored in the spring. The stop of the piezo motion at a known time is taken into account when solving the equation of motion, so that this experimental detail is not a problem to interpret Interestingly, we distinguish only three layers, which locations correspond to the three deepest layer over the nine layers seen at v = 0.075 ± 0.010 nm/s (Figure 1(b)). Such decrease of the range of the structural force with the approach velocity has already been observed for apolar molecules [START_REF] Bureau | Drainage of a nanoconfined simple fluid: Rate effects on squeeze-out dynamics[END_REF][START_REF] Bureau | Rate Effects on Layering of a Confined Linear Alkane[END_REF]. To interpret this phenomenon, the authors proposed a scenario in which the rapidly squeeze-out film is quenched in a metastable disordered state. This aspect is beyond the scope of the present work but it is of interest to note that this phenomenon can also be observed with ionic liquids.

Our determination of a negative slip length of [C 2 C 1 Im][NTf 2 ] on mica is consistent with recent measurements made by Garcia et al. [START_REF] Garcia | Nanomechanics of ionic liquids at dielectric and metallic interfaces[END_REF]. They use a modified SFA, which is able to superpose a tiny oscillation to the steady motion, and so allows to distinguish the conservative (in-phase) and dissipative (out-of-phase) responses. For

[C 4 C 1 Im][PF 6 ]
between Pyrex or Platinum surfaces, they concluded that the steady flow obeys classical hydrodynamics with a no-slip boundary condition on top of a few nanometers thick film of ordered molecules. Intriguingly, they found that the oscillating flow simultaneously obeys to classical hydrodynamics with no-slip boundary condition close to the solid surface. This observation suggests that the ordered layers of molecules can participate or not to the flow, depending on the nature of the excitation.

Furthermore, it is worth noting that the extracted value of b = 4.2±0.7 nm seems to match with the liquid thickness below which a stick-slip, solid-like response is observed in friction measurements [START_REF] Perkin | Layering and shear properties of an ionic liquid, 1-ethyl-3-methylimidazolium ethylsulfate, confined to nano-films between mica surfaces[END_REF][START_REF] Smith | Quantized friction across ionic liquid thin films[END_REF]. In these experiments, the top lens is moved laterally and the liquid is sheared at strains comparable with our steady approaches in the normal direction. These two independent observations are thus consistent with the picture of ordered layers of molecules that don't flow when subjected to a large scale excitation.

As we will see in section III B, our system clearly exhibits an electrostatic interaction at very low velocity. Anticipating this result, we analyzed the data at high velocity while incorporating the static contribution determined at very low velocity (i.e. with A = 0.22 mN/m and λ s = 8.0 nm). At such high velocity, the electrostatic force is much smaller than the total force (see Figure 2(b), and considering a purely viscous response is a very good approximation. However, given we know that there is an electrostatic force, we could have expected some electrokinetic effects, especially when probing confinement within the range λ s = 8.0 nm of the electrostatic force. In the end we found no measurable deviation from classical hydrodynamics, considering the IL as a neutral, Newtonian liquid of constant, bulk viscosity.

Finally, one can ask why we do not observe elastohydrodynamic effects [START_REF] Leroy | Hydrodynamic interactions for the measurement of thin film elastic properties[END_REF][START_REF] Leroy | Hydrodynamic Interaction between a Spherical Particle and an Elastic Surface: A Gentle Probe for Soft Thin Films[END_REF][START_REF] Villey | Effect of Surface Elasticity on the Rheology of Nanometric Liquids[END_REF][START_REF] Wang | Out-of-Contact Elastohydrodynamic Deformation due to Lubrication Forces[END_REF][START_REF] Wang | Elastic deformation during dynamic force measurements in viscous fluids[END_REF][START_REF] Wang | Elastic deformation of soft coatings due to lubrication forces[END_REF]. Indeed, although the model assumes undeformable solids, the materials used experimentally are not infinitely stiff. We have two symmetric lenses which are stacks composed of 3.6083 µm of mica / ∼ 50 nm of silver / ∼ 10 µm of EPON supported by a centimetric glass cylinder. We can estimate the order of magnitude of the indentation δ of the solid, by simply considering a rigid sphere approaching a soft, semi-infinite medium. For undeformable solids and the simple case of a no-slip boundary condition, the sphere exerts the Reynolds force 6πηR 2 Ḋ/D on a radius of action roughly given by the range √ RD of the pressure distribution. Taking these scaling as first estimates for the soft case, we can write Hooke's law [START_REF] Leroy | Hydrodynamic interactions for the measurement of thin film elastic properties[END_REF]:

δ √ RD ∼ 1 E 6πηR 2 Ḋ/D πRD , ( 2
)
where E is the effective elastic modulus, given by E/(1 -ν 2 ) for a semi-infinite medium, where E is the Young's modulus and ν the Poisson's ratio. As expected, the indentation increases with the viscosity, the velocity and the softness of the substrate. The indentation also increases when the confinement increases, and we can define a typical size D c such that δ ∼ D. This lengthscale corresponds to a situation where the actual liquid film thickness is doubled in comparison to the rigid case, and provides a scale above which elasticity cannot be neglected. From equation 2 we finally get the scaling:

D c ∼ ηR 3/2 Ḋ E 2/5
(3)

Our lenses are not semi-infinite, but are stratified with typical Young's moduli of ∼ 200 GPa for the mica, ∼ 50 GPa for the glass or the silver, 1 GPa for the EPON [START_REF] Wang | Out-of-Contact Elastohydrodynamic Deformation due to Lubrication Forces[END_REF]. In first approximation, we take the effective modulus E ∼ 10 GPa extracted from JKR tests on similar systems [START_REF] Horn | Measurement of the Deformation and Adhesion of Solids in Contact[END_REF], even if it has been shown that effective considerations in contact deformation cannot be extended to elastohydrodynamic deformation [START_REF] Wang | Out-of-Contact Elastohydrodynamic Deformation due to Lubrication Forces[END_REF][START_REF] Wang | Elastic deformation of soft coatings due to lubrication forces[END_REF]. For the maximum relative velocity Ḋ ∼ 10 nm/s reached in our experiment, we finally obtain D c ∼ 1 nm. Qualitatively, this very rough estimation means that elastic deformation could play a role when the two surfaces are closer than a few nm. Additional observations support the scenario of negligible elastohydrodynamics. First, no significant deformation of the lenses is observed when looking at the FECO pattern. Secondly, similar experiments have been performed with the same system and parameters, except that the liquid is a PDMS polymer melt [START_REF] Horn | Molecular organization and viscosity of a thin film of molten polymer between two surfaces as probed by force measurements[END_REF][START_REF] Horn | A reevaluation of forces measured across thin polymer films: Nonequilibrium and pinning effects[END_REF], and reported no significant elastohydrodynamic effect. Finally, we reproduced these measurements in our set-up with a similar PDMS polymer melt (Xiameter PMX 200 Silicone Fluid 50 cSt), leading to the same conclusion. In Figure 3(a) are shown in black the measurements performed at a very low velocity. The hydrodynamic zero is kept to its value b = 4.2 nm extracted before, but in this regime the trajectory cannot be fitted with a purely viscous force (see blue line which corresponds to A = 0). At intermediate distances (7 nm D 25 nm), a significant repulsion is not captured by a purely viscous model. On the contrary, the trajectory can be very well-fitted when considering a static, exponentially decaying contribution, interpreted as an electrostatic repulsion (see the orange curve). The best fit is obtained with a velocity v = 0.24 ± 0.01 nm/s (see the orange dashed line), a decay length λ s = 8.0 ± 0.5 nm and an amplitude A = 0.22 ± 0.02 mN/m. The trajectories computed for three different values of b ∈ {3.5 nm; 4.2 nm; 4.9 nm} (orange curves in Figure 3(a) and zoom in inset) show that the hydrodynamic zero b modifies very weakly the parameters fitted for the static contribution.

The same data are presented in terms of force in Figure 3(b) on a semi-log scale. At large distances (D 25 nm), the main contribution is the viscous force. The fit in this region is not perfect, as the semi-log scales emphasizes the imperfections due to thermal drifts. At intermediate distances (7 nm D 25 nm), the total force varies over one order of magnitude with a shape compatible with an exponential behavior. The electrostatic force is the dominant contribution, as the viscous force is one order of magnitude smaller. Additionally, the shape of the viscous force is significantly different than the total force, and it is clear that it is not possible to fit the data with a purely viscous force, for example by playing with the hydrodynamic zero b that basically translates the curve along the horizontal axis. Lastly, the region of very strong confinement (D 7 nm) is not fitted by the model; unsurprising because the model does not include the structural or Van der Waals forces.

To have a good determination of the static contribution, it is important to highlight that the initial distance should be large enough: D 0 λ s . If not, the values of v and λ s cannot be determined independently from the fit. This is a well-known experimental requirement, which has to be fulfilled even when neglecting the viscous effects.

Finally, the screening length λ s = 8.0 ± 0.5 nm obtained here is in reasonable agreement with the values previously reported for this system: 6.7 ± 0.1 nm [28] and 7.1 ± 0.7 nm [START_REF] Smith | The Electrostatic Screening Length in Concentrated Electrolytes Increases with Concentration[END_REF]. This might at first appear surprising given that, in the earlier studies, the viscous contribution was not formally taken into account. In the next section we rationalise why the electrostatic decay lengths match nonetheless.

C. Intermediate velocities: why ignoring viscosity can be taking it into account The gray dashed line is a linear fit a large distances, which doesn't take into account the strong curvature observed initially (zoom in inset) and provides the velocity 0.99 ± 0.04 nm/s. (b) Same measurements, but with data presented in terms of force, in semi-log scales. The black curve is obtained from the difference between the black curve and the orange dashed line in (a).

The total force (in orange), viscous contribution (in blue) and electrostatic contribution (in red) are deduced from the fit with non-zero A, whereas the force profile in gray line is computed from the gray dashed line in (a). Inset: same data plotted in linear scales. The vertical dashed line in gray indicates the intersection between the viscous and electrostatic contributions, and the second gray dashed line is a guide for the eyes to illustrate the linearity of the viscous component at large distances.

We now return to the analyses made in previous studies [START_REF] Gebbie | Ionic liquids behave as dilute electrolyte solutions[END_REF][START_REF] Espinosa-Marzal | Ionic Liquids Confined in Hydrophilic Nanocontacts: Structure and Lubricity in the Presence of Water[END_REF][START_REF] Cheng | Characterizing the Influence of Water on Charging and Layering at Electrified Ionic-Liquid/Solid Interfaces[END_REF][START_REF] Gebbie | Long-range electrostatic screening in ionic liquids[END_REF][START_REF] Smith | The Electrostatic Screening Length in Concentrated Electrolytes Increases with Concentration[END_REF][START_REF] Espinosa-Marzal | Microslips to "Avalanches" in Confined, Molecular Layers of Ionic Liquids[END_REF]. In those studies, comparable measurements were performed with ionic liquids of similar or higher viscosity and using velocities close to 1 nm/s. Viscous forces were not explicitly included in fitting the data; instead, in those experiments, the velocity was varied over a small range to ensure the force and decay was not velocity-dependent. This means that the motion imposed on the top lens (i.e the piezo velocity) was simply deduced from a linear fit of the trajectory at large distances. The previous fitting procedure appears particularly inappropriate when looking at the starting curvature of the trajectory, which we are now able to resolve clearly (inset of Figure 4(a)), which is a signature of the presence of a non-negligible viscous contribution. Nonetheless, we now make such a linear fit with our present data (dashed line in gray), to quantify the error committed when not taking the viscous contribution into account. The fitted velocity is v = 0.99 ± 0.04 nm/s, slightly smaller but close to the velocity extracted with our full model. A force is then deduced from the difference between the trajectory and this base line (gray line in Figure 4(b) and is finally fitted with a simple exponential law to extract a screening length. As we see in Figure 4(b) it turns out that this ill-defined force is well superimposed on the electrostatic contribution obtained from our model (red curve). As a result, the same value for the screening length is given by the two methods. This surprising matching of the electrostatic screening length, both with and without accounting for the viscous contribution, can be interpreted qualitatively by inspecting the two contributions to the total force on linear scales (inset of Figure 4(b). At distances larger than ∼ 11 nm (indicated by a vertical dashed line in gray), the viscous force is linear to very good approximation (except, of course, for the curvature at D 60 nm). At distances shorter than ∼ 11 nm, the electrostatic part is bigger than the viscous one and varies more slowly, thus the departure from linearity of the viscous force is not substantial. We find, therefore, that calculating a force from a direct linear fit of the trajectory is an uncontrolled way to subtract the viscous contribution to the interaction. Neglecting the viscous effects turns out -in this and similar cases -to take them into account.

One can ask why the viscous force, that depends on Ḋ/D, is in fact found to be linear with D. To discuss this point, we solved numerically the equation of motion for four approach velocities v ∈ {0.1 nm/s; 1 nm/s; 10 nm/s; 100 nm/s} (the other parameters being kept constant as b = 4.2 nm, λ s = 7.5 ± 0.5 nm and A = 0.18 ± 0.02 mN/m), and we computed the corresponding electrostatic and viscous forces (Figure 5(a)). To compare the dissipative forces over 3 decades of approach velocities, it is convenient to plot the ratio of the viscous force F hydro to the approach velocity v, that can be interpreted as a friction coefficient (Figure 5(b)). For all the trajectories, i.e. for the complete range of velocities explored experimentally, we observe a significant curvature at large distances caused by the sudden start-up of the piezo, a strong curvature at short distances due to the close vicinity of the solid surfaces, and a reasonably linear regime at intermediate distances. Such quasi-linear regime at intermediate distances is not intrinsic to the mathematical form of the viscous force, but results from the competition between the expansion of the piezo on the top lens that tends to accelerate the squeeze-out of the liquid and the compression of the spring on the bottom lens that tends to decelerate it. Interestingly, the "friction coefficient" is found to decrease with the velocity. This is because the higher the velocity, the greater the difference between the top lens velocity v and the maximum relative velocity Ḋ at intermediate distances. The two surfaces are in contact before a steady-state at Ḋ = v can be reached, thus the viscous force F hydro increases slower than the velocity v. Discrepancy from linearity is clearly visible for the lowest velocity, and is due to the electrostatic contribution that becomes dominant for D 30 nm (Figure 5(a)). The effect of viscosity is tiny at very low velocity, and it is of course this regime that should be preferred to have the best determination of the static contribution.

IV. CONCLUSION

In this paper we tackled the question of the contribution of the viscous force to the long-range interactions previously observed for nanoconfined ionic liquids. Using [C 2 C 1 Im][NTf 2 ] confined between mica surfaces in a SFB, we explored the effect of approach velocity over three orders of magnitude and down to 0.07 nm/s, and we systematically analyzed the data by solving numerically an equation of motion that combines viscous and electrostatic forces. This systematic study allowed us to disentangle the dynamic and static contributions to the interaction, leading to the following conclusions. First, no significant deviation from classical hydrodynamics is found for this nanoconfined ionic liquid at high velocity, and an excellent agreement is achieved when incorporating a negative slip boundary condition with a negative slip length b = 4.2 ± 0.7 nm. In other words, a nanometric film of liquid in the vicinity of the solid surface does not participate in the flow. We note that this shift of the hydrodynamic zero is in the range of the structural force and corresponds to the confinement at which stick-slip behavior is typically observed in friction experiments.

All these independent observations point toward the presence of ordered, solid-like layers of liquid molecules in the vicinity of the solid surface. Then, we showed that an equilibrium, exponentially decaying force has to be invoked in order to describe quantitatively the behavior at low velocity. The screening length λ s = 8.0 ± 0.5 nm extracted from our fitting procedure is in agreement with the values previously reported for this system. This appears surprising given that in previous studies the viscous force was not taken into account explicitly. We interpreted this surprising matching by looking at an approach at intermediate velocity, comparable to the velocities used in previous studies.

We showed that analyzing the data as if there is no viscous force is an (uncontrolled) way to subtract the viscous contribution to the total force, and thus led to a correct determination of the screening length.

Such long-range electrostatic interactions have also been reported using colloid probe AFM [START_REF] Hjalmarsson | Switchable long-range double layer force observed in a protic ionic liquid[END_REF] . Similarly, one can ask what is the contribution of viscous forces in these observations. The ratio of the viscous force to the electrostatic force expressed in this paper varies as R Ḋ. SFB typically uses a radius ∼ 1 cm and a velocity ∼ 1 nm/s, whereas colloid probe AFM typically uses a radius ∼ 10 µm and a velocity ∼ 100 nm/s. For a given system, the ratio of the viscous force to the electrostatic force for colloid probe AFM is typically ten times smaller than for SFB. Viscous effects are thus expected to have less impact in colloid probe AFM experiments although the difference is not large and similar viscous effects will often be significant in AFM studies. In general it is necessary to consider, for any scenario involving viscous drainage in combination with equilibrium surface forces, the contribution of each to the total measured force. In our study, all the trajectories recorded while exploring 3 orders of magnitude of approach velocity have been fitted quantitatively and consistently by a model that simply combines independent viscous and electrostatic forces. With the current parameters it was not found necessary to consider electrokinetic effects. It would be interesting to investigate more the possibility of such phenomena, for example by exploring higher velocities or by simultaneously applying electric fields across the liquid gap.

FIG. 1

 1 FIG. 1. a) Schematic showing the principle of SFB measurements. b) Force profile measured during approach of two mica surfaces separated by [C2C1Im][NTf2] at v = 0.075±0.010 nm/s. At short distances, below ca. 7-8 nm one observes a structural force with characteristic steps due to successive squeeze-out of ordered layers (also called solvation force). The configuration with i = 6 cations along the gap is indicated by the green dashed line and is illustrated in the inset. At larger distances a long-range force is found, and the aim of this paper is to determine the static and dynamic contributions to it. c) Chemical structure of [C2C1Im][NTf2]. d) Location of the layers in the structural force region: the slope gives the average layer thickness.

Figure 2 (

 2 Figure2(a) shows in black the trajectory recorded when the top lens is moved at a high velocity. At such high velocity the viscous force is expected to dominate and the data could be fitted by assuming a purely viscous force (A = 0), with v = 95 ± 5 nm/s and b = 4.2 ± 0.7 nm. In this regime, the viscous damping is very clearly significant:

FIG. 2 .

 2 FIG. 2. (a)Trajectory recorded when the top lens is moving at v = 95 ± 5 nm/s (piezo motion indicated by the dashed line in orange). The black curve is the experimental signal, the orange (resp. green) curve is the fit by the model with b = 4.2 ± 0.7 nm (resp. b = 0), λs = 8.0 nm and A = 0.22 mN/m. These values of λs and A are justified by fits at lower velocity as described later; as noted in the text a fit with A = 0 would be indistinguishable at this velocity. (b) Same measurements, but with data presented in terms of force, in semi-log scales. The black curve is obtained from the difference between the black curve and the orange dashed line in (a). The total force (in orange), viscous contribution (in blue) and electrostatic contribution (in red) are deduced from the fit. In both graphs, the vertical line in gray indicates the time and distance at which the piezo reaches its maximum deformation and the top lens stops.

FIG. 3 .

 3 FIG. 3. (a)Trajectory recorded when the top lens is moving at v = 0.24 ± 0.01 nm/s (piezo motion indicated by the dashed line in orange). The black curve is the experimental signal, the blue curve is an attempt to with a single, viscous force (A = 0), and the orange curve is the fit by the model with λs = 8.0 ± 0.5 nm and A = 0.22 ± 0.02 mN/m and three different values b ∈ {3.5 nm; 4.2 nm; 4.9 nm} (zoom in inset). (b) Same measurements, but with data presented in terms of force, in semi-log scales. The black curve is obtained from the difference between the black curve and the orange dashed line in (a). The total force (in orange), viscous contribution (when fitting with non-zero A; in blue) and electrostatic contribution (in red) are deduced from the fit.

10

 10

Figure 4 (

 4 Figure 4(a) shows the data obtained when moving the top lens at an intermediate velocity. Again, it is not possible to fit the trajectory at intermediate distances when considering a purely viscous force with b = 4.2 nm (blue curve, corresponding to A = 0), even if the discrepancy is less pronounced than at very low velocity. A fit can be achieved when including an electrostatic contribution (orange curve), with v = 1.13 ± 0.04 nm/s and parameters for the static force which are consistent with the values extracted at very low velocities: λ s = 7.5 ± 0.5 nm and A = 0.18 ± 0.02 mN/m. In this regime, the viscous force is never negligible, whatever the distance considered. At intermediate distances, the viscous and electrostatic forces have the same orders of magnitude and contribute equally to the total force (Figure 4(b)).

FIG. 4 .

 4 FIG. 4. (a)Trajectory recorded when the top lens is moving at v = 1.13 ± 0.04 nm/s (piezo motion indicated by the dashed line in orange). The black curve is the experimental signal, the blue curve is an attempt to fit with a single, viscous force (A = 0) with b = 4.2 nm, and the orange curve is the fit by the model with λs = 7.5 ± 0.5 nm and A = 0.18 ± 0.02 mN/m.

  FIG. 5. (a) Electrostatic and viscous forces calculated numerically for four approaches at v ∈ {0.1 nm/s; 1 nm/s; 10 nm/s; 100 nm/s}, with b = 4.2 nm, λs = 7.5 ± 0.5 nm and A = 0.18 ± 0.02 mN/m. (b) Ratioof the viscous force to the approach velocity as a function of the distance between the surfaces, for the same trajectories.

ACKNOWLEDGMENTS

We acknowledge Joelle Frechette and Yumo Wang for helping us to perform a comparative measurement with a PDMS melt. S.P. and R.L are supported The Leverhulme Trust (RPG-2015-328) and the ERC (under Starting Grant No. 676861, LIQUISWITCH). R.L is supported by the EPA Cephalosporin Junior Research Fellowship and Linacre College (University of Oxford). SP is grateful for research leave enabled by the Philip Leverhulme Prize.


Gineth Saracco

Mohamed Rouai

Konstantinos Chalikakis

Abdelilah Dekayir

Helene Miche

Adriano Mayer

Khaoula Qarqori

Brahim Abderahim

Christophe Emblanch

Geophysical constraints on the groundwater circulation and its renewal at the main Bittit spring, supplying Meknes and its region (Sais Basin

des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.