No Arabic abstract
The offset of high redshift star-forming galaxies in the [OIII]/H$beta$ versus [NII]/H$alpha$ (O3N2) diagram in comparison with the local star-forming galaxy sequence is now well established. The physical origin of the shift is the subject of some debate, and has important implications for metallicity measurements based on strong lines at all redshifts. To investigate the origin of the O3N2 offset, we use a sample of ~100,000 star-forming galaxies from SDSS DR12 to probe the empirical correlations between emission line diagnostics and measurable galaxy physical properties. In particular, we examine how surface density of star formation, ionization parameter, nitrogen-to-oxygen (N/O) ratio, and stellar mass drive position in two key diagnostic diagrams: O3N2 and [OIII]/H$beta$ versus [SII]/H$alpha$ (O3S2). We show that, at fixed [OIII]/H$beta$, galaxies falling closer to the high-redshift locus in O3N2 have higher star formation density, stellar mass and N/O ratios. We also find a tight correspondence in the distributions of stellar mass and N/O in the diagnostic diagrams. This relation, spanning a range of galaxy evolutionary states, suggests that the N/O-$M_{*}$ relation is more fundamental than the N/O-metallicity relation. We argue that a tight N/O-$M_{*}$ relation is well-motivated physically, and that the observed correlation of N/O with O/H in the local universe is primarily a reflection of the existence of the mass-metallicity relation. Because the mass-metallicity relation evolves much more rapidly with redshift than N/O-$M_{*}$, the N/O ratios of high redshift galaxies are significantly elevated in comparison with local galaxies with the same gas-phase O/H. The O3N2 shift and elevated N/O ratios observed in high redshift galaxies therefore come about as a natural consequence of the N/O-$M_{*}$ relation combined with the evolution of the mass-metallicity relation.
We investigate the Tully-Fisher Relation (TFR) for a morphologically and kine- matically diverse sample of galaxies from the SAMI Galaxy Survey using 2 dimensional spatially resolved Halpha velocity maps and find a well defined relation across the stellar mass range of 8.0 < log(M*) < 11.5. We use an adaptation of kinemetry to parametrise the kinematic Halpha asymmetry of all galaxies in the sample, and find a correlation between scatter (i.e. residuals off the TFR) and asymmetry. This effect is pronounced at low stellar mass, corresponding to the inverse relationship between stellar mass and kinematic asymmetry found in previous work. For galaxies with log(M*) < 9.5, 25 +/- 3% are scattered below the root mean square (RMS) of the TFR, whereas for galaxies with log(M*) > 9.5 the fraction is 10 +/- 1% We use simulated slits to directly compare our results with those from long slit spectroscopy and find that aligning slits with the photometric, rather than the kinematic, position angle, increases global scatter below the TFR. Further, kinematic asymmetry is correlated with misalignment between the photometric and kinematic position angles. This work demonstrates the value of 2D spatially resolved kinematics for accurate TFR studies; integral field spectroscopy reduces the underestimation of rotation velocity that can occur from slit positioning off the kinematic axis.
Emission line diagnostic diagrams probing the ionization sources in galaxies, such as the Baldwin-Phillips-Terlevich (BPT) diagram, have been used extensively to distinguish AGN from purely star-forming galaxies. Yet, they remain poorly understood at higher redshifts. We shed light on this issue with an empirical approach based on a z~0 reference sample built from ~300,000 SDSS galaxies, from which we mimic selection effects due to typical emission line detection limits at higher redshift. We combine this low-redshift reference sample with a simple prescription for luminosity evolution of the global galaxy population to predict the loci of high-redshift galaxies on the BPT and Mass-Excitation (MEx) diagnostic diagrams. The predicted bivariate distributions agree remarkably well with direct observations of galaxies out to z~1.5, including the observed stellar mass-metallicity (MZ) relation evolution. As a result, we infer that high-redshift star-forming galaxies are consistent with having normal ISM properties out to z~1.5, after accounting for selection effects and line luminosity evolution. Namely, their optical line ratios and gas-phase metallicities are comparable to that of low-redshift galaxies with equivalent emission-line luminosities. In contrast, AGN narrow-line regions may show a shift toward lower metallicities at higher redshift. While a physical evolution of the ISM conditions is not ruled out for purely star-forming galaxies, and may be more important starting at z>2, we find that reliably quantifying this evolution is hindered by selections effects. The recipes provided here may serve as a basis for future studies toward this goal. Code to predict the loci of galaxies on the BPT and MEx diagnostic diagrams, and the MZ relation as a function of emission line luminosity limits, is made publicly available.
The age distributions of stellar cluster populations have long been proposed to probe the recent formation history of the host galaxy. However, progress is hampered by the limited understanding of cluster disruption by evaporation and tidal shocks. We study the age distributions of clusters in smoothed particle hydrodynamics simulations of isolated disc galaxies, which include a self-consistent, physical model for the formation and dynamical evolution of the cluster population and account for the variation of cluster disruption in time and space. We show that the downward slope of the cluster age distribution due to disruption cannot be reproduced with a single functional form, because the disruption rate exhibits systematic trends with cluster age (the `cruel cradle effect). This problem is resolved by using the median cluster age to trace cluster disruption. Across 120 independent galaxy snapshots and simulated cluster populations, we perform two-dimensional power law fits of the median cluster age to various macroscopic physical quantities and find that it scales as $t_{rm med}propto Sigma^{-0.51pm0.03}sigma_{rm 1D}^{-0.85pm0.10}M_{rm min}^gamma$, for the gas surface density $Sigma$, gas velocity dispersion $sigma_{rm 1D}$, and minimum cluster mass $M_{rm min}$. This scaling accurately describes observed cluster populations and indicates disruption by impulsive tidal shocks from the interstellar medium. The term $M_{rm min}^gamma$ provides a model-independent way to measure the mass dependence of the cluster disruption time $gamma$. Finally, the ensemble-average cluster lifetime depends on the gas density less strongly than the instantaneous disruption time of single clusters. These results reflect the variation of cluster disruption in time and space. We provide quantitative ways of accounting for these physics in cluster population studies.
We have obtained structural parameters of about 340,000 galaxies from the Kilo Degree Survey (KiDS) in 153 square degrees of data release 1, 2 and 3. We have performed a seeing convolved 2D single Sersic fit to the galaxy images in the 4 photometric bands (u, g, r, i) observed by KiDS, by selecting high signal-to-noise ratio (S/N > 50) systems in every bands. We have classified galaxies as spheroids and disc-dominated by combining their spectral energy distribution properties and their Sersic index. Using photometric redshifts derived from a machine learning technique, we have determined the evolution of the effective radius, Re and stellar mass, mst, versus redshift, for both mass complete samples of spheroids and disc-dominated galaxies up to z ~ 0.6. Our results show a significant evolution of the structural quantities at intermediate redshift for the massive spheroids ($mbox{Log} M_*/M_odot>11$, Chabrier IMF), while almost no evolution has found for less massive ones ($mbox{Log} M_*/M_odot < 11$). On the other hand, disc dominated systems show a milder evolution in the less massive systems ($mbox{Log} M_*/M_odot < 11$) and possibly no evolution of the more massive systems. These trends are generally consistent with predictions from hydrodynamical simulations and independent datasets out to redshift z ~ 0.6, although in some cases the scatter of the data is large to drive final conclusions. These results, based on 1/10 of the expected KiDS area, reinforce precedent finding based on smaller statistical samples and show the route toward more accurate results, expected with the the next survey releases.
Using observations made with MOSFIRE on Keck I as part of the ZFIRE survey, we present the stellar mass Tully-Fisher relation at 2.0 < z < 2.5. The sample was drawn from a stellar mass limited, Ks-band selected catalog from ZFOURGE over the CANDELS area in the COSMOS field. We model the shear of the Halpha emission line to derive rotational velocities at 2.2X the scale radius of an exponential disk (V2.2). We correct for the blurring effect of a two-dimensional PSF and the fact that the MOSFIRE PSF is better approximated by a Moffat than a Gaussian, which is more typically assumed for natural seeing. We find for the Tully-Fisher relation at 2.0 < z < 2.5 that logV2.2 =(2.18 +/- 0.051)+(0.193 +/- 0.108)(logM/Msun - 10) and infer an evolution of the zeropoint of Delta M/Msun = -0.25 +/- 0.16 dex or Delta M/Msun = -0.39 +/- 0.21 dex compared to z = 0 when adopting a fixed slope of 0.29 or 1/4.5, respectively. We also derive the alternative kinematic estimator S0.5, with a best-fit relation logS0.5 =(2.06 +/- 0.032)+(0.211 +/- 0.086)(logM/Msun - 10), and infer an evolution of Delta M/Msun= -0.45 +/- 0.13 dex compared to z < 1.2 if we adopt a fixed slope. We investigate and review various systematics, ranging from PSF effects, projection effects, systematics related to stellar mass derivation, selection biases and slope. We find that discrepancies between the various literature values are reduced when taking these into account. Our observations correspond well with the gradual evolution predicted by semi-analytic models.