History of Star Formation and Chemical Enrichment in the Milky Way Disk

Based on a physical treatment of the star formation law similar to that given by Efstathiou, we have improved our two-component chemical evolution model for the Milky Way disk. Two gas infall rates are compared, one exponential, one Gaussian. It is shown that the star formation law adopted in this paper depends more strongly on the gas surface density than that in Chang et al. It has large effects on the history of star formation and gas evolution of the whole disk. In the solar neighborhood, the history of chemical evolution and star formation is not sensitive to whether the infall rate is Gaussian or exponential. For the same infall time scale, both forms predict the same behavior for the current properties of the Galactic disk. The model predictions do depend on whether or not the infall time scale varies with the radius, but current available observations cannot decide which case is the more realistic. Our results also show that it would be inadequate to describe the gradient evolution along the Gala     

银盘的径向金融丰度梯度

详细综棕了银盘（包括HII区，早B型星，行星状星云和疏散星团）径向元素丰度梯度的观测结果，分析了丰度梯度的空间和时间变化的情况，指出根据目前的观测结果，还很难确定在银盘的演化历史中径向元素丰度梯度是逐渐变平缓还是逐渐变陡，比较了目前各种化学演化模型对径向丰度梯度演化的预测结果，初步探讨了丰度梯度可能的产生机制及影响其演化的各种重要物理过程。     

Infall of gas in galaxies and triggered star formation

While the importance of merging, accretion, and infall processesin determining galactic evolution is well established boththeoretically and observationally, details on how such processesare taking place nowadays even in our own Galaxy are stillrelatively poorly known, especially due to large remaininguncertainties on the location and origin of high velocity clouds.In this paper we focus on the possible role that galacticoutflows and gas infall may have on directly triggering starformation in the halo and in galactic disks. While compellingevidence has been accumulating in recent years suggesting thatsome level of star formation directly triggered by outflows isvery likely to exist in the halo of some galaxies, the evidencefor star formation dynamically triggered by infall is far moreelusive due to confusion with other, more efficient large-scalestar forming mechanisms operating in the galactic disk. Despite ofincreasingly realistic simulations of the gas circulation betwenthe gas and the halo and of high velocity cloud impacts ongalactic disks, the efficiency of star formation directlytriggered by such impacts remains an open question.     

Evolution of the Radial Abundance Gradient and Cold Gas of the Galactic Disk

In order to understand the forming mechanism of the radial abun- dance gradient of the Galactic disk and the evolution of cold gas, we have con- structed a chemical evolution model of the Galactic disk, in which the star for- mation law concerned with molecular hydrogens is adopted, and the evolution of mass surface density is calculated for the molecular and atomic hydrogens separately, then the model predictions and the observed radial distributions of some physical quantities are compared. The result indicates that the model prediction is sensitive to the adopted infall timescale, the model which adopts the star formation law concerned with the molecular hydrogens can agree well with the major observed properties of the Galactic disk, especially can obtain naturally the radial oxygen abundance gradient of the Galactic disk, and the radial surface density profile of cold gas. The assumption of instantaneous or non-instantaneous recycling approximation has a small effect on the evolution of cold gas, especially in the case of rather low gas density.     

Chemical evolution models of local dwarf spheroidal galaxies

We calculate chemical evolution models for four dwarf spheroidal (dSph) satellites of the Milky Way (Carina, Ursa Minor, Leo I and Leo II) for which reliable non-parametric star formation histories have been derived. In this way, the independently-obtained star formation histories are used to constrain the evolution of the systems we are treating. This allows us to obtain robust inferences on the history of such crucial parameters of galactic evolution as gas infall, gas outflows and global metallicities for these systems. We can then trace the metallicity and abundance ratios of the stars formed, the gas present at any time within the systems and the details of gas ejection, of relevance to enrichment of the galaxies environment. We find that galaxies showing one single burst of star formation (Ursa Minor and Leo II) require a dark halo slightly larger that the current estimates for their tidal radii, or the presence of a metal-rich selective wind that might carry away much of the energy output of their supernovae before this might have interacted and heated the gas content, for the gas to be retained until the observed stellar populations have formed. Systems showing extended star formation histories (Carina and Leo I), however, are consistent with the idea that their tidally-limited dark haloes provide the necessary gravitational potential wells to retain their gas. The complex time structure of the star formation in these systems remains difficult to understand. Observations of detailed abundance ratios for Ursa Minor strongly suggest that the star formation history of this galaxy might in fact resemble the complex picture presented by Carina or Leo I, but localized at a very early epoch.     

