Source author record

Hyesung Kang

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

35 published item(s)

preprint2021arXiv

Effects of Forcing on Shocks and Energy Dissipation in Interstellar and Intracluster Turbulences

Observations indicate that turbulence in the interstellar medium (ISM) is supersonic ($M_{\rm turb}\gg1$) and strongly magnetized ($β\sim0.01-1$), while in the intracluster medium (ICM) it is subsonic ($M_{\rm turb}\lesssim1$) and weakly magnetized ($β\sim100$). Here, $M_{\rm turb}$ is the turbulent Mach number and $β$ is the plasma beta. We study the properties of shocks induced in these disparate environments, including the distribution of the shock Mach number, $M_s$, and the dissipation of the turbulent energy at shocks, through numerical simulations using a high-order accurate code based on the WENO scheme. In particular, we investigate the effects of different modes of the forcing that drives turbulence: solenoidal, compressive, and a mixture of the two. In the ISM turbulence, while the density distribution looks different with different forcings, the velocity power spectrum, $P_v$, on small scales exhibits only weak dependence. Hence, the statistics of shocks depend weakly on forcing either. In the ISM models with $M_{\rm turb}\approx10$ and $β\sim0.1$, the fraction of the turbulent energy dissipated at shocks is estimated to be $\sim15~\%$, not sensitive to the forcing mode. In contrast, in the ICM turbulence, $P_v$ as well as the density distribution show strong dependence on forcing. The frequency and average Mach number of shocks are greater for compressive forcing than for solenoidal forcing, so is the energy dissipation. The fraction of ensuing shock dissipation is in the range of $\sim10-35~\%$ in the ICM models with $M_{\rm turb}\approx0.5$ and $β\sim10^6$. The rest of the turbulent energy should be dissipated through turbulent cascade.

preprint2021arXiv

Electron Preacceleration at Weak Quasi-perpendicular Intracluster Shocks: Effects of Preexisting Nonthermal Electrons

Radio relics in the outskirts of galaxy clusters imply the diffusive shock acceleration (DSA) of electrons at merger-driven shocks with Mach number $M_{s}\lesssim3-4$ in the intracluster medium (ICM). Recent studies have suggested that electron preacceleration and injection, prerequisite steps for DSA, could occur at supercritical shocks with $M_{s}\gtrsim2.3$ in the ICM, thanks to the generation of multiscale waves by microinstabilities such as the Alfvén ion cyclotron (AIC) instability, the electron firehose instability (EFI), and the whistler instability (WI). On the other hand, some relics are observed to have subcritical shocks with $M_{s}\lesssim2.3$, leaving DSA at such weak shocks as an outstanding problem. Reacceleration of preexisting nonthermal electrons has been contemplated as one of possible solutions for that puzzle. To explore this idea, we perform Particle-in-Cell (PIC) simulations for weak quasi-perpendicular shocks in high-$β$ ($β=P_{\rm gas}/P_{B}$) plasmas with power-law suprathermal electrons in addition to Maxwellian thermal electrons. We find that suprathermal electrons enhance the excitation of electron-scale waves via the EFI and WI. However, they do not affect the ion reflection and the ensuing generation of ion-scale waves via the AIC instability. The presence of ion-scale waves is the key for the preacceleration of electrons up to the injection momentum, thus the shock criticality condition for electron injection to DSA is preserved. Based on the results, we conclude that preexisting nonthermal electrons in the preshock region alone would not resolve the issue of electron preacceleration at subcritical ICM shocks.

preprint2020arXiv

Gamma-ray and Neutrino Emissions due to Cosmic-Ray Protons Accelerated at Intracluster Shocks in Galaxy Clusters

We examine the cosmic-ray protons (CRp) accelerated at collisionless shocks in galaxy clusters using cosmological structure formation simulations. We find that in the intracluster medium (ICM) within the virial radius of simulated clusters, only $\sim7$\% of shock kinetic energy flux is dissipated by the shocks that are expected to accelerate CRp, that is, supercritical, quasi-parallel ($Q_\parallel$) shocks with sonic Mach number $M_s\ge2.25$. The rest is dissipated at subcritical shocks and quasi-perpendicular shocks, both of which may not accelerate CRp. Adopting the diffusive shock acceleration (DSA) model recently presented in Ryu et al. (2019), we quantify the DSA of CRp in simulated clusters. The average fraction of the shock kinetic energy transferred to CRp via DSA is assessed at $\sim(1-2)\times10^{-4}$. We also examine the energization of CRp through reacceleration using a model based on the test-particle solution. Assuming that the ICM plasma passes through shocks three times on average through the history of the universe and that CRp are reaccelerated only at supercritical $Q_\parallel$-shocks, the CRp spectrum flattens by $\sim0.05-0.1$ in slope and the total amount of CRp energy increases by $\sim40-80$\% from reacceleration. We then estimate diffuse $γ$-ray and neutrino emissions, resulting from inelastic collisions between CRp and thermal protons. The predicted $γ$-ray emissions from simulated clusters lie mostly below the upper limits set by Fermi-LAT for observed clusters. The neutrino fluxes towards nearby clusters would be $\lesssim10^{-4}$ of the IceCube flux at $E_ν=1$ PeV and $\lesssim10^{-6}$ of the atmospheric neutrino flux in the energy range of $E_ν\leq1$ TeV.

preprint2020arXiv

Limiting the shock acceleration of cosmic-ray protons in the ICM

