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

Kohji Tomisaka

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

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

Catalog footprint

What is connected

28works
6topics
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

28 published item(s)

preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mushroom-instability-driven Magnetic Reconnections in Collisionless Relativistic Jets

We study the kinetic plasma dynamics in collisionless relativistic jets with velocity shear, by carrying out particle-in-cell simulations in the transverse plane of a jet. It is discovered that intermittent magnetic reconnections (MRs) are driven by Mushroom instability (MI), which is an important kinetic-scale plasma instability in the plasma shear-flows with relativistic bulk speed. We refer to this sequence of kinetic plasma phenomena as "MI-driven MR". The MI-driven MRs intermittently occur with moving the location of the reconnection points from the vicinity of the initial velocity-shear surface towards the center of the jet. As a consequence, the number density of high energy electrons, which are accelerated by MI-driven MRs, increases with time in the region inside the initial velocity-shear surface with accompanying the generation and subsequent amplification of magnetic fields by MI. The maximum Lorentz factor of electrons increases with initial bulk Lorentz factor of the jet. A possible relation of MI-driven MR to the bright synchrotron emission in jet-spine of active galactic nucleus jets is also discussed.

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

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
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
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Polarization Structure of Filamentary Clouds

Filaments are considered to be basic structures and molecular clouds consist of filaments. Filaments are often observed as extending in the direction perpendicular to the interstellar magnetic field. The structure of filaments has been studied based on a magnetohydrostatic equilibrium model (Tomisaka 2014). Here, we simulate the expected polarization pattern for isothermal magnetohydrostatic filaments. The filament exhibits a polarization pattern in which the magnetic field is apparently perpendicular to the filament when observed from the direction perpendicular to the magnetic field. When the line-of-sight is parallel to the global magnetic field, the observed polarization pattern is dependent on the center-to-surface density ratio for the filament and the concentration of the gas mass toward the central magnetic flux tube. Filaments with low center-to-surface density ratios have an insignificant degree of polarization when observed from the direction parallel to the global magnetic field. However, models with a large center-to-surface density ratio have polarization patterns that indicate the filament is perpendicularly threaded by the magnetic field. When mass is heavily concentrated at the central magnetic flux tube, which can be realized by the ambipolar diffusion process, the polarization pattern is similar to that expected for a low center-to-surface density contrast.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Structure and Stability of Filamentary Clouds Supported by Lateral Magnetic Field

We have constructed two types of analytical models for an isothermal filamentary cloud supported mainly by magnetic tension. The first one describes an isolated cloud while the second considers filamentary clouds spaced periodically. Both the models assume that the filamentary clouds are highly flattened. The former is proved to be the asymptotic limit of the latter in which each filamentary cloud is much thinner than the distance to the neighboring filaments. We show that these models reproduce main features of the 2D equilibrium model of Tomisaka (2014) for filamentary cloud threaded by perpendicular magnetic field. It is also shown that the critical mass to flux ratio is $ M /Φ= (2 π\sqrt{G}) ^{-1} $, where $ M $, $ Φ$ and $ G $ denote the cloud mass, the total magnetic flux of the cloud, and the gravitational constant, respectively. This upper bound coincides with that for an axisymmetric cloud supported by poloidal magnetic fields. We applied the variational principle for studying the Jeans instability of the first model. Our model cloud is unstable against fragmentation as well as the filamentary clouds threaded by longitudinal magnetic field. The fastest growing mode has a wavelength several times longer than the cloud diameter. The second model describes quasi-static evolution of filamentary molecular cloud by ambipolar diffusion.

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetohydrostatic Equilibrium Structure and Mass of Filamentary Isothermal Cloud Threaded by Lateral Magnetic Field

