Context. This is the second work dedicated to the observed parallelism between galaxy clusters and early-type galaxies. The focus is on the distribution of these systems in the scaling relations (SRs) observed when effective radii, effective surface brightness, total luminosities and velocity dispersions are mutually correlated. Aims. Using the data of the Illustris simulation we try to speculate on the origin of the observed SRs. Methods. We compare the observational SRs extracted from the database of the WIde-field Nearby Galaxy-cluster Survey (WINGS) with the relevant parameters coming from the Illustris simulations. Then we use the simulated data at different redshift to infer the evolution of the SRs. Results. The comparison demonstrate that galaxy clusters (GCs) at z~0 follow the same log(L)-log(sigma) relation of early-type galaxies (ETGs) and that both in the log(Ie)-log(Re) and log(Re)-log(M*) planes the distribution of GCs is along the sequence defined by the brightest and massive early-type galaxies (BCGs). The Illustris simulation reproduces the tails of the massive galaxies visible both in the log(Ie)-log(Re) and log(Re)-log(M*) planes, but fail to give the correct estimate of the effective radii of the dwarf galaxies that appear too large and those of GCs that are too small. The evolution of the SRs up to z=4 permits to reveal the complex evolutionary paths of galaxies in the SRs and indicate that the line marking the Zone of Exclusion (ZoE), visible both in the log(Ie)-log(Re) and log(Re)-log(M*) planes, is the trend followed by virialized and passively evolving systems. Conclusions. We speculate that the observed SRs originate from the intersection of the virial theorem and a relation L=L_0 x sigma^beta where the luminosities depend on the star formation history.

Context. This is the third study of a series dedicated to the observed parallelism of properties between Galaxy Clusters and Groups(GCGs) and early-type galaxies (ETGs). Aims. Here we investigate the physical origin of the Mass-Radius Relation (MRR). Methods. Having collected literature data on masses and radii for objects going from Globular Clusters (GCs) to ETGs and GCGs, we set up the MR-plane and compare the observed distribution with the MRR predicted by theoretical models both for the monolithic and hierarchical scenarios. Results. We argue that the distributions of stellar systems in the MR-plane is due to complementary mechanisms: (i) on one hand, as shown in paper II, the relation of the virial equilibrium does intersect with a relation that provides the total luminosity as a function of the star formation history; (ii) on the other hand, the locus predicted for the collapse of systems should be convolved with the statistical expectation for the maximum mass of the halos at each comsic epoch. This second aspect provides a natural boundary limit explaining either the curved distribution observed in the MR-plane and the existence of a zone of avoidance. Conclusions. The distribution of stellar systems in the MR-plane is the result of two combined evolution, that of the stellar component and that of the halo component.

Convection is one of the fundamental energy transport processes in physics and astrophysics, and its description is central to allstellar models. In the context of stellar astrophysics, the mixing length theory is the most successful approximation to handle theconvection zones inside the stars because of its simplicity and rapidity. The price to pay is the mixing length parameter that isintroduced to derive the velocity of convective elements, the temperature gradients in the convective regions and finally the energy flux carried by convection. The mixing length is a free parameter that needs to be calibrated on observational data. Pasetto et al. (2014) have proposed a new theory that determines all the properties of convective regions and the convective transport of energy with no need for a free parameter. In this study, we aim to discuss the merits of this new approach and the limits of its applicability in comparison with the mixing length theory. We present an analytical and numerical investigation of the main physical assumptions made by Pasetto et al. (2014) and compare them with the counterparts of the mixing length theory. We also present here the homogeneous isotropic limit of the Pasetto et al. (2014) theory and discuss some numerical examples to address and clarify misconceptions often associated with the new formalism.

We have analyzed the parallelism between the properties of galaxy clusters and early-type galaxies (ETGs) by looking at the similarity between their light profiles. We find that the equivalent luminosity profiles of all these systems in the \vfilt\ band, once normalized to the effective radius \re\ and shifted in surface brightness, can be fitted by the S\'ersic's law \Sers\ and superposed with a small scatter ($\le0.3$ mag). By grouping objects in different classes of luminosity, the average profile of each class slightly deviates from the other only in the inner and outer regions (outside $0.1\leq r/R_e\leq 3$), but the range of values of $n$ remains ample for the members of each class, indicating that objects with similar luminosity have quite different shapes. The "Illustris" simulation reproduces quite well the luminosity profiles of ETGs, with the exception of in the inner and outer regions where feedback from supernovae and active galactic nuclei, wet and dry mergers, are at work. The total mass and luminosity of galaxy clusters as well as their light profiles are not well reproduced. By exploiting simulations we have followed the variation of the effective half-light and half-mass radius of ETGs up to $z=0.8$, noting that progenitors are not necessarily smaller in size than current objects. We have also analyzed the projected dark+baryonic and dark-only mass profiles discovering that after a normalization to the half-mass radius, they can be well superposed and fitted by the S\'ersic's law.

We investigate the color-magnitude diagram (CMD) of cluster galaxies in the hierarchical $\Lambda$-CDM cosmological scenario using both single stellar populations and simple galaxy models. First, we analyze the effect of bursts and mergers and companion chemical pollution and rejuvenation of the stellar content on the integrated light emitted by galaxies. The dispersion of the galaxy magnitudes and colors on the $M_V-(B-V)$ plane is mainly due to mixing of ages and metallicities of the stellar populations, with mergers weighting more than bursts of similar mass fractions. The analysis is made using the Monte-Carlo technique applied to ideal model galaxies reduced to single stellar populations with galaxy-size mass to evaluate mass, age and metallicity of each object. We show that separately determining the contributions by bursts and mergers leads to a better understanding of observed properties of CMD of cluster galaxies. Then we repeat the analysis using suitable chemo-photometric models of galaxies whose mass is derived from the cosmological predictions of the galaxy content of typical clusters. Using the halo mass function and the Monte-Carlo technique, we derive the formation redshift of each galaxy and its photometric history. These are used to simulate the CMD of the cluster galaxies. The main conclusion is that most massive galaxies have acquired the red color they show today in very early epochs and remained the same ever since. The simulations nicely reproduce the Red Sequence, the Green Valley and the Blue Cloud, the three main regions of the CMD in which galaxies crowd.

We aim to investigate the connections existing between the density profiles of the stellar populations used to define a gravitationally bound stellar system and their star formation history: we do this by developing a general framework accounting for both classical stellar population theory and classical stellar dynamics. We extend the work of Pasetto et al. (2012) on a single composite-stellar population (CSP) to multiple CSPs, including also a phase-space description of the CSP concept. In this framework, we use the concept of distribution function to define the CSP in terms of mass, metallicity, and phase-space in a suitable space of existence $\mathbb{E}$ of the CSP. We introduce the concept of foliation of $\mathbb{E}$ to describe formally any CSP as sum of disjointed Simple Stellar Populations (SSP), with the aim to offer a more general formal setting to cast the equations of stellar populations theory and stellar dynamics theory. In doing so, we allow the CSP to be object of dissipation processes thus developing its dynamics in a general non-Hamiltonian framework. Furthermore, we investigate the necessary and sufficient condition to realize a multiple CSP consistent with its mass-metallicity and phase-space distribution function over its temporal evolution, for a collisionless CSP. Finally, analytical and numerical examples show the potential of the result obtained.

We present a new Galaxy population synthesis Model (GalMod). GalMod is a star-count model featuring an asymmetric bar/bulge as well as spiral arms and related extinction. The model, initially introduced in Pasetto et al. (2016b), has been here completed with a central bar, a new bulge description, new disk vertical profiles and several new bolometric corrections. The model can generate synthetic mock catalogs of visible portions of the Milky Way (MW), external galaxies like M31, or N-body simulation initial conditions. At any given time, e.g., a chosen age of the Galaxy, the model contains a sum of discrete stellar populations, namely bulge/bar, disk, halo. These populations are in turn the sum of different components: the disk is the sum of spiral arms, thin disks, a thick disk, and various gas components, while the halo is the sum of a stellar component, a hot coronal gas, and a dark matter component. The Galactic potential is computed from these population density profiles and used to generate detailed kinematics by considering up to the first four moments of the collisionless Boltzmann equation. The same density profiles are then used to define the observed color-magnitude diagrams in a user-defined field of view from an arbitrary solar location. Several photometric systems have been included and made available on-line and no limits on the size of the field of view are imposed thus allowing full sky simulations, too. Finally, we model the extinction adopting a dust model with advanced ray-tracing solutions. The model's web page (and tutorial) can be accessed at www.GalMod.org.

We investigate the evolution of the cosmic star formation rate density (SFRD) from redshift z=20 to z=0 and compare it with the observational one by Madau and Dickinson derived from recent compilations of UV and IR data. The theoretical SFRD(z) and its evolution are obtained using a simple model which folds together the star formation histories of prototype galaxies designed to represent real objects of different morphological type along the Hubble sequence and the hierarchical growing of structures under the action of gravity from small perturbations to large scale objects in \Lambda-CDM cosmogony, i.e. the number density of dark matter halos N(M,z). Although the overall model is very simple and easy to set up, it provides results that well mimic those obtained from large scale N-body simulations of great complexity. The simplicity of our approach allows us to test different assumptions for the star formation law in galaxies, the effects of energy feedback from stars to interstellar gas and the efficiency of galactic winds, and also the effect of N(M,z). The result of our analysis is that in the framework of the hierarchical assembly of galaxies the so-called time-delayed star formation under plain assumptions mainly for the energy feedback and galactic winds can reproduce the observational SFRD(z).

The aim of this work is to show that the origin of the Fundamental Plane (FP) relation for early-type galaxies (ETGs) can be traced back to the existence of a fine-tuning between the average star formation rate $<SFR>$ of galaxies and their structural and dynamical characteristics. To get such result it is necessary to imagine the existence of two distinct "virtual planes" for each galaxy in the $\log(R_e)-\log(I_e)-\log(\sigma)$ space. The first one (named Virial Plane VP) represents the total galaxy mass using the scalar Virial Theorem and the mass-to-light ratio $M/L$, while the second plane comes from an expression of the total galaxy luminosity as a function of the mean star formation rate $<SFR>$ and the velocity dispersion $\sigma$, through a relation $L=L'_0 \sigma^{-2}$ (named here pseudo-Faber-Jackson (PFJ)) which is a mathematical convenient way for expressing the independency of light from the virial equilibrium. Its validity can be connected to the mutual correlation $L\sim\sigma\sqrt{<SFR>}$ observed for all ETGs. A posteriori it is possible to see that this approach permits to explain the observed properties of the FP (tilt and scatter) and the Zone of Exclusions (ZOE) visible in the FP projections. Furthermore, the link between the properties of the FP and the SFR of galaxies provides a new idea of the star formation, as a phenomenon driven by the initial conditions of proto-galaxies and regulated across the whole cosmic history by the variation of the main galaxy parameters (mass, luminosity, structural shape and velocity dispersion).

We present a new theoretical population synthesis model (the Galaxy Model) to examine and deal with large amounts of data from surveys of the Milky Way and to decipher the present and past structure and history of our own Galaxy. We assume the Galaxy to consist of a superposition of many composite stellar populations belonging to the thin and thick disks, the stellar halo and the bulge, and to be surrounded by a single dark matter halo component. A global model for the Milky Way's gravitational potential is built up self-consistently with the density profiles from the Poisson equation. In turn, these density profiles are used to generate synthetic probability distribution functions (PDFs) for the distribution of stars in colour-magnitude diagrams (CMDs). Finally, the gravitational potential is used to constrain the stellar kinematics by means of the moment method on a (perturbed)-distribution function. Spiral arms perturb the axisymmetric disk distribution functions in the linear response framework of density-wave theory where we present an analytical formula of the so-called `reduction factor' using Hypergeometric functions. Finally, we consider an analytical non-axisymmetric model of extinction and an algorithm based on the concept of probability distribution function to handle colour magnitude diagrams with a large number of stars. A genetic algorithm is presented to investigate both the photometric and kinematic parameter space. This galaxy model represents the natural framework to reconstruct the structure of the Milky Way from the heterogeneous data set of surveys such as Gaia-ESO, SEGUE, APOGEE2, RAVE and the Gaia mission.

