Search

The Kepler high-precision planetary sample has revealed a radius valley, separating compact super-Earths from sub-Neptunes with lower density. Super-Earths are generally assumed to be rocky planets that were probably born in-situ, while the composition and origin of sub-Neptunes remains debated. To provide more constraints on the formation history and composition, based on the planetary sample of Kepler multiple planet systems, we derive the distributions of orbital period ratios of sub-Neptune and super-Earth planet pairs and calculate the normalised fraction of near-first-order mean motion resonances. Using synthetic planetary systems generated by the Generation III Bern Model, we also obtain theoretical predictions of period ratio distributions of planet pairs of different compositions and origins. We find that actual Kepler sub-Neptune pairs show a normalised fraction smaller (larger) than the model predictions for water-rich (water-poor) pairs with confidence levels of about two sigma. The derived normalised fraction of actual Kepler Super-Earth pairs is generally consistent with that of water-poor model planet pairs but significantly smaller than that of synthetic water-rich planet pairs. Based on the distributions of orbital period ratios, we conclude that orbital migration has been more important for sub-Neptunes than for super-Earths, suggesting a partial ex situ formation of the former and an origin of the radius valley caused in part by distinct formation pathways. However, the model comparisons also show that sub-Neptunes in actual Kepler multiple systems are not likely to be all water-rich/ex situ planets but a mixture of the two (in situ/ex situ) pathways. Whereas, Kepler super-Earth planets are predominantly composed by of water-poor planets that were born inside the ice line, likely through a series of giant impacts without large scale migration., Comment: Accepted for publication in A&A; 14 pages, 6 figures in the main text, 7 figures in Appendix

Accurately characterizing intrinsic stellar photometric noise induced by stellar astrophysics, such as stellar activity, granulation, and oscillations, is of crucial importance for detecting transiting exoplanets. In this study, we investigate the relation between the intrinsic stellar photometric noise, as quantified by the Kepler rrmsCDPP measurement, and the level of stellar chromospheric activity, as indicated by the S-index of Ca II HK lines derived from the LAMOST spectra. Our results reveal a clear positive correlation between S-index and rrmsCDPP, and the correlation becomes more significant at higher activity levels and on longer timescales. We have therefore built an empirical relation between rrmsCDPP and S-index as well as Teff, logg, [Fe/H], and apparent magnitude with the XGBoost regression algorithm, using the LAMOST-Kepler common star sample as the training set. This method achieves a precision of ~20 ppm for inferring the intrinsic noise from the S-index and other stellar labels on a 6-hour integration duration. We have applied this empirical relation to the full LAMOST DR7 spectra database, and obtained the intrinsic noise predictions for 1,358,275 stars. The resultant catalog is publicly available and expected to be valuable for optimizing target selection for future exoplanet-hunting space missions, such as the Earth 2.0 mission., Comment: Accepted for publication in ApJS

The unexpected discovery of hot Jupiters challenged the classical theory of planet formation inspired by our solar system. Until now, the origin and evolution of hot Jupiters are still uncertain. Determining their age distribution and temporal evolution can provide more clues into the mechanism of their formation and subsequent evolution. Using a sample of 383 giant planets around Sun-like stars collected from the kinematic catalogs of the Planets Across Space and Time (PAST) project, we find that hot Jupiters are preferentially hosted by relatively younger stars in the Galactic thin disk. We subsequently find that the frequency of hot Jupiters declines with age. In contrast, the frequency of warm/cold Jupiters shows no significant dependence on age. Such a trend is expected from the tidal evolution of hot Jupiters' orbits, and our result offers supporting evidence using a large sample. We also perform a joint analysis on the planet frequencies in the stellar age-metallicity plane. The result suggests that the frequencies of hot Jupiters and warm/cold Jupiters, after removing the age dependence are both correlated with stellar metallicities. Moreover, we show that the above correlations can explain the bulk of the discrepancy in hot Jupiter frequencies inferred from the transit and radial velocity (RV) surveys, given that RV targets tend to be more metal-rich and younger than transits., Comment: Published in PNAS; 7 pages, 5 figures in the main text; 17 pages, 29 figures, 5 tables in the supporting information

One of the fundamental questions in astronomy is how planetary systems form and evolve. Measuring the planetary occurrence and architecture as a function of time directly addresses this question. In the fourth paper of the Planets Across Space and Time (PAST) series, we investigate the occurrence and architecture of Kepler planetary systems as a function of kinematic age by using the LAMOST-Gaia-Kepler sample. To isolate the age effect, other stellar properties (e.g., metallicity) have been controlled. We find the following results. (1) The fraction of stars with Kepler-like planets ($F_{\text{Kep}}$) is about 50% for all stars; no significant trend is found between $F_{\text{Kep}}$ and age. (2) The average planet multiplicity ($\bar{N}_p$) exhibits a decreasing trend (~2$\sigma$ significance) with age. It decreases from $\bar{N}_p$~3 for stars younger than 1 Gyr to $\bar{N}_p$~1.8 for stars about 8 Gyr. (3) The number of planets per star ($\eta=F_{\text{Kep}}\times\bar{N}_p$) also shows a decreasing trend (~2-3$\sigma$ significance). It decreases from $\eta$~1.6-1.7 for young stars to $\eta$~1.0 for old stars. (4) The mutual orbital inclination of the planets ($\sigma_{i,k}$) increases from $1.2^{+1.4}_{-0.5}$ to $3.5^{+8.1}_{-2.3}$ as stars aging from 0.5 to 8 Gyr with a best fit of $\log{\sigma_{i,k}}=0.2+0.4\times\log{\frac{\text{Age}}{\text{1Gyr}}}$. Interestingly, the Solar System also fits such a trend. The nearly independence of $F_{\text{Kep}}$~50% on age implies that planet formation is robust and stable across the Galaxy history. The age dependence of $\bar{N}_p$ and $\sigma_{i,k}$ demonstrates planetary architecture is evolving, and planetary systems generally become dynamically hotter with fewer planets as they age., Comment: 27 pages, 20 figures, 4tables, accepted for publication in AJ

