Approximate third integrals for axisymmetric potentials using local Stäckel fits

We use a set of Stäckel potentials to obtain a local approximation for an effective third integral in axisymmetric systems. We present a study on the feasibility and effectiveness of this approach. We apply it to three trial potentials of various flattenings, corresponding to nearly ellipsoidal, discy and boxy density isophotes. In all three cases, a good fit to the potential requires only a small set of Stäckel potentials, and the associated Stäckel third integral provides a very satisfactory, yet analytically simple, approximation to the trial potentials effective third integral.     

Integrable models of galactic discs with double nuclei

We introduce a new class of 2D mass models, whose potentials are of Stäckel form in elliptic coordinates. Our model galaxies have two separate strong cusps that form double nuclei. The potential and surface density distributions are locally axisymmetric near the nuclei and become highly non-axisymmetric outside the nucleus. The surface density diverges toward the cuspy nuclei with the law     Our model is sustained by four general types of regular orbits: butterfly , nucleophilic banana , horseshoe and aligned loop orbits. Horseshoes and nucleophilic bananas support the existence of cuspy regions. Butterflies and aligned loops control the non-axisymmetric shape of outer regions. Without any need for central black holes, our distributed mass models resemble the nuclei of M31 and NGC 4486B. It is also shown that the self-gravity of the stellar disc can prevent the double nucleus to collapse.     

Exact solutions for the gravitational potential of a family of heterogeneous spheroids

An infinite family of heterogeneous spheroids has been found for which exact or closed-form solutions for the Newtonian gravitational potential can be given. The family includes both axisymmetric spheroids and spheroids where the matter density varies harmonically with the azimuthal angle. For the axisymmetric family of spheroids, which have no azimuthal dependence of the density, the potential external to the spheroid is of the same form as the potential exterior to a spheroidal homoeoid. It is therefore constant on the surface of the spheroid and on all external spheroidal surfaces confocal with it. The potential is also constant on all internal confocal spheroidal surfaces, with the value on each confocal surface dependent on the density distribution chosen. However, the density is not constant on either concentric or confocal spheroids. These solutions can be considered to be generalizations of analogous spherical solutions given in a companion paper by Conway. For the classical solutions for homogeneous spheroids, the surface is not equipotential, and these are not included within the new family of solutions, except in the special case of a homogeneous sphere.     

Numerical simulations of the formation and chemical evolution of galaxies

We present the results of a set of three-dimensional SPH-Treecode simulations which model the formation and early evolution of disc galaxies, including the generation and return of heavy elements to the interstellar medium by star formation. Starting from simple initial conditions which are given by a uniform density sphere of gas which is embedded in a dark matter halo and in solid-body rotation, we are able to form realistic disc galaxies, and find that an exponential gas disc is quickly formed. Star formation within this exponential disc naturally leads to the formation of abundance gradients which are in broad agreement with those observed, although they are slightly shallower than some observations.
We investigate the systematic effects of variation of mass, rotation and star formation parameters on the abundance gradients. We find that the abundance gradients are most sensitive to changes in the star formation parameters or rotation. Including a critical-density cut-off in the star formation law causes abundance gradients to be steepened.
Analysis of gas flows within the models shows radial flows which are a function of angle of azimuth around the galaxies, with alternating inward and outward flows. This motion is linked to the presence of a bar, whose strength is related to the amount of star formation in the models, and there is a gentle drift of mass inwards. The shallow abundance gradients may be linked to these radial flows.     

Alignment and precession of a black hole with a warped accretion disc

We consider the shape of an accretion disc whose outer regions are misaligned with the spin axis of a central black hole and calculate the steady state form of the warped disc in the case where the viscosity and surface densities are power laws in the distance from the central black hole. We discuss the shape of the resulting disc in both the frame of the black hole and that of the outer disc. We note that some parts of the disc and also any companion star maybe shadowed from the central regions by the warp. We compute the torque on the black hole caused by the Lense–Thirring precession, and hence compute the alignment and precession time-scales. We generalize the case with viscosity and hence surface density independent of radius to more realistic density distributions for which the surface density is a decreasing function of radius. We find that the alignment time-scale does not change greatly but the precession time-scale is more sensitive. We also determine the effect on this time-scale if we truncate the disc. For a given truncation radius, the time-scales are less affected for more sharply falling density distributions.     

