Papers for Monday, Mar 23 2026

We derive cosmological parameter constraints from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) galaxy clustering data, based on a joint full-shape analysis of the power spectrum multipoles and the bispectrum monopole using the ShapeFit framework. This is the follow-up of our previous work, in which we obtained for the first time constraints on the ShapeFit parameters using the bispectrum of DESI DR1. Here we present the first ShapeFit cosmological inference results using the bispectrum of DESI DR1. We recover values for the matter density parameter and Hubble constant of respectively $\Omega_m=0.310\pm0.012$ and $H_0=[68.92\pm0.97]\,\mathrm{km\, s^{-1} Mpc^{-1}}$, consistent with previous results from the full DESI DR1 dataset that did not use the bispectrum signal. The inclusion of the bispectrum significantly tightens the constraints on the amplitude of fluctuations, reducing the error-bars in $\ln(A_s\times10^{10})$ by approximately 20\%, compared to using the power spectrum alone. We also explore extended cosmological models by performing fits for the evolving dark energy equation of state $w_0w_a$, and the sum of neutrino masses $\sum m_\nu$. In these cases, we obtain constraints slightly larger than the ones from previous works from the DESI collaboration, due to not combining the full-shape results with other probes in all tracers. We find no strong evidence of deviations from standard $\Lambda$CDM, with the dark energy equation-of-state remaining within 2$\sigma$ from a cosmological constant $\Lambda$, and the neutrino mass being consistent with the normal hierarchy, $\sum m_\nu<0.1\,[eV]$ at 95\% confidence limit. These constraints are broadly consistent with other DESI DR1 analyses, thus validating the robustness of the ShapeFit compression approach and the inclusion of the bispectrum for cosmological inference.

Self-interacting dark matter (SIDM) models feature short-range interactions between dark matter (DM) particles that lead to larger diversity in the inner parts of galactic rotation curves and potentially unique gravitational lensing signatures. Satellite galaxies and dark subhalos provide a valuable testing ground for such models. We develop a simulation framework to explore subhalo evolution and its gravothermal collapse for velocity- and angle-dependent self-interacting cross section in these SIDM models. Our results are essential for testing these models. We perform N-body simulations, treating the host halo analytically and modelling the scattering-induced subhalo-halo interaction process using virtual host particles, a central innovation of our work. We use the Eddington inversion method to accurately model the local velocity distribution in the halo. Our approach is significantly less computationally expensive than simulations with a fully resolved host, while incorporating tidal stripping and tidal heating. We test both isotropic and forward-dominated self-scattering, which represent limiting cases for the angular dependence of the self-interaction cross section. Environmental effects, especially the scattering-induced subhalo-halo interaction, have a strong impact on the subhalo evolution and drive a complex structural evolution. As a result, SIDM subhalos have a larger range of central densities and density profile slopes compared to collisionless DM. Our cost-efficient simulation framework enables modelling of SIDM subhalos in realistic environments. Our results highlight the necessity of accurately modelling the scattering-induced subhalo-halo interaction to predict SIDM subhalo density profiles. For the SIDM models we investigate, the enhanced diversity in the mass profiles of subhalos would leave an observable imprint on strong lensing systems and satellite galaxies.

CONTEXT. One of the most difficult and unexplored regions of the Milky Way is the highly extincted in-plane central region within the Galactic coordinates $10^\circ \lesssim |\ell |\lesssim 30^\circ $, $|b|\lesssim 3^\circ $, where we have the long-bar and 3 kpc arm with intermediate-age stellar population, whose morphological properties are still unclear. AIMS. We aim to advance our knowledge of the morphology of these two components. METHODS. We examined star counts of bright M giants in WISE-4.6$\mu $m and its distribution of distances derived from spectroscopic parallaxes with APOGEE-DR17. We also examined the distribution of distances of young OGLE-O-rich Mira variable stars, and reviewed the literature on red clump distance determination within that area. RESULTS. We corroborate the asymmetry between positive and negative longitudes in in-plane regions, thus confirming the necessity to include a long bar. We obtain an average angle between the major axis of the long bar and the line Sun-Galactic centre of $\alpha =27.4^\circ \pm 1.5^\circ $, aligned with the triaxial bulge and a semi-major-axis length $\approx 4$ kpc. The tips of the long bar are in contact with the elliptical 3 kpc arm, with the major axis again aligned with the bulge and the long bar's major axes, whose tangential lines of sight correspond to $\ell =-22^\circ $ and $\ell=+27^\circ $. In the range of 50 degrees in the sky between these two longitudes, the stellar near 3 kpc arm is clearly detected at heliocentric distances around 5 kpc, and the stellar far 3 kpc arm is tentatively detected at heliocentric distances of 9-12 kpc.

