Papers

We study the conservative and dissipative forces on a Newtonian binary interacting with a nonrelativistic ultralight scalar medium. We derive integral expressions for the instantaneous force and the orbit-averaged energy flux generated by wake perturbations of a general background field. For a homogeneous coherent background, the conservative force suffers from a well-known infrared divergence. We trace this divergence to the failure of perturbation theory about a constant scalar state and identify the gravitational Bohr scale at which the homogeneous approximation breaks down. In the long-wavelength regime, the leading dissipation into the medium is quadrupolar, and we obtain an explicit expression for the leading energy flux from an eccentric binary. In the short-wavelength regime, we identify a local ”hard” region in momentum space that produces a universal Coulomb logarithm, while the nonlogarithmic contribution remains orbit dependent. We then promote the background to a stochastic ensemble with a general velocity distribution. Three independent length scales, together with a derived geometric-mean scale, produce eight distinct scale hierarchies, for which we characterize the ensemble-mean response. Finally, we derive the two-point correlation function of the orbit-averaged energy flux and evaluate it in ”coherent-response” regimes. Wave interference produces order-unity density fluctuations, so the flux in an individual realization can differ substantially from its ensemble mean.

We assess the potential for the NASA Explorer SPHEREx mission in adding to the Milky Way census of Wolf-Rayet (WR) stars by extracting 0.75–5micron spectrophotometry of 30 candidate WR stars previously identified using Gaia XP spectra. We confirm an additional Galactic 13 WR (8 WN and 5 WC) stars, identify 3 probable [WC]-type central stars of Planetary Nebulae, plus 2 He emission line sources. SPHEREx has the potential to detect additional relatively isolated WR stars on the far size of the Milky Way by combining Machine Learning techniques with a judicious choice of candidate selection criteria. Indeed SPHEREx is sufficiently sensitive to identify bright Magellanic Cloud Wolf-Rayet stars.

In this paper, we study color kinematics duality and related double-copy and BCJ structures in the cosmological Grassmannian framework. The simplicity of these correlators in this framework makes it possible to test the double copy construction beyond flat space. We show that the four-point gluon correlator admits a representation whose kinematic numerators satisfy an algebraically exact Jacobi identity, $n_s$ + $n_t$ + $n_u$ = 0. Replacing color factors with kinematic numerators gives a graviton correlator that reproduces the correct factorization residues and flat-space limit, which matches with the BCFW result up to terms regular in all exchange channels. In addition, we construct a differential bi-adjoint scalar matrix whose null vector gives rise to a differential BCJ relation among color-ordered gluon correlators. We further explore an s-copy construction for general integer-spin correlators and show that for odd spin the corresponding numerators satisfy a Jacobi identity. Finally, we develop a differential representation for scalar four-point correlators with gluon exchange. These results establish algebraic and differential aspects of graviton and gluon correlators in curved spacetime.

The cancellation of renormalon ambiguities in the operator product expansion is usually formulated as a matching between perturbative Borel ambiguities and prescription dependence in non-perturbative condensates. We give this cancellation a contour realization in the two-dimensional Gross–Neveu model at large $N$. The exact order-$1/N$ fermion self-energy in the massive vacuum is a finite real integral; ambiguities arise only after its large-momentum behavior is decomposed into products of perturbative Wilson coefficients and the matrix elements of the associated operators. A reduced integral isolates the leading renormalon and reveals a two-dimensional contour geometry in which the perturbative series is generated at the boundary corner where the exchanged momentum reaches the external hard scale, while its Stokes discontinuity is controlled by a distinct boundary critical point at zero exchanged momentum. After renormalizing the ultraviolet divergences of the condensate scale integral, we derive its compensating lateral ambiguity without fixing it by matching to the perturbative sector, distinguish the renormalon saddle from the apparent IR Landau pole, and obtain the reduced contour variables from the large-$N$ auxiliary-field path integral. The renormalon saddle is a boundary critical point of this reduced mode-space integral. The perturbative Wilson coefficient and the renormalized condensate thus define complementary relative cycles of the same reduced integrand whose complex tails cancel, reconstructing the original real integration cycle.

