Papers

The Haldane-Shastry spin chain is a prototypical integrable model with long-range interactions, notable for hosting quasiparticles with fractional statistics and serving as a discrete analogue of a conformal field theory. Its remarkable simplicity is closely tied to a full Yangian spin symmetry. While the highest-weight states for this symmetry are known explicitly, a systematic treatment of the descendant states, needed for the computation of various physical quantities, has remained incomplete. In this work, we provide a detailed construction of these descendants in terms of the algebraic Bethe ansatz following recent work of Ferrando et al. In the limit of extreme twist, it includes the Gelfand-Tsetlin basis. As an application, we derive product and determinant formulae for norms and overlaps of these states.

In this work, we investigate the canonical quantization of topological solitons appearing in braneworld scenarios. In particular, we focus on Brane-Skyrmions, topological field configurations analogous to standard Skyrmions, which emerge as solutions of the Dirac-Nambu-Goto action supplemented by an induced curvature term. By quantizing the (iso)spin collective coordinates of the Brane-Skyrmion, we obtain a Hamiltonian that we solve perturbatively via an expansion in powers of $J^2$, in contrast to the standard Skyrme model. Furthermore, we implement a Physics-Informed Neural Network (PINN) to determine the soliton profile that minimizes the energy, consistently incorporating the backreaction from the quantized spin degrees of freedom. We conclude with a discussion of the potential applications of this framework to the description of hadronic spectra. Our results highlight both the theoretical potential of brane-defect models and the growing role of neural network methods in theoretical physics.

We investigate a physically motivated model of dynamical dark energy arising from the non-perturbative topological structure of the Quantum Chromodynamics (QCD) vacuum. The model introduces no new fundamental field or propagating degree of freedom: the dark energy (DE) density emerges as a global vacuum response to an expanding spacetime. We develop the first comprehensive cosmological implementation of this QCD-DE scenario and confront it with current observations, including Planck, ACT and SPT-3G cosmic microwave background data, DESI DR2 baryon acoustic oscillation measurements, and Type Ia supernova samples from Pantheon+ and DES-Dovekie. We compare the model with $\Lambda\mathrm{CDM}$ and $w_0w_a\mathrm{CDM}$ cosmologies. The model provides an excellent fit to the data and reproduces the late-time DE evolution preferred by DESI. The model naturally predicts effective phantom crossing behaviour at intermediate redshifts ($z\sim0.67$) while avoiding the instabilities associated with phantom scalar fields. Using goodness-of-fit statistics and Bayesian model-selection tools, including Akaike and Deviance Information Criteria and Bayesian evidence estimated from Markov-Chain Monte Carlo chains, we find that the QCD-induced model is consistently favoured over $\Lambda\mathrm{CDM}$ for the full combination of early and late-time datasets. Unlike the conventional descriptions of dynamical DE, support for QCD-DE in Bayesian evidence remains more consistent across datasets, suggesting that a physically motivated departure from a cosmological constant may provide a more economical description of the expansion history preferred by current observations.

We study the ultraviolet (UV) limit of finite-volume expectation values of vertex operators in the sine-Gordon model using the kink nonlinear integral equation (NLIE) description of the conformal limit. By analysing the integrable formulation of vacuum expectation values in the small-volume regime, we conjecture an explicit analytic expression for the leading asymptotic term in the small-volume expansion, formulated in terms of kink functions. This establishes a direct connection between the integrable finite-volume description and the expected conformal asymptotics determined by the 3-point functions of the underlying conformal field theory (CFT). The proposed formula is tested against the analytic expression known from complex Liouville conformal field theory using high-precision numerics, showing agreement to at least 19 significant digits.

We revisit the Coleman–Weinberg monopole problem introduced by Kiselev, where radiative symmetry breaking makes the broken vacuum metastable. We construct the associated static monopole–critical-bubble configuration in the full coupled radial Higgs–gauge system and show that it is a saddle of the static energy functional. The metastable monopole and monopole–critical-bubble branches are characterized by their profiles, energies, and radial Hessian spectra. The monopole–bubble solution carries a negative radial mode, while the metastable monopole remains locally stable until its lowest radial Hessian eigenvalue approaches zero. The resulting branch structure gives a direct static picture of how Coleman–Weinberg monopoles lose metastability, with critical rescaled scalar mass parameter \(\mu_c=0.064352(1)\).

