Search

Earth and other rocky bodies in the inner Solar System are significantly depleted in carbon compared to the Sun and interstellar medium (ISM) dust. Observations indicate that over half of carbon in the ISM and comets is in refractory forms, like amorphous hydrocarbons and complex organics, which can be building blocks of rocky bodies. While amorphous hydrocarbons are destroyed by photolysis and oxidation, radial transport of solid particles can limit carbon depletion, except when complex organics, which are less refractory, are the main carbon source. We aim to identify conditions for severe carbon depletion in the inner Solar System by introducing more realistic factors: differences in stickiness between icy and silicate particles, and high-temperature regions in the disk's upper optically-thin layer, which were not considered in previous studies. We perform a 3D Monte Carlo simulation of radial drift and turbulent diffusion in a steady accretion disk, incorporating ice evaporation/re-condensation, photolysis/oxidation of hydrocarbons in the upper layer, and pyrolysis of complex organics. Our results show that the carbon fraction drops by two orders of magnitude inside the snow line under two conditions: i) silicate particles are much less sticky than icy particles, leading to a rapid decline in icy pebble flux while silicates accumulate inside the snow line, and ii) high-temperature regions in the disk's upper layer stir silicate particles into UV-exposed areas. These conditions reproduce carbon depletion patterns consistent with observations and allow for diverse carbon fractions in rocky bodies. This diversity may explain the wide variation of metals in white dwarf photospheres and suggest different surface environments for rocky planets in habitable zones., Comment: 17 pages, 17 figures, accepted for publication in A&A

We present near-infrared high-dispersion spectroscopic observations of Europa using the Infrared Camera and Spectrograph (IRCS) onboard the Subaru Telescope, seeking direct evidence of water plumes on Europa and exploring spatial variations in plume activity. Using high spectral/spatial resolution and sensitivity of Subaru/IRCS, our observations have enabled a spatially resolved search for water plumes on Europa. Within our detection limits and time of observation, we found no evidence for the presence of water emission. For a rotation temperature of 50 K, we derived an upper limit on the H$_{2}$O abundance of 9.46$\times$10$^{19}$ - 5.92$\times$10$^{20}$ m$^{-2}$ in each divided slit area and 4.61$\times$10$^{19}$ m$^{-2}$ in the entire area covered by the slit. This upper limit lies below the inferred water abundance from previous UV observations by the Hubble Space Telescope (HST), while being less sensitive by a factor of three compared to the Keck telescope and by one order of magnitude or more than the James Webb Space Telescope (JWST) observations. Our results align with previous studies and demonstrate that using Subaru/IRCS is an effective strategy for searching for water plumes on Europa with high spatial resolution. Continued observations across different surface areas and orbital phases are essential to fully characterize Europa's plume activity and complement upcoming space missions., Comment: 8 pages, 7 figures. Accepted for publication in PASJ

Estimates of the frequency of planetary systems in the Milky Way are observationally limited by the low-mass planet regime. Nevertheless, substantial evidence for systems with undetectably low planetary masses now exist in the form of main-sequence stars which host debris discs, as well as metal-polluted white dwarfs. Further, low-mass sections of star formation regions impose upper bounds on protoplanetary disc masses, limiting the capacity for terrestrial or larger planets to form. Here, we use planetary population synthesis calculations to investigate the conditions that allow planetary systems to form only minor planets and smaller detritus. We simulate the accretional, collisional and migratory growth of $10^{17}$ kg embryonic seeds and then quantify which configurations with *entirely* sub-Earth-mass bodies ($\lesssim 10^{24}$ kg) survive. We find that substantial regions of the initial parameter space allow for sub-terrestrial configurations to form, with the success rate most closely tied to the initial dust mass. Total dust mass budgets of up to $10^2 M_{\oplus}$ within 10 au can be insufficiently high to form terrestrial or giant planets, resulting in systems with only minor planets. Consequently, the prevalence of planetary systems throughout the Milky Way might be higher than what is typically assumed, and minor planet-only systems may help inform the currently uncertain correspondence between planet-hosting white dwarfs and metal-polluted white dwarfs., Comment: Accepted for publication in MNRAS

Complex organic molecules (COMs) have been widely observed in molecular clouds and protostellar environments. One of the formation mechanisms of COMs is radical reactions on the icy grain surface driven by UV irradiation. While many experiments have reported that various COMs can be synthesized under such ice conditions, the majority of the reaction processes are unclear. Complementary numerical simulations are necessary to unveil the synthetic process behind the formation of COMs. In this study, we develop a chemical reaction simulation using a Monte Carlo method. To explore the complex reaction network of COM synthesis, the model was designed to eliminate the need to prepare reaction pathways and to keep computational costs low. With this simulation, we investigate the chemical reactions occurring on icy dust surfaces during and after UV irradiation, assuming a protoplanetary disk environment. We aim to reveal the types of organic molecules produced in a disk and the formation mechanisms of COMs, in particular, amino acids and sugars. The results show that photodissociation and subsequent radical-radical reactions cause random rearrangement of the covalent bonds in the initial molecules composed of methanol, formaldehyde, ammonia, and water. Consequently, highly complex molecules such as amino acids and sugars were produced in a wide range of the initial conditions. We found that the final abundances of amino acids and sugars have extremely similar dependence on the atomic ratios of the initial molecules, which peak at C/H~0.1-0.3 and O/H~0.3-0.5, although the amino acids abundance is usually more than ten times higher than that of sugars. To understand this dependence, a semi-analytical formula was derived. Additionally, parameter surveys have suggested that the decomposition reactions of amino acids and sugars undergo a rapid transition within the threshold of a given parameter., Comment: 20 pages, 20 figures, accepted for publication in Astronomy and Astrophysics (A&A)

Observations by the Voyager and Cassini spacecrafts have revealed various striking features of the gap structure in Saturn's ring, such as the density waves, sharp edge, and vertical wall structure. In order to explain these features in a single simulation, we perform a high-resolution (N~10^6-10^7) global full N-body simulation of gap formation by an embedded satellite considering gravitational interactions and inelastic collisions among all ring particles and the satellite, while these features have been mostly investigated separately with different theoretical approaches: the streamline models, 1D diffusion models, and local N-body simulation. As a first attempt of a series of papers, we here focus on the gap formation by separating satellite migration with fixing the satellite orbit in a Keplerian circular orbit. We reveal how the striking gap features - the density waves, sharp edge, and vertical wall structure - are simultaneously formed by an interplay of the satellite-ring and ring particle-particle interactions. In particular, we propose a new mechanism to quantitatively explain the creation of the vertical wall structure at the gap edge. Inelastic collisions between ring particles damp their eccentricity excited by the satellite's perturbations to enhance the surface density at the gap edge, making its sharp edges more pronounced. We find the eccentricity damping process inevitably raises the vertical wall structures the most effectively in the second epicycle waves. Particle-particle collisions generally convert their lateral epicyclic motion into vertical motion. Because the excited epicyclic motion is the greatest near the ring edge and the epicycle motions are aligned in the first waves, the conversion is the most efficient in the gap edge of the second waves and the wall height is scaled by the satellite Hill radius, which are consistent with the observations., Comment: 27 pages, 23 figures, accepted for publication in Icarus

