Source author record

Simon J. Lilly

Simon J. Lilly 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

30works
4topics
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

30 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 "proto-groups" 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 "modified accretion disk" 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'$ (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 "intrinsic scatter", to try to account for any non-ergodicity introduced from various sources. We analyze both an "integrated" 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.

preprint2016arXiv

Constraining the lifetime and opening angle of quasars using fluorescent Lyman $α$ emission: the case of Q0420-388

A toy model is developed to understand how the spatial distribution of fluorescent emitters in the vicinity of bright quasars could be affected by the geometry of the quasar bi-conical radiation field and by its lifetime. The model is then applied to the distribution of high equivalent width Lyman $α$ emitters (with rest-frame equivalent widths above 100 A, threshold used in e.g. Trainor & Steidel, 2013) identified in a deep narrow-band 36x36 arcmin$^2$ image centered on the luminous quasar Q0420-388. These emitters are found to the edge of the field and show some evidence of an azimuthal asymmetry on the sky of the type expected if the quasar is radiating in a bipolar cone. If these sources are being fluorescently illuminated by the quasar, the two most distant objects require a lifetime of at least 15 Myr for an opening angle of 60 degrees or more, increasing to more than 40 Myr if the opening angle is reduced to a minimum 30 degrees. However, few of the other expected signatures of boosted fluorescence are seen at the current survey limits, e.g. a fall off in Lyman $α$ brightness, or equivalent width, with distance. Furthermore, to have most of the Lyman $α$ emission of the two distant sources to be fluorescently boosted would require the quasar to have been significantly brighter in the past. This suggests that these particular sources may not be fluorescent, invalidating the above lifetime constraints. This would cast doubt on the use of this relatively low equivalent width threshold and thus on the lifetime analysis in Trainor & Steidel (2013).

preprint2016arXiv

Large scale structure around a z=2.1 cluster

The most prodigious starburst galaxies are absent in massive galaxy clusters today, but their connection with large scale environments is less clear at $z\gtrsim2$. We present a search of large scale structure around a galaxy cluster core at $z=2.095$ using a set of spectroscopically confirmed galaxies. We find that both color-selected star-forming galaxies (SFGs) and dusty star-forming galaxies (DSFGs) show significant overdensities around the $z=2.095$ cluster. A total of 8 DSFGs (including 3 X-ray luminous active galactic nuclei, AGNs) and 34 SFGs are found within a 10 arcmin radius (corresponds to $\sim$15 cMpc at $z\sim2.1$) from the cluster center and within a redshift range of $Δz=0.02$, which leads to galaxy overdensities of $δ_{\rm DSFG}\sim12.3$ and $δ_{\rm SFG}\sim2.8$. The cluster core and the extended DSFG- and SFG-rich structure together demonstrate an active cluster formation phase, in which the cluster is accreting a significant amount of material from large scale structure while the more mature core may begin to virialize. Our finding of this DSFG-rich structure, along with a number of other protoclusters with excess DSFGs and AGNs found to date, suggest that the overdensities of these rare sources indeed trace significant mass overdensities. However, it remains puzzling how these intense star formers are triggered concurrently. Although an increased probability of galaxy interactions and/or enhanced gas supply can trigger the excess of DSFGs, our stacking analysis based on 850 $μ$m images and morphological analysis based on rest-frame optical imaging do not show such enhancements of merger fraction and gas content in this structure.

preprint2016arXiv

Minor Mergers or Progenitor Bias? The Stellar Ages of Small and Large Quenched Early-Type Galaxies

