Modelling Shear Flows with SPH and Grid Based Methods

Given the importance of shear flows for astrophysical gas dynamics, we study the evolution of the Kelvin-Helmholtz instability (KHI) analytically and numerically. We derive the dispersion relation for the two-dimensional KHI including viscous dissipation. The resulting expression for the growth rate is then used to estimate the intrinsic viscosity of four numerical schemes depending on code-specific as well as on physical parameters. Our set of numerical schemes includes the Tree-SPH code VINE, an alternative SPH formulation developed by Price (2008), and the finite-volume grid codes FLASH and PLUTO. In the first part, we explicitly demonstrate the effect of dissipation-inhibiting mechanisms such as the Balsara viscosity on the evolution of the KHI. With VINE, increasing density contrasts lead to a continuously increasing suppression of the KHI (with complete suppression from a contrast of 6:1 or higher). The alternative SPH formulation including an artificial thermal conductivity reproduces the analytically expected growth rates up to a density contrast of 10:1. The second part addresses the shear flow evolution with FLASH and PLUTO. Both codes result in a consistent non-viscous evolution (in the equal as well as in the different density case) in agreement with the analytical prediction. The viscous evolution studied with FLASH shows minor deviations from the analytical prediction.Comment: 16 pages, 17 figure

Critical properties of spherically symmetric accretion in a fractal medium

Spherically symmetric transonic accretion of a fractal medium has been studied in both the stationary and the dynamic regimes. The stationary transonic solution is greatly sensitive to infinitesimal deviations in the outer boundary condition, but the flow becomes transonic and stable, when its evolution is followed through time. The evolution towards transonicity is more pronounced for a fractal medium than what is it for a continuum. The dynamic approach also shows that there is a remarkable closeness between an equation of motion for a perturbation in the flow, and the metric of an analogue acoustic black hole. The stationary inflow solutions of a fractal medium are as much stable under the influence of linearised perturbations, as they are for the fluid continuum.Comment: 9 pages, 4 figures. Accepted for publication in MNRAS. The definitive version is available at http://www.blackwell-synergy.co

The PN.S Elliptical Galaxy Survey: the dark matter in NGC 4494

We present new Planetary Nebula Spectrograph observations of the ordinary elliptical galaxy NGC 4494, resulting in positions and velocities of 255 PNe out to 7 effective radii (25 kpc). We also present new wide-field surface photometry from MMT/Megacam, and long-slit stellar kinematics from VLT/FORS2. The spatial and kinematical distributions of the PNe agree with the field stars in the region of overlap. The mean rotation is relatively low, with a possible kinematic axis twist outside 1 Re. The velocity dispersion profile declines with radius, though not very steeply, down to ~70 km/s at the last data point. We have constructed spherical dynamical models of the system, including Jeans analyses with multi-component LCDM-motivated galaxies as well as logarithmic potentials. These models include special attention to orbital anisotropy, which we constrain using fourth-order velocity moments. Given several different sets of modelling methods and assumptions, we find consistent results for the mass profile within the radial range constrained by the data. Some dark matter (DM) is required by the data; our best-fit solution has a radially anisotropic stellar halo, a plausible stellar mass-to-light ratio, and a DM halo with an unexpectedly low central density. We find that this result does not substantially change with a flattened axisymmetric model. Taken together with other results for galaxy halo masses, we find suggestions for a puzzling pattern wherein most intermediate-luminosity galaxies have very low concentration halos, while some high-mass ellipticals have very high concentrations. We discuss some possible implications of these results for DM and galaxy formation.Comment: 29 pages, 17 figures. MNRAS, accepte

Kinematic properties of early-type galaxy haloes using planetary nebulae

We present new planetary nebulae (PNe) positions, radial velocities, and magnitudes for 6 early-type galaxies obtained with the Planetary Nebulae Spectrograph, their two-dimensional velocity and velocity dispersion fields. We extend this study to include an additional 10 early-type galaxies with PNe radial velocity measurements available from the literature, to obtain a broader description of the outer-halo kinematics in early-type galaxies. These data extend the information derived from stellar kinematics to typically up to ~8 Re. The combination of photometry, stellar and PNe kinematics shows: i) good agreement between the PNe number density and the stellar surface brightness in the region where the two data sets overlap; ii) good agreement between PNe and stellar kinematics; iii) that the mean rms velocity profiles fall into two groups: with of the galaxies characterized by slowly decreasing profiles and the remainder having steeply falling profiles; iv) a larger variety of velocity dispersion profiles; v) that twists and misalignments in the velocity fields are more frequent at large radii, including some fast rotators; vi) that outer haloes are characterised by more complex radial profiles of the specific angular momentum-related lambda_R parameter than observed within 1Re; vii) that many objects are more rotationally dominated at large radii than in their central parts; and viii) that the halo kinematics are correlated with other galaxy properties, such as total luminosity, isophotal shape, total stellar mass, V/sigma, and alpha parameter, with a clear separation between fast and slow rotators.Comment: 36 pages, 21 figures, revised version for MNRA

Constraining Dark Matter Properties with Gamma-Rays from the Galactic Center with Fermi-LAT

We study the capabilities of the Fermi-LAT instrument on board of the Fermi mission to constrain particle dark matter properties, as annihilation cross section, mass and branching ratio into dominant annihilation channels, with gamma-ray observations from the galactic center. Besides the prompt gamma-ray flux, we also take into account the contribution from the electrons/positrons produced in dark matter annihilations to the gamma-ray signal via inverse Compton scattering off the interstellar photon background, which turns out to be crucial in the case of dark matter annihilations into mu+mu- and e+e- pairs. We study the signal dependence on different parameters like the region of observation, the density profile, the assumptions for the dark matter model and the uncertainties in the propagation model. We also show the effect of the inclusion of a 20% systematic uncertainty in the gamma-ray background. If Fermi-LAT is able to distinguish a possible dark matter signal from the large gamma-ray background, we show that for dark matter masses below ~200 GeV, Fermi-LAT will likely be able to determine dark matter properties with good accuracy.Comment: 38 pages, 13 figures, 4 tables; to match published versio