Search

The swing amplification is one of the mechanisms for spiral arm formation and gravitational scattering of stellar orbits by a high-density region in a disc galaxy is considered as its elementary process. During the swing amplification, the epicycle phases of stars are synchronized. We previously showed that gravitational scattering by a softened point mass can synchronize the epicycle phases. Here, we expand upon our previous work to consider gravitational scattering by a finite-sized wake, which we model by using a prolate body. We numerically simulate the stellar motion under the influence of the wake gravity, and we investigate the dependence of the orbital evolution on the wake properties. We find that phase synchronization is more effective for a wake with a pitch angle around π/2, and that the phase synchronization depends on the wake shape. Even for realistic pitch angles i ≲ π/4 the phase synchronization can be more effective than that by the softened point mass. Increasing the impact parameter of stars and decreasing the wake pitch angle are shown to increase the axis ratio at which synchronization occurs. Linear trailing high-density patterns form around the wake and enlarge the wake itself.

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., Accepted by Publications of the Astronomical Society of Japan, 9 pages, 7 figures, 2 tables

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., 16 pages, 14 figures. Accepted for publication in the Astrophysical Journal

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., 11 pages, 8 figures. Accepted for publication in ApJ. See also our previous paper arXiv:2205.13768

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}$)., 20 pages, 5 figures, Published in PASJ

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., 16 pages, 15 figures. Accepted for publication in ApJ

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., Comment: Accepted to ApJ

The Sun is thought to have been formed within a star cluster. The coexistence of 26Al-rich and 26Al-poor calcium–aluminum-rich inclusions indicates that a direct injection of 26Al-rich materials from a nearby core-collapse supernova would be expected to occur in the first 105 years of the existence of the Solar System. Therefore, at least one core-collapse supernova ought to 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 of experiencing at least one core-collapse supernova within the finite duration of star formation within the birth cluster. We find that the number of stars in the birth cluster may be significantly greater than that previously considered, depending on the duration of star formation.

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., 11 pages, 8 figures. Accepted for publication in A&A. See also arXiv:2210.00913 and arXiv:2205.13768

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., 59 pages, 32 figures, 12 tables, accepted for publication in PASJ, all RVs will be available online as supplementary after the publication

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., Comment: 8 pages, 6 figures, accepted for publication in ApJL

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., 31 pages, 16 figures, accepted for publication in AJ

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., 13 pages, 10 figures, accepted for publication in ApJ

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 $\alpha = 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 $\tau_{GI}$ at $r_{snow}$ is, $\tau_{GI} \simeq$ 70-80 Myr for $\alpha = 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 $\alpha = 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)., Comment: 17 pages, 7 figures, accepted by ApJ

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., Comment: 7 pages, 5 figures, 1 table, accepted for publication in A&A, a 5-minute video summary can be found here: https://exoplanet-talks.org/talk/401

In the standard formation scenario of planetary systems, planets form from a protoplanetary disk that consists of gas and dust. The scenario can be divided into three stages: (1) formation of planetesimals from dust, (2) formation of protoplanets from planetesimals, and (3) formation of planets from protoplanets. In stage (1), planetesimals form from dust through coagulation of dust grains and/or some instability of a dust layer. Planetesimals grow by mutual collisions to protoplanets or planetary embryos through runaway and oligarchic growth in stage (2). The final stage (3) of terrestrial planet formation is giant impacts among protoplanets while sweeping residual planetesimals. In the present paper, we review the elementary processes of terrestrial planet formation and discuss the extension of the standard scenario.

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., Comment: accepted for publication in ApJ

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., 15 pages, 4 figures, Accepted for publication in ApJ

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., Comment: 9 pages, 7 figures, accepted for publication in MNRAS

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., Comment: 16 pages, 10 figures, accepted for publication in ApJ

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., 11 pages, 8 figures. Accepted for publication in MNRAS

The Infrared Doppler (IRD) instrument is a fiber-fed high-resolution NIR spectrometer for the Subaru telescope covering the Y,J,H-bands simultaneously with a maximum spectral resolution of 70,000. The main purpose of IRD is a search for Earth-mass planets around nearby M-dwarfs by precise radial velocity measurements, as well as a spectroscopic characterization of exoplanet atmospheres. We report the current status of the instrument, which is undergoing commissioning at the Subaru Telescope, and the first light observation successfully done in August 2017. The general description of the instrument will be given including spectrometer optics, fiber injection system, cryogenic system, scrambler, and laser frequency comb. A large strategic survey mainly focused on late-type M-dwarfs is planned to start from 2019.

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

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., 5 pages, 4 figures, accepted for publication in A&A Letters

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 $\alpha\lesssim10^{-5}$. If the dust-to-gas ratio is increased twice, the gravitational instability occurs for $\alpha\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., Comment: 17 pages, 11 figures, accepted for publication in ApJ

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., 9 pages, 9 figures, accepted for publication in A&A

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., 13 pages, 5 figures, accepted for publication in ApJ Letters

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 $\alpha$ is less than about $4 \times 10^{-3}$. The upper limit of $\alpha$ 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., Comment: 38 pages, 14 figures, accepted for publication in ApJ

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)., 31 pages, 11 figures, accepted for publication in ApJ

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 $\gamma$ 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($\iota$ 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 $\iota$ Dra b, HD 5608 b, and HD 14067 b., Comment: 33 pases, 9 figures, accepted to ApJ

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., 23 pages, 10 figures, accepted for publication in ApJ

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., 11 pages, 9 figures, accepted for publication in ApJ