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

Tsuyoshi Inoue

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

ResearcherUnclaimed source record
astro-ph.HEastro-ph.GAastro-ph.SRphysics.plasm-phphysics.flu-dynastro-phastro-ph.IMphysics.comp-ph

Catalog footprint

What is connected

32works
8topics
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

32 published item(s)

preprint2022arXiv

An ALMA study of the massive molecular clump N159W-North in the Large Magellanic Cloud: A possible gas flow penetrating one of the most massive protocluster systems in the Local Group

Massive dense clumps in the Large Magellanic Cloud can be an important laboratory to explore the formation of populous clusters. We report multiscale ALMA observations of the N159W-North clump, which is the most CO-intense region in the galaxy. High-resolution CO isotope and 1.3 mm continuum observations with an angular resolution of $\sim$0."25($\sim$0.07 pc) revealed more than five protostellar sources with CO outflows within the main ridge clump. One of the thermal continuum sources, MMS-2, shows especially massive/dense nature whose total H$_2$ mass and peak column density are $\sim$10$^{4}$ $M_{\odot}$ and $\sim$10$^{24}$ cm$^{-2}$, respectively, and harbors massive ($\sim$100 $M_{\odot}$) starless core candidates identified as its internal substructures. The main ridge containing this source can be categorized as one of the most massive protocluster systems in the Local Group. The CO high-resolution observations found several distinct filamentary clouds extending southward from the star-forming spots. The CO (1-0) data set with a larger field of view reveals a conical-shaped, $\sim$30 pc long complex extending toward the northern direction. These features indicate that a large-scale gas compression event may have produced the massive star-forming complex. Based on the striking similarity between the N159W-North complex and the previously reported other two high-mass star-forming clouds in the nearby regions, we propose a $"$teardrops inflow model$"$ that explains the synchronized, extreme star formation across $>$50 pc, including one of the most massive protocluster clumps in the Local Group.

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

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

On the Potential of Faraday Tomography to Identify Shock Structures in Supernova Remnants

Knowledge about the magnetic fields in supernova remnants (SNRs) is of paramount importance for constraining Galactic cosmic ray acceleration models. It could also indirectly provide information on the interstellar magnetic fields. In this paper, we predict the Faraday dispersion functions (FDFs) of SNRs for the first time. For this study, we use the results of three dimensional (3D) ideal magnetohydrodynamic (MHD) simulations of SNRs expanding into a weak, regular magnetic field. We present the intrinsic FDFs of the shocked region of SNRs for different viewing angles. We find that the FDFs are generally Faraday complex, which implies that conventional rotation measure study is not sufficient to obtain the information on the magnetic fields in the shocked region and Faraday tomography is necessary. We also show that the FDF allows to derive the physical-depth distribution of polarization intensity when the line of sight is parallel to the initial magnetic field orientation. Furthermore, we demonstrate that the location of contact discontinuity can be identified from the radial profile of the width of the FDF with the accuracy of 0.1-0.2 pc.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

X-ray Line Diagnostics of Ion Temperature at Cosmic-Ray Accelerating Collisionless Shocks

A novel collisionless shock jump condition is suggested by modeling the entropy production at the shock transition region. We also calculate downstream developments of the atomic ionization balance and the ion temperature relaxation in supernova remnants (SNRs). The injection process and subsequent acceleration of cosmic-rays (CRs) in the SNR shocks are closely related to the formation process of the collisionless shocks. The formation of the shock is caused by wave-particle interactions. Since the wave-particle interactions result in energy exchanges between electromagnetic fields and charged particles, the randomization of particles associated with the shock transition may occur with the rate given by the scalar product of the electric field and current. We find that order-of-magnitude estimates of the randomization with reasonable strength of the electromagnetic fields in the SNR constrain the amount of the CR nuclei and ion temperatures. The constrained amount of the CR nuclei can be sufficient to explain the Galactic CRs. The ion temperature becomes significantly lower than in the case of no CRs. To distinguish the case without CRs, we perform synthetic observations of atomic line emissions from the downstream region of the SNR RCW~86. Future observations by {\it XRISM} and {\it Athena} can distinguish whether the SNR shock accelerates the CRs or not from the ion temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

FOREST unbiased Galactic plane imaging survey with the Nobeyama 45 m telescope (FUGIN). VII. molecular fraction of HI clouds