When spatial sections of spacetime are topologically nontrivial, the symmetry arguments usually used to choose the Bunch-Davies vacuum are no longer guaranteed to apply. We study this problem for a scalar field on de Sitter spacetime with toroidal spatial sections. The Hadamard conditions, which ensure the same short-distance structure as in flat space, are local and naturally insensitive to the compactification. By contrast, we find that the infrared sector is affected by the toroidal identifications that break global de Sitter invariance, thereby removing the usual argument for discarding zero-mode contributions. We treat the remaining zero-mode freedom as an effective topology-induced source and analyze its effects on an inflationary de Sitter expansion within a semiclassical setting. We show that, for a massless scalar, the zero-mode sector can contain a contribution to the stress-energy tensor with equation of state $P=-\rho$, behaving as a cosmological constant. Thus, the source driving an early-Universe de Sitter phase need not be determined solely by local terms in the action, but can also depend on the infrared structure of the quantum state. Nontrivial topology, as a concrete realization of de Sitter symmetry violation, turns the otherwise-unfixed infrared freedom in vacuum selection to a contribution directly relevant to the semiclassical expansion. This contribution could potentially be applicable to the current accelerated expansion of the Universe.

Pair-density-wave (PDW) superconductors admit half-quantum vortices that can bind Majorana zero modes, but whether such defects are energetically stabilized in microscopic models remains unclear. We address this question in a spinless honeycomb-lattice model with a PDW phase. Microscopic Hartree-Fock calculations determine the superfluid stiffness $\rho$ and PDW relative-phase stiffness $\kappa$, and show that $\kappa/\rho$ approaches unity near the continuous PDW–Dirac-semimetal transition, strongly reducing the long-wavelength cost of vortex fractionalization. Combining this with microscopically extracted Ginzburg–Landau couplings, we show that vortex-core energetics favor fractionalization, producing a strongly enhanced splitting scale and a field-driven full-vortex–to–half-vortex lattice transition. We determine the resulting two-flavor half-vortex lattice structure and construct the associated Majorana lattice. Majorana hybridization produces geometry- and flux-dependent bands with Dirac nodes and zero-energy Fermi lines. Together, these results establish a microscopic route from PDW superconductivity to field-induced half-vortex and Majorana lattices.

We prove that the regularized, stabilized sandwiched Rényi divergence of finite-dimensional quantum channels converges to their regularized relative entropy as the Rényi order tends to one. The key tool is the channel hockey-stick divergence: Gour's Stinespring approximation bound and a Schatten norm estimate amplify an asymptotic bound below one into exponential decay at higher threshold rates. For channel pairs with finite max-relative entropy, known operational connections then give exponential strong converses for parallel and adaptive discrimination, a sharp zero–one testing law, and the subchannel asymptotic equipartition property.

