Recent IoA Publications

arXiv:2509.13394v1 Announce Type: new Abstract: Precision cosmology with Type Ia supernovae (SNe Ia) requires robust quality control of large, heterogeneous datasets. Current data processing often relies on manual, subjective rejection of photometric data, a practice that is not scalable for forthcoming surveys like the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). We present a Bayesian framework that automates this step by integrating anomaly detection directly into the light curve fitting process. While the framework is model-agnostic and compatible with any spectral energy distribution (SED) model, here we demonstrate its application with the SALT3 model, implemented fully on GPU using \texttt{JAX-bandflux} for computational efficiency. Our method models the probability of each photometric measurement being anomalous with respect to the model, simultaneously inferring the physical supernova parameters and each data point's posterior probability of being a contaminant. Applying this framework to the Hawaii Supernova Flows dataset, we demonstrate its three principal capabilities: (i) robust mitigation of isolated outliers; (ii) automated identification and rejection of entirely corrupted bandpasses; and (iii) preservation of valid data by flagging only specific anomalous points within otherwise usable filters. We also find that contaminants are systematically brighter and bluer, which if uncorrected could bias subsequent cosmological inference.

arXiv:2507.02284v2 Announce Type: replace Abstract: The Milky Way Survey of the Dark Energy Spectroscopic Instrument (DESI) has so far observed three classical dwarf spheroidal galaxies (dSphs): Draco, Sextans and Ursa Minor. Based on the observed line-of-sight velocities and metallicities of their member stars, we apply the axisymmetric Jeans Anisotropic Multi-Gaussian Expansion modeling (JAM) approach to recover their inner dark matter distributions. In particular, both the traditional single-population Jeans model and the multiple population chemodynamical model are adopted. With the chemodynamical model, we divide member stars of each dSph into metal-rich and metal-poor populations. The metal-rich populations are more centrally concentrated and dynamically colder, featuring lower velocity dispersion profiles than the metal-poor populations. We find a diversity of the inner density slopes $\gamma$ of dark matter halos, with the best constraints by single-population or chemodynamical models consistent with each other. The inner density slopes are $0.71^{+0.34}_{-0.35}$, $0.26^{+0.22}_{-0.12}$ and $0.33^{+0.20}_{-0.16}$ for Draco, Sextans and Ursa Minor, respectively. We also present the measured astrophysical J and D factors of the three dSphs. Our results indicate that the study of the dark matter content of dSphs through stellar kinematics is still subject to uncertainties behind both the methodology and the observed data, through comparisons with previous measurements and data sets.

arXiv:2509.14156v1 Announce Type: new Abstract: In multiple-planet systems, gravitational interactions of exoplanets could lead to transit timing variations (TTVs), whose amplitude becomes significantly enhanced when planets are in or near mean-motion resonances (MMRs). In cases where both TTVs and radial velocity (RV) measurements are available, combined analysis can break degeneracies and provide robust planetary and system characterization, even detecting non-transiting planets. In this context, HIP 41378 hosts five confirmed transiting planets with periods ranging from 15 to over 542 days, providing a unique dynamical laboratory for investigating wide multi-planet systems analogous to the Solar System. In this study, we present an intensive space-based photometric follow-up of HIP 41378, combining 15 new CHEOPS observations with eight TESS sectors, alongside data from K2, Spitzer, HST, and HARPS. We dynamically modeled the TTVs and RV signals of the two inner sub-Neptunes via N-body integration. These planets, HIP 41378 b ($P_{b}$ = 15.57 days) and HIP 41378 c ($P_{c}$ = 31.71 days), are close to ($\Delta\sim1.8$ %) a 2:1 period commensurability. We report a clear detection of TTVs with amplitudes of 20 mins for planet b and greater than 3 hrs for planet c. We dynamically confirm the planetary nature of HIP 41378 g, a non-transiting planet with a period of about 64 days and a mass of about 7 $M_{\oplus}$, close to a 2:1 commensurability with planet c, suggesting a possible MMR chain in the inner system. Our precise determination of the masses, eccentricities, and radii of HIP 41378 b and c enabled us to investigate their possible volatile-rich compositions. Finally, by leveraging on the last TESS sectors we constrained the period of HIP 41378 d to three possible aliases ($P_{d} =$ 278, 371, and 1113 days) suggesting that the system could be placed in a double quasi resonant chain, highlighting its complex dynamical architecture.