Orbital eccentricity is one of the basic planetary properties, whose distribution may shed light on the history of planet formation and evolution. Here, in a series of works on Planetary Orbit Eccentricity Trends (dubbed POET), we study the distribution of planetary eccentricities and their dependence on stellar/planetary properties. In this paper, the first work of the POET series, we investigate whether and how the eccentricities of small planets depend on stellar metallicities (e.g., [Fe/H]). Previous studies on giant planets have found a significant correlation between planetary eccentricities and their host metallicities. Nevertheless, whether such a correlation exists in small planets (e.g. super-Earth and sub-Neptune) remains unclear. Here, benefiting from the large and homogeneous LAMOST-Gaia-Kepler sample, we characterize the eccentricity distributions of 244 (286) small planets in single (multiple) transiting systems with the transit duration ratio method. We confirm the eccentricity-metallicity trend that eccentricities of single small planets increase with stellar metallicities. Interestingly, a similar trend between eccentricity and metallicity is also found in the radial velocity (RV) sample. We also found that the mutual inclination of multiple transiting systems increases with metallicity, which predicts a moderate eccentricity-metallicity rising trend. Our results of the correlation between eccentricity (inclination) and metallicity for small planet support the core accretion model for planet formation, and they could be footprints of self (and/or external) excitation processes during the history of planet formation and evolution., Comment: 21 pages, 15 figures

We propose to develop a wide-field and ultra-high-precision photometric survey mission, temporarily named "Earth 2.0 (ET)". This mission is designed to measure, for the first time, the occurrence rate and the orbital distributions of Earth-sized planets. ET consists of seven 30cm telescopes, to be launched to the Earth-Sun's L2 point. Six of these are transit telescopes with a field of view of 500 square degrees. Staring in the direction that encompasses the original Kepler field for four continuous years, this monitoring will return tens of thousands of transiting planets, including the elusive Earth twins orbiting solar-type stars. The seventh telescope is a 30cm microlensing telescope that will monitor an area of 4 square degrees toward the galactic bulge. This, combined with simultaneous ground-based KMTNet observations, will measure masses for hundreds of long-period and free-floating planets. Together, the transit and the microlensing telescopes will revolutionize our understandings of terrestrial planets across a large swath of orbital distances and free space. In addition, the survey data will also facilitate studies in the fields of asteroseismology, Galactic archeology, time-domain sciences, and black holes in binaries., Comment: 116 pages,79 figures

The radius valley, a dip in the radius distribution of exoplanets at ~1.9 Earth radii separates compact rocky Super-Earths and Sub-Neptunes with lower density. Various hypotheses have been put forward to explain the radius valley. Characterizing the radius valley morphology and its correlation to stellar properties will provide crucial observation constraints on its origin mechanism and deepen the understanding of planet formation and evolution. In this paper, the third part of the Planets Across the Space and Time (PAST) series, using the LAMOST-Gaia-Kepler catalog, we perform a systematical investigation into how the radius valley morphology varies in the Galactic context, i.e., thin/thick galactic disks, stellar age and metallicity abundance ([Fe/H] and [alpha/Fe]). We find that (1) The valley becomes more prominent with the increase of both age and [Fe/H]. (2) The number ratio of super-Earths to sub-Neptunes monotonically increases with age but decreases with [Fe/H] and [alpha/Fe]. (3) The average radius of planets above the valley (2.1-6 Earth radii) decreases with age but increases with [Fe/H]. (4) In contrast, the average radius of planets below the valley (R < 1.7 Earth radii) is broadly independent on age and metallicity. Our results demonstrate that the valley morphology as well as the whole planetary radius distribution evolves on a long timescale of giga-years, and metallicities (not only Fe but also other metal elements, e.g., Mg, Si, Ca, Ti) play important roles in planet formation and in the long term planetary evolution., Comment: Accepted for pubilication in AJ, 20 Pages, 10 figures, 2 Tables (Appendix: 13 Figures)

Although the sample of exoplanets in binaries has been greatly expanded, the sample heterogeneity and observational bias are obstacles toward a clear figure of exoplanet demographics in the binary environment. To overcome the obstacles, we conduct a statistical study that focuses on S-type planetary systems detected by the Radial Velocity method. We try to account for observational biases by estimating, from available RV data, planet detection efficiencies for each individual system. Our main results are as follows: (1) Single (multiple) planetary systems are mostly found in close (wide) binaries with separation aB<(>)100-300 AU. (2) In binaries, single and multiple-planet systems are similar in 1-D distributions of mass and period as well as eccentricity (in contrast to the "eccentricity dichotomy" found in single-star systems) but different in the 2-D period-mass diagram. Specifically, there is a rectangular-shaped gap in the period-mass diagram of single-planet systems but not for multiples. This gap also depends on binary separation and is more prominent in close binaries. (3) There is a rising upper envelope in the period-mass diagram for planets in wide binaries as well as in single stars but not in close binaries. More specifically, there is a population of massive short period planets in close binaries but almost absent in wide binaries or single stars. We suggest that enhanced planetary migration, collision, and /or ejection in close binaries could be the potential underlying explanation for these three features., Comment: 24 pages, 15 figures, accepted to AJ

The Kepler telescope has discovered over 4,000 planets (candidates) by searching ? 200,000 stars over a wide range of distance (order of kpc) in our Galaxy. Characterizing the kinematic properties (e.g., Galactic component membership and kinematic age) of these Kepler targets (including the planet (candidate) hosts) is the first step towards studying Kepler planets in the Galactic context, which will reveal fresh insights into planet formation and evolution. In this paper, the second part of the Planets Across the Space and Time (PAST) series, by combining the data from LAMOST and Gaia and then applying the revised kinematic methods from PAST I, we present a catalog of kinematic properties(i.e., Galactic positions, velocities, and the relative membership probabilities among the thin disk, thick disk, Hercules stream, and the halo) as well as other basic stellar parameters for 35,835 Kepler stars. Further analyses of the LAMOST-Gaia-Kepler catalog demonstrate that our derived kinematic age reveals the expected stellar activity-age trend. Furthermore, we find that the fraction of thin(thick) disk stars increases (decreases) with the transiting planet multiplicity (Np = 0, 1, 2 and 3+) and the kinematic age decreases with Np, which could be a consequence of the dynamical evolution of planetary architecture with age. The LAMOST-Gaia-Kepler catalog will be useful for future studies on the correlations between the exoplanet distributions and the stellar Galactic environments as well as ages., Comment: Accepted for publication in Astronomical Journal, 19 pages, 9 figures, 5 tables. The full data of the LAMOST-Gaia-Kepler kinematic catalog (Table 3) will be published with the journal. If you need to use the data, please contact us via the email (jwxie@nju.edu.cn, dg1826003@smail.nju.edu.cn)

