Phase mixing in MOND

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.Comment: 10 pages, 3 figures. To appear in the Proceedings of the International Workshop "Collective Phenomena in Macroscopic Systems", G. Bertin, R. Pozzoli, M. Rome, and K.R. Sreenivasan, eds., World Scientific, Singapor

N-body simulations in modified Newtonian dynamics

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.Comment: 4 pages, no figures. To appear in EAS Publication Series (Proceedings of Symposium 7 of the JENAM 2008, Vienna

Decoupled and inhomogeneous gas flows in S0 galaxies

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.Comment: 32 pages LaTex file, plus 5 .ps figures and macro aasms4.sty -- Accepted on Ap

Cooling flows and quasars: different aspects of the same phenomenon? I. Concepts

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.Comment: 10 pages, ApJ LaTeX, plus 5 .eps figures and TeX-macro aasms4.sty - revised version - in press on ApJ Letter

Radiative feedback from massive black holes in elliptical galaxies. AGN flaring and central starburst fueled by recycled gas

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]Comment: ApJ resubmitted. 48 pages, 14 figures (some of them new, bitmapped, low resolution). New references added, typos correcte

Modelling elliptical galaxies: phase-space constraints on two-component (gamma1,gamma2) models

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.Comment: 29 pages, LaTex, plus 5 .eps figures and macro aaspp4.sty - accepted by ApJ, main journa

Formation of central massive objects via tidal compression

For a density that is not too sharply peaked towards the center, the local tidal field becomes compressive in all three directions. Available gas can then collapse and form a cluster of stars in the center, including or even being dominated by a central black hole. We show that for a wide range of (deprojected) Sersic profiles in a spherical potential, the tidal forces are compressive within a region which encloses most of the corresponding light of observed nuclear clusters in both late-type and early-type galaxies. In such models, tidal forces become disruptive nearly everywhere for relatively large Sersic indices n >= 3.5. We also show that the mass of a central massive object (CMO) required to remove all radial compressive tidal forces scales linearly with the mass of the host galaxy. If CMOs formed in (progenitor) galaxies with n ~ 1, we predict a mass fraction of ~ 0.1-0.5%, consistent with observations of nuclear clusters and super-massive black holes. While we find that tidal compression possibly drives the formation of CMOs in galaxies, beyond the central regions and on larger scales in clusters disruptive tidal forces might contribute to prevent gas from cooling.Comment: 19 pages, 4 figures. Accepted for publication in ApJ. High resolution version available at http://www-obs.univ-lyon1.fr/labo/perso/eric.emsellem/preprint

Galactic fountains and gas accretion

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.Comment: 2 pages, 1 figure, to appear in "Hunting for the Dark: The Hidden Side of Galaxy Formation", Malta, 19-23 Oct. 2009, eds. V.P. Debattista & C.C. Popescu, AIP Conf. Se

Galactic cannibalism in the galaxy cluster C0337-2522 at z=0.59

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].Comment: 30 pages, 7 figures, MNRAS (accepted

Fountain-driven gas accretion by the Milky Way

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.Comment: 2 pages, 1 figure. To appear in the proceedings of the conference "Assembling the Puzzle of the Milky Way", Le Grand-Bornand 17-22 April 2011, European Physical Journal, editors C. Reyl\'e, A. Robin and M. Schulthei