BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Tomoaki Matsumoto

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.GAastro-phastro-ph.EPastro-ph.COastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

42works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

42 published item(s)

preprint2022arXiv

Cluster Formation in GGD12-15: Infall Motion with Rotation of the Natal Clump

We report results of observations of the GGD12-15 region, where cluster formation is actively taking place, using various molecular emission lines. The C18O (J= 1-0) emission line reveals a massive clump in the region with a mass of ~2800 Mo distributed over ~2 pc. The distribution of the C18O(J= 3-2) emission is similar to that of a star cluster forming therein, with an elliptical shape of ~1 pc in size. A bipolar molecular outflow driven by IRS 9Mc, a constituent star of the cluster, is blowing in a direction perpendicular to the elongated C18O (J= 3-2) distribution, covering the entire clump. There is a massive core with a radius of 0.3 pc and a mass of 530 Mo in the center of the clump. There are two velocity components around the core, which are prominent in a position-velocity (PV) diagram along the major axis of the clump. In addition, a PV diagram along the minor axis of the clump, which is parallel to the outflow, shows a velocity gradient opposite to that of the outflow. The velocity structure strongly indicates the infalling motion of the clump. Comparison of the observational data with a simple model of infalling oblate clumps indicates that the clump is undergoing gravitational contraction with rotation.

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Circumbinary Disks of the Protostellar Binary Systems in the L1551 Region

We report ALMA Cycle 4 observations of the Class I binary protostellar system L1551 IRS 5 in the 0.9-mm continuum emission, C18O (J=3-2), OCS (J=28-27), and four other Band 7 lines. At ~0.07" (= 10 au) resolution in the 0.9 mm emission, two circumstellar disks (CSDs) associated with the binary protostars are separated from the circumbinary disk (CBD). The CBD is resolved into two spiral arms, one connecting to the CSD around the northern binary source, Source N, and the other to Source S. As compared to the CBD in the neighboring protobinary system L1551 NE, the CBD in L1551 IRS 5 is more compact (r ~150 au) and the m=1 mode of the spirals found in L1551 NE is less obvious in L1551 IRS 5. Furthermore, the dust and molecular-line brightness temperatures of CSDs and CBD reach >260 K and >100 K, respectively, in L1551 IRS 5, much hotter than those in L1551 NE. The gas motions in the spiral arms are characterized by rotation and expansion. Furthermore, the transitions from the CBD to the CSD rotations at around the L2 and L3 Lagrangian points and gas motions around the L1 point are identified. Our numerical simulations reproduce the observed two spiral arms and expanding gas motion as a result of gravitational torques from the binary, transitions from the CBD to the CSD rotations, and the gas motion around the L1 point. The higher temperature in L1551 IRS 5 likely reflects the inferred FU-Ori event.

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Distortion of Magnetic Fields in BHR 71

The magnetic field structure of a star-forming Bok globule BHR 71 was determined based on near-infrared polarimetric observations of background stars. The magnetic field in BHR 71 was mapped from 25 stars. By using a simple 2D parabolic function, the plane-of-sky magnetic axis of the core was found to be $θ_{\rm mag} = 125^{\circ} \pm 11^{\circ}$. The plane-of-sky mean magnetic field strength of BHR 71 was found to be $B_{\rm pos} = 8.8 - 15.0$ $μ$G, indicating that the BHR 71 core is magnetically supercritical with $λ= 1.44 - 2.43$. Taking into account the effect of thermal/turbulent pressure and the plane-of-sky magnetic field component, the critical mass of BHR 71 was $M_{\rm cr} = 14.5-18.7$ M$_{\odot}$, which is consistent with the observed core mass of $M_{\rm core} \approx 14.7$ M$_{\odot}$ (Yang et al. 2017). We conclude that BHR 71 is in a condition close to a kinematically critical state, and the magnetic field direction lies close to the plane of sky. Since BHR 71 is a star-forming core, a significantly subcritical condition (i.e., the magnetic field direction deviating from the plane of sky) is unlikely, and collapsed from a condition close to a kinematically critical state. There are two possible scenarios to explain the curved magnetic fields of BHR 71, one is an hourglass-like field structure due to mass accumulation and the other is the Inoue \& Fukui (2013) mechanism, which proposes the interaction of the core with a shock wave to create curved magnetic fields wrapping around the core.

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Star cluster formation and cloud dispersal by radiative feedback: dependence on metallicity and compactness