arXiv:2509.13411v1 Announce Type: new Abstract: We present results from a deep, coordinated $XMM$-$Newton$ + $NuSTAR$ observation of the type 1 Seyfert PG 1426+015, a source of particular interest as the most massive reverberation-mapped black hole to date ($\log [M_{\rm{BH}}/M_{\odot}]$ = $9.01^{+0.11}_{-0.16}$). The high-resolution RGS data confirm the 'bare' nature of the source, showing no evidence for absorption beyond the Galactic column, while the broadband spectrum unambiguously reveals the presence of relativistic reflection from the innermost accretion disc (in the form of a relativistically broadened iron emission and associated Compton reflection hump) as well as confirming the presence of the strong soft excess reported previously. We explore whether relativistic reflection can successfully account for the soft excess along with the higher-energy reflection features, utilizing the two most-commonly used reflection codes (REFLIONX, XILLVER). Ultimately we find that both models are able to successfully reproduce the soft excess, though in the case of the XILLVER model this is contingent on reducing the strength of the O VIII line included in the model, as otherwise this feature prevents the model from reproducing the data. The reflection models that successfully reproduce the broadband data imply a relatively high density for the accretion disc of $\log [n_{\rm{e}} / \rm{cm}^{-3}] \sim 18$, consistent with the loose anti-correlation seen from other AGN in the $\log [n_{\rm{e}} / \rm{cm}^{-3}]$ vs $\log[m_{\rm{BH}} \dot{m}^2]$ plane, as well as a moderate-to-high black hole spin of $a^* \gtrsim 0.7$. This preliminary spin constraint is strongly dependent on the assumption that the soft excess is dominated by relativistic reflection.

arXiv:2505.23889v2 Announce Type: replace Abstract: Active galactic nuclei (AGNs) consist of a central supermassive black hole (SMBH) embedded in a region with both high gas and stellar densities: the gas is present as a thin accretion disc that fuels the central SMBH, while the stars form a dense, roughly isotropic nuclear star cluster. The binaries present in such a cluster could be considered naturally as triples, with the SMBH as a third object, and their dynamics also depend on the interaction with the gas-rich disc. In this paper, we study the evolution of such a binary on an inclined orbit with respect to the disc. The binary experiences both eccentricity excitation via the von Zeipel-Lidov-Kozai (ZLK) effect and drag forces from each time it penetrates the disc. We find that, as the outer orbital inclination decreases, the evolution of inner orbital separation can transition from a regime of gradual hardening to a regime of rapid softening. As such binaries grow wider, their minimum pericentre distances (during ZLK oscillations) decrease. We show that a simple geometric condition, modulated by the complex ZLK evolution, dictates whether a binary expands or contracts due to the interactions with the AGN disc. Our results suggest that the interaction with gas-rich accretion disc could enhance the rate of stellar mergers and formation of gravitational wave sources, as well as other transients. The treatment introduced here is general and could apply, with the proper modifications, to hierarchical triples in other gas-rich systems.

arXiv:2509.14206v1 Announce Type: new Abstract: As old stellar systems, globular clusters (GCs) are key fossil tracers of galaxy formation and interaction histories. This paper is part of the LEWIS project, an integral-field spectroscopic survey of ultra-diffuse galaxies (UDGs) in the Hydra I cluster. We use MUSE spectroscopy and new VIRCAM $H$-band imaging data to study the GC populations and dark matter content in four dwarf galaxies. We retrieved line-of-sight velocities for all sources in the observed MUSE fields. Since the spectroscopic measurements are limited to relatively bright sources, we developed a multi-band photometric procedure to identify additional GC candidates too faint for spectroscopic confirmation. GC candidates were selected using a combination of photometric properties and morphometric criteria. Additionally, the $H$-band observations were used to constrain the stellar masses of the studied galaxies. Based on the spectroscopic classification, we confirm one GC in UDG3, two in UDG7, and four in UDG11, while UDG9 has no spectroscopically confirmed bright GCs. We identify four intra-cluster GCs in the vicinity of UDG3 and UDG11, and one ultra-compact dwarf with a radial velocity only $\Delta v = -85 \pm 10\mathrm{km\ s^{-1}}$ relative to UDG7, suggesting it may be bound to it. Considering completeness corrections and accounting for possible contamination, from photometry we estimate that the number of GCs ranges between 0 and $\sim40$ for the investigated UDGs. Their specific frequencies suggest that three out of four UDGs are either GC-rich, similar to those in the Coma cluster, or belong to an intermediate population as seen in the Perseus cluster. Dark matter content estimates, inferred from GC counts and stellar mass, indicate that these galaxies are dark-matter dominated, with dynamical-to-stellar mass ratios of $M_{\mathrm{dyn}} / M_\star \sim 10-1000$.

