Showing 1–50 of 393 results for author: Davies, B

Characterisation of the Warm-Jupiter TOI-1130 system with CHEOPS and photo-dynamical approach

Authors: L. Borsato, D. Degen, A. Leleu, M. J. Hooton, J. A. Egger, A. Bekkelien, A. Brandeker, A. Collier Cameron, M. N. Günther, V. Nascimbeni, C. M. Persson, A. Bonfanti, T. G. Wilson, A. C. M. Correia, T. Zingales, T. Guillot, A. H. M. J. Triaud, G. Piotto, D. Gandolfi, L. Abe, Y. Alibert, R. Alonso, T. Bárczy, D. Barrado Navascues, S. C. C. Barros , et al. (71 additional authors not shown)

Abstract: Among the thousands of exoplanets discovered to date, approximately a few hundred gas giants on short-period orbits are classified as "lonely" and only a few are in a multi-planet system with a smaller companion on a close orbit. The processes that formed multi-planet systems hosting gas giants on close orbits are poorly understood, and only a few examples of this kind of system have been observed… ▽ More Among the thousands of exoplanets discovered to date, approximately a few hundred gas giants on short-period orbits are classified as "lonely" and only a few are in a multi-planet system with a smaller companion on a close orbit. The processes that formed multi-planet systems hosting gas giants on close orbits are poorly understood, and only a few examples of this kind of system have been observed and well characterised. Within the contest of multi-planet system hosting gas-giant on short orbits, we characterise TOI-1130 system by measuring masses and orbital parameters. This is a 2-transiting planet system with a Jupiter-like planet (c) on a 8.35 days orbit and a Neptune-like planet (b) on an inner (4.07 days) orbit. Both planets show strong anti-correlated transit timing variations (TTVs). Furthermore, radial velocity (RV) analysis showed an additional linear trend, a possible hint of a non-transiting candidate planet on a far outer orbit. Since 2019, extensive transit and radial velocity observations of the TOI-1130 have been acquired using TESS and various ground-based facilities. We present a new photo-dynamical analysis of all available transit and RV data, with the addition of new CHEOPS and ASTEP+ data that achieve the best precision to date on the planetary radii and masses and on the timings of each transit. We were able to model interior structure of planet b constraining the presence of a gaseous envelope of H/He, while it was not possible to assess the possible water content. Furthermore, we analysed the resonant state of the two transiting planets, and we found that they lie just outside the resonant region. This could be the result of the tidal evolution that the system underwent. We obtained both masses of the planets with a precision less than 1.5%, and radii with a precision of about 1% and 3% for planet b and c, respectively. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 22 pages, 10 figures, Accepted for publication by Astronomy and Astrophysics

Unveiling the internal structure and formation history of the three planets transiting HIP 29442 (TOI-469) with CHEOPS

Authors: J. A. Egger, H. P. Osborn, D. Kubyshkina, C. Mordasini, Y. Alibert, M. N. Günther, M. Lendl, A. Brandeker, A. Heitzmann, A. Leleu, M. Damasso, A. Bonfanti, T. G. Wilson, S. G. Sousa, J. Haldemann, L. Delrez, M. J. Hooton, T. Zingales, R. Luque, R. Alonso, J. Asquier, T. Bárczy, D. Barrado Navascues, S. C. C. Barros, W. Baumjohann , et al. (69 additional authors not shown)

Abstract: Multiplanetary systems spanning the radius valley are ideal testing grounds for exploring the proposed explanations for the observed bimodality in the radius distribution of close-in exoplanets. One such system is HIP 29442 (TOI-469), an evolved K0V star hosting two super-Earths and a sub-Neptune. We observe HIP 29442 with CHEOPS for a total of 9.6 days, which we model jointly with 2 sectors of TE… ▽ More Multiplanetary systems spanning the radius valley are ideal testing grounds for exploring the proposed explanations for the observed bimodality in the radius distribution of close-in exoplanets. One such system is HIP 29442 (TOI-469), an evolved K0V star hosting two super-Earths and a sub-Neptune. We observe HIP 29442 with CHEOPS for a total of 9.6 days, which we model jointly with 2 sectors of TESS data to derive planetary radii of $3.410\pm0.046$, $1.551\pm0.045$ and $1.538\pm0.049$ R$_\oplus$ for planets b, c and d, which orbit HIP 29442 with periods of 13.6, 3.5 and 6.4 days. For planet d, this value deviates by more than 3 sigma from the median value reported in the discovery paper, leading us to conclude that caution is required when using TESS photometry to determine the radii of small planets with low per-transit S/N and large gaps between observations. Given the high precision of these new radii, combining them with published RVs from ESPRESSO and HIRES provides us with ideal conditions to investigate the internal structure and formation pathways of the planets in the system. We introduce the publicly available code plaNETic, a fast and robust neural network-based Bayesian internal structure modelling framework. We then apply hydrodynamic models to explore the upper atmospheric properties of these inferred structures. Finally, we identify planetary system analogues in a synthetic population generated with the Bern model for planet formation and evolution. Based on this analysis, we find that the planets likely formed on opposing sides of the water iceline from a protoplanetary disk with an intermediate solid mass. We finally report that the observed parameters of the HIP 29442 system are compatible with both a scenario where the second peak in the bimodal radius distribution corresponds to sub-Neptunes with a pure H/He envelope as well as a scenario with water-rich sub-Neptunes. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 30 pages, 17 figures, accepted for publication in A&A

The Predicament of Absorption-dominated Reionization II: Observational Estimate of the Clumping Factor at the End of Reionization

Authors: Frederick B. Davies, Sarah E. I. Bosman, Steven R. Furlanetto

Abstract: The history of reionization reflects the cumulative injection of ionizing photons by sources and the absorption of ionizing photons by sinks. The latter process is traditionally described in terms of a "clumping factor" which encodes the average quadratic increase in the recombination rate of dense gas within the cosmic web. The recent measurement of a short mean free path of ionizing photons from… ▽ More The history of reionization reflects the cumulative injection of ionizing photons by sources and the absorption of ionizing photons by sinks. The latter process is traditionally described in terms of a "clumping factor" which encodes the average quadratic increase in the recombination rate of dense gas within the cosmic web. The recent measurement of a short mean free path of ionizing photons from stacked quasar spectra at $z\simeq6$ has placed the importance of sinks under increased scrutiny, but its connection to the recombination rate is not immediately obvious. Here we present analytic arguments to connect the clumping factor to the mean free path by invoking ionization equilibrium within the ionized phase of the intergalactic medium at the end of (and after) reionization. We find that the latest mean free path and hydrogen photoionization rate measurements at $z=5-6$ imply a global clumping factor $C\approx12$, much higher than previous determinations from radiation-hydrodynamic simulations of the reionization process. Similar values of $C$ are also derived when applying the same procedure to observations at $2<z<5$. Compared to the traditional assumption of $C=3$, high-redshift galaxies must produce roughly twice as many ionizing photons ($\approx3$ photons per baryon) to reionize the universe by $z\sim6$. This additional requirement on the ionizing photon budget may help to reconcile the reionization history with JWST observations that suggest a far greater output of ionizing photons by the most distant galaxy populations. △ Less

Submitted 11 July, 2024; v1 submitted 26 June, 2024; originally announced June 2024.

Comments: 11 pages, 5+1 figures. Typos fixed, discussion improved, appendix added

Report number: NORDITA 2024-021

Quantifying the Precision of IGM Damping Wing Measurements Towards Quasars

Authors: Timo Kist, Joseph F. Hennawi, Frederick B. Davies

Abstract: We investigate the precision with which the Lyman-$α$ damping wing signature imprinted on the spectra of high-redshift quasars (QSOs) by the foreground neutral intergalactic medium (IGM) can measure the history of cosmic reionization. We leverage a novel inference pipeline based on a generative probabilistic model for the entire spectrum (both red- and blueward of the Lyman-$α$ line), accounting f… ▽ More We investigate the precision with which the Lyman-$α$ damping wing signature imprinted on the spectra of high-redshift quasars (QSOs) by the foreground neutral intergalactic medium (IGM) can measure the history of cosmic reionization. We leverage a novel inference pipeline based on a generative probabilistic model for the entire spectrum (both red- and blueward of the Lyman-$α$ line), accounting for all relevant sources of uncertainty - the stochasticity caused by patchy reionization, the impact of the quasar's ionizing radiation on the IGM, it's unknown intrinsic spectrum, and spectral noise. Performing fast JAX-based Hamiltonian Monte-Carlo (HMC) parameter inference, we precisely measure the underlying global IGM neutral fraction as well as the lifetime of the quasar. Running a battery of tests on over a thousand mocks, we find optimal precision when running the pipeline with a six parameter PCA continuum model (five coefficients and a normalization) on $\mathrm{S}/\mathrm{N} \sim 10$ spectra, binned to a $\sim 500\,\mathrm{km}/\mathrm{s}$ velocity pixel scale, and extending at least out to the C IV $λ\,1549\,\text{Å}$ emission line. After marginalizing out nuisance parameters associated with the quasar continuum, a single spectrum constrains the IGM neutral fraction to $28.0_{-8.8}^{+8.2}\,\%$ and the quasar lifetime to $0.80_{-0.55}^{+0.22}\,\mathrm{dex}$, improving notably towards spectra with a stronger IGM damping wing imprint. Higher precision can be achieved by averaging over statistical quasar samples. We identify two primary sources of uncertainty that contribute approximately equally to the total error budget: the uncertain quasar continuum model and the stochastic distribution of neutral regions arising from both the reionization topology and the location of the quasar's ionization front. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Submitted to MNRAS

Precisely Measuring the Cosmic Reionization History from IGM Damping Wings Towards Quasars

Authors: Joseph F. Hennawi, Timo Kist, Frederick B. Davies, John Tamanas

