Source author record

Michele Fumagalli

Michele Fumagalli 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

67works
5topics
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

67 published item(s)

preprint2024arXiv

ALMA-LEGUS I: The Influence of Galaxy Morphology on Molecular Cloud Properties

We present a comparative study of the molecular gas in two galaxies from the LEGUS sample: barred spiral NGC 1313 and flocculent spiral NGC 7793. These two galaxies have similar masses, metallicities, and star formation rates, but NGC 1313 is forming significantly more massive star clusters than NGC 7793, especially young massive clusters (<10 Myr, >10^4 Msol). Using ALMA CO(2-1) observations of the two galaxies with the same sensitivities and resolutions of 13 pc, we directly compare the molecular gas in these two similar galaxies to determine the physical conditions responsible for their large disparity in cluster formation. By fitting size-linewidth relations for the clouds in each galaxy, we find that NGC 1313 has a higher intercept than NGC 7793, implying that its clouds have higher kinetic energies at a given size scale. NGC 1313 also has more clouds near virial equilibrium than NGC 7793, which may be connected to its higher rate of massive cluster formation. However, these virially bound clouds do not show a stronger correlation with young clusters than that of the general cloud population. We find surprisingly small differences between the distributions of molecular cloud populations in the two galaxies, though the largest of those differences are that NGC 1313 has higher surface densities and lower free-fall times.

preprint2024arXiv

ALMA-LEGUS II: The Influence of Sub-Galactic Environment on Molecular Cloud Properties

We compare the molecular cloud properties in sub-galactic regions of two galaxies, barred spiral NGC 1313, which is forming many massive clusters, and flocculent spiral NGC 7793, which is forming significantly fewer massive clusters despite having a similar star formation rate to NGC 1313. We find that there are larger variations in cloud properties between different regions within each galaxy than there are between the galaxies on a global scale, especially for NGC 1313. There are higher masses, linewidths, pressures, and virial parameters in the arms of NGC 1313 and center of NGC 7793 than in the interarm and outer regions of the galaxies. The massive cluster formation of NGC 1313 may be driven by its greater variation in environments, allowing more clouds with the necessary conditions to arise, although no one parameter seems primarily responsible for the difference in star formation. Meanwhile NGC 7793 has clouds that are as massive and have as much kinetic energy as clouds in the arms of NGC 1313, but have densities and pressures more similar to the interarm regions and so are less inclined to collapse and form stars. The cloud properties in NGC 1313 and NGC 7793 suggest that spiral arms, bars, interarm regions, and flocculent spirals each represent distinct environments with regard to molecular cloud populations. We see surprisingly little difference in surface densities between the regions, suggesting that the differences in surface densities frequently seen between arm and interarm regions of lower-resolution studies are indicative of the sparsity of molecular clouds, rather than differences in their true surface density.

preprint2024arXiv

Probing coherence in metal absorption towards multiple images of strong gravitationally lensed quasars

We present a tomographic analysis of metal absorption lines arising from the circumgalactic medium (CGM) of galaxies at z~0.5-2, using Multi Unit Spectroscopic Explorer (MUSE) observations of two background quasars at z~2.2 and 2.8, which are two of the few currently known quasars with multiple images due to strong gravitational lensing by galaxy clusters at z~0.6 and 0.5, respectively. The angular separations between different pairs of quasar multiple images enable us to probe the absorption over transverse physical separations of ~0.4-150 kpc, which are based on strong lensing models exploiting MUSE observations. The fractional difference in rest-frame equivalent width (Delta Wr) of MgII, FeII, CIV absorption increases on average with physical separation, indicating that the metal-enriched gaseous structures become less coherent with distance, with a likely coherence length scale of ~10 kpc. However, Delta Wr for all the ions vary considerably over ~0.08-0.9, indicating a clumpy CGM over the full range of length scales probed. At the same time, paired MgII absorption is detected across ~100-150 kpc at similar line-of-sight velocities, which could be probing cool gas clouds within the same halo. No significant dependence of Delta Wr is found on the equivalent width and redshift of the absorbing gas and on the galaxy environment associated with the absorption. The high-ionization gas phase traced by CIV shows a higher degree of coherence than the low-ionization gas phase traced by MgII, with ~90 percent of CIV systems exhibiting Delta Wr <=0.5 at separations <=10 kpc compared to ~50 percent of MgII systems.

preprint2022arXiv

Constraining the Size of the Circumgalactic Medium Using the Transverse Autocorrelation Function of C IV Absorbers in Paired Quasar Spectra

The circumgalactic medium (CGM) plays a vital role in the formation and evolution of galaxies, acting as a lifeline between galaxies and the surrounding intergalactic medium (IGM). In this study we leverage a unique sample of quasar pairs to investigate the properties of the CGM with absorption line tomography. We present a new sample of medium resolution Keck/ESI, Magellan/MagE, and VLT/XSHOOTER spectra of 29 quasar pairs at redshift $2 < z < 3$. We supplement the sample with additional spectra of 32 pairs from the literature, creating a catalog of 61 quasar pairs with angular separations between 1.7" and 132.9" and projected physical separations ($r_\perp$) between 14 kpc and 887 kpc. We construct a catalog of 906 metal-line absorption doublets of C IV ($λλ1548, 1550$) with equivalent widths ranging from 6 mÅ $\leq W_{r, 1550} \leq 2053$ mÅ. The best fit linear model to the log-space equivalent width frequency distribution ($\log f(W_r) = m\log(W_{r}) + b$) of the sample yields coefficients of $m=-1.44\pm0.16$ and $b=-0.43\pm0.16$. To constrain the projected extent of C IV, we calculate the transverse autocorrelation function. The flattening of the autocorrelation function at low $r_\perp$ provides a lower limit for the coherence length of the metal enriched CGM - on the order of 200 $h^{-1}$ comoving kpc. This physical size constraint allows us to refine our understanding of the metals in the CGM, where the extent of C IV in the CGM depends on gas flows, feedback, timescale of metal injection and mixing, and the mass of the host galaxies.

preprint2022arXiv

Discovery of three new near-pristine absorption clouds at $z=2.6$-4.4

We report the discovery of three new "near-pristine" Lyman Limit Systems (LLSs), with metallicities ~1/1000 solar, at redshifts 2.6, 3.8 and 4.0, with a targeted survey at the Keck Observatory. High resolution echelle spectra of eight candidates yielded precise column densities of hydrogen and weak, but clearly detected, metal lines in seven LLSs; we previously reported the one remaining, apparently metal-free LLS, to have metallicity <1/10000 solar. Robust photoionisation modelling provides metallicities [Si/H] = -3.05 to -2.94, with 0.26 dex uncertainties (95% confidence) for three LLSs, and [Si/H] >~ -2.5 for the remaining four. Previous simulations suggest that near-pristine LLSs could be the remnants of PopIII supernovae, so comparing their detailed metal abundances with nucleosynthetic yields from supernovae models is an important goal. Unfortunately, at most two different metals were detected in each new system, despite their neutral hydrogen column densities (10^{19.2-19.4} cm^{-2}) being two orders of magnitude larger than the two previous, serendipitously discovered near-pristine LLSs. Nevertheless, the success of this first targeted survey for near-pristine systems demonstrates the prospect that a much larger, future survey could identify clear observational signatures of PopIII stars. With a well-understood selection function, such a survey would also yield the number density of near-pristine absorbers which, via comparison to future simulations, could reveal the origin(s) of these rare systems.

preprint2022arXiv

MUSE sneaks a peek at extreme ram-pressure stripping events -- V. Towards a complete view of the galaxy cluster A1367

We present an analysis of the kinematics and ionization conditions in a sample composed of seven star-forming galaxies undergoing ram-pressure stripping in the A1367 cluster, and the galaxy ESO137-001 in the Norma cluster. MUSE observations of two new galaxies in this sample, CGCG097-073 and CGCG097-079, are also presented. This sample is characterized by homogeneous integral field spectroscopy with MUSE and by a consistent selection based on the presence of ionised gas tails. The ratio [OI]/H$α$ is consistently elevated in the tails of these objects compared to what observed in unperturbed galaxy disks, an ubiquitous feature which we attribute to shocks or turbulent phenomena in the stripped gas. Compact star-forming regions are observed in only $\approx$ 50% of the tails, implying that specific (currently unknown) conditions are needed to trigger star formation inside the stripped gas. Focusing on the interface regions between the interstellar and intracluster medium, we observe different line ratios that we associate to different stages of the stripping process, with galaxies at an early stage of perturbation showing more prominent signatures of elevated star formation. Our analysis thus demonstrates the power of a well selected and homogeneous sample to infer general properties arising from ram-pressure stripping inside local clusters.

preprint2022arXiv

Oxygen-enhanced extremely metal-poor DLAs: A signpost of the first stars?

We present precise abundance determinations of two near-pristine damped Ly$α$ systems (DLAs) to assess the nature of the [O/Fe] ratio at [Fe/H] < -3 (i.e. <1/1000 of the solar metallicity). Prior observations indicate that the [O/Fe] ratio is consistent with a constant value, [O/Fe] ~ +0.4, when -3 < [Fe/H] < -2, but this ratio may increase when [Fe/H] < -3. In this paper, we test this picture by reporting new, high-precision [O/Fe] abundances in two of the most metal-poor DLAs currently known. We derive values of [O/Fe] = +0.50 +/- 0.10 and [O/Fe] = +0.62 +/- 0.05 for these two z ~ 3 near-pristine gas clouds. These results strengthen the idea that the [O/Fe] abundances of the most metal-poor DLAs are elevated compared to DLAs with [Fe/H] > -3. We compare the observed abundance pattern of the latter system to the nucleosynthetic yields of Population III supernovae (SNe), and find that the enrichment can be described by a (19-25) M$_{\odot}$ Population III SN that underwent a (0.9-2.4)$\times 10^{51}$ erg explosion. These high-precision measurements showcase the behaviour of [O/Fe] in the most metal-poor environments. Future high-precision measurements in new systems will contribute to a firm detection of the relationship between [O/Fe] and [Fe/H]. These data will reveal whether we are witnessing a chemical signature of enrichment from Population III stars and allow us to rule out contamination from Population II stars.

preprint2022arXiv

Probing the parameters of the intergalactic medium using quasars

