Long-term evolution of a planetesimal swarm in the vicinity of a protoplanet

Many models of planet formation involve scenarios in which one or a few large protoplanets interact with a swarm of much smaller planetesimals. In such scenarios, three-body perturbations by the protoplanet as well as mutual collisions and gravitational interactions between the swarm bodies are important in determining the velocity distribution of the swarm. We are developing a model to examine the effects of these processes on the evolution of a planetesimal swarm. The model consists of a combination of numerical integrations of the gravitational influence of one (or a few) massive protoplanets on swarm bodies together with a statistical treatment of the interactions between the planetesimals. Integrating the planetesimal orbits allows us to take into account effects that are difficult to model analytically or statistically, such as three-body collision cross-sections and resonant perturbations by the protoplanet, while using a statistical treatment for the particle-particle interactions allows us to use a large enough sample to obtain meaningful results

Collision probabilities in the presence of nebular gas drag

We are developing a model to determine what fraction of the planetesimals would have hit a protoplanet on their sunward journey as opposed to having a close approach and passing into an inferior orbit. The model involves direct numerical integration of restricted-three-body orbits using a predictor-corrector integrator. A simple gas drag law with a v(exp 2) dependence was also included in the equations of motion. Runs of 100 to 500 particles were already performed, while some future runs may require several times this number in order to get good impact statistics. All planetesimals start in superior orbits with semi-major axes 5 to 10 R(sub H) from the protoplanets, where R(sub H) is the protoplanet's Hill Sphere radius. The orbit is followed until the planetesimal passed into an inferior orbit at least 10 R(sub H) from the protoplanet. This process typically requires 10(exp 4) to 10(exp 5) orbits

Comparing HARPS and Kepler surveys: The alignment of multiple-planet systems

Aims. We study a subset of the planetary population characterized both by HARPS and Kepler surveys. We compare the statistical properties of planets in systems with m.sin i >5-10 M_Earth and R>2 R_Earth. If we assume that the underlying population has the same characteristics, the different detection sensitivity to the orbital inclination relative to the line of sight allows us to probe the planets' mutual inclination. Methods. We considered the frequency of systems with one, two and three planets as dictated by HARPS data. We used Kepler's planetary period and host mass and radii distributions (corrected from detection bias) to model planetary systems in a simple yet physically plausible way. We then varied the mutual inclination between planets in a system according to different prescriptions (completely aligned, Rayleigh distributions and isotropic) and compared the transit frequencies with one, two or three planets with those measured by Kepler. Results. The results show that the two datasets are compatible, a remarkable result especially because there are no tunable knobs other than the assumed inclination distribution. For m.sin i cutoffs of 7-10 M_Earth, which are those expected to correspond to the radius cutoff of 2 R_Earth, we conclude that the results are better described by a Rayleigh distribution with mode of 1 deg or smaller. We show that the best-fit scenario only becomes a Rayleigh distribution with mode of 5 deg if we assume a rather extreme mass-radius relationship for the planetary population. Conclusions. These results have important consequences for our understanding of the role of several proposed formation and evolution mechanisms. They confirm that planets are likely to have been formed in a disk and show that most planetary systems evolve quietly without strong angular momentum exchanges (abridged).Comment: 10 pages, 6 figures, 4 tables, accepted for publication in Astronomy & Astrophysic

Detection and measurement of planetary systems with GAIA

We use detailed numerical simulations and the $\upsilon$ Andromedae, planetary system as a template to evaluate the capability of the ESA Cornerstone Mission GAIA in detecting and measuring multiple planets around solar-type stars in the neighborhood of the Solar System. For the outer two planets of the $\upsilon$ Andromedae, system, GAIA high-precision global astrometric measurements would provide estimates of the full set of orbital elements and masses accurate to better than 1--10%, and would be capable of addressing the coplanarity issue by determining the true geometry of the system with uncertainties of order of a few degrees. Finally, we discuss the generalization to a variety of configurations of potential planetary systems in the solar neighborhood for which GAIA could provide accurate measurements of unique value for the science of extra-solar planets.Comment: 4 pages, 2 pictures, accepted for publication in A&A Letter

The Exoplanet Eccentricity Distribution from Kepler Planet Candidates

The eccentricity distribution of exoplanets is known from radial velocity surveys to be divergent from circular orbits beyond 0.1 AU. This is particularly the case for large planets where the radial velocity technique is most sensitive. The eccentricity of planetary orbits can have a large effect on the transit probability and subsequently the planet yield of transit surveys. The Kepler mission is the first transit survey that probes deep enough into period-space to allow this effect to be seen via the variation in transit durations. We use the Kepler planet candidates to show that the eccentricity distribution is consistent with that found from radial velocity surveys to a high degree of confidence. We further show that the mean eccentricity of the Kepler candidates decreases with decreasing planet size indicating that smaller planets are preferentially found in low-eccentricity orbits.Comment: 6 pages, 4 figures, accepted for publication in MNRA

