Glasses and spin glasses, random graphs, neural networks, localization, stochastic processes.
Characterizing crystalline energy landscapes is essential to predicting thermodynamic stability, electronic structure, and functional behavior. While machine learning (ML) enables rapid property predictions, the "black-box" nature of most models limits their utility for generating new scientific insights. Here, we introduce Kolmogorov-Arnold Networks (KANs) as an interpretable framework to bridge this gap. Unlike conventional neural networks with fixed activation functions, KANs employ learnable functions that reveal underlying physical relationships. We developed the Element-Weighted KAN, a composition-only model that achieves state-of-the-art accuracy in predicting formation energy, band gap, and work function across large-scale datasets. Crucially, without any explicit physical constraints, KANs uncover interpretable chemical trends aligned with the periodic table and quantum mechanical principles through embedding analysis, correlation studies, and principal component analysis. These results demonstrate that KANs provide a powerful framework with high predictive performance and scientific interpretability, establishing a new paradigm for transparent, chemistry-based materials informatics.
We present a systematic numerical investigation of the "entanglement geometry gravity" chain in random tensor networks (RTN) established by the ER EPR conjecture and Jacobson's thermodynamic derivation. First, we verify the kinematic foundation: the entanglement first law $δ\langle K\rangle=δS$ (slope=1.000), the encoding of geometry by mutual information (correlation=0.92), and the locality of holographic perturbations (3.3x). We also confirm that gravitational dynamics (JT gravity) does not emerge, identifying a sharp kinematics-dynamics boundary. Second, and more importantly, we discover that many-body localization (MBL) is the mechanism that protects emergent holographic geometry from thermalization. Replacing Haar-random evolution (geometry lifetime $t\sim6$) with an XXZ Hamiltonian plus on-site disorder, we observe a finite-size crossover at disorder strength $W_c\approx10-12$ above which mutual-information-lattice correlations persist indefinitely ($r>0.5$ for $t>50$). We map the full parameter space: the optimal regime is a near-Ising anisotropy $Δ\approx50$ with $W=30$ yielding $r=0.779\pm0.002$ (confirmed by a fine scan over $Δ\in[30,70]$); only holographic (RTN) initial states sustain geometry, while product, Néel, and Bell-pair states do not. MBL preserves the spatial structure of entanglement (adjacent/non-adjacent MI ratio ~2.6-4.2x vs. 1.0x in the thermal phase), rather than its total amount. A comparison with classical cellular automata reveals that MBL uniquely breaks the entanglement-structure trade-off imposed by quantum monogamy: classical systems achieve spatial structure only at the cost of negligible mutual information, while MBL sustains both.
The Bott index has become an indispensable tool to probe the topology of quantum matter, particularly in systems lacking translational symmetry. Constructed from a plaquette operator, it retains the phase information while discarding the amplitude. Here we introduce and develop the Bott metric, which captures this complementary amplitude information and provides a measure of the underlying quantum metric of the system. We show that, in the thermodynamic limit, the Bott metric converges to the trace of the integrated quantum metric. Our framework provides a new route to reveal the quantum metric structure in non-periodic systems, which we illustrate using representative examples ranging from disordered to amorphous models. More broadly, our definition of the Bott metric unifies the notion of topological invariants and quantum metric under the same overarching plaquette operator construction.
We introduce a solvable spherical model of coupled oscillators with fully random interactions and distributed natural frequencies. Using the dynamical mean-field theory, we derive self-consistent equations for the steady-state response and correlation functions. We show that any finite width of the natural-frequency distribution suppresses the finite-temperature spin-glass transition, because the resulting low-frequency singularity of the correlation function is incompatible with the spherical constraint. At zero temperature, however, a spin-glass phase persists for arbitrary frequency dispersion. This residual zero-temperature glassiness is likely a special feature of the spherical dynamics and would be destroyed by local nonlinearities. The model thus provides a solvable oscillator framework for studying how nonequilibrium perturbations suppress finite-temperature glassy freezing.
We analyze the one-pass stochastic gradient descent dynamics of a two-layer neural network with quadratic activations in a teacher--student framework. In the high-dimensional regime, where the input dimension $N$ and the number of samples $M$ diverge at fixed ratio $α= M/N$, and for finite hidden widths $(p,p^*)$ of the student and teacher, respectively, we study the low-dimensional ordinary differential equations that govern the evolution of the student--teacher and student--student overlap matrices. We show that overparameterization ($p>p^*$) only modestly accelerates escape from a plateau of poor generalization by modifying the prefactor of the exponential decay of the loss. We then examine how unconstrained weight norms introduce a continuous rotational symmetry that results in a nontrivial manifold of zero-loss solutions for $p>1$. From this manifold the dynamics consistently selects the closest solution to the random initialization, as enforced by a conserved quantity in the ODEs governing the evolution of the overlaps. Finally, a Hessian analysis of the population-loss landscape confirms that the plateau and the solution manifold correspond to saddles with at least one negative eigenvalue and to marginal minima in the population-loss geometry, respectively.
We investigate the time- and momentum-resolved dynamics of matter waves undergoing elastic scattering from a disordered potential in the presence of spatially uniform SU(2) gauge fields. We derive the disorder-averaged density matrix as a function of time and momentum within the weak-localization regime. By accurately approximating the frequency dependence of the ladder and maximally crossed diagram series beyond the diffusive approximation, we describe short-time spin-momentum dynamics on timescales comparable to the scattering mean free time, for arbitrary strengths of the SU(2) gauge fields and disorder. We also present a cubic equation that determines the spin isotropization time, which gives accurate asymptotic forms in the limits where the spin-orbit length is much longer (Dyakonov-Perel spin relaxation regime) or much shorter than the scattering mean free path, as well as in the SU(2)-symmetric (persistent spin helix) limit. In comparison with numerical calculations, we reproduce both the relaxation of the momentum distribution and the transient backscattering peak with a momentum offset coexisting with the robust coherent backscattering dip, supporting the reliability of our calculations.
We investigate dense higher-order associative memories in the high storage regime when the stored patterns are biased, namely when the entries of the patterns are not symmetrically distributed around zero. In this setting, the standard Hebbian prescription must be modified by recentering and rescaling the pattern entries, and an additional term must be introduced in the Hamiltonian to enforce consistency between the average activity of the network and that of the stored patterns. As a first step, we perform a signal-to-noise analysis in the zero-temperature limit and show that the bias reduces the effective storage capacity through a multiplicative correction factor (1-b^2)^P, while preserving the superlinear scaling with the system size. We then derive the quenched statistical pressure within the Replica Symmetric framework by means of Guerra's interpolation method and obtain the corresponding self consistency equations for the relevant order parameters. The analytical treatment confirms the heuristic prediction of the signal-to-noise argument, showing that the same bias dependent renormalization naturally emerges in the variance of the cross-talk noise. Finally, we discuss the resulting phase behavior of the model and its implications for retrieval performance in the model.
We study Anderson localization in a one-dimensional disordered system with long-range correlated hopping decaying as $1/r^{a}$ with complex hopping amplitudes that break time-reversal symmetry in a tunable fashion by varying their argument. We find analytically a corelation-induced algebraic localization that is robust to a finite strength of the time-reversal-symmetry-breaking parameter, beyond which all states delocalize. This establishes a localization--delocalization transition driven by the interplay between long-ranged correlated hopping and time-reversal symmetry breaking. In addition to obtaining the static localization phase diagram, we also investigate the dynamical phase diagram through the lens of wavepacket spreading. We find that the growth in time of the mean-squared displacement of a wavepacket, which is subdiffusive for the time-reversal symmetric case, becomes diffusive for any finite value of the time-reversal-symmetry-breaking parameter.
We demonstrate that wave amplification enables even weak nonlinearities to reshape linear wave-packet transport in nonreciprocal systems. We study the dynamics of bulk Gaussian wave packets in the Hatano--Nelson model with onsite cubic nonlinearity. We show that the interplay between nonlinearity and amplification generates growing frequency shifts that drive the wave packet through three successive dynamical regimes: an early nonlinear-skin regime with coherent propagation, an intermediate wave-mixing regime driven by mode resonances, and a self-trapping regime in which part of the packet localizes while the remainder ballistically spreads along the system favored direction. The crossover time scales are set by the width and average spacing of the eigen-frequency spectrum. Crucially, within the nonlinear-skin regime, we derive analytical predictions for the wave-packet dynamics and show that nonlinearity couples amplification, dispersion, and nonreciprocity, thereby modifying the magnitude of the wave-packet acceleration and introducing an explicit time dependence into its evolution. Focusing nonlinearities suppress the acceleration and cause it to decrease in time, whereas defocusing nonlinearities enhance it and cause it to increase. We further show that nonlinear interactions typically break down the wave packet before the non-Hermitian jump can occur. Our results provide a route toward accurate control of waves in nonreciprocal metamaterials.
We investigate the localization properties of a generalized SSH model. Numerical and analytical results indicate the emergence of extended states protected by unbounded hopping in this model. Moreover, this protection effect is disrupted by the appearance of generalized incommensurate zeros, causing the extended phase in the system to transition into a multifractal phase. However, at the boundaries of the phase region, we still observe the existence of extended states. These extended states coincide with multifractality-enriched mobility edges, separating the multifractal phase from the localized phase. Further analysis reveals that this extended states originates from the band edge states of SSH model. In addition, these isolated extended states also influence eigenstates with nearby energies, giving rise to an anomalous extended-to-multifractal critical transition. These findings not only enrich the behavioral repertoire of eigenstates at critical points, but also offer new insights for further understanding Anderson localization and the induction of multifractal phases.
We investigate a classical formulation of the Quantum Approximate Optimization Algorithm (QAOA), realized as a Hamiltonian dynamical system of classical kicked tops, which we call the Variational Iterative Rotation Algorithm (VIRAL). The variational parameters are the transverse and longitudinal rotation angles at each of the p layers of the circuit. We find that VIRAL outperforms QAOA on the canonical Sherrington-Kirkpatrick spin-glass benchmark at all circuit depths, with the energy density converging to the ground state value linearly in 1/p. For large circuit depths, the optimized dynamics follows a Floquet protocol in which a pitchfork bifurcation destabilizes the equatorial fixed point and drives the spins toward polar Ising configurations. Our results demonstrate that the effectiveness of QAOA-like protocols derives primarily from their underlying iterative rotation structure, and that a classical implementation of it outperforms its quantum counterpart. We further elucidate its efficiency by reducing the many-body classical evolution to an effective Landau-Lifshitz dynamics for a single spin in a stochastic magnetic field. In this picture, the covariance matrix of the effective field reveals a nearly rank-one structure in which a single mode dominates the stochastic dynamics. In contrast, quantum fluctuations make the noise covariance of the effective quantum model of higher rank, hampering the control of the system. We propose nanometer-scale magnetic tunnel junctions as a natural physical platform for implementing VIRAL, where spin rotations can be realized using magnetic fields and spin torques.
Recently, purpose-built analog hardware that can efficiently minimize the Ising energy and thereby solve a variety of combinatorial optimization problems has been receiving widespread attention. In this work, we show how multidimensional, vectorial degrees of freedom, that are either naturally present or can be artificially created in such hardware, could strengthen the capability to find optimal solutions to optimization problems. In order to achieve this, we introduce a simple model of soft vector spins that should be implementable on a variety of analog hardware platforms as well as three different dimensional annealing methods which harness the enlarged phase space of the vectorial degrees of freedom to minimize the Ising energy. We perform simulations on different benchmark problems and show that for all dimensional annealing methods we tested, vectorial degrees of freedom improve upon one-dimensional degrees of freedom when it comes to finding the ground state of the Ising model. In particular, we find that this advantage becomes most pronounced for $d \gtrsim 3$ dimensional degrees of freedom, with diminishing returns as the dimension is increased further. Our results could inspire new analog optimization hardware and algorithms that explicitly incorporate the advantage of vectorial degrees of freedom.
We introduce a systematic method for expanding general spin-glass Hamiltonians in terms of Mattis interactions, providing a novel perspective for understanding the fundamental differences between short-range Edwards-Anderson (EA) and mean-field Sherrington-Kirkpatrick (SK) spin glasses. By iteratively extracting patterns from the coupling matrix, we expand the original spin-glass system into a Hopfield-like model (a series of Mattis interactions) plus a residual system. Our analysis reveals profound distinctions between EA and SK models: while EA models in two and three dimensions break into isolated subconnected sections after expansion, the SK model exhibits remarkable self-similar behavior, with the residual system preserving the mean-field structure and Gaussian statistics throughout the expansion process. This self-similarity manifests in exponential decay of residual matrix norms and expansion coefficients, reflecting the inherent mean-field nature of the SK model. Furthermore, we demonstrate that pattern expansion can identify ultra-low energy excitations in EA models, revealing excitations with energies that decrease rapidly with expansion step. Through connected component analysis, we quantify the size-energy relationship of these independent excitation clusters, opening new avenues for understanding the low-energy landscape of spin glasses and providing insights into the nature of metastable states.
This study investigates the spatial confinement of topological $π$-modes in one-dimensional chiral-symmetric systems. In conventional periodic and quasiperiodic structures, edge-mode wave functions inevitably penetrate the bulk. To suppress this, inverse design of a potential sequence is performed using a generative model under a global topological constraint. The generated sequence reveals a characteristic structure consisting of a topological boundary layer and a macroscopic S-dense domain, leading to enhanced confinement ($ξ=0.85$) while preserving topology. Based on the physical principle extracted from this result, a minimal heterostructure composed of only two S-blocks is manually constructed, which further reduces the localization length to $ξ=0.75$. These results provide a compact design principle for strongly localized topological states.
We compare the spectrum and the localisation properties of the eigenmodes of the Laplacian and the adjacency matrix of 2D random geometric graphs, using numerical diagonalization of these matrices for different system sizes and connectivities. For sufficiently large ensembles of systems, we evaluate the spectrum, the probability distribution of the participation ratio and the relation between participation ratios and eigenvalues. While all eigenmodes of the adjacency matrix are localised for sufficiently large system sizes, the Laplacian matrix always leads to a small proportion of system-spanning modes due to a conservation law, and therefore to power-law tails in the probability distribution of the participation ratio and its relation to the eigenvalues. By disentangling the effects of finite system size, of mean degree, of component size distribution, and of network motifs, we provide a thorough understanding of the data.
We investigate the parameter space of transformer models trained on protein sequence data using a statistical mechanics framework, sampling the loss landscape at varying temperatures by Langevin dynamics to characterize the low-loss manifold and understand the mechanisms underlying the superior performance of transformers in protein structure prediction. We find that, at variance with feedforward networks, the lack of a first--order--like transition in the loss of the transformer produces a range of intermediate temperatures with good learning properties. We show that the parameters of most layers are highly conserved at these temperatures if the dimension of the embedding is optimal, and we provide an operative way to find this dimension. Finally, we show that the attention matrix is more predictive of the contact maps of the protein at higher temperatures and for higher dimensions of the embedding than those optimal for learning.
We study the number entropy and quasiparticle width in one-dimensional quasiperiodic many-body localized (MBL) systems and observe slow dynamics that have previously been investigated in detail only in random systems. In contrast, quasiperiodic systems exhibit more structured growth of both observables. We identify the modulation of the Rabi oscillation amplitude of single-particle hoppings as the mechanism underlying the slow growth even deep in the MBL regime. This quantum amplitude modulation and associated beats arise from the interaction between single-particle hopping processes at different positions in the chain. Interestingly, this mechanism is not weakened by increasing the distance between particles and is generic to many-body quantum systems. We develop an analytical model based on the aforementioned mechanism that explains the observed dynamics at all accessible timescales and provides a microscopic picture of the slow dynamics in the MBL regime. Our results are consistent with the stability of the MBL phase in the thermodynamic limit.
Randomly-assembled dynamical systems are theoretically predicted to be unstable upon crossing a critical threshold of complexity, as first shown by May. Yet, empirical complex systems exhibit remarkable stability, indicating the presence of additional mechanisms playing a stabilizing role. The relation between complexity and stability is typically assessed by assuming fixed interactions, whereas real systems often evolve in intrinsically time-dependent states. To understand how this affects stability, we linearize a general non-autonomous dynamics around a reference operating state and model the resulting parameters as stochastic processes, which represent the minimal extension of static random interactions to time-varying ones. We derive exact stability bounds that generalize complexity-stability theory to dynamically varying systems. Notably, we find that temporal variability allows systems to remain stable even when their instantaneous Jacobian would predict instability. We compare our results against a non-linear neural network model, where our theory applies exactly, and the generalized Lotka-Volterra equations, where we numerically find that time-varying interactions systematically postpone the onset of replica-symmetry breaking. Overall, our results indicate that temporal variability systematically improves stability, demonstrating a general mechanism by which complex systems can violate classical complexity-stability bounds.
The generation of complex states like multifractal critical states has been an outstanding challenge in both classical and quantum physics. Here we propose a general framework, termed the imaginary gauge phase imprint, allowing to engineer rigorous wavefunctions in any-dimensional non-Hermitian lattices. Using this method, we uncover a novel phase with exact critical wavefunctions in one (and two) dimension, dubbed the skin critical phase (SCP). Unlike conventional critical phases with overall uniform density distributions and non-Hermitian skin effect with eigenstate accumulation at open boundaries, the SCP is marked by a macroscopically multifractal distribution with all critical eigenstates sharing an identical profile and always accumulating at specific bulk interfaces under periodic boundary condition, which become topology-dependent boundary or interface skin modes under open boundary condition. We also show the ballistic dynamics in the SCP, in contrast to the diffusive behaviour in conventional critical phases. Moreover, we validate our method by imprinting configurable wavefunctions in higher dimensions, including complex fractal states with Sierpinski-carpet and Koch-snowflake profiles in non-fractal lattices and Moire states in non-Moire lattices. Our work not only offers fresh insights into fractal phenomena and critical phases, but also provides a rigorous paradigm for wave manipulations in engineered non-Hermitian systems.
The bulk-edge correspondence principle, a cornerstone of topological physics, ensures that first-order topological systems host robust chiral edge states in two dimension. This was later extended to higher-order phases, where second-order topological insulators exhibit localized, topologically protected corner states. While the transition between these distinct phases has been demonstrated in periodic systems, its existence in disordered platforms remains an open question. Here, we demonstrate a controllable topological phase transition between a second-order topological phase and a first-order topological phase in a two-dimensional photonic alloy. By tuning the magnetic doping concentration - implemented by attaching permanent magnets randomly to nonmagnetized yttrium iron garnet rods in an alternately magnetized honeycomb lattice with C3 rotational symmetry - we flexibly control the system's topology. At zero doping, we observe higher-order corner states, confirmed by a trivial Chern number and non-zero bulk polarizations of 1/3. As doping concentration increases, these corner states progressively merge with the bulk states, culminating in the closure of the bulk transmission gap. After the bulk transmission gap reopens with further increased doping, the system transitions to a first-order topological phase, characterized by a nontrivial Chern number of -1 and the emergence of a chiral edge state. This transition is reversible, providing a highly tunable and experimentally simple platform for flexibly switching between localized corner states and delocalized chiral edge states within a single photonic system.