Papers

Ultra-hot Jupiters (UHJs) are the most extreme class of close-in gas giants. We combine 13 transit observations of KELT-9b from six high-resolution spectrographs spanning 316-960 nm to conduct a census of atomic and ionised species in its atmosphere using cross-correlation. We introduce an agnostic strategy to mitigate velocity-axis correlation bias in parameter uncertainties. A forward model of the cross-correlation function (CCF) and a Monte Carlo bootstrap of flux uncertainties measure the CCF-derived velocity profiles of individual species and assess their time dependence. The stellar reference frame is determined, yielding Vsys = -17.99 +/- 0.04 km s^-1 and v sin i = 110.54 +/- 0.05 km s^-1. All previously reported species are recovered except Tb II; the combined analysis detects 29 species, including C I, Si I, K I, V II, Co I, Zr II, La II, Pr II, and Nd II. All 29 recovered species exceed the conventional 5 sigma threshold in the bootstrap Mahalanobis significance test. The collective Doppler shifts show a blueshift of Delta Voffset = -8.02 +/- 0.17 km s^-1 relative to the stellar rest frame, indicating a net day-to-night wind, with no explicit time dependence when all species are homogenised. These results provide a catalogue of detections and retrieved Doppler parameters for a companion study of KELT-9b.

Asymmetric stellar explosions impart an impulse, or natal kick, to their compact remnants by linear momentum conservation. Black hole (BH) natal kicks are often assumed to be weaker than neutron star kicks because of greater mass accretion, and they are harder to constrain observationally because isolated BHs are electromagnetically faint. Using 3D general-relativistic magnetohydrodynamic simulations lasting up to ~30 s, we show that BHs formed from the collapse of rapidly rotating massive stars (collapsars) threaded by strong large-scale magnetic fields can acquire large natal kicks 10^2–10^3 km/s. Asymmetric electromagnetic outflows, magnetic-flux eruptions, and the gravitational pull of aspherical ejecta shape the kick magnitude, and their relative contributions depend on the progenitor structure, magnetic-flux history, and BH spin. More rapidly spinning BHs receive stronger, more nearly spin-aligned kicks, primarily through the gravitational pull of asymmetric jet-driven ejecta. Delayed transitions to the magnetically arrested state produce more asymmetric outflows and can generate even stronger recoils. Because the large-scale magnetic flux required in our models is also a key ingredient of relativistic gamma-ray burst jets and their associated energetic, jet-driven supernovae, natal kicks may be a natural consequence of magnetized collapsars. Such kicks could substantially alter BH retention in dense stellar environments, binary survival, spin-orbit misalignment, and the viability of magnetized collapsars in producing highly spinning BHs within the pair-instability mass gap.

We establish an obstruction to dimension-independent local color-kinematics duality for the one-loop six-gluon amplitude in pure Yang-Mills theory. The maximal and fivefold cuts admit a polynomial solution, but no allowed contact correction reproduces a required fourfold cut. We consider parity-even Lorentz-polynomial numerators, multilinear in arbitrary transverse external polarizations, with conventional quadratic propagators, global graph identities, and pointwise formal-color cut equality. No dimension-specific Gram identities are imposed. In the sector containing one polarization inner product, the higher cuts determine the complete remaining contact freedom. Every such correction obeys an exact coefficient relation that the required physical residual violates, independently of the chosen higher-cut solution. The violation is proportional to the transverse internal vector-state count and is nonzero at every nonzero value of its algebraic continuation. Momentum homogeneity extends the degree-six contradiction to all finite polynomial degrees. An all-multiplicity polynomial family retaining these representation conditions is therefore excluded.

We show that a formal deformation quantisation of degree-2 differential graded symplectic (QP) manifolds gives rise to the complete classical kinematics and dynamics of NS-NS supergravity, or, equivalently, its duality-covariant formulation as double field theory (DFT). Our construction extends the standard tensor hierarchy, encoding gauge parameters and their redundancies, to a complex encoding the physical fields, field strengths, and Bianchi identities. We present two realisations of such a complex: 1) a differential graded Lie algebra built from the algebra of functions on the QP-manifold; 2) a cochain complex that reproduces the tensor hierarchy representations conjectured in the literature and can be interpreted as a classical BV complex of the linearised theory. The deformation that gives rise to the dilaton comes from a normal-ordered graded Moyal–Weyl star product. This naturally extends the duality structure group from $\mathrm{O}(D,D)$ to $\mathrm{O}(D,D)\times\mathbb{R}^+$. Finally, our framework provides a direct route to constructing the DFT action via local double Lorentz invariance, and offers a transparent algebraic foundation for curvature tensors in generalised Cartan geometry.