We study non-stabilizerness (magic) in the ground state of $(2+1)$-dimensional $\mathrm{SU}(2)$ Hamiltonian lattice gauge theory with matter, formulated in the dressed-site basis in the hardcore-gluon truncation and restricted to the zero baryon-number sector. Using matrix product states, we compute three facets of magic: the second-order stabilizer Rényi entropy (SRE) $M_2$, its non-local component $M_2^{\rm NL}$, estimated by the SRE $M_2(\tilde\psi)$ of the computational-basis Schmidt representative $\tilde\psi$, and a lower bound in terms of the anti-flatness $F$ of the entanglement spectrum. We also prove a stronger form of the sandwich relation: $-\log_2(1-4F)\le M_2^{\rm NL}\le M_2$; the lower bound rests on a stronger inequality that we obtain for arbitrary Schmidt bases and rank, with an optimal constant. Combined with recent results, this yields the unconditional bracket $\max\{-\log_2(1-4F),\,M_2(\tilde\psi)-4\}\le M_2^{\rm NL}\le M_2(\tilde\psi)$. We emphasize a structural distinction that makes the non-local quantities the physically preferred diagnostics: whereas the full SRE depends on the (non-unique) encoding of the gauge-invariant local Hilbert space into qubits, the non-local magic, its estimate $M_2(\tilde\psi)$, and the anti-flatness are invariant under site-local re-encodings and are therefore intrinsic to the state and bipartition. By varying the gauge coupling for two values of the fermion mass on lattices up to $8\times 8$ ($\sim 288$ qubits) with bond dimension up to $128$, we find that the non-local magic, as estimated by $M_2(\tilde\psi)$, furnishes a sharper and more bond-dimension-friendly probe of the gauge-matter delocalization crossover compared to the full SRE or the gauge-invariant entanglement entropy, retaining a clear signal at bond dimensions well below those needed to converge the ground state itself.

The ${L}_\Lambda w_{1+\infty}$ algebra is the cosmological-constant deformation of the celestial $Lw_{1+\infty}$ symmetry algebra of gravity, providing an infinite-dimensional extension of the three-dimensional conformal algebra at the AdS$_4$ boundary. We construct the Hamiltonian charges associated with these symmetries from boundary ambitwistor-space data, using LeBrun's Heaven on Earth correspondence, and translate them to spacetime via the Penrose transform. We relate the resulting expressions to light-ray operators in the boundary CFT$_3$, including ANEC operators, their conformal descendants, and their commutators.

As part of our quest for high-redshift radio galaxies (HzRGs), we present the results obtained from optical spectroscopy leading to the confirmation of one HzRG candidate selected using the Lyman break technique. The optical emission, associated with the radio source TXS 2354+015, exhibits a prominent drop in the r band, which is characteristic of sources at z~5. The optical spectrum shows a bright asymmetric emission line, with a large equivalent width of ~900 Angstrom. Its identification with the Ly$\alpha$ line is confirmed by the detection of the NV$\lambda$1640, leading to a redshift estimate of z=4.946. Its broadband radio spectrum is well reproduced using a power law with a slope of 0.94. The radio luminosity at the rest frame frequency of 500 MHz is 6.2e29 W/Hz. This result makes this source the most powerful radio galaxy currently known.

One of the primary applications of quantum mechanics is leveraging coherence to make incredibly sensitive measurements. One way to push the boundary of this sensitivity even further is to initially prepare a system in an entangled state, though the fragility of these states limits their usefulness. In this work, we propose a new approach by which sensitivity comes during the preparation of an entangled state, rather than first preparing the state and then using it for sensing. An adiabatic state preparation process is then accelerated by universal counterdiabatic driving, which is more robust to uncertainty in microscopic parameters than standard shortcuts to adiabaticity. Using a toy Lipkin-Meshkov-Glick model, we show that this approach can also lead to sensitivity beyond the Standard Quantum Limit.

Pair-density-wave (PDW) order is a superconducting state with a spatially modulated order parameter and is generally subleading to uniform superconductivity within the conventional weak-coupling BCS paradigm. Here, using a twisted bilayer checkerboard-lattice model, we show that momentum-space nonsymmorphic symmetry provides a generic mechanism that can instead promote PDW order to a leading weak-coupling instability. Combining a controlled Wilsonian renormalization-group analysis at charge neutrality with a finite-doping Bethe–Salpeter analysis, we determine the leading ordering tendencies among various competing instabilities. At charge neutrality, two symmetry-related quadratic band touchings support spin-singlet and spin-triplet PDW instabilities for different attractive interactions, whereas repulsive interactions favor quantum Hall insulators and finite-$\mathbf Q$ density-wave orders. Upon doping, the nonsymmorphic symmetry, together with time reversal, protects the Cooper logarithm in the finite-momentum $\mathbf Q$ pairing channel at generic fillings, allowing PDW order to remain a leading instability over the broad doping range studied. This mechanism extends naturally to moiré Chern bands: by breaking time-reversal symmetry while preserving the nonsymmorphic symmetry and inversion, we obtain a fully gapped topological spin-triplet PDW state with BdG Chern number $C_{\text{BdG}}=8$. Our results establish momentum-space nonsymmorphic symmetry as a general weak-coupling route to PDW superconductivity, including topological PDW states.

