Source author record

Ji-Lin Zhou

Ji-Lin Zhou 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

38works
5topics
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

38 published item(s)

preprint2023arXiv

Correcting Stellar Flare Frequency Distributions Detected by TESS and Kepler

The habitability of planets is closely connected with the stellar activity, mainly the frequency of flares and the distribution of flare energy. Kepler and TESS find many flaring stars are detected via precise time-domain photometric data, and the frequency and energy distribution of stellar flares on different types of stars are studied statistically. However, the completeness and observational bias of detected flare events from different missions (e.g. Kepler and TESS) vary a lot. We use a unified data processing and detection method for flares events based on the light curve from Kepler and TESS. Then we perform injection and recovery tests in the original light curve of each star for each flare event to correct the completeness and energy of flares. Three samples of flaring stars are selected from Kepler and TESS, with rotating periods from 1 to $\sim$ 5 days. Adopting a hot-blackbody assumption, our results show that the cumulative flare frequency distributions (FFDs) of the same stars in Kepler and TESS bands tend to be consistent after correction, revealing a more natural flaring frequency and energy distribution. Our results also extend the low-energy limit in cumulative FFD fitting to $10^{31.5-33}$ erg on different types of stars. For solar-type stars, the average power-law index of cumulative FFD ($α_{\rm cum}$) is $-0.84$, which indicates that low-energy flares contribute less to the total flare energy. With a piecewise correlation between $α_{\rm cum}$ and $T_{\rm eff}$, $α_{\rm cum}$ first rises with $T_{\rm eff}$ from M2 to K1 stars, then slightly decreases for stars hotter than K1.

preprint2022arXiv

Detecting and Monitoring Tidal Dissipation of Hot Jupiters in the Era of SiTian

Transit Timing Variation (TTV) of hot Jupiters provides direct observational evidence of planet tidal dissipation. Detecting tidal dissipation through TTV needs high precision transit timings and long timing baselines. In this work, we predict and discuss the potential scientific contribution of SiTian Survey in detecting and analyzing exoplanet TTV. We develop a tidal dissipation detection pipeline for SiTian Survey that aims at time-domain astronomy with 72 1-meter optical telescopes. The pipeline includes the modules of light curve deblending, transit timing obtaining, and TTV modeling. SiTian is capable to detect more than 25,000 exoplanets among which we expect $\sim$50 sources showing evidence of tidal dissipation. We present detection and analysis of tidal dissipating targets, based on simulated SiTian light curves of XO-3b and WASP-161b. The transit light curve modeling gives consistent results within 1$σ$ to input values of simulated light curves. Also, the parameter uncertainties predicted by Monte-Carlo Markov Chain are consistent with the distribution obtained from simulating and modeling the light curve 1000 times. The timing precision of SiTian observations is $\sim$ 0.5 minutes with one transit visit. We show that differences between TTV origins, e.g., tidal dissipation, apsidal precession, multiple planets, would be significant, considering the timing precision and baseline. The detection rate of tidal dissipating hot Jupiters would answer a crucial question of whether the planet migrates at an early formation stage or random stages due to perturbations, e.g., planet scattering, secular interaction. SiTian identified targets would be constructive given that the sample would extend tenfold.

preprint2022arXiv

Revisiting the Full Sets of Orbital Parameters for the XO-3 System: No evidence for Temporal Variation of the Spin-Orbit Angle

We present 12 new transit light curves and 16 new out-of-transit radial velocity measurements for the XO-3 system. By modelling our newly collected measurements together with archival photometric and Doppler velocimetric data, we confirmed the unusual configuration of the XO-3 system, which contains a massive planet ($M_P=11.92^{+0.59}_{-0.63} M_J$) on a relatively eccentric ($e=0.2853^{+0.0027}_{-0.0026}$) and short-period ($3.19152 \pm 0.00145\,$day) orbit around a massive star ($M_*=1.219^{+0.090}_{-0.095} M_{\odot}$). Furthermore, we find no strong evidence for a temporal change of either $V\sin i_{*}$ (and by extension, the stellar spin vector of XO-3), or the transit profile (and thus orbital angular momentum vector of XO-3b). We conclude that the discrepancy in previous Rossiter-McLaughlin measurements ($70.0^{\circ} \pm 15.0^{\circ}$ (Hebrard et al. 2008); $37.3^{\circ} \pm 3.7^{\circ}$ (Winn et al. 2009); $37.3^{\circ} \pm 3.0^{\circ}$ (Hirano et al. 2011)) may have stemmed from systematic noise sources.

preprint2022arXiv

Twinkle -- a small satellite spectroscopy mission for the next phase of exoplanet science

