Papers

Sliding friction lowers the jamming packing fraction of dense suspensions and, when activated by stress, drives non-inertial shear thickening. Additional rolling resistance lowers the jamming packing fraction further and can represent effects of particle roughness or angularity. How it changes particle motion on the approach to jamming remains unclear. Using stress-controlled two-dimensional discrete-element simulations, we compare a sliding-only suspension with systems having either uniform or surface-varying rolling friction. At fixed packing fraction, adding rolling friction changes continuous shear thickening into discontinuous shear thickening. At the same distance from the stress-dependent jamming point, $\Delta\phi=\phi-\phi_m(\sigma)$, however, the flow curves nearly collapse, showing that the rheological effect arises largely from the shift in $\phi_m$. Guided by this collapse, we compare particle kinematics in the high-stress thickened state at matched $\Delta\phi$, revealing differences hidden by the similar bulk response. Translational velocity correlations extend over several particle diameters, whereas rotational correlations remain local. Rolling friction promotes co-rotation at contact in place of strong counter-rotation and suppresses rotational relative to translational fluctuations. Rotation nevertheless becomes increasingly important near jamming in every case. The suspension with surface-varying rolling friction follows the behavior of a uniform system with $\mu_r\approx0.3$ because the coefficients sampled at contacts lie well below the surface average value, here $\mu_r\approx0.5$. Thus, $\Delta\phi$ largely organizes the shear-thickening rheology, but not the particle kinematics. These retain a distinct signature of the rolling constraint that must be considered when rolling friction is used to model rough or angular particles.

Complex metallic alloys exhibit rich disorder-driven electronic, magnetic, and vibrational behavior arising from strong chemical disorder, leading to unconventional structure-property relationships with applications in extreme-environment materials, catalysis, spintronics, and thermoelectrics. Despite their technological relevance, the microscopic nature of charge redistribution and electrostatic fluctuations in chemically disordered alloys remains incompletely understood. Here, we develop a theoretical framework that uncovers universal statistical trends of disorder-driven charge transfer and Madelung-field fluctuations. Our analytical formalism demonstrates that local charge transfer and electrostatic potentials exhibit Gaussian-like statistics and universal linear charge-potential (qV) correlations emerging directly from the underlying disorder landscape. We identify the physical origin of these correlations in the interplay between electronic screening and impurity scattering and show how their statistical properties depend on carrier density, disorder strength, and compositional complexity. We further derive scaling relations governing qV trends across binary and multicomponent alloys, including high-entropy materials. Large-supercell density-functional theory (DFT) calculations show good quantitative agreement with the predicted statistical behavior for representative metallic alloys. Our results provide a computationally efficient framework for incorporating disorder-driven electrostatic fluctuations into effective-medium electronic-structure theories and establish a statistical-physics foundation for the predictive design of complex alloys.

Freestanding thin films of quantum materials naturally develop sub-micrometer, nonuniform strain fields that strongly affect their electronic, magnetic, and structural properties in both equilibrium and nonequilibrium conditions. However, current methods to quantitatively resolve these mesoscopic features are primarily restricted to scanning probes, making it challenging for dynamical measurements such as single-shot imaging and femtosecond microscopy. Here, we present a new computational framework which we denote as bend contour electron tomography (BCET), which efficiently converts bend contours in transmission electron microscope images into quantitative two-dimensional maps of strain and topography. By iteratively minimizing a designed loss function between experimental and simulated bend contour images, BCET retrieves both surface morphology and in-plane strain tensor fields without requiring scanning or diffraction mapping. We applied BCET to freestanding SrTiO$_3$ thin films, demonstrating the successful reconstruction of the local strain distribution and curvature field with high fidelity. Our approach provides a quantitative framework for characterizing mesoscale structures in freestanding films using wide-field imaging, opening new avenues to investigate how spatial inhomogeneity governs phase transitions and nonequilibrium dynamics in two-dimensional quantum materials.

