Source author record

Masahiro N. Machida

Masahiro N. Machida 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

59works
7topics
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

59 published item(s)

preprint2022arXiv

Delivery of gas onto the circumplanetary disk of giant planets: Planetary-mass dependence of the source region of accreting gas and mass accretion rate

Gas accretion onto the circumplanetary disks and the source region of accreting gas are important to reveal dust accretion that leads to satellite formation around giant planets. We performed local three-dimensional high-resolution hydrodynamic simulations of isothermal and inviscid gas flow around a planet to investigate planetary-mass dependence of gas accretion band width and gas accretion rate onto circumplanetary disks. We examined cases with various planetary masses corresponding to M_p=0.05-1M_{Jup} at 5.2 au, where M_{Jup} is the current Jovian mass. We found that the radial width of the gas accretion band is proportional to M_p^{1/6} for the low-mass regime with M_p < 0.2 M_{Jup} while it is proportional to M_p for the high-mass regime with M_p > 0.2M_{Jup}. We found that the ratio of the mass accretion rate onto the circumplanetary disk to that into the Hill sphere is about 0.4 regardless of planetary mass for the cases we examined. Combining our results with the gap model obtained from global hydrodynamic simulations, we derive semi-analytical formulae of mass accretion rate onto circumplanetary disks. We found that the mass dependence of our three-dimensional accretion rates is the same as the previously-obtained two-dimensional case, although the qualitative behavior of accretion flow onto the CPD is quite different between the two cases.

preprint2022arXiv

Dust coagulation and fragmentation in a collapsing cloud core and their influence on non-ideal magnetohydrodynamic effects

We determine the time evolution of the dust particle size distribution during the collapse of a cloud core, accounting for both dust coagulation and dust fragmentation, to investigate the influence of dust growth on non-ideal magnetohydrodynamic effects.The density evolution of the collapsing core is given by a one-zone model. We assume two types of dust model: dust composed only of silicate (silicate dust) and dust with a surface covered by $\mathrm{H_{2}O}$ ice ($\mathrm{H_{2}O}$ ice dust). When only considering collisional coagulation, the non-ideal magnetohydrodynamic effects are not effective in the high-density region for both the silicate and $\mathrm{H_{2}O}$ ice dust cases. This is because dust coagulation reduces the abundance of small dust particles, resulting in less efficient adsorption of charged particles on the dust surface. For the silicate dust case, when collisional fragmentation is included, the non-ideal magnetohydrodynamic effects do apply at a high density of $n_{\mathrm{H}}>10^{12} \ \mathrm{cm^{-3}}$ because of the abundant production of small dust particles. On the other hand, for the $\mathrm{H_{2}O}$ ice dust case, the production of small dust particles due to fragmentation is not efficient. Therefore, for the $\mathrm{H_{2}O}$ ice dust case, non-ideal magnetohydrodynamic effects apply only in the range $n_{\mathrm{H}}\gtrsim 10^{14} \ \mathrm{cm^{-3}}$, even when collisional fragmentation is considered. Our results suggest that it is necessary to consider both dust collisional coagulation and fragmentation to activate non-ideal magnetohydrodynamic effects, which should play a significant role in the star and disk formation processes.

preprint2022arXiv

Environmental Effects of Star-Forming Cores on Mass Accretion Rate

We calculate the evolution of cloud cores embedded in different envelopes to investigate environmental effects on the mass accretion rate onto protostars. As the initial state, we neglect the magnetic field and cloud rotation, and adopt star-forming cores composed of two parts: a centrally condensed core and an outer envelope. The inner core has a critical Bonnor-Ebert density profile and is enclosed by the outer envelope. We prepare 15 star-forming cores with different outer envelope densities and gravitational radii, within which the gas flows into the collapsing core, and calculate their evolution until $\sim 2 \times10^5$ yr after protostar formation. The mass accretion rate decreases as the core is depleted when the outer envelope density is low. In contrast, the mass accretion rate is temporarily enhanced when the outer envelope density is high and the resultant protostellar mass exceeds the initial mass of the centrally condensed core. Some recent observations indicate that the mass of prestellar cores is too small to reproduce the stellar mass distribution. Our simulations show that the mass inflow from outside the core contributes greatly to protostellar mass growth when the core is embedded in a high-density envelope, which could explain the recent observations.

preprint2022arXiv

Exponentially amplified magnetic field eliminates disk fragmentation around the Population III protostar

One critical remaining issue to unclear the initial mass function of the first (Population III) stars is the final fate of secondary protostars formed in the accretion disk, specifically whether they merge or survive. We focus on the magnetic effects on the first star formation under the cosmological magnetic field. We perform a suite of ideal magnetohydrodynamic simulations until 1000 years after the first protostar formation. Instead of the sink particle technique, we employ a stiff equation of state approach to represent the magnetic field structure connecting to protostars. Ten years after the first protostar formation in the cloud initialized with $B_0 = 10^{-20}$ G at $n_0 = 10^4\,{\rm cm^{-3}}$, the magnetic field strength around protostars amplifies from pico- to kilo-gauss, which is the same strength as the present-day star. The magnetic field rapidly winds up since the gas in the vicinity of the protostar ($\leq\!10$ au) has undergone several tens orbital rotations in the first decade after protostar formation. As the mass accretion progresses, the vital magnetic field region extends outward, and the magnetic braking eliminates fragmentation of the disk that would form in the unmagnetized model. On the other hand, assuming a gas cloud with small angular momentum, this amplification might not work because the rotation would be slower. However, disk fragmentation would not occur in that case. We conclude that the exponential amplification of the cosmological magnetic field strength, about $10^{-18}$ G, eliminates disk fragmentation around the Population III protostars.

preprint2022arXiv

Implementation of dust particles in three-dimensional magnetohydrodynamics simulation: Dust dynamics in a collapsing cloud core

The aim of this study is to examine dust dynamics on a large scale and investigate the coupling of dust with gas fluid in the star formation process. We propose a method for calculating the dust trajectory in a gravitationally collapsing cloud, where the dust grains are treated as Lagrangian particles and are assumed to be neutral. We perform the dust trajectory calculations in combination with non-ideal magnetohydrodynamics simulation. Our simulation shows that dust particles with a size of $\le 10\,{\rm μm}$ are coupled with gas in a star-forming cloud core. We investigate the time evolution of the dust-to-gas mass ratio and the Stokes number, which is defined as the stopping time normalized by the freefall time-scale, and show that large dust grains ($\gtrsim 100\,{\rm μm}$) have a large Stokes number (close to unity) and tend to concentrate in the central region (i.e., protostar and rotationally supported disk) faster than do small grains ($\lesssim 10\,{\rm μm}$). Thus, large grains significantly increase the dust-to-gas mass ratio around and inside the disk. We also confirm that the dust trajectory calculations, which trace the physical quantities of each dust particle, reproduce previously reported results obtained using the Eulerian approach.

preprint2022arXiv

The First Detection of a Protostellar CO Outflow in the Small Magellanic Cloud with ALMA

Protostellar outflows are one of the most outstanding features of star formation. Observational studies over the last several decades have successfully demonstrated that outflows are ubiquitously associated with low- and high-mass protostars in the solar-metallicity Galactic conditions. However, the environmental dependence of protostellar outflow properties is still poorly understood, particularly in the low-metallicity regime. Here we report the first detection of a molecular outflow in the Small Magellanic Cloud with 0.2 $Z_{\odot}$, using Atacama Large Millimeter/submillimeter Array observations at a spatial resolution of 0.1 pc toward the massive protostar Y246. The bipolar outflow is nicely illustrated by high-velocity wings of CO(3-2) emission at $\gtrsim$15 km s$^{-1}$. The evaluated properties of the outflow (momentum, mechanical force, etc.) are consistent with those of the Galactic counterparts. Our results suggest that the molecular outflows, i.e., the guidepost of the disk accretion at the small scale, might be universally associated with protostars across the metallicity range of $\sim$0.2-1 $Z_{\odot}$.

preprint2022arXiv

Twisted magnetic field in star formation processes of L1521 F revealed by submillimeter dual band polarimetry using James Clerk Maxwell Telescope

