BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Eiichiro Kokubo

Eiichiro Kokubo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.GAastro-phastro-ph.COcond-mat.soft

Catalog footprint

What is connected

41works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

41 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 "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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 '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.

0 citations0 reviewsNo code linkedOpen access
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 "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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galactic Spiral Arms by Swing Amplification

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Planetesimal Formation by Gravitational Instability of a Porous-Dust Disk

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Swing Amplification of Galactic Spiral Arms: Phase Synchronization of Stellar Epicycle Motion

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).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamics of Self-Gravity Wakes in Dense Planetary Rings I. Pitch Angle

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effect of Stellar Encounters on Comet Cloud Formation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pitch Angle of Galactic Spiral Arms

One of the key parameters that characterize spiral arms in disk galaxies is a pitch angle that measures the inclination of a spiral arm to the direction of galactic rotation. The pitch angle differs from galaxy to galaxy, which suggests that the rotation law of galactic disks determines it. In order to investigate the relation between the pitch angle of spiral arms and the shear rate of galactic differential rotation, we perform local $N$-body simulations of pure stellar disks. We find that the pitch angle increases with the epicycle frequency and decreases with the shear rate and obtain the fitting formula. This dependence is explained by the swing amplification mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Terrestrial Planet Formation at Home and Abroad

We review the state of the field of terrestrial planet formation with the goal of understanding the formation of the inner Solar System and low-mass exoplanets. We review the dynamics and timescales of accretion from planetesimals to planetary embryos and from embryos to terrestrial planets. We discuss radial mixing and water delivery, planetary spins and the importance of parameters regarding the disk and embryo properties. Next, we connect accretion models to exoplanets. We first explain why the observed hot Super Earths probably formed by in situ accretion or inward migration. We show how terrestrial planet formation is altered in systems with gas giants by the mechanisms of giant planet migration and dynamical instabilities. Standard models of terrestrial accretion fail to reproduce the inner Solar System. The "Grand Tack" model solves this problem using ideas first developed to explain the giant exoplanets. Finally, we discuss whether most terrestrial planet systems form in the same way as ours, and highlight the key ingredients missing in the current generation of simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamics and Accretion of Planetesimals

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Substellar Companions to Seven Evolved Intermediate-Mass Stars

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tidal disruption of satellites and formation of narrow rings

In this paper we investigate the formation of narrow planetary rings such as those found around Uranus and Saturn through the tidal disruption of a weak, gravitationally bound satellite that migrates within its Roche limit. Using $N$-body simulations, we study the behaviour of rubble piles placed on circular orbits at different distances from a central planet. We consider both homogeneous satellites and differentiated bodies containing a denser core. We show that the Roche limit for a rubble pile is closer to the planet than for a fluid body of the same mean density. The Roche limit for a differentiated body is also closer to the planet than for a homogeneous satellite of the same mean density. Within its Roche limit, a homogeneous satellite totally disrupts and forms a narrow ring. The initial stages of the disruption are similar to the evolution of a viscous fluid ellipsoid, which can be computed semi-analytically. On the other hand, when a differentiated satellite is just within the Roche limit only the mantle is disrupted. This process is similar to Roche-lobe overflow in interacting binary stars and produces two narrow rings on either side of a remnant satellite. We argue that the Uranian rings, and possibly their shepherd satellites, could have been formed through the tidal disruption of a number of protomoons that were formed inside the corotation radius of Uranus and migrated slowly inwards as a result of tidal interaction with the planet.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Possible Substellar Companion to the Intermediate-mass Giant HD 175679

We report the discovery of a substellar companion around the intermediatemass giant HD 175679. Precise radial velocity data of the star from Xinglong Station and Okayama Astrophysical Observatory (OAO) revealed a Keplerian velocity variation with an orbital period of 1366.8 \pm 5.7 days, a semiamplitude of 380.2 \pm 3.2m s.1, and an eccentricity of 0.378 \pm 0.008. Adopting a stellar mass of 2.7 \pm 0.3 M\odot, we obtain the minimum mass of the HD 175679 b is 37.3 \pm 2.8 MJ, and the semimajor axis is 3.36 \pm 0.12 AU. This discovery is the second brown dwarf companion candidate from a joint planet-search program between China and Japan.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Formation of a Propeller Structure by a Moonlet in a Dense Planetary Ring

