The OpenUniverse2024 simulation suite is a cross-collaboration effort to produce matched simulated imaging for multiple surveys as they would observe a common simulated sky. Both the simulated data and associated tools used to produce it are intended to uniquely enable a wide range of studies to maximize the science potential of the next generation of cosmological surveys. We have produced simulated imaging for approximately 70 deg^2 of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Wide-Fast-Deep survey and the Nancy Grace Roman Space Telescope High-Latitude Wide-Area Survey, as well as overlapping versions of the ELAIS-S1 Deep-Drilling Field for LSST and the High-Latitude Time-Domain Survey for Roman. OpenUniverse2024 includes i) an early version of the updated extragalactic model called Diffsky, which substantially improves the realism of optical and infrared photometry of objects, compared to previous versions of these models; ii) updated transient models that extend through the wavelength range probed by Roman and Rubin; and iii) improved survey, telescope, and instrument realism based on up-to-date survey plans and known properties of the instruments. It is built on a new and updated suite of simulation tools that improves the ease of consistently simulating multiple observatories viewing the same sky. The approximately 400 TB of synthetic survey imaging and simulated universe catalogs are publicly available, and we preview some scientific uses of the simulations.

The shapes of dark matter haloes are sensitive to both cosmology and baryon physics, but are difficult to measure observationally. A promising way to constrain them is to use the positions of satellite galaxies as tracers of the underlying dark matter, but there are typically too few galaxies per halo for reliable shape estimates, resulting in biased shapes. We present a method to model sampling noise to correct for the shape bias. We compare our predicted median shape bias with that obtained from the FLAMINGO suite of simulations and find reasonable agreement. We check that our results are robust to resolution effects and baryonic feedback. We also explore the validity of our bias correction at various redshifts and we discuss how our method might be applied to observations in the future. We show that median projected halo axis ratios are on average biased low by 0.31 when they are traced by only 5 satellites. Using the satellite galaxies, the projected host halo axis ratio can be corrected with a residual bias of 0.1, by accounting for sampling bias. Hence, about two-thirds of the projected axis ratio bias can be explained by sampling noise. This enables the statistical measurement of halo shapes at lower masses than previously possible. Our method will also allow improved estimates of halo shapes in cosmological simulations using fewer particles than currently required.

ABSTRACT The high abundance of disc galaxies without a large central bulge challenges predictions of current hydrodynamic simulations of galaxy formation. We aim to shed light on the formation of these objects by studying the redshift and mass dependence of their intrinsic 3D shape distributions in the COSMOS galaxy survey below redshift z = 1.0. This distribution is inferred from the observed distribution of 2D shapes, using a reconstruction method which we test using hydrodynamic simulations. Our tests reveal a moderate bias for the inferred average disc circularity and relative thickness, but a large bias on the dispersion of these quantities. Applying the reconstruction method on COSMOS data, we find variations of the average disc circularity and relative thickness with redshift of around ∼1 per cent and ∼10 per cent, respectively, which is comparable to the error estimates on these quantities. The average relative disc thickness shows a significant mass dependence which can be accounted for by the scaling of disc radius with galaxy mass. We conclude that our data provides no evidence for a strong dependence of the average circularity and absolute thickness of disc-dominated galaxies on redshift and mass that is significant with respect to the statistical uncertainties in our analysis. These findings are expected in the absence of disruptive merging or feedback events that would affect galaxy shapes. They hence support a scenario where present-day discs form early ( z > 1.0) and subsequently undergo a tranquil evolution in isolation. However, more data and a better understanding of systematics are needed to reaffirm our results.

