Showing 1–50 of 91 results for author: Dehnen, W

Comment on 'a slightly oblate dark matter halo revealed by a retrograde precessing Galactic disk warp' by Huang et al

Authors: Walter Dehnen, Ralph Schönrich, Ronald Drimmel, Žofia Chrobáková, Eloisa Poggio, Marcin Semczuk

Abstract: Huang et al. (2024) measured the derivative of the phase $φ_{\mathrm{w}}$ of the Galactic warp traced by classical Cepheids with respect to their age $τ$ and interpreted it as the warp precession rate $ω\equiv\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}t=-\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}τ$. This interpretation is unfounded: young stars follow trajectories close to those of their parental gas and trace t… ▽ More Huang et al. (2024) measured the derivative of the phase $φ_{\mathrm{w}}$ of the Galactic warp traced by classical Cepheids with respect to their age $τ$ and interpreted it as the warp precession rate $ω\equiv\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}t=-\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}τ$. This interpretation is unfounded: young stars follow trajectories close to those of their parental gas and trace the instantaneous gas warp, not its shape at their time of birth: $φ_{\mathrm{w}}$ should hardly depend on Cepheid age. We show that the measured $\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}τ>0$ is consistent with an omitted-variable bias from neglecting the natural twist $\mathrm{d}φ_{\mathrm{w}}/\mathrm{d}R$ of the warp and the $R$-$τ$ correlation for Cepheids (originating from the Galactic metallicity gradient and the Cepheid metallicity-age correlation). △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 3 pages, 2 figures, Comment on arXiv:2407.00319

Eccentricity evolution of PTA sources from cosmological initial conditions

Authors: Federica Fastidio, Alessia Gualandris, Alberto Sesana, Elisa Bortolas, Walter Dehnen

Abstract: Recent results from pulsar timing arrays (PTAs) show evidence for a gravitational wave background (GWB) consistent with a population of unresolved supermassive black hole (SMBH) binaries (BHBs). While the data do not yet constrain the slope of the spectrum, this appears to flatten at the lowest frequencies, deviating from the power-law shape expected for circular binaries evolving solely due to gr… ▽ More Recent results from pulsar timing arrays (PTAs) show evidence for a gravitational wave background (GWB) consistent with a population of unresolved supermassive black hole (SMBH) binaries (BHBs). While the data do not yet constrain the slope of the spectrum, this appears to flatten at the lowest frequencies, deviating from the power-law shape expected for circular binaries evolving solely due to gravitational wave (GW) emission. Interestingly, such flattening can be explained with a population of eccentric rather than circular binaries. The eccentricity of BHBs is notoriously difficult to predict based simply on the parameters of the host galaxies and the initial galactic orbit, as it is subject to stochastic effects. We study the evolution of the eccentricity of BHBs formed in galactic mergers with cosmological initial conditions from pairing to coalescence, with a focus on potential PTA sources. We select galactic mergers from the IllustrisTNG100-1 simulation and re-simulate them at high resolution with the N-body code Griffin down to binary separations of the order of a parsec. We then estimate coalescence timescales with a semi-analytical model of the evolution under the effects of GW emission and stellar hardening. We find that most mergers in IllustrisTNG100-1 occur on highly eccentric orbits, and that the eccentricity of BHBs at binary formation correlates with the initial eccentricity of the merger, if this is no larger than approximately 0.9. For extremely eccentric mergers, the binaries tend to form with modest eccentricities. We discuss the implications of these results on the interpretation of the observed GWB. △ Less

Submitted 11 June, 2024; v1 submitted 4 June, 2024; originally announced June 2024.

Comments: 10 pages, 11 figures, Accepted for publication in MNRAS

A Multi-Resolution Method for Modelling Galaxy and Massive Black Hole Mergers

Authors: Kate Attard, Alessia Gualandris, Justin I. Read, Walter Dehnen

Abstract: The coalescence of the most massive black hole (MBH) binaries releases gravitational waves (GWs) within the detectable frequency range of Pulsar Timing Arrays (PTAs) $(10^{-9} - 10^{-6})$ Hz. The incoherent superposition of GWs from MBH mergers, the stochastic Gravitational Wave Background (GWB), can provide unique information on MBH parameters and the large-scale structure of the Universe. The re… ▽ More The coalescence of the most massive black hole (MBH) binaries releases gravitational waves (GWs) within the detectable frequency range of Pulsar Timing Arrays (PTAs) $(10^{-9} - 10^{-6})$ Hz. The incoherent superposition of GWs from MBH mergers, the stochastic Gravitational Wave Background (GWB), can provide unique information on MBH parameters and the large-scale structure of the Universe. The recent evidence for a GWB reported by the PTAs opens an exciting new window onto MBHs and their host galaxies. However, the astrophysical interpretation of the GWB requires accurate estimations of MBH merger timescales for a statistically representative sample of galaxy mergers. This is numerically challenging; a high numerical resolution is required to avoid spurious relaxation and stochastic effects whilst a large number of simulations is needed to sample a cosmologically representative volume. Here, we present a new multi-mass modelling method to increase the central resolution of a galaxy model at a fixed particle number. We follow mergers of galaxies hosting central MBHs with the Fast Multiple Method code Griffin at two reference resolutions and with two refinement schemes. We show that both refinement schemes are effective at increasing central resolution, reducing spurious relaxation and stochastic effects. A particle number of $N\geq 10^{6}$ within a radius of 5 times the sphere of influence of the MBHs is required to reduce numerical scatter in the binary eccentricity and the coalescence timescale to <30$\%$; a resolution that can only be reached at present with the mass refinement scheme. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 13 pages, 18 figures, accepted for publication in MNRAS

Multiple timestep reversible $N$-body integrators for close encounters in planetary systems

Authors: David M. Hernandez, Walter Dehnen

Abstract: We present new almost time-reversible integrators for solution of planetary systems consisting of "planets" and a dominant mass ("star"). The algorithms can be considered adaptive generalizations of the Wisdom--Holman method, in which all pairs of planets can be assigned timesteps. These timesteps, along with the global timestep, can be adapted time-reversibly, often at no appreciable additional c… ▽ More We present new almost time-reversible integrators for solution of planetary systems consisting of "planets" and a dominant mass ("star"). The algorithms can be considered adaptive generalizations of the Wisdom--Holman method, in which all pairs of planets can be assigned timesteps. These timesteps, along with the global timestep, can be adapted time-reversibly, often at no appreciable additional compute cost, without sacrificing any of the long-term error benefits of the Wisdom--Holman method. The method can also be considered a simpler and more flexible version of the \texttt{SYMBA} symplectic code. We perform tests on several challenging problems with close encounters and find the reversible algorithms are up to $2.6$ times faster than a code based on \texttt{SYMBA}. The codes presented here are available on Github. We also find adapting a global timestep reversibly and discretely must be done in block-synchronized manner or similar. △ Less

