Dans ce travail, les équations moyennes des mouvements planétaires sont calculées à un ordre élevé grâce à une programmation en calcul formel dédié. Le système résultant comprend plus de 150 000 termes polynomiaux et fournit une très bonne approximation de l'évolution à long terme du système solaire. Le système d'équations est développé à ordre 2 dans les masses et l'ordre 5 en excentricité et inclinaison. le système de degré 3 est intégré analytiquement au premier ordre ce qui fournit une solution de plus de 25 000 termes. Le problème des petits diviseurs séculaires est discuté et une liste de petits diviseurs de grande amplitude dans la solution est donnée, le terme principal étant lié à l'argument $ g_1-g_5 + (s_2-s_1) $ qui intervient dans l'excentricité de Mercure et de Jupiter, et dans l'inclinaison de Mercure et de Vénus. Il est démontré que la présence de ces diviseurs compromet grandement la construction d'une solution analytique, sans tenir compte du fait que la même solution de degré 5 comprend plus de 3 000 000 termes. Les équations séculaires sont ensuite intégrées numériquement d'une façon très efficace sur plus de 1 million d'années, avec un pas de 500 ans pour toutes les planètes, après avoir rajouté la contribution moyenne de la relativité générale et de la Lune. Par comparaison avec les éphémérides DE102 sur plus de 3000 ans, il est démontré que les équations séculaires représentent très bien l'évolution à long terme du système solaire.

Cette thèse traite de la modélisation des marées fluides des planètes telluriques du système solaire et des systèmes exoplanétaires.En premier lieu, nous examinons la réponse de marée des couches atmosphériques, soumises au potentiel de marée gravifique et au forçage thermique de l’étoile hôte du système. Nous proposons un nouveau modèle global prenant en compte les processus dissipatifs avec un refroidissement newtonien, modèle à partir duquel nous traitons la dynamique des ondes de marées engendrées par ces forçages, et quantifions leur dissipation, le nombre de Love et le couple de marée exercé sur la couche atmosphérique en fonction de la fréquence de forçage. Ceci nous permet d'étudier l'ensemble des configurations possibles depuis les planètes au voisinage de la synchronisation telles que Vénus jusqu'aux rotateurs rapides tels que la Terre.En second lieu, nous développons une approche similaire pour les océans de planètes de type terrestre, où la friction visqueuse effective de la topographie est prise en compte, à partir de laquelle nous quantifions la réponse de marée d’un océan global potentiellement profond et sa dépendance à la fréquence d’excitation. Dans ce cadre, et ce grâce à des modèles locaux, nous caractérisons de manière détaillée les propriétés des spectres en fréquence de la dissipation engendrée par les ondes de marées au sein des couches fluides planétaires (et stellaires) en fonction des paramètres structurels et dynamiques de ces dernières (rotation, stratification, viscosité et diffusivité thermique). / This thesis deals with the modeling of fluid tides in terrestrial planets of the Solar system and exoplanetary systems.First, we examine the tidal response of atmospheric layers, submitted to the tidal gravitational potential and the thermal forcing of the host star of the system. We propose a new global model taking into account dissipative processes with a Newtonian cooling, model that we use to treat the dynamics of tidal waves generated by these forcings, and to quantify their dissipation, the Love number and the tidal torque exerted on the atmospheric layer as a function of the forcing frequency. This allows us to study possible configurations from planets close to synchronization such as Venus to rapid rotators such as the Earth.Second, we develop a similar approach for the oceans of terrestrial planets where the action of topography is taken into account thanks to an effective viscous friction. From this modeling, we quantify the tidal response of a potentially deep global ocean and its dependence of the tidal frequency. In this framework, and by using local models, we characterize in detail the properties of the frequency spectra of dissipation generated by tidal waves within fluid planetary (and stellar) layers as functions of the structural and dynamical parameters of these latters (rotation, stratification, viscosity and thermal diffusivity).

