Source author record

Nader Haghighipour

Nader Haghighipour appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2022arXiv

Building Terrestrial Planets: Why results of perfect-merging simulations are not quantitatively reliable approximations to accurate modeling of terrestrial planet formation

Although it is accepted that perfect-merging is not a realistic outcome of collisions, some researchers state that perfect-merging simulations can still be considered as quantitatively reliable representations of the final stage of terrestrial planet formation. Citing the work of Kokubo & Genda [ApJL, 714L, 21], they argue that the differences between the final planets in simulations with perfect-merging and those where collisions are resolved accurately are small, and it is, therefore, justified to use perfect-merging results as an acceptable approximation to realistic simulations. In this paper, we show that this argument does not stand. We demonstrate that when the mass lost during collisions is taken into account, the final masses of the planets will be so different from those obtained from perfect-merging that the latter cannot be used as a valid approximation. We carried out a large number of SPH simulations of embryo-embryo collisions and determined the amount of the mass and water lost in each impact. We applied the results to collisions in a typical perfect-merging simulation and showed that even when the mass-loss in each collision is as small as 10%, perfect-merging can, on average, overestimate the masses of the final planets by $\sim 35\%$ and their water-contents by more than 18%. Our analysis demonstrates that, while perfect-merging simulations are still a powerful tool in proving concepts, they cannot be used to make predictions, draw quantitative conclusions (especially about the past history of a planetary system) and serve as a valid approximation to, or in lieu of the simulations in which collisions are resolved accurately.

preprint2020arXiv

Detailed Calculations of the Efficiency of Planetesimal Accretion in the Core-Accretion Model

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.

preprint2020arXiv

TOI-1338: TESS' First Transiting Circumbinary Planet

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.

preprint2019arXiv

Orbital Stability of Circumstellar Planets in Binary Systems

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.

preprint2016arXiv

A Catalog of Kepler Habitable Zone Exoplanet Candidates

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.

preprint2016arXiv

Kepler-1647b: the largest and longest-period Kepler transiting circumbinary planet

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.

preprint2016arXiv

Potential Jupiter-Family Comet Contamination of the Main Asteroid Belt

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.

preprint2016arXiv

Triggering Sublimation-Driven Activity of Main Belt Comets

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.

preprint2015arXiv

Evolution of circumbinary planets around eccentric binaries: The case of Kepler-34

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.

preprint2015arXiv

Formation of Terrestrial Planets in Disks with Different Surface Density Profiles

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.

preprint2015arXiv

Sublimation-Driven Activity in Main-Belt Comet 313P/Gibbs

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.

preprint2014arXiv

Calculating the Habitable Zone of Multiple Star Systems (http://astro.twam.info/hz)

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.

preprint2014arXiv

KIC 9632895 - The 10th Kepler Transiting Circumbinary Planet

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.

preprint2014arXiv

Modelling Circumbinary Planets: The case of Kepler-38

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.

preprint2014arXiv

Predicting a third planet in the Kepler-47 circumbinary system

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.

preprint2014arXiv

Properties and evolution of NEO families created by tidal disruption at Earth

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τ_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τ_{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.

preprint2013arXiv

Detection of Earth-mass and Super-Earth Trojan Planets Using Transit Timing Variation Method

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.

preprint2013arXiv

Kepler's Unparalleled Exploration of the Time Dimension

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.

preprint2013arXiv

Main-Belt Comet P/2012 T1 (PANSTARRS)

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.

preprint2013arXiv

The Formation and Dynamics of Super-Earth Planets

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.

preprint2012arXiv

A planetary system around the nearby M dwarf GJ 667C with at least one super-Earth in its habitable zone

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.

preprint2012arXiv

Can GJ 876 host four planets in resonance?

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.

preprint2012arXiv

Circumstellar Habitable Zones of Binary Star Systems in the Solar Neighborhood

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.

preprint2012arXiv

Detectability of Earth-like Planets in Circumstellar Habitable Zones of Binary Star Systems with Sun-like Components

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.

preprint2012arXiv

Discovery of Main-Belt Comet P/2006 VW139 by Pan-STARRS1

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.

preprint2012arXiv

GJ 581 update: Additional Evidence for a Super-Earth in the Habitable Zone

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 χ_ν^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σ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 χ_ν^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

preprint2012arXiv

High-Mass, Four-Planet Models for HR 8799: Constraining the Orbital Inclination and Age of the 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° and 13°, and 41, 46, and 31 Myr at 18°, 23°, and 30°, 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.

preprint2012arXiv

Kepler-47: A Transiting Circumbinary Multi-Planet System

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.

preprint2012arXiv

Observational and Dynamical Characterization of Main-Belt Comet P/2010 R2 (La Sagra)

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.

preprint2012arXiv

On The Effect of Giant Planets on the Scattering of Parent Bodies of Iron Meteorite from the Terrestrial Planet Region into the Asteroid Belt: A Concept Study

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.

preprint2012arXiv

On the HU Aquarii planetary system hypothesis

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.

preprint2012arXiv

The Lick-Carnegie Survey: A New Two-Planet System Around the Star HD 207832

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.

preprint2012arXiv

The Neptune-Sized Circumbinary Planet Kepler-38b

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.

preprint2011arXiv

New Light-Travel Time Models and Orbital Stability Study of the Proposed Planetary System HU Aquarii

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 $χ_{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 $χ_{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.

preprint2011arXiv

On the Detection of (Habitable) Super-Earths Around Low-Mass Stars Using Kepler and Transit Timing Variation Method

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.

preprint2011arXiv

Optical and Dynamical Characterization of Comet-Like Main-Belt Asteroid (596) Scheila

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.

preprint2011arXiv

Super-Earths: A New Class of Planetary Bodies

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.

preprint2010arXiv

A Survey of M Stars in the Field of View of Kepler Space Telescope

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.

preprint2010arXiv

Implications of the TTV-Detection of Close-In Terrestrial Planets Around M Stars for Their Origin and Dynamical Evolution

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.

preprint2010arXiv

The Lick-Carnegie Exoplanet Survey: A 3.1 M_Earth Planet in the Habitable Zone of the Nearby M3V Star Gliese 581

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.

preprint2010arXiv

The Lick-Carnegie Exoplanet Survey: A Saturn-Mass Planet in the Habitable Zone of the Nearby M4V Star HIP 57050

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.

preprint2010arXiv

The Lick-Carnegie Exoplanet Survey: A Uranus-mass Fourth Planet for GJ 876 in an Extrasolar Laplace Configuration

(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, $φ_{\rm Laplace}=λ_c-3λ_b+2λ_e$, for the Laplace resonance librates with an amplitude of $Δφ_{\rm Laplace}=40\pm13^{\circ}$ about $φ_{\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 $φ_{\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.''

preprint2010arXiv

The Lick-Carnegie Survey: Four New Exoplanet Candidates

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.

preprint2009arXiv

Dynamical Constraints on the Origin of Main Belt Comets

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.

preprint2009arXiv

Dynamics, Origin, and Activation of Main Belt Comets

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.

preprint2009arXiv

Planetary Dynamics and Habitable Planet Formation In Binary Star Systems

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.

preprint2009arXiv

The Stability and Prospects of the Detection of Terrestrial/Habitable Planets in Multiplanet and Multiple Star Systems

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.

preprint2007arXiv

Extreme Habitability: Formation of Habitable Planets in Systems with Close-in Giant Planets and/or Stellar Companions

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.

preprint1999arXiv

Partial Averaging Near a Resonance in Planetary Dynamics

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.