We study star cluster formation in various environments with different metallicities and column densities by performing a suite of three-dimensional radiation hydrodynamics simulations. We find that the photoionization feedback from massive stars controls the star formation efficiency (SFE) in a star-forming cloud, and its impact sensitively depends on the gas metallicity $Z$ and initial cloud surface density $Σ$. At $Z=1~Z_{\odot}$, SFE increases as a power law from 0.03 at $Σ= 10~M_{\odot}{\rm pc^{-2}}$ to 0.3 at $Σ= 300~M_{\odot}{\rm pc^{-2}}$. In low-metallicity cases $10^{-2}- 10^{-1} Z_{\odot}$, star clusters form from atomic warm gases because the molecule formation time is not short enough with respect to the cooling or dynamical time. In addition, the whole cloud is disrupted more easily by expanding H{\sc ii} bubbles which have higher temperature owing to less efficient cooling. With smaller dust attenuation, the ionizing radiation feedback from nearby massive stars is stronger and terminate star formation in dense clumps. These effects result in inefficient star formation in low-metallicity environments: the SFE drops by a factor of $\sim 3$ at $Z=10^{-2}~Z_{\odot}$ compared to the results for $Z=1~Z_{\odot}$, regardless of $Σ$. Newborn star clusters are also gravitationally less bound. We further develop a new semi-analytical model that can reproduce the simulation results well, particularly the observed dependencies of the SFEs on the cloud surface densities and metallicities.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Birth of a Massive First-Star Binary

We study the formation of massive Population III binary stars using a newly developed radiation hydrodynamics code with the adaptive mesh refinement and adaptive ray-tracing methods. We follow the evolution of a typical primordial star-forming cloud obtained from a cosmological hydrodynamics simulation. Several protostars form as a result of disk fragmentation and grow in mass by the gas accretion, which is finally quenched by the radiation feedback from the protostars. Our code enables us, for the first time, to consider the feedback by both the ionizing and dissociating radiation from the multiple protostars, which is essential for self-consistently determining their final masses. At the final step of the simulation, we observe a very wide ($\gtrsim 10^4\,\mathrm{au}$) binary stellar system consisting of $60$ and $70\,M_\odot$ stars. One of the member stars also has two smaller mass ($10\,M_\odot$) companion stars orbiting at $200$ and $800\,\mathrm{au}$, making up a mini-triplet system. Our results suggest that massive binary or multiple systems are common among Population III stars.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Distortion of Magnetic Fields in a Starless Core VI: Application of Flux Freezing Model and Core Formation of FeSt 1-457

Observational data for the hourglass-like magnetic field toward the starless dense core FeSt 1-457 were compared with a flux freezing magnetic field model (Myers et al. 2018). Fitting of the observed plane-of-sky magnetic field using the flux freezing model gave a residual angle dispersion comparable with the results based on a simple three-dimensional parabolic model. The best-fit parameters for the flux freezing model were a line-of-sight magnetic inclination angle of $γ_{\rm mag} = 35^{\circ} \pm 15^{\circ}$ and a core center to ambient (background) density contrast of $ρ_{\rm c} / ρ_{\rm bkg} = 75$. The initial density for core formation ($ρ_0$) was estimated to be $ρ_{\rm c} / 75 = 4670$ cm$^{-3}$, which is about one order of magnitude higher than the expected density ($\sim 300$ cm$^{-3}$) for the inter-clump medium of the Pipe Nebula. FeSt 1-457 is likely to have been formed from the accumulation of relatively dense gas, and the relatively dense background column density of $A_V \simeq 5$ mag supports this scenario. The initial radius (core formation radius) $R_0$ and the initial magnetic field strength $B_0$ were obtained to be 0.15 pc ($1.64 R$) and $10.8-14.6$ $μ$G, respectively. We found that the initial density $ρ_0$ is consistent with the mean density of the nearly critical magnetized filament with magnetic field strength $B_0$ and radius $R_0$. The relatively dense initial condition for core formation can be naturally understood if the origin of the core is the fragmentation of magnetized filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Discovery of Infalling Motion with Rotation of the Cluster-Forming Clump S235AB and Its Implication to the Clump Structures

