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

Kengo Tomida

Kengo Tomida 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-ph.EPastro-ph.HEastro-ph.IMphysics.comp-ph

Catalog footprint

What is connected

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

34 published item(s)

preprint2022arXiv

A VLA View of the Flared, Asymmetric Disk Around the Class 0 Protostar L1527 IRS

We present high resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm. We test the data against models including gap features in the intensity profile, and though we cannot rule such models out, they do not provide a statistically significant improvement in the quality of fit to the data. From these fits, we can however place constraints on allowed properties of any gaps that could be present in the true, underlying intensity profile. The physical nature of the asymmetry is difficult to associate with physical features due to the edge-on nature of the disk, but could be related to spiral arms or asymmetries seen in other imaging of more face-on disks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Interaction of a Relativistic Magnetized Collisionless Shock with a Dense Clump

The interactions between a relativistic magnetized collisionless shock and dense clumps have been expected to play a crucial role on the magnetic field amplification and cosmic-ray acceleration. We investigate this process by two-dimensional Particle-In-Cell (PIC) simulations for the first time, where the clump size is much larger than the gyroradius of downstream particles. We also perform relativistic magnetohydrodynamic (MHD) simulations for the same condition to see the kinetic effects. We find that particles escape from the shocked clump along magnetic field lines in the PIC simulations, so that the vorticity is lower than that in the MHD simulations. Moreover, in both the PIC and MHD simulations, the shocked clump quickly decelerates because of relativistic effects. Owing to the escape and the deceleration, the shocked clump cannot amplify the downstream magnetic field in relativistic collisionless shocks. This large-scale PIC simulation opens a new window to understand large-scale behaviors in collisionless plasma systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nature of supersonic turbulence and density distribution function in the multiphase interstellar medium

Supersonic flows in the interstellar medium (ISM) are believed to be a key driver of the molecular cloud formation and evolution. Among molecular clouds' properties, the ratio between the solenoidal and compressive modes of turbulence plays important roles in determining the star formation efficiency. We use numerical simulations of supersonic converging flows of the warm neutral medium (WNM) resolving the thermal instability to calculate the early phase of molecular cloud formation, and investigate the turbulence structure and the density probability distribution function (density PDF) of the multiphase ISM. We find that both the solenoidal and compressive modes have their power spectrum similar to the Kolmogorov spectrum. The solenoidal (compressive) modes account for >~80% (<~20%) of the total turbulence power. When we consider both the cold neutral medium (CNM) and the thermally unstable neutral medium (UNM) up to T <~ 400 K, the density PDF follows the log-normal distribution whose width sigma_s is well explained by the known relation from the isothermal turbulence as sigma_s = ln(1 + b^2 * M^2) (where b is the parameter representing the turbulence mode ratio and M is the turbulent Mach number). The density PDF of the CNM component alone (T <= 50 K), however, exhibits a narrower sigma_s by a factor of ~ 2. These results suggest that observational estimations of b based on the CNM density PDF requires the internal turbulence within each CNM clump but not the inter-clump relative velocity, the latter of which is instead powered by the WNM/UNM turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Osaka Feedback Model II: Modeling Supernova Feedback Based on High-Resolution Simulations

Feedback from supernovae (SNe) is an essential mechanism that self-regulates the growth of galaxies, and a better model of SN feedback is still needed in galaxy formation simulations. In the first part of this paper, using an Eulerian hydrodynamic code Athena++, we find universal scaling relations for the time evolution of momentum and radius for a superbubble, when the momentum and time are scaled by those at the shell-formation time. In the second part of this paper, we develop an SN feedback model based on the Athena++ simulation results utilizing Voronoi tessellation around each star particle, and implement it into the GADGET3-Osaka smoothed particle hydrodynamic code. Our feedback model was demonstrated to be isotropic and conservative in terms of energy and momentum. We examined the mass/energy/metal loading factors and find that our stochastic thermal feedback model produced galactic outflow that carries metals high above the galactic plane but with weak suppression of star formation. Additional mechanical feedback further suppressed star formation and brought the simulation results in better agreement with the observations of the Kennicutt--Schmidt relation, with all the results being within the uncertainties of observed data. We argue that both thermal and mechanical feedback are necessary for the SN feedback model of galaxy evolution when an individual SN bubble is unresolved.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Universal Properties of Dense Clumps in Magnetized Molecular Clouds Formed through Shock Compression of Two-phase Atomic Gases

