Showing 1–50 of 135 results for author: Girish

Geophysical Observations of the 24 September 2023 OSIRIS-REx Sample Return Capsule Re-Entry

Authors: Elizabeth A. Silber, Daniel C. Bowman, Chris G. Carr, David P. Eisenberg, Brian R. Elbing, Benjamin Fernando, Milton A. Garcés, Robert Haaser, Siddharth Krishnamoorthy, Charles A. Langston, Yasuhiro Nishikawa, Jeremy Webster, Jacob F. Anderson, Stephen Arrowsmith, Sonia Bazargan, Luke Beardslee, Brant Beck, Jordan W. Bishop, Philip Blom, Grant Bracht, David L. Chichester, Anthony Christe, Kenneth Cummins, James Cutts, Lisa Danielson , et al. (57 additional authors not shown)

Abstract: Sample Return Capsules (SRCs) entering Earth's atmosphere at hypervelocity from interplanetary space are a valuable resource for studying meteor phenomena. The 24 September 2023 arrival of the OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer) SRC provided an unprecedented chance for geophysical observations of a well-characterized source with kn… ▽ More Sample Return Capsules (SRCs) entering Earth's atmosphere at hypervelocity from interplanetary space are a valuable resource for studying meteor phenomena. The 24 September 2023 arrival of the OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer) SRC provided an unprecedented chance for geophysical observations of a well-characterized source with known parameters, including timing and trajectory. A collaborative effort involving researchers from 16 institutions executed a carefully planned geophysical observational campaign at strategically chosen locations, deploying over 400 ground-based sensors encompassing infrasound, seismic, distributed acoustic sensing (DAS), and GPS technologies. Additionally, balloons equipped with infrasound sensors were launched to capture signals at higher altitudes. This campaign (the largest of its kind so far) yielded a wealth of invaluable data anticipated to fuel scientific inquiry for years to come. The success of the observational campaign is evidenced by the near-universal detection of signals across instruments, both proximal and distal. This paper presents a comprehensive overview of the collective scientific effort, field deployment, and preliminary findings. The early findings have the potential to inform future space missions and terrestrial campaigns, contributing to our understanding of meteoroid interactions with planetary atmospheres. Furthermore, the dataset collected during this campaign will improve entry and propagation models as well as augment the study of atmospheric dynamics and shock phenomena generated by meteoroids and similar sources. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 87 pages, 14 figures

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

Impact of extragalactic point sources on the low-frequency sky spectrum and cosmic dawn global 21-cm measurements

Authors: Shikhar Mittal, Girish Kulkarni, Dominic Anstey, Eloy de Lera Acedo

Abstract: Contribution of resolved and unresolved extragalactic point sources to the low-frequency sky spectrum is a potentially non-negligible part of the astrophysical foregrounds for cosmic dawn 21-cm experiments. The clustering of such point sources on the sky, combined with the frequency-dependence of the antenna beam, can also make this contribution chromatic. By combining low-frequency measurements o… ▽ More Contribution of resolved and unresolved extragalactic point sources to the low-frequency sky spectrum is a potentially non-negligible part of the astrophysical foregrounds for cosmic dawn 21-cm experiments. The clustering of such point sources on the sky, combined with the frequency-dependence of the antenna beam, can also make this contribution chromatic. By combining low-frequency measurements of the luminosity function and the angular correlation function of extragalactic point sources, we develop a model for the contribution of these sources to the low-frequency sky spectrum. Using this model, we find that the contribution of sources with flux density $>10^{-6}\,$Jy to the sky-averaged spectrum is smooth and of the order of a few kelvins at 50 - $200\,$MHz. We combine this model with measurements of the galactic foreground spectrum and convolve the result with the beam of the conical log-spiral antenna planned as part of the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) project. We find that the contribution of point sources to resultant spectrum is $\sim0.4\%$ of the total foregrounds, but still larger by at least an order of magnitude than the standard predictions for the cosmological 21-cm signal. As a result, not accounting for the point-source contribution leads to a systematic bias in 21-cm signal recovery. We show, however, that in the REACH case, this reconstruction bias can be removed by modelling the point-source contribution as a power law with a running spectral index. We make our code publicly available as a Python package labelled epspy. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 13 pages, 14 figures, and 1 appendix. Submitted to MNRAS. Comments are welcome

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

Evidence of Pop~III stars' chemical signature in neutral gas at z~6. A study based on the E-XQR-30 spectroscopic sample

Authors: Alessio Sodini, Valentina D'Odorico, Stefania Salvadori, Irene Vanni, Manuela Bischetti, Guido Cupani, Rebecca Davies, George D. Becker, Eduardo Bañados, Sarah Bosman, Frederick Davies, Emanuele Paolo Farina, Andrea Ferrara, Laura Keating, Girish Kulkarni, Samuel Lai, Emma Ryan-Weber, Alma Maria Sebastian, Fabian Walter

Abstract: This study explores the metal enrichment signatures attributed to the first generation of stars (PopIII) in the Universe, focusing on the E-XQR-30 sample. We aim to identify traces of Pop III metal enrichment by analyzing neutral gas in the interstellar medium of primordial galaxies and their satellite clumps, detected in absorption. To chase the chemical signature of PopIII stars, we studied meta… ▽ More This study explores the metal enrichment signatures attributed to the first generation of stars (PopIII) in the Universe, focusing on the E-XQR-30 sample. We aim to identify traces of Pop III metal enrichment by analyzing neutral gas in the interstellar medium of primordial galaxies and their satellite clumps, detected in absorption. To chase the chemical signature of PopIII stars, we studied metal absorption systems in the E-XQR-30 sample, selected through the detection of the OI absorption line at 1302A. The OI line is a reliable tracer of HI and allowed us to overcome the challenges posed by the Lyman-$α$ forest's increasing saturation at redshifts above $\sim5$ to identify Damped Lyman-$α$ systems (DLA). We detected and analyzed 29 OI systems at $z\geq5.4$, differentiating between proximate DLAs (PDLA) and intervening DLAs. Voigt function fits were applied to obtain ionic column densities, and relative chemical abundances were determined for 28 systems. These were then compared with the predictions of theoretical models. Our findings expand the study of OI systems at $z\geq5.4$ fourfold. No systematic differences were observed in the average chemical abundances between PDLAs and intervening DLAs. The chemical abundances in our sample align with literature systems at $z>4.5$, suggesting a similar enrichment pattern for this class of absorption systems. A comparison between these DLA-analogues at $4.5<z<6.5$ with a sample of very metal-poor DLAs at $2<z<4.5$ shows in general similar average values for the relative abundances, with the exception of [C/O], [Si/Fe] and [Si/O] which are significantly larger for the high-$z$ sample. Furthermore, the dispersion of the measurements significantly increases in the high-redshift bin. This increase is predicted by the theoretical models and indicates a potential retention of PopIII signatures in the probed gas. (Abridged) △ Less

Submitted 4 June, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: Accepted by A&A

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

E-XQR-30: The evolution of MgII, CII and OI across 2<z<6

Authors: Alma Maria Sebastian, Emma Ryan-Weber, Rebecca L. Davies, George D. Becker, Laura C. Keating, Valentina D'Odorico, Romain A. Meyer, Sarah E. I. Bosman, Guido Cupani, Girish Kulkarni, Martin G. Haehnelt, Samuel Lai, Anna-Christina Eilers, Manuela Bischetti, Simona Gallerani

Abstract: Intervening metal absorbers in quasar spectra at $z > 6$ can be used as probes to study the chemical enrichment of the Universe during the Epoch of Reionization (EoR). This work presents the comoving line densities ($dn/dX$) of low ionisation absorbers, namely, Mg II (2796Å), C II (1334Å) and O I (1302Å) across $2 <z < 6$ using the E-XQR-30 metal absorber catalog prepared from 42 XSHOOTER quasar s… ▽ More Intervening metal absorbers in quasar spectra at $z > 6$ can be used as probes to study the chemical enrichment of the Universe during the Epoch of Reionization (EoR). This work presents the comoving line densities ($dn/dX$) of low ionisation absorbers, namely, Mg II (2796Å), C II (1334Å) and O I (1302Å) across $2 <z < 6$ using the E-XQR-30 metal absorber catalog prepared from 42 XSHOOTER quasar spectra at $5.8 < z < 6.6$. Here, we analyse 280 Mg II ($1.9 < z < 6.4$), 22 C II ($5.2 < z < 6.4$) and 10 O I ($5.3 < z < 6.4$) intervening absorbers, thereby building up on previous studies with improved sensitivity of 50% completeness at an equivalent width of $W > 0.03$Å. For the first time, we present the comoving line densities of 131 weak ($W < 0.3$Å) intervening Mg II absorbers at $1.9 < z < 6.4$ which exhibit constant evolution with redshift similar to medium ($0.3 < W < 1.0$Å) absorbers. However, the cosmic mass density of Mg II - dominated by strong Mg II systems - traces the evolution of global star formation history from redshift 1.9 to 5.5. E-XQR-30 also increases the absorption path length by a factor of 50% for C II and O I whose line densities show a rising trend towards $z > 5$, in agreement with previous works. In the context of a decline in metal enrichment of the Universe at $z > 5$, the overall evolution in the incidence rates of absorption systems can be explained by a weak - possibly soft fluctuating - UV background. Our results, thereby, provide evidence for a late reionization continuing to occur in metal-enriched and therefore, biased regions in the Universe. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: 20 pages, 21 figures

The dynamics of self-gravity wakes in the Mimas 5:3 bending wave: modifying the linear theory

Authors: Daniel D. Sega, Glen. Stewart, Josh E. Colwell, Girish M. Duvvuri, Richard Jerousek, Larry Esposito

