Using wide-field $^{13}$CO ($J=1-0$) data taken with the Nobeyama 45-m telescope, we investigate cloud structures of the infrared dark cloud complex in M17 with SCIMES. In total, we identified 118 clouds that contain 11 large clouds with radii larger than 1 pc. The clouds are mainly distributed in the two representative velocity ranges of 10 $-$ 20 km s$^{-1}$ and 30 $-$ 40 km s$^{-1}$. By comparing with the ATLASGAL catalog, we found that the majority of the $^{13}$CO clouds with 10 $-$ 20 km s$^{-1}$ and 30 $-$ 40 km s$^{-1}$ are likely located at distances of 2 kpc (Sagittarius arm) and 3 kpc (Scutum arm), respectively. Analyzing the spatial configuration of the identified clouds and their velocity structures, we attempt to reveal the origin of the cloud structure in this region. Here we discuss three possibilities: (1) overlapping with different velocities, (2) cloud oscillation, and (3) cloud-cloud collision. From the position-velocity diagrams, we found spatially-extended faint emission between $\sim$ 20 km s$^{-1}$ and $\sim$ 35 km s$^{-1}$, which is mainly distributed in the spatially-overlapped areas of the clouds. We also found that in some areas where clouds with different velocities overlapped, the magnetic field orientation changes abruptly. The distribution of the diffuse emission in the position-position-velocity space and the bending magnetic fields appear to favor the cloud-cloud collision scenario compared to other scenarios. In the cloud-cloud collision scenario, we propose that two $\sim$35 km s$^{-1}$ foreground clouds are colliding with clouds at $\sim$20 km s$^{-1}$ with a relative velocity of 15 km s$^{-1}$. These clouds may be substructures of two larger clouds having velocities of $\sim$ 35 km s$^{-1}$ ($\gtrsim 10^3 $ M$_{\odot}$) and $\sim$ 20 km s$^{-1}$ ($\gtrsim 10^4 $ M$_{\odot}$), respectively.

Observations of large-scale radio emissions prove the existence of shock accelerated cosmic-ray electrons in galaxy clusters, while the lack of detected $\gamma$-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.

We propose semi-analytic models for the electron momentum distribution in weak shocks that accounts for both in situ acceleration and re-acceleration through diffusive shock acceleration (DSA). In the former case, a small fraction of incoming electrons are assumed to be reflected at the shock ramp and pre-accelerated to the so-called injection momentum, $p_{\rm inj}$, above which particles can diffuse across the shock transition and participate in the DSA process. This leads to the DSA power-law distribution extending from the smallest momentum of reflected electrons, $p_{\rm ref}$, all the way to the cutoff momentum, $p_{\rm eq}$, constrained by radiative cooling. In the later case, fossil electrons, specified by a power-law spectrum with a cutoff, are assumed to be re-accelerated also from $p_{\rm ref}$ up to $p_{\rm eq}$ via DSA. We then show that, in the in situ acceleration model, the amplitude of radio synchrotron emission depends strongly on the shock Mach number, whereas it varies rather weakly in the re-acceleration model. Considering rather turbulent nature of shocks in the intracluster medium, such extreme dependence for the in situ acceleration might not be compatible with relatively smooth surface brightness of observed radio relics.

The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

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.

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 $\gamma$-ray and neutrino emissions, resulting from inelastic collisions between CRp and thermal protons. The predicted $\gamma$-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_{\nu}=1$ PeV and $\lesssim10^{-6}$ of the atmospheric neutrino flux in the energy range of $E_{\nu}\leq1$ TeV.

