Researcher profile

Kengo Tomida

Kengo Tomida contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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.

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&#39; 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.

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.

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.

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.

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 $&#34;$starless$&#34;$ 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 $&#34;$starless$&#34;$ cores in Taurus has blueshifted/redshifted wings in the $^{12}$CO (2-1) line, indicating that there is deeply embedded object driving molecular outflow. The observed velocities and sizes of the possible outflow lobes are 2-4 km s$^{-1}$, and $\sim$2 $\times$10$^3$ au, respectively, and the dynamical time is calculated to be $\sim$10$^3$ yr. In addition to this, the core is one of the strongest N$_2$D$^{+}$ (3-2) emitters in our sample. All of the observed signatures do not conflict with any of the theoretical predictions about the first hydrostatic core so far, and thus MC35-mm is unique as the only first-core candidate in the Taurus molecular cloud.

preprint2020arXiv

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

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

preprint2020arXiv

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.

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

preprint2020arXiv

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.

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.

preprint2020arXiv

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

We report results of 0.05&#34;-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 &#34;hot-disk&#34; 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.

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 &#34;task list&#34; 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.

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.

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.