Source author record

Dongsu Ryu

Dongsu Ryu 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

51works
7topics
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

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

preprint2016arXiv

Fast Radio Bursts as Probes of Magnetic Fields in the Intergalactic Medium

We examine the proposal that the dispersion measures (DMs) and Faraday rotation measures (RMs) of extragalactic linearly-polarized fast radio bursts (FRBs) can be used to probe the intergalactic magnetic field (IGMF) in filaments of galaxies. The DM through the cosmic web is dominated by contributions from the warm-hot intergalactic medium (WHIM) in filaments and from the gas in voids. On the other hand, RM is induced mostly by the hot medium in galaxy clusters, and only a fraction of it is produced in the WHIM. We show that if one excludes FRBs whose sightlines pass through galaxy clusters, the line-of-sight strength of the IGMF in filaments, $B_{||}$, is approximately $C(\langle 1+z \rangle/f_{DM})(RM/DM)$, where $C$ is a known constant. Here, the redshift of the FRB is not required to be known; $f_{DM}$ is the fraction of total DM due to the WHIM, while $\langle 1+z \rangle$ is the redshift of interevening gas weighted by the WHIM gas density, both of which can be evaluated for a given cosmology model solely from the DM of an FRB. Using data on structure formation simulations and a model IGMF, we show that $C(\langle 1+z \rangle/f_{DM})(RM/DM)$ closely reproduces the density-weighted line-of-sight strength of the IGMF in filaments of the large-scale structure.

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

SKA Deep Polarization and Cosmic Magnetism

Deep surveys with the SKA1-MID array offer for the first time the opportunity to systematically explore the polarization properties of the microJy source population. Our knowledge of the polarized sky approaching these levels is still very limited. In total intensity the population will be dominated by star-forming and normal galaxies to intermediate redshifts ($z \sim1-2$), and low-luminosity AGN to high redshift. The polarized emission from these objects is a powerful probe of their intrinsic magnetic fields and of their magnetic environments. For redshift of order 1 and above the broad bandwidth of the mid-bands span the Faraday thick and thin regimes allowing study of the intrinsic polarization properties of these objects as well as depolarization from embedded and foreground plasmas. The deep field polarization images will provide Rotation Measures data with very high solid angle density allowing a sensitive statistical analysis of the angular variation of RM on critical arc-minute scales from a magnetic component of Large Scale Structure of the Universe.

preprint2015arXiv

Vorticity, Shocks and Magnetic Fields in Subsonic, ICM-like Turbulence

We analyze high resolution simulations of compressible, MHD turbulence with properties resembling conditions in galaxy clusters. The flow is driven to turbulence Mach number $\mathcal{M}_t \sim 1/2$ in an isothermal medium with an initially very weak, uniform seed magnetic field ($β= P_g/P_B = 10^6$). Since cluster turbulence is likely to result from a mix of sheared (solenoidal) and compressive forcing processes, we examine the distinct turbulence properties for both cases. In one set of simulations velocity forcing is entirely solenoidal ($\nabla\cdot δ{\vec u} = 0$), while in the other it is entirely compressive ($\nabla\times δ{\vec u} = 0$). Both cases develop a mixture of solenoidal and compressive turbulent motions, since each generates the other. The development of compressive turbulent motions leads to shocks, even when the turbulence is solenoidally forced and subsonic. Shocks, in turn, produce and amplify vorticity, which is especially important in compressively forced turbulence. To clarify those processes we include a pair of appendices that look in detail at vorticity evolution in association with shocks. From our simulation analyses we find that magnetic fields amplified to near saturation levels in predominantly solenoidal turbulence can actually enhance vorticity on small scales by concentrating and stabilizing shear. The properties, evolution rates and relative contributions of the kinetic and magnetic turbulent elements depend strongly on the character of the forcing. Specifically, shocks are stronger, but vorticity evolution and magnetic field amplification are slower and weaker when the turbulence is compressively forced. We identify a simple relation to estimate characteristic shock strengths in terms of the turbulence Mach number and the character of the forcing. Our results will be helpful in understanding flow motions in galaxy clusters.

preprint2014arXiv

Exploring the Intergalactic Magnetic Field by Means of Faraday Tomography