The construction of a reliable renormalization group (RG) flow for discrete gravity models that preserves their underlying symmetry group, typically $U(N)$ or $O(N)$, remains an open problem. For random matrix models, which are the focus of this paper, this symmetry is intrinsically tied to the interactions encoding the random geometry of two-dimensional quantum Euclidean spacetime. We develop a novel approach based on the introduction of a partial matrix-valued intermediate field. In the large-$N$ limit, measure concentration strongly suppresses fluctuations of its singular values, allowing it to play the role of a self-consistent background field. This provides the basis for a relational RG in which the notion of scale is dynamically induced by the effective Gaussian measure in the basis where the intermediate field is diagonal. Our construction preserves the symmetry of the original model and admits a well-defined continuum limit. We show that the resulting infrared theory is described by a three-dimensional non-local Euclidean field theory with a non-trivial Wilson-Fisher-like fixed point and a single relevant direction. Remarkably, the associated critical exponent exactly matches the standard double-scaling exponent. We finally discuss extensions of the background-field approach to other discrete gravity models, such as random tensor models, and to different symmetry groups, as well as connections with more formal RG frameworks and information geometry.

From interstitial space to luminal cavities, fluid pressure and flow can remodel, reshape and even redefine a biological tissue. Fluids can either govern or react to mechanical interactions between cells. However, measuring flows at cellular scales is difficult, which makes it challenging to understand tissue hydraulics. Here, we develop a theoretical approach that captures cellular mechanics and fluid flow in one framework. We find that hydraulics can drastically influence tissue behavior. Hydraulic coupling between cell shape and size governs a tissue's response to osmotic shock, while tuning a tissue's permeabilities can channel fluid either between or across cell membranes. In active tissues, hydraulics can suppress cell mobility to the point of fracture, where we discover a hydraulic ratchet that drives fluid out of cells to generate small luminal spaces. We find experimental evidence that hydraulics can suppress cell motion in early stage zebrafish embryos injected with a thickening agent, which indicates that hydraulics may generally govern the behaviors of many multicellular systems.

We study the imaginary-time evolution of a domain-wall state constrained to return to itself in a staggered free-fermion chain. The resulting space-time profile separates into frozen, liquid, and gas phases divided by sharp boundaries known as arctic curves. The spectral gap produces an incompressible half-filled gas phase bounded by an inner arctic curve, in addition to the outer frozen-liquid boundary. We determine both arctic curves, the thermodynamic return amplitude, and the complete equal-time correlation kernel. Correlations decay algebraically in the liquid regions and exponentially in the gas, while both arctic boundaries arise as caustics of free quasiparticle trajectories. The return amplitude and correlation kernel are respectively controlled by the determinant and inverse of the same block-Toeplitz operator. We obtain an exact matrix Wiener-Hopf factorization of this operator by reducing the problem to a scalar Riemann-Hilbert problem on an elliptic spectral curve. The factorization also yields the exact thermodynamic return amplitude: its logarithm consists of a quadratic term with explicit gap dependence and a bounded periodic theta-function correction. The resulting frozen-liquid-gas structure is a continuous-time free-fermion counterpart of that found in doubly periodic dimer models.

Engineering emergent ferromagnetism at correlated-oxide interfaces offers a powerful route to creating collective states that do not exist in the parent materials. Here, we show that interfacial valence reconstruction in NdNiO$_3$/CaMnO$_3$ superlattices generates dual-sublattice ferromagnetism involving both Mn and Ni. Depth-resolved standing-wave X-ray photoelectron spectroscopy reveals enhanced Mn$^{3+}$ character on the CaMnO$_3$ side of the interface and enhanced Ni$^{2+}$ character on the NdNiO$_3$ side, establishing the configurations required for Mn$^{4+}$-O-Mn$^{3+}$ double exchange and Ni$^{2+}$-O-Mn$^{4+}$ superexchange, respectively. At low temperature, element-specific XMCD reveals ferromagnetic responses from both sublattices, with the Ni response concentrated predominantly in the Ni$^{2+}$-derived spectral component that SW-XPS independently shows to be enhanced at the interface. Together, these results show that cooperative double exchange within CaMnO$_3$ and superexchange across the interface couple the Mn and Ni sublattices, providing a general strategy for engineering interfacial ferromagnetism in correlated oxides.

We show that Landau-Zener tunneling in a small-gap insulator produces a singular DC current–voltage relation, and correspondingly, strong high-harmonic generation in AC transport. Remarkably, high-harmonic responses decrease subexponentially with harmonic order, which is parametrically slower than in ordinary conductors and PN diodes. We further derive a scaling law for the current amplitude at frequency $n\omega$ produced by the applied electric field $E$ at fundamental frequency $\omega$. Our work offers a promising route to THz generation based on frequency multiplication in solids.