In the study of duality cascades in three-dimensional gauge theories without FI parameters, an important role is played by a fundamental domain whose vertices correspond to brane configurations with vanishing relative ranks. Through the Fermi gas formalism, such brane configurations are known to be represented by factorized quantum curves. In this paper, we show that this factorized description extends naturally to quantum curves associated with del Pezzo geometries possessing exceptional Weyl-group symmetries in the presence of FI parameters. We find that the vertices of the fundamental domain corresponding to the weights of minuscule representations are realized as factorized quantum curves built from canonical operators interpreted as 5-branes dressed with FI parameters. This reveals an unexpected relation between factorized quantum curves and minuscule representations, providing a physical realization of these vertices in terms of “extremal” brane configurations.

We show that a general nonlinearly perturbed Maxwell theory of electromagnetism possesses three universal fundamental properties: (i) A finite-energy electric point charge. (ii) Exclusion of finite-energy magnetic monopoles and dually charged dyons. (iii) Arbitrary smallness of the effective radius of a point electric charge and the associated local undetectedness of the electric charge and energy. In particular, this last property offers a classical explanation for the invisibility of the electron, as a point electric charge, in accordance with the smallness of its effective radius. This nonlinear theory of electromagnetism has the feature that it minimally perturbs the Maxwell theory with a nonlinearity profile that is as general as possible such that the three universal properties stated above are all maintained.

We explore the possibility of the microscopic black hole formation at the Large Hadron Collider (LHC) in the context of the ADD model, for the centre-of-mass energy $\sqrt{s} = 14~\mathrm{TeV}$ with an integrated luminosity of $349.4~\mathrm{fb}^{-1}$ and provide constraints on the black hole mass, $M_{\mathrm{B}}$ by taking into account the effects of loss during the formation process of black holes through the parameter $\zeta$. Our analysis reveals that for $\zeta = 0$, black holes with $M_{\mathrm{B}} \leq 11.83~\mathrm{TeV}$ are disfavored in the case of three extra dimensions ($\mathcal{D}$), for the reduced Planck scale $\Lambda_{\mathcal{D}}$) of about a few TeV. The corresponding values for $\Lambda_{\mathcal{D}} = 9~\mathrm{TeV}$ turn out to be about $10.33~\mathrm{TeV}$. A significant reduction in the aforementioned limits is bserved while the loss gets higher, e.g. for $\zeta = 0.35$, $M_{\mathrm{B}}$ reduces to $7.65~(6.82)~\mathrm{TeV}$ at $\Lambda_{\mathcal{D}} = 1~(9)~\mathrm{TeV}$. These limits change to $12.03~(10.88)~\mathrm{TeV}$ and $7.80~ (7.03)~\mathrm{TeV}$ respectively for $\zeta = 0~(0.35)$ and $\Lambda_{\mathcal{D}} = 1~ (9)~\mathrm{TeV}$ at 95% CL in case $\mathcal{D}$ is raised to seven.

We investigate the impact of higher-curvature corrections on photon propagation within an effective field theory framework and their observational consequences in strong gravitational fields. We consider polarization-dependent modifications to photon trajectories in static and spherically symmetric spacetimes, focusing on Schwarzschild and Reissner–Nordström black hole backgrounds. Using the geometrical optics approximation, we derive the effective metrics governing photon propagation and study the resulting polarization-dependent shifts of the photon sphere. We compute the corresponding quasinormal modes in the eikonal limit and analyze their polarization dependence. We further investigate gravitational lensing, focusing on polarization-dependent corrections to the deflection angle in both weak- and strong-field regimes. In the strong-deflection regime, we find that even perturbatively small EFT corrections modify the coefficient of the logarithmically divergent part of the deflection angle, resulting in a potentially observable difference from the uncorrected case. This suggests that strong gravitational lensing may provide a sensitive probe of small higher-curvature corrections. While extracting EFT information directly from QNM frequencies is more subtle, QNMs may provide complementary information to gravitational lensing in future studies. Our results establish a framework for probing higher-curvature effects through polarization-dependent strong-field observables.

