Researcher profile

Dongsu Ryu

Dongsu Ryu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 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.

preprint2021arXiv

Exact solution of one dimensional relativistic jet with relativistic equation of state

We study the evolution of one-dimensional relativistic jets, using the exact solution of the Riemann problem for relativistic flows. For this purpose, we solve equations for the ideal special relativistic fluid composed of dissimilar particles in flat space-time and the thermodynamics of fluid is governed by a relativistic equation of state. We obtain the exact solution of jets impinging on denser ambient media. The time variation of the cross-section of the jet-head is modeled and incorporated. We present the initial condition that gives rise to a reverse shock. If the jet-head cross-section increases in time, the jet propagation speed slows down significantly and the reverse-shock may recede opposite to the propagation direction of the jet. We show that the composition of jet and ambient medium can affect the jet solution significantly. For instance, the propagation speed depends on the composition and is maximum for a pair-dominated jet, rather than a pure electron-positron or electron-proton jet. The propagation direction of the reverse-shock may also strongly depend on the composition of the jet.

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.

preprint2011arXiv

Cluster Turbulence: Simulation Insights

Cluster media are dynamical, not static; observational evidence suggests they are turbulent. High-resolution simulations of the intracluster media (ICMs) and of idealized, similar media help us understand the complex physics and astrophysics involved. We present a brief overview of the physics behind ICM turbulence and outline the processes that control its development. High-resolution, compressible, isothermal MHD simulations are used to illustrate important dynamical properties of turbulence that develops in media with initially very weak magnetic fields. The simulations follow the growth of magnetic fields and reproduce the characteristics of turbulence. These results are also compared with full cluster simulations that have examined the properties of ICM turbulence.

preprint2010arXiv

Quasi-Spherical, Time-Dependent Viscous Accretion Flow: One-Dimensional Results

We investigated the instability of advective accretion flow as a consequence of angular momentum transfer in one-dimensional, quasi-spherical transonic accretion flow around a non-rotating black hole. The code is designed to include the effects of viscosity; the hydrodynamics component preserves angular momentum strictly with Lagrangian and remap method in absence of viscosity, while the viscosity component updates viscous angular momentum transfer through the implicit method. We performed two tests to demonstrate the suitability of the code for accretion study. First, we simulated the inviscid, low angular momentum, transonic accretion flow with shocks around a black hole, and then the subsonic, self-similar ADAF solution around a Newtonian object. Both simulations fitted the corresponding analytical curves extremely well. We then simulated a rotating, viscous, transonic fluid with shocks. We showed that for low viscosity parameter, stable shocks at larger distance are possible. For higher viscosity parameter, more efficient angular momentum transfer in the post-shock disk makes the shock structure oscillatory. Moreover, as the shock drifts to larger distances, a secondary inner shock develops. We showed that the inner shock is the direct consequence of expansion of the outer shock, as well as creation of regions with $\partial l / \partial r < 0$ due to more efficient angular momentum transfer near the inner sonic point. We showed that all disk parameters, including emissivity, oscillate with the same period as that of the shock oscillation. Our simulation may have implication for low frequency QPOs, e.g., GRO J1655-40 and XTE J1550-564.

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.

preprint2008arXiv

Effects of Fluid Composition on Spherical Flows around Black Holes

Steady, spherically symmetric, adiabatic accretion and wind flows around non-rotating black holes were studied for fully ionized, multi-component fluids, which are described by a relativistic equation of state (EoS). We showed that the polytropic index depends on the temperature as well as on the composition of fluids, so the composition is important to the solutions of the flows. We demonstrated that fluids with different composition can produce dramatically different solutions, even if they have the same sonic point, or they start with the same specific energy or the same temperature. Then, we pointed that the Coulomb relaxation times can be longer than the dynamical time in the problem considered here, and discussed the implication.

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.