Understanding the initial conditions of star formation requires both observational studies and theoretical works taking into account the magnetic field, which plays an important role in star formation processes. Herein, we study the young nearby dense cloud core L1521 F ($n$(H$_2$) $\sim 10^{4-6}$ cm$^{-3}$) in the Taurus Molecular Cloud. This dense core hosts a 0.2 $M_\odot$ protostar, categorized as a Very Low Luminosity Objects with complex velocity structures, particularly in the vicinity of the protostar. To trace the magnetic field within the dense core, we conducted high sensitivity submillimeter polarimetry of the dust continuum at $λ$= 850 $μ$m and 450 $μ$m using the POL-2 polarimeter situated in front of the SCUBA-2 submillimeter bolometer camera on James Clerk Maxwell Tetescope. This was compared with millimeter polarimetry taken at $λ$= 3.3 mm with ALMA. The magnetic field was detected at $λ$= 850 $μ$m in the peripheral region, which is threaded in a north-south direction, while the central region traced at $λ$= 450 $μ$m shows a magnetic field with an east-west direction, i.e., orthogonal to that of the peripheral region. Magnetic field strengths are estimated to be $\sim$70 $μ$G and 200 $μ$G in the peripheral- and central-regions, respectively, using the Davis-Chandrasekhar-Fermi method. The resulting mass-to-flux ratio of 3 times larger than that of magnetically critical state for both regions indicates that L1521 F is magnetically supercritical, i.e., gravitational forces dominate over magnetic turbulence forces. Combining observational data with MHD simulations, detailed parameters of the morphological properties of this puzzling object are derived for the first time.

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.

preprint2020arXiv

A low-velocity bipolar outflow from a deeply embedded object in Taurus revealed by the Atacama Compact Array

The first hydrostatic core, the first quasi-hydrostatic object formed during the star formation process, is still the observational missing link between the prestellar and protostellar phases, mainly due to its short lifetime. Although we have not established a clear method to identify this rare object, recent theoretical studies predict that the first core has millimeter continuum emission and low-velocity outflow with a wide opening angle. An extensive continuum/outflow survey toward a large number of $"$starless$"$ cores in nearby star-forming regions works as a pathfinder. We observed 32 prestellar cores in Taurus with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in 1.3 mm continuum and molecular lines using the Atacama Large Millimeter/submillimeter Array$-$Atacama Compact Array (ALMA$-$ACA) stand-alone mode. Among the targets, MC35-mm centered at one of the densest $"$starless$"$ cores in Taurus has blueshifted/redshifted wings in the $^{12}$CO (2-1) line, indicating that there is deeply embedded object driving molecular outflow. The observed velocities and sizes of the possible outflow lobes are 2-4 km s$^{-1}$, and $\sim$2 $\times$10$^3$ au, respectively, and the dynamical time is calculated to be $\sim$10$^3$ yr. In addition to this, the core is one of the strongest N$_2$D$^{+}$ (3-2) emitters in our sample. All of the observed signatures do not conflict with any of the theoretical predictions about the first hydrostatic core so far, and thus MC35-mm is unique as the only first-core candidate in the Taurus molecular cloud.

preprint2020arXiv

Different Modes of Star Formation II: Gas Accretion Phase of Initially Subcritical Star-Forming Clouds

The accretion phase of star formation is investigated in magnetically-dominated clouds that have an initial subcritical mass-to-flux ratio. We employ nonideal magnetohydrodynamic simulations that include ambipolar diffusion and ohmic dissipation. During the early prestellar phase the mass-to-flux ratio rises toward the critical value for collapse, and during this time the angular momentum of the cloud core is reduced significantly by magnetic braking. Once a protostar is formed in the core, the accretion phase is characterized by the presence of a small amount of angular momentum but a large amount of magnetic flux in the near-protostellar environment. The low angular momentum leads to a very small (or even nonexistent) disk and weak outflow, while the large magnetic flux can lead to an interchange instability that rapidly removes flux from the central region. The effective magnetic braking in the early collapse phase can even lead to a counter-rotating disk and outflow, in which the rotation direction of the disk and outflow is opposite to that of the infalling envelope. The solutions with a counter-rotating disk, tiny disk, or nonexistent disk (direct collapse) are unique outcomes that are realized in collapse from magnetically-dominated clouds with an initial subcritical mass-to-flux ratio.

preprint2020arXiv

Disk structure around the Class I protostar L1489 IRS revealed by ALMA: a warped disk system

We have observed the Class I protostar L1489 IRS with the Atacama Millimeter/submillimeter Array (ALMA) in Band 6. The C$^{18}$O $J=$2-1 line emission shows flattened and non-axisymmetric structures in the same direction as its velocity gradient due to rotation. We discovered that the C$^{18}$O emission shows dips at a radius of ~200-300 au while the 1.3 mm continuum emission extends smoothly up to r~400 au. At the radius of the C$^{18}$O dips, the rotational axis of the outer portion appears to be tilted by ~15 degrees from that of the inner component. Both the inner and outer components with respect to the C$^{18}$O dips exhibit the $r^{-0.5}$ Keplerian rotation profiles until r~600 au. These results not only indicate that a Keplerian disk extends up to ~600 au but also that the disk is warped. We constructed a three dimensional warped disk model rotating at the Keplerian velocity, and demonstrated that the warped disk model reproduces main observed features in the velocity channel maps and the PV diagrams. Such a warped disk system can form by mass accretion from a misaligned envelope. We also discuss a possible disk evolution scenario based on comparisons of disk radii and masses between Class I and Class II sources.

preprint2020arXiv

FRagmentation and Evolution of dense cores Judged by ALMA (FREJA). I (Overview). Inner $\sim$1000 au structures of prestellar/protostellar cores in Taurus

We have performed survey-type observations in 1 mm continuum and molecular lines toward dense cores (32 prestellar + 7 protostellar) with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in the Taurus molecular clouds using the Atacama Large Millimeter/submillimeter Array-Atacama Compact Array (ALMA-ACA) stand-alone mode with an angular resolution of 6.$''$5 ($\sim$900 au). The primary purpose of this study is to investigate the innermost part of dense cores toward understanding the initial condition of star formation. In the protostellar cores, contributions from protostellar disks dominate the observed continuum flux with a range of 35-90% except for the very low-luminosity object. For the prestellar cores, we have successfully confirmed continuum emission from dense gas with a density of $\gtrsim$3 $\times$10$^5$ cm$^{-3}$ toward approximately one-third of the targets. Thanks to the lower spatial frequency coverage with the ACA-7 m array, the detection rate is significantly higher than that of the previous surveys, which have 0 or 1 continuum detected sources among large number of starless samples using the ALMA Main array. The statistical counting method tells us that the lifetime of the prestellar cores until protostar formation therein approaches the free-fall time as the density increases. Among the prestellar cores, at least two targets have possible internal substructures, which are detected in continuum emission with the size scale of $\sim$1000 au if we consider the molecular line (C$^{18}$O and N$_2$D$^{+}$) distributions. These results suggest that small-scale fragmentation/coalescence processes occur in a region smaller than 0.1 pc, which may determine the final core mass associated with individual protostar formation before starting the dynamical collapse of the core with central density of $\sim$(0.3-1) $\times$ 10$^6$ cm$^{-3}$.

preprint2020arXiv

Magnetic Field Structure of Orion Source I

We observed polarization of the SiO rotational transitions from Orion Source I (SrcI) to probe the magnetic field in bipolar outflows from this high mass protostar. Both 43 GHz $J$=1-0 and 86 GHz $J$=2-1 lines were mapped with $\sim$20 AU resolution, using the Very Large Array (VLA) and Atacama Large Millimeter/Submillimeter Array (ALMA), respectively. The $^{28}$SiO transitions in the ground vibrational state are a mixture of thermal and maser emission. Comparison of the polarization position angles in the $J$=1-0 and $J$=2-1 transitions allows us to set an upper limit on possible Faraday rotation of $10^{4}$ radians m$^{-2}$, which would twist the $J$=2-1 position angles typically by less than 10 degrees. The smooth, systematic polarization structure in the outflow lobes suggests a well ordered magnetic field on scales of a few hundred AU. The uniformity of the polarization suggests a field strength of $\sim$30 milli-Gauss. It is strong enough to shape the bipolar outflow and possibly lead to sub-Keplerian rotation of gas at the base of the outflow. The strikingly high fractional linear polarizations of 80-90% in the $^{28}$SiO $v$=0 masers require anisotropic pumping. We measured circular polarizations of 60% toward the strongest maser feature in the $v$=0 $J$=1-0 peak. Anisotropic resonant scattering (ARS) is likely to be responsible for this circular polarization. We also present maps of the $^{29}$SiO $v$=0 $J$=2-1 maser and several other SiO transitions at higher vibrational levels and isotopologues.

