Recent observations from NASA's Kepler mission detected the first planets in circumbinary orbits. The question we try to answer is where these planets formed in the circumbinary disk and how far inside they migrated to reach their present location. We investigate the first and more delicate phase of planet formation when planetesimals accumulate to form planetary embryos. We use the hydrodynamical code FARGO to study the evolution of the disk and of a test population of planetesimals embedded in it. With this hybrid hydrodynamical--N--body code we can properly account for the gas drag force on the planetesimals and for the gravitational force of the disk on them. The numerical simulations show that the gravity of the eccentric disk on the planetesimal swarm excites their eccentricities to values much larger than those induced by the binary perturbations only within 10 AU from the stars. Moreover, the disk gravity prevents a full alignment of the planetesimal pericenters. Both these effects lead to large impact velocities, beyond the critical value for erosion. Planetesimals accumulation in circumbinary disks appears to be prevented close to the stellar pair by the gravitational perturbations of the circumbinary disk. The observed planets possibly formed in the outer regions of the disk and then migrated inside by tidal interaction with the disk.

Discs in binaries have a complex behavior because of the perturbations of the companion star. Planet formation in binary-star systems both depend on the companion star parameters and on the properties of the circumstellar disc. An eccentric disc may increase the impact velocity of planetesimals and therefore jeopardize the accumulation process. We model the evolution of discs in close binaries including the effects of self-gravity and adopting different prescriptions to model the disc's radiative properties. We focus on the dynamical properties and evolutionary tracks of the discs. We use the hydrodynamical code FARGO and we include in the energy equation heating and cooling effects. Radiative discs have a lower disc eccentricity compared to locally isothermal discs with same temperature profile. As a consequence, we do not observe the formation of an internal elliptical low density region as in locally isothermal disc models. However, the disc eccentricity depends on the disc mass through the opacities. Akin to locally isothermal disc models, self-gravity forces the disc's longitude of pericenter to librate about a fixed orientation with respect to the binary apsidal line ($\pi$). The disc's radiative properties play an important role in the evolution of discs in binaries. A radiative disc has an overall shape and internal structure that are significantly different compared to a locally isothermal disc with same temperature profile. This is an important finding both for describing the evolutionary track of the disc during its progressive mass loss, and for planet formation since the internal structure of the disc is relevant for planetesimals growth in binary systems. The non-symmetrical distribution of mass in these discs causes large eccentricities for planetesimals that may affect their growth.

Hydrodynamical simulations of two giant planets embedded in a gaseous disk have shown that in case of a smooth convergent migration they end up trapped into a mean motion resonance. These findings have led to the conviction that the onset of dynamical instability causing close encounters between the planets can occur only after the dissipation of the gas when the eccentricity damping is over. We show that a system of three giant planets may undergo planet-planet scattering when the gaseous disk, with density values comparable to that of the Minimum Mass Solar Nebula, is still interacting with the planets. The hydrodynamical code FARGO--2D--1D is used to model the evolution ofthe disk and planets, modified to properly handle close encounters between the massive bodies. Our simulations predict a variety of different outcomes of the scattering phase which includes orbital exchange, planet merging and scattering of a planet in a hyperbolic orbit. This implies thatthe final fate of a multiplanet system under the action of the disk torques is not necessarily a packed resonant configuration.

We study the evolution of circumstellar massive disks around the primary star of a binary system focusing on the computation of disk eccentricity. In particular, we concentrate on its dependence on the binary eccentricity. Self-gravity is included in our numerical simulations. Our standard model assumes a semimajor axis for the binary of 30 AU, the most probable value according to the present binary statistics.

We explore planet formation in binary systems around the central star where the protoplanetary disk plane is highly inclined with respect to the companion star orbit. This might be the most frequent scenario for binary separations larger than 40 AU, according to Hale (1994). We focus on planetesimal accretion and compute average impact velocities in the habitable region and up to 6 AU from the primary.

