Source author record

Shigeru Ida

Shigeru Ida 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

64works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

64 published item(s)

preprint2022arXiv

A "no-drift" runaway pile-up of pebbles in protoplanetary disks II. Characteristics of the resulting planetesimal belt

Forming planetesimals from pebbles is a major challenge in our current understanding of planet formation. In a protoplanetary disk, pebbles drift inward near the disk midplane via gas drag and they may enter a dead zone. In this context, we identified that the backreaction of the drag of pebbles onto the gas could lead to a runaway pile-up of pebbles, the so-called no-drift mechanism. We improve upon the previous study of the no-drift mechanism by investigating the nature and characteristics of the resultant planetesimal belt. We performed 1D diffusion-advection simulations of drifting pebbles in the outer region of a dead zone by including the backreaction to the radial drift of pebbles and including planetesimal formation via the streaming instability. We considered the parameters that regulate gas accretion and vertical stirring of pebbles in the disk midplane. In this study, the pebble-to-gas mass flux ($F_{\rm p/g}$) was fixed as a parameter. We find that planetesimals initially form within a narrow ring whose width expands as accumulating pebbles radially diffuse over time. The system finally reaches a steady-state where the width of the planetesimal belt no longer changes. A non-negligible total mass of planetesimals (more than one Earth mass) is formed for a disk having $F_{\rm p/g} \gtrsim 0.1$ for more than $\sim 10-100$ kyr with nominal parameters: a gas mass flux of $\gtrsim10^{-8} {\rm M}_\oplus$/yr, $τ_{\rm s} \simeq 0.01-0.1$, $α_{\rm mid} \lesssim 10^{-4}$, and $α_{\rm acc} \simeq 10^{-3}-10^{-2}$ at $r \lesssim 10$ au, where $r$, $τ_{\rm s}$, $α_{\rm mid}$, and $α_{\rm acc}$ are the heliocentric distance, the Stokes number, and the parameters in a dead zone controlling the efficiencies of vertical turbulent diffusion of pebbles (i.e., scale height of pebbles) and gas accretion of the $α$-disk (i.e., gas surface density), respectively.

preprint2022arXiv

ALMA High-resolution Multiband Analysis for the Protoplanetary Disk around TW Hya

We present a high-resolution (2.5 au) multiband analysis of the protoplanetary disk around TW Hya using ALMA long baseline data at Bands 3, 4, 6, and 7. We aim to reconstruct a high-sensitivity millimeter continuum image and revisit the spectral index distribution. The imaging is performed by combining new ALMA data at Bands 4 and 6 with available archive data. Two methods are employed to reconstruct the images; multi-frequency synthesis (MFS) and the fiducial image-oriented method, where each band is imaged separately and the frequency dependence is fitted pixel by pixel. We find that the MFS imaging with the second order of Taylor expansion can reproduce the frequency dependence of the continuum emission between Bands 3 and 7 in a manner consistent with previous studies and is a reasonable method to reconstruct the spectral index map. The image-oriented method provides a spectral index map consistent with the MFS imaging, but with a two times lower resolution. Mock observations of an intensity model were conducted to validate the images from the two methods. We find that the MFS imaging provides a high-resolution spectral index distribution with an uncertainty of $<10$~\%. Using the submillimeter spectrum reproduced from our MFS images, we directly calculated the optical depth, power-law index of the dust opacity coefficient ($β$), and dust temperature. The derived parameters are consistent with previous works, and the enhancement of $β$ within the intensity gaps is also confirmed, supporting a deficit of millimeter-sized grains within the gaps.

preprint2022arXiv

Monte Carlo Simulation of Dust Particles in a Protoplanetary Disk: Crystalline to Amorphous Silicate Ratio in Comets

Observationally inferred crystalline abundance in silicates in comets, which should have been formed in the outer region of a protoplanetary disk, is relatively high (~ 10-60%), although crystalline silicates would be formed by annealing of amorphous precursors in the disk inner region. In order to quantitatively address this puzzle, we have performed Monte Carlo simulation of advection/diffusion of silicate particles in a turbulent disk, in the setting based on pebble accretion model: pebbles consisting of many small amorphous silicates embedded in icy mantle are formed in the disk outer region, silicate particles are released at the snow line, crystalline silicate particles are produced at the annealing line, the silicate particles diffused beyond the snow line, and they eventually stick to drifting pebbles to come back to the snow line. In a simple case without the sticking and with a steady pebble flux, we show through the simulations and analytical arguments that crystalline components in silicate materials beyond the snow line is robustly and uniformly ~ 5%. On the other hand, in a more realistic case with the sticking and with a decaying pebble flux, the crystalline abundance is raised up to ~ 20-25%, depending on the ratio of decay and diffusion timescales. This abundance is consistent with the observations. In this investigation, we assume a simple steady accretion disk. The simulations coupled with the disk evolution is needed for more detailed comparison with observed data.

preprint2022arXiv

Monte Carlo Simulation of Sugar Synthesis on Icy Dust Particles Intermittently Irradiated by UV in a Protoplanetary Disk

Context. While synthesis of organic molecules in molecular clouds or protoplanetary disks is complex, observations of interstellar grains, analyses of carbonaceous chondrites, and UV photochemistry experiments are rapidly developing and providing constraints on and clues to the complex organic molecule synthesis in space. It motivates us to construct a theoretical synthesis model. Aims. We develop a new code to simulate global reaction sequences of organic molecules to apply it for sugar synthesis by intermittent UV irradiation on the surface of icy particles in a protoplanetary disk. Here we show the first results of our new simulation. Methods. We apply a Monte Carlo method to select reaction sequences from all possible reactions, using the graph-theoretic matrix model for chemical reactions and modeling reactions on the icy particles during UV irradiation. Results. We here obtain the results consistent with the organic molecules in carbonaceous chondrites and obtained by the experiments, however, through a different pathway from the conventional formose reactions previously suggested. During UV irradiation, loosely-bonded O-rich large molecules are continuously created and destroyed. After UV irradiation is turned off, the ribose abundance rapidly increases, through the decomposition of the large molecules with break-ups of O-O bonds and replacements of C-OH by C-H to reach O/C = 1 for sugars. The sugar abundance is regulated mostly by the total atomic ratio H/O of starting materials, but not by their specific molecule forms. Deoxyribose is simultaneously synthesized, and most of the molecules end up with complex C-rich molecules.

preprint2021arXiv

A "no-drift" runaway pile-up of pebbles in protoplanetary disks in which midplane turbulence increases with radius

A notable challenge of planet formation is to find a path to directly form planetesimals from small particles. We aim to understand how drifting pebbles pile up in a protoplanetary disk with a non-uniform turbulence structure. We consider a disk structure in which the midplane turbulence viscosity is increasing with radius in protoplanetary disks as in the outer region of a dead zone. We perform 1D diffusion-advection simulations of pebbles that include back-reaction (the inertia) to radial drift and vertical/radial diffusion of pebbles for a given pebble-to-gas mass flux. We report a new mechanism, the "no-drift" runaway pile-up, leading to a runaway accumulation of pebbles in disks, thus favoring the formation of planetesimals by streaming and/or gravitational instabilities. This occurs when pebbles drifting in from the outer disk and entering a dead zone experience a decrease in vertical turbulence. The scale height of the pebble subdisk then decreases, and for small enough values of the turbulence in the dead zone and high values of the pebble to gas flux ratio, the back-reaction of pebbles on gas leads to a significant decrease in their drift velocity and thus their progressive accumulation. This occurs when the ratio of the flux of pebbles to that of the gas is large enough so that the effect dominates over any Kelvin-Helmholtz shear instability. This process is independent of the existence of a pressure bump.

preprint2021arXiv

Planetesimal formation around the snow line. II. Dust or pebbles?

Around the snow line, icy pebbles and silicate dust may locally pile-up and form icy and rocky planetesimals via streaming instability and/or gravitational instability. We perform 1D diffusion-advection simulations that include the back-reaction to radial drift and diffusion of icy pebbles and silicate dust, ice sublimation, release of silicate dust, and their recycling through recondensation and sticking onto pebbles outside the snow line. We use a realistic description of the scale height of silicate dust obtained from Ida et al. and that of pebbles including the effects of a Kelvin-Helmholtz instability. We study the dependence of solid pile-up on distinct effective viscous parameters for turbulent diffusions in the radial and vertical directions ($α_{\rm Dr}$ and $α_{\rm Dz}$) and for the gas accretion to the star ($α_{\rm acc}$) as well as that on the pebble-to-gas mass flux ($F_{\rm p/g}$). We derive the sublimation width of drifting icy pebbles which is a critical parameter to characterize the pile-up of silicate dust and pebbles around the snow line. We identify a parameter space (in the $F_{\rm p/g}-α_{\rm acc}-α_{\rm Dz}(=α_{\rm Dr})$ space) where pebbles no longer drift inward to reach the snow line due to the back-reaction that slows down radial velocity of pebbles. We show that the pile-up of solids around the snow line occurs in a broader range of parameters for $α_{\rm acc}=10^{-3}$ than for $α_{\rm acc}=10^{-2}$. Above a critical $F_{\rm p/g}$ value, the runaway pile-up of silicate dust inside the snow line is favored for $α_{\rm Dr}/α_{\rm acc} \ll 1$, while that of pebbles outside the snow line is favored for $α_{\rm Dr}/α_{\rm acc} \sim 1$. Our results imply that a distinct evolutionary path could produce a diversity of outcomes in terms of planetesimal formation around the snow line.

preprint2020arXiv

A new and simple prescription for planet orbital migration and eccentricity damping by planet-disc interactions based on dynamical friction

During planet formation gravitational interaction between a planetary embryo and the protoplanetary gas disc causes orbital migration of the planetary embryo, which plays an important role in shaping the final planetary system. While migration sometimes occurs in the supersonic regime, wherein the relative velocity between the planetary embryo and the gas is higher than the sound speed, migration prescriptions proposed thus far describing the planet-disc interaction force and the timescales of orbital change in the supersonic regime are inconsistent with one another. Here we discuss the details of existing prescriptions in the literature and derive a new simple and intuitive formulation for planet-disc interactions based on dynamical friction that can be applied in both supersonic and subsonic cases. While the existing prescriptions assume particular disc models, ours include the explicit dependence on the disc parameters; hence it can be applied to discs with any radial surface density and temperature dependence (except for the local variations with radial scales less than the disc scale height). Our prescription will reduce the uncertainty originating from different literature formulations of planet migration and will be an important tool to study planet accretion processes, especially when studying the formation of close-in low-mass planets that are commonly found in exoplanetary systems.

