Source author record

D. N. C. Lin

D. N. C. Lin 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

34works
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

34 published item(s)

preprint2020arXiv

A Shining Death of Unequal Supermassive Black Hole Binaries

In the $Λ$CDM scenario, small galaxies merge to produce larger entities. Since supermassive black holes (SMBHs) are found in galaxies of all sizes, SMBH binaries (SMBHBs) are generally expected to form during the amalgamation of galaxies. It is unclear what fraction of these binaries could eventually merge, but a general consensus is that initially the orbital decay is mediated by the surrounding gas and stars. In this Letter, we show that in active galactic nulcei (AGNs) the radiation field also causes the orbits of the accreting SMBHs to shrink. The corresponding mechanism, known as the "Poynting-Robertson drag" (PR drag), takes effect on a well-defined timescale $CT_{\rm Sal}$, where $T_{\rm Sal}$ is the Salpeter timescale of the AGN, presumably coincide with the primary SMBH, and $C=4ξ^{-1}ε^{-1}q^{1/3}(1+q)^{2/3}(1-ε)$ is a constant determined by the radiative efficiency $ε$, the mass ratio $q$ of the two black holes, and a parameter $ξ$ characterizing the size of the circum-secondary accretion disk. We find that when $q\lesssim$a few$\times10^{-5}$, the PR drag is more efficient in shrinking the binary than many other mechanisms, such as dynamical friction and type-I migration. Our finding points to a possible new channel for the coalescence of unequal SMBHBs and the clearing of intermediate-massive black holes in AGNs.

preprint2020arXiv

Departure from the Exact Location of Mean Motion Resonances Induced by the Gas Disk in the Systems Observed by Kepler

The statistical results of transiting planets show that there are two peaks around 1.5 and 2.0 in the distribution of orbital period ratios. A large number of planet pairs are found near the exact location of mean motion resonances (MMRs). In this work, we find out that the depletion and structures of gas disk play crucial roles in driving planet pairs out of exact location of MMRs. Under such scenario, planet pairs are trapped into exact MMRs during orbital migration firstly and keep migrating in a same pace. The eccentricities can be excited. Due to the existence of gas disk, eccentricities can be damped leading to the change of orbital period. It will make planet pairs depart from the exact location of MMRs. With depletion timescales larger than 1 Myr, near MMRs configurations are formed easily. Planet pairs have higher possibilities to escape from MMRs with higher disk aspect ratio. Additionally, with weaker corotation torque, planet pairs can depart farther from exact location of MMRs. The final location of the innermost planets in systems are directly related to the transition radius from optically thick region to inner optically thin disk. While the transition radius is smaller than 0.2 AU at the late stage of star evolution process, the innermost planets can reach around 10 days. Our formation scenario is a possible mechanism to explain the formation of near MMRs configuration with the innermost planet farther than 0.1 AU.

preprint2016arXiv

An alternative model for the origin of gaps in circumstellar disks

Motivated by recent observational and numerical studies suggesting that collapsing protostellar cores may be replenished from the local environment, we explore the evolution of protostellar cores submerged in the external counter-rotating environment. These models predict the formation of counter-rotating disks with a deep gap in the gas surface density separating the inner disk (corotating with the star) and the outer counter-rotating disk. The properties of these gaps are compared to those of planet-bearing gaps that form in disks hosting giant planets. We employ numerical hydrodynamics simulations of collapsing cores that are replenished from the local counter-rotating environment, as well as numerical hydrodynamic simulations of isolated disks hosting giant planets, to derive the properties of the gaps that form in both cases. Our numerical simulations demonstrate that counter-rotating disks can form for a wide range of mass and angular momentum available in the local environment. The gap that separates both disks has a depletion factor smaller than 1%, can be located at a distance from ten to over a hundred AU from the star, and can propagate inward with velocity ranging from 1 AU/Myr to >100 AU/Myr. Unlike our previous conclusion, the gap can therefore be a long-lived phenomenon, comparable in some cases to the lifetime of the disk itself. For a proper choice of the planetary mass, the viscous α-parameter and the disk mass, the planet-bearing gaps and the gaps in counter-rotating disks may show a remarkable similarity in the gas density profile and depletion factor, which may complicate their observational differentiation.

preprint2015arXiv

Giant Impact: An Efficient Mechanism for the Devolatilization of Super-Earths

Mini-Neptunes and volatile-poor super-Earths coexist on adjacent orbits in proximity to host stars such as Kepler-36 and Kepler-11. Several post-formation processes have been proposed for explaining the origin of the compositional diversity: the mass loss via stellar XUV irradiation, degassing of accreted material, and in-situ accumulation of the disk gas. Close-in planets are also likely to experience giant impacts during the advanced stage of planet formation. This study examines the possibility of transforming volatile-rich super-Earths / mini-Neptunes into volatile-depleted super-Earths through giant impacts. We present the results of three-dimensional giant impact simulations in the accretionary and disruptive regimes. Target planets are modeled with a three-layered structure composed of an iron core, silicate mantle and hydrogen/helium envelope. In the disruptive case, the giant impact can remove most of the H/He atmosphere immediately and homogenize the refractory material in the planetary interior. In the accretionary case, the planet can retain more than half of the gaseous envelope, while a compositional gradient suppresses efficient heat transfer as its interior undergoes double-diffusive convection. After the giant impact, a hot and inflated planet cools and contracts slowly. The extended atmosphere enhances the mass loss via both a Parker wind induced by thermal pressure and hydrodynamic escape driven by the stellar XUV irradiation. As a result, the entire gaseous envelope is expected to be lost due to the combination of those processes in both cases. We propose that Kepler-36b may have been significantly devolatilized by giant impacts, while a substantial fraction of Kepler-36c's atmosphere may remain intact. Furthermore, the stochastic nature of giant impacts may account for the large dispersion in the mass--radius relationship of close-in super-Earths and mini-Neptunes.

