高能暂现源 · 每日文献综述

We investigate gravitational-wave birefringence from pulsar glitches in dynamical Chern-Simons gravity, a parity-violating extension of general relativity motivated by string theory and quantum gravity. Using the Hartle-Thorne slow-rotation formalism to model the neutron star background and a perturbative treatment of the Chern-Simons coupling $α$, we derive the modified Regge-Wheeler equation governing axial gravitational perturbations and compute the resulting polarization-dependent phase shift. The parity-violating interaction produces a fractional asymmetry between the right- and left-handed circular polarizations that scales as $α^2$, grows linearly with frequency and rotation rate, and falls as the fourth power of the stellar radius. For a constant-density interior model we obtain an analytic matching solution for the background scalar field, and show that the interior curvature factor vanishes identically--the interior solution being conformally flat--so that the scalar dipole is sourced entirely in the vacuum exterior; a centrally condensed equation of state can only increase the predicted signal. We further extend the analysis to a two-fluid model incorporating differential rotation between the neutron superfluid and the charged component; the resulting correction is bounded by the fractional glitch amplitude itself, independently of the equation of state, and is negligible for typical pulsar parameters. Expressed through the dimensionless quantity $ζ\propto αM/R_\star^{3}$, the fractional polarization asymmetry for millisecond pulsars reaches $\sim 10^{-3}$ at kilohertz frequencies at the boundary of perturbative validity--four orders of magnitude above the sensitivity of current ground-based detectors to strain asymmetries. However, realistic bounds are limited by the signal-to-noise ratio required to resolve a ratio of gravitational wave strains.

We present a study of the interstellar molecular and atomic gas associated with the gamma-ray supernova remnant (SNR) RCW 103 using the Mopra $^{12}$CO($J$ = 1-0), Atacama Large Millimeter/submillimeter Array $^{13}$CO($J$ = 1-0), and the Australia Telescope Compact Array & Parkes HI data. We find that CO clouds in the velocity range from $-58.7$ to $-43.5$ km s$^{-1}$ show a clear spatial correspondence with the X-ray shell of RCW 103, particularly extending from the northwestern to the southeastern regions. The position-velocity diagram of the CO emission reveals an expanding gas motion with an expansion velocity of $ΔV \sim$$7.5$ km s$^{-1}$, which may have been produced by shock waves and/or stellar winds from the progenitor system. We estimate the total cosmic-ray proton energy to be $\sim$$3.0^{+1.1}_{-0.6} \times 10^{47}$ erg, adopting a total interstellar proton density of $\sim$$810$ cm$^{-3}$. This value is more than an order of magnitude lower than those inferred for other gamma-ray emitting SNRs of comparable age. Given that a substantial amount of dense gas remains within the interior of RCW 103, where thermal X-ray emission is dominant, our results possibly suggest that a significant fraction of the shock energy may have been converted into thermal energy via gas heating, thereby reducing the efficiency of cosmic-ray acceleration.

Stellar flares have long served as stellar magnetic activity tracers. The flare waiting time, defined as the interval between two consecutive flares, provides a valuable diagnostic for probing underlying mechanisms of energy storage and release in stellar atmospheres. In this work, utilizing flaring M dwarfs observed by the Kepler satellite, we establish a simple yet effective activity proxy, i.e., median flare waiting time ($t_{\rm{w, med}}$). Our results show that the $t_{\rm{w, med}}$ can trace long-term activity levels similar to the flare rate. However, $t_{\rm{w, med}}$ corresponding to different waiting time percentiles may encode richer physical insights than flare rate. In addition, for the first time we construct a clear relation between $t_{\rm{w, med}}$ and stellar rotation period, which is quite similar to the canonical activity--rotation relation. More intriguingly, this relation exhibits a more notable supersaturation effect (i.e., below a critical rotation period, $t_{\rm{w, med}}$ begins to increase instead of keeping constant) compared to other activity proxies. The filling factor--rotation period relation favors poleward migration of active regions as the explanation for supersaturation, rather than coronal stripping. With the dramatic increase in stellar flares detected by missions like TESS and the upcoming Earth 2.0 satellite, $t_{\rm{w, med}}$ will become a powerful diagnostic for probing stellar magnetic activity and underlying physics.

Cluster emission is an exotic decay mode between alpha-decay and fission, in which a parent nucleus emits a cluster of nucleons heavier than an alpha-particle, but lighter than what is usually considered a fission fragment. The properties of cluster emission were investigated by analyzing five high-energy events detected in a spectrum of a mass A = 230 beam produced at ISOLDE (CERN). Under the assumption that these events were caused by cluster emission, the most likely parent-cluster pair responsible for the five high-energy events, was found to be $^{230}Ra$ emitting $^{22}O$, with a branching ratio of $(4.3$ +\- $1.9) \times 10^{-9}$. Four analytical formulas were used to estimate the partial half-lives of cluster emission for a group of neutron-rich nuclei. The decay rate was calculated for a selection of cluster nuclei for each parent isotope. The rates of all decay channels of cluster emission per parent nucleus were then included in calculations of the r-process nucleosynthesis in a neutron star merger in order to study the possible impact of cluster emission on the r-process nuclear production. The resulting isotopic abundance distributions were compared to those calculated for a case in which cluster emission was not considered. It was found that the inclusion of cluster emission decay rates from the simple analytical formulas available nowadays does not influence the results of the r-process nucleosynthesis.

Geomagnetic indices quantify the disturbances of the Earth's magnetic field. They form an important space weather and space climate record that tracks solar-terrestrial effects on the Earth. Although geomagnetic activity has its origins in the photospheric magnetic fields, the two have so far not been directly related. We develop a photospheric index of geomagnetic activity (PIGA) that measures the geoeffectiveness of the Sun based on photospheric magnetic field as presented in synoptic magnetograms. We construct the PIGA from synoptic magnetograms of Wilcox Solar Observatory by relating the equatorial and axial dipole components of the photospheric magnetic field to the geomagnetic Kp index by means of linear regression. PIGA captures well the evolution of the geomagnetic activity over the past 50 years at the resolution of solar rotation, having linear correlation of 0.71 with the geomagnetic Kp index. PIGA makes it straightforward to measure the geoeffectiveness of the Sun using the dominant structure of the global solar magnetic field. PIGA can for example be used to study the geoeffectiveness of the Sun using historical simulations or the predicted solar magnetic field. PIGA also allows to use observed geomagnetic activity as an additional constraint for the early magnetic field reconstructions of the solar magnetic field.

Axion-like particles (ALPs) bound in the solar gravitational field form the so-called ALP solar-basin. Since the two-photon decay of non-relativistic particles is approximately isotropic, this population can be searched for using observations in the anti-solar direction. In this work, we propose a search strategy for narrow decay-line signals from the ALP solar basin using \textit{XMM-Newton} blank-sky observations (XMM-BSOs) stacked spectra data taken in directions opposite to the Sun. By jointly fitting the signal and background model, we obtain limits on $g_{aγγ}^{95}$ in the mass range $m_a=1.4\text{--}16~{\rm keV}$, with typical sensitivities of $g_{aγγ}\sim10^{-10}\text{--}10^{-11}~{\rm GeV}^{-1}$. We have implemented the first anti-solar search for the solar basin, demonstrating that this strategy can exploit the stacked exposure of a large number of X-ray observations and provide a scalable analysis framework for future searches.

We reconstruct the late-time expansion history and dark energy dynamics using available cosmological data. We consider dark energy as an effective minimally coupled canonical scalar field without assuming a specific form for its potential. For reconstruction, we use Gaussian Process (GP) regression with a joint dataset consisting of 32 cosmic chronometer measurements (CC32), DESI DR2 baryon acoustic oscillation data, three Type Ia supernova compilations (Pantheon+, Union3, and DES Y5), and compressed CMB distance priors. From the reconstructed Hubble parameter and its derivatives, we obtain the scalar field kinetic and potential energy densities, the equation of state $w(z)$, and the dimensionless slope parameter $|λ(z)|$ of the scalar field potential. The reconstructed potential is nearly constant at late times, with $\tilde{V}(0) \simeq 0.68$--$0.70ρ_{c,0}$, close to the dark-energy density in flat $Λ$CDM. The kinetic term remains small over $0 \lesssim z \lesssim 2.5$, while $w(z)$ remains close to $-1$ within the uncertainties. Using only CC32 and DESI DR2, we find a mild ($1σ$) hint of a phantom-divide crossing near $z\sim0.5$. However, this feature depends on the supernova compilation and cannot be regarded as a firm conclusion with current data. The slope parameter $|λ(z)|$ shows mild evolution, with present-day central values around $0.8$--$1.0$, but large uncertainties due to its dependence on higher-order derivatives of the expansion history. The curvature parameter $Γ(z)$ requires even higher-order derivatives and is not constrained by current data; therefore, we do not report its reconstruction. Allowing small spatial curvature, $Ω_k=0$--$0.02$, changes the results only slightly and remains within the existing uncertainty bands.

We analyze quasi-periodic oscillation (QPO) candidates in the 0.1-100 GeV Fermi-LAT light curves of 41 TeV blazars with a conservative, multi-stage timing analysis combining FFT power spectra, red-noise Monte Carlo simulations, false-discovery-rate control, bootstrap period uncertainties, split-sample cross-validation, Lomb-Scargle and WWZ diagnostics, a detection-only sensitivity test, observingwindow checks, and injection-recovery experiments. A first-pass FFT search flags 16 sources above a raw 95% level and six strong candidates after data-quality and split-sample filtering, none satisfying a fully automated multi-method tier. Source-level significance is quantified with forward simulations that preserve the observed sampling and non-Gaussian flux distribution, treating the red-noise slope as a nuisance parameter and never assuming independence of periodogram powers. At simulation-selected reference slopes, only two of the 41 sources reach a raw source-level p < 0.05,only J1555.7+1111 reaches Benjamini-Hochberg significance, and none survives Benjamini-Yekutieli correction. No source survives the conservative slope envelope. Even J1555.7+1111, which the first-pass search flagged as significant but did not rank among its strong candidates, has a global p-value that rises to 7.5*10^-3 under the least favorable slope (~ 0.31 after the full-sample search), so it is not a robust detection. With substantial period uncertainties, a weak observing window, and limited injection-recovery sensitivity, we establish no robust gamma-ray QPO detection: a sampling-faithful red-noise treatment substantially contracts the candidate set, providing an explicit assessment of QPO credibility.