With a focus on off-the-shelf components, Twinkle is the first in a series of cost competitive small satellites managed and financed by Blue Skies Space Ltd. The satellite is based on a high-heritage Airbus platform that will carry a 0.45 m telescope and a spectrometer which will provide simultaneous wavelength coverage from 0.5-4.5 $\rm{μm}$. The spacecraft prime is Airbus Stevenage while the telescope is being developed by Airbus Toulouse and the spectrometer by ABB Canada. Scheduled to begin scientific operations in 2025, Twinkle will sit in a thermally-stable, sun-synchronous, low-Earth orbit. The mission has a designed operation lifetime of at least seven years and, during the first three years of operation, will conduct two large-scale survey programmes: one focused on Solar System objects and the other dedicated to extrasolar targets. Here we present an overview of the architecture of the mission, refinements in the design approach, and some of the key science themes of the extrasolar survey.

preprint2020arXiv

Occurrence and Architecture of Kepler Planetary Systems as Functions of Stellar Mass and Effective Temperature

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 $Δ$AIC> 30, though the trend in planet multiplicity is somewhat tentative with $Δ$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.

preprint2016arXiv

Analytical formulation of lunar cratering asymmetries

We formulate the lunar cratering distribution and verify the cratering asymmetries generated by the main-belt asteroids (MBAs) as well as the near-Earth objects (NEOs). Based on a planar model that excludes the terrestrial and lunar gravitations on the impactors and assuming the impactor encounter speed with Earth $v_{\rm{enc}}$ is higher than the lunar orbital speed $v_{\rm{M}}$, we rigorously integrated the lunar cratering distribution, and derived its approximation to the first order of $v_{\rm{M}}/v_{\rm{enc}}$. Numerical simulations of lunar bombardment by the MBAs during the late heavy bombardment were performed with an Earth-Moon distance $a_{\rm{M}}$ = 20--60 Earth radii in five cases. The analytical model directly proves the existence of a leading/trailing asymmetry and the absence of near/far asymmetry. The approximate form of the leading/trailing asymmetry is $(1 + A_1 \cosβ)$, which decreases as the apex distance $β$ increases. The numerical simulations show evidence of a pole/equator asymmetry as well as the leading/trailing asymmetry, and the former is empirically described as $(1 + A_2 \cos2φ)$, which decreases as the latitude modulus $|φ|$ increases. The amplitudes $A_{1,2}$ are reliable measurements of asymmetries. Our analysis explicitly indicates the quantitative relations between cratering distribution and bombardment conditions (impactor properties and the lunar orbital status) like $A_1 \propto v_{\rm{M}}/v_{\rm{enc}}$, resulting in a method for reproducing the bombardment conditions through measuring the asymmetry. Mutual confirmation between analytical model and numerical simulations is found in terms of the cratering distribution and its variation with $a_{\rm{M}}$. Estimates of $A_1$ for crater density distributions generated by the MBAs and the NEOs are 0.101--0.159 and 0.117, respectively.

preprint2016arXiv

Detecting planet pairs in mean motion resonances via astrometry method

GAIA leads us to step into a new era with a high astrometry precision of 10 uas. Under such a precision, astrometry will play important roles in detecting and characterizing exoplanets. Specially, we can identify planet pairs in mean motion resonances(MMRs) via astrometry, which constrains the formation and evolution of planetary systems. In accordance with observations, we consider two Jupiters or two super-Earths systems in 1:2, 2:3 and 3:4 MMRs. Our simulations show the false alarm probabilities(FAPs) of a third planet are extremely small while the real two planets can be good fitted with signal-to-noise ratio(SNR)> 3. The probability of reconstructing a resonant system is related with the eccentricities and resonance intensity. Generally, when SNR >= 10, if eccentricities of both planets are larger than 0.01 and the resonance is quite strong, the probabilities to reconstruct the planet pair in MMRs >= 80%. Jupiter pairs in MMRs are reconstructed more easily than super-Earth pairs with similar SNR when we consider the dynamical stability. FAPs are also calculated when we detect planet pairs in or near MMRs. FAPs for 1:2 MMR are largest, i.e., FAPs > 15% when SNR <= 10. Extrapolating from the Kepler planet pairs near MMRs and assuming SNR to be 3, we will discover and reconstruct a few tens of Jupiter pairs and hundreds of super-Earth pairs in 2:3 and 1:2 MMRs within 30 pc. We also compare the differences between even and uneven data cadence and find that planets are better measured with more uniform phase coverage.

preprint2016arXiv

Exoplanet Orbital Eccentricities Derived From LAMOST-Kepler Analysis