preprint2015arXiv

On the IMF in a Triggered Star Formation Context

The origin of the stellar initial mass function (IMF) is a fundamental issue in the theory of star formation. It is generally fit with a composite power law. Some clues on the progenitors can be found in dense starless cores that have a core mass function (CMF) with a similar shape. In the low-mass end, these mass functions increase with mass, albeit the sample may be somewhat incomplete; in the high-mass end, the mass functions decrease with mass. There is an offset in the turn-over mass between the two mass distributions. The stellar mass for the IMF peak is lower than the corresponding core mass for the CMF peak in the Pipe Nebula by about a factor of three. Smaller offsets are found between the IMF and the CMFs in other nebulae. We suggest that the offset is likely induced during a starburst episode of global star formation which is triggered by the formation of a few O/B stars in the multi-phase media, which naturally emerged through the onset of thermal instability in the cloud-core formation process. We consider the scenario that the ignition of a few massive stars photoionizes the warm medium between the cores, increases the external pressure, reduces their Bonnor?Ebert mass, and triggers the collapse of some previously stable cores. We quantitatively reproduce the IMF in the low-mass end with the assumption of additional rotational fragmentation.

preprint2014arXiv

Constraining the Oblateness of Kepler Planets

We use Kepler short cadence light curves to constrain the oblateness of planet candidates in the Kepler sample. The transits of rapidly rotating planets that are deformed in shape will lead to distortions in the ingress and egress of their light curves. We report the first tentative detection of an oblate planet outside of the solar system, measuring an oblateness of $0.22 \pm 0.11$ for the 18 $M_J$ mass brown dwarf Kepler 39b (KOI-423.01). We also provide constraints on the oblateness of the planets (candidates) HAT-P-7b, KOI-686.01, and KOI-197.01 to be < 0.067, < 0.251, and < 0.186, respectively. Using the Q'-values from Jupiter and Saturn, we expect tidal synchronization for the spins of HAT-P-7b, KOI-686.01 and KOI-197.01, and for their rotational oblateness signatures to be undetectable in the current data. The potentially large oblateness of KOI-423.01 (Kepler 39b) suggests that the Q'-value of the brown dwarf needs to be two orders of magnitude larger than that of the solar system gas giants to avoid being tidally spun-down.

preprint2014arXiv

Embryo impacts and gas giant mergers II: Diversity of Hot Jupiters' internal structure

We consider the origin of compact, short-period, Jupiter-mass planets. We propose that their diverse structure is caused by giant impacts of embryos and super-Earths or mergers with other gas giants during the formation and evolution of these hot Jupiters. Through a series of numerical simulations, we show that typical head-on collisions generally lead to total coalescence of impinging gas giants. Although extremely energetic collisions can disintegrate the envelope of gas giants, these events seldom occur. During oblique and moderately energetic collisions, the merger products retain higher fraction of the colliders' cores than their envelopes. They can also deposit considerable amount of spin angular momentum to the gas giants and desynchronize their spins from their orbital mean motion. We find that the oblateness of gas giants can be used to infer the impact history. Subsequent dissipation of stellar tide inside the planets' envelope can lead to runaway inflation and potentially a substantial loss of gas through Roche-lobe overflow. The impact of super-Earths on parabolic orbits can also enlarge gas giant planets' envelope and elevates their tidal dissipation rate over $\sim $ 100 Myr time scale. Since giant impacts occur stochastically with a range of impactor sizes and energies, their diverse outcomes may account for the dispersion in the mass-radius relationship of hot Jupiters.

preprint2014arXiv

Planet Population Synthesis

With the increasing number of exoplanets discovered, statistical properties of the population as a whole become unique constraints on planet formation models provided a link between the description of the detailed processes playing a role in this formation and the observed population can be established. Planet population synthesis provides such a link. The approach allows to study how different physical models of individual processes (e.g., proto-planetary disc structure and evolution, planetesimal formation, gas accretion, migration, etc.) affect the overall properties of the population of emerging planets. By necessity, planet population synthesis relies on simplified descriptions of complex processes. These descriptions can be obtained from more detailed specialised simulations of these processes. The objective of this chapter is twofold: 1) provide an overview of the physics entering in the two main approaches to planet population synthesis and 2) present some of the results achieved as well as illustrate how it can be used to extract constraints on the models and to help interpret observations.

preprint2014arXiv

The effect of external environment on the evolution of protostellar disks