We report the discovery of infalling motion with rotation of S235AB the massive cluster-forming clump (~10^3 Mo) in the S235 region. Our C18O observations with the 45m telescope at the Nobeyama Radio Observatory have revealed the elliptical shape of the clump. Position-velocity (PV) diagram taken along its major axis exhibits two well-defined peaks symmetrically located with respect to the clump center, which is similar to that found for a dynamically infalling envelope with rotation around a single protostar modeled by N. Ohashi and his collaborators, indicating that the cluster-forming clump is also collapsing by the self-gravity toward the clump center. With analogue to Ohashi's model, we made a simple model of an infalling, rotating clump to fit the observed data. Based on the inferred model parameters as well as results of earlier observations and simulations in the literature, we discuss structures of the clump such as the relation among the global mass infall rate (~10^-3 Mo/yr), formation of a compact core (with a mass and size of ~4 Mo and < 0.1 pc) at the center, and a massive star (~11 Mo) forming in the core.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation of the Unequal-Mass Binary Protostars in L1551 NE by Rotationally-Driven Fragmentation

We present observations at 7 mm that fully resolve the two circumstellar disks, and a reanalyses of archival observations at 3.5 cm that resolve along their major axes the two ionized jets, of the class I binary protostellar system L1551 NE. We show that the two circumstellar disks are better fit by a shallow inner and steep outer power-law than a truncated power-law. The two disks have very different transition radii between their inner and outer regions of $\sim$18.6 AU and $\sim$8.9 AU respectively. Assuming that they are intrinsically circular and geometrically thin, we find that the two circumstellar disks are parallel with each other and orthogonal in projection to their respective ionized jets. Furthermore, the two disks are closely aligned if not parallel with their circumbinary disk. Over an interval of $\sim$10 yr, source B (possessing the circumsecondary disk) has moved northwards with respect to and likely away from source A, indicating an orbital motion in the same direction as the rotational motion of their circumbinary disk. All the aforementioned elements therefore share the same axis for their angular momentum, indicating that L1551 NE is a product of rotationally-driven fragmentation of its parental core. Assuming a circular orbit, the relative disk sizes are compatible with theoretical predictions for tidal truncation by a binary system having a mass ratio of $\sim$0.2, in agreement with the reported relative separations of the two protostars from the center of their circumbinary disk. The transition radii of both disks, however, are a factor of $\gtrsim$1.5 smaller than their predicted tidally-truncated radii.

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rotationally-Driven Fragmentation for the Formation of the Binary Protostellar System L1551 IRS 5

Either bulk rotation or local turbulence is widely invoked to drive fragmentation in collapsing cores so as to produce multiple star systems. Even when the two mechanisms predict different manners in which the stellar spins and orbits are aligned, subsequent internal or external interactions can drive multiple systems towards or away from alignment thus masking their formation process. Here, we demonstrate that the geometrical and dynamical relationship between the binary system and its surrounding bulk envelope provide the crucial distinction between fragmentation models. We find that the circumstellar disks of the binary protostellar system L1551 IRS 5 are closely parallel not just with each other but also with their surrounding flattened envelope. Measurements of the relative proper motion of the binary components spanning nearly 30 yr indicate an orbital motion in the same sense as the envelope rotation. Eliminating orbital solutions whereby the circumstellar disks would be tidally truncated to sizes smaller than are observed, the remaining solutions favor a circular or low-eccentricity orbit tilted by up to $\sim$25$^\circ$ from the circumstellar disks. Turbulence-driven fragmentation can generate local angular momentum to produce a coplanar binary system, but which bears no particular relationship with its surrounding envelope. Instead, the observed properties conform with predictions for rotationally-driven fragmentation. If the fragments were produced at different heights or on opposite sides of the midplane in the flattened central region of a rotating core, the resulting protostars would then exhibit circumstellar disks parallel with the surrounding envelope but tilted from the orbital plane as is observed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An origin of arc structures deeply embedded in dense molecular cloud cores

