Graphene nanoribbons grown in hBN stacks for high-performance electronics.

Van der Waals encapsulation of two-dimensional materials in hexagonal boron nitride (hBN) stacks is a promising way to create ultrahigh-performance electronic devices1-4. However, contemporary approaches for achieving van der Waals encapsulation, which involve artificial layer stacking using mechanical transfer techniques, are difficult to control, prone to contamination and unscalable. Here we report the transfer-free direct growth of high-quality graphene nanoribbons (GNRs) in hBN stacks. The as-grown embedded GNRs exhibit highly desirable features being ultralong (up to 0.25 mm), ultranarrow (<5 nm) and homochiral with zigzag edges. Our atomistic simulations show that the mechanism underlying the embedded growth involves ultralow GNR friction when sliding between AA'-stacked hBN layers. Using the grown structures, we demonstrate the transfer-free fabrication of embedded GNR field-effect devices that exhibit excellent performance at room temperature with mobilities of up to 4,600 cm2 V-1 s-1 and on-off ratios of up to 106. This paves the way for the bottom-up fabrication of high-performance electronic devices based on embedded layered materials.

Cumulative polarization in conductive interfacial ferroelectrics.

Ferroelectricity in atomically thin bilayer structures has been recently predicted1 and measured2-4 in two-dimensional materials with hexagonal non-centrosymmetric unit-cells. The crystal symmetry translates lateral shifts between parallel two-dimensional layers to sign changes in their out-of-plane electric polarization, a mechanism termed 'slide-tronics'4. These observations have been restricted to switching between only two polarization states under low charge carrier densities5-12, limiting the practical application of the revealed phenomena13. To overcome these issues, one should explore the nature of polarization in multi-layered van der Waals stacks, how it is governed by intra- and interlayer charge redistribution and to what extent it survives the addition of mobile charge carriers14. To explore these questions, we conduct surface potential measurements of parallel WSe2 and MoS2 multi-layers with aligned and anti-aligned configurations of the polar interfaces. We find evenly spaced, nearly decoupled potential steps, indicating highly confined interfacial electric fields that provide a means to design multi-state 'ladder-ferroelectrics'. Furthermore, we find that the internal polarization remains notable on electrostatic doping of mobile charge carrier densities as high as 1013 cm-2, with substantial in-plane conductivity. Using density functional theory calculations, we trace the extra charge redistribution in real and momentum spaces and identify an eventual doping-induced depolarization mechanism.

Theory and Simulations of Ionic Liquids in Nanoconfinement.

Room-temperature ionic liquids (RTILs) have exciting properties such as nonvolatility, large electrochemical windows, and remarkable variety, drawing much interest in energy storage, gating, electrocatalysis, tunable lubrication, and other applications. Confined RTILs appear in various situations, for instance, in pores of nanostructured electrodes of supercapacitors and batteries, as such electrodes increase the contact area with RTILs and enhance the total capacitance and stored energy, between crossed cylinders in surface force balance experiments, between a tip and a sample in atomic force microscopy, and between sliding surfaces in tribology experiments, where RTILs act as lubricants. The properties and functioning of RTILs in confinement, especially nanoconfinement, result in fascinating structural and dynamic phenomena, including layering, overscreening and crowding, nanoscale capillary freezing, quantized and electrotunable friction, and superionic state. This review offers a comprehensive analysis of the fundamental physical phenomena controlling the properties of such systems and the current state-of-the-art theoretical and simulation approaches developed for their description. We discuss these approaches sequentially by increasing atomistic complexity, paying particular attention to new physical phenomena emerging in nanoscale confinement. This review covers theoretical models, most of which are based on mapping the problems on pertinent statistical mechanics models with exact analytical solutions, allowing systematic analysis and new physical insights to develop more easily. We also describe a classical density functional theory, which offers a reliable and computationally inexpensive tool to account for some microscopic details and correlations that simplified models often fail to consider. Molecular simulations play a vital role in studying confined ionic liquids, enabling deep microscopic insights otherwise unavailable to researchers. We describe the basics of various simulation approaches and discuss their challenges and applicability to specific problems, focusing on RTIL structure in cylindrical and slit confinement and how it relates to friction and capacitive and dynamic properties of confined ions.