The G-dwarf problem in the solar neighbourhood: a Lagrangian picture

Recent data on the empirical metallicity distribution of G dwarfs in the disk solar neighbourhood are fitted in two different ways. We use an extended Poisson distribution in the limit where the probability of star formation is small, and a Gauss distribution in the limit where a large number of physical variables is required to determine stellar metal abundance. Both are found to reproduce the data at the same (acceptable) extent, with a slight preference for the former. The emprirical, differential metallicity distribution of G dwarfs in the disk solar neighbourhood is compared with its theoretical counterpart, in the picture of a closed, comoving model of chemical evolution. The limits of the currently used infall models are discussed and a scenario of galactic formation and evolution is presented. The Galactic history is thought as made of two main phases: contraction (which produces the extended component) and equilibrium (which gives the disk). In this view, the stars observed within the solar cylinder did not necessarily arise from the primordial gas which later collapsed into the disk solar neighbourhood. It is found that the G-dwarf problem is strongly alleviated, with the possible exception of the low-metallicity and high-metallicity tail of the distribution. The best choice of parameters implies: (i) a metal yield in the contraction phase which is larger by a factor of about 5 with respect to the equilibrium phase; (ii) a model halo mass fraction of about 0.3; (iii) a model disk mass fraction of about 0.6. It provides additional support to the idea of a generalized Schmidt star formation law, which is different in different phases of evolution. The model, cumulative, G-dwarf metallicity distribution in the disk solar neighbourhood is found to predict too may low-metallicity stars with respect to its empirical counterpart, related to a Poissonian or Gaussian fit. The main resons for the occurrence of a G-dwarf problem are discussed. Finally, a stochastic process of star formation, related to a Poisson distribution, is briefly outlined.     

Numerical simulations of the formation and chemical evolution of galaxies

We present the results of a set of three-dimensional SPH-Treecode simulations which model the formation and early evolution of disc galaxies, including the generation and return of heavy elements to the interstellar medium by star formation. Starting from simple initial conditions which are given by a uniform density sphere of gas which is embedded in a dark matter halo and in solid-body rotation, we are able to form realistic disc galaxies, and find that an exponential gas disc is quickly formed. Star formation within this exponential disc naturally leads to the formation of abundance gradients which are in broad agreement with those observed, although they are slightly shallower than some observations.
We investigate the systematic effects of variation of mass, rotation and star formation parameters on the abundance gradients. We find that the abundance gradients are most sensitive to changes in the star formation parameters or rotation. Including a critical-density cut-off in the star formation law causes abundance gradients to be steepened.
Analysis of gas flows within the models shows radial flows which are a function of angle of azimuth around the galaxies, with alternating inward and outward flows. This motion is linked to the presence of a bar, whose strength is related to the amount of star formation in the models, and there is a gentle drift of mass inwards. The shallow abundance gradients may be linked to these radial flows.     

The G-dwarf problem in the Galactic spheroid

