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

The origin of the 511 keV line signal in the Milky Way remains unresolved despite decades of observations. The measured flux of $\sim 3 \times 10^{-3}\,\mathrm{ph\,cm^{-2}\,s^{-1}}$ suggests a steady-state positron injection rate of $\sim 10^{43}$-$10^{44}\,\mathrm{e^+\,s^{-1}}$. One proposed contributor to this signal is stellar flaring activity since high energy Solar flares are known to produce positrons and associated annihilation radiation. We estimate the quasi-persistent 511 keV luminosity expected from flaring stellar populations and estimate the Galactic contribution. We constrain the 511 keV fluxes for different sources in the Galaxy, and in particular globular clusters. Using Solar flare observations as a calibration baseline, we construct a hierarchical Bayesian model to link flare energy to 511 keV luminosity. We further use flare frequency-energy distributions to estimate the time-averaged positron output of stellar populations. The resulting predictions are compared with INTEGRAL/SPI observations using spatial and population-based models. We find that stellar flares fall short by several orders of magnitude in explaining the Galactic positron annihilation rate. Reproducing $\sim10\%$ of the observed luminosity in the Galactic bulge would require unphysically large maximum flare energies per star reaching up to $E_\mathrm{max} \gtrsim 10^{37-39}\,\mathrm{erg}$. Spatial modeling further shows that stellar flare scenarios cannot reproduce the observed 511 keV morphology. Stellar flares cannot be the dominant source of Galactic positrons. Although previous studies have shown that the measured 511 keV morphology is broadly consistent with old stellar populations, our results exclude normal stellar flaring activity as the underlying source for this emission.

Next generation Very Long Baseline Interferometers (VLBI) will provide dense sampling of the Fourier space together with high signal to noise ratios allowing to reliably observe and image faint jet structure in M 87 at mm-wavelength. The proposed next generation Event Horizon Telescope (ngEHT) and next generation Very Large Array (ngVLA) offers the unique capability to simultaneously resolve and image the accretion flow around the supermassive black hole in M 87 together with the jet launching and acceleration zone. In order to explore these capabilities and to provide theoretical expectations we perform general relativistic magnetohydrodynamic simulations of accretion on to black holes and jet launching. M 87 has been the target for multiple observations across the entire electromagnetic spectrum. Among these VLBI observations provide unique capability to resolve the jet structure down to several gravitational radii. In this work we provide possible observable signatures which will allow us to distinguish between different electron heating models and particle distributions. We use general relativistic magnetohydrodynamics and simulate the accretion of the magnetised plasma onto Kerr-black holes in 3D. The multi-frequency radiative signatures of these simulations are computed taking different electron heating and distribution functions into account. The results of our simulations show that with a dynamical range of $1\times 10^4$ and a frequency range from 86 GHz to 345 GHz observations with future VLBI arrays have the potential to tell turbulent and magnetic reconnection electron heating and the electron distribution function apart.

We investigate how core-collapse supernova (CCSN) rates trace the star-formation rate densities (SFRDs) over the redshift range $0 \le z \le 5$. For this we use new high-redshift results from the James Webb Space Telescope Advanced Deep Extragalactic Survey (JADES) Transient Survey (JTS, see the companion paper by DeCoursey et al. 2026), together with published CCSN rates. Using the observed CCSN rates to constrain the CCSN production efficiency relating SFRDs to CCSN rates, we examine how the inferred connection between star formation rates and CCSN production efficiency depends on the stellar initial mass function (IMF) and the adopted CCSN progenitor mass range. We find that the observed CCSN rates are consistent with dust extinction-corrected UV+IR based SFRDs for plausible CCSN progenitor masses. Using the observed CCSN rates to directly reconstruct the cosmic star-formation history, we recover a peak at z $\sim2$, in agreement with galaxy luminosity-based determinations. Allowing the IMF to evolve with redshift has only a modest impact when SFRD estimates are treated consistently, indicating that CCSN rates are not as sensitive to the change of IMF as might be assumed. Adopting higher SFRDs that include a dust-obscured population of faint millimeter sources implies a substantial and increasing fraction of missing, dust-obscured CCSNe at higher redshifts. Although the inferred fraction of CCSNe missed by the surveys depends on the adopted CCSN production efficiency, we find an increasing fraction of supernovae missed due to obscuration, rising from modest values at low redshift to a peak at z $\sim2$, and remaining substantial toward z $\sim5$.

Cold Dark Matter predicts a population of low-mass halos which are sensitive to its fundamental nature and the primordial power spectrum, yet remain undetected. Although elusive, their discovery may be possible thanks to wave-optics lensing of gravitational waves (GWs) by the superposition of many halos along the line of sight. We study the statistical properties of stochastic diffractive lensing, which imprints correlated fluctuations on the amplitude and phase of the original waveform. The stochastic distortions can be described by an orthogonal basis that captures the dominant ''tones'' associated with the dark matter properties, or dark timbre, which is not degenerate with binary source parameters. LISA is most sensitive to halos of $O(10\text{--}10^4\,M_\odot)$, and because the imprint recurs in every source, stacking $\sim(50,500)$ loud binaries could confirm them at the $(2,5)σ$ level ($\sim0.2$ to $\gtrsim4σ$ for realistic merger rates and different concentration estimations). The per-event signal is only $O(10^{-3})$ in cold dark matter, demanding major advances in waveform accuracy and data analysis. Even short of that reach, stochastic diffraction places stringent bounds on models that enhance small-scale structure, such as axion miniclusters and primordial black holes.

We study the distribution of cosmic rays (CRs) in the vicinity of their accelerators, assuming that the transport of CRs in the interstellar medium surrounding the accelerators is described by the superdiffusion, beyond the normal diffusion. We find that the superdiffusivity, which is characterized by the superdiffusion parameter $α$, impacts significantly the distribution of CRs. For impulsive injection, the CR distribution behaves a constant radial profile, except with a power-law tail for $α< 2$ or a Gaussian tail for $α=2$ at large distance, $r$, from the accelerators. For stationary injection, the radial profile of CR protons tends to being proportional to $r^{α- 3}$, while that of CR electrons can deviate from the $r^{α- 3}$ profile, due to their severe energy losses. We also compute the $γ$-ray emission, produced by the interactions of CRs with ambient gas and radiation fields, within 100 pc regions around the accelerators. We find that by investigating the $γ$-ray morphology, we can distinguish the superdiffusion from the normal diffusion with present and next-generation imaging air Cherenkov telescopes.

We study the propagation of high-frequency electromagnetic waves coupled to an axionic scalar field in curved spacetime. By applying the effective action approach, we go beyond the standard geometric optics limit and derive spinoptics equations for the axion--Maxwell theory that are valid for arbitrary curved spacetimes and arbitrary axion profiles. These equations allow us to analyze helicity-dependent corrections to photon trajectories arising from the photon's interactions with both the axion field and spacetime curvature. We present representative examples to demonstrate how light trajectories deviate from null geodesics due to spinoptics effects in the presence of an axion field.

N63A is one of the brightest supernova remnants (SNRs) in the Large Magellanic Cloud and provides an excellent laboratory for studying SNR evolution in an H II region. We present a detailed X-ray analysis of N63A using Chandra and XMM-Newton observations, combining both imaging and high-resolution spectroscopy. Our global spectral fitting reveals three distinct thermal plasma components with temperatures $kT$$\sim$0.3, 0.7, and 1.5 keV. Spatially-resolved spectroscopy with Chandra shows substantial variations in physical parameters across the remnant, with the eastern lobes of the optical nebula exhibiting significantly higher ionization timescales than the western region. These multi-temperature plasma components indicate that N63A is evolving in a cloudy interstellar medium: the low-temperature component arises from evaporated dense cloudlets that are engulfed by the blast wave, the intermediate-temperature component represents shocked inter-cloud medium, and the highest-temperature component likely results from the shocks reflected by the clouds. We estimate a Sedov age of approximately 4.5 kyr and an explosion energy of about $4\times10^{51}$ erg. Comparison of the observed abundance ratios and masses of metal elements with the nucleosynthesis models constrains the progenitor mass to approximately $20 M_{\odot}$, supporting a single-star core-collapse origin for N63A.

We revisit the cross-correlation between the unresolved $γ$-ray background and galaxy shapes to constrain the annihilation cross section of particle dark matter. Our analysis uses $γ$-ray photons from 14 years of observations with the Fermi Large Area Telescope (LAT), together with galaxy shape catalogs from the Dark Energy Survey Year 3 (DES Y3) and the Dark Energy Camera All Data Everywhere (DECADE) project, enabling us to probe cosmological large-scale signals over a common sky area of $\sim 12{,}000\,\mathrm{deg}^2$ shared by the $γ$-ray and galaxy data sets. In order to better access signals from large-scale structure, we employ a Fourier-space estimator for the cross-correlation in contrast to the previous DES Y3 analysis. We find that our measurements are consistent with a null detection in a model-independent $χ^2$ test, while template-based analyses yield signals at the $\sim 3σ$ level. Our null results exclude an enhanced annihilation cross section for wino-like dark matter with a mass of $2-3$ TeV under a modest substructure boost factor of $\sim 30$ in Milky Way-sized halos. For larger boost factors of $\sim 100$, the constraints become significantly stronger and exclude the canonical thermal annihilation cross section $\langle σv \rangle = 3 \times 10^{-26}\,\mathrm{cm}^3/\mathrm{s}$ for a $7-40$ GeV dark matter particle annihilating into $b\bar{b}$ or $τ^{+}τ^{-}$. The template-based analysis favors a power-law $γ$-ray energy dependence of the cross-correlation, but also indicates deviations from that expected based on the mean intensity of the unresolved $γ$-ray background around 100 GeV. We further consider decaying dark matter scenarios and derive $2σ$ lower limits on the particle lifetime of $\sim 10^{26}-10^{27}\,\mathrm{s}$, depending on the decay channel.

The search for the axion like particles (APLs)-one of the deepest puzzles in modern cosmology-may hold the key to understanding dark matter and dark energy. In this work, we introduce a setup taking advantage of the quantum metrology techniques to constrain the hypothetical mass and coupling constants of APLs by employing exotic pseudoscalar spin-spin interactions between fermions mediated by ALPs. A key advantage of our approach is the exploitation of position-dependent spin sensor results to the high sensitivity of the probe. To simultaneously investigate the axion mass and coupling constants, we invoke a multi-probe detection strategy. Through this strategy, we circumvent the singularity of the quantum Fisher information matrix as an ultimate upper bound on the sensitivity of any probe. For the axion masses in the range of $m_a {\le} 10^{-3} \text{eV}$, this setup can exclude values of axion coupling constants $g^e_pg^n_p$ down to $10^{-7}$.