Structural superlubricity and ultralow friction across the length scales.

Structural superlubricity, a state of ultralow friction and wear between crystalline surfaces, is a fundamental phenomenon in modern tribology that defines a new approach to lubrication. Early measurements involved nanometre-scale contacts between layered materials, but recent experimental advances have extended its applicability to the micrometre scale. This is an important step towards practical utilization of structural superlubricity in future technological applications, such as durable nano- and micro-electromechanical devices, hard drives, mobile frictionless connectors, and mechanical bearings operating under extreme conditions. Here we provide an overview of the field, including its birth and main achievements, the current state of the art and the challenges to fulfilling its potential.

Orientational Ordering in Nano-confined Polar Liquids.

Water and other polar liquids exhibit nanoscale structuring near charged interfaces. When a polar liquid is confined between two charged surfaces, the interfacial solvent layers begin to overlap, resulting in solvation forces. Here, we perform molecular dynamics simulations of polar liquids with different dielectric constants and molecular shapes and sizes confined between charged surfaces, demonstrating strong orientational ordering in the nanoconfined liquids. To rationalize the observed structures, we apply a coarse-grained continuum theory that captures the orientational ordering and solvation forces of those liquids. Our findings reveal the subtle behavior of different nanoconfined polar liquids and establish a simple law for the decay length of the interfacial orientations of the solvents, which depends on their molecular size and polarity. These insights shed light on the nature of solvation forces, which are important in colloid and membrane science, scanning probe microscopy, and nano-electrochemistry.

Electrotunable friction with ionic liquid lubricants.

Room-temperature ionic liquids and their mixtures with organic solvents as lubricants open a route to control lubricity at the nanoscale via electrical polarization of the sliding surfaces. Electronanotribology is an emerging field that has a potential to realize in situ control of friction-that is, turning the friction on and off on demand. However, fulfilling its promise needs more research. Here we provide an overview of this emerging research area, from its birth to the current state, reviewing the main achievements in non-equilibrium molecular dynamics simulations and experiments using atomic force microscopes and surface force apparatus. We also present a discussion of the challenges that need to be solved for future applications of electrotunable friction.

Ionotronics for reverse actuation.

In the midst of an ongoing energy crisis, the search for new methods of energy harvesting has never been more important. Here we explore, analyse and discuss principles of ionotronic reverse-actuator devices based on the effect of double-layer charging. The designs that we consider in this paper operate based on a common principle - using external mechanical work, which would otherwise be wasted, to produce changes in the contact area of electrode and electrolyte, translated into the time variation of the double-layer capacitance. Periodic variation of capacitance, when connected to a reference voltage source, produces alternating electric current through a load. This concept is not new and in some forms was realised in the early works of Boland, Krupenkin and several papers of our group. The goal of the present paper is to build a comprehensive analytical platform for a description of operation of such devices in terms of materials, generated power as a function of the frequency of variation of applied force, electrical load, and other factors; the understanding of which allows us to optimise these systems and navigate their construction. The first design, discussed in the paper, is based on flat electrodes. It is the simplest one and, as such, helps elucidate some key factors determining power generation. While being easy to realise experimentally, it generates relatively low power, even when optimised. The second design, based on microporous electrodes is more sophisticated and allows a much larger power harvest. The results are also compared to the recently proposed capacitive rotor device. The developed theory is set to capture the key factors that determine the functioning of the considered reverse actuators. The structures under study are matched to fit into the sole of a shoe and produce power from walking and running. However, they can also be scaled-up to larger operating systems and various external loads.

Electric-field frictional effects in confined zwitterionic molecules.