Abstract: We introduce a new approach for analyzing the IGM damping wings imprinted on the proximity zones of quasars in the epoch of reionization (EoR). Whereas past work has typically forgone the additional constraining power afforded by the blue side continuum ($1216\,Å \lesssim λ\lesssim 1280\,Å$) and/or opted not to model the large correlated IGM transmission fluctuations in the proximity zone (… ▽ More We introduce a new approach for analyzing the IGM damping wings imprinted on the proximity zones of quasars in the epoch of reionization (EoR). Whereas past work has typically forgone the additional constraining power afforded by the blue side continuum ($1216\,Å \lesssim λ\lesssim 1280\,Å$) and/or opted not to model the large correlated IGM transmission fluctuations in the proximity zone ($λ\lesssim 1216\,Å$), we construct a generative probabilistic model for the entire spectrum accounting for all sources of error - the stochasticity induced by patchy reionization, the impact of the quasar's ionizing radiation on the IGM, the unknown intrinsic spectrum of the quasar, and spectral noise. This principled Bayesian method allows us to marginalize out nuisance parameters associated with the quasar's radiation and its unknown intrinsic spectrum to precisely measure the IGM neutral fraction, $\langle x_{\rm HI}\rangle$. A key element of our analysis is the use of dimensionality reduction (DR) to describe the intrinsic quasar spectrum via a small number of nuisance parameters. Using a large sample of $15,559$ SDSS/BOSS quasars at $z \gtrsim 2.15$ we trained and quantified the performance of six distinct DR methods, and find that a six parameter PCA model (five coefficients plus a normalization) performs best, with complex machine learning approaches providing no advantage. By conducting statistical inference on 100 realistic mock EoR quasar spectra, we demonstrate the reliability of the credibility contours that we obtain on $\langle x_{\rm HI}\rangle$ and the quasar lifetime, $t_{\rm Q}$. The new method introduced here will transform IGM damping wings into a precision probe of reionization, on the same solid methodological and statistical footing as other precision cosmological measurements. △ Less

Submitted 21 June, 2024; v1 submitted 17 June, 2024; originally announced June 2024.

Comments: Submitted to MNRAS

Optical and near-infrared spectroscopy of quasars at $z>6.5$: public data release and composite spectrum

Authors: Silvia Onorato, Joseph F. Hennawi, Jan-Torge Schindler, Jinyi Yang, Feige Wang, Aaron J. Barth, Eduardo Bañados, Anna-Christina Eilers, Sarah E. I. Bosman, Frederick B. Davies, Bram P. Venemans, Chiara Mazzucchelli, Silvia Belladitta, Fabio Vito, Emanuele Paolo Farina, Irham T. Andika, Xiaohui Fan, Fabian Walter, Roberto Decarli, Masafusa Onoue, Riccardo Nanni

Abstract: We present optical and near-infrared (NIR) spectroscopic observations for a sample of $45$ quasars at $6.50 < z \leq 7.64$ with absolute magnitudes at $1450$ Å in the range $-28.82 \leq M_{1450} \leq -24.13$ and their composite spectrum. The median redshift and $M_{1450}$ of the quasars in the sample are $z_{\rm{median}}=6.71$ and $M_{1450,\rm{median}} \simeq -26.1$, respectively. The NIR spectra… ▽ More We present optical and near-infrared (NIR) spectroscopic observations for a sample of $45$ quasars at $6.50 < z \leq 7.64$ with absolute magnitudes at $1450$ Å in the range $-28.82 \leq M_{1450} \leq -24.13$ and their composite spectrum. The median redshift and $M_{1450}$ of the quasars in the sample are $z_{\rm{median}}=6.71$ and $M_{1450,\rm{median}} \simeq -26.1$, respectively. The NIR spectra are taken with echelle spectrographs, complemented with additional data from optical long slit instruments, and then reduced consistently using the open-source Python-based spectroscopic data reduction pipeline PypeIt. The median value of the mean signal-to-noise ratios of the spectra in J, H, and K band (median $\langle \rm{SNR}_λ \rangle$) is: median $\langle \rm{SNR}_{J} \rangle=9.7$, median $\langle \rm{SNR}_{H} \rangle=10.3$, and median $\langle \rm{SNR}_{K} \rangle=11.7$; demonstrating the good data quality. This work presents the largest medium/moderate-resolution sample of quasars at $z>6.5$ from ground-based instruments. Its homogeneity and reproducibility make it ideally suited for several scientific goals, i.e., the study of the quasar proximity zones and damping wings, the Ly$α$ forest, the intergalactic medium's metal content, as well as other properties such as the distribution of SMBH masses and Eddington ratios. Our composite spectrum is compared to others at both high and low-$z$ from the literature, showing differences in the strengths of many emission lines, probably due to differences in luminosity among the samples, but a consistent continuum slope, which proves that the same spectral features are preserved in quasars at different redshift ranges. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 18 pages, 10 figures, submitted to MNRAS

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

CHEOPS in-flight performance: A comprehensive look at the first 3.5 years of operations

Authors: A. Fortier, A. E. Simon, C. Broeg, G. Olofsson, A. Deline, T. G. Wilson, P. F. L. Maxted, A. Brandeker, A. Collier Cameron, M. Beck, A. Bekkelien, N. Billot, A. Bonfanti, G. Bruno, J. Cabrera, L. Delrez, B. -O. Demory, D. Futyan, H. -G. Florén, M. N. Günther, A. Heitzmann, S. Hoyer, K. G. Isaak, S. G. Sousa, M. Stalport , et al. (106 additional authors not shown)

Abstract: CHEOPS is a space telescope specifically designed to monitor transiting exoplanets orbiting bright stars. In September 2023, CHEOPS completed its nominal mission and remains in excellent operational conditions. The mission has been extended until the end of 2026. Scientific and instrumental data have been collected throughout in-orbit commissioning and nominal operations, enabling a comprehensive… ▽ More CHEOPS is a space telescope specifically designed to monitor transiting exoplanets orbiting bright stars. In September 2023, CHEOPS completed its nominal mission and remains in excellent operational conditions. The mission has been extended until the end of 2026. Scientific and instrumental data have been collected throughout in-orbit commissioning and nominal operations, enabling a comprehensive analysis of the mission's performance. In this article, we present the results of this analysis with a twofold goal. First, we aim to inform the scientific community about the present status of the mission and what can be expected as the instrument ages. Secondly, we intend for this publication to serve as a legacy document for future missions, providing insights and lessons learned from the successful operation of CHEOPS. To evaluate the instrument performance in flight, we developed a comprehensive monitoring and characterisation programme. It consists of dedicated observations that allow us to characterise the instrument's response. In addition to the standard collection of nominal science and housekeeping data, these observations provide input for detecting, modelling, and correcting instrument systematics, discovering and addressing anomalies, and comparing the instrument's actual performance with expectations. The precision of the CHEOPS measurements has enabled the mission objectives to be met and exceeded. Careful modelling of the instrumental systematics allows the data quality to be significantly improved during the light curve analysis phase, resulting in more precise scientific measurements. CHEOPS is compliant with the driving scientific requirements of the mission. Although visible, the ageing of the instrument has not affected the mission's performance. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: Accepted for publication in Astronomy and Astrophysics

HIP 41378 observed by CHEOPS: Where is planet d?

Authors: S. Sulis, L. Borsato, S. Grouffal, H. P. Osborn, A. Santerne, A. Brandeker, M. N. Günther, A. Heitzmann, M. Lendl, M. Fridlund, D. Gandolfi, Y. Alibert, R. Alonso, T. Bárczy, D. Barrado Navascues, S. C. Barros, W. Baumjohann, T. Beck, W. Benz, M. Bergomi, N. Billot, A. Bonfanti, C. Broeg, A. Collier Cameron, C. Corral van Damme , et al. (62 additional authors not shown)

Abstract: HIP 41378 d is a long-period planet that has only been observed to transit twice, three years apart, with K2. According to stability considerations and a partial detection of the Rossiter-McLaughlin effect, $P_\mathrm{d} = 278.36$ d has been determined to be the most likely orbital period. We targeted HIP 41378 d with CHEOPS at the predicted transit timing based on $P_\mathrm{d}= 278.36$ d, but th… ▽ More HIP 41378 d is a long-period planet that has only been observed to transit twice, three years apart, with K2. According to stability considerations and a partial detection of the Rossiter-McLaughlin effect, $P_\mathrm{d} = 278.36$ d has been determined to be the most likely orbital period. We targeted HIP 41378 d with CHEOPS at the predicted transit timing based on $P_\mathrm{d}= 278.36$ d, but the observations show no transit. We find that large ($>22.4$ hours) transit timing variations (TTVs) could explain this non-detection during the CHEOPS observation window. We also investigated the possibility of an incorrect orbital solution, which would have major implications for our knowledge of this system. If $P_\mathrm{d} \neq 278.36$ d, the periods that minimize the eccentricity would be $101.22$ d and $371.14$ d. The shortest orbital period will be tested by TESS, which will observe HIP 41378 in Sector 88 starting in January 2025. Our study shows the importance of a mission like CHEOPS, which today is the only mission able to make long observations (i.e., from space) to track the ephemeris of long-period planets possibly affected by large TTVs. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 686, L18 (2024)

Photo-dynamical characterisation of the TOI-178 resonant chain

Authors: A. Leleu, J. -B. Delisle, L. Delrez, E. M. Bryant, A. Brandeker, H. P. Osborn, N. Hara, T. G. Wilson, N. Billot, M. Lendl, D. Ehrenreich, H. Chakraborty, M. N. Günther, M. J. Hooton, Y. Alibert, R. Alonso, D. R. Alves, D. R. Anderson, I. Apergis, D. Armstrong, T. Bárczy, D. Barrado Navascues, S. C. C. Barros, M. P. Battley, W. Baumjohann , et al. (82 additional authors not shown)

Abstract: The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with radii ranging from 1.2 to 2.9 earth radius and orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. The fine-tuning and fragility of such orbital configurations ensure that no significant scattering or collision ev… ▽ More The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with radii ranging from 1.2 to 2.9 earth radius and orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. The fine-tuning and fragility of such orbital configurations ensure that no significant scattering or collision event has taken place since the formation and migration of the planets in the protoplanetary disc, hence providing important anchors for planet formation models. We aim to improve the characterisation of the architecture of this key system, and in particular the masses and radii of its planets. In addition, since this system is one of the few resonant chains that can be characterised by both photometry and radial velocities, we aim to use it as a test bench for the robustness of the planetary mass determination with each technique. We perform a global analysis of all available photometry and radial velocity. We also try different sets of priors on the masses and eccentricity, as well as different stellar activity models, to study their effects on the masses estimated by each method. We show how stellar activity is preventing us from obtaining a robust mass estimation for the three outer planets using radial velocity data alone. We also show that our joint photo-dynamical and radial velocity analysis resulted in a robust mass determination for planets c to g, with precision of 12% for the mass of planet c, and better than 10% for planets d to g. The new precisions on the radii range from 2 to 3%. The understanding of this synergy between photometric and radial velocity measurements will be valuable during the PLATO mission. We also show that TOI-178 is indeed currently locked in the resonant configuration, librating around an equilibrium of the chain. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Damping Wing-Like Features in the Stacked Ly$α$ Forest: Potential Neutral Hydrogen Islands at $z<6$

Authors: Yongda Zhu, George D. Becker, Sarah E. I. Bosman, Christopher Cain, Laura C. Keating, Fahad Nasir, Valentina D'Odorico, Eduardo Bañados, Fuyan Bian, Manuela Bischetti, James S. Bolton, Huanqing Chen, Anson D'Aloisio, Frederick B. Davies, Rebecca L. Davies, Anna-Christina Eilers, Xiaohui Fan, Prakash Gaikwad, Bradley Greig, Martin G. Haehnelt, Girish Kulkarni, Samuel Lai, Ewald Puchwein, Yuxiang Qin, Emma V. Ryan-Weber , et al. (6 additional authors not shown)

Abstract: Recent quasar absorption line observations suggest that reionization may end as late as $z \approx 5.3$. As a means to search for large neutral hydrogen islands at $z<6$, we revisit long dark gaps in the Ly$β$ forest in VLT/X-Shooter and Keck/ESI quasar spectra. We stack the Ly$α$ forest corresponding to both edges of these Ly$β$ dark gaps and identify a damping wing-like extended absorption profi… ▽ More Recent quasar absorption line observations suggest that reionization may end as late as $z \approx 5.3$. As a means to search for large neutral hydrogen islands at $z<6$, we revisit long dark gaps in the Ly$β$ forest in VLT/X-Shooter and Keck/ESI quasar spectra. We stack the Ly$α$ forest corresponding to both edges of these Ly$β$ dark gaps and identify a damping wing-like extended absorption profile. The average redshift of the stacked forest is $z=5.8$. By comparing these observations with reionization simulations, we infer that such a damping wing-like feature can be naturally explained if these gaps are at least partially created by neutral islands. Conversely, simulated dark gaps lacking neutral hydrogen struggle to replicate the observed damping wing features. Furthermore, this damping wing-like profile implies that the volume-averaged neutral hydrogen fraction must be $\langle x_{\rm HI} \rangle \geq 6.1 \pm 3.9\%$ at $z = 5.8$. Our results offer robust evidence that reionization extends below $z=6$. △ Less

Submitted 28 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 8 pages, 5 figures, 1 table; accepted for publication in MNRAS Letters

Damping wings in the Lyman-α forest: a model-independent measurement of the neutral fraction at 5.4<z<6.1

Authors: Benedetta Spina, Sarah E. I. Bosman, Frederick B. Davies, Prakash Gaikwad, Yongda Zhu

Abstract: Recent observations have positioned the endpoint of the Epoch of Reionisation (EoR) at redshift $z \sim 5.3$. However, observations of the Lyman-$α$ forest have not yet been able to discern whether reionisation occurred slowly and late, with substantial neutral hydrogen persisting at redshift $\sim 6$, or rapidly and earlier, with the apparent late end driven by the fluctuating UV background. Gunn… ▽ More Recent observations have positioned the endpoint of the Epoch of Reionisation (EoR) at redshift $z \sim 5.3$. However, observations of the Lyman-$α$ forest have not yet been able to discern whether reionisation occurred slowly and late, with substantial neutral hydrogen persisting at redshift $\sim 6$, or rapidly and earlier, with the apparent late end driven by the fluctuating UV background. Gunn-Peterson (GP) absorption troughs are solid indicators that reionisation is not complete until $z=5.3$, but whether they contain significantly neutral gas has not yet been proven. We aim to answer this question by directly measuring, for the first time, the neutral hydrogen fraction ($x_\mathrm{HI}$) at the end of the EoR ($5 \lesssim z \lesssim 6$) in high-redshift quasars spectra. For high neutral fractions $x_\mathrm{HI}\gtrsim0.1$, GP troughs exhibit damping wing (DW) absorption extending over $1000$ km s$^{-1}$ beyond the troughs. While conclusively detected in Lyman-$α$ emission lines of quasars at $z\geq7$, DWs are challenging to observe in the general Lyman-$α$ forest due to absorption complexities and small-scale stochastic transmission features. We report the first successful identification of the stochastic DW signal adjacent to GP troughs at redshifts $z=5.6$ through careful stacking of the dark gaps in Lyman-$α$ forest. We use the signal to present a measurement of the corresponding global $x_\mathrm{HI}=0.19\pm0.07$ $(_{-0.16}^{+0.11})$ at $1σ$ $(2σ)$ at $z=5.6$ and a limit $x_\mathrm{HI}<0.44$ at $z=5.9$. The detection of this signal demonstrates the existence of substantially neutral islands near the conclusion of the EoR, unequivocally signaling a late-and-slow reionization scenario. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 7 pages, 4 figures, submitted to A&A Letters

A Spatially Resolved [CII] Survey of 31 $z\sim7$ Massive Galaxies Hosting Luminous Quasars

Authors: Feige Wang, Jinyi Yang, Xiaohui Fan, Bram Venemans, Roberto Decarli, Eduardo Bañados, Fabian Walter, Aaron J. Barth, Fuyan Bian, Frederick B. Davies, Anna-Christina Eilers, Emanuele Paolo Farina, Joseph F. Hennawi, Jiang-Tao Li, Chiara Mazzucchelli, Ran Wang, Xue-Bing Wu, Minghao Yue

Abstract: The [CII] 158 $μ$m emission line and the underlying far-infrared (FIR) dust continuum are important tracers for studying star formation and kinematic properties of early galaxies. We present a survey of the [CII] emission lines and FIR continua of 31 luminous quasars at $z>6.5$ using the Atacama Large Millimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) at sub-arcsec resoluti… ▽ More The [CII] 158 $μ$m emission line and the underlying far-infrared (FIR) dust continuum are important tracers for studying star formation and kinematic properties of early galaxies. We present a survey of the [CII] emission lines and FIR continua of 31 luminous quasars at $z>6.5$ using the Atacama Large Millimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) at sub-arcsec resolution. This survey more than doubles the number of quasars with [CII] and FIR observations at these redshifts and enables statistical studies of quasar host galaxies deep into the epoch of reionization. We detect [CII] emission in 27 quasar hosts with a luminosity range of $L_{\rm [CII]}=(0.3-5.5)\times10^9~L_\odot$ and detect the FIR continuum of 28 quasar hosts with a luminosity range of $L_{\rm FIR}=(0.5-13.0)\times10^{12}~L_\odot$. Both $L_{\rm [CII]}$ and $L_{\rm FIR}$ are correlated ($ρ\simeq0.4$) with the quasar bolometric luminosity, albeit with substantial scatter. The quasar hosts detected by ALMA are clearly resolved with a median diameter of $\sim$5 kpc. About 40% of the quasar host galaxies show a velocity gradient in [CII] emission, while the rest show either dispersion-dominated or disturbed kinematics. Basic estimates of the dynamical masses of the rotation-dominated host galaxies yield $M_{\rm dyn}=(0.1-7.5)\times10^{11}~M_\odot$. Considering our findings alongside those of literature studies, we found that the ratio between $M_{\rm BH}$ and $M_{\rm dyn}$ is about ten times higher than that of local $M_{\rm BH}-M_{\rm dyn}$ relation on average but with substantial scatter (the ratio difference ranging from $\sim$0.6 to 60) and large uncertainties. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: accepted for publication in ApJ

IGM damping wing constraints on the tail end of reionisation from the enlarged XQR-30 sample

Authors: Bradley Greig, Andrei Mesinger, Eduardo Bañados, George D. Becker, Sarah E. I. Bosman, Huanqing Chen, Frederick B. Davies, Valentina D'Odorico, Anna-Christina Eilers, Simona Gallerani, Martin G. Haehnelt, Laura Keating, Samuel Lai, Yuxiang Qin, Emma Ryan-Weber, Sindhu Satyavolu, Feige Wang, Jinyi Yang, Yongda Zhu