We present here the first stellar models on the Hertzsprung-Russell diagram (HRD), in which convection is treated according to the novel scale-free convection theory (SFC theory) by Pasetto et al. (2014). The aim is to compare the results of the new theory with those from the classical, calibrated mixing-length (ML) theory to examine differences and similarities. We integrate the equations describing the structure of the atmosphere from the stellar surface down to a few percent of the stellar mass using both ML theory and SFC theory. The key temperature over pressure gradients, the energy fluxes, and the extension of the convective zones are compared in both theories. The analysis is first made for the Sun and then extended to other stars of different mass and evolutionary stage. The results are adequate: the SFC theory yields convective zones, temperature gradients of the ambient and of the convective element, and energy fluxes that are very similar to those derived from the "calibrated" MT theory for main sequence stars. We conclude that the old scale dependent ML theory can now be replaced with a self-consistent scale-free theory able to predict correct results, one which is simpler and more physically grounded than the ML theory. Fundamentally, the SFC theory offers a deeper insight of the underlying physics than numerical simulations.

In this paper we explore the possibility that isolated CO-WDs with mass smaller than the Chandrasekhar limit may undergo nuclear runaway and SNa explosion, triggered by the energy produced by under-barrier pycno-nuclear reactions between carbon and light elements. Such reactions would be due to left over impurities of the light elements, which would remain inactive until the WDs transit from the liquid to the solid state. We devise a simple formulation for the coulombian potential and the local density in a ionic lattice affected by impurities and introduce it in the known rates of pycno-nuclear reactions for multi-component plasmas. Our semi-analytical results indicate that the energy generated by these pycno-nuclear reactions exceeds the WD luminosity and provides enough energy to elementary cells of matter to balance the energy cost for C-ignition at much younger ages than the age of the Universe, even for WDs with masses as low as $\simeq 0.85\, M_\odot$. A thermonuclear runaway may thus be triggered in isolated WDs. The explosion would occur from few hundred thousand to a few million years after the WD formation in the mass interval $0.85 - 1.2\, M_\odot$.