We theoretically explore the effect of a transverse electric field on the frictional response of a bi-layer of packed zwitterionic molecules. The dipole-moment reorientation promoted by the electric field can lead to either stick-slip or smooth sliding dynamics, with average shear stress values varying over a wide range. A structure-property relation is revealed by investigating the array of molecules and their mutual orientation and interlocking. Moreover, the thermal friction enhancement previously observed in these molecules is shown to be suppressed by the electric field, recovering the expected thermolubricity at large-enough fields. The same holds for other basic tribological quantities, such as the external load, which can influence friction in opposite ways depending on the strength of the applied electric field. Our findings open a route for the reversible control of friction forces via electric polarization of the sliding surface.

Anisotropic Interlayer Force Field for Group-VI Transition Metal Dichalcogenides.

An anisotropic interlayer force field that describes the interlayer interactions in homogeneous and heterogeneous interfaces of group-VI transition metal dichalcogenides (MX2, where M = Mo, W, and X = S, Se) is presented. The force field is benchmarked against density functional theory calculations for bilayer systems within the Heyd-Scuseria-Ernzerhof hybrid density functional approximation, augmented by a nonlocal many-body dispersion treatment of long-range correlation. The parametrization yields good agreement with the reference calculations of binding energy curves and sliding potential energy surfaces. It is found to be transferable to transition metal dichalcogenide (TMD) junctions outside of the training set that contain the same atom types. Calculated bulk moduli agree with most previous dispersion-corrected density functional theory predictions, which underestimate the available experimental values. Calculated phonon spectra of the various junctions under consideration demonstrate the importance of appropriately treating the anisotropic nature of the layered interfaces. Considering our previous parametrization for MoS2, the anisotropic interlayer potential enables accurate and efficient large-scale simulations of the dynamical, tribological, and thermal transport properties of a large set of homogeneous and heterogeneous TMD interfaces.

Load-induced dynamical transitions at graphene interfaces.

The structural superlubricity (SSL), a state of near-zero friction between two contacted solid surfaces, has been attracting rapidly increasing research interest since it was realized in microscale graphite in 2012. An obvious question concerns the implications of SSL for micro- and nanoscale devices such as actuators. The simplest actuators are based on the application of a normal load; here we show that this leads to remarkable dynamical phenomena in microscale graphite mesas. Under an increasing normal load, we observe mechanical instabilities leading to dynamical states, the first where the loaded mesa suddenly ejects a thin flake and the second characterized by peculiar oscillations, during which a flake repeatedly pops out of the mesa and retracts back. The measured ejection speeds are extraordinarily high (maximum of 294 m/s), and correspond to ultrahigh accelerations (maximum of 1.1×1010 m/s2). These observations are rationalized using a simple model, which takes into account SSL of graphite contacts and sample microstructure and considers a competition between the elastic and interfacial energies that defines the dynamical phase diagram of the system. Analyzing the observed flake ejection and oscillations, we conclude that our system exhibits a high speed in SSL, a low friction coefficient of 3.6×10-6, and a high quality factor of 1.3×107 compared with what has been reported in literature. Our experimental discoveries and theoretical findings suggest a route for development of SSL-based devices such as high-frequency oscillators with ultrahigh quality factors and optomechanical switches, where retractable or oscillating mirrors are required.

Velocity Dependence of Moiré Friction.

Friction force microscopy experiments on moiré superstructures of graphene-coated platinum surfaces demonstrate that in addition to atomic stick-slip dynamics, a new dominant energy dissipation route emerges. The underlying mechanism, revealed by atomistic molecular dynamics simulations, is related to moiré ridge elastic deformations and subsequent relaxation due to the action of the pushing tip. The measured frictional velocity dependence displays two distinct regimes: (i) at low velocities, the friction force is small and nearly constant; and (ii) above some threshold, friction increases logarithmically with velocity. The threshold velocity, separating the two frictional regimes, decreases with increasing normal load and moiré superstructure period. Based on the measurements and simulation results, a phenomenological model is derived, allowing us to calculate friction under a wide range of room temperature experimental conditions (sliding velocities of 1-104 nm/s and a broad range of normal loads) and providing excellent agreement with experimental observations.

Stick-Slip Dynamics of Moiré Superstructures in Polycrystalline 2D Material Interfaces.