Abstract: The attenuation of Ly$α$ photons by neutral hydrogen in the intergalactic medium (IGM) at $z\gtrsim5$ continues to be a powerful probe for studying the epoch of reionisation. Given a framework to estimate the intrinsic (true) Ly$α$ emission of high-$z$ sources, one can infer the ionisation state of the IGM during reionisation. In this work, we use the enlarged XQR-30 sample of 42 high-resolution a… ▽ More The attenuation of Ly$α$ photons by neutral hydrogen in the intergalactic medium (IGM) at $z\gtrsim5$ continues to be a powerful probe for studying the epoch of reionisation. Given a framework to estimate the intrinsic (true) Ly$α$ emission of high-$z$ sources, one can infer the ionisation state of the IGM during reionisation. In this work, we use the enlarged XQR-30 sample of 42 high-resolution and high-SNR QSO spectra between $5.8\lesssim\,z\lesssim\,6.6$ obtained with VLT/X-Shooter to place constraints on the IGM neutral fraction. This is achieved using our existing Bayesian QSO reconstruction framework which accounts for uncertainties such as the: (i) posterior distribution of predicted intrinsic Ly$α$ emission profiles (obtained via covariance matrix reconstruction of the Ly$α$ and N V emission lines from unattenuated high-ionisation emission line profiles; C IV, Si IV + O IV] and C III]) and (ii) distribution of ionised regions within the IGM using synthetic damping wing profiles drawn from a $1.6^3$ Gpc$^3$ reionisation simulation. Following careful quality control, we used 23 of the 42 available QSOs to obtain constraints/limits on the IGM neutral fraction during the tail-end of reionisation. Our median and 68th percentile constraints on the IGM neutral fraction are: $0.20\substack{+0.14\\-0.12}$ and $0.29\substack{+0.14\\-0.13}$ at $z = 6.15$~and 6.35. Further, we also report 68th percentile upper-limits of $\bar{x}_{\mathrm{H\,{\scriptscriptstyle I}}} < 0.21$, 0.20, 0.21 and 0.18 at $z = 5.8, 5.95, 6.05$~and 6.55. These results imply reionisation is still ongoing at $5.8\lesssim\,z\lesssim\,6.55$, consistent with previous results from XQR-30 (dark fraction and Ly$α$ forest) along with other observational probes considered in the literature. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 21 pages, 9 figures and 2 tables. Accepted for publication in MNRAS

The formation of transiting circumplanetary debris discs from the disruption of satellite systems during planet-planet scattering

Authors: Alexander J. Mustill, Melvyn B. Davies, Matthew A. Kenworthy

Abstract: Several stars show deep transits consistent with discs of roughly 1 Solar radius seen at moderate inclinations, likely surrounding planets on eccentric orbits. We show that this configuration arises naturally as a result of planet-planet scattering when the planets possess satellite systems. Planet-planet scattering explains the orbital eccentricities of the discs' host bodies, while the close enc… ▽ More Several stars show deep transits consistent with discs of roughly 1 Solar radius seen at moderate inclinations, likely surrounding planets on eccentric orbits. We show that this configuration arises naturally as a result of planet-planet scattering when the planets possess satellite systems. Planet-planet scattering explains the orbital eccentricities of the discs' host bodies, while the close encounters during scattering lead to the exchange of satellites between planets and/or their destabilisation. This leads to collisions between satellites and their tidal disruption close to the planet. Both of these events lead to large quantities of debris being produced, which in time will settle into a disc such as those observed. The mass of debris required is comparable to a Ceres-sized satellite. Through N-body simulations of planets with clones of the Galilean satellite system undergoing scattering, we show that 90 percent of planets undergoing scattering will possess debris from satellite destruction. Extrapolating to smaller numbers of satellites suggests that tens of percent of such planets should still possess circumplanetary debris discs. The debris trails arising from these events are often tilted at tens of degrees to the planetary orbit, consistent with the inclinations of the observed discs. Disruption of satellite systems during scattering thus simultaneously explains the existence of debris, the tilt of the discs, and the eccentricity of the planets they orbit. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: Accepted to MNRAS. Github repository at https://github.com/AJMustill/moons

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 530, Issue 4, June 2024, Pages 3606-3616

Characterisation of the TOI-421 planetary system using CHEOPS, TESS, and archival radial velocity data

Authors: A. F. Krenn, D. Kubyshkina, L. Fossati, J. A. Egger, A. Bonfanti, A. Deline, D. Ehrenreich, M. Beck, W. Benz, J. Cabrera, T. G. Wilson, A. Leleu, S. G. Sousa, V. Adibekyan, A. C. M. Correira, Y. Alibert, L. Delrez, M. Lendl, J. A. Patel, J. Venturini, R. Alonso, G. Anglada, J. Asquier, T. Bárczy, D. Barrado Navascues , et al. (66 additional authors not shown)

Abstract: The TOI-421 planetary system contains two sub-Neptune-type planets and is a prime target to study the formation and evolution of planets and their atmospheres. The inner planet is especially interesting as the existence of a hydrogen-dominated atmosphere at its orbital separation cannot be explained by current formation models without previous orbital migration. We jointly analysed photometric dat… ▽ More The TOI-421 planetary system contains two sub-Neptune-type planets and is a prime target to study the formation and evolution of planets and their atmospheres. The inner planet is especially interesting as the existence of a hydrogen-dominated atmosphere at its orbital separation cannot be explained by current formation models without previous orbital migration. We jointly analysed photometric data of three TESS sectors and six CHEOPS visits as well as 156 radial velocity data points to retrieve improved planetary parameters. We also searched for TTVs and modelled the interior structure of the planets. Finally, we simulated the evolution of the primordial H-He atmospheres of the planets using two different modelling frameworks. We determine the planetary radii and masses of TOI-421 b and c to be $R_{\rm b} = 2.64 \pm 0.08 \, R_{\oplus}$, $M_{\rm b} = 6.7 \pm 0.6 \, M_{\oplus}$, $R_{\rm c} = 5.09 \pm 0.07 \, R_{\oplus}$, and $M_{\rm c} = 14.1 \pm 1.4 \, M_{\oplus}$. We do not detect any statistically significant TTV signals. Assuming the presence of a hydrogen-dominated atmosphere, the interior structure modelling results in both planets having extensive envelopes. While the modelling of the atmospheric evolution predicts for TOI-421 b to have lost any primordial atmosphere that it could have accreted at its current orbital position, TOI-421 c could have started out with an initial atmospheric mass fraction somewhere between 10 and 35%. We conclude that the low observed mean density of TOI-421 b can only be explained by either a bias in the measured planetary parameters (e.g. driven by high-altitude clouds) and/or in the context of orbital migration. We also find that the results of atmospheric evolution models are strongly dependent on the employed planetary structure model. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 21 pages, 12 figures, accepted for publication in A&A

Deep Learning the Intergalactic Medium using Lyman-alpha Forest at $ 4 \leq z \leq 5$

Authors: Fahad Nasir, Prakash Gaikwad, Frederick B. Davies, James S. Bolton, Ewald Puchwein, Sarah E. I. Bosman

Abstract: Unveiling the thermal history of the intergalactic medium (IGM) at $4 \leq z \leq 5$ holds the potential to reveal early onset HeII reionization or lingering thermal fluctuations from HI reionization. We set out to reconstruct the IGM gas properties along simulated Lyman-alpha forest data on pixel-by-pixel basis, employing deep Bayesian neural networks. Our approach leverages the Sherwood-Relics s… ▽ More Unveiling the thermal history of the intergalactic medium (IGM) at $4 \leq z \leq 5$ holds the potential to reveal early onset HeII reionization or lingering thermal fluctuations from HI reionization. We set out to reconstruct the IGM gas properties along simulated Lyman-alpha forest data on pixel-by-pixel basis, employing deep Bayesian neural networks. Our approach leverages the Sherwood-Relics simulation suite, consisting of diverse thermal histories, to generate mock spectra. Our convolutional and residual networks with likelihood metric predicts the Ly$α$ optical depth-weighted density or temperature for each pixel in the Ly$α$ forest skewer. We find that our network can successfully reproduce IGM conditions with high fidelity across range of instrumental signal-to-noise. These predictions are subsequently translated into the temperature-density plane, facilitating the derivation of reliable constraints on thermal parameters. This allows us to estimate temperature at mean cosmic density, $T_{\rm 0}$ with one sigma confidence $δT_{\rm 0} \sim 1000{\rm K}$ using only one $20$Mpc/h sightline ($Δz\simeq 0.04$) with a typical reionization history. Existing studies utilize redshift pathlength comparable to $Δz\simeq 4$ for similar constraints. We can also provide more stringent constraints on the slope ($1σ$ confidence interval $δ{\rm γ} \lesssim 0.1$) of the IGM temperature-density relation as compared to other traditional approaches. We test the reconstruction on a single high signal-to-noise observed spectrum ($20$ Mpc/h segment), and recover thermal parameters consistent with current measurements. This machine learning approach has the potential to provide accurate yet robust measurements of IGM thermal history at the redshifts in question. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 17 pages, Submitted to MNRAS

Blind QSO reconstruction challenge: Exploring methods to reconstruct the Ly$α$ emission line of QSOs

Authors: Bradley Greig, Sarah E. I. Bosman, Frederick B. Davies, Dominika Ďurovčíková, Hassan Fathivavsari, Bin Liu, Romain A. Meyer, Zechang Sun, Valentina D'Odorico, Simona Gallerani, Andrei Mesinger, Yuan-Sen Ting

Abstract: Reconstructing the intrinsic Ly$α$ line flux from high-$z$ QSOs can place constraints on the neutral hydrogen content of the intergalactic medium during reionisation. There are now $\gtrsim10$ different Ly$α$ reconstruction pipelines using different methodologies to predict the Ly$α$ line flux from correlations with the spectral information redward of Ly$α$. However, there have been few attempts t… ▽ More Reconstructing the intrinsic Ly$α$ line flux from high-$z$ QSOs can place constraints on the neutral hydrogen content of the intergalactic medium during reionisation. There are now $\gtrsim10$ different Ly$α$ reconstruction pipelines using different methodologies to predict the Ly$α$ line flux from correlations with the spectral information redward of Ly$α$. However, there have been few attempts to directly compare the performance of these pipelines. Therefore, we devised a blind QSO challenge to compare these reconstruction pipelines on a uniform set of objects. Each author was provided de-identified, observed rest-frame QSO spectra with spectral information only redward of 1260Å rest-frame to ensure unbiased reconstruction. We constructed two samples of 30 QSOs, from X-Shooter and SDSS both spanning $3.5<z<4.5$. Importantly, the purpose of this comparison study was not to champion a single, best performing reconstruction pipeline but rather to explore the relative performance of these pipelines over a range of QSOs with broad observational characteristics to infer general trends. In summary, we find machine learning approaches in general provide the strongest ``best guesses" but underestimate the accompanying statistical uncertainty, although these can be recalibrated, whilst pipelines that decompose the spectral information, for example principal component or factor analysis generally perform better at predicting the Ly$α$ profile. Further, we found that reconstruction pipelines trained on SDSS QSOs performed similarly on average for both the X-Shooter and SDSS samples indicating no discernible biases owing to differences in the observational characteristics of the training set or QSO being reconstructed, although the recovered distributions of reconstructions for X-Shooter were broader likely due to an increased fraction of outliers. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: 30 pages, 16 figures, 4 tables. Submitted to MNRAS, comments welcome

Detailed cool star flare morphology with CHEOPS and TESS

Authors: G. Bruno, I. Pagano, G. Scandariato, H. -G. Florén, A. Brandeker, G. Olofsson, P. F. L. Maxted, A. Fortier, S. G. Sousa, S. Sulis, V. Van Grootel, Z. Garai, A. Boldog, L. Kriskovics, M. Gy. Szabó, D. Gandolfi, Y. Alibert, R. Alonso, T. Bárczy, D. Barrado Navascues, S. C. C. Barros, W. Baumjohann, M. Beck, T. Beck, W. Benz , et al. (57 additional authors not shown)

Abstract: Context. White-light stellar flares are proxies for some of the most energetic types of flares, but their triggering mechanism is still poorly understood. As they are associated with strong X and UV emission, their study is particularly relevant to estimate the amount of high-energy irradiation onto the atmospheres of exoplanets, especially those in their stars' habitable zone. Aims. We used the h… ▽ More Context. White-light stellar flares are proxies for some of the most energetic types of flares, but their triggering mechanism is still poorly understood. As they are associated with strong X and UV emission, their study is particularly relevant to estimate the amount of high-energy irradiation onto the atmospheres of exoplanets, especially those in their stars' habitable zone. Aims. We used the high-cadence, high-photometric capabilities of the CHEOPS and TESS space telescopes to study the detailed morphology of white-light flares occurring in a sample of 130 late-K and M stars, and compared our findings with results obtained at a lower cadence. We developed dedicated software for this purpose. Results. Multi-peak flares represent a significant percentage ($\gtrsim 30$\%) of the detected outburst events. Our findings suggest that high-impulse flares are more frequent than suspected from lower-cadence data, so that the most impactful flux levels that hit close-in exoplanets might be more time-limited than expected. We found significant differences in the duration distributions of single-peak and complex flare components, but not in their peak luminosity. A statistical analysis of the flare parameter distributions provides marginal support for their description with a log-normal instead of a power-law function, leaving the door open to several flare formation scenarios. We tentatively confirmed previous results about quasi-periodic pulsations in high-cadence photometry, report the possible detection of a pre-flare dip, and did not find hints of photometric variability due to an undetected flare background. Conclusions. The high-cadence study of stellar hosts might be crucial to evaluate the impact of their flares on close-in exoplanets, as their impulsive phase emission might otherwise be incorrectly estimated. Future telescopes such as PLATO and Ariel will help in this respect. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 28 pages, 25 figures, 4 tables, to be published in Astronomy & Astrophysics

Precise characterisation of HD 15337 with CHEOPS: a laboratory for planet formation and evolution

Authors: N. M. Rosário, O. D. S. Demangeon, S. C. C. Barros, D. Gandolfi, J. A. Egger, L. M. Serrano, H. P. Osborn, M. Beck, W. Benz, H. -G. Florén, P. Guterman, T. G. Wilson, Y. Alibert, L. Fossati, M. J. Hooton, L. Delrez, N. C. Santos, S. G. Sousa, A. Bonfanti, S. Salmon, V. Adibekyan, A. Nigioni, J. Venturini, R. Alonso, G. Anglada , et al. (68 additional authors not shown)