We present a bright galaxy sample with accurate and precise photometric redshifts (photo-zs), selected using ugriZYJHK(s) photometry from the Kilo-Degree Survey (KiDS) Data Release 4. The highly pure and complete dataset is flux-limited at r < 20 mag, covers similar to 1000 deg(2), and contains about 1 million galaxies after artifact masking. We exploit the overlap with Galaxy And Mass Assembly spectroscopy as calibration to determine photo-zs with the supervised machine learning neural network algorithm implemented in the ANNz2 software. The photo-zs have a mean error of vertical bar delta z vertical bar similar to 5 x 10(-4) and low scatter (scaled mean absolute deviation of similar to 0:018(1 + z)); they are both practically independent of the r-band magnitude and photo-z at 0:05 < z(phot) < 0:5. Combined with the 9-band photometry, these allow us to estimate robust absolute magnitudes and stellar masses for the full sample. As a demonstration of the usefulness of these data, we split the dataset into red and blue galaxies, used them as lenses, and measured the weak gravitational lensing signal around them for five stellar mass bins. We fit a halo model to these high-precision measurements to constrain the stellar-mass-halo-mass relations for blue and red galaxies. We find that for high stellar mass (M-* > 5 x 10(11) M-circle dot), the red galaxies occupy dark matter halos that are much more massive than those occupied by blue galaxies with the same stellar mass.

Context. Feedback processes play a fundamental role in the regulation of the star formation (SF) activity in galaxies and, in particular, in the quenching of early-type galaxies (ETGs) as has been inferred by observational and numerical studies of Λ-CDM models. At z = 0, ETGs exhibit well-known fundamental scaling relations, but the connection between scaling relations and the physical processes shaping ETG evolution remains unknown. Aims. This work aims to study the impact of the energetic feedback due to active galactic nuclei (AGN) on the formation and evolution of ETGs. We focus on assessing the impact of AGN feedback on the evolution of the mass–plane and the fundamental plane (FP; defined using mass surface density) as well as on morphology, kinematics, and stellar age across the FP. Methods. The Horizon-AGN and Horizon-noAGN cosmological hydrodynamical simulations were performed with identical initial conditions, including the same physical processes except for the activation of the AGN feedback in the former. We selected a sample of central ETGs from both simulations using the same criteria and exhaustively studied their SF activity, kinematics, and scaling relations for z ≤ 3. Results. We find that Horizon-AGN ETGs identified at z = 0 follow the observed fundamental scaling relations (mass–plane, FP, and mass–size relation) and qualitatively reproduce kinematic features albeit conserving a rotational inner component with a mass fraction regulated by the AGN feedback. We discover that AGN feedback seems to be required to reproduce the bimodality in the spin parameter distribution reported by observational works and the mass–size relation; more massive galaxies have older stellar populations, larger sizes, and are slower rotators. We study the evolution of the fundamental relations with redshift, finding a mild evolution of the mass–plane of Horizon-AGN ETGs for z < 1, whereas a stronger change is detected for z > 1. The ETGs in Horizon-noAGN show a strong systematic redshift evolution of the mass–plane. The FP of Horizon-AGN ETGs agrees with observations at z = 0. When AGN feedback is switched off, a fraction of galaxies depart from the expected FP at all analysed redshifts owing to the presence of a few extended galaxies with an excess of stellar surface density. We find that AGN feedback regulates the SF activity as a function of stellar mass and redshift being able to reproduce the observed relations. Our results show the impact of AGN feedback on the mass-to-light ratio (M/L) and its relation with the tilt of the luminosity FP (L-FP; defined using the averaged surface brightness). Overall, AGN feedback has an impact on the regulation of the SF activity, size, stellar surface density, stellar ages, rotation, and masses of ETGs that is reflected on the fundamental relations, particularly on the FP. We detect a dependence of the FP on stellar age and galaxy morphology that evolves with redshfit. The characteristics of the galaxy distribution on the FP according to these properties change drastically by z ∼ 1 in Horizon-AGN and hence this feature could provide further insight into the action of AGN feedback.