Using the fast radio bursts (FRBs) localized with the CHIME/FRB Outriggers between February and October 2025, we conduct a systematic search for dwarf host galaxies. We identify and characterize eight host galaxies for one repeating and seven apparently non-repeating FRBs, newly classifying five hosts as dwarfs with stellar masses of $\log(M_\ast/M_\odot) \approx 8.1-9.1$. We measure their gas-phase metallicities and, when combined with hosts from the literature, finding that repeating and non-repeating FRBs have statistically distinct metallicity distributions. Non-repeating FRBs inhabit relatively metal-rich environments ($12 + \rm log( O/H) \gtrsim 8.35$) across the full observed mass range, while repeating FRBs have significant preference for dwarf galaxies with $12 + \rm log( O/H) \lesssim 8.35$. Almost all repeaters associated with persistent radio sources (PRSs) are in dwarfs with $12 + \rm log( O/H) \lesssim 8.15$, suggesting a connection between repeaters in metal-poor dwarf environments and PRS formation. We also calculate the host dispersion measure ($\mathrm{DM}_{\rm host}$), Faraday rotation measure ($\mathrm{RM}_{\rm host}$), and scattering timescale ($\tau_{\rm host,1~GHz}$) for 54, 37, and 24 FRBs, respectively, including both our sample and FRBs in the literature. We find no significant ($>3\sigma$) correlation between host metallicity and $\mathrm{DM}_{\rm host}$, $\mathrm{RM}_{\rm host}$, or $\tau_{\rm host,1~GHz}$. The distinct host metallicities of apparently non-repeating and repeating FRBs may reflect metallicity-dependent massive-star evolution, with repeating FRBs and PRSs requiring more extreme magnetar birth spins than the predominantly non-repeating FRB population.

The orbital Hall effect (OHE) underpins the emerging field of orbitronics. A complete quantum-mechanical evaluation of the usual orbital-current operator must retain all matrix elements of the position operator, including its band-diagonal differential part; intrinsic calculations based on this full evaluation can differ by orders of magnitude from the conventional truncation. A relaxation-time estimate of disorder effects on the full current has been reported, but a microscopic treatment of extrinsic scattering has remained unavailable because the position operator requires wavevector-off-diagonal density-matrix elements. Here we develop such a theory by solving the quantum kinetic equation for the wavevector-off-diagonal density matrix, including band-structure, field-corrected side-jump, and skew-scattering terms within the non-crossing approximation. We apply the theory to a massive Dirac cone and to the same model with a particle-hole-symmetry-breaking quadratic term. Disorder-generated contributions of order $\tau^0$ are generically comparable to the intrinsic current in the metallic regime and cannot be separated from it by simple disorder-strength scaling. For the bare massive Dirac cone the surviving non-crossing extrinsic correction cancels the conventional intrinsic value and leaves a total current exactly twice the quantum correction, whereas particle-hole asymmetry activates conventional skew scattering for disorder with a non-zero third moment, which dominates in sufficiently clean samples. The wavevector-off-diagonal construction is independent of the two-band model and provides a general route to disorder corrections for position-dependent observables in multiband solids.