We investigate the formation of molecular clouds from atomic gas by using three-dimensional magnetohydrodynamical simulations,including non-equilibrium chemical reactions, heating/cooling processes, and self-gravity by changing the collision speed $V_0$ and the angle $θ$ between the magnetic field and colliding flow. We found that the efficiency of the dense gas formation depends on $θ$. For small $θ$, anisotropic super-Alfvénic turbulence delays the formation of gravitationally unstable clumps. An increase in $θ$ develops shock-amplified magnetic fields along which the gas is accumulated, making prominent filamentary structures. We further investigate the statistical properties of dense clumps identified with different density thresholds. The statistical properties of the dense clumps with lower densities depend on $V_0$ and $θ$ because their properties are inherited from the global turbulence structure of the molecular clouds. By contrast, denser clumps appear to have asymptotic universal statistical properties, which do not depend on the properties of the colliding flow significantly. The internal velocity dispersions approach subsonic and plasma $β$ becomes order of unity. We develop an analytic formula of the virial parameter which reproduces the simulation results reasonably well. This property may be one of the reasons for the universality of the initial mass function of stars.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Optab: Public code for generating gas opacity tables for radiation hydrodynamics simulations

We have developed a public code, Optab, that outputs Rosseland, Planck, and two-temperature Planck mean gas opacity tables for radiation hydrodynamics simulations in astrophysics. The code is developed for modern high-performance computing, being written in Fortran 90 and using Message Passing Interface and Hierarchical Data Format, Version 5. The purpose of this work is to provide a platform on which users can generate opacity tables for their own research purposes. Therefore, the code has been designed so that a user can easily modify, change, or add opacity sources in addition to those already implemented, which include bremsstrahlung, photoionization, Rayleigh scattering, line absorption, and collision-induced absorption. In this paper, we provide details of the opacity calculations in our code and present validation tests to evaluate the performance of our code.

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

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

Formation and Evolution of Disks around Young Stellar Objects

Recent observations have suggested that circumstellar disks may commonly form around young stellar objects. Although the formation of circumstellar disks can be a natural result of the conservation of angular momentum in the parent cloud, theoretical studies instead show disk formation to be difficult from dense molecular cores magnetized to a realistic level, owing to efficient magnetic braking that transports a large fraction of the angular momentum away from the circumstellar region. We review recent progress in the formation and early evolution of disks around young stellar objects of both low-mass and high-mass, with an emphasis on mechanisms that may bridge the gap between observation and theory, including non-ideal MHD effects and asymmetric perturbations in the collapsing core (e.g., magnetic field misalignment and turbulence). We also address the associated processes of outflow launching and the formation of multiple systems, and discuss possible implications in properties of protoplanetary disks.

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

GW Ori: Interactions Between a Triple-star System and its Circumtriple Disk in Action

GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Protostellar collapse: regulation of the angular momentum and onset of an ionic precursor

Through the magnetic braking and the launching of protostellar outflows, magnetic fields play a major role in the regulation of angular momentum in star formation, which directly impacts the formation and evolution of protoplanetary disks and binary systems. The aim of this paper is to quantify those phenomena in the presence of non-ideal magnetohydrodynamics effects, namely the Ohmic and ambipola r diffusion. We perform three-dimensional simulations of protostellar collapses varying the mass of the prestellar dense core, the thermal support (the $α$ ratio) and the dust grain size-distribu tion. The mass mostly influences the magnetic braking in the pseudo-disk, while the thermal support impacts the accretion rate and hence the properties of the disk. Removing the grains smaller than 0. 1 $μ$m in the Mathis, Rumpl, Nordsieck (MRN) distribution enhances the ambipolar diffusion coefficient. Similarly to previous studies, we find that this change in the distribution reduces the magnet ic braking with an impact on the disk. The outflow is also significantly weakened. In either case, the magnetic braking largely dominates the outflow as a process to remove the angular momentum from t he disk. Finally, we report a large ionic precursor to the outflow with velocities of several km s$^{-1}$, which may be observable.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Salt, Hot Water, and Silicon Compounds Tracing Massive Twin Disks