Unveiling the intergalactic magnetic field (IGMF) in filaments of galaxies is a very important and challenging subject in modern astronomy. In order to probe the IGMF from rotation measures (RMs) of extragalactic radio sources, we need to separate RMs due to other origins such as the source, intervening galaxies, and our Galaxy. In this paper, we discuss observational strategies for the separation by means of Faraday tomography (Faraday RM Synthesis). We consider an observation of a single radio source such as a radio galaxy or a quasar viewed through the Galaxy and the cosmic web. We then compare the observation with another observation of a neighbor source with a small angular separation. Our simulations with simple models of the sources suggest that it would be not easy to detect the RM due to the IGMF of order ~ 1 rad/m/m, an expected value for the IGMF through a single filament. Contrary to it, we find that the RM of at least ~10 rad/m/m could be detected with the SKA or its pathfinders/precursors, if we achieve selections of ideal sources. These results would be improved if we incorporate decomposition techniques such as RMCLEAN and QU-fitting. We discuss feasibility of the strategies for cases with complex Galactic emissions as well as with effects of observational noise and radio frequency interferences.

preprint2014arXiv

Faraday dispersion functions of galaxies

The Faraday dispersion function (FDF), which can be derived from an observed polarization spec- trum by Faraday rotation measure synthesis, is a profile of polarized emissions as a function of Faraday depth. We study intrinsic FDFs along sight lines through face-on, Milky-Way-like galaxies by means of a sophisticated galactic model incorporating 3D MHD turbulence, and investigate how much the FDF contains information intrinsically. Since the FDF reflects distributions of thermal and cosmic- ray electrons as well as magnetic fields, it has been expected that the FDF could be a new probe to examine internal structures of galaxies. We, however, find that an intrinsic FDF along a sight line through a galaxy is very complicated, depending significantly on actual configurations of turbulence. We perform 800 realizations of turbulence, and find no universal shape of the FDF even if we fix the global parameters of the model. We calculate the probability distribution functions of the standard deviation, skewness, and kurtosis of FDFs and compare them for models with different global pa- rameters. Our models predict that the presence of vertical magnetic fields and large scale-height of cosmic-ray electrons tend to make the standard deviation relatively large. Contrastingly, differences in skewness and kurtosis are relatively less significant.

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

Magnetic Field and Faraday Rotation Measure in the Turbulent Warm Ionized Medium

Wu et al. (2009) suggested an empirical relation between the magnetic field strength along the line of sight (LOS) and the dispersion of Faraday rotation measure (RM) distribution in turbulent media with root-mean-square sonic Mach number $M_s \simeq 1$. In this paper, we extend the work by incorporating the Mach number dependence. Media with $\sim 0.5 < M_s < \sim 2$ are considered to cover the Mach number range of the warm ionized medium (WIM) of our Galaxy. Three-dimensional, magnetohydrodynamic isothermal turbulence simulations with solenoidal forcing are used. We suggest a new relation among the LOS magnetic field strength, the dispersion of RM distribution, and the Mach number, which approximately represents the relation for Alfvénic perturbations. In addition, we suggest a relation between the Mach number and the dispersion of log-normal distribution of emission measure (EM), which is basically the relation for the Mach number and the density dispersion. The relations could be used for a quick and rough estimation of the LOS magnetic field strength in the turbulent WIM.

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.

preprint2014arXiv

Statistical Techniques for Detecting the Intergalactic Magnetic Field from Large Samples of Extragalactic Faraday Rotation Data

Rotation measure (RM) grids of extragalactic radio sources have been widely used for studying cosmic magnetism. But their potential for exploring the intergalactic magnetic field (IGMF) in filaments of galaxies is unclear, since other Faraday-rotation media such as the radio source itself, intervening galaxies, and the interstellar medium of our Galaxy are all significant contributors. We study statistical techniques for discriminating the Faraday rotation of filaments from other sources of Faraday rotation in future large-scale surveys of radio polarization. We consider a 30 degree times 30 degree field-of-view toward the south Galactic pole, while varying the number of sources detected in both present and future observations. We select sources located at high redshifts and toward which depolarization and optical absorption systems are not observed, so as to reduce the RM contributions from the sources and intervening galaxies. It is found that a high-pass filter can satisfactorily reduce the RM contribution from the Galaxy, since the angular scale of this component toward high Galactic latitudes would be much larger than that expected for the IGMF. Present observations do not yet provide a sufficient source density to be able to estimate the RM of filaments. However, from the proposed approach with forthcoming surveys, we predict significant residuals of RM that should be ascribable to filaments. The predicted structure of the IGMF down to scales of 0.1 degree should be observable with data from the SKA, if we achieve selections of sources toward which sight lines do not contain intervening galaxies and RM errors are less than a few rad m-2.

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