preprint2020arXiv

Orbital evolution of Saturn's satellites due to the interaction between the moons and massive Saturn's rings

Saturn's mid-sized moons (satellites) have a puzzling orbital configuration with trapping in mean-motion resonances with every other pairs (Mimas-Tethys 4:2 and Enceladus-Dione 2:1). To reproduce their current orbital configuration on the basis of Crida & Charnoz's model of satellite formation from a hypothetical ancient massive rings, adjacent pairs must pass 1st-order mean-motion resonances without being trapped. The trapping could be avoided by fast orbital migration and/or excitation of the satellite's eccentricity caused by gravitational interactions between the satellites and the rings (the disk), which are still unknown. In our research, we investigate the satellite orbital evolution due to interactions with the disk through full N-body simulations. We performed global high-resolution N-body simulations of a self-gravitating particle disk interacting with a single satellite. We used $N \sim 10^5$ particles for the disk. Gravitational forces of all the particles and their inelastic collisions are taken into account. As a result, dense short-wavelength wake structure is created by the disk self-gravity and global spiral arms with $m \sim$ a few is induced by the satellite. The self-gravity wakes regulate the orbital evolution of the satellite, which has been considered as a disk spreading mechanism but not as a driver for the orbital evolution. The self-gravity wake torque to the satellite is so effective that the satellite migration is much faster than that was predicted with the spiral arms torque. It provides a possible model to avoid the resonance capture of adjacent satellite pairs and establish the current orbital configuration of Saturn's mid-sized satellites.

preprint2020arXiv

Planetesimal formation around the snow line: I. Monte Carlo simulations of silicate dust pile-up in a turbulent disk

Context: The formation of rocky planetesimals is a long-standing problem in planet formation theory. One of the possibilities is that it results from gravitational instability as a result of pile-up of small silicate dust particles released from sublimating icy pebbles that pass the snow line. Aims: We want to understand and quantify the role of the water snow line for the formation of rock-rich and ice-rich planetesimals. In this paper, we focus on the formation of rock-rich planetesimals. A companion paper examines the combined formation of both rock-rich and ice-rich planetesimals. Methods: We develop a new Monte Carlo code to calculate the radial evolution of silicate particles in a turbulent accretion disk, accounting for the back-reaction (i.e., inertia) of the particles on their radial drift velocity and diffusion. Results depend in particular on the particle injection width (determined from the radial sublimation width of icy pebbles), the pebble scale height and the pebble mass flux through the disk. The scale height evolution of the silicate particles, which is the most important factor for the runaway pile-up, is automatically calculated in this Lagrange method. Results: From the numerical results, we derive semi-analytical relations for the scale height of the silicate dust particles and the particles-to-gas density ratio at the midplane, as functions of a pebble-to-gas mass flux ratio and the $α$ parameters for disk gas accretion and vertical/radial diffusion. We find that the runaway pile-up of the silicate particles (formation of rocky planetesimals) occurs if the pebble-to-gas mass flux ratio is $> [(α_{Dz}/α_{acc})/3 \times 10^{-2}]^{1/2}$ where $α_{Dz}$ and $α_{acc}$ are the $α$ parameters for vertical turbulent diffusion and disk gas accretion.

preprint2020arXiv

Uranian Satellite Formation by Evolution of a Water Vapor Disk Generated by a Giant Impact

The ice-giant planet Uranus likely underwent a giant impact, given that its spin axis is tilted by 98 degrees. That its satellite system is equally inclined and prograde suggests that it was formed as a consequence of the impact. However, the disks predicted by the impact simulations generally have sizes one order smaller and masses two orders larger than those of the observed system at present. Here we show, by means of a theoretical model, that the Uranian satellite formation is regulated by the evolution of the impact-generated disk. Because the vaporization temperature of water ice is low and both Uranus and the impactor are assumed to be ice-dominated, we can conclude that the impact-generated disk has mostly vaporized. We predict that the disk lost a significant amount of water vapour mass and spread to the levels of the current system until the disk cooled down enough for ice condensation and accretion of icy particles to begin. From the predicted distribution of condensed ices, our N-body simulation is able to reproduce the observed mass-orbit configuration of Uranian satellites. This scenario contrasts with the giant-impact model for the Earth's Moon, in which about half of the compact, impact-generated, solid or liquid disk is immediately incorporated into the Moon on impact.

preprint2019arXiv

The Galilean Satellites Formed Slowly from Pebbles

It is generally accepted that the four major (Galilean) satellites formed out of the gas disk that accompanied Jupiter's formation. However, understanding the specifics of the formation process is challenging as both small particles (pebbles) as well as the satellites are subject to fast migration processes. Here, we hypothesize a new scenario for the origin of the Galilean system, based on the capture of several planetesimal seeds and subsequent slow accretion of pebbles. To halt migration, we invoke an inner disk truncation radius, and other parameters are tuned for the model to match physical, dynamical, compositional, and structural constraints. In our scenario it is natural that Ganymede's mass is determined by pebble isolation. Our slow-pebble-accretion scenario then reproduces the following characteristics: (1) the mass of all the Galilean satellites; (2) the orbits of Io, Europa, and Ganymede captured in mutual 2:1 mean motion resonances; (3) the ice mass fractions of all the Galilean satellites; (4) the unique ice-rock partially differentiated Callisto and the complete differentiation of the other satellites. Our scenario is unique to simultaneously reproduce these disparate properties.

preprint2016arXiv

A Pair of Giant Planets around the Evolved Intermediate-Mass Star HD 47366: Multiple Circular Orbits or a Mutually Retrograde Configuration

We report the detection of a double planetary system around the evolved intermediate-mass star HD 47366 from precise radial-velocity measurements at Okayama Astrophysical Observatory, Xinglong Station, and Australian Astronomical Observatory. The star is a K1 giant with a mass of 1.81+-0.13M_sun, a radius of 7.30+-0.33R_sun, and solar metallicity. The planetary system is composed of two giant planets with minimum mass of 1.75^{+0.20}_{-0.17}Mjup and 1.86^{+0.16}_{-0.15}Mjup, orbital period of 363.3^{+2.5}_{-2.4} d and 684.7^{+5.0}_{-4.9} d, and eccentricity of 0.089^{+0.079}_{-0.060} and 0.278^{+0.067}_{-0.094}, respectively, which are derived by a double Keplerian orbital fit to the radial-velocity data. The system adds to the population of multi-giant-planet systems with relatively small orbital separations, which are preferentially found around evolved intermediate-mass stars. Dynamical stability analysis for the system revealed, however, that the best-fit orbits are unstable in the case of a prograde configuration. The system could be stable if the planets were in 2:1 mean-motion resonance, but this is less likely considering the observed period ratio and eccentricity. A present possible scenario for the system is that both of the planets have nearly circular orbits, namely the eccentricity of the outer planet is less than ~0.15, which is just within 1.4sigma of the best-fit value, or the planets are in a mutually retrograde configuration with a mutual orbital inclination larger than 160 degree.

preprint2016arXiv

ALMA Observations of a Gap and a Ring in the Protoplanetary Disk around TW Hya

We report the first detection of a gap and a ring in 336 GHz dust continuum emission from the protoplanetary disk around TW Hya, using the Atacama Large Millimeter/Submillimeter Array (ALMA). The gap and ring are located at around 25 and 41 au from the central star, respectively, and are associated with the CO snow line at ~30 au. The gap has a radial width of less than 15 au and a mass deficit of more than 23%, taking into account that the observations are limited to an angular resolution of ~15 au. In addition, the 13CO and C18O J = 3 - 2 lines show a decrement in CO line emission throughout the disk, down to ~10 au, indicating a freeze-out of gas-phase CO onto grain surfaces and possible subsequent surface reactions to form larger molecules. The observed gap could be caused by gravitational interaction between the disk gas and a planet with a mass less than super-Neptune (2M_{Neptune}), or could be the result of the destruction of large dust aggregates due to the sintering of CO ice.

preprint2016arXiv

Dust and gas density evolution at a radial pressure bump in protoplanetary disks

We investigate the simultaneous evolution of dust and gas density profiles at a radial pressure bump located in a protoplanetary disk. If dust particles are treated as test particles, a radial pressure bump traps dust particles that drift radially inward. As the dust particles become more concentrated at the gas pressure bump, however, the drag force from dust to gas (back-reaction), which is ignored in a test-particle approach, deforms the pressure bump. We find that the pressure bump is completely deformed by the back-reaction when the dust-to-gas mass ratio reaches $\sim 1$ for a slower bump restoration. The direct gravitational instability of dust particles is inhibited by the bump destruction. In the dust-enriched region, the radial pressure support becomes $\sim 10-100$ times lower than the global value set initially. Although the pressure bump is a favorable place for streaming instability (SI), the flattened pressure gradient inhibits SI from forming large particle clumps corresponding to $100-1000$ km sized bodies, which has been previously proposed. If SI occurs there, the dust clumps formed would be $10-100$ times smaller, that is, of about $1 - 100$ km.

preprint2016arXiv

Free-floating planets from core accretion theory: microlensing predictions

We calculate the microlensing event rate and typical time-scales for the free-floating planet (FFP) population that is predicted by the core accretion theory of planet formation. The event rate is found to be ~$1.8\times 10^{-3}$ of that for the stellar population. While the stellar microlensing event time-scale peaks at around 20 days, the median time-scale for FFP events (~0.1 day) is much shorter. Our values for the event rate and the median time-scale are significantly smaller than those required to explain the \cite{Sum+11} result, by factors of ~13 and ~16, respectively. The inclusion of planets at wide separations does not change the results significantly. This discrepancy may be too significant for standard versions of both the core accretion theory and the gravitational instability model to explain satisfactorily. Therefore, either a modification to the planet formation theory is required, or other explanations to the excess of short-time-scale microlensing events are needed. Our predictions can be tested by ongoing microlensing experiment such as KMTNet, and by future satellite missions such as WFIRST and Euclid.