We report results of 0.05"-resolution observations toward the O-type proto-binary system IRAS 16547-4247 with the Atacama Large Millimeter/submillimeter Array (ALMA). We present dynamical and chemical structures of the circumbinary disk, circumstellar disks, outflows and jets, illustrated by multi-wavelength continuum and various molecular lines. In particular, we detect sodium chloride, silicon compounds, and vibrationally-excited water lines as probes of the individual protostellar disks at a scale of 100 au. These are complementary to typical hot-core molecules tracing the circumbinary structures on a 1000-au scale. The H2O line tracing inner-disks has an upper-state energy of Eu/k>3000K, indicating a high temperature of the disks. On the other hand, despite the detected transitions of NaCl, SiO, and SiS not necessarily having high upper-state energies, they are enhanced only in the vicinity of the protostars. We interpret that these molecules are the products of dust destruction, which only happens in the inner disks. This is the second detection of alkali metal halide in protostellar systems after the case of the disk of Orion Source I, and also one of few massive protostellar disks associated with high-energy transition water and silicon compounds. These new results suggest these "hot-disk" lines may be common in innermost disks around massive protostars, and have great potential for future research of massive star formation. We also tentatively find that the twin disks are counter-rotating, which might give a hint of the origin of the massive proto-binary system IRAS 16547-4247.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers

The design and implementation of a new framework for adaptive mesh refinement (AMR) calculations is described. It is intended primarily for applications in astrophysical fluid dynamics, but its flexible and modular design enables its use for a wide variety of physics. The framework works with both uniform and nonuniform grids in Cartesian and curvilinear coordinate systems. It adopts a dynamic execution model based on a simple design called a "task list" that improves parallel performance by overlapping communication and computation, simplifies the inclusion of a diverse range of physics, and even enables multiphysics models involving different physics in different regions of the calculation. We describe physics modules implemented in this framework for both non-relativistic and relativistic magnetohydrodynamics (MHD). These modules adopt mature and robust algorithms originally developed for the Athena MHD code and incorporate new extensions: support for curvilinear coordinates, higher-order time integrators, more realistic physics such as a general equation of state, and diffusion terms that can be integrated with super-time-stepping algorithms. The modules show excellent performance and scaling, with well over 80% parallel efficiency on over half a million threads. The source code has been made publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Transition region from turbulent to dead zone in protoplanetary disks: local shearing box simulations

The dynamical evolution of protoplanetary disks is of key interest for building a comprehensive theory of planet formation and to explain the observational properties of these objects. Using the magnetohydrodynamics code Athena++, with an isothermal shearing box setup, we study the boundary between the active and dead zone, where the accretion rate changes and mass can accumulate. We quantify how the turbulence level is affected by the presence of a non uniform ohmic resistivity in the radial-x direction that leads to a region of inhibited turbulence (or dead zone). Comparing the turbulent activity to that of ideal simulations, the turbulence inhibited area shows density fluctuations and magnetic activity at its boundaries, driven by energy injection from the active (ideal) zone boundaries. We find magnetic dissipation to be significantly stronger in the ideal regions, and the turbulence penetration through the boundary of the dead zone is determined by the value of the resistivity itself, through the ohmic dissipation process, though the thickness of the transition does not play a significant role in changing the dissipation. We investigate the 1D spectra along the shearing direction: magnetic spectra appear flat at large scales both in ideal as well as resistive simulations, though a Kolmogorov scaling over more than one decade persists in the dead zone, suggesting the turbulent cascade is determined by the hydrodynamics of the system: MRI dynamo action is inhibited where sufficiently high resistivity is present.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The detection of dust gap-ring structure in the outer region of the CR Cha protoplanetary disk

We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Binary Stellar Mergers with Marginally-Bound Ejecta: Excretion Disks, Inflated Envelopes, Outflows, and their Luminous Transients

