We present results of a detailed study of the rate of the accretion of planetesimals by a growing proto-Jupiter in the core-accretion model. Using a newly developed code, we accurately combine a detailed three-body trajectory calculation with gas drag experienced during the passage of planetesimals in the protoplanet's envelope. We find that the motion of planetesimals is excited to the extent that encounters with the proto-planetary envelope become so fast that ram pressure breaks up the planetesimals in most encounters. As a result, the accretion rate is largely independent of the planetesimal size and composition. For the case we explored of a planet forming at 5.2 AU from the Sun in a disk with a solid surface density of 6 g/cm^2 (Lozovsky et al. 2017) the accretion rate we compute differs in several respects from that assumed by those authors. We find that only 4-5 M_Earth is accreted in the first 1.5x10^6 years before the onset of rapid gas accretion. Most of the mass, some 10 M_Earth, is accreted simultaneously with this rapid gas accretion. In addition, we find that the mass accretion rate remains small, but non-zero for at least a million years after this point, and an additional 0.3-0.4 M_Earth is accreted during that time. This late accretion, together with a rapid infall of gas could lead to the accreted material being mixed throughout the outer regions, and may account for the enhancement of high-Z material in Jupiter's envelope.

We report the detection of the first circumbinary planet found by TESS. The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30-minute cadence and in sectors 4 through 12 at two-minute cadence. It consists of two stars with masses of 1.1 MSun and 0.3 MSun on a slightly eccentric (0.16), 14.6-day orbit, producing prominent primary eclipses and shallow secondary eclipses. The planet has a radius of ~6.9 REarth and was observed to make three transits across the primary star of roughly equal depths (~0.2%) but different durations -- a common signature of transiting circumbinary planets. Its orbit is nearly circular (e ~ 0.09) with an orbital period of 95.2 days. The orbital planes of the binary and the planet are aligned to within ~1 degree. To obtain a complete solution for the system, we combined the TESS photometry with existing ground-based radial-velocity observations in a numerical photometric-dynamical model. The system demonstrates the discovery potential of TESS for circumbinary planets, and provides further understanding of the formation and evolution of planets orbiting close binary stars.

Planets that orbit only one of the stars in stellar binary systems (i.e., circumstellar) are dynamically constrained to a limited range of orbital parameters and thus understanding conditions on their stability is of great importance in exoplanet searches. We perform $\sim$700 million N-body simulations to identify how stability regions depend on properties of the binary, as well as, the starting planetary inclination and mean longitude relative to the binary orbit. Moreover, we provide grid interpolation maps and lookup tables for the community to use our results. Through Monte-Carlo methods we determine that planets with a semimajor axis $a_p$ $\lesssim$8\% of the binary semimajor axis $a_{bin}$ will likely be stable given the known distribution of binary star parameters. This estimate varies in the Lidov-Kozai regime or for retrograde orbits to 4\% or 10\% of $a_{bin}$, respectively. Our method to quickly determine the circumstellar stability limit is important for interpreting observations of binaries using direct imaging with JWST, photometry with TESS, or even astrometry with Gaia.

To date a dozen transiting "Tatooines" or circumbinary planets (CBPs) have been discovered, by eye, in the data from the Kepler mission; by contrast, thousands of confirmed circumstellar planets orbiting around single stars have been detected using automated algorithms. Automated detection of CBPs is challenging because their transits are strongly aperiodic with irregular profiles. Here, we describe an efficient and automated technique for detecting circumbinary planets that transit their binary hosts in Kepler light curves. Our method accounts for large transit timing and duration variations (TTVs and TDVs), induced by binary reflex motion, in two ways: 1) We directly correct for large-scale TTVs and TDVs in the light curves by using Keplerian models to approximate binary and CBP orbits; and 2) We allow additional aperiodicities on the corrected light curves by employing the Quasi-periodic Automated Transit Search algorithm (QATS). We demonstrate that our method dramatically improves detection significance using simulated data and two previously identified CBP systems, Kepler-35 and Kepler-64.

We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets.

Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the "middle planet"), which has an orbital period of 187.4 days and a radius of about $7\,R_{\oplus}$. The presence of the middle planet allows us to place much better constraints on the masses of all three planets, where the $1\sigma$ ranges are less than $26\,M_{\oplus}$, between $7-43\,M_{\oplus}$, and between $2-5\,M_{\oplus}$ for the inner, middle, and outer planets, respectively. The middle and outer planets have low bulk densities, with $\rho_{\rm middle} < 0.68$ g cm$^{-3}$ and $\rho_{\rm outer} < 0.26$ g cm$^{-3}$ at the $1\sigma$ level. The two outer planets are "tightly packed," assuming the nominal masses, meaning no other planet could stably orbit between them. All of the orbits have low eccentricities and are nearly coplanar, disfavoring violent scattering scenarios and suggesting gentle migration in the protoplanetary disk.

While recently discovered exotic new planet-types have both challenged our imaginations and broadened our knowledge of planetary system workings, perhaps the most compelling objective of exoplanet science is to detect and characterize habitable and possibly inhabited worlds orbiting in other star systems. For the foreseeable future, characterizations of extrasolar planets will be made via remote sensing of planetary spectroscopic and temporal signals, along with careful fitting of this data to advanced models of planets and their atmospheres. Terrestrial planets are small and significantly more challenging to observe compared to their larger gaseous brethren; however observatories coming on-line in the coming decade will begin to allow their characterization. Still, it is not enough to invest only in observational endeavors. Comprehensive modeling of planetary atmospheres is required in order to fully understand what it is that our grand telescopes see in the night-sky. In our quest to characterize habitable, and possibly inhabited worlds, 3D general circulation models (GCMs) should be used to evaluate potential climate states and their associated temporal and spatial dependent observable signals. 3D models allow for coupled, self-consistent, multi-dimensional simulations, which can realistically simulate the climates of terrestrial extrasolar planets. A complete theoretical understanding of terrestrial exoplanetary atmospheres, gained through comprehensive 3D modeling, is critical for interpreting spectra of exoplanets taken from current and planned instruments, and is critical for designing future missions that aim to measure spectra of potentially habitable exoplanets as one of their key science goals. We recommend continued institutional support for 3D GCM modeling teams that focus on planetary and exoplanetary applications.

The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the physical and chemical conditions in protoplanetary disks using infrared spectroscopy and (2)~studying planet-disk interactions using imaging and spectro-astrometry. Our current knowledge is based on a limited sample of targets, representing the brightest, most extreme cases, and thus almost certainly represents an incomplete understanding. ELTs will play a transformational role in this arena, thanks to the high spatial and spectral resolution data they will deliver. We recommend a key science program to conduct a volume-limited survey of high-resolution spectroscopy and high-contrast imaging of the nearest protoplanetary disks that would result in an unbiased, holistic picture of planet formation as it occurs.

We present a community-led assessment of the solar system investigations achievable with NASA's next-generation space telescope, the Wide Field InfraRed Survey Telescope (WFIRST). WFIRST will provide imaging, spectroscopic, and coronagraphic capabilities from 0.43-2.0 $\mu$m and will be a potential contemporary and eventual successor to JWST. Surveys of irregular satellites and minor bodies are where WFIRST will excel with its 0.28 deg$^2$ field of view Wide Field Instrument (WFI). Potential ground-breaking discoveries from WFIRST could include detection of the first minor bodies orbiting in the Inner Oort Cloud, identification of additional Earth Trojan asteroids, and the discovery and characterization of asteroid binary systems similar to Ida/Dactyl. Additional investigations into asteroids, giant planet satellites, Trojan asteroids, Centaurs, Kuiper Belt Objects, and comets are presented. Previous use of astrophysics assets for solar system science and synergies between WFIRST, LSST, JWST, and the proposed NEOCam mission are discussed. We also present the case for implementation of moving target tracking, a feature that will benefit from the heritage of JWST and enable a broader range of solar system observations.

The stability limit for circumbinary planets (CBPs) is not well defined and can depend on initial parameters defining either the planetary orbit or the inner binary orbit. We expand on the work of Holman & Wiegert (1999, AJ 117, 621) to develop numerical tools for quick, easy, and accurate determination of the stability limit. The results of our simulations, as well as our numerical tools, are available to the community through $\texttt{Zenodo}$ and $\texttt{GitHub}$, respectively. We employ a grid interpolation method based on $\sim$150 million full N-body simulations of initially circular, coplanar systems and compare to the 9 known Kepler CBP systems. Using a formalism from planet packing studies, we find that 55% of the Kepler CBP systems allow for an additional equal-mass planet to potentially exist on an interior orbit relative to the observed planet. Therefore, we do $\textit{not}$ find strong evidence for a pile-up in the Kepler CBP systems and more detections are needed to adequately characterize the formation mechanisms for the CBP population. Observations from the Transiting Exoplanet Survey Satellite are expected to substantially increase the number of detections using the unique geometry of CBP systems, where multiple transits can occur during a single conjunction.

It has been suggested that the comet-like activity of Main Belt Comets is due to the sublimation of sub-surface water-ice that is exposed when these objects are impacted by meter-sized bodies. We recently examined this scenario and showed that such impacts can in fact excavate ice and present a plausible mechanism for triggering the activation of MBCs (Haghighipour et al. 2016). However, because the purpose of that study was to prove the concept and identify the most viable ice-longevity model, the porosity of the object and the loss of ice due to the heat of impact were ignored. In this paper, we extend our impact simulations to porous materials and account for the loss of ice due to an impact. We show that for a porous MBC, impact craters are deeper, reaching to approximately 15 m implying that if the activation of MBCs is due to the sublimation of sub-surface ice, this ice has to be within the top 15 m of the object. Results also indicate that the loss of ice due to the heat of impact is negligible, and the re-accretion of ejected ice is small. The latter suggests that the activities of current MBCs are most probably from multiple impact sites. Our study also indicates that in order for sublimation from multiple sites to account for the observed activity of the currently known MBCs, the water content of MBCs (and their parent asteroids) needs to be larger than the values traditionally considered in models of terrestrial planet formation.

The NASA Kepler mission has discovered thousands of new planetary candidates, many of which have been confirmed through follow-up observations. A primary goal of the mission is to determine the occurrance rate of terrestrial-size planets within the Habitable Zone (HZ) of their host stars. Here we provide a list of HZ exoplanet candidates from the Kepler Data Release 24 Q1-Q17 data vetting process. This work was undertaken as part of the Kepler Habitable Zone Working Group. We use a variety of criteria regarding HZ boundaries and planetary sizes to produce complete lists of HZ candidates, including a catalog of 104 candidates within the optimistic HZ and 20 candidates with radii less than two Earth radii within the conservative HZ. We cross-match our HZ candidates with the Data Release 25 stellar properties and confirmed planet properties to provide robust stellar parameters and candidate dispositions. We also include false positive probabilities recently calculated by Morton et al. (2016) for each of the candidates within our catalogs to aid in their validation. Finally, we performed dynamical analysis simulations for multi-planet systems that contain candidates with radii less than two Earth radii as a step toward validation of those systems.

