Researcher profile

Simon J. Lilly

Simon J. Lilly contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
2topics
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

9 published item(s)

preprint2024arXiv

EIGER IV: The cool 10$^4$K circumgalactic environment of high-$z$ galaxies reveals remarkably efficient IGM enrichment

We report new observations of the cool diffuse gas around 29, $2.3<z<6.3$ galaxies, using deep JWST/NIRCam slitless grism spectroscopy around the sightline to the quasar J0100+2802. The galaxies span a stellar mass range of $7.1 \leq \log M_{*}/M_{sun} \leq 10.7$, and star-formation rates of $-0.1 < \log \; SFR/M_{sun}yr^{-1} \; <2.3$. We find galaxies for seven MgII absorption systems within 300 kpc of the quasar sightline. The MgII radial absorption profile falls off sharply with radii, with most of the absorption extending out to 2-3$R_{200}$ of the host galaxies. Six out of seven MgII absorption systems are detected around galaxies with $\log M_{*}/M_{sun} >$9. MgII absorption kinematics are shifted from the systemic redshift of host galaxies with a median absolute velocity of 135 km/s and standard deviation of 85 km/s. The high kinematic offset and large radial separation ($R> 1.3 R_{200}$), suggest that five out of the seven MgII absorption systems are gravitationally not bound to the galaxies. In contrast, most cool circumgalactic media at $z<1$ are gravitationally bound. The high incidence of unbound MgII gas in this work suggests that towards the end of reionization, galaxy halos are in a state of remarkable disequilibrium, and are highly efficient in enriching the intergalactic medium. Two strongest MgII absorption systems are detected at $z\sim$ 4.22 and 4.5, the former associated with a merging galaxy system and the latter associated with three kinematically close galaxies. Both these galaxies reside in local galaxy over-densities, indicating the presence of cool MgII absorption in two &#34;proto-groups&#34; at $z>4$.

preprint2022arXiv

The Gas-phase Metallicity Profiles of Star-forming Galaxies in the Modified Accretion Disk Framework

Simulations indicate that the inflow of gas of star-forming galaxies is almost co-planar and co-rotating with the gas disk, and that the outflow of gas driven by stellar winds and/or supernova explosions is preferentially perpendicular to the disk. This indicates that the galactic gas disk can be treated as a modified accretion disk. In this work, we focus on the metal enhancement in galactic disks in this scenario of gas accretion. Assuming that the star formation rate surface density ($Σ_{\rm SFR}$) is of exponential form, we obtain the analytic solution of gas-phase metallicity with only three free parameters: the scalelength of $Σ_{\rm SFR}$ ($h_{\rm R}$), the metallicity of the inflowing gas and the mass-loading factor defined as the wind-driven outflow rate surface density per $Σ_{\rm SFR}$. According to this simple model, the negative gradient of gas-phase metallicity is a natural consequence of the radial inflow of cold gas which is continuously enriched by in-situ star formation as it moves towards the disk center. We fit the model to the observed metallicity profiles for six nearby galaxies chosen to have well-measured metallicity profiles extending to very large radii. Our model can well characterize the overall features of the observed metallicity profiles. The observed profiles usually show a floor at the outer regions of the disk, corresponding to the metallicity of inflow gas. Furthermore, we find the $h_{\rm R}$ of $Σ_{\rm SFR}$ inferred from these fits agree well with independent estimates from $Σ_{\rm SFR}$ profiles, supporting the basic model.

preprint2022arXiv

The Origin of Exponential Star-forming Disks

The disk components of galaxies generally show an exponential profile extending over several scale lengths, both in mass and star-formation rate, but the physical origin is not well understood. We explore a physical model in which the galactic gas disk is viewed as a &#34;modified accretion disk&#34; in which coplanar gas inflow, driven by viscous stresses in the disk, provides the fuel for star formation, which progressively removes gas as it flows inwards. We show that magnetic stresses from magneto-rotational instability are the most plausible source of the required viscosity, and construct a simple physical model to explore this. A key feature is to link the magnetic field strength to the local star-formation surface density, $B_{\rm tot} \propto Σ_{\rm SFR}^α$. This provides a feed-back loop between star-formation and the flow of gas. We find that the model naturally produces stable steady-state exponential disks, as long as $α\sim$ 0.15, the value indicated from spatially-resolved observations of nearby galaxies. The disk scale-length $h_{\rm R}$ is set by the rate at which the disk is fed, by the normalization of the $B_{\rm tot}-Σ_{\rm SFR}$ relation and by the circular velocity of the halo. The angular momentum distribution of the gas and stars within the disk is a consequence of the transfer of angular momentum that is inherent to the operation of an accretion disk, rather than the initial angular momentum of the inflowing material. We suggest that magnetic stresses likely play a major role in establishing the stable exponential form of galactic disks.

preprint2021arXiv

The stellar mass versus stellar metallicity relation of star-forming galaxies at $1.6\le z\le3.0$ and implications for the evolution of the $α$-enhancement

