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

We report new optical spectroscopic results for the gamma-ray blazar PKS~0903$-$57. The optical properties of this source have historically been complicated to establish, due to a bright foreground star separated from the blazar by only 0.67$^{\prime\prime}$. During its strong high-energy (HE) and very high-energy (VHE) gamma-ray flaring period between 2020 and 2022, the Southern African Large Telescope and Very Large Telescope (VLT) took spectra of its optical counterpart. Photometric monitoring with the Rapid Eye Mount (REM) telescope complemented the observations. The observational results were inconclusive. A major improvement was achieved with a repeat of the VLT observations in 2024 using a 0.5$^{\prime\prime}$ slit under excellent seeing conditions. This led to the detection of five narrow emission features in the spectra, yielding a robust redshift measurement of $z$ = 0.2621 $\pm$ 0.0006. We also detected a broad symmetric H$α$ emission component with a full width at half maximum of 4020 $\pm$ 30 km/s in the spectra. Non-thermal jet emission was found to dominate the continuum, supporting the classification of PKS~0903$-$57 as a flat-spectrum radio quasar. Another interesting observation was a significant velocity displacement of $\sim$1500~km/s between the broad H$α$ line and its corresponding narrow line. Such a displacement may be associated with unusual dynamics in the broad-line region or accretion flow, a system containing two supermassive black holes, or a recoiling black hole produced following a supermassive black-hole merger. Building on these results, we carried out additional VLT, REM, and \emph{Swift} observations towards the end of 2024 and during 2025. We present in this paper the preliminary results from the analysis of these new data.

Electric charges carried by compact objects can affect binary dynamics and imprint characteristic signatures on gravitational-wave signals. We derive next-to-leading-order post-Newtonian frequency-domain waveforms for charged compact binaries in both the gravitational-quadrupole and electric-dipole dominated regimes, including electromagnetic corrections beyond leading-order electric-dipole radiation that reduce the degeneracy between the component charge-to-mass ratios. Using gravitational wave events from GWTC-5.0, we place constraints on the charge-to-mass ratio and mass of each component of the binary system.

When propagating through the universe, gamma rays at very-high energy (VHE, E > 100 GeV) and ultra-high energy cosmic rays (UHECRs, E > 1 EeV) can interact with the optical, infrared and microwave photon fields that permeate the universe. These interactions result in a characteristic absorption imprint in the spectra of extragalactic gamma-ray sources at VHE, and in a change in the mass-composition of UHECRs. The study of both VHE gamma rays and UHECRs therefore requires precise knowledge of the intensity of the cosmic photon fields. In this work, we explore the impact of the current uncertainties on the optical and infrared photon fields on the propagation of astroparticles. We restrict the range of available models to those best matching the recent measurements, and compare the different reconstructions resulting from these models. We find that the knowledge on the cosmic background light is no longer the dominant source of uncertainties in understanding the phenomenology of astroparticle sources in the low redshift universe (z < 0.1), enabling robust spectral and composition inference with the next generation of both gamma-ray and UHECR measurements.

The closest and brightest millisecond pulsar (MSP), PSR J0437-4715, is a high-priority pulsar for precision timing experiments, including the Parkes Pulsar Timing Array (PPTA) project, where it is used for both pulsar timing and for polarization calibration. Using approximately seven years of observations obtained with the Ultra-Wideband Low-Frequency receiver on Murriyang, CSIRO's Parkes radio telescope, we investigate the long-term spectro-polarimetric variability of PSR J0437-4715. We report the discovery of two consecutive discrete profile change events exhibiting the same phase-localised morphological variations, the first such behaviour reported in any MSP. Both events affect three phase-localised regions of the pulse profile that exhibit coupled and anti-correlated temporal evolution. We also identified, and corrected for, systematic profile variability induced by polarization calibration errors arising from long-term drift of the polarimetric properties of the receiver. The profile change exhibits frequency-dependent evolution in total intensity, together with variations in linear polarization, position angle and fractional linear polarization. These results confirm a magnetospheric origin for the events rather than propagation effects caused by the interstellar medium, and provide potential new constraints on MSP emission models.

We have developed a photometric method that produces clean hot (intrinsically blue) populations using Gaia and 2MASS photometry alone with [a] a very small (false positives < 1%) contamination from cool (intrinsically red) luminous sources, [b] independent of extinction (limited only by S/N), and [c] with only small losses (false negatives < 10%) of real OB stars. Applying the method to Gaia DR3 data we are able to recover extinguished OB stars down to E(G_BP-G_RP) = 5, producing a catalog of ~2x10^5 candidate OB stars that reaches beyond the Scutum-Centaurus arm. This method can potentially extend the sample of known OB stars by an order of magnitude, which would yield a large leap in our knowledge of the populations of compact-object progenitors, including gravitational-wave sources.

Prominences and coronal rain are two manifestations of coronal condensations whose formation mechanisms remain only partially understood. While previous studies have explored how localised steady or stochastic heating can trigger thermal instability and condensations, the role of wave driven energy transport and dissipation has not been extensively investigated. In the solar atmosphere, Alfven and kink waves are an abundant source of energy capable of contributing both to coronal heating and plasma structuring, making them compelling candidates in the formation processes of the structures we observe in the corona. We aim to explore how Alfven and kink wave energy and their dissipation influence the formation and evolution of condensations in coronal magnetic structures. In particular, we investigate each waves role independently of parameterised external heating mechanisms to isolate whether either wave can trigger condensations. We perform 2.5D simulations using the open source MPI AMRVAC code extended with the newly developed UAWSoM module. This module self consistently evolves Alfven and kink wave energy and their feedback on the plasma through wave pressure and heating terms. We examine the conditions under which condensations form, their morphology, and their dynamical evolution. Our results show that for our range of parameters, in 2.5D MHD, Alfven waves did not trigger the formation of prominences or condensations. We find that if the radius of the fine scale structures supporting the kink waves is of the order of 100 km, then kink waves can trigger coronal rain that differs considerably in its morphology from coronal rain triggered by parameterised heating functions. Through the use of UAWSoM, we have established that Alfven and kink waves leave distinct thermodynamic signatures on coronal loops, given the same wave energy injection.

We derive an extension of Chandrasekhar's pulsational stability criterion for relativistic stars with a pressure boundary condition at the surface. We then apply the new criterion to hylotropes, a model describing rapidly accreting supermassive stars with a core-envelope structure. We find that confinement typically stabilises stars, delaying the onset of the general relativistic instability when the external pressure at the surface is roughly $\sim10^{-3}$ of the stars' hydrostatic pressure scale. Beyond this value, the critical mass at marginal stability increases rapidly, exceeding the value of $\sim 10^5$ M$_{\odot}$ for helium burning, isolated $n=3$ polytropes by orders of magnitude. We discuss how the required pressure boundary conditions may be produced by the combined ram and magnetic pressure supplied by the accretion flow as it assembles the supermassive star, and provide a fitting formula for the final mass before collapse. As a rough estimate assuming equipartition between magnetic and thermal pressures in a standard $α$-disk, we find that the critical mass doubles at accretion rates of order $\dot M\sim100\,{\rm M_\odot}\,{\rm yr}^{-1}$, again increasing rapidly after this value. Our results show under what conditions, potentially realized in a subset of high redshift halos, non-rotating SMS may collapse into black holes beyond the typical scale of heavy seeds.

The interpretation of gravitational wave sources depends on the specific choice of model used to interpret signals. Previous analyses of GWTC-3 event candidates using SEOBNRv4PHM, IMRPhenomXPHM, and NRSur7dq4 have sometimes arrived at notably different conclusions about event properties. Subtle analysis settings and code differences can also produce notable differences. We revisit 60 publicly available events from the first three observing runs, released across GWTC-1, GWTC-2, GWTC-2.1, and GWTC-3, using a consistent analysis framework and three waveform models: two state-of-the-art time-domain models with higher-order modes (SEOBNRv5PHM and IMRPhenomTPHM) and one older frequency-domain model without higher-order modes (IMRPhenomPv2). The two state-of-the-art models agree well overall, although about 20\% of events have a visible difference in at least one one-dimensional marginal posterior. By contrast, the older model often arrives at qualitatively different conclusions for key events across the mass spectrum. We also identify a bug in published LVK GWTC-2.1 SEOBNRv4PHM results and in a later reprocessing of GW200105. At the adopted sampling rates, the model's Nyquist-frequency restriction excluded part of the allowed parameter space and artificially truncated the mass-ratio posterior away from equal mass for seven low-mass events. These results emphasize the need to use multiple state-of-the-art waveform models and carefully validated analysis settings to characterize uncertainty in gravitational-wave source properties.

The radio/X-ray/black-hole-mass Fundamental Plane (FP) is widely used to estimate masses of accreting intermediate-mass black hole (IMBH) candidates. Its inverse use is justified only when radio luminosity traces a compact jet core and X-rays arise from the same sub-Eddington hard-state accretion flow. We formulate an accretion-state validity criterion for inverse-FP masses using verified BH X-ray binary radio/X-ray tracks as the empirical reference. Using 424 detections from 21 BHXBs and candidates, we find that both standard and curved tracks are present and the luminosity at which radio/X-ray behavior changes is source dependent. A low X-ray Eddington ratio is therefore not, by itself, sufficient to establish inverse-FP validity. Within a scale-invariant jet/hot-flow framework, track steepening reduces the effective mass leverage of the inversion, while at much lower accretion rates uncertainties in the dominant emission mechanism and compact-core production can invalidate the underlying FP assumptions altogether. We apply the criterion to representative IMBH-relevant systems without fitting a new plane. NGC 4395 is low-Eddington but compact-radio deficient: its VLA-scale and VLBI jet-base luminosities imply formal inverse-FP mass underestimates of about 1.0 and more than 2.8 dex if the deficit is misread as a mass signal. Globular-cluster X-ray populations in NGC 1399 and NGC 4472 are source-identity limited: 71% and 67% lie below $L_{\rm Edd}(10\,M_\odot)$, and the high-luminosity tail alone does not identify an IMBH accretor. The radio non-detection of the $ω$ Centauri IMBH candidate constrains gas supply, accretion efficiency, and compact-core production rather than excluding the black hole. The resulting test defines which SKA/ngVLA-era radio cores can be converted into FP masses, and which instead constrain gas supply, source identity, or compact-core production.

The cosmological constant $Λ$, and its generalization as Dark Energy, makes up roughly 70% of the present-epoch universe and has been part of the dominant $Λ$CDM cosmological paradigm for more than a quarter of a century. It has been put on trial by recent measurements of galaxy clustering, Type Ia supernovae, and the cosmic microwave background. In this Perspective we evaluate the strength of the evidence for a constant $Λ$, versus Dynamical Dark Energy (DDE) in the light of recent developments. We also assess what the future holds both for our Universe and our understanding of its most mysterious component.

The ultraviolet continuum from a high-redshift accreting black hole is routinely used to infer its mass. GN-z11 offers a sharp test of whether that inference is secure. A thin-disk fit to its continuum gives an Eddington mass of $M_{\rm E}=1.12\times10^{7}\,\Msun$, an order of magnitude above the $\log(\Mbh/\Msun)=6.2\pm0.3$ inferred from broad N~\textsc{iv}. We show that this apparent tension is not imposed by the ultraviolet slope alone. We derive a closed-form criterion comparing the radius where observed photons are produced with the radius where a supercritical flow departs from efficient solution. For GN-z11, the criterion places the supercritical transition inside the ultraviolet-emitting region. Composite disks retaining the outer solution confirm this numerically, requiring $\Mdot_{\rm out}=3.1$--$4.4\,\Msun\,\mathrm{yr^{-1}}$ and placing $\Rsph=(1.2$--$1.8)\times10^{3}\,\rg$ inside $\Ruv=(2.6$--$2.9)\times10^{3}\,\rg$. Replacing the inner $\sim10^{3}\rg$ by a supercritical flow shifts the fitted slope by only $Δ\betaUV=+0.08$ to $+0.11$, comparable to uncertainties from standard spectral modelling. Reversing the radial ordering requires the effective transition to move outward by a factor $1.7$--$2.1$. Within our composite models, the measured slope constrains black-hole mass only to around $10^{6}\,\Msun$, an order of magnitude below the continuum Eddington mass. Continuum Eddington masses should therefore be treated as model-dependent bounds that assume global thin-disk efficiency. The same hierarchy appears in three additional JWST sources under the disk-dominated interpretation, where separation masses are only $4$--$6$\% of continuum Eddington masses. This distinction can substantially weaken the seed-mass and early-growth demands inferred from ultraviolet continua at cosmic dawn for the first massive black holes.