We derive thermal Polyakov bootstrap equations for CFT two-point functions on $S^1_\beta\times\mathbb R^{d-1}$. Combining the thermal OPE with a dispersion relation and the method of images gives a representation in blocks that satisfy Kubo-Martin-Schwinger (KMS) covariance. Their local expansions contain additional thermal blocks at classical double-twist dimensions. Matching the prescribed physical OPE requires cancellation of these generated contributions, yielding explicit linear equations for thermal OPE coefficients, including any compensating terms required by the reconstruction. We formulate the equations for bosonic scalar correlators and antiperiodic scalar structures extracted from fermion correlators, and find consistent solutions in analytically tractable thermal CFTs. The coefficient matrix, which we call the Polyakov matrix, connects the large-spin expansion to KMS consistency and reduces to minus the identity on the classical double-twist spectrum. This structure suggests a numerical approach to the three-dimensional Ising CFT that combines known low-lying vacuum data with a self-consistent large-spin tail.

Dark matter halos are expected to form with a prompt $\rho \propto r^{-3/2}$ density cusp at their centres, and in warm dark matter (WDM) cosmologies these cusps can dominate the inner structure of the low-mass subhalos that survive free-streaming suppression. We investigate whether that inner structure is visible to gravitational waves lensed in the wave-optics (WO) regime. Extending the strong lensing WO diffraction-integral framework to WDM, we generate subhalo populations with the sashimi-w semi-analytic model, assign each subhalo a prompt cusp through a cusp-halo relation, and compute the frequency-dependent amplification factor $F(f)$ across the LISA band over 500 independent realisations for $m_\chi \in \{6, 10, 20, 40, 80, 100\}$~keV. The signal is governed by the overlap between the WDM-suppressed subhalo mass function and the $m_{\rm sub} \sim 10^{4}$ to $10^{7}\,M_\odot$ range to which the LISA band is sensitive. Modulations are negligible for $m_\chi \lesssim 10$~keV, reach the percent level by 40~keV, and then saturate to a CDM-like plateau, so the observable reliably probes $m_\chi$ only over the intermediate range $\sim 20$ to 40~keV. Repeating the 40~keV ensemble on identical subhalo catalogs with the cusp amplitude set to zero, we find that prompt cusps enhance the amplitude modulation only mildly. Where cusps make up a large fraction of the subhalo mass, free-streaming has already depleted the signal-carrying mass range, and where that range is populated, the cusp lies well inside the Fresnel radius, to which the diffraction integral responds through projected enclosed mass rather than central density. A percent-level detection in a strongly lensed massive black-hole binary would therefore establish the presence of substructure on Fresnel scales without diagnosing its inner profile.

Despite the importance of reinforcement learning with verifiable rewards (RLVR), the extent to which it can learn new reasoning capabilities remains debated. Here we study the optimization landscape of RLVR on algorithmic tasks, such as iterated group and quasigroup multiplication. To this end, we map entropy-regularized RLVR over myopic tabular policies onto an energy-based (spin-glass) model over deterministic policies. This mapping upper-bounds what RLVR can achieve, and lets us rigorously characterize the landscape in this tabular setting. We show, both theoretically and experimentally, that for a wide class of models and tasks with uncorrelated inputs, this landscape is benign, containing no local minima that could trap RLVR training. Rather, the practical difficulty of these tasks appears to stem, at least in part, from issues such as diffusive barriers and gradient-estimation error in traversing the landscape. These are genuine obstacles that can prevent a solution from being found, but they are distinct from the landscape itself being rugged. We show that these obstacles can often be mitigated through the choice of entropy regulator. Consistent with this theory, we find that a transformer trained from scratch, using only last-token rewards, successfully learns an algorithmic chain of thought for iterated non-Abelian group multiplications.