Over 4,000 exoplanets have been identified and thousands of candidates are to be confirmed. The relations between the characteristics of these planetary systems and the kinematics, Galactic components, and ages of their host stars have yet to be well explored. Aiming to addressing these questions, we conduct a research project, dubbed as PAST (Planets Across Space and Time). To do this, one of the key steps is to accurately characterize the planet host stars. In this paper, the Paper I of the PAST series, we revisit the kinematic method for classification of Galactic components and extend the applicable range of velocity ellipsoid from about 100 pc to 1, 500 pc from the sun in order to cover most known planet hosts. Furthermore, we revisit the Age-Velocity dispersion Relation (AVR), which allows us to derive kinematic age with a typical uncertainty of 10-20% for an ensemble of stars. Applying the above revised methods, we present a catalog of kinematic properties (i.e. Galactic positions, velocities, the relative membership probabilities among the thin disk, thick disk, Hercules stream, and the halo) as well as other basic stellar parameters for 2,174 host stars of 2,872 planets by combining data from Gaia, LAMOST, APOGEE, RAVE, and the NASA exoplanet archive. The revised kinematic method and AVR as well as the stellar catalog of kinematic properties and ages lay foundation for future studies on exoplanets from two dimensions of space and time in the Galactic context., Comment: 26 pages, 20 figure, 10 tables, Accepted for publication in ApJ. The full data of the planet kinematic catalog (Table 10) will be published with the journal. If you need to use the data, please contact us via the email (jwxie@nju.edu.cn, dg1826003@smail.nju.edu.cn)

The Kepler space mission has detected a large number of exoplanets in multiple transiting planet systems. Previous studies found that these Kepler multiple planet systems exhibit an intra-system uniformity, namely planets in the same system have similar sizes and correlated orbital spacings. However, it is important to consider the possible role of selection effects due to observational biases. In this paper, we revisit the orbital spacing aspect of the pattern after taking observational biases into account using a forward modeling method. We find that orbital spacings, in terms of period ratios, of Kepler multiple planet systems are significantly correlated only for those tightly packed systems, and the transition from correlation to non-correlation is abrupt with a boundary at mean period ratio $\overline{PR}$ $\sim 1.5 - 1.7$. In this regard, the pattern of orbital spacing is more like a dichotomy rather than a global correlation. Furthermore, we find that such an apparent orbital spacing dichotomy could be essentially a projection of a dichotomy that related to mean motion resonance (MMR), which we dub as MMR dichotomy, and itself could be a natural result of planet migration and dynamical evolution., Comment: Accepted for publication in Astronomical Journal. 11 figures

Since September 2018, LAMOST starts a new 5-year medium-resolution spectroscopic survey (MRS) using bright/gray nights. We present the scientific goals of LAMOST-MRS and propose a near optimistic strategy of the survey. A complete footprint is also provided. Not only the regular medium-resolution survey, but also a time-domain spectroscopic survey is being conducted since 2018 and will be end in 2023. According to the detailed survey plan, we expect that LAMOST-MRS can observe about 2 million stellar spectra with ~7500 and limiting magnitude of around G=15 mag. Moreover, it will also provide about 200 thousand stars with averagely 60-epoch observations and limiting magnitude of G~14 mag. These high quality spectra will give around 20 elemental abundances, rotational velocities, emission line profiles as well as precise radial velocity with uncertainty less than 1 km/s. With these data, we expect that LAMOST can effectively leverage sciences on stellar physics, e.g. exotic binary stars, detailed observation of many types of variable stars etc., planet host stars, emission nebulae, open clusters, young pre-main-sequence stars etc., Comment: 25 pages, 10 figues. Submitted to Reseach in Astronomy and Astrophysics

The Kepler mission has discovered thousands of exoplanets around various stars with different spectral types (M, K, G, and F) and thus different masses and effective temperatures. Previous studies have shown that the planet occurrence rate, in terms of average number of planets per star, drops with increasing stellar effective temperature (Teff). In this paper, with the final Kepler Data Release (DR25) catalog, we revisit the relation between stellar effective temperature (as well as mass) and planet occurrence, but in terms of the fraction of stars with planets and the number of planets per planetary system (i.e. planet multiplicity). We find that both the fraction of stars with planets and planet multiplicity decrease with increasing stellar temperature and mass. Specifically, about 75% late-type stars (Teff<5000 K) have Kepler-like planets with an average planet multiplicity of ~2.8; while for early-type stars (Teff>6500 K), this fraction and the average multiplicity fall down to ~35% and ~1.8, respectively. The decreasing trend in the fraction of stars with planets is very significant with ${\Delta}$AIC> 30, though the trend in planet multiplicity is somewhat tentative with ${\Delta}$AIC~5. Our results also allow us to derive the dispersion of planetary orbital inclinations in relationship with stellar effective temperature. Interestingly, it is found to be similar to the well-known trend between obliquity and stellar temperature, indicating that the two trends might have a common origin., Comment: 25 pages, 21 figures, 2 tables, accepted for publication in AJ

Planetary systems can survive the stellar evolution, as evidenced by the atmospheric metal pollution and dusty disks of single white dwarfs. Recent observations show that 1 to 4 percent of single white dwarfs are accompanied by dusty disks, while the occurrence rate of metal pollution is about 25 to 50 percent . The dusty disks and metal pollution have been associated with accretion of remnant planetary systems around white dwarfs, yet the relation between these two phenomena is still unclear. Here we suggest an evolutionary scenario to link the two observational phenomena. By analyzing a sample of metal polluted white dwarfs, we find that the mass accretion rate onto the white dwarf generally follows a broken power law decay, which matches well with the theoretical prediction, if assuming dust accretion is primarily driven by Poynting-Robertson drag and the dust source is primarily delivered via dynamically falling asteroids perturbed by a Jovian planet. The presence of disks is mainly at the early stage (about 0.1 to 0.7 billion years) of the whole process of metal pollution, which is detectable until about 8 billion years, naturally explaining the fraction (about 2-16 percent) of metal-polluted white dwarfs having dusty disks. The success of this scenario also implies that the configuration of an asteroid belt with an outer gas giant might be common around stars of several solar masses., Comment: 33 pages, 3 figures, 10 supplementary figures, 3 supplementary tables. Published in Nature Astronomy