Transit Timing Observations from Kepler: VII. Confirmation of 27 planets in 13 multiplanet systems via Transit Timing Variations and orbital stability

We confirm 27 planets in 13 planetary systems by showing the existence of statistically significant anti-correlated transit timing variations (TTVs), which demonstrates that the planet candidates are in the same system, and long-term dynamical stability, which places limits on the masses of the candidates---showing that they are planetary. %This overall method of planet confirmation was first applied to \kepler systems 23 through 32. All of these newly confirmed planetary systems have orbital periods that place them near first-order mean motion resonances (MMRs), including 6 systems near the 2:1 MMR, 5 near 3:2, and one each near 4:3, 5:4, and 6:5. In addition, several unconfirmed planet candidates exist in some systems (that cannot be confirmed with this method at this time). A few of these candidates would also be near first order MMRs with either the confirmed planets or with other candidates. One system of particular interest, Kepler-56 (KOI-1241), is a pair of planets orbiting a 12th magnitude, giant star with radius over three times that of the Sun and effective temperature of 4900 K---among the largest stars known to host a transiting exoplanetary system.Comment: 12 pages, 13 figures, 5 tables. Submitted to MNRA

ϕ\phi- meson Production at RHIC energies using the PHENIX Detector

Light vector mesons are among the most informative probes to understand the strongly coupled Quark Gluon Plasma created at RHIC. The suppression of light mesons at high transverse momentum, compared to expectations from scaled $p+p$ results, reflects the properties of the strongly interacting matter formed. The $\phi$-meson is one of the probes whose systematic measurement in $p+p$, $d+Au$ and $Au+Au$ collisions can provide useful information about initial and final state effects on particle production. The mass, width and branching ratio of the $\phi$-meson decay in the di-kaon and di-electron decay channels could be modified in \au collisions due to the restoration of chiral symmetry in the QGP. The PHENIX experiment at RHIC has measured $\phi$-meson production in various systems ranging form $p+p$, $d+Au$ to $Au+Au$ collisions via both its di-electron and di-kaon decay modes. A summary of PHENIX results on invariant spectra, nuclear modification factor and elliptic flow of the $\phi$-meson are presented here

A Planetary Companion to the Nearby M4 Dwarf, Gliese 876

Doppler measurements of the M4 dwarf star, Gliese 876, taken at both Lick and Keck Observatory reveal periodic, Keplerian velocity variations with a period of 61 days. The orbital fit implies that the companion has a mass of, M = 2.1 MJUP /sin i, an orbital eccentricity of, e = 0.27+-0.03, and a semimajor axis of, a = 0.21 AU. The planet is the first found around an M dwarf, and was drawn from a survey of 24 such stars at Lick Observatory. It is the closest extrasolar planet yet found, providing opportunities for follow--up detection. The presence of a giant planet on a non-circular orbit, 0.2 AU from a 1/3 M_Sun star, presents a challenge to planet formation theory. This planet detection around an M dwarf suggests that giant planets are numerous in the Galaxy.Comment: 13 pages, 3 Figure

Last giant impact on the Neptunian system. Constraints on oligarchic masses in the trans-Saturnian region

Stochastic impacts by large bodies are, at present, the usually accepted mechanisms able to account for the obliquity of the ice giants. We attempt to set constraints on giant impacts as the cause of Neptune's current obliquity in the framework of modern theories. We also use the present orbital properties of the Neptunian irregular satellites (with the exception of Triton) to set constraints on the scenario of giant impacts at the end of Neptune formation. We model the angular momentum transfer to proto-Neptune and the impulse transfer to its irregular satellites by the last stochastic collision (GC) between the protoplanet and an oligarchic mass at the end of Neptune's formation. We obtain that an impactor mass greather than 4 Earth masses is not possible since it cannot reproduce the present rotational properties of the planet, unless the impact parameter of the collision were very small. On the other hand, if the impactor mass was greather than 1.4 Earth masses, the present Neptunian irregular satellites had to be formed or captured after the end of stochastic impacts. The upper bounds on the oligarchic masses (4 Earth masses from the obliquity of Neptune and 1.4 earth masses from the Neptunian irregular satellites) are independent of unknown parameters, such as the mass and distribution of the planetesimals, the location at which Uranus and Neptune were formed, the Solar Nebula initial surface mass density, and the growth regime. If stochastic impacts had occurred, these results should be understood as upper constraints on the oligarchic masses in the trans-Saturnian region at the end of ice planet formation and may be used to set constraints on planetary formation scenarios.Comment: Paper accepted for publication in Astronomy & Astrophysic