Source author record

Max Gronke

Max Gronke appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

29 published item(s)

preprint2025arXiv

A weak Ly$α$ halo for an extremely bright Little Red Dot: Indications of enshrouded SMBH growth

The abundant population of "Little Red Dots" (LRDs)-compact objects with red UV to optical colors and broad Balmer lines at high redshift-is unveiling new insights into the properties of early active galactic nuclei (AGN). Perhaps the most surprising features of this population are the presence of Balmer absorption and ubiquitous strong Balmer breaks. Recent models link these features to an active supermassive black hole (SMBH) cocooned in very dense gas ($N_{\rm H}\sim10^{24}\,\rm cm^{-2}$). We present a stringent test of such models using VLT/MUSE observations of A2744-45924, the most luminous LRD known to date ($L_{\rm Hα}\approx10^{44}~\rm erg\,s^{-1}$), located behind the Abell-2744 lensing cluster at $z=4.464$ ($μ=1.8$). We detect a moderately extended Ly$α$ nebula ($h\approx5.7$ pkpc), spatially offset from the point-like H$α$ seen by JWST. The Ly$α$ emission is narrow ($\rm FWHM=270\pm 15~km\,s^{-1}$), spatially offset to H$α$, and faint ($\rm Lyα=0.07Hα$) compared to Ly$α$ nebulae typically observed around quasars of similar luminosity. We detect compact N$\,$IV]$λ$1486 emission, spatially aligned with H$α$, and a spatial shift in the far-UV continuum matching the Ly$α$ offset. We discuss that H$α$ and Ly$α$ have distinct physical origins: H$α$ originates from the AGN, while Ly$α$ is powered by star formation. In the environment of A2744-45924, we identify four extended Ly$α$ halos ($Δz<0.02$, $Δr<100$ pkpc). Their Ly$α$ luminosities match expectations based on H$α$ emission, indicating no evidence for radiation from A2744-45924 affecting its surroundings. The lack of strong, compact, and broad Ly$α$ and the absence of a luminous extended halo, suggest that the UV AGN light is obscured by dense gas cloaking the SMBH with covering factor close to unity.

preprint2025arXiv

Impact of Resonance, Raman, and Thomson Scattering on Hydrogen Line Formation in Little Red Dots

Little Red Dots (LRDs) are compact sources at $z>5$ discovered through JWST spectroscopy. Their spectra exhibit broad Balmer emission lines ($\gtrsim1000\rm~km~s^{-1}$), alongside absorption features and a pronounced Balmer break -- evidence for a dense, neutral hydrogen medium with the $n=2$ state. When interpreted as arising from AGN broad-line regions, inferred black hole masses from local scaling relations exceed expectations given their stellar masses, challenging models of early black hole-galaxy co-evolution. However, radiative transfer effects in dense media may also impact the formation of hydrogen emission lines. We model three scattering processes shaping hydrogen line profiles: resonance scattering by hydrogen in the $n=2$ state, Raman scattering of UV radiation by ground-state hydrogen, and Thomson scattering by free electrons. Using 3D Monte Carlo radiative transfer simulations with multi-branching resonance transitions, we examine their imprint on line shapes and ratios. Resonance scattering produces strong deviations from Case B flux ratios, clear differences between H$α$ and H$β$, and encodes gas kinematics in line profiles but cannot broaden H$β$ due to conversion to Pa$α$. While Raman scattering can yield broad wings, scattering of UV continuum is disfavored given the absence of strong FWHM variations across transitions. Raman scattering of higher Lyman-series emission can produce H$α$/H$β$ wing width ratios of $\gtrsim1.28$, agreeing with observations. Thomson scattering can reproduce the observed $\gtrsim1000~\rm km\, s^{-1}$ wings under plausible conditions, e.g., $T_{\rm e} \sim 10^4\rm \, K$ and $N_{\rm e}\sim10^{24}\rm~cm^{-2}$ -- and lead to black hole mass overestimates by factors $\gtrsim10$. Our results provide a framework for interpreting hydrogen lines in LRDs and similar systems.

preprint2023arXiv

Extended Lyman-$α$ emission towards the SPT2349-56 protocluster at $z=4.3$

Context. Deep spectroscopic surveys with the Atacama Large Millimeter/submillimeter Array (ALMA) have revealed that some of the brightest infrared sources in the sky correspond to concentrations of dusty star-forming galaxies (DSFG) at high redshift. Among these, the SPT2349-56 protocluster system at z = 4.304 is amongst the most extreme examples due to its high source density and integrated star formation rate. Aims. We conducted a deep Lyman-$α$ line emission survey around SPT2349-56 using the Multi-Unit Spectroscopic Explorer (MUSE) at Very Large Telescope (VLT) in order to characterize this uniquely dense environment. Methods. Taking advantage of the deep three-dimensional nature of this survey, we performed a sensitive search for Lyman-$α$ emitters (LAEs) toward the core and northern extension of the protocluster, which correspond to the brightest infrared regions in this field. Using a smoothed narrowband image extracted from the MUSE datacube around the protocluster redshift, we searched for possible extended structures. Results. We identify only three LAEs at z = 4.3 in this field, in concordance with expectations for blank-fields, and an extended Lyman-$α$ structure spatially associated with core of the protocluster. All the previously-identified DSFGs in this field are undetected in Lyman-$α$ emission, consistent with the conspicuous dust obscuration in these systems. We find an extended Lyman-$α$ structure, about $60 \times 60$ kpc$^{2}$ in size, and located 56 kpc west of the protocluster core. Three DSFGs coincide spatially with the location of this structure. We conclude that either the three co-spatial DSFGs or the protocluster core itself are feeding ionizing photons to the Lyman-$α$ structure.