We propose a microscopic scheme that allows to describe ac Josephson effect in superconducting junctions in terms of dissipative quantum dynamics of subgap Andreev bound states. This approach is particularly useful for highly transparent junctions at subgap bias voltages in which case a non-trivial combination of Landau-Zener tunneling between Andreev levels and their instability due to quasiparticle escape into continuum may play an important role. In the low bias regime, we evaluate the non-equilibrium current-phase relation of superconducting junctions at arbitrary transmissions and identify a sub-Ohmic phase-dependent dissipative contribution to the current controlled by quasiparticle dynamics near the superconducting gap edge.

The Sommerfeld enhancement (SE) from long-range interactions, and the related bound state formation rate, can be important non-perturbative inputs to dark matter (DM) annihilation signals. In the presence of singular interaction potentials, the conventional boundary conditions of the Schrödinger equation fail, obfuscating the computation of SE, its physical origin, and its relation to bound states in such potentials. In this work, we clarify the origin of SE in singular potentials in a two-fold manner: using the framework of velocity power counting in non-relativistic effective field theory (NREFT), and via position-space regularization of singular potentials at short distances to compute SE and bound states. We illustrate our findings through the case of pseudoscalar mediators interacting with massive Dirac DM. We find that when such a system arises from a UV-complete theory at weak coupling, no SE is generated. However, in the case of a derivatively coupled pseudoscalar, where the interaction is described by a higher-dimension operator in an effective theory, SE (and bound states) can occur for weak couplings if there is a large hierarchy between the cutoff scale of the effective theory and the dark matter mass. This SE is sensitive to the choice of the UV completion of the potential at short distances, but a non-negligible SE can persist even when the UV physics alone would not generate any SE; we elucidate the interplay of UV and IR physics in this case.

We study the electromagnetic radiation and wave-optical imaging of an oscillating electric dipole orbiting a Kerr black hole. We derive the effective 4-current associated with a pointlike oscillating electric dipole in curved spacetime, and use black hole perturbation theory to compute the resulting radiation field at future null infinity, from first principles. We then develop a wave-optical imaging framework for a moving electromagnetic source in curved spacetime. We obtain images of an orbiting electric dipole, displaying relativistic beaming, gravitational lensing, and Einstein rings. We study polarization-dependent scattering by comparing the images produced by spinning dipoles with opposite helicities, finding a displacement that roughly decreases with the inverse of the radiation frequency, as expected for a beyond-geometric-optics effect. Our results provide a first-principles benchmark for beyond-geometric-optics descriptions of electromagnetic radiation in Kerr spacetime.

Cosmic rays can significantly alter the abundances of certain species in the upper layers of planetary atmospheres, especially in terms of their biosignatures. To fully understand the extent of this effect, it is essential to accurately model the interactions of cosmic rays with planetary magnetic fields. We used the CosmicTransmutation code to study the effect of a planetary magnetic field on the flux of galactic cosmic rays and stellar energetic particles. We found that both particle sources are significantly affected by magnetic fields, even though the effects are different due to the varying energy ranges that characterize each source. The stellar energetic particle energy flux is significantly higher for an Earth-like planet with no magnetic field, but with a magnetic field of 30\,$\mu$T or higher, the energy flux of the two sources becomes comparable. We find that the atmospheric pressures the cosmic rays reach are significantly lower than those found in previous studies using simpler models and that the effect of the atmospheric composition on this outcome is small. We found similar results using the model for the exoplanet K2-18b, even though the radius is significantly larger than Earth's. Because the two cosmic ray sources cover different energy ranges and are affected by the magnetic field in different ways as a result, both sources should be considered. Future studies should focus on combining an accurate modelling procedure of the interaction between cosmic rays and planetary magnetic fields with atmospheric chemistry models. This will require general circulation models to capture the latitudinal and longitudinal dependencies in full.