In this article, we present four different $1+1$D quantum models that realize the Yang-Lee (YL) phase transition under a deformation that preserves $PT$ symmetry. These are the antiferromagnetic Ising spin chain in transverse and longitudinal magnetic fields, the massive Schwinger model, the Blume-Capel model, and the three-state quantum clock model. Using the state-operator correspondence, we identify the YL critical point, compute the scaling dimensions of the lowest operators in each model, and find perfect agreement with the exact results for the YL criticality in two dimensions. Using bosonization for the Schwinger model and the Polyakov-Hubbard transformation for the other models, we show that in all of these quantum models the YL critical point is described, as expected, by a massless bosonic field with an $i \phi^3$ interaction. In the quantum clock model, this critical field interacts with a massive bosonic field, and we identify the massless and massive states in the Hamiltonian spectrum. In addition, we numerically compute the two-point function of $\phi$ at the Yang-Lee critical point and show that it grows with distance, in agreement with theoretical expectations.

A variational principle for quantum-state distinguishability is established in many-body systems with generalized symmetries, including noninvertible cases described by fusion categories. Standard fidelity and symmetry-resolved diagnostics emerge as coarse-grained limits of a more refined operational structure. When symmetry actions terminate at entanglement cuts, distinguishability is governed by boundary tube algebras within a symmetry-constrained measurement resource theory. The physically admissible instruments are characterized by complete positivity, entanglement-cut locality, boundary-module covariance, and sequential stability. The resulting optimal measurement structure is uniquely fixed by the center of the boundary tube algebra, $\mathcal{A}_{\mathrm{phys}} = Z\!\left(\mathrm{Tube}_{\mathcal{C}}(\mathcal{M}_A)\right)$, whose primitive dempotents define tube-sector probabilities that refine fidelity-based and symmetry-resolved descriptions. The associated tube positive-operator-valued measures (POVMs) are extremal and yield optimal one-shot hypothesis-testing distinguishability under symmetry constraints. The categorical selection rule is universal across fusion categories and independent of microscopic realization; microscopic models enter by choosing the boundary module and the protocol family used to witness a nontrivial tube fiber.

We derive the general cubic action for perturbations about a cosmological background in the quadratic subclass of degenerate higher-order scalar-tensor (DHOST) theories, focusing on scalar-scalar-tensor interactions. We adopt a fully covariant formulation and implement unitary gauge at the level of perturbations. This mixed sector provides the key ingredient for estimating the decay rate of a gravitational wave into two scalar excitations in the quadratic DHOST models of dark energy in the late-time Universe.

The (2+1)D U(1) pure gauge theory always exists in the confining phase, with strings of non-zero string tension giving a characteristic linear potential between static charges. This makes it a useful testing ground for quantum computing methods designed to study string dynamics of confining gauge theories. Here we implement a minimal U(1) quantum link model on a quantum computer with qubit degrees of freedom representing the dual height variables of the model. This facilitates an efficient realization of plaquette interactions and enables effective calculations of real-time dynamics that are inaccessible to traditional quantum Monte Carlo. A specifically tailored lattice geometry is chosen to match the heavy-hexagonal geometry of the IBM quantum hardware used here, minimizing non-adjacent qubit interactions. By performing quantum quenches from a simple initial string state, we probe the transverse quantum fluctuations of the string before it thermalizes. Our experimental results from digital quantum simulations, with up to 112 qubits, show good agreement with reference tensor-network calculations at short times and with thermal averages at long times. Near the phase transition, the quench dynamics exhibit large fluctuations of the initial string that extend across both spatial dimensions of the lattice. Nonetheless, our error-mitigated estimators from the quantum hardware also give accurate predictions in that regime, with noise-induced violations of local gauge symmetries comparable to finite-bond-dimension tensor-network results.