This paper has two parts: one about observational constraints, and the other about chemical evolution models. In the first part, the empirical differential metallicity distribution (EDMD) is deduced from two different samples involving (i) 268 K-giant bulge stars [Sadler, E.M., Rich, R.M., Terndrup, D.M., 1996. AJ 112, 171], and (ii) 149 globular clusters [Mackey, A.D., van den Bergh, S., 2005. MNRAS 360, 631], in addition to previous results (Caimmi, R., 2001b, AN 322, 241 (C01)) related to (iii) 372 solar neighbourhood halo subdwarfs [Ryan, S.G., Norris, J.E., 1991. AJ 101, 1865]. Under the assumption that each distribution is typical for the corresponding subsystem, the EDMD of the Galactic spheroid is determined by weighting the mass. The empirical age-metallicity relation (EAMR) involving absolute ages is deduced from recent results related to a homogeneous sample of globular clusters [De Angeli, F., Piotto, G., Cassisi, S., et al., 2005. AJ 130, 116]. In the second part, models of chemical evolution for the Galactic halo and bulge are computed, assuming the instantaneous recycling approximation. The EDMD data are fitted, to an acceptable extent, by simple models of chemical evolution implying both homogeneous and inhomogeneous star formation, provided that star formation is inhibited during halo formation and enhanced during bulge formation, with respect to the disk solar neighbourhood, taken to be representative of the whole disk. The initial mass function (IMF) is assumed to be a universal power law, which implies the same value of the true yield in different subsystems. The theoretical differential metallicity distribution (TDMD) is first determined for the halo and the bulge separately, and then for the Galactic spheroid by weighting the mass. The EAMR cannot be fitted into the Simple model that implies homogeneous star formation, but shows a non-monotonic trend characterized by large dispersion. On the other hand, simple models involving inhomogeneous star formation yield a theoretical age-metallicity relation (TAMR) which reproduces the data to an acceptable extent. For gas ouflow from the proto-halo, acceptable models give rise to different predictions in different alternatives. If the Galactic spheroid and disk underwent decoupled chemical evolution, i.e. no gas exchange between the related reservoirs, less than one third of the bulge mass outflowed from the proto-halo. If the Galactic spheroid and disk underwent coupled chemical evolution, i.e. some gas exchange between the related reservoirs, the existence of an unseen baryonic halo (or equivalent amount of gas lost by the Galaxy) with mass comparable to bulge mass, is necessarily needed. In this view, the outflowing proto-halo gas which remains bound to the Galaxy, produces both the bulge and the disk.     

A grid of chemical evolution models as a tool to interpret spiral and irregular galaxies data

We present a generalization of the multiphase chemical evolution model (CEM) applied to a wide set of theoretical galaxies with different masses and evolutionary rates. This generalized set of models has been computed using the so-called universal rotation curve from Persic, Salucci & Steel to calculate the radial mass distribution of 44 theoretical protogalaxies. This distribution is a fundamental input which, besides its own effect on the galaxy evolution, defines the characteristic collapse time-scale or gas infall rate on to the disc. We have adopted 10 sets of values, between 0 and 1, for the molecular cloud and star formation efficiencies, as corresponding to their probability nature, for each one of the radial distributions of total mass. Thus, we have constructed a biparametric grid of models, depending on those efficiency sets and on the rotation velocity, whose results are valid in principle for any spiral or irregular galaxy. The model results provide the time-evolution of different regions of the disc and the halo along galactocentric distance, measured by the gas (atomic and molecular) and stellar masses, the star formation rate (SFR) and chemical abundances of 14 elements, for a total of 440 models. This grid may be used to estimate the evolution of a given galaxy for which only present time information, such as radial distributions of elemental abundances, gas densities and/or star formation, which are the usual observational constraints of chemical evolution models (CEMs), is available.     

Oxygen abundance in local disk and bulge: chemical evolution with a strictly universal IMF