Abstract: The satellite Mimas launches a bending wave -- a warping of the rings that propagates radially through self-gravity -- at the 5:3 inner vertical resonance with Saturn's rings. We present a modification of the linear bending wave theory which includes the effects of satellite self-gravity wakes on the particles in the wave. We show that, when treated as rigid, these wakes generate an extra layer of… ▽ More The satellite Mimas launches a bending wave -- a warping of the rings that propagates radially through self-gravity -- at the 5:3 inner vertical resonance with Saturn's rings. We present a modification of the linear bending wave theory which includes the effects of satellite self-gravity wakes on the particles in the wave. We show that, when treated as rigid, these wakes generate an extra layer of particles whose number density is proportional to the magnitude of the slope of the warped ring. Using a ray-tracing code we compare our predictions with those of linear bending wave theory and with 60 stellar occultations observed by the Cassini Ultraviolet Imaging Spectrograph (UVIS) and find that the extra layer of particles of our perturbed bending wave model has a considerable explanatory power for the UVIS dataset. Our best model explains the most discrepant and surprising features of the Mimas 5:3 bending wave; the enhancement of the signal for the cases of occultations with high ring opening angle and the bigger-than-expected viscosity, $ν= 576 \, \mathrm{cm^2/s}$, which is more than double the viscosity computed from density waves. This shows that self-gravity wakes can be effective at transporting angular momentum in a vertically perturbed disk. Relative to neighboring density waves, we find a lower-than-expected value for the surface mass density, $σ= 36.7 \, \mathrm{g/cm^2}$, which suggests that the enhanced viscous interactions may be transporting material into the surrounding regions. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 33 pages, 14 figures

Journal ref: Icarus, Volume 413, 2024, ISSN 0019-1035,

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): Impact of Galaxies on the CGM Metal Enrichment at z > 6 Using the JWST and VLT

Authors: Siwei Zou, Zheng Cai, Feige Wang, Xiaohui Fan, Jaclyn B. Champagne, Joseph F. Hennawi, Jan-Torge Schindler, Emanuele P. Farina, Jinyi Yang, Kohei Inayoshi, Eduardo Banados, Sarah E. I. Bosman, Zihao Li, Xiaojing Lin, Yunjing Wu, Fengwu Sun, Zi-Yi Guo, Girish Kulkarni, Melanie Habouzit, Stephane Charlot, Jacopo Chevallard, Thomas Connor, Anna-Christina Eilers, Linhua Jiang, Xiangyu Jin , et al. (5 additional authors not shown)

Abstract: We characterize the multiphase circumgalactic medium and galaxy properties at z = 6.0-6.5 in four quasar fields from the James Webb Space Telescope A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE) program. We use the Very Large Telescope/X-shooter spectra of quasar J0305-3150 to identify one new metal absorber at z = 6.2713 with multiple transitions (OI, MgI, FeII and CII).… ▽ More We characterize the multiphase circumgalactic medium and galaxy properties at z = 6.0-6.5 in four quasar fields from the James Webb Space Telescope A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE) program. We use the Very Large Telescope/X-shooter spectra of quasar J0305-3150 to identify one new metal absorber at z = 6.2713 with multiple transitions (OI, MgI, FeII and CII). They are combined with the published absorbing systems in Davies et al. (2023a) at the same redshift range to form of a sample of nine metal absorbers at z = 6.03 to 6.49. We identify eight galaxies within 1000 km s$^{-1}$ and 350 kpc around the absorbing gas from the ASPIRE spectroscopic data, with their redshifts secured by [OIII]($λλ$4959, 5007) doublets and H$β$ emission lines. Our spectral energy distribution fitting indicates that the absorbing galaxies have stellar mass ranging from 10$^{7.2}$ to 10$^{8.8}M_{\odot}$ and metallicity between 0.02 and 0.4 solar. Notably, the z = 6.2713 system in the J0305-3150 field resides in a galaxy overdensity region, which contains two (tentatively) merging galaxies within 350 kpc and seven galaxies within 1 Mpc. We measure the relative abundances of $α$ elements to iron ([$α$/Fe]) and find that the CGM gas in the most overdense region exhibits a lower [$α$/Fe] ratio. Our modeling of the galaxy's chemical abundance favors a top-heavy stellar initial mass function, and hints that we may be witnessing the contribution of the first generation Population III stars to the CGM at the end of reionization epoch. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: 21 pages, 4 figures in the main text. Accepted for publication in ApJL

Chronicling the reionization history at $6\lesssim z \lesssim 7$ with emergent quasar damping wings

Authors: Dominika Ďurovčíková, Anna-Christina Eilers, Huanqing Chen, Sindhu Satyavolu, Girish Kulkarni, Robert A. Simcoe, Laura C. Keating, Martin G. Haehnelt, Eduardo Bañados

Abstract: The spectra of high-redshift ($z\gtrsim 6$) quasars contain valuable information on the progression of the Epoch of Reionization (EoR). At redshifts $z<6$, the observed Lyman-series forest shows that the intergalactic medium (IGM) is nearly ionized, while at $z>7$ the observed quasar damping wings indicate high neutral gas fractions. However, there remains a gap in neutral gas fraction constraints… ▽ More The spectra of high-redshift ($z\gtrsim 6$) quasars contain valuable information on the progression of the Epoch of Reionization (EoR). At redshifts $z<6$, the observed Lyman-series forest shows that the intergalactic medium (IGM) is nearly ionized, while at $z>7$ the observed quasar damping wings indicate high neutral gas fractions. However, there remains a gap in neutral gas fraction constraints at $6\lesssim z \lesssim 7$ where the Lyman series forest becomes saturated but damping wings have yet to fully emerge. In this work, we use a sample of 18 quasar spectra at redshifts $6.0<z<7.1$ to close this gap. We apply neural networks to reconstruct the quasars' continuum emission around the partially absorbed Lyman $α$ line to normalize their spectra, and stack these continuum-normalized spectra in three redshift bins. To increase the robustness of our results, we compare the stacks to a grid of models from two hydrodynamical simulations, ATON and CROC, and we measure the volume-averaged neutral gas fraction, $\bar{x}_{\rm HI}$, while jointly fitting for the mean quasar lifetime, $t_{\rm Q}$, for each stacked spectrum. We chronicle the evolution of neutral gas fraction using the ATON (CROC) models as follows: $\bar{x}_{\rm HI} = 0.21_{-0.07}^{+0.17}$ ($\bar{x}_{\rm HI} = 0.10_{<10^{-4}}^{+0.73}$) at $\langle z \rangle =6.10$, $\bar{x}_{\rm HI} = 0.21_{-0.07}^{+0.33}$ ($\bar{x}_{\rm HI} =0.57_{-0.47}^{+0.26}$) at $\langle z \rangle =6.46$, and $\bar{x}_{\rm HI} = 0.37_{-0.17}^{+0.17}$ ($\bar{x}_{\rm HI} =0.57_{-0.21}^{+0.26}$) at $\langle z \rangle =6.87$. At the same time we constrain the average quasar lifetime to be $t_{\rm Q} \lesssim 7\ {\rm Myr}$ across all redshift bins, in good agreement with previous studies. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: 28 pages, 14 figures

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

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

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

Submitted 13 December, 2023; originally announced December 2023.

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

The life cycle of stars and their planets from the high energy perspective

Authors: Lia Corrales, Keivan G. Stassun, Tim Cunningham, Girish Duvvuri, Jeremy J. Drake, Catherine Espaillat, Adina D. Feinstein, Elena Gallo, Hans Moritz Gunther, George W. King, Marina Kounkel, Carey M. Lisse, Rodolfo Montez Jr., David A. Principe, Jesus A. Toala, Scott J. Wolk, Raven Cilley, Tansu Daylan, Margarita Karovska, Pragati Pradhan, Peter J. Wheatley, Jun Yang

Abstract: One of the key research themes identified by the Astro2020 decadal survey is Worlds and Suns in Context. The Advanced X-ray Imaging Satellite (AXIS) is a proposed NASA APEX mission that will become the prime high-energy instrument for studying star-planet connections from birth to death. This work explores the major advances in this broad domain of research that will be enabled by the AXIS mission… ▽ More One of the key research themes identified by the Astro2020 decadal survey is Worlds and Suns in Context. The Advanced X-ray Imaging Satellite (AXIS) is a proposed NASA APEX mission that will become the prime high-energy instrument for studying star-planet connections from birth to death. This work explores the major advances in this broad domain of research that will be enabled by the AXIS mission, through X-ray observations of stars in clusters spanning a broad range of ages, flaring M-dwarf stars known to host exoplanets, and young stars exhibiting accretion interactions with their protoplanetary disks. In addition, we explore the ability of AXIS to use planetary nebulae, white dwarfs, and the Solar System to constrain important physical processes from the microscopic (e.g., charge exchange) to the macroscopic (e.g., stellar wind interactions with the surrounding interstellar medium). △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: This White Paper is part of a series commissioned for the AXIS Probe Concept Mission

Robustness of direct measurements of the mean free path of ionizing photons in the epoch of reionization

Authors: Sindhu Satyavolu, Girish Kulkarni, Laura C. Keating, Martin G. Haehnelt