We investigated the formation of arc-like structures in the infalling envelope around protostars, motivated by the recent Atacama Large Millimeter/Submillimeter Array (ALMA) observations of the high-density molecular cloud core, MC27/L1521F. We performed self-gravitational hydrodynamical numerical simulations with an adaptive mesh refinement code. A filamentary cloud with a 0.1~pc width fragments into cloud cores because of perturbations due to weak turbulence. The cloud core undergoes gravitational collapse to form multiple protostars, and gravitational torque from the orbiting protostars produces arc structures extending up to a 1000~AU scale. As well as on a spatial extent, the velocity ranges of the arc structures, $\sim0.5\,\mathrm{km\,s}^{-1}$, are in agreement with the ALMA observations. We also found that circumstellar disks are often misaligned in triple system. The misalignment is caused by the tidal interaction between the protostars when they undergo close encounters because of a highly eccentric orbit of the tight binary pair.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Star formation in turbulent molecular clouds with colliding flow

Using self-gravitational hydrodynamical numerical simulations, we investigated the evolution of high-density turbulent molecular clouds swept by a colliding flow. The interaction of shock waves due to turbulence produces networks of thin filamentary clouds with a sub-parsec width. The colliding flow accumulates the filamentary clouds into a sheet cloud and promotes active star formation for initially high-density clouds. Clouds with a colliding flow exhibit a finer filamentary network than clouds without a colliding flow. The probability distribution functions (PDFs) for the density and column density can be fitted by lognormal functions for clouds without colliding flow. When the initial turbulence is weak, the column density PDF has a power-law wing at high column densities. The colliding flow considerably deforms the PDF, such that the PDF exhibits a double peak. The stellar mass distributions reproduced here are consistent with the classical initial mass function with a power-law index of $-1.35$ when the initial clouds have a high density. The distribution of stellar velocities agrees with the gas velocity distribution, which can be fitted by Gaussian functions for clouds without colliding flow. For clouds with colliding flow, the velocity dispersion of gas tends to be larger than the stellar velocity dispersion. The signatures of colliding flows and turbulence appear in channel maps reconstructed from the simulation data. Clouds without colliding flow exhibit a cloud-scale velocity shear due to the turbulence. In contrast, clouds with colliding flow show a prominent anti-correlated distribution of thin filaments between the different velocity channels, suggesting collisions between the filamentary clouds.

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Angular Momentum Exchange by Gravitational Torques and Infall in the Circumbinary Disk of the Protostellar System L1551 NE

We report the ALMA observation of the Class I binary protostellar system L1551 NE in the 0.9-mm continuum, C18O (3-2), and 13CO (3-2) lines at a ~1.6 times higher resolution and a ~6 times higher sensitivity than those of our previous SMA observations, which revealed a r ~300 AU-scale circumbinary disk in Keplerian rotation. The 0.9-mm continuum shows two opposing U-shaped brightenings in the circumbinary disk, and exhibits a depression between the circumbinary disk and the circumstellar disk of the primary protostar. The molecular lines trace non-axisymmetric deviations from Keplerian rotation in the circumbinary disk at higher velocities relative to the systemic velocity, where our previous SMA observations could not detect the lines. In addition, we detect inward motion along the minor axis of the circumbinary disk. To explain the newly-observed features, we performed a numerical simulation of gas orbits in a Roche potential tailored to the inferred properties of L1551 NE. The observed U-shaped dust features coincide with locations where gravitational torques from the central binary system are predicted to impart angular momentum to the circumbinary disk, producing shocks and hence density enhancements seen as a pair of spiral arms. The observed inward gas motion coincides with locations where angular momentum is predicted to be lowered by the gravitational torques. The good agreement between our observation and model indicates that gravitational torques from the binary stars constitute the primary driver for exchanging angular momentum so as to permit infall through the circumbinary disk of L1551 NE.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Colliding Filaments and a Massive Dense Core in the Cygnus OB 7 Molecular Cloud

