Papers

We characterise the interplay between cation-translational and anion-rotational dynamics on the model system $\rm \alpha\text{-}Li_{2}SO_{4}$ through virtual isotopic substitution molecular dynamics (VISMD) simulations. We independently control in turn the lithium translational- and sufate rotational- diffusivity by varying the masses of lithium and oxygen, respectively, showing that both diffusivities exhibit power-law relationships with the masses of each individual component. Changing either diffusivity produces an effect on the other; i.e. there exists a bidirectional coupling between anion and cation diffusive dynamics. From our variable-temperature/variable-mass dataset we demonstrate that the apparent activation energy for lithium diffusion depends strongly on the isotope mass, decreasing with accelerated sulfate dynamics; an observation which provides clear evidence of mode-coupling between cation and anion dynamics.

Phase-field modeling has achieved great success in predicting pattern formation in materials, such as the formation of ferroelectric domains. However, because it is typically based on continuum mechanics, conventional phase-field modeling cannot be straightforwardly applied to atomic-scale pattern formation, such as dopant segregation and vacancy ordering, which are driven by chemical potentials. Here, we extend the phase-field concept to the atomic scale by formulating the free energy of atomic systems within stochastic self-consistent harmonic approximation (SSCHA) theory to allow fractional site occupations. Our methodology enables us to directly calculate the derivative of the free energy with respect to site occupation and thereby obtain the chemical potential, successfully reproducing Ag distributions in bulk Cu as well as the resulting lattice expansion. Furthermore, we applied our methodology to investigate dopant segregation around a {\Sigma}5(310)[001] Cu grain boundary doped with Ag atoms. We found that Ag atoms preferentially segregate at the vertices of the triangular motif of a grain boundary. As the number of dopants increases, excess Ag atoms segregate near the vertices and then at the bottom sites of the triangular motif. This study extends the phase-field concept to discrete atomic systems, enabling the identification of preferential dopant-segregation sites and thereby visualizing atomic-scale pattern formation.

Hyperuniformity in active matter is usually treated as a property of a prescribed density or continuum field. This view misses a basic feature of hydrodynamic active matter: an incompressible fluid does not respond equally to every microscopic force. Longitudinal forcing is absorbed into pressure, whereas transverse forcing drives flow. The relevant question is therefore not only whether particles are uniformly arranged or whether the total activity is small, but which sector of the active forcing is selected by the physical response. Here we introduce response-selected hyperuniformity, in which long-wavelength order is a property of a source-response pair. In a reversible valence-one fluid with no prescribed partners, locally neutral clusters screen the signed active-moment sector that controls transverse flow, producing a first-moment spectrum that vanishes quadratically at low wavenumber. Locally unscreened moments instead generate a nonzero infrared plateau. The resulting transverse-force spectrum has a universal crossover from fourth- to sixth-order scaling, with the crossover set by the ratio of the unscreened residual to the screened analytic contribution. Complete partner renewal preserves this normal form, establishing exchangeable multipole inheritance, while turnover tunes the residual through an independently measured local defect density. The zero-residual limit yields strictly hyperuniform velocity fluctuations; any finite residual causes defect-controlled infrared leakage and sets a finite screening length. Thus microscopic exchange need not destroy hidden hyperuniform flow order, but rare unscreened moments determine how far the quiet-flow regime survives.

Within the framework of stochastic first-passage time thermodynamics (TFPT) and power-law distributions (Mittag-Leffler and Levy), the appearance of fractional Caputo derivatives in hydrodynamic equations is substantiated. The developed framework is applied to the problem of explosive boiling of a superheated liquid: an integral formula for nonlocal entropy production is derived and the critical index of dissipation intensity is calculated, along with qualitative changes in these parameters upon introducing external pressure. It is found that dissipation at a given instant is strictly determined by the integral history of the macroscopic matter flow J(t) over the entire evolutionary interval up to the moment of phase explosion. The calculated critical dissipation index z strictly describes the kinetic catastrophe and the violation of Prigogine's minimum entropy production principle near spinodals.