arXiv:2509.13307v1 Announce Type: new Abstract: We demonstrate a GPU-accelerated nested sampling framework for efficient high-dimensional Bayesian inference in cosmology. Using JAX-based neural emulators and likelihoods for cosmic microwave background and cosmic shear analyses, our approach provides parameter constraints and direct calculation of Bayesian evidence. In the 39 dimensional $\Lambda$CDM vs $w_0w_a$ shear analysis, we produce Bayes Factors and a robust error bar in just 2 days on a single A100 GPU, without loss of accuracy. Where CPU-based nested sampling can now be outpaced by methods relying on MCMC sampling and decoupled evidence estimation, we demonstrate that with GPU acceleration nested sampling offers the necessary speed-up to put it on equal computational footing with these methods, especially where reliable model comparison is paramount. We put forward both nested and gradient-based sampling as useful tools for the modern cosmologist, where cutting-edge inference pipelines can yield orders of magnitude improvements in computation time.

arXiv:2509.12488v1 Announce Type: new Abstract: At high metallicity, a majority of massive stars have at least one close stellar companion. The evolution of such binaries is subject to strong interaction processes, heavily impacting the characteristics of their life-ending supernova and compact remnants. For the low-metallicity environments of high-redshift galaxies constraints on the multiplicity properties of massive stars over the separation range leading to binary interaction are crucially missing. Here we show that the presence of massive stars in close binaries is ubiquitous, even at low metallicity. Using the Very Large Telescope, we obtained multi-epoch radial velocity measurements of a representative sample of 139 massive O-type stars across the Small Magellanic Cloud, which has a metal content of about one fifth of the solar value. We find that 45% of them show radial velocity variations which demonstrate that they are members of close binary systems, and predominantly have orbital periods shorter than one year. Correcting for observational biases indicates that at least 70[+11:-6]% of the O stars in our sample are in close binaries, and that at least 68[+7:-8]% of all O stars interact with a companion star during their lifetime. We found no evidence supporting a statistically significant trend of the multiplicity properties with metallicity. Our results indicate that multiplicity and binary interactions govern the evolution of massive stars and determine their cosmic feedback and explosive fates.

arXiv:2509.13320v1 Announce Type: new Abstract: We present the first evolving interior structure model for sub-Neptunes that accounts for the miscibility between silicate magma and hydrogen. Silicate and hydrogen are miscible above $\sim 4000$K at pressures relevant to sub-Neptune interiors. Using the H$_2$-MgSiO$_3$ phase diagram, we self-consistently couple physics and chemistry to determine the radial extent of the fully miscible interior. Above this region lies the envelope, where hydrogen and silicates are immiscible and exist in both gaseous and melt phases. The binodal surface, representing a phase transition, provides a physically/chemically informed boundary between a planet's "interior" and "envelope". We find that young sub-Neptunes can store several tens of per cent of their hydrogen mass within their interiors. As the planet cools, its radius and the binodal surface contract, and the temperature at the binodal drops from $\sim 4000$K to $\sim 3000$K. Since the planet's interior stores hydrogen, its density is lower than that of pure-silicate. Gravitational contraction and thermal evolution lead to hydrogen exsolving from the interior into the envelope. This process slows planetary contraction compared to models without miscibility, potentially producing observable signatures in young sub-Neptune populations. At early times ($\sim 10$-$100$Myr), the high temperature at the binodal surface results in more silicate vapour in the envelope, increasing its mean molecular weight and enabling convection inhibition. After $\sim$Gyr of evolution, most hydrogen has exsolved, and the radii of miscible and immiscible models converge. However, the internal distribution of hydrogen and silicates remains distinct, with some hydrogen retained in the interior.

