Showing 1–50 of 232 results for author: Bolton, S

The Lyman-limit photon mean free path at the end of late reionisation in the Sherwood-Relics simulations

Authors: Jennifer Feron, Luke Conaboy, James S. Bolton, Emma Chapman, Martin G. Haehnelt, Laura C. Keating, Girish Kulkarni, Ewald Puchwein

Abstract: Recent evidence supporting reionisation ending at redshift z<6 includes the rapid redshift evolution of the mean free path for Lyman-limit photons through the intergalactic medium at 5<z<6. Here we investigate the mean free path predicted by the Sherwood-Relics suite of hybrid radiation hydrodynamical simulations. Simulations with comoving volumes of 40^3 h^-3 cMpc^3 (160^3 h^-3 cMpc^3), calibrate… ▽ More Recent evidence supporting reionisation ending at redshift z<6 includes the rapid redshift evolution of the mean free path for Lyman-limit photons through the intergalactic medium at 5<z<6. Here we investigate the mean free path predicted by the Sherwood-Relics suite of hybrid radiation hydrodynamical simulations. Simulations with comoving volumes of 40^3 h^-3 cMpc^3 (160^3 h^-3 cMpc^3), calibrated to match the observed Lyman-alpha forest transmission with a late end to reionisation at z<6, are consistent with recent mean free path measurements at z<5.9, and are 1.2σ(1.8σ) above the highest redshift mean free path measurement at z=5.93. The majority of the Lyman-limit opacity at the end of reionisation is attributable to highly ionised Lyman-alpha forest absorbers with neutral hydrogen column densities N_HI=10^16--10^17 cm^-2. Resolving these systems is critical for capturing the redshift evolution of the mean free path during the final stages of reionisation. After reionisation completes, overdense gas will reduce the mean free path by up to 20 per cent around haloes with masses M_h ~ 10^9--10^11 h^-1 M_sol, but during reionisation ionised bubbles will instead boost the mean free path around haloes by up to an order of magnitude when the IGM is as much as 90 per cent neutral by volume. This effect will play an important role in the visibility of Lyman-alpha emitting galaxies at z>10 discovered with JWST. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 17 pages, 14 figures, submitted to MNRAS

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 wing absorption associated with a giant Ly$α$ trough at $z < 6$: direct evidence for late-ending reionization

Authors: George D. Becker, James S. Bolton, Yongda Zhu, Seyedazim Hashemi

Abstract: Multiple observations now suggest that the hydrogen reionization may have ended well below redshift six. While there has previously been no conclusive proof of extended neutral islands in the $z < 6$ intergalactic medium, it is possible that such islands give rise to the giant Ly$α$ absorption troughs seen in the spectra of high-redshift quasars. Here we present evidence that the deepest and longe… ▽ More Multiple observations now suggest that the hydrogen reionization may have ended well below redshift six. While there has previously been no conclusive proof of extended neutral islands in the $z < 6$ intergalactic medium, it is possible that such islands give rise to the giant Ly$α$ absorption troughs seen in the spectra of high-redshift quasars. Here we present evidence that the deepest and longest-known Ly$α$ trough at $z < 6$, towards ULAS J0148+0600 (J0148), is associated with damping wing absorption. The evidence comes from a window of strong Ly$α$ transmission at the edge of the J0148 proximity zone. We show that the relatively smooth profile of this transmission window is highly unlikely to arise from resonant absorption alone, but is consistent with the presence of a damping wing. We further argue that the damping wing is unlikely to arise from a compact source due to the lack of associated metal lines, and is more likely to arise from an extended neutral island associated with the giant Ly$α$ trough. We investigate the physical conditions that may give rise to the strong transmission window, and speculate that it may signal an usually deep void, nearby ionizing sources, and/or the recent passage of an ionization front. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 17 pages, 18 figures; resubmitted to MNRAS after addressing referee comments

Properties of electrons accelerated by the Ganymede-magnetosphere interaction: survey of Juno high-latitude observations

Authors: J. Rabia, V. Hue, N. Andre, Q. Nenon, J. R. Szalay, F. Allegrini, A. H. Sulaiman, C. K. Louis, T. K. Greathouse, Y. Sarkango, D. Santos-Costa, M. Blanc, E. Penou, P. Louarn, R. W. Ebert, G. R. Gladstone, A. Mura, J. E. P. Connerney, S. J. Bolton

Abstract: The encounter between the Jovian co-rotating plasma and Ganymede gives rise to electromagnetic waves that propagate along the magnetic field lines and accelerate particles by resonant or non-resonant wave-particle interaction. They ultimately precipitate into Jupiter's atmosphere and trigger auroral emissions. In this study, we use Juno/JADE, Juno/UVS data, and magnetic field line tracing to chara… ▽ More The encounter between the Jovian co-rotating plasma and Ganymede gives rise to electromagnetic waves that propagate along the magnetic field lines and accelerate particles by resonant or non-resonant wave-particle interaction. They ultimately precipitate into Jupiter's atmosphere and trigger auroral emissions. In this study, we use Juno/JADE, Juno/UVS data, and magnetic field line tracing to characterize the properties of electrons accelerated by the Ganymede-magnetosphere interaction in the far-field region. We show that the precipitating energy flux exhibits an exponential decay as a function of downtail distance from the moon, with an e-folding value of 29°, consistent with previous UV observations from the Hubble Space Telescope (HST). We characterize the electron energy distributions and show that two distributions exist. Electrons creating the Main Alfvén Wing (MAW) spot and the auroral tail always have broadband distribution and a mean characteristic energy of 2.2 keV while in the region connected to the Transhemispheric Electron Beam (TEB) spot the electrons are distributed non-monotonically, with a higher characteristic energy above 10 keV. Based on the observation of bidirectional electron beams, we suggest that Juno was located within the acceleration region during the 11 observations reported. We thus estimate that the acceleration region is extended, at least, between an altitude of 0.5 and 1.3 Jupiter radius above the 1-bar surface. Finally, we estimate the size of the interaction region in the Ganymede orbital plane using far-field measurements. These observations provide important insights for the study of particle acceleration processes involved in moon-magnetosphere interactions. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: Accepted for publication in JGR: Space Physics

Late-end reionization with ATON-HE: towards constraints from Lyman-$α$ emitters observed with JWST

Authors: Shikhar Asthana, Martin G. Haehnelt, Girish Kulkarni, Dominique Aubert, James S. Bolton, Laura C. Keating

