BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Luca Ciotti

Luca Ciotti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-phastro-ph.HEmath-phmath.MPastro-ph.EPastro-ph.SRmath.CAnlin.CDphysics.plasm-ph

Catalog footprint

What is connected

37works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2022arXiv

A Parameter Space Exploration of High Resolution Numerically Evolved Early Type Galaxies Including AGN Feedback and Accurate Dynamical Treatment of Stellar Orbits

An extensive exploration of the model parameter space of axisymmetric Early-Type Galaxies (ETGs) hosting a central supermassive Black Hole (SMBH) is conducted by means of high resolution hydrodynamical simulations performed with our code MACER. Global properties such as 1) total SMBH accreted mass, 2) final X-ray luminosity and temperature of the X-ray emitting halos, 3) total amount of new stars formed from the cooling gas, 4) total ejected mass in form of supernovae and AGN feedback induced galactic winds, are obtained as a function of galaxy structure and internal dynamics. In addition to the galactic dark matter halo, the model galaxies are also embedded in a group/cluster dark matter halo; finally cosmological accretion is also included, with amount and time dependence derived from cosmological simulations. Angular momentum conservation leads to the formation of cold HI disks; these disks further evolve under the action of star formation induced by disk instabilities, of the associated mass discharge onto the central SMBH, and of the consequent AGN feedback. At the end of the simulations, the hot (metal enriched) gas mass is roughly $10\%$ the mass in the old stars, with twice as much having been ejected into the intergalactic medium. The cold gas disks are a $\approx$ kpc in size, and the metal rich new stars are in $0.1$ kpc disks. The masses of cold gas and new stars are roughly $0.1\%$ the mass of the old stars. Overall, the final systems appear to reproduce quite successfully the main global properties of real ETGs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the polytropic Bondi accretion in two-component galaxy models with a central massive BH

In many investigations involving accretion on a central point mass, ranging from observational studies to cosmological simulations, including semi-analytical modelling, the classical Bondi accretion theory is the standard tool widely adopted. Previous works generalised the theory to include the effects of the gravitational field of the galaxy hosting a central black hole, and of electron scattering in the optically thin limit. Here we apply this extended Bondi problem, in the general polytropic case, to a class of new two-component galaxy models recently presented. In these models, a Jaffe stellar density profile is embedded in a dark matter halo such that the total density distribution follows a $r^{-3}$ profile at large radii; the stellar dynamical quantities can be expressed in a fully analytical way. The hydrodynamical properties of the flow are set by imposing that the gas temperature at infinity is proportional to the virial temperature of the stellar component. The isothermal and adiabatic (monoatomic) cases can be solved analytically, in the other cases we explore the accretion solution numerically. As non-adiabatic accretion inevitably leads to an exchange of heat with the ambient, we also discuss some important thermodynamical properties of the polytropic Bondi accretion, and provide the expressions needed to compute the amount of heat exchanged with the environment, as a function of radius. The results can be useful for the subgrid treatment of accretion in numerical simulations, as well as for the interpretation of observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A family of Exponential Integrals suggested by Stellar Dynamics

While investigating the generalization of the Chandrasekhar (1943) dynamical friction to the case of field stars with a power-law mass spectrum and equipartition Maxwell-Boltzmann velocity distribution, a pair of 2-dimensional integrals involving the Error function occurred, with closed form solution in terms of Exponential Integrals (Ciotti 2010). Here we show that both the integrals are very special cases of the family of (real) functions $$ I(λ,μ,ν; z) :=\int_0^zx^λ\,\Enu(x^μ)\,dx= {γ\left({1+λ\overμ},z^μ\right) + z^{1+λ}\Enu(z^μ)\over 1+λ+ μ(ν-1)}, \quad μ>0,\quad z\geq 0, \eqno (1) $$ where $\Enu$ is the Exponential Integral, $γ$ is the incomplete Euler gamma function, and for existence $λ>\max \left\{-1,-1- μ(ν-1)\right\}$. Only in one of the consulted tables a related integral appears, that with some work can be reduced to eq.~(1), while computer algebra systems seem to be able to evaluate the integral in closed (and more complicated) form only provided numerical values for some of the parameters are assigned. Here we show how eq.~(1) can in fact be established by elementary methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Infra-Red Emission from Cold Gas Dusty Disks in Massive Ellipticals

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Brownian motion of supermassive black holes in galaxy cores

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Dynamical models of spheroidal multi-component stellar 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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the Role of Hot Feedback Mode in Active Galactic Nuclei Feedback in an Elliptical 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The mass discrepancy acceleration relation in early-type galaxies: extended mass profiles and the phantom menace to MOND

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fidelity and Reversibility in the Restricted Three Body Problem

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Miyamoto-Nagai discs embedded in the Binney logarithmic potential: analytical solution of the two-integrals Jeans equations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Radially anisotropic systems with $r^{-α}$ forces: equilibrium states

