Electrical charging effects on the sliding friction of a model nano-confined ionic liquid
Recent measurements suggest the possibility to exploit ionic liquids (ILs) as smart lubricants for nano-contacts, tuning their tribological and rheological properties by charging the sliding interfaces. Following our earlier theoretical study of charging effects on nanoscale confinement and squeezout of a model IL, we present here molecular dynamics simulations of the frictional and lubrication properties of that model under charging conditions. First we describe the case when two equally charged plates slide while being held together to a confinement distance of a few molecular layers. The shear sliding stress is found to rise strongly and discontinuously as the number of IL layers decreases stepwise. However the shear stress shows, within each given number of layers, only a weak dependence upon the precise value of the normal load, a result in agreement with data extracted from recent experiments. We subsequently describe the case of opposite charging of the sliding plates, and follow the shear stress when the charging is slowly and adiabatically reversed in the course of time, under fixed load. Despite the fixed load, the number and structure of the confined IL layers changes with changing charge, and that in turn drives strong friction variations. The latter involve first of all charging-induced freezing of the IL film, followed by a discharging-induced melting, both made possible by the nanoscale confinement. Another mechanism for charging-induced frictional changes is a shift of the plane of maximum shear from mid-film to the plate-film interface, and viceversa. While these occurrences and results invariably depend upon the parameters of the model IL and upon its specific interaction with the plates, the present study helps identifying a variety of possible behavior, obtained under very simple assumptions, while connecting it to an underlying equilibrium thermodynamics picture.
Ionic liquids (ILs) – organic salts that are liquid at room temperature – are
of considerable and increasing physical and technological interest.
ILs are characterized by negligible vapor pressures, high temperature stability and a wide
electrochemical window plechkova08. Moreover their physical properties can
be widely tuned by changing the molecular structure of the cation-anion pairs hallett11.
Many ILs strongly adhere to solid surfaces hayes10 and can form wear-protective
films capable of resisting much higher loads than those of molecular lubricants perkin10.
For this reason among others their use as boundary lubricants has been recently pursued in experimental
studies under nanoscale confinement, employing both the surface force apparatus (SFA) perkin10; spencer,
and atomic force microscopy (AFM) werzer12; li13; bennewitz; li14; hoth14, as well in
in a number of computer simulations at various levels of idealization padua13; federici14; fajardo15.
Our specific interest focuses here on the exploration of the frictional changes induced by electrical
charging of the plates which confine the IL, a subject only partially covered by previous theoretical
works kornyshev; fajardo15; capozza2015.
Besides lubrication, the dynamical behavior of ILs under charging is of additional interest in the field
of supercapacitors kim11; shim10.
Experimentally, there is a number of nano-frictional studies of
surfaces charged with respect to a reference electrode bennewitz; black13; li14. Not so
between two oppositely charged electrodes. Recent macroscopic friction studies have begun to appear for
IL sliding under charging conditions dold2015; yang14. Here however we shall restrict ourselves to
friction under conditions of nanoscale confinement.
At the molecular level, ILs confined between hard plates become structured in the form of layers that run parallel to the plates, not unlike other liquids but further characterized by the charge order typical of molten salts, with an alternation of positive and negative ion layers, and an interlayer separation that corresponds to the ion pair size perkin10; perkin; werzer12; hoth14; atkin07; bou-malham10; ueno10; perkin11; perkin12; li13b. Plate charging is expected to give rise to rearrangements of the IL layering and ordering, at least near the confining plate-IL interfaces, with a consequent and probable change of the lubrication properties under plate sliding. Charging-induced friction changes have already been to some extent detected and described experimentally bennewitz; li14. Not unexpectedly, they are reported to depend on the particular kind of IL and of confining surfaces chosen. For example, a strong reduction of the friction of an AFM silica colloid probe was found by negatively biasing a Au(111) surface bennewitz immersed in . However, a similarly strong friction drop occurred at positive bias upon sliding a sharp AFM tip on highly oriented pyrolytic graphite immersed in a IL li14.
This diversity of behavior, reflecting the variety of ILs and the different nature of the confining and sliding surfaces, suggests some flexibility in the corresponding theoretical modeling, aimed at a broad exploration of how some of this variety of behavior could be addressed with a handful of model parameters, rather than describing in details one or another particular case. The minimal IL model and a natural first choice is a simple molten salt, such as liquid NaCl. Fedorov and Kornyshev kornyshev did that, their model consisting of two spherical Lennard-Jones (LJ) particles with unequal radii further endowed with opposite charges. The frictional behavior of that model was recently investigated by Fajardo et al. fajardo15 who noted interesting similarities with data by Li et al. li14. However, a parallel study which we conducted of this type of model under simple confinement and squeezout capozza2015 revealed that the confined liquid film evolved far too sharply from a layered liquid for neutral plates at large spacings, to a strongly crystallized rocksalt structure solid under narrow confinement and/or plate charging. On top of that, the simple molten salt model wetted the plates much too readily, even given a reasonable choice of parameters. Both of these aspects, excessive crystallization and excessive wetting tendencies, represent a severe oversimplification by comparison with real IL lubricants whose wetting capability is generally more modest, and where, owing to far more complex molecular structures, the layering and solidification which occurs under strong confinement does not usually imply cystallization, and is replaced by glass-like disordered structures. That complex behavior is of course much better captured by fully realistic IL model simulations, of which there are several good examples in literature, padua13; federici14 but whose difficulty and intrinsic complexity make the studies much more work-intensive while at the same time emphasizing aspects that are to some extent undesirably specific to each case.
A previous study of our own IL model capozza2015 was designed to bridge the gap between these two extremes, either oversimplification
or overspecificity. We showed that a minimal modification of the charged LJ model, consisting of a neutral ”tail” rigidly
attached to one of the ions (the cation in our case, but the anion could have been equivalently chosen) could go a long
way in our desired direction, making the structural and squeezout properties
of the confined IL model considerably more realistic while still not specific. While the anion and the charged part of the
cation retain a strong Coulomb correlation, the neutral tail is generally and merely pushed out of the way, resulting
in poor tail-tail correlations. The tail-induced disorder, along with the steric hindrance they introduce, impede
overall rocksalt crystallization with drastic changes in the ILs solidification, also introducing new unexpected
structures as well as the desired glassy aspects and wettability features. These structural features were found to evolve and change under plate
charging, foreshadowing a variability that could be explored through a variation of the model parameters. At the same time, this
IL ”tailed model” (TM), to be further detailed below, is simple enough to permit a novel, quantitative calculation
and characterization of the confined IL enthalpy versus interplate distance, and – of crucial importance here – its variations
under plate charging. The minima of the enthalpy-distance curves identify stable or metastable layering states of the confined IL,
showing a perfect fit with the evolution in the structural layering oscillations of density and charge observed
in simulations of the confined IL under variable interplate separation and plate charge.
The layering transitions during squeezout, well established in real IL experiments atkin07; smith13; spencer, are
reproduced by the TM model and directly related to jumps between enthalpy minima, while these minima evolve
and cross values under closing of the gap between plates, or under variable plate charging capozza2015.
With that very useful characterization, the TM model provides a ready tool for a fresh study of some generic effects of
electrical plate charging on the nanoscale sliding friction of confined ILs.
In this work we first introduce and detail the TM model and the chosen confining geometry. We then analyze the sliding simulations
and discuss how the structural changes, induced by a negative plate charging, modify the frictional properties.
Interesting observations emerge examining the shear velocity profile inside the IL film as a function of the plate charging.
First, we observe a switching of the shear band plane from the plate-film interface, typical of low
friction, to the mid-film region, typical of larger friction. A second important effect is the
charging-induced solidification/melting process described in our previous, static study, and typically
ignored in literature. Charged plates generally increase the IL wetting, giving rise to ”electrowetting”.
But, as most ILs are not far above freezing at room temperature,
increased wetting will in turn encourage solidification, the capability to withstand load, with a large effect on friction.
Conversely, plate neutrality and partial dewetting facilitates melting of the IL film, leading to squeezout under even
The final and most specific part of this work is the investigation of the effects on friction of opposite charging (equivalent to an applied interplate voltage) of the two confining plates. A given charge state of the plates is accompanied by a strong structuring of the IL near and between the plates. A reversal of plate charges is followed by an accompanying destructuring and even melting of that layering structure, prior to the reversal of local IL ordering and eventual re-solidification. The sliding friction between plates undergoes a corresponding charge-dependent evolution as a consequence. This kind of phenomena, only poorly explored so far, constitute a fresh prediction of this study.
Ii Model, simulation details, system preparation
To pursue our desired broad scenario of possible charging effects on the lubrication properties of nano-confined ILs we
adopt the simple TM model capozza2015 which deliberately neglects the fine structure of the ions retaining
only minimal features common to most of the ILs,
such as the anion-cation asymmetry, their possible amphiphilic character and their generally irregular shape.
The TM model is a two-component charged Lennard-Jones (LJ) system where anions and cations
not only have different radii but where a tail is attached to the cation. The cation is therefore a dimer consisting of
a positively charged, small-sized LJ head, rigidly bound to an equal size, neutral particle, see Fig. 1.
The effect of the tail upon the ILs wetting and confinement properties is important, as mentioned in the
introduction and recently shown in ref. capozza2015, where the details about the adopted
potentials and the corresponding parameters are also specified and discussed. Essentially the tails give the
molecules a larger gyration ratio in the liquid state, and
prevent low temperature ordering replacing crystallization with a glass-like solidification – all features that are
common and important in real ILs.
All the molecular dynamics (MD) simulations were performed using the LAMMPS code lammps.
The long range Coulomb interactions have been treated in reciprocal space using a particle-particle particle-mesh solver.
Canonical ensemble configurations were generated by means of a Langevin thermostat
directly applied to the lubricant molecules. The plates were treated as rigid bodies,
the lower one fixed and the upper one subjected to a directed load force and driven to slide along (x,y).
Their structure consists of rigid close-packed triangular lattices of LJ
particles with nearest neighbor spacing 0.52 nm, mimicking the periodicity of a mica surface. A graphite, or graphene
plate, which unlike mica may be electrically charged li14 has a
periodicity of about 0.246 nm, which is not too different from half that value.
In most of our simulations the liquid is confined between
two identical plates with a modest surface charge density
, subsequently increased to .
As a reference, the surface charge density quoted for mica in contact with HO is -33.4 ,
a much larger value northern91.
The state of charge of a mica surface in contact with an IL will certainly differ,
but it is still generally believed that some charging remains.
Since nature as well as the algorithm used to treat the long range interaction in reciprocal space require the overall
system to be neutral, when endowing both plates with a total net nonzero charge we correspondingly
add extra ions of the opposite charge to compensate.
The lateral drive is actuated through a spring connected to the top plate and pulled at constant velocity . The same value was adopted throughout unless another value is explicitly stated. The instantaneous frictional force opposing the motion is measured through the elongation of the driving spring as , where is the coordinate of the top plate center of mass. The relevant quantity in friction is eventually the shear stress where is the area of the contact. In our case nm is a constant, so presenting or is just the same. In experiments however the contact area is not generally constant, so that the shear stress must be extracted before making comparisons.
The open geometry described in Fig. 1 has been chosen because it permits particles to be squeezed out or sucked in from two lateral IL droplets. These droplets thus serve as liquid reservoirs allowing the number of ions effectively confined inside the gap to adjust and change dynamically depending on the loading, charging, and sliding conditions, realizing a kind of effectively grand-canonical configuration. The distinguishing feature of the open boundary geometry employed here, as opposed to several previous studies which used closed periodic boundary conditions, is to address a situation much closer to real SFA or AFM experiments, permitting in particular squeezeout transitions and the consequent transverse () and in-plane reordering.
Before the sliding simulations begin, we prepared equilibrium starting configurations. In our protocol, the starting configurations were obtained by filling the open gap between initially distant plates, in any desired state of charge, with IL, and then reducing the gap width with = 0.11 m/s, generally slow enough to enact a reasonably adiabatic squeezout, as was found in Ref.capozza2015. The adiabaticity is facilitated by our chosen temperature , substantially higher the TM bulk melting temperature capozza2015, but as will be discussed below it is not automatically guaranteed, because both confinement and charging favor an increase of viscosity and a tendency to solidify. Above a minimal interplate distance and for small plate charge, the simulated IL is completely liquid, even if structured in close vicinity of the plates owing to their rigid wall nature and also to their charging. Interestingly, this near-plate local structure develops not only along , i.e., normal to the interface, but also along , parallel to the interface. Parallel IL ordering, rather more difficult to detect experimentally, is currently being pursued and detected with AFM tools segura13; elbourne15. As is further and gradually decreased, the confinement effects increase, the liquid-like diffusion diminishes and below a critical film thickness there is an effective interlocking of the structure emanating from the two facing plates, leading to a freezing of the IL into a solid-like layered arrangement with an odd number of alternating charge layers – when the two plates are neutral or identically charged. While no claim can obviously be made to represent a realistic situation, we note that IL confinement between mica plates, typically used in SFA friction experiments, is known to behave precisely like that perkin10; atkin07.
As was shown in the previous paper capozza2015 our choice of scheme and geometry permits additional thermodynamical insights. Treating the interplate separation as the external control variable, a free enthalpy can be calculated first estimating the free energy by integration of the average force measured between the plates as their distance is reduced.
where is the average z-oriented force exerted by the IL on the plates while they are a distance apart. After that the free enthalpy at given load and charge can be obtained as
Free enthalpy curves are displayed in Fig.2 for two values of normal load ,
and predict the load-induced transition layers.
The plate-plate free enthalpy generally displays
as a function of several odd- minima whose relative values depend on plate charging and load.
Once specified charge and load, one of
them is the absolute minimum, therefore the thermodynamically stable layered state, and all others are higher,
metastable states, which may still be very long lived. If the layered state is obtained by applying a force rather than by
forcing a distance , the precise number of layers attained in the course of a given simulation (or experiment) will be
history-dependent, in that more than one observed thickness may correspond to the same applied force.
This picture is quite close to what is observed in SFA and AFM perkin10; atkin2; spencer.
In addition to the alternating charge layer structure, we had previously found that the
confinement-induced solid structure of the TM model, although far from completely cystalline,
carried a second partly crystalline feature, consisting of vertical, neutral planes, inside
which the anions and cations arrange in a roughly square lattice, albeit with totally
disordered tail orientations capozza2015. Further decrease of caused this nearly solid IL
to squeeze out by successive pairs of ion planar layers, positive and negative, thus preserving local charge neutrality,
reducing the extension of vertical planes.
Starting with large and an initially fluid IL between the plates, the force resisting squeezout, initially zero, rises as drops. The IL is partly squeezed out, and the part that remains trapped acquires a structure consisting of layers (odd for equally charged plates, even for opposite charging). The squeezout transitions are clear first order transitions in this picture, the free enthalpy jumping across barriers between successive local minima. Decreasing , and thus increasing the lowest enthalpy state moves to a lower and lower number of layers. Each of these states with progressively lower, integer number of IL layers provides a starting state for our successive sliding simulations.
Upon checking, we actually found that configurations generated by fast squeezouts often possess internal layers that are only partly populated. This kind of ill-prepared states can survive as metastable for a long time. When faced with these ill-prepared states, initial particle configurations were further refined until a sufficiently reliable equilibrium particle population was reached. That allowed us to discard as a rule metastable initial configurations, and to describe the sliding behaviour of the stable, history independent ones, striving to discard metastable states as much as possible within our modest simulation durations.
Iii Simulated sliding
As mentioned earlier, there are two types of plate charging relevant for sliding friction. The first type, and the one generally realized experimentally, is single-plate charging relative to an electrode in electrical contact with the IL. In this case the aspect that matters is the IL structuring close to the charged plate, whereas the second plate, or tip, often insulating, acts as a merely mechanical probe of that structure. We replace this situation by simulating two equally charged plates, neutralized by adding the right amount of extra ions to the IL. The second type, whose experimental feasibility we ignore, but which looks even more interesting, is that of opposite plate charging. This we simulate without difficulty, restricting the charge magnitude to values small enough to correspond to an interplate voltage below 5 V.
We start off our sliding simulations with equally charged plates, endowed with a negative charge density (about 1/3 of the nominal value quoted for mica in experimental conditions) and with the IL confined between them in a variety of configurations, obtained as a result of squeezout down to different values, as detailed in the previous section, further equilibrated under an applied load , chosen to be close to . Sliding friction depends upon the IL layer number and the load .
Fig.3a illustrates two sliding simulations starting from =7 with either nN and nN. The simulated system at nN maintains layers confined in between the plates, while that at nN undergoes a relayering transition, whose kinetics seems favored by sliding.
As shown in Fig.3b, this transition corresponds to the sudden drop of the gap , associated with
the expulsion of the pair of layers, replacing the =7 metastable state with the stable =5.
This is an example of how the lateral driving can facilitate the squeezout of layers by mechanically perturbing the IL.
The close proximity of the transition can be predicted by free enthalpy curves, displayed in Fig.3c. At nN the energy barrier separating the states with and layers drops. The role of sliding appears to be kinetic, facilitating the jump to the more stable layer state.
Iv - dependent frictional shear stress
The calculated force-distance curves of Ref. capozza2015 show that a single normal load generally gives rise to several long lived states with a different number . For example nN can support capozza2015, very much as it happens in experiment hoth14; perkin; perkin10. Figure 4a shows the lateral force obtained as a function of time for different values, with nN and . Static free enthalpy curve in Fig.4b demonstrates that = 7 is stable configuration while 3 and 5 are metastable. A strong friction increase is observed with decreasing , in agreement with expectations and with SFA data perkin; bou-malham10.
The detailed load dependence of the frictional shear stress obtained for each fixed number of layers is more interesting, and is illustrated in Fig.5a.
Here the growth of shear stress with increasing load is surprisingly modest. While that might appear at first sight in disagreement with experiment,
which show a nearly linear increase of kinetic friction with load perkin, we must recall that the contact area is not a constant in SFA,
whereas shear stresses require normalizing to the area . Assuming for SFA two perpendicularly oriented
elastic cylinders in Hertzian contact, the area , where is the effective force between the two
sliders. Extracting an approximate adhesion force by linearly extrapolating the experimental kinetic friction
to zero in Fig.1 of Ref. perkin , we recover the ”experimentally derived” shear stress vs load curves in
Fig 5b. Remarkably, they now show a weaker experimental dependence on load,
actually very close to that predicted by simulations in Fig. 5a.
To understand the strong dependence of on the number of layers illustrated in Fig.4a, we now take a closer look to the lubricant shear velocity distribution during sliding. Fig.6 shows the layer-by-layer resolved ion density and velocity profiles along in the liquid for 7, 5 and 3 layers, nN and m/s in correspondence of green, blue and red symbols in Fig.5a respectively.
For the velocity profile deviates from the linear behavior (red dashed line) indicating that the layers at the boundaries are bound to the plates and tend to move at their same velocity. The plate charging promotes the IL in immediate contact to nearly solid or at least to a much higher viscosity. Moreover, for this relatively high level of charging, the density ratio of first-layer cations over the anode, and of anions over the cathode, is about 1.4, that is only slightly larger than one. kornyshev Both elements support a good grip by the plates. The three inner layers constitute by contrast the more liquid, effective shear band where the sliding is concentrated, and where the coalescence between the two solid-like structures protruding from opposite plates is easiest to fracture. By reducing the number of layers to 5 and then to 3, the shear band moves closer and closer to the boundaries and even the more solid-like structures closer to the plates must break up and move, as shown in Fig. 6d, f. As a consequence the lateral force increases. A pictorial representation of this behavior is provided by the cartoons of Fig.6. Surprising at first sight is the seemingly ”viscous” friction for , with a linear speed profile distribution between the three layers, the absence of strong stick-slip as seen in in Fig.4a and, as we shall see later in Fig.9b, a linear friction increase with speed. That behavior coexists with a very large shear stress, and a very solid structure of this extremely thin lubricant film. The evolution of layer-by-layer relative velocities shown in Fig.6 indicates that the yield stresses of the plate/IL interface and of the mid-film intra-IL interface, which for =7,5 were quite different, have now become comparable for =3. As a result all interfaces slide by about the same amount. The meandering wall static structure adopted by the IL in this pressured-and-charged regime capozza2015 now incessantly ruptures and rearranges during sliding, leading to a z-uniform velocity gradient, reminscent of fluid lubrication, despite its nature very different from that of a fluid.
The picture changes drastically by reducing the charge to symmetrically on both plates. The mid-film shear band disappears and the slippage occurs mostly at the interfaces between the IL and the plates, as shown in Fig.7. While for the top plate moved by rupturing the solid-like structures reaching across plates, at this lower charge the film behaves effectively as a solid block confined between two slippery surfaces. The shearing concentrates mostly at the two interfaces, whence the shear stress depends much more on the plate-IL interaction.
The frictional shear stress is about an order of magnitude smaller than for higher charge and displays a much weaker dependence on , as seen by comparing Fig.8b with Fig.5a. The rise of average frictional force for , displayed in Fig.8a (red line) is associated with a stick-slip regime, as shown by instantaneous spring forces in Fig.8a corresponding to the points indicated in Fig.8b. This is in contrast to the case of high charge, where friction was large but the sliding was smooth.
The essence for this difference of behavior is explained by the presence, at low but not at large plate charging,
of overscreening, a phenomenon well demonstrated in earlier work kornyshev. At the small
charging of we find that the charge density ratio of first-layer cations over
that of the plate anode has in fact risen to about 2.3, much larger than one. The overscreened plate
is covered by a dense counterion blanket, making the plate-IL adhesive grip quite effective. The IL film,
on the other hand, is solid and even harder to fracture. The plate-IL interface yields under sliding and that
V Velocity dependence of friction
After the charge, load, and layer number dependence of confined IL friction, we close the symmetric plate charging part of this study by investigating the velocity dependence. A starting consideration here is that experimental sliding velocities in SFA and AFM are very small, typically reaching , a value which is at least 6 orders of magnitude lower than velocities realistically accessible in atomic and molecular level simulations. This is a standard difficulty, well known and amply discussed vanossi13; one can in fact learn enough even by simulating at very high speeds. Essentially, friction can be either smooth and viscous, in which case the shear stress is small and proportional to velocity, or it can be stick-slip, in which case the shear stress is large and very poorly dependent on velocity. In both cases, very crudely speaking, even a fast simulation can convey the necessary information (the viscous friction coefficient in the first, the full shear stress in the second). Of course it will always be necessary to bear in mind that a) increasing velocity may turn stick-slip into smooth sliding; and b) microscopic scale stick-slip does not necessarily show up as macroscopic or mesoscopic oscillations, and is generally detectable as a velocity independent shear stress; c) quantitative aspects are controlled by parameters such as the slider’s masses, the pulling spring constant, etc., whose assumedvalues are largely arbitrary.
With these provisions we investigate the dependence of IL friction upon driving velocity ,
at fixed normal load . Fig.9 and 10 show the frictional shear
stress as a function of for equal-sign charge densities of
and respectively and for .
All results for in Fig.9 show a linear increase of shear
stress with velocity indicating a viscous-type friction. As suggested earlier, this is most
likely attributable to overscreening of both plates by the extra anions. However
vs does not extrapolate to , so that even at a finite force
is needed to set the system in motion. This is the force needed in order to rupture
the vertical IL ”solid” walls that straddle across the two charged plates (Fig.5).
Also for low charge density , the shear stress for grows proportional to velocity but here it extrapolates to zero as (see Fig.10b). This is the hallmark of true viscous friction, and it agrees with the top plate moving in a clear regime of smooth sliding. The sliding regime drastically changes for and low velocity , where a clear stick-slip regime takes over. The alternation of sticking periods followed by rapid and highly dissipative motion of the top plate is responsible for the strong increase of average lateral force at low velocities, as demonstrated in Fig.10a. By increasing , the stick-slip first turns chaotic and eventually disappears at , where smooth sliding is recovered. Particle trajectories clearly show this transition. At low the IL structure retains its worm-like meandering planes or chains, which we described in ref. capozza2015, and which were recently observed in experiments, segura13 and behaves essentially as a solid confined between two slippery planes.
Vi Friction with opposite plate charging
In this final part we study how friction force is affected by opposite charging of the two confining plates, such as one would obtain by applying a DC voltage between them. Probably not that easy to realize experimentally, this configuration is nonetheless quite interesting. Upon balancing the opposite charge on the two plates, the equilibrium layering structure now forms an even number of layers, a rearrangement which preserves local charge neutrality. In order to explore expeditely the dependence of friction upon the magnitude of the charge imbalance, we conduct a time-dependent simulation where we periodically, but very slowly and as adiabatically as possible, modulate in time the on-plate charge in the form , where the two opposite signs refer to the two opposite plates. Adiabaticity is pursued by choosing a period large enough compared to the characteristic time scale of all relevant ionic rearrangements.
The magnitude of charging used in this work, is dictated, as discussed previously capozza2015, by the requirement that, assuming an average dielectric constant of 2, the plate-plate voltage should not exceed 5V, taken as a reasonable upper limit in an experimental situation. In this respect therefore, our study does not overlap with previous work, where in a different model fajardo15 the behavior at was investigated.
Assuming a constant load nN, a charge oscillation period ns, a magnitude and a sliding velocity we obtain the spontaneous evolution of the system shown in Fig.11. What happens is, in a nutshell, that the interplate distance spontaneously changes, and so does the sliding friction as the charge changes.
The simulation begins with a charge density of on plates and the IL arranges in one single neutral layer indicated as in Fig.11, which, as shown in the static squeezout enthalpy of Fig.12, is the most stable configuration under the fixed applied loading force nN and low charge.
As the plate charge increases, the liquid in the lateral reservoirs is sucked in
the gap (despite a constant load), reaching a four layer solid state.
Actually, we expected it to go back to six layers, which is the
equilibrium free enthalpy minimum, as shown in Fig.12. It probably would
if we could simulate for much longer times, but it does not on our limited time scale,
due to tha strongly reduced mobility of ions in the near-solid film.
This strong and rigid four layers structure softens as the plate charge decreases,
eventually turning more fluid and squeezing out reaching again the single neutral layer
arrangement . These charging-induced relayering and squeezout are accompanied by
drastic change in frictional behaviour of the liquid film. The large shear stress at
high plate charge (high voltage, near 5 eV) is replaced by a nearly vanishing shear
stress at zero charge. Here therefore, contrary to the equal-charge case of the previous
sections, the thinnest film yields the lowest friction, rather than the opposite.
This is because, independently of sliding, the plate charging induces structural
transformations of freezing/melting of the the IL filmcapozza2015.
The frozen IL supports the applied load for large charge but the liquid
IL does not for zero charge. The IL
film structure at shrinks to a single, planar, square and neutral IL monolayer.
Here the neutral tails stick out of the plane, effectively lubricating the sliding.
A qualitatively similar mechanism where alkyl tails stick out of a well-packed anion-cation
plane has been proposed in a recent experimental study bennewitz
to explain the potential-dependent frictional behaviour of an AFM tip in presence of an IL.
We did change and explore a variety of parameters, including charging magnitude, applied load, and sliding speed in order to test the main aspects of these results. The processes of melting and freezing, and of expulsion and sucking always remained very asymmetrical between the fast squeezout, associated with sudden melting of the IL film upon decreasing towards zero, and the more gradual sucking-in and freezing, attained upon increasing from zero in Fig.11a. The electrically driven fast emptying but slow refilling of the interplate gap could be used, together with some spatial asymmetry to build an effective effective electro-pumping of ions in and out of the gap.
The frictional shear stress shows a very large charge dependence, as a direct consequence of the changes of structure and of gap width caused by charging. Although noisy, the large shear stress at charge maxima drops by an order of magnitude or more following squeezout at zero charge. This strong effect appears promising in view of experimental verifications, and of a possible use for the control of friction.
Lastly, we briefly dwell on questions connected with work, heat and temperature, which come about with respect both to charging-induced melting and freezing, and to frictional work and heating. In the absence of sliding, electrically induced freezing and melting would be associated with the absorption or release of latent heat. That would be large if the solid was ordered, e.g., crystalline or nearly crystalline, as is the case for example in the NaCl model of the IL capozza2015. The more the solid resembles a glass, however, the smaller the latent heat will be, owing to structural near-identity between glass and liquid. That is the case of our TM model, where indeed the latent heat of melting was very small in the simulated caloric curves capozza2015. Without sliding, the presence of a latent heat is undetectable in electrically induced freezing and melting cycles shown in Fig.13c,d. Of course, these data were obtained with a thermostat, canceling any large temperature fluctuations: despite the themostat, dips and peaks would be well visible in presence of a robust, crystal-like latent heat. Frictional heating is instead quite visible in the sliding results under cycled charging of Fig.13a,b. Despite the thermostat, the temperature is seen to oscillate, downwards when the solid confined IL lubricant goes liquid at , up when it goes back to solid at and where therefore the frictional Joule heating is much stronger, down again upon melting where the friction drops once again. Because of the presence of the thermostat we did not attempt a quantitative connection between the frictional work and the residual temperature oscillation, but the connection between the two is quite clear.
The electrical charging of plates in nearly atomic contact influences an IL when trapped in a nanoscale sized gap between them. As a result, the frictional shear stress realized upon mutual sliding of the plates is modified by the charging. The present simulation study conducted within the simple ”tailed model” previously developed for the IL explores some of the rich variety of phenomena that can take place in this context. The IL becomes structured into alternating charge layers, whose number is odd when the plates are equally charged, even when they are oppositely charged. The trapped IL film may develop a solid-like rigidity, with different and opposite frictional consequences depending on charging level and other parameters. When the plate grips the solid-like film and the shearing takes place in its middle the friction is large; but when the grip is less effective the shear concentrates at the plate-film interface and friction drops. Alternatively, at zero charging and under the same conditions the ionic liquid can melt, thus lubricating and dramatically reducing the sliding friction.
Even though the parameter-dependent nature of our results does not make a comparison with any specific experimental system particularly compelling, we believe that several elements uncovered either connect with observations already available, or else suggest pursuing newer ones. The modest dependence of friction upon load – unusual in the context of classical friction – is among the former. The possibility that electrical charging, in whichever manners that could be realized, could influence the tribology of confined IL by causing its electrically driven freezing and melting is among the latter, and deserves to be tested, supporting the possibility of an electric control of IL lubricated friction.
The authors are grateful to Daniele Passerone and Carlo Pignedoli from EMPA (Dübendorf, CH) for the computational resources and the technical assistance provided, and, at various times, for discussions with A. Kornyshev and A. Schirmeisen. Work in Trieste was partly sponsored by Sinergia Contract CRSII236287/1, and by ERC Advanced Grant 320796 - MODPHYSFRICT. This work is also supported by the COST Action MP1303 ”Understanding and Controlling Nano and Mesoscale Friction”.