In this study, we analyze molecular gas formation in neutral atomic hydrogen (HI) clouds using the latest CO data obtained from the four-beam receiver system on a 45-m telescope (FOREST) unbiased Galactic plane imaging survey with the Nobeyama 45-m telescope (FUGIN) and HI data taken from the Very Large Array (VLA) Galactic plane survey (VGPS). We applied a dendrogram algorithm to the HI data cube to identify HI clouds, and we calculated the HI mass and molecular gas mass by summing the CO line intensity within each HI cloud. On the basis of the results, we created a catalog of 5,737 identified HI clouds with local standard of rest (LSR) velocity of $V_{\rm LSR}\le -20$ km s$^{-1}$ in Galactic longitude and latitude ranges of $20^\circ \le l \le 50^\circ$ and $-1^\circ \le b \le 1^\circ$, respectively. We found that most of the HI clouds are distributed within a Galactocentric distance of 16 kpc, most of which are in the Cold Neutral Medium (CNM) phase. In addition, we determined that the high-mass end of the mass HI function is well fitted with the power-law function with an index of 2.3. Although two sequences of self-gravitating and diffuse clouds are expected to appear in the M$_{\rm tot}$-M$_{{\rm H}_2}$ diagram according to previous works based on a plane-parallel model, the observational data show only a single sequence with large scattering within these two sequences. This implies that most of the clouds are mixtures of these two types of clouds. Moreover, we suggest the following scenario of molecular gas formation: An HI-dominant cloud evolved with increasing H$_2$ mass along a path of $M_{{\rm H}_2} \propto M_{\rm tot}^2$ by collecting diffuse gas before reaching and moving along the curves of the two sequences.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Formation of Massive Star Clusters by Fast HI Gas Collision

Young massive clusters (YMCs) are dense aggregates of young stars, which are essential to galaxy evolution, owing to their ultraviolet radiation, stellar winds, and supernovae. The typical mass and radius of YMCs are M~10^4 M_sun and R~1 pc, respectively, indicating that many stars are located in a small region. The formation of YMC precursor clouds may be difficult because a very compact massive cloud should be formed before stellar feedback blows off the cloud. Recent observational studies suggest that YMCs can be formed as a consequence of the fast HI gas collision with a velocity of ~100 km s^-1, which is the typical velocity of the galaxy-galaxy interaction. In this study, we examine whether the fast HI gas collision triggers YMC formation using three-dimensional magnetohydrodynamics simulations, which include the effects of self-gravity, radiative cooling/heating, and chemistry. We demonstrate that massive gravitationally bound gas clumps with M >10^4 M_sun and L~4 pc are formed in the shock compressed region induced by the fast HI gas collision, which massive gas clumps can evolve into YMCs. Our results show that the YMC precursors are formed by the global gravitational collapse of molecular clouds, and YMCs can be formed even in low-metal environments, such as the Magellanic Clouds. Additionally, the very massive YMC precursor cloud, with M >10^5 M_sun, can be created when we consider the fast collision of HI clouds, which may explain the origin of the very massive stellar cluster R136 system in the Large Magellanic Cloud.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rapid and efficient mass collection by a supersonic cloud-cloud collision as a major mechanism of high-mass star formation

A supersonic cloud-cloud collision produces a shock-compressed layer which leads to formation of high-mass stars via gravitational instability. We carried out a detailed analysis of the layer by using the numerical simulations of magneto-hydrodynamics which deal with colliding molecular flows at a relative velocity of 20 km s$^{-1}$ (Inoue & Fukui 2013). Maximum density in the layer increases from 1000 cm$^{-3}$ to more than $10^{5}$ cm$^{-3}$ within 0.3 Myrs by compression, and the turbulence and the magnetic field in the layer are amplified by a factor of $\sim 5$, increasing the mass accretion rate by two orders of magnitude to more than $10^{-4}$ $M_{\odot}$ yr$^{-1}$. The layer becomes highly filamentary due to gas flows along the magnetic field lines, and dense cores are formed in the filaments. The massive dense cores have size and mass of 0.03 -- 0.8 pc and 8 -- 50 $M_{\odot}$ and they are usually gravitationally unstable. The mass function of the dense cores is significantly top-heavy as compared with the universal IMF, indicating that the cloud-cloud collision triggers preferentially the formation of O and early B stars. We argue that the cloud-cloud collision is a versatile mechanism which creates a variety of stellar clusters from a single O star like RCW120 and M20 to tens of O stars of a super star cluster like RCW38 and a mini-starburst W43. The core mass function predicted by the present model is consistent with the massive dense cores obtained by recent ALMA observations in RCW38 (Torii et al. 2019) and W43 (Motte et al. 2018) considering the increasing evidence for collision-triggered high-mass star formation, we argue that cloud-cloud collision is a major mechanism of high mass star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation of HI Clouds in Shock-compressed Interstellar Medium: Physical Origin of Angular Correlation Between Filamentary Structure and Magnetic Field

