Portraying the missing baryonic mass at the cosmic noon: the contribution of CUBES

The cosmological and galactic missing baryon problems are briefly reviewed with a particular focus on the contributions to the baryonic census derived from the analysis of quasar absorption spectra. The CUBES spectrograph foreseen for the ESO VLT, with its exceptional efficiency (>40\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$> 40$$\end{document} %), blue wavelength coverage (λ≃300-405\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda \simeq 300-405$$\end{document} nm) and intermediate resolution (R≃24,000\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R\simeq 24,000$$\end{document}) will allow us to tackle this issue with two approaches: using H i Lyman-α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} lines at z≃1.5-2.3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z \simeq 1.5-2.3$$\end{document}, just after the peak of star formation, and using O vi absorbers at z≃1.9-2.9\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z \simeq 1.9-2.9$$\end{document}, at the cosmic noon. In both cases, in order to derive the baryonic masses it will be necessary to acquire also higher-resolution spectra of the same target quasars to cover the region at longer wavelengths. We simulate the observations with the CUBES E2E simulator considering a sample of 40 bright quasars at redshifts zem∼2-3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z_\mathrm{em}\sim 2-3$$\end{document} observed for a total time of ∼13\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 13$$\end{document} h to reach a signal-to-noise ratio of ∼15\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 15$$\end{document} in the H i Lyman-α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}, O vi forests.


Introduction
The early universe was a highly uniform and homogeneous mix of dark and baryonic matter. The ratio between these two components as predicted by Big Bang Nucleosynthesis [10] and inferred from density fluctuations of the cosmic microwave background [19,41] 16 . This primordial mix is expected to persist as large scale structure evolves and individual gravitationally bound objects emerge during the course of cosmic time. In this hypothesis, the baryonic mass of any given cosmological structure would be On the other hand, direct measurements reveal that at z ∼ 0 a baryonic fraction of ∼ 30 − 40 % is missing on cosmological scales (e.g. [11,50]) and also galaxies are baryon deficient. The amount of missing baryons in galaxies depends on the mass/ scale of the system: while clusters are in agreement with the primordial fraction, in galaxies outside of clusters we are missing more than 60% of the baryons, rising to ≥ 90 % for gas-dominated galaxies and dwarfs 1 (e.g. [28], see Fig. 1).
Where are the missing baryons? Cosmological, hydrodynamical simulations (e.g. [7,8,14,27,53,57]) predict that, at z ∼ 0 , most of the baryons are in the form of unvirialized, diffuse gas: 40 − 50 % have been shock heated by the processes of structure formation and resides in the warm-hot intergalactic medium (WHIM) Fig. 1 Relation between baryonic mass and rotation velocity for galaxies. Red squares represent Local Group dwarf satellites. Round symbols represent rotationally supported disks. Purple squares represent the mean of many galaxy clusters. If structures possessed baryons in the cosmic fraction, they would fall on the solid line. The dashed line represents the Baryonic Tully-Fisher relation defined by the disks. Reproduced from [28] by permission of the AAS characterized by temperatures of 10 5 K < T < 10 7 K . Another 30 − 40 % is in the diffuse warm medium at T < 10 5 K. The cold gas and stars in galaxies contribute less than ∼ 20 % and the very hot gas with T > 10 7 K (the normal X-ray emitting gas, predominantly in collapsed and virialized clusters of galaxies) gives a negligible contribution. As we move back in time, the warm diffuse component becomes increasingly important and, at z > ∼ 2, it is expected to be the dominant reservoir of baryons for galaxy formation: the intergalactic medium (IGM). The WHIM is predicted to reside mainly in the filaments of the cosmic web where also galaxies lie, connecting the two scales of the missing baryon problem (Fig. 2).
The aim of this paper is: to report on the state-of-the-art of the observations that at z < 1 have investigated the contribution to the baryon fraction of the different gaseous components (Sect. 2) and, to describe what is known and what is foreseen for the future on the quest for baryons at higher redshift (Sect. 3), with a particular focus on the measurements that will be enabled by the CUBES spectrograph (Sect. 3.2).

