Circumbinary Chaos: Using Pluto's Newest Moon to Constrain the Masses of Nix & Hydra

The Pluto system provides a unique local laboratory for the study of binaries with multiple low mass companions. In this paper, we study the orbital stability of P4, the most recently discovered moon in the Pluto system. This newfound companion orbits near the plane of the Pluto-Charon binary, roughly halfway between the two minor moons Nix and Hydra. We use a suite of few body integrations to constrain the masses of Nix and Hydra, and the orbital parameters of P4. For the system to remain stable over the age of the Solar System, the masses of Nix and Hydra likely do not exceed 5e16 kg and 9e16 kg, respectively. These upper limits assume a fixed mass ratio between Nix and Hydra at the value implied by their median optical brightness. Our study finds that stability is more sensitive to their total mass and that a downward revision of Charon's eccentricity (from our adopted value of 0.0035) is unlikely to significantly affect our conclusions. Our upper limits are an order of magnitude below existing astrometric limits on the masses of Nix and Hydra. For a density at least that of ice, the albedos of Nix and Hydra would exceed 0.3. This constraint implies they are icy, as predicted by giant impact models. Even with these low masses, P4 only remains stable if its eccentricity e < 0.02. The 5:1 commensurability with Charon is particularly unstable, Combining stability constraints with the observed mean motion places the preferred orbit for P4 just exterior to the 5:1 resonance. These predictions will be tested when the New Horizons satellite visits Pluto. Based on the results for the Pluto-Charon system, we expect that circumbinary, multi-planet systems will be more widely spaced than their singleton counterparts. Further, circumbinary exoplanets close to the three-body stability boundary, such as those found by Kepler, are less likely to have other companions nearby.Comment: 12 pages, 7 figure

The Intrinsic Size of Sagittarius A* from 0.35 cm to 6 cm

We present new high-resolution observations of Sagittarius A* at wavelengths of 17.4 to 23.8 cm with the Very Large Array in A configuration with the Pie Town Very Long Baseline Array antenna. We use the measured sizes to calibrate the interstellar scattering law and find that the major axis size of the scattering law is smaller by ~6% than previous estimates. Using the new scattering law, we are able to determine the intrinsic size of Sgr A* at wavelengths from 0.35 cm to 6 cm using existing results from the VLBA. The new law increases the intrinsic size at 0.7 cm by ~20% and <5% at 0.35 cm. The intrinsic size is 13^{+7}_{-3} Schwarzschild radii at 0.35 cm and is proportional to lambda^gamma, where gamma is in the range 1.3 to 1.7.Comment: ApJL, in pres

Particle Stirring in Turbulent Gas Disks: Including Orbital Oscillations

We describe the diffusion and random velocities of solid particles due to stochastic forcing by turbulent gas. We include the orbital dynamics of Keplerian disks, both in-plane epicycles and vertical oscillations. We obtain a new result for the diffusion of solids. The Schmidt number (ratio of gas to particle diffusivity) is Sc = 1 + (Omega t_stop)^2, in terms of the particle stopping time, t_stop, and the orbital frequency, Omega. The standard result, Sc = 1 + t_stop/t_eddy, in terms of the eddy turnover time, t_eddy, is shown to be incorrect. The main difference is that Sc rises quadratically, not linearly, with stopping time. Consequently, particles larger than ~ 10 cm in protoplanetary disks will suffer less radial diffusion and will settle closer to the midplane. Such a layer of boulders would be more prone to gravitational collapse. Our predictions of RMS speeds, vertical scale height and diffusion coefficients will help interpret numerical simulations. We confirm previous results for the vertical stirring of particles (scale heights and random velocities), and add a correction for arbitrary ratios of eddy to orbital times. The particle layer becomes thinner for t_eddy > 1/Omega, with the strength of turbulent diffusion held fixed. We use two analytic techniques -- the Hinze-Tchen formalism and the Fokker-Planck equation with velocity diffusion -- with identical results when the regimes of validity overlap. We include simple physical arguments for the scaling of our results.Comment: 17 pages, 7 figures, 2 tables, accepted to Icaru

