It is well known that gas in galaxy discs is highly turbulent, but there is much debate on which mechanism can energetically maintain this turbulence. Among the possible candidates, supernova (SN) explosions are likely the primary drivers but doubts remain on whether they can be sufficient in regions of moderate star formation activity, in particular in the outer parts of discs. In this paper, we measure the SN efficiency $\eta$, namely the fraction of the total SN energy needed to sustain turbulence in galaxies, and verify that SNe can indeed be the sole driving mechanism. The key novelty of our approach is that we take into account the increased turbulence dissipation timescale associated to the flaring in outer regions of gaseous discs. We analyse the distribution and kinematics of HI and CO in 10 nearby star-forming galaxies to obtain the radial profiles of the kinetic energy per unit area, for both the atomic gas and the molecular gas. We use a theoretical model to reproduce the observed energy with the sum of turbulent energy from SNe, as inferred from the observed star formation rate (SFR) surface density, and the gas thermal energy. We find that the observed kinetic energy is remarkably well reproduced by our model across the whole extent of the galactic discs, assuming $\eta$ constant with the galactocentric radius. Taking into account the uncertainties on the SFR surface density and on the atomic gas phase, we obtain that the median SN efficiencies for our sample of galaxies are $\langle \eta_\mathrm{atom} \rangle=0.015_{-0.008}^{+0.018}$ for the atomic gas and $\langle \eta_\mathrm{mol} \rangle = 0.003_{-0.002}^{+0.006}$ for the molecular gas. We conclude that SNe alone can sustain gas turbulence in nearby galaxies with only few percent of their energy and that there is essentially no need for any further source of energy.

The evolution of the structural and kinematic properties of early-type galaxies (ETGs), their scaling relations, as well as their stellar metallicity and age contain precious information on the assembly history of these systems. We present results on the evolution of the stellar mass-velocity dispersion relation of ETGs, focusing in particular on the effects of some selection criteria used to define ETGs. We also try to shed light on the role that in-situ and ex-situ stellar populations have in massive ETGs, providing a possible explanation of the observed metallicity distributions.

We present the first analysis of the stellar content of the structures and substructures identified in the peculiar star-forming galaxy NGC5474, based on Hubble Space Telescope resolved photometry from the LEGUS survey. NGC5474 is a satellite of the giant spiral M101, and is known to have a prominent bulge that is significantly off-set from the kinematic center of the underlying HI and stellar disc. The youngest stars (age~ 100 Myr) trace a flocculent spiral pattern extending out to ~8 kpc from the center of the galaxy. On the other hand intermediate-age (age > 500 Myr) and old (age > 2 Gyr) stars dominate the off-centred bulge and a large substructure residing in the South Western part of the disc and not correlated with the spiral arms (SW over-density). The old age of the stars in the SW over-density suggests that this may be another signature of the dynamical interaction/s that have shaped this anomalous galaxy. We suggest that a fly by with M101, generally invoked as the origin of the anomalies, may not be sufficient to explain all the observations. A more local and more recent interaction may help to put all the pieces of this galactic puzzle together.

Present-day elliptical, spiral and irregular galaxies are large systems made of stars, gas and dark matter. Their properties result from a variety of physical processes that have occurred during the nearly fourteen billion years since the Big Bang. This comprehensive textbook, which bridges the gap between introductory and specialized texts, explains the key physical processes of galaxy formation, from the cosmological recombination of primordial gas to the evolution of the different galaxies that we observe in the Universe today. In a logical sequence, the book introduces cosmology, illustrates the properties of galaxies in the present-day Universe, then explains the physical processes behind galaxy formation in the cosmological context, taking into account the most recent developments in this field. The text ends on how to find distant galaxies with multi-wavelength observations, and how to extract the physical and evolutionary properties based on imaging and spectroscopic data.

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.

Several open questions on galaxy formation and evolution have their roots in the lack of a universal star formation law, that could univocally link the gas properties, e.g. its density, to the star formation rate (SFR) density. In a recent paper, we used a sample of nearby disc galaxies to infer the volumetric star formation (VSF) law, a tight correlation between the gas and the SFR volume densities derived under the assumption of hydrostatic equilibrium for the gas disc. However, due to the dearth of information about the vertical distribution of the SFR in these galaxies, we could not find a unique slope for the VSF law, but two alternative values. In this paper, we use the scale height of the SFR density distribution in our Galaxy adopting classical Cepheids (age$\lesssim 200$ Myr) as tracers of star formation. We show that this latter is fully compatible with the flaring scale height expected from gas in hydrostatic equilibrium. These scale heights allowed us to convert the observed surface densities of gas and SFR into the corresponding volume densities. Our results indicate that the VSF law $\rho_\mathrm{SFR} \propto \rho_\mathrm{gas}^\alpha$ with $\alpha \approx 2$ is valid in the Milky Way as well as in nearby disc galaxies.

Taking advantage of a Bayesian hierarchical inference formalism, we study the evolution of the observed correlation between central stellar velocity dispersion $\sigma_0$ and stellar mass $M_*$ of massive ($M_*\gtrsim 10^{10.5}\,\mathrm{M_\odot}$) early-type galaxies (ETGs) out to redshift $z\approx 2.5$. Collecting ETGs from state-of-the-art literature samples, we consider both a fiducial ($0\lesssim z\lesssim 1.2$) and an extended ($0\lesssim z\lesssim 2.5$) samples. Based on the fiducial sample, we find that at $z\lesssim 1.2$ the $M_*-\sigma_0$ relation is well represented by $\sigma_0\propto M_*^{\beta}(1+z)^{\zeta}$, with $\beta\simeq 0.23$ (at given $z$, $\sigma_0$ increases with $M_*$ as a power law with slope similar to the classical Faber-Jackson relation) and $\zeta\simeq 0.26$ (at given $M_*$, $\sigma_0$ decreases for decreasing $z$, for instance by $\simeq 20\%$ from $z=1$ to $z=0$). The analysis of the extended sample suggests that over the wider redshift range $0\lesssim z\lesssim 2.5$ the slope $\beta$ could be slightly steeper at higher $z$ ($\beta\simeq 0.2+0.2\log(1+z)$ describes the data as well as constant $\beta\simeq 0.23$) and that the redshift-dependence of the normalisation could be stronger ($\zeta\simeq 0.42$). The intrinsic scatter in $\sigma_0$ at given $M_*$ is $\simeq 0.07$ dex for the fiducial sample at $z\lesssim 1.2$, but is found to increase with redshift in the extended sample, being $\simeq 0.11$ dex at $z=2$. Our results suggest that, on average, the velocity dispersion of individual massive ($M_*\gtrsim 5\times 10^{11}M_\odot$) ETGs decreases with time while they evolve from $z\approx 1$ to $z\approx 0$.

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.