Abstract: Measurements of the mean free path of Lyman-continuum photons in the intergalactic medium during the epoch of reionization can help constrain the nature of the sources as well as sinks of hydrogen-ionizing radiation. A recent approach to this measurement has been to utilize composite spectra of multiple quasars at $z\sim 6$, and infer the mean free path after correcting the spectra for the presenc… ▽ More Measurements of the mean free path of Lyman-continuum photons in the intergalactic medium during the epoch of reionization can help constrain the nature of the sources as well as sinks of hydrogen-ionizing radiation. A recent approach to this measurement has been to utilize composite spectra of multiple quasars at $z\sim 6$, and infer the mean free path after correcting the spectra for the presence of quasar proximity zones. This has revealed not only a steep drop in the mean free path from $z=5$ to $z=6$, but also potentially a mild tension with reionization simulations. We critically examine such direct measurements of the mean free path for biases due to quasar environment, incomplete reionization, and quasar proximity zones. Using cosmological radiative transfer simulations of reionization combined with one-dimensional radiative transfer calculations of quasar proximity zones, we find that the bias in the mean free path due to overdensities around quasars is minimal at $z\sim 6$. Patchiness of reionization at this redshift also does not affect the measurements significantly. Fitting our model to the data results in a mean free path of $λ_{\mathrm{mfp}}=1.49^{+0.47}_{-0.52}$~pMpc at $z=6$, which is consistent with the recent measurements in the literature, indicating robustness with respect to the modelling of quasar proximity zones. We also compare various ways in which the mean free path has been defined in simulations before the end of reionization. Overall, our finding is that recent measurements of the mean free path appear to be robust relative to several sources of potential bias. △ Less

Submitted 12 July, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Accepted in MNRAS. Mean free path measurement revised to a higher value with narrower uncertainty using full data covariance

Radiative transfer of Lyman-$α$ photons at cosmic dawn with realistic gas physics

Authors: Shikhar Mittal, Girish Kulkarni, Thibault Garel

Abstract: The cosmic dawn 21-cm signal is enabled by Ly~$α$ photons through a process called the Wouthuysen-Field effect. An accurate model of the signal in this epoch hinges on the accuracy of the computation of the Ly~$α$ coupling, which requires one to calculate the specific intensity of UV radiation from sources such as the first stars. Most traditional calculations of the Ly~$α$ coupling assume a delta… ▽ More The cosmic dawn 21-cm signal is enabled by Ly~$α$ photons through a process called the Wouthuysen-Field effect. An accurate model of the signal in this epoch hinges on the accuracy of the computation of the Ly~$α$ coupling, which requires one to calculate the specific intensity of UV radiation from sources such as the first stars. Most traditional calculations of the Ly~$α$ coupling assume a delta-function scattering cross-section, as the resonant nature of the Ly~$α$ scattering makes an accurate radiative transfer solution computationally expensive. Attempts to improve upon this traditional approach using numerical radiative transfer have recently emerged. However, the radiative transfer computation in these treatments suffers from assumptions such as a uniform density of intergalactic gas, zero gas temperature, and absence of gas bulk motion, or numerical approximations such as core skipping. We investigate the role played by these approximations in setting the value of the Ly~$α$ coupling and the 21-cm signal at cosmic dawn. We present results of Monte Carlo radiative transfer simulations, without core skipping, and show that neglecting gas temperature in the radiative transfer significantly underestimates the scattering rate and hence the Ly~$α$ coupling and the 21-cm signal. We also discuss the effect of these processes on the 21-cm power spectrum from the cosmic dawn. This work points the way towards higher-accuracy models to enable better inferences from future measurements. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 20 pages, 14 figures, 2 appendices. Submitted to MNRAS. Comments are welcome

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

A Multi-wavelength Study of Swift J0503.7-2819: a Chimeric magnetic CV

Authors: Kala G. Pradeep, Kulinder Pal Singh, G. C. Dewangan, Elias Aydi, P. E. Barrett, D. A. H. Buckley, V. Girish, K. L. Page, S. B. Potter, E. M. Schlegel

Abstract: We present multi-wavelength temporal and spectral characteristics of a magnetic cataclysmic variable (MCV) Swift J0503.7-2819, using far ultraviolet (FUV) and X-ray data from AstroSat, supplemented with optical data from the Southern African Large Telescope and X-ray data from the XMM-Newton and Swift observatories. The X-ray modulations at 4897.6657 s and 3932.0355 s are interpreted as the orbita… ▽ More We present multi-wavelength temporal and spectral characteristics of a magnetic cataclysmic variable (MCV) Swift J0503.7-2819, using far ultraviolet (FUV) and X-ray data from AstroSat, supplemented with optical data from the Southern African Large Telescope and X-ray data from the XMM-Newton and Swift observatories. The X-ray modulations at 4897.6657 s and 3932.0355 s are interpreted as the orbital ($P_Ω$) and spin ($P_ω$) period, respectively, and are consistent with prior reports. With a spin-orbit period ratio of 0.8 and $P_Ω$ falling below the period gap (2-3 hrs) of CVs, Swift J0503.7-2819 would be the newest addition to the growing population of nearly synchronous MCVs, which we call EX Hya-like systems. Hard X-ray luminosity of $<$ $2.5\times10^{32} erg s^{-1}$, as measured with the Swift Burst Alert Telescope, identifies it to be a low-luminosity intermediate polar, similar to other EX Hya-like systems. The phenomenology of the light curves and the spectral characteristics rule out a purely disc-fed/stream-fed model and instead reveal the presence of complex accretion structures around the white dwarf. We propose a ring-like accretion flow, akin to EX Hya, using period ratio, stability arguments, and observational features. An attempt is made to differentiate between the asynchronous polar/nearly-synchronous intermediate polar nature of Swift J0503.7-2819. Further, we note that with the advent of sensitive surveys, a growing population of MCVs that exhibit characteristics of both polars and intermediate polars is beginning to be identified, likely forming a genealogical link between the two conventional classes of MCVs. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS, 17 Pages, 13 Figures

The High-Energy Spectrum of the Young Planet Host V1298 Tau

Authors: Girish M. Duvvuri, P. Wilson Cauley, Fernando Cruz Aguirre, Roy Kilgard, Kevin France, Zachory K. Berta-Thompson, J. Sebastian Pineda

Abstract: V1298 Tau is a young pre-main sequence star hosting four known exoplanets that are prime targets for transmission spectroscopy with current-generation instruments. This work pieces together observations from the NICER X-ray telescope, the Space Telescope Imaging Spectrograph and Cosmic Origins Spectrograph instruments aboard Hubble Space Telescope, and empirically informed models to create a panch… ▽ More V1298 Tau is a young pre-main sequence star hosting four known exoplanets that are prime targets for transmission spectroscopy with current-generation instruments. This work pieces together observations from the NICER X-ray telescope, the Space Telescope Imaging Spectrograph and Cosmic Origins Spectrograph instruments aboard Hubble Space Telescope, and empirically informed models to create a panchromatic spectral energy distribution for V1298 Tau spanning 1 to 100000 Angstroms. We describe the methods and assumptions used to assemble the panchromatic spectrum and show that despite this star's brightness, its high-energy spectrum is near the limit of present X-ray and ultraviolet observatories' abilities to characterize. We conclude by using the V1298 Tau spectrum as a benchmark for the activity saturation stage of high-energy radiation from solar-mass stars to compare the lifetime cumulative high-energy irradiation of the V1298 Tau planets to other planets orbiting similarly massive stars. △ Less

Submitted 29 September, 2023; originally announced October 2023.

Comments: The V1298 Tau panchromatic SED will be available via the Data Behind the Figure feature for Figure 1 of the final journal article, will be uploaded to the Exoplanet Archive as an IPAC table, and can be provided upon request

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)

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

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

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

Submitted 6 September, 2023; originally announced September 2023.

Comments: 17 pages, 13 figures, submitted to MNRAS

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

The 21-cm bispectrum from neutral hydrogen islands at z < 6

Authors: Janakee Raste, Girish Kulkarni, Catherine A. Watkinson, Laura C. Keating, Martin G. Haehnelt

Abstract: Spatial variations in the Lyman-$α$ forest opacity at $z<6$ seem to require a late end to cosmic reionization. In this picture, the universe contains neutral hydrogen 'islands' of up to 100 cMpc$/h$ in extent down to redshifts as low as $z\sim 5.3$. This delayed end to reionization also seems to be corroborated by various other observables. An implication of this scenario is that the power spectru… ▽ More Spatial variations in the Lyman-$α$ forest opacity at $z<6$ seem to require a late end to cosmic reionization. In this picture, the universe contains neutral hydrogen 'islands' of up to 100 cMpc$/h$ in extent down to redshifts as low as $z\sim 5.3$. This delayed end to reionization also seems to be corroborated by various other observables. An implication of this scenario is that the power spectrum of the cosmological 21-cm signal at $z<6$ is enhanced relative to conventional reionization models by orders of magnitude. However, these neutral hydrogen islands are also predicted to be at the locations of the deepest voids in the cosmological large-scale structure. As a result, the distribution of the 21-cm signal from them is highly non-Gaussian. We derive the 21-cm bispectrum signal from these regions using high-dynamic-range radiative transfer simulations of reionization. We find that relative to conventional models in which reionization is complete at $z>6$, our model has a significantly larger value of the 21-cm bispectrum. The neutral islands also imprint a feature in the isosceles bispectrum at a characteristic scale of $\sim 1$ cMpc$^{-1}$. We also study the 21-cm bispectrum for general triangle configuration by defining a triangle index. It should be possible to detect the 21-cm bispectrum signal at $ν\gtrsim 200$ MHz using SKA1-LOW for 1080 hours of observation, assuming optimistic foreground removal. △ Less

Submitted 5 March, 2024; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: 12 pages, 12 figures

Journal ref: MNRAS 529 (2024) 129-140

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

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

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

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

Submitted 8 August, 2023; originally announced August 2023.

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

The Variable Detection of Atmospheric Escape around the young, Hot Neptune AU Mic b

Authors: Keighley E. Rockcliffe, Elisabeth R. Newton, Allison Youngblood, Girish M. Duvvuri, Peter Plavchan, Peter Gao, Andrew W. Mann, Patrick J. Lowrance

