@article{MaitiMakwanaZhangetal.2022, author = {Maiti, Snehanshu and Makwana, Kirit and Zhang, Heshou and Yan, Huirong}, title = {Cosmic-ray transport in magnetohydrodynamic turbulence}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics / part 1}, volume = {926}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics / part 1}, number = {1}, publisher = {Institute of Physics Publ.}, address = {London}, issn = {1538-4357}, doi = {10.3847/1538-4357/ac46c8}, pages = {8}, year = {2022}, abstract = {This paper studies cosmic-ray (CR) transport in magnetohydrodynamic (MHD) turbulence. CR transport is strongly dependent on the properties of the magnetic turbulence. We perform test particle simulations to study the interactions of CR with both total MHD turbulence and decomposed MHD modes. The spatial diffusion coefficients and the pitch angle scattering diffusion coefficients are calculated from the test particle trajectories in turbulence. Our results confirm that the fast modes dominate the CR propagation, whereas Alfven and slow modes are much less efficient and have shown similar pitch-angle scattering rates. We investigate the cross field transport on large and small scales. On large/global scales, normal diffusion is observed and the diffusion coefficient is suppressed by M-A(zeta) compared to the parallel diffusion coefficients, with zeta closer to 4 in Alfven modes than that in total turbulence, as theoretically expected. For the CR transport on scales smaller than the turbulence injection scale, both the local and global magnetic reference frames are adopted. Superdiffusion is observed on such small scales in all the cases. Particularly, CR transport in Alfven modes show clear Richardson diffusion in the local reference frame. The diffusion transitions smoothly from the Richardson's one with index 1.5 to normal diffusion as the particle mean free path decreases from lambda(parallel to) >> L to lambda(parallel to) << L, where L is the injection/coherence length of turbulence. Our results have broad applications to CRs in various astrophysical environments.}, language = {en} } @article{XuYanLazarian2016, author = {Xu, Siyao and Yan, Huirong and Lazarian, A.}, title = {DAMPING OF MAGNETOHYDRODYNAMIC TURBULENCE IN PARTIALLY IONIZED PLASMA: IMPLICATIONS FOR COSMIC RAY PROPAGATION}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {826}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/0004-637X/826/2/166}, pages = {32}, year = {2016}, abstract = {We study the damping processes of both incompressible and compressible magnetohydrodynamic (MHD) turbulence in a partially ionized medium. We start from the linear analysis of MHD waves, applying both single-fluid and two-fluid treatments. The damping rates derived from the linear analysis are then used in determining the damping scales of MHD turbulence. The physical connection between the damping scale of MHD turbulence and the cutoff boundary of linear MHD waves is investigated. We find two branches of slow modes propagating in ions and neutrals, respectively, below the damping scale of slow MHD turbulence, and offer a thorough discussion of their propagation and dissipation behavior. Our analytical results are shown to be applicable in a variety of partially ionized interstellar medium (ISM) phases and the solar chromosphere. The importance of neutral viscosity in damping the Alfvenic turbulence in the interstellar warm neutral medium and the solar chromosphere is demonstrated. As a significant astrophysical utility, we introduce damping effects to the propagation of cosmic rays in partially ionized ISM. The important role of turbulence damping in both transit-time damping and gyroresonance is identified.}, 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} } @article{LiuYanWangetal.2019, author = {Liu, Ruo-Yu and Yan, Huirong and Wang, Xiang-Yu and Shao, Shi and Li, Hui}, title = {Gamma-Ray production in the extended halo of the galaxy and possible implications for the origin of galactic cosmic rays}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {871}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/1538-4357/aaf567}, pages = {11}, year = {2019}, abstract = {Various studies have implied the existence of a gaseous halo around the Galaxy extending out to similar to 100 kpc. Galactic cosmic rays (CRs) that propagate to the halo, either by diffusion or by convection with the possibly existing large-scale Galactic wind, can interact with the gas therein and produce gamma-rays via proton-proton collision. We calculate the CR distribution in the halo and the gamma-ray flux, and explore the dependence of the result on model parameters such as diffusion coefficient, CR luminosity, and CR spectral index. We find that the current measurement of isotropic gamma-ray background (IGRB) at less than or similar to TeV with the Fermi Large Area Telescope already approaches a level that can provide interesting constraints on the properties of Galactic CR (e.g., with CR luminosity L-CR <= 1041 erg s(-1)). We also discuss the possibilities of the Fermi bubble and IceCube neutrinos originating from the proton-proton collision between CRs and gas in the halo, as well as the implication of our results for the baryon budget of the hot circumgalactic medium of our Galaxy. Given that the isotropic gamma-ray background is likely to be dominated by unresolved extragalactic sources, future telescopes may extract more individual sources from the IGRB, and hence put even more stringent restrictions on the relevant quantities (such as Galactic CR luminosity and baryon budget in the halo) in the presence of a turbulent halo that we consider.}, language = {en} } @article{LiuWangXueetal.2019, author = {Liu, Ruo-Yu and Wang, Kai and Xue, Rui and Taylor, Andrew M. and Wang, Xiang-Yu and Li, Zhuo and Yan, Huirong}, title = {Hadronuclear interpretation of a high-energy neutrino event coincident with a blazar flare}, series = {Physical review : D, Particles, fields, gravitation, and cosmology}, volume = {99}, journal = {Physical review : D, Particles, fields, gravitation, and cosmology}, number = {6}, publisher = {American Physical Society}, address = {Melville}, issn = {2470-0010}, doi = {10.1103/PhysRevD.99.063008}, pages = {11}, year = {2019}, abstract = {Although many high-energy neutrinos detected by the IceCube telescope are believed to have an extraterrestrial origin, their astrophysical sources remain a mystery. Recently, an unprecedented discovery of a high-energy muon neutrino event coincident with a multiwavelength flare from a blazar, TXS 0506 + 056, shed some light on the origin of the neutrinos. It is usually believed that a blazar is produced by a relativistic jet launched from an accreting supermassive black hole (SMBH). Here, we show that the high-energy neutrino event can be interpreted by the inelastic hadronuclear interactions between the accelerated cosmic-ray protons in the relativistic jet and the dense gas clouds in the vicinity of the SMBH. Such a scenario only requires a moderate proton power in the jet, which could be much smaller than that required in the conventional hadronic model which instead calls upon the photomeson process. Meanwhile, the flux of the multiwavelength flare from the optical to gamma-ray band can be well explained by invoking a second radiation zone in the jet at a larger distance to the SMBH. In our model, the neutrino emission lasts a shorter time than the multiwavelength flare, so the neutrino event is not necessarily correlated with the flare, but it is probably accompanied by a spectrum hardening above a few giga-electron-volt (GeV).}, 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{SantosdeLimaYandeGouveiaDalPinoetal.2016, author = {Santos de Lima, Reinaldo and Yan, Huirong and de Gouveia Dal Pino, E. M. and Lazarian, A.}, title = {Limits on the ion temperature anisotropy in the turbulent intracluster medium}, series = {Monthly notices of the Royal Astronomical Society}, volume = {460}, journal = {Monthly notices of the Royal Astronomical Society}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stw1079}, pages = {2492 -- 2504}, year = {2016}, abstract = {Turbulence in the weakly collisional intracluster medium (ICM) of galaxies is able to generate strong thermal velocity anisotropies in the ions (with respect to the local magnetic field direction), if the magnetic moment of the particles is conserved in the absence of Coulomb collisions. In this scenario, the anisotropic pressure magnetohydrodynamic (AMHD) turbulence shows a very different statistical behaviour from the standard MHD one and is unable to amplify seed magnetic fields. This is in contrast to previous cosmological MHD simulations that are successful in explaining the observed magnetic fields in the ICM. On the other hand, temperature anisotropies can also drive plasma instabilities that can relax the anisotropy. This work aims to compare the relaxation rate with the growth rate of the anisotropies driven by the turbulence. We employ quasi-linear theory to estimate the ion scattering rate resulting from the parallel firehose, mirror and ion-cyclotron instabilities, for a set of plasma parameters resulting from AMHD simulations of the turbulent ICM. We show that the ICM turbulence can sustain only anisotropy levels very close to the instability thresholds. We argue that the AMHD model that bounds the anisotropies at the marginal stability levels can describe the Alfv,nic turbulence cascade in the ICM.}, language = {en} } @article{ZhaoYanLiuetal.2022, author = {Zhao, Siqi Q. and Yan, Huirong and Liu, Terry Z. and Liu, Mingzhe and Wang, Huizi}, title = {Multispacecraft analysis of the properties of magnetohydrodynamic fluctuations in Sub-Alfvenic solar wind turbulence at 1 au}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {937}, 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/ac822e}, pages = {14}, year = {2022}, abstract = {We present observations of three-dimensional magnetic power spectra in wavevector space to investigate the anisotropy and scalings of sub-Alfvenic solar wind turbulence at magnetohydrodynamic (MHD) scale using the Magnetospheric Multiscale spacecraft. The magnetic power distributions are organized in a new coordinate determined by wavevectors ((kappa) over cap) and background magnetic field ((b) over cap (0)) in Fourier space. This study utilizes two approaches to determine wavevectors: the singular value decomposition method and multispacecraft timing analysis. The combination of the two methods allows an examination of the properties of magnetic field fluctuations in terms of mode compositions without any spatiotemporal hypothesis. Observations show that fluctuations (delta B-perpendicular to 1) in the direction perpendicular to (kappa) over cap and (b) over cap (0) prominently cascade perpendicular to (b) over cap (0), and such anisotropy increases with wavenumbers. The reduced power spectra of 6.8 11 follow Goldreich-Sridhar scalings: (P) over cap (k(perpendicular to)) proportional to k(perpendicular to)(-5/3) and (P) over cap (k(parallel to)) proportional to k(parallel to)(-2). In contrast, fluctuations within the (k) over cap(b) over cap (0) plane show isotropic behaviors: perpendicular power distributions are approximately the same as parallel distributions. The reduced power spectra of fluctuations within the (k) over cap(b) over cap (0) plane follow the scalings (P) over cap (k(perpendicular to)) proportional to k(perpendicular to)(-3/2) and (P) over cap (k(parallel to)) proportional to k(parallel to)(-3/2). Comparing frequency-wavevector spectra with theoretical dispersion relations of MHD modes, we find that delta B-perpendicular to 1 are probably associated with Alfven modes. On the other hand, magnetic field fluctuations within the (k) over cap(b) over cap (0) plane more likely originate from fast modes based on their isotropic behaviors. The observations of anisotropy and scalings of different magnetic field components are consistent with the predictions of current compressible MHD theory. Moreover, for the Alfvenic component, the ratio of cascading time to the wave period is found to be a factor of a few, consistent with critical balance in the strong turbulence regime. These results are valuable for further studies of energy compositions of plasma turbulence and their effects on energetic particle transport.}, language = {en} } @article{XueLiuWangetal.2019, author = {Xue, Rui and Liu, Ruo-Yu and Wang, Xiang-Yu and Yan, Huirong and B{\"o}ttcher, Markus}, title = {On the minimum jet power of TeV BL Lac objects in the p-gamma model}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {871}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/1538-4357/aaf720}, pages = {10}, year = {2019}, abstract = {We study the requirement of the jet power in the conventional p-gamma models (photopion production and Bethe-Heitler pair production) for TeV BL Lac objects. We select a sample of TeV BL Lac objects whose spectral energy distributions are difficult to explain by the one-zone leptonic model. Based on the relation between the p-gamma interaction efficiency and the opacity of gamma gamma absorption, we find that the detection of TeV emission poses upper limits on the p-gamma interaction efficiencies in these sources and hence minimum jet powers can be derived accordingly. We find that the obtained minimum jet powers exceed the Eddington luminosity of the supermassive black holes (SMBHs). Implications for the accretion mode of the SMBHs in these BL Lac objects and the origin of their TeV emissions are discussed.}, language = {en} } @article{LiuYan2020, author = {Liu, Ruoyu and Yan, Huirong}, title = {On the unusually large spatial extent of the TeV nebula HESS J1825-137}, series = {Monthly notices of the Royal Astronomical Society}, volume = {494}, journal = {Monthly notices of the Royal Astronomical Society}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/staa911}, pages = {2618 -- 2627}, year = {2020}, abstract = {Deep observation of the High Energy Stereoscopic System (HESS) on the most extended pulsar wind nebula HESS J1825-137 reveals an enhanced energy-dependent morphology, providing useful information on the particle transport mechanism in the nebula. We find that the energy-dependent morphology is consistent with a diffusion-dominated transport of electrons/positrons. It provides an alternative possible interpretation for the unusually large spatial extent (i.e. greater than or similar to 100 pc) of the nebula, which could then be attributed to the diffusion of escaping electrons/positrons from a compact plerion. The influence of various model parameters on the energy-dependent extent of the nebula is studied in the diffusion-dominated scenario. We also show that the energy-dependent morphology of the nebula may also be used to study the spin-down history of the pulsar.}, language = {en} }