KICC papers

arXiv:2401.03712v2 Announce Type: replace Abstract: (Reduced) Using the recent cosmological high-resolution zoom-in simulations, NewHorizon and Galactica, in which the evolution of black hole spin is followed on the fly, we have tracked the cosmic history of a hundred of black holes (BHs) with a mass greater than 2x10^4 Ms. For each of them, we have studied the variations of the three dimensional angle (Psi) subtended between the BH spins and the angular momentum vectors of their host galaxies. The analysis of the individual evolution of the most massive BHs suggests that they are generally passing by three different regimes. First, for a short period after their birth, low mass BHs (10^5 Ms). In this case, the BH spins tend to be well aligned with the angular momentum of their host galaxy and this configuration is generally stable even though BH merger episodes can temporally induce misalignment. We have also derived the distributions of cos(Psi) at different redshifts and found that BHs and galaxy spins are generally aligned. Finally, based on a Monte Carlo method, we also predict statistics for the 2-d projected spin-orbit angles lambda. In particular, the distribution of lambda traces well the alignment tendency in the 3-d analysis. Such predictions provide an interesting background for future observational analyses.

arXiv:2307.12428v2 Announce Type: replace Abstract: The fragmentation of gas to form stars in molecular clouds is intrinsically linked to the turbulence within them. These internal motions are set at the birth of the cloud and may vary with galactic environment and as the cloud evolves. In this paper, we introduce a new suite of 15 high-resolution 3D molecular cloud simulations using the moving mesh code AREPO to investigate the role of different decaying turbulent modes (mixed, compressive and solenoidal) and virial ratios on the evolution of a $10^4\mathrm{M}_{\odot}$ molecular cloud. We find that diffuse regions maintain a strong relic of the initial turbulent mode, whereas the initial gravitational potential dominates dense regions. Solenoidal seeded models thus give rise to a diffuse cloud with filament-like morphology, and an excess of brown dwarf mass fragments. Compressive seeded models have an early onset of star-formation, centrally condensed morphologies and a higher accretion rate, along with overbound clouds. 3D filaments identified using DisPerSE and analyzed through a new Python toolkit we develop and make publicly available with this work called FIESTA, show no clear trend in lengths, masses and densities between initial turbulent modes. Overbound clouds, however, produce more filaments and thus have more mass in filaments. The hubs formed by converging filaments are found to favour star-formation, with surprisingly similar mass distributions independent of the number of filaments connecting the hub.

arXiv:2403.14618v1 Announce Type: new Abstract: The redshifted 21-cm signal from the Cosmic Dawn and Epoch of Reionization carries invaluable information about the cosmology and astrophysics of the early Universe. Analyzing the data from a sky-averaged 21-cm signal experiment typically involves navigating through an intricate parameter space to accurately address various factors such as foregrounds, beam uncertainties, ionospheric distortions, and receiver noise for the search of the cosmological 21-cm signal. The traditional likelihood-based sampling methods for modeling these effects could become computationally demanding for such highly complex models, which makes it infeasible to include physically motivated 21-cm signal models in the analysis. Moreover, the inference with these traditional methods is driven by the assumed functional form of the likelihood function. This work demonstrates how Simulation-Based Inference through Truncated Marginal Neural Ratio Estimation (TMNRE) can naturally handle these issues at a significantly reduced computational cost than the likelihood-based methods. We estimate the posterior distribution on our model parameters with TMNRE for simulated mock observations, composed of beam-weighted foregrounds, physically motivated 21-cm signal, and radiometric noise. We find that maximizing the information content by simultaneously analyzing the data from multiple time slices and antennas significantly improves the parameter constraints and leads to a better exploration of the cosmological signal. We discuss the application of TMNRE for the current configuration of the REACH experiment and demonstrate how it can be utilized to explore potential avenues for REACH. The method presented here can be easily generalized for any sky-averaged 21-cm signal experiment.

