Search

Submillimeter galaxies (SMGs) constitute a key population of bright star-forming galaxies at high redshift. These galaxies challenge galaxy formation models, particularly in reproducing their observed number counts and redshift distributions. Furthermore, although SMGs contribute significantly to the cosmic star formation rate density (SFRD), their precise role remains uncertain. Upcoming surveys, such as the Ultra Deep Survey with the TolTEC camera, are expected to offer valuable insights into SMG properties and their broader impact. Robust modeling of SMGs in a cosmologically representative volume is necessary to investigate their nature in preparation for next-generation submillimeter surveys. We implement and test parametric relations derived from radiative transfer calculations across three cosmological simulations: EAGLE, IllustrisTNG, and FLAMINGO. Particular emphasis is placed on the FLAMINGO due to their large volume and robust statistical sampling of SMGs. Based on the model that best reproduces observations, we forecast submillimeter fluxes within the simulations, analyze the properties of SMGs, and evaluate their evolution over cosmic time. Our results show that the FLAMINGO reproduces the observed redshift distribution and source number counts of SMGs without requiring a top-heavy initial mass function. On the other hand, the EAGLE and IllustrisTNG show a deficit of bright SMGs. We find that SMGs with S850 > 1 mJy contribute up to 27% of the SFRD at z=2.6 in the FLAMINGO, consistent with recent observations. Flux density functions reveal a rise in SMG abundance from z = 6 to 2.5, followed by a sharp decline in the number of brighter SMGs from z = 2.5 to 0. Leveraging the SMG population in FLAMINGO, we forecast that the TolTEC UDS will detect 80,000 sources over 0.8 deg^2 at 1.1 mm (at the 4{\sigma} detection limit), capturing about 50% of the cosmic SFRD at z=2.5., Comment: 19 pages, 13 figures, submitted to Astronomy & Astrophysics Journal, abstract is reduced to meet Arxiv character limit

Our concern is the data complexity of answering linear monadic datalog queries whose atoms in the rule bodies can be prefixed by operators of linear temporal logic LTL. We first observe that, for data complexity, answering any connected query with operators $\bigcirc/\bigcirc^-$ (at the next/previous moment) is either in AC0, or in $ACC0\!\setminus\!AC0$, or $NC^1$-complete, or LogSpace-hard and in NLogSpace. Then we show that the problem of deciding LogSpace-hardness of answering such queries is PSpace-complete, while checking membership in the classes AC0 and ACC0 as well as $NC^1$-completeness can be done in ExpSpace. Finally, we prove that membership in AC0 or in ACC0, $NC^1$-completeness, and LogSpace-hardness are undecidable for queries with operators $\Diamond_f/\Diamond_p$ (sometime in the future/past) provided that $NC^1 \ne NLogSpace$, and $LogSpace \ne NLogSpace$., Comment: Extended version of a paper accepted at ICDT'2025

Investigating the assembly history of the most massive and passive galaxies will enhance our understanding of why galaxies exhibit such a remarkable diversity in structure and morphology. In this paper, we simultaneously investigate the assembly history and redshift evolution of semi-analytically modelled galaxy properties of central galaxies between 0.56 < z < 4.15, alongside their connection to their halos as a function of large-scale environment. We extract sub-samples of galaxies from a mock catalogue representative for the BOSS-CMASS sample, which includes the most massive and passively evolving system known today. Utilising typical galaxy properties such as star formation rate, (g-i) colour, or cold gas-phase metallicity (Zcold), we track the redshift evolution of these properties across the main progenitor trees. We present results on galaxy and halo properties, including their growth and clustering functions. Our findings indicate that galaxies in the highest stellar and halo mass regimes are least metal-enriched (using Zcold as a proxy) and consistently exhibit significantly larger black hole masses and higher clustering amplitudes compared to sub-samples selected by e.g. colour or star formation rate. This population forms later and also retains large reservoirs of cold gas. In contrast, galaxies in the intermediate and lower stellar/halo mass regimes consume their cold gas at higher redshift and were among the earliest and quickest to assemble. We observe a clear trend where the clustering of the galaxies selected according to their Zcold-values (either low-Zcold or high-Zcold) depends on the density of their location within the large-scale environment. We assume that in particular galaxies in the low/high-Zcold sub-samples form and evolve through distinct evolutionary channels, which are predetermined by their location within the large-scale environment of the cosmic web., Comment: 21 pages, 12 figures, accepted for publication in A&A

The aim of this study is to identify quiescent galaxies in the 2-deg$^2$ COSMOS field at $z \sim 3.1$ and analyze their environment. Using data from the ODIN survey and COSMOS2020 catalog, we identify 24 massive quiescent galaxies (MQGs) with stellar masses $\geq 10^{10.6}$ and derive their star formation histories and quenching timescales using SED fitting with BAGPIPES. Voronoi-based density maps trace local and large-scale environments using Lyman-$\alpha$ Emitters and photometric galaxies. Results indicate uniformly short quenching timescales ($<$500 Myr) independent of environmental density, suggesting rapid internal mechanisms such as AGN feedback dominate over environmental factors. MQGs do not correlate with protoclusters or filaments, although some are near gas-rich filaments but show no rejuvenation. These findings suggest quenching at high redshift is driven primarily by internal processes rather than environmental interactions.

Context: Understanding galaxy evolution in dense environments, particularly proto-clusters, is crucial for studying mechanisms driving star formation and quenching. Aims: This study examines how two proto-cluster over-densities at 3 < z < 4 impact star formation rate (SFR), stellar mass, and morphology, focusing on quenched galaxies. Methods: We identified proto-cluster over-densities in the Chandra Deep Field South (CDFS) and Ultra Deep Survey (UDS) regions of the VANDELS survey. Using spectral energy distribution analysis, Bayesian methods (BEAGLE and BAGPIPES) helped derive best-fit parameters and U-V and V-J rest-frame colours (UVJ), classifying galaxies as quenched or star-forming based on UVJ diagrams and specific star formation rates (sSFR). TNG300 simulations aided interpretation. Results: Two of 13 proto-cluster over-densities host quenched galaxies with red U-V colours, low sSFR, and properties like massive passive galaxies. These quenched members are redder, older, more massive, and more compact. The highest-density peaks at z=3.55 and z=3.43 have dark matter halo masses consistent with proto-clusters and host AGNs, with five and three AGNs, respectively. Compared to field galaxies, these quenched members are in denser environments. TNG300 simulations suggest proto-clusters with quenched galaxies at high redshift evolve to contain more passive galaxies by z=1. Conclusions: The over-densities host massive quenched galaxies and AGNs in their densest peaks. Simulations reveal that sSFR for passive galaxies in proto-clusters was high at z=6, with median mass growth rates of 96% from z=6 to z=3. Conditions for mass assembly likely involve galaxy interactions and high gas accretion in dense environments. Black hole growth and AGN feedback appear to drive quenching at z=3, aligning with the properties of quenched galaxies observed in our study., Comment: 16 pages, 13 figures, 1 table. Published in A&A

As the number of gravitational-wave detections grows, the merger rate of binary black holes (BBHs) can help us to constrain their formation, the properties of their progenitors, and their birth environment. Here, we aim to address the impact of the metal-dependent star formation rate (SFR) on the BBH merger rate. To this end, we have developed a fully data-driven approach to model the metal-dependent SFR and coupled it to BBH evolution. We have adopted the most up-to-date scaling relations, based on recent observational results, and we have studied how the BBH merger rate density varies over a wide grid of galaxy and binary evolution parameters. Our results show that including a realistic metal-dependent SFR evolution yields a value of the merger rate density which is too high compared to the one inferred from GW data. Moreover, variations of the SFR in low-mass galaxies ($M_\ast \lesssim 10^8 \mathrm{M}_{\odot}$) do not contribute more than a factor $\sim 2$ to the overall merger rate density at redshift $z=0$. These results suggest that the discrepancy between the BBH merger rate density inferred from data and theoretical models is not caused by approximations in the treatment of the metal-dependent SFR, but rather stems from stellar evolution models and/or BBH formation channels., Comment: 14 pages, 9 figures. Submitted to A&A, comments are welcome

