Showing 1–50 of 286 results for author: Cristiani, S

Three super-Earths and a possible water world from TESS and ESPRESSO

Authors: M. J. Hobson, F. Bouchy, B. Lavie, C. Lovis, V. Adibekyan, C. Allende Prieto, Y. Alibert, S. C. C. Barros, A. Castro-González, S. Cristiani, V. D'Odorico, M. Damasso, P. Di Marcantonio, X. Dumusque, D. Ehrenreich, P. Figueira, R. Génova Santos, J. I. González Hernández, J. Lillo-Box, G. Lo Curto, C. J. A. P. Martins, A. Mehner, G. Micela, P. Molaro, N. J. Nunes , et al. (29 additional authors not shown)

Abstract: Since 2018, the ESPRESSO spectrograph at the VLT has been hunting for planets in the Southern skies via the RV method. One of its goals is to follow up candidate planets from transit surveys such as the TESS mission, particularly small planets. We analyzed photometry from TESS and ground-based facilities, high-resolution imaging, and RVs from ESPRESSO, HARPS, and HIRES, to confirm and characterize… ▽ More Since 2018, the ESPRESSO spectrograph at the VLT has been hunting for planets in the Southern skies via the RV method. One of its goals is to follow up candidate planets from transit surveys such as the TESS mission, particularly small planets. We analyzed photometry from TESS and ground-based facilities, high-resolution imaging, and RVs from ESPRESSO, HARPS, and HIRES, to confirm and characterize three new planets: TOI-260 b, transiting a late K-dwarf, and TOI-286 b and c, orbiting an early K-dwarf. We also update parameters for the known super-Earth TOI-134 b , hosted by an M-dwarf. TOI-260 b has a $13.475853^{+0.000013}_{-0.000011}$ d period, $4.23 \pm1.60 \mathrm{M_\oplus}$ mass and $1.71\pm0.08\mathrm{R_\oplus}$ radius. For TOI-286 b we find a $4.5117244^{+0.0000031}_{-0.0000027}$ d period, $4.53\pm0.78\mathrm{M_\oplus}$ mass and $1.42\pm0.10\mathrm{R_\oplus}$ radius; for TOI-286 c, a $39.361826^{+0.000070}_{-0.000081}$ d period, $3.72\pm2.22\mathrm{M_\oplus}$ mass and $1.88\pm 0.12\mathrm{R_\oplus}$ radius. For TOI-134 b we obtain a $1.40152604^{+0.00000074}_{-0.00000082}$ d period, $4.07\pm0.45\mathrm{M_\oplus}$ mass, and $1.63\pm0.14\mathrm{R_\oplus}$ radius. Circular models are preferred for all, although for TOI-260 b the eccentricity is not well-constrained. We compute bulk densities and place the planets in the context of composition models. TOI-260 b lies within the radius valley, and is most likely a rocky planet. However, the uncertainty on the eccentricity and thus on the mass renders its composition hard to determine. TOI-286 b and c span the radius valley, with TOI-286 b lying below it and having a likely rocky composition, while TOI-286 c is within the valley, close to the upper border, and probably has a significant water fraction. With our updated parameters for TOI-134 b, we obtain a lower density than previous findings, giving a rocky or Earth-like composition. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 61 pages (of which pp. 24-61 are appendices), 20 figures (main text). Accepted for publication in A&A

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

GHOST Commissioning Science Results III: Characterizing an iron-poor damped Lyman $α$ system

Authors: Trystyn A. M. Berg, Christian R. Hayes, Stefano Cristiani, Alan McConnachie, J. Gordon Robertson, Federico Sestito, Chris Simpson, Fletcher Waller, Timothy Chin, Adam Densmore, Ruben J. Diaz, Michael L. Edgar, Javier Fuentes Lettura, Manuel Gómez-Jiménez, Venu M. Kalari, Jon Lawrence, Steven Margheim, John Pazder, Roque Ruiz-Carmona, Ricardo Salinas, Karleyne M. G. Silva, Katherine Silversides, Kim A. Venn

Abstract: The Gemini High-resolution Optical SpecTrograph (GHOST) is a new echelle spectrograph available on the Gemini-South telescope as of Semester 2024A. We present the first high resolution spectrum of the quasar J1449-1227 (redshift z_em=3.27) using data taken during the commissioning of GHOST. The observed quasar hosts an intervening iron-poor ([Fe/H] = -2.5) damped Lyman alpha (DLA) system at redshi… ▽ More The Gemini High-resolution Optical SpecTrograph (GHOST) is a new echelle spectrograph available on the Gemini-South telescope as of Semester 2024A. We present the first high resolution spectrum of the quasar J1449-1227 (redshift z_em=3.27) using data taken during the commissioning of GHOST. The observed quasar hosts an intervening iron-poor ([Fe/H] = -2.5) damped Lyman alpha (DLA) system at redshift z=2.904. Taking advantage of the high spectral resolving power of GHOST (R~55000), we are able to accurately model the metal absorption lines of the metal-poor DLA and find a supersolar [Si/Fe], suggesting the DLA gas is in an early stage of chemical enrichment. Using simple ionization models, we find that the large range in the C IV/Si IV column density ratio of individual components within the DLA's high ionization absorption profile can be reproduced by several metal-poor Lyman limit systems surrounding the low-ionization gas of the DLA. It is possible that this metal-poor DLA resides within a complex system of metal-poor galaxies or filaments with inflowing gas. The high spectral resolution, wavelength coverage and sensitivity of GHOST makes it an ideal spectrograph for characterizing the chemistry and kinematics of quasar absorption lines. △ Less

Submitted 18 April, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS. 8 Pages, 5 figures

Parameter estimation from Ly$α$ forest in Fourier space using Information Maximising Neural Network

Authors: Soumak Maitra, Stefano Cristiani, Matteo Viel, Roberto Trotta, Guido Cupani

Abstract: We aim to present a robust parameter estimation with simulated Lya forest spectra from Sherwood-Relics simulations suite using Information Maximizing Neural Network(IMNN) to extract maximal information from Lya 1D-transmitted flux in Fourier space. We perform 1D estimations using IMNN for IGM thermal parameters $T_0$ & $γ$ at z=2-4 and cosmological parameters $σ_8$ & $n_s$ at z=3-4. We compare our… ▽ More We aim to present a robust parameter estimation with simulated Lya forest spectra from Sherwood-Relics simulations suite using Information Maximizing Neural Network(IMNN) to extract maximal information from Lya 1D-transmitted flux in Fourier space. We perform 1D estimations using IMNN for IGM thermal parameters $T_0$ & $γ$ at z=2-4 and cosmological parameters $σ_8$ & $n_s$ at z=3-4. We compare our results with estimates from power spectrum using posterior distribution from Markov Chain Monte Carlo(MCMC). We then check robustness of IMNN estimates against deviation in spectral noise levels,continuum uncertainties & instrumental smoothing effects. Using mock Lya forest sightlines from publicly available CAMELS project we also check the robustness of the trained IMNN on a different simulation. We also perform a 2D-parameter estimation for $T_0$ & HI photoionization rates $Γ_{HI}$. We obtain improved estimates of $T_0$ & $γ$ using IMNN over standard MCMC approach. These estimates are also more robust against SNR deviations at z=2 & 3. At z=4 the sensitivity to noise deviations is on par with MCMC estimates. The IMNN also provides $T_0$ and $γ$ estimates which are robust against continuum uncertainties by extracting continuum-independent small-scale information from Fourier domain. In case of $σ_8$ & $n_s$ IMNN performs on par with MCMC but still offers a significant speed boost in estimating parameters from a new dataset. The improved estimates with IMNN are seen for high instrumental-resolution(FWHM=6km/s). At medium or low resolutions IMNN performs similar to MCMC suggesting an improved extraction of small-scale information with IMNN. We also find that IMNN estimates are robust against the choice of simulation. By performing a 2D-parameter estimation for $T_0$ & $Γ_{HI}$ we also demonstrate how to take forward this approach observationally in the future. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Comments: 17 pages, 7 figures, 5 tables. Submitted to Astronomy&Astrophysics. Comments welcomed

Probing the small scale structure of the Inter-Galactic Medium with ESPRESSO: spectroscopy of the lensed QSO UM673

Authors: Stefano Cristiani, Guido Cupani, Andrea Trost, Valentina D'Odorico, Francesco Guarneri, Gaspare Lo Curto, Massimo Meneghetti, Paolo Di Marcantonio, João P. Faria, Jonay I. González Hernández, Christophe Lovis, Carlos J. A. P. Martins, Dinko Milaković, Paolo Molaro, Michael T. Murphy, Nelson J. Nunes, Francesco Pepe, Rafael Rebolo, Nuno C. Santos, Tobias M. Schmidt, Sérgio G. Sousa, Alessandro Sozzetti, María Rosa Zapatero Osorio

Abstract: The gravitationally lensed quasar J014516.6-094517 at z=2.719 has been observed with the ESPRESSO instrument at the ESO VLT to obtain high-fidelity spectra of the two images A and B with a resolving power R=70000. At the redshifts under investigation (2.1 < z < 2.7), the Lyman forests along the two sightlines are separated by sub-kiloparsec physical distances and exhibit a strong correlation. We f… ▽ More The gravitationally lensed quasar J014516.6-094517 at z=2.719 has been observed with the ESPRESSO instrument at the ESO VLT to obtain high-fidelity spectra of the two images A and B with a resolving power R=70000. At the redshifts under investigation (2.1 < z < 2.7), the Lyman forests along the two sightlines are separated by sub-kiloparsec physical distances and exhibit a strong correlation. We find that the two forests are indistinguishable at the present level of signal-to-noise ratio and do not show any global velocity shift, with the cross-correlation peaking at $Δv = 12 \pm 48$ m/s. The distribution of the difference in velocity of individual Lyman-$α$ features is compatible with a null average and a mean absolute deviation of 930 m/s. Significant differences in NHI column density are not detected, putting a limit to the RMS fluctuation in the baryon density on $\leq 1$ proper kpc scales of $Δρ/ ρ< 3$%. On the other hand, metal lines show significant differences both in velocity structure and in column density. A toy model shows that the difference in velocity of the metal features between the two sightlines is compatible with the the motions of the baryonic component associated to dark matter halos of typical mass $M\simeq 2\times 10^{10} M_\odot$, also compatible with the observed incidence of the metal systems. The present observations confirm the feasibility of the Sandage test of the cosmic redshift drift with high-fidelity spectroscopy of the Lyman forest of distant, bright quasars, but also provide an element of caution about the intrinsic noise associated to the usage of metal features for the same purpose. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 17 pages, 10 figures, 5 tables, accepted for publication in MNRAS

Fundamental physics with ESPRESSO: a new determination of the D/H ratio towards PKS1937-101

Authors: Francesco Guarneri, Luca Pasquini, Valentina D'Odorico, Stefano Cristiani, Guido Cupani, Paolo Di Marcantonio, J. I. González Hernández, C. J. A. P. Martins, Alejandro Suárez Mascareño, Dinko Milaković, Paolo Molaro, Michael T. Murphy, Nelson J. Nunes, Enric Palle, Francesco Pepe, Rafael Rebolo, Nuno C. Santos, Ricardo Génova Santos, Tobias M. Schmidt, Sérgio G. Sousa, Alessandro Sozzetti, Andrea Trost

Abstract: Primordial abundances of light elements are sensitive to the physics of the early Universe and can directly constrain cosmological quantities, such as the baryon-to-photon ratio $η_{10}$, the baryon density and the number of neutrino families. Deuterium is especially suited for these studies: its primordial abundance is sensitive and monotonically dependent on $η_{10}$, allowing an independent mea… ▽ More Primordial abundances of light elements are sensitive to the physics of the early Universe and can directly constrain cosmological quantities, such as the baryon-to-photon ratio $η_{10}$, the baryon density and the number of neutrino families. Deuterium is especially suited for these studies: its primordial abundance is sensitive and monotonically dependent on $η_{10}$, allowing an independent measurement of the cosmic baryon density that can be compared, for instance, against the Planck satellite data. The primordial deuterium abundance can be measured in high $H_I$ column density absorption systems towards distant quasars. We report here a new measurement, based on high-resolution ESPRESSO data, of the primordial $D_I$ abundance of a system at redshift $z \sim 3.572$, towards PKS1937-101. Using only ESPRESSO data, we find a D/H ratio of $2.638\pm0.128 \times 10^{-5}$, while including the available UVES data improves the precision, leading to a ratio of $2.608 \pm 0.102 \times 10^{-5}$. The results of this analysis agree with those of the most precise existing measurements. We find that the relatively low column density of this system ($\log{N_{\rm H_I}/ {\rm cm}^{-2}}\sim18 $) introduces modelling uncertainties, which become the main contributor to the error budget. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: Accepted for publication in MNRAS. 16 pages, 10 figures and 6 tables (4 in main text, 2 in Appendix). Full model is available at the end of the manuscript

TESS and ESPRESSO discover a super-Earth and a mini-Neptune orbiting the K-dwarf TOI-238

Authors: A. Suárez Mascareño, V. M. Passegger, J. I. González Hernández, D. J. Armstrong, L. D. Nielsen, C. Lovis, B. Lavie, S. G. Sousa, A. M. Silva, R. Allart, R. Rebolo, F. Pepe, N. C. Santos, S. Cristiani, A. Sozzetti, M. R. Zapatero Osorio, H. M. Tabernero, X. Dumusque, S. Udry, V. Adibekyan, C. Allende Prieto, Y. Alibert, S. C. C. Barros, F. Bouchy, A. Castro-González , et al. (31 additional authors not shown)

Abstract: The number of super-Earth and mini-Neptune planet discoveries has increased significantly in the last two decades thanks to transit and radial velocity surveys. When it is possible to apply both techniques, we can characterise the internal composition of exoplanets, which in turn provides unique insights on their architecture, formation and evolution. We performed a combined photometric and radi… ▽ More The number of super-Earth and mini-Neptune planet discoveries has increased significantly in the last two decades thanks to transit and radial velocity surveys. When it is possible to apply both techniques, we can characterise the internal composition of exoplanets, which in turn provides unique insights on their architecture, formation and evolution. We performed a combined photometric and radial velocity analysis of TOI-238 (TYC 6398-132-1), which has one short-orbit super-Earth planet candidate announced by NASA's TESS team. We aim to confirm its planetary nature using radial velocities taken with the ESPRESSO and HARPS spectrographs, to measure its mass and to detect the presence of other possible planetary companions. We carried out a joint analysis by including Gaussian processes and Keplerian orbits to account for the stellar activity and planetary signals simultaneously. We detected the signal induced by TOI-238 b in the radial velocity time-series, and the presence of a second transiting planet, TOI-238 c, whose signal appears in RV and TESS data. TOI-238 b is a planet with a radius of 1.402$^{+0.084}_{-0.086}$ R$_{\oplus}$ and a mass of 3.40$^{+0.46}_{-0.45}$ M$_{\oplus}$. It orbits at a separation of 0.02118 $\pm$ 0.00038 AU of its host star, with an orbital period of 1.2730988 $\pm$ 0.0000029 days, and has an equilibrium temperature of 1311 $\pm$ 28 K. TOI-238 c has a radius of 2.18$\pm$ 0.18 R$_{\oplus}$ and a mass of 6.7 $\pm$ 1.1 M$_{\oplus}$. It orbits at a separation of 0.0749 $\pm$ 0.0013 AU of its host star, with an orbital period of 8.465652 $\pm$ 0.000031 days, and has an equilibrium temperature of 696 $\pm$ 15 K. The mass and radius of planet b are fully consistent with an Earth-like composition, making it likely a rocky super-Earth. Planet c could be a water-rich planet or a rocky planet with a small H-He atmosphere. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: 33 pages, 31 figures, 5 tables. Accepted for publication at A&A

