Refine
Document Type
- Article (16)
- Doctoral Thesis (1)
Language
- English (17)
Is part of the Bibliography
- yes (17) (remove)
Keywords
- galaxies: evolution (17) (remove)
Institute
The work done during the PhD studies has been focused on measurements of distribution functions of rotating galaxies using integral field spectroscopy observations.
Throughout the main body of research presented here we have been using CALIFA (Calar Alto Legacy Integral Field Area) survey stellar velocity fields to obtain robust measurements of circular velocities for rotating galaxies of all morphological types. A crucial part of the work was enabled by well-defined CALIFA sample selection criteria: it enabled reconstructing sample-independent distributions of galaxy properties.
In Chapter 2, we measure the distribution in absolute magnitude - circular velocity space for a well-defined sample of 199 rotating CALIFA galaxies using their stellar kinematics. Our aim in this analysis is to avoid subjective selection criteria and to take volume and large-scale structure factors into account. Using stellar velocity fields instead of gas emission line kinematics allows including rapidly rotating early type galaxies. Our initial sample contains 277 galaxies with available stellar velocity fields and growth curve r-band photometry. After rejecting 51 velocity fields that could not be modelled due to the low number of bins, foreground contamination or significant interaction we perform Markov Chain Monte Carlo (MCMC) modelling of the velocity fields, obtaining the rotation curve and kinematic parameters and their realistic uncertainties. We perform an extinction correction and calculate the circular velocity v_circ accounting for pressure support a given galaxy has. The resulting galaxy distribution on the M_r - v_circ plane is then modelled as a mixture of two distinct populations, allowing robust and reproducible rejection of outliers, a significant fraction of which are slow rotators. The selection effects are understood well enough that the incompleteness of the sample can be corrected and the 199 galaxies can be weighted by volume and large-scale structure factors enabling us to fit a volume-corrected Tully-Fisher relation (TFR). More importantly, we also provide the volume-corrected distribution of galaxies in the M_r - v_circ plane, which can be compared with cosmological simulations. The joint distribution of the luminosity and circular velocity space densities, representative over the range of -20 > M_r > -22 mag, can place more stringent constraints on the galaxy formation and evolution scenarios than linear TFR fit parameters or the luminosity function alone.
In Chapter 3, we measure one of the marginal distributions of the M_r - v_circ distribution: the circular velocity function of rotating galaxies. The velocity function is a fundamental observable statistic of the galaxy population, being of a similar importance as the luminosity function, but much more difficult to measure. We present the first directly measured circular velocity function that is representative between 60 < v_circ < 320 km s^-1 for galaxies of all morphological types at a given rotation velocity. For the low mass galaxy population 60 < v_circ < 170 km s^-1, we use the HIPASS velocity function. For the massive galaxy population 170 < v_circ < 320 km s^-1, we use stellar circular velocities from CALIFA. The CALIFA velocity function includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies. It has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both velocity functions can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed velocity function is compared to velocity functions derived from cosmological simulations of the z = 0 galaxy population. We find that dark matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic Illustris simulations fare better, but still do not fully reproduce observations.
In Chapter 4, we present some other work done during the PhD studies, namely, a method that improves the precision of specific angular measurements by combining simultaneous Markov Chain Monte Carlo modelling of ionised gas 2D velocity fields and HI linewidths. To test the method we use a sample of 25 galaxies from the Sydney-AAO Multi-object Integral field (SAMI) survey that had matching ALFALFA HI linewidths. Such a method allows constraining the rotation curve both in the inner regions of a galaxy and in its outskirts, leading to increased precision of specific angular momentum measurements. It could be used to further constrain the observed relation between galaxy mass, specific angular momentum and morphology (Obreschkow & Glazebrook 2014).
Mathematical and computational methods are presented in the appendices.
We use a background quasar to detect the presence of circumgalactic gas around a z = 0.91 low-mass star-forming galaxy. Data from the new Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope show that the galaxy has a dust-corrected star formation rate (SFR) of 4.7 +/- 2.0. M-circle dot yr(-1), with no companion down to 0.22 M-circle dot yr(-1) (5 sigma) within 240 h(-1) kpc ("30"). Using a high-resolution spectrum of the background quasar, which is fortuitously aligned with the galaxy major axis (with an azimuth angle alpha of only 15 degrees), we find, in the gas kinematics traced by low-ionization lines, distinct signatures consistent with those expected for a "cold-flow disk" extending at least 12 kpc (3 x R-1/2). We estimate the mass accretion rate M-in to be at least two to three times larger than the SFR, using the geometric constraints from the IFU data and the H (I) column density of log N-H (I)/cm(-2) similar or equal to 20.4 obtained from a Hubble Space Telescope/COS near-UV spectrum. From a detailed analysis of the low-ionization lines (e.g., Zn II, Cr II, Ti II, MnII, Si II), the accreting material appears to be enriched to about 0.4 Z(circle dot) (albeit with large uncertainties: log Z/Z(circle dot) = -0.4 +/- 0.4), which is comparable to the galaxy metallicity (12 + log O/H = 8.7 +/- 0.2), implying a large recycling fraction from past outflows. Blueshifted Mg II and Fe II absorptions in the galaxy spectrum from the MUSE data reveal the presence of an outflow. The Mg II and Fe II absorption line ratios indicate emission infilling due to scattering processes, but the MUSE data do not show any signs of fluorescent Fe II* emission.
Observational studies have revealed that galaxy pairs tend to have lower gas-phase metallicity than isolated galaxies. This metallicity deficiency can be caused by inflows of low-metallicity gas due to the tidal forces and gravitational torques associated with galaxy mergers, diluting the metal content of the central region. In this work we demonstrate that such metallicity dilution occurs in state-of-the-art cosmological simulations of galaxy formation. We find that the dilution is typically 0.1 dex for major mergers, and is noticeable at projected separations smaller than 40 kpc. For minor mergers the metallicity dilution is still present, even though the amplitude is significantly smaller. Consistent with previous analysis of observed galaxies we find that mergers are outliers from the fundamental metallicity relation, with deviations being larger than expected for a Gaussian distribution of residuals. Our large sample of mergers within full cosmological simulations also makes it possible to estimate how the star formation rate enhancement and gas consumption timescale behave as a function of the merger mass ratio. We confirm that strong starbursts are likely to occur in major mergers, but they can also arise in minor mergers if more than two galaxies are participating in the interaction, a scenario that has largely been ignored in previous work based on idealised isolated merger simulations.
Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterise the gas in the vicinity of Milky Way and Andromeda analogues in the hestia (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The hestia suite comprise of a set of three high-resolution arepo-based simulations of the LG, run using the Auriga galaxy formation model.
For this paper, we focus only on the 𝑧 = 0 simulation datasets and generate mock skymaps along with a power spectrum analysis to show that the distributions of ions tracing low-temperature gas (H i and Si iii) are more clumpy in comparison to warmer gas tracers (O vi, O vii and O viii). We compare to the spectroscopic CGM observations of M31 and low-redshift galaxies.
hestia under-produces the column densities of the M31 observations, but the simulations are consistent with the observations of low-redshift galaxies. A possible explanation for these findings is that the spectroscopic observations of M31 are contaminated by gas residing in the CGM of the Milky Way.
We trace the specific star formation rate (sSFR) of massive star-forming galaxies (greater than or similar to 10(10)M(circle dot)) from z similar to 2 to 7. Our method is substantially different from previous analyses, as it does not rely on direct estimates of star formation rate, but on the differential evolution of the galaxy stellar mass function (SMF). We show the reliability of this approach by means of semianalytical and hydrodynamical cosmological simulations. We then apply it to real data, using the SMFs derived in the COSMOS and CANDELS fields. We find that the sSFR is proportional to (1 + z)(1.1) (+/-) (0.2) at z > 2, in agreement with other observations but in tension with the steeper evolution predicted by simulations from z similar to 4 to 2. We investigate the impact of several sources of observational bias, which, however, cannot account for this discrepancy. Although the SMF of high-redshift galaxies is still affected by significant errors, we show that future large-area surveys will substantially reduce them, making our method an effective tool to probe the massive end of the main sequence of star-forming galaxies.
Emission signatures from galactic winds provide an opportunity to directly map the outflowing gas, but this is traditionally challenging because of the low surface brightness. Using very deep observations (27 h) of the Hubble Deep Field South with the Multi Unit Spectroscopic Explorer (MUSE) instrument, we identify signatures of an outflow in both emission and absorption from a spatially resolved galaxy at z = 1.29 with a stellar mass M-star = 8 x 10(9) M-circle dot, star formation rate SFR = 77(-25)(+40) M-circle dot yr(-1), and star formation rate surface brightness Sigma(SFR) = 1.6 M-circle dot kpc(-2) within the [OII] lambda lambda 3727, 3729 half-light radius R-1/2, ([OII]) = 2.76 +/- 0.17 kpc. From a component of the strong resonant Mg II and Fe II absorptions at -350 km s(-1), we infer a mass outflow rate that is comparable to the star formation rate. We detect non-resonant Fe II* emission, at lambda 2365, lambda 2396, lambda 2612, and lambda 2626, at 1.2-2.4-1.5-2.7 x 10-(18) erg s(-1) cm(-2) respectively. The flux ratios are consistent with the expectations for optically thick gas. By combining the four non-resonant Fe II* emission lines, we spatially map the Fe II* emission from an individual galaxy for the first time. The Fe II* emission has an elliptical morphology that is roughly aligned with the galaxy minor kinematic axis, and its integrated half-light radius, R-1/2, (Fe II*) = 4.1 +/- 0.4 kpc, is 70% larger than the stellar continuum (R-1/2,(star) similar or equal to 2.34 +/- 0.17) or the [O II] nebular line. Moreover, the Fe II* emission shows a blue wing extending up to -400 km s(-1), which is more pronounced along the galaxy minor kinematic axis and reveals a C-shaped pattern in a p - v diagram along that axis. These features are consistent with a bi-conical outflow.
It is well known that satellite galaxies are not isotropically distributed among their host galaxies as suggested by most interpretations of the Λ cold dark matter (ΛCDM) model. One type of anisotropy recently detected in the Sloan Digital Sky Survey (and seen when examining the distribution of satellites in the Local Group and in the Centaurus group) is a tendency to be so-called lopsided. Namely, in pairs of galaxies (like Andromeda and the Milky Way) the satellites are more likely to inhabit the region in between the pair, rather than on opposing sides. Although recent studies found a similar set-up when comparing pairs of galaxies in ΛCDM simulations indicating that such a set-up is not inconsistent with ΛCDM, the origin has yet to be explained. Here we examine the origin of such lopsided set-ups by first identifying such distributions in pairs of galaxies in numerical cosmological simulations, and then tracking back the orbital trajectories of satellites (which at z = 0 display the effect). We report two main results: first, the lopsided distribution was stronger in the past and weakens towards z = 0. Secondly, the weakening of the signal is due to the interaction of satellite galaxies with the pair. Finally, we show that the z = 0 signal is driven primarily by satellites that are on first approach, who have yet to experience a ‘flyby’. This suggests that the signal seen in the observations is also dominated by dynamically young accretion events.
Galaxies are surrounded by massive gas reservoirs ( i.e. the circumgalactic medium; CGM) which play a key role in their evolution. The properties of the CGM, which are dependent on a variety of internal and environmental factors, are often inferred from absorption line surveys which rely on a limited number of single lines-of-sight. In this work we present an analysis of 28 galaxy haloes selected from the Auriga project, a cosmological magneto-hydrodynamical zoom-in simulation suite of isolated MilkyWay-mass galaxies, to understand the impact of CGM diversity on observational studies. Although the Auriga haloes are selected to populate a narrow range in halo mass, our work demonstrates that the CGM of L-star galaxies is extremely diverse: column densities of commonly observed species span similar to 3-4 dex and their covering fractions range from similar to 5 to 90 per cent. Despite this diversity, we identify the following correlations: 1) the covering fractions ( CF) of hydrogen and metals of the Auriga haloes positively correlate with stellar mass, 2) the CF of H I, C IV, and Si II anticorrelate with active galactic nucleus luminosity due to ionization effects, and 3) the CF of H I, C IV, and Si II positively correlate with galaxy disc fraction due to outflows populating the CGM with cool and dense gas. The Auriga sample demonstrates striking diversity within the CGM of L-star galaxies, which poses a challenge for observations reconstructing CGM characteristics from limited samples, and also indicates that long-term merger assembly history and recent star formation are not the dominant sculptors of the CGM.
Galaxies are surrounded by sizeable gas reservoirs which host a significant amount of metals: the circum-galactic medium (CGM). The CGM acts as a mediator between the galaxy and the extragalactic medium. However, our understanding of how galaxy mergers, a major evolutionary transformation, impact the CGM remains deficient. We present a theoretical study of the effect of galaxy mergers on the CGM. We use hydrodynamical cosmological zoom-in simulations of a major merger selected from the Illustris project such that the z = 0 descendant has a halo mass and stellar mass comparable to the Milky Way. To study the CGM we then re-simulated this system at a 40 times better mass resolution, and included detailed post-processing ionization modelling. Our work demonstrates the effect the merger has on the characteristic size of the CGM, its metallicity, and the predicted covering fraction of various commonly observed gas-phase species, such as H I, C IV, and O VI. We show that merger-induced outflows can increase the CGM metallicity by 0.2-0.3 dex within 0.5 Gyr post-merger. These effects last up to 6 Gyr post-merger. While the merger increases the total metal covering fractions by factors of 2-3, the covering fractions of commonly observed UV ions decrease due to the hard ionizing radiation from the active galactic nucleus, which we model explicitly. Our study of the single simulated major merger presented in this work demonstrates the significant impact that a galaxy interaction can have on the size, metallicity, and observed column densities of the CGM.
We present a study of the [O III]/[O II] ratios of star-forming galaxies drawn from Multi-Unit Spectroscopic Explorer (MUSE) data spanning a redshift range 0.28 < z < 0.85. Recently discovered Lyman continuum (LyC) emitters have extremely high oxygen line ratios: [O III]lambda 5007/[O II]lambda lambda 3726, 3729 > 4. Here we aim to understand the properties and the occurrences of galaxies with such high line ratios. Combining data from several MUSE Guaranteed Time Observing (GTO) programmes, we select a population of star-forming galaxies with bright emission lines, from which we draw 406 galaxies for our analysis based on their position in the z-dependent star formation rate (SFR) stellar mass (M*) plane. Out of this sample 15 are identified as extreme oxygen emitters based on their [O III]/[O II] ratios (3.7%) and 104 galaxies have [O III]/[O II] > 1 (26%). Our analysis shows no significant correlation between M*, SFR, and the distance from the SFR M, relation with [O III]/[O II]. We find a decrease in the fraction of galaxies with [O III]/[O II] > 1 with increasing M*, however, this is most likely a result of the relationship between [O III]/[O II] and metallicity, rather than between [O III]/[O II] and M. We draw a comparison sample of local analogues with < z > 0.03 from the Sloan Digital Sky Survey, and find similar incidence rates for this sample. In order to investigate the evolution in the fraction of high [O III]/[O II] emitters with redshift, we bin the sample into three redshift subsamples of equal number, but find no evidence for a dependence on redshift. Furthermore, we compare the observed line ratios with those predicted by nebular models with no LyC escape and find that most of the extreme oxygen emitters can be reproduced by low metallicity models. The remaining galaxies are likely LyC emitter candidates.