Abstract: Photoevaporation is a potential explanation for several features within exoplanet demographics. Atmospheric escape observed in young Neptune-sized exoplanets can provide insight into and characterize which mechanisms drive this evolution and at what times they dominate. AU Mic b is one such exoplanet, slightly larger than Neptune (4.19 Earth radii). It closely orbits a 23 Myr pre-Main Sequence M d… ▽ More Photoevaporation is a potential explanation for several features within exoplanet demographics. Atmospheric escape observed in young Neptune-sized exoplanets can provide insight into and characterize which mechanisms drive this evolution and at what times they dominate. AU Mic b is one such exoplanet, slightly larger than Neptune (4.19 Earth radii). It closely orbits a 23 Myr pre-Main Sequence M dwarf with a period of 8.46 days. We obtained two visits of AU Mic b at Lyman-alpha with HST/STIS. One flare within the first HST visit is characterized and removed from our search for a planetary transit. We present a non-detection in our first visit followed by the detection of escaping neutral hydrogen ahead of the planet in our second visit. The outflow absorbed about 30% of the star's Lyman-alpha blue-wing 2.5 hours before the planet's white-light transit. We estimate the highest velocity escaping material has a column density of 10^13.96 cm^-2 and is moving 61.26 km/s away from the host star. AU Mic b's large high energy irradiation could photoionize its escaping neutral hydrogen in 44 minutes, rendering it temporarily unobservable. Our time-variable Lyman-alpha transit ahead of AU Mic b could also be explained by an intermediate stellar wind strength from AU Mic that shapes the escaping material into a leading tail. Future Lyman-alpha observations of this system will confirm and characterize the unique variable nature of its Lyman-alpha transit, which combined with modeling will tune the importance of stellar wind and photoionization. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 24 pages, 11 figures

Journal ref: The Astronomical Journal, Volume 166, Number 2, 2023

The MUSCLES Extension for Atmospheric and Transmission Spectroscopy: UV and X-ray Host-star Observations for JWST ERS & GTO Targets

Authors: Patrick R. Behr, Kevin France, Alexander Brown, Girish Duvvuri, Jacob L. Bean, Zachory Berta-Thompson, Cynthia Froning, Yamila Miguel, J. Sebastian Pineda, David Wilson, Allison Youngblood

Abstract: X-ray through infrared spectral energy distributions (SEDs) are essential for understanding a star's effect on exoplanet atmospheric composition and evolution. We present a catalog of panchromatic SEDs, hosted on the Barbara A. Mikulski Archive for Space Telescopes (MAST), for 11 exoplanet hosting stars which have guaranteed JWST observation time as part of the ERS or GTO programs but have no prev… ▽ More X-ray through infrared spectral energy distributions (SEDs) are essential for understanding a star's effect on exoplanet atmospheric composition and evolution. We present a catalog of panchromatic SEDs, hosted on the Barbara A. Mikulski Archive for Space Telescopes (MAST), for 11 exoplanet hosting stars which have guaranteed JWST observation time as part of the ERS or GTO programs but have no previous UV characterization. The stars in this survey range from spectral type F4-M6 (0.14-1.57 M$_\odot$), rotation periods of ~4-132 days, and ages of approximately 0.5-11.4 Gyr. The SEDs are composite spectra using data from the Chandra X-ray Observatory and XMM-Newton, the Hubble Space Telescope, BT-Settl stellar atmosphere models, and scaled spectra of proxy stars of similar spectral type and activity. From our observations, we have measured a set of UV and X-ray fluxes as indicators of stellar activity level. We compare the chromospheric and coronal activity indicators of our exoplanet-hosting stars to the broader population of field stars and find that a majority of our targets have activity levels lower than the average population of cool stars in the solar neighborhood. This suggests that using SEDs of stars selected from exoplanet surveys to compute generic exoplanet atmosphere models may underestimate the typical host star's UV flux by an order of magnitude or more, and consequently, that the observed population of exoplanetary atmospheres receive lower high-energy flux levels than the typical planet in the solar neighborhood. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: 28 pages, 10 figures, 7 tables, accepted for publication in AJ

New quasar proximity zone size measurements at $z\sim 6$ using the enlarged XQR-30 sample

Authors: Sindhu Satyavolu, Anna-Christina Eilers, Girish Kulkarni, Emma Ryan-Weber, Rebecca L. Davies, George D. Becker, Sarah E. I. Bosman, Bradley Greig, Chiara Mazzucchelli, Eduardo Bañados, Manuela Bischetti, Valentina D'Odorico, Xiaohui Fan, Emanuele Paolo Farina, Martin G. Haehnelt, Laura C. Keating, Samuel Lai, Fabian Walter

Abstract: Proximity zones of high-redshift quasars are unique probes of their central supermassive black holes as well as the intergalactic medium in the last stages of reionization. We present 22 new measurements of proximity zones of quasars with redshifts between 5.8 and 6.6, using the enlarged XQR-30 sample of high-resolution, high-SNR quasar spectra. The quasars in our sample have UV magnitudes of… ▽ More Proximity zones of high-redshift quasars are unique probes of their central supermassive black holes as well as the intergalactic medium in the last stages of reionization. We present 22 new measurements of proximity zones of quasars with redshifts between 5.8 and 6.6, using the enlarged XQR-30 sample of high-resolution, high-SNR quasar spectra. The quasars in our sample have UV magnitudes of $M_{1450}\sim -27$ and black hole masses of $10^9$$\unicode{x2013}$$10^{10}$ M$_\odot$. Our inferred proximity zone sizes are 2$\unicode{x2013}$7 physical Mpc, with a typical uncertainty of less than 0.5 physical Mpc, which, for the first time, also includes uncertainty in the quasar continuum. We find that the correlation between proximity zone sizes and the quasar redshift, luminosity, or black hole mass, indicates a large diversity of quasar lifetimes. Two of our proximity zone sizes are exceptionally small. The spectrum of one of these quasars, with $z=6.02$, displays, unusually for this redshift, damping wing absorption without any detectable metal lines, which could potentially originate from the IGM. The other quasar has a high-ionization absorber $\sim$0.5 pMpc from the edge of the proximity zone. This work increases the number of proximity zone measurements available in the last stages of cosmic reionization to 87. This data will lead to better constraints on quasar lifetimes and obscuration fractions at high redshift, which in turn will help probe the seed mass and formation redshift of supermassive black holes. △ Less

Submitted 1 May, 2023; originally announced May 2023.

Comments: 16 pages, 9 figures, Accepted in MNRAS

Report number: TIFR/TH/23-5

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): JWST Reveals a Filamentary Structure around a z=6.61 Quasar

Authors: Feige Wang, Jinyi Yang, Joseph F. Hennawi, Xiaohui Fan, Fengwu Sun, Jaclyn B. Champagne, Tiago Costa, Melanie Habouzit, Ryan Endsley, Zihao Li, Xiaojing Lin, Romain A. Meyer, Jan-Torge Schindler, Yunjing Wu, Eduardo Bañados, Aaron J. Barth, Aklant K. Bhowmick, Rebekka Bieri, Laura Blecha, Sarah Bosman, Zheng Cai, Luis Colina, Thomas Connor, Frederick B. Davies, Roberto Decarli , et al. (34 additional authors not shown)

Abstract: We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) gal… ▽ More We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) galaxy redshift survey at 3-4 $μ$m among JWST Cycle-1 programs and provide extensive legacy values for studying the formation of the earliest supermassive black holes (SMBHs), the assembly of galaxies, early metal enrichment, and cosmic reionization. In this first ASPIRE paper, we report the discovery of a filamentary structure traced by the luminous quasar J0305-3150 and ten [OIII] emitters at $z=6.6$. This structure has a 3D galaxy overdensity of $δ_{\rm gal}=12.6$ over 637 cMpc$^3$, one of the most overdense structures known in the early universe, and could eventually evolve into a massive galaxy cluster. Together with existing VLT/MUSE and ALMA observations of this field, our JWST observations reveal that J0305-3150 traces a complex environment where both UV-bright and dusty galaxies are present, and indicate that the early evolution of galaxies around the quasar is not simultaneous. In addition, we discovered 31 [OIII] emitters in this field at other redshifts, $5.3<z<6.7$, with half of them situated at $z\sim5.4$ and $z\sim6.2$. This indicates that star-forming galaxies, such as [OIII] emitters, are generally clustered at high redshifts. These discoveries demonstrate the unparalleled redshift survey capabilities of NIRCam WFSS and the potential of the full ASPIRE survey dataset. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: accepted for publication in ApJL

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): A First Look at the Rest-frame Optical Spectra of $z > 6.5$ Quasars Using JWST

Authors: Jinyi Yang, Feige Wang, Xiaohui Fan, Joseph F. Hennawi, Aaron J. Barth, Eduardo Bañados, Fengwu Sun, Weizhe Liu, Zheng Cai, Linhua Jiang, Zihao Li, Masafusa Onoue, Jan-Torge Schindler, Yue Shen, Yunjing Wu, Aklant K. Bhowmick, Rebekka Bieri, Laura Blecha, Sarah Bosman, Jaclyn B. Champagne, Luis Colina, Thomas Connor, Tiago Costa, Frederick B. Davies, Roberto Decarli , et al. (31 additional authors not shown)

Abstract: Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''… ▽ More Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE)" program. Our JWST spectra cover the quasars' emission between rest frame $\sim$ 4100 and 5100 Å. The profiles of these quasars' broad H$β$ emission lines span a FWHM from 3000 to 6000 $\rm{km~s^{-1}}$. The H$β$-based virial black hole (BH) masses, ranging from 0.6 to 2.1 billion solar masses, are generally consistent with their MgII-based BH masses. The new measurements based on the more reliable H$β$ tracer thus confirm the existence of billion solar-mass BHs in the reionization epoch. In the observed [OIII] $λλ$4960,5008 doublets of these luminous quasars, broad components are more common than narrow core components ($\le~1200~\rm{km~s^{-1}}$), and only one quasar shows stronger narrow components than broad. Two quasars exhibit significantly broad and blueshifted [OIII] emission, thought to trace galactic-scale outflows, with median velocities of $-610~\rm{km~s^{-1}}$ and $-1430~\rm{km~s^{-1}}$ relative to the [CII] $158\,μ$m line. All eight quasars show strong optical FeII emission, and follow the Eigenvector 1 relations defined by low-redshift quasars. The entire ASPIRE program will eventually cover 25 quasars and provide a statistical sample for the studies of the BHs and quasar spectral properties. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 13 pages, 4 figures, accepted for publication in ApJL

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

