Search

We present the Lyman-$\alpha$ (Ly$\alpha$) and ultraviolet (UV) luminosity function (LF), in bins of redshift, of quasars selected in the Physics of the Accelerating Universe Survey (PAUS). A sample of 915 objects was selected at $2.710^{44}$ ${\rm erg\,s}$$^{-1}$. In order to obtain corrections for the LF estimation, and assess the accuracy of our selection method, we produced mock catalogs of $0

Being the most prominent HI line, Ly$\alpha$ permeates the cosmic web in emission. Despite its potential as a cosmological probe, its detection on large scales remains elusive. We present a new methodology to perform Ly$\alpha$ intensity mapping with broad-band optical images, by cross-correlating them with Ly$\alpha$ forest data using a custom one-parameter estimator. We also develop an analytical large-scale Ly$\alpha$ emission model with two parameters (average luminosity $\langle L_{\rm Ly\alpha} \rangle$ and bias $b_{\rm e}$) that respects observational constraints from QSO luminosity functions. We compute a forecast for DECaLS/BASS $g$-band images cross-correlated with DESI Ly$\alpha$ forest data, setting guidelines for reducing images into Ly$\alpha$ intensity maps. Given the transversal scales of our cross-correlation (26.4 arcmin, $\sim$33 cMpc/h), our study effectively integrates Ly$\alpha$ emission over all the cosmic volume inside the DESI footprint at $2.2 < z < 3.4$ (the $g$-band Ly$\alpha$ redshift range). Over the parameter space ($\langle L_{\rm Ly\alpha} \rangle$, $b_{\rm e}$) sampled by our forecast, we find a 3$\sigma$ of large-scale structure in Ly$\alpha$ likely, with a probability of detection of 23.95\% for DESI-DECaLS/BASS, and 54.93\% for a hypothetical DESI phase II with twice as much Ly$\alpha$ QSOs. Without a detection, we derive upper bounds on $\langle L_{\rm Ly\alpha} \rangle$ competitive with optimistic literature estimates ($2.3 \pm 1 \cdot 10^{\rm 41}$ erg/s/cMpc$^3$ for DESI, and $\sim$35\% lower for its hypothetical phase II). Extrapolation to the DESI-Rubin overlap shows that a detection of large-scale structure with Ly$\alpha$ intensity mapping using next-generation imaging surveys is certain. [abridged], Comment: 26 pages, 24 figures, submitted to MNRAS. Comments welcome!

Filaments connecting halos are a long-standing prediction of cold dark matter theories. We present a novel detection of the cosmic web emission connecting two massive quasar-host galaxies at cosmic noon in the MUSE Ultra Deep Field (MUDF) using unprecedentedly deep observations that unlock a high-definition view of the filament morphology, a measure of the transition radius between the intergalactic and circumgalactic medium, and the characterization of the surface brightness profiles along the filament and in the transverse direction. Through systematic comparisons with simulations, we validate the filaments' typical density predicted in the current cold dark matter model. Our analysis of the MUDF field, an excellent laboratory for quantitatively studying filaments in emission, opens a new avenue to understanding the cosmic web that, being a fundamental prediction of cosmology, bears key information on the essence of dark matter., Comment: submitted, comments welcome!

A population of non-stellar black holes ($\gtrsim$100 M$_{\odot}$) has been long predicted to wander the Milky Way. We aim to characterize this population by using the L-Galaxies semi-analytical model applied on top of the high resolution Millennium-II merger trees. Our results predict $\sim$10 wandering black holes with masses $\sim$2 $\times$ 10$^{3}$ M$_{\odot}$ in a typical $z$ = 0 Milky Way galaxy, accounting for $\sim$2$\%$ of the total non-stellar black hole mass budget of the galaxy. We find that the locations of these wanderers correlate with their formation scenario. While the ones concentrated at $\lesssim$1 kpc from the galactic nucleus on the disk come from past galactic mergers, the ones formed as a consequence of ejections due to gravitational recoils or the disruption of satellite galaxies are typically located at $\gtrsim$100 kpc. Such small and large distances might explain the absence of strong observational evidence for wandering black holes in the Milky Way. Our results also indicate that $\sim$67$\%$ of the wandering population is conformed by the leftovers of black hole seeds that had little to no growth since their formation. We find that wandering black holes that are leftover seeds become wanderers at an earlier time with respect to grown seeds, and also come from more metal-poor galaxies. Finally, we show that the number of wandering black holes in a Milky Way-type galaxy depends on the seeding efficiency., Comment: Accepted for publication in MNRAS, 14 pages, 9 figures