The LAMOST-Kepler (LK-) project was initiated to use the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) to make spectroscopic follow-up observations for the targets in the field of the Kepler mission. The Kepler field is divided into 14 subfields that are adapted to the LAMOST circular field with a diameter of 5 degrees. During the regular survey phase of LAMOST, the LK-project took data from 2012 June to 2017 June and covered all the 14 subfields at least twice. In particular, we describe in this paper the second Data Release of the LK-project, including all spectra acquired through 2015 May to 2017 June together with the first round observations of the LK-project from 2012 June to 2014 September. The LK-project now counts 227 870 spectra of 156 390 stars, among which we have derived atmospheric parameters (log g, T eff and [Fe/H]) and heliocentric radial velocity (RV) for 173 971 spectra of 126 172 stars. These parameters were obtained with the most recent version of the LAMOST Stellar Parameter Pipeline v 2.9.7. Nearly one half, namely 76 283 targets, are observed both by LAMOST and Kepler telescopes. These spectra, establishing a large spectroscopy library, will be useful for the entire astronomical community, particularly for planetary science and stellar variability on Kepler targets., Comment: 13 pages, 5 figures, 3 tables; ApJS in press

During the evolution of proto-planetary disc, photo-evaporations of both central and external stars play important roles. Considering the complicated radiation surroundings in the clusters, where the star formed, the proto-planetary discs survive in different lifetimes due to flyby events. In this paper, we mainly focus on the disc around a T Tauri star, which encounters with another main-sequence star with different temperatures in hyperbolic orbits with different peri-center distances, eccentricities and inclinations. We find the criterion for gap-opening due to photo-evaporation of central star after the flyby event. A gap is opened in the late stage of gas disc, and induce that the gap only influence the planet formation and migration limitedly. If the flyby orbit has a moderate value of peri-center distance, which weakly depends on the eccentricity and inclination, the external photo-evaporation lead to a maximum mass loss during the flyby event. Flyby stars in orbits with smaller eccentricities or larger inclinations induce larger mass loss. Adopting a simple multiple flyby models, we conclude that in open clusters, gas discs usually survive in typical lifetimes between 1 and 10 Myr, except there are many massive stars in dense open clusters. In globular clusters, discs disperse very quickly and hardly produce the gas giant planets. The fast-depleted discs are probably responsible for the null detection of giant planets in globular clusters., Comment: accepted by MNRAS

Stars are commonly formed in binary systems, which provide a natural laboratory for studying planet formation in extreme conditions. In our first paper (Paper I) of a series Xie et al. (2011), we have shown that the intermediate stage - from planetesimals to planetary embryos/cores - of planet formation can proceed even in highly inclined binaries. Following Paper I, here we numerically study the late stage of terrestrial planet formation, i.e., from embryos to full planets, in binary systems of various orbital configurations. We identify an orbital alignment or anti-alignment effect; namely, although an inclined binary generally misaligns the planetary orbits with respect to the spin axis of the primary host star (i.e., causing large obliquity), it could align or anti-align the planetary orbits with respect to the binary orbit. Such an orbital (anti-)alignment effect is caused by the combination of orbital differential precession and self-damping, and it is mostly significant in cases of intermediate binary separations, i.e., a_B~40-200 AU for terrestrial planet formation around 1 AU from the primary stars. In such intermediate separation binaries, somewhat contrary to intuition, the binary companion can aid planet growth by having increased the rate of collisions, forming significantly more massive but fewer planets. In the other two ends, the companion is ether too close thus plays a violently disruptive role or too wide to have significant effect on planet formation. Future observations, which can discover more planet-bearing binary star systems and constrain their masses and 3-D orbital motions will test our numerical findings., Comment: Accepted by ApJ. Main article: 9 pages and 5 figures. Appendix: 8 pages, 6 figures and 1 table

The high planetary multiplicity revealed by Kepler implies that Transit Time Variations (TTVs) are intrinsically common. The usual procedure for detecting these TTVs is biased to long-period, deep transit planets whereas most transiting planets have short periods and shallow transits. Here we introduce the Spectral Approach to TTVs technique that allows expanding the TTVs catalog towards lower TTV amplitude, shorter orbital period, and shallower transit depth. In the Spectral Approach we assume that a sinusoidal TTV exists in the data and then calculate the improvement to $\chi^2$ this model allows over that of linear ephemeris model. This enables detection of TTVs even in cases where the transits are too shallow so individual transits cannot be timed. The Spectral Approach is more sensitive due to the reduced number of free parameters in its model. Using the Spectral Approach, we: (a) detect 131 new periodic TTVs in Kepler data (an increase of ~2/3 over a previous TTV catalog); (b) Constrain the TTV periods of 34 long-period TTVs and reduce amplitude errors of known TTVs; (c) Identify cases of multi-periodic TTVs, for which absolute planetary mass determination may be possible. We further extend our analysis by using perturbation theory assuming small TTV amplitude at the detection stage, which greatly speeds up our detection (to a level of few seconds per star). Our extended TTVs sample shows no deficit of short period or low amplitude transits, in contrast to previous surveys in which the detection schemes were significantly biased against such systems., Comment: ApJS accepted. 23 pages, 40 figures