It has been suggested that the comet-like activity of Main Belt Comets are due to the sublimation of sub-surface water-ice that has been exposed as a result of their surfaces being impacted by m-sized bodies. We have examined the viability of this scenario by simulating impacts between m-sized and km-sized objects using a smooth particle hydrodynamics approach. Simulations have been carried out for different values of the impact velocity and impact angle as well as different target material and water-mass fraction. Results indicate that for the range of impact velocities corresponding to those in the asteroid belt, the depth of an impact crater is slightly larger than 10 m suggesting that if the activation of MBCs is due to the sublimation of sub-surface water-ice, this ice has to exist no deeper than a few meters from the surface. Results also show that ice-exposure occurs in the bottom and on the interior surface of impact craters as well as the surface of the target where some of the ejected icy inclusions are re-accreted. While our results demonstrate that the impact scenario is indeed a viable mechanism to expose ice and trigger the activity of MBCs, they also indicate that the activity of the current MBCs is likely due to ice sublimation from multiple impact sites and/or the water contents of these objects (and other asteroids in the outer asteroid belt) is larger than the 5% that is traditionally considered in models of terrestrial planet formation providing more ice for sublimation. We present details of our simulations and discuss their results and implications.

We report the discovery of a new Kepler transiting circumbinary planet (CBP). This latest addition to the still-small family of CBPs defies the current trend of known short-period planets orbiting near the stability limit of binary stars. Unlike the previous discoveries, the planet revolving around the eclipsing binary system Kepler-1647 has a very long orbital period (~1100 days) and was at conjunction only twice during the Kepler mission lifetime. Due to the singular configuration of the system, Kepler-1647b is not only the longest-period transiting CBP at the time of writing, but also one of the longest-period transiting planets. With a radius of 1.06+/-0.01 RJup it is also the largest CBP to date. The planet produced three transits in the light-curve of Kepler-1647 (one of them during an eclipse, creating a syzygy) and measurably perturbed the times of the stellar eclipses, allowing us to measure its mass to be 1.52+/-0.65 MJup. The planet revolves around an 11-day period eclipsing binary consisting of two Solar-mass stars on a slightly inclined, mildly eccentric (e_bin = 0.16), spin-synchronized orbit. Despite having an orbital period three times longer than Earth's, Kepler-1647b is in the conservative habitable zone of the binary star throughout its orbit.

We present the results of "snapshot" numerical integrations of test particles representing comet-like and asteroid-like objects in the inner solar system aimed at investigating the short-term dynamical evolution of objects close to the dynamical boundary between asteroids and comets as defined by the Tisserand parameter with respect to Jupiter, T_J (i.e., T_J=3). As expected, we find that T_J for individual test particles is not always a reliable indicator of initial orbit types. Furthermore, we find that a few percent of test particles with comet-like starting elements (i.e., similar to Jupiter-family comets) reach main-belt-like orbits (at least temporarily) during our 2 Myr integrations, even without the inclusion of non-gravitational forces, apparently via a combination of gravitational interactions with the terrestrial planets and temporary trapping by mean-motion resonances with Jupiter. We estimate that the fraction of real Jupiter-family comets occasionally reaching main-belt-like orbits on Myr timescales could be on the order of ~0.1-1%, although the fraction that remain on such orbits for appreciable lengths of time is certainly far lower. Thus, the number of JFC-like interlopers in the main-belt population at any given time is likely to be small, but still non-zero, a finding with significant implications for efforts to use main-belt comets to trace the primordial distribution of volatile material in the inner solar system. The test particles with comet-like starting orbital elements that transition onto main-belt-like orbits in our integrations appear to be largely prevented from reaching low eccentricity, low inclination orbits. We therefore find that low-eccentricity, low-inclination main-belt comets may provide a more reliable means for tracing the primordial ice content of the main asteroid belt than the main-belt comet population as a whole.

Spurred by the discovery of numerous exoplanets in multiple systems, binaries have become in recent years one of the main topics in planet formation research. Numerous studies have investigated to what extent the presence of a stellar companion can affect the planet formation process. Such studies have implications that can reach beyond the sole context of binaries, as they allow to test certain aspects of the planet formation scenario by submitting them to extreme environments. We review here the current understanding on this complex problem. We show in particular how each of the different stages of the planet-formation process is affected differently by binary perturbations. We focus especially on the intermediate stage of kilometre-sized planetesimal accretion, which has proven to be the most sensitive to binarity and for which the presence of some exoplanets observed in tight binaries is difficult to explain by in-situ formation following the "standard" planet-formation scenario. Some tentative solutions to this apparent paradox are presented. The last part of our review presents a thorough description of the problem of planet habitability, for which the binary environment creates a complex situation because of the presence of two irradation sources of varying distance.

We present the results of an extensive study of the final stage of terrestrial planet formation in disks with different surface density profiles and for different orbits of Jupiter and Saturn. We carried out simulations for disk densities proportional to r^-0.5, r^-1, and r^-1.5, and also for partially depleted disks as in the recent model of Mars formation by Izidoro et al (2014). The purpose of our study is to determine how the final assembly of planets and their physical properties are affected by the total mass of the disk and its radial profile. Because of the important roles of secular resonances in orbits and properties of the final planets, we studied the effects of these resonances as well. We have divided this study into two parts. In Part 1, we are interested in examining the effects of secular resonances on the formation of Mars and orbital stability of terrestrial planets. In Part 2, our goal is to determine trends that may exist between the disk surface density profile and the final properties of terrestrial planets. In the context of the depleted disk model, results show that the nu_5 resonance does not have a significant effect on the final orbits of terrestrial planets. However, nu_6 and nu_16 resonances play important roles in clearing their affected areas ensuring that no additional mass will be scattered into the accretion zone of Mars so that it can maintain its mass and orbital stability. In Part 2, our results indicate that despite some small correlations, in general, no trend seems to exist between the disk surface density profile and the mean number of the final planets, their masses, time of formation, and distances to the central star. We present the results of our simulations and discuss their implications for the formation of Mars and other terrestrial planets, as well as the physical properties of these objects such as their masses and water contents.

The existence of planets orbiting a central binary star system immediately raises questions regarding their formation and dynamical evolution. Recent discoveries of circumbinary planets by the Kepler space telescope has shown that some of these planets reside close to the dynamical stability limit where it is very difficult to form planets in situ. For binary systems with nearly circular orbits, such as Kepler-38, the observed proximity of planetary orbits to the stability limit can be understood by an evolutionary process in which planets form farther out in the disk and migrate inward to their observed position. The Kepler-34 system has a high orbital eccentricity of 0.52. Here, we analyse evolutionary scenarios for the planet observed around this system using two-dimensional hydrodynamical simulations. The highly eccentric binary opens a wide inner hole in the disk which is also eccentric, and displays a slow prograde precession. As a result of the large, eccentric inner gap, an embedded planet settles in a final equilibrium position that lies beyond the observed location of Kepler-34 b, but has the correct eccentricity. In this configuration the planetary orbit is aligned with the disk in a state of apsidal corotation.To account for the closer orbit of Kepler-34 b to the central binary, we considered a two-planet scenario and examined the evolution of the system through joint inward migration and capture into mean-motion resonances. When the inner planet orbits inside the gap of the disk, planet-planet scattering ensues. While often one object is thrown into a large, highly eccentric orbit, at times the system is left with a planet close to the observed orbit, suggesting that Kepler 34 might have had two circumbinary planets where one might have been scattered out of the system or into an orbit where it did not transit the central binary during the operation of Kepler.

We present an observational and dynamical study of newly discovered main-belt comet 313P/Gibbs. We find that the object is clearly active both in observations obtained in 2014 and in precovery observations obtained in 2003 by the Sloan Digital Sky Survey, strongly suggesting that its activity is sublimation-driven. This conclusion is supported by a photometric analysis showing an increase in the total brightness of the comet over the 2014 observing period, and dust modeling results showing that the dust emission persists over at least three months during both active periods, where we find start dates for emission no later than 2003 July 24+/-10 for the 2003 active period and 2014 July 28+/-10 for the 2014 active period. From serendipitous observations by the Subaru Telescope in 2004 when the object was apparently inactive, we estimate that the nucleus has an absolute R-band magnitude of H_R=17.1+/-0.3, corresponding to an effective nucleus radius of r_e~1.00+/-0.15 km. The object's faintness at that time means we cannot rule out the presence of activity, and so this computed radius should be considered an upper limit. We find that 313P's orbit is intrinsically chaotic, having a Lyapunov time of T_l=12000 yr and being located near two 3-body mean-motion resonances with Jupiter and Saturn, 11J-1S-5A and 10J+12S-7A, yet appears stable over >50 Myr in an apparent example of stable chaos. We furthermore find that 313P is the second main-belt comet, after P/2012 T1 (PANSTARRS), to belong to the ~155 Myr old Lixiaohua asteroid family.

We have studied the possibility that a third circumbinary planet in the Kepler-47 planetary system be the source of the single unexplained transiting event reported during the discovery of these planets. We applied the MEGNO technique to identify regions in the phase space where a third planet can maintain quasi-periodic orbits, and assessed the long-term stability of the three-planet system by integrating the entire 5 bodies (binary + planets) for 10 Myr. We identified several stable regions between the two known planets as well as a region beyond the orbit of Kepler-47c where the orbit of the third planet could be stable. To constrain the orbit of this planet, we used the measured duration of the unexplained transit event ($\sim 4.15$ hours) and compared that with the transit duration of the third planet in an ensemble of stable orbits. To remove the degeneracy among the orbits with similar transit durations, we considered the planet to be in a circular orbit and calculated its period analytically. The latter places an upper limit of 424 days on the orbital period of the third planet. Our analysis suggests that if the unexplained transit event detected during the discovery of the Kepler-47 circumbinary system is due to a planetary object, this planet will be in a low eccentricity orbit with a semimajor axis smaller than 1.24 AU. Further constraining of the mass and orbital elements of this planet requires a re-analysis of the entire currently available data including those obtained post-announcement of the discovery of this system. We present details of our methodology and discuss the implication of the results.

We present the discovery of KIC 9632895b, a 6.2 Earth-radius planet in a low-eccentricity, 240.5-day orbit about an eclipsing binary. The binary itself consists of a 0.93 and 0.194 solar mass pair of stars with an orbital period of 27.3 days. The plane of the planet's orbit is rapidly precessing, and its inclination only becomes sufficiently aligned with the primary star in the latter portion of the Kepler data. Thus three transits are present in the latter half of the light curve, but none of the three conjunctions that occurred during the first half of the light curve produced transits. The precession period is ~103 years, and during that cycle, transits are visible only ~8% of the time. This has the important implication that for every system like KIC 9632895 that we detect, there are ~12 circumbinary systems that exist but are not currently exhibiting transits. The planet's mass is too small to noticeably perturb the binary, consequently its mass is not measurable with these data; but our photodynamical model places a 1-sigma upper limit of 16 Earth masses. With a period 8.8 times that of the binary, the planet is well outside the dynamical instability zone. It does, however, lie within the habitable zone of the binary, and making it the third of ten Kepler circumbinary planets to do so.