The nearly circular (mean eccentricity <e>~0.06) and coplanar (mean mutual inclination <i>~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 (<e> ~ 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 <e>~0.3, while the multiples are on nearly circular (<e> = 0.04^{+0.03}_{-0.04}) and coplanar (<i> = 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 <e>~(1-2)x<i> 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.

preprint2016arXiv

Global Instability of Exo-Moon System Triggered by Photo-Evaporation

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.

preprint2016arXiv

Stellar Flares in the CSTAR Field: Results from the 2008 Data Set

The Chinese Small Telescope ARray (CSTAR) is the first Chinese astronomical instrument placed in Antarctica. It is a group of four identical, fully automatic $14.5\,\rm{cm}$ telescopes, with an field of view (FOV) of $20\,\rm{deg^2}$ centered on the South Celestial Pole. Placed at Antarctic Dome A, CSTAR is designed to provide high-cadence photometry for site monitoring and variable sources detection. During the 2008 observing season, CSTAR has taken high-precision photometric data for 18,145 stars around the South Celestial Pole. At $i\,=\,7.5$ and $12$, the photometric precision reaches $\sim 8$ mmag and $\sim 30$ mmag with a cadence of 20s or 30s, respectively. Using robust detection method, we have found 15 stellar flares on 13 sources, including two classified variables. We have also found a linear relation between the decay times and the total durations of the stellar flares. The details of all detected flares along with their stellar properties are presented in this work.

preprint2015arXiv

Eclipsing Binaries From the CSTAR Project at Dome A, Antarctica

The Chinese Small Telescope ARray (CSTAR) has observed an area around the Celestial South Pole at Dome A since 2008. About $20,000$ light curves in the i band were obtained lasting from March to July, 2008. The photometric precision achieves about 4 mmag at i = 7.5 and 20 mmag at i = 12 within a 30 s exposure time. These light curves are analyzed using Lomb--Scargle, Phase Dispersion Minimization, and Box Least Squares methods to search for periodic signals. False positives may appear as a variable signature caused by contaminating stars and the observation mode of CSTAR. Therefore the period and position of each variable candidate are checked to eliminate false positives. Eclipsing binaries are removed by visual inspection, frequency spectrum analysis and locally linear embedding technique. We identify 53 eclipsing binaries in the field of view of CSTAR, containing 24 detached binaries, 8 semi-detached binaries, 18 contact binaries, and 3 ellipsoidal variables. To derive the parameters of these binaries, we use the Eclipsing Binaries via Artificial Intelligence (EBAI) method. The primary and the secondary eclipse timing variations (ETVs) for semi-detached and contact systems are analyzed. Correlated primary and secondary ETVs confirmed by false alarm tests may indicate an unseen perturbing companion. Through ETV analysis, we identify two triple systems (CSTAR J084612.64-883342.9 and CSTAR J220502.55-895206.7). The orbital parameters of the third body in CSTAR J220502.55-895206.7 are derived using a simple dynamical model.

preprint2015arXiv

Eight Planets in Four Multi-planet Systems via Transit Timing Variations in 1350 Days

Analysis of the transit timing variations (TTVs) of candidate pairs near mean-motion resonances (MMRs) is an effective method to confirm planets. Hitherto, 68 planets in 34 multi-planet systems have been confirmed via TTVs. We analyze the TTVs of all candidates from the most recent Kepler data with a time span as long as about 1350 days (Q0-Q15). The anti-correlations of TTV signals and the mass upper limits of candidate pairs in the same system are calculated, using an improved method suitable for long-period TTVs. If the false alarm probability (FAP) of a candidate pair is less than 10^{-3} and the mass upper limit for each candidate is less than 13 M_J, we confirm them as planets in the same system. Finally, 8 planets in 4 multi-planet systems are confirmed via analysis of their TTVs. All of the 4 planet pairs are near first-order MMRs, including KOI-2672 near 2:1 MMR, KOI-1236, KOI-1563 and KOI-2038 near 3:2 MMR. Four planets have relatively long orbital periods (> 35 day). KOI-2672.01 has a period of 88.51658 days and a fit mass of 17 M_\oplus. It is the longest-period planet confirmed, to date, near first-order MMR via TTVs.

preprint2015arXiv

Photometric Variability in the CSTAR Field: Results From the 2008 Data Set

The Chinese Small Telescope ARray (CSTAR) is the first telescope facility built at Dome A, Antarctica. During the 2008 observing season, the installation provided long-baseline and high-cadence photometric observations in the i-band for 18,145 targets within 20 deg2 CSTAR field around the South Celestial Pole for the purpose of monitoring the astronomical observing quality of Dome A and detecting various types of photometric variability. Using sensitive and robust detection methods, we discover 274 potential variables from this data set, 83 of which are new discoveries. We characterize most of them, providing the periods, amplitudes and classes of variability. The catalog of all these variables is presented along with the discussion of their statistical properties.

preprint2014arXiv

A search for eclipsing binaries that host discs

We search for systems hosting eclipsing discs using a complete sample of eclipsing binaries (EBs); those previously identified in the third phase of the Optical Gravitational Lensing Experiment (OGLE-III). Within a subsample of 2,823 high-cadence, high-photometric precision and large eclipsing depth detached EBs previously identified in the Large Magellanic Cloud (LMC), we find that the skewness and kurtosis of the light curves magnitude distribution within the primary eclipse can distinguish EBs hosting a disc from those without. Two systems with previously identified eclipsing discs (OGLE-LMC-ECL-11893 and OGLE-LMC-ECL-17782) are identified with near zero skewness ($|S|<0.5$) and positive kurtosis. No additional eclipsing disc systems were found in the OGLE-III LMC, Small Magellanic Cloud (SMC) or Galactic Disc (GD) EB light curves. We estimate that the fraction of detached near main-sequence LMC EBs (which have a primary with an $I$-band magnitude brighter than $\simeq 19\,\rm{mag}$) that host a disc is approximately 1/1000. As circumstellar disc lifetimes are short, we expected to primarily find eclipsing discs around young stars. In addition, as there is more room for a disc in a widely separated binary and because a disk close to a luminous star would be above the dust sublimation temperature, we expected to primarily find eclipsing discs in long period binaries. However, OGLE-LMC-ECL-17782 is a 13.3 day period B star system with a transient and hot ($\sim 6000\,\rm{K}$, $\sim 0.1\,\rm{AU}$ radius) disc and Scott et al. (in prep.) estimate an age of $150\,\rm{Myr}$ for OGLE-LMC-ECL-11893. Both discs are unexpected in the EB sample and impel explanation.

preprint2014arXiv

Planetary Transit Candidates in the CSTAR Field: Analysis of the 2008 Data

The Chinese Small Telescope ARray (CSTAR) is a group of four identical, fully automated, static 14.5 cm telescopes. CSTAR is located at Dome A, Antarctica and covers 20 square degree of sky around the South Celestial Pole. The installation is designed to provide high-cadence photometry for the purpose of monitoring the quality of the astronomical observing conditions at Dome A and detecting transiting exoplanets. CSTAR has been operational since 2008, and has taken a rich and high-precision photometric data set of 10,690 stars. In the first observing season, we obtained 291,911 qualified science frames with 20-second integrations in the i-band. Photometric precision reaches about 4 mmag at 20-second cadence at i=7.5, and is about 20 mmag at i=12. Using robust detection methods, ten promising exoplanet candidates were found. Four of these were found to be giants using spectroscopic follow-up. All of these transit candidates are presented here along with the discussion of their detailed properties as well as the follow-up observations.

preprint2014arXiv

Transits of Planets with Small Intervals in Circumbinary Systems

Transit times around single stars can be described well by a linear ephemeris. However, transit times in circumbinary systems are influenced both by the gravitational perturbations and the orbital phase variations of the central binary star. Adopting a coplanar analog of Kepler-16 as an example, we find that circumbinary planets can transit the same star more than once during a single planetary orbit, a phenomenon we call "tight transits". In certain geometric, the projected orbital velocity of the planet and the secondary star can approach zero and change sign, resulting in very long transits and/or 2-3 transits during a single binary orbit. Whether tight transits are possible for a particular system depends primarily on the binary mass ratio and the orbital architecture of both the binary and the planet. We derive a time-dependent criterion to judge when tight transits are possible for any circumbinary system. These results are verified with full dynamical integrations that also reveal other tight transit characteristics, i.e., the transit durations and the intervals between tight transits. For the seven currently known circumbinary systems, we estimate these critical parameters both analytically and numerically. Due to the mutual inclination between the planet and the binary, tight transits can only occur across the less massive star B in Kepler-16, -34, -35, and -47 (for both planets). The long-term average frequency of tight transits (compared to typical transits) for Kepler-16, -34, and -35 are estimated to be several percent. Using full numerical integrations, the next tight transit for each system is predicted and the soonest example appears to be Kepler-47b and -47c, which are likely to have tight transits before 2025. These unique and valuable events often deserve special observational scrutiny.

preprint2013arXiv

A Mechanism of Exciting Planetary Inclination and Eccentricity through a Residual Gas Disk

Accordling to the theory of Kozai resonance, the initial mutual inclination between a small body and a massive planet in an outer circular orbit is as high as $\sim39.2^{\circ}$ for pumping the eccentricity of the inner small body. Here we show that, with the presence of a residual gas disk outside two planetary orbits, the inclination can be reduced as low as a few degrees. The presence of disk changes the nodal precession rates and directions of the planet orbits. At the place where the two planets achieve the same nodal processing rate, vertical secular resonance would occur so that mutual inclination of the two planets will be excited, which might trigger the Kozai resonance between the two planets further. However, in order to pump an inner Jupiter-like planet, the conditions required for the disk and the outer planet are relatively strict. We develop a set of evolution equations, which can fit the N-body simulation quite well but be integrated within a much shorter time. By scanning the parameter spaces using the evolution equations, we find that, a massive planet ($10M_J$) at 30AU with $6^o$ inclined to a massive disk ($50M_J$) can finally enter the Kozai resonance with an inner Jupiter around the snowline. And a $20^{\circ}$ inclination of the outer planet is required for flipping the inner one to a retrograde orbit. In multiple planet systems, the mechanism can happen between two nonadjacent planets, or inspire a chain reaction among more than two planets. This mechanism could be the source of the observed giant planets in moderate eccentric and inclined orbits, or hot-Jupiters in close-in, retrograde orbits after tidal damping.

preprint2013arXiv

Configurations of Bounded and Free-floating Planets in Very Young Open Clusters

Open clusters(OCs) are usually young and suitable for studying the formation and evolution of planetary systems. Hitherto, only four planets have been found with radial velocity measurements in OCs. Meanwhile, a lot of free-floating planets(FFPs) have been detected. We utilize N-body simulations to investigate the evolution and final configurations of multi-planetary systems in very young open clusters with an age <10 Myr. After an evolution of 10 Myr, 61%-72% of the planets remain bounded and more than 55% of the planetary systems will maintain their initial orbital configurations. For systems with one planet ejected, more than 25% of them have the surviving planets in misaligned orbits. In the clusters, the fraction of planetary systems with misalignment is > 6%, and only 1% have planets in retrograde orbits. We also obtain a positive correlation between the survival planet number and the distance from the cluster center $r$: planetary systems with a larger $r$ tend to be more stable. Moreover, stars with a mass >2.5 $M_\odot$ are likely unstable and lose their planets. These results are roughly consistent with current observations. Planetary systems in binaries are less stable and we achieve a rough criterion: most of the binary systems with $a_{\rm b}(1-e_{\rm b}^2)$>100 AU can keep all the initial planets survived. Finally, 80% of the FFPs are ejected out of the clusters, while the rest (~20%) still stay in host clusters and most of them are concentrated in the center (<2 pc).

preprint2013arXiv

Forever alone? Testing single eccentric planetary systems for multiple companions

Determining the orbital eccentricity of an extrasolar planet is critically important for understanding the system's dynamical environment and history. However, eccentricity is often poorly determined or entirely mischaracterized due to poor observational sampling, low signal-to-noise, and/or degeneracies with other planetary signals. Some systems previously thought to contain a single, moderate-eccentricity planet have been shown, after further monitoring, to host two planets on nearly-circular orbits. We investigate published apparent single-planet systems to see if the available data can be better fit by two lower-eccentricity planets. We identify nine promising candidate systems and perform detailed dynamical tests to confirm the stability of the potential new multiple-planet systems. Finally, we compare the expected orbits of the single- and double-planet scenarios to better inform future observations of these interesting systems.

preprint2013arXiv

Gap formation in a self-gravitating disk and the associated migration of the embedded giant planet

We present the results of our recent study on the interactions between a giant planet and a self-gravitating gas disk. We investigate how the disk's self-gravity affects the gap formation process and the migration of the giant planet. Two series of 1-D and 2-D hydrodynamic simulations are performed. We select several surface densities and focus on the gravitationally stable region. To obtain more reliable gravity torques exerted on the planet, a refined treatment of disk's gravity is adopted in the vicinity of the planet. Our results indicate that the net effect of the disk's self-gravity on the gap formation process depends on the surface density of the disk. We notice that there are two critical values, Σ_I and Σ_II. When the surface density of the disk is lower than the first one, Σ_0 < Σ_I, the effect of self-gravity suppresses the formation of a gap. When Σ_0 > Σ_I, the self-gravity of the gas tends to benefit the gap formation process and enlarge the width/depth of the gap. According to our 1-D and 2-D simulations, we estimate the first critical surface density Σ_I \approx 0.8MMSN. This effect increases until the surface density reaches the second critical value Σ_II. When Σ_0 > Σ_II, the gravitational turbulence in the disk becomes dominant and the gap formation process is suppressed again. Our 2-D simulations show that this critical surface density is around 3.5MMSN. We also study the associated orbital evolution of a giant planet. Under the effect of the disk's self-gravity, the migration rate of the giant planet increases when the disk is dominated by gravitational turbulence. We show that the migration timescale associates with the effective viscosity and can be up to 10^4 yr.

preprint2013arXiv

The effect of planet-planet scattering on the survival of exomoons

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.

preprint2013arXiv

Where to Find Habitable "Earths" in Circumbinary Systems

Hitherto, six P-type planets are found around five binary systems, i.e. Kepler-16 b, 34 b, 35 b, 38 b, 47 b, c, which are all Neptune or Jupiter-like planets. The stability of planets and the habitable zones are influenced by the gravitational and radiative perturbations of binary companions. In this Letter, we check the stability of an additional habitable Earth-mass planet in each system. Based on our simulations in 10 Myr, a habitable "Earth" is hardly stable in Kepler-16 while a stable "Earth" in Kepler-47 close to the boundaries of the habitable zone is possible. In contrast, Kepler-34, 35 and 38 seem to have high probabilities of being able to tolerate a stable "Earth" in their habitable zones. The affects of transit time variations are quite small due to the small mass of an undetected "Earth", except that of Kepler-16 b. With a time precision of 0.001 day(~ 88 s), an "Earth" in the corotational resonance with Kepler-16 b can be detected in 3 years, while habitable "Earths" in Kepler-34 and 38 systems can be detected in 10 years. Habitable "Earths" in Kepler-35 and 47 are not likely to be detected in 10 years under this precision.

preprint2012arXiv

Planetesimal Dynamics in Inclined Binary Systems: The Role of Gas-Disk Gravity

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.

preprint2012arXiv

Predicting the Configuration of Planetary System: KOI-152 Observed by Kepler

The recent Kepler discovery of KOI-152 reveals a system of three hot super-Earth candidates that are in or near a 4:2:1 mean motion resonance. It is unlikely that they formed in situ, the planets probably underwent orbital migration during the formation and evolution process. The small semimajor axes of the three planets suggest that migration stopped at the inner edge of the primordial gas disk. In this paper we focus on the influence of migration halting mechanisms, including migration "dead zones", and inner truncation by the stellar magnetic field. We show that the stellar accretion rate, stellar magnetic field and the speed of migration in the proto-planetary disk are the main factors affecting the final configuration of KOI-152. Our simulations suggest that three planets may be around a star with low star accretion rate or with high magnetic field. On the other hand, slow type I migration, which decreases to one tenth of the linear analysis results, favors forming the configuration of KOI-152. Under such formation scenario, the planets in the system are not massive enough to open gaps in the gas disk. The upper limit of the planetary masses are estimated to be about $15,~19$, and $24 M_\oplus$, respectively. Our results are also indicative of the near Laplacian configurations that are quite common in planetary systems.

preprint2012arXiv

Terrestrial Planets Formation around Circumbinary Habitable Zone: Inward Migration in the Planetesimal Swarm

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.

preprint2012arXiv

The silicate model and carbon rich model of CoRoT-7b, Kepler-9d and Kepler-10b

Possible bulk compositions of the super-Earth exoplanets, CoRoT-7b, Kepler-9d, and Kepler-10b are investigated by applying a commonly used silicate and a non-standard carbon model. Their internal structures are deduced using the suitable equation of state of the materials. The degeneracy problems of their compositions can be partly overcome, based on the fact that all three planets are extremely close to their host stars. By analyzing the numerical results, we conclude: 1) The iron core of CoRoT-7b is not more than 27% of its total mass within 1 $σ$ mass-radius error bars, so an Earth-like composition is less likely, but its carbon rich model can be compatible with an Earth-like core/mantle mass fraction; 2) Kepler-10b is more likely with a Mercury-like composition, its old age implies that its high iron content may be a result of strong solar wind or giant impact; 3) the transiting-only super-Earth Kepler-9d is also discussed. Combining its possible composition with the formation theory, we can place some constraints on its mass and bulk composition.