Submitted 5 April, 2024; v1 submitted 13 January, 2024; originally announced January 2024.

Comments: 11 pages, 7 Figures. Matches accepted MNRAS version

A twisted and precessing Cepheid warp in the outer Milky Way disc

Authors: Walter Dehnen, Marcin Semczuk, Ralph Schönrich

Abstract: We examine the Galactic warp in a sample of all classical Cepheids with Gaia DR3 radial velocity. In each radial bin, we determine (1) the inclined plane normal to the mean orbital angular momentum of the stars and (2) that best fitting their positions. We find no warping inside $R\approx 11$ kpc; for larger $R$ the disc is increasingly inclined, reaching $i\sim 3^{\circ}$ at $R \ge 14$ kpc. With… ▽ More We examine the Galactic warp in a sample of all classical Cepheids with Gaia DR3 radial velocity. In each radial bin, we determine (1) the inclined plane normal to the mean orbital angular momentum of the stars and (2) that best fitting their positions. We find no warping inside $R\approx 11$ kpc; for larger $R$ the disc is increasingly inclined, reaching $i\sim 3^{\circ}$ at $R \ge 14$ kpc. With larger $R$ the azimuth of the warp's ascending node shifts from $\varphi_{\mathrm{lon}}\approx-15^\circ$ at 11 kpc by about $14^{\circ}$/kpc in the direction of Galactic rotation, implying a leading spiral of nodes, the general behaviour of warped galaxies. From the method of fitting planes to the positions we also obtain $\dot{\varphi}_{\mathrm{lon}}$ and find prograde precession of $\dot{\varphi}_{\mathrm{lon}} \sim 12$ km/s/kpc at 12 kpc decreasing to $\sim 6$ km/s/kpc at 14 kpc and beyond. This would unwind the leading spiral of nodes in $\sim 100$ Myr, suggesting that our instantaneous measurements of $\dot{\varphi}_{\mathrm{lon}}$ reflect transient behaviour. This is consistent with existing simulations, which show oscillations in $\dot{\varphi}_{\mathrm{lon}}$ overlaying a long-term retrograde differential precession which generates the leading spiral of nodes. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: 9 pages, accepted for publication in MNRAS

Switching integrators reversibly in the astrophysical $N$-body problem

Authors: David M. Hernandez, Walter Dehnen

Abstract: We present a simple algorithm to switch between $N$-body time integrators in a reversible way. We apply it to planetary systems undergoing arbitrarily close encounters and highly eccentric orbits, but the potential applications are broader. Upgrading an ordinary non-reversible switching integrator to a reversible one is straightforward and introduces no appreciable computational burden in our test… ▽ More We present a simple algorithm to switch between $N$-body time integrators in a reversible way. We apply it to planetary systems undergoing arbitrarily close encounters and highly eccentric orbits, but the potential applications are broader. Upgrading an ordinary non-reversible switching integrator to a reversible one is straightforward and introduces no appreciable computational burden in our tests. Our method checks if the integrator during the time step violates a time-symmetric selection condition and redoes the step if necessary. In our experiments a few percent of steps would have violated the condition without our corrections. By eliminating them the algorithm avoids long-term error accumulation, of several orders magnitude in some cases. △ Less

Submitted 28 February, 2023; v1 submitted 15 January, 2023; originally announced January 2023.

Comments: 10 pages, 8 figures. Matches accepted MNRAS version

Measuring bar pattern speeds from single simulation snapshots

Authors: Walter Dehnen, Marcin Semczuk, Ralph Schoenrich

Abstract: We describe methods to measure simultaneously the orientation angle $ψ$ and pattern speed $Ω$ from single snapshots of simulated barred galaxies. Unlike previous attempts, our approach is unbiased, precise, and consistent in the sense that $ψ=\intΩ\mathrm{d}t$. It can be extended to obtain the rate and axis of rotation, i.e. the vector $\boldsymbolΩ$. We provide computer code implementing our meth… ▽ More We describe methods to measure simultaneously the orientation angle $ψ$ and pattern speed $Ω$ from single snapshots of simulated barred galaxies. Unlike previous attempts, our approach is unbiased, precise, and consistent in the sense that $ψ=\intΩ\mathrm{d}t$. It can be extended to obtain the rate and axis of rotation, i.e. the vector $\boldsymbolΩ$. We provide computer code implementing our method. △ Less

Submitted 7 November, 2022; v1 submitted 1 November, 2022; originally announced November 2022.

Comments: 7 pages, accepted at MNRAS; v2: clarified discussion of previous work

A family of potential-density pairs for galactic bars

Authors: Walter Dehnen, Hossam Aly

Abstract: We present a family of analytical potential-density pairs for barred discs, which can be combined to describe galactic bars in a realistic way, including boxy/peanut components. We illustrate this with two reasonable compound models. Computer code for the evaluation of potential, forces, density, and projected density is freely provided. We present a family of analytical potential-density pairs for barred discs, which can be combined to describe galactic bars in a realistic way, including boxy/peanut components. We illustrate this with two reasonable compound models. Computer code for the evaluation of potential, forces, density, and projected density is freely provided. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Comments: 10 pages, accepted for publication at MNRAS

Spiral arms and the angular momentum gap in Milky Way Cepheids

Authors: Marcin Semczuk, Walter Dehnen, Ralph Schoenrich, E. Athanassoula

Abstract: The angular-momentum distribution of classical Cepheids in the outer Milky Way disc is bi-modal with a gap at $L_\mathrm{gap}=2950\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}$, corresponding to $R=13\,\mathrm{kpc}$, while no similar feature has been found in the general population of disc stars. We show that star formation in multiple spiral arm segments at the same azimuth leads to such multi-mod… ▽ More The angular-momentum distribution of classical Cepheids in the outer Milky Way disc is bi-modal with a gap at $L_\mathrm{gap}=2950\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}$, corresponding to $R=13\,\mathrm{kpc}$, while no similar feature has been found in the general population of disc stars. We show that star formation in multiple spiral arm segments at the same azimuth leads to such multi-modality which quickly dissolves and only shows in young stars. Unlike other explanations, such as a 1:1 orbital resonance with the Galactic bar, this also accounts for the observed steepening of the stellar warp at $L_\mathrm{gap}$, since the adjacent spiral arms represent different parts of the warped gas disc. In this scenario the gap is clearly present only in young stars, as observed, while most purely stellar dynamical origins would affect all disc populations, including older disc stars. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 6 pages, 7 figures, submitted to MNRAS

Age-morphology dependence of the Milky Way boxy/peanut bulge seen in Mira variables

Authors: Marcin Semczuk, Walter Dehnen, Ralph Schoenrich, E. Athanassoula