We discover a population of short-period, Neptune-size planets sharing key similarities with hot Jupiters: both populations are preferentially hosted by metal-rich stars, and both are preferentially found in Kepler systems with single transiting planets. We use accurate LAMOST DR4 stellar parameters for main-sequence stars to study the distributions of short-period 1d < P < 10d Kepler planets as a function of host star metallicity. The radius distribution of planets around metal-rich stars is more "puffed up" as compared to that around metal-poor hosts. In two period-radius regimes, planets preferentially reside around metal-rich stars, while there are hardly any planets around metal-poor stars. One is the well-known hot Jupiters, and the other is a population of Neptune-size planets (2 R_Earth <~ R_p <~ 6 R_Earth), dubbed as "Hoptunes". Also like hot Jupiters, Hoptunes occur more frequently in systems with single transiting planets though the fraction of Hoptunes occurring in multiples is larger than that of hot Jupiters. About 1% of solar-type stars host "Hoptunes", and the frequencies of Hoptunes and hot Jupiters increase with consistent trends as a function of [Fe/H]. In the planet radius distribution, hot Jupiters and Hoptunes are separated by a "valley" at approximately Saturn size (in the range of 6 R_Earth <~ R_p <~ 10 R_Earth), and this "hot-Saturn valley" represents approximately an order-of-magnitude decrease in planet frequency compared to hot Jupiters and Hoptunes. The empirical "kinship" between Hoptunes and hot Jupiters suggests likely common processes (migration and/or formation) responsible for their existence., Comment: Published in PNAS. Title modified and text slightly adjusted according to PNAS style. Two ancillary files: 1) Table S1 (hoptune_tableS1.csv) with planetary and stellar parameters of the sample; 2) Table S2 (hoptune_tableS2.csv) with stellar parameters of 32556 LAMOST-Kepler targets from LAMOST DR3

Many multiple-planet systems have been found by the Kepler transit survey and various radial velocity (RV) surveys. Kepler planets show an asymmetric feature, namely, there are small but significant deficits/excesses of planet pairs with orbital period spacing slightly narrow/wide of the exact resonance, particularly near the first order mean motion resonance (MMR), such as 2:1 and 3:2 MMR. Similarly, if not exactly the same, an asymmetric feature (pileup wide of 2:1 MMR) is also seen in RV planets, but only for massive ones. We analytically and numerically study planets' orbital evolutions near and in the MMR. We find that their orbital period ratios could be asymmetrically distributed around the MMR center regardless of dissipation. In the case of no dissipation, Kepler planets' asymmetric orbital distribution could be partly reproduced for 3:2 MMR but not for 2:1 MMR, implying that dissipation might be more important to the latter. The pileup of massive RV planets just wide of 2:1 MMR is found to be consistent with the scenario that planets formed separately then migrated toward the MMR. The location of the pileup infers a K value of 1-100 on the order of magnitude for massive planets, where K is the damping rate ratio between orbital eccentricity and semimajor axis during planet migration.

Many exoplanets have been found in orbits close to their host stars and thus they are subject to the effects of photo-evaporation. Previous studies have shown that a large portion of exoplanets detected by the Kepler mission have been significantly eroded by photo-evaporation. In this paper, we numerically study the effects of photo-evaporation on the orbital evolution of a hypothesized moon system around a planet. We find that photo-evaporation is crucial to the stability of the moon system. Photo-evaporation can erode the atmosphere of the planet thus leading to significant mass loss. As the planet loses mass, its Hill radius shrinks and its moons increase their orbital semi-major axes and eccentricities. When some moons approach their critical semi-major axes, global instability of the moon system would be triggered, which usually ends up with two, one or even zero surviving moons. Some lost moons could escape from the moon system to become a new planet orbiting the star or run away further to become a free-floating object in the Galaxy. Given the destructive role of photo-evaporation, we speculate that exo-moons are less common for close-in planets ($<0.1$ AU), especially those around M-type stars, because they are more X-ray luminous and thus enhancing photo-evaporation. The lessons we learn in this study may be helpful for the target selection of on-going/future exomoon searching programs., Comment: accepted by ApJ , 14 pages, 9 figures, 2 tables

The nearly circular (mean eccentricity ~0.06) and coplanar (mean mutual inclination ~3 deg) orbits of the Solar System planets motivated Kant and Laplace to put forth the hypothesis that planets are formed in disks, which has developed into the widely accepted theory of planet formation. Surprisingly, the first several hundred extrasolar planets (mostly Jovian) discovered using the Radial Velocity (RV) technique are commonly on eccentric orbits ( ~ 0.3). This raises a fundamental question: Are the Solar System and its formation special? The Kepler mission has found thousands of transiting planets dominated by sub-Neptunes, but most of their orbital eccentricities remain unknown. By using the precise spectroscopic host star parameters from the LAMOST observations, we measure the eccentricity distributions for a large (698) and homogeneous Kepler planet sample with transit duration statistics. Nearly half of the planets are in systems with single transiting planets (singles), while the other half are multiple-transiting planets (multiples). We find an eccentricity dichotomy: on average, Kepler singles are on eccentric orbits with ~0.3, while the multiples are on nearly circular ( = 0.04^{+0.03}_{-0.04}) and coplanar ( = 1.4^{+0.8}_{-1.1} deg) orbits similar to the Solar System planets. Our results are consistent with previous studies of smaller samples and individual systems. We also show that Kepler multiples and solar system objects follow a common relation ~(1-2)x between mean eccentricities and mutual inclinations. The prevalence of circular orbits and the common relation may imply that the solar system is not so atypical in the galaxy after all., Comment: 2016 PNAS published, see http://www.pnas.org/content/early/2016/09/21/1604692113.abstract

The Kepler mission provides a wealth of multiple transiting planet systems (MTPS). The formation and evolution of multi-planet systems are likely to be influenced by companion stars given the abundance of multi stellar systems. We study the influence of stellar companions by measuring the stellar multiplicity rate of MTPS. We select 138 bright (KP < 13.5) Kepler MTPS and search for stellar companions with AO imaging data and archival radial velocity (RV) data. We obtain new AO images for 73 MTPS. Other MTPS in the sample have archival AO imaging data from the Kepler Community Follow-up Observation Program (CFOP). From these imaging data, we detect 42 stellar companions around 35 host stars. For stellar separation 1 AU < a < 100 AU, the stellar multiplicity rate is 5.2 $\pm$ 5.0% for MTPS, which is 2.8{\sigma} lower than 21.1 $\pm$ 2.8% for the control sample, i.e., the field stars in the solar neighborhood. We identify two origins for the deficit of stellar companions within 100 AU to MTPS: (1) a suppressive planet formation, and (2) the disruption of orbital coplanarity due to stellar companions. To distinguish between the two origins, we compare the stellar multiplicity rates of MTPS and single transiting planet systems (STPS). However, current data are not sufficient for this purpose. For 100 AU < a < 2000 AU, the stellar multiplicity rates are comparable for MTPS (8.0 $\pm$ 4.0%), STPS (6.4 $\pm$ 5.8%), and the control sample (12.5 $\pm$ 2.8%)., Comment: 13 pages, 3 figures, 5 tables, accepted by ApJ