Fisher Analysis on Wide-Band Polarimetry for Probing the Intergalactic Magnetic Field

We investigate the capability of ongoing radio telescopes for probing Faraday rotation measure (RM) due to the intergalactic magnetic field (IGMF) in the large-scale structure of the universe which is expected to be of order $O(1) {\rm rad/m^2}$. We consider polarization observations of a compact radio source such as quasars behind a diffuse source such as the Galaxy, and calculate Stokes parameters $Q$ and $U$ assuming a simple model of the Faraday dispersion functions with Gaussian shape. Then, we perform the Fisher analysis to estimate the expected errors in the model parameters from QU-fitting of polarization intensity, accounting for sensitivities and frequency bands of Australian Square Kilometer Array Pathfinder, Low Frequency Array, and the Giant Meterwave Radio Telescope. Finally, we examine the condition on the source intensities which are required to detect the IGMF. Our analysis indicates that the QU-fitting is promising for forthcoming wideband polarimetry to explore RM due to the IGMF in filaments of galaxies.

preprint2013arXiv

Simulated Faraday Rotation Measures toward High Galactic Latitudes

We study the Faraday rotation measure (RM) due to the Galactic magnetic field (GMF) toward high Galactic latitudes. The RM arises from the global, regular component as well as from the turbulent, random component of the GMF. We model the former based on observations and the latter using the data of magnetohydrodynamic turbulence simulations. For a large number of different GMF models, we produce mock RM maps around the Galactic poles and calculate various statistical quantities with the RM maps. We find that the observed medians of RMs toward the north and south Galactic poles, $\sim 0.0\pm 0.5 {\rm rad m^{-2}}$ and $\sim +6.3\pm 0.5 {\rm rad m^{-2}}$, are difficult to explain with any of our many alternate GMF models. The standard deviation of observed RMs, $\sim 9 {\rm rad m^{-2}}$, is clearly larger than that of simulated RMs. The second-order structure function of observed RMs is substantially larger than that of simulated RMs, especially at small angular scales. We discuss other possible contributions to RM toward high Galactic latitudes. Besides observational errors and the intrinsic RM of background radio sources against which RM is observed, we suggest that the RM due to the intergalactic magnetic field may account for a substantial fraction of the observed RM. Finally we note that reproducing the observed medians may require additional components or/and structures of the GMF that are not present in our models.

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

Simulation of radiation driven wind from disc galaxies

We present 2-D hydrodynamic simulation of rotating galactic winds driven by radiation. We study the structure and dynamics of the cool and/or warm component($T \simeq 10^4$ K) which is mixed with dust. We have taken into account the total gravity of a galactic system that consists of a disc, a bulge and a dark matter halo. We find that the combined effect of gravity and radiation pressure from a realistic disc drives the gas away to a distance of $\sim 5$ kpc in $\sim 37$ Myr for typical galactic parameters. The outflow speed increases rapidly with the disc Eddington parameter $Γ_0(=κI/(2 c G Σ)$) for $Γ_0 \ge 1.5$. We find that the rotation speed of the outflowing gas is $\lesssim 100$ km s$^{-1}$. The wind is confined in a cone which mostly consist of low angular momentum gas lifted from the central region.

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.

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

Faraday Rotation Measure due to the Intergalactic Magnetic Field II: the Cosmological Contribution