preprint2022arXiv

(Re)Solving Reionization with Lyα: How Bright Lyα Emitters account for the $z\approx2-8$ Cosmic Ionizing Background

The cosmic ionizing emissivity from star-forming galaxies has long been anchored to UV luminosity functions. Here we introduce an emissivity framework based on Ly$α$ emitters (LAEs), which naturally hones in on the subset of galaxies responsible for the ionizing background due to the intimate connection between the production and escape of Ly$α$ and LyC photons. Using constraints on the escape fractions of bright LAEs ($L_{\rm{Lyα}}>0.2 L^{*}$) at $z\approx2$ obtained from resolved Ly$α$ profiles, and arguing for their redshift-invariance, we show that: (i) quasars and LAEs together reproduce the relatively flat emissivity at $z\approx2-6$, which is non-trivial given the strong evolution in both the star-formation density and quasar number density at these epochs and (ii) LAEs produce late and rapid reionization between $z\approx6-9$ under plausible assumptions. Within this framework, the $>10\times$ rise in the UV population-averaged $f_{\rm{esc}}$ between $z\approx3-7$ naturally arises due to the same phenomena that drive the growing Ly$α$ emitter fraction with redshift. Generally, a LAE dominated emissivity yields a peak in the distribution of the ionizing budget with UV luminosity as reported in latest simulations. Using our adopted parameters ($f_{\rm{esc}}=50\%$, $ξ_{\rm{ion}}=10^{25.9}$ Hz erg$^{-1}$ for half the bright LAEs), a highly ionizing minority of galaxies with $M_{\rm UV}<-17$ accounts for the entire ionizing budget from star-forming galaxies. Rapid flashes of LyC from such rare galaxies produce a "disco" ionizing background. We conclude proposing tests to further develop our suggested Ly$α$-anchored formalism.

preprint2022arXiv

Lyα Halos around [O III]-Selected Galaxies in HETDEX

We present extended Lyman-α (Lyα) emission out to 800 kpc of 1034 [O III]-selected galaxies at redshifts 1.9<z<2.35 using the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The locations and redshifts of the galaxies are taken from the 3D-HST survey. The median-stacked surface brightness profile of Lyα emission of the [O III]-selected galaxies agrees well with that of 968 bright Lyα-emitting galaxies (LAEs) at r>40 kpc from the galaxy centers. The surface brightness in the inner parts (r<10 kpc) around the [O III]-selected galaxies, however, is ten times fainter than that of the LAEs. Our results are consistent with the notion that photons dominating the outer regions of the Lyα halos are not produced in the central galaxies but originate outside of them.

preprint2022arXiv

Radiative Turbulent Mixing Layers and the Survival of Magellanic Debris

The Magellanic Stream is sculpted by its infall through the Milky Way's circumgalactic medium, but the rates and directions of mass, momentum, and energy exchange through the Stream-halo interface are relative unknowns critical for determining the origin and fate of the Stream. Complementary to large-scale simulations of LMC-SMC interactions, we apply new insights derived from idealized, high-resolution "cloud-crushing" and radiative turbulent mixing layer simulations to the Leading Arm and Trailing Stream. Contrary to classical expectations of fast cloud breakup, we predict that the Leading Arm and much of the Trailing Stream should be surviving infall and even gaining mass due to strong radiative cooling. Provided a sufficiently supersonic tidal swing-out from the Clouds, the present-day Leading Arm could be a series of high-density clumps in the cooling tail behind the progenitor cloud. We back up our analytic framework with a suite of converged wind-tunnel simulations, finding that previous results on cloud survival and mass growth can be extended to high Mach number ($\mathcal{M}$) flows with a modified drag time $t_{drag} \propto 1 + \mathcal{M}$ and longer growth time. We also simulate the Trailing Stream; we find that the growth time is long ($\sim$ Gyrs) compared to the infall time, and approximate H$α$ emission is low on average ($\sim$few mR) but can be up to tens of mR in bright spots. Our findings also have broader extragalactic implications for e.g. galactic winds, which we discuss.

preprint2022arXiv

Survival and mass growth of cold gas in a turbulent, multiphase medium

Astrophysical gases are commonly multiphase and highly turbulent. In this work, we investigate the survival and growth of cold gas in such a turbulent, multi-phase medium using three-dimensional hydrodynamical simulations. Similar to previous work simulating coherent flow (winds), we find that cold gas survives if the cooling time of the mixed gas is shorter than the Kelvin-Helmholtz time of the cold gas clump (with some weak additional Mach number dependence). However, there are important differences. Near the survival threshold, the long-term evolution is highly stochastic, and subject to the existence of sufficiently large clumps. In a turbulent flow, the cold gas continuously fragments, enhancing its surface area. This leads to exponential mass growth, with a growth time given by the geometric mean of the cooling and the mixing time. The fragmentation process leads to a large number of small droplets which follow a scale-free $\mathrm{d} N/\mathrm{d} m \propto m^{-2}$ mass distribution, and dominate the area covering fraction. Thus, whilst survival depends on the presence of large `clouds', these in turn produce a `fog' of smaller droplets tightly coupled to the hot phase which are probed by absorption line spectroscopy. We show with the aid of Monte-Carlo simulations that the simulated mass distribution emerges naturally due to the proportional mass growth and the coagulation of droplets. We discuss the implications of our results for convergence criteria of larger scale simulations and observations of the circumgalactic medium.