We report results of molecular line observations carried out toward a massive dense core in the Cyg OB 7 molecular cloud. The core has an extraordinarily large mass ($\sim1.1 \times 10^4$ $M_\odot$) and size ($\sim2 \times 5$ pc$^2$), but there is no massive young star forming therein. We observed this core in various molecular lines such as C$^{18}$O($J=1-0$) using the 45m telescope at Nobeyama Radio Observatory. We find that the core has an elongated morphology consisting of several filaments and core-like structures. The filaments are massive ($10^2-10^3$ $M_\odot$), and they are apparently colliding against each other. Some candidates of YSOs are distributed around their intersection, suggesting that the collisions of the filaments may have influenced on their formation. To understand the formation and evolution of such colliding filaments, we performed numerical simulations using the adaptive mesh refinement (AMR) technique adopting the observed core parameters (e.g., the mass and size) as the initial conditions. Results indicate that the filaments are formed as seen in other earlier simulations for small cores in literature, but we could not reproduce the collisions of the filaments simply by assuming the large initial mass and size. We find that the collisions of the filaments occur only when there is a large velocity gradient in the initial core in a sense to compress it. We suggest that the observed core was actually compressed by an external effect, e.g., shocks of nearby supernova remnants including HB21 which has been suggested to be interacting with the Cyg OB 7 molecular cloud.

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Chemistry in the First Hydrostatic Core Stage By Adopting Three-Dimensional Radiation Hydrodynamic Simulations

We investigate molecular evolution from a molecular cloud core to a first hydrostatic core in three spatial dimensions. We perform a radiation hydrodynamic simulation in order to trace fluid parcels, in which molecular evolution is investigated, using a gas-phase and grain-surface chemical reaction network. We derive spatial distributions of molecular abundances and column densities in the core harboring the first core. We find that the total of gas and ice abundances of many species in a cold era (10 K) remain unaltered until the temperature reaches ~500 K. The gas abundances in the warm envelope and the outer layer of the first core (T < 500 K) are mainly determined via the sublimation of ice-mantle species. Above 500 K, the abundant molecules, such as H2CO, start to be destroyed, and simple molecules, such as CO, H2O and N2 are reformed. On the other hand, some molecules are effectively formed at high temperature; carbon-chains, such as C2H2 and cyanopolyynes, are formed at the temperature of >700 K. We also find that large organic molecules, such as CH3OH and HCOOCH3, are associated with the first core (r < 10 AU). Although the abundances of these molecules in the first core stage are comparable or less than in the protostellar stage (hot corino), reflecting the lower luminosity of the central object, their column densities in our model are comparable to the observed values toward the prototypical hot corino, IRAS 16293-2422. We propose that these large organic molecules can be good tracers of the first cores.

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An Implicit Scheme for Ohmic Dissipation with Adaptive Mesh Refinement

An implicit method for the ohmic dissipation is proposed. The proposed method is based on the Crank-Nicolson method and exhibits second-order accuracy in time and space. The proposed method has been implemented in the SFUMATO adaptive mesh refinement (AMR) code. The multigrid method on the grids of the AMR hierarchy converges the solution. The convergence is fast but depends on the time step, resolution, and resistivity. Test problems demonstrated that decent solutions are obtained even at the interface between fine and coarse grids. Moreover, the solution obtained by the proposed method shows good agreement with that obtained by the explicit method, which required many time steps. The present method reduces the number of time steps, and hence the computational costs, as compared with the explicit method.

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

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Direct Imaging of Bridged Twin Protoplanetary Disks in a Young Multiple Star

Studies of the structure and evolution of protoplanetary disks are important for understanding star and planet formation. Here, we present the direct image of an interacting binary protoplanetary system. Both circumprimary and circumsecondary disks are resolved in the near-infrared. There is a bridge of infrared emission connecting the two disks and a long spiral arm extending from the circumprimary disk. Numerical simulations show that the bridge corresponds to gas flow and a shock wave caused by the collision of gas rotating around the primary and secondary stars. Fresh material streams along the spiral arm, consistent with the theoretical scenarios where gas is replenished from a circummultiple reservoir.

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Protostellar collapse of magneto-turbulent cloud cores: shape during collapse and outflow formation