Abstract: We analyse the distribution of Mira variable stars in the central region of the Milky Way. We find that with increasing period, i.e. decreasing age, the Miras shift towards negative Galactic longitudes $\ell$. Comparing to a cosmological zoom simulation of a barred galaxy, we find that this shift with age can be explained by an age-morphology dependence of the boxy peanut/bulge. Owing to a combina… ▽ More We analyse the distribution of Mira variable stars in the central region of the Milky Way. We find that with increasing period, i.e. decreasing age, the Miras shift towards negative Galactic longitudes $\ell$. Comparing to a cosmological zoom simulation of a barred galaxy, we find that this shift with age can be explained by an age-morphology dependence of the boxy peanut/bulge. Owing to a combination of projection effects and the limitation of the range of Galactic longitudes, the near hump at $\ell>0$ is more truncated for younger populations, and the far hump at $\ell<0$ dominates the observed distributions. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 5 pages, 4 figures, submitted to MNRAS

Capture of interstellar objects II: by the Solar system

Authors: Walter Dehnen, Thomas O. Hands, Ralph Schönrich

Abstract: 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 cr… ▽ More 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. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 8 pages, to be published by MNRAS

Capture of interstellar objects I: the capture cross-section

Authors: Walter Dehnen, Thomas O. Hands

Abstract: 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 intera… ▽ More 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. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: 16 pages, to be published in MNRAS

The small boxy/peanut structure of the Milky Way traced by old stars

Authors: Marcin Semczuk, Walter Dehnen, Ralph Schoenrich, E. Athanassoula

