2026-02-10 | | Total: 26
We propose a geometric model predictive control framework for quantum systems subject to smoothness and state constraints. By formulating quantum state evolution intrinsically on the projective Hilbert space, we penalize covariant accelerations to generate smooth trajectories in the form of Riemannian cubics, while incorporating state-dependent constraints through potential functions. A structure-preserving variational discretization enables receding-horizon implementation, and a Lyapunov-type stability result is established for the closed-loop system. The approach is illustrated on the Bloch sphere for a two-level quantum system, providing a viable pathway toward predictive feedback control of constrained quantum dynamics.
We propose a metric-dependent geometric variant of factorization homology in conformally flat Riemannian geometry for $d \geq 2$. We introduce a symmetric monoidal category of germs of d-dimensional Riemannian manifolds and orientation-preserving conformal open embeddings, and its full monoidal subcategory generated by flat disks. A conformally flat $d$-disk algebra is a symmetric monoidal functor from this disk category to a target category; in this paper we take the target to be $\mathrm{IndHilb}$, the ind-category of Hilbert spaces, which provides a mathematical formulation of $d$-dimensional conformal field theories. The (1-categorical) left Kan extension of an $\mathrm{IndHilb}$-valued conformally flat $d$-disk algebra defines a metric-dependent invariant of conformally flat manifolds. Under suitable positivity and continuity assumptions, we prove that its value on the standard sphere $(S^d,g_{\mathrm{std}})$ reproduces the sphere partition function of the associated conformal field theory. For $d>2$, we construct nontrivial $\mathrm{IndHilb}$-valued conformally flat $d$-disk algebras from unitary representations of $\mathrm{SO}^+(d,1)$.
We determine the isometry group of the $n$-qubit state space with respect to the quantum Wasserstein distance induced by the so-called symmetric transport cost for all $n \in \mathbb{N}.$ It turns out that the isometries are precisely the Wigner symmetries, that is, the unitary or anti-unitary conjugations.
Many novel and unique physical phenomena of incommensurate systems can be illustrated and predicted using the spectra of the associated Schrödinger operators. However, the absence of periodicity in these systems poses significant challenges for obtaining the spectral information. In this paper, by embedding the system into higher dimensions together with introducing a regularization technique, we prove that the spectrum of the Schrödinger operator for the incommensurate system can be approximated by the spectra of a family of regularized Schrödinger operators, which are elliptic, retain periodicity, and enjoy favorable analytic and spectral properties. We also show the existence of Bloch-type solutions to the Schrödinger equation for the incommensurate system, which can be well approximated by the Bloch solutions to the equations associated with the regularized operators. Our analysis provides a theoretical support for understanding and computing the incommensurate systems.
We study Schrödinger operators on quantum graphs where the number of edges between points is determined by orbits of a "shift of finite type". We prove Anderson localization for these systems.
This paper concerns the derivation and validity of macroscopic descriptions of wave packets supported in the vicinity of degenerate points $(K,E)$ in the dispersion relation of tight-binding models accounting for macroscopic variations. We show that such wave packets are well approximated over long times by macroscopic models with varying orders of accuracy. Our main applications are in the analysis of single- and multilayer graphene tight-binding Hamiltonians modeling macroscopic variations such as those generated by shear or twist. Numerical simulations illustrate the theoretical findings.
Bose gas in a large cubic box with periodic boundary conditions interacting with a small potential with a positive Fourier transform. We compute the imaginary part of the phononic excitation spectrum in the lowest order of perturbation theory in thermodynamic limit at low temperatures and low momentum. Our analysis is based on perturbation theory of the standard Liouvillean. We use two approaches: the first, motivated by the standard representation of operator algebras, examines resonances near zero; the second analyzes the 2-point correlation function in the energy-momentum space.
We establish the existence of spinning $Q$-vortex solitons in a complex scalar field theory with a sextic potential on a finite domain. By reducing the governing equation to a nonlinear boundary value problem, we use variational methods to prove the existence of at least two distinct types of solutions: a ground state solution obtained via constrained minimization and an excited state of the saddle-point type obtained via the Mountain Pass Theorem. We derive bounds for the angular frequency $ω$, the wave amplitude, and the domain size $P$, and provide explicit estimates for the exponential decay of the solutions. Furthermore, we implement a spectral-Galerkin formulation to numerically compute the profiles of fundamental $Q$-vortices, illustrating the saturation behavior of the soliton's amplitude and the asymptotic dependence of the frequency on a prescribed reduced norm and vortex winding number, as well as verifying the theoretical results and visualizing the topological phase structure of the solutions.
Motivated by earlier numerical evidence for a percolation-like transition in space-time jamming, we present an analytic description of the transient dynamics of the deterministic traffic model elementary cellular automaton rule 184 (ECA184). By exploiting the deterministic structure of the dynamics, we reformulate the problem in terms of a height function constructed directly from the initial condition, and obtain an equilibrium statistical mechanics-like description over the lattice configurations. This formulation allows macroscopic observables in space-time, such as the total jam delay and jam relaxation time, as well as microscopic jam statistics, to be expressed in terms of geometric properties of the height function. We thereby derive the associated scaling forms and recover the critical exponents previously observed in numerical studies. We discuss the physical implications of this space-time geometric approach.
We establish the global well-posedness theory of small BV weak solutions to a one-dimensional compressible Navier--Stokes model for reacting gas mixtures in dynamic combustion. The unknowns of the PDE system consist of the specific volume, velocity, temperature, and mass fraction of the reactant. For initial data that are small perturbations around the constant equilibrium state $(1, 0, 1, 0)$ in the $L^1(\mathbb{R}) \cap {\rm BV}(\mathbb{R})$-norm, we establish the local-in-time existence of weak solutions via an iterative scheme, show the stability and uniqueness of local weak solutions, and prove the global-in-time existence of solutions for initial data with small BV-norm via an analysis of the Green's function of the linearised system. The large-time behaviour of the global BV weak solutions is also characterised. This work is motivated by and extends the recent global well-posedness theory for BV weak solutions to the one-dimensional isentropic Navier--Stokes and Navier--Stokes--Fourier systems developed in [T.-P. Liu, S.-H. Yu, Commun. Pure Appl. Math. 75 (2022), 223--348] and [H. Wang, S.-H. Yu, X. Zhang, Arch. Ration. Mech. Anal. 245 (2022), 375--477].
Given a Joyce structure, we show that the associated $\mathbb{C}^*$-family of non-linear connections $\mathcal{A}^ε$ can be gauged to a standard form $\mathcal{A}^{ε,\text{st}}$ by a gauge transformation $\hat{g}$, formal in $ε$. We show that the corresponding infinitesimal gauge transformation $\dot{g}=\log(\hat{g})$ has a convergent Borel transform, provided $\dot{g}$ vanishes on the base of the Joyce structure. This establishes the first step in showing that such a $\dot{g}$ is resurgent. We also use $\hat{g}$ to produce formal twistor Darboux coordinates for the complex hyperkähler structure associated to the Joyce structure, and show a similar result about convergence of the Borel transform of the formal twistor Darboux coordinates.
We prove a conjecture of Fyodorov and Keating on the supercritical moments of the partition function of the C$β$E field or equivalently the supercritical moments of moments of the characteristic polynomial of the C$β$E ensemble for general $β>0$ and general real moment exponents. Moreover, we give the first expression for all correlation functions of the $\mathsf{Sine}_β$ point process for all $β>0$. The main object behind both results is the Hua-Pickrell stochastic zeta function introduced by Li and Valkó.
The Cahn-Hilliard-Navier-Stokes (CHNS) system utilizes a diffusive phase-field for interface tracking of multi-phase fluid flows. Recently structure preserving methods for CHNS have moved into focus to construct numerical schemes that, for example, are mass conservative or obey initial bounds of the phase-field variable. In this work decoupled implicit-explicit formulations based on the Discontinuous Galerkin (DG) methodology are considered and compared to existing schemes from the literature. For the fluid flow a standard continuous Galerkin approach is applied. An adaptive conforming grid is utilized to further draw computational focus on the interface regions, while coarser meshes are utilized around pure phases. All presented methods are compared against each other in terms of bound preservation, mass conservation, and energy dissipation for different examples found in the literature, including a classical rising droplet problem.
This paper addresses uphill transport (defined as a regime in which particle flow is opposite to the prescriptions of Fick's diffusion) in drift-diffusion particle transport constrained by volume exclusion. Firstly, we show that the stationary hydrodynamic limit of a multispecies, weakly asymmetric exclusion process (SHDL) naturally predicts precisely characterized uphill regimes in the space of external drivings. Then, with specific reference to systems of oppositely charged particles, we identify well-defined model hypotheses and extensions whereby the SHDL converges to the modified Poisson-Nernst-Planck model, thus bridging the gap between exclusion-based particle models and continuum descriptions commonly used in engineering. The merits and limitations of the models in describing the particle fluxes and predicting uphill transport conditions are investigated in detail with respect to the adopted approximations and simplifications. The results demonstrate the persistence of uphill transport phenomena across modeling scales, clarify the conditions under which they occur, and suggest that uphill transport may play a significant role in nanoscale electrolytes, confined ionic and iontronic devices, and membrane-based technologies.
We study the one-dimensional Helmholtz equation with (possibly perturbed) quasiperiodic coefficients. Quasiperiodic functions are the restriction of higher dimensional periodic functions along a certain (irrational) direction. In classical settings, for real-valued frequencies, this equation is generally not well-posed: existence of solutions in L 2 is not guaranteed and uniqueness in L $\infty$ may fail. This is a well-known difficulty of Helmholtz equations, but it has never been addressed in the quasiperiodic case. We tackle this issue by using the limiting absorption principle, which consists in adding some imaginary part (also called absorption) to the frequency in order to make the equation well-posed in L 2 , and then defining the physically relevant solution by making the absorption tend to zero. In previous work, we introduced a definition of the solution of the equation with absorption based on Dirichlet-to-Neumann (DtN) boundary conditions. This approach offers two key advantages: it facilitates the limiting process and has a direct numerical counterpart. In this work, we first explain why the DtN boundary conditions have to be replaced by Robin-to-Robin boundary conditions to make the absorption go to zero. We then prove, under technical assumptions on the frequency, that the limiting absorption principle holds and we propose a numerical method to compute the physical solution.
Our recent approach to the Finkelberg-Kazhdan-Lusztig equivalence theorem centers on the construction of a fiber functor associated with the categories in the equivalence theorem, which in turn explains the underlying algebraic and analytic structure of the corresponding weak Hopf algebra in a new sense. We provide a non-technical and historical overview of the core arguments behind our proof, discuss these structural properties, and its applications to rigidity and unitarizability of braided fusion categories arising from conformal field theory. We conclude proposing some natural directions for future research.
Let $Γ=q_1\mathbb{Z}\oplus q_2 \mathbb{Z}\oplus\cdots\oplus q_d\mathbb{Z}$, with $q_j\in \mathbb{Z}^+$ for each $j\in \{1,\ldots,d\}$, and denote by $Δ$ the discrete Laplacian on $\ell^2\left( \mathbb{Z}^d\right)$. We describe various algebraic properties of the ideal of spectral invariants for the discrete Laplacian when $d=1$, including a construction of a Gröbner basis. We also present various collections of complex $Γ$-periodic potentials $V$ that are such that $Δ$ and $Δ+ V$ are Floquet isospectral. We end with a discussion of the general setting, where the $q_i$ are taken to be vectors in $\mathbb{Z}^d$.
We show that a canonical procedure of extending generalized Dynkin diagrams gives rise to families of Kac-Moody groups that satisfy homological stability. We also briefly sketch some emergent structure that appears on stabilization. Our results are illustrated for the family {E_n} which is of interest in String theory. The techniques used involve homotopy decompositions of classifying spaces of Kac-Moody groups.
This paper studies structured products of real matrices for which the top Lyapunov exponent can be accessed by reducing the dynamics to an amenable generalization of upper triangular matrices. Exploiting prescribed zero patterns (including block-triangularity and sparse decompositions, conveniently encoded by a directed sparsity graph), we obtain explicit, computable bounds and, in favorable cases, formulas for $γ_1$ by combining deterministic triangular controls with a suitable refinement of the Furstenberg--Kifer lemma for block-triangular products. The estimates apply both to tempered (possibly deterministic) sequences and to stationary ergodic random cocycles under standard integrability. We also discuss applications to perturbation models for linear systems, including low-rank updates, where the reduction converts the problem to lower-dimensional or scalar cocycles.
We discuss the classical elliptic Toda chain introduced by Krichever and the elliptic Ruijsenaars-Toda chain introduced by Adler, Shabat and Suris. It is shown that these models can be obtained as particular cases of the elliptic Ruijsenaars chain. We explain how the classical $r$-matrix structures are derived for these chains. Also, as a by-product, we prove that the elliptic Ruijsenaars-Toda chain is gauge equivalent to discrete Landau-Lifshitz model of XYZ type. The elliptic Toda chain is also gauge equivalent to XYZ chain with special values of the Casimir functions at each site.
The Lorentz mirror model provides a clean setting to study macroscopic transport generated solely by quenched environmental randomness. We introduce a hierarchical version that admits an exact recursion for the distribution of left--right crossings, and prove normal transport: the mean conductance scales as (cross-section)/(length) for all length scales if $d\ge3$. A Gaussian approximation, supported by numerics, predicts that, in the marginal case $d=2$, this scaling acquires a logarithmic correction and that the variance-to-mean ratio of conductance converges to the universal value $2/3$ (the ``$2/3$ law'') for all $d\ge2$. We conjecture that both effects persist beyond the hierarchical setting. We finally provide numerical evidence for the $2/3$ law in the original Lorentz mirror model in $d=3$, and interpret it as a universal signature of normal transport induced by random current matching. A YouTube video discussing the background and the main results of the paper is available: https://youtu.be/G1nqKd6MiXo
A systematic understanding of integrability breaking in translationally invariant spin chains with genuine three-site interactions remains lacking. In this work, we introduce and classify minimal nonintegrable spin-$1/2$ Hamiltonians, defined as models that saturate injectivity while admitting no nontrivial local conserved charges beyond the Hamiltonian. We first rigorously establish the nonintegrability of the deformed Fredkin spin chain with periodic boundary conditions by mapping it to a nearest-neighbor composite-spin representation and excluding all admissible $3$-local conserved charges. Guided by its structure, we then construct five classes of spin-$1/2$ models with genuine three-site interactions. One class is integrable, while the remaining four contain exactly two interaction terms and constitute the minimal nonintegrable three-site models. Our results delineate a sharp boundary between integrability and nonintegrability beyond the nearest-neighbor paradigm.
Adaptive mesh refinement (AMR) is widely used to efficiently resolve localized features in time-dependent partial differential equations (PDEs) by selectively refining and coarsening the mesh. However, in long-horizon simulations, repeated intergrid interpolations can introduce systematic drift in conserved quantities, especially for variational discretizations with continuous basis functions. While interpolation from parent-to-child during refinement in continuous Galerkin (CG) discretizations is naturally conservative, the standard injection-based child-to-parent coarsening interpolation is generally not. We propose a simple, scalable field-conserving coarsening operator for parallel, octree-based AMR. The method enforces discrete global conservation during coarsening by first computing field conserving coarse-element values at quadrature points and then recovering coarse nodal degrees of freedom via an $L^2$ projection (mass-matrix solve), which simultaneously controls the $L_2$ error. We evaluate the approach on mass-conserving phase-field models, including the Cahn--Hilliard and Cahn--Hilliard--Navier--Stokes systems, and compare against injection in terms of conservation error, solution quality, and computational cost.
Spin resonance phenomena are conventionally described using transition probabilities formulated in a rotating frame, whose physical meaning implicitly depends on the choice of quantum reference standard. In this Colloquium, we show that a spin in a rotating magnetic field constitutes a configuration involving two quantum descriptions that share a common quantization operator but differ in their kinematic and dynamical roles. The transition probability therefore emerges as a relational quantity between quantum reference standards rather than an intrinsic property of a single evolving spin state. By incorporating the kinematic motion of the spin vector together with the dynamical evolution, this framework restores consistent energy accounting and reveals the dual-reference structure underlying spin dynamics in rotating magnetic fields.
We study the topological fundamental groups of the locally finite infinite ordered configuration space \(Conf^{lf}_\infty(\C)\) in the plane and the homotopy quotient of $Conf^{lf}_\infty$ by the canonical action of the infinite permutation group $\Aut(\N)$: \[ H^{lf}(\infty):=π_1^{\mathrm{top}}(Conf^{lf}_\infty(\C),\widetilde{\N}), \qquad B^{lf}(\infty):=π_1^{\mathrm{top}}\!\bigl(Conf^{lf}_\infty(\C)\!/\!/\Aut(\N),[e_0,\widetilde{\N}]\bigr). \] We prove that \(H^{lf}(\infty)\) and \(B^{lf}(\infty)\) are non-discrete and complete topological groups. A main structural theorem identifies \(H^{lf}(\infty)\) with a canonical locally finite inverse-limit model built from finite pure braid groups, and we construct a complete left-invariant ultrametric compatible with the quotient topology from the loop space of $\Conf$. The direct limit of finite pure braid groups admits a dense embedding into \(H^{lf}(\infty)\), and we show that \(H^{lf}(\infty)\) is the Raĭkov completion of this subgroup. Moreover, the direct limit of finite braid groups embeds into \(B^{lf}(\infty)\) and is dense in the finitary subgroup \(B^{lf}_{\mathrm{fin}}(\infty)\subseteq B^{lf}(\infty)\).