Measurements in the FUV
In order to estimate the contribution to the local baryon budget of the gas outside galaxies, the main observational technique has been far ultra-violet (FUV) spectroscopy of bright background sources.
The H i Lyman-absorption lines, forming the so called Lymanforest observed in quasar spectra, are known to probe the neutral component of the diffuse, ionized IGM at temperature T ∼ 10 3.5 − 10 4.5 K [29]. [50] carry out a detailed study of a sample of 746 H i absorbers observed with STIS/HST from the catalogue compiled by [54]; they apply ionization corrections determined both analytically and from numerical simulations and derive a contribution to the baryon budget Ω (HI) b = 28% ± 11% 2 in agreement with simulation predictions and previous measurements (e.g. [23]).
The FUV lines due to the O vi doublet ( 1031.9, 1037.6 Å) are currently the best available tracers of the ''cooler'' portion of the shock-heated gas forming the WHIM, reaching a peak ionization fraction, f OVI ≈ 0.22 , at temperature T max = 10 5.45 K in collisional ionization equilibrium [52] and similar values in non-equilibrium cooling [18]. Surveys based on FUSE, STIS/HST and COS/HST data assuming fixed metallicity, Z∕Z ⊙ = 0.1 , and ionization correction f OVI = 0.2 , obtained a contribution to the baryon fraction of ∼ 8 − 9 % (e.g. [13,54,55]). Using simulations to determine a more reliable value of the product of metallicity and ionization fraction, [50] derive a column-density-weighted mean of f OVI (Z∕Z ⊙ ) = 0.01 , a factor of two smaller than previously assumed. This implies, Ω (OVI) b ≃ 17 ± 4 % summing up to the WHIM census in the temperature range T ∼ 10 5.3 − 10 5.7 K.
There are however two main uncertainties affecting this measurement: the value of the gas metallicity that could be larger (e.g. [24,55]), implying a lower baryon content and, in particular, the fact that O vi absorbers are subject to different physical conditions tracing both collisionally and photo-ionized gas [55]. Since the latter is also traced by the Lyman-forest, considering all O vi absorbers could amount to a certain degree of double counting.
Unfortunately, a precise determination of the metallicity of O vi absorbers is generally hampered by the scarcity of other associated metal lines. Furthermore, if the missing baryons are predominantly located in gas that is shock heated when it first accreted into galaxy potential wells, then the WHIM plasma could have a relatively low metallicity, and O vi could be difficult to detect in currently available data. An alternative route to detecting WHIM absorbers, and one that does not rely upon metal enrichment, is to detect and measure broad Lyman-absorbers (BLAs) in high-signal-to-noise ratio spectra. Absorption line width is usually denoted in terms of the corresponding Doppler parameter b ≡ FWHM∕2 √ ln2 ≡ √ 2 where is the Gaussian width. The thermal width, which contributes to the total line width, is simply a function of gas temperature T and atomic mass A: where T 5 is temperature in units of 10 5 K. BLAs are defined as Lyman-absorption lines with b ≥ 40 km s −1 and depending on the total gas column density of the WHIM absorber and its temperature, they could have column densities 12.5 ≤ log N (H i) ≤ 14 . The quantitative estimate of the baryon content of the BLAs is complicated by the effect of non-thermal line broadening processes, line blends, and noise features that can mimic broad spectral troughs. Observational studies report relatively reliable lower limits of 6 − 10 % increasing to ≈ 14 − 20 % when corrected with simulations [12,23,46,47].