We study mass loss from the outer Lagrange point (L2) in binary stellar mergers and their luminous transients by means of radiative hydrodynamical simulations. Previously, we showed that for binary mass ratios 0.06 < q < 0.8, synchronous L2 mass loss results in a radiatively inefficient, dust-forming unbound equatorial outflow. A similar outflow exists irrespective of q if the ratio of the sound speed to the orbital speed at the injection point is sufficiently large, ε= c_T/v_orb > 0.15. By contrast, for cold L2 mass-loss (ε< 0.15) from binaries with q < 0.06 or q > 0.8, the equatorial outflow instead remains marginally-bound and falls back to the binary over tens to hundreds of binary orbits, where it experiences additional tidal torqueing and shocking. As the bound gas becomes virialized with the binary, the luminosity of the system increases slowly at approximately constant photosphere radius, causing the temperature to rise. Subsequent evolution depends on the efficiency of radiative cooling. If the bound atmosphere is able to cool efficiently, as quantified by radiative diffusion time being shorter than the advection time (t_diff/t_adv < 1), then the virialized gas collapses to an excretion disk, while for t_diff/t_adv > 1 an isotropic wind is formed. Between these two extremes, an inflated envelope transports the heat generated near the binary to the surface by meridional flows. In all cases, the radiated luminosity reaches a fraction ~0.01 to 0.1 of Mdot v_orb^2/2, where Mdot is the mass outflow rate. We discuss the implications of our results for transients in the luminosity gap between classical novae and supernovae, such as V1309 Sco and V838 Mon.

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

The Origin of the Rotation Profiles in Star Forming Clouds

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

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

Cool and Luminous Transients from Mass-Losing Binary Stars

We study transients produced by equatorial disk-like outflows from catastrophically mass-losing binary stars with an asymptotic velocity and energy deposition rate near the inner edge which are proportional to the binary escape velocity v_esc. As a test case, we present the first smoothed-particle radiation-hydrodynamics calculations of the mass loss from the outer Lagrange point with realistic equation of state and opacities. The resulting spiral stream becomes unbound for binary mass ratios 0.06 < q < 0.8. For synchronous binaries with non-degenerate components, the spiral-stream arms merge at a radius of ~10a, where a is the binary semi-major axis, and the accompanying shock thermalizes about 10% of the kinetic power of the outflow. The mass-losing binary outflows produce luminosities reaching up to ~10^6 L_Sun and effective temperatures spanning 500 < T_eff < 6000 K, which is compatible with many of the class of recently-discovered red transients such as V838 Mon and V1309 Sco. Dust readily forms in the outflow, potentially in a catastrophic global cooling transition. The appearance of the transient is viewing angle-dependent due to vastly different optical depths parallel and perpendicular to the binary plane. We predict a correlation between the peak luminosity and the outflow velocity, which is roughly obeyed by the known red transients. Outflows from mass-losing binaries can produce luminous (10^5 L_Sun) and cool (T_eff < 1500 K) transients lasting a year or longer, as has potentially been detected by Spitzer surveys of nearby galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

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

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

On the role of the H2 ortho:para ratio in gravitational collapse during star formation

Hydrogen molecules (H2) come in two forms in the interstellar medium, ortho- and para-hydrogen, corresponding to the two different spin configurations of the two hydrogen atoms. The relative abundances of the two flavours in the interstellar medium are still very uncertain, and this abundance ratio has a significant impact on the thermal properties of the gas. In the context of star formation, theoretical studies have recently adopted two different strategies when considering the ortho:para ratio (OPR) of H2 molecules; the first considers the OPR to be frozen at 3:1 while the second assumes that the species are in thermal equilibrium. As the OPR potentially affects the protostellar cores which form as a result of the gravitational collapse of a dense molecular cloud, the aim of this paper is to quantify precisely what role the choice of OPR plays in the properties and evolution of the cores. We used two different ideal gas equations of state for a hydrogen and helium mix in a radiation hydrodynamics code to simulate the collapse of a dense cloud and the formation of the first and second Larson cores; the first equation of state uses a fixed OPR of 3:1 while the second assumes thermal equilibrium. Simulations using an equilibrium ratio collapse faster at early times and show noticeable oscillations around hydrostatic equilibrium, to the point where the core expands for a short time right after its formation before resuming its contraction. In the case of a fixed 3:1 OPR, the core's evolution is a lot smoother. The OPR was however found to have little impact on the size, mass and radius of the two Larson cores. We conclude that if one is solely interested in the final properties of the cores when they are formed, it does not matter which OPR is used. On the other hand, if one's focus lies primarily in the evolution of the first core, the choice of OPR becomes important.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Low-Metallicity Environments

