Papers

Multi-dimensional chaos refers to processes described by erratic functions of several dynamical variables. In this letter we analyze the string scattering amplitudes of highly-excited states and ground states. We show that the amplitudes, which depend on a scattering angle and a polarization angle, are characterized by two sets of non-intersecting curves associated with the vanishing of the derivatives with respect to the angles. We introduce the notion of the "area eigenvalue" $A_n$ associated with the $n$-th curve. We compute the spacings $\delta_{n}= A_{n+1}-A_n$ and their ratios $r_{n}=\frac{\delta_{n+1}}{\delta_n}$. We show that the distributions of the spacing ratios take the form of the RMT Gaussian $\beta$-ensembles. The curves associated with the scattering angle tend to converge to the Gaussian Orthogonal Ensemble value of $\beta=1$ and those related to the polarization angle to the Gaussian Unitary Ensemble $\beta=2$. We also compute the “areas form factor" associated with the areas and discover the regions of decline, ramp and plateau which characterize chaotic processes. The slope of the ramp seems to agree with the $\beta$ values extracted from the distribution of the spacing ratios.

The standard approach to canonical quantization encounters difficulties in dealing with perturbations that alter the kinetic structure of unperturbed theories. We show that the covariant phase space formalism provides a natural and technically efficient way to circumvent this obstruction. We illustrate the method with an exactly solvable model: a two-dimensional non-relativistic charged particle moving in a magnetic field and a harmonic confining potential, with its kinetic energy viewed as a perturbation. We quantize this model with covariant phase space formalism by constructing the solution perturbatively. We then calculate the energy spectrum and the unequal-time commutators of this model, and obtain the results that agree with the expansion of the exact theory. The procedure developed here is intended to serve as a systematic framework for the canonical quantization of more complex effective theories with higher-derivative or velocity-dependent perturbations.

In this Letter, we confront high-order post-Minkowskian (PM) predictions for generic-spin black-hole scattering with numerical-relativity (NR) simulations for the first time, targeting improvements for eccentric and precessing waveform modelling. We extract azimuthal and polar scattering angles from NR and relate them to the PM spin-kick observable. We introduce asymptotic Euler angles for unbound motion and derive their geometric relation to the scattering angles. Notably, NR exposes a strong-field precessional turning-point structure, including a polar-angle sign change absent in the perturbative PM results.

Parallel to the AdS Scale-Space construction in \cite{Yang:2026AdS}, we construct a $(d+1)$-dimensional de Sitter QFT directly from $d$-dimensional CFT data. The dS geometry is the moduli space of oriented balls, and the lifted operators are obtained by conformal-family Casimir completion. The Euclidean parent CFT gives a finite wavefunction representation, while the Minkowski parent CFT gives a double-time representation that is unitary in the parent-time polarization. This provides a CFT-based framework for reexamining puzzles of dS and double-time QFT.

In this paper, we generalize our previous work on the fermionic tidal Love numbers (TLNs) to higher-dimensional Reissner-Nordström black holes. The massless Dirac equation is solved in $D$-dimensional spacetime using ingoing Eddington coordinates and regular tetrads. After identifying the regular solution branch, we extract the fermionic TLNs from its asymptotic behavior at infinity. The resulting TLNs exhibit a rich dimension-dependent structure that generalizes the four-dimensional case. Unlike bosonic TLNs, which vanish for certain values of the total angular momentum $l$ in dimensions $D>4$, fermionic TLNs remain non-zero for all $l$ and $D \geq 4$, except for extremal black holes. Moreover, the $l$-dependence weakens as $D$ increases, disappearing entirely in the infinite-dimensional limit. These results provide new insights into black hole responses to fermionic perturbations in higher-dimensional spacetimes.

We show analytically that a dark-scalar order can preempt anisotropic $p$-wave order in a gauge-invariant hyperscaling-violating holographic superfluid. Near a point where the source-free scalar and vector instabilities coincide, their grand potentials can be compared in closed perturbative form. The homogeneous coexistence solution is a saddle, so the preferred state switches directly between an anisotropic vector-condensed phase and an isotropic dark-scalar-condensed phase across a first-order line. This line terminates where the two continuous normal-to-ordered transition lines meet, producing a bicritical point. A joint hyperscaling/Lifshitz deformation moves the bicritical point, showing that nonconformal infrared scaling can select not only the critical scale but also the spatial symmetry of the preferred ordered phase. Entanglement anisotropy follows the squared vector order and therefore jumps to zero when the dark-scalar phase takes over. The result provides a controlled dark-Higgs-motivated mechanism for ordered-to-ordered symmetry conversion at leading probe order.