Improving upon the standard evolutionary population synthesis (EPS) technique, we present spectrophotometric models of galaxies whose morphology goes from spherical structures to discs, properly accounting for the effect of dust in the interstellar medium (ISM). These models enclose three main physical components: the diffuse ISM composed by gas and dust, the complexes of molecular clouds (MCs) where active star formation occurs and the stars of any age and chemical composition. These models are based on robust evolutionary chemical models that provide the total amount of gas and stars present at any age and that are adjusted in order to match the gross properties of galaxies of different morphological type. We have employed the results for the properties of the ISM presented in Piovan, Tantalo & Chiosi (2006a) and the single stellar populations calculated by Cassar\`a et al. (2013) to derive the spectral energy distributions (SEDs) of galaxies going from pure bulge to discs passing through a number of composite systems with different combinations of the two components. The first part of the paper is devoted to recall the technical details of the method and the basic relations driving the interaction between the physical components of the galaxy. Then, the main parameters are examined and their effects on the spectral energy distribution of three prototype galaxies are highlighted. We conclude analyzing the capability of our galaxy models in reproducing the SEDs of real galaxies in the Local Universe and as a function of redshift.

Stellar convection is customarily described by Mixing-Length Theory, which makes use of the mixing-length scale to express the convective flux, velocity, and temperature gradients of the convective elements and stellar medium. The mixing-length scale is taken to be proportional to the local pressure scale height, and the proportionality factor (the mixing-length parameter) must be determined by comparing the stellar models to some calibrator, usually the Sun. No strong arguments exist to suggest that the mixing-length parameter is the same in all stars and at all evolutionary phases. The aim of this study is to present a new theory of stellar convection that does not require the mixing length parameter. We present a self-consistent analytical formulation of stellar convection that determines the properties of stellar convection as a function of the physical behaviour of the convective elements themselves and of the surrounding medium. This new theory is formulated starting from a conventional solution of the Navier-Stokes/Euler equations, i.e. the Bernoulli equation for a perfect fluid, but expressed in a non-inertial reference frame co-moving with the convective elements. In our formalism the motion of stellar convective cells inside convectively-unstable layers is fully determined by a new system of equations for convection in a non-local and time-dependent formalism. We obtain an analytical, non-local, time-dependent sub-sonic solution for the convective energy transport that does not depend on any free parameter. The theory is suitable for the outer convective zones of solar type stars and stars of all mass on the main sequence band. The predictions of the new theory are compared with those from the standard mixing-length paradigm for the most accurate calibrator, the Sun, with very satisfactory results.

We develop a simple analytical criterion to investigate the role of the environment on the onset of star formation. We will consider the main external agents that influence the star formation (i.e. ram pressure, tidal interaction, Rayleigh-Taylor and Kelvin-Helmholtz instabilities) in a spherical galaxy moving through an external environment. The theoretical framework developed here has direct applications to the cases of dwarf galaxies in galaxy clusters and dwarf galaxies orbiting our Milky Way system, as well as any primordial gas-rich cluster of stars orbiting within its host galaxy. We develop an analytic formalism to solve the fluid dynamics equations in a non-inertial reference frame mapped with spherical coordinates. The two-fluids instability at the interface between a stellar system and its surrounding hotter and less dense environment is related to the star formation processes through a set of differential equations. The solution presented here is quite general, allowing us to investigate most kinds of orbits allowed in a gravitationally bound system of stars in interaction with a major massive companion. We present an analytical criterion to elucidate the dependence of star formation in a spherical stellar system (as a dwarf galaxy or a globular cluster) on its surrounding environment useful in theoretical interpretations of numerical results as well as observational applications. We show how spherical coordinates naturally enlighten the interpretation of the two-fluids instability in a geometry that directly applies to astrophysical case. This criterion predicts the threshold value for the onset of star formation in a mass vs. size space for any orbit of interest. Moreover, we show for the first time the theoretical dependencies of the different instability phenomena acting on a system in a fully analytical way.

We investigate the influence of Asymptotic Giant Branch stars on integrated colours of star clusters of ages between ~100 Myr and a few gigayears, and composition typical for the Magellanic Clouds. We use state-of-the-art stellar evolution models that cover the full thermal pulse phase, and take into account the influence of dusty envelopes on the emerging spectra. We present an alternative approach to the usual isochrone method, and compute integrated fluxes and colours using a Monte Carlo technique that enables us to take into account statistical fluctuations due to the typical small number of cluster stars. We demonstrate how the statistical variations in the number of Asymptotic Giant Branch stars and the temperature and luminosity variations during thermal pulses fundamentally limit the accuracy of the comparison (and calibration, for population synthesis models that require a calibration of the Asymptotic Giant Branch contribution to the total luminosity) with star cluster integrated photometries. When compared to observed integrated colours of individual and stacked clusters in the Magellanic Clouds, our predictions match well most of the observations, when statistical fluctuations are taken into account, although there are discrepancies in narrow age ranges with some (but not all) set of observations.

We have employed state-of-the-art evolutionary models of low and intermediate-mass AGB stars, and included the effect of circumstellar dust shells on the spectral energy distribution (SED) of AGB stars, to revise the Padua library of isochrones (Bertelli et al. 1994). The major revision involves the thermally pulsing AGB phase, that is now taken from fully evolutionary calculations by Weiss & Ferguson (2009). Two libraries of about 600 AGB dust-enshrouded SEDs each have also been calculated, one for oxygen-rich M-stars and one for carbon-rich C-stars. Each library accounts for different values of input parameters like the optical depth {\tau}, dust composition, and temperature of the inner boundary of the dust shell. These libraries of dusty AGB spectra have been implemented into a large composite library of theoretical stellar spectra, to cover all regions of the Hertzsprung-Russell Diagram (HRD) crossed by the isochrones. With the aid of the above isochrones and libraries of stellar SEDs, we have calculated the spectro-photometric properties (SEDs, magnitudes, and colours) of single-generation stellar populations (SSPs) for six metallicities, more than fifty ages (from 3 Myr to 15 Gyr), and nine choices of the Initial Mass Function. The new isochrones and SSPs have been compared to the colour-magnitude diagrams (CMDs) of field populations in the LMC and SMC, with particular emphasis on AGB stars, and the integrated colours of star clusters in the same galaxies, using data from the SAGE (Surveying the Agents of Galaxy Evolution) catalogues. We have also examined the integrated colours of a small sample of star clusters located in the outskirts of M31. The agreement between theory and observations is generally good. In particular, the new SSPs reproduce the red tails of the AGB star distribution in the CMDs of field stars in the Magellanic Clouds.

Radial velocity surveys such as the Radial Velocity Experiment (RAVE) provide us with measurements of hundreds of thousands of nearby stars most of which belong to the Galactic thin, thick disk or halo. Ideally, to study the Galactic disks (both thin and thick) one should make use of the multi-dimensional phase-space and the whole pattern of chemical abundances of their stellar populations. In this paper, with the aid of the RAVE Survey, we study the thin and thick disks of the Milky Way, focusing on the latter. We present a technique to disentangle the stellar content of the two disks based on the kinematics and other stellar parameters such as the surface gravity of the stars. Using the Padova Galaxy Model, we checked the ability of our method to correctly isolate the thick disk component from the Galaxy mixture of stellar populations. We introduce selection criteria in order to clean the observed radial velocities from the Galactic differential rotation and to take into account the partial sky coverage of RAVE. We developed a numerical technique to statistically disentangle thin and thick disks from their mixture. We deduce the components of the solar motion relative to the Local Standard of Rest (LSR) in the radial and vertical direction, the rotational lag of the thick disk component relative to the LSR, and the square root of the absolute value of the velocity dispersion tensor for the thick disk alone. The analysis of the thin disk is presented in another paper. We find good agreement with previous independent parameter determinations. In our analysis we used photometrically determined distances. In the Appendix we show that similar values can be found for the thick disk alone as derived in the main sections of our paper even without the knowledge of photometric distances.

Aims. We study the Milky Way thin disk with the Radial Velocity Experiment (RAVE) survey. We consider the thin and thick disks as different Galactic components and present a technique to statistically disentangle the two populations. Then we focus our attention on the thin disk component. Methods. We disentangle the thin disk component from amixture of the thin and thick disks using a data set providing radial velocities, proper motions, and photometrically determined distances. Results. We present the trend of the velocity dispersions in the thin disk component of the Milky Way (MW) in the radial direction above and below the Galactic plane using data from the RAdial Velocity Experiment (RAVE). The selected sample is a limited subsample from the entire RAVE catalogue, roughly mapping up to 500 pc above and below the Galactic plane, a few degrees in azimuthal direction and covering a radial extension of 2.0 kpc around the solar position. The solar motion relative to the local standard of rest is also re-determined with the isolated thin disk component. Major results are the trend of the velocity mean and dispersion in the radial and vertical direction. In addition the azimuthal components of the solar motion relative to the local standard of rest and the velocity dispersion are discussed.

Early-type galaxies obey a narrow relation traced by their stellar content between the mass and size (Mass- Radius relation). The wealth of recently acquired observational data essentially confirms the classical relations found by Burstein, Bender, Faber, and Nolthenius, i.e. log(R_1/2) \propto log(Ms)\simeq 0.54 for high mass galaxies and log(R_1/2) \propto log(Ms) \simeq 0.3 for dwarf systems (shallower slope), where R_1/2 and Ms are the half-light radius and total mass in stars, respectively. Why do galaxies follow these characteristic trends? What can they tell us about the process of galaxy formation? We investigate the mechanisms which concur to shape the Mass-Radius relation, in order to cast light on the physical origin of its slope, its tightness, and its zero point. We perform a theoretical analysis, and couple it with the results of numerical hydrodynamical (NB-TSPH) simulations of galaxy formation, and with a simulation of the Mass-Radius plane itself. We propose a novel interpretation of the Mass-Radius relation, which we claim to be the result of two complementary mechanisms: on one hand, the result of local physical processes, which fixes the ratio between masses and radii of individual objects; on the other hand, the action of cosmological global, statistical principles, which shape the distribution of objects in the plane. We reproduce the Mass-Radius relation with a simple numerical technique based on this view.

We investigate the influence of the initial proto-galaxies over-densities and masses on their evolution, to understand whether the internal properties of the proto-galactic haloes are sufficient to account for the varied properties of the galactic populations. By means of fully hydrodynamical N-body simulations performed with the code EvoL we produce twelve self-similar models of early-type galaxies of different initial masses and over-densities, following their evolution from z \geq 20 down to z \leq 1. The simulations include radiative cooling, star formation, stellar energy feedback, a reionizing photoheating background, and chemical enrichment of the ISM. We find a strong correlation between the initial properties of the proto-haloes and their star formation histories. Massive (10^13M\odot) haloes experience a single, intense burst of star formation (with rates \geq 10^3M\odot/yr) at early epochs, consistently with observations, with a less pronounced dependence on the initial over-density; intermediate mass (10^11M\odot) haloes histories strongly depend on their initial over-density, whereas small (10^9M\odot) haloes always have fragmented histories, resulting in multiple stellar populations, due to the "galactic breathing" phenomenon. The galaxy models have morphological, structural and photometric properties comparable to real galaxies, often closely matching the observed data; even though some disagreement is still there, likely a consequence of some numerical choices. We conclude that internal properties are essentially sufficient to explain many of the observed features of early type galaxies, particularly the complicated and different star formation histories shown by haloes of very different mass. In this picture, nature seems to play the dominant role, whereas nurture has a secondary importance.

Aims. We present here a new theoretical approach to population synthesis. The aim is to predict colour magnitude diagrams (CMDs) for huge numbers of stars. With this method we generate synthetic CMDs for N-body simulations of galaxies. Sophisticated hydrodynamic N-body models of galaxies require equal quality simulations of the photometric properties of their stellar content. The only prerequisite for the method to work is very little information on the star formation and chemical enrichment histories, i.e. the age and metallicity of all star-particles as a function of time. The method takes into account the gap between the mass of real stars and that of the star-particles in N-body simulations, which best correspond to the mass of star clusters with different age and metallicity, i.e. a manifold of single stellar sopulations (SSP). Methods. The theory extends the concept of SSP to include the phase-space (position and velocity) of each star. Furthermore, it accelerates the building up of simulated CMD by using a database of theoretical SSPs that extends to all ages and metallicities of interest. Finally, it uses the concept of distribution functions to build up the CMD. The technique is independent of the mass resolution and the way the N-body simulation has been calculated. This allows us to generate CMDs for simulated stellar systems of any kind: from open clusters to globular clusters, dwarf galaxies, or spiral and elliptical galaxies. Results. The new theory is applied to an N-body simulation of a disc galaxy to test its performance and highlight its flexibility.

Dissipative phenomena occurring during the orbital evolution of a dwarf satellite galaxy around a host galaxy may leave signatures in the star formation activity and signatures in the colour magnitude diagram of the galaxy stellar content. Our goal is to reach a simple and qualitative description of the these complicated phenomena. We develop an analytical and numerical technique able to study ram pressure, Kelvin-Helmholtz instability, Rayleigh-Taylor and tidal forces acting on the star formation processes in molecular clouds. We consider it together with synthetic colour magnitude diagrams techniques. We developed a method to investigate the connections existing between gas consumption processes and star formation processes in the context of the two extended-body interaction with special attention to the dwarf galaxies dynamical regime.

Models of chemical evolution of galaxies including the dust are nowadays required to decipher the high-z universe. In a series of three papers we have tackled the problem and set a modern chemical evolution model. In the first paper (Piovan et al., 2011a) we revised the condensation coefficients for the elements that typically are present in the dust. In the second paper (Piovan et al., 2011b) we have implemented the dust into the Padova chemical model and tested it against the observational data for the Solar Neighbourhood. In this paper we extend it to the whole Disk of the Milky Way (MW). The Disk is used as a laboratory to analyze the spatial and temporal behaviour of (i) several dust grain families with the aid of which we can describe the ISM, (ii) the abundances in the gas, dust, and total ISM of the elements present in the dust and (iii) the depletion of the same elements. The temporal evolution of the dust and gas across the Disk is calculated under the effect of radial flows and a central Bar. The gradients of the abundances of C, N, O, Mg, Si, S, Ca and Fe in gas and dust across the Disk are derived as a function of time. The theoretical gradients nicely reproduce those derived from Cepheids, OB stars, Red Giants and HII regions. This provides the backbone for the companion processes of dust formation and evolution across the Disk. We examine in detail the contributions to dust by AGB stars, SNae and grain accretion in the ISM at different galacto-centric distances. Furthermore, we examine the variation of the ratio silicates/carbonaceous grains with time and position in the Disk. Finally, some hints about the depletion of the elements in regions of high and low SFR (inner and outer Disk) are presented. The results obtained make it possible to extend the model to other astrophysical situations or different theoretical models like the chemo-dynamical N-Body simulations.

The growing interest in the high-z universe, where strongly obscured objects are present, has determined an effort to improve the simulations of dust formation and evolution in galaxies. Three main basic ingredients enter the problem influencing the total dust budget and the mixture of the dust grains: the types and amounts of dust injected by AGB stars and SNae and the accretion/destruction processes of dust in the ISM. They govern the relative abundances of the gas and dust in the ISM. In this study, we focus on star-dust and present a database of condensation efficiencies for the refractory elements C, O, Mg, Si, S, Ca and Fe in AGB stars and SNae that can be easily applied to the gaseous ejecta, in order to determine the amount and kind of refractory elements locally embedded into dust and injected into the ISM. The best theoretical recipes available nowadays in literature to estimate the amount of dust produced by SNae and AGB stars have been discussed and for SNae compared to the observations to get clues on the problem. The condensation efficiencies have been analyzed in the context of a classical chemical model of dust formation and evolution in the Solar Neighbourhood and Galactic Disk. Tables of coefficients are presented for (i) AGB stars at varying the metallicity and (ii) SNae at varying the density n_H of the ISM where the SNa explosions took place. In particular, we show how the controversial CNT approximation widely adopted to form dust in SNae, still gives good results and agrees with some clues coming from the observations. A new generation of dust formation models in SNae is however required to solve some contradictions recently emerged. A simple database of condensation efficiencies is set up to be used in chemical models including the effect of dust and meant to simulate galaxies going from primordial proto-galaxies to those currently seen in the local universe.

Over the past decade a new generation of chemical models have included the dust in the treatment of the ISM. This major accomplishment has been spurred by the growing amounts of data on the highly obscured high-z Universe and the intriguing local properties of the Solar Neighbourhood (SoNE). We present here a new model able to simulate the formation and evolution of dust in the ISM. The model follows the evolution of 16 elemental species, with particular attention to those that are simultaneously present in form of gas and dust, e.g. C, N, O, Mg, Si, S, Ca and Fe. In this study we focus on the SoNe and the MW Disk as a whole which are considered as laboratories to test the physical ingredients governing the dust evolution. Infall of primordial gas, birth and death of stars, radial flows of matter between contiguous shells, presence of a central bar, star-dust emission by SNae and AGB stars, dust destruction and accretion are taken into account. The model reproduces the local depletion of the elements in the gas, and simultaneously satisfies other constraints obtained from the observations. The evolution of the element abundances in the gas and dust has been well reproduced for plausible choices of the parameters. The Mg/Si ratio, in particular, drives the formation of silicates. We show that for most of the evolution of the MW, the main process for dust enrichment is the accretion in the cold regions of the ISM. SNae dominate in the early phases of the evolution. We have examined the main factors controlling the temporal window in which SNae govern the dust budget both in low and high star forming environments. The role played by AGB stars is discussed. We find that IMFs with regular slope in the range of massive stars better reproduce the observed depletions. The results obtained for the SoNe lead us to safely extend the model.

We present ROBO, a model and its companion code for the study of the interstellar medium (ISM). The aim is to provide an accurate description of the physical evolution of the ISM and to set the ground for an ancillary tool to be inserted in NBody-Tree-SPH (NB-TSPH) simulations of large scale structures in cosmological context or of the formation and evolution of individual galaxies. The ISM model consists of gas and dust. The gas chemical composition is regulated by a network of reactions that includes a large number of species (hydrogen and deuterium based molecules, helium, and metals). New reaction rates for the charge transfer in $\mathrm H^+$ and $\mathrm H_2$ collisions are presented. The dust contains the standard mixture of carbonaceous grains (graphite grains and PAHs) and silicates of which the model follows the formation and destruction by several processes. The model takes into account an accurate treatment of the cooling process, based on several physical mechanisms, and cooling functions recently reported in the literature. The model is applied to a wide range of the input parameters and the results for important quantities describing the physical state of the gas and dust are presented. The results are organized in a database suited to the artificial neural networks (ANNs). Once trained, the ANNs yield the same results obtained by ROBO, with great accuracy. We plan to develop ANNs suitably tailored for applications to NB-TSPH simulations of cosmological structures and/or galaxies.

Modelling the complex physics of the Interstellar Medium (ISM) in the context of large-scale numerical simulations is a challenging task. A number of methods have been proposed to embed a description of the ISM into different codes. We propose a new way to achieve this task: Artificial Neural Networks (ANNs). The ANN has been trained on a pre-compiled model database, and its predictions have been compared to the expected theoretical ones, finding good agreement both in static and in dynamical tests run using the Padova Tree-SPH code \textsc{EvoL}. A neural network can reproduce the details of the interstellar gas evolution, requiring limited computational resources. We suggest that such an algorithm can replace a real-time calculation of mass elements chemical evolution in hydrodynamical codes.

We present a new study of the evolution of the Carina dwarf galaxy that includes a simultaneous derivation of its orbit and star formation history. The structure of the galaxy is constrained through orbital parameters derived from the observed distance, proper motions, radial velocity and star formation history. The different orbits admitted by the large proper motion errors are investigated in relation to the tidal force exerted by an external potential representing the Milky Way (MW). Our analysis is performed with the aid of fully consistent N-body simulations that are able to follow the dynamics and the stellar evolution of the dwarf system in order to determine self-consistently the star formation history of Carina. We find a star formation history characterized by several bursts, partially matching the observational expectation. We find also compatible results between dynamical projected quantities and the observational constraints. The possibility of a past interaction between Carina and the Magellanic Clouds is also separately considered and deemed unlikely.

One of the major challenges in modern astrophysics is to understand the origin and the evolution of galaxies, the bright, massive early type galaxies (ETGs) in particular. Therefore, these galaxies are likely to be good probes of galaxy evolution, star formation and, metal enrichment in the early Universe. In this context it is very important to set up a diagnostic tool able to combine results from chemo-dynamical N-Body-TSPH (NB-TSPH) simulations of ETGs with those of spectro-photometric population synthesis and evolution so that all key properties of galaxies can be investigated. The main goal of this paper is to provide a preliminary validation of the software package before applying it to the analysis of observational data. The galaxy models in use where calculated by the Padova group in two different cosmological scenarios: the SCDM, and the Lambda CDM. For these models, we recover their spectro-photometric evolution through the entire history of the Universe. We computed magnitudes and colors and their evolution with the redshift along with the evolutionary and cosmological corrections for the model galaxies at our disposal, and compared them with data for ETGs taken from the COSMOS and the GOODS databases. Starting from the dynamical simulations and photometric models at our disposal, we created synthetic images from which we derived the structural and morphological parameters. The theoretical results are compared with observational data of ETGs selected form the SDSS database. The simulated colors for the different cosmological scenarios follow the general trend shown by galaxies of the COSMOS and GOODS. Within the redshift range considered, all the simulated colors reproduce the observational data quite well.

Systems of shells and polar rings in early-type galaxies are considered "bona fide" tracers of mass accretion and/or mergers. Their high frequency in low density environments suggests that such episodes could drive the evolution of at least a fraction of the early-type galaxy population. Their UV emission is crucial to test whether these galaxies host ongoing/recent star formation. We used far and near ultraviolet, optical, near infrared images, HI maps, and line strength indices to investigate the nuclear and outer regions of the galaxies as well as the regions where fine structures are present. The GALEX Near (NUV) and Far UV (FUV) images of MCG-05-07-001 and NGC 1210 show complex tidal tails and debris structures. The UV morphology of both galaxies appears so different from the optical one that the early-type classification may not apply. In both GALEX bands the polar ring of MCG-05-07-001 is the dominant feature, whereas an extended tidal tail dominates the FUV bands of NGC 1210. In MCG-05-07-001 and NGC 1210 there is a strong correlation between structures detected in the FUV and NUV bands and in HI. NGC 5329 does not show evidence of shells in the UV. We try to constrain the age of the accretion episode or merger which gave rise to the shells and polar rings with the aid of composite stellar populations that take the presence of dust into account. The presence of HI in both MCG-05-07-001 and NGC 1210 argues in favour of wet mergers. Models suggest the presence of very young stellar populations in MCG-05-07-001: the observations could be explained in the framework of a conspicuous burst of star formation that occurred <=1 Gyr ago and involved a large fraction of the galaxy mass. Our models suggest that also the nuclei of NGC 1210 and NGC 5329 could have been rejuvenated by an accretion episode about 2-4 Gyr ago. (abridged)

We present EvoL, the new release of the Padova N-body code for cosmological simulations of galaxy formation and evolution. In this paper, the basic Tree + SPH code is presented and analysed, together with an overview on the software architectures. EvoL is a flexible parallel Fortran95 code, specifically designed for simulations of cosmological structure formation on cluster, galactic and sub-galactic scales. EvoL is a fully Lagrangian self-adaptive code, based on the classical Oct-tree and on the Smoothed Particle Hydrodynamics algorithm. It includes special features such as adaptive softening lengths with correcting extra-terms, and modern formulations of SPH and artificial viscosity. It is designed to be run in parallel on multiple CPUs to optimize the performance and save computational time. We describe the code in detail, and present the results of a number of standard hydrodynamical tests.

The spatial distribution of galaxies in the Local Group (LG) is the footprint of its formation mechanism and the gravitational interactions among its members and the external massive galaxies or galaxy groups. Recently, Pasetto & Chiosi (2007), using a 3D-geometrical description of the spatial distribution of all the members of the LG (not only the satellites of the MW and M31) based on present-day data on positions and distances, found that all galaxies (MW, M31, their satellites, and even the most distant objects) are confined within a slab of about 200 kpc thickness. Examining how external galaxies or groups would gravitationally affect (and eventually alter) the planar structure (and its temporal evolution) of the LG, they found that the external force field acts parallel to the plane determined by geometry and studied this with the Least Action Principle. In this paper, we have thoroughly investigated the role played by the tidal forces exerted by external galaxies or galaxy groups on the LG galaxies (the most distant dwarfs in particular) in shaping their large scale distribution. The idea based on the well known effect of tidal interactions, according to which a system of mass-points can undergo not only tidal stripping but also tidal compression and thus become flatter. Excluding the dwarf galaxies tightly bound to the MW and M31, the same tidal forces can account for the planar distribution of the remaining dwarf galaxies. We analytically recover the results of Pasetto & Chiosi (2007) and prove that a planar distribution of the LG dwarf galaxies is compatible with the external force field. We also highlight the physical cause of this result.

We present XMM-Newton X-ray observations of two shell galaxies, NGC 7070A and ESO 2400100, and far UV observations obtained with the Optical Monitor for these and for an additional shell galaxy, NGC 474, for which we also have near and far UV data from GALEX. We aim at gaining insight on the overall evolution traced by their star formation history and by their hot gas content. The X-ray and the far UV data are used to derive their X-ray spatial and spectral characteristics and their UV luminosity profiles. We use models developed ad hoc to investigate the age of the last episode of star formation from the (UV - optical) colors and line strength indices. The X-ray spatial and spectral analysis show significant differences in the two objects. A low luminosity nuclear source is the dominant component in NGC 7070A log L_X=41.7 erg s^{-1} in the 2-10 keV band. In ESO 2400100, the X-ray emission is due to a low temperature plasma with a contribution from the collective emission of individual sources. In the Optical Monitor image ESO 2400100 shows a double nucleus, one bluer than the other. This probably results from a very recent star formation event in the northern nuclear region. The extension of the UV emission is consistent with the optical extent for all galaxies, at different degrees of significance in different filters. The presence of the double nucleus, corroborated by the (UV - optical) colors and line strength indices analysis, suggests that ESO 2400100 is accreting a faint companion. We explore the evolution of the X-ray luminosity during accretion processes with time. We discuss the link between the presence of gas and age, since gas is detected either before coalescence or several Gyr (>3) after (Abridged).

Shell galaxies are widely considered the debris of recent accretion/merging episodes. Their high frequency in low density environment suggests that such episodes could be among the driver of the early-type galaxy secular evolution. We present far and near UV (FUV and NUV respectively hereafter) GALEX photometric properties of a sample of shell galaxies.

We present the GALEX far FUV and near NUV ultraviolet imaging of three nearby shell galaxies, namely NGC 2865, NGC 5018 and NGC 7135. The system of shells and fine structures visible in the optical is detected in the NUV image of NGC 2865 and in both NUV and FUV images of NGC 7135. The NUV image of NGC 5018 does not present shell structures. We detect absorption features in the nuclear region of all three galaxies. NGC 2865 has a nearly flat colour profile with (FUV-NUV)~2 throughout the whole galaxy. NGC 7135 is blue in the center (FUV-NUV)~0 and as red as (FUV-NUV)~1.5 in the outskirts, including the faint shell-like feature. We investigate the ability of the nuclear GALEX (FUV-NUV) colour to provide information about rejuvenation phenomena in the stellar populations of the shell galaxies. To this aim, we derive from theory the relationship between the Mg2, Hbeta, HgammaA, HdeltaA Lick line-strength indices and the (FUV-NUV) colour. We extend the study to a sample of early-type galaxies with emission lines in their optical spectra (Annibali et al. 2007). In the index vs.(FUV-NUV) colour diagrams, most of the galaxies are well explained by passively evolving SSPs. On the average, ages and metallicities of the galaxies in our sample estimated from optical line-strength indices are consistent with those inferred from the (FUV-NUV) colour. In general, all the colours but for (FUV-NUV) and (FUV-V), become nearly age insensitive when 1-2 Gyr have elapsed from the last star forming event. Considering composite stellar population models with a recent burst of star formation, we suggest that the position of the NGC 7135 and NGC 2865 nuclei in the (FUV-NUV)-Hbeta plane could be explained in term of a recent rejuvenation episode. (Abridged)

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.

We present the preliminary results of a two phases description of the ISM suited to numerical N Body T-SPH models of galaxy formation and evolution.

Adopting known data on positions and distances, we make use of the analytical geometry and look for the plane that minimizes the distances of all galaxies to it. A planar distribution is indeed found that, however, does not coincide with the plane found by Sawa & Fujimoto 2005. Why? The second part of this study is devoted to answer this question and to find a dynamical justification for the planar distribution. To this aim, we apply the Hamilton Method (Minimum Action) to investigate the dynamics of the two major system of the Local Group, Milky Way and Andromeda, under the action of external forces exerted by nearby galaxies or groups external to the Local Group. We find that the planar istribution is fully compatible with the minimum action and that the external force field is likely parallel to the plane. It pulls the galaxies of the Local Group without altering their planar distribution. Special care is paid to evaluate the robustness of this result. Conclusion: in this paper we have examined the spatial distribution of galaxies in the Local Group. They are confined to a plane that can be statistically and dynamically understood as the result of the Minimum Action. The planar distribution seems to be stable for a large fraction of the Hubble time. The external force field, that has likely been constant over the same time interval, does not alter the planar distribution as it is nearly parallel to it. Effects due to undetected halos of sole Dark Matter are briefly discussed. They could be a point of uncertainty of the present study.

Basing on the huge library of 1-A resolution spectra calculated by Munari et al. over a large range of logT, log g, [Fe/H] and both for solar and a-enhanced abundance ratios [a/Fe], we present theoretical absorption-line indices on the Lick system. First we derive the so-called response functions (RFs) of Tripicco & Bell for a wide range of logT, log g, [Fe/H] and [a/Fe]=+0.4 dex. The RFs are commonly used to correct indices with solar [a/Fe] ratios to indices with [a/Fe]>0. Not only the RFs vary with the type of star but also with the metallicity. Secondly, with the aid of this and the fitting functions (FFs) of Worthey et al., we derive the indices for single stellar populations and compare them with those obtained by previous authors, e.g. Tantalo & Chiosi. The new RFs not only supersede the old ones by Tripicco & Bell, but also show that Hb increases with the degree of enhancement in agreement with the results by Tantalo & Chiosi. The new indices for single stellar populations are used to derive with aid of the recursive Minimum Distance method the age, metallicity and degree of enhancement of a sample of Galactic Globular Clusters for which these key parameters have been independently derived from the Colour-Magnitude Diagram and/or spectroscopic studies. The agreement is remarkably good.

In this study, with the aid of N-Body simulations based on quasi-cosmological initial conditions, we have followed the formation and evolution of two models of early-type galaxies, from their separation from global expansion of the Universe to their collapse to virialized structures, the formation of stars and subsequent nearly passive evolution. The cosmological background we have considered is the Standard CDM. The models have significantly different nitial total mass. Particular care has been paid to the star formation process, heating and cooling of gas, and chemical enrichment. In both models star formation is completed within the first Gyrs of evolution. The structural properties of the present-day models are in good agreement with current observations. The chemical properties, mean metallicity and metallicity gradients also agree with available observational data. Finally, conspicuous galactic winds are found to occur. The models conform to the so-called revised monolithic scheme, because mergers of substructures have occurred very early in the galaxy life. Our results agree with those obtained in other similar recent studies, thus strengthening the idea that the revised monolithic scheme is the right trail to follow in the forest of galaxy formation and evolution.

In this paper, we present spectrophotometric models for galaxies of different morphological type whose spectral energy distributions (SEDs) take into account the effect of dust in absorbing UV-optical light and re-emitting it in the infrared (IR). The models contain three main components: (i) the diffuse interstellar medium (ISM) composed of gas and dust, (ii) the large complexes of molecular clouds (MCs) in which new stars are formed and (iii) the stars of any age and chemical composition. The galaxy models stand on a robust model of chemical evolution taylored in such a way to match the gross properties of galaxies of different morphological type. A simple geometrical model is adopted for each type of galaxy. The total gas and star mass provided by the chemical model are distributed over the whole volume by means of suitable density profiles, one for each component and depending on the galaxy type (spheroidal, disk and disk plus bulge). The galaxy is then split in suitable volume elements to each of which the appropriate amounts of stars, MCs and ISM are assigned. We derive the SEDs of galaxies of different morphological type. First the the basic relations driving the interaction between the physical components of the galaxy are presented. Second, the main parameters are examined and their effects on the SED of three prototype galaxies (a disk, an elliptical and a starburster) are highlighted. The final part of the paper is devoted to assess the ability of our galaxy models in reproducing the SEDs of a few real galaxies of the Local Universe.

In this paper, the first of a series of two devoted to modelling the spectra of galaxies of different morphological type in presence of dust, we present our description of the dust both in the diffuse ISM and the molecular clouds. Our model for the dust takes into account three components, i.e. graphite, silicates and polycyclic aromatic hydrocarbons. We consider and adapt to our aims two prescriptions for the size distribution of the dust grains and two models for the emission of the dusty ISM. We cross-check the emission and extinction models of the ISM by calculating the extinction curves and the emission for the typical environments of the Milky Way and the Large and Small Magellanic Clouds and by comparing the results with the observational data. The final model we have adopted is an hybrid one which stems from combining the analysis of Guhathakurta et al (1989) for the emission of graphite and silicates and Puget et al. (1985) for the PAH emission, and using the distribution law of Weingartner et al. (2001a) and the ionization model for PAHs of Weingartner et al. (2001b). We apply the model to calculate the SED of SSPs of different age and chemical composition, which may be severely affected by dust in young, massive stars while they are still embedded in their parental MCs. We use the "Ray Tracing" method to solve the problem of radiative transfer and to calculate libraries of SSP SEDs. Particular care is paid to model the contribution from PAHs, introducing different abundances of C in the population of very small carbonaceous grains (VSGs) and different ionization states in PAHs. The SEDs of young SSPs are then compared with observational data of star forming regions of four local galaxies successfully reproducing their SEDs from the UV-optical regions to the mid and far infrared.

In this study we present detailed calculations of absorption line indices on the Lick System based on the stellar models by Salasnich et al. (2000) calculated with enhanced mix of alpha-elements. Using the so-called Response Functions (RFs) of Tripicco & Bell (1995, TB95), we calculate the indices for SSPs of different age, metallicity and enhancement. We made use of the triplet Hb, Mgb and <Fe>, and Minimum-Distance Method proposed by Trager et al. (2000, TFWG00) to estimate the age, metallicity and enhancement degree for the galaxies of the Gonzalez (1993) sample, and compare the results with those TFWG00 and Thomas et al (2003). Since very large differences are found, in particular as far as the age is concerned, we analyze in a great detail all possible sources of disagreement, going from the stellar models and SSPs to many technical details of the procedure to calculate the indices, and finally the pattern of chemical elements (especially when alpha-enhanced mixtures are adopted). The key issue of the analysis is that at given metallicity Z and enhancement factor, the specific abundance ratios [Xel/Fe] adopted for some elements (e.g. O, Mg, Ti, and likely others) dominate the scene because with the TB95 RFs they may strongly affect indices like Hb and the age in turn. Finally we have drawn some remarks on the interpretation of the distribution of early-type galaxies in popular two-indices planes, like Hb vs. [MgFe]. We argue that part of the scatter along the Hb axis observed in this plane could be attributed instead of the age, the current explanation, to a spread both in the degree of enhancement and some abundance ratios. The main conclusion of this study is that deriving ages, metallicities and degree of enhancement from line indices is a cumbersome affair whose results are still uncertain.

The bottom line of this study is (i) to check whether normal (quiescent) and interacting EGs have a different behavior in the popular diagnostic planes such as Hb vs. [MgFe] (and others); (ii) to seek whether the signature of mergers and/or interactions may mirror in some specific changes of the indices that could un-equivocally hint for burts of stellar activity; (iii) to evaluate the intensity of those bursts or secondary episodes of star formation; (iv) to explore whether other alternatives can exist, i.e. distinct from obvious ones resting on large age range. The data show that both normal, field and interacting galaxies have the same scattered but smooth distribution in the Hb vs. [MgFe] diagnostic plane. There are two interpretations for this: (1) EGs span larges ranges of ages and metallicities. (2) The bulk population of stars is old but subsequent episodes of star formation scatter the object in the diagnostic planes. Both these interpretation are not able to reproduce the observed distribution. The above schemes are too demanding because of the many ad hoc ingredients that have to be introduced. Furthermore they neglect the observationally grounded hint that the stellar content of EGs is likely enhanced in alfa-elements. We propose here a new scheme, in which the bulk dispersion of galaxies in the index-index plane is caused by a different mean degree of enhancement. The possibility that EGs span large ranges of [alpha/Fe] and metallicities but narrow ranges of ages for the bulk population favors the monolithic scheme.

We investigate the bright cut-off of the [OIII]l5007 planetary nebula luminosity function (PNLF), that has been suggested as a powerful extragalactic distance indicator. Theoretical PNLFs are constructed via Monte-Carlo simulations of populations of PNe, whose individual properties are described with the aid of recent PN synthetic models (Marigo et al. 2001), coupled to a detailed photoionisation code (CLOUDY). The basic dependences of the cut-off magnitude M* are then discussed. We find that: (i) In galaxies with recent or ongoing star formation, the modelled PNLF present M* values between -4 and -5, very close to the observationally-calibrated value for the LMC. (ii) In these galaxies, the PNLF cut-off is produced by PNe with progenitor masses of about 2.5 Msun, while less massive stars give origin to fainter PNe. As a consequence M* is expected to depend strongly on the age of the last burst of star formation, dimming by as much as 5 mag as we go from young to 10-Gyr old populations. (iii) Rather than on the initial metallicity of a stellar population, M* depends on the actual [O/H] of the observed PNe, a quantity that may differ significantly from the initial value (due to dredge-up episodes), especially in young and intermediate-age PN populations. (iv) Also the transition time from the end of AGB to the PN phase, and the nuclear-burning properties (i.e. H- or He-burning) of the central stars introduce non-negligible effects on M*. The strongest indication derived from the present calculations is a serious difficulty to explain the age-invariance of the cut-off brightness over an extended interval, say from 1 to 13 Gyr, that observations of PNLFs in galaxies of late-to-early type seem to suggest.

Following Paper I, we provide extended tables of bolometric corrections, extinction coefficients, stellar isochrones, and integrated magnitudes and colours of single-burst stellar populations, for the Sloan Digital Sky Survey (SDSS) ugriz photometric system. They are tested on comparisons with DR1 data for a few stellar systems, namely the Palomar 5 and NGC 2419 globular clusters and the Draco dSph galaxy.

It is frequently debated in literature whether a "standard" Initial Mass Function (IMF) - meaning an IMF of the kind usually adopted to explain the chemical evolution in the local Solar Neighbourhood - can account for the observed metal enrichment and Iron Mass-to-Light Ratio in clusters of galaxies. We address this problem by means of straightforward estimates that should hold independently of the details of chemical evolution models. It is crucial to compute self-consistently the amount of mass and metals locked-up in stars, by accounting for the stellar mass-to-light ratio predicted by a given IMF. It becomes then clear that a "standard" Solar Neighbourhood IMF cannot provide enough metals to account for the observed chemical properties in clusters: clusters of galaxies and the local environment must be characterized by different IMFs. Alternatively, if we require the IMF to be universal, in order to explain clusters such an IMF must be much more efficient in metal production than usually estimated for the Solar Vicinity. In this case, substantial loss of metals is required from the Solar Neighbourhood and from disc galaxies in general. This "non-standard" scenario for the local chemical evolution would challenge our present understanding of the Milky Way and of disc galaxy formation.

The dwarf galaxies of the Local Group can be separated in three morphological groups: irregular, elliptical and spheroidal. As in the large galaxy clusters, there seems to be a morphology-position relationship: irregular galaxies are preferentially found in the outskirts (low density regions) of the Local Group, whereas dwarf ellipticals and spheroidals are more frequent in the central, high density regions. To cast light on the nature and origin of dwarf galaxies in the Local Group, Mayer et al.(2001) have suggested that a dwarf irregular galaxy tidally interacting with a galaxy of much larger mass may be re-shaped into a dwarf spheroidal or elliptical object. In this paper we check by means of N-Body Tree-SPH simulations whether this is possible for a selected sample of galaxies of the Local Group. Using the best data available in literature to fix the dynamical and kinematical status of a few dwarf galaxies in the Local Group, we follow the evolution of an ideal satellite, which supposedly started as an irregular object during its orbital motion around the Milky Way. We find that the tidal interactions with the Milky Way remove a large fraction of the mass of the dwarf irregular and gradually reshape it into a spherical object.

In this paper we present models for SSPs of intermediate and old ages where dust enshrouded AGB stars are introduced. To this aim, particular care is devoted to follow the evolution of the AGB stars throughout the quiet and thermally pulsing regimes, to evaluate the effect of self contamination in the outermost layers by the third dredge-up mechanism, to follow the transition from oxygen-rich to carbon-rich objects, and finally to estimate the efficiency of mass-loss by stellar winds. In addition to this, accurate physical models of the dusty shells are presented in which the re-processing of radiation from the central stars is calculated. The resulting spectral energy distribution (SED) is examined to show how important features evolve with time. The SEDs are then convolved with the IRAS filters to obtain the flux in various pass-bands for individual AGB stars of different, mass, chemical composition, and age. The comparison is made by means of SSPs along which AGB stars of the same age but different initial masses are located. The theoretical results are compared to the observational data for selected groups of stars. The same is made for the J,H,K,L pass-bands of the Johnson system. Finally, from the integrated SEDs of the SSPs, we derive the integrated Johnson J,H,K,L magnitudes and colors to be compared to infrared data for star clusters of the Magellanic Clouds.

The high metallicity of the intra-cluster medium (ICM) is generally interpreted on the base of the galactic wind scenario for elliptical galaxies. In this framework, we develop a toy-model to follow the chemical evolution of the ICM, formulated in analogy to chemical models for individual galaxies. The model computes the galaxy formation history (GFH) of cluster galaxies, connecting the final luminosity function (LF) to the corresponding metal enrichment history of the ICM. The observed LF can be reproduced with a smooth, Madau-plot like GFH peaking at z~ 1-2, plus a "burst" of formation of dwarf galaxies at high redshift. The model is used to test the response of the predicted metal content and abundance evolution of the ICM to varying input galactic models. The chemical enrichment is computed from "galactic yields" based on models of elliptical galaxies with a variable initial mass function (IMF), favouring the formation of massive stars at high redshift and/or in more massive galaxies. For a given final galactic luminosity, these model ellipticals eject into the ICM a larger quantity of gas and of metals than do standard models based on the Salpeter IMF. However, a scenario in which the IMF varies with redshift as a consequence of the effect of the the cosmic background temperature on the Jeans mass scale, appears to be too mild to account for the observed metal production in clusters. The high iron-mass-to-luminosity-ratio of the ICM can be reproduced only by assuming a more dramatic variation of the typical stellar mass, in line with other recent findings. The mass in the wind-ejected gas is predicted to exceed the mass in galaxies by a factor of 1.5-2 and to constitute roughly half of the intra-cluster gas.

Carbon stars are known to exhibit systematically redder near-infrared colours with respect to M-type stars. In the near-infrared colour-magnitude diagrams provided by the 2MASS and DENIS surveys, the LMC C-type stars draw a striking red tail, well separated from the sequences of O-rich giants. So far, this conspicuous feature has been absent from any set of available isochrones, even the few existing ones that include the TP-AGB evolution of low- and intermediate-mass stars. To investigate such issue we simulate the complete 2MASS Ks vs.(J-Ks) data towards the LMC by means of a population synthesis approach, that relies on extended libraries of published stellar evolutionary tracks, including the TP-AGB phase. The simulations provide quite a detailed description of the several vertical fingers and inclined sequences seen in 2MASS data, due to both Galactic foreground and LMC O-rich stars. Instead, as mentioned, the red tail of C-stars sets a major difficulty: we find that TP-AGB models with solar-scaled molecular opacities, the usual assumption of existing AGB calculations, do not succeed in reproducing this feature. Our tests indicate that the main reason for this failure should not be ascribed to empirical Teff - (J-K) transformations for C-type stars. Instead, the discrepancy is simply removed by adopting new evolutionary models that account for the changes in molecular opacities as AGB stars get enriched in carbon via the third dredge-up (Marigo 2002). In fact, simulations that adopt these models are able to reproduce, for the first time, the red tail of C-stars in near-infrared CMDs. Finally, we point out that these simulations also provide useful indications about the efficiency of the third dredge-up process, and the pulsation modes of long-period variables.

GAIA data will create a precise 3-dimensional map of the Galaxy, providing positional information, radial velocities, luminosity, temperature and chemical composition of a representative sample of stars. Here we present a new implementation of the Padova Galaxy Model, where kinematics simulations are included. A few examples of application to GAIA science are discussed, in particular concerning radial velocities determinations of the thin and thick disk.

We present evolutionary models of zero-metallicity very massive objects, with initial masses in the range 120 Msun -- 1000 Msun, covering their quiescent evolution up to central carbon ignition. In the attempt of exploring the possible occurrence of mass loss by stellar winds, calculations are carried out with recently-developed formalisms for the mass-loss rates driven by radiation pressure (Kudritzki 2002) and stellar rotation (Maeder & Meynet 2000).The study completes the previous analysis by Marigo et al. (2001) on the constant-mass evolution of primordial stars. Our results indicate that radiation pressure (assuming a minimum metallicity Z = 10^{-4} Zsun)is not an efficient driving force of mass loss, except for very massive stars with M >= 750 Msun. On the other hand, stellar rotation might play a crucial role in triggering powerful stellar winds, once the (Omega-Gamma)-limit is approached. However, this critical condition of intense mass loss can be maintained just for short, as the loss of angular momentum due to mass ejection quickly leads to the spinning down of the star. As by-product to the present work, the wind chemical yields from massive zero-metallicity stars are presented. The helium and metal enrichments, and the resulting Delta(Y)/Delta(Z) ratio are briefly discussed.

The color-magnitude diagrams (CMDs) of three intermediate-age LMC clusters, NGC 2173, SL556 and NGC2155 are analyzed to determine their age and metallicity basing on Padova stellar models. Synthetic CMDs are compared with cluster data. The best match is obtained using two fitting functions based on star counts in the different bins of the cluster CMD. Two different criteria are used. One of them takes into account the uncertainties in the color of the red clump stars. Given the uncertainties on the experimental values of the clusters metallicity, we provide a set of acceptable solutions. They define the correspondent values of metallicity, age, reddening and distance modulus (for the assumed IMF). The comparison with Padova models suggests for NGC 2173 a prolonged star formation (spanning a period of about 0.3 Gyr), beginning 1.7 Gyr and ending 1.4 Gyr ago. The metallicity Z is in the range 0.0016 $-$ 0.003. Contrary to what suggested for NGC 2173 a period of extended star formation was not required to fit the SL 556 and the NGC 2155 observations. For SL 556 an age of 2.0 Gyr is obtained. The metallicity value is in the range 0.002 $-$ 0.004, depending on the adopted comparison criterium. The derived age for NGC 2155 is 2.8 Gyr and its metallicity Z is in the range 0.002 $-$ 0.003. The CMD features of this cluster suggest that a more efficient overshoot should be adopted in the evolutionary models.

We discuss the evolutionary properties of primordial massive and very massive stars, supposed to have formed from metal-free gas. Stellar models are presented over a large range of initial masses (8 Msun <= Mi <= 1000 Msun), covering the hydrogen- and helium-burning phases up to the onset of carbon burning. In most cases the evolution is followed at constant mass. To estimate the possible effect of mass loss via stellar winds, recent analytic formalisms for the mass-loss rates are applied to the very massive models (Mi >= 120 Msun).

This is the first of a series of three papers devoted to the calibration of a few parameters of crucial importance in the modeling of the evolution of intermediate-mass stars, with special attention to the amount of convective core overshoot. To this end we acquired deep V and R photometry for three globular clusters of the Large Magellanic Cloud (LMC), namely NGC 2173, SL 556 and NGC 2155, in the age interval 1-3 Gyr. In this first paper, we describe the aim of the project, the VLT observations and data reduction, and we make preliminary comparisons of the color-magnitude diagrams with both Padova and Yonsei-Yale isochrones. Two following papers in this series present the results of a detailed analysis of these data, independently carried out by members of the Yale and Padova stellar evolution groups. This allows us to compare both sets of models and discuss their main differences, as well as the systematic effects that they would have to the determination of the ages and metallicities of intermediate-age single stellar populations.

A new method, AMORE - based on a genetic algorithm optimizer, is presented for the automated study of colour-magnitude diagrams. The method combines several stellar population synthesis tools developed in the last decade by or in collaboration with the Padova group. Our method is able to recover, within the uncertainties, the parameters -- distance, extinction, age, metallicity, index of a power-law initial mass function and the index of an exponential star formation rate -- from a reference synthetic stellar population. No a priori information is inserted to recover the parameters, which is done simultaneously and not one at a time. Examples are given to demonstrate and to better understand biases in the results, if one of the input parameters is deliberately set fixed to a non-optimum value.

We provide tables of theoretical isochrones in several photometric systems. To this aim, the following steps are followed: (1) First, we re-write the formalism for converting synthetic stellar spectra into tables of bolometric corrections. The resulting formulas can be applied to any photometric system, provided that the zero-points are specified by means of either ABmag, STmag, VEGAmag, or a standard star system that includes well-known spectrophotometric standards. Interstellar absorption can be considered in a self-consistent way. (2) We assemble an extended and updated library of stellar intrinsic spectra. It is mostly based on non-overshooting ATLAS9 models, suitably extended to both low and high effective temperatures. This offers an excellent coverage of the parameter space of Teff, logg, and [M/H]. We briefly discuss the main uncertainties and points still deserving more improvement. (3) From the spectral library, we derive tables of bolometric corrections for Johnson-Cousins-Glass, HST/WFPC2, HST/NICMOS, Washington, and ESO Imaging Survey systems (this latter consisting on the WFI, EMMI, and SOFI filter sets). (4) These tables are used to convert several sets of Padova isochrones into the corresponding absolute magnitudes and colours, thus providing a useful database for several astrophysical applications. All data files are made available in electronic form.

In this paper we present a quantitative analysis of the chemical abundances in "monolithic" N-body-Tree-SPH models of elliptical galaxies to determine whether the ratio [alpha/Fe] increases with the galaxy mass as suggested by observational data.

The young LMC cluster NGC 1866 represents the ideal laboratory for testing stellar evolutionary models: in particular, it can be used to discriminate among classical stellar models, in which the extension of the convective regions is fixed by the classical Schwarzschild criterion, from models with overshooting, in which an ``extra-mixing'' is considered to take place beyond the classical limit of the convective zone. Addressing this subject in a recent work, Testa et al. (1999) reached the conclusion that the classical scheme for the treatment of convection represents a good and sufficient approximation for convective interiors. Using their own data, we repeat here the analysis. First we revise the procedure followed by Testa et al. (1999) to correct the data for completeness, second we calculate new stellar models with updated physical input for both evolutionary schemes, finally we present many simulations of the colour-magnitude diagrams and luminosity functions of the cluster using the ratio of the integrated luminosity function of main sequence stars to the number of giants as the normalization factor of the simulations. We also take into account several possible physical agents that could alter the color-magnitude diagram and the luminosity function: they are unresolved binary stars, dispersion in the age, stochastic effects in the initial mass function. Their effect is analyzed separately, with the conclusion that binary stars have the largest impact. The main result of this study is that the convective overshoot hypothesis (together with a suitable percentage of unresolved binaries) is really needed to fully match the whole pattern of data. The main drawback of the classical models is that they cannot reproduce the correct ratio of main sequence to post-main sequence stars.

We develop a toy-model for the chemical evolution of the intracluster medium, polluted by the galactic winds from elliptical galaxies. The model follows the "galaxy formation history" of cluster galaxies, constrained by the observed luminosity function.

The high metallicity of the intra-cluster medium (ICM) is generally interpreted on the base of the galactic wind scenario for elliptical galaxies. In this framework, we develop a toy-model to follow the chemical evolution of the ICM, formulated in analogy to chemical models for individual galaxies. Just as the ingredients for usual models are (a) the stellar yields, amount of metals newly synthesized and re-ejected by stars; (b) the Star Formation Rate and (c) the stellar Initial Mass Function (IMF), our model for clusters involves: (a') ``galactic yields'' derived from galactic wind models of ellipticals; (b') a parametric Galactic Formation Rate; (c') a Press-Schechter-like Galactic Initial Mass Function. The model is used to test the response of the predicted metal content and abundance evolution of the ICM to varying input galactic models. The resulting luminosity function of cluster galaxies is also calculated, in order to constrain model parameters.

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.

We present extensive evolutionary models of stars with initial zero-metallicity, covering a large range of initial masses (i.e. 0.7 <= M <= 100 Msun). Calculations are carried out at constant mass, with updated input physics, and applying an overshooting scheme to convective boundaries. The nuclear network includes all the important reactions of the p-p chain, CNO-cycle and alpha-captures, and is solved by means of a suitable semi-implicit method. The evolution is followed up to the thermally pulsing AGB in the case of low- and intermediate-mass stars, or to the onset of carbon burning in massive stars. The main evolutionary features of these models are discussed, also in comparison with models of non-zero metallicity. Among several interesting aspects, particular attention has been paid to describe: i) the first synthesis of 12C inside the stars, that may suddenly trigger the CNO-cycle causing particular evolutionary features; ii) the pollution of the stellar surface by the dredge-up events, that are effective only within particular mass ranges; iii) the mass limits which conventionally define the classes of low-, intermediate-, and high-mass stars on the basis of common evolutionary properties, including the upper mass limit for the achievement of super-Eddington luminosities before C-ignition in the high-mass regime; and iv) the expected pulsational properties of zero-metallicity stars. All relevant information referring to the evolutionary tracks and isochrones is made available in computer-readable format at http://pleiadi.pd.astro.it .

We present four large sets of evolutionary tracks for stars with initial chemical compositions [Y=0.250, Z=0.008], [Y=0.273, Z=0.019], [Y=0.320, Z=0.040] and [Y=0.390, Z=0.070] and enhancement of alpha elements with respect to the solar pattern. The major improvement with respect to previous similar calculations is that we use consistent opacities - i.e. computed with the same chemical composition as adopted in the stellar models - over the whole relevant range of temperatures. For the same initial chemical compositions [Y,Z] and otherwise identical input physics we present also new evolutionary sequences with solar-scaled mixtures of abundances. Based on these stellar models we calculate the corresponding sets of isochrones both in the Johnson-Cousins UBVRIJHK and HST/WFPC2 photometric systems. Furthermore, we derive integrated magnitudes, colours and mass-to-light ratios for ideal single stellar populations with total mass equal to 1 Msun. Finally, the major changes in the tracks, isochrones, and integrated magnitudes and colours passing from solar-scaled to alpha-enhanced mixtures are briefly outlined. Retrieval of the complete data set is possible via http://pleiadi.pd.astro.it .

We develop a numerical chemical model allowing for radial flows of gas, with the aim to analyse the possible role of gas flows in the chemical evolution of the Galactic Disc. The dynamical effects of the Galactic Bar on the radial gas profile of the Disc are especially addressed.

We present new evolutionary models for zero-metallicity stars, covering a large range of initial masses (from 0.8 to 100 M_sun). Models are computed with overshooting from stellar cores and convective envelopes, and assuming mass-loss from the most massive stars. We discuss the main evolutionary features of these stars, and provide estimates of the amount of newly-synthesized elements dredged-up to the stellar surface, and possibly lost by stellar winds from the most massive stars. Full details of these models will be given in Marigo et al. (2000, in preparation).

We present a large grid of stellar evolutionary tracks, which are suitable to modelling star clusters and galaxies by means of population synthesis. The tracks are presented for the initial chemical compositions [Z=0.0004, Y=0.23], [Z=0.001, Y=0.23], [Z=0.004, Y=0.24], [Z=0.008, Y=0.25], [Z=0.019, Y=0.273] (solar composition), and [Z=0.03, Y=0.30]. They are computed with updated opacities and equation of state, and a moderate amount of convective overshoot. The range of initial masses goes from 0.15 M_sun to 7 M_sun, and the evolutionary phases extend from the zero age main sequence (ZAMS) till either the thermally pulsing AGB regime or carbon ignition. We also present an additional set of models with solar composition, computed using the classical Schwarzschild's criterion for convective boundaries. From all these tracks, we derive the theoretical isochrones in the Johnson-Cousins UBVRIJHK broad-band photometric system.

The first UV F300W and F218W WFPC2 observations of the rich galaxy cluster CL0939+4713 at z=0.41 are presented and discussed. UV/optical two-color and c-m diagrams of the sources detected in the F300W waveband are constructed. Thanks to pre-existing HST optical images of the same field a morphological classification for the majority of these objects is also provided. Moreover, taking advantage of recent redshift surveys along CL0939+4713 line of sight, separate diagrams comparing the properties of galaxies belonging to the cluster and to its close projected field are presented. Possible evolutionary effects in the UV from z~0.4 to present time are investigated by comparing the restframe [mid-UV-Optical] colors of galaxies in CL0939+4713 with balloon-borne data of the Coma cluster, as well as by resorting to suitable galaxy evolution models. Finally, current attempts to constrain the epoch of the UV-upturn onset in evolved populations by means of HST UV observations are discussed.

In this paper we analyze the line-strength indices in the Lick-system measured by Longhetti et al. (1998a, b) for a sample of 51 early-type galaxies located in low density environments (LDE) and showing signatures of fine structures and/or interactions. The sample contains 21 shell-galaxies and 30 members of interacting pairs. Firstly we perform a preliminary comparison between three different sources of calibrations of the line strength indices, namely Buzzoni et al. (1992, 1994), Worthey (1992), Worthey et al. (1994) and Idiart et al. (1995), derived from stars with different effective temperature, gravity, and metallicity. Secondly, we discuss the properties of the galaxies in our sample by comparing them both with theoretical Single Stellar Populations (SSPs) and the {\sl normal} galaxies of the Gonz\'alez (1993: G93) sample. The analysis is performed by means of several diagnostic planes. In the H$\beta$ vs. [MgFe] plane, which is perhaps best suited to infer the age of the stellar populations, the peculiar galaxies in our sample show nearly the same distribution of the {\sl normal} galaxies in the G93 sample. There is however a number of peculiar galaxies with much stronger H$\beta$. Does this mean that the scatter in the H$\beta$ vs. [MgFe] plane, of normal, shell- and pair-galaxies has a common origin, perhaps a secondary episode of star formation? We suggest that, owing to their apparent {\sl youth}, shell- and pair-galaxies should have experienced at least one interaction event after their formation. The explanation comes natural for shell- and pair-galaxies where the signatures of interactions are evident. It is more intrigued in {\sl normal} galaxies (perhaps other causes may concur).

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.

The abundance gradients and the radial gas profile of the Galactic disc are analysed by means of a model for the chemical evolution of galaxies. As one of the major uncertainties in models for galaxy evolution is the star formation (SF) process, various SF laws are considered,to assess the response of model predictions to the different assumptions. Only some SF laws are successful in reproducing the metallicity gradient, and only if combined with a suitable infall timescale increasing outward (inside-out formation scenario). Still, it is difficult to reproduce at the same time also the observed gas distribution; we therefore suggest further improvements for the models.

We present a quantitative analysis of the star formation history (SFH) of the Local Group dSph galaxy Leo I, from the information in its HST [(V-I),I] color-magnitude diagram (CMD). The method we use is based in comparing, via synthetic CMDs, the expected distribution of stars in the CMD for different evolutionary scenarios, with the observed distribution. We consider the SFH to be composed by the SFR(t), the Z(t), the IMF, and a function $\beta(f,q)$, controlling the fraction $f$ and mass ratio distribution $q$ of binary stars. The comparison between the observed CMD and the model CMDs is done through chi-square minimization of the differences in the number of stars in a set of regions of the CMD. Our solution for the SFH of Leo I defines a minimum of chi-square in a well defined position of the parameter space, and the derived SFR(t) is robust, in the sense that its main characteristics are unchanged for different combinations of the remaining parameters. However, only a narrow range of assumptions for Z(t), IMF and $\beta(f,q)$ result in a good agreement between the data and the models, namely: Z=0.0004, a Kroupa et al. (1993) IMF or slightly steeper, and a relatively large fraction of binary stars. Most star formation activity (70% to 80%) occurred between 7 and 1 Gyr ago. At 1 Gyr ago, it abruptly dropped to a negligible value, but seems to have been active until at least ~ 300 Myr ago. Our results don't unambiguously answer the question of whether Leo I began forming stars around 15 Gyr ago, but it appears that the amount of this star formation, if existing at all, would be small.

Aim of this study is to infer the age of the Galactic Disk by means of the ages of old open clusters, and comment on some recent claims that the Galactic Disk can be older than the Halo. To this purpose, we analyze the Color-Magnitude Diagrams (CMDs) of six very old clusters, namely NGC 188, NGC 6791, Collinder 261, Melotte 66, Berkeley 39 and Berkeley 17, and determine their ages. For each cluster we use the most recent photometric and spectroscopic data and metallicity estimates. The ages are derived from the isochrone fitting method using the stellar models of the Padua library (Bertelli et al . 1994, Girardi et al. 1999). We find that the ages of these clusters fall in the range 4 to 9-10 Gyr: Melotte 66 is the youngest whereas NGC 6791 and Berkeley 17 have ages of about 9-10 Gyr. Previous estimates for Berkeley 17 indicated an age as old as 12 Gyr, almost falling within the range of classical globular clusters. In our analysis, this cluster is always very old but perhaps somewhat younger than in previous studies. However we call attention on the fact that the above ages are to be taken as provisional estimates, because of the many uncertaintes still affecting stellar models in the mass range 1.0 to 1.5 Mo. Despite this drawback of extant theory of stellar structure, if NGC 6791 and Berkeley 17 set the limit to the age of the Galactic Disk, this component of the Milky Way can be as old as 9-10 Gyr, but surely younger than the Galactic Halo, at least as inferred from recent determinations of the age of globular clusters. Finally, it is worth recalling that open clusters can only provide a lower limit to the age of the Galactic Disk, while other indicators - like White Dwarfs - are perhaps more suited to this task.

We analyze the star formation properties of a sample of 21 shell galaxies and 30 early-type galaxies members of interacting pairs, located in low density environments (Longhetti et al 1998a, 1998b). The study is based on new models developed to interpret the information coming from `blue' H$\delta$/FeI, H+K(CaII) and \D4000 line-strength indices proposed by Rose (1984; 1985) and Hamilton (1985). We find that the last star forming event that occurred in the nuclear region of shell galaxies is statistically old (from 0.1 up to several Gyr) with respect to the corresponding one in the sub-sample of pair galaxies (<0.1 Gyr or even ongoing star formation). If the stellar activity is somehow related to the formation of shells, as predicted by several dynamical models of galaxy interaction, shells have to be considered long lasting structures. Since pair members show evidence of very recent star formation, we suggest that either large reservoirs of gas have to be present to maintain active star formation, if these galaxies are on periodic orbits, or most of the pair members in the present sample are experiencing unbound encounters.

