2026-08-14 | | Total: 34
We describe a mechanism of affine symmetry emergence in the maximal density regime of the prototypical model of self-organized criticality when the inverse square of the mesh of the underlying lattice is much larger than the number of random perturbation points distributed according to a prescribed probability measure supported in the interior of the ambient convex domain. Moreover, an appropriate scaling limit of the toppling function (aka odometer), which counts the number of operations per site, is a solution to a non-linear partial differential equation well known in the context of optimal transport and differential geometry, making it possible to accurately estimate the deviation of the density from its maximal value in any macroscopic window. The mechanism for the affine symmetry emergence is due to the novel empirical fact, supported in addition by inductive arguments that have recently being upgraded to a rigorous proof, that the scaling limit of the toppling function is the unique concave solution of the Monge-Ampère equation with Dirichlet boundary condition on the convex domain with the potential given by the probability measure used above as the infinite-perturbation profile.
Dyson universality typically manifests itself in local eigenvalue statistics. Here we show that its signature survives a nonlinear inverse-spectral reconstruction and reappears in the matrix geometry of the reconstructed operator. Using a dressing transformation, we map each unfolded spectrum to a deformation $f(x)$ of a fixed harmonic oscillator and represent it in the common oscillator basis by $F_{mn}=\bra m|f|n\ket$. We resolve the matrix-element weight into shells of fixed distance $d=|m-n|$, corresponding to the energy-transfer channels of the reference oscillator, and characterize the resulting distribution by distance-shell moments. Independently calibrated on Gaussian $β$-ensembles, these moments vary smoothly with $β$ and distinguish the GOE, GUE, and GSE. With this calibration fixed, applying the same diagnostic to the nontrivial zeros of the Riemann zeta function places the reconstructed operators in the GUE sector. Thus, the GUE character of the zeros is recovered not through direct statistics of the input levels, but from the distance-resolved geometry of the reconstructed operator.
This article contrasts the concepts of a proper condensate and the Onsager-Penrose criterion for Bose-Einstein condensation in the form of a debate between two proponents of the respective concepts. The prologue briefly introduces the two criteria. The epilogue contains remarks on properties of the underlying resolvent algebra that are physically relevant in this context.
We study Toeplitz determinants $\det T_n(e^f)$ for $f$ whose Fourier coefficients satisfy $f_k=O(|k|^{-1})$. This regime extends beyond $H^{1/2}$ and includes symbols with Fisher-Hartwig singularities. We develop an operator-theoretic approach based on the Baker-Campbell-Hausdorff formula that separates the quadratic term \[ \sum_{k=1}^{\infty}\min(k,n)f_kf_{-k} \] from the higher-order terms in the expansion of $\log\det T_n(e^{tf})$. We show that this quadratic term accounts for the possible growth with $n$, while every fixed higher-order coefficient remains bounded. For symbols with bounded positive and negative Fourier parts, our estimates yield two-sided bounds for the determinant after removal of the quadratic contribution. For a broader admissible class, including Fisher-Hartwig-type symbols, we obtain uniform higher-order coefficient bounds and a central limit theorem for the associated CUE linear statistics. We also obtain bounds on mixed exponential moments for CUE-derived random fields beyond the characteristic polynomial.
The most well-known beyond-Fourier heat conduction model, the Maxwell--Cattaneo--Vernotte equation is solved analytically in the presence of heat transfer boundary condition. In contrast to the corresponding Fourier problem, this boundary condition renders the underlying differential operator non-selfadjoint. With a suitable scalar product, the adjoint eigenvalue problem is established. The resulting left and right eigenfunctions constitute a biorthogonal system, allowing the expansion coefficients to be determined from arbitrary square-integrable initial conditions. This enables a convenient infinite-sum analytical solution, which is presented and thoroughly investigated for various values of the model parameters (including near-Fourier and highly hyperbolic regimes) and for two practically important initial conditions (equilibrium and flash pulse initiated). The spectral structure is analyzed in detail, including the occurrence of real, imaginary, and complex-conjugate eigenvalue roots, their asymptotic distribution, and their dependence on the dimensionless relaxation time and Biot number. We find good agreement with corresponding finite-difference numerical solutions. The comparison also illustrates the effects of spectral truncation (i.e., the Gibbs phenomenon) and numerical dissipation near propagating thermal-wave fronts. Completeness of the eigenfunction set is numerically demonstrated. The initial condition induced by the flash pulse is derived analytically. The results provide an analytical benchmark for non-selfadjoint hyperbolic heat conduction, unveil how boundary heat transfer and relaxation time influence the transition between Fourier-like diffusion and thermal-wave propagation, and open the possibility to find and investigate beyond-Fourier heat conduction via heat transfer in experiments and practical applications.
We study transparent boundary conditions (TBCs) for the time-dependent heat equation in a branching, quasi-one-dimensional domain modeled as a metric star graph. By combining the classical concept of TBCs with the theory of partial differential equations (PDEs) on networks, we derive an exact vertex condition that ensures unobstructed thermal flow across junctions. Specifically,we derive a sum rule for the diffusion coefficients that eliminates thermal backflow at the vertex.We demonstrate the validity of this analytical model numerically using a Crank-Nicolson finite-difference method, which confirms smooth, unobstructed heat propagation through the network.These results provide a practical mathematical framework for tunable control and optimization of thermal diffusion in low-dimensional structures. In particular, combining the well-known concept of the TBCs and theory for PDE on metric graphs, we propose a mathematical model providing a control tool for thermal diffusion in networks.
We construct a discrete Smorodinsky--Winternitz II superintegrable system on a triangular region of the two-dimensional square lattice. The model is built from a pair of commuting finite-difference number operators with finite spectrum together with an associated ladder-operator structure. We show that it is maximally superintegrable and that its symmetry algebra admits a Hahn-algebra presentation. The spectral problem is solved exactly in terms of bivariate orthogonal polynomials of mixed Krawtchouk and dual Hahn type associated with the factorized $A_2$-Leonard pair. Finally, we show that the continuum limit recovers the continuous Smorodinsky--Winternitz II system together with its Hermite--Laguerre eigenfunctions. We further explain how the Hahn presentation of the discrete symmetry algebra becomes singular in this limit, while the limiting algebraic structure is naturally described by the Laguerre--Heun algebra associated with Cartesian and parabolic separation of variables.
A compact formulation of electrostatics is presented. Without using constitutive relations, including the aether relations, an expression for the potential energy of a charged region under a potential function is proposed, where the charge distribution is specified by the charge potential field. Assuming that during a virtual motion of the charge, the virtual work expended by the field is equal to minus the time derivative of the potential energy, an expression for the mechanical force functional is derived. The force functional contains the action of an asymmetric active stress field D_{i}E_{j}, the skew-symmetric part of which is the mechanical couple density.
We construct a finite discrete realization of the Smorodinsky--Winternitz I superintegrable system on a triangular region of the two-dimensional square lattice. The construction is based on a pair of commuting number operators with finite spectrum together with an associated ladder-operator structure. We show that the resulting model is maximally superintegrable and that its symmetry algebra admits a Hahn-algebra presentation. Its spectral problem is solved exactly in terms of the bivariate dual Hahn polynomials of Tratnik type. Finally, we show that the continuum limit recovers the continuous Smorodinsky--Winternitz I system together with its ladder operators, eigenfunctions and symmetry algebra, thereby establishing the present construction as a genuine finite discrete realization of the continuous model.
For every Lie--Rinehart algebra, we construct an $F$-manifold algebra on the direct sum of its base algebra and module. Contrary to the assertion in \cite[Proposition 13.3.26]{LodayVallette}, the resulting structure is generally not Poisson. We determine when powers of the positive-degree ideal in the associated symmetric Poisson algebra are Poisson ideals, and relate the Leibnizator to the Lie--Rinehart differential. For a finite projective module of constant rank, the trace of the Leibnizator recovers the anchor and yields a rigidity result for injective anchors. We conclude with algebraic and geometric examples.
We construct rank-zero irregular vertex operators for the Neveu--Schwarz algebra as linear maps between irregular Verma modules of the same rank satisfying the usual superconformal commutation relations and an irregular asymptotic condition. Under the nondegeneracy assumption $Λ_{2p}\neq0$, we prove their existence and uniqueness. We then extend the decomposition theorem for the Neveu--Schwarz algebra to the irregular setting: the tensor product of the free-fermion Fock module with a rank-$p$ Neveu--Schwarz irregular Verma module decomposes into an infinite direct sum of tensor products of two Virasoro irregular Verma modules. We also prove a compatible decomposition of the irregular vertex operators into tensor products of Virasoro irregular vertex operators. Using suitable pairings and mode insertions, we derive bilinear differential equations for weighted sums of products of Virasoro irregular conformal blocks of types $(0,0,1)$ and $(0,2)$. After explicit parameter identifications and a gauge transformation in the Painlevé V case, the resulting bilinear differential operators agree with those appearing in the quantum Painlevé V and IV tau-function equations.
Uncertainty in nonlinear dynamical systems is often organized by transport structures that are not apparent from posterior geometry alone. We introduce Dynamical Uncertainty Geometry (DUG), a framework that separates three ingredients that are frequently conflated: posterior uncertainty, future observations, and model-induced transport. DUG indexes posterior credible sets and transport-conditioned subsets by a common enclosed-probability coordinate, allowing physically meaningful classes to be tracked across credibility levels rather than examined at a single threshold. The framework combines a mass-ranked credible filtration, a description of future experiments through their induced probability laws and Fisher geometry, and a labelled transport skeleton representing dynamically distinct outcomes. We establish stability results for the resulting mass-indexed persistence modules under simultaneous perturbations of posterior density and probability measure, including continuum, refinement, and finite atomic formulations. These results provide quantitative control of persistence computed from numerical approximations and weighted grids. Two benchmark problems illustrate the framework. In a short-arc orbit-determination setting, DUG identifies dynamically distinct return classes within a connected uncertainty region and highlights a difference between information-based and local Fisher-based observation-design criteria. In the Earth-Moon planar circular restricted three-body problem, DUG reveals multiple first-hit transport outcomes coexisting within a connected credible region and provides diagnostics for assessing their numerical resolution. Together, these examples show how topological summaries of posterior geometry, when coupled to transport labels and future experiments, yield a richer description of uncertainty than either posterior probabilities or dynamical classifications alone.
We prove that every PPT linear map has finite entanglement-breaking index, thereby establishing the eventual entanglement-breaking property of PPT channels in full generality. Furthermore, by utilizing completely positive maps with low entanglement dimensionality, we show that a large family of PPT maps, which strictly containing the class of 2-superpositive maps, has entanglement-breaking index bounded above by 3, uniformly in the dimension. In particular, these results provide strong evidence that the PPT-cubed conjecture may hold in full generality.
Quasi-hydrodynamics describes systems with quasi-conserved degrees of freedom, namely observables that relax on timescales that are finite but parametrically longer than microscopic relaxation times. Examples include kinetic chemistry and linear viscoelasticity. Here, we develop a rigorous effective-field-theory framework for linear quasi-hydrodynamics from kinetic-type theories. Starting from any linearized, causal kinetic-like theory endowed with slow degrees of freedom, we show that the exact dynamics of conserved and quasi-conserved observables admits a systematic expansion in the fast relaxation timescale. At zeroth order, the resulting equations form a causal, symmetric-hyperbolic theory belonging to the appropriate transient-hydrodynamic universality class, establishing Israel-Stewart-like dynamics as the universal description of slow relaxation modes. Higher-order corrections can be computed systematically and inherit universal symmetry, Onsager, positivity, and causality constraints from the underlying microscopic theory.
We identify a general mechanism for the finite-time thermalization of macroscopic observables, such as coarse-grained charge densities, in terms of elementary forms of the quantum dynamics of initial states: (1) aperiodicity, which provides a computable measure of (2) a dynamical partially ergodic exploration of the Hilbert space. Specifically, this mechanism predicts the equilibration of all (concentrated) macroscopic observables, in almost all states in an initial ensemble and almost all times within finite and longer intervals, given only the observable-independent information that the return probability of the ensemble of initial states is small over a finite time range. As a special case, it also accesses standard results on equilibration over infinitely long times in terms of (stronger versions of) the effective dimension of initial state delocalization in the energy eigenbasis. Our results incorporate macroscopic thermalization into the domain of operational quantum statistical mechanics, recently developed to provide finitely computable criteria for microscopic thermalization. We discuss an overall characterization of this approach as establishing connections between (1) the decay of a (theoretically or experimentally) computable probe indicating memorylessness, (2) a fundamental invariant mechanism in terms of the alignment of observables or states in the Hilbert space, and (3) predicting different natural forms of (classical and) quantum thermalization, most of which rigorously recover conventional eigenstate-based descriptions of infinite-time thermalization as a special case but provide stronger accessible predictions over finite observation times in the thermodynamic limit.
It has been known for a long time that the equation of motion for dissipative linear dynamical systems with constant coefficients cannot be derived from the classical principle of stationary action because the term proportional to velocity in the equation of motion leads to the time derivative of order one half in the Lagrangian. Thus, approaches utilising fractional calculus have been used; however, they suffer from deficiencies both from mathematical and physical points of view. We here present our version of such fractional calculus based approach that provides correct Euler-Lagrange and, ultimately, the Hamilton equations, energy change of the moving body, and an attempt for a geometric interpretation of how energy dissipates.
Recently we started building a bridge between two cases of Langlands duality. The latter is one of the most influential topics in mathematical research. It has many different appearances and influential subtopics. Yet there is a topic that until now seems unrelated to the Langlands program. That is the topic of invariant differential operators. That is strange since both items are deeply rooted in Harish-Chandra's representation theory of semisimple Lie groups. We started with the case of the group ~$SL(2n)$. In the present paper we deal with the group $SL(2n+1)$. The two mentioned groups are similar, but their representation theories is rather different.
The existence of global classical solutions for reaction-diffusion systems arising from chemical reaction networks remains a major open problem in the deterministic setting, especially for reactions with high polynomial growth. We prove that a suitably chosen, physically motivated transport noise yields unique strong solutions for complex balanced chemical reaction networks that are global in time with arbitrarily high probability. These solutions possess paths in $C^θ_t C^{\infty}_x$ for all $θ<1/2$ and are, in particular, classical in space. Furthermore, we show that a suitable transport noise can enhance the dissipation of spatial fluctuations at an arbitrarily prescribed exponential rate. Our proofs rely on a combination of scaling-limit arguments, maximal $L^p(L^q)$-regularity, and entropy-entropy dissipation estimates.
We construct three derived endofunctors on the derived category of graded integrable representations of the current algebra of a flat degeneration of $D(2|1,α)$, and prove that their classes in the Grothendieck group are the genus two Macdonald operators. Computing the graded $\operatorname{Ext}$ pairing explicitly, we deduce the self-adjointness of the genus two Macdonald operators and the orthogonality of genus $2$ Macdonald polynomials with respect to the certain skew-bilinear form.
The data processing inequality is a fundamental property that describes the loss of information through noisy channels. A more refined description is characterized by the strong data processing inequality (SDPI). In classical information theory, the tensorization of strong data processing inequality holds for a whole family of $f$-divergences. However, its quantum counterpart is less known. The tensorization of SDPI was shown only for some special cases previously, and the general understanding about the tensorization property of SDPI for quantum divergences remains open. In this work, we report two negative results: the tensorization property fails for certain quantum chi-square divergences, and it also does not hold for the quantum relative entropy.
We prove uniqueness in the Lorentzian Calderón problem in a semiglobal setting, comparing a potentially large perturbation of the Minkowski metric against a small one. Our proof uses a reduced amount of data, solving a formally determined version of the problem. In contrast to traditional approaches, it employs neither control-theoretic arguments nor a reduction to a geometric inverse problem via microlocal methods or high-frequency solutions. Instead, it relies on distorted plane waves and weighted $L^2$-estimates.
The main goal of tomography is the reconstruction of density function out of its line integrals (integral measurements of this density along x-ray lines). Such construction is possible and known as inverse x-ray/Radon transformations. We are interested in generalization of this problem to the case of particles. For this we need to limit the speed of these particles, otherwise the problem is reduced to the previous one. The question we discuss in this paper is how low can this speed be depending on the parameters of the studied potential. We study this problem using simple theoretical examples and Machine Learning pipeline to restore the potential.
Using patch-matrix-product methods, we study two reversible three-state interaction-round-a-face cellular automata introduced by Klobas and Prosen [J. Phys. A 55, 094003 (2022)] - the species-preserving and species-flipping rules - and a coherent quantum deformation interpolating between them. Within an explicit translationally invariant ansatz, the two classical rules share a one-parameter family of quasi-local conservation laws with an auxiliary-dimension-three realization. Every regular member is also the first centered logarithmic derivative of an analytic auxiliary-dimension-two family passing through the identity observable. On the closed span of this family and the three elementary local charges, the only nonzero Euler velocities are $\pm\sqrt{3/23}$. Consequently, $\sqrt{3/23}$ is a rigorous lower bound on the maximal Euler speed in any larger conserved sector. For the species-flipping rule, this value agrees with the reported extrapolated value $0.361$ of Klobas and Prosen, strongly supporting completeness of the known sound-active charges. The species-preserving rule admits, in addition, a staggered generating family whose logarithmic derivatives form an infinite tower of strictly local charges. In the quantum deformation, the same algebraic structures yield exact low-bond-dimension invariant states (scar candidates), exponentially long-lived quasi-local quasimodes, and diagonal quasi-local charges that produce a nonzero Mazur bound after projection away from the three elementary local charges.
We generalize the classical vector-valued tangential and normal trace theory on Lipschitz domains to the setting of non-Lipschitz $H^1$-extension domains, which includes domains with fractal boundaries such as the Koch snowflake. We define and study these operators based on the surjectivity of the trace operator for elements of $H^1(Ω)$ and the Hilbert structure of the associated trace space. The normal trace operator is defined on $\Hdiv$ in $\mathbb{R}^n$. A generalized Stokes formula allows us to introduce the tangential trace operator on $\Hcurl$ and $\Hc^1(Ω)$ in two and three dimensions. Following the approach of Buffa, Costabel and Sheen (2002) for Lipschitz domains, we define two abstract tangential boundary spaces as images of these trace operators, establish their Hilbert structure, and construct an abstract rotation operator linking them, in place of the geometric rotation based on the normal vector. We also extend Costabel's (1991) coercive bilinear form approach for Maxwell's equations to non-Lipschitz domains, yielding a theory that supports the treatment of the Hodge-Dirac operator and Green's formulas and enables the solution of boundary-value problems for the $\rotv \rotv +1$ operator within the $H^1$-extension domains framework.
We analyze the dynamics of two harmonic oscillators with nonlinear, non-Hermitian coupling between them. Specifically, we consider a synthetic antiferromagnet composed of two ferromagnetic layers with antiferromagnetic interactions between the two layers, leading to the emergence of two oscillation modes of the combined system. We analyze the response of this system to external driving fields. We determine the allowed steady states under various driving conditions. We then analyze the dynamics of small deviations away from the different steady states, which reveals the normal modes that can be probed using spectroscopic techniques in experiment. As expected, the nonlinear system exhibits multistability, in which multiple steady states can exist for the same driving conditions. We find that the dynamical response can be drastically different depending on which one of the two modes is driven and on the strength of the driving. Specifically, we can obtain level-repulsion or level-attraction patterns in the spectrum. In addition to the rich variety of spectra, we find that some steady states exhibit normal modes that contain multiple frequencies. Our results explain recent experiments on synthetic antiferromagnets and provide guidance for designing experiments that can explore previously unseen phenomena in these systems.