As hundreds of gas giant planets have been discovered, we study how these planets form and evolve in different stellar environments, specifically in multiple stellar systems. In such systems, stellar companions may have a profound influence on gas giant planet formation and evolution via several dynamical effects such as truncation and perturbation. We select 84 Kepler Objects of Interest (KOIs) with gas giant planet candidates. We obtain high-angular resolution images using telescopes with adaptive optics (AO) systems. Together with the AO data, we use archival radial velocity data and dynamical analysis to constrain the presence of stellar companions. We detect 59 stellar companions around 40 KOIs for which we develop methods of testing their physical association. These methods are based on color information and galactic stellar population statistics. We find evidence of suppressive planet formation within 20 AU by comparing stellar multiplicity. The stellar multiplicity rate for planet host stars is 0$^{+5}_{-0}$\% within 20 AU. In comparison, the stellar multiplicity rate is 18\%$\pm$2\% for the control sample, i.e., field stars in the solar neighborhood. The stellar multiplicity rate for planet host stars is 34\%$\pm$8\% for separations between 20 and 200 AU, which is higher than the control sample at 12\%$\pm$2\%. Beyond 200 AU, stellar multiplicity rates are comparable between planet host stars and the control sample. We discuss the implications of the results to gas giant planet formation and evolution., Comment: 18 pages, 5 figures, accepted by ApJ

\We present the sixth catalog of Kepler candidate planets based on nearly 4 years of high precision photometry. This catalog builds on the legacy of previous catalogs released by the Kepler project and includes 1493 new Kepler Objects of Interest (KOIs) of which 554 are planet candidates, and 131 of these candidates have best fit radii <1.5 R_earth. This brings the total number of KOIs and planet candidates to 7305 and 4173 respectively. We suspect that many of these new candidates at the low signal-to-noise limit may be false alarms created by instrumental noise, and discuss our efforts to identify such objects. We re-evaluate all previously published KOIs with orbital periods of >50 days to provide a consistently vetted sample that can be used to improve planet occurrence rate calculations. We discuss the performance of our planet detection algorithms, and the consistency of our vetting products. The full catalog is publicly available at the NASA Exoplanet Archive., Comment: 18 pages, to be published in the Astrophysical Journal Supplement Series

Almost half of the stellar systems in the solar neighborhood are made up of multiple stars. In multiple-star systems, planet formation is under the dynamical influence of stellar companions, and the planet occurrence rate is expected to be different from that for single stars. There have been numerous studies on the planet occurrence rate of single star systems. However, to fully understand planet formation, the planet occurrence rate in multiple-star systems needs to be addressed. In this work, we {{infer}} the planet occurrence rate in multiple-star systems by measuring the stellar multiplicity rate for planet host stars. For a sub-sample of 56 $Kepler$ planet host stars, we use adaptive optics (AO) imaging and the radial velocity (RV) technique to search for stellar companions. The combination of these two techniques results in high search completeness for stellar companions. We detect 59 visual stellar companions to 25 planet host stars with AO data. {{Three stellar companions are within 2$^{\prime\prime}$, and 27 within 6$^{\prime\prime}$. We also detect 2 possible stellar companions (KOI 5 and KOI 69) showing long-term RV acceleration.}} After correcting for a bias against planet detection in multiple-star systems due to flux contamination, we find that planet formation is suppressed in multiple-star systems with separations smaller than 1500 AU. Specifically, we find that compared to single star systems, planets in multiple-star systems occur $4.5\pm3.2$, $2.6\pm1.0$, and $1.7\pm0.5$ times less frequently when a stellar companion is present at a distance of 10, 100, and 1000 AU, respectively. This conclusion applies only to circumstellar planets; the planet occurrence rate for circumbinary planets requires further investigation., Comment: 17 pages, 7 figures, 5 tables, accepted by ApJ

The planet occurrence rate for multiple stars is important in two aspects. First, almost half of stellar systems in the solar neighborhood are multiple systems. Second, the comparison of the planet occurrence rate for multiple stars to that for single stars sheds light on the influence of stellar multiplicity on planet formation and evolution. We developed a method of distinguishing planet occurrence rates for single and multiple stars. From a sample of 138 bright (KP <13.5) Kepler multi-planet candidate systems, we compared the stellar multiplicity rate of these planet host stars to that of field stars. Using dynamical stability analyses and archival Doppler measurements, we find that the stellar multiplicity rate of planet host stars is significantly lower than field stars for semi-major axes less than 20 AU, suggesting that planet formation and evolution are suppressed by the presence of a close-in companion star at these separations. The influence of stellar multiplicity at larger separations is uncertain because of search in- completeness due to a limited Doppler observation time baseline and a lack of high resolution imaging observation. We calculated the planet confidence for the sample of mutlti-planet candidates, and find that the planet confidences for KOI 82.01, KOI 115.01, KOI 282.01 and KOI 1781.02 are higher than 99.7% and thus validate the planetary nature of these four planet candidates. This sample of bright Kepler multi-planet candidates with refined stellar and orbital parameters, planet confidence estimation, and nearby stellar companion identification offers a well-characterized sample for future theoretical and observational study., Comment: 52 pages, 12 figures, 5 tables, accepted to ApJ