arXiv:2403.13948v1 Announce Type: new Abstract: Arp 220 is the nearest ULIRG; it shows evidence of 100 pc-scale molecular outflows likely connected with galaxy-scale outflows traced by ionised and neutral gas. The two highly obscured nuclei of Arp 220 are the site of intense star formation, with extreme star-formation rate surface densities (~ 10^3 Msun/yr/kpc2). Despite extensive investigations searching for AGN activity in the Arp 220 nuclei, direct evidence remains elusive. We present JWST/NIRSpec IFS observations covering the 0.9 - 5.1 um wavelength range of the innermost (5''x4'', i.e. 1.8x1.5 kpc) regions of Arp 220. The primary goal is to investigate the potential presence of AGN signatures in the nuclear regions by analysing the spectra extracted from circular apertures of radius 55 pc (0.15'') around each of the two nuclei. We identify ~ 70 ionised and ~ 50 molecular emission lines in the nuclear spectra of Arp 220; we use recombination line ratios to measure optical extinctions in the range AV ~ 11 - 14 mag. High ionisation lines are not detected, except the [Mg IV] line at 4.49 um which we interpret as due to shocks rather than to AGN ionisation. We identify broadening and multiple kinematic components in the HI and H2 lines caused by outflows and shocks, with velocities up to ~ 550 km/s. Significantly higher velocities (up to ~ 900 km/s) are detected in the off-nuclear regions; however, they do not conclusively represent evidence for AGN activity. Even with the unprecedented sensitivity of JWST/NIRSpec IFS, achieving an unambiguous identification or exclusion of the presence of an AGN in the Arp 220 system remains challenging, because of its extreme dust obscuration.

arXiv:2312.00899v2 Announce Type: replace Abstract: The SPT0311-58 system resides in a massive dark matter halo at z ~ 6.9. It hosts two dusty galaxies (E and W) with a combined star formation rate of ~3500 Msun/yr. Its surrounding field exhibits an overdensity of sub-mm sources, making it a candidate proto-cluster. We use spatially-resolved spectroscopy provided by the JWST/NIRSpec Integral Field Unit (IFU) to probe a field of view (FoV) ~ 17 x 17 kpc^2 around this object. These observations have revealed ten new galaxies at z ~ 6.9, characterised by dynamical masses spanning from ~10^9 to 10^10 Msun and a range in radial velocities of ~ 1500 km/s, in addition to the already known E and W galaxies. The implied large number density, and the wide spread in velocities, indicate that SPT0311-58 is at the core of a proto-cluster, immersed in a very massive dark matter halo of ~ 6 x 10^12 Msun. Hence, it represents the most massive proto-cluster ever found at the EoR. We also study the dynamical stage of the system and find that it likely is not fully virialised. The galaxies exhibit a great diversity of properties showing a range of evolutionary stages. We derive their ongoing Ha-based unobscured SFR, and find that its contribution to the total SF varies significantly across the galaxies in the system. Their ionization conditions range from those typical of field galaxies at similar redshift recently studied with JWST to those found in more evolved objects at lower z. The metallicity spans more than 0.8 dex across the FoV, reaching nearly solar values in some cases. The detailed IFU spectroscopy of the E galaxy reveals that it is actively assembling its stellar mass, showing sub-kpc inhomogeneities, and a metallicity gradient that can be explained by accretion of low metallicity gas from the IGM. The kinematic maps indicate departures from regular rotation, high turbulence, and a possible pre-collision minor merger. (Abridged)

arXiv:2311.04270v2 Announce Type: replace Abstract: We conduct a systematic search for high-redshift galaxy overdensities at $4.9 = 7.955$ and $\left = 8.222$ (representing densities around $\sim 6$ and $\sim 12$ times that of a random volume). We estimate the total halo mass of these large-scale structures to be $11.5 \leq \mathrm{log}_{10}\left(M_{\mathrm{halo}}/M_{\odot}\right) \leq 13.4$ using an empirical stellar mass to halo mass relation, which are likely underestimates as a result of incompleteness. These protocluster candidates are expected to evolve into massive galaxy clusters with $\mathrm{log}_{10}\left(M_{\mathrm{halo}}/M_{\odot}\right) \gtrsim 14$ by $z = 0$.

arXiv:2309.00041v2 Announce Type: replace Abstract: Low-mass dwarf galaxies are expected to reside within dark matter haloes that have a pristine, `cuspy' density profile within their stellar half-light radii. This is because they form too few stars to significantly drive dark matter heating through supernova-driven outflows. Here, we study such simulated faint systems ($10^4 \leq M_{\star} \leq 2\times 10^6 \, M_\mathrm{\odot}$) drawn from high-resolution (3 pc) cosmological simulations from the `Engineering Dwarf Galaxies at the Edge of galaxy formation' (EDGE) project. We confirm that these objects have steep and rising inner dark matter density profiles at $z=0$, little affected by galaxy formation effects. But five dwarf galaxies from the suite also showcase a detectable HI reservoir ($M_{\mathrm{HI}}\approx 10^{5}-10^{6} \, M_\mathrm{\odot}$), analogous to the observed population of faint, HI-bearing dwarf galaxies. These reservoirs exhibit episodes of ordered rotation, opening windows for rotation curve analysis. Within actively star-forming dwarfs, stellar feedback easily disrupts the tenuous HI discs ($v_{\phi} \approx 10\, \mathrm{km} \, \mathrm{s}^{-1}$), making rotation short-lived ($\ll 150 \, \mathrm{Myr}$) and more challenging to interpret for dark matter inferences. In contrast, we highlight a long-lived ($\geq 500 \, \mathrm{Myr}$) and easy-to-interpret HI rotation curve extending to $\approx 2\, r_{1/2, \text{3D}}$ in a quiescent dwarf, that has not formed new stars since $z=4$. This stable gas disc is supported by an oblate dark matter halo shape that drives high-angular momentum gas flows. Our results strongly motivate further searches for HI in rotation curves in the observed population of HI-bearing low-mass dwarfs, that provide a key regime to disentangle the respective roles of dark matter microphysics and galaxy formation effects in driving dark matter heating.