arXiv:2509.12319v1 Announce Type: new Abstract: MACS J0138-2155 is the only known cluster to strongly lens two supernovae (SNe), Requiem and Encore, from the same host galaxy at z=1.949. We present seven independent mass models of the galaxy cluster built using six software packages. By conducting a blind analysis (no exchanges of results between modeling teams), we quantified uncertainties due to modeling and software. Through HST, JWST and MUSE observations, we assembled high-quality data products, including eight "gold" lensed image systems consisting of 23 images with secure spectroscopic redshifts, and one "silver" system with a likely redshift value. Restricting to the gold images, we obtain overall consistent model predictions of the positions, magnifications and time delays of SN Encore and SN Requiem images, especially for models with $\chi^2 \leq 25$. We predict the appearance of the next images of SNe Encore and Requiem with a time delay of >~3000 days and of ~3700 to 4000 days, respectively, based on a fiducial cosmological model of $H_0 = 70 {\rm\ km\ s^{-1}\ Mpc^{-1}}$ and $\Omega_{\rm m} = 0.3$. We obtain relations between $H_0$ and the time delays of SNe Encore and Requiem. In particular, for $H_0 = 73 {\rm\ km\ s^{-1}\ Mpc^{-1}}$, the four lowest $\chi^2$ models predict SN Requiem to reappear in ~Apr-Dec 2026; for $H_0 = 67 {\rm\ km\ s^{-1}\ Mpc^{-1}}$, in ~Mar-Nov 2027. Using the newly measured time delay between the two detected images of SN Encore by Pierel et al. (submitted) and our mass models, we jointly infer $H_0 = {\rm 66.9^{+11.2}_{-8.1}\ km\ s^{-1}\ Mpc^{-1}}$, where the uncertainty is dominated by that of the time delay. The long delays of the next-appearing SN Requiem and SN Encore images provide excellent opportunities to measure $H_0$ with an uncertainty of 2-3%. Our mass models form the basis for cosmological inference from this unique lens cluster with two strongly lensed SNe. (Abridged)

arXiv:2509.12301v1 Announce Type: new Abstract: Multiply-imaged supernovae (SNe) provide a novel means of constraining the Hubble constant ($H_0$). Such measurements require a combination of precise models of the lensing mass distribution and an accurate estimate of the relative time delays between arrival of the multiple images. Only two multiply-imaged SNe, Refsdal and H0pe, have enabled measurements of $H_0$ thus far. Here we detail the third such measurement for SN Encore, a $z=1.95$ SNIa discovered in JWST/NIRCam imaging. We measure the time delay, perform simulations of additional microlensing and millilensing systematics, and combine with the mass models of Suyu et al. in a double-blind analysis to obtain our $H_0$ constraint. Our final time-delay measurement is $\Delta t_{1b,1a}=-39.8_{-3.3}^{+3.9}$ days, which is combined with seven lens models weighted by the likelihood of the observed multiple image positions for a result of $H_0=66.9_{-8.1}^{+11.2} \rm{km} \rm{s}^{-1}\rm{Mpc}^{-1}$. The uncertainty on this measurement could be improved significantly if template imaging is obtained. Remarkably, a sibling to SN Encore (SN "Requiem") was discovered in the same host galaxy, making the MACS J0138.0-2155 cluster the first system known to produce more than one observed multiply-imaged SN. SN Requiem has a fourth image that is expected to appear within a few years, providing an unprecedented decade-long baseline for time-delay cosmography and an opportunity for a high-precision joint estimate of $H_0$.

arXiv:2509.11175v1 Announce Type: new Abstract: The 21-cm line of hydrogen is the most promising probe of the Dark Ages and Cosmic Dawn. We combine hydrodynamical simulations with a large-scale grid in order to calculate the effect of non-linear structure formation on the large-scale 21-cm power spectrum, focusing on redshifts $z=20-40$. As the clumping effect arises from small-scale density fluctuations, it offers a unique opportunity to probe the standard cold dark matter model in a new regime and thus potentially investigate the properties of dark matter. To this end, we also study a warm dark matter $-$ like model with a Gaussian cutoff on a scale of 50 kpc. We find that clumping has a significant impact on the large-scale 21-cm power spectrum. For example, for the Dark Ages case at $z=30$ and wavenumber $k=0.05$ Mpc$^{-1}$, small-scale clustering enhances the 21-cm power spectrum by $13\%$. Once Lyman-$\alpha$ coupling kicks in due to the first stars, the 21-cm signal strengthens, and the effect of clumping grows; it suppresses the observable power spectrum at $z=20$ by a factor of two, while the cutoff model has less than half the clumping impact. The clumping effect is significantly higher than the sensitivity of the planned Square Kilometre Array (SKA) AA$^\star$ configuration, by up to a factor of 20 for standard cold dark matter, though detection will require separation from foregrounds and from astrophysical contributions to the 21-cm power spectrum.