preprint2011arXiv

Modeling Planetary System Formation with N-Body Simulations: Role of Gas Disk and Statistics Comparing to Observations

During the late stage of planet formation when Mars-size cores appear, interactions among planetary cores can excite their orbital eccentricities, speed their merges and thus sculpture the final architecture of planet systems. This series of work contributes to the final assembling of planet systems with N-body simulations, including the type I and II migration of planets, gas accretion of massive cores in a viscous disk. In this paper, the standard formulations of type I and II migrations are adopted to investigate the formation of planet systems around solar mass stars. Statistics on the final distributions of planetary masses, semi-major axes and eccentricities are derived, which are comparable to those of the observed systems. Our simulations predict some orbital signatures of planet systems around solar mass stars: 36% of the survival planets are giant planets (Mp>10Me). Most of the massive giant planets (Mp>30Me) locate at 1-10AU. Terrestrial planets distribute more or less evenly at <1-2 AU. Planets in inner orbits (<1 AU) may accumulate at the inner edges of either the protostellar disk (3-5 days) or its MRI dead zone (30-50 days). There is a planet desert in the mass-eccecntricity diagram, i.e., lack of planets with masses 0.005 - 0.08 MJ in highly eccentric orbits (e > 0.3 - 0.4). The average eccentricity (~ 0.15) of the giant planets (Mp>10Me) are bigger than that (~ 0.05) of the terrestrial planets (Mp< 10Me). A planet system with more planets tends to have smaller planet masses and orbital eccentricities on average.