Wave excitation in three-dimensional discs by external potential

We study the excitation of density and bending waves and the associated angular momentum transfer in gaseous discs with finite thickness by a rotating external potential. The disc is assumed to be isothermal in the vertical direction and has no self-gravity. The disc perturbations are decomposed into different modes, each characterized by the azimuthal index m and the vertical index n , which specifies the nodal number of the density perturbation along the disc normal direction. The   n = 0  modes correspond to the two-dimensional density waves previously studied by Goldreich & Tremaine and others. In a three-dimensional disc, waves can be excited at both Lindblad resonances (LRs; for modes with   n = 0, 1, 2, …  ) and vertical resonances (VRs; for the   n ≥ 1  modes only). The torque on the disc is positive for waves excited at outer Lindblad/vertical resonances and negative at inner Lindblad/vertical resonances. While the   n = 0  modes are evanescent around corotation, the   n ≥ 1  modes can propagate into the corotation region where they are damped and deposit their angular momenta. We have derived analytical expressions for the amplitudes of different wave modes excited at LRs and/or VRs and the resulting torques on the disc. It is found that for   n ≥ 1  , angular momentum transfer through VRs is much more efficient than LRs. This implies that in some situations (e.g. a circumstellar disc perturbed by a planet in an inclined orbit), VRs may be an important channel of angular momentum transfer between the disc and the external potential. We have also derived new formulae for the angular momentum deposition at corotation and studied wave excitations at disc boundaries.     

Exact potential–density pairs for flattened dark haloes

Cosmological simulations suggest that dark matter haloes are not spherical, but typically moderately to strongly triaxial systems. We investigate methods to convert spherical potential–density pairs into axisymmetric ones, in which the basic characteristics of the density profile (such as the slope at small and large radii) are retained. We achieve this goal by replacing the spherical radius r by an oblate radius m in the expression of the gravitational potential  Φ( r )  .
We extend and formalize the approach pioneered by Miyamoto & Nagai to be applicable to arbitrary potential–density pairs. Unfortunately, an asymptotic study demonstrates that, at large radii, such models always show a   R ⁻³  disc superposed on a smooth roughly spherical density distribution. As a result, this recipe cannot be used to construct simple flattened potential–density pairs for dynamical systems such as dark matter haloes. Therefore, we apply a modification of our original recipe that cures the problem of the discy behaviour. An asymptotic analysis now shows that the density distribution has the desired asymptotic behaviour at large radii (if the density falls less rapidly than   r ⁻⁴  ). We also show that the flattening procedure does not alter the shape of the density distribution at small radii: while the inner density contours are flattened, the slope of the density profile is unaltered.
We apply this recipe to construct a set of flattened dark matter haloes based on the realistic spherical halo models by Dehnen & McLaughlin. This example illustrates that the method works fine for modest flattening values, whereas stronger flattening values lead to peanut-shaped density distributions.     

Power spectra from ‘spotted’ accretion discs

We are carrying out a project to calculate power spectra of variability, assuming a model of a ‘spotted’ accretion disc near a black hole. We consider relativistic effects that change photon energy and produce light‐bending and time‐delays acting on the X‐ray signal received by an observer. We assume that the life‐time and the intrinsic emissivity of individual.aring events are described in terms of a simple stochastic process. This allows us to give approximate analytical formulae and compare them with numerical computations. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Lyman α radiative transfer in a multiphase medium