Abstract: We analyse the positions of RR Lyrae stars in the central region of the Milky Way. In addition to the overall bar shape detected previously, we find evidence for a peanut shaped structure, in form of overdensities near $\ell=-2$ deg and $1$ deg at $b\sim-3.5$ deg. The corresponding physical distance between the two peaks of the peanut is $\sim0.7$ kpc, significantly shorter than that found from ne… ▽ More We analyse the positions of RR Lyrae stars in the central region of the Milky Way. In addition to the overall bar shape detected previously, we find evidence for a peanut shaped structure, in form of overdensities near $\ell=-2$ deg and $1$ deg at $b\sim-3.5$ deg. The corresponding physical distance between the two peaks of the peanut is $\sim0.7$ kpc, significantly shorter than that found from near-IR images (3.3 kpc) and red-clump stars. Qualitatively this is expected from `fractionation' of bar orbits, which we demonstrate to be matched in a simulation of an inside-out growing disc subsequently forming a bar. △ Less

Submitted 11 November, 2021; v1 submitted 22 July, 2021; originally announced July 2021.

Comments: 6 pages, 5 figures, accepted to MNRAS

EDGE: Two routes to dark matter core formation in ultra-faint dwarfs

Authors: Matthew D. A. Orkney, Justin I. Read, Martin P. Rey, Imran Nasim, Andrew Pontzen, Oscar Agertz, Stacy Y. Kim, Maxime Delorme, Walter Dehnen

Abstract: In the standard Lambda cold dark matter paradigm, pure dark matter simulations predict dwarf galaxies should inhabit dark matter haloes with a centrally diverging density `cusp'. This is in conflict with observations that typically favour a constant density `core'. We investigate this `cusp-core problem' in `ultra-faint' dwarf galaxies simulated as part of the `Engineering Dwarfs at Galaxy formati… ▽ More In the standard Lambda cold dark matter paradigm, pure dark matter simulations predict dwarf galaxies should inhabit dark matter haloes with a centrally diverging density `cusp'. This is in conflict with observations that typically favour a constant density `core'. We investigate this `cusp-core problem' in `ultra-faint' dwarf galaxies simulated as part of the `Engineering Dwarfs at Galaxy formation's Edge' (EDGE) project. We find, similarly to previous work, that gravitational potential fluctuations within the central region of the simulated dwarfs kinematically heat the dark matter particles, lowering the dwarfs' central dark matter density. However, these fluctuations are not exclusively caused by gas inflow/outflow, but also by impulsive heating from minor mergers. We use the genetic modification approach on one of our dwarf's initial conditions to show how a delayed assembly history leads to more late minor mergers and, correspondingly, more dark matter heating. This provides a mechanism by which even ultra-faint dwarfs ($M_* < 10^5\,\text{M}_{\odot}$), in which star formation was fully quenched at high redshift, can have their central dark matter density lowered over time. In contrast, we find that late major mergers can regenerate a central dark matter cusp, if the merging galaxy had sufficiently little star formation. The combination of these effects leads us to predict significant stochasticity in the central dark matter density slopes of the smallest dwarfs, driven by their unique star formation and mass assembly histories. △ Less

Submitted 7 January, 2021; originally announced January 2021.

Comments: 14 pages, 13 figures, submitted to MNRAS

Breaking Beta: A comparison of mass modelling methods for spherical systems

Authors: J. I. Read, G. A. Mamon, E. Vasiliev, L. L. Watkins, M. G. Walker, J. Penarrubia, M. Wilkinson, W. Dehnen, P. Das

Abstract: We apply four different mass modelling methods to a suite of publicly available mock data for spherical stellar systems. We focus on the recovery of the density and velocity anisotropy as a function of radius, using either line-of-sight velocity data only, or adding proper motion data. All methods perform well on isotropic and tangentially anisotropic mock data, recovering the density and velocity… ▽ More We apply four different mass modelling methods to a suite of publicly available mock data for spherical stellar systems. We focus on the recovery of the density and velocity anisotropy as a function of radius, using either line-of-sight velocity data only, or adding proper motion data. All methods perform well on isotropic and tangentially anisotropic mock data, recovering the density and velocity anisotropy within their 95% confidence intervals over the radial range 0.25 < R/Rhalf < 4, where Rhalf is the half light radius. However, radially-anisotropic mocks are more challenging. For line-of-sight data alone, only methods that use information about the shape of the velocity distribution function are able to break the degeneracy between the density profile and the velocity anisotropy to obtain an unbiased estimate of both. This shape information can be obtained through directly fitting a global phase space distribution function, by using higher order 'Virial Shape Parameters', or by assuming a Gaussian velocity distribution function locally, but projecting it self-consistently along the line of sight. Including proper motion data yields further improvements, and in this case, all methods give a good recovery of both the radial density and velocity anisotropy profiles. △ Less

Submitted 20 November, 2020; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: 17 pages, 8 figures. Accepted for publication in MNRAS

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

Authors: Imran Nasim, Alessia Gualandris, Justin I. Read, Fabio Antonini, Walter Dehnen, Maxime Delorme

Abstract: 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 ex… ▽ More 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. △ Less

Submitted 11 February, 2021; v1 submitted 9 November, 2020; originally announced November 2020.

Comments: Accepted in MNRAS

Secular Dynamics around a Supermassive Black Hole via Multipole Expansion

Authors: Jean-Baptiste Fouvry, Walter Dehnen, Scott Tremaine, Ben Bar-Or

Abstract: 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 ti… ▽ More 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. △ Less

Submitted 1 March, 2022; v1 submitted 3 November, 2020; originally announced November 2020.

Comments: 20 pages, 11 figures, submitted to AAS

Non-parametric Triaxial Deprojection of Elliptical Galaxies

Authors: Stefano de Nicola, Roberto P. Saglia, Jens Thomas, Walter Dehnen, Ralf Bender

Abstract: 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 l… ▽ More 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. △ Less

Submitted 7 July, 2020; v1 submitted 10 June, 2020; originally announced June 2020.

Comments: 26 Pages, 29 Figures, accepted by MNRAS

Journal ref: MNRAS, 2020, 3076, 3100

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

Authors: Imran Nasim, Alessia Gualandris, Justin Read, Walter Dehnen, Maxime Delorme, Fabio Antonini

Abstract: 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 mass… ▽ More 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. △ Less

Submitted 29 April, 2020; originally announced April 2020.

Comments: Submitted for publication in MNRAS

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

Authors: Tom Hands, Walter Dehnen

Abstract: 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 velo… ▽ More 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. △ Less

Submitted 6 January, 2020; v1 submitted 14 October, 2019; originally announced October 2019.

Comments: 5 pages, 3 figures, accepted for publication in MNRAS

Giant planets and brown dwarfs on wide orbits: a code comparison project

Authors: Mark Fletcher, Sergei Nayakshin, Dimitris Stamatellos, Walter Dehnen, Farzana Meru, Lucio Mayer, Hongping Deng, Ken Rice

Abstract: Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs and companion stars. Simulations show that competition between three processes -- migration, gas accretion and tidal disruption -- establishes what grows from a given seed. Here we investigate the robustness of numerical modelling of clump migration and accretion wit… ▽ More Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs and companion stars. Simulations show that competition between three processes -- migration, gas accretion and tidal disruption -- establishes what grows from a given seed. Here we investigate the robustness of numerical modelling of clump migration and accretion with the codes PHANTOM, GADGET, SPHINX, SEREN, GIZMO-MFM, SPHNG and FARGO. The test problem comprises a clump embedded in a massive disc at an initial separation of 120 AU. There is a general qualitative agreement between the codes, but the quantitative agreement in the planet migration rate ranges from $\sim 10$% to $\sim 50$%, depending on the numerical setup. We find that the artificial viscosity treatment and the sink particle prescription may account for much of the differences between the codes. In order to understand the wider implications of our work, we also attempt to reproduce the planet evolution tracks from our hydrodynamical simulations with prescriptions from three previous population synthesis studies. We find that the disagreement amongst the population synthesis models is far greater than that between our hydrodynamical simulations. The results of our code comparison project are therefore encouraging in that uncertainties in the given problem are probably dominated by the physics not yet included in the codes rather than by how hydrodynamics is modelled in them. △ Less

Submitted 23 January, 2019; originally announced January 2019.

Comments: submitted to MNRAS. This is version 2 of the paper after considering referees comments and changes

The fate of planetesimal discs in young open clusters: implications for 1I/'Oumuamua, the Kuiper belt, the Oort cloud and more

Authors: Thomas Oliver Hands, Walter Dehnen, Amery Gration, Joachim Stadel, Ben Moore

Abstract: We perform $N$-body simulations of the early phases of open cluster evolution including a large population of planetesimals, initially arranged in Kuiper-belt like discs around each star. Using a new, 4th-order and time-reversible $N$-body code on Graphics Processing Units (GPUs), we evolve the whole system under the stellar gravity, i.e. treating planetesimals as test particles, and consider two… ▽ More We perform $N$-body simulations of the early phases of open cluster evolution including a large population of planetesimals, initially arranged in Kuiper-belt like discs around each star. Using a new, 4th-order and time-reversible $N$-body code on Graphics Processing Units (GPUs), we evolve the whole system under the stellar gravity, i.e. treating planetesimals as test particles, and consider two types of initial cluster models, similar to IC348 and the Hyades, respectively. In both cases, planetesimals can be dynamically excited, transferred between stars or liberated to become free-floating (such as A/2017 U1 or 'Oumuamua) during the early cluster evolution. We find that planetesimals captured from another star are not necessarily dynamically distinct from those native to a star. After an encounter both native and captured planetesimals can exhibit aligned periastrons, qualitatively similar to that seen in the Solar system and commonly thought to be the signature of Planet 9. We discuss the implications of our results for both our Solar system and exoplanetary systems. △ Less

Submitted 12 April, 2019; v1 submitted 8 January, 2019; originally announced January 2019.

Comments: 15 pages, 7 figures, 4 tables, accepted for publication in MNRAS

How gravitational softening affects galaxy stability I. Linear mode analysis of disc galaxies

Authors: Sven De Rijcke, Jean-Baptiste Fouvry, Walter Dehnen

Abstract: Linear perturbation is used to investigate the effect of gravitational softening on the retrieved two-armed spiral eigenmodes of razor-thin stellar discs. We explore four softening kernels with different degrees of gravity bias, and with/without compact support (compact in the sense that they yield exactly Newtonian forces outside the softening kernel). These kernels are applied to two disc gala… ▽ More Linear perturbation is used to investigate the effect of gravitational softening on the retrieved two-armed spiral eigenmodes of razor-thin stellar discs. We explore four softening kernels with different degrees of gravity bias, and with/without compact support (compact in the sense that they yield exactly Newtonian forces outside the softening kernel). These kernels are applied to two disc galaxy models with well-known unsoftened unstable modes. We illustrate quantitatively the importance of a vanishing linear gravity bias to yield accurate frequency estimates of the unstable modes. As such, Plummer softening, while very popular amongst simulators, performs poorly in our tests. The best results, with excellent agreement between the softened and unsoftened mode properties, are obtained with softening kernels that have a reduced gravity bias, obtained by compensating for the sub-Newtonian forces at small interparticle distances with slightly super-Newtonian forces at radii near the softening length. We present examples of such kernels that, moreover, are analytically simple and computationally cheap. Finally, these results light the way to the construction of softening methods with even smaller gravity bias, although at the price of increasingly complex kernels. △ Less

Submitted 17 December, 2018; originally announced December 2018.

Comments: 14 pages, 7 figures, submitted to MNRAS

The hidden giant: discovery of an enormous Galactic dwarf satellite in Gaia DR2

Authors: G. Torrealba, V. Belokurov, S. E. Koposov, T. S. Li, M. G. Walker, J. L. Sanders, A. Geringer-Sameth, D. B. Zucker, K. Kuehn, N. W. Evans, W. Dehnen

Abstract: We report the discovery of a Milky Way satellite in the constellation of Antlia. The Antlia 2 dwarf galaxy is located behind the Galactic disc at a latitude of $b\sim 11^{\circ}$ and spans $1.26$ degrees, which corresponds to $\sim2.9$ kpc at its distance of 130 kpc. While similar in spatial extent to the Large Magellanic Cloud, Antlia 2 is orders of magnitude fainter at $M_V=-9$ mag, making it by… ▽ More We report the discovery of a Milky Way satellite in the constellation of Antlia. The Antlia 2 dwarf galaxy is located behind the Galactic disc at a latitude of $b\sim 11^{\circ}$ and spans $1.26$ degrees, which corresponds to $\sim2.9$ kpc at its distance of 130 kpc. While similar in spatial extent to the Large Magellanic Cloud, Antlia 2 is orders of magnitude fainter at $M_V=-9$ mag, making it by far the lowest surface brightness system known (at $\sim31.9$ mag/arcsec$^2$), $\sim100$ times more diffuse than the so-called ultra diffuse galaxies. The satellite was identified using a combination of astrometry, photometry and variability data from \textit{Gaia} Data Release 2, and its nature confirmed with deep archival DECam imaging, which revealed a conspicuous BHB signal. We have also obtained follow-up spectroscopy using AAOmega on the AAT, identifying 159 member stars, and we used them to measure the dwarf's systemic velocity, $290.9\pm0.5$km/s, its velocity dispersion, $5.7\pm1.1$ km/s, and mean metallicity, [Fe/H]$=-1.4$. From these properties we conclude that Antlia 2 inhabits one of the least dense Dark Matter (DM) halos probed to date. Dynamical modelling and tidal-disruption simulations suggest that a combination of a cored DM profile and strong tidal stripping may explain the observed properties of this satellite. The origin of this core may be consistent with aggressive feedback, or may even require alternatives to cold dark matter (such as ultra-light bosons). △ Less

Submitted 11 June, 2019; v1 submitted 9 November, 2018; originally announced November 2018.

Comments: 24 pages, 17 figures, 2 tables. Accepted in MNRAS. v2 :fixed RRL Distance, main results remain the same. v3 Accepted version. fixed Mv estimation, main results remain the same

Tidal ribbons

Authors: Walter Dehnen, Hasanuddin

Abstract: Tidal debris from Galactic satellites generally forms one-dimensional elongated streams, since nearby Galactic orbits have almost identical frequency ratios. We show that the situation is different for orbits close to the Galactic disc, whose vertical frequency $Ω_z$ is strongly amplitude dependent. As a consequence, stars stripped from a satellite obtain a range of values for $Ω_z$ and hence of f… ▽ More Tidal debris from Galactic satellites generally forms one-dimensional elongated streams, since nearby Galactic orbits have almost identical frequency ratios. We show that the situation is different for orbits close to the Galactic disc, whose vertical frequency $Ω_z$ is strongly amplitude dependent. As a consequence, stars stripped from a satellite obtain a range of values for $Ω_z$ and hence of frequency ratios, and spread into two dimensions, forming a ribbon-like structure with vertical extent comparable to that of the progenitor orbit. In integrals-of-motion space, tidal ribbons are clumps, which offers the best chance of detection and allows the determination of the Galactic potential vertically across the disc. △ Less

Submitted 28 June, 2018; v1 submitted 22 May, 2018; originally announced May 2018.

Comments: 7 pages, 6 figures, accepted for MNRAS

Feeding supermassive black holes by collisional cascades

Authors: Christian Faber, Walter Dehnen

Abstract: The processes driving gas accretion on to supermassive black holes (SMBHs) are still poorly understood. Angular momentum conservation prevents gas within $\sim10\,$pc of the black hole from reaching radii $\sim10^{-3}\,$pc where viscous accretion becomes efficient. Here we present simulations of the collapse of a clumpy shell of swept-up isothermal gas, which is assumed to have formed as a result… ▽ More The processes driving gas accretion on to supermassive black holes (SMBHs) are still poorly understood. Angular momentum conservation prevents gas within $\sim10\,$pc of the black hole from reaching radii $\sim10^{-3}\,$pc where viscous accretion becomes efficient. Here we present simulations of the collapse of a clumpy shell of swept-up isothermal gas, which is assumed to have formed as a result of feedback from a previous episode of AGN activity. The gas falls towards the SMBH forming clumps and streams, which intersect, collide, and often form a disc. These collisions promote partial cancellations of angular momenta, resulting in further infall and more collisions. This continued collisional cascade generates a tail of gas with sufficiently small angular momenta and provides a viable route for gas inflow to sub-parsec scales. The efficiency of this process hardly depends on details, such as gas temperature, initial virial ratio and power spectrum of the gas distribution, as long as it is not strongly rotating. Adding star formation to this picture might explain the near-simultaneous formation of the S-stars (from tidally disrupted binaries formed in plunging gas streams) and the sub-parsec young stellar disc around Sgr A$^{\!\star}$. △ Less

Submitted 23 April, 2018; originally announced April 2018.

Comments: 15 pages, 14 figures, accepted for publication in MNRAS

Tidal disruption of dwarf spheroidal galaxies: the strange case of Crater II

Authors: Jason L. Sanders, N. Wyn Evans, Walter Dehnen

Abstract: Dwarf spheroidal galaxies of the Local Group obey a relationship between the line-of-sight velocity dispersion and half-light radius, although there are a number of dwarfs that lie beneath this relation with suppressed velocity dispersion. The most discrepant of these (in the Milky Way) is the `feeble giant' Crater II. Using analytic arguments supported by controlled numerical simulations of tidal… ▽ More Dwarf spheroidal galaxies of the Local Group obey a relationship between the line-of-sight velocity dispersion and half-light radius, although there are a number of dwarfs that lie beneath this relation with suppressed velocity dispersion. The most discrepant of these (in the Milky Way) is the `feeble giant' Crater II. Using analytic arguments supported by controlled numerical simulations of tidally-stripped flattened two-component dwarf galaxies, we investigate interpretations of Crater II within standard galaxy formation theory. Heavy tidal disruption is necessary to explain the velocity-dispersion suppression which is plausible if the proper motion of Crater II is $(μ_{α*},μ_δ)=(-0.21\pm0.09,-0.24\pm0.09)\mathrm{mas\,yr}^{-1}$. Furthermore, we demonstrate that the velocity dispersion of tidally-disrupted systems is solely a function of the total mass loss even for weakly-embedded and flattened systems. The half-light radius evolution depends more sensitively on orbital phase and the properties of the dark matter profile. The half-light radius of weakly-embedded cusped systems rapidly decreases producing some tension with the Crater II observations. This tension is alleviated by cored dark matter profiles, in which the half-light radius can grow after tidal disruption. The evolution of flattened galaxies is characterised by two competing effects: tidal shocking makes the central regions rounder whilst tidal distortion produces a prolate tidally-locked outer envelope. After $\sim70$ per cent of the central mass is lost, tidal distortion becomes the dominant effect and the shape of the central regions of the galaxy tends to a universal prolate shape irrespective of the initial shape. △ Less