Strong electronic correlations in narrow-band systems provide a promising route to realize emergent quantum phases. While ferroelectricity in van der Waals materials is typically associated with inversion symmetry breaking driven by lattice distortions, interlayer sliding, or moiré reconstruction, the possibility of generating ferroelectricity directly from electronic interactions remains largely unexplored. Here, using molecular beam epitaxy, scanning tunneling microscopy, and ab initio calculations, we investigate two stacking geometries of bilayer 1T-TaSe$_2$, A-C and A-C$'$, formed by coupled Star-of-David charge density wave phases. We show that both stackings realize quasi-one-dimensional interacting chains, but are governed by distinct interaction mechanisms. In the A-C stacking, strong interlayer hybridization leads to dimerization and the formation of a band insulating state. In contrast, the A-C$'$ stacking is dominated by interlayer Coulomb interactions, producing a spontaneous charge imbalance between layers that gives rise to an out-of-plane ferroelectric polarization. Furthermore, we demonstrate that ferroelectric and antiferroelectric interchain configurations can be stabilized and electrically switched by an external field. Our results prove that bilayer 1T-TaSe$_2$ is a platform for interaction-driven electronic ferroelectricity, establishing an overlooked family of charge-ordered correlated states in 1T-TaSe$_2$ multilayers.

Ultrafine-grained (UFG) aluminium alloys are promising lightweight structural materials for space applications, where a high grain-boundary density can act as sinks for irradiation-induced defects. Their deployment, however, is contingent on thermal stability: aluminium components in low-Earth orbit can reach $\sim$200 $^\circ$C under solar irradiation, close to where severely deformed aluminium alloys recrystallise. Accurate, bulk-representative determination of recrystallisation onset is therefore essential, yet conventional in situ transmission electron microscopy (TEM) heating is compromised by thin-film effects, ambiguous grain-boundary contrast, and small sampling volumes. Here, a UFG AA6061 (Al-Mg-Si) alloy produced by high-pressure torsion was studied by a direct comparison of in situ TEM heating and in situ electron backscatter diffraction (EBSD) heating, complemented by differential scanning calorimetry (DSC), analytical scanning transmission electron microscopy (STEM-EDX) and microhardness. In situ EBSD sampled $\sim 10^{3}$ grains from bulk material and resolved the microstructural evolution into sequential recovery, recrystallisation and grain-growth regimes, placing the onset of instability at $\sim$198 $^\circ$C. Calorimetry, microhardness and nanoscale elemental mapping showed that grain refinement suppresses GP-zones formation and shifts precipitation to lower temperatures, with precipitation neither retarding recrystallisation nor restoring strength once the UFG structure is consumed. Revisiting the Brailsford-Bullough-Hayns sink-strength theory with a KAM-informed, temperature-dependent internal sink strength, we show that the grain-boundary sink strength collapses as recovery and recrystallisation proceed. We establish in situ EBSD heating as an in operando method for bulk-representative determination of microstructural instabilities in advanced metallic systems.

The theoretical foundation of combining external nonequilibrium thermodynamics, the nonequilibrium statistical operator method, and stochastic first passage time thermodynamics are explored. It is shown that including a random lifetime of a metastable state in the generalized distribution function allows the first passage time to be considered as a fully - fledged macroscopic coordinate. A microscopic justification for this approach is provided, and a generalized thermodynamic potential is introduced. A procedure for closing the transport equations based on on thermodynamic consistency conditions is developed. It is shown that in the generalized Maxwell - Catteneo equation, the classical relaxation time is strictly replaced by the mean first passage time, which imparts internal macroscopic nonlinearity to the system. Using renewal theory, the exact mathematical structure of transport memory kernels for non-Markovian processes is derived. A qualitative comparative analysis of various approaches to the thermodynamic description of nonequilibrium phenomena is presented.

