Binary and multiple star systems are extreme environments for the formation and long-term presence of extrasolar planets. Circumstellar planets are subject to gravitational perturbations from the distant companion star, and this interaction leads to a long-period precession of their orbits. We investigate analytical models that allow to quantify these perturbations and calculate the secular precession frequency in the dynamical model of the restricted three-body problem. These models are applied to test cases and we discuss some of their shortcomings. In addition, we introduce a modified Laplace-Lagrange model which allows to obtain better frequency estimates than the traditional model for large eccentricities of the perturber. We then generalize this model to any number of perturbers, and present an application to the four-body problem.

We investigate the collective behaviour of particle orbits in the vicinity of magnetic reconnection in Earth's magneto-tail. Various regions of different kinds of orbital stability of particle motions are found. We locate regimes of temporary capture of particle orbits in configuration space as well as locations, where strong particle accelerations take place. With this study we are able to provide a detailed map, i.e. the topology, of high and low acceleration centers close to the reconnection site. Quasi-regular and chaotic kinds of motions of elementary particles can be determined as well. The orbital stability of particle orbits is obtained by a statistical analysis of the outcome of the system of variational equations of particle orbits within the framework of particle-in-cell simulations. Using the concept of Lyapunov Characteristic Numbers to ensembles of particle orbits we introduce Lyapunov Ensemble Averages to describe the response of particle orbits to local perturbations induced by the electro-magnetic field.

In this work, we investigate the dynamical effects of a sequence of close encounters over 200 Myr varying in the interval of 10000 -- 100000 au between a binary star system and passing stars with masses ranging from 0.1$M_{\odot}$ to 10$M_{\odot}$. We focus on binaries consisting of two Sun-like stars with various orbital separations $a_{\scriptscriptstyle 0}$ from 50 au to 200 au initially on circular-planar orbits. We treat the problem statistically since each sequence is cloned 1000 times. Our study shows that orbits of binaries initially at $a_{\scriptscriptstyle 0}$ = 50 au will slightly be perturbed by each close encounter and exhibit a small deviation in eccentricity (+0.03) and in periapsis distance (+1 and -2 au) around the mean value. However increasing $a_{\scriptscriptstyle 0}$ will drastically increase these variances: up to +0.45 in eccentricity and between +63 au and -106 au in periapsis, leading to a higher rate of disrupted binaries up to 50% after the sequence of close encounters. Even though the secondary star can remain bound to the primary, $\sim$20% of the final orbits will have inclinations greater than 10$^{\circ}$. As planetary formation already takes place when stars are still members of their birth cluster, we show that the variances in eccentricity and periapsis distance of Jupiter- and Saturn-like planets will inversely decrease with $a_{\scriptscriptstyle 0}$ after successive fly-bys. This leads to higher ejection rate at $a_{\scriptscriptstyle 0}$ = 50 au but to a higher extent for Saturn-likes (60%) as those planets' apoapsis distances cross the critical stability distance for such binary separation.

Context. In tight binary star systems, tidal interactions can significantly influence the rotational and orbital evolution of both stars, and therefore their activity evolution. This can have strong effects on the atmospheric evolution of planets that are orbiting the two stars. Aims. In this paper, we aim to study the evolution of stellar rotation and of X-ray and ultraviolet (XUV) radiation in tight binary systems consisting of two solar mass stars and use our results to study planetary atmosphere evolution in the habitable zones of these systems. Methods. We have applied a rotation model developed for single stars to binary systems, taking into account the effects of tidal interactions on the rotational and orbital evolution of both stars. We used empirical rotation-activity relations to predict XUV evolution tracks for the stars, which we used to model hydrodynamic escape of hydrogen dominated atmospheres. Results. When significant, tidal interactions increase the total amount of XUV energy emitted, and in the most extreme cases by up to factor of $\sim$50. We find that in the systems that we study, habitable zone planets with masses of 1~M$_\oplus$ can lose huge hydrogen atmospheres due to the extended high levels of XUV emission, and the time that is needed to lose these atmospheres depends on the binary orbital separation.For some orbital separations, and when the stars are born as rapid rotators, it is also possible for tidal interactions to protect atmospheres from erosion by quickly spinning down the stars. For very small orbital separations, the loss of orbital angular momentum by stellar winds causes the two stars to merge. We suggest that the merging of the two stars could cause previously frozen planets to become habitable due to the habitable zone boundaries moving outwards.