preprint2011arXiv

Planet Formation In Highly Inclined Binary Systems I. Planetesimals Jump Inwards And Pile Up

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 <i_B<150^\circ$). We then find that gas drag plays a crucial role, and the results can be generally divided into two categories, i.e., the Kozai-on regime and the Kozai-off regime, depending on the specific value of $i_B$. For both regimes, a robust outcome over a wide range of parameters is that, planetesimals migrate/jump inwards and pile up, leading to a severely truncated and dense planetesimal disk around the primary. In this compact and dense disk, collision rates are high but relative velocities are low, providing conditions which are favorable for planetesimal growth, and potentially allow for the subsequent formation of planets.

preprint2010arXiv

Forming Habitable Planets around Dwarf Stars: Application to OGLE-06-109L

Dwarf stars are believed to have small protostar disk where planets may grow up. During the planet formation stage, embryos undergoing type I migration are expected to be stalled at inner edge of magnetic inactive disk ($a_{\rm crit} \sim 0.2-0.3 ~$AU). This mechanism makes the location around $a_{\rm crit}$ a sweet spot of forming planets. Especially, $a_{\rm crit}$ of dwarf stars with masses $\sim 0.5 M_\odot$ is roughly inside the habitable zone of the system. In this paper we study the formation of habitable planets due to this mechanism with a model system OGLE-06-109L. It has a $0.51 M_\odot$ dwarf star with two giant planets in 2.3 and 4.6 AU observed by microlensing. We model the embryos undergoing type I migration in the gas disk with a constant disk accretion rate ($\dot M$). Giant planets in outside orbits affect the formation of habitable planets through secular perturbations at the early stage and secular resonance at the later stage. We find that the existence and the masses of the habitable planets in OGLE-06-10L system depend on both $\dot M$ and the speed of type I migration. If planets formed earlier so that $\dot M$ is larger ($\sim 10^{-7} M_\odot $ yr$^{-1}$), terrestrial planet can not be survived unless the type I migration rate is an order of magnitude less. If planets formed later so that $\dot M$ is smaller ($\sim 10^{-8} M_\odot $ yr$^{-1}$), single and high mass terrestrial planets with high water contents ($\sim 5% $) will be formed by inward migration of outer planet cores. A slower speed migration will result in several planets by collisions of embryos, thus their water contents are low ($\sim 2%$). Mean motion resonances or apsidal resonances among planets may be observed if multiple planets survived in the inner system.