Systems composed of interacting bistable elements are commonly described by transition graphs that determine which state changes are accessible under an external drive. Under quasistatic loading, accessibility is constrained by the equilibrium structure of the system, often resulting in sparse transition networks and unreachable stable states. Here, we show that dynamic loading of dissipatively-coupled hysteron networks enhances accessibility by enabling transition pathways that are forbidden under quasistatic driving while preserving the underlying equilibrium states. In particular, we consider pulse actuation and derive a control map linking pulse amplitude and duration to state transitions. For suitable dissipative couplings, individual hysterons become independently addressable using a single scalar input, increasing transition-graph connectivity and enabling access to otherwise unreachable states. We validate the framework experimentally using pneumatic hysterons and find good agreement with theory. More generally, the framework applies to dissipatively coupled networks of bistable elements across fluidic, mechanical, and electrical domains.

The quantum dynamics of a double quantum dot two-level system, coupled to a leaking microwave resonator mode, is theoretically investigated. The double quantum dot is driven by applying a continuous field giving rise to time-dependent modulations of its energy separation between the ground levels of each dot, forming the qubit. Therefore, the cavity resonances occur when the difference between the resonator and the qubit frequencies equals the multiple of the modulation frequency, respectively. In the strong Coulomb interaction limit and the weak couplings of the combined system, i.e. the qubit and the single-mode cavity, to corresponding electronic, phononic or photonic reservoirs, we have demonstrated that the output cavity electromagnetic field is formed of a single-photon flux obeying the sub-Poissonian quantum photon statistics. This was proved via Fano's factor behaviour or by comparing the resonator's second- and third-order photon correlation functions, respectively. The phonon, the cavity mode dephasing or the environmental temperatures influence on the quantum properties of the photons are also discussed.

Coplanar magnetic spirals map locally onto uniform spin–orbit-coupled wires, but on a ring the rotating spin frame can change the electronic boundary condition. We prove, level by level in every fixed-particle-number sector, that even and odd texture windings impose boundary phases separated by half the single-electron flux period. Changing the winding by one also changes the physical pitch in inverse proportion to the circumference; a symmetric comparison removes this leading pitch correction. In the number-conserving topological phase, a global parity constraint then reverses the fermion-parity assignment of neighboring phase-winding branches. Density-matrix renormalization-group simulations show this reversal, support an Ising transition coexisting with a gapless charge mode, and relate the finite-size parity splitting to an independently calculated charge stiffness. Magnetic-texture winding thus changes the many-body spectrum through a global boundary condition that is not determined by the local band dispersion, providing a closed-geometry, number-conserving probe of topological pairing.

The Fe-Cr binary alloy is a model system for ferrtitic/martensitic steels. Cr clustering instigated loss of ductility, is one of the influential factors to deteriorate the alloys, and its severity depends on Cr concentration. This study aims at understanding (i) the nature of Cr precipitation in Fe-12Cr alloy using positron lifetime spectroscopy, transmission electron microscopy and Mössbauer spectroscopy studies and (ii) comparing the precipitation concentration and size with respect to Fe-9Cr alloy. The quenched alloys (from the high temperature FCC phase) are found to contain a high density of dislocations along with the sub-grain structure which lead to Cr rich {\sigma}-phase precipitation during post-quenching heat treatment at 748 K, while the alloy prepared with no dislocations and the sub-grain structure shows no precipitation, except at grain boundaries, consistent with observations in Fe-9Cr alloys. Selected area diffraction confirmed that the precipitates are tetragonal intermetallic {\sigma}-phase, contradicting the majority of the literature on Fe-Cr ferritic alloys which reported it as Cr rich BCC {\alpha}'-phase. Mossbauer spectroscopy studies showed signatures of Cr depletion from the rest of the matrix, indicating that the precipitates are Cr rich compared to the alloy composition. The kinetics of Cr precipitation is found to be higher in the Fe-12Cr alloy compared to the Fe-9Cr alloy. The average precipitate size is observed to be higher and the number density is lower in the Fe-12Cr alloy compared to the Fe-9Cr alloy. In addition, the experimental positron lifetime in defect-free lattice of Fe-(9, 12)Cr alloys is reported for the first time with substantial evidence.