We revisit flavor oscillations in vacuum in terms of the propagation time of intermediate states. In the limit of a long propagation time (or distance), degenerate intermediate states exhibit oscillatory behavior, as described by the Jacob–Sachs (or Grimus–Stockinger) theorem within wave-packet quantum field theory. By explicitly evaluating the relevant integrals using the saddle-point method, we derive an extended expression for the flavor-changing amplitude that remains valid even for shorter propagation times. We show that oscillations occur only when the propagation time exceeds a threshold set by the energy uncertainty of the external wave packets and by the decay width of the propagating particle. For shorter propagation, the intermediate particle behaves as a purely virtual state, in the sense that it cannot propagate over a macroscopic distance. Although a direct experimental test of the transition from virtual to real propagation is challenging, since it typically occurs at microscopic scales, our result implies that the Jacob–Sachs theorem holds to higher accuracy than previously expected, even at short propagation times. Our formalism applies not only to neutrinos but also to other propagating particles, and future improvements in energy resolution may make this threshold observable.

We formulate reheating as a constrained variational problem in the space of equation-of-state histories, rather than attempting to describe it through microscopic models. We introduce a regularized functional framework that identifies reheating histories which extremize a given cosmological observable under minimal physical assumptions. As illustrative applications, we consider prompt gravitational waves, induced gravitational waves, and primordial black holes. We find that different observables select qualitatively different regions of reheating-history space. These examples demonstrate that cosmological observables define distinct extremal directions in reheating-history space and can therefore be used to systematically explore the space of post-inflationary expansion histories.

We present the construction of the generating function of single-centered black hole index in general $\mathbb{Z}_N$ CHL models. This is done by subtracting from the index of quarter BPS dyons, described by a meromorphic Siegel modular form, the generating function for the index of two-centered black holes. We use black hole bound state metamorphosis in CHL models for the construction of the generating function of two-centered black hole index. We prove the convergence of the generating function for the cases $N=2,3$.

We study Euclidean M2-brane instantons in Freund-Rubin backgrounds $\mathrm{AdS}_{4}\times \mathrm{Y}_7$. For a seven-dimensional weak $G_{2}$ manifold $\mathrm{Y}_7$, we show that the BPS condition for an M2-brane wrapping a three-cycle $\Sigma\subset \mathrm{Y}_7$ is equivalent to the associativity condition with respect to the nearly parallel $G_{2}$-structure. When $\mathrm{Y}_7$ is Sasaki-Einstein, we identify a special class of BPS M2-branes that preserve both real internal Killing spinors and correspond to invariant three-dimensional submanifolds inheriting a Sasakian structure. We analyse the quadratic fluctuations around BPS M2-brane instantons in these backgrounds. For the special class of M2-branes in Sasaki-Einstein manifolds, the fluctuation problem reduces to transversely elliptic complexes, and the one-loop partition function can be expressed in terms of the corresponding equivariant indices. We then apply the index formula for the one-loop partition function to invariant M2-branes in $S^{7}/\mathbb{Z}_{k}$, recovering the known result for the $S^3/\mathbb{Z}_k$ instantons and discussing more general invariant BPS cycles. As a further application, we consider M2-brane instantons with $S^3$-quotient worldvolumes in the $(p,q)$-model geometry.

Understanding the time evolution of physical observables in open quantum many-body systems coupled to external environments is a natural and difficult problem, and exact results are still rare. In this work, we study this problem for Wess-Zumino-Witten (WZW) conformal field theories with Lindblad jump operators linear in Kac-Moody current modes. We investigate the exact operator dynamics generated by these Lindbladians, identifying classes of current operators whose Heisenberg equations close and can therefore be solved analytically using the underlying current algebra. In Abelian $U(1)_k$ WZW theories, this closure of operator dynamics holds for arbitrary settings of jump rates and includes exactly tractable cooling dynamics. In contrast, for non-Abelian WZW theories, exact closure occurs only for symmetric current-mode dissipation, where upward and downward current-mode transitions occur with equal rates, and even then it leads to a simple closed evolution only for a single current operator. Generic imbalances, including those needed for cooling, produce additional non-Abelian terms and prevent closure of the opeartor dynamics. Consequently, the current algebra gives rise to a broad family of exactly solvable dissipative dynamics in the Abelian setting, whereas in the non-Abelian case it singles out only a special exactly solvable dynamics corresponding to an infinite-temperature bath.