The scenario and efficiency of water transport by icy asteroids and comets are still amongst the most important unresolved questions of planetary systems. A better understanding of cometary dynamics in extrasolar systems shall provide information about cometary reservoirs and give an insight into water transport especially to planets in the habitable zone. The detection of Proxima Centauri-b (PCb), which moves in the habitable zone of this system, triggered a debate whether or not this planet can be habitable. In this work, we focus on the stability of an additional planet in the system and on water transport by minor bodies. We perform numerous N-body simulations with PCb and an outer Oort-cloud like reservoir of comets. We investigate close encounters and collisions with the planet, which are important for the transport of water. Observers found hints for a second planet with a period longer than 60 days. Our dynamical studies show that two planets in this system are stable even for a more massive second planet (~12 Earth masses). Furthermore, we perform simulations including exocomets, a second planet, and the influence of the binary Alpha Centauri. The studies on the dynamics of exocomets reveal that the outer limit for water transport is around 200 au. In addition we show that water transport would be possible from a close-in planetesimal cloud (1-4 au). From our simulations, based on typical M-star protoplanetary disks, we estimate the water mass delivered to the planets up to 51 Earth oceans.

We investigate the role of secular and mean motion resonances on the water transport from a belt of icy asteroids onto planets or embryos orbiting inside the circumprimary habitable zone (HZ) of a binary star system. In addition, the host-star has an accompanying gas giant planet. For a comparison, we perform two case studies where a secular resonance (SR) is located either inside the HZ close to 1.0 au (causing eccentric motion of a planet or embryos therein) or in the asteroid belt, beyond the snow line. In the latter case, a higher flux of icy objects moving towards the HZ is expected. Collisions between asteroids and objects in the HZ are treated analytically. Our purely dynamical study shows that the SR in the HZ boosts the water transport; however, collisions can occur at very high impact speeds. In this paper, we treat for the first time, realistic collisions using a GPU 3D-SPH code to assess the water loss in the projectile. Including the water loss into the dynamical results, we get more realistic values for the water mass fraction of the asteroid during an impact. We highlight that collisions occurring at high velocities greatly reduce the water content of the projectile and thus the amount of water transported to planets or embryos orbiting inside the HZ. Moreover, we discuss other effects that could modify our results, namely the asteroid's surface rate recession due to ice sublimation and the atmospheric drag contribution on the asteroids' mass loss.

So far, more than 130 extrasolar planets have been found in multiple stellar systems. Dynamical simulations show that the outcome of the planetary formation process can lead to different planetary architectures (i.e. location, size, mass, and water content) when the star system is single or double. In the late phase of planetary formation, when embryo-sized objects dominate the inner region of the system, asteroids are also present and can provide additional material for objects inside the habitable zone (HZ). In this study, we make a comparison of several binary star systems and aim to show how efficient they are at moving icy asteroids from beyond the snow line into orbits crossing the HZ. We also analyze the influence of secular and mean motion resonances on the water transport towards the HZ. Our study shows that small bodies also participate in bearing a non-negligible amount of water to the HZ. The proximity of a companion moving on an eccentric orbit increases the flux of asteroids to the HZ, which could result in a more efficient water transport on a short timescale, causing a heavy bombardment. In contrast to asteroids moving under the gravitational perturbations of one G-type star and a gas giant, we show that the presence of a companion star not only favors a faster depletion of our disk of planetesimals, but can also bring 4-5 times more water into the whole HZ. However, due to the secular resonance located either inside the HZ or inside the asteroid belt, impacts between icy planetesimals from the disk and big objects in the HZ can occur at high impact speed. Therefore, real collision modeling using a GPU 3D-SPH code show that in reality, the water content of the projectile is greatly reduced and therefore, also the water transported to planets or embryos initially inside the HZ.