We present a model for gas heating, driven by dynamical friction from orbiting subhalos within dark matter halos. Using data from the TNG50 simulation, we derive the subhalo mass function and calculate the dynamical friction heating rate for a wide range of halo masses and redshifts from $z = 15$ to 0. Our results show that, by converting gravitational potential energy into thermal energy, dynamical friction is an important mechanism for galaxy quenching in massive halos at low redshifts, consistent with previous studies. Additionally, we find that in the early universe at $z \sim 15$, heating rates can be comparable to the molecular hydrogen cooling rates in metal-free minihalos. This can suppress gas cooling and fragmentation and does increase the critical molecular fraction for Pop III star formation by up to one order of magnitude, thereby making Pop III star formation more difficult. In combination with the Lyman-Werner background, the dynamical friction heating mechanism favors the formation of direct-collapse black hole (DCBH) seeds in atomic cooling halos, even when the average H$_2$ fraction is $\sim 10^{-5}$ during the minihalo progenitor phase. Dynamical friction heating at a fixed host halo mass can vary by two orders of magnitude due to the scatter in the number of subhalos. To capture dynamical friction heating in simulations, it is necessary to resolve subhalos with a subhalo to host halo mass ratio $\psi \gtrsim 0.05$.

The detection of exoplanets using astrometry has long been an area of interest, but is fraught with challenges. The Gaia mission is fundamentally reshaping this field thanks to its unprecedentedly precise all-sky astrometric observations. The 2022 release of Gaia DR3 brought the first exoplanets discovered from the Gaia astrometry, including a new candidate around the bright ($V=6.6$) solar-type star HD 12800. However, two years after announcement, the Gaia exoplanet candidate was retracted. In this work we report radial velocity observations of HD 12800 acquired with the TRES spectrograph, which we began immediately after the release of Gaia DR3. Our observations failed to detect the planet candidate; nonetheless, we emphasise that the originally proposed companion would have been easily detected in our radial velocity observations. We conclude with a discussion on the role of intermediate-precision ($\approx$10 m s$^{-1}$) RV spectrographs in the follow-up of Gaia astrometric exoplanet candidates, relevant to the forthcoming release of Gaia DR4. We argue that such observations may play an important role in planet confirmation for stars between approximately $8<G<12$, likely to represent a significant fraction of Gaia exoplanet discoveries.

We apply novel developments in photoionization modeling and multi-frequency radiation hydrodynamics to the study of line driven AGN disc winds. We use a flux-averaged force multiplier approach to compute the radiation force due to lines for hydrodynamics simulations using 4 frequency bands - infrared (IR), optical (O), ultraviolet (UV) and X-rays. Though line driving is dominated by the UV, contributions from the O and X-ray bands are non-negligible and can lead to enhancements in the wind both in terms of mass flux and outflow velocity. Crucially, these effects are not captured when using a ``grey'' approach to the radiation modeling in the hydrodynamics, where frequency information is averaged over during the photoionization modeling. These results further strengthen the case for frequency dependent radiation dynamics studies for line driven winds.

Considerable observational evidence suggests that the activity of supermassive black holes in galactic nuclei is transient. The term ``active galactic nuclei archaeology'' has even been coined. This implies the possibility of reconstructing the history of activity, such as changes in the nuclear luminosity over time across various regions of the electromagnetic spectrum, by analysing how this activity manifested itself on galactic and extragalactic spatial scales. These phenomena include relic radio structures, gas clouds illuminated by the ``ionising echo'' of past activity, and Fermi/eROSITA bubbles. We provide a review of the results of galactic nucleus activity studies, focusing on its observable impact on the intergalactic medium and circumgalactic environment. Our main focus is on optical observations of ionisation cones and evidence of switching between radiative (ionisation cones) and kinetic (radio jets) modes of nuclear activity.

Vera C. Rubin Observatory's Prompt Processing system will automatically process 10 TB of raw images to produce up to 10 million transient alerts per night. We summarize how Prompt Processing meets its throughput, latency, and reliability requirements and present results from Rubin Observatory Commissioning.

Near-IR Paschen lines are potentially an excellent tracer of Type 1 AGNs that is hardly affected by dust extinction. JWST allows us, for the first time, to explore Paschen-line objects at redshift z>1. Here we present a study of 62 AGNs with broad Pa$\alpha$ and Pa$\beta$ lines at 1<z<3 using data from the JWST COSMOS-3D program. These AGNs are efficiently selected and identified using NIRCam imaging and grism slitless spectroscopic data. We separate the AGN-host emission with image decomposition and quantify dust attenuation with multi-band data. We construct a calibration sample with optical spectroscopy and use single-epoch, Mg II-based black hole masses ($M_{\mathrm{BH}}$) as an anchor to derive new, Paschen-based $M_{\mathrm{BH}}$ estimators. We obtain three sets of $M_{\mathrm{BH}}$ estimators based on Paschen line luminosities and AGN continuum luminosities at 1 and 2 $\mu$m, respectively. After dust corrections, they are well consistent with each other, and also broadly agree with previous results. With this AGN sample, we further construct the first Pa$\alpha$ and Pa$\beta$ luminosity functions (LFs) of Type 1 AGNs. The derived LFs are 3-5 times higher than those of UV/optical-selected AGNs, indicating that Paschen-selected Type 1 AGNs are more complete. In addition, the intrinsic properties of our AGNs show no dependence on dust reddening, suggesting that the observed reddening is unrelated to the central engine and is thus likely caused by line-of-sight obscuration.