Scalar field dark matter provides us with a well motivated alternative to conventional particle dark matter, especially when ultralight fields form coherent oscillations or compact self gravitating clumps. Here we develop three complementary Solar System and local Galactic level probes of such models. These probes pertain to ADAF-like flares from scalar clump encounters with Kuiper Belt Objects, atomic clock searches for oscillatory variations of fundamental constants and astrometric microlensing by compact scalar clumps. We derive simple sensitivity estimates and null detection bounds on the scalar clump fraction, clock couplings and compact lens abundance. Our results show that Gaia-like astrometry can probe compact scalar clumps at the percent level near $M_s\gtrsim10^{-2}M_\odot$, while future astrometric and clock experiments can extend the reach to lower masses and weaker couplings.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) achieves sub-microsecond timing precision for several millisecond pulsars in its pulsar timing array (PTA) with the objective of detecting and characterizing nanohertz gravitational waves. PSR J1713+0747 is one of the most precisely timed pulsars in the array, achieving sub-microsecond timing precision. However, in April 2021, PSR J1713+0747 underwent a sudden and unusual change in pulse shape that disrupted its timing stability. As PSR J1713+0747 is a key contributor to PTA sensitivity, variations in its pulse profile significantly affect the array's sensitivity to nanohertz gravitational waves. We apply frequency-dependent Gaussian component models to decompose the pulse profile and track the evolution of individual components through the event. This component-level method maintains phase-connected timing across the shape-change event. At L-band, the recovered TOAs have a median uncertainty of ~0.47 microseconds compared to ~0.69 microseconds for standard template matching. At 820 MHz, where profile evolution is stronger, the recovered TOAs have a median uncertainty of ~1.63 microseconds compared to ~0.67 microseconds for standard template matching. The recovered TOAs achieve timing uncertainties comparable to conventional template matching while allowing data affected by profile variability to be retained in PTA gravitational-wave analyses. These results represent an initial step toward profile-domain timing methods capable of accounting for pulse-profile evolution while reducing the need for additional timing model parameters.

We investigate the scenario of maximally invisible axions, namely QCD axions whose interactions with matter arise exclusively from the irreducible coupling to gluons responsible for solving the strong-CP problem. We first analyze the production of such axions in core-collapse supernovae. In particular, we derive the corresponding SN 1987A cooling constraint and compute the emission spectra for the dominant production channels, namely nucleon-nucleon bremsstrahlung and pion conversion. We then investigate the prospects for detecting maximally invisible axions in Cherenkov detectors, with the goal of establishing a robust lower bound on the overall detectability of QCD axions.

We present HAMCOR (Hamiltonian-based AGN Multi-constraint CORonal inference framework), a geometry-agnostic method for inferring the X-ray coronal structure of accreting black holes using reverberation-lag measurements. Unlike conventional template-fitting approaches, HAMCOR reframes coronal geometry inference as the ground-state selection of a physical Hamiltonian. The corona is represented as a discrete emissivity distribution over a cylindrical grid, and its geometry emerges from five competing physical constraints: magnetic coherence, lag consistency, illumination consistency, pair-production stability, and energy budget feasibility. Minimisation is performed via projected gradient descent with Armijo backtracking line search on the probability simplex. We validate HAMCOR on three synthetic geometries (lamppost, column, ring) using the same grid as the real-data fits, recovering spatial correlations rho = 0.24, 0.50, 0.12 and fractional lag errors below 24 per cent. A hyperparameter sensitivity analysis confirms robustness over more than one order of magnitude in the coupling constants. We apply HAMCOR to five sources spanning seven orders of magnitude in black hole mass: four AGN observed with XMM-Newton (Mrk 335, 1H 0707-495, IRAS 13224-3809, MCG-6-30-15) and the stellar-mass black hole binary Cyg X-1 (M_bh = 14.8 M_sun), recovering consistent extended disc-corona geometries across the full mass range. We further present a multi-epoch analysis of Mrk 335 across five XMM-Newton observations (2006-2019), revealing that the coronal centroid remains stable at (R_c, z_c) ~ (6.3, 0.5) r_g across flux states spanning a factor of ~15 in reverberation lag amplitude, arguing against a collapsing or expanding lamppost. Schwarzschild-Shapiro delay corrections amount to ~79 per cent of the flat-spacetime lag on average; the recovered spatial morphology is robust to this correction.

Along with a whole range of alternative theories of gravity, variants of the massive theory of gravity, i.e. the theory of gravity, in which the graviton has mass, have been actively discussed in recent years. Theorists have proposed versions of massive gravity theories that address the shortcomings of early versions of such theories. Astrophysicists and experimental physicists have been discussing limitations on the graviton mass from various astronomical observations. In particular, in the first LIGO paper, where the discovery of gravitational waves from binary black holes was reported, a limit on the mass of the graviton was obtained from the analysis of the profile of the gravitational wave signal. In the paper graviton mass constraints were obtained by analyzing the trajectory of a bright star in the vicinity of the center of our Galaxy, using observations from the GRAVITY and Keck groups. Briefly other astronomical ways to limit a graviton mass were discussed.

At high densities and moderate temperatures, hadronic matter is expected to undergo a first-order phase transition into a color-superconducting (CSC) state. A proto-neutron star describes the earliest evolutionary stages during the first seconds to minutes after core-collapse supernovae and therefore has the potential to assess the appearance of CSC phases at such high densities and moderate temperatures. To address this, we incorporate proto-neutron star conditions, considering neutrino-trapped and neutrino-transparent ones, into the equation of state including color-superconducting phases in a recently developed RG-consistent NJL model. Since the total baryon number of a proto-neutron star is conserved during its later evolution, tracking stellar configurations from an initial mass of the hot proto-neutron star to the final cold neutron star along isolines of baryon number allows us to investigate whether color-superconducting phases can form at any point along this trajectory. By mapping this multidimensional transition in the hot furnace of a core-collapse supernovae cooling from a neutrino-trapped birth state to a cold, neutrino-transparent final state, we reveal four distinct core evolution scenarios-our "petit four" of proto-neutron star evolution: a delayed collapse from the CSC phase to a black hole, a persistent CSC phase, a vanishing CSC phase, and a fleeting CSC phase. For our specific parameterization of the hadronic and the CSC equation of state, we find that a stable color-superconducting phase can only be sustained in the final cold neutron star for a narrow, high-mass region.

The Globular Clusters GMRT Pulsar Search (GCGPS) project, launched in mid-May 2023, has emerged as one of the most successful pulsar surveys conducted with the upgraded Giant Metrewave Radio Telescope (uGMRT), leading to the discovery of several new millisecond pulsars (MSPs). The recently concluded Phase I of the survey resulted in the discovery of seven MSPs across four globular clusters (GCs), all of which previously had no known pulsars. These discoveries have enabled the precise determination of the dispersion measures (DMs) for the respective clusters for the first time. Consistent timing follow-up revealed that out of the seven MSPs, three are confirmed binaries, while two are isolated. This paper summarises the design and the implementation of Phase I of the GCGPS project, presents the key scientific results obtained so far, and outlines the strategy and progress of the recently initiated Phase II observations.

The abundances of heavy elements observed in metal-poor, r-process-enhanced stars provide unique information about the astrophysical conditions in which the rapid neutron-capture process (r-process) occurs. We present a Bayesian framework to infer these conditions from stellar abundance patterns. Building on a large site-independent r-process survey, we use a Markov chain Monte Carlo (MCMC) sampler to fit observed abundances with weighted superpositions of abundances from nucleosynthesis calculations, each parameterized by an initial electron fraction, entropy, and expansion timescale. Applying this framework to the nearly complete r-process template star HD222925, we find that two components are required: A heavier H-component producing elements from the second to the third peak and a lighter L-component producing elements from the first to the second peak. Increasing the number of components does not significantly improve the agreement with observations. Extending the analysis to the limited-r stars HD128279 and HD122563, and to the solar r-process residuals, we find that their abundances can be reproduced by almost the same two components, albeit with different relative weights. Moreover, the lightest neutron-capture elements ($Z \lesssim 35$) in the Sun require additional contributions absent in HD222925. Residual discrepancies are concentrated in elements sensitive to observational systematics (e.g., Ag and Cd) or nuclear-physics inputs (around the third peak). Our study highlights the power of combining stellar abundances with nucleosynthesis calculations to constrain the astrophysical conditions of the r-process, and motivates further improvements in both observational data and nuclear-physics inputs.

The distortion of the spectrum of synchrotron radiation received by a distant observer due to magnetic field inhomogeneity is discussed. It is shown that for single-particle radiation, a noticeable difference in the integrated spectrum from the case of a uniform magnetic field occurs only for observation times significantly exceeding the observer's departure from the beam pattern, and only for subrelativistic particles, when the radiation arrives from regions with significantly different magnetic fields during the observation time. For shorter observation times, field inhomogeneity leads only to small-scale distortions (smaller or comparable in frequency to the gyrofrequency) and does not manifest itself in the averaged spectrum. For the case of a particle ensemble, which is of interest for astrophysical applications, the difference in the integrated spectrum becomes negligible.

We present a search for high energy gamma rays in the energy range from 1--300~GeV from the galaxy cluster PSZ2~G181.06+48.47 using 17.9~years of Fermi-LAT data. A binned likelihood analysis employing the 16-year source catalog reveals a significant $γ$-ray excess at the cluster position. Modelling the emission as a point source yields a detection with $\mathrm{TS}=19.0 (4.3 σ)$, a photon index of $Γ=3.20 \pm 0.62$, and an integrated photon flux of $(9.75\pm2.74)\times10^{-11}\,\mathrm{ph\,cm^{-2}\,s^{-1}}$, but leaves statistically significant residual emission at the cluster center. Replacing the point-source hypothesis with a RadialGaussian spatial template significantly improves the fit, with a preferred width of $σ=0.4^{\circ}$ and an extension significance of $\mathrm{TS}_{\rm ext}=36.8 (6.0σ)$. The adopted Gaussian model yields a source detection significance of $\mathrm{TS}=55.0 (7.4σ)$, a photon index of $Γ=2.62\pm0.26$, and an integrated photon flux of $(3.22\pm0.50)\times10^{-10}\,\mathrm{ph\,cm^{-2}\,s^{-1}}$, while reducing the residual emission at the cluster position to a level consistent with zero. The cluster spectral energy distribution shows significant emission only in the lowest energy interval (1.0--3.13~GeV), while all higher-energy bins are consistent with upper limits, indicating a soft $γ$-ray spectrum with no evidence for emission above $\sim5$~GeV. These results provide strong evidence that the $γ$-ray emission associated with PSZ2~G181.06+48.47 is spatially extended on a scale of $σ\approx0.4^{\circ}$ and is more consistent with diffuse intracluster emission than with a single unresolved point source.

Black holes are the sources of the strongest gravitational fields in the present-day Universe, offering unparalleled opportunities to test Einstein's theory of General Relativity in the strong-field regime. In this talk, I will examine the prospect of sending small spacecraft on an interstellar mission to the nearest black hole. While highly speculative and fraught with technical challenges, such an endeavor is not entirely beyond the realm of possibility. Although the necessary technology does not yet exist, it may become available within the next 20 to 30 years. The mission itself could span 80 to 100 years, yet the scientific return would be profound and potentially unattainable through other means.

We generate mock X-ray observations from a suite of idealized high-resolution ($\sim 4$ pc), tall-box ($\sim 2 \times 2 \times 8$ kpc$^3$) simulations of star formation driven galactic winds in an M82-like system, varying the spatial resolution as well as the strength and distribution of supernova (SN) energy injection. We compare our mock X-ray observations with deep Chandra observations of the hot plasma around M82. While the simulated total X-ray luminosity, $L_X$, increases with resolution and when SNe feedback is spatially distributed, even in the best case scenario, our simulated $L_X$ is a factor of $\sim 50-100$ lower than observed and the surface brightness profiles of X-ray emission, $S_X(z)$, fall off too quickly with distance from the galaxy. Past results were able to reproduce these observables and we discuss potential simulation differences that could explain this discrepancy. We make the first comparison of the X-ray spectrum of our simulations to observations and find that our simulated spectrum is too soft, with a deficit of hard X-ray photons at $\gtrsim 1$ keV. We discuss how physical processes missing from our simulations and prior work (e.g., thermal conduction and cosmic rays) could help resolve this discrepancy.