arXiv:2509.11055v1 Announce Type: new Abstract: It is thought that the Universe went through an early period known as the Dark Ages, during which primeval density fluctuations grew to form the first luminous objects, marking the beginning of Cosmic Dawn around 100 million years after the Big Bang. The 21-cm line of hydrogen atoms is the most promising probe of these epochs, with extensive observational efforts underway. We combine hydrodynamical simulations with a large-scale grid in order to precisely calculate the effect of non-linear structure formation on the global (sky-averaged) 21-cm radio intensity. We show that it presents a potential opportunity to probe the properties of dark matter in a new regime, corresponding to a length-scale of only 150,000 light years and a mass-scale of 20 million Solar masses. This effect can in principle be detected unambiguously during the Dark Ages, where the weak signal requires an array of global signal antennae. During Cosmic Dawn, when stellar radiation boosts the signal, a single global antenna suffices, but the clumping effect must then be separated from the effect of the stars. Our findings open new avenues for testing the nature of dark matter as well as non-standard cosmological models.

arXiv:2509.10271v1 Announce Type: new Abstract: Galaxy evolution theories require co-evolution between accreting supermassive black holes (SMBH) and galaxies to explain many properties of the local galaxy population, yet observational evidence for the mechanisms driving this co-evolution is lacking. The recent star-formation histories of the host galaxies of accreting SMBHs (Active Galactic Nuclei, AGNs) can help constrain the processes that feed SMBHs and halt star formation in galaxies, but are difficult to obtain for the most luminous AGNs (quasars). We introduce Mean-Field Independent Component Analysis (MFICA) to decompose quasar spectra and obtain recent star formation histories of their host galaxies. Applying MFICA to quasar spectra from the Sloan Digital Sky Survey (SDSS) DR7 Quasar Catalogue in the redshift range $0.16 \leq z \leq 0.76$, we find that 53 per cent of quasar host galaxies are star-forming, 17 per cent lie in the green-valley, while only 5 per cent are quiescent. This contrasts with 14, 11, and 74 per cent of a mass-matched control sample that are star-forming, green-valley, and quiescent, respectively. We find that $\sim25$ per cent of quasars are hosted by post-starburst galaxies, an excess of $28\pm1$ compared to our control sample. While the heterogeneity of recent star formation histories implies multiple SMBH feeding mechanisms, the excess of post-starburst host galaxies demonstrates the link between accreting SMBHs and a recent starburst followed by rapid quenching. Given that massive post-starburst galaxies are predominantly caused by gas-rich major mergers, our results indicate that $30-50$ per cent of quasars originate from merger-induced starbursts.

arXiv:2509.10136v1 Announce Type: new Abstract: TOI-1203 is a bright (V=8.6) G3 V star known to host a transiting warm sub-Neptune on a 25.5 d orbit. Here we report on an intensive high-precision radial velocity and photometric follow-up campaign carried out with the HARPS spectrograph and the CHEOPS space telescope. We found that TOI-1203 has an enhancement of $\alpha$ elements relative to iron of [$\alpha$/Fe]=$0.21\pm0.04$. With an age of $\sim$12.5 Gyr, TOI-1203 belongs to the old, $\alpha$-element enhanced stellar population of the galactic thick disk. We spectroscopically confirmed the planetary nature of the 25.5 d sub-Neptune TOI-1203 d, measured its mass ($M_{d}=7.39\pm0.62~M_{\oplus}$) and refined its radius ($R_{d}=2.918_{-0.045}^{+0.046}~R_{\oplus}$). We discovered the presence of an additional transiting super-Earth on a 4.2 d orbit (TOI-1203 b) with a mass of $M_{b}=3.51_{-0.32}^{+0.33}~M_{\oplus}$ and a radius of $R_{b}=1.520_{-0.046}^{+0.045}~R_{\oplus}$. We also revealed the presence of two additional low-mass planets at 13.1 d and 204.6 d (TOI-1203 c and e), with minimum masses of $5.46_{-0.50}^{+0.51}~M_{\oplus}$ and $42.10_{-1.78}^{+1.83}~M_{\oplus}$. We found that the outer planet TOI-1203 e lies on an eccentric orbit with $e_{e}=0.152\pm0.029$. We performed a stability analysis of the system confirming that there are configurations consistent with the observed parameters that are dynamically stable over billion-year timescales. While analyzing the HARPS time series, we discovered that the FWHM of the HARPS cross-correlation function shows a significant long-period signal ($\sim$615 d) that has no counterpart in the radial velocity data or in the remaining HARPS ancillary time series. We significantly detected the same signal in the FWHM of the Th-Ar calibration lines used to compute the nightly wavelength solution, and attributed this systematic effect to a long-term variation of the HARPS instrumental profile.