Le sujet de cette thèse est l'étude de la rotation à long terme des corps célestes.La première partie est consacrée à l’étude de la rotation à long terme de Cérès et Vesta, les deux corps les plus massifs de la ceinture principale d’astéroïdes. Ils sont l’objet d’étude de la sonde spatiale Dawn, qui a permis de déterminer précisément les caractéristiques physiques et de rotation nécessaires au calcul de leurs rotations. La distribution de glace sous et à la surface de Cérès dépend du mouvement de son axe de rotation par le biais de l’obliquité, inclinaison de l’équateur sur l’orbite. Les rotations de Cérès et Vesta étant rapides, l’évolution à long terme des axes de rotation de Cérès et Vesta a été obtenue à l'aide d'une intégration symplectique des équations de la rotation, où une moyenne a été réalisée sur la rotation propre rapide. La stabilité des axes de rotation de Cérès et Vesta a été étudiée en fonction des paramètres de la rotation avec un modèle séculaire semi-analytique, qui a permis de montrer que les axes de rotation ne présentaient pas de caractère chaotique.La seconde partie concerne le développement d'intégrateurs symplectiques dédiés au corps solide. L'intégration de la rotation propre d'un corps solide nécessite d’intégrer les équations issues du hamiltonien du corps solide libre. Ce hamiltonien est certes intégrable et présente une solution explicite nécessitant l’usage des fonctions elliptiques de Jacobi, cependant le coût numérique de ces fonctions est élevé. Lorsque le hamiltonien du corps solide libre est couplé avec une énergie potentielle, l’orientation du corps doit être calculée à chaque pas d’intégration, ce qui augmente le temps de calcul. Des intégrateurs symplectiques ont ainsi été précédemment proposés pour le corps solide libre. Dans ce travail, des intégrateurs spécifiques au corps solide ont été développés en utilisant les propriétés de l’algèbre de Lie du moment cinétique. / This thesis concerns the long-term rotation of celestial bodies.The first part is a study of the long-term rotation of Ceres and Vesta, the two heaviest bodies of the main asteroid belt. The spacescraft Dawn studied these two objects and determined the physical and rotational characteristics, which are necessary for the computation of their rotations. The ice distribution under and on the surface of Ceres depends on the evolution of the obliquity, which is the inclination of the equatorial plane on the orbital plane. As the rotations of Ceres and Vesta are fast, the long-term evolution of the spin axes of Ceres and Vesta was obtained by realizing a symplectic integration of the equations of the rotation averaged on the fast proper rotation. The stability of the spin axes of Ceres and Vesta was studied with respect to the parameters of the rotation with a secular and semi-analytical model, which allowed to show that the spin axes are not chaotic.The second part concerns the development of symplectic integrators dedicated to the rigid body. The integration of the proper rotation of a rigid body needs to integrate the equations given by the Hamiltonian of the free rigid body. This Hamiltonian is integrable and presents an explicit solution using the Jacobi elliptic functions. However, the numerical cost of these functions is high. When the Hamiltonian of the free rigid body is coupled to a potential energy, the orientation of the body is needed at each step, which increases the computation time. Symplectic integrators were then previously proposed for the free rigid body. In this work, symplectic integrators dedicated to the rigid body were developed using the properties of the Lie algebra of the angular momentum.