Pulsar Timing Array (PTA) experiments worldwide recently reported evidence of a nHz stochastic gravitational wave background (sGWB) compatible with the existence of slowly inspiralling massive black hole (MBH) binaries (MBHBs). The shape of the signal contains valuable information about the evolution of $z<1$ MBHs above $\rm 10^8 M_{\odot}$, suggesting a faster dynamical evolution of MBHBs towards the gravitational-wave-driven inspiral or a larger MBH growth than usually assumed. In this work, we investigate if the nHz sGWB could also provide constraints on the population of merging lower-mass MBHBs ($\rm {<} 10^7 \, M_{\odot}$) detectable by LISA. To this end, we use the $\texttt{L-Galaxies}$ semi-analytical model applied to the $\texttt{Millennium}$ suite of simulations. We generate a population of MBHs compatible simultaneously with current electromagnetic and nHz sGWB constraints by including the possibility that, in favourable environments, MBHs can accrete gas beyond the Eddington limit. The predictions of the model show that the global (integrated up to high-$z$) LISA detection rate is {\it not} significantly affected when compared to a fiducial model whose nHz sGWB signal is ${\sim}\,2$ times smaller. In both cases, the global rate yields ${\sim}\,12 \rm yr^{-1}$ and is dominated by systems of $\rm 10^{5-6} M_{\odot}$. The main differences are limited to low-$z$ ($z<3$), high-mass (${>}\rm 10^6\, M_{\odot}$) LISA MBHBs. The model compatible with the latest PTA results predicts up to ${\sim}\,1.6$ times more detections, with a rate of ${\sim}1\rm yr^{-1}$. We find that these LISA MBHB systems have 50\% probability of shining with bolometric luminosities $>10^{43}\rm erg/s$. Hence, in case PTA results are confirmed and given the current MBH modelling, our findings suggest there will be higher chances to perform multimessenger studies with LISA MBHB than previously expected., Comment: 16 pages; Submitted to A&A

We study a sample of 14 spectroscopically confirmed Ly$\alpha$ Emitters (LAEs) in the late era of reionization (at redshift $z\approx6$) based on the JWST/NIRCam imaging dataset. These LAEs with high Ly$\alpha$ luminosity of $L$(Ly$\alpha$) $\sim10^{42.4-43.4}$ erg s$^{-1}$ have been covered by the (ongoing) COSMOS-Web survey (Kartaltepe et al. 2021; Casey et al. 2022) over $0.28$ deg$^2$ in four NIRCam bands (F115W, F150W, F277W, and F444W). With JWST imaging, we determine the UV continua with $M_{\rm UV}$ ranging from ${-}20.5$ to ${-}18.5$ mag. The UV slopes have a median value of $\beta \approx-2.35$, and the steepest slope can reach $\beta<-3$. Under an excellent spatial resolution of JWST, we identify three out of the sample as potential merging/interacting systems. The 14 LAEs (and their components) are compact in morphology residing substantially below the mass-size relation of high-$z$ galaxies. We further investigate the stellar mass ($M_*$) and star-formation rates (SFRs). Most of the LAEs lie on the SFR-$M_*$ main-sequence relation while two of them featured as "little red dots" likely host active galactic nuclei (AGN), implying a ${\sim}10\%$ AGN fraction. Moreover, we reveal that a new correlation may exist between Ly$\alpha$ equivalent width and the offset between Ly$\alpha$ and UV emission ($\Delta d_{\rm Ly\alpha}$) with a median $\Delta d_{\rm Ly\alpha} \sim 1$ kpc. This could be explained by Ly$\alpha$ radiative transfer process in both ISM and CGM. The results usher a new era of detailed analysis on high-$z$ LAEs with the JWST capability., Comment: 15 pages, 5 figures, accepted by ApJL