preprint2010arXiv

From Dust To Planetesimal: The Snowball Phase ?

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.

preprint2010arXiv

On the orbital evolution of a giant planet pair embedded in a gaseous disk. I: Jupiter-Saturn configuration

We carry out a series of high resolution ($1024\times 1024$) hydrodynamical simulations to investigate the orbital evolution of Jupiter and Saturn embedded in a gaseous protostellar disk. Our work extends the results in the classical papers of Masset & Snellgrove (2001) and Morbidelli & Crida (2007) by exploring various surface density profiles ($σ$), where $σ\propto r^{-α}$. The stability of the mean motion resonances(MMRs) caused by the convergent migration of the two planets is studied as well. Our results show that:(1) The gap formation process of Saturn is greatly delayed by the tidal perturbation of Jupiter. These perturbations cause inward or outward runaway migration of Saturn, depending on the density profiles on the disk. (2) The convergent migration rate increases as $α$ increases and the type of MMRs depends on $α$ as well. When $0<α<1$, the convergent migration speed of Jupiter and Saturn is relatively slow, thus they are trapped into 2:1 MMR. When $α>4/3$, Saturn passes through the $2:1$ MMR with Jupiter and is captured into the $3:2$ MMR. (3) The $3:2$ MMR turns out to be unstable when the eccentricity of Saturn ($e_s$) increases too high. The critical value above which instability will set in is $e_s \sim 0.15$. We also observe that the two planets are trapped into $2:1$ MMR after the break of $3:2$ MMR. This process may provide useful information for the formation of orbital configuration between Jupiter and Saturn in the Solar System.