Les observations de centaines de systèmes d’exoplanètes nous ont fourni un large échantillon de configurations orbitales. Les périodes orbitales figurent parmi les données les mieux connues et les plus étonnantes. Les Super-Terres, ces planètes caractérisées par une masse entre 1 et 20 masses terrestres et une période typiquement de moins de 100 jours, sont présentes autour de la plupart des étoiles. La distribution des rapports de leurs périodes orbitales défie les astrophysiciens : pendant leur formation et migration au sein de leur disque protoplanétaire, elles devraient former des chaînes de résonances de moyen mouvement, c’est-à-dire que les rapports des périodes orbitales de planètes voisines devraient être proches de fractions simples. Toutefois, la plupart des systèmes de Super-Terres ne sont pas résonants. Dans cette thèse, je traite les aspects clés des chaînes résonantes : leur formation, leur évolution et leur stabilité. Premièrement, j’introduis les idées modernes en théorie de formation planétaire, et les méthodes utilisées dans la thèse : la mécanique Hamiltonienne, le problème planétaire et la théorie perturbative. Deuxièmement, je présente le processus de capture en résonance de moyen mouvement du premier ordre k : k − 1 par migration convergente des planètes, avec une nouvelle description analytique de l’évolution planétaire qui en suit, et je décris la dynamique résonante dans le plan orbital commun. La description analytique est confirmée par des intégrations N-corps qui incluent les interactions disque-planète. Ensuite, je me base sur des résultats existants concernant l’évolution dissipative de deux planètes en résonance qui engendre la divergence de leurs demi-grands axes. Par une approche similaire, je présente une méthode statistique qui permet de déterminer dans quelle mesure l’architecture observée d’un système de trois planètes est compatible avec une histoire dynamique résonante dissipative. Je considère par la suite la stabilité des chaînes résonantes. Des études antérieures ont montré que l’absence de systèmes exoplanétaires résonants n’est pas en contradiction avec le modèle de capture en résonance par migration dans le disque, si une phase d’instabilité est très commune après la disparition du disque. On observe un taux d’instabilité plus élevé dans les systèmes synthétiques plus compacts et peuplés par des planètes plus massives. Des simulations N-corps dédiées à l’étude de la stabilité des chaînes résonantes ont montré qu’il y a une masse planétaire maximale qui garantit la stabilité ; cette masse limite diminue si les planètes sont plus massives et/ou si la chaîne résonante est plus compacte. J’étudie la stabilité des chaînes résonantes de planètes en fonction de leur masse commune, et j’examine de façon analytique et numérique des cas spécifiques de systèmes comprenant deux ou trois planètes. Je découvre un mécanisme dynamique qui peut déclencher une excitation du système, et qui mène à une phase de rencontres proches et collisions. Ce mécanisme se généralise à différents nombres de planètes et/ou à des chaînes résonantes plus ou moins compactes, et donne une prédiction analytique de la masse critique qui est en accord qualitatif avec les expériences numériques mentionnées précédemment. Enfin, je décris un scénario dynamique qui peut expliquer la pollution des naines blanches en éléments lourds. Les systèmes planétaires compacts peuvent devenir instables pendant la phase de perte de masse qui marque la fin de l’évolution stellaire, et les impacts entre planètes génèrent des débris. En m’appuyant sur des résultats précédents, je montre que l’excentricité orbitale des débris qui résident en résonance de moyen mouvement avec une planète externe peut devenir suffisamment élevée pour que les débris soient engloutis par l’étoile, ce qui peut expliquer la pollution observée. / Observations of hundreds of exoplanetary systems have produced a huge sample of orbital configurations, and the orbital periods are one of their better constrained and most astonishing properties. A common type of exoplanets are the Super-Earths, which have a mass between 1 and 20 Earth masses and a typical period of less than 100 days. The period ratio distribution of these planets poses a challenge to astrophysicists: during their formation, still embedded in the protoplanetary disc, we expect them to form chains of mean motion resonances, where the period ratio of neighbouring planets is close to a low-integer ratio. However, most Super-Earth systems are not close to resonance. In this thesis, I discuss key dynamical aspects of resonant chains: their formation, their evolution and their stability. I first give an overview of our current understanding of planetary formation, and an introduction of the methods used in the thesis: the tools of Hamiltonian dynamics, the planetary problem and perturbation theory. Then, I present the process of capture of planets migrating in protoplanetary discs into first order k : k − 1 mean motion resonances, including a novel analytical description of the corresponding planetary evolution, and I describe the relevant aspects of resonant dynamics in the planar approximation. The analytical treatment is supported by numerical N-body simulations which include the planet-disc interactions. Next, I expand on previous results on two-planet dissipative evolution in mean motion resonance and the resulting divergence of the planets’ semi-major axes. With a similar approach, I present a statistical method which allows to determine to what extent the observed architecture of a three-planet system is compatible with a dissipative resonant dynamical history. I then address the main problem of the stability of resonant chains. Previous works have shown that the over-all lack of resonances in the exoplanet sample is not in contradiction with resonant capture, if a post-disc phase of planetary instabilities is extremely common. Higher rates of instabilities are observed in synthetic systems where planets are most massive and the configurations most compact. Specific N-body experiments on the stability of resonant chains found that there is a critical planetary mass allowed for stability, which decreases with increasing number of planets and/or increasing value of k in the chain. The origin of these instabilities was however not discussed. I study the stability of resonant chains of equal-mass planets in terms of their mass, investigating analytically and numerically specific cases of two- and threeplanet systems. I find a dynamical mechanism which can trigger an excitation of the system, leading to mutual close-encounters and collisions. This can be generalised to an arbitrary number of planets and/or value of k in the resonant chain, and gives an analytical prediction for the critical mass allowed for stability which agrees qualitatively with the aforementioned numerical experiments. Finally, I describe a dynamical scenario that can explain the pollution of White Dwarfs with heavy elements. The idea is that compact planetary systems become unstable during the mass-loss phase characterising the end of the stellar evolution, so that impacts among planets lead to the generation of collisional debris. Expanding on previous works, I show that debris residing in mean motion resonance with an outer planetary perturber can have their orbital eccentricity excited to largeenough values to be engulfed by the host star, causing the observed pollution.