Using numerical hydrodynamics simulations we studied the gravitational collapse of pre-stellar cores of sub-solar mass embedded into a low-density external environment. Four models with different magnitude and direction of rotation of the external environment with respect to the central core were studied and compared with an isolated model. We found that the infall of matter from the external environment can significantly alter the disk properties as compared to those seen in the isolated model. Depending on the magnitude and direction of rotation of the external environment, a variety of disks can form including compact (<= 200 AU) ones shrinking in size due to infall of external matter with low angular momentum, as well as extended disks forming due to infall of external matter with high angular momentum. The former are usually stable against gravitational fragmentation, while the latter are prone to fragmentation and formation of stellar systems with sub-stellar/very-low-mass companions. In the case of counterrotating external environment, very compact (< 5 AU) and short-lived (<= a few * 10^5 yr) disks can form when infalling material has low angular momentum. The most interesting case is found for the infall of counterrotating external material with high angular momentum, leading to the formation of counterrotating inner and outer disks separated by a deep gap at a few tens AU. The gap migrates inward due to accretion of the inner disk onto the protostar, turns into a central hole, and finally disappears giving way to the outer strongly gravitationally unstable disk. This model may lead to the emergence of a transient stellar system with sub-stellar/very-low-mass components counterrotating with respect to that of the star.

preprint2013arXiv

Are the Kepler Near-Resonance Planet Pairs due to Tidal Dissipation?

The multiple-planet systems discovered by the Kepler mission show an excess of planet pairs with period ratios just wide of exact commensurability for first-order resonances like 2:1 and 3:2. In principle, these planet pairs could have both resonance angles associated with the resonance librating if the orbital eccentricities are sufficiently small, because the width of first-order resonances diverges in the limit of vanishingly small eccentricity. We consider a widely-held scenario in which pairs of planets were captured into first-order resonances by migration due to planet-disk interactions, and subsequently became detached from the resonances, due to tidal dissipation in the planets. In the context of this scenario, we find a constraint on the ratio of the planet's tidal dissipation function and Love number that implies that some of the Kepler planets are likely solid. However, tides are not strong enough to move many of the planet pairs to the observed separations, suggesting that additional dissipative processes are at play.

preprint2013arXiv

Internal Gravity Waves in Massive Stars: Angular Momentum Transport

We present numerical simulations of internal gravity waves (IGW) in a star with a convective core and extended radiative envelope. We report on amplitudes, spectra, dissipation and consequent angular momentum transport by such waves. We find that these waves are generated efficiently and transport angular momentum on short timescales over large distances. We show that, as in the Earth's atmosphere, IGW drive equatorial flows which change magnitude and direction on short timescales. These results have profound consequences for the observational inferences of massive stars, as well as their long term angular momentum evolution. We suggest IGW angular momentum transport may explain many observational mysteries, such as: the misalignment of hot Jupiters around hot stars, the Be class of stars, Ni enrichment anomalies in massive stars and the non-synchronous orbits of interacting binaries.

preprint2013arXiv

On the Coagulation and Size Distribution of Pressure Confined Cores

Observations of the Pipe Nebula have led to the discovery of dense starless cores. The mass of most cores is too small for their self gravity to hold them together. Instead, they are thought to be pressure confined. The observed dense cores' mass function (CMF) matches well with the initial mass function (IMF) of stars in young clusters. Similar CMF's are observed in other star forming regions such as the Aquila Nebula, albeit with some dispersion. The shape of these CMF provides important clues to the competing physical processes which lead to star formation and its feedback on the interstellar media. In this paper, we investigate the dynamical origin of the the mass function of starless cores which are confined by a warm, less dense medium. We consider the coagulation between the cold cores and their ablation due to Kelvin-Helmholtz instability induced by their relative motion through the warm medium. We are able to reproduce the observed CMF among the starless cores in the Pipe nebula. Our results indicate that in environment similar to the Pipe nebula: 1) before the onset of their gravitational collapse, the mass distribution of the progenitor cores is similar to that of the young stars, 2) the observed CMF is a robust consequence of dynamical equilibrium between the coagulation and ablation of cores, and 3) a break in the slope of the CMF is due to the enhancement of collisional cross section and suppression of ablation for cores with masses larger than the cores' Bonnor-Ebert mass.

preprint2013arXiv

On the Tidal Dissipation of Obliquity

We investigate tidal dissipation of obliquity in hot Jupiters. Assuming an initial random orientation of obliquity and parameters relevant to the observed population, the obliquity of hot Jupiters does not evolve to purely aligned systems. In fact, the obliquity evolves to either prograde, retrograde or 90^{o} orbits where the torque due to tidal perturbations vanishes. This distribution is incompatible with observations which show that hot jupiters around cool stars are generally aligned. This calls into question the viability of tidal dissipation as the mechanism for obliquity alignment of hot Jupiters around cool stars.

preprint2013arXiv

Toward a Deterministic Model of Planetary Formation VII: Eccentricity Distribution of Gas Giants