We investigate protostellar collapse of molecular cloud cores by numerical simulations, taking into account turbulence and magnetic fields. By using the adaptive mesh refinement technique, the collapse is followed over a wide dynamic range from the scale of a turbulent cloud core to that of the first core. The cloud core is lumpy in the low density region owing to the turbulence, while it has a smooth density distribution in the dense region produced by the collapse. The shape of the dense region depends mainly on the mass of the cloud core; a massive cloud core tends to be prolate while a less massive cloud core tends to be oblate. In both cases, anisotropy of the dense region increases during the isothermal collapse. The minor axis of the dense region is always oriented parallel to the local magnetic field. All the models eventually yield spherical first cores supported mainly by the thermal pressure. Most of turbulent cloud cores exhibit protostellar outflows around the first cores. These outflows are classified into two types, bipolar and spiral flows, according to the morphology of the associated magnetic field. Bipolar flow often appears in the less massive cloud core. The rotation axis of the first core is oriented parallel to the local magnetic field for bipolar flow, while the orientation of the rotation axis from the global magnetic field depends on the magnetic field strength. In spiral flow, the rotation axis is not aligned with the local magnetic field.

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

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

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

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

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Near-Infrared Imaging Polarimetry of the NGC 2071 Star Forming Region with SIRPOL

We have conducted deep JHKs imaging polarimetry of a ~8' x 8' area of the NGC 2071 star forming region. Our polarization data have revealed various infrared reflection nebulae (IRNe) associated with the central IR young star cluster NGC2071IR and identified their illuminating sources. There are at least 4 IRNe in NGC2071IR and several additional IRNe are identified around nearby young stars in the same field-of-view. Each illuminating source coincides with a known near-IR source except for IRS3, which is only a part of IRN2 and is illuminated by the radio source 1c. Aperture polarimetry of each cluster source is used to detect unresolved circumstellar disk/outflow systems. Aperture polarimetry of the other point-like sources within the field is made in this region for the first time. The magnetic field structures (from ~1 pc down to \~0.1 pc) are derived using both aperture polarimetry of the point-like sources and imaging polarimetry of the shocked H2 emission that is seen as the dominant knotty nebulae in the Ks band image; they are both of dichroic origin and the derived field directions are consistent with each other. The magnetic field direction projected on the sky is also consistent with that inferred from the 850 micron thermal continuum emission polarimetry of the central 0.2 pc region, but running roughly perpendicular (~75 degrees) to the direction of the large scale outflow. We argue that the field strength is too weak to align the outflow in the large scale field direction via magnetic braking.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Self-gravitational Magnetohydrodynamics with Adaptive Mesh Refinement for Protostellar Collapse

A new numerical code, called SFUMATO, for solving self-gravitational magnetohydrodynamics (MHD) problems using adaptive mesh refinement (AMR) is presented. A block-structured grid is adopted as the grid of the AMR hierarchy. The total variation diminishing (TVD) cell-centered scheme is adopted as the MHD solver, with hyperbolic cleaning of divergence error of the magnetic field also implemented. The self-gravity is solved by a multigrid method composed of (1) full multigrid (FMG)-cycle on the AMR hierarchical grids, (2) V-cycle on these grids, and (3) FMG-cycle on the base grid. The multigrid method exhibits spatial second-order accuracy, fast convergence, and scalability. The numerical fluxes are conserved by using a refluxing procedure in both the MHD solver and the multigrid method. The several tests are performed indicating that the solutions are consistent with previously published results.

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

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Stability of Dynamically Collapsing Gas Sphere

We discuss stability of dynamically collapsing gas spheres. We use a similarity solution for a dynamically collapsing sphere as the unperturbed state. In the similarity solution the gas pressure is approximated by a polytrope of $ P = K ρ^γ$. We examine three types of perturbations: bar ($ \ell = 2$) mode, spin-up mode, and Ori-Piran mode. When $ γ< 1.097 $, it is unstable against bar-mode. It is unstable against spin-up mode for any $ γ$. When $ γ< 0.961 $, the similarity solution is unstable against Ori-Piran mode. The unstable mode grows in proportion to $ | t - t_0 | ^{-σ} $ while the central density increases in proportion to $ ρ_c \propto (t - t_0) ^{-2} $ in the similarity solution. The growth rate, $ σ$ is obtained numerically as a function of $ γ$ for bar mode and Ori-Piran mode. The growth rate of the bar mode is larger for a smaller $ γ$. The spin-up mode has the growth rate of $ σ= 1/3 $ for any $ γ$.

0 citations0 reviewsNo code linkedOpen access