We present a survey on exoplanetary systems of binary stars with stellar separations less than 100 au. For a sample of 11 binaries that harbour detected circumstellar giant planets we investigate the frequency of systems with secular resonances (SR) affecting the habitable zone (HZ). Secular resonances are connected to dynamically unstable or chaotic regions by enforcing highly eccentric motion. We apply a semi-analytical method to determine the locations of linear SR, which is based on finding the apsidal precession frequencies of the massive bodies. For configurations where the giant planet is located exterior to the HZ we find that there is always a SR interior to its orbit, the exact location of the SR strongly depends on the system's architecture. In systems with the giant planet interior to the HZ no SR can occur in the Newtonian framework. Taking into account the general relativistic precession of the perihelion, which increases the precession frequencies, planets with $a < 0.1$ au can cause SR in the HZ. We find two cases where the SR is located inside the HZ, and some more where it is close to the HZ. Generally, giant planets interior to the HZ are more favourable than exterior planets to avoid SR in the HZ. Around the location of the SR weaker mean-motion resonances are excited, and resonance overlap is possible. Existing analytical models are not as accurate as the semi-analytical method in locating the SR and deviate by $\sim 0.1$ au or more.

Gravitational perturbations in multi-planet systems caused by an accompanying star are the subject of this investigation. Our dynamical model is based on the binary star HD41004 AB where a giant planet orbits HD41004 A. We modify the orbital parameters of this system and analyze the motion of a hypothetical test-planet surrounding HD41004 A on an interior orbit to the detected giant planet. Our numerical computations indicate perturbations due to mean motion and secular resonances. The locations of these resonances are usually connected to high eccentricity and highly inclined motion depending strongly on the binary-planet architecture.As the positions of mean motion resonances can easily be determined, the main purpose of this study is to present a new semi-analytical method to determine the location of a secular resonance without huge computational effort.

We suggest an additional possible plasma source to explain part of the phenomena observed for the transiting hot Jupiters WASP-12b and HD 189733b in their ultraviolet (UV) light curves. In the proposed scenario, material outgasses from the molten surface of Trojan satellites on tadpole orbits near the Lagrange points L$_4$ and L$_5$. We show that the temperature at the orbital location of WASP-12b is high enough to melt the surface of rocky bodies and to form shallow lava oceans on them. In case of WASP-12b, this leads to the release of elements such as Mg and Ca, which are expected to surround the system. The predicted Mg and Ca outgassing rates from two Io-sized WASP-12b Trojans are $\approx 2.2 \times 10^{27}$ s$^{-1}$ and $\approx 2.2 \times 10^{26}$ s$^{-1}$, respectively. Trojan outgassing can lead to the apparent lack of emission in Mg{\sc ii}\,h\&k and Ca{\sc ii}\,H\&K line cores of WASP-12. For HD 189733b, the mechanism is only marginally possible due to the lower temperature. This may be one of the reasons that couldn't explain the early ingress of HD 189733b observed in the far-UV (FUV) C{\sc ii} doublet due to absence of carbon within elements outgassed by molten lava. We investigate the long-term stability region of WASP-12b and HD 189733b in case of planar and inclined motion of these satellites and show that unlike the classical exomoons orbiting the planet, Io-sized Trojans can be stable for the whole systems life time.

