Chemodynamical Modelling of Galaxy Formation and Evolution

We present our recently developed 3-dimensional chemodynamical code for galaxy evolution. This code follows the evolution of different galactic components like stars, dark matter and different components of the interstellar medium (ISM), i.e. a diffuse gaseous phase and the molecular clouds. Stars and dark matter are treated as collisionless N-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse gas and a sticky particle scheme for the molecular clouds. Additionally, the galactic components are coupled by several phase transitions like star formation, stellar death or condensation and evaporation processes within the ISM. As an example we show the dynamical and chemical evolution of a star forming dwarf galaxy with a total baryonic mass of 2 ċ 10⁹ M_⊙. After a moderate collapse phase the stars and the molecular clouds follow an exponential radial distribution, whereas the diffuse gas shows a central depression as a result of stellar feedback. The metallicities of the galactic components behave quite differently with respect to their temporal evolution as well as their radial distribution. Especially, the ISM is at no stage well mixed. This revised version was published online in September 2006 with corrections to the Cover Date.     

Evolution of a cloudy protogalaxy interacting with an early galactic wind and galaxy formation

We argue that a combined evidence from galactic and extragalactic studies suggests that a major star formation in giant galaxies is preceded by an evolutionary phase at which a strong galactic wind driven by the initial burst of star formation enriches the gaseous protogalaxy with metals and heats it up, so that the latter turns over from contraction to expansion. The result is the ejection of enriched material from the outer part of the protogalaxy into the intergalactic space, while the inner part, after a delay of about one to a few Gyr, finally contracts and cools down to form the galactic major stellar component (the hot model of galaxy formation). The paper presents a specific mechanism to produce a hot protogalaxy according to which an early galactic wind is imparting energy and momentum into a collapsing protogalaxy whose mass is contained mainly in clouds and only a small portion is in the intercloud gas that provides pressure confinement for the clouds. The model is then capable of accounting for the nearly equal mass and iron abundance in cluster giant galaxies and the intracluster gas provided the observationally plausible input parameters for giant galaxies and early galactic winds are adopted. It also predicts the formation of long-lived X-ray coronae with characteristics similar to those observed around giant ellipticals.The model specifies a characteristic length-scale that can be very naturally interpreted as a size for a stellar system to come; a very encouraging result is that it perfectly fits in with a typical size of giant ellipticals.     

Nucleosynthesis and Stellar Evolution

Two of the basic building blocks of galaxies are stars and the interstellar medium. The evolution of the abundance composition in the latter and especially the enrichment of heavy elements as a function of space and time reflects in turn the history of star formation and the lifetimes of the diverse contributing stellar objects. Therefore, the understanding of stellar evolution and its endpoints (mainly planetary nebulae, supernovae of type Ia and type II/Ib/Ic) is essential. Despite many efforts, a full and self-consistent understanding of supernovae (the main contributors to nucleosynthesis in galaxies) is not existing, yet. However, they leave fingerprints, seen either in spectra, lightcurves, radioactivities/decay gamma-rays or in galactic evolution. Here we want to address the composition of ejecta, their model uncertainties and relate them to constraints from abundance observations in galactic evolution. This revised version was published online in September 2006 with corrections to the Cover Date.     

The chemical evolution of the galactic disc from the kinematics and metallicities of proper motion stars

Members of the Galaxy components are identified according to stellar ages, metallicities and galactic orbits. The local thin disk is found to have a maximum age of 11 billion years and a small abundance scatter partially controlled by the radial gradient of abundances. Metal-rich and old metal-poor stars belong to inner galactic populations and SMRs represent the ultimate star generation in the bulge. The thick disk forms a smooth transition between the halo and thin disk.     

Mass loss of stars on the asymptotic giant branch

As low- and intermediate-mass stars reach the asymptotic giant branch (AGB), they have developed into intriguing and complex objects that are major players in the cosmic gas/dust cycle. At this stage, their appearance and evolution are strongly affected by a range of dynamical processes. Large-scale convective flows bring newly-formed chemical elements to the stellar surface and, together with pulsations, they trigger shock waves in the extended stellar atmosphere. There, massive outflows of gas and dust have their origin, which enrich the interstellar medium and, eventually, lead to a transformation of the cool luminous giants into white dwarfs. Dust grains forming in the upper atmospheric layers play a critical role in the wind acceleration process, by scattering and absorbing stellar photons and transferring their outward-directed momentum to the surrounding gas through collisions. Recent progress in high-angular-resolution instrumentation, from the visual to the radio regime, is leading to valuable new insights into the complex dynamical atmospheres of AGB stars and their wind-forming regions. Observations are revealing asymmetries and inhomogeneities in the photospheric and dust-forming layers which vary on time-scales of months, as well as more long-lived large-scale structures in the circumstellar envelopes. High-angular-resolution observations indicate at what distances from the stars dust condensation occurs, and they give information on the chemical composition and sizes of dust grains in the close vicinity of cool giants. These are essential constraints for building realistic models of wind acceleration and developing a predictive theory of mass loss for AGB stars, which is a crucial ingredient of stellar and galactic chemical evolution models. At present, it is still not fully possible to model all these phenomena from first principles, and to predict the mass-loss rate based on fundamental stellar parameters only. However, much progress has been made in recent years, which is described in this review. We complement this by discussing how observations of emission from circumstellar molecules and dust can be used to estimate the characteristics of the mass loss along the AGB, and in different environments. We also briefly touch upon the issue of binarity.     