We investigate the origin of the evolution of the population-averaged size of quenched galaxies (QGs) through a spectroscopic analysis of their stellar ages. The two most favoured scenarios for this evolution are either the size growth of individual galaxies through a sequence of dry minor merger events, or the addition of larger, newly quenched galaxies to the pre-existing population (i.e., a progenitor bias effect). We use the 20k zCOSMOS-bright spectroscopic survey to select bona fide quiescent galaxies at 0.2<z<0.8. We stack their spectra in bins of redshift, stellar mass and size to compute stellar population parameters in these bins through fits to the rest-frame optical spectra and through Lick spectral indices. We confirm a change of behaviour in the size-age relation below and above the ~10^11 MSun stellar mass scale: In our 10.5 < log M*/MSun < 11 mass bin, over the entire redshift window, the stellar populations of the largest galaxies are systematically younger than those of the smaller counterparts, pointing at progenitor bias as the main driver of the observed average size evolution at sub-10^11 MSun masses. In contrast, at higher masses, there is no clear trend in age as a function of galaxy size, supporting a substantial role of dry mergers in increasing the sizes of these most massive QGs with cosmic time. Within the errors, the [alpha/Fe] abundance ratios of QGs are (i) above-solar over the entire redshift range of our analysis, hinting at universally short timescales for the buildup of the stellar populations of QGs, and (ii) similar at all masses and sizes, suggesting similar (short) timescales for the whole QG population and strengthening the role of mergers in the buildup of the most massive QGs in the Universe.

preprint2016arXiv

Surface density effects in quenching: cause or effect?

There are very strong observed correlations between the specific star-formation rates (sSFR) of galaxies and their mean surface mass densities, Σ, as well as other aspects of their internal structure. These strong correlations have often been taken to indicate that the internal structure of a galaxy must play a major physical role, directly or indirectly, in the control of star-formation. In this paper we show by means of a very simple toy model that these correlations can arise naturally without any such physical role once the observed evolution of the size-mass relation for star-forming galaxies is taken into account. In particular, the model reproduces the sharp threshold in Σ between galaxies that are star-forming and those that are quenched, and the evolution of this threshold with redshift. Similarly, it produces iso-quenched-fraction contours in the ${f_Q(m,R_e)}$ plane that are almost exactly parallel to lines of constant Σ for centrals and shallower for satellites. It does so without any dependence on quenching on size or Σ, and without invoking any differences between centrals and satellites, beyond the different mass-dependences of their quenching laws. The toy-model also reproduces several other observations, including the sSFR gradients within galaxies and the appearance of inside-out build-up of passive galaxies. Finally, it is shown that curvature in the Main Sequence sSFR-mass relation can produce curvature in the apparent B/T ratios with mass. Our analysis therefore suggests that many of the strong correlations that are observed between galaxy structure and sSFR may well be a consequence of things unrelated to quenching and should not be taken as evidence of the physical processes that drive quenching.

preprint2016arXiv

Ubiquitous giant Ly $α$ nebulae around the brightest quasars at $z\sim3.5$ revealed with MUSE

Direct Ly $α$ imaging of intergalactic gas at $z\sim2$ has recently revealed giant cosmological structures around quasars, e.g. the Slug Nebula (Cantalupo et al. 2014). Despite their high luminosity, the detection rate of such systems in narrow-band and spectroscopic surveys is less than 10%, possibly encoding crucial information on the distribution of gas around quasars and the quasar emission properties. In this study, we use the MUSE integral-field instrument to perform a blind survey for giant Ly $α$ nebulae around 17 bright radio-quiet quasars at $3<z<4$ that does not suffer from most of the limitations of previous surveys. After data reduction and analysis performed with specifically developed tools, we found that each quasar is surrounded by giant Ly $α$ nebulae with projected sizes larger than 100 physical kpc and, in some cases, extending up to 320 kpc. The circularly averaged surface brightness profiles of the nebulae appear very similar to each other despite their different morphologies and are consistent with power laws with slopes $\approx-1.8$. The similarity between the properties of all these nebulae and the Slug Nebula suggests a similar origin for all systems and that a large fraction of gas around bright quasars could be in a relatively "cold" (T$\sim$10$^4$K) and dense phase. In addition, our results imply that such gas is ubiquitous within at least 50 kpc from bright quasars at $3<z<4$ independently of the quasar emission opening angle, or extending up to 200 kpc for quasar isotropic emission.

preprint2016arXiv

ZAP -- Enhanced PCA Sky Subtraction for Integral Field Spectroscopy

