The study, led by Lingyan Shen and Yonggui Liu with Xiahui Pan, builds on two decades of experimental work showing that the onset of frictional sliding is mediated by rupture fronts. In landmark experiments, researchers visualized the real contact area between polymethyl methacrylate (PMMA) blocks using laser illumination at angles beyond the critical angle for total internal reflection, revealing three distinct classes of detachment fronts: slow fronts, sub-Rayleigh fronts, and intersonic fronts that travel faster than shear waves but slower than longitudinal waves. Macroscopic sliding begins only after these fronts have swept across the whole interface, breaking micro-contacts one after another like a line of falling dominoes.

The new observations do not fit that sequential picture. In the experiments, two PMMA blocks forming a rough interface were loaded by a split-Hopkinson pressure bar, which delivers a well-defined compressive pulse to the trailing edge of the upper slider. Arrays of miniature piezoelectric film sensors, just 3 millimeters in active diameter with a 50-megahertz center frequency, were bonded along two orthogonal directions: one set parallel to the interface at a fixed distance of 2.5 millimeters, and another set perpendicular to it at depths ranging from 2.5 to 28 millimeters. A trigger strain gauge on the impact face defined the absolute time origin with sub-microsecond precision.

What the sensors recorded was startling. All five transducers in the parallel array detected three distinct wave packets at approximately 0.34, 0.70, and 1.00 microseconds after impact, with nearly identical amplitudes and waveforms. At the earliest of those instants, the incident compressive wave, traveling at roughly 2,520 meters per second, could have covered less than a millimeter, nowhere near the 20 millimeters separating the first sensor from the loaded edge. Simultaneous detection across the full array therefore points to a distributed, non-local source: the interface as a whole lit up at once, rather than being activated edge-first by a propagating disturbance. Control experiments, including specimen-removed electrical cross-talk tests and independent optical trigger calibration showing channel-to-channel jitter of only about 5 nanoseconds, confirmed that the synchronization was a physical reality rather than an instrumentation artifact.

The wave packets carried further signatures that set them apart from anything in the established catalog of interfacial dynamics. Their propagation speed, extracted from the trajectories of the wavefronts across the sensor arrays, matched the longitudinal wave speed of PMMA, about 2,540 meters per second experimentally. That is far above both the Rayleigh wave speed of roughly 1,237 meters per second and the shear wave speed of about 1,400 meters per second, the ceilings that classical fracture mechanics imposes on rupture fronts propagating along an interface. Moreover, the waves radiated perpendicularly into the bulk of the material, whereas rupture fronts are confined to the interfacial plane itself. The team conducted 150 experiments varying impact velocity from 5 to 15 meters per second, interface roughness, and loading geometry, and found that the characteristic arrival times, inter-packet intervals, and peak amplitudes varied by less than 4.2 percent across all conditions, a robustness that implicates collective interfacial dynamics rather than stochastic asperity-scale contact mechanics.

To probe the mechanism, the researchers built explicit dynamic finite element simulations in ABAQUS, modeling two elastic blocks pressed together along an interface with idealized asperity geometries: a smooth baseline, periodic isosceles-triangular asperities, and sharper right-triangular asperities. Using four-node quadrilateral elements with a refined mesh of 1-micrometer characteristic length and a time step of about 0.35 nanoseconds, the simulations reproduced the essential phenomenology. Von Mises stress contours revealed spherical wavelets emanating from micro-contact regions whose collective envelope formed planar wavefronts parallel to the interface, a structure that persisted across all three asperity configurations. The simulated propagation speed of 2,536 meters per second agreed with the experimental value to within 0.2 percent, and the simulated inter-packet interval of about 0.30 microseconds matched the experimental 0.30 to 0.34 microseconds, with the small offset attributable to the finite rise time of the real loading pulse and trigger uncertainty.