Abstract: We present a new suite of late-end reionization simulations performed with ATON-HE, a revised version of the GPU-based radiative transfer code ATON that includes helium. The simulations are able to reproduce the Ly$α$ flux distribution of the E-XQR-30 sample of QSO absorption spectra at $5 \lesssim z \lesssim 6.2$, and show that a large variety of reionization models are consistent with these data… ▽ More We present a new suite of late-end reionization simulations performed with ATON-HE, a revised version of the GPU-based radiative transfer code ATON that includes helium. The simulations are able to reproduce the Ly$α$ flux distribution of the E-XQR-30 sample of QSO absorption spectra at $5 \lesssim z \lesssim 6.2$, and show that a large variety of reionization models are consistent with these data. We explore a range of variations in source models and in the early-stage evolution of reionization. Our fiducial reionization history has a midpoint of reionization at $z = 6.5$, but we also explore an `Early' reionization history with a midpoint at $z = 7.5$ and an `Extremely Early' reionization history with a midpoint at $z = 9.5$. Haloes massive enough to host observed Ly$α$ emitters are highly biased. The fraction of such haloes embedded in ionized bubbles that are large enough to allow high Ly$α$ transmission becomes close to unity much before the volume filling factor of ionized regions. For our fiducial reionization history this happens at $z = 8$, probably too late to be consistent with the detection by JWST of abundant Ly$α$ emission out to $z = 11$. A reionization history in our `Early' model or perhaps even our `Extremely Early' model may be required, suggesting a Thomson scattering optical depth in tension with that reported by Planck, but consistent with recent suggestions of a significantly higher value. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 25 pages, 15 figures

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

Super-adiabatic Temperature Gradient at Jupiter's Equatorial Zone and Implications for the Water Abundance

Authors: Cheng Li, Michael Allison, Sushil Atreya, Shawn Brueshaber, Leigh N. Fletcher, Tristan Guillot, Liming Li, Jonathan Lunine, Yamila Miguel, Glenn Orton, Paul Steffes, J. Hunter Waite, Michael H. Wong, Steven Levin, Scott Bolton

Abstract: The temperature structure of a giant planet was traditionally thought to be an adiabat assuming convective mixing homogenizes entropy. The only in-situ measurement made by the Galileo Probe detected a near-adiabatic temperature structure within one of Jupiter's 5$μ$m hot spots with small but definite local departures from adiabaticity. We analyze Juno's microwave observations near Jupiter's equato… ▽ More The temperature structure of a giant planet was traditionally thought to be an adiabat assuming convective mixing homogenizes entropy. The only in-situ measurement made by the Galileo Probe detected a near-adiabatic temperature structure within one of Jupiter's 5$μ$m hot spots with small but definite local departures from adiabaticity. We analyze Juno's microwave observations near Jupiter's equator (0 ~ 5$^o$N) and find that the equatorial temperature structure is best characterized by a stable super-adiabatic temperature profile rather than an adiabatic one. Water is the only substance with sufficient abundance to alter the atmosphere's mean molecular weight and prevent dynamic instability if a super-adiabatic temperature gradient exists. Thus, from the super-adiabaticity, our results indicate a water concentration (or the oxygen to hydrogen ratio) of about 4.9 times solar with a possible range of 1.5 ~ 8.3 times solar in Jupiter's equatorial region. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Comments: 46 pages, 11 figures, accepted by Icraus

SPARCL: SPectra Analysis and Retrievable Catalog Lab

Authors: Stéphanie Juneau, Alice Jacques, Steve Pothier, Adam S. Bolton, Benjamin A. Weaver, Ragadeepika Pucha, Sean McManus, Robert Nikutta, Knut Olsen

Abstract: SPectra Analysis and Retrievable Catalog Lab (SPARCL) at NOIRLab's Astro Data Lab was created to efficiently serve large optical and infrared spectroscopic datasets. It consists of services, tools, example workflows and currently contains spectra for over 7.5 million stars, galaxies and quasars from the Sloan Digital Sky Survey (SDSS) and the Dark Energy Spectroscopic Instrument (DESI) survey. We… ▽ More SPectra Analysis and Retrievable Catalog Lab (SPARCL) at NOIRLab's Astro Data Lab was created to efficiently serve large optical and infrared spectroscopic datasets. It consists of services, tools, example workflows and currently contains spectra for over 7.5 million stars, galaxies and quasars from the Sloan Digital Sky Survey (SDSS) and the Dark Energy Spectroscopic Instrument (DESI) survey. We aim to eventually support the broad range of spectroscopic datasets that will be hosted at NOIRLab and beyond. Major elements of SPARCL include capabilities to discover and query for spectra based on parameters of interest, a fast web service that delivers desired spectra either individually or in bulk as well as documentation and example Jupyter Notebooks to empower users in their research. More information is available on the SPARCL website (https://astrosparcl.datalab.noirlab.edu). △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: 4 pages, 1 figure, Conference Proceedings for ADASS 2023 (Astronomical Data Analysis Software & Systems XXXIII)

Galaxy Formation and Symbiotic Evolution with the Inter-Galactic Medium in the Age of ELT-ANDES

Authors: Valentina D'Odorico, James S. Bolton, Lise Christensen, Annalisa De Cia, Erik Zackrisson, Aron Kordt, Luca Izzo, Jiangtao Li, Roberto Maiolino, Alessandro Marconi, Philipp Richter, Andrea Saccardi, Stefania Salvadori, Irene Vanni, Chiara Feruglio, Michele Fumagalli, Johan P. U. Fynbo, Pasquier Noterdaeme, Polychronis Papaderos, Celine Peroux, Aprajita Verma, Paolo Di Marcantonio, Livia Origlia, Alessio Zanutta

Abstract: High-resolution absorption spectroscopy toward bright background sources has had a paramount role in understanding early galaxy formation, the evolution of the intergalactic medium and the reionisation of the Universe. However, these studies are now approaching the boundaries of what can be achieved at ground-based 8-10m class telescopes. The identification of primeval systems at the highest redsh… ▽ More High-resolution absorption spectroscopy toward bright background sources has had a paramount role in understanding early galaxy formation, the evolution of the intergalactic medium and the reionisation of the Universe. However, these studies are now approaching the boundaries of what can be achieved at ground-based 8-10m class telescopes. The identification of primeval systems at the highest redshifts, within the reionisation epoch and even into the dark ages, and of the products of the first generation of stars and the chemical enrichment of the early Universe, requires observing very faint targets with a signal-to-noise ratio high enough to detect very faint spectral signatures. In this paper, we describe the giant leap forward that will be enabled by ANDES, the high-resolution spectrograph for the ELT, in these key science fields, together with a brief, non-exhaustive overview of other extragalactic research topics that will be pursued by this instrument, and its synergistic use with other facilities that will become available in the early 2030s. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 40 pages, 7 figures; submitted to Experimental Astronomy on behalf of the ANDES Science Team

Strong resemblance between surface and deep zonal winds inside Jupiter revealed by high-degree gravity moments

Authors: Hao Cao, Jeremy Bloxham, Ryan S. Park, Burkhard Militzer, Rakesh K. Yadav, Laura Kulowski, David J. Stevenson, Scott J. Bolton

Abstract: Jupiter's atmosphere-interior is a coupled fluid dynamical system strongly influenced by the rapid background rotation. While the visible atmosphere features east-west zonal winds on the order of 100 m/s (Tollefson et al. 2017), zonal flows in the dynamo region are significantly slower, on the order of 1 cm/s or less, according to the latest magnetic secular variation analysis (Bloxham et al. 2022… ▽ More Jupiter's atmosphere-interior is a coupled fluid dynamical system strongly influenced by the rapid background rotation. While the visible atmosphere features east-west zonal winds on the order of 100 m/s (Tollefson et al. 2017), zonal flows in the dynamo region are significantly slower, on the order of 1 cm/s or less, according to the latest magnetic secular variation analysis (Bloxham et al. 2022). The vertical profile of the zonal flows and the underlying mechanism remain elusive. The latest Juno radio tracking measurements afforded the derivation of Jupiter's gravity field to spherical harmonic degree 40. Here, we use the latest gravity solution to reconstruct Jupiter's deep zonal winds without a priori assumptions about their latitudinal profile. The pattern of our reconstructed deep zonal winds strongly resembles that of the surface wind within $\pm$ 35 degrees latitude from the equator, in particular the northern off-equatorial jet (NOEJ) and the southern off-equatorial jet (SOEJ) (Kulowski et al. 2021). The reconstruction features larger uncertainties in the southern hemisphere due to the north south asymmetric nature of Juno's trajectory. Amplitude of the reconstructed deep NOEJ matches that of the surface wind when the wind is truncated at a depth around 2500 km, and becomes twice that of the surface wind if the truncation depth is reduced to about 1500 km. Our analysis supports the physical picture in which prominent part of the surface zonal winds extends into Jupiter's interior significantly deeper than the water cloud layer. △ Less

Submitted 19 November, 2023; originally announced November 2023.

Comments: Accepted for publication in ApJ, 11 pages, 6 figures

Journal ref: ApJ 959 78 (2023)

Particle initialization effects on Lyman-$α$ forest statistics in cosmological SPH simulations

Authors: Nabendu Kumar Khan, Girish Kulkarni, James S. Bolton, Martin G. Haehnelt, Vid Iršič, Ewald Puchwein, Shikhar Asthana

Abstract: Confronting measurements of the Lyman-$α$ forest with cosmological hydrodynamical simulations has produced stringent constraints on models of particle dark matter and the thermal and ionization state of the intergalactic medium. We investigate the robustness of such models of the Lyman-$α$ forest, focussing on the effect of particle initial conditions on the Lyman-$α$ forest statistics in cosmolog… ▽ More Confronting measurements of the Lyman-$α$ forest with cosmological hydrodynamical simulations has produced stringent constraints on models of particle dark matter and the thermal and ionization state of the intergalactic medium. We investigate the robustness of such models of the Lyman-$α$ forest, focussing on the effect of particle initial conditions on the Lyman-$α$ forest statistics in cosmological SPH simulations. We study multiple particle initialization algorithms in simulations that are designed to be identical in other respects. In agreement with the literature, we find that the correct linear theory evolution is obtained when a glass-like configuration is used for initial unperturbed gas particle positions alongside a regular grid configuration for dark matter particles and the use of non-identical initial density perturbations for gas and dark matter. However, we report that this introduces a large scale-dependent distortion in the one-dimensional Lyman-$α$ transmission power spectrum at small scales ($k > 0.05$ s/km). The effect is close to $50\%$ at $k\sim 0.1$ s/km, and persists at higher resolution. This can severely bias inferences in parameters such as the dark matter particle mass. By considering multiple initial conditions codes and their variations, we also study the impact of a variety of other assumptions and algorithmic choices, such as adaptive softening, background radiation density, particle staggering, and perturbation theory accuracy, on the matter power spectrum, the Lyman-$α$ flux power spectrum, and the Lyman-$α$ flux PDF. This work reveals possible pathways towards more accurate theoretical models of the Lyman-$α$ forest to match the quality of upcoming measurements. △ Less

Submitted 29 February, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: 18 pages, 16 figures, accepted for publication in MNRAS, comments welcome

Unveiling Dark Matter free-streaming at the smallest scales with high redshift Lyman-alpha forest

Authors: Vid Iršič, Matteo Viel, Martin G. Haehnelt, James S. Bolton, Margherita Molaro, Ewald Puchwein, Elisa Boera, George D. Becker, Prakash Gaikwad, Laura C. Keating, Girish Kulkarni

Abstract: This study introduces novel constraints on the free-streaming of thermal relic warm dark matter (WDM) from Lyman-$α$ forest flux power spectra. Our analysis utilises a high-resolution, high-redshift sample of quasar spectra observed using the HIRES and UVES spectrographs ($z=4.2-5.0$). We employ a Bayesian inference framework and a simulation-based likelihood that encompasses various parameters in… ▽ More This study introduces novel constraints on the free-streaming of thermal relic warm dark matter (WDM) from Lyman-$α$ forest flux power spectra. Our analysis utilises a high-resolution, high-redshift sample of quasar spectra observed using the HIRES and UVES spectrographs ($z=4.2-5.0$). We employ a Bayesian inference framework and a simulation-based likelihood that encompasses various parameters including the free-streaming of dark matter, cosmological parameters, the thermal history of the intergalactic medium, and inhomogeneous reionization, to establish lower limits on the mass of a thermal relic WDM particle of $5.7\;\mathrm{keV}$ (at 95\% C.L.). This result surpasses previous limits from the Lyman-$α$ forest through reduction of the measured uncertainties due to a larger statistical sample and by measuring clustering to smaller scales ($k_{\rm max}=0.2\;\mathrm{km^{-1}\,s}$). The approximately two-fold improvement due to the expanded statistical sample suggests that the effectiveness of Lyman-$α$ forest constraints on WDM models at high redshifts are limited by the availability of high-quality quasar spectra. Restricting the analysis to comparable scales and thermal history priors as in prior studies ($k_{\rm max}<0.1\;\mathrm{km^{-1}\,s}$) lowers the bound on the WDM mass to $4.1\;\mathrm{keV}$. As the precision of the measurements increases, it becomes crucial to examine the instrumental and modelling systematics. On the modelling front, we argue that the impact of the thermal history uncertainty on the WDM particle mass constraint has diminished due to improved independent observations. At the smallest scales, the primary source of modeling systematic arises from the structure in the peculiar velocity of the intergalactic medium and inhomogeneous reionization. △ Less

Submitted 15 February, 2024; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 22 pages, 14 figures, 2 tables; published in PRD on Feb 8 2024

Journal ref: Phys. Rev. D, 109, 043511 (2024)

Investigating the characteristic shape and scatter of intergalactic damping wings during reionization

Authors: Laura C. Keating, Ewald Puchwein, James S. Bolton, Martin G. Haehnelt, Girish Kulkarni

Abstract: Ly$α$ damping wings in the spectra of bright objects at high redshift are a useful probe of the ionization state of the intergalactic medium during the reionization epoch. It has recently been noted that, despite the inhomogeneous nature of reionization, these damping wings have a characteristic shape which is a strong function of the volume-weighted average neutral hydrogen fraction of the interg… ▽ More Ly$α$ damping wings in the spectra of bright objects at high redshift are a useful probe of the ionization state of the intergalactic medium during the reionization epoch. It has recently been noted that, despite the inhomogeneous nature of reionization, these damping wings have a characteristic shape which is a strong function of the volume-weighted average neutral hydrogen fraction of the intergalactic medium. We present here a closer examination of this finding using a simulation of patchy reionization from the Sherwood-Relics simulation suite. We show that the characteristic shape and scatter of the damping wings are determined by the average neutral hydrogen density along the line of sight, weighted by its contribution to the optical depth producing the damping wing. We find that there is a redshift dependence in the characteristic shape due to the expansion of the Universe. Finally, we show that it is possible to differentiate between the shapes of damping wings in galaxies and young (or faint) quasars at different points in the reionization history at large velocity offsets from the point where the transmission first reaches zero. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 5 pages, 4 figures. Submitted to MNRAS Letters

JWST observations of galaxy damping wings during reionization interpreted with cosmological simulations

Authors: Laura C. Keating, James S. Bolton, Fergus Cullen, Martin G. Haehnelt, Ewald Puchwein, Girish Kulkarni

Abstract: Spectra of the highest redshift galaxies taken with JWST are now allowing us to see into the heart of the reionization epoch. Many of these observed galaxies exhibit strong damping wing absorption redward of their Lyman-$α$ emission. These observations have been used to measure the redshift evolution of the neutral fraction of the intergalactic medium and sizes of ionized bubbles. However, these e… ▽ More Spectra of the highest redshift galaxies taken with JWST are now allowing us to see into the heart of the reionization epoch. Many of these observed galaxies exhibit strong damping wing absorption redward of their Lyman-$α$ emission. These observations have been used to measure the redshift evolution of the neutral fraction of the intergalactic medium and sizes of ionized bubbles. However, these estimates have been made using a simple analytic model for the intergalactic damping wing. We explore the recent observations with models of inhomogeneous reionization from the Sherwood-Relics simulation suite. We carry out a comparison between the damping wings calculated from the simulations and from the analytic model. We find that although the agreement is good on the red side of the Lyman-$α$ emission, there is a discrepancy on the blue side due to residual neutral hydrogen present in the simulations, which saturates the intergalactic absorption. For this reason, we find that it is difficult to reproduce the claimed observations of large bubble sizes at z ~ 7, which are driven by a detection of transmitted flux blueward of the Lyman-$α$ emission. We suggest instead that the observations can be explained by a model with smaller ionized bubbles and larger intrinsic Lyman-$α$ emission from the host galaxy. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: 13 pages, 7 figures. Submitted to MNRAS

Spacecraft VLBI tracking to enhance stellar occultations astrometry of planetary satellites

Authors: M. Fayolle, V. Lainey, D. Dirkx, L. I. Gurvits, G. Cimo, S. J. Bolton

Abstract: Stellar occultations currently provide the most accurate ground-based measurements of the positions of natural satellites (down to a few kilometres for the Galilean moons). However, when using these observations in the calculation of satellite ephemerides, the uncertainty in the planetary ephemerides dominates the error budget of the occultation. We quantify the local refinement in the central pla… ▽ More Stellar occultations currently provide the most accurate ground-based measurements of the positions of natural satellites (down to a few kilometres for the Galilean moons). However, when using these observations in the calculation of satellite ephemerides, the uncertainty in the planetary ephemerides dominates the error budget of the occultation. We quantify the local refinement in the central planet's position achievable by performing Very Long Baseline Interferometry (VLBI) tracking of an in-system spacecraft temporally close to an occultation. We demonstrate the potential of using VLBI to enhance the science return of stellar occultations for satellite ephemerides. We identified the most promising observation and tracking opportunities offered by the Juno spacecraft around Jupiter as perfect test cases, for which we ran simulations of our VLBI experiment. VLBI tracking at Juno's perijove close to a stellar occultation locally (in time) reduces the uncertainty in Jupiter's angular position in the sky to 250-400 m. This represents up to an order of magnitude improvement with respect to current solutions and is lower than the stellar occultation error, thus allowing the moon ephemeris solution to fully benefit from the observation. Our simulations showed that the proposed tracking and observation experiment can efficiently use synergies between ground- and space-based observations to enhance the science return on both ends. The reduced error budget for stellar occultations indeed helps to improve the moons' ephemerides, which in turn benefit planetary missions and their science products, such as the recently launched JUICE and upcoming Europa Clipper missions. △ Less

Submitted 29 July, 2023; originally announced July 2023.

Highly depleted alkali metals in Jupiter's deep atmosphere

Authors: Ananyo Bhattacharya, Cheng Li, Sushil K. Atreya, Paul G. Steffes, Steven M. Levin, Scott J. Bolton, Tristan Guillot, Pranika Gupta, Andrew P. Ingersoll, Jonathan I. Lunine, Glenn S. Orton, Fabiano A. Oyafuso, J. Hunter Waite, Amadeo Belloti, Michael H. Wong

Abstract: Water and ammonia vapors are known to be the major sources of spectral absorption at pressure levels observed by the microwave radiometer (MWR) on Juno. However, the brightness temperatures and limb darkening observed by the MWR at its longest wavelength channel of 50 cm (600 MHz) in the first 9 perijove passes indicate the existence of an additional source of opacity in the deep atmosphere of Jup… ▽ More Water and ammonia vapors are known to be the major sources of spectral absorption at pressure levels observed by the microwave radiometer (MWR) on Juno. However, the brightness temperatures and limb darkening observed by the MWR at its longest wavelength channel of 50 cm (600 MHz) in the first 9 perijove passes indicate the existence of an additional source of opacity in the deep atmosphere of Jupiter (pressures beyond 100 bar). The absorption properties of ammonia and water vapor, and their relative abundances in Jupiter's atmosphere do not provide sufficient opacity in deep atmosphere to explain the 600 MHz channel observation. Here we show that free electrons due to the ionization of alkali metals, i.e. sodium, and potassium, with sub-solar metallicity [M/H] (log based 10 relative concentration to solar) in the range of [M/H] = -2 to [M/H] = -5 can provide the missing source of opacity in the deep atmosphere. If the alkali metals are not the source of additional opacity in the MWR data, then their metallicity at 1000 bars can only be even lower. The upper bound of -2 on the metallicity of the alkali metals contrasts with the other heavy elements -- C, N, S, Ar, Kr, and Xe -- which are all enriched relative to their solar abundances having a metallicity of approximately +0.5. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: This manuscript has been accepted for publication in The Astrophysical Journal Letters. The final version of the paper will be available in the published journal. This arXiv version is provided for informational purposes

The Early Data Release of the Dark Energy Spectroscopic Instrument

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering, D. M. Alexander, R. Alfarsy, C. Allende Prieto, M. Alvarez, O. Alves, A. Anand, F. Andrade-Oliveira, E. Armengaud, J. Asorey, S. Avila, A. Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, J. Bautista, J. Behera, S. F. Beltran , et al. (240 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) completed its five-month Survey Validation in May 2021. Spectra of stellar and extragalactic targets from Survey Validation constitute the first major data sample from the DESI survey. This paper describes the public release of those spectra, the catalogs of derived properties, and the intermediate data products. In total, the public release includes… ▽ More The Dark Energy Spectroscopic Instrument (DESI) completed its five-month Survey Validation in May 2021. Spectra of stellar and extragalactic targets from Survey Validation constitute the first major data sample from the DESI survey. This paper describes the public release of those spectra, the catalogs of derived properties, and the intermediate data products. In total, the public release includes good-quality spectral information from 466,447 objects targeted as part of the Milky Way Survey, 428,758 as part of the Bright Galaxy Survey, 227,318 as part of the Luminous Red Galaxy sample, 437,664 as part of the Emission Line Galaxy sample, and 76,079 as part of the Quasar sample. In addition, the release includes spectral information from 137,148 objects that expand the scope beyond the primary samples as part of a series of secondary programs. Here, we describe the spectral data, data quality, data products, Large-Scale Structure science catalogs, access to the data, and references that provide relevant background to using these spectra. △ Less

Submitted 15 June, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 43 pages, 7 figures, 17 tables, submitted to AJ, DESI EDR references added

Validation of the Scientific Program for the Dark Energy Spectroscopic Instrument

Authors: DESI Collaboration, A. G. Adame, J. Aguilar, S. Ahlen, S. Alam, G. Aldering, D. M. Alexander, R. Alfarsy, C. Allende Prieto, M. Alvarez, O. Alves, A. Anand, F. Andrade-Oliveira, E. Armengaud, J. Asorey, S. Avila, A. Aviles, S. Bailey, A. Balaguera-Antolínez, O. Ballester, C. Baltay, A. Bault, J. Bautista, J. Behera, S. F. Beltran , et al. (239 additional authors not shown)

Abstract: The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of… ▽ More The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar (MWS), bright galaxy (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the five-year program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a `One-Percent survey' conducted at the conclusion of Survey Validation covering 140 deg$^2$ using the final target selection algorithms with exposures of a depth typical of the main survey. The Survey Validation indicates that DESI will be able to complete the full 14,000 deg$^2$ program with spectroscopically-confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval $z<1.1$, 0.39% over the redshift interval $1.1<z<1.9$, and 0.46% over the redshift interval $1.9<z<3.5$. △ Less

