The Milky Way's progenitor experienced several merger events which left their imprints on the stellar halo, including the Gaia-Sausage/Enceladus. Recently, it has been proposed that this event perturbed the proto-disk and gave rise to a metal rich ([Fe/H] $>-1$) low angular momentum ($v_{\phi} < 100$ km/s) stellar population. These stars have dynamical and chemical properties different from the accreted stellar halo, but are continuous with the canonical thick disk. In this letter, we use a hydrodynamical simulation of an isolated galaxy which develops clumps that produce a bimodal thin$+$thick disk chemistry to explore whether it forms such a population. We demonstrate clump scattering forms a metal-rich, low angular momentum population, without the need for a major merger. We show that, in the simulation, these stars have chemistry, kinematics and density distribution in good agreement with those in the Milky Way.

A scenario for the formation of the bi-modality in the chemical space [$\alpha$/Fe] vs [Fe/H] of the Milky Way was recently proposed in which $\alpha$-enhanced stars are produced early and quickly in clumps. Besides accelerating the enrichment of the medium with $\alpha$-elements, these clumps scatter the old stars, converting in-plane to vertical motion, forming a geometric thick disc. In this paper, by means of a detailed analysis of the data from smooth particle hydrodynamical simulations, we investigate the geometric properties (in particular of the chemical thick disc) produced in this scenario. For mono-age populations we show that the surface radial density profiles of high-[$\alpha$/Fe] stars are well described by single exponentials, while that of low-[$\alpha$/Fe] stars require broken exponentials. This break is sharp for young populations and broadens for older ones. The position of the break does not depend significantly on age. The vertical density profiles of mono-age populations are characterized by single exponentials, which flare significantly for low-[$\alpha$/Fe] stars but only weakly (or not at all) for high-[$\alpha$/Fe] stars. For low-[$\alpha$/Fe] stars, the flaring level decreases with age, while for high-[$\alpha$/Fe] stars it weakly increases with age (although with large uncertainties). All these properties are in agreement with observational results recently reported for the Milky Way, making this a plausible scenario for the formation of the Galactic thick disc.

We introduce the study of box/peanut (B/P) bulges in the action space of the initial axisymmetric system. We explore where populations with different actions end up once a bar forms and a B/P bulge develops. We find that the density bimodality due to the B/P bulge (the X-shape) is better traced by populations with low radial, JR,0, or vertical, Jz,0, actions, or high azimuthal action, J{\phi},0. Generally populations separated by JR,0 have a greater variation in bar strength and vertical heating than those separated by Jz,0. While the bar substantially weakens the initial vertical gradient of Jz,0, it also drives a strikingly monotonic vertical profile of JR,0. We then use these results to guide us in assigning metallicity to star particles in a pure N-body model. Because stellar metallicity in unbarred galaxies depends on age as well as radial and vertical positions, the initial actions are particularly well suited for assigning metallicities. We argue that assigning metallicities based on single actions, or on positions, results in metallicity distributions inconsistent with those observed in real galaxies. We therefore use all three actions to assign metallicity to an N-body model by comparing with the actions of a star-forming, unbarred simulation. The resulting metallicity distribution is pinched on the vertical axis, has a realistic vertical gradient and has a stronger X-shape in metal-rich populations, as found in real galaxies.

This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

We show that the Weibel or currente filamentation instability can lead to the emission of circularly polarized radiation. Using particle-in-cell (PIC) simulations and a radiation post-processing numerical algorithm, we demonstrate that the level of circular polarization increases with the initial plasma magnetization, saturating at ~13% when the magnetization, given by the ratio of magnetic energy density to the electron kinetic energy density, is larger than 0.05. Furthermore, we show that this effect requires an ion-electron mass ratio greater than unity. These findings, which could also be tested in currently available laboratory conditions, show that the recent observation of circular polarization in gamma ray burst afterglows could be attributed to the presence of magnetized current filaments driven by the Weibel or the current filamentation instability.

We investigate the old problem of the fast relaxation of collisionless $N$-body systems which are collapsing or perturbed, emphasizing the importance of (non-collisional) discreteness effects. We integrate orbit ensembles in fixed external potentials, estimating the entropy of the ensemble to analyze the time evolution of the distribution function. We show that these estimates capture the correct physical behavior expected from the 2nd Law of Thermodynamics, without any spurious entropy production. For self-consistent (i.e. stationary) samples, the entropy is conserved, while for non-self-consistent samples, it increases within a few dynamical times up to a maximum where it stabilizes (even in integrable potentials). Our results shed light on the main ingredients for this fast collisionless relaxation. The fundamental ingredient is the discreteness (finite $N$) of gravitational systems in any potential. Additionally, in non-integrable potentials, the presence of chaotic orbits accelerates the entropy production. Contrary to the traditional violent relaxation scenario, our results indicate that a time-dependent potential is not necessary for this fast relaxation. For the first time, in connection with the Nyquist-Shannon theorem we determine the $N$-dependence of this discreteness-driven relaxation, deriving a typical relaxation time $T/\tau_{cr}\approx 0.1 N^{1/6}$, with slightly weaker $N$-dependencies for non-integrable potentials with substantial fractions of chaotic orbits. This timescale is much smaller than the collisional timescale even for small-$N$ systems such as open clusters and represents an upper limit for real collisionless $N$-body systems. Additionally, our results reinforce the conclusion draw in Beraldo e Silva et al. (2017) that the Vlasov equation does not provide an adequate kinetic description of the fast relaxation of collapsing collisionless $N$-body systems.

The time irreversibility and fast relaxation of collapsing $N$-body gravitating systems (as opposed to the time reversibility of the equations of motion for individual stars or particles) are traditionally attributed to information loss due to coarse-graining in the observation. We show that this subjective element is not necessary once one takes into consideration the fundamental fact that these systems are discrete, i.e. composed of a finite number $N$ of stars or particles. We show that a connection can be made between entropy estimates for discrete systems and the Nyquist-Shannon sampling criterion. Specifically, given a sample with $N$ points in a space of $d$ dimensions, the Nyquist-Shannon criterion constrains the size of the smallest structures defined by a function in the continuum that can be uniquely associated with the discrete sample. When applied to an $N$-body system, this theorem sets a lower limit to the size of phase-space structures (in the continuum) that can be resolved in the discrete data. As a consequence, the finite $N$ system tends to a uniform distribution after a relaxation time that typically scales as $N^{1/d}$. This provides an explanation for the fast achievement of a stationary state in collapsing $N$-body gravitating systems such as galaxies and star clusters, without the need to advocate for the subjective effect of coarse-graining.

We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple "recipes" for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime, in particular between 25 and 70 micron where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard "recipes" for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a "standard" constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 micron) luminosities and colors. From a Principal Component Analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within roughly 0.3 dex using only Mstar and SFR.

Strong magnetic fields in magnetospheres of neutron stars (especially magnetars) and other astrophysical objects may release their energy in violent, intense episodes of magnetic reconnection. While reconnection has been studied extensively, the extreme field strength near neutron stars introduces new effects: synchrotron cooling and electron-positron pair production. Using massively parallel particle-in-cell simulations that self-consistently incorporate these new quantum-electrodynamic effects, we investigate relativistic magnetic reconnection in the strong-field regime. We show that reconnection in this regime can efficiently convert magnetic energy to X-ray and gamma-ray radiation and thus power bright high-energy astrophysical flares. Rapid radiative cooling causes strong plasma and magnetic field compression in compact plasmoids. In the most extreme cases, the field can approach the critical quantum limit, leading to copious pair production.

Ab initio simulations of the propagation in a plasma of a soon to be available relativistic electron-positron beam or fireball beam provide an effective mean for the study of microphysics relevant to astrophysical scenarios. We show that the current filamentation instability associated with some of these scenarios reaches saturation after only 10 cm of propagation in a typical laboratory plasma with a density 10^17/cc. The different regimes of the instability, from the purely transverse to the mixed mode filamentation, can be accessed by varying the background plasma density. The instability generates large local plasma gradients, intense transverse magnetic fields, and enhanced emission of radiation. We suggest that these effects may be observed experimentally for the first time.

We present two-dimensional particle-in-cell (PIC) simulations of the fully kinetic collisionless magnetorotational instability (MRI) in weakly magnetized (high $\beta$) pair plasma. The central result of this numerical analysis is the emergence of a self-induced turbulent regime in the saturation state of the collisionless MRI, which can only be captured for large enough simulation domains. One of the underlying mechanisms for the development of this turbulent state is the drift-kink instability (DKI) of the current sheets resulting from the nonlinear evolution of the channel modes. The onset of the DKI can only be observed for simulation domain sizes exceeding several linear MRI wavelengths. The DKI, together with ensuing magnetic reconnection, activate the turbulent motion of the plasma in the late stage of the nonlinear evolution of the MRI. At steady state, the magnetic energy has an MHD-like spectrum with a slope of $k^{-5/3}$ for $k\rho<1$ and $k^{-3}$ for sub-Larmor scale ($k\rho>1$). We also examine the role of the collisionless MRI and associated magnetic reconnection in the development of pressure anisotropy. We study the stability of the system due to this pressure anisotropy, observing the development of mirror instability during the early-stage of the MRI. We further discuss the importance of magnetic reconnection for particle acceleration during the turbulence regime. In particular, consistent with reconnection studies, we show that at late times the kinetic energy presents a characteristic slope of $\epsilon^{-2}$ in the high-energy region.

The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed.

The collapse of a collisionless self-gravitating system, with the fast achievement of a quasi-stationary state, is driven by violent relaxation, with a typical particle interacting with the time-changing collective potential. It is traditionally assumed that this evolution is governed by the Vlasov-Poisson equation, in which case entropy must be conserved. We run N-body simulations of isolated self-gravitating systems, using three simulation codes: NBODY-6 (direct summation without softening), NBODY-2 (direct summation with softening) and GADGET-2 (tree code with softening), for different numbers of particles and initial conditions. At each snapshot, we estimate the Shannon entropy of the distribution function with three different techniques: Kernel, Nearest Neighbor and EnBiD. For all simulation codes and estimators, the entropy evolution converges to the same limit as N increases. During violent relaxation, the entropy has a fast increase followed by damping oscillations, indicating that violent relaxation must be described by a kinetic equation other than the Vlasov-Poisson, even for N as large as that of astronomical structures. This indicates that violent relaxation cannot be described by a time-reversible equation, shedding some light on the so-called "fundamental paradox of stellar dynamics". The long-term evolution is well described by the orbit-averaged Fokker-Planck model, with Coulomb logarithm values in the expected range 10-12. By means of NBODY-2, we also study the dependence of the 2-body relaxation time-scale on the softening length. The approach presented in the current work can potentially provide a general method for testing any kinetic equation intended to describe the macroscopic evolution of N-body systems.

In this work we investigate the far-UV to NIR shape of the dust attenuation curve of a sample of IR selected dust obscured (U)LIRGs at z$\sim$2. The spectral energy distributions (SEDs) are fitted with CIGALE, a physically-motivated spectral synthesis model based on energy balance. Its flexibility allows us to test a wide range of different analytical prescriptions for the dust attenuation curve, including the well-known Calzetti and Charlot & Fall curves, and modified versions of them. The attenuation curves computed under the assumption of our reference double power-law model are in very good agreement with those derived, in previous works, with radiative transfer (RT) SED fitting. We investigate the position of our galaxies in the IRX-$\beta$ diagram and find this to be consistent with grayer slopes, on average, in the UV. We also find evidence for a flattening of the attenuation curve in the NIR with respect to more classical Calzetti-like recipes. This larger NIR attenuation yields larger derived stellar masses from SED fitting, by a median factor of $\sim$ 1.4 and up to a factor $\sim$10 for the most extreme cases. The star formation rate appears instead to be more dependent on the total amount of attenuation in the galaxy. Our analysis highlights the need for a flexible attenuation curve when reproducing the physical properties of a large variety of objects.

Kepler-452b is currently the best example of an Earth-size planet in the habitable zone of a sun-like star, a type of planet whose number of detections is expected to increase in the future. Searching for biosignatures in the supposedly thin atmospheres of these planets is a challenging goal that requires a careful selection of the targets. Under the assumption of a rocky-dominated nature for Kepler-452b, we considered it as a test case to calculate a temperature-dependent habitability index, $h_{050}$, designed to maximize the potential presence of biosignature-producing activity (Silva et al.\ 2016). The surface temperature has been computed for a broad range of climate factors using a climate model designed for terrestrial-type exoplanets (Vladilo et al.\ 2015). After fixing the planetary data according to the experimental results (Jenkins et al.\ 2015), we changed the surface gravity, CO$_2$ abundance, surface pressure, orbital eccentricity, rotation period, axis obliquity and ocean fraction within the range of validity of our model. For most choices of parameters we find habitable solutions with $h_{050}>0.2$ only for CO$_2$ partial pressure $p_\mathrm{CO_2} \lesssim 0.04$\,bar. At this limiting value of CO$_2$ abundance the planet is still habitable if the total pressure is $p \lesssim 2$\,bar. In all cases the habitability drops for eccentricity $e \gtrsim 0.3$. Changes of rotation period and obliquity affect the habitability through their impact on the equator-pole temperature difference rather than on the mean global temperature. We calculated the variation of $h_{050}$ resulting from the luminosity evolution of the host star for a wide range of input parameters. Only a small combination of parameters yield habitability-weighted lifetimes $\gtrsim 2$\,Gyr, sufficiently long to develop atmospheric biosignatures still detectable at the present time.

By means of the updated PARSEC database of evolutionary tracks of massive stars, we compute the integrated stellar light, the ionizing photon budget and the supernova rates of young simple stellar populations (SSPs), for different metallicities and IMF upper mass limits. Using CLOUDY we compute and include in the SSP spectra the neb- ular emission contribution. We also revisit the thermal and non-thermal radio emission contribution from young stars. Using GRASIL we can thus predict the panchromatic spectrum and the main recombination lines of any type of star-forming galaxy, including the effects of dust absorption and re-emission. We check the new models against the spectral energy distributions (SEDs) of selected well-observed nearby galaxies. From the best-fit models we obtain a consistent set of star formation rate (SFR) calibrations at wavelengths ranging from ultraviolet (UV) to radio. We also provide analytical calibrations that take into account the dependence on metallcity and IMF upper mass limit of the SSPs. We show that the latter limit can be well constrained by combining information from the observed far infrared, 24 {\mu}m, 33 GHz and H{\alpha} luminosities. Another interesting property derived from the fits is that, while in a normal galaxy the attenuation in the lines is significantly higher than that in the nearby continuum, in individual star bursting regions they are similar, supporting the notion that this effect is due to an age selective extinction. Since in these conditions the Balmer decrement method may not be accurate, we provide relations to estimate the attenuation from the observed 24 {\mu}m or 33 GHz fluxes. These relations can be useful for the analysis of young high redshift galaxies.

The velocity distribution function (VDF) of dark matter (DM) halos in $\Lambda$CDM dissipationless cosmological simulations, which must be non-separable in its radial and tangential components, is still poorly known. We present the first single-parameter, non-separable, anisotropic model for the VDF in $\Lambda$CDM halos, built from an isotropic $q$-Gaussian (Tsallis) VDF of the isotropic set of dimensionless spherical velocity components (after subtraction of streaming motions), normalized by the respective velocity dispersions. We test our VDF on 90 cluster-mass halos of a dissipationless cosmological simulation. Beyond the virial radius, $r_{\rm vir}$, our model VDF adequately reproduces that measured in the simulated halos, but no $q$-Gaussian model can adequately represent the VDF within $r_{\rm vir}$, as the speed distribution function is then flatter-topped than any $q$-Gaussian can allow. Nevertheless, our VDF fits significantly better the simulations than the commonly used Maxwellian (Gaussian) distribution, at virtually all radii within $5\,r_{\rm vir}$. Within 0.4 (1) $r_{\rm vir}$, the non-Gaussianity index $q$ is (roughly) linearly related to the slope of the density profile and also to the velocity anisotropy profile. We provide a parametrization of the modulation of $q$ with radius for both the median fits and the fit of the stacked halo. At radii of a few percent of $r_{\rm vir}$, corresponding to the Solar position in the Milky Way, our best-fit VDF, although fitting better the simulations than the Gaussian one, overproduces significantly the fraction of high velocity objects, indicating that one should not blindly use these $q$-Gaussian fits to make predictions on the direct detection rate of DM particles.

A generalized equation of state corresponding to a model that includes a Chaplygin gas and a viscous term is investigated, in the context of the reconstruction program in scalar field cosmology. The corresponding inflationary model parameters can be conveniently adjusted in order to reproduce the most recent PLANCK data. The influence of the Chaplygin gas term contribution, in relation with previous models, is discussed. Exit from inflation is shown to occur quite naturally in the new model.

