BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Walter Dehnen

Walter Dehnen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.EPastro-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

9 published item(s)

preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Secular Dynamics around a Supermassive Black Hole via Multipole Expansion

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Capture of interstellar objects I: the capture cross-section

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Capture of interstellar objects II: by the Solar system

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

Massive elliptical galaxies are typically observed to have central cores in their projected radial light profiles. Such cores have long been thought to form through `binary scouring&#39; 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&#39;, `tidal deposition&#39; and `gravitational wave (GW) induced recoil&#39;. 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&#39; 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&#39; observed in the core region of A2261-BCG are best explained as stalled perturbers resulting from minor mergers, though the brightest may also represent ejected SMBHs surrounded by a stellar cloak of bound stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-parametric Triaxial Deprojection of Elliptical Galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

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

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

0 citations0 reviewsNo code linkedOpen access