Observations of large-scale radio emissions prove the existence of shock accelerated cosmic-ray electrons in galaxy clusters, while the lack of detected $γ$-rays limits the acceleration of cosmic-ray protons in galaxy clusters. This challenges our understanding of how diffusive shock acceleration works. In this work, we couple the most updated recipes for shock acceleration in the intracluster medium to state-of-the-art magneto-hydrodynamical simulations of massive galaxy clusters. Furthermore, we use passive tracer particles to follow the evolution of accelerated cosmic-rays. We show that when the interplay between magnetic field topology and the feedback from accelerated cosmic rays is taken into account, the latest developments of particle acceleration theory give results which are compatible with observational constraints.

preprint2020arXiv

Reconstruction of Radio Relics and X-ray Tails in an Off-axis Cluster Merger: Hydrodynamical Simulations of A115

Although a merging galaxy cluster is a useful laboratory to study many important astrophysical processes and fundamental physics, only limited interpretations are possible without careful analysis of the merger scenario. However, the study is demanding because a thorough comparison of multi-wavelength observations with detailed numerical simulations is required. In this paper, we present such a study for the off-axis binary merger A115. The system possesses a number of remarkable observational features, but no convincing merger scenario, explaining the shape and location of the radio relic in harmony with the orientation of the cometary X-ray tails, has been presented. Our hydrodynamical simulation, with adaptive mesh refinement, suggests that the cometary X-ray tail of A115 might be a slingshot tail and can arise $\sim0.3$ Gyrs after the impact and before the two subclusters reach their apocenters. This scenario can predict the location and orientation of the giant radio relic, which is parallel to the northern X-ray tail. In addition, our study indicates that diffusive shock acceleration alone cannot generate the observed radio power unless aided by stronger magnetic fields and/or more significant presence of fossil electrons.

preprint2019arXiv

Chandra Observations of the Spectacular A3411-12 Merger Event

We present deep Chandra observations of A3411-12, a remarkable merging cluster that hosts the most compelling evidence for electron re-acceleration at cluster shocks to date. Using the $Y_X-M$ scaling relation, we find $r_{500} \sim 1.3$ Mpc, $M_{500} = (7.1 \pm 0.7) \times 10^{14} \ M_{\rm{\odot}}$, $kT=6.5\pm 0.1$ keV, and a gas mass of $M_{\rm g,500} = (9.7 \pm 0.1) \times 10^{13} M_\odot$. The gas mass fraction within $r_{500}$ is $f_{\rm g} = 0.14 \pm 0.01$. We compute the shock strength using density jumps to conclude that the Mach number of the merging subcluster is small ($M \leq 1.15_{-0.09}^{+0.14}$). We also present pseudo-density, projected temperature, pseudo-pressure, and pseudo-entropy maps. Based on the pseudo-entropy map we conclude that the cluster is undergoing a mild merger, consistent with the small Mach number. On the other hand, radio relics extend over Mpc scale in the A3411-12 system, which strongly suggests that a population of energetic electrons already existed over extended regions of the cluster.

preprint2019arXiv

Electron Firehose Instabilities in High-$β$ ICM shocks

The preacceleration of electrons through reflection and shock drift acceleration (SDA) is essential for the diffusive shock acceleration (DSA) of nonthermal electrons in collisionless shocks. Previous studies suggested that, in weak quasi-perpendicular ($Q_\perp$) shocks in the high-$β$ ($β=P_{\rm gas}/P_{\rm B}$) intracluster medium (ICM), the temperature anisotropy due to SDA-reflected electrons can drive the electron firehose instability, which excites oblique nonpropagating waves in the shock foot. In this paper, we investigate, through a linear analysis and particle-in-cell (PIC) simulations, the firehose instabilities driven by an electron temperature anisotropy (ETAFI) and also by a drifting electron beam (EBFI) in $β\sim100$ ICM plasmas. The EBFI should be more relevant in describing the self-excitation of upstream waves in $Q_\perp$-shocks, since backstreaming electrons in the shock foot behave more like an electron beam rather than an anisotropic bi-Maxwellian population. We find that the basic properties of the two instabilities, such as the growth rate, $γ$, and the wavenumber of fast-growing oblique modes are similar in the ICM environment, with one exception; while the waves excited by the ETAFI are nonpropagating ($ω_r=0$), those excited by the EBFI have a non-zero frequency ($ω_r\neq0$). However, the frequency is small with $ω_r<γ$. Thus, we conclude that the interpretation of previous studies for the nature of upstream waves based on the ETAFI remains valid in $Q_\perp$-shocks in the ICM.

preprint2019arXiv

Turbulence Dynamo in Galaxy Clusters

The existence of microgauss magnetic fields in galaxy clusters have been established through observations of synchrotron radiation and Faraday rotation. They are conjectured to be generated via small-scale dynamo by turbulent flow motions in the intracluster medium (ICM). Some of giant radio relics, on the other hand, show the structures of synchrotron polarization vectors, organized over the scales of $\sim$ Mpc, challenging the turbulence origin of cluster magnetic fields. Unlike turbulence in the interstellar medium, turbulence in the ICM is subsonic. And it is driven sporadically in highly stratified backgrounds, when major mergers occur during the hierarchical formation of clusters. To investigate quantitatively the characteristics of turbulence dynamo in such ICM environment, we performed a set of turbulence simulations using a high-order-accurate, magnetohydrodynamic (MHD) code. We find that turbulence dynamo could generate the cluster magnetic fields up to the observed level from the primordial seed fields of $10^{-15}$ G or so within the age of the universe, if the MHD description of the ICM could be extended down to $\sim$ kpc scales. However, highly organized structures of polarization vectors, such as those observed in the Sausage relic, are difficult to be reproduced by the shock compression of turbulence-generated magnetic fields. This implies that the modeling of giant radio relics may require the pre-existing magnetic fields organized over $\sim$ Mpc scales.