preprint2016arXiv

High-Contrast Imaging of Intermediate-Mass Giants with Long-Term Radial Velocity Trends

A radial velocity (RV) survey for intermediate-mass giants has been operated for over a decade at Okayama Astrophysical Observatory (OAO). The OAO survey has revealed that some giants show long-term linear RV accelerations (RV trends), indicating the presence of outer companions. Direct imaging observations can help clarify what objects generate these RV trends. We present the results of high-contrast imaging observations or six intermediate-mass giants with long-term RV trends using the Subaru Telescope and HiCIAO camera. We detected co-moving companions to $γ$ Hya B ($0.61^{+0.12}_{-0.14} M_\odot$), HD 5608 B ($0.10 \pm 0.01 M_\odot$), and HD 109272 B ($0.28 \pm 0.06 M_\odot$). For the remaining targets($ι$ Dra, 18 Del, and HD 14067) we exclude companions more massive than 30-60 $M_\mathrm{Jup}$ at projected separations of 1arcsec-7arcsec. We examine whether these directly imaged companions or unidentified long-period companions can account for the RV trends observed around the six giants. We find that the Kozai mechanism can explain the high eccentricity of the inner planets $ι$ Dra b, HD 5608 b, and HD 14067 b.

preprint2016arXiv

On the water delivery to terrestrial embryos by ice pebble accretion

Standard accretion disk models suggest that the snow line in the solar nebula migrated interior to the Earth's orbit in a late stage of nebula evolution. In this late stage, a significant amount of ice could have been delivered to 1 AU from outer regions in the form of mm to dm-sized pebbles. This raises the question why the present Earth is so depleted of water (with the ocean mass being as small as 0.023% of the Earth mass). Here we quantify the amount of icy pebbles accreted by terrestrial embryos after the migration of the snow line assuming that no mechanism halts the pebble flow in outer disk regions. We use a simplified version of the coagulation equation to calculate the formation and radial inward drift of icy pebbles in a protoplanetary disk. The pebble accretion cross section of an embryo is calculated using analytic expressions presented by recent studies. We find that the final mass and water content of terrestrial embryos strongly depends on the radial extent of the gas disk, the strength of disk turbulence, and the time at which the snow lines arrives at 1 AU. The disk's radial extent sets the lifetime of the pebble flow, while turbulence determines the density of pebbles at the midplane where the embryos reside. We find that the final water content of the embryos falls below 0.023 wt% only if the disk is compact (< 100 AU), turbulence is strong at 1 AU, and the snow line arrives at 1 AU later than 2-4 Myr after disk formation. If the solar nebula extended to 300 AU, initially rocky embryos would have evolved into icy planets of 1-10 Earth masses unless the snow-line migration was slow. If the proto-Earth contained water of ~ 1 wt% as might be suggested by the density deficit of the Earth's outer core, the formation of the proto-Earth was possible with weaker turbulence and with earlier (> 0.5-2 Myr) snow-line migration.

preprint2016arXiv

The Giant Impact Simulations with Density Independent Smoothed Particle Hydrodynamics

At present, the giant impact (GI) is the most widely accepted model for the origin of the Moon. Most of the numerical simulations of GI have been carried out with the smoothed particle hydrodynamics (SPH) method. Recently, however, it has been pointed out that standard formulation of SPH (SSPH) has difficulties in the treatment of a contact discontinuity such as a core-mantle boundary and a free surface such as a planetary surface. This difficulty comes from the assumption of differentiability of density in SSPH. We have developed an alternative formulation of SPH, density independent SPH (DISPH), which is based on differentiability of pressure instead of density to solve the problem of a contact discontinuity. In this paper, we report the results of the GI simulations with DISPH and compare them with those obtained with SSPH. We found that the disk properties, such as mass and angular momentum produced by DISPH is different from that of SSPH. In general, the disks formed by DISPH are more compact: while formation of a smaller mass moon for low-oblique impacts is expected with DISPH, inhibition of ejection would promote formation of a larger mass moon for high-oblique impacts. Since only the improvement of core-mantle boundary significantly affects the properties of circumplanetary disks generated by GI and DISPH has not been significantly improved from SSPH for a free surface, we should be very careful when some conclusions are drawn from the numerical simulations for GI. And it is necessary to develop the numerical hydrodynamical scheme for GI that can properly treat the free surface as well as the contact discontinuity.

preprint2016arXiv

The radial dependence of pebble accretion rates: A source of diversity in planetary systems I. Analytical formulation

Context. The classical "planetesimal" accretion scenario for the formation of planets has recently evolved with the idea that "pebbles", centimeter- to meter-sized icy grains migrating in protoplanetary disks, can control planetesimal and/or planetary growth. Aims. We investigate how pebble accretion depends on disk properties and affects the formation of planetary systems Methods. We construct analytical models of pebble accretion onto planetary embryos that consistently account for the mass and orbital evolution of the pebble flow and reflect disk structure. Results. We derive simple formulas for pebble accretion rates in the so-called "settling" regime for planetary embryos with more than 100 km in size. For relatively smaller embryos or in outer disk regions, the accretion mode is 3D, meaning that the thickness of the pebble flow must be taken into account, and resulting in an accretion rate that is independent of the embryo mass. For larger embryos or in inner regions, the accretion is in a 2D mode, i.e., the pebble disk may be considered to be infinitely thin. We show that the radial dependence of the pebble accretion rata is different (even the sign of the power-law exponent changes) for different disk conditions such as the disk heating source (viscous heating or stellar irradiation), the drag law (Stokes or Epstein, and weak or strong coupling), and in the 2D or 3D accretion modes. We also discuss the effect of the sublimation/destruction of icy pebbles inside the snow line. Conclusions. Pebble accretion easily produces a large diversity of planetary systems. In other words, to correctly infer the results of planet formation through pebble accretion, detailed prescriptions of disk evolution and pebble growth, sublimation/destruction and migration are required.

preprint2015arXiv

A semi-analytical model for exploring Galilean satellites formation from a massive disk

A better knowledge of Jovian satellites' origins will bring light on the environment that surrounded Jupiter during its formation and can help us to understand the characteristics of this unique satellite system. We developed a semi-analytical model to investigate Jupiter's regular satellite formation and present the results of our population synthesis calculations. We performed simulations adopting a massive, static, low-viscosity circumplanetary disk model, in agreement with a current study of magnetorotational instability in a circum-planetary disk. We find that the high gas density leads to very rapid migration of satellitesimals due to gas drag and type II migration of satellites in a faster disk-dominated mode. A large concentration of solids, large building blocks and longer type II migration time-scales favor formation and survival of large satellites. However, bodies as massive as Ganymede and those located far away from Jupiter, such as Callisto, are difficult to form with this scenario.

preprint2015arXiv

Eccentricity Evolution Through Accretion of Protoplanets

Most of super-Earths detected by the radial velocity (RV) method have significantly smaller eccentricities than the eccentricities corresponding to velocity dispersion equal to their surface escape velocity ("escape eccentricities"). If orbital instability followed by giant impacts among protoplanets that have migrated from outer region is considered, it is usually considered that eccentricities of the merged bodies become comparable to those of orbital crossing bodies, which are excited up to their escape eccentricities by close scattering. However, the eccentricity evolution in the {\it in situ} accretion model has not been studied in detail. Here, we investigate the eccentricity evolution through {\it N}-body simulations. We have found that the merged planets tend to have much smaller eccentricities than the escape eccentricities due to very efficient collision damping. If the protoplanet orbits are initially well separated and their eccentricities are securely increased, an inner protoplanet collides at its apocenter with an outer protoplanet at its pericenter. The eccentricity of the merged body is the smallest for such configuration. Orbital inclinations are also damped by this mechanism and planets tend to share a same orbital plane, which is consistent with {\it Kepler} data. Such efficient collision damping is not found when we start calculations from densely packed orbits of the protoplanets. If the protoplanets are initially in the mean-motion resonances, which corresponds to well separated orbits, the {\it in situ} accretion model well reproduces the features of eccentricities and inclinations of multiple super-Earths/Earth systems discovered by RV and {\it Kepler} surveys.

preprint2015arXiv

Effects of Dynamical Evolution of Giant Planets on the Delivery of Atmophile Elements During Terrestrial Planet Formation

Recent observations started revealing the compositions of protostellar discs and planets beyond the Solar System. In this paper, we explore how the compositions of terrestrial planets are affected by dynamical evolution of giant planets. We estimate the initial compositions of building blocks of these rocky planets by using a simple condensation model, and numerically study the compositions of planets formed in a few different formation models of the Solar System. We find that the abundances of refractory and moderately volatile elements are nearly independent of formation models, and that all the models could reproduce the abundances of these elements of the Earth. The abundances of atmophile elements, on the other hand, depend on the scattering rate of icy planetesimals into the inner disc as well as the mixing rate of the inner planetesimal disc. For the classical formation model, neither of these mechanisms are efficient and the accretion of atmophile elements during the final assembly of terrestrial planets appears to be difficult. For the Grand Tack model, both of these mechanisms are efficient, which leads to a relatively uniform accretion of atmophile elements in the inner disc. It is also possible to have a "hybrid" scenario where the mixing is not very efficient but the scattering is efficient. The abundances of atmophile elements in this case increases with orbital radii. Such a scenario may occur in some of the extrasolar planetary systems which are not accompanied by giant planets or those without strong perturbations from giants. We also confirm that the Grand Tack scenario leads to the distribution of asteroid analogues where rocky planetesimals tend to exist interior to icy ones, and show that their overall compositions are consistent with S-type and C-type chondrites, respectively.

preprint2015arXiv

Excavation and Melting of the Hadean Continental Crust by Late Heavy Bombardment