The Cassini spacecraft discovered a propeller-shaped structure in Saturn's A ring. This propeller structure is thought to be formed by gravitational scattering of ring particles by an unseen embedded moonlet. Self-gravity wakes are prevalent in dense rings due to gravitational instability. Strong gravitational wakes affect the propeller structure. Here, we derive the condition for formation of a propeller structure by a moonlet embedded in a dense ring with gravitational wakes. We find that a propeller structure is formed when the wavelength of the gravitational wakes is smaller than the Hill radius of the moonlet. We confirm this formation condition by performing numerical simulations. This condition is consistent with observations of propeller structures in Saturn's A ring.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Origin of multiple nuclei in ultraluminous infrared galaxies

Ultraluminous infrared galaxies (ULIRGs) with multiple ($\ge 3$) nuclei are frequently observed. It has been suggested that these nuclei are produced by multiple major mergers of galaxies. The expected rate of such mergers is, however, too low to reproduce the observed number of ULIRGs with multiple nuclei. We have performed high-resolution simulations of the merging of two gas-rich disk galaxies. We found that extremely massive and compact star clusters form from the strongly disturbed gas disks after the first or second encounter between the galaxies. The mass of such clusters reaches $\sim 10^8 M_{\odot}$, and their half-mass radii are $20-30 \rm{pc}$. Since these clusters consist of young stars, they appear to be several bright cores in the galactic central region ($\sim \rm{kpc}$). The peak luminosity of these clusters reaches $\sim 10%$ of the total luminosity of the merging galaxy. These massive and compact clusters are consistent with the characteristics of the observed multiple nuclei in ULIRGs. Multiple mergers are not necessary to explain multiple nuclei in ULIRGs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Secular Gravitational Instability of a Dust Layer in Shear Turbulence

We perform a linear stability analysis of a dust layer in a turbulent gas disk. Youdin (2011) investigated the secular gravitational instability of a dust layer using hydrodynamic equations with a turbulent diffusion term. We obtain essentially the same result independently of Youdin (2011). In the present analysis, we restrict the area of interest to small dust particles, while investigating the secular gravitational instability in a more rigorous manner. We discuss the time evolution of the dust surface density distribution using a stochastic model and derive the advection-diffusion equation. The validity of the analysis by Youdin (2011) is confirmed in the strong drag limit. We demonstrate quantitatively that the finite thickness of a dust layer weakens the secular gravitational instability and that the density-dependent diffusion coefficient changes the growth rate. We apply the obtained results to the turbulence driven by the shear instability and find that the secular gravitational instability is faster than the radial drift when the gas density is three times as large as that in the minimum-mass disk model. If the dust particles are larger than chondrules, the secular gravitational instability grows within the lifetime of a protoplanetary disk.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Shock-induced star cluster formation in colliding galaxies

We studied the formation process of star clusters using high-resolution N-body/smoothed particle hydrodynamcs simulations of colliding galaxies. The total number of particles is 1.2x10^8 for our high resolution run. The gravitational softening is 5 pc and we allow gas to cool down to \sim 10 K. During the first encounter of the collision, a giant filament consists of cold and dense gas found between the progenitors by shock compression. A vigorous starburst took place in the filament, resulting in the formation of star clusters. The mass of these star clusters ranges from 10^{5-8} Msun. These star clusters formed hierarchically: at first small star clusters formed, and then they merged via gravity, resulting in larger star clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Formation of Terrestrial Planets from Protoplanets under a Realistic Accretion Condition