The ubiquity of planets and diversity of planetary systems reveal planet formation encompass many complex and competing processes. In this series of papers, we develop and upgrade a population synthesis model as a tool to identify the dominant physical effects and to calibrate the range of physical conditions. Recent planet searches leads to the discovery of many multiple-planet systems. Any theoretical models of their origins must take into account dynamical interaction between emerging protoplanets. Here, we introduce a prescription to approximate the close encounters between multiple planets. We apply this method to simulate the growth, migration, and dynamical interaction of planetary systems. Our models show that in relatively massive disks, several gas giants and rocky/icy planets emerge, migrate, and undergo dynamical instability. Secular perturbation between planets leads to orbital crossings, eccentricity excitation, and planetary ejection. In disks with modest masses, two or less gas giants form with multiple super-Earths. Orbital stability in these systems is generally maintained and they retain the kinematic structure after gas in their natal disks is depleted. These results reproduce the observed planetary mass-eccentricity and semimajor axis-eccentricity correlations. They also suggest that emerging gas giants can scatter residual cores to the outer disk regions. Subsequent in situ gas accretion onto these cores can lead to the formation of distant (> 30AU) gas giants with nearly circular orbits.

preprint2012arXiv

Internal Gravity Waves Modulate the Apparent Misalignment of Exoplanets around Hot Stars

We propose that the observed misalignment between extra-solar planets and their hot host stars can be explained by angular momentum transport within the host star. Observations have shown that this misalignment is preferentially around hot stars, which have convective cores and extended radiative envelopes. This situation is amenable to substantial angular momentum transport by internal gravity waves (IGW) generated at the convective-radiative interface. Here we present numerical simulations of this process and show that IGW can modulate the surface rotation of the star. With these two- dimensional simulations we show that IGW could explain the retrograde orbits observed in systems such as HAT-P-6 and HAT-P-7, however, extension to high obliquity objects will await future three- dimensional simulations. We note that these results also imply that individual massive stars should show temporal variations in their v sini measurements.

preprint2012arXiv

The Importance of Disk Structure in Stalling Type I Migration

As planets form they tidally interact with their natal disks. Though the tidal perturbation induced by Earth and super-Earth mass planets is generally too weak to significantly modify the structure of the disk, the interaction is potentially strong enough to cause the planets to undergo rapid type I migration. This physical process may provide a source of short-period super-Earths, though it may also pose a challenge to the emergence and retention of cores on long-period orbits with sufficient mass to evolve into gas giants. Previous numerical simulations have shown that the type I migration rate sensitively depends upon the circumstellar disk's properties, particularly the temperature and surface density gradients. Here, we derive these structure parameters for 1) a self-consistent viscous-disk model based on a constant α-prescription, 2) an irradiated disk model that takes into account heating due to the absorption of stellar photons, and 3) a layered-accretion disk model with variable α-parameter. We show that in the inner viscously-heated regions of typical protostellar disks, the horseshoe and corotation torques of super-Earths can exceed their differential Lindblad torque and cause them to undergo outward migration. However, the temperature profile due to passive stellar irradiation causes type I migration to be inwards throughout much of the disk. For disks in which there is outwards migration, we show that location and the mass range of the "planet traps" depends on some uncertain assumptions adopted for these disk models. Competing physical effects may lead to dispersion in super-Earths' mass-period distribution.

preprint2011arXiv

On the Migration of Jupiter and Saturn: Constraints from Linear Models of Secular Resonant Coupling with the Terrestrial Planets

We examine how the late divergent migration of Jupiter and Saturn may have perturbed the terrestrial planets. We identify six secular resonances between the nu_5 apsidal eigenfrequency of Jupiter and Saturn and the four eigenfrequencies of the terrestrial planets (g_{1-4}). We derive analytic upper limits on the eccentricity and orbital migration timescale of Jupiter and Saturn when these resonances were encountered to avoid perturbing the eccentricities of the terrestrial planets to values larger than the observed ones. If Jupiter and Saturn migrated with eccentricities comparable to their present day values, smooth migration with exponential timescales characteristic of planetesimal-driven migration (τ~5-10 Myr) would have perturbed the eccentricities of the terrestrial planets to values greatly exceeding the observed ones. This excitation may be mitigated if the eccentricity of Jupiter was small during the migration epoch, migration was very rapid (e.g. τ<~ 0.5 Myr perhaps via planet-planet scattering or instability-driven migration) or the observed small eccentricity amplitudes of the j=2,3 terrestrial modes result from low probability cancellation of several large amplitude contributions. Further, results of orbital integrations show that very short migration timescales (τ<0.5 Myr), characteristic of instability-driven migration, may also perturb the terrestrial planets' eccentricities by amounts comparable to their observed values. We discuss the implications of these constraints for the relative timing of terrestrial planet formation, giant planet migration, and the origin of the so-called Late Heavy Bombardment of the Moon 3.9+/-0.1 Ga ago. We suggest that the simplest way to satisfy these dynamical constraints may be for the bulk of any giant planet migration to be complete in the first 30-100 Myr of solar system history.

preprint2011arXiv

The Supernova Triggered Formation and Enrichment of Our Solar System

