What is the expected infrared output of elliptical galaxies? Here we report the latest findings obtained in this high time resolution (~10 years) and high spatial resolution (2.5 parsec at center) study. We add a set of grain physics to the MACER code, including (a) dust grains made in passive stellar evolution; (b) dust grain growth due to collision and sticking; (c) grain destruction due to thermal sputtering; (d) dust cooling of hot gas via inelastic collisions; and (e) radiation pressure on dust grains. The code improvements enable us to analyze the effects of dust on metal depletion and AGN obscuration, and also to assess the infrared output of the modeled galaxies. We simulate a representative massive elliptical galaxy of a central stellar velocity dispersion ~ 260 km/s and modest rotation. We find that: (1) the circumnuclear disk (of a size <~ 1 kpc) is dusty in its outer region where most of the metals are in dust grains, while in the inner disk most of the dust grains are destroyed by the AGN irradiation; (2) the dusty disk is optically thick to both the starlight within the disk and the radiation from the central AGN. Thus the AGN will be obscured behind the disk, and the latter is of a covering factor ~ 0.2; (3) the dust infrared emission is mainly due to the AGN irradiation. The median infrared luminosity is ~ 2e44 erg/s, and it can reach >~ 1e46 erg/s during outbursts; (4) the duty cycles of the AGN activities, star formation, and the dust infrared luminosity roughly match observations, e.g., in most of its lifetime, the simulated galaxy is a stereotypical "quiescent" elliptical galaxy with L_{IR} ~ 1e11*L_{solar}, while the star formation rate can exceed 250 M_{solar}/yr during central outbursts.

A hot plasma is the dominant phase of the interstellar medium of early-type galaxies. Its origin can reside in stellar mass losses, residual gas from the formation epoch, and accretion from outside of the galaxies. Its evolution is linked to the dynamical structure of the host galaxy, to the supernova and AGN feedback, and to (late-epoch) star formation, in a way that has yet to be fully understood. Important clues about the origin and evolution of the hot gas come from the abundances of heavy metals, that have been studied with increasing detail with XMM-Newton and Chandra. We present recent high resolution hydrodynamical simulations of the hot gas evolution that include the above processes, and where several chemical species, originating in AGB stars and supernovae of type Ia and II, have also been considered. The high resolution, of few parsecs in the central galactic region, allows us to track the metal enrichment, transportation and dilution throughout the galaxy. The comparison of model results with observed abundances reveals a good agreement for the region enriched by the AGN wind, but also discrepancies for the diffuse hot gas; the latter indicate the need for a revision of standard assumptions, and/or the importance of neglected effects as those due to the dust, and/or residual uncertainties in deriving abundances from the X-ray spectra.

The Smallest Scale of Hierarchy (SSH) survey is an ongoing strategic large program at the Large Binocular Telescope, aimed at the detection of faint stellar streams and satellites around 45 late-type dwarf galaxies located in the Local Universe within $\simeq$10 Mpc. SSH exploits the wide-field, deep photometry provided by the Large Binocular Cameras in the two wide filters $g$ and $r$. This paper describes the survey, its goals, and the observational and data reduction strategies. We present preliminary scientific results for five representative cases (UGC 12613, NGC 2366, UGC 685, NGC 5477 and UGC 4426) covering the whole distance range spanned by the SSH targets. We reach a surface brightness limit as faint as $\mu(r)\sim$ 31 mag arcsec$^{-2}$ both for targets closer than 4$-$5 Mpc, which are resolved into individual stars, and for more distant targets through the diffuse light. Our analysis reveals the presence of extended low surface brightness stellar envelopes around the dwarfs, reaching farther out than what traced by the integrated light, and as far out as, or even beyond, the observed HI disk. Stellar streams, arcs, and peculiar features are detected in some cases, indicating possible perturbation, accretion, or merging events. We also report on the discovery of an extreme case of Ultra Diffuse Galaxy ($\mu_g(0)=27.9$~mag/arcsec$^2$) in the background of one of our targets, to illustrate the power of the survey in revealing extremely low surface brightness systems.

We present a significantly improved version of our numerical code JASMINE, that can now solve the Jeans equations for axisymmetric models of stellar systems, composed of an arbitrary number of stellar populations, a Dark Matter halo, and a central Black Hole. The stellar components can have different structural (density profile, flattening, mass, scale length), dynamical (rotational support, velocity dispersion anisotropy), and population (age, metallicity, Initial Mass Function, mass-to-light ratio) properties. These models, when combined with observations, will allow to investigate important issues, such as quantifying the systematic effects of IMF variations, of mass-to-light ratio gradients, and of different stellar kinematic components (e.g. counter rotating disks, kinematically decoupled cores) on luminosity-weighted properties. The developed analytical and numerical framework aims at modeling Early-Type Galaxies, but it can also be applied to dwarf Spheroidal galaxies and Globular Clusters.

The dynamical properties of spherically symmetric galaxy models, where a Jaffe (1983) stellar density profile is embedded in a total mass density decreasing as $r^{-3}$ at large radii, are presented. The orbital structure of the stellar component is described by the Osipkov--Merritt anisotropy; the dark matter halo is isotropic, and a black hole is added at the center of the galaxy. First, the conditions for a nowhere negative and monotonically decreasing dark matter halo density profile are derived; this profile can be made asymptotically coincident with a NFW profile at the center and at large radii. Then the minimum value of the anisotropy radius for phase-space consistency is derived as a function of the galaxy parameters. The Jeans equations for the stellar component are solved analytically; the projected velocity dispersion at the center and at large radii is also obtained, for generic values of the anisotropy radius. Finally, analytical expressions for the terms entering the Virial Theorem are derived, and the fiducial anisotropy limit required to prevent the onset of Radial Orbit Instability is determined as a function of the galaxy parameters. The presented models, built following an approach already adopted in our previous works, can be a useful starting point for a more advanced modeling of the dynamics of elliptical galaxies, and can be easily implemented in numerical simulations requiring a realistic dynamical model of a galaxy.

Depending on the value of the accretion rate, black hole accretion is divided into cold and hot modes. The two modes have distinctly different physics and correspond to two feedback modes. Most previous feedback works either focus only on one mode, or the accretion physics is not always properly adopted. Here, by performing hydrodynamical numerical simulations of AGN feedback in an elliptical galaxy, we show that including both is important, and gives different results from including just one or the other. We specifically focus on the wind and radiation (but neglecting the jet) feedback in the hot mode, to explore whether this particular mode of feedback can play any important feedback role. For this aim, we have run two test models. In one model, we always adopt the cold mode no matter what value the accretion rate is; in another model, we turn off the AGN once it enters into the hot mode. We have calculated the AGN light curves, black hole growth, star formation, and AGN energy duty-cycle; and compared the results with the model in which the two modes are correctly included. Important differences are found. For example, if we were to adopt only the cold mode, the total mass of newly formed stars would become two orders of magnitude smaller, and the fraction of energy ejected within the high accretion regime (i.e., $L_{\rm BH} > 2\% L_{\rm Edd}$) would be too small to be consistent with observations. We have also investigated the respective roles of wind and radiation in the hot mode.

We expand our previous analytical and numerical studies of the family of S\'ersic models, routinely used to describe early-type galaxies and the bulges of spiral galaxies. In particular, we focus on the total energy budget, an important dynamical property that has not been discussed in detail in previous works. We use two different methods to calculate the total energy for the S\'ersic model family, resulting in two independent expressions that can be used along the entire sequence of S\'ersic models. We use these expressions to investigate whether the Spitzer concentration index is a reliable measure for the intrinsic 3D concentration of galaxies, and we conclude that it is not a very useful measure for the central concentration. The popular Third Galaxy Concentration index, on the other hand, is shown to be a reliable measure for the intrinsic 3D concentration, even though it is based on the surface brightness distribution and not on the intrinsic 3D density.