No Hadean rocks have ever been found on Earth's surface except for zircons---evidence of continental crust, suggesting that Hadean continental crust existed but later disappeared. One hypothesis for the disappearance of the continental crust is excavation/melting by the Late Heavy Bombardment (LHB), a concentration of impacts in the last phase of the Hadean eon. In this paper, we calculate the effects of LHB on Hadean continental crust in order to investigate this hypothesis. Approximating the size-frequency distribution of the impacts by a power-law scaling with an exponent α as a parameter, we have derived semi-analytical expressions for the effects of LHB impacts. We calculated the total excavation/melting volume and area affected by the LHB from two constraints of LHB on the moon, the size of the largest basin during LHB, and the density of craters larger than 20 km. We also investigated the effects of the value of α. Our results show that LHB does not excavate/melt all of Hadean continental crust directly, but over 70% of the Earth's surface area can be covered by subsequent melts in a broad range of α. If there have been no overturns of the continental crust until today, LHB could be responsible for the absence of Hadean rocks because most of Hadean continental crust is not be exposed on the Earth's surface in this case.

preprint2015arXiv

Five New Exoplanets Orbiting Three Metal-Rich, Massive Stars: Two-Planet Systems Including Long-Period Planets, and an Eccentric Planet

We report detections of new exoplanets from a radial velocity (RV) survey of metal-rich FGK stars by using three telescopes. By optimizing our RV analysis method to long time-baseline observations, we have succeeded in detecting five new Jovian-planets around three metal-rich stars HD 1605, HD 1666, and HD 67087 with the masses of $1.3 M_{\odot}$, $1.5 M_{\odot}$, and $1.4 M_{\odot}$, respectively. A K1 subgiant star HD 1605 hosts two planetary companions with the minimum masses of $ M_p \sin i = 0.96 M_{\mathrm{JUP}}$ and $3.5 M_{\mathrm{JUP}}$ in circular orbits with the planets' periods $P = 577.9$ days and $2111$ days, respectively. HD 1605 shows a significant linear trend in RVs. Such a system consisting of Jovian planets in circular orbits has rarely been found and thus HD 1605 should be an important example of a multi-planetary system that is likely unperturbed by planet-planet interactions. HD 1666 is a F7 main sequence star which hosts an eccentric and massive planet of $ M_p \sin i = 6.4 M_{\mathrm{JUP}}$ in the orbit with $a_{\rm p} = 0.94$ AU and an eccentricity $e=0.63$. Such an eccentric and massive planet can be explained as a result of planet-planet interactions among Jovian planets. While we have found the large residuals of $\mathrm{rms} = 35.6\ \mathrm{m\ s^{-1}}$, the periodogram analysis does not support any additional periodicities. Finally, HD 67087 hosts two planets of $ M_p \sin i = 3.1 M_{\mathrm{JUP}}$ and $4.9 M_{\mathrm{JUP}}$ in orbits with $P=352.2$ days and $2374$ days, and $e=0.17$ and $0.76$, respectively. Although the current RVs do not lead to accurate determinations of its orbit and mass, HD 67087 c can be one of the most eccentric planets ever discovered in multiple systems.

preprint2015arXiv

Formation of Phobos and Deimos via a Giant Impact

Although the two moons of Mars, Phobos and Deimos, have long been thought to be captured asteroids, recent observations of their compositions and orbits suggest that they may have formed from debris generated by one or more giant impacts of bodies with ~ 0.01 x target mass. Recent studies have both analytically estimated debris produced by giant impacts on Mars and numerically examined the evolution of circum-Mars debris disks. We perform a numerical study (Smoothed Particle Hydrodynamics simulation) of debris retention from giant impacts onto Mars, particularly in relation to a Borealis-scale giant impact (E ~ 3 x 10^29 J) capable of producing the Borealis basin. We find that a Borealis-scale impact is capable of producing a disk of mass ~ 5 x 10^20 kg (~ 1 - 4 % of the impactor mass), sufficient debris to form at least one of the martian moons according to recent numerical studies of martian debris disk evolution. While a Borealis-scale impact may generate sufficient debris to form both Phobos and Deimos, further studies of the debris disk evolution are necessary. Our results can serve as inputs for future studies of martian debris disk evolution.

preprint2015arXiv

Temporary Capture of Asteroids by a Planet: Dependence of Prograde/Retrograde Capture on Asteroids' Semimajor Axes

We have investigated the dependence of the prograde/retrograde temporary capture of asteroids by a planet on their original heliocentric semimajor axes through analytical arguments and numerical orbital integrations in order to discuss the origins of irregular satellites of giant planets. We found that capture is mostly retrograde for the asteroids near the planetary orbit and is prograde for those from further orbits. An analytical investigation reveals the intrinsic dynamics of these dependences and gives boundary semimajor axes for the change in prograde/retrograde capture. The numerical calculations support the idea of deriving the analytical formulae and confirm their dependence. Our numerical results show that the capture probability is much higher for bodies from the inner region than for outer ones. These results imply that retrograde irregular satellites of Jupiter are most likely to be captured bodies from the nearby orbits of Jupiter that may have the same origin as Trojan asteroids, while prograde irregular satellites originate from far inner regions such as the main-belt asteroid region.

preprint2014arXiv

A Long-period Eccentric Substellar Companion to the Evovled Intermediate-Mass Star HD 14067

We report the detection of a substellar companion orbiting an evolved intermediate-mass ($M_\star=2.4\,M_\odot$) star HD 14067 (G9 III) using precise Doppler technique. Radial velocities of this star can be well fitted either by a periodic Keplerian variation with a decreasing linear velocity trend (P=1455 days, $K_1=92.2$ m s$^{-1}$, $e=0.533$, and $\dotγ=-22.4$ m s$^{-1}$ yr$^{-1}$) or a single Keplerian orbit without linear trend (P=2850 days, $K_1=100.1$ m s$^{-1}$, and $e=0.697$). The minimum mass ($m_2\sin{i}=7.8\,M_{\rm J}$ for the model with a linear trend, or $m_2\sin{i}=9.0\,M_{\rm J}$ for the model without a linear trend) suggests a long-period giant planet around an evolved intermediate-mass star. The eccentricity of the orbit is among the highest known for planets ever detected around evolved stars.

preprint2014arXiv

On the filtering and processing of dust by planetesimals 1. Derivation of collision probabilities for non-drifting planetesimals