As a new probe, the weak lensing (WL) peak steepness statistics is sensitive to the density profile of halos that encodes important information of baryonic feedback and dark matter properties, leading to a promising means to statistically constrain these effects using WL data. In this article, we present its first application to HSC Y1 data to demonstrate the great potential of this new statistics. Within the phenomenological framework of HMcode2016 that attributes the baryonic feedback solely to the reduction of the halo concentration parameter and focusing on high peaks originated dominantly from massive clusters, our analyses by combining WL peak height and steepness statistics resulted in $S_8=0.76^{+0.08}_{-0.07}$ with the maximum-a-posteriori (MAP) of $0.79$ and low concentrations. Taking the form of the concentration-mass relation as $c(M,z)=A(1+z_{\rm f})/(1+z)$ with $z_{\rm f}$ being the formation redshift of halos with mass $M$ at redshift $z$, we obtain $A=1.93^{+1.33}_{-1.16}$ (MAP=$1.70$) in comparison with $A=3.34^{+1.52}_{-1.74}$ (MAP=3.31) from dark matter only simulated mocks. The result tends to support phenomenologically strong baryonic feedback effects at cluster scales.

High-resolution observations are expected to probe the photon-ring structure of black hole images. The critical curve is the geometrically defined limiting locus on the observer's screen toward which successive higher-order photon subrings accumulate. For a Kerr black hole, the shape of this limiting curve depends on the dimensionless spin parameter $a$ and the inclination angle $i$. We investigate whether these parameters can be determined uniquely from a connected positive-length segment of the critical curve when the segment's position and orientation on the screen are unknown. For rotating non-extremal Kerr black holes with $0<a<1$ and $0<i\leqπ/2$, we prove global uniqueness: if two such segments coincide as point sets after orientation-preserving rigid motions of the screen, then their Kerr parameter pairs coincide and the two rigid motions are identical. The proof uses an irreducible implicit polynomial to show that a segment determines the full algebraic curve containing it and then recovers the parameter pair $(a,i)$ from rigid-motion invariants of that polynomial.

Einstein-Boltzmann codes compute theoretical predictions of cosmological models and rely heavily on interpolation in their independent variables: time $τ$, wavenumber $k$ and multipole $\ell$. We give a practical summary of interpolation with Chebyshev polynomials, which converges rapidly for smooth functions and thus pairs naturally with approximation-free Einstein-Boltzmann codes. By solving the perturbations and line-of-sight integrals at Chebyshev nodes in $k$ and $\ell$, we show that Chebyshev polynomials interpolate to higher precision than traditional cubic splines from fewer explicit solutions. On a set of example spectra for matter and the cosmic microwave background (CMB), we find up to four orders of magnitude lower interpolation error using the same number of points. For a typical CMB temperature spectrum computed with interpolation in both $k$ and $\ell$, Chebyshev polynomials converge to $10^{-4}$-$10^{-5}$ relative error with only 50-80 points per variable, while cubic splines approach $10^{-4}$ error with 200 points, translating to a $2.5\times$-$4\times$ speedup. The exact improvement depends on the target function and is generally more dramatic at high precision levels. Standard Chebyshev $\ell$-interpolation needs line-of-sight integrals generalized to non-integer $\ell$, but we show a way to avoid this by rounding the nodes to integers. Chebyshev interpolation is implemented in SymBoltz, which is available at https://github.com/hersle/SymBoltz.jl.

Gaussian Processes (GPs) are widely used to model stellar variability in photometric surveys, but a statistically successful fit does not guarantee that the inferred hyperparameters correspond to physically meaningful stellar properties. This is especially important for young, active stars, whose TESS light curves contain evolving spots, harmonics, and non-sinusoidal variability. We use stellar rotation as a case study to examine when the period hyperparameter of a quasi-periodic GP can be interpreted as a physical rotation period. We introduce a regularized GP likelihood that reweights the covariance-complexity term in the marginal likelihood, reducing the tendency of unconstrained models to converge toward preferred but misleading solutions. We test this framework on 539 stars in IC2602, the Tucana-Horologium Association, Pisces-Eridanus, and GroupX, with independently reported and manually vetted rotation periods. This benchmark evaluates GP hyperparameter interpretability and automated period recovery with minimal human intervention. We compare regularized and unregularized models across several regularization strengths, λ, using literature agreement, sector-level diagnostics, and consistency across TESS sectors. Relative to the standard GP likelihood, regularization improves successful rotation-period recovery by an average of 7%. It also generally reduces the median absolute fractional deviation of periods across sectors, showing that the improvement is not limited to catastrophic failures but also mitigates smaller systematic errors. Regularization is particularly beneficial for non-sinusoidal or evolving modulation, where unconstrained GPs may recover harmonics or spurious timescales. These results show that likelihood regularization and cross-sector consistency are practical diagnostics for assessing when GP-based rotation periods are robust.

The density profiles of Dark matter (DM) halos carry imprints of the DM nature and may be constrained through the lensing effects on gravitational waves (GWs) arising from the halo gravitational potential. In this paper, we investigate GW lensing by two representative types of halo density profiles, i.e., the generalized Navarro-Frenk-White (gNFW) density profile and the Einasto density profile. Using the gravitational lensing equation, we first examine the parameter-space distribution and the imaging characteristics of both profiles under strong lensing, partitioning the parameter space into distinct regions according to the Morse index. We then conduct a detailed analysis of the modulus $\big|F\big|$ and phase $\mathrm{Arg}(F)$ of the amplification factor $F\big(w,y\big)$ at low frequency regime. Our results show that, for a fixed lens mass , increasing the gNFW slope $γ$ leads to a larger amplitude, more rapid oscillation, and more distinct wave-packet morphology in the multiple-image regime. Compared with the gNFW case, the Einasto case (with $α=0.16$, $y<0.6$, and the same $M_{200}$) produces a stronger lensing effect. Notably, the evolution of $F\big(w,y\big)$ with frequency for the Einasto profile differs from that of the gNFW case, making its behavior particularly distinctive.

Claims that neutron stars constrain quantum gravity often compare models whose new physics enters at different stages of the stellar calculation. This article provides a critical map of three representative routes: generalized-uncertainty-principle deformations of microscopic thermodynamics, noncommutative-inspired modifications of the effective matter source, and higher-curvature or scale-dependent modifications of gravitational dynamics. For each route we identify the quantity that is changed, the dimensionless parameter controlling the correction, and the consistency conditions required before masses, radii, or tidal deformabilities can be compared with observations. A deformed phase-space measure must be implemented in a thermodynamically consistent interacting equation of state; a smeared density must be completed by a conserved stress-energy tensor, generally with anisotropic pressure; and a gravitational effective field theory must remain within its derivative expansion and be applied consistently to both the equilibrium background and tidal perturbations. Simple scale estimates show that coefficients of natural size suppressed by the Planck scale are negligible in neutron stars. Observable effects consequently constrain model-dependent effective scales or enhanced couplings, not microscopic quantum gravity directly. The resulting framework is intended to separate robust phenomenological bounds from effects produced by inconsistent implementations or by degeneracies with the high-density equation of state.

Dynamical friction (DF) and gravitational-wave radiation reaction are conventionally treated as distinct dissipative mechanisms acting on a compact object inspiraling through a matter environment. We show that DF is not an additional force but a term at order $q^2ε$ in a covariant multiparameter expansion of the MiSaTaQuWa force equation that vanishes identically in the absence of a particle-induced perturbation of the environment's matter fields, where $q$ and $ε$ are respectively the mass ratio and environmental parameter. As a test of this identification, we evaluate it in weak and adiabatic field limit, solving the polar density-perturbation equation sourced by the particle's own order-$q$ vacuum metric perturbation. The reduction reproduces the Chandrasekhar-Ostriker drag force and predicts two features in strong gravity. An $\ell$-dependent splitting of the wake that only recombines into the Ostriker source in the weak-field limit, and a purely relativistic axial contribution with no Newtonian counterpart. This establishes dynamical friction as an intrinsic piece of self-force theory, extending to strong field, generic orbits, and generic environments, with direct consequences for extreme-mass-ratio inspiral waveform modeling for LISA.

The degree of spin-orbit alignment in merging binary black holes is a powerful probe of their formation history. However, assessing spin-orbit misalignment in individual events remains challenging, as most parameters commonly used to characterize precession tend to be poorly constrained. In this work, we introduce the angle between the total spin and orbital angular momentum, $θ_{LS}$, as an alternative measure of spin-orbit misalignment. Using synthetic observations, we show that $θ_{LS}$ retains more information to discriminate between the aligned- and isotropic-spins hypotheses than commonly used alternatives, including $χ_{\rm p}$. We then study this parameter on binary black hole mergers observed in the GWTC-5.0 catalog, and identify multiple events that are inconsistent with having aligned spins. In several of the examples, this is not apparent from the posteriors of the effective spin parameters ($χ_{\rm p}$ or $χ_\mathrm{eff}$) nor the individual spin-orbit tilts ($θ_1$ or $θ_2$) alone, and thus had not been previously identified. We highlight GW241127_061008 as a compelling case for a second-generation merger, with a large, tilted primary spin, unequal masses, and a large primary mass near the pair-instability-supernova mass gap. Finally, we perform population inference using $\cosθ_{LS}$ and confirm that the binary black hole population is inconsistent with either a purely isotropic or a purely aligned-spin distribution, and requires substantial spin-orbit misalignment.