Coronal X-Ray Emission from Nearby, Low-Mass, Exoplanet Host Stars Observed by the MUSCLES and Mega-MUSCLES HST Treasury Survey Projects

Authors: Alexander Brown, P. Christian Schneider, Kevin France, Cynthia S. Froning, Allison A. Youngblood, David J. Wilson, R. O. Parke Loyd, J. Sebastian Pineda, Girish M. Duvvuri, Adam F. Kowalski, Zachory K. Berta-Thompson

Abstract: The high energy X-ray and ultraviolet (UV) radiation fields of exoplanet host stars play a crucial role in controlling the atmospheric conditions and the potential habitability of exoplanets. Major surveys of the X-ray/UV emissions from late-type (K and M spectral type) exoplanet hosts have been conducted by the MUSCLES and Mega-MUSCLES Hubble Space Telescope (HST) Treasury programs. These samples… ▽ More The high energy X-ray and ultraviolet (UV) radiation fields of exoplanet host stars play a crucial role in controlling the atmospheric conditions and the potential habitability of exoplanets. Major surveys of the X-ray/UV emissions from late-type (K and M spectral type) exoplanet hosts have been conducted by the MUSCLES and Mega-MUSCLES Hubble Space Telescope (HST) Treasury programs. These samples primarily consist of relatively old, ``inactive'', low mass stars. In this paper we present results from X-ray observations of the coronal emission from these stars obtained using the Chandra X-ray Observatory, the XMM-Newton Observatory, and the Neil Gehrels Swift Observatory. The stars effectively sample the coronal activity of low-mass stars at a wide range of masses and ages. The vast majority (21 of 23) of the stars are detected and their X-ray luminosities measured. Short-term flaring variability is detected for most of the fully-convective (M $\leq$ 0.35 M$_{\odot}$) stars but not for the more massive M dwarfs during these observations. Despite this difference, the mean X-ray luminosities for these two sets of M dwarfs are similar with more massive (0.35 M$_{\odot}$ $\leq$ M $\leq$ 0.6 M$_{\odot}$) M dwarfs at $\sim$5 $\times$ 10$^{26}$ erg s$^{-1}$ compared to $\sim$2 $\times$ 10$^{26}$ erg s$^{-1}$ for fully-convective stars older than 1 Gyr. Younger, fully-convective M dwarfs have X-ray luminosities between 3 and 6 $\times$ 10$^{27}$ erg s$^{-1}$.The coronal X-ray spectra have been characterized and provide important information that is vital for the modeling of the stellar EUV spectra. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: 39 pages, 15 figures. Accepted for publication in The Astronomical Journal

Antennas for low-frequency radio telescope of SKA

Authors: Agaram Raghunathan, Keerthipriya Satish, Arasi Sathyamurthy, T. Prabu, B. S. Girish, K. S. Srivani, Shiv K. Sethi

Abstract: The low-frequency radio telescope of the Square Kilometre Array (SKA) is being built by the international radio astronomical community to (i) have orders of magnitude higher sensitivity and (ii) be able to map the sky several hundred times faster, than any other existing facilities over the frequency range of 50 - 350 MHz. The sensitivity of a radio telescope array is in general, dependent upon th… ▽ More The low-frequency radio telescope of the Square Kilometre Array (SKA) is being built by the international radio astronomical community to (i) have orders of magnitude higher sensitivity and (ii) be able to map the sky several hundred times faster, than any other existing facilities over the frequency range of 50 - 350 MHz. The sensitivity of a radio telescope array is in general, dependent upon the number of electromagnetic sensors used to receive the sky signal. The total number of them is further constrained by the effects of mutual coupling between the sensor elements, allowable grating lobes in their radiation patterns, etc. The operating frequency band is governed by the desired spatial and spectral responses, acceptable sidelobe and backlobe levels, radiation efficiency, polarization purity and calibratability of sensors' response. This paper presents a brief review of several broadband antennas considered as potential candidates by various engineering groups across the globe, for the low-frequency radio telescope of SKA covering the frequency range of 50 - 350 MHz, on the basis of their suitability for conducting primary scientific objectives. △ Less

Submitted 26 May, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 14 pages, 33 figures, JoAA - Special issue on the SKA (2023) - Accepted for publication

Progression of Digital-Receiver Architecture: From MWA to SKA1-Low,and beyond

Authors: Girish B. S., Harshavardhan Reddy S., Shiv Sethi, Srivani K. S., Abhishek R., Ajithkumar B., Sahana Bhattramakki, Kaushal Buch, Sandeep Chaudhuri, Yashwant Gupta, Kamini P. A., Sanjay Kudale, Madhavi S., Mekhala Muley, Prabu T., Raghunathan A., Shelton G. J

Abstract: Backed by advances in digital electronics, signal processing, computation, and storage technologies, aperture arrays, which had strongly influenced the design of telescopes in the early years of radio astronomy, have made a comeback. Amid all these developments, an international effort to design and build the world's largest radio telescope, the Square Kilometre Array (SKA), is ongoing. With its v… ▽ More Backed by advances in digital electronics, signal processing, computation, and storage technologies, aperture arrays, which had strongly influenced the design of telescopes in the early years of radio astronomy, have made a comeback. Amid all these developments, an international effort to design and build the world's largest radio telescope, the Square Kilometre Array (SKA), is ongoing. With its vast collecting area of 1 sq-km, the SKA is envisaged to provide unsurpassed sensitivity and leverage technological advances to implement a complex receiver to provide a large field of view through multiple beams on the sky. Many pathfinders and precursor aperture array telescopes for the SKA, operating in the frequency range of 10-300 MHz, have been constructed and operationalized to obtain valuable feedback on scientific, instrumental, and functional aspects. This review article looks explicitly into the progression of digital-receiver architecture from the Murchison Widefield Array (precursor) to the SKA1-Low. It highlights the technological advances in analog-to-digital converters (ADCs),field-programmable gate arrays (FPGAs), and central processing unit-graphics processing unit (CPU-GPU) hybrid platforms around which complex digital signal processing systems implement efficient channelizers, beamformers, and correlators. The article concludes with a preview of the design of a new generation signal processing platform based on radio frequency system-on-chip (RFSoC). △ Less

Submitted 17 January, 2023; originally announced January 2023.

Comments: 18 pages, 4 figures, Accepted for publication in the special issue (2023) on the SKA from the JoAA

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

Halo mass function in scale invariant models

Authors: Swati Gavas, J S Bagla, Nishikanta Khandai, Girish Kulkarni

Abstract: Sheth-Tormen mass function has been widely used to quantify the abundance of dark matter halos. It is a significant improvement over the Press-Schechter mass function as it uses ellipsoidal collapse in place of spherical collapse. Both of these mass functions can be written in a form that is universal, i.e., independent of cosmology and power spectrum when scaled in suitable variables. However, co… ▽ More Sheth-Tormen mass function has been widely used to quantify the abundance of dark matter halos. It is a significant improvement over the Press-Schechter mass function as it uses ellipsoidal collapse in place of spherical collapse. Both of these mass functions can be written in a form that is universal, i.e., independent of cosmology and power spectrum when scaled in suitable variables. However, cosmological simulations have shown that this universality is approximate. In this paper, we investigate the power spectrum dependence of halo mass function through a suite of dark-matter-only N-body simulations of seven power-law models in an Einstein-de Sitter cosmology. This choice of cosmology and a power-law power spectrum ensures the self-similar evolution of dark matter distribution, allowing us to isolate the power spectrum dependence of mass function. We find that the mass function shows a clear non-universality. We present fits for the parameters of the Sheth-Tormen mass function for a range of power-law power-spectrum indices. We find a mild evolution in the overall shape of the mass function with the epoch. Finally, we extend our result to LCDM cosmology. We show that the Sheth-Tormen mass function with parameter values derived from a matched power-law EdS cosmology provides a better fit to the LCDM mass function than the standard Sheth-Tormen mass function. Our results indicate that an improved analytical theory is required to provide better fits to the mass function. △ Less

Submitted 3 April, 2023; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 12 pages, 10 figures. This is a much expanded and upgraded version of 0908.2702. Accepted for Publication in the MNRAS

FUMES III: Ultraviolet and Optical Variability of M Dwarf Chromospheres

Authors: Girish M. Duvvuri, J. Sebastian Pineda, Zachory K. Berta-Thompson, Kevin France, Allison Youngblood