At non-time-reversal-invariant momenta (non-TRIMs), time-reversal symmetry does not constrain the linear term of the band dispersion. Whether $\nabla E$ vanishes is therefore determined entirely by the representation theory of the little group. For nondegenerate bands, $\nabla E$ is forced to zero if and only if the vector representation $\Gamma_{\mathrm{vec}}$ of the little group does not contain the trivial representation $\Gamma_1$. When $\Gamma_{\mathrm{vec}}$ does contain $\Gamma_1$, $\nabla E$ is not forced to vanish for any nondegenerate band; the classification instead depends on the multiplicity of $\Gamma_1$ in $\Gamma_{\mathrm{vec}}$. For degenerate bands, the Wigner-Eckart theorem and Clebsch–Gordan coefficients determine whether linear couplings vanish, with classification performed at the subband level. Applied to space group 225, the criterion explains why the $W$ point is critical for all nondegenerate bands, the degenerate $E$ bands are generically noncritical, and the $K$ and $U$ points host parameter-dependent criticality. Supporting phase diagrams reveal a two-tier hierarchy: symmetry enforces $\nabla E=0$, while band parameters determine higher-order character. We extend this classification to all space groups hosting non-TRIMs in the single-group limit, providing a symmetry-dictated, parameter-independent framework for engineering Van Hove singularities in three-dimensional quantum materials.

Although parity-time (PT)-symmetric Hamiltonians are often associated with real energy spectra, PT symmetry is neither a sufficient nor a necessary condition for a real spectrum. More generally, real spectra are associated with the broader class of pseudo-Hermitian Hamiltonians, of which PT-symmetric Hamiltonians constitute a simple subclass. Here, we investigate the nonequilibrium dynamics of the imbalanced-pairing Kitaev model, a prototypical pseudo-Hermitian system, under a linearly time-dependent chemical potential. The dynamics are analyzed within the biorthogonal framework using the concept of dynamical quantum phase transitions (DQPTs). We show that, under a linear ramp protocol, DQPTs occur only when the post-ramp Hamiltonian possesses a real energy spectrum. For positive values of the non-Hermiticity parameter ($\gamma>0$), where the energy spectrum remains entirely real, a ramp crossing a single quantum critical point gives rise to a single family of critical times, analogous to the Hermitian case. Furthermore, for ramps crossing two critical or exceptional points, the critical sweep velocity above which DQPTs disappear decreases as the non-Hermiticity parameter is reduced and vanishes in the staggered-pairing limit, $\gamma=-1$.

Optical nanoscopies including near-field optical microscopy and spectroscopy circumvent the diffraction limit of conventional optics thanks to the nanoscale light focus emerging at the apex of a sharp irradiated probe. However, while strong optical coupling between the apex and its dielectric environment affords both enhanced nanoscopic measurement sensitivity and potentially ultra-strong fields within a nano-gap cavity, the conditions for this coupling remain poorly quantified by prevailing analytic models. Here we present a robust formalism of probe-cavity eigenmodes that fully describes how mutual near-field interactions between probe and environment produce a composite response to external fields qualitatively distinct from that of its distinct components. This "EigenProbe" model identifies the fundamental excitations of realistic optical nanoscopies as nano-gap polaritons, which provide an elegant basis to accurately predict near-field microscopy and spectroscopy experiments especially when probe-sample interactions are non-perturbative. Through comparison to carefully controlled nanoscopies of polar phonons and molecular vibrations alike, we show how nano-gap polaritons are both realized and utilized for reliable and rapid "inversion" of local optical constants. This advance demands both a careful understanding of the probe response through quantitative calibration, and our efficient semi-analytic description of cavity eigenmode scattering at the probe apex. Our EigenProbe formalism sets the stage for maturing diverse and proliferating optical nanoscopies into precision metrologies of nano-scale optical environments, and guides future use of nano-gap cavities to manipulate local excitations of quantum materials and to achieve strong coupling over photonic emitters.