We introduce Zurich Atmosphere Purge (ZAP), an approach to sky subtraction based on principal component analysis (PCA) that we have developed for the Multi Unit Spectrographic Explorer (MUSE) integral field spectrograph. ZAP employs filtering and data segmentation to enhance the inherent capabilities of PCA for sky subtraction. Extensive testing shows that ZAP reduces sky emission residuals while robustly preserving the flux and line shapes of astronomical sources. The method works in a variety of observational situations from sparse fields with a low density of sources to filled fields in which the target source fills the field of view. With the inclusion of both of these situations the method is generally applicable to many different science cases and should also be useful for other instrumentation. ZAP is available for download at http://muse-vlt.eu/science/tools.

preprint2015arXiv

AGN evolution from a galaxy evolution viewpoint

We explore the connections between the evolving galaxy and AGN populations. We present a simple phenomenological model that links the evolving galaxy mass function and the evolving quasar luminosity function, which makes specific and testable predictions for the distribution of host galaxy masses for AGN of different luminosities. We show that the $ϕ^{*}$ normalisations of the galaxy mass function and of the AGN luminosity function closely track each other over a wide range of redshifts, implying a constant "duty cycle" of AGN activity. The strong redshift evolution in the AGN $L^*$ can be produced by either an evolution in the distribution of Eddington ratios, or in the $m_{bh}/m_{*}$ mass ratio, or both. To try to break this degeneracy we look at the distribution of AGN in the SDSS ($m_{bh},L$) plane, showing that an evolving ratio $m_{bh}/m_{*} \propto (1+z)^2$ reproduces the observed data and also reproduces the local relations which connect the black hole population with the host galaxies for both quenched and star-forming populations. We stress that observational studies that compare the masses of black holes in active galaxies at high redshift with those in quiescent galaxies locally will always see much weaker evolution. Evolution of this form would produce, or could be produced by, a redshift-independent $m_{bh} - σ$ relation and could explain why the local $m_{bh} - σ$ relation is tighter than $m_{bh} - m_{*}$ even if $σ$ is not directly linked to black hole growth. Irrespective of the evolution of $m_{bh}/m_{*}$, the model reproduces both the appearance of "downsizing" and the so-called "sub-Eddington boundary" without any mass-dependence in the evolution of black hole growth rates.

preprint2015arXiv

Quenching of Star Formation in SDSS Groups: Centrals, Satellites, and Galactic Conformity

We re-examine the fraction of low-redshift Sloan Digital Sky Survey satellites and centrals in which star formation has been quenched, using the environment quenching efficiency formalism that separates out the dependence of stellar mass. We show that the centrals of the groups containing the satellites are responding to the environment in the same way as their satellites (at least for stellar masses above 10^10.3 Msun), and that the well-known differences between satellites and the general set of centrals arise because the latter are overwhelmingly dominated by isolated galaxies. The widespread concept of "satellite quenching" as the cause of environmental effects in the galaxy population can therefore be generalized to "group quenching". We then explore the dependence of the quenching efficiency of satellites on overdensity, group-centric distance, halo mass, the stellar mass of the satellite, and the stellar mass and specific star formation rate (sSFR) of its central, trying to isolate the effect of these often interdependent variables. We emphasize the importance of the central sSFR in the quenching efficiency of the associated satellites, and develop the meaning of this "galactic conformity" effect in a probabilistic description of the quenching of galaxies. We show that conformity is strong, and that it varies strongly across parameter space. Several arguments then suggest that environmental quenching and mass quenching may be different manifestations of the same underlying process. The marked difference in the apparent mass dependencies of environment quenching and mass quenching which produces distinctive signatures in the mass functions of centrals and satellites will arise naturally, since, for satellites at least, the distributions of the environmental variables that we investigate in this work are essentially independent of the stellar mass of the satellite.

preprint2014arXiv

A Consistent Study of Metallicity Evolution at 0.8 < z < 2.6