The One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) survey is carrying out a systematic search for protoclusters during Cosmic Noon, using Ly$\alpha$-emitting galaxies (LAEs) as tracers. Once completed, ODIN aims to identify hundreds of protoclusters at redshifts of 2.4, 3.1, and 4.5 across seven extragalactic fields, covering a total area of up to 91~deg$^2$. In this work, we report strong clustering of high-redshift protoclusters through the protocluster-LAE cross-correlation function measurements of 150 protocluster candidates at $z~=~2.4$ and 3.1, identified in two ODIN fields with a total area of 13.9 deg$^2$. At $z~=~2.4$ and 3.1, respectively, the inferred protocluster biases are $6.6^{+1.3}_{-1.1}$ and $6.1^{+1.3}_{-1.1}$, corresponding to mean halo masses of $\log \langle M /M_\odot\rangle = 13.53^{+0.21}_{-0.24}$ and $12.96^{+0.28}_{-0.33}$. By the present day, these protoclusters are expected to evolve into virialized galaxy clusters with a mean mass of $\sim$ $10^{14.5}~M_\odot$. By comparing the observed number density of protoclusters to that of halos with the measured clustering strength, we find that our sample is highly complete. Finally, the similar descendant masses derived for our samples at $z=2.4$ and 3.1 assuming that the halo number density remains constant suggest that they represent similar structures observed at different cosmic epochs. As a consequence, any observed differences between the two samples can be understood as redshift evolution. The ODIN protocluster samples will thus provide valuable insights into the cosmic evolution of cluster galaxies., Comment: 12 pages, 5 figures, submitted to ApJ

The large-scale linear halo bias encodes the relation between the clustering of dark-matter (DM) halos and that of the underlying matter density field. Although the primary dependence of bias on halo mass is well understood in the context of structure formation, the physical origins of the multiple additional relations at fixed halo mass, commonly known as secondary halo bias, have not been fully elucidated. Of particular relevance is the secondary dependence on halo assembly history, known as halo assembly bias. Our goal is to determine whether the properties of the initial regions from which $z=0$ halos originate produce any secondary bias at $z=0$. By analyzing these initial dependencies in connection with halo assembly bias, we intend to provide insight on the physical origins of the effect. To this end, we select halos at $z=0$ in the IllustrisTNG DM-only simulation and trace back the positions and velocities of their DM particles to $z=12$. The resulting initial regions are characterized according to several shape-related and kinematic properties. The secondary bias signal produced by these properties at $z=0$ is measured using an object-by-object bias estimator, which offers significant analytical advantages. We show that, when split by the properties of their initial DM clouds, $z=0$ halos display significant secondary bias, clearly exceeding the amplitude of the well-known halo assembly bias signal produced by concentration and age. The maximum bias segregation is measured for cloud velocity dispersion and radial velocity, followed by cloud concentration, sphericity, ellipticity and triaxiality. We further show that both velocity dispersion and radial velocity are also the properties of the initial clouds that most strongly correlate with halo age and concentration at fixed halo mass. Our results highlight the importance of linear effects in shaping halo assembly bias., Comment: 9 pages, 5 figures. Submitted to A&A, comments welcome

We investigate the physical properties and redshift evolution of simulated galaxies residing in protoclusters at cosmic noon, to understand the influence of the environment on galaxy formation. This work is to build clear expectations for the ongoing ODIN survey, devoted to mapping large-scale structures at z=2.4, 3.1, and 4.5 using Ly$\alpha$-emitting galaxies (LAEs) as tracers. From the IllustrisTNG simulations, we define subregions centered on the most massive clusters ranked by total stellar mass at z=0 and study the properties of galaxies within, including LAEs. To model the LAE population, we take a semi-analytical approach that assigns Ly$\alpha$ luminosity and equivalent width based on the UV luminosities to galaxies in a probabilistic manner. We investigate stellar mass, star formation rate, major mergers, and specific star formation rate of the population of star-forming galaxies and LAEs in the field and protocluster environment and trace their evolution. We find that the overall shape of the UV luminosity function (LF) in simulated protocluster environments is characterized by a shallower faint-end slope and an excess on the bright end, signaling different formation histories for galaxies therein. The difference is milder for the Ly$\alpha$ LF. While protocluster galaxies follow the same SFR-$M_{\odot}$ scaling relation as average field galaxies, a larger fraction appears to have experienced major mergers in the last 200 Myr and as a result shows enhanced star formation at a ~60% level, leading to a flatter distribution in both SFR and $M_{\odot}$ relative to galaxies in the average field. We find that protocluster galaxies, including LAEs, begin to quench much earlier (z~0.8-1.6) than field galaxies (z~0.5-0.9); our result is in agreement with recent observational results and highlights the importance of large-scale environment on the overall formation history of galaxies., Comment: 18 pages, 13 figures, to be submitted to Astronomy & Astrophysics, the abstract was reduced to meet character requirements; edits include: updating coauthor list and affiliations, updating citations in Sec 2.1

We test whether Lyman alpha emitters (LAEs) and Lyman-break galaxies (LBGs) can be good tracers of high-z large-scale structures, using the Horizon Run 5 cosmological hydrodynamical simulation. We identify LAEs using the Ly{\alpha} emission line luminosity and its equivalent width, and LBGs using the broad-band magnitudes at z~2.4, 3.1, and 4.5. We first compare the spatial distributions of LAEs, LBGs, all galaxies, and dark matter around the filamentary structures defined by dark matter. The comparison shows that both LAEs and LBGs are more concentrated toward the dark matter filaments than dark matter. We also find an empirical fitting formula for the vertical density profile of filaments as a binomial power-law relation of the distance to the filaments. We then compare the spatial distributions of the samples around the filaments defined by themselves. LAEs and LBGs are again more concentrated toward their filaments than dark matter. We also find the overall consistency between filamentary structures defined by LAEs, LBGs, and dark matter, with the median spatial offsets that are smaller than the mean separation of the sample. These results support the idea that the LAEs and LBGs could be good tracers of large-scale structures of dark matter at high redshifts., Comment: 20 pages, 15 figures, 2 tables, accepted for publication in ApJ (affiliations are updated in v2)

We investigated the population of low-mass central galaxies with Mstar = $10^{9.5}-10^{10}$ Msun/h, inhabiting regions near massive groups and clusters of galaxies using the TNG300 and MDPL2-SAG simulations. We set out to study their evolutionary histories, aiming to find hints about the large-scale conformity signal they produce. We also used a control sample of central galaxies with the same stellar mass range located far away from massive structures. For both samples, we find a subpopulation of galaxies accreted by another halo in the past, but now considered central galaxies; we refer to these objects as former satellites. The number of former satellites is higher for quenched central galaxies near massive systems, with fractions of 45% and 17% in TNG300 and MDPL2-SAG. Our results in TNG300 show that former satellites pollute the sample of central galaxies because they suffered environmental processes when they were satellites hosted typically by massive dark matter halos (M200 $\geq 10^{13}$ Msun/h) since z$\lesssim$0.5. After removing former satellites, the evolutionary trends for quenched central galaxies near massive structures are fairly similar to those of the quenched control galaxies, showing small differences at low redshift. For MDPL2-SAG instead, former satellites were hosted by less massive halos, with a mean halo mass around $10^{11.4}$ Msun/h, and the evolutionary trends remain equal before and after removing former satellite galaxies. We also measured the two-halo conformity, i.e, the correlation in the sSFR between low-mass central galaxies and their neighbors at Mpc scales, and how former satellites contribute to the signal at three different redshifts: z=0, 0.3, and 1. The time evolution of the conformity signal in the simulations presents apparent contradictory results: it decreases from z=0 to z=1 in MDPL2-SAG, while it increases in TNG300 (abridged)., Comment: 17 pages, 11 figures. Published in A&A