Recent observations of neutral Galactic interstellar medium showed that filamentary structures of HI clouds are aligned with the interstellar magnetic field. Many interesting applications are proposed based on the alignment such as measurement of magnetic field strength through the Chandrasekhar-Fermi method and removal of polarized foreground dust emissions for the detection of inflationary polarized emission in the cosmic microwave background radiation. However, the physical origin of the alignment remains to be explained. To understand the alignment mechanism, we examine formation of HI clouds triggered by shock compression of diffuse warm neutral medium using three-dimensional magnetohydrodynamic simulations with the effects of optically thin cooling and heating. We show that the shock-compressed diffuse interstellar medium of density n~1 cm^-3 evolves into HI clouds with typical density n~50 cm^-3 via thermal instability driven by cooling, which is consistent with previous studies. We apply a machine vision transformation developed by Clark et al.(2014) to the resulting column density structures obtained by the simulations in order to measure angle correlation between filamentary structures of HI clouds and magnetic field. We find that the orientation of HI filaments depends on the environmental turbulent velocity field, particularly on the strength of shear strain in the direction of the magnetic field, which is controlled by the angle between the shock propagation direction and upstream magnetic field. When the strain along the magnetic field is weak, filamentary components of HI clouds basically lie perpendicular to the magnetic field. However, the filaments have come to align with the magnetic field, if we enhance the turbulent strain along the magnetic field or if we set turbulence in the preshock medium.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Diffusion of Cosmic Rays in a Multiphase Interstellar Medium Swept-Up by a Supernova Remnant Blast Wave

Supernova remnants (SNRs) are one of the most energetic astrophysical events and are thought to be the dominant source of Galactic cosmic rays (CRs). A recent report on observations from the Fermi satellite has shown a signature of pion decay in the gamma-ray spectra of SNRs. This provides strong evidence that high-energy protons are accelerated in SNRs. The actual gamma-ray emission from pion decay should depend on the diffusion of CRs in the interstellar medium. In order to quantitatively analyse the diffusion of high-energy CRs from acceleration sites, we have performed test particle numerical simulations of CR protons using a three-dimensional magnetohydrodynamics (MHD) simulation of an interstellar medium swept-up by a blast wave. We analyse the diffusion of CRs at a length scale of order a few pc in our simulated SNR, and find the diffusion of CRs is precisely described by a Bohm diffusion, which is required for efficient acceleration at least for particles with energies above 30 TeV for a realistic interstellar medium. Although we find the possibility of a superdiffusive process (travel distance proportional to t^0.75) in our simulations, its effect on CR diffusion at the length scale of the turbulence in the SNR is limited.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-mass star formation triggered by collision between CO filaments in N159 West in the Large Magellanic Cloud

We have carried out 13CO(J=2-1) observations of the active star-forming region N159 West in the LMC with ALMA. We have found that the CO distribution at a sub-pc scale is highly elongated with a small width. These elongated clouds called "filaments" show straight or curved distributions with a typical width of 0.5-1.0pc and a length of 5-10pc. All the known infrared YSOs are located toward the filaments. We have found broad CO wings of two molecular outflows toward young high-mass stars in N159W-N and N159W-S, whose dynamical timescale is ~10^4 yrs. This is the first discovery of protostellar outflow in external galaxies. For N159W-S which is located toward an intersection of two filaments we set up a hypothesis that the two filaments collided with each other ~10^5 yrs ago and triggered formation of the high-mass star having ~37 Mo. The colliding clouds show significant enhancement in linewidth in the intersection, suggesting excitation of turbulence in the shocked interface layer between them as is consistent with the magneto-hydro-dynamical numerical simulations (Inoue & Fukui 2013). This turbulence increases the mass accretion rate to ~6x10^-4 Mo yr^-1, which is required to overcome the stellar feedback to form the high-mass star.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On cosmic-ray production efficiency at supernova remnant shocks propagating into realistic diffuse interstellar medium