We present a molecular dynamics (MD) study of Al implantation in 4H-SiC to determine how implantation temperature and dose affect defect evolution and dopant activation during annealing. Simulations use the Gao-Weber potential with a Morse Al-SiC interaction reparameterized to density functional theory diffusion and kick-in/kick-out barriers. At Al concentrations above the saturation limit of ~2e20/cm3, implantation at 900 K promotes interstitial-rich planar clusters already during implantation. During annealing, these clusters trap Al and evolve into faulted interstitial loops, reducing substitutional Al incorporation. In contrast, lower implantation temperatures preserve stronger local disorder that is consumed during epitaxial regrowth, resulting in higher chemical activation within MD-accessible annealing times. Atomistic trajectories show a thermally activated transition to faulted loops once planar clusters reach about 60 interstitials. Rather than nucleating as a single coherent disk, several locally faulted regions form first and subsequently merge. The activation energy is ~1.1 eV for dislocation nucleation and ~2.1 eV for stacking-fault growth. Frank-type loops dominate at high temperature, while transient Shockley partials occur mainly at early stages and below 2000 K. As large planar defects form and dissolve, stable compensating Al-C complexes also emerge. These findings support the experimental hypothesis that secondary defects contribute to the reduced Al activation observed under supersaturation during annealing.

For “ideal anyons,” the Berry phase associated with a closed loop of an anyon around another is robust, i.e., independent of the size and shape of the loop, and directly yields the braid statistics. This is not the case for fractional quantum Hall (FQH) quasiparticles (QPs), which are charged and have finite size. The Berry phase $\Theta$ associated with a loop of a fractionally charged QP around another contains two contributions: $\Theta=\Theta_{\rm AB}+\Theta_{\rm braid}$, where $\Theta_{\rm AB}$, the Aharonov-Bohm phase, is proportional to the enclosed area, and $\Theta_{\rm braid}$ is the order-one contribution from braid statistics. The latter is obtained by determining $\Theta$ for a closed loop twice, once with and once without the other QP inside, and taking the difference. It is known that even small uncertainties in the enclosed areas can produce corrections to $\Theta_{\rm AB}$ larger than $\Theta_{\rm braid}$. We consider how the Berry phase depends on the shape of the QP, which, unlike its charge, is not topological and varies along the path in response to the local potential. Our calculations show that practically imperceptible differences in the shape of the outer QP in the two experiments can produce a non-universal value for the braid statistics, revealing an additional challenge for a precise determination of the braid statistics. We discuss these results in the context of interference experiments. We also note that the fractional phase jumps in these experiments can be understood without assuming QPs with sharply quantized fractional charges at the edges of the FQH system; instead, they directly measure the fractionally quantized vorticity of bulk QPs.

The coupling between layered atomically-thin materials mediated by van der Waals forces allows strong electronic correlations, unique topological signatures, and non-trivial spin textures, all of which play an important role in the creation of spin-, valley-, and orbitronic devices. Here, we demonstrate that twisted bilayer graphene encapsulated by transition metal dichalcogenides exhibits a generically non-radial spin texture yet hosts a purely collinear Edelstein effect despite this lack of radial symmetry. Moreover, we show how uniaxial pressure can be used to further tune the band structure and change the number of active Fermi surfaces without compromising the collinear response. Lastly, we illustrate how the quantum geometry changes in the encapsulated twisted graphene bilayer with larger pressures spreading the Berry curvature over large regions of the moiré Brillouin zone. These results illustrate how encapsulation and pressure can be used to drastically alter the topology and spin-charge interconversion processes of moiré heterostructures.

Entanglement and magic are resources that reveal complementary aspects of quantum many-body systems. Their interplay is captured by nonlocal magic, the magic that survives arbitrary local changes of basis. Yet their markedly different dynamical behavior leaves open how this irreducible component of magic spreads. Here we connect nonlocal magic to the capacity of entanglement, a tractable quantity measuring fluctuations of the entanglement Hamiltonian. This connection enables analytically controlled predictions across a wide range of many-body dynamics, from chaotic to integrable systems, which we investigate also using large-scale numerical simulations. In chaotic systems, nonlocal magic exhibits a transient buildup, with logarithmic growth in time followed by decay to a size-independent value. For integrable systems, we develop a quasiparticle picture in quantitative agreement with numerics, showing that the same initial growth instead leads to saturation at a value logarithmic in subsystem size. These contrasting behaviors have an operational consequence for entanglement embezzlement: the strongest scramblers embezzle only transiently, whereas free-fermionic dynamics can sustain universal embezzlement in the steady state.