preprint2016arXiv

Re-acceleration Model for Radio Relics with Spectral Curvature

Most of the observed features of radio {\it gischt} relics such as spectral steepening across the relic width and power-law-like integrated spectrum can be adequately explained by diffusive shock acceleration (DSA) model, in which relativistic electrons are (re-)accelerated at shock waves induced in the intracluster medium. However, the steep spectral curvature in the integrated spectrum above $\sim 2$ GHz detected in some radio relics such as the Sausage relic in cluster CIZA J2242.8+5301 may not be interpreted by simple radiative cooling of postshock electrons. In order to understand such steepening, we here consider a model in which a spherical shock sweeps through and then exits out of a finite-size cloud with fossil relativistic electrons. The ensuing integrated radio spectrum is expected to steepen much more than predicted for aging postshock electrons, since the re-acceleration stops after the cloud-crossing time. Using DSA simulations that are intended to reproduce radio observations of the Sausage relic, we show that both the integrated radio spectrum and the surface brightness profile can be fitted reasonably well, if a shock of speed, $u_s \sim 2.5-2.8\times 10^3 \kms$, and sonic Mach number, $M_s \sim 2.7-3.0$, traverses a fossil cloud for $\sim 45$ Myr and the postshock electrons cool further for another $\sim 10$ Myr. This attempt illustrates that steep curved spectra of some radio gischt relics could be modeled by adjusting the shape of the fossil electron spectrum and adopting the specific configuration of the fossil cloud.

preprint2015arXiv

Curved Radio Spectra of Weak Cluster Shocks

In order to understand certain observed features of arc-like giant radio relics such as the rareness, uniform surface brightness, and curved integrated spectra, we explore a diffusive shock acceleration (DSA) model for radio relics in which a spherical shock impinges on a magnetized cloud containing fossil relativistic electrons. Toward this end, we perform DSA simulations of spherical shocks with the parameters relevant for the Sausage radio relic in cluster CIZA J2242.8+5301, and calculate the ensuing radio synchrotron emission from re-accelerated electrons. Three types of fossil electron populations are considered: a delta-function like population with the shock injection momentum, a power-law distribution, and a power-law with an exponential cutoff. The surface brightness profile of radio-emitting postshock region and the volume-integrated radio spectrum are calculated and compared with observations. We find that the observed width of the Sausage relic can be explained reasonably well by shocks with speed $u_s \sim 3\times 10^3 \kms$ and sonic Mach number $M_s \sim 3$. These shocks produce curved radio spectra that steepen gradually over $(0.1-10) ν_{\rm br}$ with break frequency $ ν_{\rm br}\sim 1$ GHz, if the duration of electron acceleration is $\sim 60 - 80$ Myr. However, the abrupt increase of spectral index above $\sim 1.5$ GHz observed in the Sausage relic seems to indicate that additional physical processes, other than radiative losses, operate for electrons with $γ_e \gtrsim 10^4$.

preprint2015arXiv

Radio and X-ray Shocks in Clusters of galaxies

Radio relics detected in the outskirts of galaxy clusters are thought to trace radio-emitting relativistic electrons accelerated at cosmological shocks. In this study, using the cosmological hydrodynamic simulation data for the large-scale structure formation and adopting a diffusive shock acceleration (DSA) model for the production of cosmic-ray (CR) electrons, we construct mock radio and X-ray maps of simulated galaxy clusters that are projected in the sky plane. Various properties of shocks and radio relics, including the shock Mach number, radio spectral index and luminosity are extracted from the synthetic maps and compared with observations. A substantial fraction of radio and X-ray shocks identified in these maps involve multiple shock surfaces along line of sights (LoSs), and the morphology of shock distributions in the maps depends on the projection direction. Among multiple shocks in a given LoS, radio observations tend to pick up stronger shocks with flatter radio spectra, while X-ray observations preferentially select weaker shocks with larger kinetic energy flux. As a result, the shock Mach numbers and locations derived from radio and X-ray observations could differ from each other in some cases. We also find that the distributions of the spectral index and radio power of the synthetic radio relics are somewhat inconsistent with those of observed real relics; a bit more radio relics have been observed closer to the cluster core and with steeper spectral indices. We suggest the inconsistency could be explained, if very weak shocks with $M_s \la 2$ accelerate CR electrons more efficiently, compared to the DSA model adopted here.

preprint2015arXiv

Radio emission from weak spherical shocks in the outskirts of galaxy clusters

In Kang (2015) we calculated the acceleration of cosmic-ray electrons and the ensuing radio synchrotron emission at weak spherical shocks that are expected to form in the outskirts of galaxy clusters.There we demonstrated that, at decelerating spherical shocks, the volume integrated spectra of both electrons and radiation deviate significantly from the test-particle power-laws predicted for constant planar shocks, because the shock compression ratio and the flux of injected electrons decrease in time. In this study, we consider spherical blast waves propagating into a constant density core surrounded by an isothermal halo with a decreasing density profile in order to explore how the deceleration rate of the shock speed affects the radio emission from accelerated electrons. The surface brightness profile and the volume-integrated radio spectrum of the model shocks are calculated by assuming a ribbon-like shock surface on a spherical shell and the associated downstream region of relativistic electrons. If the postshock magnetic field strength is about 7 microgauss, at the shock age of ~50 Myr, the volume-integrated radio spectrum steepens gradually with the spectral index from alpha_{inj} to alpha_{inj}+0.5 over 0.1-10 GHz, where alpha_{inj} is the injection index at the shock positionexpected from the diffusive shock acceleration theory. Such gradual steepening could explain the curved radio spectrum of the radio relic in cluster A2266, which was interpreted as a broken power-law by Trasatti et al. (2014), if the relic shock is young enough so that the break frequency falls in around 1 GHz.