We present a three-dimensional visco-magnetoacoustic framework for analog gravity in magnetohydrodynamic (MHD) flows. While standard fluid dissipation typically breaks the Lorentzian signature of acoustic metrics, we demonstrate that evaluating wave perturbations in the eikonal limit alongside the Shakura-Sunyaev $α$-viscosity prescription preserves a well-defined effective spacetime geometry for the fast magnetoacoustic mode. The slow-magnetoacoustic mode and Alfvén mode do not admit a non-degenerate metric. To investigate analog horizon thermodynamics, we utilize astrophysical accretion disks as background media, specifically modeling numerical magnetized advective accretion flows around rotating black holes and analytical advection-dominated inflow-outflow solutions (ADIOS). Standard self-similar ADIOS models strictly enforce a constant Mach number, precluding horizon formation. We therefore introduce a magnetic field perturbation that breaks self-similarity, generates a dynamic Mach number, and enables the formation of a visco-magnetoacoustic horizon. By evaluating the spontaneous phonon emission at these horizons, we reveal that the analog Hawking temperature is highly sensitive to the magnetic field topology; the spatial orientation of the background magnetic gradients dictates whether the Hawking radiation is amplified or suppressed. Furthermore, we find that increasing the viscosity parameter leads to a monotonic increase in the analog Hawking temperature.

Recent IceCube observations point to the ultra-high energy neutrino flux exhibiting a very soft spectral nature beyond tens of TeV, which is unexpected from the standard cosmic-ray neutrino connection. As a potential explanation for this behaviour, we investigate neutrino self-interactions in an extra-dimensional $ {U(1)}_{L_μ-L_τ} $ gauge theory, where symmetry breaking leads to a tower of Kaluza-Klein gauge bosons in compactified four dimensions. These multiple gauge bosons may enhance the attenuation of astrophysical neutrinos as they propagate in the C$ν$B medium. Unlike single-mediator scenarios, for example, which arise from broken $ {U(1)}_{L_μ- L_τ} $ gauge symmetries in four dimensions, the presence of a Kaluza-Klein tower of mediators produces multiple closely spaced resonances whose interference gives rise to a rich energy-dependent behaviour of the scattering cross-section over a vast range of incident neutrino energies. Alongside $s$-channel resonances, off-resonant $t$- and $u$-channel contributions also become important. We explore the possibility that repeated resonant scattering between astrophysical neutrinos and those in the C$ν$B, mediated by these new multiple gauge bosons may increasingly attenuate the former at higher energies, thereby addressing the possibility of softening its spectral nature within the consideration of a single power law.

We present a multi-epoch, primarily X-ray study of the narrow-line Seyfert 1 (NLSy1) galaxy Mrk~766 by combining archival \textit{XMM-Newton} observations spanning nearly two decades together with a \textit{XRISM} observation obtained during its Performance Verification Phase. We analyzed the absorption and emission features in the soft X-ray spectra from RGS and EPIC-pn data, revealing a stratified multi-phase potentially ionised absorber with a column density varying at a timescale of $\sim$1 day. The covering fraction shows positive correlation with the continuum flux, however with significant intrinsic scatter, between the epochs. Mrk~766 exhibits a multi-component Fe\,K$α$ emission line as revealed by fits to the \textit{XRISM} data. A quasi-phenomenological model comprising a torus and a broad \textit{diskline} component distinguishes the two components and constrains the inner radius of the potential broad line emitting disk substructure between $R_{in}=40-60\,r_{\rm g}$. Analysis of multiple epochs of the Fe K$α$ complex across \textit{XMM-Newton} observations reveals that the broad component tracks the continuum flux closely and effectively reflects the trend with respect to the continuum. The narrow component is consistent with a distant emitter, potentially the torus. Quasi-simultaneous monitoring with \textit{Swift}-XRT and UVOT further indicates UV emission lagging the X-rays by $5.9^{+4.1}_{-5.7}$ hours, consistent with disk reprocessing of coronal X-ray fluctuations, and corresponds to a light travel distance of up to $10^{3}\,r_{\rm g}$. Overall, our analysis maps the extent of the dynamically evolving absorber, the iron line emitter, and the accretion disk in this highly accreting NLSy1 system.

The large scale clustering of gravitational wave (GW) sources offers an independent probe of cosmology and compact-binary astrophysics. Interpreting future GW clustering measurements requires understanding how merger rates depend on host-galaxy properties and environment. We develop a simulation-based framework to model the clustering of binary black hole (BBH) mergers and estimate the corresponding GW bias. We populate IllustrisTNG300 with GW sources using a machine-learning emulator trained on the GALAXYRATE data set, which provides population-synthesis calibrated merger rates linked to galaxy formation histories. This enables physically motivated mock GW catalogs with galaxy-specific BBH merger rates that depend on stellar mass, star formation rate, metallicity, and redshift. We compare this framework to a simpler phenomenological model in which the merger rate depends only on stellar mass. In both approaches, GW sources are more strongly biased than the overall galaxy population, reflecting the preference for BBH mergers to occur in more massive, more strongly clustered halos. The machine-learned model also produces stronger scale dependence in the GW bias than stellar-mass weighting alone, indicating that additional host properties affect GW clustering. We further find that metallicity correlates positively with both GW and galaxy bias. These results highlight the importance of physically motivated merger rate models for GW large scale structure analyses and provide a framework for interpreting future clustering measurements from experiments such as the Einstein Telescope and Cosmic Explorer.

Gravitational-wave searches for subsolar-mass compact binaries have left a gap between the planetary-mass regime accessible to long-duration searches and the higher masses probed by conventional compact-binary analyses. We close this gap by searching for the inspiral of ultra-compact binaries with chirp masses $10^{-2}$--$10^{-1}\,M_\odot$ in LIGO O4ab data. Using the statistically robust BinaryGFH-v2 method, a time-frequency technique to search for rapidly evolving inspirals lasting minutes to hours, we cover this previously challenging parameter space with a computational cost of only $\sim 350$ CPU-hours, and find no significant candidates. We thus set 90\% confidence upper limits on the binary merger-rate density, providing the first gravitational-wave constraints across this chirp-mass range. Assuming that primordial-black-hole binaries formed in the early Universe with no binary suppression ($f_{\rm sup}=1$), these results provide the first gravitational-wave constraints on the primordial-black-hole \dmh fraction in this mass range, yielding upper limits with $f_{\rm PBH}<1$ for equal-mass binaries with component masses $10^{-2}$--$10^{-1}\,M_\odot$ and asymmetric binaries with secondary masses $m_2=3\times10^{-4}$--$10^{-2}\,M_\odot$ for a primary mass $m_1=2.5\,M_\odot$.

Stars embedded in the disks of active galactic nuclei (AGN) can grow to hundreds of solar masses, and the same disks are expected to host a population of stellar-mass black holes. A star may therefore capture passing black holes, turning its interior into a potential factory for forming compact binary black hole systems and, ultimately, for the gravitational-wave events seen by our detectors. We test how readily this channel operates using three-dimensional hydrodynamic simulations of black hole--star encounters. As a proof of principle, we adopt a star-to-black-hole mass ratio of $34\!:\!1$, and find that the capture does not disrupt the structure of the star: it heats the star by $10$--$30\%$, and the star loses only $\lesssim1\%$ of its mass throughout. The captured black hole loses its orbital angular momentum rapidly to dynamical friction, sinking to the stellar center within a stellar dynamical time $\lesssim10^{4}\,$s. If the star already harbors a black hole at its center, the two form a bound binary whose gravitational-wave coalescence time falls below $10^{4}\,$yr. Depending on the geometry and speed of the encounter, the resulting orbit may be nearly circular or retain substantial eccentricity. Our calculations confirm that black hole--star encounters in AGN disks are indeed a viable channel for assembling stellar-mass black hole binaries and supplying sources for gravitational-wave detectors.

In this work we establish general constitutive equations for a relativistic plasma composed of an arbitrary number of classical, charged, and chemically non-reacting species in the presence of a background electromagnetic field. The intensity of the electromagnetic force acting on either species is assumed to be comparable with the gradients of the state variables and is thus considered as a first-order driving force for dissipation within the gradient expansion. Assuming that the species share the same temperature and velocity flow when in equilibrium, we determine the entropy production and establish constitutive equations in terms of first-order frame invariant quantities. Conditions on the corresponding transport coefficients are given in order to comply with Onsager's reciprocal relations and the second law of thermodynamics. In particular, we discuss cross effects for the binary mixture which may be relevant for astrophysical and cosmological applications. Finally, we derive a complete set of evolution equations governing the dynamics of a binary mixture propagating on an arbitrary spacetime background.

Context. Globular clusters (GCs) are fragments of the chemical history and evolution of the Milky Way. Their notable intra-cluster abundance variations have been used to characterise chemical trends and stellar populations. However, the origins of these populations and GCs in general are still largely unknown. Aims. Our aim is to explore the origin of the stellar populations belonging to 14 GCs making use of their abundances derived from the Gaia-ESO Survey (GES) spectra. In particular, we use the Mg-Al anti-correlation to separate the populations. Methods. The GES spectra were revisited to derive additional homogenous chemical abundances for Mg and Al. A Gaussian Mixture Model (GMM) was applied to the [Mg/Al] distribution for each cluster to define a threshold between the two main populations. Results. Eight clusters imply the presence of Mg-Al anti-correlations and dual populations. We find there is a regime transition near [Fe/H] approximately -1.4, where, the Mg-Al anti-correlation disappears at higher metallicity. Several clusters, regardless of their Mg-Al anti- correlation behaviour, exhibit metallicity spreads (e.g., NGC 2808, NGC 6752) and s-process variations (e.g., NGC 7078, NGC 1851), suggesting complex enrichment histories. Conclusions. The presence and spread of the Mg-Al anti-correlation depend on cluster mass and metallicity, with metal poor and/or massive clusters exhibiting the clearest anti-correlation. This work demonstrates additional abundances that can be retrieved from Gaia-ESO Survey spectra and provides further insight into the chemistry of multiple stellar populations within GCs. These findings highlight the chemical complexity of GCs and the need for large, homogeneous datasets to trace their formation histories.

We forecast constraints on the parameters of the coupled dark energy model using DESI and Euclid data with the one-loop galaxy power spectrum. We investigate the distinguishability of our model from the zero-coupling scenario at the $1σ$ level and explore how the parameter constraints depend on the fiducial coupling $β_{\rm fid}$, the fixed potential slope parameter $ν$, the maximum wavenumber $k_{\max}$, and different prior choices. We find that the inclusion of mildly non-linear scales improves the constraints on the coupling by roughly a factor of five. Then, we address the question of which is the minimum value of $β$ that can be distinguished from zero at $1σ$. We find that for our reference case $k_{\rm max}=0.2 h/$Mpc, $β=0.15$ lies roughly $1σ$ above zero. In the most optimistic case with $k_{\rm max}=0.3 h/$Mpc and including a Planck prior on $Ω_{m0}$, this value can be reduced to $0.05$. These values are substantially larger than the current constraints on $β$, but the latter have been obtained assuming $β$ to be constant from at least the decoupling epoch to today, while we only employ late-time data. We conclude therefore that only models that allow for time-varying couplings can be detected with late-time datasets.

We apply a newly proposed bound state method to revisit the computation of gravitational and electromagnetic quasinormal modes (QNMs) for Reissner-Nordström black holes (RNBHs). It is found that the method yields QNM frequencies of high accuracy for low-lying modes with overtones $n\leq3$. The accuracy degrades for higher overtones, however, a homotopy deformation to the potential enables us to compute more high overtone modes reliably. It is known that the continued fraction method for computing QNMs of extremal RNBHs differs significantly from that used for nonextremal ones. In contrast, the bound state method demonstrates advantage of simplicity: one only need to account for the definition of tortoise coordinate in the extremal case, and the rest of the computational procedure remains exactly the same as that for the nonextremal case.

For the $k$-self-similar naked singularity solutions of Einstein--scalar field equations constructed by Christodoulou, we construct a family of interior gravitational perturbations leading to trapped surface formation that remain large but finite at the threshold Hölder regularity, while converging to zero in all regularities below the threshold. This extends the interior spherically symmetric result below the threshold in \cite{Li25} to non-spherically symmetric setting. We further construct a genuinely localized family of perturbations whose angular support shrinks to a single point, with smooth convergence to zero away from that point. This exploits the additional angular freedom available beyond spherical symmetry. Moreover, when measured instead in Sobolev regularity, such genuinely localized family of perturbations still converges to zero in all regularities below the same threshold.