The simulations also exposed a remarkable discrete structure in the interfacial response. Monitoring six adjacent material points along the interface showed statistically indistinguishable time histories, confirming synchronous activation, while the normal stress and displacement traces revealed wave packets arriving at equally spaced intervals of roughly 0.3 microseconds. When the plateau amplitudes were normalized, they followed a striking quantization law: the amplitude at the nth plateau grows as 3n(n − 1), a closed-form rule indicating that the interface accumulates wave energy in discrete, integer-multipled steps. Spectral analysis of the oscillatory signal showed a dominant carrier frequency of about 34.4 megahertz, with higher harmonics, whose amplitude but never its frequency changed from packet to packet. The normal displacement evolved as a staircase of flat plateaus, each step coinciding with the arrival of an additional packet, echoing discrete force plateaus previously observed in stick-slip systems.

To explain these features from first principles, the authors propose a non-local micro-oscillation model that extends Huygens’ principle to elastodynamics. In this picture, the frictional interface behaves as a distributed secondary source: once the incident plane wave arrives, every material point on the interface acts as a point emitter of spherical wavelets, and the envelope of those wavelets is precisely the observed planar longitudinal wavefront. The model’s Green’s function formalism naturally yields the characteristic 1/r geometric decay of the radiated amplitude, which the finite element results match closely; the experimental amplitudes fall systematically lower, a discrepancy the authors attribute to viscoelastic damping in the real PMMA that the conservative theory deliberately excludes. A Taylor expansion of the discrete-time source term reproduces the step-wise amplification, with the theory predicting a proportionality coefficient of π where the simulations give 3, a modest deviation likely arising from the discrete sampling of the continuous source distribution.

Perhaps the deepest implication concerns the foundations of solid mechanics itself. Classical continuum theory rests on the principle of local action, the assumption that a material point responds only to conditions in its immediate neighborhood, which underpins the standard hyperbolic wave equation and guarantees that disturbances propagate at finite speeds. The synchronous activation of the entire interface within about 0.3 microseconds, when acoustic transit across the 40-millimeter contact would take roughly 16 microseconds, suggests instead that elastic coupling among micro-asperities establishes a correlated source density across the whole interface almost instantaneously. The authors formalize this by making the source term of the bulk wave equation satisfy an elliptic constraint, a modified Helmholtz equation with an exponential-decaying Green’s function, so that the wave operator remains hyperbolic while its source behaves non-locally. The temporal discreteness, meanwhile, challenges the assumed smoothness of classical elastodynamic solutions and points toward delay-differential and pseudodifferential operator frameworks for interface-mediated wave propagation.

The practical stakes are considerable. Because the discrete wave packets are intrinsic to the interface rather than artifacts of loading or morphology, they offer a quantitative fingerprint that could underpin non-destructive evaluation of bolted joints and friction assemblies, where hidden interfacial degradation is a persistent safety concern. The authors also point to seismic precursor monitoring as a potential application, since the same fracture-wave coupling operates at laboratory faults and may scale to natural ones. Much remains to be done: a self-consistent derivation of the 0.3-microsecond characteristic timescale from material and geometric parameters is still outstanding, and interface curvature effects await study. But the core message is already clear and likely to reverberate well beyond tribology: under impulsive loading, a frictional interface is not merely a passive plane where cracks nucleate. It is a coherent, quantized radiator, a collective oscillator whose synchronized voice challenges how physicists have modeled contact for generations.

Subject of Research: Non-local, quantized stress wave emission from frictional interfaces under dynamic shock loading

Article Title: Observation of non-local waves at frictional interface

Article References: Shen, L., Liu, Y., & Pan, X. (2026). Observation of non-local waves at frictional interface. Results in Engineering, 32, Article 113150. https://doi.org/10.1016/j.rineng.2026.113150

Image Credits: AI Generated

DOI: 10.1016/j.rineng.2026.113150