We present deep HST F555W (V) and F814W (I) observations of a central field in the Local Group dwarf spheroidal (dSph) galaxy Leo I. The resulting color-magnitude diagram (CMD) reaches I \simeq 26 and reveals the oldest ~10-15 Gyr old turnoffs. Nevertheless, a horizontal branch is not obvious in the CMD. Given the low metallicity of the galaxy, this likely indicates that the first substantial star formation in the galaxy may have been somehow delayed in Leo I in comparison with the other dSph satellites of the Milky Way. The subgiant region is well and uniformly populated from the oldest turnoffs up to the 1 Gyr old turnoff, indicating that star formation has proceeded in a continuous way, with possible variations in intensity but no big gaps between successive bursts, over the galaxy's lifetime. The structure of the red-clump of core He-burning stars is consistent with the large amount of intermediate-age population inferred from the main sequence and the subgiant region. In spite of the lack of gas in Leo I, the CMD clearly shows star formation continuing until 1 Gyr ago and possibly until a few hundred Myrs ago in the central part of the galaxy.

In this study we present a simple model of elliptical galaxies aimed at interpreting the gradients in colours and narrow band indices observed across these systems. Salient features of the model are the gradients in mass density and star formation and infall of primordial gas aimed at simulating the collapse of a galaxy into the potential well of dark matter. Adopting a multi-zone model we follow in detail the history of star formation, gas consumption, and chemical enrichment of the galaxy and also allow for the occurrence of galactic winds according to the classical supernova (and stellar winds) energy deposit. The outline of the model, the time scale of gas accretion and rate of star formation as a function of the galacto-centric distance in particular, seek to closely mimic the results from Tree-SPH dynamical models. Although same specific ingredients of the model can be questioned from many points of view (of which we are well aware) the model predictions have to be considered as a gross tool for exploring the consequences of different receipts of gas accretion and star formation in which the simple one-zone scheme is abandoned. With the aid of this model we discuss the observational data on the gradients in metallicity, colour, and narrow band indices across elliptical galaxies.