Scalar-induced gravitational waves (SIGWs) probe primordial curvature perturbations beyond cosmic microwave background scales, while primordial black holes (PBHs) probe the rare tail of the same statistics. We ask whether SIGW morphology can reveal scale-dependent primordial non-Gaussianity. We adopt the internal-leg separable kernel, $\mathcal F_{\rm NL}(\mathbf{k};\mathbf q,\mathbf{k}-\mathbf q)=\bar f_{\rm NL}f(q)f(\lVert\mathbf{k}-\mathbf q\rVert)$, and compute Gaussian, reducible, and connected contributions through $\mathcal O(\bar f_{\rm NL}^{4})$ in the second-order tensor-source approximation. We compare persistent power-law running with a UV-matched tanh weight that saturates at high scalar momentum. Across narrow, finite-width, asymmetric, and multi-slope scalar spectra, running produces more than amplitude renormalization. Power-law running shifts the scalar support sampled by non-Gaussian convolutions, generating peak displacement, asymmetric shoulders, and persistent ultraviolet deformations. The tanh template instead produces a transition-localized modification and approaches a momentum-independent ultraviolet plateau. Analytical estimates of the perturbative hierarchy, weighted slopes, log-normal saddle displacement, peak shifts, and infrared scaling explain the trends and the enhanced sensitivity of the $\mathcal O(\bar f_{\rm NL}^{4})$ sector. For localized sources with finite weighted moments, smooth running preserves the leading $k^3$ infrared behavior up to logarithmic corrections. Since the PBH mass scales as $M\propto k^{-2}$, preferential high-wavenumber weighting conditionally favors lower PBH masses in a narrow-support interpretation. Peak position, spectral curvature, shoulders, and ultraviolet slopes are diagnostics for PTAs, LISA, and third-generation observations, providing finite-order templates for primordial interactions on otherwise inaccessible scales.

Gravitational radiation from the early Universe offers a powerful window into physics beyond the Standard Model (SM). In this paper, we investigate the low-frequency tail of gravitational radiation produced by the decay of massive particles. The energy spectrum in Minkowski space-time, $dE/dω$, exhibits two distinct regimes: a flat bremsstrahlung/memory spectrum above a characteristic cutoff frequency and a suppressed quadratic frequency scaling below it. The typical rate of momentum change of the decay products determines the cutoff frequency. We calculate $dE/dω$ in two independent but equivalent ways: using Weinberg's soft graviton theorem and by solving the equations of motion. This further establishes the connection between Weinberg's soft graviton theorem and the gravitational memory formula, which we generalize to incorporate the effects of final-state scattering. We then compute the cosmological stochastic gravitational energy spectrum $h^2Ω_{\rm h}$ generated by many particle decays in the early Universe. We demonstrate that $h^2Ω_{\rm h}$ exhibits linear scaling with frequency above the cutoff frequency, transitioning to cubic scaling below it. Our findings reveal a suppression of the low-frequency signal compared to a naive linear extrapolation, with implications for detection prospects. Our results also reveal the connection between bremsstrahlung and memory in the context of cosmological stochastic gravitational energy spectra from the early Universe.

The Nancy Grace Roman Space Telescope (Roman), NASA's next flagship mission, is currently scheduled to launch in August 2026. As part of its mission, Roman will conduct the Galactic Bulge Time Domain Survey (GBTDS), which will generate $\sim$50,000 epochs of high-precision photometric and astrometric data for $\sim 10^8$ sources across 1.7 deg$^2$ in the Galactic Bulge. Roman GBTDS astrometry is comparable to Gaia Data Release 4 in terms of number of stars, astrometric precision, and time baseline, and is highly complementary in terms of wavelength and sky location. In this paper, we investigate a synthetic population of GBTDS sources to characterize the detectability of unresolved astrometric binaries, in particular those with compact object companions. Assuming the occurrence rate of black holes (BHs) and neutron stars (NSs) in AU-scale orbits around stars is $10^{-7}$ and $10^{-6}$, respectively, and that Roman achieves an astrometric precision of $1\%$ of a pixel, $\mathcal{O}(10)$ BH+star and $\mathcal{O}(10)$ NS+star detached binaries will be detectable. The BHs will have median mass measurement uncertainties of $\sim 25\%$, increasing the existing sample of detached astrometric BH binaries by a factor of three. Together with the $\mathcal{O}(10^2)$ isolated BHs expected to be discovered by microlensing in the Roman GBTDS and an additional $\mathcal{O}(10)$ detached BH binaries in Gaia DR4, this will provide a representative view of the quiescent Galactic stellar-mass BH population.

Magnetic activity in cool stars governs the thermal structure and high-energy radiative output of their outer atmospheres, thereby influencing stellar evolution, stellar winds, and the atmospheres of orbiting (exo)planets. A long-standing question in stellar astrophysics is whether the mechanisms responsible for atmospheric heating are universal across the Sun and cool stars spanning different ages and activity levels. In this paper, we review recent progress in understanding magnetic heating across the Sun and Sun-like stars through empirical scaling relations between photospheric magnetic flux and radiative output from the chromosphere, transition region, and corona. Analysis of more than a decade of Sun-as-a-star observations reveals that irradiance and magnetic flux follow power-law relationships over a wide temperature range. While coronal emissions exhibit superlinear scaling with magnetic flux, chromospheric and transition-region diagnostics show weaker, sublinear dependencies. Remarkably, observations of G-type stars with ages ranging from 50 Myr to 4.5 Gyr are found to lie on extensions of the solar scaling laws, suggesting that a common magnetic-heating mechanism operates across different levels of stellar activity. We further discuss the application of these scaling relations to reconstructing stellar X-ray and ultraviolet (XUV) spectra from observed magnetic fluxes. The resulting synthetic spectra reproduce actual observations of young, active solar analogs, providing a practical tool to estimate the ionizing radiation whose extreme UV emissions cannot be directly measured.

We present a new statistical technique that utilizes the synergy of precision astrometry and time series photometry from modern space missions to obtain reliable physical parameters of pulsating white dwarf stars. We compute a grid of white dwarf structural models that span the helium-atmosphere pulsating white dwarf (DBV) instability strip, showing that mean period spacings between adjacent pulsation modes and absolute magnitudes derived from Gaia astrometry vary monotonically and in opposing directions across parameter space. While most efforts in white dwarf asteroseismology to directly fit individual pulsation periods to stellar models result in degenerate and poorly resolved solutions, the "seismic technique" produces unique and reliable seismic solutions for global parameters of mass and effective temperature when a reliable mean period spacing is detected. Our models sample various physically plausible interior chemical composition profiles based on modern evolutionary models to propagate uncertainty from the precise structures of actual stars. Once the global stellar parameters are tightly constrained, seismically resolving white dwarf interior structures becomes more computationally tractable, and degeneracies can be resolved. This new seismic technique is largely insensitive to the precise absorption line profiles interpreted by the widely used spectroscopic technique, and therefore provides complementary constraints that can be used to test spectroscopic methods. We demonstrate the method for the pulsating helium-atmosphere white dwarf WD 0158-160 observed by TESS, obtaining $T_\mathrm{eff} = 24584\pm 971$ K and $M_\star = 0.608\pm 0.013$ $M_\odot$.

We present a user-friendly web application that allows users to extract, visualise, and interact with TESS light-curves (LCs) from full-frame images. The service works for any star identified in the TIC database or using stellar coordinates. TESSExtractor performs photometry using a circular aperture centred on the object of interest, with a sky annulus for background subtraction. The application also uses cotrending basis vectors provided by TESS to correct photometry for possible systematic effects. The Lomb-Scargle Periodogram is used to obtain periods from the LC. The TESSExtractor's website offers an intuitive interface for scientific and educational purposes and is publicly available at https://www.tessextractor.app. This tool simplifies the analysis of TESS LCs and provides valuable insight for researchers seeking to extract and analyse TESS data for a variety of scientific goals (e.g., stellar rotation, activity, transient detection, asteroseismology, or simply detection of planetary transits or stellar companions, among others). TESS provides insights into our understanding of time-domain astrophysics, and TESSExtractor provides an accessible and intuitive web interface for scientific exploitation of TESS data. We demonstrate that TESSExtractor produces LCs of sufficient quality for both quantitative and morphological time-domain studies.

The energy-dependent apparent extension of the Crab Nebula, inconsistent with leptonic expectations between X-ray and TeV bands, can be naturally accounted for in a spatially resolved lepto-hadronic scenario. A transition to hadronic dominance occurs above $\sim400~\rm{TeV}$, with nearly pure hadronic emission at PeV energies. The predicted $68\%$ containment radius at $1.1~\rm{PeV}$ is $\sim1.4'$, about six times larger than in purely leptonic scenarios. This energy-dependent extension therefore provides an observable discriminator of the origin of ultra-high-energy gamma-ray emission in pulsar wind nebulae, offering a diagnostic of PeV hadronic acceleration in the Crab Nebula.

Scalar-induced gravitational waves (SIGWs) provide a direct probe of the enhanced primordial curvature perturbations that may also produce primordial black holes (PBHs). We forecast the capability of the space-based gravitational-wave observatory Taiji to search for an SIGW background generated by a broken-power-law curvature spectrum. A signal-injection study is used to validate the analysis pipeline, after which a pure-noise realization is employed to derive projected upper limits on the curvature-spectrum parameters. We translate these limits into constraints on the PBH abundance using the nonlinear compaction function, critical collapse, and the joint Gaussian distribution of the compaction amplitude and curvature at its peak. The resulting projected $95\%$ upper limits on the PBH dark-matter fraction satisfy $f_{\mathrm{PBH}}^{95\%}<1$ over PBH masses from approximately $6.2\times10^{-18}\,M_\odot$ to $1.7\times10^{-8}\,M_\odot$. We compare the forecast with representative Hawking-evaporation and microlensing bounds. In part of the asteroid-mass interval, the projected Taiji limit is more restrictive than the current Subaru Hyper Suprime-Cam (HSC) microlensing constraint.