We review some selected topics in celestial holography on the search for a celestial dual to quantum gravity in flat spacetimes. We focus on the bottom-up approach, emphasizing symmetries, key elements in celestial CFT, interplay with twistor theory, and connection to AdS/CFT.

We revisit equivariant Gromov-Witten theories on P1 and on P1 x C2. One can introduce three equivariant parameters associated to rotations of the sphere as well as the two planes. A number of points in the parameter space have known holographic duals. These include the symmetric orbifold point dual to the AdS3 x S3 x C2 string theory at string scale radius of curvature, the grand canonical Hurwitz theory and the product of two Kontsevich models. Within this framework, we discuss interpolations in the equivariant parameters. Firstly, we move between the small and large equivariant parameter regimes in Gromov-Witten theory on P1. At large equivariant parameter, the model is dominated by the pure topological gravity theories at the two fixed points while at small equivariant parameter the theory is equivalent to the grand canonical Hurwitz theory. The deformation is a solvable analogue for the interpolation in the transposition coupling in the moduli space of the AdS3/CFT2 duality. Moreover, we propose that the full equivariant correspondence between the Gromov-Witten theory on P1 x C2 and the symmetric orbifold of the equivariant plane can be embedded in string theory. On the boundary side of that correspondence, we analyze the scaling limit from the equivariant to the ordinary cohomology ring for the Hilbert scheme of points on the plane in terms of Jack symmetric polynomials. We explicitly compute several structure constants of the equivariant cohomology ring and point out their intriguing positivity and integrality properties.

Recent DESI baryon acoustic oscillation measurements, especially when combined with Type Ia supernova and CMB data, sharpen the case for possible low-redshift dynamics in the dark energy sector. We study a simple and physically transparent realization of such dynamics: a quintessence field that is Hubble frozen for most of cosmic history and starts to oscillate around its minimum recently (at a redshift $z\approx 0.1$). This late onset of oscillations can occur in a broad class of models where the quintessence potentials have a shallow slope away from the minimum and steepen near it. This class of models can improve the fit relative to $\Lambda$CDM, with $\Delta\chi^2\simeq -9$, while remaining competitive with common phenomenological dark energy parameterizations with the same number of parameters. The preference is driven mainly by the background expansion history, and near the best-fit region the resonant growth of quintessence perturbations and the associated Integrated Sachs-Wolfe (ISW) contribution remain small. More precise low-redshift distance measurements, together with late-time probes such as the ISW effect and lensing, may help distinguish this oscillating quintessence scenario from other forms of late-time dark energy dynamics.

The fermion $\nu_{R}+C\overline{\nu_{R}}^{T}$ in the Type I seesaw model, which is commonly identified with a Majorana fermion, is not a Majorana fermion and thus no neutrinoless double $\beta$ decay (in the case of $\nu_{L}-C\overline{\nu_{L}}^{T}$) due to the theorem on the absence of a Majorana-Weyl fermion in the representation of the d=4 Lorentz group. The leptogenesis is realized naturally together with the neutrinoless double $\beta$ decay by first defining C symmetric Majorana fermions in the seesaw model using an analogue of the Bogoliubov transformation.

The chance observation of a novel traveling wave in a canal led over time to the formulation of a nonlinear wave equation – the Korteweg–de Vries equation – that describes strikingly robust disturbances now called solitons. The figure of an isolated soliton corresponds to a reflectionless potential that supports a single bound state in the one-dimensional Schrödinger equation. An appropriate combination of individual solitons yields a symmetric reflectionless potential that supports multiple bound states. Thus, the KdV equation opens the path to solving the inverse scattering problem for a collection of bound states. Applied to the quarkonium spectra, this formalism allows the construction of reflectionless approximations to the confining potentials that account for the force between quarks, and to tests of the flavor-independence of the interquark interaction.