Thermoelectric effects in superconductors are generally suppressed by the approximate particle-hole symmetry of the quasiparticle spectrum, but can become pronounced near defects that host particle-hole asymmetric bound states. Here, we investigate how the local thermoelectric response is modified when multiple vortices or nonmagnetic impurities are brought into proximity. Using a lattice Bogolioubov-de Gennes approach combined with linear-response tunneling theory, we calculate the spatially resolved density of states and Seebeck coefficient in $s$- and $d$-wave superconductors. We find that the thermoelectric response can be strongly enhanced when spatially extended defect-induced states overlap. For vortices, this enhancement persists beyond the immediate core regions and originates from interference between vortex-bound states. For impurities, the thermoelectric response exhibits a comparable dependence on impurity separation in the $s$- and $d$-wave cases, although the associated spectral reconstruction is considerably more localized in the $s$-wave superconductor. Our results show that the spatial extent and interference of defect-induced quasiparticle states provide a means of controlling local thermoelectricity in inhomogeneous superconductors, with potential relevance for cryogenic thermoelectric sensing and energy conversion.

This work presents commissioning of the Hat Creek Outrigger (HCO), a dual-polarization 256-element radio interferometer that is part of the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst (CHIME/FRB) Outrigger project. Driven by a complex radio-frequency interference (RFI) environment that consistently contaminates $\sim40\%$ of HCO's usable bandwidth, we implement an improved calibration scheme using Gaussian Process Regression to recover complex gain solutions over RFI contaminated channels. To validate our method, we test the performance of the array using known transients and continuum sources over relevant timescales used for fast-transient research ($\lesssim$ seconds). We find that our updated calibration scheme results in a $\sim1.69\times~\mathrm{to}~1.87\times$ improvement in the array's point-source sensitivity, while simultaneously maintaining noise properties consistent with thermal statistics. As a result, we find that the array performs consistently within theoretical expectations across $\gtrsim80\%$ of HCO's bandpass. We further observe an improvement in the interferometric performance after applying recently developed spatial filtering techniques for RFI mitigation, which rely on accurate calibration solutions for effective removal of unwanted interference. We conclude that our approach provides a valid framework for improving calibration solutions over RFI contaminated channels for large-$N$ interferometric arrays more broadly. Our work motivates future development of more sophisticated techniques to recover astrophysical information in RFI-contaminated channels, departing from the historical practice of discarding them outright.

We study the representation theory of the vertex operator algebras (VOAs) associated with rank-one 4d $\mathcal{N} = 3$ superconformal field theories. For the $\mathbb{Z}_3$ S-fold theory, whose VOA $\mathcal{W}_{\mathbb{Z}_3}$ has central charge $c_{\mathrm{2d}} = -15$, we use the $\mathcal{N} = 1$ Lagrangian description to obtain the unflavored Schur index in terms of Dedekind eta functions, while Wilson-loop indices yield the unflavored non-vacuum characters. These characters all solve a modular linear differential equation (MLDE) whose solution space also contains a logarithmic character. Combining flavored MLDEs from null states with Zhu's associative algebra and a free-field realization, we study four highest-weight modules of $\mathcal{W}_{\mathbb{Z}_3}$ and their flavored characters in closed form. A parallel analysis applies to the $\mathcal{N} = 3$ theories obtained by gauging a discrete $\mathbb{Z}_n$ flavor subgroup of $\mathcal{N} = 4$ $U(1)$ and $SU(2)$ super-Yang–Mills, for which we also obtain closed-form Schur indices and a new free-field realization of the VOA of the $\mathbb{Z}_4$ quotient.