We investigate the Faraday rotation measure (RM) due to the intergalactic magnetic field (IGMF) through the cosmic web up to cosmological distances, using a model IGMF based on turbulence dynamo in the large-scale structure of the universe. By stacking the IGMF and gas density data up to redshift $z=5$ and taking account of the redshift distribution of polarized background radio sources against which the RM is measured, we simulate the sky map of the RM. The contribution from galaxy clusters is subtracted from the map, based on several different criteria of X-ray brightness and temperature. Our findings are as follows. The distribution of RM for radio sources of different redshifts shows that the root-mean-square (rms) value increases with redshift and saturates for $z \ga 1$. The saturated value is RM$_{\rm rms} \approx$ several ${\rm rad m^{-2}}$. The probability distribution function of $|{\rm RM}|$ follows the lognormal distribution. The power spectrum has a broad plateau over the angular scale of $\sim 1 - 0.1^\circ$ with a peak around $\sim 0.15^\circ$. The second-order structure function has a flat profile in the angular separation of $\ga 0.2^\circ$. Our results could provide useful insights for surveys to explore the IGMF with the Square Kilometer Array (SKA) and upcoming SKA pathfinders.

preprint2011arXiv

Magnetic fields in the Large-Scale Structure of the Universe

Magnetic fields appear to be ubiquitous in astrophysical environments. Their existence in the intracluster medium is established through observations of synchrotron emission and Faraday rotation. On the other hand, the nature of magnetic fields outside of clusters, where observations are scarce and controversial, remains largely unknown. In this chapter, we review recent developments in our understanding of the nature and origin of intergalactic magnetic fields, and in particular, intercluster fields. A plausible scenario for the origin of galactic and intergalactic magnetic fields is for seed fields, created in the early universe, to be amplified by turbulent flows induced during the formation of the large scale structure. We present several mechanisms for the generation of seed fields both before and after recombination. We then discuss the evolution and role of magnetic fields during the formation of the first starts. We describe the turbulent amplification of seed fields during the formation of large scale structure and the nature of the magnetic fields that arise. Finally, we discuss implications of intergalactic magnetic fields.

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.

preprint2011arXiv

The First Magnetic Fields

We review current ideas on the origin of galactic and extragalactic magnetic fields. We begin by summarizing observations of magnetic fields at cosmological redshifts and on cosmological scales. These observations translate into constraints on the strength and scale magnetic fields must have during the early stages of galaxy formation in order to seed the galactic dynamo. We examine mechanisms for the generation of magnetic fields that operate prior during inflation and during subsequent phase transitions such as electroweak symmetry breaking and the quark-hadron phase transition. The implications of strong primordial magnetic fields for the reionization epoch as well as the first generation of stars is discussed in detail. The exotic, early-Universe mechanisms are contrasted with astrophysical processes that generate fields after recombination. For example, a Biermann-type battery can operate in a proto-galaxy during the early stages of structure formation. Moreover, magnetic fields in either an early generation of stars or active galactic nuclei can be dispersed into the intergalactic medium.

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.

preprint2010arXiv

Faraday Rotation Measure due to the Intergalactic Magnetic Field

Studying the nature and origin of the intergalactic magnetic field (IGMF) is an outstanding problem of cosmology. Measuring Faraday rotation would be a promising method to explore the IGMF in the large-scale structure (LSS) of the universe. We investigated the Faraday rotation measure (RM) due to the IGMF in filaments of galaxies using simulations for cosmological structure formation. We employed a model IGMF based on turbulence dynamo in the LSS of the universe; it has an average strength of $< B > \sim 10$ nG and a coherence length of several $\times\ 100\ h^{-1}$ kpc in filaments. With the coherence length smaller than path length, the inducement of RM would be a random walk process, and we found that the resultant RM is dominantly contributed by the density peak along line of sight. The rms of RM through filaments at the present universe was predicted to be $\sim 1\ {\rm rad\ m^{-2}}$. In addition, we predicted that the probability distribution function of $|{\rm RM}|$ through filaments follows the log-normal distribution, and the power spectrum of RM in the local universe peaks at a scale of $\sim 1\ h^{-1}$ Mpc. Our prediction of RM could be tested with future instruments.

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

Characteristic Lengths of Magnetic Field in Magnetohydrodynamic Turbulence