Abstract: We obtained ultraviolet and optical spectra for 9 M~dwarfs across a range of rotation periods to determine whether they showed stochastic intrinsic variability distinguishable from flares. The ultraviolet spectra were observed during the Far Ultraviolet M~Dwarf Evolution Survey \emph{Hubble~Space~Telescope} program using the Space Telescope Imaging Spectrograph. The optical observations were taken… ▽ More We obtained ultraviolet and optical spectra for 9 M~dwarfs across a range of rotation periods to determine whether they showed stochastic intrinsic variability distinguishable from flares. The ultraviolet spectra were observed during the Far Ultraviolet M~Dwarf Evolution Survey \emph{Hubble~Space~Telescope} program using the Space Telescope Imaging Spectrograph. The optical observations were taken from the Apache Point Observatory 3.5-meter telescope using the Dual Imaging Spectrograph and from the Gemini South Observatory using the Gemini Multi-Object Spectrograph. We used the optical spectra to measure multiple chromospheric lines: the Balmer series from H$α$ to H$10$ and the \ion{Ca}{2}~H and K lines. We find that after excising flares, these lines vary on the order of $1-20\%$ at minute-cadence over the course of an hour. The absolute amplitude of variability was greater for the faster rotating M~dwarfs in our sample. Among the 5 stars for which measured the weaker Balmer lines, we note a tentative trend that the fractional amplitude of the variability increases for higher order Balmer lines. We measured the integrated flux of multiple ultraviolet emission features formed in the transition region: the \ion{N}{5}, \ion{Si}{4} and \ion{C}{4} resonance line doublets, and the \ion{C}{2} and \ion{He}{2} multiplets. The signal-to-noise (S/N) ratio of the UV data was too low for us to detect non-flare variability at the same scale and time cadence as the optical. We consider multiple mechanisms for the observed stochastic variability and propose both observational and theoretical avenues of investigation to determine the physical causes of intrinsic variability in the chromospheres of M~dwarfs. △ Less

Submitted 17 October, 2022; originally announced October 2022.

X-ray Emission from the Exoplanet Hosting LTT 1445 Triple Star System

Authors: Alexander Brown, Cynthia S. Froning, Allison Youngblood, Kevin France, David J. Wilson, Girish Duvvuri, Yamila Miguel, Hannah Diamond-Lowe

Abstract: JWST will be able to observe the atmospheres of rocky planets transiting nearby M dwarfs. The M dwarf triple star system LTT 1445, at a distance of 6.86 pc, hosts some of the nearest rocky terrestrial planets. These planets most likely orbit the M 3.5V star LTT 1445A. During a 28.6 ksec Chandra ACIS-S3 observation we have i) spatially resolved and detected all three stars in the LTT 1445 system, i… ▽ More JWST will be able to observe the atmospheres of rocky planets transiting nearby M dwarfs. The M dwarf triple star system LTT 1445, at a distance of 6.86 pc, hosts some of the nearest rocky terrestrial planets. These planets most likely orbit the M 3.5V star LTT 1445A. During a 28.6 ksec Chandra ACIS-S3 observation we have i) spatially resolved and detected all three stars in the LTT 1445 system, ii) measured the X-ray luminosity of the individual stars, including LTT 1445A, for the first time, iii) studied the flux variability of the X-ray sources and found strong variability from the A and C components, and iv) investigated how the coronal luminosities, temperatures and volume emission measures vary at different activity levels. Combining these X-ray data with upcoming HST ultraviolet observations will allow a differential emission measure (DEM) estimation of the star's EUV spectrum, thereby facilitating modeling of the rocky planets' atmospheres. △ Less

Submitted 22 September, 2022; originally announced September 2022.

Comments: 18 pages, 4 figures. Accepted to The Astronomical Journal (2022 September 1)

The need for obscured supermassive black hole growth to explain quasar proximity zones in the epoch of reionization

Authors: Sindhu Satyavolu, Girish Kulkarni, Laura C. Keating, Martin G. Haehnelt

Abstract: Proximity zones of high-redshift quasars are unique probes of supermassive black hole formation, but simultaneously explaining proximity zone sizes and black hole masses has proved to be challenging. We study the robustness of some of the assumptions that are usually made to infer quasar lifetimes from proximity zone sizes. We show that small proximity zones can be readily explained by quasars tha… ▽ More Proximity zones of high-redshift quasars are unique probes of supermassive black hole formation, but simultaneously explaining proximity zone sizes and black hole masses has proved to be challenging. We study the robustness of some of the assumptions that are usually made to infer quasar lifetimes from proximity zone sizes. We show that small proximity zones can be readily explained by quasars that vary in brightness with a short duty cycle of $f_\mathrm{duty}\sim 0.1$ and short bright periods of $t_\mathrm{on}\sim 10^4$ yr, even for long lifetimes. We further show that reconciling this with black hole mass estimates requires the black hole to continue to grow and accrete during its obscured phase. The consequent obscured fractions of $\gtrsim$ 0.7 or higher are consistent with low-redshift measurements and models of black hole accretion. Such short duty cycles and long obscured phases are also consistent with observations of large proximity zones, thus providing a simple, unified model for proximity zones of all sizes. The large dynamic range of our simulation, and its calibration to the Lyman-$α$ forest, allows us to investigate the influence of the large-scale topology of reionization and the quasar's host halo mass on proximity zones. We find that incomplete reionization can impede the growth of proximity zones and make them smaller up to 30%, but the quasar host halo mass only affects proximity zones weakly and indirectly. Our work suggests that high-redshift proximity zones can be an effective tool to study quasar variability and black hole growth. △ Less

Submitted 17 March, 2023; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 19 pages, 17 figures, Accepted for publication in MNRAS

Report number: TIFR/TH/22-38

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/

High-Performance Computing for SKA Transient Search: Use of FPGA based Accelerators -- a brief review

Authors: R. Aafreen, R. Abhishek, B. Ajithkumar, Arunkumar M. Vaidyanathan, Indrajit. V. Barve, Sahana Bhattramakki, Shashank Bhat, B. S. Girish, Atul Ghalame, Y. Gupta, Harshal G. Hayatnagarkar, P. A. Kamini, A. Karastergiou, L. Levin, S. Madhavi, M. Mekhala, M. Mickaliger, V. Mugundhan, Arun Naidu, J. Oppermann, B. Arul Pandian, N. Patra, A. Raghunathan, Jayanta Roy, Shiv Sethi , et al. (12 additional authors not shown)

Abstract: This paper presents the High-Performance computing efforts with FPGA for the accelerated pulsar/transient search for the SKA. Case studies are presented from within SKA and pathfinder telescopes highlighting future opportunities. It reviews the scenario that has shifted from offline processing of the radio telescope data to digitizing several hundreds/thousands of antenna outputs over huge bandwid… ▽ More This paper presents the High-Performance computing efforts with FPGA for the accelerated pulsar/transient search for the SKA. Case studies are presented from within SKA and pathfinder telescopes highlighting future opportunities. It reviews the scenario that has shifted from offline processing of the radio telescope data to digitizing several hundreds/thousands of antenna outputs over huge bandwidths, forming several 100s of beams, and processing the data in the SKA real-time pulsar search pipelines. A brief account of the different architectures of the accelerators, primarily the new generation Field Programmable Gate Array-based accelerators, showing their critical roles to achieve high-performance computing and in handling the enormous data volume problems of the SKA is presented here. It also presents the power-performance efficiency of this emerging technology and presents potential future scenarios. △ Less

Submitted 17 January, 2023; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: Accepted for publication in the JoAA special issue on SKA (2023)

AU Microscopii in the FUV: Observations in Quiescence, During Flares, and Implications for AU Mic b and c

Authors: Adina D. Feinstein, Kevin France, Allison Youngblood, Girish M. Duvvuri, DJ Teal, P. Wilson Cauley, Darryl Z. Seligman, Eric Gaidos, Eliza M. R. Kempton, Jacob L. Bean, Hannah Diamond-Lowe, Elisabeth Newton, Sivan Ginzburg, Peter Plavchan, Peter Gao, Hilke Schlichting

Abstract: High energy X-ray and ultraviolet (UV) radiation from young stars impacts planetary atmospheric chemistry and mass loss. The active $\sim 22$ Myr M dwarf AU Mic hosts two exoplanets orbiting interior to its debris disk. Therefore, this system provides a unique opportunity to quantify the effects of stellar XUV irradiation on planetary atmospheres as a function of both age and orbital separation. I… ▽ More High energy X-ray and ultraviolet (UV) radiation from young stars impacts planetary atmospheric chemistry and mass loss. The active $\sim 22$ Myr M dwarf AU Mic hosts two exoplanets orbiting interior to its debris disk. Therefore, this system provides a unique opportunity to quantify the effects of stellar XUV irradiation on planetary atmospheres as a function of both age and orbital separation. In this paper we present over 5 hours of Far-UV (FUV) observations of AU Mic taken with the Cosmic Origins Spectrograph (COS; 1070-1360 Angstrom) on the Hubble Space Telescope (HST). We provide an itemization of $120$ emission features in the HST/COS FUV spectrum and quantify the flux contributions from formation temperatures ranging from $10^4-10^7$ K. We detect 13flares in the FUV white-light curve with energies ranging from $10^{29} - 10^{31}$ ergs. The majority of the energy in each of these flares is released from the transition region between the chromosphere and the corona. There is a 100$\times$ increase in flux at continuum wavelengths $λ< 1100$ Angstrom in each flare which may be caused by thermal Bremsstrahlung emission. We calculate that the baseline atmospheric mass-loss rate for AU Mic b is $\sim 10^8$ g s$^{-1}$, although this rate can be as high as $\sim 10^{14}$ g s$^{-1}$ during flares with $L_\textrm{flare} \simeq 10^{33}$ erg s$^{-1}$. Finally, we model the transmission spectra for AU Mic b and c with a new panchromatic spectrum of AU Mic c and motivate future JWST observations of these planets. △ Less

Submitted 12 August, 2022; v1 submitted 19 May, 2022; originally announced May 2022.

Comments: 29 pages, 12 figures, 3 tables, Accepted to AJ

Long Dark Gaps in the Ly$β$ Forest at $z<6$: Evidence of Ultra Late Reionization from XQR-30 Spectra

Authors: Yongda Zhu, George D. Becker, Sarah E. I. Bosman, Laura C. Keating, Valentina D'Odorico, Rebecca L. Davies, Holly M. Christenson, Eduardo Bañados, Fuyan Bian, Manuela Bischetti, Huanqing Chen, Frederick B. Davies, Anna-Christina Eilers, Xiaohui Fan, Prakash Gaikwad, Bradley Greig, Martin G. Haehnelt, Girish Kulkarni, Samuel Lai, Andrea Pallottini, Yuxiang Qin, Emma Ryan-Weber, Fabian Walter, Feige Wang, Jinyi Yang