Submitted 26 February, 2018; originally announced February 2018.

Comments: submitted to MNRAS

Warp, Waves, and Wrinkles in the Milky Way

Authors: Ralph Schönrich, Walter Dehnen

Abstract: We derive unbiased distance estimates for the Gaia-TGAS dataset by correcting for the bias due to the distance dependence of the selection function, which we measure directly from the data. From these distances and proper motions, we estimate the vertical and azimuthal velocities, $W$ and $V_φ$, and angular momentum $L_z$ for stars in the Galactic centre and anti-centre directions. The resulting m… ▽ More We derive unbiased distance estimates for the Gaia-TGAS dataset by correcting for the bias due to the distance dependence of the selection function, which we measure directly from the data. From these distances and proper motions, we estimate the vertical and azimuthal velocities, $W$ and $V_φ$, and angular momentum $L_z$ for stars in the Galactic centre and anti-centre directions. The resulting mean vertical motion $\overline{W}$ shows a linear increase with both $V_φ$ and $L_z$ at $10 σ$ significance. Such a trend is expected from and consistent with the known Galactic warp. This signal extends to stars with guiding centre radii $R_g<R_0$, placing the onset of the warp at $R\lesssim7{\rm kpc}$. At equally high significance, we detect a previously unknown wave-like pattern of $\overline{W}$ over guiding centre $R_g$ with amplitude $\sim1{\rm kms}^{-1}$ and wavelength $\sim2.5{\rm kpc}$. This pattern is present in both the centre and anti-centre directions, consistent with a winding (corrugated) warp or bending wave, likely related to known features in the outer disc (TriAnd and Monoceros over-densities), and may be caused by the interaction with the Sgr dwarf galaxy $\sim1{\rm Gyr}$ ago. The only significant deviation from this simple fit is a stream-like feature near $R_g\sim9{\rm kpc}$ ($|L_z|\sim2150{\rm kpckms}^{-1}$). △ Less