preprint2022arXiv

The Lyman Alpha Spectral Database (LASD)

Lyman $α$ (Ly$α$) emission from star-forming galaxies is an important tool to study a large range of astrophysical questions: it has the potential to carry information about the source galaxy, its nearby circumgalactic medium, and also the surrounding intergalactic medium. Much observational and theoretical work has therefore focused on understanding the details of this emission line. These efforts have been hampered, however, by an absence of spectroscopic reference samples that can be used both as comparisons for observational studies and as critical tests for theoretical work. For this reason we have compiled a large sample of Ly$α$ spectra, at both low and high redshift, and created a publicly available online database, at http://lasd.lyman-alpha.com. The Lyman alpha Spectral Database (LASD) hosts these spectra, as well as large set of spectral and kinematic quantities that have been homogeneously measured for the entire sample. These can easily be viewed online and downloaded in tabular form. The LASD has the capacity for users to easily upload their own Ly$α$ spectra, and all the same spectral measurements will be made, reported, and ingested into the database. We actively invite the community to do so, and the LASD is intended to be a long-term community resource. In this paper we present the design of the database as well as descriptions of the underlying algorithms and the initial Ly$α$ emitter samples that are in the database.

preprint2022arXiv

Tracing Lyman-alpha and Lyman Continuum Escape in Galaxies with Mg II Emission

Star-forming galaxies are considered the likeliest source of the H I ionizing Lyman Continuum (LyC) photons that reionized the intergalactic medium at high redshifts. However, above z >~ 6, the neutral intergalactic medium prevents direct observations of LyC. Therefore, recent years have seen the development of indirect indicators for LyC that can be calibrated at lower redshifts and applied in the Epoch of Reionization. Emission from Mg II \ly\ly 2796, 2803 doublet has been proposed as a promising LyC proxy. In this paper, we present new Hubble Space Telescope/Cosmic Origins Spectrograph observations for 8 LyC emitter candidates, selected to have strong Mg II emission lines. We securely detect LyC emission in 50% (4/8) galaxies with 2$σ$ significance. This high detection rate suggests that strong Mg II emitters might be more likely to leak LyC than similar galaxies without strong Mg II. Using photoionization models, we constrain the escape fraction of Mg II as ~ 15 -- 60%. We confirm that the escape fraction of Mg II correlates tightly with that of Lyman-alpha (LyA), which we interpret as an indication that the escape fraction of both species is controlled by resonant scattering in the same low column density gas. Furthermore, we show that the combination of the Mg II emission and dust attenuation can be used to estimate the escape fraction of LyC statistically. These findings confirm that Mg II emission can be adopted to estimate the escape fraction of LyA and LyC in local star-forming galaxies and may serve as a useful indirect indicator at the Epoch of Reionization.

preprint2022arXiv

Where Outflows Meet Inflows: Gas Kinematics in SSA22 Lyman-$α$ Blob 2 Decoded by Advanced Radiative Transfer Modelling

