Source author record

Toru Tsuribe

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

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

9 published item(s)

preprint2016arXiv

Circumstellar Disks of the Most Vigorously Accreting Young Stars

Young stellar objects (YSOs) may not accumulate their mass steadily, as was previously thought, but in a series of violent events manifesting themselves as sharp stellar brightening. These events can be caused by fragmentation due to gravitational instabilities in massive gaseous disks surrounding young stars, followed by migration of dense gaseous clumps onto the star. We report our high angular resolution, coronagraphic near-infrared polarization imaging observations using the High Contrast Instrument for the Subaru Next Generation Adaptive Optics (HiCIAO) of the Subaru 8.2 m Telescope, towards four YSOs which are undergoing luminous accretion outbursts. The obtained infrared images have verified the presence of several hundred AUs scale arms and arcs surrounding these YSOs. In addition, our hydrodynamics simulations and radiative transfer models further demonstrate that these observed structures can indeed be explained by strong gravitational instabilities occurring at the beginning of the disk formation phase. The effect of those tempestuous episodes of disk evolution on star and planet formation remains to be understood.

preprint2016arXiv

Evolution of binary seeds in collapsing protostellar gas clouds

We perform three dimensional smoothed particle hydrodynamics (SPH) simulations of gas accretion onto the seeds of binary stars to investigate their short-term evolution. Our simulation setup is more realistic compared to the previous works by taking into account of dynamically evolving envelope with non-uniform distribution of gas density and angular momentum of accreting flow. Our initial condition includes a seed binary and a surrounding gas envelope, modelling the phase of core collapse of gas cloud when the fragmentation has already occurred. We assume that the seed binary has no eccentricity and no growth by gas accretion. The envelope is assumed to be an isothermal gas with no self-gravity. We run multiple simulations with different values of initial mass ratio $q_0$ (the ratio of secondary over primary mass) and gas temperature, and find a critical value of $q_{\rm c} = 0.25$ which distinguishes the later evolution of mass ratio $q$ as a function of time. If $q_0$ >~ $q_{\rm c}$, the secondary seed grows faster and $q$ increases monotonically towards unity. If $q_0$ <~ $q_{\rm c}$, on the other hand, the primary seed grows faster and $q$ is lower than $q_0$ at the end of the simulation. Based on our numerical results, we analytically calculate the long-term evolution of the seed binary including the growth of binary by gas accretion. We find that the seed binary with $q_0$ >~ $q_{\rm c}$ evolves towards an equal-mass binary star, and that with $q_0$ <~ $q_{\rm c}$ evolves to a binary with an extreme value of $q$. Binary separation is a monotonically increasing function of time for any $q_0$, suggesting that the binary growth by accretion does not lead to the formation of close binaries.

preprint2014arXiv

The Role for the Inner Disk in Mass Accretion to the Star in the Early Phase of Star Formation

A physical mechanism that drives FU Orionis-type outbursts is reconsidered. We study the effect of inner part of a circumstellar disk covering a region from near the central star to the radius of approximately $5$ AU (hereafter, the inner disk). Using the fluctuated mass accretion rate onto the inner disk $\dot{M}_{\rm out}$, we consider the viscous evolution of the inner disk and the time variability of the mass accretion rate onto the central star $\dot{M}_{\rm in}$ by means of numerical calculation of an unsteady viscous accretion disk in a one-dimensional axisymmetric model. First, we calculate the evolution of the inner disk assuming an oscillating $\dot{M}_{\rm out}$. It is shown that the time variability of $\dot{M}_{\rm in}$ does not coincide with $\dot{M}_{\rm out}$ due to viscous diffusion. Second, we investigate the properties of spontaneous outbursts with temporally constant $\dot{M}_{\rm out}$. Outburst occur only in a limited range of mass accretion rates onto the inner disk $10^{-10}<\dot{M}_{\rm out}< 3\times 10^{-6}~{\rm M}_{\odot} {\rm yr}^{-1}$ due to gravo-magneto limit cycle (GML). Finally, we discuss the case with a combination of episodic $\dot{M}_{\rm out}$ and accretion outbursts cause by the GML in the inner disk. The GML can drive accretion outbursts onto the star even for the case of fluctuating $\dot{M}_{\rm out}$, although fluctuations of $\dot{M}$ decay during transmitting the inner disk inwards. We newly identified two modes of outburst which are spontaneous one and stimulated one. In a stimulated mode of outburst, $\dot{M}_{\rm out}$ does appear directly in $\dot{M}_{\rm in}$ (the latter defining the stellar accretion luminosity). In a spontaneous mode of outburst, $\dot{M}_{\rm out}$ appears as the interval between outbursts.