Second-order nonlinear optical processes in van der Waals (vdW) semiconductors offer a compelling route toward compact, integrable photon-conversion platforms. Among III-VI vdW semiconductors, GaS is particularly attractive owing to its wide bandgap suppressing two-photon absorption under near-infrared laser excitation. However, bulk GaS typically crystallizes in the centrosymmetric $\beta$ phase, which eliminates second-order nonlinearity and severely limits its application in nonlinear photonics. Here, we demonstrate that GaS nanobelts synthesized via self-catalyzed vapor-liquid-solid growth predominantly crystallize in non-centrosymmetric $\gamma$-phase stacking. This behavior originates from edge-selective growth kinetics at the Ga catalyst interface, which stabilizes the $\gamma$ phase and enables deterministic in-plane dipole moment alignment. The GaS nanobelts exhibit strong second-harmonic generation (SHG) with intensities comparable to those of GaSe, a widely used nonlinear optical material. Furthermore, we integrate the nanobelts onto SiN waveguides and demonstrate efficient on-chip SHG and sum-frequency generation. These results establish $\gamma$-GaS nanobelts as a transferable one-dimensional nonlinear materials well suited for on-chip photonic integration and indicate their strong potential for nonlinear optical applications.

Mechanical instabilities provide a general design principle for shaping developing organs and engineering soft materials. However, in slender structures, simple elastic buckling tends to erase rather than preserve shape complexity: structures relax to the simplest possible shape, erasing finer detail. Living systems nonetheless build complex, reproducible morphologies from continually remodeling material, while remaining robust to noise arising across scales. Using analytical theory and numerical simulations of a minimal model of growing visco-elasto-plastic rods, we show that remodeling plays two opposing roles: At low plasticity, patterns coarsen through elastic relaxation, while high plasticity converts fluctuations into geometric disorder. This sets a trade-off between shape complexity and reproducibility with an optimal intermediate plasticity, which protects initial patterns. Growth breaks this trade-off by suppressing both failure modes, enabling complex shapes to be reproducibly generated. Our results identify remodeling and growth rates as two knobs governing whether an encoded pattern is remembered, degraded, or transformed.

Physical reservoirs, neuromorphic devices, and wave-mediated systems often possess memory, feedback, and rich state-dependent dynamics, but these properties do not by themselves establish autonomous computation. Here we develop a closure criterion for autonomous physical computation, motivated by the wave–particle walker. We formulate the walker as a stroboscopic reservoir with state, where the wave field stores an exponentially decaying trace of previous droplet impacts and guides future motion through local slope coupling. This model separates physical writing, storage, reading, feedback, and externally triggered erasure. We then define computation as robust coarse-grained transition preservation: a physical map implements an abstract transition only when a coarse-graining satisfies compositionality, with abstract states realized by separated physical basins and transitions stable under noise. Autonomous physical computation requires a further closure condition: an internal physical readout state must select the next physical operation. This criterion classifies the wave–particle walker as a wave-memory machine with genuine Turing-like primitives, but not as a closed autonomous physical computer, because the erasing phase shift is externally imposed. The framework turns this distinction into a design principle: memory becomes autonomous computation when physical readout basins are coupled back to operation selection.

X-ray absorption spectroscopy (XAS) is central to understanding the local electronic and atomic structure of materials, yet most published spectra remain inaccessible to data-driven analysis because they are embedded in figures and described through fragmented textual context in the literature. Here, we use multimodal (image and text) literature mining to transform this dispersed knowledge into an AI-ready experimental data resource. We developed a scalable spectroscopy data digitization pipeline that identifies XAS figures in full-text articles, digitizes spectral curves, and links each spectrum to accompanying metadata on the measured edge and material. Applying this pipeline to the battery literature produced an open dataset of 13,740 XAS spectra, spanning 66 absorbing elements and diverse battery chemistries, with expert validation confirming accurate extraction of spectral and metadata information. By converting literature-embedded spectra into structured numerical data, this dataset provides a foundation for large-scale XAS analysis, cross-laboratory comparison, high-throughput characterization, and autonomous discovery of advanced materials.