We present new spectroscopic observations of Lyman-$α$ (Ly$α$) Blob 2 ($z \sim$ 3.1). We observed extended Ly$α$ emission in three distinct regions, where the highest Ly$α$ surface brightness (SB) center is far away from the known continuum sources. We searched through the MOSFIRE slits that cover the high Ly$α$ SB regions, but were unable to detect any significant nebular emission near the highest SB center. We further mapped the flux ratio of the blue peak to the red peak and found it is anti-correlated with Ly$α$ SB with a power-law index of $\sim$ -0.4. We used radiative transfer models with both multiphase, clumpy and shell geometries and successfully reproduced the diverse Ly$α$ morphologies. We found that most spectra suggest outflow-dominated kinematics, while 4/15 spectra imply inflows. A significant correlation exists between parameter pairs, and the multiphase, clumpy model may alleviate previously reported discrepancies. We also modeled Ly$α$ spectra at different positions simultaneously and found that the variation of the inferred clump outflow velocities can be approximately explained by line-of-sight projection effects. Our results support the `central powering + scattering' scenario, i.e. the Ly$α$ photons are generated by a central powering source and then scatter with outflowing, multiphase HI gas while propagating outwards. The infalling of cool gas near the blob outskirts shapes the observed blue-dominated Ly$α$ profiles, but its energy contribution to the total Ly$α$ luminosity is less than 10%, i.e. minor compared to the photo-ionization by star-forming galaxies and/or AGNs.

preprint2021arXiv

Can galaxy evolution mimic cosmic reionization?

Lyman-$α$ (Ly$α$) emitting galaxies are powerful tools to probe the late stages of cosmic reionization. The observed sudden drop in Ly$α$ fraction at $z>6$ is often interpreted as a sign of reionization, since the intergalactic medium (IGM) is more neutral and opaque to Ly$α$ photons. Crucially, this interpretation of the observations is only valid under the assumption that galaxies themselves experience a minimal evolution at these epochs. By modelling Ly$α$ radiative transfer effects in and around galaxies, we examine whether a change in the galactic properties can reproduce the observed drop in the Ly$α$ fraction. We find that an increase in the galactic neutral hydrogen content or a reduction in the outflow velocity toward higher redshift both lead to a lower Ly$α$ escape fraction, and can thus mimic an increasing neutral fraction of the IGM. We furthermore find that this change in galactic properties leads to systematically different Ly$α$ spectra which can be used to differentiate the two competing effects. Using the CANDELSz7 survey measurements which indicate slightly broader lines at $z\sim 6$, we find that the scenario of a mere increase in the galactic column density towards higher $z$ is highly unlikely. We also show that a decrease in outflow velocity is not ruled out by existing data but leads to more prominent blue peaks at $z>6$. Our results caution the use of Ly$α$ observations to estimate the IGM neutral fraction without accounting for the potential change in the galactic properties, e.g., by mapping out the evolution of Ly$α$ spectral characteristics.

preprint2021arXiv

Deciphering the Lyman-$α$ Emission Line: Towards the Understanding of Galactic Properties Extracted from Ly$α$ Spectra via Radiative Transfer Modeling

Existing ubiquitously in the Universe with the highest luminosity, the Lyman-$α$ emission line encodes abundant physical information about the gaseous medium it interacts with. Nevertheless, the resonant nature of Ly$α$ complicates the radiative transfer (RT) modeling of the line profile, making the extraction of physical properties of the surrounding gaseous medium notoriously difficult. In this paper, we revisit the problem of deciphering the Ly$α$ emission line with RT modeling. We reveal intrinsic parameter degeneracies in the widely-used shell model in the optically thick regime for both static and outflowing cases, which suggest the limitations of the model. We have also explored the connection between the more physically realistic multiphase, clumpy model and the shell model. We find that the parameters of a ``very clumpy'' slab model and the shell model have the following correspondences: (1) the total column density of the clumpy slab model is equal to the HI column density of the shell model; (2) the effective temperature of the clumpy slab model, which incorporates the clump velocity dispersion, is equal to the effective temperature of the shell model; (3) the average radial clump outflow velocity is equal to the shell expansion velocity; (4) large intrinsic line widths are required in the shell model to reproduce the wings of the clumpy slab models; (5) adding another phase of hot inter-clump medium will increase peak separation, and the fitted shell expansion velocity lies between the outflow velocities of two phases of gas. Our results provide a viable solution to the major discrepancies associated with Ly$α$ fitting reported in previous literature, and emphasize the importance of utilizing information from additional observations to break the intrinsic degeneracies as well as interpreting the model parameters in a more physically realistic context.

preprint2021arXiv

Revisiting the Gas Kinematics in SSA22 Lyman-$α$ Blob 1 with Radiative Transfer Modeling in a Multiphase, Clumpy Medium

We present new observations of Lyman-$α$ (Ly$α$) Blob 1 (LAB1) in the SSA22 protocluster region ($z$ = 3.09) using the Keck Cosmic Web Imager (KCWI) and the Keck Multi-object Spectrometer for Infrared Exploration (MOSFIRE). By applying matched filtering to the KCWI datacube, we have created a narrow-band Ly$α$ image and identified several prominent features. By comparing the spatial distributions and intensities of Ly$α$ and H$β$, we find that recombination of photo-ionized HI gas followed by resonant scattering is sufficient to explain all the observed Ly$α$/H$β$ ratios. We further decode the spatially-resolved Ly$α$ profiles using both moment maps and Monte-Carlo radiative transfer (MCRT) modeling. By fitting a set of multiphase, 'clumpy' models to the observed Ly$α$ profiles, we are able to reasonably constrain many parameters, namely the HI number density in the inter-clump medium (ICM), the cloud volume filling factor, the random velocity and outflow velocity of the clumps, the HI outflow velocity of the ICM and the local systemic redshift. Our model has successfully reproduced the diverse Ly$α$ morphologies at different locations, and the main results are: (1) The observed Ly$α$ spectra require relatively few clumps per line-of-sight as they have significant fluxes at the line center; (2) The velocity dispersion of the clumps yields a significant broadening of the spectra as observed; (3) The clump bulk outflow can also cause additional broadening if the HI in the ICM is optically thick; (4) The HI in the ICM is responsible for the absorption feature close to the Ly$α$ line center.

preprint2020arXiv

Measuring the properties of reionised bubbles with resolved Lyman alpha spectra

Identifying and characterising reionised bubbles enables us to track both their size distribution, which depends on the primary ionising sources, and the relationship between reionisation and galaxy evolution. We demonstrate that spectrally resolved $z\gtrsim6$ Lyman-alpha (Ly$α$) emission can constrain properties of reionised regions. Specifically, the distant from a source to a neutral region sets the minimum observable Ly$α$ velocity offset from systemic. Detection of flux on the blue side of the Ly$α$ resonance implies the source resides in a large, sufficiently ionised region that photons can escape without significant resonant absorption, and thus constrains both the sizes of and the residual neutral fractions within ionised bubbles. We estimate the extent of the region around galaxies which is optically thin to blue Ly$α$ photons, analogous to quasar proximity zones, as a function of the source's ionising photon output and surrounding gas density. This optically thin region is typically $\lesssim 0.3$ pMpc in radius (allowing transmission of flux $\gtrsim -250$ km s$^{-1}$), $\lesssim 20$% of the distance to the neutral region. In a proof-of-concept, we demonstrate the $z\approx6.6$ galaxy COLA1 -- with a blue Ly$α$ peak -- likely resides in an ionised region $>0.7$ pMpc, with residual neutral fraction $<10^{-5.5}$. To ionise its own proximity zone we infer COLA1 has a high ionising photon escape fraction ($f_{\mathrm{esc}}>0.50$), relatively steep UV slope ($β< -1.79$), and low line-of-sight gas density ($\sim0.5\times$ the cosmic mean), suggesting it is a rare, underdense line-of-sight.

preprint2020arXiv

Spectral shapes of the Lyman-alpha emission from galaxies: I. blueshifted emission and intrinsic invariance with redshift

We demonstrate the redshift-evolution of the spectral profile of H i Lyman-alpha (Lyα) emission from star-forming galaxies. In this first study we pay special attention to the contribution of blueshifted emission. At redshift z = 2.9-6.6, we compile spectra of a sample of 229 Lyα-selected galaxies identified with the Multi-Unit Spectroscopic Explorer at the Very Large Telescope, while at low-z (< 0.44) we use a sample of 74 ultraviolet-selected galaxies observed with the Cosmic Origin Spectrograph onboard the Hubble Space Telescope. At low-z, where absorption from the intergalactic medium (IGM) is negligible, we show that the ratio of Lyα luminosity bluewards and redwards of line center (L_B/R) increases rapidly with increasing equivalent width (W Lyα). This correlation does not, however, emerge at z = 3-4, and we use bootstrap simulations to demonstrate that trends in L B/R should be suppressed by variations in IGM absorption. Our main result is that the observed blueshifted contribution evolves rapidly downwards with increasing redshift: L_B/R~30 % at z~0, but drops to 15 % at z~3, and to below 3 % by z~6. Applying further simulations of the IGM absorption to the unabsorbed COS spectrum, we demonstrate that this decrease in the blue-wing contribution can be entirely attributed to the thickening of intervening Lyα absorbing systems, with no need for additional H i opacity from local structure, companion galaxies, or cosmic infall. We discuss our results in light of the numerical radiative transfer simulations, the evolving total Lyα and ionizing output of galaxies, and the utility of resolved Lyα spectra in the reionization epoch.

preprint2019arXiv

Resolved Lyman-α properties of a luminous Lyman-break galaxy in a large ionised bubble at z = 6.53

The observed properties of the Lyman-$α$ (Ly$α$) emission line are a powerful probe of neutral gas in and around galaxies. We present spatially resolved Ly$α$ spectroscopy with VLT/MUSE targeting VR7, a UV-luminous galaxy at $z=6.532$ with moderate Ly$α$ equivalent width (EW$_0\approx38$ Å). These data are combined with deep resolved [CII]$_{\rm 158 μm}$ spectroscopy obtained with ALMA and UV imaging from {\it HST} and we also detect UV continuum with MUSE. Ly$α$ emission is clearly detected with S/N $\approx40$ and FWHM of 374 km s$^{-1}$. Ly$α$ and [CII] are similarly extended beyond the UV, with effective radius r$_{\rm eff} = 2.1\pm0.2$ kpc for a single exponential model or r$_{\rm eff, Lyα, halo} = 3.45^{+1.08}_{-0.87}$ kpc when measured jointly with the UV continuum. The Ly$α$ profile is broader and redshifted with respect to the [CII] line (by 213 km s$^{-1}$), but there are spatial variations that are qualitatively similar in both lines and coincide with resolved UV components. This suggests that the emission originates from two components with plausibly different HI column densities. We place VR7 in the context of other galaxies at similar and lower redshift. The Ly$α$ halo scale length is similar at different redshifts and velocity shifts with respect to the systemic are typically smaller. Overall, we find little indications of a more neutral vicinity at higher redshift. This means that the local ($\sim 10$ kpc) neutral gas conditions that determine the observed Ly$α$ properties in VR7 resemble the conditions in post-re-ionisation galaxies.

preprint2019arXiv

Three Lyman-alpha emitting galaxies within a quasar proximity zone at z=5.8

Quasar proximity zones at $z>5.5$ correspond to over-dense and over-ionized environments. Galaxies found inside proximity zones can therefore display features which would otherwise be masked by absorption in the IGM. We demonstrate the utility of this quasar-galaxy synergy by reporting the discovery of the first three `proximate Lyman-$α$ emitters' (LAEs) within the proximity zone of quasar J0836 at $z=5.802$ (\textit{Aerith A, B} and \textit{C}). \textit{Aerith A}, located behind the quasar with an impact parameter $D_\perp = 278$ pkpc, provides the first detection of a Lyman-$α$ transverse proximity effect. We model the transmission and show it constrains the onset of J0836's quasar phase to $0.2 \text{Myr}<t<20\text{Myr}$ in the past. The second object, \textit{Aerith B} at a distance $D=750$ pkpc from the quasar, displays a bright, broad double-peaked \lal emission line. Based on relations calibrated at $z\leq3$, the peak separation implies a low ionizing $f_{\text{esc}} \lesssim 1\%$, the most direct such constraint on a reionization-era galaxy. We fit the Ly-$α$ line with an outflowing shell model, finding a completely typical central density $\text{log N}_{\text{HI}}/\text{cm}^{-2} = 19.3_{-0.2}^{+0.8}$, outflow velocity $v=16_{-11}^{+4}$ km s$^{-1}$, and gas temperature $\text{log} T/\text{K} = 3.8_{-0.7}^{+0.8}$ compared to $2<z<3$ analogue LAEs. Finally, we detect an emission line at $λ=8177$ Å in object \textit{Aerith C} which, if it is \lal at $z=5.726$, would correspond closely with the end of the quasar's proximity zone ($Δz<0.02$ from the boundary) and suggests the quasar influences the IGM up to $\sim85$ cMpc away, making it the largest quasar proximity zone. Via the analyses conducted here, we illustrate how proximate LAEs offer unique insight into the ionizing properties of both quasars and galaxies during the epoch of reionization.

