2026-03-23 | | Total: 61
The scale-dependent bias in halo and galaxy power spectra is a key signature of local primordial non-Gaussianity (local PNG), with PNG sensitivity scaling as $b_φ/b_1$ -- the ratio of their responses to long-wavelength primordial potential $b_φ$ and late-time density fluctuations $b_1$. For number density fluctuations, these responses are closely tied by the universality relation, limiting the achievable ratio. We show that size density fluctuations strongly violate this relation, thus evading the limit. For galaxy-mass halos, sizes have a vanishingly small density response but a sizable, negative local PNG response, implying an effective $b_φ/b_1$ that is large in magnitude and opposite in sign to that of number counts. This makes galaxy sizes complementary probes of local PNG from the same galaxy sample, without any sample split. For a DESI-like survey, a multi-tracer analysis combining galaxy numbers and sizes improves the local-PNG detection significance by a factor of $\sim\!3.6$. Due to the sign flip, the number-size cross power spectrum further provides a handle on systematics in the event of a detection.
A new measurement of the Hubble constant $H_0$ is presented using the statistical dark siren method applied to a sample of seven well-localized gravitational-wave (GW) events from the fourth LIGO-Virgo-KAGRA (LVK) observing run and ten additional events from the first three runs. Galaxy catalogs from the DESI Legacy Imaging Survey (LS) are combined with a deep learning model to compute photometric redshift probability density functions. We extend our previous analysis by including the events GW230731_215307 and GW230927_153832, using sky maps from the fourth Gravitational-Wave Transient Catalog (GWTC-4), and introducing key methodological improvements: $r$-band luminosity weighting of host galaxies; an extended GW likelihood that incorporates information from the binary black hole component masses; and a consistent treatment of selection effects that accounts for the incompleteness of the magnitude-limited LS galaxy catalog. Using a total of 17 well-localized dark sirens (seven from the first part of the fourth observing run, O4a), we obtain $H_0 = 78.8^{+14.6}_{-12.2}$ km/s/Mpc without luminosity weighting and $H_0 = 78.2^{+12.0}_{-11.0}$ km/s/Mpc when applying $r$-band luminosity weighting. Finally, we combine the luminosity-weighted dark siren sample with the bright siren GW170817, including constraints on the jet viewing angle and corrections for the host galaxy peculiar velocity, to obtain a final constraint of $H_0 = 69.9^{+4.1}_{-4.0}$ km/s/Mpc, representing an improvement of approximately 11% in the uncertainty relative to the GW170817-only result.
Satellite galaxies in the Local Group tend to be distributed in thin, planar configurations, with many sharing coherent orbital motion. Galaxy formation simulations in $Λ$CDM have historically struggled to produce similar structures, leading to the so-called "planes of satellites problem". In this work, we investigate whether the emergence of such structures is connected to the mass of a major merger at $z\sim2$, analogous to the Gaia-Sausage-Enceladus (GSE) event in the Milky Way. We use the VINTERGATAN-GM suite of high-resolution zoom-in simulations, comprising five realizations of the same Milky Way-mass halo generated through targeted genetic modifications of a GSE progenitor. The GSE-like merger mass ratio is systematically varied from 1:10 to 1:2.1, while keeping the final dynamical mass and large-scale environment fixed. We find a clear and consistent trend: more massive GSE-like mergers lead to satellite populations that are both more planar and more kinematically coherent. In particular, simulations with merger mass ratios larger than 1:6 develop Kinematic Persistent Planes (KPPs), in which at least 40% of satellites co-orbit around a common axis over extended periods, comparable to those observed in the Milky Way. These structures arise when sufficiently massive mergers, accreted along the direction of maximum compression of the Lagrangian volume, produce flattened host halos with anisotropic velocity dispersions aligned with the merger direction. The merger aligns the host halo's minor axis with the direction of flattening of the surrounding cosmic web, and planes of satellites then emerge through two complementary processes: (i) satellites preferentially infall along the host's equatorial plane, and (ii) anisotropic dynamical friction in the non-spherical halo gradually reshapes their orbits toward this plane, generating coherent and long-lived planar configurations.
Maps of the sky in millimeter wavelengths contain rich information on cosmology through anisotropies of the cosmic microwave background (CMB). Creating multifrequency sky maps of anisotropies in the $I$, $Q$, and $U$ Stokes parameters is one of the first steps of CMB cosmology analyses. In this work, we describe the production and validation of a set of sky maps from the South Pole Telescope's third-generation camera, SPT-3G. The maps are from data taken in frequency bands centered at 95, 150, and 220 GHz and taken during the first two years, 2019 and 2020, of the SPT-3G Main survey, which covers $4\%$ of the sky. We applied high-pass filters to time series of individual detectors and binned the filtered time series samples into map pixels. After that, we calibrated and cleaned the maps to reduce known systematic errors. In addition, we searched for other systematic errors through null tests and mitigated a significant systematic error detected therein. The white noise levels of the full-depth maps of the $I$ Stokes parameter are $5.4$, $4.4$, and $16.2$ $\mathrm{μK}$-$\mathrm{arcmin}$ in the 95, 150, and 220 GHz bands, respectively, and $8.4$, $6.6$, and $25.8$ $\mathrm{μK}$-$\mathrm{arcmin}$ for $Q/U$. These maps are the deepest to date used for measurements of mid-to-high-$\ell$ primary temperature and $E$-mode polarization CMB anisotropies, and reconstructions of the CMB gravitational lensing potential. We make these maps and supporting data products publicly accessible.
We compute the full-sky angular power spectrum and bispectrum, along with their Fisher matrices, to forecast constraints on cosmological parameters for the BINGO and SKA1-MID Band 2 radio telescopes. This represents the first forecast analysis using the full-sky relativistic bispectrum in redshift space for these surveys. Our results show that the second-order velocity contribution, often neglected under the Limber approximation, accounts for approximately $24\%$ of the total signal at low redshifts, indicating that it must be included for accurate modeling. Using these forecasts, we find that while the bispectrum provides constraints comparable to the angular power spectrum for $Λ$CDM and ${\rm w}$CDM models, it becomes a powerful probe of dynamical dark energy. Restricting the analysis to linear scales, we show that the inclusion of the bispectrum yields a substantial improvement in the determination of the Chevallier-Polarski-Linder (CPL) parameters. In particular, the joint analysis of the bispectrum, power spectrum, and Planck CMB data improves constraints on ${\rm w}_0$ and ${\rm w}_a$ by over $70\%$, and the Hubble parameter $h$ by approximately $60\%$. These results underscore the importance of relativistic bispectrum for breaking parameter degeneracies and probing the nature of dark energy with upcoming large-scale structure surveys.
The Potential Field Source Surface (PFSS) model is the most used approach for extrapolating the global coronal magnetic field, offering efficiency and strong performance at large scales. However, PFSS assumes a potential coronal field, so it cannot account for distortions from electric currents. More advanced methods, such as nonlinear force-free field (NLFFF) models, can represent these effects but are much more computationally intensive. Recent observational techniques also allow 3D reconstruction of coronal loops, which trace solar magnetic field geometry. This work develops a new approach that constrains the PFSS model using 3D coronal loop information, improving agreement with observations while keeping efficiency. The model is based on PFSS field constraints from photospheric data but allows magnetic field deviations from the potential state within loop-influenced regions, maintaining control over divergence and force-freeness. We adapted NLFFF optimization to the PFSS framework, enabling multiple physical constraints. Our functional includes up to three terms: divergence-free, loop geometry, and force-free. The resulting Python algorithm was tested with synthetic loops, using Carrington rotation 2284 as the lower boundary. This method yields magnetic field solutions that better match the geometry of included loops and controls divergence and force-freeness. Our results show that 3D coronal loop information can be incorporated into PFSS, largely preserving computational efficiency even with many loops. This approach lets PFSS better reflect observed coronal structures without significant computational cost.
The effects of the polarization characteristics of beam trains in optical long-baseline interferometers are well known and have led to difficulties in measuring the spatial coherence of astronomical sources in the past. This has been overcome by designing symmetrical optical trains. With the advent of interferometers using large telescopes, observations of faint sources with high degrees of polarization have become even more possible. As in the radio domain, where radiation processes usually lead to high polarization rates, a description of coherence for polarized or unpolarized sources observed with non-polarization neutral interferometers is necessary. A theory of optical long-baseline interferometry fully taking into account the polarization characteristics of beam trains and those of the sources is presented in this paper, building on concepts developed for radio aperture synthesis. The concept of generalized Mueller matrix is introduced for the case of multi-aperture interferometers leading to a simple matrix relationship between the observed Stokes visibilities, as they are disturbed by the instrument polarization characteristics, and the object Stokes visibilities. This relationship is applied to the case of single-mode interferometers. The formalism also shows that classical complex visibilities (squared moduli, phases and closure phases) need to be debiased from polarization crosstalk, even when the source is not polarized as in this case ghost polarized visibilities are created.
The orbital and eccentricity evolution for compact object binaries through gravitational wave emission first derived by Peters and Mathews are used extensively throughout the gravitational wave community for calculating the orbital evolution and merger time of compact binaries. While improved calculations of the binary merger time have been the focus of several investigations since, the orbital evolution has not received the same attention. As the equations lack a closed form solution, a numerical integrator is required, but standard methods typically break when the point of merger is overstepped. We present a rewrite of Peters' equations in $\ln$-space, which allows common numerical solvers to converge. This leads to a more numerically robust and computationally efficient method for evolving compact binaries due to gravitational wave emission, reducing the number of function evaluations by 60\% to 70\% in our tests.
In this Letter, we present early emission data analysis of the tidal disruption event TDE2025aarm, including optical, UV and X-ray data. At a redshift of z = 0.01368, TDE2025aarm is the second closest TDE ever discovered, offering an unprecedented opportunity to study such phenomena in great details. We observed TDE2025aarm in optical with the Liverpool Telescope for a total of three epochs, and complemented our dataset with ancillary spectroscopic and photometric data. The early optical spectra are characterized by a blue-continuum and helium, hydrogen and possibly Bowen lines typical of H+He events. The optical light curves peak at M_g ~ -18.63 mag and are well described by fallback of a M_star ~ 0.16 M_sun star onto a M_BH ~ 2x10^{7} M_sun black hole. We report Swift-XRT detection in the 0.3-10 keV range, with a total flux of F_X ~ 1.42x10^{-14} erg s-1 cm-2, fitted by a black-body with kT ~ 0.39 keV. This makes TDE2025aarm a new event among optical/UV bright TDEs detected in soft X-rays. Our analysis suggests that the early emission from TDE2025aarm is powered by circularization shocks, and that the delayed accretion scenario best describes the observed features.
RS Canum Venaticorum (RS CVn) close binaries, characterized by tidal locking, rapid rotations, and strong magnetic fields, are ideal laboratories for high-resolution radio observations to probe emission processes, magnetic field configurations, and interaction activity. Despite their importance, only a few RS CVn sources have been explored by polarimetric observations of very long baseline interferometry (VLBI). To expand the effort, we have analyzed the existing Very Long Baseline Array (VLBA) astrometric data for the RS CVn binary FF Ursae Majoris (FF UMa). In the 5GHz VLBA experiments conducted between 2021 and 2024, both total intensity and circularly polarized emission were clearly detected at six of seven epochs. The consistently high brightness temperatures (10^7 K) and the moderate fractional circular polarization (10%-30%) over about three years indicate that the radio emission is mainly produced by gyrosynchrotron radiation from mildly relativistic electrons in the highly-ordered magnetic field. The radio luminosities are also comparable to those of previously studied powerful RS CVn binaries and show a significant anti-correlation with fractional circular polarization. A mean centroid offset of 13.4 +/- 3.1 solar radii between the Stokes I and V emission was found across multiple epochs, indicating a possible additional contribution from the secondary star via a magnetically active corona, a giant magnetic loop, or significant interaction activity with the primary star in the quiescent state.
At a distance of 5.1 pc, the 70 Oph AB binary star system is one of the most favorable targets for future direct imaging and astrometry missions surveying mature, terrestrial planets. We present new radial velocities (RVs) obtained with the Planet Finder Spectrograph (PFS) on the 6.5\,m Magellan II Clay Telescope in Chile. We collected 499 measurements of 70 Oph A and 334 measurements of 70 Oph B during 2023--2025. Combining these data with decades of archival RVs and astrometry, we derive an updated orbital solution for the binary and dynamical masses of $0.88 \pm 0.004\,M_\odot$ and $0.73 \pm 0.003\,M_\odot$ for the primary and secondary components, respectively. We find that the long-term RV variability of both components is consistent with stellar activity modulated by rotation periods, and we detect no coherent planetary signals in either component. We place upper limits on any planets orbiting in the plane of the binary. The 27 yr RV baseline for 70 Oph A excludes Jupiter-mass planets interior to 5 au and reaches a sensitivity of $0.3\,M_{\rm Jup}$ at 1 au or $0.5\,M_{\rm Jup}$ at 2 au. For 70 Oph B, with PFS data we rule out planets more massive than $0.25$--$0.3\,M_{\rm Jup}$ inside 0.5 au. We show that stable S-type orbits around 70 Oph A extend to $\sim2.5$ au, covering the habitable zone. Thus, Saturn-mass planets or smaller on stable orbits in the habitable zone of 70 Oph A are allowed. Overall, our results provide important guidance for future planet searches around this stellar system.
The detection and characterization of potentially habitable exoplanets is one of the chief goals of astrophysics for the coming decades. Imaging in reflected light is well suited for characterizing Earth-like planets, as much can be learned about these planets in this wavelength range (i.e., ~0.3-2 μm). Several studies have been conducted to determine the abilities and limitations of reflectance spectroscopy, but most previous studies assumed a homogeneous atmospheric and surface composition. Here we investigate how heterogeneities in the atmosphere and surface of an Earth-like planet impact retrieval results. We extend the ExoReL retrieval framework to include a step function for retrieving wavelength varying surface albedo. We then use it to retrieve on visible-to-near-infrared spectra of realistic 3D Earth models with different surface features in view and varying cloud types/distributions synthesized with the Planetary Spectrum Generator. Including the ability to fit for wavelength dependent albedo mitigates degeneracies that arise when using 1D models to analyze 3D planets, and we recover an Earth-like planet in all cases. We detect surface albedo steps at ~0.7 and ~1.1 μm despite clouds, both when significant lands are in view and when the spectra are averaged to account for a longer integration time. Our findings support the application of the vegetation red edge as a biosignature in the context of the Habitable Worlds Observatory. This study highlights the importance of considering a range of-particularly wavelength-dependent-surface albedos when using reflectance spectroscopy to characterize Earth-like exoplanets.
We present the first spatially resolved characterisation of the cool intragroup medium (IGrM) in a spectroscopically confirmed galaxy group at z=1.167. Using 30 independent sightlines towards the gravitationally lensed galaxy SGAS J0033+02, we combine background light from an extended gravitational arc and various sources in the field to map the distribution and kinematics of diffuse, metal-enriched gas pertaining to the group. We detect prominent MgII, FeII, CaII, and MgI absorption extending up to 62 kpc from a massive star-forming spiral galaxy and its interacting companion. Together with four other members, these form a compact group with a virial radius of 313 kpc. Down-the-barrel, blueshifted absorption indicates outflows. The distribution and two-dimensional kinematics of this gas suggest the influence of tidal stripping and star formation-driven winds. Intervening absorption across the field partly traces internal galaxy motions. A simple superposition of individual discs cannot reproduce the velocity field at large impact parameters or in counter-rotating regions, while a global IGrM halo with a rotational velocity of ~130 km/s provides a good match. Beyond individual galaxy envelopes, the data are consistent with a group-scale structure that co-rotates in concert with the galaxies. Assuming dynamical equilibrium, we estimate a total (cool+warm+hot) gas mass of 1.3-2.5x10^11 Msol, with large systematic uncertainties, corresponding to approximately 50% of all baryons, within one-quarter of the group's virial radius. These results point to a multiphase IGrM in which cool (~10^4 K) clouds are embedded within a dynamically coherent, group-wide halo. The gas appears gravitationally bound to the group rather than reaccreting onto individual galaxies.
In this work, we used numerical integration of the 4-body problem to study 3-body resonance chains (two planets and an asteroid in the innermost orbit) as a possible mechanism for white dwarf pollution. Two 3-body resonance chains were selected for study: the 6:3:2 and the 4:2:1. Asteroids in both a dynamically colder initial orbit in the 6:3:2 resonance and hotter initial orbits in both resonances were studied. An asteroid had up to a 1.08% chance of being delivered to the stellar Roche Limit of the white dwarf. This probability was strongly linearly correlated with the mass of the inner planet but was not correlated with the mass of the outer planet for both colder and hotter orbits. Average dynamical lifetimes ranged from 23 kyr to 1137 kyr for the dynamically colder orbit and from 12.9 kyr to 89.2 kyr and 10.8 kyr to 793.4 kyr for the dynamically hotter orbits in the 6:3:2 and 4:2:1 resonances, respectively. Average dynamical lifetime was exponentially anticorrelated with the outer planet mass and usually with the inner planet mass except in one case. The hotter 4:2:1 resonance delivered 1.1 times more asteroids to the stellar Roche Limit than the hotter 6:3:2 resonance. The hotter 6:3:2 resonance delivered 1.2 times more asteroids to the stellar Roche Limit than the colder 6:3:2 resonance. A typical accretion rate for a white dwarf star of 10^8 grams s^-1 could be explained by the accretion of an equivalent mass of one of our simulated asteroids every 13.8 Myr.
Dramatically increasing data volumes are forcing astronomers to adopt automated methods for the identification and classification of astronomical objects. Although deep-learning models are often well-suited to this task, obtaining a measure of uncertainty on their predictions is challenging. Here we consider the suitability of Monte Carlo conformal prediction (MCCP) set size and confidence as measures of model uncertainty for the astronomical classification of radio galaxies. We demonstrate this approach using model predictions from a pre-trained radio galaxy foundation model, fine-tuned on a smaller set of labelled radio galaxies. We calibrate the MCCP by obtaining annotator-derived soft label distributions, i.e. probability distributions over classes instead of single class assignments, for each of these labelled radio galaxies and compare the resulting set sizes and confidence scores to predictive entropy measures for each galaxy obtained using a supervised Bayesian deep-learning model trained using Hamiltonian Monte Carlo (HMC). The comparison reveals only a weak correlation between the measures.
Water and methanol are key components of interstellar ices and gas in star- and planet-forming regions, but direct observations of water in low-mass protostars are challenging due to atmospheric absorption. We present high-resolution (R = 70,500) mid-infrared spectroscopy of the Class I protostar SVS13-A with EXES on board SOFIA at 26 $μ$m, targeting both H$_2$O and CH$_3$OH absorption lines. Several lines of each species are detected, tracing warm gas with rotational temperatures of $\sim$140--170 K. Remarkably, the methanol column density is a factor of $\sim$4 higher than that of water, well above typical interstellar ice ratios ($<$10\%). Comparison with previous millimeter observations indicates that absorption and emission probe distinct regions, with the mid-IR lines likely tracing cooler gas along the line of sight. The surprising observed CH$_3$OH/H$_2$O ratio may reflect selective sublimation due to the distribution of binding energies or ice stratification in the inner envelope. These observations probe the inner regions of the protostar, where planets are expected to form and inherit the chemical composition of their natal environment, providing a direct link between ice sublimation and gas-phase chemistry. Our results represent the first high-spectral-resolution mid-infrared view of both water and methanol toward a low-mass protostar, offering a unique window into the chemical composition of the innermost envelope and planet-forming region, and highlighting the diagnostic power of high-resolution mid-infrared spectroscopy to uncover hidden chemical layers and the ice-to-gas transition in embedded protostars.
Observations of the binary system DX Cha (HD 104237) reveal a compact, asymmetric ring structure with a radius of 0.43\,au. This ring is just outside the binary orbit, which has semi-major axis $a_{\rm b} = 0.22$\,au and eccentricity $e_{\rm b} = 0.665$; placing the ring at $\approx 1.2$ times the binary apocenter distance. The inner regions of circumbinary disks, $\approx 2-3\,a_{\rm b}$, are typically evacuated by strong gravitational torques from the binary, resulting in a deep gap between the binary and the disk. Accordingly, previous numerical simulations of DX Cha have found an eccentric inner cavity with almost no material inside $\approx 1$\,au, and we find similar results when making the same assumption that the circumbinary disk orbits in the same direction as the binary. However, the disk can exist much closer to the binary if it is retrograde. For DX Cha we find that the inner edge of a retrograde disk occurs at $\approx 2a_{\rm b}$, and moreover takes the form of one or two arcs, in agreement with observations. We therefore suggest that the circumbinary disk in the DX Cha system could be orbiting retrograde to the binary star system in the center. We conclude that compact circumbinary disks observed in young stellar systems are important targets for future observations; if the disks are prograde then their properties are likely to be significantly different from current estimates, while if they are retrograde then this will have profound implications for our understanding of star and planet formation.
Active small bodies in extrasolar systems, the extrasolar analogues of Solar System comets, provide insights into the orbital evolution and physical processes shaping planetary systems. Since the discovery of exocomets around $β$ Pictoris, these small, icy bodies have shown the potential to become key probes for understanding planetary formation and migration. This review presents an overview of current observational techniques used to detect exocomets, focusing on individual systems and large-scale searches. We discuss photometric methods that identify exocomet transits through asymmetric light curves and spectroscopic techniques revealing cometary gases via time-variable absorption lines. Despite progress, significant open questions remain regarding the physical properties, occurrence rates, and similarities between exocomets and their Solar System counterparts. This review explores future opportunities in observational exocomet research, highlighting advancements required to further our understanding of these active small bodies and their role in the context of planetary system evolution.
Ground Level Enhancements (GLEs) probe the earliest, highest-energy solar energetic particles and thus provide a unique window onto particle release and transport from the low corona to 1 AU. We present a uniform, event-resolved analysis of the early anisotropy for ten well-observed GLEs, combining consistently reconstructed neutron-monitor pitch-angle distributions (PADs) with Parker-spiral footpoint mapping. We find a clear, monotonic decline of initial anisotropy with increasing magnetic connection angle: well-connected events exhibit strong, persistent forward-directed beams, while poorly connected events show systematically weaker and more rapidly decaying anisotropies. This relationship holds across a wide range of flare classes and CME speeds, demonstrating that magnetic connectivity and interplanetary transport, rather than eruption magnitude, dominate the directional properties of the earliest relativistic arrivals at Earth. A principal component analysis was applied to time-resolved spectral and angular parameters to separate source-driven changes from transport effects. By explicitly identifying and removing secondary sunward (back-scattered) components-attributable to scattering and reflection from solar-wind structures and transient interplanetary features-from the PAD fits, we isolate the intrinsic relaxation of the primary forward beam and show that apparent departures from simple exponential decay are frequently attributable to reflected or delayed populations rather than prolonged source injection. The empirical anisotropy--connection-angle relation reported here provides an event-resolved, quantitative benchmark that constrains focused-transport and shock-acceleration models and offers immediate operational value: rapid footpoint estimates can meaningfully limit expected initial beaming and directional radiation risk.
Hydrogen-rich supernovae can efficiently accelerate particles when the expanding ejecta interact with the surrounding circumstellar medium (CSM), producing high-energy (TeV--PeV) neutrinos. In this work we investigate the nearby SN~2024ggi, whose proximity and clear signatures of ejecta--CSM interaction make it a promising candidate for studying high-energy ($ν$) emission. We apply a new semi-analytical model that consistently links the electromagnetic and neutrino emission components, allowing us to constrain the main explosion parameters, including the kinetic energy, ejecta mass, progenitor radius, and nickel yield. The predicted high-energy ($ν$) fluence at Earth peaks at TeV energies and remains below the sensitivity of current detectors. However, the modeling establishes a robust framework for interpreting future signals from nearby interacting supernovae and fine-tuning observational strategies for next-generation multi-messenger facilities such as IceCube-Gen2 and KM3NeT/ARCA.
Measurements of the shear induced by weak gravitational lensing around galaxy cluster lines of sight are the gold standard for calibrating cluster observable-mass relations, thereby enabling a robust and precise inference of cosmological parameters. The weak-lensing mass bias is the systematic offset between the true halo mass and the mass that is inferred from the lensing data using an imperfect model for the halo mass distribution. We study the impact of cosmology on the lensing mass bias to inform future cosmological analyses of galaxy clusters. We create synthetic lensing shear maps for 115,920 projections of clusters with $M_{200\mathrm c}>1.56\times10^{14}\,h^{-1}M_\odot$ in a suite of Magneticum simulations. The simulation boxes are $896\,h^{-1}$Mpc on a side and are set up with 15 different combinations of the cosmological parameters $Ω_\mathrm{m}$, $Ω_\mathrm{b}$, $σ_8$, and $H_0$. Assuming a Navarro-Frenk-White profile, we extract weak-lensing mass measurements and quantify their bias $b_\mathrm{WL}$ with respect to the true halo mass. To investigate the impact of baryonic effects, we perform the analysis on gravity-only simulations and on their full-physics hydrodynamical counterparts. We confirm that assuming a fixed halo concentration or a fixed concentration-mass relation leads to cosmology-dependent changes of the mass bias. We report changes of up to $Δ\ln b_\mathrm{WL}=0.030$ with respect to the bias obtained at the fiducial WMAP7 cosmology. Adopting a model for the concentration that also depends on cosmology absorbs the changes in halo profiles and we recover essentially constant values for the mass bias. Our analysis of hydrodynamical simulations suggests that future, more accurate models will also need to explicitly account for the strength of baryonic effects.
BL Lacertae (BL Lac), the archetype blazar of its subclass and one of the most studied blazars in the last decades, has gone through a series of major multi-wavelength outbursts since 2020, resulting in its highest recorded $γ$-ray flare up to date between September and November 2022 together with those from August 2021 and October 2024. We characterise the $γ$-ray and multi-wavelength emission and spectral energy distribution (SED) of BL Lac, as well as their evolution during the major and extended $γ$-ray and multi-wavelength flare occurring between September and November 2022. We evaluate the variability of the flare, with focus on the nights of October 20 and November 13, when clear intranight very-high-energy (VHE, $E>100$ GeV) $γ$-ray variability is observed. We model the $γ$-ray and broadband SEDs during periods of stable emission identified with a Bayesian block analysis, interpreting their evolution of the flare from the variability of the relativistic particles and physical parameters of the jet. The VHE emission shows an average flux of 0.23 Crab Units (C.U.) above 200 GeV during this flare and a variability amplitude of more than a factor 10. Intranight doubling-flux variations as fast as $\sim$8 minutes are observed during the nights of October 20 and November 13, when maximum fluxes of 4.4 C.U. above 100 GeV and 2.8 C.U. above 200 GeV are reached. The spectral analysis reveals a transition of the X-ray emission from the high- to the low-energy SED peak, and a shift of the $γ$-ray peak towards higher energies. The broadband emission was interpreted within a leptonic two-zone model in which intranight variability is explained as magnetic reconnection in a compact region closely oriented with the line of sight while weekly-scale variations can be explained as variations of the electron distributions and the injection of accelerated particles.
It has been an unsolved question what leads a supernova remnant (SNR) to a thermal composite rather than a typical shell-like morphology, and what causes recombining plasma inside it. With the 13-ks observation of the Following-up X-ray Telescope onboard the Einstein Probe, we give an overall X-ray picture of W28, one of the prototypical thermal composite SNRs. The observation revealed a shell-like structure west of W28 in radio, optical, and X-ray images, which may revise the known extent of the SNR to $72'\times45'$. Spectral analysis explicitly maps that the special relationship where the plasma experiences recombination in the interior of the remnant, spatially coincident with H$α$ emissions, while in the other regions, the plasma is ionization-dominated. We found that W28 is generally isobaric from its center to the newly discovered shell, and it is even isothermal with a temperature of $\sim0.6$-0.7 keV in the center before the cooling of the plasma. Saturated thermal conduction and cloud evaporation may cool down the plasma within $\sim3$ kyr, the estimated recombination timescale. We revised the SNR dynamical age to $\sim8$ kyr, much younger than previous estimates. The complex structure and complex ionization state distribution may suggest that centrally filled and shell-like morphologies coexist in W28. This state may depend on the environment in which the SNR evolves.
Gravitational-wave observations of binary black hole (BH) mergers provide a novel avenue for testing massive-star evolution and the resulting BH mass spectrum. Recent population analyses under the hierarchical-merger hypothesis have offered evidence for the BH mass gap and inferred its lower edge to $\sim 44 - 68$ M$_\odot$. Motivated by these findings, we compute low-metallicity ($Z=10^{-5}$) helium star models with MESA and systematically explore the effect of uncertainties in the $^{12}$C$(α, γ)^{16}$O and $^{16}$O+$^{16}$O reaction rates on the final fate. Varying the $^{12}$C$(α, γ)^{16}$O reaction rate by $-3 σ$ to $+3σ$, we find that the predicted BH mass gap shifts from $\sim104 - 184$ M$_\odot$ to $\sim45 - 135$ M$_\odot$. In contrast, scaling the $^{16}$O+$^{16}$O reaction rate by global factors of 0.1, 1, and 10 has only a modest effect on the lower edge of the BH mass gap (less than 5 M$_\odot$), and shifts the upper edge by more than 10 M$_\odot$. Using the predictions of our models together with the literature estimates for the lower edge of the BH mass gap, we constrain the astrophysical S factor of $^{12}$C$(α, γ)^{16}$O reaction at 300 keV of $S_{300} \simeq$ 137.6 - 263.4 keV barn.
We aim to determine the most effective approach for estimating uncertainties in quasar photo-$z$ and to evaluate the ability of different models to reconstruct the true redshift distribution under varying data quality. We use photometric magnitudes from the Kilo-Degree Survey Data Release 5 and spectroscopically confirmed quasars from the Dark Energy Spectroscopic Instrument Data Release 1. We compare artificial neural networks (ANNs), Mixture Density Networks (MDNs), and Bayesian Neural Networks (BNNs), both latter combined with Gaussian Mixture Model (GMM) outputs. To assess robustness to observational limitations, we construct four test sets covering all combinations of sources fainter than those in the training sample and missing photometric bands. ANNs show substantial deviations in reconstructing the redshift distribution. MDNs require at least two Gaussian components to achieve accurate reconstruction, with the three-component MDN providing the best performance in this class. BNNs improve results for sources fainter than the training range, yielding a negative log-likelihood (NLL) gain of $0.11$, but reduce performance for brighter data by $0.07$ NLL. Reconstruction remains feasible for either fainter data or missing magnitudes individually; however, their combination leads to pronounced deviations. Unsupervised clustering identifies two dominant degenerate solutions at redshift pairs of $(1.2, 2.3)$ and $(1.6, 2.5)$. Accurate uncertainty modeling is essential for reliable reconstruction of the redshift distribution directly from photo-$z$. BNNs are particularly beneficial for out-of-distribution inference, although at the expense of reduced accuracy for brighter sources. Our methodology enables the identification and removal of degenerate photo-$z$ estimates unsuitable for tomographic analyses.