We address an issue that arises when self-consistently flattened dynamical stellar systems are constructed by adopting a distribution function (DF) that depends on the action integrals. The velocity distribution at points on the symmetry axis is controlled by the the dependence of the DF on just one action, while at points off the symmetry axis two actions are involved. Consequently, the physical requirement that the velocity distribution evolves continuously in the neighbourhood of the symmetry axis restricts the functional forms of acceptable DFs. An algorithm for conforming to this restriction is presented and used to construct a variety of flattened models.

Dwarf spheroidal galaxies (dSphs) appear to be some of the most dark matter dominated objects in the Universe. Their dynamical masses are commonly derived using the kinematics of stars under the assumption of equilibrium. However, these objects are satellites of massive galaxies (e.g.\ the Milky Way) and thus can be influenced by their tidal fields. We investigate the implication of the assumption of equilibrium focusing on the Sculptor dSph by means of ad-hoc $N$-body simulations tuned to reproduce the observed properties of Sculptor following the evolution along some observationally motivated orbits in the Milky Way gravitational field. For this purpose, we used state-of-the-art spectroscopic and photometric samples of Sculptor's stars. We found that the stellar component of the simulated object is not directly influenced by the tidal field, while $\approx 30\%-60\%$ the mass of the more diffuse DM halo is stripped. We conclude that, considering the most recent estimate of the Sculptor proper motion, the system is not affected by the tides and the stellar kinematics represents a robust tracer of the internal dynamics. In the simulations that match the observed properties of Sculptor, the present-day dark-to-luminous mass ratio is $\approx 6$ within the stellar half-light radius ($\approx0.3$ kpc) and $>50$ within the maximum radius of the analysed dataset ($\approx1.5^\circ\approx2$ kpc).

A new family of self-consistent DF-based models of stellar systems is explored. The stellar component of the models is described by a distribution function (DF) depending on the action integrals, previously used to model the Fornax dwarf spheroidal galaxy (dSph). The stellar component may cohabit with either a dark halo, also described by a DF, or with a massive central black hole. In all cases we solve for the model's self-consistent potential. Focussing on spherically symmetric models, we show how the stellar observables vary with the anisotropy prescribed by the DF, with the dominance and nature of the dark halo, and with the mass of the black hole. We show that precise fits to the observed surface brightness profiles of four globular clusters can be obtained for a wide range of prescribed velocity anisotropies. We also obtain precise fits to the observed projected densities of four dSphs. Finally, we present a three-component model of the Scupltor dSph with distinct DFs for the red and blue horizontal branch stars and the dark matter halo.

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.

We present new dynamical models of dwarf spheroidal galaxies (dSphs) in which both the stellar component and the dark halo are described by analytic distribution functions that depend on the action integrals. In their most general form these distribution functions can represent axisymmetric and possibly rotating stellar systems. Here, as a first application, we model the Fornax dSph, limiting ourselves, for simplicity, to the non rotating, spherical case. The models are compared with state-of-the-art spectroscopic and photometric observations of Fornax, exploiting the knowledge of the line-of-sight velocity distribution of the models and accounting for the foreground contamination from the Milky Way. The model that best fits the structural and kinematic properties of Fornax has a cored dark halo, with core size $r_{\rm c}\simeq1.03$ kpc. The dark-to-luminous mass ratio is $(M_{\rm dm}/M_{\star})|_{R_{\rm eff}}\simeq9.6$ within the effective radius $R_{\rm eff} \simeq 0.62\,$kpc and $(M_{\rm dm}/M_{\star})|_{3 {\rm kpc}} \simeq 144$ within 3 kpc. The stellar velocity distribution is isotropic almost over the full radial range covered by the spectroscopic data and slightly radially anisotropic in the outskirts of the stellar distribution. The dark-matter annihilation $J$-factor and decay $D$-factor are, respectively, $\log_{10}(J$ $[$GeV$^2$ cm$^{-5}])\simeq18.34$ and $\log_{10}(D$ $[$GeV cm$^{-2}])\simeq18.55$, for integration angle $\theta = 0.5^{\circ}$. This cored halo model of Fornax is preferred, with high statistical significance, to both models with a Navarro, Frenk and White dark halo and simple mass-follows-light models.

We study the statistical properties of mergers between central and satellite galaxies in galaxy clusters in the redshift range $0<z<1$, using a sample of dark-matter only cosmological N-body simulations from Le SBARBINE dataset. Using a spherical overdensity algorithm to identify dark-matter haloes, we construct halo merger trees for different values of the over-density $\Delta_c$. While the virial overdensity definition allows us to probe the accretion of satellites at the cluster virial radius $r_{vir}$, higher overdensities probe satellite mergers in the central region of the cluster, down to $\approx 0.06 r_{vir}$, which can be considered a proxy for the accretion of satellite galaxies onto central galaxies. We find that the characteristic merger mass ratio increases for increasing values of $\Delta_c$: more than $60\%$ of the mass accreted by central galaxies since $z\approx 1$ comes from major mergers. The orbits of satellites accreting onto central galaxies tend to be more tangential and more bound than orbits of haloes accreting at the virial radius. The obtained distributions of merger mass ratios and orbital parameters are useful to model the evolution of the high-mass end of the galaxy scaling relations without resorting to hydrodynamic cosmological simulations.

We carry out a systematic investigation of the total mass density profile of massive (Mstar~3e11 Msun) early-type galaxies and its dependence on redshift, specifically in the range 0<z<1. We start from a large sample of SDSS early-type galaxies with stellar masses and effective radii measured assuming two different profiles, de Vaucouleurs and S\'{e}rsic. We assign dark matter haloes to galaxies via abundance matching relations with standard LCDM profiles and concentrations. We then compute the total, mass-weighted density slope at the effective radius gamma', and study its redshift dependence at fixed stellar mass. We find that a necessary condition to induce an increasingly flatter gamma' at higher redshifts, as suggested by current strong lensing data, is to allow the intrinsic stellar profile of massive galaxies to be S\'{e}rsic and the input S\'{e}rsic index n to vary with redshift approximately as n(z)~(1+z)^(-1). This conclusion holds irrespective of the input Mstar-Mhalo relation, the assumed stellar initial mass function, or even the chosen level of adiabatic contraction in the model. Secondary contributors to the observed redshift evolution of gamma' may come from an increased contribution at higher redshifts of adiabatic contraction and/or bottom-light stellar initial mass functions. The strong lensing selection effects we have simulated seem not to contribute to this effect. A steadily increasing S\'{e}rsic index with cosmic time is supported by independent observations, though it is not yet clear whether cosmological hierarchical models (e.g., mergers) are capable of reproducing such a fast and sharp evolution.