A new frictional mechanism, based on collective stick-slip motion of moiré superstructures across polycrystalline two-dimensional material interfaces, is predicted. The dissipative stick-slip behavior originates from an energetic bistability between low- and high-commensurability configurations of large-scale moiré superstructures. When the grain boundary separates between grains of small and large interfacial twist angle, the corresponding moiré periods are significantly different, resulting in forbidden grain boundary crossing of the moiré superstructures during shear induced motion. For small twist angle grains, where the moiré periods are much larger than the lattice constant, this results in multiple reflections of collective surface waves between the surrounding grain boundaries. In combination with the individual grain boundary dislocation snap-through buckling mechanism dominating at the low normal load regime, the friction exhibits nonmonotonic behavior with the normal load. While the discovered phenomenon is demonstrated for h-BN/graphene polycrystalline junctions, it is expected to be of general nature and occur in many other large-scale layered material interfaces.

Parity-Dependent Moiré Superlattices in Graphene/h-BN Heterostructures: A Route to Mechanomutable Metamaterials.

The superlattice of alternating graphene/h-BN few-layered heterostructures is found to exhibit strong dependence on the parity of the number of layers within the stack. Odd-parity systems show a unique flamingolike pattern, whereas their even-parity counterparts exhibit regular hexagonal or rectangular superlattices. When the alternating stack consists of 7 layers or more, the flamingo pattern becomes favorable, regardless of parity. Notably, the out-of-plane corrugation of the system strongly depends on the shape of the superstructure resulting in significant parity dependence of its mechanical properties. The predicted phenomenon originates in an intricate competition between moiré patterns developing at the interface of consecutive layers. This mechanism is of general nature and is expected to occur in other alternating stacks of closely matched rigid layered materials as demonstrated for homogeneous alternating junctions of twisted graphene and h-BN. Our findings thus allow for the rational design of mechanomutable metamaterials based on van der Waals heterostructures.

Structural effects in nanotribology of nanoscale films of ionic liquids confined between metallic surfaces.

Room Temperature Ionic Liquids (RTILs) attract significant interest in nanotribology. However, their microscopic lubrication mechanism is still under debate. Here, using non-equilibrium molecular dynamics simulations, we investigate the lubrication performance of ultra-thin (<2 nm) films of [C2MIM]+ [NTf2]- confined between plane-parallel neutral surfaces of Au(111) or Au(100). We find that films consisting of tri-layers or bilayers, form ordered structures with a flat orientation of the imidazolium rings with respect to the gold surface plane. Tri-layers are unstable against loads >0.5 GPa, while bi-layers sustain pressures in the 1-2 GPa range. The compression of these films results in monolayers that can sustain loads of several GPa without significant loss in their lubrication performance. Surprisingly, in such ultra-thin films the imidazolium rings show higher orientational in-plane disorder, with and the rings adopting a tilted orientation with respect to the gold surface. The friction force and friction coefficient of the monolayers depends strongly on the structure of the gold plates, with the friction coefficient being four times higher for monolayers confined between Au(100) surfaces than for more compact Au(111) surfaces. We show that the general behaviour described here is independent of whether the metallic surfaces are modelled as polarizable or non-polarizable surfaces and speculate on the nature of this unexpected conclusion.

Controllable Thermal Conductivity in Twisted Homogeneous Interfaces of Graphene and Hexagonal Boron Nitride.

Thermal conductivity of homogeneous twisted stacks of graphite is found to strongly depend on the misfit angle. The underlying mechanism relies on the angle dependence of phonon-phonon couplings across the twisted interface. Excellent agreement between the calculated thermal conductivity of narrow graphitic stacks and corresponding experimental results indicates the validity of the predictions. This is attributed to the accuracy of interlayer interaction descriptions obtained by the dedicated registry-dependent interlayer potential used. Similar results for h-BN stacks indicate overall higher conductivity and reduced misfit angle variation. This opens the way for the design of tunable heterogeneous junctions with controllable heat-transport properties ranging from substrate-isolation to efficient heat evacuation.