We investigate whether the recently reported high-energy nuclear recoil event by the LUX-ZEPLIN (LZ) collaboration can be interpreted as a signal of freeze-in dark matter (DM) in a non-standard cosmological history. We consider DM production from the Standard Model thermal bath during an era of modified expansion preceding Big Bang nucleosynthesis. Requiring consistency with the observed DM abundance, cosmological and collider constraints, we identify viable parameter regions that can accommodate the LZ event. Our results demonstrate the potential of direct detection experiments to probe non-standard cosmological histories and, more broadly, alternative theories of gravity.

In this manuscript, we develop topologically charged static and spherically symmetric wormhole (WH) configurations in Einstein gravity by employing the Einasto dark matter profile to construct an analytic shape function. The resulting spacetime satisfies the fundamental traversability conditions, while the matter sector exhibits localized violations of energy conditions near the throat, indicating confined exotic matter content. We analyze equilibrium through the TOV framework, revealing a non-equilibrium force competition in which hydrostatic and anisotropic contributions self-organize to sustain static configurations within specific parameter domains. Stability is examined via the anisotropy parameter and further extended to linearized radial perturbations of the surrounding shell, uncovering parameter-dependent stability windows that characterize the system's nonlinear response. A detailed dynamical analysis shows that the monopole charge parameter significantly restructures the effective potential landscape, modifying curvature and inducing qualitative transitions in particle motion. The coupled interplay between monopole charge, angular momentum, and particle energy governs a nonlinear transition from tightly wound quasi-bound states to unbounded scattering trajectories. This transition reflects an emergent dynamical phase structure in the WH spacetime. Furthermore, the complexity factor exhibits strong localization near the throat and vanishes asymptotically, indicating that structural complexity is confined to the inner nonlinear regime, whereas a volume integral quantifier is employed to estimate the total exotic matter content required to sustain the configuration.

FOSSIL (FTS fOr CMB Spectral diStortIon expLoration) is a proposed ESA M8 mission tailored to measure the spectral distortions of the Cosmic Microwave Background (CMB) with a sensitivity three orders of magnitude beyond the COBE/FIRAS legacy measurement. Achieving this sensitivity demands an extraordinarily challenging cryogenic architecture: the scientific instrument must be maintained at the lowest achievable temperature (4.5 K), while the detector focal plane assembly operates at 50 mK. This paper presents the thermal architecture of the FOSSIL payload and instrument, from the spacecraft service module at 293 K down to the sub-kelvin detector stage. The thermal design draws on heritage from the Planck and ARIEL missions and relies on a staged passive cooling chain comprising a multi-layer insulation blanket, three V-groove radiators (operating at approximately 130 K, 90 K, and 50 K), and a 25 K actively cooled shield fed by an ESA-provided 4 K cryocooler chain. Sub-kelvin temperatures are achieved via a multi-stage adiabatic demagnetisation refrigerator (ADR) system developed for NewAthena/X-IFU, providing continuous cooling at 1.8 K and 350 mK, and 50 mK with an 80% duty cycle. The Focal Plane Assembly (FPA), which houses four Kinetic Inductance Detector (KID) arrays at 50 mK, is thermally isolated from the 4.5~K bench via a carbon-fibre reinforced polymer (CFRP) hexapod structure with staged heat interception. We present the steady-state thermal budget across all stages, demonstrating comfortable margins at every temperature level, and discuss key thermal design drivers on the Blackbody Internal Reference (BBIR).

In this paper, we derive a new class of analytic black hole solutions within the framework of Hoyle--Narlikar gravity, where the black hole is surrounded by an electric charge acting as the source. After obtaining the corresponding metric, we investigate the thermodynamic properties of the resulting black hole solutions, including the Hawking temperature and entropy. Subsequently, we compute the specific heat and Gibbs free energy in order to analyse the thermodynamic stability of the black hole in Hoyle--Narlikar gravity. An interesting result of our analysis is that the black hole becomes thermodynamically stable for large values of the horizon radius within this framework.

Identifying a shock from a moving solar brightness front requires relating the observed emission to changes in the plasma state. The Riemann framework describes shocks, rarefactions and other waves as parts of a pattern connecting different states. We present the Riemann Map Operator (RMO), an inverse diagnostic framework implemented as a web application to help observers identify magnetohydrodynamic (MHD) shock families from incomplete measurements. RMO tests candidate local connections against conservation laws, entropy requirements and characteristic conditions, while retaining observational uncertainties and model assumptions. We demonstrate the approach with five solar examples and a complementary in-situ comparison. An event combining extreme-ultraviolet (EUV) imaging and spectroscopy admits both fast and slow shock connections within the tested conditions. For an EUV-radio event, we derive an upstream-flow constraint that would exclude ordinary slow shocks within the adopted model, provided both diagnostics sample the same front. A flare-loop benchmark reproduces the published spectral profiles and recovers slow-shock ordering in the isothermal limit. The Solar Orbiter comparison recovers an Alfvenic relation between measured velocity and magnetic-field changes. Together, the examples show how RMO connects observations to admissible MHD interpretations and identifies the additional measurements needed to distinguish them.

We present the first results of two Spanish TAC programs using the WEAVE LIFU that have been obtaining data since October 2023. The first one is obtaining high-resolution multi-epoch spectroscopy of the nine stellar clusters that contain a minimum number of 3 O-type stars (and many more B stars) within the field of view of the WEAVE LIFU and that are accessible from La Palma with the purpose of studying the spectroscopic multiplicity of the OB stars. The second program is obtaining single-shot low-resolution spectroscopy of as many stellar clusters with OB stars as possible with the purpose of characterizing the OB-star population in those regions that are too crowded to be properly studied with the guaranteed time SCIP program. That program has already observed 54 clusters in a wide range of environments, from H II regions to gas-free clusters, and distances, from clusters in the solar neighbourhood to (nearby) extragalactic systems. We also present the tools developed for this program, which include a data cube processing and visualisation utility and UNWIND, a code that fits and subtracts the ISM signature in the observed spectra.

We uncover the quantum geometric structure of graviton Wigner functions and stress-energy tensors via quantum field theory in curved spacetime, revealing the phase-space geometry of spacetime and relativistic quantum-state manifolds. We show that the polarization mode expansion of the underlying graviton field operator is fully captured by quantum geometry, as encoded in the spinor-helicity formalism, thereby establishing the geometric framework from the outset. Applying this to a weak gravitational background, we demonstrate the roles of the quantum metric and connection in describing graviton transport beyond the chiral vortical effect. We also clarify the role of the quantum metric in normalizing the perfect-fluid graviton stress tensor within this approach. This work paves the path for explorations of multistate Hilbert-space geometry in high-energy and gravitational physics. In addition, this framework naturally encompasses lower-spin excitations, allowing for a unified quantum geometric treatment of bosonic and fermionic many-body systems in condensed matter and particle physics.

Black holes in asymptotically anti-de Sitter (AdS) spacetime whose extended-phase-space thermodynamics reproduces that of a van der Waals fluid exactly form a distinctive class, the van der Waals black holes (VBHs). The timelike geodesics of VBHs are solved, their periodic orbits are catalogued, and the gravitational waves emitted along them are computed, so as to trace how the thermodynamic pressure P and the molecular-volume parameter b shape the orbital dynamics. Integrating the geodesic equations numerically yields the energy, angular momentum and radius of the innermost stable circular orbit (ISCO) and maps out the family of bound orbits. To delimit the parameter window in which the effective potential stays well behaved, we introduce a critical angular momentum Ls. Within the orbit-classification scheme, we compute the energies corresponding to rational q for several parameter sets and display the resulting orbits. The waveforms produced along these orbits are generated within the "Kludge" formalism, which makes it possible to quantify the imprint of b on the wave amplitude and on the period. Both the orbital architecture and the gravitational-wave signature turn out to be strongly sensitive to b. The dynamical behavior of VBHs and Schwarzschild-AdS (SAdS) black holes differs markedly, providing a new perspective for probing black hole thermodynamics via gravitational wave observations.

Lunar laser ranging (LLR) provides a sensitive test of the equivalence principle (EP) through a possible synodic perturbation in the Earth--Moon distance. In this work, we develop an independent LLR data-reduction framework and investigate this signature in the post-fit O--C residuals obtained with the EPM21, INPOP21, and DE430 planetary and lunar ephemerides. The same observation model and parameter-adjustment procedure are applied to the three ephemerides, allowing their residual-based results to be compared within a consistent framework. To account for the nonuniform distribution of LLR observations over the synodic angle $D$, the residuals are fitted with cosine and sine harmonics up to the third order. Signal-injection tests are also performed to quantify the attenuation of a possible $\cos D$ signal during the parameter adjustment. After correcting for signal attenuation, solar radiation pressure, and the thermal response of the lunar retroreflectors, the preferred EPM21, INPOP21, and DE430 solutions yield $Δ(m_g/m_i)_{\rm EM}=(0.900\pm3.183)\times10^{-14}$, $(0.496\pm3.624)\times10^{-14}$, and $(2.669\pm7.048)\times10^{-14}$, respectively. The EPM21 and INPOP21 results are based on observations from 1970 to 2024, whereas the DE430 solution was obtained using observations up to 2016. All three estimates are consistent with zero within their uncertainties. This work demonstrates that residual-based analysis can provide an independent and transparent EP test, while also offering guidance for assessing sampling and ephemeris-dependent effects in future high-precision LLR studies.

We present a deterministic machine-learning method for denoising candidate signals from all-sky searches for continuous gravitational waves from rotating neutron stars. Using the time-domain $F$-statistic search pipeline, we post-process the resulting $F(f,\dot{f})$ candidate patterns with iterative $α$-(de)blending, a deterministic diffusion-type generative model implemented as a U-Net, well suited to the non-Gaussian, correlated noise of the $F$-statistic output. Two models, trained on data based on simulated 6- and 12-day time-domain segments with software-injected signals similar to the LVK hardware injections, are tested by comparing deblended images to the library of expected signal patterns via the Structural Similarity Index Measure. In general the method recovers the correct sky-position-dependent pattern for signal-to-noise ratios $ρ\gtrsim 4$, moderately below typical all-sky detection thresholds, with recovery depending strongly on pattern morphology, demonstrating that deterministic diffusion-based denoising may serve as a consistency/veto tool ahead of the semi-coherent coincidence stage. We also discuss limitations and possible improvements of the method.

We explore the cosmological dynamics and formation of structures within the Tsallis-Cirto modified gravity framework with a model parameter $β$ that is inspired by nonextensive statistics. By applying the gravity-thermodynamics conjecture, we modify the Friedmann equations and the evolution of the Hubble parameter. The impact of the Tsallis-Cirto entropy on different cosmographic parameters and the development of the linear matter overdensities have been analyzed by formulating perturbed field equations using the spherical collapse approach in a flat Friedmann-Lemaître-Robertson-Walker background. We present a novel and established diagnostic approach to distinguish among distinct cosmological models compared to flat and non-flat $Λ$CDM scenarios, discovering that the Tsallis-Cirto entropy-inspired modified cosmology ($β\ne 1$) successfully withstands all tests, contradicting both the flat and non-flat $Λ$CDM models. This model also meets the conditions for the Universe to reach thermodynamic equilibrium in the far future. Additionally, we investigate the halo mass function and cluster number counts within this modified gravity framework. All findings are compared to the fiducial $Λ$CDM profile, indicating that the added entropic correction affects the history of expansion, the growth rate of structures, and the dark matter halo abundance. We find that the more massive structures are less abundant and develop at later times, which aligns with the hierarchical model of formation of large-scale structures.