Hydrogen Lyman α (Lyα) is our primary emission-line window into high-redshift galaxies. Despite an extensive literature, Lyα radiative transfer in the most realistic case of a dusty, multiphase medium has received surprisingly little detailed theoretical attention. We investigate Lyα resonant scattering through an ensemble of dusty, moving, optically thick gas clumps. We treat each clump as a scattering particle and use Monte Carlo simulations of surface scattering to quantify continuum and Lyα surface scattering angles, absorption probabilities, and frequency redistribution, as a function of the gas dust content. This atomistic approach speeds up the simulations by many orders of magnitude, making possible calculations which are otherwise intractable. Our fitting formulae can be readily adapted for fast radiative transfer in numerical simulations. With these surface scattering results, we develop an analytic framework for estimating escape fractions and line widths as a function of gas geometry, motion, and dust content. Our simple analytic model shows good agreement with full Monte Carlo simulations. We show that the key geometric parameter is the average number of surface scatters for escape in the absence of absorption,     , and we provide fitting formulae for several geometries of astrophysical interest. We consider the following two interesting applications. (i) Equivalent widths ( EWs ). Lyα can preferentially escape from a dusty multiphase interstellar medium if most of the dust lies in cold neutral clouds, which Lyα photons cannot penetrate. This might explain the anomalously high EWs sometimes seen in high-redshift/submillimetre sources. (ii) Multiphase galactic outflows . We show the characteristic profile is asymmetric with a broad red tail, and relate the profile features to the outflow speed and gas geometry. Many future applications are envisaged.     

Modelling star formation and feedback in simulations of interacting galaxies

We discuss a heuristic model to implement star formation and feedback in hydrodynamical simulations of galaxy formation and evolution. In this model, gas is allowed to cool radiatively and to form stars at a rate given by a simple Schmidt-type law. We assume that supernova feedback results in turbulent motions of gas below resolved scales, a process that can pressurize the diffuse gaseous medium effectively, even if it lacks substantial thermal support. Ignoring the complicated detailed physics of the feedback processes, we try to describe their net effect on the interstellar medium with a fiducial second reservoir of internal energy, which accounts for the kinetic energy content of the gas on unresolved scales. Applying the model to three-dimensional, fully self-consistent models of isolated disc galaxies, we show that the resulting feedback loop can be modelled with smoothed particle hydrodynamics such that converged results can be reached with moderate numerical resolution. With an appropriate choice of the free parameters, Kennicutt's phenomenological star formation law can be reproduced over many orders of magnitude in gas surface density. We also apply the model to mergers of equal-mass disc galaxies, typically resulting in strong nuclear starbursts. Confirming previous findings, the presence of a bulge can delay the onset of the starburst from the first encounter of the galaxies until their final coalescence. The final density profiles of the merger remnants are consistent with de Vaucouleurs profiles, except for the innermost region, where the newly created stars give rise to a luminous core with stellar densities that may be in excess of those observed in the cores of most elliptical galaxies. By comparing the isophotal shapes of collisionless and dissipative merger simulations we show that dissipation leads to isophotes that are more discy than those of corresponding collisionless simulations.     

The Monitor project: JW 380 – a 0.26-, 0.15-M⊙, pre-main-sequence eclipsing binary in the Orion nebula cluster

We present a general recipe for constructing N -body realizations of galaxies comprising near spherical and disc components. First, an exact spherical distribution function for the spheroids (halo and bulge) is determined, such that it is in equilibrium with the gravitational monopole of the disc components. Second, an N -body realization of this model is adapted to the full disc potential by growing the latter adiabatically from its monopole. Finally, the disc is sampled with particles drawn from an appropriate distribution function, avoiding local-Maxwellian approximations. We performed test simulations and find that the halo and bulge radial density profile very closely match their target model, while they become slightly oblate due to the added disc gravity. Our findings suggest that vertical thickening of the initially thin disc is caused predominantly by spiral and bar instabilities, which also result in a radial re-distribution of matter, rather than scattering off interloping massive halo particles.     

Exact density–potential pairs from complex-shifted axisymmetric systems