We investigate the dynamics of supermassive black holes (SMBHs) in galactic cores by means of a semi-analytic model based on the Langevin equation, including dynamical friction and stochastic noise accounting for the gravitational interactions with stars. The model is validated against direct $N$-body simulations of intermediate-mass black holes in stellar clusters where a realistic number of particles is accessible. For the galactic case, we find that the SMBH experiences a Brownian-like motion with a typical displacement from the geometric center of the Galaxy of a few parsecs, for system parameters compatible with M87. \keywords{stellar dynamics, black hole physics, methods: n-body simulations, methods: statistical.

We present new Hubble Space Telescope (HST) imaging of a stream-like system associated with the dwarf galaxy DDO 68, located in the Lynx-Cancer Void at a distance of D$\sim$12.65 Mpc from us. The stream, previously identified in deep Large Binocular Telescope images as a diffuse low surface brightness structure, is resolved into individual stars in the F606W (broad V) and F814W ($\sim$I) images acquired with the Wide Field Camera 3. The resulting V, I color-magnitude diagram (CMD) of the resolved stars is dominated by old (age$\gtrsim$1-2 Gyr) red giant branch (RGB) stars. From the observed RGB tip, we conclude that the stream is at the same distance as DDO 68, confirming the physical association with it. A synthetic CMD analysis indicates that the large majority of the star formation activity in the stream occurred at epochs earlier than $\sim$1 Gyr ago, and that the star formation at epochs more recent than $\sim$500 Myr ago is compatible with zero. The total stellar mass of the stream is $\sim10^{6} M_{\odot}$, about 1/100 of that of DDO~68. This is a striking example of hierarchical merging in action at the dwarf galaxy scales.

The S\'ersic or $R^{1/m}$ model has become the de facto standard model to describe the surface brightness profiles of early-type galaxies and the bulges of spiral galaxies. The photometric, intrinsic, and dynamical properties of this model have been investigated, but mainly for fairly large S\'ersic indices $m$. For small values of $m$, appropriate for low-mass and dwarf ellipticals, a detailed investigation of these properties is still lacking. In this study, we used a combination of numerical and analytical techniques to investigate the S\'ersic model over the entire range of S\'ersic parameters, focussing on the small $m$ regime, where a number of interesting and surprising properties are found. For all values $m<1$, the model is characterised by a finite central luminosity density, and for $m<\tfrac12$, even a central depression in the luminosity density profile. This behaviour translates to the dynamical properties: we show that all S\'ersic models with $m \geqslant\tfrac12$ can be supported by an isotropic velocity dispersion tensor, and that these isotropic models are stable to both radial and non-radial perturbations. The models with $m < \tfrac12$, on the other hand, cannot be supported by an isotropic velocity dispersion tensor.

We add a suite of chemical abundances to the MACER (Massive AGN Controlled Ellipticals Resolved) 2D code, by solving 12 additional continuity equations for H, He, C, N, O, Ne, Mg, Si, S, Ca, Fe and Ni respectively with sources from AGB stars and supernovae of type Ia and II with metal yields based on standard stellar physics. New stars, formed in Toomre unstable circumnuclear disks (of a size $\lesssim150$ parsec), are assumed to have a top-heavy initial mass function with a power index of 1.65. The metal dilution effects due to cosmic accretion are also included. With the high resolution of few parsecs in central regions, resolved black hole accretion and AGN feedback, we can track the metal enrichment, transportation and dilution throughout the modeled massive elliptical galaxy of velocity dispersion $\sim280$ km/s. We retrieve the chemical composition of the BAL winds launched by the central AGN, synthesize the X-ray features of the hot ISM, and find that (1) the simulated metallicity in the BAL winds could be up to $\sim 8 Z_\odot$, while that of the hot ISM in the host galaxy is $\sim 2.3 Z_\odot$, matching well with SDSS observations of BLR gas; (2) the X-ray emitting hot gas is metal enriched with a typical value $\sim2.5 Z_\odot$; (3) the circumunuclear cold gas disk, where the metals are condensed, further enriched and recycled, plays a critical role in the metal enrichment; (4) the black hole accretion rate $\dot{M}_{\rm BH}$ linearly correlates with the star formation rate $\dot{M}_\star^+$ in the circumnuclear disk, i.e, $\dot{M}_\star^+ \sim 7.7 \dot{M}_{\rm BH}$, but lagged in time by roughly $10^6$ years.

Based on our previous modeling of AGN feedback in isolated elliptical galaxies (Gan et al. 2014) using the MACER (Massive AGN Controlled Ellipticals Resolved) code, we extend and improve the model to include rotation, to limit star formation to regions of high density and low temperature, to facilitate angular momentum transfer via the Toomre instability in gaseous disks, and to improve the treatment of hot mode (low accretion rate) AGN feedback. The model galaxy now has an extended dark matter profile that matches with standard observations, but has a resolution of parsecs in the inner region and resolves the Bondi radius. We find that the results agree reasonably well with a panoply of observations: (1) both AGN activity and star formation are primarily in central cold gaseous disks, are bursty and mainly driven by the Toomre instability; (2)AGN duty cycle agrees well with the Soltan argument, i.e., the AGN spends most of its lifetime when it is in low luminosity (half of time with $L/L_{Edd}<7\times10^{-5}$), while emitting most of its energy when it is in high luminosity (half of radiant energy emitted with $L/L_{Edd}>0.06$); (3) the total star formation is $\sim$ few percents of the initial stellar mass, occurring in the bursts that would be associated with the observed E+A phenomenon. Most of the star formation occurs in the circumnuclear disk of a size <1 kpc, which is in agreement with recent observations; (4) the ISM X-ray luminosity varies within a reasonable range (median $L_{\rm X, ISM}=9.1\times10^{39}$ erg/s) in agreement with observations.

We present chemical abundances and radial velocities of six HII regions in the extremely metal-poor star-forming dwarf galaxy DDO 68. They are derived from deep spectra in the wavelength range 3500 - 10,000 {\AA}, acquired with the Multi Object Double Spectrograph (MODS) at the Large Binocular Telescope (LBT). In the three regions where the [O III]$\lambda$4363 {\AA} line was detected, we inferred the abundance of He, N, O, Ne, Ar, and S through the "direct" method. We also derived the oxygen abundances of all the six regions adopting indirect method calibrations. We confirm that DDO 68 is an extremely metal-poor galaxy, and a strong outlier in the luminosity - metallicity relation defined by star-forming galaxies. With the direct-method we find indeed an oxygen abundance of 12+log(O/H)=7.14$\pm$0.07 in the northernmost region of the galaxy and, although with large uncertainties, an even lower 12+log(O/H)=6.96$\pm$0.09 in the "tail". This is, at face value, the most metal-poor direct abundance detection of any galaxy known. We derive a radial oxygen gradient of -0.06$\pm$0.03 dex/kpc (or -0.30 dex $R_{25}^{-1}$) with the direct method, and a steeper gradient of -0.12$\pm$0.03 dex/kpc (or -0.59 dex $R_{25}^{-1}$) from the indirect method. For the $\alpha$-element to oxygen ratios we obtain values in agreement with those found in other metal-poor star-forming dwarfs. For nitrogen, instead, we infer much higher values, leading to log(N/O)$\sim-1.4$, at variance with the suggested existence of a tight plateau at $-1.6$ in extremely metal poor dwarfs. The derived helium mass fraction ranges from Y=0.240$\pm$0.005 to Y=0.25$\pm$0.02, compatible with standard big bang nucleosynthesis. Finally, we measured HII region radial velocities in the range 479$-$522 km/s from the tail to the head of the "comet", consistent with the rotation derived in the HI.

Understanding how galaxies maintain the inefficiency of star formation with physically self-consistent models is a central problem for galaxy evolution. Although numerous theoretical models have been proposed in recent decades, the debate still exists. By means of high-resolution two-dimensional hydrodynamical simulations, we study the three feedback effects (the stellar wind heating, SNe feedback, and AGN feedback) in suppressing star formation activities on the evolution of early-type galaxies with different stellar masses. AGN feedback models are updated based on \citet{Yuan2018}. The gas sources comes exclusively from the mass losses of dying low-mass stars for most of our models. We find that SNe feedback can keep star formation at a significantly low level for low mass elliptical galaxies for a cosmological evolution time. For the high mass galaxies, AGN feedback can efficiently offset the radiative cooling and thus regulate the star formation activities. Such a suppression of star formation is extremely efficient in the inner region of the galaxies. AGB heating cannot account for this suppression for low and high mass galaxies. The X-ray temperature $T_{\rm X}$ and luminosity $L_{\rm X}$ of hot plasma can be in agreement with the observed data with the inclusion of effective feedback processes. These results thus suggest that we can use $T_{\rm X}$ and $L_{\rm X}$ to probe the role of different feedback processes. The inclusion of additional gas sources can make the mass scale between SNe and AGN feedback dominating in suppressing star formation decrease to an observationally inferred value of a few $10^{10}~M_{\odot}$.

The fully analytical solution for isothermal Bondi accretion on a black hole (MBH) at the center of two-component Jaffe (1983) galaxy models is presented. In a previous work we provided the analytical expressions for the critical accretion parameter and the radial profile of the Mach number in the case of accretion on a MBH at the center of a spherically symmetric one-component Jaffe galaxy model. Here we apply this solution to galaxy models where both the stellar and total mass density distributions are described by the Jaffe profile, with different scale-lengths and masses, and to which a central MBH is added. For such galaxy models all the relevant stellar dynamical properties can also be derived analytically (Ciotti & Ziaee Lorzad 2018). In these new models the hydrodynamical and stellar dynamical properties are linked by imposing that the gas temperature is proportional to the virial temperature of the galaxy stellar component. The formulae that are provided allow to evaluate all flow properties, and are then useful for estimates of the scale-radius and the mass flow rate when modeling accretion on massive black holes at the center of galaxies. As an application, we quantify the departure from the true mass accretion rate of estimates obtained using the gas properties at various distances from the MBH, under the hypothesis of classical Bondi accretion.

The fully analytical solution for isothermal Bondi accretion on a black hole (MBH) at the center of JJ two-component Jaffe (1983) galaxy models is presented. In JJ models the stellar and total mass density distributions are described by the Jaffe profile, with different scale-lengths and masses, and to which a central MBH is added; all the relevant stellar dynamical properties can also be derived analytically. In these new accretion solutions the hydrodynamical and stellar dynamical properties are linked by imposing that the gas temperature is proportional to the virial temperature of the stellar component. The formulae that are provided allow to evaluate all flow properties, and are then useful for estimates of the accretion radius and the mass flow rate when modeling accretion on MBHs at the center of galaxies.

This is the second paper of our series of works of studying the effects of active galactic nuclei (AGN) feedback on the cosmological evolution of an isolated elliptical galaxy by performing two-dimensional high-resolution hydrodynamical numerical simulations. In these simulations, the inner boundary is chosen so that the Bondi radius is resolved. Physical processes like star formation, SNe Ia and II are taken into account. Compared to previous works, the main improvements is that we adopt the most updated AGN physics, which is described in detail in the first paper of this series (Yuan et al. 2018, Paper I). These improvements include the discrimination of the two accretion modes of the central AGN and the most updated descriptions of the wind and radiation in the two modes. In Paper I, we consider the case that the specific angular momentum of the gas in the galaxy is very low. In this paper, we consider the case that the specific angular momentum of the gas is high. In the galactic scale, we adopt the gravitational torques raised due to non-axisymmetric structure in the galaxy as the mechanism of the transfer of angular momentum of gas, as proposed in some recent works. Since our simulations are axisymmetric, we make use of a parameterized prescription to mimic this mechanism. Same as Paper I, we investigate the AGN light curve, typical AGN lifetime, growth of the black hole mass, AGN duty-cycle, star formation, and the X-ray surface brightness of the galaxy. Special attention is paid to the effects of specific angular momentum of the galaxy on these properties. We find that some results are qualitatively similar to those shown in Paper I, while some results such as star formation and black hole growth do show a significant difference due to the mass concentration in the galactic disk as a consequence of galactic rotation.

This paper summarizes our recent works of studying AGN feedback in an isolated elliptical galaxy by performing high-resolution hydrodynamical numerical simulations. Bondi radius is resolved and the mass accretion rate of the black hole is calculated. The most updated AGN physics, namely the discrimination of cold and hot accretion modes and the exact descriptions of the AGN radiation and wind for a given accretion rate are adopted and their interaction with the gas in the host galaxy is calculated. Physical processes such as star formation and SNe feedback are taken into account. Consistent with observation, we find the AGN spends most of the time in the low-luminosity regime. AGN feedback overall suppresses the star formation; but depending on location in the galaxy and time, it can also enhance it. The light curve of specific star formation rate is not synchronous with the AGN light curve. We find that wind usually plays a dominant role in controlling the AGN luminosity and star formation, but radiation also cannot be neglected.

A recent determination of the relationships between the X-ray luminosity of the ISM (Lx) and the stellar and total mass, for a sample of nearby early-type galaxies (ETGs), is used to investigate the origin of the hot gas, via a comparison with the results of hydrodynamical simulations of the ISM evolution for a large set of isolated ETGs. After the epoch of major galaxy formation (after z~2), the ISM is replenished by stellar mass losses and SN ejecta, at the rate predicted by stellar evolution, and is depleted by star formation; it is heated by the thermalization of stellar motions, SNe explosions and the mechanical (from winds) and radiative AGN feedback. The models agree well with the observed relations, even for the largely different Lx values at the same mass, thanks to the sensitivity of the gas flow to many galaxy properties; this holds for models including AGN feedback, and those without. Therefore, the mass input from the stellar population is able to account for a major part of the observed Lx; and AGN feedback, while very important to maintain massive ETGs in a time-averaged quasi-steady state, keeping low star formation and the black hole mass, does not dramatically alter the gas content originating in stellar recycled material. These conclusions are based on theoretical predictions for the stellar population contributions in mass and energy, and on a self-consistent modeling of AGN feedback.

Dynamical properties of spherically symmetric galaxy models where both the stellar and total mass density distributions are described by the Jaffe (1983) profile (with different scale-lenghts and masses), are presented. The orbital structure of the stellar component is described by Osipkov--Merritt anisotropy, and a black hole (BH) is added at the center of the galaxy; the dark matter halo is isotropic. First, the conditions required to have a nowhere negative and monothonically decreasing dark matter halo density profile, are derived. We then show that the phase-space distribution function can be recovered by using the Lambert-Euler $W$ function, while in absence of the central BH only elementary functions appears in the integrand of the inversion formula. The minimum value of the anisotropy radius for consistency is derived in terms of the galaxy parameters. The Jeans equations for the stellar component are solved analytically, and the projected velocity dispersion at the center and at large radii are also obtained analytically for generic values of the anisotropy radius. Finally, the relevant global quantities entering the Virial Theorem are computed analytically, and the fiducial anisotropy limit required to prevent the onset of Radial Orbit Instability is determined as a function of the galaxy parameters. The presented models, even though highly idealized, represent a substantial generalization of the models presentd in Ciotti et al. (2009), and can be useful as starting point for more advanced modeling the dynamics and the mass distribution of elliptical galaxies.

One of the most active fields of research of modern-day astrophysics is that of massive black hole formation and co-evolution with the host galaxy. In these investigations, ranging from cosmological simulations, to semi-analytical modeling, to observational studies, the Bondi solution for accretion on a central point mass is widely adopted. In this work we generalize the classical Bondi accretion theory to take into account the effects of the gravitational potential of the host galaxy, and of radiation pressure in the optically thin limit. Then, we present the fully analytical solution, in terms of the Lambert-Euler $W$-function, for isothermal accretion in Jaffe and Hernquist galaxies with a central black hole. The flow structure is found to be sensitive to the shape of the mass profile of the host galaxy. These results and the formulae that are provided, mostly important the one for the critical accretion parameter, allow for a direct evaluation of all flow properties, and are then useful for the above mentioned studies. As an application, we examine the departure from the true mass accretion rate of estimates obtained using the gas properties at various distances from the black hole, under the hypothesis of classical Bondi accretion. An overestimate is obtained from regions close to the black hole, and an underestimate outside a few Bondi radii; the exact position of the transition between the two kinds of departure depends on the galaxy model.

We continue to investigate the dynamics of collisionless systems of particles interacting via additive $r^{-\alpha}$ interparticle forces. Here we focus on the dependence of the radial-orbit instability on the force exponent $\alpha$. By means of direct $N$-body simulations we study the stability of equilibrium radially anisotropic Osipkov-Merritt spherical models with Hernquist density profile and with $1\leq\alpha<3$. We determine, as a function of $\alpha$, the minimum value for stability of the anisotropy radius $r_{as}$ and of the maximum value of the associated stability indicator $\xi_s$. We find that, for decreasing $\alpha$, $r_{as}$ decreases and $\xi_s$ increases, i.e. longer-range forces are more robust against radial-orbit instability. The isotropic systems are found to be stable for all the explored values of $\alpha$. The end products of unstable systems are all markedly triaxial with minor-to-major axial ratio $>0.3$, so they are never flatter than an E7 system.

We report on five compact, extremely young (<10Myr) and blue (\beta_UV<-2.5, F_\lambda =\lambda^\beta) objects observed with VLT/MUSE at redshift 3.1169, 3.235, in addition to three objects at z=6.145. These sources are magnified by the Hubble Frontier Field galaxy clusters MACS~J0416 and AS1063. Their de-lensed half light radii (Re) are between 16 to 140pc, the stellar masses are ~1-20 X 10^6 Msun, the magnitudes are m_uv=28.8 - 31.4 (-17<Muv<-15) and specific star formation rates can be as large as ~800Gyr^-1. Multiple images of these systems are widely separated in the sky (up to 50'') and individually magnified by factors 3-40. Remarkably, the inferred physical properties of two objects are similar to those expected in some globular cluster formation scenarios, representing the best candidate proto-globular clusters (proto-GC) discovered so far. Rest-frame optical high dispersion spectroscopy of one of them at z=3.1169 yields a velocity dispersion \sigma_v~20km/s, implying a dynamical mass dominated by the stellar mass. Another object at z=6.145, with de-lensed Muv ~ -15.3 (m_uv ~ 31.4), shows a stellar mass and a star-formation rate surface density consistent with the values expected from popular GC formation scenarios. An additional star-forming region at z=6.145, with de-lensed m_uv ~ 32, a stellar mass of 0.5 X 10^6 Msun and a star formation rate of 0.06 Msun/yr is also identified. These objects currently represent the faintest spectroscopically confirmed star-forming systems at z>3, elusive even in the deepest blank fields. We discuss how proto-GCs might contribute to the ionization budget of the universe and augment Lya visibility during reionization. This work underlines the crucial role of JWST in characterizing the rest-frame optical and near-infrared properties of such low-luminosity high-z objects.

We present 2D hydrodynamical simulations for the evolution of early-type galaxies containing central massive black holes (MBHs), starting at age 2 Gyr. The code contains accurate and physically consistent radiative and mechanical AGN wind feedback, with parsec-scale central resolution. Mass input comes from stellar evolution; energy input includes Type Ia and II supernova and stellar heating; star-formation is included. Realistic, axisymmetric dynamical models for the galaxies are built solving the Jeans' equations. The lowest mass models (Mstar = 8 10^{10}Msun) develop global outflows sustained by SNIa's heating, ending with a significantly lower amount of hot gas and new stars. In more massive models, nuclear outbursts last to the present epoch, with large and frequent fluctuations in nuclear emission and from the gas (Lx). Each burst last ~ 10^{7.5} yr, during which (for r < 2-3 kpc) cold, inflowing, and hot, outflowing gas phases coexist. The Lx-T relation for the gas matches that of local galaxies. AGN activity causes positive feedback for star formation. Roughly half of the total mass loss is recycled into new stars (DeltaMstar), just ~ 3% of it is accreted on the MBH, the remainder being ejected from the galaxy. The ratio between the mass of gas expelled to that in to new stars, the load factor, is ~0.6. Rounder galaxies shapes lead to larger final MBH masses, DeltaMstar, and Lx. Almost all the time is spent at very low nuclear luminosities, yet one quarter of the total energy is emitted at an Eddington ratio > 0.1. The duty-cycle of AGN activity approximates 4% (Abridged).

We present new photometry of the dwarf irregular galaxy DDO 68, one of the most metal-poor and least massive dwarfs, located in the Lynx-Cancer Void. The images were acquired with the Large Binocular Telescope in the g and r passbands, and show unequivocally that DDO 68 has previously unknown stellar streams related to the accretion of at least two smaller companions: a flea with smaller fleas biting it, to put it in Jonathan Swift's words. Our data provide direct observational evidence of multiple merging occurring at very low galactic mass scales. We present the results of an N-body simulation of the interaction of three dwarf galaxies which reproduces well the main morphological features of DDO 68.

The dark matter (DM) haloes around spiral galaxies appear to conspire with their baryonic content: empirically, significant amounts of DM are inferred only below a universal characteristic acceleration scale. Moreover, the discrepancy between the baryonic and dynamical mass, which is usually interpreted as the presence of DM, follows a very tight mass discrepancy acceleration (MDA) relation. Its universality, and its tightness in spiral galaxies, poses a challenge for the DM interpretation and was used to argue in favour of MOdified Newtonian Dynamics (MOND). Here, we test whether or not this applies to early-type galaxies. We use the dynamical models of fast-rotator early-type galaxies by Cappellari et al. based on ATLAS$^{3D}$ and SLUGGS data, which was the first homogenous study of this kind, reaching ~4 $R_e$, where DM begins to dominate the total mass budget. We find the early-type galaxies to follow an MDA relation similar to spiral galaxies, but systematically offset. Also, while the slopes of the mass density profiles inferred from galaxy dynamics show consistency with those expected from their stellar content assuming MOND, some profiles of individual galaxies show discrepancies.

Accretion onto central massive black holes in galaxies is often modelled with the Bondi solution. In this paper we study a generalization of the classical Bondi accretion theory, considering the additional effects of the gravitational potential of the host galaxy, and of electron scattering in the optically thin limit. We provide a general analysis of the bias in the estimates of the Bondi radius and mass accretion rate, when adopting as fiducial values for the density and temperature at infinity the values of these quantities measured at finite distance from the central black hole. We also give general formulae to compute the correction terms of the critical accretion parameter in relevant asymptotic regimes. A full analytical discussion is presented in the case of an Hernquist galaxy, when the problem reduces to the discussion of a cubic equation, therefore allowing for more than one critical point in the accretion structure. The results are useful for observational works (especially in the case of low-luminosity systems), as well as for numerical simulations, where accretion rates are usually defined in terms of the gas properties near the

AGN feedback from supermassive black holes (SMBHs) at the center of early type galaxies is commonly invoked as the explanation for the quenching of star formation in these systems. The situation is complicated by the significant amount of mass injected in the galaxy by the evolving stellar population over cosmological times. In absence of feedback, this mass would lead to unobserved galactic cooling flows, and to SMBHs two orders of magnitude more massive than observed. By using high-resolution 2D hydrodynamical simulations with radiative transport and star formation in state-of-the-art galaxy models, we show how the intermittent AGN feedback is highly structured on spatial and temporal scales, and how its effects are not only negative (shutting down the recurrent cooling episodes of the ISM), but also positive, inducing star formation in the inner regions of the host galaxy.

Early-type galaxies (ETGs) host a hot ISM produced mainly by stellar winds, and heated by Type Ia supernovae and the thermalization of stellar motions. High resolution 2D hydrodynamical simulations showed that ordered rotation in the stellar component results in the formation of a centrifugally supported cold equatorial disc. In a recent numerical investigation we found that subsequent generations of stars are formed in this cold disc; this process consumes most of the cold gas, leaving at the present epoch cold masses comparable to those observed. Most of the new stellar mass formed a few Gyrs ago, and resides in a disc.

We use the Reversibility Error Method and the Fidelity to analyze the global effects of a small perturbation in a non-integrable system. Both methods have already been proposed and used in the literature but the aim of this paper is to compare them in a physically significant example adding some considerations on the equivalence, observed in this case, between round-off and random perturbations. As a paradigmatic example we adopt the restricted planar circular three body problem. The cumulative effect of random perturbations or round-off leads to a divergence of the perturbed orbit from the reference one. Rather than computing the distance of the perturbed orbit from the reference one, after a given number n of iterations, a procedure we name the Forward Error Method (FEM), we measure the distance of the reversed orbit (n periods forward and backward) from the initial point. This approach, that we name Reversibility Error Method (REM), does not require the computation of the unperturbed map. The loss of memory of the perturbed map is quantified by the Fidelity decay rate whose computation requires a statistical average over an invariant region.

We continue the study of collisionless systems governed by additive $r^{-\alpha}$ interparticle forces by focusing on the influence of the force exponent $\alpha$ on radial orbital anisotropy. In this preparatory work we construct the radially anisotropic Osipkov-Merritt phase-space distribution functions for self-consistent spherical Hernquist models with $r^{-\alpha}$ forces and $1\leq\alpha<3$. The resulting systems are isotropic at the center and increasingly dominated by radial orbits at radii larger than the anisotropy radius $r_a$. For radially anisotropic models we determine the minimum value of the anisotropy radius $r_{ac}$ as a function of $\alpha$ for phase-space consistency (such that the phase-space distribution function is nowhere negative for $r_a\geq r_{ac}$). We find that $r_{ac}$ decreases for decreasing $\alpha$, and that the amount of kinetic energy that can be stored in the radial direction relative to that stored in the tangential directions for marginally consistent models increases for decreasing $\alpha$. In particular, we find that isotropic systems are consistent in the explored range of $\alpha$. By means of direct $N$-body simulations we finally verify that the isotropic systems are also stable.

High resolution 2D hydrodynamical simulations describing the evolution of the hot ISM in axisymmetric two-component models of early-type galaxies well reproduced the observed trends of the X-ray luminosity ($L_\mathrm{x}$) and temperature ($T_\mathrm{x}$) with galaxy shape and rotation, however they also revealed the formation of an exceedingly massive cooled gas disc in rotating systems. In a follow-up of this study, here we investigate the effects of star formation in the disc, including the consequent injection of mass, momentum and energy in the pre-existing interstellar medium. It is found that subsequent generations of stars originate one after the other in the equatorial region; the mean age of the new stars is $> 5$ Gyr, and the adopted recipe for star formation can reproduce the empirical Kennicutt-Schmidt relation. The results of the previous investigation without star formation, concerning $L_\mathrm{x}$ and $T_\mathrm{x}$ of the hot gas, and their trends with galactic shape and rotation, are confirmed. At the same time, the consumption of most of the cold gas disc into new stars leads to more realistic final systems, whose cold gas mass and star formation rate agree well with those observed in the local universe. In particular, our models could explain the observation of kinematically aligned gas in massive, fast-rotating early-type galaxies.

In the last years, stellar migration in galactic discs has been the subject of several investigations. However, its impact on the chemical evolution of the Milky Way still needs to be fully quantified. In this paper, we aim at imposing some constraints on the significance of this phenomenon by considering its influence on the chemical evolution of the Milky Way thin disc. We do not investigate the physical mechanisms underlying the migration of stars. Rather, we introduce a simple, heuristic treatment of stellar migration in a detailed chemical evolution model for the thin disc of the Milky Way, which already includes radial gas flows and reproduces several observational constraints for the solar vicinity and the whole Galactic disc. When stellar migration is implemented according to the results of chemo-dynamical simulations by Minchev et. al. (2013) and finite stellar velocities of 1 km s$^{-1}$ are taken into account, the high-metallicity tail of the metallicity distribution function of long-lived thin-disc stars is well reproduced. By exploring the velocity space, we find that the migrating stars must travel with velocities in the range 0.5 -2 km s$^{-1}$ to properly reproduce the high-metallicity tail of the metallicity distribution. We confirm previous findings by other authors that the observed spread in the age-metallicity relation of solar neighbourhood stars can be explained by the presence of stars which originated at different Galactocentric distances, and we conclude that the chemical properties of stars currently observed in the solar vicinity do suggest that stellar migration is present to some extent.

We present the analytical solution of the two-integrals Jeans equations for Miyamoto-Nagai discs embedded in Binney logarithmic dark matter haloes. The equations can be solved (both with standard methods and with the Residue Theorem) for arbitrary choices of the parameters, thus providing a very flexible two-component galaxy model, ranging from flattened discs to spherical systems. A particularly interesting case is obtained when the dark matter halo reduces to the Singular Isothermal Sphere. Azimuthal motions are separated in the ordered and velocity dispersion components by using the Satoh decomposition. The obtained formulae can be used in numerical simulations of galactic gas flows, for testing codes of stellar dynamics, and to study the dependence of the stellar velocity dispersion and of the asymmetric drift in the equatorial plane as a function of disc and halo flattenings. Here, we estimate the inflow radial velocities of the interstellar medium, expected by the mixing of the stellar mass losses of the lagging stars in the disc with a pre-existing gas in circular orbit.

We present an approach to the design of distribution functions that depend on the phase-space coordinates through the action integrals. The approach makes it easy to construct a dynamical model of a given stellar component. We illustrate the approach by deriving distribution functions that self-consistently generate several popular stellar systems, including the Hernquist, Jaffe, and Navarro, Frenk and White models. We focus on non-rotating spherical systems, but extension to flattened and rotating systems is trivial. Our distribution functions are easily added to each other and to previously published distribution functions for discs to create self-consistent multi-component galaxies. The models this approach makes possible should prove valuable both for the interpretation of observational data and for exploring the non-equilibrium dynamics of galaxies via N-body simulation.

We present an investigation about the shape of the initial mass function (IMF) of early-type galaxies (ETGs), based on a joint lensing and dynamical analysis, and on stellar population synthesis models, for a sample of 55 lens ETGs identified by the Sloan Lens ACS (SLACS) Survey. We construct axisymmetric dynamical models based on the Jeans equations which allow for orbital anisotropy and include a dark matter halo. The models reproduce in detail the observed \textit{HST} photometry and are constrained by the total projected mass within the Einstein radius and the stellar velocity dispersion ($\sigma$) within the SDSS fibers. Comparing the dynamically-derived stellar mass-to-light ratios $(M_*/L)_{\rm dyn}$, obtained for an assumed halo slope $\rho_{\rm h}\propto r^{-1}$, to the stellar population ones $(M_*/L)_{\rm pop}$, derived from full-spectrum fitting and assuming a Salpeter IMF, we infer the mass normalization of the IMF. Our results confirm the previous analysis by the SLACS team that the mass normalization of the IMF of high $\sigma$ galaxies is consistent on average with a Salpeter slope. Our study allows for a fully consistent study of the trend between IMF and $\sigma$ for both the SLACS and \ATLAS samples, which explore quite different $\sigma$ ranges. The two samples are highly complementary, the first being essentially $\sigma$ selected, and the latter volume-limited and nearly mass selected. We find that the two samples merge smoothly into a single trend of the form $\log\alpha =(0.38\pm0.04)\times\log(\sigma_{\rm e}/200\,\mathrm{km~s}^{-1})+(-0.06\pm0.01)$, where $\alpha=(M_*/L)_{\rm dyn}/(M_*/L)_{\rm pop}$ and $\sigma_{\rm e}$ is the luminosity averaged $\sigma$ within one effective radius $R_{\rm e}$. This is consistent with a systematic variation of the IMF normalization from Kroupa to Salpeter in the interval $\sigma_{\rm e}\approx90-270\,\mathrm{km~s}^{-1}$.

We use the Chandrasekhar formalism and direct N-body simulations to study the effect of dynamical friction on a test object only slightly more massive than the field stars, orbiting a spherically symmetric background of particles with a mass spectrum. The main goal is to verify whether the dynamical friction time (t_DF) develops a non-monotonic radial-dependence that could explain the bimodality of the Blue Straggler radial distributions observed in globular clusters. In these systems, in fact, relaxation effects lead to a mass and velocity radial segregation of the different mass components, so that mass-spectrum effects on t_DF are expected to be dependent on radius. We find that, in spite of the presence of different masses, t_DF is always a monotonic function of radius, at all evolutionary times and independently of the initial concentration of the simulated cluster. This because the radial dependence of t_DF is largely dominated by the total mass density profile of the background stars (which is monotonically decreasing with radius). Hence, a progressive temporal erosion of the BSS population at larger and larger distances from the cluster center remains the simplest and the most likely explanation of the shape of the observed BSS radial distributions, as suggested in previous works. We also confirm the theoretical expectation that approximating a multi-mass globular cluster as made of (averaged) equal-mass stars can lead to significant overestimates of t_DF within the half-mass radius.

By means of high resolution 2D hydrodynamical simulations, we study the evolution of the hot ISM for a large set of early-type galaxy models, characterized by various degrees of flattening and internal rotation. The galaxies are described by state-of-the-art axisymmetric two-component models, tailored to reproduce real systems; the dark matter haloes follow the Navarro-Frenk-White or the Einasto profile. The gas is produced by the evolving stars, and heated by Type Ia SNe. We find that, in general, the rotation field of the ISM in rotating galaxies is very similar to that of the stars, with a consequent negligible heating contribution from thermalization of the ordered motions. The relative importance of flattening and rotation in determining the final X-ray luminosity $L_x$ and temperature $T_x$ of the hot haloes is a function of the galactic mass. Flattening and rotation in low mass galaxies favour the establishment of global winds, with the consequent reduction of $L_x$. In medium-to-high mass galaxies, flattening and rotation are not sufficient to induce global winds, however, in the rotating models the nature of the gas flows is deeply affected by conservation of angular momentum, resulting in a reduction of both $L_x$ and $T_x$.

Based on two-dimensional high resolution hydrodynamic numerical simulation, we study the mechanical and radiative feedback effects from the central AGN on the cosmological evolution of an isolated elliptical galaxy. Physical processes such as star formation and supernovae are considered. The inner boundary of the simulation domain is carefully chosen so that the fiducial Bondi radius is resolved and the accretion rate of the black hole is determined self-consistently. In analogy to previous works, we assume that the specific angular momentum of the galaxy is low. It is well-known that when the accretion rates are high and low, the central AGNs will be in cold and hot accretion modes, which correspond to the radiative and kinetic feedback modes, respectively. The emitted spectrum from the hot accretion flows is harder than that from the cold accretion flows, which could result in a higher Compton temperature accompanied by a more efficient radiative heating, according to previous theoretical works. Such a difference of the Compton temperature between the two feedback modes, the focus of this study, has been neglected in previous works. Significant differences in the kinetic feedback mode are found as a result of the stronger Compton heating and accretion becomes more chaotic. More importantly, if we constrain models to correctly predict black hole growth and AGN duty cycle after cosmological evolution, we find that the favored model parameters are constrained: mechanical feedback efficiency diminishes with decreasing luminosity (the maximum efficiency being $\simeq 10^{-3.5}$) and X-ray Compton temperature increases with decreasing luminosity, although models with fixed mechanical efficiency and Compton temperature can be found that are satisfactory as well. We conclude that radiative feedback in the kinetic mode is much more important than previously thought.

We study the effects of a detailed dust treatment on the properties and evolution of early-type galaxies containing central black holes, as determined by AGN feedback. We find that during cooling flow episodes, radiation pressure on the dust in and interior to infalling shells of cold gas can greatly impact the amount of gas able to be accreted and therefore the frequency of AGN bursts. However, the overall hydrodynamic evolution of all models, including mass budget, is relatively robust to the assumptions on dust. We find that IR re-emission from hot dust can dominate the bolometric luminosity of the galaxy during the early stages of an AGN burst, reaching values in excess of $10^{46}$ erg/s. The AGN-emitted UV is largely absorbed, but the optical depth in the IR does not exceed unity, so the radiation momentum input never exceeds $L_{\rm BH}/c$. We constrain the viability of our models by comparing the AGN duty cycle, broadband luminosities, dust mass, black hole mass, and other model predictions to current observations. These constraints force us to models wherein the dust to metals ratios are $\simeq 1%$ of the Galactic value, and only models with a dynamic dust to gas ratio are able to produce both quiescent galaxies consistent with observations and high obscured fractions during AGN "on" phases. During AGN outbursts, we predict that a large fraction of the FIR luminosity can be attributed to warm dust emission ($\simeq100$ K) from dense dusty gas within $\leq 1$ kpc reradiating the AGN UV emission.

Early-type galaxies (ETGs) are observed to be more compact, on average, at $z \gtrsim 2$ than at $z\simeq 0$, at fixed stellar mass. Recent observational works suggest that such size evolution could reflect the similar evolution of the host dark matter halo density as a function of the time of galaxy quenching. We explore this hypothesis by studying the distribution of halo central velocity dispersion ($\sigma_0$) and half-mass radius ($r_{\rm h}$) as functions of halo mass $M$ and redshift $z$, in a cosmological $\Lambda$-CDM $N$-body simulation. In the range $0\lesssim z\lesssim 2.5$, we find $\sigma_0\propto M^{0.31-0.37}$ and $r_{\rm h}\propto M^{0.28-0.32}$, close to the values expected for homologous virialized systems. At fixed $M$ in the range $10^{11} M_\odot \lesssim M\lesssim 5.5 \times 10^{14} M_\odot$ we find $\sigma_0\propto(1+z)^{0.35}$ and $r_{\rm h}\propto(1+z)^{-0.7}$. We show that such evolution of the halo scaling laws is driven by individual haloes growing in mass following the evolutionary tracks $\sigma_0\propto M^{0.2}$ and $r_{\rm h}\propto M^{0.6}$, consistent with simple dissipationless merging models in which the encounter orbital energy is accounted for. We compare the $N$-body data with ETGs observed at $0\lesssim z\lesssim3$ by populating the haloes with a stellar component under simple but justified assumptions: the resulting galaxies evolve consistently with the observed ETGs up to $z \simeq 2$, but the model has difficulty reproducing the fast evolution observed at $z\gtrsim 2$. We conclude that a substantial fraction of the size evolution of ETGs can be ascribed to a systematic dependence on redshift of the dark matter haloes structural properties.

Observational and numerical studies gave hints that the hot gaseous haloes of ETGs may be sensitive to the galaxy internal kinematics. By using high resolution 2D hydro simulations, and realistic two-component (stars plus dark matter) axisymmetric galaxy models, we study the evolution of the hot haloes in a suite of flat ETGs of fixed mass distribution, but with variable amounts of azimuthal velocity dispersion and rotational support, including the possibility of a counter-rotating inner stellar disc. The hot halo is fed by stellar mass losses and heated by SNIa explosions and thermalization of stellar motions. We measure the value of the thermalization parameter gamma (the ratio between the heating due to the relative velocity between the stellar streaming and the ISM bulk flow, and the heating attainable by complete thermalization of the stellar streaming motions). We find that 1) the X-ray emission and the average temperature are larger in fully velocity dispersion supported systems; 2) 0.1<gamma<0.2 for isotropic rotators (with a trend for being larger for lower dark mass models); 3) systems that are isotropic rotators at large radii with an inner counter-rotating disc, or fully velocity dispersion supported systems with an inner rotating disc, have gamma=1, again with a trend to increase for lower dark mass contents. We also find that the lower X-ray luminosities of isotropic rotators cannot be explained just by their low gamma, but are due to the complicated flow structure, consequence of the angular momentum stored at large radii. X-ray emission weighted temperatures and luminosities nicely match observed values; the X-ray isophotes are boxy in case of significant galaxy rotation. Overall, it is found that rotation has an important role to explain the observational result that more rotationally supported ETGs on average show a lower X-ray emission [abridged].