Astronomical spectra, which encode rich astrophysical and chemical information, are fundamental to understanding celestial objects and universal laws. The advent of large-scale spectroscopic surveys, generating tens of millions of spectra, presents significant challenges for efficient data processing and analysis. To address these challenges, we develop an AI-powered platform (named ``SpecZoo'') for spectral visualization and analysis. This platform integrates modern information technology and machine learning to lower the barrier to spectral data utilization and enhance research efficiency. Its core functionalities include interactive visualization, automated spectral classification, physical parameter measurement, spectral annotation, and multi-band/multi-modal data fusion, all supported by flexible user and data management systems. It has become an essential tool for the National Astronomical Data Center, directly supporting spectral data processing and research for major projects including LAMOST, SDSS, DESI, and so on. Furthermore, the platform demonstrates strong potential for science-education integration, providing a novel resource for cultivating talent in astronomy and data science.

ATLAS J1138-5139 is a newly detected ultra-compact double white dwarf (DWD) system which is composed of a $1.02\,M_{\odot}$ carbon-oxygen white dwarf (CO WD) and a $0.24\,M_{\odot}$ helium (He) WD with an orbital period of about 27.68 min, making it one of the shortest-period DWD systems known. The future evolution and final fate of this system remain unexplored. In this work, we investigate the evolution of ATLAS J1138-5139 with the one-dimensional stellar evolution code Modules for Experiments in Stellar Astrophysics (MESA). We find that ATLAS J1138-5139 will evolve into an AM Canum Venaticorum (AM CVn) system in about \sim 6.3 Myr. Afterwards, the transferred material from the He WD companion start to build up to form a He shell near the surface of the CO WD. This accumulated He-shell masses can be up to approximately $0.12\,M_{\odot}$, which is likely to trigger a double-detonation (DDet) explosion of the CO WD. We therefore expect that ATLAS J1138-5139 will likely explode as a type Ia supernova eventually through the DDet explosion mechanism. Moreover, our calculations show that ATLAS J1138-5139 will be a promising target for gravitational-wave (GW) detection by future detectors like LISA, Tianqin and Taiji.

The Framework for Unified and Robust data Analysis with JAX (Furax) is an open-source Python framework for modeling data acquisition systems and solving inverse problems in astrophysics and cosmology. Built on JAX, Furax provides composable building blocks in the form of general-purpose and domain-specific linear operators, along with preconditioners and solvers for their numerical inversion. Domain-specific tools are provided for astrophysical and cosmic microwave background (CMB) data analysis$-$including map-making, instrument modeling, and astrophysical component separation$-$with a modular architecture designed to extend to other fields.

The Next-Generation Atmospheric Cherenkov Telescope Array (NG-ACTA) is proposed as a prospective infrastructure for very high energy (VHE) gamma-ray astronomy, consisting of a mixed-aperture array of 88 telescopes with a maximum array diameter of 10 km. The array adopts a three-tier configuration of 30 m large-aperture Large Size Telescopes (LSTs), 12 m medium-aperture Medium Size Telescopes (MSTs), and 6 m small-aperture Small Size Telescopes (SSTs), enabling continuous gamma-ray detection across the full energy band from 20 GeV to 100 TeV. With core advantages of an ultra-low detection threshold ($\leq20$ GeV), ultra-high angular resolution ($\leq0.04^\circ$), ultra-large effective area ($\geq1\times10^5$ m$^2$), extreme cosmic ray background rejection (proton rejection efficiency $\geq99.99\%$), and rapid transient response ($\leq100$ ns trigger latency), NG-ACTA targets the most cutting-edge and transformative fundamental scientific topics in modern astrophysics and particle physics, including VHE gamma-ray astronomy, cosmic ray origin, multi-messenger astronomy, and dark matter as well as new physics tests. The array's scientific goals cover five core fields: particle astrophysics, VHE gamma-ray astronomy, cosmic ray physics, multi-messenger astronomy, and new physics exploration, with six hierarchical and mutually supportive scientific objectives from Galactic to extragalactic sources, steady to transient objects, and conventional objects to dark matter. A comprehensive comparison with international under-construction facilities (e.g., CTAO-North, CTAO-South) and Chinese facilities (e.g., LACT) demonstrates that NG-ACTA leads the world in low-energy threshold, baseline length, background suppression, and multi-messenger rapid response capabilities.

