Errata for MNRAS 330, 821 (2002, astro-ph/0111084). A notation error in Eq. 7 and some typos in Table 3 corrected.

In Smoothed Particles Hydrodynamics (SPH) codes with a large number of particles, star formation as well as gas and metal restitution from dying stars can be treated statistically. This approach allows to include detailed chemical evolution and gas re-ejection with minor computational effort. Here we report on a new statistical algorithm for star formation and chemical evolution, especially conceived for SPH simulations with large numbers of particles, and for parallel SPH codes. For the sake of illustration, we present also two astrophysical simulations obtained with this algorithm, implemented into the Tree-SPH code by Lia & Carraro (2000). In the first one, we follow the formation of an individual disc-like galaxy, predict the final structure and metallicity evolution, and test resolution effects. In the second one we simulate the formation and evolution of a cluster of galaxies, to demonstrate the capabilities of the algorithm in investigating the chemo-dynamical evolution of galaxies and of the intergalactic medium in a cosmological context.

We present a new statistical algorithm of Chemical Evolution implemented into the Tree-SPH code developed by Lia & Carraro (2000 MNRAS 314, 145).

The aim of this paper is to cast light on the formation and evolution of elliptical galaxies by means of N-body/hydro-dynamical simulations that include star formation, feed-back and chemical evolution. Particular attention is paid to the case of dwarf spheroidals of the Local Group which, thanks to their proximity and modern ground-based and space instrumentation, can be resolved into single stars so that independent determinations of their age and star formation history can be derived. Dwarf galaxies are known to exhibit complicated histories of star formation ranging from a single very old episode to a series of bursts over most of the Hubble time. We start from virialized haloes of dark matter, and follow the infall of gas into the potential wells and the formation of stars. We find that in objects of the same total mass, different star formation histories are possible, if the collapse phase started at different initial densities. We predict the final structure of dwarf spheroidal galaxies, their kinematics, their large scale distribution of gas and stars, and their detailed histories of the star formation and metal enrichment. Using a population synthesis technique, star formation and metal enrichment rates are then adopted to generate the present color-magnitude diagrams of the stellar populations hosted by dwarf spheroidal galaxies. varying the cosmological parameters $H_0$ and $q_0$. compared with the observational ones for some dwarf spheroidals of the Local Group.

Using N-body/hydrodynamical simulations which include prescriptions for Star Formation, Feed-Back and Chemical Evolution, we explore the interaction between baryons and Dark Matter at galactic scale. The N-body simulations are performed using a Tree-SPH code that follows the evolution of individual DM halos inside which stars form from cooling gas, and evolve delivering in the interstellar medium mass, metals, and energy. We examine the formation and evolution of a giant and a dwarf elliptical galaxy of total mass 10^12 and 10^9 solar masses, respectively. Starting from an initial density profile like the universal Navarro profile in the inner region, baryons sink towards the center due to cooling energy losses. At the end of the collapse, the innermost part (1/20 of the halo size) of the galaxy is baryon-dominated, whereas the outer regions are DM dominated. The Star Formation proceeds at much faster speed in the giant galaxy where a spheroid of 8 \times 10^10 solar masses is formed in 2 Gyr, with respect to the dwarf galaxy where the spheroid of 2 \times 10^7 solar masses is formed in 4 Gyr. For the two objects the final distribution of stars is well fitted by a Hernquist profile with effective radii of r_e = 30 kpc, and 2.8 kpc, respectively. The dark-to-luminous tarnsition radius r_IBD occurs roughly at 1 r_e, as in real ellipticals. The DM halo density evolution is non-adiabatic and does not lead to a core radius.

By means of N-body/hydrodynamical simulations we model the evolution of a primordial 10^{8} solar masses density peak which ends up in an object closely resembling the present day globular cluster Omega Centauri. We succeed to reproduce the main features of the cluster, namely the structure, kinematics and metallicity distribution. We suggest that Omega Centauri might be a cosmological dwarf elliptical, formed at high redshift, evolved in isolation and self-enriched, and eventually fallen inside the potential well of the Milky Way, in agreement with the Searle-Zinn (1978) paradigm for galactic globular clusters formation. We finally suggest that Omega Centauri is probably surrounded by an extended Dark Matter (DM) halo, for which no observational evidence is at present available. We expect that signatures, if any, of the DM halo can be found in the kinematics of stars outside about 20 arcmin.