This paper has two parts: one about observational constraints related to the empirical differential oxygen abundance distribution (EDOD), and the other about inhomogeneous models of chemical evolution, in particular the theoretical differential oxygen abundance distribution (TDOD). In the first part, the EDOD is deduced from subsamples related to two different samples involving (i) N=532 solar neighbourhood (SN) stars within the range, −1.5<[Fe/H]<0.5, for which the oxygen abundance has been determined both in presence and in absence of the local thermodynamical equilibrium (LTE) approximation (Ramirez et al. in Astron. Astrophys. 465:271, 2007); and (ii) N=64 SN thick disk, SN thin disk, and bulge K-giant stars within the range, −1.7<[Fe/H]<0.5, for which the oxygen abundance has been determined (Melendez et al. in Astron. Astrophys. 484:L21, 2008). A comparison is made with previous results implying use of [O/H]–[Fe/H] empirical relations (Caimmi in Astron. Nachr. 322:241, 2001b; New Astron. 12:289, 2007) related to (iii) 372 SN halo subdwarfs (Ryan and Norris in Astron. J. 101:1865, 1991); and (iv) 268 K-giant bulge stars (Sadler et al. in Astron. J. 112:171, 1996). The EDOD of the SN thick + thin disk is determined by weighting the mass, for assumed SN thick to thin disk mass ratio within the range, 0.1–0.9. In the second part, inhomogeneous models of chemical evolution for the SN thick disk, the SN thin disk, the SN thick + thin disk, the SN halo, and the bulge, are computed assuming the instantaneous recycling approximation. The EDOD data are fitted, to an acceptable extent, by their TDOD counterparts with the exception of the thin or thick + thin disk, where two additional restrictions are needed: (i) still undetected, low-oxygen abundance thin disk stars exist, and (ii) a single oxygen overabundant star is removed from a thin disk subsample. In any case, the (assumed power-law) stellar initial mass function (IMF) is universal but gas can be inhibited from, or enhanced in, forming stars at different rates with respect to a selected reference case. Models involving a strictly universal IMF (i.e. gas neither inhibited from, nor enhanced in, forming stars with respect to a selected reference case) can also reproduce the data to an acceptable extent. Our main conclusions are (1) different models are necessary to fit the (incomplete) halo sample, which is consistent with the idea of two distinct halo components: an inner, flattened halo in slow prograde rotation, and an outer, spherical halo in net retrograde rotation (Carollo et al. in Nature 450:1020, 2007); (2) the oxygen enrichment within the inner SN halo, the SN thick disk, and the bulge, was similar and coeval within the same metallicity range, as inferred from observations (Prochaska et al. in Astron. J. 120:2513, 2000); (3) the fit to thin disk data implies an oxygen abundance range similar to its thick disk counterpart, with the extension of conclusion (2) to the thin disk, and the evolution of the thick + thin disk as a whole (Haywood in Mon. Not. R. Astron. Soc. 388:1175, 2008) cannot be excluded; (4) leaving outside the outer halo, a fit to the data related to different environments is provided by models with a strictly universal IMF but different probabilities of a region being active, which implies different global efficiencies of the star formation rate; (5) a special case of stellar migration across the disk can be described by models with enhanced star formation, where a fraction of currently observed SN stars were born in situ and a comparable fraction is due to the net effect of stellar migration, according to recent results based on high-resolution N-body + smooth particle hydrodynamics simulations (Roškar et al. in Astrophys. J. Lett. 684:L79, 2008).     

Chemical Evolution of the Magellanic Clouds

Previous models for the chemical evolution of the Magellanic Clouds have assumed either a steepened IMF compared to the solar neighbourhood or preferential expulsion of oxygen and α-particle elements by selective galactic winds. These assumptions were largely motivated by a belief that the O/Fe ratio in the Clouds is substantially lower than in the Galaxy, but the difference appears to have been exaggerated: Galactic supergiants have a similar O/Fe ratio as Cloud supergiants, there is no corresponding effect in Mg and other α-elements and a combination of data from planetary nebulae, H II regions and supernova remnants indicates an O/Fe ratio more or less equal to solar. Consequently new analytical models for the chemical evolution of the Magellanic Clouds have been developed, assuming chemical yields and time delays identical to those we previously assumed for the solar neighbourhood, but assuming (in addition to infall) non-selective galactic winds and burst-like modes of star formation represented by discontinuous variations in the star formation rate per unit gas mass. We find adequate agreement with age-metallicity relations and element:element ratios within their substantial uncertainties, whereas our LMC model turns out to give an excellent fit to the anomalous Galactic halo stars discovered by Nissen and Schuster (1997). It also gives an enhanced SNIa/SNII ratio compared to the solar neighbourhood, due to the assumption that the SFR has declined in the past 1 to 2 Gyr. This revised version was published online in July 2006 with corrections to the Cover Date.     

Cold gas accretion in galaxies

Evidence for the accretion of cold gas in galaxies has been rapidly accumulating in the past years. HI observations of galaxies and their environment have brought to light new facts and phenomena which are evidence of ongoing or recent accretion: (1) A large number of galaxies are accompanied by gas-rich dwarfs or are surrounded by HI cloud complexes, tails and filaments. This suggests ongoing minor mergers and recent arrival of external gas. It may be regarded, therefore, as direct evidence of cold gas accretion in the local universe. It is probably the same kind of phenomenon of material infall as the stellar streams observed in the halos of our galaxy and M 31. (2) Considerable amounts of extra-planar HI have been found in nearby spiral galaxies. While a large fraction of this gas is undoubtedly produced by galactic fountains, it is likely that a part of it is of extragalactic origin. Also the Milky Way has extra-planar gas complexes: the Intermediate- and High-Velocity Clouds (IVCs and HVCs). (3) Spirals are known to have extended and warped outer layers of HI. It is not clear how these have formed, and how and for how long the warps can be sustained. Gas infall has been proposed as the origin. (4) The majority of galactic disks are lopsided in their morphology as well as in their kinematics. Also here recent accretion has been advocated as a possible cause. In our view, accretion takes place both through the arrival and merging of gas-rich satellites and through gas infall from the intergalactic medium (IGM). The new gas could be added to the halo or be deposited in the outer parts of galaxies and form reservoirs for replenishing the inner parts and feeding star formation. The infall may have observable effects on the disk such as bursts of star formation and lopsidedness. We infer a mean “visible” accretion rate of cold gas in galaxies of at least . In order to reach the accretion rates needed to sustain the observed star formation (), additional infall of large amounts of gas from the IGM seems to be required.     