Observations of compact binary coalescences (CBCs) by the LIGO-Virgo-KAGRA collaboration resulted in 90 events in the first three observing runs. We search the data near the loudest of these detections for so-called "gravitational glints," by performing Bayesian model comparison. Gravitational glints are gravitational-wave tails caused by interactions between a signal and a spacetime perturber and are theorized to present as echoes of the primary signal. We do not find convincing evidence of gravitational glints in binary black hole or intermediate-mass black hole binary events with a signal-to-noise ratio (SNR) of at least 12. We also found that standard matched-filter CBC searches with template banks that do not contain gravitational glints would likely not have missed signals containing a gravitational glint. We therefore estimate upper limits on the probability of a glint as a function of its relative amplitude and use these estimates to constrain the properties of perturbers in the Universe. With increased sensitivity in the fourth observing run, which will increase the number and SNR of detectable events, we remain optimistic that glints will soon be detectable.

We develop a general Fokker-Planck framework describing the dynamical evolution of Keplerian binaries subjected to stochastic perturbations. The formalism provides an algorithmic way to obtain the Fokker-Planck drift and diffusion coefficients of any set of orbital variables given the statistics of the perturbations. We apply the method to three physically distinct regimes: adiabatic tidal perturbations, white-noise tidal perturbations, and impulsive encounters with a third body of arbitrary density profile. In each regime we provide explicit drift and diffusion coefficients for all six orbital elements, derive the associated evolution timescale, and obtain analytic steady-state distribution functions. Our results extend previous treatments by including the evolution of the binary's orientation, retaining the complete tensor structure of tidal correlators, treating non-pointlike perturbers, and resolving the exact geometry of impulsive encounters. The latter correction leads to a steady-state eccentricity distribution that is slightly sub-thermal. We also show how these equations can be applied directly in several astrophysical scenarios, including binaries perturbed by dark matter subhaloes, ultralight dark matter, and the interstellar medium. This work delivers both a complete mathematical framework and a practical toolkit for stochastic binary evolution, providing ready-to-evaluate equations to be applied directly to binary population data.

Massive astrophysical objects can source huge static configurations of a scalar field. When such an object ends up forming a black hole, for instance, via a core-collapse supernova, the scalar field loses its source abruptly; then the static configuration becomes dynamical and propagates away in a burst, a "scalar tsunami". These bursts accumulate over cosmological history, forming a relic stochastic diffuse scalar background peaked in the $1-10^3~\text{Hz}$ range. We propose this as a novel mechanism for the generation of such a background, compute its spectrum, and compare it with the sensitivity of future experiments. We show how this extends the experimental sensitivity to scalar masses $m_φ\lesssim 10^{-13}~\text{eV}$, ten orders of magnitude larger than those accessible via individual transient events previously considered.

We establish a universal connection between the classical energy conditions and directly observable black hole properties. By imposing the energy conditions locally at the photon sphere, we derive analytic and testable bounds on the shadow size, Lyapunov exponent, Lyapunov time, and photon-ring time delay that apply to all static, spherically symmetric black holes in General Relativity. Expressed equivalently in terms of either the spacetime metric functions or the local components of the stress-energy tensor, these bounds provide complementary tests of strong-field gravity that are independent of any specific black hole solution or matter model. They are directly applicable to current and forthcoming horizon-scale observations, constrain both admissible black hole spacetimes and their underlying matter sources, and establish the classical energy conditions as observational probes of strong-field gravity. Any observational violation of these bounds would signal either a violation of the classical energy conditions by the matter fields or a breakdown of General Relativity in the strong-field regime.

Kinetic instabilities, notably the whistler heat-flux instability (WHFI), are known to suppress thermal transport significantly in the moderate- to high-$β$ plasmas relevant to many astrophysical systems. This paper explores WHFI-regulated heat transport in a new regime: ultra-high-$β$ plasmas with $β_{e} \gtrsim L_{\mathrm{T}}/ρ_e$. Extrapolating previous theories of the WHFI to ultra-high-$β$ plasmas, we propose that the magnetic energy in unstable whistler fluctuations becomes comparable to that of the background magnetic field at saturation. We corroborate this hypothesis using 1D3V and 2D3V kinetic simulations using the particle-in-cell code OSIRIS. We find that, in ultra-high-$β$ plasmas, the heat flux is localised and no longer regulated primarily by resonant pitch-angle scattering of electrons; instead, thermal energy is transported predominantly by advection at the whistler phase velocity. Heat-flux suppression is observed in 1D3V and 2D3V simulations; however, we show that the saturation of the WHFI and the regulation of heat flux are sensitive to dimensionality in the ultra-high-$β$ regime. The amplitude and phase velocity of the heat-flux-regulating whistler waves scale differently with $β_e$, yielding parallel heat fluxes, normalised to the free-streaming value, of $q_{e\parallel} / q_\mathrm{fs} \approx 4.7 β_{e}^{-1}$ and $q_{e\parallel} / q_\mathrm{fs} \approx 0.3 β_{e}^{-1/2}$ in 2D3V and 1D3V simulations, respectively. We perform 2D3V simulations with background magnetic fields inclined to the temperature gradient, showing cross-field heat transport remains negligible. We develop a heuristic theory from kinetic equations that explains these phenomena. Our work extends our understanding of how the WHFI modifies thermal transport to regimes applicable to high-energy-density physics and the reionised intergalactic medium.

We study the problem of measuring the anisotropy of high-energy cosmic rays across all angular scales. The limited field of view and non-uniform exposure of ground-based cosmic-ray experiments reduce their sensitivity to real anisotropy of cosmic-ray arrival directions. A widely used signature of anisotropy -- a dipole of the flux expansion over the right ascension -- provides very limited understanding of the underlying physics of cosmic-ray origin and propagation. In this study we test two other methods: the angular power spectrum and the autocorrelation function for sensitivity to possible medium- and small-scale anisotropies of cosmic-ray flux. We find that the autocorrelation function is the most sensitive estimator for an underlying physical anisotropy in most of the models tested, while the angular power spectrum can provide additional knowledge about cosmic-ray flux properties, when the anisotropy is strong enough. As a test of our findings, we apply these methods to $10\%$ sample of the KASCADE experiment public data. Namely, we consider all-particle set of events and sets of individual mass groups classified by a convolutional neural network. We find an indication of anisotropy at $> 2.5 σ$ level at $\sim 10^\circ$ angular scale for the iron nuclei mass group at $E \gtrsim 20$ PeV with both the angular power spectrum and the autocorrelation methods.

In recent years, NICER data has been extensively used to determine neutron star radius and mass using pulse profile modeling. The pulse profile modeling is implemented with the {\tt X-PSI} package and best-fit parameters are obtained using Bayesian inference. Using synthetic data, we demonstrate the application of frequentist inference to determine the neutron star radius, where the nuisance parameters are treated using profile likelihood. We find that the profile likelihood technique can recover the true radius to $< 1σ$. Its precision is also comparable to that of Bayesian analysis while being computationally much faster. Therefore, this work serves as a proof-of-principle application of frequentist inference to determine neutron star radius using pulse profile modeling and complements the Bayesian inference technique used. We have also made our analysis codes for frequentist inference using {\tt X-PSI} publicly available.

Binary black hole mergers without electromagnetic counterparts are expected to dominate O5 gravitational-wave catalogs. Recent CHIMERA~2.0 forecasts typically fix $Ω_m$ to a CMB-informed value and use spectroscopic hosts when available, but the corresponding sensitivity of the $Ξ$ posterior has not been assessed for pure dark sirens on the public O5 mock catalog. We analyze 300 O5-sensitivity mock events without galaxy catalogs, varying $Ω_m$ over $[0.20,\,0.35]$ (including Planck $0.315$), and compare fixed-background inference with joint $(H_0,\,Ω_m,\,Ξ)$ inference. The marginalized $Ξ$ posterior is $0.9783 \pm 0.3548$ and shows no change across this interval. Only GW luminosity distances enter the analysis, so the likelihood constrains $Ξ\, D_L^{\rm EM}(H_0,\,Ω_m)$; when $Ω_m$ is changed, $H_0$ shifts to compensate and the $Ξ$ marginal remains unchanged. Joint inference gives $Ξ^{\rm joint} = 0.9550 \pm 0.3710$, with $|ρ| \lesssim 0.05$ for the $Ω_m$--$Ξ$ and $H_0$--$Ξ$ pairs, whereas $ρ_{H_0 Ω_m} \simeq -0.4$ and the $H_0$ median moves by $\simeq 4.2\,{\rm km\,s^{-1}\,Mpc^{-1}}$ over the adopted $Ω_m$ range. Galaxy-catalog analyses on the same events at fixed $Ω_m = 0.3$ reach $\sim 7.5\%$ precision on $Ξ$, compared with $\pm 36.3\%$ here; the larger uncertainty is driven mainly by missing host redshifts. Sub-percent $Ξ$ tests will therefore still require measured redshifts even if dark sirens dominate the detection rate.

We calculate high-order late-time tails of the retarded Green function of the Teukolsky equation for linear field perturbations of (subextremal) Kerr spacetime. We calculate these tails at a fixed spheroidal harmonic $\ell$ and azimuthal number $m$ up to the first three orders for the field point: at finite radius (away from the event horizon) for large Boyer-Lindquist time $t$; along the future event horizon $\mathscr{H}^+$ for large ingoing Eddington-Finkelstein coordinate $v$; and along future null infinity $\mathscr{I}^+$ for large outgoing Eddington-Finkelstein coordinate $u$. We obtain the tail powers for generic integer field spin $s$ and the tail coefficients specifically for gravitational ($s=-2$) perturbations. Our asymptotics include the known leading power-law (generic) tails, respectively,$t^{-2\ell-3}$, $e^{imΩ_H v}v^{-2\ell-3-b}$ (where $b=1$ for $s>0, m=0$ and $b=0$ otherwise, and where $Ω_H$ is the angular velocity of the event horizon) and $u^{-\ell+s-2}$, as well as their higher-order logarithmic corrections: $t^{-2\ell-5}\ln t$, $e^{imΩ_H v}v^{-2\ell-5-b}\ln v$ and $u^{-\ell+s-3}\ln u$ (as well as $u^{-\ell+s-4}\ln^2 u$). Since we obtain the high-order expansions for modes for generic $\ell$ and $m$, we can readily infer the explicit expansions of the {\it full} retarded Green function for $s=-2$ (and its decay powers for generic integer $s$). We obtain the late-time asymptotics from small-frequency expansions of the Fourier modes of the retarded Green function in the frequency domain. Accordingly, we also provide small-frequency expansions of various quantities of interest in the scattering theory. We also attach two notebooks which provide expansions for specific values of $s$: one notebook provides them to the first three leading orders for generic $\ell$ and the other one to arbitrary order for specific values of $\ell$.

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-$β$, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (ICWs). Prior work studying this stream (Bowen et al., 2025) showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process, and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions (VDFs), from the ion electrostatic analyzer (SPANi) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via non-coherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier, and heating rates for RH via $\parallel$-ICWs do not peak in the same region of velocity-space as the empirical measurements, nor reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas, and shows evidence of a Fokker-Planck like diffusive process in the near-Sun solar wind.

We develop a kinematic framework that relates the monopole, dipole, and quadrupole of the luminosity distance to an ellipsoidal peculiar velocity field describing the dynamics of the Laniakea supercluster. By properly accounting for the transformations between the CMB and Laniakea reference frames and selecting Type Ia supernovae within the volume associated with the superstructure, we show that luminosity distance multipoles encode the expansion and shear of the local velocity field. This allows the ellipsoidal kinematics of Laniakea to be inferred directly from supernova observations. Our results provide a physically motivated interpretation of luminosity distance anisotropies within the volume dynamically dominated by Laniakea as signatures of its large-scale kinematics. The framework developed here suggests that future supernova observations can enable the independent reconstruction of large-scale structures in the local Universe