Halo assembly bias describes the dependence of halo clustering on properties beyond halo mass and provides a means of connecting halo assembly history to the large-scale matter distribution. We develop and test a relative clustering methodology for measuring differences in the large-scale bias of two halo populations using the ratio of their cluster--galaxy pair counts. In the linear regime, this ratio can be related to the difference in cluster bias, while the galaxy autocorrelation independently constrains the galaxy bias. We apply the method to a volume-limited X-ray selected galaxy clusters from the eROSITA eRASS1 catalogue, dividing them into high- and low-temperature subsamples with matched distributions in X-ray luminosity, redshift, and exposure. We first verify the validity of the linear approximation, finding $ξ(r)\ll1.0$ over $4$--$10\,h^{-1}\mathrm{Mpc}$. The joint clustering analysis yields $b_{\rm g}=1.610^{+0.018}_{-0.019}$ and a bias difference $Δb=0.8\pm1.1$. DES Y3 weak gravitational lensing measurements provide an independent test of the halo masses, yielding $\log(M_{200\mathrm{m}}/h^{-1}M_\odot)=14.33^{+0.07}_{-0.07}$ and $14.33^{+0.10}_{-0.09}$ for the high- and low-temperature samples, respectively, with a mass difference of $Δ\log M_{200\mathrm{m}}=-0.01\pm0.12$. We therefore find no statistically significant evidence for assembly bias at the $1σ$ level, but our measurement is consistent with the bias difference measured for analogous halo populations in the FLAMINGO simulation. The result therefore does not rule out assembly bias at the level expected from the simulation, while demonstrating that the relative clustering methodology provides a framework for future assembly-bias measurements with larger X-ray cluster and galaxy surveys.

Observations of the Large Scale Structure sit on the light cone. Cosmological statistics must then be projected from 3D to angular space. For the bispectrum and higher orders, the polyspectra crucially depend on the angles between Fourier modes in addition to their moduli. This leads to delicate cancellations in the squeezed limit and make Limber's approximation fail. This is problematic for predicting the Super-Sample Covariance (SSC), with no known full-sky computation of the angle-dependent terms. Here we develop a new analytical approach to the angular bispectrum, writing Fourier vectors as derivatives that we project on the sphere and tackle with spin-weighted spherical harmonics. We find bispectrum equations free of cancellations and amenable to beyond-Limber numerical methods. Limber's approximation becomes well-behaved and works in the squeezed limit and any other configuration. We then derive the consequence for the covariance between 1-point (e.g. cluster or peak counts) and 2-point statistics (e.g. galaxy clustering, weak lensing or 21cm), which relates to the squeezed bispectrum. The results split into intra-survey covariance and SSC, where we find significant differences with past literature. Our results are in line with the 2+1 geometry of the light cone, having e.g. temporal derivatives in addition to spatial ones. We also find that the tidal term is degenerate with the simpler growth-only term. The total SSC contains a wealth of new terms including wide-angle effects and responses to the background velocity. However we show numerically that it is well approximated by a single term with a density-only response. This ultimately makes our results simpler than past literature. This full-sky result should be more representative than the past flat-sky limit for large survey areas, thus being the one relevant for coming observations with e.g. Euclid, LSST and Roman.

Dark-matter self-interactions can drive subhalos through gravothermal core collapse, but conventional lensing observables vary continuously with uncertain inner profiles. We show that gravitational-wave image multiplicity provides a discrete alternative. Near a macro-critical minimum, increasing subhalo central convergence opens two caustics in sequence, converting one detectable waveform into three and then five stationary images. Using weighted SASHIMI subhalos and an ET+CE selection model, we forecast almost 2 systems with at least three copies detectable in five years if dense remnants survive, compared with the prediction of 0.36 in 5yrs for standard CDM. The much more diagnostic second case is rare: the corresponding 5yrs yield with at least four detectable copies is 0.149, versus $1.55\times10^{-6}$ for CDM. According to the Poisson statistics, this gives only a 13.8\% probability of at least one event detected in 5 years. Thus this is a rare-event, exposure-limited test rather than a guaranteed detection: a resolved second pair would be difficult to reconcile with the fiducial single-halo CDM population, whereas a five-year null result would not by itself constrain remnant survival strongly.

We study slowly rotating, nonlinearly charged black holes in arbitrary spacetime dimension $d \geq 4$, using a generalized Lense--Thirring ansatz. Working to linear order in the rotation parameters, we solve the Einstein--nonlinear electrodynamics (NLE) field equations for restricted theories ${\cal L}({\cal S})$ depending on single invariant. We provide an extension of the four-dimensional "No Go theorem" proven in ArXiv:2203.01919 and demonstrate the uniqueness of Maxwell theory in $d \geq 4$. Similar to four dimensions, it is proven that the standard rotation generating trick fails to produce full rotating solutions for any nontrivial NLE in higher dimensions. We identify a family of NLE theories admitting slowly rotating solutions with vector potentials encoding rotation effect in Maxwell-like manner; in four dimensions this reproduces the previously found RegMax theory, while in $d\geq4$ we obtain the corresponding electric potential in closed form and the associated Lagrangian implicitly as a function of the radial coordinate.

We present a new approach to constraining the delay time distribution (DTD) governing r- and s-process element production. Rather than using metallicity as a proxy for time, we work directly with stellar ages at solar metallicity, where chemical equilibrium models are most applicable. We analyze two complementary samples: more than $10,000$ red clump (RC) thin disk stars at solar metallicity from LAMOST, and solar twins with high-precision abundances; both have well-determined stellar ages $τ$. We fit chemical-equilibrium models with power-law DTDs to the evolution of [Eu/Mg]$(τ)$ and [Ba/Mg]$(τ)$ in each sample, taking europium as a canonical r-process element and barium as an s-process element, each referenced to magnesium as a near-instantaneous, CCSNe-dominated element. To validate the method, we first apply it to [Fe/Mg]$(τ)$ and recover that $\approx\!70\%$ of Fe comes from a delayed channel with a declining DTD $\propto t^{-1.2}$, in agreement with well-established results for iron production in thermonuclear supernovae. It also shows that about $60\%$ of Eu and Ba forming today originate from strongly-delayed processes, and the inferred delayed-production rates of both Eu and Ba rise with delay time, opposite to that of thermonuclear supernovae. These results provide direct temporal constraints on r- and s-process enrichment, complementary to metallicity-based chemical-evolution studies. Because we model abundance referenced to Mg, and the same trend is recovered across stars of different guiding-center radii, this appears unaffected by radial migration and chemical evolution varying with Galactocentric radius. If confirmed, it points to delayed neutron-capture production beyond NS--NS mergers.

Binary stars are the progenitors of exotic objects, such as supernovae and gravitational wave sources. However, some systems merge earlier when the primary expands into a giant star at the end of its lifetime on the main sequence and causes its companion to spiral in. TYC 4144-329-2 is a candidate post-merger first ascent giant star with a circumstellar disk seen in the infrared (IR). We detect weak X-ray and far ultraviolet (FUV) emission from TYC 4144-329-2, which we interpret as a signature of weak coronal activity. We also find variability in Ha line profiles and in the optical light curves and we suggest an intermediate geometry where the line-of-sight passes through the upper layers of a flared disk with a time variable column density. While the X-rays must be coronal, the Ha line profiles point to ongoing and variable accretion. We suggest that TYC 4144-329-2's merger was more recent than other stars in this class and that it did not yet have time to develop a deep convection zone. If this scenario holds, TYC 4144-329-2 would be a unique probe of the earliest stages in the evolution of post-merger systems.

The paper discusses the fraction of extragalactic cosmic iron nuclei with energies of 10^19 eV that reached Earth from distant sources without fragmentation in intergalactic space. It is shown that the fraction of particles with energies E = 10^19 eV is significant (about 10%), even if the redshifts of the sources reach 0.2, and not only for nearby sources.

The kinetic Sunyaev-Zeldovich (kSZ) effect is a sensitive probe of the free-electron momentum field, tapping into spatial correlations relevant for both large-scale cosmology and smaller-scale astrophysics. Precision kSZ measurements will provide a robust channel for constraining inflationary physics and equivalence principle violation, while simultaneously characterizing the small-scale gas distribution and associated baryonic feedback. However, robust kSZ measurement for cosmological applications is challenging, as the primary source of large-scale information is convolved with unknown small-scale dynamics in a degenerate combination. In this paper we address this challenge using the Effective Field Theory (EFT) framework that provides a principled way of treating the effects of small-scale dynamics on large scale-structure. Extracting the kSZ signal can be reformulated as a measurement of the one-loop galaxy-galaxy-electron momentum bispectrum. We compute here a pure dark matter version of this bispectrum in EFT and compare it to simulations. We derive the EFT counterterms that capture small-scale backreaction and find that two of them acquire a transverse direction. One transverse counterterm gets generated entirely by the vorticity of the velocity field. Our computation of the bispectrum dipole agrees with simulation data to $5\%$ up to $k_{\rm max}=0.23~iM$, which is significantly larger than the reach of tree level theory, which breaks down beyond k_max = 0.23 h/Mpc and k_max=0.07 h/Mpc (at redshift $z=0.5$). For the bispectrum computed with the projected momentum field, relevant to measuring the kSZ field, we find a similar $k_{\rm max}$, and make the first detection of the transverse counterterms from simulations. Our results suggest that one-loop bispectrum computations can play an important role in advancing kSZ cosmology.

We study a curvature-coupled dark energy model that can modify cosmological evolution both before recombination and during the late-time accelerated era. The model belongs to the class of scalar-tensor theories, in which a quintessence field is non-minimally coupled to the Ricci scalar. We specify the model through a shifted quartic coupling, $f(\varphi)=α(\varphi^2-\varphi_{\rm today}^2)^2$, and an inverse power-law potential, $V(\varphi)=Λ\varphi^{-σ}$, where $\varphi_{\rm today}$ is a constant fixed by requiring the effective Planck mass to recover its present-day normalization. We implement the model in a modified version of $\mathtt{hi\_class}$ and compute its background and linear cosmological predictions. This specific form of non-minimal coupling allows an effective crossing of the phantom divide at late times while naturally suppressing deviations from standard gravity today and satisfying local gravity constraints. At the same time, the scalar field can modify the expansion history before recombination, shifting the acoustic scale in the direction required to alleviate the $H_0$ tension, while remaining dynamically relevant as dark energy at low redshift. For the parameter choices studied here, we find tens-of-percent deviations from $Λ$CDM in the expansion history, matter clustering, and metric-potential spectra. The modified evolution of the gravitational potentials leaves characteristic signatures in relativistic observables, with weak-lensing power suppressed by $\mathcal{O}(20$-$40\%)$ at low multipoles and order-unity changes in the late Integrated Sachs-Wolfe signal. Together, these results reveal a broad, scale-dependent phenomenology that motivates both a full parameter-space analysis and an extension of the present linear treatment to a dedicated non-linear $N$-body implementation.

The neutron skin, defined by the difference between neutron and proton root-mean-square radii, is a characteristic manifestation of isospin asymmetry and an important probe of the isovector nuclear interaction. This focused review examines how neutron skins influence particle emission and collective dynamics in heavy-ion collisions, from the Fermi-energy regime to ultra-relativistic energies. By modifying the initial neutron and proton density profiles, the neutron skin affects the isospin composition and geometry of the participant region, pre-equilibrium emission, particle production, fragment formation, and collective flow. We review neutron-to-proton and $\rm{t}/^3\rm{He}$ yield ratios, light clusters, pion ratios, bremsstrahlung photons, isoscaling and fragment momentum distributions, and neutron-proton differential flow and momentum observables, emphasizing their interplay with the symmetry energy and transport dynamics. At high energies, neutron skins also modify the initial geometry, eccentricities, multiplicities, and anisotropic flows in isobar and heavy-nucleus collisions. We discuss the challenge of disentangling these effects from deformation, surface diffuseness, shell structure, clustering, and model dependence. Broader connections to parity-violating electron scattering, dipole responses, coherent elastic neutrino-nucleus scattering, SRC-induced proton skins in momentum space, and neutron-star observables are also explored. Finally, we highlight opportunities from radioactive beams, improved collision experiments, microscopic many-body and transport calculations, and Bayesian inference. Combining multiple reaction systems and observables with complementary nuclear-structure and astrophysical information will be essential for quantitatively constraining neutron skins and the density dependence of the symmetry energy.

We report the first detection of the neutral methyl radical CH$_3$ in a highly evolved C-rich object, SMP LMC 011, using JWST MIRI/MRS. CH$_3$ is well fitted by an excitation temperature of $T_{\rm ex} \simeq$~190\,K and a column density of $N_{\rm tot}(\mathrm{CH_3}) \simeq 5.6\times10^{17}$ cm$^{-2}$. We also report the non-detection of ethane (C$_2$H$_6$), consistent with CH$_3$ reacting preferentially with unsaturated radicals rather than self-recombination, which may be inefficient in SMP LMC 011. Combined with the exceptionally large benzene column density of SMP LMC 011, this supports a methyl-addition route from benzene towards alkyl-substituted aromatics. We propose that the high CH$_3$ abundance is driven by the erosion of hydrogenated amorphous carbon (HAC) dust grains in the dense warm torus by UV photons from the central star and/or by shocks, which release CH$_3$ directly into the gas phase. These results establish CH$_3$ as a key reactive intermediate in the formation of complex hydrocarbons in C-rich circumstellar environments. Chemical models that incorporate methyl-addition reactions are needed for a better understanding of PAH formation pathways in evolved stars. The potential detection of alkyl-substituted aromatics such as toluene (C$_7$H$_8$) and ethylbenzene (C$_8$H$_{10}$) in C-rich evolved stars would provide direct confirmation of CH$_3$-driven aromatic growth in circumstellar environments.