Interesting dynamics at high mutual inclination in the framework of the Kozai problem with an eccentric perturber

We study the dynamics of the 3-D three-body problem of a small body moving under the attractions of a star and a giant planet which orbits the star on a much wider and elliptic orbit. In particular, we focus on the influence of an eccentric orbit of the outer perturber on the dynamics of a small highly inclined inner body. Our analytical study of the secular perturbations relies on the classical octupole hamiltonian expansion (third-order theory in the ratio of the semi-major axes), as third-order terms are needed to consider the secular variations of the outer perturber and potential secular resonances between the arguments of the pericenter and/or longitudes of the node of both bodies. Short-period averaging and node reduction (Laplace plane) reduce the problem to two degrees of freedom. The four-dimensional dynamics is analyzed through representative planes which identify the main equilibria of the problem. As in the circular problem (i.e. perturber on a circular orbit), the "Kozai-bifurcated" equilibria play a major role in the dynamics of an inner body on quasi-circular orbit: its eccentricity variations are very limited for mutual inclination between the orbital planes smaller than ~40^{\deg}, while they become large and chaotic for higher mutual inclination. Particular attention is also given to a region around 35^{\deg} of mutual inclination, detected numerically by Funk et al. (2011) and consisting of long-time stable and particularly low eccentric orbits of the small body. Using a 12th-order Hamiltonian expansion in eccentricities and inclinations, in particular its action-angle formulation obtained by Lie transforms in Libert & Henrard (2008), we show that this region presents an equality of two fundamental frequencies and can be regarded as a secular resonance. Our results also apply to binary star systems where a planet is revolving around one of the two stars.Comment: 12 pages, 9 figures, accepted for publication in MNRA

On the Formation of Planetesimals via Secular Gravitational Instabilities with Turbulent Stirring

We study the gravitational instability (GI) of small solids in a gas disk as a mechanism to form planetesimals. Dissipation from gas drag introduces secular GI, which proceeds even when standard GI criteria for a critical density or Toomre's $Q$ predict stability. We include the stabilizing effects of turbulent diffusion, which suppresses small scale GI. The radially wide rings that do collapse contain up to $\sim 0.1$ Earth masses of solids. Subsequent fragmentation of the ring (not modeled here) would produce a clan of chemically homogenous planetesimals. Particle radial drift time scales (and, to a lesser extent, disk lifetimes and sizes) restrict the viability of secular GI to disks with weak turbulent diffusion, characterized by $\alpha \lesssim 10^{-4}$. Thus midplane dead zones are a preferred environment. Large solids with radii $\gtrsim 10$ cm collapse most rapidly because they partially decouple from the gas disk. Smaller solids, even below $\sim$ mm-sizes could collapse if particle-driven turbulence is weakened by either localized pressure maxima or super-Solar metallicity. Comparison with simulations that include particle clumping by the streaming instability shows that our linear model underpredicts rapid, small scale gravitational collapse. Thus the inclusion of more detailed gas dynamics promotes the formation of planetesimals. We discuss relevant constraints from Solar System and accretion disk observations.Comment: Accepted for publication in the Astrophysical Journal; 20 pages, 10 figure

Parallel electric field generation by Alfven wave turbulence