Transmission spectroscopy during planetary transits, which is based on the measurements of the variations of planet-to-star radius ratio as a function of wavelength, is a powerful technique to study the atmospheric properties of transiting planets. One of the main limitation of this technique is the effects of stellar activity, which up until now, have been taken into account only by assessing the effect of non-occulted stellar spots on the estimates of planet-to-star radius ratio. In this paper, we study, for the first time, the impact of the occultation of a stellar spot and plage on the transmission spectra of transiting exoplanets. We simulated this effect by generating a large number of transit light curves for different transiting planets, stellar spectral types, and for different wavelengths. Results of our simulations indicate that the anomalies inside the transit light curve can lead to a significant underestimation or overestimation of the planet-to-star radius ratio as a function of wavelength. At short wavelengths, the effect can reach to a difference of up to 10% in the planet-to-star radius ratio, mimicking the signature of light scattering in the planetary atmosphere. Atmospheric scattering has been proposed to interpret the increasing slopes of transmission spectra toward blue for exoplanets HD 189733b and GJ 3470b. Here we show that these signatures can be alternatively interpreted by the occultation of stellar plages. Results also suggest that the best strategy to identify and quantify the effects of stellar activities on the transmission spectrum of a planet is to perform several observations during the transit epoch at the same wavelength. This will allow for identifying the possible variations in transit depth as a function of time due to stellar activity variability.

