Cosmological Shock Waves in the Large Scale Structure of the Universe: Non-gravitational Effects

Cosmological shock waves result from supersonic flow motions induced by hierarchical clustering of nonlinear structures in the universe. These shocks govern the nature of cosmic plasma through thermalization of gas and acceleration of nonthermal, cosmic-ray (CR) particles. We study the statistics and energetics of shocks formed in cosmological simulations of a concordance $\Lambda$CDM universe, with a special emphasis on the effects of non-gravitational processes such as radiative cooling, photoionization/heating, and galactic superwind feedbacks. Adopting an improved model for gas thermalization and CR acceleration efficiencies based on nonlinear diffusive shock acceleration calculations, we then estimate the gas thermal energy and the CR energy dissipated at shocks through the history of the universe. Since shocks can serve as sites for generation of vorticity, we also examine the vorticity that should have been generated mostly at curved shocks in cosmological simulations. We find that the dynamics and energetics of shocks are governed primarily by the gravity of matter, so other non-gravitational processes do not affect significantly the global energy dissipation and vorticity generation at cosmological shocks. Our results reinforce scenarios in which the intracluster medium and warm-hot intergalactic medium contain energetically significant populations of nonthermal particles and turbulent flow motions.Comment: Submitted to ApJ. Pdf with full resolution figures can be downloaded from http://canopus.cnu.ac.kr/ryu/krco.pd

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$+\Lambda$ 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.Comment: Princeton University Observatory, in ApJ press, figs can be ftp'ed from ftp://astro.princeton.edu/cen/XRAY

A Comparison of Cosmological Hydrodynamic Codes

We present a detailed comparison of the simulation results of various cosmological hydrodynamic codes. Starting with identical initial conditions based on the Cold Dark Matter scenario for the growth of structure, we integrate from redshift $z=20$ to $z=0$ to determine the physical state within a representative volume of size $L^3$ where $L=64 h^{-1} {\rm Mpc}$. Five independent codes are compared: three of them Eulerian mesh based and two variants of the Smooth Particle Hydrodynamics "SPH" Lagrangian approach. The Eulerian codes were run at $N^3=(32^3,~64^3,~128^3,~{\rm and},~256^3)$ cells, the SPH codes at $N^3= 32^3$ and $64^3$ particles. Results were then rebinned to a $16^3$ grid with the expectation that the rebinned data should converge, by all techniques, to a common and correct result as $N \rightarrow \infty$. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainties in even primitive quantities such as $\langle T \rangle$, $\langle \rho^2\rangle^{1/2}$ persists at the 3\%-17\% level after completion of very large simulations. The two SPH codes and the two shock capturing Eulerian codes achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high density, high temperature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by $\rho^2$) is only 4.5\%. Overall the comparison allows us to better estimate errors, it points to ways of improving this current generation of hydrodynamic codes and of suiting their use to problems which exploit their individually best features.Comment: 20p plaintex to appear in The Astrophysical Journal on July 20, 199

Nonthermal Radiation from Type Ia Supernova Remnants

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

Cosmic Ray Electrons in Groups and Clusters of Galaxies: Primary and Secondary Populations from a Numerical Cosmological Simulation

We study the generation and distribution of high energy electrons in cosmic environment and their observational consequences by carrying out the first cosmological simulation that includes directly cosmic ray (CR) particles. Starting from cosmological initial conditions we follow the evolution of primary and secondary electrons (CRE), CR ions (CRI) and a passive magnetic field. CRIs and primary CREs are injected and accelerated at large scale structure shocks. Secondary CREs are continuously generated through inelastic p-p collisions. We include spatial transport, adiabatic expansion/compression, Coulomb collisions, bremsstrahlung, synchrotron (SE)and inverse Compton (IC) emission. We find that, from the perspective of cosmic shock energy and acceleration efficiency, the few detections of hard X-ray radiation excess could be explained in the framework of IC emission of primary CREs in clusters undergoing high accretion/merger phase. Instead, IC emission from both primary and secondary CREs accounts at most for a small fraction of the radiation excesses detected in the extreme-UV (except for the Coma cluster as reported by Bowyer et al.1999). Next, we calculate the SE after normalizing the magnetic field so that for a Coma-like cluster ^1/2~3 \muG. Our results indicate that the SE from secondary CREs reproduces several general properties of radio halos, including the recently found P_1.4GHz vs T relation, the morphology and polarization of the emitting region and, to some extent, the spectral index. Moreover, SE from primary CREs turns out sufficient to power extended regions resembling radio relics observed at the outskirts of clusters. Again we find striking resemblance between morphology, polarization and spectral index of our synthetic maps and those reported in the literature.Comment: emulateapj, 27 pages, 10 figures, 5 tables; ApJ in pres

105110^{51} Ergs: The Evolution of Shell Supernova Remnants