In a previous paper (Gayon & Bois 2008a), we have shown the general efficiency of retrograde resonances for stabilizing compact planetary systems. Such retrograde resonances can be found when two-planets of a three-body planetary system are both in mean motion resonance and revolve in opposite directions. For a particular two-planet system, we have also obtained a new orbital fit involving such a counter-revolving configuration and consistent with the observational data. In the present paper, we analytically investigate the three-body problem in this particular case of retrograde resonances. We therefore define a new set of canonical variables allowing to express correctly the resonance angles and obtain the Hamiltonian of a system harboring planets revolving in opposite directions. The acquiring of an analytical "rail" may notably contribute to a deeper understanding of our numerical investigation and provides the major structures related to the stability properties. A comparison between our analytical and numerical results is also carried out.

Recent studies have shown that alpha Centauri B might be, from an observational point of view, an ideal candidate for the detection of an Earth-like planet in or near its habitable zone (0.5-0.9AU). We study here if such habitable planets can form, by numerically investigating the planet-formation stage which is probably the most sensitive to binarity effects: the mutual accretion of km-sized planetesimals. Using a state-of-the-art algorithm for computing the impact velocities within a test planetesimal population, we find that planetesimal growth is only possible, although marginally, in the innermost part of the HZ around 0.5AU. Beyond this point, the combination of secular perturbations by the binary companion and gas drag drive the mutual velocities beyond the erosion limit. Impact velocities might later decrease during the gas removal phase, but this probably happens too late for preventing most km-sized objects to be removed by inward drift, thus preventing accretion from starting anew. A more promising hypothesis is that the binary formed in a crowded cluster, where it might have been wider in its initial stages, when planetary formation was ongoing. We explore this scenario and find that a starting separation roughly 15 AU wider, or an eccentricity 2.5 times lower than the present ones are required to have an accretion-friendly environment in the whole HZ.

The five planets discovered around the main-sequence star 55 Cnc may represent a case of stable chaos. By using both the Frequency Map Analysis and MEGNO we find that about 15 % of the systems that can be build from the nominal orbital elements of the system are highly chaotic. However, in spite of the fast diffusion rate in the phase space, the planetary system is not destabilized over 400 Myr and close encounters between the planets are avoided.

We numerically explore planet formation around alpha Cen A by focusing on the crucial planetesimals-to-embryos phase. Our code computes the relative velocity distribution, and thus the accretion vs. fragmentation trend, of planetesimal populations having any given size distribution. This is a critical aspect of planet formation in binaries since the pericenter alignment of planetesimal orbits due to the gravitational perturbations of the companion star and to gas friction strongly depends on size. We find that, for the nominal case of a MMSN gas disc, the region beyond 0.5AU from the primary is hostile to planetesimal accretion. In this area, impact velocities between different-size bodies are increased, by the differential orbital phasing, to values too high to allow mutual accretion. For any realistic size distribution for the planetesimal population, this accretion-inhibiting effect is the dominant collision outcome and the accretion process is halted. Results are robust with respect to the profile and density of the gas disc: except for an unrealistic almost gas-free case, the inner accretion safe area never extends beyond 0.75AU. We conclude that planet formation is very difficult in the terrestrial region around alpha Cen A, unless it started from fast-formed very large (>30km) planetesimals. Notwithstanding these unlikely initial conditions, the only possible explanation for the presence of planets around 1 AU from the star would be the hypothetical outward migration of planets formed closer to the star or a different orbital configuration in the binary's early history. Our conclusions differ from those of several studies focusing on the later embryos-to-planets stage, confirming that the planetesimals-to-embryos phase is more affected by binary perturbations.

Numerical simulations of planetesimal accretion in circumprimary and circumbinary orbits are described. The secular perturbations by the com- panion star and gas drag are included in our models. We derive limits on the parameters of the binary system for which accretion and then planetary forma- tion are possible. In the circumbinary case we also outline the radial distance from the baricenter of the stars beyond which accumulation always occurs. Hy- drodynamical simulations are also presented to validate our N-body approach based on the axisymmetric approximation for the gas of the disk.