Measurements in the soft X-rays
Observations in the FUV probe WHIM temperatures up to maximum T ∼ 10 5.7 K. The hot phase of the WHIM can be investigated through absorptions of ionized metals in the soft X-rays [3]. In this energetic regime the most promising ionic transitions are O vii (21.60 Å) and O viii (18.97 Å) since they are strong transitions of an abundant element with peak collisional ionization equilibrium (CIE) abundances at temperatures of (12) × 10 6 K. Other metal transitions that could be detected in absorption in the X-rays (note that this list is by no means complete) are N vii 24.78 Å, Ne ix 13.45 Å, and the Lyman-line of O vii at 18.63 Å.
The mostly used tracers of this WHIM phase are the O vii He absorption lines; these are predicted to be relatively narrow (with a Doppler parameter thermal component b th ≈ 20 − 46 km s −1 ), extremely shallow (rest-frame equivalent widths, EW ≤ 10 mÅ) and rare [8,14]. Current X-ray spectrometers, XMM-Newton and Chandra, do not resolve such lines and thus need a signal-to-noise ratio per resolution element greater than or equal to ∼ 20 in the continuum to detect them at a significance level exceeding ∼ 3 . Although substantial observational efforts have been carried out in the last 10 years to target the rare brightest sources at z ≥ 0.3 , no firm, high-significance detections have been reported up to now (see [36] for a recent review). Even the most recent claim of a double detection of O vii at z = 0.3551 and z = 0.4339 [35] has been subject to revisions that decreased its reliability [21,34]. The full exploitation of the soft X-ray metal absorption lines to investigate the nature of the WHIM will have to wait for the next generation of instruments, in particular the X-ray integral field unit spectrometer (XIFU: [1]) on board the Athena observatory 3 , which is due for launch in the early 2030s.

Accounting for baryons in galaxy halos
The baryon mass estimates described in the previous sections are based on surveys of absorption lines tracing different environments: from the diffuse IGM, to filaments in the cosmic web and the outskirts of galaxies.
Combining absorption line studies with the identification of galaxies in the field of the background source allows us to investigate the galactic halo properties. The gaseous halo extending out to the virial radius and beyond -dubbed circumgalactic medium (CGM) -is where a large fraction of the baryons missing from galaxies is expected to reside (e.g. [7]). (Fig. 2).
Several studies have shown that the CGM is a major reservoir of heavy elements with a mass rivaling and possibly exceeding that within galaxies [39,56,59]. Metals have to be transported from the ISM to ∼ 100 kpc distances in the CGM and to attain this goal several mechanisms have been proposed including galactic winds, active galactic nucleus (AGN) feedback, accretion, tidal stripping, and ram pressure. The detailed study of the metallicity and ionization state of the CGM is of great importance to establish which one of those mechanisms was more relevant for the enrichment. Furthermore, heavy element ions are critical to By visual inspection, we see that the diffuse IGM is ubiquitous, whereas the WHIM is concentrated near knots and filaments. Adapted from [27] establish which ionization corrections have to be applied to the measured neutral hydrogen column density (N(H i)) to infer the total gaseous halo mass.
An estimate of the baryon content of the CGM at z ∼ 0 has been carried out in particular in the context of the COS-Halos survey [60], which investigated the halo properties of a sample of L ≈ L * galaxies whose angular offsets from quasar sightlines and photometric redshifts implied impact parameters ( R < 160 kpc) well inside their virial radii.
Thanks to the detection in COS spectra of low and intermediate-ionization absorption lines of heavy elements, together with the corresponding H iabsorption lines, it was possible to assess the physical conditions and mass of highly ionized, cool ( T ≈ 10 4 K) CGM gas observed within 160 kpc of 33 low-redshift, L ≈ L * galaxies [59]. The derived contribution of the cool CGM to the baryonic budget ( 0.16 M halo ) ranges between 25% and 45% (see Fig. 3). Note that this determination is subject to significant uncertainties [44]. On the one hand, several systems had an uncertain determination of the H icolumn density due to the lack of spectral coverage of the whole Lyman series: the precise measurement of N(H i) is critical for the determination of the halo mass. On the other hand, the very simplified standard photoionization models applied to determine the ionization Fig. 3 Baryonic budget of the fiducial COS-Halos galaxy, at L ≈ L * , represented as a bar chart showing the most massive baryonic components of the galaxy. The solid filled bars are lower limits to the fraction each component will contribute, while the hatched area above the solid bars shows potential additional contributions allowed by the data. The sum of these components is given by the brown bar, illustrating that galaxies have anywhere between 45% and 100% of their baryons relative to the cosmological fraction. Adapted from [59] and reproduced by permission of the AAS correction are known to fail when a wider set of metal ions is considered, suggesting a complex density structure and/or additional ionization mechanisms. More sophisticated models should be developed to reproduce the ionization structure observed in the CGM.
As described in Sect. 2.1, the O vi absorption doublet is the most commonly used tracer of the WHIM gas, being one of the few ions easily observed at rest wavelength r > 912 Å with ionization potential (IP) > ∼ 100 eV. However, discriminating between photoionized and collisionally ionized O vi has proved to be difficult (e.g. [58]). The COS Absorption Survey of Baryon Harbors (CASBaH) obtained high spectral resolution COS/HST and STIS/HST spectra of nine FUV-bright quasars at 0.92 < z em < 1.48 without known Lyman Limit systems, in order to access transitions at r < 912 Å. The quasar spectra are complemented by a survey for galaxies in the fields surrounding the nine sightlines, for a total of ∼ 5900 galaxies whithin ≈ 10 comoving Mpc [42]. A goal of the CASBaH survey was to investigate the presence of the Ne viii 770.41, 780.32 Å doublet in the CGM of galaxies close to the line of sight. With IP = 207 eV, Ne viii can persist in gas with 10 5 < T < 10 6 K, and thereby complements O vi with the potential to break degeneracies in ionization models. [4] detected Ne viii associated with 9 galaxies of the 29 that they could consider for their experiment. They estimated the mass in Ne viii-traced gas to be M gas (Ne viii ) ≥ 10 9.5 M ⊙ (Z∕Z ⊙ ) −1 , or 6%-20% of the expected baryonic mass depending on the Ne viii absorbers' metallicity. This estimate assumes the peak of the ionization fraction, Ne viii/Ne, under collisional ionization, or T ∼ 10 5.8 K. Interestingly, this is similar to the virial temperature for the median halo mass and redshift ( ⟨z⟩ = 0.68 ) of their sample, implying that the Ne viii-bearing gas is plausibly collisionally ionized near this temperature.