Systematic variations in the line strength indices H_beta, Mg2, and <Fe> are observed across elliptical galaxies and from one object to another. Furthermore arguments are given for an enhancement of Mg (\alfa in general) with respect to Fe toward the center of these galaxies. In this study we have investigated the ability of the H_beta, Mg2 and <Fe> diagnostics to assess the metallicity, [Mg/Fe] ratios, and ages of elliptical galaxies. To this aim, first we derive basic calibrations for the variations dH_beta, dMg2 and d<Fe> as a function of Dlog(t), Dlog(Z/Zo), and D[Mg/Fe]. Second, we analyze how the difference dH_beta, dMg2, and d<Fe> between the external and central values translates into D[Mg/Fe], Dlog(Z/Zo), and Dlog(t). Third, we explore the variation from galaxy to galaxy of the nuclear values of H_beta, Mg2, and <Fe> limited to a sub-sample of the Gonzales (1993) list. The differences dH_beta, dMg2, and d<Fe> are converted into the differences Dlog(t), Dlog(Z/Zo), and D[Mg/Fe]. Various correlations among the age, metallicity, and enhancement variations are explored. In particular we thoroughly examine the relationship Dlog(t)-Mv. It is found that a sort of age limit is likely to exist in this plane, traced by galaxies with mild or no sign of rejuvenation. In these objects, the duration of the star forming activity is likely to have increased at decreasing galactic mass. The results are discussed in the context of current theories of galaxy formation and evolution. i.e. merger and isolation. Finally, the suggestion is advanced that models with time and space dependent IMF might be able to alleviate some of the difficulties encountered by the standard SN-driven galactic wind model.

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.