We present the Lyman-$a$ (Lya) Luminosity Function (LF) at $2.0543.5. This work serves as a proof-of-concept for the results that can be obtained with the upcoming data releases of the J-PAS survey., Comment: 25 pages, 15 figures, Accepted for publication in A&A

The gravitational wave (GW) antenna LISA will detect the signal from coalescing massive black hole binaries (MBHBs) of $\rm 10^4\,{-}\,10^7\, M_{\odot}$, providing clues on their formation and growth along cosmic history. Some of these events will be localized with a precision of several to less than a deg$^2$, enabling the possible identification of their host galaxy. This work explores the properties of the host galaxies of LISA MBHBs below $z\,{\lesssim}\,3$. We generate a simulated lightcone by using the semi-analytical model $\mathrm{\texttt{L-Galaxies}}$ applied on the merger trees of the high-resolution N-body cosmological simulation $\mathrm{\texttt{Millennium-II}}$. The model shows that LISA MBHBs are expected to be found in optically dim ($r\,{>}\,20$), star-forming ($\rm sSFR\,{>}\,10^{-10}\, \rm yr^{-1}$), gas-rich ($f_{\rm gas}\,{>}\,0.6$) and disc-dominated ($\rm B/T\,{<}\,0.7$) \textit{low-mass galaxies} of stellar masses $10^8\,{-}\,10^9 M_{\odot}$. However, these properties are indistinguishable from those of galaxies harboring single massive black holes with comparable mass, making difficult the selection of LISA hosts among the whole population of low-mass galaxies. Motivated by this, we explore the possibility of using merger signatures to select LISA hosts. We find that 40-80% of the galaxies housing LISA MBHBs display merger features related to the interaction which brought the secondary MBH to the galaxy. Despite this, around 60% of dwarf galaxies placed in the surroundings of the LISA hosts will show such kind of features as well, challenging the unequivocal detection of LISA hosts through the search of merger signatures. Consequently, the detection of an electromagnetic transient associated with the MBHB merger will be vital to pinpoint the star-forming dwarf galaxy where these binary systems evolve and coalesce., Comment: 17 pages, 14 figures, Submitted to A&A

We present a study of seven spectroscopically confirmed (Ly$\alpha$ emitting) galaxies at redshift $z\simeq6$ using the $JWST$/NIRCam imaging data. These galaxies, with a wide range of Ly$\alpha$ luminosities, were recently observed in a series of NIRCam broad- and medium-bands. We constrain the rest-frame UV/optical continua and measure the H$\alpha$ line emission of the galaxies using the combination of the $JWST$/NIRCam and archival $HST$/WFC3 infrared photometry. We further estimate their escape fractions of Ly$\alpha$ photons ($f_{\rm esc}^{\rm Ly\alpha}$) and the production efficiency of ionizing photons ($\xi_{\rm ion}$). Among the sample, 6/7 galaxies have Ly$\alpha$ escape fractions of ${\lesssim}10\%$, which might be the status for most of star-forming galaxies at $z\simeq6$. One UV-faint Ly$\alpha$ galaxy with an extremely blue UV slope owns a large value of $f_{\rm esc}^{\rm Ly\alpha}$ reaching ${\simeq}50\%$. These galaxies spread a broad range of $\xi_{\rm ion}$ over log$_{10}$ $\xi_{\rm ion, 0}$ (Hz erg$^{-1}$) $\sim25.0-26.5$. We find that UV-fainter galaxies with bluer UV continuum slopes likely have higher escape fractions of Ly$\alpha$ photons. We also find that galaxies with higher Ly$\alpha$ line emission tend to produce ionizing photons more efficiently. The most Ly$\alpha$-luminous galaxy in the sample has a very high $\xi_{\rm ion, 0}$ of log$_{10}$ $\xi_{\rm ion, 0}$ (Hz erg$^{-1}$) $>26$. Our results support that Ly$\alpha$ galaxies may have served as an important contributor to the cosmic reionization. Blue and bright Ly$\alpha$ galaxies are excellent targets for $JWST$ follow-up spectroscopic observations., Comment: Accepted by ApJL; an improved version