The mass returned to the ambient medium by aging stellar populations over cosmological times sums up to a significant fraction (20% - 30% or more) of their initial mass. This continuous mass injection plays a fundamental role in phenomena such as galaxy formation and evolution, fueling of supermassive black holes in galaxies and the consequent (negative and positive) feedback phenomena, and the origin of multiple stellar populations in globular clusters. In numerical simulations the calculation of the mass return can be time consuming, since it requires at each time step the evaluation of a convolution integral over the whole star formation history, so the computational time increases quadratically with the number of time-steps. The situation can be especially critical in hydrodynamical simulations, where different grid points are characterized by different star formation histories, and the gas cooling and heating times are shorter by orders of magnitude than the characteristic stellar lifetimes. In this paper we present a fast and accurate method to compute the mass return from stellar populations undergoing arbitrarily complicated star formation histories. At each time-step the mass return is calculated from its value at the previous time, and the star formation rate over the last time-step only. Therefore in the new scheme there is no need to store the whole star formation history, and the computational time increases linearly with the number of time-steps.

We present a detailed diagnostic study of the observed temperatures of the hot X-ray coronae of early-type galaxies. By extending the investigation carried out in Pellegrini (2011) with spherical models, we focus on the dependence of the energy budget and temperature of the hot gas on the galaxy structure and internal stellar kinematics. By solving the Jeans equations we construct realistic axisymmetric three-component galaxy models (stars, dark matter halo, central black hole) with different degrees of flattening and rotational support. The kinematical fields are projected along different lines of sight, and the aperture velocity dispersion is computed within a fraction of the circularized effective radius. The model parameters are chosen so that the models resemble real ETGs and lie on the Faber-Jackson and Size-Luminosity relations. For these models we compute T_* (the stellar heating contribution to the gas injection temperature) and T_gm (the temperature equivalent of the energy required for the gas escape). In particular, different degrees of thermalisation of the ordered rotational field of the galaxy are considered. We find that T_* and T_gm can vary only mildly due to a pure change of shape. Galaxy rotation instead, when not thermalised, can lead to a large decrease of T_*; this effect can be larger in flatter galaxies that can be more rotationally supported. Recent temperature measurements T_x, obtained with Chandra, are larger than, but close to, the T_* values of the models, and show a possible trend for a lower T_x in flatter and more rotationally supported galaxies; this trend can be explained by the lack of thermalisation of the whole stellar kinetic energy. Flat and rotating galaxies also show lower L_x values, and then a lower gas content, but this is unlikely to be due to the small variation of T_gm found here for them.

