Abstract
Feedback from active galactic nuclei (AGN) has become a major component in simulations of galaxy evolution, in particular for massive galaxies. AGN jets have been shown to provide a large ...amount of energy and are capable of quenching cooling flows. Their impact on the host galaxy, however, is still not understood. Subgrid models of AGN activity in a galaxy evolution context so far have been mostly focused on the quenching of star formation. To shed more light on the actual physics of the 'radio mode' part of AGN activity, we have performed simulations of the interaction of a powerful AGN jet with the massive gaseous disc () of a high-redshift galaxy. We spatially resolve both the jet and the clumpy, multi-phase interstellar medium (ISM) and include an explicit star formation model in the simulation. Following the system over more than 107 yr, we find that the jet activity excavates the central region, but overall causes a significant change to the shape of the density probability distribution function and hence the star formation rate due to the formation of a blast wave with strong compression and cooling in the ISM. This results in a ring- or disc-shaped population of young stars. At later times, the increase in star formation rate also occurs in the disc regions further out since the jet cocoon pressurizes the ISM. The total mass of the additionally formed stars may be up to for one duty cycle. We discuss the details of this jet-induced star formation (positive feedback) and its potential consequences for galaxy evolution and observable signatures.
Abstract
We present the results of a study investigating the dust attenuation law at z ≃ 5, based on synthetic spectral energy distributions (SEDs) calculated for a sample of N = 498 galaxies drawn ...from the First Billion Years (FiBY) simulation project. The simulated galaxies at z ≃ 5, which have M
1500 ≤ −18.0 and
$7.5 \le \rm {log(M/M}_{{\odot }}\rm {)} \le 10.2$
, display a mass-dependent α-enhancement, with a median value of
$\alpha /{\rm {Fe}}_{z=5} \simeq 4 \times \alpha /{\rm {Fe}}_{Z_{{\odot }}}$
. The median Fe/H ratio of the simulated galaxies is 0.14 ± 0.05 which produces steep intrinsic ultraviolet (UV) continuum slopes; 〈β
i
〉 = −2.4 ± 0.05. Using a set of simple dust attenuation models, in which the wavelength-dependent attenuation is assumed to be of the form A(λ) ∝ λ
n
, we explore the parameter values which best reproduce the observed z = 5 luminosity function (LF) and colour–magnitude relation (CMR). We find that a simple model in which the absolute UV attenuation is a linearly increasing function of log stellar mass (A
1500 = 0.5 × log(M/M⊙) − 3.3), and the dust attenuation slope (n) is within the range −0.7 ≤ n ≤ −0.3, can successfully reproduce the LF and CMR over a wide range of stellar population synthesis model assumptions, including the effects of massive binaries. This range of attenuation curves is consistent with a power-law fit to the Calzetti attenuation law in the UV (n = −0.55). In contrast, curves as steep as the Small Magellanic Cloud extinction curve (n = −1.24) are formally ruled out. Finally, we show that our models are consistent with recent 1.3 mm Atacama Large Millimeter Array observations of the Hubble Ultra Deep Field, and predict the form of the z ≃ 5 infrared excess (IRX)–β relation.
Abstract
We present the results of a new study of the relationship between infrared excess (IRX ≡ LIR/LUV), ultraviolet (UV) spectral slope (β) and stellar mass at redshifts 2 < z < 3, based on a ...deep Atacama Large Millimeter Array (ALMA) 1.3-mm continuum mosaic of the Hubble Ultra Deep Field. Excluding the most heavily obscured sources, we use a stacking analysis to show that z ≃ 2.5 star-forming galaxies in the mass range $9.25\le \log (M_{\ast }/\rm M_{{\odot }}) \le 10.75$ are fully consistent with the IRX–β relation expected for a relatively grey attenuation curve, similar to the commonly adopted Calzetti law. Based on a large, mass-complete sample of 2 ≤ z ≤ 3 star-forming galaxies drawn from multiple surveys, we proceed to derive a new empirical relationship between β and stellar mass, making it possible to predict UV attenuation (A1600) and IRX as a function of stellar mass, for any assumed attenuation law. Once again, we find that z ≃ 2.5 star-forming galaxies follow A1600–M* and IRX–M* relations consistent with a relatively grey attenuation law, and find no compelling evidence that star-forming galaxies at this epoch follow a reddening law as steep as the Small Magellanic Cloud (SMC) extinction curve. In fact, we use a simple simulation to demonstrate that previous determinations of the IRX–β relation may have been biased towards low values of IRX at red values of β, mimicking the signature expected for an SMC-like dust law. We show that this provides a plausible mechanism for reconciling apparently contradictory results in the literature and that, based on typical measurement uncertainties, stellar mass provides a cleaner prediction of UV attenuation than β. Although the situation at lower stellar masses remains uncertain, we conclude that for 2 < z < 3 star-forming galaxies with $\log (M_{\ast }/\rm M_{{\odot }}) \ge 9.75$, both the IRX–β and IRX–M* relations are well described by a Calzetti-like attenuation law.
ABSTRACT
Many astrophysical systems can only be accurately modelled when the behaviour of their baryonic gas components is well understood. The residual distribution (RD) family of partial ...differential equation (PDE) solvers produce approximate solutions to the corresponding fluid equations. We present a new implementation of the RD method. The solver efficiently calculates the evolution of the fluid, with up to second order accuracy in both time and space, across an unstructured triangulation, in both 2D and 3D. We implement a novel variable time stepping routine, which applies a drifting mechanism to greatly improve the computational efficiency of the method. We conduct extensive testing of the new implementation, demonstrating its innate ability to resolve complex fluid structures, even at very low resolution. We can resolve complex structures with as few as 3–5 resolution elements, demonstrated by Kelvin–Helmholtz and Sedov blast tests. We also note that we find cold cloud destruction time scales consistent with those predicted by a typical PPE solver, albeit the exact evolution shows small differences. The code includes three residual calculation modes, the LDA, N, and blended schemes, tailored for scenarios from smooth flows (LDA), to extreme shocks (N), and both (blended). We compare our RD solver results to state-of-the-art solvers used in other astrophysical codes, demonstrating the competitiveness of the new approach, particularly at low resolution. This is of particular interest in large scale astrophysical simulations, where important structures, such as star forming gas clouds, are often resolved by small numbers of fluid elements.
In order to specify cosmologically motivated initial conditions for major galaxy mergers (mass ratios ≤ 4:1) that are supposed to explain the formation of elliptical galaxies we study the orbital ...parameters of major mergers of cold dark matter halos using a high-resolution cosmological simulation. Almost half of all encounters are nearly parabolic with eccentricities $e \approx 1$ and no correlations between the halo spin planes or the orbital planes. The pericentric argument ω shows no correlation with the other orbital parameters and is distributed randomly. In addition we find that 50% of typical pericenter distances are larger than half the halo's virial radii which is much larger than typically assumed in numerical simulations of galaxy mergers. In contrast to the usual assumption made in semi-analytic models of galaxy formation the circularities of major mergers are found to be not randomly distributed but to peak around a value of $\epsilon \approx 0.5$. Additionally all results are independent of the minimum progenitor mass and major merger definitions (i.e. mass ratios ≤ 4:1; 3:1; 2:1).
We present a simulation of the cosmic evolution of the atomic and molecular phases of the cold hydrogen gas in about 3 X 107 galaxies, obtained by postprocessing the virtual galaxy catalog produced ...by De Lucia & Blaizot on the Millennium Simulation of cosmic structure. Our method uses a set of physical prescriptions to assign neutral atomic hydrogen (H I) and molecular hydrogen (H2) to galaxies, based on their total cold gas masses and a few additional galaxy properties. These prescriptions are specially designed for large cosmological simulations, where, given current computational limitations, individual galaxies can only be represented by simplistic model objects with a few global properties. Our recipes allow us to (1) split total cold gas masses between H I, H2, and helium, (2) assign realistic sizes to both the H I and H2 disks, and (3) evaluate the corresponding velocity profiles and shapes of the characteristic radio emission lines. The results presented in this paper include the local H I and H2 mass functions, the CO luminosity function, the cold gas mass-diameter relation, and the Tully-Fisher relation (TFR), which all match recent observational data from the local universe. We also present high-redshift predictions of cold gas diameters and the TFR, both of which appear to evolve markedly with redshift.
We investigate the role that dry mergers play in the build-up of massive galaxies within the cold dark matter paradigm. Implementing an empirical shut-off mass scale for star formation, we find a ...nearly constant dry merger rate of ∼6 × 10−5 Mpc−3 Gyr−1 at z≤ 1 and a steep decline at larger z. Less than half of these mergers are between two galaxies that are morphologically classified as early-types, and the other half is mostly between an early- and late-type galaxy. Latter are prime candidates for the origin of tidal features around red elliptical galaxies. The introduction of a transition mass scale for star formation has a strong impact on the evolution of galaxies, allowing them to grow above a characteristic mass scale of M*,c∼ 6.3 × 1010 M⊙ by mergers only. As a consequence of this transition, we find that around M*,c, the fraction of 1:1 mergers is enhanced with respect to unequal mass major mergers. This suggests that it is possible to detect the existence of a transition mass scale by measuring the relative contribution of equal mass mergers to unequal mass mergers as a function of galaxy mass. The evolution of the high-mass end of the luminosity function is mainly driven by dry mergers at low z. We however find that only 10–20 per cent of galaxies more massive than M*,c experience dry major mergers within their last Gyr at any given redshift z≤ 1.
A deep ALMA image of the Hubble Ultra Deep Field Dunlop, J. S; McLure, R. J; Biggs, A. D ...
Monthly notices of the Royal Astronomical Society,
04/2017, Letnik:
466, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Abstract
We present the results of the first, deep Atacama Large Millimeter Array (ALMA) imaging covering the full ≃4.5 arcmin2 of the Hubble Ultra Deep Field (HUDF) imaged with Wide Field Camera ...3/IR on HST. Using a 45-pointing mosaic, we have obtained a homogeneous 1.3-mm image reaching σ1.3 ≃ 35 μJy, at a resolution of ≃0.7 arcsec. From an initial list of ≃50 > 3.5σ peaks, a rigorous analysis confirms 16 sources with S
1.3 > 120 μJy. All of these have secure galaxy counterparts with robust redshifts (〈z〉 = 2.15). Due to the unparalleled supporting data, the physical properties of the ALMA sources are well constrained, including their stellar masses (M
*) and UV+FIR star formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-redshift Universe; indeed at z ≥ 2 our ALMA sample contains seven of the nine galaxies in the HUDF with M
* ≥ 2 × 1010 M⊙, and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S
1.3 ≃ 10 μJy. We find strong evidence for a steep star-forming ‘main sequence’ at z ≃ 2, with SFR ∝M
* and a mean specific SFR ≃ 2.2 Gyr−1. Moreover, we find that ≃85 per cent of total star formation at z ≃ 2 is enshrouded in dust, with ≃65 per cent of all star formation at this epoch occurring in high-mass galaxies (M
* > 2 × 1010 M⊙), for which the average obscured:unobscured SF ratio is ≃200. Finally, we revisit the cosmic evolution of SFR density; we find this peaks at z ≃ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured at z ≃ 4.
We investigate the stellar composition of bulges and elliptical galaxies as predicted by the cold dark matter paradigm using semi-analytical modelling. We argue that spheroid stars are built up of ...two main components, merger and quiescent, according to the origin of the stars. The merger component is formed during major mergers by gas driven to the centre, while the quiescent component is formed in gaseous discs and added later to the spheroid during major mergers. Galaxies more massive than M
C= 3 × 1010 M⊙ have on average only a 15 per cent merger component in their spheroids, while smaller galaxies can have up to 30 per cent. The merger component increases with redshift due to mergers involving more gas. However, we do not find mergers with gas fraction above ∼40 per cent of the remnants mass. Generally, the gas fraction is a decreasing function of the redshift at which the merger occurs and the mass of the remnant, with more massive remnants having smaller gas fraction and hence smaller merger components. This trend is independent of the environment of the galaxy with the only impact of the environment being that galaxies less massive than M
C have slightly larger merger components in dense environments. The fraction of stars in bulges for galaxies more massive than M
C is larger than 50 per cent. We find that the majority of stars in galaxies more massive than M
C reside within bulges and ellipticals independent of redshift and that the fraction increases with redshift. The most massive galaxies at each redshift are elliptical galaxies.