The compact multi-planet system GJ 9827 revisited with ESPRESSO

Authors: V. M. Passegger, A. Suárez Mascareño, R. Allart, J. I. González Hernández, C. Lovis, B. Lavie, A. M. Silva, H. M. Müller, H. M. Tabernero, S. Cristiani, F. Pepe, R. Rebolo, N. C. Santos, V. Adibekyan, Y. Alibert, C. Allende Prieto, S. C. C. Barros, F. Bouchy, A. Castro-González, V. D'Odorico, X. Dumusque, P. Di Marcantonio, D. Ehrenreich, P. Figueira, R. Génova Santos , et al. (14 additional authors not shown)

Abstract: GJ 9827 is a bright, nearby K7V star orbited by two super-Earths and one mini-Neptune on close-in orbits. The system was first discovered using K2 data and then further characterized by other spectroscopic and photometric instruments. Previous literature studies provide several mass measurements for the three planets, however, with large variations and uncertainties. To better constrain the planet… ▽ More GJ 9827 is a bright, nearby K7V star orbited by two super-Earths and one mini-Neptune on close-in orbits. The system was first discovered using K2 data and then further characterized by other spectroscopic and photometric instruments. Previous literature studies provide several mass measurements for the three planets, however, with large variations and uncertainties. To better constrain the planetary masses, we added high-precision radial velocity measurements from ESPRESSO to published datasets from HARPS, HARPS-N, and HIRES and we performed a Gaussian process analysis combining radial velocity and photometric datasets from K2 and TESS. This method allowed us to model the stellar activity signal and derive precise planetary parameters. We determined planetary masses of $M_b = 4.28_{-0.33}^{+0.35}$ M${_\oplus}$, $M_c = 1.86_{-0.39}^{+0.37}$ M${_\oplus}$, and $M_d = 3.02_{-0.57}^{+0.58}$ M${_\oplus}$, and orbital periods of $1.208974 \pm 0.000001$ days for planet b, $3.648103_{-0.000010}^{+0.000013}$ days for planet c, and $6.201812 \pm 0.000009$ days for planet d. We compared our results to literature values and found that our derived uncertainties for the planetary mass, period, and radial velocity amplitude are smaller than the previously determined uncertainties. We modeled the interior composition of the three planets using the machine-learning-based tool ExoMDN and conclude that GJ 9827 b and c have an Earth-like composition, whereas GJ 9827 d has an hydrogen envelope, which, together with its density, places it in the mini-Neptune regime. △ Less

Submitted 16 January, 2024; v1 submitted 11 January, 2024; originally announced January 2024.

Comments: accepted in A&A

HYPERION. Coevolution of supermassive black holes and galaxies at $z>6$ and the build-up of massive galaxies

Authors: R. Tripodi, C. Feruglio, F. Fiore, L. Zappacosta, E. Piconcelli, M. Bischetti, A. Bongiorno, S. Carniani, F. Civano, C. -C. Chen, S. Cristiani, G. Cupani, F. Di Mascia, V. D'Odorico, X. Fan, A. Ferrara, S. Gallerani, M. Ginolfi, R. Maiolino, V. Mainieri, A. Marconi, I. Saccheo, F. Salvestrini, A. Tortosa, R. Valiante

Abstract: We used low- to high-frequency ALMA observations to investigate the cold gas and dust in ten QSOs at $z\gtrsim 6$. Our analysis of the CO(6-5) and CO(7-6) emission lines in the selected QSOs provided insights into their molecular gas masses, which average around $10^{10}\ \rm M_\odot$, consistent with typical values for high-redshift QSOs. Proprietary and archival ALMA observations in bands 8 and… ▽ More We used low- to high-frequency ALMA observations to investigate the cold gas and dust in ten QSOs at $z\gtrsim 6$. Our analysis of the CO(6-5) and CO(7-6) emission lines in the selected QSOs provided insights into their molecular gas masses, which average around $10^{10}\ \rm M_\odot$, consistent with typical values for high-redshift QSOs. Proprietary and archival ALMA observations in bands 8 and 9 enabled precise constraints on the dust properties and star formation rate (SFR) of four QSOs in our sample for the first time. The examination of the redshift distribution of dust temperatures revealed a general trend of increasing $T_{\rm dust}$ with redshift, which agrees with theoretical expectations. We computed a mean cold dust spectral energy distribution considering all ten QSOs. This offers a comprehensive view of the dust properties of high-$z$ QSOs. The QSOs marked by a more intense growth of the supermassive black hole (HYPERION QSOs) showed lower dust masses and higher gas-to-dust ratios on average, but their $\rm H_2$ gas reservoirs are consistent with those of other QSOs at the same redshift. The observed high SFR in our sample yields high SF efficiencies and thus very short gas depletion timescales ($τ_{\rm dep}\sim 10^{-2}$ Gyr). Beyond supporting the paradigm that high-$z$ QSOs reside in highly star-forming galaxies, our findings portrayed an interesting evolutionary path at $z>6$. Our study suggests that they are undergoing rapid galaxy growth that might be regulated by strong outflows. Their inferred evolutionary path shows a convergence toward the massive end of the local relation, which supports the idea that they are candidate progenitors of local massive galaxies. The observed pathway involves intense BH growth followed by substantial galaxy growth, in contrast with a symbiotic growth scenario. The abstract has been shortened (full version in the article). △ Less

Submitted 28 June, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 29 pages; 6 tables; 18 figures. Accepted by A&A. A section about SF efficiency has been added compared to the previous version

Boost recall in QSO selection from highly imbalanced photometric datasets

Authors: Giorgio Calderone, Francesco Guarneri, Matteo Porru, Stefano Cristiani, Andrea Grazian, Luciano Nicastro, Manuela Bischetti, Konstantina Boutsia, Guido Cupani, Valentina D'Odorico, Chiara Feruglio, Fabio Fontanot

Abstract: Context. The identification of bright QSOs is of great importance to probe the intergalactic medium and address open questions in cosmology. Several approaches have been adopted to find such sources in currently available photometric surveys, including machine learning methods. However, the rarity of bright QSOs at high redshifts compared to contaminating sources (such as stars and galaxies) makes… ▽ More Context. The identification of bright QSOs is of great importance to probe the intergalactic medium and address open questions in cosmology. Several approaches have been adopted to find such sources in currently available photometric surveys, including machine learning methods. However, the rarity of bright QSOs at high redshifts compared to contaminating sources (such as stars and galaxies) makes the selection of reliable candidates a difficult task, especially when high completeness is required. Aims. We present a novel technique to boost recall (i.e., completeness within the considered sample) in the selection of QSOs from photometric datasets dominated by stars, galaxies, and low-z QSOs (imbalanced datasets). Methods. Our method operates by iteratively removing sources whose probability of belonging to a noninteresting class exceeds a user-defined threshold, until the remaining dataset contains mainly high-z QSOs. Any existing machine learning method can be used as underlying classifier, provided it allows for a classification probability to be estimated. We applied the method to a dataset obtained by cross-matching PanSTARRS1, Gaia, and WISE, and identified the high-z QSO candidates using both our method and its direct multi-label counterpart. Results. We ran several tests by randomly choosing the training and test datasets, and achieved significant improvements in recall which increased from 50% to 85% for QSOs with z>2.5, and from 70% to 90% for QSOs with z>3. Also, we identified a sample of 3098 new QSO candidates on a sample of 2.6x10^6 sources with no known classification. We obtained follow-up spectroscopy for 121 candidates, confirming 107 new QSOs with z>2.5. Finally, a comparison of our candidates with those selected by an independent method shows that the two samples overlap by more than 90% and that both methods are capable of achieving a high level of completeness. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: Accepted for publication on A&A

Cosmology and fundamental physics with the ELT-ANDES spectrograph

Authors: C. J. A. P. Martins, R. Cooke, J. Liske, M. T. Murphy, P. Noterdaeme, T. M. Schmidt, J. S. Alcaniz, C. S. Alves, S. Balashev, S. Cristiani, P. Di Marcantonio, R. Génova Santos, R. S. Gonçalves, J. I. González Hernández, R. Maiolino, A. Marconi, C. M. J. Marques, M. A. F. Melo e Sousa, N. J. Nunes, L. Origlia, C. Péroux, S. Vinzl, A. Zanutta