We present the correlations between stellar mass, star formation rate (SFR) and [NII]/Ha flux ratio as indicator of gas-phase metallicity for a sample of 222 galaxies at 0.8 < z < 2.6 and log(M*/Msun)=9.0-11.5 from the LUCI, SINS/zC-SINF and KMOS3D surveys. This sample provides a unique analysis of the mass-metallicity relation (MZR) over an extended redshift range using consistent data analysis techniques and strong-line metallicity indicator. We find a constant slope at the low-mass end of the relation and can fully describe its redshift evolution through the evolution of the characteristic turnover mass where the relation begins to flatten at the asymptotic metallicity. At fixed mass and redshift, our data do not show a correlation between the [NII]/Ha ratio and SFR, which disagrees with the 0.2-0.3dex offset in [NII]/Ha predicted by the "fundamental relation" between stellar mass, SFR and metallicity discussed in recent literature. However, the overall evolution towards lower [NII]/Ha at earlier times does broadly agree with these predictions.

preprint2014arXiv

Mass and Environment as Drivers of Galaxy Evolution III: The constancy of the faint-end slope and the merging of galaxies

We explore using our continuity approach the underlying connections between the evolution of the faint-end slope of the stellar mass function of star-forming galaxies, the logarithmic slope of the sSFR-mass relation and the merging of galaxies. We derive analytically the consequences of the observed constancy of the faint-end slope since redshifts of at least z ~ 2. If the logarithmic slope of the sSFR-mass relation is negative, then the faint-end slope should quickly diverge due to the differential mass increase of galaxies on the star-forming main sequence, and this will also quickly destroy the Schechter form of the mass function. This problem can be solved by removing low mass galaxies by merging them into more massive galaxies. We quantify this process by introducing the specific merger mass rate (sMMR) as the specific rate of mass added to a given galaxy through mergers. For a modest negative value of the logarithmic slope of the sSFR-mass relation of beta ~ -0.1, an average sMMR ~ 0.1sSFR across the population is required to keep the faint-end slope constant with epoch, as observed. This in turn implies a merger rate of ~ 0.2sSFR for major mergers, which is consistent with the available observational estimates. More negative values of beta require higher sMMR and higher merger rates, and the steepening of the mass function becomes impossible to control for beta < ~ -0.6, for an observed value of the faint-end slope of ~ -1.4. The close link that is required between the in situ sSFR and the sMMR probably arises because both are closely linked to the buildup of dark matter haloes. These new findings further develop the formalism for the evolving galaxy population that we introduced earlier, and show how striking symmetries in the galaxy population can emerge as the result of deep links between the physical processes involved.

preprint2013arXiv

A comparative analysis of virial black-hole mass estimates of moderate-luminosity active galactic nuclei using Subaru/FMOS

We present an analysis of broad emission lines observed in moderate-luminosity active galactic nuclei (AGNs), typical of those found in X-ray surveys of deep fields, with the aim to test the validity of single-epoch virial black hole mass estimates. We have acquired near-infrared (NIR) spectra of AGNs up to z ~ 1.8 in the COSMOS and Extended Chandra Deep Field-South Survey, with the Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope. These low-resolution NIR spectra provide a significant detection of the broad Halpha line that has been shown to be a reliable probe of black hole mass at low redshift. Our sample has existing optical spectroscopy which provides a detection of MgII, a broad emission line typically used for black hole mass estimation at z > 1. We carry out a spectral-line fitting procedure using both Halpha and MgII to determine the virial velocity of gas in the broad line region, the monochromatic continuum luminosity at 3000 A, and the total Halpha line luminosity. With a sample of 43 AGNs spanning a range of two decades in luminosity (i.e., L ~ 10^44-46 ergs/s), we find a tight correlation between the continuum and line luminosity with a distribution characterized by <log(L_3000/L_Halpha)> = 1.52 and a dispersion sigma = 0.16. There is also a close one-to-one relationship between the FWHM of Halpha and of MgII up to 10000 km/s with a dispersion of 0.14 in the distribution of the logarithm of their ratios. Both of these then lead to there being very good agreement between Halpha- and MgII-based masses over a wide range in black hole mass (i.e., M_BH ~ 10^7-9 M_sun). We do find a small offset in MgII-based masses, relative to those based on Halpha, of +0.17 dex and a dispersion sigma = 0.32. In general, these results demonstrate that local scaling relations, using MgII or Halpha, are applicable for AGN at moderate luminosities and up to z ~ 2.