In a previous paper, the complex-shift method has been applied to self-gravitating spherical systems, producing new analytical axisymmetric density–potential pairs. We now extend the treatment to the Miyamoto–Nagai disc and the Binney logarithmic halo, and we study the resulting axisymmetric and triaxial analytical density–potential pairs; we also show how to obtain the surface density of shifted systems from the complex shift of the surface density of the parent model. In particular, the systems obtained from Miyamoto–Nagai discs can be used to describe disc galaxies with a peanut-shaped bulge or with a central triaxial bar, depending on the direction of the shift vector. By using a constructive method that can be applied to generic axisymmetric systems, we finally show that the Miyamoto–Nagai and the Satoh discs, and the Binney logarithmic halo cannot be obtained from the complex shift of any spherical parent distribution. As a by-product of this study, we also found two new generating functions in closed form for even and odd Legendre polynomials, respectively.     

The effect of Poisson noise on SPH calculations

A simple algorithm is presented which generates a lattice-like, regular spacing of smoothed particle hydrodynamics (SPH) particles in discs, with any prescribed density gradient. Using this for comparison we demonstrate the effect of Poisson noise on SPH estimates of density, pressure and viscous forces when calculated using randomly distributed particles. The standard deviation of density and pressure is typically increased by greater than an order of magnitude. In a disc with a Keplerian velocity profile, the effectiveness of the Balsara switch in reducing the shear component of SPH artificial viscosity is greatly enhanced when the particles are properly spaced, reducing the magnitude of viscosity by two orders of magnitude. Noise problems are exacerbated, not removed, by increasing the numbers of SPH particles, if the number of neighbours used is kept constant. However, comparison of the evolution of a disc created using perfectly spaced particles and a disc with identical parameters but randomly placed particles, reveals very similar results. Although there are subtle differences in the evolution, and the smooth disc takes longer to begin developing structure, in both cases the identical number of objects is created by gravitational collapse. SPH disc simulations do not depend on initial density disturbances to evolve objects by gravitational collapse, which gives added credence to the validity of the results. It also appears that complicated disc settling procedures are unnecessary.     

Self-consistent response of a galactic disc to vertical perturbations

We study the self-consistent, linear response of a galactic disc to vertical perturbations, as induced, say, by a tidal interaction. We calculate the self-gravitational potential corresponding to a non-axisymmetric, self-consistent density response of the disc using the Green's function method. The response potential is shown to oppose the perturbation potential because the self-gravity of the disc resists the imposed potential, and this resistance is stronger in the inner parts of a galactic disc. For the   m = 1  azimuthal wavenumber, the disc response opposes the imposed perturbation up to a radius that spans a range of 4–6 disc scalelengths, so that the disc shows a net warp only beyond this region. This physically explains the well known but so far unexplained observation that warps typically set in beyond this range of radii. We show that the inclusion of a dark matter halo in the calculation only marginally changes (by ∼10 per cent) the radius for the onset of warps. For perturbations with higher azimuthal wavenumbers, the net signature of the vertical perturbations can only be seen at larger radii – for example, beyond 7 exponential disc scalelengths for   m = 10  . Also, for the high- m cases, the magnitude of the negative disc response due to the disc self-gravity is much smaller. This is shown to result in corrugations of the mid-plane density, which explains the puzzling scalloping with   m = 10  detected in H  i in the outermost regions ∼30 kpc in the Galaxy.     

Measuring angular diameters of extended sources

When measuring diameters of partially resolved sources like planetary nebulae, H  ii regions or galaxies, often a technique called Gaussian deconvolution is used. This technique yields a Gaussian diameter, which subsequently has to be multiplied by a conversion factor to obtain the true angular diameter of the source. This conversion factor is a function of the FWHM of the beam or point spread function, and also depends on the intrinsic surface brightness distribution of the source.
In this paper, conversion factors are presented for a number of simple geometries: a circular constant surface brightness disc and a spherical constant emissivity shell, using a range of values for the inner radius. Also, more realistic geometries are studied, based on a spherically symmetric photoionization model of a planetary nebula. This enables a study of optical depth effects, a comparison between images in various emission lines, and the use of power-law density distributions. It is found that the conversion factor depends quite critically on the intrinsic surface brightness distribution, which is usually unknown. The uncertainty is particularly large if extended regions of low surface brightness are present in the nebula. In such cases the use of Gaussian or second-moment deconvolution is not recommended.
As an alternative, a new algorithm is presented which allows the determination of the intrinsic FWHM of the source using only the observed surface brightness distribution and the FWHM of the beam. Hence no assumptions concerning the intrinsic surface brightness distribution are needed. Tests show that this implicit deconvolution method works well in realistic conditions, even when the signal-to-noise ratio is low, provided that the beamsize is less than roughly 2/3 of the observed FWHM and the beam profile can be approximated by a Gaussian. A code implementing this algorithm is available.     