We explore how the co-evolution of massive black holes (MBHs) and galaxies is affected by environmental effects, addressing in particular MBHs hosted in the central galaxies of clusters (we will refer to these galaxies in general as 'CGs'). Recently the sample of MBHs in CGs with dynamically measured masses has increased, and it has been suggested that these MBH masses (M_BH) deviate from the expected correlations with velocity dispersion (sigma) and mass of the bulge (M_bulge) of the host galaxy: MBHs in CGs appear to be `over-massive'. This discrepancy is more pronounced when considering the M_BH-sigma relation than the M_BH-M_bulge one. We show that this behavior stems from a combination of two natural factors, (i) that CGs experience more mergers involving spheroidal galaxies and their MBHs, and (ii) that such mergers are preferentially gas-poor. We use a combination of analytical and semi-analytical models to investigate the MBH-galaxy co-evolution in different environments and find that the combination of these two factors explains the trends observed in current data-sets.

In Newtonian gravity the final states of cold dissipationless collapses are characterized by several structural and dynamical properties remarkably similar to those of observed elliptical galaxies. Are these properties a peculiarity of the Newtonian force or a more general feature of long-range forces? We study this problem by means of $N-$body simulations of dissipationless collapse of systems of particles interacting via additive $r^{-\alpha}$ forces. We find that most of the results holding in Newtonian gravity are also valid for $\alpha\neq2$. In particular the end products are triaxial and never flatter than an E7 system, their surface density profiles are well described by the S\'ersic law, the global density slope-anisotropy inequality is obeyed, the differential energy distribution is an exponential over a large range of energies (for $\alpha\geq1$), and the pseudo phase-space density is a power law of radius. In addition, we show that the process of virialization takes longer (in units of the system's dynamical time) for decreasing values of $\alpha$, and becomes infinite for $\alpha=-1$ (the harmonic oscillator). This is in agreement with the results of deep-MOND collapses (qualitatively corresponding to $\alpha=1$) and it is due to the fact the force becomes more and more similar to the $\alpha=-1$ case, where as well known no relaxation can happen and the system oscillates forever.

Elliptical galaxies have hot coronae with X-ray luminosities and mean gas temperatures that span over wide ranges. This variation can be partially due to the energy budget of the hot gas, that depends on the host galaxy structure and internal kinematics. With the aid of realistic axisymmetric galaxy models, we performed a diagnostic study focussed on the effects of galaxy flattening and rotational support on the hot gas temperature.

With 2D hydrodynamical simulations, we study the evolution of the hot gas flows in early-type galaxies, focussing on the effects of galaxy rotation on the thermal and dynamical status of the ISM. The galaxy is modelled as a two-component axisymmetric system (stars and dark matter), with a variable amount of azimuthal velocity dispersion and rotational support; the presence of a counter rotating stellar disk is also considered. It is found that the ISM of the rotationally supported (isotropic) model is more prone to thermal instabilities than the fully velocity dispersion counterpart, while its ISM temperature and X-ray luminosity are lower. The model with counter rotation shows an intermediate behaviour.

To facilitate the study of black hole fueling, star formation, and feedback in galaxies, we outline a method for treating the radial forces on interstellar gas due to absorption of photons by dust grains. The method gives the correct behavior in all of the relevant limits (dominated by the central point source; dominated by the distributed isotropic source; optically thin; optically thick to UV/optical; optically thick to IR) and reasonably interpolates between the limits when necessary. The method is explicitly energy conserving so that UV/optical photons that are absorbed are not lost, but are rather redistributed to the IR where they may scatter out of the galaxy. We implement the radiative transfer algorithm in a two-dimensional hydrodynamical code designed to study feedback processes in the context of early-type galaxies. We find that the dynamics and final state of simulations are measurably but only moderately affected by radiative forces on dust, even when assumptions about the dust-to-gas ratio are varied from zero to a value appropriate for the Milky Way. In simulations with high gas densities designed to mimic ULIRGs with a star formation rate of several hundred solar masses per year, dust makes a more substantial contribution to the dynamics and outcome of the simulation. We find that, despite the large opacity of dust to UV radiation, the momentum input to the flow from radiation very rarely exceeds L/c due to two factors: the low opacity of dust to the re-radiated IR and the tendency for dust to be destroyed by sputtering in hot gas environments. We also develop a simplification of our radiative transfer algorithm that respects the essential physics but is much easier to implement and requires a fraction of the computational cost.

We study mass models that correspond to MOND (triaxial) potentials for which the Hamilton-Jacobi equation separates in ellipsoidal coordinates. The problem is first discussed in the simpler case of deep-MOND systems, and then generalized to the full MOND regime. We prove that the Kuzmin property for Newtonian gravity still holds, i.e., that the density distribution of separable potentials is fully determined once the density profile along the minor axis is assigned. At variance with the Newtonian case, the fact that a positive density along the minor axis leads to a positive density everywhere remains unproven. We also prove that (i) all regular separable models in MOND have a vanishing density at the origin, so that they would correspond to centrally dark-matter dominated systems in Newtonian gravity; (ii) triaxial separable potentials regular at large radii and associated with finite total mass leads to density distributions that at large radii are not spherical and decline as ln(r)/r^5; (iii) when the triaxial potentials admit a genuine Frobenius expansion with exponent 0<epsilon<1, the density distributions become spherical at large radii, with the profile ln(r)/r^(3+2epsilon). After presenting a suite of positive density distributions associated with MOND separable potentials, we also consider the important family of (non-separable) triaxial potentials V_1 introduced by de Zeeuw and Pfenniger, and we show that, as already known for Newtonian gravity, they obey the Kuzmin property also in MOND. The ordinary differential equation relating their potential and density along the z-axis is an Abel equation of the second kind that, in the oblate case, can be explicitly reduced to canonical form.

We extend the black hole (BH) feedback models of Ciotti, Ostriker, and Proga to two dimensions. In this paper, we focus on identifying the differences between the one-dimensional and two-dimensional hydrodynamical simulations. We examine a normal, isolated $L_*$ galaxy subject to the cooling flow instability of gas in the inner regions. Allowance is made for subsequent star formation, Type Ia and Type II supernovae, radiation pressure, and inflow to the central BH from mildly rotating galactic gas which is being replenished as a normal consequence of stellar evolution. The central BH accretes some of the infalling gas and expels a conical wind with mass, momentum, and energy flux derived from both observational and theoretical studies. The galaxy is assumed to have low specific angular momentum in analogy with the existing one-dimensional case in order to isolate the effect of dimensionality. The code then tracks the interaction of the outflowing radiation and winds with the galactic gas and their effects on regulating the accretion. After matching physical modeling to the extent possible between the one-dimensional and two-dimensional treatments, we find essentially similar results in terms of BH growth and duty cycle (fraction of the time above a given fraction of the Eddington luminosity). In the two-dimensional calculations, the cool shells forming at 0.1--1 kpc from the center are Rayleigh--Taylor unstable to fragmentation, leading to a somewhat higher accretion rate, less effective feedback, and a more irregular pattern of bursting compared to the one-dimensional case.

We investigate how environmental effects by gas stripping alter the growth of a super massive black hole (SMBH) and its host galaxy evolution, by means of 1D hydrodynamical simulations that include both mechanical and radiative AGN feedback effects. By changing the truncation radius of the gas distribution (R_t), beyond which gas stripping is assumed to be effective, we simulate possible environments for satellite and central galaxies in galaxy clusters and groups. The continuous escape of gas outside the truncation radius strongly suppresses star formation, while the growth of the SMBH is less affected by gas stripping because the SMBH accretion is primarily ruled by the density of the central region. As we allow for increasing environmental effects - the truncation radius decreasing from about 410 to 50 kpc - we find that the final SMBH mass declines from about 10^9 to 8 x 10^8 Msol, but the outflowing mass is roughly constant at about 2 x 10^10 Msol. There are larger change in the mass of stars formed, which declines from about 2 x 10^10 to 2 x 10^9 Msol, and the final thermal X-ray gas, which declines from about 10^9 to 5 x 10^8 Msol, with increasing environmental stripping. Most dramatic is the decline in the total time that the objects would be seen as quasars, which declines from 52 Myr (for R_t = 377 kpc) to 7.9 Myr (for R_t = 51 kpc). The typical case might be interpreted as a red and dead galaxy having episodic cooling flows followed by AGN feedback effects resulting in temporary transitions of the overall galaxy color from red to green or to blue, with (cluster) central galaxies spending a much larger fraction of their time in the elevated state than do satellite galaxies.(Abridged)

Accretion of fresh gas at a rate of ~ 1 M_{sun} yr^{-1} is necessary in star-forming disc galaxies, such as the Milky Way, in order to sustain their star-formation rates. In this work we present the results of a new hydrodynamic simulation supporting the scenario in which the gas required for star formation is drawn from the hot corona that surrounds the star-forming disc. In particular, the cooling of this hot gas and its accretion on to the disc are caused by the passage of cold galactic fountain clouds through the corona.

The ISM evolution of elliptical galaxies experiencing feedback from accretion onto a central black hole was studied recently with high-resolution 1D hydrodynamical simulations including radiative heating and pressure effects, a RIAF-like radiative efficiency, mechanical input from AGN winds, and accretion-driven starbursts. Here we focus on the observational properties of the models in the X-ray band (nuclear luminosity; hot ISM luminosity and temperature; temperature and brightness profiles during quiescence and during outbursts). The nuclear bursts last for ~10^7 yr, with a duty-cycle of a few X (10^-3-10^-2); the present epoch bolometric nuclear emission is very sub-Eddington. The ISM thermal luminosity \lx oscillates in phase with the nuclear one; this helps reproduce statistically the observed large \lx variation. In quiescence the temperature profile has a negative gradient; thanks to past outbursts, the brightness profile lacks the steep shape typical of inflowing models. Outbursts produce disturbances in these profiles. Most significantly, a hot bubble from shocked hot gas is inflated at the galaxy center; the bubble would be conical in shape, and show radio emission. The ISM resumes a smooth appearance on a time-scale of ~200 Myr; the duty-cycle of perturbances in the ISM is of the order of 5-10%. From the present analysis, additional input physics is important in the ISM-black hole coevolution, to fully account for the properties of real galaxies, as a confining external medium and a jet. The jet will reduce further the mass available for accretion (and then the Eddington ratio $l$), and may help, together with an external pressure, to produce flat or positive temperature gradient profiles (observed in high density environments). Alternatively, $l$ can be reduced if the switch from high to low radiative efficiency takes place at a larger $l$ than routinely assumed.

In galaxies like the Milky Way, cold (~ 10^4 K) gas ejected from the disc by stellar activity (the so-called galactic-fountain gas) is expected to interact with the virial-temperature (~ 10^6 K) gas of the corona. The associated transfer of momentum between cold and hot gas has important consequences for the dynamics of both gas phases. We quantify the effects of such an interaction using hydrodynamical simulations of cold clouds travelling through a hot medium at different relative velocities. Our main finding is that there is a velocity threshold between clouds and corona, of about 75 km/s, below which the hot gas ceases to absorb momentum from the cold clouds. It follows that in a disc galaxy like the Milky Way a static corona would be rapidly accelerated: the corona is expected to rotate and to lag, in the inner regions, by ~ 80-120 km/s with respect to the cold disc. We also show how the existence of this velocity threshold can explain the observed kinematics of the cold extra-planar gas.

The importance of feedback (radiative and mechanical) from massive black holes at the centers of elliptical galaxies is not in doubt, given the well established relation among black hole mass and galaxy optical luminosity. Here, with the aid of high-resolution hydrodynamical simulations, we discuss how this feedback affects the hot ISM of isolated elliptical galaxies of different mass. The cooling and heating functions include photoionization plus Compton heating, the radiative transport equations are solved, and the mechanical feedback due to the nuclear wind is also described on a physical basis; star formation is considered. In the medium-high mass galaxies the resulting evolution is highly unsteady. At early times major accretion episodes caused by cooling flows in the recycled gas produced by stellar evolution trigger AGN flaring: relaxation instabilities occur so that duty cycles are small enough to account for the very small fraction of massive ellipticals observed to be in the QSO-phase, when the accretion luminosity approaches the Eddington luminosity. At low redshift all models are characterized by smooth, very sub-Eddington mass accretion rates. The mass accumulated by the central black hole is limited to range observed today, even though the mass lost by the evolving stellar population is roughly two order of magnitude larger than the black hole masses observed in elliptical galaxies.

The process of relaxation of a system of particles interacting with long-range forces is relevant to many areas of Physics. For obvious reasons, in Stellar Dynamics much attention has been paid to the case of 1/r^2 force law. However, recently the interest in alternative gravities emerged, and significant differences with respect to Newtonian gravity have been found in relaxation phenomena. Here we begin to explore this matter further, by using a numerical model of spherical shells interacting with an 1/r^alpha force law obeying the superposition principle. We find that the virialization and phase-mixing times depend on the exponent alpha, with small values of alpha corresponding to longer relaxation times, similarly to what happens when comparing for N-body simulations in classical gravity and in Modified Newtonian Dynamics.

The stability of radially anisotropic spherical stellar systems in modified Newtonian dynamics (MOND) is explored by means of numerical simulations performed with the N-body code N-MODY. We find that Osipkov-Merritt MOND models require for stability larger minimum anisotropy radius than equivalent Newtonian systems (ENSs) with dark matter, and also than purely baryonic Newtonian models with the same density profile. The maximum value for stability of the Fridman-Polyachenko-Shukhman parameter in MOND models is lower than in ENSs, but higher than in Newtonian models with no dark matter. We conclude that MOND systems are substantially more prone to radial-orbit instability than ENSs with dark matter, while they are able to support a larger amount of kinetic energy stored in radial orbits than purely baryonic Newtonian systems. An explanation of these results is attempted, and their relevance to the MOND interpretation of the observed kinematics of globular clusters, dwarf spheroidal and elliptical galaxies is briefly discussed.

By means of a one-zone evolutionary model, we study the co-evolution of supermassive black holes and their host galaxies, as a function of the accretion radiative efficiency, dark matter content, and cosmological infall of gas. In particular, the radiation feedback is computed by using the self-regulated Bondi accretion. The models are characterized by strong oscillations when the galaxy is in the AGN state with a high accretion luminosity. We found that these one-zone models are able to reproduce two important phases of galaxy evolution, namely an obscured-cold phase when the bulk of star formation and black hole accretion occur, and the following quiescent hot phase in which accretion remains highly sub-Eddington. A Compton-thick phase is also found in almost all models, associated with the cold phase. An exploration of the parameter space reveals that the closest agreement with the present-day Magorrian relation is obtained, independently of the dark matter halo mass, for galaxies with a low-mass seed black hole, and the accretion radiative efficiency ~0.1.

The deposition of mechanical feedback from a supermassive black hole (SMBH) in an active galactic nucleus (AGN) into the surrounding galaxy occurs via broad-line winds which must carry mass and radial momentum as well as energy. The effect can be summarized by the dimensionless parameter $\eta=dot{M_outflow}/dot{M_accretion}= (2 \epsilon_w c^2)/v_w^2$ where ($\epslion_w \equiv dot{E}_w/(dot{M_accretion} c^2)$) is the efficiency by which accreted matter is turned into wind energy in the disc surrounding the central SMBH. The outflowing mass and omentum are proportional to $\eta$, and many prior treatments have essentially assumed that $\eta=0$. We perform one- and two-dimensional simulations and find that the growth of the central SMBH is very sensitive to the inclusion of the mass and momentum driving but is insensitive to the assumed mechanical efficiency. For example in representative calculations, the omission of momentum and mass feedback leads to an hundred fold increase in the mass of the SMBH to over $10^{10} \Msun$. When allowance is made for momentum driving, the final SMBH mass is much lower and the wind efficiencies which lead to the most observationally acceptable results are relatively low with $\epsilon_w \lesssim 10^{-4}$.

Following the seminal result of An & Evans, known as the central density slope-anisotropy theorem, successive investigations unexpectedly revealed that the density slope-anisotropy inequality holds not only at the center, but at all radii in a very large class of spherical systems whenever the phase-space distribution function is positive. In this paper we derive a criterion that holds for all spherical systems in which the augmented density is a separable function of radius and potential: this new finding allows to unify all the previous results in a very elegant way, and opens the way for more general investigations. As a first application, we prove that the global density slope-anisotropy inequality is also satisfied by all the explored additional families of multi-component stellar systems. The present results, and the absence of known counter-examples, lead us to conjecture that the global density slope-anisotropy inequality could actually be a universal property of spherical systems with positive distribution function.

We find, from high-resolution hydro simulations, that winds from AGN effectively heat the inner parts (~100 pc) of elliptical galaxies, reducing infall to the central SMBH; and radiative (photoionization and X-ray) heating reduces cooling flows at the kpc scale. Including both types of feedback with (peak) efficiencies of 3 10^{-4} < epsilon_mech < 10^{-3} and of epsilon_rad ~10^{-1.3} respectively, produces systems having duty-cycles, central SMBH masses, X-ray luminosities, optical light profiles, and E+A spectra in accord with the broad suite of modern observations of massive elliptical systems. Our main conclusion is that mechanical feedback (including all three of energy, momentum and mass) is necessary but the efficiency, based on several independent arguments must be a factor of 10 lower than is commonly assumed. Bursts are frequent at z>1 and decline in frequency towards the present epoch as energy and metal rich gas are expelled from the galaxies into the surrounding medium. For a representative galaxy of final stellar mass ~3 10^{11} Msun, roughly 3 10^{10} Msun of recycled gas has been added to the ISM since z~2 and, of that, roughly 63% has been expelled from the galaxy, 19% has been converted into new metal rich stars in the central few hundred parsecs, and 2% has been added to the central SMBH, with the remaining 16% in the form hot X-ray emitting ISM. The bursts occupy a total time of ~170 Myr, which is roughly 1.4% of the available time. Of this time, the central SMBH would be seen as an UV or optical source for ~45% and ~71$% of the time, respectively. Restricting to the last 8.5 Gyr, the burst occupy ~44 Myr, corresponding to a fiducial duty-cycle of ~5 10^{-3}.

By using high-resolution 1D hydrodynamical simulations, we investigate the effects of purely mechanical feedback from super massive black holes (SMBHs) in the evolution of elliptical galaxies for a broad range of feedback efficiencies and compare the results to four major observational constraints. In particular, we focus on 1) the central black hole to stellar mass ratio of the host galaxy, 2) the lifetime of the luminous quasar phase, 3) the mass of stars formed in the host galaxy within the last Gyr, and 4) the X-ray luminosity of the hot diffuse gas. As a result, we try to pin down the most successful range of mechanical feedback efficiencies. We find that while low feedback efficiencies result in too much growth of the SMBH, high efficiencies totally blow out the hot interstellar gas, and the models are characterized by very low thermal X-ray luminosity well below the observed range. The net lifetime of the quasar phase is strongly coupled to the mass ratio between SMBH and its host galaxy, while the X-ray luminosity is generally correlated to the recent star formation within the last Gyr. When considering the popularly adopted model of the constant feedback efficiency, the feedback energy deposited into the ambient medium should be more than 0.01% of the SMBH accretion energy to be consistent with the SMBH mass to stellar mass ratio in the local universe. Yet, the X-ray luminosity of the hot gas favors about 0.005% of the accretion energy as the mechanical AGN feedback energy. We conclude that the purely mechanical feedback mode is unlikely to be simultaneously compatible with all four observable tests, even allowing a broad range of feedback efficiencies, and that including both radiative and mechanical feedback together may be a solution to comply the observational constraints. [abridged]

The analytical generalization of the classical dynamical friction formula (derived under the assumption that all the field particles have the same mass) to the case in which the masses of the field particles are distributed with a mass spectrum is presented. Two extreme cases are considered: in the first, energy equipartition is assumed, in the second all the field particles have the same (Maxwellian) velocity distribution. Three different mass spectra are studied in detail, namely the exponential, discrete (two components), and power-law cases. It is found that the dynamical friction deceleration can be significantly stronger than in the equivalent classical case, with the largest differences (up to a factor of 10 or more in extreme cases) arising for test particle velocities comparable to the mass-averaged velocity dispersion of the field particles. The present results are relevant to our understanding of the dynamical evolution of globular clusters, in particular in the modelization of mass segregation and sedimentation of Blue Straggler stars and Neutron stars, and for the study of binary black holes in galactic nuclei.

Starting from the central density slope-anisotropy theorem of An and Evans (2006), recent investigations have shown that the involved density slope-anisotropy inequality holds not only at the center, but at all radii (i.e. globally) in a very large class of spherical systems with positive phase-space distribution function. Here we present some additional analytical cases that further extend the validity of the global density slope-anisotropy inequality. These new results, several numerical evidences, and the absence of known counter-examples, lead us to conjecture that the global density slope-anisotropy inequality could actually be a universal property of spherical systems with positive distribution function.

It is argued that galaxies like ours sustain their star formation by transferring gas from an extensive corona to the star-forming disc. The transfer is effected by the galactic fountain -- cool clouds that are shot up from the plane to kiloparsec heights above the plane. The Kelvin-Helmholtz instability strips gas from these clouds. If the pressure and the the metallicity of the corona are high enough, the stripped gas causes a similar mass of coronal gas to condense in the cloud's wake. Hydrodynamical simulations of cloud-corona interaction are presented. These confirm the existence of a critical ablation rate above which the corona is condensed, and imply that for the likely parameters of the Galactic corona this rate lies near the actual ablation rate of clouds. In external galaxies trails of HI behind individual clouds will not be detectable, although the integrated emission from all such trails should be significant. Parts of the trails of the clouds that make up the Galaxy's fountain should be observable and may account for features in targeted 21-cm observations of individual high-velocity clouds and surveys of Galactic HI emission. Taken in conjunction with the known decline in the availability of cold infall with increasing cosmic time and halo mass, the proposed mechanism offers a promising explanation of the division of galaxies between the blue cloud to the red sequence in the colour-luminosity plane.

Star-forming disc galaxies such as the Milky Way need to accrete $\gsim$ 1 $M_{\odot}$ of gas each year to sustain their star formation. This gas accretion is likely to come from the cooling of the hot corona, however it is still not clear how this process can take place. We present simulations supporting the idea that this cooling and the subsequent accretion are caused by the passage of cold galactic-fountain clouds through the hot corona. The Kelvin-Helmholtz instability strips gas from these clouds and the stripped gas causes coronal gas to condense in the cloud's wake. For likely parameters of the Galactic corona and of typical fountain clouds we obtain a global accretion rate of the order of that required to feed the star formation.

Both radiative and mechanical feedback from Active Galactic Nuclei have been found to be important for the evolution of elliptical galaxies. We compute how a shock may be driven from a central black hole into the gaseous envelope of an elliptical galaxy by mechanical as well as radiative feedback (in the form of nuclear winds) using high resolution 1-D hydrodynamical simulations. We calculate the synchrotron emission from the electron cosmic rays accelerated by the shocks (not the jets) in the central part of elliptical galaxies, and we also study the synchrotron spectrum's evolution using the standard diffusive shock acceleration mechanism, which is routinely applied to the scaled volume case of supernova remnants. We find good agreement quantitatively between the synchrotron radio emission produced via this mechanism with extant observations of elliptical galaxies which are undergoing outbursts. Additionally, we also find that synchrotron optical and X-ray emission can co-exist inside elliptical galaxies during a certain phase of evolution subsequent to central outbursts. In fact, our calculations predict a synchrotron luminosity of $\sim 1.3\times 10^6 L_{\odot}$ at the frequency 5 GHz (radio band), of $\sim 1.1\times 10^6 L_{\odot}$ at $4.3\times10^{14}$ Hz (R band, corresponding to the absolute magnitude -10.4), and of $\sim 1.5\times 10^{7} L_{\odot}$ at $2.4\times10^{17}$ Hz (soft X-ray, 0.5 -- 2.0 keV band).

The kinematics of the extra-planar neutral and ionised gas in disc galaxies shows a systematic decline of the rotational velocity with height from the plane (vertical gradient). This feature is not expected for a barotropic gas, whilst it is well reproduced by baroclinic fluid homogeneous models. The problem with the latter is that they require gas temperatures (above $10^5$ K) much higher than the temperatures of the cold and warm components of the extra-planar gas layer. In this paper, we attempt to overcome this problem by describing the extra-planar gas as a system of gas clouds obeying the Jeans equations. In particular, we consider models having the observed extra-planar gas distribution and gravitational potential of the disc galaxy NGC 891: for each model we construct pseudo-data cubes and we compare them with the HI data cube of NGC 891. In all cases the rotational velocity gradients are in qualitative agreement with the observations, but the synthetic and the observed data cubes of NGC 891 show systematic differences that cannot be accommodated by any of the explored models. We conclude that the extra-planar gas in disc galaxies cannot be satisfactorily described by a stationary Jeans-like system of gas clouds.

General criteria to check the positivity of the distribution function (phase-space consistency) of stellar systems of assigned density and anisotropy profile are useful starting points in Jeans-based modeling. Here we substantially extend previous results, and we present the inversion formula and the analytical necessary and sufficient conditions for phase-space consistency of the family of multi-component Cuddeford spherical systems: the distribution function of each density component of these systems is defined as the sum of an arbitrary number of Cuddeford distribution functions with arbitrary values of the anisotropy radius, but identical angular momentum exponent. The radial trend of anisotropy that can be realized by these models is therefore very general. As a surprising by-product of our study, we found that the ``central cusp-anisotropy theorem'' (a necessary condition for consistency relating the values of the central density slope and of the anisotropy parameter) holds not only at the center, but at all radii in consistent multi-component generalized Cuddeford systems. This last result suggests that the so--called mass--anisotropy degeneracy could be less severe than what is sometimes feared.

The effects of dry and wet merging on the Scaling Laws (SLs) of elliptical galaxies (Es) are discussed. It is found that the galaxy SLs, possibly established at high redshift by the fast collapse of gas-rich and clumpy stellar distributions in preexisting dark matter halos following the cosmological SLs, are compatible with a (small) number of galaxy mergers at lower redshift.

The passively evolving stellar population in elliptical galaxies (Es) provides a continuous source of fuel for accretion on the central supermassive black hole (SMBH), which is 1) extended over the entire galaxy life (but declining with cosmic time), 2) linearly proportional to the stellar mass of the host spheroid, 3) summing up to a total gas mass that is >100 times larger than the currently observed SMBH masses, 4) available independently of merging events. The main results of numerical simulations of Es with central SMBH, in which a physically based implementation of radiative and mechanical feedback effects is considered, are presented.

The centers of elliptical galaxies host supermassive black holes that significantly affect the surrounding interstellar medium through feedback resulting from the accretion process. The evolution of this gas and of the nuclear emission during the galaxies' lifetime has been studied recently with high-resolution hydrodynamical simulations. These included gas cooling and heating specific for an average AGN spectral energy distribution, a radiative efficiency declining at low mass accretion rates, and mechanical coupling between the hot gas and AGN winds. Here we present a short summary of the observational properties resulting from the simulations, focussing on 1) the nuclear luminosity; 2) the global luminosity and temperature of the hot gas; 3) its temperature profile and X-ray brightness profile. These properties are compared with those of galaxies of the local universe, pointing out the successes of the adopted feedback and the needs for new input in the simulations.

The importance of the radiative feedback from SMBHs at the centers of elliptical galaxies is not in doubt, given the well established relations among electromagnetic output, black hole mass and galaxy optical luminosity. In addition, feedback due to mechanical and thermal deposition of energy from jets and winds emitted by the accretion disk around the central SMBH is also expected to occur. In this paper we improve and extend the accretion and feedback physics explored in our previous papers to include also a physically motivated mechanical feedback. We study the evolution of an isolated elliptical galaxy with the aid of a high-resolution 1-D hydrodynamical code, where the cooling and heating functions include photoionization and Compton effects, and restricting to models which include only radiative or only mechanical feedback. We confirm that for Eddington ratios above 0.01 both the accretion and radiative output are forced by feedback effects to be in burst mode, so that strong intermittencies are expected at early times, while at low redshift the explored models are characterized by smooth, very sub-Eddington mass accretion rates punctuated by rare outbursts. However, the explored models always fail some observational tests. If we assume the high mechanical efficiency of 10^{-2.3}, we find that most of the gas is ejected from the galaxy, the resulting X-ray luminosity is far less than is typically observed and little SMBH growth occurs. But models with low enough mechanical efficiency to accomodate satisfactory SMBH growth tend to allow too strong cooling flows and leave galaxies at z=0 with E+A spectra more frequently than is observed. We conclude that both types of feedback are required. Models with combined feedback are explored in a forthcoming paper [abridged]

We apply the joint lensing and dynamics code for the analysis of early-type galaxies, CAULDRON, to a rotating N-body stellar system with dark matter halo which significantly violates the two major assumptions of the method, i.e. axial symmetry supported by a two-integral distribution function. The goal is to study how CAULDRON performs in an extreme case, and to determine which galaxy properties can still be robustly recovered. Three data sets, corresponding to orthogonal lines of sight, are generated from the N-body system and analysed with the identical procedure followed in the study of real lens galaxies, adopting an axisymmetric power-law total density distribution. We find that several global properties of the N-body system are recovered with remarkable accuracy, despite the fact that the adopted power-law model is too simple to account for the lack of symmetry of the true density distribution. In particular, the logarithmic slope of the total density distribution is robustly recovered to within less than 10 per cent (with the exception of the ill-constrained very inner regions), the inferred angle-averaged radial profile of the total mass closely follows the true distribution, and the dark matter fraction of the system (inside the effective radius) is correctly determined within ~ 10 per cent of the total mass. Unless the line of sight direction is almost parallel to the total angular momentum vector of the system, reliably recovered quantities also include the angular momentum, the V/sigma ratio, and the anisotropy parameter delta. We conclude that the CAULDRON code can be safely and effectively applied to real early-type lens galaxies, providing reliable information also for systems that depart significantly from the method's assumptions.

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

We describe some results obtained with N-MODY, a code for N-body simulations of collisionless stellar systems in modified Newtonian dynamics (MOND). We found that a few fundamental dynamical processes are profoundly different in MOND and in Newtonian gravity with dark matter. In particular, violent relaxation, phase mixing and galaxy merging take significantly longer in MOND than in Newtonian gravity, while dynamical friction is more effective in a MOND system than in an equivalent Newtonian system with dark matter.

It is well known that the density and anisotropy profile in the inner regions of a stellar system with positive phase-space distribution function are not fully independent. Here we study the interplay between density profile and orbital anisotropy at large radii in physically admissible (consistent) stellar systems. The analysis is carried out by using two-component (n-gamma,gamma_1) spherical self-consistent galaxy models, in which one density distribution follows a generalized gamma profile with external logarithmic slope n, and the other a standard gamma_1 profile (with external slope 4). The necessary and sufficient conditions for phase-space consistency are determined analytically, also in presence of a dominant massive central black hole, and the analytical phase-space distribution function of (n-gamma,1) models, and of n-gamma models with a central black hole, is derived for gamma=0,1,2. It is found that the density slope in the external regions of a stellar system can play an important role in determining the amount of admissible anisotropy: in particular, for fixed density slopes in the central regions, systems with a steeper external density profile can support more radial anisotropy than externally flatter models. This is quantified by an inequality formally identical to the ``cusp slope-central anisotropy'' theorem (An & Evans 2006), relating at all radii (and not just at the center) the density logarithmic slope and the anisotropy indicator in all Osipkov-Merritt systems [Abridged].

Dynamical properties of two-component galaxy models whose stellar density distribution is described by a gamma-model while the total density distribution has a pure r^(-2) profile, are presented. The orbital structure of the stellar component is described by Osipkov-Merritt anisotropy, while the dark matter halo is isotropic. After a description of minimum halo models, the positivity of the phase-space density (the model consistency) is investigated, and necessary and sufficient conditions for consistency are obtained analytically as a function of the stellar inner density slope gamma and anisotropy radius. The explicit phase-space distribution function is recovered for integer values of gamma, and it is shown that while models with gamma>4/17 are consistent when the anisotropy radius is larger than a critical value (dependent on gamma), the gamma=0 models are unphysical even in the fully isotropic case. The Jeans equations for the stellar component are then solved analytically; in addition, the projected velocity dispersion at the center and at large radii are also obtained analytically for generic values of the anisotropy radius, and it is found that they are given by remarkably simple expressions. The presented models, even though highly idealized, can be useful as starting point for more advanced modeling of the mass distribution of elliptical galaxies in studies combining stellar dynamics and gravitational lensing.

After the discovery that supermassive black holes (SMBHs) are ubiquitous at the center of stellar spheroids and that their mass Mbh, in the range 10^6-10^9 Msun, is tightly related to global properties of the host stellar system, the idea of the co-evolution of elliptical galaxies and of their SMBHs has become a central topic of modern astrophysics. Here, I summarize some consequences that can be derived from the galaxy Scaling Laws (SLs) and present a coherent scenario for the formation and evolution of elliptical galaxies and their central SMBHs, focusing in particular on the establishment and maintenance of their SLs. In particular, after a first observationally based part, the discussion focuses on the physical interpretation of the Fundamental Plane. Then, two important processes in principle able to destroy the galaxy and SMBH SLs, namely galaxy merging and cooling flows, are analyzed. Arguments supporting the necessity to clearly distinguish between the origin and maintenance of the different SLs, and the unavoidable occurrence of SMBH feedback on the galaxy ISM in the late stages of galaxy evolution (when elliptical galaxies are sometimes considered as ``dead, red objects''), are then presented. At the end of the paper I will discuss some implications of the recent discovery of super-dense ellipticals in the distant Universe. In particular, I will argue that, if confirmed, these new observations would lead to the conclusion that at early epochs a relation between the stellar mass of the galaxy and the mass of the central SMBH should hold, consistent with the present day Magorrian relation, while the proportionality coefficient between Mbh and the scale of velocity dispersion of the hosting spheroids should be significantly smaller than that at the present epoch.

Estimating the temperature and metal abundance of the intracluster and the intragroup media is crucial to determine their global metal content and to determine fundamental cosmological parameters. When a spatially resolved temperature or abundance profile cannot be recovered from observations (e.g., for distant objects), or deprojection is difficult (e.g., due to a significant non-spherical shape), only global average temperature and abundance are derived. After introducing a general technique to build hydrostatic gaseous distributions of prescribed density profile in potential wells of any shape, we compute the global mass weighted and emission weighted temperature and abundance for a large set of barotropic equilibria and an observationally motivated abundance gradient. We also compute the spectroscopic-like temperature that is recovered from a single temperature fit of observed spectra. The derived emission weighted abundance and temperatures are higher by 50% to 100% than the corresponding mass weighted quantities, with overestimates that increase with the gas mean temperature. Spectroscopic temperatures are intermediate between mass and luminosity weighted temperatures. Dark matter flattening does not lead to significant differences in the values of the average temperatures or abundances with respect to the corresponding spherical case (except for extreme cases).

In a previous paper the complex-shift method has been applied to self-gravitating spherical systems, producing new analytical axisymmetric density-potential pairs. We now extend the treatement to the Miyamoto-Nagai disc and to the Binney logarithmic halo, and we study the resulting axisymmetric and triaxial analytical density-potential pairs; we also show how to obtain the surface density of shifted systems from the complex-shift of the surface density of the parent model. In particular, the systems obtained from Miyamoto-Nagai discs can be used to describe disc galaxies with a peanut-shaped bulge or with a central triaxial bar, depending on the direction of the shift vector. By using a constructive method that can be applied to generic axisymmetric systems, we finally show that the Miyamoto-Nagai and the Satoh discs, and the Binney logarithmic halo, cannot be obtained from the complex-shift of any spherical parent distribution. As a by-product of this study we also found two new generating functions in closed form for even and odd Legendre polynomials, respectively.

We have tested a previous analytical estimate of the dynamical friction timescale in Modified Newtonian Dynamics (MOND) with fully non-linear N-body simulations. The simulations confirm that the dynamical friction timescale is significantly shorter in MOND than in equivalent Newtonian systems, i.e. systems with the same phase-space distribution of baryons and additional dark matter. An apparent conflict between this result and the long timescales determined for bars to slow and mergers to be completed in previous N-body simulations of MOND systems is explained. The confirmation of the short dynamical-friction timescale in MOND underlines the challenge that the Fornax dwarf spheroidal poses to the viability of MOND.

We computed the chemical evolution of spiral bulges hosting Seyfert nuclei, based on updated chemical and spectro-photometrical evolution models for the bulge of our Galaxy, made predictions about other quantities measured in Seyferts, and modeled the photometry of local bulges. The chemical evolution model contains detailed calculations of the Galactic potential and of the feedback from the central supermassive black hole, and the spectro-photometric model covers a wide range of stellar ages and metallicities. We followed the evolution of bulges in the mass range 10^9 - 10^{11} Msun by scaling the star formation efficiency and the bulge scalelength as in the inverse-wind scenario for elliptical galaxies, and considering an Eddington limited accretion onto the central supermassive black hole. We successfully reproduced the observed black hole-host bulge mass relation. The observed nuclear bolometric luminosity is reproduced only at high redshift or for the most massive bulges; in the other cases, at z = 0 a rejuvenation mechanism is necessary. The black hole feedback is in most cases not significant in triggering the galactic wind. The observed high star formation rates and metal overabundances are easily achieved, as well as the constancy of chemical abundances with redshift and the bulge present-day colours. Those results are not affected if we vary the index of the stellar IMF from x=0.95 to x=1.35; a steeper IMF is instead required in order to reproduce the colour-magnitude relation and the present K-band luminosity of the bulge.

We computed the chemical evolution of Seyfert galaxies, residing in spiral bulges, based on an updated model for the Milky Way bulge with updated calculations of the Galactic potential and of the feedback from the central supermassive black hole (BH) in a spherical approximation. We followed the evolution of bulges of masses $2\times 10^{9}-10^{11}M_{\odot}$ by scaling the star-formation efficiency and the bulge scalelenght as in the inverse-wind scenario for ellipticals. We successfully reproduced the observed relation between the BH mass and that of the host bulge, and the observed peak nuclear bolometric luminosity. The observed metal overabundances are easily achieved, as well as the constancy of chemical abundances with the redshift.

We present the results of N-body simulations of dissipationless galaxy merging in Modified Newtonian Dynamics (MOND). For comparison, we also studied Newtonian merging between galaxies embedded in dark matter halos, with internal dynamics equivalent to the MOND systems. We found that the merging timescales are significantly longer in MOND than in Newtonian gravity with dark matter, suggesting that observational evidence of rapid merging could be difficult to explain in MOND. However, when two galaxies eventually merge, the MOND merging end-product is hardly distinguishable from the final stellar distribution of an equivalent Newtonian merger with dark matter.

The importance of the radiative feedback from massive black holes at the centers of elliptical galaxies is not in doubt, given the well established relations among electromagnetic output, black hole mass and galaxy optical luminosity. We show how this AGN radiative output affects the hot ISM of an isolated elliptical galaxy with the aid of a high-resolution hydrodynamical code, where the cooling and heating functions include photoionization plus Compton heating. We find that radiative heating is a key factor in the self-regulated coevolution of massive black holes and their host galaxies and that 1) the mass accumulated by the central black hole is limited by feedback to the range observed today, and 2) relaxation instabilities occur so that duty cycles are small enough (~0.03) to account for the very small fraction of massive ellipticals observed to be in the "on" -QSO- phase, when the accretion luminosity approaches the Eddington luminosity. The duty cycle of the hot bubbles inflated at the galaxy center during major accretion episodes is of the order of 0.1-0.4. Major accretion episodes caused by cooling flows in the recycled gas produced by normal stellar evolution trigger nuclear starbursts coincident with AGN flaring. During such episodes the central sources are often obscured; but overall, in the bursting phase (1<z<3), the duty cycle of the black hole in its "on" phase is of the order of percents and it is unobscured approximately one-third of the time. Mechanical energy output from non-relativistic gas winds integrates to 2.3 10^{59} erg, with most of it caused by broadline AGN outflows. [abridged]

We investigate the possibility of discriminating between Modified Newtonian Dynamics (MOND) and Newtonian gravity with dark matter, by studying the vertical dynamics of disk galaxies. We consider models with the same circular velocity in the equatorial plane (purely baryonic disks in MOND and the same disks in Newtonian gravity embedded in spherical dark matter haloes), and we construct their intrinsic and projected kinematical fields by solving the Jeans equations under the assumption of a two-integral distribution function. We found that the vertical velocity dispersion of deep-MOND disks can be much larger than in the equivalent spherical Newtonian models. However, in the more realistic case of high-surface density disks this effect is significantly reduced, casting doubts on the possibility of discriminating between MOND and Newtonian gravity with dark matter by using current observations.

Dissipationless collapses in Modified Newtonian Dynamics (MOND) have been studied by using our MOND particle-mesh N-body code, finding that the projected density profiles of the final virialized systems are well described by Sersic profiles with index m<4 (down to m~2 for a deep-MOND collapse). The simulations provided also strong evidence that phase mixing is much less effective in MOND than in Newtonian gravity. Here we describe "ad hoc" numerical simulations with the force angular components frozen to zero, thus producing radial collapses. Our previous findings are confirmed, indicating that possible differences in radial orbit instability under Newtonian and MOND gravity are not relevant in the present context.

Dissipationless collapses in Modified Newtonian Dynamics (MOND) are studied by using a new particle-mesh N-body code based on our numerical MOND potential solver. We found that low surface-density end-products have shallower inner density profile, flatter radial velocity-dispersion profile, and more radially anisotropic orbital distribution than high surface-density end-products. The projected density profiles of the final virialized systems are well described by Sersic profiles with index m~4, down to m~2 for a deep-MOND collapse. Consistently with observations of elliptical galaxies, the MOND end-products, if interpreted in the context of Newtonian gravity, would appear to have little or no dark matter within the effective radius. However, we found impossible (under the assumption of constant mass-to-light ratio) to simultaneously place the resulting systems on the observed Kormendy, Faber-Jackson and Fundamental Plane relations of elliptical galaxies. Finally, the simulations provide strong evidence that phase mixing is less effective in MOND than in Newtonian gravity.

We show how the complex-shift method introduced by Appell in gravity to the case of a point mass (and applied among others in electrodynamics by Newman, Carter, Lynden-Bell, and Kaiser to determine remarkable properties of the electromagnetic field of rotating charged configurations), can be extended to obtain new and explicit density-potential pairs for self-gravitating systems departing significantly from spherical symmetry. The rotational properties of two axysimmetric baroclinic gaseous configurations derived with the proposed method are illustrated.

With the aid of a simple yet robust approach we investigate the influence of dissipationless and dissipative merging on galaxy structure, and the consequent effects on the scaling laws followed by elliptical galaxies. Our results suggest that ellipticals cannot be originated by parabolic merging of low mass spheroids only, even in presence of substantial gas dissipation. However, we also found that scaling laws such as the Faber-Jackson, Kormendy, Fundamental Plane, and the Mbh -sigma relations, when considered over the whole mass range spanned by ellipticals in the local universe, are robust against merging. We conclude that galaxy scaling laws, possibly established at high redshift by the fast collapse in pre-existing dark matter halos of gas rich and clumpy stellar distributions, are compatible with a (small) number of galaxy mergers at lower redshift.

NGC3379 is a well studied nearby elliptical for which optical investigations have claimed a little dark matter content, or even no dark matter. Recently, its total mass profile M(r) has been derived by exploiting Chandra observations of its extended and X-ray emitting interstellar medium, based on the hypothesis of hydrostatic equilibrium for the hot gas. The resulting total mass within the effective radius Re has been claimed to be a few times larger than found by optical studies. Here we show that part of the discrepancy can be due to an underestimate of the optically derived mass, and the remaining discrepancy of a factor of ~2 can be explained by deviations from hydrostatic equilibrium of the hot gas. By using hydrodynamical simulations tailored to reproduce the observed hot gas properties of NGC3379, and by assuming as input for the simulations the total mass profile derived optically, we show that i) the hot gas at the present time has X-ray properties consistent with those observed only if it is outflowing over most of the galactic body, and ii) an overestimate of M of the same size found in the recent X-ray analysis is recovered when assuming hydrostatic equilibrium. We also show that the hot gas is outflowing even for a dark matter fraction within Re as large as derived with the standard X-ray procedure based on the hydrostatic equilibrium assumption, which shows the unapplicability of the method for this galaxy. Finally, we find that the whole range of dark mass amount and distribution allowed for by optical studies is compatible with a hot gas flow with the observed X-ray properties.

By means of high-resolution N-body simulations we revisited the dissipationless collapse scenario for galaxy formation. We considered both single-component collapses and collapses of a cold stellar distribution in a live dark matter halo. Single-component collapses lead to stellar systems whose projected profiles are fitted very well by the Sersic R^(1/m) law with 3.6 < m < 8. The stellar end-products of collapses in a dark matter halo are still well described by the R^(1/m) law, but with 1.9 < m < 12, where the lowest m values are obtained when the halo is dominant. In all the explored cases the profiles at small radii deviate from their global best-fit R^(1/m) model, being significantly flatter. The break-radius values are comparable with those measured in `core' elliptical galaxies, and are directly related to the coldness of the initial conditions. The dissipationless collapse of initially cold stellar distributions in pre-existing dark matter haloes may thus have a role in determining the observed weak homology of elliptical galaxies.