Submitted 19 June, 2018; v1 submitted 18 December, 2017; originally announced December 2017.

Comments: 17 pages, 15 figures, MNRAS in press, some additional tests

Towards time symmetric N-body integration

Authors: Walter Dehnen

Abstract: Computational efficiency demands discretised, hierarchically organised, and individually adaptive time-step sizes (known as the block-step scheme) for the time integration of N-body models. However, most existing N-body codes adapt individual step sizes in a way that violates time symmetry (and symplecticity), resulting in artificial secular dissipation (and often secular growth of energy errors).… ▽ More Computational efficiency demands discretised, hierarchically organised, and individually adaptive time-step sizes (known as the block-step scheme) for the time integration of N-body models. However, most existing N-body codes adapt individual step sizes in a way that violates time symmetry (and symplecticity), resulting in artificial secular dissipation (and often secular growth of energy errors). Using single-orbit integrations, I investigate various possibilities to reduce or eliminate irreversibility from the time stepping scheme. Significant improvements over the standard approach are possible at little extra effort. However, in order to reduce irreversible step-size changes to negligible amounts, such as suitable for long-term integrations of planetary systems, more computational effort is needed, while exact time reversibility appears elusive for discretised individual step sizes. △ Less

Submitted 27 July, 2017; originally announced July 2017.

Comments: 13 pages, 16 figures, accepted for publication in MNRAS

A study of symplectic integrators for planetary system problems: error analysis and comparisons

Authors: David M. Hernandez, Walter Dehnen

Abstract: The symplectic Wisdom-Holman map revolutionized long-term integrations of planetary systems. There is freedom in such methods of how to split the Hamiltonian and which coordinate system to employ, and several options have been proposed in the literature. These choices lead to different integration errors, which we study analytically and numerically. The Wisdom-Holman method in Jacobi coordinates a… ▽ More The symplectic Wisdom-Holman map revolutionized long-term integrations of planetary systems. There is freedom in such methods of how to split the Hamiltonian and which coordinate system to employ, and several options have been proposed in the literature. These choices lead to different integration errors, which we study analytically and numerically. The Wisdom-Holman method in Jacobi coordinates and the method of Hernandez, H16, compare favorably and avoid problems of some of the other maps, such as incorrect center-of-mass position or truncation errors even in the one-planet case. We use H16 to compute the evolution of Pluto's orbital elements over 500 million years in a new calculation. △ Less

Submitted 1 March, 2017; v1 submitted 15 December, 2016; originally announced December 2016.

Comments: 24 pages, 9 figures, 1 table. Matches accepted MNRAS version

Journal ref: MNRAS (2017) 468 (3): 2614-2636

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

Authors: Alessia Gualandris, Justin I. Read, Walter Dehnen, Elisa Bortolas

Abstract: 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 mor… ▽ More 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. △ Less

Submitted 29 September, 2016; originally announced September 2016.

Comments: 11 pages, 10 figures, accepted for publication in MNRAS

Symplectic fourth-order maps for the collisional N-body problem

Authors: Walter Dehnen, David M. Hernandez

Abstract: We study analytically and experimentally certain symplectic and time-reversible N-body integrators which employ a Kepler solver for each pair-wise interaction, including the method of Hernandez & Bertschinger (2015). Owing to the Kepler solver, these methods treat close two-body interactions correctly, while close three-body encounters contribute to the truncation error at second order and above.… ▽ More We study analytically and experimentally certain symplectic and time-reversible N-body integrators which employ a Kepler solver for each pair-wise interaction, including the method of Hernandez & Bertschinger (2015). Owing to the Kepler solver, these methods treat close two-body interactions correctly, while close three-body encounters contribute to the truncation error at second order and above. The second-order errors can be corrected to obtain a fourth-order scheme with little computational overhead. We generalise this map to an integrator which employs a Kepler solver only for selected interactions and yet retains fourth-order accuracy without backward steps. In this case, however, two-body encounters not treated via a Kepler solver contribute to the truncation error. △ Less

Submitted 29 September, 2016; originally announced September 2016.

Comments: 17 pages, re-submitted to MNRAS

Journal ref: MNRAS (2017) 465 (1): 1201-1217

Indirect Dark Matter Detection for Flattened Dwarf Galaxies

Authors: Jason L. Sanders, N. Wyn Evans, Alex Geringer-Sameth, Walter Dehnen

Abstract: 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 effe… ▽ More 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) △ Less

Submitted 16 September, 2016; v1 submitted 19 April, 2016; originally announced April 2016.

Comments: 22 pages, 10 figures, matches accepted version

Journal ref: Phys. Rev. D 94, 063521 (2016)

On the Alignment of the Stress Tensor in Galaxies

Authors: N. W. Evans, J. L. Sanders, A. A. Williams, J. An, D. Lynden-Bell, W. Dehnen