Secular and mean motion resonances (MMR) are effective perturbations for shaping planetary systems. In binary star systems, they play a key role during the early and late phases of planetary formation, as well as for the dynamical stability of a planetary system. In this study, we aim to correlate the presence of orbital resonances with the rate of icy asteroids crossing the habitable zone (HZ), from a circumprimary disk of planetesimals in various binary star systems. We modelled a belt of small bodies in the inner and outer regions, interior and exterior to the orbit of a gas giant planet, respectively. The planetesimals are equally placed around a primary G-type star and move under the gravitational influence of the two stars and the gas giant. We numerically integrated the system for 50 Myr, considering various parameters for the secondary star. Its stellar type varies from a M- to F-type; its semimajor axis is either 50 au or 100 au, and its eccentricity is either 0.1 or 0.3. For comparison, we also varied the gas giant's orbital and physical parameters. Our simulations highlight that a disk of planetesimals will suffer from perturbations owing to a perturbed gas giant, mean motion, and secular resonances. We show that a secular resonance -- with location and width varying according to the secondary star's characteristics -- can exist in the icy asteroid belt region and overlap with MMRs which have an impact on the dynamical lifetime of the disk. In addition, we point out that in any case, the 2:1 MMR, the 5:3 MMR and the secular resonance are powerful perturbations for the flux of icy asteroids towards the HZ and the transport of water therein.

So far, multiple stellar systems harbor more than 130 extra solar planets. Dynamical simulations show that the outcome of planetary formation process can lead to various planetary architecture (i.e. location, size, mass and water content) when the star system is single or double. In the late phase of planetary formation, when embryo-sized objects dominate the inner region of the system, asteroids are also present and can provide additional material for objects inside the habitable zone (hereafter HZ). In this study, we make a comparison of several binary star systems and their efficiency to move icy asteroids from beyond the snow-line into orbits crossing the HZ. We modeled a belt of 10000 asteroids (remnants from the late phase of planetary formation process) beyond the snow-line. The planetesimals are placed randomly around the primary star and move under the gravitational influence of the two stars and a gas giant. As the planetesimals do not interact with each other, we divided the belt into 100 subrings which were separately integrated. In this statistical study, several double star configurations with a G-type star as primary are investigated. Our results show that small bodies also participate in bearing a non-negligible amount of water to the HZ. The proximity of a companion moving on an eccentric orbit increases the flux of asteroids to the HZ, which could result into a more efficient water transport on a short timescale, causing a heavy bombardment. In contrast to asteroids moving under the gravitational perturbations of one G-type star and a gas giant, we show that the presence of a companion star can not only favor a faster depletion of our disk of planetesimals but can also bring 4 -- 5 times more water into the whole HZ.

As part of the national scientific network 'Pathways to Habitable Worlds' the delivery of water onto terrestrial planets is a key question since water is essential for the development of life as we know it. After summarizing the state of the art we show some first results of the transport of water in the early Solar System for scattered main belt objects. Hereby we investigate the questions whether planetesimals and planetesimal fragments which have gained considerable inclination due to the strong dynamical interactions in the main belt region around 2 AU can be efficient water transporting vessels. The Hungaria asteroid group is the best example that such scenarios are realistic. Assuming that the gas giants and the terrestrial planets are already formed, we monitor the collisions of scattered small bodies containing water (in the order of a few percent) with the terrestrial planets. Thus we are able to give a first estimate concerning the respective contribution of such bodies to the actual water content in the crust of the Earth.

Locating planets in HabitableZones (HZs) around other stars is a growing field in contemporary astronomy. Since a large percentage of all G-M stars in the solar neighbourhood are expected to be part of binary or multiple stellar systems, investigations of whether habitable planets are likely to be discovered in such environments are of prime interest to the scientific community. As current exoplanet statistics predicts that the chances are higher to find new worlds in systems that are already known to have planets, we examine four known extrasolar planetary systems in tight binaries in order to determine their capacity to host additional habitable terrestrial planets. Those systems are Gliese 86, gamma Cephei, HD 41004 and HD 196885. In the case of gamma Cephei, our results suggest that only the M dwarf companion could host additional potentially habitable worlds. Neither could we identify stable, potentially habitable regions around HD 196885 A. HD 196885 B can be considered a slightly more promising target in the search for Earth-twins. Gliese 86 A turned out to be a very good candidate, assuming that the systems history has not been excessively violent. For HD 41004, we have identified admissible stable orbits for habitable planets, but those strongly depend on the parameters of the system. A more detailed investigation shows that for some initial conditions stable planetary motion is possible in the HZ of HD 41004 A. In spite of the massive companion HD 41004 Bb, we found that HD 41004 B, too, could host additional habitable worlds.