Studies of intrinsic alignment effects mostly focus on the correlations between shapes of galaxies with each other, or with the underlying density field of the large scale structure of the universe. Lately, the correlation between shapes of galaxies and the large-scale velocity field has been proposed as an additional probe of the large scale structure. We use a Fisher forecast to make a prediction for the detectability of this velocity-shape correlation with a combination of redshifts and shapes from the 4MOST+LSST surveys, and radial velocity reconstruction from the Simons Observatory. The signal-to-noise ratio for the velocity-shape (dipole) correlation is 23, relative to 44 for the galaxy density-shape (monopole) correlation and for a maximum wavenumber of $0.2\: \mathrm{Mpc^{-1}}$. Increasing signal-to-noise for higher values of the maximum wavenumbers (respectively, 56 and 69, for a maximum wavenumber of $1\: \mathrm{Mpc^{-1}}$) indicate potential gains in the nonlinear regime. Encouraged by these predictions, we discuss two possible applications for the velocity-shape correlation. Measuring the velocity-shape correlation could improve the mitigation of selection effects induced by intrinsic alignments on galaxy clustering. We also find that velocity-shape measurements could potentially aid in determining the scale-dependence of intrinsic alignments when multiple shape measurements of the same galaxies are provided.

Galaxies and clusters embedded in the large-scale structure of the Universe are observed to align in preferential directions. Galaxy alignment has been established as a potential probe for cosmological information, but the application of cluster alignments for these purposes remains unexplored. Clusters are observed to have a higher alignment amplitude than galaxies, but because galaxies are much more numerous, the trade-off in detectability between the two signals remains unclear. We present forecasts comparing cluster and galaxy alignments for two extragalactic survey set-ups: a currently-available low redshift survey (SDSS) and an upcoming higher redshift survey (LSST). For SDSS, we rely on the publicly available redMaPPer catalogue to describe the cluster sample. For LSST, we perform estimations of the expected number counts while we extrapolate the alignment measurements from SDSS. Clusters in SDSS have typically higher alignment signal-to-noise than galaxies. For LSST, the cluster alignment signals quickly wash out with redshift due to a relatively low number count and a decreasing alignment amplitude. Nevertheless, a potential strong-suit of clusters is in their interplay with weak lensing: intrinsic alignments can be more easily isolated for clusters than for galaxies. The signal-to-noise of cluster alignment can in general be improved by isolating close pairs along the line of sight.

We present a bright galaxy sample with accurate and precise photometric redshifts (photo-zs), selected usingugriZYJHKsphotometry from the Kilo-Degree Survey (KiDS) Data Release 4. The highly pure and complete dataset is flux-limited atr < 20 mag, covers ∼1000 deg2, and contains about 1 million galaxies after artifact masking. We exploit the overlap with Galaxy And Mass Assembly spectroscopy as calibration to determine photo-zs with the supervised machine learning neural network algorithm implemented in the ANNz2 software. The photo-zs have a mean error of |⟨δz⟩|∼5 × 10−4and low scatter (scaled mean absolute deviation of ∼0.018(1 + z)); they are both practically independent of ther-band magnitude and photo-zat 0.05 < zphot < 0.5. Combined with the 9-band photometry, these allow us to estimate robust absolute magnitudes and stellar masses for the full sample. As a demonstration of the usefulness of these data, we split the dataset into red and blue galaxies, used them as lenses, and measured the weak gravitational lensing signal around them for five stellar mass bins. We fit a halo model to these high-precision measurements to constrain the stellar-mass–halo-mass relations for blue and red galaxies. We find that for high stellar mass (M⋆ > 5 × 1011 M⊙), the red galaxies occupy dark matter halos that are much more massive than those occupied by blue galaxies with the same stellar mass.