We continue the study of collisionless systems governed by additive $r^{-α}$ interparticle forces by focusing on the influence of the force exponent $α$ 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^{-α}$ forces and $1\leqα<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 $α$ 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 $α$, 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 $α$. In particular, we find that isotropic systems are consistent in the explored range of $α$. By means of direct $N$-body simulations we finally verify that the isotropic systems are also stable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

X-ray haloes and star formation in early-type galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Action-based distribution functions for spheroidal galaxy components

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

AGN feedback in an isolated elliptical galaxy: the effect of strong radiative feedback in the kinetic mode

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamical friction in multi-component evolving globular clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Grain Physics and IR Dust Emission in AGN Environments

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The effects of galaxy shape and rotation on the X-ray haloes of early-type galaxies - II. Numerical simulations

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$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The effects of stellar dynamics on the X-ray emission of flat early-type galaxies

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].

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The imprint of dark matter haloes on the size and velocity dispersion evolution of early-type galaxies

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 ($σ_0$) and half-mass radius ($r_{\rm h}$) as functions of halo mass $M$ and redshift $z$, in a cosmological $Λ$-CDM $N$-body simulation. In the range $0\lesssim z\lesssim 2.5$, we find $σ_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 $σ_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 $σ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The stellar initial mass function of early type galaxies from low to high stellar velocity dispersion: homogeneous analysis of ATLAS$^{\rm 3D}$ and Sloan Lens ACS galaxies

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 ($σ$) within the SDSS fibers. Comparing the dynamically-derived stellar mass-to-light ratios $(M_*/L)_{\rm dyn}$, obtained for an assumed halo slope $ρ_{\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 $σ$ galaxies is consistent on average with a Salpeter slope. Our study allows for a fully consistent study of the trend between IMF and $σ$ for both the SLACS and \ATLAS samples, which explore quite different $σ$ ranges. The two samples are highly complementary, the first being essentially $σ$ 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α=(0.38\pm0.04)\times\log(σ_{\rm e}/200\,\mathrm{km~s}^{-1})+(-0.06\pm0.01)$, where $α=(M_*/L)_{\rm dyn}/(M_*/L)_{\rm pop}$ and $σ_{\rm e}$ is the luminosity averaged $σ$ 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 $σ_{\rm e}\approx90-270\,\mathrm{km~s}^{-1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Disk dynamics and the X-ray emission of S0 and flat early-type galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Massive black holes in central cluster galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relaxation of N-body systems with additive r^-α interparticle forces

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^{-α}$ forces. We find that most of the results holding in Newtonian gravity are also valid for $α\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érsic law, the global density slope-anisotropy inequality is obeyed, the differential energy distribution is an exponential over a large range of energies (for $α\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 $α$, and becomes infinite for $α=-1$ (the harmonic oscillator). This is in agreement with the results of deep-MOND collapses (qualitatively corresponding to $α=1$) and it is due to the fact the force becomes more and more similar to the $α=-1$ case, where as well known no relaxation can happen and the system oscillates forever.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The effects of flattening and rotation on the temperature of the X-ray halos of elliptical galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Effects of Galaxy Shape and Rotation on the X-ray Haloes of Early-Type Galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Environmental effects on the growth of super massive black holes and AGN feedback

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)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radial-orbit instability 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamical Friction from field particles with a mass spectrum

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Feedback from central black holes in elliptical galaxies. II: Can purely mechanical energy feedback models work?

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]

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

How general is the global density slope-anisotropy inequality?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Momentum Driving: which physical processes dominate AGN feedback?

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 $η=dot{M_outflow}/dot{M_accretion}= (2 ε_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 $η$, and many prior treatments have essentially assumed that $η=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 $ε_w \lesssim 10^{-4}$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the global density slope-anisotropy inequality

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Consistency criteria for generalized Cuddeford systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Crash-testing the CAULDRON code for joint lensing and dynamics analysis of early-type galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Galaxy merging and the Fundamental Plane of elliptical galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

What does the local black hole mass distribution tell us about the evolution of the quasar luminosity function?

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.

0 citations0 reviewsNo code linkedOpen access