We measure the relationship between stellar mass and stellar metallicity, the stellar mass--metallicity relation (MZR), for 1336 star-forming galaxies at $1.6\le z\le3.0$ (<z>=2.2) using rest-frame far-ultraviolet spectra from the zCOSMOS-deep survey. High signal-to-noise composite spectra containing stellar absorption features are fit with population synthesis model spectra of a range of metallicity. We find stellar metallicities, which mostly reflect iron abundances, scaling as $(Z_{Fe,\ast}/Z_{Fe,\odot})=-(0.81\pm0.01)+(0.32+0.03)\log(M_\ast/10^{10}M_\odot)$ across the mass range of $10^9\lesssim M_\ast/M_\odot\lesssim10^{11}$, being $\approx6\times$ lower than seen locally at the same masses. The instantaneous oxygen-to-iron ratio ($α$-enhancement) inferred using the gas-phase oxygen MZRs, is on average found to be [O/Fe]$\approx0.47$, being higher than the local [O/Fe]$\approx0$. The observed changes in [O/Fe] and [Fe/H] are reproduced in simple flow-through gas-regulator models with steady star-formation histories (SFHs) that follow the evolving main sequence. Our models show that the [O/Fe] is determined almost entirely by the instantaneous specific star formation rate alone while being independent of the SFHs, mass, and the gas-regulation characteristics of the systems. We find that the locations of $\sim10^{10}M_\odot$ galaxies at z~2 in the [O/Fe]--metallicity planes are in remarkable agreement with the sequence of low-metallicity thick-disk stars in our Galaxy. This manifests a beautiful concordance between the results of Galactic archaeology and observations of high-redshift Milky Way progenitors. However, there remains a question of how and when the old metal-rich, low-$α$/Fe stars seen in the bulge had formed by z~2 because such a stellar population is not seen in our data and difficult to explain in the context of our models.

preprint2020arXiv

Evidence for a highly opaque large-scale galaxy void at the end of reionization

We present evidence that a region of high effective Ly$α$ optical depth at $z\sim5.7$ is associated with an underdense region at the tail end of cosmic reionization. We carried out a survey of Lyman-break Galaxies (LBGs) using Subaru Hyper Suprime-Cam in the field of the $z=5.98$ quasar J0148+0600, whose spectrum presents an unusually long ($\sim160 \mathrm{cMpc}$) and opaque ($τ\gtrsim7$) Ly$α$ trough at $5.5\le z\le 5.9$. LBG candidates were selected to lie within the redshift range of the trough, and the projected number densities were measured within 90~cMpc of the quasar sightline. The region within $8&#39;$ (or $\approx 19~\mathrm{cMpc}$) of the quasar position is the most underdense of the whole field. The significance of the presence of the void is estimated to be $99\%$. This is consistent with the significant deficit of Ly$α$ emitters (LAEs) at $z=5.72$ reported by Becker et al. and suggests that the paucity of LAEs is not purely due to the removal of the Ly$α$ emission by the high opacity but reflects a real coherent underdensity of galaxies across the entire redshift range of the trough. These observations are consistent with scenarios in which large optical depth fluctuations arise due to fluctuations in the galaxy-dominant UV background or due to residual neutral islands that are expected from reionization that is completed at redshifts as low as $z\lesssim5.5$.

preprint2020arXiv

From haloes to galaxies -- II. The fundamental relations in star formation and quenching

Star formation and quenching are two of the most important processes in galaxy formation and evolution. We explore in the local Universe the interrelationships among key integrated galaxy properties, including stellar mass $M_*$, star formation rate (SFR), specific SFR (sSFR), molecular gas mass $M_{\rm H_2}$, star formation efficiency (SFE) of the molecular gas and molecular gas to stellar mass ratio $μ$. We aim to identify the most fundamental scaling relations among these key galaxy properties and their interrelationships. We show the integrated $M_{\rm H_2}$-SFR, SFR-$M_*$ and $M_{\rm H_2}$-$M_*$ relation can be simply transformed from the $μ$-sSFR, SFE-$μ$ and SFE-sSFR relation, respectively. The transformation, in principle, can increase or decrease the scatter of each relation. Interestingly, we find the latter three relations all have significantly smaller scatter than the former three corresponding relations. We show the probability to achieve the observed small scatter by accident is extremely close to zero. This suggests that the smaller scatters of the latter three relations are driven by a more fundamental physical connection among these quantities. We then show the large scatters in the former relations are due to their systematic dependence on other galaxy properties, and on star formation and quenching process. We propose the sSFR-$μ$-SFE relation as the Fundamental Formation Relation (FFR), which governs the star formation and quenching process, and provides a simple framework to study galaxy evolution. Other scaling relations, including integrated Kennicutt-Schmidt law, star-forming main sequence and molecular gas main sequence, can all be derived from the FFR.

preprint2020arXiv

LATIS: The Ly$α$ Tomography IMACS Survey

