The Nuker model for galactic nuclei

The Nuker profile, characterised by an inner and outer power-law profile smoothly merged around a break radius, is a very popular model to describe the surface brightness profile of galactic nuclei. A disadvantage of this model for dynamical studies is that the spatial density distribution that corresponds to this surface brightness profile cannot be written in terms of elementary or regular special functions. We derive a compact and elegant analytical expression for the density of the Nuker model, based the Mellin integral transform method. We use this expression to discuss the general behaviour and asymptotic expansion of the density. We also discuss the special subclass of Nuker models with an infinitely sharp break and demonstrate that these models are always characterised by non-monotonous and hence unphysical density profile. We extend our study to the dynamical structure of spherical isotropic galactic nuclei with a Nuker surface brightness profile. Based on this analysis, we extend and refine the classification of spherical isotropic galactic nuclei introduced by Tremaine et al. (1994, AJ, 107, 634). We demonstrate that both the inner density slope and the sharpness of the break between the inner and outer profiles critically determine the consistency and stability of the Nuker models

Stellar systems following the R^1/m luminosity law, IV : the total energy and the central concentration of galaxies

We expand our previous analytical and numerical studies of the family of Sérsic models, which are routinely used to describe early-type galaxies and the bulges of spiral galaxies. In particular, we focus on the total energy budget, an important dynamical property that has not been discussed in detail in previous works. We use two different methods to calculate the total energy for the Sérsic model family that result in two independent expressions that can be used along the entire sequence of Sérsic models. We use these expressions to investigate whether the Spitzer concentration index is a reliable measure for the intrinsic 3D concentration of galaxies, and we conclude that it is not a very useful measure for the central concentration. The popular Third Galaxy Concentration index, on the other hand, is shown to be a reliable measure for the intrinsic 3D concentration, even though it is based on the surface brightness distribution and not on the intrinsic 3D density

Stellar systems following the R1/mR^{1/m} luminosity law. III. Photometric, intrinsic, and dynamical properties for all S\'ersic indices

The S\'ersic or $R^{1/m}$ model has become the de facto standard model to describe the surface brightness profiles of early-type galaxies and the bulges of spiral galaxies. The photometric, intrinsic, and dynamical properties of this model have been investigated, but mainly for fairly large S\'ersic indices $m$. For small values of $m$, appropriate for low-mass and dwarf ellipticals, a detailed investigation of these properties is still lacking. In this study, we used a combination of numerical and analytical techniques to investigate the S\'ersic model over the entire range of S\'ersic parameters, focussing on the small $m$ regime, where a number of interesting and surprising properties are found. For all values $m<1$, the model is characterised by a finite central luminosity density, and for $m<\tfrac12$, even a central depression in the luminosity density profile. This behaviour translates to the dynamical properties: we show that all S\'ersic models with $m \geqslant\tfrac12$ can be supported by an isotropic velocity dispersion tensor, and that these isotropic models are stable to both radial and non-radial perturbations. The models with $m < \tfrac12$, on the other hand, cannot be supported by an isotropic velocity dispersion tensor.Comment: 10 pages, 5 figures, accepted for publication in A&

The Failure of Monte Carlo Radiative Transfer at Medium to High Optical Depths

Computer simulations of photon transport through an absorbing and/or scattering medium form an important research tool in astrophysics. Nearly all software codes performing such simulations for three-dimensional geometries employ the Monte Carlo radiative transfer method, including various forms of biasing to accelerate the calculations. Because of the probabilistic nature of the Monte Carlo technique, the outputs are inherently noisy, but it is often assumed that the average values provide the physically correct result. We show that this assumption is not always justified. Specifically, we study the intensity of radiation penetrating an infinite, uniform slab of material that absorbs and scatters the radiation with equal probability. The basic Monte Carlo radiative transfer method, without any biasing mechanisms, starts to break down for transverse optical depths above ~20 because so few of the simulated photon packets reach the other side of the slab. When including biasing techniques such as absorption/scattering splitting and path length stretching, the simulated photon packets do reach the other side of the slab but the biased weights do not necessarily add up to the correct solution. While the noise levels seem to be acceptable, the average values sometimes severely underestimate the correct solution. Detecting these anomalies requires the judicious application of statistical tests, similar to those used in the field of nuclear particle transport, possibly in combination with convergence tests employing consecutively larger numbers of photon packets. In any case, for transverse optical depths above ~75 the Monte Carlo methods used in our study fail to solve the one-dimensional slab problem, implying the need for approximations such as a modified random walk.Comment: Accepted for publication in the ApJ; 13 pages, 6 figure

Dynamical models for dusty disk galaxies

Disk galaxies contain a large amount of interstellar dust, which affects the projection of kinematic quantities. We investigate in detail the effects of dust extinction on the mean projected velocity and the projected velocity dispersion. We use our results to construct a general strategy to determine the dynamical structure of disk galaxies, with the aim to constrain their mass distribution and dynamical history.Comment: to be published in the proceedings of "Galaxy Disks and Disk Galaxies", Funes J.G. and Corsini E.M. eds., ASP Conference Serie

Tracing the relation between black holes and dark haloes