In the framework of turbulence dynamo, flow motions amplify a weak seed magnetic field through the stretching of field lines. Although the amplification process has been a topic of active research, less attention has been paid to the length scales of magnetic field. In this paper, we described a numerical study on characteristic lengths of magnetic field in magnetohydrodynamic turbulence. We considered the case of very weak or zero mean magnetic field, which is applicable to the turbulence in the intergalactic space. Our findings are as follows. (1) At saturation, the peak of magnetic field spectrum occurs at $\sim L_0/2$, where $L_0$ is the energy injection scale, while the most energy containing scale is $\sim L_0/5$. The peak scale of spectrum of projected, two-dimensional field is $\sim L_0$. (2) During the stage of magnetic field amplification, the energy equipartition scale shows a power-law increase of $\sim t^{1.5}$, while the integral and curvature scales show a linear increase. The equipartition, integral, and curvature scales saturate at $\sim L_0$, $\sim 0.3L_0$, and $\sim 0.15L_0$, respectively. (3) The coherence length of magnetic field defined in the Faraday rotation measure (RM) due to the intergalactic magnetic field (IGMF) is related to the integral scale. We presented a formula that expresses the standard deviation of RM, $σ_{RM}$, in terms of the integral scale and rms strength of the IGMF, and estimated that $σ_{RM}$ would be $\sim 100$ and $\sim$ a few rad m$^{-2}$ for clusters and filaments, respectively.

preprint2009arXiv

Estimation of Magnetic Field Strength in the Turbulent Warm Ionized Medium

We studied Faraday rotation measure (RM) in turbulent media with the rms Mach number of unity, using isothermal, magnetohydrodynamic turbulence simulations. Four cases with different values of initial plasma beta were considered. Our main findings are as follows. (1) There is no strong correlation between the fluctuations of magnetic field strength and gas density. So the magnetic field strength estimated with RM/DM (DM is the dispersion measure) correctly represents the true mean strength of the magnetic field along the line of sight. (2) The frequency distribution of RMs is well fitted to the Gaussian. In addition, there is a good relation between the width of the distribution of RM/$\bar{\rm RM}$ ($\bar{\rm RM}$ is the average value of RMs) and the strength of the regular field along the line of sight; the width is narrower, if the field strength is stronger. We discussed the implications of our findings in the warm ionized medium where the Mach number of turbulent motions is around unity.

preprint2009arXiv

Hydrodynamic Simulations of Oscillating Shock Waves in a Sub-Keplerian Accretion Flow Around Black Holes

We study the accretion processes on a black hole by numerical simulation. We use a grid based finite difference code for this purpose. We scan the parameter space spanned by the specific energy and the angular momentum and compare the time-dependent solutions with those obtained from theoretical considerations. We found several important results (a) The time dependent flow behaves close to a constant height model flow in the pre-shock region and a flow with vertical equilibrium in the post-shock region. (c) The infall time scale in the post-shock region is several times higher than the free-fall time scale. (b) There are two discontinuities in the flow, one being just outside of the inner sonic point. Turbulence plays a major role in determining the locations of these discontinuities. (d) The two discontinuities oscillate with two different frequencies and behave as a coupled harmonic oscillator. A Fourier analysis of the variation of the outer shock location indicates higher power at the lower frequency and lower power at the higher frequency. The opposite is true when the analysis of the inner shock is made. These behaviours will have implications in the spectral and timing properties of black hole candidates.

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.

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.

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.

preprint1999arXiv

On the Exchange of Kinetic and Magnetic Energy Between Clouds and the Interstellar Medium

We investigate, through 2D MHD numerical simulations, the interaction of a uniform magnetic field oblique to a moving interstellar cloud. In particular we explore the transformation of cloud kinetic energy into magnetic energy as a result of field line stretching. Some previous simulations have emphasized the possible dynamical importance of a ``magnetic shield'' formed around clouds when the magnetic field is perpendicular to the cloud motion (Jones et al. 1996, Miniati et al. 1998). It was not clear, however, how dependent those findings were to the assumed field configuration and cloud properties. To expand our understanding of this effect, we examine several new cases by varing the magnetic field orientation angle with respect to the cloud motion (θ), the cloud-background density contrast, and the cloud Mach number. We show that in 2D and with θlarge enough, the magnetic field tension can become dominant in the dynamics of the motion of high density contrast, low Mach number clouds. In such cases a significant fraction of cloud kinetic energy can be transformed into magnetic energy with the magnetic pressure at the cloud nose exceeding the ram pressure of the impinging flow. We derive a characteristic timescale for this process of energy ``conversion''. We find also that unless the cloud motion is highly aligned to the magnetic field, reconnection through tearing mode instabilities in the cloud wake limit the formation of a strong flux rope feature following the cloud. Finally we attempt to interpret some observational properties of the magnetic field in view of our results.

