2026-05-11 | | Total: 23
We studied the constrained Hamiltonian formulation of a supersymmetric Korteweg-de Vries (KdV) equation, which is observed to be a constrained system similar to its classical version. We found a nontrivial Lagrangian description, where we select $a=2$ for the free parameter $a$ in the supersymmetric extension. The corresponding degenerate Lagrangian requires an exclusive consideration and the utilization of the Dirac-Bergmann algorithm. We explicitly determined the full set of primary and secondary constraints and constructed the total Hamiltonian governing the dynamics of the system. In this analysis, in addition to a nontrivial constraint involving the fermionic fields, the consistency conditions give rise to a nonlocal contribution to the Hamiltonian density. This highlights a distinctive feature of this supersymmetric extension. We showed that the resulting Hamilton equations of motion reproduce the supersymmetric KdV system in the component form. Finally, we derived a compact superspace representation of the Hamiltonian and demonstrated its consistency with the component-level formulation.
In recent work, we initiated a research program aimed at the systematic investigation of quantum superintegrable systems describing the interaction of two non-relativistic spin-$1/2$ particles in three-dimensional Euclidean space. In that study, we classified all such superintegrable systems admitting additional first-order scalar integrals of motion. In the present paper, we continue this program by focusing on systems that admit additional pseudo-scalar integrals of motion. Starting from the most general rotationally invariant Hamiltonian for two interacting spin-$1/2$ particles, we construct the most general first-order pseudo-scalar operator in the form of a matrix polynomial in the momenta. Imposing the commutativity of this operator with the Hamiltonian leads to a system of determining equations. By solving these equations, we obtain a complete classification of such superintegrable systems and determine the corresponding pseudo-scalar integrals of motion. The resulting classification provides new families of superintegrable systems with spin-dependent interactions. These systems enrich the class of integrable models relevant to nucleon--nucleon interactions and contribute to the broader program of classifying superintegrable quantum systems with spin. For selected cases, we further construct the associated polynomial symmetry algebras generated by the integrals of motion, providing additional insight into the algebraic structure of the systems.
We prove a volume-uniform effective-Hamiltonian theorem for bounded finite-range quantum spin systems on possibly infinite lattices. For any finite target region, we construct an energy-truncated Hamiltonian and prove a volume-uniform spectral-overlap bound controlling the leakage of its low-energy spectral subspace into the high-energy spectral subspace of the original Hamiltonian. The bound may contain non-exponential spectral-window terms, but its cutoff-dependent remainder decays exponentially in the cutoff. In finite volume, this yields stability of low-lying eigenvalues, with eigenvalue errors controlled by the exponentially small cutoff-dependent remainder. In infinite volume, we prove the corresponding spectral-overlap estimate in the GNS representation of an infinite-volume ground state. Thus, for bounded finite-range interactions, we extend and strengthen the effective-Hamiltonian mechanism of Arad, Kuwahara, and Landau by replacing the finite-volume operator-norm formulation with a volume-uniform spectral-overlap formulation applicable in the thermodynamic limit.
We consider an open quantum system governed $N$-body Lindblad equation and study mean-field limits in this setting. We prove that the $N$-particle dynamics converges, in the sense of quantum relative entropy, to the tensorized solution of the limiting nonlinear equation. More precisely, we establish explicit bounds of order $1/N$ on the relative entropy between the $N$-particle density operator and the corresponding product state, thereby providing a quantitative propagation of chaos.
We review how the local particle number cutoff introduced in [11] is used to build trial states for the dilute Bose gas that capture the substantial correlation structure of the ground state in the thermodynamic limit. In particular, we provide a simplified derivation of the Lee-Huang-Yang correction as an upper bound for the ground state energy.
The eigenstate thermalization hypothesis provides a framework for understanding thermalization in isolated quantum many-body systems by characterizing statistical properties of local observables in energy eigenstates. Here we demonstrate that distributions of matrix elements in macroscopic systems may depend not only on the macrostate parameters, such as the densities of local conserved charges, but generally also on the properties of ensembles used in sampling eigenstates. To this end, we depart from the conventional analysis of microcanonical windows and consider statistical ensembles with an adjustable scale parameter prescribing the magnitude of charge fluctuations. We specifically consider an integrable field theory that permits efficient numerical sampling of matrix elements and reliable extrapolation to the thermodynamic limit. Moreover, in this system, we identify a class of states that enables explicit closed-form computation of the suppression rate of matrix elements. Our findings reveal an underlying multiscale structure of matrix elements captured by a non-analytic fluctuation-scale dependence of algebraic exponents governing their statistical properties.
We consider the phase ordering problem for the low-temperature Ising dynamics initialized from a biased and disordered initialization. Work of Fontes, Schonmann, Sidoravicius (2002) showed that at zero-temperature, Ising Glauber dynamics on $\mathbb Z^d$ for $d\ge 2$ initialized from i.i.d. spins on each vertex that are $+1$ with sufficiently large probability, absorbs into the all-plus configuration quickly. We prove that analogous behavior holds throughout the low-temperature regime of the Ising model in two dimensions. Namely, there exists $p_0 <1$ such that Ising Glauber dynamics initialized from i.i.d. spins that are $+1$ with probability $p>p_0$, run at any low temperature $β>β_c$ converges rapidly to the plus phase measure $π^+$. The result is proved using a spacetime multiscale coupling valid in any $d\ge 2$, that boosts a uniform-in-$β$ quasi-polynomial bound on the mixing time of Ising dynamics with plus boundary conditions, into rapid phase ordering from biased initializations with no boundary conditions.
We prove that the KSGNS construction can be viewed as an endofunctor on a category whose objects are positive $C^*$-correspondences from a fixed $C^*$-algebra and morphisms are given by intertwiners which account for automorphisms of the fixed $C^*$-algebra. Using this perspective, we provide a functorial perspective for strict positive equivariant $C^*$-correspondences of $C^*$-dynamical systems and show every strict positive equivariant $C^*$-correspondence of $C^*$-dynamical systems unitarily uniquely dilates under the KSGNS construction to an equivariant $C^*$-correspondence of the dynamical systems.
We present a comprehensive Lefschetz thimble analysis of the one-dimensional Salpeter equation under a strong electric field. By treating the non-local square-root operator within the framework of algebraic analysis, we construct the full solution space, which includes relativistic generalizations of the Airy Ai and Bi functions and their negative-energy counterparts. Through a direct comparison with the Dirac and Klein-Gordon equations, we provide a geometric explanation for the absence of Klein paradox and the Schwinger effect in the Salpeter equation. Furthermore, our findings establish a unified geometric interpretation of the Schwinger effect across different relativistic wave equations.
Relaxation to equilibrium of a drifted Brownian motion is quantified by a probability density function, whose main (multiplicative) entry is an inferred Feynman-Kac kernel of the Schrödinger semigroup operator. Although seemingly devoid of a natural probabilistic significance (except for its explicit path integral definition), the pertinent kernel relaxes to equilibrium as well. The implicit Feynman-Kac potential ${\cal{V}}(x)$, continuous, confining and bounded from below, may take negative values. If positive, ${\cal{V}}(x)$ can be interpreted as the killing rate of the decaying diffusion process. In case of relaxing F-K kernels the killing effects are tamed (often overcompensated). The taming inavoidably appears in conjunction with the existence of the negativity subdomains of ${\cal{V}}(x)$ in $R$. If locally ${\cal{V}}(x) < 0$, its sign inversion $- {\cal{V}}(x)$ can be interpreted as the branching (cloning, alternatively bifurcation) rate in the course of the other wise free random motion. The arising killed diffusion processes with branching, we interpret as the possible path-wise background of tamed (relaxing) Feynman-Kac diffusions. We present acomputer-assisted path-wise arguments, towards a consistency of the killing/branching taming scenario, for a number of nonlinear model systems in one space dimension. Special attention is paid to Feynman-Kac potential shapes, presumed to be in the double well form, where an analytic access to eigenvalues and eigenfunctions is scarce beyond the semiclassical regime.
We study the classification of the gapped phases or massive renormalization group (RG) flows dual to the massless RG flows under changing the sign of the coupling constants. Whereas our formalism is based on combining Higgs- or Nambu-Goldstone-type arguments with Cardy's smeared boundary conformal field theories (SBCFTs), several puzzling structures arise. More specifically, the established Higgs or Nambu-Goldstone type arguments on the duality imply that the natural basis for the gapped states should be constructed from a set of smeared Ishibashi states, which are unphysical in boundary critical phenomena. Hence, the module of the gapped phases can be outside of that of boundary critical phenomena, whereas one can still calculate characterizing quantities by applying SBCFTs to the models. For example, we demonstrate that the massive RG flow dual to the massless RG flow from the tricritical Ising model to the Ising model, one of the simplest massless RG flows, has this unusual structure. This can be regarded as a quantum field-theoretic analogue of order-disorder coexistence in lattice models. More generally, the resultant gapped phases usually spontaneously break non-group-like symmetry (or noninvertible symmetry). Our work provides systematic quantum field theoretic descriptions of such unusual phases with spontaneous symmetry breaking of non-group-like (or noninvertible) symmetries.
We provide a rigorous construction of I.M. Singer's universal connection, a natural connection on a bundle of paths associated to any manifold, using the theory of diffeology. Furthermore, we generalize the universal connection to the diffeological setting, which enables the reconstruction of diffeological principal bundles with connections from their holonomy representations. We show that any two diffeological bundle-connection pairs with conjugate holonomy representations must be equivalent in a certain sense. These constructions are functorial in that, ultimately, our results can be summarized as an equivalence of categories between the so-called holonomy category and the category of diffeological bundle-connection pairs.
We prove a quantum ergodicity theorem for sequences of closed hyperbolic surfaces converging to the Poincaré disc in the Benjamini-Schramm sense. Assuming a uniform lower bound on the injectivity radius and a spectral gap, we establish vanishing of quantum variance on fixed spectral windows for a class of observables that contains differential operators and finite-propagation smooth operators. This generalises a result of Le Masson and Sahlsten from scalar observables to both local and 'pseudolocal' operator settings.
We study homeomorphisms of the circle that are smooth diffeomorphisms away from a finite set of $n$ points. These "broken diffeomorphisms" do not form a Lie group, but instead naturally assemble into a Lie groupoid. We construct an explicit nontrivial $n$-dimensional central extension of this groupoid, which restricts to the classical Virasoro group when confined to smooth diffeomorphisms. We further describe the associated "broken Virasoro" algebroid, defined as a nontrivial $n$-dimensional central extension of the Lie algebroid of vector fields on the circle that are smooth except at $n$ points. This construction generalizes the Virasoro algebra. As a byproduct, we analyze a related setting on an interval: we construct a nontrivial central extension of the Lie algebra of vector fields vanishing at the endpoints, together with the corresponding central extension of the group of diffeomorphisms fixing the endpoints. We also describe the associated Lie algebroid and groupoid obtained by allowing the endpoints to vary.
The $n$-species particle-exchange process (PEP($n$)) is an exclusion process in which particles of $n$ different species exchange positions on neighbouring sites with rates chosen such that the invariant measure on the discrete torus is a product measure. We address the large-scale hydrodynamic behaviour of this process which yields a system of $n$ coupled inviscid Burgers equations. This system of conservation laws is shown to admit Riemann invariants for arbitrary $n$ from which explicit solutions of the Riemann problem in terms of shock waves and rarefaction fans are obtained. We also introduce the open PEP($n$), in which particles are exchanged with boundary reservoirs. For a distinguished manifold of boundary rates, we prove that the invariant measure is the same product measure as in the periodic system. The hydrodynamic description in terms of Riemann invariants is used to derive the stationary phase diagram explicitly in terms of microscopic boundary rates. In the generic case, the steady state exhibits $2n+1$ phases, with boundary-induced phase transitions analogous to those of the single-species asymmetric simple exclusion process.
Motivated by the average partition function of c free bosons $($Afhkami-Jeddi et al. \cite{Afhk2021}$)$ and the average of the genus 1 partition function over the Narain moduli space $($Maloney-Witten \cite{Witten2020}$)$, we investigate ratios of theta functions. In this paper, we completely classify the minimizers (or maximizers) for ratios of theta and Epstein zeta functions. We find that the hexagonal lattice plays a pivotal role there. These results have direct applications in conformal and Liouville field theory via partition functions. Additionally, they yield the minima of differences of theta and Epstein zeta functions, which have implications for the mathematics of crystallization and interacting particle theory (\cite{Bet2016,Bet2019AMP}).
Single-Hessian Gaussian wavepacket dynamics (GWD) significantly reduces the computational burden of Heller's local harmonic GWD, while maintaining comparable accuracy in approximating vibronic spectra. Here, we provide a new, symplectic derivation of the equations of motion of single-Hessian GWD and show that, unlike the local harmonic version, this method conserves the non-canonical symplectic structure on the manifold of Gaussian wavepackets and$-$for bounded dynamics in smooth potentials$-$avoids the drift of energy. Our numerical results suggest that, despite being much more efficient than the local harmonic variant, the single-Hessian GWD exhibits the same $\mathcal{O}(\hbar)$ asymptotic error in averages of observables. To further accelerate numerical simulations, we implement high-order time-stepping geometric integrators that are time-reversible and conserve the norm and symplectic structure exactly, regardless of the time step. In addition, we present explicit expressions for the exact evolution of the width of a single-Hessian Gaussian wavepacket in a general potential, as well as for the exact evolution of the whole wavepacket in a global harmonic potential. Using on-the-fly ab initio Gaussian wavepacket dynamics on the first excited-state surface of ammonia, we numerically confirm the conservation of geometric properties by these integrators and demonstrate that high-order integrators can enhance both accuracy and computational efficiency. We also compute the photoelectron spectrum of the difluorocarbene anion and the absorption spectrum of methylamine, and find that, in comparison with experiment, single-Hessian GWD outperforms global harmonic models and matches the accuracy of local harmonic GWD. Finally, we identify which spectral features are sensitive to the choice of reference Hessian.
Certain Feynman integrals can be expressed as periods of differential forms on Calabi--Yau manifolds. We provide a mathematical proof of a result of Duhr and Maggio on the modularity of the two-dimensional conformal traintrack integral. Our approach is based on a factorization of the associated Picard-Fuchs system into a tensor product of Gauss hypergeometric systems via a gauge transformation due to Clingher, Doran and Malmendier.
We discuss the problem of counting certain LEGO structures, primarily those comprising parallel $w \times 1$ tiles. These can be combined, as a single LEGO structure, by interlocking the tiles. %Alternatively, if the interlocking condition is relaxed, so that tiles can also be placed end-to-end, a greater number of possible configurations results. We also study the historically earlier problem of counting the number of ways to combine $2 \times 4$ LEGO tiles, which in this case gives a 3-dimensional structure. In all cases the number of configurations is dominated by an exponential growth term, $μ^n$ where $n$ is the number of tiles. We present an algorithm for counting these various LEGO configurations, and use the data to estimate the asymptotics. We analyse the data so generated, and conjecture that, for the two-dimensional structures, the number of possible configurations grows like $A(w)μ(w)^n/n,$ and we give numerical estimates for $A(w)$ and $μ(w)$ for $w < 11,$ while for the three-dimensional structure the number of possible configurations is conjectured to grow like $Aμ^n/n^{3/2},$ where $μ= 117.25 \pm 0.05.$ We also study the sequences that arise when we fix the number of tiles $n,$ and vary the tile size $w.$ We prove that the sequences are polynomials of degree $n-1,$ and we give these explicitly for $n=1 \ldots 14.$
We consider the fluctuations of the free energy in generalized Sherrington-Kirkpatrick models and the log likelihood ratio of spiked Wigner models in the high temperature/subcritical regime. We prove that the limiting laws of the fluctuations are Gaussian under suitable assumptions, and the result is universal in the sense that it does not depend on the distribution of the disorder or the prior except that the means and the variances of the limiting laws depend on a few parameters of the model. The proof is based on the multigraph expansion that provides a unified approach to analyze both models.
We introduce a formalism for conformal field theory in four dimensions: a symplectic bi-Grassmannian representation of CFT$_4$ Wightman correlators. Working in Klein space with off-shell spinor-helicity variables, we show that correlators of $Δ= 2$ scalars and symmetric-traceless conserved currents are encoded by integrals over a pair of $n$-planes in a $2n$-dimensional symplectic vector space. These planes are constrained to be mutually symplectically orthogonal and aligned with the external kinematics. Conformal invariance, momentum conservation, and little-group covariance all follow geometrically from this structure. We derive all two- and three-point functions involving scalars, fermions, conserved currents, and stress tensors. As a non-trivial test, we show that the construction reproduces the full set of independent conformally invariant structures of $\langle JJJ\rangle$ and $\langle TTT\rangle$ in CFT$_4$. The resulting expressions are considerably more compact than their momentum-space counterparts. They also make manifest the double copy between Yang--Mills $\langle JJJ \rangle$ and Einstein-gravity $\langle TTT \rangle$. We further present a helicity-basis reformulation that makes the GL(1,R) and SL(2,R) weights of individual helicity components explicit. This basis also provides a natural starting point for a twistor-space formulation of the correlators.
We develop a variational thermodynamic framework for statistical systems governed by a self-referential nonlinear operator Omega characterized by structural exponents alpha > 0, beta >= 0, a symmetric kernel K, and a self-coupling constant kappa >= 0. The central object is the self-consistency entropy S[Psi] = -D_KL(Psi || Omega Psi), which vanishes at the fixed points of Omega and serves as a natural Lyapunov functional. Within the local kernel (mean-field) approximation, minimization of the free energy F = U - T S admits the Tsallis q-exponential distribution as an equilibrium state, with the entropic index q = alpha + beta emerging directly from the fixed-point structure of the operator rather than being postulated. The framework yields a consistent thermodynamic description, including a generalized equation of state PV = (2 - q) T, response functions, and a critical temperature associated with spontaneous symmetry breaking. The relation q = alpha + beta connects independently measurable structural exponents of the feedback mechanism to the observed tail index, providing a parameter-free criterion that distinguishes this approach from superstatistics, constrained entropy maximization, and q-deformed formalisms. This work establishes an operator-theoretic foundation for nonextensive statistical mechanics in which nonlinear self-referential feedback naturally generates Tsallis statistics in the mean-field limit.
We solve an equation of motion for the axial-vector field under boundary conditions of the bulk-to-boundary propagator in the framework of the hard-wall model of AdS/QCD. The equation is reduced to the form of a homogeneous ordinary differential equation with varying coefficients. We solve conjugate equations and find fundamental solutions. This allows us to establish main relations and necessity conditions. The integral equation has both Volterra and Fredholm terms and was solved by the iteration method. We apply a new scheme to solve the equation, since all linearly independent necessary conditions for the existence of a solution were defined and used to establish a sufficient condition for Fredholm solvability of the problem.