The study of black hole perturbations typically follows two main approaches: the direct perturbation of the metric, or the perturbation within the Newman-Penrose (NP) or Geroch-Held-Penrose (GHP) formalism. In the latter case, a reconstruction procedure, such as the Chrzanowski-Cohen-Kegeles (CCK) method based on the Debye (Hertz) potential, is required to obtain the corresponding metric perturbation. However, the reconstructed metric is then expressed in the radiation gauge, which is not always optimal. In this paper, we analyze stationary and axially symmetric perturbations of the Kerr black hole within both frameworks. Focusing on the vacuum part of the spacetime (outside the sources), we derive an explicit gauge transformation between the averaged radiation gauge and the gauge in which the metric takes its standard Weyl-Lewis-Papapetrou (WLP) form, and we express the linearized WLP metric functions directly in terms of the Debye potential. We further discuss perturbations of the Kerr black hole towards general type D spacetimes, and analyze the mass and angular momentum perturbations in more detail. Finally, we illustrate the procedure on two examples beyond the type D class: a perturbation of the Schwarzschild black hole by a thin disk, and a perturbation of the Kerr black hole by a rotating particle.

We investigate the thermodynamic properties, Joule-Thomson expansion, and microscopic structure of a dyonic Reissner-Nordström-AdS black hole in nonlinear electrodynamics surrounded by quintessence. Within the extended phase-space formalism, analytical expressions for the mass, Hawking temperature, equation of state, heat capacity, Gibbs free energy, Joule-Thomson coefficient, and Ruppeiner curvature are derived. The effects of the nonlinear-electrodynamics parameter $ω$, the quintessence parameters $(c,ε)$, the dyonic charge $Q$, and the pressure $P$ are systematically investigated. The temperature profiles and equation of state exhibit a Van der Waals--like phase structure, with critical points strongly modified by the quintessence and nonlinear electrodynamivs contributions. The heat capacity and Gibbs free energy further reveal the corresponding thermodynamic stability and phase-transition behaviour. The Joule-Thomson analysis shows distinct cooling and heating regions separated by inversion points. In particular, quintessence generally lowers the critical and inversion temperatures, while increasing the charge and varying the nonlinear parameter significantly deform the inversion structure. The Ruppeiner curvature exhibits pronounced peaks and sign changes near critical regions, indicating strong thermodynamic correlations and alternating attractive and repulsive effective microscopic interactions. The limiting cases $ω\sim 1$, $c\sim 0$, and $Q\sim 0$ consistently recover the corresponding Maxwell, no-quintessence, and neutral black hole behaviours. These results demonstrate that the combined effects of nonlinear electrodynamics, dyonic charge, quintessence, and AdS pressure substantially reshape both the macroscopic phase structure and the effective microscopic interactions of the black hole.

The baryon acoustic oscillations (BAO) is a vital probe to measure cosmological distances and constrain dark energy. The ongoing and near-future space-based telescopes, including the Chinese Space Station Survey Telescope (CSST), will perform galaxy redshift surveys using slitless spectroscopy. Due to the relatively low spectral resolution and observational effects, there will be uncertainties and interloper contamination in measured redshifts. In this study, we consider the Gaussian-distributed redshift uncertainty and [O III]-H$β$ interlopers in the CSST-like slitless spectroscopic redshift (spec-z) samples, which are simulated using the approximate $N$-body code FastPM. We study the effect of spec-z errors on BAO before and after density field reconstruction. Both the redshift uncertainty and interlopers can damp the BAO signal, and the latter can further induce oscillations in the broadband shape of the power spectrum. We model the interloper effect on the propagator, i.e. the cross-correlation between the observed and initial density fields, as well as on the BAO power spectrum. We study the influence of the spec-z errors on the fitted anisotropic BAO parameters $α_{\perp}$ and $α_{\parallel}$. The resulting systematic bias on $α_{\perp}$ and $α_{\parallel}$ is mild, mostly within $0.1$ per cent and $0.2$ per cent, respectively. Even with the spec-z errors, BAO reconstruction can still significantly reduce the systematic bias and statistical errors on the $α$ parameters. Furthermore, we vary the redshift uncertainty and study the reconstruction performance.

We present the first self-consistent theoretical study of the competing $^{14}\mathrm{C}(n,γ)^{15}\mathrm{C}$ and $^{14}\mathrm{C}(p,γ)^{15}\mathrm{N}$ reactions within the same modified potential cluster model (MPCM). For the $^{14}$C$(p,γ_{0})^{15}$N reaction, total cross sections, astrophysical $S$ factors, and reaction rates are calculated using interaction potentials constrained by the available scattering and bound-state data. The astrophysical $S$-factor is estimated as $S(0)=4.5(1)$~keV$\cdot \text{b}$. Combining these results with our recent MPCM calculations for $^{14}\mathrm{C}(n,γ)^{15}\mathrm{C}$, we determine the transition temperature at which proton capture overtakes neutron capture in the production of $^{15}\mathrm{N}$. The self-consistent comparison predicts a transition temperature $T_9^{\rm c.p.}=2.5$ under Maxwell--Boltzmann statistics, significantly higher than previous estimates. The analysis is extended to Tsallis statistics, demonstrating that deviations from thermal equilibrium produce substantial shifts of the transition temperature. These results provide improved nuclear-physics input for astrophysical nucleosynthesis calculations.

In both collisional and collisionless plasmas the presence of a broad range of electromagnetic and plasma fluctuations provides anomalous electric fields that can be important for the dynamical evolution of the system as it is the case of magnetic reconnection, plasma turbulence and dynamo theory. In the context of plasma turbulence at scales larger than the ion's inertial length, the plasma satisfies the frozen-in condition, and the anomalous electric fields are produced by correlations between turbulent velocity and magnetic fields. Conversely, for collisionless plasmas and at kinetic scales, the total electric field has additional contributions that arise from kinetic phenomena, namely, charge separation, ambipolar electric fields and electron inertia, and each of these terms presents an anomalous counterpart. In this work we characterize the anomalous electric fields. We present a framework that can explain partial balance between the anomalous resistivity and anomalous transport, and we use it to study anomalous electric fields at kinetic scales. We establish how the different contributions to the anomalous electric field couple to the large-scale electric fields and show that the anomalous terms act back-reacting on magnetic-field organization and they may contribute to effective turbulent magnetic diffusion or shielding. Finally, we test two Sub-Grid-Scale models including anisotropic transport coefficients and show that although these models partially recover the spectral information, the phase coherence is not entirely recovered by these models suggesting a more complex non-linear contribution of the anomalous terms to the resolved scales.

The existence of supermassive black holes with masses $M\gtrsim10^9\,M_{\odot}$ and large dimensionless spins $χ\sim0.9-0.99$ at high redshift remains a challenge to our understanding of the early Universe. In this work, we study the adiabatic co-evolution of a Kerr black hole seed surrounded by two ultralight scalar dark matter clouds occupying different bound states, and show that this configuration allows the black hole to grow into the supermassive mass range while imprinting a characteristic final spin. The evolution proceeds through two stages. During the first stage, a spherical cloud described by the $\ell=0$ mode is completely depleted through a runaway dark-to-black accretion mechanism on a timescale of hundreds of millions of years for boson masses $μ\sim10^{-18}-10^{-17}\,\mathrm{eV}$. Since the accreted material does not carry angular momentum, the black hole spin is universally driven to $χ\simeq0$, independently of its initial spin. Throughout this stage, the second cloud, described by the $\ell=m=1$ mode, remains in the superradiant regime with negligible evolution. However, once the first stage is completed, this cloud transitions to the accreting regime, rapidly transferring both mass and angular momentum to the black hole. Starting from $χ\simeq0$, the black hole spin increases until the evolution self-consistently saturates close to the threshold $χ_{\rm sat}$, defined by the condition $Ω_H(χ_{\rm sat})=μ$, on an e-folding timescale of thousands of years, orders of magnitude shorter than the first stage. This final saturation spin is largely independent of both the initial black hole spin and the mass of the secondary cloud, providing a spin-imprinting mechanism in which the primordial spin is first erased by spherical accretion and then reset to a value determined only by the boson mass and the final black hole mass.

Magnetic fields play a fundamental role in the dynamics and radiation of jets associated with active galactic nuclei (AGNs), where they accelerate and collimate the relativistic plasma and produce the observed synchrotron emission. Faraday tomography is a powerful technique for probing the three-dimensional distribution of magnetic fields and synchrotron-emitting plasma along the line of sight, by reconstructing the Faraday Dispersion Function (FDF) from the observed polarization spectrum. To explore what aspects of AGN-jet structure can be revealed by Faraday tomography, we construct a simple three-dimensional model of an AGN jet that captures its essential features, and compute the resulting FDF for a range of model parameters. The jet is represented as a cylindrical region of uniform thermal-electron density threaded by a coherent helical magnetic field, to which a random turbulent component can optionally be added. The parameters varied include the inclination angle θ between the jet axis and the line of sight, the wavenumber of the helical field, and the amplitude of the random magnetic-field component. By systematically varying these parameters, we examine how the underlying magnetic geometry and turbulence are encoded in the observable Faraday and polarization structures. The main observational signatures identified by our model include top-bottom asymmetry across the jet, curved or double peaked FDF profiles and increasingly fragmented Faraday depth structure along longer lines of sight, providing potential diagnostics of large scale helical magnetic fields and small scale turbulent components. These signatures providing a physical framework for interpreting future polarimetric observations of AGN jets from instruments such as LOFAR, MeerKAT, and the VLA.