We associate a gravitational Hilbert space $\mathbf{H}_\sigma$ to any closed compact $(d-1)$-manifold $\sigma$ with real metric. A quantum state $\mathcal{J}(\sigma)$ is a $d$-manifold bounded by $\sigma$ and equipped with elliptic data. An inner product is defined by gluing states pairwise across $\sigma$ and evaluated by viewing the resulting closed $d$-manifold as a boundary condition on the gravitational path integral (GPI) over $(d+1)$-manifolds. If $\sigma$ is nonempty and the GPI is dominated by a single $(d+1)$-manifold $M$ in the $G_N\to 0$ limit, then $M$ contains a Lorentzian CRT fixed-point set, providing $\mathcal{J}(\sigma)$ with a classical spacetime interpretation. Conversely, given a finite Lorentzian domain with edge $\sigma$, a state $\mathcal{J}(\sigma)$ may be associated to it by deforming its initial data off the real Lorentzian section and retaining only elliptic data. This establishes a broad correspondence between non-asymptotic spacetime regions and quantum states. Assuming that $\mathbf{H}_\sigma$ factorizes over connected components of $\sigma$, our framework admits operators and partial traces. This allows us to explore the information-theoretic structure of the states we define. As an example, we construct a family of states by deforming partial Cauchy slices $\Sigma$ that straddle a two-sided black hole; $\sigma$ consists of two spheres. We construct the reduced state on one sphere and find that its Rényi entropies are positive, monotonic, and sensitive to all aspects of $\Sigma$ and its complex deformation. The von Neumann entropy, however, is controlled only by the maximin surface in the causal domain of $\Sigma$, independently of other parameters, so long as the complex deformation does not vanish. Our proposal may thus explain the efficacy of tensor network toy models of holography while transcending their limitations.

We construct Floquet codes from the Majorana XYZ subsystem code with local realizations both in qubits and directly in microscopic Majorana modes with lattice size $L\times L$. For a three-step cycle, odd $L$ supports one static logical qubit, whereas even $L$ supports two. For even $L$ with $L=4n+2$, both logical qubits admit time-independent Pauli representatives. For $L=4n$, by contrast, one Pauli of the second logical qubit remains fixed throughout the cycle, while every representative of its conjugate must evolve through the measurement cycle. This distinction follows from a parity-dependent algebraic obstruction and disappears when the cycle is reduced to two steps, which restores a fully static logical pair. Thus, the same encoded logical degree of freedom can be switched between static and partially dynamical forms by the measurement schedule. With one open direction, suitable protocols can be implemented using only local Majorana parity measurements. To our knowledge, this is the first Floquet-code family with both a local qubit representation and a direct microscopic Majorana realization.

We construct and study lattice realizations of finite 2-group symmetries in ${2+1}$d quantum lattice systems with finite-dimensional tensor-product Hilbert spaces. We focus on two broad classes of 2-groups with 0-form symmetry group $G$ and 1-form symmetry group $A$: split 2-groups with trivial Postnikov class ${[\beta]\in\mathcal{H}^3(G,A)}$, and central 2-groups with trivial action ${\rho\colon G\to\text{Aut}(A)}$. In both cases, we construct symmetry operators on the full tensor-product Hilbert space that become 2-group symmetry operators when restricted to the topological subspace of the lattice $A$ 1-form symmetry. While the lattice split 2-group symmetry operators are onsite, the lattice central 2-group symmetry operators are not, and can only be made onsite after introducing ancillae. We extensively explore various manifestations of $\rho$ and $[\beta]$ for these lattice 2-group symmetry operators and demonstrate their agreement with expectations from quantum field theory. These manifestations arise in the transformation of operators carrying symmetry charge, the structure of lattice 2-group symmetry defects, and the dual fusion 2-category symmetries obtained by gauging the lattice 2-group symmetries. We further propose families of local symmetric Hamiltonians for both classes of lattice 2-group symmetries and identify exactly solvable limits lying in phases with spontaneous 2-group symmetry breaking and nontrivial symmetry-enriched topological order. In one such limit, the gauged Hamiltonians are exactly solvable lattice realizations of the corresponding 2-group gauge theories, whose ground-state degeneracies we calculate.