arXiv:2509.10455v1 Announce Type: new Abstract: There is no consensus on how baryon feedback shapes the underlying matter distribution from either simulations or observations. We confront the uncertain landscape by jointly analyzing new measurements of the gas distribution around groups and clusters -- DESI+ACT kinetic Sunyaev-Zel'dovich (kSZ) effect profiles and eROSITA X-ray gas masses -- with mean halo masses characterized by galaxy-galaxy lensing. Across a wide range of halo masses ($M_{500}=10^{13-14}M_\odot$) and redshifts ($08 \sigma$ discrepant. The FLAMINGO simulation variant with the most gas expulsion, and therefore the most suppression of the matter power spectrum relative to a dark matter only simulation, provides a good description of how much gas is expelled and how far it extends; the enhanced gas depletion is achieved by more powerful but less frequent AGN outbursts. Joint kSZ, X-ray, and lensing measurements form a consistent picture of gas expulsion beyond several $R_{500}$, implying a more suppressed matter power spectrum than predicted by most recent simulations. Complementary observables and next-generation simulations are critical to understanding the physical mechanism behind this extreme gas expulsion and mapping its impact on the large-scale matter distribution.

arXiv:2503.15317v2 Announce Type: replace Abstract: Galaxy major mergers are a key pathway to trigger AGN. We present the first detection of major mergers in the Euclid Deep Fields and analyse their connection with AGN. We constructed a stellar-mass-complete ($M_*>10^{9.8}\,M_{\odot}$) sample of galaxies from the first quick data release (Q1), in the redshift range z=0.5-2. We selected AGN using X-ray data, optical spectroscopy, mid-infrared colours, and processing \IE observations with an image decomposition algorithm. We used CNNs trained on cosmological simulations to classify galaxies as mergers and non-mergers. We found a larger fraction of AGN in mergers compared to the non-merger controls for all AGN selections, with AGN excess factors ranging from 2 to 6. Likewise, a generally larger merger fraction ($f_{merg}$) is seen in active galaxies than in the non-active controls. We analysed $f_{merg}$ as a function of the AGN bolometric luminosity ($L_{bol}$) and the contribution of the point-source to the total galaxy light in the \IE-band ($f_{PSF}$) as a proxy for the relative AGN contribution fraction. We uncovered a rising $f_{merg}$, with increasing $f_{PSF}$ up to $f_{PSF}=0.55$, after which we observed a decreasing trend. We then derived the point-source luminosity ($L_{PSF}$) and showed that $f_{merg}$ monotonically increases as a function of $L_{PSF}$ at z<0.9, with $f_{merg}>$50% for $L_{PSF}>2\,10^{43}$ erg/s. At z>0.9, $f_{merg}$ rises as a function of $L_{PSF}$, though mergers do not dominate until $L_{PSF}=10^{45}$ erg/s. For X-ray and spectroscopic AGN, we computed $L_{bol}$, which has a positive correlation with $f_{merg}$ for X-ray AGN, while shows a less pronounced trend for spectroscopic AGN due to the smaller sample size. At $L_{bol}>10^{45}$ erg/s, AGN mostly reside in mergers. We concluded that mergers are strongly linked to the most powerful, dust-obscured AGN, associated with rapid supermassive black hole growth.