preprint2014arXiv

Injection of $κ$-like Suprathermal Particles into Diffusive Shock Acceleration

We consider a phenomenological model for the thermal leakage injection in the diffusive shock acceleration (DSA) process, in which suprathermal protons and electrons near the shock transition zone are assumed to have the so-called $κ$-distributions produced by interactions of background thermal particles with pre-existing and/or self-excited plasma/MHD waves or turbulence. The $κ$-distribution has a power-law tail, instead of an exponential cutoff, well above the thermal peak momentum. So there are a larger number of potential seed particles with momentum, above that required for participation in the DSA process. As a result, the injection fraction for the $κ$-distribution depends on the shock Mach number much less severely compared to that for the Maxwellian distribution. Thus, the existence of $κ$-like suprathermal tails at shocks would ease the problem of extremely low injection fractions, especially for electrons and especially at weak shocks such as those found in the intracluster medium. We suggest that the injection fraction for protons ranges $10^{-4}-10^{-3}$ for a $κ$-distribution with $10 < κ_p < 30$ at quasi-parallel shocks, while the injection fraction for electrons becomes $10^{-6}-10^{-5}$ for a $κ$-distribution with $κ_e < 2$ at quasi-perpendicular shocks. For such $κ$ values the ratio of cosmic ray electrons to protons naturally becomes $K_{e/p}\sim 10^{-3}-10^{-2}$, which is required to explain the observed ratio for Galactic cosmic rays.

preprint2014arXiv

Shock Waves and Cosmic Ray Acceleration in the Outskirts of Galaxy Clusters

The outskirts of galaxy clusters are continuously disturbed by mergers and gas infall along filaments, which in turn induce turbulent flow motions and shock waves. We examine the properties of shocks that form within $r_{200}$ in sample galaxy clusters from structure formation simulations. While most of these shocks are weak and inefficient accelerators of cosmic rays (CRs), there are a number of strong, energetic shocks which can produce large amounts of CR protons via diffusive shock acceleration. We show that the energetic shocks reside mostly in the outskirts and a substantial fraction of them are induced by infall of the warm-hot intergalactic medium from filaments. As a result, the radial profile of the CR pressure in the intracluster medium is expected to be broad, dropping off more slowly than that of the gas pressure, and might be even temporarily inverted, peaking in the outskirts. The volume-integrated momentum spectrum of CR protons inside $r_{200}$ has the power-law slope of $4.25 - 4.5$, indicating that the average Mach number of the shocks of main CR production is in the range of $\left< M_s \right>_{\rm CR} \approx 3 - 4$. We suggest that some radio relics with relatively flat radio spectrum could be explained by primary electrons accelerated by energetic infall shocks with $M_s >~ 3$ induced in the cluster outskirts.

preprint2013arXiv

A comparison study of CORSIKA and COSMOS simulations for extensive air showers

Cosmic rays with energy exceeding ~ 10^{18} eV are referred to as ultra-high energy cosmic rays (UHECRs). Monte Carlo codes for extensive air shower (EAS) simulate the development of EASs initiated by UHECRs in the Earth's atmosphere. Experiments to detect UHECRs utilize EAS simulations to estimate their energy, arrival direction, and composition. In this paper, we compare EAS simulations with two different codes, CORSIKA and COSMOS, presenting quantities including the longitudinal distribution of particles, depth of shower maximum, kinetic energy distribution of particle at the ground, and energy deposited to the air. We then discuss implications of our results to UHECR experiments.

preprint2013arXiv

Nonthermal Radiation from Supernova Remnants: Effects of Magnetic Field Amplification and Particle Escape

We explore nonlinear effects of wave-particle interactions on the diffusive shock acceleration (DSA) process in Type Ia-like, SNR blast waves, by implementing phenomenological models for magnetic field amplification, Alfv'enic drift, and particle escape in time-dependent numerical simulations of nonlinear DSA. For typical SNR parameters the CR protons can be accelerated to PeV energies only if the region of amplified field ahead of the shock is extensive enough to contain the diffusion lengths of the particles of interest. Even with the help of Alfv'enic drift, it remains somewhat challenging to construct a nonlinear DSA model for SNRs in which order of 10 % of the supernova explosion energy is converted to the CR energy and the magnetic field is amplified by a factor of 10 or so in the shock precursor, while, at the same time, the energy spectrum of PeV protons is steeper than E^{-2}. To explore the influence of these physical effects on observed SNR emissions, we also compute resulting radio-to-gamma-ray spectra. Nonthermal emission spectra, especially in X-ray and gamma-ray bands,depend on the time dependent evolution of CR injection process, magnetic field amplification, and particle escape, as well as the shock dynamic evolution. This result comes from the fact that the high energy end of the CR spectrum is composed of the particles that are injected in the very early stages of blast wave evolution. Thus it is crucial to understand better the plasma wave-particle interactions associated with collisionless shocks in detail modeling of nonthermal radiation from SNRs.

preprint2012arXiv

Diffusive Shock Acceleration at Cosmological Shock Waves