Asteroseismology of member pulsators provides a robust physical constraint on cluster parameters by linking internal stellar structures to the global properties of the host cluster. However, the parameters of NGC 1647 remains poorly constrained due to limited investigation, a situation that cluster asteroseismology can significantly refine. In this study, we identified 271 high confidential cluster members in NGC 1647, using HDBSCAN clustering with radial-velocity validation. Its initial age is determined in the range of 1250-280 Myr, derived from isochrone fitting based on multi-survey metallicities extinction-corrected Gaia photometry. Among the members, we found 96 periodic variables from TESS and K2 photometry, including nine p-mode pulsators (five {\delta} Sct and four hybrid {\delta} Sct-{\gamma} Dor stars). Assuming a common cluster age and initial chemical composition, joint asteroseismic modeling is performed based on measured large frequency separations and individual mode frequencies. This yields a metallicity of [Fe/H] = -0.08+0.04-0.01, well consistent with the spectroscopic determinations, and a seismic age of 178+11-9 Myr, more precise than isochrone-based estimates. This work shows the diagnostic potential of {\delta} Sct asteroseismology in young open clusters and establishes a high-precision benchmark for future studies of NGC 1647 and other open clusters.

We investigate the morphological properties of large-scale structure in the Universe and the physical processes that modify the excursion-set morphology of the three-dimensional matter density field. Using the Quijote N-body simulation suite, we study how an initially Gaussian random matter density field is altered by non-linear gravitational evolution, redshift-space distortions, and massive neutrino free-streaming. To quantify these effects, we employ a comprehensive set of morphological descriptors, including Minkowski Functionals, Betti numbers, Minkowski Tensors, and local measures of the size and shape of connected components and cavities. We find that gravitational evolution, on quasi-linear scales $R_G \sim 10 h^{-1} \mathrm{Mpc}$, strongly skews the one-point distribution and slightly smooths the field via the merging of critical points, with a more pronounced effect for minima and wall saddle points than for peaks. Redshift-space distortions produce the strongest morphological signal, generating pronounced anisotropies that are robustly captured by Minkowski Tensors and local shape measures, arising from both coherent large-scale flows and non-linear Finger-of-God effects. In contrast, massive neutrinos induce an approximately isotropic suppression of small-scale structure, slightly reducing the amplitudes of the Minkowski Functionals while leaving individual shape measures largely unchanged. We further explore the sensitivity of these statistics to variations in cosmological parameters $\Omega_m$, $n_s$, and $\sigma_8$, finding that they probe strongly degenerate combinations of $\Omega_m$ and $n_s$, while also exhibiting sensitivity to $\sigma_8$ through the non-Gaussianity of the evolved density field.

To assess the velocity stability of CO spectral lines in the Milky Way Imaging Scroll Painting (MWISP) survey, we employ a cross-correlation method to measure velocity shifts across $\sim$ 10,000 CO spectra from six reference sources observed over the ten-year duration of the survey. The standard deviations ($\sigma$) of these measured velocity shifts range from 0.03 to 0.23 km s$^{-1}$ for their $^{12}$CO lines and 0.03 to 0.16 km s$^{-1}$ for $^{13}$CO lines. We find that larger shifts are associated with broader linewidths, more pronounced differences between monthly and long-term variations, and a stronger correlation between velocity shifts of $^{12}$CO and $^{13}$CO lines. By examining the relation of velocity shifts with the Azimuth-Elevation of the telescope, as well as the velocity fields of these extended sources, we find that the velocity shifts exhibit systematic changes across different Azimuth-Elevation ranges. The patterns and amplitudes of these changes vary among sources and are closely linked to the extended velocity fields of sources. This indicates that the increased velocity shifts are primarily caused by pointing errors, which are more sensitive to reference sources with higher velocity gradients. We also provide high signal-to-noise, velocity-aligned template spectra for reference sources.

Current pulsar timing models face challenges when applied to binary pulsars with wide orbits and low orbital eccentricities. The conventional \texttt{DD} model accurately characterizes the orbits of such systems, but it suffers from strong correlations between the time of periastron passage ($T_0$) and the longitude of periastron ($\omega$). The ELL1 model avoids these parameter correlations, yet fails due to the limitations of its first-order low-eccentricity approximation. Recent enhancements to the ELL1 model (dubbed ELL1+ model) have incorporated higher-order terms but retain the low-eccentricity approximation. In this study, we propose a further improved model, ELL1R, which eliminates reliance on the low-eccentricity approximation through rigorous calculation of the Römer delay. This modification can avoid strong parameter correlations in the DD model, and it can be used in systems with mild eccentricity $0.01\lesssim e\lesssim0.1$ where the ELL1+ model can not. Using the ELL1R model, we present the first phase-coherent timing solutions for three binary pulsars: PSR~J1851--0108 (orbital period: 228 days), PSR~J1910+0423 (886 days), and PSR~J1923+2022 (777 days). Validation against the DD and ELL1+ models confirms that ELL1R yields consistent timing results while integrating the advantages of the two models. Our analysis further indicates that all three pulsars are mildly recycled. The companions of PSRs J1910+0423 and J1923+2022 are likely white dwarfs, whereas the nature of PSR J1851--0108's companion remains unknown.