We investigate the enrichment of the pre-solar cloud core with short lived radionuclides (SLRs), especially 26Al. The homogeneity and the surprisingly small spread in the ratio 26Al/27Al observed in the overwhelming majority of calcium-aluminium-rich inclusions (CAIs) in a vast variety of primitive chondritic meteorites places strong constraints on the formation of the the solar system. Freshly synthesized radioactive 26Al has to be included and well mixed within 20kyr. After discussing various scenarios including X-winds, AGB stars and Wolf-Rayet stars, we come to the conclusion that triggering the collapse of a cold cloud core by a nearby supernova is the most promising scenario. We then narrow down the vast parameter space by considering the pre-explosion survivability of such a clump as well as the cross-section necessary for sufficient enrichment. We employ numerical simulations to address the mixing of the radioactively enriched SN gas with the pre-existing gas and the forced collapse within 20kyr. We show that a cold clump of 10Msun at a distance of 5pc can be sufficiently enriched in 26Al and triggered into collapse fast enough - within 18kyr after encountering the supernova shock - for a range of different metallicities and progenitor masses, even if the enriched material is assumed to be distributed homogeneously in the entire supernova bubble. In summary, we envision an environment for the birth place of the Solar System 4.567Gyr ago similar to the situation of the pillars in M16 nowadays, where molecular cloud cores adjacent to an HII region will be hit by a supernova explosion in the future. We show that the triggered collapse and formation of the Solar System as well as the required enrichment with radioactive 26Al are possible in this scenario.

preprint2010arXiv

A Population of Very-Hot Super-Earths in Multiple-Planet Systems Should be Uncovered by Kepler

We simulate a Kepler-like observation of a theoretical exoplanet population and we show that the observed orbital period distribution of the Kepler giant planet candidates is best matched by an average stellar specific dissipation function Q_* in the interval 10^6 ~< Q_* ~< 10^7. In that situation, the few super-Earths that are driven to orbital periods P < 1 day by dynamical interactions in multiple-planet systems will survive tidal disruption for a significant fraction of the main-sequence lifetimes of their stellar hosts. Consequently, though these very-hot super-Earths are not characteristic of the overall super-Earth population, their substantial transit probability implies that they should be significant contributors to the full super-Earth population uncovered by Kepler. As a result, the CoRoT-7 system may be the first representative of a population of very-hot super-Earths that we suggest should be found in multiple-planet systems preferentially orbiting the least-dissipative stellar hosts in the Kepler sample.

preprint2010arXiv

Binaries migrating in a gaseous disk: Where are the Galactic center binaries?

The massive stars in the Galactic center inner arcsecond share analogous properties with the so-called Hot Jupiters. Most of these young stars have highly eccentric orbits, and were probably not formed in-situ. It has been proposed that these stars acquired their current orbits from the tidal disruption of compact massive binaries scattered toward the proximity of the central supermassive black hole. Assuming a binary star formed in a thin gaseous disk beyond 0.1 pc from the central object, we investigate the relevance of disk-satellite interactions to harden the binding energy of the binary, and to drive its inward migration. A massive, equal-mass binary star is found to become more tightly wound as it migrates inwards toward the central black hole. The migration timescale is very similar to that of a single-star satellite of the same mass. The binary's hardening is caused by the formation of spiral tails lagging the stars inside the binary's Hill radius. We show that the hardening timescale is mostly determined by the mass of gas inside the binary's Hill radius, and that it is much shorter than the migration timescale. We discuss some implications of the binary's hardening process. When the more massive (primary) components of close binaries eject most their mass through supernova explosion, their secondary stars may attain a range of eccentricities and inclinations. Such processes may provide an alternative unified scenario for the origin of the kinematic properties of the central cluster and S-stars in the Galactic center as well as the high velocity stars in the Galactic halo.

preprint2010arXiv

Dusty Disks around White Dwarfs I: Origin of Debris Disks

A significant fraction of the mature FGK stars have cool dusty disks at least an orders of magnitudes brighter than the solar system's outer zodiacal light. Since such dusts must be continually replenished, they are generally assumed to be the collisional fragments of residual planetesimals analogous to the Kuiper Belt objects. At least 10% of solar type stars also bear gas giant planets. The fraction of stars with known gas giants or detectable debris disks (or both) appears to increase with the stellar mass. Here, we examine the dynamical evolution of systems of long-period gas giant planets and residual planetesimals as their host stars evolve off the main sequence, lose mass, and form planetary nebula around remnant white dwarf cores. The orbits of distant gas giant planets and super-km-size planetesimals expand adiabatically. During the most intense AGB mass loss phase, sub-meter-size particles migrate toward their host stars due to the strong hydrodynamical drag by the intense stellar wind. Along their migration paths, gas giant planets capture and sweep up sub-km-size planetesimals onto their mean-motion resonances. These planetesimals also acquire modest eccentricities which are determined by the mass of the perturbing planets, the rate and speed of stellar mass loss. The swept-up planetesimals undergo disruptive collisions which lead to the production of grains with an extended size range. The radiation drag on these particles is ineffective against the planets' resonant barrier and they form 30-to-150-AU-sizes rings which can effective reprocess the stellar irradiation in the form of FIR continuum. We identify the recently discovered dust ring around the white dwarf WD 2226-210 at the center of the Helix nebula as a prototype of such disks and suggest such rings may be common.

preprint2010arXiv

Embryo impacts and gas giant mergers I: Dichotomy of Jupiter and Saturn's core mass