We study the stimulation of superconductivity by an external electromagnetic radiation — the Eliashberg effect — in a disordered $d$-wave superconductor. We use the self-consistent Keldysh–Nambu quasiclassical formalism augmented by an inelastic relaxation term to compute the static correction to the pairing amplitude by treating elastic scattering within the Born approximation with the full second-order impurity vertex retained. We find that the clean $d$-wave superconductor shows no enhancement at any frequency: without a momentum relaxation channel a spatially uniform drive produces no absorption, only pair breaking. Elastic impurity scattering restores absorption and with it activates response in the Eliashberg channel. We find enhancement of the pairing amplitude which is confined to a bounded window in frequency, disorder strength, inelastic scattering rate and gap magnitude. Both edges of the window scale with the gap but are set by the competition between quasiparticle redistribution, heating and pair breaking rather than by a feature of the spectrum. Furthermore, in contrast to the $s$-wave when the threshold is pinned at $2\Delta$, the upper edge is not pinned to the maximum gap and crosses it as the system approaches the critical temperature.

Universal statistical laws are a hallmark of quantum chaos. A natural route to probing them in quantum circuits is offered by the statistics of output probabilities, which for sufficiently deep chaotic circuits approach the Porter–Thomas distribution. Here we show, theoretically and experimentally, that the approach to this limit is itself universal, and is organized into geometry-dependent universality classes. At shallow depths, growing only logarithmically with system size, we derive universal output-probability distributions that differ markedly from Porter–Thomas statistics, and we test them on a programmable superconducting processor. Two-dimensional lattices, open one-dimensional chains, and periodic rings exhibit distinct universal forms, retaining signatures of geometry and boundary topology within the same logarithmic-depth regime. Within each class, the distributions depend on only two parameters: an effective many-body Thouless scale and the accumulated number of incoherent errors. Extracting these directly from measured bitstring frequencies provides access to scrambling, entanglement growth of the target state, and global fidelity, without reconstructing or classically simulating the implemented circuit. We also identify a noise-driven transition between output statistics retaining system-wide quantum correlations and those reproducible by independent small patches. Our results establish universality beyond Porter–Thomas statistics as a quantitative tool for characterizing noisy quantum processors.

Software underlies virtually all modern research, yet citation practices for software remain inconsistent and often inadequate. This short primer offers practical guidance in two directions: how to properly cite the software you use in your research, and how to make your own software easy for others to cite. We describe a seven-step workflow for citing software and a five-step workflow for making your software citeable, along with the community tools and resources that support each step. We also address common questions about what software, and how much, to cite. The guidance is field-independent; we use astronomy as a running example throughout, and readers in other disciplines can substitute the equivalent tools and conventions of their own community. Our goal is to lower the practical barrier to good software citation practice and to support a research culture in which software receives the scholarly credit it deserves.

We find the shape of the modes containing most of the entanglement between two disjoint spherical regions in the vacuum of a free massless real scalar field in 1+3 Minkowski spacetime. Using the Klco-Beck-Savage method we approximate the field and momentum profile of the vacuum most entangled modes (MEMs) between two spheres. We define a non-perturbative mode-swap entanglement harvesting protocol where two spacelike separated probes extract the maximal amount of entanglement between two regions. Our results are a first step in the direction of understanding the local entanglement structure of quantum fields, pinpointing the exact degrees of freedom that encode accessible quantum correlations in QFT.

We study magnetic systems featuring competing antiferromagnetic and spin-nematic phases, described by dipolar primary and quadrupolar secondary order parameters. While the antiferromagnetic phase breaks time-reversal and spin-rotational symmetries, only spin-rotational symmetry is broken in the nematic phase, which can be understood as a spin-vestigial phase. We develop a classical continuum field theory with independent vector and tensor fields and use mean-field and renormalization group methods to analyze its phases and finite-temperature transitions. We determine the fixed-point structure using an $\epsilon$ expansion about the upper critical dimension of six and track the resulting fixed points to lower dimensions using a perturbative fixed-dimension renormalization group approach. The multicritical fixed point governing the meeting of the antiferromagnetic, nematic, and paramagnetic phases occurs at imaginary coupling throughout the dimensions accessible to our analysis, indicating a first-order transition through the triple point. Away from the triple point, a direct paramagnetic-to-antiferromagnetic transition can be continuous when the tensor mass is sufficiently large, without an intervening vestigial phase. We argue that it is governed by the cubic universality class in $d=3$ and the Ising universality class in $d=2$ for three-component order parameters. The theory also supports a two-step transition with an intermediate spin-nematic phase. The paramagnetic-to-nematic transition is governed by the four-state Potts universality class in $d=2$ and is first order in $d=3$, while the nematic-to-antiferromagnetic transition is generically expected to be continuous and of Ising type in both dimensions. Our results provide a field-theoretical framework for understanding the finite-temperature transitions observed in the candidate Kitaev material Na$_2$Co$_2$TeO$_6$.