Abstract: We show that, provided the principal axes of the second velocity moment tensor of a stellar population are generally unequal and are oriented perpendicular to a set of orthogonal surfaces at each point, then those surfaces must be confocal quadric surfaces and the potential must be separable or Stackel. This is true under the mild assumption that the even part of the distribution function is invar… ▽ More We show that, provided the principal axes of the second velocity moment tensor of a stellar population are generally unequal and are oriented perpendicular to a set of orthogonal surfaces at each point, then those surfaces must be confocal quadric surfaces and the potential must be separable or Stackel. This is true under the mild assumption that the even part of the distribution function is invariant under time reversal $v_i \rightarrow -v_i$ of each velocity component. In particular, if the second velocity moment tensor is everywhere exactly aligned in spherical polar coordinates, then the potential must be of Stackel form (excepting degenerate cases where two or more of the semiaxes of ellipsoid are everywhere the same). The theorem also has consequences for alignment in cylindrical polar coordinates, which is used in the popular Jeans Anisotropic Models (JAM). We analyse data on the radial velocities and proper motions of a sample of $\sim 7400$ stars in the stellar halo of the Milky Way. We provide the distributions of the tilt angles or misalignments from the spherical polar coordinate systems. We show that in this sample the misalignment is always small (usually within $3^\circ$) for radii between 7 and 12 kpc. The velocity anisotropy is very radially biased ($β\approx 0.7$), and almost invariant across the volume in our study. Finally, we construct a triaxial stellar halo in a triaxial NFW dark matter halo using a made-to-measure method. Despite the triaxiality of the potential, the velocity ellipsoid of the stellar halo is nearly spherically aligned within $\sim6^\circ$ for large regions of space, particularly outside the scale radius of the stellar halo. We conclude that the second velocity moment ellipsoid can be close to spherically aligned for a much wider class of potentials than the strong constraints that arise from exact alignment might suggest. △ Less

Submitted 24 November, 2015; v1 submitted 29 September, 2015; originally announced September 2015.

Comments: MNRAS, in press (Section 2, which comprises the analysis of the data, has been largely re-done from the first version)

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

Authors: Hossam Aly, Walter Dehnen, Chris Nixon, Andrew King

Abstract: 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 fou… ▽ More 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. △ Less

Submitted 19 January, 2015; originally announced January 2015.

Comments: 12 pages, 9 figures, Accepted for publication in MNRAS

Understanding the assembly of Kepler's compact planetary systems

Authors: T. O. Hands, R. D. Alexander, W. Dehnen

Abstract: The Kepler mission has recently discovered a number of exoplanetary systems, such as Kepler-11 and Kepler-32, in which ensembles of several planets are found in very closely packed orbits (often within a few percent of an AU of one another). These compact configurations present a challenge for traditional planet formation and migration scenarios. We present a dynamical study of the assembly of the… ▽ More The Kepler mission has recently discovered a number of exoplanetary systems, such as Kepler-11 and Kepler-32, in which ensembles of several planets are found in very closely packed orbits (often within a few percent of an AU of one another). These compact configurations present a challenge for traditional planet formation and migration scenarios. We present a dynamical study of the assembly of these systems, using an N-body method which incorporates a parametrized model of planet migration in a turbulent protoplanetary disc. We explore a wide parameter space, and find that under suitable conditions it is possible to form compact, close-packed planetary systems via traditional disc-driven migration. We find that simultaneous migration of multiple planets is a viable mechanism for the assembly of tightly-packed planetary systems, as long as the disc provides significant eccentricity damping and the level of turbulence in the disc is modest. We discuss the implications of our preferred parameters for the protoplanetary discs in which these systems formed, and comment on the occurrence and significance of mean-motion resonances in our simulations. △ Less

Submitted 1 September, 2014; originally announced September 2014.

Comments: 12 pages, 4 figures, 2 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

A fast multipole method for stellar dynamics

Authors: Walter Dehnen

Abstract: 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 fr… ▽ More 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$. △ Less

Submitted 9 May, 2014; originally announced May 2014.

Comments: 21 pages, 15 figures, accepted for publication in Journal for Computational Astrophysics and Cosmology

Black hole foraging: feedback drives feeding

Authors: Walter Dehnen, Andrew King

Abstract: 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 s… ▽ More 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. △ Less

Submitted 8 October, 2013; originally announced October 2013.

Comments: accepted for publication in ApJ Letters, 5 pages, 1 figure

Constraints on the Galactic Bar from the Hercules stream as traced with RAVE across the Galaxy

Authors: T. Antoja, A. Helmi, W. Dehnen, O. Bienayme, J. Bland-Hawthorn, B. Famaey, K. Freeman, B. K. Gibson, G. Gilmore, E. K. Grebel, G. Kordopatis, A. Kunder, I. Minchev, U. Munari, J. Navarro, Q. Parker, W. A. Reid, G. Seabroke, A. Siebert, M. Steinmetz, F. Watson, R. F. G. Wyse, T. Zwitter

Abstract: Non-axisymmetries in the Galactic potential (spiral arms and bar) induce kinematic groups such as the Hercules stream. Assuming that Hercules is caused by the effects of the Outer Lindblad Resonance of the Galactic bar, we model analytically its properties as a function of position in the Galaxy and its dependence on the bar's pattern speed and orientation. Using data from the RAVE survey we find… ▽ More Non-axisymmetries in the Galactic potential (spiral arms and bar) induce kinematic groups such as the Hercules stream. Assuming that Hercules is caused by the effects of the Outer Lindblad Resonance of the Galactic bar, we model analytically its properties as a function of position in the Galaxy and its dependence on the bar's pattern speed and orientation. Using data from the RAVE survey we find that the azimuthal velocity of the Hercules structure decreases as a function of Galactocentric radius, in a manner consistent with our analytical model. This allows us to obtain new estimates of the parameters of the Milky Way's bar. The combined likelihood function of the bar's pattern speed and angle has its maximum for a pattern speed of $Ω_b=(1.89 \pm 0.08) \times Ω_0$ where $Ω_0$ is the local circular frequency. Assuming a Solar radius of 8.05 kpc and a local circular velocity of 238 km/s, this corresponds to $Ω_b=56\pm 2$ km/s/kpc. On the other hand, the bar's orientation $φ_b$ cannot be constrained with the available data. In fact, the likelihood function shows that a tight correlation exists between the pattern speed and the orientation, implying that a better description of our best fit results is given by the linear relation $Ω_b/Ω_0=1.91 + 0.0044[φ_b(°)-48]$, with standard deviation of 0.02. For example, for an angle of $φ_b=30°$ the pattern speed is $54.0\pm0.5$ km/s/kpc. These results are not very sensitive to the other Galactic parameters such as the circular velocity curve or the peculiar motion of the Sun, and are robust to biases in distance. △ Less

Submitted 28 January, 2014; v1 submitted 17 September, 2013; originally announced September 2013.

Comments: 16 pages, 16 figures, 3 tables, 2 appendix (2 tables, 2 figures), accepted for publication in A&A

Gravitational Collapse in One Dimension

Authors: A. E. Schulz, Walter Dehnen, Gerard Jungman, Scott Tremaine

Abstract: 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 prof… ▽ More 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. △ Less

Submitted 1 June, 2012; originally announced June 2012.

Comments: 14 Pages, 13 Figures. Submitted to MNRAS