星系中原子和分子气体的半解析模型

星际气体是星系中重子物质的重要组成部分,其中的分子气体（主要是分子氢H₂）以及原子气体（主要是中性氢HI）对于星系中发生的各个物理过程至关重要。本文在前人的星系形成和演化的半解析模型基础上,加入了描述星系盘中分子气体和原子气体成分的物理模型,来研究分子气体和原子气体对于星系形成和演化所起的作用。我们主要使用了马普天体物理所Munich Group的L-Galaxies半解析星系形成模型,并借鉴了星系化学演化模型的方法,把半解析模型中的每一个星系盘分成了多个同心圆圈,然后在每个圈中分别追踪气体下落、分子气体和原子气体转化、恒星形成、金属增丰、超新星爆发加热冷气体等发生在星系盘上的物理过程,并且每个同心圈都是独立演化的。在我们的模型中,一个基本假设是每个时间步内气体都是以指数形式下落到星系盘上,并且直接叠加在已有的气体径向面密度轮廓之上,其中指数盘的标长r_d正比于星系所在暗物质晕的维里半径r_vir与旋转参量λ的乘积。我们的模型使用了两种描述分子气体形成的模型:一种是基于Krumholz等人解析模型的结果,其中分子气体的比例与局域气体面密度以及局域气体金属丰度相关;另一种是分子气体比例与星际压强相关的模型,根据Obreschkow等人的近似,分子气体的比例与气体面密度以及恒星质量面密度相关。由于恒星形成过程发生在星际巨分子云之中,并且根据Leroy等人的观测结果,恒星形成率面密度近似正比于分子气体的面密度,因此我们在模型中使用了与分子气体面密度相关的恒星形成规律。     

The chemo-dynamical evolution of a disk galaxy

I present a model for the formation and evolution of a massive disk galaxy, within a growing dark halo whose mass evolves according to cosmological simulations of structure formation. The galactic evolution is simulated with a new three-dimensional chemo-dynamical code, including dark matter, stars and a multi-phase ISM. We follow the evolution from redshift z= 4.85 until the present epoch. The energy release by massive stars and supernovae prevents a rapid collapse of the baryonic matter and delays the maximum star formation until redshift z ≈ 1. The galaxy forms radially from inside-out and vertically from top-to-bottom. Correspondingly, the inner halo is the oldest component, followed by the outer halo, the bar/bulge, the thick and the thin disk. The bulge in the model consists of at least two stellar subpopulations, an early collapse population and a population that formed later in the bar. This revised version was published online in August 2006 with corrections to the Cover Date.     

Spectroscopic evidence for protostellar infall

The observational evidence for infall associated with star formation is discussed. Whilst spectral energy distributions of young protostellar objects are consistent with infall, the best direct evidence comes from millimetre and sub-millimetre spectral line observations. Considerations of the formation of the line profiles and the chemical effects of gas-grain interactions suggest that there is only a very short ‘window’ in the evolutionary track of a protostellar object during which infall is directly observable. This may explain why so few infall candidates have been detected. It is argued that self-consistent models of the dynamical and chemical evolution of collapsing cores, coupled to multiple high resolution line observations, will provide definitive evidence for the presence of infall in these objects.     

Galaxy formation and evolution: low-surface-brightness galaxies