We study the dynamics of Jupiter Trojans in the early phase of the Solar system while the outer planets migrated due to their interaction with the planetesimal disk.

We numerically investigate the possibility of planetesimal accretion in circumbinary disks, under the coupled influence of both stars' secular perturbations and friction due to the gaseous component of the protoplanetary disk. We focus on one crucial parameter: the distribution of encounter velocities between planetesimals in the 0.5 to 100km size range. An extended range of binary systems with differing orbital parameters is explored. The resulting encounter velocities are compared to the threshold velocities below which the net outcome of a collision is accumulation into a larger body instead of mass erosion. For each binary configuration, we derive the critical radial distance from the binary barycenter beyond which planetesimal accretion is possible. This critical radial distance is smallest for equal-mass binaries on almost circular orbits. It shifts to larger values for increasing eccentricities and decreasing mass ratio. The importance of the planetesimals' orbital alignments of planetesimals due to gas drag effects is discussed.

We study to which extent planetesimal accretion is affected by the perturbing presence of a compagnon star. We concentrate on one crucial parameter: the distribution of encounter velocities within the planetesimal swarm. We numerically explore the evolution of this parameter taking into account the secular perturbations of the binary and friction due to the very likely presence of gas in the disk. $<\Delta v>$ maps are derived, for planetesimals of different sizes, for a total of 120 binary configurations (eccentricity eb and separation ab). We identify for each case 3 different accreting behaviours. 1) In regions where no significant dV increase is observed, ``classical'' (i.e., single-star) accretion is possible. 2) In regions where $dV>v_{ero}$, the threshold velocity above which all impacts are eroding, no accretion is possible and planet growth is stopped. 3) In between these 2 limiting behaviours, a large fraction of binary configurations leads to significant dV increase, but still below the erosion threshold. In this intermediate case, planetesimal growth can occur, but proceeds slower than in the single-star case, possibly following the so-called type II runaway groth mode.

We investigate classical planetesimal accretion in a binary star system of separation ab<50AU by numerical simulations, with particular focus on the region at a distance of 1 AU from the primary. The planetesimals orbit the primary, are perturbed by the companion and are in addition subjected to a gas drag force. We concentrate on the problem of relative velocities dv among planetesimals of different sizes. For various stellar mass ratios and binary orbital parameters we determine regions where dv exceed planetesimal escape velocities v_esc (thus preventing runaway accretion) or even the threshold velocity v_ero for which erosion dominates accretion. Gaseous friction has two crucial effects on the velocity distribution: it damps secular perturbations by forcing periastron alignment of orbits, but at the same time the size--dependence of this orbital alignment induces a significant dv increase between bodies of different sizes. This differential phasing effect proves very efficient and almost always increases dv to values preventing runaway accretion, except in a narrow domain of almost circular companion orbits. The erosion threshold dv>v_ero is reached in a wide (ab,eb) space for small (<10km) planetesimals, but in a much more limited region for bigger ~50km objects. In the intermediate v_esc<dv < v_ero domain, a possible growth mode would be the type II runaway growth identified by Kortenkmap et al.(2001)

We present the data acquisition strategy and characterization procedures for the Canada-France Ecliptic Plane Survey (CFEPS), a sub-component of the Canada-France-Hawaii Telescope Legacy Survey. The survey began in early 2003 and as of summer 2005 has covered 430 square degrees of sky within a few degrees of the ecliptic. Moving objects beyond the orbit of Uranus are detected to a magnitude limit of $m_R$=23 -- 24 (depending on the image quality). To track as large a sample as possible and avoid introducing followup bias, we have developed a multi-epoch observing strategy that is spread over several years. We present the evolution of the uncertainties in ephemeris position and orbital elements as the objects progress through the epochs. We then present a small 10-object sample that was tracked in this manner as part of a preliminary survey starting a year before the main CFEPS project. We describe the CFEPS survey simulator, to be released in 2006, which allows theoretical models of the Kuiper Belt to be compared with the survey discoveries since CFEPS has a well-documented pointing history with characterized detection efficiencies as a function of magnitude and rate of motion on the sky. Using the pre-survey objects we illustrate the usage of the simulator in modeling the classical Kuiper Belt.