We present a conjecture for the double-logarithmic behavior of the logarithm of large-charge correlators in the null limit in planar $\mathcal{N} = 4$ Super-Yang-Mills theory. Generalizing earlier results for four- and five-point functions, our proposal predicts this behavior for $n$-point functions to all loops in terms of the tilted cusp anomalous dimension. In the dual amplitude description, it reproduces the expected small-mass behavior of massive amplitudes in the equal-mass limit.

Lagrangian formulations for (ir)reducible integer higher-spin massless and massive Poincare group representations subject to Young tableau with $k$ columns $Y[\hat{s}_1,\hat{s}_2,...,\hat{s}_k]$ in $d$-dimensional Minkowski space-time are firstly presented. The particles are described in a metric-like formulation by tensor fields with $k$ groups of antisymmetric Lorentz indices $\Phi_{\mu^1[{\hat{s}_1}],\mu^2[{\hat{s}_2}],..., \mu^k[{\hat{s}_k}]}$ by means of the BRST procedure with complete, $Q$, and incomplete, $Q_c$, BRST operators. Starting from a description of bosonic mixed-antisymmetric higher-spin fields in terms of an auxiliary Fock space associated with a special Poincare module, we realize a conversion of the initial operator constraint system into a system of first-class operator constraints. To this aim, we find, in first time, by means of Verma module the auxiliary representations of the constraint subalgebra, to be isomorphic due to Howe duality to $so(k,k)$ algebra, and containing the subsystem of second-class operators in terms of new oscillator variables forming the Fock module. An unconstrained (with $Q$) and constrained (with $Q_c$ and BRST invariant algebraic constraints) gauge Lagrangian formulations with equivalent dynamics, but different configuration spaces are found. Concept of consistent cubic interactions are suggested within $Q_c$-approach together with evaluation the structure of cubic vertices for $p\geq 3$ copies of interacting massless mixed-antisymmetric fields.

Motivated by the recent results reported by the CONUS$+$ collaboration, in which coherent elastic neutrino-nucleus scattering (CE$\nu \mathcal{N}$S) with reactor antineutrinos was observed for the first time, we perform a statistical analysis to constrain possible low-energy scenarios of physics beyond the Standard Model (BSM). The models considered include effective non-standard vector interactions of neutrinos (NSI) and generalized neutrino interactions (NGI) by light vector bosons within the $E_6$ and $U(1)_{L_e-L_\mu}$ frameworks.

Microscopic models of neutron-star matter have been widely used in astrophysical applications. The focus of attention has been on densities up to the maximal densities reached in stable neutron stars. The possibility that the underlying model assumptions may have important implications at higher densities has not been addressed. Here, we show that the behaviour at higher densities is strongly constrained by requiring a causal, stable, and thermodynamically consistent extension to the perturbative-QCD regime. We explicitly reveal what that behaviour must be and provide a tool for constructing and visualizing such extensions. We find that purely hadronic models trusted up to the maximal central density often require radically different behaviour at higher densities from that assumed in the original model, while models with additional degrees of freedom fare better. Our analysis disfavors purely nucleonic models for describing all stable neutron stars and supports the appearance of some type of additional degrees of freedom in stable massive neutron stars.

This note is a detailed review of the geometry behind the Batalin-Vilkovisky formalism and how it fits into the framework of the quantum odd symplectic category and the odd quantization functor.

We use supersymmetric localization for 3d $\mathcal{N}=2$ SCFTs to establish a relation of the form $Z_{S^1\times S^2} \sim |Z_{S^3_b}|^2$ between the superconformal index and the squashed three-sphere partition function. This applies to general 3d $\mathcal{N}=2$ SCFTs and is derived using a saddle point approximation in the Cardy-like limit of small $S^1$ radius and large squashing parameter. We also show a similar relation between the topologically twisted index and the squashed sphere partition function. In the context of holography our results lead to a relation between the partition function of supersymmetric asymptotically AdS$_4$ black holes and the large $N$ limit of $Z_{S^3_b}$ of the dual SCFT. Since the latter encodes the gauged supergravity prepotential, this result is akin to the OSV conjecture for asymptotically flat black holes. We confirm this relation in detail for SCFTs arising from M2-branes using recent large $N$ results from supersymmetric localization. We also briefly discuss a similar relation for 5d SCFTs and its implications for asymptotically AdS$_6$ black holes.