Following on from Paper I in our series (Xie 2013), we report the confirmation by Transit Timing Variations (TTVs) of a further 30 planets in 15 multiple planet systems, using the publicly available Kepler light curves (Q0-Q16). All of these fifteen pairs are near first-order Mean Motion Resonances (MMR), showing sinusoidal TTVs consistent with theoretically predicted periods, which demonstrate they are orbiting and interacting in the same systems. Although individual masses cannot be accurately extracted based only on TTVs (because of the well known degeneracy between mass and eccentricity), the measured TTV phases and amplitudes can still place relatively tight constraints on their mass ratios and upper limits on their masses, which confirm their planetary nature. Some of these systems (KOI-274, KOI-285, KOI-370 and KOI-2672) are relatively bright and thus suitable for further follow-up observations., Comment: Accepted to ApJS (with assigned Kepler numbers and newly revised stellar properties from Huber et al. 2013) . 13 pages, 7 figures, 3 tables. Text follows the same style and has substantial overlap with Paper I in our series (Xie 2012, arXiv:1208.3312)

A transiting planet exhibits sinusoidal transit-time-variations (TTVs) if perturbed by a companion near a mean-motion-resonance (MMR). We search for sinusoidal TTVs in more than 2600 Kepler candidates, using the publicly available Kepler light-curves (Q0-Q12). We find that the TTV fractions rise strikingly with the transit multiplicity. Systems where four or more planets transit enjoy roughly five times higher TTV fraction than those where a single planet transits, and about twice higher than those for doubles and triples. In contrast, models in which all transiting planets arise from similar dynamical configurations predict comparable TTV fractions among these different systems. One simple explanation for our results is that there are at least two different classes of Kepler systems, one closely packed and one more sparsely populated., Comment: 7 pages, 5 figures, 1 table, accepted to ApJ

Compared to the giant planets in the solar system, exoplanets have many remarkable properties such as the prevalence of giant planets on eccentric orbits and the presence of hot Jupiters. Planet-planet scattering (PPS) between giant planets is a possible mechanism in interpreting above and other observed properties. If the observed giant planet architectures are indeed the outcomes of PPS, such drastic dynamical process must affect their primordial moon systems. In this Letter, we discuss the effect of the PPS on the survival of their regular moons. From the viewpoint of observations, some preliminary conclusions are drawn from the simulations. 1. PPS is a destructive process to the moon systems, single planets on eccentric orbits are not the ideal moon-search targets. 2. If hot Jupiters formed through PPS, their original moons have little chance to survive. 3. Planets in multiple systems with small eccentricities are more likely holding their primordial moons. 4. Compared to the lower-mass planets, the massive ones in multiple systems may not be the preferred moon-search targets if the system underwent a PPS history., Comment: published at ApJL

According to the core accretion theory, circumbinary embryos can form only beyond a critical semimajor axis (CSMA). However, due to the relatively high density of solid materials in the inner disk, significant amount of small planetesimals must exist in the inner zone when embryos were forming outside this CSMA. So embryos migration induced by the planetesimal swarm is possible after the gas disk depletion. Through numerical simulations, we found (i) the scattering-driven inward migration of embryos is robust, planets can form in the habitable zone if we adopt a mass distribution of MMSN-like disk; (ii) the total mass of the planetesimals in the inner region and continuous embryo-embryo scattering are two key factors that cause significant embryo migrations; (iii) the scattering-driven migration of embryos is a natural water-deliver mechanism. We propose that planet detections should focus on the close binary with its habitable zone near CSMA., Comment: Accepted for publication at ApJL

We extract Transit Timing Variation (TTV) signals for 12 pairs of transiting planet candidates that are near first-order Mean Motion Resonances (MMR), using publicly available Kepler light curves (Q0-Q14). These pairs show significant sinusoidal TTVs with theoretically predicted periods, which demonstrate these planet candidates are orbiting and interacting in the same system. Although individual masses cannot be accurately extracted based only on TTVs because of the well known degeneracy between mass and eccentricity, TTV phases and amplitudes can still place upper limits on the masses of the candidates, confirming their planetary nature. Furthermore, the mass ratios of these planet pairs can be relatively tight constrained using these TTVs. The planetary pair in Kepler-82 (KOI-880) seems to have a particularly high mass ratio and density ratio, which might indicate very different internal compositions of these two planets. Some of these newly confirmed planets are also near MMR with other candidates in the system, forming unique resonance chains, e.g., Kepler-80 (KOI-500)., Comment: Accepted to ApJS. 17 pages, 11 figures, 3 tables. KOI 869 is replaced with KOI 2038. Kepler numbers are assigned

The metallic gas associated with the Beta Pic debris disk is not believed to be primordial, but arises from the destruction of dust grains. Recent observations have shown that carbon and oxygen in this gas are exceptionally overabundant compared to other elements, by some 400 times. We study the origin of this enrichment under two opposing hypothesis, preferential production, where the gas is produced with the observed unusual abundance (as may happen if gas is produced by photo-desorption from C/O-rich icy grains), and preferential depletion, where the gas evolves to the observed state from an original solar abundance (if outgassing occurs under high-speed collisions) under a number of dynamical processes. We include ... ... We find ... ..., Comment: accepted to ApJ, revised text is highlighted in red color, the old figure 2 is removed, other figures are revised