Interior to the gaseous envelopes of Saturn, Uranus, and Neptune, there are high-density cores with masses larger than 10 Earth masses. According to the conventional sequential accretion hypothesis, such massive cores are needed for the onset of efficient accretion of their gaseous envelopes. However, Jupiter's gaseous envelope is more massive and core may be less massive than those of Saturn. In order to account for this structural diversity and the super-solar metallicity in the envelope of Jupiter and Saturn, we investigate the possibility that they may have either merged with other gas giants or consumed several Earth-mass proto-planetary embryos during or after the rapid accretion of their envelope. In general, impinging sub-Earth-mass planetesimals disintegrate in gas giants' envelopes deposit heavy elements well outside the cores and locally suppress the convection. Consequently, their fragments sediment to promote the growth of cores. Through a series of numerical simulations, we show that it is possible for colliding super-Earth-mass embryos to reach the cores of gas giants. Direct parabolic collisions also lead to the coalescence of gas giants and merging of their cores. In these cases, the energy released from the impact leads to vigorous convective motion throughout the envelope and the erosion of the cores. This dichotomy contributes to the observed dispersion in the internal structure and atmospheric composition between Jupiter and Saturn and other gas giant planets and elsewhere.

preprint2010arXiv

Toward a Deterministic Model of Planetary Formation VI: Dynamical Interaction and Coagulation of Multiple Rocky Embryos and Super-Earth Systems around Solar Type Stars

Radial velocity and transit surveys indicate that solar-type stars bear super-Earths, with mass and period up to ~ 20 M_E and a few months, are more common than those with Jupiter-mass gas giants. In many cases, these super-Earths are members of multiple-planet systems in which their mutual dynamical interaction has influenced their formation and evolution. In this paper, we modify an existing numerical population synthesis scheme to take into account protoplanetary embryos' interaction with their evolving natal gaseous disk, as well as their close scatterings and resonant interaction with each other. We show that it is possible for a group of compact embryos to emerge interior to the ice line, grow, migrate, and congregate into closely-packed convoys which stall in the proximity of their host stars. After the disk-gas depletion, they undergo orbit crossing, close scattering, and giant impacts to form multiple rocky Earths or super-Earths in non-resonant orbits around ~ 0.1AU with moderate eccentricities of ~0.01-0.1. We suggest that most refractory super-Earths with period in the range of a few days to weeks may have formed through this process. These super-Earths differ from Neptune-like ice giants by their compact sizes and lack of a substantial gaseous envelope.

preprint2010arXiv

Type I planet migration in nearly laminar disks - long term behavior

We carry out 2-D high resolution numerical simulations of type I planet migration with different disk viscosities. We find that the planet migration is strongly dependent on disk viscosities. Two kinds of density wave damping mechanisms are discussed. Accordingly, the angular momentum transport can be either viscosity dominated or shock dominated, depending on the disk viscosities. The long term migration behavior is different as well. Influences of the Rossby vortex instability on planet migration are also discussed. In addition, we investigate very weak shock generation in inviscid disks by small mass planets and compare the results with prior analytic results.

preprint2009arXiv

Protoplanetary migration in turbulent isothermal disks

In order to reproduce the statistical properties of the observed exoplanets, population synthesis models have shown that the migration of protoplanets should be significantly slowed down, and that processes stalling migration should be at work. Much current theoretical efforts have thus been dedicated to find physical effects that slow down, halt or even reverse migration. Many of these studies rely on the horseshoe drag, whose long term-evolution (saturated or not) is intimately related to the disk viscosity in laminar disk models. We investigate how the horseshoe drag exerted on a low-mass planet is altered by a more realistic treatment of the turbulence in protoplanetary disks. Two-dimensional hydrodynamic simulations are performed with a turbulence model that reproduces the main turbulence properties of three-dimensional magnetohydrodynamic calculations. We find that the horseshoe drag can remain unsaturated on the long term, depending on the turbulence strength. We show that the desaturation of the horseshoe drag by turbulence can be modeled by vortensity diffusion across the time-averaged planet's horseshoe region. At low-turbulence, the running-time averaged torque is in good agreement with the total torque obtained for an equivalent laminar model, with a similar vortensity's diffusion coefficient. At high-turbulence, differences arise due to the time-evolution of the averaged density profile with turbulence.

preprint2008arXiv

The Signature of the Ice Line and Modest Type I Migration in the Observed Exoplanet Mass-Semimajor Axis Distribution

Existing exoplanet radial velocity surveys are complete in the planetary mass-semimajor axis (Mp-a) plane over the range 0.1 AU < a < 2.0 AU where Mp >~ 100 M_Earth. We marginalize over mass in this complete domain of parameter space and demonstrate that the observed semimajor axis distribution is inconsistent with models of planet formation that use the full Type I migration rate derived from a linear theory and that do not include the effect of the ice line on the disk surface density profile. However, the efficiency of Type I migration can be suppressed by both nonlinear feedback and the barriers introduced by local maxima in the disk pressure distribution, and we confirm that the synthesized Mp-a distribution is compatible with the observed data if we account for both retention of protoplanetary embryos near the ice line and an order-of-magnitude reduction in the efficiency of Type I migration. The validity of these assumption can be checked because they also predict a population of short-period rocky planets with a range of masses comparable to that of the Earth as well as a "desert" in the Mp-a distribution centered around Mp ~ 30-50 M_Earth and a < 1 AU. We show that the expected "desert" in the Mp-a plane will be discernible by a radial velocity survey with 1 m/s precision and n ~ 700 radial velocity observations of program stars.