We study energy extraction via magnetic reconnection from a rotating dyonic black hole in four-dimensional $N=2$, $U(1)^2$ gauged supergravity. Using the Comisso-Asenjo mechanism in the ZAMO frame, we derive the asymptotic hydrodynamic energy per unit enthalpy $\epsilon^{\infty}_\pm$ and determine when reconnection outflows attained negative energy at infinity. By varying the spin $a$, electric and magnetic charges, NUT parameter $N_g$, and gauge coupling $g$, we compute the cutoff magnetization $\sigma_0^{\rm cutoff}$ and map the region of parameter space that admits $\epsilon^{\infty}_-<0$. We find that $\sigma_0^{\rm cutoff}$ and the very existence of Comisso-Asenjo extraction are tightly controlled by $g$ and the dyonic charges: increasing $g$ or pushing the charges toward extremality raises $\sigma_0^{\rm cutoff}$ and shrinks the CA-active part of the ergoregion. Unlike Kerr, the spin enters through the normalization factor $\Xi$, and the quartic horizon function $\Delta_g$, so geometric effects from the AdS/NUT deformation dominate the usual frame-dragging enhancement. As a result, the extracted power and efficiency are non-monotonic in $a$ and peak at intermediate spin ($a\sim0.8$); near-extremal rotation is not required for high efficiency, provided $g$ is small and $Q$ is moderate. Efficient extraction further demands extreme magnetization and nearly radial outflows, confining the active reconnection layer to a thin shell, just outside the horizon.

Primordial Black Holes (PBHs) have emerged as a leading non-particulate candidate for dark matter and a unique cosmological probe, a paradigm shift accelerated by the detection of anomalous binary mergers by the LIGO-Virgo-KAGRA (LVK) collaboration. While the literature is rich with phenomenological constraints, the fundamental quantum and relativistic underpinnings governing PBH genesis and evolution often receive comparatively less emphasis. This review aims to bridge that gap by systematically detailing the physics of PBH formation and their subsequent evolutionary trajectory. We critically examine the hydrodynamic complexity of the early universe, establishing the relativistic thresholds for collapse, the non-linear race against sound in the primordial plasma, and the rigorous mathematical utility of the compaction function. Furthermore, by incorporating the dynamic nature of FLRW backgrounds, higher curvature corrections, and quantum backreaction via the memory burden effect, we challenge the standard hawking evaporation and show that extreme-curvature environments halt evaporation entirely, leaving Planck-scale relics that evade current extragalactic bounds. Finally, we map the multimessenger observational landscape, highlighting how the imminent search for sub-solar mass inspirals by next-generation gravitational wave observatories such as the Einstein Telescope and Cosmic Explorer could yield smoking-gun evidence for the PBH paradigm, ultimately transforming these primordial relics into unparalleled laboratories for high-energy physics.

We introduce a variational method to solve continuum quantum models with discrete tensor network techniques. The method leverages wavelet matrix product states (wMPS): matrix product states built on top of sufficiently regular ($N\geq 6$) Daubechies scaling functions. These states live in the continuum field theory Fock space, have finite energy density, and can be optimized with standard algorithms, without restriction to free theories. Further, exploiting the multi-resolution analysis built into wavelets, and its quantum circuit description, we can iteratively refine wMPS to obtain accurate approximations at arbitrarily fine length-scales. We showcase the efficiency of the method on the Lieb-Liniger model, computing energy density and correlation functions.