By now, observations of exoplanets have found more than 50 binary star systems hosting 71 planets. We expect these numbers to increase as more than 70% of the main sequence stars in the solar neighborhood are members of binary or multiple systems. The planetary motion in such systems depends strongly on both the parameters of the stellar system (stellar separation and eccentricity) and the architecture of the planetary system (number of planets and their orbital behaviour). In case a terrestrial planet moves in the so-called habitable zone (HZ) of its host star, the habitability of this planet depends on many parameters. A crucial factor is certainly the amount of water. We investigate in this work the transport of water from beyond the snow-line to the HZ in a binary star system and compare it to a single star system.

From the numerous detected planets outside the Solar system, no terrestrial planet comparable to our Earth has been discovered so far. The search for an Exo-Earth is certainly a big challenge which may require the detections of planetary systems resembling our Solar system in order to find life like on Earth. However, even if we find Solar system analogues, it is not certain that a planet in Earth position will have similar circumstances as those of Earth. Small changes in the architecture of the giant planets can lead to orbital perturbations which may change the conditions of habitability for a terrestrial planet in the habitable zone (HZ). We present a numerical investigation where we first study the motion of test-planets in a particular Jupiter-Saturn configuration for which we can expect strong gravitational perturbations on the motion at Earth position according to a previous work. In this study, we show that these strong perturbations can be reduced significantly by the neighboring planets of Earth. In the second part of our study we investigate the motion of test-planets in inclined Jupiter-Saturn systems where we analyze changes in the dynamical behavior of the inner planetary system. Moderate values of inclination seem to counteract the perturbations in the HZ while high inclinations induce more chaos in this region. Finally, we carry out a stability study of the actual orbits of Venus, Earth and Mars moving in the inclined Jupiter-Saturn systems for which we used the Solar system parameters. This study shows that the three terrestrial planets will only move in low-eccentric orbits if Saturn's inclination is <=10{\deg}. Therefore, it seems that it is advantageous for the habitability of Earth when all planets move nearly in the same plane.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Context. In binary star systems, the winds from the two components impact each other, leading to strong shocks and regions of enhanced density and temperature. Potentially habitable circumbinary planets must continually be exposed to these interactions regions. Aims. We study, for the first time, the interactions between winds from low-mass stars in a binary system, to show the wind conditions seen by potentially habitable circumbinary planets. Methods. We use the advanced 3D numerical hydrodynamic code Nurgush to model the wind interactions of two identical winds from two solar mass stars with circular orbits and a binary separation of 0.5 AU. As input into this model, we use a 1D hydrodynamic simulation of the solar wind, run using the Versatile Advection Code. We derive the locations of stable and habitable orbits in this system to explore what wind conditions potentially habitable planets will be exposed to during their orbits. Results. Our wind interaction simulations result in the formation of two strong shock waves separated by a region of enhanced density and temperature. The wind-wind interaction region has a spiral shape due to Coriolis forces generated by the orbital motions of the two stars. The stable and habitable zone in this system extends from approximately 1.4 AU to 2.4 AU. (TRUNCATED)

Determining planetary habitability is a complex matter, as the interplay between a planet's physical and atmospheric properties with stellar insolation has to be studied in a self consistent manner. Standardized atmospheric models for Earth-like planets exist and are commonly accepted as a reference for estimates of Habitable Zones. In order to define Habitable Zone boundaries, circular orbital configurations around main sequence stars are generally assumed. In gravitationally interacting multibody systems, such as double stars, however, planetary orbits are forcibly becoming non circular with time. Especially in binary star systems even relatively small changes in a planet's orbit can have a large impact on habitability. Hence, we argue that a minimum model for calculating Habitable Zones in binary star systems has to include dynamical interactions.