Context. Circumstellar disks are known to contain a significant mass in dust ranging from micron to centimeter size. Meteorites are evidence that individual grains of those sizes were collected and assembled into planetesimals in the young solar system. Aims. We assess the efficiency of dust collection of a swarm of non-drifting planetesimals {\rev with radii ranging from 1 to $10^3$\,km and beyond. Methods. We calculate the collision probability of dust drifting in the disk due to gas drag by planetesimal accounting for several regimes depending on the size of the planetesimal, dust, and orbital distance: the geometric, Safronov, settling, and three-body regimes. We also include a hydrodynamical regime to account for the fact that small grains tend to be carried by the gas flow around planetesimals. Results. We provide expressions for the collision probability of dust by planetesimals and for the filtering efficiency by a swarm of planetesimals. For standard turbulence conditions (i.e., a turbulence parameter $α=10^{-2}$), filtering is found to be inefficient, meaning that when crossing a minimum-mass solar nebula (MMSN) belt of planetesimals extending between 0.1 AU and 35 AU most dust particles are eventually accreted by the central star rather than colliding with planetesimals. However, if the disk is weakly turbulent ($α=10^{-4}$) filtering becomes efficient in two regimes: (i) when planetesimals are all smaller than about 10 km in size, in which case collisions mostly take place in the geometric regime; and (ii) when planetary embryos larger than about 1000 km in size dominate the distribution, have a scale height smaller than one tenth of the gas scale height, and dust is of millimeter size or larger in which case most collisions take place in the settling regime. These two regimes have very different properties: we find that the local filtering efficiency $x_{filter,MMSN}$ scales with $r^{-7/4}$ (where $r$ is the orbital distance) in the geometric regime, but with $r^{-1/4}$ to $r^{1/4}$ in the settling regime. This implies that the filtering of dust by small planetesimals should occur close to the central star and with a short spread in orbital distances. On the other hand, the filtering by embryos in the settling regime is expected to be more gradual and determined by the extent of the disk of embryos. Dust particles much smaller than millimeter size tend only to be captured by the smallest planetesimals because they otherwise move on gas streamlines and their collisions take place in the hydrodynamical regime. Conclusions. Our results hint at an inside-out formation of planetesimals in the infant solar system because small planetesimals in the geometrical limit can filter dust much more efficiently close to the central star. However, even a fully-formed belt of planetesimals such as the MMSN only marginally captures inward-drifting dust and this seems to imply that dust in the protosolar disk has been filtered by planetesimals even smaller than 1 km (not included in this study) or that it has been assembled into planetesimals by other mechanisms (e.g., orderly growth, capture into vortexes). Further refinement of our work concerns, among other things: a quantitative description of the transition region between the hydro and settling regimes; an assessment of the role of disk turbulence for collisions, in particular in the hydro regime; and the coupling of our model to a planetesimal formation model.

preprint2014arXiv

The Transiting Exoplanet Survey Satellite

The Transiting Exoplanet Survey Satellite (TESS) will search for planets transiting bright and nearby stars. TESS has been selected by NASA for launch in 2017 as an Astrophysics Explorer mission. The spacecraft will be placed into a highly elliptical 13.7-day orbit around the Earth. During its two-year mission, TESS will employ four wide-field optical CCD cameras to monitor at least 200,000 main-sequence dwarf stars with I = 4-13 for temporary drops in brightness caused by planetary transits. Each star will be observed for an interval ranging from one month to one year, depending mainly on the star's ecliptic latitude. The longest observing intervals will be for stars near the ecliptic poles, which are the optimal locations for follow-up observations with the James Webb Space Telescope. Brightness measurements of preselected target stars will be recorded every 2 min, and full frame images will be recorded every 30 min. TESS stars will be 10-100 times brighter than those surveyed by the pioneering Kepler mission. This will make TESS planets easier to characterize with follow-up observations. TESS is expected to find more than a thousand planets smaller than Neptune, including dozens that are comparable in size to the Earth. Public data releases will occur every four months, inviting immediate community-wide efforts to study the new planets. The TESS legacy will be a catalog of the nearest and brightest stars hosting transiting planets, which will endure as highly favorable targets for detailed investigations.

preprint2013arXiv

Do Giant Planets Survive Type II Migration?

Planetary migration is one of the most serious problems to systematically understand the observations of exoplanets. We clarify that the theoretically predicted type II migration is too fast, as well as type I migration, by developing detailed analytical arguments in which the timescale of type II migration is compared with the disk lifetime. In the disk-dominated regime, the type II migration timescale is characterized by a {\it local} viscous diffusion timescale, while the disk lifetime characterized by a {\it global} diffusion timescale that is much longer than the local one. Even in the planet-dominated regime where the inertia of the planet mass reduces the migration speed, the timescale is still shorter than the disk lifetime except in the final disk evolution stage where the total disk mass decays below the planet mass. This suggests that most giant planets plunge into the central stars within the disk lifetime, and it contradicts the exoplanet observations that gas giants are piled up at $r \ga 1$AU. We examine additional processes that may arise in protoplanetary disks: dead zones, photoevaporation of gas, and gas flow across a gap formed by a type II migrator. Although they make the type II migration timescale closer to the disk lifetime, we show that none of them can act as an effective barrier for rapid type II migration with the current knowledge of these processes. We point out that gas flow across a gap and the fraction of the flow accreted onto the planets are uncertain and they may have a potential to solve the problem. Much more detailed investigation for each process may be needed to explain the observed distribution of gas giants in extrasolar planetary systems.

preprint2013arXiv

Effects of Dynamical Evolution of Giant Planets on Survival of Terrestrial Planets

The orbital distributions of currently observed extrasolar giant planets allow marginally stable orbits for hypothetical, terrestrial planets. In this paper, we propose that many of these systems may not have additional planets on these "stable" orbits, since past dynamical instability among giant planets could have removed them. We numerically investigate the effects of early evolution of multiple giant planets on the orbital stability of the inner, sub-Neptune-like planets which are modeled as test particles, and determine their dynamically unstable region. Previous studies have shown that the majority of such test particles are ejected out of the system as a result of close encounters with giant planets. Here, we show that secular perturbations from giant planets can remove test particles at least down to 10 times smaller than their minimum pericenter distance. Our results indicate that, unless the dynamical instability among giant planets is either absent or quiet like planet-planet collisions, most test particles down to ~0.1 AU within the orbits of giant planets at a few AU may be gone. In fact, out of ~30 % of survived test particles, about three quarters belong to the planet-planet collision cases. We find a good agreement between our numerical results and the secular theory, and present a semi-analytical formula which estimates the dynamically unstable region of the test particles just from the evolution of giant planets. Finally, our numerical results agree well with the observations, and also predict the existence of hot rocky planets in eccentric giant planet systems.

preprint2013arXiv

Planetary Companions to Three Evolved Intermediate-Mass Stars: HD 2952, HD 120084, and omega Serpentis

We report the detections of planetary companions orbiting around three evolved intermediate-mass stars from precise radial velocity measurements at Okayama Astrophysical Observatory. HD 2952 (K0III, 2.5 M_sun) and omega Ser (G8III, 2.2 M_sun) host a relatively low mass planet with minimum mass of m_2sin i=1.6 M_J and 1.7 M_J in nearly circular orbits with period of P=312 and 277 d, respectively. HD 120084 (G7 III, 2.4 M_sun) hosts an eccentric planet with m_2sin i=4.5 M_J in an orbit with P=2082 d and eccentricity of e=0.66. The planet has one of the largest eccentricities among those ever discovered around evolved intermediate-mass stars, almost all of which have eccentricity smaller than 0.4. We also show that radial velocity variations of stellar oscillations for G giants can be averaged out below a level of a few m/s at least in timescale of a week by high cadence observations, which enables us to detect a super-Earth and a Neptune-mass planet in short-period orbits even around such giant stars.

preprint2013arXiv

The Gemini NICI Planet-Finding Campaign: The Frequency of Giant Planets around Young B and A Stars

We have carried out high contrast imaging of 70 young, nearby B and A stars to search for brown dwarf and planetary companions as part of the Gemini NICI Planet-Finding Campaign. Our survey represents the largest, deepest survey for planets around high-mass stars (~1.5-2.5 M_sun) conducted to date and includes the planet hosts beta Pic and Fomalhaut. We obtained follow-up astrometry of all candidate companions within 400 AU projected separation for stars in uncrowded fields and identified new low-mass companions to HD 1160 and HIP 79797. We have found that the previously known young brown dwarf companion to HIP 79797 is itself a tight (3 AU) binary, composed of brown dwarfs with masses 58 (+21, -20) M_Jup and 55 (+20, -19) M_Jup, making this system one of the rare substellar binaries in orbit around a star. Considering the contrast limits of our NICI data and the fact that we did not detect any planets, we use high-fidelity Monte Carlo simulations to show that fewer than 20% of 2 M_sun stars can have giant planets greater than 4 M_Jup between 59 and 460 AU at 95% confidence, and fewer than 10% of these stars can have a planet more massive than 10 M_Jup between 38 and 650 AU. Overall, we find that large-separation giant planets are not common around B and A stars: fewer than 10% of B and A stars can have an analog to the HR 8799 b (7 M_Jup, 68 AU) planet at 95% confidence. We also describe a new Bayesian technique for determining the ages of field B and A stars from photometry and theoretical isochrones. Our method produces more plausible ages for high-mass stars than previous age-dating techniques, which tend to underestimate stellar ages and their uncertainties.

preprint2012arXiv

A Common Proper Motion Stellar Companion to HAT-P-7

We report that HAT-P-7 has a common proper motion stellar companion. The companion is located at $\sim3.9$ arcsec to the east and estimated as an M5.5V dwarf based on its colors. We also confirm the presence of the third companion, which was first reported by Winn et al. (2009), based on long-term radial velocity measurements. We revisit the migration mechanism of HAT-P-7b given the presence of those companions, and propose sequential Kozai migration as a likely scenario in this system. This scenario may explain the reason for an outlier in the discussion of the spin-orbit alignment timescale for HAT-P-7b by Albrecht et al. (2012).

preprint2012arXiv

A Double Planetary System around the Evolved Intermediate-Mass Star HD 4732

We report the detection of a double planetary system orbiting around the evolved intermediate-mass star HD 4732 from precise Doppler measurements at Okayama Astrophysical Observatory (OAO) and Anglo-Australian Observatory (AAO). The star is a K0 subgiant with a mass of 1.7 M_sun and solar metallicity. The planetary system is composed of two giant planets with minimum mass of msini=2.4 M_J, orbital period of 360.2 d and 2732 d, and eccentricity of 0.13 and 0.23, respectively. Based on dynamical stability analysis for the system, we set the upper limit on the mass of the planets to be about 28 M_J (i>5 deg) in the case of coplanar prograde configuration.

preprint2012arXiv

Dynamics and Accretion of Planetesimals

We review the basic dynamics and accretion of planetesimals by showing N-body simulations. The orbits of planetesimals evolve through two-body gravitational relaxation: viscous stirring increases the random velocity and dynamical friction realizes the equiparation of the random energy. In the early stage of planetesimal accretion the growth mode of planetesimals is runaway growth where larger planetesimals grow faster than smaller ones. When a protoplanet (runaway-growing planetesimal) exceeds a critical mass the growth mode shifts to oligarchic growth where similar-sized protoplanets grow keeping a certain orbital separation. The final stage of terrestrial planet formation is collision among protoplanets known as giant impacts. We also summarize the dynamical effects of disk gas on planets and the core accretion model for formation of gas giants and discuss the diversity of planetary systems.

preprint2012arXiv

Migration rates of planets due to scattering of planetesimals

Planets migrate due to the recoil they experience from scattering solid (planetesimal) bodies. To first order, the torques exerted by the interior and exterior disks cancel, analogous to the cancellation of the torques from the gravitational interaction with the gas (type I migration). Assuming the dispersion-dominated regime and power-laws characterized by indices α and β for the surface density and eccentricity profiles, we calculate the net torque on the planet. We consider both distant encounters and close (orbit-crossing) encounters. We find that the close and distant encounter torques have opposite signs with respect to their α and β dependences; and that the torque is especially sensitive to the eccentricity gradient (β). Compared to type-I migration due to excitation of density waves, the planetesimal-driven migration rate is generally lower due to the lower surface density of solids in gas-rich disk, although this may be partially or fully offset when their eccentricity and inclination are small. Allowing for the feedback of the planet on the planetesimal disk through viscous stirring, we find that under certain conditions a self-regulated migration scenario emerges, in which the planet migrates at a steady pace that approaches the rate corresponding to the one-sided torque. If the local planetesimal disk mass to planet mass ratio is low, however, migration stalls. We quantify the boundaries separating the three migration regimes.

preprint2012arXiv

Minimum Dust Abundances for Planetesimal Formation via Secular Gravitational Instabilities

We estimate minimum dust abundances required for secular gravitational instability (SGI) to operate at the midplane dust layer of protoplanetary disks. For SGI to be a viable process, the growth time of the instability T_grow must be shorter than the radial drift time of the dust T_drift. The growth time depends on the turbulent diffusion parameter alpha, because the modes with short wavelengths are stabilized by turbulent diffusion. Assuming that turbulence is excited via the Kelvin-Helmholtz or streaming instabilities in the dust layer, and that its strength is controlled by the energy supply rate from dust accretion, we estimate the diffusion parameter and the growth time of the instability. The condition T_grow < T_drift requires that the dust abundance must be greater than a critical abundance Z_min, which is a function of the Toomre parameter Q_g and aspect ratio h_g / r of the gas disk. For a wide range of parameter space, the required dust abundance is less than 0.1. A slight increase in dust abundance opens a possible route for the dust to directly collapse to planetesimals.

preprint2012arXiv

N-body Simulations of Satellite Formation around Giant Planets: Origin of Orbital Configuration of the Galilean Moons

As the number of discovered extrasolar planets has been increasing, diversity of planetary systems requires studies of new formation scenarios. It is important to study satellite formation in circumplanetary disks, which is often viewed as analogous to formation of rocky planets in protoplanetary disks. We investigated satellite formation from satellitesimals around giant planets through N-body simulations that include gravitational interactions with a circumplanetary gas disk. Our main aim is to reproduce the observable properties of the Galilean satellites around Jupiter through numerical simulations, as previous N-body simulations have not explained the origin of the resonant configuration. We performed accretion simulations based on the work of Sasaki et al. (2010), in which an inner cavity is added to the model of Canup & Ward (2002, 2006). We found that several satellites are formed and captured in mutual mean motion resonances outside the disk inner edge and are stable after rapid disk gas dissipation, which explains the characteristics of the Galilean satellites. In addition, owing to the existence of the disk edge, a radial compositional gradient of the Galilean satellites can also be reproduced. An additional objective of this study is to discuss orbital properties of formed satellites for a wide range of conditions by considering large uncertainties in model parameters. Through numerical experiments and semianalytical arguments, we determined that if the inner edge of a disk is introduced, a Galilean-like configuration in which several satellites are captured into a 2:1 resonance outside the disk inner cavity is almost universal. In fact, such a configuration is produced even for a massive disk and rapid type I migration. This result implies the inevitability of a Galilean satellite formation in addition to providing theoretical predictions for extrasolar satellites.

preprint2012arXiv

Rapid Formation of Saturn after Jupiter Completion

We have investigated Saturn's core formation at a radial pressure maximum in a protoplanetary disk, which is created by gap opening by Jupiter. A core formed via planetesimal accretion induces the fragmentation of surrounding planetesimals, which generally inhibits further growth of the core by removal of the resulting fragments due to radial drift caused by gas drag. However, the emergence of the pressure maximum halts the drift of the fragments, while their orbital eccentricities and inclinations are efficiently damped by gas drag. As a result, the core of Saturn rapidly grows via accretion of the fragments near the pressure maximum. We have found that in the minimum-mass solar nebula, kilometer sized planetesimals can produce a core exceeding 10 Earth masses within two million years. Since Jupiter may not have undergone significant type II inward migration, it is likely that Jupiter's formation was completed when the local disk mass has already decayed to a value comparable to or less than Jovian mass. The expected rapid growth of Saturn's core on a timescale comparable to or shorter than observationally inferred disk lifetime enables Saturn to acquire the current amount of envelope gas before the disk gas is completely depleted. The high heat energy release rate onto the core surface due to the rapid accretion of the fragments delays onset of runaway gas accretion until the core mass becomes somewhat larger than that of Jupiter, which is consistent with the estimate based on interior modeling. Therefore, the rapid formation of Saturn induced by gap opening of Jupiter can account for the formation of multiple gas giants (Jupiter and Saturn) without significant inward migration and larger core mass of Saturn than that of Jupiter.

preprint2012arXiv

Substellar Companions to Seven Evolved Intermediate-Mass Stars

We report the detections of substellar companions orbiting around seven evolved intermediate-mass stars from precise Doppler measurements at Okayama Astrophysical Observatory. o UMa (G4 II-III) is a giant with a mass of 3.1 M_sun and hosts a planet with minimum mass of m_2sini=4.1 M_J in an orbit with a period P=1630 d and an eccentricity e=0.13. This is the first planet candidate (< 13 M_J) ever discovered around stars more massive than 3 M_sun. o CrB (K0 III) is a 2.1 M_sun giant and has a planet of m_2sini=1.5 M_J in a 187.8 d orbit with e=0.19. This is one of the least massive planets ever discovered around ~2 M_sun stars. HD 5608 (K0 IV) is an 1.6 M_sun subgiant hosting a planet of m_2sini=1.4 M_J in a 793 d orbit with e=0.19. The star also exhibits a linear velocity trend suggesting the existence of an outer, more massive companion. 75 Cet (G3 III:) is a 2.5 M_sun giant hosting a planet of m_2sini=3.0 M_J in a 692 d orbit with e=0.12. The star also shows possible additional periodicity of about 200 d and 1880 d with velocity amplitude of ~7--10 m/s, although these are not significant at this stage. nu Oph (K0 III) is a 3.0 M_sun giant and has two brown-dwarf companions of m_2sini= 24 M_J and 27 M_J, in orbits with P=530.3 d and 3190 d, and e=0.126 and 0.17, respectively, which were independently announced by Quirrenbach et al. (2011). The ratio of the periods is close to 1:6, suggesting that the companions are in mean motion resonance. We also independently confirmed planets around k CrB (K0 III-IV) and HD 210702 (K1 IV), which had been announced by Johnson et al. (2008) and Johnson et al. (2007a), respectively. All of the orbital parameters we obtained are consistent with the previous results.

preprint2012arXiv

The Orbital Stability of Planets Trapped in the First-Order Mean-Motion Resonances

Many extrasolar planetary systems containing multiple super-Earths have been discovered. N-body simulations taking into account standard type-I planetary migration suggest that protoplanets are captured into mean-motion resonant orbits near the inner disk edge at which the migration is halted. Previous N-body simulations suggested that orbital stability of the resonant systems depends on number of the captured planets. In the unstable case, through close scattering and merging between planets, non-resonant multiple systems are finally formed. In this paper, we investigate the critical number of the resonantly trapped planets beyond which orbital instability occurs after disk gas depletion. We find that when the total number of planets ($N$) is larger than the critical number ($N_{\rm crit}$), crossing time that is a timescale of initiation of the orbital instability is similar to non-resonant cases, while the orbital instability never occurs within the orbital calculation time ($10^8$ Kepler time) for $N\leq N_{\rm crit}$. Thus, the transition of crossing time across the critical number is drastic. When all the planets are trapped in 7:6 resonance of adjacent pairs, $N_{\rm crit} = 4$. We examine the dependence of the critical number of 4:3, 6:5 and 8:7 resonance by changing the orbital separation in mutual Hill radii and planetary mass. The critical number increases with increasing the orbital separation in mutual Hill radii with fixed planetary mass and increases with increasing planetary mass with fixed the orbital separation in mutual Hill radii. We also calculate the case of a system which is not composed of the same resonance. The sharp transition of the stability can be responsible for the diversity of multiple super-Earths (non-resonant or resonant), that is being revealed by $Kepler$ mission.

preprint2011arXiv

Evolution of Snow Line in Optically Thick Protoplanetary Disks: Effects of Water Ice Opacity and Dust Grain Size

Evolution of a snow line in an optically-thick protoplanetary disk is investigated with numerical simulations. The ice-condensing region in the disk is obtained by calculating the temperature and the density with the 1+1D approach. The snow line migrates as the mass accretion rate (\dot{M}) in the disk decreases with time. Calculations are carried out from an early phase with high disk accretion rates (\dot{M} \sim 10^{-7} M_sun/yr) to a later phase with low disk accretion rates (\dot{M} \sim 10^{-12} M_sun/yr) using the same numerical method. It is found that the snow line moves inward for \dot{M} > 10^{-10} M_sun/yr, while it gradually moves outward in the later evolution phase with \dot{M} < 10^{-10} M_sun/yr. In addition to the silicate opacity, the ice opacity is taken into consideration. In the inward migration phase, the additional ice opacity increases the distance of the snow line from the central star by a factor of 1.3 for dust grains < 10 micro-meter in size and 1.6 for > 100 micro-meter. It is inevitable that the snow line comes inside the Earth's orbit in the course of the disk evolution, if the alpha viscosity parameter is in a range 0.001-0.1, the dust-to-gas mass ratio is higher than a tenth of the solar abundance value, and the dust grains are smaller than 1 mm. The formation of water-devoid planetesimals in the terrestrial planet region seems to be difficult throughout the disk evolution, which imposes a new challenge to planet formation theory.

preprint2011arXiv

Induced Turbulence and the Density Structure of the Dust Layer in a Protoplanetary Disk

We study the turbulence induced in the dust layer of a protoplanetary disk based on the energetics of dust accretion due to gas drag. We estimate turbulence strength from the energy supplied by dust accretion, using the radial drift velocity of the dust particles in a laminar disk. Our estimate of the turbulence strength agrees with previous analytical and numerical research on the turbulence induced by Kelvin-Helmholtz and/or streaming instabilities for particles whose stopping time is less than the Keplerian time. For such small particles, the strongest turbulence is expected to occur when the dust-to-gas ratio of the disk is ~C_eff^(1/2) (h_g / r) ~ 10^(-2), where C_eff ~ 0.2 represents the energy supply efficiency to turbulence and h_g / r ~ 5 x 10^(-2) is the aspect ratio of the gas disk. The maximum viscosity parameter is alpha_max ~ C_eff T_s (h_g / r)^2 ~ 10^(-4) T_s, where T_s (<1) is the non-dimensional stopping time of the dust particles. Modification in the dust-to-gas ratio from the standard value, 10^(-2), by any process, results in weaker turbulence and a thinner dust layer, and consequently may accelerate the growth process of the dust particles.

preprint2011arXiv

On the Interaction between a Protoplanetary Disk and a Planet in an Eccentric Orbit: Application of Dynamical Friction

We present a new analytic approach to the disk-planet interaction that is especially useful for planets with eccentricity larger than the disk aspect ratio. We make use of the dynamical friction formula to calculate the force exerted on the planet by the disk, and the force is averaged over the period of the planet. The resulting migration and eccentricity damping timescale agrees very well with the previous works in which the planet eccentricity is moderately larger than the disk aspect ratio. The advantage of this approach is that it is possible to apply this formulation to arbitrary large eccentricity. We have found that the timescale of the orbital evolution depends largely on the adopted disk model in the case of highly eccentric planets. We discuss the possible implication of our results to the theory of planet formation.

preprint2011arXiv

Planetesimal Formation at the Boundary Between Steady Super/Sub-Keplerian Flow Created by Inhomogeneous Growth of Magnetorotational Instability

We have studied formation of planetesimals at a radial pressure bump in a protoplanetary disk created by radially inhomogeneous magnetorotational instability (MRI), through three-dimensional resistive MHD simulations including dust particles. In our previous papers, we showed that the inhomogeneous MRI developing in non-uniform structure of magnetic field or magnetic resistivity can transform the local gas flow in the disk to a quasi-steady state with local rigid rotation that is no more unstable against the MRI. Since the outer part of the rigid rotation is super-Keplerian flow, a quasi-static pressure bump is created and dust concentration is expected there. In this paper, we perform simulations of the same systems, adding dust particles that suffer gas drag and modulate gas flow via the back-reaction of the gas drag (dust drag). We use O(10^7) super-particles, each of which represents many dust particles with sizes of centimeter to meter. We have found that the dust drag suppresses turbulent motion to decrease the velocity dispersion of the dust particles while it broadens the dust concentrated regions to limit peaky dust concentration, compared with the simulation without the dust drag. We found that reduction in the velocity dispersion) is dominated over the suppression in particle concentration. For meter-size particles with the friction time ~1/Omega, where Omega is Keplerian frequency, the gravitational instability of the dust particles that may lead to planetesimal formation is expected. For such a situation, we further introduced the self-gravity of dust particles to the simulation to demonstrate that several gravitationally bound clumps are actually formed. Through analytical arguments, we found that the planetesimal formation from meter-sized dust particles can be possible at ~5AU, if dust spatial density is a few times larger than that in the minimum mass solar nebula.