We present a temperature-dependent infrared study of intrinsic and Zn-doped GaAs within a generalized dispersion framework describing coupled phonon and free-carrier excitations. Reflectance, transmittance, and ellipsometric data measured from 300 to 440 K are fitted simultaneously using Gaussian-broadened distributions for single-phonon modes and triangle-delta-triangle representations for multiphonon absorption bands. The formalism provides a self-consistent description of lattice and carrier contributions to the dielectric response and reproduces the evolution of phonon-plasmon coupling with doping and temperature. In heavily doped samples, a pronounced Fano-type asymmetry of the TO-phonon resonance is observed and attributed to interference between phonon excitations and the Drude background of light and heavy holes. The resulting optical constants remain physically consistent across the far- and mid-infrared range and demonstrate the applicability of the framework to coupled bound and free excitations in doped semiconductors. The approach provides a quantitative tool for modeling infrared optical properties of GaAs and related photonic materials.

We investigate the Uhlenbeck-Ford (UF) model as a reference system for free-energy calculations in two-dimensional (2D) fluids. The 2D virial coefficients are computed exactly up to tenth order and combined with molecular simulation data to construct highly accurate numerical representations of the equation of state and the excess Helmholtz free energy. We then determine the phase diagram of the model in order to establish the thermodynamic stability limits of the fluid phase and thereby identify the range of applicability of the UF model as a fluid reference system. In the course of this analysis, we identify the solid, hexatic, and fluid phases, and show that the fluid remains the only thermodynamically stable phase, independent of density, for scaling parameters up to $p\lesssim 70$. Finally, we demonstrate the practical applicability of the 2D UF model as a reference system through thermodynamic integration calculations of the free energy of a two-dimensional Lennard-Jones fluid.

We present a method that enables efficient simulations of coarse-grained electrolyte solutions inside a narrow slit pore with an anisotropic dielectric permittivity tensor. The electrostatic equations for polarizable surfaces are solved using a 2D periodic Green's function method combined with a slab-corrected anisotropic 3D Ewald summation. We apply this approach in Monte Carlo simulations to study 1:1 electrolytes confined between both polarizable dielectric and metallic surfaces. Our results show that dielectric anisotropy aggressively reshapes the double-layer structure. While a coordinate stretching transformation demonstrates that individual ion-image interactions depend strictly on the bulk-like parallel permittivity, the suppression of the perpendicular permittivity dramatically amplifies direct in-plane ion-ion correlations. Under strong anisotropy, these lateral correlations dominate the thermodynamics completely, rendering the structural profiles of mutually opposing dielectric and metallic boundaries practically identical by forcing the smaller cations directly into the contact plane of the larger anions.

The discovery of scaling laws has motivated training neural networks on ever increasing quantities of data. This is typically done with a constant decoupled weight decay which causes the network weights to shrink steadily over the course of training. Taking inspiration from the Robbins–Monro conditions, we propose to scale weight decay by the fraction of the peak learning rate $\eta/\eta_{\max}$. We prove that this scaled weight decay preserves the asymptotic stationarity guarantees of the corresponding unregularized methods for both stochastic gradient descent and the non-Euclidean spectral optimizer Muon, thereby avoiding the additional asymptotic bias introduced by constant decoupled weight decay. This retains the stability benefits of weight decay without changing the asymptotic optimization target. Using a steady-state analysis, we explain why under standard weight decay the weight norm shrinks steadily as training proceeds, whereas under scaled weight decay it settles to a roughly constant value. When applied to the training of mixture-of-experts models, Muon with scaled weight decay (Muon-SW) consistently outpaces Muon with identical hyperparameters, reaching the same validation loss $\mathbf{30\%}$ faster at our largest scale across models from $72 - 930$ million parameters trained at $\sim 600$ tokens per active parameter. If this trend continues to hold, the method promises to substantially accelerate the pre-training of frontier models while requiring only a few lines of code to implement.