We reexamine nonlinear diffusive shock acceleration (DSA) at cosmological shocks in the large scale structure of the Universe, incorporating wave-particle interactions that are expected to operate in collisionless shocks. Adopting simple phenomenological models for magnetic field amplification (MFA) by cosmic-ray (CR) streaming instabilities and Alfv'enic drift, we perform kinetic DSA simulations for a wide range of sonic and Alfv'enic Mach numbers and evaluate the CR injection fraction and acceleration efficiency. In our DSA model the CR acceleration efficiency is determined mainly by the sonic Mach number Ms, while the MFA factor depends on the Alfv'enic Mach number and the degree of shock modification by CRs. We show that at strong CR modified shocks, if scattering centers drift with an effective Alfv'en speed in the amplified magnetic field, the CR energy spectrum is steepened and the acceleration efficiency is reduced significantly, compared to the cases without such effects. As a result, the postshock CR pressure saturates roughly at ~ 20 % of the shock ram pressure for strong shocks with Ms>~ 10. In the test-particle regime (Ms<~ 3), it is expected that the magnetic field is not amplified and the Alfv'enic drift effects are insignificant, although relevant plasma physical processes at low Mach number shocks remain largely uncertain.

preprint2012arXiv

Diffusive Shock Acceleration Simulations of Radio Relics

Recent radio observations have identified a class of structures, so-called radio relics, in clusters of galaxies. The radio emission from these sources is interpreted as synchrotron radiation from GeV electrons gyrating in microG-level magnetic fields. Radio relics, located mostly in the outskirts of clusters, seem to associate with shock waves, especially those developed during mergers. In fact, they seem to be good structures to identify and probe such shocks in intracluster media (ICMs), provided we understand the electron acceleration and re-acceleration at those shocks. In this paper, we describe time-dependent simulations for diffusive shock acceleration at weak shocks that are expected to be found in ICMs. Freshly injected as well as pre-existing populations of cosmic-ray (CR) electrons are considered, and energy losses via synchrotron and inverse Compton are included. We then compare the synchrotron flux and spectral distributions estimated from the simulations with those in two well-observed radio relics in CIZA J2242.8+5301 and ZwCl0008.8+5215. Considering that the CR electron injection is rather inefficient at weak shocks with Mach number M <~ a few, the existence of radio relics could indicate the pre-existing population of low-energy CR electrons in ICMs. The implication of our results on the merger shock scenario of radio relics is discussed.

preprint2012arXiv

Diffusive shock acceleration with magnetic field amplification and Alfvenic drift

We explore how wave-particle interactions affect diffusive shock acceleration (DSA) at astrophysical shocks by performing time-dependent kinetic simulations, in which phenomenological models for magnetic field amplification (MFA), Alfvenic drift, thermal leakage injection, Bohm-like diffusion, and a free escape boundary are implemented. If the injection fraction of cosmic-ray (CR) particles is greater than 2x10^{-4}, for the shock parameters relevant for young supernova remnants, DSA is efficient enough to develop a significant shock precursor due to CR feedback, and magnetic field can be amplified up to a factor of 20 via CR streaming instability in the upstream region. If scattering centers drift with Alfven speed in the amplified magnetic field, the CR energy spectrum can be steepened significantly and the acceleration efficiency is reduced. Nonlinear DSA with self-consistent MFA and Alfvenic drift predicts that the postshock CR pressure saturates roughly at 10 % of the shock ram pressure for strong shocks with a sonic Mach number ranging 20< M_s< 100. Since the amplified magnetic field follows the flow modification in the precursor, the low energy end of the particle spectrum is softened much more than the high energy end. As a result, the concave curvature in the energy spectra does not disappear entirely even with the help of Alfvenic drift. For shocks with a moderate Alfven Mach number (M_A<10), the accelerated CR spectrum can become as steep as E^{-2.1}-E^{-2.3}, which is more consistent with the observed CR spectrum and gamma-ray photon spectrum of several young supernova remnants.

preprint2011arXiv

Comparison of CORSIKA and COSMOS simulations

Ultra-high-energy cosmic rays (UHECRs) refer to cosmic rays with energy above 10^{18} eV. UHECR experiments utilize simulations of extensive air shower to estimate the properties of UHECRs. The Telescope Array (TA) experiment employs the Monte Carlo codes of CORSIKA and COSMOS to obtain EAS simulations. In this paper, we compare the results of the simulations obtained from CORSIKA and COSMOS and report differences between them in terms of the longitudinal distribution, Xmax-value, calorimetric energy, and energy spectrum at ground.

preprint2011arXiv

Cosmic Ray Spectrum in Supernova Remnant Shocks

We performed kinetic simulations of diffusive shock acceleration in Type Ia supernova remnants (SNRs) expanding into a uniform interstellar medium (ISM). The preshock gas temperature is the primary parameter that governs the cosmic ray (CR) acceleration, while magnetic field strength and CR injection rate are secondary parameters. SNRs in the hot ISM, with an injection fraction smaller than 10^{-4}, are inefficient accelerators with less than 10 % energy getting converted to CRs. The shock structure is almost test-particle like and the ensuing CR spectrum can be steeper than E^{-2}. Although the particles can be accelerated to the knee energy of 10^{15.5}Z eV with amplified magnetic fields in the precursor, Alfv'enic drift of scattering centers softens the source spectrum as steep as E^{-2.1} and reduces the CR acceleration efficiency.

preprint2011arXiv

Nonthermal Radiation from Type Ia Supernova Remnants