We exploit the first \gaia data release to study the properties of the Galactic stellar halo as traced by RR Lyrae. We demonstrate that it is possible to select a pure sample of RR Lyrae using only photometric information available in the Gaia+2MASS catalogue. The final sample contains about 21600 RR Lyrae covering an unprecedented fraction ($\sim60\%$) of the volume of the Galactic inner halo ($\text{R}<28$ kpc). We study the morphology of the stellar halo by analysing the RR Lyrae distribution with parametric and non-parametric techniques. Taking advantage of the uniform all-sky coverage, we test halo models more sophisticated than usually considered in the literature, such as those with varying flattening, tilt and/or offset of the halo with respect to the Galactic disc. A consistent picture emerges: the inner halo is well reproduced by a smooth distribution of stars settled on triaxial ellipsoids. The minor axis is perpendicular to the Milky Ways disc, while the major axis is misaligned by $\sim20^{\circ}$ from the Galactic Y axis. The elongation along the major axis is mild ($\text{p}=1.27$), and the vertical flattening is shown to evolve from a squashed state with $\text{q}\approx0.57$ in the centre to a more spherical $\text{q}\approx0.75$ at the outer edge of our dataset. The density slope is well approximated by a single power-law with exponent $\alpha=-2.96$. Within the range probed, we see no significant evidence for a change of the radial density slope, out of the plane tilt or an offset of the halo with respect to the Galaxy's centre.

We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and S\'{e}rsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a S\'{e}rsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations.

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.

Exploiting the data of the Grism Lens-Amplified Survey from Space (GLASS), we characterize the spatial distribution of star formation in 76 high star forming galaxies in 10 clusters at 0.3< z <0.7. All these galaxies are likely restricted to first infall. In a companion paper we contrast the properties of field and cluster galaxies, whereas here we correlate the properties of H{\alpha} emitters to a number of tracers of the cluster environment to investigate its role in driving galaxy transformations. H{\alpha} emitters are found in the clusters out to 0.5 virial radii, the maximum radius covered by GLASS. The peak of the H{\alpha} emission is offset with respect to the peak of the UV-continuum. We decompose this offsets into a radial and tangential component. The radial compo- nent points away from the cluster center in 60% of the cases, with 95% confidence. The decompositions agree with cosmological simulations, i.e. the H{\alpha} emission offset correlates with galaxy velocity and ram-pressure stripping signatures. Trends between H{\alpha} emitter properties and surface mass density distributions and X-ray emissions emerge only for unrelaxed clusters. The lack of strong correlations with the global environment does not allow us to identify a unique environmental effect originating from the cluster center. In contrast, correla- tions between H{\alpha} morphology and local number density emerge. We conclude that local effects, uncorrelated to the cluster-centric radius, play a more important role in shaping galaxy properties.

Central galaxies (CGs) in galaxy groups and clusters are believed to form and assemble a good portion of their stellar mass at early times, but they also accrete significant mass at late times via galactic cannibalism, that is merging with companion group or cluster galaxies that experience dynamical friction against the common host dark-matter halo. The effect of these mergers on the structure and kinematics of the CG depends not only on the properties of the accreted satellites, but also on the orbital parameters of the encounters. Here we present the results of numerical simulations aimed at estimating the distribution of merging orbital parameters of satellites cannibalized by the CGs in groups and clusters. As a consequence of dynamical friction, the satellites' orbits evolve losing energy and angular momentum, with no clear trend towards orbit circularization. The distributions of the orbital parameters of the central-satellite encounters are markedly different from the distributions found for halo-halo mergers in cosmological simulations. The orbits of satellites accreted by the CGs are on average less bound and less eccentric than those of cosmological halo-halo encounters. We provide fits to the distributions of the central-satellite merging orbital parameters that can be used to study the merger-driven evolution of the scaling relations of CGs.