Strongly disordered solids exhibit distinctive properties at the nanoscale, where conventional continuum theory breaks down. We show that their disorder-averaged linear response can be described by a modified continuum theory in which the response coefficients are local but depend nonlocally on the structural properties. The stability criterion requires the response operator to be positive semidefinite, which naturally leads to a correlated Wishart disorder. In the limit of strong disorder, the resulting theory reduces the long-wavelength response to a scalar disorder-induced contrast field. In elasticity, this field describes the formation of a stiffened interphase around rigid nanoparticles and boundaries, whose characteristic extent is set by the nonaffine length. The predictions are confirmed by molecular-dynamics simulations of a Lennard-Jones glass and a model polymer. Direct calculations of the nonlocal elastic kernels show that their spatial range is much shorter than the extent of the stiffened interphase, demonstrating that the latter does not require long-range constitutive nonlocality.

We study the Berry curvature effects of the Bogoliubov quasiparticles in chiral superconducting rhombohedral graphene. Using a two-band Bogoliubov-de Gennes Hamiltonian to describe the superconducting quasiparticles, we calculate the momentum-space Berry curvature, orbital magnetic moment, and anomalous thermal Hall, spin Nernst, and orbital Nernst transport of chiral $p$-wave superconducting states. We investigate the impact of normal-state band warping in rhombohedral graphene, which can generate Bogoliubov Fermi surfaces for quasiparticle excitations. We find that Bogoliubov Fermi surfaces qualitatively modify the anomalous transport responses, inducing a deviation of the thermal Hall conductivity from the quantized value and strongly enhancing the spin and orbital Nernst responses.

We consider the spatial spread of epidemic into an unstable, healthy state. When the transmission rate depends on the fraction of infected, the propagating pulse of infection can be either a pulled front or a pushed front. We determined the phase space of parameters for the pulled and pushed regions both by numerically solving the spatial SIR partial differential equations and by a theoretical analysis of the front propagation phenomenon. We found both a continuous and a discontinuous transition between the pulled and pushed front solutions; in the latter transition, the front speed undergoes a jump as a certain parameter crosses the critical threshold. The behavior of the propagating pulses near the transitions has been analyzed and a good agreement between the theory and numerics has been observed. Finally, a bistable region where both pushed and pulled propagating pulses can be realized has been discovered.

Phase-change memory (PCM) is a mature technology for fast, scalable, non-volatile data storage, with applications spanning embedded memory, as well as in-memory and neuromorphic computing. PCM predominantly relies on chalcogenide alloys, with $\mathrm{Ge_2Sb_2Te_5}$ (GST) as the industry standard. Yet in these alloys, the individual Ge, Sb, and Te atoms redistribute upon cycling, causing stochastic operation and ultimately device failure. To address this issue, elemental antimony was proposed as a PCM material, but it exhibits a metastable amorphous state that prevents reliable data retention. Moreover, tellurium and antimony can contaminate complementary metal-oxide-semiconductor (CMOS) production lines or act as unintended dopants, restricting manufacturing of PCM to dedicated fabs. Here we introduce elemental germanium (Ge) as a CMOS-native phase-change material that overcomes these fundamental limitations. In a vertical PCM cell architecture, Ge enables sub-nanosecond crystallization (240 ps, 40 times faster than GST), non-volatile data storage with excellent thermal stability ($>$110 {\deg}C for 10 years vs. $\sim$87 {\deg}C for GST), and a resistance drift coefficient approximately 60% lower than in GST. These results establish pure Ge, a standard semiconductor, as an alternative to chalcogenide phase-change materials, achieving superior performance in key metrics and enabling phase-change memory to be fabricated in standard semiconductor facilities.