The Ly$\alpha$ emission is an important tracer of neutral gas in a circum-galactic medium (CGM) around high-z QSOs. The origin of Lya emission around QSOs is still under debate which has significant implications for galaxy formation and evolution. In this paper, we study Ly$\alpha$ nebulae around two high redshift QSOs, SDSS J141935.58+525710.7 at $z=3.218$ (hereafter QSO1) and SDSS J141813.40+525240.4 at $z=3.287$ (hereafter QSO2), from the miniJPAS survey within the AEGIS field. Using the contiguous narrow-band (NB) images from the miniJPAS survey and SDSS spectra, we analyzed their morphology, nature, and origin. We report the serendipitous detection of double-core Ly\al\ morphology around two QSOs which is rarely seen among other QSOs. The separations of the two Ly\al~cores are 11.07 $\pm$ 2.26 kpcs (1.47 $\pm$ 0.3$^{\prime\prime}$) and 9.73 $\pm$ 1.55 kpcs (1.31 $\pm$ 0.21$^{\prime\prime}$) with Ly$\alpha$~line luminosities of $\sim$ 3.35 $\times 10^{44}$ erg s $^{-1} $ and $\sim$ 6.99 $\times$ 10$^{44}$ erg s $^{-1}$ for QSO1 and QSO2, respectively. The miniJPAS NB images show evidence of extended Ly$\alpha$ and CIV morphology for both QSOs and extended HeII morphology for QSO1. These two QSOs may be potential candidates for the new enormous Lyman alpha nebula (ELAN) found from the miniJPAS survey due to their extended morphology in the shallow depth and relatively high Ly$\alpha$ luminosities. We suggest that galactic outflows are the major powering mechanism for the double-core Ly$\alpha$ morphology. Considering the relatively shallow exposures of miniJPAS, the objects found here could be the tip of the iceberg of a promising number of such objects that will be uncovered in the upcoming full J-PAS survey and deep IFU observations with 8-10m telescopes will be essential for constraining the underlying physical mechanism that is responsible for the double-cored morphology., Comment: Five figures and two tables (Accepted in Astronomy and Astrophysics)