Cette thèse contient deux parties distinctes, reliées par le thème de l’étude géométrique du problème à trois corps. La première partie présente un point de vue sur les enjeux et les perspectives liés à la diffusion des mathématiques, et illustre ce point de vue à l’aide de deux projets de diffusion « grand public » : une exposition virtuelle autour de la mécanique céleste et du problème à trois corps, et un duo de contes mathématiques pour enfants, l’un sur la forme de la lune, et l’autre sur l’enlacement de courbes fermées. La présentation de ces projets est suivie d’une analyse a priori et d’une étude des observations recueillies lors de différentes expérimentations auprès de publics variés. La deuxième partie est consacrée à l’étude – théorique et numérique – de l’enlacement des trajectoires de quelques systèmes dynamiques sur la 3-sphère, et en particulier de certaines instances du problème à trois corps. On y présente d’abord le problème à trois corps restreint, plan, circulaire, en s’intéressant tout particulièrement au cas où une des deux primaires disparait. On se ramène ainsi à un flot sur la 3-shpère dont on connaît explicitement des sections de Birkhoff en disque ou en anneau, et on met en lumière des éléments qui tendent à montrer le caractère lévogyre de ce flot. On explore ensuite, à l’aide de simulations numériques, la possibilité que le système reste lévogyre sur un domaine assez éloigné de ce cas dégénéré. Enfin, on s’intéresse aux flots sur la 3-sphère qui admettent une section de Birkhoff en disque et on traduit la notion d’enlacement de mesures invariantes pour le flot en termes d’enroulement de mesures invariantes pour le difféomorphisme de premier retour. / This thesis contains two distinct parts, connected by the subject of the geometric study of the three body problem.The first part presents a point of view about the stakes and prospects of the popularization of mathematics, and it illustrates this point of view with two projects of popularization for a general public : a virtual exhibition about celestial mechanics and the three body problem, and a pair of mathematical tales for children, one about the shape of the moon, and the other about the linking number of two closed curves. The presentation of these projects is followed by an initial analysis and by a study of the observations collected during different experimentations towards various publics. The second part is devoted to the theoretical and computational study of the linking number of trajectories from a few dynamical systems on the 3-sphere, and in particular from some cases of the restricted three body problem. We first present the planar, circular, restricted three body problem, with a particular attention to the case where one of the two heavy bodies vanishes. We thus restrict ourselves to a flow on the 3-shpere for which disk-like or annular-like Birkhoff sections are explicitely known, and we bring to light evidences of the right-handedness of this flow. Then we investigate, with the help of computer simulations, the possibility for the system to stay right-handed over a domain rather distant from this degenerate case. Finally, we consider the flows on the 3-sphere which admit a disk-like Birkhoff section, and we translate the notion of linking for measures that are invariant by a flow into the notion of winding for measures that are invariant by the first return map on the disk.

