Simulating the formation and evolution of galaxies: Multi-phase description of the interstellar medium, star formation, and energy feedback
TL;DR: In this paper, a multi-phase representation of the ISM in NB-TSPH simulations of galaxy formation and evolution with particular attention to the case of early-type galaxies is presented.
Abstract: We present a multi-phase representation of the ISM in NB-TSPH simulations of galaxy formation and evolution with particular attention to the case of early-type galaxies. Cold gas clouds are described by the so-called sticky particles algorithm. They can freely move throughout the hot ISM medium; stars form within these clouds and the mass exchange among the three baryonic phases (hot gas, cold clouds, stars) is governed by radiative and Compton cooling and energy feedback by supernova (SN) explosions, stellar winds, and UV radiation. We also consider thermal conduction, cloud-cloud collisions, and chemical enrichment. Our model agrees with and improves upon previous studies on the same subject. The results for the star formation rate are very promising and agree with recent observational data on early-type galaxies. These models lend further support to the revised monolithic scheme of galaxy formation, which has recently been also strengthened by high redshift data leading to the so-called downsizing and top-down scenarios.
Citations
More filters
TL;DR: In this paper , the authors derived the mutual relationship among the stellar mass, effective radius, velocity dispersion, and luminosity of ETGs as a function of the ϵ-factor and calculated the coefficients of the FP.
Abstract: We have recently suggested that the combination of the scalar virial theorem ($M_s \sim R_e \sigma^2$) and the $L=L'_0 \sigma^\beta(t)$ law, with L'_0 and $\beta$ changing from galaxy to galaxy (and with time), can provide a new set of equations valid for investigating the evolution of early-type galaxies (ETGs) (Donofrio&Chiosi, 2022). These equations are able to account for the tilt of the Fundamental Plane (FP) and to explain the observed distributions of ETGs in all its projections. In this paper we analyze the advantages offered by those equations, derive the $\beta$ and $L'_0$ parameters for real and simulated galaxies, and demonstrate that, according to the value of $\beta$, galaxies can move only along some permitted directions in the FP projections. Then, we show that simple galaxy models that grow in mass by infall of gas and form stars with a star formation rate depending on the stellar velocity dispersion nicely reproduce the observed distributions of ETGs in the FP projections and yield $\beta$s that agree with the measured ones. We derive the mutual relationships among the stellar mass, effective radius, velocity dispersion, and luminosity of ETGs as a function of $\beta$ and calculate the coefficients of the FP. Then, using the simple infall models, we show that the star formation history of ETGs is compatible with the $\sigma$-dependent star formation rate, and that both positive and negative values of $\beta$ are possible in a standard theory of galaxy evolution. The parameter $\beta(t)$ offers a new view of the evolution of ETGs. In brief, i) it gives a coherent interpretation of the FP and of the motions of galaxies in its projections; ii) it is the fingerprint of their evolution; iii) it measures the degree of virialization of ETGs; iv) and finally it allows us to infer their evolution in the near past.
1 citations
References
More filters
TL;DR: In this paper, the evolutionary significance of the observed luminosity function for main-sequence stars in the solar neighborhood is discussed and it is shown that stars move off the main sequence after burning about 10 per cent of their hydrogen mass and that stars have been created at a uniform rate in a solar neighborhood for the last five billion years.
Abstract: The evolutionary significance of the observed luminosity function for main-sequence stars in the solar neighborhood is discussed. The hypothesis is made that stars move off the main sequence after burning about 10 per cent of their hydrogen mass and that stars have been created at a uniform rate in the solar neighborhood for the last five billion years. Using this hypothesis and the observed luminosity function, the rate of star creation as a function of stellar mass is calculated. The total number and mass of stars which have moved off the main sequence is found to be comparable with the total number of white dwarfs and with the total mass of all fainter main-sequence stars, respectively.
8,607 citations
TL;DR: GADGET-2 as mentioned in this paper is a massively parallel tree-SPH code, capable of following a collisionless fluid with the N-body method, and an ideal gas by means of smoothed particle hydrodynamics.
Abstract: We discuss the cosmological simulation code GADGET-2, a new massively parallel TreeSPH code, capable of following a collisionless fluid with the N-body method, and an ideal gas by means of smoothed particle hydrodynamics (SPH). Our implementation of SPH manifestly conserves energy and entropy in regions free of dissipation, while allowing for fully adaptive smoothing lengths. Gravitational forces are computed with a hierarchical multipole expansion, which can optionally be applied in the form of a TreePM algorithm, where only short-range forces are computed with the ‘tree’ method while long-range forces are determined with Fourier techniques. Time integration is based on a quasi-symplectic scheme where long-range and short-range forces can be integrated with different time-steps. Individual and adaptive short-range time-steps may also be employed. The domain decomposition used in the parallelization algorithm is based on a space-filling curve, resulting in high flexibility and tree force errors that do not depend on the way the domains are cut. The code is efficient in terms of memory consumption and required communication bandwidth. It has been used to compute the first cosmological N-body simulation with more than 10 10 dark matter particles, reaching a homogeneous spatial dynamic range of 10 5 per dimension in a three-dimensional box. It has also been used to carry out very large cosmological SPH simulations that account for radiative cooling and star formation, reaching total particle numbers of more than 250 million. We present the algorithms used by the code and discuss their accuracy and performance using a number of test problems. GADGET-2 is publicly released to the research community. Ke yw ords: methods: numerical ‐ galaxies: interactions ‐ dark matter.
6,196 citations
TL;DR: A finite-size particle scheme for the numerical solution of two-and three-dimensional gas dynamical problems of astronomical interest is described and tested in this article, which is then applied to the fission problem for optically thick protostars.
Abstract: A finite-size particle scheme for the numerical solution of two- and three-dimensional gas dynamical problems of astronomical interest is described and tested. The scheme is then applied to the fission problem for optically thick protostars. Results are given, showing the evolution of one such protostar from an initial state as a single, rotating star to a final state as a triple system whose components contain 60% of the original mass. The decisiveness of this numerical test of the fission hypothesis and its relevance to observed binaries are briefly discussed.
5,508 citations
TL;DR: A novel method of directly calculating the force on N bodies that grows only as N log N is described, using a tree-structured hierarchical subdivision of space into cubic cells, each is recursively divided into eight subcells whenever more than one particle is found to occupy the same cell.
Abstract: Until recently the gravitational N-body problem has been modelled numerically either by direct integration, in which the computation needed increases as N2, or by an iterative potential method in which the number of operations grows as N log N. Here we describe a novel method of directly calculating the force on N bodies that grows only as N log N. The technique uses a tree-structured hierarchical subdivision of space into cubic cells, each of which is recursively divided into eight subcells whenever more than one particle is found to occupy the same cell. This tree is constructed anew at every time step, avoiding ambiguity and tangling. Advantages over potential-solving codes are: accurate local interactions; freedom from geometrical assumptions and restrictions; and applicability to a wide class of systems, including (proto-)planetary, stellar, galactic and cosmological ones. Advantages over previous hierarchical tree-codes include simplicity and the possibility of rigorous analysis of error. Although we concentrate here on stellar dynamical applications, our techniques of efficiently handling a large number of long-range interactions and concentrating computational effort where most needed have potential applications in other areas of astrophysics as well.
3,750 citations
TL;DR: The first extensive catalog of galactic embedded clusters is compiled, finding that the embedded cluster birthrate exceeds that of visible open clusters by an order of magnitude or more indicating a high infant mortality rate for protocluster systems.
Abstract: ▪ Abstract Stellar clusters are born embedded within giant molecular clouds (GMCs) and during their formation and early evolution are often only visible at infrared wavelengths, being heavily obscured by dust. Over the past 15 years advances in infrared detection capabilities have enabled the first systematic studies of embedded clusters in galactic molecular clouds. In this article we review the current state of empirical knowledge concerning these extremely young protocluster systems. From a survey of the literature we compile the first extensive catalog of galactic embedded clusters. We use the catalog to construct the mass function and estimate the birthrate for embedded clusters within ∼2 kpc of the sun. We find that the embedded cluster birthrate exceeds that of visible open clusters by an order of magnitude or more indicating a high infant mortality rate for protocluster systems. Less than 4–7% of embedded clusters survive emergence from molecular clouds to become bound clusters of Pleiades age. Th...
2,949 citations