One of the unique aspects of Earth is that it has a fractionally large Moon, which is thought to have formed from a Moon-forming disk generated by a giant impact. The Moon stabilizes the Earth's spin axis at least by several degrees and contributes to Earth's stable climate. Given that impacts are common during planet formation, exomoons, which are moons around planets in extrasolar systems, should be common as well, but no exomoon has been confirmed. Here we propose that an initially vapor-rich moon-forming disk is not capable of forming a large moon because growing moonlets, which are building blocks of a moon, experience strong gas drag and quickly fall toward the planet. Our impact simulations show that terrestrial and icy planets that are larger than $\sim 1.3-1.6 R_\oplus$ produce entirely vapor disks, which fail to form a large moon. This indicates that (1) our model supports the Moon-formation models that produce vapor-poor disks and (2) rocky and icy exoplanets whose radii are smaller than $\sim 1.6 R_\oplus$ are ideal candidates for hosting fractionally large exomoons., Comment: This paper was published in 2022

Dust particles in protoplanetary disks experience various chemical reactions under different physicochemical conditions through their accretion and diffusion, which results in the radial chemical gradient of dust. We performed three-dimensional Monte Carlo simulations to evaluate the dust trajectories and the progress of fictitious irreversible reactions, of which kinetics is expressed by the Johnson-Mehl-Avrami equation. The distribution of the highest temperature that each particle experiences before the degree of reaction exceeds a certain level shows the log-normal distribution, and its mode temperature was used as the effective reaction temperature. Semi-analytical prediction formulas of the effective reaction temperature and its dispersion were derived by comparing a reaction timescale with a diffusive transport timescale of dust as a function of the reaction parameters and the disk parameters. The formulas reproduce the numerical results of the effective reaction temperatures and their dispersions within 5.5 and 24 %, respectively, in a wide temperature range (200-1400 K). We applied the formulas for the crystallization of amorphous silicate dust and its oxygen isotope exchange with the H2O vapor based on the experimentally determined kinetics. For sub-micron sized amorphous forsterite dust, the predicted effective reaction temperature for the oxygen isotope exchange was lower than that of crystallization without overlap even considering their dispersions. This suggests that the amorphous silicate dust in the protosolar disk could exchange their oxygen isotopes efficiently with the 16O-poor H2O vapor, resulting in the distinct oxygen isotope compositions from the Sun., Comment: 28 pages, 7 figures, accepted for publication in ApJ

In this study, we revisit 32 planetary systems around evolved stars observed within the framework of the Okayama Planet Search Program and its collaborative framework of the EAPS-Net to search for additional companions and investigate the properties of stars and giant planets in multiple-planet systems. With our latest radial velocities obtained from Okayama Astrophysical Observatory (OAO), we confirm an additional giant planet in the wide orbit of 75 Cet system ($P_{\rm{c}} = 2051.62_{-40.47}^{+45.98}\ \rm{d}$, $M_{\rm{c}}\sin i=0.912_{-0.090}^{+0.088}\ M_{\rm{J}}$, and $a_{\rm{c}}=3.929_{-0.058}^{+0.052}\ \rm{au}$), along with five stars exhibiting long-term radial velocity accelerations, which indicates massive companions in the wide orbits. We have also found that the radial velocity variations of several planet-harboring stars may indicate additional planet candidates, stellar activities, or other understudied sources. These stars include $\epsilon$ Tau, 11 Com, 24 Boo, 41 Lyn, 14 And, HD 32518, and $\omega$ Ser. We further constrain the orbital configuration of the HD 5608, HD 14067, HD 120084, and HD 175679 systems by combining radial velocities with astrometry, as their host central stars exhibit significant astrometric accelerations. For other systems, we simply refine their orbital parameters. Moreover, our study indicates that the OPSP planet-harboring stars are more metal-poor compared to the currently known planet-harboring stars, and this is likely due to the $B-V$ color upper limit at 1.0 for star selection in the beginning of the survey. Finally, by investigating the less-massive giant planets ($< 5 M_{\rm{J}}$) around currently known planet-harboring evolved stars, we have found that metallicity positively correlates with the multiplicity and total planet mass of the system, which can be evidence for the core-accretion planet formation model., Comment: 49 figures, 4 tables, accepted by PASJ, RV data will be available online as supplementary after the publication

The suite of over 60 known planetary debris discs which orbit white dwarfs, along with detections of multiple minor planets in these systems, motivate investigations about the migration properties of planetesimals embedded within the discs. Here, we determine whether any of the migration regimes which are common in (pre-)main-sequence protoplanetary discs, debris discs and ring systems could be active and important in white dwarf discs. We investigate both dust-dominated and gas-dominated regions, and quantitatively demonstrate that Type I and Type II migration, as well as their particulate disc analogues, are too slow to be relevant in white dwarf discs. However, we find that the analogue of Type III migration for particulate discs may be rapid in the dusty regions of asteroid- or moon-generated ($>10^{18}$ kg) white dwarf discs, where a planetesimal exterior to its Roche radius may migrate across the entire disc within its lifetime. This result holds over a wide range of disc boundaries, both within and exterior to $1R_{\odot}$, and such that the probability of migration occurring increases with higher disc masses., Comment: Accepted for publication in MNRAS; corrected typos in Table 1 and References section

We investigate the spin state of a protoplanet during the pebble accretion influenced by the gas flow in the gravitational potential of the protoplanet and how it depends on the planetary mass, the headwind speed, the distance from the host star, and the pebble size. We perform nonisothermal three-dimensional hydrodynamical simulations in a local frame to obtain the gas flow around the planet. We then numerically integrate three-dimensional orbits of pebbles under the obtained gas flow. Finally, assuming uniform spatial distribution of incoming pebbles, we calculate net spin by summing up specific angular momentum that individual pebbles transfer to the protoplanet at impacts. We find that a protoplanet with the envelope acquires prograde net spin rotation regardless of the planetary mass, the pebble size, and the headwind speed of the gas. This is because accreting pebbles are dragged by the envelope that commonly has prograde rotation. As the planetary mass or orbital radius increases, the envelope is thicker and the prograde rotation is faster, resulting in faster net prograde spin. When the dimensionless thermal mass of the planet, $m = R_{\mathrm{Bondi}} / H$, where $R_{\mathrm{Bondi}}$ and $H$ are the Bondi radius and the disk gas scale height, is larger than a certain critical mass ($m \gtrsim 0.3$ at $0.1 \, \mathrm{au}$ or $m \gtrsim 0.1$ at $1 \, \mathrm{au}$), the spin rotation exceeds the breakup one. The predicted spin frequency reaches the breakup one at the planetary mass $m_{\mathrm{iso,rot}} \sim 0.1 \, (a / 1 \, \mathrm{au})^{-1/2}$ (where $a$ is the orbital radius), suggesting that the protoplanet cannot grow beyond $m_{\mathrm{iso,rot}}$. It is consistent with the Earth's current mass and could help the formation of the Moon by a giant impact on fast-spinning proto-Earth., Comment: 14 pages, 10 figures, Accepted for publication in Astronomy and Astrophysics (A&A)