Model for Bundling of Keratin Intermediate Filaments.

Keratin intermediate filaments form dynamic intracellular networks, which span the entire cytoplasm and provide mechanical strength to the cell. The mechanical resilience of the keratin intermediate filament network itself is determined by filament bundling. The bundling process can be reproduced in artificial conditions in the absence of any specific cross-linking proteins, which suggests that it is driven by generic physical forces acting between filaments. Here, we suggest a detailed model for bundling of keratin intermediate filaments based on interfilament electrostatic and hydrophobic interactions. It predicts that the process is limited by an optimal bundle thickness, which is determined by the electric charge of the filaments, the number of hydrophobic residues in the constituent keratin polypeptides, and the extent to which the electrolyte ions are excluded from the bundle interior. We evaluate the kinetics of the bundling process by considering the energy barrier a filament has to overcome for joining a bundle.

Origin of Friction in Superlubric Graphite Contacts.

More than thirty years ago, it was theoretically predicted that friction for incommensurate contacts between atomically smooth, infinite, crystalline materials (e.g., graphite, MoS_{2}) is vanishing in the low speed limit, and this corresponding state was called structural superlubricity (SSL). However, experimental validation of this prediction has met challenges, since real contacts always have a finite size, and the overall friction arises not only from the atoms located within the contact area, but also from those at the contact edges which can contribute a finite amount of friction even when the incommensurate area does not. Here, we report, using a novel method, the decoupling of these contributions for the first time. The results obtained from nanoscale to microscale incommensurate contacts of graphite under ambient conditions verify that the average frictional contribution of an inner atom is no more than 10^{-4} that of an atom at the edge. Correspondingly, the total friction force is dominated by friction between the contact edges for contacts up to 10 µm in lateral size. We discuss the physical mechanisms of friction observed in SSL contacts, and provide guidelines for the rational design of large-scale SSL contacts.

Robust microscale superlubricity in graphite/hexagonal boron nitride layered heterojunctions.

Structural superlubricity is a fascinating tribological phenomenon, in which the lateral interactions between two incommensurate contacting surfaces are effectively cancelled resulting in ultralow sliding friction. Here we report the experimental realization of robust superlubricity in microscale monocrystalline heterojunctions, which constitutes an important step towards the macroscopic scale-up of superlubricity. The results for interfaces between graphite and hexagonal boron nitride clearly demonstrate that structural superlubricity persists even when the aligned contact sustains external loads under ambient conditions. The observed frictional anisotropy in the heterojunctions is found to be orders of magnitude smaller than that measured for their homogeneous counterparts. Atomistic simulations reveal that the underlying frictional mechanisms in the two cases originate from completely different dynamical regimes. Our results are expected to be of a general nature and should be applicable to other van der Waals heterostructures.

Negative Friction Coefficients in Superlubric Graphite-Hexagonal Boron Nitride Heterojunctions.

Negative friction coefficients, where friction is reduced upon increasing normal load, are predicted for superlubric graphite-hexagonal boron nitride heterojunctions. The origin of this counterintuitive behavior lies in the load-induced suppression of the moiré superstructure out-of-plane distortions leading to a less dissipative interfacial dynamics. Thermally induced enhancement of the out-of-plane fluctuations leads to an unusual increase of friction with temperature. The highlighted frictional mechanism is of a general nature and is expected to appear in many layered material heterojunctions.

Load and Velocity Dependence of Friction Mediated by Dynamics of Interfacial Contacts.

Studying the frictional properties of interfaces with dynamic chemical bonds advances understanding of the mechanism underlying rate and state laws, and offers new pathways for the rational control of frictional response. In this work, we revisit the load dependence of interfacial chemical-bond-induced (ICBI) friction experimentally and find that the velocity dependence of friction can be reversed by changing the normal load. We propose a theoretical model, whose analytical solution allows us to interpret the experimental data on timescales and length scales that are relevant to experimental conditions. Our work provides a promising avenue for exploring the dynamics of ICBI friction.