The non-uniform distribution of gas and protostars in molecular clouds is caused by combinations of various physical processes that are difficult to separate. We explore this non-uniform distribution in the M17 molecular cloud complex that hosts massive star formation activity using the $^{12}$CO ($J=1-0$) and $^{13}$CO ($J=1-0$) emission lines obtained with the Nobeyama 45m telescope. Differences in clump properties such as mass, size, and gravitational boundedness reflect the different evolutionary stages of the M17-H{\scriptsize II} and M17-IRDC clouds. Clumps in the M17-H{\scriptsize II} cloud are denser, more compact, and more gravitationally bound than those in M17-IRDC. While M17-H{\scriptsize II} hosts a large fraction of very dense gas (27\%) that has column density larger than the threshold of $\sim$ 1 g cm$^{-2}$ theoretically predicted for massive star formation, this very dense gas is deficient in M17-IRDC (0.46\%). Our HCO$^+$ ($J=1-0$) and HCN ($J=1-0$) observations with the TRAO 14m telescope, {\nlqb trace all gas with column density higher than $3\times 10^{22}$ cm$^{-2}$}, confirm the deficiency of high density ($\gtrsim 10^5$ cm$^{-3}$) gas in M17-IRDC. Although M17-IRDC is massive enough to potentially form massive stars, its deficiency of very dense gas and gravitationally bound clumps can explain the current lack of massive star formation.

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-$\beta$ ($\beta=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 $\beta\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, $\gamma$, 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 ($\omega_r=0$), those excited by the EBFI have a non-zero frequency ($\omega_r\neq0$). However, the frequency is small with $\omega_r<\gamma$. 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.

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.

A high concentration of ultra-high-energy cosmic ray (UHECR) events, called a hotspot, was reported by the Telescope Array (TA) experiment, but its origin still remains unsolved. One of the obstacles is that there is no astronomical object, which could be the source, behind the TA hotpot. In an effort to understand the origin of the TA hotspot, we suggested a model based on the magnetized cosmic web structure. The UHECRs were produced from sources in the Virgo cluster and were initially confined by cluster magnetic fields for a certain period. Next, some of them preferentially escaped to and propagated along filaments. Eventually, they were scattered by filament magnetic fields, and come to us. To examine the model, we followed the propagation trajectories of UHE protons in a simulated universe with clusters, filaments, and voids, by employing a number of models for cosmic magnetic fields. In this study, we present some of the initial results, such as the ratio between the particles directly escaping from the clusters to the voids and particles escaping from the clusters to the filaments. We also discuss the feasibility of our model for the origin of the hotspot by examining the trajectories of the UHE protons.

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.

Low sonic Mach number shocks form in the intracluster medium (ICM) during the formation of the large-scale structure of the universe. Nonthermal cosmic-ray (CR) protons are expected to be accelerated via diffusive shock acceleration (DSA) in those ICM shocks, although observational evidence for the $\gamma$-ray emission of hadronic origin from galaxy clusters has yet to be established. Considering the results obtained from recent plasma simulations, we improve the analytic test-particle DSA model for weak quasi-parallel ($Q_\parallel$) shocks, previously suggested by \citet{kang2010}. In the model CR spectrum, the transition from the postshock thermal to CR populations occurs at the injection momentum, $p_{\rm inj}$, above which protons can undergo the full DSA process. As the shock energy is transferred to CR protons, the postshock gas temperature should decrease accordingly and the subshock strength weakens due to the dynamical feed of the CR pressure to the shock structure. This results in the reduction of the injection fraction, although the postshock CR pressure approaches an asymptotic value when the CR spectrum extends to the relativistic regime. Our new DSA model self-consistently accounts for such behaviors and adopts better estimations for $p_{\rm inj}$. With our model DSA spectrum, the CR acceleration efficiency ranges $\eta\sim10^{-3}-0.01$ for supercritical, $Q_\parallel$-shocks with sonic Mach number $2.25\lesssim M_{\rm s}\lesssim5$ in the ICM. Based on \citet{ha2018b}, on the other hand, we argue that proton acceleration would be negligible in subcritical shocks with $M_{\rm s}<2.25$.

"TRAO FUNS" is a project to survey Gould Belt's clouds in molecular lines. This paper presents its first results on the central region of the California molecular cloud, L1478. We performed On-The-Fly mapping observations using the Taedeok Radio Astronomy Observatory (TRAO) 14m single dish telescope equipped with a 16 multi-beam array covering $\sim$1.0 square degree area of this region using C$^{18}$O (1-0) mainly tracing low density cloud and about 460 square arcminute area using N$_{2}$H$^{+}$ (1-0) mainly tracing dense cores. CS (2-1) and SO $(3_{2}-2_{1})$ were also used simultaneously to map $\sim$440 square arcminute area of this region. We identified 10 filaments by applying the dendrogram technique to the C$^{18}$O data-cube and 8 dense N$_{2}$H$^{+}$ cores by using {\sc FellWalker}. Basic physical properties of filaments such as mass, length, width, velocity field, and velocity dispersion are derived. It is found that L1478 consists of several filaments with slightly different velocities. Especially the filaments which are supercritical are found to contain dense cores detected in N$_{2}$H$^{+}$. Comparison of non-thermal velocity dispersions derived from C$^{18}$O and N$_{2}$H$^{+}$ for the filaments and dense cores indicates that some of dense cores share similar kinematics with those of the surrounding filaments while several dense cores have different kinematics with those of their filaments. This suggests that the formation mechanism of dense cores and filaments can be different in individual filaments depending on their morphologies and environments.

Galaxy clusters are the largest and most massive bound objects resulting from cosmic hierarchical structure formation. Baryons account for somewhat more than 10% of that mass, with roughly 90% of the baryonic matter distributed throughout the clusters as hot ($T>1$ keV), high-$\beta$, very weakly collisional plasma; the so-called "intracluster medium" (ICM). Cluster mergers, close gravitational encounters and accretion, along with violent feedback from galaxies and relativistic jets from active galactic nuclei, drive winds, gravity waves, turbulence and shocks within the ICM. Those dynamics, in turn, generate cluster-scale magnetic fields and accelerate and mediate the transport of high-energy charged particles. Kinetic-scale, collective plasma processes define the basic character and fundamental signatures of these ICM phenomena, which are observed primarily by X-ray and radio astronomers.

Giant radio relics in the outskirts of galaxy clusters are known to be lit up by the relativistic electrons produced via diffusive shock acceleration (DSA) in shocks with low sonic Mach numbers, $M_{\rm s}\lesssim3$. The particle acceleration at these collisionless shocks critically depends on the kinetic plasma processes that govern the injection to DSA. Here, we study the preacceleration of suprathermal electrons in weak, quasi-perpendicular ($Q_\perp$) shocks in the hot, high-$\beta$ ($\beta = P_{\rm gas}/P_{\rm B}$) intracluster medium (ICM) through two-dimensional particle-in-cell simulations. \citet{guo2014a,guo2014b} showed that in high-$\beta$ $Q_\perp$-shocks, some of incoming electrons could be reflected upstream and gain energy via shock drift acceleration (SDA). The temperature anisotropy due to the SDA-energized electrons then induces the electron firehose instability (EFI), and oblique waves are generated, leading to a Fermi-like process and multiple cycles of SDA in the preshock region. We find that such electron preacceleration is effective only in shocks above a critical Mach number $M_{\rm ef}^*\approx2.3$. This means that in ICM plasmas, $Q_\perp$-shocks with $M_{\rm s}\lesssim2.3$ may not efficiently accelerate electrons. We also find that even in $Q_\perp$-shocks with $M_{\rm s}\gtrsim2.3$, electrons may not reach high enough energies to be injected to the full Fermi-I process of DSA, because long-wavelength waves are not developed via the EFI alone. Our results indicate that additional electron preaccelerations are required for DSA in ICM shocks, and the presence of fossil relativistic electrons in the shock upstream region may be necessary to explain observed radio relics.

In a growing number of galaxy clusters diffuse extended radio sources have been found. These sources are not directly associated with individual cluster galaxies. The radio emission reveal the presence of cosmic rays and magnetic fields in the intracluster medium (ICM). We classify diffuse cluster radio sources into radio halos, cluster radio shocks (relics), and revived AGN fossil plasma sources. Radio halo sources can be further divided into giant halos, mini-halos, and possible `intermediate' sources. Halos are generally positioned at cluster center and their brightness approximately follows the distribution of the thermal ICM. Cluster radio shocks (relics) are polarized sources mostly found in the cluster's periphery. They trace merger induced shock waves. Revived fossil plasma sources are characterized by their radio steep-spectra and often irregular morphologies. In this review we give an overview of the properties of diffuse cluster radio sources, with an emphasis on recent observational results. We discuss the resulting implications for the underlying physical acceleration processes that operate in the ICM, the role of relativistic fossil plasma, and the properties of ICM shocks and magnetic fields. We also compile an updated list of diffuse cluster radio sources which will be available on-line http://galaxyclusters.com. We end this review with a discussion on the detection of diffuse radio emission from the cosmic web.

Ultra-high-energy cosmic rays (UHECRs) are known to come from outside of our Galaxy, but their origin still remains unknown. The Telescope Array (TA) experiment recently identified a high concentration in the arrival directions of UHECRs with energies above $5.7 \times 10^{19} eV$, called hotspot. We here report the presence of filaments of galaxies, connected to the Virgo Cluster, on the sky around the hotspot, and a statistically significant correlation between hotspot events and the filaments. With 5-year TA data, the maximum significance of binomial statistics for the correlation is estimated to be 6.1 $\sigma$ at correlation angle 3.4 degree. The probability that the above significance appears by chance is $\sim 2.0 \times 10^{-8}$ (5.6 $\sigma$). Based on this finding, we suggest a model for the origin of TA hotspot UHECRs; they are produced at sources in the Virgo Cluster, and escape to and propagate along filaments, before they are scattered toward us. This picture requires the filament magnetic fields of strength $\gtrsim 20$ nG, which need to be confirmed in future observations.

Star formation takes place in the dense gas phase, and therefore a simple dense gas and star formation rate relation has been proposed. With the advent of multi-beam receivers, new observations show that the deviation from linear relations is possible. In addition, different dense gas tracers might also change significantly the measurement of dense gas mass and subsequently the relation between star formation rate and dense gas mass. We report the preliminary results the DEnse GAs in MAssive star-forming regions in the Milky Way (DEGAMA) survey that observed the dense gas toward a suit of well-characterized massive star forming regions in the Milky Way. Using the resulting maps of HCO$^{+}$ 1--0, HCN 1--0, CS 2--1, we discuss the current understanding of the dense gas phase where star formation takes place.

Galaxy clusters are known to host many active galaxies (AGNs) with radio jets, which could expand to form radio bubbles with relativistic electrons in the intracluster medium (ICM). It has been suggested that fossil relativistic electrons contained in remnant bubbles from extinct radio galaxies can be re-accelerated to radio-emitting energies by merger-driven shocks via diffusive shock acceleration (DSA), leading to the birth of radio relics detected in clusters. In this study we assume that such bubble consist primarily of thermal gas entrained from the surrounding medium and dynamically-insignificant amounts of relativistic electrons. We also consider several realistic models for magnetic fields in the cluster outskirts, including the ICM field that scales with the gas density as $B_{\rm ICM}\propto n_{\rm ICM}^{0.5}$. Then we perform time-dependent DSA simulations of a spherical shock that runs into a lower-density but higher-temperature bubble with the ratio $n_{\rm b}/n_{\rm ICM} \approx T_{\rm ICM}/T_{\rm b}\approx 0.5$. We find that inside the bubble the shock speed increases by about 20 \%, but the Mach number decreases by about 15\% in the case under consideration. In this re-acceleration model, the observed properties of a radio relic such as radio flux, spectral index, and integrated spectrum would be governed mainly by the presence of seed relativistic electrons and the magnetic field profile as well as shock dynamics. Thus it is crucial to understand how fossil electrons are deposited by AGNs in the ICM and how the downstream magnetic field evolves behind the shock in detailed modeling of radio relics.

Stellar halos observed around normal galaxies are extended and faint stellar structures formed by debris of tidally disrupted dwarf galaxies accreted over time by the host galaxy. Around dwarf galaxies, these stellar halos may not exist if all the accreted satellites are dark halos without stars. However, if a stellar halo is found in sufficiently small mass dwarfs, the whole stellar halo is composed of tidal debris of fossil galaxies, and we refer to it as ghostly halo. Fossil galaxies are called so because they formed most of their stars before the epoch of reionization, and have been identified as the ultra-faint dwarf galaxies found around the Milky Way and M31. In this paper we carry out semi-analytical simulations to characterize the sizes and stellar masses of ghostly stellar halos in dwarf galaxies as a function of their dark matter halo mass. By comparing the models to observations of six isolated dwarf galaxies in the Local Group showing evidence of extended stellar halos, we are able to constrain the star formation efficiency in fossil galaxies. We find that at redshift z~6, dark matter halos in the mass range 10^7-10^9 M_sun have a mean star formation efficiency f_* = M_*/M_{dm} ~ 0.05%-0.1% nearly constant as a function of the dark matter halo mass. There is a tentative indication of a sharp increase of f_* at M_{dm} ~ 10^6-10^7 M_sun, but it is very uncertain and based on only one dwarf galaxy (Leo~T).

Despite progress in understanding radio relics, there are still open questions regarding the underlying particle acceleration mechanisms. In this paper we present deep 1--4 GHz VLA observations of CIZA\,J2242.8+5301 ($z=0.1921$), a double radio relic cluster characterized by small projection on the plane of the sky. Our VLA observations reveal, for the first time, the complex morphology of the diffuse sources and the filamentary structure of the northern relic. We discover new faint diffuse radio emission extending north of the main northern relic. Our Mach number estimates for the northern and southern relics, based on the radio spectral index map obtained using the VLA observations and existing LOFAR and GMRT data, are consistent with previous radio and X-ray studies ($\mathcal{M}_{\rm RN}=2.58\pm0.17$ and $\mathcal{M}_{\rm RS}=2.10\pm0.08$). However, color-color diagrams and modelings suggest a flatter injection spectral index than the one obtained from the spectral index map, indicating that projection effects might be not entirely negligible. The southern relic consists of five "arms". Embedded in it, we find a tailed radio galaxy which seems to be connected to the relic. A spectral index flattening, where the radio tail connects to the relic, is also measured. We propose that the southern relic may trace AGN fossil electrons that are re-accelerated at a shock, with an estimated strength of $\mathcal{M}=2.4$. High-resolution mapping of other tailed radio galaxies also supports a scenario where AGN fossil electrons are revived by the merger event and could be related to the formation of some diffuse cluster radio emission.

Collisionless shocks with low sonic Mach numbers, $M_{\rm s} \lesssim 4$, are expected to accelerate cosmic ray (CR) protons via diffusive shock acceleration (DSA) in the intracluster medium (ICM). However, observational evidence for CR protons in the ICM has yet to be established. Performing particle-in-cell simulations, we study the injection of protons into DSA and the early development of a nonthermal particle population in weak shocks in high $\beta$ ($\approx 100$) plasmas. Reflection of incident protons, self-excitation of plasma waves via CR-driven instabilities, and multiple cycles of shock drift acceleration are essential to the early acceleration of CR protons in supercritical quasi-parallel shocks. We find that only in ICM shocks with $M_{\rm s} \gtrsim M_{\rm s}^*\approx 2.25$, a sufficient fraction of incoming protons are reflected by the overshoot in the shock electric potential and magnetic mirror at locally perpendicular magnetic fields, leading to efficient excitation of magnetic waves via CR streaming instabilities and the injection into the DSA process. Since a significant fraction of ICM shocks have $M_{\rm s} < M_{\rm s}^*$, CR proton acceleration in the ICM might be less efficient than previously expected. This may explain why the diffuse gamma-ray emission from galaxy clusters due to proton-proton collisions has not been detected so far.

Massive stars blow powerful stellar winds throughout their evolutionary stages from the main sequence to Wolf-Rayet phases. The amount of mechanical energy deposited in the interstellar medium by the wind from a massive star can be comparable to the explosion energy of a core-collapse supernova that detonates at the end of its life. In this study, we estimate the kinetic energy deposition by massive stars in our Galaxy by considering the integrated Galactic initial mass function and modeling the stellar wind luminosity. The mass loss rate and terminal velocity of stellar winds during the main sequence, red supergiant, and Wolf-Rayet stages are estimated by adopting theoretical calculations and observational data published in the literature. We find that the total stellar wind luminosity due to all massive stars in the Galaxy is about $L_w\approx 1.1\times 10^{41}$ erg/s, which is about 1/4 of the power of supernova explosions, $L_{SN} \approx 4.8\times 10^{41}$ erg/s. If we assume that $\sim 1-10$ % of the wind luminosity could be converted to Galactic cosmic rays (GCRs) through collisonless shocks such as termination shocks in stellar bubbles and superbubbles, colliding-wind shocks in binaries, and bow-shocks of massive runaway stars, stellar winds might be expected to make a significant contribution to GCR production, though lower than that of supernova remnants.

Non-detection of $\gamma$-ray emission from galaxy clusters has challenged diffusive shock acceleration (DSA) of cosmic-ray (CR) protons at weak collisionless shocks that are expected to form in the intracluster medium. As an effort to address this problem, we here explore possible roles of Alfv\'en waves self-excited via resonant streaming instability during the CR acceleration at parallel shocks. The mean drift of Alfv\'en waves may either increase or decrease the scattering center compression ratio, depending on the postshock cross-helicity, leading to either flatter or steeper CR spectra. We first examine such effects at planar shocks, based on the transport of Alfv\'en waves in the small amplitude limit. For the shock parameters relevant to cluster shocks, Alfv\'enic drift flattens the CR spectrum slightly, resulting in a small increase of the CR acceleration efficiency, $\eta$. We then consider two additional, physically motivated cases: (1) postshock waves are isotropized via MHD and plasma processes across the shock transition and (2) postshock waves contain only forward waves propagating along with the flow due to a possible gradient of CR pressure behind the shock. In these cases, Alfv\'enic drift could reduce $\eta$ by as much as a factor of 5 for weak cluster shocks. For the canonical parameters adopted here, we suggest $\eta\sim10^{-4}-10^{-2}$ for shocks with sonic Mach number $M_{\rm s}\approx2-3$. The possible reduction of $\eta$ may help ease the tension between non-detection of $\gamma$-rays from galaxy clusters and DSA predictions.

X-ray shocks and radio relics detected in the cluster outskirts are commonly interpreted as shocks induced by mergers of sub-clumps. We study the properties of merger shocks in merging galaxy clusters, using a set of cosmological simulations for the large-scale structure formation of the universe. As a representative case, we here focus on the simulated clusters that undergo almost head-on collisions with mass ratio $\sim2$. Due to the turbulent nature of the intracluster medium, shock surfaces are not smooth, but composed of shocks with different Mach numbers. As the merger shocks expand outward from the core to the outskirts, the average Mach number, $\left<M_s\right>$, increases in time. We suggest that the shocks propagating along the merger axis could be manifested as X-ray shocks and/or radio relics. The kinetic energy through the shocks, $F_\phi$, peaks at $\sim1$ Gyr after their initial launching, or at $\sim1-2$ Mpc from the core. Because of the Mach number dependent model adopted here for the cosmic ray (CR) acceleration efficiency, their CR-energy-weighted Mach number is higher with $\left< M_s \right>_{\rm CR}\sim3-4$, compared to the kinetic-energy-weighted Mach number, $\left<M_s\right>_{\phi}\sim2-3$. Most energetic shocks are to be found ahead of the lighter dark matter (DM) clump, while the heavier DM clump is located in the opposite side of clusters. Although our study is limited to the merger case considered, the results such as the means and variations of shock properties and their time evolution could be compared with the observed characteristics of merger shocks, constraining interpretations of relevant observations.

According to structure formation simulations, weak shocks with typical Mach number, $M_{\rm s}\lesssim 3$, are expected to form in merging galaxy clusters. The presence of such shocks has been indicated by X-ray and radio observations of many merging clusters. In particular, diffuse radio sources known as radio relics could be explained by synchrotron-emitting electrons accelerated via diffusive shock acceleration (Fermi I) at quasi-perpendicular shocks. Here we also consider possible roles of stochastic acceleration (Fermi II) by compressive MHD turbulence downstream of the shock. Then we explore a puzzling discrepancy that for some radio relics, the shock Mach number inferred from the radio spectral index is substantially larger than that estimated from X-ray observations. This problem could be understood, if shock surfaces associated with radio relics consist of multiple shocks with different strengths.In that case, X-ray observations tend to pick up the part of shocks with lower Mach numbers and higher kinetic energy flux, while radio emissions come preferentially from the part of shocks with higher Mach numbers and higher cosmic ray (CR) production. We also show that the Fermi I reacceleration model with preexisting fossil electrons supplemented by Fermi II acceleration due to postshock turbulence could reproduce observed profiles of radio flux densities and integrated radio spectra of two giant radio relics. This study demonstrates the CR electrons can be accelerated at collisionless shocks in galaxy clusters just like supernova remnant shock in the interstellar medium and interplanetary shocks in the solar wind.

Cosmological hydrodynamic simulations have demonstrated that shock waves could be produced in the intergalactic medium by supersonic flow motions during the course of hierarchical clustering of the large-scale-structure in the Universe. Similar to interplanetary shocks and supernova remnants (SNRs), these structure formation shocks can accelerate cosmic ray (CR) protons and electrons via diffusive shock acceleration. External accretion shocks, which form in the outermost surfaces of nonlinear structures, are as strong as SNR shocks and could be potential accelerations sites for high energy CR protons up to $10^{18}$ eV. But it could be difficult to detect their signatures due to extremely low kinetic energy flux associated with those accretion shocks. On the other hand, radiative features of internal shocks in the hot intracluster medium have been identified as temperature and density discontinuities in X-ray observations and diffuse radio emission from accelerated CR electrons. However, the non-detection of gamma-ray emission from galaxy clusters due to $\pi^0$ decay still remains to be an outstanding problem.

We explore the shock acceleration model for giant radio relics, in which relativistic electrons are accelerated via diffusive shock acceleration (DSA) by merger-driven shocks in the outskirts of galaxy clusters. In addition to DSA, turbulent acceleration by compressive MHD mode downstream of the shock is included as well as energy losses of postshock electrons by Coulomb scattering, synchrotron emission, and inverse Compton scattering off the cosmic background radiation. Considering that only a small fraction of merging clusters host radio relics, we favor the reacceleration scenario in which radio relics are generated preferentially when shocks encounter the regions containing low-energy ($\gamma_{\rm e} \lesssim 300$) cosmic ray electrons (CRe). We perform time-dependent DSA simulations of spherically expanding shocks with physical parameters relevant for the Sausage radio relic, and calculate the radio synchrotron emission from the accelerated CRe. We find that significant level of postshock turbulent acceleration is required in order to reproduce broad profiles of the observed radio flux densities of the Sausage relic. Moreover, the spectral curvature in the observed integrated radio spectrum can be explained, if the putative shock should have swept up and exited out of the preshock region of fossil CRe about 10~Myr ago.

Although many of the observed properties of giant radio relics detected in the outskirts of galaxy clusters can be explained by relativistic electrons accelerated at merger-driven shocks, significant puzzles remain. In the case of the so-called Toothbrush relic, the shock Mach number estimated from X-ray observations ($M_{\rm X}\approx1.2-1.5$) is substantially weaker than that inferred from the radio spectral index ($M_{\rm rad}\approx2.8$).Toward understanding such a discrepancy, we here consider the following diffusive shock acceleration (DSA) models:(1) weak-shock models with $M_{\rm s}\lesssim 2$ and a preexisting population of cosmic-ray electrons (CRe) with a flat energy spectrum,and (2) strong-shock models with $M_{\rm s}\approx3$ and either shock-generated suprathermal electrons or preexisting fossil CRe. We calculate the synchrotron emission from the accelerated CRe, following the time evolution of the electron DSA, and subsequent radiative cooling and postshock turbulent acceleration (TA). We find that both models could reproduce reasonably well the observed integrated radio spectrum of the Toothbrush relic, but the observed broad transverse profile requires the stochastic acceleration by downstream turbulence, which we label "turbulent acceleration" or TA to distinguish it from DSA. Moreover, to account for the almost uniform radio spectral index profile along the length of the relic, the weak-shock models require a preshock region over 400~kpc with a uniform population of preexisting CRe with a high cutoff energy ($\gtrsim 40$ GeV). Due to the short cooling time, it is challenging to explain the origin of such energetic electrons. Therefore, we suggest the strong-shock models with low-energy seed CRe ($\lesssim 150$~MeV) are preferred for the radio observations of this relic.

On the largest scales, the Universe consists of voids and filaments making up the cosmic web. Galaxy clusters are located at the knots in this web, at the intersection of filaments. Clusters grow through accretion from these large-scale filaments and by mergers with other clusters and groups. In a growing number of galaxy clusters, elongated Mpc-size radio sources have been found, so-called radio relics. These relics are thought to trace relativistic electrons in the intracluster plasma accelerated by low-Mach number collisionless shocks generated by cluster-cluster merger events. A long-standing problem is how low-Mach number shocks can accelerate electrons so efficiently to explain the observed radio relics. Here we report on the discovery of a direct connection between a radio relic and a radio galaxy in the merging galaxy cluster Abell 3411-3412. This discovery indicates that fossil relativistic electrons from active galactic nuclei are re-accelerated at cluster shocks. It also implies that radio galaxies play an important role in governing the non-thermal component of the intracluster medium in merging clusters.

The Sausage radio relic is the arc-like radio structure in the cluster CIZA J2242.8+5301, whose observed properties can be best understood by synchrotron emission from relativistic electrons accelerated at a merger-driven shock. However, there remain a few puzzles that cannot be explained by the shock acceleration model with only in-situ injection. In particular, the Mach number inferred from the observed radio spectral index, $M_{\rm radio}\approx 4.6$, while the Mach number estimated from X-ray observations, $M_{\rm X-ray}\approx 2.7$. In an attempt to resolve such a discrepancy, here we consider the re-acceleration model in which a shock of $M_s\approx 3$ sweeps through the intracluster gas with a pre-existing population of relativistic electrons. We find that observed brightness profiles at multi frequencies provide strong constraints on the spectral shape of pre-existing electrons. The models with a power-law momentum spectrum with the slope, $s\approx 4.1$, and the cutoff Lorentz factor, $\gamma_{e,c}\approx 3-5\times 10^4$ can reproduce reasonably well the observed spatial profiles of radio fluxes and integrated radio spectrum of the Sausage relic. The possible origins of such relativistic electrons in the intracluster medium remain to be investigated further.

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.

The Toothbrush radio relic associated the merging cluster 1RXS J060303.3 is presumed to be produced by relativistic electrons accelerated at merger-driven shocks. Since the shock Mach number inferred from the observed radio spectral index, $M_{radio}\approx 2.8$, is larger than that estimated from X-ray observations, $M_{X-ray}\lesssim 1.5$, we consider the re-acceleration model in which a weak shock of $M_s\approx 1.2-1.5$ sweeps through the intracluster plasma with a preshock population of relativistic electrons. We find the models with a power-law momentum spectrum with the slope, $s\approx 4.6$, and the cutoff Lorentz factor, $\gamma_{e,c}\approx 7-8\times 10^4$ can reproduce reasonably well the observed profiles of radio fluxes and integrated radio spectrum of the head portion of the Toothbrush relic. This study confirms the strong connection between the ubiquitous presence of fossil relativistic plasma originated from AGNs and the shock-acceleration model of radio relics in the intracluster medium.

The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected.

We present a simultaneous single-dish survey of 22 GHz water maser and 44 GHz and 95 GHz class I methanol masers toward 77 6.7 GHz class II methanol maser sources, which were selected from the Arecibo methanol maser Galactic plane survey (AMGPS) catalog.Water maser emission is detected in 39 (51%) sources, of which 15 are new detections. Methanol maser emission at 44 GHz and 95 GHz is found in 25 (32%) and 19 (25%) sources, of which 21 and 13 sources are newly detected, respectively. We find 4 high-velocity (> 30 km/s) water maser sources, including 3 dominant blue- or redshifted outflows.The 95 GHz masers always appear with the 44 GHz maser emission. They are strongly correlated with 44 GHz masers in velocity, flux density, and luminosity, while they are not correlated with either water or 6.7 GHz class II methanol masers. The average peak flux density ratio of 95 GHz to 44 GHz masers is close to unity, which is two times higher than previous estimates. The flux densities of class I methanol masers are more closely correlated with the associated BGPS core mass than those of water or class II methanol masers. Using the large velocity gradient (LVG) model and assuming unsaturated class I methanol maser emission, we derive the fractional abundance of methanol to be in a range of 4.2*10^-8 to 2.3*10^-6, with a median value of 3.3\pm2.7*10^-7.

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.

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) \nu_{\rm br}$ with break frequency $ \nu_{\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 $\gamma_e \gtrsim 10^4$.

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

We study the evolution of the energy spectrum of cosmic-ray electrons accelerated at spherically expanding shocks with low Mach numbers and the ensuing spectral signatures imprinted in radio synchrotron emission. Time-dependent simulations of diffusive shock acceleration (DSA) of electrons in the test-particle limit have been performed for spherical shocks with parameters relevant for typical shocks in the intracluster medium. The electron and radiation spectra at the shock location can be described properly by the test-particle DSA predictions with instantaneous shock parameters. However, the volume integrated spectra of both electrons and radiation deviate significantly from the test-particle power-laws, because the shock compression ratio and the flux of injected electrons at the shock gradually decrease as the shock slows down in time.So one needs to be cautious about interpreting observed radio spectra of evolving shocks based on simple DSA models in the test-particle regime.

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 $\kappa$-distributions produced by interactions of background thermal particles with pre-existing and/or self-excited plasma/MHD waves or turbulence. The $\kappa$-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 $\kappa$-distribution depends on the shock Mach number much less severely compared to that for the Maxwellian distribution. Thus, the existence of $\kappa$-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 $\kappa$-distribution with $10 < \kappa_p < 30$ at quasi-parallel shocks, while the injection fraction for electrons becomes $10^{-6}-10^{-5}$ for a $\kappa$-distribution with $\kappa_e < 2$ at quasi-perpendicular shocks. For such $\kappa$ 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.

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.

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

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.

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.

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

We present the results of 12CO J = 1-0 line observations of eleven Galactic supernova remnants (SNRs) obtained using the Seoul Radio Astronomy Observatory (SRAO) 6-m radio telescope. The observation was made as a part of the SRAO CO survey of SNRs between l = 70 and 190 deg, which is intended to identify SNRs interacting with molecular clouds. The mapping areas for the individual SNRs are determined to cover their full extent in the radio continuum. We used halfbeam grid spacing (60") for 9 SNRs and full-beam grid spacing (120") for the rest. We detected CO emission towards most of the remnants. In six SNRs, molecular clouds showed a good spatial relation with their radio morphology, although no direct evidence for the interaction was detected. Two SNRs are particularly interesting: G85.4+0.7, where there is a filamentary molecular cloud along the radio shell, and 3C434.1, where a large molecular cloud appears to block the western half of the remnant. We briefly summarize the results obtained for individual SNRs.

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.

We report a multi-epoch, simultaneous 22 GHz H2O and 44 GHz class I CH3OH maser line survey towards 180 intermediate-mass young stellar objects, including 14 Class 0, 19 Class I objects, and 147 Herbig Ae/Be stars. We detected H2O and CH3OH maser emission towards 16 (9%) and 10 (6%) sources with one new H2O and six new CH3OH maser sources. The detection rates of both masers rapidly decrease as the central (proto)stars evolve, which is contrary to the trends in high-mass star-forming regions. This suggests that the excitations of the two masers are closely related to the evolutionary stage of the central (proto)stars and the circumstellar environments. H2O maser velocities deviate on average 9 km s^-1 from the ambient gas velocities whereas CH3OH maser velocities match quite well with the ambient gas velocities. For both maser emissions, large velocity differences (|v_{H2O} - v_{sys} | > 10 km s^-1 and |v_{CH3OH} - v_{sys}| > 1 km s^-1) are mostly confined to Class 0 objects. The formation and disappearance of H2O masers is frequent and their integrated intensities change by up to two orders of magnitude. In contrast, CH3OH maser lines usually show no significant change in intensity, shape, or velocity. This is consistent with the previous suggestion that H2O maser emission originates from the base of an outflow while 44 GHz class I CH3OH maser emission arises from the interaction region of the outflow with the ambient gas. The isotropic maser luminosities are well correlated with the bolometric luminosities of the central objects. The fitted relations are L_{H2O} = 1.71 * 10^{-9} (L_{bol})^{0.97} and L_{CH3OH} = 1.71 * 10^{-10} (L_{bol})^{1.22}.

[...] We present results for the statistics of thermal gas and the shock wave properties for a large volume simulated with three different cosmological numerical codes: the Eulerian total variations diminishing code TVD, the Eulerian piecewise parabolic method-based code ENZO, and the Lagrangian smoothed-particle hydrodynamics code GADGET. Starting from a shared set of initial conditions, we present convergence tests for a cosmological volume of side-length 100 Mpc/h, studying in detail the morphological and statistical properties of the thermal gas as a function of mass and spatial resolution in all codes. By applying shock finding methods to each code, we measure the statistics of shock waves and the related cosmic ray acceleration efficiencies, within the sample of simulations and for the results of the different approaches. We discuss the regimes of uncertainties and disagreement among codes, with a particular focus on the results at the scale of galaxy clusters. Even if the bulk of thermal and shock properties are reasonably in agreement among the three codes, yet some significant differences exist (especially between Eulerian methods and smoothed particle hydrodynamics). In particular, we report: a) differences of huge factors (10-100) in the values of average gas density, temperature, entropy, Mach number and shock thermal energy flux in the most rarefied regions of the simulations between grid and SPH methods; b) the hint of an entropy core inside clusters simulated in grid codes; c) significantly different phase diagrams of shocked cells in grid codes compared to SPH; d) sizable differences in the morphologies of accretion shocks between grid and SPH methods.

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

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.

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.

We calculate the energy spectra of cosmic ray (CR) protons and electrons at a plane shock with quasi-parallel magnetic fields, using time-dependent, diffusive shock acceleration (DSA) simulations, including energy losses via synchrotron emission and Inverse Compton (IC) scattering. A thermal leakage injection model and a Bohm type diffusion coefficient are adopted. The electron spectrum at the shock becomes steady after the DSA energy gains balance the synchrotron/IC losses, and it cuts off at the equilibrium momentum p_{eq}. In the postshock region the cutoff momentum of the electron spectrum decreases with the distance from the shock due to the energy losses and the thickness of the spatial distribution of electrons scales as p^{-1}. Thus the slope of the downstream integrated spectrum steepens by one power of p for p_{br}<p<p_{eq}, where the break momentum decrease with the shock age as p_{br}\propto t^{-1}. In a CR modified shock, both the proton and electron spectrum exhibit a concave curvature and deviate from the canonical test-particle power-law, and the upstream integrated electron spectrum could dominate over the downstream integrated spectrum near the cutoff momentum. Thus the spectral shape near the cutoff of X-ray synchrotron emission could reveal a signature of nonlinear DSA.

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\'enic 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.

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, \xi_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, \xi_c < 10^{-3} for shocks with M < 5, while \xi_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.

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

We perform kinetic simulations of diffusive shock acceleration (DSA) in Type Ia supernova remnants (SNRs) expanding into a uniform interstellar medium (ISM). Bohm-like diffusion assumed, and simple models for Alfvenic drift and dissipation are adopted. Phenomenological models for thermal leakage injection are considered as well. We find that the preshock gas temperature is the primary parameter that governs the cosmic ray (CR) acceleration efficiency and energy spectrum, while the CR injection rate is a secondary parameter. For SNRs in the warm ISM, if the injection fraction is larger than 10^{-4}, the DSA is efficient enough to convert more than 20 % of the SN explosion energy into CRs and the accelerated CR spectrum exhibits a concave curvature flattening to E^{-1.6}. Such a flat source spectrum near the knee energy, however, may not be reconciled with the CR spectrum observed at Earth. On the other hand, SNRs in the hot ISM, with an injection fraction smaller than 10^{-4}, are inefficient accelerators with less than 10 % of the explosion energy getting converted to CRs. Also the shock structure is almost test-particle like and the ensuing CR spectrum can be steeper than E^{-2}. With amplified magnetic field strength of order of 30 microG, Alfven waves generated by the streaming instability may drift upstream fast enough to make the modified test-particle power-law as steep as E^{-2.3}, which is more consistent with the observed CR spectrum.

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 $<\theta > \sim 15^{\circ}$ and the median value of $\tilde \theta \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.

Ultrahigh energy cosmic rays (UHECRs) are believed to originate from astrophysical sources, which should trace the large scale structure (LSS) of the universe. On the other hand, the magnetic field in the intergalactic space (IGMF), which also traces the LSS of the universe, deflects the trajectories of the charged UHECRs and spoils the positional correlation of the observed UHECR events with their true sources. To explore this problem, we studied a simulation of the propagation of UHE protons through the magnetized LSS of the universe, reported earlier in Das et al. (2008), in which the IGMF was estimated based on a turbulence dynamo model (Ryu et al. 2008). Hypothetical sources were placed inside clusters and groups of galaxies in the simulated universe, while observers were located inside groups of galaxies that have similar properties as the Local Group. We calculated the statistics of the angular distance between the arrival directions of simulated UHE proton events and the positions of candidate sources in our simulation. We compared the statistics from our simulation with those calculated with the Auger data. We discussed the implication of our works on the nature of the sources of UHECRs.

For over four decades, synchrotron-radiating sources have played a series of pathfinding roles in the study of galaxy clusters and large scale structure. Such sources are uniquely sensitive to the turbulence and shock structures of large-scale environments, and their cosmic rays and magnetic fields often play important dynamic and thermodynamic roles. They provide essential complements to studies at other wavebands. Over the next decade, they will fill essential gaps in both cluster astrophysics and the cosmological growth of structure in the universe, especially where the signatures of shocks and turbulence, or even the underlying thermal plasma itself, are otherwise undetectable. Simultaneously, synchrotron studies offer a unique tool for exploring the fundamental question of the origins of cosmic magnetic fields. This work will be based on the new generation of m/cm-wave radio telescopes now in construction, as well as major advances in the sophistication of 3-D MHD simulations.

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\'en 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\'enic 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.

This paper reports the discovery of a Very Low Luminosity Object (VeLLO) in the "starless" dense core L328, using the Spitzer Space Telescope and ground based observations from near-infrared to millimeter wavelengths. The Spitzer 8 micron image indicates that L328 consists of three subcores of which the smallest one may harbor a source, L328-IRS while two other subcores remain starless. L328-IRS is a Class 0 protostar according to its bolometric temperature (44 K) and the high fraction ~72 % of its luminosity emitted at sub-millimeter wavelengths. Its inferred "internal luminosity" (0.04 - 0.06 Lsun) using a radiative transfer model under the most plausible assumption of its distance as 200 pc is much fainter than for a typical protostar, and even fainter than other VeLLOs studied previously. Note, however, that its inferred luminosity may be uncertain by a factor of 2-3 if we consider two extreme values of the distance of L328-IRS (125 or 310 pc). Low angular resolution observations of CO do not show any clear evidence of a molecular outflow activity. But broad line widths toward L328, and Spitzer and near-infrared images showing nebulosity possibly tracing an outflow cavity, strongly suggest the existence of outflow activity. Provided that an envelope of at most ~0.1 Msunis the only mass accretion reservoir for L328-IRS, and the star formation efficiency is close to the canonical value ~30%, L328-IRS has not yet accreted more than 0.05 Msun. At the assumed distance of 200 pc, L328-IRS is destined to be a brown dwarf.

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.

The nature and origin of turbulence and magnetic fields in the intergalactic space are important problems that are yet to be understood. We propose a scenario in which turbulent flow motions are induced via the cascade of the vorticity generated at cosmological shocks during the formation of the large scale structure. The turbulence in turn amplifies weak seed magnetic fields of any origin. Supercomputer simulations show that the turbulence is subsonic inside clusters/groups of galaxies, whereas it is transonic or mildly supersonic in filaments. Based on a turbulence dynamo model, we then estimate that the average magnetic field strength would be a few microgauss inside clusters/groups, approximately 0.1 microgauss around clusters/groups, and approximately 10 nanogauss in filaments. Our model presents a physical mechanism that transfers the gravitation energy to the turbulence and magnetic field energies in the large scale structure of the universe.

If ultra-high-energy cosmic rays (UHECRs) originate from extragalactic sources, understanding the propagation of charged particles through the magnetized large scale structure (LSS) of the universe is crucial in the search for the astrophysical accelerators. Based on a novel model of the turbulence dynamo, we estimate the intergalactic magnetic fields (IGMFs) in cosmological simulations of the formation of the LSS. Under the premise that the sources of UHECRs are strongly associated with the LSS, we consider a model in which protons with E >10^{19} eV are injected by sources that represent active galactic nuclei located inside clusters of galaxies. With the model IGMFs, we then follow the trajectories of the protons, while taking into account the energy losses due to interactions with the cosmic background radiation. For observers located inside groups of galaxies like ours, about 70% and 35% of UHECR events above 60 EeV arrive within ~15 degree and ~5 degree, respectively, of the source position with time delays of less than ~10^7 yr. This implies that the arrival direction of super-GZK protons might exhibit a correlation with the distribution of cosmological sources on the sky. In this model, nearby sources (within 10 - 20 Mpc) should contribute significantly to the particle flux above ~10^{20} eV.

We study the evolution of multiple supernova (SN) explosions inside a pre-exiting cavity blown by winds from massive progenitor stars. Hydrodynamic simulations in one-dimensional spherical geometry, including radiative cooling and thermal conduction, are carried out to follow first the development of the wind-blown bubble during the main sequence and then the evolution of the SN-driven bubble. We find the size and mass of the SN-driven bubble shell depend on the structure of the pre-existing wind bubble as well as the SN explosion energy E_{SN} (= N_{SN} 10^{51} ergs). The hot cavity inside the bubble is 2-3 times bigger in volume and hotter than that of a bubble created by SNe exploded in a uniform interstellar medium (ISM). For an association with 10 massive stars in the average ISM, the SN-driven shell has an outer radius of R_{ss} ~ (85 pc) N_{SN}^{0.1} and a mass of M_{ss} ~ (10^{4.8} Msun) N_{SN}^{0.3}at 10^6 years after the explosion. By that time most of the explosion energy is lost via radiative cooling, while ~10% remains as kinetic energy and ~10% as thermal energy. We also calculate the total integrated spectrum of diffuse radiation emitted by the shock-heated gas of the SN bubble. For the models with 0.1 solar metalicity, the radiative energy loss is smaller and the fraction of non-ionizing photons is larger, compared to those with solar metalicity. We conclude the photoionization/heating by diffuse radiation is the most dominant form of feedback from SN explosions into the surrounding medium.

The propagation of UHECRs is affected by the intergalactic magnetic fields that were produced during the course of the large scale structure formation of the universe. We adopt a novel model where the large scale extragalactic magnetic fields (EGMF) are estimated from local dynamic properties of the gas flows in hydrodynamic simulations of a concordance LambdaCDM universe. With the model magnetic fields, we calculate the deflection angle, time delay and energy spectrum of protons with E > 10^{19} eV that are injected at cosmological sources and then travel through the large scale structure of the universe, losing the energy due to interactions with the cosmic background radiation. Implications of this study on the origin of UHECRs are discussed.

We present results of time evolution of oblique MHD plane shocks including diffusive cosmic ray acceleration with backreaction on the plasma flows. The simulations include self-consistent effects of finite Alfven wave propagation and dissipation. From the computed cosmic ray particle phase space distributions we calculate expected leptonic and hadronic emissions resulting from interactions between the cosmic rays, magnetic fields, the thermal particle population and relevant astrophysical photon fields.

Using an improved version of the previously introduced CRASH (Cosmic Ray Acceleration SHock) code, we have calculated the time evolution of cosmic-ray (CR) modified quasi-parallel plane shocks for Bohm-like diffusion, including self-consistent models of Alfven wave drift and dissipation, along with thermal leakage injection of CRs. The new simulations follow evolution of the CR distribution to much higher energies than our previous study, providing a better examination of evolutionary and asymptotic behaviors. The postshock CR pressure becomes constant after quick initial adjustment, since the evolution of the CR partial pressure expressed in terms of a momentum similarity variable is self-similar. The shock precursor, which scales as the diffusion length of the highest energy CRs, subsequently broadens approximately linearly with time, independent of diffusion model, so long as CRs continue to be accelerated to ever-higher energies. This means the nonlinear shock structure can be described approximately in terms of the similarity variable, x/(u_s t), where u_s is the shock speed once the postshock pressure reaches an approximate time asymptotic state. As before, the shock Mach number is the key parameter determining the evolution and the CR acceleration efficiency, although finite Alfven wave drift and wave energy dissipation in the shock precursor reduce the effective velocity change experienced by CRs, so reduce acceleration efficiency noticeably, thus, providing a second important parameter at low and moderate Mach numbers.

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.

We have calculated the cosmic ray (CR) acceleration at young remnants from Type Ia supernovae expanding into a uniform interstellar medium (ISM). Adopting quasi-parallel magnetic fields, gasdynamic equations and the diffusion convection equation for the particle distribution function are solved in a comoving spherical grid which expands with the shock. Bohm-type diffusion due to self-excited Alfven waves, drift and dissipation of these waves in the precursor and thermal leakage injection were included. With magnetic fields amplified by the CR streaming instability, the particle energy can reach up to 10^{16}Z eV at young supernova remnants (SNRs) of several thousand years old. The fraction of the explosion energy transferred to the CR component asymptotes to 40-50 % by that time. For a typical SNR in a warm ISM, the accelerated CR energy spectrum should exhibit a concave curvature with the power-law slope flattening from 2 to 1.6 at E>0.1 TeV.

Thermal-gravitational instability (TGI) is studied in the protogalactic environment. We extend our previous work, where we found that dense clumps first form out of hot background gas by thermal instability and later a small fraction of them grow to virialized clouds of mass M_c >~ 6X10^6 M_sun by gravitational infall and merging. But these clouds have large angular momentum, so they would be difficult, if not impossible, to further evolve into globular clusters. In this paper, through three-dimensional hydrodynamic simulations in a uniformly rotating frame, we explore if the Coriolis force due to rotation in protogalactic disk regions can hinder binary merging and reduce angular momentum of the clouds formed. With rotation comparable to the Galactic rotation at the Solar circle, the Coriolis force is smaller than the pressure force during the early thermal instability stage. So the properties of clumps formed by thermal instability are not affected noticeably by rotation, except increased angular momentum. However, during later stage the Coriolis force becomes dominant over the gravity, and hence the further growth to gravitationally bound clouds by gravitational infall and merging is prohibited. Our results show that the Coriolis force effectively destroys the picture of cloud formation via TGI, rather than alleviate the problem of large angular momentum.

We have developed a cosmic ray (CR) shock code in one dimensional spherical geometry with which the particle distribution, the gas flow and their nonlinear interaction can be followed numerically in a frame comoving with an expanding shock. In order to accommodate a very wide dynamic range of diffusion length scales in the CR shock problem, we have incorporated subzone shock tracking and adaptive mesh refinement techniques. We find the spatial grid resolution required for numerical convergence is less stringent in this code compared to typical, fixed-grid Eulerian codes. The improved convergence behavior derives from maintaining the shock discontinuity inside the same grid zone in the comoving code. That feature improves numerical estimates of the compression rate experienced by CRs crossing the subshock compared to codes that allow the subshock to drift on the grid. Using this code with a Bohm-like diffusion model we have calculated the CR acceleration and the nonlinear feedback at supernova remnant shocks during the Sedov-Taylor stage. Similarly to plane-parallel shocks, with an adopted thermal leakage injection model, about 10^{-3} of the particles that pass through the shock and up to 60 % of the explosion energy are transferred to the CR component. These results are in good agreement with previous nonlinear spherical CR shock calculations of Berezhko and collaborators.

We have developed a new, very efficient numerical scheme to solve the CR diffusion convection equation that can be applied to the study of the nonlinear time evolution of CR modified shocks for arbitrary spatial diffusion properties. The efficiency of the scheme derives from its use of coarse-grained finite momentum volumes. This approach has enabled us, using $\sim 10 - 20$ momentum bins spanning nine orders of magnitude in momentum, to carry out simulations that agree well with results from simulations of modified shocks carried out with our conventional finite difference scheme requiring more than an order of magnitude more momentum points. The coarse-grained, CGMV scheme reduces execution times by a factor approximately half the ratio of momentum bins used in the two methods. Depending on the momentum dependence of the diffusion, additional economies in required spatial and time resolution can be utilized in the CGMV scheme, as well. These allow a computational speed-up of at least an order of magnitude in some cases.

We study thermal-gravitational instability in simplified models for protogalactic halos using three-dimensional hydrodynamic simulations. The simulations followed the evolution of gas with radiative cooling down to T = 10^4 K, background heating, and self-gravity. Then cooled and condensed clouds were identified and their physical properties were examined in detail. During early stage clouds start to form around initial density peaks by thermal instability. Small clouds appear first and they are pressure-bound. Subsequently, the clouds grow through compression by the background pressure as well as gravitational infall. During late stage cloud-cloud collisions become important, and clouds grow mostly through gravitational merging. Gravitationally bound clouds with mass M_c > ~6 X 10^6 Msun are found in the late stage. They are approximately in virial equilibrium and have radius R_c = \~150 - 200 pc. Those clouds have gained angular momentum through tidal torque as well as merging, so they have large angular momentum with the spin parameter <lambda_s> ~ 0.3. The clouds formed in a denser background tend to have smaller spin parameters. We discuss briefly the implications of our results on the formation of protoglobular cluster clouds in protogalactic halos. (abridged)

We have calculated the evolution of cosmic ray (CR) modified astrophysical shocks for a wide range of shock Mach numbers and shock speeds through numerical simulations of diffusive shock acceleration (DSA) in 1D quasi- parallel plane shocks. The simulations include thermal leakage injection of seed CRs, as well as pre-existing, upstream CR populations. Bohm-like diffusion is assumed. We model shocks similar to those expected around cosmic structure pancakes as well as other accretion shocks driven by flows with upstream gas temperatures in the range $T_0=10^4-10^{7.6}$K and shock Mach numbers spanning $M_s=2.4-133$. We show that CR modified shocks evolve to time-asymptotic states by the time injected particles are accelerated to moderately relativistic energies ($p/mc \gsim 1$), and that two shocks with the same Mach number, but with different shock speeds, evolve qualitatively similarly when the results are presented in terms of a characteristic diffusion length and diffusion time. For these models the time asymptotic value for the CR acceleration efficiency is controlled mainly by shock Mach number. The modeled high Mach number shocks all evolve towards efficiencies $\sim 50$%, regardless of the upstream CR pressure. On the other hand, the upstream CR pressure increases the overall CR energy in moderate strength shocks ($M_s \sim {\rm a few}$). (abridged)

Cosmological shock waves play a major role in heating baryonic gas in the large scale structure of the universe. In this contribution we study the shock-heated gas in connection with shocks themselves. The distributions of shock speed and temperature of shock-hated gas both should reflect the depth of gravitational potential of associated nonlinear structures, so their morphology. From their distributions we find, in addition to hot gas in and around clusters/groups and warm-hot intergalactic medium (WHIM) with T=10^5-10^7 K mostly in filaments, there is a significant amount of low temperature WHIM with T < 10^5 K distributed mostly as sheet-like structures. Not only the WHIM with T=10^5-10^7 K but also the WHIM with T < 10^5 K make up a significant fraction of gas mass, implying the low temperature WHIM could be important in resolving the missing baryon problem. The shock-heated gas in filaments and sheets are manifested best through emissions and absorptions in soft X-ray and far UV. We confirm that the WHIM with T=10^5-10^7 K makes significant contributions to the soft X-ray background, absorptions of highly ionized species such as OVII and OVIII in AGN spectra, and line emissions from OVII and OVIII ions, as pointed by previous studies. However, the WHIM with T < 10^5 K is the major contributor to absorptions of lower ionized species such as OV and OVI, because these photo-ionized ions are most abundant in sheets of low density and temperature. (abridged)

We study the properties of cosmological shock waves identified in high-resolution, N-body/hydrodynamic simulations of a $\Lambda$CDM universe and their role on thermalization of gas and acceleration of nonthermal, cosmic ray (CR) particles. External shocks form around sheets, filaments and knots of mass distribution when the gas in void regions accretes onto them. Within those nonlinear structures, internal shocks are produced by infall of previously shocked gas to filaments and knots, and during subclump mergers, as well as by chaotic flow motions. Due to the low temperature of the accreting gas, the Mach number of external shocks is high, extending up to $M\sim 100$ or higher. In contrast, internal shocks have mostly low Mach numbers. For all shocks of $M\ge1.5$ the mean distance between shock surfaces over the entire computed volume is $\sim4 h^{-1}$ Mpc at present, or $\sim 1 h^{-1}$ Mpc for internal shocks within nonlinear structures. Identified external shocks are more extensive, with their surface area $\sim2$ times larger than that of identified internal shocks at present. However, especially because of higher preshock densities, but also due to higher shock speeds, internal shocks dissipate more energy. Hence, the internal shocks are mainly responsible for gas thermalization as well as CR acceleration. In fact, internal shocks with $2 \la M \la 4$ contribute $\sim 1/2$ of the total dissipation. Using a nonlinear diffusive shock acceleration model for CR protons, we estimate the ratio of CR energy to gas thermal energy dissipated at cosmological shock waves to be $\sim1/2$ through the history of the universe. Our result supports scenarios in which the intracluster medium contains energetically significant populations of CRs.

Nonthermal particles can be produced due to incomplete thermalization at collisionless shocks and further accelerated to very high energies via diffusive shock acceleration. In a previous study we explored the cosmic ray (CR) acceleration at cosmic shocks through numerical simulations of CR modified, quasi-parallel shocks in 1D plane-parallel geometry with the physical parameters relevant for the shocks emerging in the large scale structure formation of the universe (Kang & Jones 2002). In order to consider the CR acceleration at shocks with a broader range of physical properties, in this contribution we present additional simulations. For a given Mach number, shocks with higher speeds accelerate CRs faster with a greater number of particles. However, two shocks with a same Mach number but with different shock speeds evolve qualitatively similarly when the results are presented in terms of diffusion length and time scales. Therefore, the time asymptotic value for the fraction of shock kinetic energy transferred to CRs is mainly controlled by shock Mach number rather than shock speed. We present the ``CR energy ratio'' for a wide range of shock parameters and for $\epsilon=0.2-0.3$ at terminal time of our simulations. We suggest that these values can be considered as time-asymptotic values for the CR acceleration efficiency, since the time-dependent evolution of CR modified shocks has become approximately self-similar before the terminal time.

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.

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

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.

We have performed two-dimensional hydrodynamic simulations of the collapse of isolated axisymmetric clouds condensing via radiative cooling in a primordial background gas. In order to study the development of the so-called ``shape-instability'', we have considered two types of axisymmetric clouds, oblate and prolate clouds of various sizes and with axial ratios of $0.5 \leq {R_{\rm c,R}} /{R_{\rm c,z}} \leq 2$. We find that the degree of oblateness or prolateness is enhanced during the initial cooling phase. But it can be reversed later, if the initial contrast in cooling times between the cloud gas and the background gas is much greater than one. In such cases an oblate cloud collapses to a structure composed of an outer thin disk and a central prolate component. A prolate cloud, on the other hand, becomes a thin cigar-shape structure with a central dense oblate component. The reversal of shape in the central part of the cooled clouds is due to supersonic motions either along the disk plane in the case of oblate clouds or along the symmetry axis in the case of prolate clouds. For a background gas of $T_h=1.7\times 10^6$K and $n_h=0.1 \cm3$ in a protogalactic halo environment, the mean density of the cloud gas that has cooled to $10^4$K increases to $100 n_h$ or so, in our simulations where nonequilibrium cooling is adopted and the background gas cools too. The spherical Jeans mass of such gas is estimated to be about $M_J \sim 5\times10^{7}\Msun$. In order for cloud mass to exceed the Jeans mass and at the same time in order for the thermal instability to operate, the initial cloud size should be around $1 - 1.5 l_{\rm cool}$ where $l_{\rm cool}$ is the cooling length.

A numerical scheme that incorporates a thermal leakage injection model into a combined gas dynamics and cosmic ray (CR, hereafter) diffusion-convection code has been developed. The particle injection is followed numerically by filtering the diffusive flux of suprathermal particles across the shock to the upstream region according to a velocity-dependent transparency function that controls the fraction of leaking particles. We have studied CR injection and acceleration efficiencies during the evolution of CR modified planar shocks for a wide range of initial shock Mach numbers, $M_0$, assuming a Bohm-like diffusion coefficient. The injection process is very efficient when the subshock is strong, leading to fast and significant modification of the shock structure. As the CR pressure increases, the subshock weakens and the injection rate decreases accordingly, so that the subshock does not disappear. Although some fraction of the particles injected early in the evolution continue to be accelerated to ever higher energies, the postshock CR pressure reaches an approximate time-asymptotic value due to a balance between fresh injection/acceleration and advection/diffusion of the CR particles away from the shock. We conclude that the injection rates in strong parallel shocks are sufficient to lead to rapid nonlinear modifications to the shock structures and that self-consistent injection and time-dependent simulations are crucial to understanding the non-linear evolution of CR modified shocks.

Shocks are a ubiquitous consequence of cosmic structure formation, and they play an essential role in heating galaxy cluster media. Virtually all of the gas in clusters has been processed by one or more shocks of at least moderate strength. These are collisionless shocks, so likely sites for diffusive shock acceleration of high energy particles. We have carried out numerical simulations of cosmic structure formation that directly include acceleration and transport of nonthermal protons, as well as primary and secondary electrons. Nonthermal emissions have also been computed from the resulting particle spatial and energy distributions. Here we outline some of our current findings, showing that nonthermal protons may contribute a significant pressure in cluster media, and that expected radio, X-ray and $\gamma$-ray emissions from these populations should be important cluster diagnostics.

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 <B^2>^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.

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.

Diffusive shock acceleration (DSA) is now widely accepted as the model to explain the production of cosmic rays (CRs) in a wide range of astrophysical environments. Despite initial successes of the theory in explaining the energetics and the spectrum of CRs accelerated by supernova remnants, there still remain some unresolved issues such as particle injection out of the thermal plasma at shocks, CR diffusion due to the self-generated MHD waves and yet-to-be-detected gamma-ray emission due to the ionic CRs. Recent technical advancements to resolve these issues are reviewed.

We have developed a novel computer code designed to follow the evolution of cosmic-ray modified shocks, including the full momentum dependence of the particles for a realistic diffusion coefficient model. In this form the problem is technically very difficult, because one needs to cover a wide range of diffusive scales, beginning with those slightly larger than the physical shock thickness. With most finite difference schemes for Euler's equations the numerical shock thickness is at least one zone across, so this provides a lower bound on the physical scale for diffusive transport computation. Our code uses sub-zone shock tracking and multi-level adaptive mesh refinement to provide enhanced spatial resolution around shocks at modest cost compared to the coarse grid and vastly improved cost effectiveness compared to a uniform, highly refined grid. We present and discuss the implications from our initial results.

One of the key questions to understanding the efficiency of diffusive shock acceleration of the cosmic rays (CRs) is the injection process from thermal particles. A self-consistent injection model based on the interactions of the suprathermal particles with self-generated magneto-hydrodynamic waves has been developed recently by Malkov (1998, Phys. Rev. E 58, 4911). By adopting this analytic solution, a numerical treatment of the plasma-physical injection model at a strong quasi-parallel shock has been devised and incorporated into the combined gas dynamics and the CR diffusion-convection code. In order to investigate self-consistently the injection and acceleration efficiencies, we have applied this code to the CR modified shocks of both high and low Mach numbers (M=30 and M=2.24) with a Bohm type diffusion model. We find the injection process is self-regulated in such a way that the injection rate reaches and stays at a nearly stable value after quick initial adjustment. For both shocks about 0.1% of the incoming thermal particles are injected into the CRs. For the weak shock, the shock has reached a steady state within our integration time and about 10% of the total available shock energy is transfered into the CR energy density. The strong shock has achieved a higher acceleration efficiency of about 20% by the end of our simulation, but has not yet reached a steady-state.

Many models of globular cluster formation assume the presence of cold dense clouds in early universe. Here we re-examine the Fall & Rees (1985) model for formation of proto-globular cluster clouds (PGCCs) via thermal instabilities in a protogalactic halo. We first argue, based on the previous study by others, that under the protogalactic environments only nonlinear density inhomogeneities can condense into PGCCs. We then carry out numerical simulations of the collapse of overdense clouds in one-dimensional spherical geometry, including self-gravity and radiative cooling down to T=10^4 K. Since imprinting of Jeans mass at 10^4 K is essential to this model, here we focus on the cases where external UV background radiation prevents the formation of H2 molecules and so prevent the cloud from cooling below 10^4 K. The quantitative results from these simulations can be summarized as follows: 1) Perturbations smaller than M_min ~ (10^{5.6} M_sun) (n_h/0.05 cm3)^{-2} cool isobarically, while perturbations larger than M_max ~ (10^8 M_sun) (n_h/0.05 cm3)^{-2} cool isochorically. On the other hand, intermediate size perturbations (M_min< M_pgcc < M_max) are compressed supersonically. 2) For supersonically collapsing clouds, the density compression factor after they cool to T_c=10^4 K range 10^{2.5}-10^6. 3) For supersonically collapsing clouds the Jeans mass can be reduced to as small as 10^{5.5} M_sun (n_h/0.05 cm3)^{-1/2} at the maximum compression. 4) The density profile of simulated PGCCs can be approximated by a constant core with a halo of rho ~ r^{-2} rather than a singular isothermal sphere.

We review the possible roles of large scale shocks as particle accelerators in clusters of galaxies. Recent observational and theoretical work has suggested that high energy charged particles may constitute a substantial pressure component in clusters. If true that would alter the expected dynamical evolution of clusters and increase the dynamical masses consistent with hydrostatic equilibrium. Moderately strong shocks are probably common in clusters, through the actions of several agents. The most obvious of these agents include winds from galaxies undergoing intense episodes of starbursts, active galaxies and cosmic inflows, such as accretion and cluster mergers. We describe our own work derived from simulations of large scale structure formation, in which we have, for the first time, explicitly included passive components of high energy particles. We find, indeed that shocks associated with these large scale flows can lead to nonthermal particle pressures big enough to influence cluster dynamics. These same simulations allow us also to compute nonthermal emissions from the clusters. Here we present resulting predictions of gamma-ray fluxes.

We investigate the dynamical importance of a newly recognized possible source of significant feedback generated during structure formation; namely cosmic ray (CR) pressure. We present evidence for the existence of numerous shocks in the hot gas of galaxy clusters (GCs). We employ for the first time an explicit numerical treatment of CR acceleration and transport in hydro simulations of structure formation. According to our results, CRs provide an important fraction of the total pressure inside GCs, up to several tenths. This was true even at high redshift (z=2), meaning that such non-thermal component could affect the evolution of 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.

A cosmological hydrodynamic code is described, which includes a routine to compute cosmic ray acceleration and transport in a simplified way. The routine was designed to follow explicitly diffusive acceleration at shocks, and second-order Fermi acceleration and adiabatic loss in smooth flows. Synchrotron cooling of the electron population can also be followed. The updated code is intended to be used to study the properties of nonthermal synchrotron emission and inverse Compton scattering from electron cosmic rays in clusters of galaxies, in addition to the properties of thermal bremsstrahlung emission from hot gas. The results of a test simulation using a grid of $128^3$ cells are presented, where cosmic rays and magnetic field have been treated passively and synchrotron cooling of cosmic ray electrons has not been included.

The properties of cosmic shock waves are studied through numerical simulations in two cosmological scenarios (SCDM and LCDM). The scaling relations for the average radius and velocity associated with the accretion shocks are somewhat different, yet qualitatively similar to the self similar solutions for a flat Omega_M=1 universe. The energy supplied by infalling gas at accretion shock waves is large enough to sustain production of abundant cosmic ray populations if a viable acceleration mechanism can take place there. Finally, in addition to shocks created by the encounter of the merging ICMs of two clusters of galaxies, accretion shocks associated with the merging clusters generate strong ``relic'' shocks which propagate through the ICM producing additional heating of the ICM, and associated CR acceleration.

The theory of diffusive particle acceleration explains the spectral properties of the cosmic rays below energies of approx. 10^6 GeV as produced at strong shocks in supernova remnants (SNR's). To supply the observed flux of cosmic rays, a significant fraction of the energy released by a supernova has to be transfered to cosmic rays. The key to the question of the efficiency of SNR's in producing cosmic rays is the injection process from thermal energies. A self-consistent model has to take into account the interaction of the accelerated particles with magneto-hydrodynamic waves, which generate the particle diffusion, a requisite for the acceleration process. Such a nonlinear model of the turbulent background plasma has been developed recently (Malkov, 1998, Phys. Rev. E 58, 4911). We use this model for the first numerical treatment of the gas dynamics and the diffusion-convection equation at a quasi-parallel strong shock, which incorporates a plasma-physical injection model to investigate the cosmic ray production.

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 \mu 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.

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

The recent discoveries of several reliable events of high energy cosmic rays at an energy above 10^20 eV raise questions about their path through the nearby universe. The two analyses of the Haverah Park data set and the Akeno data set appear to have an inconsistent pattern. Both data sets showed some measure of a correlation with the supergalactic plane, the locus of cosmologically nearby galaxies, radio galaxies and clusters of galaxies. We need a reasonable model of the true intergalactic magnetic field. Using recent cosmological simulations of structure formation in the universe, we estimate the magnetic fields which correspond to the upper limits in the Rotation Measure to distant radio sources. We estimate that the magnetic field strength in supergalactic sheets and filaments may be in the range of 0.1 to 1 muG. If such strengths are realized inside our Local Supercluster, this opens up the possibility to focus charged particles in the direction perpendicular to the supergalatic plane, analoguously but in the opposite direction to solar wind modulation. If focusing exists, it means that for all particles captured into the sheets, the dilution with distance d is 1/d instead of 1/d^2. This means in effect, that we may see sources to much larger distances than expected sofar. This effect is relevant only for energies for which the possible distances are smaller than the void scale of the cosmological galaxy distribution, in the range possibly up to 100 Mpc.

The simulations of the formation of cosmological structure allows to determine the spatial inhomogeneity of cosmic magnetic fields. Such simulations, however, do not give an absolute number for the strength of the magnetic field due to insufficient spatial resolution. Combining these simulations with observations of the Rotation Measure to distant radio sources allows then to deduce upper limits for the strength of the magnetic field. These upper limits are of order 0.2 - 2 muG along the filaments and sheets of the galaxy distribution. In one case, the sheet outside the Coma cluster, there is a definitive estimate of the strength of the magnetic field consistent with this range. Such estimates are almost three orders of magnitude higher than hitherto assumed usually. High energy cosmic ray particles can be either focussed or strongly scattered in such magnetic filaments and sheets, depending on the initial transverse momentum. The cosmological background in radio and X-ray wavelengths will have contributions from these intergalactic filaments and sheets, should the magnetic fields really be as high as 0.2 - 2 muG.

In the standard theory of the large scale structure formation, matter accretes onto high density perturbations via gravitational instability. Collisionless dark matter forms caustics around such structures, while collisional baryonic matter forms accretion shocks which then halt and heat the infalling gas. Here, we discuss the characteristics, roles, and observational consequences of these accretion shocks.

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

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 $\Omega_o<1$, and the flat universe with $\Omega_{\Lambda}=1-\Omega_o$. A density parameter in the range $0.1\le \Omega_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 $\Omega_o$, but only by up to $\sim24\%$ in the models with $0.1\le \Omega_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 $\Omega_o$ by up to $\sim41\%$ and $\sim65\%$, respectively. So, it can provide a better

We report simulations of diffusive particle acceleration in oblique magnetohydrodynamical (MHD) shocks. These calculations are based on extension to oblique shocks of a numerical model for ``thermal leakage'' injection of particles at low energy into the cosmic-ray population. That technique, incorporated into a fully dynamical diffusion-convection formalism, was recently introduced for parallel shocks by Kang \& Jones (1995). Here, we have compared results of time dependent numerical simulations using our technique with Monte Carlo simulations by Ellison, Baring \& Jones 1995 and with {\it in situ} observations from the Ulysses spacecraft of oblique interplanetary shocks discussed by Baring \etal (1995). Through the success of these comparisons we have demonstrated that our {diffusion-convection} method and injection techniques provide a practical tool to capture essential physics of the injection process and particle acceleration at oblique MHD shocks. In addition to the diffusion-convection simulations, we have included time dependent two-fluid simulations for a couple of the shocks to demonstrate the basic validity of that formalism in the oblique shock context. Using simple models for the two-fluid closure parameters based on test-particle considerations, we find good agreement with the dynamical properties of the more detailed diffusion-convection results. We emphasize, however, that such two-fluid results can be sensitive to the properties of these closure parameters when the flows are not truly steady. Furthermore, we emphasize through example how the validity of the two-fluid formalism does not necessarily mean that {\it steady-state} two-fluid models provide a reliable tool for predicting the efficiency of particle acceleration in real shocks.

We have examined consequences of strong tidal encounters between a neutron star and a normal star using SPH as a possible formation mechanism of isolated recycled pulsars in globular clusters. We have made a number of SPH simulations for close encounters between a main-sequence star of mass ranging from 0.2 to 0.7 Solar masses represented by an n=3/2 polytrope and a neutron star represented by a point mass. The outcomes of the first encounters are found to be dependent only on the dimensionless parameter eta' = (m/(m+M))^(1/2) (r_min/R_MS)^(3/2) (m/M)^(1/6), where m and M are the masses of the main-sequence star and the neutron star, respectively, r_min the minimum separation between two stars, and R_MS the size of the main-sequence star. The material from the (at least partially) disrupted star forms a disk around the neutron star. If all material in the disk is to be acctreted onto the neutron star's surface, the mass of the disk is enough to spin up the neutron star to spin period of 1 ms.

A star around a massive black hole can be disrupted tidally by the gravity of the black hole. Then, its debris may form a precessing stream which may even collide with itself. In order to understand the dynamical effects of the stream-stream collision on the eventual accretion of the stellar debris onto the black hole, we have studied how gas flow behaves when the outgoing stream collides supersonically with the incoming stream. We have investigated the problem analytically with one-dimensional plane-parallel streams and numerically with more realistic three-dimensional streams. A shock formed around the contact surface converts the bulk of the orbital streaming kinetic energy into thermal energy. In three-dimensional simulations, the accumulated hot post-shock gas then expands adiabatically and drives another shock into the low density ambient region. Through this expansion, thermal energy is converted back to the kinetic energy associated with the expanding motion. Thus, in the end, only a small fraction of the orbital kinetic energy is actually converted to the thermal energy, while most of it is transferred to the kinetic energy of the expanding gas. Nevertheless the collision is effective in circularizing the debris orbit, because the shock efficiently transforms the ordered motion of the streams into the expanding motion in directions perpendicular to the streams. The circularization efficiency decreases, if two colliding streams have a large ratio of cross sections and a large density contrast. But even in such cases, the main shock extends beyond the overlapping contact surface and the high pressure region behind the shock keeps the stream of the larger cross section from passing freely. Thus the stream-stream collisions are still expected to circularize the stellar debris rather efficiently, unless the ratio

Three-dimensional hydrodynamic simulations of large scale structure in the Universe have shown that accretion shocks form during the gravitational collapse of one-dimensional caustics, and that clusters of galaxies formed at intersections of the caustics are surrounded by these accretion shocks. Estimated speed and curvature radius of the shocks are 1000-3000 \kms and about 5 Mpc, respectively, in the $\Omega=1$ CDM universe. Assuming that energetic protons are accelerated by these accretion shocks via the first-order Fermi process and modeling particle transport around the shocks through Bohm diffusion, we suggest that protons can be accelerated up to the {\it Greisen cutoff energy} near $6\times 10^{19}$ eV, provided the mean magnetic field strength in the region around the shocks is at least of order a microgauss. We have also estimated the proton flux at earth from the Virgo cluster. Assuming a few (1-10) \% of the ram pressure of the infalling matter would be transferred to the cosmic-rays, the estimated flux for $ E \sim 10^{19}$eV is consistent with observations, so that such clusters could be plausible sources of the UHE CRs.

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.

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

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.

We exmaine the distribution of hot gas in a standard CDM model of the universe using high resolution hydrodynamic simulations. Adopting standard parameters determined from COBE and light element nucleosynthesis, $\sigma_8=1.05$, $\Omega_b=0.06$ and assuming $h=0.5$, we find the X-ray emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution and estimated sizes as well as the evolution of these quantities with redshift. This standard CDM model, normalized to COBE, produces approximately 5 times too much emission from clusters having $L_x>10^{43}$erg/s, a not unexpected result. If all other parameters were unchanged, we would expect adequate agreement for $\sigma_8=0.6$. This provides a new and independent argument for lower small scale power than standard CDM at the $8h^{-1}$Mpc scale. The background radiation field at 1keV due to clusters in this model is approximately $1/3$ of the observed background which, after correction for numerical effects, again indicates approximately 5 times too much emission and the appropriateness of $\sigma_8=0.6$. If we had used the observed ratio of gas to total mass in clusters, rather than basing the mean density on light element nucleosynthesis, then the computed luminosity of each cluster would have increased still further, by a factor of approximately ten. Examining the ratio of gas to total mass in the clusters (which we find to be anti-biased by a factor of approximately 0.6), normalized to $\Omega_b h^2=0.015$, and comparing to observations, we conclude, in agreement with S. White, that the cluster observations argue for an open universe.