We present a simple method, based on the deformation of spherically symmetric potentials, to construct explicit axisymmetric and triaxial MOND density-potential pairs. General guidelines to the choice of suitable deformations, so that the resulting density distribution is nowhere negative, are presented. This flexible method offers for the first time the possibility to study the MOND gravitational field for sufficiently general and realistic density distributions without resorting to sophisticated numerical codes. The technique is illustrated by constructing the MOND density-potential pair for a triaxial galaxy model that, in the absence of deformation, reduces to the Hernquist sphere. Such analytical solutions are also relevant to test and validate numerical codes. Here we present a new numerical potential solver designed to solve the MOND field equation for arbitrary density distributions: the code is tested with excellent results against the analytic MOND triaxial Hernquist model and the MOND razor-thin Kuzmin disk, and a simple application is finally presented.

We show that fluid stationary models are able to reproduce the observed, negative vertical gradient of the rotation velocity of the extra-planar gas in spiral galaxies. We have constructed models based on the simple condition that the pressure of the medium does not depend on density alone (baroclinic instead of barotropic solutions: isodensity and isothermal surfaces do not coincide). As an illustration, we have successfully applied our method to reproduce the observed velocity gradient of the lagging gaseous halo of NGC 891. The fluid stationary models discussed here can describe a hot homogeneous medium as well as a "gas" made of discrete, cold HI clouds with an isotropic velocity dispersion distribution. Although the method presented here generates a density and velocity field consistent with observational constraints, the stability of these configurations remains an open question.