We investigate the effects of gas-disk gravity on the planetesimal dynamics in inclined binary systems, where the circumprimary disk plane is tilted by a significant angle ($i_B$) with respect to the binary disk plane. Our focus is on the Lidov-Kozai mechanism and the evolution of planetesimal eccentricity and inclination. Using both analytical and numerical methods, we find that, on one hand, the disk gravity generally narrows down the Kozai-on region, i.e., the Lidov-Kozai effect can be suppressed in certain parts of (or even the whole of) the disk, depending on various parameters. In the Kozai-off region, planetesimals would move on orbits close to the mid-plane of gas-disk, with the relative angle ($i^{'}$) following a small amplitude periodical oscillation. On the other hand, when we include the effects of disk gravity, we find that the Lidov-Kozai effect can operate even at arbitrarily low inclinations ($i_B$), although lower $i_B$ leads to a smaller Kozai-on region. Furthermore, in the Kozai-on region, most planetesimals' eccentricities can be excited to extremely high values ($\sim 1$), and such extreme high eccentricities usually accompany orbital flipping, i.e., planetesimal orbit flips back and forth between anterograde and retrograde. Once a planetesimal reaches very high orbital eccentricity, gas drag damping will shrink the planetesimal orbit, forming a "hot planetesimal" on a near circular orbit very close to the primary star. Such a mechanism, if replacing the planetesimals and gas drag damping with Jupiters and tidal damping respectively, may lead to frequent production of hot-Jupiters., Comment: 11 pages, 10 figures, Accepted to ApJ

Most detected planet-bearing binaries are in wide orbits, for which a high inclination, $i_B$, between the binary orbital plane and the plane of the planetary disk around the primary is likely to be common. In this paper, we investigate the intermediate stages - from planetesimals to planetary embryos/cores - of planet formation in such highly inclined cases. Our focus is on the effects of gas drag on the planetesimals' orbital evolution, in particular on the evolution of the planetesimals' semimajor axis distribution and their mutual relative velocities. We first demonstrate that a non-evolving axisymmetric disk model is a good approximation for studying the effects of gas drag on a planetesimal in the highly inclined case ($30^\circ

We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup., Comment: 47 pages, 11 figures. Accepted for publication in ApJ

The standard model of planet formation considers an initial phase in which planetesimals form from a dust disk, followed by a phase of mutual planetesimal-planetesimal collisions, leading eventually to the formation of planetary embryos. However, there is a potential transition phase (which we call the "snowball phase"), between the formation of the first planetesimals and the onset of mutual collisions amongst them, which has often been either ignored or underestimated in previous studies. In this snowball phase, isolated planetesimals move on Keplerian orbits and grow solely via the direct accretion of sub-cm sized dust entrained with the gas in the protoplanetary disk. Using a simplified model in which planetesimals are progressively produced from the dust, we consider the expected sizes to which the planetesimals can grow before mutual collisions commence and derive the dependence of this size on a number of critical parameters, including the degree of disk turbulence, the planetesimal size at birth and the rate of planetesimal creation. For systems in which turbulence is weak and the planetesimals are created at a low rate and with relatively small birth size, we show that the snowball growth phase can be very important, allowing planetesimals to grow by a factor of 10^6 in mass before mutual collisions take over. In such cases, the snowball growth phase can be the dominant mode to transfer mass from the dust to planetesimals. Moreover, such growth can take place within the typical lifetime of a protoplanetary gas disk. A noteworthy result is that ... ...(see the paper). For the specific case of close binaries such as Alpha Centauri ... ... (see the paper). From a more general perspective, these preliminary results suggest that an efficient snowball growth phase provides a large amount of "room at the bottom" for theories of planet formation., Comment: Accepted for publication in the Astrophysical Journal. 15 pages, 4 figures, 1 table

Stellar perturbations affect planet-formation in binary systems. Recent studies show that the planet-formation stage of mutual accretion of km-sized planetesimals is most sensitive to binary effects. In this paper, the condition for planetesimal accretion is investigated around Alpha CenB, which is believed to be an ideal candidate for detection of an Earth-like planet in or near its habitable zone(0.5-0.9 AU). A simplified scaling method is developed to estimate the accretion timescale of the planetesimals embedded in a protoplanetary disk. Twenty-four cases with different binary inclinations(i_B=0, 0.1, 1.0, and 10 deg), gas densities(0.3,1,and 3 times of the Minimum Mass of Solar Nebula, MMSN hereafter), and with and without gas depletion, are simulated. We find: (1)...(2)...(3)...(4)...(see the paper for details). In other words, our results suggest that formation of Earth-like planets through accretion of km-sized planetesimals is possible in Alpha CenB, while formation of gaseous giant planets is not favorable., Comment: 16 pages 7 figures, ApJ 708, 1566

Recent observations show that planet can reside in close binary systems with stellar separation of only about 20 AU. However, planet formation in such close binary systems is a challenge to current theory. One of the major theoretical problems occurs in the intermediate stage-planetesimals accretion into planetary embryos-during which the companion's perturbations can stir up the relative velocites(dV) of planetesimals and thus slow down or even cease their growth. However, all previous studies assumed a 2-dimentional (2D) disk and a coplanar binary orbit. Extending previous studies by including a 3D gas disk and an inclined binary orbit with small relative inclination of i_B=0.1-5 deg, we numerically investigate the conditions for planetesimal accretion at 1-2 AU, an extension of the habitable zone(1-1.3 AU), around alpha Centauri A in this paper. Inclusion of the binary inclination leads to: (1) differential orbital phasing is realized in the 3D space, and thus different-sized bodies are separated from each other; (2) total impact rate becomes lower, and impacts mainly occur between similar-sized bodies; (3) accretion is more favored, but the balance between accretion and erosion remains uncertain, and the "possible accretion region" extends up to 2AU when assuming an optimistic Q*(critical specific energy that leads to catastrophic fragmentation); and (4) impact velocities (dV) are significantly reduced but still much larger than their escape velocities, which infers that planetesimals grow by means of type II runaway mode. As a conclusion, inclusion of a small binary inclination is a promising mechanism that favors accretion, opening a possibility that planet formation in close binary systems can go through the difficult stage of planetesimals accretion into planetary embryos., Comment: 9 pages, 6 figures, ApJ, 698, 2066

Currently, one of major problems concerning planet formation theory in close binary systems is, the strong perturbation from the companion star can increase relative velocities ($\triangle V$) of planetesimals around the primary and thus hinder their growth. According to previous studies, while gas drag can reduce the $\triangle V$ between bodies of the same sizes by forcing orbital alignment to planetesimals, it increases the $\triangle V$ among bodies of different sizes. In this paper, focusing on the $\gamma$ Cephei binary system, we propose a mechanism that can overcome this difficulty. We show that in a dissipating gas disk (with a typical dissipating timescale of $\sim 10^5-10^6$ years), all the planetesimals eventually converge towards the same forced orbits regardless of their sizes, leading to much lower impact velocities among them. These $\triangle V$ decrease processes progressively increase net mass accretion and even trigger runaway growth for large bodies (radius $>15$ km). The effect of size distribution of planetesimals is discussed, and found to be one of the dominant factors that determine the outcome of collisional evolution. Anyway, it can be concluded that by including the gas dissipation in the early stage of disk evolution, the conditions for planetesimal accretion become much better, and the process from planetesimal to planet-embryo can be carried out in close binary systems like $\gamma$ Cephei., Comment: 27 pages, 9 figures, accepted by ApJ