Using three-dimensional magnetohydrodynamics simulations, we show that the efficiency of cosmic-ray (CR) production at supernova remnants (SNRs) is over-predicted if it could be estimated based on proper motion measurements of H$α$ filaments in combination with shock-jump conditions. Density fluctuations of upstream medium make shock waves rippled and oblique almost everywhere. The kinetic energy of the shock wave is transferred into that of downstream turbulence as well as thermal energy which is related to the shock velocity component normal to the shock surface. Our synthetic observation shows that the CR acceleration efficiency as estimated from a lower downstream plasma temperature, is overestimated by 10-40%, because rippled shock does not immediately dissipate all upstream kinetic energy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Formation and Destruction of Molecular Clouds and Galactic Star Formation

We describe an overall picture of galactic-scale star formation. Recent high-resolution magneto-hydrodynamical simulations of two-fluid dynamics with cooling/heating and thermal conduction have shown that the formation of molecular clouds requires multiple episodes of supersonic compression. This finding enables us to create a scenario in which molecular clouds form in interacting shells or bubbles on a galactic scale. First we estimate the ensemble-averaged growth rate of molecular clouds over a timescale larger than a million years. Next we perform radiation hydrodynamics simulations to evaluate the destruction rate of magnetized molecular clouds by the stellar FUV radiation. We also investigate the resultant star formation efficiency within a cloud which amounts to a low value (a few percent) if we adopt the power-law exponent -2.5 for the mass distribution of stars in the cloud. We finally describe the time evolution of the mass function of molecular clouds over a long timescale (>1Myr) and discuss the steady state exponent of the power-law slope in various environments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Turbulent Reconnection in Relativistic Plasmas And Effects of Compressibility

We report turbulence effects on magnetic reconnection in relativistic plasmas using 3-dimensional relativistic resistive magnetohydrodynamics simulations. We found reconnection rate became independent of the plasma resistivity due to turbulence effects similarly to non-relativistic cases. We also found compressible turbulence effects modified the turbulent reconnection rate predicted in non-relativistic incompressible plasmas; The reconnection rate saturates and even decays as the injected velocity approaches to the Alfvén velocity. Our results indicate the compressibility cannot be neglected when compressible component becomes about half of incompressible mode occurring when the Alfvén Mach number reaches about $0.3$. The obtained maximum reconnection rate is around $0.05$ to $0.1$, which will be able to reach around $0.1$ to $0.2$ if injection scales are comparable to the sheet length.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evaporation of grain-surface species by shock waves in proto-planetary disk

Recent ALMA (Atacama Large Millimeter/submillimeter Array) observations of young protostellar objects detected warm SO emission, which could be associated with a forming protostellar disk. In order to investigate if such warm gas can be produced by accretion shock onto the forming disk, we calculate the sputtering and thermal desorption of various grain surface species in one dimensional shock waves. We find that thermal desorption is much more efficient than the sputtering in the post-shock region. While H$_{2}$O can be thermally desorbed, if the accretion velocity is larger than 8 km s$^{-1}$ with the pre-shock gas number density of 10$^{9}$ cm$^{-3}$, SO is desorbed, if the accretion velocity $\gtrsim$ 2 km s$^{-1}$ and $\gtrsim$ 4km s$^{-1}$, with the pre-shock density of 10$^{9}$ cm$^{-3}$ and 10$^{8}$ cm$^{-3}$, respectively. We also find that the column density of hydrogen nuclei in warm post-shock gas is $N_{\rm warm} \sim 10^{21}$ cm$^{-2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Thermal instability and multi-phase interstellar medium in the first galaxies

We examine the linear stability and nonlinear growth of the thermal instability in isobarically contracting gas with various metallicities and FUV field strengths. When the H2 cooling is suppressed by FUV fields (G_0>10^-3) or the metallicity is high enough (Z/Zs>10^-3), the interstellar medium is thermally unstable in the temperature range 100-7000 K owing to the cooling by CII and OI fine-structure lines. In this case, a bi-phasic medium with a bimodal density probability distribution is formed as a consequence of the thermal instability. The characteristic scales of the thermal instability become smaller with increasing metallicity. Comparisons of the nonlinear simulations with different resolution indicates that the maximum scale of the thermal instability should be resolved with more than 60 cells to follow runaway cooling driven by the thermal instability. Under sufficiently weak FUV fields and with low metallcity, the density range of the thermal instability shrinks owing to dominance of H2 cooling. As the FUV intensity is reduced, bi-phasic structure becomes less remarkable and disappears eventually. Our basic results suggest that, in early galaxies, i) the thermal instability has little effect for the medium with Z/Zs<10^-4, ii) fragmentation by the thermal instability could determine mass spectrum of star clusters for 10^-4<Z/Zs<0.04, and iii) thermally bistable turbulent interstellar medium like our galaxy becomes ubiquitous for Z/Zs>0.04, although the threshold metallicity depends on the conditions such as thermal pressure, FUV strength and redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Critical Magnetic Field Strength for Suppression of the Richtmyer-Meshkov Instability in Plasmas