We investigate the fully relativistic dynamical tidal response of neutron stars up to second order in the frequency. Combining the worldline effective field theory for extended gravitating bodies with perturbation theory of relativistic stellar models, we derive the tidal deformation induced by an external time-dependent field, including a universal logarithmic running term. In the effective theory, we work in dimensional regularization and, through a consistent matching procedure, obtain for the first time the complete leading-order dynamical tidal corrections to both the conservative dynamics and the gravitational-wave signal of compact binaries, including the scheme-dependent finite terms in addition to the running. We show that, in the relativistic regime, dynamical effects cannot be fully captured by mode excitations alone. The magnitude of the additional contribution depends on the stellar compactness, the equation of state, and the running term. Dynamical Love numbers are significantly enhanced with respect to their static counterparts for relatively small compactness. As a result, although they formally enter the gravitational-wave phase at eighth post-Newtonian order, dynamical tidal effects yield a nonnegligible contribution during the late inspiral. Using a Fisher-matrix analysis, we show that third-generation detectors such as the Einstein Telescope could measure dynamical Love numbers for a range of neutron-star masses and equations of state. Conversely, neglecting these effects can lead to significant biases in the inference of static Love numbers, and hence on the nuclear equation of state. Our results highlight the importance of dynamical tidal effects for high-precision gravitational-wave modeling with future detectors.

We provide a moduli stabilisation mechanism for realising anisotropic string compactifications with one or two large extra dimensions, corresponding to the ADD and Dark Dimension scenarios. This is achieved within the type IIB Large Volume Scenario, where an exponentially large Calabi-Yau volume in string units can naturally generate a parametrically low Kaluza-Klein scale. Anisotropy is realised by considering a Calabi-Yau threefold which is a K3 fibration over a $\mathbb{P}^1$ base. The volume of the 4D K3 fibre is stabilised at relatively small values by perturbative corrections to the effective action, in particular string loops and higher-derivative effects, leaving an exponentially large volume of the 2D $\mathbb{P}^1$ base. We argue that complex structure moduli stabilisation can dynamically deform the $\mathbb{P}^1$ base, corresponding to a Tyurin degeneration limit where the internal geometry effectively develops a single large 1D cycle. Within a unified description, the ADD case is instead recovered as a symmetric alternative limit. The potential can feature either a dS vacuum or a quintessence runaway, although in both cases some degree of tuning is required to match the observed cosmological constant scale. We also present an explicit Calabi-Yau orientifold example with consistent brane setup, tadpole cancellation and moduli stabilisation. We analyse the resulting moduli spectrum and associated phenomenological constraints, including supersymmetry breaking, cosmological moduli overproduction and fifth force bounds.

We investigate top-down derivations of effective descriptions of radiation from gravitational systems such as binaries. With a specified cutoff prescription, one can derive worldline effective field theories, although the associated cutoff dependence can complicate their formulation. We investigate a related effective approach in which the dynamics are parameterized by an action on the boundaries of cavities enclosing the individual bodies. We present examples of such cavity descriptions for black holes and for simple models of modified black hole behavior. We also show how cavity effective descriptions connect to observable quantities including Love numbers and aspects of wave profiles, such as the accumulated phase shifts of emitted signals. A primary motivation is to develop a systematic framework for inferring how modifications of classical black hole behavior affect gravitational wave signals. Such modifications may be motivated by the requirement that black hole evolution be consistent with quantum mechanics, or by other models of new black hole behavior. Accumulated phase shifts during the inspiral have, in particular, been argued to provide sensitivity to small new effects. To illustrate the basic principles clearly, we develop the formalism primarily for scalar radiation, and then outline its extension to gravitational perturbations.