The contribution of ordered rotation to the observed tilt and thickness of the Fundamental Plane of elliptical galaxies is studied by means of oblate, two-integrals cuspy galaxy models with adjustable flattening, variable amount of ordered rotational support, and possible presence of a dark matter halo and of a central super-massive black hole. We find that, when restricting the measure of the velocity dispersion to the central galactic regions, rotation has a negligible effect, and so cannot be responsible of the observed tilt. However, streaming velocity effect can be significant when observing small and rotationally supported galaxies through large (relative) aperture (as for example in Fundamental Plane studies at high redshift), and can lead to unrealistically low mass-to-light ratios. The effect of a central supermassive black hole on the kinematical fields, and the models position in the v/sigma-ellipticity plane are also discussed.

We have performed VLT spectroscopy of an almost complete sample of 18 early-type galaxies with 0.88 < z < 1.3 plus two at z=0.67, selected from the K20 survey, and derived the velocity dispersion for 15+2 of them. By combining these data with HST and VLT images, we study the Fundamental Plane (FP), the Faber-Jackson and the Kormendy relations at z~1, and compare them with the local ones. The FP at z~1 keeps a remarkably small scatter, and shows both an offset and a rotation, which we interpret in terms of evolution of the mass-to-light ratio, and possibly of the size. We show evidence that the evolution rate depends on galaxy mass, being faster for less massive galaxies. We discuss the possible factors driving the evolution of spheroids and compare our results with the predictions of the hierachical models of galaxy formation.

We start from a study of the density-potential relation for classical homeoids in terms of an asymptotic expansion for small deviations from spherical symmetry. We then show that such expansion is a useful device that allows us to construct a variety of exact density-potential pairs with spheroidal, toroidal, or triaxial shapes for which the deviation from spherical symmetry is finite. As concrete analytical applications, we describe: (1) The construction of a family of toroidal axisymmetric density-potential pairs one of which is associated with a perfectly flat rotation curve (for a member of this family, the supporting two-integral phase-space distribution function is obtained in closed form); (2) The determination of the aperture velocity dispersion in a wide class of two-integral axisymmetric models not stratified on homeoids with central black hole, which may be useful for the discussion of the dynamical contributions to the characteristics of the Fundamental Plane of early-type galaxies; and (3) For such class of models, the construction of the v/sigma-ellipticity relation, often considered to assess the role of rotation in the structure of elliptical galaxies.

We discuss the importance of feedback via photoionization and Compton heating on the co-evolution of massive black holes (MBHs) at the center of spheroidal galaxies and their stellar and gaseous components. We first assess the energetics of the radiative feedback from a typical quasar on the ambient interstellar medium. We then demonstrate that the observed Mbh-sigma relation could be established following the conversion of most of the gas of an elliptical progenitor into stars, specifically when the gas-to-stars mass ratio in the central regions has dropped to a low level ~0.01 or less, so that gas cooling is no longer able to keep up with the radiative heating by the growing central MBH. A considerable amount of the remaining gas will be expelled and both MBH accretion and star formation will proceed at significantly reduced rates thereafter, in agreement with observations of present day ellipticals. We find further support for this scenario by evolving over an equivalent Hubble time a simple, physically based toy model that additionally takes into account the mass and energy return for the spheroid evolving stellar population, a physical ingredient often neglected in similar approaches.

We discuss the role of feedback via photoionization and Compton heating in the co-evolution of massive black holes at the center of spheroidal galaxies and their stellar and gaseous components. We first assess the energetics of the radiative feedback from a typical quasar on the ambient interstellar gas. We then demonstrate that the observed Mbh-sigma relation could be established at a relatively early epoch in galactic evolution when the formation of the stellar bulge was almost completed and the gas-to-stars mass ratio was reduced to a low level ~0.01 such that cooling could not keep up with radiative heating. A considerable amount of gas was expelled at that time and black hole accretion proceeded at a much lower rate thereafter.

There is abundant evidence that heating processes in the central regions of elliptical galaxies has both prevented large-scale cooling flows and assisted in the expulsion of metal rich gas. We now know that each such spheroidal system harbors in its core a massive black hole weighing approximately 0.13% of the mass in stars and also know that energy was emitted by each of these black holes with an efficiency exceeding 10% of its rest mass. Since, if only 0.5% of that radiant energy were intercepted by the ambient gas, its thermal state would be drastically altered, it is worth examining in detail the interaction between the out-flowing radiation and the equilibrium or inflowing gas. On the basis of detailed hydrodynamic computations we find that relaxation oscillations are to be expected with the radiative feedback quite capable of regulating both the growth of the central black hole and also the density and thermal state of the gas in the galaxy. Mechanical input of energy by jets may assist or dominate over these radiative effects. We propose specific observational tests to identify systems which have experienced strong bursts of radiative heating from their central black holes.

We show how to construct families of stationary hydrodynamical configurations that reproduce the observed vertical gradient of the rotation velocity of the extra-planar gas in spiral galaxies. We then present a simple model for the lagging halo of the spiral galaxy NGC 891, which is in agreement with the HI observations. Our method is based on well known properties of baroclinic solutions, and it is an elementary application of a much more general and flexible method.

We present a summary of the results obtained with a time-dependent, one-zone toy model aimed at exploring the importance of radiative feedback on the co-evolution of massive black holes (MBHs) at the center of stellar spheroids and their stellar and gaseous components. We consider cosmological infall of gas as well as the mass and energy return for the evolving stellar population. The AGN radiative heating and cooling are described by assuming photoionization equilibrium of a plasma interacting with the average quasar SED. Our results nicely support a new scenario in which the AGN accretion phase characterized by a very short duty-cycle (and now common in the Universe) is due to radiative feedback. The establishment of this phase is recorded as a fossil in the Magorrian and Mbh-sigma relations.

Using N-body simulations we show that the process of formation of the brightest cluster galaxy through dissipationless galactic cannibalism can affect the inner cluster dark matter density profile. In particular, we use as realistic test case the dynamical evolution of the galaxy cluster C0337-2522 at redshift z=0.59, hosting in its centre a group of five elliptical galaxies which are likely to be the progenitor of a central giant elliptical. After the formation of the brightest cluster galaxy, the inner cluster dark matter density profile is significantly flatter (logarithmic slope 0.49<beta<0.90) than the original cusp (beta=1), as a consequence of dynamical friction heating of the massive galaxies against the diffuse cluster dark matter. In our simulations we have assumed that the cluster galaxies are made of stars only. We also show that the presence of galactic dark matter haloes can steepen the cluster central density profile. We conclude that galactic cannibalism could be a viable physical mechanism to reconcile - at least at the cluster scale - the flat dark matter haloes inferred observationally in some galaxy clusters with the steep haloes predicted by cosmological simulations.

By means of simple numerical models we discuss whether "collisionless stellar evaporation" from cluster elliptical galaxies could be an effective mechanism for the production of intracluster stellar populations. The effectiveness of this mechanism is due to the fact that, for realistic galaxy and cluster models, the galaxy oscillation periods near equilibrium configurations in the cluster tidal field are of the same order of stellar orbital times in the external parts of the galaxies themselves. With the aid of Monte-Carlo simulations we explore the evolution of stellar orbits in oscillating galaxies placed near different equilibrium positions. We found that, over an Hubble time, the main effect is a substantial expansion of the galactic outskirts, particularly affecting the galaxy at the cluster center and those orbiting near the cluster core radius: overall, approximately the 10% of the galaxy mass is affected. Thus, the proposed mechanism could be of some importance in the shaping of the halo of cD galaxies and in making easier ``galaxy harassment'' in the formation of the intracluster stellar population.

A naive extension to MOND of the standard computation of the two-body relaxation time Tb implies that Tb is comparable to the crossing time regardless of the number N of stars in the system. This computation is questionable in view of the non-linearity of MOND's field equation. A non-standard approach to the calculation of Tb is developed that can be extended to MOND whenever discreteness noise generates force fluctuations that are small compared to the mean-field force. It is shown that this approach yields standard Newtonian results for systems in which the mean density profile is either plane-parallel or spherical. In the plane-parallel case we find that in the deep-MOND regime Tbb scales with N as in the Newtonian case, but is shorter by the square of the factor by which MOND enhances the gravitational force over its Newtonian value for the same system. Application of these results to dwarf galaxies and groups and clusters of galaxies reveals that in MOND luminosity segregation should be far advanced in groups and clusters of galaxies, two body relaxation should have substantially modified the density profiles of galaxy groups, while objects with masses in excess of ~10M_sun should have spiralled to the centres of dwarf galaxies.(Abridged)

