Source author record

Yusuke Tsukamoto

Yusuke Tsukamoto 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

17works
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

17 published item(s)

preprint2024arXiv

Formation of unipolar outflow and $\textit{protostellar rocket effect}$ in magnetized turbulent molecular cloud cores

Observed protostellar outflows exhibit a variety of asymmetrical features, including remarkable unipolar outflows and bending outflows. Revealing the formation and early evolution of such asymmetrical protostellar outflows, especially the unipolar outflows, is essential for a better understanding of the star and planet formation because they can dramatically change the mass accretion and angular momentum transport to the protostars and protoplanetary disks. Here, we perform the three-dimensional non-ideal magnetohydrodynamics simulations to investigate the formation and early evolution of the asymmetrical protostellar outflows in magnetized turbulent isolated molecular cloud cores. We find, for the first time to our knowledge, that the unipolar outflow forms even in the single low-mass protostellar system. The results show that the unipolar outflow is driven in the weakly magnetized cloud cores with the dimensionless mass-to-flux ratios of $μ=8$ and $16$. Furthermore, we find the $\textit{protostellar rocket effect}$ of the unipolar outflow, which is similar to the launch and propulsion of a rocket. The unipolar outflow ejects the protostellar system from the central dense region to the outer region of the parent cloud core, and the ram pressure caused by its ejection suppresses the driving of additional new outflows. In contrast, the bending bipolar outflow is driven in the moderately magnetized cloud core with $μ=4$. The ratio of the magnetic to turbulent energies of a parent cloud core may play a key role in the formation of asymmetrical protostellar outflows.

preprint2022arXiv

B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main

We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.

preprint2022arXiv

Calibration and performance of the readout system based on switched capacitor arrays for the Large-Sized Telescope of the Cherenkov Telescope Array

The Cherenkov Telescope Array (CTA) is the next-generation ground-based very-high-energy gamma-ray observatory. The Large-Sized Telescope (LST) of CTA is designed to detect gamma rays between 20 GeV and a few TeV with a 23-meter diameter mirror. We have developed the focal plane camera of the first LST, which has 1855 photomultiplier tubes (PMTs) and the readout system which samples a PMT waveform at GHz with switched capacitor arrays, Domino Ring Sampler ver4 (DRS4). To measure the precise pulse charge and arrival time of Cherenkov signals, we developed a method to calibrate the output voltage of DRS4 and the sampling time interval, as well as an analysis method to correct the spike noise of DRS4. Since the first LST was inaugurated in 2018, we have performed the commissioning tests and calibrated the camera. We characterised the camera in terms of the charge pedestal under various conditions of the night sky background, the charge resolution of each pixel, the charge uniformity of the whole camera, and the time resolutions with a test pulse and calibration laser.

preprint2022arXiv

Impact of dust size distribution including large dust grains on magnetic resistivity: an analytical approach

This paper investigates the impact of dust size distribution on magnetic resistivity. In particular, we focus on its impact when the maximum dust size significantly increases from sub-micron. The first half of the paper describes our calculation method for magnetic resistivity based on the model of \citet{1987ApJ...320..803D} and shows that the method reproduces the results of a more realistic chemical reaction network calculations reasonably well. Then, we describe the results of the resistivity calculations for dust distributions with large maximum dust grains. Our results show that resistivity tends to decrease with dust growth, which is particularly true when the dust size power exponent $q$ is $q=2.5$. On the other hand, the decrease is less pronounced when the dust size power exponent $q$ is $q=3.5$, i.e., when the small dust is also responsible for the dust cross-section. Our results suggest that detailed dust coagulation and fragmentation processes play a vital role in the magnetic resistivities in protostar formation.

preprint2020arXiv

Star-disk alignment in the protoplanetary disks: SPH simulation of the collapse of turbulent molecular cloud cores

We perform a series of three-dimensional smoothed particle hydrodynamics (SPH) simulations to study the evolution of the angle between the protostellar spin and the protoplanetary disk rotation axes (the star-disk angle $ψ_{\rm sd}$) in turbulent molecular cloud cores. While $ψ_{\rm sd}$ at the protostar formation epoch exhibits broad distribution up to $\sim 130^{\circ}$, $ψ_{\rm sd}$ decreases ($\lesssim 20^{\circ}$) in a timescale of $\sim 10^{4}$ yr. This timescale of the star-disk alignment, $t_{\rm alignment}$, corresponds basically to the mass doubling time of the central protostar, in which the protostar forgets its initial spin direction due to the mass accretion from the disk. Values of $ψ_{\rm sd}$ both at $t=10^2$ yr and $t=10^5$ yr after the protostar formation are independent of the ratios of thermal and turbulent energies to gravitational energy of the initial cloud cores: $α=E_{\rm thermal}/|E_{\rm gravity}|$ and $γ_{\rm turb}=E_{\rm turbulence}/|E_{\rm gravity}|$. We also find that a warped disk is possibly formed by the turbulent accretion flow from the circumstellar envelope.