arXiv:2403.10072v1 Announce Type: new Abstract: Intervening metal absorbers in quasar spectra at $z > 6$ can be used as probes to study the chemical enrichment of the Universe during the Epoch of Reionization (EoR). This work presents the comoving line densities ($dn/dX$) of low ionisation absorbers, namely, Mg II (2796\r{A}), C II (1334\r{A}) and O I (1302\r{A}) across $2 0.03$\r{A}. For the first time, we present the comoving line densities of 131 weak ($W 5$, in agreement with previous works. In the context of a decline in metal enrichment of the Universe at $z > 5$, the overall evolution in the incidence rates of absorption systems can be explained by a weak - possibly soft fluctuating - UV background. Our results, thereby, provide evidence for a late reionization continuing to occur in metal-enriched and therefore, biased regions in the Universe.

arXiv:2403.08431v1 Announce Type: new Abstract: We present the first statistical investigation of spatially resolved emission-line properties in a sample of 63 low-mass galaxies at $4\leq z

arXiv:2403.08241v1 Announce Type: new Abstract: We measure the three-dimensional power spectrum (P3D) of the transmitted flux in the Lyman-$\alpha$ (Ly-$\alpha$) forest using the complete extended Baryon Oscillation Spectroscopic Survey data release 16 (eBOSS DR16). This sample consists of 205,012 quasar spectra in the redshift range 2

arXiv:2306.11801v2 Announce Type: replace Abstract: We investigate the incidence and properties of ionized gas outflows in a sample of 52 galaxies with stellar mass between $10^7$ M$_{\odot}$ and $10^9$ M$_{\odot}$ observed with ultra-deep JWST/NIRSpec MSA spectroscopy as part of the JWST Advanced Deep Extragalactic Survey (JADES). The high-spectral resolution (R2700) NIRSpec observations allowed us to identify for the first time the signature of outflows in the rest-frame optical nebular lines in low-mass galaxies at $z>3$. The incidence fraction of ionized outflows, traced by broad components, is about 25-40$\%$ depending on the intensity of the emission lines. The low incidence fraction might be due to both the sensitivity limit and the fact that outflows are not isotropic but have a limited opening angle which results in a detection only when this is directed toward our line of sight. Evidence for outflows increases slightly with stellar mass and star-formation rate. The median velocity and mass loading factor (i.e., the ratio between mass outflow rate and star formation rate) of the outflowing ionized gas are 350 km s$^{-1}$ and $\eta=2.0^{+1.6}_{-1.5}$, respectively. These are 1.5 and 100 times higher, respectively than the typical values observed in local dwarf galaxies. These outflows are able to escape the gravitational potential of the galaxy and enrich the circum-galactic medium and, potentially, the inter-galactic medium. Our results indicate that outflows can significantly impact the star formation activity in low-mass galaxies within the first 2 Gyr of the Universe.

arXiv:2403.09589v1 Announce Type: new Abstract: We introduce cosmocnc, a Python package for computing the number count likelihood of galaxy cluster catalogues in a fast, flexible and accurate way. cosmocnc offers three types of likelihoods: an unbinned, a binned, and an extreme value likelihood. It also supports the addition of stacked cluster data, which is modelled consistently with the cluster catalogue. The unbinned likelihood, which is the main focus of the code, can take an arbitrary number of mass observables as input and deal with several complexities in the data, such as variations in the properties of the cluster observable across the survey footprint, the possibility of different clusters having measurements for different combinations of mass observables, redshift measurement uncertainties, and the presence on unconfirmed detections in the catalogue. If there are more than one mass observables, the unbinned likelihood is computed with the backward convolutional approach, a novel approach that is first implemented in cosmocnc. After developing the likelihood formalism and describing its implementation, we validate the code with synthetic Simons-Observatory-like catalogues, finding excellent agreement between their properties and cosmocnc's predictions and obtaining constraints on cosmological and scaling relation parameters featuring negligible biases. cosmocnc is publicly available at github.com/inigozubeldia/cosmocnc.

arXiv:2403.09538v1 Announce Type: new Abstract: The detection of blue-shifted absorption lines likely associated with ionized Iron K-shell transitions in the X-ray spectra of many Active Galactic Nuclei (AGN) suggests the presence of a highly ionized gas outflowing with mildly relativistic velocities (0.03c-0.6c), named Ultra-Fast Outflow (UFO). Within the SUBWAYS project we characterized these winds starting from a sample of 22 radio-quiet quasars at 0.1

arXiv:2403.07871v1 Announce Type: new Abstract: We present the first fully simulation-based hierarchical analysis of the light curves of a population of low-redshift type Ia supernovae (SNae Ia). Our hardware-accelerated forward model, released in the Python package slicsim, includes stochastic variations of each SN's spectral flux distribution (based on the pre-trained BayeSN model), extinction from dust in the host and in the Milky Way, redshift, and realistic instrumental noise. By utilising truncated marginal neural ratio estimation (TMNRE), a neural network-enabled simulation-based inference technique, we implicitly marginalise over 4000 latent variables (for a set of $\approx 100$ SNae Ia) to efficiently infer SN Ia absolute magnitudes and host-galaxy dust properties at the population level while also constraining the parameters of individual objects. Amortisation of the inference procedure allows us to obtain coverage guarantees for our results through Bayesian validation and frequentist calibration. Furthermore, we show a detailed comparison to full likelihood-based inference, implemented through Hamiltonian Monte Carlo, on simulated data and then apply TMNRE to the light curves of 86 SNae Ia from the Carnegie Supernova Project, deriving marginal posteriors in excellent agreement with previous work. Given its ability to accommodate arbitrarily complex extensions to the forward model -- e.g. different populations based on host properties, redshift evolution, complicated photometric redshift estimates, selection effects, and non-Ia contamination -- without significant modifications to the inference procedure, TMNRE has the potential to become the tool of choice for cosmological parameter inference from future, large SN Ia samples.

arXiv:2403.07414v1 Announce Type: new Abstract: Age-dating and weighting stellar populations in galaxies at various cosmic epochs are essential steps to study galaxy formation through cosmic times. Evolutionary population synthesis models with different input physics are used towards this aim. In particular, the contribution from the thermally pulsing asymptotic-giant-branch (TP-AGB) stellar phase, which peaks for intermediate-age 0.6-2 Gyr systems, has been debated upon for decades. Here we report the detection of strong cool star signatures in the rest-frame near-infrared spectra of three young (~1 Gyr), massive (~10^10 Msun) quiescent galaxies at large look-back time, z=1-2, using JWST/NIRSpec. The co-existence of oxygen- and carbon-type absorption features, spectral edges and features from rare species such as Vanadium, and possibly Zirconium, reveal a strong contribution from TP-AGB stars. Population synthesis models with significant TP-AGB contribution reproduce the observations considerably better than those with weak TP-AGB, which are those commonly used. These findings call for revisions of published stellar population fitting results, pointing to lower masses and younger ages, with additional implications on cosmic dust production and chemical enrichment. These results will stimulate new generations of improved models informed by these and future observations.

arXiv:2403.05620v1 Announce Type: new Abstract: We introduce the Hawai`i Supernova Flows project and present summary statistics of the first 1218 astronomical transients observed, 669 of which are spectroscopically classified Type Ia Supernovae (SNe Ia). Our project is designed to obtain systematics-limited distances to SNe Ia while consuming minimal dedicated observational resources. This growing sample will provide increasing resolution into peculiar velocities as a function of position on the sky and redshift, allowing us to more accurately map the structure of dark matter. This can be used to derive cosmological parameters such as $\sigma_8$ and can be compared with large scale flow maps from other methods such as luminosity-line width or luminosity-velocity dispersion correlations in galaxies. Additionally, our photometry will provide a valuable test bed for analyses of SNe Ia incorporating near-infrared data. In this survey paper, we describe the methodology used to select targets, collect and reduce data, and calculate distances.

arXiv:2403.05506v1 Announce Type: new Abstract: The Near-infrared Spectrograph (NIRSpec) on the James Webb Space Telescope is uniquely suited to studying galaxies in the distant Universe with its combination of multi-object capabilities and sensitivity over a large range in wavelength (0.6-5.3 microns). Here we present the NIRSpec Wide survey, part of the NIRSpec Instrument Science Team's Guaranteed Time Observations, using NIRSpec's microshutter array to obtain spectra of more than 3200 galaxies at $z>1$ at both low- and high-resolution ($R\approx100$ and 2700) for a total of 105 hours. With 31 pointings covering $\approx$320 arcmin$^2$ across the five CANDELS fields with exquisite ancillary photometry from the Hubble Space Telescope, the NIRSpec Wide survey represents a fast and efficient way of using JWST to probe galaxies in the early Universe. Pointing centers are determined to maximize the observability of the rarest, high-value sources. Subsequently, the microshutter configurations are optimized to observe the maximum number of "census" galaxies with a selection function based primarily on HST/F160W magnitude, photometric/slitless grism redshift, and predicted \ha\ flux tracing the bulk of the galaxy population at cosmic noon ($z_{\rm med}=2.0$). We present details on the survey strategy, the target selection, an outline of the motivating science cases, and discuss upcoming public data releases to the community.

arXiv:2403.04901v1 Announce Type: new Abstract: We present the largest catalog to-date of star clusters and compact associations in nearby galaxies. We have performed a V-band-selected census of clusters across the 38 spiral galaxies of the PHANGS-HST Treasury Survey, and measured integrated, aperture-corrected NUV-U-B-V-I photometry. This work has resulted in uniform catalogs that contain $\sim$20,000 clusters and compact associations which have passed human inspection and morphological classification, and a larger sample of $\sim$100,000 classified by neural network models. Here, we report on the observed properties of these samples, and demonstrate that tremendous insight can be gained from just the observed properties of clusters, even in the absence of their transformation into physical quantities. In particular, we show the utility of the UBVI color-color diagram, and the three principal features revealed by the PHANGS-HST cluster sample: the young cluster locus, the middle-age plume, and the old globular cluster clump. We present an atlas of maps of the 2D spatial distribution of clusters and compact associations in the context of the molecular clouds from PHANGS-ALMA. We explore new ways of understanding this large dataset in a multi-scale context by bringing together once-separate techniques for the characterization of clusters (color-color diagrams and spatial distributions) and their parent galaxies (galaxy morphology and location relative to the galaxy main sequence). A companion paper presents the physical properties: ages, masses, and dust reddenings derived using improved spectral energy distribution (SED) fitting techniques.

arXiv:2303.17676v2 Announce Type: replace Abstract: Very metal-poor stars ($\rm[Fe/H] $ 16) and obtain 38,000/41,000 additional metal-poor candidates with purity 29\%/52\%, respectively. We make our metal-poor star catalogs publicly available, for further exploration of the metal-poor Milky Way.

arXiv:2403.03977v1 Announce Type: new Abstract: The study of gas-phase metallicity and its spatial distribution at high redshift is crucial to understand the processes that shaped the growth and evolution of galaxies in the early Universe. Here we study the spatially resolved metallicity in three systems at $z\sim6-8$, namely A2744-YD4, BDF-3299, and COSMOS24108, with JWST NIRSpec IFU low-resolution ($R$ $\sim$ 100) spectroscopic observations. These are among the highest-$z$ sources in which metallicity gradients have been probed so far. Each of these systems hosts several spatial components in the process of merging within a few kiloparsecs, identified from the rest-frame UV and optical stellar continuum and ionised gas emission line maps. The sources have heterogeneous properties, with stellar masses log($M_*/M_\odot) \sim 7.6-9.3$, star formation rates (SFRs) $\sim1-15$ $M_\odot$ yr$^{-1}$, and gas-phase metallicities 12+log(O/H) $\sim 7.7-8.3$, which exhibit a large scatter within each system. Their properties are generally consistent with those of the highest-redshift samples to date ($z\sim3-10$), though the sources in A2744-YD4 and COSMOS24108 are at the high end of the mass-metallicity relation (MZR) defined by the $z\sim3-10$ sources. Moreover, the targets in this work follow the predicted slope of the MZR at $z\sim 6-8$ from most cosmological simulations. The gas-phase metallicity gradients are consistent with being flat in the main sources of each system. Flat metallicity gradients are thought to arise from gas mixing processes on galaxy scales, such as mergers or galactic outflows and SN winds driven by intense stellar feedback, which wash out any gradient formed in the galaxy. The existence of flat gradients at $z\sim6-8$ sets also important constraints on cosmological simulations, whose predictions on the cosmic evolution of metallicity gradients often differ significantly, especially at high redshift.