We present new velocity dispersion measurements for a set of 12 spiral galaxies and use them to derive a more accurate V_c - sigma relation which holds for a wide morphological range of galaxies. Combined with the M_BH - sigma relation, this relation can be used as a tool to estimate supermassive black hole (SMBH) masses by means of the asymptotic circular velocity. Together with the Tully-Fisher relation, it serves as a constraint for galaxy formation and evolution models.Comment: 2 pages, 2 figures, to appear in Proc. IAU Symp. 220, "Dark Matter in Galaxies" eds. S. Ryder, D.J. Pisano, M. Walker, & K. Freeman (San Francisco: ASP

SKIRT: the design of a suite of input models for Monte Carlo radiative transfer simulations

The Monte Carlo method is the most popular technique to perform radiative transfer simulations in a general 3D geometry. The algorithms behind and acceleration techniques for Monte Carlo radiative transfer are discussed extensively in the literature, and many different Monte Carlo codes are publicly available. On the contrary, the design of a suite of components that can be used for the distribution of sources and sinks in radiative transfer codes has received very little attention. The availability of such models, with different degrees of complexity, has many benefits. For example, they can serve as toy models to test new physical ingredients, or as parameterised models for inverse radiative transfer fitting. For 3D Monte Carlo codes, this requires algorithms to efficiently generate random positions from 3D density distributions. We describe the design of a flexible suite of components for the Monte Carlo radiative transfer code SKIRT. The design is based on a combination of basic building blocks (which can be either analytical toy models or numerical models defined on grids or a set of particles) and the extensive use of decorators that combine and alter these building blocks to more complex structures. For a number of decorators, e.g. those that add spiral structure or clumpiness, we provide a detailed description of the algorithms that can be used to generate random positions. Advantages of this decorator-based design include code transparency, the avoidance of code duplication, and an increase in code maintainability. Moreover, since decorators can be chained without problems, very complex models can easily be constructed out of simple building blocks. Finally, based on a number of test simulations, we demonstrate that our design using customised random position generators is superior to a simpler design based on a generic black-box random position generator.Comment: 15 pages, 4 figures, accepted for publication in Astronomy and Computin

The nature of the UV halo around the spiral galaxy NGC 3628

Thanks to deep UV observations with GALEX and Swift, diffuse UV haloes have recently been discovered around galaxies. Based on UV-optical colours, it has been advocated that the UV haloes around spiral galaxies are due to UV radiation emitted from the disc and scattered off dust grains at high latitudes. Detailed UV radiative transfer models that take into account scattering and absorption can explain the morphology of the UV haloes, and they require the presence of an additional thick dust disc next the to traditional thin disc for half of the galaxies in their sample. We test whether such an additional thick dust disc agrees with the observed infrared emission in NGC 3628, an edge-on galaxy with a clear signature of a thick dust disc. We extend the far-ultraviolet radiative transfer models to full-scale panchromatic models. Our model, which contains no fine-tuning, can almost perfectly reproduce the observed spectral energy distribution from UV to mm wavelengths. These results corroborate the interpretation of the extended UV emission in NGC 3628 as scattering off dust grains, and hence of the presence of a substantial amount of diffuse extra-planar dust. A significant caveat, however, is the geometrical simplicity and non-uniqueness of our model: other models with a different geometrical setting could lead to a similar spectral energy distribution. More detailed radiative transfer simulations that compare the model results to images from UV to submm wavelengths are a way to break this degeneracy, as are UV polarisation measurements.Comment: 6 pages, 2 figures, accepted for publication in Astronomy & Astrophysic

Black hole mass measurements using ionized gas discs: systematic dust effects

Using detailed Monte Carlo radiative transfer simulations in realistic models for galactic nuclei, we investigate the influence of interstellar dust in ionized gas discs on the rotation curves and the resulting black hole mass measurements. We find that absorption and scattering by interstellar dust leaves the shape of the rotation curves basically unaltered, but slightly decreases the central slope of the rotation curves. As a result, the "observed" black hole masses are systematically underestimated by some 10 to 20% for realistic optical depths. We therefore argue that the systematic effect of dust attenuation should be taken into account when estimating SMBH masses using ionized gas kinematics.Comment: 4 pages, 2 figures, to appear in "Observational Evidence for Black Holes in the Universe", AIP Conference Proceeding

Optical depth in polarised Monte Carlo radiative transfer

Context: The Monte Carlo method is the most widely used method to solve radiative transfer problems in astronomy, especially in a fully general 3D geometry. A crucial concept in any Monte Carlo radiative transfer code is the random generation of the next interaction location. In polarised Monte Carlo radiative transfer with aligned non-spherical grains, the nature of dichroism complicates the concept of optical depth. Aims: We investigate, in detail, the relation between optical depth and the optical properties and density of the attenuating medium in polarised Monte Carlo radiative transfer codes that take dichroic extinction into account. Methods: Based on solutions for the radiative transfer equation, we discuss the optical depth scale in polarised radiative transfer with spheroidal grains. We compare the dichroic optical depth to the extinction and total optical depth scale. Results: In a dichroic medium, the optical depth is not equal to the usual extinction optical depth, nor to the total optical depth. For representative values of the optical properties of dust grains, the dichroic optical depth can differ from the extinction or total optical depth by several tens of percent. A closed expression for the dichroic optical depth cannot be given, but it can be derived efficiently through an algorithm that is based on the analytical result corresponding to elongated grains with a uniform grain alignment. Conclusions: Optical depth is more complex in dichroic media than in systems without dichroic attenuation, and this complexity needs to be considered when generating random free path lengths in Monte Carlo radiative transfer simulations. There is no benefit in using approximations instead of the dichroic optical depth