Chemodynamical gas flow cycles and their influence on the chemical evolution of dwarf irregular galaxies

Here we investigate an exemplary chemodynamical evolutionary simulation of a dwarf irregular galaxy. By means of this model we demonstrate the existence of three gas mixing cycles: 1) An inner local cycle mixing the metals produced in stars locally, and 2) an outer galactic cycle on which hot gas is driven out of the galaxy by multiple supernovae type II and mixes on a short timescale with the available cold gas. 3) Only a small fraction of the metals leaves the galactic gravitational field and follows the global cycle with the intergalactic matter. The large-scale mixing results in a temporary depletion of supernova ejected metals. We will discuss this delayed recycling and its influence on the chemical evolution, especially on the nitrogen over oxygen ratio which is increased temporarily. The results presented here are also relevant for less sophisticated analytical approaches and chemical evolutionary models of galaxies which have to parameterize the metal loss through outflow. This revised version was published online in August 2006 with corrections to the Cover Date.     

A galactic model with a pulsating active nucleus

The chemical evolution of the Galaxy with a pulsating active nucleus is investigated. The surface densities of gas, stellar remnants, stars and chemical species such as helium and heavy elements inZ6 are calculated as functions of the position in the Galaxy and of the evolutional time of the Galaxy. According to this model, the entire luminosity of the galactic disk becomes almost constant at some 2×10⁹ yr after the galactic formation, but the nuclear bulge, whose dimensions gradually diminishes, becomes more and more luminous with time. On the other hand, the abundance depletion of helium and heavy elements appears in the inner region of the disk after some 6×10⁹ yr of the galactic formation. It also becomes clear that the activity for the nucleosynthesis in the nucleus is limited only in the early history of the Galaxy and has been reduced rapidly with time. Using this model, we can account for the observed phenomena such as the smooth dependence of the elemental abundance in the halo population on the distance from the galactic center, the high abundance of heavy elements in quasar spectra and etc.     

Constraining the star formation histories of spiral bulges

Stellar populations in spiral bulges are investigated using the Lick system of spectral indices. Long-slit spectroscopic observations of line strengths and kinematics made along the minor axes of four spiral bulges are reported. Comparisons are made between central line strengths in spiral bulges and those in other morphological types [elliptical, spheroidal (Sph) and S0]. The bulges investigated are found to have central line strengths comparable to those of single stellar populations of approximately solar abundance or above. Negative radial gradients are observed in line strengths, similar to those exhibited by elliptical galaxies. The bulge data are also consistent with correlations between Mg₂, Mg₂ gradient and central velocity dispersion observed in elliptical galaxies. In contrast to elliptical galaxies, central line strengths lie within the loci defining the range of 〈Fe〉 and Mg₂ achieved by Worthey's solar abundance ratio, single stellar populations (SSPs). The implication of solar abundance ratios indicates significant differences in the star formation histories of spiral bulges and elliptical galaxies. A 'single zone with infall' model of galactic chemical evolution, using Worthey's SSPs, is used to constrain the possible star formation histories of our sample. We show that the 〈Fe〉, Mg₂ and H β line strengths observed in these bulges cannot be reproduced using primordial collapse models of formation but can be reproduced by models with extended infall of gas and star formation (2–17 Gyr) in the region modelled. One galaxy (NGC 5689) shows a central population with a luminosity-weighted average age of ∼5 Gyr, supporting the idea of extended star formation. Kinematic substructure, possibly associated with a central spike in metallicity, is observed at the centre of the Sa galaxy NGC 3623.     

Evolution of the Galactic Halo and Disk

Correlations between stellar kinematics and chemical abundances are fossil evidence for evolutionary connections between Galactic structural components. Extensive stellar surveys show that the only tolerably clear distinction between galactic components appears in the distributions of specific angular momentum. Here the stellar metal-poor halo and the metal-rich bulge are indistinguishable from each other, as are the thick disk and the old disk. Each pair is very distinct from the other. This leads to an evolutionary model in which the metal-poor stellar halo evolves into the inner bulge, while the thick disk is a precursor to the thin disk. These evolutionary sequences are distinct. The galaxy is made of two discrete 'populations', one of low and one of high angular momentum. Some (minor?) complexity is added to this picture by the debris of late and continuing mergers, which will be especially important in the outer stellar halo.     