{This work aims to investigate the spectral structure of the parallel electric field generated by strong anisotropic and balanced Alfvenic turbulence in relation with the problem of electron acceleration from the thermal population in solar flare plasma conditions.} {We consider anisotropic Alfvenic fluctuations in the presence of a strong background magnetic field. Exploiting this anisotropy, a set of reduced equations governing non-linear, two-fluid plasma dynamics is derived. The low-$\beta$ limit of this model is used to follow the turbulent cascade of the energy resulting from the non-linear interaction between kinetic Alfven waves, from the large magnetohydrodynamics (MHD) scales with $k_{\perp}\rho_{s}\ll 1$ down to the small "kinetic" scales with $k_{\perp}\rho_{s} \gg 1$, $\rho_{s}$ being the ion sound gyroradius.} {Scaling relations are obtained for the magnitude of the turbulent electromagnetic fluctuations, as a function of $k_{\perp}$ and $k_{\parallel}$, showing that the electric field develops a component parallel to the magnetic field at large MHD scales.} {The spectrum we derive for the parallel electric field fluctuations can be effectively used to model stochastic resonant acceleration and heating of electrons by Alfven waves in solar flare plasma conditions

Adding particle collisions to the formation of asteroids and Kuiper belt objects via streaming instabilities

Modelling the formation of super-km-sized planetesimals by gravitational collapse of regions overdense in small particles requires numerical algorithms capable of handling simultaneously hydrodynamics, particle dynamics and particle collisions. While the initial phases of radial contraction are dictated by drag forces and gravity, particle collisions become gradually more significant as filaments contract beyond Roche density. Here we present a new numerical algorithm for treating momentum and energy exchange in collisions between numerical superparticles representing a high number of physical particles. We adopt a Monte Carlo approach where superparticle pairs in a grid cell collide statistically on the physical collision time-scale. Collisions occur by enlarging particles until they touch and solving for the collision outcome, accounting for energy dissipation in inelastic collisions. We demonstrate that superparticle collisions can be consistently implemented at a modest computational cost. In protoplanetary disc turbulence driven by the streaming instability, we argue that the relative Keplerian shear velocity should be subtracted during the collision calculation. If it is not subtracted, density inhomogeneities are too rapidly diffused away, as bloated particles exaggerate collision speeds. Local particle densities reach several thousand times the mid-plane gas density. We find efficient formation of gravitationally bound clumps, with a range of masses corresponding to contracted radii from 100 to 400 km when applied to the asteroid belt and 150 to 730 km when applied to the Kuiper belt, extrapolated using a constant self-gravity parameter. The smaller planetesimals are not observed at low resolution, but the masses of the largest planetesimals are relatively independent of resolution and treatment of collisions.Comment: Version accepted for publication in A&

Astrophysical turbulence modeling

The role of turbulence in various astrophysical settings is reviewed. Among the differences to laboratory and atmospheric turbulence we highlight the ubiquitous presence of magnetic fields that are generally produced and maintained by dynamo action. The extreme temperature and density contrasts and stratifications are emphasized in connection with turbulence in the interstellar medium and in stars with outer convection zones, respectively. In many cases turbulence plays an essential role in facilitating enhanced transport of mass, momentum, energy, and magnetic fields in terms of the corresponding coarse-grained mean fields. Those transport properties are usually strongly modified by anisotropies and often completely new effects emerge in such a description that have no correspondence in terms of the original (non coarse-grained) fields.Comment: 88 pages, 26 figures, published in Reports on Progress in Physic

Against all odds? Forming the planet of the HD196885 binary

HD196885Ab is the most "extreme" planet-in-a-binary discovered to date, whose orbit places it at the limit for orbital stability. The presence of a planet in such a highly perturbed region poses a clear challenge to planet-formation scenarios. We investigate this issue by focusing on the planet-formation stage that is arguably the most sensitive to binary perturbations: the mutual accretion of kilometre-sized planetesimals. To this effect we numerically estimate the impact velocities $dv$ amongst a population of circumprimary planetesimals. We find that most of the circumprimary disc is strongly hostile to planetesimal accretion, especially the region around 2.6AU (the planet's location) where binary perturbations induce planetesimal-shattering $dv$ of more than 1km/s. Possible solutions to the paradox of having a planet in such accretion-hostile regions are 1) that initial planetesimals were very big, at least 250km, 2) that the binary had an initial orbit at least twice the present one, and was later compacted due to early stellar encounters, 3) that planetesimals did not grow by mutual impacts but by sweeping of dust (the "snowball" growth mode identified by Xie et al., 2010b), or 4) that HD196885Ab was formed not by core-accretion but by the concurent disc instability mechanism. All of these 4 scenarios remain however highly conjectural.Comment: accepted for publication by Celestial Mechanics and Dynamical Astronomy (Special issue on EXOPLANETS