Researcher profile

J. Xavier Prochaska

J. Xavier Prochaska contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
49works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2026arXiv

Featureless stars: Flux Calibration for Extremely Large Telescopes

The spectrophotometric flux calibration of recent spectroscopic surveys has reached a limiting systematic precision of approximately 1-3 percent, and is often biased near the wavelengths associated with H I Balmer absorption. As we prepare for the next generation of imaging and spectroscopic surveys, and high-precision cosmology experiments, we must find a way to address this systematic. Towards this goal, we have identified a global network of 29 bright (G < 17.5) featureless white dwarf stars that have a spectral energy distribution consistent with an almost pure blackbody form over the entire optical and near-infrared wavelength range. Based on this sample, we have computed the systematic uncertainty and AB magnitude offsets associated with Gaia, SDSS, SMSS, PanSTARRS, DES, and 2MASS, and we have also checked the consistency of our objects with both GALEX and WISE. The magnitude range of the featureless stars reported here are ideally suited to observations taken with the forthcoming generation of extremely large telescopes, as well as calibrating the survey data acquired by the Rubin, Euclid and Roman observatories. Finally, all of the high-precision spectrophotometric standard stars reported here have been included in the latest release of the PypeIt data reduction pipeline.

preprint2026arXiv

The Second CHIME/FRB Catalog of Fast Radio Bursts