The X-ray emitting interstellar medium of early type galaxies is often used as a tool to determine their total mass M and stellar orbital anisotropy beta profiles, based on the hypothesis of hydrostatic equilibrium for the hot gas. Here we investigate the effects that deviations from equilibrium have on M and beta estimates, by using simple analytical calculations and hydrodynamical simulations representative of gas rich galaxies. We show that the deviations of the X-ray determined beta^{est} and M^{est} from the true values are linked by a remarkably simple relation; in particular, M is underestimated if beta is overestimated, and more radially anisotropic orbital distributions than true are deduced in presence of gas infall velocities of the order of the local stellar velocity dispersion (as are likely in the central regions of galactic cooling flows). The results of this analysis are applied to the most thoroughly investigated bright elliptical, NGC4472. We show that beta^{est} recently derived from X-rays corresponds to a galaxy that is unstable by radial orbit instability. Then, assuming as true beta and M the optically derived values, we show that the differences beta^{est}-beta and M^{est}-M agree with the predictions found here in the case of lack of hydrostatic equilibrium, which points to the latter as a possible explanation for the discrepancies (abridged).

In standard treatments the Eddington luminosity is calculated by assuming that the electron-photon cross section is well described by the Thomson cross section which is gray (frequency independent). Here we discuss some consequence of the introduction of free-free opacity in the Eddington luminosity computation: in particular, due to the dependence of free-free emission on the square of the gas density, it follows that the associated absorption cross section increases linearly with the gas density, so that in high density environments Eddington luminosity is correspondingly reduced. We present a summary of an ongoing exploration of the parameter space of the problem, and we conclude that Eddington luminosity in high density environments can be lowered by a factor of ten or more, making it considerably easier for black holes to accelerate and eject ambient gas.

We present preliminary results on a new family of distribution functions that are able to generate axisymmetric, truncated (i.e., finite size) stellar dynamical models characterized by toroidal shapes. The relevant distribution functions generalize those that are known to describe polytropic spheres, for which all the dynamical and structural properties of the system can be expressed in explicit form as elementary functions of the system gravitational potential. The model construction is then completed by a numerical study of the associated Poisson equation. We note that our axisymmetric models can also include the presence of an external gravitational field, such as that produced by a massive disk or by a central mass concentration (e.g., a supermassive black hole).

Recent observations by Chandra and XMM-Newton demonstrate that the central gas in "cooling flow" galaxy clusters has a mass cooling rate that decreases rapidly with decreasing temperature. This contrasts the predictions of a steady state cooling flow model. On the basis of these observational results, the gas can be in a steady state only if a steady temperature dependent heating mechanism is present; alternatively the gas could be in an unsteady state, i.e., heated intermittently. Intermittent heating can be produced by accretion on the supermassive black hole residing in the central cluster galaxy, via Compton heating. This mechanism can be effective provided that the radiation temperature of the emitted spectrum is higher than the gas temperature. Here we explore whether this heating mechanism can be at the origin of the enigmatic behavior of the hot gas in the central regions of ``cooling flow'' clusters. Although several characteristics of Compton heating appear attractive in this respect, we find that the fraction of absorbed heating for realistic gas and radiation temperatures falls short by two orders of magnitude of the required heating.

By means of high-resolution N-body simulations in a LambdaCDM cosmology, we verify that scaling relations similar to those observed for nearby galaxy clusters are also defined by their dark matter hosts; the slopes, however, are not the same. We then show that the scaling relations of galaxy clusters can be explained as the result of the cosmological collapse of density fluctuations at the appropriate scales, plus a systematic trend of the M/L ratio with cluster mass. The empirical fact that the exponent of the Faber-Jackson relation of elliptical galaxies is significantly different (higher) than that of clusters, force us to conclude that the galaxy scaling laws might derive from the cosmological collapse of density fluctuations at the epoch when galactic scales became non-linear, plus modifications afterward due to early-time dissipative merging.

The X-ray cluster C0337-2522 at redshift $z=0.59$ hosts in its core a group of five elliptical galaxies. Using N-body simulations we show that a multiple merging event among the five galaxies is expected to take place in the next few Gyrs, forming a central brightest cluster galaxy. We also find indications that dynamical friction heating associated with this event is likely to modify the central slope of the cluster dark matter density profile.

Like early-type galaxies, also nearby galaxy clusters define a Fundamental Plane, a luminosity-radius, and a luminosity-velocity dispersion relations, whose physical origin is still unclear. By means of high resolution N--body simulations of massive dark matter halos in a Lambda-CDM cosmology, we find that scaling relations similar to those observed for galaxy clusters are already defined by their dark matter hosts. The slopes however are not the same, and among the various possibilities in principle able to bring the simulated and the observed scaling relations in mutual agreement, we show that the preferred solution is a luminosity dependent mass-to-light ratio (M/L ~ L^0.3), that well corresponds to what inferred observationally. We then show that at galactic scales there is a conflict between the cosmological predictions of structure formation, the observed trend of the mass-to-light ratio in ellipticals, and the slope of their luminosity-velocity dispersion relation (that significantly differs from the analogous one followed by clusters). The conclusion is that the scaling laws of elliptical galaxies might be the combined result of the cosmological collapse of density fluctuations at the epoch when galactic scales became non-linear, plus important modifications afterward due to early-time dissipative merging. Finally, we briefly discuss the possible evolution of the cluster scaling relations with redshift.

According to the galactic cannibalism model, cD galaxies are formed in the center of galaxy clusters by merging of massive galaxies and accretion of smaller stellar systems: however, observational examples of the initial phases of this process are lacking. We have identified a strong candidate for this early stage of cD galaxy formation: a group of five elliptical galaxies in the core of the X-ray cluster C0337-2522 at redshift z=0.59. With the aid of numerical simulations, in which the galaxies are represented by N-body systems, we study their dynamical evolution up to z=0; the cluster dark matter distribution is also described as a N-body system. We find that a multiple merging event in the considered group of galaxies will take place before z=0 and that the merger remnant preserves the Fundamental Plane and the Faber-Jackson relations, while its behavior with respect to the Mbh-sigma relation is quite sensitive to the details of black hole merging [abridged].

We discuss a simple model for the growth of supermassive black holes (BHs) at the center of spheroidal stellar systems. In particular, we assess the hypotheses that (1) star formation in spheroids and BH fueling are proportional to one another, and (2) the BH accretion luminosity stays near the Eddington limit during luminous quasar phases. With the aid of this simple model, we are able to interpret many properties of the QSO luminosity function, including the puzzling steep decline of the characteristic luminosity from redshift z=2 to to z=0: indeed the residual star formation in spheroidal systems is today limited to a small number of bulges, characterized by stellar velocity dispersions a factor of 2-3 smaller those of the elliptical galaxies hosting QSOs at z > 2. A simple consequence of our hypotheses is that the redshift evolution of the QSO emissivity and of the star formation history in spheroids should be roughly parallel. We find this result to be broadly consistent with our knowledge of the evolution of both the global star formation rate, and of the evolution of the QSO emissivity, but we identify interesting discrepancies at both low and high redshifts, to which we offer tentative solutions. Finally, our hypotheses allow us to present a robust method to derive the duty cycle of QSO activity, based on the observed QSO luminosity function, and on the present-day relation between the masses of supermassive BHs and those of their spheroidal host stellar systems. The duty cycle is found to be substantially less than unity, with characteristic values in the range (3-6)x10^(-3), and we compute that the average bolometric radiative efficiency is epsilon=0.07. Finally, we find that the growth in mass of individual black holes at high redshift (z>2) can be dominated by mergers, and is therefore not necessarily limited by accretion.

We study the contribution of projection effects to the intrinsic thickness of the Fundamental Plane (FP) of elliptical galaxies. The Monte-Carlo mapping technique between model properties and observed quantities, introduced by Bertin, Ciotti, and Del Principe (2002), is extended to oblate, two-integrals galaxy models, with non-homologous density profiles, adjustable flattening, variable amount of ordered rotational support, and for which all the relevant projected dynamical quantities can be expressed in fully analytical way. In agreement with previous works, it is found that projection effects move models not exactly parallel to the edge-on FP, by an amount that can be as large as the observed FP thickness. The statistical contribution of projection effects to the FP thickness is however marginal, and the estimated physical FP rms thickness is ~90% of the observed one (when corrected for measurement errors).

We explore the effects of dissipationless merging on the Fundamental Plane of elliptical galaxies by using a N-body code based on a new, high performance numerical scheme (Dehnen 2002). We investigate the extreme cases of galaxy growth by equal mass merging and accretion of small stellar systems; in a subset of simulations we also consider the presence of dark matter halos around the merging galaxies. Curiously, we found that the Fundamental Plane is preserved by major merging, while in the accretion scenario its edge--on thickness is only marginally reproduced, with substantial thickening in the case of merging with low angular momentum. We also found that both the Faber-Jackson and Kormendy relations are not reproduced by the simulations. Finally, we discuss the implications of our results for the origin of the Mbh-sigma and Magorrian relations. We found that dissipationless merging is unable to reproduce the Mbh-sigma relation, if Mbhs add linearly (while the Magorrian relation is nicely reproduced); on the contrary a black hole merging with substantial emission of gravitational waves reproduces the Mbh-sigma relation but fails at reproducing the Magorrian relation. [Abridged]

We describe a particular aspect of the effects of the parent cluster tidal field (CTF) on stellar orbits inside cluster Elliptical galaxies. In particular we discuss, with the aid of a simple numerical model, the possibility that collisionless stellar evaporation from elliptical galaxies is an effective mechanism for the production of the recently discovered intracluster stellar populations. A preliminary investigation, based on very idealized galaxy density profiles (Ferrers density distributions), showed that over an Hubble time, the amount of stars lost by a representative galaxy may sum up to the 10% of the initial galaxy mass, a fraction in interesting agreement with observational data. The effectiveness of this mechanism is due to the fact that the galaxy oscillation periods near equilibrium configurations in the CTF are comparable to stellar orbital times in the external galaxy regions. Here we extend our previous study to more realistic galaxy density profiles, in particular by adopting a triaxial Hernquist model.

With the aid of simple numerical models, we discuss a particular aspect of the interaction between stellar orbital periods inside elliptical galaxies (Es) and the parent cluster tidal field (CTF), i.e., the possibility that collisionless stellar evaporation from Es is an effective mechanism for the production of the recently discovered intracluster stellar populations (ISP). These very preliminary investigations, based on idealized galaxy density profiles (such as Ferrers density distributions) show that, over an Hubble time, the amount of stars lost by a representative galaxy may sum up to the 10% of the initial galaxy mass, a fraction in interesting agreement with observational data. The effectiveness of this mechanism is due to the fact that the galaxy oscillation periods near its equilibrium configurations in the CTF are of the same order of stellar orbital times in the external galaxy regions.

Pure N-body high-resolution re-simulations of 13 massive dark matter halos in a Lambda-CDM cosmology are presented. The resulting sample is used to investigate the physical origin of the Fundamental Plane, ``Faber-Jackson'', and ``Kormendy'' relations observed for nearby galaxy clusters. In particular, we focus on the role of dissipationless hierarchical merging in establishing or modifying these relations. Contrarily to what found in the case of the Faber-Jackson and Kormendy relations for galaxies (see Londrillo, Nipoti and Ciotti, this conference), dissipationless merging on cluster scales produces scaling relations remarkably similar to the observed ones. This suggests that gas dissipation plays a minor role in the formation of galaxy clusters, and the hierarchical merger scenario is not at odd with the observed regularity of these systems.

A new, momentum preserving fast Poisson solver for N-body systems sharing effective O(N) computational complexity, has been recently developed by Dehnen (2000, 2002). We have implemented the proposed algorithms in a Fortran-90 code, and parallelized it by a domain decomposition using the MPI routines. The code has been applied to intensive numerical investigations of galaxy mergers, in particular focusing on the possible origin of some of the observed scaling relations of elliptical galaxies. We found that the Fundamental Plane is preserved by an equal mass merging hierarchy, while it is not in a scenario where galaxies grow by accretion of smaller stellar systems. In addition, both the Faber-Jackson and Kormendy relations are not reproduced by our simulations.

We start by studying a small set of objects characterized by photometric profiles that have been pointed out to deviate significantly from the standard R^{1/4} law. For these objects we confirm that a generic R^{1/n} law, with n a free parameter, can provide superior fits (the best-fit value of n can be lower than 2.5 or higher than 10), better than those that can be obtained by a pure R^{1/4} law, by an R^{1/4}+exponential model, and by other dynamically justified self--consistent models. Therefore, strictly speaking, elliptical galaxies should not be considered homologous dynamical systems. Still, a case for "weak homology", useful for the interpretation of the Fundamental Plane of elliptical galaxies, could be made if the best-fit parameter n, as often reported, correlates with galaxy luminosity L, provided the underlying dynamical structure also follows a systematic trend with luminosity. We demonstrate that this statement may be true even in the presence of significant scatter in the correlation n(L). Preliminary indications provided by a set of "data points" associated with a sample of 14 galaxies suggest that neither the strict homology nor the constant stellar mass--to--light solution are a satisfactory explanation of the observed Fundamental Plane (abridged).

The existence of the Fundamental Plane (FP) imposes strong constraints on the structure and dynamics of elliptical galaxies, and thus contains important information on the processes of their formation and evolution. Here we focus on the relations between the FP thinness and tilt and the amount of radial orbital anisotropy. By using N-body simulations of galaxy models characterized by observationally motivated density profiles, and also allowing for the presence of live, massive dark matter halos, we explore the impact of radial orbital anisotropy and instability on the FP properties. The numerical results confirm a previous semi--analytical finding: the requirement of stability matches almost exactly the thinness of the FP. In other words, galaxy models that are radially anisotropic enough to be found outside the observed FP (with their isotropic parent models lying on the FP) are unstable, and their end--products fall back on the FP itself. We also find that a systematic increase of radial orbit anisotropy with galaxy luminosity cannot explain by itself the whole tilt of the FP, becoming the galaxy models unstable at moderately high luminosities: at variance with the previous case their end--products are found well outside the FP itself (abridged).

We present preliminary results of numerical simulations of dissipationless merging of stellar systems, aimed at exploring the consequences of merging between gas free, spheroidal systems. In particular, we study the dynamical and structural characteristics of hierarchical merging between equal mass stellar systems, and we compare the properties of the end-products with the most important structural and dynamical scaling relations obeyed by spheroids. In the explored hierarchy of four successive mergings we find that the FP tilt is marginally conserved, but both the Faber-Jackson and Kormendy relations are not conserved.

We present a robust method to derive the duty cycle of QSO activity based on the empirical QSO luminosity function and on the present-day linear relation between the masses of supermassive black holes and those of their spheroidal host stellar systems. It is found that the duty cycle is substantially less than unity, with characteristic values in the range $3-6\times 10^{-3}$. Finally, we tested the expectation that the QSO luminosity evolution and the star formation history should be roughly parallel, as a consequence of the above--mentioned relation between BH and galaxy masses.

Recovering the true average abundance of the intracluster medium (ICM) is crucial for estimates of its global metal content, which in turn is linked to its past evolution and to the star formation history of the stellar component of the cluster. We analyze here how the presence of abundance gradients affects commonly adopted estimates of the average abundance, assuming various plausible ICM density and temperature profiles. We find that, adopting the observed abundance gradients, the true average mass weighted abundance is less than (although not largely deviating from) the commonly used emission weighted abundance.

Most elliptical galaxies contain central black holes, whose masses scale proportionally to the observed central velocity dispersions of the host galaxies according to the so-called Mbh-sigma relation. Here we discuss some consequences that can be derived by combining the Mbh-sigma relation with the scaling relation describing the Fundamental Plane of elliptical galaxies. In particular, the possibility of substantial dissipationless merging in the formation and evolution of elliptical galaxies is discussed. Enforcing the merger end-products to satisfy the two scaling relations mentioned above, a major role of dissipation in galaxy formation is strongly suggested by our analysis. Moreover, we show that existing observational data may shed some light on the complex process of black hole merging.

Most elliptical galaxies contain central black holes (BHs), and most also contain significant amounts of hot gas capable of accreting on to the central BH due to cooling times short compared to the Hubble time. Why therefore do we not see AGNs at the center of most elliptical galaxies rather than in only (at most) a few per cent of them? We propose here the simple idea that feedback from accretion events heats the ambient gas retarding subsequent infall. In this context, we present a class of 1D hydrodynamical evolutionary sequences for the gas flows in elliptical galaxies with a massive central BH. The resulting evolution is characterized by strong oscillations, in which very fast and energetic bursts of the BH are followed by longer periods in which the X-ray galaxy emission comes from the coronal gas. We also allow for departures from spherical symmetry by examining scenarios in which the central engine is either an ADAF or a more conventional accretion disk that is optically thick except for a polar region. In all cases the duty cycle (fraction of the time that the system will be seen as an AGN) is quite small and in the range 10^{-4} - 10^{-3}. Thus, for any reasonable value of the efficiency, the presence of a massive BH at the center of a galaxy seems to be incompatible with the presence of a long-lived cooling flow.

In this paper we describe some analytical properties of the R^{1/m} law proposed by Sersic (1968) to categorize the photometric profiles of elliptical galaxies. In particular, we present the full asymptotic expansion for the dimensionless scale factor b(m) that is introduced when referring the profile to the standard effective radius. Surprisingly, our asymptotic analysis turns out to be useful even for values of m as low as unity, thus providing a unified analytical tool for observational and theoretical investigations based on the R^{1/m} law for the entire range of interesting photometric profiles, from spiral to elliptical galaxies.