The final stage of terrestrial planet formation is known as the giant impact stage where protoplanets collide with one another to form planets. So far this stage has been mainly investigated by N-body simulations with an assumption of perfect accretion in which all collisions lead to accretion. However, this assumption breaks for collisions with high velocity and/or a large impact parameter. We derive an accretion condition for protoplanet collisions in terms of impact velocity and angle and masses of colliding bodies, from the results of numerical collision experiments. For the first time, we adopt this realistic accretion condition in N-body simulations of terrestrial planet formation from protoplanets. We compare the results with those with perfect accretion and show how the accretion condition affects terrestrial planet formation. We find that in the realistic accretion model, about half of collisions do not lead to accretion. However, the final number, mass, orbital elements, and even growth timescale of planets are barely affected by the accretion condition. For the standard protoplanetary disk model, typically two Earth-sized planets form in the terrestrial planet region over about 100 M years in both realistic and perfect accretion models. We also find that for the realistic accretion model, the spin angular velocity is about 30% smaller than that for the perfect accretion model that is as large as the critical spin angular velocity for rotational instability. The spin angular velocity and obliquity obey Gaussian and isotropic distributions, respectively, independently of the accretion condition.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gas Accretion onto a Protoplanet and Formation of a Gas Giant Planet

We investigate gas accretion onto a protoplanet, by considering the thermal effect of gas in three-dimensional hydrodynamical simulations, in which the wide region from a protoplanetary gas disk to a Jovian radius planet is resolved using the nested-grid method. We estimate the mass accretion rate and growth timescale of gas giant planets. The mass accretion rate increases with protoplanet mass for M_p < M_cri, while it becomes saturated or decreases for M_p > M_cri, where M_cri = 0.036 M_Jup (a_p/1AU)^0.75, and M_Jup and a_p are the Jovian mass and the orbital radius, respectively. The growth timescale of a gas giant planet or the timescale of the gas accretion onto the protoplanet is about 10^5 yr, that is two orders of magnitude shorter than the growth timescale of the solid core. The thermal effects barely affect the mass accretion rate because the gravitational energy dominates the thermal energy around the protoplanet. The mass accretion rate obtained in our local simulations agrees quantitatively well with those obtained in global simulations with coarser spatial resolution. The mass accretion rate is mainly determined by the protoplanet mass and the property of the protoplanetary disk. We find that the mass accretion rate is correctly calculated when the Hill or Bondi radius is sufficiently resolved. Using the oligarchic growth of protoplanets, we discuss the formation timescale of gas giant planets.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

N-Body Simulation of Planetesimal Formation through Gravitational Instability of a Dust Layer in Laminar Gas Disk

We investigate the formation process of planetesimals from the dust layer by the gravitational instability in the gas disk using local $N$-body simulations. The gas is modeled as a background laminar flow. We study the formation process of planetesimals and its dependence on the strength of the gas drag. Our simulation results show that the formation process is divided into three stages qualitatively: the formation of wake-like density structures, the creation of planetesimal seeds, and their collisional growth. The linear analysis of the dissipative gravitational instability shows that the dust layer is secularly unstable although Toomre's $Q$ value is larger than unity. However, in the initial stage, the growth time of the gravitational instability is longer than that of the dust sedimentation and the decrease in the velocity dispersion. Thus, the velocity dispersion decreases and the disk shrinks vertically. As the velocity dispersion becomes sufficiently small, the gravitational instability finally becomes dominant. Then wake-like density structures are formed by the gravitational instability. These structures fragment into planetesimal seeds. The seeds grow rapidly owing to mutual collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

N-Body Simulation of Planetesimal Formation through Gravitational Instability and Coagulation. II. Accretion Model

The gravitational instability of a dust layer is one of the scenarios for planetesimal formation. If the density of a dust layer becomes sufficiently high as a result of the sedimentation of dust grains toward the midplane of a protoplanetary disk, the layer becomes gravitationally unstable and spontaneously fragments into planetesimals. Using a shearing box method, we performed local $N$-body simulations of gravitational instability of a dust layer and subsequent coagulation without gas and investigated the basic formation process of planetesimals. In this paper, we adopted the accretion model as a collision model. A gravitationally bound pair of particles is replaced by a single particle with the total mass of the pair. This accretion model enables us to perform long-term and large-scale calculations. We confirmed that the formation process of planetesimals is the same as that in the previous paper with the rubble pile models. The formation process is divided into three stages: the formation of non-axisymmetric structures, the creation of planetesimal seeds, and their collisional growth. We investigated the dependence of the planetesimal mass on the simulation domain size. We found that the mean mass of planetesimals formed in simulations is proportional to $L_y^{3/2}$, where $L_y$ is the size of the computational domain in the direction of rotation. However, the mean mass of planetesimals is independent of $L_x$, where $L_x$ is the size of the computational domain in the radial direction if $L_x$ is sufficiently large. We presented the estimation formula of the planetesimal mass taking into account the simulation domain size.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Toward First-Principle Simulations of Galaxy Formation: II. Shock-Induced Starburst at a Collision Interface During the First Encounter of Interacting Galaxies