We identify a mechanism in which soft nonunitary fluctuations of the pair condensate can suppress phase coherence in a two-dimensional Rashba superconductor proximate to easy-plane ferromagnetic order. This spin polarization couples linearly to the supercurrent through a Lifshitz invariant, acting as an emergent gauge field for the superconducting $\mathrm{U}(1)$ phase. Integrating out soft magnetic fluctuations renormalizes the effective superfluid stiffness and drives the system into a "smectic" superconducting state, characterized by a finite-momentum helical phase texture and a long-wavelength stiffness that vanishes in one direction. Consequently, vortices no longer have a logarithmically divergent energy cost, and the quasi-long-range superconducting order is destroyed at any finite temperature. Weak lattice anisotropy pins the ferromagnetic order, cutting off the smectic-like behaviour at long wavelengths and thereby restoring a finite effective stiffness. Nonetheless, near the onset of magnetic order, the enhanced magnetic susceptibility leads to a vortex-deconfined regime for arbitrarily weak magnetoelectric coupling. The results provide insight into how intertwined superconducting and easy-plane ferromagnetic orders can reshape phase coherence, with relevance to rhombohedral graphene proximitized by a transition-metal dichalcogenide substrate.

We introduce a bootstrap method for deriving bounds on disorder-averaged observables in classical statistical systems with quenched disorder. The three main ingredients of the bootstrap are the positivity of ensembles of probability measures, disorder-averaged equations of motion, and explicit expectation values of random coupling variables. Together, these ingredients lead to a linear programming problem over the space of disorder averages. By further employing the positivity of the two-replica Gram matrix, we also formulate a semidefinite programming problem that produces bounds on two-replica correlators. We demonstrate the method using two examples: the two-dimensional random site-diluted Ising model and the one-dimensional disordered contact process. For the latter example, the bootstrap bounds exhibit a region of multiple kinks that we relate to the Griffiths region, whose finite-size scaling analysis yields a typical correlation length exponent in good numerical agreement with the known value.

Inspired by recent ARPES studies on multilayer ($n\geq3$ layers) cuprate superconductors, we use the unrestricted Hartree-Fock approximation to explore the ground-state phase diagram of two coupled, inequivalently doped square-lattice Hubbard layers. In the decoupled-layer limit, the lightly hole-doped ground state is typically an incommensurate spin-stripe state. However, with sufficiently strong interlayer coupling, stripe order is destabilized relative to a commensurate Néel-ordered metal. The resulting state exhibits a reconstructed Fermi surface with hole pockets centered at $(\pm\pi/2,\pm\pi/2)$ that are similar in character to those seen in experiments. Our results illustrate the qualitatively new physics that can arise from interlayer coupling in multilayer cuprates.

We present a bootstrap framework that yields rigorous two-sided bounds on ground-state observables of quantum systems with quenched disorder. The key idea is to extend the operator algebra to include classical disorder variables and to encode the disorder distribution through a set of moments. The method directly bounds disorder-averaged expectation values for finite or infinite lattice systems and for arbitrary disorder distributions. In contrast to conventional approaches, we do not need to sample disorder realizations, which enables us to make thorough use of symmetries that hold only on average rather than for each disorder realization. Additionally, we show that multi-replica operators can be used to access disorder fluctuations such as the variance of an observable. We apply the method to the one-dimensional random transverse-field Ising model (RTFIM) with discrete, uniform, and Gaussian disorder, obtaining tight bounds on the disorder-averaged ground-state energy density and short-range spin-spin correlators, as well as bounds on the disorder variance of the energy density. Additionally, we show that the difference between the upper and lower bounds on the disorder-averaged ground-state energy density, which we term the bootstrap gap, can be used to map out the phase diagram of the RTFIM. The bootstrap gap is small and flat within the paramagnetic and ferromagnetic phases at weak disorder, grows upon entering the Griffiths regions, and peaks along the critical line.