With the discovery of hundreds of exoplanets and a potentially huge number of Earth-like planets waiting to be discovered, the conditions for their habitability have become a focal point in exoplanetary research. The classical picture of habitable zones primarily relies on the stellar flux allowing liquid water to exist on the surface of an Earth-like planet with a suitable atmosphere. However, numerous further stellar and planetary properties constrain habitability. Apart from "geophysical" processes depending on the internal structure and composition of a planet, a complex array of astrophysical factors additionally determine habitability. Among these, variable stellar UV, EUV, and X-ray radiation, stellar and interplanetary magnetic fields, ionized winds, and energetic particles control the constitution of upper planetary atmospheres and their physical and chemical evolution. Short- and long-term stellar variability necessitates full time-dependent studies to understand planetary habitability at any point in time. Furthermore, dynamical effects in planetary systems and transport of water to Earth-like planets set fundamentally important constraints. We will review these astrophysical conditions for habitability under the crucial aspects of the long-term evolution of stellar properties, the consequent extreme conditions in the early evolutionary phase of planetary systems, and the important interplay between properties of the host star and its planets.

We investigate the influence of impacts of large planetesimals and small planetary embryos on the early Martian surface on the hydrodynamic escape of an early steam atmosphere that is exposed to the high soft X-ray and EUV flux of the young Sun. Impact statistics in terms of number, masses, velocities, and angles of asteroid impacts onto the early Mars are determined via n-body integrations. Based on these statistics, smoothed particle hydrodynamics (SPH) simulations result in estimates of energy transfer into the planetary surface material and according surface heating. For the estimation of the atmospheric escape rates we applied a soft X-ray and EUV absorption model and a 1-D upper atmosphere hydrodynamic model to a magma ocean-related catastrophically outgassed steam atmosphere with surface pressure values of 52 bar H2O and 11 bar CO2. The estimated impact rates and energy deposition onto an early Martian surface can account for substantial heating. The energy influx and conversion rate into internal energy is most likely sufficient to keep a shallow magma ocean liquid for an extended period of time. Higher surface temperatures keep the outgassed steam atmosphere longer in vapor form and therefore enhance its escape to space within about 0.6 Myr after its formation.

Given the considerable percentage of stars that are members of binaries or stellar multiples in the Solar neighborhood, it is expected that many of these binaries host planets, possibly even habitable ones. The discovery of a terrestrial planet in the alpha Centauri system supports this notion. Due to the potentially strong gravitational interaction that an Earth-like planet may experience in such systems, classical approaches to determining habitable zones, especially in close S-Type binary systems, can be rather inaccurate. Recent progress in this field, however, allows to identify regions around the star permitting permanent habitability. While the discovery of alpha Cen Bb has shown that terrestrial planets can be detected in solar-type binary stars using current observational facilities, it remains to be shown whether this is also the case for Earth analogues in habitable zones. We provide analytical expressions for the maximum and RMS values of radial velocity and astrometric signals, as well as transit probabilities of terrestrial planets in such systems, showing that the dynamical interaction of the second star with the planet may indeed facilitate the planets detection. As an example, we discuss the detectability of additional Earth-like planets in the averaged, extended, and permanent habitable zones around both stars of the alpha Centauri system.

Binary and multiple systems constitute more than half of the total stellar population in the Solar neighborhood (Kiseleva-Eggleton and Eggleton 2001). Their frequent occurrence as well as the fact that more than 70 (Schneider et al. 2011) planets have already been discovered in such configurations - most noteably the telluric companion of alpha Centauri B (Dumusque et al. 2012) - make them interesting targets in the search for habitable worlds. Recent studies (Eggl et al. 2012b, Forgan 2012) have shown, that despite the variations in gravitational and radiative environment, there are indeed circumstellar regions where planets can stay within habitable insolation limits on secular dynamical timescales. In this article we provide habitable zones for 19 near S-Type binary systems from the Hipparchos and WDS catalogues with semimajor axes between 1 and 100 AU. Hereby, we accounted for the combined dynamical and radiative influence of the second star on the Earth-like planet. Out of the 19 systems presented, 17 offer dynamically stable habitable zones around at least one component. The 17 potentially habitable systems contain 5 F, 3 G, 7 K and 16 M class stars. As their proximity to the Solar System (d < 31 pc) makes the selected binary stars exquisite targets for observational campaigns, we offer estimates on radial velocity, astrometric and transit signatures produced by habitable Earth-like planets in eccentric circumstellar orbits.

