@article{WarrenEllisonBarkovetal.2017, author = {Warren, Donald C. and Ellison, Donald C. and Barkov, Maxim V. and Nagataki, Shigehiro}, title = {Nonlinear Particle Acceleration and Thermal Particles in GRB Afterglows}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {835}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {2}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/1538-4357/aa56c3}, pages = {15}, year = {2017}, abstract = {The standard model for GRB afterglow emission treats the accelerated electron population as a simple power law, N(E) proportional to E-p for p greater than or similar to 2. However, in standard Fermi shock acceleration, a substantial fraction of the swept-up particles do not enter the acceleration process at all. Additionally, if acceleration is efficient, then the nonlinear back-reaction of accelerated particles on the shock structure modifies the shape of the nonthermal tail of the particle spectra. Both of these modifications to the standard synchrotron afterglow impact the luminosity, spectra, and temporal variation of the afterglow. To examine the effects of including thermal particles and nonlinear particle acceleration on afterglow emission, we follow a hydrodynamical model for an afterglow jet and simulate acceleration at numerous points during the evolution. When thermal particles are included, we find that the electron population is at no time well fitted by a single power law, though the highest-energy electrons are; if the acceleration is efficient, then the power-law region is even smaller. Our model predicts hard-soft-hard spectral evolution at X-ray energies, as well as an uncoupled X-ray and optical light curve. Additionally, we show that including emission from thermal particles has drastic effects (increases by factors of 100 and 30, respectively) on the observed flux at optical and GeV energies. This enhancement of GeV emission makes afterglow detections by future gamma-ray observatories, such as CTA, very likely.}, language = {en} } @article{ZhangYan2018, author = {Zhang, Heshou and Yan, Huirong}, title = {Polarization of submillimetre lines from interstellar medium}, series = {Monthly notices of the Royal Astronomical Society}, volume = {475}, journal = {Monthly notices of the Royal Astronomical Society}, number = {2}, publisher = {Oxford University Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stx3164}, pages = {2415 -- 2420}, year = {2018}, abstract = {Magnetic fields play important roles in many astrophysical processes. However, there is no universal diagnostic for the magnetic fields in the interstellar medium (ISM) and each magnetic tracer has its limitation. Any new detection method is thus valuable. Theoretical studies have shown that submillimetre fine-structure lines are polarized due to atomic alignment by ultraviolet photon-excitation, which opens up a new avenue to probe interstellar magnetic fields. We will, for the first time, perform synthetic observations on the simulated three-dimensional ISM to demonstrate the measurability of the polarization of submillimetre atomic lines. The maximum polarization for different absorption and emission lines expected from various sources, including star-forming regions are provided. Our results demonstrate that the polarization of submillimetre atomic lines is a powerful magnetic tracer and add great value to the observational studies of the submilimetre astronomy.}, language = {en} } @article{LebigaSantosLimaYan2018, author = {Lebiga, O. and Santos-Lima, Reinaldo and Yan, Huirong}, title = {Kinetic-MHD simulations of gyroresonance instability driven by CR pressure anisotropy}, series = {Monthly notices of the Royal Astronomical Society}, volume = {476}, journal = {Monthly notices of the Royal Astronomical Society}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/sty309}, pages = {2779 -- 2791}, year = {2018}, abstract = {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.}, language = {en} } @article{VafinLazarFichtneretal.2018, author = {Vafin, Sergei and Lazar, M. and Fichtner, H. and Schlickeiser, R. and Drillisch, M.}, title = {Solar wind temperature anisotropy constraints from streaming instabilities}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {613}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201731852}, pages = {5}, year = {2018}, abstract = {Due to the relatively low rate of particle-particle collisions in the solar wind, kinetic instabilities (e.g., the mirror and firehose) play an important role in regulating large deviations from temperature isotropy. These instabilities operate in the high beta(parallel to) > 1 plasmas, and cannot explain the other limits of the temperature anisotropy reported by observations in the low beta beta(parallel to) < 1 regimes. However, the instability conditions are drastically modified in the presence of streaming (or counterstreaming) components, which are ubiquitous in space plasmas. These effects have been analyzed for the solar wind conditions in a large interval of heliospheric distances, 0.3-2.5 AU. It was found that proton counter-streams are much more crucial for plasma stability than electron ones. Moreover, new instability thresholds can potentially explain all observed bounds on the temperature anisotropy, and also the level of differential streaming in the solar wind.}, language = {en} } @article{BroseSushchPohletal.2019, author = {Brose, Robert and Sushch, Iuri and Pohl, Martin and Luken, K. J. and Filipovic, M. D. and Lin, R.}, title = {Nonthermal emission from the reverse shock of the youngest galactic supernova remnant G1.9+0.3}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {627}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201834430}, pages = {9}, year = {2019}, abstract = {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.}, language = {en} } @article{RuedigerKuekerKapylaetal.2019, author = {R{\"u}diger, G{\"u}nther and K{\"u}ker, Manfred and Kapyla, P. J. and Strassmeier, Klaus G.}, title = {Antisolar differential rotation of slowly rotating cool stars}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {630}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201935280}, pages = {9}, year = {2019}, abstract = {Rotating stellar convection transports angular momentum towards the equator, generating the characteristic equatorial acceleration of the solar rotation while the radial flux of angular momentum is always inwards. New numerical box simulations for the meridional cross-correlation < u(theta)u(phi)>, however, reveal the angular momentum transport towards the poles for slow rotation and towards the equator for fast rotation. The explanation is that for slow rotation a negative radial gradient of the angular velocity always appears, which in combination with a so-far neglected rotation-induced off-diagonal eddy viscosity term nu(perpendicular to) provides "antisolar rotation" laws with a decelerated equator Similarly, the simulations provided positive values for the rotation-induced correlation < u(r)u(theta)>, which is relevant for the resulting latitudinal temperature profiles (cool or warm poles) for slow rotation and negative values for fast rotation. Observations of the differential rotation of slowly rotating stars will therefore lead to a better understanding of the actual stress-strain relation, the heat transport, and the underlying model of the rotating convection.}, language = {en} } @article{PrattBusseMuelleretal.2017, author = {Pratt, Jane and Busse, Angela and Mueller, W-C and Watkins, Nikolas W. and Chapman, Sandra C.}, title = {Extreme-value statistics from Lagrangian convex hull analysis for homogeneous turbulent Boussinesq convection and MHD convection}, series = {New journal of physics : the open-access journal for physics}, volume = {19}, journal = {New journal of physics : the open-access journal for physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/aa6fe8}, pages = {18}, year = {2017}, abstract = {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.}, language = {en} } @article{KobzarNiemiecPohletal.2017, author = {Kobzar, Oleh and Niemiec, Jacek and Pohl, Martin and Bohdan, Artem}, title = {Spatio-temporal evolution of the non-resonant instability in shock precursors of young supernova remnants}, series = {Monthly notices of the Royal Astronomical Society}, volume = {469}, journal = {Monthly notices of the Royal Astronomical Society}, publisher = {Oxford Univ. Press}, address = {Oxford}, organization = {ANTARES Collaboration;H E S S Collaboration}, issn = {0035-8711}, doi = {10.1093/mnras/stx1201}, pages = {4985 -- 4998}, year = {2017}, abstract = {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.}, language = {en} } @article{VafinRiazantsevaYoon2017, author = {Vafin, Sergei and Riazantseva, M. and Yoon, P. H.}, title = {Kinetic Features in the Ion Flux Spectrum}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {850}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/1538-4357/aa9519}, pages = {7}, year = {2017}, abstract = {An interesting feature of solar wind fluctuations is the occasional presence of a well-pronounced peak near the spectral knee. These peaks are well investigated in the context of magnetic field fluctuations in the magnetosheath and they are typically related to kinetic plasma instabilities. Recently, similar peaks were observed in the spectrum of ion flux fluctuations of the solar wind and magnetosheath. In this paper, we propose a simple analytical model to describe such peaks in the ion flux spectrum based on the linear theory of plasma fluctuations. We compare our predictions with a sample observation in the solar wind. For the given observation, the peak requires similar to 10 minutes to grow up to the observed level that agrees with the quasi-linear relaxation time. Moreover, our model well reproduces the form of the measured peak in the ion flux spectrum. The observed lifetime of the peak is about 50 minutes, which is relatively close to the nonlinear Landau damping time of 30-40 minutes. Overall, our model proposes a plausible scenario explaining the observation.}, language = {en} } @article{GeHeYan2016, author = {Ge, J. X. and He, J. H. and Yan, Huirong}, title = {Effects of turbulent dust grain motion to interstellar chemistry}, series = {Monthly notices of the Royal Astronomical Society}, volume = {455}, journal = {Monthly notices of the Royal Astronomical Society}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stv2560}, pages = {3570 -- 3587}, year = {2016}, abstract = {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.}, language = {en} }