We study the phenomenology of the topological portal between QCD and pseudo-Nambu–Goldstone dark matter proposed in~\cite{Davighi:2024zip}. We construct a weakly coupled ultraviolet completion in which a vector mediator gauges a baryonic current of QCD for the light quark flavours, linking it to dark pions parametrizing the coset $SU(2)_D/U(1)_D$. As a concrete anomaly-free realisation, we gauge the leptophobic combination $B_1+B_2-2B_3$. This completion allows us to calculate thermal coannihilation beyond the regime of validity of the pion effective theory, including resonant mediator effects, and to test the resulting thermal target with boosted-dijet and monojet searches at the LHC, with LEP and bottomonium data, and with flavour observables. These collider signatures follow directly from the symmetry structure required by matching onto the topological portal. At low momentum transfer, the topological operator instead controls the decay of the heavier dark pion, which is long-lived. A future high-statistics $Z$ factory, such as FCC-ee, can probe the low-mass region where hadron-collider sensitivity deteriorates, providing complementary tests of the topological portal.

Mock catalogues for the Roman Galaxy Redshift Survey benefit from galaxy populations whose physical properties and emission-line observables are predicted within a single galaxy-formation framework. We present an emulator-assisted calibration of the semi-analytic model Galacticus for this purpose. Because direct posterior exploration with Galacticus is computationally prohibitive, we train Gaussian-process emulators to predict the observables entering the calibration likelihood. We calibrate the model simultaneously to stellar mass functions, the star-formation-rate function, size–mass relations, the black-hole mass–velocity-dispersion relation, the mass–metallicity relation, and H-alpha luminosity functions. Combining these observables breaks parameter degeneracies present in individual fits and can reveal tensions where different observables favour different regions of parameter space. The resulting model provides a broadly successful simultaneous description of the calibration data, including the abundance and redshift evolution of H-alpha emitters relevant to the Roman Galaxy Redshift Survey. A direct Galacticus calculation at the maximum a posteriori point gives similarly good agreement with the calibration data, providing a calibrated galaxy-formation model suitable for constructing physically consistent mock catalogues for Roman.

We develop an Efficient Wake procedure for computing the distribution of hadrons originating from jet wakes in heavy ion collisions - the hydrodynamic response of a droplet of quark-gluon plasma to the energy and momentum deposited in it by high-energy partons propagating through it. The procedure employs the linearity of linearized hydrodynamics and takes account of the effects of both longitudinal expansion and transverse radial flow on the hydrodynamic evolution of the wakes and on the resulting particle production at the freezeout hypersurface. It makes repeated use of template solutions to linearized hydrodynamics in a Bjorken flow background with no transverse flow, templates that need only be computed once, and uses suitable rotations and boosts to map fluctuations from these templates to fluctuations at a point on the freezeout hypersurface in a way that incorporates the effects of the radial flow. We benchmark this procedure by comparing its results to results obtained from full $(3+1)$-dimensional nonlinear hydrodynamics calculations, find reasonable agreement, and find that our Efficient Wake procedure yields a much better description of the distribution of hadrons originating from jet wakes than does the older oversimplified procedure employed in the Hybrid Model. And, the Efficient Wake procedure is computationally efficient: it is at least tens of thousands of times faster than full nonlinear hydrodynamics calculations. Hence, we anticipate that when our new procedure is implemented in Monte Carlo analyses of jets in heavy ion collisions, for example in the Hybrid Model, it will greatly improve the description of the soft component of many jet and jet substructure observables as compared to experimental data.

The concept of quantum geometry has recently led to reinvigorated insights in a wide range of fields including physical responses, superconductivity, and optical transitions, with effects most pronounced in systems with nearly flat dispersion. Here, we show that it plays an essential role in charge density wave formation (CDW) – an important physical phenomenon that is responsible for driving various sharp changes in material transport properties including metal-insulator transitions. We derive an effective Ginzburg-Landau theory including uncharted contributions and, as a highlight, discover a general criterion for both CDW formation and commensurability transitions where underlying electron-phonon interactions manifest purely as electronic quantum geometric enhancements/suppressions. We benchmark our framework in a class of transition-metal dichalcogenides and resolve a longstanding puzzle where well-established purely kinetic CDW criteria fail in describing the correct ordering wavevector. Besides rendering robust criteria and fundamental insights that are immediately relevant to several experimental charge ordering systems, our theory can also be applied directly to other phonon-mediated phases such as superconductivity, and can be used as an important tool to explore the interplay between various such states. More generally, our framework provides a recipe for investigating the role of quantum geometry in phase transitions.