Next-generation gravitational-wave observatories like Cosmic Explorer will provide an unprecedented opportunity to probe the cosmological stochastic gravitational-wave background, providing insight into the early Universe that is inaccessible with other means. We investigate the prospects for measuring a Gaussian cosmological gravitational-wave background in Cosmic Explorer in the presence of a distinguishable, non-Gaussian foreground of binary black hole mergers, using a Bayesian search that simultaneously fits binary black hole and background parameters. We successfully recover a Gaussian cosmological background with a signal-to-noise ratio of 5 in one day of simulated data in Cosmic Explorer without contamination from unresolved signals or data exclusion. However, when introducing waveform systematic errors to our analysis, our inference on the stochastic background parameters becomes biased and can falsely infer the presence of a background when there is none. The bias is driven by signal mismatches at low frequencies (<20 Hz), especially when the recovery waveform is missing physical effects included in the injected signals. While we propose a metric for diagnosing if measurements will be strongly biased by waveform systematics, our results suggest improvements in waveform modeling or other mitigation techniques will be required to accurately recover a primordial stochastic gravitational-wave background in next-generation detectors.

We present the first application of neural-network-based baryon acoustic oscillation (BAO) reconstruction to real galaxy survey data, restoring the acoustic signature damped by nonlinear structure growth. {\texttt{NERV}} ({\bf N}eural-network {\bf E}nhanced {\bf R}econstruction of the Uni{\bf V}erse) augments standard reconstruction with a convolutional neural network trained on cubic $N$-body simulations, and explicitly accounts for realistic observational effects including the curved-sky geometry, the redshift-dependent selection function, and finite survey boundaries, by tessellating the survey volume into local patches. We validate the method on the \textsc{MultiDark-Patchy} mock catalogs, recovering unbiased BAO dilation parameters. Applied to the BOSS DR12 galaxy sample, NERV improves the precision of the BAO distance measurements significantly. These results establish neural reconstruction as a practical component of BAO analyses for ongoing surveys such as DESI, with the potential to substantially tighten constraints on the cosmic expansion history and the nature of dark energy.

We confront a class of dark-energy equations of state emerging from group field theory (GFT) quantum gravity with DESI Data Release 2 baryon acoustic oscillations and Pantheon+ type-Ia supernovae. We introduce sampling parametrisations that replace microscopic initial-condition parameters by combinations more directly measured by background probes. The GFT solutions separate into logarithmic, power-law and oscillatory branches, determined by the microscopic interaction parameter $m$. The logarithmic branch is constrained to lie extremely close to a cosmological constant, while the power-law branch permits a small phantom deviation. Without perturbative-theory priors, oscillatory solutions can reproduce the mild preference of the distance data for a dip in $w(z)$ near $z\simeq0.5$--$1$. Profile-likelihood constraints favour $m\sim-2$ and $m\sim-5$ from BAO and supernovae, shifting towards $m\sim-3.5$ when CMB information is included. Conservative perturbative priors strongly suppress these deviations from $Λ$CDM. The quantum-gravity scale $z_q$, related to the average number of quantum gravity atoms, remains unconstrained, although its role in the time evolution makes higher-redshift observations a promising route to probing it. We further find that strong projection effects highlight the importance of performing likelihood profiling alongside our marginal posterior constraints. Our results provide a first direct test of GFT-motivated dynamical dark energy and demonstrate the potential for cosmological observations to inform quantum-gravity model building.

We find the necessary and sufficient conditions for the gravitational path integral (GPI) to admit an interpretation in terms of (i) an inner product in the Hilbert space of open/closed universes or (ii) a statistical average of boundary observables. The conditions needed for (i) are necessary but not sufficient for (ii). Requiring the stronger condition (ii) places additional constraints on all wormhole amplitudes, giving a concrete diagnostic of whether the GPI computes an average or merely a pseudo-average. To exemplify the difference between the two constraints, we present a class of gravitational toy models that satisfy (i) but fail the stronger condition (ii); to emphasize the power of the stronger condition, we initiate a bootstrap study of statistical ensembles based on condition (ii), finding that the resulting constraints substantially sharpen the bounds on the moments of the ensemble obtained by imposing (i). More broadly, both of these positivity conditions impose nontrivial constraints on wormhole amplitudes and thereby restrict which gravitational effective field theories can be embedded in consistent theories of quantum gravity.

We introduce a model of dark matter charged under a new confining gauge group akin to the Standard Model's quantum chromodynamics (QCD). The fermionic sector contains both heavy and light quark-like fermions charged under the dark QCD. As the Universe cools, the gauge force confines and dark radiation transitions from free particles to low-mass dark pions. This form of dark radiation behaves as an imperfect fluid with nonzero viscosity. If the confinement occurs at the MeV scale, the model can alleviate the tension between early- and late-time measurements of the Hubble parameter, with a Bayesian posterior of $H_0 = 71.60\pm 0.63$ km/s/Mpc and a $3σ$ residual tension across Planck, ACT, BBN, Pantheon+, DESI DR2, and SH0ES. In the late Universe, the heavy dark matter can evade existing constraints on dark matter self-interactions while containing a rich phenomenology that may be accessible in future observations.

The "Little red dots" (LRDs) are a population of accreting supermassive black holes (SMBHs) in the early Universe which often exhibit undermassive or even undetectable stellar hosts. Their early emergence, high space density, and extremely large black-hole-to-stellar mass ratios pose a serious challenge to conventional seeding scenarios that rely on baryon for both the formation and growth of black holes. Here we demonstrate that the above anomalies can be naturally resolved if dark matter is self-interacting. We apply a fully relativistic, non-equilibrium halo-evolution model, first developed in our earlier work, to trace the complete gravothermal evolution of self-interacting dark matter (SIDM) halos, from the initial collapse into BH seeds to the subsequent accretion of dark matter. We find that in highly concentrated halos assembled before reionization, gravothermal collapse efficiently produces stellar-mass black-hole seeds within a few hundred million years. Remarkably, and contrary to standard expectations for dark-matter accretion, heat conduction in SIDM then sustains a prolonged super-Bondi inflow that drives these seeds to supermassive scale by the LRD epoch, without baryonic assistance. The halo conditions required for completing these processes, together with the probability of avoiding major mergers that disrupt gravothermal evolution, result in an SMBH population consistent with the observed abundance and redshift distribution of LRDs. Our findings establish a pathway in which SMBHs are seeded and assembled primarily from dark matter, well before substantial galaxies form around them, thereby offering both a compelling physical explanation for LRDs and a new observational probe of dark-matter microphysics.

Kehle and Unger have recently shown that, under suitable conditions, a massless charged classical scalar field can produce an extremal black hole in finite time. This result appears to violate the third law of thermodynamics (TLT), since extremal black holes have zero temperature. However, because the TLT is intrinsically a quantum principle, a self-consistent analysis of this problem must account for quantum effects. In particular, we point out that the electric field of the system may spontaneously polarize the vacuum, thereby creating massless charged pairs that will discharge the black hole. This mechanism guarantees the {\it quantum} validity of the TLT in the regime $|eQ|\geq\hbar/2$ (here ${Q,e}$ denote, respectively, the electric charges of the black hole and the scalar field). This leaves open the question of what happens for lower values of $|eQ|/\hbar$. Recently, Schneider derived a classical bound showing that extremal black holes cannot be formed à la Kehle-Unger in the small-charge regime $|eQ|\leq\hbar/3$. His analysis leaves open the possibility that a stronger bound may exist. We conjecture here that this bound can indeed be strengthened such that classical considerations would prevent extremal (zero-temperature) charged Reissner-Nordström black holes from forming dynamically in the regime $|eQ|<\hbar/2$, where quantum polarization effects are too weak to preserve the validity of the TLT.

We present simulations of the microwave sky at 2.5 arcsecond resolution over 100 square degrees, generated from existing full-sky, lower-resolution simulations, and use them to demonstrate extragalactic foreground removal for a CMB-HD survey. Our cleaning method detects and removes the cosmic infrared background and radio galaxies, yielding source catalogs that are 95% complete down to flux limits of 0.008 and 0.04 mJy at 90 and 148 GHz, respectively. It also identifies and removes galaxy clusters via the thermal Sunyaev-Zel'dovich effect, producing a cluster sample that is nearly complete above $M_{500c} = 5 \times 10^{13} M_\odot$. After cleaning, the residual foreground-plus-noise power spectrum of the coadded 90 and 148 GHz temperature map is 50% higher than previous idealized estimates, increasing cosmological parameter uncertainties for an 11-parameter $Λ$CDM + $N_\mathrm{eff}$ + $\sum m_ν$ + $T_\mathrm{AGN}$ + $A_\mathrm{kSZ}$ + $n_\mathrm{kSZ}$ model by less than 7%. The small impact on cosmological parameters reflects the strong constraining power of CMB polarization, the temperature-polarization cross spectra, and CMB lensing spectra reconstructed from polarization-only estimators, all of which are minimally affected by extragalactic foregrounds. In particular, the survey remains a sensitive probe of light thermal relic particles, achieving $σ(N_\mathrm{eff}) = 0.0154$, which can exclude any new species ($ΔN_\mathrm{eff} \ge 0.027$) with at least 90% confidence. Our simulation and foreground-removal codes are publicly available and should aid the development of analysis pipelines for ultradeep, ultrahigh-resolution microwave surveys.

Nonextensive Bose and Fermi state is well studied with counting near black-hole horizons in a one-function, nonextremal, static, spherically symmetric class. The calculation is organized around a $q$-deformed grand ensemble whose total partition function is a $q$-product of single-level partitions. This identifies the bosonic occupation-number sum as the primary level partition and allows us to separate two operations that are often conflated: replacing that sum by a geometric-series form and applying the ordinary reciprocal-log identity to a $q$-logarithm. A finite-energy tunneling prescription yields a leading near-horizon state-counting term controlled by the surface gravity and a horizon-response factor. Entropy matching then gives an energy scale proportional to the Hawking temperature; the geometry factor cancels from normalized nonextensive corrections and from the Bose--Fermi ratio. We derive the first-order coefficients analytically and evaluate the unexpanded $q$-integrals numerically. The linear approximation is accurate to better than one percent when $q$ lies within roughly two percent of unity. All prescriptions recover the extensive limit, but the exact bosonic level sum and the geometric approximations give opposite signs for the leading deformation of the conventional $7/8$ Bose--Fermi relation. The result isolates a robust geometry-factorized structure and quantifies a distinct $q$-algebra systematic uncertainty.

Null-congruences encode the causal structure of spacetime, describing a phase space of radiation channels, whose expansion scalar $θ$ satisfies the Raychaudhuri equation. In spherical symmetry, the expansion scalar of wave fronts of area $A=4πr^2$ satisfies $θ= A^\prime/A = 2/r$. About a mass $M$, this carries an encoding area $A_E = λ\varphi \ell_p^2 = λR_g r \le A$ in a UV-IR consistent coupling $\ell_p^2= G\hbar/c^3$ to its Compton phase $\varphi = (mc/\hbar)r \propto 1/\hbar$. Here, $λ= 8π$ is fixed by the Newtonian limit at large distances and the saturation limit $A_E=A$ at horizon about $M$ based on the non-local extension of the equivalence principle in Rindler spacetime, without using the Einstein equations. In spherical symmetry, this recovers the Schwarzschild metric from a geometric embedding of $p=A_E/A$, which is directly observable in gravitational lensing. The result shows a non-thermal origin of the scaling dimension of two of gravitational phase space, previously derived based on the second law of thermodynamics by Bekenstein (1973), with possible implications for cosmological spacetime.