We present calculations of expected continuum emissions from Sedov-Taylor phase Type Ia supernova remnants (SNRs), using the energy spectra of cosmic ray (CR) electrons and protons from nonlinear diffusive shock acceleration (DSA) simulations. A new, general-purpose radiative process code, Cosmicp, was employed to calculate the radiation expected from CR electrons and protons and their secondary products. These radio, X-ray and gamma-ray emissions are generally consistent with current observations of Type Ia SNRs. The emissions from electrons in these models dominate the radio through X-ray bands. Decays of π^0 s from p-p collisions mostly dominate the gamma-ray range, although for a hot, low density ISM case (n_{ISM}=0.003 cm^{-3}), the pion decay contribution is reduced sufficiently to reveal the inverse Compton contribution to TeV gamma-rays. In addition, we present simple scalings for the contributing emission processes to allow a crude exploration of model parameter space, enabling these results to be used more broadly. We also discuss the radial surface brightness profiles expected for these model SNRs in the X-ray and gamma-ray bands.

preprint2011arXiv

Re-acceleration of Nonthermal Particles at Weak Cosmological Shock Waves

We examine diffusive shock acceleration (DSA) of the pre-exisiting as well as freshly injected populations of nonthermal, cosmic-ray (CR) particles at weak cosmological shocks. Assuming simple models for thermal leakage injection and Alfvénic drift, we derive analytic, time-dependent solutions for the two populations of CRs accelerated in the test-particle regime. We then compare them with the results from kinetic DSA simulations for shock waves that are expected to form in intracluster media and cluster outskirts in the course of large-scale structure formation. We show that the test-particle solutions provide a good approximation for the pressure and spectrum of CRs accelerated at these weak shocks. Since the injection is extremely inefficient at weak shocks, the pre-existing CR population dominates over the injected population. If the pressure due to pre-existing CR protons is about 5 % of the gas thermal pressure in the upstream flow, the downstream CR pressure can absorb typically a few to 10 % of the shock ram pressure at shocks with the Mach number $M \la 3$. Yet, the re-acceleration of CR electrons can result in a substantial synchrotron emission behind the shock. The enhancement in synchrotron radiation across the shock is estimated to be about a few to several for $M \sim 1.5$ and $10^2-10^3$ for $M \sim 3$, depending on the detail model parameters. The implication of our findings for observed bright radio relics is discussed.

preprint2010arXiv

Comparison of Different Methods for Nonlinear Diffusive Shock Acceleration

We provide a both qualitative and quantitative comparison among different approaches aimed to solve the problem of non-linear diffusive acceleration of particles at shocks. In particular, we show that state-of-the-art models (numerical, Monte Carlo and semi-analytical), even if based on different physical assumptions and implementations, for typical environmental parameters lead to very consistent results in terms of shock hydrodynamics, cosmic ray spectrum and also escaping flux spectrum and anisotropy. Strong points and limits of each approach are also discussed, as a function of the problem one wants to study.

preprint2010arXiv

Diffusive Shock Acceleration in Test-Particle Regime

We examine the test-particle solution for diffusive shock acceleration, based on simple models for thermal leakage injection and Alfv'enic drift. The critical injection rate, ξ_c, above which the cosmic ray (CR) pressure becomes dynamically significant, depends mainly on the sonic shock Mach number, M, and preshock gas temperature, T_1. In the hot-phase interstellar medium (ISM) and intracluster medium, ξ_c < 10^{-3} for shocks with M < 5, while ξ_c ~ 10^{-4}(T_1/10^6 K)^{1/2} for shocks with M > 10. For T_1=10^6 K, for example, the test-particle solution would be valid if the injection momentum, p_{inj} > 3.8 p_{th}. This leads to the postshock CR pressure less than 10% of the shock ram pressure. If the Alfv'en speed is comparable to the sound speed in the preshock flow, as in the hot-phase ISM, the power-law slope of CR spectrum can be significantly softer than the canonical test-particle slope. Then the CR spectrum at the shock can be approximated by the revised test-particle power-law with an exponential cutoff at the highest accelerated momentum, p_{max}(t). An analytic form of the exponential cutoff is also suggested.

preprint2009arXiv

Intergalactic Magnetic Field and Arrival Direction of Ultra-High-Energy Protons

We studied how the intergalactic magnetic field (IGMF) affects the propagation of super-GZK protons that originate from extragalactic sources within the local GZK sphere. Toward this end, we set up hypothetical sources of ultra-high-energy cosmic-rays (UHECRs), virtual observers, and the magnetized cosmic web in a model universe constructed from cosmological structure formation simulations. We then arranged a set of reference objects mimicking active galactic nuclei (AGNs) in the local universe, with which correlations of simulated UHECR events are analyzed. With our model IGMF, the deflection angle between the arrival direction of super-GZK protons and the sky position of their actual sources is quite large with the mean value of $<θ> \sim 15^{\circ}$ and the median value of $\tilde θ\sim 7 - 10^{\circ}$. On the other hand, the separation angle between the arrival direction and the sky position of nearest reference objects is substantially smaller with $<S > \sim 3.5 - 4^{\circ}$, which is similar to the mean angular distance in the sky to nearest neighbors among the reference objects. This is a direct consequence of our model that the sources, observers, reference objects, and the IGMF all trace the matter distribution of the universe. The result implies that extragalactic objects lying closest to the arrival direction of UHECRs are not necessary their actual sources. With our model for the distribution of reference objects, the fraction of super-GZK proton events, whose closest AGNs are true sources, is less than 1/3. We discussed implications of our findings for correlation studies of real UHECR events.

preprint2009arXiv

Self-Similar Evolution of Cosmic-Ray Modified Shocks: The Cosmic-Ray Spectrum