Grain boundary engineering is a major avenue for tailoring the mechanical behavior of polycrystalline materials. Grain boundary complexions, including amorphous intergranular films, are classically accessed through thermal driving forces and solute segregation. Here, we discover that plastic deformation can drive amorphous complexion transitions at room temperature in a chemically primed nanocrystalline binary CuZr alloy. High-resolution and four-dimensional scanning transmission electron microscopy reveal that the amorphous complexions preferentially emerge at incoherent twin boundaries. Spatially-resolved electron pair distribution function analysis at the atomic scale, the local-order characterization of amorphous complexions, demonstrates short-range and medium-range order gradients from crystal-templated interfaces to a metallic-glass-like core. We thus uncover a novel amorphous complexion transformation-induced plasticity mechanism that concurrently increases the yield strength, fracture strain, and tensile toughness about a factor of two relative to a designed reference material. Our findings establish mechanical deformation as a non-thermal pathway to trigger amorphous interfacial states for enhancing damage tolerance in nanostructured metals.

Anomalies of generalized symmetries are closely related to the statistics of topological excitations. Truncating a symmetry transformation to a finite region produces an operator that creates symmetry defects on its boundary; we call this a symmetry patch operator. Such operators resemble hopping operators that create topological excitations at the boundary of their support, but this geometric similarity does not justify identifying the two. For finite Abelian groups in one spatial dimension, we prove that if the same family of local operators defines both valid hopping operators and compatible symmetry patches, the type-III component of the anomaly is $-2$ times that of the statistics. We establish a one-to-one correspondence between statistics and anomaly by treating hopping operators as symmetric operators that commute with the global symmetry, following the description of symmetry through the algebra of local symmetric operators. We explain this correspondence through gauging by coupling boundary matter to a Dijkgraaf-Witten gauge field in one higher spatial dimension. Hopping operators act on the gauge-invariant degrees of freedom. After fixing the background gauge field, the residual gauge transformations become the global symmetry, and the gauge invariance of hopping operators becomes the requirement that they commute with this symmetry.

Quantum many-body dynamics is difficult to simulate because the global wavefunction develops extensive complexity, whereas local observables may retain only an intrinsic, system-size-independent complexity. This distinction is especially relevant in thermalizing systems, where local subsystems lose microscopic information and relax towards equilibrium. We introduce a variational time-evolution principle that exploits this separation by replacing global-state fidelity with a cost function defined on local reduced density matrices. The resulting locally optimized matrix-product-state (LOMPS) algorithm retains coherent short-time dynamics while using the simplification generated by thermalization at later times, and admits closed-form equations of motion analogous to the time-dependent variational principle. The same principle yields a quantum-classical algorithm in which a quantum processor evaluates the local cost and a noise-robust classical procedure performs the optimization. Proof-of-principle implementations on Quantinuum H2 and IBM Heron processors recover the characteristic local dynamics.

We study the Anderson model on hyperbolic lattices. The graph is the Cartesian product of the Euclidean lattice with the regular hyperbolic square lattice having five squares at each vertex. We prove the existence of absolutely continuous spectrum and the absence of singular spectrum on a deterministic set of positive measure at weak disorder, and for $0<\lambda\le2$ when the single-site density is bounded below on $[-1,1]$. Truncated Cauchy distributions and controlled perturbations give purely absolutely continuous spectrum on intervals.

Active nematics can spontaneously develop flow beyond a critical activity in confined geometries. We show analytically that the onset of flow is governed by two competing instabilities: a long-wavelength mode leading to unidirectional flow and a finite-wavelength instability producing transverse rolls. A reduced description reveals how activity, flow alignment, and nematic elasticity control the competition between these modes. We identify regimes in which the finite-wavelength instability has a lower critical activity than the long-wavelength instability, causing vortices to emerge before unidirectional flow. This establishes wavelength selection as an intrinsic feature of the onset of spontaneous flow in active nematics.