Herschel observation has recently revealed that interstellar molecular clouds consist of many filaments. Polarization observations in optical and infrared wavelengths indicate that the magnetic field often runs perpendicular to the filament. In this paper, the magnetohydrostatic configuration of isothermal gas is studied, in which the thermal pressure and the Lorentz force are balanced against the self-gravity and the magnetic field is globally perpendicular to the axis of the filament. The model is controlled by three parameters: center-to-surface density ratio (ρ_c/ρ_s), plasma βof surrounding interstellar gas (β_0) and the radius of the hypothetical parent cloud normalized by the scale-height (R'_0), although there remains a freedom how the mass is distributed against the magnetic flux (mass loading). In the case that R'_0 is small enough, the magnetic field plays a role in confining the gas. However, the magnetic field generally has an effect in supporting the cloud. There is a maximum line-mass (mass per unit length) above which the cloud is not supported against the gravity. Compared with the maximum line-mass of non-magnetized cloud (2c_s^2/G, where c_s and G represent respectively the isothermal sound speed and the gravitational constant), that of the magnetized filament is larger than the non-magnetized one. The maximum line-mass is numerically obtained as λ_max \simeq 0.24 Φ_cl/G^1/2 + 1.66 c_s^2/G, where Φ_cl represents one half of the magnetic flux threading the filament per unit length. The maximum mass of the filamentary cloud is shown to be significantly affected by the magnetic field when the magnetic flux per unit length exceeds Φ_cl \gtrsim 3 pc μG (c_s/190 m s^-1)^2.

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

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.

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

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

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

Observational Identification of First Cores: Non-LTE Radiative Transfer Simulation

A first core is a first hydrostatic object formed in the course of dynamical contraction of a molecular cloud core. Since the inflow pattern changes drastically before and after the first core formation, it is regarded as a milestone in the star formation process. In order to identify the first core from a mapping observation, the features expected for the first core are studied for CS rotation transitions at radio wavelengths. The non-LTE radiation transfer is calculated for the results of radiation magnetohydrodynamical simulations of the contraction of the magnetized molecular cloud core in rotation (Tomida et al. 2010a). We use the Monte-Carlo method to solve the non-LTE radiation transfer in a nested grid hierarchy. In the first core phase, an outflow arises from the vicinity of the first core due to the twisted magnetic field amplified by the rotation motion of the contracting gas disk. The disk and outflow system has several characteristic observational features: (i) in the pole-on view, relatively opaque lines indicate asymmetry in the emission lines in which the blue side is stronger than the red side (an infall signature of the envelope); (ii) in the edge-on view, the disk has a signature of simultaneous rotation and infall, i.e., the integrated intensity of the reaching side is brighter than that of the receding side and the gradient in the intensity-weighted velocity is larger in the reaching side; (iii) the observed outflow indicates rotation around the rotation axis. The size of the outflow gives the approximate age after the first core is formed, since the outflow is not expected for the earlier runaway isothermal collapse phase.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Origin of Molecular Outflow Determined from Thermal Dust Polarization

The observational expectation of polarization measurements of thermal dust radiation is investigated to find information on molecular outflows based on magnetohydrodynamical (MHD) and radiation transfer simulations. There are two major proposed models for the driving of molecular outflows: (1) molecular gas is accelerated by a magnetic pressure gradient or magnetocentrifugal wind mechanism before the magnetic field and molecular gas are decoupled, (2) the linear momentum of a highly collimated jet is transferred to the ambient molecular gas. In order to distinguish between these two models, it is crucial to observe the configuration of the magnetic field. An observation of a toroidal magnetic field is strong evidence that the first of the models is appropriate. In this paper, we calculated the polarization distribution of thermal dust radiation due to the alignment of dust grains along the magnetic field using molecular outflow data calculated by two-dimensional axisymmetric MHD simulations. An asymmetric distribution around the z-axis is characteristic for magnetic fields composed of both poloidal and toroidal components. We determined that the outflow has a low polarization degree compared with the envelope and that the envelope and outflow have different polarization directions (B-vector), namely, the magnetic field within the envelope is parallel to the global magnetic field lines while the magnetic field of the outflow is perpendicular to it. Thus we have demonstrated that the point-symmetric (rather than axisymmetric) distributions of low polarization regions indicate that molecular outflows are likely to be magnetically driven. Observations of this polarization distribution with tools such as ALMA would confirm the origin of the molecular outflow.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Origin of multiple nuclei in ultraluminous infrared galaxies

Ultraluminous infrared galaxies (ULIRGs) with multiple ($\ge 3$) nuclei are frequently observed. It has been suggested that these nuclei are produced by multiple major mergers of galaxies. The expected rate of such mergers is, however, too low to reproduce the observed number of ULIRGs with multiple nuclei. We have performed high-resolution simulations of the merging of two gas-rich disk galaxies. We found that extremely massive and compact star clusters form from the strongly disturbed gas disks after the first or second encounter between the galaxies. The mass of such clusters reaches $\sim 10^8 M_{\odot}$, and their half-mass radii are $20-30 \rm{pc}$. Since these clusters consist of young stars, they appear to be several bright cores in the galactic central region ($\sim \rm{kpc}$). The peak luminosity of these clusters reaches $\sim 10%$ of the total luminosity of the merging galaxy. These massive and compact clusters are consistent with the characteristics of the observed multiple nuclei in ULIRGs. Multiple mergers are not necessary to explain multiple nuclei in ULIRGs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Shock-induced star cluster formation in colliding galaxies

We studied the formation process of star clusters using high-resolution N-body/smoothed particle hydrodynamcs simulations of colliding galaxies. The total number of particles is 1.2x10^8 for our high resolution run. The gravitational softening is 5 pc and we allow gas to cool down to \sim 10 K. During the first encounter of the collision, a giant filament consists of cold and dense gas found between the progenitors by shock compression. A vigorous starburst took place in the filament, resulting in the formation of star clusters. The mass of these star clusters ranges from 10^{5-8} Msun. These star clusters formed hierarchically: at first small star clusters formed, and then they merged via gravity, resulting in larger star clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Generation of Seed Magnetic Fields in Primordial Supernova Remnants

Origin of the magnetic field ubiquitous in the Universe is studied based on the Biermann mechanism, which is expected to work in the non-barotropic region. We perform a series of two-dimensional MHD simulations of the first generation supernova remnant (SNR) expanding in the inhomogeneous interstellar matter (ISM) and study the Biermann mechanism working in the interior of the SNR. Especially, we pay attention to the relaxation process of electron and ion temperatures via the Coulomb interaction. In the early SNR in which the electron temperature is much lower than the ion temperature, the Biermann mechanism is ineffective, since the gradient of electron pressure is small. Magnetic fields begin to be generated just behind the shock front when the electron temperature is sufficiently relaxed. Assuming the explosion energy of 10^52 erg, the total magnetic energy generated reaches about 10^26 erg and does not depend strongly on the parameters of either SNR or ISM. Analytic expression to estimate the magnetic total energy is presented and it is shown this agrees well with our numerical results. Finally we evaluate the expected amplitude of magnetic fields in protogalaxies as ~10^-19 G, which is sufficient for seed fields of the observed galactic magnetic fields.

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
preprint2009arXiv

Formation of Galactic Center Magnetic Loops

A survey for the molecular clouds in the Galaxy with NANTEN mm telescope has discovered molecular loops in the Galactic center region. The loops show monotonic gradients of the line of sight velocity along the loops and the large velocity dispersions towards their foot points. It is suggested that these loops are explained in terms of the buoyant rise of magnetic loops due to the Parker instability. We have carried out global three-dimensional magneto-hydrodynamic simulations of the gas disk in the Galactic center. The gravitational potential is approximated by the axisymmetric potential proposed by Miyamoto & Nagai (1975). At the initial state, we assume a warm (~ 10^4 K) gas torus threaded by azimuthal magnetic fields. Self-gravity and radiative cooling of the gas are ignored. We found that buoyantly rising magnetic loops are formed above the differentially rotating, magnetically turbulent disk. By analyzing the results of global MHD simulations, we have identified individual loops, about 180 in the upper half of the disk, and studied their statistical properties such as their length, width, height, and velocity distributions along the loops. Typical length and height of a loop are 1kpc and 200pc, respectively. The line of sight velocity changes linearly along a loop and shows large dispersions around the foot-points. Numerical results indicate that loops emerge preferentially from the region where magnetic pressure is large. We argue that these properties are consistent with those of the molecular loops discovered by NANTEN.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Toward First-Principle Simulations of Galaxy Formation: II. Shock-Induced Starburst at a Collision Interface During the First Encounter of Interacting Galaxies

We investigated the evolution of interacting disk galaxies using high-resolution $N$-body/SPH simulations, taking into account the multiphase nature of the interstellar medium (ISM). In our high-resolution simulations, a large-scale starburst occurred naturally at the collision interface between two gas disks at the first encounter, resulting in the formation of star clusters. This is consistent with observations of interacting galaxies. The probability distribution function (PDF) of gas density showed clear change during the galaxy-galaxy encounter. The compression of gas at the collision interface between the gas disks first appears as an excess at $n_{\rm H} \sim 10{\rm cm^{-3}}$ in the PDF, and then the excess moves to higher densities ($n_{\rm H} \gtrsim 100{\rm cm^{-3}}$) in a few times $10^7$ years where starburst takes place. After the starburst, the PDF goes back to the quasi-steady state. These results give a simple picture of starburst phenomena in galaxy-galaxy encounters.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Toward first-principle simulations of galaxy formation: I. How should we choose star formation criteria in high-resolution simulations of disk galaxies?

We performed 3-dimensional N-body/SPH simulations to study how mass resolution and other model parameters such as the star formation efficiency parameter, C* and the threshold density, nth affect structures of the galactic gaseous/stellar disk in a static galactic potential. We employ 10^6 - 10^7 particles to resolve a cold and dense (T < 100 K & n_H > 100 cm^{-3}) phase. We found that structures of the ISM and the distribution of young stars are sensitive to the assumed nth. High-nth models with nth = 100 cm^{-3} yield clumpy multi-phase features in the ISM. Young stars are distributed in a thin disk of which half-mass scale height is 10 - 30 pc. In low-nth models with nth = 0.1 cm^{-3}, the stellar disk is found to be several times thicker, and the gas disk appears smoother than the high-nth models. A high-resolution simulation with high-nth is necessary to reproduce the complex structure of the gas disk. The global properties of the model galaxies in low-nth models, such as star formation histories, are similar to those in the high-nth models when we tune the value of C* so that they reproduce the observed relation between surface gas density and surface star formation rate density. We however emphasize that high-nth models automatically reproduce the relation, regardless of the values of C*. The ISM structure, phase distribution, and distributions of young star forming region are quite similar between two runs with values of C* which differ by a factor of 15. We also found that the timescale of the flow from n_H ~1 cm^{-3} to n_H > 100 cm^{-3} is about 5 times as long as the local dynamical time and is independent of the value of C*. The use of a high-nth criterion for star formation in high-resolution simulations makes numerical models fairy insensitive to the modelling of star formation. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Accretion in Gravitationally Contracting Clouds

Accretion flow in a contracting magnetized isothermal cloud was studied using magnetohydrodynamical simulations and a nested grid technique. First, the interstellar magnetized cloud experiences a ``runaway collapse'' phase, in which the central density increases drastically within a finite time scale. Finally, it enters an accretion phase, in which inflowing matter accretes onto a central high-density disk or a new-born star. We found that the accretion rate reaches (4 -- 40) $\times c_s^3/G$, where $c_s$ and $G$ represent the isothermal sound speed and the gravitational constant, respectively. This is much larger than the standard accretion rate of $0.975c_s^3/G$ for a hydrostatic isothermal spherical cloud (Shu 1977, AAA19.065.044). Due to the effect of an extra infall velocity achieved in the runaway phase ($\sim 2 c_s$), the accretion rate is boosted. This rate declines with time in contrast to Shu's solution, but keeps $\gtsim 2.5 c_s^3/G$. The observed gas infall rate around proto-stars such as L1551 IRS 5 and HL Tau is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Collapse and Fragmentation of Cylindrical Magnetized Clouds. II. Simulation with Nested Grid Scheme

Fragmentation process in a cylindrical magnetized cloud is studied with the nested grid method. The nested grid scheme use 15 levels of grids with different spatial resolution overlaid subsequently, which enables us to trace the evolution from the molecular cloud density $\sim 100 {\rm cm}^{-3}$ to that of the protostellar disk $\sim 10^{10} {\rm cm} ^{-3}$ or more. Fluctuation with small amplitude grows by the gravita- tional instability. It forms a disk perpendicular to the magnetic fields which runs in the direction parallel to the major axis of the cloud. Matter accrets on to the disk mainly flowing along the magnetic fields and this makes the column density increase. The radial inflow, whose velocity is slower than that perpendicular to the disk, is driven by the increase of the gravity. While the equation of state is isothermal and magnetic fields are perfectly coupled with the matter, which is realized in the density range of $ρ\lesssim 10^{10} {\rm cm}^{-3}$, never stops the contraction. The structure of the contracting disk reaches that of a singular solution as the density and the column density obey $ρ(r)\propto r^{-2}$ and $σ(r) \propto r^{-1}$, respectively. The magnetic field strength on the mid-plane is proportional to $ρ(r)^{1/2}$ and further that of the center ($B_c$) evolves as proportional to the square root of the gas density ($\propto ρ_c^{1/2}$). It is shown that isothermal clouds experience ``run-away'' collapses. The evolution after the equation of state becomes hard is also discussed.

0 citations0 reviewsNo code linkedOpen access