This paper reports on a workshop hosted by the University of Minnesota, March 23-26, 1997. It addressed fundamental dynamical issues associated with the evolution of shell supernova remnants and the relationships between supernova remnants and their environments. The workshop considered, in addition to classical shell SNRs, dynamical issues involving X-ray filled composite remnants and pulsar driven shells, such as that in the Crab Nebula. Approximately 75 participants with wide ranging interests attended the workshop. An even larger community helped through extensive on-line debates prior to the meeting. Each of the several sessions, organized mostly around chronological labels, also addressed some underlying, general physical themes: How are SNR dynamics and structures modified by the character of the CSM and the ISM and vice versa? How are magnetic fields generated in SNRs and how do magnetic fields influence SNRs? Where and how are cosmic-rays (electrons and ions) produced in SNRs and how does their presence influence or reveal SNR dynamics? How does SNR blast energy partition into various components over time and what controls conversion between components? In lieu of a proceedings volume, we present here a synopsis of the workshop in the form of brief summaries of the workshop sessions. The sharpest impressions from the workshop were the crucial and under-appreciated roles that environments have on SNR appearance and dynamics and the critical need for broad-based studies to understand these beautiful, but enigmatic objects. \\Comment: 54 pages text, no figures, Latex (aasms4.sty). submitted to the PAS

Properties of Cosmic Shock Waves in Large Scale Structure Formation

We have examined the properties of shock waves in simulations of large scale structure formation for two cosmological scenarios (a SCDM and a LCDM with Omega =1). Large-scale shocks result from accretion onto sheets, filaments and Galaxy Clusters (GCs) on a scale of circa 5 Mpc/h in both cases. Energetic motions, both residual of past accretion history and due to current asymmetric inflow along filaments, generate additional, common shocks on a scale of about 1 Mpc/h, which penetrate deep inside GCs. Also collisions between substructures inside GCs form merger shocks. Consequently, the topology of the shocks is very complex and highly connected. During cosmic evolution the comoving shock surface density decreases, reflecting the ongoing structure merger process in both scenarios. Accretion shocks have very high Mach numbers (10-10^3), when photo-heating of the pre-shock gas is not included. The typical shock speed is of order v_{sh}(z) =H(z)lambda_{NL}(z), with lambda_{NL}(z) the wavelength scale of the nonlinear perturbation at the given epoch. However, the Mach number for shocks occuring within clusters is usually smaller (3-10), due to the fact that the intracluster gas is already hot. Statistical fits of shock speed around GCs as a function of GCs temperature give power-law's in accord with 1-D predictions. However, a very different result is obtained for fits of the shock radius, reflecting the very complex shock structures forming in 3-D simulations. The in-flowing kinetic energy across such shocks, giving the power available for cosmic-ray acceleration, is comparable to the cluster X-ray luminosity emitted from a central region of radius 0.5 Mpc/h. Considering their large size and long lifetimes, those shocks are potentially interesting sites for cosmic-ray acceleration, if modest magnetic fields exist within them.Comment: 20 Pages, 11 figures, ApJ in press. Complete set of full resolution figures available at http://www.msi.umn.edu:80/Projects/twj/figures.tar.g

Shock Waves in the Large-Scale Structure of the Universe

Cosmological shock waves are induced during hierarchical formation of large-scale structure in the universe. Like most astrophysical shocks, they are collisionless, since they form in the tenuous intergalactic medium through electromagnetic viscosities. The gravitational energy released during structure formation is transferred by these shocks to the intergalactic gas as heat, cosmic-rays, turbulence, and magnetic fields. Here we briefly describe the properties and consequences of the shock waves in the context of the large-scale structure of the universe.Comment: Submitted to Astrophysics and Space Science (Special Issue for the proceedings of International Conference on HEDP/HEDLA-08). Pdf with full resolution Figure 1 can be downloaded from http://canopus.cnu.ac.kr/ryu/rk.pd

Cosmic Ray Protons Accelerated at Cosmological Shocks and Their Impact on Groups and Clusters of Galaxies

We investigate the production of cosmic ray (CR) protons at cosmological shocks by performing, for the first time, numerical simulations of large scale structure formation that include directly the acceleration, transport and energy losses of the high energy particles. CRs are injected at shocks according to the thermal leakage model and, thereafter, accelerated to a power-law distribution as indicated by the test particle limit of the diffusive shock acceleration theory. The evolution of the CR protons accounts for losses due to adiabatic expansion/compression, Coulomb collisions and inelastic p-p scattering. Our results suggest that CR protons produced at shocks formed in association with the process of large scale structure formation could amount to a substantial fraction of the total pressure in the intra-cluster medium. Their presence should be easily revealed by GLAST through detection of gamma-ray flux from the decay of neutral pions produced in inelastic p-p collisions of such CR protons with nuclei of the intra-cluster gas. This measurement will allow a direct determination of the CR pressure contribution in the intra-cluster medium. We also find that the spatial distribution of CR is typically more irregular than that of the thermal gas because it is more influenced by the underlying distribution of shocks. This feature is reflected in the appearance of our gamma-ray synthetic images. Finally, the average CR pressure distribution appears statistically slightly more extended than the thermal pressure.Comment: 16 pages, 7 figures, ApJ in pres