Is the missing baryon problem solved?
Recently, a new astronomical phenomenon has been discovered: the fast radio bursts (FRB), radio flashes of unknown physical nature with duration of milliseconds. The dispersion measure (DM) of a short and coherent burst of radiation from a distant source offers an integral constraint on the electron distribution along the path between the source and the Earth. This includes the IGM, our Galaxy, our Local Group, the galaxy hosting the FRB, and the baryons residing in other galactic haloes near the sightline (e.g. [30,43,49]).
Thanks to the repeatability of some of these FRB it has been possible to localize them and associate them with galactic counterparts (e.g. [9]). [25] have recently exploited the dispersion observed in a sample of localized FRB at z ∼ 0.1 − 0.5 to measure the electron column density and account for every ionised baryon along the line of sight. This independent measurement of the baryon content of the Universe with a golden sample of five FRBs is consistent with the CMB and Big Bang Nucleosynthesis values, substantially solving the missing baryon problem at low redshift.
This result has established the total amount of baryons, confirming that they are distributed in several environments. However, with the present sample of localized FRB it is not possible to distinguish the share of the different gaseous components. In this respect, quasar absorption lines studies will continue to play a key role, in particular, they could help understand whether a substantial fraction of the mass at high redshift is also contained in diffuse halo components.

Baryons in the high redshift IGM/CGM
At z ∼ 3 most of the baryons ( ∼ 90 %) are in the diffuse IGM which is traced by the Lyman-forest [16,45]. Also the diffuse gas in the halo of galaxies is expected to contribute significantly to the total budget, hosting the (metal enriched) gas expelled from galaxies by stellar and AGN feedback processes and also the gas accreting from the IGM to fuel the star formation (e.g. [38]). The contribution of the warm-hot gas to the baryon census at z ∼ 2.5 was determined in [51] by studying the O vi absorbers along 5 lines of sight observed with HIRES at Keck. They found that these systems arise at the interface between the IGM and galaxies, in the CGM, and that they contribute ≥ 5 % to Ω b . A very similar experiment was carried out by [6] with 2 lines of sight observed with UVES/VLT, probing O vi absorbers at z ∼ 2 . These authors found that gas traced by O vi represent ∼ 20 % of Ω b , implying a large uncertainty in this determination. [22] investigated the baryonic content of the CGM as traced by absorbers optically thick at the Lyman limit (characterized by column densities N(H i ) ≥ 10 17 cm −2 ) at z ∼ 2 − 3 . They determined the metallicity of the absorbers, which is critical to estimate the ionized fraction of the gas (since the ionized gas is only observed via the metal lines) and obtained that highly ionized gas in the LL > 1 absorbers is very likely the second largest contributor after the Lyman-forest to the baryon budget at z ∼ 2 − 3 and also a substantial contributor to the metal budget. This work was based on 22 quasar spectra obtained with the HIRES/Keck spectrograph and selected from the 1st data release of the KODIAQ survey [37]. The considered spectra are characterized by slightly different resolutions and signal-to-noise ratios in the O vi forest varying from 6 to 54. To our knowledge there are no other recent studies that dealt with the determination of the baryonic content of the z > ∼ 2 CGM and we think that further investigations are needed based on homogeneous and unbiased samples of quasar spectra. We note that a good coverage and, in particular, a large efficiency in the UV region of the spectrum is essential to cover more transitions in the Lyman series and determine with greater accuracy the H icolumn density, from which the halo mass is derived.