preprint2020arXiv

Testing Disk Identification Methods Through Numerical Simulations of Protostellar Evolution

We test whether the radii of circumstellar disks can be reliably determined in observations applying the results of a numerical simulation. Firstly, we execute a core collapse simulation which starts from a rotating magnetized spherical core, and continue the calculation until the protostellar mass reaches 0.5 Msun. Then, for each set of simulation data, we calculate the radiative transfer to generate the data cube for the synthetic observation. The spatial and velocity resolutions of the synthetic observation are 0.15 arcsec (20 au) and 0.1 km/s, respectively. We define seven different disk radii. Four radii are estimated from the synthetic observation, using the continuum image, continuum visibility, C18O channel map, and C18O position velocity (PV) diagram. The other three radii are taken from the simulation and use the disk rotation, infall motion, and density contrast around the protostar to identify the disk. Finally, we compare the disk radii estimated from the systemic observation with those from the simulation. We find that the disk radius defined using the PV diagram can reliably trace the Keplerian disk when the protostellar mass is larger than M_*>~ 0.2 Msun independent of the inclination angle to the line of sight. In addition, the PV diagram provides an accurate estimate of the central stellar mass through the whole protostellar evolution. The simulation also indicates that the circumstellar disk is massive enough to be gravitationally unstable through the evolution. Such an unstable disk can show either a circular or spiral morphology on a similar timescale.

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.

preprint2020arXiv

Twin Jets and Close Binary Formation

The formation of a close binary system is investigated using a three-dimensional resistive magnetohydrodynamics simulation. Starting from a prestellar cloud, the cloud evolution is calculated until about 400 yr after protostar formation. Fragmentation occurs in the gravitationally collapsing cloud and two fragments evolve into protostars. The protostars orbit each other and a protobinary system appears. A wide-angle low-velocity outflow emerges from the circumbinary streams that encloses two protostars, while each protostar episodically drives high-velocity jets. Thus, the two high-velocity jets are surrounded by the low-velocity circumbinary outflow. The speed of the jets exceeds $\gtrsim 100 km s^{-1}$. Although the jets have a collimated structure, they are swung back on the small scale and are tangled at the large scale due to the binary orbital motion. A circumstellar disk also appears around each protostar. In the early main accretion phase, the binary orbit is complicated, while the binary separation is within $<30$ au. For the first time, all the characteristics of protobinary systems recently observed with large telescopes are reproduced in a numerical simulation.

preprint2019arXiv

Misalignment of Magnetic Fields, Outflows and Discs in Star-forming Clouds

Using resistive magnetohydrodynamics simulations, the propagation of protostellar jets, the formation of circumstellar discs and the configuration of magnetic fields are investigated from the prestellar cloud phase until $\sim$500 yr after protostar formation. As the initial state, we prepare magnetized rotating clouds, in which the rotation axis is misaligned with the global magnetic field by an angle $θ_0$. We calculate the cloud evolution for nine models with different $θ_0 (=$ 0, 5, 10, 30, 45, 60, 80, 85, 90$^\circ$). Our simulations show that there is no significant difference in the physical quantities of the protostellar jet, such as the mass and momentum, among the models except for the model with $θ_0=90^\circ$. On the other hand, the directions of the jet, disc normal and magnetic field are never aligned with each other during the early phase of star formation except for the model with $θ_0=0^\circ$. Even when the rotation axis of the prestellar cloud is slightly inclined to the global magnetic field, the directions of the jet, disc normal and local magnetic field differ considerably, and they randomly change over time. Our results indicate that it is very difficult to extract any information from the observations of the directions of the outflow, disc and magnetic field at the scale of $\lesssim 1000$ au. Thus, we cannot use such observations to derive any restrictions on the star formation process.

preprint2019arXiv

Protostellar Evolution in Serpens Main: Possible Origin of Disk-Size Diversity

We have observed the submillimeter continuum condensations SMM2, SMM4, SMM9, and SMM11 in the star forming cluster Serpens Main using the Atacama Large Millimeter/submillimeter Array during Cycle 3 in the 1.3 mm continuum, 12CO J=2-1, SO J_N=6_5-5_4, and C18O J=2-1 lines at an angular resolution of ~0.55 (240 au). Sixteen sources have been detected in the 1.3 mm continuum, which can be classified into three groups. Group 1 consists of six sources showing extended continuum emission and bipolar/monopolar 12CO outflows. Although all the Group 1 members are classified as Class 0 protostars, our observations suggest evolutionary trends among them in terms of 12CO outflow dynamical time, SO emission distribution, C18O fractional abundance, and continuum morphology. Group 2 consists of four sources associated with a continuum filamentary structure and no 12CO outflows. Central densities estimated from the 1.3 mm continuum intensity suggest that they are prestellar sources in a marginally Jeans unstable state. Group 3 consists of six Spitzer sources showing point-like 1.3 mm continuum emission and clumpy 12CO outflows. These features of Group 3 suggest envelope dissipation, preventing disk growth from the present size, r <~ 60 au. The Group 3 members are protostars that may be precursors to the T Tauri stars associated with small disks at tens-au radii identified in recent surveys.

preprint2016arXiv

ALMA band 8 continuum emission from Orion Source I

We have measured continuum flux densities of a high-mass protostar candidate, a radio source I in the Orion KL region (Orion Source I) using the Atacama Large Millimeter/Submillimeter Array (ALMA) at band 8 with an angular resolution of 0.1". The continuum emission at 430, 460, and 490 GHz associated with Source I shows an elongated structure along the northwest-southeast direction perpendicular to the so-called low-velocity bipolar outflow. The deconvolved size of the continuum source, 90 au times 20 au, is consistent with those reported previously at other millimeter/submillimeter wavelength. The flux density can be well fitted to the optically thick black-body spectral energy distribution (SED), and the brightness temperature is evaluated to be 700-800 K. It is much lower than that in the case of proton-electron or H- free-free radiations. Our data are consistent with the latest ALMA results by Plambeck & Wright (2016), in which the continuum emission have been proposed to arise from the edge-on circumstellar disk via thermal dust emission, unless the continuum source consists of an unresolved structure with the smaller beam filling factor.

preprint2016arXiv

Revealing a detailed mass distribution of a high-density core MC27/L1521F in Taurus with ALMA

We present the results of ALMA observations of dust continuum emission and molecular rotational lines toward a dense core MC27 (aka L1521F) in Taurus, which is considered to be at a very early stage of star formation. The detailed column density distribution on size scales from a few tens AU to ~10,000 AU scale are revealed by combining the ALMA (12 m array + 7 m array) data with the published/unpublished single-dish data. The high angular resolution observations at 0.87 mm with a synthesized beam size of ~0."74 x 0."32 reveal that a protostellar source, MMS-1, is not spatially resolved and lacks associated gas emission, while a starless high-density core, MMS-2, has substructures both in dust and molecular emission. The averaged radial column density distribution of the inner part of MC27/L1521F (r $\lesssim$ 3000 AU) is N(H2) ~r$^{-0.4}$, clearly flatter than that of the outer part, ~r$^{-1.0}$. The complex velocity/spatial structure obtained with previous ALMA observations is located inside the inner flatter region, which may reflect the dynamical status of the dense core.

preprint2016arXiv

The Origin of the Rotation Profiles in Star Forming Clouds

Angular momentum distribution and its redistribution are of crucial importance in formation and evolution of circumstellar disks. Many molecular line observations toward young stellar objects indicate that radial distributions of the specific angular momentum $j$ are more or less universal. In small scales, typically $R\lesssim$ 100 AU, the specific angular momenta distribute like $j\propto r^{1/2}$, indicating existence of rotationally supported disk. In outer regions, $R\gtrsim$ 5000 AU, $j$ increases as the radius increases and the slope is steeper than unity, which is supposed to reflect the original angular momentum distributions in the maternal molecular clouds. And lastly there is a connecting region, 100 AU $\lesssim R \lesssim$ 5000 AU, in which $j$-distribution looks almost flat. While this is often interpreted as a consequence of conservation of the specific angular momentum, it actually is insufficient and requires a stronger condition that the initial distribution of $j$ must be spatially uniform. However, this requirement is unrealistic and inconsistent with observations. In this work, we propose a simple alternative explanation; the flat $j$ profile is produced by strong prolongation due to the large velocity gradient in the accreting flow no matter what the initial $j$ distribution is. We provide a simple analytic model for gravitational collapse of molecular clouds. This model can be used to estimate ages of protostars based solely on the observed rotation profile. We demonstrate its validity in comparison with hydrodynamic simulations, and apply the model to young stellar objects such as L1527 IRS, TMC-1A and B335.