preprint2013arXiv

Gas-regulation of galaxies: the evolution of the cosmic sSFR, the metallicity-mass-SFR relation and the stellar content of haloes

A very simple physical model of galaxies, in which the formation of stars is instantaneously regulated by the mass of gas in a reservoir, links together three different aspects of the evolving galaxy population:(a) the cosmic time evolution of the specific star-formation rate sSFR relative to the growth of haloes, (b) the gas-phase metallicities across the galaxy population and over cosmic time, and (c) the ratio of the stellar to dark matter mass of haloes. If the SFR efficiency and wind mass loading are constant, the sSFR is set to the specific accretion rate of the galaxy: more realistic situations lead to an sSFR which is perturbed from this identity. The metallicity is set by the instantaneous operation of the regulator system rather than by the past history of the system. The regulator system naturally produces a Z(mstar, SFR) relation, with SFR as a second parameter in the mass-metallicity relation. This will be the same at all epochs unless the efficiency and mass-loading change with time, naturally producing a so-called "fundamental metallicity relation". The observed Z(mstar) relation of SDSS galaxies implies a strong dependence of stellar mass on halo mass that reconciles the different faint end slopes of the stellar and halo mass-functions in standard ΛCDM. It also boosts the sSFR relative to the specific accretion rate and produces a different dependence on mass, both of which are observed. The derived Z(mstar, SFR) relation for the regulator system is fit to published Z(mstar, SFR) data. The fitted efficiency is consistent with observed molecular gas-depletion timescales in galaxies while the fitted mass-loading is also plausible. The model also successfully reproduces the mass-metallicity relation of star-forming galaxies at z ~ 2.

preprint2013arXiv

The SINS/zC-SINF Survey of z~2 Galaxy Kinematics: The Nature of Dispersion Dominated Galaxies

