1,228 research outputs found
Mapping the three-body system - decay time and reversibility
In this paper we carry out a quantitative analysis of the three-body systems
and map them as a function of decaying time and intial conguration, look at
this problem as an example of a simple deterministic system, and ask to what
extent the orbits are really predictable. We have investigated the behavior of
about 200 000 general Newtonian three body systems using the simplest initial
conditions. Within our resolution these cover all the possible states where the
objects are initially at rest and have no angular momentum. We have determined
the decay time-scales of the triple systems and show that the distribution of
this parameter is fractal in appearance. Some areas that appear stable on large
scales exhibit very narrow strips of instability and the overall pattern,
dominated by resonances, reminds us of a traditional Maasai warrior shield.
Also an attempt is made to recover the original starting conguration of the
three bodies by backward integration. We find there are instances where the
evolution to the future and to the past lead to different orbits, in spite of
time symmetric initial conditions. This implies that even in simple
deterministic systems there exists an Arrow of Time.Comment: 8 pages, 9 figures. Accepted for publication in MNRAS. Includes
low-resolution figures. High-resolution figures are available as PNG
A Hybrid N-Body Code Incorporating Algorithmic Regularization and Post-Newtonian Forces
We describe a novel N-body code designed for simulations of the central
regions of galaxies containing massive black holes. The code incorporates
Mikkola's 'algorithmic' chain regularization scheme including post-Newtonian
terms up to PN2.5 order. Stars moving beyond the chain are advanced using a
fourth-order integrator with forces computed on a GRAPE board. Performance
tests confirm that the hybrid code achieves better energy conservation, in less
elapsed time, than the standard scheme and that it reproduces the orbits of
stars tightly bound to the black hole with high precision. The hybrid code is
applied to two sample problems: the effect of finite-N gravitational
fluctuations on the orbits of the S-stars; and inspiral of an intermediate-mass
black hole into the galactic center.Comment: 12 pages, 15 figures, accepted for publication in MNRA
Asteroids in the Inner Solar System I - Existence
Ensembles of in-plane and inclined orbits in the vicinity of the Lagrange
points of the terrestrial planets are integrated for up to 100 million years.
The integrations incorporate the gravitational effects of Sun and the eight
planets (Pluto is neglected). Mercury is the least likely planet, as it is
unable to retain tadpole orbits over 100 million year timescales. Both Venus
and the Earth are much more promising, as they possess rich families of stable
tadpole and horseshoe orbits. Our survey of Trojans in the orbital plane of
Venus is undertaken for 25 million years. Some 40% of the survivors are on
tadpole orbits. For the Earth, the integrations are pursued for 50 million
years. The stable zones in the orbital plane are larger for the Earth than for
Venus, but fewer of the survivors are tadpoles. Both Venus and the Earth also
have regions in which inclined test particles can endure near the Lagrange
points. For Venus, only test particles close to the orbital plane are stable.
For the Earth, there are two bands of stability, one at low inclinations (i <
16 degrees) and one at moderate inclinations (between 24 degrees and 34
degrees). The inclined test particles that evade close encounters are primarily
moving on tadpole orbits. Our survey of in-plane test particles near the
Martian Lagrange points shows no survivors after 60 million years. Low
inclination test particles do not persist, as their inclinations are quickly
increased until the effects of a secular resonance with Jupiter cause
de-stabilisation. Numerical integrations of inclined test particles for
timespans of 25 million years show stable zones for inclinations between 14 and
40 degrees.Comment: 20 pages, 21 figures, Monthly Notices (in press
Chaos in the one-dimensional gravitational three-body problem
We have investigated the appearance of chaos in the 1-dimensional Newtonian
gravitational three-body system (three masses on a line with pairwise
potential). We have concentrated in particular on how the behavior changes when
the relative masses of the three bodies change (with negative total energy).
For two mass choices we have calculated 18000 full orbits (with initial states
on a lattice on the Poincar\'e section) and obtained dwell time
distributions. For 105 mass choices we have calculated Poincar\'e maps for
starting points. Our results show that the Poincar\'e section
(and hence the phase space) divides into three well defined regions with orbits
of different characteristics: 1) There is a region of fast scattering, with a
minimum of pairwise collisions and smooth dependence on initial values. 2) In
the chaotic scattering region the interaction times are longer, and both the
interaction time and the final state depend sensitively on the starting point
on the Poincar\'e section. For both 1) and 2) the initial and final states
consists of a binary + single particle. 3) The third region consists of
quasiperiodic orbits where the three masses are bound together forever. At the
center of the quasiperiodic region there is the periodic Schubart orbit, whose
stability turns out to correlate strongly with the global behavior.Comment: 24 pages of text (REVTEX 3.0) + 21 pages of figures. Figures are only
available in paper form, ask for a preprint from the author
Collisional dynamics around binary black holes in galactic centers
We follow the sinking of two massive black holes in a spherical stellar
system where the black holes become bound under the influence of dynamical
friction. Once bound, the binary hardens by three-body encounters with
surrounding stars. We find that the binary wanders inside the core, providing
an enhanced supply of reaction partners for the hardening. The binary evolves
into a highly eccentric orbit leading to coalescence well beyond a Hubble time.
These are the first results from a hybrid ``self consistent field'' (SCF) and
direct Aarseth N-body integrator (NBODY6), which combines the advantages of the
direct force calculation with the efficiency of the field method. The code is
designed for use on parallel architectures and is therefore applicable to
collisional N-body integrations with extraordinarily large particle numbers (>
10^5). This creates the possibility of simulating the dynamics of both globular
clusters with realistic collisional relaxation and stellar systems surrounding
supermassive black holes in galactic nuclei.Comment: 38 pages, 13 figures, submitted to ApJ, accepted, revised text and
added figure
Long-Term Evolution of Massive Black Hole Binaries. III. Binary Evolution in Collisional Nuclei
[Abridged] In galactic nuclei with sufficiently short relaxation times,
binary supermassive black holes can evolve beyond their stalling radii via
continued interaction with stars. We study this "collisional" evolutionary
regime using both fully self-consistent N-body integrations and approximate
Fokker-Planck models. The N-body integrations employ particle numbers up to
0.26M and a direct-summation potential solver; close interactions involving the
binary are treated using a new implementation of the Mikkola-Aarseth chain
regularization algorithm. Even at these large values of N, two-body scattering
occurs at high enough rates in the simulations that they can not be simply
scaled to the large-N regime of real galaxies. The Fokker-Planck model is used
to bridge this gap; it includes, for the first time, binary-induced changes in
the stellar density and potential. The Fokker-Planck model is shown to
accurately reproduce the results of the N-body integrations, and is then
extended to the much larger N regime of real galaxies. Analytic expressions are
derived that accurately reproduce the time dependence of the binary semi-major
axis as predicted by the Fokker-Planck model. Gravitational wave coalescence is
shown to occur in <10 Gyr in nuclei with velocity dispersions below about 80
km/s. Formation of a core results from a competition between ejection of stars
by the binary and re-supply of depleted orbits via two-body scattering. Mass
deficits as large as ~4 times the binary mass are produced before coalescence.
After the two black holes coalesce, a Bahcall-Wolf cusp appears around the
single hole in one relaxation time, resulting in a nuclear density profile
consisting of a flat core with an inner, compact cluster, similar to what is
observed at the centers of low-luminosity spheroids.Comment: 21 page
Dynamics of two planets in co-orbital motion
We study the stability regions and families of periodic orbits of two planets
locked in a co-orbital configuration. We consider different ratios of planetary
masses and orbital eccentricities, also we assume that both planets share the
same orbital plane. Initially we perform numerical simulations over a grid of
osculating initial conditions to map the regions of stable/chaotic motion and
identify equilibrium solutions. These results are later analyzed in more detail
using a semi-analytical model. Apart from the well known quasi-satellite (QS)
orbits and the classical equilibrium Lagrangian points L4 and L5, we also find
a new regime of asymmetric periodic solutions. For low eccentricities these are
located at , where \sigma is
the difference in mean longitudes and \Delta\omega is the difference in
longitudes of pericenter. The position of these Anti-Lagrangian solutions
changes with the mass ratio and the orbital eccentricities, and are found for
eccentricities as high as ~ 0.7. Finally, we also applied a slow mass variation
to one of the planets, and analyzed its effect on an initially asymmetric
periodic orbit. We found that the resonant solution is preserved as long as the
mass variation is adiabatic, with practically no change in the equilibrium
values of the angles.Comment: 9 pages, 11 figure
Interaction of massive black hole binaries with their stellar environment: II. Loss-cone depletion and binary orbital decay
We study the long-term evolution of massive black hole binaries (MBHBs) at
the centers of galaxies using detailed scattering experiments to solve the full
three-body problem. Ambient stars drawn from a isotropic Maxwellian
distribution unbound to the binary are ejected by the gravitational slingshot.
We construct a minimal, hybrid model for the depletion of the loss cone and the
orbital decay of the binary, and show that secondary slingshots - stars
returning on small impact parameter orbits to have a second super-elastic
scattering with the MBHB - may considerably help the shrinking of the pair in
the case of large binary mass ratios. In the absence of loss-cone refilling by
two-body relaxation or other processes, the mass ejected before the stalling of
a MBHB is half the binary reduced mass. About 50% of the ejected stars are
expelled ejected in a "burst" lasting ~1E4 yrs M_6^1/4, where M_6 is the binary
mass in units of 1E6 Msun. The loss cone is completely emptied in a few bulge
crossing timescales, 1E7 yrs M_6^1/4. Even in the absence of two-body
relaxation or gas dynamical processes, unequal mass and/or eccentric binaries
with M_6 >0.1 can shrink to the gravitational wave emission regime in less than
a Hubble time, and are therefore "safe" targets for the planned Laser
Interferometer Space Antenna (LISA).Comment: Minor revision. 10 pages, 7 figures, ApJ in pres
- âŠ