preprint2015arXiv

Accretion Phase of Star Formation in Clouds with Different Metallicities

The main accretion phase of star formation is investigated in clouds with different metallicities in the range of 0 \le Z \le Z_\odot, resolving the protostellar radius. Starting from a near-equilibrium prestellar cloud, we calculate the cloud evolution up to \sim100 yr after the first protostar formation. The star formation process considerably differs between clouds with lower (Z \le 10^-4 Z_\odot) and higher (Z > 10^-4 Z_\odot) metallicities. Fragmentation frequently occurs and many protostars appear without forming a stable circumstellar disc in lower-metallicity clouds. In these clouds, although protostars mutually interact and some are ejected from the cloud centre, many remain as a small stellar cluster. In contrast, higher-metallicity clouds produce a single protostar surrounded by a nearly stable rotation-supported disc. In these clouds, although fragmentation occasionally occurs in the disc, the fragments migrate inwards and finally fall onto the central protostar. The difference in cloud evolution is due to different thermal evolutions and mass accretion rates. The thermal evolution of the cloud determines the emergence and lifetime of the first core. The first core develops prior to the protostar formation in higher-metallicity clouds, whereas no (obvious) first core appears in lower-metallicity clouds. The first core evolves into a circumstellar disc with a spiral pattern, which effectively transfers the angular momentum outwards and suppresses frequent fragmentation. In lower-metallicity clouds, the higher mass accretion rate increases the disc surface density within a very short time, rendering the disc unstable to self-gravity and inducing vigorous fragmentation.

preprint2015arXiv

ALMA Observations of the Transition from Infall Motion to Keplerian Rotation around the Late-phase Protostar TMC-1A

We have observed the Class I protostar TMC-1A with Atacama Millimeter/submillimeter Array (ALMA) in 12CO and C18O (J=2-1), and 1.3-mm dust continuum emission. Continuum emission with a deconvolved size of 0.50"x0.37", perpendicular to the 12CO outflow, is detected. It most likely traces a circumstellar disk around TMC-1A, as previously reported. In contrast, a more extended structure is detected in C18O although it is still elongated with a deconvolved size of 3.3"x2.2", indicating that C18O traces mainly a flattened envelope surrounding the disk and the central protostar. C18O shows a clear velocity gradient perpendicular to the outflow at higher velocities, indicative of rotation, while an additional velocity gradient along the outflow is found at lower velocities. The radial profile of the rotational velocity is analyzed in detail, finding that it is given as a power-law \propto r^{-a} with an index of ~0.5 at higher velocities. This suggests that the rotation at higher velocities can be explained as Keplerian rotation orbiting a protostar with a dynamical mass of 0.68 Mo (inclination-corrected). The additional velocity gradient of C18O along the outflow is considered to be mainly infall motions in the envelope. Position-Velocity diagrams made from models consisting of an infalling envelope and a Keplerian disk are compared with the observations, revealing that the observed infall velocity is ~0.3 times smaller than the free fall velocity yielded by the dynamical mass of the protostar. Magnetic fields could be responsible for the slow infall velocity. A possible scenario of Keplerian disk formation is discussed.

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.

preprint2015arXiv

Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Non-Ideal Magnetohydrodynamic Effects and Early Formation of Circumstellar Disks

The transport of angular momentum by magnetic fields is a crucial physical process in formation and evolution of stars and disks. Because the ionization degree in star forming clouds is extremely low, non-ideal magnetohydrodynamic (MHD) effects such as ambipolar diffusion and Ohmic dissipation work strongly during protostellar collapse. These effects have significant impacts in the early phase of star formation as they redistribute magnetic flux and suppress angular momentum transport by magnetic fields. We perform three-dimensional nested-grid radiation magnetohydrodynamic (RMHD) simulations including Ohmic dissipation and ambipolar diffusion. Without these effects, magnetic fields transport angular momentum so efficiently that no rotationally supported disk is formed even after the second collapse. Ohmic dissipation works only in a relatively high density region within the first core and suppresses angular momentum transport, enabling formation of a very small rotationally supported disk after the second collapse. With both Ohmic dissipation and ambipolar diffusion, these effects work effectively in almost the entire region within the first core and significant magnetic flux loss occurs. As a result, a rotationally supported disk is formed even before a protostellar core forms. The size of the disk is still small, about 5 AU at the end of the first core phase, but this disk will grow later as gas accretion continues. Thus the non-ideal MHD effects can resolve the so-called magnetic braking catastrophe while maintaining the disk size small in the early phase, which is implied from recent interferometric observations.

preprint2014arXiv

Accretion of Solid Materials onto Circumplanetary Disks from Protoplanetary Disks

We investigate accretion of solid materials onto circumplanetary disks from heliocentric orbits rotating in protoplanetary disks, which is a key process for the formation of regular satellite systems. In the late stage of gas-capturing phase of giant planet formation, the accreting gas from protoplanetary disks forms circumplanetary disks. Since the accretion flow toward the circumplanetary disks affects the particle motion through gas drag force, we use hydrodynamic simulation data for the gas drag term to calculate the motion of solid materials. We consider wide range of size for the solid particles ($10^{-2}$-$10^6$m), and find that the accretion efficiency of the solid particles peaks around 10m-sized particles because energy dissipation of drag with circum-planetary disk gas in this size regime is most effective. The efficiency for particles larger than 10m size becomes lower because gas drag becomes less effective. For particles smaller than 10m, the efficiency is lower because the particles are strongly coupled with the back-ground gas flow, which prevent particles from accretion. We also find that the distance from the planet where the particles are captured by the circumplanetary disks is in a narrow range and well described as a function of the particle size.

preprint2014arXiv

ALMA observations of a high-density core in Taurus: dynamical gas interaction at the possible site of a multiple star formation

Starless dense cores eventually collapse dynamically, which forms protostars inside them, and the physical properties of the cores determine the nature of the forming protostars. We report ALMA observations of dust continuum emission and molecular rotational lines toward MC27 or L1521F, which is considered to be very close to the first protostellar core phase. We found a few starless high-density cores, one of which has a very high density of ~10^7 cm^-3, within a region of several hundred AU around a very low-luminosity protostar detected by Spitzer. A very compact bipolar outflow with a dynamical timescale of a few hundred years was found toward the protostar. The molecular line observation shows several cores with an arc-like structure, possibly due to the dynamical gas interaction. These complex structures revealed in the present observations suggest that the initial condition of star formation is highly dynamical in nature, which is considered to be a key factor in understanding fundamental issues of star formation such as the formation of multiple stars and the origin of the initial mass function of stars.

preprint2014arXiv

ALMA Observations of Infalling Flows toward the Keplerian Disk around the Class I Protostar L1489 IRS

We have conducted ALMA observations in the 1.3 mm continuum and 12CO (2-1), C18O (2-1) and SO (5_6-4_5) lines toward L1489 IRS, a Class I protostar surrounded by a Keplerian disk and an infalling envelope. The Keplerian disk is clearly identified in the 12CO and C18O emission, and its outer radius (~700 AU) and mass (~0.005 Msun) are comparable to those of disks around T Tauri stars. The protostellar mass is estimated to be 1.6 Msun with the inclination angle of 66 deg. In addition to the Keplerian disk, there are blueshifted and redshifted off-axis protrusions seen in the C18O emission pointing toward the north and the south, respectively, adjunct to the middle part of the Keplerian disk. The shape and kinematics of these protrusions can be interpreted as streams of infalling flows with a conserved angular momentum following parabolic trajectories toward the Keplerian disk, and the mass infalling rate is estimated to be ~5E-7 Msun/yr. The specific angular momentum of the infalling flows (~2.5E-3 km/s*pc) is comparable to that at the outer radius of the Keplerian disk (~4.8E-3 km/s*pc). The SO emission is elongated along the disk major axis and exhibits a linear velocity gradient along the axis, which is interpreted as that the SO emission primarily traces a ring region in the flared Keplerian disk at radii of ~250-390 AU. The local enhancement of the SO abundance in the ring region can be due to the accretion shocks at the centrifugal radius where the infalling flows fall onto the disk. Our ALMA observations unveiled both the Keplerian disk and the infalling gas onto the disk, and the disk can further grow by accreting material and angular momenta from the infalling gas.

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.