To understand the formation and evolution of massive cosmic structures, studying them at high redshift, in the epoch when they formed the majority of their mass is essential. The One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) survey is undertaking the widest-area narrowband program to date, to use Ly$\alpha$-emitting galaxies (LAEs) to trace the large-scale structure (LSS) of the Universe on the scale of 10 - 100 cMpc at three cosmic epochs. In this work, we present results at $z$ = 3.1 based on early ODIN data in the COSMOS field. We identify and characterize protoclusters and cosmic filaments using multiple methods and discuss their strengths and weaknesses. We then compare our observations against the IllustrisTNG suite of cosmological hydrodynamical simulations. The two are in excellent agreement, with a similar number and angular size of structures identified above a specified density threshold. We are able to recover the simulated protoclusters with $\log$(M$_{z=0}$/$M_\odot$) $\gtrsim$ 14.4 in $\sim$ 60% of the cases. With these objects we show that the descendant masses of the protoclusters in our sample can be estimated purely based on our 2D measurements, finding a median $z$ = 0 mass of $\sim10^{14.5}$M$_\odot$. The lack of information on the radial extent of each protocluster introduces a $\sim$0.4 dex uncertainty in its descendant mass. Finally, we show that the recovery of the cosmic web in the vicinity of protoclusters is both efficient and accurate. The similarity of our observations and the simulations imply that our structure selection is likewise robust and efficient, demonstrating that LAEs are reliable tracers of the LSS., Comment: 29 pages, 20 figures; Accepted for publication in ApJ

Definite descriptions, such as 'the General Chair of KR 2024', are a semantically transparent device for object identification in knowledge representation. In first-order modal logic, definite descriptions have been widely investigated for their non-rigidity, which allows them to designate different objects (or none at all) at different states. We propose expressive modal description logics with non-rigid definite descriptions and names, and investigate decidability and complexity of the satisfaction problem. We first systematically link satisfiability for the one-variable fragment of first-order modal logic with counting to our modal description logics. Then, we prove a promising NEXPTIME-completeness result for concept satisfiability for the fundamental epistemic multi-agent logic $\mathbf{S5}^{n}$ and its neighbours, and show that some expressive logics that are undecidable with constant domain become decidable (but Ackermann-hard) with expanding domains. Finally, we conduct a fine-grained analysis of decidability of temporal logics.

Observations suggest that red central galaxies align closely with their group galaxies and the large-scale environment.This finding was also replicated in simulations, which added information about the alignment of the stars that form the galaxies with the dark matter in the halo they inhabit. These results were obtained for the present universe. Our study aims to build upon previous findings by examining the evolution of central galaxy alignment with their environment, as well as the alignment between their stellar and dark matter components. Based on previous studies, in this work, we describe the evolution of the alignment of bright central galaxies over time and try to understand the process leading to the current observed alignment. By employing the merger trees from the simulation, we track the alignment evolution of the central galaxy sample at z=0 used in a previous study, whose results correspond to the observations. In particular, we exploit the anisotropic correlation function to study the alignment of the central galaxies with the environment and the probability distribution of the angle between the axes of the shape tensor calculated for each component to deepen the analysis of the stellar and dark matter components. A description was given of the evolution of alignment in bright central galaxies with a focus on the distinctions between red and blue galaxies. Furthermore, it was found that the alignment of the dark matter halo differs from that of the stellar material within it. According to the findings, the assembly process and mergers influenced the evolution of alignment, Comment: 11 pages, 13 figuras, Accepted in A&A

Hierarchical black hole (BH) mergers are one of the most straightforward mechanisms to produce BHs inside and above the pair-instability mass gap. Here, we investigate the impact of globular cluster (GC) evolution on hierarchical mergers, and we account for the uncertainties related to BH mass pairing functions on the predicted primary BH mass, mass ratio and spin distribution. We find that the evolution of the host GC quenches the hierarchical BH assembly already at the third generation, mainly due to cluster expansion powered by a central BH sub-system. Hierarchical mergers match the primary BH mass distribution from GW events for $m_1 > 50 \, \mathrm{M_{\odot}}$, regardless of the assumed BH pairing function. At lower masses, however, different pairing functions lead to dramatically different predictions on the primary BH mass merger rate density. We find that the primary BH mass distribution evolves with redshift, with a larger contribution from mergers with $m_1 \geq 30 \, \mathrm{M_{\odot}}$ for $z\geq{}2$. Finally, we calculate the mixing fraction of BBHs from GCs and isolated binary systems. Our predictions are very sensitive to the spins, which favor a large fraction ($>0.6$) of BBHs born in GCs, in order to reproduce misaligned spin observations., Comment: 18 pages, 15 figures, 1 table. Comments welcome

Lyman-Alpha Emitting galaxies (LAEs) are typically young, low-mass, star-forming galaxies with little extinction from interstellar dust. Their low dust attenuation allows their Ly$\alpha$ emission to shine brightly in spectroscopic and photometric observations, providing an observational window into the high-redshift universe. Narrowband surveys reveal large, uniform samples of LAEs at specific redshifts that probe large scale structure and the temporal evolution of galaxy properties. The One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) utilizes three custom-made narrowband filters on the Dark Energy Camera (DECam) to discover LAEs at three equally spaced periods in cosmological history. In this paper, we introduce the hybrid-weighted double-broadband continuum estimation technique, which yields improved estimation of Ly$\alpha$ equivalent widths. Using this method, we discover 6032, 5691, and 4066 LAE candidates at $z =$ 2.4, 3.1, and 4.5 in the extended COSMOS field ($\sim$9 deg$^2$). We find that [O II] emitters are a minimal contaminant in our LAE samples, but that interloping Green Pea-like [O III] emitters are important for our redshift 4.5 sample. We introduce an innovative method for identifying [O II] and [O III] emitters via a combination of narrowband excess and galaxy colors, enabling their study as separate classes of objects. We present scaled median stacked SEDs for each galaxy sample, revealing the overall success of our selection methods. We also calculate rest-frame Ly$\alpha$ equivalent widths for our LAE samples and find that the EW distributions are best fit by exponential functions with scale lengths of $w_0$ = 53 $\pm$ 1, 65 $\pm$ 1, and 59 $\pm$ 1 Angstroms, respectively., Comment: Accepted to ApJ, 24 pages, 14 figures

We describe the survey design and science goals for ODIN (One-hundred-deg^2 DECam Imaging in Narrowbands), a NOIRLab survey using the Dark Energy Camera (DECam) to obtain deep (AB~25.7) narrow-band images over an unprecedented area of sky. The three custom-built narrow-band filters, N419, N501, and N673, have central wavelengths of 419, 501, and 673 nm and respective full-widthat-half-maxima of 7.2, 7.4, and 9.8 nm, corresponding to Lya at z=2.4, 3.1, and 4.5 and cosmic times of 2.8, 2.1, and 1.4 Gyr, respectively. When combined with even deeper, public broad-band data from Hyper Suprime-Cam, DECam, and in the future, LSST, the ODIN narrow-band images will enable the selection of over 100,000 Lya-emitting (LAE) galaxies at these epochs. ODIN-selected LAEs will identify protoclusters as galaxy overdensities, and the deep narrow-band images enable detection of highly extended Lya blobs (LABs). Primary science goals include measuring the clustering strength and dark matter halo connection of LAEs, LABs, and protoclusters, and their respective relationship to filaments in the cosmic web. The three epochs allow the redshift evolution of these properties to be determined during the period known as Cosmic Noon, where star formation was at its peak. The two narrow-band filter wavelengths are designed to enable interloper rejection and further scientific studies by revealing [O II] and [O III] at z=0.34, Lya and He II 1640 at z=3.1, and Lyman continuum plus Lya at z=4.5. Ancillary science includes similar studies of the lower-redshift emission-line galaxy samples and investigations of nearby star-forming galaxies resolved into numerous [O III] and [S II] emitting regions., Comment: 21 pages, 8 figures, submitted to Astrophysical Journal