preprint1999arXiv

The MHD Kelvin-Helmholtz Instability III: The Role of Sheared Magnetic Field in Planar Flows

We have carried out simulations of the nonlinear evolution of the magnetohydrodynamic (MHD) Kelvin-Helmholtz (KH) instability for compressible fluids in $2\frac{1}{2}$-dimensions, extending our previous work by Frank et al (1996) and Jones \etal (1997). In the present work we have simulated flows in the x-y plane in which a ``sheared'' magnetic field of uniform strength ``smoothly'' rotates across a thin velocity shear layer from the z direction to the x direction, aligned with the flow field. We focus on dynamical evolution of fluid features, kinetic energy dissipation, and mixing of the fluid between the two layers, considering their dependence on magnetic field strength for this geometry. The introduction of magnetic shear can allow a Cat's Eye-like vortex to form, even when the field is stronger than the nominal linear instability limit given above. For strong fields that vortex is asymmetric with respect to the preliminary shear layer, however, so the subsequent dissipation is enhanced over the uniform field cases of comparable field strength. In fact, so long as the magnetic field achieves some level of dynamical importance during an eddy turnover time, the asymmetries introduced through the magnetic shear will increase flow complexity, and, with that, dissipation and mixing. The degree of the fluid mixing between the two layers is strongly influenced by the magnetic field strength. Mixing of the fluid is most effective when the vortex is disrupted by magnetic tension during transient reconnection, through local chaotic behavior that follows.

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.

preprint1998arXiv

Hydrodynamics of Cloud Collisions in 2D: The Fate of Clouds in a Multi-phase Medium

We have studied head-on collisions between equal-mass, mildly supersonic (Mach number 1.5) HI clouds, in a standard Two-phase ISM (T_cl = 74 K, n_cl = 22 cm^-3, χ= 100). We explore the role of various factors, including the radiative cooling parameter η= t_rad/t_coll (t_coll=R_c/v_c), evolutionary modifications on the cloud structure (by colliding clouds ``evolved'' through independent motion within the intercloud medium (WIM)), and the symmetry of the problem (by colliding initially identical clouds, evolved to different ages before impact). The presence of bow shocks and ram pressure from material in the cloud wake, developed during such evolution through the WIM, significantly alters these interactions with respect to the standard case of non-evolved clouds. In general, in our adiabatic collisions the clouds are disrupted and convert their gas into a few low density contrast clumps. By contrast, for symmetric radiative cases we find that the two clouds coalesce, with almost all the initial kinetic energy radiated away. On the other hand, for both adiabatic and radiative collisions, asymmetric collisions have a much greater tendency to disrupt the two clouds. Fragmentation of the clouds may occur, and instabilities are in general enhanced. In addition, radiative cooling is less efficient in our asymmetric interactions, so that those parts of the clouds that initially seem to merge are more likely to re-expand and fade into the WIM. Since the majority of real cloud collisions should be asymmetric for one reason or another, we conclude that most gasdynamical diffuse cloud collisions will be disruptive, at least in the absence of significant self-gravity or of a significant magnetic field.

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

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

Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for Multi-Dimensional Flow

We present for astrophysical use a multi-dimensional numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference method on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. Multiple spatial dimensions are treated through a Strang-type operator splitting. The constraint of a divergence-free field is enforced exactly by calculating a correction via a gauge transformation in each time step. Results from two-dimensional shock tube tests show that the code captures correctly discontinuities in all three MHD waves families as well as contact discontinuities. The numerical viscosities and resistivity in the code, which are useful in order to understand simulations involving turbulent flows, are estimated through the decay of two-dimensional linear waves. Finally, the robustness of the code in two-dimensions is demonstrated through calculations of the Kelvin-Helmholtz instability and the Orszag-Tang vortex.