Among many physical processes involved in star formation, radiation transfer is one of the key processes since it dominantly controls the thermodynamics. Because metallicities control opacities, they are one of the important environmental parameters which affect star formation processes. In this work, I investigate protostellar collapse in solar-metallicity and low-metallicity ($Z = 0.1 Z_\odot$) environments using 3D radiation hydrodynamic and magnetohydrodynamic simulations. Because radiation cooling is more effective in the low-metallicity environments, first cores are colder and have lower entropies. As a result, first cores are smaller, less massive and have shorter lifetimes in the low-metallicity clouds. Therefore, first cores would be less likely to be found in low-metallicity star forming clouds. This also implies that first cores tend to be more gravitationally unstable and susceptible to fragmentation. The evolution and structure of protostellar cores formed after the second collapse weakly depend on metallicities in the spherical and magnetized models despite the large difference in the metallicities. Because this is due to the change of the heat capacity by dissociation and ionization of hydrogen, it is a general consequence of the second collapse as long as the effects of radiation cooling are not very large during the second collapse. On the other hand, the effects of different metallicities are more significant in the rotating models without magnetic fields, because they evolve slower than other models and therefore more affected by radiation cooling.

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

Explicit-Implicit Scheme for Relativistic Radiation Hydrodynamics

We propose an explicit-implicit scheme for numerically solving Special Relativistic Radiation Hydrodynamic (RRHD) equations, which ensures a conservation of total energy and momentum (matter and radiation). In our scheme, 0th and 1st moment equations of the radiation transfer equation are numerically solved without employing a flux-limited diffusion (FLD) approximation. For an hyperbolic term, of which the time scale is the light crossing time when the flow velocity is comparable to the speed of light, is explicitly solved using an approximate Riemann solver. Source terms describing an exchange of energy and momentum between the matter and the radiation via the gas-radiation interaction are implicitly integrated using an iteration method. The implicit scheme allows us to relax the Courant-Friedrichs-Lewy condition in optically thick media, where heating/cooling and scattering timescales could be much shorter than the dynamical timescale. We show that our numerical code can pass test problems of one- and two-dimensional radiation energy transport, and one-dimensional radiation hydrodynamics. Our newly developed scheme could be useful for a number of relativistic astrophysical problems. We also discuss how to extend our explicit-implicit scheme to the relativistic radiation magnetohydrodynamics.

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
preprint2012arXiv

Spatially Resolving Substructures within the Massive Envelope around an Intermediate-mass Protostar: MMS 6/OMC-3

With the Submillimeter Array, the brightest (sub)millimeter continuum source in the OMC-2/3 region, MMS 6, has been observed in the 850 um continuum emission with approximately 10 times better angular resolution than previous studies (~0.3"; ~120 AU at Orion). The deconvolved size, the mass, and the column density of MMS 6-main are estimated to be 0.32"x0.29" (132 AUx120 AU), 0.29 Mo, and 2.1x10^{25} cm^{-2}, respectively. The estimated extremely high mean number density, 1.5x10^{10} cm^{-3}, suggests that MMS 6-main is likely optically thick at 850 um. We compare our observational data with three theoretical core models: prestellar core, protostellar core + disk-like structure, and first adiabatic core. These comparisons clearly show that the observational data cannot be modeled as a simple prestellar core with a gas temperature of 20 K. A self-luminous source is necessary to explain the observed flux density in the (sub)millimeter wavelengths. Our recent detection of a very compact and energetic outflow in the CO (3-2) and HCN (4-3) lines, supports the presence of a protostar. We suggest that MMS 6 is one of the first cases of an intermediate mass protostellar core at an extremely young stage. In addition to the MMS 6-main peak, we have also spatially resolved a number of spiky structures and sub-clumps, distributed over the central 1000 AU. The masses of these sub-clumps are estimated to be 0.066-0.073 Mo, which are on the order of brown dwarf masses. Higher angular resolution and higher sensitivity observations with ALMA and EVLA will reveal the origin and nature of these structures such as whether they are originated from fragmentations, spiral arms, or inhomogeneity within the disk-like structures/envelope.

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