Researcher profile

S. J. Aarseth

S. J. Aarseth contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
6works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

6 published item(s)

preprint2016arXiv

Dynamical impact of the Planet Nine scenario: N-body experiments

The Planet Nine hypothesis has now enough constraints to deserve further attention in the form of detailed numerical experiments. The results of such studies can help us improve our understanding of the dynamical effects of such a hypothetical object on the extreme trans-Neptunian objects or ETNOs and perhaps provide additional constraints on the orbit of Planet Nine itself. Here, we present the results of direct N-body calculations including the latest data available on the Planet Nine conjecture. The present-day orbits of the six ETNOs originally linked to the hypothesis are evolved backwards in time and into the future under some plausible incarnations of the hypothesis to investigate if the values of several orbital elements, including the argument of perihelion, remain confined to relatively narrow ranges. We find that a nominal Planet Nine can keep the orbits of (90377) Sedna and 2012 VP113 relatively well confined in orbital parameter space for hundreds of Myr, but it may make the orbits of 2004 VN112, 2007 TG422 and 2013 RF98 very unstable on time-scales of dozens of Myr, turning them retrograde and eventually triggering their ejection from the Solar system. Far more stable orbital evolution is found with slightly modified orbits for Planet Nine.

preprint2015arXiv

Chasing the Chelyabinsk asteroid N-body style

On 2013 February 15 a small asteroid rammed against the atmosphere above the region of Chelyabinsk in Russia, producing the most powerful superbolide since the Tunguska event in 1908. Lacking proper astrometric observations, the pre-impact orbit of this object has been determined using videos, satellite images, and pure geometry. Unfortunately, more than two years after the event, the published estimates vary so much that there is no clear orbital solution that could be used to investigate the origin of the impactor and the existence of dynamically, or perhaps even genetically, related asteroids. Here, we revisit this topic using a full N-body approach. A robust statistical test is applied to published solutions to discard those unable to produce a virtual impact at the observed time (03:20:20.8 s UTC). The same N-body methodology and the latest ephemerides are used to compute a new orbital solution: a=1.6247 AU, e=0.5318, i=3.9750 degrees, Omega=326.4607 degrees and omega=109.7012 degrees. This new solution --which has an impact probability > 0.99999 and uncertainties in time and space of 0.2 s and 6 km, respectively-- is utilized to explore the past orbital evolution of the impactor as well as the presence of near-Earth objects moving in similar paths. A dynamical link between asteroid 2011 EO40 and the Chelyabinsk impactor is confirmed. Alternative orbital solutions are extensively explored.

preprint2015arXiv

Dynamical evolution of massive black holes in galactic-scale N-body simulations - introducing the regularized tree code "rVINE"

We present a hybrid code combining the OpenMP-parallel tree code VINE with an algorithmic chain regularization scheme. The new code, called "rVINE", aims to significantly improve the accuracy of close encounters of massive bodies with supermassive black holes in galaxy-scale numerical simulations. We demonstrate the capabilities of the code by studying two test problems, the sinking of a single massive black hole to the centre of a gas-free galaxy due to dynamical friction and the hardening of a supermassive black hole binary due to close stellar encounters. We show that results obtained with rVINE compare well with NBODY7 for problems with particle numbers that can be simulated with NBODY7. In particular, in both NBODY7 and rVINE we find a clear N-dependence of the binary hardening rate, a low binary eccentricity and moderate eccentricity evolution, as well as the conversion of the galaxy's inner density profile from a cusp to a a core via the ejection of stars at high velocity. The much larger number of particles that can be handled by rVINE will open up exciting opportunities to model stellar dynamics close to SMBHs much more accurately in a realistic galactic context. This will help to remedy the inherent limitations of commonly used tree solvers to follow the correct dynamical evolution of black holes in galaxy scale simulations.

preprint2014arXiv

Flipping minor bodies: what comet 96P/Machholz 1 can tell us about the orbital evolution of extreme trans-Neptunian objects and the production of near-Earth objects on retrograde orbits

Nearly all known extreme trans-Neptunian objects (ETNOs) have argument of perihelion close to 0 degrees. An existing observational bias strongly favours the detection of ETNOs with arguments of perihelion close to 0 degrees and 180 degrees yet no objects have been found at 180 degrees. No plausible explanation has been offered so far to account for this unusual pattern. Here, we study the dynamical evolution of comet 96P/Machholz 1, a bizarre near-Earth object (NEO) that may provide the key to explain the puzzling clustering of orbits around argument of perihelion close to 0 degrees recently found for the population of ETNOs. Comet 96P/Machholz 1 is currently locked in a Kozai resonance with Jupiter such that the value of its argument of perihelion is always close to 0 degrees at its shortest possible perihelion (highest eccentricity and lowest inclination) and about 180 degrees near its shortest aphelion (longest perihelion distance, lowest eccentricity and highest inclination). If this object is a dynamical analogue (albeit limited) of the known ETNOs, this implies that massive perturbers must keep them confined in orbital parameter space. Besides, its future dynamical evolution displays orbital flips when its eccentricity is excited to a high value and its orbit turns over by nearly 180 degrees, rolling over its major axis. This unusual behaviour, that is preserved when post-Newtonian terms are included in the numerical integrations, may also help understand the production of NEOs on retrograde orbits.

preprint1998arXiv

Basic N-Body Modelling of the Evolution of Globular Clusters. I. Time Scaling

We consider the use of N-body simulations for studying the evolution of rich star clusters (i.e. globular clusters). The dynamical processes included in this study are restricted to gravitational (point-mass) interactions, the steady tidal field of a galaxy, and instantaneous mass loss resulting from stellar evolution. With evolution driven by these mechanisms, it is known that clusters fall roughly into two broad classes: those that dissipate promptly in the tidal field, as a result of mass loss, and those that survive long enough for their evolution to become dominated by two-body relaxation. The time scales of the processes we consider scale in different ways with the number of stars in the simulation, and the main aim of the paper is to suggest how the scaling of a simulation should be done so that the results are representative of the evolution of a `real' cluster. We investigate three different ways of scaling time. One of these is appropriate to the first type of cluster, i.e. those that dissipate rapidly, and similarly a second scaling is appropriate only to the second (relaxation-dominated) type. We also develop a hybrid scaling which is a satisfactory compromise for both types of cluster. Finally we present evidence that the widely used Fokker-Planck method produces models which are in good agreement with N-body models of those clusters which are relaxation dominated, at least for N-body models with several thousand particles, but that the Fokker-Planck models evolve too fast for clusters which dissipate promptly.

preprint1997arXiv

Hierarchical Systems in Open Clusters

In this paper we study the formation, evolution and disruption of hierarchical systems in open clusters. With this purpose, N-body simulations of star clusters containing an initial population of binaries have been carried out using Aarseth's NBODY4 and NBODY5 codes. Stable triples may form from strong interactions of two binaries in which the widest pair is disrupted. The most frequent type of hierarchical systems found in the cluster models are triples in which the outer star is single, but in some cases the outer body is also a binary, giving a hierarchical quadruple. Many triple systems are non-coplanar and the presence of even a very distant and small outer companion may affect the orbital parameters of the inner binary, including a possible mechanism of significant shrinkage if the binary experiences a weak tidal dissipation. The main features of these systems are analyzed in order to derive general properties which can be checked by observations. The formation of exotic objects, such as blue stragglers and white dwarf binaries, inside hierarchical triple systems is particularly interesting.