We study the genesis and evolution of super-massive black hole (SMBH) seeds through different formation channels, from PopIII remnants to massive seeds, modeled within the L-Galaxies semi-analytic code. We run the model on the Millennium-II simulation (MR-II) merger trees, as their halo-mass resolution (M_{vir,res}~10^7 Msun h^-1) allows to study in a cosmological volume (L_{box=100 Mpc h^-1) the evolution of atomic-cooling halos (T_{vir}>10^4 K) where intermediate-mass and heavy seeds are expected to form. We track the formation of these seeds according to spatial variations of the chemical and radiative feedback of star formation. Not being able to resolve the first mini-halos (T_{vir}~10^3 K), we inherit evolved PopIII remnants in a sub-grid fashion, using the results of the GQd model. We also include the formation of heavy seeds in gas-rich massive mergers, who are very rare in the MR-II volume. The descendants of light seeds numerically prevail among our SMBHs population at all masses and z. Heavier seeds form in dense environments where close neighbors provide the required UV illumination. Overall, our model produces a z=0 SMBHs population whose statistical properties meet current constraints. We find that the BH occupation fraction highly depends on the seeding efficiency and that the scaling relation between BH and stellar mass, in the dwarf-mass regime, is flatter than in the high-mass range. Finally, a fraction of BHs hosted in local dwarf galaxies never grow since they form at z>6., Comment: 23 pages, 13 figures, accepted in MNRAS

We make use of z = 0 samples of strongly barred and unbarred disc galaxies from the TNG100 and TNG50 cosmological hydrodynamical simulations to assess the performance of the simple disc instability criterion proposed by Efstathiou, Lake & Negroponte (1982) (ELN-criterion). We find that strongly barred galaxies generally assemble earlier, are more star-dominated in their central regions, and have more massive and more compact discs than unbarred galaxies. The ELN-criterion successfully identifies ~75% and ~80% of the strongly barred and the unbarred galaxies, respectively. Strongly barred galaxies that the criterion fails to identify tend to have more extended discs, higher spin values and bars that assembled later than is typical for the bulk of the barred population. The bars in many of these cases appear to be produced by an interaction with a close neighbour (i.e. to be externally triggered) rather than to result from secular growth in the disc. On the other hand, we find that unbarred galaxies misclassified as barred by the ELN-criterion typically have stellar discs similar to those of barred galaxies, although more extended in the vertical direction and less star-dominated in their central regions, possibly reflecting later formation times. In addition, the bulge component of these galaxies is significantly more prominent at early times than in the strongly barred sample. Thus, the ELN-criterion robustly identifies secular bar instabilities in most simulated disc galaxies, but additional environmental criteria are needed to account for interaction-induced bar formation., Comment: 15 pages, 12 Figures, Accepted by MNRAS

We present the photometric determination of the bright-end (L_Lya>10^43.5 erg/s) of the Lya luminosity function (LF) within four redshifts windows in the interval 2.210^43.3 erg/s) Lya-emitting candidates. We test our selection with two spectroscopic follow-up programs at the GTC telescope, confirming ~89% of the targets as line-emitting sources, with ~64% being genuine z~2.2 QSOs. We extend the 2.22 AGN/QSOs (likely dominant in our samples) are described by a structurally different LF than z>2 star-forming LAEs, namely with L^*_QSOs ~ 100 L^*_LAEs and Phi^*_QSOs ~ 10^-3 Phi^*_LAEs. Finally, our method identifies very efficiently as high-z line-emitters sources without previous spectroscopic confirmation, currently classified as stars (~2000 objects in each redshift bin, on average). Assuming a large predominance of Lya-emitting AGN/QSOs in our samples, this supports the scenario by which these are the most abundant class of z>2 Lya emitters at L_Lya>10^43.3 erg/s., Comment: 27 pages (6 in appendix), 21 figures (5 in appendix), submitted to A&A

We study the mass assembly and spin evolution of supermassive black holes (BHs) across cosmic time as well as the impact of gravitational recoil on the population of nuclear and wandering black holes (wBHs) by using the semi-analytical model LGalaxies run on top of Millennium merger trees. We track spin changes that BHs experience during both coalescence events and gas accretion phases. For the latter, we assume that spin changes are coupled with the bulge assembly. This assumption leads to predictions for the median spin values of z=0 BHs that depend on whether they are hosted by pseudobulges, classical bulges or ellipticals, being a ~ 0.9, 0.7 and 0.4, respectively. The outcomes of the model display a good consistency with z<4 quasar luminosity functions and the z=0 BH mass function, spin values and black hole-bulge mass correlation. Regarding the wBHs, we assume that they can originate from both the disruption of satellite galaxies (orphan wBH) and ejections due to gravitational recoils (ejected wBH). The model points to a number density of wBHs that increases with decreasing redshift, although this population is always ~2 dex smaller than the one of nuclear black holes. At all redshifts, wBHs are typically hosted in M_{halo}>10^{13} M_{\odot} and M_{stellar}>10^{10} M_{\odot}, being orphan wBHs the dominant type. Besides, independently of redshift and halo mass, \textit{ejected} wBHs inhabit the central regions (<0.3R_{200}) of the host DM halo, while orphan wBHs linger at larger scales (>0.5R_{200}). Finally, we find that gravitational recoils cause a progressive depletion of nuclear BHs with decreasing redshift and stellar mass. Moreover, ejection events lead to changes in the predicted local black hole-bulge relation, in particular for BHs in pseudobulges, for which the relation is flattened at M_{bulge}>10^{10.2} M{\odot} and the scatter increase up to ~3dex., Comment: 30 pages, 17 figures, accepted by MNRAS