In an effort to derive temperature based criteria of habitability for multicellular life, we investigated the thermal limits of terrestrial poikilotherms, i.e. organisms whose body temperature and the functioning of all vital processes is directly affected by the ambient temperature. Multicellular poikilotherms are the most common and evolutionarily ancient form of complex life on earth. The thermal limits for their active metabolism and reproduction are bracketed by the temperature interval 0C<T<50C. The same interval applies to the photosynthetic production of oxygen, an essential ingredient of complex life, and for the generation of atmospheric biosignatures. Analysis of the main mechanisms responsible for the thermal thresholds of terrestrial life suggests that the same mechanisms would apply to other forms of chemical life. We propose a habitability index for complex life, h050, representing the mean orbital fraction of planetary surface that satisfies the temperature limits 0C<T<50C. With the aid of a climate model tailored for the calculation of the surface temperature of Earth-like planets, we calculated h050 as a function of planet insolation S, and atmospheric columnar mass Natm, for a few earth-like atmospheric compositions. By displaying h050 as a function of S and Natm, we built up an atmospheric mass habitable zone (AMHZ) for complex life. At variance with the classic habitable zone, the inner edge of the complex life HZ is not affected by the uncertainties inherent to the calculation of the runaway greenhouse limit. The complex life HZ is significantly narrower than the HZ of dry planets. Our calculations illustrate how changes in ambient conditions dependent on S and Natm, such as temperature excursions and surface dose of secondary particles of cosmic rays, may influence the type of life potentially present at different epochs of planetary evolution inside the AMHZ.

Collisionless shocks occur in various fields of physics. In the context of space and astrophysics they have been investigated for many decades. However, a thorough understanding of shock formation and particle acceleration is still missing. Collisionless shocks can be distinguished into electromagnetic and electrostatic shocks. Electromagnetic shocks are of importance mainly in astrophysical environments and they are mediated by the Weibel or filamentation instability. In such shocks, charged particles gain energy by diffusive shock acceleration. Electrostatic shocks are characterized by a strong electrostatic field, which leads to electron trapping. Ions are accelerated by reflection from the electrostatic potential. Shock formation and particle acceleration will be discussed in theory and simulations.

The microphysics of relativistic collisionless sheared flows is investigated in a configuration consisting of a globally neutral, relativistic $e^-e^+$ beam streaming through a hollow plasma/dielectric channel. We show through multidimensional PIC simulations that this scenario excites the Mushroom instability (MI), a transverse shear instability on the electron-scale, when there is no overlap (no contact) between the $e^-e^+$ beam and the walls of the hollow plasma channel. The onset of the MI leads to the conversion of the beam's kinetic energy into magnetic (and electric) field energy, effectively slowing down a globally neutral body in the absence of contact. The collisionless shear physics explored in this configuration may operate in astrophysical environments, particularly in highly relativistic and supersonic settings where macroscopic shear processes are stable.

In this paper we show a plausible mechanism that could lead to the formation of the Dark Lanes in Lunar Swirls, and the electromagnetic shielding of the lunar surface that results in the preservation of the white colour of the lunar regolith. We present the results of a fully self-consistent 2 and 3 dimensional particle-in-cell simulations of mini-magnetospheres that form above the lunar surface and show that they are consistent with the formation of `lunar swirls' such as the archetypal formation Reiner Gamma. The simulations show how the microphysics of the deflection/shielding of plasma operates from a kinetic-scale cavity, and show that this interaction leads to a footprint with sharp features that could be the mechanism behind the generation of `dark lanes'. The physics of mini-magnetospheres is described and shown to be controlled by space-charge fields arising due to the magnetized electrons and unmagnetized ions. A comparison between model and observation is shown for a number of key plasma parameters.

Electron-scale surface waves are shown to be unstable in the transverse plane of a shear flow in an initially unmagnetized plasma, unlike in the (magneto)hydrodynamics case. It is found that these unstable modes have a higher growth rate than the closely related electron-scale Kelvin-Helmholtz instability in relativistic shears. Multidimensional particle-in-cell simulations verify the analytic results and further reveal the emergence of mushroom-like electron density structures in the nonlinear phase of the instability, similar to those observed in the Rayleigh Taylor instability despite the great disparity in scales and different underlying physics. Macroscopic ($\gg c/\omega_{pe}$) fields are shown to be generated by these microscopic shear instabilities, which are relevant for particle acceleration, radiation emission and to seed MHD processes at long time-scales.

We introduce a novel Earth-like planet surface temperature model (ESTM) for habitability studies based on the spatial-temporal distribution of planetary surface temperatures. The ESTM adopts a surface Energy Balance Model complemented by: radiative-convective atmospheric column calculations, a set of physically-based parameterizations of meridional transport, and descriptions of surface and cloud properties more refined than in standard EBMs. The parameterization is valid for rotating terrestrial planets with shallow atmospheres and moderate values of axis obliquity (epsilon >= 45^o). Comparison with a 3D model of atmospheric dynamics from the literature shows that the equator-to-pole temperature differences predicted by the two models agree within ~5K when the rotation rate, insolation, surface pressure and planet radius are varied in the intervals 0.5 <= Omega/Omega_o <= 2, 0.75 <= S/S_o <= 1.25, 0.3 <= p/(1 bar) <= 10, and 0.5 <= R/R_o <= 2, respectively. The ESTM has an extremely low computational cost and can be used when the planetary parameters are scarcely known (as for most exoplanets) and/or whenever many runs for different parameter configurations are needed. Model simulations of a test-case exoplanet (Kepler-62e) indicate that an uncertainty in surface pressure within the range expected for terrestrial planets may impact the mean temperature by ~60 K. Within the limits of validity of the ESTM, the impact of surface pressure is larger than that predicted by uncertainties in rotation rate, axis obliquity, and ocean fractions. We discuss the possibility of performing a statistical ranking of planetary habitability taking advantage of the flexibility of the ESTM.

We analyse the full shock formation process in electron-ion plasmas in theory and simulations. It is accepted that electromagnetic shocks in initially unmagnetised relativistic plasmas are triggered by the filamentation instability. However, the transition from the first unstable phase to the quasi-steady shock is still missing. We derive a theoretical model for the shock formation time, taking into account the filament merging in the non-linear phase of the filamentation instability. This process is much slower than in electron-positron pair shocks, so that the shock formation is longer by a factor proportional to sqrt(m_i/m_e) ln(m_i/m_e).

We present a backward approach for the interpretation of the evolution of the near-infrared and the far-infrared luminosity functions across the redshift range 0<z<3. In our method, late-type galaxies are treated by means of a parametric phenomenological method based on PEP/HerMES data up to z~4, whereas spheroids are described by means of a physically motivated backward model. The spectral evolution of spheroids is modelled by means of a single-mass model, associated to a present-day elliptical with K-band luminosity comparable to the break of the local early-type luminosity function. The formation of proto-spheroids is assumed to occurr across the redshift range 1< z < 5. The key parameter is represented by the redshift z_0.5 at which half proto-spheroids are already formed. A statistical study indicates for this parameter values between z_0.5=1.5 and z_0.5=3. We assume as fiducial value z_0.5~2, and show that this assumption allows us to describe accourately the redshift distributions and the source counts. By assuming z_0.5 ~ 2 at the far-IR flux limit of the PEP-COSMOS survey, the PEP-selected sources observed at z>2 can be explained as progenitors of local spheroids caught during their formation. We also test the effects of mass downsizing by dividing the spheroids into three populations of different present-day stellar masses. The results obtained in this case confirm the validity of our approach, i.e. that the bulk of proto-spheroids can be modelled by means of a single model which describes the evolution of galaxies at the break of the present-day early type K-band LF.

We complement our previous analysis of a sample of z~1-2 luminous and ultra-luminous infrared galaxies ((U)LIRGs), by adding deep VLA radio observations at 1.4 GHz to a large data-set from the far-UV to the sub-mm, including Spitzer and Herschel data. Given the relatively small number of (U)LIRGs in our sample with high S/N radio data, and to extend our study to a different family of galaxies, we also include 6 well sampled near IR-selected BzK galaxies at z~1.5. From our analysis based on the radiative transfer spectral synthesis code GRASIL, we find that, while the IR luminosity may be a biased tracer of the star formation rate (SFR) depending on the age of stars dominating the dust heating, the inclusion of the radio flux offers significantly tighter constraints on SFR. Our predicted SFRs are in good agreement with the estimates based on rest-frame radio luminosity and the Bell (2003) calibration. The extensive spectro-photometric coverage of our sample allows us to set important constraints on the SF history of individual objects. For essentially all galaxies we find evidence for a rather continuous SFR and a peak epoch of SF preceding that of the observation by a few Gyrs. This seems to correspond to a formation redshift of z~5-6. We finally show that our physical analysis may affect the interpretation of the SFR-M* diagram, by possibly shifting, with respect to previous works, the position of the most dust obscured objects to higher M* and lower SFRs.

We propose an association between the phase-space mixing level of a self-gravitating system and the indistinguishability of its constituents (stars or dark matter particles). This represents a refinement in the study of systems exhibiting incomplete violent relaxation. Within a combinatorial analysis similar to that of Lynden-Bell, we make use of this association to obtain a distribution function that deviates from the Maxwell-Boltzmann distribution, increasing its slope for high energies. Considering the smallness of the occupation numbers for large distances from the center of the system, we apply a correction to Stirling's approximation which increases the distribution slope also for low energies. The distribution function thus obtained presents some resemblance to the "S" shape of distributions associated with cuspy density profiles (as compared to the distribution function obtained from the Einasto profile), although it is not quite able to produce sharp cusps. We also argue how the association between mixing level and indistinguishability can provide a physical meaning to the assumption of particle-permutation symmetry in the N-particle distribution function, when it is used to derive the one-particle Vlasov equation, which raises doubts about the validity of this equation during violent relaxation.

A new magnetic field generation mechanism in electrostatic shocks is found, which can produce fields with magnetic energy density as high as 0.01 of the kinetic energy density of the flows on time scales $ \tilde \, 10^4 \, {\omega}_{pe}^{-1}$. Electron trapping during the shock formation process creates a strong temperature anisotropy in the distribution function, giving rise to the pure Weibel instability. The generated magnetic field is well-confined to the downstream region of the electrostatic shock. The shock formation process is not modified and the features of the shock front responsible for ion acceleration, which are currently probed in laser-plasma laboratory experiments, are maintained. However, such a strong magnetic field determines the particle trajectories downstream and has the potential to modify the signatures of the collisionless shock.

We preform hybrid simulations of super Alfvenic quasi-parallel shock, driven by a Coronal Mass Ejection (CME), propagating in the Outer Coronal or Solar Wind at distances of between 3 to 6 solar radii. The hybrid treatment of the problem enable the study of the shock propagation on the ion time scale, preserving ion kinetics and allowing for a self consistent treatment of the shock propagation and particle acceleration. The CME plasma drags the embedded magnetic field lines stretching from the sun, and propagates out into interplanetary space at a greater velocity than the in-situ solar wind, driving the shock, and producing very energetic particles. Our results show electromagnetic Alfven waves are generated at the shock front. The waves propagate upstream of the shock and are produced by the counter streaming ions of the solar wind plasma being reflected at the shock. A significant fraction of the particles are accelerated in two distinct phases: first, particles drift from the shock and are accelerated in the upstream region and, second, particles arriving at the shock get trapped, and are accelerated at the shock front. A fraction of the particles diffused back to the shock, which is consistent with the Fermi acceleration mechanism.