Radial Migration Does Little for Galactic Disc Thickening

Authors: I. Minchev, B. Famaey, A. C. Quillen, W. Dehnen, M. Martig, A. Siebert

Abstract: Non-axisymmetric components, such as spirals and central bars, play a major role in shaping galactic discs. An important aspect of the disc secular evolution driven by these perturbers is the radial migration of stars. It has been suggested recently that migration can populate a thick-disc component from inner-disc stars with high vertical energies. Since this has never been demonstrated in simula… ▽ More Non-axisymmetric components, such as spirals and central bars, play a major role in shaping galactic discs. An important aspect of the disc secular evolution driven by these perturbers is the radial migration of stars. It has been suggested recently that migration can populate a thick-disc component from inner-disc stars with high vertical energies. Since this has never been demonstrated in simulations, we study in detail the effect of radial migration on the disc velocity dispersion and disc thickness, by separating simulated stars into migrators and non-migrators. We apply this method to three isolated barred Tree-SPH N-body galaxies with strong radial migration. Contrary to expectations, we find that as stellar samples migrate, on the average, their velocity dispersion change (by as much as 50%) in such a way as to approximately match the non-migrating population at the radius at which they arrive. We show that, in fact, migrators suppress heating in parts of the disc. To confirm the validity of our findings, we also apply our technique to three cosmological re-simulations, which use a completely different simulation scheme and, remarkably, find very similar results. We believe the inability of migration to thicken discs is a fundamental property of internal disc evolution, irrespective of the migration mechanism at work. We explain this with the approximate conservation of the (average) vertical and radial actions rather than the energy. This "action mixing" can be used to constrain the migration rate in the Milky Way: estimates of the average vertical action in observations for different populations of stars should reveal flattening with radius for older groups of stars. △ Less

Submitted 15 October, 2012; v1 submitted 29 May, 2012; originally announced May 2012.

Comments: 18 Pages, 13 Figures, A&A Accepted with minor revisions: a paragraph added in Sec.10 discussing the effect of resolution. Higher resolution version can be found at www.ivanminchev.com

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

Authors: David R. Cole, Walter Dehnen, Justin I. Read, Mark I. Wilkinson

Abstract: 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 tha… ▽ More 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. △ Less

Submitted 29 May, 2012; originally announced May 2012.

Comments: 12 pages, 8 figures, 3 tables, submitted to MNRAS

Improving convergence in smoothed particle hydrodynamics simulations without pairing instability

Authors: Walter Dehnen, Hossam Aly

Abstract: 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 insuffic… ▽ More 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. △ Less

Submitted 1 June, 2012; v1 submitted 11 April, 2012; originally announced April 2012.

Comments: substantially revised version, accepted for publication in MNRAS, 15 pages, 13 figures

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

Authors: Ivan Minchev, Benoit Famaey, Alice C. Quillen, Walter Dehnen

Abstract: 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 pa… ▽ More 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. △ Less

Submitted 7 November, 2011; v1 submitted 1 November, 2011; originally announced November 2011.

Comments: 6 pages, 3 figures, to appear in the proceedings of "Assembling the Puzzle of the Milky Way", Le Grand Bornand (Apr. 17-22, 2011), C. Reyle, A. Robin, M. Schultheis (eds.)

Weakening dark-matter cusps by clumpy baryonic infall

Authors: David Cole, Walter Dehnen, Mark Wilkinson

Abstract: 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 rem… ▽ More 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. △ Less

Submitted 20 May, 2011; originally announced May 2011.

Comments: 17 pages, 13 figures, accepted for publication in MNRAS

N-body simulations of gravitational dynamics

Authors: Walter Dehnen, Justin Read

Abstract: 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 imp… ▽ More 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. △ Less

Submitted 5 May, 2011; originally announced May 2011.

Comments: invited review (28 pages), to appear in European Physics Journal Plus

Journal ref: 2011, Eur. Phys. J. Plus 126, 55

Inviscid SPH

Authors: Lee Cullen, Walter Dehnen

Abstract: 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 subj… ▽ More 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. △ Less

Submitted 8 June, 2010; originally announced June 2010.

Comments: 14 pages (15 in arXiv), 15 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.408:669-683,2010

On the isolated dwarf galaxies: from cuspy to flat dark matter density profiles and metallicity gradients

Authors: S. Pasetto, E. K. Grebel, P. Berczik, R. Spurzem, W. Dehnen

Abstract: The chemodynamical evolution of spherical multi-component self-gravitating models for isolated dwarf galaxies is studied. We compare their evolution with and without feedback effects from star formation processes. We find that initially cuspy dark matter profiles flatten with time as a result of star formation, without any special tuning conditions. Thus the seemingly flattened profiles found in… ▽ More The chemodynamical evolution of spherical multi-component self-gravitating models for isolated dwarf galaxies is studied. We compare their evolution with and without feedback effects from star formation processes. We find that initially cuspy dark matter profiles flatten with time as a result of star formation, without any special tuning conditions. Thus the seemingly flattened profiles found in many dwarfs do not contradict the cuspy profiles predicted by cosmological models. We also calculate the chemical evolution of stars and gas, to permit comparisons with observational data. △ Less

Submitted 4 February, 2010; originally announced February 2010.

Comments: Accepted the 28/01/10

Local Kinematics and the Local Standard of Rest

Authors: R. Schoenrich, J. Binney, W. Dehnen

Abstract: We re-examine the stellar kinematics of the Solar neighbourhood in terms of the velocity of the Sun with respect to the local standard of rest. We show that the classical determination of its component V_sun in the direction of Galactic rotation via Stroemberg's relation is undermined by the metallicity gradient in the disc, which introduces a correlation between the colour of a group of stars a… ▽ More We re-examine the stellar kinematics of the Solar neighbourhood in terms of the velocity of the Sun with respect to the local standard of rest. We show that the classical determination of its component V_sun in the direction of Galactic rotation via Stroemberg's relation is undermined by the metallicity gradient in the disc, which introduces a correlation between the colour of a group of stars and the radial gradients of its properties. Comparing the local stellar kinematics to a chemodynamical model which accounts for these effects, we obtain (U,V,W)_sun = (11.1 +/- 0.74, 12.24 +/- 0.47, 7.25 +/-0.37) km/s, with additional systematic uncertainties of ~ (1,2,0.5) km/s. In particular, V_sun is 7 km/s larger than previously estimated. The new values of solar motion are extremely insensitive to the metallicity gradient within the disc. △ Less

Submitted 18 December, 2009; originally announced December 2009.

Comments: 5 pages, submitted to MNRAS

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

Authors: Walter Dehnen

Abstract: 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. Th… ▽ More 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 △ Less

Submitted 12 February, 2009; originally announced February 2009.

Comments: 9 pages, 6 figures, MNRAS, in press

Journal ref: Mon.Not.Roy.Astron.Soc.395:1079-1086,2009