Dwarf Irregular galaxies (dIrrs) are the smallest stellar systems with extended HI discs. The study of the kinematics of such discs is a powerful tool to estimate the total matter distribution at these very small scales. In this work, we study the HI kinematics of 17 galaxies extracted from the `Local Irregulars That Trace Luminosity Extremes, The HI Nearby Galaxy Survey' (LITTLE THINGS). Our approach differs significantly from previous studies in that we directly fit 3D models (two spatial dimensions plus one spectral dimension) using the software $^\text{3D}$BAROLO, fully exploiting the information in the HI datacubes. For each galaxy we derive the geometric parameters of the HI disc (inclination and position angle), the radial distribution of the surface density, the velocity-dispersion ($\sigma_v$) profile and the rotation curve. The circular velocity (V$_{\text{c}}$), which traces directly the galactic potential, is then obtained by correcting the rotation curve for the asymmetric drift. As an initial application, we show that these dIrrs lie on a baryonic Tully-Fisher relation in excellent agreement with that seen on larger scales. The final products of this work are high-quality, ready-to-use kinematic data ($\textrm{V}_\textrm{c}$ and $\sigma_v$) that we make publicly available. These can be used to perform dynamical studies and improve our understanding of these low-mass galaxies.

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 stellar initial mass function (IMF) of early-type galaxies is the combination of the IMF of the stellar population formed in-situ and that of accreted stellar populations. Using as an observable the effective IMF $\alpha_{IMF}$, defined as the ratio between the true stellar mass of a galaxy and the stellar mass inferred assuming a Salpeter IMF, we present a theoretical model for its evolution as a result of dry mergers. We use a simple dry merger evolution model, based on cosmological $N$-body simulations, together with empirically motivated prescriptions for the IMF to make predictions for how the effective IMF of massive early-type galaxies changes from $z=2$ to $z=0$. We find that the IMF normalization of individual galaxies becomes lighter with time. At fixed velocity dispersion, $\alpha_{IMF}$ is predicted to be constant with redshift. Current constraints on the evolution of the IMF are in slight tension with this prediction, even though systematic uncertainties prevent a conclusive statement. The correlation of $\alpha_{IMF}$ with stellar mass becomes shallower with time, while the correlation between $\alpha_{IMF}$ and velocity dispersion is mostly preserved by dry mergers. We also find that dry mergers can mix the dependence of the IMF on stellar mass and velocity dispersion, making it challenging to infer, from $z=0$ observations of global galactic properties, what is the quantity that is originally coupled with the IMF.

The surface-brightness profiles of galaxies I(R) and the density profiles of dark-matter halos rho(r) are well represented by the same analytic function, named after either S\'ersic, I~exp[-(R/R*)^(1/m)], or Einasto, rho~[exp[-(r/r*)^alpha], where R* and r* are characteristic radii. Systems with high S\'ersic index m (or low Einasto index alpha) have steep central profiles and shallow outer profiles, while systems with low m (or high alpha) have shallow central profiles and steep profiles in the outskirts. We present the results of idealized numerical experiments which suggest that the origin of these profiles can be traced back to the initial density fluctuation field: high-alpha (low-m) systems form in smooth regions via few mergers, while low-alpha (high-m) systems form in clumpy regions via several mergers.

Estimates of the mass distribution and dark-matter (DM) content of dwarf spheroidal galaxies (dSphs) are usually derived under the assumption that the effect of the tidal field of the host galaxy is negligible over the radial extent probed by kinematic data-sets. We assess the implications of this assumption in the specific case of the Fornax dSph by means of N-body simulations of a satellite orbiting around the Milky Way. We consider observationally-motivated orbits and we tailor the initial distributions of the satellite's stars and DM to match, at the end of the simulations, the observed structure and kinematics of Fornax. In all our simulations the present-day observable properties of Fornax are not significantly influenced by tidal effects. The DM component is altered by the interaction with the Galactic field (up to 20% of the DM mass within 1.6 kpc is lost), but the structure and kinematics of the stellar component are only mildly affected even in the more eccentric orbit (more than 99% of the stellar particles remain bound to the dwarf). In the simulations that successfully reproduce Fornax's observables, the dark-to-luminous mass ratio within 1.6 kpc is in the range 5-6, and up to 16-18 if measured within 3 kpc.

Clusters of galaxies are embedded in halos of optically thin, gravitationally stratified, weakly magnetized plasma at the system's virial temperature. Due to radiative cooling and anisotropic heat conduction, such intracluster medium (ICM) is subject to local instabilities, which are combinations of the thermal, magnetothermal and heat-flux-driven buoyancy instabilities. If the ICM rotates significantly, its stability properties are substantially modified and, in particular, also the magnetorotational instability (MRI) can play an important role. We study simple models of rotating cool-core clusters and we demonstrate that the MRI can be the dominant instability over significant portions of the clusters, with possible implications for the dynamics and evolution of the cool cores. Our results give further motivation for measuring the rotation of the ICM with future X-ray missions such as ASTRO-H and ATHENA.

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.

The surface-brightness profiles of galaxies are well described by the S\'ersic law: systems with high S\'ersic index m have steep central profiles and shallow outer profiles, while systems with low m have shallow central profiles and steep outer profiles. R. Cen (2014, ApJL, 790, L24) has conjectured that these profiles arise naturally in the standard cosmological model with initial density fluctuations represented by a Gaussian random field (GRF). We explore and confirm this hypothesis with N-body simulations of dissipationless collapses in which the initial conditions are generated from GRFs with different power spectra. The numerical results show that GRFs with more power on small scales lead to systems with higher m. In our purely dissipationless simulations the S\'ersic index is in the range 2<m<6.5. It follows that systems with S\'ersic index as low as m=2 can be produced by coherent dissipationless collapse, while high-m systems can be obtained if the assembly history is characterized by several mergers. As expected, dissipative processes appear to be required to obtain exponential profiles (m=1).

We investigate the cosmic evolution of the internal structure of massive early-type galaxies over half of the age of the Universe. We perform a joint lensing and stellar dynamics analysis of a sample of 81 strong lenses from the SL2S and SLACS surveys and combine the results with a hierarchical Bayesian inference method to measure the distribution of dark matter mass and stellar IMF across the population of massive early-type galaxies. Lensing selection effects are taken into account. We find that the dark matter mass projected within the inner 5 kpc increases for increasing redshift, decreases for increasing stellar mass density, but is roughly constant along the evolutionary tracks of early-type galaxies. The average dark matter slope is consistent with that of an NFW profile, but is not well constrained. The stellar IMF normalization is close to a Salpeter IMF at $\log{M_*} = 11.5$ and scales strongly with increasing stellar mass. No dependence of the IMF on redshift or stellar mass density is detected. The anti-correlation between dark matter mass and stellar mass density supports the idea of mergers being more frequent in more massive dark matter halos.

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.

Central Galaxies (CGs) in massive halos live in unique environments with formation histories closely linked to that of the host halo. In local clusters they have larger sizes ($R_e$) and lower velocity dispersions (sigma) at fixed stellar mass M_star, and much larger R_e at a fixed $\sigma$ than field and satellite galaxies (non-CGs). Using spectroscopic observations of group galaxies selected from the COSMOS survey, we compare the dynamical scaling relations of early-type CGs and non-CGs at z~0.6, to distinguish possible mechanisms that produce the required evolution. CGs are systematically offset towards larger R_e at fixed $\sigma$ compared to non-CGs with similar M_star. The CG R_e-M_star relation also shows differences, primarily driven by a sub-population (~15%) of galaxies with large $R_e$, while the M_star-sigma relations are indistinguishable. These results are accentuated when double Sersic profiles, which better fit light in the outer regions of galaxies, are adopted. They suggest that even group-scale CGs can develop extended components by these redshifts that can increase total $R_e$ and M_star estimates by factors of ~2. To probe the evolutionary link between our sample and cluster CGs, we also analyze two cluster samples at z~0.6 and z~0. We find similar results for the more massive halos at comparable z, but much more distinct CG scaling relations at low-z. Thus, the rapid, late-time accretion of outer components, perhaps via the stripping and accretion of satellites, would appear to be a key feature that distinguishes the evolutionary history of CGs.

Simulations of the clustering of cold dark matter yield dark-matter halos that have central density cusps, but observations of totally dark-matter dominated dwarf spheroidal galaxies imply that they do not have cuspy central density profiles. We use analytic calculations and numerical modelling to argue that whenever stars form, central density cusps are likely to be erased. Gas that accumulates in the potential well of an initially cuspy dark-matter halo settles into a disc. Eventually the surface density of the gas exceeds the threshold for fragmentation into self-gravitating clouds. The clouds are massive enough to transfer energy to the dark-matter particles via dynamical friction on a short time-scale. The halo's central cusp is heated to form a core with central logarithmic density slope gamma=0 before stellar feedback makes its impact. Since star formation is an inefficient process, the clouds are disrupted by feedback when only a small fraction of their mass has been converted to stars, and the dark matter dominates the final mass distribution.

A long-standing question is whether radiative cooling can lead to local condensations of cold gas in the hot atmospheres of galaxies and galaxy clusters. We address this problem by studying the nature of local instabilities in rotating, stratified, weakly magnetized, optically thin plasmas in the presence of radiative cooling and anisotropic thermal conduction. For both axisymmetric and non-axisymmetric linear perturbations we provide the general equations that can be applied locally to specific systems to establish whether they are unstable and, in case of instability, to determine the kind of evolution (monotonically growing or over-stable) and the growth rates of unstable modes. We present results for models of rotating plasmas representative of Milky Way-like galaxy coronae and cool-cores of galaxy clusters. It is shown that the unstable modes arise from a combination of thermal, magnetothermal, magnetorotational and heat-flux-driven buoyancy instabilities. Local condensation of cold clouds tends to be hampered in cluster cool cores, while it is possible under certain conditions in rotating galactic coronae. If the magnetic field is sufficiently weak the magnetorotational instability is dominant even in these pressure supported systems.

Several studies have suggested that the observed size evolution of massive early-type galaxies (ETGs) can be explained as a combination of dry mergers and progenitor bias, at least since z\sim1. In this paper we carry out a new test of the dry-merger scenario based on recent lensing measurements of the evolution of the mass density profile of ETGs. We construct a theoretical model for the joint evolution of the size and mass density profile slope \gamma' driven by dry mergers occurring at rates given by cosmological simulations. Such dry-merger model predicts a strong decrease of \gamma' with cosmic time, inconsistent with the almost constant \gamma' inferred from observations in the redshift range 0<z<1. We then show with a simple toy model that a modest amount of cold gas in the mergers -- consistent with the upper limits on recent star formation in ETGs -- is sufficient to reconcile the model with measurements of \gamma'. By fitting for the amount of gas accreted during mergers, we find that models with dissipation are consistent with observations of the evolution in both size and density slope, if \sim4% of the total final stellar mass arises from the gas accreted since z\sim1. Purely dry merger models are ruled out at > 99% CL. We thus suggest a scenario where the outer regions of massive ETGs grow by accretion of stars and dark matter, while small amounts of dissipation and nuclear star formation conspire to keep the mass density profile constant and approximately isothermal.

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.

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 optical and near infrared spectroscopy obtained at Keck, VLT, and Gemini for a sample of 36 secure strong gravitational lens systems and 17 candidates identified as part of the SL2S survey. The deflectors are massive early-type galaxies in the redshift range z_d=0.2-0.8, while the lensed sources are at z_s=1-3.5. We combine this data with photometric and lensing measurements presented in the companion paper III and with lenses from the SLACS and LSD surveys to investigate the cosmic evolution of the internal structure of massive early-type galaxies over half the age of the universe. We study the dependence of the slope of the total mass density profile \gamma' (\rho(r)\propto r^{-\gamma'}) on stellar mass, size, and redshift. We find that two parameters are sufficent to determine \gamma' with less than 6% residual scatter. At fixed redshift, \gamma' depends solely on the surface stellar mass density \partial \gamma'/ \partial \Sigma_*=0.38\pm 0.07, i.e. galaxies with denser stars also have steeper slopes. At fixed M_* and R_{eff}, \gamma' depends on redshift, in the sense that galaxies at a lower redshift have steeper slopes (\partial \gamma' / \partial z = -0.31\pm 0.10). However, the mean redshift evolution of \gamma' for an individual galaxy is consistent with zero d\gamma'/dz=-0.10\pm0.12. This result is obtained by combining our measured dependencies of \gamma' on z,M_*,R_{eff} with the evolution of the R_{eff}-M_* taken from the literature, and is broadly consistent with current models of the formation and evolution of massive early-type galaxies. Detailed quantitative comparisons of our results with theory will provide qualitatively new information on the detailed physical processes at work.

We study simple models of massive galaxy clusters in which the intracluster medium (ICM) rotates differentially in equilibrium in the cluster gravitational potential. We obtain the X-ray surface brightness maps, evaluating the isophote flattening due to the gas rotation. Using a set of different rotation laws, we put constraint on the amplitude of the rotation velocity, finding that rotation curves with peak velocity up to \sim 600 km s^-1 are consistent with the ellipticity profiles of observed clusters. We convolve each of our models with the instrument response of the X-ray Calorimeter Spectrometer on board the ASTRO-H to calculate the simulated X-ray spectra at different distance from the X-ray centre. We demonstrate that such an instrument will allow us to measure rotation of the ICM in massive clusters, even with rotation velocities as low as \sim 100 km s^-1

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.

Clusters of galaxies are excellent locations to probe the distribution of baryons and dark matter (DM) over a wide range of scales. We study a sample of seven massive, relaxed galaxy clusters with centrally-located brightest cluster galaxies (BCGs) at z=0.2-0.3. Using the observational tools of strong and weak gravitational lensing, combined with resolved stellar kinematics within the BCG, we measure the total radial density profile, comprising both dark and baryonic matter, over scales of ~3-3000 kpc. Lensing-derived mass profiles typically agree with independent X-ray estimates within ~15%, suggesting that departures from hydrostatic equilibrium are small and that the clusters in our sample (except A383) are not strongly elongated along the line of sight. The inner logarithmic slope gamma_tot of the total density profile measured over r/r200=0.003-0.03, where rho_tot ~ r^(-gamma_tot), is found to be nearly universal, with a mean <gamma_tot> = 1.16 +- 0.05 (random) +0.05-0.07 (systematic) and an intrinsic scatter of < 0.13 (68% confidence). This is further supported by the very homogeneous shape of the observed velocity dispersion profiles, obtained via Keck spectroscopy, which are mutually consistent after a simple scaling. Remarkably, this slope agrees closely with numerical simulations that contain only dark matter, despite the significant contribution of stellar mass on the scales we probe. The Navarro-Frenk-White profile characteristic of collisionless cold dark matter is a better description of the total mass density at radii >~ 5-10 kpc than that of dark matter alone. Hydrodynamical simulations that include baryons, cooling, and feedback currently provide a poorer match. We discuss the significance of our findings for understanding the assembly of BCGs and cluster cores, particularly the influence of baryons on the inner DM halo. [abridged]

A few years ago we started an observational campaign aimed at the thorough study of the massive and remote globular cluster NGC2419. We have used the collected data, e.g., to test alternative theories of gravitation, to constrain the stellar M/L ratio by direct analysis of the observed luminosity function, and to search for Dark Matter within the cluster. Here we present some recent results about (a) the peculiar abundance pattern that we observed in a sample of cluster giants, and (b) newly found photometric evidence for the presence of multiple populations in the cluster. In particular, from new deep and accurate uVI LBT photometry, we find that the color spread on the Red Giant Branch is significantly larger than the observational errors both in V-I and u-V, and that the stars lying to the blue of the RGB ridge line are more concentrated toward the center of the cluster than those lying to the red of the ridge line.

We use recently published measurements of the kinematics, surface brightness and stellar mass-to-light ratio of the globular cluster NGC 2419 to examine the possibility that this Galactic halo satellite is embedded in a low-mass dark matter halo. NGC 2419 is a promising target for such a study, since its extreme Galactocentric distance and large mass would have greatly facilitated the retention of dark matter. A Markov-Chain Monte Carlo approach is used to investigate composite dynamical models containing a stellar and a dark matter component. We find that it is unlikely that a significant amount of dark matter (less than approx. 6% of the luminous mass inside the tidal limit of the cluster) can be present if the stars follow an anisotropic Michie model and the dark matter a double power law model. However, we find that more general models, derived using a new technique we have developed to compute non-parametric solutions to the spherical Jeans equation, suggest the presence of a significant dark matter fraction (approximately twice the stellar mass). Thus the presence of a dark matter halo around NGC 2419 cannot be fully ruled out at present, yet any dark matter within the 10 arcmin visible extent of the cluster must be highly concentrated and cannot exceed 1.1x10^6 Solar masses (99% confidence), in stark contrast to expectations for a plausible progenitor halo of this structure.

The thermal stability of a weakly magnetized, rotating, stratified, optically thin plasma is studied by means of linear-perturbation analysis. We derive dispersion relations and criteria for stability against axisymmetric perturbations that generalize previous results on either non-rotating or unmagnetized fluids. The implications for the hot atmospheres of galaxies and galaxy clusters are discussed.

Early-type galaxies (ETGs) are observed to be more compact at z>2 than in the local Universe. Remarkably, much of this size evolution appears to take place in a short (1.8 Gyr) time span between z=2.2 and z=1.3, which poses a serious challenge to hierarchical galaxy formation models where mergers occurring on a similar timescale are the main mechanism for galaxy growth. We compute the merger-driven redshift evolution of stellar mass Mstar\propto(1+z)^aM, half-mass radius Re\propto(1+z)^aR and velocity-dispersion sigma0\propto(1+z)^asigma predicted by concordance Lambda cold dark matter for a typical massive ETG in the redshift range z=1.3-2.2. Neglecting dissipative processes, and thus maximizing evolution in surface density, we find -1.5<aM<-0.6, -1.9<aR<-0.7 and 0.06<asigma<0.22, under the assumption that the accreted satellites are spheroids. It follows that the predicted z=2.2 progenitors of z=1.3 ETGs are significantly less compact (on average a factor of 2 larger Re at given Mstar) than the quiescent galaxies observed at z>2. Furthermore, we find that the scatter introduced in the size-mass correlation by the predicted merger-driven growth is difficult to reconcile with the tightness of the observed scaling law. We conclude that - barring unknown systematics or selection biases in the current measurements - minor and major mergers with spheroids are not sufficient to explain the observed size growth of ETGs within the standard model.

We present the analysis of a large sample of early-type galaxies (ETGs) at 0<z<3 aimed at tracing the cosmic evolution of their size and compare it with a model of pure dissipationless (dry) merging in the LambdaCDM framework. The effective radius R_e depends on stellar mass M as R_e(M) \propto M}^{alpha} with alpha ~ 0.5 at all redshifts. The redshift evolution of the mass- or SDSS-normalized size can be reproduced as \propto (1+z)^beta with beta ~ -1, with the most massive ETGs possibly showing the fastest evolutionary rate (beta ~ -1.4). This size evolution slows down significantly to beta ~ -0.6 if the ETGs at z>2 are removed from the sample, suggesting an accelerated increase of the typical sizes at z>2, especially for the ETGs with the largest masses. A pure dry merging LambdaCDM model is marginally consistent with the average size evolution at 0<z<1.7, but predicts descendants too compact for z>2 progenitor ETGs. This opens the crucial question on what physical mechanism can explain the accelerated evolution at z>2, or whether an unclear observational bias is partly responsible for that.

We present the results of a detailed analysis of the projected velocity dispersion of the globular cluster Palomar 14 performed using recent high-resolution spectroscopic data and extensive Monte Carlo simulations. The comparison between the data and a set of dynamical models (differing in fraction of binaries, degree of anisotropy, mass-to-light ratio M/L, cluster orbit and theory of gravity) shows that the observed velocity dispersion of this stellar system is well reproduced by Newtonian models with a fraction of binaries f_b<30% and a M/L compatible with the predictions of stellar evolution models. Instead, models computed with a large fraction of binaries systematically overestimate the cluster velocity dispersion. We also show that, across the parameter space sampled by our simulations, models based on the Modified Newtonian Dynamics theory can be reconciled with observations only assuming values of M/L lower than those predicted by stellar evolution models under standard assumptions.

We present an analysis of the globular cluster NGC 2419, using a polytropic model in Modified Newtonian Dynamics (MOND) to reproduce recently published high quality data of the structure and kinematics of the system. We show that a specific MOND polytropic model of NGC 2419 suggested by a previous study can be completely ruled out by the data. Furthermore, the highest likelihood fit polytrope in MOND is a substantially worse model (by a factor of approximately 5000) than a Newtonian Michie model studied in an earlier contribution by Ibata et al. (2011). We conclude that the structure and dynamics of NGC 2419 favor Newtonian dynamics and do indeed challenge the MOND theory.

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.

Dry (i.e. dissipationless) merging has been proposed as the main driver of the observed size evolution of early-type galaxies (ETGs). The actual role of this mechanism is questioned by the tightness of the local stellar mass-size relation of ETGs. Combining this observed scaling law with simple merging models, which should bracket cosmologically motivated merging histories, we draw the following conclusions: 1) local massive ETGs can have assembled at most ~45% of their stellar mass via dry mergers; 2) extreme fine tuning is required for this to be the case.

We present the analysis of a kinematic data set of stars in the globular cluster NGC 2419, taken with Keck/DEIMOS. Combined with a reanalysis of deep HST and Subaru imaging data, which provide an accurate luminosity profile of the cluster, we investigate the validity of a large set of dynamical models of the system, which are checked for stability via N-body simulations. We find that isotropic models in either Newtonian or Modified Newtonian Dynamics (MOND) are ruled out with extremely high confidence. However, a simple Michie model in Newtonian gravity with anisotropic velocity dispersion provides an excellent representation of the luminosity profile and kinematics. In contrast, with MOND we find that Michie models that reproduce the luminosity profile either over-predict the velocity dispersion on the outskirts of the cluster if the mass to light ratio is kept at astrophysically-motivated values, or else they under-predict the central velocity dispersion if the mass to light ratio is taken to be very small. We find that the best Michie model in MOND is a factor of 10000 less likely than the Newtonian model that best fits the system. A likelihood ratio of 350 is found when we investigate more general models by solving the Jeans equation with a Markov-Chain Monte Carlo scheme. We verified with N-body simulations that these results are not significantly different when the MOND external field effect is accounted for. If the assumptions that the cluster is in dynamical equilibrium, spherical, not on a peculiar orbit, and possesses a single dynamical tracer population of constant M/L are correct, we conclude that the present observations provide a very severe challenge for MOND. [abridged]

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

The thermal stability of rotating, stratified, unmagnetized atmospheres is studied by means of linear-perturbation analysis, finding stability, overstability or instability, depending on the properties of the gas distribution, but also on the nature of the perturbations. In the relevant case of distributions with outward-increasing specific entropy and angular momentum, axisymmetric perturbations grow exponentially, unless their wavelength is short enough that they are damped by thermal conduction; non-axisymmetric perturbations typically undergo overstable oscillations in the limit of zero conductivity, but are effectively stabilized by thermal conduction, provided rotation is differential. To the extent that the studied models are representative of the poorly constrained hot atmospheres of disc galaxies, these results imply that blob-like, cool overdensities are unlikely to grow in galactic coronae, suggesting an external origin for the high-velocity clouds of the Milky Way.

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.

A substantial fraction of the baryons of disk galaxies like the Milky Way is expected to reside in coronae of gas at the virial temperature. This is the only realistic reservoir of gas available to feed star formation in the disks. The way this feeding occurs depends crucially on whether galactic coronae can fragment into cool clouds via thermal instability. Here I summarize arguments suggesting that galactic coronae are not prone to thermal instability, and I briefly discuss the implications for galaxy-formation models and for the origin of the high-velocity clouds of the Milky Way.

The characteristic size of early-type galaxies (ETGs) of given stellar mass is observed to increase significantly with cosmic time, from redshift z>2 to the present. A popular explanation for this size evolution is that ETGs grow through dissipationless ("dry") mergers, thus becoming less compact. Combining N-body simulations with up-to-date scaling relations of local ETGs, we show that such an explanation is problematic, because dry mergers do not decrease the galaxy stellar-mass surface-density enough to explain the observed size evolution, and also introduce substantial scatter in the scaling relations. Based on our set of simulations, we estimate that major and minor dry mergers increase half-light radius and projected velocity dispersion with stellar mass (M) as M^(1.09+/-0.29) and M^(0.07+/-0.11), respectively. This implies that: 1) if the high-z ETGs are indeed as dense as estimated, they cannot evolve into present-day ETGs via dry mergers; 2) present-day ETGs cannot have assembled more than ~45% of their stellar mass via dry mergers. Alternatively, dry mergers could be reconciled with the observations if there was extreme fine tuning between merger history and galaxy properties, at variance with our assumptions. Full cosmological simulations will be needed to evaluate whether this fine-tuned solution is acceptable.

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.

If the dark matter sector in the universe is composed by metastable particles, galaxies and galaxy clusters are expected to undergo significant secular evolution from high to low redshift. We show that the decay of dark matter, with a lifetime compatible with cosmological constraints, can be at the origin of the observed evolution of the Tully-Fisher relation of disk galaxies and alleviate the problem of the size-evolution of elliptical galaxies, while being consistent with the current observational constraints on the gas fraction of clusters of galaxies.

We test the modified Newtonian dynamics (MOND) theory with the velocity-dispersion profiles of Galactic globular clusters populating the outermost region of the Milky Way halo, where the Galactic acceleration is lower than the characteristic MOND acceleration a_0. For this purpose, we constructed self-consistent, spherical models of stellar systems in MOND, which are the analogues of the Newtonian King models. The models are spatially limited, reproduce well the surface-brightness profiles of globular clusters, and have velocity-dispersion profiles that differ remarkably in shape from the corresponding Newtonian models. We present dynamical models of six globular clusters, which can be used to efficiently test MOND with the available observing facilities. A comparison with recent spectroscopic data obtained for NGC2419 suggests that the kinematics of this cluster might be hard to explain in MOND.

The most evident features of colour-magnitude diagrams of galaxies are the red sequence of quiescent galaxies, extending up to the brightest elliptical galaxies, and the blue cloud of star-forming galaxies, which is truncated at a luminosity L~L*. The truncation of the blue cloud indicates that in the most massive systems star formation must be quenched. For this to happen the virial-temperature galactic gas must be kept hot and any accreted cold gas must be heated. The elimination of accreted cold gas can be due to thermal evaporation by the hot interstellar medium, which in turn is prevented from cooling by feedback from active galactic nuclei.

Dissipationless (gas-free or "dry") mergers have been suggested to play a major role in the formation and evolution of early-type galaxies, particularly in growing their mass and size without altering their stellar populations. We perform a new test of the dry merger hypothesis by comparing N-body simulations of realistic systems to empirical constraints provided by recent studies of lens early-type galaxies. We find that major and minor dry mergers: i) preserve the nearly isothermal structure of early-type galaxies within the observed scatter; ii) do not change more than the observed scatter the ratio between total mass M and "virial" mass R_e*sigma/2G (where R_e is the half-light radius and sigma the projected velocity dispersion); iii) increase strongly galaxy sizes [as M^(0.85+/-0.17)] and weakly velocity dispersions [as M^(0.06+/-0.08)] with mass, thus moving galaxies away from the local observed M-R_e and M-sigma relations; iv) introduce substantial scatter in the M-R_e and M-sigma relations. Our findings imply that, unless there is a high degree of fine tuning of the mix of progenitors and types of interactions, present-day massive early-type galaxies cannot have assembled more than ~50% of their mass, and increased their size by more than a factor ~1.8, via dry merging.

There are two families of luminous elliptical galaxies: cusp galaxies, with steep central surface-brightness profiles, and core galaxies, whose surface-brightness profiles have flat central cores. Thermal evaporation of accreted cold gas by the hot interstellar medium may be at the origin of this cusp-core dichotomy: in less massive (hot-gas poor) galaxies central cores are likely to be refilled by central starbursts following cold gas infall, while in more massive (hot-gas rich) galaxies most cold gas is eliminated and central cores survive. This scenario is consistent with the observation that cusp and core galaxies differ systematically in terms of optical luminosity, X-ray gas content, age of the central stellar population, and properties of the active galactic nucleus.

We examine the proposal that the HI "high-velocity" clouds (HVCs) surrounding the Milky Way and other disc galaxies form by condensation of the hot galactic corona via thermal instability. Under the assumption that the galactic corona is well represented by a non-rotating, stratified atmosphere, we find that for this formation mechanism to work the corona must have an almost perfectly flat entropy profile. In all other cases the growth of thermal perturbations is suppressed by a combination of buoyancy and thermal conduction. Even if the entropy profile were nearly flat, cold clouds with sizes smaller than 10 kpc could form in the corona of the Milky Way only at radii larger than 100 kpc, in contradiction with the determined distances of the largest HVC complexes. Clouds with sizes of a few kpc can form in the inner halo only in low-mass systems. We conclude that unless even slow rotation qualitatively changes the dynamics of a corona, thermal instability is unlikely to be a viable mechanism for formation of cold clouds around disc 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.

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.

Massive early-type galaxies are observed to lie on the Mass Plane (MP), a two-dimensional manifold in the space of effective radius R_e, projected mass M_p (measured via strong gravitational lensing) and projected velocity dispersion sigma within R_e/2. The MP is less `tilted' than the Fundamental Plane, and the two have comparable associated scatter. This means that c_e2=2*G*M_p/(R_e*sigma^2) is a nearly universal constant in the range sigma=175-400 km/s. This finding can be used to constrain the mass distribution and internal dynamics of early-type galaxies. We find that a relatively wide class of spherical galaxy models has values of c_e2 in the observed range, because c_e2 is not very strongly sensitive to the mass distribution and orbital anisotropy. If the total mass distribution is isothermal, a broad range of stellar luminosity profile and anisotropy is consistent with the observations, while NFW dark-matter halos require more fine tuning of the stellar mass fraction, luminosity profile and anisotropy. If future data can cover a broader range of masses, the MP could be seen to be tilted and the value of any such tilt would provide a discriminant between models for the total mass-density profile of the galaxies. [Abridged]

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 describe the numerical code N-MODY, a parallel particle-mesh code for collisionless N-body simulations in modified Newtonian dynamics (MOND). N-MODY is based on a numerical potential solver in spherical coordinates that solves the non-linear MOND field equation, and is ideally suited to simulate isolated stellar systems. N-MODY can be used also to compute the MOND potential of arbitrary static density distributions. A few applications of N-MODY indicate that some astrophysically relevant dynamical processes are profoundly different in MOND and in Newtonian gravity with dark matter.