Keywords: friction, rupture fronts, stress waves, frictional interfaces, PMMA, split-Hopkinson pressure bar, finite element simulation, non-locality, Huygens principle, acoustic emission, earthquake nucleation, non-destructive evaluation

The research team, led by Kunchen He, Haichun Hao and Bengao Yang, together with colleagues from institutions including Sichuan University’s State Key Laboratory of Intelligent Construction and Healthy Operation and Maintenance of Deep Underground Engineering, started from an awkward truth about existing coring technology. Pressure-holding coring tools date back to the 1960s and were refined through the Deep Sea Drilling Project, the Ocean Drilling Program and later the European HYACE and HYACINTH systems, which added the Pressure Core Analysis and Transfer System. More recently, a so-called five-preservation concept has been proposed to maintain pressure, temperature, mass, moisture and light simultaneously. But comprehensive fidelity comes at a price: costs balloon and success rates fall. Because prior work had already shown that in weakly water-bearing formations the dominant deep factors are temperature and pressure, and that these two factors interact competitively rather than additively, the researchers set out to determine exactly when each one rules the behavior of deep coal.

The experimental material came from the No. 8 cross-heading in the southern part of a colliery in the Qinshui coalfield, sampled at a vertical depth of 568.2 meters within the No. 3 coal seam of the Lower Permian Shanxi Formation. The seam sits in a weak aquifer horizon, so the natural moisture content of the fresh samples was only 2.67 to 3.18 percent, allowing the team to justify excluding pore water pressure and focusing on the coupling of confining pressure and temperature. The coal is a bright, hard black anthracite with an original gas pressure of 2.5 to 2.75 megapascals and gas contents between 20.62 and 27.49 cubic meters per ton. Proximate and elemental analyses confirmed the material’s high maturity: fixed-carbon contents near 81 percent, very low to low total sulfur of 0.41 to 0.45 percent, and average maximum vitrinite reflectance between 3.15 and 3.17 percent.

Because deep anthracite is notoriously brittle and riddled with endogenetic fissures, only 27 standard cylindrical specimens of 50 millimeters in diameter and 100 millimeters in height could be prepared from the block samples, which were sealed in bubble wrap, taped, and packed in foam-lined wooden boxes on site. Physical properties varied considerably, with densities between 1.38 and 1.5 grams per cubic centimeter and P-wave velocities spanning 1.19 to 2.43 kilometers per second, a fluctuation of up to 51.2 percent. To prevent this heterogeneity from contaminating the results, the team applied a local outlier factor algorithm using density and wave velocity as characteristic variables. The filter retained 19 statistically similar specimens, cutting the mean square deviation of density by 58.4 percent and that of wave velocity by 60.6 percent, a crucial step for a study that hinges on detecting subtle temperature-driven changes.

The core of the work was a two-factor, four-level orthogonal scheme spanning 16 experimental groups, with temperatures of 20, 50, 80 and 95 degrees Celsius crossed with confining pressures of 1, 22, 44 and 55 megapascals, chosen using formation gradients of 30 degrees Celsius per kilometer and 22 megapascals per kilometer to represent depths down to 2500 meters. Testing used the RTRX-140-65 GCTS rock mechanics system, capable of 1000 kilonewtons of axial load, 70 megapascals of confining pressure and heating to 140 degrees Celsius, with a Micro-II acoustic emission system tracking crack growth in real time through six sensors. A loading path of temperature first, confining pressure second, was adopted to protect the apparatus, with a gentle heating rate of 0.5 degrees Celsius per minute, followed by a two-hour temperature-pressure retention experiment and then real-time strain-controlled triaxial loading at 0.05 percent per minute.