As parameters are varied, the set of saddle points contributing to a (path) integral may change discontinuously, leading to a corresponding change in the asymptotic expansion of the integral. This behavior is known as the Stokes phenomenon. We explore this phenomenon in the context of the analytic continuation problem relating the AdS/CFT and dS/CFT correspondences. In this paper, we study these correspondences for three-dimensional pure gravity and two-dimensional Liouville theory, using independent calculations in bulk minisuperspace and in the boundary Liouville zero-mode. In the bulk, the dS contour selects a single saddle and yields the tunneling wave function. Upon continuation to AdS, the contour instead selects an infinite family of saddles. The boundary calculation independently reproduces the same Stokes structure at leading semiclassical order, providing a nontrivial holographic consistency check. Our construction also offers a contour prescription for the conformal factor problem in Euclidean AdS$_3$ quantum gravity within minisuperspace.

We introduce $Krylov$ $edge$ $spectroscopy$, a many-body operator-space protocol for detecting and classifying one-dimensional symmetry-protected topological phases from local boundary dynamics. A Hermitian boundary operator generates a semi-infinite Krylov hopping chain whose boundary weight obeys an exact zero-frequency normalizability criterion. Open-periodic, boundary-bulk, and symmetry-preserving boundary-perturbation tests identify protected boundary memory, while a finite-depth leakage residual certifies when an explicitly reconstructed operator is already near zero frequency. Classification minimizes the normalized commutator over symmetry-resolved boundary operators in fixed charge sectors. For bosonic $\mathbb{Z}_N \times \mathbb{Z}_N$ phases, the recovered endpoint charge gives the cohomology label. The method requires no explicit ground-state wavefunction, entanglement spectrum, or guessed dressed edge or string operator. For finite-range Hamiltonians, locality eliminates full-system exact diagonalization and guarantees thermodynamic-limit convergence at fixed Krylov depth. We demonstrate the protocol in cluster, clock, Haldane, and trivial spin-1 chains. In the exactly solvable cluster chain, the transition between the gapped topological and trivial phases manifests as a localization-delocalization transition of the Krylov edge mode on the Krylov chain. This transition occurs precisely at the bulk gap closing, and its localization-length exponent coincides with the Ising correlation-length exponent. Through operator Krylov dynamics, our work turns local boundary evolution into a direct spectroscopy of many-body topology.

We present the data reduction and transient detection pipeline for the Wide-field Infrared Transient Explorer (WINTER) surveyor and report its on-sky performance. The WINTER camera utilizes cost-effective InGaAs sensors as alternatives to traditional IR sensors, and is mounted on a dedicated 1-m robotic telescope at Palomar Observatory. The WINTER camera has six detectors producing a combined field-of-view of 1.2 sq. deg. equipped with y, J, and shortened-H bands. WINTER saw first light in June 2023 and has been operating robotically since. The WINTER data processing pipeline ($\texttt{winterdrp}$) has been implemented within the broader framework $\texttt{mirar}$: a modular, open-source $\texttt{python}$ package developed for realtime processing of images from time-domain surveys. $\texttt{winterdrp}$ performs end-to-end data processing implementing data reduction and image subtraction to go from raw dithered WINTER images to transient alerts in the $\texttt{avro}$ format, which are then sent to $\texttt{SkyPortal}$ for vetting and follow-up. During a year of observations in 2024, WINTER achieved J-band median 5-$\sigma$ depths ranging from $18.1-18.8$ mag (AB) on its six detectors in 960 second integrations as part of its survey, with an astrometric accuracy of $\approx0.2$ arcsec (a fifth of a pixel) and a detector-performance limited photometric accuracy ranging from $\approx0.09-0.18$ mag for its six detectors. We present early science results from WINTER, which include the identification of a stellar merger in M31, dust-enshrouded outbursting young stellar objects and classical novae in the Galactic plane, NIR followup of known supernovae, and multi-messenger follow-up of neutrinos, gravitational waves, fast X-ray transients and gamma-ray bursts.

