Researcher profile

Eiichiro Kokubo

Eiichiro Kokubo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 published item(s)

preprint2026arXiv

Unraveling Year-Long Radial Velocity Variations in Red Clump Region -- I: Comprehensive analysis of a K0 Giant star, 2 Draconis

Slow-rotating evolved stars frequently exhibit radial velocity (RV) variations on annual timescales, complicated by instrumental systematics and aliasing in the one-year regime. Here we investigate the origin of the near-yearly periodicity in 2 Dra, a star located in the red-clump region, assessing possible causes between stellar activity, instrumental profile (IP) effects, sampling alias, and planetary companions. We applied two independent approaches: (1) constraining diagnostic signals and performing a correlation analysis ($r$) between period-confined signals, and (2) evaluating phase stability by partitioning Keplerian fits. These methods enabled us to examine the physical connections and phase coherence among stellar activity indicators, RV measurements, and IP diagnostics. Our analysis suggests a stellar rotation period of $\simeq270\text{--}320$\,d for 2~Dra. The 340-d RV signal does not appear to originate from stellar activity in this chromospherically quiet star ($|r| \lesssim 0.33$), nor from instrumental systematics near the annual period ($|r| \lesssim 0.1$). This conclusion is supported by contrasting phase behavior: the RV and stellar activity phases remain stable, whereas the IP phases do not. We therefore propose that the 340-d variation likely arises from either small-amplitude intrinsic variability or a tentative gas giant companion with potential weak activity-induced modulation. The case of 2~Dra provides a framework for distinguishing the origins of $\sim$1-yr RV variations in other evolved stars.

preprint2022arXiv

A Super-Earth Orbiting Near the Inner Edge of the Habitable Zone around the M4.5-dwarf Ross 508

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.

preprint2022arXiv

Impacts of viscous dissipation on collisional growth and fragmentation of dust aggregates

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.

preprint2022arXiv

Orbital structure of planetary systems formed by giant impacts: stellar mass dependence

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.

preprint2022arXiv

Planetesimal Dynamics in the Presence of a Giant Planet II: Dependence on Planet Mass and Eccentricity

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.

preprint2022arXiv

Rapid-then-slow migration reproduces mass distribution of TRAPPIST-1 system

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.

preprint2021arXiv

Elemental abundances of nearby M dwarfs based on high-resolution near-infrared spectra obtained by the Subaru/IRD survey: Proof of concept

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.

preprint2021arXiv

Regular Radial Velocity Variations in Nine G- and K-type Giant Stars: Eight Planets and One Planet Candidate

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 &#34;planet desert&#34; 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.

preprint2020arXiv

Coherent Stellar Motion in Galactic Spiral Arms by Swing Amplification

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.

preprint2019arXiv

Hyperbolic Orbits in the Solar System: Interstellar Origin or Perturbed Oort Cloud Comets?

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 &#39;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.

preprint2019arXiv

Planet Formation around Super Massive Black Holes in the Active Galactic Nuclei

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 &#34;snow line&#34;, located several parsecs from SMBHs in low luminosity active galactic nuclei (AGNs). In contrast to protoplanetary disks, the &#34;radial drift barrier&#34; 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&#39;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&#39;s lifetime ($\sim 10^8$ yr) for low mass ($M_{BH} \sim 10^6 M_\odot$) SMBHs.