Material strength effects have been recently shown to be significant in giant impacts even at scales of planetary collisions. Despite this, their effects are often neglected in numerical giant impact simulations. We present an implementation of a basic strength model (pressure dependent shear strength) in the massively parallel smoothed particle hydrodynamics code pkdgrav3. The model includes elastic deviatoric stresses, plasticity with pressure-dependent yield strength, and thermal softening, and is fully integrated into the GPU-accelerated framework introduced in Paper I, preserving its scalability and performance characteristics. We validate the implementation against laboratory experiments of granular cliff collapse and our simulation results are in excellent agreement. We then determine the catastrophic disruption threshold, $Q_{RD}^*$, over a wide mass range of the colliding bodies using simulations performed both with and without material strength. Consistent with prior work, we find that strength substantially increases $Q_{RD}^*$ in the low-mass regime, while convergence toward the fluid limit occurs only near $R_{C1} \sim 10^7$ m ($\sim 0.7,M_\oplus$), well above the often assumed $\sim 100$ km size limit. Entropy production and remnant morphology likewise remain sensitive to rheology at intermediate masses. Performance measurements show that including strength introduces only modest computational overhead while maintaining favorable scaling, thereby enabling realistic solid mechanics in large-scale impact simulations.

We present the HI galaxy observation results of the FATHOMER (FAst neuTral HydrOgen intensity Mapping ExpeRiment), a pilot drift scan survey by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The survey comprises 28 hours of observations over 7 nights in 2021, covering a $60\, °^2$ sky area in the frequency range 1.05-1.45 GHz. The HI galaxies are identified using both a matched-filtering algorithm and the SoFiA source-finding pipeline, which yield consistent detections. We derive the velocity width ($W_{50}$), flux density, and HI mass for detected galaxies. A total of 702 galaxies are identified with HI mass above $10^{6.2}\,{M_\odot}$, signal-to-noise ratio greater than 5, and redshift $z < 0.09$. Among these, 331 are previously known from the ALFALFA survey. Of the newly detected sources, 9 have spectroscopic confirmation from SDSS, 285 are matched to SDSS or DESI photometric data, and 77 lack optical counterparts--possible candidates for dark or faint galaxies. Comparison with ALFALFA shows that FAST enables detection of galaxies at higher redshifts and with lower HI fluxes, despite the radio frequency interference (RFI) and partial data masking. A preliminary HI mass function analysis reveals a higher characteristic mass and steeper low-mass slope than ALFALFA, indicating FAST's enhanced sensitivity to massive and distant HI systems. These results demonstrate FAST's strong potential for future deep HI surveys and highlight the importance of improved RFI mitigation and completeness correction.

Recent transmission spectra of the temperate sub-Neptune K2-18 b obtained with JWST have attracted significant attention. Debates have quickly arisen over the interpretation of the spectral data, particularly the recent MIRI observation where dimethyl sulfide (DMS) and dimethyl disulfide (DMDS) are claimed. Here we revisit K2-18 b as a case study to examine several key questions that are also broadly relevant to the temperate sub-Neptune population: i) Can the low water abundance be reconciled with water clouds driven by orbital eccentricity? ii) Are the observed and non-observed atmospheric compositions mutually consistent? iii) Is it kinetically possible to produce DMS under sub-Neptune conditions? To address these questions, we couple climate and photochemical models to obtain self-consistent climate-photochemistry states for K2-18 b with a moderate orbital eccentricity of 0.2, as suggested by radial-velocity measurements. In addition, we present new laboratory measurements of DMS and DMDS infrared opacities by HFML-FELIX and compile updated C$_2$H$_6$ (ethane) opacities that include weak overtone bands. Our results support the interpretation of a sub-Neptune scenario without invoking DMS, and we do not find strong evidence for a water-rich interior.

Aims. In a first study, we characterised the properties of the gas component in the circumstellar envelopes surrounding a sample of 29 AGB stars with UV excesses. Now we intend to complement this information with an analysis of the dust component and compare the estimated parameters with those previously inferred from larger samples of AGB stars. Methods. We modelled the spectral energy distributions of the sample using dust radiative transfer models. In some cases, we complemented the analysis with Herschel/PACS radial surface brightness profiles. Results. We derived mass-loss rates and gas-to-dust ratios, which are in the typical ranges for AGB stars. We found that the stellar and mass-loss parameters follow similar trends than those presented in the literature. There is an anticorrelation between the gas-to-dust ratio and the UV emission, although it is weaker than its correlations with pulsation and mass-loss. We also estimated the dust attenuation produced by the dust at UV wavelengths and describe its effects on the intrinsic UV emission. Conclusions. Stellar and mass-loss parameters of UV emitting AGB stars follow similar trends as found for larger samples of AGB stars. High-angular resolution observations are required to explore the dust forming regions and identify the presence of stellar companions. Circumstellar dust attenuation might play a dominant role in the observed UV emission, and needs to be accounted to estimate the intrinsic UV emission.