The circumgalactic medium (CGM) -- the multiphase gas reservoirs surrounding galaxies -- remains the least understood component of the baryon cycle governing galaxy growth, despite its central role in the Astro2020 Decadal Survey's priorities. Existing constraints come almost exclusively from pencil-beam absorption spectroscopy, leaving the spatial structure, kinematics, and phase interactions of CGM gas fundamentally unmapped. We present Ardua, a mission concept for NASA's ASTRA Initiative that combines wide-field far-ultraviolet spectroscopy with a Line Emission Mapper (LEM)-derived X-ray microcalorimeter instrument to obtain the first comprehensive emission maps spanning the full CGM temperature range, including cool neutral gas, ionized warm-hot phase gas, and the volume-filling hot corona. By observing more than 50 nearby galaxies comprehensively in the UV and X-ray, Ardua will test competing galaxy formation models, resolve multiphase gas flows and feedback-driven outflows, and extend baryon-cycle science to the intergalactic medium and the environments of exoplanet-hosting stars. Beyond its core CGM/IGM program, Ardua's wide-field, high-sensitivity instruments are designed to serve as a flexible community resource, supporting guest-investigator science across astrophysics. No planned or approved mission is designed to deliver this combined UV/X-ray survey capability.

Spacetime torsion provides a natural extension of general relativity and may lead to black hole solutions that differ significantly from their Einsteinian counterparts. We investigate a class of Reissner-Nordström-like black holes in Poincaré gauge gravity, where the effective charge is generated entirely by spacetime torsion instead of an electromagnetic field. Within the physically relevant torsion sector, the spacetime exhibits a single-horizon structure, free from the inner horizons and extremal states characteristic of charged black holes. Using the sixth-order Wentzel-Kramers-Brillouin (WKB) approximation, Leaver's continued-fraction method, and the eikonal correspondence between quasinormal modes and unstable null geodesics, we study scalar perturbations and spin-2 test fields on the torsion-modified background. We find that increasing torsion decreases both the oscillation frequencies and damping rates, leading to longer-lived ringdown signals. We further compare the model predictions with Event Horizon Telescope observations of Sgr A* and M87*, along with representative LIGO-Virgo-KAGRA (LVK) ringdown scales, to derive constraints on the torsion parameter via a profile-$χ^{2}$ analysis supplemented by Monte Carlo sampling. Although the resulting bounds remain consistent with the Schwarzschild limit within current observational uncertainties, our results show that spacetime torsion leaves correlated imprints on both black hole shadows and ringdown observables.

Although currently poorly constrained, cosmological ultra-large scales are expected to provide formidable tests not only of General Relativity, but also of the content of $Λ$CDM and the Early Universe. However, in this regime, cosmic variance plays a major role in limiting sensitivity, and controlling systematic errors becomes a crucial aspect in preserving the limited information content of current and future observations. Multi-tracer analyses represent a useful technique that simultaneously allows to limit the impact of cosmic variance and mitigate the presence of systematics. In this sense, gravitational waves might represent the perfect alternative tracer of the large-scale structure, since their detection suffers from a set of uncertainties completely different from that of traditional large-scale structure surveys. In this work, we provide a concrete example of how gravitational wave clustering mitigates the presence of systematics, and facilitate the discovery of New Physics signatures. Specifically, we focus on systematics that degrade the constraining power on local Primordial non-Gaussianity for future galaxy surveys; and show how catalogs of gravitational wave events detected by third-generation observatories reduce the impact of systematics while being able to maintain flexibility in the statistical analysis. Additionally, we release a new version of the Multi_CLASS code, which now provides an enhanced level of customization of the different tracers and is compatible with the latest releases of CLASS.

Close white dwarf (WD) + main-sequence (MS) binaries are end products of mass transfer (MT) that occurred when the WD progenitor was a giant, making their population demographics a powerful probe of binary evolution. Several recent works have constructed samples of WD+MS binaries with well-understood selection functions using data from wide-field surveys. These include (a) au-scale astrometric binaries from Gaia that can be shown to contain a WD on dynamical grounds, (b) au-scale astrometric binaries in which a hot WD is detected through a GALEX UV excess, and (c) close binaries discovered through eclipses. Together, these samples probe outcomes of both stable MT and common-envelope evolution, and interactions on both the red giant branch (RGB) and asymptotic giant branch (AGB). We forward-model the three observed samples simultaneously. This approach produces robust constraints on uncertain binary evolution parameters because binaries removed from one population are predicted to appear in another. Our modeling includes a realistic initial binary population and treatments of the selection effects affecting all samples. We find that MT from AGB donors requires a critical accretor-to-donor mass ratio of $\sim0.4$ and is more stable than MT from RGB donors, for which we find a critical ratio of $\gtrsim0.65$. A common-envelope efficiency of $αλ\sim0.3$ matches the relative numbers of close and wide systems and the period distribution of close systems. Most stable-MT products in the sample, including those with RGB donors, retain nonzero eccentricities ($\simeq0.1$). The model does not fully reproduce the mass distribution of main-sequence stars in post-common-envelope binaries, which exhibits a cliff well below the fully convective limit, pointing to missing physics that should be explored in future work.

Recent studies have suggested that galaxy mergers may help drive the escape of ionizing photons from galaxies. However, the ionization properties of merging galaxies have not yet been systematically studied across the full merger sequence in a large, well-defined sample. Here, we investigate the impact of mergers on the line ratio $\rm{O}_{32}$ tracing the ionization state and study its evolution as a function of merger progression. We identify 7641 galaxy mergers with robust ($\rm{SNR}>3$) emission line fluxes from the Sloan Digital Sky Survey Data Release 17, using spectroscopic pair analysis and visual classifications from the participatory science experiment "Cosmic Disco: Characterizing Galaxy Collisions". We find a continuous increase in $\rm{O}_{32}$ as a function of merger progression traced by pair separation and visual merger stages, with median and 90th percentile peaking post-coalescence. $\rm{O}_{32}$ is significantly enhanced in mergers past the first pericenter passage as compared to a control sample of isolated galaxies. We investigate the properties of galaxy mergers with $\rm{O}_{32}>3$ and find that these are mainly low-mass ($\widetilde{\rm{M}_{*}}=6\times 10^{8}\,\rm{M}{_\odot}$) mergers close to coalescence. In this population, large $\rm{O}_{32}$ values are due to mergers driving a simultaneous increase in specific star formation rate and decrease in metallicity. Thus, galaxy mergers enhance the ionization state of galaxies and likely facilitate the escape of ionizing radiation from galaxies, even at low stellar masses. Finally, we infer that Green Pea galaxies, a population characterized by high $\rm{O}_{32}$ , are low-mass galaxy mergers at coalescence, a scenario that reconciles findings from environmental and interferometric 21cm studies.

We present updated bounds on ultra-light axions (ULAs) as a subcomponent of dark matter, derived from the full-shape analysis of the DESI Data Release 1 galaxy power spectra combined with Cosmic Microwave Background (CMB) data from ACT and Planck. We focus on the mass window $10^{-32}\,\mathrm{eV}\leq m_a \leq 10^{-24}\,\mathrm{eV}$, employing state-of-the-art analysis methods rooted in the Effective Field Theory of Large Scale Structure. For the smallest masses, our joint analysis with DESI improves over CMB-only constraints by more than a factor of 2, establishing the most stringent limits to date. For instance, for $m_a \sim 10^{-29}\,\mathrm{eV}$ ULAs are constrained to be a fraction as small as $0.3\%$ of the total matter energy density. The DESI Luminous Red Galaxy sample shows a mild preference for an ULA subcomponent with $m_a \approx 10^{-26}\,\mathrm{eV}$, mirroring previous hints from BOSS, but this preference vanishes upon combination with CMB data. Probing the largest masses, $m_a\gtrsim10^{-25}\,\mathrm{eV}$, in future studies will benefit from extending the data analysis to smaller scales, both for galaxy clustering and CMB lensing, but will also require concurrent improvements in the theoretical modeling.

We investigate the long-wavelength evolution of linear perturbations in a homogeneous and anisotropic background with a scalar field coupled to a vector field. Using the spatial gradient expansion in the uniform-$\mathcal{N}$ gauge in which the number of $e$-folds is unperturbed, we derive the complete set of superhorizon solutions and establish their correspondence with infinitesimal variations of the homogeneous anisotropic background. This extends the separate-universe picture, previously known for isotropic FLRW cosmology, to anisotropic spacetimes despite the mixing of scalar, vector, and tensor perturbations induced by broken rotational symmetry. We show that the long-wavelength equations form a self-consistent system and identify a conserved quantity that generalizes the conserved Wronskian of isotropic cosmology. Unlike the isotropic case, the superhorizon modes sourcing the curvature perturbation are governed by three independent conserved channels associated with the scalar field, the background shear, and the gauge-field tilt, together with an additional dynamical shear contribution originating from the anisotropic geometry. This reveals that the evolution of curvature perturbations around anisotropic background is intrinsically richer than in isotropic multi-field models. Our formulation provides a practical prescription for computing the final curvature perturbation directly from horizon-crossing fluctuations, thereby establishing the anisotropic generalization of the $δN$ formalism. We further derive an explicit relation between curvature perturbations and primordial gravitational waves, demonstrating how anisotropic expansion couples scalar and tensor sectors on superhorizon scales. Our framework provides a practical basis for predicting statistical anisotropies in primordial scalar and tensor perturbations.

We study the motion of spinning test particles in Schwarzschild spacetime within the Mathisson--Papapetrou--Dixon pole--dipole approximation, imposing the Tulczyjew--Dixon spin supplementary condition. Restricting to equatorial orbits with the particle spin aligned with the orbital angular momentum, and retaining terms through linear order in the specific spin $s$, we derive the spin-corrected radial potential, circular-orbit conditions, bound periodic trajectories, epicyclic frequencies, and Lyapunov exponents of unstable circular orbits. The spin--curvature coupling shifts the circular-orbit energy and angular momentum and moves the innermost stable circular orbit to $r_{\rm ISCO}=6M-2\sqrt{2/3}\,s+\mathcal{O}(s^2)$ in the sign convention adopted here. We construct bound periodic orbits using the Levin--Perez-Giz zoom--whirl taxonomy and show how the particle spin deforms the corresponding energy--angular-momentum map. We then obtain the coordinate-time azimuthal and radial epicyclic frequencies and use them as kinematical inputs for relativistic-precession and resonance prescriptions for quasi-periodic oscillations. Finally, we relate the Lyapunov exponent of unstable circular orbits to the local separatrix structure governing near-homoclinic zoom--whirl motion. The resulting formulation provides a compact analytic connection between linear-in-spin MPD dynamics, periodic-orbit taxonomy, epicyclic-frequency shifts, and transient strong-field phenomenology in a nonrotating black-hole background. Also, we study the gravitational waveforms from the periodic orbits of a massive spinning particle around a black hole, presenting those associated with extreme mass-ratio inspirals involving a stellar-mass compact spinning object orbiting a supermassive black hole.

