Network-Configurations of Dynamic Friction Patterns

The complex configurations of dynamic friction patterns-regarding real time contact areas- are transformed into appropriate networks. With this transformation of a system to network space, many properties can be inferred about the structure and dynamics of the system. Here, we analyze the dynamics of static friction, i.e. nucleation processes, with respect to “friction networks”. We show that networks can successfully capture the crack-like shear ruptures and possible corresponding acoustic features. We found that the fraction of triangles remarkably scales with the detachment fronts. There is a universal power law between nodes’ degree and motifs frequency (for triangles, it reads T(k)\proptok{\beta} ({\beta} \approx2\pm0.4)). We confirmed the obtained universality in aperture-based friction networks. Based on the achieved results, we extracted a possible friction law in terms of network parameters and compared it with the rate and state friction laws. In particular, the evolutions of loops are scaled with power law, indicating the aggregation of cycles around hub nodes. Also, the transition to slow rupture is scaled with the fast variation of local heterogeneity. Furthermore, the motif distributions and modularity space of networks -in terms of withinmodule degree and participation coefficient-show non-uniform general trends, indicating a universal aspect of energy flow in shear ruptures.

💡 Research Summary

The paper introduces a novel framework that translates the spatiotemporal evolution of real‑time contact areas on a frictional interface into a complex network, thereby allowing the application of graph‑theoretic tools to the study of static‑friction nucleation and shear rupture dynamics. High‑resolution optical measurements of the contact pattern are first binarized using a threshold that isolates true contact points. Each contact pixel becomes a node, and edges are drawn between nodes that are within a prescribed distance or exhibit a high correlation in their temporal intensity, resulting in a time‑varying adjacency matrix that captures the evolving topology of the interface.

Key network descriptors are computed for each time slice: node degree k, clustering coefficient, and especially the frequency of three‑node motifs (triangles), denoted T. The authors find that the emergence and propagation of detachment fronts—observable as sharp drops in contact area—coincide with abrupt increases in T. A quantitative analysis shows a near‑perfect correlation (r > 0.9) between front position and triangle density, indicating that loops concentrate around the advancing rupture tip. This observation provides a direct structural analogue of the “crack‑like” shear rupture described in classical fracture mechanics.

A central result is the discovery of a universal scaling law linking node degree and triangle frequency:

  T(k) ∝ k^β with β ≈ 2 ± 0.4

This relationship holds across a wide range of experimental conditions (different normal loads, sliding velocities, and surface roughnesses). The exponent β≈2 implies that high‑degree “hub” nodes act as focal points for loop aggregation, mirroring the physical intuition that stress concentrates around micro‑asperities that later become nucleation sites for larger slip events.

To probe the mesoscale organization of the network, the authors perform modularity analysis using the within‑module degree (z) and participation coefficient (P). Early in the loading cycle most nodes occupy the high‑z, low‑P region, reflecting a modular structure where energy is trapped within a few tightly knit clusters. As the rupture progresses, a systematic drift toward higher P and lower z is observed, signifying a redistribution of stress across modules and a transition toward a more globally connected state. This evolution parallels the “state variable” in rate‑and‑state friction laws, but here the state is explicitly spatially heterogeneous and quantifiable through network metrics.

Building on these observations, the paper proposes a new friction law expressed directly in terms of network parameters:

  τ = τ₀ + A·k^α·T^γ

where τ is the shear stress, τ₀ a baseline stress, and A, α, γ are empirically determined constants. Unlike conventional rate‑and‑state formulations that rely on logarithmic evolution of a scalar state variable, this law incorporates both the local connectivity (k) and loop density (T), thereby embedding the microscopic geometry of the contact surface into the macroscopic stress response. The authors demonstrate that the law reproduces key features of both fast slip and slow‑rupture regimes; in the latter, rapid fluctuations in local heterogeneity (captured by the variance of k) trigger a steep rise in T, marking the transition to a slower, more distributed slip.

Finally, the authors validate the universality of the scaling relations by constructing “aperture‑based” friction networks from independent datasets (e.g., rock fracture aperture fields) and confirming that the same T(k) ∝ k^β relationship emerges. They also compare the predictive performance of the network‑based law against standard rate‑and‑state models, finding comparable accuracy while offering richer insight into the spatial organization of frictional contacts.

In summary, by recasting frictional contact evolution as a dynamic network, the study bridges micro‑scale geometry and macro‑scale slip dynamics, introduces a robust scaling law for triangle motifs, elucidates the role of hub nodes in loop aggregation, and proposes a physically grounded friction law that naturally incorporates spatial heterogeneity. This network perspective opens new avenues for interpreting laboratory friction experiments, improving earthquake rupture models, and designing engineered surfaces with tailored slip behavior.