Submitted 12 January, 2024; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: 42 pages, 18 figures, accepted by AJ

The Io, Europa and Ganymede auroral footprints at Jupiter in the ultraviolet: positions and equatorial lead angles

Authors: Vincent Hue, Randy Gladstone, Corentin K. Louis, Thomas K. Greathouse, Bertrand Bonfond, Jamey R. Szalay, Alessandro Moirano, Rohini S. Giles, Joshua A. Kammer, Masafumi Imai, Alessandro Mura, Maarten H. Versteeg, George Clark, Jean-Claude Gérard, Denis C. Grodent, Jonas Rabia, Ali H. Sulaiman, Scott J. Bolton, John E. P. Connerney

Abstract: Jupiter's satellite auroral footprints are a consequence of the interaction between the Jovian magnetic field with co-rotating iogenic plasma and the Galilean moons. The disturbances created near the moons propagate as Alfvén waves along the magnetic field lines. The position of the moons is therefore "Alfvénically" connected to their respective auroral footprint. The angular separation from the i… ▽ More Jupiter's satellite auroral footprints are a consequence of the interaction between the Jovian magnetic field with co-rotating iogenic plasma and the Galilean moons. The disturbances created near the moons propagate as Alfvén waves along the magnetic field lines. The position of the moons is therefore "Alfvénically" connected to their respective auroral footprint. The angular separation from the instantaneous magnetic footprint can be estimated by the so-called lead angle. That lead angle varies periodically as a function of orbital longitude, since the time for the Alfvén waves to reach the Jovian ionosphere varies accordingly. Using spectral images of the Main Alfvén Wing auroral spots collected by Juno-UVS during the first forty-three orbits, this work provides the first empirical model of the Io, Europa and Ganymede equatorial lead angles for the northern and southern hemispheres. Alfvén travel times between the three innermost Galilean moons to Jupiter's northern and southern hemispheres are estimated from the lead angle measurements. We also demonstrate the accuracy of the mapping from the Juno magnetic field reference model (JRM33) at the completion of the prime mission for M-shells extending to at least 15RJ . Finally, we shows how the added knowledge of the lead angle can improve the interpretation of the moon-induced decametric emissions. △ Less

Submitted 28 April, 2023; originally announced April 2023.

Comments: 20 pages, 8 figures, Accepted for publication in Journal of Geophysical Research: Space Physics on 20 April 2023

Measuring the photo-ionization rate, neutral fraction and mean free path of HI ionizing photons at $4.9 \leq z \leq 6.0$ from a large sample of XShooter and ESI spectra

Authors: Prakash Gaikwad, Martin G. Haehnelt, Frederick B. Davies, Sarah E. I. Bosman, Margherita Molaro, Girish Kulkarni, Valentina D'Odorico, George D. Becker, Rebecca L. Davies, Fahad Nasir, James S. Bolton, Laura C. Keating, Vid Iršič, Ewald Puchwein, Yongda Zhu, Shikhar Asthana, Jinyi Yang, Samuel Lai, Anna-Christina Eilers

Abstract: We measure the mean free path ($λ_{\rm mfp,HI}$), photo-ionization rate ($\langle Γ_{\rm HI} \rangle$) and neutral fraction ($\langle f_{\rm HI} \rangle$) of hydrogen in 12 redshift bins at $4.85<z<6.05$ from a large sample of moderate resolution XShooter and ESI QSO absorption spectra. The fluctuations in ionizing radiation field are modeled by post-processing simulations from the Sherwood suite… ▽ More We measure the mean free path ($λ_{\rm mfp,HI}$), photo-ionization rate ($\langle Γ_{\rm HI} \rangle$) and neutral fraction ($\langle f_{\rm HI} \rangle$) of hydrogen in 12 redshift bins at $4.85<z<6.05$ from a large sample of moderate resolution XShooter and ESI QSO absorption spectra. The fluctuations in ionizing radiation field are modeled by post-processing simulations from the Sherwood suite using our new code "EXtended reionization based on the Code for Ionization and Temperature Evolution" (EX-CITE). EX-CITE uses efficient Octree summation for computing intergalactic medium attenuation and can generate large number of high resolution $Γ_{\rm HI}$ fluctuation models. Our simulation with EX-CITE shows remarkable agreement with simulations performed with the radiative transfer code Aton and can recover the simulated parameters within $1σ$ uncertainty. We measure the three parameters by forward-modeling the Ly$α$ forest and comparing the effective optical depth ($τ_{\rm eff, HI}$) distribution in simulations and observations. The final uncertainties in our measured parameters account for the uncertainties due to thermal parameters, modeling parameters, observational systematics and cosmic variance. Our best fit parameters show significant evolution with redshift such that $λ_{\rm mfp,HI}$ and $\langle f_{\rm HI} \rangle$ decreases and increases by a factor $\sim 6$ and $\sim 10^{4}$, respectively from $z \sim 5$ to $z \sim 6$. By comparing our $λ_{\rm mfp,HI}$, $\langle Γ_{\rm HI} \rangle$ and $\langle f_{\rm HI} \rangle$ evolution with that in state-of-the-art Aton radiative transfer simulations and the Thesan and CoDa-III simulations, we find that our best fit parameter evolution is consistent with a model in which reionization completes by $z \sim 5.2$. △ Less

Submitted 23 August, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: 30 pages (+14 pages appendices), 14 figures (+18 figures appendices); Accepted for publication in MNRAS; Main results are summarized in Fig. 10, Fig. 11 and Table 3