We have calculated the coherence and detectable lifetimes of synthetic near-Earth object (NEO) families created by catastrophic disruption of a progenitor as it suffers a very close Earth approach. The closest or slowest approaches yield the most violent `s-class' disruption events. We found that the average slope of the absolute magnitude (H) distribution, $N(H)\propto10^{(0.55\pm0.04)\,H}$, for the fragments in the s-class families is steeper than the slope of the NEO population \citep{mainzer2011} in the same size range. The families remain coherent as statistically significant clusters of orbits within the NEO population for an average of $\bar\tau_c = (14.7\pm0.6)\times10^3$ years after disruption. The s-class families are detectable with the techniques developed by \citet{fu2005} and \citet{Schunova2012} for an average duration ($\bar\tau_{det}$) ranging from about 2,000 to about 12,000 years for progenitors in the absolute magnitude ($H_p$) range from 20 to 13 corresponding to diameters in the range from about 0.5 to 10$\km$ respectively. The short detectability lifetime explains why zero NEO families have been discovered to-date. Nonetheless, every tidal disruption event of a progenitor with D$>0.5\km$ is capable of producing several million fragments in the $1\meter$ to $10\meter$ diameter range that can contribute to temporary local density enhancements of small NEOs in Earth's vicinity. We expect that there are about 1,200 objects in the steady state NEO population in this size range due to tidal disruption assuming that one $1\km$ diameter NEO tidally disrupts at Earth every 2,500 years. These objects may be suitable targets for asteroid retrieval missions due to their Earth-like orbits with corresponding low $v_{\infty}$. The fragments from the tidal disruptions at Earth have $\sim5\times$ the collision probability with Earth compared to the background NEO population.

Recently, a number of planets orbiting binary stars have been discovered by the Kepler space telescope. In a few systems the planets reside close to the dynamical stability limit. Due to the difficulty of forming planets in such close orbits, it is believed that they have formed further out in the disk and migrated to their present locations. Our goal is to construct more realistic models of planet migration in circumbinary disks, and to determine the final position of these planets. In our work, we focus on the system Kepler-38. The evolution of the circumbinary disk is studied using two-dimensional hydrodynamical simulations. We study locally isothermal disks as well as more realistic models with viscous heating, radiative cooling from the disk surfaces, and radiative diffusion in the disk mid plane. After the disk has been brought into equilibrium, a 115 Earth-mass planet is embedded and its evolution is followed. In all cases the planets stop inward migration near the inner edge of the disk. In isothermal disks with a typical disk scale height of H/r = 0.05, the final outcome agrees very well with the observed location of planet Kepler-38b. For the radiative models, the disk thickness and location of the inner edge is determined by the mass in the system. For surface densities in the order of 3000 g/cm^2 at 1 AU, the inner gap lies close to the binary and planets stop in the region between the 5:1 and 4:1 mean-motion resonances with the binary. A model with a disk with approximately a quarter of the mass yields a final position very close to the observed one. For planets migrating in circumbinary disks, the final position is dictated by the structure of the disk. Knowing the observed orbits of circumbinary planets, radiative disk simulations with embedded planets can provide important information on the physical state of the system during the final stages of its evolution.

We have developed a comprehensive methodology and an interactive website for calculating the habitable zone (HZ) of multiple star systems. Using the concept of spectral weight factor, as introduced in our previous studies of the calculations of HZ in and around binary star systems, we calculate the contribution of each star (based on its spectral energy distribution) to the total flux received at the top of the atmosphere of an Earth-like planet, and use the models of the HZ of the Sun to determine the boundaries of the HZ in multiple star systems. Our interactive website for carrying out these calculations is publicly available at http://astro.twam.info/hz . We discuss the details of our methodology and present its application to some of the multiple star systems detected by the Kepler space telescope. We also present the instructions for using our interactive website, and demonstrate its capabilities by calculating the HZ for two interesting analytical solutions of the three-body problem.

Models of terrestrial planet formation for our solar system have been successful in producing planets with masses and orbits similar to those of Venus and Earth. However, these models have generally failed to produce Mars-sized objects around 1.5 AU. The body that is usually formed around Mars' semimajor axis is, in general, much more massive than Mars. Only when Jupiter and Saturn are assumed to have initially very eccentric orbits (e $\sim$ 0.1), which seems fairly unlikely for the solar system, or alternately, if the protoplanetary disk is truncated at 1.0 AU, simulations have been able to produce Mars-like bodies in the correct location. In this paper, we examine an alternative scenario for the formation of Mars in which a local depletion in the density of the protosolar nebula results in a non-uniform formation of planetary embryos and ultimately the formation of Mars-sized planets around 1.5 AU. We have carried out extensive numerical simulations of the formation of terrestrial planets in such a disk for different scales of the local density depletion, and for different orbital configurations of the giant planets. Our simulations point to the possibility of the formation of Mars-sized bodies around 1.5 AU, specifically when the scale of the disk local mass-depletion is moderately high (50-75%) and Jupiter and Saturn are initially in their current orbits. In these systems, Mars-analogs are formed from the protoplanetary materials that originate in the regions of disk interior or exterior to the local mass-depletion. Results also indicate that Earth-sized planets can form around 1 AU with a substantial amount of water accreted via primitive water-rich planetesimals and planetary embryos. We present the results of our study and discuss their implications for the formation of terrestrial planets in our solar system.

We have developed a comprehensive methodology for calculating the circumbinary HZ in planet-hosting P-type binary star systems. We present a general formalism for determining the contribution of each star of the binary to the total flux received at the top of the atmosphere of an Earth-like planet, and use the Sun's HZ to calculate the locations of the inner and outer boundaries of the HZ around a binary star system. We apply our calculations to the Kepler's currently known circumbinary planetary system and show the combined stellar flux that determines the boundaries of their HZs. We also show that the HZ in P-type systems is dynamic and depending on the luminosity of the binary stars, their spectral types, and the binary eccentricity, its boundaries vary as the stars of the binary undergo their orbital motion. We present the details of our calculations and discuss the implications of the results.

We have developed a comprehensive methodology for calculating the boundaries of the habitable zone (HZ) of planet-hosting S-type binary star systems. Our approach is general and takes into account the contribution of both stars to the location and extent of the Binary HZ with different stellar spectral types. We have studied how the binary eccentricity and stellar energy distribution affect the extent of the habitable zone. Results indicate that in binaries where the combination of mass-ratio and orbital eccentricity allows planet formation around a star of the system to proceed successfully, the effect of a less luminous secondary on the location of the primary's habitable zone is generally negligible. However, when the secondary is more luminous, it can influence the extent of the HZ. We present the details of the derivations of our methodology and discuss its application to the Binary HZ around the primary and secondary main sequence stars of an FF, MM, and FM binary, as well as two known planet-hosting binaries alpha Cen AB and HD 196886.

We show that the Kepler spacecraft in two-reaction wheel mode of operation is very well suited for the study of eclipsing binary star systems. Continued observations of the Kepler field will provide the most enduring and long-term valuable science. It will enable the discovery and characterization of eclipsing binaries with periods greater than 1 year - these are the most important, yet least understood binaries for habitable-zone planet background considerations. The continued mission will also enable the investigation of hierarchical multiple systems (discovered through eclipse timing variations), and provide drastically improved orbital parameters for circumbinary planetary systems.

In general, in the studies of transit light-curves and the Rossiter-McLaughlin (RM), the contribution of the planet's gravitational microlensing is neglected. Theoretical studies, have, however shown that the planet's microlensing can affect the transit light-curve and in some extreme cases cause the transit depth to vanish. In this letter, we present the results of our quantitative analysis of microlening on the RM effect. Results indicate that for massive planets in on long period orbits, the planet's microlensing will have considerable contribution to the star's RV measurements. We present the details of our study, and discuss our analysis and results.

We have carried out an extensive study of the possibility of the detection of Earth-mass and super-Earth Trojan planets using transit timing variation method with the Kepler space telescope. We have considered a system consisting of a transiting Jovian-type planet in a short period orbit, and determined the induced variations in its transit timing due to an Earth-mass/super-Earth Trojan planet. We mapped a large section of the phase space around the 1:1 mean-motion resonance and identified regions corresponding to several other mean-motion resonances where the orbit of the planet would be stable. We calculated TTVs for different values of the mass and orbital elements of the transiting and perturbing bodies as well as the mass of central star, and identified orbital configurations of these objects (ranges of their orbital elements and masses) for which the resulted TTVs would be within the range of the variations of the transit timing of Kepler's planetary candidates. Results of our study indicate that in general, the amplitudes of the TTVs fall within the detectable range of timing precision obtained from the Kepler's long-cadence data, and depending on the parameters of the system, their magnitudes may become as large as a few hours. The probability of detection is higher for super-Earth Trojans with slightly eccentric orbits around short-period Jovian-type planets with masses slightly smaller than Jupiter. We present the details of our study and discuss the implications of its results.

Stellar activity features such as spots can create complications in determining planetary parameters through spectroscopic and photometric observations. The overlap of a transiting planet and a stellar spot, for instance, can produce anomalies in the transit light curves that may lead to inaccurate estimation of the transit duration, depth and timing. For instance such inaccuracies can affect the precise derivation of the planet's radius. In this paper, we present the results of a quantitative study on the effects of stellar spots on high precision transit light curves. We show that spot anomalies can lead to the estimate of a planet radius that is 4% smaller than the real value. The effects on the transit duration can also be of the order of 4%, longer or shorter. Depending on the size and distribution of spots, anomalies can also produce transit timing variations with significant amplitudes. For instance, TTVs with signal amplitudes of 200 seconds can be produced when the spot is completely dark and as large as the largest Sun spot. Our study also indicates that the minimum size of a stellar spot with detectable affects on the high precision transit light curve is around 0.03 time the stellar radius for typical Kepler telescope precision. We also show that the strategy of including more free parameters (such as transit depth and duration) in the fitting procedure to measure the transit time of each individual transit will not produce accurate results in the case of active stars.

Super-Earths, objects slightly larger than Earth and slightly smaller than Uranus, have found a special place in exoplanetary science. As a new class of planetary bodies, these objects have challenged models of planet formation at both ends of the spectrum and have triggered a great deal of research on the composition and interior dynamics of rocky planets in connection to their masses and radii. Being relatively easier to detect than an Earth-sized planet at 1 AU around a G star, super-Earths have become the focus of worldwide observational campaigns to search for habitable planets. With a range of masses that allows these objects to retain moderate atmospheres and perhaps even plate tectonics, super-Earths may be habitable if they maintain long-term orbits in the habitable zones of their host stars. Given that in the past two years a few such potentially habitable super-Earths have in fact been discovered, it is necessary to develop a deep understanding of the formation and dynamical evolution of these objects. This article reviews the current state of research on the formation of super-Earths and discusses different models of their formation and dynamical evolution.

We present initial results from observations and numerical analyses aimed at characterizing main-belt comet P/2012 T1 (PANSTARRS). Optical monitoring observations were made between October 2012 and February 2013 using the University of Hawaii 2.2 m telescope, the Keck I telescope, the Baade and Clay Magellan telescopes, Faulkes Telescope South, the Perkins Telescope at Lowell Observatory, and the Southern Astrophysical Research (SOAR) telescope. The object's intrinsic brightness approximately doubles from the time of its discovery in early October until mid-November and then decreases by ~60% between late December and early February, similar to photometric behavior exhibited by several other main-belt comets and unlike that exhibited by disrupted asteroid (596) Scheila. We also used Keck to conduct spectroscopic searches for CN emission as well as absorption at 0.7 microns that could indicate the presence of hydrated minerals, finding an upper limit CN production rate of QCN<1.5x10^23 mol/s, from which we infer a water production rate of QH2O<5x10^25 mol/s, and no evidence of the presence of hydrated minerals. Numerical simulations indicate that P/2012 T1 is largely dynamically stable for >100 Myr and is unlikely to be a recently implanted interloper from the outer solar system, while a search for potential asteroid family associations reveal that it is dynamically linked to the ~155 Myr-old Lixiaohua asteroid family.

One of the most important subjects of debate in the formation of the solar system is the origin of Earth's water. Comets have long been considered as the most likely source of the delivery of water to Earth. However, elemental and isotopic arguments suggest a very small contribution from these objects. Other sources have also been proposed, among which, local adsorption of water vapor onto dust grains in the primordial nebula and delivery through planetesimals and planetary embryos have become more prominent. However, no sole source of water provides a satisfactory explanation for Earth's water as a whole. In view of that, using numerical simulations, we have developed a compound model incorporating both the principal endogenous and exogenous theories, and investigating their implications for terrestrial planet formation and water-delivery. Comets are also considered in the final analysis, as it is likely that at least some of Earth's water has cometary origin. We analyze our results comparing two different water distribution models, and complement our study using D/H ratio, finding possible relative contributions from each source, focusing on planets formed in the habitable zone. We find that the compound model play an important role by showing more advantage in the amount and time of water-delivery in Earth-like planets.

This paper presents new observations of the planet-hosting, visual binary GJ 86 (HR 637) using the Hubble Space Telescope. Ultraviolet and optical imaging with WFC3 confirms the stellar companion is a degenerate star and indicates the binary semimajor axis is larger than previous estimates, with a > 28 AU. Optical STIS spectroscopy of the secondary reveals a helium-rich white dwarf with C2 absorption bands and Teff = 8180 K, thus making the binary system rather similar to Procyon. Based on the 10.8 pc distance, the companion has 0.59 Msun and descended from a main-sequence A star of 1.9 Msun with an original orbital separation a > 14 AU. If the giant planet is coplanar with the binary, the mass of GJ 86Ab is between 4.4 and 4.7 MJup. The similarity of GJ 86 and Procyon prompted a re-analysis of the white dwarf in the latter system, with the tentative conclusion that Procyon hosts a planetesimal population. The periastron distance in Procyon is 20% smaller than in alpha Cen AB, but the metal-enriched atmosphere of Procyon B indicates that the planet formation process minimally attained 25 km bodies, if not small planets as in alpha Cen.

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.

We present an improved version of SOAP (Boisse et al. 2012) named "SOAP-T", which can generate the radial velocity variations and light-curves for systems consisting of a rotating spotted star with a transiting planet. This tool can be used to study the anomalies inside transit light-curves and the Rossiter-McLaughlin effect, to better constrain the orbital configuration and properties of planetary systems and active zones of their host stars. Tests of the code are presented to illustrate its performance and to validate its capability when compared with analytical models and real data. Finally, we apply SOAP-T to the active star, HAT-P-11, observed by the NASA Kepler space telescope and use this system to discuss the capability of this tool in analyzing light-curves for the cases where the transiting planet overlaps with the star's spots.

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.

We report the detection of Kepler-47, a system consisting of two planets orbiting around an eclipsing pair of stars. The inner and outer planets have radii 3.0 and 4.6 times that of the Earth, respectively. The binary star consists of a Sun-like star and a companion roughly one-third its size, orbiting each other every 7.45 days. With an orbital period of 49.5 days, eighteen transits of the inner planet have been observed, allowing a detailed characterization of its orbit and those of the stars. The outer planet's orbital period is 303.2 days, and although the planet is not Earth-like, it resides within the classical "habitable zone", where liquid water could exist on an Earth-like planet. With its two known planets, Kepler-47 establishes that close binary stars can host complete planetary systems.

We discuss the discovery and characterization of the circumbinary planet Kepler-38b. The stellar binary is single-lined, with a period of 18.8 days, and consists of a moderately evolved main-sequence star (M_A = 0.949 +/- 0.059 solar masses and R_A = 1.757 +/- 0.034 solar radii) paired with a low-mass star (M_B = 0.249 +/- 0.010 solar masses and R_B = 0.2724 +/- 0.0053 solar radii) in a mildly eccentric (e=0.103) orbit. A total of eight transits due to a circumbinary planet crossing the primary star were identified in the Kepler light curve (using Kepler Quarters 1 through 11), from which a planetary period of 105.595 +/- 0.053 days can be established. A photometric dynamical model fit to the radial velocity curve and Kepler light curve yields a planetary radius of 4.35 +/- 0.11 Earth radii, or 1.12 +/- 0.03 Neptune radii. Since the planet is not sufficiently massive to observably alter the orbit of the binary from Keplerian motion, we can only place an upper limit on the mass of the planet of 122 Earth masses (7.11 Neptune masses or 0.384 Jupiter masses) at 95% confidence. This upper limit should decrease as more Kepler data become available.

We present an analysis of the significantly expanded HARPS 2011 radial velocity data set for GJ 581 that was presented by Forveille et al. (2011). Our analysis reaches substantially different conclusions regarding the evidence for a Super-Earth-mass planet in the star's Habitable Zone. We were able to reproduce their reported \chi_{\nu}^2 and RMS values only after removing some outliers from their models and refitting the trimmed down RV set. A suite of 4000 N-body simulations of their Keplerian model all resulted in unstable systems and revealed that their reported 3.6\sigma detection of e=0.32 for the eccentricity of GJ 581e is manifestly incompatible with the system's dynamical stability. Furthermore, their Keplerian model, when integrated only over the time baseline of the observations, significantly increases the \chi_{\nu}^2 and demonstrates the need for including non-Keplerian orbital precession when modeling this system. We find that a four-planet model with all of the planets on circular or nearly circular orbits provides both an excellent self-consistent fit to their RV data and also results in a very stable configuration. The periodogram of the residuals to a 4-planet all-circular-orbit model reveals significant peaks that suggest one or more additional planets in this system. We conclude that the present 240-point HARPS data set, when analyzed in its entirety, and modeled with fully self-consistent stable orbits, by and of itself does offer significant support for a fifth signal in the data with a period near 32 days. This signal has a False Alarm Probability of <4% and is consistent with a planet of minimum mass of 2.2 Earth-masses, orbiting squarely in the star's Habitable Zone at 0.13 AU, where liquid water on planetary surfaces is a distinct possibility

Keck/HIRES precision radial velocities of HD 207832 indicate the presence of two Jovian-type planetary companions in Keplerian orbits around this G star. The planets have minimum masses of 0.56 and 0.73 Jupiter-masses with orbital periods of ~162 and ~1156 days, and eccentricities of 0.13 and 0.27, respectively. Stromgren b and y photometry reveals a clear stellar rotation signature of the host star with a period of 17.8 days, well separated from the period of the radial velocity variations, reinforcing their Keplerian origin. The values of the semimajor axes of the planets suggest that these objects have migrated from the region of giant planet formation to closer orbits. In order to examine the possibility of the existence of additional (small) planets in the system, we studied the orbital stability of hypothetical terrestrial-sized objects in the region between the two planets and interior to the orbit of the inner body. Results indicated that stable orbits exist only in a small region interior to planet b. However, the current observational data offer no evidence for the existence of additional objects in this system.

Debates regarding the age and inclination of the planetary system orbiting HR 8799, and the release of additional astrometric data following the discovery of the fourth planet prompted us to examine the possibility of constraining these two quantities by studying the long-term stability of this system at different orbital inclinations and in its high-mass configuration (7-10-10-10 MJup). We carried out ~1.5 million N-body integrations for different combinations of orbital elements of the four planets. The most dynamically stable combinations survived less than ~5 Myr at inclinations of 0{\deg} and 13{\deg}, and 41, 46, and 31 Myr at 18{\deg}, 23{\deg}, and 30{\deg}, respectively. Given such short lifetimes and the location of the system on the age-luminosity diagram for low-mass objects, the most reasonable conclusion of our study is that the planetary masses are less than 7-10-10-10 MJup and the system is quite young. Two trends to note from our work are as follows. (1) In the most stable systems, the higher the inclination, the more the coordinates for planets b and c diverge from the oldest archival astrometric data (released after we completed our N-body integrations), suggesting that either these planets are in eccentric orbits or have lower orbital inclinations than that of planet d. (2) The most stable systems place planet e closer to the central star than is observed, supporting the conclusion that the planets are more massive and the system is young. We present the details of our simulations and discuss the implications of the results.

In this work, we investigate the eclipse timing of the polar binary HU Aquarii that has been observed for almost two decades. Recently, Qian et al. attributed large (O-C) deviations between the eclipse ephemeris and observations to a compact system of two massive jovian companions. We improve the Keplerian, kinematic model of the Light Travel Time (LTT) effect and re-analyse the whole currently available data set. We add almost 60 new, yet unpublished, mostly precision light curves obtained using the time high-resolution photo-polarimeter OPTIMA, as well as photometric observations performed at the MONET/N, PIRATE and TCS telescopes. We determine new mid--egress times with a mean uncertainty at the level of 1 second or better. We claim that because the observations that currently exist in the literature are non-homogeneous with respect to spectral windows (ultraviolet, X-ray, visual, polarimetric mode) and the reported mid--egress measurements errors, they may introduce systematics that affect orbital fits. Indeed, we find that the published data, when taken literally, cannot be explained by any unique solution. Many qualitatively different and best-fit 2-planet configurations, including self-consistent, Newtonian N-body solutions may be able to explain the data. However, using high resolution, precision OPTIMA light curves, we find that the (O-C) deviations are best explained by the presence of a single circumbinary companion orbiting at a distance of ~4.5 AU with a small eccentricity and having ~7 Jupiter-masses. This object could be the next circumbinary planet detected from the ground, similar to the announced companions around close binaries HW Vir, NN Ser, UZ For, DP Leo or SZ Her, and planets of this type around Kepler-16, Kepler-34 and Kepler-35.

Prior to the detection of its outermost Uranus-mass object, it had been suggested that GJ 876 could host an Earth-sized planet in a 15-day orbit. Observation, however, did not support this idea, but instead revealed evidence for the existence of a larger body in a $\sim$125-day orbit, near a three-body resonance with the two giant planets of this system. In this paper, we present a detailed analysis of the dynamics of the four-planet system of GJ 876, and examine the possibility of the existence of other planetary objects interior to its outermost body. We have developed a numerical scheme that enables us to search the orbital parameter-space very effectively and, in a short time, identify regions where an object may be stable. We present details of this integration method and discuss its application to the GJ 876 four-planet system. The results of our initial analysis suggested possible stable orbits at regions exterior to the orbit of the outermost planet and also indicated that an island of stability may exist in and around the 15-day orbit. However, examining the long-term stability of an object in that region by direct integration revealed that the 15-day orbit becomes unstable and that the system of GJ 876 is most likely dynamically full. We present the results of our study and discuss their implications for the formation and final orbital architecture of this system.

In a recent study, Lee et al. presented new photometric follow-up timing observations of the semi-detached binary system SZ Herculis and proposed the existence of two hierarchical cirumbinary companions. Based on the light-travel time effect, the two low-mass M-dwarf companions are found to orbit the binary pair on moderate to high eccentric orbits. The derived periods of these two companions are close to a 2:1 mean-motion orbital resonance. We have studied the stability of the system using the osculating orbital elements as presented by Lee et al. Results indicate an orbit-crossing architecture exhibiting short-term dynamical instabilities leading to the escape of one of the proposed companions. We have examined the system's underlying model parameter-space by following a Monte Carlo approach and found an improved fit to the timing data. A study of the stability of our best-fitting orbits also indicates that the proposed system is generally unstable. If the observed anomalous timing variations of the binary period is due to additional circumbinary companions, then the resulting system should exhibit a long-term stable orbital configuration much different from the orbits suggested by Lee et al. We, therefore, suggest that based on Newtonian-dynamical considerations, the proposed quadruple system cannot exist. To uncover the true nature of the observed period variations of this system, we recommend future photometric follow-up observations that could further constrain eclipse-timing variations and/or refine light-travel time models.

We present a model-independent technique for calculating the time of mid-transits. This technique, named "barycenter method", uses the light-curve's symmetry to determine the transit timing by calculating the transit light-curve barycenter. Unlike the other methods of calculating mid-transit timing, this technique does not depend on the parameters of the system and central star. We demonstrate the capabilities of the barycenter method by applying this technique to some known transiting systems including several \emph{Kepler} confirmed planets. Results indicate that for complete and symmetric transit lightcurves, the barycenter method achieves the same precision as other techniques, but with fewer assumptions and much faster. Among the transiting systems studied with the barycenter method, we focus in particular on LHS 6343C, a brown dwarf that transits a member of an M+M binary system, LHS 6343AB. We present the results of our analysis, which can be used to set an upper limit on the period and mass of a possible second small perturber.

We present observations of comet-like main-belt object P/2010 R2 (La Sagra) obtained by Pan-STARRS 1 and the Faulkes Telescope-North on Haleakala in Hawaii, the University of Hawaii 2.2 m, Gemini-North, and Keck I telescopes on Mauna Kea, the Danish 1.54 m telescope at La Silla, and the Isaac Newton Telescope on La Palma. An antisolar dust tail is observed from August 2010 through February 2011, while a dust trail aligned with the object's orbit plane is also observed from December 2010 through August 2011. Assuming typical phase darkening behavior, P/La Sagra is seen to increase in brightness by >1 mag between August 2010 and December 2010, suggesting that dust production is ongoing over this period. These results strongly suggest that the observed activity is cometary in nature (i.e., driven by the sublimation of volatile material), and that P/La Sagra is therefore the most recent main-belt comet to be discovered. We find an approximate absolute magnitude for the nucleus of H_R=17.9+/-0.2 mag, corresponding to a nucleus radius of ~0.7 km, assuming an albedo of p=0.05. Using optical spectroscopy, we find no evidence of sublimation products (i.e., gas emission), finding an upper limit CN production rate of Q_CN<6x10^23 mol/s, from which we infer an H2O production rate of Q_H2O<10^26 mol/s. Numerical simulations indicate that P/La Sagra is dynamically stable for >100 Myr, suggesting that it is likely native to its current location and that its composition is likely representative of other objects in the same region of the main belt, though the relatively close proximity of the 13:6 mean-motion resonance with Jupiter and the (3,-2,-1) three-body mean-motion resonance with Jupiter and Saturn mean that dynamical instability on larger timescales cannot be ruled out.

In their model for the origin of the parent bodies of iron meteorites, Bottke et al proposed differentiated planetesimals that were formed in the region of 1-2 AU during the first 1.5 Myr, as the parent bodies, and suggested that these objects and their fragments were scattered into the asteroid belt as a result of interactions with planetary embryos. Although viable, this model does not include the effect of a giant planet that might have existed or been growing in the outer regions. We present the results of a concept study where we have examined the effect of a planetary body in the orbit of Jupiter on the early scattering of planetesimals from terrestrial region into the asteroid belt. We integrated the orbits of a large battery of planetesimals in a disk of planetary embryos, and studied their evolutions for different values of the mass of the planet. Results indicate that when the mass of the planet is smaller than 10 Earth-masses, its effects on the interactions among planetesimals and planetary embryos is negligible. However, when the planet mass is between 10 and 50 Earth-masses, simulations point to a transitional regime with ~50 Earth-mass being the value for which the perturbing effect of the planet can no longer be ignored. Simulations also show that further increase of the mass of the planet strongly reduces the efficiency of the scattering of planetesimals from the terrestrial planet region into the asteroid belt. We present the results of our simulations and discuss their possible implications for the time of giant planet formation.

Main belt asteroid (300163) 2006 VW139 (later designated P/2006 VW139) was discovered to exhibit comet-like activity by the Pan-STARRS1 survey telescope using automated point-spread-function analyses performed by PS1's Moving Object Processing System. Deep follow-up observations show both a short (\sim 10") antisolar dust tail and a longer (\sim 60") dust trail aligned with the object's orbit plane, similar to the morphology observed for another main-belt comet, P/2010 R2 (La Sagra), and other well-established comets, implying the action of a long-lived, sublimation-driven emission event. Photometry showing the brightness of the near-nucleus coma remaining constant over \sim 30 days provides further evidence for this object's cometary nature, suggesting it is in fact a main-belt comet, and not a disrupted asteroid. A spectroscopic search for CN emission was unsuccessful, though we find an upper limit CN production rate of Q_CN < 1.3x10^24 mol/s, from which we infer a water production rate of Q_H2O < 10^26 mol/s. We also find an approximately linear optical spectral slope of 7.2%/1000A, similar to other cometary dust comae. Numerical simulations indicate that P/2006 VW139 is dynamically stable for > 100 Myr, while a search for a potential asteroid family around the object reveals a cluster of 24 asteroids within a cutoff distance of 68 m/s. At 70 m/s, this cluster merges with the Themis family, suggesting that it could be similar to the Beagle family to which another main-belt comet, 133P/Elst-Pizarro, belongs.

We re-analyze 4 years of HARPS spectra of the nearby M1.5 dwarf GJ 667C available through the ESO public archive. The new radial velocity (RV) measurements were obtained using a new data analysis technique that derives the Doppler measurement and other instrumental effects using a least-squares approach. Combining these new 143 measurements with 41 additional RVs from the Magellan/PFS and Keck/HIRES spectrometers, reveals 3 additional signals beyond the previously reported 7.2-day candidate, with periods of 28 days, 75 days, and a secular trend consistent with the presence of a gas giant (Period sim 10 years). The 28-day signal implies a planet candidate with a minimum mass of 4.5 Mearth orbiting well within the canonical definition of the star's liquid water habitable zone, this is, the region around the star at which an Earth-like planet could sustain liquid water on its surface. Still, the ultimate water supporting capability of this candidate depends on properties that are unknown such as its albedo, atmospheric composition and interior dynamics. The 75-day signal is less certain, being significantly affected by aliasing interactions among a potential 91-day signal, and the likely rotation period of the star at 105 days detected in two activity indices. GJ 667C is the common proper motion companion to the GJ 667AB binary, which is metal poor compared to the Sun. The presence of a super-Earth in the habitable zone of a metal poor M dwarf in a triple star system, supports the evidence that such worlds should be ubiquitous in the Galaxy.

In this work we propose a new orbital architecture for the two proposed circumbinary planets around the polar eclipsing binary HU Aquarii. We base the new two-planet, light-travel time model on the result of a Monte Carlo simulation driving a least-squares Levenberg-Marquardt minimisation algorithm on the observed eclipse egress times. Our best-fitting model with $\chi_{r}^2=1.43$ resulted in high final eccentricities for the two companions leading to an unstable orbital configuration. From a large ensemble of initial guesses we examined the distribution of final eccentricities and semi-major axes for different $\chi_{r}^2$ parameter intervals and encountered qualitatively a second population of best-fitting parameters. The main characteristic of this population is described by low-eccentric orbits favouring long-term orbital stability of the system. We present our best-fitting model candidate for the proposed two-planet system and demonstrate orbital stability over one million years using numerical integrations.

We present observations and a dynamical analysis of the comet-like main-belt object, (596) Scheila. V-band photometry obtained on UT 2010 December 12 indicates that Scheila's dust cloud has a scattering cross-section ~1.4 times larger than that of the nucleus, corresponding to a dust mass of M_d~3x10^7 kg. V-R color measurements indicate that both the nucleus and dust are redder than the Sun, with no significant color differences between the dust cloud's northern and southern plumes. We also undertake an ultimately unsuccessful search for CN emission, where we find CN and H2O production rates of Q(CN) < 9x10^23 s^-1 and Q(H2O) < 10^27 s^{-1}. Numerical simulations indicate that Scheila is dynamically stable for >100 Myr, suggesting that it is likely native to its current location. We also find that it does not belong to a dynamical asteroid family of any significance. We consider sublimation-driven scenarios that could produce the appearance of multiple plumes of dust emission, but reject them as being physically implausible. Instead, we concur with previous studies that the unusual morphology of Scheila's dust cloud is most simply explained by a single oblique impact, meaning this object is likely not a main-belt comet, but is instead the second disrupted asteroid after P/2010 A2 (LINEAR) to be discovered.

Super-Earths, a class of planetary bodies with masses ranging from a few Earth-masses to slightly smaller than Uranus, have recently found a special place in the exoplanetary science. Being slightly larger than a typical terrestrial planet, super-Earths may have physical and dynamical characteristics similar to those of Earth whereas unlike terrestrial planets, they are relatively easier to detect. Because of their sizes, super-Earths can maintain moderate atmospheres and possibly dynamic interiors with plate tectonics. They also seem to be more common around low-mass stars where the habitable zone is in closer distances. This article presents a review of the current state of research on super-Earths, and discusses the models of the formation, dynamical evolution, and possible habitability of these objects. Given the recent advances in detection techniques, the detectability of super-Earths is also discussed, and a review of the prospects of their detection in the habitable zones of low-mass stars is presented.

We present the results of an extensive study of the detectability of Earth-sized planets and super-Earths in the habitable zones of cool and low-mass stars using transit timing variation method. We have considered a system consisting of a star, a transiting giant planet, and a terrestrial-class perturber, and calculated TTVs for different values of the parameters of the system. To identify ranges of the parameters for which these variations would be detectable by Kepler, we considered the analysis presented by Ford et al. (2011, ArXiv:1102.0544) and assumed that a peak-to-peak variation of 20 seconds would be within the range of the photometric sensitivity of this telescope. We carried out simulations for resonant and non-resonant orbits, and identified ranges of the semimajor axes and eccentricities of the transiting and perturbing bodies for which an Earth-sized planet or a super-Earth in the habitable zone of a low-mass star would produce such TTVs. Results of our simulations indicate that in general, outer perturbers near first- and second-order resonances show a higher prospect for detection. Inner perturbers are potentially detectable only when near 1:2 and 1:3 mean-motion resonances. For a typical M star with a Jupiter-mass transiting planet, for instance, an Earth-mass perturber in the habitable zone can produce detectable TTVs when the orbit of the transiting planet is between 15 and 80 days. We present the details of our simulations and discuss the implication of the results for the detection of terrestrial planets around different low-mass stars.

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 has been shown that an Earth-size planet or a super-Earth, in resonance with a transiting Jupiter-like body in a short-period orbit around an M star, can create detectable TTV signals (Kirste \& Haghighipour, 2011). Given the low masses of M stars and their circumstellar disks, it is expected that such a transiting giant planet to have formed at large distances and migrated to its close-in orbit. That implies, if such systems are discovered around M stars, the terrestrial planet had to form during the migration of the giant planet. The formation of this object may be either in-situ (in a close-in orbit) followed by its capture in resonance, or the object is formed at larger distances where it was subsequently captured in a resonance with the migrating giant planet. We have investigated these two scenarios by simulating the dynamics of a disk of protoplanetary embryos and the formation of terrestrial planets during the migration of a Jupiter-like planet around an M star. Results suggest that unless the migration of the giant planet is very slow (slower than 1E-7 AU/year), it is unlikely that the close-in terrestrial planet is formed in-situ. If a terrestrial planet is detected in a mean-motion resonance with a close-in giant planet around an M star, the terrestrial planet was most likley formed at large distances and carried to its close-in resonant orbit by the migrating giant body.

M dwarfs constitute more than 70% of the stars in the solar neighborhood. They are cooler and smaller than Sun-like stars and have less-massive disks which suggests that planets around these stars are more likely to be Neptune-size or smaller. The transit depths and transit times of planets around M stars are large and well-matched to the Kepler temporal resolution. As a result, M stars have been of particular interest for searching for planets in both radial velocity and transit photometry surveys. We have recently started a project on searching for possible planet-hosting M stars in the publicly available data from Kepler space telescope. We have used four criteria, namely, the magnitude, proper motion, H-Ks and J-H colors, and searched for M stars in Q0 and Q1 data sets. We have been able to find 108 M stars among which 54 had not been previously identified among Kepler's targets. We discuss the details of our selection process and present the results.

We present new precise HIRES radial velocity (RV) data sets of five nearby stars obtained at Keck Observatory. HD 31253, HD 218566, HD 177830, HD 99492 and HD 74156 are host stars of spectral classes F through K and show radial velocity variations consistent with new or additional planetary companions in Keplerian motion. The orbital parameters of the candidate planets in the five planetary systems span minimum masses of M sin i = 27.43 M_{earth} to M sin i = 8.28 M_{jup}, periods of 17.05 to 4696.95 days and eccentricities ranging from circular to extremely eccentric (e ~ 0.63). The 5th star, HD 74156, was known to have both a 52-day and a 2500-day planet, and was claimed to also harbor a 3rd planet at 336d, in apparent support of the "Packed Planetary System" hypothesis. Our greatly expanded data set for HD 74156 provides strong confirmation of both the 52-day and 2500-d planets, but strongly contradicts the existence of a 336-day planet, and offers no significant evidence for any other planets in the system.

We present 11 years of HIRES precision radial velocities (RV) of the nearby M3V star Gliese 581, combining our data set of 122 precision RVs with an existing published 4.3-year set of 119 HARPS precision RVs. The velocity set now indicates 6 companions in Keplerian motion around this star. Differential photometry indicates a likely stellar rotation period of ~94 days and reveals no significant periodic variability at any of the Keplerian periods, supporting planetary orbital motion as the cause of all the radial velocity variations. The combined data set strongly confirms the 5.37-day, 12.9-day, 3.15-day, and 67-day planets previously announced by Bonfils et al. (2005), Udry et al. (2007), and Mayor et al (2009). The observations also indicate a 5th planet in the system, GJ 581f, a minimum-mass 7.0 M_Earth planet orbiting in a 0.758 AU orbit of period 433 days and a 6th planet, GJ 581g, a minimum-mass 3.1 M_Earth planet orbiting at 0.146 AU with a period of 36.6 days. The estimated equilibrium temperature of GJ 581g is 228 K, placing it squarely in the middle of the habitable zone of the star and offering a very compelling case for a potentially habitable planet around a very nearby star. That a system harboring a potentially habitable planet has been found this nearby, and this soon in the relatively early history of precision RV surveys, indicates that eta_Earth, the fraction of stars with potentially habitable planets, is likely to be substantial. This detection, coupled with statistics of the incompleteness of present-day precision RV surveys for volume-limited samples of stars in the immediate solar neighborhood suggests that eta_Earth could well be on the order of a few tens of percent. If the local stellar neighborhood is a representative sample of the galaxy as a whole, our Milky Way could be teeming with potentially habitable planets.

(Abreviated) Continued radial velocity monitoring of the nearby M4V red dwarf star GJ~876 with Keck/HIRES has revealed the presence of a Uranus-mass fourth planetary companion in the system. The new planet has a mean period of $P_e=126.6$ days (over the 12.6-year baseline of the radial velocity observations), and a minimum mass of $m_e\sin{i_e}=12.9\pm 1.7\,M_{\oplus}$. Self-consistent, N-body fits to the radial velocity data set show that the four-planet system has an invariable plane with an inclination relative to the plane of the sky of $i=59.5^{\circ}$. The fit is not significantly improved by the introduction of a mutual inclination between the planets ``b'' and ``c,'' but the new data do confirm a non-zero eccentricity, $e_d=0.207\pm0.055$ for the innermost planet, ``d.'' In our best-fit coplanar model, the mass of the new component is $m_e=14.6\pm1.7\,M_{\oplus}$. Our best-fitting model places the new planet in a 3-body resonance with the previously known giant planets (which have mean periods of $P_c=30.4$ and $P_b=61.1$ days). The critical argument, $\varphi_{\rm Laplace}=\lambda_c-3\lambda_b+2\lambda_e$, for the Laplace resonance librates with an amplitude of $\Delta\varphi_{\rm Laplace}=40\pm13^{\circ}$ about $\varphi_{\rm Laplace}=0^{\circ}$. Numerical integration indicates that the four-planet system is stable for at least a billion years (at least for the coplanar cases). This resonant configuration of three giant planets orbiting an M-dwarf primary differs from the well-known Laplace configuration of the three inner Galilean satellites of Jupiter, which are executing very small librations about $\varphi_{\rm Laplace}=180^{\circ}$, and which never experience triple conjunctions. The GJ~876 system, by contrast, comes close to a triple conjunction between the outer three planets once per every orbit of the outer planet, ``e.''

We analyze the dynamics of individual kilometer-size planetesimals in circumstellar orbits of a tight binary system. We include both the gravitational perturbations of the secondary star and a non-linear gas drag stemming from an eccentric gas disk with a finite precession rate. We consider several precession rates and eccentricities for the gas, and compare the results with a static disk in circular orbit. The disk precession introduces three main differences with respect to the classical static case: (i) The equilibrium secular solutions generated by the gas drag are no longer fixed points in the averaged system, but limit cycles with frequency equal to the precession rate of the gas. The amplitude of the cycle is inversely dependent on the body size, reaching negligible values for $\sim 50$ km size planetesimals. (ii) The maximum final eccentricity attainable by small bodies is restricted to the interval between the gas eccentricity and the forced eccentricity, and apsidal alignment is no longer guaranteed for planetesimals strongly coupled with the gas. (iii) The characteristic timescales of orbital decay and secular evolution decrease significantly with increasing precession rates, with values up to two orders of magnitude smaller than for static disks. Finally, we apply this analysis to the $\gamma$-Cephei system and estimate impact velocities for different size bodies and values of the gas eccentricity. For high disk eccentricities, we find that the disk precession decreases the velocity dispersion between different size planetesimals, thus contributing to accretional collisions in the outer parts of the disk. The opposite occurs for almost circular gas disks, where precession generates an increase in the relative velocities.

Precision radial velocities from Keck/HIRES reveal a Saturn-mass planet orbiting the nearby M4V star HIP 57050. The planet has a minimum mass of 0.3 Jupiter-mass, an orbital period of 41.4 days, and an orbital eccentricity of 0.31. V-band photometry reveals a clear stellar rotation signature of the host star with a period of 98 days, well separated from the period of the radial velocity variations and reinforcing a Keplerian origin for the observed velocity variations. The orbital period of this planet corresponds to an orbit in the habitable zone of HIP 57050, with an expected planetary temperature of approximately 230 K. The star has a metallicity of [Fe/H] = 0.32+/-0.06 dex, of order twice solar and among the highest metallicity stars in the immediate solar neighborhood. This newly discovered planet provides further support that the well-known planet-metallicity correlation for F, G, and K stars also extends down into the M-dwarf regime. The a priori geometric probability for transits of this planet is only about 1%. However, the expected eclipse depth is ~7%, considerably larger than that yet observed for any transiting planet. Though long on the odds, such a transit is worth pursuing as it would allow for high quality studies of the atmosphere via transmission spectroscopy with HST. At the expected planetary effective temperature, the atmosphere may contain water clouds.

The discovery of Main Belt Comets (MBCs) has raised many questions regarding the origin and activation mechanism of these objects. Results of a study of the dynamics of these bodies suggest that MBCs were formed in-situ as the remnants of the break-up of large icy asteroids. Simulations show that similar to the asteroids in the main belt, MBCs with orbital eccentricities smaller than 0.2 and inclinations lower than 25 degrees have stable orbits implying that many MBCs with initially larger eccentricities and inclinations might have been scattered to other regions of the asteroid belt. Among scattered MBCs, approximately 20 percent reach the region of terrestrial planets where they might have contributed to the accumulation of water on Earth. Simulations also show that collisions among MBCs and small objects could have played an important role in triggering the cometary activity of these bodies. Such collisions might have exposed sub-surface water ice which sublimated and created thin atmospheres and tails around MBCs. This paper discusses the results of numerical studies of the dynamics of MBCs and their implications for the origin of these objects. The results of a large numerical modeling of the collisions of m-sized bodies with km-sized asteroids in the outer part of the asteroid belt are also presented and the viability of the collision-triggering activation scenario is discussed.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

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.

In an effort to understand the origin of the Main Belt Comets (MBCs) 7968 Elst-Pizzaro, 118401, and P/2005 U1, the dynamics of these three icy asteroids and a large number of hypothetical MBCs were studied. Results of extensive numerical integrations of these objects suggest that these MBCs were formed in-place through the collisional break up of a larger precursor body. Simulations point specifically to the Themis family of asteroids as the origin of these objects and rule out the possibility of a cometary origin (i.e. inward scattering of comets from outer solar system and their primordial capture in the asteroid belt). Results also indicate that while 7968 Elst-Pizzaro and 118401 maintain their orbits for 1 Gyr, P/2005 U1 diffuses chaotically in eccentricity and becomes unstable in ~20 Myr. The latter suggest that this MBC used to be a member of the Themis family and is now escaping away. Numerical integrations of the orbits of hypothetical MBCs in the vicinity of the Themis family show a clustering of stable orbits (with eccentricities smaller than 0.2 and inclinations less than 25 deg.) suggesting that many more MBCs may exist in the vicinity of this family (although they might have not been activated yet). The details of the results of simulations and the constraints on the models of the formation and origins of MBCs are presented, and their implications for the detection of more of these objects are discussed.

Given the tendency of planets to form in multiples, and the observational evidence in support of the existence of potential planet-hosting stars in binaries or clusters, it is expected that extrasolar terrestrial planes are more likely to be found in multiple body systems. This paper discusses the prospects of the detection of terrestrial/habitable planets in multibody systems by presenting the results of a study of the long-term stability of these objects in systems with multiple giant planets (particularly those in eccentric and/or in mean-motion resonant orbits), systems with close-in Jupiter-like bodies, and systems of binary stars. The results of simulations show that while short-period terrestrial-class objects that are captured in near mean-motion resonances with migrating giant planets are potentially detectable via transit photometry or the measurement of the variations of the transit-timing due to their close-in Jovian-mass planetary companions, the prospect of the detection of habitable planets with radial velocity technique is higher in systems with multiple giant planets outside the habitable zone and binary systems with moderately separated stellar companions.

A survey of currently known planet-hosting stars indicates that approximately 25% of extrasolar planetary systems are within dual-star environments. Several of these systems contain stellar companions on moderately close orbits, implying that studies of the formation and dynamical evolution of giant and terrestrial planets, in and around binary star systems have now found realistic grounds. With the recent launch of the space telescope CoRoT, and the launch of NASA's Kepler satellite in 2009, the number of such dynamically complex systems will soon increase and many more of their diverse and interesting dynamical characteristics will soon be discovered. It is therefore, both timely and necessary, to obtain a deep understanding of the history and current status of research on planets in binary star systems. This chapter will serve this purpose by reviewing the models of the formation of giant and terrestrial planets in dual-star environments, and by presenting results of the studies of their dynamical evolution and habitability, as well as the mechanisms of delivery of water and other volatiles to their terrestrial-class objects. In this chapter, the reader is presented with a comprehensive, yet relatively less technical approach to the study of planets in and around binary stars, and with discussions on the differences between dynamical characteristics of these systems and planetary systems around single stars.

We present a study of the newly-discovered main-belt comet P/2008 R1 (Garradd), an object with the dynamical characteristics of an asteroid and the physical characteristics of a comet. Photometry sets a limit to the effective radius of the nucleus at r_e < 0.7 km (red geometric albedo 0.05 assumed). The coma shows a secular fading in our data caused by the escape of dust particles from the near-nucleus environment. The optical reflection spectrum is a nearly neutral continuum devoid of gaseous emission lines, from which we derive a limit to the cyanide (CN) radical production rate of Q_CN <1.4e23/s and infer a mass loss rate <1.5 kg/s at the time of our observations. Unlike the first-reported main-belt comets, P/2008 R1 is not dynamically stable. The nearby 8:3 mean-motion resonance with Jupiter induces dynamical instability on timescales 20 to 30 Myr. Hence, we conclude that P/2008 R1 has recently arrived from a more stable source elsewhere. The high Tisserand parameter of the orbit (in fact, with T_J = 3.216 it is the highest of any comet) points to a source in the asteroid belt itself, instead of in the Kuiper belt (putative source of the Jupiter family comets). We infer that P/2008 R1 is an icy body from the outer asteroid belt in which sublimation has been triggered by rising temperatures resulting from a decreasing perihelion distance.

We present visible and near-infrared reflectance spectra and interpreted surface mineralogy for asteroid 10537 (1991 RY16). The spectrum of this object is without precedent amongst the Main Belt asteroids. A unique absorption band centered at 0.63 microns could be attributed to one of several mineralogies. Pronounced 1- and 2-micron absorption bands suggest that the composition of 10537 is a mixture of pyroxenes and olivine and that it originated from a parent body that was partially or fully differentiated. The closest available analog is the large Main Belt asteroid 349 Dembowska but 10537 may be an isolated fragment from a completely eroded parent body.

An interesting feature of the giant planets of our solar system is the existence of regions around these objects where no irregular satellites are observed. Surveys have shown that, around Jupiter, such a region extends from the outermost regular satellite Callisto, to the vicinity of Themisto, the innermost irregular satellite. To understand the reason for the existence of such a satellite-void region, we have studied the dynamical evolution of Jovian irregulars by numerically integrating the orbits of several hundred test particles, distributed in a region between 30 and 80 Jupiter-radii, for different values of their semimajor axes, orbital eccentricities, and inclinations. As expected, our simulations indicate that objects in or close to the influence zones of the Galilean satellites become unstable because of interactions with Ganymede and Callisto. However, these perturbations cannot account for the lack of irregular satellites in the entire region between Callisto and Themisto. It is suggested that at distances between 60 and 80 Jupiter-radii, Ganymede and Callisto may have long-term perturbative effects, which may require the integrations to be extended to times much longer than 10 Myr. The interactions of irregular satellites with protosatellites of Jupiter at the time of the formation of Jovian regulars may also be a destabilizing mechanism in this region. We present the results of our numerical simulations and discuss their applicability to similar satellite void-regions around other giant planets.

Understanding the formation and dynamical evolution of habitable planets in extrasolar planetary systems is a challenging task. In this respect, systems with multiple giant planets and/or multiple stars present special complications. The formation of habitable planets in these environments is strongly affected by the dynamics of their giant planets and/or their stellar companions. These objects have profound effects on the structure of the disk of planetesimals and protoplanetary objects in which terrestrial-class planets are formed. To what extent the current theories of planet formation can be applied to such "extreme" planetary systems depends on the dynamical characteristics of their planets and/or their binary stars. In this paper, I present the results of a study of the possibility of the existence of Earth-like objects in systems with multiple giant planets (namely Upsilon Andromedae, 47 UMa, GJ 876, and 55 Cnc) and discuss the dynamics of the newly discovered Neptune-size object in 55 Cnc system. I will also review habitable planet formation in binary systems and present the results of a systematic search of the parameter-space for which Earth-like objects can form and maintain long-term stable orbits in the habitable zones of binary stars.

With more than 260 extrasolar planetary systems discovered to-date, the search for habitable planets has found new grounds. Unlike our solar system, the stars of many of these planets are hosts to eccentric or close-in giant bodies. Several of these stars are also members of moderately close ($<$40 AU) binary or multi-star systems. The formation of terrestrial objects in these "extreme" environments is strongly affected by the dynamics of their giant planets and/or their stellar companions. These objects have profound effects on the chemical structure of the disk of planetesimals and the radial mixing of these bodies in the terrestrial regions of their host stars. For many years, it was believed that such effects would be so destructive that binary stars and also systems with close-in giant planets would not be able to form and harbor habitable bodies. Recent simulations have, however, proven otherwise. I will review the results of the simulations of the formation and long-term stability of Earth-like objects in the habitable zones of such "extreme" planetary systems, and discuss the possibility of the formation of terrestrial planets, with significant amounts of water, in systems with hot Jupiters, and also around the primaries of moderately eccentric close binary stars.

The search for habitable planets like Earth around other stars fulfils an ancient imperative to understand our origins and place in the cosmos. The past decade has seen the discovery of hundreds of planets, but nearly all are gas giants like Jupiter and Saturn. Recent advances in instrumentation and new missions are extending searches to planets the size of the Earth, but closer to their host stars. There are several possible ways such planets could form, and future observations will soon test those theories. Many of these planets we discover may be quite unlike Earth in their surface temperature and composition, but their study will nonetheless inform us about the process of planet formation and the frequency of Earth-like planets around other stars.

A survey of currently known extrasolar planets indicates that close to 20% of their hosting stars are members of binary systems. While the majority of these binaries are wide (i.e., with separations between 250 and 6500 AU), the detection of Jovian-type planets in the three binaries of Gamma Cephei (separation of 18.5 AU), GL 86 (separation of 21 AU), and HD 41004 (separation of 23 AU) have brought to the forefront questions on the formation of giant planets and the possibility of the existence of smaller bodies in moderately close binary star systems. This paper discusses the late stage of the formation of habitable planets in binary systems that host Jovian-type bodies, and reviews the effects of the binary companion on the formation of Earth-like planets in the system's habitable zone. The results of a large survey of the parameter-space of binary-planetary systems in search of regions where habitable planets can form and have long-term stable orbits are also presented.

One of the most surprising discoveries of extrasolar planets is the detection of planets in moderately close binary star systems. The Jovian-type planets in the two binaries of Gamma Cephei and GJ 86 have brought to the forefront questions on the formation of giant planets and the possibility of the existence of smaller bodies in such dynamically complex environments. The diverse dynamical characteristics of these objects have made scientists wonder to what extent the current theories of planet formation can be applied to binaries and multiple star systems. At present, the sensitivity of the detection techniques does not allow routine discovery of Earth-sized bodies in binary systems. However, with the advancement of new techniques, and with the recent launch of CoRoT and the launch of Kepler in late 2008, the detection of more planets (possibly terrestrial-class objects) in such systems is on the horizon. Theoretical studies and numerical modeling of terrestrial and habitable planet formation are, therefore, necessary to gain fundamental insights into the prospects for life in such systems and have great strategic impact on NASA science and missions.

All four giant planets in the Solar system possess irregular satellites, characterized by large, highly eccentric and/or inclined orbits that are distinct from the nearly circular, uninclined orbits of the regular satellites. This difference can be traced directly to different modes of formation. Whereas the regular satellites grew by accretion within circumplanetary disks the irregular satellites were captured from initially heliocentric orbits at an early epoch. Recently, powerful survey observations have greatly increased the number of known irregular satellites, permitting a fresh look at the group properties of these objects and motivating a re-examination of the mechanisms of capture. None of the suggested mechanisms, including gas-drag, pull-down, and three-body capture, convincingly fit the group characteristics of the irregular satellites. The sources of the satellites also remain unidentified.

Recent radial velocity observations have indicated that Jovian-type planets can exist in moderately close binary star systems. Numerical simulations of the dynamical stability of terrestrial-class planets in such environments have shown that, in addition to their giant planets, these systems can also harbor Earth-like objects. In this paper, we study the late stage of terrestrial planet formation in such binary-planetary systems, and present the results of the simulations of the formation of Earth-like bodies in their habitable zones. We consider a circumprimary disk of Moon- to Mars-sized objects and numerically integrate the orbits of these bodies at the presence of the Jovian-type planet of the system and for different values of the mass, semimajor axis, and orbital eccentricity of the secondary star. Results indicate that, Earth-like objects, with substantial amounts of water, can form in the habitable zone of the primary star. Simulations also indicate that, by transferring angular momentum from the secondary star to protoplanetary objects, the giant planet of the system plays a key role in the radial mixing of these bodies and the water contents of the final terrestrial planets. We will discuss the results of our simulation and show that the formation of habitable planets in binary-planetary systems is more probable in binaries with moderate to large perihelia.

It has been suggested that the long-lived residual radial velocity variations observed in the precision radial velocity measurements of the primary of Gamma Cephei (HR8974, HD222404, HIP116727) are likely due to a Jupiter-like planet around this star (Hatzes et al, 2003). In this paper, the orbital dynamics of this plant is studied and also the possibility of the existence of a hypothetical Earth-like planet in the habitable zone of its central star is discussed. Simulations, which have been carried out for different values of the eccentricity and semimajor axis of the binary, as well as the orbital inclination of its Jupiter-like planet, expand on previous studies of this system and indicate that, for the values of the binary eccentricity smaller than 0.5, and for all values of the orbital inclination of the Jupiter-like planet ranging from 0 to 40 degrees, the orbit of this planet is stable. For larger values of the binary eccentricity, the system becomes gradually unstable. Integrations also indicate that, within this range of orbital parameters, a hypothetical Earth-like planet can have a long-term stable orbit only at distances of 0.3 to 0.8 AU from the primary star. The habitable zone of the primary, at a range of approximately 3.1 to 3.8 AU, is, however, unstable.

A summary of the results of a numerical study of the growth of solid particles in the vicinity of an azimuthally symmetric density enhancement of a protostellar disk are presented. The effects of gas drag and pressure gradients on the rate of growth of dust particles and their settling on the midplane of the nebula are also discussed.

We present results of extensive numerical studies of the stability of non-interacting particles in the planetary systems of stars Upsilon Andromedae, GJ 876, 47 UMa, and 55 Cancri. We discuss the possibility of the existence of islands of stability and/or instability at different regions in these multi-body systems, and their probable correspondence to certain mean-motion resonances. The results of our study may be applied to questions concerning the stability of terrestrial planets in these systems and also the trapping of particles in resonances with the planets.

Asymptotic properties of electromagnetic waves are studied within the context of Friedmann-Robertson-Walker (FRW) cosmology. Electromagnetic fields are considered as small perturbations on the background spacetime and Maxwell's equations are solved for all three cases of flat, closed and open FRW universes. The asymptotic character of these solutions are investigated and their relevance to the problem of cosmological tails of electromagnetic waves is discussed.

Precision radial velocity measurements of the Gamma Cephei (HR8974) binary system suggest the existence of a planetary companion with a minimum mass of 1.7 Jupiter-mass on an elliptical orbit with a ~2.14 AU semimajor axis and 0.12 eccentricity (hatzes et al. 2003). I present in this paper a summary of the results of an extensive numerical study of the orbital stability of this three-body system for different values of the semimajor axis and orbital eccentricity of the binary, and also the orbital inclination of the planet. Numerical integrations indicate that the system is stable for the planet's orbital inclination ranging from 0 to 60 degrees, and for the binary's orbital eccentricity less than 0.5. The results also indicate that for large values of the inclination, the system may be locked in a Kozai resonance.

It has recently been stated that the warp structure observed around the star Beta Pictoris may be due to four planets embedded in its debris disk (Wahhaj et al., 2003). It, therefore, becomes important to investigate for what range of parameters, and for how long such a multibody system will be dynamically stable. We present the results of the numerical integration of the suggested planetary system for different values of the mass and radii of the planets, and their relative positions and velocities. We also present a statistical analysis of the results in an effort to understand the relation between the different regions of the parameter-space of the system and the duration of the orbital stability of the embedded planets.

The results of a numerical study of the growth of solid particles, ranging from 1 micron to 1 millimeter in size, in the vicinity of an azimuthally symmetric density enhancement of a protostellar disk are presented. It is shown that the combined effect of gas drag and pressure gradients, which causes solid objects to rapidly migrate toward the location of the local maximum density, can also enhance the rate of growth of dust particles to larger objects. The results of numerical simulations of such growth processes are presented and the effects of the changes of the physical properties of particles on the rates of their growth are also discussed.

We study the motions of small solids, ranging from micron-sized dust grains to meter-sized objects, in the vicinity of local pressure enhancements of a gaseous nebula. Integrating numerically, we show that as a result of the combined effect of gas drag and pressure gradients, solids tend to accumulate at the locations where the pressure of the gas maximizes. The rate of migration of solids varies with their sizes and densities and also with the physical properties of the gas. The results of our numerical simulations indicate that such migrations are most rapid for meter-sized objects. The applicability of the results to the enhancement of the collision and coagulation of solids and also to the growth-rate of planetesimals is discussed.

The results of an extensive numerical study of the periodic orbits of planar, elliptic restricted three-body planetary systems consisting of a star, an inner massive planet and an outer mass-less body in the external 1:2 mean-motion resonance are presented. Using the method of differential continuation, the locations of the resonant periodic orbits of such systems are identified and through an extensive study of their phase-parameter space, it is found that the majority of the resonant periodic orbits are unstable. For certain values of the mass and the orbital eccentricity of the inner planet, however, stable periodic orbits can be found. The applicability of such studies to the 1:2 resonance of the extrasolar planetary system GJ876 is also discussed.

In this paper, we present preliminary results on the stability of massless particles in two and three-planet systems. The results of our study may be used to address questions concerning the stability of terrestrial planets in these systems and also the trapping of particles in resonances with the planets. The possibility of the existence of islands of stability and/or instability at different regions in multi-body systems and their probable correspondence to certain resonances are also discussed.

The results of an extensive numerical study of the orbital dynamics of small bodies ranging from micron-sized dust grains to 1 km objects subject to gas drag and also the gravitational attraction of a non-uniform gaseous nebula are presented. The results indicate that it is possible for small bodies to migrate rapidly toward the locations of the maxima of the gas density where the probabilities of collisions and coagulations are enhanced.

We study the motions of small solids, ranging from micron-sized dust grains to 100-m objects, in the vicinity of a local density enhancement of an isothermal gaseous solar nebula. Being interested in possible application of the results to the formation of clumps and spiral arms in a circumstellar disk, we numerically integrate the equations of motion of such solids and study their migration for different values of their sizes and masses and also for different physical properties of the gas, such as its density and temperature. We show that, considering the drag force of the gas and also the gravitational attraction of the nebula, it is possible for solids, within a certain range of size and mass, to migrate rapidly (i.e. within ~1000 years) toward the location of a local maximum density where collisions and coagulation may result in an accelerated rate of planetesimal formation.

The existence of Uranus and Neptune presents severe difficulties for the core accretion model for the formation of ice giant planets. We suggest an alternative mechanism, namely disk instability leading to the formation of gas giant protoplanets, coagulation and settling of dust grains to form ice/rock cores at their centers, and photoevaporation of their gaseous envelopes by a nearby OB star, as a possible means of forming ice giant planets.

Based on the value of the orbital eccentricity of a particle and also its proximity to the exact resonant orbit in a three-body system, the Pendulum Approximation (Dermott & Murray 1983) or the Second Fundamental Model of Resonance (Andoyer 1903; Henrard & Lemaitre 1983) are commonly used to study the motion of that particle near its resonance state. In this paper, we present the method of partial averaging as an analytical approach to study the dynamical evolution of a body near a resonance. To focus attention on the capabilities of this technique, a restricted, circular and planar three-body system is considered and the dynamics of its outer planet while captured in a resonance with the inner body is studied. It is shown that the first-order partially averaged system resembles a mathematical pendulum whose librational motion can be viewed as a geometrical interpretation of the resonance capture phenomenon. The driving force of this pendulum corresponds to the gravitational attraction of the inner body and its contribution, at different resonant states, is shown to be proportional to e^s, where s is the order of the resonance and e is the orbital eccentricity of the outer planet. As examples of such systems, the cases of (1:1), (1:2) and (1:3) resonances are discussed and the results are compared with known planetary systems such as the Sun-Jupiter-Trojan asteroids.

Following the general numerical analysis of Melita and Woolfson (1996), I showed in a recent paper that a restricted, planar, circular planetary system consisting of Sun, Jupiter and Saturn would be captured in a near (2:1) resonance when one would allow for frictional dissipation due to interplanetary medium (Haghighipour, 1998). In order to analytically explain this resonance phenomenon, the method of partial averaging near a resonance was utilized and the dynamics of the first-order partially averaged system at resonance was studied. Although in this manner, the finding that resonance lock occurs for all initial relative positions of Jupiter and Saturn was confirmed, the first-order partially averaged system at resonance did not provide a complete picture of the evolutionary dynamics of the system and the similarity between the dynamical behavior of the averaged system and the main planetary system held only for short time intervals. To overcome these limitations, the method of partial averaging near a resonance is extended to the second order of perturbation in this paper and a complete picture of dynamical behavior of the system at resonance is presented. I show in this study that the dynamics of the second-order partially averaged system at resonance resembles the dynamical evolution of the main system during the resonance lock in general, and I present analytical explanations for the evolution of the orbital elements of the main system while captured in resonance.

A restricted planar circular three-body system, consisting of the Sun and two planets, is studied as a simple model for a planetary system. The mass of the inner planet is considered to be larger and the system is assumed to be moving in a uniform interplanetary medium with constant density. Numerical integrations of this system indicate a resonance capture when the dynamical friction of the interplanetary medium is taken into account. As a result of this resonance trapping, the ratio of orbital periods of the two planets becomes nearly commensurate and the eccentricity and semimajor axis of the orbit of the outer planet and also its angular momentum and total energy become constant. It appears from the numerical work that the resulting commensurability and also the resonant values of the orbital elements of the outer planet are essentially independent of the initial relative positions of the two bodies. The results of numerical integrations of this system are presented and the first-order partially averaged equations are studied in order to elucidate the behavior of the system while captured in resonance.