In colour-magnitude diagrams most galaxies fall in either the ``blue cloud'' or the ``red sequence'', with the red sequence extending to significantly brighter magnitudes than the blue cloud. The bright-end of the red sequence comprises elliptical galaxies (Es) with boxy isophotes and luminosity profiles with shallow central cores, while fainter Es have disky isophotes and power-law inner surface-brightness (SB) profiles. An analysis of published data reveals that the centres of galaxies with power-law central SB profiles have younger stellar populations than the centres of cored galaxies. We argue that thermal evaporation of cold gas by virial-temperature gas plays an important role in determining these phenomena. In less massive galaxies, thermal evaporation is not very efficient, so significant amounts of cold gas can reach the galaxy centre and fill a central core with newly formed stars, consistent with the young stellar ages of the cusps of Es with power-law SB profiles. In more massive galaxies, cold gas is evaporated within a dynamical time, so star formation is inhibited, and a core in the stellar density profile produced by dissipationless dynamics cannot be refilled. The different observed properties of AGN in higher-mass and lower-mass ellipticals are also explained because in the former the central black holes invariably accrete hot gas, while in the latter they typically accrete cold gas. An important consequence of our results is that at the present time there cannot be blue, star-forming galaxies in the most massive galactic halos, consistent with the observed truncation of the blue cloud at L*. [abridged]

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.

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 analysed the long-term variability of four microquasars (GRS 1915+105, Cyg X-1, Cyg X-3, and Sco X-1) in radio and X rays. The results of our analysis indicate the existence of two distinct modes of energy output, which we refer to as the `coupled' mode and the `flaring' mode. The coupled mode is responsible for mildly fluctuating, flat-spectrum radio emission, coupled with the X-ray emission; the flaring mode produces powerful, steep-spectrum radio flares, with no significant counterpart in X rays. We find that the fraction of time spent by a typical microquasar in the flaring mode is similar to the fraction of quasars that are radio-loud. This is consistent with the hypothesis that radio-loudness of quasars is a function of the epoch at which the source is observed.

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.