preprint2007arXiv

Hypervelocity binary stars: smoking gun of massive binary black holes

The hypervelocity stars recently found in the Galactic halo are expelled from the Galactic center through interactions between binary stars and the central massive black hole or between single stars and a hypothetical massive binary black hole. In this paper, we demonstrate that binary stars can be ejected out of the Galactic center with velocities up to 10^3 km/s, while preserving their integrity, through interactions with a massive binary black hole. Binary stars are unlikely to attain such high velocities via scattering by a single massive black hole or through any other mechanisms. Based on the above theoretical prediction, we propose a search for binary systems among the hypervelocity stars. Discovery of hypervelocity binary stars, even one, is a definitive evidence of the existence of a massive binary black hole in the Galactic center.

preprint2007arXiv

On the origin of kinematic distribution of the sub-parsec young stars in the Galactic center

Within a half-parsec from the Galactic center (GC), there is a population of coeval young stars which appear to reside in a coherent disk. Surrounding this dynamically-cool stellar system, there is a population of stars with a similar age and much larger eccentricities and inclinations relative to the disk. We propose a hypothesis for the origin of this dynamical dichotomy. Without specifying any specific mechanism, we consider the possibility that both stellar populations were formed within a disk some 6 Myr ago. But this orderly structure was dynamically perturbed outside-in by an intruding object with a mass ~10^4 Msun, which may be an intermediate-mass black hole (IMBH) or a dark stellar cluster hosting an IMBH. We suggest that the perturber migrated inward to ~0.15-0.3pc from the GC under the action of dynamical friction. Along the way, it captured many stars in the outer disk region into its mean-motion resonance, forced them to migrate with it, closely encountered with them, and induced the growth of their eccentricity and inclination. But stars in the inner regions of the disk retain their initial coplanar structure. We predict that some of the inclined and eccentric stars surrounding the disk may have similar Galactocentric semimajor axis. Future precision determination of their kinematic distribution of these stars will not only provide a test for this hypothesis but also evidences for the presence of an IMBH or a dark cluster at the immediate proximity of the massive black hole at the GC. (abridged)

preprint2007arXiv

Oxygen abundance of open cluster dwarfs

We present oxygen abundances of dwarfs in the young open cluster IC 4665 deduced from the OI $λ$7774 triplet lines and of dwarfs in the open cluster Pleiades derived from the [OI] $λ$6300 forbidden line. Stellar parameters and oxygen abundances were derived using the spectroscopic synthesis tool SME (Spectroscopy Made Easy). We find a dramatic increase in the upper boundary of the OI triplet abundances with decreasing temperature in the dwarfs of IC 4665, consistent with the trend found by Schuler et al. in the open clusters Pleiades and M 34, and to a less extent in the cool dwarfs of Hyades (Schuler et al. 2006a) and UMa (King & Schuler 2005). By contrast, oxygen abundances derived from the [OI] $λ$6300 forbidden line for stars in Pleiades and Hyades (Schuler et al. 2006b) are constant within the errors. Possible mechanisms that may lead a varying oxygen triplet line abundance are examined, including systematic errors in the stellar parameter determinations, the NLTE effects, surface activities and granulation. The age-related effects stellar surface activities (especially the chromospheric activities) are suggested by our analysis to blame for the large spreads of oxygen triplet line abundances.

preprint2007arXiv

Tidal Barrier and the Asymptotic Mass of Proto Gas-Giant Planets

Extrasolar planets found with radial velocity surveys have masses ranging from several Earth to several Jupiter masses. While mass accretion onto protoplanetary cores in weak-line T-Tauri disks may eventually be quenched by a global depletion of gas, such a mechanism is unlikely to have stalled the growth of some known planetary systems which contain relatively low-mass and close-in planets along with more massive and longer period companions. Here, we suggest a potential solution for this conundrum. In general, supersonic infall of surrounding gas onto a protoplanet is only possible interior to both of its Bondi and Roche radii. At a critical mass, a protoplanet's Bondi and Roche radii are equal to the disk thickness. Above this mass, the protoplanets' tidal perturbation induces the formation of a gap. Although the disk gas may continue to diffuse into the gap, the azimuthal flux across the protoplanets' Roche lobe is quenched. Using two different schemes, we present the results of numerical simulations and analysis to show that the accretion rate increases rapidly with the ratio of the protoplanet's Roche to Bondi radii or equivalently to the disk thickness. In regions with low geometric aspect ratios, gas accretion is quenched with relatively low protoplanetary masses. This effect is important for determining the gas-giant planets' mass function, the distribution of their masses within multiple planet systems around solar type stars, and for suppressing the emergence of gas-giants around low mass stars.

preprint2007arXiv

Tidal dissipation in rotating solar-type stars

