2026-11-06 | | Total: 43
We introduce the notion of an oriented $\vee$-system which gives a class of solutions of the oriented associativity equations generalising Veselov's $\vee$-system solutions of WDVV equations. An oriented $\vee$-system is a special configuration of pairs of a vector and a covector in a linear space. We show that the class of oriented $\vee$-systems is closed under taking subsystems and restrictions. We provide examples which include a multiparameter family related to imprimitive complex reflection groups.
Overcompressive shocks in a $2\times2$ hyperbolic system of conservation laws describing small-amplitude longitudinal-torsional waves in an elastic rod are investigated. A viscous regularization method with an anisotropic dissipative matrix is applied to the system. A qualitative analysis of phase portraits has revealed that a two-parameter family of heteroclinic trajectories, representing the structure of an overcompressive shock, exists in one of six domains in the parameter plane. It is established that the model admits overcompressive shocks, which belong to the first type according to the classification of Liu.
We prove multifractal scaling of measures related to the eigenfunctions of Schrödinger operators with point-like bump potentials, in the thermodynamic limit, in dimension $d\geq 5$. This is consistent with a longstanding expectation in the physics literature in these dimensions. The fractal exponents are computed explicitly as a function of the parameter defining the moments. Our result requires a weak disorder condition on the location of the scatterers and a strong coupling condition on the sequence of eigenvalues.
Tumour evolution and treatment resistance are driven by the interplay between clonal heterogeneity, selective pressures, and therapeutic interventions. In this work, we propose an integrated mathematical and statistical framework to investigate clonal selection, predict radioresistance, and optimize radiotherapy protocols. We develop a multi-clonal dynamical model incorporating proliferation, resource competition, epigenetic instability, and clone-specific radiosensitivity. Through the statistical analysis of 255,150 in silico simulations, we identify key biological determinants of treatment response and characterize the effects of evolutionary competition and radiotherapy on clonal composition. Our results show that evolutionary dynamics favour highly proliferative clones, whereas radiotherapy selects radioresistant populations and reduces tumour heterogeneity. We further investigate how tumour characteristics influence resistance trajectories and identify trade-offs between tumour control, toxicity, and resistance selection across different treatment schedules. Optimal protocols are determined for distinct patient profiles, accounting for dose, fractionation, and timing. Finally, we demonstrate how treatment history affects the efficacy of second-line therapies. Overall, our findings highlight the potential of combining mechanistic modelling and statistical inference to understand cancer evolution and support personalized radiotherapy design.
Topological recursion is a universal procedure that generates an infinite family of multi-differentials with different interpretations in various areas of mathematics. Recently, a new understanding of topological recursion has emerged from $x$-$y$ duality, which relates two such families through explicit formulae. We study $x$-$y$ duality, its parametrised version, and symplectic duality in the setting of equivariant topological recursion. We explicitly analyse the pole structure of the resulting dual families and provide recursive formulae for the different types of polar parts. In contrast to classical topological recursion, the dual family has additional poles at fixed points of the group action, and at anti-diagonals determined by the group action. We show that the partition function $Z_P(E)=\int_{H_N} e^{-N\,\mathrm{Tr}(EΦ^2+P(Φ))}dΦ$ of the Kontsevich matrix model with arbitrary potential $P$, where $dΦ$ is the Lebesgue measure on Hermitian $N\times N$ matrices and $E$ a positive matrix, admits a Fock space realisation related to the character expansion of the IZHC integral. The Fock space framework provides two mutually $X$-$Y$ dual families of multi-differentials relative to an implicitly defined spectral curve $(X,Y)$. The solution $Y(X)$ succeeds for the quartic Kontsevich model with potential $P(Φ)=\fracλ{4}Φ^4$. We identify a network of duality transformations and symplectic shifts which relate the quartic Kontsevich model with the $\mathbb{Z}/2\mathbb{Z}$-equivariant spectral curve. The proposed construction yields a conjectural all-genus solution to the quartic Kontsevich model, which reproduces recursive structures previously established for genus $g\leq 1$.
Non-local couplings offer a means of engineering the dynamics and topology of discrete lattices beyond the constraints of nearest-neighbour (NN) interactions. Here, we extend the mass-spring Su-Schrieffer-Heeger (SSH) model to incorporate beyond-nearest-neighbour (BNN) connections that preserve the chiral symmetry of the diatomic unit cell, investigating how their range and strength shape the bulk response. We show that the accessible winding numbers are bounded by the lattice connectivity, while the spring stiffnesses select the specific topological phase. We further investigate the challenges posed by non-local interfaces, where the extended connectivity inherently introduces local boundary perturbations. The dynamic freedom afforded by BNN couplings is exploited by examining interfaces between dynamically and topologically dissimilar bulk media, revealing a band-gap-induced asymmetry in the localisation of interface states. We show that the interfacial configuration provides a means of controlling the existence and character of localised interface modes, with changes to the interface enabling the tuning of mode frequencies, localisation and hybridisation with the bulk bands.
We describe a symplectic reduction of the Zhukovskiy integrable case of motion of a rigid body about a fixed point with a gyroscope in $\mathbb R^n$ to an integrable magnetic flow on a sphere $S^{n-1}$. We also present Zhukovskiy's geometric interpretation of motion of the Euler top with a gyroscope, aka the Zhukovskiy--Volterra gyroscope, and perform a symplectic reduction of this system to an integrable magnetic flow on a $2$-dimensional sphere.
We investigate the collisional quantum dynamics of an open system in contact with a thermal bath. We propose a particle-bath interaction model based on the Weyl quantization principle and the Wigner phase space description of quantum evolution. By factorizing the reduced Hamiltonian of the open system via supersymmetric quantum mechanics techniques, we construct system-dependent creation and annihilation operators for arbitrary potential landscapes. This allows us to extend the classical Fokker--Planck relaxation mechanism to incorporate Lindblad quantum stochastic jumps within the Wigner phase-space formalism. Through a semiclassical expansion of the Moyal star product to second order, we obtain a generalized Wigner--Lindblad kinetic model featuring position-dependent friction and diffusion terms. We establish the mathematical well-posedness of the Lindblad problem and we validate our approach through numerical simulations. Our model resolves unphysical long-time artifacts of standard quantum Fokker--Planck models in quartic double-well traps, accurately captures dissipative quantum state localization in 2D molecular potentials (4D phase space), and effectively evaluates electron transport and trapping times in InAs/InP resonant tunneling diodes.
We provide an explicit scheme to realize a quantum memory in any 2D non-chiral topological order. For this, we construct syndrome extraction circuits based on the Turaev-Viro path integral as well as a classical algorithm for decoding on the fly. We prove that our scheme has a threshold under a local stochastic circuit-level noise model and discuss how to generalize it to implement topological logical gates.
Sampling experiments are a leading approach to demonstrating quantum computational advantage, yet their robustness to experimental noise is not completely understood. We develop a common framework for classical simulation of boson and fermion sampling under thermal loss, including pure loss. For Fock-state boson sampling, Gaussian boson sampling with fixed uniform squeezing, and fermion sampling with magic inputs, we obtain polynomial-time sampling at inverse-polynomial total variation error, for any fixed nonzero overall loss rate, and with high probability over Haar-random interferometers. A milder flatness assumption on the ideal output distribution suffices for quasipolynomial-time simulation without Haar randomness. With constant thermal loss per layer, the same flatness assumption gives polynomial-time simulation at logarithmic depth. Our bounds for Fock-state boson sampling cover both dilute and saturated regimes and become stronger near full input filling: when only a bounded number of input modes are empty, the noisy output for Haar-random interferometers converges, as the number of modes grows, to a distribution independent of the interferometer.
A signal band-limited to one rectangular band and observed through a finite window is the classical concentration problem of Slepian, Landau, Pollak, and Widom. The eigenvalues of the operator cluster near one and near zero, separated by a transition called the plunge. The plunge width grows like the logarithm of the window length, with a universal coefficient of one over pi squared set by the band-edge jump. This paper follows the clusters and the plunge as the band is deformed. An elementwise power of the kernel smooths the band edge, and the plunge gives way to a power-law edge with exponent one less than the power. Disjoint copies of the band widen the plunge in proportion to the number of band edges, as the multiple-interval theorem of Landau and Widom implies, so the coefficient is a per-edge constant. Shrinking each band to a line gives the classical finite-rank operator. Overlapping bands give a multilevel staircase, and a self-similar band gives a cluster hierarchy governed by the Cantor dimension. At the smooth edge, rescaled eigenvalues do not collapse onto one curve. The edge descriptor passes from a constant to an exponent to an integer to a fractal dimension. For the single band, classical results place all eigenfilter zeros on the unit circle, with more of them inside the band as the order rises, so the dominant eigenvectors act as lowpass filters.
Building on our previous work, in which geometric Gibbs measures were introduced to complete the space of classical Gibbs measures under the Teichmüller metric, we establish the uniqueness of geometric Gibbs measures in this paper. This result completes the construction of a Teichmüller complex manifold structure on the space of all geometric Gibbs measures. Using the welding technique, we uncover a connection between geometric Gibbs measures and harmonic measures on dynamical quasicircles. This connection allows us to establish asymptotic conformal rigidity for these dynamical quasicircles. Finally, we establish a foundation for the systematic study of the ergodic properties of geometric Gibbs measures.
The geometry of spherical spin glasses is relatively well understood, thanks in large part to the Kac-Rice formula. However, this formula does not directly apply to discrete spin systems, and much less is known mathematically about the geometry of Ising spin glasses. In this paper, we present a simple method for counting local optima and gapped states in a broad class of discrete spin glasses. To illustrate its scope, we compute the annealed complexities ("configurational entropies," in the physics language) of mixed Ising spin glasses with potentials, such as external fields, and of the Ghatak-Sherrington model. These computations recover several existing results in the mathematics and physics literature as special cases. We can also compute energy-constrained complexity, finding a critical energy threshold analogous to the spherical case, but for gapped states of pure Ising spin glasses. Finally, we use our method to establish a form of disorder persistence for gapped states, identifying a transition between $3$-spin and $4$-spin glasses.
We prove a logarithmic central limit theorem for a tagged particle in infinite-dimensional Dyson Brownian motion at inverse temperature $β=2$, starting from the sine$_2$ point process. The fluctuation is of order $\sqrt{\log t}$, which is subdiffusive and substantially slower than the polynomial fluctuations of classical short-range interacting particle systems, reflecting the rigidity induced by long-range logarithmic repulsion. We further identify the scaling limit as the singular Gaussian field corresponding to the zero-Hurst parameter, which provides a realisation of the singular Gaussian field as a scaling limit of a tagged particle raised by Dürr, Goldstein and Lebowitz (Comm. Pure Appl. Math. 38 (1985)).
We study the long-time, large-distance asymptotic behaviour of the two-point correlation function of fields of the impenetrable Bose gas in equilibrium states characterized by a quasi-momentum distribution (or filling fraction). This includes, in particular, the cases of Gibbs and generalized Gibbs ensembles. The overall structure of the asymptotic expansion depends on the number of Fermi points, those points for which the quasi-momentum distribution takes the value $1/2$. The leading asymptotic behaviour is determined by a relaxation rate, a logarithmic correction term and an amplitude for which we provide simple and explicit expressions. In the case of two Fermi points the first sub-leading correction is calculated as well.
We consider a prototypical class of dynamics of $N\gg 1$ degrees of freedom subject to two competing terms, one favoring relaxation to a trivial equilibrium and one taking the form of a Gaussian force field promoting the emergence of many local equilibria. The latter is a non-gradient field, i.e. it includes a non-reciprocal component with a tunable strength. For this model, we compute the quenched complexity, i.e. the exponential growth rate of the typical number of equilibria, resolved by stability index. We compare it with the annealed complexity, i.e. the exponential growth rate of the average number of equilibria. We show that the two coincide only in the fully non-reciprocal limit: for any finite degree of reciprocity, typical realizations contain exponentially fewer equilibria than predicted by the average count. Increasing reciprocity enhances the discrepancy between the two. The quenched complexity depends explicitly on the covariance of the random force field and on the degree of non-reciprocity: the universality observed at the annealed level breaks down. The average counting of equilibria correctly identifies the location of the topology trivialization transition, but it does not capture the typical organization of equilibria near this transition, nor the stability index of the dominant equilibria.
We present a systematic framework for gauging continuous symmetries, and show that the Schur-Weyl duality arises naturally as a realization of this construction with open boundary conditions. Restricting to compact Lie group symmetries, we construct dual models for one-dimensional quantum lattice models exhibiting those symmetries. We demonstrate that matrix product operators enable the explicit construction of intertwiners between the dual theories and their (non-invertible) symmetries. With closed boundary conditions, we can relate the full spectrum of these dual theories by considering all possible twisted boundary conditions.We illustrate how the Peter-Weyl transform appears when the duality acts on the boundary Hilbert space of a $G$-symmetric model. The presented framework generalizes previous work on dualities for unitary fusion category symmetric models. It exemplifies how dualities for both finite and continuous symmetries rely on the relation between a bond algebra and a representation of the symmetry in the lattice as double centralizers. We argue that dual non-invertible spontaneous symmetry breaking is the mechanism that underlies symmetric tensor network algorithms. We illustrate these results for $U(1)$- and $SU(2)$-symmetric quantum lattice models.
This is the first of a series of papers which will be devoted to a study of the index of mesoscopic Wilson--Dirac operators on closed Riemannian manifolds. The aim of the series is to discretise the index of the Dirac operator. On flat tori, this is achieved by the Wilson--Dirac operator of lattice gauge theory on a regular lattice. On a general closed Riemannian manifold, however, no regular lattice is available. As an intermediate step between the continuum Dirac operator and its discretisations, we introduce mesoscopic Wilson--Dirac operators on any closed Riemannian manifold, defined by integrating covariant differences over distances comparable to a fixed length scale. In this paper, we define the index of the mesoscopic Wilson--Dirac operators via spectral flow and prove that it agrees with that of the continuum Dirac operator at all sufficiently small scales. The proof rests on an elliptic estimate uniform in the scale and a Rellich-type compactness theorem. In the second paper, we shall pass from mesoscopic operators to discrete Wilson--Dirac operators.
We develop a Pohlmeyer-type reduction for classical non-relativistic strings on string Newton-Cartan backgrounds. For strings on a two-sphere, we identify the invariant reduced field, show that its dynamics is integrable, and construct a Lax representation. We relate this intrinsic formulation to the non-relativistic expansion of relativistic strings and obtain a dictionary between non-relativistic string solutions and their reduced counterparts. The construction explicitly determines which relativistic solitonic configurations survive the non-relativistic limit and which decouple. Moreover, it provides a framework for applying integrability methods to Newton-Cartan string dynamics.
Thermodynamic curvature is widely used as a diagnostic of interactions and phase transitions, yet its interpretation is usually tested only in the thermodynamic limit. Small spin clusters are physically relevant in their own right, since a macroscopic sample of independent molecular trimers or a few-qubit quantum simulator is governed exactly by the few-spin partition function. They also provide the only setting in which the geometry can be obtained in closed form together with its full dependence on the interaction topology. We study the Fisher-Rao (Ruppeiner) geometry of three Ising spins with arbitrary pairwise couplings in a longitudinal field. This cluster interpolates continuously between an open chain and a geometrically frustrated triangle. For the symmetric triangle we obtain the scalar curvature in closed form. It vanishes on the noninteracting line and on a line of fixed field-to-temperature ratio that does not depend on the coupling, which coincides with the inflection point of the single-spin susceptibility. The near-zero behaviour is therefore exactly linear and reflects analyticity rather than critical scaling. In the frustrated antiferromagnetic regime the curvature remains finite and negative down to zero temperature. A general low-temperature formula relates this plateau to the change in magnetization content between the degenerate ground manifold and its first excitation. Any bond anisotropy destroys the plateau below a crossover temperature that depends only logarithmically on the anisotropy, so the curvature is an extremely sensitive detector of exact ground-state degeneracy. On the ferromagnetic side the zero-temperature curvature is universal for all couplings, whereas additional zero lines in anisotropic clusters track field-induced crossings of excited levels.
We prove exponentially slow thermalization for the quantum random energy model (QREM) under Davies dynamics. The QREM is the $n$-qubit Hamiltonian $H=ΓT+U$, where $U$ assigns independent Gaussian energies with variance equal to the system size to computational-basis configurations and $T$ is the transverse field generated by single-qubit Pauli $X$ operators. Let $β_c=\sqrt{2\log 2}$ be the freezing inverse temperature of the classical random energy model. For every inverse temperature $β>β_c$ and every transverse-field strength $0\leqΓ<Γ_c(β)$, where $e^{ββ_c}=2\cosh(βΓ_c(β))$, we show that the mixing time of the dissipative Davies dynamics associated with single-qubit Pauli system--bath couplings grows exponentially with the system size, almost surely; the bound is unchanged if the standard zero-Bohr-frequency dephasing terms are included. The result is uniform over all Davies rates satisfying detailed balance and a polynomial upper bound and does not assume primitivity. In the spin-glass phase, the Gibbs weight is carried by a sparse family of deep classical configurations that are separated by an extensive Hamming distance. A localization theorem for the QREM [Manai, Warzel, Comm. Math. Phys. 402, 1259 (2023)] transfers this geometry to the corresponding low-energy quantum eigenstates. Consequently, the Gibbs weight of the intermediate spectral sector is exponentially small, while single-qubit bath operators have exponentially small matrix elements between the two macroscopic sides of the bottleneck. A Davies-specific flux estimate then converts these two facts into an exponential lower bound on the mixing time. This gives a rigorous open-system dynamical manifestation of the frozen, localized low-energy structure of the QREM throughout its equilibrium spin-glass phase below the quantum transition line.
In this paper, we construct the complex Ginzburg-Landau dynamic associated with the harmonic oscillator, with a random force which is white in time and covers the full subcritical range in space. We are able to renormalise the non-linearity using the same constant as the one used to renormalise the problem associated to the Laplacian. By an extension of the problem to complex times, we are able to prove that the solution of the complex Ginzburg-Landau is equal in law to the solution of the purely dissipative case, i.e. the associated complex $Φ^4$ model, but along a complex axis of time.
The square-lattice Ising model is exactly solvable in zero magnetic field, but the scope of exact solutions at nonzero field remains unresolved. We show that, at every fixed dimensionless field $H=h/(k_BT)\ne 0$, its free energy and magnetization obey no algebraic differential equation in the coupling $K=J/(k_BT)$. This excludes Onsager-type formulas and broad classes of special-function expressions. Previous work left open a countable set of fields with algebraic-integer fugacities. We eliminate these exceptions using Galois conjugation and the large-order growth of droplets of reversed spins. At conjugate fugacities outside the unit circle, nearly square droplets produce coefficient growth incompatible with an algebraic differential equation. We prove that cancellations cannot remove this growth in the full pressure series. The free-energy result also holds for the Lieb lattice and square lattices with equal-length chain decorations and uniform coupling and field. In every dimension $d\ge2$, we additionally resolve fields for which $e^{-2|H|}$ has a real algebraic conjugate greater than one.
Ostracism was the practice in use in Athens and other ancient Greek cities to send into exile citizens whose influence was deemed to be dangerous. Quite similarly, splay ostracism is here the approximation that bans splay distortion from thin films of ferroelectric nematic liquid crystals on a curved shell in space. This reduces to a mere geometric problem the quest for equilibria of director vortexes and the walls isolating them, as splay ostracism requires the director field lines to be everywhere orthogonal to systems of geodesics on the shell. Consequences of this theorem are explored for shells with Gaussian curvature of either sign. The resulting quasi static dynamics reveals both cases where walls escape from the shell and cases where they are trapped on it, with possibly different reflections on the vortex dynamics.
We study the Boltzmann equation for a gas on the three-dimensional torus. Its particles collide with one another and with a thermal reservoir whose temperature varies in space. The reservoir heats the gas in some regions and cools it in others, driving a nonequilibrium energy current. For the isotropic cutoff kernel, with Maxwell molecules or hard potentials, we prove global well-posedness and exponential convergence to the same stationary state for initial data with a common Gaussian upper bound and Gaussian moment bound. These bounds determine how weak self collisions must be relative to bath collisions. The estimates are uniform over bath profiles between fixed positive temperature bounds. For any nonconstant reservoir profile on the three-dimensional torus, this stationary state cannot be a local Maxwellian. When particles collide only with the bath, relative entropy to the stationary state cannot increase. Once self collisions are added, arbitrarily small perturbations can initially increase it, even though they later return exponentially.