A growing number of debris discs have been detected around metal-polluted white dwarfs. They are thought to be originated from tidally disrupted exoplanetary bodies and responsible for metal accretion onto host WDs. To explain (1) the observationally inferred accretion rate higher than that induced by Poynting-Robertson drag, $\dot{M}_{\rm PR}$, and (2) refractory-rich photosphere composition indicating the accretion of terrestrial rocky materials, previous studies proposed runaway accretion of silicate particles due to gas drag by the increasing silicate vapor produced by the sublimation of the particles. Because re-condensation of the vapor diffused beyond the sublimation line was neglected, we revisit this problem by one-dimensional advection/diffusion simulation that consistently incorporates silicate sublimation/condensation and back-reaction to particle drift due to gas drag in the solid-rich disc. We find that the silicate vapor density in the region overlapping the solid particles follows the saturating vapor pressure and that no runaway accretion occurs if the re-condensation is included. This always limits the accretion rate from mono-compositional silicate discs to $\dot{M}_{\rm PR}$ in the equilibrium state. Alternatively, by performing additional simulations that couple the volatile gas (e.g., water vapor), we demonstrate that the volatile gas enhances the silicate accretion to $>\dot{M}_{\rm PR}$ through gas drag. The refractory-rich accretion is simultaneously reproduced when the initial volatile fraction of disc is $\lesssim 10$ wt\% because of the suppression of volatile accretion due to the efficient back-reaction of solid to gas. The discs originating from C-type asteroid analogs might be a possible clue to the high-$\dot{M}$ puzzle., Comment: 20 pages, 14 figures, accepted for publication in MNRAS

We report the discovery of a triple-giant-planet system around an evolved star HD 184010 (HR 7421, HIP 96016). This discovery is based on observations from Okayama Planet Search Program, a precise radial velocity survey, undertaken at Okayama Astrophysical Observatory between 2004 April and 2021 June. The star is K0 type and located at beginning of the red-giant branch. It has a mass of $1.35_{-0.21}^{+0.19} M_{\odot}$, a radius of $4.86_{-0.49}^{+0.55} R_{\odot}$, and a surface gravity $\log g$ of $3.18_{-0.07}^{+0.08}$. The planetary system is composed of three giant planets in a compact configuration: The planets have minimum masses of $M_{\rm{b}}\sin i = 0.31_{-0.04}^{+0.03} M_{\rm{J}}$, $M_{\rm{c}}\sin i = 0.30_{-0.05}^{+0.04} M_{\rm{J}}$, and $M_{\rm{d}}\sin i = 0.45_{-0.06}^{+0.04} M_{\rm{J}}$, and orbital periods of $P_{\rm{b}}=286.6_{-0.7}^{+2.4}\ \rm{d}$, $P_{\rm{c}}=484.3_{-3.5}^{+5.5}\ \rm{d}$, and $P_{\rm{d}}=836.4_{-8.4}^{+8.4}\ \rm{d}$, respectively, which are derived from a triple Keplerian orbital fit to three sets of radial velocity data. The ratio of orbital periods are close to $P_{\rm{d}}:P_{\rm{c}}:P_{\rm{b}} \sim 21:12:7$, which means the period ratios between neighboring planets are both lower than $2:1$. The dynamical stability analysis reveals that the planets should have near-circular orbits. The system could remain stable over 1 Gyr, initialized from co-planar orbits, low eccentricities ($e=0.05$), and planet masses equal to the minimum mass derived from the best-fit circular orbit fitting. Besides, the planets are not likely in mean motion resonance. HD 184010 system is unique: it is the first system discovered to have a highly evolved star ($\log g < 3.5$ cgs) and more than two giant planets all with intermediate orbital periods ($10^2\ \rm{d} < P < 10^3\ \rm{d}$)., Comment: 20 pages, 5 figures, Published in PASJ

