Hydrogen-Helium Mixtures at High Pressure

The properties of hydrogen-helium mixtures at high pressure are crucial to address important questions about the interior of Giant planets e.g. whether Jupiter has a rocky core and did it emerge via core accretion? Using path integral Monte Carlo simulations, we study the properties of these mixtures as a function of temperature, density and composition. The equation of state is calculated and compared to chemical models. We probe the accuracy of the ideal mixing approximation commonly used in such models. Finally, we discuss the structure of the liquid in terms of pair correlation functions.Comment: Proceedings article of the 5th Conference on Cryocrystals and Quantum Crystals in Wroclaw, Poland, submitted to J. Low. Temp. Phys. (2004

The distribution of transit durations for Kepler planet candidates and implications for their orbital eccentricities

‘In these times, during the rise in the popularity of institutional repositories, the Society does not forbid authors from depositing their work in such repositories. However, the AAS regards the deposit of scholarly work in such repositories to be a decision of the individual scholar, as long as the individual's actions respect the diligence of the journals and their reviewers.’ Original article can be found at : http://iopscience.iop.org/ Copyright American Astronomical SocietyDoppler planet searches have discovered that giant planets follow orbits with a wide range of orbital eccentricities, revolutionizing theories of planet formation. The discovery of hundreds of exoplanet candidates by NASA's Kepler mission enables astronomers to characterize the eccentricity distribution of small exoplanets. Measuring the eccentricity of individual planets is only practical in favorable cases that are amenable to complementary techniques (e.g., radial velocities, transit timing variations, occultation photometry). Yet even in the absence of individual eccentricities, it is possible to study the distribution of eccentricities based on the distribution of transit durations (relative to the maximum transit duration for a circular orbit). We analyze the transit duration distribution of Kepler planet candidates. We find that for host stars with T > 5100 K we cannot invert this to infer the eccentricity distribution at this time due to uncertainties and possible systematics in the host star densities. With this limitation in mind, we compare the observed transit duration distribution with models to rule out extreme distributions. If we assume a Rayleigh eccentricity distribution for Kepler planet candidates, then we find best fits with a mean eccentricity of 0.1-0.25 for host stars with T ≤ 5100 K. We compare the transit duration distribution for different subsets of Kepler planet candidates and discuss tentative trends with planetary radius and multiplicity. High-precision spectroscopic follow-up observations for a large sample of host stars will be required to confirm which trends are real and which are the results of systematic errors in stellar radii. Finally, we identify planet candidates that must be eccentric or have a significantly underestimated stellar radius.Peer reviewedFinal Accepted Versio

Orbital structure of the GJ876 extrasolar planetary system, based on the latest Keck and HARPS radial velocity data

We use full available array of radial velocity data, including recently published HARPS and Keck observatory sets, to characterize the orbital configuration of the planetary system orbiting GJ876. First, we propose and describe in detail a fast method to fit perturbed orbital configuration, based on the integration of the sensitivity equations inferred by the equations of the original $N$-body problem. Further, we find that it is unsatisfactory to treat the available radial velocity data for GJ876 in the traditional white noise model, because the actual noise appears autocorrelated (and demonstrates non-white frequency spectrum). The time scale of this correlation is about a few days, and the contribution of the correlated noise is about 2 m/s (i.e., similar to the level of internal errors in the Keck data). We propose a variation of the maximum-likelihood algorithm to estimate the orbital configuration of the system, taking into account the red noise effects. We show, in particular, that the non-zero orbital eccentricity of the innermost planet \emph{d}, obtained in previous studies, is likely a result of misinterpreted red noise in the data. In addition to offsets in some orbital parameters, the red noise also makes the fit uncertainties systematically underestimated (while they are treated in the traditional white noise model). Also, we show that the orbital eccentricity of the outermost planet is actually ill-determined, although bounded by $\sim 0.2$. Finally, we investigate possible orbital non-coplanarity of the system, and limit the mutual inclination between the planets \emph{b} and \emph{c} orbits by $5^\circ-15^\circ$, depending on the angular position of the mutual orbital nodes.Comment: 36 pages, 11 figures, 3 tables; Accepted to Celestial Mechanics and Dynamical Astronom

Search for Small Trans-Neptunian Objects by the TAOS Project

The Taiwan-America Occultation Survey (TAOS) aims to determine the number of small icy bodies in the outer reach of the Solar System by means of stellar occultation. An array of 4 robotic small (D=0.5 m), wide-field (f/1.9) telescopes have been installed at Lulin Observatory in Taiwan to simultaneously monitor some thousand of stars for such rare occultation events. Because a typical occultation event by a TNO a few km across will last for only a fraction of a second, fast photometry is necessary. A special CCD readout scheme has been devised to allow for stellar photometry taken a few times per second. Effective analysis pipelines have been developed to process stellar light curves and to correlate any possible flux changes among all telescopes. A few billion photometric measurements have been collected since the routine survey began in early 2005. Our preliminary result of a very low detection rate suggests a deficit of small TNOs down to a few km size, consistent with the extrapolation of some recent studies of larger (30--100 km) TNOs.Comment: 4 pages, 3 figures, IAU Symposium 23

Origin and Evolution of Saturn's Ring System

The origin and long-term evolution of Saturn's rings is still an unsolved problem in modern planetary science. In this chapter we review the current state of our knowledge on this long-standing question for the main rings (A, Cassini Division, B, C), the F Ring, and the diffuse rings (E and G). During the Voyager era, models of evolutionary processes affecting the rings on long time scales (erosion, viscous spreading, accretion, ballistic transport, etc.) had suggested that Saturn's rings are not older than 100 My. In addition, Saturn's large system of diffuse rings has been thought to be the result of material loss from one or more of Saturn's satellites. In the Cassini era, high spatial and spectral resolution data have allowed progress to be made on some of these questions. Discoveries such as the ''propellers'' in the A ring, the shape of ring-embedded moonlets, the clumps in the F Ring, and Enceladus' plume provide new constraints on evolutionary processes in Saturn's rings. At the same time, advances in numerical simulations over the last 20 years have opened the way to realistic models of the rings's fine scale structure, and progress in our understanding of the formation of the Solar System provides a better-defined historical context in which to understand ring formation. All these elements have important implications for the origin and long-term evolution of Saturn's rings. They strengthen the idea that Saturn's rings are very dynamical and rapidly evolving, while new arguments suggest that the rings could be older than previously believed, provided that they are regularly renewed. Key evolutionary processes, timescales and possible scenarios for the rings's origin are reviewed in the light of tComment: Chapter 17 of the book ''Saturn After Cassini-Huygens'' Saturn from Cassini-Huygens, Dougherty, M.K.; Esposito, L.W.; Krimigis, S.M. (Ed.) (2009) 537-57

Collisional Velocities and Rates in Resonant Planetesimal Belts

We consider a belt of small bodies around a star, captured in one of the external or 1:1 mean-motion resonances with a massive perturber. The objects in the belt collide with each other. Combining methods of celestial mechanics and statistical physics, we calculate mean collisional velocities and collisional rates, averaged over the belt. The results are compared to collisional velocities and rates in a similar, but non-resonant belt, as predicted by the particle-in-a-box method. It is found that the effect of the resonant lock on the velocities is rather small, while on the rates more substantial. The collisional rates between objects in an external resonance are by about a factor of two higher than those in a similar belt of objects not locked in a resonance. For Trojans under the same conditions, the collisional rates may be enhanced by up to an order of magnitude. Our results imply, in particular, shorter collisional lifetimes of resonant Kuiper belt objects in the solar system and higher efficiency of dust production by resonant planetesimals in debris disks around other stars.Comment: 31 pages, 11 figures (some of them heavily compressed to fit into arxiv-maximum filesize), accepted for publication at "Celestial Mechanics and Dynamical Astronomy

Forced oscillations in a hydrodynamical accretion disk and QPOs

This is the second of a series of papers aimed to look for an explanation on the generation of high frequency quasi-periodic oscillations (QPOs) in accretion disks around neutron star, black hole, and white dwarf binaries. The model is inspired by the general idea of a resonance mechanism in the accretion disk oscillations as was already pointed out by Abramowicz & Klu{\'z}niak (\cite{Abramowicz2001}). In a first paper (P\'etri \cite{Petri2005a}, paper I), we showed that a rotating misaligned magnetic field of a neutron star gives rise to some resonances close to the inner edge of the accretion disk. In this second paper, we suggest that this process does also exist for an asymmetry in the gravitational potential of the compact object. We prove that the same physics applies, at least in the linear stage of the response to the disturbance in the system. This kind of asymmetry is well suited for neutron stars or white dwarfs possessing an inhomogeneous interior allowing for a deviation from a perfectly spherically symmetric gravitational field. We show by a linear analysis that the disk initially in a cylindrically symmetric stationary state is subject to three kinds of resonances: a corotation resonance, a Lindblad resonance due to a driven force and a parametric sonance. The highest kHz QPOs are then interpreted as the orbital frequency of the disk at locations where the response to the resonances are maximal. It is also found that strong gravity is not required to excite the resonances.Comment: Accepte

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

Kepler-22b: A 2.4 Earth-radius Planet in the Habitable Zone of a Sun-like Star

A search of the time-series photometry from NASA's Kepler spacecraft reveals a transiting planet candidate orbiting the 11th magnitude G5 dwarf KIC 10593626 with a period of 290 days. The characteristics of the host star are well constrained by high-resolution spectroscopy combined with an asteroseismic analysis of the Kepler photometry, leading to an estimated mass and radius of 0.970 +/- 0.060 MSun and 0.979 +/- 0.020 RSun. The depth of 492 +/- 10ppm for the three observed transits yields a radius of 2.38 +/- 0.13 REarth for the planet. The system passes a battery of tests for false positives, including reconnaissance spectroscopy, high-resolution imaging, and centroid motion. A full BLENDER analysis provides further validation of the planet interpretation by showing that contamination of the target by an eclipsing system would rarely mimic the observed shape of the transits. The final validation of the planet is provided by 16 radial velocities obtained with HIRES on Keck 1 over a one year span. Although the velocities do not lead to a reliable orbit and mass determination, they are able to constrain the mass to a 3{\sigma} upper limit of 124 MEarth, safely in the regime of planetary masses, thus earning the designation Kepler-22b. The radiative equilibrium temperature is 262K for a planet in Kepler-22b's orbit. Although there is no evidence that Kepler-22b is a rocky planet, it is the first confirmed planet with a measured radius to orbit in the Habitable Zone of any star other than the Sun.Comment: Accepted to Ap

Resonances of low orders in the planetary system of HD37124

The full set of published radial velocity data (52 measurements from Keck + 58 ones from ELODIE + 17 ones from CORALIE) for the star HD37124 is analysed. Two families of dynamically stable high-eccentricity orbital solutions for the planetary system are found. In the first one, the outer planets c and d are trapped in the 2/1 mean-motion resonance. The second family of solutions corresponds to the 5/2 mean-motion resonance between these planets. In both families, the planets are locked in (or close to) an apsidal corotation resonance. In the case of the 2/1 MMR, it is an asymmetric apsidal corotation (with the difference between the longitudes of periastra $\Delta\omega\sim 60^\circ$), whereas in the case of the 5/2 MMR it is a symmetric antialigned one ($\Delta\omega = 180^\circ$). It remains also possible that the two outer planets are not trapped in an orbital resonance. Then their orbital eccentricities should be relatively small (less than, say, 0.15) and the ratio of their orbital periods is unlikely to exceed $2.3-2.5$.Comment: 28 pages, 10 figures, 3 tables; Accepted to Celestial Mechanics and Dynamical Astronom