We use kinetic simulations of diffusive shock acceleration (DSA) to study the time-dependent evolution of plane, quasi-parallel, cosmic-ray (CR) modified shocks. Thermal leakage injection of low energy CRs and finite Alfvén wave propagation and dissipation are included. Bohm diffusion as well as the diffusion with the power-law momentum dependence are modeled. As long as the acceleration time scale to relativistic energies is much shorter than the dynamical evolution time scale of the shocks, the precursor and subshock transition approach the time-asymptotic state, which depends on the shock sonic and Alfvénic Mach numbers and the CR injection efficiency. For the diffusion models we employ, the shock precursor structure evolves in an approximately self-similar fashion, depending only on the similarity variable, x/(u_s t). During this self-similar stage, the CR distribution at the subshock maintains a characteristic form as it evolves: the sum of two power-laws with the slopes determined by the subshock and total compression ratios with an exponential cutoff at the highest accelerated momentum, p_{max}(t). Based on the results of the DSA simulations spanning a range of Mach numbers, we suggest functional forms for the shock structure parameters, from which the aforementioned form of CR spectrum can be constructed. These analytic forms may represent approximate solutions to the DSA problem for astrophysical shocks during the self-similar evolutionary stage as well as during the steady-state stage if p_{max} is fixed.

preprint2002arXiv

Acceleration of Cosmic Rays at Large Scale Cosmic Shocks in the Universe

Cosmological hydrodynamic simulations of large scale structure in the universe have shown that accretion shocks and merger shocks form due to flow motions associated with the gravitational collapse of nonlinear structures. Estimated speed and curvature radius of these shocks could be as large as a few 1000 km/s and several Mpc, respectively. According to the diffusive shock acceleration theory, populations of cosmic-ray particles can be injected and accelerated to very high energy by astrophysical shocks in tenuous plasmas. In order to explore the cosmic ray acceleration at the cosmic shocks, we have performed nonlinear numerical simulations of cosmic ray (CR) modified shocks with the newly developed CRASH (Cosmic Ray Amr SHock) numerical code. We adopted the Bohm diffusion model for CRs, based on the hypothesis that strong Alfvén waves are self-generated by streaming CRs. The shock formation simulation includes a plasma-physics-based ``injection'' model that transfers a small proportion of the thermal proton flux through the shock into low energy CRs for acceleration there. We found that, for strong accretion shocks, CRs can absorb most of shock kinetic energy and the accretion shock speed is reduced up to 20 %, compared to pure gas dynamic shocks. For merger shocks with small Mach numbers, however, the energy transfer to CRs is only about 10-20 % with an associated CR particle fraction of $10^{-3}$. Nonlinear feedback due to the CR pressure is insignificant in the latter shocks. Although detailed results depend on models for the particle diffusion and injection, these calculations show that cosmic shocks in large scale structure could provide acceleration sites of extragalactic cosmic rays of the highest energy.

preprint2002arXiv

Clusters of Galaxies: Shock Waves and Cosmic Rays

Recent observations of galaxy clusters in radio and X-ray indicate that cosmic rays and magnetic fields may be energetically important in the intracluster medium. According to the estimates based on theses observational studies, the combined pressure of these two components of the intracluster medium may range between $10% \sim 100 %$ of gas pressure, although their total energy is probably time dependent. Hence, these non-thermal components may have influenced the formation and evolution of cosmic structures, and may provide unique and vital diagnostic information through various radiations emitted via their interactions with surrounding matter and cosmic background photons. We suggest that shock waves associated with cosmic structures, along with individual sources such as active galactic nuclei and radio galaxies, supply the cosmic rays and magnetic fields to the intracluster medium and to surrounding large scale structures. In order to study 1) the properties of cosmic shock waves emerging during the large scale structure formation of the universe, and 2) the dynamical influence of cosmic rays, which were ejected by AGN-like sources into the intracluster medium, on structure formation, we have performed two sets of N-body/hydrodynamic simulations of cosmic structure formation. In this contribution, we report the preliminary results of these simulations.

preprint2002arXiv

Cosmic Rays in Clusters of Galaxies

We argue that clusters of galaxies have an intergalactic medium, which is permeated by strong magnetic fields and also has a contribution of pressure from cosmic rays. These two components of total pressure are probably highly time dependent, and range probably between 1/10 of the gas pressure up to equipartition between gas pressure and the sum of the two other components. Radio galaxies are likely to provide the main source for both magnetic fields and cosmic rays. In this concept it becomes easy to understand the occasional mismatch between the total mass inferred from the assumption of hydrostatic equilibrium derived purely from gas, and the total mass derived from lensing data. We also suggest that the structure and topology of the magnetic field may be highly inhomogeneous - at least over a certain range of scales, and may contain long twisted filaments of strong magnetic fields, as on the Sun. The analogy with the interstellar medium may be fruitful to explore further, where we do not know where magnetic fields come from, but suspect that the cosmic rays derive from supernova explosions. In such an analogy it becomes useful to refer to "radio galaxy explosions" in clusters of galaxies. A full scale exploration of all the implications, especially of the notion that occasionally complete equipartition may be reached, is a task for the future.

preprint1998arXiv

Cosmic Magnetic Fields in Large Scale Filaments and Sheets

We consider the possibility that cosmic magnetic field, instead of being uniformly distributed, is strongly correlated with the large scale structure of the universe. Then, the observed rotational measure of extra-galactic radio sources would be caused mostly by the clumpy magnetic field in cosmological filaments/sheets rather than by a uniform magnetic field, which was often assumed in previous studies. As a model for the inhomogeneity of the cosmological magnetic field, we adopt a cosmological hydrodynamic simulation, where the field is passively included, and can approximately represent the real field distribution with an arbitrary normalization for the field strength. Then, we derive an upper limit of the magnetic field strength by comparing the observed limit of rotational measure with the rotational measure expected from the magnetic field geometry in the simulated model universe. The resulting upper limit to the magnetic field in filaments and sheets is ${\bar B}_{fs} \la 1 μG$ which is $\sim10^3$ times higher than the previously quoted values. This value is close to, but larger than, the equipartition magnetic field strength in filaments and sheets. The amplification mechanism of the magnetic field to the above strength is uncertain. The implications of such a strength of the cosmic magnetic field are discussed.