Collisionless shocks in plasmas play an important role in space physics (Earth's bow shock) and astrophysics (supernova remnants, relativistic jets, gamma-ray bursts, high energy cosmic rays). While the formation of a fluid shock through the steepening of a large amplitude sound wave has been understood for long, there is currently no detailed picture of the mechanism responsible for the formation of a collisionless shock. We unravel the physical mechanism at work and show that an electromagnetic Weibel shock always forms when two relativistic collisionless, initially unmagnetized, plasma shells encounter. The predicted shock formation time is in good agreement with 2D and 3D particle-in-cell simulations of counterstreaming pair plasmas. By predicting the shock formation time, experimental setups aiming at producing such shocks can be optimised to favourable conditions.

In this paper we explore the effectiveness of an artificial mini-magnetosphere as a potential radiation shelter for long term human space missions. Our study includes the differences that the plasma environment makes to the efficiency of the shielding from the high energy charged particle component of solar and cosmic rays, which radically alters the power requirements. The incoming electrostatic charges are shielded by fields supported by the self captured environmental plasma of the solar wind, potentially augmented with additional density. The artificial magnetic field generated on board acts as the means of confinement and control. Evidence for similar behaviour of electromagnetic fields and ionised particles in interplanetary space can be gained by the example of the enhanced shielding effectiveness of naturally occurring "mini-magnetospheres" on the moon. The shielding effect of surface magnetic fields of the order of ~100s nanoTesla is sufficient to provide effective shielding from solar proton bombardment that culminate in visible discolouration of the lunar regolith known as "lunar swirls". Supporting evidence comes from theory, laboratory experiments and computer simulations that have been obtained on this topic. The result of this work is, hopefully, to provide the tools for a more realistic estimation of the resources versus effectiveness and risk that spacecraft engineers need to work with in designing radiation protection for long-duration human space missions.

Strong shear flow regions found in astrophysical jets are shown to be important dissipation regions, where the shear flow kinetic energy is converted into electric and magnetic field energy via shear instabilities. The emergence of these self-consistent fields make shear flows significant sites for radiation emission and particle acceleration. We focus on electron-scale instabilities, namely the collisionless, unmagnetized Kelvin-Helmholtz instability (KHI) and a large-scale dc magnetic field generation mechanism on the electron scales. We show that these processes are important candidates to generate magnetic fields in the presence of strong velocity shears, which may naturally originate in energetic matter outburst of active galactic nuclei and gamma-ray bursters. We show that the KHI is robust to density jumps between shearing flows, thus operating in various scenarios with different density contrasts. Multidimensional particle-in-cell (PIC) simulations of the KHI, performed with OSIRIS, reveal the emergence of a strong and large-scale dc magnetic field component, which is not captured by the standard linear fluid theory. This dc component arises from kinetic effects associated with the thermal expansion of electrons of one flow into the other across the shear layer, whilst ions remain unperturbed due to their inertia. The electron expansion forms dc current sheets, which induce a dc magnetic field. Our results indicate that most of the electromagnetic energy developed in the KHI is stored in the dc component, reaching values of equipartition on the order of $10^{-3}$ in the electron time-scale, and persists longer than the proton time-scale. Particle scattering/acceleration in the self generated fields of these shear flow instabilities is also analyzed.

Using cosmological galaxy simulations from the MaGICC project, we study the evolution of the stellar masses, star formation rates and gas phase abundances of star forming galaxies. We derive the stellar masses and star formation rates using observational relations based on spectral energy distributions by applying the new radiative transfer code GRASIL-3D to our simulated galaxies. The simulations match well the evolution of the stellar mass-halo mass relation, have a star forming main sequence that maintains a constant slope out to redshift z $\sim$ 2, and populate projections of the stellar mass - star formation - metallicity plane, similar to observed star forming disc galaxies. We discuss small differences between these projections in observational data and in simulations, and the possible causes for the discrepancies. The light-weighted stellar masses are in good agreement with the simulation values, the differences between the two varying between 0.06 dex and 0.20 dex. We also find a good agreement between the star formation rate tracer and the true (time-averaged) simulation star formation rates. Regardless if we use mass- or light-weighted quantities, our simulations indicate that bursty star formation cycles can account for the scatter in the star forming main sequence.

We introduce a new model for the spectral energy distribution of galaxies, GRASIL-3D, which includes a careful modelling of the dust component of the interstellar medium. GRASIL-3D is an entirely new model based on the formalism of an existing and widely applied spectrophotometric model, GRASIL, but specifically designed to be interfaced with galaxies with any arbitrarily given geometry, such as galaxies calculated by theoretical hydrodynamical galaxy formation codes. GRASIL-3D is designed to separately treat radiative transfer in molecular clouds and in the diffuse cirrus component. The code has a general applicability to the outputs of simulated galaxies, either from Lagrangian or Eulerian hydrodynamic codes. As an application, the new model has been interfaced to the P-DEVA and GASOLINE smoothed-particle hydrodynamic codes, and has been used to calculate the spectral energy distribution for a variety of simulated galaxies from UV to sub-millimeter wavelengths, whose comparison with observational data gives encouraging results. In addition, GRASIL-3D allows 2D images of such galaxies to be obtained, at several angles and in different bands.

We use deep Herschel PACS and SPIRE observations in GOODSS, GOODSN and COSMOS to estimate the average dust mass (Mdust) of galaxies on a redshift-stellar mass (Mstar)-SFR grid. We study the scaling relations between Mdust, Mstar and SFR at z<=2.5. No clear evolution of Mdust is observed at fixed SFR and Mstar. We find a tight correlation between SFR and Mdust, likely a consequence of the Schmidt-Kennicutt (S-K) law. The Mstar-Mdust correlation observed by previous works flattens or sometimes disappears when fixing the SFR. Most of it likely derives from the combination of the Mdust-SFR and Mstar-SFR correlations. We then investigate the gas content as inferred by converting Mdust by assuming that the dust/gas ratio scales linearly with the gas metallicity. All galaxies in the sample follow, within uncertainties, the same SFR-Mgas relation (integrated S-K law), which broadly agrees with CO-based results for the bulk of the population, despite the completely different approaches. The majority of galaxies at z~2 form stars with an efficiency (SFE=SFR/Mgas) ~5 times higher than at z~0. It is not clear what fraction of such variation is an intrinsic redshift evolution and what fraction arises from selection effects. The gas fraction (fgas) decreases with Mstar and increases with SFR, and does not evolve with z at fixed Mstar and SFR. We explain these trends by introducing a universal relation between fgas, Mstar and SFR, non-evolving out to z~2.5. Galaxies move across this relation as their gas content evolves in time. We use the 3D fundamental fgas-Mstar-SFR relation and the redshift evolution of the Main Sequence to estimate the evolution of fgas in the average population of galaxies as a function of z and Mstar, and we find evidence a downsizing scenario.

We use the stacked gravitational lensing mass profile of four high-mass (M >~10^{15}Msun) galaxy clusters around z ~ 0.3 from Umetsu et al. to fit density profiles of phenomenological [Navarro-Frenk-White (NFW), Einasto, S\'ersic, Stadel, Baltz-Marshall-Oguri (BMO) and Hernquist] and theoretical (non-singular Isothermal Sphere, DARKexp and Kang & He) models of the dark matter distribution. We account for large-scale structure effects, including a 2-halo term in the analysis. We find that the BMO model provides the best fit to the data as measured by the reduced chi^2. It is followed by the Stadel profile, the generalized NFW profile with a free inner slope and by the Einasto profile. The NFW model provides the best fit if we neglect the 2-halo term, in agreement with results from Umetsu et al. Among the theoretical profiles, the DARKexp model with a single form parameter has the best performance, very close to that of the BMO profile. This may indicate a connection between this theoretical model and the phenomenology of dark matter halos, shedding light on the dynamical basis of empirical profiles which emerge from numerical simulations.

The effect of the scalar spectral index on inflationary super-Hubble waves is to amplify/damp large wavelengths according to whether the spectrum is red ($n_{s}<1$) or blue ($n_{s}>1$). As a consequence, the large-scale temperature correlation function will unavoidably change sign at some angle if our spectrum is red, while it will always be positive if it is blue. We show that this inflationary filtering property also affects our estimates of the size of the homogeneous patch of the universe through the Grishchuk-Zel'dovich effect. Using the recent quadrupole measurement of ESA's Planck mission, we find that the homogeneous patch of universe is at least 87 times bigger than our visible universe if we accept Planck's best fit value $n_{s}=0.9624$. An independent estimation of the size of the universe could be used to independently constrain $n_{s}$, thus narrowing the space of inflationary models.

Collisionless shocks are ubiquitous in astrophysics and in the lab. Recent numerical simulations and experiments have shown how they can arise from the encounter of two collisionless plasma shells. When the shells interpenetrate, the overlapping region turns unstable, triggering the shock formation. As a first step towards a microscopic understanding of the process, we analyze here in detail the initial instability phase. On the one hand, 2D relativistic PIC simulations are performed where two symmetric initially cold pair plasmas collide. On the other hand, the instabilities at work are analyzed, as well as the field at saturation and the seed field which gets amplified. For mildly relativistic motions and onward, Weibel modes govern the linear phase. We derive an expression for the duration of the linear phase in good agreement with the simulations. This saturation time constitutes indeed a lower-bound for the shock formation time.

As a contribution to the study of the habitability of extrasolar planets, we implemented a 1-D Energy Balance Model (EBM), the simplest seasonal model of planetary climate, with new prescriptions for most physical quantities. Here we apply our EBM to investigate the surface habitability of planets with an Earth-like atmospheric composition but different levels of surface pressure. The habitability, defined as the mean fraction of the planet's surface on which liquid water could exist, is estimated from the pressure-dependent liquid water temperature range, taking into account seasonal and latitudinal variations of surface temperature. By running several thousands of EBM simulations we generated a map of the habitable zone (HZ) in the plane of the orbital semi-major axis, a, and surface pressure, p, for planets in circular orbits around a Sun-like star. As pressure increases, the HZ becomes broader, with an increase of 0.25 AU in its radial extent from p=1/3 bar to p=3 bar. At low pressure, the habitability is low and varies with a; at high pressure, the habitability is high and relatively constant inside the HZ. We interpret these results in terms of the pressure dependence of the greenhouse effect, the effciency of horizontal heat transport, and the extent of the liquid water temperature range. Within the limits discussed in the paper, the results can be extended to planets in eccentric orbits around non-solar type stars. The main characteristics of the pressure-dependent HZ are modestly affected by variations of planetary properties, particularly at high pressure.

We combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $\mu$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-the-art treatment of dust extinction and reprocessing. We find that all of our galaxies appear to require massive populations of old (>1 Gyr) stars and, at the same time, to host a moderate ongoing activity of SF (SFR < 100 M$_{\odot}$/yr). The bulk of the stars appear to have been formed a few Gyr before the observation in essentially all cases. Only five galaxies of the sample require a recent starburst superimposed on a quiescent star formation history (SFH). We also find discrepancies between our results and those based on optical-only SED fitting for the same objects; by fitting their observed Spectral Energy Distributions with our physical model we find higher extinctions (by $\Delta$A_{V} ~ 0.81 and 1.14) and higher stellar masses (by $\Delta$Log(M*) ~ 0.16 and 0.36 dex) for z~1 and z~2 (U)LIRGs, respectively. The stellar mass difference is larger for the most dust obscured objects. We also find lower SFRs than those computed from L_{IR} using the Kennicutt relation due to the significant contribution to the dust heating by intermediate-age stellar populations through 'cirrus' emission (~73% and ~66% of total L_{IR} for z~1 and z~2 (U)LIRGs, respectively).

The generation of DC magnetic fields in unmagnetized plasmas with velocity shear is predicted for non relativistic and relativistic scenarios either due to thermal effects or due to the onset of the Kelvin-Helmholtz instability (KHI). A kinetic model describes the growth and the saturation of the DC field. The predictions of the theory are confirmed by multidimensional particle-in-cell simulations, demonstrating the formation of long lived magnetic fields ($t \sim 100s \omega_{pi}^{-1}$) along the full longitudinal extent of the shear layer, with transverse width on the electron length scale ($\sqrt{\gamma_0}c/\omega_{pe}$), reaching magnitudes $eB_{\mathrm{DC}}/m_ec\omega_{pe}\sim \beta_0\sqrt{\gamma_0}$.

We assess the impact of non-thermally shock-accelerated particles on the magnetohydrodynamic (MHD) jump conditions of relativistic shocks. The adiabatic constant is calculated directly from first principle particle-in-cell simulation data, enabling a semi-kinetic approach to improve the standard fluid model and allowing for an identification of the key parameters that define the shock structure. We find that the evolving upstream parameters have a stronger impact than the corrections due to non-thermal particles. We find that the decrease of the upstream bulk speed yields deviations from the standard MHD model up to 10%. Furthermore, we obtain a quantitative definition of the shock transition region from our analysis. For Weibel-mediated shocks the inclusion of a magnetic field in the MHD conservation equations is addressed for the first time.

We empirically test the relation between the SFR(LIR) derived from the infrared luminosity, LIR, and the SFR(Ha) derived from the Ha emission line luminosity using simple conversion relations. We use a sample of 474 galaxies at z = 0.06 - 0.46 with both Ha detection (from 20k zCOSMOS survey) and new far-IR Herschel data (100 and 160 {\mu}m). We derive SFR(Ha) from the Ha extinction corrected emission line luminosity. We find a very clear trend between E(B - V) and LIR that allows to estimate extinction values for each galaxy even if the Ha emission line measurement is not reliable. We calculate the LIR by integrating from 8 up to 1000 {\mu}m the SED that is best fitting our data. We compare SFR(Ha) with the SFR(LIR). We find a very good agreement between the two SFR estimates, with a slope of m = 1.01 \pm 0.03 in the SFR(LIR) vs SFR(Ha) diagram, a normalization constant of a = -0.08 \pm 0.03 and a dispersion of sigma = 0.28 dex.We study the effect of some intrinsic properties of the galaxies in the SFR(LIR)-SFR(Ha) relation, such as the redshift, the mass, the SSFR or the metallicity. The metallicity is the parameter that affects most the SFR comparison. The mean ratio of the two SFR estimators log[SFR(LIR)/SFR(Ha)] varies by approx. 0.6 dex from metal-poor to metal-rich galaxies (8.1 < log(O/H) + 12 < 9.2). This effect is consistent with the prediction of a theoretical model for the dust evolution in spiral galaxies. Considering different morphological types, we find a very good agreement between the two SFR indicators for the Sa, Sb and Sc morphologically classified galaxies, both in slope and normalization. For the Sd, irregular sample (Sd/Irr), the formal best-fit slope becomes much steeper (m = 1.62 \pm 0.43), but it is still consistent with 1 at the 1.5 sigma level, because of the reduced statistics of this sub-sample.

In this paper we present in-situ satellite data, theory and laboratory validation that show how small scale collisionless shocks and mini-magnetospheres can form on the electron inertial scale length. The resulting retardation and deflection of the solar wind ions could be responsible for the unusual "lunar swirl" patterns seen on the surface of the Moon.

We investigate the acceleration of light particles in perpendicular shocks for plasmas consisting of a mixture of leptonic and hadronic particles. Starting from the full set of conservation equations for the mixed plasma constituents, we generalize the magneto-hydrodynamical jump conditions for a multi-component plasma, including information about the specific adiabatic constants for the different species. The impact of deviations from the standard model of an ideal gas is compared in theory and particle-in-cell simulations, showing that the standard-MHD model is a good approximation. The simulations of shocks in electron-positron-ion plasmas are for the first time multi-dimensional, transverse effects are small in this configuration and 1D simulations are a good representation if the initial magnetization is chosen high. 1D runs with a mass ratio of 1836 are performed, which identify the Larmor frequency \omega_{ci} as the dominant frequency that determines the shock physics in mixed component plasmas. The maximum energy in the non-thermal tail of the particle spectra evolves in time according to a power-law proportional to t^\alpha with \alpha in the range 1/3 < \alpha < 1, depending on the initial parameters. A connection is made with transport theoretical models by Drury (1983) and Gargate & Spitkovsky (2011), which predict an acceleration time proportional to \gamma and the theory for small wavelength scattering by Kirk & Reville (2010), which predicts a behavior rather as proportional to \gamma^2. Furthermore, we compare different magnetic field orientations with B_0 inside and out of the plane, observing qualitatively different particle spectra than in pure electron-ion shocks.

Excesses on positron and electron fluxes measured by ATIC, and the PAMELA and Fermi--LAT telescopes can be explained by dark matter annihilation in our Galaxy. However, this requires large boosts on the dark matter annihilation rate. There are many possible enhancement mechanisms, such as the Sommerfeld effect or the existence of dark matter clumps in our halo. If enhancements on the dark matter annihilation cross section are taking place, the dark matter annihilation in the core of the Earth should also be enhanced. Here we use recent results from the IceCube 40-string configuration to probe generic enhancement scenarios. We present results as a function of the dark matter-proton interaction cross section, $\sigma_{\chi p}$ weighted by the branching fraction into neutrinos, $f_{\nu\bar{\nu}}$, as a function of a generic boost factor, $B_F$, which parametrizes the expected enhancement of the annihilation rate. We find that dark matter models which require annihilation enhancements of $\mathcal{O}(100)$ or more and that annihilate significantly into neutrinos are excluded as the explanation for these excesses. We also determine the boost range that can be probed by the full IceCube telescope.

The formation of non-relativistic collisionless shocks in laboratory with ultrahigh intensity lasers is studied via \emph{ab initio} multi-dimensional particle-in-cell simulations. The microphysics behind shock formation and dissipation, and the detailed shock structure are analyzed, illustrating that the Weibel instability plays a crucial role in the generation of strong subequipartition magnetic fields that isotropize the incoming flow and lead to the formation of a collisionless shock, similarly to what occurs in astrophysical scenarios. The possibility of generating such collisionless shocks in laboratory opens the way to the direct study of the physics associated with astrophysical shocks.

A serious concern for semi-analytical galaxy formation models, aiming to simulate multi-wavelength surveys and to thoroughly explore the model parameter space, is the extremely time consuming numerical solution of the radiative transfer of stellar radiation through dusty media. To overcome this problem, we have implemented an artificial neural network algorithm in the radiative transfer code GRASIL, in order to significantly speed up the computation of the infrared SED. The ANN we have implemented is of general use, in that its input neurons are defined as those quantities effectively determining the shape of the IR SED. Therefore, the training of the ANN can be performed with any model and then applied to other models. We made a blind test to check the algorithm, by applying a net trained with a standard chemical evolution model (i.e. CHE_EVO) to a mock catalogue extracted from the SAM MORGANA, and compared galaxy counts and evolution of the luminosity functions in several near-IR to sub-mm bands, and also the spectral differences for a large subset of randomly extracted models. The ANN is able to excellently approximate the full computation, but with a gain in CPU time by $\sim 2$ orders of magnitude. It is only advisable that the training covers reasonably well the range of values of the input neurons in the application. Indeed in the sub-mm at high redshift, a tiny fraction of models with some sensible input neurons out of the range of the trained net cause wrong answer by the ANN. These are extreme starbursting models with high optical depths, favorably selected by sub-mm observations, and difficult to predict a priori.

The cosmic star formation rate (CSFR), is an important clue to investigate the history of the assembly and evolution of galaxies. Here, we develop a method to study the CSFR from a purely theoretical point of view. Starting from detailed models of chemical evolution, we obtain the histories of star formation of galaxies of different morphological types. These histories are then used to determine the luminosity functions of the same galaxies by means of a spectro-photometric code. We obtain the CSFR under different hypothesis. First, we study the hypothesis of a pure luminosity evolution scenario, in which all galaxies are supposed to form at the same redshift and then evolve only in luminosity. Then we consider scenarios in which the number density or the slope of the LFs are assumed to vary with redshift. After comparison with available data we conclude that a pure luminosity evolution does not provide a good fit to the data, especially at very high redshift, although many uncertainties are still present in the data. On the other hand, a variation in the number density of ellipticals and spirals as a function of redshift can provide a better fit to the observed CSFR. We also explore cases of variable slope of the LFs with redshift and variations of number density and slope at the same time. We cannot find any of those cases which can improve the fit to the data respect to the solely number density variation. Finally, we compute the evolution of the average cosmic metallicity in galaxies with redshift.

Collisionless plasma instabilities are fundamental in magnetic field generation in astrophysical scenarios, but their role has been addressed in scenarios where velocity shear is absent. In this work we show that velocity shears must be considered when studying realistic astrophysical scenarios, since these trigger the collisionless Kelvin-Helmholtz instability (KHI). We present the first self-consistent three-dimensional (3D) particle-in-cell (PIC) simulations of the KHI in conditions relevant for unmagnetized relativistic outflows with velocity shear, such as active galactic nuclei (AGN) and gamma-ray bursts (GRBs). We show the generation of a strong large-scale DC magnetic field, which extends over the entire shear-surface, reaching thicknesses of a few tens of electron skin depths, and persisting on time-scales much longer than the electron time scale. This DC magnetic field is not captured by MHD models since it arises from intrinsically kinetic effects. Our results indicate that the KHI can generate intense magnetic fields yielding equipartition values up to \epsilon_B/\epsilon_p ~ 10^-3 in the electron time-scale. The KHI-induced magnetic fields have a characteristic structure that will lead to a distinct radiation signature, and can seed the turbulent dynamo amplification process. The dynamics of the KHI are relevant for non-thermal radiation modeling and can also have a strong impact on the formation of relativistic shocks in presence of velocity shears.

The spectral energy distribution of galaxies is a complex function of the star formation history and geometrical arrangement of stars and gas in galaxies. The computation of the radiative transfer of stellar radiation through the dust distribution is time-consuming. This aspect becomes unacceptable in particular when dealing with the predictions by semi-analytical galaxy formation models populating cosmological volumes, to be then compared with multi-wavelength surveys. Mainly for this aim, we have implemented an artificial neural network algorithm into the spectro-photometric and radiative transfer code GRASIL in order to compute the spectral energy distribution of galaxies in a short computing time. This allows to avoid the adoption of empirical templates that may have nothing to do with the mock galaxies output by models. The ANN has been implemented to compute the dust emission spectrum (the bottleneck of the computation), and separately for the star-forming molecular clouds and the diffuse dust (due to their different properties and dependencies). We have defined the input neurons effectively determining their emission, which means this implementation has a general applicability and is not linked to a particular galaxy formation model. We have trained the net for the disc and spherical geometries, and tested its performance to reproduce the SED of disc and starburst galaxies, as well as for a semi-analytical model for spheroidal galaxies. We have checked that for this model both the SEDs and the galaxy counts in the Herschel bands obtained with the ANN approximation are almost superimposed to the same quantities obtained with the full GRASIL. We conclude that this method appears robust and advantageous, and will present the application to a more complex SAM in another paper.

The presence of dust strongly affects the way we see galaxies and also the chemical abundances we measure in gas. It is therefore important to study he chemical evolution of galaxies by taking into account dust evolution. We aim at performing a detailed study of abundance ratios of high redshift objects and their dust properties. We focus on Lyman-Break galaxies (LBGs) and Quasar (QSO) hosts and likely progenitors of low- and high-mass present-day elliptical galaxies, respectively. We have adopted a chemical evolution model for elliptical galaxies taking account the dust production from low and intermediate mass stars, supernovae Ia, supernovae II, QSOs and both dust destruction and accretion processes. By means of such a model we have followed the chemical evolution of ellipticals of different baryonic masses. Our model complies with chemical downsizing. We made predictions for the abundance ratios versus metallicity trends for models of differing masses that can be used to constrain the star formation rate, initial mass function and dust mass in observed galaxies. We predict the existence of a high redshift dust mass-stellar mass relationship. We have found a good agreement with the properties of LBGs if we assume that they formed at redshift z=2-4. In particular, a non-negligible amount of dust is needed to explain the observed abundance pattern. We studied the QSO SDSS J114816, one of the most distant QSO ever observed (z=6.4), and we have been able to reproduce the amount of dust measured in this object. The dust is clearly due to the production from supernovae and the most massive AGB stars as well as from the grain growth in the interstellar medium. The QSO dust is likely to dominate only in the very central regions of the galaxies and during the early development of the galactic wind.

We present the analysis of Spitzer-IRS spectra of four early-type galaxies, NGC 1297, NGC 5044, NGC 6868, and NGC 7079, all classified as LINERs in the optical bands. Their IRS spectra present the full series of H2 rotational emission lines in the range 5--38 microns, atomic lines, and prominent PAH features. We investigate the nature and origin of the PAH emission, characterized by unusually low 6 -- 9/11.3 microns inter-band ratios. After the subtraction of a passive early type galaxy template, we find that the 7 -- 9 microns spectral region requires dust features not normally present in star forming galaxies. Each spectrum is then analyzed with the aim of identifying their components and origin. In contrast to normal star forming galaxies, where cationic PAH emission prevails, our 6--14 microns spectra seem to be dominated by large and neutral PAH emission, responsible for the low 6 -- 9/11.3 microns ratios, plus two broad dust emission features peaking at 8.2 microns and 12 microns. Theses broad components, observed until now mainly in evolved carbon stars and usually attributed to pristine material, contribute approximately 30-50% of the total PAH flux in the 6--14 microns region. We propose that the PAH molecules in our ETGs arise from fresh carbonaceous material which is continuously released by a population of carbon stars, formed in a rejuvenation episode which occurred within the last few Gyr. The analysis of the MIR spectra allows us to infer that, in order to maintain the peculiar size and charge distributions biased to large and neutral PAHs, this material must be shocked, and excited by the weak UV interstellar radiation field of our ETG.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

We use deep observations taken with the Photodetector Array Camera and Spectrometer (PACS), on board the Herschel satellite as part of the PACS evolutionary probe (PEP) guaranteed project along with submm ground-based observations to measure the dust mass of a sample of high-z submillimeter galaxies (SMGs). We investigate their dust content relative to their stellar and gas masses, and compare them with local star-forming galaxies. High-z SMGs are dust rich, i.e. they have higher dust-to-stellar mass ratios compared to local spiral galaxies (by a factor of 30) and also compared to local ultraluminous infrared galaxies (ULIRGs, by a factor of 6). This indicates that the large masses of gas typically hosted in SMGs have already been highly enriched with metals and dust. Indeed, for those SMGs whose gas mass is measured, we infer dust-to-gas ratios similar or higher than local spirals and ULIRGs. However, similarly to other strongly star-forming galaxies in the local Universe and at high-z, SMGs are characterized by gas metalicities lower (by a factor of a few) than local spirals, as inferred from their optical nebular lines, which are generally ascribed to infall of metal-poor gas. This is in contrast with the large dust content inferred from the far-IR and submm data. In short, the metalicity inferred from the dust mass is much higher (by more than an order of magnitude) than that inferred from the optical nebular lines. We discuss the possible explanations of this discrepancy and the possible implications for the investigation of the metalicity evolution at high-z.

A non-resonant instability for the amplification of the interstellar magnetic field in young Supernova Remnant (SNR) shocks was predicted by Bell (2004), and is thought to be relevant for the acceleration of cosmic ray (CR) particles. For this instability, the CRs streaming ahead of SNR shock fronts drive electromagnetic waves with wavelengths much shorter than the typical CR Larmor radius, by inducing a current parallel to the background magnetic field. We explore the nonlinear regime of the non-resonant mode using Particle-in-Cell (PIC) hybrid simulations, with kinetic ions and fluid electrons, and analyze the saturation mechanism for realistic CR and background plasma parameters. In the linear regime, the observed growth rates and wavelengths match the theoretical predictions; the nonlinear stage of the instability shows a strong reaction of both the background plasma and the CR particles, with the saturation level of the magnetic field varying with the CR parameters. The simulations with CR-to-background density ratios of n_CR/n_b=10^(-5) reveal the highest magnetic field saturation levels, with energy also being transferred to the background plasma and to the perpendicular velocity components of the CR particles. The results show that amplification factors >10 for the magnetic field can be achieved, and suggest that this instability is important for the generation of magnetic field turbulence, and for the acceleration of CR particles.

We use a model for the evolution of galaxies in the far-IR based on the LambdaCDM cosmology to make detailed predictions for upcoming cosmological surveys with the Herschel Space Observatory. We use the combined GALFORM semi-analytical galaxy formation model and GRASIL spectrophotometric code to compute galaxy SEDs including the reprocessing of radiation by dust. The model, which is the same as that in Baugh et al. (2005), assumes two different IMFs: a normal solar neighbourhood IMF for quiescent star formation in disks, and a very top-heavy IMF in starbursts triggered by galaxy mergers. We have shown previously that the top-heavy IMF appears necessary to explain the number counts and redshifts of faint sub-mm galaxies. In this paper, we present predictions for galaxy luminosity functions, number counts and redshift distributions in the Herschel imaging bands. We find that source confusion will be a serious problem in the deepest planned surveys. We also show predictions for physical properties such as star formation rates and stellar, gas and halo masses, together with fluxes at other wavelengths (from the far-UV to the radio) relevant for multi-wavelength follow-up observations. We investigate what fraction of the total IR emission from dust and of the high-mass star formation over the history of the Universe should be resolved by planned surveys with Herschel, and find a fraction ~30-50%, depending on confusion. Finally, we show that galaxies in Herschel surveys should be significantly clustered.

We study in detail the mid-infrared Spitzer-IRS spectrum of M 87 in the range 5 to 20 micron. Thanks to the high sensitivity of our Spitzer-IRS spectra we can disentangle the stellar and nuclear components of this active galaxy. To this end we have properly subtracted from the M 87 spectrum, the contribution of the underlying stellar continuum, derived from passive Virgo galaxies in our sample. The residual is a clear power-law, without any additional thermal component, with a zero point consistent with that obtained by high spatial resolution, ground based observations. The residual is independent of the adopted passive template. This indicates that the 10 micron silicate emission shown in spectra of M 87 can be entirely accounted for by the underlying old stellar population, leaving little room for a possible torus contribution. The MIR power-law has a slope alpha ~ 0.77-0.82 (S$_\nu\propto\nu^{-\alpha}$), consistent with optically thin synchrotron emission.

It has recently been shown that galaxy formation models within the LambdaCDM cosmology predict that, compared to the observed population, small galaxies (with stellar masses < 10^{11} M_sun) form too early, are too passive since z ~ 3 and host too old stellar populations at z=0. We then expect an overproduction of small galaxies at z > 4 that should be visible as an excess of faint Lyman-break galaxies. To check whether this excess is present, we use the MORGANA galaxy formation model and GRASIL spectro-photometric + radiative transfer code to generate mock catalogues of deep fields observed with HST-ACS. We add observational noise and the effect of Lyman-alpha emission, and perform color-color selections to identify Lyman-break galaxies. The resulting mock candidates have plausible properties that closely resemble those of observed galaxies. We are able to reproduce the evolution of the bright tail of the luminosity function of Lyman-break galaxies (with a possible underestimate of the number of the brightest i-dropouts), but uncertainties and degeneracies in dust absorption parameters do not allow to give strong constraints to the model. Besides, our model shows a clear excess with respect to observations of faint Lyman-break galaxies, especially of z_{850} ~ 27 V-dropouts at z ~ 5. We quantify the properties of these "excess" galaxies and discuss the implications: these galaxies are hosted in dark matter halos with circular velocities in excess of 100 km s^{-1}, and their suppression may require a deep re-thinking of stellar feedback processes taking place in galaxy formation.

We introduce a new technique based on artificial neural networks which allows us to make accurate predictions for the spectral energy distributions (SEDs) of large samples of galaxies, at wavelengths ranging from the far-ultra-violet to the sub-millimetre and radio. The neural net is trained to reproduce the SEDs predicted by a hybrid code comprised of the GALFORM semi-analytical model of galaxy formation, which predicts the full star formation and galaxy merger histories, and the GRASIL spectro-photometric code, which carries out a self-consistent calculation of the SED, including absorption and emission of radiation by dust. Using a small number of galaxy properties predicted by GALFORM, the method reproduces the luminosities of galaxies in the majority of cases to within 10% of those computed directly using GRASIL. The method performs best in the sub-mm and reasonably well in the mid-infrared and the far-ultra-violet. The luminosity error introduced by the method has negligible impact on predicted statistical distributions, such as luminosity functions or colour distributions of galaxies. We use the neural net to predict the overlap between galaxies selected in the rest-frame UV and in the observer-frame sub-mm at z=2. We find that around half of the galaxies with a 850um flux above 5 mJy should have optical magnitudes brighter than R_AB < 25 mag. However, only 1% of the galaxies selected in the rest-frame UV down to R_AB < 25 mag should have 850um fluxes brighter than 5 mJy. Our technique will allow the generation of wide-angle mock catalogues of galaxies selected at rest-frame UV or mid- and far-infrared wavelengths.

Ab-initio numerical study of collisionless shocks in electron-ion unmagnetized plasmas is performed with fully relativistic particle in cell simulations. The main properties of the shock are shown, focusing on the implications for particle acceleration. Results from previous works with a distinct numerical framework are recovered, including the shock structure and the overall acceleration features. Particle tracking is then used to analyze in detail the particle dynamics and the acceleration process. We observe an energy growth in time that can be reproduced by a Fermi-like mechanism with a reduced number of scatterings, in which the time between collisions increases as the particle gains energy, and the average acceleration efficiency is not ideal. The in depth analysis of the underlying physics is relevant to understand the generation of high energy cosmic rays, the impact on the astrophysical shock dynamics, and the consequent emission of radiation.

We fit the near-infrared to radio spectral energy distributions of a sample of 30 luminous and ultra-luminous infrared galaxies with models that include both starburst and AGN components. The aim of the work was to determine important physical parameters for this kind of objects such as the optical depth towards the luminosity source, the star formation rate, the star formation efficiency and the AGN fraction. We found that although about half of our sample have best-fit models that include an AGN component, only 30 % have an AGN which accounts for more than 10 % of the infrared luminosity whereas all have an energetically dominant starburst. Our models also determine the mass of dense molecular gas. Assuming that this mass is that traced by the HCN molecule, we reproduce the observed linear relation between HCN luminosity and infrared luminosity found by Gao and Solomon (2004). However, our derived conversion factor between HCN luminosity and the mass of dense molecular gas is a factor of 2 smaller than that assumed by these authors. Finally, we find that the star formation efficiency falls as the starburst ages.

We present photometric evolution models of galaxies, in which, in addition to the stellar component, the effects of an evolving dusty interstellar medium have been included with particular care. Starting from the work of Calura, Pipino & Matteucci (2008), in which chemical evolution models have been used to study the evolution of both the gas and dust components of the interstellar medium in the solar neighbourhood, elliptical and irregular galaxies, it has been possible to combine these models with a spectrophotometric stellar code that includes dust reprocessing (GRASIL) (Silva et al. 1998) to analyse the evolution of the spectral energy distributions (SED) of these galaxies. We test our models against observed SEDs both in the local universe and at high redshift and use them to predict how the percentage of reprocessed starlight evolves for each type of galaxy. The importance of following the dust evolution is investigated by comparing our results with those obtained by adopting simple assumptions to treat this component.

The treatment of dust attenuation is crucial in order to compare the predictions of galaxy formation models with multiwavelength observations. Most past studies have either used simple analytic prescriptions or else full radiative transfer (RT) calculations. Here, we couple star formation histories and morphologies predicted by the semi-analytic galaxy formation model MORGANA with RT calculations from the spectrophotometric and dust code GRASIL to create a library of galaxy SEDs from the UV/optical through the far Infrared, and compare the predictions of the RT calculations with analytic prescriptions. We consider a low and high redshift sample, as well as an additional library constructed with empirical, non-cosmological star formation histories and simple (pure bulge or disc) morphologies. Based on these libraries, we derive fitting formulae for the effective dust optical depth as a function of galaxy physical properties such as metallicity, gas mass, and radius. We show that such fitting formulae can predict the V-band optical depth with a scatter smaller than 0.4 dex for both the low and high redshift samples, but that there is a large galaxy-to-galaxy scatter in the shapes of attenuation curves, probably due to geometrical variations, which our simple recipe does not capture well. However, our new recipe provides a better approximation to the GRASIL results at optical wavelength than standard analytic prescriptions from the literature, particularly at high redshift.

We use 16 micron, Spitzer-IRS, blue peakup photometry of 50 early-type galaxies in the Coma cluster to define the mid-infrared colour-magnitude relation. We compare with recent simple stellar population models that include the mid-infrared emission from the extended, dusty envelopes of evolved stars. The Ks-[16] colour in these models is very sensitive to the relative population of dusty Asymptotic Giant Branch (AGB) stars. We find that the passively evolving early-type galaxies define a sequence of approximately constant age (~10 Gyr) with varying metallicity. Several galaxies that lie on the optical/near-infrared colour-magnitude relation do not lie on the mid-infrared relation. This illustrates the sensitivity of the Ks-[16] colour to age. The fact that a colour-magnitude relation is seen in the mid-infrared underlines the extremely passive nature of the majority (68%) of early-type galaxies in the Coma cluster. The corollary of this is that 32% of the early-type galaxies in our sample are not `passive', insofar as they are either significantly younger than 10 Gyr or they have had some rejuvenation episode within the last few Gyr.

We present initial results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity relation at z>3 by means of deep near-IR spectroscopy. Gas metallicities are measured, for an initial sample of nine star forming galaxies at z~3.5, by means of optical nebular lines redshifted into the near-IR. Stellar masses are accurately determined by using Spitzer-IRAC data, which sample the rest-frame near-IR stellar light in these distant galaxies. When compared with previous surveys, the mass-metallicity relation inferred at z~3.5 shows an evolution much stronger than observed at lower redshifts. The evolution is prominent even in massive galaxies, indicating that z~3 is an epoch of major action in terms of star formation and metal enrichment also for massive systems. There are also indications that the metallicity evolution of low mass galaxies is stronger relative to high mass systems, an effect which can be considered the chemical version of the galaxy downsizing. The mass-metallicity relation observed at z~3.5 is difficult to reconcile with the predictions of some hierarchical evolutionary models. Such discrepancies suggest that at z>3 galaxies are assembled mostly with relatively un-evolved sub-units, i.e. small galaxies with low star formation efficiency. The bulk of the star formation and metallicity evolution probably occurs once small galaxies are already assembled into bigger systems.

In a previous work, three bright MIR/radio sources were discovered in the nuclear region of NGC1365. We here confirm that these sources are young and massive ``embedded'' clusters, and derive their parameters, such as extinction, age and mass. Using ISAAC and VISIR at the VLT we obtained NIR and MIR maps and LR spectra. The dataset is first interpreted by comparing the observations with images and spectra of the close-by young cluster R136 in the LMC and then by using model predictions for both the nebular emission lines and the spectral energy distribution of the sources. We produce maps of the region containing the three sources in the R, J, Ks, L' bands and at 12.8micro. We also provide spectra in K, L and N. The spectral energy distribution of the three sources rises with wavelength. Emission lines from ionised hydrogen and molecular hydrogen are detected, as well as PAH emission. Conspicuous [NeII]12.8 line is also present, while neither the [ArIII] 8.9 nor the [SIV] 10.4 lines are detected. This provides a stringent constraint on the age of the sources: they are relatively evolved young clusters (6-8 Myr). Owing to their ionising photon rates and ages, they must be extremely massive clusters (around 10^7 solar masses). Their MIR spectral energy distribution suggests the presence of two components: (1) an optically thin component, with a continuum comparable to that of R136, and (2) an optically thick component which might be related to subsequent or on-going episodes of star formation. These sources are good candidates for evolving according to a bi-modal hydrodynamical regime, in which matter is trapped at the centre of a compact and massive cluster.

We present a large set of theoretical isochrones, whose distinctive features mostly reside on the greatly improved treatment of the thermally pulsing asymptotic giant branch (TP-AGB) phase. Essentially, we have coupled the TP-AGB tracks described in Paper I, at their stages of pre-flash quiescent H-shell burning, with the evolutionary tracks for the previous evolutionary phases from Girardi et al. (2000). Theoretical isochrones for any intermediate value of age and metallicity are then derived by interpolation in the grids. We take care that the isochrones keep, to a good level of detail, the several peculiarities present in these TP-AGB tracks. Theoretical isochrones are then converted to about 20 different photometric systems -- including traditional ground-based systems, and those of recent major wide-field surveys such as SDSS, OGLE, DENIS, 2MASS, UKIDSS, etc., -- by means of synthetic photometry applied to an updated library of stellar spectra, suitably extended to include C-type stars. Finally, we correct the predicted photometry by the effect of circumstellar dust during the mass-losing stages of the AGB evolution, which allows us to improve the results for the optical-to-infrared systems, and to simulate mid- and far-IR systems such as those of Spitzer and AKARI. Access to the data is provided both via a web repository of static tables (http://stev.oapd.inaf.it/dustyAGB07 and CDS), and via an interactive web interface (http://stev.oapd.inaf.it/cmd) that provides tables for any intermediate value of age and metallicity, for several photometric systems, and for different choices of dust properties.

We present predictions for the evolution of the galaxy luminosity function, number counts and redshift distributions in the IR based on the Lambda-CDM cosmological model. We use the combined GALFORM semi-analytical galaxy formation model and GRASIL spectrophotometric code to compute galaxy SEDs including the reprocessing of radiation by dust. The model, which is the same as that in Baugh et al (2005), assumes two different IMFs: a normal solar neighbourhood IMF for quiescent star formation in disks, and a very top-heavy IMF in starbursts triggered by galaxy mergers. We have shown previously that the top-heavy IMF seems to be necessary to explain the number counts of faint sub-mm galaxies. We compare the model with observational data from the SPITZER Space Telescope, with the model parameters fixed at values chosen before SPITZER data became available. We find that the model matches the observed evolution in the IR remarkably well over the whole range of wavelengths probed by SPITZER. In particular, the SPITZER data show that there is strong evolution in the mid-IR galaxy luminosity function over the redshift range z ~ 0-2, and this is reproduced by our model without requiring any adjustment of parameters. On the other hand, a model with a normal IMF in starbursts predicts far too little evolution in the mid-IR luminosity function, and is therefore excluded.

We present preliminary results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity relation at z~3 by means of deep near-IR spectroscopy. Gas metallicities and stellar masses are measured for an initial sample of nine star forming galaxies at z~3.3. When compared with previous surveys, the mass-metallicity relation inferred at z~3.3 shows an evolution significantly stronger than observed at lower redshifts. There are also some indications that the metallicity evolution of low mass galaxies is stronger relative to high mass systems, an effect which can be considered as the chemical version of the galaxy downsizing. The mass-metallicity relation observed at z~3.3 is difficult to reconcile with the predictions of some hierarchical evolutionary models. We shortly discuss the possible implications of such discrepancies.

We fit the near-infrared to radio spectral energy distributions of 30 luminous and ultra-luminous infrared galaxies with pure starburst models or models that include both starburst and AGN components to determine important physical parameters for this population of objects. In particular we constrain the optical depth towards the luminosity source, the star formation rate, the star formation efficiency and the AGN fraction. We find that although about half of our sample have best-fit models that include an AGN component, only 30% have an AGN which accounts for more than 10% of the infrared luminosity, whereas all have an energetically dominant starburst. Our derived AGN fractions are generally in good agreement other measurements based in the mid-infrared line ratios measured by Spitzer IRS, but lower than those derived from PAH equivalent widths or the mid-infrared spectral slope. Our models determine the mass of dense molecular gas via the extinction required to reproduce the SED. Assuming that this mass is that traced by HCN, we reproduce the observed linear relation between HCN and infrared luminosities found by Gao & Solomon. We also find that the star formation efficiency, defined as the current star formation rate per unit of dense molecular gas mass, is enhanced in the ULIRGs phase. If the evolution of ULIRGs includes a phase in which an AGN contributes an important fraction to the infrared luminosity, this phase should last an order of magnitude less time than the starburst phase. Because the mass of dense molecular gas which we derive is consistent with observations of the HCN molecule,it should be possible to estimate the mass of dense, star-forming molecular gas in such objects when molecular line data are not available.

As a constraint for new starburst/AGN models of IRAS bright galaxies we determine the radio spectra of 31 luminous and ultraluminous IRAS galaxies (LIRGs/ULIRGs). We construct the radio spectra using both new and archival data. From our sample of radio spectra we find that very few have a straight power-law slope. Although some sources show a flattening of the radio spectral slope at high frequencies the average spectrum shows a steepening of the radio spectrum from 1.4 to 22.5 GHz. This is unexpected because in sources with high rates of star formation we expect flat spectrum, free-free emission to make a significant contribution to the radio flux at higher radio frequencies. Despite this trend the radio spectral indices between 8.4 and 22.5 GHz are flatter for sources with higher values of the FIR-radio flux density ratio q, when this is calculated at 8.4 GHz. Therefore, sources that are deficient in radio emission relative to FIR emission (presumably younger sources) have a larger thermal component to their radio emission. However, we find no correlation between the radio spectral index between 1.4 and 4.8 GHz and q at 8.4 GHz. Because the low frequency spectral index is affected by free-free absorption, and this is a function of source size for a given mass of ionized gas, this is evidence that the ionized gas in ULIRGs shows a range of densities. The youngest LIRGs and ULIRGs are characterized by a larger contribution to their high-frequency radio spectra from free-free emission. However, the youngest sources are not those that have the greatest free-free absorption at low radio frequencies. The sources in which the effects of free-free absorption are strongest are instead the most compact sources. Although these have the warmest FIR colours, they are not necessarily the youngest sources.

In order to gain insight into the physical mechanisms leading to the formation of stars and their assembly in galaxies, we compare the predictions of the MOdel for the Rise of GAlaxies aNd Active nuclei (MORGANA) to the properties of K- and 850 micron-selected galaxies (such as number counts, redshift distributions and luminosity functions) by combining MORGANA with the spectrophotometric model GRASIL. We find that it is possible to reproduce the K- and 850 micron-band datasets at the same time and with a standard Salpeter IMF, and ascribe this success to our improved modeling of cooling in DM halos. We then predict that massively star-forming discs are common at z~2 and dominate the star-formation rate, but most of them merge with other galaxies within ~100 Myr. Our preferred model produces an overabundance of bright galaxies at z<1; this overabundance might be connected to the build-up of the diffuse stellar component in galaxy clusters, as suggested by Monaco et al. (2006), but a naive implementation of the mechanism suggested in that paper does not produce a sufficient slow-down of the evolution of these objects. Moreover, our model over-predicts the number of 10^{10}-10^{11} M_sun galaxies at z~1; this is a common behavior of theoretical models as shown by Fontana et al. (2006). These findings show that, while the overall build-up of the stellar mass is correctly reproduced by galaxy formation models, the ``downsizing'' trend of galaxies is not fully reproduced yet. This hints to some missing feedback mechanism in order to reproduce at the same time the formation of both the massive and the small galaxies.

We present the redshift distribution of the SHADES galaxy population based on the rest-frame radio-mm-FIR colours of 120 robustly detected 850um sources in the Lockman Hole East (LH) and Subaru XMM-Newton Deep Field (SXDF). The redshift distribution derived from the full SED information is shown to be narrower than that determined from the radio-submm spectral index, as more photometric bands contribute to a higher redshift accuracy. The redshift distribution of sources derived from at least two photometric bands peaks at z ~ 2.4 and has a near-Gaussian distribution, with 50 per cent (interquartile range) of sources at z=1.8-3.1. We find a statistically-significant difference between the measured redshift distributions in the two fields; the SXDF peaking at a slightly lower redshift (median z ~ 2.2) than the LH (median z ~ 2.7), which we attribute to the noise-properties of the radio observations. We demonstrate however that there could also be field-to-field variations that are consistent with the measured differences in the redshift distributions, and hence, that the incomplete area observed by SHADES with SCUBA, despite being the largest sub-mm survey to date, may still be too small to fully characterize the bright sub-mm galaxy population. Finally we present a brief comparison with the predicted, or assumed, redshift distributions of sub-mm galaxy formation and evolution models, and we derive the contribution of these SHADES sources and the general sub-mm galaxy population to the star formation-rate density at different epochs.

We present a quantitative analysis of the astrophysical and cosmological information that can be extracted from the many important wide-area, shallow surveys that will be carried out in the next few years. Our calculations combine the predictions of the physical model by Granato et al. (2004) for the formation and evolution of spheroidal galaxies with up-to-date phenomenological models for the evolution of starburst and normal late-type galaxies and of radio sources. We compute the expected number counts and the redshift distributions of these source populations separately and then focus on proto-spheroidal galaxies. For the latter objects we predict the counts and redshift distributions of strongly lensed sources at 250, 350, 500, and 850 micron, the angular correlation function of sources detected in the surveys considered, the angular power spectra due to clustering of sources below the detection limit in Herschel and Planck surveys. An optimal survey for selecting strongly lensed proto-spheroidal galaxies is described, and it is shown how they can be easily distinguished from the other source populations. We also discuss the detectability of the imprints of the 1-halo and 2-halo regimes on angular correlation functions and clustering power spectra, as well as the constraints on cosmological parameters that can be obtained from the determinations of these quantities. The novel data relevant to derive the first sub-millimeter estimates of the local luminosity functions of starburst and late-type galaxies, and the constraints on the properties of rare source populations, such as blazars, are also briefly described.

The mid infrared emission of early type galaxies traces the presence of intermediate age stellar populations as well as even tiny amounts of ongoing star formation. Here we discuss high S/N Spitzer IRS spectra of a sample of Virgo early type galaxies, with particular reference to NGC 4435. We show that, by combining mid infrared spectroscopic observations with existing broad band fluxes, it is possible to obtain a very clean picture of the nuclear activity in this galaxy.

We are performing a systematic study of the properties of early-type galaxies in the mid infrared spectral region with the Spitzer Space Telescope. We present here high S/N Spitzer IRS spectra of 17 Virgo early-type galaxies. Thirteen objects of the sample (76%) show a pronounced broad feature (above 10 microns) which is spatially extended and likely of stellar origin. We argue that this feature is (mostly) due to silicate emission from circumstellar envelopes of asymptotic giant branch (AGB) stars. The remaining 4 objects, namely NGC 4486, NGC 4636, NGC 4550 and NGC 4435, are characterized by various levels and type of activity.

We analyse the attenuation properties of a sample of UV selected galaxies, with the use of the spectrophotometric model Grasil. In particular, we focus on the relation between dust attenuation and the reddening in the UV spectral region. We show that a realistic modelling of geometrical distribution of dust and of the different population of stars can explain the UV reddening of normal spiral galaxies also with a standard Milky Way dust. Our results clearly underline that it is fundamental to take into account that younger stars suffer a higher attenuation than older stars (the age-dependent extinction) because stars are born in more-than-average dusty environments. In this work we also find that the concentration of young stars on the galactic plane of spirals has a relevant impact on the expected UV colours, impact that has not been explored before this paper. Finally, we discuss the role of IMF in shaping the relation between UV reddening and dust attenuation, and we show that a Kroupa IMF is more consistent with observed data than the classical Salpeter IMF.

The ~800 optically unseen (R>25.5) 24mum-selected sources in the complete Spitzer First Look Survey sample (Fadda et al. 2006) with F[24mum]>0.35 mJy are found to be very strongly clustered. If, as indicated by several lines of circumstantial evidence, they are ultraluminous far-IR galaxies at z ~ [1.6-2.7], the amplitude of their spatial correlation function is very high. The associated comoving clustering length is estimated to be r_0=14.0_{-2.4}^{+2.1} Mpc, value which puts these sources amongst the most strongly clustered populations of our known universe. Their 8mum-24mum colours suggest that the AGN contribution dominates above F[24mum] ~ 0.8 mJy, consistent with earlier analyses. The properties of these objects (number counts, redshift distribution, clustering amplitude) are fully consistent with those of proto-spheroidal galaxies in the process of forming most of their stars and of growing their active nucleus, as described by the Granato et al. (2004) model. In particular, the inferred space density of such galaxies at z ~ 2 is much higher than what expected from most semi-analytic models. Matches of the observed projected correlation function w(\theta) with models derived within the so-called Halo Occupation Scenario show that these sources have to be hosted by haloes more massive than ~10^{13.4} M_\odot. This value is significantly higher than that for the typical galactic haloes hosting massive elliptical galaxies, suggesting a duration of the starburst phase of massive high-redshift dusty galaxies of T_B ~ 0.5 Gyr.

Many problems in astrophysics involve relativistic outflows. The plasma dynamics in these scenarios is critical to determine the conditions for the self-consistent evolution of the fields and particle acceleration. Advances in computer power now allow for kinetic plasma simulations, based on the particle-in-cell (PIC) paradigm, capable of providing information about the role of plasma instabilities in relativistic outflows. A discussion of the key issues associated with PIC simulations is presented, along with some the most important results and open questions, with a particular emphasis on the long time evolution of the filamentation, or Weibel, instability, and on the possible collisionless mechanisms for particle acceleration arising in the collision of relativistic plasma shells.

We present a population synthesis study of NGC4435, an early-type Virgo galaxy interacting with NGC4438. We combine new spectroscopic observations obtained with the Spitzer Space Telescope IRS instrument with IRAC archival data and broad band data from the literature. The IRS spectrum shows prominent PAH features, low ionization emission lines and H_2 rotational lines arising from the dusty circumnuclear disk characterizing this galaxy. The central SED, from X-ray to radio, is well fitted by a model of an exponential burst superimposed on an old simple stellar population. From the lack of high excitation nebular lines, the [NeIII]15.5/[NeII]12.8 ratio, the temperature of molecular hydrogen, and the fit to the full X-ray to radio SED we argue that the present activity of the galaxy is driven by star formation alone. The AGN contribution to the ionizing flux is constrained to be less than 2%. The age of the burst is found to be around 190 Myr and it is fully consistent with the notion that the star formation process has been triggered by the interaction with NGC4438. The mass involved in the rejuvenation episode turns out to be less than 1.5% of the stellar galaxy mass sampled in a 5" central aperture. This is enough to render NGC4435 closely similar to a typical interacting early-type galaxy with inverted CaII[H+K] lines that will later turn into a typical cluster E+A galaxy and enforces the notion that these objects are the result of a recent rejuvenation episode rather than a genuine delayed formation.

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

The dynamics of two counter-streaming electron-positron-ion unmagnetized plasma shells with zero net charge is analyzed in the context of magnetic field generation in GRB internal shocks due to the Weibel instability. The effects of large thermal motion of plasma particles, arbitrary mixture of plasma species and space charge effects are taken into account. We show that, although thermal effects slow down the instability, baryon loading leads to a non-negligible growth rate even for large temperatures and different shell velocities, thus guaranteeing the robustness and the occurrence of the Weibel instability for a wide range of scenarios.

We show that our previously proposed anti-hierarchical baryon collapse scenario for the joint evolution of black holes and host galaxies predicts quasar luminosity functions at redshifts 1.5<z<6 and local properties in nice agreement with observations. In our model the quasar activity marks and originates the transition between an earlier phase of violent and heavily dust-enshrouded starburst activity promoting rapid black hole growth, and a later phase of almost passive evolution; the former is traced by the submillimeter-selected sources, while the latter accounts for the high number density of massive galaxies at substantial redshifts z>1.5, the population of Extremely Red Objects, and the properties of local ellipticals.

We have obtained high S/N Spitzer IRS spectra of 17 Virgo early-type galaxies that lie on the colour-magnitude relation of passively evolving galaxies in the cluster. To flux calibrate these extended sources we have devised a new procedure that allows us to obtain the intrinsic spectral energy distribution and to disentangle resolved and unresolved emission within the same object. Thirteen objects of the sample (76%) show a pronounced broad silicate feature (above 10micron) which is spatially extended and likely of stellar origin, in agreement with model predictions. The other 4 objects (24%) are characterized by different levels of activity. In NGC 4486 (M87) the line emission and the broad silicate emission are evidently unresolved and, given also the typical shape of the continuum, they likely originate in the nuclear torus. NGC 4636 show emission lines superimposed to extended silicate emission (i.e. likely of stellar origin, pushing the percentage of galaxies with silicate emission to 82%). Finally NGC 4550 and NGC 4435 are characterized by PAH and line emission, arising from a central unresolved region.

We present a hierarchical galaxy formation model which can account for the number counts of sources detected through their emission at sub-millimetre wavelengths. The first stage in our approach is an ab initio calculation of the star formation histories for a representative sample of galaxies, which is carried out using the semi-analytical galaxy formation model GALFORM. These star formation histories are then input into the spectro-photometric code GRASIL, to produce a spectral energy distribution for each galaxy. Dust extinction and emission are treated self consistently in our model, without having to resort to ad-hoc assumptions about the amount of attenuation by dust or the temperature at which the dust radiates. We argue that it is necessary to modify the form of the stellar initial mass function in starbursts in order to match the observed number of sub-mm sources, if we are to retain the previous good matches enjoyed between observations and model predictions in the local universe. We also list some other observational tests that have been passed by our model.

We show that the ABC scenario we proposed for the co-evolution of spheroids and QSOs predicts accretion rates and masses of supermassive black holes in sub-mm galaxies in keeping with recent X-ray determinations. These masses are well below the local values, and those predicted by alternative models. The observed column densities may be mostly due to ISM in the galaxy. The contribution of the associated nuclear activity to the X-ray background is likely negligible, while they may contribute a sizeable fraction $\sim 10 %$ to hard-X cumulative counts at the faintest observed fluxes.

At variance with most semi-analytic models, in the Anti-hierarchical Baryon Collapse scenario (Granato et al. 2001, 2004) the main driver of the galaxy formation and evolution is not the merging sequence but are baryon processes. This approach emphasizes, still in the framework of the hierarchical clustering paradigm for dark matter halos, feedback processes from supernova explosions and from active nuclei, that tie together star formation in spheroidal galaxies and the growth of black holes at their centers. We review some recent results showing the remarkably successful predictive power of this scenario, which allows us to account for the evolution with cosmic time of a broad variety of properties of galaxies and active nuclei, which proved to be very challenging for competing models.

We present high signal to noise ratio Spitzer Infrared Spectrograph observations of 17 Virgo early-type galaxies. The galaxies were selected from those that define the colour-magnitude relation of the cluster, with the aim of detecting the silicate emission of their dusty, mass-losing evolved stars. To flux calibrate these extended sources we have devised a new procedure that allows us to obtain the intrinsic spectral energy distribution and to disentangle resolved and unresolved emission within the same object. We have found that thirteen objects of the sample (76%) are passively evolving galaxies with a pronounced broad silicate feature which is spatially extended and likely of stellar origin, in agreement with model predictions. The other 4 objects (24%) are characterized by different levels of activity. In NGC 4486 (M 87) the line emission and the broad silicate emission are evidently unresolved and, given also the typical shape of the continuum, they likely originate in the nuclear torus. NGC 4636 shows emission lines superimposed on extended (i.e. stellar) silicate emission, thus pushing the percentage of galaxies with silicate emission to 82%. Finally, NGC 4550 and NGC 4435 are characterized by polycyclic aromatic hydrocarbon (PAH) and line emission, arising from a central unresolved region. A more detailed analysis of our sample, with updated models, will be presented in a forthcoming paper.

The kinetic theory of collisionless electrostatic shocks resulting from the collision of plasma slabs with different temperatures and densities is presented. The theoretical results are confirmed by self-consistent particle-in-cell simulations, revealing the formation and stable propagation of electrostatic shocks with very high Mach numbers ($M \gg 10$), well above the predictions of the classical theories for electrostatic shocks.

The origin of the micro-Gauss magnetic fields in galaxy clusters is one of the outstanding problem of modern cosmology. We have performed three-dimensional particle-in-cell simulations of the nonrelativistic Weibel instability in an electron-proton plasma, in conditions typical of cosmological shocks. These simulations indicate that cluster fields could have been produced by shocks propagating through the intergalactic medium during the formation of large-scale structure or by shocks within the cluster. The strengths of the shock-generated fields range from tens of nano-Gauss in the intercluster medium to a few micro-Gauss inside galaxy clusters.

In this paper we analyze the evolution of actively star forming galaxies in the mid-infrared (MIR). This spectral region, characterized by continuum emission by hot dust and by the presence of strong emission features generally ascribed to polycyclic aromatic hydrocarbon (PAH) molecules, is the most strongly affected by the heating processes associated with star formation and/or active galactic nuclei (AGN). Following the detailed observational characterization of galaxies in the MIR by ISO, we have updated the modelling of this spectral region in our spectro-photometric model GRASIL (Silva et al. 1998). In the diffuse component we have updated the treatment of PAHs according to the model by Li & Draine (2001). As for the dense phase of the ISM associated with the star forming regions, the molecular clouds, we strongly decrease the abundance of PAHs as compared to that in the cirrus, basing on the observational evidences of the lack or weakness of PAH bands close to the newly formed stars, possibly due to the destruction of the molecules in strong UV fields. The robustness of the model is checked by fitting near infrared to radio broad band spectra and the corresponding detailed MIR spectra of a large sample of galaxies (Lu et al. 2003), at once. With this model, we have analyzed the larger sample of actively star forming galaxies by Dale et al. (2000). We show that the observed trends of galaxies in the ISO-IRAS-Radio color-color plots can be interpreted in terms of different evolutionary phases of star formation activity, and the consequent different dominance in the spectral energy distribution (SED) of the diffuse or dense phase of the ISM.

We present new results from a multi-wavelength model of galaxy formation, which combines a semi-analytical treatment of the formation of galaxies within the CDM framework with a sophisticated treatment of absorption and emission of radiation by dust. We find that the model, which incorporates a top-heavy IMF in bursts, agrees well with the evolution of the rest-frame far-UV luminosity function over the range z=0-6, with the IR number counts in all bands measured by SPITZER, and with the observed evolution of the mid-IR luminosity function for z=0-2.

We summarize our physical model for the early co-evolution of spheroidal galaxies and of active nuclei at their centers. Our predictions are in excellent agreement with a number of observables which proved to be extremely challenging for all the current semi-analytic models, including the sub-mm counts and the corresponding redshift distributions, and the epoch-dependent K-band luminosity function of spheroidal galaxies. Also, the black hole mass function and the relationship between the black hole mass and the velocity dispersion in the galaxy are nicely reproduced. The mild AGN activity revealed by X-ray observations of SCUBA sources is in keeping with our scenario, and testify the build up of SMBH triggered by intense star formation.

The possible contributions of the various classes of extragalactic sources (including, in addition to the canonical radio sources, GHz Peaked Spectrum sources, advection-dominated sources, starburst galaxies, high-redshift proto-spheroidal galaxies) to the arcminute scale fluctuations measured by the CBI, BIMA, and ACBAR experiments are discussed. At 30 GHz, fluctuations due to radio sources undetected by ancillary low-frequency surveys may be higher than estimated by the CBI and BIMA groups. High-redshift dusty galaxies, whose fluctuations may be strongly enhanced by the effect of clustering, could contribute to the BIMA excess signal, and dominate at 150 GHz (the ACBAR frequency). Moreover, in the present data situation, the dust emission of these high-redshift sources set an unavoidable limit to the detection of primordial CMB anisotropies at high multipoles, even at frequencies as low as $\simeq 30$ GHz. It is concluded that the possibility that the excess power at high multipoles is dominated by unsubtracted extragalactic sources cannot be ruled out. On the other hand, there is room for a contribution from the Sunyaev-Zeldovich effect within clusters of galaxies, with a density fluctuation amplitude parameter $\sigma_8$ consistent with the values preferred by current data.

We present here a model for simulating the panchromatic spectral energy distribution of galaxies, which aims to be a complete tool to study the complex multi-wavelength picture of the universe. The model take into account all important components that concur to the SED of galaxies at wavelengths from X-rays to the radio. We review the modeling of each component and provide several applications, interpreting observations of galaxy of different types at all the wavelengths.

We present predictions for the abundance of sub-mm galaxies (SMGs) and Lyman-break galaxies (LBGs) in the $\Lambda$CDM cosmology. A key feature of our model is the self-consistent calculation of the absorption and emission of radiation by dust. The new model successfully matches the LBG luminosity function, as well reproducing the properties of the local galaxy population in the optical and IR. The model can also explain the observed galaxy number counts at $850\mum$, but only if we assume a top-heavy IMF for the stars formed in bursts. The predicted redshift distribution of SMGs depends relatively little on their flux over the range 1-$10\mjy$, with a median value of $z\approx 2.0$ at a flux of $5\mjy$, in very good agreement with the recent measurement by Chapman et al The counts of SMGs are predicted to be dominated by ongoing starbursts. However, in the model these bursts are responsible for making only a few per cent of the stellar mass locked up in massive ellipticals at the present day.

Granato et al (2004) have elaborated a physically grounded model exploiting the mutual feedback between star-forming spheroidal galaxies and the active nuclei growing in their cores to overcome, in the framework of the hierarchical clustering scenario for galaxy formation, one of the main challenges facing such scenario, the fact that massive spheroidal galaxies appear to have formed much earlier and faster than predicted by previous semi-analytical models. After having assessed the values of the two parameters that control the effect of the complex radiative transfer processes on the time-dependent SEDs we have compared the model predictions with a variety of IR to mm data. Our results support a rather strict continuity between objects where stars formed and evolved massive early-type galaxies, indicating that large spheroidal galaxies formed most of their stars when they were already assembled as single objects. The model is successful in reproducing the observed z distribution of K\le20 galaxies at z>1, in contrast with both the classical monolithic and the semi-analytic models, the ratio of star-forming to passively evolving spheroids and the counts and z distributions of EROs. The model also favourably compares with the ISOCAM 6.7 mu counts, with the corresponding z distribution, and with IRAC counts, which probe primarily the passive evolution phase, and with the submm SCUBA and MAMBO data, probing the active star-formation phase. The observed fraction of 24mu selected sources with no detectable emission in either the 8mu or R band nicely corresponds to the predicted surface density of star-forming spheroids with 8mu fluxes below the detection limit. Predictions for the z distributions of 24mu sources detected by MIPS surveys are pointed out. [Abridged]

Although the primary goal of ESA's Planck mission is to produce high resolution maps of the temperature and polarization anisotropies of the Cosmic Microwave Background (CMB), its high-sensitivity all-sky surveys of extragalactic sources at 9 frequencies in the range 30--860 GHz will constitute a major aspect of its science products. In particular, Planck surveys will provide key information on several highly interesting radio source populations, such as Flat Spectrum Radio Quasars, BL Lac objects, and, especially, extreme GHz Peaked Spectrum sources, thought to correspond to the very earliest phases of the evolution of radio sources. Above 100 GHz, Planck will provide the first all-sky surveys, that are expected to supply rich samples of highly gravitationally amplified dusty proto-galaxies and large samples of candidate proto-clusters at z~2-3, thus shedding light on the evolution of large scale structure across the cosmic epoch when dark energy should start dominating the cosmic dynamics.

We present detailed predictions of the contributions of the various source populations to the counts at frequencies of tens of GHz. New evolutionary models are worked out for flat-spectrum radio quasars, BL Lac objects, and steep-spectrum sources. Source populations characterized by spectra peaking at high radio frequencies, such as extreme GPS sources, ADAF/ADIOS sources and early phases of gamma-ray burst afterglows are also dealt with. The counts of different populations of star-forming galaxies (normal spirals, starbursts, high-z galaxies detected by SCUBA and MAMBO surveys, interpreted as proto-spheroidal galaxies) are estimated taking into account both synchrotron and free-free emission, and dust re-radiation. Our analysis is completed by updated counts of Sunyaev-Zeldovich effects in clusters of galaxies and by a preliminary estimate of galactic-scale Sunyaev-Zeldovich signals associated to proto-galactic plasma.

We present a physical model for the coevolution of massive spheroidal galaxies and active nuclei at their centers. Supernova heating is increasingly effective in slowing down the star formation and in driving gas outflows in smaller and smaller dark matter halos. Thus the more massive protogalaxies virializing at early times are the sites of faster star formation. The correspondingly higher radiation drag causes a faster angular momentum loss by the gas and induces a larger accretion rate onto the central black hole. In turn, the kinetic energy of the outflows powered by the active nuclei can unbind the residual gas in a time shorter for larger halos. The model accounts for a broad variety of dynamical, photometric and metallicity properties of early-type galaxies, for the M_BH -- \sigma relation and for the local supermassive black-hole mass function.

We investigate the long-time evolution of magnetic fields generated by the two-stream instability at ultra- and sub-relativistic astrophysical collisionless shocks. Based on 3D PIC simulation results, we introduce a 2D toy model of interacting current filaments. Within the framework of this model, we demonstrate that the field correlation scale in the region far downstream the shock grows nearly as the light crossing time, lambda(t) ~ ct, thus making the diffusive field dissipation inefficient. The obtained theoretical scaling is tested using numerical PIC simulations. This result extends our understanding of the structure of collisionless shocks in gamma-ray bursts and other astrophysical objects.

We estimate the contribution of AGNs and of their host galaxies to the infrared background. We use the luminosity function and evolution of AGNs recently determined by the hard X-ray surveys, and new Spectral Energy Distributions connecting the X-ray and the infrared emission, divided in intervals of absorption. These two ingredients allow us to determine the contribution of AGNs to the infrared background by using mostly observed quantities, with only minor assumptions. We obtain that AGN emission contributes little to the infrared background ($<$5% over most of the infrared bands), implying that the latter is dominated by star formation. However, AGN host galaxies may contribute significantly to the infrared background, and more specifically 10--20% in the 1--20$\mu$m range and $\sim$5% at $\lambda<60\mu m$. We also give the contribution of AGNs and of their host galaxies to the source number counts in various infrared bands, focusing on those which will be observed with Spitzer. We also report a significant discrepancy between the expected contribution of AGN hosts to the submm background and bright submm number counts with the observational constraints. We discuss the causes and implications of this discrepancy and the possible effects on the Spitzer far-IR bands.

We present detailed predictions for the confusion noise due to extragalactic sources in the far-IR/(sub)-millimeter channels of ESA/ISO, NASA/Spitzer, ESA/Herschel and ESA/Planck satellites, including the contribution from clustering of unresolved SCUBA galaxies. Clustering is found to increase the confusion noise, compared to the case of purely Poisson fluctuations, by 10-15% for the lowest frequency (i.e. lowest angular resolution) Spitzer and Herschel channels, by 25-35% for the 175 micron ISOPHOT channel, and to dominate in the case of Planck/HFI channels at nu>143GHz. Although our calculations make use of a specific evolutionary model (Granato et al. 2004), the results are strongly constrained by the observed counts and by data on the redshift distribution of SCUBA sources, and therefore are not expected to be heavily model dependent. The main uncertainty arises from the poor observational definition of the source clustering properties. Two models have been used for the latter: a power-law with constant slope and a redshift-independent comoving correlation length,r_0, and the standard theoretical model for clustering evolution in a LambdaCDM universe, with a redshift-dependent bias factor. In both cases, the clustering amplitude has been normalized to yield a unit angular correlation function at theta_0=1-2 arcsec for 850 micron sources fainter than 2 mJy, consistent with the results by Peacock et al. (2000). This normalization yields, for the first model, r_0=8.3$ Mpc/h, and, for the second model, an effective mass of dark matter haloes in which these sources reside of M_halo=1.8*10^{13} M_sun/h. These results are consistent with independent estimates for SCUBA galaxies and for other, likely related, sources.

We have performed 23 GHz VLA observations of 7 compact, luminous infrared galaxies, selected to have evidence of starburst activity. New and published multi-frequency data are combined to obtain the spectral energy distributions of all 7 galaxies from the near-infrared to the radio (at 1.4 GHz). These SEDs are compared with new models, for dust enshrouded galaxies, which account for both starburst and AGN components. In all 7 galaxies the starburst provides the dominant contribution to the infrared luminosity; in 4 sources no contribution from an AGN is required. Although AGN may contribute up to 50 percent of the total far--infrared emission, the starbursts always dominate in the radio. The SEDs of most of our sources are best fit with a very high optical depth of (>=50) at 1 micron. The scatter in the far-infrared/radio correlation, found among luminous IRAS sources, is due mainly to the different evolutionary status of their starburst components. The short time-scale of the star formation process amplifies the delay between the far-infrared and radio emission. This becomes more evident at low radio frequencies (below about 1 GHz) where synchrotron radiation is the dominant process. In the far-infrared (at wavelengths shorter than 100 micron) an additional source of scatter is provided by AGN, where present. AGN may be detected in the near-infrared by the absence of the knee, typical of stellar photospheres. However, near-infrared data alone cannot constrain the level at which AGN contribute because the interpretation of their observed properties, in this wave-band, depends strongly on model parameters.

Granato et al(2004) have elaborated a physically grounded model exploiting the mutual feedback between star-forming spheroidal galaxies and the active nuclei growing in their cores to overcome, within the hierarchical clustering scenario for galaxy formation, one of the main challenges facing such scenario, the fact that massive spheroidal galaxies appear to have formed earlier and faster than predicted by previous models. Adopting the choice by Granato et al for the parameters governing the history of the SF,of chemical abundances and of the gas and dust content of galaxies, we are left with only two parameters affecting the time and mass dependent SED of spheroidal galaxies. After complementing the model with a simple description of evolutionary properties of starburst, normal late-type galaxies and AGNs we have successfully compared the model with a broad variety of observational data, deep K-band, ISO, IRAS, SCUBA, radio counts, the corresponding redshift distributions, the IR background spectrum, and also with data for EROs. We also present detailed predictions for the GOODS and SWIRE surveys with the Spitzer Space Telescope. We find that the GOODS deep survey at 24$\mu$m and the SWIRE surveys at 70 and 160$\mu$m are likely to be severely confusion limited. The GOODS surveys in the IRAC channels are expected to resolve most of the background, to explore the full passive evolution phase of spheroidal galaxies and most of their active star-forming phase, detecting galaxies up to z\simeq 4. A substantial number of high z star-forming spheroidal galaxies should also be detected by the 24\mum SWIRE and GOODS surveys, while the 70 and 160\mum will be particularly useful to study the evolution of such galaxies in the range 1 \lsim z \lsim 2.[abridged]

The mutual feedback between star formation and nuclear activity in large spheroidal galaxies may be a key ingredient to overcome several difficulties plaguing current semi-analytic models for galaxy formation. We discuss some observational implications of the model by Granato et al. (2003) for the co-evolution of galaxies and active nuclei at their centers and stress the potential of the forthcoming surveys of the Sunyaev-Zeldovich effect on arcminute scales, down to $\mu$K levels, to investigate the early galaxy formation phases, difficult to access by other means.

The Sunyaev-Zeldovich (SZ) effect has proven to be an extremely powerful tool to study the physical and evolutionary properties of rich clusters of galaxies. Upcoming SZ experiments, with their much improved sensitivity and angular resolution, will provide unique information also on phases of galaxy evolution characterized by the presence of large amounts of hot proto-galactic gas. We present a preliminary analysis of the SZ signals that can be expected at the collapse of the proto-galaxy when, according to the standard scenario, the gas is heated at its virial temperature, and during episodes of strong energy injections from the active nucleus. The contributions of such signals to excess power on arc-minute scales recently found by CBI and BIMA experiments are briefly discussed.

Charged plasma particles form an anisotropic counter-streaming distribution at the front of a shock. In the near- and ultra-relativistic regimes, Weibel (or two-stream) instability produces near-equipartition, chaotic, small-scale magnetic fields. The fields introduce effective collisions and "thermalize" the plasma particles via pitch-angle scattering. The properties of jitter radiation emitted by accelerated electrons from these small-scale fields are markedly different from synchrotron spectra. Here we theoretical summary and the results of recent numerical 3D PIC plasma kinetic simulations. The relation of the obtained results to the theory of gamma-ray bursts is outlined.

We summarize our detailed, physically grounded, model for the early co-evolution of spheroidal galaxies and of active nuclei at their centers (astro-ph/0307202). Our predictions are excellent agreement with observations for a number of observables which proved to be extremely challenging for all the current semi-analytic models, including the sub-mm counts and the corresponding redshift distributions, and the epoch-dependent K-band luminosity function of spheroidal galaxies. Also, the black hole mass function and the relationship between the black hole mass and the velocity dispersion in the galaxy are nicely reproduced.

We present a physically motivated model for the early co-evolution of massive spheroidal galaxies and active nuclei at their centers. Within dark matter halos, forming at the rate predicted by the canonical hierarchical clustering scenario, the gas evolution is controlled by gravity, radiative cooling, and heating by feedback from supernovae and from the growing active nucleus. Supernova heating is increasingly effective with decreasing binding energy in slowing down the star formation and in driving gas outflows. The more massive proto-galaxies virializing at earlier times are thus the sites of the faster star-formation. The correspondingly higher radiation drag fastens the angular momentum loss by the gas, resulting in a larger accretion rate onto the central black-hole. In turn, the kinetic energy carried by outflows driven by active nuclei can unbind the residual gas, thus halting both the star formation and the black-hole growth, in a time again shorter for larger halos. For the most massive galaxies the gas unbinding time is short enough for the bulk of the star-formation to be completed before type Ia supernovae can substantially increase the $Fe$ abundance of the interstellar medium, thus accounting for the $\alpha$-enhancement seen in the largest galaxies. The feedback from supernovae and from the active nucleus also determines the relationship between the black-hole mass and the mass, or the velocity dispersion, of the host galaxy, as well as the black-hole mass function. These and other model predictions are in excellent agreement with observations for a number of observables which proved to be extremely challenging for all the current semi-analytic models, including the sub-mm counts and the corresponding redshift distributions, and the epoch-dependent K-band luminosity function of spheroidal galaxies.

We present the first three-dimensional fully kinetic electromagnetic relativistic particle-in-cell simulations of the collision of two interpenetrating plasma shells. The highly accurate plasma-kinetic "particle-in-cell" (with the total of $10^8$ particles) parallel code OSIRIS has been used. Our simulations show: (i) the generation of long-lived near-equipartition (electro)magnetic fields, (ii) non-thermal particle acceleration, and (iii) short-scale to long-scale magnetic field evolution, in the collision region. Our results provide new insights into the magnetic field generation and particle acceleration in relativistic and sub-relativistic colliding streams of particles, which are present in gamma-ray bursters, supernova remnants, relativistic jets, pulsar winds, etc..

We present a model for nebular emission in star forming galaxies, which takes into account the effects of dust reprocessing. The nebular emissions have been computed with CLOUDY and then included into GRASIL, our spectrophotometric code specifically developed for dusty galaxies. The interface between nebular emission and population synthesis is based on a set of pre-computed HII region emission models covering a wide range of physical quantities. Concerning the extinction properties of normal star forming galaxies, we are able to interpret the observed lack of correlation between the attenuation measured at Halpha and in the UV band as a consequence of age selective extinction. We also find that, for these galaxies with modest SFR, the ratio FIR/UV provides the best constraints on the UV attenuation. Our model also allows to deal with different SFR estimators in a consistent way, from the UV to radio wavelengths, and to discuss the uncertainties arising from the different physical conditions encountered in star forming galaxies. We provide our best estimates of SFR/luminosity calibrations, together with their expected range of variation. It results that SFR derived through Halpha, even when corrected for extinction using the Balmer decrement, is affected by important uncertainties due to age selective extinction. Another remarkable result is that SFR from UV luminosity corrected by means of the ratio FIR/UV has a small uncertainty. Finally, our model provides a calibration of SFR from radio luminosity; we are also able to reproduce the observed FIR/radio ratio.

The trend towards ever larger telescopes and more advanced adaptive optics systems is driving the need for deformable mirrors with a large number of low cost actuators. Liquid mirrors have long been recognized a potential low cost alternative to conventional solid mirrors. By using a water or oil based ferrofluid we are able to benefit from a stronger magnetic response than is found in magnetic liquid metal amalgams and avoid the difficulty of passing a uniform current through a liquid. Depositing a thin silver colloid known as a metal liquid-like film (MELLF) on the ferrofluid surface solves the problem of low reflectivity of pure ferrofluids. This combination provides a liquid optical surface that can be precisely shaped in a magnetic field. We present experimental results obtained with a prototype deformable liquid mirror based on this combination.

In view of the extensive evidence of tight inter-relationships between spheroidal galaxies (and galactic bulges) with massive black holes hosted at their centers, a consistent model must deal jointly with the evolution of the two components. We describe one such model, which successfully accounts for the local luminosity function of spheroidal galaxies, for their photometric and chemical properties, for deep galaxy counts in different wavebands, including those in the (sub)-mm region which proved to be critical for current semi-analytic models stemming from the standard hierarchical clustering picture, for clustering properties of SCUBA galaxies, of EROs, and of LBGs, as well as for the local mass function of massive black holes and for quasar evolution. Predictions that can be tested by surveys carried out by SIRTF are presented.

Our multi-wavelength model GRASIL for the SED of galaxies is described, in particular the recent extension to the radio and X-ray range. With our model we can study different aspects of galaxy evolution by exploiting all available spectral observations, where different emission components dominate.

We revisit the nature of the FIR/Radio correlation by means of the most recent models for star forming galaxies. We model the IR emission with our population synthesis code, GRASIL (Silva et al. 1998). As for the radio emission, we revisit the simple model of Condon & Yin (1990). We find that a tightFIR/Radio correlation is natural when the synchrotron mechanism dominates over the inverse Compton, and the electrons cooling time is shorter than the fading time of the supernova rate. Observations indicate that both these conditions are met in star forming galaxies. However since the radio non thermal emission is delayed, deviations are expected both in the early phases of a starburst, when the radio thermal component dominates, and in the post-starburst phase, when the bulk of the NT component originates from less massive stars. This delay allows the analysis of obscured starbursts with a time resolution of a few tens of Myrs, unreachable with other star formation indicators. We suggest to complement the analysis of the deviations from the FIR/Radio correlation with the radio slope to obtain characteristic parameters of the burst. The analysis of a sample of compact ULIRGs shows that they are intense but transient starbursts, to which one should not apply usual SF indicators devised for constant SF rates. We also discuss the possibility of using the q- radio slope diagram to asses the presence of obscured AGN. A firm prediction of the models is an apparent radio excess during the post-starburst phase, which seems to be typical of a class of star forming galaxies in rich cluster cores. We discuss how deviations from the correlation, due to the evolutionary status of the starburst, affect the technique of photometric redshift determination widely used for high-z sources.

Surveys at mm wavelengths emphasize rare classes of extragalactic radio sources characterized by spectra keeping flat or inverted up to high frequencies (such as blazars, GPS sources, advection dominated sources, or even gamma-ray burst afterglows), which are are difficult to single out at lower frequencies where the counts are dominated by far more numerous populations. Such surveys may lead to the discovery of populations with extreme properties, such as sources of the general GPS or blazar type but peaking at mm wavelengths, whose existence is hinted by currently available data. At relatively lower flux levels mm counts are dominated by high-redshift dusty galaxies, such as those detected by SCUBA and MAMBO surveys. These data have profound implications for our understanding of the formation and early evolution of galaxies. High sensitivity surveys at cm/mm wavelengths may also detect galaxy-scale Sunyaev-Zeldovich effects due to gas heated by central active nuclei.

Direct and indirect pieces of observational evidence point to a strong connection between high-redshift quasars and their host galaxies. In the framework of a model where the shining of the quasar is the episode that stops the formation of the galactic spheroid inside a virialized halo, it has been proven possible to explain the submillimetre source counts together with their related statistics and the local luminosity function of spheroidal galaxies. The time delay between the virialization and the quasar manifestation required to fit the counts is short and incresing with decresing the host galaxy mass. In this paper we compute the detailed chemical evolution of gas and stars inside virialized haloes in the framework of the same model, taking into account the combined effects of cooling and stellar feedback. Under the assumption of negligible angular momentum, we are able to reproduce the main observed chemical properties of local ellipticals. In particular, by using the same duration of the bursts which are required in order to fit the submillimetre source counts, we recover the observed increase of the Mg/Fe ratio with galactic mass. Since for the most massive objects the assumed duration of the burst is Tburst < 0.6 Gyr, we end up with a picture for elliptical galaxy formation in which massive spheroids complete their assembly at early times, thus resembling a monolithic collapse, whereas smaller galaxies are allowed for a more prolonged star formation, thus allowing for a more complicated evolutionary history.

We present results from a deep mid-infrared survey of the Hubble Deep Field South (HDF-S) region performed at 7 and 15 micron with the CAM instrument on board the Infrared Space Observatory (ISO). The final map in each band was constructed by the coaddition of four independent rasters, registered using bright sources securely detected in all rasters, with the absolute astrometry being defined by a radio source detected at both 7 and 15 micron. We sought detections of bright sources in a circular region of radius 2.5 arcmin at the centre of each map, in a manner that simulations indicated would produce highly reliable and complete source catalogues using simple selection criteria. Merging source lists in the two bands yielded a catalogue of 35 distinct sources, which we calibrated photometrically using photospheric models of late-type stars detected in our data. We present extragalactic source count results in both bands, and discuss the constraints they impose on models of galaxy evolution models, given the volume of space sampled by this galaxy population.

We present results from a deep mid-IR survey of the Hubble Deep Field South (HDF-S) region performed at 7 and 15um with the CAM instrument on board ISO. We found reliable optical/near-IR associations for 32 of the 35 sources detected in this field by Oliver et al. (2002, Paper I): eight of them were identified as stars, one is definitely an AGN, a second seems likely to be an AGN, too, while the remaining 22 appear to be normal spiral or starburst galaxies. Using model spectral energy distributions (SEDs) of similar galaxies, we compare methods for estimating the star formation rates (SFRs) in these objects, finding that an estimator based on integrated (3-1000um) IR luminosity reproduces the model SFRs best. Applying this estimator to model fits to the SEDs of our 22 spiral and starburst galaxies, we find that they are forming stars at rates of ~1-100 M_sol/yr, with a median value of ~40M_sol/yr, assuming an Einstein - de Sitter universe with a Hubble constant of 50 km/s/Mpc, and star formation taking place according to a Salpeter (1955) IMF across the mass range 0.1-100M_sol. We split the redshift range 0.0<z<0.6 into two equal-volume bins to compute raw estimates of the star formation rate density contributed by these sources, assuming the same cosmology and IMF as above and computing errors based on estimated uncertainties in the SFRs of individual galaxies. We compare these results with other estimates of the SFR density made with the same assumptions, showing them to be consistent with the results of Flores et al. (1999) from their ISO survey of the CFRS 1415+52 field. However, the relatively small volume of our survey means that our SFR density estimates suffer from a large sampling variance, implying that our results, by themselves, do not place tight constraints on the global mean SFR density.

In view of the extensive evidence of a tight inter-relationship between spheroidal galaxies (and galactic bulges) and massive black holes hosted at their centers, a consistent model must deal jointly with the evolution of the two components. We describe one viable model, which successfully accounts for the local luminosity function of spheroidal galaxies, their photometric and chemical properties, deep galaxy counts in different wavebands, including those in the (sub)-mm region which proved to be critical for current semi-analytic models stemming from the standard hierarchical clustering picture, clustering properties of SCUBA galaxies, of EROs, and of LBGs, as well as for the local mass function of massive black holes and for quasar evolution. Predictions that can be tested by surveys carried out by SIRTF are presented.

In view of the extensive evidence of tight inter-relationships between spheroidal galaxies (and galactic bulges) with massive black holes hosted at their centers, a consistent model must deal jointly with the evolution of the two components. We describe one such model, which successfully accounts for the local luminosity function of spheroidal galaxies, for their photometric and chemical properties, for deep galaxy counts in different wavebands, including those in the (sub)-mm region which proved to be critical for current semi-analytic models stemming from the standard hierarchical clustering picture, for clustering properties of SCUBA galaxies, of EROs, and of LBGs, as well as for the local mass function of massive black holes and for quasar evolution. Predictions that can be tested by surveys carried out by SIRTF are presented.

We show that the predictions of the recent physical model by Granato et al. (2001) are fully consistent with the available estimates of the clustering properties of galaxies at (sub)-mm wavelengths. The model is then exploited to compute the expected power spectrum of fluctuations due to clustering at the frequencies of the High Frequency Instrument (HFI) of ESA's {\sc Planck} satellite. The clustering signal is found to be detectable in regions of low interstellar dust emission. A further distinctive prediction of the adopted model is a remarkably high fraction of gravitationally lensed sources at bright mm/sub-mm fluxes. In fact, according to the adopted model, only thanks to gravitational lensing, a significant number of high-$z$ forming spheroidal galaxies will be detectable by large area, shallow surveys at mm/sub-mm wavelengths, such as those carried out by P{\sc lanck}/HFI. Allowing for other source populations, we find that the fraction of gravitationally lensed sources at bright mm/sub-mm fluxes is expected to be up to $\simeq 40%$.

We present our model to interpret the SED of galaxies. The model for the UV to sub-mm SED is already well established (Silva et al 1998). We remind here its main features and show some applications. Recently we have extended the model to the radio range (Bressan et al 2001), and we have started to include the X-ray emission from the stellar component.

Star forming galaxies exhibit a variety of physical conditions, from quiescent normal spirals to the most powerful dusty starbursts. In order to study these complex systems, we need a suitable tool to analyze the information coming from observations at all wavelengths. We present a new spectro-photometric model which considers in a consistent way starlight as reprocessed by gas and dust. We discuss preliminary results to interpret some observed properties of VLIRGs.

We summarize recent advances in the determination of the cosmic history of star formation and other properties of high-z galaxies, and the relevance of this information in our understanding of the formation of structures. We emphasize the importance of dust reprocessing in the high--z universe, as demonstrated in particular by IR and sub-mm data. This demand a panchromatic approach to observations and suitable modelling tools. We spend also some words on expectations from future instruments.

The effects of the highly damped modes in the energy and reaction rates in a plasma are discussed. These modes, with wavenumbers $k \gg k_{D}$, even being only weakly excited, with less than $k_{B}T$ per mode, make a significant contribution to the energy and screening in a plasma. When the de Broglie wavelength is much less than the distance of closest approach of thermal electrons, a classical analysis of the plasma can

We examine the tight correlation observed between the far infrared (FIR) and radio emission in normal and star burst galaxies. We show that significant deviations from the average relation are to be expected in young star burst or post star burst galaxies, due to the different fading times of FIR and Radio power.

Direct and indirect observational evidence leads to the conclusion that high redshift QSOs did shine in the core of early type proto-galaxies during their main episode of star formation. Exploting this fact, we derive the rate of formation of this kind of stellar systems at high redshift by using the QSO Luminosity Function. The elemental proportions in elliptical galaxies, the descendents of the QSO hosts, suggest that the star formation was more rapid in more massive objects. We show that this is expected to occur in Dark Matter haloes, when the processes of cooling and heating is considered. This is also confirmed by comparing the observed sub-mm counts to those derived by coupling the formation rate and the star formation rate of the spheroidal galaxies with a detailed model for their SED evolution. In this scenario SCUBA galaxies and Lyman Break Galaxies are early type proto-galaxies forming the bulk of their stars before the onset of QSO activity.

We analyze a morphologically-selected complete sample of 52 late-type (spiral and irregular) galaxies in the Hubble Deep Field North with total K-magnitudes brighter than K=20.47. This sample exploits in particular the ultimate imaging quality achieved by HST in this field, allowing us to clearly disentangle the galaxy morphologies, based on accurate profiles of the surface brightness distributions. Our purpose was to investigate systematic differences between the two classes, as for colours, redshift distributions and ages of the dominant stellar populations. This sample appears to miss significantly galaxies above z=1.4 (in a similar way as an early-type galaxy sample previously studied by us), a fact which may be explained as a global decline of the underlying mass function for galaxies at these high redshifts. Differences between early and late-types are apparent -particularly in the colour distributions and the evolutionary star-formation (SF) rates per unit volume-, although the complication in spectro-photometric modelling introduced by dust-extinction in the gas-rich systems prevents us to reach conclusive results on the single sources. However, we find that an integrated quantity like the comoving star-formation rate density (SFR(t)) as a function of redshift is much less affected by these uncertainties: by combining this with the previously studied early-type galaxy sample, we find a shallower dependence of SFR(t) on z between z=0.2 and z=1.5 than found by Lilly et al. (1995).

Recently (Granato, Lacey, Silva et al. 2000, astro-ph/0001308) we have combined our spectrophotometric galaxy evolution code which includes dust reprocessing (GRASIL, Silva et al. 1998) with semi-analytical galaxy formation models (GALFORM, Cole et al. 1999). One of the most characteristic features of the former is that the dust is divided in two main phases: molecular cloud complexes, where stars are assumed to be born, and the diffuse interstellar medium. As a consequence, stellar populations of different ages have different geometrical relationships with the two phases, which is essential in understanding several observed properties of galaxies, in particular those undergoing major episodes of star formation at any redshift. Indeed, our merged GRASIL+GALFORM model reproduces fairly well the SEDs of normal spirals and starbursts from the far-UV to the sub-mm and their internal extinction properties. In particular in the model the observed starburst attenuation law (Calzetti 1999) is accounted for as an effect of geometry of stars and dust, and has nothing to do with the optical properties of dust grains.

We present a new evolutionary model for the far-UV to sub-mm properties of the galaxy population. This combines a semi-analytic galaxy formation model based on hierarchical clustering (GALFORM) with a spectro-photometric code which includes dust reprocessing (GRASIL). The former provides the star formation and metal enrichment histories, together with the gas mass and various geometrical parameters, for a representative sample of galaxies formed in different density environments. These quantities allow us to model the SEDs of galaxies, taking into account stellar emission and also dust extinction and re-emission. Two phases are considered for the dust: molecular cloud complexes, where stars are assumed to be born, and the diffuse interstellar medium. The model includes both galaxies forming stars quiescently in disks, and starbursts triggered by galaxy mergers. We test our models against the observed spectro-photometric properties of galaxies in the local Universe. The models reproduce fairly well the SEDs of normal spirals and starbursts, and their internal extinction properties. The starbursts follow the observed relationship between the FIR to UV luminosity ratio and the slope of the UV continuum. They also reproduce the observed starburst attenuation law (Calzetti et al 99). This result is remarkable, because we use a dust mixture which reproduces the Milky Way extinction law. It suggests that the observed attenuation law is related to the geometry of the stars and dust. We compute galaxy luminosity functions over our wide range of wavelengths, which turn out to be in good agreement with observational data. The UV continuum turns out to be a poor star formation indicator for our models, whilst the infrared luminosity is much more reliable.

We have explored, in the general framework of the ``young source'' scenario, evolutionary models for GHz Peaked Spectrum (GPS) galaxies and quasars which reproduce the observed counts, redshift and peak frequency distributions of currently available samples. Substantially different cosmological evolution properties are found for the two populations: the quasar luminosity function must evolve strongly up to $z\sim 1$, while the data on galaxies may be consistent with no evolution. The models show that GPS sources (mostly quasars) may comprise quite a significant fraction of bright ($S> 1 $Jy) radio sources at $\nu \geq 30 $GHz if the intrinsic distribution of peak frequencies extends up to $\sim 1000 $GHz. In any case, however, their fraction decreases rapidly with decreasing flux and their contribution to small scale fluctuations in the frequency range covered by the forthcoming space missions MAP and Planck Surveyor is expected to be minor.

We present the basic features and preliminary results of the interface between our spectro-photometric model GRASIL (that calculates galactic SED from the UV to the sub-mm with a detailed computation of dust extinction and thermal reemission) with the semi-analytical galaxy formation model GALFORM (that computes galaxy formation and evolution in the hierarchical scenario, providing the star formation history as an input to our model). With these two models we are able to synthetize simulated samples of a few thousands galaxies suited for statistical studies of galaxy properties to investigate on galaxy formation and evolution. We find good agreement with the available data of SED and luminosity functions.

We summarize our modelling of galaxy photometric evolution (GRASIL code). By including the effects of dust grains and PAHs molecules in a two phases clumpy medium, where clumps are associated to star forming regions, we reproduce the observed UV to radio SEDs of galaxies with star formation rates from zero to several hundreds $M_\odot/yr$. GRASIL is a powerful tool to investigate the star formation, the initial mass function and the supernovae rate in nearby starbursts and normal galaxies, as well as to predict the evolution of luminosity functions of different types of galaxies at wavelengths covering six decades. It may be interfaced with any device providing the star formation and metallicity histories of a galaxy. As an application, we have investigated the properties of early-type galaxies in the HDF, tracking the contribution of this population to the cosmic star-formation history, which has a broad peak between z=1.5 and 4. To explain the absence of objects at z $\gsim$ 1.3, we suggest a sequence of dust-enshrouded merging-driven starbursts in the first few Gyrs of galaxies lifetime. We are at present working on a complementary sample of late type objects selected in a similar way.

The derivation of the ponderomotive force due to neutrinos by Hardy and Melrose, Phys.Rev.D v.54, p.6491 (1996) is based on a flawed analysis of the ponderomotive force concept. Their conclusions also contain an erroneous physical assumption related to the neutrino emission in supernovae. A correct analysis shows the importance of the ponderomotive force due to neutrinos in type II supernovae explosions.

We present models of photometric evolution of galaxies in which the effects of a dusty interstellar medium have been included with particular care. A chemical evolution code follows the star formation rate, the gas fraction and the metallicity, basic ingredients for the stellar population synthesis. The latter is performed with a grid of integrated spectra of simple stellar populations (SSP) of different ages and metallicities, in which the effects of dusty envelopes around asymptotic giant branch (AGB) stars are included. The residual fraction of gas in the galaxy is divided into two phases: the star forming molecular clouds and the cirrus. The relative amount is a model parameter. The molecular gas is sub--divided into clouds of given mass and radius: it is supposed that each SSP is born within the cloud and progressively escapes it. The emitted spectrum of the star forming molecular clouds is computed with a radiative transfer code. The cirrus emission 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 PAHs are included. The model is compared and calibrated with available data of normal and starburst galaxies in the local universe, in particular new broad--band and spectroscopic ISO observations. It will be a powerful tool to investigate the star formation, the IMF, supernovae rate in nearby starbursts and normal galaxies, as well as to predict the evolution of luminosity functions of different types of galaxies at wavelengths covering four decades.

We have investigated the properties of a complete K-band selected sample of 35 elliptical and S0 galaxies brighter than K=20.15 in the Hubble Deep Field, as representative of the field galaxy population. This sample has been derived from deep K-band image by the KPNO-IRIM camera, by applying a rigorous morphological classification scheme based on quantitative analyses of the surface brightness profiles. The broad-band spectra of the sample galaxies allow us to date their dominant stellar populations. The majority of bright early-types in this field are found at redshifts z<1.3 to display colors indicative of a fairly wide range of ages (typically 1.5 to 3 Gyrs). We find that the major episodes of star-formation building up typical M_star galaxies have taken place during a wide redshift interval 1<z<4 for q_0=0.5, which becomes 1<z<3 for q_0=0.15. Our estimated galactic masses, for a Salpeter IMF, are found in the range from a few billion to a few 10^{11} solar masses already at z=1. So the bright end of the E/S0 population is mostly in place by that cosmic epoch, with space densities, masses and luminosities consistent with those of the local field E/S0 population. What distinguishes the present sample is a remarkable absence of objects at z>1.3, which should be detectable during the luminous star-formation phase expected to happen at these redshifts. Obvious solutions are a) that the merging events imply perturbed morphologies which prevent selecting them by our morphological classification filter, or b) that a dust-polluted ISM obscures the (either continuous or episodic) events of star-formation. We conclude that the likely solution is a combination thereof, i.e. a set of dust-enshrouded merging-driven starbursts occurring during the first few Gyrs of the galaxy's lifetime.

We have investigated the spectro-photometric properties of the Asymptotic Giant Branch (AGB) stars and their contribution to the integrated infrared emission in simple stellar populations (SSP). Adopting analytical relations describing the evolution of these stars in the HR diagram and empirical relations for the mass-loss rate and the wind terminal velocity, we were able to model the effects of the dusty envelope around these stars, with a minimal number of parameters. We computed isochrones at different age and initial metal content. We compare our models with existing infrared colors of M giants and Mira stars and with IRAS PSC data. Contrary to previous models, in the new isochrones the mass-loss rate, which establishes the duration of the AGB phase, also determines the spectral properties of the stars. The contribution of these stars to the integrated light of the population is thus obtained in a consistent way. We find that the emission in the mid infrared is about one order of magnitude larger when dust is taken into account in an intermediate age population, irrespective of the particular mixture adopted. The dependence of the integrated colors on the metallicity and age is discussed, with particular emphasis on the problem of age-metallicity degeneracy. We show that, contrary to the case of optical or near infrared colors, the adoption of a suitable pass-band in the mid infrared allows a fair separation of the two effects. We suggest intermediate redshift elliptical galaxies as possible targets of this method of solving the age-metallicity dilemma. The new SSP models constitute a first step in a more extended study aimed at modelling the spectral properties of the galaxies from the ultraviolet to the far infrared.

Our present understanding of the extragalactic source counts and background radiation at infrared and sub-mm wavelengths is reviewed. Available count data are used to constrain evolutionary models of galaxies and Active Nuclei. The CIRB, on the other hand, provides crucial information on the integrated past IR emissivity, including sources so faint to be never accessible. The near-IR and sub-mm are suited to search for the CIRB. These cosmological windows are ideal to detect redshifted photons from the two most prominent broad emission features in galaxy spectra: the stellar photospheric 1 micron peak and the one at 100 micron due to dust re-radiation. The recently claimed detection of an isotropic diffuse component at 100-200 microns would support the evidence for strong cosmic evolution from $z=0$ to $z\sim 2$ of faint gas-rich late-type galaxies, as inferred from direct long-$\lambda$ counts. Furthermore, an equally intense CIRB flux measured in the wavelength range 200 to 500 $\mu m$ may be in support of models envisaging a dust-enshrouded phase during formation of spheroidal galaxies (e.g. Franceschini et al., 1994): the background spectrum at such long wavelengths would indicate an high redshift ($z>2-3$) for this active phase of star formation, implying a relevant constraint on structure formation scenarios. We argue that, given the difficulties for sky surveys in this spectral domain, only sporadic tests of these ideas will be achieved in the next several years. Only a mission like FIRST, in combination with large mm arrays and optical telescopes, will allow to thoroughly investigate, via photometric and spectroscopic surveys, these IR-bright early phases in the evolution of galaxies and AGNs.

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

We have investigated the properties of a bright K-band selected sample of early-type galaxies in the Hubble Deep Field, as representative of the field galaxy population. This dataset is unique as for the morphological information on faint high-z sources, and for complete photometric and spectroscopic coverage. The majority of bright early-type galaxies in this field are found at redshifts $z \leq 1.3$ to share common properties with those of high-z cluster samples, as for the age and mass of the dominant stellar population -- which are found to be as old as 3-5 Gyr and as massive as $10^{11} M_\odot$ already at $z\simeq 1$. There is no evidence of a steepening of the mass function from $z$=0 to $z$=1, as inferred by some authors from analyses of optically-selected samples and favoured by hierarchical clustering models forming most of the E/S0s at $z<1$. What distinguishes this sample is a remarkable absence of objects at $z>1.3$, which would be expected as clearly detectable above the flux limits, given the aged properties of the lower redshift counterparts. So, something hide them at high redshifts. Merging could be an explanation, but it would require that already formed and old stellar systems would assemble on short timescales of 1 Gyr or less. Only a negligible fraction of the galactic mass in young stars has to be added, a rather contrived situation from the dynamical point of view. An alternative interpretation could be that a dust-polluted ISM obscures the first 3 to 4 Gyr of the evolution of these massive systems, after which a galactic wind and gas consumption makes them transparent. The presence of dust would have relevant (and testable) implications for the evaluation of the global energetics from galaxy formation and for the visibility of the early phases by current and future infrared