Effects of ram pressure on the gas distribution and star formation in the Large Magellanic Cloud

We use high-resolution N -body/smoothed particle hydrodynamics (SPH) simulations to study the hydrodynamical interaction between the Large Magellanic Cloud (LMC) and the hot halo of the Milky Way. We investigate whether ram pressure acting on the satellite's interstellar medium can explain the peculiarities observed in the H  i distribution and the location of the recent star formation activity.
Due to the present nearly edge-on orientation of the disc with respect to the orbital motion, compression at the leading edge can explain the high density region observed in H  i at the south-east border. In the case of a face-on disc (according to Mastropietro the LMC was moving almost face-on before the last perigalactic passage), ram pressure directed perpendicular to the disc produces a clumpy structure characterized by voids and high density filaments that resemble those observed by the Parkes H  i survey. As a consequence of the very recent edge-on motion, the Hα emission is mainly concentrated on the eastern side where 30 Doradus and most of the supergiant shells are located, although some Hα complexes form a patchy distribution on the entire disc. In this scenario, only the youngest stellar complexes show a progression in age along the leading border of the disc.     

Damped Lyα systems at high redshift and models of protogalactic discs

We employ observationally determined intrinsic velocity widths and column densities of damped Lyman alpha (Lyα) systems at high redshift to investigate the distribution of baryons in protogalaxies within the context of a standard cold dark matter (CDM) model. We proceed under the assumption that damped Lyα systems represent a population of cold, rotationally supported, protogalactic discs, and that the abundance of dark matter haloes is well approximated by a CDM model with critical density and vanishing cosmological constant. Using conditional cross-sections to observe a damped system with a given velocity width and column density, we compare observationally inferred velocity width and column density distributions to the corresponding theoretically determined distributions for a variety of disc parameters and CDM normalizations. In general, we find that the observations cannot be reproduced by the models for most disc parameters and CDM normalizations. Whereas the column density distribution favours small discs with large neutral gas fraction, the velocity width distribution favours large and thick discs with small neutral gas fraction. The possible resolutions of this problem in the context of this CDM model may be (1) an increased contribution of rapidly rotating discs within massive dark matter haloes to damped Lyα absorption, or (2) the abandoning of simple disc models within this CDM model for damped Lyα systems at high redshift. Here the first possibility may be achieved by supposing that damped Lyα system formation occurs only in haloes with fairly large circular velocities, and the second possibility may result from a large contribution of mergers and double discs to damped Lyα absorption at high redshift.     

Multimass spherical structure models for N-body simulations

We present a simple and efficient method to set up spherical structure models for N -body simulations with a multimass technique. This technique reduces by a substantial factor the computer run time needed in order to resolve a given scale as compared to single-mass models. It therefore allows to resolve smaller scales in N -body simulations for a given computer run time. Here, we present several models with an effective resolution of up to  1.68 × 10⁹  particles within their virial radius which are stable over cosmologically relevant time-scales. As an application, we confirm the theoretical prediction by Dehnen that in mergers of collisionless structures like dark matter haloes always the cusp of the steepest progenitor is preserved. We model each merger progenitor with an effective number of particles of approximately 10⁸ particles. We also find that in a core–core merger the central density approximately doubles whereas in the cusp–cusp case the central density only increases by approximately 50 per cent. This may suggest that the central regions of flat structures are better protected and get less energy input through the merger process.