Obtaining accurate fluxes for faint sources from the ground at near- and mid-infrared wavelengths is challenging because of rapidly varying atmospheric absorption, and is often limited by the lack of a nearby spectro-photometric calibrator. We present the STellar Absolute Reference Spectroscopic Flux Library (STARSFLUX), an all-sky catalogue of 64,484 calibrator spectra spanning 0.3--30 $μ$m, based on the Mid-infrared stellar Diameters and Fluxes compilation Catalogue (MDFC). STARSFLUX combines Gaia DR3 stellar parameters, multi-band photometry from space- and ground-based surveys, and synthetic NewEra PHOENIX atmosphere models. For each star, we fit the observed SED with an interpolated model spectrum, angular diameter, and extinction to produce an absolutely flux-calibrated spectrum. The derived angular diameters can also be used to calibrate interferometric observations. We validate the catalogue using three independent tests. First, STARSFLUX angular radii agree closely with Gaia DR3 radii, with a median $|ΔR|/R_{\rm Gaia,DR3} \simeq 4.8$ per cent. Second, integrated L-band (2.8--4.2 $μ$m) fluxes for 12 stars in common with the Cohen infrared standards agree with Cohen values with a mean absolute percentage difference of $4.7 \pm 2.6$ per cent. Third, comparison with Gaia DR3 BP/RP spectra shows that the near-UV/visible/near-IR absolute spectro-photometric fluxes agree to approximately 3 per cent. The FITS spectra are available at https://home.strw.leidenuniv.nl/~gamez/ and will be submitted to VizieR.

Flux ratio anomalies in multiply-imaged quasar lenses are a long-standing issue. Using a classical system B1422+231 as a case study, we investigate how typical non-clumpy perturbations beyond elliptical shapes -- multipoles $m_3, m_4$ and radial variations in $q, φ_q$ -- can account for the observed image positions and flux ratios under different observational precisions. We extract these perturbations from a pre-selected strong-lensing galaxy sample from the TNG100 simulation. Smooth macroscopic models (SIE+$γ$, EPL+$γ$) are then fitted to the observed image positions alone and to both positions and flux ratios, with and without including the extracted perturbations. With astrometric uncertainty of $σ_{p}=10$ mas, both macro-models alone can already successfully fit image positions within $3σ_{p}$. At $σ_{p}=2$ mas, however, 'astrometric anomalies' appear if smooth macro-models alone are adopted. In this case, adding the extracted perturbations can explain the anomalous image positions. When both positions and flux ratios are adopted, the SIE+$γ$ model family already shows 'flux ratio anomalies' at photometric uncertainty $σ_{f} \le 10\%$ (keeping $σ_{p}=10$ mas). When EPL+$γ$ is used, the smooth model alone can simultaneously fit both positions and flux ratios with $σ_{f}=10\%, 5\%$, but not with $σ_{f}=2\%$, where 'flux ratio anomalies' appear. Adding all four types of extracted perturbations can rescue the macro-models and explain the observed anomalous flux ratios. We present important lessons learned regarding model flexibility and degeneracy.

This paper investigates a warm Dirac-Born-Infeld (DBI) inflationary model with nonminimal derivative coupling (NMDC) to gravity, where the inflaton kinetic term interacts with the Einstein tensor, thereby improving the effective gravitational friction. This model seamlessly integrates the noncanonical DBI kinetic structure, the NMDC-induced gravitational friction, and thermal dissipation. We formulate the background evolution equations along with the corresponding slow-roll stability conditions, leading to analytic results for the scalar spectral index $n_s$ and the tensor-to-scalar ratio $R$. By applying these results to power-law potentials with $n=2$ and $n=4$, the model parameter space is constrained using the \emph{Planck} 2018 data. The findings indicate that the interaction between NMDC-induced gravitational friction and thermal dissipation effectively modulates $n_s$ and significantly expands the viable parameter space. In the $(n_s,R)$ plane, the predictions for $N=50$ fall within the 95\% confidence-level region, while those for $N=60$ extend into the 68\% confidence-level region and approach the observationally preferred central values. For the representative parameter choices examined, the tensor-to-scalar ratio is notably suppressed, generally within the range $10^{-8}\lesssim R\lesssim10^{-5}$. Moreover, the combined damping mechanism relaxes the slow-roll condition related to $η$ and limits the inflaton field excursion, thus addressing the $η$ problem without incurring super-Planckian field variations. These results indicate that warm DBI inflation with NMDC offers a theoretically coherent and observationally viable inflationary model, showcasing the complementary effects of thermal dissipation and enhanced gravitational friction in the context of modified gravity.

Mutual coupling between antennas has emerged as the dominant direction-dependent corruption in dense aperture arrays, imprinting pronounced sub-MHz spatial and spectral structure that compromises the time-gating and foreground-separation strategies used to isolate the faint 21-cm signal. In this work, we introduce \textit{Direct Primary Beam Correction}, a domain-agnostic framework for reconstructing the far-field radiation pattern relative to an arbitrary reference via a regularised, direction-weighted linear inversion of stacked Jones matrices, thereby enabling the removal of direction-dependent distortions such as mutual coupling. Using full-wave electromagnetic simulations of SKA-Low, we demonstrate that this framework reconstructs the radiation pattern down to the numerical noise floor within a suitably conditioned field of view, with the reconstruction accuracy governed by the regularised inversion and the fidelity of the underlying beam model. Applying the framework to a simulated 4-hour observation of the EoR0 field in the $122$--$134$~MHz band, we identify two principal implications for 21-cm power-spectrum analysis. First, restricting the correction to the main lobe and near sidelobes is inadequate: chromatic grating lobe contributions, whether left insufficiently or entirely uncorrected, continue to contaminate the EoR window. Second, mutual-coupling-induced contamination is temporally coherent and, being anchored to the fixed array geometry, does not average down across snapshots as the EoR field is tracked. Direct primary beam correction, therefore, provides a computationally efficient means of mitigating mutual coupling; however, robust recovery of the EoR window necessitates either full-sky correction or explicit separation of main-beam and sidelobe contributions prior to power-spectrum estimation.

The origin of magnetic switchbacks-large-amplitude, spherically polarized magnetic-field fluctuations with local polarity inversions shown to be nearly ubiquitous in the inner heliosphere by Parker Solar Probe-is presently one of the most outstanding open questions in solar and heliospheric science. The occurrence of these structures in the young solar wind is a topic of active research, focusing not only on their characteristics and evolution but also on unraveling the puzzle of their origin. We first discuss the potential influence on switchback (SB) formation of large-scale coronal dynamics and restructuring processes that shape and regulate the extended solar corona. Following that, we review the dynamics, physical properties, occurrence rate, and energetics of numerous possible precursor small-scale dynamic events that occur throughout the solar atmosphere, emphasizing their possible roles in creating SBs. Finally, we discuss some recent studies attempting to connect in situ observations of SBs and SB patches with remote-sensing observations, and the related challenges in identifying solar wind source regions precisely and linking in situ observations to transient solar phenomena. We also clarify the terminology used by the solar community to describe small-scale solar phenomena.

Efficient, differentiable Poisson solvers are a key component of modern particle--mesh simulations and field-level inference pipelines. FFT-based solvers are extremely effective on fixed Cartesian meshes, but they impose global all-to-all communication and rely on symmetries that are lost in adaptive or non-Cartesian coordinates. In this work, we present a JAX-native geometric multigrid framework for particle--mesh gravity and argue that multigrid plays two complementary roles: on fixed meshes it can be a competitive, communication-avoiding alternative to FFTs, while on moving meshes it becomes the enabling solver. For static FastPM evolution, warm-started Chebyshev multigrid acts as a defect-correction method, exploiting temporal coherence between time steps to reduce the number of V-cycles required for field-level accuracy. At large mesh sizes this reduces memory pressure and yields comparable or faster wall-clock performance than distributed FFTs, with up to a factor of two reduction in total GPU time at fixed final mesh size. We then embed the same solver in a differentiable moving-mesh particle--mesh method, where adaptive coordinate deformation produces a variable-coefficient curvilinear Poisson equation that cannot be solved by ordinary FFT diagonalization. The resulting method concentrates force resolution in nonlinear structures while retaining a regular, JAX-compilable, automatically differentiable array workflow. These results suggest geometric multigrid can be a practical bridge between fast fixed-grid PM methods and differentiable adaptive-force cosmological simulations.

We make use of the free data formalism developed in \cite{CK25,CK26} to provide a direct construction of short-pulse type Cauchy data. The construction of the spacelike short-pulse follows from the local existence result established in \cite{CK26}. The forward integration construction in \cite{CK26} allows us to show that such data can be extended to a set of asymptotically flat Cauchy data. This greatly extends the result of Li--Yu \cite{LiYu}.

We present the analysis of four cool and strongly metal-polluted H-atmosphere white dwarfs observed with X-shooter. We compared their atmospheric parameters obtained from spectroscopy, photometry and a hybrid method, finding a difference of up to $\simeq160$ K in $T_\mathrm{eff}$, and $\simeq0.3$ dex in $\log g$, between the three analyses. We adopt the $T_\mathrm{eff}$ and $\log g$ from the photometric analysis, and subsequently measured the metal abundances of the photospheric elements from spectroscopic modelling, analysing their compositions. We identified from five to eleven unique metals in the photospheres of the four white dwarfs, with total accretion rates ranging from $10^{8}$ to $10^{9}$ $\mathrm{g~s^{-1}}$. The compositional analysis of WD J035826.49$+$215726.16 suggests an accreted planetesimal akin to a core-rich differentiated body, with a core mass fraction of 70 per cent, placing it among the most core-rich objects known to be accreted by white dwarfs. The parent body accreted by WD J042643.98$-$415341.44 shows an enhancement in Na compared to that of the Earth, making it most similar to primitive chondrites. The photosphere of WD J013222.88$+$052923.71 is greatly depleted in core-like material, and its composition resembles that of pure crust/mantle material. WD J232428.21$-$021643.65 has accreted the most Fe-rich planetesimal, with an Fe mass fraction of 67 per cent. We compare these results to other published studies and conclude that these white dwarfs are among the most heavily polluted cool DAZs studied to date, increasing the sample of cool H-dominated white dwarfs with five or more metals by 50 per cent.

Massive fields exchanged during inflation source cosmological collider features in the correlators of primordial curvature fluctuations. While their structure is fixed by the symmetries of inflation, the amplitude of such signatures is a priori unknown, though bounded by physical principles such as unitarity. We assess the observational viability of the scalar cosmological collider by combining precise numerical predictions derived using CosmoFlow with observational data from Planck bispectrum reconstructions. Working in the limit of weak quadratic mixing between the Goldstone and scalar sectors (as assumed by all previous searches), we find that perturbativity restricts the collider signals to be unobservably small, except for double- and triple-exchange configurations at low masses. Larger bispectra can be generated by increasing the mixing to non-perturbative (but still unitary) values: however, this distorts the signal, damping its amplitude and shifting the collider oscillations. Furthermore, strong mixing breaks the factorization of the bispectrum into a free amplitude and a fixed shape, and complicates both theoretical predictions and observational searches. Despite the challenges, we argue that strong mixing analyses are necessary to obtain meaningful constraints on the scalar cosmological collider from current and upcoming datasets.

We study the thermodynamic properties and radiation characteristics of a regular compact object obtained by applying the Simpson-Visser regularisation to the Schwarzschild black hole in modified gravity. The resulting SV-MOG spacetime, whose lapse function involves both the MOG coupling parameter and the black bounce parameter, smoothly interpolates between a regular black hole, a one-way wormhole, and a traversable wormhole depending on the parameter. We derive the Hawking temperature and heat capacity for both the black hole and wormhole branches, identifying second-order phase transitions signaled by sign changes in the heat capacity. Quantum gravitational corrections to the entropy are incorporated via logarithmic terms parameterized by coefficients, and we show that deviations from the Bekenstein-Hawking area law become significant at small horizon radii. For the radiation sector, we compute the electromagnetic flux, effective disk temperature, and spectral luminosity of geometrically thin accretion disks surrounding both black hole and wormhole configurations. Our results demonstrate that increasing the modified parameter enlarges the event horizon and enhances the emission, while increasing the black bounce parameter suppresses the horizon and softens the spectral profile, providing observational signatures that may distinguish SV-MOG compact objects from their Schwarzschild and pure SV counterparts.