The star formation rate in disc galaxies: thresholds and dependence on gas amount

We reassess the applicability of the Toomre criterion in galactic discs and we study the local star formation law in 16 disc galaxies for which abundance gradients are published. The data we use consist of stellar light profiles, atomic and molecular gas (deduced from CO with a metallicity-dependent conversion factor), star formation rates (from Hα emissivities), metallicities, dispersion velocities and rotation curves. We show that the Toomre criterion applies successfully to the case of the Milky Way disc, but it has limited success with the data of our sample; depending on whether or not the stellar component is included in the stability analysis, we find average values for the threshold ratio of the gas surface density to the critical surface density in the range 0.5–0.7. We also test various star formation laws proposed in the literature, i.e. either the simple Schmidt law or modifications of it, that take into account dynamical factors. We find only small differences among them as far as the overall fit to our data is concerned; in particular, we find that all three star formation laws (with parameters derived from the fits to our data) match observations in the Milky Way disc particularly well. In all cases we find that the exponent n of our best-fitting star formation rate has slightly higher values than in other recent works and we suggest several reasons that may cause that discrepancy.     

A cosmological study of the star formation history in the solar neighbourhood

We use a cosmological galactic evolutionary approach to model the Milky Way. A detailed treatment of the mass aggregation and dynamical history of the growing dark halo is included, together with a self-consistent physical treatment for the star formation processes within the growing galactic disc. This allows us to calculate the temporal evolution of star and gas surface densities at all galactic radii, in particular, the star formation history (SFH) at the solar radius. A large range of cosmological mass aggregation histories (MAHs) is capable of producing a galaxy with the present-day properties of the Milky Way. The resulting SFHs for the solar neighbourhood bracket the available observational data for this feature, the most probable MAH yielding the optimal comparison with these observations. We also find that the rotation curve for our Galaxy implies the presence of a constant density core in its dark-matter halo.     

Warped galactic disks and optical projection effects

Direct and indirect evidence for global warping in galactic disks is reviewed. The existing data suggest that inversion symmetric warping is common in spiral galaxies, that it can occur throughout the inner and outer disk regions, and that it can involve stellar, dust, and gas components. The amplitude of the warping tends to be proportional to the radial distance from the galactic center and varies from instances wherein the disk is essentially flat to cases wherein thez-axis displacement is at least 30% of the galactocentric distance. Projection effects that are inherent features of warped optical disks are discussed.     

Testing model predictions of the cold dark matter cosmology for the sizes, colours, morphologies and luminosities of galaxies with the SDSS

The huge size and uniformity of the Sloan Digital Sky Survey (SDSS) make possible an exacting test of current models of galaxy formation. We compare the predictions of the galform semi-analytical galaxy formation model for the luminosities, morphologies, colours and scalelengths of local galaxies. galform models the luminosity and size of the disc and bulge components of a galaxy, and so we can compute quantities which can be compared directly with SDSS observations, such as the Petrosian magnitude and the Sérsic index. We test the predictions of two published models set in the cold dark matter cosmology: the Baugh et al. model, which assumes a top-heavy initial mass function (IMF) in starbursts and superwind feedback, and the Bower et al. model, which uses active galactic nucleus feedback and a standard IMF. The Bower et al. model better reproduces the overall shape of the luminosity function, the morphology–luminosity relation and the colour bimodality observed in the SDSS data, but gives a poor match to the size–luminosity relation. The Baugh et al. model successfully predicts the size–luminosity relation for late-type galaxies. Both models fail to reproduce the sizes of bright early-type galaxies. These problems highlight the need to understand better both the role of feedback processes in determining galaxy sizes, in particular the treatment of the angular momentum of gas reheated by supernovae, and the sizes of the stellar spheroids formed by galaxy mergers and disc instabilities.     

Cosmic Rays and Galactic Winds

Repeated SN-explosion provide large amounts of thermal energy as well as energetic particles through a 1. order Fermi-process. Both effects together with the generation of Alfvén-waves are considered to drive a large scale outflow from a galaxy. These so-called galactic windstransport stellar material enriched by heavy elements into the intergalactic space explaining also the large amount of metals found inthe intergalactic gas. The present contribution is focused on time-dependenteffects which originate from galactic winds driven by a star burst activity. Shock waves travelling through the galactic wind and radiative cooling within the expanding plasma lead to complex flow structures. Depending e.g. on theSFR of the galaxy galactic winds can remove almost all ISM into the galactic halo and therefore cease a subsequent star formation.     

Photoionization modelling of planetary nebulae – II. Galactic bulge nebulae, a comparison with literature results