We study the cosmological build-up of pseudobulges using the LGalaxies semi-analytical model for galaxy formation with a new approach for following separately the assembly of classical bulges and pseudobulges. Classical bulges are assumed to be the result of violent processes (i.e., mergers and starbursts), while the formation of pseudobulges is connected to the secular growth of disks. We apply the model to both the Millennium and the Millennium II simulations, in order to study our results across a wide range of stellar masses ($10^{7} - 10^{11.5} M_{\odot}$). We find that $z=0$ pseudobulges mainly reside in galaxies of $ \rm M_{stellar} \sim 10^{10} - 10^{10.5} M_{\odot}$ ($\rm M_{halo} \sim 10^{11.5}-10^{12} M_{\odot}$) and we recover structural properties of these objects (e.g., sizes and bulge-to-total ratios) that are in good agreement with observational results. Tracing their formation history, we find that pseudobulges assembled in galaxies with a very quiet merger history, as opposed to the host galaxies of classical bulges. Regarding the bulge structure, we find that $\sim$ 30\% of the galaxies with a predominant pseudobulge feature a composite structure, hosting both a pseudo and a classical bulge component. The classical component typically constitutes $\sim$10\% of the total bulge galaxy mass. When looking at the properties of the host galaxies, we find that $z = 0$ pseudobulges are hosted by main sequence galaxies, characterized by a stellar population which is generally younger compared to the one of the hosts of classical bulges., Comment: Submitted to MNRAS

The emergence of a large-scale stellar bar is one of the most striking features in disc galaxies. By means of state-of-the-art cosmological zoom-in simulations, we study the formation and evolution of bars in Milky Way-like galaxies in a fully cosmological context, including the physics of gas dissipation, star formation, and supernova feedback. Our goal is to characterise the actual trigger of the non-axisymmetric perturbation that leads to the strong bar observable in the simulations at z=0, discriminating between an internal/secular versus an external/tidal origin. To this aim, we run a suite of cosmological zoom-in simulations altering the original history of galaxy-satellite interactions at a time when the main galaxy, though already bar-unstable, does not feature any non-axisymmetric structure yet. We find that the main effect of a late minor merger and of a close fly-by is to delay the time of bar formation and those two dynamical events are not directly responsible for the development of the bar and do not alter significantly its global properties (e.g. its final extension). We conclude that, once the disc has grown to a mass large enough to sustain global non-axisymmetric modes, then bar formation is inevitable., Comment: 15 pages, 6 Figure; Accepted for publication in MNRAS - Full-resolution images are available at http://www.dfm.uninsubria.it/fh/FHpages/Research.html

We analyse the output of the hi-res cosmological zoom-in simulation ErisBH to study self-consistently the formation of a strong stellar bar in a Milky Way-type galaxy and its effect on the galactic structure, on the central gas distribution and on star formation. The simulation includes radiative cooling, star formation, SN feedback and a central massive black hole which is undergoing gas accretion and is heating the surroundings via thermal AGN feedback. A large central region in the ErisBH disk becomes bar-unstable after z~1.4, but a clear bar-like structure starts to grow significantly only after z~0.4, possibly triggered by the interaction with a massive satellite. At z~0.1 the bar reaches its maximum radial extent of l~2.2 kpc. As the bar grows, it becomes prone to buckling instability, which we quantify based on the anisotropy of the stellar velocity dispersion. The actual buckling event is observable at z~0.1, resulting in the formation of a boxy-peanut bulge clearly discernible in the edge-on view of the galaxy at z=0. The bar in ErisBH does not dissolve during the formation of the bulge but remains strongly non-axisymmetric down to the resolution limit of ~100 pc at z=0. During its early growth, the bar exerts a strong torque on the gas within its extent and drives gas inflows that enhance the nuclear star formation on sub-kpc scales. Later on the infalling gas is nearly all consumed into stars and, to a lesser extent, accreted onto the central black hole, leaving behind a gas-depleted region within the central ~2 kpc. Observations would more likely identify a prominent, large-scale bar at the stage when the galactic central region has already been quenched. Bar-driven quenching may play an important role in disk-dominated galaxies at all redshift. [Abridged], Comment: 13 pages, 12 figures, MNRAS submitted