BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Walter Dehnen

Walter Dehnen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-phastro-ph.IMastro-ph.EPphysics.comp-phastro-ph.HEphysics.flu-dyn

Catalog footprint

What is connected

25works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2026arXiv

Converging on the eccentricity of massive black hole binaries in galactic mergers

Binaries of supermassive black holes (MBHBs) represent the primary sources of the gravitational wave background (GWB) detectable by Pulsar Timing Arrays (PTAs). The eccentricity with which binaries form in galactic mergers is the key parameter determining their evolutionary timescale from pairing to coalescence. However, accurately determining the binary eccentricity at formation is difficult in simulations due to stochastic effects. We present a numerical study of the formation and evolution of MBHBs that are potential PTA sources. We simulate mergers of equal-mass galaxies on different initial orbits and follow the dynamics of the MBHBs through the hardening phase. We find that low-resolution simulations are affected by stochasticity due to torques from the stellar distribution acting at pericentre passages. The dispersion in binary eccentricity decreases with increasing central resolution, as expected for a Poisson process. We provide a fitting formula for the resolution requirement of an N-body simulation of MBHB formation and evolution as a function of the initial eccentricity of the merger, e0, and the required accuracy in the binary eccentricity, eb. We find that binaries experience a torque at first pericentre that is approximately independent of initial eccentricity, producing a general trend in which the binary eccentricity decreases above sufficiently large initial orbital eccentricities. While this behaviour is generic, the precise cross-over eccentricity (e0~0.97 in our models) and the sharpness of the drop-off depend on the galaxy initial conditions. We provide a fitting formula for eb(e0) that can be used in semi-analytical models to determine the merger timescales of MBHBs as well as the amplitude and slope of the GWB.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Secular Dynamics around a Supermassive Black Hole via Multipole Expansion

In galactic nuclei, the gravitational potential is dominated by the central supermassive black hole, so stars follow quasi-Keplerian orbits. These orbits are distorted by gravitational forces from other stars, leading to long-term orbital relaxation. The direct numerical study of these processes is challenging because the fast orbital motion imposed by the central black hole requires very small timesteps. An alternative approach, pioneered by Gauss, is to use the secular approximation of smearing out the $N$ stars over their Keplerian orbits, using $K$ nodes along each orbit. In this study we propose three novel improvements to this method. First, we re-formulate the discretisation of the rates of change of the variables describing the orbital states to ensure that all conservation laws are exactly satisfied. Second, we replace the pairwise sum over nodes by a multipole expansion up to order $\ell_{\mathrm{max}}$, reducing the overall computational costs from $O(N^2K^2)$ to $O(NK\ell_{\mathrm{max}}^2)$. Finally, we show that the averaged dynamical system is equivalent to $2N$ interacting unit spin vectors and provide two time integrators: a second-order symplectic scheme and a fourth-order Lie-group Runge-Kutta method, both of which are straightforward to generalize to higher order. These new simulations recover the diffusion coefficients of stellar eccentricities obtained through analytical calculations of the secular dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Capture of interstellar objects I: the capture cross-section