The low-frequency power spectra of the X-ray and radio emission from four microquasars suggest that two distinct modes of energy output are at work: (i) the `coupled' mode in which the X-ray and radio luminosities are closely coupled and vary only weakly, and (ii) the `flaring' mode, which dramatically boosts the radio luminosity but makes no impact on the X-ray luminosity. The systems are in the flaring mode only a few percent of the time. However, flares completely dominate the power spectrum of radio emission, with the consequence that sources in which the flaring mode occurs, such as GRS 1915+105 and Cyg X-3, have radio power spectra that lie more than an order of magnitude above the corresponding X-ray power spectra. Of the four microquasars for which we have examined data, in only one, Cyg X-1, is the flaring mode seemingly inactive. While Cyg X-1 is a black-hole candidate, one of the three flaring sources, Sco X-1, is a neutron star. Consequently, it is likely that both modes are driven by the accretion disk rather than black-hole spin. Radio imaging strongly suggests that the flaring mode involves relativistic jets. A typical microquasar is in the flaring mode a few percent of the time, which is similar to the fraction of quasars that are radio loud. Thus there may be no essential difference between radio-loud and radio-quiet quasars; radio loudness may simply be a function of the epoch at which the source is observed.

There is increasing evidence that AGN mechanical feedback is important in the energetics of cooling flows in galaxies and galaxy clusters. We investigate the implications of the variability of AGN mechanical luminosity L_m on observations of cooling flows and radio galaxies in general. It is natural to assume that l=ln(L_m/L_x) is a Gaussian process. Then L_m will be log-normally distributed at fixed cooling luminosity L_x, and the variance in a measure of L_m will increase with the time-resolution of the measure. We test the consistency of these predictions with existing data. These tests hinge on the power spectrum of l(t). Monitoring of Seyfert galaxies combined with estimates of the duty cycle of quasars imply flicker noise spectra, similar to those of microquasars. We combine a sample of sources in cooling flows that have cavities with the assumption that the average mechanical luminosity of the AGN equals L_x. Given that the mechanical luminosities are characterized by flicker noise, we find that their spectral amplitudes lie between the estimated amplitudes of quasars and the measured values for the radio luminosities of microquasars. The model together with the observation that powerful radio galaxies lie within a narrow range in optical luminosity, predicts the luminosity function of radio galaxies, in agreement with observations. Forthcoming radio surveys will test the prediction that the luminosity function turns over at about the smallest luminosities so far probed. [Abridged]

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.

We determine the critical size l_crit of a filament of cold (T~10^4 K) gas that is in radiative equilibrium with X-ray emitting gas at temperatures T_out~10^6 - 10^8 K. Filaments smaller than l_crit will be rapidly evaporated, while longer ones will induce the condensation of the ambient medium. At fixed pressure P, l_crit increases as T_out^(11/4), while at fixed T_out it scales as 1/P. It scales as f^(1/2), where f is the factor by which the magnetic field depresses the thermal conductivity below Spitzer's benchmark value. For plausible values of f, l_crit is similar to the lengths of observed filaments. In a cluster such as Perseus, the value of l_crit increases by over an order of magnitude between the centre and a radius of 100 kpc. If the spectrum of seed filament lengths l is strongly falling with l, as is natural, then these results explain why filaments are only seen within a few kiloparsecs of the centres of clusters, and are not seen in clusters that have no cooling flow. We calculate the differential emission measure as a function of temperature for the interface between filaments and ambient gas of various temperatures. We discuss the implications of our results for the origin of the galaxy luminosity function.

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.

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

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.

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.