We present a combined tomographic weak gravitational lensing analysis of the Kilo Degree Survey (KV450) and the Dark Energy Survey (DES-Y1). We homogenize the analysis of these two public cosmic shear datasets by adopting consistent priors and modeling of nonlinear scales, and determine new redshift distributions for DES-Y1 based on deep public spectroscopic surveys. Adopting these revised redshifts results in a $0.8\sigma$ reduction in the DES-inferred value for $S_8$, which decreases to a $0.5\sigma$ reduction when including a systematic redshift calibration error model from mock DES data based on the MICE2 simulation. The combined KV450 + DES-Y1 constraint on $S_8 = 0.762^{+0.025}_{-0.024}$ is in tension with the Planck 2018 constraint from the cosmic microwave background at the level of $2.5\sigma$. This result highlights the importance of developing methods to provide accurate redshift calibration for current and future weak lensing surveys.

Context . Accurate model predictions including the physics of baryons are required to make the most of the upcoming large cosmological surveys devoted to gravitational lensing. The advent of hydrodynamical cosmological simulations enables such predictions on sufficiently sizeable volumes. Aims . Lensing quantities (deflection, shear, convergence) and their statistics (convergence power spectrum, shear correlation functions, galaxy-galaxy lensing) are computed in the past lightcone built in the Horizon-AGN hydrodynamical cosmological simulation, which implements our best knowledge on baryonic physics at the galaxy scale in order to mimic galaxy populations over cosmic time. Methods . Lensing quantities are generated over a one square degree field of view by performing multiple-lens plane ray-tracing through the lightcone, taking full advantage of the 1 kpc resolution and splitting the line of sight over 500 planes all the way to redshift z ∼ 7. Two methods are explored (standard projection of particles with adaptive smoothing, and integration of the acceleration field) to ensure a good implementation. The focus is on small scales where baryons matter most. Results . Standard cosmic shear statistics are affected at the 10% level by the baryonic component for angular scales below a few arcminutes. The galaxy-galaxy lensing signal, or galaxy-shear correlation function, is consistent with measurements for the redshift z ∼ 0.5 massive galaxy population. At higher redshift z ≳ 1, the effect of magnification bias on this correlation is relevant for separations greater than 1 Mpc. Conclusions . This work is pivotal for all current and upcoming weak-lensing surveys and represents a first step towards building a full end-to-end generation of lensed mock images from large cosmological hydrodynamical simulations.

Context. Accurate model predictions including the physics of baryons are required to make the most of the upcoming large cosmological surveys devoted to gravitational lensing. The advent of hydrodynamical cosmological simulations enables such predictions on sufficiently sizeable volumes. Aims. Lensing quantities (deflection, shear, convergence) and their statistics (convergence power spectrum, shear correlation functions, galaxy-galaxy lensing) are computed in the past lightcone built in the Horizon-AGN hydrodynamical cosmological simulation, which implements our best knowledge on baryonic physics at the galaxy scale in order to mimic galaxy populations over cosmic time. Methods. Lensing quantities are generated over a one square degree field of view by performing multiple-lens plane ray-tracing through the lightcone, taking full advantage of the 1 kpc resolution and splitting the line of sight over 500 planes all the way to red-shift z similar to 7. Two methods are explored (standard projection of particles with adaptive smoothing, and integration of the acceleration field) to ensure a good implementation. The focus is on small scales where baryons matter most. Results. Standard cosmic shear statistics are affected at the 10% level by the baryonic component for angular scales below a few arcminutes. The galaxy-galaxy lensing signal, or galaxy-shear correlation function, is consistent with measurements for the redshift z greater than or similar to 0.5 massive galaxy population. At higher redshift z greater than or similar to 1, the effect of magnification bias on this correlation is relevant for separations greater than 1 Mpc. Conclusions. This work is pivotal for all current and upcoming weak-lensing surveys and represents a first step towards building a full end-to-end generation of lensed mock images from large cosmological hydrodynamical simulations.