The large-scale distribution of baryons is sensitive to gravitational collapse, mergers, and galactic feedback. Known as the Cosmic Web, its large-scale structure (LSS) can be classified as halos, filaments, and voids. Fast Radio Bursts (FRBs) are extragalactic sources that undergo dispersion along their propagation paths. They provide insight into ionised matter along their sightlines via their dispersion measures (DMs), and have been investigated as probes of the LSS baryon fraction, the diffuse baryon distribution, and of cosmological parameters. We use the cosmological simulation IllustrisTNG to study FRB DMs accumulated while traversing different types of LSS. We combine methods for deriving electron density, classifying LSS, and tracing FRB sightlines. We identify halos, filaments, voids, and collapsed structures along random sightlines and calculate their DM contributions. We analyse the redshift-evolving cosmological DM components of the Cosmic Web. We find that the filamentary contribution dominates, increasing from ~71% to ~80% on average for FRBs originating at z=0.1 vs z=5, while the halo contribution falls, and the void contribution remains consistent to within ~1%. The majority of DM variance originates from halos and filaments, potentially making void-only sightlines more precise probes of cosmological parameters. We find that, on average, an FRB originating at z=1 will intersect ~1.8 foreground collapsed structures, increasing to ~12.4 structures for a z=5 FRB. The impact parameters between our sightlines and TNG structures of any mass appear consistent with those reported for likely galaxy-intersecting FRBs. However, we measure lower average accumulated DMs from these structures than the $\sim90\;{\rm pc\;cm^{-3}}$ DM excesses reported for these literature FRBs, indicating some DM may arise beyond the structures themselves., Comment: 19 pages, 10 figures, accepted for publication by Astronomy & Astrophysics on 14/09/2023

Definite descriptions, such as 'the smallest planet in the Solar System', have been recently recognised as semantically transparent devices for object identification in knowledge representation formalisms. Along with individual names, they have been introduced also in the context of description logic languages, enriching the expressivity of standard nominal constructors. Moreover, in the first-order modal logic literature, definite descriptions have been widely investigated for their non-rigid behaviour, which allows them to denote different objects at different states. In this direction, we introduce epistemic and temporal extensions of standard description logics, with nominals and the universal role, additionally equipped with definite descriptions constructors. Regarding names and descriptions, in these languages we allow for: possible lack of denotation, ensured by partial models, coming from free logic semantics as a generalisation of the classical ones; and non-rigid designation features, obtained by assigning to terms distinct values across states, as opposed to the standard rigidity condition on individual expressions. In the absence of the rigid designator assumption, we show that the satisfiability problem for epistemic free description logics is NExpTime-complete, while satisfiability for temporal free description logics over linear time structures is undecidable.

We perform a detailed study of the cosmological bias of gravitational gave (GW) events produced by binary black hole mergers (BBHM). We start from a BBHM distribution modeled inside the EAGLE hydrodyamical simulation using the population synthesis code MOBSE. We then compare our findings with predictions from different Halo Occupation Distribution (HOD) prescriptions and find overall agreement, provided that the modeled properties of host galaxies and halos in the semi-analytical treatment match those in the simulations. By highlighting the sources of these discrepancies, we provide the stepping stone to build future more robust models that prevent the shortcoming of both simulation-based and analytical models. Finally, we train a neural network to build a simulation-based HOD and perform feature importance analysis to gain intuition on which host halo/galaxy parameters are the most relevant in determining the actual distribution and power spectrum of BBHM. We find that the distribution of BBHM in a galaxy does not only depend on its size, star formation rate and metallicity, but also by its kinetic state.

Galactic binary neutron stars (BNSs) are a unique laboratory to probe the evolution of BNSs and their progenitors. Here, we use a new version of the population synthesis code SEVN to evolve the population of Galactic BNSs, by modeling the spin up and down of pulsars self-consistently. We analyze the merger rate $\mathcal{R}_{\rm MW}$, orbital period $P_{\rm orb}$, eccentricity $e$, spin period $P$, and spin period derivative $\dot{P}$ of the BNS population. Values of the common envelope parameter $\alpha=1 - 3$ and an accurate model of the Milky Way star formation history best reproduce the BNS merger rate in our Galaxy ($\mathcal{R}_{\rm MW}\approx{}30$ Myr$^{-1}$). We apply radio-selection effects to our simulated BNSs and compare them to the observed population. Using a Dirichlet process Gaussian mixture method, we evaluate the four-dimensional likelihood in the $(P_{\rm orb}, e, P, \dot{P})$ space, by comparing our radio-selected simulated pulsars against Galactic BNSs. Our analysis favours an uniform initial distribution for both the magnetic field ($10^{10-13}$ G) and the spin period ($10-100$ ms). The implementation of radio selection effects is critical to match not only the spin period and period derivative, but also the orbital period and eccentricity of Galactic BNSs. According to our fiducial model, the Square Kilometre Array will detect $\sim 20$ new BNSs in the Milky Way., Comment: 20 pages, 8 figures, 8 tables, comments welcome

In this work, we aim at investigating the morphology evolution of Milky Way mass-like dark matter haloes selected from the CIELO and IllustrisTNG Projects. The connection between halo shapes and their environment has been studied in previous works at z=0 but their connection remains yet to be fully understood. We focus on the evolution across cosmic time of the halo shapes and the relation with the infalling material, using hydrodynamical simulations. Our findings show that haloes tend to be more triaxial at earlier times as a consequence of stronger accretion in the direction of the filaments. As the haloes evolve towards a dominant isotropic accretion mode and relaxation, their shape at 20 percent of the virial mass becomes more spherical. In agreement with previous results, baryons have an important effect within the inner regions of the haloes, driving them from triaxial to rounder shapes. We also find a correlation between the strength of the quadrupole infalling mode and the degree of ellipticity of the haloes: as the filament strength decreases steadily with redshift, the haloes became more spherical and less elliptical.

While many Lyman-alpha Blobs (LABs) are found in and around several well-known protoclusters at high redshift, how they trace the underlying large-scale structure is still poorly understood. In this work, we utilize 5,352 Lyman-alpha emitters (LAEs) and 129 LABs at z=3.1 identified over a $\sim$ 9.5 sq. degree area in early data from the ongoing One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) survey to investigate this question. Using LAEs as tracers of the underlying matter distribution, we identify overdense structures as galaxy groups, protoclusters, and filaments of the cosmic web. We find that LABs preferentially reside in regions of higher-than-average density and are located in closer proximity to overdense structures, which represent the sites of protoclusters and their substructures. Moreover, protoclusters hosting one or more LABs tend to have a higher descendant mass than those which do not. Blobs are also strongly associated with filaments of the cosmic web, with $\sim$ 70% of the population being within a projected distance of 2.4 pMpc from a filament. We show that the proximity of LABs to protoclusters is naturally explained by their association with filaments as large cosmic structures are where many filaments converge. The contiguous wide-field coverage of the ODIN survey allows us for the first time to firmly establish a connection between LABs as a population and their environment., Comment: 25 pages, 17 figures; accepted for publication in ApJ

Satellite galaxies undergo a variety of physical processes when they are accreted by groups and clusters, often resulting in the loss of baryonic and dark matter (DM) mass. In this work, we evaluate the predictions from the IllustrisTNG hydrodynamical simulation regarding the evolution of the matter content of satellites, focusing on a population that are accreted at $z>1$ and retain their identity as satellites down to $z=0$. At fixed host halo mass, the amount of DM and stellar mass stripped depends mostly on the pericentric distance, $d_{\rm peri}$, here normalised by host halo virial radius. The closest encounters result in significant loss of DM, with subhaloes retaining between 20 and a few per cent of their $z=1$ mass. At fixed $d_{\rm peri}$, DM mass stripping seems more severe in lower mass haloes. Conversely, the average satellite in higher mass haloes has its stellar mass growth halted earlier, having lost a higher fraction of stellar mass by $z=0$. We also show that mass stripping has a strong impact on the quenched fractions. The IllustrisTNG boxes are qualitatively consistent in these predictions, with quantitative differences mostly originating from the distinct subhalo mass ranges covered by the boxes. Finally, we have identified DM-deficient systems in all TNG boxes. These objects are preferentially found in massive clusters ($M_{\rm host } \gtrsim 10^{13}$ M$_\odot$), had very close encounters with their central galaxies ($d_{\rm peri}\simeq0.05\, R_{\rm vir}$), and were accreted at high redshift ($z_{\rm infall} \gtrsim 1.4$), reinforcing the notion that tidal stripping is responsible for their remarkable lack of DM., Comment: 16 pages, 12 figures. Submitted to MNRAS, comments welcome