preprint1996arXiv

Accreting Matter around Clusters of Galaxies: One-Dimensional Considerations

During the formation of the large scale structure of the Universe, matter accretes onto high density peaks. Accreting collisionless dark matter (DM) forms caustics around them, while accreting collisional baryonic matter (BM) forms accretion shocks. The properties of the accreting matter depend upon the power spectrum of the initial perturbations on a given scale as well as the background expansion in a given cosmological model. In this paper, we have calculated the accretion of DM particles in one-dimensional spherical geometry under various cosmological models including the Einstein-de Sitter universe, the open universe with $Ω_o<1$, and the flat universe with $Ω_Λ=1-Ω_o$. A density parameter in the range $0.1\le Ω_o \le 1$ has been considered. The initial perturbation characterized by a point mass at the origin has been considered. Since the accretion shock of BM is expected to form close to the first caustic of DM, the properties of the accreting BM are common with those of the DM. Hence, the accretion calculations with DM particles have been used to find the position and velocity of the accretion shock and the cluster mass inside it. The average temperature of BM has been estimated by adopting simplifying assumptions. The velocity of the accreting BM around clusters of a given temperature is smaller in a universe with smaller $Ω_o$, but only by up to $\sim24\%$ in the models with $0.1\le Ω_o \le 1$. Thus, it would be difficult to use that quantity to discriminate among the cosmological models. However, the accretion velocity around clusters of a given mass or a given radius depends more sensitively on the cosmological models. It is smaller in a universe with smaller $Ω_o$ by up to $\sim41\%$ and $\sim65\%$, respectively. So, it can provide a better

preprint1996arXiv

Contributions to the Cosmic Ray Flux above the Ankle: Clusters of Galaxies

Motivated by the suggestion of Kang, Ryu \& Jones (1996) that particles can be accelerated to high energies via diffusive shock acceleration process at the accretion shocks formed by the infalling flow toward the clusters of galaxies, we have calculated the expected particle flux from a cosmological ensemble of clusters. We use the observed temperature distribution of local clusters and assume a simple power-law evolutionary model for the comoving density of the clusters. The shock parameters such as the shock radius and velocity are deduced from the ICM temperature using the self-similar solutions for secondary infall onto the clusters. The magnetic field strength is assume to be in equipartition with the postshock thermal energy behind the accretion shock. We also assume that the injected energy spectrum is a power-law with the exponential cutoff at the maximum energy which is calculated from the condition that the energy gain rate for diffusive shock acceleration is balanced by the loss rate due to the interactions with the cosmic background radiation. In contrast to the earlier paper we have adopted here the description of the cosmic ray diffusion by Jokipii (1987) which leads to considerably higher particle energies. Finally the injected particle spectrum at the clusters is integrated over the cosmological distance to earth by considering the energy loss due to the interactions with the cosmic background radiation. Our calculations show that the expected spectrum of high-energy protons from the cosmological ensemble of the cluster accretion shocks could match well the observed cosmic ray spectrum near $10^{19}$eV with reasonable parameters and models if about $10^{-4}$ of the infalling kinetic energy can be injected into the intergalactic space as the high energy particles.

preprint1995arXiv

Background X-ray Emission from Hot Gas in CDM and CDM+Lambda Universes: Spectral Signatures

We present a new treatment of two popular models for the growth of structure, examining the X-ray emission from hot gas with allowance for spectral line emission from various atomic species, primarily ``metals". The X-ray emission from the bright cluster sources is not significantly changed from prior work and, as noted earlier, shows the CDM$+Λ$ model (LCDM) to be consistent but the standard, COBE normalized model (SCDM) to be inconsistent with existing observations --- after allowance for still the considerable numerical modelling uncertainties. But we find one important new result. Radiation in the softer band 0.5-1.0keV is predominantly emitted by gas far from cluster centers (hence ``background"). This background emission dominates over the cluster emission below 1keV and observations of it should show clear spectral signatures indicating its origin. In particular the ``iron blend" should be seen prominantly in this spectral bin from cosmic background hot gas at high galactic latitudes and should show shadowing against the SMC indicating its extragalactic origin. Certain OVII lines also provide a signature of this gas which emits a spectrum characteristic of $10^{6.6\pm 0.6}$K gas. Recent ASCA observations of the X-ray background tentatively indicate the presence of component with exactly the spectral features we predict here.

preprint1995arXiv

Diffusive Shock Acceleration Simulations: Comparison with Particle Methods and Bow Shock Measurements

Direct comparisons of diffusive particle acceleration numerical simulations have been made against Monte Carlo and hybrid plasma simulations by Ellison {\it et. al.} (1993) and against observations at the earth's bow shock presented by Ellison {\it et. al.} (1990). Toward this end we have introduced a new numerical scheme for injection of cosmic-ray particles out of the thermal plasma, modeled by way of the diffusive scattering process itself; that is, the diffusion and acceleration across the shock front of particles out of the suprathermal tail of the Maxwellian distribution. We find that all of these computational methods (diffusion-advection, two-fluid, Monte Carlo and hybrid) are in substantial agreement on the issues they can simultaneously address, so that the essential physics of diffusive particle acceleration is adequately contained within each. This is despite the fact that each makes what appear to be very different assumptions or approximations.