This work investigates the alignment of galactic spins with the cosmic web across cosmic time using the cosmological hydrodynamical simulation Horizon-AGN. The cosmic web structure is extracted via the persistent skeleton as implemented in the DISPERSE algorithm. It is found that the spin of low-mass galaxies is more likely to be aligned with the filaments of the cosmic web and to lie within the plane of the walls while more massive galaxies tend to have a spin perpendicular to the axis of the filaments and to the walls. The mass transition is detected with a significance of 9σ. This galactic alignment is consistent with the alignment of the spin of dark haloes found in pure dark matter simulations and with predictions from (anisotropic) tidal torque theory. However, unlike haloes, the alignment of low-mass galaxies is weak and disappears at low redshifts while the orthogonal spin orientation of massive galaxies is strong and increases with time, probably as a result of mergers. At fixed mass, alignments are correlated with galaxy morphology: the high-redshift alignment is dominated by spiral galaxies while elliptical centrals are mainly responsible for the perpendicular signal. The two regimes probed in this work induce competing galactic alignment signals for weak lensing, with opposite redshift and luminosity evolution. Understanding the details of these intrinsic alignments will be key to exploit future major cosmic shear surveys like Euclid or LSST.

Using the cosmological hydrodynamics simulation Horizon-AGN, we investigate the spatial distribution of satellite galaxies relative to their central counterpart in the redshift range between 0.3 and 0.8. We find that, on average, these satellites tend to be located on the galactic plane of the central object. This effect is detected for central galaxies with a stellar mass larger than 10^10 solar masses and found to be strongest for red passive galaxies, while blue galaxies exhibit a weaker trend. For galaxies with a minor axis parallel to the direction of the nearest filament, we find that the coplanarity is stronger in the vicinity of the central galaxy, and decreases when moving towards the outskirts of the host halo. By contrast, the spatial distribution of satellite galaxies relative to their closest filament follows the opposite trend: their tendency to align with them dominates at large distances from the central galaxy, and fades away in its vicinity. Relying on mock catalogs of galaxies in that redshift range, we show that massive red centrals with a spin perpendicular to their filament also have corotating satellites well aligned with both the galactic plane and the filament. On the other hand, lower-mass blue centrals with a spin parallel to their filament have satellites flowing straight along this filament, and hence orthogonally to their galactic plane. The orbit of these satellites is then progressively bent towards a better alignment with the galactic plane as they penetrate the central region of their host halo. The kinematics previously described are consistent with satellite infall and spin build-up via quasi-polar flows, followed by a re-orientation of the spin of massive red galaxies through mergers.