Motivated by observational results, we use IllustrisTNG hydrodynamical numerical simulations to study the alignment of the central galaxies in groups with the surrounding structures. This approach allows us to analyse galaxy and group properties not available in observations. To perform this analysis, we use a modified version of the two-point cross-correlation function and a measure of the angle between the semi-major axes of the central galaxies and the larger structures. Overall, our results reproduce observational ones, as we find large-scale anisotropy, which is dominated by the red central galaxies. In addition, the latter is noticeably more aligned with their group than the blue ones. In contrast to the observations, we find a strong dependence of the anisotropy on the central galaxy with mass, probably associated with the inability of observational methods to determine them. This result allows us to link the alignment to the process of halo assembly and the well-known dependence of halo anisotropy on mass. When we include the dark matter distribution in our analysis, we conclude that the galaxy alignment found in simulations (and observations) can be explained by a combination of physical processes at different scales: the central galaxy aligns with the dark matter halo it inhabits, and this, in turn, aligns with the surrounding structures at large scales., Comment: 9 pages, 11 figures, Accepted by MNRAS

Linear Temporal Logic (LTL) is the de-facto standard temporal logic for system specification, whose foundational properties have been studied for over five decades. Safety and cosafety properties define notable fragments of LTL, where a prefix of a trace suffices to establish whether a formula is true or not over that trace. In this paper, we study the complexity of the problems of satisfiability, validity, and realizability over infinite and finite traces for the safety and cosafety fragments of LTL. As for satisfiability and validity over infinite traces, we prove that the majority of the fragments have the same complexity as full LTL, that is, they are PSPACE-complete. The picture is radically different for realizability: we find fragments with the same expressive power whose complexity varies from 2EXPTIME-complete (as full LTL) to EXPTIME-complete. Notably, for all cosafety fragments, the complexity of the three problems does not change passing from infinite to finite traces, while for all safety fragments the complexity of satisfiability (resp., realizability) over finite traces drops to NP-complete (resp., ${\Pi}^P_2$-complete).

The merger rate density evolution of binary compact objects and the properties of their host galaxies carry crucial information to understand the sources of gravitational waves. Here, we present galaxyRate, a new code that estimates the merger rate density of binary compact objects and the properties of their host galaxies, based on observational scaling relations. We generate our synthetic galaxies according to the galaxy stellar mass function. We estimate the metallicity according to both the mass-metallicity relation (MZR) and the fundamental metallicity relation (FMR). Also, we take into account galaxy-galaxy mergers and the evolution of the galaxy properties from the formation to the merger of the binary compact object. We find that the merger rate density changes dramatically depending on the choice of the star-forming galaxy main sequence, especially in the case of binary black holes (BBHs) and black hole neutron star systems (BHNSs). The slope of the merger rate density of BBHs and BHNSs is steeper if we assume the MZR with respect to the FMR, because the latter predicts a shallower decrease of metallicity with redshift. In contrast, binary neutron stars (BNSs) are only mildly affected by both the galaxy main sequence and metallicity relation. Overall, BBHs and BHNSs tend to form in low-mass metal-poor galaxies and merge in high-mass metal-rich galaxies, while BNSs form and merge in massive galaxies. We predict that passive galaxies host at least ~5-10%, ~15-25%, and ~15-35% of all BNS, BHNS and BBH mergers in the local Universe., Comment: 21 pages, 22 figures (including appendices), 3 tables, published in MNRAS

In this work, we combine ellipticity and major axis position angle measurements from the Sloan Digital Sky Server Data Release 16 (SDSS DR16) with the group finder algorithm of Rodriguez \& Merch\'an to determine the alignment of the central galaxies with the surrounding structures and satellite galaxies lying in their group. We use two independent methods: A modified version of the two-point cross-correlation function and the angle between the central galaxy orientation and the satellite galaxies relative position. The first method allows us to study the inner and outer regions of the cluster, while the second method provides information within the halos. Our results show that central galaxies present anysotropy in the correlation function up to $\sim 10 h^{-1}Mpc$, which becomes $\sim$10\% stronger for the brightest ones ($^{0.1}M_{r}<-21.5$). When we split the galaxy sample by colour, we find that red central galaxies are the main contributors to this anisotropy. We also show that this behaviour does not depend on the group mass or central galaxy ellipticity. Finally, our results are in agreement with previous findings, showing that the two-point cross-correlation function is a good tracer of the galaxy alignments using all galaxies and not only those of the group to which it belongs. In addition, this feature allows us to explore the behaviour of the alignment on larger scales., Comment: 11 pages, 10 figures, Approved for publication in MNRAS

Formalisms based on temporal logics interpreted over finite strict linear orders, known in the literature as finite traces, have been used for temporal specification in automated planning, process modelling, (runtime) verification and synthesis of programs, as well as in knowledge representation and reasoning. In this paper, we focus on first-order temporal logic on finite traces. We first investigate preservation of equivalences and satisfiability of formulas between finite and infinite traces, by providing a set of semantic and syntactic conditions to guarantee when the distinction between reasoning in the two cases can be blurred. Moreover, we show that the satisfiability problem on finite traces for several decidable fragments of first-order temporal logic is ExpSpace-complete, as in the infinite trace case, while it decreases to NExpTime when finite traces bounded in the number of instants are considered. This leads also to new complexity results for temporal description logics over finite traces. Finally, we investigate applications to planning and verification, in particular by establishing connections with the notions of insensitivity to infiniteness and safety from the literature.

Elucidating the connection between the properties of galaxies and the properties of their hosting haloes is a key element in galaxy formation. When the spatial distribution of objects is also taken under consideration, it becomes very relevant for cosmological measurements. In this paper, we use machine learning techniques to analyse these intricate relations in the IllustrisTNG300 magnetohydrodynamical simulation, predicting baryonic properties from halo properties. We employ four different algorithms: extremely randomized trees, K-nearest neighbours, light gradient boosting machine, and neural networks, along with a unique and powerful combination of the results from all four approaches. Overall, the different algorithms produce consistent results in terms of predicting galaxy properties from a set of input halo properties that include halo mass, concentration, spin, and halo overdensity. For stellar mass, the Pearson correlation coefficient is 0.98, dropping down to 0.7-0.8 for specific star formation rate (sSFR), colour, and size. In addition, we apply, for the first time in this context, an existing data augmentation method, synthetic minority over-sampling technique for regression with Gaussian noise (SMOGN), designed to alleviate the problem of imbalanced data sets, showing that it improves the overall shape of the predicted distributions and the scatter in the halo-galaxy relations. We also demonstrate that our predictions are good enough to reproduce the power spectra of multiple galaxy populations, defined in terms of stellar mass, sSFR, colour, and size with high accuracy. Our results align with previous reports suggesting that certain galaxy properties cannot be reproduced using halo features alone., Comment: Matches published version; very minor changes wrt V1

The detection of electromagnetic radiation (EM) accompanying the gravitational wave (GW) signal from the binary neutron star (BNS) merger GW170817 has revealed that these systems constitute at least a fraction of the progenitors of short gamma-ray bursts (SGRBs). As gravitational wave detectors keep pushing their detection horizons, it is important to assess coupled GW/EM probabilities, and how to maximize observational prospects. Here we perform population synthesis calculations of BNS evolution with the code MOBSE, and seed the binaries in galaxies at three representative redshifts (z=0.01,0.1,1) of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the astrophysical parameters of GRB170817A as a prototype, we numerically compute the broadband lightcurves of jets from BNS mergers, with the afterglow brightness depending on the local medium density at the merger sites. We perform Monte Carlo simulations of the resulting EM population assuming either a random viewing angle with respect to the jet, or a jet aligned with the orbital angular momentum of the binary, which biases the viewing angle probability for GW-triggered events. We find that ~70-80% of BNSs from z=0.01 should be detectable in gamma-rays. The afterglow detection probabilities of GW-triggered BNS mergers vary between ~0.3-0.7%, with higher values for jets aligned with the BNS angular momentum, and are comparable across the high and low-energy bands, unlike gamma-ray-triggered events (cosmological SGRBs) which are significantly brighter at higher energies. We further quantify observational biases with respect to host galaxy masses., Comment: 15 pages, 12 figures, submitted to MNRAS