The dark energy equation of state (EoS) and a possible dark-sector interaction are degenerate at the level of background expansion: the same expansion history may be interpreted as time-varying dark energy, energy exchange with dark matter, or a mixture of both. We introduce a data-driven, nonparametric expansion--growth framework that breaks this degeneracy by simultaneously reconstructing the dark energy EoS and the dark-sector coupling. Allowing the interacting matter density to evolve freely, we incorporate the coupling into the linear matter perturbation equation and reconstruct the expansion and growth histories using Gaussian processes, with hyperparameters marginalized in a Bayesian treatment. Applying this method to Pantheon+, cosmic chronometers, BAO measurements including DESI DR2, and redshift-space-distortion data, we infer both $w_{\rm de}(z)$ and the interaction history over $0\lesssim z\lesssim 2$, without assuming either a parametric EoS or a prescribed interaction form. We find no statistically significant evidence for either a nonzero interaction or dark energy dynamics: the reconstructed coupling and dark energy EoS remain consistent with the $Λ$CDM limit. Our results establish an expansion--growth consistency test for coupled dark-sector physics and provide a model-independent route to distinguish genuine dark energy dynamics from effective dark-sector energy exchange.

A. A. Friedmann (04.06.1888 -- 16.09.1925) proposed the first physical cosmological models in 1920s. Despite the fact that Friedmann's works were very famous soon after their publication, the study of dynamic models of the Universe in the USSR was actually banned in the 30s - 50s of the last century, and Soviet philosophers and propagandists wrote that models of the evolving Universe were invented by Lemaitre on the demand of the Roman Pope, because according to Soviet philosophers such the birth of the universe is very similar to the divine creation of the world described in the Bible. Thus, in the USSR, Friedman's works were in oblivion in since 1930s until 1960s.

Analytic infinite derivative gravity theories provide a renormalizable and ghost-free description of gravity around covariantly constant backgrounds. These theories can have non-singular bouncing Universe solutions. In this paper we aim to address a question whether it is possible to realize a bouncing solution without the presence of a ghost or a tachyon instability in this framework. We perform a detailed analysis of degrees of freedom in $(1+3)$ formalism around Minkowski and de Sitter space-times. As a result it becomes clear that on a very general basis one cannot construct an instability free bounce without a negative cosmological constant. An analysis of known bouncing solutions in this model shows that an analyticity of higher derivative form factors in combination with solutions parameters result in the presence of a ghost radiation. Being motivated by the idea of resolving the cosmological singularity problem we proceed by analyzing scalar and tensor modes anyway. Scalar modes appear to not influence Cosmic Microwave Background observations at all, while tensor modes spectrum is computed and the corresponding implications are discussed.

We solve the semiclassical Einstein equations for an evaporating Schwarzschild black hole in closed form, obtaining its backreacted geometry in 3+1 dimensions. This is achieved by a reformulation of the Hadamard renormalization prescription, which admits explicit stress-energy tensors for scalar quantum fields in 3+1-dimensional curved spacetimes without the standard mode-sum construction. The Unruh-like states we construct reproduce the Hawking flux at future null infinity and the associated ingoing negative-energy flux at the horizon, yielding the first analytical quasi-stationary solution describing black-hole evaporation in four-dimensional semiclassical gravity. The resulting geometry exhibits a timelike apparent horizon, providing a new analytical setting to study the causal structure of black-hole evaporation and its connection with the information-loss issue. Our results open a new route to semiclassical backreaction in four dimensions, extending the level of analytical control previously available only in two-dimensional effective models.

We present a photometric catalogue of the Local Group dwarf galaxy NGC 6822 based on deep JWST observations obtained with the Near-Infrared Camera (NIRCam) and Mid-Infrared Instrument (MIRI). Point-spread-function photometry and band matching were performed with StarbugII. The resulting catalogue contains 864,114 NIRCam point sources and 17,235 MIRI detections, with 10,079 detected in both instruments. Blackbody fitting yields effective temperatures and bolometric luminosities for 119,621 stars, providing a detailed view of the resolved stellar content. Candidate evolved stars were selected from NIRCam-MIRI colour-magnitude diagrams, with the sample refined by removing resolved contaminants and excluding young stellar objects through spectral energy distribution (SED) fitting. The final sample of 1226 evolved stars was analysed using the Grid of Red Supergiant and Asymptotic Giant Branch Models (GRAMS), from which dust-production rates were obtained and carbon-rich or oxygen-rich asymptotic giant branch (AGB) classifications assigned via a likelihood-weighted comparison across the full model set. Across the JWST fields, evolved stars return 5.6 x 10^-7 Msun yr^-1 of dust to the interstellar medium, with oxygen-rich AGB stars contributing 60% and carbon-rich AGB stars 35%, despite the low metallicity of NGC 6822. The unexpectedly high oxygen-rich contribution indicates the presence of intermediate-mass AGB stars undergoing hot-bottom burning, in line with the recent star-formation history of NGC 6822. Dust-producing AGB stars exhibit a centrally concentrated carbon-rich population and a more extended oxygen-rich population. We also identify JWST colour relations that provide robust photometric estimators of dust-production rates for evolved stars.

The majority of the Solar System planets are sources of bright radio emission, driven by energetic electrons trapped within each planet's magnetic field. Detection of this emission from exoplanets provides a unique opportunity to characterise their magnetic fields, which is key to determining the atmospheric evolution of exoplanets. However, a conclusive detection of radio emission from an exoplanet remains at large, primarily due to a lack of sensitivity at low radio frequencies. On the other hand, planet-like radio signatures have been detected on objects called ultracool dwarfs (UCDs) for over two decades. UCDs are of comparable sizes to Jupiter, but are more massive. They also possess similar interior structures to Jupiter, the region where magnetic fields are generated. Therefore, UCDs are ideal targets to study to advance our understanding of how magnetic fields manifest at planetary scales. In this Chapter, we outline the revolutionary role that the Square Kilometre Array will play in the study of exoplanets and UCDs. We anticipate that it will facilitate the first detection of radio emission from giant exoplanets with strong magnetic fields, and will deliver thousands of detections of UCDs within a few hundred parsecs. Combined with very long baseline interferometry, we also expect that astrometric monitoring will enable the detection of planets of a few Earth masses orbiting nearby radio-emitting UCDs. These findings will open a new window into how planets form and evolve in extrasolar systems.