Accurate cosmology from upcoming weak lensing surveys relies on knowledge of the total matter power spectrum at per cent level at scales k < 10 hMpc(-1), for which modelling the impact of baryonic physics is crucial. We compare measurements of the total matter power spectrum from the Horizon cosmological hydrodynamical simulations: a dark-matteronly run, one with full baryonic physics, and another lacking active galactic nucleus (AGN) feedback. Baryons cause a suppression of power at k similar or equal to 10 h Mpc(-1) of <15 per cent at z = 0, and an enhancement of a factor of a few at smaller scales due to the more efficient cooling and star formation. The results are sensitive to the presence of the highest mass haloes in the simulation and the distribution of dark matter is also impacted up to a few per cent. The redshift evolution of the effect is non-monotonic throughout z = 0-5 due to an interplay between AGN feedback and gas pressure, and the growth of structure. We investigate the effectiveness of an analytic `baryonic correction model' in describing our results. We require a different redshift evolution and propose an alternative fitting function with four free parameters that reproduces our results within 5 per cent. Compared to other simulations, we find the impact of baryonic processes on the total matter power spectrum to be smaller at z = 0. Correspondingly, our results suggest that AGN feedback is not strong enough in the simulation. Total matter power spectra from the Horizon simulations are made publicly available at https://www.horizon-simulation.org/catalogues.html.

The observed massive end of the galaxy stellar mass function is steeper than its predicted dark matter halo counterpart in the standard $\Lambda $CDM paradigm. In this paper, we investigate the impact of active galactic nuclei (AGN) feedback on star formation in massive galaxies. We isolate the impact of AGNs by comparing two simulations from the HORIZON suite, which are identical except that one also includes super massive black holes (SMBH), and related feedback models. This allows us to cross-identify individual galaxies between simulations and quantify the effect of AGN feedback on their properties, including stellar mass and gas outflows. We find that massive galaxies ($ \rm M_{*} \geq 10^{11} M_\odot $) are quenched by AGN feedback to the extent that their stellar masses decrease by up to 80% at $z=0$. SMBHs affect their host halo through a combination of outflows that reduce their baryonic mass, particularly for galaxies in the mass range $ \rm 10^9 M_\odot \leq M_{*} \leq 10^{11} M_\odot $, and a disruption of central gas inflows, which limits in-situ star formation. As a result, net gas inflows onto massive galaxies, $ \rm M_{*} \geq 10^{11} M_\odot $, drop by up to 70%. We measure a redshift evolution in the stellar mass ratio of twin galaxies with and without AGN feedback, with galaxies of a given stellar mass showing stronger signs of quenching earlier on. This evolution is driven by a progressive flattening of the $\rm M_{\rm SMBH}-M_* $ relation with redshift, particularly for galaxies with $\rm M_{*} \leq 10^{10} M_\odot $. $\rm M_{\rm SMBH}/M_*$ ratios decrease over time, as falling average gas densities in galaxies curb SMBH growth.

The HSC collaboration team: S. Abe, H. Aihara, M. Akiyama, K. Aoki, N. Arimoto, N. A. Bahcall, S. J. Bickerton, J. Bosch, K. Bundy†(3), C. W. Chen, M. Chiba†(4), T. Chiba, N. E. Chisari, J. Coupon, M. Doi, M. Enoki S. Foucaud, M. Fukugita, H. Furusawa†(5), T. Futamase, R. Goto, T. Goto, J. E. Greene, J. E. Gunn†(6), T. Hamana†(5), T. Hashimoto, M. Hayashi, Y. Higuchi, C. Hikage, J. C. Hill, P. T. P. Ho, B. C. Hsieh, K. Y. Huang†(7), H. Ikeda, M. Imanishi, N. Inada, A. K. Inoue, W.-H. Ip, T. Ito, K. Iwasawa, M. Iye, H. Y. Jian, Y. Kakazu, H. Karoji, N. Kashikawa, N. Katayama, T. Kawaguchi, S. Kawanomoto, I. Kayo, T. Kitayama, G. R. Knapp, T. Kodama, K. Kohno, M. Koike, E. Kokubo, M. Kokubo, Y. Komiyama, A. Konno, Y. Koyama, C. N. Lackner, D. Lang, A. Leauthaud†(3), M. J. Lehner, K.Y. Lin, L. Lin, Y.-T. Lin†(7), C. P. Loomis, R. H. Lupton†(6), P. S. Lykawka, K. Maeda, R. Mandelbaum†(22), Y. Matsuda, K. Matsuoka, Y. Matsuoka, S. Mineo, T. Minezaki, H. Miyatake, R. Momose, A. More, S. More, T. J. Moriya, T. Morokuma†(2), H. Murayama, K. Nagamine, T. Nagao†(23), S. Nagataki, Y. Naito, K. Nakajima, F. Nakata, H. Nakaya, T. Namikawa, C.-C. Ngeow, T. Nishimichi, H. Nishioka, A. J. Nishizawa†(3), K. Nomoto, M. Oguri†(3), A. Oka, N. Okabe, S. Okamoto, S. Okamura, J. Okumura, S. Okumura, Y. Okura, Y. Ono M. Onodera, K. Ota, M. Ouchi†(2), S. Oyabu, P. A. Price, R. Quimby, C. E. Rusu, S. Saito, T. Saito, Y. Saitou, M. Sato, T. Shibuya, K. Shimasaku†(2), A. Shimono, S. Shinogi, M. Shirasaki, J. D. Silverman, D. N. Spergel, M. A. Strauss†(6), H. Sugai, N. Sugiyama, D. Suto, Y. Suto, K. Tadaki, M. Takada†(3), R. Takahashi, S. Takahashi, T. Takata, T. T. Takeuchi, N. Tamura, M. Tanaka, M. Tanaka†(3), M. Tanaka, Y. Taniguchi, A. Taruya, T. Terai, Y. Terashima, N. Tominaga, J. Toshikawa, T. Totani, M. Tsai, E. L. Turner, Y. Ueda, K. Umetsu, Y. Urata†(1), Y. Utsumi, B. Vulcani, K. Wada, S.-Y. Wang, W.-H. Wang, T. Yamada, Y. Yamada, K. Yamamoto, H. Yamanoi, C.-H. Yan, N. Yasuda†(3), A. Yonehara, F. Yoshida†(5), N. Yoshida, M. Yoshikawa, S. Yuma (1) NCU, Taiwan (2) Tokyo (3) Kavli IPMU (4) Tohoku (5) NAOJ

Context. One of the most striking discoveries of the INTEGRAL observatory is the existence of a previously unknown population of X-ray sources in the inner arms of the Galaxy. The investigations of the optical/NIR counterparts of some of them have provided evidence that they are highly absorbed high mass X-ray binaries hosting supergiants. Aims. We aim to identify the optical/NIR counterpart of one of the newly discovered INTEGRAL sources, IGR J16283-4838, and determine the nature of this system. Methods. We present optical and NIR observations of the field of IGR J16283-4838, and use the astrometry and photometry of the sources within it to identify its counterpart. We obtain its NIR spectrum, and its optical/NIR spectral energy distribution by means of broadband photometry. We search for the intrinsic polarization of its light, and its short and long-term photometric variability. Results. We demonstrate that this source is a highly absorbed HMXB located beyond the Galactic center, and that it may be surrounded by a variable circumstellar medium.

Starting jet airflow is investigated in a channel with a pair of consecutive slitted constrictions approximating the true and false vocal folds in the human larynx. The flow is visualized using the Schlieren optical technique and simulated by solving the Navier-Stokes equations for an incompressible two-dimensional viscous flow. Laboratory and numerical experiments show the spontaneous formation of three different classes of vortex dipolar structures in several regions of the laryngeal profile under conditions that may be assimilated to those of voice onset.

In this work, we investigate the nature of the host galaxies of long gamma-ray bursts (LGRBs) using a galaxy catalogue constructed from the Millennium Simulation. We developed an LGRB synthetic model based on the hypothesis that these events originate at the end of the life of massive stars following the collapsar model, with the possibility of including a constraint on the metallicity of the progenitor star. A complete observability pipeline was designed to calculate a probability estimation for a galaxy to be observationally identified as a host for LGRBs detected by present observational facilities. This new tool allows us to build an observable host galaxy catalogue which is required to reproduce the current stellar mass distribution of observed hosts. This observability pipeline predicts that the minimum mass for the progenitor stars should be similar to 75M(circle dot) in order to be able to reproduce BATSE observations. Systems in our observable catalogue are able to reproduce the observed properties of host galaxies, namely stellar masses, colours, luminosity, star formation activity and metallicities as a function of redshift. At z > 2, our model predicts that the observable host galaxies would be very similar to the global galaxy population. We found that similar to 88 per cent of the observable host galaxies with mean gas metallicity lower than 0.6Z(circle dot) have stellar masses in the range 10(8.5)-10(10.3) M-circle dot, in excellent agreement with observations. Interestingly in our model, observable host galaxies remain mainly within this mass range regardless of redshift, since lower stellar mass systems would have a low probability of being observed while more massive ones would be too metal-rich. Observable host galaxies are predicted to preferentially inhabit dark matter haloes in the range 10(11)-10(11.5) M-circle dot, with a weak dependence on redshift. They are also found to preferentially map different density environments at different stages of evolution of the Universe. At high redshifts, the observable host galaxies are predicted to be located in similar environments as the global galaxy population but have a slightly higher probability to have a close companion.