Possible evidence for a large-scale enhancement in the Lyman-$α$ forest power spectrum at redshift $\mathbf{\textit{z}\geq 4}$

Authors: Margherita Molaro, Vid Iršič, James S. Bolton, Maggie Lieu, Laura C. Keating, Ewald Puchwein, Martin G. Haehnelt, Matteo Viel

Abstract: Inhomogeneous reionization enhances the 1D Lyman-$α$ forest power spectrum on large scales at redshifts $z\geq4$. This is due to coherent fluctuations in the ionized hydrogen fraction that arise from large-scale variations in the post-reionization gas temperature, which fade as the gas cools. It is therefore possible to use these relic fluctuations to constrain inhomogeneous reionization with the… ▽ More Inhomogeneous reionization enhances the 1D Lyman-$α$ forest power spectrum on large scales at redshifts $z\geq4$. This is due to coherent fluctuations in the ionized hydrogen fraction that arise from large-scale variations in the post-reionization gas temperature, which fade as the gas cools. It is therefore possible to use these relic fluctuations to constrain inhomogeneous reionization with the power spectrum at wavenumbers $\log_{10}(k/{\rm km^{-1}\,s})\lesssim -1.5$. We use the Sherwood-Relics suite of hybrid radiation hydrodynamical simulations to perform a first analysis of new Lyman-$α$ forest power spectrum measurements at $4.0\leq z \leq 4.6$. These data extend to wavenumbers $\log_{10}(k/{\rm km^{-1}\,s})\simeq -3$, with a relative uncertainty of $10$--$20$ per cent in each wavenumber bin. Our analysis returns a $2.7σ$ preference for an enhancement in the Lyman-$α$ forest power spectrum at large scales, in excess of that expected for a spatially uniform ultraviolet background. This large-scale enhancement could be a signature of inhomogeneous reionization, although the statistical precision of these data is not yet sufficient for obtaining a robust detection of the relic post-reionization fluctuations. We show that future power spectrum measurements with relative uncertainties of $\lesssim 2.5$ per cent should provide unambiguous evidence for an enhancement in the power spectrum on large scales. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: Accepted by MNRAS, 13 pages, 8 figures

Enhanced C$_2$H$_2$ absorption within Jupiter's southern auroral oval from Juno UVS observations

Authors: Rohini S. Giles, Vincent Hue, Thomas K. Greathouse, G. Randall Gladstone, Joshua A. Kammer, Maarten H. Versteeg, Bertrand Bonfond, Denis G. Grodent, Jean-Claude Gérard, James A. Sinclair, Scott J. Bolton, Steven M. Levin

Abstract: Reflected sunlight observations from the Ultraviolet Spectrograph (UVS) on the Juno spacecraft were used to study the distribution of acetylene (C$_2$H$_2$) at Jupiter's south pole. We find that the shape of the C$_2$H$_2$ absorption feature varies significantly across the polar region, and this can be used to infer spatial variability in the C$_2$H$_2$ abundance. There is a localized region of en… ▽ More Reflected sunlight observations from the Ultraviolet Spectrograph (UVS) on the Juno spacecraft were used to study the distribution of acetylene (C$_2$H$_2$) at Jupiter's south pole. We find that the shape of the C$_2$H$_2$ absorption feature varies significantly across the polar region, and this can be used to infer spatial variability in the C$_2$H$_2$ abundance. There is a localized region of enhanced C$_2$H$_2$ absorption which coincides with the location of Jupiter's southern polar aurora; the C$_2$H$_2$ abundance poleward of the auroral oval is a factor of 3 higher than adjacent quiescent, non-auroral longitudes. This builds on previous infrared studies which found enhanced C$_2$H$_2$ abundances within the northern auroral oval. This suggests that Jupiter's upper-atmosphere chemistry is being strongly influenced by the influx of charged auroral particles and demonstrates the necessity of developing ion-neutral photochemical models of Jupiter's polar regions. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: Accepted in JGR: Planets

Jupiter's interior from Juno: Equation-of-state uncertainties and dilute core extent

Authors: S. Howard, T. Guillot, M. Bazot, Y. Miguel, D. J. Stevenson, E. Galanti, Y. Kaspi, W. B. Hubbard, B. Militzer, R. Helled, N. Nettelmann, B. Idini, S. Bolton

Abstract: The Juno mission has provided measurements of Jupiter s gravity field with an outstanding level of accuracy, leading to better constraints on the interior of the planet. Improving our knowledge of the internal structure of Jupiter is key to understanding its formation and evolution but is also important in the framework of exoplanet exploration. In this study, we investigated the differences betwe… ▽ More The Juno mission has provided measurements of Jupiter s gravity field with an outstanding level of accuracy, leading to better constraints on the interior of the planet. Improving our knowledge of the internal structure of Jupiter is key to understanding its formation and evolution but is also important in the framework of exoplanet exploration. In this study, we investigated the differences between the state-of-the-art equations of state and their impact on the properties of interior models. Accounting for uncertainty on the hydrogen and helium equation of state, we assessed the span of the interior features of Jupiter. We carried out an extensive exploration of the parameter space and studied a wide range of interior models using Markov chain Monte Carlo (MCMC) simulations. To consider the uncertainty on the equation of state, we allowed for modifications of the equation of state in our calculations. Our models harbour a dilute core and indicate that Jupiter s internal entropy is higher than what is usually assumed from the Galileo probe measurements. We obtain solutions with extended dilute cores, but contrary to other recent interior models of Jupiter, we also obtain models with small dilute cores. The dilute cores in such solutions extend to 20% of Jupiter s mass, leading to better agreement with formation evolution models. We conclude that the equations of state used in Jupiter models have a crucial effect on the inferred structure and composition. Further explorations of the behaviour of hydrogen helium mixtures at the pressure and temperature conditions in Jupiter will help to constrain the interior of the planet, and therefore its origin. △ Less

Submitted 17 February, 2023; originally announced February 2023.

Journal ref: A&A 672, A33 (2023)

Astro Data Lab Spectral Viewer Requirements for Wide-Area Spectroscopic Surveys

Authors: Leah M. Fulmer, Stephanie Juneau, Catherine Merrill, Adam S. Bolton, David L. Nidever, Robert Nikutta, Stephen T. Ridgway, Knut A. G. Olsen, Benjamin A. Weaver

Abstract: The Astro Data Lab is preparing to host large spectroscopic datasets such as a copy of the Dark Energy Spectroscopic Instrument (DESI) survey, which is projected to include approximately 40 million spectra of galaxies and quasars as well as over 10 million spectra of stars by 2026. Currently, we serve DR16 spectra from the Sloan Digital Sky Survey (SDSS), including Baryon Oscillation Spectroscopic… ▽ More The Astro Data Lab is preparing to host large spectroscopic datasets such as a copy of the Dark Energy Spectroscopic Instrument (DESI) survey, which is projected to include approximately 40 million spectra of galaxies and quasars as well as over 10 million spectra of stars by 2026. Currently, we serve DR16 spectra from the Sloan Digital Sky Survey (SDSS), including Baryon Oscillation Spectroscopic Survey (BOSS), and Extended BOSS (eBOSS) spectra. A spectral viewer tool allows users to visually and interactively inspect spectra. Given the large size of these spectroscopic datasets, a typical use case might consist of a selection or query for a subset of objects of interest (e.g., a subsample of stars or galaxies or quasars), followed by visual inspection of the selected spectra. It is anticipated that in some cases, users will want to go through a long list of spectra (e.g., thousands) quickly while looking for specific features. This document contains a description of the requirements for such a spectral viewer tool to be incorporated within the Astro Data Lab environment at NSF's NOIRLab. For each object, the spectral viewer will display the observed spectrum and, if available, the noise spectrum, sky spectrum, and best-fit template spectrum. Users will be able to control the display interactively after they launch the tool as part of their Data Lab workflow. The primary objective will be to support the visualization of spectroscopic datasets hosted at the Astro Data Lab but this requirements document could be a useful reference or inspiration for other applications and/or other datasets in the astronomy community. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: Spectral Viewer requirements document for astrophysical software, 27 pages

Effect of a magnetosphere compression on Jovian radio emissions: in situ case study using Juno data

Authors: C. K. Louis, C. M. Jackman, G. Hospodarsky, A. O'Kane Hackett, E. Devon-Hurley, P. Zarka, W. S. Kurth, R. W. Ebert, D. M. Weigt, A. R. Fogg, J. E. Waters, S. Mc Entee, J. E. P. Connerney, P. Louarn, S. Levin, S. J. Bolton

Abstract: During its 53-day polar orbit around Jupiter, Juno often crosses the boundaries of the Jovian magnetosphere (namely the magnetopause and bow shock). From the boundary locations, the upstream solar wind dynamic pressure can be inferred, which in turn illustrates the state of compression or relaxation of the system. The aim of this study is to examine Jovian radio emissions during magnetospheric com… ▽ More During its 53-day polar orbit around Jupiter, Juno often crosses the boundaries of the Jovian magnetosphere (namely the magnetopause and bow shock). From the boundary locations, the upstream solar wind dynamic pressure can be inferred, which in turn illustrates the state of compression or relaxation of the system. The aim of this study is to examine Jovian radio emissions during magnetospheric compressions, in order to determine the relationship between the solar wind and Jovian radio emissions. In this paper, we give a complete list of bow shock and magnetopause crossings (from June 2016 to August 2022), along with some extra informations (e.g. solar wind dynamic pressure and position of the standoff distances inferred from Joy et al. (2002)). We then select two compression events that occur in succession (inferred from magnetopause crossings) and we present a case study of the response of the Jovian radio emissions. We demonstrate that magnetospheric compressions lead to the activation of new radio sources. Newly activated broadband kilometric emissions are observed almost simultaneously to compression of the magnetosphere, with sources covering a large range of longitudes. Decametric emission sources are seen to be activated more than one rotation later only at specific longitudes and dusk local times. Finally, the activation of narrowband kilometric radiation is not observed during the compression phase, but when the magnetosphere is in its expansion phase. △ Less

Submitted 10 August, 2023; v1 submitted 7 December, 2022; originally announced December 2022.

Planetary Exploration Horizon 2061 Report, Chapter 4: From planetary exploration goals to technology requirements

Authors: Jérémie Lasue, Pierre Bousquet, Michel Blanc, Nicolas André, Pierre Beck, Gilles Berger, Scott Bolton, Emma Bunce, Baptiste Chide, Bernard Foing, Heidi Hammel, Emmanuel Lellouch, Lea Griton, Ralph Mcnutt, Sylvestre Maurice, Olivier Mousis, Merav Opher, Christophe Sotin, Dave Senske, Linda Spilker, Pierre Vernazza, Qiugang Zong

Abstract: This chapter reviews for each province and destination of the Solar System the representative space missions that will have to be designed and implemented by 2061 to address the six key science questions about the diversity, origins, workings and habitability of planetary systems (described in chapter 1) and to perform the critical observations that have been described in chapters 3 and partly 2.… ▽ More This chapter reviews for each province and destination of the Solar System the representative space missions that will have to be designed and implemented by 2061 to address the six key science questions about the diversity, origins, workings and habitability of planetary systems (described in chapter 1) and to perform the critical observations that have been described in chapters 3 and partly 2. It derives from this set of future representative missions, some of which will have to be flown during the 2041-2061 period, the critical technologies and supporting infrastructures that will be needed to fly these challenging missions, thus laying the foundation for the description of technologies and infrastructures for the future of planetary exploration that is given in chapters 5 and 6, respectively. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: 83 pages, 23 figures, 7 tables, Horizon 2061 is a science-driven, foresight exercise, for future scientific investigations

Reconstructing and Classifying SDSS DR16 Galaxy Spectra with Machine-Learning and Dimensionality Reduction Algorithms