We have constructed photoionization models of five galactic bulge planetary nebulae using our automatic method, which enables a fully self-consistent determination of the physical parameters of a planetary nebula. The models are constrained using the spectrum, the IRAS and radio fluxes and the angular diameter of the nebula. We also conducted a literature search for physical parameters determined with classical methods for these nebulae. Comparison of the distance-independent physical parameters with published data shows that the stellar temperatures generally are in good agreement and can be considered reliable. The literature data for the electron temperature, electron density and also for the abundances show a large spread, indicating that the use of line diagnostics is not reliable and that the accuracy of these methods needs to be improved. Comparison of the various abundance determinations indicates that the uncertainty in the electron temperature is the main source of uncertainty in the abundance determination. The stellar magnitudes predicted by the photoionization models are in good agreement with observed values.     

银盘的径向金融丰度梯度

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

Galactic mergers and gravitationally unbound populations

Motivated by the observations on the intracluster light and intergalactic stellar populations, N -body simulations are used to model the galactic merging events as a goal to investigate the production and distribution of gravitational unbound populations (GUPs). Both the parabolic and hyperbolic mergers are considered, and each category includes six models with different relative orientations between two galaxies. Our results show that there are more (about a factor of 2) GUPs after a hyperbolic merging event than after a parabolic one. In general, depending on the relative orientation and also on the relative velocity of the two galaxies in a merging pair, a head-on collision of a galaxy pair would only make a tiny fraction (less than 1 per cent) of the initial stellar mass luminous GUP, but a considerable fraction (8–14 per cent) of the dark matter becomes dark GUP.     

Gas flow in a two component galactic disk

Numerical simulations of two-component (stars + gas) self-gravitating galactic disks show that the interstellar gas can significantly affect the dynamical evolution of the disk even if its mass fraction (relative to the total galaxy mass) is as low as several percent. Aided by efficient energy dissipation, the gas becomes gravitationally unstable onlocal scale and forms massive clumps. Gravitational scattering of stars by these clumps leads to suppression of bar instability usually seen in heavy stellar disks. In this case, gas inflow towards the galactic center is driven by dynamical friction which gas clumps suffer instead of bar forcing.     

Large-scale flow of interstellar gas in galactic spiral waves: Effects of thermal balance and self-gravitation

Using the recent observational data on atomic and molecular hydrogen in the Galaxy, we analyse the dynamics of the interstellar gas in a spiral density wave. Within the framework of Marochniket al.'s (1972) model of the galactic spiral structure, the gas flow is obtained, with self-gravitation and thermal processes taken into account.It is shown that: (1) the self-gravitation of gas does not practically affect the galactic shock if the dominant contribution into the gas density comes from atomic hydrogen; (2) the effects of self-gravitation could be essential for both the gas flow and the stellar spiral wave only if the density contribution of H₂ exceeded several times that ofHi, with molecular hydrogen as a continuous medium having the isothermal equation of state; (3) however, regardless of the estimates of H₂ abundance in the Galaxy, its reaction to the density wave is weak, since it forms a collisionless system not dragged by the intercloud gas.It has been found that, if we allow for thermal processes in the interstellar medium, new types of gas flow can develop alongside with a previously-known continuous flow and galactic shock. They are: (1) galactic shock with the phase transition leading to the formation of dense cold clouds; (2) a three-phase flow where hot cavities and dense cold clouds coexist with an initial, moderately dense and cold phase; (3) an accretion wave which is a specific type of nonlinear wave with an amplitude of 11/2 orders of magnitude larger than that of the isothermal galactic shock appearing under the same conditions, but without heating and cooling.     

Gas-driven evolution of stellar orbits in barred galaxies

We carry out a detailed orbit analysis of gravitational potentials selected at different times from an evolving self-consistent model galaxy consisting of a two-component disc (stars+gas) and a live halo. The results are compared with a pure stellar model, subject to nearly identical initial conditions, which are chosen so as to make the models develop a large-scale stellar bar. The bars are also subject to hose-pipe (buckling) instability which modifies the vertical structure of the disc. The diverging morphological evolution of both models is explained in terms of gas radial inflow, the resulting change in the gravitational potential at smaller radii, and the subsequent modification of the main families of orbits, both in and out of the disc plane.   We find that dynamical instabilities become milder in the presence of the gas component, and that the stability of planar and 3D stellar orbits is strongly affected by the related changes in the potential — both are destabilized, with the gas accumulation at the centre. This is reflected in the overall lower amplitude of the bar mode and in the substantial weakening of the bar, which appears to be a gradual process. The vertical buckling of the bar is much less pronounced and the characteristic peanut shape of the galactic bulge almost disappears when there is a substantial gas inflow towards the centre. Milder instability results in a smaller bulge, the basic parameters of which are in agreement with observations. We also find that the overall evolution in the model with a gas component is accelerated because of the larger central mass concentration and the resulting decrease in the characteristic dynamical time.