Given a quiver Q with gauge dimension $\bf v$ and framing dimension $\bf w$, one can define the extended quiver variety $\widetilde{\mathcal M}(\mathbf v,\mathbf w)$, which is a smooth family of deformations of the Nakajima quiver variety $\mathcal M(\mathbf v,\mathbf w)$. In this paper we discuss two vertex algebras which chiralize the geometry $\widetilde{\mathcal M}(\mathbf v,\mathbf w)$. We construct a sheaf of $\hbar$-adic vertex superalgebras $\mathscr D^{\mathrm{ch}}_{\widetilde{\mathcal M}(\mathbf v,\mathbf w),\hbar}$ on $\widetilde{\mathcal M}(\mathbf v,\mathbf w)$ which quantizes the jet bundle of $\widetilde{\mathcal M}(\mathbf v,\mathbf w)$, and define a vertex algebra $\mathsf D^{\mathrm{ch}}(\widetilde{\mathcal M}(\mathbf v,\mathbf w))$ to be the $\hbar=1$ specialization of the $\mathbb C^{\times}$-finite part of the vector space of global sections $\Gamma(\widetilde{\mathcal M}(\mathbf v,\mathbf w), \mathscr D^{\mathrm{ch}}_{\widetilde{\mathcal M}(\mathbf v,\mathbf w),\hbar})$. We define another vertex superalgebra $\mathcal V(\mathbf v,\mathbf w)$ by BRST reduction of the tensor product of the $\beta\gamma bc$-system and Heisenberg VOA associated to the quiver Q, and show that there exists a natural vertex superalgebra map from $\mathcal V(\mathbf v,\mathbf w)$ to $\mathsf D^{\mathrm{ch}}(\widetilde{\mathcal M}(\mathbf v,\mathbf w))$. Under certain technical assumptions, we prove that the negative degree BRST cohomologies of the tensor product of $\beta\gamma bc$-systems and Heisenberg VOA associated to the quiver Q are zero, and under stronger assumptions, that the aforementioned vertex superalgebra map is injective. Physically, the vertex superalgebra $\mathcal V(\mathbf v,\mathbf w)$ is closely related to the boundary VOA of the H-twisted 3D $\mathcal N=4$ quiver gauge theory associated to the quiver Q with gauge and framing dimension vectors $\bf v$ and $\bf w$.

We initiate the study of nonlinear effects in the ringdown phase of black hole mergers using the effective theory of black hole dynamics in the large-D limit. This framework offers several advantages: the quasinormal mode spectrum, including nonlinear corrections, is analytically tractable; numerical simulations of collisions are computationally inexpensive; and the extraction and analysis of the ringdown signal are clean and controlled. As a proof of concept, we derive analytic expressions for the third-order response of a static black hole driven by a single quasinormal mode, and apply them to study the ringdown following head-on collisions of non-spinning black holes across a range of velocities and mass ratios. We find that including nonlinear effects, up to quadratic and cubic order, improves the accuracy of quasinormal-mode modelling of black hole relaxation by several orders of magnitude. The results also show a clear growth in the strength of nonlinear effects as the collision velocity increases.

Chaotic dynamics play an important role in a number of physical systems. One of the qualitative hallmarks of this behavior is the appearance of a sufficiently "complex" spectrum of energy levels. This also makes it challenging to directly verify the onset of chaos in interacting quantum field theories. We present a class of 4D superconformal field theories (SCFTs) given by orbifolds of 4D $\mathcal{N} = 4$ Super Yang–Mills theory in which operator mixing in a controlled subsector is described by an effective spin chain in one spatial dimension with nearest neighbor interactions tuned by the marginal couplings of the SCFT. Tuning the marginal couplings results in a chaotic spectrum, while generically the spin chain exhibits Anderson localization. We diagnose the onset of chaos by analyzing the statistical distribution of eigenvalues of the dilatation operator, in particular properties such as eigenvalue level repulsion, spectral rigidity, and the spectral form factor. We also show that other diagnostics such as Krylov complexity sometimes do not faithfully capture this information. This structure defines a chaotic billiard in the target space of the stringy realization. We also comment on the large $N$ holographic dual description, where the controlled single spin chain approximation must be supplemented by multi-trace dynamics, i.e., the splitting and joining of multiple spin chains.

We obtain a new electrovacuum and related background spacetime that contains a cosmological horizon supported entirely by the electromagnetic field. The local solution belongs to the Kundt class of type D, but is globally distinct. We further construct an exact solution describing a Kerr black hole in this background. We study the global structure including horizons and singularities, and derive the first law of black hole thermodynamics. The emergence of a cosmological horizon in Einstein-Maxwell gravity without invoking a positive cosmological constant or dark energy is tantalizing, and may provide a new avenue for exploring cosmological and astrophysical phenomena related to black holes and the late-time cosmology.

In arXiv:2604.19861, we gave general criteria for when one zero-mean Gaussian state can be excited out of another in a (generalized) free field theory. Here we extend this analysis to the case of nonzero mean, i.e., to Gaussian states with vacuum expectation values (VEVs). We prove that excitability is possible exactly when (i) the connected two-point functions satisfy criteria like those in arXiv:2604.19861, and (ii) the difference of the VEVs is bounded relative to the two-point functions. As an application, we give an explicit computation showing that in anti-de Sitter spacetime, a VEV shift can be excited from the Klein-Gordon vacuum if and only if its boundary extrapolation can be excited from the vacuum of the dual conformal field theory.