We explore the phenomenological viability of non-minimally coupled Weyl connection gravity by confronting its static, spherically symmetric black hole solutions with classical weak-field Solar System tests and stellar-orbit observations near Sgr A*. In this geometric framework, non-metricity is encoded via a Weyl vector field, giving rise to two distinct families of Schwarzschild-like vacuum solutions characterized by a free parameter \(ω\) with dimensions of length. We compute the corrections to four classical observables---gravitational redshift, Mercury's perihelion advance, light deflection, and radar echo delay---and derive stringent lower bounds on \(ω\) using current observational data. For Solution I (purely radial Weyl vector), the leading metric corrections scale as \(1/ω\), yielding bounds as strong as \(ω\gtrsim 10^{30}\) m from perihelion precession. For Solution II (time-radial Weyl vector), the corrections scale as \(1/ω^2\), resulting in weaker constraints, with \(ω\gtrsim 10^{20}\) m from the same test. This marked difference arises from the distinct behavior of the linear corrections in each solution: while Solution I exhibits an unsuppressed linear term \(2r/ω\) that dominates in the Solar System regime, Solution II features a linear term suppressed by an additional factor of \(M/ω\), making the quadratic term \(-r^2/4ω^2\) dominant throughout. We then analyze stellar orbits near the Galactic Center, finding that current observations of the S2 star provide complementary constraints on both solutions. (...)

Galactic binaries are expected to be the most numerous long-lived sources in the Laser Interferometer Space Antenna (LISA) data stream and are commonly modeled as quasi-circular, nearly monochromatic systems. In this work we present \texttt{eGB-multi}, a fast time-domain eccentric Galactic-binary waveform package for LISA. The package combines a post-Newtonian-accurate quasi-Keplerian eccentric source model with a time-domain LISA response, evaluating the signal on the retarded one-way links and then constructing time-delay interferometry observables. Its modular implementation allows the user to switch between different levels of orbital description, ranging from Newtonian closed eccentric orbits and relativistically precessing eccentric orbits incorporating the 1PN-accurate periastron advance, to precessing and shrinking eccentric orbits that additionally account for gravitational radiation reaction at leading quadrupolar order. The framework also provides an automatic evolution mode, in which the appropriate level of orbital dynamics is selected based on a pre-defined mismatch tolerance. We employ this framework to quantify the regimes in which orbital eccentricity renders the commonly adopted quasi-circular approximation for Galactic binaries inadequate.

Bondi showed in 1957 that bodies of opposite active gravitational mass self-accelerate: the negative chases the positive it repels, and the pair runs off together. We evolve this "Bondi dipole" in $3+1$ numerical relativity: two complex scalars share identical Klein-Gordon dynamics; only the phantom enters Einstein's equations with a minus sign -- inertial and passive masses positive, active mass negative. Across a matrix of thirty-four evolutions -- runaway pairs, controls, and parameter scans -- a mass-matched pair released at rest accelerates as a unit. The midpoint moves $3.00\pm0.01$ by $t=200$ with the separation held to 1%, and reaches a speed of $0.056c$ by $t=400$ with the acceleration steady to 2%; the total signed momentum holds at zero to $\lesssim 1$%. The force is gravity on both of its axes: $a\propto d^{-2.03\pm0.01}$ over $d=8$ to $20$, $a\propto M^{0.97\pm0.06}$ over a factor $2.5$ in mass, and $a\,d^2/\bar{M}=1$ within 2.4% on the equal-mass ladder. Swapping the sectors inverts the acceleration to two parts in $10^5$; gauge, solver-depth and mesh variations move the drift by $\lesssim 0.01$% and box doubling by 4%; same-sign control pairs hold their centroids to $\lesssim 8\times10^{-4}$ even while merging. The runaway carries no detectable gravitational radiation: the signed dipole cannot radiate, the quadrupole's $\ddot{Q}$ is constant, and the measured $\ell=2$ amplitude falls as $r^{-4.8}$ -- near zone, not flux. The phantom star is, to our knowledge, the first asymptotically flat body of negative ADM mass evolved in numerical relativity; alone it survives to $t=1000$, slowly relaxing outward.

In low-regularity spacetimes and Lorentzian length spaces, achronal causal curves play the role of null (pre-)geodesics. Because of the lack of a geodesic equation, they do not come with a canonical parametrization. In this context, we discuss a notion of affine parametrization via limits of affinely parametrized timelike geodesics. However, we point out a major drawback: an example where this approximation procedures gives a non-unique result, in a way that even completeness or incompleteness of the limit null geodesic is not well-defined. This example involves discontinuous gluing of two Lorentzian metrics across a null hypersurface to give a well-behaved Lorentzian length space and as such is, just like the parametrization problem itself, manifestly Lorentzian. In view of this, we explore possibilities for a Penrose-type singularity theorem using timelike geodesics.

Recently initiated and upcoming spectroscopic surveys, such as DESI and WST, will provide more than $10^6$ high-$z$ quasar spectra, enabling dense sky coverage by the Ly$α$ forest. This is expected to establish the three-dimensional (3D) Ly$α$ forest power spectrum, $P_{\rm 3D,α}$, as a probe of the thermal and ionization history of the intergalactic medium (IGM). To exploit this opportunity, we quantify the imprints of reionization on the post-reionization IGM using high-fidelity numerical models. We use the Sherwood and Sherwood-Relics cosmological hydrodynamical simulations with box sizes up to $160\,h^{-1}\,\rm cMpc$ to investigate the impact of box size, mass resolution, and extracted grid resolution on $P_{\rm 3D,α}$ over $2.4\leq z\leq4.8$. After applying a Zel'dovich control variate correction, simulation volume has a modest impact over most scales and orientations, whereas degrading the mass resolution produces differences of up to $\sim13\%$. Insufficient resolution of the grid used for the optical-depth calculation can artificially enhance small-scale power by up to $\sim35\%$. The timing of $\mathrm{H\,I}$ reionization leaves only a percent-level imprint on $P_{\rm 3D,α}$ at $z=2.4$, whereas varying the photoheating rate by a factor of two changes the large-scale power by $\sim4-8\%$. Our results demonstrate that numerical effects can be comparable to, or exceed, the relic astrophysical signatures encoded in $P_{\rm 3D,α}$, making numerical convergence essential for interpreting precise Ly$α$ forest measurements. The strongest astrophysical imprint arises from spatially inhomogeneous $\mathrm{H\,I}$ reionization, which enhances the large-scale power by up to $\sim70\%$ at $z=4.2$. This highlights the potential of post-reionization Ly$α$ forest measurements as a probe of the thermal history and spatial morphology of cosmic reionization.