We present a catalog of 4539 fast radio bursts (FRBs) observed with the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope between 25 July 2018 and 15 September 2023. These bursts originate from 3641 unique sources, including 981 bursts from 83 known repeating sources. For each FRB, the catalog provides a $O(10&#39;)$ estimate of sky location along with corresponding measurements of cumulative exposure time and survey sensitivity over the observing period. It includes a total-intensity dynamic spectrum between 400 and 800 MHz at 0.983 ms resolution. From this spectrum, we constrain a model of the burst morphology and measure key parameters such as arrival time, intrinsic temporal width, dispersion measure, scattering time, and flux density. This second catalog includes all FRBs from the first catalog, with every event reprocessed using a uniform and improved analysis framework. We show that previously published inferences remain valid under the updated measurements. We assess consistency of the detection rate across observational parameters, present initial distributions of burst properties, and outline ongoing and future studies that will use this catalog to investigate the nature of FRBs and their utility as astrophysical and cosmological probes.

preprint2025arXiv

Localisation and host galaxy identification of new Fast Radio Bursts with MeerKAT

Accurately localising fast radio bursts (FRBs) is essential for understanding their birth environments and for their use as cosmological probes. Recent advances in radio interferometry, particularly with MeerKAT, have enabled the localisation of individual bursts with arcsecond precision. In this work, we present the localisation of 15 apparently non-repeating FRBs detected with MeerKAT. Two of the FRBs, discovered in 2022, were localised in 8 second images from the projects which MeerTRAP was commensal to, while eight were localised using the transient buffer (TB) pipeline, and another one through SeeKAT, all with arcsecond precision. Four additional FRBs lacked TB triggers and sufficient signal, limiting their localisation only to arcminute precision. For eight of the FRBs in our sample, we identify host galaxies with greater than 90% confidence, and one with 80% confidence, while two FRBs have ambiguous associations. We measured spectroscopic redshifts for six host galaxies, ranging from 0.33 to 0.85, demonstrating MeerKAT&#39;s sensitivity to high redshift FRBs. We modelled the spectral energy distributions of host galaxies with sufficient photometric coverage to derive their stellar population and star formation properties. This work represents one of the largest uniform samples of well-localised distant FRBs to date, laying the groundwork for using MeerKAT FRBs as cosmological probes and understand how FRB hosts evolve at high redshift.

preprint2023arXiv

SDSS DR17: The Cosmic Slime Value Added Catalog

The &#34;cosmic web&#34;, the filamentary large-scale structure in a cold dark matter Universe, is readily apparent via galaxy tracers in spectroscopic surveys. However, the underlying dark matter structure is as of yet unobservable and mapping the diffuse gas permeating it lies beyond practical observational capabilities. A recently developed technique, inspired by the growth and movement of Physarum polycephalum &#34;slime mold&#34;, has been used to map the cosmic web of a low redshift sub-sample of the SDSS spectroscopic galaxy catalog. This model, the Monte Carlo Physarum Machine (MCPM) was shown to promisingly reconstruct the cosmic web. Here, we improve the formalism used in calibrating the MCPM to better recreate the Bolshoi-Planck cosmological simulation&#39;s density distributions and apply them to a significantly larger cosmological volume than previous works using the Sloan Digital Sky Survey (SDSS, $z < 0.1$) and the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) Luminous Red Galaxy (LRG, $z \lesssim 0.5$) spectroscopic catalogs. We present the &#34;Cosmic Slime Value Added Catalog&#34; which provides estimates for the cosmic overdensity for the sample of galaxies probed spectroscopically by the above SDSS surveys. In addition, we provide the fully reconstructed 3D density cubes of these volumes. These data products were released as part of Sloan Digital Sky Survey Data Release 17 and are publicly available. We present the input catalogs and the methodology for constructing these data products. We also highlight exciting potential applications to galaxy evolution, cosmology, the intergalactic and circumgalactic medium, and transient phenomenon localization.

preprint2023arXiv

Second Data Release of the COSMOS Lyman-alpha Mapping and Tomographic Observation: The First 3D Maps of the Detailed Cosmic Web at 2.05<z<2.55

We present the second data release of the COSMOS Lyman-Alpha Mapping And Tomography Observations (CLAMATO) Survey conducted with the LRIS spectrograph on the Keck-I telescope. This project used Lyman-alpha forest absorption in the spectra of faint star forming galaxies and quasars at z ~ 2-3 to trace neutral hydrogen in the intergalactic medium. In particular, we use 320 objects over a footprint of ~0.2 deg^2 to reconstruct the absorption field at 2.05 < z < 2.55 at ~2 h^{-1}Mpc resolution. We apply a Wiener filtering technique to the observed data to reconstruct three dimensional maps of the field over a volume of 4.1 x 10^5 comoving cubic Mpc. In addition to the filtered flux maps, for the first time we infer the underlying dark matter field through a forward modeling framework from a joint likelihood of galaxy and Lyman-alpha forest data, finding clear examples of the detailed cosmic web consisting of cosmic voids, sheets, filaments, and nodes. In addition to traditional figures, we present a number of interactive three dimensional models to allow exploration of the data and qualitative comparisons to known galaxy surveys. We find that our inferred over-densities are consistent with those found from galaxy fields. Our reduced spectra, extracted Lyman-alpha forest pixel data, and reconstructed tomographic maps are available publicly at https://doi.org/10.5281/zenodo.7524313

preprint2022arXiv

A fast radio burst progenitor born in a galaxy merger

We report a Giant Metrewave Radio Telescope HI 21cm mapping study of the neutral atomic hydrogen (HI) in the host galaxy of the fast radio burst (FRB) FRB20180916B at $z \approx 0.03399$. We find that the FRB host has an HI mass of $\rm M_{HI} = (2.74 \pm 0.33) \times 10^9 \ M_\odot$ and a high HI-to-stellar mass ratio, $\approx 1.3$. The FRB host is thus a gas-rich but near-quiescent galaxy, that is likely to have acquired a significant mass of HI in the recent past. The HI distribution is disturbed, with extended HI 21cm emission detected in a north-eastern tail, a counter-tail towards the south, an HI hole between the galaxy centre and the FRB location, and a high HI column density measured close to the FRB position. The FRB host is part of a group with four companions detected in their HI 21cm emission, the nearest of which is only 22~kpc from the FRB location. The gas-richness and disturbed HI distribution indicate that the FRB host has recently undergone a minor merger, which increased its HI mass, disturbed the HI in the galaxy disk, and compressed the HI near the FRB location to increase its surface density. We propose that this merger caused the burst of star-formation in the outskirts of the galaxy that gave rise to the FRB progenitor. The evidence for a minor merger is consistent with scenarios in which the FRB progenitor is a massive star, formed due to the merger event.

preprint2022arXiv

A Massive, Dusty HI-Absorption-Selected Galaxy at $z \approx 2.46$ Identified in a CO Emission Survey

We report a NOrthern Extended Millimeter Array (NOEMA) and Atacama Large Millimeter/submillimeter Array (ALMA) search for redshifted CO emission from the galaxies associated with seven high-metallicity ([M/H] $\geq -1.03$) damped Lyman-$α$ absorbers (DLAs) at $z\approx1.64-2.51$. Our observations yielded one new detection of CO(3-2) emission from a galaxy at $z=2.4604$ using NOEMA, associated with the $z=2.4628$ DLA towards QSO B0201+365. Including previous searches, our search results in detection rates of CO emission of $\approx56^{+38}_{-24}$ % and $\approx11^{+26}_{-9}$ %, respectively, in the fields of DLAs with ${\rm [M/H]}>-0.3$ and ${\rm [M/H]}<-0.3$. Further, the HI-selected galaxies associated with five DLAs with [M/H] $>-0.3$ all have high molecular gas masses, $\gtrsim5\times10^{10}\ {\rm M}_\odot$. This indicates that the highest-metallicity DLAs at $z\approx2$ are associated with the most massive galaxies. The newly-identified $z\approx2.4604$ HI-selected galaxy, DLA0201+365g, has an impact parameter of $\approx7$ kpc to the QSO sightline, and an implied molecular gas mass of $(5.04\pm0.78) \times10^{10}\times(α_{\rm CO}/4.36)\times(r_{31}/0.55)\ {\rm M}_\odot$. Archival Hubble Space Telescope Wide Field and Planetary Camera 2 imaging covering the rest-frame near-ultraviolet (NUV) and far-ultraviolet (FUV) emission from this galaxy yield non-detections of rest-frame NUV and FUV emission, and a $5σ$ upper limit of 2.3 M$_\odot$ yr$^{-1}$ on the unobscured star formation rate (SFR). The low NUV-based SFR estimate, despite the very high molecular gas mass, indicates that DLA0201+365g either is a very dusty galaxy, or has a molecular gas depletion time that is around two orders of magnitude larger than that of star-forming galaxies at similar redshifts.

preprint2022arXiv

Constraining the Cosmic Baryon Distribution with Fast Radio Burst Foreground Mapping

The dispersion measures (DM) of fast radio bursts (FRBs) encode the integrated electron density along the line-of-sight, which is typically dominated by the intergalactic medium (IGM) contribution in the case of extragalactic FRBs. In this paper, we show that incorporating wide-field spectroscopic galaxy survey data in the foreground of localized FRBs can significantly improve constraints on the partition of diffuse cosmic baryons. Using mock DMs and realistic lightcone galaxy catalogs derived from the Millennium simulation, we define spectroscopic surveys that can be carried out with 4m and 8m-class wide field spectroscopic facilities. On these simulated surveys, we carry out Bayesian density reconstructions in order to estimate the foreground matter density field. In comparison with the `true&#39; matter density field, we show that these can help reduce the uncertainties in the foreground structures by $\sim 2-3\times$ compared to cosmic variance. We calculate the Fisher matrix to forecast that $N=30\: (96)$ localized FRBs should be able to constrain the diffuse cosmic baryon fraction to $\sim 10\%\: (\sim 5\%) $, and parameters governing the size and baryon fraction of galaxy circumgalactic halos to within $\sim 20-25\%\: (\sim 8-12\%)$. From the Fisher analysis, we show that the foreground data increases the sensitivity of localized FRBs toward our parameters of interest by $\sim 25\times$. We briefly introduce FLIMFLAM, an ongoing galaxy redshift survey that aims to obtain foreground data on $\sim 30$ localized FRB fields.

preprint2022arXiv

Deep Learning of DESI Mock Spectra to Find Damped Lyα Systems

We have updated and applied a convolutional neural network (CNN) machine learning model to discover and characterize damped Ly$α$ systems (DLAs) based on Dark Energy Spectroscopic Instrument (DESI) mock spectra. We have optimized the training process and constructed a CNN model that yields a DLA classification accuracy above 99$\%$ for spectra which have signal-to-noise (S/N) above 5 per pixel. Classification accuracy is the rate of correct classifications. This accuracy remains above 97$\%$ for lower signal-to-noise (S/N) $\approx1$ spectra. This CNN model provides estimations for redshift and HI column density with standard deviations of 0.002 and 0.17 dex for spectra with S/N above 3 per pixel. Also, this DLA finder is able to identify overlapping DLAs and sub-DLAs. Further, the impact of different DLA catalogs on the measurement of Baryon Acoustic Oscillation (BAO) is investigated. The cosmological fitting parameter result for BAO has less than $0.61\%$ difference compared to analysis of the mock results with perfect knowledge of DLAs. This difference is lower than the statistical error for the first year estimated from the mock spectra: above $1.7\%$. We also compared the performance of CNN and Gaussian Process (GP) model. Our improved CNN model has moderately 14$\%$ higher purity and 7$\%$ higher completeness than an older version of GP code, for S/N $>$ 3. Both codes provide good DLA redshift estimates, but the GP produces a better column density estimate by $24\%$ less standard deviation. A credible DLA catalog for DESI main survey can be provided by combining these two algorithms.

preprint2022arXiv

First discoveries and localisations of Fast Radio Bursts with MeerTRAP: a real-time, commensal MeerKAT survey

We report on the discovery and localization of fast radio bursts (FRBs) from the MeerTRAP project, a commensal fast radio transient-detection programme at MeerKAT in South Africa. Our hybrid approach combines a coherent search with an average field-of-view of 0.4 $\rm deg^{2}$ with an incoherent search utilizing a field-of-view of $\sim$1.27 $\rm deg^{2}$ (both at 1284~MHz). Here, we present results on the first three FRBs: FRB 20200413A (DM=1990.05 pc cm$^{-3}$), FRB 20200915A (DM=740.65 pc cm$^{-3}$), and FRB 20201123A (DM=433.55 pc cm$^{-3}$). FRB 20200413A was discovered only in the incoherent beam. FRB 20200915A (also discovered only in the incoherent beam) shows speckled emission in the dynamic spectrum which cannot be explained by interstellar scintillation in our Galaxy or plasma lensing, and might be intrinsic to the source. FRB 20201123A shows a faint post-cursor burst about 200 ms after the main burst and warrants further follow-up to confirm whether it is a repeating FRB. FRB 20201123A also exhibits significant temporal broadening consistent with scattering by a turbulent medium. The broadening exceeds that predicted for medium along the sightline through our Galaxy. We associate this scattering with the turbulent medium in the environment of the FRB in the host galaxy. Within the approximately $1&#39;$ localization region of FRB 20201123A, we identify one luminous galaxy ($r \approx 15.67$; J173438.35$-$504550.4) that dominates the posterior probability for a host association. The galaxy&#39;s measured properties are consistent with other FRB hosts with secure associations.

preprint2022arXiv

HeII Ly$α$ Transmission Spikes and Absorption Troughs in Eight High-resolution Spectra Probing the End of HeII Reionization

We present statistics of HeII Lya transmission spikes and large-scale absorption troughs using archival high-resolution ($R=λ/Δλ\simeq 12,500$-$18,000$) far-UV spectra of eight HeII-transparent quasars obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. The sample covers the redshift range 2.5<z<3.8, thereby probing the rapidly evolving HeII absorption at the end of the HeII reionization epoch. The measured lengths of the troughs decrease dramatically from L>100cMpc at z>3 to L~30cMpc at z~2.7, signaling a significant progression of HeII reionization at these redshifts. Furthermore, unexpectedly long L~65cMpc troughs detected at z~2.9 suggest that the UV background fluctuates at larger scales than predicted by current models. By comparing the measured incidence of transmission spikes to predictions from forward-modeled mock spectra created from the outputs of a (146cMpc)^3 optically thin Nyx hydrodynamical simulation employing different UV background models, we infer the redshift evolution of the HeII photoionization rate $Γ_\mathrm{He\,II}(z)$. The photoionization rate decreases with increasing redshift from $\simeq 4.6\times 10^{-15}\mathrm{\,s^{-1}}$ at z~2.6 to $\simeq 1.2 \times 10^{-15}\mathrm{\,s^{-1}}$ at z~3.2, in agreement with previous inferences from the HeII effective optical depth, and following expected trends of current models of a fluctuating HeII-ionizing background.

preprint2022arXiv

Jansky Very Large Array detections of CO(1-0) emission in HI-absorption-selected galaxies at $z \gtrsim 2$

We report a Karl G. Jansky Very Large Array search for redshifted CO(1-0) emission from three HI-absorption-selected galaxies at $z \approx 2$, identified earlier in their CO(3-2) or CO(4-3) emission. We detect CO(1-0) emission from DLA B1228-113 at $z\approx2.1933$ and DLA J0918+1636 at $z\approx2.5848$; these are the first detections of CO(1-0) emission in high-$z$ HI-selected galaxies. We obtain high molecular gas masses, $\rm M_{mol}\approx10^{11}\times(α_{\rm CO}/4.36)\ M_\odot$, for the two objects with CO(1-0) detections, which are a factor of $\approx1.5-2$ lower than earlier estimates. We determine the excitation of the mid$-J$ CO rotational levels relative to the $J=1$ level, r$_{ J1}$, in HI-selected galaxies for the first time, obtaining r$_{\rm 31}=1.00\pm0.20$ and r$_{\rm 41}=1.03\pm0.23$ for DLA J0918+1636, and r$_{\rm 31}=0.86\pm0.21$ for DLA B1228-113. These values are consistent with thermal excitation of the $J=3,4$ levels. The excitation of the $J=3$ level in the HI-selected galaxies is similar to that seen in massive main-sequence and sub-mm galaxies at $z\gtrsim2$, but higher than that in main-sequence galaxies at $z\approx1.5$; the higher excitation of the galaxies at $z\gtrsim2$ is likely to be due to their higher star-formation rate (SFR) surface density. We use Hubble Space Telescope Wide Field Camera 3 imaging to detect the rest-frame near-ultraviolet emission of DLA B1228-113, obtaining an NUV SFR of $4.44\pm0.47$ M$_{\odot}$ yr$^{-1}$, significantly lower than that obtained from the total infrared luminosity, indicating significant dust extinction in the $z\approx2.1933$ galaxy.

preprint2022arXiv

Monte Carlo Physarum Machine: Characteristics of Pattern Formation in Continuous Stochastic Transport Networks

We present Monte Carlo Physarum Machine: a computational model suitable for reconstructing continuous transport networks from sparse 2D and 3D data. MCPM is a probabilistic generalization of Jones&#39;s 2010 agent-based model for simulating the growth of Physarum polycephalum slime mold. We compare MCPM to Jones&#39;s work on theoretical grounds, and describe a task-specific variant designed for reconstructing the large-scale distribution of gas and dark matter in the Universe known as the Cosmic web. To analyze the new model, we first explore MCPM&#39;s self-patterning behavior, showing a wide range of continuous network-like morphologies -- called &#34;polyphorms&#34; -- that the model produces from geometrically intuitive parameters. Applying MCPM to both simulated and observational cosmological datasets, we then evaluate its ability to produce consistent 3D density maps of the Cosmic web. Finally, we examine other possible tasks where MCPM could be useful, along with several examples of fitting to domain-specific data as proofs of concept.

preprint2021arXiv

A Long Stream of Metal-Poor Cool Gas around a Massive Starburst Galaxy at z = 2.67

We present the first detailed dissection of the circumgalactic medium (CGM) of massive starburst galaxies at z > 2. Our target is a submillimeter galaxy (SMG) at z = 2.674 that has a star formation rate of 1200 $M_\odot$/yr and a molecular gas reservoir of $1.3\times10^{11} M_\odot$. We characterize its CGM with two background QSOs at impact parameters of 93 kpc and 176 kpc. We detect strong HI and metal-line absorption near the redshift of the SMG towards both QSOs, each consisting of three main subsystems spanning over 1500 km/s. The absorbers show remarkable kinematic and metallicity coherence across a separation of 86 kpc. In particular, the cool gas in the CGM of the SMG exhibits high HI column densities ($\log N_{\rm HI}/{\rm cm}^{-2} = 20.2, 18.6$), low metallicities ([M/H] $\approx$ -2.0), and similar radial velocities ($\approx$ -300 km/s). While the HI column densities match previous results on the CGM around QSOs at z > 2, the metallicities are lower by more than an order of magnitude, making it an outlier in the line width$-$metallicity relation of damped Ly$α$ absorbers. The large physical extent, the velocity coherence, the high surface density, and the low metallicity are all consistent with the cool, inflowing, and near-pristine gas streams predicted to penetrate hot massive halos at z > 1.5. We estimate a total gas accretion rate of ~100 $M_\odot$/yr from three such streams, which falls short of the star formation rate but is consistent with simulations. At this rate, it takes about a gigayear to acquire the molecular gas reservoir of the central starburst.

preprint2021arXiv

A Multiwavelength Study of ELAN Environments (AMUSE$^2$): Detection of a dusty star-forming galaxy within the enormous Lyman $α$ nebula at $z=2.3$ sheds light on its origin

We present ALMA observations on and around the radio-quiet quasar UM287 at $z=2.28$. Together with a companion quasar, UM287 is believed to play a major role in powering the surrounding enormous Ly$α$ nebula (ELAN), dubbed the Slug ELAN, that has an end-to-end size of 450 physical kpc. In addition to the quasars, we detect a new dusty star-forming galaxy (DSFG), dubbed the Slug-DSFG, in 2 mm continuum with a single emission line consistent with CO(4-3). The Slug-DSFG sits at a projected distance of 100 kpc south-east from UM287, with a systemic velocity difference of $-360\pm30$ km s$^{-1}$ with respect to UM287, suggesting it being a possible contributor to the powering of the Slug ELAN. With careful modeling of SED and dynamical analyses it is found that the Slug-DSFG and UM287 appear low in both gas fraction and gas-to-dust ratio, suggesting environmental effects due to the host massive halo. In addition, our Keck long-slit spectra reveal significant Ly$α$ emissions from the Slug-DSFG, as well as a Ly$α$ tail that starts at the location and velocity of the Slug-DSFG and extends towards the south, with a projected length of about 100 kpc. Supported by various analytical estimates we propose that the Ly$α$ tail is a result of the Slug-DSFG experiencing ram pressure stripping. The gas mass stripped is estimated to be about 10$^9$ M$_\odot$, contributing to the dense warm/cool gas reservoir that is believed to help power the exceptional Ly$α$ luminosity.

preprint2021arXiv

A Multiwavelength Study of ELAN Environments (AMUSE$^2$): Ubiquitous dusty star-forming galaxies associated with enormous Ly$α$ nebulae on megaparsec scales

We have been undertaking a systematic survey at 850 $μ$m based on a sample of four prototypical $z\sim2-3$ enormous Ly$α$ nebulae (ELANe) as well as their megaparsec-scale (Mpc-scale) environments to study the physical connections between ELANe and their coeval dusty submillimeter galaxies (SMGs). By analysing the SCUBA-2 data with self-consistent Monte Carlo simulations to construct the number counts, here, we report on the overabundance of 850 $μ$m-selected submillimeter sources around all the four ELANe, by a factor of 3.6$\pm$0.6 (weighted average) compared to the blank fields. This suggests that the excessive number of submillimeter sources are likely to be part of the Mpc-scale environment around the ELANe, corroborating the co-evolution scenario for SMGs and quasars; this is a process which may be more commonly observed in the ELAN fields. If the current form of the underlying count models continues toward the fainter end, our results would suggest an excess of the 850 $μ$m extragalactic background light by a factor of between 2-10, an indication of significant background light fluctuations on the survey scales. Finally, by assuming that all the excessive submillimeter sources are associated with their corresponding ELAN environments, we estimate the SFR densities of each ELAN field, as well as a weighted average of $Σ$SFR=1200$\pm$300 $M_\odot$ yr$^{-1}$ Mpc$^{-3}$, consistent with that found in the vicinity of other quasar systems or proto-clusters at similar redshifts; in addition, it is a factor of about 300 greater than the cosmic mean.

preprint2021arXiv

A Multiwavelength Study of ELAN Environments (AMUSE$^2$). Mass budget, satellites spin alignment and gas infall in a massive $z\sim3$ quasar host halo

The systematic targeting of extended Ly$α$ emission around high-redshift quasars resulted in the discovery of rare and bright Enormous Ly$α$ Nebulae (ELANe) associated with multiple active galactic nuclei (AGN). We here initiate &#34;a multiwavelength study of ELAN environments&#34; (AMUSE$^2$) focusing on the ELAN around the $z\sim3$ quasar SDSS J1040+1020, a.k.a. the Fabulous ELAN. We report on VLT/HAWK-I, APEX/LABOCA, JCMT/SCUBA-2, SMA/850$μ$m, ALMA/CO(5-4) and 2mm observations and compare them to previously published VLT/MUSE data. The continuum and line detections enable a first estimate of the star-formation rates, dust, stellar and molecular gas masses in four objects associated with the ELAN (three AGNs and one Ly$α$ emitter), confirming that the quasar host is the most star-forming (${\rm SFR}\sim500$ M$_\odot$ yr$^{-1}$) and massive galaxy ($M_{\rm star}\sim10^{11}$ M$_{\odot}$) in the system, and thus can be assumed as central. All four embedded objects have similar molecular gas reservoirs ($M_{\rm H_2}\sim10^{10}$ M$_{\odot}$), resulting in short depletion time scales. This fact together with the estimated total dark-matter halo mass, $M_{\rm DM}=(0.8-2)\times10^{13}$ M$_{\odot}$, implies that this ELAN will evolve into a giant elliptical galaxy. Consistently, the constraint on the baryonic mass budget for the whole system indicates that the majority of baryons should reside in a massive warm-hot reservoir (up to $10^{12}$ M$_{\odot}$), needed to complete the baryons count. Additionally, we discuss signatures of gas infall on the compact objects as traced by Ly$α$ radiative transfer effects and the evidence for the alignment between the satellites&#39; spins and their directions to the central.

preprint2021arXiv

Characterizing the FRB host galaxy population and its connection to transients in the local and extragalactic Universe

We present the localization and host galaxies of one repeating and two apparently non-repeating Fast Radio Bursts. FRB20180301A was detected and localized with the Karl G. Jansky Very Large Array to a star-forming galaxy at $z=0.3304$. FRB20191228A, and FRB20200906A were detected and localized by the Australian Square Kilometre Array Pathfinder to host galaxies at $z=0.2430$ and $z=0.3688$, respectively. We combine these with 13 other well-localized FRBs in the literature, and analyze the host galaxy properties. We find no significant differences in the host properties of repeating and apparently non-repeating FRBs. FRB hosts are moderately star-forming, with masses slightly offset from the star-forming main-sequence. Star formation and low-ionization nuclear emission-line region (LINER) emission are major sources of ionization in FRB host galaxies, with the former dominant in repeating FRB hosts. FRB hosts do not track stellar mass and star formation as seen in field galaxies (more than 95% confidence). FRBs are rare in massive red galaxies, suggesting that progenitor formation channels are not solely dominated by delayed channels which lag star formation by Gigayears. The global properties of FRB hosts are indistinguishable from core-collapse supernovae (CCSNe) and short gamma-ray bursts (SGRBs) hosts, and the spatial offset (from galaxy centers) of FRBs is mostly inconsistent with that of the Galactic neutron star population (95% confidence). The spatial offsets of FRBs (normalized to the galaxy effective radius) also differ from those of globular clusters (GCs) in late- and early-type galaxies with 95% confidence.

preprint2021arXiv

On the Environments of Giant Radio Galaxies

We test the hypothesis that environments play a key role in enabling the growth of enormous radio structures spanning more than 700 kpc, an extreme population of radio galaxies called giant radio galaxies (GRGs). To achieve this, we explore (1) the relationships between the occurrence of GRGs and the surface number density of surrounding galaxies, including satellite galaxies and galaxies from neighboring halos, as well as (2) the GRG locations towards large-scale structures. The analysis is done by making use of a homogeneous sample of 110 GRGs detected from the LOFAR Two-metre Sky Survey in combination with photometric galaxies from the DESI Legacy Imaging Surveys and a large-scale filament catalog from the Sloan Digital Sky Survey. Our results show that the properties of galaxies around GRGs are similar with that around the two control samples, consisting of galaxies with optical colors and luminosity matched to the properties of the GRG host galaxies. Additionally, the properties of surrounding galaxies depend on neither their relative positions to the radio jet/lobe structures nor the sizes of GRGs. We also find that the locations of GRGs and the control samples with respect to the nearby large-scale structures are consistent with each other. These results demonstrate that there is no correlation between the GRG properties and their environments traced by stars, indicating that external galaxy environments are not the primary cause of the large sizes of the radio structures. Finally, regarding radio feedback, we show that the fraction of blue satellites does not correlate with the GRG properties, suggesting that the current epoch of radio jets have minimal influence on the nature of their surrounding galaxies.

preprint2020arXiv

3D Distribution Map of HI Gas and Galaxies Around an Enormous Ly$α$ Nebula and Three QSOs at $z=2.3$ Revealed by the HI Tomographic Mapping Technique

We present an IGM HI tomography map in a survey volume of $16 \times 19 \times 131 \ h^{-3} {\rm comoving \ Mpc}^{3}$ (cMpc$^3$) centered at MAMMOTH-1 nebula and three neighbouring quasars at $z=2.3$. MAMMOTH-1 nebula is an enormous Ly$α$ nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over $1\ h^{-1}$ cMpc with not fully clear physical origin. Here we investigate the HI-gas distribution around MAMMOTH1-QSO with the ELAN and three neighbouring type-I quasars, making the IGM HI tomography map with a spatial resolution of $2.6\ h^{-1}$ cMpc. Our HI tomography map is reconstructed with HI Ly$α$ forest absorption of bright background objects at $z=2.4-2.9$: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of HI flux overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with significantly weak HI absorption. This suggests that MAMMOTH1-QSO has a proximity zone where quasar illuminates and photo-ionizes the surrounding HI gas and suppresses HI absorption, and that the ELAN is probably a photo-ionized cloud embedded in the cosmic web. The HI radial profile of MAMMOTH1-QSO is very similar to those of three neighbouring type-I quasars at $z=2.3$, which is compatible with the AGN unification model. We compare the distributions of the HI absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and HI gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.

preprint2020arXiv

A Cold, Massive, Rotating Disk Galaxy 1.5 Billion Years after the Big Bang

Massive disk galaxies like the Milky Way are expected to form at late times in traditional models of galaxy formation, but recent numerical simulations suggest that such galaxies could form as early as a billion years after the Big Bang through the accretion of cold material and mergers. Observationally, it has been difficult to identify disk galaxies in emission at high redshift, in order to discern between competing models of galaxy formation. Here we report imaging, with a resolution of about 1.3 kiloparsecs, of the 158-micrometre emission line from singly ionized carbon, the far-infrared dust continuum and the near-ultraviolet continuum emission from a galaxy at a redshift of 4.2603, identified by detecting its absorption of quasar light. These observations show that the emission arises from gas inside a cold, dusty, rotating disk with a rotational velocity of 272 kilometres per second. The detection of emission from carbon monoxide in the galaxy yields a molecular mass that is consistent with the estimate from the ionized carbon emission of about 72 billion solar masses. The existence of such a massive, rotationally supported, cold disk galaxy when the Universe was only 1.5 billion years old favours formation through either cold-mode accretion or mergers, although its large rotational velocity and large content of cold gas remain challenging to reproduce with most numerical simulations.

preprint2020arXiv

A Data-Driven Technique Using Millisecond Transients to Measure the Milky Way Halo

We introduce a new technique to constrain the line-of-sight integrated electron density of our Galactic halo $\text{DM}_\text{MW,halo}$ through analysis of the observed dispersion measure distributions of pulsars $\text{DM}_\text{pulsar}$ and fast radio bursts $\text{DM}_\text{FRB}$. We model these distributions, correcting for the Galactic interstellar medium, with kernel density estimation---well-suited to the small data regime---to find lower/upper bounds to the corrected $\text{DM}_\text{pulsar}$/$\text{DM}_\text{FRB}$ distributions: $\max[\text{DM}_\text{pulsar}] \approx 7\pm2 \text{ (stat)} \pm 9 \text{ (sys) pc cm}^{-3}$ and $\min[\text{DM}_\text{FRB}] \approx 63^{+27}_{-21} \text{ (stat)} \pm 9 \text{ (sys) pc cm}^{-3}$. Using bootstrap resampling to estimate uncertainties, we set conservative limits on the Galactic halo dispersion measure $-2 < \text{DM}_\text{MW,halo} < 123 \text{pc cm}^{-3}$ (95\% c.l.). The upper limit is especially conservative because it may include a non-negligible contribution from the FRB host galaxies and a non-zero contribution from the cosmic web. It strongly disfavors models where the Galaxy has retained the majority of its baryons with a density profile tracking the presumed dark matter density profile. Last, we perform Monte Carlo simulations of larger FRB samples to validate our technique and assess the sensitivity of ongoing and future surveys. We recover bounds of several tens $\text{pc cm}^{-3}$ which may be sufficient to test whether the Galaxy has retained a majority of its baryonic mass. We estimate that a sample of several thousand FRBs will significantly tighten constraints on $\text{DM}_\text{MW,halo}$ and offer a valuable complement to other analyses.

preprint2020arXiv

A metal-poor damped Ly-alpha system at redshift 6.4

We identify a strong Ly-alpha damping wing profile in the spectrum of the quasar P183+05 at z=6.4386. Given the detection of several narrow metal absorption lines at z=6.40392, the most likely explanation for the absorption profile is that it is due to a damped Ly-alpha system. However, in order to match the data a contribution of an intergalactic medium 5-38% neutral or additional weaker absorbers near the quasar is also required. The absorption system presented here is the most distant damped Ly-alpha system currently known. We estimate an HI column density ($10^{20.68\pm0.25}\,$cm$^{-2}$), metallicity ([O/H]$=-2.92\pm 0.32$), and relative chemical abundances of a system consistent with a low-mass galaxy during the first Gyr of the universe. This object is among the most metal-poor damped Ly-alpha systems known and, even though it is observed only ~850 Myr after the big bang, its relative abundances do not show signatures of chemical enrichment by Population III stars.

preprint2020arXiv

A search for supernova-like optical counterparts to ASKAP-localised Fast Radio Bursts

Fast radio bursts (FRBs) are millisecond-scale radio pulses, which originate in distant galaxies and are produced by unknown sources. The mystery remains partially because of the typical difficulty in localising FRBs to host galaxies. Accurate localisations delivered by the Commensal Real-time ASKAP Fast Transients (CRAFT) survey now provide an opportunity to study the host galaxies and potential transient counterparts of FRBs at a large range of wavelengths. In this work, we investigate whether the first three FRBs accurately localised by CRAFT have supernova-like transient counterparts. We obtained two sets of imaging epochs with the Very Large Telescope for three host galaxies, one soon after the burst detection and one several months later. After subtracting these images no optical counterparts were identified in the associated FRB host galaxies, so we instead place limits on the brightness of any potential optical transients. A Monte Carlo approach, in which supernova light curves were modelled and their base properties randomised, was used to estimate the probability of a supernova associated with each FRB going undetected. We conclude that Type Ia and IIn supernovae are unlikely to accompany every apparently non-repeating FRB.

preprint2020arXiv

ALMA [CII] 158um imaging of an HI-Selected Major Merger at z~4

We present high spatial-resolution (~2kpc) Atacama Large Millimeter/submillimeter Array (ALMA) observations of [CII] 158um and dust-continuum emission from a galaxy at z=3.7978 selected by its strong HI absorption (a damped Lya absorber, DLA) against a background QSO. Our ALMA images reveal a pair of star-forming galaxies separated by ~6kpc (projected) undergoing a major merger. Between these galaxies is a third emission component with highly elevated (2x) [CII] 158um emission relative to the dust continuum, which is likely to arise from stripped gas associated with the merger. This merger of two otherwise-normal galaxies is not accompanied by enhanced star-formation, contrary to mergers detected in most luminosity-selected samples. The DLA associated with the merger exhibits extreme kinematics, with a velocity width for the low-ionization metal lines of Dv90~470km/s that spans the velocity spread revealed in the [CII] 158um emission. We propose that DLAs with high Dv90 values are a signpost of major mergers in normal galaxies at high redshifts, and use the distribution of the velocity widths of metal lines in high-z DLAs to provide a rough estimate the fraction of z>3 galaxies that are undergoing a major merger.

preprint2020arXiv

Confronting the Magnetar Interpretation of Fast Radio Bursts Through Their Host Galaxy Demographics

We explore the prompt magnetar progenitor scenario in the context of fast radio burst (FRB) host galaxies demographics and offset distributions. Magnetars are neutron stars with strong magnetic fields on the order of $10^{15}$ G with a short decay lifetime of less than $10^4$ years. Due to their extremely short lifetimes, magnetars should follow the demographics of galaxies according to their current star-formation rate (SFR). Moreover, we hypothesize that magnetars should follow the SFR profile within galaxies, which we assume to follow an exponential profile. We construct a simple model for the host galaxies of magnetars assuming these events track SFR in all galaxies and compare it to observed properties from a sample of \nsecure\ secure FRB hosts. We find the distribution of observed SFRs is inconsistent with the model at $>95\%$ c.l. The offset distribution is consistent with this scenario; however, this could be due to the limited sample size and the seeing limited estimates for the effective radii of the FRB host galaxies. Despite the recent association of an FRB with a magnetar in the Milky Way, magnetars may not be the only source of FRBs in the universe, yet any other successful model must account for the demographics of the FRB host in SFR and their observed galactocentric offsets.

preprint2020arXiv

Constraining Magnetic Fields in the Circumgalactic Medium

We study the properties of magnetic fields in the circumgalactic medium (CGM) of $z<1$ galaxies by correlating Faraday rotation measures of $\sim 1,000$ high-redshift radio sources with the foreground galaxy number density estimated from the DESI Legacy Imaging Surveys. This method enables us to extract signals of rotation measures contributed by intervening gas around multiple galaxies. Our results show that there is no detectable correlation between the distribution of rotation measures and the number of foreground galaxies, contrary to several previous results. Utilizing the non-detection signals, we estimate $3σ$ upper limits to the rotation measures from the CGM of $\sim 20 \rm \ rad/m^{2}$ within~50 kpc and $\sim 10 \rm \ rad / m^{2}$ at separations of $100$ kpc. By adopting a column density distribution of ionized gas obtained from absorption line measurements, we further estimate the strengths of coherent magnetic fields parallel to the line of sight of $<\rm 2 \ μG$ in the CGM. We show that the estimated upper limits of rotation measures and magnetic field strengths are sufficient to constrain outputs of recent galaxy magneto-hydrodynamic simulations. Finally, we discuss possible causes for the inconsistency between our results and previous works.

preprint2020arXiv

Deep Hubble Space Telescope Imaging on the Extended Ly$α$ Emission of a QSO at $z = 2.19$ with Damped Lyman Alpha System as a Natural Coronagraph

Recent observations suggest that proximate damped Ly$α$ (PDLA) systems can be used to study the host galaxies of Quasi-stellar objects (QSOs), because the PDLAs can block the bright point-spread-function (PSF) from central QSOs. Using six-orbits of narrowband imaging with $\it{HST}$/WFC3, we present the first high resolution narrowband image of the Ly$α$ emission in the PDLA trough of the QSO SDSSJ115432.67-021537.9 . We detect one major component and one minor component in the narrowband imaging. Combining the $\it{HST}$/WFC3 imaging with deep Magellan/MagE spectra, we measure that the Ly$α$ flux F$_{\rm{Lyα}}$ = 1.56$\pm0.10 \times10^{-16}$ $\rm{erg\ s^{-1}\ cm^{-2}}$, which is among the luminous ($\approx$ 2.7 L$^{\star}_{\rm{Lyα}}$) Ly$α$ emitters at $\it{z} =$ 2.19. The $\rm{Lyα}$-based star formation rate (SFR) is $\sim$ 7 $\textrm{M_{\sun} \ yr^{-1}}$. These observational results favor that the star formation from the host galaxy could be the main mechanism to power the Ly$α$ emission. This new method sheds new light on the study of the kinematic structure and the spatial distribution of the extended Ly$α$ emitting regions around the QSO host.

preprint2020arXiv

Effective Opacity of the Intergalactic Medium from Galaxy Spectra Analysis

We measure the effective opacity ($τ_{eff}$) of the Intergalactic Medium (IGM) from the composite spectra of 281 Lyman-Break Galaxies (LBGs) in the redshift range $2 \lesssim z \lesssim 3$. Our spectra are taken from the COSMOS Lyman-Alpha Mapping And Tomographic Observations (CLAMATO) survey derived from the Low Resolution Imaging Spectrometer (LRIS) on the W.M. Keck I telescope. We generate composite spectra in two redshift intervals and fit them with spectral energy distribution (SED) models composed of simple stellar populations. Extrapolating these SED models into the Ly$α$ forest, we measure the effective Ly$α$ opacity ($τ_{eff}$) in the $2.02 \leq z \leq 2.44$ range. At $z = 2.22$, we estimate $τ_{eff} = 0.159 \pm 0.001$ from a power-law fit to the data. These measurements are consistent with estimates from quasar analyses at $z<2.5$ indicating that the systematic errors associated with normalizing quasar continua are not substantial. We provide a Gaussian Processes model of our results and previous $τ_{eff}$ measurements that describes the steep redshift evolution in $τ_{eff}$ from $z = 1.5 - 4$.

preprint2020arXiv

First Constraints on Compact Dark Matter from Fast Radio Burst Microstructure

Despite existing constraints, it remains possible that up to $35\%$ of all dark matter is comprised of compact objects, such as the black holes in the 10-100\,M$_\odot$ range whose existence has been confirmed by LIGO. The strong gravitational lensing of transients such as FRBs and GRBs has been suggested as a more sensitive probe for compact dark matter than intensity fluctuations observed in microlensing experiments. Recently ASKAP has reported burst substructure down to $15μ$s timescales in FRBs in the redshift range $0.3-0.5$. We investigate here the implications of this for the detectability of compact dark matter by FRBs. We find that a sample size of $\sim130$ FRBs would be required to constrain compact dark matter to less than the existing 35$\%$ limit with 95$\%$ confidence, if it were distributed along $\gtrsim 1\,$Gpc-long FRB sightlines through the cosmic web. Conversely, existing constraints on the fraction of compact dark matter permit as many as 1 in $\approx 40$ of all $z \lesssim 0.4$ FRBs to exhibit micro-lensed burst structure. Approximately $170$ FRBs intercepting halos within $\sim 50\,$kpc would be required to exclude the fraction of compact dark matter in each intercepted halo to a similar level. Furthermore, we consider the cumulative effects of lensing of the FRB signal by a macroscopic dark matter distribution. We conclude that lensing from a uniform distribution of compact objects is likely not observable, but suggest that FRBs may set meaningful limits on power-law distributions of dark matter.

preprint2020arXiv

Fully probabilistic quasar continua predictions near Lyman-α with conditional neural spline flows

Measurement of the red damping wing of neutral hydrogen in quasar spectra provides a probe of the epoch of reionization in the early Universe. Such quantification requires precise and unbiased estimates of the intrinsic continua near Lyman-$α$ (Ly$α$), a challenging task given the highly variable Ly$α$ emission profiles of quasars. Here, we introduce a fully probabilistic approach to intrinsic continua prediction. We frame the problem as a conditional density estimation task and explicitly model the distribution over plausible blue-side continua ($1190\ \unicode{xC5} \leq λ_{\text{rest}} < 1290\ \unicode{xC5}$) conditional on the red-side spectrum ($1290\ \unicode{xC5} \leq λ_{\text{rest}} < 2900\ \unicode{xC5}$) using normalizing flows. Our approach achieves state-of-the-art precision and accuracy, allows for sampling one thousand plausible continua in less than a tenth of a second, and can natively provide confidence intervals on the blue-side continua via Monte Carlo sampling. We measure the damping wing effect in two $z>7$ quasars and estimate the volume-averaged neutral fraction of hydrogen from each, finding $\bar{x}_\text{HI}=0.304 \pm 0.042$ for ULAS J1120+0641 ($z=7.09$) and $\bar{x}_\text{HI}=0.384 \pm 0.133$ for ULAS J1342+0928 ($z=7.54$).

preprint2020arXiv

High time resolution and polarisation properties of ASKAP-localised fast radio bursts

Combining high time and frequency resolution full-polarisation spectra of Fast Radio Bursts (FRBs) with knowledge of their host galaxy properties provides an opportunity to study both the emission mechanism generating them and the impact of their propagation through their local environment, host galaxy, and the intergalactic medium. The Australian Square Kilometre Array Pathfinder (ASKAP) telescope has provided the first ensemble of bursts with this information. In this paper, we present the high time and spectral resolution, full polarisation observations of five localised FRBs to complement the results published for the previously studied ASKAP FRB~181112. We find that every FRB is highly polarised, with polarisation fractions ranging from 80 -- 100\%, and that they are generally dominated by linear polarisation. While some FRBs in our sample exhibit properties associated with an emerging archetype (i.e., repeating or apparently non-repeating), others exhibit characteristic features of both, implying the existence of a continuum of FRB properties. When examined at high time resolution, we find that all FRBs in our sample have evidence for multiple sub-components and for scattering at a level greater than expected from the Milky Way. We find no correlation between the diverse range of FRB properties (e.g., scattering time, intrinsic width, and rotation measure) and any global property of their host galaxy. The most heavily scattered bursts reside in the outskirts of their host galaxies, suggesting that the source-local environment rather than the host interstellar medium is likely the dominant origin of the scattering in our sample.

preprint2020arXiv

Modelling the prompt optical emission of GRB 180325A: the evolution of a spike from the optical to gamma-rays

The transition from prompt to the afterglow emission is one of the most exciting and least understood phases in gamma-ray bursts (GRBs). Correlations among optical, X-ray and gamma-ray emission in GRBs have been explored, to attempt to answer whether the earliest optical emission comes from internal and/or external shocks. We present optical photometric observations of GRB 180325A collected with the TAROT and RATIR ground-based telescopes. These observations show two strong optical flashes with separate peaks at $\sim50\;$s and $\sim120\;$s, followed by a temporally extended optical emission. We also present X-rays and gamma-ray observations of GRB 180325A, detected by the Burst Alert Telescope (BAT) and X-ray Telescope (XRT), on the Neil Gehrels Swift observatory, which both observed a narrow flash at $\sim80\;$s. We show that the prompt gamma- and X-ray early emission shares similar temporal and spectral features consistent with internal dissipation within the relativistic outflow (e.g. by internal shocks or magnetic reconnection), while the early optical flashes are likely generated by the reverse shock that decelerates the ejecta as it sweeps up the external medium.

preprint2020arXiv

Multi-phase outflows in post starburst E+A galaxies -- II. A direct connection between the neutral and ionized outflow phases

Post starburst E+A galaxies are systems that hosted a powerful starburst that was quenched abruptly. Simulations suggest that these systems provide the missing link between major merger ULIRGs and red and dead ellipticals, where AGN feedback is responsible for the expulsion or destruction of the molecular gas. However, many details remain unresolved and little is known about AGN-driven winds in this short-lived phase. We present spatially-resolved IFU spectroscopy with MUSE/VLT of SDSS J124754.95-033738.6, a post starburst E+A galaxy with a recent starburst that started 70 Myrs ago and ended 30 Myrs ago, with a peak SFR of $\sim 150\,\mathrm{M_{\odot}/yr}$. We detect disturbed gas throughout the entire field of view, suggesting triggering by a minor merger. We detect fast-moving multi-phased gas clouds, embedded in a double-cone face-on outflow, which are traced by ionized emission lines and neutral NaID \emph{emission} and absorption lines. We find remarkable similarities between the kinematics, spatial extents, and line luminosities of the ionized and neutral gas phases, and propose a model in which they are part of the same outflowing clouds, which are exposed to both stellar and AGN radiation. Our photoionization model provides consistent ionized line ratios, NaID absorption optical depths and EWs, and dust reddening. Using the model, we estimate, for the first time, the neutral-to-ionized gas mass ratio (about 20), the sodium neutral fraction, and the size of the outflowing clouds. This is one of the best ever observed direct connections between the neutral and ionized outflow phases in AGN.

preprint2020arXiv

Polyphorm: Structural Analysis of Cosmological Datasets via Interactive Physarum Polycephalum Visualization

This paper introduces Polyphorm, an interactive visualization and model fitting tool that provides a novel approach for investigating cosmological datasets. Through a fast computational simulation method inspired by the behavior of Physarum polycephalum, an unicellular slime mold organism that efficiently forages for nutrients, astrophysicists are able to extrapolate from sparse datasets, such as galaxy maps archived in the Sloan Digital Sky Survey, and then use these extrapolations to inform analyses of a wide range of other data, such as spectroscopic observations captured by the Hubble Space Telescope. Researchers can interactively update the simulation by adjusting model parameters, and then investigate the resulting visual output to form hypotheses about the data. We describe details of Polyphorm&#39;s simulation model and its interaction and visualization modalities, and we evaluate Polyphorm through three scientific use cases that demonstrate the effectiveness of our approach.

preprint2020arXiv

Project AMIGA: The Circumgalactic Medium of Andromeda

Project AMIGA (Absorption Maps In the Gas of Andromeda) is a large ultraviolet Hubble Space Telescope program, which has assembled a sample of 43 QSOs that pierce the circumgalactic medium (CGM) of Andromeda (M31) from R=25 to 569 kpc (25 of them probing gas from 25 kpc to about the virial radius-Rvir = 300 kpc-of M31). Our large sample provides an unparalleled look at the physical conditions and distribution of metals in the CGM of a single galaxy using ions that probe a wide range of gas phases (Si II, Si III, Si IV, C II, C IV, and O VI, the latter being from the Far Ultraviolet Spectroscopic Explorer). We find that Si III and O VI have near unity covering factor maintained all the way out to 1.2Rvir and 1.9Rvir, respectively. We show that Si III is the dominant ion over Si II and Si IV at any R. While we do not find that the properties of the CGM of M31 depend strongly on the azimuth, we show that they change remarkably around 0.3-0.5Rvir, conveying that the inner regions of the CGM of M31 are more dynamic and have more complicated multi-phase gas-structures than at R>0.5Rvir. We estimate the metal mass of the CGM within Rvir as probed by Si II, Si III, and Si IV is 2x10^7 Msun and by O VI is >8x10^7 Msun, while the baryon mass of the 10^4-10^5.5 K gas is ~4x10^10 (Z/0.3 Zsun)^(-1) Msun within Rvir. We show that different zoom-in cosmological simulations of L* galaxies better reproduce the column density profile of O VI with R than Si III or the other studied ions. We find that observations of the M31 CGM and zoom-in simulations of L* galaxies have both lower ions showing higher column density dispersion and dependence on R than higher ions, indicating that the higher ionization structures are larger and/or more broadly distributed.

preprint2020arXiv

PypeIt: The Python Spectroscopic Data Reduction Pipeline

PypeIt is a Python package for semi-automated reduction of astronomical, spectroscopic data. Its algorithms build on decades-long development of previous data reduction pipelines by the developers (Bernstein, Burles, & Prochaska, 2015; Bochanski et al., 2009). The reduction procedure -- including a complete list of the input parameters and available functionality -- is provided as online documentation hosted by Read the Docs, which is regularly updated. (https://pypeit.readthedocs.io/en/latest/). Release v1.0.3 serves the following spectrographs: Gemini/GNIRS, Gemini/GMOS, Gemini/FLAMINGOS 2, Lick/Kast, Magellan/MagE, Magellan/Fire, MDM/OSMOS, Keck/DEIMOS (600ZD, 830G, 1200G), Keck/LRIS, Keck/MOSFIRE (J and Y gratings tested), Keck/NIRES, Keck/NIRSPEC (low-dispersion), LBT/Luci-I, Luci-II, LBT/MODS (beta), NOT/ALFOSC (grism4), VLT/X-Shooter (VIS, NIR), VLT/FORS2 (300I, 300V), WHT/ISIS.

preprint2020arXiv

Revealing the Dark Threads of the Cosmic Web

Modern cosmology predicts that matter in our Universe has assembled today into a vast network of filamentary structures colloquially termed the Cosmic Web. Because this matter is either electromagnetically invisible (i.e., dark) or too diffuse to image in emission, tests of this cosmic web paradigm are limited. Wide-field surveys do reveal web-like structures in the galaxy distribution, but these luminous galaxies represent less than 10% of baryonic matter. Statistics of absorption by the intergalactic medium (IGM) via spectroscopy of distant quasars support the model yet have not conclusively tied the diffuse IGM to the web. Here, we report on a new method inspired by the Physarum polycephalum slime mold that is able to infer the density field of the Cosmic Web from galaxy surveys. Applying our technique to galaxy and absorption-line surveys of the local Universe, we demonstrate that the bulk of the IGM indeed resides in the Cosmic Web. From the outskirts of Cosmic Web filaments, at approximately the cosmic mean matter density (rho_m) and approx. 5 virial radii from nearby galaxies, we detect an increasing H I absorption signature towards higher densities and the circumgalactic medium, to approx. 200 rho_m. However, the absorption is suppressed within the densest environments, suggesting shock-heating and ionization deep within filaments and/or feedback processes within galaxies.

preprint2020arXiv

Spectropolarimetric analysis of FRB 181112 at microsecond resolution: Implications for Fast Radio Burst emission mechanism

We have developed a new coherent dedispersion mode to study the emission of Fast Radio Bursts that trigger the voltage capture capability of the Australian SKA Pathfinder (ASKAP) interferometer. In principle the mode can probe emission timescales down to 3 ns with full polarimetric information preserved. Enabled by the new capability, here we present a spectropolarimetric analysis of FRB 181112 detected by ASKAP, localized to a galaxy at redshift 0.47. At microsecond time resolution the burst is resolved into four narrow pulses with a rise time of just $15 μ$s for the brightest. The pulses have a diversity of morphology, but do not show evidence for temporal broadening by turbulent plasma along the line of sight, nor is there any evidence for periodicity in their arrival times. The pulses are highly polarized (up to 95%), with the polarization position angle varying both between and within pulses. The pulses have apparent rotation measures that vary by $15\pm 2\, {\rm rad \,m^{-2}}$ and apparent dispersion measures that vary by $0.041\pm 0.004\,{\rm pc\,cm^{-3}}$. Conversion between linear and circular polarization is observed across the brightest pulse. We conclude that the FRB 181112 pulses are most consistent with being a direct manifestation of the emission process or the result of propagation through a relativistic plasma close to the source. This demonstrates that our method, which facilitates high-time-resolution polarimetric observations of FRBs, can be used to study not only burst emission processes, but also a diversity of propagation effects present on the gigaparsec paths they traverse.

preprint2020arXiv

Statistical correlation between the distribution of Ly$α$ emitters and IGM HI at $z\sim2.2$ mapped by Subaru/Hyper Suprime-Cam

The correlation between neutral Hydrogen (HI) in the intergalactic medium (IGM) and galaxies now attracts great interests. We select four fields which include several coherently strong Ly$α$ absorption systems at $z\sim2.2$ detected by using background quasars from the whole SDSS/(e)BOSS database. Deep narrow-band and $g$-band imaging are performed using the Hyper Suprime-Cam on the Subaru Telescope. We select out 2,642 Ly$α$ emitter (LAE) candidates at $z=2.177\pm0.023$ down to the Ly$α$ luminosity of $L_{\text{Ly}α}\approx 2 \times 10^{42} {\rm erg~s}^{-1}$ to construct the galaxy overdensity maps, covering an effective area of 5.39 deg$^2$. Combining the sample with the Ly$α$ absorption estimated from 64 (e)BOSS quasar spectra, we find a moderate to strong correlation between the LAE overdensity $δ_{\rm LAE}$ and the effective optical depth $τ_{\rm LoS}$ in line-of-sights, with $P$-value$=0.09\%$ ($<0.01\%$) when the field that contains a significant quasar overdensity is in(ex)cluded. The cross-correlation analysis also clearly suggests that up to $4\pm1$ pMpc, LAEs tend to cluster in the regions rich in HI gas, indicated by the high $τ_{\rm LoS}$, and avoid the low $τ_{\rm LoS}$ region where the HI gas is deficient. By averaging the $τ_{\rm LoS}$ as a function of the projected distance ($d$) to LAEs, we find a $30\%$ excess signal at $2σ$ level at $d<200$ pkpc, indicating the dense HI in circumgalactic medium, and a tentative excess at $400<d<600$ pkpc in IGM regime, corroborating the cross-correlation signal detected at about $0.5$ pMpc. These statistical analyses indicate that galaxy$-$IGM HI correlations exist on scales ranging from several hundred pkpc to several pMpc at $z\sim2.2$.

preprint2020arXiv

Testing Galaxy Formation Simulations with Damped Lyman-$α$ Abundance and Metallicity Evolution

We examine the properties of damped Lyman-$α$ absorbers (DLAs) emerging from a single set of cosmological initial conditions in two state-of-the-art cosmological hydrodynamic simulations: {\sc Simba} and {\sc Technicolor Dawn}. The former includes star formation and black hole feedback treatments that yield a good match with low-redshift galaxy properties, while the latter uses multi-frequency radiative transfer to model an inhomogeneous ultraviolet background (UVB) self-consistently and is calibrated to match the Thomson scattering optical depth, UVB amplitude, and Ly-$α$ forest mean transmission at $z>5$. Both simulations are in reasonable agreement with the measured stellar mass and star formation rate functions at $z\geq 3$, and both reproduce the observed neutral hydrogen cosmological mass density, $Ω_{\rm HI}(z)$. However, the DLA abundance and metallicity distribution are sensitive to the galactic outflows&#39; feedback and the UVB amplitude. Adopting a strong UVB and/or slow outflows under-produces the observed DLA abundance, but yields broad agreement with the observed DLA metallicity distribution. By contrast, faster outflows eject metals to larger distances, yielding more metal-rich DLAs whose observational selection may be more sensitive to dust bias. The DLA metallicity distribution in models adopting an ${\rm H}_2$-regulated star formation recipe includes a tail extending to $[M/H] \ll -3$, lower than any DLA observed to date, owing to curtailed star formation in low-metallicity galaxies. Our results show that DLA observations play an imporant role in constraining key physical ingredients in galaxy formation models, complementing traditional ensemble statistics such as the stellar mass and star formation rate functions.

preprint2020arXiv

The host galaxies and progenitors of Fast Radio Bursts localized with the Australian Square Kilometre Array Pathfinder

The Australian SKA Pathfinder (ASKAP) telescope has started to localize Fast Radio Bursts (FRBs) to arcsecond accuracy from the detection of a single pulse, allowing their host galaxies to be reliably identified. We discuss the global properties of the host galaxies of the first four FRBs localized by ASKAP, which lie in the redshift range $0.11<z<0.48$. All four are massive galaxies (log( $M_{*}/ M_{\odot}$) $\sim 9.4 -10.4$) with modest star-formation rates of up to $2M_{\odot}$yr$^{-1}$ -- very different to the host galaxy of the first repeating FRB 121102, which is a dwarf galaxy with a high specific star-formation rate. The FRBs localized by ASKAP typically lie in the outskirts of their host galaxies, which appears to rule out FRB progenitor models that invoke active galactic nuclei (AGN) or free-floating cosmic strings. The stellar population seen in these host galaxies also disfavors models in which all FRBs arise from young magnetars produced by superluminous supernovae (SLSNe), as proposed for the progenitor of FRB 121102. A range of other progenitor models (including compact-object mergers and magnetars arising from normal core-collapse supernovae) remain plausible.

preprint2020arXiv

The PHLEK Survey: A New Determination of the Primordial Helium Abundance

We present Keck NIRSPEC and Keck NIRES spectroscopy of sixteen metal-poor galaxies that have pre-existing optical observations. The near-infrared (NIR) spectroscopy specifically targets the HeI 10830A emission line, due to its sensitivity to the physical conditions of the gas in HII regions. We use these NIR observations, combined with optical spectroscopy, to determine the helium abundance of sixteen galaxies across a metallicity range 12 + log (O/H) = 7.13 - 8.00. This data set is combined with two other samples where metallicity and helium abundance measurements can be secured: star-forming galaxies selected from the Sloan Digital Sky Survey spectroscopic database and existing low-metallicity systems in the literature. We calculate a linear fit to these measurements, accounting for intrinsic scatter, and report a new determination of the primordial helium number abundance, y_P = 0.0805+/-0.0017, which corresponds to a primordial helium mass fraction Y_P = 0.2436+/-0.0040. Using our determination of the primordial helium abundance in combination with the latest primordial deuterium measurement, (D/H)_P x 10^5 =2.527+/-0.030, we place a bound on the baryon density Omega_bh^2 = 0.0215+/-0.0005 and the effective number of neutrino species Neff = 2.85+/-0.28. These values are in 1.3-sigma agreement with those deduced from the Planck satellite observations of the temperature fluctuations imprinted on the Cosmic Microwave Background.

preprint2020arXiv

The Third Data Release of the KODIAQ Survey

We present and make publicly available the third data release (DR3) of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ) survey. KODIAQ DR3 consists of a fully-reduced sample of 727 quasars at 0.1 < z < 6.4 observed with ESI at moderate resolution (4000 < R < 10000). DR3 contains 872 spectra available in flux calibrated form, representing a sum total exposure time of approximately 2.8 megaseconds. These co-added spectra arise from a total of 2753 individual exposures of quasars taken from the Keck Observatory Archive (KOA) in raw form and uniformly processed using a data reduction package made available through the XIDL distribution. DR3 is publicly available to the community, housed as a higher level science product at the KOA and in the igmspec database.

preprint2019arXiv

Discovery of a Rare Late-Type, Low-Mass Wolf-Rayet Star in the LMC

We report the serendipitous discovery of an object, UVQS J060819.93-715737.4, with a spectrum dominated by extremely intense, narrow C II emission lines. The spectrum is similar to those of the very rare, late-type [WC11] low-mass Wolf-Rayet stars. Despite the recognition of these stars as a distinct class decades ago, there remains barely a handful of Galactic members, all of which are also planetary-nebula central stars. Although no obvious surrounding nebulosity is present in J0608, [O II], [N II], and [S II] emission suggest the presence of an inconspicuous, low-excitation nebula. There is low-amplitude incoherent photometric variability on timescales of days to years, as well as numerous prominent P Cygni profiles, implying mass loss. There are indications of a binary companion. The star is located on the outskirts of the LMC, and the observed radial velocity (~250 km/s) and proper motion strongly suggest membership. If indeed an LMC member, this is the first extragalactic late [WC] star, and the first with an accurately determined luminosity, as the Galactic examples are too distant for precise parallax determinations. A high-quality, broad-coverage spectrum of the prototype of the late [WC] class, CPD -56 8032, is also presented. We discuss different excitation mechanisms capable of producing the great strength of the C II emission. Numerous autoionizing levels of C II are definitely populated by processes other than dielectronic recombination. Despite the spectacular emission spectra, observational selection makes objects such as these difficult to discover. Members of the [WC11] class may in fact be considerably more common than the handful of previously known late [WC] stars.

preprint2019arXiv

Evolution of the Cool Gas in the Circumgalactic Medium (CGM) of Massive Halos -- A Keck Cosmic Web Imager (KCWI) Survey of Ly$α$ Emission around QSOs at $z\approx2$

Motivated by the recent discovery of the near-ubiquity of Ly$α$ emission around $z \gtrsim 3$ QSOs, we performed a systematic study of QSO circumgalactic Ly$α$ emission at $z\approx2$, utilizing the unique capability of the Keck Cosmic Web Imager (KCWI) -- a new wide-field, blue sensitive integral-field spectrograph (IFU). In this paper, we present KCWI observations on a sample of 16 ultraluminous Type-I QSOs at $z=2.1-2.3$ with ionizing luminosities of $L_{\rm{ν_{\rm{LL}}}}=10^{31.1-32.3}$ erg s$^{-1}$ Hz$^{-1}$. We found that 14 out of 16 QSOs are associated with Ly$α$ nebulae with projected linear-sizes larger than 50 physical kpc (pkpc). Among them, four nebulae have enormous Ly$α$ emission with the Ly$α$ surface brightness $SB_{\rm{Lyα}} >10^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ on the $>100$ kpc scale, extending beyond the field of view of KCWI. Our KCWI observations reveal that most $z\approx2$ QSO nebulae have a more irregular morphology compared to those at $z\gtrsim3$. In turn, we measure that the circularly-averaged surface brightness (SB) at $z\approx2$ is 0.4 dex fainter than the redshift-corrected, median SB at $z\gtrsim 3$. The Ly$α$ SB profile (SB$_{\rm{Lyα}}$) of QSOs at $z\approx 2$ can be described by a power law of SB$_{\rm{Lyα},z\approx2.3}= 3.7\times10^{-17}\times(r/40)^{-1.8}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, with the slope similar to that at $z\gtrsim3$. The observed lower redshift-corrected, circularly-averaged SB may be mainly due to the lower covering factor of cool gas clouds in massive halos at $z\approx2$.

preprint2019arXiv

GRB 180620A: Evidence for late-time energy injection

The early optical emission of gamma-ray bursts gives an opportunity to understand the central engine and first stages of these events. About 30\% of GRBs present flares whose origin is still a subject of discussion. We present optical photometry of GRB 180620A with the COATLI telescope and RATIR instrument. COATLI started to observe from the end of prompt emission at $T+39.3$~s and RATIR from $T+121.4$~s. We supplement the optical data with the X-ray light curve from \emph{Swift}/XRT. %The optical and X-ray light curves show very unusual behavior with features clearly beyond the standard fireball model. We observe an optical flare from $T+110$ to $T+550$~s, with a temporal index decay $α_\mathrm{O,decay}=1.32\pm 0.01$, and a $Δt/t=1.63$, which we interpret as the signature of a reverse shock component. After the initial normal decay the light curves show a long plateau from $T+500$ to $T+7800$~s both in X-rays and the optical before decaying again after an achromatic jet break at $T+7800$~s. Fluctuations are seen during the plateau phase in the optical. Adding to the complexity of GRB afterglows, the plateau phase (typically associated with the coasting phase of the jet) is seen in this object after the ``normal&#39;&#39; decay phase (emitted during the deceleration phase of the jet) and the jet break phase occurs directly after the plateau. We suggest that this sequence of events can be explained by a rapid deceleration of the jet with $t_d\lesssim 40$ s due to the high density of the environment ($\approx 100$ cm$^{-3}$) followed by reactivation of the central engine which causes the flare and powers the plateau phase.

preprint2019arXiv

MUSE Analysis of Gas around Galaxies (MAGG) -- I: Survey design and the environment of a near pristine gas cloud at z~3.5

We present the design, methods, and first results of the MUSE Analysis of Gas around Galaxies (MAGG) survey, a large programme on the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope (VLT) which targets 28 z > 3.2 quasars to investigate the connection between optically-thick gas and galaxies at z~3-4. MAGG maps the environment of 52 strong absorption line systems at z > 3, providing the first statistical sample of galaxies associated with gas-rich structures in the early Universe. In this paper, we study the galaxy population around a very metal poor gas cloud at z~3.5 towards the quasar J124957.23-015928.8. We detect three Lyman alpha emitters within <200km/s of the cloud redshift, at projected separations <185 kpc (physical). The presence of star-forming galaxies near a very metal-poor cloud indicates that metal enrichment is still spatially inhomogeneous at this redshift. Based on its very low metallicity and the presence of nearby galaxies, we propose that the most likely scenario for this LLS is that it lies within a filament which may be accreting onto a nearby galaxy. Taken together with the small number of other LLSs studied with MUSE, the observations to date show a range of different environments near strong absorption systems. The full MAGG survey will significantly expand this sample and enable a statistical analysis of the link between gas and galaxies to pin down the origin of these diverse environments at z~3-4.

preprint2019arXiv

The low density and magnetization of a massive galaxy halo exposed by a fast radio burst

Present-day galaxies are surrounded by cool and enriched halo gas extending to hundreds of kiloparsecs. This halo gas is thought to be the dominant reservoir of material available to fuel future star formation, but direct constraints on its mass and physical properties have been difficult to obtain. We report the detection of a fast radio burst (FRB 181112) with arcsecond precision, which passes through the halo of a foreground galaxy. Analysis of the burst shows the halo gas has low net magnetization and turbulence. Our results imply predominantly diffuse gas in massive galactic halos, even those hosting active supermassive black holes, contrary to some previous results.