Recent high-spatial-resolution observations have revealed dust substructures in protoplanetary disks such as rings and gaps, which do not always correlate with gas. Because radial gas flow induced by low-mass, non-gas-gap-opening planets could affect the radial drift of dust, it potentially forms these dust substructures in disks. We investigate the potential of gas flow induced by low-mass planets to sculpt the rings and gaps in the dust profiles. We first perform three-dimensional hydrodynamical simulations, which resolve the local gas flow past a planet. We then calculate the trajectories of dust influenced by the planet-induced gas flow. Finally, we compute the steady-state dust surface density by incorporating the influences of the planet-induced gas flow into a one-dimensional dust advection-diffusion model. The outflow of the gas toward the outside of the planetary orbit inhibits the radial drift of dust, leading to dust accumulation (the dust ring). The outflow toward the inside of the planetary orbit enhances the inward drift of dust, causing dust depletion around the planetary orbit (the dust gap). Under weak turbulence ($\alpha_{\rm diff}\lesssim10^{-4}$, where $\alpha_{\rm diff}$ is the turbulence strength parameter), the gas flow induced by the planet with $\gtrsim1\,M_{\oplus}$ (Earth mass) generates the dust ring and gap in the distribution of small dust grains ($\lesssim1$ cm) with the radial extent of $\sim1\text{--}10$ times gas scale height around the planetary orbit without creating a gas gap and pressure bump. The gas flow induced by low-mass, non-gas-gap-opening planets can be considered a possible origin of the observed dust substructures in disks. Our results may be helpful to explain the disks whose dust substructures were found not to correlate with those of the gas., Comment: 25 pages, 20 figures, Accepted for publication in Astronomy and Astrophysics (A&A)

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., Comment: 14 pages, accepted for publication in Astronomy and Astrophysics

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, $\tau_{\rm s} \simeq 0.01-0.1$, $\alpha_{\rm mid} \lesssim 10^{-4}$, and $\alpha_{\rm acc} \simeq 10^{-3}-10^{-2}$ at $r \lesssim 10$ au, where $r$, $\tau_{\rm s}$, $\alpha_{\rm mid}$, and $\alpha_{\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 $\alpha$-disk (i.e., gas surface density), respectively., Comment: 13 pages, 7 Figures, accepted for publication in Astronomy & Astrophysics (A&A)

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 ($\beta$), and dust temperature. The derived parameters are consistent with previous works, and the enhancement of $\beta$ within the intensity gaps is also confirmed, supporting a deficit of millimeter-sized grains within the gaps., Comment: 17pages, 12 figures, Accepted for publication in The Astrophysical Journal

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., Comment: Accepted for publication in the Astrophysical Journal

We present molecular line observations of 13CO and C18O J=3-2, CN N = 3 - 2, and CS J=7-6 lines in the protoplanetary disk around TW Hya at a high spatial resolution of ~9 au (angular resolution of 0.15''), using the Atacama Large Millimeter/Submillimeter Array. A possible gas gap is found in the deprojected radial intensity profile of the integrated C18O line around a disk radius of ~58 au, slightly beyond the location of the au-scale dust clump at ~52 au, which resembles predictions from hydrodynamic simulations of planet-disk interaction. In addition, we construct models for the physical and chemical structure of the TW Hya disk, taking account of the dust surface density profile obtained from high spatial resolution dust continuum observations. As a result, the observed flat radial profile of the CN line intensities is reproduced due to a high dust-to-gas surface density ratio inside ~20 au. Meanwhile, the CO isotopologue line intensities trace high temperature gas and increase rapidly inside a disk radius of ~30 au. A model with either CO gas depletion or depletion of gas-phase oxygen elemental abundance is required to reproduce the relatively weak CO isotopologue line intensities observed in the outer disk, consistent with previous atomic and molecular line observations towards the TW Hya disk. {Further observations of line emission of carbon-bearing species, such as atomic carbon and HCN, with high spatial resolution would help to better constrain the distribution of elemental carbon abundance in the disk gas., Comment: 22 pages, 7 figures, Accepted for publication in ApJ

Aims. The connection between initial disc conditions and final orbital and physical properties of planets is not well-understood. In this paper, we numerically study the formation of planetary systems via pebble accretion and investigate the effects of disc properties such as masses, dissipation timescales, and metallicities on planet formation outcomes. Methods. We improved the N-body code SyMBA that was modified by taking account of new planet-disc interaction models and type II migration. We adopted the 'two-alpha disc' model to mimic the effects of both the standard disc turbulence and the mass accretion driven by the magnetic disc wind. Results. We successfully reproduced the overall distribution trends of semi-major axes, eccentricities, and planetary masses of extrasolar giant planets. We find that, when planet formation happens fast enough, giant planets are fully grown (Jupiter mass or higher) and are distributed widely across the disc. On the other hand, when planet formation is limited by the disc's dissipation, discs generally form low-mass cold Jupiters (CJs). Our simulations also naturally explain why hot Jupiters (HJs) tend to be alone and how the observed eccentricity-metallicity trends arise. The low-metallicity discs tend to form nearly circular and coplanar HJs in situ, because planet formation is slower than high-metallicity discs, and thus protoplanetary cores migrate significantly before gas accretion. The high-metallicity discs, on the other hand, generate HJs in situ or via tidal circularisation of eccentric orbits. Both pathways usually involve dynamical instabilities, and thus HJs tend to have broader eccentricity and inclination distributions. When giant planets with very wide orbits ('super-cold Jupiters') are formed, we find that they often belong to metal-rich stars, have eccentric orbits, and tend to have (~80%) companions interior to their orbits., Comment: 29 pages, 16 figures, 2 tables, A&A (in press)

We aim to understand the effect of stellar evolution on the evolution of protoplanetary disks. We focus in particular on the disk evolution around intermediate-mass (IM) stars, which evolve more rapidly than low-mass ones. We numerically solve the long-term evolution of disks around 0.5-5 solar-mass stars considering viscous accretion and photoevaporation (PE) driven by stellar far-ultraviolet (FUV), extreme-ultraviolet (EUV), and X-ray emission. We also take stellar evolution into account and consider the time evolution of the PE rate. We find that the FUV, EUV, and X-ray luminosities of IM stars evolve by orders of magnitude within a few Myr along with the time evolution of stellar structure, stellar effective temperature, or accretion rate. Therefore, the PE rate also evolves with time by orders of magnitude, and we conclude that stellar evolution is crucial for the disk evolution around IM stars., Comment: 21 pages, 15 figures, published in ApJ. Evolutionary models of young 0.5-5 Msun stars are available at https://zenodo.org/record/4595807

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., Comment: 7 pages, 3 Figures, 1 Figure in Appendix; Accepted for publication in Astronomy & Astrophysics Letters (A&A Letters)

Uranus and Neptune are the archetypes of "ice giants", a class of planets that may be among the most common in the Galaxy. They hold the keys to understand the atmospheric dynamics and structure of planets with hydrogen atmospheres inside and outside the solar system; however, they are also the last unexplored planets of the Solar System. Their atmospheres are active and storms are believed to be fueled by methane condensation which is both extremely abundant and occurs at low optical depth. This means that mapping temperature and methane abundance as a function of position and depth will inform us on how convection organizes in an atmosphere with no surface and condensates that are heavier than the surrounding air, a general feature of giant planets. Owing to the spatial and temporal variability of these atmospheres, an orbiter is required. A probe would provide a reference atmospheric profile to lift ambiguities inherent to remote observations. It would also measure the abundances of noble gases which can be used to reconstruct the history of planet formation in the Solar System. Finally, mapping the planets' gravity and magnetic fields will be essential to constrain their global composition, atmospheric dynamics, structure and evolution. An exploration of Uranus or Neptune will be essential to understand these planets and will also be key to constrain and analyze data obtained at Jupiter, Saturn, and for numerous exoplanets with hydrogen atmospheres., Comment: arXiv admin note: substantial text overlap with arXiv:1908.02092

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 ($\alpha_{\rm Dr}$ and $\alpha_{\rm Dz}$) and for the gas accretion to the star ($\alpha_{\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}-\alpha_{\rm acc}-\alpha_{\rm Dz}(=\alpha_{\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 $\alpha_{\rm acc}=10^{-3}$ than for $\alpha_{\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 $\alpha_{\rm Dr}/\alpha_{\rm acc} \ll 1$, while that of pebbles outside the snow line is favored for $\alpha_{\rm Dr}/\alpha_{\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., Comment: 18 pages, 6 Figures, 3 Figures in Appendix, accepted for publication in Astronomy & Astrophysics (A&A)

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 $\alpha$ 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 $> [(\alpha_{Dz}/\alpha_{acc})/3 \times 10^{-2}]^{1/2}$ where $\alpha_{Dz}$ and $\alpha_{acc}$ are the $\alpha$ parameters for vertical turbulent diffusion and disk gas accretion., Comment: 15 pages, accepted for publication in Astronomy & Astrophysics

We investigate the impact of pre-main sequence stellar luminosity evolution on the thermal and chemical properties of disc midplanes. We create template disc models exemplifying initial conditions for giant planet formation for a variety of stellar masses and ages. These models include the 2D physical structure of gas as well as 1D chemical structure in the disc midplane. The disc temperature profiles are calculated using fully physically consistent radiative transfer models for stars between 0.5 and 3 Msun and ages up to 10 Myr. The resulting temperature profiles are used to determine how the chemical conditions in the mid-plane change over time. We therefore obtain gas and ice-phase abundances of the main carbon and oxygen carrier species. While the temperature profiles produced are not markedly different for the stars of different masses at early stages (<1 Myr), they start to diverge significantly beyond 2 Myr. Discs around stars with mass >1.5 Msun become warmer over time as the stellar luminosity increases, whereas low-mass stars decrease in luminosity leading to cooler discs. This has an observable effect on the location of the CO snowline, which is located >200 au in most models for a 3 Msun star, but is always within 80 au for 0.5 Msun star. The chemical compositions calculated show that a well defined stellar mass and age range exists in which high C/O gas giants can form. In the case of the exoplanet HR8799b, our models show it must have formed before the star was 1 Myr old.

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., Comment: (6 pages, 4 figures, Accepted in A&A)

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., Comment: 10 pages, 1 figure, accepted for publication in MNRAS; typos corrected, the reference list was completed

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., Comment: Published in Nature Astronomy, 18 pages, 2 figures

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., Comment: 25 pages, 11 figures, 6 tables, accepted for publication in ApJ

It is important to clarify where and when rocky and icy planetesimals are formed in a viscously evolving disk. We wish to understand how local runaway pile-up of solids occurs inside or outside the snow line. We assume an icy pebble contains micron-sized silicate grains that are uniformly mixed with ice and are released during the ice sublimation. Using a local one-dimensional code, we solve the radial drift and the turbulent diffusion of solids and the water vapor, taking account of their sublimation/condensation around the snow line. We systematically investigate effects of back-reactions of the solids to gas on the radial drift and diffusion of solids, scale height evolution of the released silicate particles, and possible difference in effective viscous parameters between that for turbulent diffusion ($\alpha_{\rm tur}$) and that for the gas accretion rate onto the central star ($\alpha_{\rm acc}$). We study the dependence on the ratio of the solid mass flux to the gas ($F_{\rm p/g}$). We show that the favorable locations for the pile-up of silicate grains and icy pebbles are the regions in the proximity of the water snow line inside and outside it, respectively. We found that runaway pile-ups occur when both the back-reactions for radial drift and diffusion are included. In the case with only the back-reaction for the radial drift, no runaway pile-up is found except for extremely high pebble flux, while the condition of streaming instability can be satisfied for relatively large $F_{\rm p/g}$ as found in the past literatures. If the back-reactions for radial diffusion is considered, the runaway pile-up occurs for reasonable value of pebble flux. The runaway pile-up of silicate grains that would lead to formation of rocky planetesimals occurs for $\alpha_{\rm tur} \ll \alpha_{\rm acc}$, while the runaway pile-up of icy pebbles is favored for $\alpha_{\rm tur} \sim \alpha_{\rm acc}$., Comment: Accepted for publication in Astronomy & Astrophysics (A&A); 13pages, 8 figures

Planets in the "Habitable Zones" around M-type stars are important targets for characterization in future observations. Due to tidal-locking in synchronous spin-orbit rotations, the planets tend to have a hot dayside and a cold nightside. On the cold nightside, water vapor transferred from the dayside can be frozen in ("cold trap") or the major atmospheric constituent could also condense ("atmospheric collapse") if the atmosphere is so thin that the heat re-distribution is not efficient, in the case of a single M-type star. Motivated by the abundance of binary star systems, we investigate the effects of irradiation from a G-type companion star on the climate of a tidally locked planet around an M-type star using the 2D energy balance model. We find that the irradiation from the G-type star is more effective at warming up the nightside of the planet than the dayside. This contributes to the prevention of the irreversible trapping of water and atmosphere on the cold nightside, broadening the parameter space where tidally locked planets can maintain surface liquid water. Tidally locked ocean planets with < ~ 0.3 bar atmospheres or land planets with < ~ 3 bar atmospheres can realize temperate climate with surface liquid water when they are also irradiated by a companion star with a separation of 1 - 4 au. We also demonstrate that planets with given properties can be in the Earth-like temperate climate regime or in a completely frozen state under the same total irradiation., Comment: 28 pages, 15 figures, accepted for publication in The Astrophysical Journal

We report the detection of an excess in dust continuum emission at 233~GHz (1.3~mm in wavelength) in the protoplanetary disk around TW~Hya revealed through high-sensitivity observations at $\sim$3~au resolution with the Atacama Large Millimeter/submillimeter Array (ALMA). The sensitivity of the 233~GHz image has been improved by a factor of 3 with regard to that of our previous cycle 3 observations. The overall structure is mostly axisymmetric, and there are apparent gaps at 25 and 41 au as previously reported. The most remarkable new finding is a few au-scale excess emission in the south-west part of the protoplanetary disk. The excess emission is located at 52 au from the disk center and is 1.5 times brighter than the surrounding protoplanetary disk at a significance of 12$\sigma$. We performed a visibility fitting to the extracted emission after subtracting the axisymmetric protoplanetary disk emission and found that the inferred size and the total flux density of the excess emission are 4.4$\times$1.0~au and 250~$\mu$Jy, respectively. The dust mass of the excess emission corresponds to 0.03~$M_\oplus$ if a dust temperature of 18~K is assumed. Since the excess emission can also be marginally identified in the Band 7 image at almost the same position, the feature is unlikely to be a background source. The excess emission can be explained by a dust clump accumulated in a small elongated vortex or a massive circumplanetary disk around a Neptune mass forming-planet., Comment: 9pages, 5 figures, Accepted for publication in The Astrophysical Journal Letters

Understanding the fate of planetary systems through white dwarfs which accrete debris crucially relies on tracing the orbital and physical properties of exo-asteroids during the giant branch phase of stellar evolution. Giant branch luminosities exceed the Sun's by over three orders of magnitude, leading to significantly enhanced Yarkovsky and YORP effects on minor planets. Here, we place bounds on Yarkovsky-induced differential migration between asteroids and planets during giant branch mass loss by modelling one exo-Neptune with inner and outer exo-Kuiper belts. In our bounding models, the asteroids move too quickly past the planet to be diverted from their eventual fate, which can range from: (i) populating the outer regions of systems out to 10^4-10^5 au, (ii) being engulfed within the host star, or (iii) experiencing Yarkovsky-induced orbital inclination flipping without any Yarkovsky-induced semimajor axis drift. In these violent limiting cases, temporary resonant trapping of asteroids with radii of under about 10 km by the planet is insignificant, and capture within the planet's Hill sphere requires fine-tuned dissipation. The wide variety of outcomes presented here demonstrates the need to employ sophisticated structure and radiative exo-asteroid models in future studies. Determining where metal-polluting asteroids reside around a white dwarf depends on understanding extreme Yarkovsky physics., Comment: Accepted for publication in MNRAS

The three-dimensional structure of the gas flow around a planet is thought to influence the accretion of both gas and solid materials. In particular, the outflow in the mid-plane region may prevent the accretion of the solid materials and delay the formation of super-Earths' cores. However, it is not yet understood how the nature of the flow field and outflow speed change as a function of the planetary mass. In this study, we investigate the dependence of gas flow around a planet embedded in a protoplanetary disc on the planetary mass. Assuming an isothermal, inviscid gas disc, we perform three-dimensional hydrodynamical simulations on the spherical polar grid, which has a planet located at its centre. We find that gas enters the Bondi or Hill sphere at high latitudes and exits through the mid-plane region of the disc regardless of the assumed dimensionless planetary mass $m=R_{\rm Bondi}/H$, where $R_{\rm Bondi}$ and $H$ are the Bondi radius of the planet and disc scale height, respectively. The altitude from where gas predominantly enters the envelope varies with the planetary mass. The outflow speed can be expressed as $|u_{\rm out}|=\sqrt{3/2}mc_{\rm s}$ $(R_{\rm Bondi}\leq R_{\rm Hill})$ or $|u_{\rm out}|=\sqrt{3/2}(m/3)^{1/3} c_{\rm s}$ ($R_{\rm Bondi}\geq R_{\rm Hill}$), where $c_{\rm s}$ is the isothermal sound speed and $R_{\rm Hill}$ is the Hill radius. The outflow around a planet may reduce the accretion of dust and pebbles onto the planet when $m\gtrsim\sqrt{\rm St}$, where St is the Stokes number. Our results suggest that the flow around proto-cores of super-Earths may delay their growth and, consequently, help them to avoid runaway gas accretion within the lifetime of the gas disc., Comment: 16 pages, 14 figures, Accepted for publication in Astronomy and Astrophysics (A&A)

The Earth's ocean mass is only 2.3 x 10^{-4} of the whole planet mass. Even including water in the interior, it would be at most 10^{-3}-10^{-2}. Ancient Mars may have had a similar or slightly smaller water fraction. It is important to clarify the water delivery mechanism to rocky planets in habitable zones in exoplanetary systems, as well as that to the Earth and Mars. Here, we consider water delivery to planets by icy pebbles after the snowline inwardly passes the planetary orbits and derive the water mass fraction (f_{water}) of the final planet as a function of disk parameters and discuss the parameters that reproduce f_{water} comparable to that inferred for the Earth and ancient Mars. We calculate the growth of icy pebbles and their radial drift with a 1D model, and accretion of icy pebbles onto planets, by simultaneously solving the snowline migration and the disk dissipation, to evaluate f_{water} of the planets. We find that f_{water} is regulated by the total mass (M_{res}) of icy dust materials preserved in the outer disk regions at the timing (t = t_{snow}) of the snowline passage of the planetary orbit. Because M_{res} decays rapidly after the pebble formation front reaches the disk outer edge (at t = t_{pff}), f_{water} is sensitive to the ratio t_{snow}/t_{pff}, which is determined by the disk parameters. We find t_{snow}/t_{pff} < 10 or > 10 is important. Deriving an analytical formula for f_{water} that reproduces the numerical results, we find that f_{water} of a rocky planet near 1 au is ~ 10^{-4}-10^{-2}, in the disks with initial disk size ~ 30-50 au and the initial disk mass accretion rate ~ (10^{-8}-10^{-7}) M_sun/y. Because these disks may be median or slightly compact/massive disks, the water fraction of rocky planets in habitable zones may be often similar to that of the Earth, if the icy pebble accretion is responsible for the water delivery., Comment: published in A&A in 2019; a corrected version from the published version; in the previous replacement, some figures were not shown. This is the version with complete figures

We compare the planet-to-star mass-ratio distribution measured by gravitational microlensing to core accretion theory predictions from population synthesis models. The core accretion theory's runaway gas accretion process predicts a dearth of intermediate-mass giant planets that is not seen in the microlensing results. In particular, the models predict $\sim10\,\times$ fewer planets at mass ratios of $10^{-4} \leq q \leq 4 \times 10^{-4}$ than inferred from microlensing observations. This tension implies that gas giant formation may involve processes that have hitherto been overlooked by existing core accretion models or that the planet-forming environment varies considerably as a function of host-star mass. Variation from the usual assumptions for the protoplanetary disk viscosity and thickness could reduce this discrepancy, but such changes might conflict with microlensing results at larger or smaller mass ratios, or with other observations. The resolution of this discrepancy may have important implications for planetary habitability because it has been suggested that the runaway gas accretion process may have triggered the delivery of water to our inner solar system. So, an understanding of giant planet formation may help us to determine the occurrence rate of habitable planets., Comment: 12 pages, 2 figures, 1 table, accepted for publication in ApJL

Planetary migration is a major challenge for planet formation theories. The speed of Type I migration is proportional to the mass of a protoplanet, while the final decade of growth of a pebble-accreting planetary core takes place at a rate that scales with the mass to the two-thirds power. This results in planetary growth tracks (i.e., the evolution of a protoplanet's mass versus its distance from the star) that become increasingly horizontal (migration-dominated) with rising mass of the protoplanet. It has been shown recently that the migration torque on a protoplanet is reduced proportional to the relative height of the gas gap carved by the growing planet. Here we show from 1-D simulations of planet-disc interaction that the mass at which a planet carves a 50% gap is approximately 2.3 times the pebble isolation mass. Our measurements of the pebble isolation mass from 1-D simulations match published 3-D results relatively well, except at very low viscosities where the 3-D pebble isolation mass is significantly higher, possibly due to gap edge instabilities not captured in 1-D. The pebble isolation mass demarks the transition from pebble accretion to gas accretion. Gas accretion to form gas-giant planets therefore takes place over a few astronomical units of migration after reaching first the pebble isolation mass and, shortly after, the 50% gap mass. Our results demonstrate how planetary growth can outperform migration, both during core accretion and during gas accretion, even when the Stokes number of the pebbles is small, St~0.01, and the pebble-to-gas flux ratio in the protoplanetary disc is in the nominal range of 0.01-0.02. We find that planetary growth is very rapid in the first million years of the protoplanetary disc and that the probability for forming gas-giant planets increases with the initial size of the protoplanetary disc and with decreasing turbulent diffusion., Comment: Accepted for publication in Astronomy & Astrophysics

The formation and orbital evolution of Saturn's inner mid-sized moons are still debated. The most puzzling aspects are 1) how the Tethys-Dione pair and the Mimas-Enceladus pair passed through their strong 3:2 mean-motion resonances during the tidal orbital evolution, and 2) the current strong heat flow from Enceladus, which is a few orders of magnitude higher than the tidal energy dissipation caused by the present orbital eccentricity of Enceladus. Here we perform N-body simulations of the moons' orbital evolution from various initial conditions -- assuming that the moons were formed from Saturn's hypothetical massive ring -- and investigate possible paths to solve the above difficulties. If the influence of the rings is neglected, we find that the Tethys-Dione pair cannot avoid becoming trapped in the mean-motion resonances as they recede from Saturn, and that the Tethys-Enceladus pair cannot avoid collisions after the resonance trapping, in case Saturn's quality factor is smaller than 15,000. These findings are inconsistent with the current orbital configuration. However, taking into account both the eccentricity excitation and the orbital expansion caused by the ring torque, we find that these resonance captures are avoided. With the high eccentricity pumped up by the torque, Enceladus passes through the resonances with Tethys, and the Dione-Tethys pair passes through their 2:1 and possibly the 3:2 resonance as well. After Enceladus resides beyond the 2:1 resonance with the outer ring edge, the eccentricity can be tidally damped. While this is a promising path of evolution, in most runs, Enceladus collides with Tethys by the excited eccentricity. There is a hint that a ring mass decrease could avoid the collision. The tidal heat due to the eccentricity excitation by the ring torque is stored in the moons and slowly radiated away. It may account for the current high heat flow., Comment: 37 pages, 7 figures, accepted for publication in Icarus

We want to investigate how planet formation is imprinted on stellar surface composition using up-to-date stellar evolution models. We simulate the evolution of pre-main-sequence stars as a function of the efficiency of heat injection during accretion, the deuterium mass fraction, and the stellar mass. For simplicity, we assume that planet formation leads to the late accretion of zero-metallicity gas, diluting the surface stellar composition as a function of the mass of the stellar outer convective zone. We adopt $150\,{\mathrm{M}_\oplus}(M_\star/\mathrm{M}_\odot)(Z/\mathrm{Z}_\odot)$ as an uncertain but plausible estimate of the mass of heavy elements that is not accreted by stars with giant planets, including our Sun. By combining our stellar evolution models to these estimates, we evaluate the consequences of planet formation on stellar surface composition. We show that after the first $\sim0.1$ Myr, the evolution of the convective zone follows classical evolutionary tracks within a factor of two in age. We find that planet formation should lead to a scatter in stellar surface composition that is larger for high-mass stars than for low-mass stars. We predict a spread in [Fe/H] of approximately $0.02$ dex for stars with $T_\mathrm{eff}\sim 5500\,$K, marginally compatible with differences in metallicities observed in some binary stars with planets. Stars with $T_\mathrm{eff}\geq 7000\,$K may show much larger [Fe/H] deficits, by 0.6 dex or more, compatible with the existence of refractory-poor $\lambda$ Boo stars. We also find that planet formation may explain the lack of refractory elements seen in the Sun as compared to solar twins, but only if the ice-to-rock ratio in the solar-system planets is less than $\approx0.4$ and planet formation began less than $\approx1.3$ Myr after the beginning of the formation of the Sun. (abbreviated), Comment: Accepted for publicatoin in A&A. 18 pages, 14 figures

The mass and semimajor axis distribution of gas giants in exoplanetary systems obtained by radial velocity surveys shows that super-jupiter-mass planets are piled up at > 1 au, while jupiter/sub-jupiter-mass planets are broadly distributed from ~0.03 au to beyond 1 au. This feature has not been explained by theoretical predictions. In order to reconcile this inconsistency, we investigate evolution of gas giants with a new type II migration formula by Kanagawa et al. (2018), by comparing the migration, growth timescales of gas giants, and disk lifetime and by population synthesis simulation. While the classical migration model assumes that a gas giant opens up a clear gap in the protoplanetary disk and the planet migration tied to the disk gas accretion, recent high-resolution simulations show that the migration of gap-opening planets is decoupled from the disk gas accretion and Kanagawa et al. (2018) proposed that type II migration speed is no other than type I migration speed with the reduced disk gas surface density in the gap. We show that with this new formula, type II migration is significantly reduced for super-jupiter-mass planets, if the disk accretion is driven by the disk wind as suggested by recent MHD simulations. Population synthesis simulations show that super-jupiter-mass planets remain at > 1 au without any additional ingredient such as disk photoevaporation. Therefore, the mystery of the pile-up of gas giants at > 1 au would be theoretically solved, if the new formula is confirmed and wind-driven disk accretion dominates., Comment: Accepted for publication in ApJ; the typos in Eq.(31) were corrected

Dynamical models of planet formation coupled with cosmochemical data from martian meteorites show that Mars' isotopic composition is distinct from that of Earth. Reconciliation of formation models with meteorite data require that Mars grew further from the Sun than its present position. Here, we evaluate this compositional difference in more detail by comparing output from two $N$-body planet formation models. The first of these planet formation models simulates what is termed the "Classical" case wherein Jupiter and Saturn are kept in their current orbits. We compare these results with another model based on the "Grand Tack", in which Jupiter and Saturn migrate through the primordial asteroid belt. Our estimate of the average fraction of chondrite assembled into Earth and Mars assumes that the initial solid disk consists of only sources of enstatite chondrite composition in the inner region, and ordinary chondrite in the outer region. Results of these analyses show that both models tend to yield Earth and Mars analogues whose accretion zones overlap. The Classical case fares better in forming Mars with its documented composition (29% to 68% enstatite chondrite plus 32% to 67% ordinary chondrite) though the Mars analogues are generally too massive. However, if we include the restriction of mass on the Mars analogues, the Classical model does not work better. We also further calculate the isotopic composition of $^{17} \rm O$, $^{50} \rm Ti$, $^{54} \rm Cr$, $^{142} \rm Nd$, $^{64} \rm Ni$, and $^{92} \rm Mo$ in the martian mantle from the Grand Tack simulations. We find that it is possible to match the calculated isotopic composition of all the above elements in Mars' mantle with their measured values, but the resulting uncertainties are too large to place good restriction on the early dynamical evolution and birth place of Mars., Comment: 14 pages, 8 figures, presented in the 2017 DPS meeting, 2018 Solar system symposium in Sapporo and 2018 AOGS annual meeting, Accepted for publishing in A&A

We report the detection of planets around two evolved giant stars from radial velocity measurements at Okayama Astrophysical observatory. 24 Boo (G3IV) has a mass of $0.99\,M_{\odot}$, a radius of $10.64\,R_{\odot}$, and a metallicity of ${\rm [Fe/H]}=-0.77$. The star hosts one planet with a minimum mass of $0.91\,M_{\rm Jup}$ and an orbital period of $30.35{\rm d}$. The planet has one of the shortest orbital periods among those ever found around evolved stars by radial-veloocity methods. The stellar radial velocities show additional periodicity with $150{\rm d}$, which are probably attributed to stellar activity. The star is one of the lowest-metallicity stars orbited by planets currently known. $\gamma$ Lib (K0III) is also a metal-poor giant with a mass of $1.47\,M_{\odot}$, a radius of $11.1\,R_{\odot}$, and ${\rm [Fe/H]}=-0.30$. The star hosts two planets with minimum masses of $1.02M_{\rm Jup}$ and $4.58\,M_{\rm Jup}$, and periods of $415{\rm d}$ and $964{\rm d}$, respectively. The star has the second lowest metallicity among the giant stars hosting more than two planets. Dynamical stability analysis for the $\gamma$ Lib system sets a minimum orbital inclination angle to be about $70^{\circ}$ and suggests that the planets are in 7:3 mean-motion resonance, though the current best-fitted orbits to the radial-velocity data are not totally regular., Comment: 29 pages, 19 figures, accepted for publication in PASJ

The icy satellites around Jupiter are considered to have formed in a circumplanetary disk. While previous models focused on the formation of satellites starting from satellitesimals, the question of how satellitesimals form from smaller dust particles has not been addressed so far. In this work, we study the possibility that satellitesimals form in situ in a circumplanetary disk. We calculate the radial distribution of the surface density and representative size of icy dust particles that grow by colliding with each other and drift toward the central planet in a steady circumplanetary disk with a continuous supply of gas and dust from the parent protoplanetary disk. The radial drift barrier is overcome if the ratio of the dust to gas accretion rates onto the circumplanetary disk, $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$, is high and the strength of turbulence, $\alpha$, is not too low. The collision velocity is lower than the critical velocity of fragmentation when $\alpha$ is low. Taken together, we find that the conditions for satellitesimal formation via dust coagulation are given by $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}\ge1$ and $10^{-4}\le\alpha<10^{-2}$. The former condition is generally difficult to achieve, suggesting that the in-situ satellitesimal formation via particle sticking is viable only under an extreme condition. We also show that neither satellitesimal formation via the collisional growth of porous aggregates nor via streaming instability is viable as long as $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$ is low., Comment: 12 pages, 8 figures, 1 table, accepted for publication in ApJ

We have investigated the formation of a circumstellar wide-orbit gas giant planet in a multiple stellar system. We consider a model of orbital circularization for the core of a giant planet after it is scattered from an inner disk region by a more massive planet, which was proposed by Kikuchi et al (2014). We extend their model for single star systems to binary (multiple) star systems, by taking into account tidal truncation of the protoplanetary gas disk by a binary companion. As an example, we consider a wide-orbit gas giant in a hierarchical triple system, HD131399Ab. The best-fit orbit of the planet is that with semimajor axis $\sim 80$ au and eccentricity $\sim 0.35$. As the binary separation is $\sim 350$ au, it is very close to the stability limit, which is puzzling. With the original core location $\sim 20$-30 au, the core (planet) mass $\sim 50 M_{\rm E}$ and the disk truncation radius $\sim 150$ au, our model reproduces the best-fit orbit of HD131399Ab. We find that the orbit after the circularization is usually close to the stability limit against the perturbations from the binary companion, because the scattered core accretes gas from the truncated disk. Our conclusion can also be applied to wider or more compact binary systems if the separation is not too large and another planet with $> \sim$ 20-30 Earth masses that scattered the core existed in inner region of the system., Comment: published in AJ

Context. Planet formation with pebbles has been proposed to solve a couple of long-standing issues in the classical formation model. Some sophisticated simulations have been done to confirm the efficiency of pebble accretion. However, there has not been any global N-body simulations that compare the outcomes of planet formation via pebble accretion with observed extrasolar planetary systems. Aims. In this paper, we study the effects of a range of initial parameters of planet formation via pebble accretion, and present the first results of our simulations. Methods. We incorporate the pebble accretion model by Ida et al. (2016) in the N-body code SyMBA (Duncan et al. 1998), along with the effects of gas accretion, eccentricity and inclination damping and planet migration in the disc. Results. We confirm that pebble accretion leads to a variety of planetary systems, but have difficulty in reproducing observed properties of exoplanetary systems, such as planetary mass, semimajor axis, and eccentricity distributions. The main reason behind this is a too-efficient type I migration, which sensitively depends on the disc model. However, our simulations also lead to a few interesting predictions. First, we find that formation efficiencies of planets depend on the stellar metallicities, not only for giant planets, but also for Earths (Es) and Super-Earths (SEs). The dependency for Es/SEs is subtle. Although higher metallicity environments lead to faster formation of a larger number of Es/SEs, they also tend to be lost later via dynamical instability. Second, our results indicate that a wide range of bulk densities observed for Es and SEs is a natural consequence of dynamical evolution of planetary systems. Third, the ejection trend of our simulations suggest that one free-floating E/SE may be expected for two smaller-mass planets., Comment: Accepted for publication in A&A, 21 pages, 15 figures

Recent theoretical work has shown that the pre-main-sequence (PMS) evolution of stars is much more complex than previously envisioned. Instead of the traditional steady, one-dimensional solution, accretion may be episodic and not necessarily symmetrical, thereby affecting the energy deposited inside the star and its interior structure. Given this new framework, we want to understand what controls the evolution of accreting stars. We use the MESA stellar evolution code with various sets of conditions. In particular, we account for the (unknown) efficiency of accretion in burying gravitational energy into the protostar through a parameter, $\xi$, and we vary the amount of deuterium present. We confirm the findings of previous works that the evolution changes significantly with the amount of energy that is lost during accretion. We find that deuterium burning also regulates the PMS evolution. In the low-entropy accretion scenario, the evolutionary tracks in the H-R diagram are significantly different from the classical tracks and are sensitive to the deuterium content. A comparison of theoretical evolutionary tracks and observations allows us to exclude some cold accretion models ($\xi\sim 0$) with low deuterium abundances. We confirm that the luminosity spread seen in clusters can be explained by models with a somewhat inefficient injection of accretion heat. The resulting evolutionary tracks then become sensitive to the accretion heat efficiency, initial core entropy, and deuterium content. In this context, we predict that clusters with a higher D/H ratio should have less scatter in luminosity than clusters with a smaller D/H. Future work on this issue should include radiation-hydrodynamic simulations to determine the efficiency of accretion heating and further observations to investigate the deuterium content in star-forming regions. (abbrev.), Comment: Published in A&A. 16 pages, 14 figures

We have investigated the probability of temporary capture of asteroids in eccentric orbits by a planet in a circular or an eccentric orbit through analytical and numerical calculations. We found that in the limit of the circular orbit, the capture probability is $\sim 0.1\%$ of encounters to the planet's Hill sphere, independent of planetary mass and semimajor axis. In general, the temporary capture becomes more difficult as the planet's eccentricity ($e_{\rm p}$) increases. We found that the capture probability is almost independent of $e_{\rm p}$ until a critical value ($e_{\rm p}^{\rm c}$) that is given by $\simeq$ 5 times the Hill radius scaled by the planet's semimajor axis. For $e_{\rm p} > e_{\rm p}^{\rm c}$, the probability decreases approximately in proportion to $e_{\rm p}^{-1}$. The current orbital eccentricity of Mars is several times larger than $e_{\rm p}^{\rm c}$. However, since the range of secular change in Martian eccentricity overlaps $e_{\rm p}^{\rm c}$, the capture of minor bodies by the past Mars is not ruled out., Comment: Accepted for publication in AJ, 14 pages and 6 figures

Content: For up to a few millions of years, pebbles must provide a quasi-steady inflow of solids from the outer parts of protoplanetary disks to their inner regions. Aims: We wish to understand how a significant fraction of the pebbles grows into planetesimals instead of being lost to the host star. Methods:We examined analytically how the inward flow of pebbles is affected by the snow line and under which conditions dust-rich (rocky) planetesimals form. When calculating the inward drift of solids that is due to gas drag, we included the back-reaction of the gas to the motion of the solids. Results: We show that in low-viscosity protoplanetary disks (with a monotonous surface density similar to that of the minimum-mass solar nebula), the flow of pebbles does not usually reach the required surface density to form planetesimals by streaming instability. We show, however, that if the pebble-to-gas-mass flux exceeds a critical value, no steady solution can be found for the solid-to-gas ratio. This is particularly important for low-viscosity disks (alpha < 10^(-3)) where we show that inside of the snow line, silicate-dust grains ejected from sublimating pebbles can accumulate, eventually leading to the formation of dust-rich planetesimals directly by gravitational instability. Conclusions: This formation of dust-rich planetesimals may occur for extended periods of time, while the snow line sweeps from several au to inside of 1 au. The rock-to-ice ratio may thus be globally significantly higher in planetesimals and planets than in the central star., Comment: 5 pages, 3 figures; accepted for publication in Astronomy and Astrophysics