preprint2010arXiv

On the orbital evolution of a giant planet pair embedded in a gaseous disk. II. A Saturn-Jupiter configuration

We carry out a series of high-resolution (1024 X 1024) hydrodynamic simulations to investigate the orbital evolution of a Saturn-Jupiter pair embedded in a gaseous disk. This work extends the results of our previous work by exploring a different orbital configuration---Jupiter lies outside Saturn (q<1, where q= M_i/M_o is the mass ratio of the inner planet and the outer one). We focus on the effects of different initial separations (d) between the two planets and the various surface density profiles of the disk, where σ\propto r^{-α}. We also compare the results of different orbital configurations of the planet pair. Our results show that: (1) when the initial separation is relatively large(d>d_{iLr}, where d_{iLr} is the distance between Jupiter and its first inner Lindblad resonance), the two planets undergo divergent migration. However, the inward migration of Saturn could be halted when Jupiter compresses the inner disk in which Saturn is embedded. (2) Convergent migration occurs when the initial separation is smaller (d<d_{iLr}) and the density slope of the disk is nearly flat (α<1/2). Saturn is then forced by Jupiter to migrate inward when the two planets are trapped into mean motion resonances (MMRs), and Saturn may get very close to the central star. (3) In the case of q<1, the eccentricity of Saturn could be excited to a very high value (e_{S}~0.4-0.5) by the MMRs and the system could maintain stability. These results explain the formation of MMRs in the exoplanet systems where the outer planet is more massive than the inner one. It also helps us to understand the origin of the "hot Jupiter/Saturn" undergoing high eccentric orbit.