We present a tool for the automated generation of Universal Feynman Output (UFO) files, used as input by many event generators, for generic Beyond the Standard Model (BSM) theories and their low-energy Effective Field Theory (EFT) descriptions. Building on the Matchete package, which provides a convenient interface for defining BSM models and deriving their EFTs via matching, our extension obtains the Lagrangian in the broken phase of the electroweak symmetry, performs automatic gauge fixing, computes all Feynman rules, imposes flavor symmetries, and exports the model in the UFO format. This significantly streamlines BSM and EFT analyses, as we illustrate with an EFT validity study for a multi-channel leptoquark search. As an application, we provide a complete UFO implementation of the Standard Model EFT up to dimension eight using the $M_W$ electroweak input scheme.

The Vertical Shear Instability (VSI) is known to create turbulence and strong vertical mixing in protoplanetary disks if thermal relaxation is sufficiently fast. In simulations where this condition is met, VSI can loft dust particles to large altitudes, creating a vertically extended appearance in mock millimeter-wavelength observations, which is inconsistent with the morphology of the majority of observed class II disks. We present simulations of protoplanetary disks with VSI that are consistent with the observed thin-disk geometries, while maintaining the commonly observed bowl-shaped morphology in scattered light images at micrometer wavelength. We show that this outcome arises naturally when the effects of dust coagulation, sedimentation, and dust-gas thermal accommodation are taken into account. Sedimentation-driven coagulation removes large amounts of dust from the disk atmosphere, in the process slowing down the dust-driven cooling of the gas. At the same time, a dense midplane layer of millimeter-sized grains forms, which exerts aerodynamic drag on the gas. This results in the termination of the VSI's corrugation mode. Only weak VSI activity remains in the upper and lower hemispheres. These processes occur on the typical dust growth timescale and suppress strong VSI-induced turbulence within a few hundred thousand years. VSI could thus generally be restricted to the class I evolutionary stages of protoplanetary disks.

Given any four-point scattering amplitude, the partial wave expansion resolves its residues into contributions from individual spin-$s$ internal particles, where each coefficient in the expansion must be non-negative due to unitarity. In this paper, we introduce a partial wave expansion tailored to "half-ladder" Feynman diagrams at generic multiplicity $n$ and spacetime dimension $d$, with scalars on the external lines and massive bosons on the internal lines. Derived by enforcing that the half-ladder residue is an eigenvector of each internal particle's angular momentum, each basis element is written as a sum over products of special functions related to the Gegenbauer polynomials, and we give a set of simple graphical rules for explicitly building any of them at generic $n$ and $d$. To put the formalism into action, we apply our technology to deformed "string theory" amplitudes, defined a priori in any dimension and with arbitrary Regge intercept $\alpha_0$. By constructing multi-particle positivity tests for the half-ladder and "twisted" half-ladder diagrams, we show that the only amplitudes consistent with unitarity are those belonging to open bosonic string theory ($\alpha_0 = -1$) below the critical dimension $d \leq 26$. In particular, we confirm the expectation that $Z$-theory ($\alpha_0 = 0$) fails to define a unitary string theory at finite $\alpha'$. Finally, for the open bosonic string, we show how knowledge of the partial wave expansion allows us to resolve the contributions to the half-ladder diagram from individual degenerate states at a given level and spin in a purely on-shell analysis. As an ancillary file, we include a Python script combwaves.py which implements the partial wave expansion given a generic (twisted) half-ladder residue.