We introduce LATIS, the Ly$α$ Tomography IMACS Survey, a spectroscopic survey at Magellan designed to map the z=2.2-2.8 intergalactic medium (IGM) in three dimensions by observing the Ly$α$ forest in the spectra of galaxies and QSOs. Within an area of 1.7 deg${}^2$, we will observe approximately half of $\gtrsim L^*$ galaxies at z=2.2-3.2 for typically 12 hours, providing a dense network of sightlines piercing the IGM with an average transverse separation of 2.5 $h^{-1}$ comoving Mpc (1 physical Mpc). At these scales, the opacity of the IGM is expected to be closely related to the dark matter density, and LATIS will therefore map the density field in the $z \sim 2.5$ universe at $\sim$Mpc resolution over the largest volume to date. Ultimately LATIS will produce approximately 3800 spectra of z=2.2-3.2 galaxies that probe the IGM within a volume of $4 \times 10^6 h^{-3}$ Mpc${}^3$, large enough to contain a representative sample of structures from protoclusters to large voids. Observations are already complete over one-third of the survey area. In this paper, we describe the survey design and execution. We present the largest IGM tomographic maps at comparable resolution yet made. We show that the recovered matter overdensities are broadly consistent with cosmological expectations based on realistic mock surveys, that they correspond to galaxy overdensities, and that we can recover structures identified using other tracers. LATIS is conducted in Canada-France-Hawaii Telescope Legacy Survey fields, including COSMOS. Coupling the LATIS tomographic maps with the rich data sets collected in these fields will enable novel studies of environment-dependent galaxy evolution and the galaxy-IGM connection at cosmic noon.

preprint2020arXiv

The Variability of Star Formation Rate in Galaxies: II. Power Spectrum Distribution on the Main Sequence

We constrain the temporal power spectrum of the sSFR(t) of star-forming galaxies, using a well-defined sample of Main Sequence galaxies from MaNGA and our earlier measurements of the ratio of the SFR averaged within the last 5 Myr to that averaged over the last 800 Myr. We explore the assumptions of stationarity and ergodicity that are implicit in this approach. We assume a single power-law form of the PSD but introduce an additional free parameter, the &#34;intrinsic scatter&#34;, to try to account for any non-ergodicity introduced from various sources. We analyze both an &#34;integrated&#34; sample consisting of global measurements of all of the galaxies, and also 25 sub-samples obtained by considering five radial regions and five bins of integrated stellar mass. Assuming that any intrinsic scatter is not the dominant contribution to the Main Sequence dispersion of galaxies, we find that the PSDs have slopes between 1.0 and 2.0, indicating that the power (per log interval of frequency) is mostly contributed by longer timescale variations. We find a correlation between the returned PSDs and the inferred gas depletion times ($τ_{\rm dep,eff}$) obtained from application of the extended Schmidt Law, in that regions with shorter gas depletion times show larger integrated power and flatter PSD. Intriguingly, it is found that shifting the PSDs by the inferred $τ_{\rm dep,eff}$ causes all of the 25 PSDs to closely overlap, at least in that region where the PSD is best constrained and least affected by uncertainties about any intrinsic scatter. A possible explanation of these results is the dynamical response of the gas regulator system of Lilly et al. 2013 to a uniform time-varying inflow, as previously proposed in Wang et al. 2019.

preprint2020arXiv

The Variability of the Star Formation Rate in Galaxies: I. Star Formation Histories Traced by EW(H$α$) and EW(H$δ_A$)

To investigate the variability of the star formation rate (SFR) of galaxies, we define a star formation change parameter, SFR$_{\rm 5Myr}$/SFR$_{\rm 800Myr}$ which is the ratio of the SFR averaged within the last 5 Myr to the SFR averaged within the last 800 Myr. We show that this parameter can be determined from a combination of H$α$ emission and H$δ$ absorption, plus the 4000 A break, with an uncertainty of $\sim$0.07 dex for star-forming galaxies. We then apply this estimator to MaNGA galaxies, both globally within Re and within radial annuli. We find that galaxies with higher global SFR$_{\rm 5Myr}$/SFR$_{\rm 800Myr}$ appear to have higher SFR$_{\rm 5Myr}$/SFR$_{\rm 800Myr}$ at all galactic radii, i.e. that galaxies with a recent temporal enhancement in overall SFR have enhanced star formation at all galactic radii. The dispersion of the SFR$_{\rm 5Myr}$/SFR$_{\rm 800Myr}$ at a given relative galactic radius and a given stellar mass decreases with the (indirectly inferred) gas depletion time: locations with short gas depletion time appear to undergo bigger variations in their star-formation rates on Gyr or less timescales. In Wang et al. (2019) we showed that the dispersion in star-formation rate surface densities $Σ_{\rm SFR}$ in the galaxy population appears to be inversely correlated with the inferred gas depletion timescale and interpreted this in terms of the dynamical response of a gas-regulator system to changes in the gas inflow rate. In this paper, we can now prove directly with SFR$_{\rm 5Myr}$/SFR$_{\rm 800Myr}$ that these effects are indeed due to genuine temporal variations in the SFR of individual galaxies on timescales between $10^7$ and $10^9$ years rather than possibly reflecting intrinsic, non-temporal, differences between different galaxies.