Papers

The Collins-Soper kernel is calculated from a vacuum soft function using space-like Wilson lines with complex-directional vectors on the Euclidean lattice. Our pure gauge calculations with this method achieve high statistical precision in computing the soft function, whose rapidity dependence is well described by Collins-Soper evolution across a wide range of rapidity differences. The extracted kernel contains errors comparable to those achieved in state-of-the-art lattice calculations based on hadronic observables, but exhibits saturated behavior at large transverse Wilson-line separations.

We prove a `nonlinear graviton theorem' for higher-spin self-dual gravity. We consider small deformations of the complex structure of the non-projective twistor space that are bounded in a specified region near the origin and investigate the space $M_{HS}$ of holomorphically embedded complex planes $\mathbb{C}^2$ that intersect the origin. We show that this space is an infinite dimensional complex manifold with a canonical projection onto a four-dimensional holomorphic self-dual spacetime $\mathcal{M}$, and discuss the geometry induced on this new higher-spin space. Solutions of higher-spin self-dual gravity are then obtained by choosing an embedding of spacetime $\mathcal{M}$ into higher-spin space $M_{HS}$, with higher-spin symmetries arising from the different choices of embedding. Integrability of the theory is manifested in the form of a Lax pair for the system that we present. We conjecture that chiral higher-spin gravity can similarly be realized by considering deformations that are unconstrained at the origin.

We study the dynamics of first-order inverse phase transitions (driven by superheating) at strong coupling, focusing on the top-down Witten-Sakai-Sugimoto model for holographic QCD. Two cases are considered: the deconfinement transition in the unflavored version of the model and a chiral symmetry-restoring transition occurring in the deconfined phase of the full theory. In both cases, we imagine driving the system into a metastable phase at high temperature, inducing the nucleation of bubbles of the stable phase. For both classes of transitions, we find the corresponding Euclidean bounce solutions and compute the bubble nucleation rates and the relevant transition parameters. For the deconfinement transition, the large jump in the number of degrees of freedom between the two phases suggests that the bubble wall velocity is parametrically small; we provide a rough estimate of it near the critical temperature. In the case of the chiral transition, instead, we compute the bubble wall velocity and the friction force exerted on the bubbles employing motivated ansatze and approximations for the steady-state configurations.

By incorporating leading $c\,$-dependent corrections to the Carroll transformations, we introduce the “post-Carroll transformations”. We demonstrate that these transformations are consistent with post-Carrollian mechanics \cite{Najafizadeh:2025ksm}; furthermore, they give rise to the so-called “post-Carroll algebra”. We show that, unlike the Carroll algebra, this new structure allows for a central charge in higher dimensions; we refer to it as the “Carroll-Bargmann algebra”. To construct conformal extensions, we first build the conformal extension of the post-Carroll algebra and study field theories invariant under this symmetry. We then construct the conformal extension of the Carroll-Bargmann algebra, referred to as the “Carroll-Schrödinger algebra”, and demonstrate that it precisely matches the symmetry algebra of the higher-dimensional Carroll-Schrödinger theory \cite{Najafizadeh:2024imn}. Finally, we derive the general form of two-point functions in a post-Carrollian CFT, which in $1+1$ dimensions exhibits both electric and magnetic sectors, while in higher dimensions only the magnetic sector survives.

We show that Gauss-Bonnet black branes in five-dimensional anti-de Sitter gravity are unstable when the Gauss-Bonnet coupling falls outside the range allowed by the conformal collider bounds. The unstable modes and the boundary causality violating modes are connected by a phase rotation of complex boundary momentum.

We study four-point functions of identical scalar operators in conformal field theories with AdS duals in large-$N$ expansion. We analyze the appearance of higher-trace operators in theories dual to bulk $\phi^3$ and $\phi^4$ interactions, focusing on how these operators are required by crossing symmetry. We compute part of the OPE data associated with these operators. We also introduce a diagrammatic framework for organizing the different terms in the conformal block expansion within the large-$N$ expansion. This framework refines the use of crossing symmetry by allowing it to be applied to individual diagrammatic topologies, rather than only to the full correlator. It further separates different contributions to the OPE data by associating them with different cut diagrams. In this language, the emergence of higher-trace operators and their relation to lower-trace OPE data become more manifest.

We study diffractive photo- and electroproduction of the $S$-wave heavy quarkonia $J/\psi$, $\psi(2S)$, and $\Upsilon(nS)$ at energies relevant for the Electron-Ion Collider (EIC). The production amplitude is evaluated while retaining the full transverse-momentum ($\ell_t$) dependence of the hard two-gluon kernel, that is, without expanding the impact-parameter Bessel kernel as is done in the small-size color-dipole limit. The quarkonia light-cone wave functions are built from Cornell-potential solutions of the Schrödinger equation, normalized to the measured leptonic widths, and combined with a modern collinear gluon distribution. After benchmarking the framework against the full set of HERA charmonium cross-section ratio data, we provide a consistent set of bottomonium cross-section ratio predictions in EIC kinematics. We find that the full $\ell_t$-resolved treatment systematically improves the description of the radially excited states relative to the leading dipole limit, and we identify the kinematic windows where this difference is largest.

Color graphs and their subgraphs, referred to as bubble graphs, correspond bijectively to the simplicial complexes of pseudomanifolds and their subsimplices, respectively. In this paper, we introduce matrix representations for colored graphs and their associated bubble graphs. By using this correspondence, we define simplicial-complex matrices and subsimplex matrices that encode the simplicial complexes of pseudomanifolds and their subsimplices. Moreover, we formulate mutation and crossover operations on colored graphs. Through the established correspondence among simplicial complexes, colored graphs, and simplicial-complex matrices, we extend these operations to simplicial complexes and simplicial-complex matrices. We further implement an algorithm generating simplicial-complex matrices and a genetic algorithm performing mutation and crossover of them to produce pseudomanifolds exhibiting diverse topologies. In addition, we implement procedures for decomposing the generated simplicial-complex matrices into simplex matrices, reconstructing the simplicial complexes of the associated pseudomanifolds from this information, and computing geometric quantities such as the volume, circumcenter, and dual-simplex volume of each simplex.

The path integral of the Dirac fermion with vector and axial gauge backgrounds is analyzed near the infrared limit in the presence of residual irrelevant current-current interaction. After integrating out fermions, a semiclassical low-energy effective action is obtained, written in terms of currents. Its expression is found to correspond to the hydrodynamic action previously proposed for perfect barotropic fluids with anomalies at zero temperature. This approach also leads to two further hydrodynamic actions to be associated, respectively, with the Weyl fermion, and the Dirac fermion having independent vector and axial currents. These actions feature four- and five-dimensional bulk-boundary terms, owing to anomaly inflow, which are identified as being the so-called transgression forms. These are generalizations of Chern–Simons forms that involve two gauge fields: the dynamical field and the background field. The path-integral argument provides a “microscopic” explanation for several ingredients of the action formulation of hydrodynamics that are necessary to incorporate anomalies. It also clarifies the infrared reduction required to pass from the effective field theory to a local hydrodynamic description. This reduction is implemented by considering restricted variations of the action, familiar from hydrodynamics, which at the same time lead to four-dimensional equations of motion from the five-dimensional transgression terms.

We show that the moment of the transfer matrix of the double scaled SYK model is written as an overlap between the stationary state of ASEP (asymmetric simple exclusion process) and a product state. We argue that this overlap is an analogue of the strange correlator appearing in the correspondence between the Levin-Wen string-net model and the Turaev-Viro state sum.

Spontaneous CP violation (SCPV) provides an attractive solution to the strong CP problem. However, SCPV after inflation suffers from the formation of CP domain walls, requiring the maximal temperature of the Universe to lie below the CP-breaking scale. In the present work, we then propose a dynamical mechanism that removes this cosmological constraint without introducing permanent explicit CPV. We consider a new scalar field that acquires a large field value with a nontrivial phase in the early Universe and induces a transient bias among degenerate CP vacua through a higher-dimensional interaction with a CP-breaking scalar field. This bias triggers the decay of CP domain walls after they form. As the new scalar field evolves toward the origin, the bias disappears, leaving the low-energy CP structure intact. We derive the conditions for successful domain wall decay and identify the viable parameter space. Furthermore, we point out that the coherent oscillation of the new scalar field naturally survives as dark matter, linking the resolution of the CP domain wall problem to the origin of dark matter.

In this work, we construct the simplest pseudo-Hermitian quantum field theory that is supersymmetric. This is the pseudo-Hermitian Wess-Zumino model in the sense that it contains a pair of symplectic fermions (anti-commuting scalar fields) that satisfy the Klein-Gordon equation and a spin-half boson that satisfies the Dirac equation. The conventional spin-statistics theorem is circumvented through the use of pseudo-Hermitian conjugation to define field adjoints. To make the supersymmetry manifest, we formulate the pseudo-Hermitian Wess-Zumino model using the superfield formalism. These superfields are Grassmann-odd so it is not possible to construct non-vanishing cubic interactions using only these superfields. We show that this problem can be resolved by coupling the pseudo-Hermitian Wess-Zumino model with the Hermitian Wess-Zumino model while preserving supersymmetry.

If entanglement wedge reconstruction is exact, then (under certain assumptions) the area term in the RT formula is a c-number, indicating that the choice of a bulk quantum state does not influence the geometry. Recently Cao, Cheng, Karthikeyan, Li, and Preskill considered a generic perturbation away from exact entanglement wedge reconstruction. The optimal reconstruction was defined; based on this, an effective area function that depends nontrivially on the quantum state was defined and its properties were analyzed. Here we make one aspect of this picture more quantitative, by showing that if as expected the area term in the RT formula is of order 1/G while the bulk entropy is of order 1, then the corrections to entanglement wedge reconstruction are exponentially small (in G) relative to corrections to the area function. In the framework under discussion, there is an area function but no area operator; we discuss to what extent this is the expected behavior in holography.

We present a holographic dual of a two dimensional conformal field theory with non-hermitian but Parity-Time (PT) symmetric boundary conditions, by applying the AdS/BCFT duality and by introducing an imaginary valued scalar field localized on an end-of-the-world brane. We find that as we increase the strength of the non-hermitian PT symmetric interactions, the system experiences a spontaneous PT symmetry breaking. We also consider its Wick rotated setup as a new quantum quenched state and show that its growth of entanglement entropy can be larger than the standard results obtained from standard Cardy states.

The complex Sachdev-Ye-Kitaev (SYK) model, featuring fermions with all-to-all interactions, serves as a dual paradigm for understanding non-Fermi liquid behavior and the holographic nature of charged black holes. Two defining characteristics of the standard SYK model are its maximal chaos (Lyapunov exponent $\lambda_{\mathrm{L}}=2\pi T$ at temperature $T$), and its finite zero-temperature residual entropy. While previous studies have largely focused on couplings drawn from a zero-mean Gaussian distribution, we investigate a generalized model with a finite mean-to-standard-deviation ratio, $g\equiv J_{0}/\delta J$ of the coupling distribution in order to get deeper insight into the evolution of chaos. We find that increasing $g$ yields the following effects: (i) The system remains a fast scrambler with $\lambda_{\mathrm{L}}=A~T$, but with a suppressed coefficient $A<2\pi$. (ii) In the limit $g\to \infty$, out-of-time-ordered correlators (OTOCs) no longer exhibit exponential growth with $\lambda_{\mathrm{L}}\simeq 0$. (iii) The spectral correlations indicative of late-time chaos maintain Wigner-Dyson level spacing statistics for all values of $g$. (iv) The system preserves a finite residual entropy, albeit with reduced magnitude, for all $g$ values. We conclude that in this generalized SYK model, there is a chaotic to non-chaotic crossover. Moreover different measures of chaos decouple, demonstrating that the presence of finite residual entropy does not strictly imply maximal chaos.

We construct exact, positive-energy, normalizable wave-packet solutions of the Dirac equation in the axisymmetric potential $V=-\,v_0/\rho$ – to our knowledge, the first such solutions in any external potential. Remarkably, one family comprises only elementary functions whose longitudinal profiles reproduce the free-Schrödinger Hermite–Gauss wave packets in the nonrelativistic limit. All packets share two striking features: (i) a probability density that is pointwise decoupled from spin orientation – despite the inherent spin-orbit coupling of the Dirac equation – and (ii) a complete freezing of their time evolution at the critical coupling $v_0\to\hbar c/2$. We also present a simple scheme that maps solutions of the 2D Helmholtz equation to further exact Dirac wave packets.

We study the semiclassical limit of SL_2-opers on the four-punctured sphere in Nekrasov-Rosly-Shatashvili Darboux coordinates. Using exact Wentzel-Kramers-Brillouin (WKB) connection formulae, we express the trace coordinates of the corresponding SL_2(C) character variety as finite Laurent sums of Voros exponentials. Tropicalizing these formulae and the NRS relations gives a chamberwise integer affine linear system for the leading logarithms of the NRS coordinates. In chambers where this system is unimodular and the selected cycles form a primitive symplectic pair, the leading asymptotics agree, up to flavor-period shifts, with Seiberg-Witten periods of the N=2 SU(2) theory with N_f=4 fundamental hypermultiplets. We verify this mechanism in a sample chamber and in the weak-coupling degeneration. No global coordinate-independent recovery theorem is claimed; non-unimodular or degenerate chambers are treated as limitations of the chosen NRS chart. In the weak-coupling degeneration, we show that the NRS chart can be chosen compatibly with the plumbing limit so that the resulting chamber is unimodular and non-degenerate away from tropical walls.

We study two-point functions of half-BPS operators in maximally supersymmetric Yang-Mills theory continued to $d=4-2\,\varepsilon$ dimensions. Using supersymmetric localization on $S^d$, we derive perturbative matrix-model expressions for the $\varepsilon$-expansion of these correlators and obtain all-loop results at leading order in $\varepsilon$ in the planar limit, with extensions to finite-$N$ corrections and higher-charge operators. We compare the localization results with direct perturbative computations in flat space. At order $\varepsilon$ the two descriptions agree perfectly, while at higher orders our construction fails to reproduce the perturbative data due to the breaking of conformal symmetry away from four dimensions. Nevertheless, in the case of the dimension-two operator we conjecture an all-loop formula at order $\varepsilon^2$ by exploiting the precise form of the mismatch.

We investigate the thermodynamic stability and the stochastic thermal fluctuations of the warped black hole solution in three-dimensional General Massive Gravity. We demonstrate that the black hole is thermodynamically unstable and identify the nontrivial Davies phase-transition curves from the behavior of its admissible heat capacities. Going beyond the classical stability analysis, we study thermal fluctuations within a modified finite-time nonequilibrium extension of Ruppeiner's Hessian-based fluctuation theory. For a class of isentropic and isoenergetic processes, we derive exact on-shell angular momentum trajectories in the thermodynamic state space and compute the corresponding thermodynamic lengths. These quantities characterize relaxation processes between macrostates and provide an estimate of the associated relaxation times. Furthermore, we show that the thermodynamic geodesic equations do not admit constant-angular-momentum solutions, suggesting a continuous change of the black hole's angular momentum. Our results consistently reproduce the warped AdS$_3$ black hole limit of Topological Massive Gravity.

The Ashtekar-Lewandowski (AL) volume operator of loop quantum gravity is central to the Hamiltonian constraint, but its vertex action is usually obtained from dense spectral decompositions of finite recoupling matrices, obstructing numerical analysis on large kinematical Hilbert spaces or high-valence vertices. We formulate a matrix free action of the $SU(2)$ AL vertex volume operator in standard recoupling basis, making use of the Brunnemann-Thiemann expression for the oriented AL volume density $Q_{v}$ whose matrix elements can be generated locally from recoupling theory without forming the full matrix. Based on the Balakrishnan-Stieltjes representation of $(Q_{v}^{2})^{1/4}$ we approximate the volume by shifted-resolvent quadrature (SRQ). The resulting action uses only repeated applications of $Q_{v}$ and shifted positive linear solves, making it compatible with multi-shift Krylov methods. We prove exact preservation of the volume kernel, provide operator-norm and residual error estimates, discuss sector-wise scaling bounds, and validate the method on an embedded $K_{5}$ graph at small spin cutoffs against exact dense local-block operators. Numerical simulations show rapid convergence of vertex expectation values, controlled dependence on bound parameters, and exact preservation of zero-volume modes. We further demonstrate matrix free Monte Carlo estimates at doubled-spin cutoff $2j=250000$ beyond dense materialisation, and show that SRQ can be combined with stochastic Lanczos quadrature to estimate fixed-sector volume spectral measures without dense volume matrices.

The sign problem is a notorious obstacle for classically simulating quantum theories with fermions. We propose an effective field theory method for analyzing the sign problem. At high temperatures, a $d$+1 dimensional field theory reduces to a bosonic $d$-dimensional theory; the phase of the Pfaffian in the higher dimensional theory is encoded in an operator in the lower dimensional theory. We apply this framework to the D0-brane/BFSS matrix quantum mechanics, where the phase becomes an operator in a bosonic multi-matrix integral. Our results show that the continuum theory has a sign problem that persists in the large-$N$ 't Hooft regime. However, detecting the sign problem involves going to 10-loop order in the high-temperature expansion. This delayed onset follows from the fact that the Pfaffian phase transforms as an $O(9)$ pseudoscalar. Furthermore, the relevant diagrams give a numerically small prefactor. Consequently, ignoring the sign problem leads to a relatively small fractional error in thermodynamic quantities for temperatures $T \gtrsim \lambda^{1/3}$. However, at stronger coupling in the 't Hooft regime, the sign problem may become more severe. Finally, we initiate the application of this framework to higher-dimensional maximally supersymmetric Yang-Mills theories.

Type II string theory in the absence of orientifolds requires the target spacetime $X$ to admit a spin structure. We show that this familiar requirement arises directly from the consistency of the worldsheet GSO projection. The obstruction is a mixed global anomaly between $(-1)^{\text{F}_L}$ and the target space background data, detected by the spin bordism group $\Omega_3^{\mathrm{spin}}(B\mathbb{Z}_2\times X)$. We compute the relevant mixed bordism group and identify the bordism class of the GSO-projected worldsheet theory. For smooth target spaces, vanishing of the anomaly reduces to the condition that $X$ is spin, while for general orbifolds $[\hat{X}/G]$, $\hat{X}$ has to carry a $G$-equivariant spin structure. We also classify all possible theta angles in the worldsheet theory and show that they correspond to all possible continuous and discrete background fields of the target space theory visible in string perturbation theory.

While the Standard Model (SM) of particle physics contains the most precise set of predictions ever devised by humanity, that precision comes at a cost. The strange nature of the Higgs particle requires its parameters to be tuned so precisely that if the SM is indeed the true description of reality, one is forced to wonder how such a miracle as galactic structure and life could occur. Instead, we search in this work for a natural explanation. The concept of naturalness is comprehensively explored, and a new tuning measure proposed, with an aim to place it on well-defined Bayesian footing. We then turn this measure on to the analysis of a class of intriguing new physics - Composite Higgs models. These effective models are the result of a plethora of underlying theories, and they allow the production of a naturally light Higgs particle, appearing as the SM Higgs at low energy. We establish the background required to appreciate the N-site 4D Composite Higgs model, and subsequently focus on the simplest incarnations of this class. A global fit is performed on the Minimal 4D Composite Higgs model (M4DCHM), with strong exclusion bounds placed on collider search channels. We analyse any improvement in tuning that could be gained from several extensions to this model. The Leptonic M4DCHM is explored, with a composite tau lepton embedded in various representations. The possibility of a dark matter candidate existing in the Next-to-Minimal 4DCHM is considered. Ultimately, we are able to define what, if any, benefit to naturalness can come to the Composite Higgs sector by introducing these extensions.

We study two-point functions of protected scalar operators in N=4 SYM, focusing on their transcendentality properties in dimensional reduction, where quantum corrections are subleading in the regulator. We compute the correlators explicitly through two loops and operators up to classical dimension 10, for all trace structures. The one-loop correction is universal. At two loops, we find a controlled partial breaking of uniform transcendentality for higher-dimensional operators, which can be cancelled by suitable combinations of correlators in a fully predictable way. A main result is a complete planar extrapolation for two-loop correlators at arbitrary dimension and trace structure, whose dependence is entirely controlled by the number of stress-tensor multiplet factors in the operator. The perturbative results agree with localization predictions in all cases where comparisons are possible.

We consider the problem of identifying a topological order based on bulk entanglement of the ground-state wavefunction. Previous work showed that some universal information can be extracted from multi-entropy measures, a class of multipartite entanglement measures obtained by applying permutation operators exchanging the degrees of freedom between different replicas of the wavefunction. It remains an open question to what extent such entanglement measures can be used to extract any universal information from the ground state. Here we show that the topological partition function $Z(M)$ of a manifold $M$ can be extracted provided that $M$ satisfies a topological condition which we term “local achirality". We show that locally-achiral manifolds can be used to extract universal properties of 2+1d topological phases that go beyond the $S$ and $T$ matrices. As a first step towards classifying locally-achiral manifolds, we show that, in four dimensions, such manifolds have vanishing Pontryagin number. We relate this property to the existence of beyond-cohomology time-reversal symmetry protected topological order (T-SPT) in four dimensions. Finally, we present an entanglement measure that detects this nontrivial T-SPT.

We study signatures of quantum chaos in finite-loop truncations of the planar dilatation operator in the $\mathfrak{su}(2)$ sector of $\mathcal N=4$ super Yang-Mills and its $\beta$-deformation. These truncations define holographically motivated long-range deformations of the nearest-neighbour XXX spin chain. At one-loop the model is integrable, while the all-loop planar theory is expected to again be integrable. Finite-loop truncations therefore provide a natural setting for investigating how chaotic behaviour emerges between these two integrable limits. We analyse this question using spectral statistics, eigenvector diagnostics and spread complexity. We find that the two- and four-loop truncations develop GOE-like level statistics at sufficiently large coupling but with features characteristic of weak integrability breaking. The integrability breaking at four-loops is weaker than at two-loops and the critical coupling at which chaos occurs is larger, at least for long spin chains. The three-loop truncation does not show the same onset of chaos in the range studied. Eigenvector diagnostics show that the corresponding eigenstates remain less random than GOE vectors, indicating weak ergodicity and multifractality. Finally, we can identify signatures of the eigenvalue and eigenvector chaos in the Krylov-space data. Namely, we demonstrate a correlation of the level spacing statistics with the peak of spread complexity and disorder on the Krylov chain. The delocalisation of the initial state in the Hamiltonian eigenbasis is shown to strongly affect the saturation of complexity. Our results suggest that finite-loop dilatation operators are not generic long-range spin chain Hamiltonians, but already display patterns consistent with the restoration of integrability in the all-loop planar theory.

We elucidate the physical nature of ghosts above the multi-particle threshold by contrasting them with unstable particles in quantum field theory. We first consider the asymptotic formulation, where ordinary positive-norm one-particle states can be unstable and decay, whereas ghosts survive asymptotically without decaying, yet admit no particle interpretation due to interference with the multi-particle component which masks the negative-norm one-particle state. This distinction originates from two different analytic structures of the dressed propagator, whose complex conjugate poles lie in the first or second Riemann sheet in the ghost or ordinary case, respectively. Ghost resonances are, in principle, phenomenologically distinguishable from ordinary ones, being narrower and exhibiting weaker interference between positive- and negative-energy peaks. We then formulate the quantum field theory in a finite interval of time and, working within a suitable approximation for the dressed propagator, find that finite-time effects amplify differences in the resonant behavior and give rise to new features, such as higher peaks in ghost resonances. Distinct temporal regimes are also identified: for times shorter than the inverse width, an approximate free-particle description is valid, whereas at later times interactions and interference effects dominate, leading to decay or multi-particle masking. Complex poles in the dressed propagator emerge only at late times and become complex-conjugate pairs asymptotically, determining the asymptotic dynamics. This study supports the absence of freely propagating ghost particles in the asymptotic limit.

A long-standing conjecture claims that Virasoro conformal blocks exponentiate in the semiclassical limit $c \to \infty$ with $h/c$ finite. However, this has been proven only for four-point blocks on the sphere and one-point blocks on the torus. Here we extend the proof to general conformal blocks for higher-point functions and higher-genus backgrounds in arbitrary channels. The statement is to be understood at the level of a formal power series. Our proof builds upon a novel extension of the oscillator method for the computation of conformal blocks to cases where three internal lines meet at a vertex. This extension also gives a new constructive method to compute global conformal blocks in 2d CFTs at general genus.

This paper generalizes the method used in the previous article 2512.09930 to black $p$-brane thermodynamics in arbitrary dimensions containing black holes and strings as special cases: thermodynamic quantities can be derived without constructing the corresponding black $p$-brane solutions. We further extend the discussion to black holes or $p$-branes involving a general scalar coset.

We derive an analytical approximation for the graviton propagator from Asymptotic Safety. We find neither extra poles nor indications of unitarity or causality violations in the spin-two sector. Our results strengthen the case that Asymptotic Safety does not introduce new degrees of freedom, and thus propagates the same field content as General Relativity. We also identify the underlying mechanism: the residues of spurious poles in finite-order derivative expansions approach zero as the order is increased.