We study the capture of interstellar objects (ISOs) by a planet-star binary with mass ratio $q\ll1$, semi-major axis $a_p$, orbital speed $v_c$, and eccentricity $e_p$. Very close (slingshot) and wide encounters with the planet are amenable to analytic treatment, while numerically obtained capture cross-sections $σ$ closely follow the analytical results even in the intermediate regime. Wide interactions can only generate energy changes $ΔE\lesssim qv_c^2$, when $σ\propto v_\infty^{-2} |\lnΔE|^{2/3}$ (with $v_\infty$ the ISO's incoming speed far away from the binary), which is slightly enhanced for $e_p>0$. Energy changes $ΔE\gtrsim qv_c^2$, on the other hand, require close interactions when $σ\propto (v_\infty ΔE)^{-2}$ hardly depending on $e_p$. Finally, at $ΔE\gtrsim v_c^2$, the cross-section drops to zero, depending on the planet's radius $R_p$ through the Safronov number $Θ=qa_p/R_p$. We also derive the cross-sections for collisions of ISOs with planets or moons.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Capture of interstellar objects II: by the Solar system

Capture of interstellar objects (ISOs) into the Solar system is dominated by ISOs with asymptotic incoming speeds $v_\infty<4\,$km\,s$^{-1}$. The capture rate is proportional to the ISO phase-space density in the Solar vicinity and does not vary along the Sun's Galactic orbit, i.e.\ is not enhanced during a passage through a cloud of ISOs (in contrast to previous suggestions). Most bound orbits crossing those of Jupiter and Saturn are fully mixed with unbound phase space, implying that they hold the same ISO phase-space density. Assuming an interstellar number density $n_{iso}\sim0.1\,$au$^{-3}$, we estimate that in 1000 years the planets capture $\sim2$ ISOs (while $\sim17$ fall into the Sun), resulting in a population of $\sim8$ captured ISOs within 5\,au of the Sun at any time, less than the number of visiting ISOs passing through the same volume on hyperbolic orbits. In terms of phase-space volume, capture onto and ejection from the Solar system are equal, such that on average ISOs will not remain captive at $a\lesssim2000\,$au for extensive periods.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Formation of the largest galactic cores through binary scouring and gravitational wave recoil

Massive elliptical galaxies are typically observed to have central cores in their projected radial light profiles. Such cores have long been thought to form through `binary scouring' as supermassive black holes (SMBHs), brought in through mergers, form a hard binary and eject stars from the galactic centre. However, the most massive cores, like the ~3kpc core in A2261-BCG, remain challenging to explain in this way. In this paper, we run a suite of dry galaxy merger simulations to explore three different scenarios for central core formation in massive elliptical galaxies: `binary scouring', `tidal deposition' and `gravitational wave (GW) induced recoil'. Using the Griffin code, we self-consistently model the stars, dark matter and SMBHs in our merging galaxies, following the SMBH dynamics through to the formation of a hard binary. We find that we can only explain the large surface brightness core of A2261-BCG with a combination of a major merger that produces a small ~1kpc core through binary scouring, followed by the subsequent GW recoil of its SMBH that acts to grow the core size. We show that this same model can also explain the bright `knots' observed in the core region of A2261-BCG. Key predictions of this scenario are an offset SMBH surrounded by a compact cluster of bound stars and a non-divergent central density profile. We show that the bright `knots' observed in the core region of A2261-BCG are best explained as stalled perturbers resulting from minor mergers, though the brightest may also represent ejected SMBHs surrounded by a stellar cloak of bound stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Capture of interstellar objects: a source of long-period comets

We simulate the passage through the Sun-Jupiter system of interstellar objects (ISOs) similar to 1I/`Oumuamua or 2I/Borisov. Capture of such objects is rare and overwhelmingly from low incoming speeds onto orbits akin to those of known long-period comets. This suggests that some of these comets could be of extra-solar origin, in particular inactive ones. Assuming ISOs follow the local stellar velocity distribution, we infer a volume capture rate of $0.051\,\mathrm{au}^3 \mathrm{yr}^{-1}$. Current estimates for orbital lifetimes and space densities then imply steady-state captured populations of $\sim10^2$ comets and $\sim10^5$ `Oumuamua-like rocks, of which 0.033% are within 6au at any time.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Defeating stochasticity: coalescence timescales of massive black holes in galaxy mergers

The coalescence of massive black hole binaries (BHBs) in galactic mergers is the primary source of gravitational waves (GWs) at low frequencies. Current estimates of GW detection rates for the Laser Interferometer Space Antenna and the Pulsar Timing Array vary by three orders of magnitude. To understand this variation, we simulate the merger of equal-mass, eccentric, galaxy pairs with central massive black holes and shallow inner density cusps. We model the formation and hardening of a central BHB using the Fast Multiple Method as a force solver, which features a $O(N)$ scaling with the number $N$ of particles and obtains results equivalent to direct-summation simulations. At $N \sim 5\times 10^5$, typical for contemporary studies, the eccentricity of the BHBs can vary significantly for different random realisations of the same initial condition, resulting in a substantial variation of the merger timescale. This scatter owes to the stochasticity of stellar encounters with the BHB and decreases with increasing $N$. We estimate that $N \sim 10^7$ within the stellar half-light radius suffices to reduce the scatter in the merger timescale to $\sim 10$\%. Our results suggest that at least some of the uncertainty in low-frequency GW rates owes to insufficient numerical resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-parametric Triaxial Deprojection of Elliptical Galaxies

We present a grid-based non-parametric approach to obtain a triaxial three-dimensional luminosity density from its surface brightness distribution. Triaxial deprojection is highly degenerate and our approach illustrates the resulting difficulties. Fortunately, for massive elliptical galaxies, many deprojections for a particular line of sight can be discarded, because their projection along other lines of sight does not resemble elliptical galaxies. The near-elliptical isophotes of these objects imply near ellipsoidal intrinsic shapes. In fact, deprojection is unique for densities distributed on ellipsoidal shells. The constrained non-parametric deprojection method we present here relaxes this constraint and assumes that the contours of the luminosity density are boxy/discy ellipsoids with radially varying axis ratios. With this approach we are able to reconstruct the intrinsic triaxial densities of our test models, including one drawn from an $N$-body simulation. The method also allows to compare the relative likelihood of deprojections at different viewing angles. We show that the viewing orientations of individual galaxies with nearly ellipsoidal isophotal shapes can be constrained from photometric data alone.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Collisionless loss-cone refilling: there is no final parsec problem

Coalescing massive black hole binaries, formed during galaxy mergers, are expected to be a primary source of low frequency gravitational waves. Yet in isolated gas-free spherical stellar systems, the hardening of the binary stalls at parsec-scale separations owing to the inefficiency of relaxation-driven loss-cone refilling. Repopulation via collisionless orbit diffusion in triaxial systems is more efficient, but published simulation results are contradictory. While sustained hardening has been reported in simulations of galaxy mergers with $N \sim 10^6$ stars and in early simulations of rotating models, in isolated non-rotating triaxial models the hardening rate continues to fall with increasing N, a signature of spurious two-body relaxation. We present a novel approach for studying loss cone repopulation in galactic nuclei. Since loss cone repopulation in triaxial systems owes to orbit diffusion, it is a purely collisionless phenomenon and can be studied with an approximated force calculation technique, provided the force errors are well behaved and sufficiently small. We achieve this using an accurate fast multipole method and define a proxy for the hardening rate that depends only on stellar angular momenta. We find that the loss cone is efficiently replenished even in very mildly triaxial models (with axis ratios 1 : 0.9 : 0.8). Such triaxiality is unavoidable following galactic mergers and can drive binaries into the gravitational wave regime. We conclude that there is no final parsec problem.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Indirect Dark Matter Detection for Flattened Dwarf Galaxies

Gamma-ray experiments seeking to detect evidence of dark matter annihilation in dwarf spheroidal galaxies require knowledge of the distribution of dark matter within these systems. We analyze the effects of flattening on the annihilation (J) and decay (D) factors of dwarf spheroidal galaxies with both analytic and numerical methods. Flattening has two consequences: first, there is a geometric effect as the squeezing (or stretching) of the dark matter distribution enhances (or diminishes) the J-factor; second, the line of sight velocity dispersion of stars must hold up the flattened baryonic component in the flattened dark matter halo. We provide analytic formulae and a simple numerical approach to estimate the correction to the J- and D-factors required over simple spherical modeling. The formulae are validated with a series of equilibrium models of flattened stellar distributions embedded in flattened dark-matter distributions. We compute corrections to the J- and D-factors for the Milky Way dwarf spheroidal galaxies under the assumption that they are all prolate or all oblate and find that the hierarchy of J-factors for the dwarf spheroidals is slightly altered (typical correction factors for an ellipticity of $0.4$ are $0.75$ for the oblate case and $1.6$ for the prolate case). We demonstrate that spherical estimates of the D-factors are very insensitive to the flattening and introduce uncertainties significantly less than the uncertainties in the D-factors from the other observables for all the dwarf spheroidals (for example, ${}^{+10\mathrm{\,percent}}_{-3\mathrm{\,percent}}$ for a typical ellipticity of $0.4$). We conclude by investigating the spread in correction factors produced by triaxial figures and provide uncertainties in the J-factors for the dwarf spheroidals using different physically-motivated assumptions for their intrinsic shape and axis alignments. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Misaligned gas discs around eccentric black-hole binaries and implications for the final-parsec problem

We investigate the evolution of low mass (Md /Mb = 0.005) misaligned gaseous discs around eccentric supermassive black hole (SMBH) binaries. These are expected to form from randomly oriented accretion events onto a SMBH binary formed in a galaxy merger. When expanding the interaction terms between the binary and a circular ring to quadrupole order and averaging over the binary orbit, we expect four non-precessing disc orientations: aligned or counter-aligned with the binary, or polar orbits around the binary eccentricity vector with either sense of rotation. All other orientations precess around either of these, with the polar precession dominating for high eccentricity. These expectations are borne out by smoothed particle hydrodynamics simulations of initially misaligned viscous circumbinary discs, resulting in the formation of polar rings around highly eccentric binaries in contrast to the co-planar discs around circular binaries. Moreover, we observe disc tearing and violent interactions between differentially precessing rings in the disc significantly disrupting the disc structure and causing gas to fall onto the binary with little angular momentum. While accretion from a polar disc may not promote SMBH binary coalescence (solving the `final-parsec problem'), ejection of this infalling low-angular momentum material via gravitational slingshot is a possible mechanism to reduce the binary separation. Moreover, this process acts on dynamical rather than viscous time scales, and so is much faster.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A fast multipole method for stellar dynamics

The approximate computation of all gravitational forces between $N$ interacting particles via the fast multipole method (FMM) can be made as accurate as direct summation, but requires less than $\mathcal{O}(N)$ operations. FMM groups particles into spatially bounded cells and uses cell-cell interactions to approximate the force at any position within the sink cell by a Taylor expansion obtained from the multipole expansion of the source cell. By employing a novel estimate for the errors incurred in this process, I minimise the computational effort required for a given accuracy and obtain a well-behaved distribution of force errors. For relative force errors of $\sim10^{-7}$, the computational costs exhibit an empirical scaling of $\propto N^{0.87}$. My implementation (running on a 16 core node) out-performs a GPU-based direct summation with comparable force errors for $N\gtrsim10^5$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Black hole foraging: feedback drives feeding

We suggest a new picture of supermassive black hole (SMBH) growth in galaxy centers. Momentum-driven feedback from an accreting hole gives significant orbital energy but little angular momentum to the surrounding gas. Once central accretion drops, the feedback weakens and swept-up gas falls back towards the SMBH on near-parabolic orbits. These intersect near the black hole with partially opposed specific angular momenta, causing further infall and ultimately the formation of a small-scale accretion disk. The feeding rates into the disk typically exceed Eddington by factors of a few, growing the hole on the Salpeter timescale and stimulating further feedback. Natural consequences of this picture include (i) the formation and maintenance of a roughly toroidal distribution of obscuring matter near the hole; (ii) random orientations of successive accretion disk episodes; (iii) the possibility of rapid SMBH growth; (iv) tidal disruption of stars and close binaries formed from infalling gas, resulting in visible flares and ejection of hypervelocity stars; (v) super-solar abundances of the matter accreting on to the SMBH; and (vi) a lower central dark-matter density, and hence annihilation signal, than adiabatic SMBH growth implies. We also suggest a simple sub-grid recipe for implementing this process in numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational Collapse in One Dimension

We simulate the evolution of one-dimensional gravitating collisionless systems from non- equilibrium initial conditions, similar to the conditions that lead to the formation of dark- matter halos in three dimensions. As in the case of 3D halo formation we find that initially cold, nearly homogeneous particle distributions collapse to approach a final equilibrium state with a universal density profile. At small radii, this attractor exhibits a power-law behavior in density, ρ(x) \propto |x|^(-γ_crit), γ_crit \simeq 0.47, slightly but significantly shallower than the value γ = 1/2 suggested previously. This state develops from the initial conditions through a process of phase mixing and violent relaxation. This process preserves the energy ranks of particles. By warming the initial conditions, we illustrate a cross-over from this power-law final state to a final state containing a homogeneous core. We further show that inhomogeneous but cold power-law initial conditions, with initial exponent γ_i > γ_crit, do not evolve toward the attractor but reach a final state that retains their original power-law behavior in the interior of the profile, indicating a bifurcation in the final state as a function of the initial exponent. Our results rely on a high-fidelity event-driven simulation technique.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Improving convergence in smoothed particle hydrodynamics simulations without pairing instability

The numerical convergence of smoothed particle hydrodynamics (SPH) can be severely restricted by random force errors induced by particle disorder, especially in shear flows, which are ubiquitous in astrophysics. The increase in the number NH of neighbours when switching to more extended smoothing kernels at fixed resolution (using an appropriate definition for the SPH resolution scale) is insufficient to combat these errors. Consequently, trading resolution for better convergence is necessary, but for traditional smoothing kernels this option is limited by the pairing (or clumping) instability. Therefore, we investigate the suitability of the Wendland functions as smoothing kernels and compare them with the traditional B-splines. Linear stability analysis in three dimensions and test simulations demonstrate that the Wendland kernels avoid the pairing instability for all NH, despite having vanishing derivative at the origin (disproving traditional ideas about the origin of this instability; instead, we uncover a relation with the kernel Fourier transform and give an explanation in terms of the SPH density estimator). The Wendland kernels are computationally more convenient than the higher-order B-splines, allowing large NH and hence better numerical convergence (note that computational costs rise sub-linear with NH). Our analysis also shows that at low NH the quartic spline kernel with NH ~= 60 obtains much better convergence then the standard cubic spline.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The mass distribution of the Fornax dSph: constraints from its globular cluster distribution

Uniquely among the dwarf spheroidal (dSph) satellite galaxies of the Milky Way, Fornax hosts globular clusters. It remains a puzzle as to why dynamical friction has not yet dragged any of Fornax's five globular clusters to the centre, and also why there is no evidence that any similar star cluster has been in the past (for Fornax or any other dSph). We set up a suite of 2800 N-body simulations that sample the full range of globular-cluster orbits and mass models consistent with all existing observational constraints for Fornax. In agreement with previous work, we find that if Fornax has a large dark-matter core then its globular clusters remain close to their currently observed locations for long times. Furthermore, we find previously unreported behaviour for clusters that start inside the core region. These are pushed out of the core and gain orbital energy, a process we call 'dynamical buoyancy'. Thus a cored mass distribution in Fornax will naturally lead to a shell-like globular cluster distribution near the core radius, independent of the initial conditions. By contrast, CDM-type cusped mass distributions lead to the rapid infall of at least one cluster within Δt = 1-2Gyr, except when picking unlikely initial conditions for the cluster orbits (\sim 2% probability), and almost all clusters within Δt = 10Gyr. Alternatively, if Fornax has only a weakly cusped mass distribution, dynamical friction is much reduced. While over Δt = 10Gyr this still leads to the infall of 1-4 clusters from their present orbits, the infall of any cluster within Δt = 1-2Gyr is much less likely (with probability 0-70%, depending on Δt and the strength of the cusp). Such a solution to the timing problem requires that in the past the globular clusters were somewhat further from Fornax than today; they most likely did not form within Fornax, but were accreted.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modeling disc non-axisymmetries: multiple patterns, radial migration, and thick discs

Disc non-axisymmetrc components, such as spirals and central bars, are nowadays known to play an important role in shaping galactic discs. Here we use Tree-SPH N-body simulations to examine the effect of these perturbers on two aspects: the occurrence of multiple patterns in discs and the effects of radial migration on disc thickening. We find that, in addition to a central bar, multiple spiral patterns and lopsided modes develop in all models. Interaction among these asymmetric features results in a large scale stellar migration. However, we show that, despite the strong radial mixing, discs cannot be thickened sufficiently to match observed thick discs. We relate this to the adiabatic cooling as stars migrate radially outwards. We also find that the bulge contribution to a thick-disc component for an Sa-type galaxy at ~2.5 disc scale-lengths is less than 1% and zero in the case of a Milky Way-like, Sb-type. Our findings cast doubt on the plausibility of thick disc formation via stellar radial migration.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

N-body simulations of gravitational dynamics

We describe the astrophysical and numerical basis of N-body simulations, both of collisional stellar systems (dense star clusters and galactic centres) and collisionless stellar dynamics (galaxies and large-scale structure). We explain and discuss the state-of-the-art algorithms used for these quite different regimes, attempt to give a fair critique, and point out possible directions of future improvement and development. We briefly touch upon the history of N-body simulations and their most important results.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Weakening dark-matter cusps by clumpy baryonic infall

We consider the infall of a massive clump into a dark-matter halo as a simple and extreme model for the effect of baryonic physics (neglected in gravity-only simulations of large-scale structure formation) on the dark-matter. We find that such an infalling clump is extremely efficient in altering the structure of the halo and reducing its central density: a clump of 1% the mass of the halo can remove about twice its own mass from the inner halo and transform a cusp into a core or weaker cusp. If the clump is subsequently removed, mimicking a galactic wind, the central halo density is further reduced and the mass removed from the inner halo doubled. Lighter clumps are even more efficient: the ratio of removed mass to clump mass increases slightly towards smaller clump masses. This process is the more efficient the more radially anisotropic the initial dark-matter velocities. While such a clumpy infall may be somewhat unrealistic, it demonstrates that the baryons need to transfer only a small fraction of their initial energy to the dark matter via dynamical friction to explain the discrepancy between predicted dark-matter density profiles and those inferred from observations of dark-matter dominated galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Inviscid SPH

In smooth-particle hydrodynamics (SPH), artificial viscosity is necessary for the correct treatment of shocks, but often generates unwanted dissipation away from shocks. We present a novel method of controlling the amount of artificial viscosity, which uses the total time derivative of the velocity divergence as shock indicator and aims at completely eliminating viscosity away from shocks. We subject the new scheme to numerous tests and find that the method works at least as well as any previous technique in the strong-shock regime, but becomes virtually inviscid away from shocks, while still maintaining particle order. In particular sound waves or oscillations of gas spheres are hardly damped over many periods.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Tailoring triaxial N-body models via a novel made-to-measure method

The made-to-measure N-body method (Syer & Tremaine 1996) slowly adapts the particle weights of an N-body model, whilst integrating the trajectories in an assumed static potential, until some constraints are satisfied, such as optimal fits to observational data. I propose a novel technique for this adaption procedure, which overcomes several limitations and shortcomings of the original method. The capability of the new technique is demonstrated by generating realistic N-body equilibrium models for dark-matter haloes with prescribed density profile, triaxial shape, and slowly outwardly growing radial velocity anisotropy

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

The Distribution of Nearby Stars in Velocity Space Inferred from Hipparcos Data

(abridged) The velocity distribution f(v) of nearby stars is estimated, via a maximum- likelihood algorithm, from the positions and tangential velocities of a kinematically unbiased sample of 14369 stars observed by the HIPPARCOS satellite. f(v) shows rich structure in the radial and azimuthal motions, v_R and v_phi, but not in the vertical velocity, v_z: there are four prominent and many smaller maxima, many of which correspond to well known moving groups. While samples of early-type stars are dominated by these maxima, also up to 25% of red main-sequence stars are associated with them. These moving groups are responsible for the vertex deviation measured even for samples of late-type stars; they appear more frequently for ever redder samples; and as a whole they follow an asymmetric-drift relation, in the sense that those only present in red samples predominantly have large |v_R| and lag in v_phi w.r.t. the local standard of rest (LSR). The question arise, how these old moving groups got on their eccentric orbits. A plausible mechanism, known from solar system dynamics, which is able to manage a shift in orbit space involves locking into an orbital resonance. Apart from these moving groups, there is a smooth background distribution, akin to Schwarzschild's ellipsoidal model, with axis ratio of about 1:0.6:0.35 in v_R, v_phi, and v_z. The contours are aligned with the $v_r$ direction, but not w.r.t. the v_phi and v_z axes: the mean v_z increases for stars rotating faster than the LSR. This effect can be explained by the stellar warp of the Galactic disk. If this explanation is correct, the warp's inner edge must not be within the solar circle, while its pattern rotates with frequency of about 13 km/s/kpc or more retrograde w.r.t. the stellar orbits.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Local stellar kinematics from Hipparcos data

(shortened) From a kinematically unbiased subsample of the Hipparcos catalogue we have redetermined as a function of colour the kinematics of main-sequence stars. The stars' mean heliocentric velocity nicely follows the asymmetric drift relation, except for stars blueward of B-V=0.1. Extrapolating to zero dispersion yields for the velocity of the Sun w.r.t. the LSR in km/s: U_0=10.00+/-0.36 (radially inwards), V_0=5.23+/-0.62 (in direction of galactic rotation), and W_0=7.17+/-0.38 (vertically upwards). A plot of velocity dispersion vs. colour beautifully shows Parenago's discontinuity: the dispersion is constant for B-V>0.62 and decreases towards bluer colour. We determine the velocity-dispersion tensor sigma^2_ij as function of B-V. The mixed moments involving vertical motion are zero within the errors, while sigma^2_xy is non-zero at about (10km/s)^2 independent of colour. The resulting vertex deviations are about 20 deg for early-type stars and 10+/-4 deg for old-disc stars. The persistence of the vertex deviation to late-type stars implies that the Galactic potential is significantly non-axisymmetric at the solar radius. If spiral arms are responsible for this, they cannot be tightly wound. Except for stars bluer than B-V=0.1 the ratios of the principal velocity dispersions are 2.2 : 1.4 :1, while the absolute values increase with colour from sigma_1=20km/s at B-V=0.2 to sigma_1=38km/s at Parenago's discontinuity and beyond. These ratios imply significant heating of the disc by spiral structure and that R_0/R_d=3 to 3.5, where R_d is the scale length of the disc.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

The outer rotation curve of the Milky Way

A straightforward determination of the circular-speed curve vc(R) of the Milky Way suggests that near the Sun, vc starts to rise approximately linearly with R. If this result were correct, the Galactic mass density would have to be independent of radius at R ~> R0. We show that the apparent linear rise in v_c arises naturally if the true circular-speed curve is about constant or gently falling at R0 < R ~< 2 R0, but most tracers that appear to be at R ~> 1.25 R0 are actually concentrated into a ring of radius ~1.6 R0.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Modelling Galaxies with f(E,Lz); a Black Hole in M32

A technique for the construction of axisymmetric distribution functions for individual galaxies is presented. It starts from the observed surface bright- ness distribution, which is deprojected to gain the axisymmetric luminosity density, from which follows the stars' gravitational potential. After adding dark mass components, such as a central black hole, the two-integral distribu- tion function (2I-DF) f(E,Lz), which depends only on the classical integrals of motion in an axisymmetric potential, is constructed using the Richardson- Lucy algorithm. This algorithm proved to be very efficient in finding f(E,Lz) provided the integral equation to be solved has been properly modified. Once the 2I-\df\ is constructed, its kinematics can be computed and compared with those observed. Many discrepancies may be remedied by altering the assumed inclination angle, mass-to-light ratio, dark components, and odd part of the 2I-DF. Remaining discrepancies may indicate, that the distribution function depends on the non-classical third integral, or is non-axisymmetric. The method has been applied to the nearby elliptical galaxy M32. A 2I-DF with ~55 degrees inclination and a central black hole (or other compact dark mass inside ~1pc) of 1.6-2*10^6 Msun fits the high-spatial-resolution kinema- tic data of van der Marel et al. remarkably well. 2I-DFs with a significantly less or more massive central dark mass or with edge-on inclination can be ruled out for M32. Predictions are made for HST-observations: spectroscopy using its smallest square aperture of 0.09"*0.09" should yield a non-gaussian central velocity profile with broad wings, true and gaussian-fit velocity dispersion of 150-170km/s and 120-130km/s, respectively.

0 citations0 reviewsNo code linkedOpen access