Intensity mapping (IM) represents an innovative and potent probe to cosmology. One of its prime applications is to measure the Baryonic Acoustic Oscillations (BAO) in the late universe. We study the BAO measurement by IM in configuration space using simulations, focusing on the impact of the telescope beam and foreground removal effects. Three types of correlation functions are applied to measure BAO, including the radial correlation function, multipole correlation function, and wedge correlation function. We check our pipeline against a set of IM mock catalogs, finding good agreement with the numerical results. We use the mock catalogs to look for the parameter choices that optimize the BAO constraint for the correlation function estimators. With the optimal settings, our pipeline is utilized to forecast the BAO constraint for the 21 cm IM experiments: BINGO, MeerKAT, and SKA-mid. We find that for the low redshift experiments BINGO and MeerKAT, the wedge correlation function achieves the tightest constraint for both the transverse and radial BAO. For SKA-mid, the radial correlation function and wedge correlation function deliver the tightest constraint for the radial and transverse BAO, respectively.

Variability is one of the classic features of active galactic nuclei (AGNs). The normalized structure function was applied to distinguish variability samples from OVRO, ASAS-SN and Fermi. A power-law function model was selected to fit the structure functions of samples of three bands. We present the available samples of three bands, and by integrating two parameters, we obtain ideal discrimination results for three bands. Meanwhile, the differences between BL Lacs and FSRQs of Fermi and non-Fermi samples are well verified. The results show that the improved structure function can effectively distinguish samples of radio, optical, and gamma-ray. Additionally, BL Lacs and FSRQs in both Fermi and non-Fermi samples can be distinguished. The conclusion obtained through the distinction of structural functions in different bands supports that the variability in the three bands are caused by different physical mechanisms respectively: the samples in the optical band are radio quiet AGNs, and their variability is mainly caused by the fluctuations of the accretion disk, and the samples of radio band and gamma-ray band are radio loud AGNs whose variability is mainly caused by relativistic jet radiation. This conclusion conforms to the unified standard interpretation of variability about AGNs. Using these two parameters, we verify that there is no fundamental difference between Fermi and non-Fermi BL Lacs, while significant differences exist between FSRQs. However, the power exponent of the two can well distinguish BL Lacs.

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.

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$(\alpha, \gamma)^{16}$O and $^{16}$O+$^{16}$O reaction rates on the final fate. Varying the $^{12}$C$(\alpha, \gamma)^{16}$O reaction rate by $-3 \sigma$ to $+3\sigma$, 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$(\alpha, \gamma)^{16}$O reaction at 300 keV of $S_{300} \simeq$ 137.6 - 263.4 keV barn.

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 $\gamma$-ray flare up to date between September and November 2022 together with those from August 2021 and October 2024. We characterise the $\gamma$-ray and multi-wavelength emission and spectral energy distribution (SED) of BL Lac, as well as their evolution during the major and extended $\gamma$-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) $\gamma$-ray variability is observed. We model the $\gamma$-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 $\gamma$-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.

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 $\Omega_\mathrm{m}$, $\Omega_\mathrm{b}$, $\sigma_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 $\Delta\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.

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 ($\nu$) 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 ($\nu$) 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.

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.

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 $\beta$ 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.

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.

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 $\mu$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.

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.

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.

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.

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 {\mu}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 {\mu}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.

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.

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.

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.

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.

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 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.

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 $\Lambda$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.

Satellite galaxies in the Local Group tend to be distributed in thin, planar configurations, with many sharing coherent orbital motion. Galaxy formation simulations in $\Lambda$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.

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.

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_\phi/b_1$ -- the ratio of their responses to long-wavelength primordial potential $b_\phi$ 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_\phi/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.

Measurements of Waves instrument onboard the Juno spacecraft suggest that narrowband kilometric radiation (nKOM; 20-141 kHz) and narrowband low-frequency radiation (nLF; 5-70 kHz) are generated within the plasma near the Io plasma torus (IPT) in low-latitude regions. While these emissions are thought to result from the conversion of the natural modes of the plasma into escaping radio waves, at either the fundamental or the first harmonic of the plasma frequency, there is no consensus on the specific mechanism involved. Using the electron density and the magnetic field data from the Jovian Auroral Distribution Experiment (JADE) and the FluxGate Magnetometer (FGM), we determine the range of frequencies accessible to different wave modes during Juno's crossing of the plasma disk environment. We classify the observed nKOM and nLF according to their propagation modes: trapped (Z-mode or Whistler), escaping (X-mode or O-mode), and undetermined (either trapped or escaping). We apply the 3D numerical modeling method that was developed in Boudouma et al. (2024) to the escaping and undetermined nKOM and nLF observations, deriving macroscopic constraints on the generation mechanisms, wave modes, characteristic frequencies, beaming and source locations. Our results support the interpretation in which high-latitude nKOM is consistent with O-mode, while low-latitude is rather X-mode. Both nKOM and nLF appear to be generated near the fundamental of the plasma frequency, but only nLF shows compatibility with emission near the first harmonic, suggesting the possible coexistence of both linear and nonlinear generation mechanisms.