We continue our series of papers on intergalactic medium (IGM) tracers using quasi-stellar objects (QSOs), having examined gamma-ray bursts (GRBs) and blazars in earlier studies. We have estimated the IGM properties of hydrogen column density (Nhxigm), temperature and metallicity using XMM-Newton QSO spectra over a large redshift range, with a collisional ionisation equilibrium (CIE) model for the ionised plasma. The Nhxigm parameter results were robust with respect to intrinsic power laws, spectral counts, reflection hump and soft excess features. There is scope for a luminosity bias given both luminosity and Nhxigm scale with redshift, but we find this unlikely given the consistent IGM parameter results across the other tracer types reviewed. The impact of intervening high column density absorbers was found to be minimal. The Nhxigm from the QSO sample scales as (1 + z)^1.5+/-0.2. The mean hydrogen density at z = 0 is n0 = (2.8 +/- 0.3) x 10^-7 cm^-3, the mean IGM temperature over the full redshift range is log(T/K) = 6.5+/-0.1, and the mean metallicity is [X/H] = -1.3+/-0.1(Z 0.05). Aggregating with our previous GRB and blazar tracers, we conclude that we have provided evidence of the IGM contributing substantially and consistently to the total X-ray absorption seen in the spectra. These results are based on the necessarily simplistic slab model used for the IGM, due to the inability of current X-ray data to constrain the IGM redshift distribution

preprint2021arXiv

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

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

preprint2021arXiv

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

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

preprint2021arXiv

Probing the physical properties of the intergalactic medium using gamma-ray bursts

We use Gamma-ray burst (GRB) spectra total continuum absorption to estimate the key intergalactic medium (IGM) properties of hydrogen column density ($\mathit{N}\textsc{hxigm}$), metallicity, temperature and ionisation parameter over a redshift range of $1.6 \leq z \leq 6.3$, using photo-ionisation (PIE) and collisional ionisation equilibrium (CIE) models for the ionised plasma. We use more realistic host metallicity, dust corrected where available, in generating the host absorption model, assuming that the host intrinsic hydrogen column density is equal to the measured ionisation corrected intrinsic neutral column from UV spectra ($\textit{N}\textsc{h}\/\ \textsc{i,ic}$). We find that the IGM property results are similar, regardless of whether the model assumes all PIE or CIE. The $\mathit{N}\textsc{hxigm}$ scales as $(1 + z)^{1.0\/\ -\/\ 1.9}$, with equivalent hydrogen mean density at $z = 0$ of $n_0 = 1.8^{+1.5}_{-1.2} \times 10^{-7}$ cm$^{-3}$. The metallicity ranges from $\sim0.1Z\sun$ at $z \sim 2$ to $\sim0.001Z\sun$ at redshift $z > 4$. The PIE model implies a less rapid decline in average metallicity with redshift compared to CIE. Under CIE, the temperature ranges between $5.0 <$ log$(T/$K$)<\/\ 7.1$. For PIE the ionisation parameter ranges between $0.1 <$ log$(ξ) < 2.9$. Using our model, we conclude that the IGM contributes substantially to the total absorption seen in GRB spectra and that this contribution rises with redshift, explaining why the hydrogen column density inferred from X-rays is substantially in excess of the intrinsic host contribution measured in UV.

preprint2021arXiv

Sub-damped Lyman alpha systems in the XQ-100 survey II -- Chemical evolution at 2.4<z<4.3

We present the measured gas-phase metal column densities in 155 sub-damped Lyman alpha systems (subDLAs) with the aim to investigate the contribution of subDLAs to the chemical evolution of the Universe. The sample was identified within the absorber-blind XQ-100 quasar spectroscopic survey over the redshift range 2.4<=z<=4.3. Using all available column densities of the ionic species investigated (mainly CIV, SiII, MgII, SiIV, AlII, FeII, CII, and OI; in order of decreasing detection frequency), we estimate the ionization-corrected gas-phase metallicity of each system using Markov Chain Monte Carlo techniques to explore a large grid of Cloudy ionization models. Without accounting for ionization and dust depletion effects, we find that the HI-weighted gas-phase metallicity evolution of subDLAs are consistent with damped Lyman alpha systems (DLAs). When ionization corrections are included, subDLAs are systematically more metal-poor than DLAs (between ~0.5 sigma and ~3 sigma significance) by up to ~1.0 dex over the redshift range 3<=z<=4.3. The correlation of gas-phase [Si/Fe] with metallicity in subDLAs appears to be consistent with that of DLAs, suggesting that the two classes of absorbers have a similar relative dust depletion pattern. As previously seen for Lyman limit systems, the gas-phase [C/O] in subDLAs remains constantly solar for all metallicities indicating that both subDLAs and Lyman limit systems could trace carbon-rich ejecta, potentially in circumgalactic environments.

preprint2020arXiv

A bound on the 12C/13C ratio in near-pristine gas with ESPRESSO

Using science verification observations obtained with ESPRESSO at the Very Large Telescope (VLT) in 4UT mode, we report the first bound on the carbon isotope ratio 12C/13C of a quiescent, near-pristine damped Ly-alpha (DLA) system at z=2.34. We recover a limit log10(12C/13C) > +0.37 (2 sigma). We use the abundance pattern of this DLA, combined with a stochastic chemical enrichment model, to infer the properties of the enriching stars, finding the total gas mass of this system to be log10(M_gas/M_sun)=6.3+1.4-0.9 and the total stellar mass to be log10(M_*/M_sun)=4.8+/-1.3. The current observations disfavour enrichment by metal-poor Asymptotic Giant Branch (AGB) stars with masses <2.4 Msun, limiting the epoch at which this DLA formed most of its enriching stars. Our modelling suggests that this DLA formed very few stars until >1 Gyr after the cosmic reionization of hydrogen and, despite its very low metallicity (~1/1000 of solar), this DLA appears to have formed most of its stars in the past few hundred Myr. Combining the inferred star formation history with evidence that some of the most metal-poor DLAs display an elevated [C/O] ratio at redshift z<3, we suggest that very metal-poor DLAs may have been affected by reionization quenching. Finally, given the simplicity and quiescence of the absorption features associated with the DLA studied here, we use these ESPRESSO data to place a bound on the possible variability of the fine-structure constant, Delta alpha/alpha=(-1.2 +/- 1.1)x10^-5.

preprint2020arXiv

A limit on Planck-scale froth with ESPRESSO