We calculate the excitation and dissipation of low-frequency tidal oscillations in uniformly rotating solar-type stars. For tidal frequencies smaller than twice the spin frequency, inertial waves are excited in the convective envelope and are dissipated by turbulent viscosity. Enhanced dissipation occurs over the entire frequency range rather than in a series of very narrow resonant peaks, and is relatively insensitive to the effective viscosity. Hough waves are excited at the base of the convective zone and propagate into the radiative interior. We calculate the associated dissipation rate under the assumption that they do not reflect coherently from the center of the star. Tidal dissipation in a model based on the present Sun is significantly enhanced through the inclusion of the Coriolis force but may still fall short of that required to explain the circularization of close binary stars. However, the dependence of the results on the spin frequency, tidal frequency, and stellar model indicate that a more detailed evolutionary study including inertial and Hough waves is required. We also discuss the case of higher tidal frequencies appropriate to stars with very close planetary companions. The survival of even the closest hot Jupiters can be plausibly explained provided that the Hough waves they generate are not damped at the center of the star. We argue that this is the case because the tide excited by a hot Jupiter in the present Sun would marginally fail to achieve nonlinearity. As conditions at the center of the star evolve, nonlinearity may set in at a critical age, resulting in a relatively rapid inspiral of the hot Jupiter.

preprint2005arXiv

Towards a deterministic model of planetary formation. III. Mass distribution of short-period planets around stars of various masses

The origin of a recently discovered close-in Neptune-mass planet around GJ436 poses a challenge to the current theories of planet formation. Based on the sequential accretion hypothesis and the standard theory of gap formation and orbital migration, we show that around M dwarf stars, close-in Neptune-mass ice-giant planets may be relatively common, while close-in Jupiter-mass gas-giant planets are relatively rare. The mass distribution of close-in planets generally has two peaks at about Neptune mass and Jupiter mass. The lower-mass peak takes the maximum frequency for M dwarfs. Around more massive solar-type stars (G dwarfs), the higher-mass peak is much more pronounced. These are because planets tend to undergo type II migration after fully accreting gas around G dwarfs while they tend to migrate faster than gas accretion around M stars. Close-in Neptune-mass planets may also exist around G dwarfs, though they tend to be mostly composed of silicates and iron cores and their frequency is expected to be much smaller than that of Neptune-mass planets around M dwarfs and that of gas giants around G dwarfs. We also show that the conditions for planets' migration due to their tidal interaction with the disk and the stellar-mass dependence in the disk-mass distribution can be calibrated by the mass distribution of short-period planets around host stars with various masses.

preprint2003arXiv

On the Radii of Extrasolar Giant Planets

We have computed evolutionary models for extrasolar planets which range in mass from 0.1 to 3.0 Jovian Masses, and which range in equilibrium temperature from 113 K to 2000 K. We present four sequences of models, designed to show the structural effects of a solid core and of internal heating due to the conversion of kinetic to thermal energy at pressures of tens of bars. The model planetary radii are intended for comparisons with radii derived from observations of transiting extrasolar planets. To provide such comparisons, we expect that of order 10 transiting planets with orbital periods less than 200 days can be detected around bright (V<10) main-sequence stars for which accurate, well-sampled radial velocity measurements can be readily accumulated. Through these observations, structural properties of the planets will be derivable, particularly for low-mass, high-temperature planets. Implications regarding the transiting companion to OGLE-TR-56 recently announced by Konacki et al. are discussed. With regard to the confirmed transiting planet, HD 209458b, we find, in accordance with other recent calculations, that models without internal heating predict a radius that is 35 percent smaller than the observed radius. We explore the possibility that HD 209458b owes its large size to dissipation of energy arising from ongoing tidal circularization of the orbit. We show that residual scatter in the current radial velocity data set for HD 209458b is consistent with the presence of an as-of-yet undetected second companion, and that further radial velocity monitoring of the star is indicated.

preprint2002arXiv

Radial Flow of Dust Particles in Accretion Disks

We study the radial migration of dust particles in accreting protostellar disks analogous to the primordial solar nebula. This study takes account of the two dimensional (radial and normal) structure of the disk gas, including the effects of the variation in the gas velocity as a function of distance from the midplane. It is shown that the dust component of disks accretes slower than the gas component. At high altitude from the disk midplane, the gas rotates faster than particles because of the inward pressure gradient force, and its drag force causes particles to move outward in the radial direction. Viscous torque induces the gas within a scale height from the disk midplane to flow outward, carrying small (size < 100 micron at 10 AU) particles with it. Only particles at intermediate altitude or with sufficiently large sizes (> 1 mm at 10 AU) move inward. When the particles' radial velocities are averaged over the entire vertical direction, particles have a net inward flux. At large distances from the central star, particles migrate inward with a velocity much faster than the gas accretion velocity. However, their inward velocity is reduced below that of the gas in the inner regions of the disk. The rate of velocity decrease is a function of the particles' size. While larger particles retain fast accretion velocity until they approach closer to the star, 10 micron particles have slower velocity than the gas in the most part of the disk (r < 100 AU). This differential migration of particles causes the size fractionation. Dust disks composed mostly of small particles (size < 10 micron) accrete slower than gas disks, resulting in the increase in the dust-gas ratio during the gas accretion phase.