Future measurements
Several current and future instruments are being designed with the study of the highredshift (z > ∼ 2) CGM in mind. One of the cases developed for the ELT-MOSAIC instrument ( [48], see Fig. 4 and note that the goal is to extend the visible coverage of MOSAIC down to 390 nm) is to use high-redshift Lyman-break galaxies as background sources to reconstruct the 3D density field of the IGM and probe the CGM of galaxies at z ∼ 3 [20]. The Keck Cosmic Web Imager (KCWI) instrument [32], commissioned in 2018, has partly been designed to address this topic and will certainly make significant contributions, but its blueward coverage stops at 350 nm. Observations in the ground-UV can play an exciting role in this field, considering also the fact that, with the present design, the ELT will be almost blind at < ∼ 400 nm. The CUBES spectrograph, foreseen for the ESO VLT, has been designed to be extremely efficient ( > 40 %) in the wavelength range ∼ 300 − 405 nm with a resolving power R ≃ 24, 000 [61] and with the possibility to have a fiberlink connection with the UVES/ VLT spectrograph to cover simultaneously the region at > 420 nm. These characteristics, makes CUBES an ideal instrument to carry out further investigations of the cosmological baryonic mass content and of the contribution to the baryonic budget of the CGM in redshift regimes not particularly explored before, as we describe in the next section.