We investigate in detail the hypothesis that low-surface-brightness galaxies (LSBs) differ from ordinary galaxies simply because they form in haloes with large spin parameters. We compute star formation rates using the Schmidt law, assuming the same gas infall dependence on surface density as used in models of the Milky Way. We build stellar population models, predicting colours, spectra and chemical abundances. We compare our predictions with observed values of metallicity and colours for LSBs, and find excellent agreement with all observables. In particular, integrated colours, colour gradients, surface brightness and metallicity match very well to the observed values of LSBs for models with ages larger than 7 Gyr and high values (λ > 0.05) for the spin parameter of the haloes. We also compute the global star formation rate (SFR) in the Universe due to LSBs, and show that it has a flatter evolution with redshift than the corresponding SFR for normal discs. We furthermore compare the evolution in redshift of [ Zn / H ] for our models to those observed in damped Lyman α systems by Pettini et al. and show that damped Lyman α system abundances are consistent with the predicted abundances at different radii for LSBs. Finally, we show how the required late redshift of collapse of the halo may constrain the power spectrum of fluctuations.     

The G‐dwarf problem in the solar neighbourhood: II. Halo subdwarfs

The oxygen abundance distribution in solar neighbourhood halo subdwarfs is deduced, using two alternative, known empirical relationships, involving the presence or the absence of [O/Fe] plateau for low [Fe/H] values, from a sample of 372 kinematically selected halo stars, for which the iron abundance distribution has been determined by Ryan & Norris (1991). The data are interpreted by a simple, either homogeneous or inhomogeneous model of chemical evolution, using an updated value of the solar oxygen abundance. The effect of changing the solar oxygen abundance, the power‐law exponent in the initial mass function, and the rate of oxygen nucleosyntesis, keeping the remaining input parameters unchanged, is investigated, and a theorem is stated. In all cases, part of the gas must necessarily be inhibited from forming stars, and no disk contamination has to be advocated for fitting the empirical oxygen abundance distribution in halo subdwarfs of the solar neighbourhood (EGD). Then a theorem is stated, which allows a one‐to‐one correspondence between simple, homogeneous models with and without inhibited gas, related to same independent parameters of chemical evolution, except lower stellar mass limit, real yield, and inhibition parameter. The mutual correlations between the latter parameters are also specified. In addition the starting point, and the point related to the first step, of the theoretical distribution of oxygen abundance (TGD) predicted by simple, inhomogeneous models, is calculated analytically. The mean oxygen abundance of the total and only inhibited gas, respectively, are also determined. Following the idea of a universal, initial mass function (IMF), a power‐law with both an exponent p = 2.9, which is acceptably close to Scalo IMF for m ≳ m_⊙, and an exponent p = 2.35, i.e. Salpeter IMF, have been considered. In general, both the age‐metallicity relationship and the empirical distribution of oxygen abundance in G dwarfs of the disk solar neighbourhood, are fitted by power‐law IMF exponents in the range 2.35 ≤ p ≤ 2.9. Acceptable models predict about 15% of the total mass in form of long‐lived stars and remnants, at the end of halo evolution, with a mean gas oxygen abundance which is substantially lower than the mean bulge and initial disk oxygen abundance. To avoid this discrepancy, either the existence of a still undetected, baryonic dark halo with about 15% of the total mass, or an equal amount of gas loss during bulge and disk formation, is necessary. The latter alternative implies a lower stellar mass limit close to 0.2 m_⊙, which is related to a power‐law IMF exponent close to 2.77. Acceptable models also imply a rapid halo formation, mainly during the first step, Δt = 0.5 Gyr, followed by a period (three steps) where small changes occur. Accordingly, statistical fluctuations are found to produce only minor effects on the evolution.     

The Chemical Evolution of the Galactic Bulge

In this paper we review the chemical evolution models for the Galactic bulge: in particular, we discuss the predictions of models as compared with the available abundance data and infer the mechanism as well as the time scale for the formation of the Galactic bulge. We show that good chemical evolution models reproducing the observed metallicity distribution of stars in the bulge predict that the [α/Fe] >0 over most of the metallicity range. This is a very important constraint indicating that the bulge of our Galaxy formed at the same time and even faster than the inner Galactic halo. We also discuss predictions for the evolution of light elements such as D and ⁷Li and conclude that the D astration should be maximum due to the high star formation rate required for the bulge whereas the evolution of the abundance of Li should be similar to that observed in the solar neighbourhood, but with an higher Li abundance in the interstellar medium at the present time. This revised version was published online in July 2006 with corrections to the Cover Date.