Aiming at ontology-based data access to temporal data, we design two-dimensional temporal ontology and query languages by combining logics from the (extended) DL-Lite family with linear temporal logic LTL over discrete time (Z,<). Our main concern is first-order rewritability of ontology-mediated queries (OMQs) that consist of a 2D ontology and a positive temporal instance query. Our target languages for FO-rewritings are two-sorted FO(<) - first-order logic with sorts for time instants ordered by the built-in precedence relation < and for the domain of individuals - its extension FOE with the standard congruence predicates t \equiv 0 mod n, for any fixed n > 1, and FO(RPR) that admits relational primitive recursion. In terms of circuit complexity, FOE- and FO(RPR)-rewritability guarantee answering OMQs in uniform AC0 and NC1, respectively. We proceed in three steps. First, we define a hierarchy of 2D DL-Lite/LTL ontology languages and investigate the FO-rewritability of OMQs with atomic queries by constructing projections onto 1D LTL OMQs and employing recent results on the FO-rewritability of propositional LTL OMQs. As the projections involve deciding consistency of ontologies and data, we also consider the consistency problem for our languages. While the undecidability of consistency for 2D ontology languages with expressive Boolean role inclusions might be expected, we also show that, rather surprisingly, the restriction to Krom and Horn role inclusions leads to decidability (and ExpSpace-completeness), even if one admits full Booleans on concepts. As a final step, we lift some of the rewritability results for atomic OMQs to OMQs with expressive positive temporal instance queries. The lifting results are based on an in-depth study of the canonical models and only concern Horn ontologies., Comment: Accepted for JAIR

The environment within dark matter haloes can quench the star formation of galaxies. However, environmental effects beyond the virial radius of haloes ($\gtrsim$ 1 Mpc) are less evident. An example is the debated correlation between colour or star formation in central galaxies and neighbour galaxies in adjacent haloes at large separations of several Mpc, referred to as two-halo galactic conformity. We use two galaxy catalogues generated from different versions of the semi-analytic model SAG applied to the MDPL2 cosmological simulation and the IllustrisTNG300 cosmological hydrodynamical simulation to study the two-halo conformity by measuring the quenched fraction of neighbouring galaxies as a function of the real-space distance from central galaxies. We find that low-mass central galaxies in the vicinity of massive systems ($M_{\rm 200c}$ $\geq$ 10$^{13}$ $h^{-1}~\rm M_{\odot}$) out to 5 $h^{-1}$ Mpc are preferentially quenched compared to other central galaxies at fixed stellar mass $M_{\star}$ or fixed host halo mass $M_{\rm 200c}$ at $z$ ~ 0. In all the galaxies catalogues is consistent that the low-mass ($M_{\star} < 10^{10}$ $h^{-1}~\rm M_{\odot}$ or $M_{\rm 200c} < 10^{11.8}$ $h^{-1}~\rm M_{\odot}$) central galaxies in the vicinity of clusters and, especially, groups of galaxies mostly produce the two-halo galactic conformity. On average, the quenched low-mass central galaxies are much closer to massive haloes than star-forming central galaxies of the same mass (by a factor of ~5). Our results agree with other works regarding the environmental influence of massive haloes that can extend beyond the virial radius and affect nearby low-mass central galaxies., Comment: 15 pages and 10 figures without appendix. Accepted for publication in MNRAS

We study the clustering of Gravitational Wave (GW) merger events and Supernovae IA (SN), as cosmic tracers in Luminosity Distance Space. We modify the publicly available CAMB code to numerically evaluate auto- and cross- power spectra for the different sources, including Luminosity Distance Space distortion effects generated by peculiar velocities and lensing convergence. We perform a multitracer Fisher analysis to forecast expected constraints on cosmological and GW bias coefficients, using outputs from hydrodinamical N-body simulations to determine the bias fiducial model and considering future observations from the Vera Rubin Observatory and Einstein Telescope (ET), both single and in a 3 detector network configuration. We find that adding SN to the GW merger dataset considerably improves the forecast, mostly by breaking significant parameter degeneracies, with final constraints comparable to those obtainable from a Euclid-like survey. GW merger bias is forecasted to be detectable with good significance even in the single ET case., Comment: Version accepted for publication on JCAP, 25 pages, 10 figures

The growing population of binary black holes (BBHs) observed by gravitational wave detectors is a potential Rosetta stone for understanding their formation channels. Here, we use an upgraded version of our semi-analytic codes {\sc fastcluster} and {\sc cosmo$\mathcal{R}$ate} to investigate the cosmic evolution of four different BBH populations: isolated BBHs and dynamically formed BBHs in nuclear star clusters (NSCs), globular clusters (GCs), and young star clusters (YSCs). With our approach, we can study different channels assuming the same stellar and binary input physics. We find that the merger rate density of BBHs in GCs and NSCs is barely affected by stellar metallicity ($Z$), while the rate of isolated BBHs changes wildly with $Z$. BBHs in YSCs behave in an intermediate way between isolated and GC/NSC BBHs. The local merger rate density of Nth-generation black holes (BHs), obtained by summing up hierarchical mergers in GCs, NSCs and YSCs, ranges from $\sim{1}$ to $\sim{4}$ Gpc$^{-3}$ yr$^{-1}$ and is mostly sensitive to the spin parameter. We find that the mass function of primary BHs evolves with redshift in GCs and NSCs, becoming more top-heavy at higher $z$. In contrast, the primary BH mass function almost does not change with redshift in YSCs and in the field. This signature of the BH mass function has relevant implications for Einstein Telescope and Cosmic Explorer. Finally, our analysis suggests that multiple channels contribute to the BBH population of the second gravitational-wave transient catalog., Comment: 20 pages, 14 figures, 3 tables, accepted for publication in MNRAS

Definite descriptions are phrases of the form 'the $x$ such that $\varphi$', used to refer to single entities in a context. They are often more meaningful to users than individual names alone, in particular when modelling or querying data over ontologies. We investigate free description logics with both individual names and definite descriptions as terms of the language, while also accounting for their possible lack of denotation. We focus on the extensions of $\mathcal{ALC}$ and, respectively, $\mathcal{EL}$ with nominals, the universal role, and definite descriptions. We show that standard reasoning in these extensions is not harder than in the original languages, and we characterise the expressive power of concepts relative to first-order formulas using a suitable notion of bisimulation. Moreover, we lay the foundations for automated support for definite descriptions generation by studying the complexity of deciding the existence of definite descriptions for an individual under an ontology. Finally, we provide a polynomial-time reduction of reasoning in other free description logic languages based on dual-domain semantics to the case of partial interpretations.

Star formation in molecular clouds is clumpy, hierarchically subclustered. Fractal structure also emerges in hydro-dynamical simulations of star-forming clouds. Simulating the formation of realistic star clusters with hydro-dynamical simulations is a computational challenge, considering that only the statistically averaged results of large batches of simulations are reliable, due to the chaotic nature of the gravitational N-body problem. While large sets of initial conditions for N-body runs can be produced by hydro-dynamical simulations of star formation, this is prohibitively expensive in terms of computational time. Here we address this issue by introducing a new technique for generating many sets of new initial conditions from a given set of star masses, positions and velocities from a hydro-dynamical simulation. We use hierarchical clustering in phase space to learn a tree representation of the spatial and kinematic relations between stars. This constitutes the basis for the random generation of new sets of stars which share the same clustering structure of the original ones but have individually different masses, positions, and velocities. We apply this method to the output of a number of hydro-dynamical star-formation simulations, comparing the generated initial conditions to the original ones through a series of quantitative tests, including comparing mass and velocity distributions and fractal dimension. Finally, we evolve both the original and the generated star clusters using a direct N-body code, obtaining a qualitatively similar evolution., Comment: Updated version of the manuscript: "Seeing the forest for the trees: hierarchical generative models for star clusters from hydro-dynamical simulations". 15 pages, 15 figures, 2 tables. Comments welcome