Abstract: We aim to constrain the internal structure and composition of HD 15337 b and c, two short-period planets situated on opposite sides of the radius valley, using new transit photometry and radial velocity data. We acquire 6 new transit visits with the CHaracterising ExOPlanet Satellite (CHEOPS) and 32 new radial velocity measurements from the High Accuracy Radial Velocity Planet Searcher (HARPS) to… ▽ More We aim to constrain the internal structure and composition of HD 15337 b and c, two short-period planets situated on opposite sides of the radius valley, using new transit photometry and radial velocity data. We acquire 6 new transit visits with the CHaracterising ExOPlanet Satellite (CHEOPS) and 32 new radial velocity measurements from the High Accuracy Radial Velocity Planet Searcher (HARPS) to improve the accuracy of the mass and radius estimates for both planets. We reanalyse light curves from TESS sectors 3 and 4 and analyse new data from sector 30, correcting for long-term stellar activity. Subsequently, we perform a joint fit of the TESS and CHEOPS light curves, and all available RV data from HARPS and the Planet Finder Spectrograph (PFS). Our model fits the planetary signals, the stellar activity signal and the instrumental decorrelation model for the CHEOPS data simultaneously. The stellar activity was modelled using a Gaussian-process regression on both the RV and activity indicators. We finally employ a Bayesian retrieval code to determine the internal composition and structure of the planets. We derive updated and highly precise parameters for the HD 15337 system. Our improved precision on the planetary parameters makes HD 15337 b one of the most precisely characterised rocky exoplanets, with radius and mass measurements achieving a precision better than 2\% and 7\%, respectively. We are able to improve the precision of the radius measurement of HD 15337 c to 3\%. Our results imply that the composition of HD 15337 b is predominantly rocky, while HD 15337 c exhibits a gas envelope with a mass of at least $0.01\ M_\oplus$.Our results lay the groundwork for future studies, which can further unravel the atmospheric evolution of these exoplanets and give new insights into their composition and formation history and the causes behind the radius gap. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 17 pages, including appendix

Journal ref: A&A 686, A282 (2024)

The tidal deformation and atmosphere of WASP-12b from its phase curve

Authors: B. Akinsanmi, S. C. C. Barros, M. Lendl, L. Carone, P. E. Cubillos, A. Bekkelien, A. Fortier, H. -G. Florén, A. Collier Cameron, G. Boué, G. Bruno, B. -O. Demory, A. Brandeker, S. G. Sousa, T. G. Wilson, A. Deline, A. Bonfanti, G. Scandariato, M. J. Hooton, A. C. M. Correia, O. D. S. Demangeon, A. M. S. Smith, V. Singh, Y. Alibert, R. Alonso , et al. (63 additional authors not shown)

Abstract: Ultra-hot Jupiters present a unique opportunity to understand the physics and chemistry of planets at extreme conditions. WASP-12b stands out as an archetype of this class of exoplanets. We performed comprehensive analyses of the transits, occultations, and phase curves of WASP-12b by combining new CHEOPS observations with previous TESS and Spitzer data to measure the planet's tidal deformation, a… ▽ More Ultra-hot Jupiters present a unique opportunity to understand the physics and chemistry of planets at extreme conditions. WASP-12b stands out as an archetype of this class of exoplanets. We performed comprehensive analyses of the transits, occultations, and phase curves of WASP-12b by combining new CHEOPS observations with previous TESS and Spitzer data to measure the planet's tidal deformation, atmospheric properties, and orbital decay rate. The planet was modeled as a triaxial ellipsoid parameterized by the second-order fluid Love number, $h_2$, which quantifies its radial deformation and provides insight into the interior structure. We measured the tidal deformation of WASP-12b and estimated a Love number of $h_2=1.55_{-0.49}^{+0.45}$ (at 3.2$σ$) from its phase curve. We measured occultation depths of $333\pm24$ppm and $493\pm29$ppm in the CHEOPS and TESS bands, respectively, while the dayside emission spectrum indicates that CHEOPS and TESS probe similar pressure levels in the atmosphere at a temperature of 2900K. We also estimated low geometric albedos of $0.086\pm0.017$ and $0.01\pm0.023$ in the CHEOPS and TESS passbands, respectively, suggesting the absence of reflective clouds in the dayside of the WASP-12b. The CHEOPS occultations do not show strong evidence for variability in the dayside atmosphere of the planet. Finally, we refine the orbital decay rate by 12% to a value of -30.23$\pm$0.82 ms/yr. WASP-12b becomes the second exoplanet, after WASP-103b, for which the Love number has been measured (at 3$sigma$) from the effect of tidal deformation in the light curve. However, constraining the core mass fraction of the planet requires measuring $h_2$ with a higher precision. This can be achieved with high signal-to-noise observations with JWST since the phase curve amplitude, and consequently the induced tidal deformation effect, is higher in the infrared. △ Less

Submitted 20 February, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: accepted for publication in A&A

A Massive Protocluster Anchored by a Luminous Quasar at $z=6.63$

Authors: Feige Wang, Jinyi Yang, Joseph F. Hennawi, Xiaohui Fan, Minghao Yue, Eduardo Bañados, Shane Bechtel, Fuyan Bian, Sarah Bosman, Jaclyn B. Champagne, Frederick B. Davies, Roberto Decarli, Emanuele Paolo Farina, Chiara Mazzucchelli, Bram Venemans, Fabian Walter

Abstract: Protoclusters, the progenitors of galaxy clusters, trace large scale structures in the early Universe and are important to our understanding of structure formation and galaxy evolution. To date, only a handful of protoclusters have been identified in the Epoch of Reionization (EoR). As one of the rarest populations in the early Universe, distant quasars that host active supermassive black holes ar… ▽ More Protoclusters, the progenitors of galaxy clusters, trace large scale structures in the early Universe and are important to our understanding of structure formation and galaxy evolution. To date, only a handful of protoclusters have been identified in the Epoch of Reionization (EoR). As one of the rarest populations in the early Universe, distant quasars that host active supermassive black holes are thought to reside in the most massive dark matter halos at that cosmic epoch, and could thus potentially pinpoint some of the earliest protoclusters. In this letter, we report the discovery of a massive protocluster around a luminous quasar at $z=6.63$. This protocluster is anchored by the quasar, and includes three [CII] emitters at $z\sim6.63$, 12 spectroscopically confirmed Ly$α$ emitters (LAEs) at $6.54<z\le6.64$, and a large number of narrow-band imaging selected LAE candidates at the same redshift. This structure has an overall overdensity of $δ=3.3^{+1.1}_{-0.9}$ within $\sim35\times74$ cMpc$^2$ on the sky and an extreme overdensity of $δ>30$ in its central region (i.e., $R\lesssim2$ cMpc). We estimate that this protocluster will collapse into a galaxy cluster with a mass of $6.9^{+1.2}_{-1.4}\times10^{15}~M_\odot$ at the current epoch, more massive than the most massive clusters known in the local Universe such as Coma. In the quasar vicinity, we discover a double-peaked LAE which implies that the quasar has a UV lifetime greater than 0.8 Myrs and has already ionized its surrounding intergalactic medium. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL

Discovery of two warm mini-Neptunes with contrasting densities orbiting the young K3V star TOI-815

Authors: Angelica Psaridi, Hugh Osborn, François Bouchy, Monika Lendl, Léna Parc, Nicolas Billot, Christopher Broeg, Sérgio G. Sousa, Vardan Adibekyan, Omar Attia, Andrea Bonfanti, Hritam Chakraborty, Karen A. Collins, Jeanne Davoult, Elisa Delgado-Mena, Nolan Grieves, Tristan Guillot, Alexis Heitzmann, Ravit Helled, Coel Hellier, Jon M. Jenkins, Henrik Knierim, Andreas Krenn, JackJ. Lissauer, Rafael Luque , et al. (108 additional authors not shown)

Abstract: We present the discovery and characterization of two warm mini-Neptunes transiting the K3V star TOI-815 in a K-M binary system. Analysis of the spectra and rotation period reveal it to be a young star with an age of $200^{+400}_{-200}$Myr. TOI-815b has a 11.2-day period and a radius of 2.94$\pm$0.05$\it{R_{\rm\mathrm{\oplus}}}$ with transits observed by TESS, CHEOPS, ASTEP, and LCOGT. The outer pl… ▽ More We present the discovery and characterization of two warm mini-Neptunes transiting the K3V star TOI-815 in a K-M binary system. Analysis of the spectra and rotation period reveal it to be a young star with an age of $200^{+400}_{-200}$Myr. TOI-815b has a 11.2-day period and a radius of 2.94$\pm$0.05$\it{R_{\rm\mathrm{\oplus}}}$ with transits observed by TESS, CHEOPS, ASTEP, and LCOGT. The outer planet, TOI-815c, has a radius of 2.62$\pm$0.10$\it{R_{\rm\mathrm{\oplus}}}$, based on observations of three non-consecutive transits with TESS, while targeted CHEOPS photometry and radial velocity follow-up with ESPRESSO were required to confirm the 35-day period. ESPRESSO confirmed the planetary nature of both planets and measured masses of 7.6$\pm$1.5 $\it{M_{\rm \mathrm{\oplus}}}$ ($ρ_\mathrm{P}$=1.64$^{+0.33}_{-0.31}$gcm$^{-3}$) and 23.5$\pm$2.4$\it{M_{\rm\mathrm{\oplus}}}$ ($ρ_\mathrm{P}$=7.2$^{+1.1}_{-1.0}$gcm$^{-3}$) respectively. Thus, the planets have very different masses, unlike the usual similarity of masses in compact multi-planet systems. Moreover, our statistical analysis of mini-Neptunes orbiting FGK stars suggests that weakly irradiated planets tend to have higher bulk densities compared to those suffering strong irradiation. This could be ascribed to their cooler atmospheres, which are more compressed and denser. Internal structure modeling of TOI-815b suggests it likely has a H-He atmosphere constituting a few percent of the total planet mass, or higher if the planet is assumed to have no water. In contrast, the measured mass and radius of TOI-815c can be explained without invoking any atmosphere, challenging planetary formation theories. Finally, we infer from our measurements that the star is viewed close to pole-on, which implies a spin-orbit misalignment at the 3$σ$ level. △ Less

Submitted 30 January, 2024; v1 submitted 28 January, 2024; originally announced January 2024.

Comments: 24 pages, 27 figures, 6 tables

Recognizing Blazars Using Radio Morphology from the VLA Sky Survey

Authors: Zhang-Liang Xie, Eduardo Banados, Silvia Belladitta, Chiara Mazzucchelli, Jan-Torge Schindler, Frederick B. Davies, Bram P. Venemans

Abstract: Blazars are radio-loud Active Galactic Nuclei (AGN) whose jets have a very small angle to our line of sight. Observationally, the radio emission are mostly compact or a compact-core with a 1-sided jet. With 2.5$^{\prime\prime}$ resolution at 3 GHz, the Very Large Array Sky Survey (VLASS) enables us to resolve the structure of some blazar candidates in the sky north of Decl. $-40$ deg. We introduce… ▽ More Blazars are radio-loud Active Galactic Nuclei (AGN) whose jets have a very small angle to our line of sight. Observationally, the radio emission are mostly compact or a compact-core with a 1-sided jet. With 2.5$^{\prime\prime}$ resolution at 3 GHz, the Very Large Array Sky Survey (VLASS) enables us to resolve the structure of some blazar candidates in the sky north of Decl. $-40$ deg. We introduce an algorithm to classify radio sources as either blazar-like or non-blazar-like based on their morphology in the VLASS images. We apply our algorithm to three existing catalogs, including one of known blazars (Roma-BzCAT) and two of blazar candidates identified by WISE colors and radio emission (WIBRaLS, KDEBLLACS). We show that in all three catalogs, there are objects with morphology inconsistent with being blazars. Considering all the catalogs, more than 12% of the candidates are unlikely to be blazars, based on this analysis. Notably, we show that 3% of the Roma-BzCAT "confirmed" blazars could be a misclassification based on their VLASS morphology. The resulting table with all sources and their radio morphological classification is available online. △ Less

Submitted 9 January, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 21 pages, 11 figures, accepted for publication in The Astrophysical Journal

The effect of winds on atmospheric layers of red supergiants II. Modelling VLTI/GRAVITY and MATISSE observations of AH Sco, KW Sgr, V602 Car, CK Car and V460 Car

Authors: G. González-Torà, M. Wittkowski, B. Davies, B. Plez

Abstract: Mass loss plays a crucial role in the lives of massive stars, especially as the star leaves the main sequence and evolves to the red supergiant (RSG) phase. However, the physical processes that trigger mass-loss events in RSGs are not well understood. Recently, we showed that adding a static wind to a MARCS atmospheric model can accurately reproduce observed extensions in the atmospheres of RSGs.… ▽ More Mass loss plays a crucial role in the lives of massive stars, especially as the star leaves the main sequence and evolves to the red supergiant (RSG) phase. However, the physical processes that trigger mass-loss events in RSGs are not well understood. Recently, we showed that adding a static wind to a MARCS atmospheric model can accurately reproduce observed extensions in the atmospheres of RSGs. In this work, we compute synthetic observables that match new interferometric data of the RSGs AH Sco, KW Sgr, V602 Car, CK Car, and V460 Car obtained with VLTI/MATISSE and GRAVITY. We computed model spectra and visibilities and found the best-fit model, mass-loss rate, and best-fit angular Rosseland diameter for the observations. We matched our model to the data, covering a wavelength range of $1.8-5.0\,μ$m, which corresponds to the $K$, $L,$ and $M$ bands. Our models reproduce the spectro-interferometric data over this wide wavelength range, including extended atmospheric layers of CO, H$_2$O, and SiO. We obtain Rosseland angular diameters between $3.0<θ_{\mathrm{Ross}}<5.5$ mas and mass-loss rates of $-6.5<\log \dot{M}/M_{\odot}\mathrm{yr}^{-1}<-4$. The partial pressure of SiO relative to the gas pressure and the SiO 4.0$\,μ$m line intensity increase between 2 and 3 stellar radii. The relative intensity depends on the luminosity used for our models, since the more luminous models have a higher mass-loss rate. This work further demonstrates that our MARCS+wind model can reproduce the observed physical extension of RSG atmospheres for several spectral diagnostics spanning a broad wavelength range. We reproduce both spectra and visibilities of newly obtained data as well as provide temperature and density stratifications that are consistent with the observations. With the MATISSE data, we newly include the extension of SiO layers as a precursor of silicate dust. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: Accepted for publication in A&A

Journal ref: 2024A&A...683A..19G

The EBLM Project XI. Mass, radius and effective temperature measurements for 23 M-dwarf companions to solar-type stars observed with CHEOPS

Authors: M. I. Swayne, P. F. L. Maxted, A. H. M. J. Triaud, S. G. Sousa, A. Deline, D. Ehrenreich, S. Hoyer, G. Olofsson, I. Boisse, A. Duck, S. Gill, D. Martin, J. McCormac, C. M. Persson, A. Santerne, D. Sebastian, M. R. Standing, L. Acuña, Y. Alibert, R. Alonso, G. Anglada, T. Bárczy, D. Barrado Navascues, S. C. C. Barros, W. Baumjohann , et al. (82 additional authors not shown)