We present for the first time the exact solution of the massive Klein-Gordon equation in the Schwarzschild star (perfect-fluid, uniform-density, spherically-symmetric star), including the non-minimal curvature-scalar coupling. The solution is expressed in terms of the general Heun function. A geometry-induced algebraic coordinate transformation reveals a hidden Fuchsian structure that underlies the exact solvability. Known leading- and next-to-leading-order results are recovered in the low-compactness limit. In the Buchdahl limit, we derive a regularity condition for static modes and describe analytically the divergence in amplitude and oscillation wave vector of dynamic modes as they approach the pressure singularity at the center of the star.

Theories with massless particles contain $t$-channel (forward scattering) singularities that cause standard fixed order expressions for partial-wave amplitudes to be ill-defined. This presents an obstruction to systematically improvable partial-wave unitarity bounds. In this work, we study the construction of partial-wave amplitudes in a modified perturbation theory that incorporates long-range interactions focusing on the role of off-shell Coulomb modes. We find that there exists a universal description of the forward scattering region that renders the amplitudes renormalization scale independent. The resulting partial-wave amplitudes become well defined single-scale objects without spurious dependence on the infrared regulator, and we present a practical method for their computation order-by-order in perturbation theory.

We present an algorithm that constructs Calabi-Yau threefold orientifolds and their $F$-theory uplifts to elliptically-fibered Calabi-Yau fourfolds, embedded in toric varieties at codimension one and two respectively. The resulting Calabi-Yau fourfolds arise from triangulations of $6d$ reflexive polytopes – which our method constructs from orientifold data – and are smooth away from isolated terminal singularities. For many of our fourfolds, the construction of the mirror manifold is immediate, enabling the computation of fourfold periods, and thus the seven-brane superpotential. We present multiple examples that demonstrate these capabilities. Our algorithms work with $\mathtt{CYTools}$ and are available through a GitHub repository.

We initiate a bootstrap study of pion form factors in large $N$ QCD. We consider the mixed system of the vector-current two-point function, the pion vector form factor, and the pion scattering amplitude in the chiral limit. At large $N$ these observables are meromorphic, with spectral data constrained by unitarity, crossing symmetry, and Regge boundedness. We obtain bounds of two kinds. The first are rigorous and universal: from analyticity, unitarity and the asymptotic Brodsky-Farrar scaling, we constrain low-energy form-factor coefficients. The second are more phenomenological, of the Shifman-Vainshtein-Zakharov type: feeding in the perturbative ultraviolet behavior at a finite scale lets us bound the pion decay constant, convert a large $N$ lattice measurement into a lower bound on the scale at which asymptotic freedom sets in, and constrain the pion charge radius. Combining these inputs, the space of allowed chiral Lagrangians shrinks toward the region where large $N$ QCD is expected to sit. Our results illustrate how local gauge-invariant probes provide a canonical bridge between the hadronic bootstrap and the microscopic QCD Lagrangian.

In equilibrium, topological defect insertions in quantum and classical partition functions provide non-perturbative probes of phase transitions beyond local observables. In non-equilibrium physics, the spectral form factor provides a minimal probe of universal quantum dynamics and is a product of two partition functions at imaginary inverse temperature. We define the topological spectral form factor (TopSFF) by inserting topological defects acting non-trivially on the doubled partition functions, producing mismatched spacetime world-sheet topologies. For the minimal $\mathbb{Z}_2$ spatially extended defect, implemented by the global swap operator, we derive an exact mapping of the TopSFF of a generic $(1+1)$D many-body chaotic system to an emergent $(3+1)$D non-Hermitian single-particle problem describing a temporal domain wall (tDW). We show analytically that the effective tDW dynamics undergoes a $\mathcal{PT}$ symmetry breaking transition at finite interaction strength $\epsilon_{\mathrm{EP}}$: below $\epsilon_{\mathrm{EP}}$, the leading modes are polarized into Gaussian or non-Gaussian tDW sectors and the TopSFF varies monotonically and exponentially with system size; above $\epsilon_{\mathrm{EP}}$, the tDW sectors hybridize and the TopSFF oscillates with system size; at the exceptional point $\epsilon_{\mathrm{EP}}$, Jordan non-diagonality produces a linear-in-system-size enhancement. For minimal temporally extended topological defects associated with global swap or time translation, we derive exact universal TopSFF free-energy scaling forms. We verify these signatures numerically in independent models.