Engineering interactions that stabilize non-Abelian fractional quantum Hall phases is a central challenge in strongly correlated topological matter and quantum simulation. We introduce a differentiable framework for inverse Hamiltonian design, in which Haldane pseudopotentials are optimized by gradient-based exact diagonalization to stabilize target fractional quantum Hall phases. In spherical geometry, the Haldane pseudopotentials are treated as variational parameters and optimized in a JAX-based exact-diagonalization framework. By directly maximizing the overlap between the many-body ground state and the Moore-Read state, we obtain a robust pseudopotential profile that has Pfaffian overlaps exceeding $99\%$ for systems up to $N_e=12$, substantially improving over conventional Coulomb interactions. Analyses of the neutral excitation spectrum and orbital entanglement spectrum further confirm that the optimized interaction stabilizes the Pfaffian topological phase. Our results demonstrate that essential features of the three-body Pfaffian parent Hamiltonian can be effectively encoded in a suitably designed two-body interaction. Furthermore, they identify a nearly universal exponentially decaying pseudopotential profile that stabilizes the Pfaffian phase and establishes a general framework toward engineering non-Abelian topological order in quantum simulation.

We present a tunable mesoscale model to provide a basis for future rheological calculations of of bacterial biofilms, explicitly incorporating reversible crosslinking within the extracellular polymeric substance (EPS) matrix. Using a Dissipative Particle Dynamics framework combined with a Gillespie-inspired algorithm, bonds between polymers and bacteria dynamically form and break, capturing the intrinsically evolving nature of the network. We show that biofilm structure is governed by a competition between polymer-polymer and polymer-bacteria crosslinks, controlled by binding energy, linker availability, and bond stiffness and provide a minimal model that helps to understand the competition between both species.

Pair-density-wave (PDW) superconductivity is a state in which the superconducting order parameter modulates at a finite wavevector. Using large-scale density matrix renormalization group, we study the doped altermagnetic Mott insulator in the checkerboard $t$-$J$ model, where altermagnetic exchange anisotropy is encoded microscopically through anisotropic ferromagnetic next-nearest-neighbor exchange. By mapping the ground-state phase diagram as a function of doping and altermagnetic anisotropy, mainly on six-leg cylinders, we identify a transition from a uniform $d$-wave superconducting regime with charge modulation to a PDW regime coexisting with stripe order. In the PDW regime, we report an unconventional wave-vector locking $\mathbf Q_{\mathrm{PDW}}\approx 2\mathbf Q_{\mathrm{Stripe}}$ along the cylinder direction, in contrast to the conventional relation. Pair correlations reveal a two-scale structure, consisting of short-distance local $d$-wave pairing and long-distance finite-momentum PDW correlations. A symmetry-based Ginzburg–Landau analysis is presented for the observed locking. Our results identify altermagnetism as a strong-coupling, microscopically grounded route to finite-momentum superconductivity in doped Mott insulators.

We solve the closed-system problem of cooling $M$ qubits, selected from $N$ identical thermal qubits, to the lowest common local temperature allowed by unitarity. The optimal protocol consists of two conceptually distinct steps. First, a passive rearrangement assigns the largest eigenvalues of the initial state to target sectors of lowest Hamming weight, thereby minimizing the total target energy. Second, a target-only complex-Hadamard transformation within each fixed-Hamming-weight subspace equalizes the one-qubit target marginals without changing any target-sector probability or the total energy. Consequently, imposing a common local temperature costs neither cooling depth nor additional work: the constrained optimum coincides with the unconstrained passive minimum for every $N>M$ and every initial temperature. The complex-Hadamard correction may nevertheless be costly at the circuit level. We therefore derive an exact arithmetic criterion for when the same optimum can be attained by a temperature-independent computational-basis permutation alone, and exhaustively classify the resulting finite-size islands of feasibility for $M+2\leq N\leq 128$. At isolated temperatures, further optimal permutations can arise through numerical cancellations between different thermal eigenvalue shells. These alternative realizations may reduce implementation complexity, but they cannot improve the cooling curve already attained by the universal protocol. We also derive the exact cooling curve, prove that at least two ancillary qubits are necessary and sufficient for nontrivial cooling, and show that joint many-target cooling can strictly outperform parallel single-target strategies.