We analyze the spectra, spatial distributions and kinematics of Ha, [NII] and [SII] emission in a sample of 42, z~2.2 UV/optically selected star forming galaxies (SFGs) from the SINS & zC-SINF surveys, 35 of which were observed in the adaptive optics mode of SINFONI. This is supplemented by kinematic data from 48 z~1-2.5 galaxies from the literature. We find that the kinematic classification of the high-z SFGs as `dispersion dominated' or `rotation dominated' correlates most strongly with their intrinsic sizes. Smaller galaxies are more likely `dispersion-dominated' for two main reasons: 1) The rotation velocity scales linearly with galaxy size but intrinsic velocity dispersion does not depend on size, and as such, their ratio is systematically lower for smaller galaxies, and 2) Beam smearing strongly decreases large-scale velocity gradients and increases observed dispersion much more for galaxies with sizes at or below the resolution. Dispersion dominated SFGs may thus have intrinsic properties similar to `rotation dominated' SFGs, but are primarily more compact, lower mass, less metal enriched and may have higher gas fractions, plausibly because they represent an earlier evolutionary state.

preprint2012arXiv

Detection of dark galaxies and circum-galactic filaments fluorescently illuminated by a quasar at z=2.4

A deep narrow-band survey for Ly-alpha emission carried out on the VLT-FORS2 has revealed 98 Ly-alpha candidates down to a flux limit of 4.e-18 erg/s/cm^2 in a volume of 5500 comoving Mpc^3 at z=2.4 centered on the hyperluminous quasar HE0109-3518. The properties of the detected sources in terms of their i) equivalent width distribution, ii) luminosity function, and iii) the average luminosity versus projected distance from the quasar, all suggest that a large fraction of these objects have been fluorescently "illuminated" by HE0109-3518. This conclusion is supported by comparison with detailed radiative transfer simulations of the effects of the quasar illumination. 18 objects have a rest-frame Equivalent Width (EW0) larger than 240A, the expected limit for Ly-alpha emission powered by Population II star formation and 12 sources among these do not have any continuum counterpart in a deep V-band imaging of the same field. For these, a stacking analysis indicates EW0>800A, effectively ruling out Ly-alpha powered by internal star formation. These sources are thus the best candidates so far for proto-galactic clouds or "dark" galaxies at high-redshift, whose existence has recently been suggested by several theoretical studies. Assuming they are mostly ionized by the quasar radiation, we estimate that their gas masses would be about 10^9 Msun implying that their star formation efficiencies (SFE) are less than 1.e-11 yr^-1 several times below the SFE of the most gas-rich dwarf galaxies locally, and five hundred times lower than typical massive star-forming galaxies at z~2. We have also discovered extended, filamentary gas, also likely illuminated by the quasar, around some of the brightest continuum-detected sources with EW0>240A. This emission is compatible with the expectations for circum-galactic cold streams but other origins, including tidal stripping, are also possible.

preprint2012arXiv

Mass and Environment as Drivers of Galaxy Evolution II: The quenching of satellite galaxies as the origin of environmental effects

We extend the phenomenological study of the evolving galaxy population of Peng et al (2010) to the central/satellite dichotomy in Yang et al. SDSS groups. We find that satellite galaxies are responsible for all the environmental effects in our earlier work. The fraction of centrals that are red does not depend on their environment but only on their stellar masses, whereas that of the satellites depends on both. We define a relative satellite quenching efficiency, which is the fraction of blue centrals that are quenched upon becoming the satellite of another galaxy. This is shown to be independent of stellar mass, but to depend strongly on local overdensity. The red fraction of satellites correlate much better with the local overdensity, a measure of location within the group, than with the richness of the group, i.e. dark matter halo mass. This, and the fact that satellite quenching depends on local density and not on either the stellar mass of the galaxy or the halo mass gives clues as to the nature of the satellite-quenching process. We furthermore show that the action of mass-quenching on satellite galaxies is also independent of the DM mass of the parent halo. We then apply the Peng et al (2010) approach to predict the mass functions of central and satellite galaxies, split into passive and active galaxies, and show that these match very well the observed mass functions from SDSS, further strengthening the validity of this phenomenological approach. We highlight the fact that the observed M* is the same for the star-forming centrals and satellites and the observed M* for the star-forming satellites is independent of halo mass above 10^12M\odot, which emphasizes the universality of the mass-quenching process that we identified in Peng et al (2010). Post-quenching merging modifies the mass function of the central galaxies but can increase the mass of typical centrals by only about 25%.

preprint2012arXiv

Shocked Superwinds from the z~2 Clumpy Star-forming Galaxy, ZC406690

We have obtained high-resolution data of the z 2 ring-like, clumpy star-forming galaxy (SFG) ZC406690 using the VLT/SINFONI with AO (in K-band) and in seeing-limited mode (in H- and J-band). Our data includes all of the main strong optical emission lines: [OII], [OIII], Ha, Hb, [NII] and [SII]. We find broad, blueshifted Ha and [OIII] emission line wings in the spectra of the galaxy's massive, star-forming clumps (sigma \sim 85 km s^-1) and even broader wings (up to 70% of the total Ha flux, with sigma \sim 290 km s^-1) in regions spatially offset from the clumps by \sim 2 kpc. The broad emission likely originates from large-scale outflows with mass outflow rates from individual clumps that are 1-8x the SFR of the clumps. Based on emission line ratio diagnostics ([NII]/Ha and [SII]/Ha) and photoionization and shock models, we find that the emission from the clumps is due to a combination of photoionization from the star-forming regions and shocks generated in the outflowing component, with 5-30% of the emission deriving from shocks. In terms of the ionization parameter (6x10^7-10^8 cm/s, based on both the SFR and the O32 ratio), density (local electron densities of 300-1800 cm^-3 in and around the clumps, and ionized gas column densities of 1200-8000 Msol/pc^2), and SFR (10-40 Msol/yr), these clumps more closely resemble nuclear starburst regions of local ULIRGs and dwarf irregulars than HII regions in local galaxies. However, the star-forming clumps are not located in the nucleus as in local starburst galaxies but instead are situated in a ring several kpc from the center of their high-redshift host galaxy, and have an overall disk-like morphology. The two brightest clumps are quite different in terms of their internal properties, energetics and relative ages, and thus we are given a glimpse at two different stages in the formation and evolution of rapidly star-forming giant clumps at high-z.

preprint2012arXiv

The SINS/zC-SINF survey of z~2 galaxy kinematics: Outflow properties

Based on SINFONI Ha, [NII] and [SII] AO data of 30 z \sim 2 star-forming galaxies (SFGs) from the SINS and zcSINF surveys, we find a strong correlation of the Ha broad flux fraction with the star formation surface density of the galaxy, with an apparent threshold for strong outflows occurring at 1 Msun yr^-1 kpc^-2. Above this threshold, we find that SFGs with logm_\ast>10 have similar or perhaps greater wind mass loading factors (eta = Mdotout/SFR) and faster outflow velocities than lower mass SFGs. This trend suggests that the majority of outflowing gas at z \sim 2 may derive from high-mass SFGs, and that the z \sim 2 mass-metallicity relation is driven more by dilution of enriched gas in the galaxy gas reservoir than by the efficiency of outflows. The mass loading factor is also correlated with the SFR and inclination, such that more star-forming and face-on galaxies launch more powerful outflows. For galaxies that have evidence for strong outflows, we find that the broad emission is spatially extended to at least the half-light radius (\sim a few kpc). We propose that the observed threshold for strong outflows and the observed mass loading of these winds can be explained by a simple model wherein break-out of winds is governed by pressure balance in the disk. Using the ratio of the [SII] doublet in a broad and narrow component, we find that outflowing gas has a density of \sim10-100 cm^-3, significantly less than that of the star forming gas (600 cm^-3).

preprint2010arXiv

The bulk of the black hole growth since z~1 occurs in a secular universe: No major merger-AGN connection

What is the relevance of major mergers and interactions as triggering mechanisms for active galactic nuclei (AGN) activity? To answer this longstanding question, we analyze 140 XMM-selected AGN host galaxies and a matched control sample of 1264 inactive galaxies over z~0.3-1.0 and log(M_*/M_sun)<11.7 with high-resolution HST/ACS imaging from the COSMOS field. The visual analysis of their morphologies by 10 independent human classifiers yields a measure of the fraction of distorted morphologies in the AGN and control samples, i.e. quantifying the signature of recent mergers which might potentially be responsible for fueling/triggering the AGN. We find that (1) the vast majority (>85%) of the AGN host galaxies do not show strong distortions, and (2) there is no significant difference in the distortion fractions between active and inactive galaxies. Our findings provide the best direct evidence that, since z~1, the bulk of black hole accretion has not been triggered by major galaxy mergers, therefore arguing that the alternative mechanisms, i.e., secular processes and minor interactions, are the leading triggers for the episodes of major black hole growth. We also exclude an alternative interpretation of our results: a significant time lag between merging and the observability of the AGN phase could wash out the most significant merging signatures, explaining the lack of enhancement of strong distortions on the AGN hosts. We show that this alternative scenario is unlikely due to: (1) recent major mergers being ruled out for the majority of sources due to the high fraction of disk-hosted AGN, (2) the lack of a significant X-ray signal in merging inactive galaxies as a signature of a potential buried AGN, and (3) the low levels of soft X-ray obscuration for AGN hosted by interacting galaxies, in contrast to model predictions.

preprint2009arXiv

Photo-z Performance for Precision Cosmology

Current and future weak lensing surveys will rely on photometrically estimated redshifts of very large numbers of galaxies. In this paper, we address several different aspects of the demanding photo-z performance that will be required for future experiments, such as the proposed ESA Euclid mission. It is first shown that the proposed all-sky near-infrared photometry from Euclid, in combination with anticipated ground-based photometry (e.g. PanStarrs-2 or DES) should yield the required precision in individual photo-z of sigma(z) < 0.05(1+z) at I_AB < 24.5. Simple a priori rejection schemes based on the photometry alone can be tuned to recognise objects with wildly discrepant photo-z and to reduce the outlier fraction to < 0.25% with only modest loss of otherwise usable objects. Turning to the more challenging problem of determining the mean redshift <z> of a set of galaxies to a precision of 0.002(1+z) we argue that, for many different reasons, this is best accomplished by relying on the photo-z themselves rather than on the direct measurement of <z> from spectroscopic redshifts of a representative subset of the galaxies. A simple adaptive scheme based on the statistical properties of the photo-z likelihood functions is shown to meet this stringent systematic requirement. We also examine the effect of an imprecise correction for Galactic extinction and the effects of contamination by fainter over-lapping objects in photo-z determination. The overall conclusion of this work is that the acquisition of photometrically estimated redshifts with the precision required for Euclid, or other similar experiments, will be challenging but possible. (abridged)

preprint2008arXiv

The COSMOS AGN Spectroscopic Survey I: XMM Counterparts

We present optical spectroscopy for an X-ray and optical flux-limited sample of 677 XMM-Newton selected targets covering the 2 deg^2 COSMOS field, with a yield of 485 high-confidence redshifts. The majority of the spectra were obtained over three seasons (2005-2007) with the IMACS instrument on the Magellan (Baade) telescope. We also include in the sample previously published Sloan Digital Sky Survey spectra and supplemental observations with MMT/Hectospec. We detail the observations and classification analyses. The survey is 90% complete to flux limits of f_{0.5-10 keV}>8 x 10^-16 erg cm^-2 s^-1 and i_AB+<22, where over 90% of targets have high-confidence redshifts. Making simple corrections for incompleteness due to redshift and spectral type allows for a description of the complete population to $i_AB+<23. The corrected sample includes 57% broad emission line (Type 1, unobscured) AGN at 0.13<z<4.26, 25% narrow emission line (Type 2, obscured) AGN at 0.07<z<1.29, and 18% absorption line (host-dominated, obscured) AGN at 0<z<1.22 (excluding the stars that made up 4% of the X-ray targets). We show that the survey's limits in X-ray and optical flux include nearly all X-ray AGN (defined by L_{0.5-10 keV}>3 x 10^42 erg s^-1) to z<1, of both optically obscured and unobscured types. We find statistically significant evidence that the obscured to unobscured AGN ratio at z<1 increases with redshift and decreases with luminosity.

preprint2007arXiv

Strong MgII systems in quasar and gamma-ray burst spectra

The incidence of strong MgII systems in gamma-ray burst (GRB) spectra is a few times higher than in quasar (QSO) spectra. We investigate several possible explanations for this effect, including: dust obscuration bias, clustering of the absorbers, different beam sizes of the sources, multiband magnification bias of GRBs, association of the absorbers with the GRB event or the circumburst environment. We find that: i) the incidence rate of MgII systems in QSO spectra could be underestimated by a factor 1.3-2 due to dust obscuration; ii) the equivalent-width distribution of the MgII absorbers along GRBs is consistent with that observed along QSOs thus suggesting that the absorbers are more extended than the beam sizes of the sources; iii) on average, GRB afterglows showing more than one MgII system are a factor of 1.7 brighter than the others, suggesting a lensing origin of the observed discrepancy; iv) gravitational lensing (in different forms, from galaxy lensing to microlensing) can bias high the counts of MgII systems along GRBs if the luminosity functions of the prompt gamma-ray emission and of the optical afterglows have a mean faint-end slope approaching -5/3 -- -2; v) some of the absorbers can be associated with the circumburst environment or produced by supernova remnants unrelated to the GRB event itself but lying in the same star-forming region. With the possible exception of magnification bias, it is unlikely that one of these effects on its own can fully account for the observed counts. However, the combined action of some of them can substantially reduce the statistical significance of the discrepancy.