Halo assembly bias is the secondary dependence of the clustering of dark-matter haloes on their assembly histories at fixed halo mass. This established dependence is expected to manifest itself on the clustering of the galaxy population, a potential effect commonly known as galaxy assembly bias. Using the IllustrisTNG300 magnetohydrodynamical simulation, we analyse the dependence of the properties and clustering of galaxies on the shape of the specific mass accretion history of their hosting haloes (sMAH). We first show that several halo and galaxy properties strongly correlate with the slope of the sMAH ($\beta$) at fixed halo mass. Namely, haloes with increasingly steeper $\beta$ increment their halo masses faster at early times, and their hosted galaxies present larger stellar-to-halo mass ratios, lose their gas faster, reach the peak of their star formation histories at higher redshift, and become quenched earlier. We also demonstrate that $\beta$ is more directly connected to these key galaxy formation properties than other broadly employed halo proxies, such as formation time. Finally, we measure the secondary dependence of galaxy clustering on $\beta$ at fixed halo mass as a function of redshift. By tracing back the evolution of individual haloes, we show that the amplitude of the galaxy assembly bias signal for the progenitors of $z=0$ galaxies increases with redshift, reaching a factor of 2 at $z = 1$ for haloes of $M_\mathrm{halo}=10^{11.5}-10^{12}$ $h^{-1}\mathrm{M}_\odot$. The measurement of the evolution of assembly bias along the merger tree provides a new theoretical perspective to the study of secondary bias. Our findings, which show a tight relationship between halo accretion and both the clustering and the observational properties of the galaxy population, have also important implications for the generation of mock catalogues for upcoming cosmological surveys., Comment: 10 pages, 6 figures. Submitted to MNRAS

We explore hierarchical black hole (BH) mergers in nuclear star clusters (NSCs), globular clusters (GCs) and young star clusters (YSCs), accounting for both original and dynamically assembled binary BHs (BBHs). We find that the median mass of both first- and nth-generation dynamical mergers is larger in GCs and YSCs with respect to NSCs, because the lighter BHs are ejected by supernova kicks from the lower-mass clusters. Also, first- and nth-generation BH masses are strongly affected by the metallicity of the progenitor stars: the median mass of the primary BH of a nth-generation merger is $\sim{}24-38$ M$_\odot$ ($\sim{}9-15$ M$_\odot$) in metal-poor (metal-rich) NSCs. The maximum BH mass mainly depends on the escape velocity: BHs with mass up to several thousand M$_\odot$ form in NSCs, while YSCs and GCs host BHs with mass up to several hundred M$_\odot$. Furthermore, we calculate the fraction of mergers with at least one component in the pair-instability mass gap ($f_{\rm PI}$) and in the intermediate-mass BH regime ($f_{\rm IMBH}$). In the fiducial model for dynamical BBHs with metallicity $Z=0.002$, we find $f_{\rm PI}\approx{}0.05$, $0.02$ and $0.007$ ($f_{\rm IMBH}\approx{}0.01$, $0.002$ and $0.001$) in NSCs, GCs and YSCs, respectively. Both $f_{\rm PI}$ and $f_{\rm IMBH}$ drop by at least one order of magnitude at solar metallicity. Finally, we investigate the formation of GW190521 by assuming that it is either a nearly equal-mass BBH or an intermediate-mass ratio inspiral., Comment: 20 pages, 13 figures, accepted for publication in MNRAS

With the recent release of the second gravitational-wave transient catalogue (GWTC-2), which introduced dozens of new detections, we are at a turning point of gravitational wave astronomy, as we are now able to directly infer constraints on the astrophysical population of compact objects. Here, we tackle the burning issue of understanding the origin of binary black hole (BBH) mergers. To this effect, we make use of state-of-the-art population synthesis and N-body simulations, to represent two distinct formation channels: BBHs formed in the field (isolated channel) and in young star clusters (dynamical channel). We then use a Bayesian hierarchical approach to infer the distribution of the mixing fraction $f$, with $f=0$ ($f=1$) in the pure dynamical (isolated) channel. %that controls the proportion of isolated and dynamical BBHs. We explore the effects of additional hyper-parameters of the model, such as the spread in metallicity $\sigma_{\text{Z}}$ and the parameter $\sigma_{\text{sp}}$, describing the distribution of spin magnitudes. We find that the dynamical model is slightly favoured with a median value of $f=0.26$, when $\sigma_{\text{sp}}=0.1$ and $\sigma_{\text{Z}}=0.4$. Models with higher spin magnitudes tend to strongly favour dynamically formed BBHs ($f\le{}0.1$ if $\sigma_{\text{sp}}=0.3$). Furthermore, we show that hyper-parameters controlling the rates of the model, such as $\sigma_{\rm Z}$, have a large impact on the inference of the mixing fraction, which rises from $0.18$ to $0.43$ when we increase $\sigma_{\text{Z}}$ from 0.2 to 0.6, for a fixed value of $\sigma_{\text{sp}}=0.1$. Finally, our current set of observations is better described by a combination of both formation channels, as a pure dynamical scenario is excluded at the $99\%$ credible interval, except when the spin magnitude is high., Comment: 13 pages, 10 figures, 2 tables, published in MNRAS

We study the intrinsic large-scale distribution and evolution of seven ionized metals in the IllustrisTNG magneto-hydrodynamical cosmological simulation. We focus on the fractions of C\,\textsc{ii}, C\,\textsc{iv}, Mg\,\textsc{ii}, N\,\textsc{v}, Ne\,\textsc{viii}, O\,\textsc{vi}, and Si\,\textsc{iv} in different cosmic web structures (filaments, haloes, and voids) and gas phases (warm-hot intergalactic medium WHIM, hot, diffuse, and condensed gas) from $z=6$ to $z=0$. Our analysis provides a new perspective to the study of the distribution and evolution of baryons across cosmic time while offering new hints in the context of the well-known missing baryons problem. The cosmic web components are here identified using the local comoving dark matter density, which provides a simple but effective way of mapping baryons on large scales. Our results show that C\,\textsc{ii} and Mg\,\textsc{ii} are mostly located in condensed gas inside haloes in high-density and low-temperature star-forming regions ($\rho_{\rm gas}/\bar{\rho}_{\rm bar}\gtrsim10^3$, and ${\rm T}\lesssim10^{5}$~K). C\,\textsc{iv} and Si\,\textsc{iv} present similar evolution of their mass fractions in haloes and filaments across cosmic time. In particular, their mass budgets in haloes in condensed phase ($\rho_{\rm gas}/\bar{\rho}_{\rm bar}\gtrsim10^3$, and ${\rm T}\lesssim10^{5}$~K) are driven by gas cooling and star formation with a peak at $z\sim2$. Finally, our results confirm that O\,\textsc{vi}, Ne\,\textsc{viii}, and N\,\textsc{v} are good tracers of warm/hot and low-density gas at low redshift ($\rho_{\rm gas}/\bar{\rho}_{\rm bar}\lesssim10^3$, and ${\rm T}\gtrsim10^{5}$~K), regions that are likely to contain most of the missing baryons in the local Universe., Comment: 15 pages, 9 figures. Submitted to MNRAS

We use the IllustrisTNG100 hydrodynamical simulation to study the dependence of the galaxy two-point correlation function on a broad range of secondary subhalo and galactic properties. We construct galaxy mock catalogues adopting a standard sub-halo abundance matching scheme coupled with a secondary assignment between galaxy colour or specific star formation rate and the following subhalo properties: starvation redshift z$_{\rm starve}$, concentration at infall, overdensity $\delta_R^{\rm env}$, tidal anisotropy $\alpha_R$, and tidal overdensity $\delta_R$. The last two quantities allow us to fully characterise the tidal field of our subhaloes, acting as mediators between their internal and large-scale properties. The resulting mock catalogues overall return good agreement with the IllustrisTNG100 measurements. The accuracy of each model strongly depends on the correlation between the secondary galaxy and subhalo properties employed. Among all the subhalo proxies tested, we find that z$_{\rm starve}$ and $c_{\rm infall}$ are the ones that best trace the large-scale structure, producing robust clustering predictions for different samples of red/blue and quenched/star-forming galaxies., Comment: 12 pages, 8 figures