We compute spectral functions of graviton modes in Lorentzian quantum gravity, interpolating between classical general relativity and an asymptotically safe ultraviolet fixed point. Using functional renormalisation adapted for theories in Lorentzian signature, and enhanced by new symmetry conditions to account for underlying Ward identities, we derive and solve flow equations directly for the Källén-Lehmann representation of propagators. Consistent results are found for several sets of renormalisation conditions yielding normalisable spectral functions for the graviton and the scalar graviton mode, in agreement with effective theory in the infrared. We further calculate the full quantum effective action to quadratic order in curvature, extract graviton-induced form factors, and discuss implications for unitarity of quantum gravity.

Ultra-cold bosonic systems can be tuned to exhibit quantum phase transitions. For example, the Rabi-coupled bosonic system exhibits ferromagnetic and paramagnetic phases, whereas the spin-orbit-coupled system exhibits exciting phases such as supersolidity. The physics of these phases and phase transitions is very rich. It is an important topic of research to probe these phases and phase transitions using various tools in many-body physics. The operator product expansion (OPE) provides one such tool. It expresses the product of two separated operators as a series expansion of local operators. In this article, we will derive the OPE of two operators $\psi^\dagger_\sigma(\vec r)$ and $\psi_{\sigma'}(\vec r')$. More specifically, we look for the contact density term, which controls many of the universal physics of the underlying bosonic system.

We study the closed-string tree-level Virasoro-Shapiro amplitude using the dispersive S-matrix bootstrap. For the ten-dimensional maximally supersymmetric four-point amplitude, we impose analyticity, crossing symmetry, partial-wave unitarity, and Regge boundedness. With the massless graviton pole kept explicitly, the resulting dispersion relations and crossing null constraints give numerical bounds on the leading low-energy coefficients normalized by the gravitational coupling, which agree at our reference truncation with earlier bounds from smeared functionals. We then introduce a Virasoro-inspired ansatz, which becomes a set of nonlinear relations among Wilson coefficients and shrinks the allowed region toward the Virasoro-Shapiro trajectory. Finally, we study a gravity-pole-subtracted setup, where the regular part of the amplitude has a well-defined forward limit. In this stripped problem, on the slice of the leading coefficient that places the first Virasoro-Shapiro resonance at the cutoff, the nonlinear constraints reduce the allowed region to a small island containing the Virasoro-Shapiro point. Imposing the ansatz to all orders, together with spectral assumptions not used in the numerics, we show analytically that positivity selects the Virasoro-Shapiro spectrum.

Black hole static Love numbers vanish, but their dynamical counterparts do not. We present the entire scheme-independent dynamical response $\bar{F}_{\ell,s}$ of a Schwarzschild black hole in closed form, to all orders in $G$, and for every spin $s$ and multipole $\ell$. The result is $\bar{F}_{\ell,s}/4\pi R_S^{2\ell+1}=\Phi_{\ell,s}(\bar{y})-\tfrac12\eta\,\Phi_{\ell,s}'(\bar{y})$ with $\bar{y}=-\tfrac12\eta^2\tau$ and $\eta=i\omega R_S$. Here $\Phi_{\ell,s}$ is simply the leading-log solution to the renormalization group equation, but lifting the running logarithm to $\tau=\log(R_S/R)-2\sum_{k\ge2}\zeta_k\,\eta^{k-1}$ resums it to all orders. This tower of Riemann zeta values is the Newtonian phase in disguise: it originates from the long-range $1/r$ potential in the far zone, and is universal across multipole and spin. Our result exhibits a factorization pinned to three ingredients: the strength of horizon absorption, the running rate in the near zone, and the dressed logarithm in the far zone. Using shell effective field theory, we independently verify our formula for scalar, electromagnetic, and gravitational perturbations, reaching $\mathcal O(G^{15})$.