We show that tree-level gluon correlators in $\textrm{AdS}_4$ admit a natural expansion in terms of flat-space scattering amplitudes at all multiplicities. In particular, every $n$-point correlator can be decomposed into a sum over energy poles whose residues are flat-space amplitudes. The $n$-point amplitude encodes the flat-space limit while curvature corrections are captured by lower-point amplitudes with merged external data. The merging of external polarizations is recursively defined via an AdS analogue of the Berends-Giele currents, giving rise to all-multiplicity formulae which we verify against Feynman diagram computations up to five points. Crucially, our approach works at the level of full correlators rather than individual diagrams, providing an elegant and transparent alternative to conventional approaches for computing correlators in anti-de Sitter space.

While QCD fundamentally only depends on the values of the strong coupling and the quark masses, it exhibits a rich nonperturbative structure at low energies, where composite fields emerge as the relevant degrees of freedom. In this work, we present a first-principles framework that captures the transition from fundamental QCD to its low-energy sector in vacuum. It builds on the dynamical hadronization technique within the functional renormalization group approach to two-flavor QCD. In this framework, the low-energy constants relevant for effective models, including effective masses and coupling strengths, naturally emerge from the underlying renormalization group flow without introducing free parameters beyond those of QCD itself. We investigate the dynamical emergence of the pion, the $\sigma$-meson and the scalar diquark in both imaginary and real time, and determine a set of QCD low-energy constants which can be used to fix the free parameters of models of dense quark matter with a two-flavor color superconducting phase. In particular, this includes previously unknown properties of the scalar diquark. Our results provide important microscopic input for constraining color superconducting phases, which are expected to play a key role in our understanding of dense neutron star matter.

Andreev reflection converts an incoming fermion into an outgoing hole and is usually tied to a superconducting interface. We show that an analogous charge-conjugating boundary condition emerges from a purely lattice duality defect. Starting from a Majorana representation of the transverse-field Ising chain, we construct a folded lattice model in which a boundary Majorana impurity implements a one-site translation of a staggered Majorana chain. In the continuum, this translation becomes a chiral fermion-parity defect: it flips the sign of the only left-moving Majorana mode while leaving the right-moving mode unchanged. When the two Majorana modes are recombined into a complex fermion in the folded geometry, this sign flip becomes the Andreev-like boundary condition. Our lattice formulation gives a microscopic interpretation of the Emery–Kivelson boundary of the two-channel Kondo problem and of Maldacena–Ludwig monopole scattering, while identifying the boundary as the interface between a Kitaev-chain SPT phase and a gapless chain. The same Majorana translation defect also provides a lattice realization of an axial $U(1)_A$-symmetric charge-flip boundary.

We investigate test particle dynamics and gravitational wave (GW) phenomenology in an exact spherically symmetric vacuum solution of Freund - Nambu scalar - tensor gravity. This framework generalizes the Janis - Newman - Winicour (JNW) naked singularity via a geometric non - linear coupling $q$ and a direct scalar - particle coupling $g_s$. We demonstrate that these parameters systematically modify the Innermost Stable Circular Orbit (ISCO) - which shifts inward for $g_s > 0$ - and the Marginally Bound Orbit (MBO). Furthermore, we classify bound periodic trajectories to isolate extreme zoom - whirl orbits exhibiting intense periapsis precession. By applying the Numerical Kludge method to Extreme Mass - Ratio Inspirals (EMRIs), we reveal that scalar - tensor corrections induce a macroscopic temporal dephasing in high - frequency GW bursts, even when the orbit's spatial topology is preserved. These unique phase shifts offer a robust diagnostic signature for future space-based observatories like LISA to probe the strong - field regime and constrain scalar - tensor extensions of general relativity.

We study the decay of de Sitter space into black holes whose horizons are dressed by BPS vortices of a $\mathrm{CP}^1$ action. The process is mediated by a regular Euclidean instanton obtained as a vortex-deformed generalization of the Nariai instanton. Its Euclidean geometry has the form $S^2\times\Sigma$, where $\Sigma$ is a compact surface whose geometry is shaped by the vortex configuration. The resulting decay rates are controlled by a discrete topological charge, showing that matter vortices open a new topologically organized family of decay channels for de Sitter space.