The results revealed a stark antagonism. At a confining pressure of just 1 megapascal, temperature ran riot: peak strength fell by 14.73 to 69.65 percent and elastic modulus by 5.5 to 32.6 percent as temperature rose, with the decline accelerating, so that the strength loss jumped from 15.95 percent between 50 and 80 degrees Celsius to 57.65 percent between 80 and 95 degrees Celsius. Stress-strain curves at 80 degrees Celsius and above developed periodic fluctuations, the first acoustic emission event arrived progressively earlier, from 5.41 to 1.4 minutes, and the fraction of shear cracks climbed 2.8 times to 51.62 percent, signaling a shift from brittle toward ductile failure. But once the confining pressure reached 22 megapascals, the thermal fingerprint vanished almost entirely. Strength and modulus values stayed within the normal fluctuation ranges across all temperatures, acoustic emission characteristics became nearly indistinguishable, and shear crack proportions settled between 26.38 and 32.82 percent regardless of temperature.

The mechanism behind this tug-of-war turned out to be largely a battle between internal pore pressure and external confinement. During low-pressure, high-temperature retention tests, the researchers observed mass loss of 0.55 percent, a collapse of moisture content by 91.79 percent, a 22.44 percent drop in volatile components, and thermal shrink film inflated by escaping colorless gas. Drawing on the classic firedamp drainage observation that anthracite desorbs roughly 0.8 percent of its methane per degree Celsius, and applying the ideal gas law, the team estimated that pore pressure could rise more than six-fold as temperature climbs from 20 to 50 degrees Celsius. Above roughly 60 degrees Celsius, desorbing methane and water vapor generate pressures strong enough to carve gas-erosion crack networks along natural weak directions, explaining the earlier crack initiation, larger fissures and erratic strain paths. At 22 megapascals or more, high confining pressure suppresses molecular activity, inhibits methane desorption and moisture evaporation, raises effective stress, and effectively strangles the gas-erosion process before it begins.

To quantify exactly where control flips from one field to the other, the team devised high-temperature confrontation tests in which samples were heated first and then subjected to gradually increasing confining pressure while deformation was tracked. By analyzing strain rate responses, they defined the controlling confining pressure, the pressure that fully cancels the thermal damage of a given temperature, and the controlling temperature, the temperature that overwhelms a given pressure. The numbers were striking: at 95 degrees Celsius, confining pressures up to 2 megapascals were powerless, with the absolute volumetric strain rate rising 4.5-fold, but at 15 megapascals the sample stabilized with deformation below instrument accuracy. The fitted relationships, a quadratic curve for controlling confining pressure and a linear one for controlling temperature, each achieved a correlation coefficient of 0.99.

Plotting these two curves together produced the study’s headline deliverable: a two-dimensional optimization atlas that divides the deep sampling environment into three preservation zones. Above the quadratic boundary, pressure governs and temperature damage is negligible, so a coring tool needs to maintain only the pressure field. Below the linear boundary, temperature dominates and only the temperature field must be preserved. Between the curves, both fields matter and dual preservation is required. For engineers designing the next generation of fidelity coring tools, the atlas promises a direct path to lower costs, higher success rates and scientifically valid cores, because it replaces blanket five-field preservation with a tailored strategy grounded in the genuine in-situ conditions of each target depth. As China’s mines push past 900 meters into first-level high-temperature zones, that kind of targeted efficiency may determine whether the deep coal frontier can be opened both safely and economically.

Subject of Research: An experimental and theoretical study of the antagonistic temperature-pressure mechanisms controlling the physical and mechanical properties of deep coal and the resulting selection strategy for fidelity coring.

Article Title: Selection strategy for preserving physical fields in deep coal fidelity coring

Article References: He, K., Hao, H., Yang, B., Xie, J., Xu, L., Duan, H., & Gao, M. (2026). Selection strategy for preserving physical fields in deep coal fidelity coring. Results in Engineering, 32, Article 113000. https://doi.org/10.1016/j.rineng.2026.113000

Image Credits: AI Generated

DOI: Not provided

Keywords: deep coal, fidelity coring, temperature-pressure coupling, thermal damage, confining pressure, acoustic emission, methane desorption, pore pressure, triaxial testing, anthracite, pressure preservation, coring technology