preprint2016arXiv

Estimates of cluster masses in screened modified gravity

We use cosmological hydrodynamical simulations to study the effect of screened modified gravity models on the mass estimates of galaxy clusters. In particular, we focus on two novel aspects: (i) we study modified gravity models in which baryons and dark matter are coupled with different strengths to the scalar field, and, (ii) we put the simulation results into the greater context of a general screened-modified gravity parametrization. We have compared the mass of clusters inferred via lensing versus the mass inferred via kinematical measurements as a probe of violations of the equivalence principle at Mpc scales. We find that estimates of cluster masses via X-ray observations is mainly sensitive to the coupling between the scalar degree of freedom and baryons - while the kinematical mass is mainly sensitive to the coupling to dark matter. Therefore, the relation between the two mass estimates is a probe of a possible non-universal coupling between the scalar field, the standard model fields, and dark matter. Finally, we used observational data of kinetic, thermal and lensing masses to place constraints on deviations from general relativity on cluster scales for a general parametrization of screened modified gravity theories which contains $f(R)$ and Symmetron models. We find that while the kinematic mass can be used to place competitive constraints, using thermal measurements is challenging as a potential non-thermal contribution is degenerate with the imprint of modified gravity.

preprint2016arXiv

From Mirrors to Windows: Lyman-Alpha Radiative Transfer in a Very Clumpy Medium