We report improved global-minimum configurations for the classical Thomson problem of $N=60$, $61$, $92$, and $99$ repulsive Coulomb charges confined to a disk. By combining the quenched molecular dynamics (QMD) method with the fixed-border heuristic introduced by Amore and Zarate, we systematically obtain configurations with energies $E_{\mathrm{QMD}}(60)=2159.3584240930$, $E_{\mathrm{QMD}}(61)=2237.19264190$, $E_{\mathrm{QMD}}(92)=5358.35353314$, and $E_{\mathrm{QMD}}(99)=6254.83029083$, which improve upon the previously best-known values. For $N=60$, the Voronoi diagram of our configuration differs from the one reported earlier. The symmetries for $N=61$ $(C_{2})$ and $N=99$ $(D_{1})$ are confirmed and are consistent with the known symmetry patterns for these configurations. For $N=92$, whereas the previously reported Voronoi diagram has lower symmetry, our solution exhibits a clear $D_{1}$ (axial) symmetry of defects. The results are highly reproducible across multiple independent runs, providing strong evidence that these configurations are robust global-minimum candidates.

Hole spin qubits in semiconductor quantum dots (QDs) enable high-fidelity all-electric control, but conventional electric dipole spin resonance (EDSR) can require substantial rf drive power at the low magnetic fields that are favorable for qubit coherence and readout. In planar Ge hole spin qubits, this can reach -27 dBm at the device, posing challenges for scalable architectures due to heating and crosstalk. Here, we demonstrate a flopping-mode (FM) qubit in Ge, where a single spin is delocalized in a double QD, combining first-order protection against charge noise with exceptionally efficient electric driving. By mapping out coherence sweet-spots as a function of magnetic field orientation we achieve $T_2^*= 1.4\mu\mathrm{s}$, $T_2^{\mathrm{Hahn}}= 11.5 \mu\mathrm{s}$, $T^{\phi, \mathrm{CPMG32}}_2= 130 \mu\mathrm{s}$, and $T_1= 226 \mu \mathrm{s}$, and a single-qubit gate fidelity of up to 99.76$\%$ for a gate time $t_{X\pi} = 88$ ns. Importantly, these results are obtained at a nearly in-plane magnetic field of 5 mT using only -52 dBm drive power at the device. We further find that qubit relaxation in this regime is consistent with a two-photon Orbach process, providing a route for further optimization. Our results demonstrate that FM-EDSR supports ultra-low-power, high-fidelity single-qubit operations, improvements that could benefit scalable hole-spin-based architectures and hybrid spin-photon interfaces.

Distinguishing altermagnetism from conventional ferromagnetism and antiferromagnetism typically relies on momentum-space probes. Here, we show that the real-space spreading of a localized spin excitation provides a distinctive dynamical fingerprint of altermagnetic order. Using linear spin-wave theory on a two-dimensional checkerboard lattice with nearest-neighbor and next-nearest-neighbor couplings J_1 and J_2, we demonstrate that finite J_2 produces direction-dependent magnon group velocities and splits the two magnon branches. The resulting propagation remains closer to that of an antiferromagnet than a ferromagnet, retaining an approximately circular outer wavefront, while the direction-dependent magnon velocities produce a pronounced cross-like spatial structure inside this front. We further show that changing the sign of J_2 to enter the unfrustrated regime (J_2<0) increases the characteristic propagation velocities and interchanges the diagonal directions of enhanced propagation, corresponding to a pi/2 rotation of the anisotropic spatial pattern. These results establish spacetime dynamics as a complementary probe of altermagnetic magnons, providing a route to identifying unconventional magnetic order through real-space propagation patterns in solid-state and synthetic quantum systems.

Group interactions are widespread, and higher-order extensions of familiar models display collective phenomena absent from their pairwise baselines, which are routinely offered as evidence of a distinct higher-order physics. We ask whether that claim survives the test that gives new a precise meaning in statistical physics: that of universality. Revisiting canonical higher-order models, we argue that what classifies collective behavior are the infrared ingredients that survive at long scales. Arity is not a universality label. Beyond universality, we examine two further questions: whether higher-order structure is a fact about the system or a choice of description, and what data can and cannot tell us about interaction order. Higher-order descriptions remain indispensable when they expose the organizing structure, provide a better mechanistic language, or improve prediction. We close with what should be measured before new phenomena can be claimed, and where higher-order structure already earns its place.