We numerically study perturbations of spatially homogeneous and orthogonal solutions to the Einstein-Euler equations with positive cosmological constant and linear equation of state $p=Kρ$ for $K \in (\frac{1}{3},1)$. Recent numerical and analytic work has demonstrated that the fluid density gradient in these perturbed spacetimes can form steep gradients and blow up as future timelike infinity is approached, a phenomenon now known as the `Rendall instability'. In this article, we study the onset of this instability outside of one spatial dimension for the first time. Our results provide compelling evidence that the Rendall instability causes a forward turbulent cascade and $H^{s}$ norm inflation in the energy density of the fluid from arbitrarily small $U(1)$-symmetric perturbations. This is, to the best of our knowledge, the first example of small data turbulence for solutions of the field equations with a positive cosmological constant. Moreover, we find that this instability can drive structure formation in the late universe.

Characterising post-common envelope binaries (PCEBs) containing a white dwarf and a main-sequence companion is essential for improving theories of binary evolution. This paper presents the first direct spectroscopic confirmations of the white dwarf components in four PCEB systems within 20 pc of the Sun: G 203-47, GJ 207.1, LHS 1817, and Wolf 1130. To detect the white dwarfs we obtained near-UV spectroscopy from STIS on the Hubble Space Telescope, fitting with white dwarf models and M dwarf proxy spectra. We provide estimates of the white dwarf effective temperatures, which range from approximately 5300 K to 6300 K. We compare these parameters to those determined from modelling with photometry alone, and find a 5 - 8 per cent discrepancy, due to emission features. Notably, 27 years after its initial detection, we confirm the presence of a white dwarf in G 203-47, which is the ninth closest white dwarf to the Sun. Using Swift XRT data, we find that despite the 14.9-day orbital period of G 203-47, it is not tidally locked, possessing a rotation period likely exceeding 100 days, and making it a rare example of a long-period PCEB formed via a brief common envelope interaction. We update the local white dwarf space density to (5.2 $\pm $0.4) $\times$ 10$^{-3}$ pc$^{-3}$, and compare our results to models from the Binary Populations and Spectral Synthesis (BPASS) framework, finding a good agreement with the predicted and observed numbers of PCEBs within 20 pc.

High-redshift binary black hole (BBH) events are promising candidates for primordial black holes (PBHs) detectable by next-generation gravitational wave (GW) detectors. A redshifted mass distribution of detected PBH candidates can be obtained from GW observations, from which the underlying PBH mass function can be reconstructed. In this work, we develop a framework that applies the gradient-descent method to the observed redshifted mass distribution and reconstructs the PBH mass function and, subsequently, the primordial power spectrum (PPS) on small scales. As an illustrative application, we analyze BBH events in the LIGO--Virgo--KAGRA (LVK) catalogs under a specified PBH selection criterion. We find a regularization-stable candidate bump-like enhancement of order $\mathcal{O}(10^{-2})$ in the reconstructed PPS, centered around $k_{\mathrm{peak}}\simeq 5.7\times 10^5~\mathrm{Mpc}^{-1}$ under the adopted assumptions. Our results demonstrate the feasibility of reconstructing the small-scale PPS from high-redshift BBH observations with next-generation GW detectors.

Recent analyses of Proca theories with non-minimal curvature couplings have uncovered an additional scalar degree of freedom with an identically vanishing propagation speed, $c_s^2=0$, signaling a scale-dependent strong-coupling problem. In this {\it Letter}, we show that an additional derivative interaction can lift this zero-speed degeneracy and restore a non-degenerate quadratic dynamics for the scalar perturbations. In an open region of parameter space, the scalar kinetic matrix is positive definite and the high-frequency propagation speeds are real and positive. At quadratic order, this removes the specific signature of strong coupling associated with the $c_s^2=0$ mode. We also uncover a non-uniform limit as the temporal vector condensate approaches zero: at fixed nonzero $A_0$, the formal extreme-ultraviolet regime develops a ghost, while the kinetic structure of the exactly vanishing-condensate branch is different. The momentum scale at which this ghost appears is pushed toward increasingly high values as $A_0\to0$. Whether the ghost scale ultimately lies above the EFT cutoff depends on an independent determination of the EFT cutoff. Finally, we identify a stable de Sitter fixed point with $A_0\neq0$, surrounded by a finite region in which the scalar no-ghost and high-frequency stability conditions remain satisfied.

We investigate the thermodynamic modifications of the Reissner-Nordstroem anti-de Sitter (RN AdS) black hole induced by a minimal length scale,which naturally emerges in $κ$-deformed space-time. By constructing the modified metric via phase-space commutation relations,we derive the deformed Hawking temperature, entropy, and enthalpy. We analyze the thermal stability through the heat capacity and study the P-V criticality, revealing that the black hole undergoes a small-to-large phase transition analogous to the Van der Waals system, albeit with a critical ratio slightly lowered by non-commutativity. Furthermore, we examine the Joule-Thomson expansion and find that the non-commutative (NC) parameter expands the cooling region in the temperature-pressure plane. Our results demonstrate that while the overall thermodynamic analogy with the Van der Waals fluid persists, the minimal length scale systematically deform the coexistence region and inversion curves, offering potential observational signatures for quantum gravity.

Recent Gaia studies have identified numerous open clusters (OCs) & tidal tail catalogues, enabling systematic searches for white dwarfs (WDs) associated with clusters & their extended structures. We compile a literature-based sample of OC-WD pairs to validate WD membership in cluster cores & tails, investigate the initial-final mass relation (IFMR), identify WDs formed through non-canonical evolution, and interpret the observed WD populations using a grid of N-body simulations. We combine Gaia DR3 cluster & tail catalogues with UV-IR photometry to analyse the OC-WD pairs. WD masses, cooling ages, radii, temperatures & luminosities are estimated using colour-magnitude diagrams & spectral energy distributions. These observations are interpreted in the context of N-body simulations. We identify 235 OC-WD pairs (99 in tails) in 80 clusters. More than 28% of the pairs are likely spurious, with contamination substantially higher in the tails (>48%) than in the cluster cores (>13%), indicating significant field-star contamination in current Gaia-based catalogues. The Pleiade tails also show severe contamination by old WDs. Simulations predict that the fraction of core WDs increases with cluster age, reaching >10%, whereas the observed fractions remain systematically lower, consistent with the WD deficit problem. Despite the high contamination rate, most tail WDs (~83%) are consistent with having been born inside the tidal radius. We also identify 63 candidate binary-origin WDs & 47 new IFMR candidates. WDs provide a powerful probe of contamination in cluster and tail catalogues and place important constraints on cluster detection methods & N-body simulations. Resolving the WD deficit and improving membership validation will require improved observations, membership methods, WD physics, and spectroscopic follow-up, enabling stronger constraints on dynamical cluster evolution & the WD IFMR.

We investigate the presence of compact, long-range, and vacuumless regimes in models described by real scalar fields in 2D flat and 5D warped geometries. In the two-dimensional case, we study in detail several models that support localized solutions. When embedded into a warped $\text{AdS}_5$ background, these fields generate distinct kinds of brane profiles.

We determine the quantum corrections to transport arising from fluctuations of the near-AdS${}_2$ throat of near-extremal black branes in holographic models with momentum relaxation. By computing the shear viscosity and electrical conductivity at both zero and finite chemical potential, we uncover a universal low-temperature enhancement of transport generated by Schwarzian quantum fluctuations. Specifically, transport coefficients extracted from retarded Green's functions increase throughout the regime $C T \ll 1$. For operators with Schwarzian scaling dimension $Δ>1$, this enhancement is preceded by a universal minimum at a characteristic temperature $T_{\rm min} \propto C^{-1}$, leading to a non-monotonic temperature dependence strikingly similar to that observed in many correlated materials. In contrast, for the case $Δ=1$, relevant for the electrical conductivity, the transport coefficient evolves monotonically toward a constant value. Our results identify universal signatures of near-horizon quantum gravity in the transport properties of holographic quantum matter.

We revisit holographic electrical conductivity in the regime of very low temperatures where quantum fluctuations in the throat of the near-extremal, asymptotically AdS$_4$ dual black brane are strongly coupled. We treat the fluctuations via an effective two-dimensional action capturing the effects of the Schwarzian modes at the scale $1/C$. Our main result is a non-monotonic temperature dependence: the quantum-corrected conductivity decreases as the temperature is lowered, reaches a minimum at $(CT_{\mathrm{min}}\simeq 0.023)$, and grows as $(1/\sqrt{CT})$ for $(CT\ll 1)$. We further estimate the one-loop contribution to the conductivity directly from the four-dimensional gravitational path integral and find qualitative agreement with the effective two-dimensional description in the regime $(CT\gg 1)$. We find that quantum corrections substantially modify the low-temperature behavior of the holographic conductivity in both approaches.

The unification of quantum mechanics and general relativity remains one of the major open problems of theoretical physics. The Anti-de Sitter/Conformal Field Theory (AdS/CFT) correspondence provides a valuable theoretical framework for this effort via a holographic duality between a theory of quantum gravity in asymptotically AdS spacetime and a conformal quantum field theory on the lower-dimensional boundary. Here, we implement a toy model of this duality called the HaPPY code, a quantum error-correcting code in the form of a tensor network with hyperbolic entanglement patterns, on a trapped-ion quantum computer. We present the first experimental confirmation of the Faulkner-Lewkowycz-Maldacena formula in this model - a key test of the holographic correspondence. We then enrich it with non-stabilizerness, or magic, and observe entropic precursors expected of emergent gravity. Finally, we present and measure a code construction whose entropic behavior is reminiscent of a highly quantum wormhole. Our experiments illustrate how quantum computers can serve as testbeds for modeling the emergence of spacetime.

Quantum theory challenges the view that individual measurement outcomes are predefined and independent of the measurement context. Yet the quantum state itself -- the catalogue of probabilities for all possible measurements -- is usually assumed to be well defined. We argue that this assumption tacitly relies on measurements being performed relative to ideal, infinitely-resourceful reference frames. We show that, when measurements are made relative to non-ideal quantum reference frames, the probabilities themselves become indefinite: even in the limit of arbitrarily large number of runs, the relative frequencies may remain uncertain. The uncertainty is irreducible in a quantum-mechanical sense, as we show by proving a Bell-type theorem for relative frequencies. We further propose a quantum-optical implementation of these relational measurements based on pulsed homodyne detection. Our findings motivate an extension of the notion of the quantum state to regimes constrained by finite resources. We expect them to be especially relevant at the interface between quantum theory and general relativity, where the resources and information available in a bounded region of spacetime are fundamentally limited.

Light but Massive Relics (LiMRs), cosmological relics that are relativistic at some point in the observable early universe but are non-relativistic today, are a generic prediction of many theories of physics beyond the Standard Model, and are a target of high interest for current and upcoming cosmological surveys. The enormous variety of scenarios that can give rise to LiMRs also gives rise to a wide range of possible phase space distributions for these relics, making it valuable to have a model-agnostic quantification of LiMR signatures. Toward that end we identify three independent physical quantities that govern LiMRs' dominant effects on cosmological observables, and introduce a phenomenological family of "Squished Relics" (SquiRels) that allow us to efficiently and flexibly search for the imprint of LiMRs with different phase space distributions in data. Our primary focus here is on LiMRs that transition from radiation to matter during the Cosmic Microwave Background epoch. Using this framework we explicitly demonstrate that current Planck observations are not sensitive to the shape of the LiMR's phase space distribution function, which allows us to provide a model-insensitive and readily reinterpretable limit on LiMRs that become non-relativistic at redshifts z_{NR}<10^5. We further demonstrate that upcoming and future ground-based observatories will be able to distinguish LiMRs from massless dark radiation species, and could begin to provide information about the shape of its distribution.

The snowlines of major volatiles in protoplanetary disks play a pivotal role in dust evolutions and volatile delivery to nascent planetary systems. In this paper, we report the Atacama Large Millimeter/submillimeter Array Band 1 ($\approx7.5\,\mathrm{mm}$) observations of methanol (CH$_3$OH) emission lines in the disk around the FU-Ori type star V883 Ori, where accretion outburst heats the disk and the majority of ices has sublimated. We detect three CH$_3$OH emission lines at an angular resolution of $\approx0.\!\!^{\prime\prime}2$. The stacked CH$_3$OH image exhibits a centrally-peaked morphology in contrast to the previous (sub-)mm observations that show a central depression. By fitting radially-resolved line profiles, we derive the radial intensity profile of the CH$_3$OH emission where we find a steep increase at $\lesssim40\,\mathrm{au}$. The column density of CH$_3$OH reaches at least $\sim10^{19}\mathrm{-}10^{20}\,\mathrm{cm^{-2}}$ at $\sim10\,\mathrm{au}$. This provides direct evidence that a significant amount of warm gaseous methanol is present in the innermost region of the disk where its emission has been suppressed in previous (sub-)mm observations due to the optically thick dust emission. The steep increase in the intensity profile may indicate that the CH$_3$OH snowline in the midplane is located at $\sim30\mathrm{-}55\,\mathrm{au}$, or that the CH$_3$OH emission traces the temperature structure given that the emission is likely optically thick. Our results demonstrate the capability and significance of (sub-)cm observations in probing the innermost opaque region of disks, paving the way for the future observations with upcoming facilities.

We analyse the initial-value problem of metric higher-derivative effective theories of gravity. We show that any such theory whose characteristic velocities are independent of derivatives of the metric is intrinsically weakly hyperbolic, independently of the gauge fixing. To show this, we identify the spin-$2$ physical sector directly from the characteristic equation; this can be done without introducing an order-reduced formulation, which greatly simplifies the computation. In this sector, every metric theory with more than two derivatives in the equations of motion contains a weakly hyperbolic block. Since this obstruction is physical, no choice of gauge or constraint addition can remove it, providing a structural explanation for the failure of strong hyperbolicity in this broad class of theories.

We investigate linear and nonlinear primordial perturbations in an anisotropic Dirac-Born-Infeld (DBI) inflationary model with a non-vacuum initial state. Using the Arnowitt-Deser- Misner (ADM) formalism, we expand the action up to second and third order in the curvature perturbation and derive the corresponding scalar and tensor power spectra, as well as the bispectrum and the equilateral non-linearity parameter \(f_{NL}^{\mathrm{equil}}\). The effects of anisotropic corrections and non-Bunch-Davies (non-BD) initial conditions are incorporated through the slow-roll sector and Bogoliubov coefficients. For the numerical analysis, we consider an intermediate expansion scenario together with a phenomenological ansatz for the excited-state occupation number \(N_k\). By comparing the model predictions with recent observational datasets, including Planck2018 TT, TE, EE + lowE + lensing + BK18 + BAO and DESI+CMB+DESY5 data, we identify observationally viable regions in the parameter space of the model. Our analysis indicates that the anisotropic DBI scenario with non-vacuum initial conditions can remain compatible with current constraints on the scalar spectral index, tensor-to-scalar ratio, and equilateral non-Gaussianity for suitable ranges of the anisotropy parameter \(c\) and the initial-state parameter \(N_{k,0}\).

Cosmic voids are underdense regions that can provide an effective large-scale environment with a de Sitter-like gravitational behavior. Motivated by recent black-hole solutions embedded in void density profiles, we construct a symmetric thin-shell wormhole by gluing two copies of the positive-lapse region of a black hole inside a cosmic void. The surface stresses are obtained from the Darmois--Israel junction conditions, and the corresponding null, weak, dominant, and strong energy-condition combinations are written directly in terms of the void mass function and density profile.We further develop the thermodynamics of the static shell, deriving a first law that relates the shell entropy to the black-hole and cosmological-like horizon entropies.We then formulate the radial dynamics of the throat through an effective potential and study the local stability of static configurations when the exotic matter on the shell obeys either a generalized cosmic Chaplygin gas or a modified cosmic Chaplygin gas equation of state. In both models the Chaplygin parameter $B$ is fixed by the static junction condition, so that the remaining stability test is governed by the void geometry and by the equation-of-state parameters. Numerical results reveal that GCCG-supported configurations are generically unstable, whereas MCCG-supported wormholes can be stable for a sufficiently large linear term in the equation of state. The resulting framework connects the de Sitter-like structure of cosmic voids with the standard thin-shell wormhole formalism and provides a starting point for identifying stable or unstable wormhole configurations located between the black-hole and cosmological-like horizons of the void spacetime.

We investigate the modified thermodynamic properties of a deformed Schwarzschild black hole arising from new noncommutative corrections to the tetrad fields obtained through the Seiberg-Witten map. By deriving the corrected event horizon radius, we evaluate the corresponding thermodynamic quantities, including the Hawking temperature, entropy, heat capacity, Helmholtz free energy, and pressure. We also examine the evaporation process of the black hole within this noncommutative framework. The behavior of the heat capacity reveals a rich thermodynamic structure consisting of three distinct phases: two unstable phases with negative heat capacity, corresponding to a large black hole and a small black hole, separated by an intermediate stable phase with positive heat capacity. Our analysis also shows that the Schwarzschild black hole in noncommutative geometry possesses a finite lifetime that exceeds that of its classical counterpart, owing to the effective contribution of noncommutativity to the black hole mass. The noncommutative correction naturally introduces a fundamental length scale, Θ=1,6.10^-35 m, which is of the order of the Planck length.

In this paper we present new spectropolarimetric and velocimetric observations of the prototypical classical T Tauri star TW Hya obtained with SPIRou at the Canada-France-Hawaii Telescope, expanding our previous monitoring over two new seasons (2024 and 2025). We confirm that the large-scale magnetic field of TW Hya varied with time, and find that it showed fluctuations on a timescale of about a year in addition to the longer term variations outlined in the previous study. Using Zeeman-Doppler imaging, we obtain that the large-scale field of TW Hya mostly consisted of a poloidal dipole of mean polar strength 0.83 kG, inclined at an average 17degr to the rotation axis. We also find that the radial velocities of TW Hya, once fully filtered from telluric contamination, were dominated by rotational modulation induced by activity, with residuals of 32 m/s rms. No signal from a putative close-in planet is found, with an upper limit on the planet mass ranging from 0.33 to 0.98 Mjup for distances of 0.053 to 0.41 au from the central star. Emission lines indicate that the mass accretion rate was equal to $10^{-8.33\pm0.20}$ Msun/yr on average, with peak-to-peak fluctuations by a factor of ~5 from season to season. This confirms that accretion onto TW Hya is unstable, with the magnetospheric gap carved by the large-scale field at the center of the disk extending on average no further than 33-40% of the corotation radius where the disk Keplerian angular velocity equals the rotation rate at the stellar surface.

We study timelike entanglement entropy and timelike subregion complexity in localized black holes with asymptotic AdS3*S3*T4 geometry, focusing on the black-pole solution. Unlike the BTZ solution, the black pole exhibits a nontrivial dependence on the internal sphere through the functions $K_y(r,θ)$ and $G(r,θ)$. Both observables are constructed from spacelike and timelike Lorentzian branches, but they probe the geometry in different ways: timelike entanglement yields a complex lifted area, while timelike complexity gives a real, finite renormalized volume. We employ a localized timelike prescription in which the branch profile is built at an angular label $θ_0$ and subsequently lifted over the physical internal angle $θ$. In the large-$r$ regime, the leading angular dependence drops out, recovering the expected short-interval behaviour. In the exact black-pole geometry, the temporal families become non-monotonic, making a fixed-boundary-interval selection essential. As the boundary interval increases, the selected branches move inward and become sensitive to the localized cap-horizon transition region. These results demonstrate that timelike Lorentzian observables probe localized-geometry effects that are absent in BTZ and in the leading large-$r$ description.

The solar chromosphere is permeated by complex magnetic fields that guide plasma flows and energy into the corona. This work presents a detailed analysis of a unique, high-resolution observation of a giant chromospheric spiral structure that emerges due to a large magnetic pore, captured by the Swedish 1-m Solar Telescope (SST). A comprehensive data analysis pipeline is developed to automatically detect the edges of 2255 plasma flows (loops) that constitute the spiral, and these are used to extract the kinematics of flows propagating along the magnetic field. The analysis reveals three primary insights into the spiral's physics. First, magnetic curvature is correlated with oscillatory flow dynamics, i.e. regions of high loop curvature exhibit a statistically significant excess of higher-order oscillation modes compared to straighter loops; it is also correlated with higher intensity and longer periods. Second, spatial distribution of oscillation period shows an inverse trend, decreasing from $\sim$3.5 minutes in the pore to $\sim$3 minutes in the outer spiral arms. This is interpreted as a signature of the overlying trans-equatorial quadrupolar coronal loop system compressing the pore's field lines into a near-horizontal orientation, producing a period gradient that challenges the standard expanding canopy model. Finally, the emission signature confirms that oscillating threads represent localised channels of brightness that lie within cooler, absorbing loop material. This study provides the first statistical analysis of oscillatory flows in a large-scale spiral, probing energy flow through the chromosphere through curved magnetic structure.

Brown dwarfs are massive analogues of extrasolar giant planets. Compared to exoplanets whose observations are generally limited by the presence of their bright host star, brown dwarfs are ideal targets for studying substellar atmospheric physics, chemistry and dynamics. Previous observations and simulations of their atmospheres suggest preferential cloud formation around the equator, associated with an equator-pole thermal gradient. Here we show that this atmospheric structure should induce latitudinal chemical variations detectable by the Doppler effect. We introduce a new method - Differential Molecular Rotational Broadening - which consists in comparing the apparent rotational broadening of individual molecules from high-resolution spectra. Application of this approach to VLT-CRIRES observations for different molecules (CO, H2O, CH4 and NH3) in the atmosphere of the brown dwarf DENIS J0255-4700 confirms the existence of latitudinal chemical variations. Our data suggest a depletion of CH4 and NH3 at low latitudes, consistent with an equatorial cloud belt. Our method could be applied to multiple brown dwarfs and exoplanets to map their atmospheres and to study various atmospheric processes.

HD 98800 is a young hierarchical quadruple system composed of the tight binaries AaAb and BaBb on a wide, highly inclined outer orbit. The B subsystem hosts a circumbinary disc in a polar configuration, and the geometry of the system offers a rare opportunity to observe the passage of the disc around BaBb in front of AaAb. We aim to update the overall orbital solution to offer precise time windows of the main occultation features by combining the orbit with the state-of-the-art knowledge of the disc's structure. We combine new and published radial-velocity measurements and multi-wavelength astrometry for the outer AB orbit and both inner subsystems, AaAb and BaBb, in a joint orbital fit. The revised solution is consistent with previous dynamical masses and improved the outer orbit, reducing the uncertainties in the period and periastron epoch by about a factor of two. It also narrows the predicted crossing-phase windows to 5--15 days at the 1σ level, improving the timing of ingress, egress, and the first cavity-crossing predictions. These results provide a more accurate timing framework for future observations of the occultation, although the predicted epochs remain model-dependent because of uncertainties in the disc structure.

We consider simple one dimensional models of acoustic black holes formed by Bose-Einstein condensates where the flow is stepwise homogeneous. We shall concentrate on the case in which an extended sonic region is present, i.e. a region where the flow velocity equals exactly the speed of sound. These models can be thought of as acoustic black holes with a thick horizon. Particular attention will be devoted to the number of created particles and to the density-density correlation functions, these latter being at present time the experimental tool to identify Hawking like radiation in acoustic black holes.

We derive the Belinfante--Rosenfeld symmetrization procedure from the metric-affine conservation law by means of an affine lift of flat-spacetime field theories. Following minimal coupling to an independent affine connection, variation of the action with respect to the connection defines the hypermomentum current. Taking the flat-spacetime limit of the resulting conservation law, we recover the Belinfante--Rosenfeld relation, with the divergence of the hypermomentum current reproducing the Belinfante improvement term. We thoroughly study the form of the couplings with this property and their physical significance. The construction is illustrated for Quantum Electrodynamics and Quantum Chromodynamics.

We investigate the thermodynamic and optical properties of Kerr--Bertotti--Robinson black holes, namely rotating black holes immersed in an external Bertotti--Robinson electromagnetic background. In the fixed-$a$ ensemble, we derive the horizon mass relation, the Hawking temperature, the entropy, the Helmholtz-type free energy, the heat capacity, and the extremal remnant configuration. These quantities reduce smoothly to their Kerr counterparts as $B\to0$. In the weak-field regime, the leading thermodynamic corrections arise at order $B^2$; the extremal radius is shifted at this order, whereas the remnant mass receives its first correction only at order $B^4$. We also introduce a formal AdS-like thermodynamic interpretation of the Bertotti--Robinson scale, treating the associated pressure as an effective response variable rather than a genuine cosmological pressure. Because the spacetime is not asymptotically flat, we further compute the finite-radius Komar mass and the Komar charge associated with the horizon generator. Using the Hamilton--Jacobi formalism, we derive the separated null-geodesic potentials, the impact parameters of spherical photon orbits, and the celestial coordinates of the shadow boundary for a finite-distance observer. We then characterize the photon-region boundaries, ergosphere thickness, photon--ergosphere gap, shadow area, and magnetic shadow susceptibility. Within the perturbative regime considered, the Bertotti--Robinson background decreases the averaged ergosphere thickness and shadow area, increases the photon--ergosphere gap, and produces a negative shadow susceptibility whose magnitude is enhanced by rotation.

SKAO surveys will provide an unprecedented window into the large-scale structure of the universe through HI 21cm galaxy and intensity mapping surveys, and radio continuum surveys. We present forecasts for the cosmological constraining power of ``$N\times2\mathrm{pt}$'' analyses -- which combine galaxy clustering, galaxy weak lensing, galaxy-galaxy lensing signals and 21cm Intensity Maps. By assuming cosmology surveys from an SKA-Mid AA4, we show that such an $N\times2\mathrm{pt}$ analysis will be able to deliver measurements of $\sim1\%$ precision on $Λ$CDM cosmological parameters. We also explore dynamical dark energy in the $w_0, w_a$ model, the sum of neutrino masses $\rm M_ν$, and the background curvature $Ω_{\rm k}$.

The intrinsic alignment (IA) of galaxies and the dark matter haloes that host them is one of the leading astrophysical systematics for weak-lensing cosmology, yet how the IA signal depends on the large-scale environment in which haloes reside is not yet fully characterised. We use the high-resolution $N$-body simulations of the Dark Quest suite to measure the environmental dependence of the IA of dark matter haloes over the redshift range $z=0.1$--$1.5$. We quantify each halo's environment through the overdensity $δ_8$, defined from the number of neighbouring haloes within $8\, h^{-1}Mpc$, and we isolate the environmental effect from its degeneracy with the halo-mass dependence by comparing the most overdense and most underdense haloes constructed to share the same halo-mass distribution. We find that haloes in underdense environments exhibit systematically larger IA amplitudes $A_{\rm IA}$ than haloes of the same mass in overdense environments, by a factor of $\sim1.5$--$1.8$, and that this trend persists across the mass and redshift ranges probed, strengthening towards low redshift. Using an orientation-only (unit-ellipticity) estimator, we further show that this environmental contrast is driven by a combination of two effects: haloes in underdense regions are both intrinsically less spherical and more strongly aligned with the large-scale tidal field than their overdense counterparts of the same mass. These results indicate that the large-scale environment is a non-negligible variable in modelling halo and galaxy alignments, and may be a particularly important factor for beyond-two-point weak-lensing analyses.

It is widely believed that global symmetries cannot exist in a consistent theory of quantum gravity. A prominent mechanism underlying this expectation is provided by gravitational instantons or Euclidean wormholes, whose contributions to the Euclidean path integral generates local symmetry-breaking operators in the effective action with couplings parametrized by the $α$-parameters $\overrightarrowα=(α_1^{},α_2^{},\cdots)$. The appearance of these $α$-parameters is often taken as evidence that wormholes explicitly break global symmetries. In this paper, we argue that this conclusion is premature. The fate of the symmetry is determined by the ensemble average over these $α$-parameters. We show that, under broad situations, this average can be sharply dominated by the symmetric point $\overrightarrowα=0$, while contributions from other symmetry-breaking critical points $\overrightarrowα\neq 0$ are typically suppressed by a doubly exponential factor $\exp(-e^{2S_{\rm ins}})$, where $S_{\rm ins}$ is the instanton action. In particular, standard $\mathrm{U}(1)^{}$ Peccei--Quinn models fall into this class, implying that the conventional axion quality problem does not arise in the wormhole-induced effective theory. Our analysis is formulated for general $\mathrm{U}(1)^{}$ $p$-form global symmetries.