During a Hubble time, cluster galaxies may undergo several mutual encounters close enough to gravitationally perturb their hot, X-ray emitting gas flows. We ran several 2D, time dependent hydrodynamical models to investigate the effects of such perturbations on the gas flow inside elliptical galaxies. In particular, we studied in detail the modifications occurring in the scenario proposed by D'Ercole et al. (1989), in which the galactic interstellar medium produced by the aging galactic stellar population, is heated by SNIa at a decreasing rate. We find that, although the tidal interaction in our models lasts less than 1 Gyr, its effect extends over several Gyrs. The tidally induced turbulent flows create dense filaments which cool quickly and accrete onto the galactic center, producing large spikes in the global Lx. Once this mechanism starts, it is fed by gravity and amplified by SNIa. In cooling flow models without supernovae the amplitude of the Lx fluctuations due to the tidal interaction is substantially reduced. We conclude that, if SNIa significantly contribute to the energetics of the gas flows in ellipticals, then the observed spread in the Lx-Lb diagram may be caused, at least in part, by this mechanism. On the contrary, tidal interactions cannot be responsible for the observed spread if the pure cooling flow scenario applies (abridged).

In this paper we describe some analytical properties of the R^{1/m} law proposed by Sersic to categorize the photometric profiles of elliptical galaxies. In particular, we present the full asymptotic expansion for the dimensionless scale factor b(m) that is introduced when referring the profile to the standard effective radius. Surprisingly, our asymptotic analysis turns out to be useful even for values of m as low as unity, thus providing a unified analytical tool for observational and theoretical investigations based on the R^{1/m} law for the entire range of interesting photometric profiles, from spiral to elliptical galaxies.

There are observational and theoretical indications that both the visible (stars) and the dark matter density distributions in elliptical galaxies increase significantly up to the galactic center. I present here some analytical results obtained with the aid of self-consistent, spherically symmetric two component galaxy models. These results suggest the possibility that this similar behavior could be a direct consequence of the structural and dynamical constraints imposed by the request of positivity of the phase-space distribution function of each density component.

In the context of the study of the properties of the mutual mass distribution of the bright and dark matter in elliptical galaxies, present a family of two-component, spherical, self-consistent galaxy models, where one density distribution follows a gamma_1 profile, and the other a gamma_2 profile [(gamma_1,gamma_2) models], with different total masses and ``core'' radii. A variable amount of Osipkov-Merritt (radial) orbital anisotropy is allowed in both components. For these models, I derive analytically the necessary and sufficient conditions that the model parameters must satisfy in order to correspond to a physical system. Moreover, the possibility of adding a black hole at the center of radially anisotropic gamma models is discussed, determining analytically a lower limit of the anisotropy radius as a function of gamma. The analytical phase-space distribution function for (1,0) models is presented, together with the solution of the Jeans equations and the quantities entering the scalar virial theorem. It is proved that a globally isotropic gamma=1 component is consistent for any mass and core radius of the superimposed gamma=0 model; on the contrary, only a maximum value of the core radius is allowed for the gamma=0 model when a gamma=1 density distribution is added. The combined effects of mass concentration and orbital anisotropy are investigated, and an interesting behavior of the distribution function of the anisotropic gamma=0 component is found: there exists a region in the parameter space where a sufficient amount of anisotropy results in a consistent model, while the structurally identical but isotropic model would be inconsistent.

In this paper I present three new results of astronomical interest concerning the theory of Abel inversion. 1) I show that in the case of a spatial emissivity that is constant on toroidal surfaces and projected along the symmetry axis perpendicular to the torus' equatorial plane, it is possible to invert the projection integral. From the surface (i.e. projected) brightness profile one then formally recovers the original spatial distribution as a function of the toroidal radius. 2) By applying the above-described inversion formula, I show that if the projected profile is described by a truncated off-center gaussian, the functional form of the related spatial emissivity is very simple and - most important - nowhere negative for any value of the gaussian parameters, a property which is not guaranteed - in general - by Abel inversion. 3) Finally, I show how a generic multimodal centrally symmetric brightness distribution can be deprojected using a sum of truncated off-center gaussians, recovering the spatial emissivity as a sum of nowhere negative toroidal distributions.

The properties of the analitycal phase--space distribution function of two-component spherical self-consistent galaxy models, where one density distribution follows the Hernquist profile, and the other a gamma=0 model, with different total masses and core radii (H0 models), presented in Ciotti (1998), are here summarized. A variable amount of radial Osipkov-Merritt orbital anisotropy is allowed in both components. The necessary and sufficient conditions that the model parameters must satisfy in order to correspond to a model where each one of the two distinct components has a positive DF (the so-called model consistency) are analytically derived, together with some results on the more general problem of the consistency of two-component gamma1+gamma2 models. The possibility to add in a consistent way a black hole at the center of radially anisotropic gamma models is also discussed. In particular, it is proved that a globally isotropic Hernquist component is consistent for any mass and core radius of the superimposed gamma=0 halo; on the contrary, only a maximum value of the core radius is allowed to the gamma=0 component when a Hernquist halo is added. The combined effect of halo concentration and orbital anisotropy is successively investigated.

We present the results of a new set of hydrodynamical simulations of hot gas flows in model elliptical galaxies with the following characteristics: the spatial luminosity distribution approaches a power law form at small radii, in accordance with the results of recent ground based observations, and with the Hubble Space Telescope; the dark matter has a peaked profile too, as indicated by high resolution numerical simulations of dissipationless collapse; the dark to luminous mass ratio spans a large range of values, including low values found by optical studies confined to within two effective radii; finally, the type Ia supernova rate is that given by the latest estimates of optical surveys, or zero, as suggested by the iron abundances recently measured in the hot gas. We find that the resulting gas flows are strongly decoupled: an inflow develops in the central region of the galaxies, while the external parts are still degassing, i.e., the flows are mostly partial winds. This behavior can be explained in terms of the local energy balance of the hot gas. A large spread in the X-ray luminosity Lx at fixed optical luminosity Lb can be produced as in previous simulations that used King models plus massive quasi-isothermal dark halos, and higher supernova rates; the key factor causing large Lx variations is now the size of the central inflow region. The highest Lx observed correspond to global inflows. Finally, non negligible amounts of cold gas can be produced by the partial winds; this could be an explanation for the possible discovery of cold matter at the center of elliptical galaxies, an alternative to the presence of a steady state cooling flow.

We investigate the motion, near the equilibrium configurations, of an initially spinless rigid body subject to an external tidal field. Two cases are considered: when the center of mass of the body is at rest at the equilibrium point of the field generated by a generic mass distribution, and when it is placed on a circular orbit subject to a spherically symmetric potential. A complete analysis of the equilibrium configurations is carried out for both cases. First, we derive the conditions for the stable equilibria, and then we analyze the frequencies of oscillations around the equilibrium positions. In view of these results, we consider the problem of alignment of galaxies in clusters. After estimating the period of the oscillations induced on the galaxies by the tidal field of the cluster, we discuss the possible effect of resonances between stellar orbits inside the galaxy and the oscillations of the galaxy as a whole; this may be a mechanism responsible for producing an intracluster stellar population.

A recent analysis of the "Einstein" sample of early-type galaxies has revealed that at any fixed optical luminosity Lb S0 galaxies have lower mean X-ray luminosity Lx per unit Lb than ellipticals. Following a previous analytical investigation of this problem (Ciotti & Pellegrini 1996), we have performed 2D numerical simulations of the gas flows inside S0 galaxies in order to ascertain the effectiveness of rotation and/or galaxy flattening in reducing the Lx/Lb ratio. The flow in models without SNIa heating is considerably ordered, and essentially all the gas lost by the stars is cooled and accumulated in the galaxy center. If rotation is present, the cold material settles in a disk on the galactic equatorial plane. Models with a time decreasing SNIa heating host gas flows that can be much more complex. After an initial wind phase, gas flows in energetically strongly bound galaxies tend to reverse to inflows. This occurs in the polar regions, while the disk is still in the outflow phase. In this phase of strong decoupling, cold filaments are created at the interface between inflowing and outflowing gas. Models with more realistic values of the dynamical quantities are preferentially found in the wind phase with respect to their spherical counterparts of equal Lb. The resulting Lx of this class of models is lower than in spherical models with the same Lb and SNIa heating. At variance with cooling flow models, rotation is shown to have only a marginal effect in this reduction, while the flattening is one of the driving parameters for such underluminosity, in accordance with the analytical investigation.

We present a new class of solutions for the gas flows in elliptical galaxies containing massive central black holes (BH). Modified King model galaxies are assumed. Two source terms operate: mass loss from evolving stars, and a secularly declining heating by SNIa. Relevant atomic physical processes are modeled in detail. Like the previous models investigated by Ciotti et al. (1991), these new models first evolve through three consecutive evolutionary stages: wind, outflow, and inflow. At this point the presence of the BH alters dramatically the subsequent evolution, because the energy emitted by the BH can heat the surrounding gas to above virial temperatures, causing the formation of a hot expanding central bubble. Short and strong nuclear bursts of radiation are followed by longer periods during which the X-ray galaxy emission comes from the coronal gas (Lx). The range and approximate distribution spanned by Lx are found to be in accordance with observations of X-ray early type galaxies. Moreover, although high accretion rates occur during bursting phases when the central BH has a luminosity characteristic of QSOs, the total mass accreted is very small when compared to that predicted by stationary cooling-flow solutions and computed masses are in accord with putative BH nuclear masses. In the bursting phases Lx is low and the surface brightness profile is very low compared to pre-burst or to cooling flow models. We propose that these new models, while solving some long-standing problems of the cooling flow scenario, can provide a unified description of QSO-like objects and X-ray emitting elliptical galaxies, these being the same objects observed at two different evolutionary phases.

In the first contribution the basic problems posed by the existence of the Fundamental Plane (FP), and its relations with the Virial Theorem are discussed. Various possibilities are presented that can produce the observed uniform departure from homology and/or from a constant stellar mass-to-light ratio. The role of orbital anisotropy and its relation with the FP thickness are also discussed. In the second contribution the observational evidence for the existence of a FP for galaxy clusters, analogous to that of ellipticals, is presented and discussed. In the third contribution the dynamical properties of $R^{1/m}$ spherical galaxies with variable orbital anisotropy are discussed. The stability of the models against radial--orbit instability is studied, and the consequences on the FP tilt and thickness are discussed. An innovative tool for the construction of Evolutionary Synthesis models of Stellar Populations is described in the fourth contribution. As an illustrative example, the evolution of the stellar mass-to-light ratios has been computed and the implications on the FP of elliptical galaxies are briefly discussed. Finally, in the fifth contribution, with the aid of the first large sample of X-ray emitting early-type galaxies with known kinematics, the effect of rotation on the X-ray emission, particularly with regard to the X-ray underluminosity of flat systems is discussed. The role of structural parameters as the ellipticity, the isophotal shape parameter $a_4/a$, and the inner slope of the surface brightness profile recently measured by HST are also investigated.

Long slit spectroscopic data for 7 early-type galaxies with X-ray emission have been analyzed to derive velocity dispersion and radial velocity profiles. Major axis rotation curves out to $R\sim R_e$ are presented. Adding these new data to those available in the literature, we have built a sample of early-type galaxies with detected X-ray emission and known kinematics (central velocity dispersion $\sigma_0$ and maximum rotational velocity $v_r$). Using this sample we investigate from an observational point of view the role of rotation and flattening on the X-ray emission, particularly with regard to the X-ray underluminosity of flat systems. The trend between the X-ray to optical ratio $L_X/L_B$, a measure of the hot gas content of the galaxies, and $v_r/\sigma_0$ is L-shaped, with the X-ray brightest objects confined at $v_r/\sigma_0 <0.4$, both for Es and S0s. Neither for low or intermediate, nor for high $L_X/L_B$, there is any clear correlation between X-ray emission and rotational properties. The trend between $L_X/L_B$ and the ellipticity $\epsilon$ is also L-shaped: there are no high $L_X/L_B$ objects with high $\epsilon$. The existence of a relation between $v_r/\sigma_0$ and $\epsilon$ makes it difficult to assess whether rotation or flattening is at the basis of the trends found. Our observational findings are finally discussed in connection with the effects that rotation and flattening are predicted to have in the cooling flow and in the wind/outflow/inflow scenarios.

We study the dynamical properties of spherical galaxies with surface luminosity profile described by the R^{1/m}-law, in which a variable degree of orbital anisotropy is allowed. The parameter m for the present models covers the range [1,10]. For these models we study the self-consistently generated phase-space distribution function (DF), and we derive - as a function of m - the minimum value of the anisotropy radius for the model consistency (i.e., in order to have a nowhere negative DF). Then we study the region in the parameter space where the R^{1/m} models are likely to be stable against radial-orbit instability, and we compare its size with that of the larger region corresponding to the consistency requirement. For stable anisotropic models the spatial and projected velocity dispersion profiles are obtained solving the Jeans equation, and compared to those of the globally isotropic case, already discussed in Paper I. The relevance of the results in connection with the Fundamental Plane (FP) of elliptical galaxies is pointed out: the effect on the projected velocity dispersion due to the maximum orbital anisotropy allowed by the stability requirement is well within the FP thickness, and so no fine-tuning for anisotropy is required. Finally, the Velocity Profiles are constructed as function of the projected radius and for various degrees of anisotropy, and their deviations from a gaussian discussed.

The analytical phase-space distribution function (DF) of spherical self--consistent galaxy (or cluster) models, embedded in a dark matter halo, where both density distributions follow the Hernquist profile, with different total masses and core radii (hereafter called HH models), is presented. The concentration and the amount of the stellar and dark matter distributions are described by four parameters: the mass and core radius of the {\it reference} component, and two dimensionless parameters describing the mass and core radius of the {\it halo} component. A variable amount of orbital anisotropy is allowed in both components, following the widely used parameterization of Osipkov-Merritt. An important case is obtained for a null core radius of the halo, corresponding to the presence of a central black hole (BH). It is proved that globally isotropic HH models are consistent for any mass ratio and core radii ratio, even in the case of a central BH. In this last case the analytical expression for a lower limit of the anisotropy radius of the host system as a function of the BH mass is given. In the particular case of global isotropy the stability of HH models is proved, and the stability of radially anisotropic HH models is briefly discussed. The expression derived for the DF is useful for understanding the relations between anisotropy, density shape and external potential well in a consistent stellar system, and to produce initial conditions for N-body simulations of two-component galaxies or galaxy clusters.

In this paper we explore several structural and dynamical effects on the projected velocity dispersion as possible causes of the fundamental plane (FP) tilt of elliptical galaxies. Specifically, we determine the size of the systematic trend along the FP in the orbital radial anisotropy, in the dark matter (DM) content and distribution relative to the bright matter, and in the shape of the light profile that would be needed to produce the tilt, under the assumption of a constant stellar mass to light ratio. Spherical, non rotating, two--components models are constructed, where the light profiles resemble the $R^{1/4}$ law. For the investigated models anisotropy cannot play a major role in causing the tilt, while a systematic increase in the DM content and/or concentration may formally produce it. Also a suitable variation of the shape of the light profile can produce the desired effect, and there may be some observational hints supporting this possibility. However, fine tuning is always required in order to produce the tilt, while preserving the {\it tightness} of the galaxies distribution about the FP.

The global energy budget of the hot gas flows in early-type galaxies is used to explain the recent observational results that S0 galaxies have lower mean X-ray luminosity $L_X$ per unit optical luminosity $L_B$ than do ellipticals. This could be explained with a higher heat input or with a lower gravitational energy (all per unit gas mass), at fixed $L_B$. The effect of the flattening of the mass distribution is investigated with galaxy models described by the Miyamoto-Nagai potential--density pair, to which a dark matter halo of various shapes is added. For these models the analytical expressions of the gravitational energy, and the stellar kinetic energy associated with various relative amounts of random motions and rotational streaming, are derived. It is found that rotation cannot produce a change in the flow phase of the hot gas, independently of the galaxy shape and the presence of dark matter. The effect of flattening instead can be substantial in reducing the binding energy of the hot gas. Thus S0s and possibly non spherical Es are less able to retain a significant halo of hot gas than rounder Es of the same $\lb$.

We address the problem of deriving, from the observed projected metallicity gradients, the intrinsic metallicity distribution of elliptical galaxies as a function of their integrals of motion. The method is illustrated by an application to anisotropic spherical Hernquist models. We also compare the derived metallicity distribution with those expected from two very simple models of galaxy formation and find that the more dissipative scheme agrees better with the typical metallicity distribution of ellipticals.

In this paper the concept of tomography of a collisionless stellar system of general shape is introduced, and a generalization of the Projected Virial Theorem is obtained. Applying the tomographic procedure we then derive a new family of virial equations which coincides with the already known ones for spherically symmetric systems. This result is obtained without any use of explicit expressions for the line-of-sight velocity dispersion, or spherical coordinate system.

We investigate two possible effects of the tidal field induced by a spherical cluster on its elliptical galaxy members: the modification of the ellipticity of a spherical galaxy and the isophotal alignment in the cluster radial direction of a misaligned prolate galaxy. Numerical N-body simulations have been performed for radial and circular galactic orbits. The properties of the stars' zero--velocity surfaces in the perturbed galaxies are explored briefly, and the adiabaticity of the galaxy to the external field is discussed. For a choice of parameters characteristic of rich clusters we find that the induced ellipticity on a spherical galaxy is below or close to the detectability level. But we find that the tidal torque can result in significant isophotal alignment of the galaxies' major axis with the cluster radial direction if the galaxy is outside the cluster core radius. The time required for the alignment is very short compared with the Hubble time. A significant increase in the ellipticity of the outer isophotes of the prolate model is also found, but with no observable isophotal twisting. Our main prediction is an alignment segregation of the elliptical galaxy population according to whether their orbits lie mostly outside or inside the cluster core radius. These results also suggest that galactic alignment in rich clusters is not incompatible with a bottom-up galaxy formation scenario.