Many researchers have proposed their own potential-energy functions for the study of diatomic molecules, some of which have been proved to be equivalent. In this paper we suggest the application of a few simple rules to determine the equivalence between molecular potentials. The resulting systematic procedure is illustrated by means of several examples.

A classic theorem of Weyl (1921) states that a Weyl metric -- a natural generalisation of a pseudo-Riemannian metric -- is uniquely determined by its conformal and projective structures (i.e. by its conformal structure and its set of unparametrised geodesics). An equivalent formulation of Weyl's result is that a torsion-free linear connection compatible with a pseudo-Riemannian conformal structure is uniquely determined by its projective structure. We discuss analogous results for suitably defined notions of conformal structure for Galilei and Carroll geometry, i.e. for spacetime geometries arising as the `non-relativistic' and `ultra-relativistic' limits of Lorentzian geometry.

We identify the quadratic symmetry algebra of the two-dimensional Smorodinsky--Winternitz II system with a Laguerre-type confluent Heun algebra. The system is separable in Cartesian and parabolic coordinates. The complementary Cartesian separation operator \[ Y=\partial_y^2-\omega^2y^2+\frac{1/4-c^2}{y^2} \] is of Laguerre type, while the parabolic integral \(W=L_2\) is its algebraic Heun partner. With \(Z=[Y,W]\), the defining relations are \[ [Y,Z]=16\omega^2W-2bY,\qquad [W,Z]=6Y^2-4HY+2bW+8\omega^2(1-c^2), \] where \(H\) is central. This gives a direct superintegrable realization of the Laguerre--Heun algebra.

We show that axisymmetric Willmore surfaces admit a first-order formulation obtained by combining two independent first integrals. If $\rho$ denotes the distance from the axis of revolution and $\Psi=\sin\psi$, where $\psi$ is the tangent angle of the generating curve, then the profile satisfies \begin{equation*} \frac{\left[\Psi(\rho\Psi'-\Psi)^2+2(\rho\Psi'-\Psi)+2C_1\rho\right]^2}{1-\Psi^2} +\left[(\rho\Psi'-\Psi)^2-2\right]^2=C_2, \end{equation*} where $C_1$ and $C_2$ are constants of integration and the prime denotes differentiation with respect to $\rho$. This equation reduces the axisymmetric Willmore equation to a first-order ordinary differential equation and provides a convenient classification scheme for Willmore surfaces of revolution. The sphere and the Clifford torus are discussed as elementary checks of the formulation.

The aim of this paper is to revisit, in Souriau's book "Structure of Dynamical Systems", the chapter devoted to the geometric quantization where the justifications of important results and formulae are not given and are difficult to prove. After recalling the coadjoint orbit method and its application to the relativistic particle with spin, we state and prove two keystone theorems that allow to equip its prequantum manifold with the symplectic and contact structures. We apply them to the relativistic particle with spin 1/2, leading to the Dirac equation and the conservation of the probability current. We propose also an identity of conservation of the spin current and, invoking the Kaluza-Klein theory, a systematic construction of the symmetries of charge conjugation, parity transformation and time reversal which seems to us more convenient and readable than that of the classical presentations.

In this work we study the Schwinger-Dyson equations and saddle point equations of matrix integrals that come from type $(0,1)$ random fuzzy geometries coupled to fermions or bosons. Such random fuzzy geometries are bi-tracial Hermitian matrix ensembles with a determinant contribution in the integrand. We derive the Schwinger-Dyson equations using complex analytic techniques from the saddle point equation. For arbitrary potentials with either bosonic or fermionic contributions, their Schwinger-Dyson equations can be solved iteratively. For both the Gaussian models with either one boson or fermion we rigorously derive the formula for the free energy and first moment in terms of elliptic integrals. In the bosonic case this solution is closely related to the Hoppe model and the three-colour model.

We review the kinetic theory of one-dimensional nonlinear oscillator chains, of which the most famous example is the Fermi-Pasta-Ulam-Tsingou equation. We provide detailed, though not rigorous, accounts of the microscopic to mesoscopic, and mesoscopic to macroscopic limits: derivation of the kinetic wave equation and hydrodynamic limit. We also present the state of the art of the mathematical theory, including proofs. We discuss the connection to two famous problems of Mathematical Physics: the Fermi-Pasta-Ulam-Tsingou paradox, and the derivation of Fourier's law. Finally, many open problems and possible directions for future research are proposed.

We study the exact spectrum of the XX spin chain with constrained non-diagonal boundary fields, which can be analyzed by solving the associated Bethe Ansatz equations. In these equations, the number of Bethe roots has a definite parity, and all Bethe roots are located at the zeros of a unary function. We investigate the possible positions of the Bethe roots. Based on numerical observations, we analyze the Bethe root configurations for the ground state and the first excited state. Our results show that elementary excitations are characterized by the cooperative change of a pair of Bethe roots. Furthermore, we obtain an analytical expression for the ground state energy in the thermodynamic limit.

We investigate the Cauchy problem for the focusing modified Korteweg--de Vries (mKdV) equation with finite-genus algebro-geometric quasi-periodic initial data. By applying the nonlinear steepest-descent method of Deift--Zhou to the associated Riemann--Hilbert (RH) problem, we derive the long-time asymptotics of the solution in the critical regime where complex stationary phase points coalesce with the endpoints of the finite-genus branch cuts. The collision is resolved via a local Painlevé XXXIV parametrix, and the discrete spectrum (breathers) is incorporated into the analysis. The resulting expansion is valid uniformly up to an error of order $\mathcal{O}(t^{-1/2})$. In this critical region, the leading-order term comprises the finite-genus algebro-geometric background together with breathers, whose parameters are slowly modulated by the background solution.

We study $O(2)$-equivariant spectral flow for Dirac operators on a finite warped cylinder equipped with fixed admissible regularized APS boundary conditions. The twisting bundle is a real higher-rank orthogonal bundle, and reflection symmetry is implemented by a fiber involution. After complexifying the twisting bundle and diagonalizing the orthogonal twist, the Dirac equation decomposes into a scalar Fourier-mode radial equation, with moving rotating blocks and stationary neutral blocks. After regrouping conjugate and reflection-paired blocks, the crossing contributions define real $RO(O(2))$-classes. Consequently, we obtain an explicit blockwise formula for the $RO(O(2))$-valued spectral flow of the resulting regularized APS family. Under the standard self-adjoint Fredholm, endpoint-invertibility, and regular-crossing hypotheses, together with a fixed neutral-sector convention, this formula is obtained by assembling the local crossing contributions of the separated blocks. It refines ordinary integer-valued spectral flow and shows explicitly how the dimension map $RO(O(2))\to\mathbb Z$ loses representation-theoretic information. We also discuss the rank-three case, including the role of the fixed neutral sector, and the corresponding endpoint $\eta$/APS index interpretation.

The deterministic incompressible Navier-Stokes equations are physically incomplete: any viscous fluid at finite temperature must exhibit thermal fluctuations whose form is dictated by the fluctuation-dissipation relation. We formulate the stochastic Navier-Stokes equations with the kinematically selected deformation Laplacian on compact Riemannian manifolds with strictly negative Ricci curvature. The fluctuation-dissipation relation, derived from a topological (Poincaré lemma) argument, uniquely determines the noise from the viscous operator. For the spectrally truncated system, we prove that the unique stationary distribution is the Gibbs measure (Gaussian in the mode amplitudes, because the nonlinear convective terms preserve energy), and that convergence to equilibrium is exponentially fast with rate at least $2\nu\lambda_\Def$, where $\nu$ is the kinematic viscosity and $\lambda_\Def$ is the spectral gap of the deformation Laplacian. The spectral gap satisfies $\lambda_\Def \geq \kappa^2$ when $\Ric \leq -\kappa^2 g$, and is independent of the volume of the domain. On flat space, the analogous thermalisation rate vanishes in the infinite-volume limit. The equilibrium velocity-velocity correlation function decays exponentially in geodesic distance, in contrast to the algebraic decay on flat space. These results provide a rigorous statistical-mechanical foundation for viscous fluids on negatively curved manifolds and illustrate how the geometry of the domain controls not only the deterministic dynamics but also the approach to thermal equilibrium.

In the representation theory of Lorentzian orthogonal groups there are well known arguments as to why the parity operator $\mathcal{P}$ and the time reversal operator $\mathcal{T}$ should be realized as linear and anti-linear operators respectively (Wigner 1932). In this paper it is shown that the only operators satisfying the requisite properties for the spinor representation of the de Sitter group lead to fermion self-couplings which are necessarily parity odd, ruling out standard Dirac mass terms from any de Sitter symmetric theory of spin-$\tfrac{1}{2}$ fermions.

Systems moving in low Reynolds number fluid regimes are known to be governed by a ``motility map'' which linearly relates their shape change rates to they body frame velocity moving through the fluid. A consequence of this is ``Purcell's Scallop Theorem'' -- a locomotion system that undergoes shape changes that follow the same path forward and backward in time (reciprocal body deformations) cannot achieve net displacement, regardless of pacing of those this http URL show that linear-in-velocity motility maps extend to any power law viscosity (a.k.a. Ostwald--de Waele fluid), and therefore to many biological fluids in intermediate shear ranges. We also show that the linear-in-velocity property can be violated in Carreau-Yasuda fluids to produce net motion using an ``inchworm'' model consisting of two unequal masses with unequal drag coefficients performing reciprocal motions. Interestingly, the direction of motion can be switched by changing speeds. Our results show that the linear motility map of geometric mechaincs can be used to analyze and design locomotion in power-law fluids, and that some nonlinear drag relationships such as Carreau-Yasuda can be exploited to generate net locomotion in seeming violation of the ``scallop theorem''.

In this work we study the general relativistic massless Boltzmann equation on Friedmann-Lemaître-Robertson-Walker spacetimes with spatial topology $\mathbb{T}^3$ in the linear and decelerated expanding regimes, where the scale factor is $t^{\mathfrak{q}}$ with $\mathfrak{q}\in [0,1]$. The massless Boltzmann equation on these backgrounds admits non-stationary Maxwell-Jüttner equilibria of the form $\exp(- |t^{2\mathfrak{q}}p|)$. For $0 \leq \mathfrak{q} \leq 1$, we prove future global-in-time existence and uniqueness of small perturbations of these equilibria in the case of hard ball interaction without symmetry assumptions. For $0\leq \mathfrak{q} < 1/3$, we prove that the perturbation -- measured in a suitable $L^2_p$ based energy norm -- decays at the superpolynomial time-decay rate of $t^{-3\mathfrak{q}}\exp(-t^{1-3\mathfrak{q}})$, whereas for $1/3< \mathfrak{q} \leq 1$ we obtain the polynomial time-decay rate of $t^{-3\mathfrak{q}}$. In the borderline case $\mathfrak{q}=1/3$, we show the time-decay of $t^{-3\mathfrak{q} -c}$ with a uniform constant $c>0$. Finally, for $\frac{1}{3}< \mathfrak{q}\leq 1$, we prove future global-in-time existence and uniqueness of small perturbations of the vacuum solution on $\mathbb{T}^3$.

Stochastic inflation is a powerful technique for calculating the probability distribution function (PDF) of large inflationary perturbations, which may collapse to form Primordial Black Holes. The PDF, $P({\cal N})$, of the stochastic number of e-folds, ${\cal N}$, satisfies an adjoint Fokker-Planck Equation. We develop a new self-contained eigenvalue technique which can be used to determine $P({\cal N})$. First we apply this method to the simple case of quantum diffusion along a flat potential without any classical drift. We recover the expression for the PDF that has previously been found using characteristic functions, with an exponential tail. We also identify an intermediate regime between the peak and the exponential tail of the PDF, which has not been emphasized in earlier studies, where it exhibits a power-law behaviour, $P({\cal N}) \propto {\cal N}^{-3/2}$. Finally we apply the method to constant drift inflation, in the narrow- and broad-well limits. In the narrow-well limit, there is an analytic solution and the PDF is similar to the drift-free case, with a mildly suppressed tail. In the broad-well limit, determining the full set of eigenvalues and eigenfunctions requires a piecewise construction of the spectrum, and the broad-well PDF is qualitatively different, with an enhanced peak and a strongly suppressed tail.

The causal propagator (or Pauli-Jordan function), which multiplied by $i$ is the spacetime commutator of the field $[\phi(x),\phi(x')]$, plays an essential role in scalar quantum field theory. We discuss the role of the causal propagator and its spectrum in recent developments in defining quantum field theory in a more explicitly covariant manner, as well as in causal set theory. We then present a conjecture for its asymptotic spectral density in a free theory, and give examples that lend evidence to the conjectured scaling. Our work has implications for Lorentzian spectral geometry in much the same way as Weyl's asymptotic law has for Riemannian spectral geometry.

This work deals with the ergodicity for the stochastic dissipative Hasegawa-Mima equations in a bounded domain and perturbed by additive noise. The uniqueness of the invariant measure is obtained by the asymptotic coupling technique.

We study the fluctuations of logarithmic Green singularities in the spherical ensemble, viewed as a random discretization of the two-sphere. Smooth observables exhibit the usual Sobolev or Gaussian free field fluctuations, whereas logarithmic singularities live on a larger logarithmic scale and asymptotically decouple in high-dimension, producing an explicit white-noise limit. The result gives precise asymptotics for logarithmic potentials and characteristic polynomials, with constants expressed through chordal geometry on the sphere.

We develop a deterministic framework showing that a power-law form of semi-uniform localization of eigenfunctions (SULE) imposes strong structural and dynamical constraints on lattice operators. In particular, we prove that power-law SULE yields geometric constraints on localization centers, quantitative bounds on eigenfunction correlators, and power-law localization in the sense of finite moments of the position operator. Conversely, suitable bounds on eigenfunction correlators imply a corresponding form of power-law SULE, establishing a close connection between these notions. This highlights the role of power-law SULE as a structural mechanism governing localization beyond the exponential regime, including features typically associated with Anderson-type models. Our results reveal that power-law localization is intrinsically geometric: the spatial distribution of localization centers directly influences eigenfunction correlators and transport properties. As an application, we obtain power-law localization for long-range lattice operators with Stark-type potentials of sublinear growth whose spectral regime exhibits asymptotically collapsing spectral gaps and quasi-resonant structures, without relying on perturbative methods. Applications to long-range random operators are also discussed.

We combine Reimann's spectral typicality theorem -- a modern formulation of quantum ergodicity -- with the framework of Interval Quantum Mechanics (IQM). In IQM, quantum states are represented not by points but by \emph{quantum parcels}: weak open convex sets of density matrices defined by finitely many expectation intervals. Such parcels are the exact mathematical representation of the epistemic knowledge obtained from finite-precision measurements of macroscopic observables. We prove that for a single parcel in which every state has large effective dimension (a condition that ensures thermalization), the expectation interval of any bounded observable becomes concentrated around the microcanonical value for most late times. The asymptotic bound depends only on the minimal effective dimension within the parcel, not on its detailed shape. For a double parcel \((O_1,O_2)\) with both components contained in an energy shell, separated by a conserved quantity \(Q^*\) that is supported on the range of the measurement projector, we show that the expectation intervals of both parcels become concentrated near the microcanonical values of bounded observables, the separation is preserved exactly, and the updated double parcel after a fuzzy measurement remains valid.

We consider natural Hamiltonian systems with potentials that are $C^0$ or $C^1$ on a hypersurface and $C^{\infty}$-smooth in the complement and introduce and study corresponding notions of their integrabilty and non-integrability. As a motivating example, we derive and analyze models of bungee jumping. We provide prototype examples of the Liuoville-Arnol'd theorem for $C^0$ and $C^1$ Hamiltonians.

A spin-orbit Hamiltonian with an effective gauge structure carries two distinct loop objects that are routinely conflated: an energy-independent Wilson holonomy, which organizes interference and internal spin transport, and an energy-dependent monodromy, which quantizes the spectrum. We show that cleanly separating these objects supplies a precise, computable bridge between the loop/holonomy representation of gauge theories and condensed-matter spin-orbit transport. The construction maps a spin-orbit Hamiltonian to an effective $U(1)$ plus internal non-Abelian connection, reduces it to a first-order transport problem, and reads physical predictions from holonomy, monodromy, curvature, and eigenphase data. Two rings make the separation explicit. For a Dirac (graphene) ring with Rashba coupling and Aharonov-Bohm flux, the total holonomy factorizes exactly into a commuting $U(1)$ flux phase times an internal spin/pseudospin holonomy, and the spectrum follows from a holonomy-eigenvalue condition. For a Rashba-Dresselhaus ring, the internal $SU(2)$ transport is genuinely non-Abelian away from the $\alpha=\pm\beta$ pure-gauge locus, where curvature controls path ordering; spectral quantization then requires an explicit first-order reduction obtained by phase-space doubling of the second-order Schrödinger problem. A non-Abelian Stokes formulation and Magnus expansion serve as ordering diagnostics rather than spectral tools. Spin-network ideas enter only as historical geometric motivation, not as a dynamical import into spintronics.

We study the Riemann problem for the compressible Euler equations with a stationary coupling interface across which a discontinuity in the heat flux is prescribed. This coupling gives rise to non-conservative effects and models heat addition mechanisms such as condensation-induced waves. Without imposing restrictions on sonic states, we analyze the problem in all Mach number regimes. Lax weak entropy solutions are constructed via half-Riemann problems, and we show that non-uniqueness occurs for a large class of initial data. To address this, we introduce an admissibility criterion derived from the evolutionarity criterion, and we characterize the full structure of admissible Riemann solutions. Our analysis establishes local existence of admissible Riemann solutions provided the heat flux jump is sufficiently small, while also identifying families of initial data for which admissible Riemann solutions cannot exist for any fixed, nonzero heat flux jump. Numerical experiments are included to illustrate the theoretical findings.

The Laplace-Beltrami formalism, in which the Ricci tensor in the Einstein field equations (EFEs) is formulated at leading-order in terms of the partial-differential Laplace-Beltrami operator, was previously applied to coalescing compact binaries (CCBs) generating gravitational waves (GWs). Supposing that the CCB is an effective singular body -- a hollow mass-shell -- that follows a Kerr metric Ansatz, the EFEs were approached variationally such that the Ansatz geometric signature dictates the energetic output via $G_{\mu\nu}=8\pi GT_{\mu\nu}$. For the CCB mass-shell representation, the generated GW energy is treated as radiated surface energy via $E:=T_{00}V$. This surface energy yielded a close approximation to the cataloged GW coalescence energy, as previously shown in past comparisons. Given this success, it is logical to ask whether the Laplace-Beltrami formalism can be applied to other general relativistic systems, whether ``simple" or ``perturbative", beyond CCBs. This heuristic work focuses broadly on the EFEs themselves under the Laplace-Beltrami formalism, considering all differential orders up to second-order. This namely includes a deeper analysis on the variational methodology employed on the EFEs in the second-order sector, utilized in previous works, and the benchmark analysis of the lower first- and zeroth-order terms. This all-order report utilizes representative examples and select metric Ansätze to explore the formalism's practicality and its limitations; this is shown that the first-order decomposition showcases heuristically the mechanics of vector and scalar fields upon a curved spacetime.

We consider time-dependent Kohn-Sham equations in dimension $3$ with a fractional dispersion relation $(1-\Delta)^s$, $s\in(0,\frac32)$, and a class of interaction terms including, in particular, external potentials, internal potentials associated to Hartree-type non-linearities, and exchange terms described by energy subcritical pure-power non-linearities. We prove the local existence of weak solutions in $H^s$ using an approximation procedure regularizing the non-linearities. Assuming that the interaction energies can be controlled by the kinetic energy, we show that the solutions can be extended to global solutions using energy estimates. If $s\in[1,\frac32)$, we establish in addition the well-posedness of the time-dependent Kohn-Sham equations using Strichartz estimates.

We investigate the ill-posedness mechanism of the hard-sphere Boltzmann equation in $H_{x}^{s}$ Sobolev space. Via a direct construction, we prove a strong-weak type ill-posedness result in the low-regularity regime $s<1$, establishing a sharp threshold in connection to the local $s>1$ well-posedness result [11]. Instead of originating from the large-velocity growth of the collision kernel, this illposedness is generated by the loss term and dispersive effects. Consequently, we prove a dispersion-driven nonlinear instability mechanism for the hard-sphere Boltzmann equation, and provide a capstone of the ill-posedness series [18,20].

We consider integrable hierarchies such as KP, modified KP, 2D Toda lattice, BKP (small and large), DKP, Pfaff-Toda and their multi-component generalizations. We work in the framework of the bilinear formalism in which the universal dependent variable is a tau-function satisfying bilinear equations of the Hirota-Miwa type. Our principal interest in this paper is the dispersionless versions of the hierarchies. In the limit of zero dispersion the main object is an $F$-function, which is the limit of properly re-scaled logarithm of the tau-function. We show that in all the cases there exists an algebraic curve built into the structure of the hierarchy. We call it the {\it dynamical curve}. For the KP, modified KP and Toda lattice hierarchies, as well as for their multi-component generalizations, the curve is rational (of genus 0) and can be uniformized by rational or trigonometric functions. For hierarchies of the Pfaff type (DKP and Pfaff-Toda) the dynamical curve is in general a smooth elliptic curve (of genus 1), with its modular parameter being a dynamical variable. It is also shown that the large BKP hierarchy admits two different dispersionless versions. In one of them the dynamical curve degenerates to a rational curve while in the other one it remains to be elliptic. We show that a reformulation of the hierarchies based on uniformization of the dynamical curves by elliptic (or trigonometric) functions makes their structure nice and clear, especially in the multi-component case.

We define the analogue of Lie-Rinehart algebras over $C^\infty$-rings. We show that given a Poisson $C^\infty$-ring $\mathcal{A}$ its module $\Omega_{\mathcal{A}}^{1}$ of $C^\infty$-Kähler differentials is (part of) a Lie-Rinehart algebra. Conversely, given a Lie-Rinehart algebra $\mathcal{M} \xrightarrow{\rho} C^\infty\mathrm{Der}(\mathcal{A})$ over a $C^\infty$-ring $\mathcal{A}$, there is a natural Poisson bracket on the $C^\infty$-ring $\mathcal{F}(\mathcal{M})$ associated with the $\mathcal{A}$-module $\mathcal{M}$ (the $C^\infty$-ring analogue of an $\mathcal{A}$-algebra freely generated by the module $\mathcal{M}$). In the case where $\mathcal{A}$ is the $C^\infty$-ring of smooth functions on a manifold $M$ and $\mathcal{M}$ is the module $\Gamma(E)$ of sections of a Lie algebroid $E \to M$, the $C^\infty$-ring $\mathcal{F}(\Gamma(E))$ is the ring of functions $C^\infty(E^\vee)$ on the total space of the vector bundle $E^\vee \to M$ dual to the vector bundle $E$.

Starting from Maxwell's and linear momentum balance equations, we derive a ferrofluid model using the generalized Onsager's principle. Guided by a discrete perturbation estimate, we design and analyze families of Galerkin schemes that converge to sufficiently regular solutions and derive error estimates. Finally, we numerically explore the model with our proposed method.

Preconditioning is a fundamental technique for accelerating classical linear system solvers, and understanding when its benefits persist in quantum linear system (QLS) solvers is important for assessing the practical resource requirements of quantum linear algebra. In QLS algorithms, however, the potential advantage of preconditioning may be offset by the normalization overhead incurred by composing separate block-encodings of the system matrix and the preconditioner, as observed in recent work. This limitation leaves open whether additional algebraic structure can make preconditioning effective in quantum access models. Motivated by this question, we show that Pauli-structured representations of both the system matrix and the preconditioner allow the preconditioned operator to be accessed through regrouped Pauli expansions. In this setting, algebraic regrouping of Pauli products can reduce the Pauli coefficient weight of the preconditioned operator, thereby altering the normalization parameters relevant to quantum algorithms. We derive explicit size and coefficient-weight bounds for the regrouped Pauli representations, and we trace their consequences for both direct block-encoding constructions and randomized Pauli linear system solvers. These results identify when Pauli-structured preconditioning can reduce the effective complexity parameters of QLS algorithms, rather than merely improving the classical condition number. Numerical experiments on a finite-dimensional synthetic benchmark show reductions in norm-aware direct block-encoding diagnostics and in the randomized QLS per-sample depth proxy.

In this article we study periodic orbits of an electron attracted by a proton subject to Lorentz, electric, and Euler forces where each of them is allowed to depend periodically on time. This setup is motivated by the elliptic restricted three-body-problem where the Lorentz force corresponds to Coriolis force, the Coulomb force is replaced by the gravitational force, and the electric force of an external source is a combination of centrifugal forces and gravitational forces of other bodies. This is a singular version of a Euler-Hamilton system as discussed in [FW26b]. The singularity is due to collisions of the electron with the proton, respectively of two masses. Due to the possibility of collisions this problem has to be regularized. We show how periodic collisional solutions of this problem can be detected variationally in a non-local Lagrangian setup as well as in a non-local Hamiltonian setup.

Continual learning struggles to balance retaining past knowledge with absorbing new tasks. Stefan-CL elegantly resolves this stability-plasticity dilemma through the physics of melting. It frames consolidated knowledge as a protected "solid" and unused capacity as an adaptable "liquid." As the network learns, this boundary expands, governed by a "latent heat" tuning dial. By mathematically freezing the learned interior, Stefan-CL cuts forgetting to near zero, matching memory-heavy baselines without storing raw data, forging a beautiful, physics-grounded path for AI.

We found sets of exact analytic quasi-stationary states of a massive scalar field in a dyonic Kerr-Sen black hole~(DKSBH) background in the maximally extended spacetime region. A central novelty is the use of horizon-regular ingoing Eddington-Finkelstein coordinates, which enables a direct and unambiguous imposition of the ingoing boundary condition at the horizon. The exact radial solutions are in the form of confluent Heun functions. Imposing regularity at spatial infinity enforces a series truncation condition, yielding an exact quantization of the quasi-stationary frequencies. The spectrum exhibits a rich multi-branch structure, which we show splits into two distinct classes: modes that are insensitive to the black hole spin and charges and modes that explicitly depend on them. We uncover a clear asymmetry between co-rotating and counter-rotating configurations, driven by the spin-angular momentum coupling, as well as a systematic shift of the spectrum induced by electric and magnetic charges. The physical branches exhibit a universal behavior: modes with positive real frequency possess positive imaginary parts and therefore grow exponentially in time, whereas modes with negative real frequency are damped and decay. This suggests that positive-energy excitations in the region behind the outer horizon including the inner region of the inner horizon which contains the closed-timelike-curve, exponentially destabilize the background spacetime, supporting Hawking's chronology protection conjecture. In addition, the purely imaginary modes contain no oscillatory component and hence do not propagate through the spacetime, preventing traveling excitations along closed timelike curves and remaining consistent with the conjecture.

We introduce a novel exclusion process with two conservation laws, mass and energy, designed to mimic the essential features of continuous systems like interacting oscillators within the framework of interacting particle systems. This distinguishes our model from conventional multi-species processes where only particle numbers are conserved. As a basis for our fluctuation analysis, we first show that applying nonlinear fluctuating hydrodynamics (NFH) to this model reveals a wide variety of universality classes depending on the parameter choices. The main objective of this work is to study the stationary fluctuations of these conserved quantities. For a suitable choice of parameters, we rigorously show that the fluctuation fields converge to uncoupled stochastic Burgers equations (SBE) in the scaling limit. The proof relies on the second-order Boltzmann-Gibbs principle that we establish for this model, along with the spectral gap estimate and the equivalence of ensembles. Of independent interest is our general proof of the diagonalizability of the Jacobian matrix for the macroscopic current with distinct real eigenvalues. While this property is often taken as given in the physics literature, we establish it rigorously for multi-component systems even when the eigenvectors cannot be explicitly computed, offering a firm mathematical foundation for a broad class of models.

We derive a fully 3-dimensional Summation-By-Parts scheme for a class of linear wave equations on hyperboloidal slices that meet future null infinity on a Minkowski background. The scheme is derived in spherical polar coordinates, with a major strength being that it is provably stable and allows having grid points at the origin and on the $z$-axis, despite coordinate singularities, and at infinity, by introducing compactification followed by rescaling. Reducing it to the standard Cauchy problem, or on finite spacelike slices with an outer boundary, will follow a similar procedure. Interesting relations are obtained between the rescaling and compactification factors that simplify the equations, and the conditions on constraint addition terms are discovered to maintain symmetric hyperbolicity. Numerical implementation is achieved using finite-difference methods at second-order accuracy, which can be generalized to higher-order or spectral accuracies as well. Dissipation operators are given a more abstract treatment, which makes it possible to define them everywhere in the domain, including at the boundary points, in curvilinear coordinates, such that they satisfy the dissipative property (DP) in our energy norms. These generalizations reduce to the well-known Kreiss-Oliger dissipation operators whenever defined on a Cartesian grid in the bulk and satisfy the DP in the standard $L^2$-norms. We also propose new norm convergence tests that produce more accurate outputs. Promising results are obtained, giving hope for application to fully nonlinear systems, like the Einstein Field Equations, and extracting the resulting gravitational waves free of systematic errors or gauge ambiguities.

We construct a cohesive $\infty$-topos $\mathbf{H}_{\mathbb{Q}}$ equipped with a \emph{quantum modality} -- an idempotent product-preserving comonad $Q^{\diamond}$ with right adjoint $Q_{\bullet}$ satisfying the Beck--Chevalley compatibility conditions with the cohesive structure $(\Pi,\flat,\sharp)$. The model is the functor $\infty$-topos $\operatorname{Fun}(\mathbf{C}^{*}\mathbf{Alg}_{\mathrm{fd}},\; \mathbf{H}_{\mathrm{sm}})$, where $\mathbf{H}_{\mathrm{sm}}$ is the smooth cohesive $\infty$-topos and $\mathbf{C}^{*}\mathbf{Alg}_{\mathrm{fd}}$ is the category of finite-dimensional $C^{*}$-algebras with centre-preserving $*$-homomorphisms. Cohesion is lifted pointwise from $\mathbf{H}_{\mathrm{sm}}$; the quantum comonad is precomposition with the centre functor. We endow the topos with the Day convolution monoidal structure $\otimes_{\mathrm{Day}}$ induced by the tensor product of $C^{*}$-algebras and prove that $Q^{\diamond}$ is a strong monoidal comonad. The category of $Q^{\diamond}$-coalgebras is equivalent, via Gelfand duality, to the topos $\operatorname{Fun}(\mathbf{FinSet}^{\mathrm{op}},\mathbf{H}_{\mathrm{sm}})$ of discrete classical field theories. The comonad is interpreted as decoherence. This yields a cohesive linear $\infty$-topos in which the cartesian linear-logic structure degenerates, while the Day convolution provides a non-degenerate affine model of multiplicative intuitionistic linear logic. We also prove a synthetic no-cloning theorem and discuss the limits of the centre modality for representing quantum channels. This work provides the first rigorous instance of the cohesive linear framework and settles the open problem of finding a concrete model for cohesive linear homotopy type theory.

We investigate a three-dimensional ghostly Hamiltonian realisation of the fully degenerate resonant sixth-order Pais-Uhlenbeck oscillator. On the classical level, the phase-space flow is non-diagonalisable and decomposes into two complex-conjugate Jordan chains of length three, explaining the appearance of oscillatory solutions with secular terms. Upon quantisation, we construct intertwining operators whose quadratic combinations generate a hidden spectrum-generating $\mathfrak{u}(2,1)$-algebra. The associated descendant spaces are finite-dimensional invariant subspaces carrying non-trivial Jordan structure. Although these spaces admit a natural decomposition into irreducible modules of a distinguished $\mathfrak{sl}_2$-subalgebra, this decomposition does not in general coincide with the Jordan decomposition of the Hamiltonian. We further derive a tri-Hamiltonian formulation from Lie point symmetries of the classical flow and show that the corresponding Hamiltonians are naturally encoded by the same hidden algebra. Nevertheless, unlike in the non-resonant case, no positive-definite linear combination of them generates the same dynamics. Finally, we analyse the common centraliser of the tri-Hamiltonian family in $U(\mathfrak u(2,1))$, showing that the natural higher-order candidate $Q$ is reducible and yields no independent classical or quantum integral. The model thus provides a resonant higher-derivative system in which hidden $\mathfrak{u}(2,1)$ symmetry, classical and quantum Jordan structures, and multi-Hamiltonian geometry coexist.

Reduced-order models (ROMs) are widely employed to describe complex system dynamics when simulations with full-order models (FOMs) are computationally prohibitive. This study presents POD-AS-PRS, a novel model-reduction framework based on the active subspaces (AS) technique, which performs dimensionality reduction in both the state and parameter spaces, enabling efficient and high-fidelity approximations of quantities of interest (QoI). The approach employs proper orthogonal decomposition (POD) to extract low-dimensional coefficients from CFD snapshots, which are inputs to a residual neural network (ResNet) with linear layers to learn their nonlinear mapping to QoI. Reverse-mode automatic differentiation (AD) is utilized to compute gradients with respect to the coefficients, enabling AS analysis to identify influential modes by shifting the analysis to the POD coefficient space, thereby achieving a dual-stage dimensionality reduction driven by QoI sensitivity rather than modal energy. A surrogate model is subsequently constructed using a polynomial response surface (PRS) based on AS-derived active variables, retaining only the highly influential POD coefficients to ensure accurate and efficient QoI reconstruction. The framework is validated on periodic and chaotic bluff-body flows, demonstrating high accuracy with few influential parameters, while AD-based gradients achieve a two-order-of-magnitude speed-up over finite-difference approximations. Sensitivity analysis further reveals that the influential coefficients are not necessarily proportional to modal energy, highlighting the critical flow structures. Consequently, POD-AS-PRS identifies a low-dimensional manifold of sensitivity-dominant parameters that govern the QoI, elucidating the essential flow structures and their coupling with control parameters, thereby enabling efficient and accurate QoI reconstruction.

The Neural Spectral Element Method (NSEM) evaluates each network only at fixed Legendre-Gauss-Lobatto quadrature nodes and replaces all derivative calls with precomputed spectral differentiation matrices. The resulting deterministic loss enables limited-memory BFGS (L-BFGS) to reach residuals of 10^-9 to 10^-10. A Kosloff-Tal-Ezer coordinate map resolves electrochemical boundary layers, while a mesh-free neural mortar framework couples multi-element domains. On the four-example Poisson-Nernst-Planck (PNP) benchmark of Huang and co-workers, NSEM attains 10^-4 to 10^-7 relative pointwise error with two orders of magnitude fewer collocation points than the adaptive-resampling PINN baseline. Both a tanh multilayer perceptron (MLP) and a basis-aligned Legendre Kolmogorov-Arnold Network (KAN) backbone attain spectral accuracy within the same NSEM infrastructure, with the KAN requiring roughly half the Adam steps to enter the L-BFGS basin of attraction on the 1D PNP benchmark.

Density cube theory extends the canonical quantum density matrix $\rho_{ij}$ with the addition of an extra index to $\rho_{ijk}$. The elements of the density cube with two different indices, $\rho_{iij}$ and $\rho_{ijj}$, correspond to the real and imaginary parts of the off-diagonal element $\rho_{ij}$ of the density matrix and describe double-path interference, while those with three different indices describe non-canonical triple-path interference. In this letter, we propose an equation of motion for the density cube, obtained from the quantization of ternary Nambu dynamics, and find that pairs of triple-path interferences oscillate into each other.

We introduce the notion of timelike ideal boundary of a Lorentzian length space as the set of asymptotic classes of future or past-directed timelike geodesic rays, a construction complementary to the causal boundary in the sense of Geroch-Kronheimer-Penrose and akin to the concept of ideal boundary of a metric space. We endow such a timelike ideal boundary with a natural cone topology and an angular metric, and establish upper curvature bounds for the resulting metric space. Finally, we consider generalized cones as a model and study the relation between the timelike ideal boundary and both the metric ideal boundary of the fiber and the asymptotic behaviour of the warping function.

In this paper we define and compare several new Quillen model structures which present the homotopy theory of algebraic quantum field theories. In this way, we expand foundational work of Benini et al. by providing a richer framework to detect and treat homotopical phenomena in quantum field theory. Our main technical tool is a new extension model structure on operadic algebras which is constructed via (right) Bousfield localization. We expect that this tool is useful in other contexts.

The inverse source problem for the Helmholtz equation poses significant challenges, particularly when sources exhibit complex or discontinuous geometries. Traditional numerical methods suffer from prohibitive computational costs, while machine learning-based approaches such as Physics-Informed Neural Networks (PINNs) and the Random Feature Method (RFM) -- though computationally efficient for inverse problems -- lack the intrinsic machinery to handle the sharp morphological features in such singular problems, leading to inaccurate solutions. To address this issue, we propose the Morphology-Adaptive Random Feature Method (MA-RFM), a novel two-phase framework that adaptively locates critical regions and adds morphology activation functions for tackling the multi-frequency inverse source problem with complex geometry. Our framework recasts the ill-posed inverse problem into a well-posed, strictly convex optimization problem by reformulating the governing Helmholtz equation as a Tikhonov-regularized integral equation via its fundamental solution. In the first stage, the Integral Adaptive RFM (IA-RFM), employs an adaptive algorithm to rapidly localize the source support, thereby reducing computational overhead and accelerating convergence. In the second stage, posterior geometric information is progressively integrated into the solver via hybrid basis functions, enabling a precise reconstruction of complex morphologies. The MA-RFM extends the capabilities of RFM to handle PDEs with singular solutions while preserving its mesh-free efficiency. We demonstrate the superior performance of our approach through ample challenging 2D and 3D benchmark problems, even under limited and noisy measurement conditions, highlighting its robustness and accuracy in reconstructing complex and disjoint sources.

Surface holonomy plays a central role in higher gauge theory, bundle gerbes and the geometric formulation of Wess--Zumino terms in string theory. In this work, we consider the relation between surface holonomy and nonabelian multiplicative integration on surfaces. In this framework, we interpret the local Stokes law as a curvature obstruction law for higher holonomy and investigate its consequences in the abelian setting. We derive a global three-dimensional Stokes relation and show that it reproduces the familiar Wess-Zumino phase formula. In particular, the phase difference between two surfaces with common boundary is governed by the integral of the corresponding $3$-form curvature over an interpolating three-manifold. These results provide a geometric interpretation of multiplicative integration on surfaces in terms of surface holonomy and clarify its relationship with the classical theory of bundle gerbes and Wess-Zumino terms. We conclude by discussing possible extensions to nonabelian higher gauge theories and their relation to Wilson surface operators and generalized symmetries.

We consider the problem of reconstruction of an $n\times n$ matrix with coefficients depending rationally on $x\in \mathbb P^1$ from the data of: (a) its characteristic polynomial and (b) a line bundle of degree $g+n-1$, with $g$ the geometric genus of the spectral curve, represented by a choice of $g+n+1$ points forming a (non-positive) divisor of the given degree. We thus provide a reconstruction formula that does not involve transcendental functions; this includes formulas for the spectral projectors and for the change of line bundle, thus integrating the isospectral flows. The formula is a single residue formula which depends rationally on the coordinates of the points involved, the coefficients of the spectral curve, and the position of the finite poles of $L$. We also discuss the canonical bi-differential associated with the Lax matrix and its relationship with other bi-differentials that appear in Topological Recursion and integrable systems.

We study the ground-state of a Fermi gas with short range attrative interactions in one or two dimensions. N fermions are placed in a confining potential, and interact with each other through a negative potential, whose range is larger than the typical distance between particles. We show the convergence of the ground state energy of the Hamiltonian to a Thomas-Fermi energy in the large N limit. Furthermore, we prove convergence of the ground states, in the sense of their Husimi functions. These results extend to the case of a repulsive interaction of positive Fourier transform.

Although irregular vectors for the Virasoro algebra are widely used in modern mathematical physics, a rigorous existence and uniqueness theorem in arbitrary rank has not been available in the literature. In this paper, we develop an algebraic framework, based on Virasoro differential operators on the parameter space, which gives such a theorem for arbitrary integer and half-integer ranks. A key ingredient is the construction of a canonical operator \(L_*\) from the coefficient matrix of the vector-field part of a truncated Virasoro realization. This operator closes the recursive system by isolating the derivative with respect to the highest irregular parameter. Using this mechanism, we prove the existence and uniqueness of formal irregular vectors of arbitrary integer rank. We then construct the truncated Virasoro vector fields required in the half-integer rank setting and prove the existence and uniqueness of the corresponding half-integer rank formal irregular vectors. We also prove that, after a scalar gauge normalization, the canonical solutions satisfy the full lower Virasoro deformation equations. These results provide an algebraic foundation for the rigorous construction of irregular conformal blocks built from higher-rank irregular vectors.

We study the Klein-Gordon equation with general interaction terms, which may be linear or nonlinear, and space-time dependent. We initiate the study of such equations with large (non-radial) data. We prove that global solutions are asymptotically given by a free wave and a weakly localized part. The proof is based on constructing in a new way the Free Channel Wave Operator, and further tools from the recent works \cite{Liu-Sof1,Liu-Sof2,SW2020,SW2022}. This work generalizes the results of part of \cite{Liu-Sof1,Liu-Sof2} on the Schrödinger equation to arbitrary dimension, and non-radial data.

In this paper we introduce a multiparameter version of the quantum universal enveloping superalgebras introduced by Yamane in [H. Yamane, "Quantized enveloping algebras associated to simple Lie superalgebras and their universal $R$-matrices", Publ. Res. Inst. Math. Sci. 30 (1994), no. 1, 15-87]. For these objects we consider: - (1) their deformations by twist and by 2-cocycle (both of "toral type"); in particular, we prove that this family is stable under both types of deformations; - (2) their semiclassical limits, which are multiparameter Lie superbialgebras; - (3) the deformations by twist and by 2-cocycle (of "toral type") of these multiparameter Lie superbialgebras: in particular, we prove that this family is stable under these deformations, and that "quantization commutes with deformation".

Considered in this work is the Yang-Mills field in an extremal Reissner-Nordström black hole, a physically motivated mathematical model introduced by Bizoń and Kahl. The kink is a fundamental, strongly unstable stationary solution in this non-perturbative, variable coefficients model, with a polynomial tail and no explicit form. In this paper, we introduce and extend several virial techniques, adapt them to the inhomogeneous medium setting, and construct a finite codimensional manifold of the energy space where the kink is asymptotically stable. In particular, we handle, using virial techniques, the emergence of a weak threshold resonance in the description of the stable manifold.

We approach the problem of constructing a quantum analogue of the immensely fruitful classical transport cost theory of Monge from a new angle. Going back to the original motivations, by which the transport is a bilinear function of a mass distribution (without loss of generality a probability density) and a transport plan (a stochastic kernel), we explore the quantum version where the mass distribution is generalised to a density matrix, and the transport plan to a completely positive and trace preserving map. % These two data are naturally integrated into their Jordan product, which is called state over time (``stote''), and the transport cost is postulated to be a linear function of it. We explore the properties of this transport cost, as well as the optimal transport cost between two given states (simply the minimum cost over all suitable transport plans). After that, we analyse in considerable detail the case of unitary invariant cost, for which we can calculate many costs analytically. These findings suggest that our quantum transport cost is qualitatively different from Monge's classical transport.

The space discreteness hypothesis asserts that the nature of space at short distances is radically different from that at large distances. Based on the Bronstein inequality, here, we use a totally disconnected topological space $\mathcal{X}$ as a model for the space. However, we consider the time as a real variable. In this framework, the formalism of Dirac-von Neumann can be used. This discreteness hypothesis implies that given two different points in space, there is no continuous curve (a world line) joining them. Consequently, this hypothesis is not compatible with the theory of relativity. We propose $\mathbb{R}\times(\mathbb{R}\times\mathcal{X})^{3}$ as a model of a space-time. For simplicity, we work out our models using $\mathbb{R}\times(\mathbb{R}\times\mathcal{X})$ as the configuration space. Quantum mechanics (QM), in the sense of Dirac-von Neumann, on the Hilbert space $L^{2}(\mathbb{R}\times\mathcal{X})$ is a non-local theory: the Hamiltonians are non-local operators, and thus, spooky action at a distance is allowed. The paradigm asserting that the universe is non-locally real implies that the proposed version of QM admits realism. This version of QM can be specialized to standard QM by using Hamiltonians acting on wavefunctions supported on the region $\mathbb{R}\times\mathbb{R}$. We apply the developed formalism to the measurement problem. We propose a new mechanism for the collapse of the wavefunction. The mechanism resembles the one proposed by Ghirardi, Ramini, and Weber, but there are significant differences. The most important feature is that the Schrödinger equation describes the dynamics at all times, even at the moment of measurement. We also discuss a model for the two-slit experiment, where bright and dark states of light (proposed recently) naturally occur.

For site percolation on the triangular lattice, we define two lattice fields that form a logarithmic pair in the sense of conformal field theory. We show that, at the critical point, their two- and three-point correlation functions have well-defined scaling limits, whose structure agrees with that predicted for logarithmic field theories. One of the two fields can be identified with the percolation analog of the Ising energy field, while the other is related to the percolation four-arm event.

In quantum gravity, the gravitational path integral involves a sum over topologies, representing the joining and splitting of multiple universes. To account for topology change, one is led to allow the creation and annihilation of closed and open universes in a framework often called third quantization or universe field theory. We argue that since topology change in gravity is a rare event, its contribution to exponentially late-time physics is universally described by a Poisson process. This universal Poisson process is responsible for the plateau in the multi-boundary generalization of the spectral form factor at late times. In the Fock space of closed baby universes, this universality implies that the statistics of the total number of baby universes is captured by a coherent state. Allowing for the creation of asymptotic open universes calls for a noncommutative generalization of a Poisson process. We propose such an operator algebraic framework, called Poissonization, which takes as input the observable algebra and a (unnormalized) state of a quantum system and outputs a von Neumann algebra of a many-body theory represented on its symmetric Fock space. Physically, Poissonization is a generalization of the coherent state vacua of bipartite quantum systems or matrix quantum mechanics. The multi-boundary correlators of the Marolf-Maxfield toy model of baby universes, and the sum over bordisms of closed-open 2D topological quantum field theory, are entirely captured by Poissonization. In the Jackiew-Teitelboim gravity, the universal Poisson process reproduces the late-time correlators in the $\tau$-scaling limit.

We study the stability of ground states in the Edwards-Anderson Ising spin glass in dimensions two and higher against perturbations of a single coupling. After reviewing the concepts of critical droplets, flexibilities and metastates, we show that, in any dimension, a certain kind of critical droplet with space-filling (i.e., positive spatial density) boundary does not exist in ground states generated by coupling-independent boundary conditions. Using this we show that if incongruent ground states exist in any dimension, the variance of their energy difference restricted to finite volumes scales proportionally to the volume. This in turn is used to prove that a metastate generated by (e.g.) periodic boundary conditions is unique and supported on a single pair of spin-reversed ground states in two dimensions. We further show that any excitation above an infinite-volume ground state with a positive-density interface must have the property that its energy difference with the ground state in a restricted volume diverges as the square root of the restricted volume.

Strongly nonlinear dynamics, from fluid turbulence to quantum chromodynamics, have long constituted some of the most challenging problems in theoretical physics. This review describes a unified theoretical framework, the loop space calculus, which offers an analytical approach to these problems. The central idea is a shift in perspective from pointwise fields to integrated loop observables, a transformation that recasts the governing nonlinear equations into a universal linear diffusion equation in the space of loops. This framework, supported by recent mathematical analysis, is analytically solvable and yields an exact, parameter-free solution for decaying hydrodynamic turbulence--the Euler ensemble--which is shown to be dual to a solvable string theory. The theory's predictions include: (i) the unification of spatial and temporal scaling laws, governed by two related, infinite spectra of intermittency and decay exponents derived from the nontrivial zeros of the Riemann zeta function; (ii) a first-order phase transition in magnetohydrodynamic (MHD) turbulence; and (iii) the formation of quantized, concentric shells in passive scalar mixing. The appearance of identical mathematical structures as solutions to the turbulent regime of Yang-Mills gradient flow points to the broad applicability of this approach. The framework also yields a new type of analytic Hodge-dual matrix surface that solves the Yang-Mills fixed-point loop equation by harmonic map, opening the way for a geometric formulation of QCD string theory.

In this paper, we study the quantization of classical observables (i.e., functions) on a Kähler manifold $X$ as differential operators acting on holomorphic sections of tensor powers $L^{\otimes k}$ of the pre-quantum line bundle $L$. We prove two global results as follows. (1). For a general smooth function $f \in C^\infty(X)$, we construct higher order generalizations of Kostant-Souriau's pre-quantum differential operators using our Fedosov-type constructions of Bargmann-Fock sheaves in previous works. We prove that these differential operators are asymptotic to the Berezin-Toeplitz operators $T_{f,k}$ acting on the Hilbert space $H^0(X, L^{\otimes k})$ as $k \to \infty$. (2). If a smooth function $f \in C^\infty(X)$ is furthermore the symbol of a level $k$ quantizable function , then we prove that the associated Berezin-Toeplitz operator $T_{f,k}$ is a holomorphic differential operator. Conversely, Berezin-Toeplitz operators that are holomorphic differential operators all arise in this way. This gives a complete characterization of when Berezin-Toeplitz operators are holomorphic differential operators. To prove these results, we establish new orthogonality relations which generalize the classical Tuynman's Lemma, and employ various differential-geometric and analytic technqiues such as Hörmander's estimates.

We provide a numerical method for computing solutions to a free boundary problem arising from the equilibrium state of a floating drop. This numerical method is based on a Newton's method for the underlying nonlinear boundary value problems, and at each iterative step a Chebyshev spectral collocation method is employed. The problems considered here are those that can be described by using generating curves, and include problems in $\mathbb{R}^2$ and $\mathbb{R}^3$. The resulting nine-dimensional space of physical parameters is explored, and examples are given that highlight the potential energy of centrally located drops, wall-bound drops, and asymmetrical configurations in $\mathbb{R}^2$. Non-uniqueness of solutions to the corresponding Euler-Lagrange equations is displayed, and also strong evidence of non-uniqueness of energy minimizers is given.

Over an arbitrary commutative ring $R$, we develop a theory of quantum cellular automata. We then use algebraic K-theory to construct a space $\mathbf{Q}(X)$ of quantum cellular automata (QCA) on a given metric space $X$. In most cases of interest, $\pi_0 \mathbf{Q}(X)$ classifies QCA up to quantum circuits and stabilization. Notably, the QCA spaces are related by homotopy equivalences $\mathbf{Q}(*) \simeq \Omega^n \mathbf{Q}(\mathbb{Z}^n)$ for all $n$, which shows that the classification of QCA on Euclidean lattices is given by an $\Omega$-spectrum indexed by the dimension $n$. As a corollary, we also obtain a non-connective delooping of the K-theory of Azumaya $R$-algebras, which may be of independent interest. We also include a section leading to the $\Omega$-spectrum for QCA over $C^*$-algebras with unitary circuits.

Hamiltonians ${\cal H}^{a}_k$ of new integrable systems associated with the integer rays $(1,a)$ (commutative subalgebras) of Ding-Iohara-Miki (DIM) algebra in the $N$-body representation are closely related to commuting twisted Cherednik Hamiltonians $\mathfrak{C}_i^{(a)}$, ${\cal H}^{a}_k = \sum_{i=1}^N (\mathfrak{C}_i^{(a)})^k$. Moreover, symmetric combinations of eigenfunctions in the twisted Cherednik system were recently shown to produce the DIM Hamiltonian eigenstates. We explicitly construct these twisted Cherednik eigenfunctions recurrently by action of some (creation and permutation) operations. It resembles of a far-going generalization of Kirillov-Noumi operators, but exact relation remains to be specified.

Efficient classical simulation has matured to a critical component of the quantum computing stack, driving hardware validation, algorithm design, benchmarking, and the study of structured quantum dynamics. Lie-algebraic simulation ($\mathfrak{g}$-sim) offers a compelling approach: it represents Heisenberg-picture dynamics in the adjoint space whose dimension is set by the dynamical Lie algebra (DLA) governing the circuit, enabling efficient simulation of expectation values whenever the DLA grows only polynomially with system size. Despite this promise, existing applications of $\mathfrak{g}$-sim have been confined to free-fermionic settings. It has therefore remained unclear if the method can be applied to other structured circuit families, especially when their generators have large Pauli expansions, and hence whether Lie-algebraic simulability presents a genuinely broader paradigm than free fermions. In this work, we resolve this question by identifying additional non-trivial families of polynomial-dimensional DLAs and introducing symmetry-adapted bases that make the required adjoint-space preprocessing tractable. In particular, we develop an explicit Pauli orbit representation for permutation-equivariant dynamics, enabling efficient processing of cubic-dimensional algebras despite exponential Pauli support, and a modified generalized Gell--Mann representation for bounded Hamming-weight ($U(1)$-equivariant) dynamics, yielding polynomial simulation costs on fixed excitation sectors. Together with streamlined routines for free-fermionic algebras, these constructions significantly broaden the practical scope of $\mathfrak{g}$-sim as a unifying simulation tool for structured quantum circuits. Numerical benchmarks confirm favorable preprocessing scaling and validate large-scale proof-of-concept simulations beyond the reach of state-vector simulation.

We develop a mesoscopic formulation of equilibrium statistical mechanics based on coarse-grained occupation-number sectors of one-particle phase space. A mesoscopic partition function is constructed by averaging the microscopic Hamiltonian over configurations compatible with a given occupation profile. The construction converges to the canonical Gibbs partition function in the fine-graining limit and remains compatible with interacting many-body systems. Within this framework, thermodynamic extensivity is shown to be equivalent to asymptotic factorisation of the mesoscopic partition function, while residual inter-cell correlations generate subextensive corrections. The resulting formalism provides a mathematically consistent bridge between microscopic Gibbs theory and mesoscopic thermodynamics.

On metric measure spaces with sub-Gaussian heat kernel behavior in small time, we obtain a sufficient condition to solve Wick renormalized stochastic quantization equations with polynomial interaction. Given the power of the nonlinearity, the local solution condition depends on the Hausdorff dimension $d_h$, the walk dimension $d_w$, and the maximal spatial Hölder regularity of the heat kernel $\Theta$. A slightly more restrictive condition based on the same parameters is required for a global solution. For all global solutions, we construct an invariant measure for the Markov process defined by the solution. Our results apply to many rough spaces such as Barlow--Kigami type fractals as well as their Cartesian products and open up the possibility of making rigorous various structures in quantum field theory and statistical mechanics in non-integer dimensions. In the process, we build entirely from the short-time heat semigroup the necessary analytic framework that accommodates the issues which come with allowing rough local geometry.

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.

We investigate involutive, non-combinatorial solutions of the braid equation, viewing them as special deformations of the permutation map. Utilizing these solutions, we identify the associated quantum algebra and introduce it as the $\mathfrak{gl}_{k,m}$ Yangian. This newly derived Yangian is distinct from the known Yangian of the general linear Lie superalgebra; crucially, as a Hopf algebra, it possesses the standard tensor product algebra structure. The underlying algebra $\mathfrak{gl}_{k,m}$ is also introduced as a novel structure and constitutes a subalgebra of the Yangian. We then construct specific highest-weight modules of $\mathfrak{gl}_{k,m}$ that simultaneously yield the eigenstates of spin-chain-like ``Hamiltonians'', which are defined as the sum of the generators of the $A$-type braid group. Furthermore, we study the highest-weight representations and the corresponding combinatorial bases for $\mathfrak{gl}_{1,1}$, explicitly linking them to specific shapes of Young tableaux.

Standard quantum mechanics is an idealisation based on infinite-precision objects: point states, exact probabilities, and sharp measurements. Yet every real experiment has finite resolution, and for macroscopic systems we never have access to the microscopic state. Following Heisenberg's call for a theory built only on observable quantities, and von Neumann's insight that a complete description of a macroscopic system is neither possible nor necessary, we elevate the macroscopic state to a fundamental concept. We introduce Interval Quantum Mechanics (IQM), in which the state of a quantum system is never a point but a quantum parcel -- a convex weak open set of density matrices defined by finitely many open expectation intervals, representing exactly the set of microscopic states compatible with a finite set of macroscopic measurements. We show that unitary evolution lifts to a deterministic flow on parcels, and that a fuzzy measurement is represented by a volume-contracting update, strictly increasing the geometric information defined as the inverse of the Hilbert-Schmidt volume. By introducing a second impossible set we obtain a double-parcel whose geometric information increases monotonically -- dissolving the entropy stagnation problem, since von Neumann entropy is defined on point states that no finite-precision experiment can certify. The framework reformulates foundational puzzles without additional interpretational assumptions: wave-particle duality becomes a smooth trade-off; Schroedinger's cat is described by a parcel containing many microscopically distinct states; and spooky action at a distance disappears, replaced by a purely epistemic geometric update. All empirical predictions of standard quantum mechanics are recovered in the infinite-precision limit, which is never physically attained.

The incompressible Navier-Stokes equations contain viscous dissipation but no thermal noise. I show, using a topological argument based on Poincaré's lemma, that the fluctuation-dissipation relation for the full nonlinear dynamics can be derived without the linearisation or structural assumptions that all previous derivations require. The nonlinear convective term is Hamiltonian (energy-preserving and phase-space-volume-preserving) and drops out of the Fokker-Planck equilibrium condition exactly, so the noise derived from linearised fluctuations near equilibrium is in fact exact for the full nonlinear system. This result proves, rather than assumes, the reversible/irreversible decomposition that the GENERIC framework postulates, provided Poincaré's lemma holds on the phase space. The resulting stochastic system, with a physical molecular-scale spectral cutoff, is trivially globally well-posed: a finite-dimensional stochastic differential equation with non-degenerate noise and a confining Lyapunov function. It has a unique Gibbs equilibrium and converges to it exponentially. The difficulty of the Clay Millennium Prize Problem arises entirely from two idealisations, zero temperature and infinite spectral resolution, neither of which is satisfied by any physical fluid.

We study sharp weighted Sobolev-type inequalities of the form\[ \int_{0}^{1}|u(x)|\rho(x) \,\mathrm{d}x \leqslant \Lambda \Bigl(\int_{0}^{1}|u^{(k)}(x)|^2 \,\mathrm{d}x \Bigr)^{1/2}, \qquad u\in H_0^k(0,1), \]where $\rho$ is a non-negative weight. We characterize the minimizers and identify the optimal constant $\Lambda(k,\rho)$ by showing that every minimizer has a constant sign and therefore solves a nonlinear eigenvalue problem of polyharmonic type. This yields an explicit characterization of extremizers for a broad class of weights. Moreover, we even provide with a an explicit computation of the optimal constant in terms of the weight function. The new weighted estimates turn to be very useful and, among other applications, allow us to recover several previous sharp estimates and Hardy type inequalities on finite intervals.