The critical strength of a magnetic field required for the suppression of the Richtmyer-Meshkov instability (RMI) is investigated numerically by using a two-dimensional single-mode analysis. For the cases of MHD parallel shocks, the RMI can be stabilized as a result of the extraction of vorticity from the interface. A useful formula describing a critical condition for MHD RMI has been introduced, and which is successfully confirmed by the direct numerical simulations. The critical field strength is found to be largely depending on the Mach number of the incident shock. If the shock is strong enough, even low-$β$ plasmas can be subject to the growth of the RMI.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Formation of Massive Molecular Cloud Cores by Cloud-cloud Collision

Recent observations of molecular clouds around rich massive star clusters including NGC3603, Westerlund 2, and M20 revealed that the formation of massive stars could be triggered by a cloud-cloud collision. By using three-dimensional, isothermal, magnetohydrodynamics simulations with the effect of self-gravity, we demonstrate that massive, gravitationally unstable, molecular cloud cores are formed behind the strong shock waves induced by the cloud-cloud collision. We find that the massive molecular cloud cores have large effective Jeans mass owing to the enhancement of the magnetic field strength by shock compression and turbulence in the compressed layer. Our results predict that massive molecular cloud cores formed by the cloud-cloud collision are filamentary and threaded by magnetic fields perpendicular to the filament.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Origin of Radially Aligned Magnetic Fields in Young Supernova Remnants

It has been suggested by radio observations of polarized synchrotron emissions that downstream magnetic field in some young supernova remnants are oriented radially. We study magnetic field distribution of turbulent supernova remnant driven by the Richtmyer-Meshkov instability -- in other words, the effect of rippled shock -- by using three-dimensional magnetohydrodynamics simulations. We find that the induced turbulence has radially biased anisotropic velocity dispersion that leads to a selective amplification of the radial component of the magnetic field. The Richtmyer-Meshkov instability is induced by the interaction between the shock and upstream density fluctuations. Future high-resolution polarization observations can distinguish the following candidates responsible for the upstream density fluctuations: (i) inhomogeneity caused by the cascade of large-scale turbulence in the ISM so-called the big-power-law-in-the-sky, (ii) structures generated by the Drury instability in the cosmic-ray modified shock, and (iii) fluctuations induced by the non-linear feedback of the cosmic-ray streaming instability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermal Instability behind a Shock Wave in HI and Molecular Clouds

We performed one-dimensional hydrodynamic simulations with detailed cooling, heating and chemical processes to examine the thermal stability of shocked gas in cold neutral medium (CNM) and molecular clouds. We find that both CNM and molecular clouds can be thermally unstable in the cooling layer behind the shock wave. The characteristic wavelength of the thermal instability ranges from 10^-5 pc to 0.1 pc in the CNM, and from 10^-7 pc to 0.1 pc in the molecular clouds. This coincides with the size of observed tiny scale structures in the CNM and molecular clouds, indicating that the thermal instability in the shocked gas could be a formation mechanism of these tiny structures in the interstellar medium. We have also calculated the e-folding number of the thermal instability to estimate the amplification of the density fluctuation in the shocked gas. Density perturbations in the CNM grow by a factor of exp(5)~150, whereas the perturbations in the molecular clouds grow only by a factor of a few behind a high Mach number shock. The amplification factor is larger at lower densities and higher velocities. Formation of very small-scale structures by thermal instability in shocked gas is more effective in lower densities.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Enhanced Dissipation Rate of Magnetic Field in Striped Pulsar Winds by the Effect of Turbulence