New detailed stellar yields of several elemental species are derived for massive stars in a wide range of masses (from 6 to 120 Msol) and metallicities (Z= 0.0004, 0.004, 0.008, 0.02, 0.05). Our calculations are based on the Padova evolutionary tracks and take into account recent results on stellar evolution, such as overshooting and quiescent mass-loss, paying major attention to the effects of the initial chemical composition of the star. We finally include modern results on explosive nucleosynthesis in SNae by Woosley & Weaver 1995. The issue of the chemical yields of Very Massive Objects (from 120 to 1000 Msol) is also addressed. Our grid of stellar yields for massive stars is complementary to the results by Marigo et al. (1996, 1997) on the evolution and nucleosynthesis of low and intermediate mass stars, also based on the Padova evolutionary tracks. Altogether, they represent a complete set of stellar yields of unprecedented homogeneity and self-consistency. Our new stellar yields are inserted in a code for the chemical evolution of the Galactic disc with infall of primordial gas, according to the formulation first suggested by Talbot & Arnett (1971, 1973, 1975) and Chiosi (1980). As a first application, the code is used to develop a model of the chemical evolution of the Solar Vicinity, with a detailed comparison to the available observational constraints.

In this paper we present new models of massive stars based on recent advancements in the theory of diffusive mixing and a new empirical formulation of the mass-loss rates of red supergiant stars. We compute two sets of stellar models of massive stars with initial chemical composition [Z=0.008, Y=0.25] and [Z=0.020, Y=0.28]. Mass loss by stellar wind is also taken into account according to the empirical relationship by de Jager et al. (1988). Despite the new mixing prescription, these models share the same problems of older models in literature as far as the interpretation of the observational distribution of stars across the HRD is concerned. Examining possible causes of the failure, we find that the adopted rate of mass loss for the red supergiant stages under-estimates the observational values by a large factor. Revising the whole problem, first we adopt the recent formulation by Feast (1991), and secondly we take also into account the possibility that the dust to gas ratio varies with the stellar luminosity. Stellar models are then calculated with the new prescription for the mass-loss rates. The models now possess very extended loops in the HRD and are able to match the distribution of stars across the HRD from the earliest to the latest spectral types both in the Milky Way, LMC and SMC. Finally, because the surface H-abundance is close to the limit adopted to start the Wolf-Rayet phase (WNL type), we suggest that a new channel is possible for the formation of low luminosity WNL stars, i.e. by progenitors whose mass can be as low as 20 \msol, that have evolved horizontally across the HRD following the blue-red-blue scheme and suffering large mass loss during the red supergiant stages.