The LUX-ZEPLIN (LZ) Collaboration has reported a nuclear-recoil candidate at $248$ keV, motivating interpretations in terms of endothermic dark matter (DM). We investigate this signal in Higgs-coupled minimal dark matter (HC-MDM), where the same electroweak interaction responsible for the terrestrial recoil also induces DM capture in the Sun and a potentially observable high-energy neutrino signal. We independently determine the mass splittings required to reproduce the LZ event and perform an improved calculation of Solar capture and post-capture evolution. Our treatment includes natural isotope mixtures, exact inelastic kinematics, nuclear form factors, finite-temperature nuclear motion near kinematic closure, orbital cooling through loop-induced elastic scattering, and an explicit capture--annihilation equilibrium test. The resulting Solar predictions are numerically stable and robust against the astrophysical and nuclear variations considered. For the two thermal benchmarks within the released IceCube DM mass grid, $3_M2_D$ and $5_M4_D$, the predicted annihilation rates exceed the mass-matched $W^+W^-$ IceCube response by factors of approximately $2.2\times10^3$ and $91$, placing both standard-halo solutions under strong tension within the $W^+W^-$ response mapping used here. For the heavier representations, whose masses lie above the published $10$ TeV IceCube grid, we construct a response-based extrapolation to higher masses. The $22$ TeV $7_M6_D$ benchmark remains in strong tension with IceCube even under our conservative response-loss extrapolation, while the $48$ TeV $9_M8_D$ benchmark also remains above the extrapolated response, although much closer to the sensitivity boundary. The $82$ TeV $11_M10_D$ benchmark is response dependent, whereas the $130$ TeV $13_M12_D$ benchmark remains below the extrapolated IceCube sensitivity.

Population III (Pop III) stars are expected to produce extremely hard ionizing spectra, yet direct evidence for their signature--strong narrow He II $\lambda1640$ emission with weak or absent metal lines--remains elusive. Although Pop III formation peaks at $z\gtrsim10$, models predict that nearly pristine gas pockets can survive to later times, making intermediate-redshift searches a probe of metal mixing across cosmic time. We search 109,545 high-confidence Ly$α$-emitting galaxies at $1.9<z<2.3$ in the Hobby--Eberly Telescope Dark Energy Experiment's database for systems with strong He II $\lambda1640$ and no detected UV metal lines, recovering eight candidates and a first-order, tentative comoving number density of $\sim30$ Gpc$^{-3}$. The eight-object stack yields a rest-frame He II equivalent width of $28.7\pm6.7$ Angstrom, He II/Ly$α=0.354\pm0.094$, and a resolution-corrected He II FWHM of $449\pm105$ km s$^{-1}$. We find no significant N V, C IV, or O III] emission; using the He II width to set the scale, we derive $3σ$ upper limits of N V/He II $<0.167$, C IV/He II $<0.128$, and O III]/He II $<0.140$. This combination of strong He II, elevated He II/Ly$α$, moderate line width, and weak metal lines is difficult to reproduce with metal-enriched stellar populations, AGN narrow-line regions, shocks, or classical Wolf--Rayet features. We interpret the sample as candidates for very metal-poor or Pop III-like ionizing sources at $z\sim2$. Definitive confirmation requires deeper rest-UV spectroscopy to measure He II and metal-line limits in individual objects, followed by rest-optical spectroscopy of H$β$, [O III] $\lambda4959,\lambda5007$, [O II] $\lambda3727$, and H$α$ to constrain gas-phase metallicity, ionization conditions, and AGN activity.

The Kepler mission has significantly advanced our understanding of exoplanet demographics. However, Kepler planet candidates showing centroid shifts are often classified as false-positives and thus excluded from demographic analyses. Some of these systems may be erroneous false positives, resulting from unresolved binary stars with planets orbiting one of the stars, thereby introducing biases into past exoplanet demographic studies. Here, we present updated Kepler centroid offset values derived from Gaia Data Release 3 (DR3) and analyze KOI-1623.01, a Kepler centroid offset false positive (COFP) with a clear transit signal and an orbital period of 111 days. High-contrast imaging with Keck/NIRC2 reveals that KOI-1623A has a close companion with a separation of $0.704 \pm 0.002$ arcsec, consisting of a slightly evolved G-type primary and a K-dwarf secondary, while Gaia detects a tertiary companion with common proper motion at a separation of $\sim 8.4$ arcsec. By combining transit photometry, Kepler pixel level data, direct imaging, and Gaia DR3, we find that the transit signal is consistent with a $10.02 ^{+0.44}_{-0.31} R_{\oplus}$ giant planet located in the habitable zone of the secondary. We furthermore use resolved companions from Gaia to estimate that another $\approx 70$ out of $515$ of Kepler's centroid false positives were misclassified and instead are real planets that belong to binary systems.

We study the relaxation of a complex massive scalar field around a supermassive black hole, modeling the near-horizon dynamics of Fuzzy Dark Matter (FDM) with spatial granularity. The field is initialized as a broadband, anisotropic multi-mode random field in momentum space with various spectral widths $σ$, decomposed into spherical harmonics up to $\ell_{\max}=100$, and numerically evolved on a finite domain using a specialized radial basis adapted to the background spatial geometry. We track the evolution of each multipole and calculate its Noether flux across both the event horizon and the outer boundary. Since angular momentum barriers suppress the horizon absorption of high-$\ell$ modes, while high-frequency radial components and fine-scale angular structures dissipate via outward radiation and horizon accretion, the Schwarzschild spacetime acts as an effective low-pass filter. Regardless of the initial spectral width, all configurations converge toward a universal late-time relaxation state dominated almost exclusively by low multipoles ($\ell \le 2$) and small radial wavenumbers, with the initial bandwidth $σ$ determining the timescale of global charge depletion. For concrete physical estimates, we set the boson mass to $m_b=10^{-22}\,\mathrm{eV}/c^2$ and the black hole mass to $M_{\mathrm{BH}}=6.5\times10^9 M_\odot$.

We present near-infrared spectroscopic observations of L1415-IRS, a Class I type young stellar object that underwent a factor of $>$30 photometric outburst in the early 2000s. Ours is the first reported spectroscopy of this deeply embedded source. Both components of L1415-IRS present continuum plus $H_2$ emission. The southern component has additional weak HI Br$γ$ and $CO$ emission. Based upon our spectrum and evidence of a declining lightcurve, we hypothesize that either: the outburst was not large enough and therefore hot enough to produce the FU Ori type absorption spectrum commonly interpreted as an accretion disk atmosphere, or the source has already recovered from FU Ori-type equatorial accretion to re-establish T Tauri-type magnetospheric accretion.

A narrow gamma-ray line would be a distinctive signature of dark matter, but its interpretation depends as much on its spatial distribution as on its energy. We reanalyze the reported 43 GeV feature in Fermi-LAT observations of Virgo, Fornax, and Ophiuchus. The excess is reproduced in the same public data and survives event-type partitions, a test using data withheld from the original selection, and detector-coordinate permutation tests. Realistic hadronic and inverse-Compton emission is too broad to explain it. A spatial--spectral likelihood analysis restricted to Good Time Intervals and including diffuse backgrounds and catalogued sources finds a broad line-like component. Uniform surface brightness and substructure-enhanced annihilation are preferred, whereas smooth NFW annihilation and a central point source are unlikely explanations. An independent annular analysis agrees: the events associated with the excess extend through much of the virial region. Calibration over 20--70 GeV and five morphologies gives a conditional global significance of approximately $2.7σ$, peaking at 44.5 GeV for uniform brightness. This result is conditional on the preselected three-cluster sample. The broad morphology is the main obstacle to a dark-matter interpretation. Such an interpretation would require photons to be a leading annihilation channel, annihilation to remain efficient in cold subhalos, and subhalos to supply nearly all the cluster luminosity. Yet a boost consistent with $Λ$CDM drives the required two-photon cross section into the range most tightly constrained by Galactic-halo line limits. The most likely explanation is a chance fluctuation amplified by analysis and target-selection effects, possibly compounded by residual Galactic diffuse or unresolved-source mismodeling; if celestial, its morphology makes a simple dark-matter interpretation unlikely.

In this literature, we derive an approximate black hole metric within the \( f(R, G, T) \) gravity framework, where the black hole is surrounded by an anisotropic fluid as the matter source. This approach allows us to study the influence of anisotropic pressures on the black hole structure. We then investigate the thermodynamic properties of the approximate solution, considering various types of matter fields. We analyze the thermodynamic quantities such as Wald's entropy, Hawking temperature, and specific heat for each case. After examining the specific heat, we find that the black holes are thermodynamically stable, with positive specific heat values indicating stability. This result suggests that the black holes under these conditions do not undergo thermal instability, ensuring their long-term viability.

Context. Accurate stellar spectra are essential to derive stellar properties. Traditional model atmospheres often oversimplify phenomena that break spherical symmetry, such as rotational deformation or multiplicity. The Spectroscopic PAtch Model for Massive Stars (SPAMMS) accounts for these effects, but its applicability has been limited by the spectral types covered by its model atmosphere grids. Aims. We aim to extend the parameter space of the model atmosphere grids available to SPAMMS, enabling spectral synthesis across a broader range of stellar types. Methods. We computed specific intensities, $I\left(λ,μ\right)$, for $101$ emergent angles using PRISMAS (Pipeline of Radiative Intensity Synthesis for Meshed Atmospheric Surfaces), and pre-computed LTE and non-LTE atmospheres from two ATLAS9-Kurucz grids and the TLUSTY-based OSTAR2002 and BSTAR2006 models. Results. The intensity grids cover effective temperatures from $3500$ to $55000\,\mathrm{K}$ and surface gravities from $0.0$ to $5.0\,\mathrm{dex}$, spanning O- to K-type stars. They include metallicities from $0$ to $30\,\mathrm{Z_\odot}$ and microturbulent velocities of $1$, $3$, $5$, and $10\,\mathrm{km\,s^{-1}}$. The spectral range extends from $3000$ to $9000\,\mathring{\mathrm{A}}$ with $Δλ=0.01\,\mathring{\mathrm{A}}$. As a proof of concept, we modelled with SPAMMS a rapidly rotating B-type star and an eclipsing Algol-type binary. Conclusions. The new LTE-Kurucz and NLTE-TLUSTY grids substantially expand the parameter space accessible to SPAMMS. The code can now generate synthetic spectra for a broader range of stellar types and geometries, including rapidly rotating and multiple systems -- providing a more comprehensive framework for modelling non-spherical stellar surfaces.

The cosmic web governs the growth of dark matter halos through anisotropic matter flows, yet the role of local velocity-field dynamics in halo mass assembly remains poorly understood. The velocity field carries non-linear information about matter transport that is not fully captured by density alone, particularly within filaments. We investigate this connection through velocity divergence and vorticity using dark matter halos in the IllustrisTNG TNG50-1-Dark simulation at $z=0$. Cosmic filaments are extracted with the 1-DREAM framework, and we statistically study how halo mass relates to divergence, vorticity, and radial distance from filament spines across three overdensity regimes. The velocity field correlates strongly with the cosmic web, with converging flows tracing filaments and vorticity arising primarily in overdense, non-linear regions. Within filaments, both divergence and vorticity show enhanced kinematic activity at intermediate distances from the spine, with amplitudes increasing with overdensity. We confirm a clear mass segregation, with more massive halos preferentially located near filament centers. Most notably, halos with masses $M \geq 10^{12} M_\odot$ are found almost exclusively in regions of low absolute divergence and low vorticity, while highly dynamical environments are dominated by low-mass halos. We further find that halos residing in kinematically calm regions assembled later and maintain higher recent accretion rates than those in dynamically active environments, linking present-day velocity-field properties to halo growth histories. These results demonstrate that the local kinematic state of the cosmic flow encodes information about halo mass and assembly beyond that contained in density alone, providing a complementary dynamical description of halo growth within the cosmic web.

Photometric galaxy clustering is a major probe of large-scale structure whose interpretation relies on accurate characterization of the galaxy redshift distribution, $n(z)$. Standard analyses assume the $n(z)$ of a selected lens sample to be isotropic, whereas observational systematics can imprint spatial variations. We quantify how this anisotropy biases angular galaxy clustering, internal consistency, and cosmological inference. We develop an analytical formalism for clustering under spatially varying $n(z)$ and forward-model realistic observational systematics, galaxy selection, and tomographic binning. For a fiducial magnitude-limited sample, modeling $w(θ)$ with a single global $\bar n(z)$ underestimates the clustering signal by up to $\sim10\%$ at small scales in the highest-redshift bin, decreasing towards lower redshift. In a $2\times2$pt analysis, cosmological parameter shifts remain below $1σ$ for KiDS-like and DES-like setups. For an LSST-like setup retaining KiDS-like selection and observing-condition variation, $σ_8$ is biased low by $1.7σ$ and $Ω_m$ high by $1.4σ$, while the high-redshift galaxy bias shifts by up to $\sim4σ$. The anisotropy also affects whether the probes are consistent enough to combine: a calibrated posterior-predictive internal-consistency test finds the contaminated angular clustering improbable given the joint cosmic-shear and galaxy--galaxy-lensing constraints (median calibrated $\tilde p\approx8\times10^{-4}$ under LSST-like noise). Anisotropic redshift distributions are therefore a potentially significant systematic for future photometric clustering analyses. Our forward model can be used to mitigate it given a realistic selection function and galaxy population model.