preprint2013arXiv

Growth of a Protostar and a Young Circumstellar Disk with High Mass Accretion Rate onto the Disk

The growing process of both a young protostar and a circumstellar disk is investigated. Viscous evolution of a disk around a single star is considered with a model where a disk increases its mass by dynamically accreting envelope and simultaneously loses its mass via viscous accretion onto the central star. We focus on the circumstellar disk with high mass accretion rate onto the disk $\dot{M}=8.512c_{\rm s}^3/G$ as a result of dynamical collapse of rotating molecular cloud core. We study the origin of the surface density distribution and the origin of the disk-to-star mass ratio by means of numerical calculations of unsteady viscous accretion disk in one-dimensional axisymmetric model. It is shown that the radial profiles of the surface density $Σ$, azimuthal velocity $v_ϕ$, and mass accretion rate $\dot{M}$ in the inner region approach to the quasi-steady state. Profile of the surface density distribution in the quasi-steady state is determined as a result of angular momentum transport rather than its original distribution of angular momentum in the cloud core. It is also shown that the disk mass becomes larger than the central star in the long time limit as long as temporary constant mass flux onto the disk is assumed. After the mass infall rate onto the disk declines owing to the depletion of the parent cloud core, the disk-to-star mass ratio $M_{\rm disk}/M_*$ decreases. The disk-to-star mass ratio becomes smaller than unity after $t> 10^5 \rm yr$ and $t>10^6 \rm yr$ from the beginning of the accretion phase in the case with $α_0 =1 {\rm and} 0.1$, respectively, where $α_0 $ is the constant part of viscous parameter. In the case with $α_0 \leq 10^{-2}$, $M_{\rm disk}/M_*$ is still larger than unity at $2 \rm Myr$ from the beginning of the accretion phase.

preprint2013arXiv

Kelvin-Helmholtz Instabilities in Multi-Sized Dust Layers

We examine the effect of the dust size distribution on Kelvin-Helmholtz instabilities in the protoplanetary disk with dust sedimentation. With newly taking into account the dust size distribution, the growth rate of the Kelvin-Helmholtz instability is calculated using the linear stability analysis with the dust density distribution consistent with sedimentation. Dust abundance required for gravitational instabilities before the Kelvin-Helmholtz instability is derived from the linear stability analysis, and it is found that the required dust abundance significantly coincides with that estimated from the Richardson number. It is also found that when the dust size distribution is taken into account, the critical Richardson number for the onset of the Kelvin-Helmholtz instability tends to increase with dust abundance. This result is different from that in the case without the dust size distribution.

preprint2011arXiv

Gravitational Fragmentation of Expanding Shells. I. Linear Analysis

We perform a linear perturbation analysis of expanding shells driven by expansions of HII regions. The ambient gas is assumed to be uniform. As an unperturbed state, we develop a semi-analytic method for deriving the time evolution of the density profile across the thickness. It is found that the time evolution of the density profile can be divided into three evolutionary phases, deceleration-dominated, intermediate, and self-gravity-dominated phases. The density peak moves relatively from the shock front to the contact discontinuity as the shell expands. We perform a linear analysis taking into account the asymmetric density profile obtained by the semi-analytic method, and imposing the boundary conditions for the shock front and the contact discontinuity while the evolutionary effect of the shell is neglected. It is found that the growth rate is enhanced compared with the previous studies based on the thin-shell approximation. This is due to the boundary effect of the contact discontinuity and asymmetric density profile that were not taken into account in previous works.

preprint2011arXiv

Gravitational Fragmentation of Expanding Shells. II. Three-dimensional Simulations

