Pervasive orientational and directional locking at geometrically heterogeneous sliding interfaces

Evidence of Directional Structural Superlubricity and Lévy Flights in a van der Waals Heterostructure

Structural superlubricity is a special frictionless contact in which two crystals are in incommensurate arrangement such that relative in-plane translation is associated with vanishing energy barrier crossing. So far, it has been realized in multilayer graphene and other van der Waals (2D crystals with hexagonal or triangular crystalline symmetries, leading to isotropic frictionless contacts. Directional structural superlubricity, to date unrealized in 2D systems, is possible when the reciprocal lattices of the two crystals coincide in one direction only. Here, directional structural superlubricity a α-bismuthene/graphite van der Waals system is evidenced, manifested by spontaneous hopping of the islands over hundreds of nanometers at room temperature, resolved by low-energy electron microscopy and supported by registry simulations. Statistical analysis of individual and collective α-bismuthene islands populations reveal a heavy-tailed distribution of the hopping lengths and sticking times indicative of Lévy flight dynamics, largely unobserved in condensed-matter systems.

Tuning the stability of a model quasicrystal and its approximants with a periodic substrate

Quasicrystals and their periodic approximants are complex crystalline phases. They have now been observed in many metallic alloys, soft matter systems, and particle simulations. In recent experiments of thin-film perovskites on solid substrates, the type of complex phase was found to change depending on thermodynamic conditions and the type of substrate used. Here, we investigate the effect of a substrate on the relative thermodynamic stability of a two-dimensional model quasicrystal and its approximants. Our simulation model is particles interacting via the Lennard-Jones-Gauss potential. Our numerical methods are molecular dynamics simulations and free energy calculations that take into account phason flips explicitly. For substrates interacting weakly with the particles, we observe an incommensurate-commensurate transition, in which a continuous series of quasicrystal approximants locks into a small number of approximants. Interestingly, we observe that the 3/2 approximant exhibits phason mode fluctuations in thermodynamic equilibrium. Such fluctuations are reminiscent of random tiling and a phenomenon usually associated only with quasiperiodic order. For stronger substrates, we find an enhancement of the stability of the dodecagonal quasicrystal and variants of square lattices. We explain all observed phenomena by the interplay of the model system with the substrate. Our results demonstrate that designing novel complex periodic and quasiperiodic structures by choice of suitable substrates is a promising strategy.

Twisting Dynamics of Large Lattice-Mismatch van der Waals Heterostructures

van der Waals (vdW) homo/heterostructures are ideal systems for studying interfacial tribological properties such as structural superlubricity. Previous studies concentrated on the mechanism of translational motion in vdW interfaces. However, detailed mechanisms and general properties of the rotational motion are barely explored. Here, we combine experiments and simulations to reveal the twisting dynamics of the MoS₂/graphite heterostructure. Unlike the translational friction falling into the superlubricity regime with no twist angle dependence, the dynamic rotational resistances highly depend on twist angles. Our results show that the periodic rotational resistance force originates from structural potential energy changes during the twisting. The structural potential energy of MoS₂/graphite heterostructure increases monotonically from 0° to 30° twist angles, and the estimated relative energy barrier is (1.43 ± 0.36) × 10^-3 J/m². The formation of Moiré superstructures in the graphene layer is the key to controlling the structural potential energy of the MoS₂/graphene heterostructure. Our results suggest that in twisting 2D heterostructures, even if the interface sliding friction is negligible, the evolving potential energy change results in a nonvanishing rotational resistance force. The structural change of the heterostructure can be an additional pathway for energy dissipation in the rotational motion, further enhancing the rotational friction force.

Clogging, dynamics, and reentrant fluid for active matter on periodic substrates

We examine the collective states of run-and-tumble active matter disks driven over a periodic obstacle array. When the drive is applied along a symmetry direction of the array, we find a clog-free uniform liquid state for low activity, while at higher activity, the density becomes increasingly heterogeneous and an active clogged state emerges in which the mobility is strongly reduced. For driving along nonsymmetry or incommensurate directions, there are two different clogging behaviors consisting of a drive-dependent clogged state in the low activity thermal limit and a drive-independent clogged state at high activity. These regimes are separated by a uniform flowing liquid at intermediate activity. There is a critical activity level above which the thermal clogged state does not occur, as well as an optimal activity level that maximizes the disk mobility. Thermal clogged states are dependent on the driving direction while active clogged states are not. In the low activity regime, diluting the obstacles produces a monotonic increase in the mobility; however, for large activities, the mobility is more robust against obstacle dilution. We also examine the velocity-force curves for driving along nonsymmetry directions and find that they are linear when the activity is low or intermediate but become nonlinear at high activity and show behavior similar to that found for the plastic depinning of solids. At higher drives, the active clustering is lost. For low activity, we also find a reentrant fluid phase, where the system transitions from a high mobility fluid at low drives to a clogged state at higher drives and then back into another fluid phase at very high drives. We map the regions in which the thermally clogged, partially clogged, active uniform fluid, clustered fluid, active clogged, and directionally locked states occur as a function of disk density, drift force, and activity.