preprint2014arXiv

Formation of a Keplerian disk in the infalling envelope around L1527 IRS: transformation from infalling motions to Kepler motions

We report Atacama Large Millimeter/submillimeter Array (ALMA) cycle 0 observations of C$^{18}$O ($J=2-1$), SO ($J_N= 6_5-5_4$) and 1.3mm dust continuum toward L1527 IRS, a class 0 solar-type protostar surrounded by an infalling and rotating envelope. C$^{18}$O emission shows strong redshifted absorption against the bright continuum emission associated with L1527 IRS, strongly suggesting infall motions in the C$^{18}$O envelope. The C$^{18}$O envelope also rotates with a velocity mostly proportional to $r^{-1}$, where $r$ is the radius, while the rotation profile at the innermost radius (54 AU) may be shallower than $r^{-1}$, suggestive of formation of a Keplerian disk around the central protostar of 0.3 Mo in dynamical mass. SO emission arising from the inner part of the C$^{18}$O envelope also shows rotation in the same direction as the C$^{18}$O envelope. The rotation is, however, rigid-body like which is very different from the differential rotation shown by C$^{18}$O. In order to explain the line profiles and the position-velocity (PV) diagrams of C$^{18}$O and SO observed, simple models composed of an infalling envelope surrounding a Keplerian disk of 54 AU in radius orbiting a star of 0.3 Mo are examined. It is found that in order to reproduce characteristic features of the observed line profiles and PV diagrams, the infall velocity in the model has to be smaller than the free-fall velocity yielded by a star of 0.3 Mo. Possible reasons for the reduced infall velocities are discussed.

preprint2014arXiv

Protostellar Jets Enclosed by Low-velocity Outflows

A protostellar jet and outflow are calculated for \sim 270 yr following the protostar formation using a three dimensional magnetohydrodynamics simulation, in which both the protostar and its parent cloud are spatially resolved. A high-velocity (\sim100km/s) jet with good collimation is driven near the disk's inner edge, while a low-velocity (<10km/s) outflow with a wide opening angle appears in the outer-disk region. The high-velocity jet propagates into the low-velocity outflow, forming a nested velocity structure in which a narrow high-velocity flow is enclosed by a wide low-velocity flow. The low-velocity outflow is in a nearly steady state, while the high-velocity jet appears intermittently. The time-variability of the jet is related to the episodic accretion from the disk onto the protostar, which is caused by gravitational instability and magnetic effects such as magnetic braking and magnetorotational instability. Although the high-velocity jet has a large kinetic energy, the mass and momentum of the jet are much smaller than those of the low-velocity outflow. A large fraction of the infalling gas is ejected by the low-velocity outflow. Thus, the low-velocity outflow actually has a more significant effect than the high-velocity jet in the very early phase of the star formation.

preprint2013arXiv

A Semi-Analytical Description for the Formation and Gravitational Evolution of Protoplanetary Disks

We investigate the formation process of self-gravitating protoplanetary disks in unmagnetized molecular clouds. The angular momentum is redistributed by the action of gravitational torques in the massive disk during its early formation. We develop a simplified one-dimensional accretion disk model that takes into account the infall of gas from the envelope onto the disk and the transfer of angular momentum in the disk with an effective viscosity. First we evaluate the gas accretion rate from the cloud core onto the disk by approximately estimating the effects of gas pressure and gravity acting on the cloud core. We formulate the effective viscosity as a function of the Toomre Q parameter that measures the local gravitational stability of the rotating thin disk. We use a function for viscosity that changes sensitively with Q when the disk is gravitationally unstable. We find a strong self-regulation mechanism in the disk evolution. During the formation stage of protoplanetary disks, the evolution of the surface density does not depend on the other details of the modeling of effective viscosity, such as the prefactor of the viscosity coefficient. Next, to verify our model, we compare the time evolution of the disk calculated with our formulation with that of three-dimensional hydrodynamical simulations. The structures of the resultant disks from the one-dimensional accretion disk model agree well with those of the three-dimensional simulations. Our model is a useful tool for the further modeling of chemistry, radiative transfer, and planet formation in protoplanetary disks.

preprint2013arXiv

Conditions for Circumstellar Disk Formation: Effects of Initial Conditions and Sink Treatment

The formation of a circumstellar disk in collapsing cloud cores is investigated with three-dimensional magnetohydrodynamic simulations. We prepare four types of initial cloud having different density profiles and calculate their evolution with or without a sink. To investigate the effect of magnetic dissipation on disk formation, the Ohmic dissipation is considered in some models. Calculations show that disk formation is very sensitive to both the initial cloud configuration and the sink treatment. The disk size considerably differs in clouds with different density profiles even when the initial clouds have almost the same mass-to-flux ratio. Only a very small disk (\sim 10 AU in size) appears in clouds with a uniform density profile, whereas a large disk (\sim 100 AU in size) forms in clouds with a Bonnor-Ebert density profile. In addition, a large sink accretion radius numerically impedes disk formation during the main accretion phase and tends to foster the misleading notion that disk formation is completely suppressed by magnetic braking. The protostellar outflow is also greatly affected by the sink properties. A sink accretion radius of \lesssim 1 AU and sink threshold density of \gtrsim 10^ 13 cm^-3 are necessary for investigating disk formation during the main accretion phase.

preprint2013arXiv

Evolution of Protostellar Outflow around Low-mass Protostar

The evolution of protostellar outflow is investigated with resistive magneto-hydrodynamic nested-grid simulations that cover a wide range of spatial scales (\sim 1AU - 1pc). We follow cloud evolution from the pre-stellar core stage until the infalling envelope dissipates long after the protostar formation. We also calculate protostellar evolution to derive protostellar luminosity with time-dependent mass accretion through a circumstellar disk. The protostellar outflow is driven by the first core prior to protostar formation and is directly driven by the circumstellar disk after protostar formation. The opening angle of the outflow is large in the Class 0 stage. A large fraction of the cloud mass is ejected in this stage, which reduces the star formation efficiency to \sim 50%. After the outflow breaks out from the natal cloud, the outflow collimation is gradually improved in the Class I stage. The head of the outflow travels more than \sim 10^5AU in \sim 10^5 yr. The outflow momentum, energy and mass derived in our calculations agree well with observations. In addition, our simulations show the same correlations among outflow momentum flux, protostellar luminosity and envelope mass as those in observations. These correlations differ between Class 0 and I stages, which is explained by different evolutionary stages of the outflow; in the Class 0 stage, the outflow is powered by the accreting mass and acquires its momentum from the infalling envelope; in the Class I stage, the outflow enters the momentum-driven snow-plough phase. Our results suggest that protostellar outflow should determine the final stellar mass and significantly affect the early evolution of low-mass protostars.

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.

preprint2013arXiv

The Formation of Population III Stars in Gas Accretion Stage: Effects of Magnetic Fields

The formation of Population III stars is investigated using resistive magnetohydrodynamic simulations. Starting from a magnetized primordial prestellar cloud, we calculate the cloud evolution several hundreds of years after first protostar formation, resolving the protostellar radius. When the natal minihalo field strength is weaker than B \lesssim 10^-13 (n/1 cm^-3)^-2/3 G (n is the hydrogen number density), magnetic effects can be ignored. In this case, fragmentation occurs frequently and a stellar cluster forms, in which stellar mergers and mass exchange between protostars contribute to the mass growth of these protostars. During the early gas accretion phase, the most massive protostar remains near the cloud centre, whereas some of the less massive protostars are ejected. The magnetic field significantly affects Population III star formation when B_amb \gtrsim 10^-12 (n/1 cm^-3)^-2/3 G. In this case, because the angular momentum around the protostar is effectively transferred by both magnetic braking and protostellar jets, the gas falls directly onto the protostar without forming a disk, and only a single massive star forms. In addition, a massive binary stellar system appears when B_amb \sim 10^-12 (n/1cm^-3)^-2/3 G. Therefore, the magnetic field determines the end result of the formation process (cluster, binary or single star) for Population III stars. Moreover, no persistent circumstellar disk appears around the protostar regardless of the magnetic field strength, which may influence the further evolution of Population III stars.