We investigated the evolution of interacting disk galaxies using high-resolution $N$-body/SPH simulations, taking into account the multiphase nature of the interstellar medium (ISM). In our high-resolution simulations, a large-scale starburst occurred naturally at the collision interface between two gas disks at the first encounter, resulting in the formation of star clusters. This is consistent with observations of interacting galaxies. The probability distribution function (PDF) of gas density showed clear change during the galaxy-galaxy encounter. The compression of gas at the collision interface between the gas disks first appears as an excess at $n_{\rm H} \sim 10{\rm cm^{-3}}$ in the PDF, and then the excess moves to higher densities ($n_{\rm H} \gtrsim 100{\rm cm^{-3}}$) in a few times $10^7$ years where starburst takes place. After the starburst, the PDF goes back to the quasi-steady state. These results give a simple picture of starburst phenomena in galaxy-galaxy encounters.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Toward first-principle simulations of galaxy formation: I. How should we choose star formation criteria in high-resolution simulations of disk galaxies?

We performed 3-dimensional N-body/SPH simulations to study how mass resolution and other model parameters such as the star formation efficiency parameter, C* and the threshold density, nth affect structures of the galactic gaseous/stellar disk in a static galactic potential. We employ 10^6 - 10^7 particles to resolve a cold and dense (T < 100 K & n_H > 100 cm^{-3}) phase. We found that structures of the ISM and the distribution of young stars are sensitive to the assumed nth. High-nth models with nth = 100 cm^{-3} yield clumpy multi-phase features in the ISM. Young stars are distributed in a thin disk of which half-mass scale height is 10 - 30 pc. In low-nth models with nth = 0.1 cm^{-3}, the stellar disk is found to be several times thicker, and the gas disk appears smoother than the high-nth models. A high-resolution simulation with high-nth is necessary to reproduce the complex structure of the gas disk. The global properties of the model galaxies in low-nth models, such as star formation histories, are similar to those in the high-nth models when we tune the value of C* so that they reproduce the observed relation between surface gas density and surface star formation rate density. We however emphasize that high-nth models automatically reproduce the relation, regardless of the values of C*. The ISM structure, phase distribution, and distributions of young star forming region are quite similar between two runs with values of C* which differ by a factor of 15. We also found that the timescale of the flow from n_H ~1 cm^{-3} to n_H > 100 cm^{-3} is about 5 times as long as the local dynamical time and is independent of the value of C*. The use of a high-nth criterion for star formation in high-resolution simulations makes numerical models fairy insensitive to the modelling of star formation. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The formation of Kuiper-belt Binaries through Exchange Reactions

Recent observations have revealed an unexpectedly high binary fraction among the Trans-Neptunian Objects (TNOs) that populate the Kuiper-belt. The discovered binaries have four characteristics they comprise a few percent of the TNOs, the mass ratio of their components is close to unity, their internal orbits are highly eccentric, and the orbits are more than 100 times wider than the primary's radius. In contrast, theories of binary asteroid formation tend to produce close, circular binaries. Therefore, a new approach is required to explain the unique characteristics of the TNO binaries. Two models have been proposed. Both, however, require extreme assumptions on the size distribution of TNOs. Here we show a mechanism which is guaranteed to produces binaries of the required type during the early TNO growth phase, based on only one plausible assumption, namely that initially TNOs were formed through gravitational instabilities of the protoplanetary dust layer.

0 citations0 reviewsNo code linkedOpen access