Authors: Felix Pat, Stéphanie Juneau, Vanessa Böhm, Ragadeepika Pucha, A. G. Kim, A. S. Bolton, Cleo Lepart, Dylan Green, Adam D. Myers

Abstract: Optical spectra of galaxies and quasars from large cosmological surveys are used to measure redshifts and infer distances. They are also rich with information on the intrinsic properties of these astronomical objects. However, their physical interpretation can be challenging due to the substantial number of degrees of freedom, various sources of noise, and degeneracies between physical parameters… ▽ More Optical spectra of galaxies and quasars from large cosmological surveys are used to measure redshifts and infer distances. They are also rich with information on the intrinsic properties of these astronomical objects. However, their physical interpretation can be challenging due to the substantial number of degrees of freedom, various sources of noise, and degeneracies between physical parameters that cause similar spectral characteristics. To gain deeper insights into these degeneracies, we apply two unsupervised machine learning frameworks to a sample from the Sloan Digital Sky Survey data release 16 (SDSS DR16). The first framework is a Probabilistic Auto-Encoder (PAE), a two-stage deep learning framework consisting of a data compression stage from 1000 elements to 10 parameters and a density estimation stage. The second framework is a Uniform Manifold Approximation and Projection (UMAP), which we apply to both the uncompressed and compressed data. Exploring across regions on the compressed data UMAP, we construct sequences of stacked spectra which show a gradual transition from star-forming galaxies with narrow emission lines and blue spectra to passive galaxies with absorption lines and red spectra. Focusing on galaxies with broad emission lines produced by quasars, we find a sequence with varying levels of obscuration caused by cosmic dust. The experiments we present here inform future applications of neural networks and dimensionality reduction algorithms for large astronomical spectroscopic surveys. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: ASP Conference Series, Compendium of Undergraduate Research in Astronomy and Space Science (accepted), 24 pages, 14 figures

Probing quasar lifetimes with proximate $21$-centimetre absorption in the diffuse intergalactic medium at redshifts $z\geq 6$

Authors: Tomáš Šoltinský, James S. Bolton, Margherita Molaro, Nina Hatch, Martin G. Haehnelt, Laura C. Keating, Girish Kulkarni, Ewald Puchwein

Abstract: Enhanced ionizing radiation in close proximity to redshift $z\gtrsim 6$ quasars creates short windows of intergalactic Ly$α$ transmission blueward of the quasar Ly$α$ emission lines. The majority of these Ly$α$ near-zones are consistent with quasars that have optically/UV bright lifetimes of $t_{\rm Q}\sim 10^{5}-10^{7}\rm\,yr$. However, lifetimes as short as $t_{\rm Q}\lesssim 10^{4}\rm\,yr$ appe… ▽ More Enhanced ionizing radiation in close proximity to redshift $z\gtrsim 6$ quasars creates short windows of intergalactic Ly$α$ transmission blueward of the quasar Ly$α$ emission lines. The majority of these Ly$α$ near-zones are consistent with quasars that have optically/UV bright lifetimes of $t_{\rm Q}\sim 10^{5}-10^{7}\rm\,yr$. However, lifetimes as short as $t_{\rm Q}\lesssim 10^{4}\rm\,yr$ appear to be required by the smallest Ly$α$ near-zones. These short lifetimes present an apparent challenge for the growth of $\sim 10^{9}\rm\,M_{\odot}$ black holes at $z\gtrsim 6$. Accretion over longer timescales is only possible if black holes grow primarily in an obscured phase, or if the quasars are variable on timescales comparable to the equilibriation time for ionized hydrogen. Distinguishing between very young quasars and older quasars that have experienced episodic accretion with Ly$α$ absorption alone is challenging, however. We therefore predict the signature of proximate 21-cm absorption around $z\gtrsim 6$ radio-loud quasars. For modest pre-heating of intergalactic hydrogen by the X-ray background, where the spin temperature $T_{\rm S} \lesssim 10^{2}\rm\,K$ prior to any quasar heating, we find proximate 21-cm absorption should be observable in the spectra of radio-loud quasars. The extent of the proximate 21-cm absorption is sensitive to the integrated lifetime of the quasar. Evidence for proximate 21-cm absorption from the diffuse intergalactic medium within $2-3\rm\,pMpc$ of a (radio-loud) quasar would be consistent with a short quasar lifetime, $t_{\rm Q}\lesssim 10^{5}\rm\,yr$, and would provide a complementary constraint on models for high redshift black hole growth. △ Less

Submitted 19 December, 2022; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 19 pages, 12 figures, 2 tables. Accepted for publication in MNRAS

LensWatch: I. Resolved HST Observations and Constraints on the Strongly-Lensed Type Ia Supernova 2022qmx ("SN Zwicky")

Authors: J. D. R. Pierel, N. Arendse, S. Ertl, X. Huang, L. A. Moustakas, S. Schuldt, A. J. Shajib, Y. Shu, S. Birrer, M. Bronikowski, J. Hjorth, S. H. Suyu, S. Agarwal, A. Agnello, A. S. Bolton, S. Chakrabarti, C. Cold, F. Courbin, J. M. Della Costa, S. Dhawan, M. Engesser, O. D. Fox, C. Gall, S. Gomez, A. Goobar , et al. (17 additional authors not shown)

Abstract: Supernovae (SNe) that have been multiply-imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply-image… ▽ More Supernovae (SNe) that have been multiply-imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply-imaged Type Ia SN 2022qmx (aka, "SN Zwicky"; Goobar et al. 2022) at z = 0.3544. SN Zwicky was discovered by the Zwicky Transient Facility (ZTF) in spatially unresolved data. Here we present follow-up Hubble Space Telescope observations of SN Zwicky, the first from the multi-cycle "LensWatch" program (www.lenswatch.org). We measure photometry for each of the four images of SN Zwicky, which are resolved in three WFC3/UVIS filters (F475W, F625W, F814W) but unresolved with WFC3/IR F160W, and produce an analysis of the lensing system using a variety of independent lens modeling methods. We find consistency between time delays estimated with the single epoch of HST photometry and the lens model predictions constrained through the multiple image positions, with both inferring time delays of <1 day. Our lens models converge to an Einstein radius of (0.168+0.009-0.005)", the smallest yet seen in a lensed SN. The "standard candle" nature of SN Zwicky provides magnification estimates independent of the lens modeling that are brighter by ~1.5 mag and ~0.8 mag for two of the four images, suggesting significant microlensing and/or additional substructure beyond the flexibility of our image-position mass models. △ Less

Submitted 8 November, 2022; v1 submitted 7 November, 2022; originally announced November 2022.

The Spectroscopic Data Processing Pipeline for the Dark Energy Spectroscopic Instrument

Authors: J. Guy, S. Bailey, A. Kremin, Shadab Alam, D. M. Alexander, C. Allende Prieto, S. BenZvi, A. S. Bolton, D. Brooks, E. Chaussidon, A. P. Cooper, K. Dawson, A. de la Macorra, A. Dey, Biprateep Dey, G. Dhungana, D. J. Eisenstein, A. Font-Ribera, J. E. Forero-Romero, E. Gaztañaga, S. Gontcho A Gontcho, D. Green, K. Honscheid, M. Ishak, R. Kehoe , et al. (33 additional authors not shown)

Abstract: We describe the spectroscopic data processing pipeline of the Dark Energy Spectroscopic Instrument (DESI), which is conducting a redshift survey of about 40 million galaxies and quasars using a purpose-built instrument on the 4-m Mayall Telescope at Kitt Peak National Observatory. The main goal of DESI is to measure with unprecedented precision the expansion history of the Universe with the Baryon… ▽ More We describe the spectroscopic data processing pipeline of the Dark Energy Spectroscopic Instrument (DESI), which is conducting a redshift survey of about 40 million galaxies and quasars using a purpose-built instrument on the 4-m Mayall Telescope at Kitt Peak National Observatory. The main goal of DESI is to measure with unprecedented precision the expansion history of the Universe with the Baryon Acoustic Oscillation technique and the growth rate of structure with Redshift Space Distortions. Ten spectrographs with three cameras each disperse the light from 5000 fibers onto 30 CCDs, covering the near UV to near infrared (3600 to 9800 Angstrom) with a spectral resolution ranging from 2000 to 5000. The DESI data pipeline generates wavelength- and flux-calibrated spectra of all the targets, along with spectroscopic classifications and redshift measurements. Fully processed data from each night are typically available to the DESI collaboration the following morning. We give details about the pipeline's algorithms, and provide performance results on the stability of the optics, the quality of the sky background subtraction, and the precision and accuracy of the instrumental calibration. This pipeline has been used to process the DESI Survey Validation data set, and has exceeded the project's requirements for redshift performance, with high efficiency and a purity greater than 99 percent for all target classes. △ Less

Submitted 9 January, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: AJ, revised version, 55 pages, 55 figures, 4 tables

Observation of Turbulent Magnetohydrodynamic Cascade in the Jovian Magnetosheath

Authors: N. Andrés, R. Bandyopadhyay, D. J. McComas, J. R. Szalay, F. Allegrini, R. W. Ebert, D. J. Gershman, J. E. P. Connerney, S. J. Bolton

Abstract: We present the first estimation of the energy cascade rate in Jupiter's magnetosheath (MS). We use in-situ observations from the Jovian Auroral Distributions Experiment (JADE) and the magnetometer investigation (MAG) instruments onboard the Juno spacecraft, in concert with two recent compressible models to investigate the cascade rate in the magnetohydrodynamic (MHD) scales. While a high level of… ▽ More We present the first estimation of the energy cascade rate in Jupiter's magnetosheath (MS). We use in-situ observations from the Jovian Auroral Distributions Experiment (JADE) and the magnetometer investigation (MAG) instruments onboard the Juno spacecraft, in concert with two recent compressible models to investigate the cascade rate in the magnetohydrodynamic (MHD) scales. While a high level of compressible density fluctuations is observed in the Jovian MS, a constant energy flux exists in the MHD inertial range. The compressible isothermal and polytropic energy cascade rates increase in the MHD range when density fluctuations are present. We find that the energy cascade rate in Jupiter's magnetosheath is at least two orders of magnitude (100 times) smaller than the corresponding typical value in the Earth's magnetosheath. △ Less

Submitted 31 January, 2023; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: To be published in The Astrophysical Journal

The MegaMapper: A Stage-5 Spectroscopic Instrument Concept for the Study of Inflation and Dark Energy

Authors: David J. Schlegel, Juna A. Kollmeier, Greg Aldering, Stephen Bailey, Charles Baltay, Christopher Bebek, Segev BenZvi, Robert Besuner, Guillermo Blanc, Adam S. Bolton, Ana Bonaca, Mohamed Bouri, David Brooks, Elizabeth Buckley-Geer, Zheng Cai, Jeffrey Crane, Regina Demina, Joseph DeRose, Arjun Dey, Peter Doel, Xiaohui Fan, Simone Ferraro, Douglas Finkbeiner, Andreu Font-Ribera, Satya Gontcho A Gontcho , et al. (64 additional authors not shown)

Abstract: In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2<z<5$. In order to achieve path-breaking results with a mid-scale investment, the MegaMapper combines existing technologies for critical path elements and pushes innovative development in other design areas. To this… ▽ More In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2<z<5$. In order to achieve path-breaking results with a mid-scale investment, the MegaMapper combines existing technologies for critical path elements and pushes innovative development in other design areas. To this aim, we envision a 6.5-m Magellan-like telescope, with a newly designed wide field, coupled with DESI spectrographs, and small-pitch robots to achieve multiplexing of at least 26,000. This will match the expected achievable target density in the redshift range of interest and provide a 10x capability over the existing state-of the art, without a 10x increase in project budget. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Comments: Contributed White Paper to Snowmass 2021. arXiv admin note: substantial text overlap with arXiv:1907.11171. text overlap with arXiv:2209.03585