We describe a new implementation of a parallel Tree-SPH code with the aim to simulate Galaxy Formation and Evolution. The code has been parallelized using SHMEM, a Cray proprietary library to handle communications between the 256 processors of the Silicon Graphics T3E massively parallel supercomputer hosted by the Cineca Supercomputing Center (Bologna, Italy). The code combines the Smoothed Particle Hydrodynamics (SPH) method to solve hydro-dynamical equations with the popular Barnes and Hut (1986) tree-code to perform gravity calculation with a NlogN scaling, and it is based on the scalar Tree-SPH code developed by Carraro et al(1998)[MNRAS 297, 1021]. Parallelization is achieved distributing particles along processors according to a work-load criterion. Benchmarks, in terms of load-balance and scalability, of the code are analyzed and critically discussed against the adiabatic collapse of an isothermal gas sphere test using 20,000 particles on 8 processors. The code results balanced at more that 95% level. Increasing the number of processors, the load-balance slightly worsens. The deviation from perfect scalability at increasing number of processors is almost negligible up to 32 processors. Finally we present a simulation of the formation of an X-ray galaxy cluster in a flat cold dark matter cosmology, using 200,000 particles and 32 processors, and compare our results with Evrard (1988) P3M-SPH simulations. Additionaly we have incorporated radiative cooling, star formation, feed-back from SNae of type II and Ia, stellar winds and UV flux from massive stars, and an algorithm to follow the chemical enrichment of the inter-stellar medium. Simulations with some of these ingredients are also presented.

In this paper we present a critical discussion of the algorithms commonly used in N-body simulations of Galaxy Formation to deal with the energy equation governing heating and cooling, to model star formation and the star formation rate, and to account for energy feed-back from stars. First, we propose our technique for solving the energy equation in presence of heating and cooling, which includes some difference with respect to the standard semi-implicit technique. Second, we examine the current criteria for the onset of the star formation activity. We suggest a new approach, in which star formaiton is let depend on the total mass density - baryonic (gas and stars) and dark matter - of the system and on the metal-dependent cooling efficiency. Third, we check and discuss the separate effects of energy (and mass) feed-back from several sources - namely supernovae, stellar winds from massive stars, and UV flux from the same objects. All the simulations are performed in the framework of the formation and evolution of a disk galaxy. We show that the inclusion of these physical phenomena has a signigicant impact on the evolution of the galaxy model.

We describe a new implementation of a parallel Tree-SPH code with the aim to simulate Galaxy Formation and Evolution. The code has been parallelized using SHMEM, a Cray proprietary library to handle communications between the 256 processors of the Silicon Graphics T3E massively parallel supercomputer hosted by the Cineca Super-computing Center (Bologna, Italy). The code combines the Smoothed Particle Hydrodynamics (SPH) method to solve hydro-dynamical equations with the popular Barnes and Hut (1986) tree-code to perform gravity calculation with a $N \times logN$ scaling, and it is based on the scalar Tree-SPH code developed by Carraro et al (1998)[MNRAS 297, 1021]. Parallelization is achieved distributing particles along processors according to a work-load criterium. Benchmarks, in terms of load-balance and scalability, of the code are analised and critically discussed against the adiabatic collapse of an isothermal gas sphere test using $2 \times 10^{4}$ particles on 8 processors. The code results balanced at more than 95% level. Increasing the number of processors, the load balance sligthly worsens. The deviation from perfect scalability at increasing number of processors is negligible up to 64 processors. Additionally we have incorporated radiative cooling, star formation, feed-back and an algorithm to follow the chemical enrichment of the interstellar medium.

We present preliminary results on galactic Dark Matter (DM), halo structure, and galaxy evolution. We show how during the first Gyr of the evolution of a $10^{10} M_{\odot}$ dwarf elliptica feed-back from stars (SN\ae~ and stellar winds) leads to an extended constant density isothermal core with radius of 0.15 the virial radius $R_{200}$. We also present first results on galaxy merging as a possibile scenario to form ellipticals, studying in particular how the details of the merging evolution vary as a function of the mass ratio of the interacting galaxies.

In this report we describe a parallel implementation of a Tree-SPH code realized using the SHMEM libraries in the Cray T3E supercomputer at CINECA. We show the result of a 3D test to check the code performances against its scalar version. Finally we compare the load balancing and scalability of the code with PTreeSPH (Dav\'e et al 1997), the only other parallel Tree-SPH code present in the literature.

In this paper we report on PD-SPH the new tree-sph code developed in Padua. The main features of the code are described and the results of a new and independent series of 1-D and 3-D tests are shown. The paper is mainly dedicated to the presentation of the code and to the critical discussion of its performances. In particular great attention is devoted to the convergency analysis. The code is highly adaptive in space and time by means of individual smoothing lengths and individual time steps. At present it contains both dark and baryonic matter, this latter in form of gas and stars, cooling, thermal conduction, star formation, and feed-back from Type I and II supernovae, stellar winds, and ultraviolet flux from massive stars, and finally chemical enrichment. New cooling rates that depend on the metal abundance of the interstellar medium are employed, and the differences with respect to the standard ones are outlined. Finally, we show the simulation of the dynamical and chemical evolution of a disk-like galaxy with and without feed-back. The code is suitably designed to study in a global fashion the problem of formation and evolution of elliptical galaxies, and in particular to feed a spectro-photometric code from which the integrated spectra, magnitudes, and colors (together with their spatial gradients) can be derived.