In this paper we focus on the TP-AGB evolution of intermediate-mass stars experiencing envelope burning (M=4-5Mo). The approach we have adopted is a semi-analytical one as it combines analytical relationships derived from complete models of TP-AGB stars with sole envelope models in which the physical structure is calculated from the photosphere down to the core. Our efforts are directed to analyse the effects produced by envelope burning, such as: i) the energy contribution which may drive significant deviations from the standard core mass-luminosity relationship; and ii) the changes in the surface chemical composition due to nuclear burning via the CNO cycle. Evolutionary models for stars with initial mass of 4.0, 4.5, 5.0 Mo and two choices of the initial chemical composition ([Y=0.28,Z=0.02] and [Y=0.25,Z=0.008]) are calculated from the first thermal pulse till the complete ejection of the envelope. We find that massive TP-AGB stars can rapidly reach high luminosities (-6>Mbol>-7), without exceeding, however, the classical limit to the AGB luminosity of Mbol~-7.1 corresponding to the Chandrasekhar value of the core mass. No carbon stars brighter than Mbol~-6.5 are predicted to form (the alternative of a possible transition from M-star to C-star during the final pulses is also explored), in agreement with observations which indicate that most of the very luminous AGB stars are oxygen-rich. Finally, new chemical yields from stars in the mass range 4-5 Mo are presented, so as to extend the sets of stellar yields from low-mass stars already calculated by Marigo et al. (1996). For each CNO element we give both the secondary and the primary components.