Cosmic-ray (CR) antinuclei, especially antideuteron $\overline{\rm D}$ and antihelium-3 nuclei ${}^3\overline{\rm He}$, are among the most promising messengers for indirect dark matter (DM) searches. This is because secondary production in CR interactions with the interstellar medium is strongly suppressed at kinetic energies $K\simeq (0.1 - 1)$ GeV/$n$, typically one to two orders of magnitude below fluxes expected in standard DM scenarios. From the theoretical side, the formation of $\overline{\rm D}$ and ${}^3\overline{\rm He}$ is governed by coalescence, whose dynamics cannot yet be reliably derived from first principles. Phenomenological approaches therefore introduce effective coalescence parameters, possibly dependent on collision energy and production environment (hadronic versus electroweak). In this work we show, for the first time, that a common set of physically motivated coalescence models can simultaneously reproduce collider data in two qualitatively different regimes: ALICE measurements of (anti)deuteron production in $pp$ collisions at $\sqrt{s}=(0.9 - 13)$ TeV and the ALEPH $\overline{\rm D}$ multiplicity in hadronic $Z$ decays at $\sqrt{s}=m_Z$. We test both simple event-by-event prescriptions based on a relative-momentum cutoff, finding a preferred coalescence scale $p_{\rm coal}\simeq 0.2$ GeV, and quantum-mechanical models in the Wigner formalism. In the latter, a Gaussian bound-state wavefunction gives a best-fit momentum width, corresponding to $\delta\simeq 1.7$ fm, while a parameter-free implementation using the Argonne $v_{18}$ wavefunction (constrained by proton-neutron scattering data) agrees with ALICE spectra at the $\sim 25\%$ level. Overall, our results support an approximately universal coalescence description across energies and production environments, strengthening the theoretical basis for interpreting upcoming CR antinuclei searches.

We study a recursive Penrose process and the energy extraction for the decay of electrically charged particles in a Reissner-Nordström black hole spacetime with anti-de Sitter (AdS) asymptotics, incorporating the backreaction on the black hole's mass and charge. A recursive process requires that the decay products are confined in a finite region so that the emitted particles bounce back for further decay. In AdS spacetimes, the confinement arises naturally. Outgoing particles encounter a turning point and are reflected. One may impose a mirror at finite radius, but in AdS, backreaction makes these two confinement methods equivalent. Let $Q_n$ be the black hole charge after $n$ decays, and define $n_{\rm c}$ as the index for which the black hole's charge is zero, $Q_{n_{\rm c}}=0$. For $n_{\rm c}$ integer the black hole's charge decreases and reaches exactly zero after a finite number of decays, terminating the process. However, the last particle turns back, and encountering zero charge, falls into the hole. The final state is a charged black hole whose charge equals the sum of the original black hole and the initial particle charges. For $n_{\rm c}$ noninteger, the black hole charge decreases and can be arbitrarily small, but is never zero. The last allowed decay occurs at $n=n_{c}^-$, where $n=n_{c}^-$ is the greatest integer less than $n_{\rm c}$. Any further decay invalidates the approximations, the particles would carry a charge comparable to the black hole mass, transforming the problem into a two-body problem. The would-be subsequent decay would violate cosmic censorship and the process terminates before any inconsistency arises. In the integer and noninteger cases, the system yields a finite energy gain. Backreaction ensures that the process extracts a finite amount of energy. No black hole bomb occurs, the system works at most as an energy factory.

We perform forecasts for the Hubble parameter H_0 and for the parameter Xi_0 that describes modified gravitational-wave propagation, using information from the binary neutron star (BNS) mass function, for Einstein Telescope (ET), taken either in the triangle or in the ``2L'' configuration, as well as for detector network made by ET together with a 40-km Cosmic Explorer (CE). We restrict ourselves to BNSs with a large signal-to-noise ratio, SNR>50, which still give O(10^3) events yr-1, and we perform a full joint cosmology-population Bayesian inference. We find that, for ET in isolation, the two ET configurations perform comparably, yielding uncertainties of 12% and 11% on H_0 for the triangular and 2L designs, respectively, and 18% uncertainty on Xi_0 in both cases. For networks including ET and CE, we can constrain H_0 and Xi_0 to precisions of 9% and 6%, respectively. These results should be taken as a very conservative estimate of third-generation detectors' capabilities as a consequence of the high SNR cut. We project the constraints on the Lambda CDM expansion history and find that ET alone (triangular and 2L configurations) achieves its best precision on H(z) at z=0.23 and z=0.28, yielding a 10% and 6% precision, respectively. When CE is added to the network, the precision improves to 4% and 3% at z=0.37 and z=0.38, respectively.

In this work, we investigate observable signatures of a magnetically charged Anti-de Sitter black hole in string-inspired Euler-Heisenberg theory. We analyze photon trajectories, the photon sphere, and the resulting black hole shadow. We derive the photon sphere and shadow radii and show that both deviate from the Schwarzschild and Schwarzschild-AdS cases. In particular, the radii decrease monotonically as the magnetic charge parameter $Q_m$ increases, indicating that magnetic charge modifies light propagation near the black hole. We also study neutral and charged particle motion and compute the corresponding epicyclic frequencies. Using the effective potential method, we obtain the specific energy and angular momentum for stable circular orbits and determine the innermost stable circular orbit (ISCO). The presence of $Q_m$ shifts the ISCO radius and alters the orbital structure. The radial, vertical, and orbital frequencies show clear deviations from the Schwarzschild case. Finally, we confront the model with twin-peak quasi-periodic oscillation (QPO) data from stellar-mass, intermediate-mass, and supermassive black hole candidates. A two-dimensional Delta chi-square analysis in the ($r$, $Q_m$) space shows that the best fit corresponds to $Q_m=0$, although finite values remain allowed within confidence levels. At the 1 sigma level, we obtain an upper bound $Q_m/M$ less than about 0.2. These results indicate that while magnetic charge produces measurable theoretical deviations, current QPO data place only moderate constraints on its magnitude.

