TY - JOUR A1 - Berenstein, Igal A1 - Beta, Carsten T1 - Flow-induced control of chemical turbulence JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - We report spatiotemporal chaos in the Oregonator model of the Belousov-Zhabotinsky reaction. Spatiotemporal chaos spontaneously develops in a regime, where the underlying local dynamics show stable limit cycle oscillations (diffusion-induced turbulence). We show that spatiotemporal chaos can be suppressed by a unidirectional flow in the system. With increasing flow velocity, we observe a transition scenario from spatiotemporal chaos via a regime of travelling waves to a stationary steady state. At large flow velocities, we recover the known regime of flow distributed oscillations. KW - chaos KW - chemical equilibrium KW - chemically reactive flow KW - reaction kinetics theory KW - spatiotemporal phenomena KW - turbulence Y1 - 2011 U6 - https://doi.org/10.1063/1.3656248 SN - 0021-9606 VL - 135 IS - 16 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Brose, Robert A1 - Sushch, Iuri A1 - Pohl, Martin A1 - Luken, K. J. A1 - Filipovic, M. D. A1 - Lin, R. T1 - Nonthermal emission from the reverse shock of the youngest galactic supernova remnant G1.9+0.3 JF - Astronomy and astrophysics : an international weekly journal N2 - Context. The youngest Galactic supernova remnant G1.9+0.3 is an interesting target for next-generation gamma-ray observatories. So far, the remnant is only detected in the radio and the X-ray bands, but its young age of approximate to 100 yr and inferred shock speed of approximate to 14 000 km s(-1) could make it an efficient particle accelerator. Aims. We aim to model the observed radio and X-ray spectra together with the morphology of the remnant. At the same time, we aim to estimate the gamma-ray flux from the source and evaluate the prospects of its detection with future gamma-ray experiments. Methods. We performed spherical symmetric 1D simulations with the RATPaC code, in which we simultaneously solved the transport equation for cosmic rays, the transport equation for magnetic turbulence, and the hydro-dynamical equations for the gas flow. Separately computed distributions of the particles accelerated at the forward and the reverse shock were then used to calculate the spectra of synchrotron, inverse Compton, and pion-decay radiation from the source. Results. The emission from G1.9+0.3 can be self-consistently explained within the test-particle limit. We find that the X-ray flux is dominated by emission from the forward shock while most of the radio emission originates near the reverse shock, which makes G1.9+0.3 the first remnant with nonthermal radiation detected from the reverse shock. The flux of very-high-energy gamma-ray emission from G1.9+0.3 is expected to be close to the sensitivity threshold of the Cherenkov Telescope Array. The limited time available to grow large-scale turbulence limits the maximum energy of particles to values below 100 TeV, hence G1.9+0.3 is not a PeVatron. KW - acceleration of particles KW - turbulence KW - ISM: supernova remnants KW - gamma rays: ISM Y1 - 2019 U6 - https://doi.org/10.1051/0004-6361/201834430 SN - 1432-0746 VL - 627 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Brose, Robert A1 - Telezhinsky, Igor O. A1 - Pohl, Martin T1 - Transport of magnetic turbulence in supernova remnants JF - Physical review letters N2 - Context. Supernova remnants are known as sources of Galactic cosmic rays for their nonthermal emission of radio waves, X-rays, and gamma rays. However, the observed soft broken power-law spectra are hard to reproduce within standard acceleration theory based on the assumption of Bohm diffusion and steady-state calculations. Aims. We point out that a time-dependent treatment of the acceleration process together with a self-consistent treatment of the scattering turbulence amplification is necessary. Methods. We numerically solve the coupled system of transport equations for cosmic rays and isotropic Alfvenic turbulence. The equations are coupled through the growth rate of turbulence determined by the cosmic-ray gradient and the spatial diffusion coefficient of cosmic rays determined by the energy density of the turbulence. The system is solved on a comoving expanding grid extending upstream for dozens of shock radii, allowing for the self-consistent study of cosmic-ray diffusion in the vicinity of their acceleration site. The transport equation for cosmic rays is solved in a test-particle approach. Results. We demonstrate that the system is typically not in a steady state. In fact, even after several thousand years of evolution, no equilibrium situation is reached. The resulting time-dependent particle spectra strongly differ from those derived assuming a steady state and Bohm diffusion. Our results indicate that proper accounting for the evolution of the scattering turbulence and hence the particle diffusion coefficient is crucial for the formation of the observed soft spectra. In any case, the need to continuously develop magnetic turbulence upstream of the shock introduces nonlinearity in addition to that imposed by cosmic-ray feedback. KW - ISM: supernova remnants KW - acceleration of particles KW - turbulence Y1 - 2016 U6 - https://doi.org/10.1051/0004-6361/201527345 SN - 1432-0746 VL - 593 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Ge, J. X. A1 - He, J. H. A1 - Yan, Huirong T1 - Effects of turbulent dust grain motion to interstellar chemistry JF - Monthly notices of the Royal Astronomical Society N2 - Theoretical studies have revealed that dust grains are usually moving fast through the turbulent interstellar gas, which could have significant effects upon interstellar chemistry by modifying grain accretion. This effect is investigated in this work on the basis of numerical gas-grain chemical modelling. Major features of the grain motion effect in the typical environment of dark clouds (DC) can be summarized as follows: (1) decrease of gas-phase (both neutral and ionic) abundances and increase of surface abundances by up to 2-3 orders of magnitude; (2) shifts of the existing chemical jumps to earlier evolution ages for gas-phase species and to later ages for surface species by factors of about 10; (3) a few exceptional cases in which some species turn out to be insensitive to this effect and some other species can show opposite behaviours too. These effects usually begin to emerge from a typical DC model age of about 10(5) yr. The grain motion in a typical cold neutral medium (CNM) can help overcome the Coulomb repulsive barrier to enable effective accretion of cations on to positively charged grains. As a result, the grain motion greatly enhances the abundances of some gas-phase and surface species by factors up to 2-6 or more orders of magnitude in the CNM model. The grain motion effect in a typical molecular cloud (MC) is intermediate between that of the DC and CNM models, but with weaker strength. The grain motion is found to be important to consider in chemical simulations of typical interstellar medium. KW - astrochemistry KW - turbulence KW - ISM: abundances KW - ISM: clouds KW - dust, extinction KW - ISM: molecules Y1 - 2016 U6 - https://doi.org/10.1093/mnras/stv2560 SN - 0035-8711 SN - 1365-2966 VL - 455 SP - 3570 EP - 3587 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Kobzar, Oleh A1 - Niemiec, Jacek A1 - Pohl, Martin A1 - Bohdan, Artem T1 - Spatio-temporal evolution of the non-resonant instability in shock precursors of young supernova remnants JF - Monthly notices of the Royal Astronomical Society N2 - A non-resonant cosmic ray (CR) current-driven instability may operate in the shock precursors of young supernova remnants and be responsible for magnetic-field amplification, plasma heating and turbulence. Earlier simulations demonstrated magnetic-field amplification, and in kinetic studies a reduction of the relative drift between CRs and thermal plasma was observed as backreaction. However, all published simulations used periodic boundary conditions, which do not account for mass conservation in decelerating flows and only allow the temporal development to be studied. Here we report results of fully kinetic particle-in-cell simulations with open boundaries that permit inflow of plasma on one side of the simulation box and outflow at the other end, hence allowing an investigation of both the temporal and the spatial development of the instability. Magnetic-field amplification proceeds as in studies with periodic boundaries and, observed here for the first time, the reduction of relative drifts causes the formation of a shock-like compression structure at which a fraction of the plasma ions are reflected. Turbulent electric field generated by the non-resonant instability inelastically scatters CRs, modifying and anisotropizing their energy distribution. Spatial CR scattering is compatible with Bohm diffusion. Electromagnetic turbulence leads to significant non-adiabatic heating of the background plasma maintaining bulk equipartition between ions and electrons. The highest temperatures are reached at sites of large-amplitude electrostatic fields. Ion spectra show supra-thermal tails resulting from stochastic scattering in the turbulent electric field. Together, these modifications in the plasma flow will affect the properties of the shock and particle acceleration there. KW - acceleration of particles KW - shock waves KW - turbulence KW - methods: numerical KW - cosmic rays KW - ISM: supernova remnants Y1 - 2017 U6 - https://doi.org/10.1093/mnras/stx1201 SN - 0035-8711 SN - 1365-2966 VL - 469 SP - 4985 EP - 4998 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Lebiga, O. A1 - Santos-Lima, Reinaldo A1 - Yan, Huirong T1 - Kinetic-MHD simulations of gyroresonance instability driven by CR pressure anisotropy JF - Monthly notices of the Royal Astronomical Society N2 - The transport of cosmic rays (CRs) is crucial for the understanding of almost all high-energy phenomena. Both pre-existing large-scale magnetohydrodynamic (MHD) turbulence and locally generated turbulence through plasma instabilities are important for the CR propagation in astrophysical media. The potential role of the resonant instability triggered by CR pressure anisotropy to regulate the parallel spatial diffusion of low-energy CRs (less than or similar to 100 GeV) in the interstellar and intracluster medium of galaxies has been shown in previous theoretical works. This work aims to study the gyroresonance instability via direct numerical simulations, in order to access quantitatively the wave-particle scattering rates. For this, we employ a 1D PIC-MHD code to follow the growth and saturation of the gyroresonance instability. We extract from the simulations the pitch-angle diffusion coefficient D-mu mu produced by the instability during the linear and saturation phases, and a very good agreement (within a factor of 3) is found with the values predicted by the quasi-linear theory (QLT). Our results support the applicability of the QLT for modelling the scattering of low-energy CRs by the gyroresonance instability in the complex interplay between this instability and the large-scale MHD turbulence. KW - MHD KW - plasmas KW - turbulence KW - cosmic rays Y1 - 2018 U6 - https://doi.org/10.1093/mnras/sty309 SN - 0035-8711 SN - 1365-2966 VL - 476 IS - 2 SP - 2779 EP - 2791 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Mizuno, Yosuke A1 - Pohl, Martin A1 - Niemiec, Jacek A1 - Zhang, Bing A1 - Nishikawa, Ken-Ichi A1 - Hardee, Philip E. T1 - Magnetic-field amplification by turbulence in a relativistic shockpropagating through an inhomogeneous medium JF - The astrophysical journal : an international review of spectroscopy and astronomical physics N2 - We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the timescale of magnetic-field growth depends on the shock strength. KW - gamma-ray burst: general KW - magnetohydrodynamics (MHD) KW - methods: numerical KW - relativistic processes KW - shock waves KW - turbulence Y1 - 2011 U6 - https://doi.org/10.1088/0004-637X/726/2/62 SN - 0004-637X VL - 726 IS - 2 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Mizuno, Yosuke A1 - Pohl, Martin A1 - Niemiec, Jacek A1 - Zhang, Bing A1 - Nishikawa, Ken-Ichi A1 - Hardee, Philip E. T1 - Magnetic field amplification and saturation in turbulence behind a relativistic shock JF - Monthly notices of the Royal Astronomical Society N2 - We have investigated via 2D relativistic magnetohydrodynamic simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the post-shock region, magnetic field is strongly amplified by turbulent motions triggered by pre-shock density inhomogeneities. Using a long-simulation box we have followed the magnetic field amplification until it is fully developed and saturated. The turbulent velocity is subrelativistic even for a strong shock. Magnetic field amplification is controlled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic field amplification and saturation depend on the initial strength and direction of the magnetic field in the pre-shock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the post-shock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the post-shock region. If the initial magnetic field in the pre-shock medium is strong, the post-shock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the post-shock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei. KW - MHD KW - relativistic processes KW - shock waves KW - turbulence KW - methods: numerical KW - gamma-ray burst: general Y1 - 2014 U6 - https://doi.org/10.1093/mnras/stu196 SN - 0035-8711 SN - 1365-2966 VL - 439 IS - 4 SP - 3490 EP - 3503 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Pohl, Manuela A1 - Wilhelm, Alina A1 - Telezhinsky, Igor O. T1 - Reacceleration of electrons in supernova remnants JF - Astronomy and astrophysics : an international weekly journal N2 - Context. radio spectra of many shell-type supernova remnants show deviations from those expected on theoretical grounds. Aims. In this paper we determine the effect of stochastic reacceleration on the spectra of electrons in the GeV band and at lower energies, and we investigate whether reacceleration can explain the observed variation in radio spectral indices. Methods. We explicitely calculated the momentum diffusion coefficient for 3 types of turbulence expected downstream of the forward shock: fast-mode waves, small-scale non-resonant modes, and large-scale modes arising from turbulent dynamo activity. After noting that low-energy particles are efficiently coupled to the quasi-thermal plasma, a simplified cosmic-ray transport equation can be formulated and is numerically solved. Results. Only fast-mode waves can provide momentum diffusion fast enough to significantly modify the spectra of particles. Using a synchrotron emissivity that accurately reflects a highly turbulent magnetic field, we calculated the radio spectral index and find that soft spectra with index a alpha less than or similar to -0.6 can be maintained over more than 2 decades in radio frequency, even if the electrons experience reacceleration for only one acceleration time. A spectral hardening is possible but considerably more frequency-dependent. The spectral modification imposed by stochastic reacceleration downstream of the forward shock depends only weakly on the initial spectrum provided by, e.g., diffusive shock acceleration at the shock itself. KW - acceleration of particles KW - turbulence KW - cosmic rays KW - ISM: supernova remnants Y1 - 2015 U6 - https://doi.org/10.1051/0004-6361/201425027 SN - 0004-6361 SN - 1432-0746 VL - 574 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Pratt, Jane A1 - Busse, Angela A1 - Mueller, W-C A1 - Watkins, Nikolas W. A1 - Chapman, Sandra C. T1 - Extreme-value statistics from Lagrangian convex hull analysis for homogeneous turbulent Boussinesq convection and MHD convection JF - New journal of physics : the open-access journal for physics N2 - We investigate the utility of the convex hull of many Lagrangian tracers to analyze transport properties of turbulent flows with different anisotropy. In direct numerical simulations of statistically homogeneous and stationary Navier-Stokes turbulence, neutral fluid Boussinesq convection, and MHD Boussinesq convection a comparison with Lagrangian pair dispersion shows that convex hull statistics capture the asymptotic dispersive behavior of a large group of passive tracer particles. Moreover, convex hull analysis provides additional information on the sub-ensemble of tracers that on average disperse most efficiently in the form of extreme value statistics and flow anisotropy via the geometric properties of the convex hulls. We use the convex hull surface geometry to examine the anisotropy that occurs in turbulent convection. Applying extreme value theory, we show that the maximal square extensions of convex hull vertices are well described by a classic extreme value distribution, the Gumbel distribution. During turbulent convection, intermittent convective plumes grow and accelerate the dispersion of Lagrangian tracers. Convex hull analysis yields information that supplements standard Lagrangian analysis of coherent turbulent structures and their influence on the global statistics of the flow. KW - turbulence KW - magnetohydrodynamics KW - Lagrangian statistics KW - magnetoconvection KW - turbulent transport Y1 - 2017 U6 - https://doi.org/10.1088/1367-2630/aa6fe8 SN - 1367-2630 VL - 19 PB - IOP Publ. Ltd. CY - Bristol ER -