Abstract: Observations of low-mass stars have frequently shown a disagreement between observed stellar radii and radii predicted by theoretical stellar structure models. This ``radius inflation'' problem could have an impact on both stellar and exoplanetary science. We present the final results of our observation programme with the CHEOPS satellite to obtain high-precision light curves of eclipsing binaries… ▽ More Observations of low-mass stars have frequently shown a disagreement between observed stellar radii and radii predicted by theoretical stellar structure models. This ``radius inflation'' problem could have an impact on both stellar and exoplanetary science. We present the final results of our observation programme with the CHEOPS satellite to obtain high-precision light curves of eclipsing binaries with low mass stellar companions (EBLMs). Combined with the spectroscopic orbits of the solar-type companion, we can derive the masses, radii and effective temperatures of 23 M-dwarf stars. We use the PYCHEOPS data analysis software to analyse their primary and secondary occultations. For all but one target, we also perform analyses with TESS light curves for comparison. We have assessed the impact of starspot-induced variation on our derived parameters and account for this in our radius and effective temperature uncertainties using simulated light curves. We observe trends for inflation with both metallicity and orbital separation. We also observe a strong trend in the difference between theoretical and observational effective temperatures with metallicity. There is no such trend with orbital separation. These results are not consistent with the idea that observed inflation in stellar radius combines with lower effective temperature to preserve the luminosity predicted by low-mass stellar models. Our EBLM systems are high-quality and homogeneous measurements that can be used in further studies into radius inflation. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 21 pages, 10 figures, accepted for publication in MNRAS, Supplementary material provided as ancillary files

Constraints on the Evolution of the Ionizing Background and Ionizing Photon Mean Free Path at the End of Reionization

Authors: Frederick B. Davies, Sarah E. I. Bosman, Prakash Gaikwad, Fahad Nasir, Joseph F. Hennawi, George D. Becker, Martin G. Haehnelt, Valentina D'Odorico, Manuela Bischetti, Anna-Christina Eilers, Laura C. Keating, Girish Kulkarni, Samuel Lai, Chiara Mazzucchelli, Yuxiang Qin, Sindhu Satyavolu, Feige Wang, Jinyi Yang, Yongda Zhu

Abstract: The variations in Ly$α$ forest opacity observed at $z>5.3$ between lines of sight to different background quasars are too strong to be caused by fluctuations in the density field alone. The leading hypothesis for the cause of this excess variance is a late, ongoing reionization process at redshifts below six. Another model proposes strong ionizing background fluctuations coupled to a short, spatia… ▽ More The variations in Ly$α$ forest opacity observed at $z>5.3$ between lines of sight to different background quasars are too strong to be caused by fluctuations in the density field alone. The leading hypothesis for the cause of this excess variance is a late, ongoing reionization process at redshifts below six. Another model proposes strong ionizing background fluctuations coupled to a short, spatially varying mean free path of ionizing photons, without explicitly invoking incomplete reionization. With recent observations suggesting a short mean free path at $z\sim6$, and a dramatic improvement in $z>5$ Ly$α$ forest data quality, we revisit this latter possibility. Here we apply the likelihood-free inference technique of approximate Bayesian computation to jointly constrain the hydrogen photoionization rate $Γ_{\rm HI}$ and the mean free path of ionizing photons $λ_{\rm mfp}$ from the effective optical depth distributions at $z=5.0$-$6.1$ from XQR-30. We find that the observations are well-described by fluctuating mean free path models with average mean free paths that are consistent with the steep trend implied by independent measurements at $z\sim5$-$6$, with a concomitant rapid evolution of the photoionization rate. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 23 pages, 13 figures, resubmitted to ApJ after referee's comments

Signatures of Small-scale Structure of the Pre-reionization Intergalactic Medium in $z\gtrsim7$ Quasar Proximity Zones

Authors: Frederick B. Davies, Joseph F. Hennawi

Abstract: The small-scale structure of baryons in the intergalactic medium is intimately linked to their past thermal history. Prior to the $\gtrsim10^4$ K photoheating during the epoch of reionization, cold baryons may have closely traced the clumpy cosmic web of dark matter down to scales as low as $\lesssim1$ comoving kpc, depending on the degree of heating by the X-ray background. After the passage of t… ▽ More The small-scale structure of baryons in the intergalactic medium is intimately linked to their past thermal history. Prior to the $\gtrsim10^4$ K photoheating during the epoch of reionization, cold baryons may have closely traced the clumpy cosmic web of dark matter down to scales as low as $\lesssim1$ comoving kpc, depending on the degree of heating by the X-ray background. After the passage of the ionization front, this clumpy structure can persist for $\sim10^{8}$ years. The strong Ly$α$ damping wings detected towards a few of the highest redshift quasars, in addition to their smaller-than-expected Ly$α$-transmissive proximity zones, suggest that they have ionized and heated the foreground intergalactic medium less than $10^7$ years ago. Signatures of the pre-reionization small-scale structure should thus persist in their intergalactic surroundings. Here we explore how the persistence of this clumpy structure can affect the statistics of Ly$α$ transmission inside the transparent proximity zones of $z\gtrsim7$ quasars by post-processing a suite of small-volume hydrodynamical simulations with 1D ionizing radiative transfer. We find that the Ly$α$ flux power spectrum and flux PDF statistics of ten $z=7.5$ proximity zones, with realistic observational parameters, could distinguish the gaseous structure of a $T_{\rm IGM}\sim2$ K CDM model from warm dark matter models with particle masses $m_{\rm WDM}>10$ keV and X-ray heated models with $f_{\rm X}f_{\rm abs}>0.1$ ($T_{\rm IGM}(z=7.5)\gtrsim275$ K) at the $2σ$ level. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 18 pages, 17 figures, submitted to ApJ

Possible Contamination of the Intergalactic Medium Damping Wing in ULAS J1342+0928 by Proximate Damped Ly$α$ Absorption

Authors: Frederick B. Davies, Eduardo Bañados, Joseph F. Hennawi, Sarah E. I. Bosman

Abstract: The red damping wing from neutral hydrogen in the intergalactic medium is a smoking-gun signal of ongoing reionization. One potential contaminant of the intergalactic damping wing signal is dense gas associated with foreground galaxies, which can give rise to proximate damped Ly$α$ absorbers. The Ly$α$ imprint of such absorbers on background quasars is indistinguishable from the intergalactic medi… ▽ More The red damping wing from neutral hydrogen in the intergalactic medium is a smoking-gun signal of ongoing reionization. One potential contaminant of the intergalactic damping wing signal is dense gas associated with foreground galaxies, which can give rise to proximate damped Ly$α$ absorbers. The Ly$α$ imprint of such absorbers on background quasars is indistinguishable from the intergalactic medium within the uncertainty of the intrinsic quasar continuum, and their abundance at $z\gtrsim7$ is unknown. Here we show that the complex of low-ionization metal absorption systems recently discovered by deep JWST/NIRSpec observations in the foreground of the $z=7.54$ quasar ULAS~J1342$+$0928 can potentially reproduce the quasar's spectral profile close to rest-frame Ly$α$ without invoking a substantial contribution from the intergalactic medium, but only if the absorbing gas is extremely metal-poor ($[{\rm O}/{\rm H}]\sim-3.5$). Such a low oxygen abundance has never been observed in a damped Ly$α$ absorber at any redshift, but this possibility still complicates the interpretation of the spectrum. Our analysis highlights the need for deep spectroscopy of high-redshift quasars with JWST or ELT to "purify" damping wing quasar samples, an exercise which is impossible for much fainter objects like galaxies. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 10 pages, 5 figures, submitted to ApJL

The impact of varying inhomogeneous reionization histories on metrics of Ly$α$ opacity

Authors: Caitlin C. Doughty, Joseph F. Hennawi, Jose Oñorbe, Frederick B. Davies, Zarija Lukić

Abstract: The epoch of hydrogen reionization is complete by $z=5$, but its progression at higher redshifts is uncertain. Measurements of Ly$α$ forest opacity show large scatter at $z<6$, suggestive of spatial fluctuations in neutral fraction ($x_\mathrm{HI}$), temperature, or ionizing background, either individually or in combination. However, these effects are degenerate, necessitating modeling these physi… ▽ More The epoch of hydrogen reionization is complete by $z=5$, but its progression at higher redshifts is uncertain. Measurements of Ly$α$ forest opacity show large scatter at $z<6$, suggestive of spatial fluctuations in neutral fraction ($x_\mathrm{HI}$), temperature, or ionizing background, either individually or in combination. However, these effects are degenerate, necessitating modeling these physics in tandem in order to properly interpret the observations. We begin this process by developing a framework for modeling the reionization history and associated temperature fluctuations, with the intention of incorporating ionizing background fluctuations at a later time. To do this, we generate several reionization histories using semi-numerical code AMBER, selecting histories with volume-weighted neutral fractions that adhere to the observed CMB optical depth and dark pixel fractions. Implementing these histories in the \texttt{Nyx} cosmological hydrodynamics code, we examine the evolution of gas within the simulation, and the associated metrics of the Ly$α$ forest opacity. We find that the pressure smoothing scale within the IGM is strongly correlated with the adiabatic index of the temperature-density relation. We find that while models with 20,000 K photoheating at reionization are better able to reproduce the shape of the observed $z=5$ 1D flux power spectrum than those with 10,000 K, they fail to match the highest wavenumbers. The simulated autocorrelation function and optical depth distributions are systematically low and narrow, respectively, compared to the observed values, but are in better agreement when the reionization history is longer in duration, more symmetric in its distribution of reionization redshifts, or if there are remaining neutral regions at $z<6$. The systematically low variance likely requires the addition of a fluctuating UVB. △ Less

Submitted 3 December, 2023; originally announced December 2023.

Comments: To be submitted to MNRAS. Comments welcome

A resonant sextuplet of sub-Neptunes transiting the bright star HD 110067

Authors: R. Luque, H. P. Osborn, A. Leleu, E. Pallé, A. Bonfanti, O. Barragán, T. G. Wilson, C. Broeg, A. Collier Cameron, M. Lendl, P. F. L. Maxted, Y. Alibert, D. Gandolfi, J. -B. Delisle, M. J. Hooton, J. A. Egger, G. Nowak, M. Lafarga, D. Rapetti, J. D. Twicken, J. C. Morales, I. Carleo, J. Orell-Miquel, V. Adibekyan, R. Alonso , et al. (127 additional authors not shown)

Abstract: Planets with radii between that of the Earth and Neptune (hereafter referred to as sub-Neptunes) are found in close-in orbits around more than half of all Sun-like stars. Yet, their composition, formation, and evolution remain poorly understood. The study of multi-planetary systems offers an opportunity to investigate the outcomes of planet formation and evolution while controlling for initial con… ▽ More Planets with radii between that of the Earth and Neptune (hereafter referred to as sub-Neptunes) are found in close-in orbits around more than half of all Sun-like stars. Yet, their composition, formation, and evolution remain poorly understood. The study of multi-planetary systems offers an opportunity to investigate the outcomes of planet formation and evolution while controlling for initial conditions and environment. Those in resonance (with their orbital periods related by a ratio of small integers) are particularly valuable because they imply a system architecture practically unchanged since its birth. Here, we present the observations of six transiting planets around the bright nearby star HD 110067. We find that the planets follow a chain of resonant orbits. A dynamical study of the innermost planet triplet allowed the prediction and later confirmation of the orbits of the rest of the planets in the system. The six planets are found to be sub-Neptunes with radii ranging from 1.94 to 2.85 Re. Three of the planets have measured masses, yielding low bulk densities that suggest the presence of large hydrogen-dominated atmospheres. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: Published in Nature on November 30, 2023. Supplementary Information can be found in the online version of the paper in the journal

Journal ref: Nature 623, 932-937 (2023)

Characterising TOI-732 b and c: New insights into the M-dwarf radius and density valley

Authors: A. Bonfanti, M. Brady, T. G. Wilson, J. Venturini, J. A. Egger, A. Brandeker, S. G. Sousa, M. Lendl, A. E. Simon, D. Queloz, G. Olofsson, V. Adibekyan, Y. Alibert, L. Fossati, M. J. Hooton, D. Kubyshkina, R. Luque, F. Murgas, A. J. Mustill, N. C. Santos, V. Van Grootel, R. Alonso, J. Asquier, T. Bandy, T. Bárczy , et al. (66 additional authors not shown)

Abstract: TOI-732 is an M dwarf hosting two transiting planets that are located on the two opposite sides of the radius valley. By doubling the number of available space-based observations and increasing the number of radial velocity (RV) measurements, we aim at refining the parameters of TOI-732 b and c. We also use the results to study the slope of the radius valley and the density valley for a well-chara… ▽ More TOI-732 is an M dwarf hosting two transiting planets that are located on the two opposite sides of the radius valley. By doubling the number of available space-based observations and increasing the number of radial velocity (RV) measurements, we aim at refining the parameters of TOI-732 b and c. We also use the results to study the slope of the radius valley and the density valley for a well-characterised sample of M-dwarf exoplanets. We performed a global MCMC analysis by jointly modelling ground-based light curves and CHEOPS and TESS observations, along with RV time series both taken from the literature and obtained with the MAROON-X spectrograph. The slopes of the M-dwarf valleys were quantified via a Support Vector Machine (SVM) procedure. TOI-732 b is an ultrashort-period planet ($P\sim0.77$ d) with a radius $R_b=1.325_{-0.058}^{+0.057}$ $R_{\oplus}$ and a mass $M_b=2.46\pm0.19$ $M_{\oplus}$ (mean density $ρ_b=5.8_{-0.8}^{+1.0}$ g cm$^{-3}$), while the outer planet at $P\sim12.25$ d has $R_c=2.39_{-0.11}^{+0.10}$ $R_{\oplus}$, $M_c=8.04_{-0.48}^{+0.50}$ $M_{\oplus}$, and thus $ρ_c=3.24_{-0.43}^{+0.55}$ g cm$^{-3}$. Also taking into account our interior structure calculations, TOI-732 b is a super-Earth and TOI-732 c is a mini-Neptune. Following the SVM approach, we quantified $\mathrm{d}\log{R_{p,{\mathrm{valley}}}}/\mathrm{d}\log{P}=-0.065_{-0.013}^{+0.024}$, which is flatter than for Sun-like stars. In line with former analyses, we note that the radius valley for M-dwarf planets is more densely populated, and we further quantify the slope of the density valley as $\mathrm{d}\log{\hatρ_{\mathrm{valley}}}/\mathrm{d}\log{P}=-0.02_{-0.04}^{+0.12}$. Compared to FGK stars, the weaker dependence of the position of the radius valley on the orbital period might indicate that the formation shapes the radius valley around M dwarfs more strongly than the evolution mechanisms. △ Less