Cette thèse considère quatre systèmes physiques complexes pour lesquels il est exceptionnellement possible d’identifier des variables lentes qui contrôlent l'évolution à temps long du système complet. La séparation d'échelle de temps entre ces variables lentes et les autres variables permet d'utiliser la technique de moyennisation stochastique pour obtenir une dynamique effective pour les variables lentes. Cette thèse considère la possibilité de prédire les événements rares dans le système solaire. Nous avons étudié deux types d’événements rares. Le premier est un renversement possible de l'axe de rotation de la Terre en l'absence des effets de marée de la Lune. Le second est la désintégration de l'ensemble du système solaire interne suite à une instabilité dans l'orbite de Mercure. Pour chacun des deux problèmes, il existe des variables lentes non triviales, qui ne sont pas données par des variables physiques naturelles. La moyennisation stochastique a permis de découvrir le mécanisme physique qui conduit à ces événements rares et de donner, par une approche purement théorique, l'ordre de grandeur de la probabilité de ces phénomènes. Nous avons également montré que la déstabilisation de Mercure sur un temps inférieur à l'âge du système solaire obéit à un mécanisme d'instanton bien décrit par la théorie des grandes déviations. Le travail effectué dans cette thèse ouvre donc un nouveau champ d'action pour l'utilisation d'algorithmes de calcul d'événements rares. Nous avons utilisé pour la première fois les théorèmes de moyennisation stochastique dans le cadre de la mécanique céleste pour quantifier l'effet stochastique des astéroïdes sur la trajectoire des planètes. Enfin, une partie du travail porte sur un problème de turbulence géophysique: dans l'atmosphère de Jupiter, on peut observer des structures zonales (jets) à grande échelles évoluant beaucoup plus lentement que les tourbillons environnants. Nous montrons qu'il est pour la première fois possible d'obtenir explicitement le profil de ces jets par moyennisation des degrés de liberté turbulents rapides. / The present thesis describes four complex dynamical systems. In each system, the long-term behavior is controlled by a few number of slow variables that can be clearly identified. We show that in the limit of a large timescale separation between the slow variables and the other variables, stochastic averaging can be performed and leads to an effective dynamics for the set of slow variables. This thesis also deals with rare events predictions in the solar system. We consider two possible rare events. The first one is a very large variation of the spin axis orientation of a Moonless Earth. The second one is the disintegration of the inner solar system because of an instability in Mercury’s orbit. Both systems are controlled by non-trivial slow variables that are not given by simple physical quantities. Stochastic averaging has led to the discovery of the mechanism leading to those rare events and gives theoretical bases to compute the rare events probabilities. We also show that Mercury’s short-term destabilizations (compared to the age of the solar system) follow an instanton mechanism, and can be predicted using large deviation theory. The special algorithms devoted to the computation of rare event probabilities can thus find surprising applications in the field of celestial mechanics. We have used for the first time stochastic averaging in the field of celestial mechanics to give a relevant orders of magnitude for the long-term perturbation of planetary orbits by asteroids. A part of the work is about geophysical fluid mechanics. In Jupiter atmosphere, large scale structures (jets) can be observed, the typical time of evolution of which is much larger than that of the surrounding turbulence. We show for the first time that the mean wind velocity can be obtained explicitly by averaging the fast turbulent degrees of freedom.