preprint2012arXiv

Exploring Magnetic Field Structure in Star-Forming Cores with Polarization of Thermal Dust Emission

The configuration and evolution of the magnetic field in star-forming cores are investigated in order to directly compare simulations and observations. We prepare four different initial clouds having different magnetic field strengths and rotation rates, in which magnetic field lines are aligned/misaligned with the rotation axis. First, we calculate the evolution of such clouds from the prestellar stage until long after protostar formation. Then, we calculate the polarization of thermal dust emission expected from the simulation data. We create polarization maps with arbitrary viewing angles and compare them with observations. Using this procedure, we confirmed that the polarization distribution projected on the celestial plane strongly depends on the viewing angle of the cloud. Thus, by comparing the observations with the polarization map predicted by the simulations, we can roughly determine the angle between the direction of the global magnetic field and the line of sight. The configuration of the polarization vectors also depends on the viewing angle. We find that an hourglass configuration of magnetic field lines is not always realized in a collapsing cloud when the global magnetic field is misaligned with the cloud rotation axis. Depending on the viewing angle, an S-shaped configuration of the magnetic field (or the polarization vectors) appears early in the protostellar accretion phase. This indicates that not only the magnetic field but also the cloud rotation affects the dynamical evolution of such a cloud. In addition, by comparing the simulated polarization with actual observations, we can estimate properties of the host cloud such as the evolutionary stage, magnetic field strength, and rotation rate.

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.

preprint2012arXiv

Magnetic Field in The Isolated Massive Dense Clump IRAS 20126+4104

We measured polarized dust emission at 350um towards the high-mass star forming massive dense clump IRAS 20126+4104 using the SHARC II Polarimeter, SHARP, at the Caltech Submillimeter Observatory. Most of the observed magnetic field vectors agree well with magnetic field vectors obtained from a numerical simulation for the case when the global magnetic field lines are inclined with respect to the rotation axis of the dense clump. The results of the numerical simulation show that rotation plays an important role on the evolution of the massive dense clump and its magnetic field. The direction of the cold CO 1-0 bipolar outflow is parallel to the observed magnetic field within the dense clump as well as the global magnetic field, as inferred from optical polarimetry data, indicating that the magnetic field also plays a critical role in an early stage of massive star formation. The large-scale Keplerian disk of the massive (proto)star rotates in almost opposite sense to the clump's envelope. The observed magnetic field morphology and the counter-rotating feature of the massive dense clump system provide hints to constrain the role of magnetic fields in the process of high mass star formation.

preprint2012arXiv

Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Protostellar Core Formation

We report the first three-dimensional radiation magnetohydrodynamic (RMHD) simulations of protostellar collapse with and without Ohmic dissipation.We take into account many physical processes required to study star formation processes, including a realistic equation of state. We follow the evolution from molecular cloud cores until protostellar cores are formed with sufficiently high resolutions without introducing a sink particle. The physical processes involved in the simulations and adopted numerical methods are described in detail.We can calculate only about one year after the formation of the protostellar cores with our direct three-dimensional RMHD simulations because of the extremely short timescale in the deep interior of the formed protostellar cores, but successfully describe the early phase of star formation processes. The thermal evolution and the structure of the first and second (protostellar) cores are consistent with previous one-dimensional simulations using full radiation transfer, but differ considerably from preceding multi-dimensional studies with the barotropic approximation. The protostellar cores evolve virtually spherically symmetric in the ideal MHD models because of efficient angular momentum transport by magnetic fields, but Ohmic dissipation enables the formation of the circumstellar disks in the vicinity of the protostellar cores as in previous MHD studies with the barotropic approximation. The formed disks are still small (less than 0.35 AU) because we simulate only the earliest evolution. We also confirm that two different types of outflows are naturally launched by magnetic fields from the first cores and protostellar cores in the resistive MHD models.

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

Distribution of Accreting Gas and Angular Momentum onto Circumplanetary Disks

We investigate gas accretion flow onto a circumplanetary disk from a protoplanetary disk in detail by using high-resolution three-dimensional nested-grid hydrodynamic simulations, in order to provide a basis of formation processes of satellites around giant planets. Based on detailed analyses of gas accretion flow, we find that most of gas accretion onto circumplanetary disks occurs nearly vertically toward the disk surface from high altitude, which generates a shock surface at several scale heights of the circumplanetary disk. The gas that has passed through the shock surface moves inward because its specific angular momentum is smaller than that of the local Keplerian rotation, while gas near the midplane in the protoplanetary disk cannot accrete to the circumplanetary disk. Gas near the midplane within the planet's Hill sphere spirals outward and escapes from the Hill sphere through the two Lagrangian points L$_1$ and L$_2$. We also analyze fluxes of accreting mass and angular momentum in detail and find that the distributions of the fluxes onto the disk surface are well described by power-law functions and that a large fraction of gas accretion occurs at the outer region of the disk, i.e., at about 0.1 times the Hill radius. The nature of power-law functions indicates that, other than the outer edge, there is no specific radius where gas accretion is concentrated. These source functions of mass and angular momentum in the circumplanetary disk would provide us with useful constraints on the structure and evolution of the circumplanetary disk, which is important for satellite formation.

preprint2011arXiv

Effect of Magnetic Braking on the Circumstellar Disk Formation in a Strongly Magnetized Cloud

Using resistive magnetohydrodynamics simulation, we investigate circumstellar disk formation in a strongly magnetized cloud. As the initial state, an isolated cloud core embedded in a low-density interstellar medium with a uniform magnetic field is adopted. The cloud evolution is calculated until almost all gas inside the initial cloud falls onto either the circumstellar disk or a protostar, and a part of the gas is ejected into the interstellar medium by the protostellar outflow driven by the circumstellar disk. In the early main accretion phase, the disk size is limited to \sim 10 AU because the angular momentum of the circumstellar disk is effectively transferred by both magnetic braking and the protostellar outflow. In the later main accretion phase, however, the circumstellar disk grows rapidly and exceeds 100 AU by the end of the main accretion phase. This rapid growth of the circumstellar disk is caused by the depletion of the infalling envelope, while magnetic braking is effective when the infalling envelope is more massive than the circumstellar disk. The infalling envelope cannot brake the circumstellar disk when the latter is more massive than the former. In addition, the protostellar outflow weakens and disappears in the later main accretion phase, because the outflow is powered by gas accretion onto the circumstellar disk. Although the circumstellar disk formed in a magnetized cloud is considerably smaller than that in an unmagnetized cloud, a circumstellar disk exceeding 100 AU can form even in a strongly magnetized cloud.

preprint2011arXiv

Exposed Long-lifetime First-core: A New Model of First Cores Based on Radiation Hydrodynamics

A first adiabatic core is a transient object formed in the early phase of star formation. The observation of a first core is believed to be difficult because of its short lifetime and low luminosity. On the basis of radiation hydrodynamic simulations, we propose a novel theoretical model of first cores, Exposed Long-lifetime First core (ELF). In the very low-mass molecular core, the first core evolves slowly and lives longer than 10,000 years because the accretion rate is considerably low. The evolution of ELFs is different from that of ordinary first cores because radiation cooling has a significant effect there. We also carry out radiation transfer calculation of dust-continuum emission from ELFs to predict their observational properties. ELFs have slightly fainter but similar SEDs to ordinary first cores in radio wavelengths, therefore they can be observed. Although the probabilities that such low mass cores become gravitationally unstable and start to collapse are low, we still can expect that a considerable number of ELFs can be formed because there are many low-mass molecular cloud cores in star-forming regions that can be progenitors of ELFs.

preprint2011arXiv

Impact of Protostellar Outflow on Star Formation: Effects of Initial Cloud Mass