Submitted 30 November, 2023; v1 submitted 21 November, 2023; originally announced November 2023.

Comments: 28 pages (17 in the main text), 18 figures (9 in the main text), 11 tables (7 in the main text). Accepted for publication in A&A

CHEOPS observations of KELT-20 b/MASCARA-2 b: An aligned orbit and signs of variability from a reflective dayside

Authors: V. Singh, G. Scandariato, A. M. S. Smith, P. E. Cubillos, M. Lendl, N. Billot, A. Fortier, D. Queloz, S. G. Sousa, Sz. Csizmadia, A. Brandeker, L. Carone, T. G. Wilson, B. Akinsanmi, J. A. Patel, A. Krenn, O. D. S. Demangeon, G. Bruno, I. Pagano, M. J. Hooton, J. Cabrera, N. C. Santos, Y. Alibert, R. Alonso, J. Asquier , et al. (65 additional authors not shown)

Abstract: Occultations are windows of opportunity to indirectly peek into the dayside atmosphere of exoplanets. High-precision transit events provide information on the spin-orbit alignment of exoplanets around fast-rotating hosts. We aim to precisely measure the planetary radius and geometric albedo of the ultra-hot Jupiter (UHJ) KELT-20 b as well as the system's spin-orbit alignment. We obtained optical h… ▽ More Occultations are windows of opportunity to indirectly peek into the dayside atmosphere of exoplanets. High-precision transit events provide information on the spin-orbit alignment of exoplanets around fast-rotating hosts. We aim to precisely measure the planetary radius and geometric albedo of the ultra-hot Jupiter (UHJ) KELT-20 b as well as the system's spin-orbit alignment. We obtained optical high-precision transits and occultations of KELT-20 b using CHEOPS observations in conjunction with the simultaneous TESS observations. We interpreted the occultation measurements together with archival infrared observations to measure the planetary geometric albedo and dayside temperatures. We further used the host star's gravity-darkened nature to measure the system's obliquity. We present a time-averaged precise occultation depth of 82(6) ppm measured with seven CHEOPS visits and 131(+8/-7) ppm from the analysis of all available TESS photometry. Using these measurements, we precisely constrain the geometric albedo of KELT-20 b to 0.26(0.04) and the brightness temperature of the dayside hemisphere to 2566(+77/-80) K. Assuming Lambertian scattering law, we constrain the Bond albedo to 0.36(+0.04/-0.05) along with a minimal heat transfer to the night side. Furthermore, using five transit observations we provide stricter constraints of 3.9(1.1) degrees on the sky-projected obliquity of the system. The aligned orbit of KELT-20 b is in contrast to previous CHEOPS studies that have found strongly inclined orbits for planets orbiting other A-type stars. The comparably high planetary geometric albedo of KELT-20 b corroborates a known trend of strongly irradiated planets being more reflective. Finally, we tentatively detect signs of temporal variability in the occultation depths, which might indicate variable cloud cover advecting onto the planetary day side. △ Less

Submitted 29 November, 2023; v1 submitted 6 November, 2023; originally announced November 2023.

Comments: 27 pages, 15 figures, Accepted for publication in Astronomy & Astrophysics

No random transits in CHEOPS observations of HD 139139

Authors: R. Alonso, S. Hoyer, M. Deleuil, A. E. Simon, M. Beck, W. Benz, H. -G. Florén, P. Guterman, L. Borsato, A. Brandeker, D. Gandolfi, T. G. Wilson, T. Zingales, Y. Alibert, G. Anglada, T. Bárczy, D. Barrado Navascues, S. C. C. Barros, W. Baumjohann, T. Beck, N. Billot, X. Bonfils, Ch. Broeg, S. Charnoz, A. Collier Cameron , et al. (56 additional authors not shown)

Abstract: HD 139139 (a.k.a. 'The Random Transiter') is a star that exhibited enigmatic transit-like features with no apparent periodicity in K2 data. The shallow depth of the events ($\sim$200 ppm -- equivalent to transiting objects with radii of $\sim$1.5 R$_\oplus$ in front of a Sun-like star), and their non-periodicity, constitutes a challenge for the photometric follow-up of this star. The goal of this… ▽ More HD 139139 (a.k.a. 'The Random Transiter') is a star that exhibited enigmatic transit-like features with no apparent periodicity in K2 data. The shallow depth of the events ($\sim$200 ppm -- equivalent to transiting objects with radii of $\sim$1.5 R$_\oplus$ in front of a Sun-like star), and their non-periodicity, constitutes a challenge for the photometric follow-up of this star. The goal of this study is to confirm with independent measurements the presence of shallow, non-periodic transit-like features on this object. We performed observations with CHEOPS, for a total accumulated time of 12.75 d, distributed in visits of roughly 20 h in two observing campaigns in years 2021 and 2022. The precision of the data is sufficient to detect 150 ppm features with durations longer than 1.5 h. We use the duration and times of the events seen in the K2 curve to estimate how many should have been detected in our campaigns, under the assumption that their behaviour during the CHEOPS observations would be the same as in the K2 data of 2017. We do not detect events with depths larger than 150 ppm in our data set. If the frequency, depth, and duration of the events were the same as in the K2 campaign, we estimate the probability of having missed all events due to our limited observing window would be 4.8 %. We suggest three different scenarios to explain our results: 1) Our observing window was not long enough, and the events were missed with the estimated 4.8 % probability. 2) The events recorded in the K2 observations were time critical, and the mechanism producing them was either not active in the 2021 and 2022 campaigns or created shallower events under our detectability level. 3) The enigmatic events in the K2 data are the result of an unidentified and infrequent instrumental noise in the original data set or its data treatment. △ Less

Submitted 25 October, 2023; v1 submitted 16 October, 2023; originally announced October 2023.

Comments: 7 pages, 3 figures, accepted for publication in A&A. Language-corrected version

Journal ref: A&A 680, A78 (2023)

JWST reveals a luminous infrared source at the position of the failed supernova candidate N6946-BH1

Authors: Emma R. Beasor, Griffin Hosseinzadeh, Nathan Smith, Ben Davies, Jacob E. Jencson, Jeniveve Pearson, David J. Sand

Abstract: N6946-BH1 is the first plausible candidate for a failed supernova (SN), a peculiar event in which a massive star disappears without the expected bright SN, accompanied by collapse into a black hole (BH). Following a luminous outburst in 2009, the source experienced a significant decline in optical brightness, while maintaining a persistent infrared (IR) presence. While it was proposed to be a pote… ▽ More N6946-BH1 is the first plausible candidate for a failed supernova (SN), a peculiar event in which a massive star disappears without the expected bright SN, accompanied by collapse into a black hole (BH). Following a luminous outburst in 2009, the source experienced a significant decline in optical brightness, while maintaining a persistent infrared (IR) presence. While it was proposed to be a potential failed SN, such behavior has been observed in SN impostor events in nearby galaxies. Here, we present late-time observations of BH1, taken 14 years after disappearance, using JWST's NIRCam and MIRI instruments to probe a never-before-observed region of the object's spectral energy distribution. We show for the first time that all previous observations of BH1 (pre- and post-disappearance) are actually a blend of at least 3 sources. In the near-IR, BH1 is notably fainter than the progenitor but retains similar brightness to its state in 2017. In the mid-IR, the flux appears to have brightened compared to the inferred fluxes from the best-fitting progenitor model. The total luminosity of the source is between 13 - 25% that of the progenitor. We also show that the IR SED appears consistent with PAH features that arise when dust is illuminated by near-ultraviolet radiation. At present, the interpretation of N6946-BH1 remains uncertain. The observations match expectations for a stellar merger, but theoretical ambiguity in the failed SN hypothesis makes it hard to dismiss. △ Less

Submitted 22 January, 2024; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: accepted to ApJ

Constraining the reflective properties of WASP-178b using Cheops photometry

Authors: I. Pagano, G. Scandariato, V. Singh, M. Lendl, D. Queloz, A. E. Simon, S. G. Sousa, A. Brandeker, A. Collier Cameron, S. Sulis, V. Van Grootel, T. G. Wilson, Y. Alibert, R. Alonso, G. Anglada, T. Bárczy, D. Barrado Navascues, S. C. C. Barros, W. Baumjohann, M. Beck, T. Beck, W. Benz, N. Billot, X. Bonfils, L. Borsato , et al. (57 additional authors not shown)

Abstract: Multiwavelength photometry of the secondary eclipses of extrasolar planets is able to disentangle the reflected and thermally emitted light radiated from the planetary dayside. This leads to the measurement of the planetary geometric albedo $A_g$, which is an indicator of the presence of clouds in the atmosphere, and the recirculation efficiency $ε$, which quantifies the energy transport within th… ▽ More Multiwavelength photometry of the secondary eclipses of extrasolar planets is able to disentangle the reflected and thermally emitted light radiated from the planetary dayside. This leads to the measurement of the planetary geometric albedo $A_g$, which is an indicator of the presence of clouds in the atmosphere, and the recirculation efficiency $ε$, which quantifies the energy transport within the atmosphere. In this work we aim to measure $A_g$ and $ε$ for the planet WASP-178 b, a highly irradiated giant planet with an estimated equilibrium temperature of 2450 K.} We analyzed archival spectra and the light curves collected by Cheops and Tess to characterize the host WASP-178, refine the ephemeris of the system and measure the eclipse depth in the passbands of the two respective telescopes. We measured a marginally significant eclipse depth of 70$\pm$40 ppm in the Tess passband and statistically significant depth of 70$\pm$20 ppm in the Cheops passband. Combining the eclipse depth measurement in the Cheops (lambda_eff=6300 AA) and Tess (lambda_eff=8000 AA) passbands we constrained the dayside brightness temperature of WASP-178 b in the 2250-2800 K interval. The geometric albedo 0.1<$\rm A_g$<0.35 is in general agreement with the picture of poorly reflective giant planets, while the recirculation efficiency $ε>$0.7 makes WASP-178 b an interesting laboratory to test the current heat recirculation models. △ Less

Submitted 16 September, 2023; originally announced September 2023.

Comments: Accepted by Astronomy and Astrophysics on 31/08/2023

Forecasting constraints on the high-z IGM thermal state from the Lyman-$α$ forest flux auto-correlation function

Authors: Molly Wolfson, Joseph F. Hennawi, Frederick B. Davies, Zarija Lukić, Jose Oñorbe

Abstract: The auto-correlation function of the Lyman-$α$ (Ly$α$) forest flux from high-z quasars can statistically probe all scales of the intergalactic medium (IGM) just after the epoch of reionization. The thermal state of the IGM, which is determined by the physics of reionization, sets the amount of small-scale power seen in the \lya forest. To study the sensitivity of the auto-correlation function to t… ▽ More The auto-correlation function of the Lyman-$α$ (Ly$α$) forest flux from high-z quasars can statistically probe all scales of the intergalactic medium (IGM) just after the epoch of reionization. The thermal state of the IGM, which is determined by the physics of reionization, sets the amount of small-scale power seen in the \lya forest. To study the sensitivity of the auto-correlation function to the thermal state of the IGM, we compute the auto-correlation function from cosmological hydrodynamical simulations with semi-numerical models of the thermal state of the IGM. We create mock data sets of 20 quasars to forecast constraints on $T_0$ and $γ$, which characterize a tight temperature-density relation in the IGM, at $5.4 \leq z \leq 6$. At $z = 5.4$ we find that an ideal data set constrains $T_0$ to 29\% and $γ$ to 9\%. In addition, we investigate four realistic reionization scenarios that combine temperature and ultra-violet background (UVB) fluctuations at $z = 5.8$. We find that, when using mock data generated from a model that includes temperature and UVB fluctuations, we can rule out a model with no temperature or UVB fluctuations at $>1σ$ level 50.5\% of the time. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 18 pages, 13 figures, comments welcome. arXiv admin note: text overlap with arXiv:2208.09013

Measurements of the $z > 5$ Lyman-$α$ forest flux auto-correlation functions from the extended XQR-30 data set

Authors: Molly Wolfson, Joseph F. Hennawi, Sarah E. I. Bosman, Frederick B. Davies, Zarija Lukić, George D. Becker, Huanqing Chen, Guido Cupani, Valentina D'Odorico, Anna-Christina Eilers, Martin G. Haehnelt, Laura C. Keating, Girish Kulkarni, Samuel Lai, Andrei Mesinger, Fabian Walter, Yongda Zhu

Abstract: Recently, the Lyman-$α$ (Ly$α$) forest flux auto-correlation function has been shown to be sensitive to the mean free path of hydrogen-ionizing photons, $λ_{\text{mfp}}$, for simulations at $z \geq 5.4$. Measuring $λ_{\text{mfp}}$ at these redshifts will give vital information on the ending of reionization. Here we present the first observational measurements of the Ly$α$ forest flux auto-correlat… ▽ More Recently, the Lyman-$α$ (Ly$α$) forest flux auto-correlation function has been shown to be sensitive to the mean free path of hydrogen-ionizing photons, $λ_{\text{mfp}}$, for simulations at $z \geq 5.4$. Measuring $λ_{\text{mfp}}$ at these redshifts will give vital information on the ending of reionization. Here we present the first observational measurements of the Ly$α$ forest flux auto-correlation functions in ten redshift bins from $5.1 \leq z \leq 6.0$. We use a sample of 35 quasar sightlines at $z > 5.7$ from the extended XQR-30 data set, this data has signal-to-noise ratios of $> 20$ per spectral pixel. We carefully account for systematic errors in continuum reconstruction, instrumentation, and contamination by damped Ly$α$ systems. With these measurements, we introduce software tools to generate auto-correlation function measurements from any simulation. For an initial comparison, we show our auto-correlation measurements with simulation models for recently measured $λ_{\text{mfp}}$ values and find good agreements. Further work in modeling and understanding the covariance matrices of the data is necessary to get robust measurements of $λ_{\text{mfp}}$ from this data. △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: 17 pages, 13 figures, submitted to MNRAS

The relationship between IGM Lyman-alpha opacity and galaxy density near the end of reionization

Authors: Holly M. Christenson, George D. Becker, Anson D'Aloisio, Frederick B. Davies, Yongda Zhu, Elisa Boera, Fahad Nasir, Steven R. Furlanetto, Matthew A. Malkan