Quantum critical theories with continuously varying critical exponents remain challenging to access experimentally. Here, we propose quantum-dot architectures for realizing the quantum Ashkin-Teller and XYZ/eight-vertex models using resonator-mediated and direct Coulomb interactions, respectively. We focus on the Ashkin-Teller and eight-vertex critical lines, which are connected by a non-local duality relating local order parameters to topological string operators. The resulting platforms provide microscopic control over the parameters explicitly controlling critical exponents, with the Coulomb-based implementation enabling stronger interaction regimes. We demonstrate that continuously varying critical behavior can already be resolved in short chains. For experimentally realistic parameters, the accessible interaction range is expected to produce changes in the critical exponents of order $10\%$ for the Ashkin-Teller model and substantially larger changes for the eight-vertex model. In both cases, the predicted variations lie well above the estimated measurement uncertainties.

Active Brownian particles confined to rigid substrates are known to accumulate near rigid boundaries and, under suitable conditions, undergo motility-induced phase separation (MIPS). A particularly intriguing manifestation of this behavior is the formation of self-sustained vortices around circular obstacles, which act as localized nucleation sites for particle aggregation. While several dynamical properties of such vortices have been previously characterized, their behavior in the thermodynamic limit remains largely unexplored. Here, we investigate how the mass and spatial extent of a dry active-matter vortex scale with system size. Using numerical simulations of repulsive active Brownian Particles interacting with a fixed circular obstacle, we measure the vortex mass, mean radius, and maximum radius as functions of the global area fraction, obstacle size, and system size. We find two distinct scaling regimes. At low densities, the vortex remains localized and its characteristic properties saturate as the system size increases. Above a critical density, however, the vortex mass grows extensively with the total number of particles, while its spatial dimensions scale linearly with the system size, indicating the emergence of an obstacle-stabilized condensed state.

We study semiflexible active polymers confined within a soft, deformable droplet suspended in a fluid using dissipative particle dynamics (DPD) simulations. Extensile and contractile force dipoles along the polymer backbone generate distinct hydrodynamic flow fields that mediate effective interactions between polymer segments. Extensile activity promotes parallel alignment and lateral attraction, whereas contractile activity favors predominantly perpendicular organization, leading to qualitatively different collective behavior. The activity-generated flows couple strongly to the deformable interface, selectively enhancing low-order spherical-harmonic modes associated with long-wavelength droplet deformations. Activity also drives the interfacial relaxation away from passive capillary behavior, with extensile and contractile droplets exhibiting distinct mode-dependent dynamics. These differences in internal organization and interfacial fluctuations strongly influence droplet motility. Extensile activity produces increasingly persistent droplet motion with increasing polymer number, whereas contractile activity can sustain long-lived near-ballistic motion whose duration depends sensitively on polymer number and activity strength. Our results demonstrate how active hydrodynamics governs the interplay between internal structure, interfacial fluctuations, and emergent motility in confined active matter systems.

Detecting phases of matter in general relies on identifying the correct order parameter - a task that remains notoriously difficult for unknown transitions and traditionally is guided by physical intuition and educated guess. Neural networks have recently offered an alternative route by locating phase transitions in known models without any a priori physical knowledge. Yet these approaches remain black boxes and only identify phases without elucidating their properties. Moreover, they often struggle when confronted with realistic, noisy experimental data, which constitute the ultimate testbed for automated methods in physics. Here, we bridge these perspectives by introducing TetrisCNN, a convolutional architecture with parallel branches of differently shaped filters, reminiscent of Tetris blocks, that learns sparse, interpretable latent representations directly in terms of spin correlators. Applied to experimental snapshots of two-dimensional Ising and XY quantum simulators measured in multiple bases, the network not only detects phase transitions and crossovers but also expresses its latent representation and decision boundaries as symbolic formulas built from experimentally measurable spin correlators. This framework opens the way to integrating interpretable neural networks with quantum simulators to uncover and understand new phases of matter.

We study the large-$N$ asymptotics of a constrained spherical integral for non-Hermitian random matrices, in which the norms and mutual scalar product of two vectors are fixed. In the delocalized regime, the asymptotics are governed by the non-Hermitian transforms $\mathcal R_1$ and $\mathcal R_2$. At saddle-point level, the constrained integral exhibits a transition to a localized regime controlled by the largest singular value of a shifted matrix and by the overlap of its associated left and right singular vectors. Motivated by the Hermitian spherical-integral mechanism and by the Coulomb-gas picture, we formulate conjectures for one-eigenvalue large deviations and boundary fluctuations. Throughout, our analysis is carried out in the spirit of mathematical physics.