preprint2020arXiv

The Effect of Misalignment between Rotation Axis and Magnetic Field on Circumstellar Disk

The formation of circumstellar disks is investigated using three-dimensional resistive magnetohydrodynamic simulations, in which the initial prestellar cloud has a misaligned rotation axis with respect to the magnetic field. We examine the effects of (i) the initial angle difference between the global magnetic field and the cloud rotation axis ($θ_0$) and (ii) the ratio of the thermal to gravitational energy ($α_0$). We study $16$ models in total and calculate the cloud evolution until $\sim \! 5000$ yr after protostar formation. Our simulation results indicate that an initial non-zero $θ_0$ ($> 0$) promotes the disk formation but tends to suppress the outflow driving, for models that are moderately gravitationally unstable, $α_0 \lesssim 1$. In these models, a large-sized rotationally-supported disk forms and a weak outflow appears, in contrast to a smaller disk and strong outflow in the aligned case ($θ_0 = 0$). Furthermore, we find that when the initial cloud is highly unstable with small $α_0$, the initial angle difference $θ_0$ does not significantly affect the disk formation and outflow driving.

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.

preprint2016arXiv

A revised condition for self-gravitational fragmentation of protoplanetary disks

Fragmentation of protoplanetary disks due to gravitational instabilities is a candidate of a formation mechanism of binary stars, brown dwarfs, and gaseous giant planets. The condition for the fragmentation has been thought that the disk cooling timescale is comparable to its dynamical timescale. However, some numerical simulations suggest that the fragmentation does not occur even if the cooling time is small enough, or the fragmentation can occur even when the cooling is inefficient. To reveal a realistic condition for fragmentation of self-gravitating disks, we perform two-dimensional numerical simulations that take into account the effect of the irradiation of the central star and radiation cooling of the disk, and precisely investigate the structure of the spiral arms formed in the protoplanetary disks. We show that the Toomre Q parameter in the spiral arms is an essential parameter for fragmentation. The spiral arms fragment only when Q < 0.6 in the spiral arms. We have further confirmed that this fragmentation condition observed in the numerical simulations can be obtained from the linear stability analysis for the self-gravitating spiral arms. These results indicate that the process of fragmentation of protoplanetary disks is divided into two stages: formation of the spiral arms in the disks; and fragmentation of the spiral arm. Our work reduces the condition for the fragmentation of the protoplanetary disks to the condition of the formation of the spiral arm that satisfies Q < 0.6.

preprint2016arXiv

Chemistry in a forming protoplanetary disk: main accretion phase

We investigate the chemistry in a radiation-hydrodynamics model of star-forming core which evolves from a cold ($\sim 10$ K) prestellar core to the main accretion phase in $\sim 10^5$ yr. A rotationally-supported gravitationally unstable disk is formed around a protostar. We extract the temporal variation of physical parameters in $\sim 1.5 \times 10^3$ SPH particles which end up in the disk, and perform post-processing calculations of the gas-grain chemistry adopting a three-phase model. Inside the disk, the SPH particles migrate both inward and outward. Since a significant fraction of volatiles such as CO can be trapped in the water-dominant ice in the three-phase model, the ice mantle composition depends not only on the current position in the disk but also on whether the dust grain has ever experienced higher temperatures than the water sublimation temperature. Stable molecules such as H$_2$O, CH$_4$, NH$_3$ and CH$_3$OH are already abundant at the onset of gravitational collapse and simply sublimated as the fluid parcels migrate inside the water snow line. On the other hand, various molecules such as carbon chains and complex organic molecules (COMs) are formed in the disk. COMs abundance sensitively depends on the outcomes of photodissociation and diffusion rates of photofragments in bulk ice mantle. As for S-bearing species, H$_2$S ice is abundant in the collapse phase. In the warm regions in the disk, H$_2$S is sublimated to be destroyed, while SO, H$_2$CS, OCS and SO$_2$ become abundant.

preprint2016arXiv

Magnetic field and early evolution of circumstellar disks