Lyman-Alpha (Ly$α$) is the strongest emission line in the Universe and is frequently used to detect and study the most distant galaxies. Because Lya is a resonant line, photons typically scatter prior to escaping; this scattering process complicates the interpretation of Ly$α$ spectra, but also encodes a wealth of information about the structure and kinematics of neutral gas in the galaxy. Modeling the Ly$α$ line therefore allows us to study tiny-scale features of the gas, even in the most distant galaxies. Curiously, observed Ly$α$ spectra can be modeled successfully with very simple, homogeneous geometries (such as an expanding, spherical shell), whereas more realistic, multiphase geometries often fail to reproduce the observed spectra. This seems paradoxical since the gas in galaxies is known to be multiphase. In this Letter, we show that spectra emerging from extremely clumpy geometries with many clouds along the line of sight converge to the predictions from simplified, homogeneous models. We suggest that this resolves the apparent discrepancy, and may provide a way to study the gas structure in galaxies on scales far smaller than can be probed in either cosmological simulations or direct (i.e., spatially-resolved) observations.

preprint2016arXiv

Green Pea Galaxies Reveal Secrets of Ly$α$ Escape

We analyze archival Ly$α$ spectra of 12 "Green Pea" galaxies observed with the Hubble Space Telescope, model their Ly$α$ profiles with radiative transfer models, and explore the dependence of Ly$α$ escape fraction on various properties. Green Pea galaxies are nearby compact starburst galaxies with [OIII]$λ$5007 equivalent widths of hundreds of Å. All 12 Green Pea galaxies in our sample show Ly$α$ lines in emission, with a Ly$α$ equivalent width distribution similar to high redshift Ly$α$ emitters. Combining the optical and UV spectra of Green Pea galaxies, we estimate their Ly$α$ escape fractions and find correlations between Ly$α$ escape fraction and kinematic features of Ly$α$ profiles. The escape fraction of Ly$α$ in these galaxies ranges from 1.4% to 67%. We also find that the Ly$α$ escape fraction depends strongly on metallicity and moderately on dust extinction. We compare their high-quality Ly$α$ profiles with single HI shell radiative transfer models and find that the Ly$α$ escape fraction anti-correlates with the derived HI column densities. Single shell models fit most Ly$α$ profiles well, but not the ones with highest escape fractions of Ly$α$. Our results suggest that low HI column density and low metallicity are essential for Ly$α$ escape, and make a galaxy a Ly$α$ emitter.

preprint2016arXiv

Ly$α$ Spectra from Multiphase Outflows, and their Connection to Shell Models