Star formation efficiency controlled by the protostellar outflow in a single cloud core is investigated by three-dimensional resistive MHD simulations. Starting from the prestellar cloud core, the star formation process is calculated until the end of the main accretion phase. In the calculations, the mass of the prestellar cloud is parameterized. During the star formation, the protostellar outflow is driven by the circumstellar disk. The outflow extends also in the transverse direction until its width becomes comparable to the initial cloud scale, and thus, the outflow has a wide opening angle of >40 degrees. As a result, the protostellar outflow sweeps up a large fraction of the infalling material and ejects it into the interstellar space. The outflow can eject at most over half of the host cloud mass, significantly decreasing star formation efficiency. The outflow power is stronger in clouds with a greater initial mass. Thus, the protostellar outflow effectively suppresses star formation efficiency in a massive cloud. The outflow weakens significantly and disappears in several free-fall timescales of the initial cloud after the cloud begins to collapse. The natal prestellar core influences the lifetime and size of the outflow. At the end of the main accretion phase, a massive circumstellar disk comparable in mass to the protostar remains. Calculations show that typically, ~30% of the initial cloud mass is converted into the protostar and ~20% remains in the circumstellar disk, while ~40% is ejected into the interstellar space by the protostellar outflow. Therefore, a single cloud core typically has a star formation efficiency of 30-50%.

preprint2011arXiv

Recurrent Planet Formation and Intermittent Protostellar Outflows Induced by Episodic Mass Accretion

The formation and evolution of a circumstellar disk in magnetized cloud cores is investigated from prestellar core stage until sim 10^4 yr after protostar formation. In the circumstellar disk, fragmentation first occurs due to gravitational instability in a magnetically inactive region, and substellar-mass objects appear. The substellar-mass objects lose their orbital angular momenta by gravitational interaction with the massive circumstellar disk and finally fall onto the protostar. After this fall, the circumstellar disk increases its mass by mass accretion and again induces fragmentation. The formation and falling of substellar-mass objects are repeated in the circumstellar disk until the end of the main accretion phase. In this process, the mass of fragments remain small, because the circumstellar disk loses its mass by fragmentation and subsequent falling of fragments before it becomes very massive. In addition, when fragments orbit near the protostar, they disturb the inner disk region and promote mass accretion onto the protostar. The orbital motion of substellar-mass objects clearly synchronizes with the time variation of the accretion luminosity of the protostar. Moreover, as the objects fall, the protostar shows a strong brightening for a short duration. The intermittent protostellar outflows are also driven by the circumstellar disk whose magnetic field lines are highly tangled owing to the orbital motion of fragments. The time-variable protostellar luminosity and intermittent outflows may be a clue for detecting planetary-mass objects in the circumstellar disk.

preprint2010arXiv

Emergence of Protoplanetary Disks and Successive Formation of Gaseous Planets by Gravitational Instability

We use resistive magnetohydrodynamical simulations with the nested grid technique to study the formation of protoplanetary disks around protostars from molecular cloud cores that provide the realistic environments for planet formation. We find that gaseous planetary-mass objects are formed in the early evolutionary phase by gravitational instability in regions that are decoupled from the magnetic field and surrounded by the injection points of the magnetohydrodynamical outflows during the formation phase of protoplanetary disks. Magnetic decoupling enables massive disks to form and these are subject to gravitational instability, even at ~ 10 AU. The frequent formation of planetary-mass objects in the disk suggests the possibility of constructing a hybrid planet formation scenario, where the rocky planets form later under the influence of the giant planets in the protoplanetary disk.

preprint2010arXiv

Formation Process of the Circumstellar Disk: Long-term Simulations in the Main Accretion Phase of Star Formation

The formation and evolution of the circumstellar disk in unmagnetized molecular clouds is investigated using three-dimensional hydrodynamic simulations from the prestellar core until the end of the main accretion phase. In collapsing clouds, the first (adiabatic) core with a size of ~10AU forms prior to the formation of the protostar. At its formation, the first core has a thick disk-like structure, and is mainly supported by the thermal pressure. After the protostar formation, it decreases the thickness gradually, and becomes supported by the centrifugal force. We found that the first core is a precursor of the circumstellar disk. This indicates that the circumstellar disk is formed before the protostar formation with a size of ~10AU, which means that no protoplanetary disk smaller than <10AU exists. Reflecting the thermodynamics of the collapsing gas, at the protostar formation epoch, the circumstellar disk has a mass of ~0.01-0.1 solar mass, while the protostar has a mass of ~10^-3 solar mass. Thus, just after the protostar formation, the circumstellar disk is about 10-100 times more massive than the protostar. Even in the main accretion phase that lasts for ~10^5yr, the circumstellar disk mass dominates the protostellar mass. Such a massive disk is unstable to gravitational instability, and tends to show fragmentation. Our calculations indicate that the planet or brown-dwarf mass object may form in the circumstellar disk in the main accretion phase. In addition, the mass accretion rate onto the protostar shows strong time variability that is caused by the perturbation of proto-planets and/or the spiral arms in the circumstellar disk. Such variability provides a useful signature for detecting the planet-sized companion in the circumstellar disk around very young protostars.

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.

preprint2010arXiv

Population III Star Formation in Magnetized Primordial Clouds

The evolution of primordial collapsing clouds and formation of proto-Population III stars are investigated using three-dimensional ideal MHD simulation. We calculated the evolution of magnetized primordial clouds from the prestellar stage until the epoch after the proto-Population III star formation, spatially resolving both parsec-scale clouds and sub-AU scale protostars. The formation process of proto-population III star is characterized by the ratio of rotational to magnetic energy of the parent cloud. When the rotational energy is larger than the magnetic energy, fragmentation occurs in the collapsing primordial cloud before the proto-Population III star formation and binary or multiple system appears. Instead, when the magnetic energy is larger than the rotational energy, strong jet with >100km s^-1 is driven by circumstellar disk around the proto-population III star without fragmentation. Thus, even in the early universe, the magnetic field plays an important role in the star formation process.

preprint2010arXiv

Radiation Magnetohydrodynamics Simulation of Proto-Stellar Collapse: Two-Component Molecular Outflow

We perform a three-dimensional nested-grid radiation magneto-hydrodynamics (RMHD) simulation with self-gravity to study the early phase of the low-mass star formation process from a rotating molecular cloud core to a first adiabatic core just before the second collapse begins. Radiation transfer is handled with the flux-limited diffusion approximation, operator-splitting and implicit time-integrator. In the RMHD simulation, the outer region of the first core attains a higher entropy and the size of first core is larger than that in the magnetohydrodynamics simulations with the barotropic approximation. Bipolar molecular outflow consisting of two components is driven by magnetic Lorentz force via different mechanisms, and shock heating by the outflow is observed. Using the RMHD simulation we can predict and interpret the observed properties of star-forming clouds, first cores and outflows with millimeter/submillimeter radio interferometers, especially the Atacama Large Millimeter/submillimeter Array (ALMA).

preprint2010arXiv

Recent Developments in Simulations of Low-mass Star Formation

In star forming regions, we can observe different evolutionary stages of various objects and phenomena such as molecular clouds, protostellar jets and outflows, circumstellar disks, and protostars. However, it is difficult to directly observe the star formation process itself, because it is veiled by the dense infalling envelope. Numerical simulations can unveil the star formation process in the collapsing gas cloud. Recently, some studies showed protostar formation from the prestellar core stage, in which both molecular clouds and protostars are resolved with sufficient spatial resolution. These simulations showed fragmentation and binary formation, outflow and jet driving, and circumstellar disk formation in the collapsing gas clouds. In addition, the angular momentum transfer and dissipation process of the magnetic field in the star formation process were investigated. In this paper, I briefly review recent developments in numerical simulations of low-mass star formation.

preprint2010arXiv

The Origin and Formation of the Circumstellar Disk

The formation and evolution of the circumstellar disk in the collapsing molecular cloud is investigated from the prestellar stage resolving both the molecular cloud core and the protostar itself. In the collapsing cloud, the first adiabatic core appears prior to the protostar formation. Reflecting the thermodynamics of the collapsing gas, the first core is much more massive than the protostar. When the molecular cloud has no angular momentum, the first core falls onto the protostar and disappears a few years after the protostar formation. On the other hand, when the molecular cloud has an angular momentum, the first core does not disappear even after the protostar formation, and directly evolves into the circumstellar disk with a Keplerian rotation. There are two paths for the formation of the circumstellar disk. When the initial cloud has a considerably small rotational energy, two nested disks appear just after the protostar formation. During the early main accretion phase, the inner disk increases its size and merges with the outer disk (i.e. first core) to form a single circumstellar disk with a Keplerian rotation. On the other hand, when the molecular cloud has a rotational energy comparable to observations, a single Keplerian disk that corresponds to the first core already exists prior to the protostar formation. In such a cloud, the first core density gradually increases, maintaining the Keplerian rotation and forms the protostar inside it. Thus, the protostar is born in the Keplerian disk. In other words, a massive disk already exists before the protostar formation. In each case, the protostar at its formation is already surrounded by a massive circumstellar disk. The circumstellar disk is about 10-100 times more massive than the protostar in the main accretion disk. Such disks are favourable sites for the formation of binary companions and gas-giant planets.

