2026-10-16 | | Total: 34
I investigate the average and variance of species abundance distributions for fixed-point solutions to generalized Lotka-Volterra ecosystems. Two different analytical methods are constructed. First, I build up on the perturbative Random Matrix Theory approach developed in my earlier work [Phys. Rev. E. 113, L062202 (2026)], and formalize this method diagrammatically. The abundance average is given by a sum of planar diagrams with non-crossing vertex pairings which can be encoded as Dyck paths. This allows me to express the combinatorial factor multiplying perturbative terms at each order as Catalan and Narayana numbers, following which the abundance average is given by the latter's generating functions. In this way, the infinite perturbation expansion is exactly resummed in a closed-form formula. I have not been able to find a similar resummation for the abundance variance, however diagrammatic Feynman rules are expressed, following which the perturbative series can be calculated to any desirable order. Second, I do the same calculations using the cavity method, which gives closed-form formulas for both the abundance average and variance. For the average, the cavity formula gives the exact same result as the perturbative approach when the average interspecies interaction vanishes. When this is not the case, the two results are functionally different, however they are hardly distinguishable in their parametric range of validity. Comparison with numerical data indicate that the cavity result is more accurate over a larger parameter range.
We announce the solution of constructing perturbative quantum BCOV invariants on arbitrary Calabi-Yau manifolds. In this note, we outline the construction and present an explicit solution. The quantization includes all descendant fields and all loop orders, and expresses the quantum BCOV interaction by an explicit integration over the moduli space of marked Riemann surfaces. This canonical solution provides a mathematical construction of all higher genus BCOV invariants on compact Calabi-Yau manifolds, which are conjecturally equivalent to all higher genus Gromov-Witten invariants on the mirror.
We prove a version of topological recursion for the correlation differentials of colored HOMFLY-PT polynomials of torus knots on the corresponding augmentation curves, following the proposal of Gu, Jockers, Klemm, and Soroush. To this end, we explicitly construct the spectral curve from the so-called auxiliary curve of Brini, Eynard, and Mariño and prove that it coincides with the augmentation curve for every torus knot with coprime $(r,s)$. Finally, our version of topological recursion proves the higher-genus closed-string prediction of the Aganagic--Vafa conjecture for all torus knot augmentation curves.
We study dynamical conservation laws and the propagation of thermodynamic structure in translation-invariant lattice quantum spin systems. Using Lieb-Robinson bounds [CMP 28 (1972), 251-257] in arbitrary dimension, and Araki's bounds [CMP 14 (1969), 120-157] in dimension one, we prove conservation of the specific entropy, specific energy, and specific relative entropy under suitable assumptions on the interactions. The transported local Hamiltonians admit surface energies that vanish per unit volume in the van Hove limit, and local observables remain a core for the transported generators at every finite time. Weak Gibbsianity is preserved for sufficiently small times in arbitrary dimension and for all times in dimension one. These results provide the analytical foundation for the structural theory of approach to equilibrium developed in a companion paper, and establish partial results toward the conjectures formulated there.
We continue the development, started in [AHP 25 (2023), 715-749], of a structural theory of approach to equilibrium in translation-invariant quantum spin systems. Dynamical Steady States (DSS) are defined as w*-limit points of the time-averaged evolution of an initial state. We show that if a DSS is an equilibrium state of a sufficiently regular interaction, then the specific energy of this interaction is necessarily a constant of the motion. Two distinct mechanisms lead to this conclusion, one based on dynamical surface energies and the other on conservation of specific relative entropy. We conjecture that these properties hold for a natural class of interactions. For a DSS admitting a physical equilibrium description, the associated equilibrium states are characterized by a dynamical Gibbs variational principle: they are precisely the minimizers of the specific relative entropy from the initial state over the physical equilibrium states. The achieved minimum equals the increase of specific entropy. The structural theory also yields characterizations of approach to thermal equilibrium, based on two model-dependent questions: restoration of a sufficiently local thermodynamic description at large times and classification of the corresponding local constants of motion. Once locality is restored, the constants of motion determine the possible asymptotic equilibrium states. In particular, if finite-range locality is restored and the energy is the only finite-range constant of the motion, then every DSS is a thermal equilibrium state at the temperature fixed by energy conservation. The finite-time dynamical input underlying this theory is established in a companion paper.
We study the two lowest Bloch bands of the infinite Chern band model in the semiclassical limit. When the Berry curvature parameter is an odd multiple of \(π\), we derive an explicit asymptotic formula for their exponentially small splitting, uniformly over the Brillouin torus. The proof combines Agmon estimates and WKB approximations for a matrix Schrödinger operator whose lowest potential eigenvalue has two wells. Geometric properties of the Agmon distance allow us to compute the leading tunneling coefficients by the method of stationary phase. We provide numerical computations to illustrate the asymptotic formula and support the conjecture that the two bands remain separated for other values of the Berry curvature parameter.
We prove Poisson local eigenvalue statistics for one-dimensional Schrödinger operators with nondegenerate Bernoulli potentials at any fixed nonzero coupling. After unfolding by the integrated density of states (IDS), the local eigenvalue process converges to a Poisson point process of unit intensity at every energy in the interior of the almost-sure spectrum. The theorem also holds for bounded nonconstant iid potentials, provided the IDS is continuous and strictly increasing near the reference energy. If the IDS also has a finite positive derivative at that energy, the usual linear normalization gives the same limit. The proof uses a phase coordinate with exact IDS drift to obtain a uniform exponential crossing law and joint Poisson counts in disjoint energy windows.
We give a new proof that translation-invariant free Fermion frustration-free models have zero Hall conductance. Our proof is analytic in nature and it uses the bulk-edge correspondence. The central argument, which holds beyond the translation-invariant situation, is that in frustration-free models, edge modes cannot merge into the ground state band. We give an example of a gapped frustration-free model that is nonetheless gapless on the half-plane.
Jellium, a system of electrons in a uniform neutralizing background, and the uniform electron gas (UEG), a system of electrons with constant density and no background, are two standard models of the electron gas. We prove that the ground-state energy per unit volume of Jellium, as defined by Lieb and Narnhofer, coincides with the UEG energy defined by Lewin, Lieb and Seiringer, i.e. \[ e^\hbar_{\rm Jel}(ρ_0)=e^\hbar_{\rm UEG}(ρ_0). \] This equality was previously established in the classical case by Lewin, Lieb and Seiringer and independently by Cotar and Petrache.
On a compact Riemann surface of genus $g \geq 2$, we study Simpson's stratification of the de Rham moduli space of rank-2 flat connections. Upon representing flat $\mathrm{SL}_2(\mathbb{C})$-connections by branched projective structures, or equivalently by opers with apparent singularities, the stratum containing a flat connection is determined by the minimal number of apparent singularities needed to represent it. We construct explicit holomorphic families in which two simple apparent singularities collide and disappear, so that the family crosses into a smaller stratum; this realizes Simpson's closure relation between strata near generic points. The analogous families of Higgs bundles are described by Baker--Akhiezer divisors on the spectral curve in which two points converge to a pair of points exchanged by the spectral involution. We show that, before the collision, the two apparent singularities are joined by twin paths bounding a bubble, so that these families arise from the isomonodromic bubbling surgery discussed by Calsamiglia, Deroin and Francaviglia.
The classical capacity of a noisy optical channel is the maximum rate at which information can be transmitted reliably, given a constraint on the mean input energy. For the ubiquitous bosonic Gaussian channels, which include standard models of loss and noise in fibre and free-space optical links, determining this rate requires addressing whether entanglement between channel uses or non-Gaussian signals can outperform Gaussian encodings. More precisely, the question is whether the energy-constrained Holevo information is additive across arbitrarily many channel uses and coincides with its Gaussian counterpart. We answer these fundamental questions in the affirmative for every single-mode bosonic Gaussian channel at every input energy. Namely, we prove that independently prepared Gaussian states with Gaussian-distributed random displacements achieve the capacity, including in the low-energy regime where squeezing the signals competes with encoding information in their displacement. Our proof connects an entanglement inequality to a sharp, additive lower bound on output entropy. We then characterise communication above capacity: under a strict total photon-number cutoff, the probability of correct decoding decreases exponentially with the number of channel uses, and we derive explicit bounds on the decay exponent. When only the mean photon number is constrained, a simple energy-rescaling formula determines exactly the entire asymptotic rate--error trade-off: the capacity at mean photon number $N$ and allowed average decoding error $ε$ equals the vanishing-error capacity at mean photon number $N/(1-ε)$. Together, these results identify the ultimate limits of classical communication for some of the most fundamental models of optical communication.
We introduce a quadratic ancestor inequality for weighted spanning trees on finite undirected graphs and give two applications. First, for every finite $n$-vertex transitive graph, the expected distance in its uniform spanning tree (UST) between two independent uniform vertices is at least $\sqrt{n}-1$, resolving a conjecture in [Benjamini--Kozma, Comm. Math. Phys. 259 (2005), 257--286]. Second, we prove that the growth exponent of three-dimensional loop-erased random walk (LERW) is at least 3/2.
Compactified Imaginary Liouville Theory (CILT for short) is a probabilistically constructed Conformal Field Theory (CFT), expected to have connections to minimal models, Coulomb gas integrals, and loop models. In this paper, we study the quantization of Segal's semigroup given by CILT. Each element $\mathbf T_f$ in the quantized semigroup corresponds to a conformal annulus: the unit disk $\mathbb D$ with $f(\mathbb D^\circ)$ removed, where $f:\mathbb D\to f(\mathbb D)\subset \mathbb D^\circ$ is a contracting biholomorphism. This semigroup encodes the local conformal structure of the CFT, in particular the Virasoro symmetry with central charge $<1$ for CILT. We first relate the quantized Segal's semigroup to Segal's amplitudes of annuli and recover the Virasoro generators $\mathbf L_n,\widetilde{\mathbf L}_n$ ($n\geq0$) of CILT as generators of this semigroup. Compared to the real Liouville CFT, the main difficulty here is the absence of positivity, which we address through estimates for the imaginary Gaussian Multiplicative Chaos and non-self-adjoint operators. Using Segal's gluing axiom as in the real case, we differentiate the amplitude and use conformal covariance to derive the Ward identities in terms of all Virasoro operators $\mathbf L_n,\widetilde{\mathbf L}_n$ ($n\in\mathbb Z$). We also outline the parallel theory for Boundary CILT (BCILT), concerning Segal's semigroup of semiannuli. These results provide essential tools for the conformal bootstrap program of (B)CILT, boundary BPZ equations, boundary structure constants, and the global construction of conformal blocks.
We establish local nonpersistence of full-dimensional invariant tori with polynomially decaying actions for the periodic quintic nonlinear Schrödinger equation \[ iu_t-u_{xx}+T_vu+μ|u|^4u=0,\qquad x\in\mathbb T:=\mathbb R/(2π\mathbb Z). \] For every $κ\geq1$, we construct a fixed bounded real Fourier multiplier $T_ve^{ijx}=v_je^{ijx}$ whose frequencies $λ_j=j^2+v_j$ satisfy the Bourgain-type Diophantine condition \[ \operatorname{dist}\left(\sum_j\ell_jλ_j,\mathbb Z\right) \geqγ_κ\prod_j(1+\ell_j^2\langle j\rangle^{2κ})^{-1} \] for some $γ_κ>0$ and every nonzero finitely supported integer sequence $\ell$, where $\langle j\rangle=\max\{1,|j|\}$. For this multiplier, every $3<M<4κ+\frac12$, every $A>0$, and all sufficiently small $μA^2>0$, the closed relative action tube of radius $3/10$ around $I_j=A\langle j\rangle^{-M}$ contains no full-dimensional invariant torus in $H^1(\mathbb T)$ carrying a complete flow, with arbitrary frequencies allowed. This advances the study of the polynomial-decay problem identified by Kuksin as a main open problem in KAM theory for one-dimensional Hamiltonian PDEs [Amer. Math. Soc. Transl. Ser. 2, 212 (2004), 237--258, Problem 7.1]. Bourgain discussed a related question for short-range Hamiltonian systems [Russian Math. Surveys 59 (2004), 231--246, Section 3]. The proof combines simultaneous Diophantine and weighted homological estimates, an exact common-phase transformation, and one normal-form step on four adjacent Fourier coordinates. Uniform phase rotation forces an action excursion, while overlapping parameter intervals cover all sufficiently small positive perturbations. The mechanism suggests applications to other Hamiltonian PDEs.
A fusion ring is said to be vanishing if one of its multiplication matrices has determinant zero. We classify all vanishing fusion rings of rank four. We then provide a more explicit classification in the integral case, which allows us to derive several new consequences: there are infinitely many nonisomorphic simple integral vanishing fusion rings of rank four and no simple integral fusion ring of rank four is either $1$-Frobenius or Burnside.
For any translation invariant generalized complex structure J on a torus we construct the derived category D(J) generalizing the derived category of coherent sheaves on a complex torus and the affine/convergent Fukaya category of a symplectic torus. The category is defined in terms of a non-commutative curved deformation of the Dolbeault algebra of a complex torus. We show that this category is independent of the chosen deformation path, which we prove using non-commutative Fourier-Mukai duality. We also prove that the action of O(2n,2n; Z) on the space of generalized complex structures leaves the category D(J) invariant. In particular, T-duality acts trivially on D(J). We construct concrete objects associated to affine Lagrangians in the derived category associated to a symplectic torus and show that their spaces of morphisms are given by Floer cohomology. We show that our construction is consistent with the homological mirror symmetry conjecture (including its non-commutative version).
We investigate a single-component Abelian-Higgs model with an eighth-order scalar potential whose vacuum manifold consists of two disconnected circles. The model supports two species of vortices and domain walls interpolating between the two vacuum components. We find type-1.5-like behavior characterized by nonmonotonic intervortex interactions and bound states of spatially separated vortices. Inner-vacuum vortices form clusters through short-range repulsion and long-range attraction, whereas outer-vacuum vortices form vortex-filled bubbles in which a closed domain wall confines multiple inner-vacuum-vortex cores. We construct the domain wall analytically, derive the asymptotic vortex-wall interactions, and analyze the fluctuation spectra of coincident multivortices to identify their splitting instabilities.
We study the stability of equal wealth shares in a spatial model combining local sharing, global redistribution, and random nonreciprocal interactions. On a two-dimensional square lattice in the infinite-system limit, the effect of local sharing depends on whether the influence of one site on another tends to have the same or opposite sign as the reverse influence. When these influences are nonnegatively correlated, diffusion reduces the redistribution required for stability, but a positive rate remains necessary. Negative correlations instead allow stability without redistribution, with the largest stability margin at finite diffusion and a return toward marginal stability at strong diffusion. This nonmonotonicity reflects a competition between spatial smoothing and the weakening of the stabilizing contribution of reciprocal anticorrelations. Equal shares can even remain stable under wealth-concentrating transfers, but only within an intermediate range of diffusion strengths.
We study the von Neumann operator-space entanglement entropy (OSEE) of initial product operators evolved under random gaussian dynamics, as modeled by a Brownian SYK$_2$ model. By exploiting a mapping to a domain-wall melting problem, we derive an exact expression for the time-dependent bipartite OSEE in the thermodynamic limit, and obtain an analytic characterization of the corresponding gaussian Page curve emerging at late times. Our results build on recent work on the entanglement dynamics in the all-to-all quantum simple exclusion process, leveraging tools from RMT, including the Coulomb-gas approach and a mapping to the free Jacobi process. Finally, based on our results, we discuss the possibility of relating the OSEE with non-gaussianity measures in the Heisenberg picture.
Algebraic entropy measures the exponential growth of the degrees of iterates of rational self-maps. It is introduced by Bellon and Viallet in \cite{BellonandViallet}. We study degree growth for rational maps arising from recurrences with the Laurent property, with particular emphasis on Gale--Robinson recurrences. We develop a direct method for determining the degree growth from the recurrence relation itself, reducing the problem to the study of associated max-linear, or tropical recurrences. We first obtain a closed formula for the degrees of iterates for a broad class of homogeneous quadratic Laurent recurrences. We then apply this result to the two-term Gale--Robinson family and derive an explicit quadratic growth formula with a periodic correction term. For three-term Gale--Robinson recurrences, we express the algebraic entropy in terms of the associated tropical growth rate sequence. Based on explicit computations, we formulate a conjecture for the general three-term family predicting quadratic growth with a periodic correction. We prove the conjecture for the infinite family with parameters $p=1$, $q=2$, and $r\geq3$.
We propose a general statistical framework for characterizing erratic functions through the statistics of their extrema, with the aim of identifying statistical signatures associated with chaotic processes. It is based on analyzing the extrema points of erratic functions associated with these processes. One is instructed to determine the probability distribution functions ( PDF)s of the locations of the extrema points, the spacings between the locations, in particular the nearest neighbors, and the ratios between adjacent spacings. We then fit these PDFs with PDFs known to relate to chaotic behavior like of random matrix theory. The extrema considered establish a cloud of points in general in $D$ dimensions for erratic function of several variables. For certain systems one finds also extrema curves. We further determine the ``thermodynamic" description of the erratic functions by computing the Shannon or differential entropy and by equating the PDFs to the Gibbs one. With this map we determine the ``energy function" and for certain cases also the ``temperature". We analyze three benchmark spacing PDFs: the Bernoulli shift map, the Poisson and the Gaussian beta ensemble and several others. We want to emphasize that the map to thermodynamics can be applied to variety of chaotic processes including non-physics ones. We distinguish marginal spacing distributions from joint distributions of consecutive spacings. The latter probe correlations and statistical organization that are invisible in one-spacing PDFs and may provide stronger signatures for distinguishing structured chaotic-like behavior from uncorrelated irregularity.
We study the connection between the $SU(N)_{K}$ Wess-Zumino-Witten (WZW) model and $SU(N)_{K'}$ WZW models, enforcing the unbroken $\mathbb{Z}_{N}$ symmetries. The fundamental underlying structures are the (half-)integer spin simple currents in the topological quantum field theories and anomaly matching conditions in conformal field theories. For the (half-)integer spin simple current condition of the folded models $SU(N)_{K}\boxtimes_{\mathbb{Z}_{N}} \overline{SU(N)_{K'}}$ or $SU(N)_{K}\boxtimes_{\mathbb{Z}_{N}} SU(N)_{K'}$, we obtain the generalized pairing condition $K= \pm R^{2}K' (\text{mod}.N)$ with $\text{gcd}(R,N)=1$. On the other hand, the anomaly matching condition results in the weaker condition $K= \pm RUK' (\text{mod}.N)$ with $\text{gcd}(R,N)=\text{gcd}(U,N)=1$. Surprisingly, the anomaly matching without the pairing condition results in the intrinsic mismatch of the preservation law of the $\mathbb{Z}_{N}$ defects and symmetries. The pairing condition canonically relates symmetries and defects, as in category theory, whereas it is unclear whether it can be a guiding principle in prohibiting the massless renormalization group flows.
Time-domain boundary measurements of diffuse light are temporal point-spread functions whose peak position is used as a robust scalar datum in diffuse optical and fluorescence imaging. We consider the exact response generated by a single fluorescent point target in a homogeneous half-space with a Robin boundary condition. The response is a temporal convolution of two diffusion Green-function factors containing a complementary-error-function boundary correction. For every positive diffusion, absorption, depth, and source--target and detector--target distance parameter, and for every nonnegative Robin parameter, we prove that the normalized exact response is a generalized gamma convolution density. Its extension by zero to the negative half-line is bell-shaped. A differentiable small-time expansion, established by sectorial analyticity and Cauchy estimates, together with the sign-change structure of the first two derivatives, excludes tangential critical points and yields a unique peak with strict increase before it and strict decrease after it. The argument is exact and requires no large-distance, source--detector balance, symmetry, restricted-time, or small-parameter assumption. The same conclusion holds after convolution with an exponential fluorescence-lifetime density for every positive lifetime. The key ingredients are an explicit Stieltjes representation of the Robin factor, hyperbolic complete monotonicity, closure of the generalized gamma convolution class, and bell-shape theory. A derivation of the exact half-space model and the differentiated endpoint asymptotics are included in two appendices.
We consider exponential last passage percolation in half-space with a deterministic first row. This model corresponds to a version of TASEP in half-space in which particles enter at the origin and are restricted by a moving wall. We determine the possible limiting distributions of rescaled last passage times depending on the boundary strength and the choice of the wall. As key inputs to our analysis, we establish several results for point-to-point half-space last passage times, including upper tail bounds, weak convergence, and bounds on the probabilities that geodesics touch the diagonal for certain start and end points.
We develop an operadic theory of anti-splitting and anti-replication for binary operads. For a binary operad $\opd$, we define its anti-bisuccessor $\ABSu(\opd)$ and anti-trisuccessor $\ATSu(\opd)$, with the property that an operation of $\opd$ splits into a negative sum of two or three operations, yielding canonical operad morphisms $\opd \to\ABSu(\opd)$ and $\opd\to\ATSu(\opd)$. For binary quadratic operads, the process gives Manin black products, \[ \ABSu(\opd) \cong \antiprelie \bullet \opd, \qquad \ATSu(\opd)\cong \antipostlie\bullet\opd. \] Dually, we introduce the anti-duplicator \(\ADup(\mathcal P)\) and anti-triplicator \(\ATrip(\mathcal P)\) and prove \[ \ADup(\opd^!)\cong \ABSu(\opd)^!, \qquad \ATrip(\opd^!)\cong \ATSu(\opd)^!. \] Consequently, the anti-duplicator gives a new way to calculate Manin white products. %, in particular, $\ADup(\opd)\cong \antiperm\circ\opd$. The framework recovers known anti-splittings, including anti-dendriform, anti-pre-Lie, and anti-pre-Poisson structures, in a uniform operadic form and produces further examples and dual pairs.