Abstract: Observed scatter in the Lyman-alpha opacity of quasar sightlines at $z<6$ has motivated measurements of the correlation between Ly$α$ opacity and galaxy density, as models that predict this scatter make strong and sometimes opposite predictions for how they should be related. Our previous work associated two highly opaque Ly$α$ troughs at $z\sim5.7$ with a deficit of Lyman-$α$ emitting galaxies (L… ▽ More Observed scatter in the Lyman-alpha opacity of quasar sightlines at $z<6$ has motivated measurements of the correlation between Ly$α$ opacity and galaxy density, as models that predict this scatter make strong and sometimes opposite predictions for how they should be related. Our previous work associated two highly opaque Ly$α$ troughs at $z\sim5.7$ with a deficit of Lyman-$α$ emitting galaxies (LAEs). In this work, we survey two of the most highly transmissive lines of sight at this redshift, towards the $z=6.02$ quasar SDSS J1306+0356 and the $z=6.17$ quasar PSO J359-06. We find that both fields are underdense in LAEs within 10 $h^{-1}$ Mpc of the quasar sightline, somewhat less extensive than underdensities associated with Ly$α$ troughs. We combine our observations with three additional fields from the literature, and find that while fields with extreme opacities are generally underdense, moderate opacities span a wider density range. The results at high opacities are consistent with models that invoke UV background fluctuations and/or late reionization to explain the observed scatter in IGM Ly$α$ opacities. There is tension at low opacities, however, as the models tend to associate lower IGM Ly$α$ opacities with higher densities. Although the number of fields surveyed is still small, the low-opacity results may support a scenario in which the ionizing background in low-density regions increases more rapidly than some models suggest after becoming ionized. Elevated gas temperatures from recent reionization may also be making these regions more transparent. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Comments: 20 pages, 17 figures. Accepted to ApJ

First measurement of the Mg II forest correlation function in the Epoch of Reionization

Authors: Suk Sien Tie, Joseph F. Hennawi, Feige Wang, Silvia Onorato, Jinyi Yang, Eduardo Bañados, Frederick B. Davies, Jose Oñorbe

Abstract: In the process of producing the roughly three ionizing photons per atom required to reionize the IGM, the same massive stars explode and eject metals into their surroundings. While the overly sensitive Lya transition makes Gunn-Peterson absorption of background quasar light an ineffective probe of reionization at z > 6, strong low-ionization transitions like the MgII doublet will give rise to a de… ▽ More In the process of producing the roughly three ionizing photons per atom required to reionize the IGM, the same massive stars explode and eject metals into their surroundings. While the overly sensitive Lya transition makes Gunn-Peterson absorption of background quasar light an ineffective probe of reionization at z > 6, strong low-ionization transitions like the MgII doublet will give rise to a detectable "metal-line forest", if metals pollute the neutral IGM. We measure the auto-correlation of the MgII forest transmission using a sample of ten ground based z >= 6.80 quasar spectra probing the redshift range 5.96 < z_MgII < 7.42 (z_MgII,median = 6.47). The correlation function exhibits strong small-scale clustering and a pronounced peak at the doublet velocity (768 km/s) arising from strong absorbers in the CGM of galaxies. After these strong absorbers are identified and masked the signal is consistent with noise. Our measurements are compared to a suite of models generated by combining a large hydrodynamical simulation with a semi-numerical reionization topology, assuming a simple uniform enrichment model. We obtain a 95% credibility upper limit of [Mg/H] < -3.73 at z_MgII,median = 6.47, assuming uninformative priors on [Mg/H] and the IGM neutral fraction x_HI. Splitting the data into low-z (5.96 < z_MgII < 6.47; z_MgII,median = 6.235) and high-z (6.47 < z_MgII < 7.42; z_MgII,median = 6.72) subsamples again yields null-detections and 95% upper limits of [Mg/H] < -3.75 and [Mg/H] < -3.45, respectively. These first measurements set the stage for an approved JWST Cycle 2 program (GO 3526) targeting a similar number of quasars that will be an order of magnitude more sensitive, making the Mgii forest an emerging powerful tool to deliver precision constraints on the reionization and enrichment history of the Universe. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: Submitted to MNRAS; 23 pages, 21 figures, 1 table

Refining the properties of the TOI-178 system with CHEOPS and TESS

Authors: L. Delrez, A. Leleu, A. Brandeker, M. Gillon, M. J. Hooton, A. Collier Cameron, A. Deline, A. Fortier, D. Queloz, A. Bonfanti, V. Van Grootel, T. G. Wilson, J. A. Egger, Y. Alibert, R. Alonso, G. Anglada, J. Asquier, T. Bárczy, D. Barrado y Navascues, S. C. C. Barros, W. Baumjohann, M. Beck, T. Beck, W. Benz, N. Billot , et al. (62 additional authors not shown)

Abstract: The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. Mass estimates derived from a preliminary radial velocity (RV) dataset suggest that the planetary densities do not decrease in a monotonic way with the orbital d… ▽ More The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. Mass estimates derived from a preliminary radial velocity (RV) dataset suggest that the planetary densities do not decrease in a monotonic way with the orbital distance to the star, contrary to what one would expect based on simple formation and evolution models. To improve the characterisation of this key system and prepare for future studies (in particular with JWST), we perform a detailed photometric study based on 40 new CHEOPS visits, one new TESS sector, as well as previously published CHEOPS, TESS, and NGTS data. First we perform a global analysis of the 100 transits contained in our data to refine the transit parameters of the six planets and study their transit timing variations (TTVs). We then use our extensive dataset to place constraints on the radii and orbital periods of potential additional transiting planets in the system. Our analysis significantly refines the transit parameters of the six planets, most notably their radii, for which we now obtain relative precisions $\lesssim$3%, with the exception of the smallest planet $b$ for which the precision is 5.1%. Combined with the RV mass estimates, the measured TTVs allow us to constrain the eccentricities of planets $c$ to $g$, which are found to be all below 0.02, as expected from stability requirements. Taken alone, the TTVs also suggest a higher mass for planet $d$ than the one estimated from the RVs, which had been found to yield a surprisingly low density for this planet. However, the masses derived from the current TTV dataset are very prior-dependent and further observations, over a longer temporal baseline, are needed to deepen our understanding of this iconic planetary system. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 20 pages, 5 figures, 9 tables. Accepted for publication in A&A

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

CHEOPS and TESS view of the ultra-short period super-Earth TOI-561 b

Authors: J. A. Patel, J. A. Egger, T. G. Wilson, V. Bourrier, L. Carone, M. Beck, D. Ehrenreich, S. G. Sousa, W. Benz, A. Brandeker, A. Deline, Y. Alibert, K. W. F. Lam, M. Lendl, R. Alonso, G. Anglada, T. Bárczy, D. Barrado, S. C. C. Barros, W. Baumjohann, T. Beck, N. Billot, X. Bonfils, C. Broeg, M. -D. Busch , et al. (53 additional authors not shown)

Abstract: Ultra-short period planets (USPs) are a unique class of super-Earths with an orbital period of less than a day and hence subject to intense radiation from their host star. While most of them are consistent with bare rocks, some show evidence of a heavyweight envelope, which could be a water layer or a secondary metal-rich atmosphere sustained by an outgassing surface. Much remains to be learned ab… ▽ More Ultra-short period planets (USPs) are a unique class of super-Earths with an orbital period of less than a day and hence subject to intense radiation from their host star. While most of them are consistent with bare rocks, some show evidence of a heavyweight envelope, which could be a water layer or a secondary metal-rich atmosphere sustained by an outgassing surface. Much remains to be learned about the nature of USPs. The prime goal of the present work is to study the bulk planetary properties and atmosphere of TOI-561b, through the study of its transits and occultations. We obtained ultra-precise transit photometry of TOI-561b with CHEOPS and performed a joint analysis of this data with four TESS sectors. Our analysis of TOI-561b transit photometry put strong constraints on its properties, especially on its radius, Rp=1.42 +/- 0.02 R_Earth (at ~2% error). The internal structure modelling of the planet shows that the observations are consistent with negligible H/He atmosphere, however requiring other lighter materials, in addition to pure iron core and silicate mantle to explain the observed density. We find that this can be explained by the inclusion of a water layer in our model. We searched for variability in the measured Rp/R* over time to trace changes in the structure of the planetary envelope but none found within the data precision. In addition to the transit event, we tentatively detect occultation signal in the TESS data with an eclipse depth of ~27 +/- 11 ppm. Using the models of outgassed atmospheres from the literature we find that the thermal emission from the planet can mostly explain the observation. Based on this, we predict that NIR/MIR observations with JWST should be able to detect silicate species in the atmosphere of the planet. This could also reveal important clues about the planetary interior and help disentangle planet formation and evolution models. △ Less

Submitted 16 August, 2023; originally announced August 2023.

Comments: 17 pages, 10 + 3 figures, 4 tables, accepted for publication in A&A (abstract abbreviated)

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

Probing Ultra-late Reionization: Direct Measurements of the Mean Free Path over $5<z<6$

Authors: Yongda Zhu, George D. Becker, Holly M. Christenson, Anson D'Aloisio, Sarah E. I. Bosman, Tom Bakx, Valentina D'Odorico, Manuela Bischetti, Christopher Cain, Frederick B. Davies, Rebecca L. Davies, Anna-Christina Eilers, Xiaohui Fan, Prakash Gaikwad, Martin G. Haehnelt, Laura C. Keating, Girish Kulkarni, Samuel Lai, Hai-Xia Ma, Andrei Mesinger, Yuxiang Qin, Sindhu Satyavolu, Tsutomu T. Takeuchi, Hideki Umehata, Jinyi Yang

Abstract: The mean free path of ionizing photons, $λ_{\rm mfp}$, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements of $λ_{\rm mfp}$ from QSO spectra over the redshift range $5<z<6$, including the first measurements at $z\simeq5.3$ and 5.6. Our sample includes data from the XQR-30 VLT large program, as well as new Keck/ESI o… ▽ More The mean free path of ionizing photons, $λ_{\rm mfp}$, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements of $λ_{\rm mfp}$ from QSO spectra over the redshift range $5<z<6$, including the first measurements at $z\simeq5.3$ and 5.6. Our sample includes data from the XQR-30 VLT large program, as well as new Keck/ESI observations of QSOs near $z \sim 5.5$, for which we also acquire new [C II] 158$μ$m redshifts with ALMA. By measuring the Lyman continuum transmission profile in stacked QSO spectra, we find $λ_{\rm mfp} = 9.33_{-1.80}^{+2.06}$, $5.40_{-1.40}^{+1.47}$, $3.31_{-1.34}^{+2.74}$, and $0.81_{-0.48}^{+0.73}$ pMpc at $z=5.08$, 5.31, 5.65, and 5.93, respectively. Our results demonstrate that $λ_{\rm mfp}$ increases steadily and rapidly with time over $5<z<6$. Notably, we find that $λ_{\rm mfp}$ deviates significantly from predictions based on a fully ionized and relaxed IGM as late as $z=5.3$. By comparing our results to model predictions and indirect $λ_{\rm mfp}$ constraints based on IGM Ly$α$ opacity, we find that the $λ_{\rm mfp}$ evolution is consistent with scenarios wherein the IGM is still undergoing reionization and/or retains large fluctuations in the ionizing UV background well below redshift six. △ Less

Submitted 8 August, 2023; originally announced August 2023.

Comments: 17 pages, 6 figures, 2 tables; accepted for publication in ApJ

Investigating the visible phase-curve variability of 55 Cnc e

Authors: E. A. Meier Valdés, B. M. Morris, B. -O. Demory, A. Brandeker, D. Kitzmann, W. Benz, A. Deline, H. -G. Florén, S. G. Sousa, V. Bourrier, V. Singh, K. Heng, A. Strugarek, D. J. Bower, N. Jäggi, L. Carone, M. Lendl, K. Jones, A. V. Oza, O. D. S. Demangeon, Y. Alibert, R. Alonso, G. Anglada, J. Asquier, T. Bárczy , et al. (65 additional authors not shown)

Abstract: 55 Cnc e is an ultra-short period super-Earth transiting a Sun-like star. Previous observations in the optical range detected a time-variable flux modulation that is phased with the planetary orbital period, whose amplitude is too large to be explained by reflected light and thermal emission alone. The goal of the study is to investigate the origin of the variability and timescale of the phase-cur… ▽ More 55 Cnc e is an ultra-short period super-Earth transiting a Sun-like star. Previous observations in the optical range detected a time-variable flux modulation that is phased with the planetary orbital period, whose amplitude is too large to be explained by reflected light and thermal emission alone. The goal of the study is to investigate the origin of the variability and timescale of the phase-curve modulation in 55 Cnc e. To this end, we used the CHaracterising ExOPlanet Satellite (CHEOPS), whose exquisite photometric precision provides an opportunity to characterise minute changes in the phase curve from one orbit to the next. CHEOPS observed 29 individual visits of 55 Cnc e between March 2020 and February 2022. Based on these observations, we investigated the different processes that could be at the origin of the observed modulation. In particular, we built a toy model to assess whether a circumstellar torus of dust driven by radiation pressure and gravity might match the observed flux variability timescale. We find that the phase-curve amplitude and peak offset of 55 Cnc e do vary between visits. The sublimation timescales of selected dust species reveal that silicates expected in an Earth-like mantle would not survive long enough to explain the observed phase-curve modulation. We find that silicon carbide, quartz, and graphite are plausible candidates for the circumstellar torus composition because their sublimation timescales are long. The extensive CHEOPS observations confirm that the phase-curve amplitude and offset vary in time.We find that dust could provide the grey opacity source required to match the observations. However, the data at hand do not provide evidence that circumstellar material with a variable grain mass per unit area causes the observed variability. Future observations with the James Webb Space Telescope promise exciting insights into this iconic super-Earth. △ Less

Submitted 27 July, 2023; v1 submitted 12 July, 2023; originally announced July 2023.

Comments: 27 pages, 22 figures. Accepted for publication on A&A

Journal ref: A&A 677, A112 (2023)

Making hot Jupiters in stellar clusters II: efficient formation in binary systems

Authors: Daohai Li, Alexander J. Mustill, Melvyn B. Davies, Yan-Xiang Gong

Abstract: Observations suggested that the occurrence rate of hot Jupiters (HJs) in open clusters is largely consistent with the field ($\sim1\%$) but in the binary-rich cluster M67, the rate is $\sim5\%$. How does the cluster environment boost HJ formation via the high-eccentricity tidal migration initiated by the extreme-amplitude von Zeipel-Lidov-Kozai (XZKL) mechanism forced by a companion star? Our anal… ▽ More Observations suggested that the occurrence rate of hot Jupiters (HJs) in open clusters is largely consistent with the field ($\sim1\%$) but in the binary-rich cluster M67, the rate is $\sim5\%$. How does the cluster environment boost HJ formation via the high-eccentricity tidal migration initiated by the extreme-amplitude von Zeipel-Lidov-Kozai (XZKL) mechanism forced by a companion star? Our analytical treatment shows that the cluster's collective gravitational potential alters the companion's orbit slowly, which may render the star-planet-companion configuration XZKL-favourable, a phenomenon only possible for very wide binaries. We have also performed direct Gyr $N$-body simulations of the star cluster evolution and XZKL of planets' orbit around member stars. We find that an initially-single star may acquire a companion star via stellar scattering and the companion may enable XZKL in the planets' orbit. Planets around an initially-binary star may also be XZKL-activated by the companion. In both scenarios, the companion's orbit has likely been significantly changed by star scattering and the cluster potential before XZKL occurs in the planets' orbits. Across different cluster models, 0.8\%-3\% of the planets orbiting initially-single stars have experienced XZKL while the fraction is 2\%-26\% for initially-binary stars. Notably, the ejection fraction is similar to or appreciably smaller than XZKL. Around a star that is binary at 1 Gyr, 13\%-32\% of its planets have undergone XZKL, and combined with single stars, the overall XZKL fraction is 3\%-21\%, most affected by the cluster binarity. If 10\% of the stars in M67 host a giant planet, our model predicts an HJ occurrence rate of $\sim1\%$. We suggest that HJ surveys target old, high-binarity, not-too-dense open clusters and prioritise wide binaries to maximise HJ yield. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 12 figures, 2 tables, submitted to MNRAS

TESS and CHEOPS Discover Two Warm Sub-Neptunes Transiting the Bright K-dwarf HD 15906

Authors: Amy Tuson, Didier Queloz, Hugh P. Osborn, Thomas G. Wilson, Matthew J. Hooton, Mathias Beck, Monika Lendl, Göran Olofsson, Andrea Fortier, Andrea Bonfanti, Alexis Brandeker, Lars A. Buchhave, Andrew Collier Cameron, David R. Ciardi, Karen A. Collins, Davide Gandolfi, Zoltan Garai, Steven Giacalone, João Gomes da Silva, Steve B. Howell, Jayshil A. Patel, Carina M. Persson, Luisa M. Serrano, Sérgio G. Sousa, Solène Ulmer-Moll , et al. (97 additional authors not shown)

Abstract: We report the discovery of two warm sub-Neptunes transiting the bright (G = 9.5 mag) K-dwarf HD 15906 (TOI 461, TIC 4646810). This star was observed by the Transiting Exoplanet Survey Satellite (TESS) in sectors 4 and 31, revealing two small transiting planets. The inner planet, HD 15906 b, was detected with an unambiguous period but the outer planet, HD 15906 c, showed only two transits separated… ▽ More We report the discovery of two warm sub-Neptunes transiting the bright (G = 9.5 mag) K-dwarf HD 15906 (TOI 461, TIC 4646810). This star was observed by the Transiting Exoplanet Survey Satellite (TESS) in sectors 4 and 31, revealing two small transiting planets. The inner planet, HD 15906 b, was detected with an unambiguous period but the outer planet, HD 15906 c, showed only two transits separated by $\sim$ 734 days, leading to 36 possible values of its period. We performed follow-up observations with the CHaracterising ExOPlanet Satellite (CHEOPS) to confirm the true period of HD 15906 c and improve the radius precision of the two planets. From TESS, CHEOPS and additional ground-based photometry, we find that HD 15906 b has a radius of 2.24 $\pm$ 0.08 R$_\oplus$ and a period of 10.924709 $\pm$ 0.000032 days, whilst HD 15906 c has a radius of 2.93$^{+0.07}_{-0.06}$ R$_\oplus$ and a period of 21.583298$^{+0.000052}_{-0.000055}$ days. Assuming zero bond albedo and full day-night heat redistribution, the inner and outer planet have equilibrium temperatures of 668 $\pm$ 13 K and 532 $\pm$ 10 K, respectively. The HD 15906 system has become one of only six multiplanet systems with two warm ($\lesssim$ 700 K) sub-Neptune sized planets transiting a bright star (G $\leq$ 10 mag). It is an excellent target for detailed characterisation studies to constrain the composition of sub-Neptune planets and test theories of planet formation and evolution. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 30 pages, 20 figures, 11 tables (including appendix). Published in MNRAS

Journal ref: MNRAS, volume 523, issue 2, pp 3090-3118 (2023)

Refined parameters of the HD 22946 planetary system and the true orbital period of planet d

Authors: Z. Garai, H. P. Osborn, D. Gandolfi, A. Brandeker, S. G. Sousa, M. Lendl, A. Bekkelien, C. Broeg, A. Collier Cameron, J. A. Egger, M. J. Hooton, Y. Alibert, L. Delrez, L. Fossati, S. Salmon, T. G. Wilson, A. Bonfanti, A. Tuson, S. Ulmer-Moll, L. M. Serrano, L. Borsato, R. Alonso, G. Anglada, J. Asquier, D. Barrado y Navascues , et al. (63 additional authors not shown)

Abstract: Multi-planet systems are important sources of information regarding the evolution of planets. However, the long-period planets in these systems often escape detection. HD 22946 is a bright star around which 3 transiting planets were identified via TESS photometry, but the true orbital period of the outermost planet d was unknown until now. We aim to use CHEOPS to uncover the true orbital period of… ▽ More Multi-planet systems are important sources of information regarding the evolution of planets. However, the long-period planets in these systems often escape detection. HD 22946 is a bright star around which 3 transiting planets were identified via TESS photometry, but the true orbital period of the outermost planet d was unknown until now. We aim to use CHEOPS to uncover the true orbital period of HD 22946d and to refine the orbital and planetary properties of the system, especially the radii of the planets. We used the available TESS photometry of HD 22946 and observed several transits of the planets b, c, and d using CHEOPS. We identified 2 transits of planet d in the TESS photometry, calculated the most probable period aliases based on these data, and then scheduled CHEOPS observations. The photometric data were supplemented with ESPRESSO radial velocity data. Finally, a combined model was fitted to the entire dataset. We successfully determined the true orbital period of the planet d to be 47.42489 $\pm$ 0.00011 d, and derived precise radii of the planets in the system, namely 1.362 $\pm$ 0.040 R$_\oplus$, 2.328 $\pm$ 0.039 R$_\oplus$, and 2.607 $\pm$ 0.060 R$_\oplus$ for planets b, c, and d, respectively. Due to the low number of radial velocities, we were only able to determine 3$σ$ upper limits for these respective planet masses, which are 13.71 M$_\oplus$, 9.72 M$_\oplus$, and 26.57 M$_\oplus$. We estimated that another 48 ESPRESSO radial velocities are needed to measure the predicted masses of all planets in HD 22946. Planet c appears to be a promising target for future atmospheric characterisation. We can also conclude that planet d, as a warm sub-Neptune, is very interesting because there are only a few similar confirmed exoplanets to date. Such objects are worth investigating in the near future, for example in terms of their composition and internal structure. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Journal ref: A&A 674, A44 (2023)

Two Warm Neptunes transiting HIP 9618 revealed by TESS & Cheops

Authors: Hugh P. Osborn, Grzegorz Nowak, Guillaume Hébrard, Thomas Masseron, J. Lillo-Box, Enric Pallé, Anja Bekkelien, Hans-Gustav Florén, Pascal Guterman, Attila E. Simon, V. Adibekyan, Allyson Bieryla, Luca Borsato, Alexis Brandeker, David R. Ciardi, Andrew Collier Cameron, Karen A. Collins, Jo A. Egger, Davide Gandolfi, Matthew J. Hooton, David W. Latham, Monika Lendl, Elisabeth C. Matthews, Amy Tuson, Solène Ulmer-Moll , et al. (104 additional authors not shown)

Abstract: HIP 9618 (HD 12572, TOI-1471, TIC 306263608) is a bright ($G=9.0$ mag) solar analogue. TESS photometry revealed the star to have two candidate planets with radii of $3.9 \pm 0.044$ $R_\oplus$ (HIP 9618 b) and $3.343 \pm 0.039$ $R_\oplus$ (HIP 9618 c). While the 20.77291 day period of HIP 9618 b was measured unambiguously, HIP 9618 c showed only two transits separated by a 680-day gap in the time s… ▽ More HIP 9618 (HD 12572, TOI-1471, TIC 306263608) is a bright ($G=9.0$ mag) solar analogue. TESS photometry revealed the star to have two candidate planets with radii of $3.9 \pm 0.044$ $R_\oplus$ (HIP 9618 b) and $3.343 \pm 0.039$ $R_\oplus$ (HIP 9618 c). While the 20.77291 day period of HIP 9618 b was measured unambiguously, HIP 9618 c showed only two transits separated by a 680-day gap in the time series, leaving many possibilities for the period. To solve this issue, CHEOPS performed targeted photometry of period aliases to attempt to recover the true period of planet c, and successfully determined the true period to be 52.56349 d. High-resolution spectroscopy with HARPS-N, SOPHIE and CAFE revealed a mass of $10.0 \pm 3.1 M_\oplus$ for HIP 9618 b, which, according to our interior structure models, corresponds to a $6.8\pm1.4\%$ gas fraction. HIP 9618 c appears to have a lower mass than HIP 9618 b, with a 3-sigma upper limit of $< 18M_\oplus$. Follow-up and archival RV measurements also reveal a clear long-term trend which, when combined with imaging and astrometric information, reveal a low-mass companion ($0.08^{+0.12}_{-0.05} M_\odot$) orbiting at $26^{+19}_{-11}$ au. This detection makes HIP 9618 one of only five bright ($K<8$ mag) transiting multi-planet systems known to host a planet with $P>50$ d, opening the door for the atmospheric characterisation of warm ($T_{\rm eq}<750$ K) sub-Neptunes. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 19 pages, 16 figures, 9 tables. Accepted at MNRAS. CHEOPS, RV and ground-based photometric data is available on CDS at https://cdsarc.cds.unistra.fr/viz-bin/cat/J/MNRAS/523/3069

Journal ref: MNRAS, Vol. 523, 2023, issue 2, pp 3069-3089

TOI-5678 b: A 48-day transiting Neptune-mass planet characterized with CHEOPS and HARPS

Authors: S. Ulmer-Moll, H. P. Osborn, A. Tuson, J. A. Egger, M. Lendl, P. Maxted, A. Bekkelien, A. E. Simon, G. Olofsson, V. Adibekyan, Y. Alibert, A. Bonfanti, F. Bouchy, A. Brandeker, M. Fridlund, D. Gandolfi, C. Mordasini, C. M. Persson, S. Salmon, L. M. Serrano, S. G. Sousa, T. G. Wilson, M. Rieder, J. Hasiba, J. Asquier , et al. (70 additional authors not shown)

Abstract: A large sample of long-period giant planets has been discovered thanks to long-term radial velocity surveys, but only a few dozen of these planets have a precise radius measurement. Transiting gas giants are crucial targets for the study of atmospheric composition across a wide range of equilibrium temperatures and for shedding light on the formation and evolution of planetary systems. Indeed, com… ▽ More A large sample of long-period giant planets has been discovered thanks to long-term radial velocity surveys, but only a few dozen of these planets have a precise radius measurement. Transiting gas giants are crucial targets for the study of atmospheric composition across a wide range of equilibrium temperatures and for shedding light on the formation and evolution of planetary systems. Indeed, compared to hot Jupiters, the atmospheric properties and orbital parameters of cooler gas giants are unaltered by intense stellar irradiation and tidal effects. We identify long-period planets in the Transiting Exoplanet Survey Satellite (TESS) data as duo-transit events. To solve the orbital periods of TESS duo-transit candidates, we use the CHaracterising ExOPlanet Satellite (CHEOPS) to observe the highest-probability period aliases in order to discard or confirm a transit event at a given period. We also collect spectroscopic observations with CORALIE and HARPS in order to confirm the planetary nature and measure the mass of the candidates. We report the discovery of a warm transiting Neptune-mass planet orbiting TOI-5678. After four non-detections corresponding to possible periods, CHEOPS detected a transit event matching a unique period alias. Joint modeling reveals that TOI-5678 hosts a 47.73 day period planet. TOI-5678 b has a mass of 20 (+-4) Me and a radius of 4.91 (+-0.08 Re) . Using interior structure modeling, we find that TOI-5678 b is composed of a low-mass core surrounded by a large H/He layer with a mass of 3.2 (+1.7, -1.3) Me. TOI-5678 b is part of a growing sample of well-characterized transiting gas giants receiving moderate amounts of stellar insolation (11 Se). Precise density measurement gives us insight into their interior composition, and the objects orbiting bright stars are suitable targets to study the atmospheric composition of cooler gas giants. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 17 pages, 10 figures, accepted to A&A

Convergence of small scale Ly$α$ structure at high-$z$ under different reionization scenarios

Authors: Caitlin C. Doughty, Joseph F. Hennawi, Frederick B. Davies, Zarija Lukić, Jose Oñorbe

Abstract: The Ly$α$ forest (LAF) at $z>5$ probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-$z$, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting f… ▽ More The Ly$α$ forest (LAF) at $z>5$ probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-$z$, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting for a late end to reionization, small structures down to the sub-kpc level may survive to later times than conventionally thought due to the reduced time for pressure smoothing to impact the gas, further tightening simulation resolution requirements. We perform a suite of simulations using the Eulerian cosmological hydrodynamics code Nyx, spanning domain sizes of 1.25-10 $h^{-1}$ Mpc and 5-80 $h^{-1}$ kpc cells, and explore the interaction of these variables with the timing of reionization on the properties of the matter distribution and the simulated LAF at $z=5.5$. In observable Ly$α$ power, convergence within 10% is achieved for $k< 0.1$ s/km, but larger $k$ shows deviation of up to 20 percent. While a later reionization retains more small structure in the density field, because of the greater thermal broadening there is little difference in the convergence of LAF power between early ($z=9$) and later ($z=6$) reionizations. We conclude that at $z\sim5.5$, resolutions of 10 kpc are necessary for convergence of LAF power at $k<0.1$ s/km, while higher-$k$ modes require higher resolution, and that the timing of reionization does not significantly impact convergence given realistic photoheating. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: Submitted to MNRAS. Comments welcome