The gradients in line strength indices such as Mg2 and <Fe> (Carollo & Danziger 1994a,b; Carollo et al. 1993), and their increase toward massive system (Worthey et al. 1992, Faber et al. 1992, Davies et al. 1993, Weiss et al. 1995) observed in elliptical galaxies are customarily considered as indicators of gradients in the chemical abundances of Mg and Fe and arguments are given for an enhancement of Mg ($\alpha$-elements in general) with respect to Fe toward the center of these galaxies or going from dwarf to massive ellipticals. In this paper, we present a detailed analysis of the whole problem aimed at understanding how the indices Mg2 and <Fe> depend on the the chemical abundances and age of the stellar mix of a galaxy. We find that Mg2 and <Fe> do not simply correlate with the abundances of Mg and Fe and the ratio Mg/Fe, but depend also on the underlaying distribution of stars in different metallicity bins N(Z), or equivalently on the past history of star formation. In practice, inferring the abundance of Mg and Fe and their ratio is hampered by the need of some information about N(Z). Finally, the observational gradients in Mg2 and <Fe> do not automatically imply gradients in chemical abundances.

In this paper, basic observational properties of elliptical galaxies such as the integrated spectra, the chemical abundances and the enhancement of $\alpha$-elements inferred from broad-band colors and line strength indices $Mg_2$ and $< Fe >$ (and their gradients), the color-magnitude relation, the UV fluxes, and the mass-luminosity ratios, are examined in the light of current theoretical interpretations, and attention is called on several points of internal contradiction. Indeed existing models for the formation and evolution of elliptical galaxies are not able to simultaneously account for all of the above observational features. We suggest that the new initial mass function (IMF) by Padoan et al. (1997), which depends on the temperature, density, and velocity dispersion of the medium in which stars are formed, may constitute a break-through in the difficulties in question. Models of elliptical galaxies incorporating the new IMF (varying with time and position inside a galaxy) are presented and discussed at some extent. As a result of the changing IMF, the enhancement of $\alpha$-elements and tilt of the Fundamental Plane are easily explained leaving unaltered the interpretation of the remaining properties under examination. Finally, some implications concerning the relative proportions of visible stars, collapsed remnants (baryonic dark matter), and gas left over by the star forming process are examined.

We present the first results of a detailed modeling of chemical and photometric evolution of galaxies including the effects of a dusty interstellar medium. A chemical evolution code follows the SF rate, the gas fraction and the metallicity, basic ingredients for the stellar population synthesis. The latter is performed with a grid of integrated spectra for SSP of different ages and metallicities, in which the effects of dusty envelopes around AGB stellar are included. The residual fraction of gas in the galaxy is divided into two phases: the star forming molecular clouds and the diffuse medium. The molecular gas is sub-divided into clouds of given mass and radius: it is supposed that any SSP is borne within the cloud and progressively escapes it. The emitted spectrum of the star forming molecular clouds is computed with a radiative transfer code. The diffuse dust emission (cirrus) is derived by describing the galaxy as an axially symmetric system, in which the local dust emissivity is consistently calculated as a function of the local field intensity due to the stellar component. Effects of very small grains, subject to temperature fluctuations, as well as PAH are included. The model is compared and calibrated with available data of starburst galaxies in the local universe, in particular new broad-band and spectroscopic ISO observations. It will be a fundamental tool to extract information on global evolution of the stellar component and of ISM of galaxies from observations covering four decades in $\lambda$. The combination of observations from HST, Keck, ISO and ground-based optical, IR and sub-mm telescopes already provided several objects at relevant $z$ with spectral information on this large $\lambda$ range. Their number will increase after the completion of ISO surveys, and will burst when SIRFT, FIRST, Planck

The line strength indices Hb, Mg2, and <Fe>, together with their gradients and dependence on the galaxy luminosity (mass), observed in elliptical galaxies are customarily considered as reliable indicators of systematic differences in age and abundances of Mg and Fe. In this paper, we address the question whether or not the indices Mg2 and <Fe> (and their gradients) are real indicators of chemical abundances and enhancement of these or other effects must be taken into account. We show that this is not the case because the above indices are severely affected by the unknown relative percentage of stars as function of the metallicity. In order to cast the problem, first we provide basic calibrations for the variations dHb, dMg2, and d<Fe> as a function of the age Dlog(t) (in Gyr), metallicity Dlog(Z/Z_o), and D[Mg/Fe]. Second, we analyze the implications of the gradients in Mg2 and <Fe> observed across these systems. Finally, the above calibration is used to explore the variations from galaxy to galaxy of the nuclear values of Hb, Mg2, and <Fe>. The differences dHb, dMg2, and d<Fe> are converted to differences Dlog(t), Dlog(Z/Z_o), and D[Mg/Fe]. We thoroughly examine the relationship Dlog(t)-Mv, Dlog(Z/Z_o)-Mv, and D[Mg/Fe]-Mv, and advance the suggestion that the duration of the star forming period gets longer or the age of the last episode of stellar activity gets closer to the present at decreasing galaxy mass. This result is discussed in the context of current theories of galaxy formation and evolution. i.e. merger and isolation. In brief, we conclude that none of these can explain the results of our analysis, and suggest that the kind of time and space dependent IMF proposed by Padoan et al. (1997) and the associated models of elliptical galaxies elaborated by Chiosi et al. (1997) should be at work.

In this paper, we discuss a method to conduct a quantitative study of the star formation history (SFH) of Local Group (LG) galaxies using (HST) data. This method has proven to be successful in the analysis of the SFH of the same kind of galaxies using ground-based observations. It is based on the comparison of observed CMDs with a set of model CMDs. The latter are computed assuming different evolutionary scenarios, and include a detailed simulation of observational effects. HST CMDs are ~3 mags deeper than typical ground-based CMDs, allowing the observation, for all LG galaxies, of a part of the CMD that up till now had remained accessible only for the very nearest galaxies. A very important feature that will become accessible is the HB+red-clump. The distribution of stars along this structure is quite sensitive to age and metallicity and should provide a very important improvement in the time resolution of the SFH for stars older than ~2-3 Gyr. We show and discuss four model CMDs which would be comparable with CMDs from deep HST observations. These model CMDs represent the following evolutionary scenarios corresponding to a wide range of dwarf galaxy sub-types from dI to dE: A) a constant SFR from 15Gyr ago to the present time; B) as A), but with the SFR stopped 0.5 Gyr ago; C) a constant SFR in the age range 10-9Gyr and D) as C) but in the age range 15-12 Gyr. In all four cases a range of metallicity from Z=0.0001 to Z=0.004 has been assumed. The present analysis is just a first qualitative approach to what one may expect to find in the CMDs of LG galaxies. However a complete set of model CMDs must be computed to analize the data for each galaxy, using the crowding effects derived for that particular galaxy.

Observations of (B-V) colour distributions towards the galactic poles are compared with those obtained from synthetic colour-magnitude diagrams to determine the major constituents in the disc and spheroid. The disc is described with four stellar sub-populations: the young, intermediate, old, and thick disc populations, which have respectively scale heights of 100 pc, 250 pc, 0.5 kpc, and 1.0 kpc. The spheroid is described with stellar contributions from the bulge and halo. The bulge is not well constrained with the data analyzed in this study. A non-flattened power-law describes the observed distributions at fainter magnitudes better than a deprojected R^{1/4}-law. Details about the age, metallicity, and normalizations are listed in Table 1. The star counts and the colour distributions from the stars in the intermediate fields towards the galactic anti-centre are well described with the stellar populations mentioned above. Arguments are given that the actual solar offset is about 15 pc north from the galactic plane.

In this paper we address the question whether age and metallicity effects can be disentangled with the aid of the broad-band colours and spectral indices from absorption feature strengths, so that the age of elliptical galaxies can be inferred. The observational data under examination are the indices Hbeta and [MgFe], and the velocity dispersion Sigma for the sample of galaxies of Gonzales (1993), supplemented by the ultra-violet data, i.e. the colour (1550-V), of Burstein et al. (1988). The analysis is performed with the aid of chemo-spectro-photometric models of elliptical galaxies with infall of primordial gas (aimed at simulating the collapse phase of galaxy formation) and the occurrence of galactic winds. The galaxy models are from Tantalo et al. (1995).

In this paper we present new chemo-spectro-photometric models of elliptical galaxies in which infall of primordial gas is allowed to occur. They aim to simulate the collapse of a galaxy made of two components, i.e. luminous material and dark matter. The mass of the dark component is assumed to be constant in time, whereas that of the luminous material is supposed to accrete at a suitable rate. They also include the effect of galactic winds powered by supernova explosions and stellar winds from massive, early-type stars. The models are constrained to match a number of properties of elliptical galaxies, i.e. the slope and mean colours of the colour-magnitude relation (CMR), V versus (V--K), the UV excess as measured by the colour (1550--V) together with the overall shape of the integrated spectral energy distribution (ISED) in the ultraviolet, the relation between the Mg2 index and (1550--V), the mass to blue luminosity ratio M/Lb as a function of the B luminosity, and finally the broad-band colours (U--B), (B--V), (V--I), (V--K), etc.

The stellar content of the Im galaxy NGC 2366 is discussed on the basis of CCD BVR photometry. The three brightest blue and red stars have been used to estimate its distance, obtaining a balue of 2.9 Mpc. The spatial distribution of the young stellar population is discussed in the light of the integrated color indices and the color-magnitude diagrams of different zones of the galaxy. A generalized star formation burst seems to have taken place about 50 Myr ago. The youngest stars are preferentially formed in the South-West part of the bar, where the giant HII complex NGC 2363 is located, being younger and bluer. The bar seems to play a role favouring star formation in one of its extremes. Self-propagation however, does not seem to be triggering star formation at large scale. A small region, populated by very young stars has also been found at the East of the galaxy.