arXiv:2508.00774v2 Announce Type: replace Abstract: We present a novel formation channel for supermassive black hole (SMBH) binaries in the early Universe, driven by primordial black holes (PBHs). Using high-resolution hydrodynamical simulations, we explore the role of massive PBHs ($m_{BH} \sim 10^6 M_\odot$) in catalyzing the formation of direct-collapse black holes (DCBHs), providing a natural in situ pathway for binary SMBH formation. PBHs enhance local overdensities, accelerate structure formation, and exert thermal feedback on the surrounding medium via accretion. Lyman-Werner (LW) radiation from accreting PBHs suppresses H$2$ cooling, shifting the dominant gas coolant to atomic hydrogen. When combined with significant baryon-dark matter streaming velocities ($v_{b\chi} \gtrsim 0.8 \sigma_{b\chi}$, where $\sigma_{b\chi}$ is the root-mean-square streaming velocity), these effects facilitate the formation of dense, gravitationally unstable, atomically cooling gas clouds in the PBH's wake. These clouds exhibit sustained high inflow rates ($\dot{M}_{infall} \gtrsim 0.01 - 0.1 M_\odot yr^{-1}$), providing ideal conditions for DCBH formation from rapidly growing supermassive stars of $\sim 10^5 M_\odot$ at redshifts $z \sim 20 - 10$. The resulting systems form SMBH binaries with initial mass ratios $q \sim O(0.1)$ and separations of $\sim 10$ pc. Such PBH-DCBH binaries provide testable predictions for JWST and ALMA, potentially explaining select high-$z$ sources such as the Little Red Dots (LRDs), and represent gravitational-wave sources for future missions like LISA and TianQin-bridging early-Universe black hole physics, multi-messenger astronomy, and dark matter theory.

arXiv:2507.12116v2 Announce Type: replace Abstract: We present two extensive sets of 3500+1000 simulations of dark matter haloes on the past light cone, and two corresponding sets of simulated (`mock') galaxy catalogues that represent the Euclid spectroscopic sample. The simulations were produced with the latest version of the PINOCCHIO code, and provide the largest, public set of simulated skies. Mock galaxy catalogues were obtained by populating haloes with galaxies using an halo occupation distribution (HOD) model extracted from the Flagship galaxy catalogue provided by Euclid Collaboration. The Geppetto set of 3500 simulated skies was obtained by tiling a 1.2 Gpc/h box to cover a light-cone whose sky footprint is a circle of 30 deg radius, for an area of 2763 deg$^2$ and a minimum halo mass of $1.5\times10^{11}$ Msun/h. The relatively small box size makes this set unfit for measuring very large scales. The EuclidLargeBox set consists of 1000 simulations of 3.38 Gpc/h, with the same mass resolution and a footprint that covers half of the sky, excluding the Milky Way zone of avoidance. From this we produced a set of 1000 EuclidLargeMocks on the 30 deg radius footprint, whose comoving volume is fully contained in the simulation box. We validated the two sets of catalogues by analysing number densities, power spectra, and 2-point correlation functions, showing that the Flagship spectroscopic catalogue is consistent with being one of the realisations of the simulated sets, although we noticed small deviations limited to the quadrupole at k>0.2 h/Mpc. We show cosmological parameter inference from these catalogues and demonstrate that using one realisation of EuclidLargeMocks in place of the Flagship mock produces the same posteriors, to within the expected shift given by sample variance. These simulated skies will be used for the galaxy clustering analysis of Euclid's Data Release 1 (DR1).

arXiv:2410.21365v2 Announce Type: replace Abstract: The stellar halo of the Milky Way comprises an abundance of chemical signatures from accretion events and \textit{in-situ} evolution, that form an interweaving tapestry in kinematic space. To untangle this, we consider the mixtures of chemical information, in a given region of integral of motion space, as a variant of the blind source separation problem using non-negative matrix factorisation (NMF). Specifically, we examine the variation in [Fe/H], [Mg/Fe], and [Al/Fe] distributions of APOGEE DR17 stars across the $(E,L_z)$ plane of the halo. When 2 components are prescribed, the NMF algorithm splits stellar halo into low- and high-energy components in the $(E,L_z)$ plane which approximately correspond to the accreted and \textit{in-situ} halo respectively. We use these components to define a boundary between the \textit{in-situ} and the accreted stellar halo, and calculate their fractional contribution to the stellar halo as a function of energy, galactocentric spherical radius ($r$), height ($z$), and galactocentric cylindrical radius ($R$). Using a stellar halo defined by kinematic cuts, we derive a boundary in $(E,L_z)$ space where the halo transitions from \textit{in-situ} dominated to accretion dominated. Spatially, we find that this transition happens at $(r,z,R) \approx (8.7, 3.0, 8.1)$ kpc. We find that between 34\% to 53\% of the stellar halo's content is of accreted origin. Upon prescribing more components to the NMF model, we find evidence for overlapping chemical evolution sequences. We examine features within these components that resemble known substructures in the halo, such as \textit{Eos} and \textit{Aurora}.