We perform Lyman-$α$ (Ly$α$) Monte-Carlo radiative transfer calculations on a suite of $2500$ models of multiphase, outflowing media, which are characterized by $14$ parameters. We focus on the Ly$α$ spectra emerging from these media, and investigate which properties are dominant in shaping the emerging Ly$α$ profile. Multiphase models give rise to a wide variety of emerging spectra, including single, double and triple peaked spectra. We find that the dominant parameters in shaping the spectra include (i) the cloud covering factor, $f_c$, in agreement with earlier studies, and (ii) the temperature and number density of residual HI in the hot ionized medium. We attempt to reproduce spectra emerging from multiphase models with `shell models' which are commonly used to fit observed Ly$α$ spectra, and investigate the connection between shell-model parameters and the physical parameters of the clumpy media. In shell models, the neutral hydrogen content of the shell is one of the key parameters controlling Ly$α$ radiative transfer. Because Ly$α$ spectra emerging from multi-phase media depend much less on the neutral hydrogen content of the clumps, the shell model parameters such as HI column density (but also shell velocity and dust content) are generally not well matched to the associated physical parameters of clumpy media.

preprint2016arXiv

Lyman Alpha Signatures from Direct Collapse Black Holes

`Direct collapse black holes' (DCBHs) provide possible seeds for supermassive black holes that exist at redshifts as high as z~7. We study Lyman Alpha (Lya) radiative transfer through simplified representations of the DCBH-scenario. We find that gravitational heating of the collapsing cloud gives rise to a Lya cooling luminosity of up to ~ 1e38(M_gas/1e6 Msun)^2 erg/s. The Lya production rate can be significantly larger during the final stages of collapse, but collisional deexcitation efficiently suppresses the emerging Lya flux. Photoionization by a central source boosts the Lya luminosity to L~1e43(M_BH/1e6 M_sun) erg/s during specific evolutionary stages of the cloud, where M_BH denotes the mass of the black hole powering this source. We predict that the width and velocity off-set of the Lya spectral line range from a few tens to few thousands km/s, depending sensitively on the evolutionary state of the cloud. We also compare our predictions to observations of CR7 (Sobral et al. 2015), a luminous Lya emitter at z~7, which is potentially associated with a DCBH. If CR7 is powered by a black hole, then its Lya flux alone requires that M_BH> 1e7 M_sun, which exceeds the mass of DCBHs when they first form. The observed width of the Lya spectrum favors the presence of only a low column density of hydrogen, log [N_HI/cm^-2]~19-20. The shape of the Lya spectrum indicates that this gas is outflowing. These requirements imply that if CR7 harbors a DCBH, then the physical conditions that enabled its formation have been mostly erased, which is in agreement with theoretical expectations.

preprint2016arXiv

The LyA-LyC Connection: Evidence for an Enhanced Contribution of UV-faint Galaxies to Cosmic Reionization