The recent observation of a single event at nuclear recoil energy $248\pm23\pm23$ keV by the LZ collaboration has prompted dark matter (DM) model builders to explore new physics explanations. While most of the existing literature invokes an endothermic scattering by a TeV-scale DM to explain the event, the favored DM interaction strength is in tension with the indirect detection constraints for most of the models. To evade this, this work explores exothermic DM scattering in a hadrophilic $q_1-q_2$ gauge extension setup with a new $Z'$ boson as a mediator to the dark sector. This can naturally escape the indirect searches in the absence of tree-level coupling of the mediator to neutrinos, $W,τ$, and $b$ quark. We consider freeze-in production at a low reheating temperature to populate DM in the early universe, with a comparatively smaller coupling than the thermal freeze-out scenario. We identify the $Z'$ parameter space that can explain the observed relic density and LZ event simultaneously. We also discuss the relevant constraints coming from flavor mixing and FCNC.

In this paper, we derive the complete transformations of a generic first order perturbative metric in Schwarzschild spacetime to the Newman-Unti (NU) gauge in series expansion near null infinity. This allows us to determine the asymptotic shear, the mass aspect and the angular momentum aspect of the perturbative fields. As a direct application, we derive the corresponding asymptotic NU data for quasinormal modes (QNMs) of a Schwarzschild black hole. The total energy and angular momentum of QNMs obtained from the asymptotic NU data vanish but classical supertranslation charges are non-trivial, which can be applied to fix the BMS frame at first perturbative order.

Many well-motivated dark matter (DM) scenarios, including dark photon DM, axion-like particles, sterile neutrinos, and primordial black holes, can deposit energy into ordinary matter and significantly raise the thermal energy of astrophysical gas. We present a simple, model-independent implementation of DM-induced heating in the RAMSES hydrodynamic simulation code, parameterized by an effective local heat-transfer rate $Γ_{\rm heat}$. Using simulations of an isolated gas-rich dwarf galaxy with properties similar to WLM, we find that DM heating substantially alters the thermal structure of the circumgalactic medium (CGM) while leaving the star formation rate and the dense ISM nearly unchanged. For $Γ_{\rm heat} = 10^{-27}$ to $10^{-25}\ {\rm s}^{-1}$, the neutral hydrogen column at 10 kpc is suppressed by 0.6 to 3.9 dex relative to the control simulation. Similarly, the C II and Si II column densities decrease by 4.1 and 5.5 dex, and by 2.2 and 3.0 dex for C IV and Si IV. These signatures are accessible to quasar absorption-line spectroscopy, so CGM observations of nearby field dwarf galaxies provide a promising avenue for constraining dark matter models that inject heat into ordinary matter.

Thermal conduction can act as a heating mechanism for cool-core clusters, yet the same transport must erode any temperature discontinuity it crosses. Ophiuchus permits both effects to be tested in the same atmosphere. Its temperature rises from about $1$ keV in the innermost core to $9$ keV at $r\sim30$ kpc, making conductive heating especially favorable, while two sharply resolved cold fronts independently constrain the same transport. \textcolor{black}{We construct a steady radial energy balance from deprojected \textit{Chandra} density and temperature profiles, with XRISM-based limits on turbulent heating and a subsonic inflow contribution.} Closing the balance with conduction requires only sub-Spitzer coefficients --- \Ftwentyfive\ at 25 kpc and \Finner\ and \Fouter\ at the inner and outer cold fronts --- so conduction can energetically supply the missing heat. Yet the same coefficients, applied normal to the cold-front surfaces, would broaden the fronts to their observed width limits within \Tinner\ and \Touter\ Myr, short compared to the characteristic sloshing timescales that generate and sustain such interfaces. For a representative 200 Myr residence time, the radial conductivity demanded by the energy balance exceeds the cross-front conductivity permitted by front survival by roughly a factor of three, with the required anisotropy exceeding unity in 99.9\% and 98.8\% of accepted models. The result leaves two broad regimes. Thermal conduction is either subdominant in Ophiuchus or contributes substantially but is strongly anisotropic, with heat transport suppressed across the cold fronts relative to the radial direction. Together they disfavor locally isotropic conduction as a dominant heating mechanism in cool cores generally.

We present FOSSIL, a proposed space mission which is designed to deliver unprecedented absolute spectroscopic measurements of the sky intensity over a broad frequency range (50GHz to 2THz) to accurately measure chemical-potential mu-type, Compton y-type, and Compton relativistic y-type CMB spectral distortions. FOSSIL will reach monopole sensitivities representing a gain of roughly three orders of magnitude over COBE/FIRAS. These capabilities will allow us to constrain primordial density fluctuations on scales k~(10 to 10^4)Mpc^-1, test inflationary models, probe dark matter annihilation/decay/interactions, constrain primordial black hole seed scenarios, and pin down feedback models in galaxy evolution and structure formation by measuring the cosmic thermal energy content. Furthermore, its 130 spectral channels will enable legacy astrophysics, including precise characterisation of the Cosmic Infrared Background, line intensity mapping during cosmic noon, and detailed Galactic dust mapping.

Massive stars, including Wolf-Rayet stars, are predominantly found in binary systems. These systems are known to emit both thermal emission from stellar winds and, occasionally, non-thermal emission produced by relativistic electrons accelerated in the wind-wind interaction region. We intend to provide the most complete census of radio emission from WR stars in the Galactic plane, using the SARAO MeerKAT Galactic Plane Survey (SMGPS) complemented by the MeerKAT Galactic Center Survey. Our main motivation is to identify hints of synchrotron radio emission indicative of particle-accelerating colliding-wind binaries (PACWBs). We compiled an input catalogue of 428 WR stars positionally covered by the SMGPS and the MGCS. Using the survey data, we measured the radio emission at 1.3 GHz for detected WR stars. We also measured the upper limits for objects located in sufficiently low radio background regions. For detected objects, we searched for a radio excess by comparing the measured fluxes to two different evaluators of the thermal emission from massive star winds. We detected 23 targets and determined the upper limits for 279 WR stars. Among the detected objects, 15 display a (significant or potential) radio excess that cannot be explained by unresolved circumstellar emission of any kind. After removing already known PACWBs, we report the identification of 12 potential new PACWB candidates. Our study has led to the compilation of the most extensive catalogue of Galactic WR radio emission to date, based on a homogeneous dataset covering the Galactic plane. The low detection rate indicates either a low occurrence rate of synchrotron emission or substantial attenuation by turnover processes that are clearly dominated by free-free absorption from the WR wind material. Our results open the door to dedicated follow-up observations aimed at ascertaining the nature of the identified radio excesses.

The frequency spectrum of the cosmic microwave background is the most precise blackbody ever measured in nature, with deviations constrained at the level of almost one part per million from the COBE satellite. Nevertheless, departures away from a perfect blackbody are present in standard $Λ$CDM cosmology, lurking just beneath the surface of our current observational bounds. These spectral distortions provide invaluable information on our thermal history in both the post- and pre-recombination epochs, allowing us to peer beyond last scattering and into the primordial Universe. Here, we present an overview of the underlying physics responsible for generating CMB spectral distortions at all epochs. As an illustration of this rich physics, we review a comprehensive set of mechanisms capable of generating distortions both within and beyond the standard $Λ$CDM paradigm. We also discuss the information that can be gleaned by going beyond the monopole (sky-averaged) spectrum and exploiting the spatial information present in anisotropic spectral distortions. To supplement our discussion of the diverse science of spectral distortions, we provide an overview of the upcoming and proposed experimental landscape. We highlight that the combination of the TMS, COSMO, and BISOU experiments will provide the first discovery of a monopole $y$-type distortion within the coming decade. From space, the proposed FOSSIL experiment is forecasted to improve upon the original COBE/FIRAS measurement by roughly three orders of magnitude in sensitivity, bringing with it the detection of the $Λ$CDM $μ$-type distortion sourced by the dissipation of small scale acoustic modes in the pre-recombination plasma. With transformational measurements on the horizon, CMB spectral distortions offer a uniquely sensitive probe of the thermal history of the Universe at redshifts $z \lesssim 2 \times 10^6$.

Pulsations in $δ$ Scuti stars exhibit a wide range of excited frequencies and amplitudes, but the underlying principles governing their amplitude distribution remain unclear. We revisit this problem using a global, energy-based perspective. For a complete sample of 3200 bright $δ$ Scuti stars observed by TESS, integrated power spectral density (PSD) of the Fourier spectrum allows us to measure the overall pulsation signal. From this, we derive a proxy for the total pulsation energy. We uncover a bimodality among higher-frequency delta Scuti stars in their integrated PSD, possibly suggesting the presence of two distinct pulsation regimes. In addition, we identify an apparent upper limit to the total pulsation energy, indicating a constrained energy budget for these pulsations. This feature is not evident when stars are characterized solely by their dominant mode amplitude.

We present the first application of neural field-level baryon-acoustic oscillation (BAO) reconstruction to real spectroscopic survey data. We develop Linear Field Transformer (LiFT), a 3D vision transformer that takes as input the observed galaxy field, its standard reconstruction, and a set of context channels encoding local line of sight, survey coverage, and redshift, and learns to correct standard reconstruction toward the linear density field. We construct a forward-modeling pipeline to produce mock lightcones similar to the DESI Data Release 1 (DR1) luminous red galaxy (LRG) sample, and train LiFT on these. We validate LiFT on held-out simulations, as well as on additional mocks which differ in gravity solver, halo finder, HOD, cosmology, and fiber-assignment history, as well as on mocks analyzed with a distorted distance-redshift relation; ultimately, we find unbiased dilation parameters with consistently tighter constraints than standard reconstruction. Applied to the DESI DR1 LRGs, LiFT improves errors on $α_{\rm iso}$ by 13%, 24%, and 30% and on $α_{\rm AP}$ by 5%, 22%, and 30% relative to the DESI DR1 standard reconstruction analysis in the LRG1, LRG2, and LRG3 bins respectively; meanwhile, our DR1 central values remain consistent with DR1 and DR2. This equates to a factor of 1.2, 1.7 and 2.0 increase in Figure of Merit (or effective survey volume) if one were to only use standard reconstruction. Ultimately, these results establish LiFT as a validated, survey-ready tool for current and upcoming galaxy surveys.

The LUX-ZEPLIN (LZ) experiment has reported a high-energy nuclear recoil candidate, LZ230616, at $E_R = 248 \pm 23_{\rm stat} \pm 23_{\rm sys}\,$keV, which is difficult to reconcile with elastic scattering of halo dark matter (DM). Endothermic inelastic scattering offers a natural explanation, but it rests on a nucleon coupling that is off-diagonal rather than diagonal. We show that this feature arises automatically within a radiative Dirac neutrino mass framework. The Standard Model is extended by three right-handed neutrinos (RHNs), a pair of vector-like neutral fermions for each generations, and a scalar sector comprising an inert doublet and a real singlet, governed by a $\mathcal{Z}_2 \times \mathcal{Z}_4$ symmetry. The $\mathcal{Z}_2$ symmetry stabilizes the DM, while the $\mathcal{Z}_4$ symmetry forbids the tree-level Dirac Yukawa coupling together with all renormalizable Majorana mass terms, thereby allowing the generation of Dirac neutrino masses at the one-loop level through a softly broken scalar trilinear coupling `$κ$'. Since the neutral dark-sector fields can be expressed in terms of real scalar mass eigenstates, the $Z$ boson couples to them purely off-diagonally and elastic $Z$-mediated scattering is absent identically rather than simply suppressed. The singlet--doublet mixing induced by `$κ$' governs both the inelastic rate and the one-loop Dirac neutrino mass, which therefore vanish as $κ\to 0$. We find that DM masses in the few hundred GeV to TeV range, with splittings of $\mathcal{O}(340-360)\,$ keV, simultaneously reproduce the observed relic abundance, account for the LZ event and yield neutrino masses of the correct order, while respecting elastic direct-detection limits and all relevant theoretical and experimental constraints. A substantial part of the surviving parameter space lies within reach of DARWIN.