An eV-scale QCD axion comprising the observed dark matter (DM) abundance is expected to generate a photon line at infrared energies that would be observable or near-observable in data collected by the James Webb Space Telescope (JWST), as might more general axion-like particles (ALPs) over a broader range of masses and couplings. This has motivated a number of efforts to either forecast JWST sensitivities to a QCD axion or realize them through analyses of publicly available datasets. At present, no consensus exists; leading analyses disagree by as much as an order of magnitude in terms of axion-coupling sensitivity, implying an orders-of-magnitude discrepancy in raw flux density sensitivity, and consistency between the analyses and prior forecasts is unclear. We address these outstanding discrepancies with a bespoke data reduction and flexible nonparametric inference procedure that lead to well-controlled and robust limits on the decay of eV-scale axion DM, consistent with previously forecasted sensitivities. We further demonstrate that the strongest previously claimed sensitivities exceed those attainable by any analysis of the datasets from which they were derived. We exclude QCD axion DM for masses between $500\,\mathrm{meV}$ and $2.5\,\mathrm{eV}$ using NIRSpec data, while setting limits on ALP DM complementary to other astrophysical constraints at masses between $100\,\mathrm{meV}$ and $500\,\mathrm{meV}$. However, we find the sensitivities to be systematically limited, and therefore unlikely to be improved upon by ongoing data collection or re-analysis unless instrumental modeling and data reduction pipelines improve considerably.

In a broad class of inflationary models, the leading non-Gaussian bispectrum arises from one-loop rather than tree-level processes. However, realistic one-loop contributions remain largely unexplored for phenomenologically relevant masses beyond the simplest bubble diagrams. This is insufficient for observational purposes because, for generic masses and couplings, triangle contributions need not be suppressed relative to their bubble counterparts. We present a numerical method for evaluating scalar one-loop in-in diagrams contributing to the inflationary bispectrum at fully general external momenta and masses. First, the Witten-Feynman parameterization reduces de Sitter loop integrals to generalized Euler-Mellin integrals governed by Symanzik graph polynomials. Unfortunately, the exponents of the polynomials in the integrand become complex for sufficiently heavy masses, leading to sign problems when evaluating the integral using standard numerical techniques. To remedy this problem, we introduce the reduced Schwinger method, which evaluates a highly oscillatory subintegral analytically to obtain a Gauss hypergeometric kernel, leaving the remaining integrals to standard numerical quadrature. We validate the algorithm by reproducing known analytic results for the tree-level bispectrum in terms of ${}_3F_2$ functions, and then apply it to the one-loop bubble and triangle contributions at general kinematics. This yields the first direct numerical evaluation of the complete one-loop bispectrum, valid across the full kinematic range. Using these numerical results, we construct bispectrum templates and compare them with CMB data.

Within the Hot Big Bang picture, as usually interpreted in the standard LCDM framework and its hierarchical scenario for elliptical galaxy assembly, the cosmic microwave background (CMB) is associated with photon decoupling near z = 1100, and the Planck angular power spectra are taken to constrain primordial 10^{-5}-level fluctuations that seeded later structure formation. We show that, in addition to the LCDM model, two cosmological models based on Milgromian gravitation are also consistent with the Planck CMB power spectra, such that the spectra do not uniquely favor dark matter. We then argue that the published CMB spectra are not direct observables, but reconstructed quantities obtained after foreground subtraction and corrections for weak lensing, photon-electron scattering, and gravitational energy redshifts along the line of sight. Because these steps depend on the assumed growth history of structure, the resulting spectra depend on the cosmological model adopted in the data reduction. We also note that the reported correlation between the hemispherical CMB asymmetry and the excess on-sky distribution of elliptical galaxies may point to matter inhomogeneities spanning large portions of the observable universe. Motivated by observational tensions in the hierarchical assembly of elliptical galaxies, we review the consequences of elliptical-galaxy formation on a downsizing timescale, which implies a significant extragalactic foreground released at 15< z < 20. Under conservative assumptions, this foreground contributes at least 1.4 percent of the observed CMB energy density and has not been included in CMB analyses. It follows that the reported 10^{-5}-level fluctuations is introduced by overly aggressive foreground cleaning.