With more and more extrasolar planets discovered in and around binary star systems, questions concerning the determination of the classical Habitable Zone arise. Do the radiative and gravitational perturbations of the second star influence the extent of the Habitable Zone significantly, or is it sufficient to consider the host-star only? In this article we investigate the implications of stellar companions with different spectral types on the insolation a terrestrial planet receives orbiting a Sun-like primary. We present time independent analytical estimates and compare these to insolation statistics gained via high precision numerical orbit calculations. Results suggest a strong dependence of permanent habitability on the binary's eccentricity, as well as a possible extension of Habitable Zones towards the secondary in close binary systems.

In close eclipsing binaries, measurements of the variations in binary's eclipse timing may be used to infer information about the existence of circumbinary objects. To determine the possibility of the detection of such variations with CoRoT and Kepler space telescopes, we have carried out an extensive study of the dynamics of a binary star system with a circumbinary planet, and calculated its eclipse timing variations (ETV) for different values of the mass-ratio and orbital elements of the binary and the perturbing body. Here, we present the results of our study and assess the detectability of the planet by comparing the resulting values of ETVs with the temporal sensitivity of CoRoT and Kepler. Results point to extended regions in the parameter-space where the perturbation of a planet may become large enough to create measurable variations in the eclipse timing of the secondary star. Many of these variations point to potentially detectable ETVs and the possible existence of Jovian-type planets.

It is known since the seminal study of Laskar (1989) that the inner planetary system is chaotic with respect to its orbits and even escapes are not impossible, although in time scales of billions of years. The aim of this investigation is to locate the orbits of Venus and Earth in phase space, respectively to see how close their orbits are to chaotic motion which would lead to unstable orbits for the inner planets on much shorter time scales. Therefore we did numerical experiments in different dynamical models with different initial conditions -- on one hand the couple Venus-Earth was set close to different mean motion resonances (MMR), and on the other hand Venus' orbital eccentricity (or inclination) was set to values as large as e = 0.36 (i = 40deg). The couple Venus-Earth is almost exactly in the 13:8 mean motion resonance. The stronger acting 8:5 MMR inside, and the 5:3 MMR outside the 13:8 resonance are within a small shift in the Earth's semimajor axis (only 1.5 percent). Especially Mercury is strongly affected by relatively small changes in eccentricity and/or inclination of Venus in these resonances. Even escapes for the innermost planet are possible which may happen quite rapidly.

Whether binaries can harbor potentially habitable planets depends on several factors including the physical properties and the orbital characteristics of the binary system. While the former determines the location of the habitable zone (HZ), the latter affects the dynamics of the material from which terrestrial planets are formed (i.e., planetesimals and planetary embryos), and drives the final architecture of the planets assembly. In order for a habitable planet to form in a binary star system, these two factors have to work in harmony. That is, the orbital dynamics of the two stars and their interactions with the planet-forming material have to allow terrestrial planet formation in the habitable zone, and ensure that the orbit of a potentially habitable planet will be stable for long times. We have organized this chapter with the same order in mind. We begin by presenting a general discussion on the motion of planets in binary stars and their stability. We then discuss the stability of terrestrial planets, and the formation of potentially habitable planets in a binary-planetary system.

We investigate the stability regions of hypothetical terrestrial planets around the Lagrangian equilibrium points L4 and L5 in some specific extrasolar planetary systems. The problem of their stability can be treated in the framework of the restricted three body problem where the host star and a massive Jupiter-like planet are the primary bodies and the terrestrial planet is regarded as being massless. From these theoretical investigations one cannot determine the extension of the stable zones around the equilibrium points. Using numerical experiments we determined their largeness for three test systems chosen from the table of the know extrasolar planets, where a giant planet is moving close to the so-called habitable zone around the host star in low eccentric orbits. The results show the dependence of the size and structure of this region, which shrinks significantly with the eccentricity of the known gas giant.