A Spectroscopic Road Map for Cosmic Frontier: DESI, DESI-II, Stage-5

Authors: David J. Schlegel, Simone Ferraro, Greg Aldering, Charles Baltay, Segev BenZvi, Robert Besuner, Guillermo A. Blanc, Adam S. Bolton, Ana Bonaca, David Brooks, Elizabeth Buckley-Geer, Zheng Cai, Joseph DeRose, Arjun Dey, Peter Doel, Alex Drlica-Wagner, Xiaohui Fan, Gaston Gutierrez, Daniel Green, Julien Guy, Dragan Huterer, Leopoldo Infante, Patrick Jelinsky, Dionysios Karagiannis, Stephen M. Kent , et al. (40 additional authors not shown)

Abstract: In this white paper, we present an experimental road map for spectroscopic experiments beyond DESI. DESI will be a transformative cosmological survey in the 2020s, mapping 40 million galaxies and quasars and capturing a significant fraction of the available linear modes up to z=1.2. DESI-II will pilot observations of galaxies both at much higher densities and extending to higher redshifts. A Stage… ▽ More In this white paper, we present an experimental road map for spectroscopic experiments beyond DESI. DESI will be a transformative cosmological survey in the 2020s, mapping 40 million galaxies and quasars and capturing a significant fraction of the available linear modes up to z=1.2. DESI-II will pilot observations of galaxies both at much higher densities and extending to higher redshifts. A Stage-5 experiment would build out those high-density and high-redshift observations, mapping hundreds of millions of stars and galaxies in three dimensions, to address the problems of inflation, dark energy, light relativistic species, and dark matter. These spectroscopic data will also complement the next generation of weak lensing, line intensity mapping and CMB experiments and allow them to reach their full potential. △ Less

Submitted 8 September, 2022; originally announced September 2022.

Comments: Contribution to Snowmass 2021

Alternating north-south brightness ratio of Ganymede's auroral ovals: Hubble Space Telescope observations around the Juno PJ34 flyby

Authors: Joachim Saur, Stefan Duling, Alexandre Wennmacher, Clarissa Willmes, Lorenz Roth, Darrell F. Strobel, Frédéric Allegrini, Fran Bagenal, Scott J. Bolton, Bertrand Bonfond, George Clark, Randy Gladstone, T. K. Greathouse, Denis C. Grodent, Candice J. Hansen, W. S. Kurth, Glenn S. Orton, Kurt D. Retherford, Abigail M. Rymer, Ali H. Sulaiman

Abstract: We report results of Hubble Space Telescope observations from Ganymede's orbitally trailing side which were taken around the flyby of the Juno spacecraft on June 7, 2021. We find that Ganymede's northern and southern auroral ovals alternate in brightness such that the oval facing Jupiter's magnetospheric plasma sheet is brighter than the other one. This suggests that the generator that powers Gany… ▽ More We report results of Hubble Space Telescope observations from Ganymede's orbitally trailing side which were taken around the flyby of the Juno spacecraft on June 7, 2021. We find that Ganymede's northern and southern auroral ovals alternate in brightness such that the oval facing Jupiter's magnetospheric plasma sheet is brighter than the other one. This suggests that the generator that powers Ganymede's aurora is the momentum of the Jovian plasma sheet north and south of Ganymede's magnetosphere. Magnetic coupling of Ganymede to the plasma sheet above and below the moon causes asymmetric magnetic stresses and electromagnetic energy fluxes ultimately powering the auroral acceleration process. No clear statistically significant time variability of the auroral emission on short time scales of 100s could be resolved. We show that electron energy fluxes of several tens of mW m$^{-2}$ are required for its OI 1356 Å$\;$ emission making Ganymede a very poor auroral emitter. △ Less

Submitted 18 August, 2022; originally announced August 2022.

Comments: Accepted for publication with Geophys. Res. Let

The Sherwood-Relics simulations: overview and impact of patchy reionization and pressure smoothing on the intergalactic medium

Authors: Ewald Puchwein, James S. Bolton, Laura C. Keating, Margherita Molaro, Prakash Gaikwad, Girish Kulkarni, Martin G. Haehnelt, Vid Iršič, Tomáš Šoltinský, Matteo Viel, Dominique Aubert, George D. Becker, Avery Meiksin

Abstract: We present the Sherwood-Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for trea… ▽ More We present the Sherwood-Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for treating radiative transfer effects. This scheme follows the spatial variations in the ionizing radiation field, as well as the associated fluctuations in IGM temperature and pressure smoothing. It is computationally much cheaper than full radiation hydrodynamics simulations and circumvents the difficult task of calibrating a galaxy formation model to observational constraints on cosmic reionization. Using this hybrid technique, we study the spatial fluctuations in IGM properties that are seeded by patchy cosmic reionization. We investigate the relevant physical processes and assess their impact on the z > 4 Lyman-alpha forest. Our main findings are: (i) Consistent with previous studies patchy reionization causes large scale temperature fluctuations that persist well after the end of reionization, (ii) these increase the Lyman-alpha forest flux power spectrum on large scales, and (iii) result in a spatially varying pressure smoothing that correlates well with the local reionization redshift. (iv) Structures evaporated or puffed up by photoheating cause notable features in the Lyman-alpha forest, such as flat-bottom or double-dip absorption profiles. △ Less

Submitted 28 February, 2023; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: Version accepted for publication in MNRAS, some discussions added, data added to Fig. 8. 24 pages, 15 figures. Images and movies may be viewed at https://www.nottingham.ac.uk/astronomy/sherwood-relics/

Hints of dark photon dark matter from observations and hydrodynamical simulations of the low-redshift Lyman-$α$ forest

Authors: James S. Bolton, Andrea Caputo, Hongwan Liu, Matteo Viel