In this work, we combine size and stellar mass measurements from the Sloan Digital Sky Server (SDSS) with the group finder algorithm of Rodriguez \& Merch\'an in order to determine the stellar and halo mass -- size relations of central and satellite galaxies separately. We show that, while central and satellite galaxies display similar stellar mass -- size relations, their halo mass -- size relations differ significantly. As expected, more massive haloes tend to host larger central galaxies. However, the size of satellite galaxies depends only slightly on halo virial mass. We show that these results are compatible with a remarkably simple model in which the size of central and satellite galaxies scales as the cubic root of their host halo mass, with the normalization for satellites being $\sim$ 30 \% smaller than that for central galaxies, which can be attributed to tidal stripping. We further check that our measurements are in excellent agreement with predictions from the IllustrisTNG hydrodynamical simulation. In the second part of this paper, we analyse how the clustering properties of central and satellite galaxies depend on their size. We demonstrate that, independently of the stellar mass threshold adopted, smaller galaxies are more tightly clustered than larger galaxies when either the entire sample or only satellites are considered. The opposite trend is observed on large scales when the size split is performed for the central galaxies alone. Our results place significant constraints for halo-galaxy connection models that link galaxy size with the properties of their hosting haloes., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS

We evaluate the redshift distribution of binary black hole (BBH), black hole - neutron star binary (BHNS) and binary neutron star (BNS) mergers, exploring the main sources of uncertainty: star formation rate (SFR) density, metallicity evolution, common envelope, mass transfer via Roche lobe overflow, natal kicks, core-collapse supernova model and initial mass function. Among binary evolution processes, uncertainties on common envelope ejection have a major impact: the local merger rate density of BNSs varies from $\sim{}10^3$ to $\sim{}20$ Gpc$^{-3}$ yr$^{-1}$ if we change the common envelope efficiency parameter from $\alpha_{\rm CE}=7$ to 0.5, while the local merger rates of BBHs and BHNSs vary by a factor of $\sim{}2-3$. The BBH merger rate changes by one order of magnitude, when $1 \sigma$ uncertainties on metallicity evolution are taken into account. In contrast, the BNS merger rate is almost insensitive to metallicity. Hence, BNSs are the ideal test bed to put constraints on uncertain binary evolution processes, such as common envelope and natal kicks. Only models assuming values of $\alpha_{\rm CE}\gtrsim{}2$ and moderately low natal kicks (depending on the ejected mass and the SN mechanism), result in a local BNS merger rate density within the 90% credible interval inferred from the second gravitational-wave transient catalogue., Comment: 14 pages, 12 figures, 2 tables, accepted for publication in MNRAS

We explore the potential of the kinetic Sunyaev-Zel'dovich (kSZ) effect as the cornerstone of a future observational probe for halo spin bias, the secondary dependence of halo clustering on halo spin at fixed halo mass. Using the IllustrisTNG magneto-hydrodynamical cosmological simulation, we measure both the kSZ and the thermal SZ (tSZ) effects produced by the baryonic content of more than 50,000 haloes within the halo mass range $11 < \log_{10} ({\rm M_{vir}}/ h^{-1} {\rm M_{\odot}}) \lesssim 14.5$. First, we confirm that the magnitude of both effects depends strongly on the total gas and virial mass of the haloes, and that the integrated kSZ signal displays a significant correlation with the angular momentum of the intra-halo gas, particularly for massive haloes. Second, we show that both the integrated kSZ signal and the ratio of the integrated kSZ and tSZ signals trace total halo spin, even though significant scatter exists. Finally, we demonstrate that, in the absence of observational and instrumental uncertainties, these SZ-related statistics can be used to recover most of the underlying IllustrisTNG halo spin bias signal. Our analysis represents the first attempt to develop a future observational probe for halo spin bias, bringing forward alternative routes for measuring the secondary bias effects., Comment: 17 pages, 10 figures, submitted to MNRAS

This paper investigates the feasibility of automated reasoning over temporal DL-Lite (TDL-Lite) knowledge bases (KBs). We test the usage of off-the-shelf LTL reasoners to check satisfiability of TDL-Lite KBs. In particular, we test the robustness and the scalability of reasoners when dealing with TDL-Lite TBoxes paired with a temporal ABox. We conduct various experiments to analyse the performance of different reasoners by randomly generating TDL-Lite KBs and then measuring the running time and the size of the translations. Furthermore, in an effort to make the usage of TDL-Lite KBs a reality, we present a fully fledged tool with a graphical interface to design them. Our interface is based on conceptual modelling principles and it is integrated with our translation tool and a temporal reasoner.

In this work we aim at investigating the effects of baryons on the dark matter (DM) haloes structure, focusing on the correlation between the presence and importance of stellar discs and the halo shapes. We study the properties of a subsample of DM haloes from Fenix and EAGLE cosmological simulations. We inspect the central regions of haloes in the mass range [10.9 - 992.3] x 10^{10} Msun at z=0, comparing the hydrodynamic runs and their dark matter only (DMo) counterparts. Our results indicate that baryons have a significant impact on the shape of the inner halo, mainly within ~ 20 percent of the virial radius. We find haloes to be more spherical when hosting baryons. While the impact of baryons depends on the mass of the haloes, we also find a trend with morphology which suggests that the way baryons are assembled is also relevant in agreement with previous works. Our findings also indicate that disc galaxies preferentially form in haloes whose DMo counterparts were originally more spherical and with stronger velocity anisotropy. The presence of baryons alter the orbital structure of the DM particles of the haloes, which show a decrease in their velocity anisotropy, towards more tangentially biased orbits. This relative decrease is weaker in the case of disc-dominated galaxies. Our results point out to a cosmological connection between the final morphology of galaxies and the intrinsic properties of their DM haloes, which gets reinforce by the growth of the discs., Comment: 13 pages, 13 figures, published after revision in MNRAS

Hierarchical mergers are one of the distinctive signatures of binary black hole (BBH) formation through dynamical evolution. Here, we present a fast semi-analytic approach to simulate hierarchical mergers in nuclear star clusters (NSCs), globular clusters (GCs) and young star clusters (YSCs). Hierarchical mergers are more common in NSCs than they are in both GCs and YSCs, because of the different escape velocity. The mass distribution of hierarchical BBHs strongly depends on the properties of first-generation BBHs, such as their progenitor's metallicity. In our fiducial model, we form black holes (BHs) with masses up to $\sim{}10^3$ M$_\odot$ in NSCs and up to $\sim{}10^2$ M$_\odot$ in both GCs and YSCs. When escape velocities in excess of 100 km~s$^{-1}$ are considered, BHs with mass $>10^3$ M$_\odot$ are allowed to form in NSCs. Hierarchical mergers lead to the formation of BHs in the pair instability mass gap and intermediate-mass BHs, but only in metal-poor environments. The local BBH merger rate in our models ranges from $\sim{}10$ to $\sim{} 60$ Gpc$^{-3}$ yr$^{-1}$; hierarchical BBHs in NSCs account for $\sim{}10^{-2}- 0.2$ Gpc$^{-3}$ yr$^{-1}$, with a strong upper limit of $\sim{}10$ Gpc$^{-3}$ yr$^{-1}$. When comparing our models with the second gravitational-wave transient catalog, we find that multiple formation channels are favored to reproduce the observed BBH population., Comment: 36 pages, 12 figures, 2 tables, accepted for publication in Symmetry (MDPI), special issue "Probing the Universe with Gravitational Waves", editors: Jos\'e A. de Freitas Pacheco and Pasquale Bosso