2026-09-25 | | Total: 37
The Anderson tight-binding model is a fundamental model of quantum transport and localization in disordered media. Completing a problem left open by Bourgain and Kenig, this paper proves Anderson localization near the bottom of the spectrum for the lattice Anderson model with Bernoulli potential, in any dimension $d\ge 2$. The proof uses the multiscale framework of Fröhlich-Spencer and Bourgain-Kenig, and the main new ingredient is a probabilistic discrete unique continuation principle (PDUC) for the discrete Schrödinger equation. This PDUC is established via a bootstrap argument and a key probabilistic lemma proved by adaptively revealing the random potential.
This paper systematically investigates admissible configurations of the interaction between a rarefaction wave (RW) and a dispersive shock wave (DSW) by constructing a two-step piecewise constant initial condition in the Gerdjikov--Ivanov (GI) equation with non-convex hydrodynamics. Owing to the two-valued mapping associated with the Riemann invariants, the GI equation also supports contact DSWs (CDSWs), which are included in the interaction analysis. Using Whitham modulation theory, we elucidate the self-similar descriptions of the incoming waves, derive the boundary matching conditions, and obtain analytical expressions of the full evolution before, during, and after the interaction through generalized hodograph transformation. The results show that the interaction of a RW with either a DSW or a CDSW can exhibit two distinct scenarios: in one case, the two waves exchange spectral parameters after the interaction and fully separate, propagating independently with modified phases; in the other case, the wave structures merge upon interaction and subsequently evolve into either a soliton train or a small-amplitude harmonic wave for a classical DSW, and into an algebraic soliton train for a CDSW. The theoretical results are in good agreement with numerical simulations. This work can provide theoretical support for wave interaction phenomena in non-convex dispersive systems.
We consider a hierarchy of models that can be obtained from the original Nyzhnyk system by imposing conditions on parameters, introducing potentials or pseudopotentials, performing limiting processes with respect to a scaling parameter, applying differential substitutions and interpreting parts of independent and/or dependent variables as complex or real. Therefore, among the Nyzhnyk models, one can distinguish between symmetric and asymmetric, dispersive and dispersionless, standard and modified, as well as real, complex, mixed and specific models. One can also distinguish single partial differential equations or systems of such equations, as well as linear or nonlinear Lax representations in the dispersive or dispersionless cases, respectively. For each specified model, we compute the maximal Lie invariance pseudoalgebra and, in the symmetric case, find the point- and contact-symmetry pseudogroups using the megaideal-based version of the algebraic method. It is shown that relations between these pseudoalgebras and between these pseudogroups are induced by relations between the corresponding models. Defining (finite-dimensional) subalgebras are singled out in all these pseudoalgebras in the symmetric and specific cases. Based on the established correspondences between the dispersionless and dispersive models, we completely classify one- and two-dimensional subalgebras of the maximal Lie invariance pseudoalgebra of the (dispersive symmetric potential) Nyzhnyk equation and one-dimensional subalgebras of the maximal Lie invariance pseudoalgebra of its linear Lax representation.
The winding number is a simple topological invariant. In the case of chiral symmetry it characterises gapped phases of Fermions. We study statistical properties of this topological index or invariant in a chiral symplectic setting using Random Matrix Theory. We consider ensembles of Hamilton matrices in the symmetry class CII (quaternionic matrices) according to the classification in the tenfold way. We set up a parametric random matrix model and derive expressions for parametric correlations of the winding number density as well as for the discrete winding number distribution. We found a super-universality in the limit of large matrix dimensions for the one- and two-point correlators for the bulk of parameters, meaning that the results agree up to rescaling with those of the class AIII (complex matrices). In this context we employ a new method of unfolding and discover the Gaussian behaviour of the winding number distribution.
Wetting transitions on textured substrates play a central role in the design of functional surfaces, but resolving their transition mechanisms requires efficient exploration of complex energy landscapes. In particular, saddle points are essential for revealing the connectivity among metastable states and transition pathways. In this work, we develop a diffuse-domain accelerated saddle dynamics (DD-ASD) framework for mass-constrained wetting transitions on complex textured substrates. The diffuse-domain formulation embeds complex solid geometries into a Cartesian grid, avoiding body-fitted mesh construction in repeated saddle-point searches. An orthogonal projection is applied consistently to the phase-field state and unstable-direction dynamics, preserving the prescribed droplet mass at the discrete level. Momentum acceleration is incorporated into the saddle dynamics to improve the efficiency of high-index saddle searches. Under the stated local spectral and exact-eigenspace assumptions, we establish a local convergence rate of $1-\mathcal{O}(1/\sqrtκ)$ on the effective mass-conserving subspace. Numerical experiments demonstrate the effectiveness of the proposed method in resolving wetting solution landscapes, transition pathways, and energy barriers on textured substrates. We further extend the framework to fully three-dimensional wetting-landscape computations.
Equivalence groups associated with the field equations of arbitrary motions of inhomogeneous hyperelastic solids are investigated. The explicit forms of the infinitesimal generators, namely the components of the isovector field associated with the equivalence groups are derived. Additional constraints on the isovector components are established to ensure that the transformed equations are supported by a strain energy function. Possible mappings between the equations of materials with different properties are examined. Not only the transformations between homogeneous and inhomogeneous hyperelasticity equations, but also the transformations between the equations and strain energy functions of isotropic and anisotropic materials are investigated. The results are illustrated with several examples.
The random hopping Hamiltonian is a toy model for a disorder-driven topological phase transitions in one-dimensional chiral Hamiltonians. It has a vanishing Lyapunov exponent and a Dyson peak in the density of states at zero energy. This work shows how the quantum dynamics leads to a logarithmic growth of the moments of the position operator, similar as in classical Sinai diffusion. The mechanism at the origin of this phenomenon is that the eigenfunction of the smallest eigenvalue is well-approximated by the exponential of a classical random walk, leading to two essentially independent subexponential localization centers.
We introduce a bi-oriented interlaced matrix product on $2\times 2$ matrices over a unital, possibly non-associative algebra, extending a matrix mnemonic for Cayley--Dickson doubling. The opposite orientations $X_R$ and $X_L$ coincide precisely when the underlying algebra is commutative. For the canonical position and momentum operators $Q,P$ on $L^2(\mathbb{R})$, a matrix $D$ and its companion $\widetilde{D}$ satisfy, on the common Schwartz domain, $\widetilde{D} X_R D =\widetilde{D} X_L D =2H Id_2$, where $H=(P^2+Q^2)/2$. We study the associated interlaced eigenvalue equations for operator-valued columns. The combinations $Φ=A+iB$ and $Ψ=A-iB$ decouple the equations into two-sided Sylvester equations. We construct explicit rank-one generalized eigenoperators from position and momentum eigenstates. Their distributional kernels belong to the doubled rigged Hilbert space $S(\mathbb{R}^2)\subset L^2(\mathbb{R}^2) \subset S'(\mathbb{R}^2)$; the corresponding position and momentum operators act on different tensor factors and therefore commute. The rank-one families do not exhaust the distributional solutions: higher-order generalized solutions are also exhibited. We distinguish the algebraic identities on the Schwartz core from questions concerning closed-operator realizations.
For the neutral periodic one-dimensional Riesz gas with pair potential locally $-|x|^a$, $0<a<1$, we prove a particle-displacement variance bound of order $β^{-1}$, uniformly in the number of particles. Log-concavity also gives exponential displacement tails. Every stationary periodic thermodynamic limit is simple, has intensity one, and admits a stationary ordered matching to the unit lattice with the same bounds. Each limit satisfies the canonical Gibbs equations for the full-line interaction, with an ordinary symmetric spatial principal value for the exterior potential. The matching implies uniformly bounded interval number variance and a positive limiting second moment of the reciprocal-lattice Fourier average at sufficiently low temperature. The main estimate compares the inverse random Hessian, through deterministic electrical flows, to a transient long-range network.
We describe a family of convex uniformly continuous functions $E^λ_F$, $λ>0$, on the set of states of a bipartite quantum system (consisting of finite-dimensional or infinite-dimensional subsystems), which monotonically increase and converge pointwise to the Entanglement of Formation (EoF) as $λ\to0$. These functions are "nonselective" entanglement monotones defined by the way close to the construction of the Moreau-Yosida regularization (the Moreau envelope) of a convex function on a convex set used in the modern convex analysis. So, we call the functions $E^λ_F$ the Moreau-Yosida approximations of the EoF and describe their equivalent definitions and basic properties. The semicontinuity bounds for the EoF (presented in [Lob.J.Math., 46(6), 2632-2658]) allow us to obtain easily computable upper bounds on the difference $E_F(ρ)-E^λ_F(ρ)$ for a given state $ρ$. These bounds give easily computable bounds on the rate of uniform convergence of the function $E^λ_F$ to the EoF as $λ\to0^+$ on the sets of states with bounded rank/energy of one of the marginal states. We also discuss sufficient conditions for the coincidence of $E_F(ρ)$ and $E^λ_F(ρ)$ at a given state $ρ$ for all $λ$ small enough and consider several classes of states for which such coincidence takes place. The conjectured selective LOCC-monotonicity of the functions $E^λ_F$ and a possible way to prove it are briefly discussed. Secondary parts of the article are written with the help of ChatGPT-5.6.
Non-Abelian quantum codes connect quantum error correction, phases of matter, and computational resources. In this work, we develop a general framework for constructing non-Abelian quantum low-density parity-check (qLDPC) codes by gauging sheaf codes via cup products and use it to obtain families with constant encoding rate and linear distance. We resolve the coupled logical constraints using explicit representatives to characterize the full gauged code space. We provide a fundamental treatment of code distance based on the general Knill--Laflamme condition and combine expansion with cleaning to establish protection against arbitrary low-weight errors. We further construct an almost-good family whose entire code space exhibits long-range magic. Gauging and ungauging also enable logical Clifford measurements that prepare encoded magic states. These results extend good qLDPC codes beyond the Pauli stabilizer setting and provide a concrete foundation for exploring non-Abelian phases beyond geometric locality and pursuing the no low-energy trivial magic conjecture.
We determine the exact universal high-temperature separability threshold for Pauli Hamiltonians of bounded degree $Δ\ge2$. If every coefficient in the Hamiltonian has magnitude at most one and each term has overlapping support with at most $Δ$ other terms, the Gibbs state is a mixture of product Pauli eigenstates whenever \[ β\le z_Δ:= \operatorname{arctanh}\left[\max_{0\le x \le 1}x\left(\frac{1-x}{1+x}\right)^{Δ-1} \right]. \] For every $β>z_Δ$, a finite commuting Hamiltonian with maximum overlap degree at most $Δ$ has an entangled Gibbs state. At any fixed $β<z_Δ$ strictly below the threshold, a classical polynomial-time algorithm produces samples from a distribution over product Pauli eigenstates approximating the Gibbs state in trace distance.
We explicitly construct the additional integrals of motion, ensuring maximal superintegrability of the discrete-time rational Ruijsenaars-Schneider model. Using them, we investigate the algebraic aspects of superintegrability in both continuous- and discrete-time settings. In particular, we determine the complete structures of the polynomial symmetry algebras associated with both the rational Ruijsenaars-Schneider model and its discretization. We demonstrate that discretization leads to a nontrivial deformation of the continuous symmetry algebra with respect to the discretization parameter, thereby extending recent analogous results from the rational Calogero-Moser system to its relativistic generalization.
The classical Becker-Döring equations describe the formation of clusters by aggregation and fragmentation of monomers. If the total amount of mass is supercritical, larger and larger clusters are formed, leading to an asymptotic loss of mass in the long time limit. The two-component Becker-Döring system arises when clusters are built from two types of monomers. Here, we study the natural extension of the one-component system, where no energy or entropy is leaving or entering the system - the so called detailed balance assumption. We rigorously prove that all initial conditions admit a solution minimising the relative entropy as time approaches infinity. This shows, that the long time limit in the weak* topology is selected through the initial Type I and Type II masses. Furthermore, the proportion of lost Type I mass is determined by the limit point via the mixing ratio of Type I and Type II monomers that maximises the binding energy. The main difficulty of the two-component system is that no relative entropy is weak* continuous, which is crucial for the classical argument [2]. Instead, our proof is based on a discrete two-dimensional logarithmic Sobolev inequality, that bounds the entropy dissipation on the correct timescale. Our approach also improves current assumptions to determine the long time behaviour for the one-component system.
The two-component Becker-Döring model is a two-dimensional coagulation-fragmentation equation that describes the dynamics of clusters built from two different types of monomers growing and shrinking through monomer interaction. We show the existence and mass conservation of solutions and give a partial result regarding uniqueness. Furthermore, we prove under the assumption of detailed balance the entropy equation, followed by a detailed study of the equilibrium states and relative entropies. We can determine the minimum of all relative entropies under sequences of fixed Type I and Type II masses and fully describe the minimising sequences in the weak* topology.
For independent bond percolation on $\mathbb{Z}^d$ with parameter $p$, let $p_c^b(d)$ be the critical probability. We prove that for each $d \geq 9$, there is $ε_d>0$ such that for each $p \in (p_c^b(d), p_c^b(d)+ε_d)$, the complement of the infinite open cluster stochastically dominates a supercritical site percolation on $\mathbb{Z}^d$. This improves the previous results by Grimmett, Holroyd and Kozma 2014, and Bock, Damron, Newman and Sidoravicius 2020. Numerical estimates of $p_c^b(d)$ and $p_c^s(d)$ (the site critical probability) suggest that a similar stochastic domination result holds for all $d \geq 4$. For the Ising model on $\mathbb{Z}^d$ with inverse temperature $β$, let $β_c(d)$ be the critical inverse temperature. We prove that for each $d \geq 8$, there is $ε_d>0$ such that for each $β\in [0,β_c(d)+ε_d)$, the $+$ spins under the minus phase stochastically dominate a supercritical site percolation on $\mathbb{Z}^d$. This improves the previous result of Aizenman, Bricmont and Lebowitz 1987. Numerical estimates of $β_c(d)$ and $p_c^s(d)$ suggest that a similar stochastic domination result holds for all $d \geq 5$.
We study the $L^p$-boundedness of Isozaki--Kitada type modified wave operators associated with Schrödinger operators $P=-\partial_x^2 + V$ in one dimension for potentials $V$ including long-range potentials such as negative Coulomb-like ones $V(x) = -(1+|x|^2)^{-μ/2}$ for $μ>0$. For such potentials, if $ μ\in (1, 2)$, we prove that the middle and the high energy parts are bounded on $L^p (\mathbb{R})$ for any $p \in [1, \infty]$ but the low energy part is unbounded except for $p=2$. Moreover, if $μ\in (0, 1]$, we show that even the middle energy part is unbounded except for $p=2$. Actually, regarding the $L^p$-boundedness in the middle and the high energy regimes, we give a complete classification of slowly decaying potentials without imposing any sign condition.
We study the phase boundary of the finite-lifespan frog model. For rate-one continuous-time simple random walk on an infinite, connected, locally finite transitive graph of superlinear growth, we prove that the critical-density curve is continuous and strictly decreasing, with $-\logλ_c$ locally bi-Lipschitz. The inverse critical-lifespan curve has the corresponding regularity wherever it is finite. Whenever every positive density has finite critical lifespan, this resolves a conjecture of Angel, de la Riva, Hermon, and Shi on the regularity of the critical-parameter curves; in particular, it does so on nonamenable and superlinear polynomial-growth graphs. We also establish a small-lifespan scaling limit and extend critical-curve regularity and sharpness to a class of long-range frog models. The main tool is a local-to-global principle for activation processes generated by independent finite rooted ranges: local one-hit comparison implies comparison of global reachability and survival.
Let $\mathbf{c}=(λ,α,μ,ν)\in\mathbb{R}^{4}$ denote the four running couplings of four-dimensional nonplanar Euclidean scalar $\varphi^{4}$ theory, and let $B(\mathbf{c})=\sum_{p\ge2}\sum_{|\mathbf{m}|=p}β(\mathbf{m})\mathbf{c}^{\mathbf{m}}$ be its formal beta-function series. Assuming the factorial coefficient bounds $|β(\mathbf{m})|\le (p-1)!C^{p-1}$ for $|\mathbf{m}|=p$, we construct an explicit summation rule $\mathcal{S}$, given by a cut Borel--Laplace transform whose cut is placed by a scale-covariant real-analytic gauge. We prove that it assigns a well-defined real value to $B$ on an explicit nonempty open domain, in fact on all of $\mathbb{R}^{4}$. The summed function $\mathcal{S}B$ is linear and real, real-analytic on $\mathbb{R}^{4}\setminus\{0\}$, $C^\infty$ on $\mathbb{R}^{4}$ with Taylor series $B$ at the origin, and Gevrey-$1$ asymptotic to $B$ to all orders with an explicit remainder bound. At optimal truncation the remainder is $O(e^{-η/(C\|\mathbf{c}\|_\infty)})$ for any $η<1$. We show that $\mathcal{S}B$ reproduces the classical Borel sum and the ordinary sum of a convergent series up to errors of the same exponentially small order. Exact reproduction is impossible for any linear rule defined on the whole admissible class: regularity forces all weights to equal $1$, leading to divergence. Finally, we prove optimality under the stated hypothesis. Admissible series may diverge at every $\mathbf{c}\ne0$, and their Borel transforms may have the circle $|τ|=1/(C\|\mathbf{c}\|_\infty)$ as a natural boundary, so classical Borel-type summation need not be available on the class. Moreover, there exist two functions indistinguishable at the level of the coefficient bounds yet differing by exactly $(2π/C)e^{-1/(C\|\mathbf{c}\|_\infty)}$. All constants are explicit and all arguments are self-contained.
We study graded arrays of highly dispersive subwavelength resonators through their volume-potential matrix reduction. The reduction gives a finite-dimensional nonlinear resonance problem whose diagonal entries encode the local dispersive response and whose off-diagonal entries encode volume-mediated coupling between distinct resonators. We show that monotone material grading turns this matrix into an analytically controllable frequency sorter. When graded detuning dominates coupling in the reduced characteristic matrix, the characteristic vectors localize near individual resonators; under an explicit half-lattice-site localization condition, this localization becomes a monotone frequency-position map. We also prove a constructive inverse-design statement: prescribed monotone target frequencies can be realized by choosing the local dispersive parameters, with robustness under sufficiently small grading errors. Numerical experiments, whose parameters and algorithms are recorded in an appendix, illustrate localized resonances, designed frequency-position maps, disorder stability, spacing dependence, and a resolved physical Lippmann-Schwinger validation beyond the leading logarithmic quadrature model.
Quantum incommensurate systems have attracted widespread interest due to their unique physical properties. Related studies have made notable progress in recent years. To gain deeper insight, it is both significant and challenging to study a broader range of physical observables for these systems, which requires a comprehensive understanding of their spectral properties and wavefunction behavior. Based on the regularized model recently proposed, this work introduces a regularization framework, rendering physical observables for incommensurate systems mathematically well-defined and computationally accessible with theoretical guarantees. Based on this framework, taking the density of states and the electron density as representative examples, we establish their mathematical characterization, derive their planewave approximations, and provide a rigorous convergence analysis. Numerical experiments validate the effectiveness of our method, demonstrating its practical applicability for generic 1D and 2D incommensurate systems.
A quantum measurement transforms a quantum statistical experiment into a classical one, but different measurements generally yield experiments with different statistical information. We introduce quantum-to-classical sufficiency: relative to a prescribed measurement class, a measurement is sufficient if its induced classical experiment Blackwell-dominates those from every other measurement in the class. This separates two components of quantum inference: reducing the measurement class through symmetry or decision-theoretic arguments, and identifying a canonical measurement in the reduced class by an experiment-level sufficiency condition. We establish a factorization criterion where the parameter enters the group-averaged experiment only through weights of an orthogonal block decomposition. The block-label measurement is then sufficient relative to all invariant measurements. For tensor-product models under unitary conjugation, this identifies weak Schur sampling as the sufficient measurement for spectral inference. Combined with risk-preserving symmetrization, it reduces invariant Bayes and minimax problems over all measurements exactly to classical decision problems based on Young diagrams. We derive sharp first-order asymptotic Bayes and minimax risks for smooth spectral functionals. We establish an analogous reduction for a noncompact group action in a continuous-variable setting. We consider displaced thermal states and functionals of the thermal parameter with an unknown displacement. Within displacement-invariant measurements, total residual photon number measurement emerges as the sufficient measurement. Combined with Hunt-Stein reduction, this yields exact finite-sample reductions of Bayes-minimax and minimax problems over all measurements to classical decision problems based on the total residual photon count, which follows a negative binomial distribution.
We establish a pointwise lower bound for the heat kernel of the elliptic operator $ Λ=(1+|x|^α)Δ+b|x|^{α-2}x\cdot\nabla-|x|^β$, where $d\geq3$, $α>2$, $β>α-2$, and $b\in\mathbb R$. For every $τ>0$, we prove that $$ \begin{aligned} p(t,x,y) &\geq C_τe^{λ_0t} \left(\frac{1+|y|^α}{1+|x|^α}\right)^{\frac{b}{2α}} \frac{(|x||y|)^{-\frac{d-1}{2}-\frac{β-α}{4}}}{1+|y|^α}\\ &\quad\times\exp \left[ -\int_1^{|x|}\sqrt{\frac{s^β}{1+s^α}}\,\mathrm ds -\int_1^{|y|}\sqrt{\frac{s^β}{1+s^α}}\,\mathrm ds \right] \end{aligned} $$ for all $t\geqτ$ and $|x|,|y|\geq1$, where $λ_0<0$ is the largest eigenvalue of $Λ$ and $C_τ>0$ is independent of $t,x,y$. The proof applies the classical Davies - Simon argument in a weighted symmetric setting.
We study large-overtone quasi-normal modes of a massive scalar field in a four-dimensional Schwarzschild de-Sitter black hole in the presence of a wall placed at the static sphere. Dirichlet boundary condition is imposed at the wall. At leading order the quasi-normal modes with large overtones scale as $n$ for $n\gg 1$ and are inversely proportional to the bouncing time. We analyze correction terms up to $n^{-1}$ order in the asymptotic expansion.
This data article provides the datasets, symbolic derivations, and scripts used to reproduce master diagrams, stationary-phase spectra, windowed first-order coherence functions, and quantum-noise stability maps for strongly chirped dissipative solitons of the cubic-quintic complex Ginzburg-Landau equation in normal and anomalous group-delay dispersion regimes. The repository includes node-regularized normal-dispersion spectra and energies; small-parameter expansions of the branch roots; cavity-map gain-loss update relations; Airy uniformization at the normal-dispersion spectral edge; anomalous-dispersion spectra and coherence calculations; and processed tables for plotting and stability analysis. OriginLab projects are accompanied by open-format .csv/.txt numerical tables to support reuse without proprietary plotting software. Data and code repository: https://doi.org/10.5281/zenodo.22690899.