Abstract: Recent work has suggested that an additional $\lesssim 6.9 \rm{\, eV}$ per baryon of heating in the intergalactic medium is needed to reconcile hydrodynamical simulations with Lyman-$α$ forest absorption line widths at redshift $z\simeq 0.1$. Resonant conversion of dark photon dark matter into low frequency photons is a viable source of such heating. We perform the first hydrodynamical simulations… ▽ More Recent work has suggested that an additional $\lesssim 6.9 \rm{\, eV}$ per baryon of heating in the intergalactic medium is needed to reconcile hydrodynamical simulations with Lyman-$α$ forest absorption line widths at redshift $z\simeq 0.1$. Resonant conversion of dark photon dark matter into low frequency photons is a viable source of such heating. We perform the first hydrodynamical simulations including dark photon heating and show that dark photons with mass $m_{A'}\sim 8\times 10^{-14}\rm\,eV\,c^{-2}$ and kinetic mixing $ε\sim 5\times 10^{-15}$ can alleviate the heating excess. A prediction of this model is a non-standard thermal history for underdense gas at $z \gtrsim 3$. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 5+1 pages, 2+1 figures

Jupiter's Temperature Structure: A Reassessment of the Voyager Radio Occultation Measurements

Authors: Pranika Gupta, Sushil K. Atreya, Paul G. Steffes, Leigh N. Fletcher, Tristan Guillot, Michael D. Allison, Scott J. Bolton, Ravit Helled, Steven Levin, Cheng Li, Jonathan I. Lunine, Yamila Miguel, Glenn S. Orton, J. Hunter Waite, Paul Withers

Abstract: The thermal structure of planetary atmospheres is an essential input for predicting and retrieving the distribution of gases and aerosols, as well as the bulk chemical abundances. In the case of Jupiter, the temperature at a reference level - generally taken at 1 bar - serves as the anchor in models used to derive the planet's interior structure and composition. Most models assume the temperature… ▽ More The thermal structure of planetary atmospheres is an essential input for predicting and retrieving the distribution of gases and aerosols, as well as the bulk chemical abundances. In the case of Jupiter, the temperature at a reference level - generally taken at 1 bar - serves as the anchor in models used to derive the planet's interior structure and composition. Most models assume the temperature measured by the Galileo probe (Seiff et al. 1998). However, those data correspond to a single location, an unusually clear, dry region, affected by local atmospheric dynamics. On the other hand, the Voyager radio occultation observations cover a wider range of latitudes, longitudes, and times (Lindal et al. 1981). The Voyager retrievals were based on atmospheric composition and radio refractivity data that require updating and were never properly tabulated: the few existing tabulations are incomplete and ambiguous. Here, we present a systematic electronic digitization of all available temperature profiles from Voyager, followed by their reanalysis, employing currently accepted values of the abundances and radio refractivities of atmospheric species. We find the corrected temperature at the 1 bar level to be up to 4 K greater than previously published values, i.e., 170.3{\pm}3.8 K at 12°S (Voyager 1 ingress) and 167.3{\pm}3.8 K at 0°N (Voyager 1 egress). This is to be compared with the Galileo probe value of 166.1{\pm}0.8 K at the edge of an unusual feature at 6.57°N. Altogether, this suggests that Jupiter's tropospheric temperatures may vary spatially by up to 7 K between 7°N and 12°S. △ Less

Submitted 14 July, 2022; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: 11 pages, 3 figures

Journal ref: Planet. Sci. J. 3 159 (2022)

Jupiter's inhomogeneous envelope

Authors: Y. Miguel, M. Bazot, T. Guillot, S. Howard, E. Galanti, Y. Kaspi, W. B. Hubbard, B. Militzer, R. Helled, S. K. Atreya, J. E. P. Connerney, D. Durante, L. Kulowski, J. I. Lunine, D. Stevenson, S. Bolton

Abstract: While Jupiter's massive gas envelope consists mainly of hydrogen and helium, the key to understanding Jupiter's formation and evolution lies in the distribution of the remaining (heavy) elements. Before the Juno mission, the lack of high-precision gravity harmonics precluded the use of statistical analyses in a robust determination of the heavy-elements distribution in Jupiter's envelope. In this… ▽ More While Jupiter's massive gas envelope consists mainly of hydrogen and helium, the key to understanding Jupiter's formation and evolution lies in the distribution of the remaining (heavy) elements. Before the Juno mission, the lack of high-precision gravity harmonics precluded the use of statistical analyses in a robust determination of the heavy-elements distribution in Jupiter's envelope. In this paper, we assemble the most comprehensive and diverse collection of Jupiter interior models to date and use it to study the distribution of heavy elements in the planet's envelope. We apply a Bayesian statistical approach to our interior model calculations, reproducing the Juno gravitational and atmospheric measurements and constraints from the deep zonal flows. Our results show that the gravity constraints lead to a deep entropy of Jupiter corresponding to a 1 bar temperature 5-15 K higher than traditionally assumed. We also find that uncertainties in the equation of state are crucial when determining the amount of heavy elements in Jupiter's interior. Our models put an upper limit to the inner compact core of Jupiter of 7 Earth masses, independently on the structure model (with or without dilute core) and the equation of state considered. Furthermore, we robustly demonstrate that Jupiter's envelope is inhomogenous, with a heavy-element enrichment in the interior relative to the outer envelope. This implies that heavy element enrichment continued through the gas accretion phase, with important implications for the formation of giant planets in our solar system and beyond. △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 662, A18 (2022)

Revelations on Jupiter's Formation, Evolution and Interior: Challenges from Juno Results

Authors: Ravit Helled, David J. Stevenson, Jonathan I. Lunine, Scott J. Bolton, Nadine Nettelmann, Sushil Atreya, Tristan Guillot, Burkhard Militzer, Yamila Miguel, William B. Hubbard

Abstract: The Juno mission has revolutionized and challenged our understanding of Jupiter. As Juno transitioned to its extended mission, we review the major findings of Jupiter's internal structure relevant to understanding Jupiter's formation and evolution. Results from Juno's investigation of Jupiter's interior structure imply that the planet has compositional gradients and is accordingly non-adiabatic, w… ▽ More The Juno mission has revolutionized and challenged our understanding of Jupiter. As Juno transitioned to its extended mission, we review the major findings of Jupiter's internal structure relevant to understanding Jupiter's formation and evolution. Results from Juno's investigation of Jupiter's interior structure imply that the planet has compositional gradients and is accordingly non-adiabatic, with a complex internal structure. These new results imply that current models of Jupiter's formation and evolution require a revision. In this paper, we discuss potential formation and evolution paths that can lead to an internal structure model consistent with Juno data, and the constraints they provide. We note that standard core accretion formation models, including the heavy-element enrichment during planetary growth is consistent with an interior that is inhomogeneous with composition gradients in its deep interior. However, such formation models typically predict that this region, which could be interpreted as a primordial dilute core, is confined to about 10% of Jupiter's total mass. In contrast, structure models that fit Juno data imply that this region contains 30% of the mass or more. One way to explain the origin of this extended region is by invoking a relatively long (about 2 Myrs) formation phase where the growing planet accretes gas and planetesimals delaying the runaway gas accretion. Alternatively, Jupiter's fuzzy core could be a result of a giant impact or convection post-formation. These novel scenarios require somewhat special and specific conditions. Clarity on the plausibility of such conditions could come from future high-resolution observations of planet-forming regions around other stars, from the observed and modeled architectures of extrasolar systems with giant planets, and future Juno data obtained during its extended mission. △ Less

Submitted 21 February, 2022; originally announced February 2022.

Comments: invited review for Icarus, accepted for publication

GIGA-Lens: Fast Bayesian Inference for Strong Gravitational Lens Modeling

Authors: A. Gu, X. Huang, W. Sheu, G. Aldering, A. S. Bolton, K. Boone, A. Dey, A. Filipp, E. Jullo, S. Perlmutter, D. Rubin, E. F. Schlafly, D. J. Schlegel, Y. Shu, S. H. Suyu

Abstract: We present GIGA-Lens: a gradient-informed, GPU-accelerated Bayesian framework for modeling strong gravitational lensing systems, implemented in TensorFlow and JAX. The three components, optimization using multi-start gradient descent, posterior covariance estimation with variational inference, and sampling via Hamiltonian Monte Carlo, all take advantage of gradient information through automatic di… ▽ More We present GIGA-Lens: a gradient-informed, GPU-accelerated Bayesian framework for modeling strong gravitational lensing systems, implemented in TensorFlow and JAX. The three components, optimization using multi-start gradient descent, posterior covariance estimation with variational inference, and sampling via Hamiltonian Monte Carlo, all take advantage of gradient information through automatic differentiation and massive parallelization on graphics processing units (GPUs). We test our pipeline on a large set of simulated systems and demonstrate in detail its high level of performance. The average time to model a single system on four Nvidia A100 GPUs is 105 seconds. The robustness, speed, and scalability offered by this framework make it possible to model the large number of strong lenses found in current surveys and present a very promising prospect for the modeling of $\mathcal{O}(10^5)$ lensing systems expected to be discovered in the era of the Vera C. Rubin Observatory, Euclid, and the Nancy Grace Roman Space Telescope. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 23 pages, 13 figures, 2 tables. Submitted to ApJ

The Seventeenth Data Release of the Sloan Digital Sky Surveys: Complete Release of MaNGA, MaStar and APOGEE-2 Data

Authors: Abdurro'uf, Katherine Accetta, Conny Aerts, Victor Silva Aguirre, Romina Ahumada, Nikhil Ajgaonkar, N. Filiz Ak, Shadab Alam, Carlos Allende Prieto, Andres Almeida, Friedrich Anders, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Erik Aquino-Ortiz, Alfonso Aragon-Salamanca, Maria Argudo-Fernandez, Metin Ata, Marie Aubert, Vladimir Avila-Reese, Carles Badenes, Rodolfo H. Barba, Kat Barger, Jorge K. Barrera-Ballesteros, Rachael L. Beaton , et al. (316 additional authors not shown)

Abstract: This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies… ▽ More This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys. △ Less

Submitted 13 January, 2022; v1 submitted 3 December, 2021; originally announced December 2021.

Comments: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

Limits on non-canonical heating and turbulence in the intergalactic medium from the low redshift Lyman-alpha forest

Authors: James S. Bolton, Prakash Gaikwad, Martin G. Haehnelt, Tae-Sun Kim, Fahad Nasir, Ewald Puchwein, Matteo Viel, Bart P. Wakker

Abstract: We examine the column density distribution function (CDDF) and Doppler parameter distribution from hydrodynamical simulations and Cosmic Origins Spectrograph (COS) observations of the Lyman-alpha forest at redshift $0\leq z\leq 0.2$. Allowing for a factor of two uncertainty in the metagalactic HI photoionisation rate, our hydrodynamical simulations are in good agreement ($1$-$1.5σ$) with the shape… ▽ More We examine the column density distribution function (CDDF) and Doppler parameter distribution from hydrodynamical simulations and Cosmic Origins Spectrograph (COS) observations of the Lyman-alpha forest at redshift $0\leq z\leq 0.2$. Allowing for a factor of two uncertainty in the metagalactic HI photoionisation rate, our hydrodynamical simulations are in good agreement ($1$-$1.5σ$) with the shape and amplitude of the observed CDDF at HI column densities $10^{13.3}\rm\,cm^{-2}\leq N_{\rm HI}\leq 10^{14.5}\rm\,cm^{-2}$. However, the Doppler widths of the simulated lines remain too narrow with respect to the COS data. We argue that invoking AGN feedback does not resolve this discrepancy. We also disfavour enhanced photoheating rates as a potential solution, as this requires an unphysically hard UV background spectrum. If instead appealing to a non-canonical source of heating, an additional specific heat injection of $u \lesssim 6.9\rm\,eV\,m_{\rm p}^{-1}$ is required at $z\lesssim 2.5$ for gas that has $N_{\rm HI}\simeq 10^{13.5}\rm\,cm^{-2}$ by $z=0.1$. Alternatively, there may be an unresolved line of sight turbulent velocity component of $v_{\rm turb}\lesssim 8.5\rm\,km\,s^{-1}(N_{\rm HI}/10^{13.5}\rm\,cm^{-2})^{0.21}$ for the coldest gas in the diffuse IGM. △ Less

Submitted 25 March, 2022; v1 submitted 18 November, 2021; originally announced November 2021.

Comments: 21 pages, 13 figures, MNRAS accepted following revision. The accepted version includes an extra panel in figure 2 and some additional clarifications on the line broadening within the models. Conclusions unchanged

Theory of Figures to the 7th order and the interiors of Jupiter and Saturn

Authors: N. Nettelmann, N. Movshovitz, D. Ni, J. J. Fortney, E. Galanti, Y. Kaspi, R. Helled, C. R. Mankovich, S. Bolton

Abstract: Interior modeling of Jupiter and Saturn has advanced to a state where thousands of models are generated that cover the uncertainty space of many parameters. This approach demands a fast method of computing their gravity field and shape. Moreover, the Cassini mission at Saturn and the ongoing Juno mission delivered gravitational harmonics up to J12. Here, we report the expansion of the Theory of Fi… ▽ More Interior modeling of Jupiter and Saturn has advanced to a state where thousands of models are generated that cover the uncertainty space of many parameters. This approach demands a fast method of computing their gravity field and shape. Moreover, the Cassini mission at Saturn and the ongoing Juno mission delivered gravitational harmonics up to J12. Here, we report the expansion of the Theory of Figures, which is a fast method for gravity field and shape computation, to the 7th-order (ToF7), which allows for computation of up to J14. We apply three different codes to compare the accuracy using polytropic models. We apply ToF7 to Jupiter and Saturn interior models in conjunction with CMS-19 H/He-EOS. For Jupiter, we find that J6 is best matched by a transition from He-depleted to He-enriched envelope at 2-2.5 Mbar. However, the atmospheric metallicity reaches 1xtimes solar only if the adiabat is perturbed toward lower densities, or if the surface temperature is enhanced by ~14 K from the Galileo value. Our Saturn models imply a largely homogeneous-in-Z envelope at 1.5-4x solar atop a small core. Perturbing the adiabat yields metallicity profiles with extended, heavy-element enriched deep interior (diffuse core) out to 0.4 RSat, as for Jupiter. Classical models with compact, dilute, or no core are possible as long as the deep interior is enriched in heavy-elements. Including a thermal wind fitted to the observed wind speeds, representative Jupiter and Saturn models are consistent with all observed Jn values. △ Less

Submitted 28 October, 2021; originally announced October 2021.

Comments: accepted to the Planetary Science Journal

Jupiter's Temperate Belt/Zone Contrasts Revealed at Depth by Juno Microwave Observations

Authors: L. N. Fletcher, F. A. Oyafuso, M. Allison, A. Ingersoll, L. Li, Y. Kaspi, E. Galanti, M. H. Wong, G. S. Orton, K. Duer, Z. Zhang, C. Li, T. Guillot, S. M. Levin, S. Bolton

Abstract: Juno Microwave Radiometer (MWR) observations of Jupiter's mid-latitudes reveal a strong correlation between brightness temperature contrasts and zonal winds, confirming that the banded structure extends throughout the troposphere. However, the microwave brightness gradient is observed to change sign with depth: the belts are microwave-bright in the $p<5$ bar range and microwave-dark in the $p>10$… ▽ More Juno Microwave Radiometer (MWR) observations of Jupiter's mid-latitudes reveal a strong correlation between brightness temperature contrasts and zonal winds, confirming that the banded structure extends throughout the troposphere. However, the microwave brightness gradient is observed to change sign with depth: the belts are microwave-bright in the $p<5$ bar range and microwave-dark in the $p>10$ bar range. The transition level (which we call the jovicline) is evident in the MWR 11.5 cm channel, which samples the 5-14 bar range when using the limb-darkening at all emission angles. The transition is located between 4 and 10 bars, and implies that belts change with depth from being NH$_3$-depleted to NH$_3$-enriched, or from physically-warm to physically-cool, or more likely a combination of both. The change in character occurs near the statically stable layer associated with water condensation. The implications of the transition are discussed in terms of ammonia redistribution via meridional circulation cells with opposing flows above and below the water condensation layer, and in terms of the `mushball' precipitation model, which predicts steeper vertical ammonia gradients in the belts versus the zones. We show via the moist thermal wind equation that both the temperature and ammonia interpretations can lead to vertical shear on the zonal winds, but the shear is $\sim50\times$ weaker if only NH$_3$ gradients are considered. Conversely, if MWR observations are associated with kinetic temperature gradients then it would produce zonal winds that increase in strength down to the jovicline, consistent with Galileo probe measurements; then decay slowly at higher pressures. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: Published in JGR: Planets

Journal ref: Journal of Geophysical Research: Planets, 126, e2021JE006858

Evidence for multiple Ferrel-like cells on Jupiter

Authors: Keren Duer, Nimrod Gavriel, Eli Galanti, Yohai Kaspi, Leigh N. Fletcher, Tristan Guillot, Scott J. Bolton, Steven M. Levin, Sushil K. Atreya, Davide Grassi, Andrew P. Ingersoll, Cheng Li, Liming Li, Jonathan I. Lunine, Glenn S. Orton, Fabiano A. Oyafuso, J. Hunter Waite Jr

Abstract: Jupiter's atmosphere is dominated by multiple jet streams which are strongly tied to its 3D atmospheric circulation. Lacking a rigid bottom boundary, several models exist for how the meridional circulation extends into the planetary interior. Here we show, collecting evidence from multiple instruments of the Juno mission, the existence of mid-latitudinal meridional circulation cells which are driv… ▽ More Jupiter's atmosphere is dominated by multiple jet streams which are strongly tied to its 3D atmospheric circulation. Lacking a rigid bottom boundary, several models exist for how the meridional circulation extends into the planetary interior. Here we show, collecting evidence from multiple instruments of the Juno mission, the existence of mid-latitudinal meridional circulation cells which are driven by turbulence, similar to the Ferrel cells on Earth. Different than Earth, which contains only one such cell in each hemisphere, the larger, faster rotating Jupiter can incorporate multiple cells. The cells form regions of upwelling and downwelling, which we show are clearly evident in Juno's microwave data between latitude 60S and 60N. The existence of these cells is confirmed by reproducing the ammonia observations using a simplistic model. This study solves a long-standing puzzle regarding the nature of Jupiter's sub-cloud dynamics and provides evidence for 8 cells in each Jovian hemisphere. △ Less

Submitted 14 October, 2021; originally announced October 2021.

Comments: 13 pages, 5 figures, Supporting information. Accepted to Geophysical Research Letters

The effect of inhomogeneous reionisation on the Lyman-$α$ forest power spectrum at redshift $z>4$: implications for thermal parameter recovery

Authors: Margherita Molaro, Vid Iršič, James S. Bolton, Laura C. Keating, Ewald Puchwein, Prakash Gaikwad, Martin G. Haehnelt, Girish Kulkarni, Matteo Viel

Abstract: We use the Sherwood-Relics suite of hybrid hydrodynamical and radiative transfer simulations to model the effect of inhomogeneous reionisation on the 1D power spectrum of the \Lya forest transmitted flux at redshifts $4.2\leq z \leq 5$. Relative to models that assume a homogeneous UV background, reionisation suppresses the power spectrum at small scales, $k \sim 0.1\rm\,km^{-1}\,s$ by $\sim 10$ pe… ▽ More We use the Sherwood-Relics suite of hybrid hydrodynamical and radiative transfer simulations to model the effect of inhomogeneous reionisation on the 1D power spectrum of the \Lya forest transmitted flux at redshifts $4.2\leq z \leq 5$. Relative to models that assume a homogeneous UV background, reionisation suppresses the power spectrum at small scales, $k \sim 0.1\rm\,km^{-1}\,s$ by $\sim 10$ per cent because of spatial variations in the thermal broadening kernel and the divergent peculiar velocity field associated with over-pressurised intergalactic gas. On larger scales, $k<0.03\rm\,km^{-1}\,s$, the power spectrum is instead enhanced by $10$--$50$ per cent by large scale spatial variations in the neutral hydrogen fraction. The effect of inhomogeneous reionisation must therefore be accounted for in analyses of forthcoming high precision measurements. We provide a correction for the Ly$α$ forest power spectrum at $4.1\leq z \leq 5.4$ that can be easily applied within other parameter inference frameworks using similar reionisation models. We perform a Bayesian analysis of mock data to assess the extent of systematic biases that may arise in measurements of the intergalactic medium if ignoring this correction. At the scales probed by current high resolution Ly$α$ forest data at $z>4$, $0.006 \rm \,km^{-1}\,s\leq k \leq 0.2 \rm\, km^{-1}\,s$, we find inhomogeneous reionisation does not introduce any significant bias in thermal parameter recovery for the current measurement uncertainties of $\sim 10$ per cent. However, for $5$ per cent uncertainties, $\sim 1σ$ shifts between the estimated and true parameters occur. △ Less

Submitted 5 January, 2022; v1 submitted 14 September, 2021; originally announced September 2021.

Comments: Accepted by MNRAS (21 November 2021), 20 pages, 11 figures

Journal ref: MNRAS, Volume 509, Issue 4, February 2022, Pages 6119-6137, Published 26 November 2021

Meridional variations of C$_2$H$_2$ in Jupiter's stratosphere from Juno UVS observations

Authors: Rohini S. Giles, Thomas K. Greathouse, Vincent Hue, G. Randall Gladstone, Henrik Melin, Leigh N. Fletcher, Patrick G. J. Irwin, Joshua A. Kammer, Maarten H. Versteeg, Bertrand Bonfond, Denis C. Grodent, Scott J. Bolton, Steven M. Levin

Abstract: The UVS instrument on the Juno mission records far-ultraviolet reflected sunlight from Jupiter. These spectra are sensitive to the abundances of chemical species in the upper atmosphere and to the distribution of the stratospheric haze layer. We combine observations from the first 30 perijoves of the mission in order to study the meridional distribution of acetylene (C$_2$H$_2$) in Jupiter's strat… ▽ More The UVS instrument on the Juno mission records far-ultraviolet reflected sunlight from Jupiter. These spectra are sensitive to the abundances of chemical species in the upper atmosphere and to the distribution of the stratospheric haze layer. We combine observations from the first 30 perijoves of the mission in order to study the meridional distribution of acetylene (C$_2$H$_2$) in Jupiter's stratosphere. We find that the abundance of C$_2$H$_2$ decreases towards the poles by a factor of 2-4, in agreement with previous analyses of mid-infrared spectra. This result is expected from insolation rates: near the equator, the UV solar flux is higher, allowing more C$_2$H$_2$ to be generated from the UV photolysis of CH$_4$. The decrease in abundance towards the poles suggests that horizontal mixing rates are not rapid enough to homogenize the latitudinal distribution. △ Less

Submitted 23 July, 2021; originally announced July 2021.

Comments: Accepted in JGR Planets. 19 pages, 3 figures

Searching for the shadows of giants II: the effect of local ionisation on the Lyman-$α$ absorption signatures of protoclusters at redshift $z\sim2.4$

Authors: Joel S. A. Miller, James S. Bolton, Nina Hatch

Abstract: Local variations in the intergalactic medium (IGM) neutral hydrogen fraction will affect the Ly-$α$ absorption signature of protoclusters identified in tomographic surveys. Using the IllustrisTNG simulations, we investigate how the AGN proximity effect and hot, collisionally ionised gas arising from gravitational infall and black hole feedback changes the Ly-$α$ absorption associated with… ▽ More Local variations in the intergalactic medium (IGM) neutral hydrogen fraction will affect the Ly-$α$ absorption signature of protoclusters identified in tomographic surveys. Using the IllustrisTNG simulations, we investigate how the AGN proximity effect and hot, collisionally ionised gas arising from gravitational infall and black hole feedback changes the Ly-$α$ absorption associated with $M_{z=0}\simeq10^{14}\,M_\odot$ protoclusters at $z\simeq2.4$. We find that protocluster galaxy overdensities exhibit a weak anti-correlation with Ly-$α$ transmission in IGM transmission maps, but local HI ionisation enhancements due to hot $T>10^{6}\rm\,K$ gas or nearby AGN can disrupt this relationship within individual protoclusters. On average, however, we find that strong reductions in the IGM neutral fraction are limited to within $\lesssim 5h^{-1}\,\textrm{cMpc}$ of the dark matter haloes. Local ionisation enhancements will therefore have a minimal impact on the completeness of protocluster identification in tomographic surveys if smoothing Ly-$α$ transmission maps over scales of $\sim4 h^{-1}\,\textrm{cMpc}$, as is typically done in observations. However, if calibrating the relationship between the matter density and Ly-$α$ transmission in tomographic maps using simple analytical models for the Ly-$α$ forest opacity, the presence of hot gas around haloes can still result in systematically lower estimates of $M_{z=0}$ for the most massive protoclusters. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 14 pages (including 1 appendix), 5 figures (plus 2 in the appendix), 1 table. Accepted for publication in MNRAS

The detectability of strong 21 centimetre forest absorbers from the diffuse intergalactic medium in late reionisation models

Authors: Tomáš Šoltinský, James S. Bolton, Nina Hatch, Martin G. Haehnelt, Laura C. Keating, Girish Kulkarni, Ewald Puchwein, Jonathan Chardin, Dominique Aubert

Abstract: A late end to reionisation at redshift $z\simeq 5.3$ is consistent with observed spatial variations in the Ly$α$ forest transmission and the deficit of Ly$α$ emitting galaxies around extended Ly$α$ absorption troughs at $z=5.5$. In this model, large islands of neutral hydrogen should persist in the diffuse intergalactic medium (IGM) until $z\simeq 6$. We use a novel, hybrid approach that combines… ▽ More A late end to reionisation at redshift $z\simeq 5.3$ is consistent with observed spatial variations in the Ly$α$ forest transmission and the deficit of Ly$α$ emitting galaxies around extended Ly$α$ absorption troughs at $z=5.5$. In this model, large islands of neutral hydrogen should persist in the diffuse intergalactic medium (IGM) until $z\simeq 6$. We use a novel, hybrid approach that combines high resolution cosmological hydrodynamical simulations with radiative transfer to predict the incidence of strong 21 cm forest absorbers with optical depths $τ_{21}>10^{-2}$ from the diffuse IGM in these late reionisation models. We include the effect of redshift space distortions on the simulated 21 cm forest spectra, and treat the highly uncertain heating of the pre-reionisation IGM by soft X-rays as a free parameter. For a model with only modest IGM pre-heating, such that average gas kinetic temperatures in the diffuse IGM remain below $T_{\rm K}\simeq 10^{2} \rm\, K$, we find that strong 21 cm forest absorption lines should persist until $z=6$. For a sample of $\sim 10$ sufficiently radio loud background sources, a null-detection of 21 cm forest absorbers at $z\simeq 6$ with SKA1-low or possibly LOFAR should provide an informative lower limit on the still largely unconstrained soft X-ray background at high redshift and the temperature of the pre-reionisation IGM. △ Less

Submitted 24 June, 2021; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: 18 pages, 11 figures, 4 tables. Accepted for publication in MNRAS

The mean free path of ionizing photons at 5 < z < 6: evidence for rapid evolution near reionization

Authors: George D. Becker, Anson D'Aloisio, Holly M. Christenson, Yongda Zhu, Gábor Worseck, James S. Bolton

Abstract: The mean free path of ionizing photons, $λ_{\rm mfp}$, is a key factor in the photoionization of the intergalactic medium (IGM). At $z \gtrsim 5$, however, $λ_{\rm mfp}$ may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of $λ_{\rm mfp}$ that address this bias and extend up to $z \sim 6$ for the first time. Our measurements… ▽ More The mean free path of ionizing photons, $λ_{\rm mfp}$, is a key factor in the photoionization of the intergalactic medium (IGM). At $z \gtrsim 5$, however, $λ_{\rm mfp}$ may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of $λ_{\rm mfp}$ that address this bias and extend up to $z \sim 6$ for the first time. Our measurements at $z \sim 5$ are based on data from the Giant Gemini GMOS survey and new Keck LRIS observations of low-luminosity QSOs. At $z \sim 6$ we use QSO spectra from Keck ESI and VLT X-Shooter. We measure $λ_{\rm mfp} = 9.09^{+1.62}_{-1.28}$ proper Mpc and $0.75^{+0.65}_{-0.45}$ proper Mpc (68% confidence) at $z = 5.1$ and 6.0, respectively. The results at $z = 5.1$ are consistent with existing measurements, suggesting that bias from the proximity effect is minor at this redshift. At $z = 6.0$, however, we find that neglecting the proximity effect biases the result high by a factor of two or more. Our measurement at $z = 6.0$ falls well below extrapolations from lower redshifts, indicating rapid evolution in $λ_{\rm mfp}$ over $5 < z < 6$. This evolution disfavors models in which reionization ended early enough that the IGM had time to fully relax hydrodynamically by $z = 6$, but is qualitatively consistent with models wherein reionization completed at $z = 6$ or even significantly later. Our mean free path results are most consistent with late reionization models wherein the IGM is still 20% neutral at $z=6$, although our measurement at $z = 6.0$ is even lower than these models prefer. △ Less

Submitted 15 September, 2021; v1 submitted 30 March, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS