Phase mixing in MOND

Dissipationless collapses in Modified Newtonian Dynamics (MOND) have been studied by using our MOND particle-mesh N-body code, finding that the projected density profiles of the final virialized systems are well described by Sersic profiles with index m<4 (down to m~2 for a deep-MOND collapse). The simulations provided also strong evidence that phase mixing is much less effective in MOND than in Newtonian gravity. Here we describe "ad hoc" numerical simulations with the force angular components frozen to zero, thus producing radial collapses. Our previous findings are confirmed, indicating that possible differences in radial orbit instability under Newtonian and MOND gravity are not relevant in the present context.Comment: 10 pages, 3 figures. To appear in the Proceedings of the International Workshop "Collective Phenomena in Macroscopic Systems", G. Bertin, R. Pozzoli, M. Rome, and K.R. Sreenivasan, eds., World Scientific, Singapor

ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics

We present a new numerical code, ECHO, based on an Eulerian Conservative High Order scheme for time dependent three-dimensional general relativistic magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at providing a shock-capturing conservative method able to work at an arbitrary level of formal accuracy (for smooth flows), where the other existing GRMHD and GRMD schemes yield an overall second order at most. Moreover, our goal is to present a general framework, based on the 3+1 Eulerian formalism, allowing for different sets of equations, different algorithms, and working in a generic space-time metric, so that ECHO may be easily coupled to any solver for Einstein's equations. Various high order reconstruction methods are implemented and a two-wave approximate Riemann solver is used. The induction equation is treated by adopting the Upwind Constrained Transport (UCT) procedures, appropriate to preserve the divergence-free condition of the magnetic field in shock-capturing methods. The limiting case of magnetodynamics (also known as force-free degenerate electrodynamics) is implemented by simply replacing the fluid velocity with the electromagnetic drift velocity and by neglecting the matter contribution to the stress tensor. ECHO is particularly accurate, efficient, versatile, and robust. It has been tested against several astrophysical applications, including a novel test on the propagation of large amplitude circularly polarized Alfven waves. In particular, we show that reconstruction based on a Monotonicity Preserving filter applied to a fixed 5-point stencil gives highly accurate results for smooth solutions, both in flat and curved metric (up to the nominal fifth order), while at the same time providing sharp profiles in tests involving discontinuities.Comment: 20 pages, revised version submitted to A&

Galactic fountains and gas accretion

Star-forming disc galaxies such as the Milky Way need to accrete \gsim 1 $M_{\odot}$ of gas each year to sustain their star formation. This gas accretion is likely to come from the cooling of the hot corona, however it is still not clear how this process can take place. We present simulations supporting the idea that this cooling and the subsequent accretion are caused by the passage of cold galactic-fountain clouds through the hot corona. The Kelvin-Helmholtz instability strips gas from these clouds and the stripped gas causes coronal gas to condense in the cloud's wake. For likely parameters of the Galactic corona and of typical fountain clouds we obtain a global accretion rate of the order of that required to feed the star formation.Comment: 2 pages, 1 figure, to appear in "Hunting for the Dark: The Hidden Side of Galaxy Formation", Malta, 19-23 Oct. 2009, eds. V.P. Debattista & C.C. Popescu, AIP Conf. Se

Self-consistent stellar dynamical tori

We present preliminary results on a new family of distribution functions that are able to generate axisymmetric, truncated (i.e., finite size) stellar dynamical models characterized by toroidal shapes. The relevant distribution functions generalize those that are known to describe polytropic spheres, for which all the dynamical and structural properties of the system can be expressed in explicit form as elementary functions of the system gravitational potential. The model construction is then completed by a numerical study of the associated Poisson equation. We note that our axisymmetric models can also include the presence of an external gravitational field, such as that produced by a massive disk or by a central mass concentration (e.g., a supermassive black hole).Comment: 4 pages, 2 figures, to appear in "Plasmas in the Laboratory and in the Universe: new insights and new challenges", G. Bertin, D. Farina, R. Pozzoli eds., AIP Conference Proceeding

Fountain-driven gas accretion by the Milky Way

Accretion of fresh gas at a rate of ~ 1 M_{sun} yr^{-1} is necessary in star-forming disc galaxies, such as the Milky Way, in order to sustain their star-formation rates. In this work we present the results of a new hydrodynamic simulation supporting the scenario in which the gas required for star formation is drawn from the hot corona that surrounds the star-forming disc. In particular, the cooling of this hot gas and its accretion on to the disc are caused by the passage of cold galactic fountain clouds through the corona.Comment: 2 pages, 1 figure. To appear in the proceedings of the conference "Assembling the Puzzle of the Milky Way", Le Grand-Bornand 17-22 April 2011, European Physical Journal, editors C. Reyl\'e, A. Robin and M. Schulthei

Dissipationless collapse, weak homology and central cores of elliptical galaxies

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.Comment: Accepted for publication in MNRAS (11 pages, 10 figures

Dissipationless collapses in MOND

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

Three-dimensional magnetic reconnection simulations using the Eulerian Conservative High Order (ECHO) code

Magnetic reconnection and shear driven instabilities are pervasive phenomena in the heliosphere and in astrophysical plasmas in general. Magnetic reconnection and Kelvin-Helmholtz-like instabilities require the use of high-order numerical approximations to study their linear and non-linear evolution. At the same time, in compressible MHD the dynamical activity following reconnection processes leads to formation of discontinuous modes which should be treated by shock-capturing numerical schemes. For this purpose we have designed an Eulerian Conservative High Order (ECHO) code in which, i) explicit diffusivity is taken into account, ii) high-order numerical approximations of flux derivatives are included and iii) shock-capturing algorithms are employed in managing flux discontinuities. This code has been applied successfully in studying the linear and non-linear 3D evolution of the tearing instability and in following the 3D evolution of a current sheet embedded in a sheared flow