Ultralight dark matter could be an axial vector field, a simple possibility which motivates new experiments. Couplings to fermions, through an axial vector current or a dark electric dipole moment operator, lead to enhanced spin torques compared to axion dark matter. An axial vector also has an analogue of the axion-photon coupling, though its consistent realization requires a photon mass. Under this "axial-photon" coupling, the effect of a background magnetic field is suppressed, and the strongest experimental probes involve polarimetry, electric fields in superconducting cavities, and the cosmic microwave background. Since a massive axial vector is dual to a massive two-form, our results also apply to "Kalb-Ramond" dark matter motivated by string theory.

We study the magnetic deformation of the $(2+1)$d Ising CFT, also known as Ising Field Theory (IFT), on a spatial sphere. We use the Fuzzy Sphere (FS) Ising setup and compare it to the Truncated Conformal Space Approach (TCSA). Universal, infinite volume, IFT quantities are extracted by modeling the effects of curvature, as expected from a local effective theory on the sphere. We find evidence that IFT contains a bound-state with a small binding energy. Locality also allows us to extract the same IFT quantities from different angular momentum sectors, providing additional consistency checks.

We identify the precise F-theory backgrounds dual to the heterotic string theory compactified on $T^2$ with general Wilson lines turned on for one of the two $E_8$'s. The elliptic curve describing the backgrounds is expressed in terms of orthogonal modular forms and is determined so that it admits a nontrivial scaling limit yielding a rational elliptic surface. We confirm that the curve reduces to the Seiberg-Witten curve for the E-string theory in this limit. We conjecture that our result gives an explicit inverse period map for algebraic K3 surfaces polarized by the $U\oplus E_8(-1)$ lattice.

Recent developments in the Swampland program motivate the existence of a mesoscopic Dark Dimension of size $0.1-10~\mu$m with bulk Majorana fermions in the $1-10$ keV range. Motivated by this, we derive terrestrial constraints on Majorana bulk neutrinos as a function of the compactification radius. After determining the mass spectrum and mixing structure, we confront the model with Daya Bay neutrino-oscillation data, the KATRIN beta-decay bound, and the KamLAND-Zen neutrinoless double beta decay limit. For the latter, we obtain a closed-form expression for the effective Majorana mass and show that the Kaluza-Klein tower efficiently screens neutrinoless double beta decay when the nuclear momentum exceeds the compactification and Majorana scales. In the Dark-Dimension window, beta decay provides the strongest constraint, pushing the allowed Yukawa couplings down to $\mathcal{O}(10^{-3})$. Oscillation and neutrinoless double beta decay searches remain complementary, probing smaller and larger Majorana masses, respectively. For completeness, outside the Dark-Dimension assumptions, we investigate how hierarchical Majorana masses can weaken the constraints.

We derive an entropy formula for recently proposed tensor network models which prepare diffeomorphism invariant states of topological quantum field theories with non-compact and/or continuous gauge groups. We show that our entropy formula generalizes the notion of “topological entanglement entropy” to incorporate the infinite number of particle-like excitations in such theories. When our networks are endowed with gauge group $\mathrm{SL}(2,\mathbb{R})$, we can interpret them as models of three-dimensional gravity with small Newton's constant and possibly non-invertible metrics.

Leakage, occurring when a qubit undetectably exits the computational subspace, poses a significant challenge for quantum error correction by inducing correlated errors in space and time. These correlations reduce both the code threshold and the effective code distance. To address this challenge, we introduce two zero-overhead walking circuits for the honeycomb Floquet code (hFC), termed the swirling and sliding circuits, which periodically remove leakage while dynamically protecting logical qubits via a schedule of anticommuting two-body measurements. Unlike previous dynamic circuits for the hFC, ours preserves the distance to Pauli errors. We further study their performance under leakage and find that, for two leakage noise models, they may correct as many leakage errors as Pauli errors, improving on the more widely known walking surface code. Through numerical simulation, we observe that finite-error-rate performance under leakage is strongly influenced by entropic effects, producing a pronounced waterfall regime in which the logical error rate decreases much more rapidly with physical error rate than the expected distance-limited scaling. In this regime, even when asymptotic predictions suggest otherwise, the hFC can outperform the same-distance walking surface code in error rates and sub-threshold logical error scaling.