preprint2011arXiv

The Gemini NICI Planet-Finding Campaign : Discovery of a Substellar L Dwarf Companion to the Nearby Young M Dwarf CD-35 2722

We present the discovery of a wide (67 AU) substellar companion to the nearby (21 pc) young solar-metallicity M1 dwarf CD-35 2722, a member of the ~100 Myr AB Doradus association. Two epochs of astrometry from the NICI Planet-Finding Campaign confirm that CD-35 2722 B is physically associated with the primary star. Near-IR spectra indicate a spectral type of L4\pm1 with a moderately low surface gravity, making it one of the coolest young companions found to date. The absorption lines and near-IR continuum shape of CD-35 2722 B agree especially well the dusty field L4.5 dwarf 2MASS J22244381-0158521, while the near-IR colors and absolute magnitudes match those of the 5 Myr old L4 planetary-mass companion, 1RXS J160929.1-210524 b. Overall, CD-35 2722 B appears to be an intermediate-age benchmark for L-dwarfs, with a less peaked H-band continuum than the youngest objects and near-IR absorption lines comparable to field objects. We fit Ames-Dusty model atmospheres to the near-IR spectra and find T=1700-1900 K and log(g) =4.5\pm0.5. The spectra also show that the radial velocities of components A and B agree to within \pm10 km/s, further confirming their physical association. Using the age and bolometric luminosity of CD-35 2722 B, we derive a mass of 31\pm8 Mjup from the Lyon/Dusty evolutionary models. Altogether, young late-M to mid-L type companions appear to be over-luminous for their near-IR spectral type compared to field objects, in contrast to the under-luminosity of young late-L and early-T dwarfs.