preprint2010arXiv

Planetesimal Accretion in Binary Systems: Could Planets Form Around Alpha Centauri B ?

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.

preprint2009arXiv

Dynamics and Eccentricity Formation of Planets in OGLE-06-109L System

Recent observation of microlensing technique reveals two giant planets at 2.3 AU and 4.6 AU around the star OGLE-06-109L. The eccentricity of the outer planet (ec) is estimated to be 0.11(+0.17,-0.04), comparable to that of Saturn (0.01-0.09). The similarities between the OGLE-06-109L system and the solar system indicate that they may have passed through similar histories during their formation stage. In this paper we investigate the dynamics and formation of the orbital architecture in the OGLE-06-109L system. For the present two planets with their nominal locations, the secular motions are stable as long as their eccentricities (eb; ec) fulfill eb^2 + ec^2 < = 0.3^2. Earth-size bodies might be formed and are stable in the habitable zone (0.25AU-0.36AU) of the system. Three possible scenarios may be accounted for formation of eb and ec: (i) convergent migration of two planets and the 3:1 MMR trapping; (ii) planetary scattering; (iii) divergent migration and the 3:1 MMR crossing. As we showed that the probability for the two giant planets in 3:1 MMR is low (~3%), scenario (i) is less likely. According to models (ii) and (iii), the final eccentricity of inner planet (eb) may oscillate between [0-0.06], comparable to that of Jupiter (0.03-0.06). An inspection of eb, ec secular motion may be helpful to understand which model is really responsible for the eccentricity formation.

preprint2009arXiv

Formation and tidal evolution of hot super-Earths in multiple planetary systems

Hot super-Earths are exoplanets with masses < 10 Earth masses and orbital periods < 20 days. Around 8 hot super-Earths have been discovered in the neighborhood of solar system. In this lecture, we review the mechanisms for the formation of hot super-Earths, dynamical effects that play important roles in sculpting the architecture of the multiple planetary systems. Two example systems (HD 40307 and GJ 436) are presented to show the formation and evolution of hot super-Earths or Neptunes.

preprint2009arXiv

Planetesimal Accretion in Binary Systems: Role of the Companion's Orbital Inclination

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.

preprint2008arXiv

Eccentricity modulation of a close-in planet by a companion - application to GJ 436 system

GJ 436b is a Neptune-size planet with 23.2 Earth masses in an elliptical orbit of period 2.64 days and eccentricity 0.16. With a typical tidal dissipation factor (Q' ~ 10^6) as that of a giant planet with convective envelope, its orbital circularization timescale under internal tidal dissipation is around 1 Gyr, at least two times less than the stellar age (>3 Gyr). A plausible mechanism is that the eccentricity of GJ 436b is modulated by a planetary companion due to their mutual perturbation. Here we investigate this possibility from the dynamical viewpoint. A general method is given to predict the possible locations of the dynamically coupled companions, including in nearby/distance non-resonant or mean motion resonance orbits with the first planet. Applying the method to GJ 436 system, we find it is very unlikely that the eccentricity of GJ 436b is maintained at the present location by a nearby/distance companion through secular perturbation or mean motion resonance. In fact, in all these simulated cases, GJ 436b will undergo eccentricity damp and orbital decay, leaving the present location within the stellar age. However, these results do not rule out the possible existence of planet companions in nearby/distance orbits, although they are not able to maintain the eccentricity of GJ 436b.

preprint2008arXiv

Migration and Final Location of Hot Super Earths in the Presence of Gas Giants

Based on the conventional sequential-accretion paradigm, we have proposed that, during the migration of first-born gas giants outside the orbits of planetary embryos, super Earth planets will form inside the 2:1 resonance location by sweeping of mean motion resonances (Zhou et al. 2005). In this paper, we study the subsequent evolution of a super Earth (m_1) under the effects of tidal dissipation and perturbation from a first-born gas giant (m_2) in an outside orbit. Secular perturbation and mean motion resonances (especially 2:1 and 5:2 resonances) between m_1 and m_2 excite the eccentricity of m_1, which causes the migration of m_1 and results in a hot super Earth. The calculated final location of the hot super Earth is independent of the tidal energy dissipation factor Q'. The study of migration history of a Hot Super Earth is useful to reveal its Q' value and to predict its final location in the presence of one or more hot gas giants. When this investigation is applied to the GJ876 system, it correctly reproduces the observed location of GJ876d around 0.02AU.