Stellar variability has become the dominant limitation to achieving the radial velocity (RV) precision required for the detection and characterization of Earth-like exoplanets. While significant progress has been made in mitigating the effects of oscillations and magnetic activity, convective granulation and its interaction with stellar magnetic fields remain among the least understood sources of RV variability. To address these challenges, we organized the splinter session `Connecting Granulation and Magnetic Activity in Radial Velocities: The Next Breakthrough for High-Precision Spectroscopy' at Cool Stars 23. The session brought together researchers working on observations, numerical simulations, and data-driven techniques to discuss the current understanding of granulation-driven RV signals and identify the most promising directions for future progress. Through invited and contributed talks, followed by community discussions, participants emphasized the importance of combining physically motivated models with data-driven approaches, developing standardized benchmark datasets, and validating simulations against high-quality observations across a range of stellar types. The discussions also highlighted the need for coordinated observing strategies, improved characterization of individual spectral lines, and physically informed line-by-line analyses to disentangle convective and magnetic signals. This contribution summarises the scientific discussions and community perspectives that emerged during the session and outlines the key challenges that must be addressed to reach the sub-40 cm/s precision required for the next generation of RV planet searches.

Warm dark matter (WDM) has served as a valuable benchmark for constraining small-scale structure in the past decades. In this note, I examine what it would take for dark matter to be warm in the literal sense assumed by that benchmark, i.e., a thermal relic that decoupled while relativistic. Satisfying current constraints on the WDM mass, which approach the $\sim$10 keV scale, requires one of three possibilities. One possibility is that there were $\sim 10^4$ relativistic degrees of freedom in the thermal bath at the time of decoupling, which is far beyond what is available in the Standard Model or plausible extensions of it. An alternative is that there was a period of early matter domination whose entropy injection diluted the relic density. However, in this scenario, the perturbations would have evolved through an expansion history that was different from the radiation-dominated one assumed in deriving WDM transfer functions. The third possibility is a dark sector that was never in thermal contact with the Standard Model and was simply born colder via asymmetric reheating. All of these possibilities rely on strong coincidences, where physics that has nothing to do with WDM happens to provide the exact initial conditions assumed in a WDM cosmology. Meanwhile, WDM is often used as a proxy for models with self-consistent thermal histories that generically predict the suppression of structure formation on small scales in a way that is both quantitatively and qualitatively different from WDM. I therefore advocate for a transition to more expressive parameterizations and simulation-based methods in order to extract more useful information from the wealth of upcoming data.

Recent advances in cosmological measurements such as DESI and ACT may indicate an increase in the preferred value of the spectral tilt $n_s$ such that it disfavors rather popular models such as the Higgs or Starobinski inflation models. We argue that it actually means that the data is now sensitive enough to quantum corrections beyond simple tree-level models. Resurrecting the old theme of Coleman-Weinberg effective potential, we analyze the predictions of such models, as well as likelihood analysis, leading to a more established and interesting predictive framework.

Early observations with the BICEP Array 150 GHz (BA2-150) and 220/270 GHz (BA3-220/270) receivers revealed detector noise equivalent temperatures (NETs) higher than expected, together with substantial detector-to-detector and module-to-module scatter. Noise measurements acquired with multiplexing off and high frequency sampling demonstrate that this excess originates from elevated high-frequency detector noise that aliases into the science band during time-division multiplexing. We show that the excess high-frequency noise is correlated with anomalously large logarithmic TES transition slopes, α, resulting in elevated electrothermal loop gain and operation near the detector stability boundary. Measurements of α indicate values substantially larger than expected, consistent with the sharper superconducting transitions introduced by the inverted TES fabrication process adopted for BA2-150 and BA3-220/270 detectors. Operational mitigation strategies were investigated through both increased multiplexing rates and elevated focal-plane operating temperatures. Faster multiplexing reduces aliasing by shifting the multiplexing Nyquist frequency beyond the excess noise roll-off, while elevated bath temperatures reduce TES electrical power and loop gain, improving detector stability and reducing NET by approximately 10%. These results demonstrate the importance of balancing TES responsivity, electrothermal stability, and multiplexed readout performance in next-generation CMB polarimeters.

Main sequence stars of spectral types late F through M exhibit systematic variability in photometric light curves, particularly when they are young. Rotational modulation of starspots manifests as quasi-sinusoidal variability, which enables the measurement of rotation periods. Variability can also be stochastic, as in stellar flaring. However, since measurements of stochastic processes depend on the time of observation, they are typically noisier. Considering that different manifestations of variability have unique observational nuances, models that naturally unify these are incredibly useful for stellar characterization. Towards this goal, we have developed EncoTESS: a Time Series Foundation Model (TSFM) trained on a subset of TESS 2-min light curves. EncoTESS is specifically designed to handle the observational noise, heteroskedastic measurements, irregular sampling, and large data gaps common to TESS data. It is also ~1% of the size of a typical literature TSFM, so can be run easily on a modern laptop. EncoTESS encodes light curves into a fixed-size latent parameter space, which can be used to infer physical stellar properties and recovers light curve summary statistics well. EncoTESS outperforms rotation period and variability amplitude as age indicators for stars that have not converged onto the slow rotator sequence yet; broadly these include K and M stars less than ~100 Myr, and M stars less than ~1 Gyr. We focus on age inference as an application of EncoTESS in this work, but other downstream tasks such as stellar classification could also be explored. The architecture of EncoTESS enables its future extension to TESS light curves of all cadences, and additional surveys such as Kepler and the upcoming PLATO mission. The core EncoTESS framework and library of encodings produced for the stars used in this work are publicly available at https://github.com/philvanlane/encotess.

The recent JWST discovery of a numerous population of broad-line emitters (BLEs) at $z>4$ has reopened the question of whether sources other than faint galaxies contributed significantly to cosmic reionization. To assess this contribution, we present a systematic census of compact, blue broad-line emitters at $4\leq z\leq7$ selected from the DJA, designed to isolate sources in which the UV continuum and the broad-line emission originate from the same compact physical region. Starting from a parent sample of 4145 galaxies with NIRSpec/PRISM spectroscopy and NIRCam imaging, we apply criteria on broad Ha emission (FWHM$>2000$km s-1, detected at S/N$\geq5$, complete at $L_{Ha,broad}\geq10^{42.7}$erg s-1), blue continuum slopes ($β_{UV}\leq-1.5$, $β_{opt}<0.5$), and morphological compactness in both the rest-frame UV and optical, returning a final sample of 20 sources. This constitutes $\sim20%$ of the PRISM-selected broad-line sample, which is dominated by Little Red Dots (LRDs). We find that the compact blue BLEs display hot ionizing continua, high-ionization UV lines (CIV, NIV], HeII), and elevated Lya emitter fractions ($X_{Lya}=50-67%$) relative to star-forming galaxies and LRDs. Beyond their similar lack of X-ray emission, their optical spectra closely resemble those of LRDs. This suggests that compact blue BLEs extend this population towards lower column density envelopes, but with similar engines, rather than representing a physically distinct population. Using Sirocco to model physical configurations that reproduce the observed spectra, we find a high ionizing photon production efficiency (log$ξ_{ion}\sim25.4$Hz erg-1) and a weighted average escape fraction $fesc\sim0.2$ for the bluest sources, and $\sim0.01$ for LRDs. This implies that broad-line sources do not dominate reionization globally, but given their luminosity, their contribution can dominate up to few Mpc scales.

If the Universe possesses a compact spatial topology with characteristic scale not much larger than the diameter of the observable Universe, then the associated symmetry transformations imprint specific statistical correlations on density perturbations. We derive the covariance between perturbation modes in the full topology volume and compute the induced covariance between modes in the observable subvolume. We then use this covariance to estimate the available information about topology contained in the full set of observable linear perturbations. For two sample Euclidean examples - the cubic three-torus ($E_1$) and a three-torus with a quarter turn ($E_3$) - the three-dimensional density field carries enough information to detect topology on scales up to 25 percent larger than those accessible from the two-dimensional cosmic microwave background anisotropies alone. Future probes of the cosmic density field using deep galaxy surveys or neutral-hydrogen 21-cm intensity mapping, combined with computationally intensive searches over the large parameter space of admissible topologies, thus offer a plausible avenue to extending the range of detectable cosmic topologies.

The large-scale clustering of galaxies depends not only on their internal properties but also on their location within the cosmic web, which defines the anisotropic environments where galaxies form and evolve. Filaments serve as bridges through which galaxies and dark matter flow from low-density regions toward the highest-density nodes. In this work, we aim to map galaxy bias, which quantifies how these tracers follow the underlying dark matter density field, through the cosmic web, with a special emphasis on filaments. We look for dependencies on filament properties, such as length and density, and characterize the spatial variations of the large-scale bias along the filamentary spine. We applied the DisPerSE algorithm to identify the cosmic web in the TNG300 volume of the IllustrisTNG simulation. To measure large-scale galaxy bias, we utilized an object-by-object estimator, which provides advantages over standard estimators. We find that galaxies in node outskirts exhibit the highest large-scale bias values, reaching up to $\sim4$ times the values expected from theoretical models based on halo mass alone. Low-mass red galaxies in filaments and filament outskirts also display enhanced bias, which is strongly reduced after excluding galaxies close to nodes. Together, these results suggest that proximity to massive nodes plays a central role in shaping environmental secondary bias. Furthermore, galaxy bias decreases with filament length, from mean values of $1.4$ for short filaments to $-0.5$ for longer structures. This relation persists at fixed halo mass and galaxy color, and is not primarily driven by the average local galaxy density of the filament. Finally, short filaments exhibit an approximately uniform longitudinal bias profile, whereas long filaments show an increase in normalized bias from $\sim0.85$ near the saddle point to $\sim1.06$ close to the node.

We investigate the impact of scalar-tensor theories on the 21cm brightness temperature. We consider a specific model in Horndeski scalar-tensor theory in which dark matter and dark energy in the framework of shift-symmetric generalized Galileon theories was described and an effective unified cosmic fluid was obtained. We calculate the Hubble parameter of the model and we plot the curve of the evolution of the 21-cm temperature in a specific range of the redshift parameter $z$ and compare it to the curve corresponding to $Λ$CDM paradigm.

We develop a convolutional neural network classifier that performs spectroscopic identification of stellar youth (<40 Myr) in APOGEE spectra. This classifier is sensitive to several youth-related features, including rotational broadening, which is common in younger stars, and to several discrete lines that appear to be indicative of a very spotted photosphere. The model is successful at identifying youth across a wide range of stars, achieving its strongest performance at discriminating young M and K dwarfs, while maintaining useful discriminatory power for hotter stars as well. This work enables more robust separation of pre-main-sequence stars from more evolved sources in the field even when they have comparable Teff and log(g), providing a reliable means to substantially reduce contamination in photometrically selected YSO candidates. In addition to constructing a classifier for APOGEE spectra, we also perform classification of young stars in optical BOSS spectra.

We introduce the concept of the optimal redshift for dark energy, a statistically motivated redshift at which departures of the dark-energy equation of state (EOS) from the cosmological-constant value are tested most effectively. Within a generalized CPL parameterization, we derive an analytic expression for the optimal redshift by maximizing the separation of the EOS from $w=-1$ relative to the corresponding uncertainty. We establish the optimal redshift relationship to the phantom-crossing and pivot redshifts. As an illustration, we apply the formalism to the dataset combination of DESI DR2 BAO measurements, the DES Year-6 independent BAO measurement, and the recalibrated DES-Dovekie supernova sample. Adopting the null hypothesis $\mathcal{H}_0:w(a_{\rm opt})=-1$, we find a tension of ~$2.8σ$ with the cosmological-constant prediction at the optimal redshift, compared with ~$2.6σ$ at the pivot redshift, despite the latter having a smaller EOS uncertainty. This behavior reflects the specific design of the optimal redshift to maximize the statistical significance of a departure of the EOS from the -1 value of the $Λ$CDM model. We further provide statistical and geometrical interpretations of the optimum. In the CPL parameter space, the pivot corresponds to the projection that minimizes the variance of the equation of state, whereas the optimum maximizes its squared-distance from the cosmological-constant value over its variance. While the optimum is dataset and parameterization dependent, the underlying optimization principle and definition of the optimal redshift can be extended beyond the CPL framework. Future applications to DESI, Rubin LSST, Euclid, Roman Space Telescope, and other Stage-IV dark-energy surveys appear particularly promising.

We study static and spherically symmetric traversable wormholes in $f(T,τ)$ gravity, where the torsion scalar $T$ is coupled to the trace $τ$ of the matter energy--momentum tensor. We consider the linear model $f(T,τ)=T+βτ$ with an anisotropic fluid and the mean-pressure matter Lagrangian $\Lm=\Pmean=(p_r+2p_t)/3$. The field equations are obtained for the Morris--Thorne geometry without fixing the redshift or shape function at the outset. For a constant redshift function, the energy-condition problem takes a simple form. On the branch $β>8π$ and for $b(r)>0$, the energy density together with the null, weak, and strong energy conditions is satisfied throughout the spacetime if and only if $r b(r)$ is non-increasing. The same condition also implies asymptotic flatness, $b(r)<r$ outside the throat, and $b'(r_0)\leq -1$. The allowed geometries can therefore be written as $b(r)=r_0^2 h(r)/r$, where $h(r_0)=1$ and $h(r)$ is positive and non-increasing. For the representative family $b(r)=r_0(r_0/r)^n$, the null, weak, and strong energy conditions hold for $n\geq1$, while the dominant energy condition requires $n\geq3(β-2π)/(β-6π)$. We also separate the physical matter from the effective source and show how the trace coupling allows the physical matter to remain non-exotic although the effective source violates the null energy condition. Finally, we examine the marginal case $b(r)=r_0^2/r$ with a non-constant redshift function. A decreasing redshift function can improve the tangential null energy condition at the throat, but this improvement cannot be maintained throughout an asymptotically flat exterior. These results show that the matter--torsion coupling can support a broad class of traversable wormholes without requiring exotic physical matter.

Previous work has shown that first-order Palatini gravity may be cast as a non-conservative Herglotz Lagrangian field theory, which excludes the conformal mode from its space of dynamical variables. In the present work, we analyse the role of the $SO(1,3)$ gauge connection in the construction of this action-dependent theory. The connection may be decomposed into `shape' and `scale' pieces, where only the former is retained within the Herglotz theory, while the latter is incorporated into the action-dependent sector. It is shown that the choice of connection decomposition is a space that can be smoothly parameterised. Such an observation leads to the finding that there exists a `moduli space of frictional field theories', in which each point corresponds to a distinct splitting of the shape and scale connection dynamics. Movement within the moduli space transitions between Herglotz Lagrangians whose on-shell dynamics all reproduce first-order Palatini gravity. The distinction lies in how the reproduction of the scale dynamics of the original theory is partitioned between action dependence and dynamical torsion.

The emerging picture of radio haloes and mini-haloes in galaxy clusters is more blurred and far less distinct than once thought. While both trace cosmic ray electrons and magnetic fields in the intracluster medium (ICM), our knowledge of the mechanisms powering these sources, their properties and their connection to cluster dynamics, is incomplete. We present new deep MeerKAT L-band (856$-$1712\,MHz) observations of the sloshing cool-core cluster 2A~0335+096, performed as part of the ``MeerKAT-meets-LOFAR'' mini-halo census, as well as reprocessed LOFAR 120$-$168\,MHz LOFAR Two-metre Sky Survey data. We present highly sensitive radio continuum images, from which we perform a spectral analysis and a comparison of the thermal and non-thermal properties using reprocessed archival \textit{XMM-Newton} data. Our MeerKAT data reveal that the known mini-halo is embedded in a large-scale radio halo spanning up to 510\,kpc in this sloshing cluster, although only the inner mini-halo is detected by LOFAR. We find a typical steep spectrum of $α= -1.20 \pm 0.09$ for the mini-halo; the non-detection of the larger-scale halo by LOFAR indicates a spectral index no steeper than $-1.2$. Our results are broadly consistent with an interpretation of turbulent (re-)acceleration powered by sloshing within the cluster potential well. This work emphasises the need for using deep, high-dynamic multi-frequency radio observations in conjunction with high-quality X-ray observations to complete our understanding of the connection between radio halo properties and cluster dynamics.

We reconstruct the primordial power spectrum (PPS) across the full range of cosmological scales accessible to Cosmic Microwave Background (CMB) observations. Using the Modified Richardson--Lucy algorithm, we perform free-form reconstructions through deconvolution of Planck PR3 and PR4, Atacama Cosmology Telescope (ACT) DR6, and South Pole Telescope 3G (SPT-3G) D1 data. Across different regularization schemes, we find no evidence for significant deviations from a power-law primordial spectrum. The reconstructed spectra show a strong correlation between Planck and ACT, even at the level of localized features over their overlapping range, $0.07 < k < 0.14\,\mathrm{Mpc}^{-1}$, demonstrating consistency between the two observations. Using a complementary parametric Bayesian reconstruction, we find that the previously discussed preference for a blueward tilt in the ACT data at small scales is preferred only at the $1σ$ level.

Quantum simulators based on Bose-Einstein condensates (BECs) provide a powerful platform to study relativistic quantum-field phenomena in the lab and in particular emulate relativistic analog particle detectors. However, such detectors typically couple to excitations beyond the acoustic phonon regime, so that the nonlinear Bogoliubov dispersion introduces corrections that cause the effective relativistic description to break down. We propose to overcome this limitation using an intrinsically momentum-selective Unruh-DeWitt detector model in which a tweezer-trapped atom couples to the condensate through a multimode optical cavity. Although the atom is effectively pointlike, the cavity-mediated interaction kernel causes it to sample density fluctuations over a controllable finite region of the BEC equivalent to a finite-size detector. The effective detector size is tunable via the cavity mode structure determining a controllable momentum cutoff that filters out nonphononic excitations. We analytically calculate the response of such a detector for uniform linear acceleration and uniform circular motion and show that this momentum selectivity can enhance the signal at experimentally accessible accelerations. Compatible with existing cavity-QED platforms our scheme thus provides a versatile setting for analog quantum-field measurements or relativistic quantum-information protocols.

A cosmological network of axion strings may exist in the Universe today. If axion-like particles couple to electromagnetism, such a network induces spatially varying birefringence in the polarization of the cosmic microwave background (CMB), which can be probed by current and next-generation CMB experiments. We calibrate a loop-crossing model against a large-scale adaptive-mesh-refinement (AMR) simulation of axion-string network dynamics in the early Universe and use the calibrated model to predict CMB birefringence from recombination to today. We find that the non-detection of anisotropic birefringence in CMB observations places a strong upper bound on the electromagnetic anomaly coefficient $\mathcal{A}$ that enters the axion-photon coupling $g_{aγγ} = - \mathcal{A} α_\mathrm{em} / πf_a$. A joint analysis of available anisotropic birefringence measurements constrains $|\mathcal{A}| < 0.24$ at 95% C.L., which is independent of the Peccei-Quinn scale $f_a$, assuming that the axions are hyperlight so that the network survives until today. This limit strongly restricts the high-energy embedding of hyperlight axions, excluding the minimal Grand Unified Theory prediction for the electromagnetic anomaly coefficient at high significance. In addition, we discuss the implications of an axion-string origin for the recently reported evidence of isotropic birefringence.

We extend the earlier analysis of cosmic censorship and the Weak Gravity Conjecture (WGC) in Einstein--Maxwell-dilaton (EMD) theory by performing a complete higher-order Sorce--Wald gedankenexperiment on the static charged Gibbons--Maeda--Garfinkle--Horowitz--Strominger (GMG) black hole. In the probe limit, an extremal GMG black hole can apparently be overcharged by a test particle, seemingly violating the weak cosmic censorship conjecture (WCCC), which is nothing but first-order variational identity in the Iyer--Wald formalism. This window can be closed by either incorporating backreaction via the Hoop Conjecture or further including second-order variational identities. In this paper, we explicitly compute third-order variational identities. The extremality parameter eta remains non-negative to the third perturbative orders confirming WCCC protection, regardless of the apparent failure of all-order generalization, model-independent proof of Lu, et al.. We further incorporate the Swampland Distance Conjecture (SDC): the resulting infinite tower of light charged states discharges any would-be overcharged configuration on a timescale tau much less than 2M, exponentially faster than a horizon can be destroyed, reinforcing WCCC in dynamical regimes and strongly relaxing the original WGC charge-to-mass bound.

In this work we introduce exact cosmological solutions in General Relativity in (3+1) dimensions which take into account both the usual fluid component as well as a superfluid component describing superfluid Pionic vortices. Such superfluid vortices are source of inhomogeneities which can be treated exactly, no perturbation theory is required. These analytic solutions are relevant to obtain a non-perturbative description of inhomogeneities in cosmology. We discuss the peculiar effects of these hadronic inhomogeneities on the propagation of electromagnetic fields on such cosmological backgrounds analyzing possible observable effects. Finally, the properties that the effective superfluid satisfy for such solutions to exist are discussed, together with possible extensions of this family of solutions.

Symmetric teleparallel gravity provides an alternative description of gravitation in which non-metricity replaces curvature and torsion. Its extension through $f(Q)$ gravity offers a different geometric description of the late-time expansion of the Universe and its accelerated phase. In this work, we investigate two $f(Q)$ models, a normalized power-law model and a square-root exponential model, and test their ability to describe the late-time expansion history. We constrain the model parameters through Markov chain Monte Carlo analyses using Cosmic Chronometer measurements, DESI DR2 baryon acoustic oscillations, strong-lensing time-delay observations, and three Type Ia supernova compilations, Pantheon$^+$, Union 3.0, and DES Y5. We compare both models with the flat $Λ$CDM model using the minimum $χ^2$, Akaike information criterion, and Bayesian information criterion. The square-root exponential model provides a better statistical fit than $Λ$CDM for the combinations of Cosmic Chronometer, DESI DR2, and strong-lensing time-delay data with Pantheon$^+$ and Union 3.0, with improvements in both the goodness of fit and information criteria. The normalized power-law model remains statistically competitive with $Λ$CDM for the supernova-inclusive combinations, although the information criteria do not favor its additional parameter. We also determine the transition redshift from cosmic deceleration to acceleration for both models, obtaining consistent values across the different dataset combinations. The transition redshifts agree with observational estimates of the cosmic acceleration epoch. Overall, our results support $f(Q)$ gravity as a viable alternative to $Λ$CDM for explaining the late-time accelerated expansion of the Universe without requiring a cosmological constant.

Recent advancements in direct detection (DD) experiments stimulate the investigation of the interactions of dark matter (DM) with nucleons and nuclei. In the framework of Composite Higgs Models (CHMs) the lightest Dirac composite particle (LDCP) can be stable composing a significant fraction ξof the observed DM relic abundance. We consider the elastic scattering of the LDCP on nucleons as well as on xenon (Xe) and argon (Ar) nuclei within the CHMs in which the LDCP magnetic moment, its mass and its coupling to the Higgs doublet are suppressed by an approximate U(1) symmetry. The LDCP with non-zero magnetic dipole moment can result in a substantial enhancement of the differential event rate in the DD experiments at low recoil energies of nuclei. Assuming ξ\ge 0.1, we identify the region of the parameter space where such enhancement may be potentially observable. In addition we specify some observables that can be useful in discriminating between the DM fermions with non-zero magnetic moment and other types of dark matter particles which don't have similar electromagnetic properties.

We characterize all maximally charged black hole spacetimes in $3+1$ dimensions. More precisely, given any initial data set $(M,g,k)$ saturating the mass-charge inequality that satisfies the charged dominant energy condition, we show that $(M,g,k)$ must arise from an isometric embedding into a Majumdar-Papapetrou spacetime.

Few planetary systems form in isolation. Rather, interactions in their birth environments may sculpt their architectures or drive them into unstable configurations. Here we study the role of environmental perturbations in triggering dynamical instabilities leading to hot Jupiter formation in dense clusters, focusing on whether the elevated occurrence rate of hot Jupiters in the open cluster M67 can be explained through flyby-induced high-eccentricity migration. We develop a hybrid method for modelling the secular and tidal evolution of planetary systems under external perturbations by passing stars using a combination of analytic and numerical approaches. We evolve 10,000 realizations each of systems with either a single planet, two planets, or a planet and a stellar companion for $\sim$4 Gyr in an M67-like cluster, comparing outcome statistics against a control sample without flybys. In single- and two-planet systems, the rate of flyby-induced hot Jupiter formation is negligible. However, in systems with an initially isotropically oriented stellar companion, the cluster environment boosts hot Jupiter formation by a factor of $\sim$2, accompanied by a factor of $\sim$3 increase in the rate of planets undergoing tidal disruption. Based on observationally motivated estimates of primordial populations, our hot Jupiter yields are consistent with both field and M67 occurrence rates, provided the primordial wide binary fraction among solar-type stars is close to $\sim$50 per cent, binary-driven high-eccentricity migration is the dominant formation pathway, and a substantial fraction of tidally disrupted systems survive as hot Jupiters through partial mass loss.

We establish a uniqueness theorem in four-dimensional spacetime, demonstrating that a nonlinear spinor system naturally constrains the background geometry to the Kastor-Traschen spacetime, which describes dynamical, cosmic multi-black-holes. By analyzing a nonlinear field equation where the super-covariant derivative of a Dirac spinor is sourced by intrinsic fermion currents, we evaluate its integrability under two physical requirements: a phase-locking condition enforcing a real hermitian product $\overlineψ_R ψ_L$, and a pure-electric source condition. We show that the field equation uniquely compels the spinor to be an eigenvector of the super-covariant Dirac operator, with the eigenvalue dynamically identified as the Hubble constant of the de Sitter background. These results provide a novel geometric mechanism wherein both the cosmological constant and dynamical black hole configurations spontaneously emerge from fermion condensates.

Anti-ultralocality refers to the growth of spatial gradient terms relative to velocity terms in the coupled Einstein--scalar field equations. It is a characteristic feature of decelerated expansion before the onset of inflation. Previous numerical relativity studies have shown that anti-ultralocality prevents the onset of inflation in models with power-law inflaton potentials. In this paper, we show that models with plateau-shaped inflaton potentials, which are considered to be the simplest way to generate a tensor-to-scalar ratio below current observational upper limits, are especially vulnerable to anti-ultralocality effects. The reasons are the flatness of the plateau and the energy density gap of $\sim 10$ orders of magnitude between the Planck density and the plateau potential energy. To study the problem, we develop a protocol for assessing the viability of inflationary models in general, and we apply it to a plateau potential using a previously validated numerical relativity code. We find that, starting from generic initial conditions, the growth of gradient terms in the Einstein equations relative to non-gradient terms either prevents inflation from lasting for enough $e$-folds or triggers a phase of quantum runaway. We show that the fine-tuning of initial conditions necessary to avoid these issues becomes more severe as the energy scale of inflation is made smaller, disfavoring common approaches for reducing the tensor-to-scalar ratio.

Massive hot stars are among the major ionizing sources in the Universe. The ionizing flux of a star strongly depends on its temperature, which is best measured spectroscopically using lines of different ionization stages of the same element. For the quantitative spectral analysis of O-type stars, helium lines are commonly employed, and for B-type stars, lines of silicon and magnesium are used. At low metallicity, many of the diagnostic metal lines disappear. We conduct a systematic theoretical analysis of stellar atmosphere models of OB stars to study the effect of metallicity on key temperature diagnostic lines and how this impacts spectral classification. We computed large grids of state-of-the-art non-LTE stellar atmosphere models at four different metallicities, ranging from solar to 1/31 solar. We presented equivalent widths of selected diagnostic lines in the temperature-gravity plane and investigated their dependence on metallicity. Using Galactic stellar templates of O- and B-type stars, we established equivalent width ratios for spectral classification and applied them to our models. For the hottest stars in our model grids, He I lines become detectable at lower metallicity due to the reduced back-warming in the atmosphere, opening new possibilities for their classification. Metal lines used for the analysis of B-type stars are mostly absent in the spectra of stars with metallicity below that of the Small Magellanic Cloud, except Si lines, challenging current classification schemes. Our stellar atmosphere grids reveal that the decrease in metallicity causes many metal lines to disappear, leaving different sensitivities of He I lines the only temperature indicators for B-type stars in the optical. The ionizing flux of OB-type stars can be partially trapped within the stellar wind. As a result, at lower metallicity, a larger fraction of stars can emit hard ionizing radiation.

As cosmic microwave background (CMB) polarization experiments, including BICEP3, BICEP Array, and future BICEP experiments, achieve ever-deeper polarization maps in search of primordial B-modes sourced from inflation, constraining instrumental systematics below statistical uncertainties becomes progressively more challenging. Since polarimetry in the BICEP telescopes is performed by pair-differencing co-located, orthogonally polarized detectors, differential beam response leads to temperature-to-polarization ($T \rightarrow P$) leakage, introducing a potential systematic bias on the inferred tensor-to-scalar ratio $r$. To mitigate this leakage, the lowest-order beam mismatch modes are filtered out of the CMB polarization maps through deprojection; however, residual undeprojected modes remain. To quantify this residual contamination, we perform dedicated in situ far-field beam measurements of the BICEP receivers during austral-summer calibration campaigns. We quantify the systematic impact of the undeprojected residuals with a specialized set of timestream simulations based on the measured per-detector beams. These "beam measurement-informed simulations" yield an estimate of the false polarized signal sourced by the undeprojected residuals. We summarize the beam measurements relevant to the BK24 data release and present preliminary residual-leakage results for BICEP3 at 95 GHz. For BICEP3 over 2016-2024, deprojecting all six standard templates together with readout-crosstalk templates and their radially smoothed counterparts reduces the equivalent-$r$ leakage amplitude from $ρ=(4.5\pm0.7)\times10^{-3}$ to $(1.12\pm0.06)\times10^{-3}$. We further describe an ongoing program to extend the deprojection basis beyond its historical six modes, guided by a forward optical model that relates candidate leakage modes to perturbations of physical instrument parameters.

We re-examine the problem of mimicking the standard cosmological model, characterized by $j=1$ (where $j$ is the cosmographic jerk parameter), within the context of the first connection branch of $f(Q)$ gravity. While this problem has previously been addressed via reconstruction techniques$-$yielding the analytic form $f(Q)=-2Λ+αQ + β\sqrt{-Q}$ with $Q=-6H^2$ $-$here we tackle this problem from two distinct perspectives, applying either a cosmographic closure strategy or an auxiliary-variable hierarchy approach to the traditional dynamical systems formulation. Remarkably, the cosmographic closure renders the dynamical system integrable for $Λ$CDM-mimicking $f(Q)$ models. Consequently, we obtain closed-form analytical solutions for all relevant cosmological quantities at both the background and linear perturbation levels, with the perturbative solutions elegantly expressed in terms of generalized Heun functions. Furthermore, we analyze the structural stability of the $Λ$CDM-mimicking phase space against small kinematic deviations from $j=1$, demonstrating the robustness of these solutions. Our approach provides a systematic route to studying $Λ$CDM-mimicking dynamics without requiring a closed-form analytic reconstruction of the underlying action. Finally, we utilize this framework to establish a direct comparison between $Λ$CDM-mimicking dynamics in the first connection branch of $f(Q)$ gravity and those in $f(R)$ gravity.

The $n$-point single-minus graviton amplitudes are nonzero in a special kinematic region that preserves some supersymmetry when embedded in supergravity. We compute these amplitudes within the BFSS matrix theory by relating them to an index problem considered by A. Sen. More precisely, the single-minus amplitudes are obtained from an index that counts BPS states in the Coulomb branch of four-dimensional ${\cal N}=4$ super Yang-Mills. We explicitly check this relation in a specific kinematic region by computing the index using BPS wall-crossing formulae and matching it to the gravity result. We further employ BPS wall-crossing to show that the index satisfies a tower of $w_{1+\infty}$ soft theorems in any kinematic regime, again matching gravity. Finally, we show that for certain kinematics, the single-minus amplitudes can be realized by scattering a graviton off a plane wave.