The orbits in the HD 160691 planetary system at first appeared highly unstable, but using the MEGNO and FLI techniques of global dynamics analysis in the orbital parameter space we have found a stabilizing mechanism that could be the key to its existence. In order to be dynamically stable, the HD 160691 planetary system has to satisfy the following conditions: (1) a 2:1 mean motion resonance, combined with (2) an apsidal secular resonance in (3) a configuration $P_{c}(ap) - S - P_{b}(ap)$ where the two apsidal lines are anti-aligned, and (4) specific conditions on the respective sizes of the eccentricities (high eccentricity for the outer orbit is in particular the most probable necessary condition). More generally, in this original orbital topology, where the resonance variables $\theta_{1}$ and $\theta_{3}$ librate about $180^{\circ}$ while $\theta_{2}$ librates about $0^{\circ}$, the HD 160691 system and its mechanism have revealed aspects of the 2:1 orbital resonances that have not been observed nor analyzed before. The present topology combined with the 2:1 resonance is indeed more wide-ranging than the particular case of the HD 160691 planetary system. It is a new theoretical possibility suitable for a stable regime despite relatively small semi-major axes with respect to the important masses in interactions.

Using numerical methods we thoroughly investigate the dynamical stability in the region between the two planets found in HD 74156. The two planets with minimum masses 1.56 M_JUP (HD 74156b) and 7.5 M_JUP (HD 74156c), semimajor axes 0.276 AU and 3.47 AU move on quite eccentric orbits (e=0.649 and 0.395). There is a region between 0.7 and 1.4 AU which may host additional planets which we checked via numerical integrations using different dynamical models. Besides the orbital evolution of several thousands of massless regarded planets in a three-dimensional restricted 4-body problem (host star, two planets + massless bodies) we also have undertaken test computation for the orbital evolution for fictive planets with masses of 0.1, 0.3 and 1 M_JUP in the region between HD74156b and HD74156c. For direct numerical integrations up to 10^7 years we used the Lie-integrator, a method with adaptive stepsize; additionally we used the Fast Lyapunov Indicators as tool for detecting chaotic motion in this region. We emphasize the important role of the inner resonances (with the outer planet) and the outer resonances (with the inner planet) with test bodies located inside the resonances. In these two "resonance" regions almost no orbits survive. The region between the 1:5 outer resonance (0.8 AU) and the 5:1 inner resonance (1.3 AU), just in the right position for habitability, is also very unstable probably due to three-body-resonances acting there. Our results do not strictly "forbid" planets to move there, but the existence of a planet on a stable orbit between 0.8 and 1.3 AU is unlikely.

The recently discovered planetary system in the binary GamCep was studied concerning its dynamical evolution. We confirm that the orbital parameters found by the observers are in a stable configuration. The primary aim of this study was to find stable planetary orbits in a habitable region in this system, which consists of a double star (a=21.36 AU) and a relatively close (a=2.15 AU) massive (1.7 Mjup sin i) planet. We did straightforward numerical integrations of the equations of motion in different dynamical models and determined the stability regions for a fictitious massless planet in the interval of the semimajor axis 0.5 AU < a < 1.85 AU around the more massive primary. To confirm the results we used the Fast Lyapunov Indicators (FLI) in separate computations, which are a common tool for determining the chaoticity of an orbit. Both results are in good agreement and unveiled a small island of stable motions close to 1 AU up to an inclination of about 15 deg (which corresponds to the 3:1 mean motion resonance between the two planets). Additionally we computed the orbits of earthlike planets (up to 90 earthmasses) in the small stable island and found out, that there exists a small window of stable orbits on the inner edge of the habitable zone in GamCep even for massive planets.