-
Dynamics of Binary Planets within Star Clusters
Authors:
Yukun Huang,
Wei Zhu,
Eiichiro Kokubo
Abstract:
We develop analytical tools and perform three-body simulations to investigate the orbital evolution and dynamical stability of binary planets within star clusters. Our analytical results show that the orbital stability of a planetary-mass binary against passing stars is mainly related to its orbital period. We find that critical flybys, defined as stellar encounters with energy kicks comparable to…
▽ More
We develop analytical tools and perform three-body simulations to investigate the orbital evolution and dynamical stability of binary planets within star clusters. Our analytical results show that the orbital stability of a planetary-mass binary against passing stars is mainly related to its orbital period. We find that critical flybys, defined as stellar encounters with energy kicks comparable to the binary binding energy, can efficiently produce a wide range of semimajor axes ($a$) and eccentricities ($e$) from a dominant population of primordially tight JuMBOs. Applying our results to the recently discovered Jupiter-Mass Binary Objects (JuMBOs) by the James Webb Space Telescope (JWST), our simulations suggest that to match the observed $\sim$9\% wide binary fraction, an initial semimajor axis of $a_0 = 10$--20~au and a density-weighted residence time of $χ\gtrsim 10^4$~Myr~pc$^{-3}$ are favored. These results imply that the JWST JuMBOs probably formed as tight binaries near the cluster core.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
-
Gliese 12 b: A temperate Earth-sized planet at 12 pc ideal for atmospheric transmission spectroscopy
Authors:
M. Kuzuhara,
A. Fukui,
J. H. Livingston,
J. A. Caballero,
J. P. de Leon,
T. Hirano,
Y. Kasagi,
F. Murgas,
N. Narita,
M. Omiya,
Jaume Orell-Miquel,
E. Palle,
Q. Changeat,
E. Esparza-Borges,
H. Harakawa,
C. Hellier,
Yasunori Hori,
Kai Ikuta,
H. T. Ishikawa,
T. Kodama,
T. Kotani,
T. Kudo,
J. C. Morales,
M. Mori,
E. Nagel
, et al. (81 additional authors not shown)
Abstract:
Recent discoveries of Earth-sized planets transiting nearby M dwarfs have made it possible to characterize the atmospheres of terrestrial planets via follow-up spectroscopic observations. However, the number of such planets receiving low insolation is still small, limiting our ability to understand the diversity of the atmospheric composition and climates of temperate terrestrial planets. We repor…
▽ More
Recent discoveries of Earth-sized planets transiting nearby M dwarfs have made it possible to characterize the atmospheres of terrestrial planets via follow-up spectroscopic observations. However, the number of such planets receiving low insolation is still small, limiting our ability to understand the diversity of the atmospheric composition and climates of temperate terrestrial planets. We report the discovery of an Earth-sized planet transiting the nearby (12 pc) inactive M3.0 dwarf Gliese 12 (TOI-6251) with an orbital period ($P_{\rm{orb}}$) of 12.76 days. The planet, Gliese 12b, was initially identified as a candidate with an ambiguous $P_{\rm{orb}}$ from TESS data. We confirmed the transit signal and $P_{\rm{orb}}$ using ground-based photometry with MuSCAT2 and MuSCAT3, and validated the planetary nature of the signal using high-resolution images from Gemini/NIRI and Keck/NIRC2 as well as radial velocity (RV) measurements from the InfraRed Doppler instrument on the Subaru 8.2 m telescope and from CARMENES on the CAHA 3.5 m telescope. X-ray observations with XMM-Newton showed the host star is inactive, with an X-ray-to-bolometric luminosity ratio of $\log L_{\rm X}/L_{\rm bol} \approx -5.7$. Joint analysis of the light curves and RV measurements revealed that Gliese 12b has a radius of 0.96 $\pm$ 0.05 $R_\oplus$, a 3$σ$ mass upper limit of 3.9 $M_\oplus$, and an equilibrium temperature of 315 $\pm$ 6 K assuming zero albedo. The transmission spectroscopy metric (TSM) value of Gliese 12b is close to the TSM values of the TRAPPIST-1 planets, adding Gliese 12b to the small list of potentially terrestrial, temperate planets amenable to atmospheric characterization with JWST.
△ Less
Submitted 23 May, 2024;
originally announced May 2024.
-
Global N-body Simulation of Gap Edge Structures Created by Perturbations from a Small Satellite Embedded in Saturn's Rings
Authors:
Naoya Torii,
Shigeru Ida,
Eiichiro Kokubo,
Shugo Michikoshi
Abstract:
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 i…
▽ More
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.
△ Less
Submitted 5 March, 2024;
originally announced March 2024.
-
Two long-period giant planets around two giant stars: HD 112570 and HD 154391
Authors:
Guang-Yao Xiao,
Huan-Yu Teng,
Jianzhao Zhou,
Bun'ei Sato,
Yu-Juan Liu,
Shaolan Bi,
Takuya Takarada,
Masayuki Kuzuhara,
Marc Hon,
Liang Wang,
Masashi Omiya,
Hiroki Harakawa,
Fei Zhao,
Gang Zhao,
Eiji Kambe,
Hideyuki Izumiura,
Hiroyasu Ando,
Kunio Noguchi,
Wei Wang,
Meng Zhai,
Nan Song,
Chengqun Yang,
Tanda Li,
Timothy D. Brandt,
Michitoshi Yoshida
, et al. (2 additional authors not shown)
Abstract:
We present the discoveries of two giant planets orbiting the red giant branch (RGB) star HD 112570 and the red clump (RC) star HD 154391, based on the radial velocity (RV) measurements from Xinglong station and Okayama Astrophysical Observatory (OAO). Spectroscopic and asteroseismic analyses suggest that HD 112570 has a mass of $1.15\pm0.12\,M_{\odot}$, a radius of $9.85\pm0.23\,R_{\odot}$, a meta…
▽ More
We present the discoveries of two giant planets orbiting the red giant branch (RGB) star HD 112570 and the red clump (RC) star HD 154391, based on the radial velocity (RV) measurements from Xinglong station and Okayama Astrophysical Observatory (OAO). Spectroscopic and asteroseismic analyses suggest that HD 112570 has a mass of $1.15\pm0.12\,M_{\odot}$, a radius of $9.85\pm0.23\,R_{\odot}$, a metallicity [Fe/H] of $-0.46\pm0.1$ and a ${\rm log}\,g$ of $2.47\pm0.1$. With the joint analysis of RV and Hipparcos-Gaia astrometry, we obtain a dynamical mass of $M_{\rm p}={3.42}_{-0.84}^{+1.4}\ M_{\rm Jup}$, a period of $P={2615}_{-77}^{+85}$ days and a moderate eccentricity of $e={0.20}_{-0.14}^{+0.16}$ for the Jovian planet HD 112570 b. For HD 154391, it has a mass of $2.07\pm0.03\,M_{\odot}$, a radius of $8.56\pm0.05\,R_{\odot}$, a metallicity [Fe/H] of $0.07\pm0.1$ and a ${\rm log}\,g$ of $2.86\pm0.1$. The super-Jupiter HD 154391 b has a mass of $M_{\rm p}={9.1}_{-1.9}^{+2.8}\ M_{\rm Jup}$, a period of $P={5163}_{-57}^{+60}$ days and an eccentricity of $e={0.20}_{-0.04}^{+0.04}$. We found HD 154391 b has one of the longest orbital period among those ever discovered orbiting evolved stars, which may provide a valuable case in our understanding of planetary formation at wider orbits. Moreover, while a mass gap at $4\,M_{\rm Jup}$ seems to be present in the population of giant stars, there appears to be no significant differences in the distribution of metallicity among giant planets with masses above or below this threshold. Finally, The origin of the abnormal accumulation near 2 au for planets around large evolved stars ($R_{\star}>21\,R_{\odot}$), remains unclear.
△ Less
Submitted 3 December, 2023;
originally announced December 2023.
-
Formation of inner planets in the presence of a Cold Jupiter: orbital evolution and relative velocities of planetesimals
Authors:
Kangrou Guo,
Eiichiro Kokubo
Abstract:
We investigate the orbital evolution of planetesimals in the inner disk in the presence of nebula gas and a (proto-) cold Jupiter. By varying the mass, eccentricity, and semi-major axis of the planet, we study the dependence of the relative velocities of the planetesimals on these parameters. For classic small planetesimals ($10^{16}-10^{20} $g) whose mutual gravitational interaction is negligible…
▽ More
We investigate the orbital evolution of planetesimals in the inner disk in the presence of nebula gas and a (proto-) cold Jupiter. By varying the mass, eccentricity, and semi-major axis of the planet, we study the dependence of the relative velocities of the planetesimals on these parameters. For classic small planetesimals ($10^{16}-10^{20} $g) whose mutual gravitational interaction is negligible, gas drag introduces a size-dependent alignment of orbits and keeps the relative velocity low for similar-size bodies, while preventing orbital alignment for different-size planetesimals. Regardless of the location and the mass ratio of the planetesimals, increasing the mass and eccentricity or decreasing the orbital distance of the planet always leads to higher relative velocities of planetesimals. However, for massive planetesimals, the interplay of viscous stirring, gas damping, and secular perturbation results in lower velocity dispersion of equal-size planetesimals when the planet is more massive or when it is located on a closer or more eccentric orbit. The random velocities of such planetesimals remain almost unperturbed when the planet is located beyond Jupiter's current orbit, or when it is less massive or less eccentric than Jupiter. Unlike small planetesimals, such large planetesimals can grow in a runaway fashion as in the unperturbed case. Our results imply that the presence of a cold Jupiter does not impede the formation of inner rocky planets through planetesimal accretion, provided that the planetesimals are initially large.
△ Less
Submitted 28 August, 2023;
originally announced August 2023.
-
Revisiting Planetary Systems in Okayama Planet Search Program: A new long-period planet, RV astrometry joint analysis, and multiplicity-metallicity trend around evolved stars
Authors:
Huan-Yu Teng,
Bun'ei Sato,
Masayuki Kuzuhara,
Takuya Takarada,
Masashi Omiya,
Hiroki Harakawa,
Hideyuki Izumiura,
Eiji Kambe,
Mesut Yilmaz,
Ilfan Bikmaev,
Selim O. Selam,
Timothy D. Brandt,
Guang-Yao Xiao,
Michitoshi Yoshida,
Yoichi Itoh,
Hiroyasu Ando,
Eiichiro Kokubo,
Shigeru Ida
Abstract:
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 con…
▽ More
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 $ε$ Tau, 11 Com, 24 Boo, 41 Lyn, 14 And, HD 32518, and $ω$ 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.
△ Less
Submitted 10 August, 2023;
originally announced August 2023.
-
Size Dependence of the Bouncing Barrier in Protoplanetary Dust Growth
Authors:
Sota Arakawa,
Satoshi Okuzumi,
Misako Tatsuuma,
Hidekazu Tanaka,
Eiichiro Kokubo,
Daisuke Nishiura,
Mikito Furuichi,
Taishi Nakamoto
Abstract:
Understanding the collisional behavior of dust aggregates is essential in the context of planet formation. It is known that low-velocity collisions of dust aggregates result in bouncing rather than sticking when the filling factor of colliding dust aggregates is higher than a threshold value. However, a large discrepancy between numerical and experimental results on the threshold filling factor wa…
▽ More
Understanding the collisional behavior of dust aggregates is essential in the context of planet formation. It is known that low-velocity collisions of dust aggregates result in bouncing rather than sticking when the filling factor of colliding dust aggregates is higher than a threshold value. However, a large discrepancy between numerical and experimental results on the threshold filling factor was reported so far. In this study, we perform numerical simulations using soft-sphere discrete element methods and demonstrate that the sticking probability decreases with increasing aggregates radius. Our results suggest that the large discrepancy in the threshold filling factor may reflect the difference in the size of dust aggregates in earlier numerical simulations and laboratory experiments.
△ Less
Submitted 6 June, 2023;
originally announced June 2023.
-
Insights on the Sun birth environment in the context of star-cluster formation in hub-filament systems
Authors:
Doris Arzoumanian,
Sota Arakawa,
Masato I. N. Kobayashi,
Kazunari Iwasaki,
Kohei Fukada,
Shoji Mori,
Yutaka Hirai,
Masanobu Kunitomo,
M. S. Nanda Kumar,
Eiichiro Kokubo
Abstract:
Cylindrical molecular filaments are observed to be the main sites of Sun-like star formation, while massive stars form in dense hubs, at the junction of multiple filaments. The role of hub-filament configurations has not been discussed yet in relation to the birth environment of the solar system and to infer the origin of isotopic ratios of Short-Lived Radionuclides (SLR, such as $^{26}$Al) of Cal…
▽ More
Cylindrical molecular filaments are observed to be the main sites of Sun-like star formation, while massive stars form in dense hubs, at the junction of multiple filaments. The role of hub-filament configurations has not been discussed yet in relation to the birth environment of the solar system and to infer the origin of isotopic ratios of Short-Lived Radionuclides (SLR, such as $^{26}$Al) of Calcium-Aluminum-rich Inclusions (CAIs) observed in meteorites. In this work, we present simple analytical estimates of the impact of stellar feedback on the young solar system forming along a filament of a hub-filament system. We find that the host filament can shield the young solar system from the stellar feedback, both during the formation and evolution of stars (stellar outflow, wind, and radiation) and at the end of their life (supernovae). We show that the young solar system formed along a dense filament can be enriched with supernova ejecta (e.g., $^{26}$Al) during the formation timescale of CAIs. We also propose that the streamers recently observed around protostars may be channeling the SLR-rich material onto the young solar system. We conclude that considering hub-filament configurations as the birth environment of the Sun is important when deriving theoretical models explaining the observed properties of the solar system.
△ Less
Submitted 27 March, 2023;
originally announced March 2023.
-
Direct Imaging Explorations for Companions around Mid-Late M Stars from the Subaru/IRD Strategic Program
Authors:
Taichi Uyama,
Charles Beichman,
Masayuki Kuzuhara,
Markus Janson,
Takayuki Kotani,
Dimitri Mawet,
Bun'ei Sato,
Motohide Tamura,
Hiroyuki Tako Ishikawa,
Bryson Cale,
Thayne Currie,
Hiroki Harakawa,
Thomas Henning,
Teruyuki Hirano,
Klaus Hodapp,
Yasunori Hori,
Masato Ishizuka,
Shane Jacobson,
Yui Kasagi,
Eiichiro Kokubo,
Mihoko Konishi,
Tomoyuki Kudo,
Takashi Kurokawa,
Nobuhiko Kusakabe,
Jungmi Kwon
, et al. (10 additional authors not shown)
Abstract:
The Subaru telescope is currently performing a strategic program (SSP) using the high-precision near-infrared (NIR) spectrometer IRD to search for exoplanets around nearby mid/late-M~dwarfs via radial velocity (RV) monitoring. As part of the observing strategy for the exoplanet survey, signatures of massive companions such as RV trends are used to reduce the priority of those stars. However, this…
▽ More
The Subaru telescope is currently performing a strategic program (SSP) using the high-precision near-infrared (NIR) spectrometer IRD to search for exoplanets around nearby mid/late-M~dwarfs via radial velocity (RV) monitoring. As part of the observing strategy for the exoplanet survey, signatures of massive companions such as RV trends are used to reduce the priority of those stars. However, this RV information remains useful for studying the stellar multiplicity of nearby M~dwarfs. To search for companions around such ``deprioritized" M~dwarfs, we observed 14 IRD-SSP targets using Keck/NIRC2 observations with pyramid wavefront sensing at NIR wavelengths, leading to high sensitivity to substellar-mass companions within a few arcseconds. We detected two new companions (LSPM~J1002+1459~B and LSPM~J2204+1505~B) and two new candidates that are likely companions (LSPM~J0825+6902~B and LSPM~J1645+0444~B) as well as one known companion. Including two known companions resolved by the IRD fiber injection module camera, we detected seven (four new) companions at projected separations between $\sim2-20$~au in total. A comparison of the colors with the spectral library suggests that LSPM~J2204+1505~B and LSPM~J0825+6902~B are located at the boundary between late-M and early-L spectral types. Our deep high-contrast imaging for targets where no bright companions were resolved did not reveal any additional companion candidates. The NIRC2 detection limits could constrain potential substellar-mass companions ($\sim10-75\ M_{\rm Jup}$) at 10~au or further. The failure with Keck/NIRC2 around the IRD-SSP stars having significant RV trends makes these objects promising targets for further RV monitoring or deeper imaging with JWST to search for smaller-mass companions below the NIRC2 detection limits.
△ Less
Submitted 26 February, 2023;
originally announced February 2023.
-
An Earth-sized Planet around an M5 Dwarf Star at 22 pc
Authors:
Teruyuki Hirano,
Fei Dai,
John H. Livingston,
Yui Kasagi,
Norio Narita,
Hiroyuki Tako Ishikawa,
Sascha Grziwa,
Kristine W. F. Lam,
Kohei Miyakawa,
Luisa M. Serrano,
Yuji Matsumoto,
Eiichiro Kokubo,
Tadahiro Kimura,
Masahiro Ikoma,
Joshua N. Winn,
John P. Wisniewski,
Hiroki Harakawa,
Huan-Yu Teng,
William D. Cochran,
Akihiko Fukui,
Davide Gandolfi,
Eike W. Guenther,
Yasunori Hori,
Kai Ikuta,
Kiyoe Kawauchi
, et al. (24 additional authors not shown)
Abstract:
We report on the discovery of an Earth-sized transiting planet ($R_p=1.015\pm0.051\,R_\oplus$) in a $P=4.02$ day orbit around K2-415 (EPIC 211414619), an M5V star at 22 pc. The planet candidate was first identified by analyzing the light curve data by the K2 mission, and is here shown to exist in the most recent data from TESS. Combining the light curves with the data secured by our follow-up obse…
▽ More
We report on the discovery of an Earth-sized transiting planet ($R_p=1.015\pm0.051\,R_\oplus$) in a $P=4.02$ day orbit around K2-415 (EPIC 211414619), an M5V star at 22 pc. The planet candidate was first identified by analyzing the light curve data by the K2 mission, and is here shown to exist in the most recent data from TESS. Combining the light curves with the data secured by our follow-up observations including high-resolution imaging and near infrared spectroscopy with IRD, we rule out false positive scenarios, finding a low false positive probability of $2\times 10^{-4}$. Based on IRD's radial velocities of K2-415, which were sparsely taken over three years, we obtain the planet mass of $3.0\pm 2.7\,M_\oplus$ ($M_p<7.5\,M_\oplus$ at $95\,\%$ confidence) for K2-415b. Being one of the lowest mass stars ($\approx 0.16\,M_\odot$) known to host an Earth-sized transiting planet, K2-415 will be an interesting target for further follow-up observations, including additional radial velocity monitoring and transit spectroscopy.
△ Less
Submitted 1 February, 2023;
originally announced February 2023.
-
Threshold velocity for collisional growth of porous dust aggregates consisting of cohesive frictionless spheres
Authors:
Sota Arakawa,
Hidekazu Tanaka,
Eiichiro Kokubo,
Daisuke Nishiura,
Mikito Furuichi
Abstract:
Understanding the collisional outcomes of dust aggregates and dependence on material properties of the constituting particles is of great importance toward understanding planet formation. Recent numerical simulations have revealed that interparticle tangential friction plays a crucial role in energy dissipation during collisions between porous dust aggregates; however, the importance of friction o…
▽ More
Understanding the collisional outcomes of dust aggregates and dependence on material properties of the constituting particles is of great importance toward understanding planet formation. Recent numerical simulations have revealed that interparticle tangential friction plays a crucial role in energy dissipation during collisions between porous dust aggregates; however, the importance of friction on the collisional growth of dust aggregates remains poorly understood. Here we demonstrate the effects of interparticle tangential friction on the collisional growth of dust aggregates. We performed numerical simulations of collisions between equal-mass porous dust aggregates consisting of cohesive and frictionless spheres. We changed the collision velocity and impact angle systematically and calculated the collisional growth efficiency as a function of the collision velocity. We found that the threshold velocity for collisional growth decreases when dust aggregates are made of frictionless spheres as compared to frictional spheres. Our results highlight the importance of tangential interactions on the collisional behavior of dust aggregates and indicate that the predictive equation for threshold velocity should be reconstructed.
△ Less
Submitted 1 February, 2023;
originally announced February 2023.
-
On the Number of Stars in the Sun's Birth Cluster
Authors:
Sota Arakawa,
Eiichiro Kokubo
Abstract:
The Sun is thought to be formed within a star cluster. The coexistence of $^{26}{\rm Al}$-rich and $^{26}{\rm Al}$-poor calcium--aluminum-rich inclusions indicates that a direct injection of $^{26}{\rm Al}$-rich materials from a nearby core-collapse supernova should occur in the first $10^5$ years of the solar system. Therefore, at least one core-collapse supernova should occur within the duration…
▽ More
The Sun is thought to be formed within a star cluster. The coexistence of $^{26}{\rm Al}$-rich and $^{26}{\rm Al}$-poor calcium--aluminum-rich inclusions indicates that a direct injection of $^{26}{\rm Al}$-rich materials from a nearby core-collapse supernova should occur in the first $10^5$ years of the solar system. Therefore, at least one core-collapse supernova should occur within the duration of star formation in the Sun's birth cluster. Here we revisit the number of stars in the Sun's birth cluster from the point of view of the probability for acquiring at least one core-collapse supernova within the finite duration of star formation in the birth cluster. We find that the number of stars in the birth cluster can be significantly larger than that previously considered, depending on the duration of star formation.
△ Less
Submitted 28 December, 2022;
originally announced December 2022.
-
Orbital structure of planetary systems formed by giant impacts: stellar mass dependence
Authors:
Haruka Hoshino,
Eiichiro Kokubo
Abstract:
Recent exoplanet surveys revealed that for solar-type stars, close-in Super-Earths are ubiquitous and many of them are in multi-planet systems. These systems are more compact than the Solar System's terrestrial planets. However, there have been few theoretical studies on the formation of such planets around low-mass stars. In the standard model, the final stage of terrestrial planet formation is t…
▽ More
Recent exoplanet surveys revealed that for solar-type stars, close-in Super-Earths are ubiquitous and many of them are in multi-planet systems. These systems are more compact than the Solar System's terrestrial planets. However, there have been few theoretical studies on the formation of such planets around low-mass stars. In the standard model, the final stage of terrestrial planet formation is the giant impact stage, where protoplanets gravitationally scatter and collide with each other and then evolve into a stable planetary system. We investigate the effect of the stellar mass on the architecture of planetary systems formed by giant impacts. We perform {\it N}-body simulations around stars with masses of 0.1--2 times the solar mass. Using the isolation mass of protoplanets, we distribute the initial protoplanets in 0.05--0.15 au from the central star and follow the evolution for 200 million orbital periods of the innermost protoplanet. We find that for a given protoplanet system, the mass of planets increases as the stellar mass decreases, while the number of planets decreases. The eccentricity and inclination of orbits and the orbital separation of adjacent planets increase with decreasing the stellar mass. This is because as the stellar mass decreases, the relative strength of planetary scattering becomes more effective. We also discuss the properties of planets formed in the habitable zone using the minimum-mass extrasolar nebula model.
△ Less
Submitted 17 December, 2022;
originally announced December 2022.
-
A Close-in Planet Orbiting Giant Star HD 167768
Authors:
Huan-Yu Teng,
Bun'ei Sato,
Masanobu Kunitomo,
Takuya Takarada,
Masashi Omiya,
Hiroki Harakawa,
Guang-Yao Xiao,
Yu-Juan Liu,
Hideyuki Izumiura,
Eiji Kambe,
Michitoshi Yoshida,
Yoichi Itoh,
Hiroyasu Ando,
Eiichiro Kokubo,
Ida Shigeru
Abstract:
We report the detection of a giant planet orbiting a G-type giant star HD 167768 from radial velocity measurements using HIgh Dispersion Echelle Spectrograph (HIDES) at Okayama Astrophysical Observatory (OAO). HD 167768 has a mass of $1.08_{-0.12}^{+0.14} M_{\odot}$, a radius of $9.70_{-0.25}^{+0.25} R_{\odot}$, a metallicity of $\rm{[Fe/H]}=-0.67_{-0.08}^{+0.09}$, and a surface gravity of…
▽ More
We report the detection of a giant planet orbiting a G-type giant star HD 167768 from radial velocity measurements using HIgh Dispersion Echelle Spectrograph (HIDES) at Okayama Astrophysical Observatory (OAO). HD 167768 has a mass of $1.08_{-0.12}^{+0.14} M_{\odot}$, a radius of $9.70_{-0.25}^{+0.25} R_{\odot}$, a metallicity of $\rm{[Fe/H]}=-0.67_{-0.08}^{+0.09}$, and a surface gravity of $\log g = 2.50_{-0.06}^{+0.06}$. The planet orbiting the star is a warm Jupiter, having a period of $20.6532_{-0.0032}^{+0.0032}\ \rm{d}$, a minimum mass of $0.85_{-0.11}^{+0.12}\ M_{\rm{J}}$, and an orbital semimajor axis of $0.1512_{-0.0063}^{+0.0058}\ \rm{au}$. The planet has one of the shortest orbital periods among those ever found around deeply evolved stars ($\log g < 3.5$) using radial velocity methods. The equilibrium temperature of the planet is $1874\ \rm{K}$, as high as a hot Jupiter. The radial velocities show two additional regular variations at $41\ \rm{d}$ and $95\ \rm{d}$, suggesting the possibility of outer companions in the system. Follow-up monitoring will enable validation of the periodicity. We also calculated the orbital evolution of HD 167768 b and found that the planet will be engulfed within 0.15\,Gyr.
△ Less
Submitted 12 November, 2022;
originally announced November 2022.
-
Planetesimal formation by the gravitational instability of dust ring structures
Authors:
Sanemichi Z. Takahashi,
Eiichiro Kokubo,
Shu-ichiro Inutsuka
Abstract:
We investigate the gravitational instability (GI) of dust-ring structures and the formation of planetesimals by their gravitational collapse. The normalized dispersion relation of a self-gravitating ring structure includes two parameters that are related to its width and line mass (the mass per unit length). We survey these parameters and calculate the growth rate and wavenumber. Additionally, we…
▽ More
We investigate the gravitational instability (GI) of dust-ring structures and the formation of planetesimals by their gravitational collapse. The normalized dispersion relation of a self-gravitating ring structure includes two parameters that are related to its width and line mass (the mass per unit length). We survey these parameters and calculate the growth rate and wavenumber. Additionally, we investigate the planetesimal formation by growth of the GI of the ring that is formed by the growth of the secular GI of the protoplanetary disk. We adopt a massive, dust rich disk as a disk model. We find the range of radii for the fragmentation by the ring GI as a function of the width of the ring. The inner-most radius for the ring GI is smaller for the smaller ring width. We also determine the range of the initial planetesimal mass resulting from the fragmentation of the ring GI. Our results indicate that the planetesimal mass can be as large as 10^28 g at its birth after the fragmentation. It can be as low as about 10^25 g if the ring width is 0.1% of the ring radius and the lower limit increases with the ring width. Furthermore, we obtain approximate formulas for the upper and lower limits of the planetesimal mass. We predict that the planetesimals formed by the ring GI have prograde rotations because of the Coriolis force acting on the contracting dust. This is consistent with the fact that many trans-Neptunian binaries exhibit prograde rotation.
△ Less
Submitted 4 November, 2022;
originally announced November 2022.
-
Collisional growth efficiency of dust aggregates and its independence of the strength of interparticle rolling friction
Authors:
Sota Arakawa,
Hidekazu Tanaka,
Eiichiro Kokubo
Abstract:
The pairwise collisional growth of dust aggregates consisting submicron-sized grains is the first step of the planet formation, and understanding the collisional behavior of dust aggregates is therefore essential. It is known that the main energy dissipation mechanisms are the tangential frictions between particles in contact, namely, rolling, sliding, and twisting. However, there is a large uncer…
▽ More
The pairwise collisional growth of dust aggregates consisting submicron-sized grains is the first step of the planet formation, and understanding the collisional behavior of dust aggregates is therefore essential. It is known that the main energy dissipation mechanisms are the tangential frictions between particles in contact, namely, rolling, sliding, and twisting. However, there is a large uncertainty for the strength of rolling friction, and the dependence of the collisional growth condition on the strength of rolling friction was poorly understood. Here we performed numerical simulations of collisions between two equal-mass porous aggregates with various collision velocities and impact parameters, and we also changed the strength of rolling friction systematically. We found that the threshold of the collision velocity for the fragmentation of dust aggregates is nearly independent of the strength of rolling friction. This is because the total amount of the energy dissipation by the tangential frictions is nearly constant even though the strength of rolling friction is varied.
△ Less
Submitted 3 October, 2022;
originally announced October 2022.
-
A Trio of Giant Planets Orbiting Evolved Star HD 184010
Authors:
Huan-Yu Teng,
Bun'ei Sato,
Takuya Takarada,
Masashi Omiya,
Hiroki Harakawa,
Makiko Nagasawa,
Ryo Hasegawa,
Hideyuki Izumiura,
Eiji Kambe,
Michitoshi Yoshida,
Yoichi Itoh,
Hiroyasu Ando,
Eiichiro Kokubo,
Shigeru Ida
Abstract:
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…
▽ More
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}$).
△ Less
Submitted 19 September, 2022;
originally announced September 2022.
-
Planetesimal Dynamics in the Presence of a Giant Planet II: Dependence on Planet Mass and Eccentricity
Authors:
Kangrou Guo,
Eiichiro Kokubo
Abstract:
The presence of an early-formed giant planet in the protoplanetary disk has mixed influence on the growth of other planetary embryos. Gravitational perturbation from the planet can increase the relative velocities of planetesimals at the mean motion resonances to very high values and impede accretion at those locations. However, gas drag can also align the orbital pericenters of equal-size planete…
▽ More
The presence of an early-formed giant planet in the protoplanetary disk has mixed influence on the growth of other planetary embryos. Gravitational perturbation from the planet can increase the relative velocities of planetesimals at the mean motion resonances to very high values and impede accretion at those locations. However, gas drag can also align the orbital pericenters of equal-size planetesimals in certain disk locations and make them dynamically quiet and "accretion-friendly" locations for planetesimals of similar sizes. Following the previous paper, where we investigated the effect of a Jupiter-like planet on an external planetesimal disk, we generalize our findings to extrasolar planetary systems by varying the planet parameters. In particular, we focus on the dependence of the planetesimal relative velocities on the mass and eccentricity of the existing planet. We found that the velocity dispersion of identical-mass particles increases monotonically with increasing planet mass. Meanwhile, the dependence of the relative velocity between different-mass planetesimals on their mass ratio becomes weaker as the planet mass increases. While the relative velocities generally increases with increasing planet eccentricity, the velocity dispersion of smaller-mass particles ($m \lesssim 10^{18}~\rm{g}$) is almost independent of planet eccentricity owing to their strong coupling to gas. We find that the erosion limits are met for a wider range of parameters (planet mass/eccentricity, planetesimal mass ratio) when the planetesimal size decreases. Our results could provide some clues for the formation of Saturn's core as well as the architecture of some exoplanetary systems with multiple cold giant planets.
△ Less
Submitted 27 July, 2022;
originally announced July 2022.
-
Impacts of viscous dissipation on collisional growth and fragmentation of dust aggregates
Authors:
Sota Arakawa,
Hidekazu Tanaka,
Eiichiro Kokubo
Abstract:
Understanding the collisional behavior of dust aggregates consisting of submicron-sized grains is essential to unveiling how planetesimals formed in protoplanetary disks. It is known that the collisional behavior of individual dust particles strongly depends on the strength of viscous dissipation force; however, impacts of viscous dissipation on the collisional behavior of dust aggregates have not…
▽ More
Understanding the collisional behavior of dust aggregates consisting of submicron-sized grains is essential to unveiling how planetesimals formed in protoplanetary disks. It is known that the collisional behavior of individual dust particles strongly depends on the strength of viscous dissipation force; however, impacts of viscous dissipation on the collisional behavior of dust aggregates have not been studied in detail, especially for the cases of oblique collisions. Here we investigated the impacts of viscous dissipation on the collisional behavior of dust aggregates. We performed numerical simulations of collisions between two equal-mass dust aggregates with various collision velocities and impact parameters. We also changed the strength of viscous dissipation force systematically. We found that the threshold collision velocity for the fragmentation of dust aggregates barely depends on the strength of viscous dissipation force when we consider oblique collisions. In contrast, the size distribution of fragments changes significantly when the viscous dissipation force is considered. We obtained the empirical fitting formulae for the size distribution of fragments for the case of strong dissipation, which would be useful to study the evolution of size and spatial distributions of dust aggregates in protoplanetary disks.
△ Less
Submitted 27 May, 2022;
originally announced May 2022.
-
A Super-Earth Orbiting Near the Inner Edge of the Habitable Zone around the M4.5-dwarf Ross 508
Authors:
Hiroki Harakawa,
Takuya Takarada,
Yui Kasagi,
Teruyuki Hirano,
Takayuki Kotani,
Masayuki Kuzuhara,
Masashi Omiya,
Hajime Kawahara,
Akihiko Fukui,
Yasunori Hori,
Hiroyuki Tako Ishikawa,
Masahiro Ogihara,
John Livingston,
Timothy D. Brandt,
Thayne Currie,
Wako Aoki,
Charles A. Beichman,
Thomas Henning,
Klaus Hodapp,
Masato Ishizuka,
Hideyuki Izumiura,
Shane Jacobson,
Markus Janson,
Eiji Kambe,
Takanori Kodama
, et al. (24 additional authors not shown)
Abstract:
We report the near-infrared radial-velocity (RV) discovery of a super-Earth planet on a 10.77-day orbit around the M4.5 dwarf Ross 508 ($J_\mathrm{mag}=9.1$). Using precision RVs from the Subaru Telescope IRD (InfraRed Doppler) instrument, we derive a semi-amplitude of $3.92^{+0.60}_{-0.58}$ ${\rm m\,s}^{-1}$, corresponding to a planet with a minimum mass…
▽ More
We report the near-infrared radial-velocity (RV) discovery of a super-Earth planet on a 10.77-day orbit around the M4.5 dwarf Ross 508 ($J_\mathrm{mag}=9.1$). Using precision RVs from the Subaru Telescope IRD (InfraRed Doppler) instrument, we derive a semi-amplitude of $3.92^{+0.60}_{-0.58}$ ${\rm m\,s}^{-1}$, corresponding to a planet with a minimum mass $m \sin i = 4.00^{+0.53}_{-0.55}\ M_{\oplus}$. We find no evidence of significant signals at the detected period in spectroscopic stellar activity indicators or MEarth photometry. The planet, Ross 508 b, has a semimajor-axis of $0.05366^{+0.00056}_{-0.00049}$ au. This gives an orbit-averaged insolation of $\approx$1.4 times the Earth's value, placing Ross 508 b near the inner edge of its star's habitable zone. We have explored the possibility that the planet has a high eccentricity and its host is accompanied by an additional unconfirmed companion on a wide orbit. Our discovery demonstrates that the near-infrared RV search can play a crucial role to find a low-mass planet around cool M dwarfs like Ross 508.
△ Less
Submitted 24 May, 2022;
originally announced May 2022.
-
Rapid-then-slow migration reproduces mass distribution of TRAPPIST-1 system
Authors:
Masahiro Ogihara,
Eiichiro Kokubo,
Ryuunosuke Nakano,
Takeru K. Suzuki
Abstract:
The TRAPPIST-1 system is an iconic planetary system in various aspects (e.g., habitability, resonant relation, and multiplicity) and hence has attracted considerable attention. The mass distribution of the TRAPPIST-1 planets is characterized by two features: the two inner planets are large, and the masses of the four planets in the outer orbit increase with orbital distance. The origin of these fe…
▽ More
The TRAPPIST-1 system is an iconic planetary system in various aspects (e.g., habitability, resonant relation, and multiplicity) and hence has attracted considerable attention. The mass distribution of the TRAPPIST-1 planets is characterized by two features: the two inner planets are large, and the masses of the four planets in the outer orbit increase with orbital distance. The origin of these features cannot be explained by previous formation models. We investigate whether the mass distribution of the TRAPPIST-1 system can be reproduced by a planet formation model using N-body simulations. We used a gas disk evolution model around a low-mass star constructed by considering disk winds and followed the growth and orbital migration from planetary embryos with the isolation mass, which increases with orbital distance. As a result, we find that from the initial phase, planets in inner orbits undergo rapid orbital migration, and the coalescence growth near the inner disk edge is enhanced. This allows the inner planets to grow larger. Meanwhile, compared with the inner planets, planets in outer orbits migrate more slowly and do not frequently collide with neighboring planets. Therefore, the trend of increasing mass toward the outer orbit, called reversed mass ranking, is maintained. The final mass distribution approximately agrees with the two features of the mass distribution in the TRAPPIST-1 system. We discover that the mass distribution in the TRAPPIST-1 system can be reproduced when embryos experience rapid migration and become trapped near the disk inner edge, and then more massive embryos undergo slower migration. This migration transition can be achieved naturally in a disk evolution model with disk winds.
△ Less
Submitted 21 January, 2022;
originally announced January 2022.
-
Regular Radial Velocity Variations in Nine G- and K-type Giant Stars: Eight Planets and One Planet Candidate
Authors:
Huan-Yu Teng,
Bun'ei Sato,
Takuya Takarada,
Masashi Omiya,
Hiroki Harakawa,
Hideyuki Izumiura,
Eiji Kambe,
Yoichi Takeda,
Michitoshi Yoshida,
Yoichi Itoh,
Hiroyasu Ando,
Eiichiro Kokubo
Abstract:
We report the detection of radial velocity variations in nine evolved G- and K-type giant stars. The observations were conducted at Okayama Astrophysical Observatory. Planets or planet candidates can best explain these regular variations. However, a coincidence of near 280-day variability among five of them prevents us from fully ruling out stellar origins for some of the variations, since all nin…
▽ More
We report the detection of radial velocity variations in nine evolved G- and K-type giant stars. The observations were conducted at Okayama Astrophysical Observatory. Planets or planet candidates can best explain these regular variations. However, a coincidence of near 280-day variability among five of them prevents us from fully ruling out stellar origins for some of the variations, since all nine stars behave similarly in stellar properties. In the planet hypotheses to the RV variations, the planets (including one candidate) may survive close to the boundary of the so-called "planet desert" around evolved stars, having orbital periods between 255 and 555 days. Besides, they are the least-massive giant planets detected around G- and K-type giant stars, with minimum masses between 0.45$M_{\rm{J}}$ and 1.34$M_{\rm{J}}$. We further investigated other hypotheses for our detection, yet none of them can better explain regular RV variation. With our detection, it is convinced that year-long regular variation with amplitude down to 15 $\rm{m\ s^{-1}}$ for G- and K-type giant stars is detectable. Moreover, we performed simulations to further confirm the detectability of planets around these stars. Finally, we explored giant planets around intermediate-mass stars, and likewise found a 4 Jupiter mass gap (e.g. Santo et al. 2017), which is probably a boundary of the giant planet population.
△ Less
Submitted 14 December, 2021;
originally announced December 2021.
-
Elemental abundances of nearby M dwarfs based on high-resolution near-infrared spectra obtained by the Subaru/IRD survey: Proof of concept
Authors:
Hiroyuki Tako Ishikawa,
Wako Aoki,
Teruyuki Hirano,
Takayuki Kotani,
Masayuki Kuzuhara,
Masashi Omiya,
Yasunori Hori,
Eiichiro Kokubo,
Tomoyuki Kudo,
Takashi Kurokawa,
Nobuhiko Kusakabe,
Norio Narita,
Jun Nishikawa,
Masahiro Ogihara,
Akitoshi Ueda,
Thayne Currie,
Thomas Henning,
Yui Kasagi,
Jared R. Kolecki,
Jungmi Kwon,
Masahiro N. Machida,
Michael W. McElwain,
Takao Nakagawa,
Sebastien Vievard,
Ji Wang
, et al. (2 additional authors not shown)
Abstract:
Detailed chemical analyses of M dwarfs are scarce but necessary to constrain the formation environment and internal structure of planets being found around them. We present elemental abundances of 13 M dwarfs (2900 < Teff < 3500 K) observed in the Subaru/IRD planet search project. They are mid-to-late M dwarfs whose abundance of individual elements has not been well studied. We use the high-resolu…
▽ More
Detailed chemical analyses of M dwarfs are scarce but necessary to constrain the formation environment and internal structure of planets being found around them. We present elemental abundances of 13 M dwarfs (2900 < Teff < 3500 K) observed in the Subaru/IRD planet search project. They are mid-to-late M dwarfs whose abundance of individual elements has not been well studied. We use the high-resolution (~70,000) near-infrared (970-1750 nm) spectra to measure the abundances of Na, Mg, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Sr by the line-by-line analysis based on model atmospheres, with typical errors ranging from 0.2 dex for [Fe/H] to 0.3-0.4 dex for other [X/H]. We measure radial velocities from the spectra and combine them with Gaia astrometry to calculate the Galactocentric space velocities UVW. The resulting [Fe/H] values agree with previous estimates based on medium-resolution K-band spectroscopy, showing a wide distribution of metallicity (-0.6 < [Fe/H] < +0.4). The abundance ratios of individual elements [X/Fe] are generally aligned with the solar values in all targets. While the [X/Fe] distributions are comparable to those of nearby FGK stars, most of which belong to the thin disk population, the most metal-poor object, GJ 699, could be a thick disk star. The UVW velocities also support this. The results raise the prospect that near-infrared spectra of M dwarfs obtained in the planet search projects can be used to grasp the trend of elemental abundances and Galactic stellar population of nearby M dwarfs.
△ Less
Submitted 30 November, 2021;
originally announced December 2021.
-
Early Initiation of Inner Solar System Formation at Dead-Zone Inner Edge
Authors:
Takahiro Ueda,
Masahiro Ogihara,
Eiichiro Kokubo,
Satoshi Okuzumi
Abstract:
The inner solar system possesses a unique orbital structure in which there are no planets inside the Mercury orbit and the mass is concentrated around the Venus and Earth orbits. The origins of these features still remain unclear. We propose a novel concept that the building blocks of the inner solar system formed at the dead-zone inner edge in the early phase of the protosolar disk evolution, whe…
▽ More
The inner solar system possesses a unique orbital structure in which there are no planets inside the Mercury orbit and the mass is concentrated around the Venus and Earth orbits. The origins of these features still remain unclear. We propose a novel concept that the building blocks of the inner solar system formed at the dead-zone inner edge in the early phase of the protosolar disk evolution, where the disk is effectively heated by the disk accretion. First, we compute the dust evolution in a gas disk with a dead zone and obtain the spatial distribution of rocky planetesimals. The disk is allowed to evolve both by a viscous diffusion and magnetically-driven winds. We find that the rocky planetesimals are formed in concentrations around $\sim$ 1 au with a total mass comparable to the mass of the current inner solar system in the early phase of the disk evolution within $\lesssim0.1$ Myr. Based on the planetesimal distribution and the gas disk structure, we subsequently perform \textit{N}-body simulations of protoplanets to investigate the dynamical configuration of the planetary system. We find that the protoplanets can grow into planets without significant orbital migration because of the rapid clearing of the inner disk by the magnetically-driven disk winds. Our model can explain the origins of the orbital structure of the inner solar system. Several other features such as the rocky composition can also be explained by the early formation of rocky planetesimals.
△ Less
Submitted 13 October, 2021;
originally announced October 2021.
-
Size evolution of close-in super-Earths through giant impacts and photoevaporation
Authors:
Yuji Matsumoto,
Eiichiro Kokubo,
Pin-Gao Gu,
Kenji Kurosaki
Abstract:
The Kepler transit survey with follow-up spectroscopic observations has discovered numerous super-Earth sized planets and revealed intriguing features of their sizes, orbital periods, and their relations between adjacent planets. For the first time, we investigate the size evolution of planets via both giant impacts and photoevaporation to compare with these observed features. We calculate the siz…
▽ More
The Kepler transit survey with follow-up spectroscopic observations has discovered numerous super-Earth sized planets and revealed intriguing features of their sizes, orbital periods, and their relations between adjacent planets. For the first time, we investigate the size evolution of planets via both giant impacts and photoevaporation to compare with these observed features. We calculate the size of a protoplanet, which is the sum of its core and envelope sizes, by analytical models. $N$-body simulations are performed to evolve planet sizes during the giant impact phase with envelope stripping via impact shocks. We consider the initial radial profile of the core mass and the initial envelope mass fractions as parameters. Inner planets can lose their whole envelopes via giant impacts, while outer planets can keep their initial envelopes since they do not experience giant impacts. Photoevaporation is simulated to evolve planet sizes afterward. Our results suggest that the period-radius distribution of the observed planets would be reproduced if we perform simulations in which the initial radial profile of the core mass follows a wide range of power-law distributions and the initial envelope mass fractions are $\sim0.1$. Moreover, our model shows that the adjacent planetary pairs have similar sizes and regular spacings, with slight differences from detailed observational results such as the radius gap.
△ Less
Submitted 28 September, 2021;
originally announced September 2021.
-
Merging Criteria for Planetesimal Collisions
Authors:
Takashi Shibata,
Eiichiro Kokubo,
Natsuki Hosono
Abstract:
In the standard scenario of planet formation, terrestrial planets, ice giants, and cores of gas giants are formed by the accumulation of planetesimals. However, there are few N-body simulation studies of planetesimal accretion that correctly take into account the merging condition of planetesimals. In order to investigate a realistic accretion process of planetesimals, it is necessary to clarify t…
▽ More
In the standard scenario of planet formation, terrestrial planets, ice giants, and cores of gas giants are formed by the accumulation of planetesimals. However, there are few N-body simulation studies of planetesimal accretion that correctly take into account the merging condition of planetesimals. In order to investigate a realistic accretion process of planetesimals, it is necessary to clarify the merging criteria of planetesimals at collision. We perform numerical collision experiments using smoothed particle hydrodynamics and obtain the merging criteria as a function of planetesimal mass and impact parameters for undifferentiated rocky and icy planetesimals and differentiated icy planetesimals. We vary the total mass of colliding planetesimals, their mass ratios, and the impact angle and obtain the critical impact velocity as the merging criteria distinguishing merging from hit-and-run collision. We find that the critical impact velocity normalized by the two-body surface escape velocity decreases with increasing impact angle. The critical impact velocity does not depend on the total mass, while it has a weak positive dependence on the mass ratio. These results barely depend on the composition and internal structure of the planetesimals.
△ Less
Submitted 19 August, 2021;
originally announced August 2021.
-
Planetesimal Dynamics in the Presence of a Giant Planet
Authors:
Kangrou Guo,
Eiichiro Kokubo
Abstract:
Standard models of planet formation explain how planets form in axisymmetric, unperturbed disks in single star systems. However, it is possible that giant planets could have already formed when other planetary embryos start to grow. We investigate the dynamics of planetesimals under the perturbation of a giant planet in a gaseous disk. Our aim is to understand the effect of the planet's perturbati…
▽ More
Standard models of planet formation explain how planets form in axisymmetric, unperturbed disks in single star systems. However, it is possible that giant planets could have already formed when other planetary embryos start to grow. We investigate the dynamics of planetesimals under the perturbation of a giant planet in a gaseous disk. Our aim is to understand the effect of the planet's perturbation on the formation of giant planet cores outside the orbit of the planet. We calculate the orbital evolution of planetesimals ranging from $10^{13}$ to $10^{20}$g, with a Jupiter-mass planet located at 5.2 au. We find orbital alignment of planetesimals distributed in $\simeq 9$-15 au, except for the mean motion resonance (MMR) locations. The degree of alignment increases with increasing distance from the planet and decreasing planetesimal mass. Aligned orbits lead to low encounter velocity and thus faster growth. The typical velocity dispersion for identical-mass planetesimals is $\sim \mathcal{O}(10)$ $\rm{m\text{ } s}^{-1}$ except for the MMR locations. The relative velocity decreases with increasing distance from the planet and decreasing mass ratio of planetesimals. When the eccentricity vectors of planetesimals reach equilibrium under the gas drag and secular perturbation, the relative velocity becomes lower when the masses of two planetesimals are both on the larger end of the mass spectrum. Our results show that with a giant planet embedded in the disk, the growth of another planetary core outside the planet orbit might be accelerated in certain locations.
△ Less
Submitted 11 June, 2021;
originally announced June 2021.
-
Elementary Process of Galactic Spiral Arm Formation: Phase Synchronization of Epicyclic Motion by Gravitational Scattering
Authors:
Yuki Yoshida,
Eiichiro Kokubo
Abstract:
Swing amplification is a model of spiral arm formation in disk galaxies. Previous $N$-body simulations show that the epicycle phases of stars in spiral arms are synchronized. However, the elementary process of the phase synchronization is not well understood. In order to investigate phase synchronization, we investigate the orbital evolution of stars due to gravitational scattering by a perturber…
▽ More
Swing amplification is a model of spiral arm formation in disk galaxies. Previous $N$-body simulations show that the epicycle phases of stars in spiral arms are synchronized. However, the elementary process of the phase synchronization is not well understood. In order to investigate phase synchronization, we investigate the orbital evolution of stars due to gravitational scattering by a perturber under the epicycle approximation and its dependence on orbital elements and a disk parameter. We find that gravitational scattering by the perturber can cause phase synchronization of stellar orbits. The epicycle phases are better synchronized for smaller initial epicycle amplitudes of stars and larger shear rates of galactic disks. The vertical motion of stars does not affect the phase synchronization. The phase synchronization forms trailing dense regions, which may correspond to spiral arms.
△ Less
Submitted 6 May, 2021; v1 submitted 8 April, 2021;
originally announced April 2021.
-
Formation of "Blanets" from Dust Grains around the Supermassive Black Holes in Galaxies
Authors:
Keiichi Wada,
Yusuke Tsukamoto,
Eiichiro Kokubo
Abstract:
In Wada, Tsukamoto, and Kokubo (2019), we proposed for the first time that a new class of planets, "blanets" (i.e., black hole planets), can be formed around supermassive black holes (SMBHs) in the galactic center. Here, we investigate the dust coagulation processes and physical conditions of the blanet formation outside the snow line ($r_{snow} \sim $ several parsecs) in more detail, especially c…
▽ More
In Wada, Tsukamoto, and Kokubo (2019), we proposed for the first time that a new class of planets, "blanets" (i.e., black hole planets), can be formed around supermassive black holes (SMBHs) in the galactic center. Here, we investigate the dust coagulation processes and physical conditions of the blanet formation outside the snow line ($r_{snow} \sim $ several parsecs) in more detail, especially considering the effect of the radial advection of the dust aggregates. We found that a dimensionless parameter $α= v_t^2/c_s^2$, where $v_t$ is the turbulent velocity and $c_s$ is the sound velocity, describing the turbulent viscosity should be smaller than 0.04 in the circumnuclear disk, to prevent the destruction of the aggregates due to collisions. The formation timescale of blanets $τ_{GI}$ at $r_{snow}$ is, $τ_{GI} \simeq$ 70-80 Myr for $α= 0.01-0.04$ and $M_{BH} = 10^6 M_\odot$. The mass of the blanets ranges from $\sim 20 M_E$ to $3000 M_E$ in $r < 4$ pc for $α= 0.02$ ($M_E$ is the Earth mass), which is in contrast with $ 4 M_E-6 M_E$ for the case without the radial advection. Our results suggest that \textit{blanets} could be formed around relatively low-luminosity AGNs ($L_{bol} \sim 10^{42}$ erg s$^{-1}$) during their lifetime ($\lesssim 10^8$ yr).
△ Less
Submitted 15 December, 2020; v1 submitted 29 July, 2020;
originally announced July 2020.
-
Ejection of close-in super-Earths around low-mass stars in the giant impact stage
Authors:
Yuji Matsumoto,
Pin-Gao Gu,
Eiichiro Kokubo,
Shoichi Oshino,
Masashi Omiya
Abstract:
Earth-sized planets were observed in close-in orbits around M dwarfs. While more and more planets are expected to be uncovered around M dwarfs, theories of their formation and dynamical evolution are still in their infancy. We investigate the giant impact growth of protoplanets, which includes strong scattering around low-mass stars. The aim is to clarify whether strong scattering around low-mass…
▽ More
Earth-sized planets were observed in close-in orbits around M dwarfs. While more and more planets are expected to be uncovered around M dwarfs, theories of their formation and dynamical evolution are still in their infancy. We investigate the giant impact growth of protoplanets, which includes strong scattering around low-mass stars. The aim is to clarify whether strong scattering around low-mass stars affects the orbital and mass distributions of the planets. We perform $N$-body simulation of protoplanets by systematically surveying the parameter space of the stellar mass and surface density of protoplanets. We find that protoplanets are often ejected after twice or three times close-scattering around late M dwarfs. The ejection sets the upper limit of the largest planet mass. Adopting the surface density scaling linearly with the stellar mass, we find that as the stellar mass decreases less massive planets are formed in orbits with higher eccentricities and inclinations. Under this scaling, we also find that a few close-in protoplanets are generally ejected. The ejection of protoplanets plays an important role in the mass distribution of super-Earths around late M dwarfs. The mass relation of observed close-in super-Earths and their central star mass is well reproduced by ejection.
△ Less
Submitted 28 July, 2020;
originally announced July 2020.
-
Coherent Stellar Motion in Galactic Spiral Arms by Swing Amplification
Authors:
Shugo Michikoshi,
Eiichiro Kokubo
Abstract:
We perform local $N$-body simulations of disk galaxies and investigate the evolution of spiral arms. We calculate the time autocorrelation of the surface density of spiral arms and find that the typical evolution timescale is described by the epicycle period. We investigate the distribution of the orbital elements of stars and find that in spiral arms the epicycle motions of stars are in phase whi…
▽ More
We perform local $N$-body simulations of disk galaxies and investigate the evolution of spiral arms. We calculate the time autocorrelation of the surface density of spiral arms and find that the typical evolution timescale is described by the epicycle period. We investigate the distribution of the orbital elements of stars and find that in spiral arms the epicycle motions of stars are in phase while the spatial distribution of the guiding center is nearly uniform. These facts clearly show that the phase synchronization of the epicycle motion takes place, which is theoretically predicted by the swing amplification.
△ Less
Submitted 15 May, 2020;
originally announced May 2020.
-
Hyperbolic Orbits in the Solar System: Interstellar Origin or Perturbed Oort Cloud Comets?
Authors:
Arika Higuchi,
Eiichiro Kokubo
Abstract:
We study the dynamical properties of objects in hyperbolic orbits passing through the inner Solar system in the context of two different potential sources: interstellar space and the Oort cloud. We analytically derive the probability distributions of eccentricity, $e$, and perihelion distance, $q$, for each source and estimate the numbers of objects produced per unit of time as a function of these…
▽ More
We study the dynamical properties of objects in hyperbolic orbits passing through the inner Solar system in the context of two different potential sources: interstellar space and the Oort cloud. We analytically derive the probability distributions of eccentricity, $e$, and perihelion distance, $q$, for each source and estimate the numbers of objects produced per unit of time as a function of these quantities. By comparing the numbers from the two sources, we assess which origin is more likely for a hyperbolic object having a given eccentricity and perihelion distance. We find that the likelihood that a given hyperbolic object is of interstellar origin increases with decreasing eccentricity and perihelion. Conversely, the likelihood that a hyperbolic object has been scattered from the Oort cloud by a passing star increases with decreasing eccentricity and increasing perihelion. By carefully considering their orbital elements, we conclude that both 1I/2017 U1 'Oumuamua ($e\simeq$ 1.2 and $q\simeq$ 0.26 au) and 2I/2019 Q4 Borisov ($e\simeq$ 3.3 and $q\simeq$ 2 au) are most likely of interstellar origin, not scattered from the Oort cloud. However, we also find that Oort cloud objects can be scattered into hyperbolic orbits like those of the two known examples, by sub-stellar and even sub-Jovian mass perturbers. This highlights the need for better characterization of the low mass end of the free-floating brown dwarf and planet population.
△ Less
Submitted 11 November, 2019;
originally announced November 2019.
-
Planet Formation around Super Massive Black Holes in the Active Galactic Nuclei
Authors:
Keiichi Wada,
Yusuke Tsukamoto,
Eiichiro Kokubo
Abstract:
As a natural consequence of the elementary processes of dust growth, we discovered that a new class of planets can be formed around supermassive black holes (SMBHs). We investigated a growth path from sub-micron sized icy dust monomers to Earth-sized bodies outside the "snow line", located several parsecs from SMBHs in low luminosity active galactic nuclei (AGNs). In contrast to protoplanetary dis…
▽ More
As a natural consequence of the elementary processes of dust growth, we discovered that a new class of planets can be formed around supermassive black holes (SMBHs). We investigated a growth path from sub-micron sized icy dust monomers to Earth-sized bodies outside the "snow line", located several parsecs from SMBHs in low luminosity active galactic nuclei (AGNs). In contrast to protoplanetary disks, the "radial drift barrier" does not prevent the formation of planetesimals. In the early phase of the evolution, low collision velocity between dust particles promotes sticking; therefore, the internal density of the dust aggregates decreases with growth. When the porous aggregate's size reaches 0.1--1 cm, the collisional compression becomes effective, and the decrease in internal density stops. Once 10--100 m sized aggregates are formed, they are decoupled from gas turbulence, and the aggregate layer becomes gravitationally unstable, leading to the formation of planets by the fragmentation of the layer, with ten times the mass of the earth. The growth time scale depends on the turbulent strength of the circumnuclear disk and the black hole mass $M_{BH}$, and it is comparable to the AGN's lifetime ($\sim 10^8$ yr) for low mass ($M_{BH} \sim 10^6 M_\odot$) SMBHs.
△ Less
Submitted 9 October, 2019; v1 submitted 15 September, 2019;
originally announced September 2019.
-
Global N-Body Simulation of Galactic Spiral Arms
Authors:
Shugo Michikoshi,
Eiichiro Kokubo
Abstract:
The origin of galactic spiral arms is one of fundamental problems in astrophysics. Based on the local analysis Toomre (1981) proposed the swing amplification mechanism in which the self-gravity forms spiral arms as leading waves of stars rotate to trailing ones due to galactic shear. The structure of spiral arms is characterized by their number and pitch angle. We perform global $N$-body simulatio…
▽ More
The origin of galactic spiral arms is one of fundamental problems in astrophysics. Based on the local analysis Toomre (1981) proposed the swing amplification mechanism in which the self-gravity forms spiral arms as leading waves of stars rotate to trailing ones due to galactic shear. The structure of spiral arms is characterized by their number and pitch angle. We perform global $N$-body simulations of spiral galaxies to investigate the dependence of the spiral structure on disk parameters and compare the simulation results with the swing amplification model. We find that the spiral structure in the $N$-body simulations agrees well with that predicted by the swing amplification for the wide range of parameters. The pitch angle decreases with increasing the shear rate and is independent of the disk mass fraction. The number of spiral arms decreases with both increasing the shear rate and the disk mass fraction. If the disk mass fraction is fixed, the pitch angle increases with the number of spiral arms.
△ Less
Submitted 24 August, 2018;
originally announced August 2018.
-
Planets around the evolved stars 24 Booties and $γ$ Libra: A 30d-period planet and a double giant-planet system in possible 7:3 MMR
Authors:
Takuya Takarada,
Bun'ei Sato,
Masashi Omiya,
Hiroki Harakawa,
Makiko Nagasawa,
Hideyuki Izumiura,
Eiji Kambe,
Yoichi Takeda,
Michitoshi Yoshida,
Yoichi Itoh,
Hiroyasu Ando,
Eiichiro Kokubo,
Shigeru Ida
Abstract:
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 shor…
▽ More
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. $γ$ 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 $γ$ 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.
△ Less
Submitted 11 April, 2018;
originally announced April 2018.
-
Formation of the terrestrial planets in the solar system around 1 au via radial concentration of planetesimals
Authors:
Masahiro Ogihara,
Eiichiro Kokubo,
Takeru K. Suzuki,
Alessandro Morbidelli
Abstract:
No planets exist inside the orbit of Mercury and the terrestrial planets of the solar system exhibit a localized configuration. According to thermal structure calculation of protoplanetary disks, a silicate condensation line (~ 1300 K) is located around 0.1 au from the Sun except for the early phase of disk evolution, and planetesimals could have formed inside the orbit of Mercury. A recent study…
▽ More
No planets exist inside the orbit of Mercury and the terrestrial planets of the solar system exhibit a localized configuration. According to thermal structure calculation of protoplanetary disks, a silicate condensation line (~ 1300 K) is located around 0.1 au from the Sun except for the early phase of disk evolution, and planetesimals could have formed inside the orbit of Mercury. A recent study of disk evolution that includes magnetically driven disk winds showed that the gas disk obtains a positive surface density slope inside ~ 1 au from the central star. In a region with positive midplane pressure gradient, planetesimals undergo outward radial drift. We investigate the radial drift of planetesimals and type I migration of planetary embryos in a disk that viscously evolves with magnetically driven disk winds. We show a case in which no planets remain in the close-in region. Radial drifts of planetesimals are simulated using a recent disk evolution model that includes effects of disk winds. The late stage of planet formation is also examined by performing N-body simulations of planetary embryos. We demonstrate that in the middle stage of disk evolution, planetesimals can undergo convergent radial drift in a magnetorotational instability (MRI)-inactive disk, in which the pressure maximum is created, and accumulate in a narrow ring-like region with an inner edge at ~ 0.7 au from the Sun. We also show that planetary embryos that may grow from the narrow planetesimal ring do not exhibit significant type I migration in the late stage of disk evolution. The origin of the localized configuration of the terrestrial planets of the solar system, in particular the deficit of close-in planets, can be explained by the convergent radial drift of planetesimals in disks with a positive pressure gradient in the close-in region.
△ Less
Submitted 6 April, 2018;
originally announced April 2018.
-
Formation of close-in super-Earths in evolving protoplanetary disks due to disk winds
Authors:
Masahiro Ogihara,
Eiichiro Kokubo,
Takeru K. Suzuki,
Alessandro Morbidelli
Abstract:
Planets with masses larger than about 0.1 Earth-masses undergo rapid inward migration (type I migration) in a standard protoplanetary disk. Recent magnetohydrodynamical simulations revealed the presence of magnetically driven disk winds, which would alter the disk profile and the type I migration in the close-in region. We investigate orbital evolution of planetary embryos in disks that viscously…
▽ More
Planets with masses larger than about 0.1 Earth-masses undergo rapid inward migration (type I migration) in a standard protoplanetary disk. Recent magnetohydrodynamical simulations revealed the presence of magnetically driven disk winds, which would alter the disk profile and the type I migration in the close-in region. We investigate orbital evolution of planetary embryos in disks that viscously evolve under the effects of disk winds. The aim is to discuss effects of altered disk profiles on type I migration. In addition, we aim to examine whether observed distributions of close-in super-Earths can be reproduced by simulations that include effects of disk winds. We perform N-body simulations of super-Earth formation from planetary embryos, in which a recent model for disk evolution is used. We explore a wide range of parameters and draw general trends. We also carry out N-body simulations of close-in super-Earth formation from embryos in such disks under various conditions. We find that the type I migration is significantly suppressed in many cases. Even in cases in which inward migration occurs, the migration timescale is lengthened to 1 Myr, which mitigates the type I migration problem. This is because the gas surface density is decreased and has a flatter profile in the close-in region due to disk winds. We find that when the type I migration is significantly suppressed, planets undergo late orbital instability during the gas depletion, leading to a non-resonant configuration. We also find that observed distributions of close-in super-Earths (e.g., period ratio, mass ratio) can be reproduced. In addition, we show that in some results of simulations, systems with a chain of resonant planets, like the TRAPPIST-1 system, form.
△ Less
Submitted 3 April, 2018;
originally announced April 2018.
-
Gravitational instability of a dust layer composed of porous silicate dust aggregates in a protoplanetary disk
Authors:
Misako Tatsuuma,
Shugo Michikoshi,
Eiichiro Kokubo
Abstract:
Planetesimal formation is one of the most important unsolved problems in planet formation theory. In particular, rocky planetesimal formation is difficult because silicate dust grains are easily broken when they collide. Recently, it has been proposed that they can grow as porous aggregates when their monomer radius is smaller than $\sim$ 10 nm, which can also avoid the radial drift toward the cen…
▽ More
Planetesimal formation is one of the most important unsolved problems in planet formation theory. In particular, rocky planetesimal formation is difficult because silicate dust grains are easily broken when they collide. Recently, it has been proposed that they can grow as porous aggregates when their monomer radius is smaller than $\sim$ 10 nm, which can also avoid the radial drift toward the central star. However, the stability of a layer composed of such porous silicate dust aggregates has not been investigated. Therefore, we investigate the gravitational instability of this dust layer. To evaluate the disk stability, we calculate Toomre's stability parameter $Q$, for which we need to evaluate the equilibrium random velocity of dust aggregates. We calculate the equilibrium random velocity considering gravitational scattering and collisions between dust aggregates, drag by mean flow of gas, stirring by gas turbulence, and gravitational scattering by gas density fluctuation due to turbulence. We derive the condition of the gravitational instability using the disk mass, dust-to-gas ratio, turbulent strength, orbital radius, and dust monomer radius. We find that, for the minimum mass solar nebula model at 1 au, the dust layer becomes gravitationally unstable when the turbulent strength $α\lesssim10^{-5}$. If the dust-to-gas ratio is increased twice, the gravitational instability occurs for $α\lesssim10^{-4}$. We also find that the dust layer is more unstable in disks with larger mass, higher dust-to-gas ratio, and weaker turbulent strength, at larger orbital radius, and with a larger monomer radius.
△ Less
Submitted 8 February, 2018;
originally announced February 2018.
-
Effects of global gas flows on type I migration
Authors:
Masahiro Ogihara,
Eiichiro Kokubo,
Takeru K. Suzuki,
Alessandro Morbidelli,
Aurélien Crida
Abstract:
Magnetically-driven disk winds would alter the surface density slope of gas in the inner region of a protoplanetary disk $(r \lesssim 1 {\rm au})$. This in turn affects planet formation. Recently, the effect of disk wind torque has been considered with the suggestion that it would carve out the surface density of the disk from inside and would induce global gas flows (wind-driven accretion). We ai…
▽ More
Magnetically-driven disk winds would alter the surface density slope of gas in the inner region of a protoplanetary disk $(r \lesssim 1 {\rm au})$. This in turn affects planet formation. Recently, the effect of disk wind torque has been considered with the suggestion that it would carve out the surface density of the disk from inside and would induce global gas flows (wind-driven accretion). We aim to investigate effects of global gas flows on type I migration and also examine planet formation. A simplified approach was taken to address this issue, and N-body simulations with isolation-mass planets were also performed. In previous studies, the effect of gas flow induced by turbulence-driven accretion has been taken into account for its desaturation effect of the corotation torque. If more rapid gas flows (e.g., wind-driven accretion) are considered, the desaturation effect can be modified. In MRI-inactive disks, in which the wind-driven accretion dominates the disk evolution, the gas flow at the midplane plays an important role. If this flow is fast, the corotation torque is efficiently desaturated. Then, the fact that the surface density slope can be positive in the inner region due to the wind torque can generate an outward migration region extended to super-Earth mass planets. In this case, we observe that no planets fall onto the central star in N-body simulations with migration forces imposed to reproduce such migration pattern. We also see that super-Earth mass planets can undergo outward migration. Relatively rapid gas flows affects type I migration and thus the formation of close-in planets.
△ Less
Submitted 3 October, 2017;
originally announced October 2017.
-
Formation of Close-in Super-Earths by Giant Impacts: Effects of Initial Eccentricities and Inclinations of Protoplanets
Authors:
Yuji Matsumoto,
Eiichiro Kokubo
Abstract:
Recent observations have revealed the eccentricity and inclination distributions of close-in super-Earths. These distributions have the potential to constrain their formation processes. In the in-situ formation scenario, the eccentricities and inclinations of planets are determined by gravitational scattering and collisions between protoplanets on the giant impact stage. We investigate the effect…
▽ More
Recent observations have revealed the eccentricity and inclination distributions of close-in super-Earths. These distributions have the potential to constrain their formation processes. In the in-situ formation scenario, the eccentricities and inclinations of planets are determined by gravitational scattering and collisions between protoplanets on the giant impact stage. We investigate the effect of the initial eccentricities and inclinations of protoplanets on the formation of close-in super-Earths. We perform $N$-body simulations of protoplanets in gas-free disks, changing the initial eccentricities and inclinations systematically. We find that while the eccentricities of protoplanets are well relaxed through their evolution, the inclinations are not. When the initial inclinations are small, they are not generally pumped up since scattering is less effective and collisions occur immediately after orbital crossing. On the other hand, when the initial inclinations are large, they tend to be kept large since collisional damping is less effective. Not only the resultant inclinations of planets, but also their number, eccentricities, angular momentum deficit, and orbital separations are affected by the initial inclinations of protoplanets.
△ Less
Submitted 22 May, 2017;
originally announced May 2017.
-
Dynamics of Porous Dust Aggregates and Gravitational Instability of Their Disk
Authors:
Shugo Michikoshi,
Eiichiro Kokubo
Abstract:
We consider the dynamics of porous icy dust aggregates in a turbulent gas disk and investigate the stability of the disk. We evaluate the random velocity of porous dust aggregates by considering their self-gravity, collisions, aerodynamic drag, turbulent stirring and scattering due to gas. We extend our previous work by introducing the anisotropic velocity dispersion and the relaxation time of the…
▽ More
We consider the dynamics of porous icy dust aggregates in a turbulent gas disk and investigate the stability of the disk. We evaluate the random velocity of porous dust aggregates by considering their self-gravity, collisions, aerodynamic drag, turbulent stirring and scattering due to gas. We extend our previous work by introducing the anisotropic velocity dispersion and the relaxation time of the random velocity. We find the minimum mass solar nebular model to be gravitationally unstable if the turbulent viscosity parameter $α$ is less than about $4 \times 10^{-3}$. The upper limit of $α$ for the onset of gravitational instability is derived as a function of the disk parameters. We discuss the implications of the gravitational instability for planetesimal formation.
△ Less
Submitted 12 May, 2017;
originally announced May 2017.
-
Simulating the Smallest Ring World of Chariklo
Authors:
Shugo Michikoshi,
Eiichiro Kokubo
Abstract:
A ring system consisting of two dense narrow rings has been discovered around Centaur Chariklo. The existence of these rings around a small object poses various questions, such as their origin, stability, and lifetime. In order to understand the nature of Chariklo's rings, we perform global $N$-body simulations of the self-gravitating collisional particle rings for the first time. We find that Cha…
▽ More
A ring system consisting of two dense narrow rings has been discovered around Centaur Chariklo. The existence of these rings around a small object poses various questions, such as their origin, stability, and lifetime. In order to understand the nature of Chariklo's rings, we perform global $N$-body simulations of the self-gravitating collisional particle rings for the first time. We find that Chariklo should be denser than the ring material to avoid the rapid diffusion of the rings. If Chariklo is denser than the ring material, fine spiral structures called self-gravity wakes occur in the inner ring. These wakes accelerate the viscous spreading of the ring significantly and they typically occur on timescales of about $100\,\mathrm{years}$ for m-sized ring particles, which is considerably shorter than the timescales suggested in previous studies. The existence of these narrow rings implies smaller ring particles or the existence of shepherding satellites.
△ Less
Submitted 21 February, 2017;
originally announced February 2017.
-
Planetesimal Formation by Gravitational Instability of a Porous-Dust Disk
Authors:
Shugo Michikoshi,
Eiichiro Kokubo
Abstract:
Recently it is proposed that porous icy dust aggregates are formed by pairwise accretion of dust aggregates beyond the snowline. We calculate the equilibrium random velocity of porous dust aggregates taking into account mutual gravitational scattering, collisions, gas drag, and turbulent stirring and scattering. We find that the disk of porous dust aggregates becomes gravitationally unstable as th…
▽ More
Recently it is proposed that porous icy dust aggregates are formed by pairwise accretion of dust aggregates beyond the snowline. We calculate the equilibrium random velocity of porous dust aggregates taking into account mutual gravitational scattering, collisions, gas drag, and turbulent stirring and scattering. We find that the disk of porous dust aggregates becomes gravitationally unstable as they evolve through gravitational compression in the minimum-mass solar nebula model for a reasonable range of turbulence strength, which leads to rapid formation of planetesimals.
△ Less
Submitted 22 June, 2016;
originally announced June 2016.
-
Swing Amplification of Galactic Spiral Arms: Phase Synchronization of Stellar Epicycle Motion
Authors:
Shugo Michikoshi,
Eiichiro Kokubo
Abstract:
We revisit the swing amplification model of galactic spiral arms proposed by Toomre (1981). We describe the derivation of the perturbation equation in detail and investigate the amplification process of stellar spirals. We find that the elementary process of the swing amplification is the phase synchronization of the stellar epicycle motion. Regardless of the initial epicycle phase, the epicycle p…
▽ More
We revisit the swing amplification model of galactic spiral arms proposed by Toomre (1981). We describe the derivation of the perturbation equation in detail and investigate the amplification process of stellar spirals. We find that the elementary process of the swing amplification is the phase synchronization of the stellar epicycle motion. Regardless of the initial epicycle phase, the epicycle phases of stars in a spiral are synchronized during the amplification. Based on the phase synchronization, we explain the dependence of the pitch angle of spirals on the epicycle frequency. We find the most amplified spiral mode and calculate its pitch angle, wavelengths, and amplification factor, which are consistent with those obtained by the more rigorous model based on the Boltzmann equation by Julian and Toomre (1966).
△ Less
Submitted 11 April, 2016;
originally announced April 2016.
-
High-Contrast Imaging of Intermediate-Mass Giants with Long-Term Radial Velocity Trends
Authors:
Tsuguru Ryu,
Bun'ei Sato,
Masayuki Kuzuhara,
Norio Narita,
Yasuhiro H. Takahashi,
Taichi Uyama,
Tomoyuki Kudo,
Nobuhiko Kusakabe,
Jun Hashimoto,
Masashi Omiya,
Hiroki Harakawa,
Lyu Abe,
Hiroyasu Ando,
Wolfgang Brandner,
Timothy D. Brandt,
Joseph C. Carson,
Thayne Currie,
Sebastian Egner,
Markus Feldt,
Miwa Goto,
Carol A. Grady,
Olivier Guyon,
Yutaka Hayano,
Masahiko Hayashi,
Saeko S. Hayashi
, et al. (40 additional authors not shown)
Abstract:
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-cont…
▽ More
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.
△ Less
Submitted 4 April, 2016; v1 submitted 7 March, 2016;
originally announced March 2016.
-
Galactic Spiral Arms by Swing Amplification
Authors:
Shugo Michikoshi,
Eiichiro Kokubo
Abstract:
Based on the swing amplification model of Julian and Toomre (1966), we investigate the formation and structure of stellar spirals in disk galaxies. We calculate the pitch angle, wavelengths, and amplification factor of the most amplified mode. We also obtain the fitting formulae of these quantities as a function of the epicycle frequency and Toomre's $Q$. As the epicycle frequency increases, the p…
▽ More
Based on the swing amplification model of Julian and Toomre (1966), we investigate the formation and structure of stellar spirals in disk galaxies. We calculate the pitch angle, wavelengths, and amplification factor of the most amplified mode. We also obtain the fitting formulae of these quantities as a function of the epicycle frequency and Toomre's $Q$. As the epicycle frequency increases, the pitch angle and radial wavelength increases, while the azimuthal wavelength decreases. The pitch angle and radial wavelength increases with $Q$, while the azimuthal wavelength weakly depends on $Q$. The amplification factor decreases with $Q$ rapidly. In order to confirm the swing amplification model, we perform local $N$-body simulations. The wavelengths and pitch angle by the swing amplification model are in good agreement with those by $N$-body simulations. The dependence of the amplification factor on the epicycle frequency in $N$-body simulations is generally consistent with that in the swing amplification model. Using these results, we estimate the number of spiral arms as a function of the shear rate. The number of spiral arms increases with the shear rate if the disk to halo mass ratio is fixed.
△ Less
Submitted 2 March, 2016;
originally announced March 2016.
-
A Pair of Giant Planets around the Evolved Intermediate-Mass Star HD 47366: Multiple Circular Orbits or a Mutually Retrograde Configuration
Authors:
Bun'ei Sato,
Liang Wang,
Yu-Juan Liu,
Gang Zhao,
Masashi Omiya,
Hiroki Harakawa,
Makiko Nagasawa,
Robert A. Wittenmyer,
Paul Butler,
Nan Song,
Wei He,
Fei Zhao,
Eiji Kambe,
Kunio Noguchi,
Hiroyasu Ando,
Hideyuki Izumiura,
Norio Okada,
Michitoshi Yoshida,
Yoichi Takeda,
Yoichi Itoh,
Eiichiro Kokubo,
Shigeru Ida
Abstract:
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 plan…
▽ More
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.
△ Less
Submitted 18 January, 2016;
originally announced January 2016.
-
Dynamics of Self-Gravity Wakes in Dense Planetary Rings I. Pitch Angle
Authors:
Shugo Michikoshi,
Akihiko Fujii,
Eiichiro Kokubo,
Heikki Salo
Abstract:
We investigate the dynamics of self-gravity wakes in dense planetary rings. In particular, we examine how the pitch angle of self-gravity wakes depend on ring parameters using N-body simulations. We calculate the pitch angles using the two-dimensional autocorrelation function of the ring surface density. We obtain the pitch angles for the inner and outer parts of the autocorrelation function separ…
▽ More
We investigate the dynamics of self-gravity wakes in dense planetary rings. In particular, we examine how the pitch angle of self-gravity wakes depend on ring parameters using N-body simulations. We calculate the pitch angles using the two-dimensional autocorrelation function of the ring surface density. We obtain the pitch angles for the inner and outer parts of the autocorrelation function separately. We confirm that the pitch angles are 15 to 30 degrees for reasonable ring parameters, which are consistent with previous studies. We find that the inner pitch angle increases with the Saturnicentric distance, while it barely depends on the optical depth and the restitution coefficient of ring particles. The increase of the inner pitch angle with the Saturnicentric distance is consistent with the observations of the A ring. The outer pitch angle does not have the clear dependence on any ring parameters and is about 10 - 15 degrees. This value is consistent with the pitch angle of spiral arms in collisionless systems.
△ Less
Submitted 13 September, 2015;
originally announced September 2015.
-
Warm Debris Disks Produced by Giant Impacts During Terrestrial Planet Formation
Authors:
H. Genda,
H. Kobayashi,
E. Kokubo
Abstract:
In our solar system, Mars-sized protoplanets frequently collided with each other during the last stage of terrestrial planet formation called the giant impact stage. Giant impacts eject a large amount of material from the colliding protoplanets into the terrestrial planet region, which may form debris disks with observable infrared excesses. Indeed, tens of warm debris disks around young solar-typ…
▽ More
In our solar system, Mars-sized protoplanets frequently collided with each other during the last stage of terrestrial planet formation called the giant impact stage. Giant impacts eject a large amount of material from the colliding protoplanets into the terrestrial planet region, which may form debris disks with observable infrared excesses. Indeed, tens of warm debris disks around young solar-type stars have been observed. Here, we quantitatively estimate the total mass of ejected materials during the giant impact stages. We found that $\sim$0.4 times the Earth's mass is ejected in total throughout the giant impact stage. Ejected materials are ground down by collisional cascade until micron-sized grains are blown out by radiation pressure. The depletion timescale of these ejected materials is determined primarily by the mass of the largest body among them. We conducted high-resolution simulations of giant impacts to accurately obtain the mass of the largest ejected body. We then calculated the evolution of the debris disks produced by a series of giant impacts and depleted by collisional cascades to obtain the infrared excess evolution of the debris disks. We found that the infrared excess is almost always higher than the stellar infrared flux throughout the giant impact stage ($\sim$100 Myr) and is sometimes $\sim$10 times higher immediately after a giant impact. Therefore, giant impact stages would explain the infrared excess from most observed warm debris disks. The observed fraction of stars with warm debris disks indicates that the formation probability of our solar system-like terrestrial planets is approximately 10%.
△ Less
Submitted 5 August, 2015;
originally announced August 2015.
-
Effect of Stellar Encounters on Comet Cloud Formation
Authors:
Arika Higuchi,
Eiichiro Kokubo
Abstract:
We have investigated the effect of stellar encounters on the formation and disruption of the Oort cloud using the classical impulse approximation. We calculate the evolution of a planetesimal disk into a spherical Oort cloud due to the perturbation from passing stars for 10 Gyr. We obtain the empirical fits of the $e$-folding time for the number of Oort cloud comets using the standard exponential…
▽ More
We have investigated the effect of stellar encounters on the formation and disruption of the Oort cloud using the classical impulse approximation. We calculate the evolution of a planetesimal disk into a spherical Oort cloud due to the perturbation from passing stars for 10 Gyr. We obtain the empirical fits of the $e$-folding time for the number of Oort cloud comets using the standard exponential and Kohlrausch formulae as functions of the stellar parameters and the initial semimajor axes of planetesimals. The $e$-folding time and the evolution timescales of the orbital elements are also analytically derived. In some calculations, the effect of the Galactic tide is additionally considered. We also show the radial variations of the $e$-folding times to the Oort cloud. From these timescales, we show that if the initial planetesimal disk has the semimajor axes distribution ${\rm d}n/{\rm d}a\propto a^{-2}$, which is produced by planetary scattering (Higuchi et al. 2006), the $e$-folding time for planetesimals in the Oort cloud is $\sim$10 Gyr at any heliocentric distance $r$. This uniform $e$-folding time over the Oort cloud means that the supply of comets from the inner Oort cloud to the outer Oort cloud is sufficiently effective to keep the comet distribution as ${\rm d}n/{\rm d}r\propto r^{-2}$. We also show that the final distribution of the semimajor axes in the Oort cloud is approximately proportional to $a^{-2}$ for any initial distribution.
△ Less
Submitted 2 July, 2015;
originally announced July 2015.