Au moins 20% des étoiles de la séquence principale abritent des disques de débris, analogues à la ceinture de Kuiper. Ces disques sont la preuve que l'accumulation de solides a au moins permis la formation de corps de taille kilométrique. Il n'est donc pas surprenant que plusieurs de ces disques soient accompagnés de planètes, qui, en laissant leur empreinte dynamique sur la structure spatiale de ces disques, révèlent leur présence. Par conséquent, la détection d'un disque de débris excentrique entourant Zeta² Ret par le télescope spatial Herschel indique la présence d'un perturbateur massif dans ce système. Zeta² Ret étant un système mature, âgé de 2-3 Gyr, et en ce sens, analogue à notre propre système solaire, il offre un exemple différent d'évolution dynamique à long terme. Cette thèse comprend une modélisation détaillée de la structure du disque de débris de Zeta² Ret, ce qui conduit à des contraintes sur la masse et l'orbite du perturbateur suspecté. Cette étude révèle également que les structures excentriques dans les disques de débris peuvent survivre sur des échelles de temps Gyr.La modélisation de disques de débris peut permettre la découverte postérieure de planètes comme c'est le cas pour le système de Fomalhaut. La forme excentrique de son disque de débris fut d'abord attribuée à Fom b, un compagnon détecté près du bord interne du disque, mais qui se révèle finalement trop excentrique (e~0.6-0.9) pour lui donner sa forme, indiquant la présence d'un autre corps massif, Fom c. Le système planétaire qui en résulte est très instable, ce qui implique une diffusion récente de Fom b sur son orbite actuelle, éventuellement par Fom c. L'étude de ce scénario révèle qu'en ayant résidé dans une résonance de moyen-mouvement interne avec une Fom c excentrique et de masse comparable à Neptune ou Saturne, Fom b aurait subi une augmentation progressive de son excentricité sur des périodes comparables à l'âge du système (~440 Myr), ce qui l'aurait amenée assez proche de Fom c pour subir une diffusion récente, qui, complétée par une évolution séculaire avec Fom c, explique sa configuration orbitale actuelle. Ce mécanisme en trois étapes pourrait également avoir généré d'importantes quantités de matériel en orbites très excentriques, ce qui en retour pourrait alimenter en poussière les parties internes du système. Par conséquent, ce mécanisme pourrait aussi expliquer la présence de ceintures de poussières internes dans le système de Fomalhaut, mais aussi la découverte d'une importante population de ceintures de poussières chaudes et massives dans les systèmes âgés de plus de 100 Myr.Les systèmes planétaires découverts jusqu'ici présentent une grande variété d'architectures, et notre système solaire est loin d'être un modèle générique. Un des principaux mécanismes déterminant la morphologie d'un système planétaire est la migration planétaire. On attend d'un compagnon stellaire - ce que notre système solaire ne possède pas - qu'il affecte les conditions de migration planétaire, et conduise potentiellement à la formation de systèmes planétaires très différents. Ce phénomène est évidemment non négligeable puisque les systèmes binaires représentent au moins la moitié des systèmes stellaires. Dans les systèmes planétaires matures, la migration planétaire peut se produire suite à l'interaction avec le matériel solide et l'impact de la binarité sur cette migration tardive est exploré dans cette thèse. Un compagnon circumstellaire pourrait inverser la tendance à la migration interne des planètes dans les systèmes stellaires simples, et rapprocher ces planètes des régions perturbées par le compagnon binaire, où elles ne pourraient pas s'être formées in situ. Cela pourrait fournir une explication à la présence de planètes qui présentent des signes de migration externe vers un compagnon circumstellaire dans les systèmes de Gamma Cephei et HD 196885. / At least 20% of Main-Sequence stars are known to harbor debris disks analogs to the Kuiper Belt. These disks are proof that the accretion of solids has permitted the formation of at least km-sized bodies. It is thus not surprising that several of these disks are accompanied by planets, which may reveal themselves by setting their dynamical imprints on the spatial structure of debris disks. Therefore, the detection of an eccentric debris disk surrounding Zeta² Ret by the Herschel space telescope provides evidence for the presence of a massive perturber in this system.Zeta² Ret being a mature Gyr-old system, and in that sense, analogous to our own Solar System, it offers a different example of long-term dynamical evolution. This thesis includes a detailed modeling of the structure of the debris disk of Zeta² Ret, which leads to constraints on the mass and orbital characteristics of the putative perturber. This study also reveals that eccentric structures in debris disks can survive on Gyr timescales.Detailed modeling of the structure of debris disks can allow the posterior discovery of hidden planets as is the case for the Fomalhaut system. The eccentric shape of the debris disk observed around this star was first attributed to Fom b, a companion detected near the belt inner-edge, which revealed to be highly eccentric (e~0.6-0.9), and thus very unlikely shaping the belt. This hints at the presence of another massive body in this system, Fom c, which drives the debris disk shape. The resulting planetary system is highly unstable, which involves a recent scattering of Fom b on its current orbit, potentially with the yet undetected Fom c. This scenario is investigated in this thesis and its study reveals that by having resided in inner mean-motion resonance with a Neptune or Saturn-mass belt-shaping eccentric Fom c and therefore have suffered a gradual resonant eccentricity increase on timescales comparable to the age of the system (~440 Myr), Fom b could have been brought close enough to Fom c and suffered a recent scattering event, which, complemented by a secular evolution with Fom c, explains its current orbital configuration. This three-step scenario also implies that significant amounts of material may have been set on extremely eccentric orbits such as this of Fom b through this mechanism, which in return could feed in dust the inner parts of the system. Therefore, this mechanism may also explain the presence of inner dust belts in the Fomalhaut system, but also the discovery a significant population of very bright hot dust belts in systems older than 100 Myr.The planetary systems discovered so far exhibit a great variety of architectures, and our solar system is far from being a generic model. One of the main mechanism that determines a planetary system morphology is planetary migration. The presence of a stellar binary companion - which our solar system is deprived of - is expected to affect planetary migration conditions, and potentially lead to the formation of very different planetary systems. This phenomenon is obviously non-negligible since binary systems represent at least half of stellar systems. At late stages of planetary systems formation, planetary migration may occur as the result of interactions with remaining solid planetesimals and the impact of binarity on this planetesimal-driven migration is explored in this thesis. A stellar binary companion may in fact reverse the tendency for planets in single star systems to migrate inwards, and bring them closer to regions perturbed by the binary companion, where they could not have formed in situ. This may give an explanation for the presence of planets which present signs of outward migration towards a circumstellar companion in the Gamma Cephei and HD 196885 systems.