We explore a model of sneutrino dark matter in the Minimal R-Symmetric Supersymmetric Standard Model as a possible interpretation of the recent event at LZ with high reconstructed nuclear recoil energy. This model can account for a small mass splitting between neutral states in a technically natural way due to terms which violate lepton number; importantly, this model does not require sending masses of other SUSY particles to an extremely high energy scale. Mixing with additional singlet states suppresses the interaction strength with standard model particles which allows the sneutrino to sneak around bounds obtained from searches for annihilating dark matter captured in the sun and allows the dark matter to remain near the TeV scale and still potentially account for the LZ results. Typically, sneutrinos produce too much dark matter unless they are very light; we consider a minimal scenario where production occurs through co-annihilation as one specific mechanism where the dark matter mass is predicted to be around $360$ GeV.

In this time of uncertainty and anxiety about artificial intelligence (AI) in mathematics, it is more important than ever to reflect on how we can prepare the next generation of professional research mathematicians for success. In this article, we report, from the front lines, on how AI is affecting the long-running summer math research program at the University of Minnesota Duluth. We do not have answers to the big questions about how AI will shape the future of professional mathematics. What we can offer instead is a collection of practical observations about the role of AI in the 2026 Duluth program: how it was useful and where it presented potential pitfalls; what surprised us and how we adapted; and what went well and what we wished we had done differently.

Autonomous research agents can design experiments, evaluate results, and write reports, giving them control over both a scientific result and the evidence used to support it. This creates a risk of reward hacking: meeting the reward criteria without achieving the intended goal. We study (1) how often models reward-hack without instructions to do so, (2) how effective and detectable their methods are when hacking is allowed, and (3) how they adapt when an LLM review panel returns its decision and reasons. Across 17 language models and 38 tasks, the spontaneous reward-hacking rate is 30.5% on open-ended research-pipeline tasks and 2.9% on task-specific kernels. When hacking is allowed on tasks whose pass thresholds exceed our best compliant baselines, 505/677 attempts (74.6%) are confirmed reward hacks: they both clear the threshold and receive mechanism-verification panel confirmation of an evaluation exploit. An LLM panel reviewing only submitted code and reported scores misses 33/505 confirmed hacks (6.5%). Direct methods that achieve the highest scores are often easy to detect, while less direct methods evade more often. In a five-round loop, the number of model-task pairs with an evasion rises from 7 to 56. Among 79 pairs evaluated under two feedback conditions, cumulative evasion reaches 40.5% with detailed feedback and 20.3% with generic rejection. The detailed condition includes the review decision, reasons, and attempt history, so this comparison does not isolate the effect of explanations. These findings highlight the need for stronger defenses, including metrics kept outside the agent's control and independent recomputation on data chosen to expose likely exploits.

Most studies of prompt injection focus on generative agents, leaving their effects on models with schema-defined outputs unclear. We examine these effects in Jev, a non-generative decision model, using 510 reconstructed InjecAgent cases. Malicious content shifts action probabilities but rarely causes Jev to select the attacker's target. Override markers reduce this influence, while claims of contextual relatedness have small effects. Adaptive attacks using score feedback double the mean highest attacker-target probability found during optimization, while success on fresh validation calls rises from 1.8% to 3.5%. Exploratory analysis links these successes to small initial decision margins or greater attacker control over the observation. Together, these findings show that schema-defined outputs change but do not eliminate prompt-injection risk, highlighting the need to evaluate how untrusted content influences choices within the allowed action set.

Recently, the presence of a large shell around the bright hypergiant star VY Canis Majoris was reported in images from the SPHEREx mission. We report similar shells around other bright stars at similar locations in SPHEREx Band 1 images, suggesting that this shell is not astrophysical but instead a detector artifact.