The magnetic field plays a central role in the formation and evolution of circumstellar disks. The magnetic field connects the rapidly rotating central region with the outer envelope and extracts angular momentum from the central region during gravitational collapse of the cloud core. This process is known as magnetic braking. Both analytical and multidimensional simulations have shown that disk formation is strongly suppressed by magnetic braking in moderately magnetized cloud cores in the ideal magnetohydrodynamic limit. On the other hand, recent observations have provided growing evidence of a relatively large disk several tens of astronomical units in size existing in some Class 0 young stellar objects. This introduces a serious discrepancy between the theoretical study and observations. Various physical mechanisms have been proposed to solve the problem of catastrophic magnetic braking, such as misalignment between the magnetic field and the rotation axis, turbulence, and non-ideal effect. In this paper, we review the mechanism of magnetic braking, its effect on disk formation and early evolution, and the mechanisms that resolve the magnetic braking problem. In particular, we emphasize the importance of non-ideal effects. The combination of magnetic diffusion and thermal evolution during gravitational collapse provides a robust formation process for the circumstellar disk at the very early phase of protostar formation. The rotation induced by the Hall effect can supply a sufficient amount of angular momentum for typical circumstellar disks around T Tauri stars. By examining the combination of the suggested mechanisms, we conclude that the circumstellar disks commonly form in the very early phase of protostar formation.

preprint2015arXiv

Effects of Ohmic and ambipolar diffusion on the formation and evolution of the first cores, protostars and circumstellar discs

We investigate the formation and evolution of a first core, protostar, and circumstellar disc with a three-dimensional non-ideal (including both Ohmic and ambipolar diffusion) radiation magnetohydrodynamics simulation. We found that the magnetic flux is largely removed by magnetic diffusion in the first core phase and that the plasma $β$ of the centre of the first core becomes large, $β>10^4$. Thus, proper treatment of first core phase is crucial in investigating the formation of protostar and disc. On the other hand, in an ideal simulation, $β\sim 10$ at the centre of the first core. The simulations with magnetic diffusion show that the circumstellar disc forms at almost the same time of protostar formation even with a relatively strong initial magnetic field (the value for the initial mass-to-flux ratio of the cloud core relative to the critical value is $μ=4$). The disc has a radius of $r \sim 1$ AU at the protostar formation epoch. We confirm that the disc is rotationally supported. We also show that the disc is massive ($Q\sim 1$) and that gravitational instability may play an important role in the subsequent disc evolution.

preprint2014arXiv

Effects of radiation transfer on the structure of self-gravitating disks, their fragmentation and evolution of the fragments

We investigate structure of self-gravitating disks, their fragmentation and evolution of the fragments (the clumps) using both analytic approach and three-dimensional radiation hydrodynamics simulations starting from molecular cores. The simulations show that non-local radiative transfer determines disk temperature. We find the disk structure is well described by an analytical model of quasi-steady self-gravitating disk with radial radiative transfer. Because the radiative process is not local and radiation from the interstellar medium cannot be ignored, the local radiative cooling would not be balanced with the viscous heating in a massive disk around a low mass star. In our simulations, there are cases in which the disk does not fragment even though it satisfies the fragmentation criterion based on disk cooling time ($Q \sim 1$ and $Ωt_{\rm cool}\sim 1$). This indicates that at least the criterion is not sufficient condition for fragmentation. We determine the parameter range for the host cloud core in which disk fragmentation occurs. In addition, we show that the temperature evolution of the center of the clump is close to that of typical first cores and the minimum initial mass of clumps to be about a few Jupiter mass.

preprint2013arXiv

An explicit scheme for ohmic dissipation with smoothed particle magneto-hydrodynamics

In this paper, we present an explicit scheme for Ohmic dissipation with smoothed particle magneto-hydrodynamics (SPMHD). We propose a SPH discretization of Ohmic dissipation and solve Ohmic dissipation part of induction equation with the supertime-stepping method (STS) which allows us to take a longer time-step than Courant-Friedrich-Levy stability condition. Our scheme is second-order accurate in space and first-order accurate in time. Our numerical experiments show that optimal choice of the parameters of STS for Ohmic dissipation of SPMHD is νsts ~ 0.01 and Nsts ~ 5.

preprint2013arXiv

Formation, orbital and thermal evolution, and survival of planetary-mass clumps in the early phase of circumstellar disk evolution