Abstract: We present a new investigation of the intergalactic medium (IGM) near reionization using dark gaps in the Lyman-$β$ (Ly$β$) forest. With its lower optical depth, Ly$β$ offers a potentially more sensitive probe to any remaining neutral gas compared to commonly used Ly$α$ line. We identify dark gaps in the Ly$β$ forest using spectra of 42 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30… ▽ More We present a new investigation of the intergalactic medium (IGM) near reionization using dark gaps in the Lyman-$β$ (Ly$β$) forest. With its lower optical depth, Ly$β$ offers a potentially more sensitive probe to any remaining neutral gas compared to commonly used Ly$α$ line. We identify dark gaps in the Ly$β$ forest using spectra of 42 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30 VLT Large Programme. Approximately $40\%$ of these QSO spectra exhibit dark gaps longer than $10h^{-1}{\rm Mpc}$ at $z\simeq5.8$. By comparing the results to predictions from simulations, we find that the data are broadly consistent both with models where fluctuations in the Ly$α$ forest are caused solely by ionizing ultraviolet background (UVB) fluctuations and with models that include large neutral hydrogen patches at $z<6$ due to a late end to reionization. Of particular interest is a very long ($L=28h^{-1}{\rm Mpc}$) and dark ($τ_{\rm eff} \gtrsim 6$) gap persisting down to $z\simeq 5.5$ in the Ly$β$ forest of the $z_{\rm}=5.85$ QSO PSO J025$-$11. This gap may support late reionization models with a volume-weighted average neutral hydrogen fraction of $ \langle x_{\rm HI}\rangle \gtrsim 5\%$ by $z=5.6$. Finally, we infer constraints on $\langle x_{\rm HI}\rangle$ over $5.5 \lesssim z \lesssim 6.0$ based on the observed Ly$β$ dark gap length distribution and a conservative relationship between gap length and neutral fraction derived from simulations. We find $\langle x_{\rm HI}\rangle \leq 0.05$, 0.17, and 0.29 at $z\simeq 5.55$, 5.75, and 5.95, respectively. These constraints are consistent with models where reionization ends significantly later than $z = 6$. △ Less

Submitted 9 May, 2022; originally announced May 2022.

Comments: 17 pages, 13 figures, accepted for publication in ApJ. The Figure set and machine readable table are available at https://ydzhuastro.github.io/lyb.html

Implications of the cosmological 21-cm absorption profile for high-redshift star formation and deep JWST surveys

Authors: Shikhar Mittal, Girish Kulkarni

Abstract: Apart from its anomalously large depth, the cosmological 21-cm absorption signal measured by the EDGES collaboration also has a shape that is distinctly different from theoretical predictions. Models with non-traditional components such as super-adiabatic baryonic cooling or an excess radio background explain the depth of the observed profile, but still conspicuously fail to explain its shape. In… ▽ More Apart from its anomalously large depth, the cosmological 21-cm absorption signal measured by the EDGES collaboration also has a shape that is distinctly different from theoretical predictions. Models with non-traditional components such as super-adiabatic baryonic cooling or an excess radio background explain the depth of the observed profile, but still conspicuously fail to explain its shape. In this paper, we quantify the requirements imposed by the EDGES measurement on sources of Ly $α$ and X-ray photons in the presence of excess radio background at cosmic dawn. In extreme cases, the Ly $α$ and X-ray emissivities require to be enhanced by up to an order of magnitude relative to traditional models. Furthermore, this enhancement needs to be active only for a short duration. We find that under conventional assumptions for the cosmic star formation rate density, standard stellar populations are incapable of meeting these conditions. Only highly unusual models of massive metal-free stars seem to provide a possible mechanism. Conversely, if the sources of Ly $α$ and X-ray photons are compelled to have standard properties, the EDGES measurement puts strong demands on the cosmic star formation rate density. This provides interesting falsifiable predictions for high-redshift galaxy surveys enabled by \textit{James Webb Space Telescope} (\textit{JWST}). We derive predictions for galaxy UV luminosity functions and number densities, and show that a deep \textit{JWST} survey with a limiting UV magnitude of $m_\mathrm{UV,lim}=32$ would potentially be able to rule out the predictions enforced by the EDGES measurement. △ Less

Submitted 2 August, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: Published in MNRAS

Report number: TIFR/TH/22-14

Journal ref: MNRAS 515 (2022) 2901

The Lick Observatory Supernova Search follow-up program: photometry data release of 70 stripped-envelope supernovae

Authors: WeiKang Zheng, Benjamin E. Stahl, Thomas de Jaeger, Alexei V. Filippenko, Shan-Qin Wang, Wen-Pei Gan, Thomas G. Brink, Ivan Altunin, Raphael Baer-Way, Andrew Bigley, Kyle Blanchard, Peter K. Blanchard, James Bradley, Samantha K. Cargill, Chadwick Casper, Teagan Chapman, Vidhi Chander, Sanyum Channa, Byung Yun Choi, Nick Choksi, Matthew Chu, Kelsey I. Clubb, Daniel P. Cohen, Paul A. Dalba, Asia deGraw , et al. (63 additional authors not shown)

Abstract: We present BVRI and unfiltered Clear light curves of 70 stripped-envelope supernovae (SESNe), observed between 2003 and 2020, from the Lick Observatory Supernova Search (LOSS) follow-up program. Our SESN sample consists of 19 spectroscopically normal SNe~Ib, two peculiar SNe Ib, six SN Ibn, 14 normal SNe Ic, one peculiar SN Ic, ten SNe Ic-BL, 15 SNe IIb, one ambiguous SN IIb/Ib/c, and two superlum… ▽ More We present BVRI and unfiltered Clear light curves of 70 stripped-envelope supernovae (SESNe), observed between 2003 and 2020, from the Lick Observatory Supernova Search (LOSS) follow-up program. Our SESN sample consists of 19 spectroscopically normal SNe~Ib, two peculiar SNe Ib, six SN Ibn, 14 normal SNe Ic, one peculiar SN Ic, ten SNe Ic-BL, 15 SNe IIb, one ambiguous SN IIb/Ib/c, and two superluminous SNe. Our follow-up photometry has (on a per-SN basis) a mean coverage of 81 photometric points (median of 58 points) and a mean cadence of 3.6d (median of 1.2d). From our full sample, a subset of 38 SNe have pre-maximum coverage in at least one passband, allowing for the peak brightness of each SN in this subset to be quantitatively determined. We describe our data collection and processing techniques, with emphasis toward our automated photometry pipeline, from which we derive publicly available data products to enable and encourage further study by the community. Using these data products, we derive host-galaxy extinction values through the empirical colour evolution relationship and, for the first time, produce accurate rise-time measurements for a large sample of SESNe in both optical and infrared passbands. By modeling multiband light curves, we find that SNe Ic tend to have lower ejecta masses and lower ejecta velocities than SNe~Ib and IIb, but higher $^{56}$Ni masses. △ Less

Submitted 10 March, 2022; originally announced March 2022.

Comments: Accepted by MNRAS

On the detection of a cosmic dawn signal in the radio background

Authors: Saurabh Singh, Jishnu Nambissan T., Ravi Subrahmanyan, N. Udaya Shankar, B. S. Girish, A. Raghunathan, R. Somashekar, K. S. Srivani, Mayuri Sathyanarayana Rao

Abstract: The astrophysics of cosmic dawn, when star formation commenced in the first collapsed objects, is predicted to be revealed as spectral and spatial signatures in the cosmic radio background at long wavelengths. The sky-averaged redshifted 21-cm absorption line of neutral hydrogen is a probe of cosmic dawn. The line profile is determined by the evolving thermal state of the gas, radiation background… ▽ More The astrophysics of cosmic dawn, when star formation commenced in the first collapsed objects, is predicted to be revealed as spectral and spatial signatures in the cosmic radio background at long wavelengths. The sky-averaged redshifted 21-cm absorption line of neutral hydrogen is a probe of cosmic dawn. The line profile is determined by the evolving thermal state of the gas, radiation background, Lyman-$α$ radiation from stars scattering off cold primordial gas and the relative populations of the hyperfine spin levels in neutral hydrogen atoms. We report a radiometer measurement of the spectrum of the radio sky in the 55--85~MHz band, which shows that the profile found by Bowman et al. in data taken with the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) low-band instrument is not of astrophysical origin; their best-fitting profile is rejected with 95.3\% confidence. The profile was interpreted to be a signature of cosmic dawn; however, its amplitude was substantially higher than that predicted by standard cosmological models. Explanations for the amplitude of the profile included non-standard cosmology, additional mechanisms for cooling the baryons, perhaps via interactions with millicharged dark matter and an excess radio background at redshifts beyond 17. Our non-detection bears out earlier concerns and suggests that the profile found by Bowman et al. is not evidence for new astrophysics or non-standard cosmology. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Comments: Accepted in principle

Background of radio photons from primordial black holes

Authors: Shikhar Mittal, Girish Kulkarni