In this letter we report on turbulent acceleration of the dissipation of magnetic field in the postshock re- gion of a Poynting flux-dominated flow, such as the Crab pulsar wind nebula. We have performed two- dimensional resistive relativistic magnetohydrodynamics simulations of subsonic turbulence driven by the Richtmyer-Meshkov instability at the shock fronts of the Poynting flux-dominated flows in pulsar winds. We find that turbulence stretches current sheets which substantially enhances the dissipation of magnetic field, and that most of the initial magnetic field energy is dissipated within a few eddy-turnover times. We also develop a simple analytical model for turbulent dissipation of magnetic field that agrees well with our simulations. The analytical model indicates that the dissipation rate does not depend on resistivity even in the small resistivity limit. Our findings can possibly alleviate the σ-problem in the Crab pulsar wind nebulae.

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

Formation of Turbulent and Magnetized Molecular Clouds via Accretion Flows of HI Clouds

Using 3D MHD simulation with the effects of radiative cooling/heating, chemical reactions, and thermal conduction, we investigate the formation of molecular cloud in the ISM. We consider the formation of molecular cloud by accretion of the HI clouds as suggested in recent observations. The simulation shows that the initial HI medium is compressed and piled up behind the shock waves induced by the accretion flows. Since the initial medium is highly inhomogeneous as a consequence of the thermal instability, the formed molecular cloud becomes very turbulent owing to the development of the Richtmyer-Meshkov instability. The structure of the post shock region is composed of dense cold gas (T<100 K) and diffuse warm gas (T>1,000 K), which are spatially well mixed owing to the turbulence. Because the energy source of the turbulence is the accretion flows, the turbulence is highly anisotropic biased toward the direction of accretion flows. The kinetic energy of the turbulence dominates the thermal, magnetic, and gravitational energies in the total 10 Myr evolution. However, the kinetic energy measured by using the CO-fraction-weighted density is comparable to the other energies. This suggests that the true kinetic energy of turbulence in molecular cloud as a hole can be much larger than the kinetic energy of turbulence estimated by line-width of molecular emissions. The clumps in the molecular cloud show statistically homogeneous evolution as follows: The typical plasma beta of the clumps is roughly constant <beta>~ 0.4. The size-velocity dispersion relation show dv ~ 1.5 km s^{-1} (l/1 pc)^{0.5}, irrespective of the density. The clumps evolve toward magnetically supercritical, gravitationally unstable cores. The clumps seem to evolve into cores that satisfy the condition for fragmentation into binary. These statistical properties may provide the initial condition of star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Magnetic Field Amplification Associated with the Richtmyer-Meshkov Instability

The amplification of a magnetic field due to the Richtmyer-Meshkov instability (RMI) is investigated by two-dimensional MHD simulations. Single-mode analysis is adopted to reveal definite relation between the nonlinear evolution of RMI and the field enhancement. It is found that an ambient magnetic field is stretched by fluid motions associated with the RMI, and the strength is amplified significantly by more than two orders of magnitude. The saturation level of the field is determined by a balance between the amplified magnetic pressure and the thermal pressure after shock passage. This effective amplification can be achieved in a wide range of the conditions for the RMI such as the Mach number of an incident shock and the density ratio at a contact discontinuity. The results suggest that the RMI could be a robust mechanism of the amplification of interstellar magnetic fields and cause the origin of localized strong fields observed at the shock of supernova remnants.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Richtmyer-Meshkov Type Instability of a Current Sheet in a Relativistically Magnetized Plasma

Linear stability of a current sheet that is subject to an impulsive acceleration due to a shock passage is studied with the effect of guide magnetic field. We find that the current sheet embedded in relativistically magnetized plasma always shows a Richtmyer-Meshkov type instability, while it depends on the density structure in the Newtonian limit. The growth of the instability is expected to generate turbulence around the current sheet that can induce so-called turbulent reconnection whose rate is essentially free from plasma resistivity. Thus, the instability can be applied as a triggering mechanism of rapid magnetic energy release in variety of high-energy astrophysical phenomena such as pulsar wind nebulae, gamma-ray bursts, and active galactic nuclei, where the shock wave is supposed to play a crucial role.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A New numerical scheme for resistive relativistic MHD using method of characteristics