preprint2010arXiv

Detection of a low-eccentricity and super-massive planet to the subgiant HD 38801

We report the detection of a large mass planet orbiting around the K0 metal-rich subgiant HD38801 ($V=8.26$) by precise radial velocity (RV) measurements from the Subaru Telescope and the Keck Telescope. The star has a mass of $1.36M_{\odot}$ and metallicity of [Fe/H]= +0.26. The RV variations are consistent with a circular orbit with a period of 696.0 days and a velocity semiamplitude of 200.0\mps, which yield a minimum-mass for the companion of $10.7\mjup$ and semimajor axis of 1.71 AU. Such super-massive objects with very low-eccentricities and hundreds of days period are uncommon among the ensemble of known exoplanets.

preprint2010arXiv

Dust Concentration at the Boundary Between Steady Super/Sub-Keplerian Flow Created by Inhomogeneous Growth of MRI

How to create planetesimals from tiny dust particles in a proto-planetary disk before the dust particles spiral to the central star is one of the most challenging problems in the theory of planetary system formation. In our previous paper Kato et al. (2009), we have shown that a steady angular velocity profile that consists of both super and sub-Keplerian regions is created in the disk through non-uniform excitation of Magneto-Rotational Instability (MRI). Such non-uniform MRI excitation is reasonably expected in a part of disks with relatively low ionization degree. In this paper, we show through three-dimensional resistive MHD simulations with test particles that this radial structure of the angular velocity indeed leads to prevention of spiral-in of dust particles and furthermore to their accumulation at the boundary of super-Keplerian and sub-Keplerian regions. Treating dust particles as test particles, their motions under the influence of the non-uniform MRI through gas drag are simulated. In the most favorable cases (meter-size dust particles in the disk region with a relatively large fraction of MRI-stable region), we found that the dust concentration is peaked around the super/sub-Keplerian flow boundary and the peak dust density is 10,000 times as high as the initial value. The peak density is high enough for the subsequent gravitational instability to set in, suggesting a possible route to planetesimal formation via non-uniformly excited MRI in weakly ionized regions of a disk.

preprint2010arXiv

Eccentricity Trap: Trapping of Resonantly Interacting Planets near the Disk Inner Edge

Using orbital integration and analytical arguments, we have found a new mechanism (an "eccentricity trap") to halt type I migration of planets near the inner edge of a protoplanetary disk. Because asymmetric eccentricity damping due to disk-planet interaction on the innermost planet at the disk edge plays a crucial role in the trap, this mechanism requires continuous eccentricity excitation and hence works for a resonantly interacting convoy of planets. This trap is so strong that the edge torque exerted on the innermost planet can completely halt type I migrations of many outer planets through mutual resonant perturbations. Consequently, the convoy stays outside the disk edge, as a whole. We have derived semi-analytical formula for the condition for the eccentricity trap and predict how many planets are likely to be trapped. We found that several planets or more should be trapped by this mechanism in protoplanetary disks that have cavities. It can be responsible for the formation of non-resonant, multiple, close-in super-Earth systems extending beyond 0.1AU. Such systems are being revealed by radial velocity observations to be quite common around solar-type stars.

preprint2010arXiv

Origin of the Different Architectures of the Jovian and Saturnian Satellite Systems

The Jovian regular satellite system mainly consists of four Galilean satellites that have similar masses and are trapped in mutual mean motion resonances except for the outer satellite, Callisto. On the other hand, the Saturnian regular satellite system has only one big icy body, Titan, and a population of much smaller icy moons. We have investigated the origin of these major differences between the Jovian and Saturnian satellite systems by semi-analytically simulating the growth and orbital migration of proto-satellites in an accreting proto-satellite disk. We set up two different disk evolution/structure models that correspond to Jovian and Saturnian systems, by building upon previously developed models of an actively-supplied proto-satellite disk, the formation of gas giants, and observations of young stars. Our simulations extend previous models by including the (1) different termination timescales of gas infall onto the proto-satellite disk and (2) different evolution of a cavity in the disk, between the Jovian and Saturnian systems. We have performed Monte Carlo simulations and show that in the case of the Jovian systems, four to five similar-mass satellites are likely to remain trapped in mean motion resonances. This orbital configuration is formed by type I migration, temporal stopping of the migration near the disk inner edge, and quick truncation of gas infall caused by Jupiter opening a gap in the Solar nebula. The Saturnian systems tend to end up with one dominant body in the outer regions caused by the slower decay of gas infall associated with global depletion of the Solar nebula. The total mass and compositional zoning of the predicted Jovian and Saturnian satellite systems are consistent with the observed satellite systems.

preprint2010arXiv

Planet Migration

Planet migration is the process by which a planet's orbital radius changes in time. The main agent for causing gas giant planet migration is the gravitational interaction of the young planet with the gaseous disk from which it forms. We describe the migration rates resulting from these interactions based on a simple model for disk properties. These migration rates are higher than is reasonable for planet survival. We discuss some proposed models for which the migration rates are lower. There are major uncertainties in migration rates due to a lack of knowledge about the detailed physical properties of disks. We also describe some additional forms of migration.

preprint2010arXiv

Search for Outer Massive Bodies around Transiting Planetary Systems: Candidates of Faint Stellar Companions around HAT-P-7

We present results of direct imaging observations for HAT-P-7 taken with the Subaru HiCIAO and the Calar Alto AstraLux. Since the close-in transiting planet HAT-P-7b was reported to have a highly tilted orbit, massive bodies such as giant planets, brown dwarfs, or a binary star are expected to exist in the outer region of this system. We show that there are indeed two candidates for distant faint stellar companions around HAT-P-7. We discuss possible roles played by such companions on the orbital evolution of HAT-P-7b. We conclude that as there is a third body in the system as reported by Winn et al. (2009, ApJL, 763, L99), the Kozai migration is less likely while planet-planet scattering is possible.

preprint2010arXiv

Substellar Companions to Evolved Intermediate-Mass Stars: HD 145457 and HD 180314

We report the detections of two substellar companions orbiting around evolved intermediate-mass stars from precise Doppler measurements at Subaru Telescope and Okayama Astrophysical Observatory. HD 145457 is a K0 giant with a mass of 1.9 M_sun and has a planet of minimum mass m_2sini=2.9 M_J orbiting with period of P=176 d and eccentricity of e=0.11. HD 180314 is also a K0 giant with 2.6 M_sun and hosts a substellar companion of m_2sin i=22 M_J, which falls in brown-dwarf mass regime, in an orbit with P=396 d and e=0.26. HD 145457 b is one of the innermost planets and HD 180314 b is the seventh candidate of brown-dwarf-mass companion found around intermediate-mass evolved stars.

preprint2010arXiv

The Gemini NICI Planet-Finding Campaign

Our team is carrying out a multi-year observing program to directly image and characterize young extrasolar planets using the Near-Infrared Coronagraphic Imager (NICI) on the Gemini-South 8.1-meter telescope. NICI is the first instrument on a large telescope designed from the outset for high-contrast imaging, comprising a high-performance curvature adaptive optics system with a simultaneous dual-channel coronagraphic imager. Combined with state-of-the-art observing methods and data processing, NICI typically achieves ~2 magnitudes better contrast compared to previous ground-based or space-based programs, at separations inside of ~2 arcsec. In preparation for the Campaign, we carried out efforts to identify previously unrecognized young stars, to rigorously construct our observing strategy, and to optimize the combination of angular and spectral differential imaging. The Planet-Finding Campaign is in its second year, with first-epoch imaging of 174 stars already obtained out of a total sample of 300 stars. We describe the Campaign's goals, design, implementation, performance, and preliminary results. The NICI Campaign represents the largest and most sensitive imaging survey to date for massive (~1 Mjup) planets around other stars. Upon completion, the Campaign will establish the best measurements to date on the properties of young gas-giant planets at ~5-10 AU separations. Finally, Campaign discoveries will be well-suited to long-term orbital monitoring and detailed spectrophotometric followup with next-generation planet-finding instruments.