The escape of ionizing Lyman Continuum (LyC) photons requires the existence of low-N_HI sightlines, which also promote escape of Lyman-Alpha (Lya). We use a suite of 2500 Lya Monte-Carlo radiative transfer simulations through models of dusty, clumpy interstellar (`multiphase') media from Gronke & Dijkstra (2016), and compare the escape fractions of Lya [f_esc(Lya)] and LyC radiation [f_esc(LyC)]. We find that f_esc(LyC) and f_esc(Lya) are correlated: galaxies with a low f_esc(Lya) consistently have a low f_esc(LyC), while galaxies with a high f_esc(Lya) exhibit a large dispersion in f_esc(LyC). We argue that there is increasing observational evidence that Lya escapes more easily from UV-faint galaxies. The correlation between f_esc(LyC) and f_esc(Lya) then implies that UV-faint galaxies contribute more to the ionizing background than implied by the faint-end slope of the UV-luminosity function. In multiphase gases, the ionizing escape fraction is most strongly affected by the cloud covering factor, f_cl, which implies that f_esc(LyC) is closely connected to the observed Lya spectral line shape. Specifically, LyC emitting galaxies typically having narrower, more symmetric line profiles. This prediction is qualitatively similar to that for `shell models'.

preprint2015arXiv

A systematic study of Lyman-Alpha transfer through outflowing shells: Model parameter estimation

Outflows promote the escape of Lyman-$α$ (Ly$α$) photons from dusty interstellar media. The process of radiative transfer through interstellar outflows is often modelled by a spherically symmetric, geometrically thin shell of gas that scatters photons emitted by a central Ly$α$ source. Despite its simplified geometry, this `shell model' has been surprisingly successful at reproducing observed Ly$α$ line shapes. In this paper we perform automated line fitting on a set of noisy simulated shell model spectra, in order to determine whether degeneracies exist between the different shell model parameters. While there are some significant degeneracies, we find that most parameters are accurately recovered, especially the HI column density ($N_{\rm HI}$) and outflow velocity ($v_{\rm exp}$). This work represents an important first step in determining how the shell model parameters relate to the actual physical properties of Ly$α$ sources. To aid further exploration of the parameter space, we have made our simulated model spectra available through an interactive online tool.

preprint2015arXiv

Blue Wings and Bumps via Fermi-like Acceleration of Lyman-alpha photons across Shocks

We explore the impact of Fermi-like acceleration of Lyman-alpha (Lyα) photons across shock fronts on the observed Lyα spectral line shape. We first confirm the result of Neufeld & McKee (1988) that this mechanism gives rise to extended blue wings which may have been observed in some radio galaxies. Our Monte-Carlo radiative transfer calculations further show that in a minor modification of the shell-model, in which we add an additional static shell of hydrogen, this process can naturally explain the small blue bumps observed in a subset of Lyα emitting galaxies, which have been difficult to explain with conventional shell-models. Blue bumps can be produced with an additional column density of static hydrogen as small as $N_{HI}^{static} \ll N_{HI}^{shell}$, and typically occur at roughly the outflow velocity of the shell. In our model the spectra of so-called 'blue-bump objects' might reflect an evolutionary stage in which the outflows regulating the escape of Lyα photons are still engulfed within a static interstellar medium.

preprint2015arXiv

Connecting Faint End Slopes of the Lyman-$α$ emitter and Lyman-break Galaxy Luminosity Functions

We predict Lyman-$α$ (Ly$α$) luminosity functions (LFs) of Ly$α$-selected galaxies (Ly$α$ emitters, or LAEs) at $z=3-6$ using the phenomenological model of Dijkstra & Wyithe (2012). This model combines observed UV-LFs of Lyman-break galaxies (LBGs, or drop out galaxies), with constraints on their distribution of Ly$α$ line strengths as a function of UV-luminosity and redshift. Our analysis shows that while Ly$α$ LFs of LAEs are generally not Schechter functions, these provide a good description over the luminosity range of $\log_{10}( L_α/{\rm erg}\,{\rm s}^{-1})=41-44$. Motivated by this result, we predict Schechter function parameters at $z=3-6$. Our analysis further shows that (i) the faint end slope of the Ly$α$ LF is steeper than that of the UV-LF of Lyman-break galaxies, (with a median $α_{Lyα} < -2.0$ at $z\gtrsim 4$), and (ii) a turn-over in the Ly$α$ LF of LAEs at Ly$α$ luminosities $10^{40}$ erg s$^{-1}<L_α\lesssim 10^{41}$ erg s$^{-1}$ may signal a flattening of UV-LF of Lyman-break galaxies at $-12>M_{\rm UV}>-14$. We discuss the implications of these results - which can be tested directly with upcoming surveys - for the Epoch of Reionization.

preprint2015arXiv

Universal predictions of screened modified gravity on cluster scales

Modified gravity models require a screening mechanism to be able to evade the stringent constraints from local gravity experiments and, at the same time, give rise to observable astrophysical and cosmological signatures. Such screened modified gravity models necessarily have dynamics determined by complex nonlinear equations that usually need to be solved on a model-by-model basis to produce predictions. This makes testing them a cumbersome process. In this paper, we investigate whether there is a common signature for all the different models that is suitable to testing them on cluster scales. To do this we propose an observable related to the fifth force, which can be observationally related to the ratio of dynamical-to-lensing mass of a halo, and then show that the predictions for this observable can be rescaled to a near universal form for a large class of modified gravity models. We demonstrate this using the Hu-Sawicki $f(R)$, the Symmetron, the nDGP, and the Dilaton models, as well as unifying parametrizations. The universal form is determined by only three quantities: a strength, a mass, and a width parameter. We also show how these parameters can be derived from a specific theory. This self-similarity in the predictions can hopefully be used to search for signatures of modified gravity on cluster scales in a model-independent way.

preprint2014arXiv

Directional Lya equivalent boosting I. Spherically symmetric distributions of clumps

We quantify the directional dependence of the escape fraction of Lyman-$α$ (Ly$α$) and non-ionising UV-continuum photons from a multiphase medium, and investigate whether there exist directional enhancements in the Ly$α$ equivalent width (EW). Our multiphase medium consists of spherically symmetric distributions of cold, dusty clumps embedded within a hot dust-free medium. We focus on three models from the analysis presented by Laursen et al. (2013). We find that for a Ly$α$ and UV-continuum point source, it is possible to find an EW boost $b(θ,ϕ) > 5 \bar{b}$ in a few per cent of sight lines, where $\bar{b}$ denotes the boost averaged over all photons. For spatially extended sources this directional dependence vanishes quickly when the size of the UV emitting region exceeds the mean distance between cold dusty clumps. Our analysis suggests that directional EW boosting can occur, and that this is mostly driven by reduced escape fractions of UV photons (which gives rise to UV-continuum `shadows'), and less due to an enhanced Ly$α$ escape fraction (or beaming), in certain directions.

preprint2014arXiv

Halo velocity profiles in screened modified gravity theories

Screened modified gravity predicts potentially large signatures in the peculiar velocity field that makes it an interesting probe to test gravity on cosmological scales. We investigate the signatures induced by the Symmetron and a Chameleon $f(R)$ model in the peculiar velocity field using $N$-body simulations. By studying fifth force and halo velocity profiles we identify three general categories of effects found in screened modified gravity models: a fully screened regime where we recover $Λ$CDM to high precision, an unscreened regime where the fifth force is in full operation, and, a partially screened regime where screening occurs in the inner part of a halo, but the fifth force is active at larger radii. These three regimes can be pointed out very clearly by analyzing the deviation in the maximum cluster velocity. Observationally, the partially screened regime is of particular interest since an uniform increase of the gravitational force - as present in the unscreened regime - is degenerate with the (dynamical) halo mass estimate, and, thus, hard to detect.