We investigate the gravitational fragmentation of expanding shells driven by HII regions using the three-dimensional Lagrangian simulation codes based on the Riemann solver, called Godunov smoothed particle hydrodynamics. The ambient gas is assumed to be uniform. In order to attain high resolution to resolve the geometrically thin dense shell, we calculate not the whole but a part of the shell. We find that perturbations begin to grow earlier than the prediction of the linear analysis under the thin-shell approximation. The wavenumber of the most unstable mode is larger than that in the thin-shell linear analysis. The development of the gravitational instability is accompanied by the significant deformation of the contact discontinuity. These results are consistent with a linear analysis presented by Iwasaki et al. (2011) that have taken into account the density profile across the thickness and approximate shock and contact discontinuity boundary conditions. We derive useful analytic formulae for the fragment scale and the epoch when the gravitational instability begins to grow.

preprint2009arXiv

Fragmentation of a dynamically condensing radiative layer

In this paper, the stability of a dynamically condensing radiative gas layer is investigated by linear analysis. Our own time-dependent, self-similar solutions describing a dynamical condensing radiative gas layer are used as an unperturbed state. We consider perturbations that are both perpendicular and parallel to the direction of condensation. The transverse wave number of the perturbation is defined by $k$. For $k=0$, it is found that the condensing gas layer is unstable. However, the growth rate is too low to become nonlinear during dynamical condensation. For $k\ne0$, in general, perturbation equations for constant wave number cannot be reduced to an eigenvalue problem due to the unsteady unperturbed state. Therefore, direct numerical integration of the perturbation equations is performed. For comparison, an eigenvalue problem neglecting the time evolution of the unperturbed state is also solved and both results agree well. The gas layer is unstable for all wave numbers, and the growth rate depends a little on wave number. The behaviour of the perturbation is specified by $kL_\mathrm{cool}$ at the centre, where the cooling length, $L_\mathrm{cool}$, represents the length that a sound wave can travel during the cooling time. For $kL_\mathrm{cool}\gg1$, the perturbation grows isobarically. For $kL_\mathrm{cool}\ll1$, the perturbation grows because each part has a different collapse time without interaction. Since the growth rate is sufficiently high, it is not long before the perturbations become nonlinear during the dynamical condensation. Therefore, according to the linear analysis, the cooling layer is expected to split into fragments with various scales.

preprint2008arXiv

Gravitational Instability of Shocked Interstellar Gas Layers

In this paper we investigate gravitational instability of shocked gas layers using linear analysis. An unperturbed state is a self-gravitating isothermal layer which grows with time by the accretion of gas through shock fronts due to a cloud-cloud collision. Since the unperturbed state is not static, and cannot be described by a self-similar solution, we numerically solved the perturbation equations and directly integrated them over time. We took account of the distribution of physical quantities across the thickness. Linearized Rankine-Hugoniot relations were imposed at shock fronts as boundary conditions. The following results are found from our unsteady linear analysis: the perturbation initially evolves in oscillatory mode, and begins to grow at a certain epoch. The wavenumber of the fastest growing mode is given by $k=2\sqrt{2πGρ_\mathrm{E} {\cal M\mit}}/c_\mathrm{s}$, where $ρ_\mathrm{E}, c_\mathrm{s}$ and $\cal M\mit$ are the density of parent clouds, the sound velocity and the Mach number of the collision velocity, respectively. For this mode, the transition epoch from oscillatory to growing mode is given by $t_g = 1.2/\sqrt{2πGρ_\mathrm{E} {\cal M\mit}}$. The epoch at which the fastest growing mode becomes non-linear is given by $2.4δ_0^{-0.1}/\sqrt{2πG ρ_\mathrm{E}{\cal M\mit}}$, where $δ_0$ is the initial amplitude of the perturbation of the column density. As an application of our linear analysis, we investigate criteria for collision-induced fragmentation. Collision-induced fragmentation will occur only when parent clouds are cold, or $α_0=5c_\mathrm{s}^2 R/2G M < 1$, where $R$ and $M$ are the radius and the mass of parent clouds, respectively.