Abstract: We compute the isotropic radiation background due to Hawking emission from primordial black holes (PBHs), and examine if this background is a viable option in explaining the excess radiowave background observed by the ARCADE2 and LWA1 experiments at $\lesssim 1\,$GHz. We find that even under the extreme assumption that all of the dark matter is in the form of PBHs, the radio brightness temperature… ▽ More We compute the isotropic radiation background due to Hawking emission from primordial black holes (PBHs), and examine if this background is a viable option in explaining the excess radiowave background observed by the ARCADE2 and LWA1 experiments at $\lesssim 1\,$GHz. We find that even under the extreme assumption that all of the dark matter is in the form of PBHs, the radio brightness temperature induced by Hawking evaporation of PBHs is $\mathcal{O}(10^{-46})\,$K, highly subdominant compared to the cosmic microwave background. The main reason for this is that for PBHs in the mass range $\sim10^{12}$-$10^{14}\,$kg, which can be constrained by Hawking emission, the spectrum peaks at $10^7$ to $10^5\,$eV. As the Hawking spectrum is power law suppressed towards lower energies, negligible flux of $μ$eV photons is obtained. The peak of the Hawking spectrum shifts to lower energies for higher masses, but the number density is low and so is the specific intensity. Because Hawking emission from PBHs is thus unable to explain the observed excess radio background, we also consider the alternative possibility of radio emission from gas accretion onto supermassive PBHs. These PBHs can readily produce strong radio emission that could easily explain the ARCADE2/LWA1 excess. △ Less

Submitted 20 January, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: Published

Report number: TIFR/TH/21-16

Journal ref: MNRAS 510 (2022) 4992

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

Chasing the Tail of Cosmic Reionization with Dark Gap Statistics in the Ly$α$ Forest over $5 < z < 6$

Authors: Yongda Zhu, George D. Becker, Sarah E. I. Bosman, Laura C. Keating, Holly M. Christenson, Eduardo Bañados, Fuyan Bian, Frederick B. Davies, Valentina D'Odorico, Anna-Christina Eilers, Xiaohui Fan, Martin G. Haehnelt, Girish Kulkarni, Andrea Pallottini, Yuxiang Qin, Feige Wang, Jinyi Yang

Abstract: We present a new investigation of the intergalactic medium (IGM) near the end of reionization using "dark gaps" in the Lyman-alpha (Ly$α$) forest. Using spectra of 55 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30 VLT Large Programme, we identify gaps in the Ly$α$ forest where the transmission averaged over 1 comoving $h^{-1}\,{\rm Mpc}$ bins falls below 5%. Nine ultra-long (… ▽ More We present a new investigation of the intergalactic medium (IGM) near the end of reionization using "dark gaps" in the Lyman-alpha (Ly$α$) forest. Using spectra of 55 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30 VLT Large Programme, we identify gaps in the Ly$α$ forest where the transmission averaged over 1 comoving $h^{-1}\,{\rm Mpc}$ bins falls below 5%. Nine ultra-long ($L > 80~h^{-1}\,{\rm Mpc}$) dark gaps are identified at $z<6$. In addition, we quantify the fraction of QSO spectra exhibiting gaps longer than $30~h^{-1}\,{\rm Mpc}$, $F_{30}$, as a function of redshift. We measure $F_{30} \simeq 0.9$, 0.6, and 0.15 at $z = 6.0$, 5.8, and 5.6, respectively, with the last of these long dark gaps persisting down to $z \simeq 5.3$. Comparing our results with predictions from hydrodynamical simulations, we find that the data are consistent with models wherein reionization extends significantly below redshift six. Models wherein the IGM is essentially fully reionized that retain large-scale fluctuations in the ionizing UV background at $z \lesssim 6$ are also potentially consistent with the data. Overall, our results suggest that signature of reionization in the form of islands of neutral hydrogen and/or large-scale fluctuations in the ionizing background remain present in the IGM until at least $z \simeq 5.3$. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 30 pages, 19 figures, 3 tables, accepted for publication in ApJ

Hydrogen reionisation ends by $z=5.3$: Lyman-$α$ optical depth measured by the XQR-30 sample

Authors: Sarah E. I. Bosman, Frederick B. Davies, George D. Becker, Laura C. Keating, Rebecca L. Davies, Yongda Zhu, Anna-Christina Eilers, Valentina D'Odorico, Fuyan Bian, Manuela Bischetti, Stefano V. Cristiani, Xiaohui Fan, Emanuele P. Farina, Martin G. Haehnelt, Joseph F. Hennawi, Girish Kulkarni, Andrei Mesinger, Romain A. Meyer, Masafusa Onoue, Andrea Pallottini, Yuxiang Qin, Emma Ryan-Weber, Jan-Torge Schindler, Fabian Walter, Feige Wang , et al. (1 additional authors not shown)

Abstract: The presence of excess scatter in the Ly-$α$ forest at $z\sim 5.5$, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionisation. However, low data quality and systematic uncertainties complicate the use of Ly-$α$ transmission as a precision probe of reionisation's end stages. In this paper, we assembl… ▽ More The presence of excess scatter in the Ly-$α$ forest at $z\sim 5.5$, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionisation. However, low data quality and systematic uncertainties complicate the use of Ly-$α$ transmission as a precision probe of reionisation's end stages. In this paper, we assemble a sample of 67 quasar sightlines at $z>5.5$ with high signal-to-noise ratios of $>10$ per $\leq 15$ km s$^{-1}$ spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR $>20$ per pixel) spectra of 30 quasars at $z>5.7$. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-$α$ systems. We present improved measurements of the mean Ly-$α$ transmission over $4.9<z<6.1$. Using all known systematics in a forward modelling analysis, we find excellent agreement between the observed Ly-$α$ transmission distributions and the homogeneous-UVB simulations Sherwood and Nyx up to $z\leq5.2$ ($<1 σ$), and mild tension ($\sim 2.5 σ$) at $z=5.3$. Homogeneous UVB models are ruled out by excess Ly-$α$ transmission scatter at $z\geq5.4$ with high confidence ($>3.5 σ$). Our results indicate that reionisation-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least $z=5.3$ ($t=1.1$ Gyr after the Big Bang). This is further evidence for a late end to reionisation. △ Less

Submitted 11 April, 2022; v1 submitted 8 August, 2021; originally announced August 2021.

Comments: Accepted for publication in MNRAS. Full measurement datasets are available on the first author's website

PIPS, an advanced platform for period detection in time series -- I. Fourier-likelihood periodogram and application to RR Lyrae Stars

Authors: Yukei S. Murakami, Connor Jennings, Andrew M. Hoffman, Arjun B. Savel, James Sunseri, Raphael Baer-Way, Benjamin E. Stahl, Ivan Altunin, Nachiket Girish, Alexei V. Filippenko

Abstract: We describe the $\texttt{Period detection and Identification Pipeline Suite}$ ($\texttt{PIPS}$) -- a new, fast, and statistically robust platform for period detection and analysis of astrophysical time-series data. $\texttt{PIPS}$ is an open-source Python package that provides various pre-implemented methods and a customisable framework for automated, robust period measurements with principled unc… ▽ More We describe the $\texttt{Period detection and Identification Pipeline Suite}$ ($\texttt{PIPS}$) -- a new, fast, and statistically robust platform for period detection and analysis of astrophysical time-series data. $\texttt{PIPS}$ is an open-source Python package that provides various pre-implemented methods and a customisable framework for automated, robust period measurements with principled uncertainties and statistical significance calculations. In addition to detailing the general algorithm that underlies $\texttt{PIPS}$, this paper discusses one of $\texttt{PIPS'}$ central and novel features, the Fourier-likelihood periodogram, and compares its performance to existing methods. The resulting improved performance implies that one can construct deeper, larger, and more reliable sets of derived properties from various observations, including all-sky surveys. We present a comprehensive validation of $\texttt{PIPS}$ against artificially generated data, which demonstrates the reliable performance of our algorithm for a class of periodic variable stars (RR Lyrae stars). △ Less

Submitted 1 June, 2022; v1 submitted 29 July, 2021; originally announced July 2021.

Comments: 18 pages, 13 figures. Accepted to MNRAS

Constraining primordial black holes as dark matter using the global 21-cm signal with X-ray heating and excess radio background

Authors: Shikhar Mittal, Anupam Ray, Girish Kulkarni, Basudeb Dasgupta

Abstract: Using the global 21-cm signal measurement by the EDGES collaboration, we derive constraints on the fraction of the dark matter that is in the form of primordial black holes (PBHs) with masses in the range $10^{15}$-$10^{17}\,$g. Improving upon previous analyses, we consider the effect of the X-ray heating of the intergalactic medium on these constraints, and also use the full shape of the 21-cm ab… ▽ More Using the global 21-cm signal measurement by the EDGES collaboration, we derive constraints on the fraction of the dark matter that is in the form of primordial black holes (PBHs) with masses in the range $10^{15}$-$10^{17}\,$g. Improving upon previous analyses, we consider the effect of the X-ray heating of the intergalactic medium on these constraints, and also use the full shape of the 21-cm absorption feature in our inference. In order to account for the anomalously deep absorption amplitude, we also consider an excess radio background motivated by LWA1 and ARCADE2 observations. Because the heating rate induced by PBH evaporation evolves slowly, the data favour a scenario in which PBH-induced heating is accompanied by X-ray heating. Also, for the same reason, using the full measurement across the EDGES observation band yields much stronger constraints on PBHs than just the redshift of absorption. We find that 21-cm observations exclude $f_{\mathrm{PBH}} \gtrsim 10^{-9.7}$ at 95% CL for $M_{\mathrm{PBH}}=10^{15}\,$g. This limit weakens approximately as $M_{\mathrm{PBH}}^4$ towards higher masses, thus providing the strongest constraints on ultralight evaporating PBHs as dark matter over the entire mass range $10^{15}$-$10^{17}\,$g. Under the assumption of a simple spherical gravitational collapse based on the Press-Schechter formalism, we also derive bounds on the curvature power spectrum at extremely small scales ($k\sim 10^{15}\,$Mpc$^{-1}$). This highlights the usefulness of global 21-cm measurements, including non-detections, across wide frequency bands for probing exotic physical processes. △ Less

Submitted 15 March, 2022; v1 submitted 5 July, 2021; originally announced July 2021.

Comments: Published in JCAP

Report number: TIFR/TH/21-6

Journal ref: JCAP 03 (2022) 030