preprint2010arXiv

The Gemini NICI Planet-Finding Campaign: Discovery of a Close Substellar Companion to the Young Debris Disk Star PZ Tel

We report the discovery of a tight substellar companion to the young solar analog PZ Tel, a member of the Beta Pictoris moving group observed with high contrast adaptive optics imaging as part of the Gemini NICI Planet-Finding Campaign. The companion was detected at a projected separation of 16.4 +/- 1.0 AU (0.33 +/- 0.01") in April 2009. Second-epoch observations in May 2010 demonstrate that the companion is physically associated and shows significant orbital motion. Monte Carlo modeling constrains the orbit of PZ Tel B to eccentricities > 0.6. The near-IR colors of PZ Tel B indicate a spectral type of M7+/-2 and thus this object will be a new benchmark companion for studies of ultracool, low-gravity photospheres. Adopting an age of 12 +8 -4 Myr for the system, we estimate a mass of 36 +/- 6 Mjup based on the Lyon/DUSTY evolutionary models. PZ Tel B is one of few young substellar companions directly imaged at orbital separations similar to those of giant planets in our own solar system. Additionally, the primary star PZ Tel A shows a 70 um emission excess, evidence for a significant quantity of circumstellar dust that has not been disrupted by the orbital motion of the companion.

preprint2010arXiv

The Occurrence and Mass Distribution of Close-in Super-Earths, Neptunes, and Jupiters

The questions of how planets form and how common Earth-like planets are can be addressed by measuring the distribution of exoplanet masses and orbital periods. We report the occurrence rate of close-in planets (with orbital periods less than 50 days) based on precise Doppler measurements of 166 Sun-like stars. We measured increasing planet occurrence with decreasing planet mass (M). Extrapolation of a power law mass distribution fitted to our measurements, df/dlogM = 0.39M^-0.48, predicts that 23% of stars harbor a close-in Earth-mass planet (ranging from 0.5 to 2.0 Earth masses). Theoretical models of planet formation predict a deficit of planets in the domain from 5 to 30 Earth masses and with orbital periods less than 50 days. This region of parameter space is in fact well populated, implying that such models need substantial revision.

preprint2010arXiv

Thermal evolution and lifetime of intrinsic magnetic fields of Super Earths in habitable zones

We have numerically studied the thermal evolution of various-mass terrestrial planets in habitable zones, focusing on duration of dynamo activity to generate their intrinsic magnetic fields, which may be one of key factors in habitability on the planets. In particular, we are concerned with super-Earths, observations of which are rapidly developing. We calculated evolution of temperature distributions in planetary interior, using Vinet equations of state, Arrhenius-type formula for mantle viscosity, and the astrophysical mixing length theory for convective heat transfer modified for mantle convection. After calibrating the model with terrestrial planets in the Solar system, we apply it for 0.1--$10M_{\oplus}$ rocky planets with surface temperature of $300~\mbox{K}$ (in habitable zones) and the Earth-like compositions. With the criterion for heat flux at the CMB (core-mantle boundary), the lifetime of the magnetic fields is evaluated from the calculated thermal evolution. We found that the lifetime slowly increases with the planetary mass ($M_p$) independent of initial temperature gap at the core-mantle boundary ($ΔT_{\rm CMB}$) but beyond a critical value $M_{c,p}$ ($\sim O(1)M_{\oplus}$) it abruptly declines by the mantle viscosity enhancement due to the pressure effect. We derived $M_{c,p}$ as a function of $ΔT_{\rm CMB}$ and a rheological parameter (activation volume, $V^*$). Thus, the magnetic field lifetime of super-Earths with $M_p > M_{p,c}$ sensitively depends on $ΔT_{\rm CMB}$, which reflects planetary accretion, and $V^*$, which has uncertainty at very high pressure. More advanced high-pressure experiments and first-principle simulation as well as planetary accretion simulation are needed to discuss habitability of super-Earths.

preprint2009arXiv

A Substellar Companion in a 1.3 yr Nearly-circular Orbit of HD 16760

We report the detection of a substellar companion orbiting the G5 dwarf HD 16760 from the N2K sample. Precise Doppler measurements of the star from Subaru and Keck revealed a Keplerian velocity variation with a period of 466.47+-0.35 d, a semiamplitude of 407.71+-0.84 m/s, and an eccentricity of 0.084+-0.003. Adopting a stellar mass of 0.78+-0.05 M_Sun, we obtain a minimum mass for the companion of 13.13+-0.56 M_JUP, which is close to the planet/brown-dwarf transition, and the semimajor axis of 1.084+-0.023 AU. The nearly circular orbit despite the large mass and intermediate orbital period makes this companion unique among known substellar companions.

preprint2009arXiv

N-body simulations of planetary accretion around M dwarf stars

We have investigated planetary accretion from planetesimals in terrestrial planet regions inside the ice line around M dwarf stars through N-body simulations including tidal interactions with disk gas. Because of low luminosity of M dwarfs, habitable zones (HZs) are located in inner regions. In the close-in HZ, type-I migration and the orbital decay induced by eccentricity damping are efficient according to the high disk gas density in the small orbital radii. In the case of full efficiency of type-I migration predicted by the linear theory, we found that protoplanets that migrate to the vicinity of the host star undergo close scatterings and collisions, and 4 to 6 planets eventually remain in mutual mean motion resonances and their orbits have small eccentricities and they are stable both before and after disk gas decays. In the case of slow migration, the resonant capture is so efficient that densely-packed ~ 40 small protoplanets remain in mutual mean motion resonances. In this case, they start orbit crossing, after the disk gas decays and eccentricity damping due to tidal interaction with gas is no more effective. Through merging of the protoplanets, several planets in widely-separated non-resonant orbits with relatively large eccentricities are formed. Thus, the final orbital configurations of the terrestrial planets around M dwarfs sensitively depend on strength of type-I migration. We also found that large amount of water-ice is delivered by type-I migration from outer regions and final planets near the inner disk edge around M dwarfs are generally abundant in water-ice except for the innermost one that is shielded by the outer planets, unless type-I migration speed is reduced by a factor of more than 100 from that predicted by the linear theory.

preprint2008arXiv

Planetary Companions around Three Intermediate-Mass G and K Giants: 18 Del, xi Aql, and HD 81688

We report the detection of 3 new extrasolar planets from the precise Doppler survey of G and K giants at Okayama Astrophysical Observatory. The host stars, namely, 18 Del (G6 III), xi Aql (K0 III) and HD 81688 (K0 III-IV), are located at the clump region on the HR diagram with estimated masses of 2.1-2.3 M_solar. 18 Del b has a minimum mass of 10.3 M_Jup and resides in a nearly circular orbit with period of 993 days, which is the longest one ever discovered around evolved stars. xi Aql b and HD 81688 b have minimum masses of 2.8 and 2.7 M_Jup, and reside in nearly circular orbits with periods of 137 and 184 days, respectively, which are the shortest ones among planets around evolved stars. All of the substellar companions ever discovered around possible intermediate-mass (1.7-3.9 M_solar) clump giants have semimajor axes larger than 0.68 AU, suggesting the lack of short-period planets. Our numerical calculations suggest that Jupiter-mass planets within about 0.5 AU (even up to 1 AU depending on the metallicity and adopted models) around 2-3 M_solar stars could be engulfed by the central stars at the tip of RGB due to tidal torque from the central stars. Assuming that most of the clump giants are post-RGB stars, we can not distinguish whether the lack of short-period planets is primordial or due to engulfment by central stars. Deriving reliable mass and evolutionary status for evolved stars is highly required for further investigation of formation and evolution of planetary systems around intermediate-mass stars.

preprint2008arXiv

Planetary Companions to Evolved Intermediate-Mass Stars: 14 Andromedae, 81 Ceti, 6 Lyncis, and HD 167042

We report on the detection of four extrasolar planets orbiting evolved intermediate-mass stars from a precise Doppler survey of G and K giants at Okayama Astrophysical Observatory. All of the host stars are considered to be formerly early F-type or A-type dwarfs when they were on the main sequence. 14 And (K0 III) is a clump giant with a mass of 2.2 M_solar and has a planet of minimum mass m_2sin i=4.8 M_Jup in a nearly circular orbit with a 186 day period. This is one of the innermost planets around evolved intermediate-mass stars and such planets have only been discovered in clump giants. 81 Cet (G5 III) is a clump giant with 2.4 M_solar hosting a planet of m_2sin i=5.3 M_Jup in a 953 day orbit with an eccentricity of e=0.21. 6 Lyn (K0 IV) is a less evolved subgiant with 1.7 M_solar and has a planet of m_2sin i=2.4 M_Jup in a 899 day orbit with e=0.13. HD 167042 (K1 IV) is also a less evolved star with 1.5 M_solar hosting a planet of m_2sin i=1.6 M_Jup in a 418 day orbit with e=0.10. This planet was independently announced by Johnson et al. (2008, ApJ, 675, 784). All of the host stars have solar or sub-solar metallicity, which supports the lack of metal-rich tendency in planet-harboring giants in contrast to the case of dwarfs.

preprint2006arXiv

The Separation/Period Gap in the Distribution of Extrasolar Planets around Stars with Masses M >= 1.2 Msolar

The evidence for a shortage of exosolar planets with semimajor axes -1.1 <= log (a/AU) <= -0.2 is investigated. It is shown that this valley results from a gap in the radial distribution of planets, orbiting stars with masses M >= 1.2 Msolar (the high-mass sample, HMS). No underabundance is found for planets orbiting stars with smaller masses. The observational data also indicate that within the HMS population it is preferentially the more massive planets with M sin(i) >= 0.8 M_J that are missing. Monte-Carlo simulations of planet formation and migration are presented that reproduce the observed shortage of planets in the observed radius regime. A dependence on the disk depletion timescale tau_dep is found. The gap is more pronounced for tau_dep = 10^6 - 10^7 yrs than for tau_dep = 3*10^6 - 3*10^7 yrs. This might explain the observed trend with stellar mass if disks around stars with masses M* >= 1.2 Msolar have shorter depletion timescales than those around less massive stars. Possible reasons for such a dependence are a decrease of disk size and an increase of stellar EUV flux with stellar mass.