We investigate the loss of coherence of a two-level Unruh--DeWitt detector coupled to massless real scalar fields in $(1+3)$-dimensional de Sitter spacetime. Treating the detector--field interaction perturbatively to second order, we derive analytic expressions for the real asymptotic coherence-decay coefficients for conformally and minimally coupled scalar fields. For the conformally coupled field, the coefficient remains finite in the small-gap limit and grows linearly with the detector energy gap in the large-gap regime. For the massless minimally coupled field, the logarithmic sector of the Wightman function produces an additional positive contribution. This provides an explicit analytic extension of the previously known infrared enhancement of detector transition responses to the asymptotic coefficient governing detector coherence. The explicitly average-time-dependent sector does not contribute to the real nonzero-frequency asymptotic coefficient under the adiabatically regulated half-line prescription adopted here, when the regulator is removed at fixed average time and fixed positive detector gap. This asymptotic result does not describe the complete finite-time dynamics, which depends on the switching function and observation interval. Within the stated prescription, the coefficients satisfy $Γ_{\rm MMC}>Γ_{\rm CC}$ for $ω>0$, with $Γ_{\rm MMC}/Γ_{\rm CC} =1+(H/ω)^2$. The distinction is most pronounced in the small-gap regime, whereas the two coefficients approach one another for large detector gaps.

We discuss the effectiveness of the recently proposed quantity $J_{s}$ in suppressing the contribution of a string's self-field to its spectrum, which is used in field theory simulations to track the emission of axions by decaying cosmic strings. We compute the contribution of the self-field to the spectrum of this quantity and to the usually computed spectrum of $φ\partial_{t}α$ for an infinitely long straight string. Although we demonstrate that the $J_s$ approach is a substantial improvement, we also point out that this highly symmetric model doesn't capture the full contribution of the self-field to the spectrum emitted by strings of a more complex shape, which are typical in network simulations. We then illustrate this point numerically using a simulation of a sinusoidally perturbed straight string. For this simple configuration, we managed to separate the contribution from the self-field and the radiation to the spectrum of $J_{s}$, using the self-field subtraction method. This allows us to show that the spectrum of $J_{s}$ is still dominated by the $n=1$ mode of the string's oscillation, which can be attributed to the self-field and is largely suppressed when the self-field is removed. We also demonstrate that this mode is predominantly sourced by variations along the direction parallel to the string, which is missing in the unconnected segment model, used by arXiv:2512.13653 to claim the effectiveness of $J_{s}$ in suppressing the self-field's contribution to the spectrum.

In this paper we investigate accelerated gluodynamics for a broad intervals of temperature and acceleration. Our study is carried out within lattice simulation in the co-moving reference frame parameterized by the Rindler coordinates. We developed the renormalization prescription that allowed us to calculate renormalized Polyakov loop as a function of coordinate in the Rindler spacetime. Using the data for the renormalized local Polyakov loop, we calculated spatial dependence of the static quark free energy and effective mass of static quark. Besides the Rindler coordinates, it is believed that accelerated gluodynamics can be approximated utilizing non-accelerated gluodynamics with a properly adjusted temperature gradient in accordance with the Tolman-Ehrenfest law. We compared these approaches for the observables under study. It was found that they agree quite well close to the critical temperature and demonstrate disagreement at higher temperatures. We believe that this disagreement might be attributed to the Tolman-Ehrenfest law corrections which appear in the Rindler gluodynamics.

We explore the thermodynamic properties of static black holes in conformal gravities in arbitrary, even dimensions. In the spherically symmetric case, up to local Weyl rescalings, these spacetimes are static and unique, due to a Birkhoff's theorem that applies for all conformal gravities whose Lagrangians are complete contractions of $k$ Weyl tensors in dimension $D=2k$. The solutions can be chosen to be asymptotically locally (A)dS or flat, and the horizon can be extended to be a surface of arbitrary constant curvature, $γ$. We focus on the asymptotically AdS case, for which the cosmological constant at infinity is an integration constant, and compute the mass of the spacetime in terms of the diverse hair parameters which have a purely gravitational origin. We show that the first law of thermodynamics is satisfied, and study phase transitions in the canonical ensemble, showing the existence of zeroth-order black-hole/black-hole phase transitions for arbitrary dimension.

The Baade-Wesselink (BW) method compares the linear and angular variations of Cepheids to derive their distance. This method is limited, however, by the projection factor, which relates the observed radial velocity to the true pulsation velocity of the star. Using simultaneous observations from the CHARA Array interferometer in the K, H, and R bands with the MYSTIC, MIRC-X, and SPICA combiners, respectively, we aim to understand the physics of Cepheid atmospheres better. Applying a specific method to multi-chromatic simultaneous interferometric observations of Cepheids, we derived robust limb-darkened angular diameters that were then used to calibrate the SBCR and study the projection factor. We also developed a strategy to measure the limb-darkening of Cepheids in R, H, and K bands. These measurements were then used to constrain the geometrical component of the projection factor. From the limb-darkened angular diameter curves of three Cepheids, we decreased the scatter of the SBCR in $V-K$ colour to 0.0011 magnitude and to 0.0040 for the SBCR in $G_{BP}-G_{RP}$. These SBCRs are particularly robust because they are based on multi-chromatic diameters and are homogeneous, which previous calibrations for Cepheids were not. For the very first time, we derived the limb-darkening of Cepheids in R, H, and K bands. For $δ$ Cep, we derived an R-based projection factor of $1.275 \pm 0.051$, and the geometrical part obtained from the measured limb-darkening coefficient in R band was estimated to $1.420 \pm 0.016$ (on average), as expected from stellar static and hydrodynamical atmosphere models. The limb-darkening coefficients obtained in H and K band are consistent with models. These results demonstrate that multi-chromatic interferometry can improve the accuracy of the BW method. It is therefore essential to continue the CHARA survey of Cepheids in the coming years.

Whether classical spacetime remains welldefined at arbitrarily high curvature is a central question. In general relativity, singularities signal the breakdown of the classical description, motivating modifications of spacetime's short-distance structure. Pseudocomplex general relativity (pcGR) extends spacetime geometry to pseudo-complex coordinates, naturally introducing two metric sectors, the physical metric gμν and an auxiliary field fμν. The magnitude of the pseudo imaginary component defines the invariant acceleration scale a0, a direct consequence of the pseudo-complex structure. The simultaneity condition, requiring both idempotent sectors to satisfy the pseudo-complex Einstein equations independently, fixes the auxiliary field algebraically, with no new propagating degrees of freedom, and sets a lower bound on the lapse that regularizes curvature. The regularization is achieved through the combined effect of the lapse-gap condition eν(r) gt. a0, which ensures the radial metric component remains nondegenerate, together with the regularity conditions at the areal-radius origin, B(0) = 1 and B'(0) = 0, which emerge naturally from the pseudo-complex geometry. These conditions jointly eliminate the Schwarzschildtype curvature divergence, rather than the gap condition acting alone.

We identify the spatial projection of a conserved string charge density with a bit-thread flow. This correspondence gives bit threads a string worldsheet description and allows us to compute black hole entropy from the maximal entropy flux of a worldsheet-derived current. We first reproduce the BTZ entropy and then extend the construction to the ten-dimensional F1-NS5-P and D1-D5-P systems. Under Type IIB S-duality, the F1 current maps to the D1 current. Although the local currents, metrics, and flow norms change, the maximal physical flux remains unchanged and gives the same entropy. This provides a nontrivial check of the correspondence. Finally, we propose a D1-D5 boundary-matching conjecture. In the coarse-grained picture, each string-charge line defines an entropy-carrying tube. The boundary CFT entanglement entropy and the bulk black hole entropy both fix the product of the number of tubes and the maximum entropy carried by each tube. For a fixed tube discretization, they assign the same capacity to each tube. They therefore describe the same conserved maximal entropy flux, providing an information-flow interpretation of the equality between entanglement entropy and black hole entropy.

We investigate how neutrino-mass constraints from cosmology depend on the assumed thermal history of the universe. Photon injection after neutrino decoupling would decrease the neutrino abundance inferred from the temperature of the cosmic microwave background (CMB), potentially loosening the upper limit on their masses. We first evaluate how the cosmological neutrino-mass bound is altered by the decay of massive particles into photons and dark radiation after Big Bang nucleosynthesis. To test the robustness of this constraint more generally, we also assess the impact of varying the temperature of the cosmic neutrino background without assuming a specific physical mechanism. We perform MCMC analyses of both frameworks with primary CMB observations from Planck, CMB lensing measurements from Planck and ACT, and baryon acoustic oscillation data from DESI. For the degenerate mass ordering, the $95\%$ credible limit tightens from $\sum m_ν<0.0691$ eV in a standard thermal history to $\sum m_ν<0.0652$ eV if the massive particles decay only into photons, while the injection of dark radiation in addition to photons slightly relaxes this limit to $\sum m_ν<0.0710$ eV. The same pattern holds for the normal and inverted orderings, and the decay scenario shifts the bound on the sum of the neutrino masses by at most $0.004$ eV for fixed mass splittings. Allowing model-agnostic changes in the neutrino-to-photon ratio yields a $95\%$ credible limit of $\sum m_ν<0.0724$ eV for the degenerate ordering, indicating that the stringency of our neutrino-mass bounds is not driven by constraints on the decay scenario. We find that the neutrino temperature and the sum of the neutrino masses are positively correlated, which implies that reducing the pre-recombination radiation density will only worsen the emerging tension between cosmological bounds on the neutrino masses and the measured mass splittings.

Gravitational sums over spacetime topology suggest a baby-universe field theory whose quanta are entire universes. We develop such a description from group-averaged amplitudes that impose gravitational constraints and define physical inner products in the single-universe theory. To sharpen the contrast with ordinary QFT, we start with gravity on one-dimensional worldlines, which can arise from reducing higher-dimensional cylinders. For multiple noninteracting worldlines that pair boundary components, the baby-universe Hilbert space is the symmetric Fock space over the physical one-universe Hilbert space, and the universe field expands in the on-shell modes spanning the latter. Each mode multiplies the sum of creation and annihilation operators, so the mode coefficients commute. The mode expansion thus makes explicit the field's commutativity and its Wheeler-DeWitt linear field equation. Diagonalizing the coefficients recovers a statistical description by a Gaussian path integral over on-shell field configurations labelling superselection ($α$-)sectors. We take exploratory steps toward more general topologies. Disk-like topologies coherently shift the no-boundary state, turning on the Hartle-Hawking wavefunction as a one-point function. Even with nontrivial topologies represented by interacting Feynman diagrams, we argue that group averaging near each boundary still enforces the same linear field equation, while interactions instead enter through a non-Gaussian measure over $α$-sectors. In the topological Marolf-Maxfield model, we promote a disk-and-cylinder truncation to an exact Fock-space representation of the full theory incorporating all topologies. We also relate the truncated and exact inner products by an insertion of interactions, whose perturbative expansion produces a sum over time evolutions generated by pairs of pants, evoking an analogy to group averaging.

We employ the statefinder hierarchy along with the fractional growth parameter to investigate the recently proposed noninteracting Viaggiu holographic dark energy (VHDE) model with the future event horizon as the infrared cut-off. By adjusting different values of the VHDE model parameter $δ$, we illustrate the evolution trajectories of several important parameters, prticularly, the statefinder hierarchy $S_3^{(1)}$, $S_3^{(2)}$, $S_4^{(1)}$, $S_4^{(2)}$ and the fractional growth parameter $ε$ against the redshift $z$ and the fractional energy density of matter, $Ω_m$. Together with the growth rate of matter perturbations, we present the statefinder hierarchy that defines a composite null diagnostic (CND), which can differentiate dynamical dark energy from $Λ$CDM. The joint use of the statefinder hierarchy and the fractional growth parameter offers a powerful diagnostic tool for probing VHDE, effectively resolving the degeneracy between different parameter choices within the model. We find that employing a CND pair, rather than a single diagnostic, significantly enhances the diagnostic power for the VHDE model. Finally, we perform an $ω_d-ω_d'$ analysis and, in addition, assess the squared sound speed $v_{s}^{2}$ to determine the model's dynamical stability under small perturbations.