We report the results of our three-dimensional radiation hydrodynamics simulation of collapsing unmagnetized molecular cloud cores. We investigate the formation and evolution of the circumstellar disk and the clumps formed by disk fragmentation. Our simulation shows that disk fragmentation occurs in the early phase of circumstellar disk evolution and many clumps form. The clump can be represented by a polytrope sphere of index $n \sim 3$ and $n \gtrsim 4$ at central temperature $T_c \lesssim100$ K and $T_c \gtrsim 100$ K, respectively. We demonstrate, numerically and theoretically, that the maximum mass of the clump, beyond which it inevitably collapses, is $\sim 0.03 M_{\odot}$. The entropy of the clump increases during its evolution, implying that evolution is chiefly determined by mass accretion from the disk rather than by radiative cooling. Although most of the clumps rapidly migrate inward and finally fall onto the protostar, a few clumps remain in the disk. The central density and temperature of the surviving clump rapidly increase and the clump undergoes a second collapse within 1000 - 2000 years after its formation. In our simulation, three second cores of masses $0.2\msun$, $0.15\msun$, and $0.06\msun$ formed. These are protostars or brown dwarfs rather than protoplanets. For the clumps to survive as planetary-mass objects, the rapid mass accretion should be prevented by some mechanisms.

preprint2012arXiv

Formation and Early Evolution of Circumstellar Disks in Turbulent Molecular Cloud Cores

We investigate the formation and evolution of circumstellar disks in turbulent cloud cores until several 104 years after protostar formation using smoothed particle hydrodynamics (SPH) calculations. The formation and evolution process of circumstellar disk in turbulent cloud cores differs substantially from that in rigidly rotating cloud cores. In turbulent cloud cores, a filamentary structure appears before the protostar formation and the protostar forms in the filament. If the turbulence is initially sufficiently strong, the remaining filament twists around the protostar and directly becomes a rotation-supported disk. Upon formation, the disk orientation is generally misaligned with the angular momentum of its host cloud core and it dynamically varies during the main accretion phase, even though the turbulence is weak. This is because the angular momentum of the entire cloud core is mainly determined by the large scale velocity field whose wavelength is comparable to the cloud scale, whereas the angular momentum of the disk is determined by the local velocity field where the protostar forms and these two velocity fields do not correlate with each other. In the case of disk evolution in a binary or multiple stars, the disks are misaligned with each other at least during the main accretion phase, because there is no correlation between the velocity fields around the position where each protostar forms. In addition, each disk is also misaligned with the binary orbital plane. Such misalignment can explain the recent observations of misaligned disks and misaligned protostellar outflows.

preprint2011arXiv

Classification of the Circumstellar Disk Evolution During the Main Accretion Phase

We have carried out hydro-dynamical simulations to investigate the formation and evolution of protostar and circumstellar disks from the prestellar cloud. As the initial state, we adopt the molecular cloud core with two non-dimensional parameters representing the thermal and rotational energies. With these parameters, we make 17 models and calculate the cloud evolution ~ 10^4 years after the protostar formation. We find that early evolution of the star-disk system can be qualitatively classified into four modes: the massive disk, early fragmentation, late fragmentation, and protostar dominant modes. In the 'massive disk mode' to which the majority of models belong, the disk mass is greater than the protostellar mass for over 10^4 years and no fragmentation occurs in the circumstellar disk. The collapsing cloud shows fragmentation before the protostar formation in the 'early fragmentation mode'. The circumstellar disk shows fragmentation after the protostar formation in the 'late fragmentation mode', in which the secondary star substantially gains its mass from the circumstellar disk after fragmentation and it has a mass comparable to that of the primary star. The protostellar mass rapidly increases and exceeds the circumstellar disk mass in the 'protostar dominant mode'. This mode appears only when the initial molecular cloud core has a considerably small rotational energy. Comparison of our results with observations indicates that a majority of protostar has a fairly massive disk during the main accretion phase: the circumstellar disk mass is comparable to or more massive than the protostar. It is expected that such a massive disk promotes gas-giant formation by gravitational instability in a subsequent evolution stage.

preprint2011arXiv

The self gravity effect on the orbital stability of Twotinos

We investigate how the self gravity changes the orbital stability of Twotinos using N-body simulations in which the gravitational interaction between planetesimals are fully calculated. We show the timescale in which the half of Twotinos becomes unstable, t_{half}, obey the formula, t_{half} = 4 \times 106(Mtot0.1M\oplus)-1(mp7.6 \times 1023g)-1(hi2i1/20.002)(years), if we assume the primordial planetesimals disk have the power low surface mass density, Σ= Σ_0 \times r-2/3. Where Mtot, mp, hi2i1/2 are the total mass of large bodies of Twotinos, the maximum mass of planetesimals, and the inclination dispersion. With this formulae, we conclude the total mass of Twotinos is reduced to the order of 0.01 M\oplus by the self gravity and secular perturbation of the Planets even if there was huge mass such as several order of the earth mass in 1:2 MMR of Neptune at the early age of the solar system. These results will invoke reexamination to many previous works which tried to explain the dynamical evolution of TNOs.