Some models of quantum gravity predict that the very structure of spacetime is `frothy' due to quantum fluctuations. Although the effect is expected to be tiny, if these spacetime fluctuations grow over a large distance, the initial state of a photon, such as its energy, is gradually washed out as the photon propagates. Thus, in these models, even the most monochromatic light source would gradually disperse in energy due to spacetime fluctuations over large distances. In this paper, we use science verification observations obtained with ESPRESSO at the Very Large Telescope to place a novel bound on the growth of spacetime fluctuations. To achieve this, we directly measure the width of a narrow Fe II absorption line produced by a quiescent gas cloud at redshift z=2.34, corresponding to a comoving distance of ~5.8 Gpc. Using a heuristic model where the energy fluctuations grow as sigma_E / E = (E/E_P)^alpha, where E_P = 1.22 x 10^28 eV is the Planck energy, we rule out models with alpha < 0.634, including models where the quantum fluctuations grow as a random walk process (alpha = 0.5). Finally, we present a new formalism where the uncertainty accrued at discrete spacetime steps is drawn from a continuous distribution. We conclude, if photons take discrete steps through spacetime and accumulate Planck-scale uncertainties at each step, then our ESPRESSO observations require that the step size must be at least >10^13.2 L_P, where L_P is the Planck length.

preprint2020arXiv

Detecting neutral hydrogen at z > 3 in large spectroscopic surveys of quasars

We present a pipeline based on a random forest classifier for the identification of high column-density clouds of neutral hydrogen (i.e. the Lyman limit systems, LLSs) in absorption within large spectroscopic surveys of z>3 quasars. We test the performance of this method on mock quasar spectra that reproduce the expected data quality of the Dark Energy Spectroscopic Instrument (DESI) and the WHT Enhanced Area Velocity Explorer (WEAVE) surveys, finding >90% completeness and purity for N(HI)> 10^17.2 cm^-2 LLSs against quasars of g < 23 mag at z~3.5-3.7. After training and applying our method on 10,000 quasar spectra at z~3.5-4.0 from the Sloan Digital Sky Survey (Data Release 16), we identify ~6600 LLSs with N(HI)>10^17.5 cm^-2 between z~3.1-4.0 with a completeness and purity of >90% for the classification of LLSs. Using this sample, we measure a number of LLSs per unit redshift of 2.32 +/- 0.08 at z=[3.3,3.6]. We also present results on the performance of random forest for the measurement of the LLS redshifts and HI column densities, and for the identification of broad absorption line quasars.

preprint2020arXiv

Discovery of a damped Ly$α$ galaxy at z $\sim$ 3 towards the quasar SDSS J011852+040644

We report the detection of the host galaxy of a damped Ly$α$ system (DLA) with log N(HI) $ [\rm cm^{-2}]$ = $21.0 \pm 0.10$ at $z \approx 3.0091$ towards the background quasar SDSS J011852+040644 using the Palomar Cosmic Web Imager (PCWI) at the Hale (P200) telescope. We detect Ly$α$ emission in the dark core of the DLA trough at a 3.3$σ$ confidence level, with Ly$α$ luminosity of $L_{\rm Lyα}$ $\rm = (3.8 \pm 0.8) \times 10^{42}\ erg\ s^{-1}$, corresponding to a star formation rate of $\gtrsim 2\ \rm M_{\odot}\ yr^{-1}$ (considering a lower limit on Ly$α$ escape fraction $f_{esc}^{Lyα} \sim 2\%$) as typical for Lyman break galaxies at these redshifts. The Ly$α$ emission is blueshifted with respect to the systemic redshift derived from metal absorption lines by $281 \pm 43$ km/s. The associated galaxy is at very small impact parameter of $\lesssim 12 \rm\ kpc$ from the background quasar, which is in line with the observed anticorrelation between column density and impact parameter in spectroscopic searches tracing the large-scale environments of DLA host galaxies.

preprint2020arXiv

Interpreting Observations of Absorption Lines in the Circumgalactic Medium with a Turbulent Medium

Single-phase photoionization equilibrium (PIE) models are often used to infer the underlying physical properties of galaxy halos probed in absorption with ions at different ionization potentials. To incorporate the effects of turbulence, we use the MAIHEM code to model an isotropic turbulent medium exposed to a redshift zero metagalactic UV background, while tracking the ionizations, recombinations, and species-by-species radiative cooling for a wide range of ions. By comparing observations and simulations over a wide range of turbulent velocities, densities, and metallicity with a Markov chain Monte Carlo technique, we find that MAIHEM models provide an equally good fit to the observed low-ionization species compared to PIE models, while reproducing at the same time high-ionization species such as \ion{Si}{4} and \ion{O}{6}. By including multiple phases, MAIHEM models favor a higher metallicity ($Z/Z_\odot \approx 40\%$) for the circumgalactic medium compared to PIE models. Furthermore, all of the solutions require some amount of turbulence ($σ_{\rm 3D} \geqslant 26\ {\rm km}\ {\rm s}^{-1}$). Correlations between turbulence, metallicity, column density, and impact parameter are discussed alongside mechanisms that drive turbulence within the halo.

preprint2020arXiv

Into the Lyα jungle: exploring the circumgalactic medium of galaxies at z ~ 4-5 with MUSE

We present a study of the galaxy environment of 9 strong HI+CIV absorption line systems ($16.2<{\rm log}(N({\rm HI}))<21.2$) spanning a wide range in metallicity at $z\sim4-5$, using MUSE integral field and X-Shooter spectroscopic data collected in a $z\approx 5.26$ quasar field. We identify galaxies within a 250 kpc and $\pm1000$ km s$^{-1}$ window for 6 out of the 9 absorption systems, with 2 of the absorption line systems showing multiple associated galaxies within the MUSE field of view. The space density of Ly$α$ emitting galaxies (LAEs) around the HI and CIV systems is $\approx10-20$ times the average sky density of LAEs given the flux limit of our survey, showing a clear correlation between the absorption and galaxy populations. Further, we find that the strongest CIV systems in our sample are those that are most closely aligned with galaxies in velocity space, i.e. within velocities of $\pm500$ km s$^{-1}$. The two most metal poor systems lie in the most dense galaxy environments, implying we are potentially tracing gas that is infalling for the first time into star-forming groups at high redshift. Finally, we detect an extended Ly$α$ nebula around the $z\approx 5.26$ quasar, which extends up to $\approx50$ kpc at the surface brightness limit of $3.8 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. After scaling for surface brightness dimming, we find that this nebula is centrally brighter, having a steeper radial profile than the average for nebulae studied at $z\sim3$ and is consistent with the mild redshift evolution seen from $z\approx 2$.

preprint2019arXiv

Determining the primordial helium abundance and UV background using fluorescent emission in star-free dark matter haloes

Observational measures of the primordial helium mass fraction, $Y_{\rm P}$, are of interest for cosmology and fundamental particle physics. Current measures obtained from H II regions agree with the Standard Model prediction to approximately 1% precision, although these determinations may be affected by systematic uncertainties. This possibility can only be tested by independently measuring the helium abundance in new ways. Here, we propose a novel method to obtain a measurement of $Y_{\rm P}$ using hydrogen and helium recombination line emission from RELHICs: pristine, gas-rich but star-free low-mass dark matter haloes whose existence is predicted by hydrodynamical simulations. Although expected to be uncommon and intrinsically faint in emission, the primordial composition and simple physical properties of these objects make them an ideal laboratory to determine $Y_{\rm P}$. We present radiative transfer simulations to demonstrate the effectiveness of this approach, finding that comparing the emission in H and He lines, either via their volumetric emissivities, or integrated properties such as the surface brightness and total flux, may be used to infer $Y_{\rm P}$. Furthermore, we show that RELHICs can be used to provide an entirely novel constraint on the spectral slope of the ultraviolet background, and discuss the possibility of measuring this slope and the primordial helium abundance simultaneously.

preprint2019arXiv

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

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

preprint2016arXiv

Exploring Damped Lyman-$α$ System Host Galaxies using Gamma-ray Bursts

We present a sample of 45 Damped Lyman-$α$ system (DLA, \nhi\ $ \geq 2 \times 10^{20} {\rm cm}^{-2}$) counterparts (33 detections, 12 upper limits) which host gamma-ray bursts (GRB-DLAs) in order to investigate star-formation and metallicity within galaxies hosting DLAs. Our sample spans $z \sim 2-6$ and is nearly three times larger than any previously detected DLA counterparts survey based on quasar line-of-sight searches (QSO-DLAs). We report star formation rates (SFRs) from rest-frame UV photometry and SED modeling. We find that DLA counterpart SFRs are not correlated with either redshift or HI column density. Thanks to the combination of \hst\ and ground-based observations, we also investigate DLA host star-formation efficiency. Our GRB-DLA counterpart sample spans both higher efficiency and low efficiency star formation regions compared to the local Kennicutt-Schmidt relation, local star formation laws, and $z \sim 3$ cosmological simulations. We also compare the depletion times of our DLA hosts sample to other objects in the local Universe, our sample appears to deviate from the star formation efficiencies measured in local spiral and dwarf galaxies. Furthermore, we find similar efficiencies as local inner disks, SMC, and LBG outskirts. Finally, our enrichment time measurements show a spread of systems with under- and over-abundance of metals which may suggest that these systems had episodic star formation and a metal enrichment/depletion as a result of strong stellar feedback and/or metal inflow/outflow.

preprint2016arXiv

MUSE searches for galaxies near very metal-poor gas clouds at z~3: new constraints for cold accretion models

We report on the search for galaxies in the proximity of two very metal-poor gas clouds at z~3 towards the quasar Q0956+122. With a 5-hour MUSE integration in a ~500x500 kpc^2 region centred at the quasar position, we achieve a >80% complete spectroscopic survey of continuum-detected galaxies with m<25 mag and Lyα emitters with luminosity L>3e41 erg/s. We do not identify galaxies at the redshift of a z~3.2 Lyman limit system (LLS) with log Z/Zsun = -3.35 +/- 0.05, placing this gas cloud in the intergalactic medium or circumgalactic medium of a galaxy below our sensitivity limits. Conversely, we detect five Lyα emitters at the redshift of a pristine z~3.1 LLS with log Z/Zsun < -3.8, while ~0.4 sources were expected given the z~3 Lyα luminosity function. Both this high detection rate and the fact that at least three emitters appear aligned in projection with the LLS suggest that this pristine cloud is tracing a gas filament that is feeding one or multiple galaxies. Our observations uncover two different environments for metal-poor LLSs, implying a complex link between these absorbers and galaxy halos, which ongoing MUSE surveys will soon explore in detail. Moreover, in agreement with recent MUSE observations, we detected a ~90 kpc Lyα nebula at the quasar redshift and three Lyα emitters reminiscent of a "dark galaxy" population.

preprint2016arXiv

MUSE sneaks a peek at extreme ram-pressure stripping events. II. The physical properties of the gas tail of ESO137-001

We present a study of the physical properties of the disc and tail of ESO137-001, a galaxy suffering from extreme ram-pressure stripping during its infall into the Norma cluster. With sensitive and spatially-resolved MUSE spectroscopy, we analyse the emission line diagnostics in the tail of ESO137-001, finding high values of [NII]/H$α$ and [OI]/H$α$ that are suggestive of the presence of shocks in turbulent gas. However, the observed line ratios are not as strong as commonly seen in pure shock heating models, suggesting that other emission mechanisms may contribute to the observed emission. Indeed, part of the observed emission, particularly at close separations from the galaxy disc, may originate from recombination of photoionised gas stripped from the main body of ESO137-001. We also identify a large number of bright compact knots within in the tail, with line ratios characteristic of HII regions. These HII regions, despite residing in a stripped gas tail, have quite typical line ratios, densities, temperatures, and metallicity ($\sim0.7$ solar). The majority of these HII regions are embedded within diffuse gas from the tail, which is dynamically cool ($σ\sim 25-50\ \rm{km\ s^{-1}}$ ). This fact, together with a lack of appreciable gradients in age and metallicity, suggests that these HII regions formed in situ. While our analysis represents a first attempt to characterise the rich physics of the ESO137-001 tail, future work is needed to address the importance of other mechanisms, such as thermal conduction and magneto hydrodynamic waves, in powering the emission in the tail.

preprint2016arXiv

On the connection between the metal-enriched intergalactic medium and galaxies: an OVI-galaxy cross-correlation study at $z < 1$

We present new results on the auto- and cross-correlation functions of galaxies and OVI absorbers in a $\sim 18~\textrm{Gpc}^3$ comoving volume at $z < 1$. We use a sample of 51,296 galaxies and 140 OVI absorbers in the column density range $13 \lesssim \log N \lesssim 15$ to measure two-point correlation functions in the two dimensions transverse and orthogonal to the line-of-sight $ξ(r_{\perp}, r_{\parallel})$. We furthermore infer the corresponding 'real-space' correlation functions, $ξ(r)$, by projecting $ξ(r_{\perp}, r_{\parallel})$ along $r_{\parallel}$, and assuming a power-law form, $ξ(r) = (r / r_0)^{-γ}$. Comparing the results from the absorber-galaxy cross-correlation function, $ξ_{\textrm{ag}}$, the galaxy auto-correlation function, $ξ_{\textrm{gg}}$, and the absorber auto-correlation function, $ξ_{\textrm{aa}}$, we constrain the statistical connection between galaxies and the metal-enriched intergalactic medium as a function of star-formation activity. We also compare these results to predictions from the EAGLE cosmological hydrodynamical simulation and find a reasonable agreement. We find that: (i) OVI absorbers show very little velocity dispersion with respect to galaxies on $\sim$ Mpc scales, likely $\lesssim$ 100 \kms; (ii) OVI absorbers and galaxies may not linearly trace the same underlying distribution of matter in general. In particular, our results demonstrate that OVI absorbers are less clustered, and potentially more extended around galaxies than galaxies are around themselves; (iii) On $\gtrsim 100$ kpc scales, the likelihood of finding OVI absorbers around star-forming galaxies is similar to the likelihood of finding OVI absorbers around non star-forming galaxies (abridged)

preprint2016arXiv

Robust automatic photometry of local galaxies from SDSS Dissecting the color magnitude relation with color profiles

We present an automatic procedure to perform reliable photometry of galaxies on SDSS images. We selected a sample of 5853 galaxies in the Coma and Virgo superclusters. For each galaxy, we derive Petrosian g and i magnitudes, surface brightness profiles and color profiles. Unlike the SDSS pipeline, our procedure is not affected by the well known shredding problem and efficiently extracts Petrosian magnitudes for all galaxies. Hence we derived magnitudes even from the population of galaxies missed by the SDSS that represents 25% of all Local supercluster galaxies and ~95% of galaxies with g < 11 mag. After correcting the g and i magnitudes for Galactic and internal extinction, the blue and red sequences in the color magnitude diagram are well separated, with similar slopes. In addition, we study (i) the color-magnitude diagrams in different galaxy regions, the inner (r <= 1 kpc), intermediate (0.2R{Pet} <=r<= 0.3R{Pet}) and outer, disk-dominated (r>=0.35R{Pet})) zone, and (ii), we compute template color profiles, discussing the dependences of the templates on the galaxy masses and on their morphological type. The two analyses consistently lead to a picture where elliptical galaxies show no color gradients, irrespective of their masses. Spirals, instead, display a steeper gradient in their color profiles with increasing mass, consistently with the growing relevance of a bulge and/or a bar component above 10^{10} Modot.

preprint2016arXiv

The Cosmic Evolution of the Metallicity Distribution of Ionized Gas Traced by Lyman Limit Systems

We present the first results from our KODIAQ Z survey aimed to determine the metallicity distribution and physical properties of the z>2 partial and full Lyman limit systems (pLLSs and LLSs; 16.2<log N(HI)<19), which are probed of the interface regions between the intergalactic medium (IGM) and galaxies. We study 31 HI-selected pLLSs and LLSs at 2.3<z<3.3 observed with Keck/HIRES in absorption against background QSOs. We compare the column densities of metal-ions to HI and use photoionization models to assess the metallicity. The metallicity distribution of the pLLSs/LLSs at 2.3<z<3.3 is consistent with a unimodal distribution peaking at [X/H]=-2. The metallicity distribution of these absorbers therefore evolves markedly with z since at z<1 it is bimodal with peaks at [X/H]=-1.8 and -0.3. There is a substantial fraction (25-41%) of pLLSs/LLSs with metallicities well below those of damped Lya absorbers (DLAs) at any studied z from z<1 to z~2-4, implying reservoirs of metal-poor cool, dense gas in the IGM/galaxy interface at all z. However, the gas probed by pLLSs and LLSs is rarely pristine, with a fraction 3-18% for pLLSs/LLSs with [X/H]<-3. We find C/alpha enhancement in several pLLSs and LLSs in the metallicity range -2<[X/H]<-0.5, where C/alpha is 2-5 times larger than observed in Galactic metal-poor stars or high redshift DLAs at similar metallicities. This is likely caused by preferential ejection of carbon from metal-poor galaxies into their surroundings.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2016arXiv

The Star Formation Rate Efficiency of Neutral Atomic-dominated Hydrogen Gas in the Outskirts of Star Forming Galaxies from z~1 to z~3

Current observational evidence suggests that the star formation rate (SFR) efficiency of neutral atomic hydrogen gas measured in Damped Ly-alpha Systems (DLAs) at z~3 is more than 10 times lower than predicted by the Kennicutt-Schmidt (KS) relation. To understand the origin of this deficit, and to investigate possible evolution with redshift and galaxy properties, we measure the SFR efficiency of atomic gas at z~1, z~2, and z~3 around star-forming galaxies. We use new robust photometric redshifts in the Hubble Ultra Deep Field to create galaxy stacks in these three redshift bins, and measure the SFR efficiency by combining DLA absorber statistics with the observed rest-frame UV emission in the galaxies' outskirts. We find that the SFR efficiency of HI gas at z>1 is ~1-3% of that predicted by the KS relation. Contrary to simulations and models that predict a reduced SFR efficiency with decreasing metallicity and thus with increasing redshift, we find no significant evolution in the SFR efficiency with redshift. Our analysis instead suggests that the reduced SFR efficiency is driven by the low molecular content of this atomic-dominated phase, with metallicity playing a secondary effect in regulating the conversion between atomic and molecular gas. This interpretation is supported by the similarity between the observed SFR efficiency and that observed in local atomic-dominated gas, such as in the outskirts of local spiral galaxies and local dwarf galaxies.

preprint2015arXiv

Halpha3: an Halpha imaging survey of HI selected galaxies from ALFALFA. V: The Coma Supercluster survey completion

Neutral hydrogen represents the major observable baryonic constituent of galaxies that fuels the formation of stars through the transformation in molecular hydrogen. The emission of the hydrogen recombination line Halpha is the most direct tracer of the process that transforms gas (fuel) into stars. We continue to present Halpha3 (acronym for Halpha-alpha-alpha), an extensive Halpha+[NII] narrow-band imaging campaign of galaxies selected from the HI Arecibo Legacy Fast ALFA Survey (ALFALFA), using the instrumentation available at the San Pedro Martir observatory (Mexico). In only four years since 2011 we were able to complete in 48 nights the Halpha imaging observations of 724 galaxies in the region of the Coma supercluster 10^h < R.A. <16^h; 24^o < Dec. <28^o and 3900<cz<9000 kms^{-1}. Of these, 603 are selected from the HI Arecibo Legacy Fast ALFA Survey (ALFALFA) and constitute a 97% complete sample. They provide for the first time a complete census of the massive star formation properties of local gas-rich galaxies belonging to different environments (cluster vs filaments), morphological type (spirals vs dwarf Irr), over a wide range of stellar mass (10^{8}-10^{11.5} Modot) in the Coma Supercluster. The present Paper V provides the Halpha data and the derived star formation rates for the observed galaxies.

preprint2015arXiv

Metal-enriched, sub-kiloparsec gas clumps in the circumgalactic medium of a faint z = 2.5 galaxy

We report the serendipitous detection of a 0.2 L$^*$, Lyman-$α$ emitting galaxy at redshift 2.5 at an impact parameter of 50 kpc from a bright background QSO sightline. A high-resolution spectrum of the QSO reveals a partial Lyman-limit absorption system ($N_\mathrm{HI}=10^{16.94\pm0.10}$ cm$^{-2}$) with many associated metal absorption lines at the same redshift as the foreground galaxy. Using photoionization models that carefully treat measurement errors and marginalise over uncertainties in the shape and normalisation of the ionizing radiation spectrum, we derive the total hydrogen column density $N_\mathrm{H}=10^{19.4\pm0.3}$ cm$^{-2}$, and show that all the absorbing clouds are metal enriched, with $Z=0.1$-$0.6 Z_\odot$. These metallicities and the system's large velocity width ($436$ km$\,$s$^{-1}$) suggest the gas is produced by an outflowing wind. Using an expanding shell model we estimate a mass outflow rate of $\sim5 M_\odot\,$yr$^{-1}$. Our photoionization model yields extremely small sizes ($<$100-500 pc) for the absorbing clouds, which we argue are typical of high column density absorbers in the circumgalactic medium (CGM). Given these small sizes and extreme kinematics, it is unclear how the clumps survive in the CGM without being destroyed by hydrodynamic instabilities. The small cloud sizes imply that even state-of-the-art cosmological simulations require more than a $1000$-fold improvement in mass resolution to resolve the hydrodynamics relevant for cool gas in the CGM.

preprint2015arXiv

SLUG -- Stochastically Lighting Up Galaxies. III: A Suite of Tools for Simulated Photometry, Spectroscopy, and Bayesian Inference with Stochastic Stellar Populations

Stellar population synthesis techniques for predicting the observable light emitted by a stellar population have extensive applications in numerous areas of astronomy. However, accurate predictions for small populations of young stars, such as those found in individual star clusters, star-forming dwarf galaxies, and small segments of spiral galaxies, require that the population be treated stochastically. Conversely, accurate deductions of the properties of such objects also requires consideration of stochasticity. Here we describe a comprehensive suite of modular, open-source software tools for tackling these related problems. These include: a greatly-enhanced version of the slug code introduced by da Silva et al. (2012), which computes spectra and photometry for stochastically- or deterministically-sampled stellar populations with nearly-arbitrary star formation histories, clustering properties, and initial mass functions; cloudy_slug, a tool that automatically couples slug-computed spectra with the cloudy radiative transfer code in order to predict stochastic nebular emission; bayesphot, a general-purpose tool for performing Bayesian inference on the physical properties of stellar systems based on unresolved photometry; and cluster_slug and sfr_slug, a pair of tools that use bayesphot on a library of slug models to compute the mass, age, and extinction of mono-age star clusters, and the star formation rate of galaxies, respectively. The latter two tools make use of an extensive library of pre-computed stellar population models, which are included the software. The complete package is available at http://www.slugsps.com.

preprint2015arXiv

Star Cluster Properties in Two LEGUS Galaxies Computed with Stochastic Stellar Population Synthesis Models

We investigate a novel Bayesian analysis method, based on the Stochastically Lighting Up Galaxies (slug) code, to derive the masses, ages, and extinctions of star clusters from integrated light photometry. Unlike many analysis methods, slug correctly accounts for incomplete IMF sampling, and returns full posterior probability distributions rather than simply probability maxima. We apply our technique to 621 visually-confirmed clusters in two nearby galaxies, NGC 628 and NGC 7793, that are part of the Legacy Extragalactic UV Survey (LEGUS). LEGUS provides Hubble Space Telescope photometry in the NUV, U, B, V, and I bands. We analyze the sensitivity of the derived cluster properties to choices of prior probability distribution, evolutionary tracks, IMF, metallicity, treatment of nebular emission, and extinction curve. We find that slug's results for individual clusters are insensitive to most of these choices, but that the posterior probability distributions we derive are often quite broad, and sometimes multi-peaked and quite sensitive to the choice of priors. In contrast, the properties of the cluster population as a whole are relatively robust against all of these choices. We also compare our results from slug to those derived with a conventional non-stochastic fitting code, Yggdrasil. We show that slug's stochastic models are generally a better fit to the observations than the deterministic ones used by Yggdrasil. However, the overall properties of the cluster populations recovered by both codes are qualitatively similar.

preprint2015arXiv

The Keck+Magellan Survey for Lyman Limit Absorption III: Sample Definition and Column Density Measurements

We present an absorption-line survey of optically thick gas clouds -- Lyman Limit Systems (LLSs) -- observed at high dispersion with spectrometers on the Keck and Magellan telescopes. We measure column densities of neutral hydrogen NHI and associated metal-line transitions for 157 LLSs at z=1.76-4.39 restricted to 10^17.3 < NHI < 10^20.3. An empirical analysis of ionic ratios indicates an increasing ionization state of the gas with decreasing NHI and that the majority of LLSs are highly ionized, confirming previous expectations. The Si^+/H^0 ratio spans nearly four orders-of-magnitude, implying a large dispersion in the gas metallicity. Fewer than 5% of these LLSs have no positive detection of a metal transition; by z~3, nearly all gas that is dense enough to exhibit a very high Lyman limit opacity has previously been polluted by heavy elements. We add new measurements to the small subset of LLS (~5-10) that may have super-solar abundances. High Si^+/Fe^+ ratios suggest an alpha-enhanced medium whereas the Si^+/C^+ ratios do not exhibit the super-solar enhancement inferred previously for the Lya forest.

preprint2015arXiv

The Neutral Hydrogen Cosmological Mass Density at z=5

We present the largest homogeneous survey of $z>4.4$ damped Lyman-$α$ systems (DLAs) using the spectra of 163 QSOs that comprise the Giant Gemini GMOS (GGG) survey. With this survey we make the most precise high-redshift measurement of the cosmological mass density of neutral hydrogen, $Ω_{\rm HI}$. At such high redshift important systematic uncertainties in the identification of DLAs are produced by strong intergalactic medium absorption and QSO continuum placement. These can cause spurious DLA detections, result in real DLAs being missed, or bias the inferred DLA column density distribution. We correct for these effects using a combination of mock and higher-resolution spectra, and show that for the GGG DLA sample the uncertainties introduced are smaller than the statistical errors on $Ω_{\rm HI}$. We find $Ω_{\rm HI}=0.98^{+0.20}_{-0.18}\times10^{-3}$ at $\langle z\rangle=4.9$, assuming a 20% contribution from lower column density systems below the DLA threshold. By comparing to literature measurements at lower redshifts, we show that $Ω_{\rm HI}$ can be described by the functional form $Ω_{\rm HI}(z)\propto(1+z)^{0.4}$. This gradual decrease from $z=5$ to $0$ is consistent with the bulk of HI gas being a transitory phase fuelling star formation, which is continually replenished by more highly-ionized gas from the intergalactic medium, and from recycled galactic winds.

preprint2015arXiv

The physical properties of z>2 Lyman limit systems: new constraints for feedback and accretion models

We study the physical properties of a homogeneous sample of 157 optically-thick absorption line systems at redshifts ~1.8-4.4, selected from a high-dispersion spectroscopic survey of Lyman limit systems (LLSs). By means of multiple ionisation models and Bayesian techniques, we derive the posterior probability distribution functions for the density, metallicity, temperature, and dust content of the absorbing gas. We find that z>2 LLSs are highly ionised with ionisation parameters between -3<log U<-2, depending on the HI column density. LLSs are characterised by low temperatures (T<5x10^4 K) and reside in dust-poor environments. Between z~2.5-3.5, ~80% of the LLSs have physical densities between n(H)~10^-3.5-10^-2 cm^-3 for the assumed UV background, but we caution that a degeneracy between the ionisation parameter and the intensity of the radiation field prevents robust inference on the density and sizes of LLSs. Conversely, metallicity estimates are less sensitive to the assumptions behind ionisation corrections. LLSs at z>2 are characterised by a broad unimodal distribution over >4 orders of magnitude, with a peak at log Z/Zsun~-2. LLSs are metal poor, significantly less enriched than DLAs, with ~70% of the metallicity PDF below log Z/Zsun<-1.5. The median metallicity of super LLSs with log N(HI)>19 rapidly evolves with redshift, with a ten-fold increase between z~2.1-3.6 (~1.5 Gyr). Based on this sample, we find that LLSs at z=2.5-3.5 account for ~15% of all the metals produced by UV-selected galaxies. The implications for theories of cold gas accretion and metal ejection from galaxies are also discussed.

preprint2014arXiv

A compact, metal-rich, kpc-scale outflow in FBQS J0209-0438: Detailed diagnostics from HST/COS extreme UV observations

We present HST/COS observations of highly ionized absorption lines associated with a radio-loud QSO at $z=1.1319$. The absorption system has multiple velocity components, tracing gas that is largely outflowing from the QSO at velocities of a few 100 km s$^{-1}$. There is an unprecedented range in ionization, with detections of HI, NIII, NIV, NV, OIV, OIV*, OV, OVI, NeVIII, MgX, SV and ArVIII. We estimate the total hydrogen number density from the column density ratio N(OIV*)/N(OIV) to be $\log(n_{\textrm{H}}/\textrm{cm}^3)\sim 3$. Assuming photoionization equilibrium, we derive a distance to the absorbing complex of $2.3<R<6.0$ kpc from the centre of the QSO. A range in ionization parameter, covering $\sim 2$ orders of magnitude, suggest absorption path lengths in the range $10^{-4.5}<l_{\textrm{abs}}<1$ pc. In addition, the absorbing gas only partially covers the background emission from the QSO continuum, which suggests clouds with transverse sizes $l_{\textrm{trans}}<10^{-2.5}$ pc. Widely differing absorption path lengths, combined with covering fractions less than unity across all ions pose a challenge to models involving simple cloud geometries. These issues may be mitigated by the presence of non-equilibrium effects, together with the possibility of multiple gas temperatures. The dynamics and expected lifetimes of the gas clouds suggest that they do not originate from close to the AGN, but are instead formed close to their observed location. Their inferred distance, outflow velocities and gas densities are broadly consistent with scenarios involving gas entrainment or condensations in winds driven by either supernovae, or the supermassive black hole accretion disc. In the case of the latter, the present data most likely does not trace the bulk of the outflow by mass, which could instead manifest itself as an accompanying warm absorber, detectable in X-rays.

preprint2014arXiv

Directly imaging damped Ly-alpha galaxies at z>2. II: Imaging and spectroscopic observations of 32 quasar fields

Damped Ly-alpha absorbers (DLAs) are a well-studied class of absorption line systems, and yet the properties of their host galaxies remain largely unknown. To investigate the origin of these systems, we have conducted an imaging survey of 32 quasar fields with intervening DLAs between z~1.9-3.8, leveraging a technique that allows us to image galaxies at any small angular separation from the background quasars. In this paper, we present the properties of the targeted DLA sample, new imaging observations of the quasar fields, and the analysis of new and archival spectra of the background quasars.

preprint2014arXiv

Directly imaging damped Ly-alpha galaxies at z>2. III: The star formation rates of neutral gas reservoirs at z~2.7

We present results from a survey designed to probe the star formation properties of 32 damped Ly-alpha systems (DLAs) at z~2.7. By using the "double-DLA" technique that eliminates the glare of the bright background quasars, we directly measure the rest-frame FUV flux from DLAs and their neighbouring galaxies. At the position of the absorbing gas, we place stringent constraints on the unobscured star formation rates (SFRs) of DLAs to 2-sigma limits of <0.09-0.27 M/yr, corresponding to SFR surface densities <10^(-2.6)-10^(-1.5) M/yr/kpc^2. The implications of these limits for the star formation law, metal enrichment, and cooling rates of DLAs are examined. By studying the distribution of impact parameters as a function of SFRs for all the galaxies detected around these DLAs, we place new direct constraints on the bright end of the UV luminosity function of DLA hosts. We find that <13% of the hosts have SFR>2 M/yr at impact parameters b < (SFR/(M/yr))^(0.8)+6 kpc, differently from current samples of confirmed DLA galaxies. Our observations also disfavor a scenario in which the majority of DLAs arise from bright LBGs at distances 20 < b < 100 kpc. These new findings corroborate a picture in which DLAs do not originate from highly star forming systems that are coincident with the absorbers, and instead suggest that DLAs are associated with faint, possibly isolated, star-forming galaxies. Potential shortcomings of this scenario and future strategies for further investigation are discussed.

preprint2014arXiv

MUSE sneaks a peek at extreme ram-pressure stripping events. I. A kinematic study of the archetypal galaxy ESO137-001

We present MUSE observations of ESO137-001, a spiral galaxy infalling towards the center of the massive Norma cluster at z~0.0162. During the high-velocity encounter of ESO137-001 with the intracluster medium, a dramatic ram-pressure stripping event gives rise to an extended gaseous tail, traced by our MUSE observations to >30 kpc from the galaxy center. By studying the H-alpha surface brightness and kinematics in tandem with the stellar velocity field, we conclude that ram pressure has completely removed the interstellar medium from the outer disk, while the primary tail is still fed by gas from the inner regions. Gravitational interactions do not appear to be a primary mechanism for gas removal. The stripped gas retains the imprint of the disk rotational velocity to ~20 kpc downstream, without a significant gradient along the tail, which suggests that ESO137-001 is fast moving along a radial orbit in the plane of the sky. Conversely, beyond ~20 kpc, a greater degree of turbulence is seen, with velocity dispersion up to >100 km/s. For a model-dependent infall velocity of ~3000 km/s, we conclude that the transition from laminar to turbulent flow in the tail occurs on timescales >6.5 Myr. Our work demonstrates the terrific potential of MUSE for detailed studies of how ram-pressure stripping operates on small scales, providing a deep understanding of how galaxies interact with the dense plasma of the cluster environment.

preprint2014arXiv

SLUG - Stochastically Lighting Up Galaxies. II: Quantifying the Effects of Stochasticity on Star Formation Rate Indicators

The integrated light of a stellar population, measured through photometric filters that are sensitive to the presence of young stars, is often used to infer the star formation rate (SFR) for that population. However, these techniques rely on an assumption that star formation is a continuous process, whereas in reality stars form in discrete spatially- and temporally-correlated structures. This discreteness causes the light output to undergo significant time-dependent fluctuations, which, if not accounted for, introduce systematic errors in the inferred SFRs due to the intrinsic distribution of luminosities at any fix SFR. We use SLUG, a code that Stochastically Lights Up Galaxies, to simulate galaxies undergoing stochastic star formation. We then use these simulations to present a quantitative analysis of these effects and provide tools for calculating probability distribution functions of SFRs given a set of observations. We show that, depending on the SFR tracer used, stochastic fluctuations can produce non-trivial errors at SFRs as high as 1 M/yr, and biases >0.5 dex at the lowest SFRs. We emphasize that due to the stochastic behavior of blue SFR tracers, one cannot assign a deterministic single value to the SFR of an individual galaxy, but we suggest methods by which future analyses that rely on integrated-light indicators can properly account for these stochastic effects.

preprint2014arXiv

The Giant Gemini GMOS survey of z>4.4 quasars - I. Measuring the mean free path across cosmic time

We have obtained spectra of 163 quasars at $z_\mathrm{em}>4.4$ with the Gemini Multi Object Spectrometers on the Gemini North and South telescopes, the largest publicly available sample of high-quality, low-resolution spectra at these redshifts. From this homogeneous data set, we generated stacked quasar spectra in three redshift intervals at $z\sim 5$. We have modelled the flux below the rest-frame Lyman limit ($λ_\mathrm{r}<912$Å) to assess the mean free path $λ_\mathrm{mfp}^{912}$ of the intergalactic medium to HI-ionizing radiation. At mean redshifts $z_\mathrm{q}=4.56$, 4.86 and 5.16, we measure $λ_\mathrm{mfp}^{912}=(22.2\pm 2.3, 15.1\pm 1.8, 10.3\pm 1.6)h_{70}^{-1}$ proper Mpc with uncertainties dominated by sample variance. Combining our results with $λ_\mathrm{mfp}^{912}$ measurements from lower redshifts, the data are well modelled by a simple power-law $λ_\mathrm{mfp}^{912}=A[(1+z)/5]^η$ with $A=(37\pm 2)h_{70}^{-1}$ Mpc and $η= -5.4\pm 0.4$ between $z=2.3$ and $z=5.5$. This rapid evolution requires a physical mechanism -- beyond cosmological expansion -- which reduces the cosmic effective Lyman limit opacity. We speculate that the majority of HI Lyman limit opacity manifests in gas outside galactic dark matter haloes, tracing large-scale structures (e.g. filaments) whose average density (and consequently neutral fraction) decreases with cosmic time. Our measurements of the strongly redshift-dependent mean free path shortly after the completion of HI reionization serve as a valuable boundary condition for numerical models thereof. Having measured $λ_\mathrm{mfp}^{912}\approx 10$ Mpc at $z=5.2$, we confirm that the intergalactic medium is highly ionized by that epoch and that the redshift evolution of the mean free path does not show a break that would indicate a recent end to HI reionization.

preprint2014arXiv

The nature of massive black hole binary candidates: II. Spectral energy distribution atlas

Recoiling supermassive black holes (SMBHs) are considered one plausible physical mechanism to explain high velocity shifts between narrow and broad emission lines sometimes observed in quasar spectra. If the sphere of influence of the recoiling SMBH is such that only the accretion disc is bound, the dusty torus would be left behind, hence the SED should then present distinctive features (i.e. a mid-infrared deficit). Here we present results from fitting the Spectral Energy Distributions (SEDs) of 32 Type-1 AGN with high velocity shifts between broad and narrow lines. The aim is to find peculiar properties in the multi-wavelength SEDs of such objects by comparing their physical parameters (torus and disc luminosity, intrinsic reddening, and size of the 12$μ$m emitter) with those estimated from a control sample of $\sim1000$ \emph{typical} quasars selected from the Sloan Digital Sky Survey in the same redshift range. We find that all sources, with the possible exception of J1154+0134, analysed here present a significant amount of 12~$μ$m emission. This is in contrast with a scenario of a SMBH displaced from the center of the galaxy, as expected for an undergoing recoil event.

preprint2014arXiv

The Rapid Decline in Metallicity of Damped Ly-alpha Systems at z~5

We present evidence that the cosmological mean metallicity of neutral atomic hydrogen gas shows a sudden decrease at z>4.7 down to <Z>=-2.03(+0.09-0.11), which is 6 sigma deviant from that predicted by a linear fit to the data at lower redshifts. This measurement is made possible by the chemical abundance measurements of 8 new damped Ly-alpha (DLA) systems at z>4.7 observed with the Echellette Spectrograph and Imager on the Keck II telescope, doubling the number of measurements at z>4.7 to 16. Possible explanations for this sudden decrease in metallicity include a change in the physical processes that enrich the neutral gas within disks, or an increase of the covering factor of neutral gas outside disks due to a lower ultra-violet radiation field and higher density at high redshift. The later possibility would result in a new population of presumably lower metallicity DLAs, with an increased contribution to the DLA population at higher redshifts resulting in a reduced mean metallicity. Furthermore, we provide evidence of a possible decrease at z>4.7 in the comoving metal mass density of DLAs, which is flat out to z~4.3. Such a decrease is expected, as otherwise most of the metals from star-forming galaxies would reside in DLAs by z~6. While the metallicity is decreasing at high redshift, the contribution of DLAs to the total metal budget of the universe increases with redshift, with DLAs at z~4.3 accounting for ~20% as many metals as produced by Lyman break galaxies.

preprint2014arXiv

Towards a Unified Description of the Intergalactic Medium at Redshift z~2.5

We examine recent measurements of the z~2.5 intergalactic medium (IGM) which constrain the HI frequency distribution f(N) and the mean free path MFP to ionizing radiation. We argue that line-blending and the clustering of strong absorption-line systems have led previous authors to systematically overestimate the effective Lyman limit opacity, yielding too small of a MFP for the IGM. We further show that recently published measurements of f(N) at log NHI ~ 16 lie in strong disagreement, implying underestimated uncertainty from sample variance and/or systematics like line-saturation. Allowing for a larger uncertainty in the f(N) measurements, we provide a new cubic Hermite spline model for f(N) which reasonably satisfies all of the observational constraints under the assumption of randomly distributed absorption systems. We caution, however, that this formalism is invalid in light of absorber clustering and use a toy model to estimate the effects. Future work must properly account for the non-Poissonian nature of the IGM.

preprint2013arXiv

Confronting Simulations of Optically Thick Gas in Massive Halos with Observations at z=2-3

Cosmological hydrodynamic simulations predict the physical state of baryons in the circumgalactic medium (CGM), which can be directly tested via quasar absorption line observations. We use high resolution "zoom-in" simulations of 21 galaxies to characterize the distribution of neutral hydrogen around halos in the mass range M_vir~ 2x10^11 - 4x10^12 M_sun at z~2. We find that both the mass fraction of cool (T <= 3x10^4 K) gas and the covering fraction of optically-thick Lyman limit systems (LLSs) depend only weakly on halo mass, even around the critical value for the formation of stable virial shocks. The covering fraction of LLSs interior to the virial radius varies between f_c~0.05 - 0.2, with significant scatter among halos. Our simulations of massive halos (M_vir >= 10^12 M_sun) underpredict the covering fraction of optically-thick gas observed in the quasar CGM by a large factor. The reason for this discrepancy is unclear, but several possibilities are discussed. In the lower mass halos (M_vir >= 5x10^11 M_sun) hosting star-forming galaxies, the predicted covering factor agrees with observations, however current samples of quasar-galaxy pairs are too small for a conclusive comparison. To overcome this limitation, we propose a new observable: the small-scale auto-correlation function of optically-thick absorbers detected in the foreground of close quasar pairs. We show that this new observable can constrain the underlying dark halos hosting LLSs at z~2-3, as well as the characteristic size and covering factor of the CGM.

preprint2013arXiv

Dissecting the properties of optically-thick hydrogen at the peak of cosmic star formation history

We present results of a blind survey of Lyman limit systems (LLSs) detected in absorption against 105 quasars at z~3 using the blue sensitive MagE spectrograph at the Magellan Clay telescope. By searching for Lyman limit absorption in the wavelength range ~3000-4000A, we measure the number of LLSs per unit redshift l(z)=1.21+/-0.28 at z~2.8. Using a stacking analysis, we further estimate the mean free path of ionizing photons in the z~3 Universe lambda^912 = 100+/-29 Mpc/h_70.4. Combined with our LLS survey, we conclude that systems with log N_HI >= 17.5 cm^-2 contribute only ~40% to the observed mean free path at these redshifts. Further, with the aid of photo-ionization modeling, we infer that a population of ionized and metal poor systems is likely required to reproduce the metal line strengths observed in a composite spectrum of 20 LLSs with N_HI = 17.5-19 cm^-2 at z~2.6-3.0. Finally, with a simple toy model, we deduce that gas in the halos of galaxies can alone account for the totality of LLSs at z<~3, but a progressively higher contribution from the intergalactic medium is required beyond z~3.5. We also show how the weakly evolving number of LLSs per unit redshift at z<~3 can be modeled either by requiring that the spatial extent of the circumgalactic medium is redshift invariant in the last ~10 Gyr of cosmic evolution or by postulating that LLSs arise in halos that are rare fluctuations in the density field at each redshift.

preprint2013arXiv

Halpha3: an Halpha imaging survey of HI selected galaxies from ALFALFA. II. The star formation properties of galaxies in the Virgo cluster and surroundings

We present the analysis of Halpha3, an Halpha imaging survey of 409 galaxies selected from the HI Arecibo ALFALFA Survey in the Local Supercluster, including the Virgo cluster. We explore the relations between the stellar mass, the HI mass and the current, massive SFR of nearby galaxies in the Virgo cluster and we compare them with those of isolated galaxies in the Local Supercluster, disentangling the role of the environment in shaping the star formation properties of galaxies at the present cosmological epoch. We investigate the relationships between atomic neutral gas and newly formed stars in different environments, across many morphological types, and over a wide range of stellar masses adopting an updated calibration of the HI deficiency parameter. Studying the mean properties of late-type galaxies in the Local Supercluster, we find that galaxies in increasing local galaxy density conditions (or decreasing projected angular separation from M87) show a significant decrease in the HI content and in the mean specific star formation rate, along with a progressive reddening of their stellar populations. The gradual quenching of the star formation occurs outside-in, consistently with the predictions of the ram pressure model. Once considered as a whole, the Virgo cluster is effective in removing neutral hydrogen from galaxies, and this perturbation is strong enough to appreciably reduce the SFR of its entire galaxy population. An estimate of the present infall rate of 300-400 galaxies per Gyr in the Virgo cluster is obtained from the number of existing HI-rich late-type systems, assuming 200-300 Myr as the time scale for HI ablation. If the infall process has been acting at constant rate this would imply that the Virgo cluster has formed approximately 2 Gyr ago, consistently with the idea that Virgo is in a young state of dynamical evolution.

preprint2013arXiv

Halpha3: an Halpha imaging survey of HI selected galaxies from ALFALFA. III. Nurture shapes up the Hubble sequence in the Great Wall

We present the analysis of Halpha3, an Halpha imaging survey of galaxies selected from the HI ALFALFA Survey in the Coma Supercluster. By using the Halpha line as a tracer of the "instantaneous" star formation, complemented with optical colors from SDSS we explore the hypothesis that a morphological sequence of galaxies of progressively earlier type, lower gas-content exists in the neighborhood of the Coma cluster, with specific star formation activity decreasing with increasing local galaxy density and velocity dispersion. In the dwarf regime (8.5<\log(M*)< 9.5) we identify a 4-step sequence of galaxies with progressively redder colors, i.e. of decreasing specific star formation, from (1) HI-rich Late-Type Galaxies belonging to the blue-cloud exhibiting extended plus nuclear star formation, (2) HI-poor LTGs with nuclear star formation only, (3) HI-poor galaxies with no star formation either extended or nuclear, but with nuclear Post-Star-Burst signature,(4) Early-type Galaxies in the red-sequence, with no gas or star formation on all scales. Along this sequence the quenching of the star formation proceeds radially outside-in. The progression toward redder colors found along this "morphological" (gas content) sequence is comparable to the one obtained increasing the local galaxy density, from the cosmic filaments (1,2), to the rich clusters (2,3,4). In the dwarf regime we find evidence for evolution of HI-rich LTGs into ETGs, passing through HI-poor and PSB galaxies, driven by the environment. We identify ram-pressure as the mechanism most likely responsible for this transformation. We conclude that infall of low-mass galaxies has proceeded for the last 7.5 Gyr building up the Coma cluster at a rate of approximately 100 galaxies per Gyr.

preprint2013arXiv

Halpha3: an Halpha imaging survey of HI selected galaxies from ALFALFA. IV. Structure of galaxies in the Local and Coma Superclusters

We present the analysis of the galaxy structural parameters from Halpha3, an Halpha narrow-band imaging follow-up survey of ~800 galaxies selected from the HI ALFALFA Survey in the Local and Coma Superclusters. Taking advantage of Halpha3 which provides the complete census of the recent star-forming, HI-rich galaxies in the local universe, we aim to investigate the structural parameters of both the young (<10 Myr) and the old (>1 Gyr) stellar populations. By comparing the sizes of these stellar components we investigated the spatial scale on which galaxies are growing at the present cosmological epoch and the role of the environment in quenching the star-formation activity. We computed the concentration, asymmetry, and clumpiness structural parameters. To quantify the sizes we computed half-light radii and a new parameter dubbed EW/r. The concentration index computed in the r band depends on the stellar mass and on the Hubble type, these variables being related since most massive galaxies are bulge dominated thus most concentrated. Going toward later spirals and irregulars both the concentration index and the stellar mass decrease. Blue Compact dwarfs represent an exception since they have similar stellar mass but they are more concentrated than dwarf irregulars. The asymmetry and the clumpiness increase along the spiral sequence then they decrease going into the dwarf regime, where the light distribution is smooth and more symmetric. When measured on Halpha images, the CAS parameters do not exhibit obvious correlations with Hubble type. We found that the concentration index is the main parameter that describes the current growth of isolated galaxies but, for a fixed concentration, the stellar mass plays a second order role. At the present epoch, massive galaxies are growing inside-out, conversely the dwarfs are growing on the scale of their already assembled mass.

preprint2013arXiv

Metal abundances in the high-redshift intergalactic medium

Twenty years of high-resolution spectroscopy at the 8-10 m class telescopes have drastically expanded our view of the gas-phase metallicity in the z>2 universe. This contribution briefly summarizes how these studies reveal a widespread metal pollution in the intergalactic medium with a median abundance [C/H] ~ -3.5 at z~3 that is increasing by a factor of ~2-3 from z~4.3 to z~2.4. At the higher densities that are typical of galactic halos, observations uncover a metallicity spread of five orders of magnitude in Lyman limit systems, ranging from super-solar ([M/H]~+0.7) to pristine ([M/H]< -4) gas clouds. Finally, the neutral damped Ly-alpha systems are enriched to a median metallicity of [M/H] ~ -1.5 that slowly declines with redshift up to z ~ 4.5, at which point it appears to more rapidly evolve as one approaches the end of reionization.

preprint2013arXiv

On the Redshift of the Very High Energy Blazar 3C 66A

As a bright gamma-ray source, 3C 66A is of great interest to the high-energy astrophysics community, having a potential for placing cosmological constraints on models for the extragalactic background light (EBL) and the processes which contribute to this photon field. No firm spectroscopic redshift measurement has been possible for this blazar due to a lack of intrinsic emission and absorption features in optical spectra. We present new far-ultraviolet spectra from the Hubble Space Telescope/Cosmic Origins Spectrograph (HST/COS) of the BL Lac object 3C 66A covering the wavelength range 1132 to 1800 Angstroms . The data show a smooth continuum with intergalactic medium absorption features which can be used to place a firm lower limit on the blazar redshift of z>0.3347. An upper limit is set by statistically treating the non-detection of additional absorbers beyond z=0.3347, indicating a redshift of less than 0.41 at 99% confidence and ruling out z>0.444 at 99.9% confidence. We conclude by showing how the redshift limits derived from the COS spectra remove the potential for this gamma-ray emitting blazar to place an upper limit on the flux of the EBL using high energy data from a flare in October of 2009.

preprint2013arXiv

The First AllWISE Proper Motion Discovery: WISEA J070720.50+170532.7

While quality checking a new motion-aware co-addition of all 12.5 months of WISE data, the source WISE J070720.48+170533.0 was found to have moved 0.9 arc-sec in 6 months. Backtracking this motion allowed us to identify this source as 2MASS J07071961+1705464 and with several entries in the USNO B catalog. An astrometric fit to these archival data gives a proper motion of 1793 +/- 2 mas/yr and a parallax of 35 +/- 42 mas. Photometry from WISE, 2MASS and the POSS can be fit reasonably well by a blackbody with T = 3658 K and an angular radius of 4.36E-11 radians. No clear evidence of H_2 collision-induced absorption is seen in the near-IR. An optical spectrum shows broad deep CaH bands at 638 and 690 nm, broad deep Na D at 598.2 nm, and weak or absent TiO, indicating that this source is an ultra-subdwarf M star with a radial velocity about -21 +/- 18 km/sec relative to the Sun. Given its apparent magnitude, the distance is about 39 +/- 9 pc and the tangential velocity is probably about 330 km/sec, but a more precise parallax is needed to be certain.

preprint2013arXiv

The nature of massive black hole binary candidates: I. Spectral properties and evolution

Theoretically, bound binaries of massive black holes are expected as the natural outcome of mergers of massive galaxies. From the observational side, however, massive black hole binaries remain elusive. Velocity shifts between narrow and broad emission lines in quasar spectra are considered a promising observational tool to search for spatially unresolved, dynamically bound binaries. In this series of papers we investigate the nature of such candidates through analyses of their spectra, images and multi-wavelength spectral energy distributions. Here we investigate the properties of the optical spectra, including the evolution of the broad line profiles, of all the sources identified in our previous study. We find a diverse phenomenology of broad and narrow line luminosities, widths, shapes, ionization conditions and time variability, which we can broadly ascribe to 4 classes based on the shape of the broad line profiles: 1) Objects with bell-shaped broad lines with big velocity shifts (>1000 km/s) compared to their narrow lines show a variety of broad line widths and luminosities, modest flux variations over a few years, and no significant change in the broad line peak wavelength. 2) Objects with double-peaked broad emission lines tend to show very luminous and broadened lines, and little time variability. 3) Objects with asymmetric broad emission lines show a broad range of broad line luminosities and significant variability of the line profiles. 4) The remaining sources tend to show moderate to low broad line luminosities, and can be ascribed to diverse phenomena. We discuss the implications of our findings in the context of massive black hole binary searches.

preprint2012arXiv

65 kpc of ionized gas trailing behind NGC 4848 during its first crossing of the Coma cluster

In a 5 hour Halpha exposure of the N-W region of the Coma cluster with the 2.1m telescope at SPM (Mx) we discovered a 65 kpc cometary emission of ionized gas trailing behind the SBab galaxy NGC 4848. The tail points in the opposite direction of the cluster center, in the same direction where stripped HI has been detected in previous observations. The galaxy shows bright HII regions in an inner ring-like pattern, where the star formation takes place at the prodigious rate of 8.9 Msun/yr. From the morphology of the galaxy and of the trailing material, we infer that the galaxy is suffering from ram pressure due to its high velocity motion through the cluster IGM. We estimate that 4 x 10^9 Msun of gas is swept out from the galaxy forming the tail. Given the ambient conditions in the Coma cluster (rho = 6.3 x 10^-27 g/cm^3; sigma_vel = 940 km/s) simulations predict that the ram pressure mechanism is able to remove such an amount of gas in less than 200 Myr. This, combined with the geometry of the interaction, indicating radial infall into the cluster, leads to the conclusion that NGC 4848 is caught in its first passage through the dense cluster environment.

preprint2012arXiv

A search of CO emission lines in blazars: the low molecular gas content of BL Lac objects compared to quasars

BL Lacertae (Lac) objects that are detected at very-high energies (VHE) are of fundamental importance to study multiple astrophysical processes, including the physics of jets, the properties of the extragalactic background light and the strength of the intergalactic magnetic field. Unfortunately, since most blazars have featureless optical spectra that preclude a redshift determination, a substantial fraction of these VHE extragalactic sources cannot be used for cosmological studies. To assess whether molecular lines are a viable way to establish distances, we have undertaken a pilot program at the IRAM 30m telescope to search for CO lines in three BL Lac objects with known redshifts. We report a positive detection of M_H2 ~ 3x10^8 Msun toward 1ES 1959+650, but due to the poor quality of the baseline, this value is affected by a large systematic uncertainty. For the remaining two sources, W Comae and RGB J0710+591, we derive 3sigma upper limits at, respectively, M_H2 < 8.0x10^8 Msun and M_H2 < 1.6x10^9 Msun, assuming a line width of 150 km/s and a standard conversion factor alpha=4 M_sun/(K km/s pc^2). If these low molecular gas masses are typical for blazars, blind redshift searches in molecular lines are currently unfeasible. However, deep observations are still a promising way to obtain precise redshifts for sources whose approximate distances are known via indirect methods. Our observations further reveal a deficiency of molecular gas in BL Lac objects compared to quasars, suggesting that the host galaxies of these two types of active galactic nuclei (AGN) are not drawn from the same parent population. Future observations are needed to assess whether this discrepancy is statistically significant, but our pilot program shows how studies of the interstellar medium in AGN can provide key information to explore the connection between the active nuclei and the host galaxies.

preprint2012arXiv

Halpha3: Halpha imaging survey of HI selected galaxies from ALFALFA

We present Halpha3 (acronym for Halpha-alpha-alpha), an Halpha narrow-band imaging survey of ~400 galaxies selected from the HI Arecibo Legacy Fast ALFA Survey (ALFALFA) in the Local Supercluster, including the Virgo cluster. By using hydrogen recombination lines as a tracer of recent star formation, we aim to investigate the relationships between atomic neutral gas and newly formed stars in different environments (cluster and field), morphological types (spirals and dwarfs), and over a wide range of stellar masses (~10^7.5-10^11.5 Msun). We image in Halpha+[NII] all the galaxies that contain more than 10^7 Msun of neutral atomic hydrogen in the sky region 11^h < R.A. <16^h 4^o < Dec. <16^o; 350< cz <2000 km/s using the San Pedro Martir 2m telescope. This survey provides a complete census of the star formation in HI rich galaxies of the local universe. We present the properties of the galaxy sample, together with Halpha fluxes and equivalent widths. We find an excellent agreement between the fluxes determined from our images in apertures of 3 arcsec diameter and the fluxes derived from the SDSS spectral database. From the Halpha fluxes corrected for galactic and internal extinction and for [NII] contamination we derive the global star formation rates (SFRs).

preprint2012arXiv

On the redshift of the blazar PKS0447-439

PKS0447-439 is one of the BL Lacertae objects that have been detected at very high energy. There has been a recent report of a lower limit of z>1.246 for the redshift of this blazar, challenging the current paradigm in which very high-energy gamma-rays cannot freely propagate in the z>1 universe. In this research note, we present a new MagE/Magellan spectrum of PKS0447-439 with exquisite signal-to-noise (S/N>150 at 6500A). Our analysis confirms the presence of the previously-reported absorption line at 6280A, which we identify, however, with a known telluric absorption, invalidating the claim that this blazar lies at z>1. Since no other extragalactic spectral features are detected, we cannot establish a redshift based on our spectrum.

preprint2012arXiv

Stripped gas as fuel for newly formed HII regions in the encounter between VCC1249 and M49: a unified picture from NGVS and GUViCS

Context: We study the peculiar interacting galaxy system of VCC1249/M49 located in the core of the Virgo B subcluster. Owing to a recent interaction between the dwarf galaxy VCC1249 and the halo gas of the gE M49, neutral hydrogen has been displaced from the interstellar medium of this dwarf into the Virgo ICM. Observations also reveal multiple compact star-forming regions that are embedded in this HI cloud, with a projected separation up to 13 kpc from VCC1249 in the northwest direction. Aims: Motivated by recent NUV imaging from GUViCS of the VCC1249/M49 system that shows significant ongoing/recent star formation in the compact regions, we aim to constrain the origin of these outlying HII regions with a multi-wavelength approach. Methods: Using deep optical (u, g, i, z) imaging from NGVS and new Halpha imaging obtained at the San Pedro Martir observatory together with Keck long-slit spectroscopy, we characterize the SFR, ages, and metallicity of VCC1249 and its outlying compact regions. Moreover, we analyze the color and luminosity profile of the galaxy to investigate its recent interaction with M49. Results: Our new observations indicate that VCC1249 underwent a recent interaction with M49 in which both ram-pressure stripping and tidal interaction occured. The joint action of the two mechanisms led to the removal of the HI gas from the ISM of VCC1249, while the gravitational tides triggered the stellar tail and counter-tail of VCC1249. Our SED analysis reveals that the star formation in this galaxy was truncated around 200 Myr ago and that the outlying HII regions were born in situ about 10 Myr ago out of pre-enriched gas removed from the dwarf galaxy. These observations also reveal that interactions between central and satellite galaxies similar to the one between VCC1249/M49 may be an effective way of dispersing metals into the halos of massive galaxies.

preprint2011arXiv

Absorption line systems in simulated galaxies fed by cold streams

Hydro cosmological simulations reveal that massive galaxies at high redshift are fed by long narrow streams of merging galaxies and a smoother component of cold gas. We post-process seven high-resolution simulated galaxies with radiative transfer to study the absorption characteristics of the gas in galaxies and streams, in comparison with the statistics of observed absorption-line systems. We find that much of the stream gas is ionized by UV radiation from background and local stellar sources, but still optically thick (N_HI > 10^17 cm^-2) so that the streams appear as Lyman-limit systems (LLSs). At z>3, the fraction of neutral gas in streams becomes non-negligible, giving rise to damped Lyman-alpha (DLA) absorbers as well. The gas in the central and incoming galaxies remains mostly neutral, responsible for DLAs. Within one (two) virial radii, the covering factor of optically thick gas is <25% (10%) for LLSs and <5% (1%) for DLAs, slowly declining with time following the universal expansion. Nevertheless, galaxies and their cold streams in the studied mass range, M_vir = 10^10 - 10^12 Msun, account for >30% of the observed absorbers in the foreground of quasars, the rest possibly arising from smaller galaxies or the intergalactic medium. The mean metallicity in the streams is ~1% solar, much lower than in the galaxies. The simulated galaxies reproduce the Lyalpha-absorption equivalent widths observed around Lyman-break galaxies, but they severely underpredict the equivalent widths in metal lines, suggesting that the latter may arise from outflows. We conclude that the observed metal-poor LLSs are likely detections of the predicted cold streams. Revised analysis of the observed LLSs kinematics and simulations with more massive outflows in conjunction with the inflows may enable a clearer distinction between the signatures of the various gas modes.

preprint2011arXiv

Detection of Pristine Gas Two Billion Years after the Big Bang

In the current cosmological model, only the three lightest elements were created in the first few minutes after the Big Bang; all other elements were produced later in stars. To date, however, heavy elements have been observed in all astrophysical environments. We report the detection of two gas clouds with no discernible elements heavier than hydrogen. These systems exhibit the lowest heavy-element abundance in the early universe and thus are potential fuel for the most metal poor halo stars. The detection of deuterium in one system at the level predicted by primordial nucleosynthesis provides a direct confirmation of the standard cosmological model. The composition of these clouds further implies that the transport of heavy elements from galaxies to their surroundings is highly inhomogeneous.

preprint2011arXiv

SLUG - Stochastically Lighting Up Galaxies I: Methods and Validating Tests

The effects of stochasticity on the luminosities of stellar populations are an often neglected but crucial element for understanding populations in the low mass or low star formation rate regime. To address this issue, we present SLUG, a new code to "Stochastically Light Up Galaxies". SLUG synthesizes stellar populations using a Monte Carlo technique that treats stochastic sampling properly including the effects of clustering, the stellar initial mass function, star formation history, stellar evolution, and cluster disruption. This code produces many useful outputs, including i) catalogs of star clusters and their properties, such as their stellar initial mass distributions and their photometric properties in a variety of filters, ii) two dimensional histograms of color-magnitude diagrams of every star in the simulation, iii) and the photometric properties of field stars and the integrated photometry of the entire simulated galaxy. After presenting the SLUG algorithm in detail, we validate the code through comparisons with starburst99 in the well-sampled regime, and with observed photometry of Milky Way clusters. Finally, we demonstrate the SLUG's capabilities by presenting outputs in the stochastic regime. SLUG is publicly distributed through the website http://sites.google.com/site/runslug/.

preprint2011arXiv

Stochastic star formation and a (nearly) uniform stellar initial mass function

Recent observations indicate a lower Halpha to FUV ratio in dwarf galaxies than in brighter systems, a trend that could be explained by a truncated and/or steeper IMF in small galaxies. However, at low star formation rates (SFRs), the Halpha to FUV ratio can vary due to stochastic sampling even for a universal IMF, a hypothesis that has, prior to this work, received limited investigation. Using SLUG, a fully stochastic code for synthetic photometry in star clusters and galaxies, we compare the Halpha and FUV luminosity in a sample of ~450 nearby galaxies with models drawn from a universal Kroupa IMF and a modified IMF, the integrated galactic initial mass function (IGIMF). Once random sampling and time evolution are included, a Kroupa IMF convolved with the cluster mass function reproduces the observed Halpha distribution at all FUV luminosities, while a truncated IMF as implemented in current IGIMF models underpredicts the Halpha luminosity by more than an order of magnitude at the lowest SFRs. We conclude that the observed luminosity is the result of the joint probability distribution function of the SFR, cluster mass function, and a universal IMF, consistent with parts of the IGIMF theory, but that a truncation in the IMF in clusters is inconsistent with the observations. Future work will examine stochastic star formation and its time dependence in detail to study whether random sampling can explain other observations that suggest a varying IMF.

preprint2010arXiv

Directly imaging damped Lyman-alpha galaxies at z>2. I: Methodology and First Results

We present the methodology for, and the first results from, a new imaging program aimed at identifying and characterizing the host galaxies of damped Lyman-alpha absorbers (DLAs) at z>2. We target quasar sightlines with multiple optically-thick HI absorbers and use the higher-redshift system as a "blocking filter" (via its Lyman-limit absorption) to eliminate all far-ultraviolet (FUV) emission from the quasar. This allows us to directly image the rest-frame FUV continuum emission of the lower-redshift DLA, without any quasar contamination and with no bias towards large impact parameters. We introduce a formalism based on galaxy number counts and Bayesian statistics with which we quantify the probability that a candidate is the DLA host galaxy. This method will allow the identification of a bona fide sample of DLAs that are too faint to be spectroscopically confirmed. The same formalism can be adopted to the study of other quasar absorption line systems (e.g. MgII absorbers). We have applied this imaging technique to two QSO sightlines. For the z~2.69 DLA towards J073149+285449, a galaxy with impact parameter b=1.54"=11.89 kpc and implied star formation rate (SFR) of ~5 M/yr is identified as the most reliable candidate. In the case of the z~2.92 DLA towards J211444-005533, no likely host is found down to a 3-sigma SFR limit of 1.4 M/yr. Studying the HI column density as a function of the impact parameter, including 6 DLAs with known hosts from the literature, we find evidence that the observed HI distribution is more extended than what is generally predicted from numerical simulation.

preprint2010arXiv

SLUG: A new way to Stochastically Light Up Galaxies

We present SLUG, a new code to "Stochastically Light Up Galaxies". SLUG populates star clusters by randomly drawing stars from an initial mass function (IMF) and then following their time evolution with stellar models and an observationally-motivated prescription for cluster disruption. For a choice of star formation history, metallicity, and IMF, SLUG outputs synthetic photometry for clusters and field stars with a proper treatment of stochastic star formation. SLUG generates realistic distributions of star clusters, demonstrating the range of properties that result from finite sampling of an IMF and a random distribution of ages. The simulated data sets provide a quantitative means to address open problems in studies of star formation in galaxies and clusters, such as a test for IMF variations that are suggested by the systematic deficiency in the H-alpha/UV ratio in outer disks or in dwarf galaxies. SLUG will be made publicly available through the website http://sites.google.com/site/runslug/.

preprint2010arXiv

Testing models for molecular gas formation in galaxies: hydrostatic pressure or gas and dust shielding?

Stars in galaxies form in giant molecular clouds that coalesce when the atomic hydrogen is converted into molecules. There are currently two dominant models for what property of the galactic disk determines its molecular fraction: either hydrostatic pressure driven by the gravity of gas and stars, or a combination of gas column density and metallicity. To assess the validity of these models, we compare theoretical predictions to the observed atomic gas content of low-metallicity dwarf galaxies with high stellar densities. The extreme conditions found in these systems are optimal to distinguish the two models, otherwise degenerate in nearby spirals. Locally, on scales <100 pc, we find that the state of the interstellar medium is mostly sensitive to the gas column density and metallicity rather than hydrostatic pressure. On larger scales where the average stellar density is considerably lower, both pressure and shielding models reproduce the observations, even at low metallicity. We conclude that models based on gas and dust shielding more closely describe the process of molecular formation, especially at the high resolution that can be achieved in modern galaxy simulations or with future radio/millimeter arrays.

preprint2009arXiv

Molecular hydrogen deficiency in HI-poor galaxies and its implications for star formation

We use a sample of 47 homogeneous and high sensitivity CO images taken from the Nobeyama and BIMA surveys to demonstrate that, contrary to common belief, a significant number (~40%) of HI-deficient nearby spiral galaxies are also depleted in molecular hydrogen. While HI-deficiency by itself is not a sufficient condition for molecular gas depletion, we find that H2 reduction is associated with the removal of HI inside the galaxy optical disk. Those HI-deficient galaxies with normal H2 content have lost HI mainly from outside their optical disks, where the H2 content is low in all galaxies. This finding is consistent with theoretical models in which the molecular fraction in a galaxy is determined primarily by its gas column density. Our result is supported by indirect evidence that molecular deficient galaxies form stars at a lower rate or have dimmer far infrared fluxes than gas rich galaxies, as expected if the star formation rate is determined by the molecular hydrogen content. Our result is consistent with a scenario in which, when the atomic gas column density is lowered inside the optical disk below the critical value required to form molecular hydrogen and stars, spirals become quiescent and passive evolving systems. We speculate that this process would act on the time-scale set by the gas depletion rate and might be a first step for the transition between the blue and red sequence observed in the color-magnitude diagram.