We numerically investigate under which conditions the planet detected at 2.1 AU of Gamma-Cephei could form through the core-accretion scenario despite the perturbing presence of the highly eccentric companion star. We first show that the initial stage of runaway accretion of kilometer-sized planetesimals is possible within 2.5 AU from the central star only if large amounts of gas are present. In this case, gaseous friction induces periastron alignment of the orbits which reduces the otherwise high mutual impact velocities due to the companion's secular perturbations. The following stage of mutual accretion of large embryos is also modeled. According to our simulations, the giant impacts among the embryos always lead to a core of 10 Mearth within 10 Myr, the average lifetime of gaseous discs. However, the core always ends up within 1.5 AU from the central star. Either the core grows more quickly in the inner region of the disc, or it migrates inside by scattering the residual embryos

At present the possible existence of planets around the stars of a close binary system is still matter of debate. Can planetary bodies form in spite of the strong gravitational perturbations of the companion star? We study in this paper via numerical simulation the last stage of planetary formation, from embryos to terrestrial planets in the Alpha Cen system, the prototype of close binary systems. We find that Earth class planets can grow around Alpha Cen A on a time-scale of 50 Myr. In some of our numerical models the planets form directly in the habitable zone of the star in low eccentric orbits. In one simulation two of the final planets are in a 2:1 mean motion resonance that, however, becomes unstable after 200 Myr. During the formation process some planetary embryos fall into the stars possibly altering their metallicity.

We have performed a series of high resolution N-body experiments on a Connection Machine CM-5 in order to study the stability of collisionless self-gravitating spherical systems. We interpret our results in the framework of symplectic mechanics, which provides the definition of a new class of particular perturbations: The preserving perturbations, which are a generalization of the radial ones. Using models defined by the Ossipkov-Merritt algorithm, we show that the stability of a spherical anisotropic system is directly related to the preserving or non-preserving nature of the perturbations acting on the system. We then generalize our results to all spherical systems. Since the ``isotropic component'' of the linear variation of the distribution function cannot be used to predict the stability or instability of a spherical system, we propose a more useful stability parameter which is derived from the ``anisotropic'' component of the linear variation.

We propose a morphological multi-scale analysis of large-scale structures obtained by computer simulations and by observations. Structures are obtained at different scales by applying a wavelet transform on the observed and simulated data. Application of a segmentation algorithm allows a quantitative morphological description of the structures at each scale. The morphological parameter which we propose represents the deviation of a structure from sphericity. The dependence on scale of this simple parameter is shown to characterize in simulations the underlying cosmological model. We compare the HDM, CDM and mixed models with the CfA catalogue. This comparison favours a mixed model containing $65\%$ of CDM, $30\%$ of HDM, and $5\%$ baryonic matter.

We have performed a series of N-body experiments on Connection Machine-5 in order to simulate the formation of galaxy clusters gravitationally dominated by a massive dark background. In accordance with previous authors we find an extremely inhomogeneous evolution where subcondensations are continually formed and merged. The final distribution of galaxies is more centrally condensed than that of dark matter particles. We have analyzed the origin of this galaxy-dark matter segregation and also the origin of subclustering leading to this effect. We have then analysed, analytically and numerically, the dependence of final segregation on the physical parameters characterizing the model of protocluster. We also find that such a segregation persists even when inelastic encounters of dark halos around galaxies are taken into account. We conclude that this effect cannot be in general avoided in any hierarchical clustering scenario.