The CUBES baryon experiment
The spectral properties of CUBES allow us to foresee two experiments for the measure of the baryonic mass with quasar absorption lines. In the first one, the baryons in the diffuse IGM will be chased through the detection and analysis of H iLyman-absorption lines at z ≃ 1.5 to 2.3 (see Fig. 4). This redshift range, immediately after the era of peak star-formation in the Universe (e.g. [26]), is poorly investigated due to observational difficulties (namely the low efficiency of ground-based spectrographs in the UV) but critical to connect the low and high redshift results. The relatively low number density of lines in the Lyman-forest at these redshifts implies that, even with observations at R ∼ 24, 000 , the impact of blending with metal lines in the Lyman-forest can be minimised (which is not the case at larger redshifts). To complement the information on the neutral IGM component derived from the Fig. 4 Lines available in the ground-UV for studies of the CGM of galaxies at 1.5 < z < 2.3 (or 1.9 < z < 2.9 for O vi absorbers). They span a range of ionisation states (green = low, red = high) to study the relative fractions of cold vs. warm vs. hot gas, at redshifts where the contamination of key ions such as O vi by the Lyman-forest is less severe. CUBES observations will neatly complement future observations in the blue-visible with ELT-MOSAIC, which is targeting the CGM in galaxy halos at z ∼ 3 Lyman-forest, we will need to check also the associated metal absorption lines (in particular due to C iv and Si iv) and derive the contribution of the ionised gas. To this purpose, UVES/VLT (or ESPRESSO/VLT) spectra of the same targets will have to be obtained (or retrieved from the ESO archive) to cover the wavelength range at > 405 nm at higher resolution.
In the second experiment, we will use the important O vi absorption line at 1.9 < z < 2.9 to trace the warm-hot gas at T > 10 5 K, associated with the IGM (e.g. [11]) or with the CGM, using the approach described in [22]. Also in this case, CUBES observations will need complementary UVES (or ESPRESSO) spectra to complete the coverage of the associated H i and metal transitions.
Based on previous works, we estimated the time necessary to carry out the CUBES baryon experiments with a sample of ∼ 40 bright quasars at redshifts z em ∼ 2 − 3.
For our purpose, we want to be sensitive at 3 to H i Lyman-lines with column density N(H i) We used the CUBES E2E simulator [17] with the custom composite quasar spectrum QSO_composite_absorbed.dat with a superposed Lyman-forest, and simulated observations carried out three days from new moon (Fig. 5). Considering an ideal sample formed by 10 targets at z em = 2 , 20 at z em = 2.5 and 10 at z em = 3 , all with a conservative magnitude u(AB) = 18 , we derived that observing times of 900s ×2 , 600s ×2 and 300s ×3 are necessary at the three considered redshifts to reach the desired signal-to-noise ratio in the Lyman-and O vi forests, summing up to a total exposure time of ∼ 13 hours.
We note that in the SQUAD sample [33], compiled from UVES archival quasar spectra, there are ∼ 40 targets with 2 ≤ z em ≤ 3 and signal-to-noise ratio SNR ≥ 20 at 350nm (13 with SNR ≥ 40 ), covering the wavelength range 305-400nm, although not always completely. We plan to use these spectra (and others with similar wavelength coverage and SNR from KODIAQ, if available) for a preparatory work that will allow us to determine the optimal number and redshift distribution of targets for the CUBES experiments.
The ongoing QUBRICS survey [2,5] which identifies quasars in the southern hemisphere with i(AB) < 18 mag, will provide the bright targets needed for this case.

Conclusions
It is almost thirty years that scientists look for a definitive answer to the question "Where are the baryons?" (e.g. [16,40]). Today, we can probably state that all the baryons predicted by the CMB and Big-Bang nucleosynthesis measurements have been detected: at low redshifts thanks to the dispersion measure in FRBs [25], at z > ∼ 2 based on the Lyman-forest observations and simulations (e.g. [31]). This settled, the next step is to investigate in details how baryonic matter is distributed among the different gaseous components with the aim of understanding and constraining the mechanisms that determined the observed distribution (stellar and AGN feedback, accretion, etc.).
The CUBES spectrograph, with its exceptional efficiency in the UV and relatively high resolution, will play a critical role on these issues, in particular it will allow us Fig. 5 Output of the E2E simulator for the case of a quasar at z em = 2.5 with magnitude u(AB) = 18 observed for 600s. The upper panel shows the input spectrum (solid black line, target raw), the input spectrum with atmospheric extinction (dashed black line, target extincted), the spectrum that is seen by the detector, i.e. scaled by telescope and instrument efficiency (solid black line, on detector), the background seen on the detector (solid blue line, on detector background), the spectrum recovered by simple (non-optimal) extraction (red dots, extracted), the error on the extracted spectrum (green dots, extracted error). In the lower panel the signal-to-noise ratio per pixel is reported to dig into the complex nature of the inter and circum-galactic gas at z ∼ 1.5 − 3 . For reference, we determined that an exposure time of ∼ 13 hours would be needed with CUBES to observe a sample of ∼ 40 bright background quasars at z em ∼ 2 − 3 , reaching the necessary signal-to-noise ratio for the proposed science case. Its improved performances with respect to UVES/VLT will open up background sources 1-2 magnitudes fainter than the quasars used at present, significantly increasing the number (and the space density) of available targets and/or the possibility to obtain very high signal-to-noise ratio spectra to detect the faint lines (e.g. [15]).