We present a new numerical method of special relativistic resistive magnetohydrodynamics with scalar resistivity that can treat a range of phenomena, from nonrelativistic to relativistic (shock, contact discontinuity, and Alfvén wave). The present scheme calculates the numerical flux of fluid by using an approximate Riemann solver, and electromagnetic field by using the method of characteristics. Since this scheme uses appropriate characteristic velocities, it is capable of accurately solving problems that cannot be approximated as ideal magnetohydrodynamics and whose characteristic velocity is much lower than light velocity. The numerical results show that our scheme can solve the above problems as well as nearly ideal MHD problems. Our new scheme is particularly well suited to systems with initially weak magnetic field, and mixed phenomena of relativistic and non-relativistic velocity; for example, MRI in accretion disk, and super Alfvénic turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Three-Dimensional Simulations of MHD Turbulence Behind Relativistic Shock Waves and Their Implications for GRBs

Relativistic astrophysical phenomena such as gamma-ray bursts (GRBs) and active galactic nuclei often require long-lived strong magnetic field. Here, we report on three-dimensional special-relativistic magnetohydrodynamic (MHD) simulations to explore the amplification and decay of macroscopic turbulence dynamo excited by the so-called Richtmyer-Meshkov instability (RMI; a Rayleigh-Taylor type instability). This instability is an inevitable outcome of interactions between shock and ambient density fluctuations. We find that the magnetic energy grows exponentially in a few eddy turnover times, and then, following the decay of kinetic turbulence, decays with a temporal power-law exponent of -0.7. The magnetic-energy fraction can reach $epsilon_B \sim$ 0.1 but depends on the initial magnetic field strength. We find that the magnetic energy grows by at least two orders of magnitude compared to the magnetic energy immediately behind the shock. This minimum degree of the amplification does not depend on the amplitude of the initial density fluctuations, while the growth timescale and the maximum magnetic energy depend on the degree of inhomogeneity in the density. The transition from Kolmogorov cascade to MHD critical balance cascade occurs at $\sim$ 1/10th the initial inhomogeneity scale, which limits the maximum synchrotron polarization to less than 2%. New results include the avoidance of electron cooling with RMI turbulence, the turbulent photosphere model via RMI, and the shallow decay of the early afterglow from RMI. We also performed a simulation of freely decaying turbulence with relativistic velocity dispersion. We find that relativistic turbulence begins to decay much faster than one eddy-turnover time because of fast shock dissipation, which does not support the relativistic turbulence model by Narayan & Kumar.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Toward Understanding the Cosmic-ray Acceleration at Young Supernova Remnants Interacting with Interstellar Clouds: Possible Applications to RX J1713.7-3946

Using three-dimensional magnetohydrodynamics simulations, we investigate general properties of a blast wave shock interacting with interstellar clouds. The shock-cloud interaction generates a turbulent shell through the vorticity generations. In the turbulent shell, the magnetic field is amplified as a result of turbulent dynamo action. In the case of a young supernova remnant, the corresponding strength of the magnetic field is approximately 1 mG. The propagation speed of the shock wave is significantly stalled in the clouds because of the high density, while the shock maintains a high velocity in the diffuse surrounding. In addition, when the shock wave hits the clouds, reflection shock waves are generated that propagate back into the shocked shell. From these simulation results, many observational characteristics of a young SNR RX J1713.7-3946 that is suggested to be interacting with molecular clouds, can be explained as follows: The reflection shocks can accelerate particles in the turbulent downstream region where the magnetic field strength reaches 1mG, which causes short-time variability of synchrotron X-rays. Since the shock velocity is stalled locally in the clouds, the temperature in the shocked cloud is suppressed below 1 keV. Thus, thermal X-ray line emissions would be faint even if the SNR is interacting with molecular clouds. We also find that the photon index of the $π^0$-decay gamma rays generated by cosmic-ray protons can be 1.5, because the penetration depth of high-energy particles into the clumpy clouds depends on their energy, which is consistent with the recent gamma-ray observation by Fermi space telescope.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Two-Step Acceleration Model of Cosmic Rays at Middle-Aged SNR