We report new results obtained in our lattice studies of the $Sp(4)$ gauge theory coupled to two fundamental Dirac fermions. This theory provides a candidate for the dynamical origin of dark matter models within the strongly interacting massive particle paradigm. We employ Lüscher's formalism to analyse finite-volume energy levels and study the scattering amplitude of two pseudoscalar states in the spin-1 channel. We present our preliminary findings for a set of ensembles generated within a broad range of (Wilson) fermion masses.

This work investigates the universality of cusp formation in the shadows of compact objects. The emergence of cusps is accompanied by three interrelated phenomena: a topological charge transition, an equal-area law governing the self-intersecting structure, and universal critical scaling behavior. We demonstrate that, because these phenomena originate from the global morphology of the shadow, they are fundamentally independent of specific spacetime metric details and apply across diverse models. These features are systematically analyzed for the Kerr black hole endowed with a running Newton coupling. By extending our framework to rotating traversable wormholes, we confirm that the same universal behavior persists in more general compact objects. Our study uncovers the universality underlying cusp formation, offering a model-independent framework for identifying non-Kerr signatures in future black hole observations.

We investigate a concrete realization of the Dark Dimension scenario, where a single large extra dimension is set at sub-millimeter scales. In this framework, the Casimir energy of bulk fields accounts for the observed dark energy. Working in a 5-dimensional setup with the Standard Model confined to a 4-dimensional brane, we derive the effective action for the radion. We demonstrate that a minimal model comprising only gravity and three right-handed bulk neutrinos typically yields a negative radion potential. To realize a positive vacuum energy, we consider some extensions with additional bulk degrees of freedom. These extensions generate a sufficiently flat positive potential that allows the radion to behave as a quintessence field, evolving slowly at the sub-eV scale. Finally, we analyze the evolution of the dark-energy equation-of-state parameter and show that our model is consistent with recent DESI BAO measurements, including the distance ratios $D_H/r_d$ and $D_M/r_d$.

We investigate quarkyonic matter within a relativistic quark model by combining the dual quarkyonic picture with the quark-meson coupling (QMC) model. Using relativistic gaussian quark wavefunctions for the nucleon, we construct the quarkyonic QMC (QQMC) model and study the properties of symmetric nuclear matter and pure neutron matter. We find that the quark saturation density depends sensitively on the nucleon size parameter and that nuclear interactions quantitatively modify the high-density behavior of the equation of state (EoS) and the sound velocity. In particular, the QQMC model yields an earlier onset of quark saturation than the noninteracting gaussian quarkyonic (GQ) model, indicating that nuclear interactions enhance the stiffening of the EoS in the quarkyonic regime.

Mechanisms that dynamically relax the vacuum energy offer a concrete way to approach the cosmological constant problem, but because relaxation is not confined to the vacuum energy alone it can have consequences for the rest of low-energy physics. We explore this issue using the recently proposed 'yoga' relaxation models as an explicit framework and show how relaxation differentially suppresses 'slow' physics relative to a characteristic timescale set by the mass of the relaxon. It therefore need not alter e.g. Higgs & collider physics but can dramatically change how light scalar fields participate in cosmology. We revisit the QCD axion in this setting and show that the suppression of the axion's vacuum potential reshapes its behaviour on cosmological timescales while leaving fast, high-energy processes unaffected. The result is to alter the axion mass-coupling relation away from the standard QCD band, driving it into a regime already ruled out by observational constraints. In particular, suppression of the vacuum axion potential allows the QCD matter-induced potential to dominate even for matter densities relevant to cosmology and everyday matter, potentially driving the axion away from the CP-conserving minimum for QCD-motivated parameters. We conclude that conventional QCD axions are unlikely to remain viable in their standard form within vacuum-energy relaxation frameworks.

The complex cislunar dynamical environment poses challenges for spacecraft navigation and Space Domain Awareness (SDA) operations, where the knowledge of current and future spacecraft states is essential. Conventional Gaussian-based approaches for SDA degrade under the nonlinearities that manifest in this regime. To accurately model the underlying dynamics and characterize uncertainty, this work explores the Modified Generalized Equinoctial Orbital Elements under high-fidelity propagation for cislunar applications. The Henze-Zirkler test for multivariate normality is leveraged to evaluate uncertainty evolution across a range of orbits, demonstrating improved preservation of Gaussian behavior in cislunar space.