preprint2009arXiv

Binary Formation with Different Metallicities: Dependence on Initial Conditions

The fragmentation process in collapsing clouds with various metallicities is studied using three-dimensional nested-grid hydrodynamics. Initial clouds are specified by three parameters: cloud metallicity, initial rotation energy and initial cloud shape. For different combinations of these parameters, we calculate 480 models in total and study cloud evolution, fragmentation conditions, orbital separation and binary frequency. For the cloud to fragment during collapse, the initial angular momentum must be higher than a threshold value, which decreases with decreasing metallicity. Although the exact fragmentation conditions depend also on the initial cloud shape, this dependence is only modest. Our results indicate a higher binary frequency in lower-metallicity gas. In particular, with the same median rotation parameter as in the solar neighbourhood, a majority of stars are born as members of binary/multiple systems for < 10^-4 Z_sun. With initial mass < 0.1 M_sun, if fragments are ejected in embryo from the host clouds by multi-body interaction, they evolve to substellar-mass objects. This provides a formation channel for low-mass stars in zero- or low-metallicity environments.

preprint2009arXiv

First Direct Simulation of Brown Dwarf Formation in a Compact Cloud Core

Brown dwarf formation and star formation efficiency are studied using a nested grid simulation that covers five orders of magnitude in spatial scale (10^4 - 0.1AU). Starting with a rotating magnetized compact cloud with a mass of 0.22 M_sun, we follow the cloud evolution until the end of main accretion phase. Outflow of about 5 km/s emerges about 100 yr before the protostar formation and does not disappear until the end of the calculation. The mass accretion rate declines from 10^-6 M_sun/yr to 10^-8 - 10^-12 M_sun/yr in a short time (about 10^4 yr) after the protostar formation. This is because (1) a large fraction of mass is ejected from the host cloud by the protostellar outflow and (2) the gas escapes from the host cloud by the thermal pressure. At the end of the calculation, 74% (167 M_Jup) of the total mass (225 M_Jup) is outflowing from the protostar, in which 34% (77 M_Jup) of the total mass is ejected by the protostellar outflow with supersonic velocity and 40% (90 M_Jup) escapes with subsonic velocity. On the other hand, 20% (45 M_Jup) is converted into the protostar and 6% (13 M_Jup) remains as the circumstellar disk. Thus, the star formation efficiency is epsilon = 0.2. The resultant protostellar mass is in the mass range of brown dwarfs. Our results indicate that brown dwarfs can be formed in compact cores in the same manner as hydrogen-burning stars, and the magnetic field and protostellar outflow are essential in determining the star formation efficiency and stellar mass.

preprint2009arXiv

Protostellar Jet and Outflow in the Collapsing Cloud Core

We investigate the driving mechanism of outflows and jets in star formation process using resistive MHD nested grid simulations. We found two distinct flows in the collapsing cloud core: Low-velocity outflows (sim 5 km/s) with a wide opening angle, driven from the first adiabatic core, and high-velocity jets (sim 50 km/s) with good collimation, driven from the protostar. High-velocity jets are enclosed by low-velocity outflow. The difference in the degree of collimation between the two flows is caused by the strength of the magnetic field and configuration of the magnetic field lines. The magnetic field around an adiabatic core is strong and has an hourglass configuration. Therefore, the low-velocity outflow from the adiabatic core are driven mainly by the magnetocentrifugal mechanism and guided by the hourglass-like field lines. In contrast, the magnetic field around the protostar is weak and has a straight configuration owing to Ohmic dissipation in the high-density gas region. Therefore, high-velocity jet from the protostar are driven mainly by the magnetic pressure gradient force and guided by straight field lines. Differing depth of the gravitational potential between the adiabatic core and the protostar cause the difference of the flow speed. Low-velocity outflows correspond to the observed molecular outflows, while high-velocity jets correspond to the observed optical jets. We suggest that the protostellar outflow and the jet are driven by different cores (the first adiabatic core and protostar), rather than that the outflow being entrained by the jet.

preprint2009arXiv

Star Formation in Relic HII Regions of the First Stars: Binarity and Outflow Driving

Star formation in relic HII regions of the first stars is investigated using magneto-hydrodynamical simulations with a nested grid method that covers 10 orders of magnitude in spatial scale and 20 orders of magnitude in density contrast. Due to larger fraction of H_2 and HD molecules, its prestellar thermal evolution is considerably different from that in the first star formation. Reflecting the difference, two hydrostatic cores appear in a nested manner: a protostar is enclosed by a transient hydrostatic core, which appears during the prestellar collapse. If the initial natal core rotates fast at a rate with rotational to gravitational energy ratio β_0 > 0.01-0.1, the transient hydrostatic core fragments to \sim 10 M_\odot sub-cores at density \sim 10^9 cm^-3. With smaller rotation energy, fragmentation occurs at higher density while a single protostar forms without fragmentation if rotation is extremely slow with β_0 < 10^-6 -10^-5. If magnetic field is present, these threshold values of β_0 is boosted owing to angular momentum transport by the magnetic breaking. Magnetic field also drives the protostellar outflows. With strong magnetic field, two distinct outflows are observed: The slower one emanates from the transient hydrostatic core, while the faster one from the protostar. These flows may affect the final stellar mass by ejecting some of masses in the initial core, and also may play some role in driving and maintenance of interstellar turbulence in young galaxies.

preprint2009arXiv

The Circumbinary Outflow: A Protostellar Outflow Driven by a Circumbinary Disk

The protostellar outflows have indispensable role in the formation of single stars, because they carry off the excess angular momentum from the centre of the shrinking gas cloud, and permits further collapse to form a star. On the other hand, a significant fraction of stars is supposed to be born as binaries with circumbinary disk that are frequently observed. Here, we investigate the evolution of a magnetized rotating cloud using three-dimensional resistive MHD nested-grid code, and show that the outflow is driven by the circumbinary disk and has an important role even in the binary formation. After the adiabatic core formation in the collapsing cloud core, the magnetic flux is significantly removed from the centre of the cloud by the Ohmic dissipation. Since this removal makes the magnetic braking ineffective, the adiabatic core continuously acquires the angular momentum to induce fragmentation and subsequent binary formation. The magnetic field accumulates in the circumbinary disk where the removal and accretion of magnetic field are balanced, and finally drives circumbinary outflow. This result explains the spectacular morphology of some specific young stellar objects such as L1551 IRS5. We can infer that most of the bipolar molecular outflows observed by low density tracers (i.e., CO) would correspond to circumbinary or circum-multiple outflows found in this report, since most of the young stellar objects are supposed to be binaries or multiples.

preprint2006arXiv

Outflows driven by Giant Protoplanets

We investigate outflows driven by a giant protoplanet using three-dimensional MHD nested grid simulations. We consider a local region around the protoplanet in the protoplanetary disk, and calculate three models: (a) unmagnetized disk model, (b) magnetized disk model having magnetic field azimuthally parallel to the disk, and (c) magnetic field perpendicular to the disk. Outflows with velocities, at least, 10 km/s are driven by the protoplanets in both magnetized disk models, while outflow does not appear in unmagnetized disk model. Tube-like outflows along the azimuthal direction of the protoplanetary disk appear in model with magnetic field being parallel to the disk. In this model, the magnetically dominated regions (i.e., density gap) are clearly contrasted from other regions and spiral waves appear near the protoplanet. On the other hand, in model with magnetic field being perpendicular to the disk, outflows are driven by a protoplanet with cone-like structure just as seen in the outflow driven by a protostar. Magnetic field lines are strongly twisted near the protoplanet and the outflows have well-collimated structures in this model.These outflows can be landmarks for searching exo-protoplanets in their formation stages. Our results indicate that the accretion rate onto the protoplanet tend to have a larger value than that expected from previous hydrodynamical calculations, since a fraction of the angular momentum of circum-planetary disk is removed by outflows, enhanced non-axisymmetric patterns caused by magnetic field, and magnetic braking. Possible implications for observation are also briefly discussed.