Recent gamma-ray observations of middle-aged supernova remnants revealed a mysterious broken power-law spectrum. Using three-dimensional magnetohydrodynamics simulations, we show that the interaction between a supernova blast wave and interstellar clouds formed by thermal instability generates multiple reflected shocks. The typical Mach numbers of the reflected shocks are shown to be M ~ 2 depending on the density contrast between the diffuse intercloud gas and clouds. These secondary shocks can further energize cosmic-ray particles originally accelerated at the blast-wave shock. This "two-step" acceleration scenario reproduces the observed gamma-ray spectrum and predicts the high-energy spectral index ranging approximately from 3 to 4.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Turbulence and Magnetic Field Amplification in Supernova Remnants: Interactions Between A Strong Shock Wave and Multi-Phase Interstellar Medium

We examine MHD simulations of the propagation of a strong shock wave through the interstellar two-phase medium composed of small-scale cloudlets and diffuse warm neutral medium in two-dimensional geometry. The pre-shock two-phase medium is provided as a natural consequence of the thermal instability that is expected to be ubiquitous in the interstellar medium. We show that the shock-compressed shell becomes turbulent owing to the preshock density inhomogeneity and magnetic field amplification takes place in the shell. The maximum field strength is determined by the condition that plasma beta ~ 1, which gives the field strength on the order of 1 mG in the case of shock velocity ~ 1,000 km/s. The strongly magnetized region shows filamentary and knot-like structures in two-dimensional simulations. The spatial scale of the regions with magnetic field of 1 mG in our simulation is roughly 0.05 pc which is comparable to the spatial scale of the X-ray hot spots recently discovered in supernova remnants where the magnetic field strength is indicated to be amplified up to the order of 1 mG. This result may also suggest that the turbulent region with locally strong magnetic field is expected to be spread out in the region with frequent supernova explosions, such as in the Galactic center and starburst galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Two-Fluid MHD Simulations of Converging Hi Flows in the Interstellar Medium. II: Are Molecular Clouds Generated Directly from Warm Neutral Medium?

Formation of interstellar clouds as a consequence of thermal instability is studied using two-dimensional two-fluid magnetohydrodynamic simulations. We consider the situation of converging, supersonic flows of warm neutral medium in the interstellar medium that generate a shocked slab of thermally unstable gas in which clouds form. We found, as speculated in paper I, that in the shocked slab magnetic pressure dominates thermal pressure and the thermal instability grows in the isochorically cooling, thermally unstable slab that leads formation of HI clouds whose number density is typically n < 100 cm^-3, even if the angle between magnetic field and converging flows is small. We also found that even if there is a large dispersion of magnetic field, evolution of the shocked slab is essentially determined by the angle between the mean magnetic field and converging flows. Thus, the direct formation of molecular clouds by piling up warm neutral medium does not seem a typical molecular cloud formation process, unless the direction of supersonic converging flows is biased to the orientation of mean magnetic field by some mechanism. However, when the angle is small, the HI shell generated as a result of converging flows is massive and possibly evolves into molecular clouds, provided gas in the massive HI shell is piled up again along the magnetic field line. We expect that another subsequent shock wave can pile up again the gas of the massive shell and produce a larger cloud. We thus emphasize the importance of multiple episodes of converging flows, as a typical formation process of molecular clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Two-Fluid MHD Simulations of Converging HI Flows in the Interstellar Medium. I: Methodology and Basic Results

We develop an unconditionally stable numerical method for solving the coupling between two fluids (frictional forces/heatings, ionization, and recombination), and investigate the dynamical condensation process of thermally unstable gas that is provided by the shock waves in a weakly ionized and magnetized interstellar medium by using two-dimensional two-fluid magnetohydrodynamical simulations. If we neglect the effect of magnetic field, it is known that condensation driven by thermal instability can generate high density clouds whose physical condition corresponds to molecular clouds (precursor of molecular clouds). In this paper, we study the effect of magnetic field on the evolution of supersonic converging HI flows and focus on the case in which the orientation of magnetic field to converging flows is orthogonal. We show that the magnetic pressure gradient parallel to the flows prevents the formation of high density and high column density clouds, but instead generates fragmented, filamentary HI clouds. With this restricted geometry, magnetic field drastically diminishes the opportunity of fast molecular cloud formation directly from the warm neutral medium, in contrast to the case without magnetic field.

0 citations0 reviewsNo code linkedOpen access