Abstract: State-of-the-art 19th century spectroscopy led to the discovery of quantum mechanics, and 20th century spectroscopy led to the confirmation of quantum electrodynamics. State-of-the-art 21st century astrophysical spectrographs, especially ANDES at ESO's ELT, have another opportunity to play a key role in the search for, and characterization of, the new physics which is known to be out there, waitin… ▽ More State-of-the-art 19th century spectroscopy led to the discovery of quantum mechanics, and 20th century spectroscopy led to the confirmation of quantum electrodynamics. State-of-the-art 21st century astrophysical spectrographs, especially ANDES at ESO's ELT, have another opportunity to play a key role in the search for, and characterization of, the new physics which is known to be out there, waiting to be discovered. We rely on detailed simulations and forecast techniques to discuss four important examples of this point: big bang nucleosynthesis, the evolution of the cosmic microwave background temperature, tests of the universality of physical laws, and a real-time model-independent mapping of the expansion history of the universe (also known as the redshift drift). The last two are among the flagship science drivers for the ELT. We also highlight what is required for the ESO community to be able to play a meaningful role in 2030s fundamental cosmology and show that, even if ANDES only provides null results, such `minimum guaranteed science' will be in the form of constraints on key cosmological paradigms: these are independent from, and can be competitive with, those obtained from traditional cosmological probes. △ Less

Submitted 1 February, 2024; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 31 pages, 6 figures; submitted to Experimental Astronomy on behalf of the ANDES Science Team; v2: matches accepted version

The VANDELS ESO public spectroscopic survey: The spectroscopic measurements catalogue

Authors: M. Talia, C. Schreiber, B. Garilli, L. Pentericci, L. Pozzetti, G. Zamorani, F. Cullen, M. Moresco, A. Calabrò, M. Castellano, J. P. U. Fynbo, L. Guaita, F. Marchi, S. Mascia, R. McLure, M. Mignoli, E. Pompei, E. Vanzella, A. Bongiorno, G. Vietri, R. O. Amorín, M. Bolzonella, A. C. Carnall, A. Cimatti, G. Cresci , et al. (14 additional authors not shown)

Abstract: VANDELS is a deep spectroscopic survey, performed with the VIMOS instrument at VLT, aimed at studying in detail the physical properties of high-redshift galaxies. VANDELS targeted about 2100 sources at 1<z<6.5 in the CANDELS Chandra Deep-Field South (CDFS) and Ultra-Deep Survey (UDS) fields. In this paper we present the public release of the spectroscopic measurement catalogues from this survey, f… ▽ More VANDELS is a deep spectroscopic survey, performed with the VIMOS instrument at VLT, aimed at studying in detail the physical properties of high-redshift galaxies. VANDELS targeted about 2100 sources at 1<z<6.5 in the CANDELS Chandra Deep-Field South (CDFS) and Ultra-Deep Survey (UDS) fields. In this paper we present the public release of the spectroscopic measurement catalogues from this survey, featuring emission and absorption line centroids, fluxes, and rest-frame equivalent widths obtained through a Gaussian fit, as well as a number of atomic and molecular indices (e.g. Lick) and continuum breaks (e.g. D4000), and including a correction to be applied to the error spectra. We describe the measurement methods and the validation of the codes that were used. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: 12 pages, 9 figures; accepted for publication in Astronomy & Astrophysics. The measurement catalogues are accessible through the survey database (http://vandels.inaf.it) where all information can be queried interactively, and at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/cgi-bin/qcat?J/A+A/

Journal ref: A&A 678, A25 (2023)

On the 12C/13C isotopic ratio at the dawn of chemical evolution

Authors: P. Molaro, D. S. Aguado, E. Caffau, C. Allende Prieto, P. Bonifacio, J. I. Gonzalez Hernandez, R. Rebolo, M. R. Zapatero Osorio, S. Cristiani, F. Pepe, N. C. Santos, Y. Alibert, G. Cupani, P. Di Marcantonio, V. D'Odorico, C. Lovis, C. J. A. P. Martins, D. Milakovic, M. Murphy, N. J. Nunes, T. M. Schmidt, S. Sousa, a. Sozzetti, A. Suarez Mascareno

Abstract: The known Mega and Hyper Metal-Poor (MMP-HMP) stars with [Fe/H]<-6.0 and <-5.0, respectively, likely belong to the CEMP-no class, i.e. carbon-enhanced stars with low or absent second peak neutron capture elements. They are likely second generation stars and the few elements measurable in their atmospheres are used to infer the properties of single or very few progenitors. The high carbon abundance… ▽ More The known Mega and Hyper Metal-Poor (MMP-HMP) stars with [Fe/H]<-6.0 and <-5.0, respectively, likely belong to the CEMP-no class, i.e. carbon-enhanced stars with low or absent second peak neutron capture elements. They are likely second generation stars and the few elements measurable in their atmospheres are used to infer the properties of single or very few progenitors. The high carbon abundance in the CEMP-no stars offers a unique opportunity to measure the carbon isotopic ratio, which directly monitors the presence of mixing between the He and H-burning layers either within the star or in the progenitor(s). By means of high-resolution spectra acquired with the ESPRESSO spectrograph at the VLT we aim to derive values for the 12C/13C ratio at the lowest metallicities. A spectral synthesis technique based on the SYNTHE code and on ATLAS models is used within a Markov-chain Monte Carlo methodology to derive 12C/13C in the stellar atmospheres of five of the most metal poor stars. These are the Mega Metal-Poor giant SMS J0313-6708 ([Fe/H]<-7.1), the Hyper Metal-Poor dwarf HE1327-2326 ([Fe/H]=-5.8),the Hyper Metal-Poor giant SDSS J1313-0019 ([Fe/H] = -5.0) and the Ultra Metal-Poor subgiant HE0233-0343 ([Fe/H]=-4.7). We also revise a previous value for the Mega Metal-Poor giant SMSS J1605-1443 with ([Fe/H] = -6.2). In four stars we derive an isotopic value while for HE1327-2326 we provide a lower limit. All Measurements are in the range 39<12C/13C<100 showing that the He- and H-burning layers underwent partial mixing either in the stars or, more likely, in their progenitors. This provides evidence of a primary production of 13C at the dawn of chemical evolution. [abridged] △ Less

Submitted 2 October, 2023; v1 submitted 20 September, 2023; originally announced September 2023.

Comments: 10 pages, 6 figure, accepted A&A

Journal ref: A&A 679, A72 (2023)

A compact multi-planet system transiting HIP 29442 (TOI-469) discovered by TESS and ESPRESSO. Radial velocities lead to the detection of transits with low signal-to-noise ratio

Authors: M. Damasso, J. Rodrigues, A. Castro-González, B. Lavie, J. Davoult, M. R. Zapatero Osorio, J. Dou, S. G. Sousa, J. E. Owen, P. Sossi, V. Adibekyan, H. Osborn, Z. Leinhardt, Y. Alibert, C. Lovis, E. Delgado Mena, A. Sozzetti, S. C. C. Barros, D. Bossini, C. Ziegler, D. R. Ciardi, E. C. Matthews, P. J. Carter, J. Lillo-Box, A. Suárez Mascareño , et al. (30 additional authors not shown)

Abstract: We followed-up with ESPRESSO the K0V star HIP 29442 (TOI-469), already known to host a validated sub-Neptune companion TOI-469.01. We aim to verify the planetary nature of TOI-469.01. We modelled radial velocity and photometric time series to measure the dynamical mass, radius, and ephemeris, and to characterise the internal structure and composition of TOI-469.01. We confirmed the planetary natur… ▽ More We followed-up with ESPRESSO the K0V star HIP 29442 (TOI-469), already known to host a validated sub-Neptune companion TOI-469.01. We aim to verify the planetary nature of TOI-469.01. We modelled radial velocity and photometric time series to measure the dynamical mass, radius, and ephemeris, and to characterise the internal structure and composition of TOI-469.01. We confirmed the planetary nature of TOI-469.01. Thanks to ESPRESSO we discovered two additional close-in companions. We also detected their low signal-to-noise transit signals in the TESS light curve. HIP 29442 is a compact multi-planet system, and the three planets have orbital periods $P_{\rm orb, b}=13.63083\pm0.00003$, $P_{\rm orb, c}=3.53796\pm0.00003$, and $P_{\rm orb, d}=6.42975^{+0.00009}_{-0.00010}$ days, and we measured their masses with high precision: $m_{\rm p, b}=9.6\pm0.8~M_{\oplus}$, $m_{\rm p, c}=4.5\pm0.3~M_{\oplus}$, and $m_{\rm p, d}=5.1\pm0.4~M_{\oplus}$. We measured radii and bulk densities of all the planets (the 3$σ$ confidence intervals are shown in parenthesis): $R_{\rm p, b}=3.48^{+0.07 (+0.19)}_{-0.08 (-0.28)} ~R_{\oplus}$ and $ρ_{\rm p, b}=1.3\pm0.2 (0.3) g~cm^{-3}$; $R_{\rm p, c}=1.58^{+0.10 (+0.30)}_{-0.11 (-0.34)}~R_{\oplus}$ and $ρ_{\rm p, c}=6.3^{+1.7 (+6.0)}_{-1.3 (-2.7)} g~cm^{-3}$; $R_{\rm p, d}=1.37\pm0.11^{(+0.32)}_{(-0.43)}~R_{\oplus}$ and $ρ_{\rm p, d}=11.0^{+3.4 (+21.0)}_{-2.4 (-6.3)} g~cm^{-3}$. We used the more conservative 3$σ$ confidence intervals for the radii as input to the interior structure modelling. We find that HIP 29442 $b$ appears as a typical sub-Neptune, likely surrounded by a gas layer of pure H-He with a mass of $0.27^{+0.24}_{-0.17} M_{\oplus}$ and a thickness of $1.4\pm0.5 R_{\oplus}$. For the innermost companions HIP 29442 $c$ HIP 29442 $d$, the model supports an Earth-like composition. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: Accepted for publication on A&A

The MUSE Ultra Deep Field (MUDF). IV. A pair of X-ray weak quasars at the heart of two extended Lyα nebulae

Authors: Elisabeta Lusso, Emanuele Nardini, Michele Fumagalli, Matteo Fossati, Fabrizio Arrigoni Battaia, Mitchell Revalski, Marc Rafelski, Valentina D'Odorico, Celine Peroux, Stefano Cristiani, Pratika Dayal, Francesco Haardt, Emma K. Lofthouse

Abstract: We present the results obtained from follow-up observations of the MUSE Ultra Deep Field (MUDF) at X-ray energies with XMM-Newton. The MUDF is centred on a unique field with two bright, physically associated quasars at $z\simeq3.23$, separated by $\sim$500 kpc in projection. Both quasars are embedded within extended Ly$α$ nebulae ($\gtrsim 100~\rm kpc$ at a surface brightness flux level of… ▽ More We present the results obtained from follow-up observations of the MUSE Ultra Deep Field (MUDF) at X-ray energies with XMM-Newton. The MUDF is centred on a unique field with two bright, physically associated quasars at $z\simeq3.23$, separated by $\sim$500 kpc in projection. Both quasars are embedded within extended Ly$α$ nebulae ($\gtrsim 100~\rm kpc$ at a surface brightness flux level of $\approx 6\times 10^{-19} \rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$), whose elongated morphology is suggestive of an extended filament connecting the quasar haloes. The new X-ray observations presented here allow us to characterise the physical properties (e.g. X-ray slope, luminosities, gas column densities) in the innermost region of the MUDF quasars. We find that both quasars are X-ray underluminous compared to objects at similar ultraviolet luminosities. Based on our X-ray spectral analysis, absorbing columns of $N_H(z)\gtrsim$ 10$^{23}$ cm$^{-2}$ appear unlikely, therefore such a weakness is possibly intrinsic. When also including literature data, we do not observe any detectable trend between the area of the nebulae and nuclear luminosities at both the rest-frame 2 keV and 2500 $\rm \mathring{A}$. The area is also not correlated with the X-ray photon index nor with the integrated band flux in the hard band (2$-$10 keV). We also do not find any trend between the extended Ly$α$ emission of the nebulae and the nuclear X-ray luminosity. Finally, the properties of the MUDF quasars' nebulae are consistent with the observed relation between the Ly$α$ integrated luminosity of the nebulae and their area. Our results suggest that the quasar ionization power is not a strong driver of the morphology and size of the nebulae. △ Less

Submitted 4 September, 2023; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: 15 pages, 9 figures, reference added, published in MNRAS

Crossing the Rubicon of Reionization with z~5 QSOs

Authors: A. Grazian, K. Boutsia, E. Giallongo, S. Cristiani, F. Fontanot, M. Bischetti, A. Bongiorno, G. Calderone, G. Cupani, V. D'Odorico, C. Feruglio, F. Fiore, F. Guarneri, M. Porru, I. Saccheo

Abstract: One of the key open questions in Cosmology is the nature of the sources that completed the cosmological hydrogen Reionization at z~5.2. High-z primeval galaxies have been long considered the main drivers for Reionization, with a minor role played by high-z AGN. However, in order to confirm this scenario, it is fundamental to measure the photo-ionization rate produced by active SMBHs close to the e… ▽ More One of the key open questions in Cosmology is the nature of the sources that completed the cosmological hydrogen Reionization at z~5.2. High-z primeval galaxies have been long considered the main drivers for Reionization, with a minor role played by high-z AGN. However, in order to confirm this scenario, it is fundamental to measure the photo-ionization rate produced by active SMBHs close to the epoch of Reionization. Given the pivotal role played by spectroscopically complete observations of high-z QSOs, in this paper we present the first results of the RUBICON (Reionizing the Universe with BrIght COsmological Nuclei) survey. It consists of a color selected sample of bona-fide z~5 QSO candidates from the Hyper Suprime-Cam Subaru Strategic Survey. Our QSO candidates have been validated both by photometric redshifts based on SED fitting and by spectroscopic redshifts, confirming that they lie at 4.5<z_spec<5.2. A relatively large space density of QSOs (Phi~1.4x10^-8 cMpc^-3) is thus confirmed at z~5 and M1450~-27, consistent with a pure density evolution of the AGN luminosity function from z=4 to z=5, with a mild density evolution rate of 0.25 dex. This indicates that AGN could play a non-negligible role in the cosmic Reionization. The Rubicon of Reionization has been crossed. △ Less

Submitted 3 August, 2023; v1 submitted 23 July, 2023; originally announced July 2023.

Comments: Accepted for publication in ApJ, 25 pages, 11 figures, 3 tables, figures 1, 3, 5, and 6 resized

Advanced Data Analysis for Observational Cosmology: applications to the study of the Intergalactic Medium

Authors: Guido Cupani, Giorgio Calderone, Stefano Cristiani, Francesco Guarneri

Abstract: The analysis of absorption features along the line of sight to distant sources is an invaluable tool for observational cosmology, giving a direct insight into the physical and chemical state of the inter/circumgalactic medium. Such endeavour entails the accessibility of bright QSOs as background beacons, and the availability of software tools to extract the information in a reproducible way. In th… ▽ More The analysis of absorption features along the line of sight to distant sources is an invaluable tool for observational cosmology, giving a direct insight into the physical and chemical state of the inter/circumgalactic medium. Such endeavour entails the accessibility of bright QSOs as background beacons, and the availability of software tools to extract the information in a reproducible way. In this article, we will present the latest results we obtained in both directions within the QUBRICS project: we will describe how machine learning techniques were applied to detect hundreds of previously unknown QSOs in the southern hemisphere, and how state-of-the art software like QSFit and Astrocook was integrated in the analysis of the targets, opening up new possibilities for the next era of observations. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Comments: 4 pages, 2 figures; proceedings of ADASS XXXI, accepted by ASP Conference Series

XQR-30: the ultimate XSHOOTER quasar sample at the reionization epoch

Authors: Valentina D'Odorico, E. Banados, G. D. Becker, M. Bischetti, S. E. I. Bosman, G. Cupani, R. Davies, E. P. Farina, A. Ferrara, C. Feruglio, C. Mazzucchelli, E. Ryan-Weber, J. -T. Schindler, A. Sodini, B. P. Venemans, F. Walter, H. Chen, S. Lai, Y. Zhu, F. Bian, S. Campo, S. Carniani, S. Cristiani, F. Davies, R. Decarli , et al. (24 additional authors not shown)

Abstract: The final phase of the reionization process can be probed by rest-frame UV absorption spectra of quasars at z>6, shedding light on the properties of the diffuse intergalactic medium within the first Gyr of the Universe. The ESO Large Programme "XQR-30: the ultimate XSHOOTER legacy survey of quasars at z~5.8-6.6" dedicated ~250 hours of observations at the VLT to create a homogeneous and high-quali… ▽ More The final phase of the reionization process can be probed by rest-frame UV absorption spectra of quasars at z>6, shedding light on the properties of the diffuse intergalactic medium within the first Gyr of the Universe. The ESO Large Programme "XQR-30: the ultimate XSHOOTER legacy survey of quasars at z~5.8-6.6" dedicated ~250 hours of observations at the VLT to create a homogeneous and high-quality sample of spectra of 30 luminous quasars at z~6, covering the rest wavelength range from the Lyman limit to beyond the MgII emission. Twelve quasar spectra of similar quality from the XSHOOTER archive were added to form the enlarged XQR-30 sample, corresponding to a total of ~350 hours of on-source exposure time. The median effective resolving power of the 42 spectra is R~11400 and 9800 in the VIS and NIR arm, respectively. The signal-to-noise ratio per 10 km/s pixel ranges from ~11 to 114 at $λ\simeq 1285$ Årest frame, with a median value of ~29. We describe the observations, data reduction and analysis of the spectra, together with some first results based on the E-XQR-30 sample. New photometry in the H and K bands are provided for the XQR-30 quasars, together with composite spectra whose characteristics reflect the large absolute magnitudes of the sample. The composite and the reduced spectra are released to the community through a public repository, and will enable a range of studies addressing outstanding questions regarding the first Gyr of the Universe. △ Less

Submitted 19 June, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 21 pages, 10 figures. Final version accepted by MNRAS

Journal ref: MNRAS, Volume 523, Issue 1, pp.1399-1420 (2023)

An unusually low-density super-Earth transiting the bright early-type M-dwarf GJ 1018 (TOI-244)

Authors: A. Castro-González, O. D. S. Demangeon, J. Lillo-Box, C. Lovis, B. Lavie, V. Adibekyan, L. Acuña, M. Deleuil, A. Aguichine, M. R. Zapatero Osorio, H. M. Tabernero, J. Davoult, Y. Alibert, N. Santos, S. G. Sousa, A. Antoniadis-Karnavas, F. Borsa, J. N. Winn, C. Allende Prieto, P. Figueira, J. M. Jenkins, A. Sozzetti, M. Damasso, A. M. Silva, N. Astudillo-Defru , et al. (12 additional authors not shown)

Abstract: Small planets located at the lower mode of the bimodal radius distribution are generally assumed to be composed of iron and silicates in a proportion similar to that of the Earth. However, recent discoveries are revealing a new group of low-density planets that are inconsistent with that description. We intend to confirm and characterize the TESS planet candidate TOI-244.01, which orbits the brigh… ▽ More Small planets located at the lower mode of the bimodal radius distribution are generally assumed to be composed of iron and silicates in a proportion similar to that of the Earth. However, recent discoveries are revealing a new group of low-density planets that are inconsistent with that description. We intend to confirm and characterize the TESS planet candidate TOI-244.01, which orbits the bright ($K$ = 7.97 mag), nearby ($d$ = 22 pc), and early-type (M2.5 V) M-dwarf star GJ 1018 with an orbital period of 7.4 days. We used Markov Chain Monte Carlo methods to model 57 precise radial velocity measurements acquired by the ESPRESSO spectrograph together with TESS photometry and complementary HARPS data. We find TOI-244 b to be a super-Earth with a radius of $R_{\rm p}$ = 1.52 $\pm$ 0.12 $\rm R_{\oplus}$ and a mass of $M_{\rm p}$ = 2.68 $\pm$ 0.30 $\rm M_{\oplus}$. These values correspond to a density of $ρ$ = 4.2 $\pm$ 1.1 $\rm g \cdot cm^{-3}$, which is below what would be expected for an Earth-like composition. We find that atmospheric loss processes may have been efficient to remove a potential primordial hydrogen envelope, but high mean molecular weight volatiles such as water could have been retained. Our internal structure modeling suggests that TOI-244 b has a $479^{+128}_{-96}$ km thick hydrosphere over a 1.17 $\pm$ 0.09 $\rm R_{\oplus}$ solid structure composed of a Fe-rich core and a silicate-dominated mantle compatible with that of the Earth. On a population level, we find two tentative trends in the density-metallicity and density-insolation parameter space for the low-density super-Earths, which may hint at their composition. With a 8$\%$ precision in radius and 12$\%$ precision in mass, TOI-244 b is among the most precisely characterized super-Earths, which, together with the likely presence of an extended hydrosphere, makes it a key target for atmospheric observations. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 675, A52 (2023)

HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). A new regime for the X-ray nuclear properties of the first quasars

Authors: L. Zappacosta, E. Piconcelli, F. Fiore, I. Saccheo, R. Valiante, C. Vignali, F. Vito, M. Volonteri, M. Bischetti, A. Comastri, C. Done, M. Elvis, E. Giallongo, F. La Franca, G. Lanzuisi, M. Laurenti, G. Miniutti, A. Bongiorno, M. Brusa, F. Civano, S. Carniani, V. D'Odorico, C. Feruglio, S. Gallerani, R. Gilli , et al. (18 additional authors not shown)

Abstract: The existence of luminous quasars (QSO) at the Epoch of Reionization (EoR; i.e. z>6) powered by supermassive black holes (SMBH) with masses $\gtrsim10^9~M_\odot$ challenges models of early SMBH formation. To shed light on the nature of these sources we started a multiwavelength programme based on a sample of 18 HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). These are the luminous Q… ▽ More The existence of luminous quasars (QSO) at the Epoch of Reionization (EoR; i.e. z>6) powered by supermassive black holes (SMBH) with masses $\gtrsim10^9~M_\odot$ challenges models of early SMBH formation. To shed light on the nature of these sources we started a multiwavelength programme based on a sample of 18 HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). These are the luminous QSOs whose SMBH must have had the fastest mass growth during the Universe first Gyr. In this paper we present the HYPERION sample and report on the first of the 3 years planned observations of the 2.4 Ms XMM-Newton Multi-Year Heritage program on which HYPERION is based. The goal of this program is to accurately characterize the X-ray nuclear properties of QSOs at the EoR. Through a joint X-ray spectral analysis of 10 sources, in the rest-frame $\sim2-50$ keV range, we report a steep average photon index ($Γ\sim2.4\pm0.1$). Absorption is not required. The average $Γ$ is inconsistent at $\geq4σ$ level with the canonical 1.8-2 value measured in QSO at z<6. This spectral slope is also much steeper than that reported in lower-z QSOs with similar luminosity or accretion rate, thus suggesting a genuine redshift evolution. Alternatively, we can interpret this result as the presence of an unusually low-energy cutoff $E_{cut}\sim20$ keV on a standard $Γ=1.9$ power-law. We also report on mild indications that HYPERION QSOs show higher soft X-ray emission at 2 keV compared to the UV one at 2500A than expected by lower-z luminous AGN. We speculate that a redshift-dependent coupling between the corona and accretion disc or intrinsically different coronal properties may account for the steep spectral slopes, especially in the presence of powerful winds. The reported slopes, if confirmed at lower luminosities, may have an important impact on future X-ray AGN studies in the early Universe. △ Less

Submitted 19 July, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 21 pages (including appendix), 12 figures, 4 tables. Accepted for pubblication in A&A

Journal ref: A&A 678, A201 (2023)

Spectroscopy of QUBRICS quasar candidates: 1672 new redshifts and a Golden Sample for the Sandage Test of the Redshift Drift

Authors: Stefano Cristiani, Matteo Porru, Francesco Guarneri, Giorgio Calderone, Konstantina Boutsia, Andrea Grazian, Guido Cupani, Valentina D'Odorico, Fabio Fontanot, Carlos J. A. P. Martins, Catarina M. J. Marques, Soumak Maitra, Andrea Trost

Abstract: The QUBRICS (QUasars as BRIght beacons for Cosmology in the Southern hemisphere) survey aims at constructing a sample of the brightest quasars with z>~2.5, observable with facilities in the Southern Hemisphere. QUBRICS makes use of the available optical and IR wide-field surveys in the South and of Machine Learning techniques to produce thousands of bright quasar candidates of which only a few hun… ▽ More The QUBRICS (QUasars as BRIght beacons for Cosmology in the Southern hemisphere) survey aims at constructing a sample of the brightest quasars with z>~2.5, observable with facilities in the Southern Hemisphere. QUBRICS makes use of the available optical and IR wide-field surveys in the South and of Machine Learning techniques to produce thousands of bright quasar candidates of which only a few hundred have been confirmed with follow-up spectroscopy. Taking advantage of the recent Gaia Data Release 3, which contains 220 million low-resolution spectra, and of a newly developed spectral energy distribution fitting technique, designed to combine the photometric information with the Gaia spectroscopy, it has been possible to measure 1672 new secure redshifts of QUBRICS candidates, with a typical uncertainty $σ_z = 0.02$. This significant progress of QUBRICS brings it closer to (one of) its primary goals: providing a sample of bright quasars at redshift 2.5 < z < 5 to perform the Sandage test of the cosmological redshift drift. A Golden Sample of seven quasars is presented that makes it possible to carry out this experiment in about 1500 hours of observation in 25 years, using the ANDES spectrograph at the 39m ELT, a significant improvement with respect to previous estimates. △ Less

Submitted 1 April, 2023; originally announced April 2023.

Comments: 11 pages, 10 figures, accepted for publication in MNRAS

Spectrographs and Spectroscopists for the Sandage Test

Authors: S. Cristiani, K. Boutsia, G. Calderone, G. Cupani, V. D'Odorico, F. Fontanot, A. Grazian, F. Guarneri, C. Martins, L. Pasquini, M. Porru, E. Vanzella

Abstract: The redshift drift is a small, dynamic change in the redshift of objects following the Hubble flow. Its measurement provides a direct, real-time, model-independent mapping of the expansion rate of the Universe. It is fundamentally different from other cosmological probes: instead of mapping our (present-day) past light-cone, it directly compares different past light-cones. Being independent of any… ▽ More The redshift drift is a small, dynamic change in the redshift of objects following the Hubble flow. Its measurement provides a direct, real-time, model-independent mapping of the expansion rate of the Universe. It is fundamentally different from other cosmological probes: instead of mapping our (present-day) past light-cone, it directly compares different past light-cones. Being independent of any assumptions on gravity, geometry or clustering, it directly tests the pillars of the Lambda-CDM paradigm. Recent theoretical studies have uncovered unique synergies with other cosmological probes, including the characterization of the physical properties of dark energy. At the time of the original proposal by Sandage (1962) the expected change in the redshift of objects at cosmological distances appeared to be exceedingly small for reasonable observing times and beyond technological capabilities. In the last decades progress in the spectrographs (e.g. ESPRESSO), in the collecting area of telescopes and in the samples of cosmic beacons, enabled by new datasets and new machine-learning-based selections, have drastically changed the situation, bringing the Redshift Drift Grail within reach. As a consequence, this measurement is a flagship objective of the Extremely Large Telescope (ELT), specifically of its high-resolution spectrograph, ANDES. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: Proceedings of the Conference "Past, Present and Future of Astrophysical Spectroscopy", Trieste, 6-10 June 2022; accepted for publication; 8 pages, 4 figures

Journal ref: Mem.S.A.It., Proceedings of the Conference "Past, Present and Future of Astrophysical Spectroscopy", Trieste, 6-10 June 2022

Eddington accreting Black Holes in the Epoch of Reionization

Authors: Fabio Fontanot, Stefano Cristiani, Andrea Grazian, Francesco Haardt, Valentina D'Odorico, Konstantina Boutsia, Giorgio Calderone, Guido Cupani, Francesco Guarneri, Chiara Fiorin, Giulia Rodighiero

Abstract: The evolution of the luminosity function (LF) of Active Galactic Nuclei (AGNs) at $z \gtrsim 5$ represents a key constraint to understand their contribution to the ionizing photon budget necessary to trigger the last phase transition in the Universe, i.e. the epoch of Reionization. Recent searches for bright high-z AGNs suggest that the space densities of this population at $z>4$ has to be revised… ▽ More The evolution of the luminosity function (LF) of Active Galactic Nuclei (AGNs) at $z \gtrsim 5$ represents a key constraint to understand their contribution to the ionizing photon budget necessary to trigger the last phase transition in the Universe, i.e. the epoch of Reionization. Recent searches for bright high-z AGNs suggest that the space densities of this population at $z>4$ has to be revised upwards, and sparks new questions about their evolutionary paths. Gas accretion is the key physical mechanism to understand both the distribution of luminous sources and the growth of central Super-Massive Black Holes (SMBHs). In this work, we model the high-z AGN-LF assuming that high-z luminous AGN shine at their Eddington limit: we derive the expected evolution as a function of the ``duty-cycle'' ($f_{\rm dc}$), i.e. the fraction of life-time that a given SMBH spends accreting at the Eddington rate. Our results show that intermediate values ($f_{\rm dc} \simeq 0.1$) predict the best agreement with the ionizing background and photoionization rate, but do not provide enough ionizing photons to account for the observed evolution of the hydrogen neutral fraction. Smaller values ($f_{\rm dc} \lesssim 0.05$) are required for AGNs to be the dominant population responsible for Hydrogen reionization in the Early Universe. We then show that this low-$f_{\rm dc}$ evolution can be reconciled with the current constraints on Helium reionization, although it implies a relatively large number of inactive SMBHs at $z\gtrsim5$, in tension with SMBH growth models based on heavy seeding. △ Less

Submitted 17 January, 2023; originally announced January 2023.

Comments: 11 pages, 5 figures, MNRAS accepted

The pristine nature of SMSS 1605$-$1443 revealed by ESPRESSO

Authors: D. S. Aguado, E. Caffau, P. Molaro, C. Allende Prieto, P. Bonifacio, J. I. González Hernández, R. Rebolo, S. Salvadori, M. R. Zapatero Osorio, S. Cristiani, F. Pepe, C. Santos, G. Cupani, P. Di Marcantonio, V. D'Odorico, C. Lovis, N. J. Nunes, C. J. A. P. Martins, D. Milakovic, J. Rodrigues, T. M. Schmidt, A. Sozzetti, A. Suarez Mascareño

Abstract: SMSS J160540.18$-$144323.1 is the carbon-enhanced metal-poor (CEMP) star with the lowest iron abundance ever measured, [Fe/H]=-6.2, which was first reported with the SkyMapper telescope. The carbon abundance is A(C)~6.1 in the low-C band, as the majority of the stars in this metallicity range. Yet, constraining the isotopic ratio of key species, such as carbon, sheds light on the properties and or… ▽ More SMSS J160540.18$-$144323.1 is the carbon-enhanced metal-poor (CEMP) star with the lowest iron abundance ever measured, [Fe/H]=-6.2, which was first reported with the SkyMapper telescope. The carbon abundance is A(C)~6.1 in the low-C band, as the majority of the stars in this metallicity range. Yet, constraining the isotopic ratio of key species, such as carbon, sheds light on the properties and origin of these elusive stars. We performed high-resolution observations of SMSS1605$-$1443 with the ESPRESSO spectrograph to look for variations in the radial velocity ($v_{rad}$) with time. These data have been combined with older MIKE and UVES archival observations to enlarge the temporal baseline. The $^{12}$C/$^{13}$C isotopic ratio is also studied to explore the possibility of mass transfer from a binary companion. A cross-correlation function against a natural template was applied to detect $v_{rad}$ variability and a spectral synthesis technique was used to derive $^{12}$C/$^{13}$C in the stellar atmosphere. We confirm previous indications of binarity in SMSS1605$-$1443 and measured a lower limit $^{12}$C/$^{13}$C$>60$ at more than a 3$σ$ confidence level, proving that this system is chemically unmixed and that no mass transfer from the unseen companion has happened so far. Thus, we confirm the CEMP-no nature of SMSS1605$-$1443 and show that the pristine chemical composition of the cloud from which it formed is currently imprinted in its stellar atmosphere free of contamination. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: Published in A&A, 9 pages, 6 figures

CUBES: a UV spectrograph for the future

Authors: S. Covino, S. Cristiani, J. M. Alcala', S. H. P. Alencar, S. A. Balashev, B. Barbuy, N. Bastian, U. Battino, L. Bissell, P. Bristow, A. Calcines, G. Calderone, P. Cambianica, R. Carini, B. Carter, S. Cassisi, B. V. Castilho, G. Cescutti, N. Christlieb, R. Cirami, R. Conzelmann, I. Coretti, R. Cooke, G. Cremonese, K. Cunha , et al. (64 additional authors not shown)

Abstract: In spite of the advent of extremely large telescopes in the UV/optical/NIR range, the current generation of 8-10m facilities is likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral r… ▽ More In spite of the advent of extremely large telescopes in the UV/optical/NIR range, the current generation of 8-10m facilities is likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R>20,000, although a lower-resolution, sky-limited mode of R ~ 7,000 is also planned. CUBES will offer new possibilities in many fields of astrophysics, providing access to key lines of stellar spectra: a tremendous diversity of iron-peak and heavy elements, lighter elements (in particular Beryllium) and light-element molecules (CO, CN, OH), as well as Balmer lines and the Balmer jump (particularly important for young stellar objects). The UV range is also critical in extragalactic studies: the circumgalactic medium of distant galaxies, the contribution of different types of sources to the cosmic UV background, the measurement of H2 and primordial Deuterium in a regime of relatively transparent intergalactic medium, and follow-up of explosive transients. The CUBES project completed a Phase A conceptual design in June 2021 and has now entered the Phase B dedicated to detailed design and construction. First science operations are planned for 2028. In this paper, we briefly describe the CUBES project development and goals, the main science cases, the instrument design and the project organization and management. △ Less

Submitted 24 December, 2022; originally announced December 2022.

Comments: Proceedings for the HACK100 conference, Trieste, June 2022. arXiv admin note: substantial text overlap with arXiv:2208.01672

Two temperate Earth-mass planets orbiting the nearby star GJ1002

Authors: A. Suárez Mascareño, E. González-Álvarez, M. R. Zapatero Osorio, J. Lillo-Box, J. P. Faria, V. M. Passegger, J. I. González Hernández, P. Figueira, A. Sozzetti, R. Rebolo, F. Pepe, N. C. Santos, S. Cristiani, C. Lovis, A. M. Silva, I. Ribas, P. J. Amado, J. A. Caballero, A. Quirrenbach, A. Reiners, M. Zechmeister, V. Adibekyan, Y. Alibert, V. J. S. Béjar, S. Benatti , et al. (20 additional authors not shown)

Abstract: We report the discovery and characterisation of two Earth-mass planets orbiting in the habitable zone of the nearby M-dwarf GJ~1002 based on the analysis of the radial-velocity (RV) time series from the ESPRESSO and CARMENES spectrographs. The host star is the quiet M5.5~V star GJ~1002 (relatively faint in the optical, $V \sim 13.8$ mag, but brighter in the infrared, $J \sim 8.3$ mag), located at… ▽ More We report the discovery and characterisation of two Earth-mass planets orbiting in the habitable zone of the nearby M-dwarf GJ~1002 based on the analysis of the radial-velocity (RV) time series from the ESPRESSO and CARMENES spectrographs. The host star is the quiet M5.5~V star GJ~1002 (relatively faint in the optical, $V \sim 13.8$ mag, but brighter in the infrared, $J \sim 8.3$ mag), located at 4.84 pc from the Sun. We analyse 139 spectroscopic observations taken between 2017 and 2021. We performed a joint analysis of the time series of the RV and full-width half maximum (FWHM) of the cross-correlation function (CCF) to model the planetary and stellar signals present in the data, applying Gaussian process regression to deal with the stellar activity. We detect the signal of two planets orbiting GJ~1002. GJ~1002~b is a planet with a minimum mass $m_p \sin i $ of 1.08 $\pm$ 0.13 M$_{\oplus}$ with an orbital period of 10.3465 $\pm$ 0.0027 days at a distance of 0.0457 $\pm$ 0.0013 au from its parent star, receiving an estimated stellar flux of 0.67 $F_{\oplus}$. GJ~1002 c is a planet with a minimum mass $m_p \sin i $ of 1.36 $\pm$ 0.17 M$_{\oplus}$ with an orbital period of 21.202 $\pm$ 0.013 days at a distance of 0.0738 $\pm$ 0.0021 au from its parent star, receiving an estimated stellar flux of 0.257 $F_{\oplus}$. We also detect the rotation signature of the star, with a period of 126 $\pm$ 15 days. GJ~1002 is one of the few known nearby systems with planets that could potentially host habitable environments. The closeness of the host star to the Sun makes the angular sizes of the orbits of both planets ($\sim$ 9.7 mas and $\sim$ 15.7 mas, respectively) large enough for their atmosphere to be studied via high-contrast high-resolution spectroscopy with instruments such as the future spectrograph ANDES for the ELT or the LIFE mission. △ Less

Submitted 25 January, 2023; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 23 pages, 20 figures, 3 tables ; Replacement to fix some typos

Planetary system around LTT 1445A unveiled by ESPRESSO: Multiple planets in a triple M-dwarf system

Authors: B. Lavie, F. Bouchy, C. Lovis, M. Zapatero Osorio, A. Deline, S. Barros, P. Figueira, A. Sozzetti, J. I. Gonzalez Hernandez, J. Lillo-Box, J. Rodrigues, A. Mehner, M. Damasso, V. Adibekyan, Y. Alibert, C. Allende Prieto, S. Cristiani, V. DOdorico, P. Di Marcantonio, D. Ehrenreich, R. Genova Santos, G. Lo Curto, C. J. A. P. Martins, G. Micela, P. Molaro , et al. (10 additional authors not shown)

Abstract: We present radial velocity follow-up obtained with ESPRESSO of the M-type star LTT 1445A (TOI-455), for which a transiting planet b with an orbital period of~5.4 days was detected by TESS. We report the discovery of a second transiting planet (LTT 1445A c) and a third non-transiting candidate planet (LTT 1445A d) with orbital periods of 3.12 and 24.30 days, respectively. The host star is the main… ▽ More We present radial velocity follow-up obtained with ESPRESSO of the M-type star LTT 1445A (TOI-455), for which a transiting planet b with an orbital period of~5.4 days was detected by TESS. We report the discovery of a second transiting planet (LTT 1445A c) and a third non-transiting candidate planet (LTT 1445A d) with orbital periods of 3.12 and 24.30 days, respectively. The host star is the main component of a triple M-dwarf system at a distance of 6.9 pc. We used 84 ESPRESSO high-resolution spectra to determine accurate masses of 2.3$\pm$0.3 $\mathrm{M}_\oplus$ and 1.0$\pm$0.2 $\mathrm{M}_\oplus$ for planets b and c and a minimum mass of 2.7$\pm$0.7 $\mathrm{M}_\oplus$ for planet d. Based on its radius of 1.43$\pm0.09$ $\mathrm{R}_\oplus$ as derived from the TESS observations, LTT 1445A b has a lower density than the Earth and may therefore hold a sizeable atmosphere, which makes it a prime target for the James Webb Space Telescope. We used a Bayesian inference approach with the nested sampling algorithm and a set of models to test the robustness of the retrieved physical values of the system. There is a probability of 85$\%$ that the transit of planet c is grazing, which results in a retrieved radius with large uncertainties at 1.60$^{+0.67}_{-0.34}$ $\mathrm{R}_\oplus$. LTT 1445A d orbits the inner boundary of the habitable zone of its host star and could be a prime target for the James Webb Space Telescope. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 31 pages, 20 figures Accepted A&A

Journal ref: A&A 673, A69 (2023)

Detection of barium in the atmospheres of the ultra-hot gas giants WASP-76b and WASP-121b

Authors: T. Azevedo Silva, O. D. S. Demangeon, N. C. Santos, R. Allart, F. Borsa, E. Cristo, E. Esparza-Borges, J. V. Seidel, E. Palle, S. G. Sousa, H. M. Tabernero, M. R. Zapatero Osorio, S. Cristiani, F. Pepe, R. Rebolo, V. Adibekyan, Y. Alibert, S. C. C. Barros, F. Bouchy, V. Bourrier, G. Lo Curto, P. Di Marcantonio, V. D'Odorico, D. Ehrenreich, P. Figueira , et al. (11 additional authors not shown)

Abstract: High-resolution spectroscopy studies of ultra-hot Jupiters have been key in our understanding of exoplanet atmospheres. Observing into the atmospheres of these giant planets allows for direct constraints on their atmospheric compositions and dynamics while laying the groundwork for new research regarding their formation and evolution environments. Two of the most well-studied ultra-hot Jupiters ar… ▽ More High-resolution spectroscopy studies of ultra-hot Jupiters have been key in our understanding of exoplanet atmospheres. Observing into the atmospheres of these giant planets allows for direct constraints on their atmospheric compositions and dynamics while laying the groundwork for new research regarding their formation and evolution environments. Two of the most well-studied ultra-hot Jupiters are WASP-76b and WASP-121b, with multiple detected chemical species and strong signatures of their atmospheric dynamics. We take a new look at these two exceptional ultra-hot Jupiters by reanalyzing the transit observations taken with ESPRESSO at the Very Large Telescope and attempt to detect additional species. To extract the planetary spectra of the two targets, we corrected for the telluric absorption and removed the stellar spectrum contributions. We then exploited new synthetic templates that were specifically designed for ultra-hot Jupiters in combination with the cross-correlation technique to unveil species that remained undetected by previous analyses. We add a novel detection of Ba+ to the known atmospheric compositions of WASP-76b and WASP-121b, the heaviest species detected to date in any exoplanetary atmosphere, with additional new detections of Co and Sr+ and a tentative detection of Ti+ for WASP-121b. We also confirm the presence of Ca+, Cr, Fe, H, Li, Mg, Mn, Na, and V on both WASP-76b and WASP-121b, with the addition of Ca, Fe+, and Ni for the latter. Finally, we also confirm the clear asymmetric absorption feature of Ca+ on WASP-121b, with an excess absorption at the bluer wavelengths and an effective planet radius beyond the Roche lobe. This indicates that the signal may arise from the escape of planetary atmosphere. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: 21 pages, 16 figures. Published in A&A Volume 666, October 2022 - Letter to the Editor

Journal ref: A&A 666, L10 (2022)

Espresso observations of HE 0107$-$5240 and other CEMP-no stars with $\rm [Fe/H]\le -4.5$

Authors: D. Aguado, P. Molaro, E. Caffau, J. I. González Hernández, M. Zapatero Osorio, P. Bonifacio, C. Allende Prieto, R. Rebolo, M. Damasso, A. Suárez Mascareño, S. B. Howell, E. Furlan, S. Cristiani, G. Cupani, P. Di Marcantonio, V. D'Odorico, C. Lovis, C. J. A. P. Martins, D. Milakovic, M. T. Murphy, N. J. Nunes, F. Pepe, N. C. Santos, T. M. Schmidt, A. Sozzetti

Abstract: HE 0107$-$5240 is a hyper metal-poor star with $\rm [Fe/H]=-5.39$. We performed high-res observations with the ESPRESSO spectrograph at the VLT to constrain the kinematical properties of the binary system HE 0107$-$5240 and to probe the binarity of the sample of 8 most metal-poor stars with $\rm [Fe/H]<-4.5$. Radial velocities are obtained by using cross-correlation in the interval 4200$-$4315A, w… ▽ More HE 0107$-$5240 is a hyper metal-poor star with $\rm [Fe/H]=-5.39$. We performed high-res observations with the ESPRESSO spectrograph at the VLT to constrain the kinematical properties of the binary system HE 0107$-$5240 and to probe the binarity of the sample of 8 most metal-poor stars with $\rm [Fe/H]<-4.5$. Radial velocities are obtained by using cross-correlation in the interval 4200$-$4315A, which contains the strong CH band, against a template in an iterative process. A Bayesian method is applied to calculate the orbit by using the ESPRESSO measurements and others from the literature. A chemical analysis has also been performed in HE0107$-$5240 by means of spectral synthesis. Observations of HE 0107$-$5240 spanning more than 3 years show a monotonic decreasing trend in radial velocity at a rate of approximately by 0.5 ms$^{-1}$d$^{-1}$. The period is constrained at $P_{\rm orb} = 13009_{-1370}^{+1496}$d. New more stringent upper-limits have been found for several elements: a)[Sr/Fe] and [Ba/Fe] are lower than $-0.76$ and $+0.2$ respectively, confirming the star is a CEMP-no; b)$A(Li)< 0.5$ is well below the plateau at $A(Li)=1.1$ found in the Lower Red Giant Branch stars, suggesting Li was originally depleted; and c)the isotopic ratio $^{12}$C/$^{13}$C is 87$\pm6$ showing very low $^{13}$C in contrast to what expected from a spinstar progenitor. We confirm that HE 0107$-$5240 is a binary star with a long period of about 13000d ($\sim36$ years).The carbon isotopic ratio excludes the possibility that the companion has gone through the AGB phase and transferred mass to the currently observed star. The binarity of HE 0107$-$5240 implies some of the first generations of low-mass stars form in multiple systems and indicates that the low metallicity does not preclude the formation of binaries. Finally, a solid indication of $v_{ rad}$ variation has been found also in SMSS 1605$-$1443. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: 15 pages, 7 Figures. Accepted for publication in A&A

Journal ref: A&A 668, A86 (2022)

The probabilistic random forest applied to the QUBRICS survey: improving the selection of high-redshift quasars with synthetic data

Authors: Francesco Guarneri, Giorgio Calderone, Stefano Cristiani, Matteo Porru, Fabio Fontanot, Konstantina Boutsia, Guido Cupani, Andrea Grazian, Valentina D'Odorico, Michael T. Murphy, Angela Bongiorno, Ivano Saccheo, Luciano Nicastro

Abstract: Several recent works have focused on the search for bright, high-z quasars (QSOs) in the South. Among them, the QUasars as BRIght beacons for Cosmology in the Southern hemisphere (QUBRICS) survey has now delivered hundreds of new spectroscopically confirmed QSOs selected by means of machine learning algorithms. Building upon the results obtained by introducing the probabilistic random forest (PRF)… ▽ More Several recent works have focused on the search for bright, high-z quasars (QSOs) in the South. Among them, the QUasars as BRIght beacons for Cosmology in the Southern hemisphere (QUBRICS) survey has now delivered hundreds of new spectroscopically confirmed QSOs selected by means of machine learning algorithms. Building upon the results obtained by introducing the probabilistic random forest (PRF) for the QUBRICS selection, we explore in this work the feasibility of training the algorithm on synthetic data to improve the completeness in the higher redshift bins. We also compare the performances of the algorithm if colours are used as primary features instead of magnitudes. We generate synthetic data based on a composite QSO spectral energy distribution. We first train the PRF to identify QSOs among stars and galaxies, then separate high-z quasar from low-z contaminants. We apply the algorithm on an updated dataset, based on SkyMapper DR3, combined with Gaia eDR3, 2MASS and WISE magnitudes. We find that employing colours as features slightly improves the results with respect to the algorithm trained on magnitude data. Adding synthetic data to the training set provides significantly better results with respect to the PRF trained only on spectroscopically confirmed QSOs. We estimate, on a testing dataset, a completeness of ~86% and a contamination of ~36%. Finally, 207 PRF-selected candidates were observed: 149/207 turned out to be genuine QSOs with z > 2.5, 41 with z < 2.5, 3 galaxies and 14 stars. The result confirms the ability of the PRF to select high-z quasars in large datasets. △ Less

Submitted 30 September, 2022; v1 submitted 15 September, 2022; originally announced September 2022.

Comments: 18 pages, 15 figures. Accepted for publication (MNRAS). Removed a source from the catalogue due to uncertain classification (a star, QID = 1016233)

Automatic model-based telluric correction for the ESPRESSO data reduction software. Model description and application to radial velocity computation

Authors: R. Allart, C. Lovis, J. Faria, X. Dumusque, D. Sosnowska, P. Figueira, A. M. Silva, A. Mehner, F. Pepe, S. Cristiani, R. Rebolo, N. C. Santos, V. Adibekyan, G. Cupani, P. Di Marcantonio, V. D'Odorico, J. I. González Hernández, C. J. A. P. Martins, D. Milaković, N. J. Nunes, A. Sozzetti, A. Suárez Mascareño, H. Tabernero, M. R. Zapatero Osorio

Abstract: Ground-based high-resolution spectrographs are key instruments for several astrophysical domains. Unfortunately, the observed spectra are contaminated by the Earth's atmosphere. While different techniques exist to correct for telluric lines in exoplanet atmospheric studies, in radial velocity (RV) studies, telluric lines with an absorption depth of >2% are generally masked, which poses a problem f… ▽ More Ground-based high-resolution spectrographs are key instruments for several astrophysical domains. Unfortunately, the observed spectra are contaminated by the Earth's atmosphere. While different techniques exist to correct for telluric lines in exoplanet atmospheric studies, in radial velocity (RV) studies, telluric lines with an absorption depth of >2% are generally masked, which poses a problem for faint targets and M dwarfs as most of their RV content is present where telluric contamination is important. We propose a simple telluric model to be embedded in the ESPRESSO DRS. The goal is to provide telluric-free spectra and enable RV measurements, including spectral ranges where telluric lines fall. The model is a line-by-line radiative transfer code that assumes a single atmospheric layer. We use the sky conditions and the physical properties of the lines from HITRAN to create the telluric spectrum. A subset of selected telluric lines is used to robustly fit the spectrum through a Levenberg-Marquardt minimization algorithm. When applied to stellar spectra from A0- to M5-type stars, the residuals of the strongest H2O lines are below 2% for all spectral types, with the exception of M dwarfs, which are within the pseudo-continuum. We then determined the RVs from the telluric-corrected ESPRESSO spectra of Tau Ceti and Proxima. We created telluric-free masks and compared the obtained RVs with the DRS RVs. In the case of Tau Ceti, we identified that micro-telluric lines introduce systematics up to an amplitude of 58 cm/s and with a period of one year. For Proxima, the gain in spectral content at redder wavelengths is equivalent to a gain of 25% in photon noise. This leads to better constraints on the semi-amplitude and eccentricity of Proxima d. We showcase that our model can be applied to other molecules, and thus to other wavelength regions observed by other spectrographs, such as NIRPS. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Comments: 18 pages, 18 figures, accepted to A&A

Journal ref: A&A 666, A196 (2022)

The CUBES Science Case

Authors: Chris Evans, Stefano Cristiani, Cyrielle Opitom, Gabriele Cescutti, Valentina D'Odorico, Juan Manuel Alcalá, Silvia H. P. Alencar, Sergei Balashev, Beatriz Barbuy, Nate Bastian, Umberto Battino, Pamela Cambianica, Roberta Carini, Brad Carter, Santi Cassisi, Bruno Vaz Castilho, Norbert Christlieb, Ryan Cooke, Stefano Covino, Gabriele Cremonese, Katia Cunha, André R. da Silva, Valerio D'Elia, Annalisa De Cia, Gayandhi De Silva , et al. (29 additional authors not shown)

Abstract: We introduce the scientific motivations for the development of the Cassegrain U-Band Efficient Spectrograph (CUBES) that is now in construction for the Very Large Telescope. The assembled cases span a broad range of contemporary topics across Solar System, Galactic and extragalactic astronomy, where observations are limited by the performance of current ground-based spectrographs shortwards of 400… ▽ More We introduce the scientific motivations for the development of the Cassegrain U-Band Efficient Spectrograph (CUBES) that is now in construction for the Very Large Telescope. The assembled cases span a broad range of contemporary topics across Solar System, Galactic and extragalactic astronomy, where observations are limited by the performance of current ground-based spectrographs shortwards of 400nm. A brief background to each case is presented and specific technical requirements on the instrument design that flow-down from each case are identified. These were used as inputs to the CUBES design, that will provide a factor of ten gain in efficiency for astronomical spectroscopy over 300-405nm, at resolving powers of R~24,000 and ~7,000. We include performance estimates that demonstrate the ability of CUBES to observe sources that are up to three magnitudes fainter than currently possible at ground-ultraviolet wavelengths, and we place its predicted performance in the context of existing facillities. △ Less

Submitted 30 September, 2022; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: Accepted for publication in Experimental Astronomy

CUBES, the Cassegrain U-Band Efficient Spectrograph

Authors: S. Cristiani, J. M. Alcalá, S. H. P. Alencar, S. A. Balashev, N. Bastian, B. Barbuy, U. Battino, A. Calcines, G. Calderone, P. Cambianica, R. Carini, B. Carter, S. Cassisi, B. V. Castilho, G. Cescutti, N. Christlieb, R. Cirami, I. Coretti, R. Cooke, S. Covino, G. Cremonese, K. Cunha, G. Cupani, A. R. da Silva, V. De Caprio , et al. (52 additional authors not shown)

Abstract: In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R>20,000 (with a lowe… ▽ More In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at ground-UV wavelengths for the foreseeable future. The Cassegrain U-Band Efficient Spectrograph (CUBES) has been designed to provide high-efficiency (>40%) observations in the near UV (305-400 nm requirement, 300-420 nm goal) at a spectral resolving power of R>20,000 (with a lower-resolution, sky-limited mode of R ~ 7,000). With the design focusing on maximizing the instrument throughput (ensuring a Signal to Noise Ratio (SNR) ~20 per high-resolution element at 313 nm for U ~18.5 mag objects in 1h of observations), it will offer new possibilities in many fields of astrophysics, providing access to key lines of stellar spectra: a tremendous diversity of iron-peak and heavy elements, lighter elements (in particular Beryllium) and light-element molecules (CO, CN, OH), as well as Balmer lines and the Balmer jump (particularly important for young stellar objects). The UV range is also critical in extragalactic studies: the circumgalactic medium of distant galaxies, the contribution of different types of sources to the cosmic UV background, the measurement of H2 and primordial Deuterium in a regime of relatively transparent intergalactic medium, and follow-up of explosive transients. The CUBES project completed a Phase A conceptual design in June 2021 and has now entered the detailed design and construction phase. First science operations are planned for 2028. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: SPIE proceedings, SPIE Astronomical Telescopes + Instrumentation 2022, Montréal, Canada; 20 pages, 13 figures, 2 tables

The black hole and host galaxy growth in an isolated $z\sim 6$ QSO observed with ALMA

Authors: R. Tripodi, C. Feruglio, F. Fiore, M. Bischetti, V. D'Odorico, S. Carniani, S. Cristiani, S. Gallerani, R. Maiolino, A. Marconi, A. Pallottini, E. Piconcelli, L. Vallini, T. Zana

Abstract: The outstanding mass growth of supermassive black holes (SMBHs) at the Reionisation Epoch and how it is related to the concurrent growth of their host galaxies, poses challenges to theoretical models aimed at explaining how these systems formed in short timescales (<1 Gyr). To trace the average evolutionary paths of quasi-stellar objects (QSOs) and their host galaxies in the BH mass-host mass (… ▽ More The outstanding mass growth of supermassive black holes (SMBHs) at the Reionisation Epoch and how it is related to the concurrent growth of their host galaxies, poses challenges to theoretical models aimed at explaining how these systems formed in short timescales (<1 Gyr). To trace the average evolutionary paths of quasi-stellar objects (QSOs) and their host galaxies in the BH mass-host mass ($M_{\rm dyn}$) plane, we compare the star formation rate (SFR), derived from the accurate estimate of the dust temperature and the dust mass ($T_{\rm dust}, M_{\rm dust}$), with the BH accretion rate. To this aim, we analysed a deep, $900$ pc resolution ALMA observation of the sub-mm continuum, [CII] and H$_2$O of the $z\sim 6$ QSO J2310+1855, enabling a detailed study of dust properties and cold gas kinematics. We performed an accurate SED analysis obtaining a dust temperature of $T_{\rm dust} = 71$ K and a dust mass of $M_{\rm dust}= 4.4 \times 10^8\ \rm M_{\odot}$. The implied AGN-corrected SFR is $1240 \ \rm M_{\odot}yr^{-1}$, a factor of 2 smaller than previously reported for this QSO. We derived the best estimate of the dynamical mass $M_{\rm dyn} = 5.2\times 10^{10}\ \rm M_{\odot}$ within $r = 1.7$ kpc, based on a dynamical model of the system. We found that ${\rm SFR}/M_{\rm dyn}>\dot M_{\rm BH}/M_{\rm BH}$, suggesting that AGN feedback might be efficiently acting to slow down the SMBH accretion, while the stellar mass assembly is still vigorously taking place in the host galaxy. In addition, we were also able to detect high-velocity emission on the red and blue sides of the [CII] emission line, that traces a gaseous outflow, and for the first time, we mapped a spatially-resolved water vapour disk through the H$_2$O v=0 $3_{(2,2)}-3_{(1,3)}$ emission line detected at $ν_{\rm obs} = 274.074$ GHz, whose kinematic properties and size are broadly consistent with those of the [CII] disk. △ Less

Submitted 30 July, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: 15 pages, 10 figures, 4 tables. Accepted in A&A

Journal ref: A&A 665, A107 (2022)

Transmission spectroscopy of MASCARA-1b with ESPRESSO: Challenges of overlapping orbital and Doppler tracks

Authors: N. Casasayas-Barris, F. Borsa, E. Pallé, R. Allart, V. Bourrier, J. I. González Hernández, A. Kesseli, A. Sánchez-López, M. R. Zapatero Osorio, I. A. G. Snellen, J. Orell-Miquel, M. Stangret, E. Esparza-Borges, C. Lovis, M. Hooton, M. Lend, A. M. S. Smith, F. Pepe, R. Rebolo, S. Cristiani, N. C. Santos, V. Adibekyan, Y. Alibert, E. Cristo, O. D. S. Demangeon , et al. (10 additional authors not shown)

Abstract: Atmospheric studies at high spectral resolution have shown the presence of molecules, neutral and ionised metals, and hydrogen in the transmission spectrum of ultra-hot Jupiters, and have started to probe the dynamics of their atmospheres. We analyse the transmission spectrum of MASCARA-1b, one of the densest ultra-hot Jupiters orbiting a bright (V=8.3) star. We focus on the CaII H&K, NaI, LiI, H… ▽ More Atmospheric studies at high spectral resolution have shown the presence of molecules, neutral and ionised metals, and hydrogen in the transmission spectrum of ultra-hot Jupiters, and have started to probe the dynamics of their atmospheres. We analyse the transmission spectrum of MASCARA-1b, one of the densest ultra-hot Jupiters orbiting a bright (V=8.3) star. We focus on the CaII H&K, NaI, LiI, H$α$, and KI D1 spectral lines and on the cross-correlated FeI, FeII, CaI, YI, VI, VII, CaH, and TiO lines. For those species that are not present in the stellar spectrum, no detections are reported, but we measure upper limits with an excellent precision ($\sim10$ ppm for particular species). For those species that are present in the stellar spectrum and whose planet-occulted spectral lines induce spurious features in the planetary transmission spectrum, an accurate modelling of the Rossiter-McLaughlin effect (RM) and centre-to-limb variations (CLV) is necessary to recover possible atmospheric signals. In the case of MASCARA-1b, this is difficult due to the overlap between the radial velocities of the stellar surface regions occulted by MASCARA-1b and the orbital track along which the planet atmospheric signal is expected to be found. To try to disentangle a possible planetary signal, we compare our results with models of the RM and CLV effects, and estimate the uncertainties of our models depending on the different system parameters. Unfortunately, more precise measurements of the spin-orbit angle are necessary to better constrain the planet-occulted track and correct for the transit effects in the transmission spectrum with enough precision to be able to detect or discard possible planetary absorptions. Finally, we discuss the possibility that non-detections are related to the low absorption expected for a high surface gravity planet such as MASCARA-1b. △ Less

Submitted 19 June, 2022; originally announced June 2022.

Comments: Accepted

Fundamental physics with ESPRESSO: Constraints on Bekenstein and dark energy models from astrophysical and local probes

Authors: C. J. A. P. Martins, S. Cristiani, G. Cupani, V. D'Odorico, R. Génova Santos, A. C. O. Leite, C. M. J. Marques, D. Milaković, P. Molaro, Michael T. Murphy, N. J. Nunes, Tobias M. Schmidt, V. Adibekyan, Y. Alibert, Paolo Di Marcantonio, J. I. González Hernández, D. Mégevand, E. Palle, F. A. Pepe, N. C. Santos, S. G. Sousa, A. Sozzetti, A. Suárez Mascareño, M. R. Zapatero Osorio

Abstract: Dynamical scalar fields in an effective four-dimensional field theory are naturally expected to couple to the rest of the theory's degrees of freedom, unless some new symmetry is postulated to suppress these couplings. In particular, a coupling to the electromagnetic sector will lead to spacetime variations of the fine-structure constant, $α$. Astrophysical tests of the space-time stability of… ▽ More Dynamical scalar fields in an effective four-dimensional field theory are naturally expected to couple to the rest of the theory's degrees of freedom, unless some new symmetry is postulated to suppress these couplings. In particular, a coupling to the electromagnetic sector will lead to spacetime variations of the fine-structure constant, $α$. Astrophysical tests of the space-time stability of $α$ are therefore a powerful probe of new physics. Here we use ESPRESSO and other contemporary measurements of $α$, together with background cosmology data, local laboratory atomic clock and Weak Equivalence Principle measurements, to place stringent constraints on the simplest examples of the two broad classes of varying $α$ models: Bekenstein models and quintessence-type dark energy models, both of which are parametric extensions of the canonical $Λ$CDM model. In both cases, previously reported constraints are improved by more than a factor of ten. This improvement is largely due to the very strong local constraints, but astrophysical measurements can help to break degeneracies between cosmology and fundamental physics parameters. △ Less

Submitted 27 May, 2022; originally announced May 2022.

Comments: 8 pages, 2 figures, 2 tables; in press at Phys. Rev. D

A novel framework for semi-Bayesian radial velocities through template matching

Authors: A. M. Silva, J. P. Faria, N. C. Santos, S. G. Sousa, P. T. P. Viana, J. H. C. Martins, P. Figueira, C. Lovis, F. Pepe, S. Cristiani, R. Rebolo, R. Allart, A. Cabral, A. Mehner, A. Sozzetti, A. Suárez Mascareño, C. J. A. P. Martins, D. Ehrenreich, D. Mégevand, E. Palle, G. Lo Curto, H. M. Tabernero, J. Lillo-Box, J. I. González Hernández, M. R. Zapatero Osorio , et al. (6 additional authors not shown)

Abstract: The detection and characterization of an increasing variety of exoplanets has been in part possible thanks to the continuous development of high-resolution, stable spectrographs, and using the Doppler radial-velocity (RV) method. The Cross Correlation Function (CCF) method is one of the traditional approaches for RV extraction. More recently, template matching was introduced as an advantageous alt… ▽ More The detection and characterization of an increasing variety of exoplanets has been in part possible thanks to the continuous development of high-resolution, stable spectrographs, and using the Doppler radial-velocity (RV) method. The Cross Correlation Function (CCF) method is one of the traditional approaches for RV extraction. More recently, template matching was introduced as an advantageous alternative for M-dwarf stars. In this paper, we describe a new implementation of template matching within a semi-Bayesian framework, providing a more statistically principled characterization of the RV measurements. In this context, a common RV shift is used to describe the difference between each spectral order of a given stellar spectrum and a template built from the available observations. Posterior probability distributions are obtained for the relative RV associated with each spectrum, after marginalizing with respect to the continuum. This methodology was named S-BART: Semi-Bayesian Approach for RVs with Template-matching, and it can be applied to HARPS and ESPRESSO. The application of our method to HARPS archival observations of Barnard's star allowed us to validate our implementation against HARPS-TERRA and SERVAL. Then, we applied it to 33 ESPRESSO targets, evaluating its performance and comparing it with the CCF method. We found a decrease in the median RV scatter of \sim 10\% and \sim 4\% for M- and K-type stars, respectively. S-BART yields more precise RV estimates than the CCF method, particularly in the case of M-type stars where a median uncertainty of \sim 15 cm/s is achieved over 309 observations. Further, we estimated the nightly zero point (NZP) of ESPRESSO, finding a weighted NZP scatter below \sim 0.7 m/s. As this includes stellar variability, photon noise, and potential planetary signals, it should be taken as an upper limit of the RV precision attainable with ESPRESSO data. △ Less

Submitted 29 April, 2022; originally announced May 2022.

Comments: 17 pages, 15 Figures, Accepted for publication by Astronomy & Astrophysics (A&A), Code available in https://github.com/iastro-pt/sBART

Journal ref: A&A 663, A143 (2022)

Widespread, strong outflows in XQR-30 quasars at the Reionisation epoch

Authors: M. Bischetti, C. Feruglio, V. D'Odorico, N. Arav, E. Bañados, G. Becker, S. E. I. Bosman, S. Carniani, S. Cristiani, G. Cupani, R. Davies, A. C. Eilers, E. P. Farina, A. Ferrara, R. Maiolino, C. Mazzucchelli, A. Mesinger, R. Meyer, M. Onoue, E. Piconcelli, E. Ryan-Weber, J-T. Schindler, F. Wang, J. Yang, Y. Zhu , et al. (1 additional authors not shown)

Abstract: Luminous quasars powered by accretion onto billion solar mass black holes already exist at the epoch of Reionisation, when the Universe was 0.5-1 Gyr old. These objects likely reside in over-dense regions of the Universe, and will grow to form today's giant galaxies. How their huge black holes formed in such short times is debated, particularly as they lie above the local black hole mass-galaxy dy… ▽ More Luminous quasars powered by accretion onto billion solar mass black holes already exist at the epoch of Reionisation, when the Universe was 0.5-1 Gyr old. These objects likely reside in over-dense regions of the Universe, and will grow to form today's giant galaxies. How their huge black holes formed in such short times is debated, particularly as they lie above the local black hole mass-galaxy dynamical mass correlation, thus following the black hole-dominance growth path. It is unknown what slowed down the black hole growth, leading towards the symbiotic growth observed in the local Universe, and when this process started, although black hole feedback is a likely driver. This deadlock is due to the lack of large, homogeneous samples of high-redshift quasars with high-quality, broad-band spectroscopic information. Here we report results from a VLT/X-shooter survey of 30 quasars at redshift 5.8$\le$z$\le$6.6 (XQR-30). About 50% of their spectra reveal broad blue-shifted absorption line (BAL) throughs, tracing powerful ionised winds. The BAL fraction in z$\gtrsim$6 quasars is 2-3 times higher than in quasars at z~2-4.5. XQR-30 BAL quasars exhibit extreme outflow velocities, up to 17% of the light speed, rarely observed at lower redshift. These outflows inject large amounts of energy into the galaxy interstellar medium, which can contrast nuclear gas accretion, slowing down the black-hole growth. The star-formation rate in high-z quasar hosts is generally $>$100 M$_\odot$/yr, so these galaxies are growing at a fast rate. The BAL phase may then mark the beginning of significant feedback, acting first on black hole growth and possibly later on galaxy growth. The red optical colors of BAL quasars at z$\gtrsim$6 indeed suggest that these systems are dusty and may be caught during an initial quenching phase of obscured accretion. △ Less

Submitted 17 May, 2022; v1 submitted 29 April, 2022; originally announced May 2022.

Comments: This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature, and is available online at https://doi.org/10.1038/s41586-022-04608-1

Journal ref: Nature, 605, 244-247 (2022)

Fundamental Physics with ESPRESSO, Constraining a simple parametrisation for varying $α$

Authors: Vitor da Fonseca, Tiago Barreiro, Nelson J. Nunes, Stefano Cristiani, Guido Cupani, Valentina D'Odorico, Ricardo Génova Santos, Ana C. O. Leite, Catarina M. J. Marques, Carlos J. A. P. Martins, Dinko Milaković, Paolo Molaro, Michael T. Murphy, Tobias M. Schmidt, Manuel Abreu, Vardan Adibekyan, Alexandre Cabral, Paolo Di Marcantonio, Jonay I. González Hernández, Enric Palle, Francesco A. Pepe, Rafael Rebolo, Nuno C. Santos, Sérgio G. Sousa, Alessandro Sozzetti , et al. (2 additional authors not shown)

Abstract: The spectrograph ESPRESSO recently obtained a limit on the variation of the fine-structure constant, $α$, through measurements along the line of sight of a bright quasar with a precision of $1.36$ ppm at $1σ$ level. This imposes new constraints on cosmological models with a varying $α$. We assume such a model where the electromagnetic sector is coupled to a scalar field dark energy responsible for… ▽ More The spectrograph ESPRESSO recently obtained a limit on the variation of the fine-structure constant, $α$, through measurements along the line of sight of a bright quasar with a precision of $1.36$ ppm at $1σ$ level. This imposes new constraints on cosmological models with a varying $α$. We assume such a model where the electromagnetic sector is coupled to a scalar field dark energy responsible for the current acceleration of the Universe. We parametrise the variation of $α$ with two extra parameters, one defining the cosmological evolution of the quintessence component and the other fixing the coupling with the electromagnetic field. The objective of this work is to constrain these parameters with both astrophysical and local probes. We also carried out a comparative analysis of how each data probe may constrain our parametrisation. We performed a Bayesian analysis by comparing the predictions of the model with observations. The astrophysical datasets are composed of quasar spectra measurements, including the latest ESPRESSO data point, as well as Planck observations of the cosmic microwave background. We combined these with local results from atomic clocks and the MICROSCOPE experiment. The constraints placed on the quintessence parameter are consistent with a null variation of the field, and are therefore compatible with a $Λ$CDM cosmology. The constraints on the coupling to the electromagnetic sector are dominated by the Eötvös parameter local bound. More precise measurements with ESPRESSO will be extremely important to study the cosmological evolution of $α$ as it probes an interval of redshift not accessible to other types of observations. However, for this particular model, current available data favour a null variation of $α$ resulting mostly from the strong MICROSCOPE limits. △ Less

Submitted 14 November, 2022; v1 submitted 6 April, 2022; originally announced April 2022.

Comments: 9 pages, 7 figures, as published in A&A

Journal ref: A&A 666, A57 (2022)

The Low-Redshift Lyman Continuum Survey II: New Insights into LyC Diagnostics

Authors: Sophia R. Flury, Anne E. Jaskot, Harry C. Ferguson, Gabor Worseck, Kirill Makan, John Chisholm, Alberto Saldana-Lopez, Daniel Schaerer, Stephan McCandliss, Bingjie Wang, N. M. Ford, M. S. Oey, Timothy Heckman, Zhiyuan Ji, Mauro Giavalisco, Ricardo Amorin, Hakim Atek, Jeremy Blaizot, Sanchayeeta Borthakur, Cody Carr, Marco Castellano, Stefano Cristiani, Stephane de Barros, Mark Dickinson, Steven L. Finkelstein , et al. (22 additional authors not shown)

Abstract: The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z {\gtrsim} 6) due to intergalactic H I absorption. To identify Lyman continuum emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical t… ▽ More The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z {\gtrsim} 6) due to intergalactic H I absorption. To identify Lyman continuum emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical test of these diagnostics. While optical depth indicators based on Lyα, such as peak velocity separation and equivalent width, perform well, we also find that other diagnostics, such as the [O III]/[O II] flux ratio and star formation rate surface density, predict whether a galaxy is a LCE. The relationship between these galaxy properties and the fraction of escaping LyC flux suggests that LyC escape depends strongly on H I column density, ionization parameter, and stellar feedback. We find LCEs occupy a range of stellar masses, metallicities, star formation histories, and ionization parameters, which may indicate episodic and/or different physical causes of LyC escape. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: ApJ, accepted. 31 pages, 26 figures, 1 table

The CUBES Instrument Model and Simulation Tools. Their role in the project Phase A study

Authors: Matteo Genoni, Marco Landoni, Guido Cupani, Mariagrazia Franchini, Roberto Cirami, Alessio Zanutta, Chris Evans, Paolo Di Marcantonio, Stefano Cristiani, Andrea Trost, Sonia Zorba

Abstract: We present the simulation tools developed to aid the design phase of the Cassegrain U-Band Efficient Spectrograph (CUBES) for the Very Large Telescope (VLT), exploring aspects of the system design and evaluating the performance for different design configurations. CUBES aims to be the 'ultimate' ultraviolet (UV) instrument at the European Southern Observatory (ESO) in terms of throughput, with the… ▽ More We present the simulation tools developed to aid the design phase of the Cassegrain U-Band Efficient Spectrograph (CUBES) for the Very Large Telescope (VLT), exploring aspects of the system design and evaluating the performance for different design configurations. CUBES aims to be the 'ultimate' ultraviolet (UV) instrument at the European Southern Observatory (ESO) in terms of throughput, with the goal to cover the bluest part of the spectrum accessible from the ground (300 nm to 400 nm) with the highest possible efficiency. Here we introduce the End-to-End (E2E) and the Exposure Time Calculator (ETC) tools. The E2E simulator has been developed with different versions to meet the needs of different users, including a version that can be accessed for use by the broader scientific community using a Jupyter notebook. The E2E tool was used by the system team to help define the Phase A baseline design of the instrument, as well as in scientific evaluation of a possible low-resolution mode. The ETC is a web-based tool through which the science community are able to test a range of science cases for CUBES, demonstrating its potential to push the limiting magnitude for the detection of specific UV-features, such as abundance estimates of beryllium in main sequence stars. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: Accepted for publication in Experimental Astronomy

CUBES Phase A design overview -- The Cassegrain U-Band Efficient Spectrograph for the Very Large Telescope

Authors: Alessio Zanutta, Stefano Cristiani, David Atkinson, Veronica Baldini, Andrea Balestra, Beatriz Barbuy, Vanessa Bawden P. Macanhan, Ariadna Calcines, Giorgio Calderone, Scott Case, Bruno V. Castilho, Gabriele Cescutti, Roberto Cirami, Igor Coretti, Stefano Covino, Guido Cupani, Vincenzo De Caprio, Hans Dekker, Paolo Di Marcantonio, Valentina D'Odorico, Heitor Ernandes, Chris Evans, Tobias Feger, Carmen Feiz, Mariagrazia Franchini , et al. (29 additional authors not shown)

Abstract: We present the baseline conceptual design of the Cassegrain U-Band Efficient Spectrograph (CUBES) for the Very Large Telescope. CUBES will provide unprecedented sensitivity for spectroscopy on a 8 - 10 m class telescope in the ground ultraviolet (UV), spanning a bandwidth of > 100 nm that starts at 300 nm, the shortest wavelength accessible from the ground. The design has been optimized for end-to… ▽ More We present the baseline conceptual design of the Cassegrain U-Band Efficient Spectrograph (CUBES) for the Very Large Telescope. CUBES will provide unprecedented sensitivity for spectroscopy on a 8 - 10 m class telescope in the ground ultraviolet (UV), spanning a bandwidth of > 100 nm that starts at 300 nm, the shortest wavelength accessible from the ground. The design has been optimized for end-to-end efficiency and provides a spectral resolving power of R > 20000, that will unlock a broad range of new topics across solar system, Galactic and extraglactic astronomy. The design also features a second, lower-resolution (R \sim 7000) mode and has the option of a fiberlink to the UVES instrument for simultaneous observations at longer wavelengths. Here we present the optical, mechanical and software design of the various subsystems of the instrument after the Phase A study of the project. We discuss the expected performances for the layout choices and highlight some of the performance trade-offs considered to best meet the instrument top-level requirements. We also introduce the model-based system engineering approach used to organize and manage the project activities and interfaces, in the context that it is increasingly necessary to integrate such tools in the development of complex astronomical projects. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: Accepted for publication in Experimental Astronomy

The evolution of the Si IV content in the Universe from the epoch of reionization to cosmic noon

Authors: V. D'Odorico, K. Finlator, S. Cristiani, G. Cupani, S. Perrotta, F. Calura, M. Cènturion, G. Becker, T. A. M. Berg, S. Lopez, S. Ellison, E. Pomante

Abstract: We investigate the abundance and distribution of metals in the high-redshift intergalactic medium and circum-galactic medium through the analysis of a sample of almost 600 SiIV absorption lines detected in high and intermediate resolution spectra of 147 quasars. The evolution of the number density of SiIV lines, the column density distribution function and the cosmic mass density are studied in th… ▽ More We investigate the abundance and distribution of metals in the high-redshift intergalactic medium and circum-galactic medium through the analysis of a sample of almost 600 SiIV absorption lines detected in high and intermediate resolution spectra of 147 quasars. The evolution of the number density of SiIV lines, the column density distribution function and the cosmic mass density are studied in the redshift interval 1.7 <= z <= 6.2 and for log N(SiIV) >= 12.5. All quantities show a rapid increase between z~6 and z< 5 and then an almost constant behaviour to z~2 in very good agreement with what is already observed for CIV absorption lines. The present results are challenging for numerical simulations: when simulations reproduce our SiIV results, they tend to underpredict the properties of CIV, and when the properties of CIV are reproduced, the number of strong SiIV lines (log N(SiIV) > 14) is overpredicted. △ Less

Submitted 24 February, 2022; originally announced February 2022.

Comments: 14 pages, 8 figures, accepted for publication in MNRAS

A candidate short-period sub-Earth orbiting Proxima Centauri

Authors: J. P. Faria, A. Suárez Mascareño, P. Figueira, A. M. Silva, M. Damasso, O. Demangeon, F. Pepe, N. C. Santos, R. Rebolo, S. Cristiani, V. Adibekyan, Y. Alibert, R. Allart, S. C. C. Barros, A. Cabral, V. D'Odorico, P. Di Marcantonio, X. Dumusque, D. Ehrenreich, J. I. González Hernández, N. Hara, J. Lillo-Box, G. Lo Curto, C. Lovis, C. J. A. P. Martins , et al. (12 additional authors not shown)

Abstract: Proxima Centauri is the closest star to the Sun. This small, low-mass, mid M dwarf is known to host an Earth-mass exoplanet with an orbital period of 11.2 days within the habitable zone, as well as a long-period planet candidate with an orbital period of close to 5 years. We report on the analysis of a large set of observations taken with the ESPRESSO spectrograph at the VLT aimed at a thorough ev… ▽ More Proxima Centauri is the closest star to the Sun. This small, low-mass, mid M dwarf is known to host an Earth-mass exoplanet with an orbital period of 11.2 days within the habitable zone, as well as a long-period planet candidate with an orbital period of close to 5 years. We report on the analysis of a large set of observations taken with the ESPRESSO spectrograph at the VLT aimed at a thorough evaluation of the presence of a third low-mass planetary companion, which started emerging during a previous campaign. Radial velocities (RVs) were calculated using both a cross-correlation function (CCF) and a template matching approach. The RV analysis includes a component to model Proxima's activity using a Gaussian process (GP). We use the CCF's full width at half maximum to help constrain the GP, and we study other simultaneous observables as activity indicators in order to assess the nature of any potential RV signals. We detect a signal at 5.12 $\pm$ 0.04 days with a semi-amplitude of 39 $\pm$ 7 cm/s. The analysis of subsets of the ESPRESSO data, the activity indicators, and chromatic RVs suggest that this signal is not caused by stellar variability but instead by a planetary companion with a minimum mass of 0.26 $\pm$ 0.05 $M_\oplus$ (about twice the mass of Mars) orbiting at 0.029 au from the star. The orbital eccentricity is well constrained and compatible with a circular orbit. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: 17 pages, 13 figures

The Low-Redshift Lyman Continuum Survey I: New, Diverse Local Lyman-Continuum Emitters

Authors: Sophia R. Flury, Anne E. Jaskot, Harry C. Ferguson, Gabor Worseck, Kirill Makan, John Chisholm, Alberto Saldana-Lopez, Daniel Schaerer, Stephan McCandless, Bingjie Wang, N. M. Ford, Timothy Heckman, Zhiyuan Ji, Mauro Giavalisco, Ricardo Amorin, Hakim Atek, Jeremy Blaizot, Sanchayeeta Borthakur, Cody Carr, Marco Castellano, Stefano Cristiani, Stephane de Barros, Mark Dickinson, Steven L. Finkelstein, Brian Fleming , et al. (22 additional authors not shown)

Abstract: The origins of Lyman continuum (LyC) photons responsible for the reionization of the universe are as of yet unknown and highly contested. Detecting LyC photons from the epoch of reionization is not possible due to absorption by the intergalactic medium, which has prompted the development of several indirect diagnostics to infer the rate at which galaxies contribute LyC photons to reionize the univ… ▽ More The origins of Lyman continuum (LyC) photons responsible for the reionization of the universe are as of yet unknown and highly contested. Detecting LyC photons from the epoch of reionization is not possible due to absorption by the intergalactic medium, which has prompted the development of several indirect diagnostics to infer the rate at which galaxies contribute LyC photons to reionize the universe by studying lower-redshift analogs. We present the Low-redshift Lyman Continuum Survey (LzLCS) comprising measurements made with HST/COS for a z=0.2-0.4 sample of 66 galaxies. After careful processing of the FUV spectra, we obtain a total of 35 Lyman continuum emitters (LCEs) detected with 97.725% confidence, nearly tripling the number of known local LCEs. We estimate escape fractions from the detected LyC flux and upper limits on the undetected LyC flux, finding a range of LyC escape fractions up to 50%. Of the 35 LzLCS LCEs, 12 have LyC escape fractions greater than 5%, more than doubling the number of known local LCEs with cosmologically relevant LyC escape. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: ApJS, accepted. 28 pages, 19 figures, 8 tables. Full machine readable tables will be made available by the publisher at the time of publication

CaRM: Exploring the chromatic Rossiter-McLaughlin effect. The cases of HD 189733b and WASP-127b

Authors: E. Cristo, N. C. Santos, O. Demangeon, J. H. C. Martins, P. Figueira, N. Casasayas-Barris, M. R. Zapatero Osorio, F. Borsa, S. G. Sousa, M. Oshagh, G. Micela, H. M. Tabernero, J. V. Seidel, S. Cristiani, F. Pepe, R. Rebolo, V. Adibekyan, R. Allart, Y. Alibert, T. Azevedo Silva, V. Bourrier, A. Cabral, E. Esparza Borges, J. I. González Hernández, J. Lillo-Box , et al. (12 additional authors not shown)

Abstract: In this paper we introduce CaRM, a semi-automatic code for the retrieval of broadband transmission spectra of transiting planets through the chromatic Rossiter-McLaughlin method. We applied it to HARPS and ESPRESSO observations of two exoplanets to retrieve the transmission spectrum and we analyze its fitting transmission models. We used the strong radius dependence of the Rossiter-McLaughlin (RM)… ▽ More In this paper we introduce CaRM, a semi-automatic code for the retrieval of broadband transmission spectra of transiting planets through the chromatic Rossiter-McLaughlin method. We applied it to HARPS and ESPRESSO observations of two exoplanets to retrieve the transmission spectrum and we analyze its fitting transmission models. We used the strong radius dependence of the Rossiter-McLaughlin (RM) effect amplitude, caused by planetary companions, to measure the apparent radius change caused by the exoplanet atmosphere. In order to retrieve the transmission spectrum, the radial velocities, which were computed over wavelength bins that encompass several spectral orders, were used to simultaneously fit the Keplerian motion and the RM effect. From this, the radius ratio was computed as a function of the wavelength, which allows one to retrieve the low-resolution broadband transmission spectrum of a given exoplanet. CaRM offers the possibility to use two Rossiter-McLaughlin models taken from ARoME and PyAstronomy, associated with a Keplerian function to fit radial velocities during transit observations automatically. Furthermore it offers the possibility to use some methods that could, in theory, mitigate the effect of perturbation in the radial velocities during transits. The CaRM code allows one to retrieve the transmission spectrum of a given exoplanet using minimal user interaction. We demonstrate that it allows one to compute the low-resolution broadband transmission spectra of exoplanets observed using high-resolution spectrographs such as HARPS and ESPRESSO. △ Less

Submitted 19 January, 2022; v1 submitted 17 January, 2022; originally announced January 2022.

Journal ref: A&A 660, A52 (2022)

Fundamental physics with ESPRESSO: Precise limit on variations in the fine-structure constant towards the bright quasar HE 0515$-$4414

Authors: Michael T. Murphy, Paolo Molaro, Ana C. O. Leite, Guido Cupani, Stefano Cristiani, Valentina D'Odorico, Ricardo Génova Santos, Carlos J. A. P. Martins, Dinko Milaković, Nelson J. Nunes, Tobias M. Schmidt, Francesco A. Pepe, Rafael Rebolo, Nuno C. Santos, Sérgio G. Sousa, Maria-Rosa Zapatero Osorio, Manuel Amate, Vardan Adibekyan, Yann Alibert, Carlos Allende Prieto, Veronica Baldini, Willy Benz, François Bouchy, Alexandre Cabral, Hans Dekker , et al. (18 additional authors not shown)

Abstract: The strong intervening absorption system at redshift 1.15 towards the very bright quasar HE 0515$-$4414 is the most studied absorber for measuring possible cosmological variations in the fine-structure constant, $α$. We observed HE 0515$-$4414 for 16.1$\,$h with the Very Large Telescope and present here the first constraint on relative variations in $α$ with parts-per-million (ppm) precision from… ▽ More The strong intervening absorption system at redshift 1.15 towards the very bright quasar HE 0515$-$4414 is the most studied absorber for measuring possible cosmological variations in the fine-structure constant, $α$. We observed HE 0515$-$4414 for 16.1$\,$h with the Very Large Telescope and present here the first constraint on relative variations in $α$ with parts-per-million (ppm) precision from the new ESPRESSO spectrograph: $Δα/α= 1.3 \pm 1.3_{\rm stat} \pm 0.4_{\rm sys}\,{\rm ppm}$. The statistical uncertainty (1$σ$) is similar to the ensemble precision of previous large samples of absorbers, and derives from the high S/N achieved ($\approx$105 per 0.4$\,$km$\,$s$^{-1}$ pixel). ESPRESSO's design, and calibration of our observations with its laser frequency comb, effectively removed wavelength calibration errors from our measurement. The high resolving power of our ESPRESSO spectrum ($R=145000$) enabled the identification of very narrow components within the absorption profile, allowing a more robust analysis of $Δα/α$. The evidence for the narrow components is corroborated by their correspondence with previously detected molecular hydrogen and neutral carbon. The main remaining systematic errors arise from ambiguities in the absorption profile modelling, effects from redispersing the individual quasar exposures, and convergence of the parameter estimation algorithm. All analyses of the spectrum, including systematic error estimates, were initially blinded to avoid human biases. We make our reduced ESPRESSO spectrum of HE 0515$-$4414 publicly available for further analysis. Combining our ESPRESSO result with 28 measurements, from other spectrographs, in which wavelength calibration errors have been mitigated, yields a weighted mean $Δα/α= -0.5 \pm 0.5_{\rm stat} \pm 0.4_{\rm sys}\,$ppm at redshifts 0.6-2.4. △ Less

Submitted 10 December, 2021; originally announced December 2021.

Comments: Accepted by A&A. The ESPRESSO quasar spectrum, reduced contributing exposures, and absorption profile fits are available at https://github.com/MTMurphy77/ESPRESSO_HE0515-4414

Journal ref: A&A 658, A123 (2022)

Near-infrared spectroscopy of extreme BAL QSOs from the QUBRICS bright quasar survey

Authors: Guido Cupani, Giorgio Calderone, Pierluigi Selvelli, Stefano Cristiani, Konstantina Boutsia, Andrea Grazian, Fabio Fontanot, Francesco Guarneri, Valentina D'Odorico, Emanuele Giallongo, Nicola Menci

Abstract: We report on the spectral confirmation of 18 QSO candidates from the "QUasars as BRIght beacons for Cosmology in the Southern hemisphere'' survey (QUBRICS), previously observed in the optical band, for which we acquired new spectroscopic data in the near-infrared band with the Folded-port InfraRed Echellette spectrograph (FIRE) at the Magellan Baade telescope. In most cases, further observations w… ▽ More We report on the spectral confirmation of 18 QSO candidates from the "QUasars as BRIght beacons for Cosmology in the Southern hemisphere'' survey (QUBRICS), previously observed in the optical band, for which we acquired new spectroscopic data in the near-infrared band with the Folded-port InfraRed Echellette spectrograph (FIRE) at the Magellan Baade telescope. In most cases, further observations were prompted by the peculiar nature of the targets, whose optical spectra displayed unexpected absorption features. All candidates have been confirmed as bona fide QSOs, with average emission redshift $z\simeq 2.1$. The analysis of the emission and absorption features in the spectra, performed with Astrocook and QSFit, reveals that the large majority of these objects are broad-absorption line (BAL) QSOs, with almost half of them displaying strong Fe II absorption (typical of the so-called FeLoBAL QSOs). The detection of such a large fraction of rare objects (which are estimated to account for less than one percent of the general QSO population) is interpreted as an unexpected (yet favourable) consequence of the particular candidate selection procedure adopted within the QUBRICS survey. The measured properties of FeLoBAL QSOs observed so far provide no evidence that they are a manifestation of a particular stage in AGN evolution. In this paper we present an explorative analysis of the individual QSOs, to serve as a basis for a further, more detailed investigation. △ Less

Submitted 5 December, 2021; originally announced December 2021.

Comments: 20 pages, 7 figures. This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review. Associated data and procedures available at https://www.ict.inaf.it/index.php/31-doi/137-ds-2021-03

The Rossiter-McLaughlin effect Revolutions: An ultra-short period planet and a warm mini-Neptune on perpendicular orbits

Authors: V. Bourrier, C. Lovis, M. Cretignier, R. Allart, X. Dumusque, J. -B. Delisle, A. Deline, S. G. Sousa, V. Adibekyan, Y. Alibert, S. C. C. Barros, F. Borsa, S. Cristiani, O. Demangeon, D. Ehrenreich, P. Figueira, J. I. González Hernández, M. Lendl, J. Lillo-Box, G. Lo Curto, P. Di Marcantonio, C. J. A. P. Martins, D. Mégevand, A. Mehner, G. Micela , et al. (12 additional authors not shown)

Abstract: Comparisons of the alignment of exoplanets with a common host star can be used to distinguish among concurrent evolution scenarios. However, multi-planet systems usually host mini-Neptunes and super-Earths, whose size make orbital architecture measurements challenging. We introduce the Rossiter-McLaughlin effect Revolutions technique, which can access spin-orbit angles of small planets by exploiti… ▽ More Comparisons of the alignment of exoplanets with a common host star can be used to distinguish among concurrent evolution scenarios. However, multi-planet systems usually host mini-Neptunes and super-Earths, whose size make orbital architecture measurements challenging. We introduce the Rossiter-McLaughlin effect Revolutions technique, which can access spin-orbit angles of small planets by exploiting the full information contained in spectral transit time series. We validated the technique on published HARPS-N data of the mini-Neptune HD3167c, refining its high sky-projected spin-orbit angle (-108.9+5.4-5.5 deg), and we applied it to new ESPRESSO observations of the super-Earth HD3167b, revealing an aligned orbit (-6.6+6.6-7.9 deg). Surprisingly different variations in the contrast of the stellar lines occulted by the planets can be reconciled with a latitudinal dependence of the stellar line shape. In this scenario, a joint fit to both datasets constrains the inclination of the star (111.6+3.1-3.3 deg) and the 3D spin-orbit angles of HD3167b (29.5+7.2-9.4 deg) and HD3167c (107.7+5.1-4.9 deg). The projected spin-orbit angles do not depend on the model for the line contrast variations, and so, with a mutual inclination of 102.3+7.4-8.0 deg, we conclude that the two planets are on perpendicular orbits. This could be explained by HD3167b being strongly coupled to the star and retaining its primordial alignment, whereas HD3167c would have been brought to a nearly polar orbit via secular gravitational interactions with an outer companion. Follow-up observations and dynamical evolution simulations are required to search for this companion and explore this scenario. HD3167b is the smallest exoplanet with a confirmed spectroscopic Rossiter-McLaughlin signal. Our new technique opens the way to determining the orbital architectures of the super-Earth and Earth-sized planet populations. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 20 pages, 19 figures, accepted for publication in A&A (14 September 2021)

The space density of ultra-luminous QSOs at the end of reionization epoch by the QUBRICS Survey and the AGN contribution to the hydrogen ionizing background

Authors: A. Grazian, E. Giallongo, K. Boutsia, G. Calderone, S. Cristiani, G. Cupani, F. Fontanot, F. Guarneri, Y. Ozdalkiran

Abstract: Motivated by evidences favoring a rapid and late hydrogen reionization process completing at z~5.2-5.5 and mainly driven by rare and luminous sources, we have reassessed the estimate of the space density of ultra-luminous QSOs at z~5 in the framework of the QUBRICS survey. A ~90% complete sample of 14 spectroscopically confirmed QSOs at M1450<-28.3 and 4.5<z<5.0 has been derived in an area of 12,4… ▽ More Motivated by evidences favoring a rapid and late hydrogen reionization process completing at z~5.2-5.5 and mainly driven by rare and luminous sources, we have reassessed the estimate of the space density of ultra-luminous QSOs at z~5 in the framework of the QUBRICS survey. A ~90% complete sample of 14 spectroscopically confirmed QSOs at M1450<-28.3 and 4.5<z<5.0 has been derived in an area of 12,400 sq. deg., thanks to multi-wavelength selection and GAIA astrometry. The space density of z~5 QSOs within -29.3<M1450<-28.3 is three times higher than previous determinations. Our results suggest a steep bright-end slope for the QSO luminosity function at z~5 and a mild redshift evolution of the space density of ultra-bright QSOs (M1450~-28.5) at 3<z<5.5, in agreement with the redshift evolution of the much fainter AGN population at M1450~-23. These findings are consistent with a pure density evolution for the AGN population at z>3. Adopting our z~4 QSO luminosity function and applying a mild density evolution in redshift, a photo-ionization rate of GammaHI=0.46+0.17-0.09x10^-12 s^-1 has been obtained at z=4.75, assuming an escape fraction of ~70% and a steep faint-end slope of the AGN luminosity function. The derived photo-ionization rate is ~50-100% of the ionizing background measured at the end of the reionization epoch, suggesting that AGNs could play an important role in the cosmological reionization process. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJ. 15 pages, 3 figures, 3 tables