Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.

We propose a mechanism for the generation of a magnetic field in the early universe during the QCD crossover assuming that dark matter is made of axions. Thermoelectric fields arise at pressure gradients in the primordial plasma due to the difference in charge, energy density and equation of state between the quark and lepton components. The axion field is coupled to the EM field, so when its spatial gradient is misaligned with the thermoelectric field, an electric current is driven. Due to the finite resistivity of the plasma an electric field appears that is generally rotational. For a QCD axion mass consistent with observational constraints and a conventional efficiency for turbulent dynamo amplification --- driven by the same pressure gradients responsible for the thermoelectric fields --- a magnetic field is generated on subhorizon scales. After significant Alfv\'enic unwinding it reaches a present day strength of $B \sim 10^{-13 }$ G on a characteristic scale $L_B \sim $ 20 pc. The resulting combination of $BL_B^{1/2}$ is significantly stronger than in any astrophysical scenario, providing a clear test for the cosmological origin of the field through $\gamma$-ray observations of distant blazars. The amplitude of the pressure gradients may be inferred from the detection of concomitant gravitational waves, while several experiments are underway to confirm or rule out the existence of axions.

The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and study the turbulent dynamo mechanism in the laboratory. Here we describe the design of these experiments through simulation campaigns using FLASH, a highly capable radiation magnetohydrodynamics code that we have developed, and large-scale three-dimensional simulations on the Mira supercomputer at Argonne National Laboratory. The simulation results indicate that the experimental platform may be capable of reaching a turbulent plasma state and study dynamo amplification. We validate and compare our numerical results with a small subset of experimental data using synthetic diagnostics.

Magnetic fields are ubiquitous in the Universe. Extragalactic disks, halos and clusters have consistently been shown, via diffuse radio-synchrotron emission and Faraday rotation measurements, to exhibit magnetic field strengths ranging from a few nG to tens of $\mu$G. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are embedded, making magnetic fields essential players in the dynamics of the luminous matter. The standard theoretical model for the origin of these strong magnetic fields is through the amplification of tiny seed fields via turbulent dynamo to the level consistent with current observations. Here we demonstrate, using laser-produced colliding plasma flows, that turbulence is indeed capable of rapidly amplifying seed fields to near equipartition with the turbulent fluid motions. These results support the notion that turbulent dynamo is a viable mechanism responsible for the observed present-day magnetization of the Universe.

We introduce CRASH-AMR, a new version of the cosmological Radiative Transfer (RT) code CRASH, enabled to use refined grids. This new feature allows us to attain higher resolution in our RT simulations and thus to describe more accurately ionisation and temperature patterns in high density regions. We have tested CRASH-AMR by simulating the evolution of an ionised region produced by a single source embedded in gas at constant density, as well as by a more realistic configuration of multiple sources in an inhomogeneous density field. While we find an excellent agreement with the previous version of CRASH when the AMR feature is disabled, showing that no numerical artifact has been introduced in CRASH-AMR, when additional refinement levels are used the code can simulate more accurately the physics of ionised gas in high density regions. This result has been attained at no computational loss, as RT simulations on AMR grids with maximum resolution equivalent to that of a uniform cartesian grid can be run with a gain of up to 60% in computational time.

There is evidence that magnetized material along the line of sight to distant quasars is detectable in the polarization properties of the background sources. The polarization properties appear to be correlated with the presence of intervening MgII absorption, which is thought to arise in outflowing material from star forming galaxies. In order to investigate this further, we have obtained high spectral resolution polarization measurements, with the VLA and ATCA, of a set of 49 unresolved quasars for which we have high quality optical spectra. These enable us to produce a Faraday Depth spectrum for each source, using Rotation Measure Synthesis. Our new independent radio data confirms that interveners are strongly associated with depolarization. We characterize the complexity of the Faraday Depth spectrum using a number of parameters and show how these are related, or not, to the depolarization and to the presence of MgII absorption along the line of sight. We argue that complexity and structure in the Faraday Depth distribution likely arise from both intervening material and intrinsically to the background source and attempt to separate these. We find that the strong radio depolarization effects associated with intervening material at redshifts out to $z \approx 1$ arise from inhomogeneous Faraday screens producing a dispersion in Rotation Measure across individual sources of around 10~rad/m$^2$. This is likely produced by disordered fields with strengths of at least $3\;\mu$G.

Massive galaxy clusters are filled with a hot, turbulent and magnetized intra-cluster medium. Still forming under the action of gravitational instability, they grow in mass by accretion of supersonic flows. These flows partially dissipate into heat through a complex network of large-scale shocks [1], while residual transonic flows create giant turbulent eddies and cascades [2,3]. Turbulence heats the intra-cluster medium [4] and also amplifies magnetic energy by way of dynamo action [5-8]. However, the pattern regulating the transformation of gravitational energy into kinetic, thermal, turbulent and magnetic energies remains unknown. Here we report that the energy components of the intra-cluster medium are ordered according to a permanent hierarchy, in which the ratio of thermal to turbulent to magnetic energy densities remains virtually unaltered throughout the cluster's history, despite evolution of each individual component and the drive towards equipartition of the turbulent dynamo. This result revolves around the approximately constant efficiency of turbulence generation from the gravitational energy that is freed during mass accretion, revealed by our computational model of cosmological structure formation [3,9]. The permanent character of this hierarchy reflects yet another type of self-similarity in cosmology [10-13], while its structure, consistent with current data [14-18], encodes information about the efficiency of turbulent heating and dynamo action.

We compare DNS calculations of homogeneous isotropic turbulence with the statistical properties of intra-cluster turbulence from the Matryoshka Run (Miniati 2014) and find remarkable similarities between their inertial ranges. This allowed us to use the time dependent statistical properties of intra-cluster turbulence to evaluate dynamo action in the intra-cluster medium, based on earlier results from numerically resolved nonlinear magneto-hydrodynamic turbulent dynamo (Beresnyak 2012). We argue that this approach is necessary (a) to properly normalize dynamo action to the available intra-cluster turbulent energy and (b) to overcome the limitations of low Re affecting current numerical models of the intra-cluster medium. We find that while the properties of intra-cluster magnetic field are largely insensitive to the value and origin of the seed field, the resulting values for the Alfven speed and the outer scale of the magnetic field are consistent with current observational estimates, basically confirming the idea that magnetic field in today's galaxy clusters is a record of its past turbulent activity.

We use the Zurich ENvironmental Study (ZENS) database to investigate the environmental dependence of the merger fraction $\Gamma$ and merging galaxy properties in a sample of ~1300 group galaxies with $M>10^{9.2}M_\odot$ and 0.05<z<0.0585. In all galaxy mass bins investigated in our study, we find that $\Gamma$ decreases by a factor of ~2-3 in groups with halo masses $M_{HALO}>10^{13.5} M_\odot$ relative to less massive systems, indicating a suppression of merger activity in large potential wells. In the fiducial case of relaxed groups only, we measure a variation $\Delta\Gamma/\Delta \log (M_{HALO}) \sim - 0.07$ dex$^{-1}$, which is almost independent of galaxy mass and merger stage. At galaxy masses $>10^{10.2} M_\odot$, most mergers are dry accretions of quenched satellites onto quenched centrals, leading to a strong increase of $\Gamma$ with decreasing group-centric distance at these mass scales.Both satellite and central galaxies in these high mass mergers do not differ in color and structural properties from a control sample of nonmerging galaxies of equal mass and rank. At galaxy masses $<10^{10.2} M_\odot$, where we mostly probe satellite-satellite pairs and mergers between star-forming systems, close pairs (projected distance $<10-20$ kpc) show instead $\sim2\times$ enhanced (specific) star formation rates and $\sim1.5\times$ larger sizes than similar mass, nonmerging satellites. The increase in both size and SFR leads to similar surface star-formation densities in the merging and control-sample satellite populations.

We use the Matryoshka run to study the time dependent statistics of structure-formation driven turbulence in the intracluster medium of a 10$^{15}M_\odot$ galaxy cluster. We investigate the turbulent cascade in the inner Mpc for both compressional and incompressible velocity components. The flow maintains approximate conditions of fully developed turbulence, with departures thereof settling in about an eddy-turnover-time. Turbulent velocity dispersion remains above $700$ km s$^{-1}$ even at low mass accretion rate, with the fraction of compressional energy between 10% and 40%. Normalisation and slope of compressional turbulence is susceptible to large variations on short time scales, unlike the incompressible counterpart. A major merger occurs around redshift $z\simeq0$ and is accompanied by a long period of enhanced turbulence, ascribed to temporal clustering of mass accretion related to spatial clustering of matter. We test models of stochastic acceleration by compressional modes for the origin of diffuse radio emission in galaxy clusters. The turbulence simulation model constrains an important unknown of this complex problem and brings forth its dependence on the elusive micro-physics of the intracluster plasma. In particular, the specifics of the plasma collisionality and the dissipation physics of weak shocks affect the cascade of compressional modes with strong impact on the acceleration rates. In this context radio halos emerge as complex phenomena in which a hierarchy of processes acting on progressively smaller scales are at work. Stochastic acceleration by compressional modes implies statistical correlation of radio power and spectral index with merging cores distance, both testable in principle with radio surveys.

We present the results of a pilot XMM-$Newton$ and $Chandra$ program aimed at studying the diffuse intragroup medium (DIM) of optically-selected nearby groups from the Zurinch ENvironmental Study (ZENS) catalog. The groups are in a narrow mass range about $10^{13}M_\odot$, a mass scale at which the interplay between the DIM and the group member galaxies is still largely unprobed. X-ray emission from the DIM is detected in the energy band 0.5--2 keV with flux $\le 10^{-14}$ erg cm$^{-1}$ s$^{-1}$, which is one order of magnitude fainter than for typical ROSAT groups (RASS). For many groups we set upper limits to the X-ray luminosity, indicating that the detections are likely probing the upper envelope of the X-ray emitting groups. We find evidence for our optically selected groups to be under-luminous with respect to predictions from X-ray scaling relations. X-ray mass determinations are in best agreement with those based on the member galaxies bulge luminosity, followed by their total optical luminosity and velocity dispersion. We measure a stellar mass fraction with a median value of about 1$\%$, with a contribution from the most massive galaxies between 30 to 50 \%. Optical and X-ray data give often complementary answers concerning the dynamical state of the groups, and are essential for a complete picture of the system. Extending this pilot program to a larger sample of groups is mandatory to unveil any imprint of interaction between member galaxies and DIM in halo potentials of key importance for environmentally-driven galactic evolution.

We use ZENS data at low redshift to study the dependence of the quenched satellite fraction and of the morphological mix of these quenched satellites, on three different environmental parameters: group halo mass, halo-centric distance and large-scale structure over-density. The fraction of quenched satellites is independent of halo mass and the surrounding large-scale overdensity, but increases towards the centres of the haloes, as previously found. The morphological mix is, however, constant with radial position, indicating that the well-known morphology-density relation results from the increasing fraction of quenched galaxies towards the centres of haloes. The constancy of the morphological outcome suggests that mass-quenching and satellite quenching have the same effect on the morphologies of the galaxies. The quenched satellites have larger B/T and smaller half-light radii than the star-forming satellites. These are mostly due to differences in the disks. The bulges in quenched satellites have very similar luminosities and surface brightness profiles, and any mass growth of the bulges associated with quenching cannot greatly change these quantities. The quenched disks are fainter and have smaller scale lengths than in star-forming satellites. This can be explained either by a differential fading of the disks or if disks were generally smaller in the past, both of which are expected in an inside-out growth of disks. At least at low redshifts, the structure of massive quenched satellites is produced by processes that operate before quenching takes place. A comparison with semi-analytic models argues for a reduction in the efficiency of group halos in quenching their disk satellites and for mechanisms to increase the B/T of low mass quenched satellites. [abridged]

We study the statistical properties of turbulence driven by structure formation in a massive merging galaxy cluster at z=0. The development of turbulence is ensured as the largest eddy turnover time is much shorter than the Hubble time independent of mass and redshift. We achieve a large dynamic range of spatial scales through a novel Eulerian refinement strategy where the cluster volume is refined with progressively finer uniform nested grids during gravitational collapse. This provides an unprecedented resolution of 7.3 h^{-1} kpc across the virial volume. The probability density functions of various velocity derived quantities exhibit the features characteristic of fully developed compressible turbulence observed in dedicated periodic-box simulations. Shocks generate only 60% of the total vorticity within \rvir/3 region and 40% beyond that. We compute second and third order, longitudinal and transverse, structure functions for both solenoidal and compressional components, in the cluster core, virial region and beyond. The structure functions exhibit a well defined inertial range. The injection scale is comparable to the virial radius but increases towards the outskirts. Within \rvir/3, the spectral slope of the solenoidal component is close to Kolmogorov's, but for the compressional component is substantially steeper and close to Burgers'; the flow is mostly solenoidal and statistically rigorously consistent with fully developed, homogeneous and isotropic turbulence. Small scale anisotropy appears due to numerical artifact. Towards the virial region, the flow becomes compressional, the structure functions flatter and modest genuine anisotropy appear. In comparison, mesh adaptivity based on Lagrangian refinement and the same finest resolution, leads to lack of turbulent power on small scale and excess thereof on large scales, and unreliable density weighted structure functions.

We present photometric measurements for the galaxies - and when possible their bulges and disks - in the 0.05<z<0.0585 groups of the Zurich Environmental Study (ZENS); these measurements include (B-I) colors, color gradients and maps, color dispersions, as well as stellar masses and star-formation rates. The ZENS galaxies are classified into quenched, moderately star-forming, and strongly star-forming using a combination of spectral features and FUV-to-optical colors; this approach optimally distinguishes quenched systems from dust-reddened star-forming galaxies. The latter contribute up to 50% to the (B-I) "red sequence" at ~10^10Msun. At fixed morphological or spectral type, we find that galaxy stellar masses are largely independent of environment, and especially of halo mass. As a first utilization of our photometric database, we study, at fixed stellar mass and Hubble type, how (B-I) colors, color gradients and color dispersion of disk satellites depend on group mass (M_GROUP), group-centric distance (R/R_200) and large-scale structure overdensity. The strongest environmental trend is found for disk-dominated satellites with M_GROUP and R/R_200. At M<10^10 Msun, disk-dominated satellites are redder in the inner regions of the groups than in the outer parts. At M>10^10 Msun, these satellites have shallower color gradients in higher mass groups and in the cores of groups compared with lower mass groups and the outskirts of groups. Stellar population analyses and semi-analytic models suggest that disk-dominated satellites undergo quenching of star formation in their outer disks, on timescales ~2 Gyr, as they progressively move inside the group potential.

We describe X-ray observations with Chandra and XMM-Newton of 18 galaxy groups (M_group ~ 1-6x10^13 Msolar, z~0.05) from the Zurich Environmental Study (ZENS). We aim to establish the frequency and properties, unaffected by host galaxy dilution and obscuration, of AGNs in central and satellite galaxy members, also as a function of halo-centric distance. X-ray point-source detections are reported for 22 of 177 observed galaxies, down to a limit of f_(0.5-8 keV) ~ 5x10^-15 erg cm^-2 s^-1, corresponding to a limiting luminosity of L_(0.5-8 keV)~3x10^40 erg s^-1. With the majority of the X-ray sources attributed to AGNs of low-to-moderate levels (L/L_Edd>~10^-4), we discuss the detection rate in the context of the occupation of AGNs to halos of this mass scale and redshift, and compare the structural/morphological properties between AGN-active and non-active galaxies of different rank and location within the group halos. We see a slight tendency for AGN hosts to have either relatively brighter/denser disks (or relatively fainter/diffuse bulges) than non-active galaxies of similar mass. At galaxy mass scales <10^11 Msolar, central galaxies appear to be a factor ~4 more likely to host AGNs than satellite galaxies of similar mass. This effect, coupled with the tendency for AGNs to reside in massive galaxies, explains the (weak) trend for AGNs to be preferentially found in the inner regions of groups, with no detectable trend with halo-centric distance in the frequency of AGNs within the satellite population. Finally, our data support other analyses in finding that the rate of decline with redshift of AGN activity in groups matches that of the global AGN population, indicating that either AGNs occur preferentially in groups, or that the evolution rate is independent of halo mass. These trends are of potential importance, and require X-ray coverage of a larger sample to be solidly confirmed.

We present structural measurements for the galaxies in the 0.05<z<0.0585 groups of the Zurich Environmental Study, aimed at establishing how galaxy properties depend on four environmental parameters: group halo mass M_GROUP, group-centric distance R/R_200, ranking into central or satellite, and large-scale structure density delta_LSS. Global galaxy structure is quantified both parametrically and non-parametrically. We correct all these measurements for observational biases due to PSF blurring and surface brightness effects as a function of galaxy size, magnitude, steepness of light profile and ellipticity. Structural parameters are derived also for bulges, disks and bars. We use the galaxy bulge-to-total ratios (B/T), together with the calibrated non-parametric structural estimators, to implement a quantitative morphological classification that maximizes purity in the resulting morphological samples. We investigate how the concentration C of satellite galaxies depends on galaxy mass for each Hubble type, and on M_GROUP, R/R_200 and delta_LSS. At galaxy masses M>10^10 M_sun, the concentration of disk satellites increases with increasing stellar mass, separately within each morphological bin of B/T. The known increase in concentration with stellar mass for disk satellites is thus due, at least in part, to an increase in galaxy central stellar density at constant B/T. The correlation between concentration and galaxy stellar mass becomes progressively steeper for later morphological types. The concentration of disk satellites shows a barely significant dependence on delta_LSS or R/R_200. The strongest environmental effect is found with group mass for M>10^10 M_sun disk-dominated satellites, which are ~10% more concentrated in high mass groups than in lower mass groups.

The Zurich Environmental Study (ZENS) is based on a sample of ~1500 galaxy members of 141 groups in the mass range ~10^12.5-14.5 M_sun within the narrow redshift range 0.05<z<0.0585. ZENS adopts novel approaches, here described, to quantify four different galactic environments, namely: (1) the mass of the host group halo; (2) the projected halo-centric distance; (3) the rank of galaxies as central or satellites within their group halos; and (4) the filamentary large-scale structure (LSS) density. No self-consistent identification of a central galaxy is found in ~40% of <10^13.5 M_sun groups, from which we estimate that ~15% of groups at these masses are dynamically unrelaxed systems. Central galaxies in relaxed and unrelaxed groups have in general similar properties, suggesting that centrals are regulated by their mass and not by their environment. Centrals in relaxed groups have however ~30% larger sizes than in unrelaxed groups, possibly due accretion of small satellites in virialized group halos. At M>10^10 M_sun, satellite galaxies in relaxed and unrelaxed groups have similar size, color and (specific) star formation rate distributions; at lower galaxy masses, satellites are marginally redder in relaxed relative to unrelaxed groups, suggesting quenching of star formation in low-mass satellites by physical processes active in relaxed halos. Finally, relaxed and unrelated groups show similar stellar mass conversion efficiencies, peaking at halo masses around 10^12.5 M_sun. In the enclosed ZENS catalogue we publish all environmental diagnostics as well as the galaxy structural and photometric measurements described in companion ZENS papers II and III.

Radio quasar sightlines with strong MgII absorption lines display statistically enhanced Faraday Rotation Measures (RM) indicating the presence of additional magneto-active plasma with respect to sightlines free of such absorption. In this letter, we use multi-color optical imaging to identify the likely galaxies hosting the magneto-active plasma, and to constrain the location of the latter with respect to the putative parent halo. We find that all of the sightlines with high |RM| pass within 50 kpc of a galaxy and that the |RM| distribution for low impact parameters, D < 50 kpc, is significantly different than for larger impact parameters. In addition, we find a decrease in the ratio of the polarization at 21 cm and 1.5 cm, p_{21}/p_{1.5}, towards lower D. These two effects are most likely related, strengthen the association of excess |RM| with intervening galaxies, and suggest that intervening galaxies operate as inhomogeneous Faraday screens. These results are difficult to reconciliate with only a disk model for the magnetic field but are consistent with highly magnetized winds associated with MgII systems. We infer strong magnetic fields of a few tens of mikro Gauss, consistent with values required by the lack of evolution of the FIR-radio correlation at high redshifts. Finally, these findings lends support to the idea that the small scale helicity bottleneck of alpha-Omega galactic dynamos can be significantly alleviated via galatic winds.

The stability properties of a low density ultra relativistic pair beam produced in the intergalactic medium by multi-TeV gamma-ray photons from blazars are analyzed. The problem is relevant for probes of magnetic field in cosmic voids through gamma-ray observations. In addition, dissipation of such beams could affect considerably the thermal history of the intergalactic medium and structure formation. We use a Monte Carlo method to quantify the properties of the blazar induced electromagnetic shower, in particular the bulk Lorentz factor and the angular spread of the pair beam generated by the shower, as a function of distance from the blazar itself. We then use linear and nonlinear kinetic theory to study the stability of the pair beam against the growth of electrostatic plasma waves, employing the Monte Carlo results for our quantitative estimates. We find that the fastest growing mode, like any perturbation mode with even a very modest component perpendicular to the beam direction, cannot be described in the reactive regime. Due to the effect of non-linear Landau damping, which suppresses the growth of plasma oscillations, the beam relaxation timescale is found significantly longer than the inverse Compton loss time. Finally, density inhomogeneities associated with cosmic structure induce loss of resonance between the beam particles and plasma oscillations, strongly inhibiting their growth. We conclude that relativistic pair beams produced by blazars in the intergalactic medium are stable on timescales long compared to the electromagnetic cascade's. There appears to be little or no effect of pair-beams on the intergalactic medium.

We analyze the redshift evolution of the Rotation Measure (RM) in Taylor et al. (2009) dataset, which is based on NVSS radio data at 21 cm, and compare with results from our previous work (Kronberg et al. (2008), Bernet et al. (2008),(2010)), based on RMs determined at lower wavelengths, e.g. 6 cm. We find that, in spite of the same analysis, Taylor's dataset produces neither an increase of the RM dispersion with redshift as found in Kronberg et al. (2008), nor the correlation of RM strength with MgII absorption lines found in Bernet et al. (2008). We develop a simple model to understand the discrepancy. The model assumes that the Faraday Rotators, namely the QSO's host galaxy and the intervening MgII host galaxies along the line of sight, contain partially inhomogeneous RM screens. We find that this leads to an increasing depolarization towards longer wavelengths and to wavelength dependent RM values. In particular, due to cosmological redshift, observations at fixed wavelength of sources at different redshift are affected differently by depolarization and are sensitive to different Faraday active components. For example, at 21 cm the polarized signal is averaged out by inhomogeneous Faraday screens and the measured RM mostly reflects the Milky Way contributions for low redshift QSOs, while polarization is relatively unaffected for high redshift QSOs. Similar effects are produced by intervening galaxies acting as inhomogeneous screens. Finally, we assess the performance of Rotation Measure synthesis on our synthetic models and conclude that the study of magnetic fields in galaxies as a function of cosmic time will benefit considerably from the application of such a technique, provided enough instrumental bandwidth. For this purpose, high frequency channels appear preferable but not strictly necessary.

Miniati and Bell (2011) proposed a mechanism for the generation of magnetic seeds that is based the finite resistivity of the low temperature IGM in the high redshift universe. In this model, cosmic-ray protons generated by the first generation of galaxies, escape into the intergalactic medium carrying an electric current that induces return currents, $j_t$, and associated electric fields, $\vec E=\eta\vec j_t$ there. Because the resistivity, $\eta$, depends on the IGM temperature, which is highly inhomogeneous due to adiabatic contraction and shocks produced by structure formation, a non-vanishing curl of the electric field exists which sustains the growth of magnetic field. In this contribution we have developed an approximate numerical model for this process by implementing the source terms of the resistive mechanism in the cosmological code CHARM. Our numerical estimates substantiate the earlier analysis in Miniati and Bell (2011) which found magnetic seeds between 10$^{-18}$ and 10$^{-16}$ Gauss throughout cosmic space at redshift z~6, consistent with conservative estimates of magnetic fields in voids at z~0 from recent gamma-ray experiments.

We discuss the amplification of magnetic fields by the small-scale dynamo, a process that could efficiently produce strong magnetic fields in the first galaxies. In addition, we derive constraints on the primordial field strength from the epoch of reionization.

Observations of giant radio halos provide unambiguous evidence for the existence of cosmic ray (CR) electrons and magnetic fields in galaxy clusters. The physical mechanism generating radio halos is still heavily debated. We critically discuss the proposed models for the radio halo emission and highlight the weaknesses underlying each explanation. We present an idea how the interplay of CR propagation and turbulent advection selects a bimodal spatial CR distribution that is characteristic for the dynamical state of a cluster. As a result, strongly turbulent, merging clusters should have a more centrally concentrated CR energy density profile with respect to relaxed ones with very subsonic turbulence. This translates into a bimodality of the expected diffuse radio and gamma ray emission of clusters. Thus, the observed bimodality of cluster radio halos appears to be a natural consequence of the interplay of CR transport processes, independent of the model of radio halo formation, be it hadronic interactions of CR protons or re-acceleration of low-energy CR electrons.

Primordial magnetic fields generated in the early universe are subject of considerable investigation, and observational limits on their strength are required to constrain the theory. Due to their impact on the reionization process, the strength of primordial fields can be limited using the latest data on reionization and the observed UV-luminosity function of high-redshift galaxies. Given the steep faint-end slope of the luminosity function, faint galaxies contribute substantial ionizing photons, and the low-luminosity cutoff has an impact on the total budget thereof. Magnetic pressure from primordial fields affects such cutoff by preventing collapse in halos with mass below 10^{10} M_solar (B_0 / 3 nG)^3, with B_0 the co-moving field strength. In this letter, the implications of these effects are consistently incorporated in a simplified model for reionization, and the uncertainties due to the cosmological parameters, the reionization parameters and the observed UV luminosity function are addressed. We show that the observed ionization degree at z\sim7 leads to the strongest upper limit of B_0\lsim 2-3nG. Stronger limits could follow from measurements of high ionization degree at z>7.

Massive structures like cluster of galaxies, embedded in cosmic filaments, release enormous amount of energy through their interactions. These events are associated with production of Mpc-scale shocks and injection of considerable amount of turbulence, affecting the non-thermal energy budget of the ICM. In order to study this thoroughly, we performed a set of cosmological simulations using the hydrodynamical code Enzo. We studied the formation of clusters undergoing major mergers, the propagation of merger shocks and their interaction with the filamentary cosmic web. This interaction is shown to produce peripheral structures remarkably similar to giant radio relics observed, for example, in Abell 3376 and Abell 3667. We find a relatively long timescale (about 4 Gyr) for turbulence decay in the centre of major merging clusters. This timescale is substantially longer than typically assumed in the turbulent re-acceleration models, invoked for explaining the statistics of observed radio halos.

We present the implementation of a three-dimensional, second order accurate Godunov-type algorithm for magneto-hydrodynamic (MHD), in the adaptive-mesh-refinement (AMR) cosmological code {\tt CHARM}. The algorithm is based on the full 12-solve spatially unsplit Corner-Transport-Upwind (CTU) scheme. The fluid quantities are cell-centered and are updated using the Piecewise-Parabolic-Method (PPM), while the magnetic field variables are face-centered and are evolved through application of the Stokes theorem on cell edges via a Constrained-Transport (CT) method. The multidimensional MHD source terms required in the predictor step for high-order accuracy are applied in a simplified form which reduces their complexity in three dimensions without loss of accuracy or robustness. The algorithm is implemented on an AMR framework which requires specific synchronization steps across refinement levels. These include face-centered restriction and prolongation operations and a {\it reflux-curl} operation, which maintains a solenoidal magnetic field across refinement boundaries. The code is tested against a large suite of test problems, including convergence tests in smooth flows, shock-tube tests, classical two- and three-dimensional MHD tests, a three-dimensional shock-cloud interaction problem and the formation of a cluster of galaxies in a fully cosmological context. The magnetic field divergence is shown to remain negligible throughout.

High energy electrons and positrons from annihilating dark matter can imprint unique angular anisotropies on the diffuse gamma-ray flux by inverse Compton scattering off the interstellar radiation field. We develop a numerical tool to compute gamma-ray emission from such electrons and positrons produced in the smooth host halo and in substructure halos with masses down to 10^(-6)M_sun. We show that the angular power spectrum from inverse Compton scattering is exponentially suppressed below an angular scale determined by the diffusion length of electrons and positrons. We also find that the total flux and the shape of the angular power spectrum depends sensitively on the spatial distribution of subhalos in the Milky Way. Finally, the contribution from the smooth host halo component to the gamma-ray mean intensity is negligibly small compared to subhalos.

Deep X-ray data from CHANDRA and XMM-Newton, and GMRT radio data are presented for ZwCl 2341.1+0000, an extremely unusual and complex merging cluster of galaxies at the intersection of optical filaments. We propose that energetics of multiple mergers and accretion flows has resulted in wide-spread shocks, acceleration of cosmic ray particles and amplification of weak magnetic fields. This results in Mpc-scale peripheral radio relics and halo like non-thermal emission observed near the merging center.

We examine in detail the recent proposal that extreme Cosmic-Ray-Dominated-Regions (CRDRs) characterize the ISM of galaxies during events of high-density star formation, fundamentally altering its initial conditions (Papadopoulos 2010). Solving the coupled chemical and thermal state equations for dense UV-shielded gas reveals that the large cosmic ray energy densities in such systems (U_{CR} (few)x(10^3-10^4) U_{CR,Gal}) will indeed raise the minimum temperature of this phase (where the initial conditions of star formation are set) from ~10K (as in the Milky Way) to (50-100)K. Moreover in such extreme CRDRs the gas temperature remains fully decoupled from that of the dust, with T_{kin} >> T_{dust}, even at high densities (n(H_2)~10^5--10^6 cm^{-3}), quite unlike CRDRs in the Milky Way where T_k T_{dust} when n(H_2) >= 10^5 cm^{-3}. These dramatically different star formation initial conditions will: a) boost the Jeans mass of UV-shielded gas regions by factors of ~10--100 with respect to those in quiescent or less extreme star forming systems, and b) "erase" the so-called inflection point of the effective equation of state (EOS) of molecular gas. Both these effects occur across the entire density range of typical molecular clouds, and may represent {\it a new paradigm for all high-density star formation in the Universe}, with cosmic rays as the key driving mechanism, operating efficiently even in the high dust extinction environments of extreme starbursts...

We performed a set of cosmological simulations of major mergers in galaxy clusters to study the evolution of merger shocks and the subsequent injection of turbulence in the post-shock region and in the intra-cluster medium (ICM). The computations were done with the grid-based, adaptive mesh refinement hydro code Enzo, using an especially designed refinement criteria for refining turbulent flows in the vicinity of shocks. A substantial amount of turbulence energy is injected in the ICM due to major merger. Our simulations show that the shock launched after a major merger develops an ellipsoidal shape and gets broken by the interaction with the filamentary cosmic web around the merging cluster. The size of the post-shock region along the direction of shock propagation is about 300 kpc h^-1, and the turbulent velocity dispersion in this region is larger than 100 km s^-1. Scaling analysis of the turbulence energy with the cluster mass within our cluster sample is consistent with M^(5/3), i.e. the scaling law for the thermal energy in the self-similar cluster model. This clearly indicates the close relation between virialization and injection of turbulence in the cluster evolution. We found that the ratio of the turbulent to total pressure in the cluster core within 2 Gyr after the major merger is larger than 10%, and it takes about 4 Gyr to get relaxed, which is substantially longer than typically assumed in the turbulent re-acceleration models, invoked to explain the statistics of observed radio halos. Striking similarities in the morphology and other physical parameters between our simulations and the "symmetrical radio relics" found at the periphery of the merging cluster A3376 are finally discussed. In particular, the interaction between the merger shock and the filaments surrounding the cluster could explain the presence of "notch-like" features at the edges of the double relics.

We investigate the interplay of cosmic ray (CR) propagation and advection in galaxy clusters. Propagation in form of CR diffusion and streaming tends to drive the CR radial profiles towards being flat, with equal CR number density everywhere. Advection of CR by the turbulent gas motions tends to produce centrally enhanced profiles. We assume that the CR streaming velocity is of the order of the sound velocity. This is motivated by plasma physical arguments. The CR streaming is then usually larger than typical advection velocities and becomes comparable or lower than this only for periods with trans- and super-sonic cluster turbulence. As a consequence a bimodality of the CR spatial distribution results. Strongly turbulent, merging clusters should have a more centrally concentrated CR energy density profile with respect to relaxed ones with very subsonic turbulence. This translates into a bimodality of the expected diffuse radio and gamma-ray emission of clusters, since more centrally concentrated CR will find higher target densities for hadronic CR proton interactions, higher plasma wave energy densities for CR electron and proton re-acceleration, and stronger magnetic fields. Thus, the observed bimodality of cluster radio halos appears to be a natural consequence of the interplay of CR transport processes, independent of the model of radio halo formation, be it hadronic interactions of CR protons or re-acceleration of low-energy CR electrons. Energy dependence of the CR propagation should lead to spectral steepening of dying radio halos. Furthermore, we show that the interplay of CR diffusion with advection implies first order CR re-acceleration in the pressure-stratified atmospheres of galaxy clusters. Finally, we argue that CR streaming could be important in turbulent cool cores of galaxy clusters since it heats preferentially the central gas with highest cooling rate.

Relativistic charged particles (CR for cosmic-rays) produced by Supernova explosion of the first generation of massive stars that are responsible for the re-ionization of the universe escape into the intergalactic medium, carrying an electric current. Charge imbalance and induction give rise to a return current, $\vec j_t$, carried by the cold thermal plasma which tends to cancel the CR current. The electric field, $\vec E=\eta \vec j_t$, required to draw the collisional return current opposes the outflow of low energy cosmic rays and ohmically heats the cold plasma. Owing to inhomogeneities in the resistivity, $\eta(T)$, caused by structure in the temperature, $T$, of the intergalactic plasma, the electric field possesses a rotational component which sustains Faraday's induction. It is found that magnetic field is robustly generated throughout intergalactic space at rate of 10$^{-17}-10^{-16}$ Gauss/Gyr, until the temperature of the intergalactic medium is raised by cosmic reionization. The magnetic field may seed the subsequent growth of magnetic fields in the intergalactic environment.

We present a stable and convergent method for studying a system of gas and dust, coupled through viscous drag in both non-stiff and stiff regimes. To account for the effects of dust drag in the update of the fluid quantities, we employ a fluid description of the dust component and study the modified gas-dust hyperbolic system following the approach in Miniati & Colella (2007). In addition to two entropy waves for the gas and dust components, respectively, the extended system includes three waves driven partially by gas pressure and partially by dust drift, which, in the limit of vanishing coupling, tend to the two original acoustic waves and the unhindered dust streaming. Based on this analysis we formulate a predictor step providing first order accurate reconstruction of the time-averaged state variables at cell interfaces, whence a second order accurate estimate of the conservative fluxes can be obtained through a suitable linearized Riemann solver. The final source term update is carried out using a one-step, second order accurate, L-stable, predictor corrector asymptotic method (the alpha-QSS method suggested by Mott et. al. 2000). This procedure completely defines a two-fluid method for gas-dust system. Using the updated fluid solution allows us to then advance the individual particle solutions, including self-consistently the time evolution of the gas velocity in the estimate of the drag force. This is done with a suitable particle scheme also based on the alpha-QSS method. A set of benchmark problems shows that our method is stable and convergent. When dust is modeled as a fluid (two-fluid) second order accuracy is achieved in both stiff and non-stiff regimes, whereas when dust is modeled with particles (hybrid) second order is achieved in the non-stiff regime and first order otherwise.

We present a catalogue of MgII absorption systems obtained from high resolution UVES/VLT data of 77 QSOs in the redshift range 0.6 < z < 2.0, and down to an equivalent width W_0 > 0.1 \AA. The statistical properties of our sample are found to be in agreement with those from previous work in the literature. However, we point out that the previously observed increase with redshift of dN/dz for weak absorbers, pertains exclusively to very weak absorbers with W_0 < 0.1 \AA. Instead, dN/dz for absorbers with W_0 in the range 0.1-0.3 \AA actually decreases with redshift, similarly to the case of strong absorbers. We then use this catalogue to extend our earlier analysis of the links between the Faraday Rotation Measure of the quasars and the presence of intervening MgII absorbing systems in their spectra. In contrast to the case with strong MgII absorption systems W_0 > 0.3 \AA, the weaker systems do not contribute significantly to the observed Rotation Measure of the background quasars. This is possibly due to the higher impact parameters of the weak systems compared to strong ones, suggesting that the high column density magnetized material that is responsible for the Faraday Rotation is located within about 50 kpc of the galaxies. Finally, we show that this result also rules out the possibility that some unexpected secondary correlation between the quasar redshift and its intrinsic Rotation Measure is responsible for the association of high Rotation Measure and strong intervening MgII absorption that we have presented elsewhere, since this would have produced an equal effect for the weak absorption line systems, which exhibit a very similar distribution of quasar redshifts.

Context: Hierarchal models of large scale structure (LSS) formation predict that galaxy clusters grow via gravitational infall and mergers of (smaller) mass concentrations, such as clusters and galaxy groups. Diffuse radio emission, in the form of radio halos and relics, is found in clusters undergoing a merger, indicating that shocks or turbulence associated with the merger are capable of accelerating electrons to highly relativistic energies. Here we report on radio observations of ZwCl 2341.1+0000, a complex merging structure of galaxies located at z=0.27, using Giant Metrewave Radio Telescope (GMRT) observations. Aims: The main aim of the observations is to study the nature of the diffuse radio emission in the galaxy cluster ZwCl 2341.1+0000. Methods: We have carried out GMRT 610, 241, and 157 MHz continuum observations of ZwCl 2341.1+0000. The radio observations are combined with X-ray and optical data of the cluster. Results: The GMRT observations show the presence of a double peripheral radio relic in the cluster ZwCl 2341.1+0000. The spectral index is -0.49 \pm 0.18 for the northern relic and -0.76 \pm 0.17 for the southern relic respectively. We have derived values of 0.48-0.93 microGauss for the equipartition magnetic field strength. The relics are probably associated with an outwards traveling merger shock waves.

The origin and growth of magnetic fields in galaxies is still something of an enigma. It is generally assumed that seed fields are amplified over time through the dynamo effect, but there are few constraints on the timescale. It has recently been demonstrated that field strengths as traced by rotation measures of distant quasars are comparable to those seen today, but it was unclear whether the high fields were in the exotic environments of the quasars themselves or distributed along the line of sight. Here we demonstrate that the quasars with strong MgII absorption lines are unambiguously associated with larger rotation measures. Since MgII absorption occurs in the haloes of normal galaxies along the sightline to the quasars, this association requires that organized fields of surprisingly high strength are associated with normal galaxies when the Universe was only about one-third of its present age.

Large-scale accretion shocks around massive clusters of galaxies, generically expected in the cold dark matter scenario of cosmological structure formation, are shown to be plausible sources of the observed ultrahigh energy cosmic rays (UHECRs) by accelerating a mixture of heavy nuclei including the iron group elements. Current observations can be explained if the source composition at injection for the heavier nuclei is somewhat enhanced from simple expectations for the accreting gas. The proposed picture should be clearly testable by current and upcoming facilities in the near future through characteristic features in the UHECR spectrum, composition and anisotropy, in particular the rapid increase of the average mass composition with energy from $10^{19}$ to $10^{20}$ eV.

Faraday rotation (RM) probes of magnetic fields in the universe are sensitive to cosmological and evolutionary effects as $z$ increases beyond $\sim $1 because of the scalings of electron density and magnetic fields, and the growth in the number of expected intersections with galaxy-scale intervenors, $d$N/$dz$. In this new global analysis of an unprecedented large sample of RM's of high latitude quasars extending out to $z\sim $3.7 we find that the distribution of RM broadens with redshift in the 20 $-$ 80 rad m$^{-2}$ range range, despite the (1 +$z$)$^{-2}$ wavelength dilution expected in the observed Faraday rotation. Our results indicate that the Universe becomes increasingly ``Faraday-opaque'' to sources beyond $z \sim$ 2, that is, as $z$ increases progressively fewer sources are found with a ``small'' RM in the observer's frame. This is in contrast to sources at $z \la$1. They suggest that the environments of galaxies were significantly magnetized at high redshifts, with magnetic field strengths that were at least as strong within a few Gyr of the Big Bang as at the current epoch. We separately investigate a simple unevolving toy model in which the RM is produced by MgII absorber systems, and find that it can approximately reproduce the observed trend with redshift. An additional possibility is that the intrinsic RM associated with the radio sources was much higher in the past, and we show that this is not a trivial consequence of the higher radio luminosities of the high redshift sources.

To understand the origin of ultra-high energy cosmic rays (UHECRs, defined to be above 10^18 eV), it is required to model in a realistic way their propagation in the Universe. UHECRs can interact with low energy radio, microwave, infrared and optical photons to produce electron/positron pairs or pions. The latter decay and give rise to neutrinos and electromagnetic cascades extending down to MeV energies. In addition, deflections in cosmic magnetic fields can influence the spectrum and sky distribution of primary cosmic rays and, due to the increased propagation path length, the secondary neutrino and gamma-ray fluxes. Neutrino, gamma-ray, cosmic ray physics and extra-galactic magnetic fields are, therefore, strongly linked subjects and should be considered together in order to extract maximal information from existing and future data, like the one expected from the Auger Observatory. For that purpose, we have developed CRPropa, a publicly-available numerical package which takes into account interactions and deflections of primary UHECRs as well as propagation of secondary electromagnetic cascades and neutrinos. CRPropa allows to compute the observable properties of UHECRs and their secondaries in a variety of models for the sources and propagation of these particles. Here we present physical processes taken into account as well as benchmark examples; a detailed documentation of the code can be found on our web site.

Large-scale accretion shocks around massive clusters of galaxies, generically expected in hierarchical scenarios of cosmological structure formation, are shown to be potential sources of the observed ultrahigh energy cosmic rays (UHECRs) by accelerating a mixture of heavy nuclei including the iron group elements. Current observations can be explained if the source composition at injection for the heavier nuclei is somewhat enhanced from simple expectations for the accreting gas. The proposed picture should be testable by current and upcoming facilities in the near future through characteristic features in the UHECR spectrum, composition and anisotropy. The associated X-ray and gamma-ray signatures are also briefly discussed.

A numerical method for integrating the equations describing a dynamically coupled system made of a fluid and cosmic-rays is developed. In smooth flows the effect of CR pressure is accounted for by modification of the characteristic equations and the energy exchange between cosmic-rays and the fluid, due to diffusive processes in configuration and momentum space, is modeled with a flux conserving method. Provided the shock acceleration efficiency as a function of the upstream conditions and shock Mach number, we show that the Riemann solver can be modified to take into account the cosmic-ray mediation without having to resolve the cosmic-ray induced substructure. Shocks are advanced with Glimm's method which preserves their discontinuous character without any smearing, thus allowing to maintain self-consistency in the shock solutions. In smooth flows either Glimm's or a higher order Godunov's method can be applied, with the latter producing better results when approximations are introduced in the Riemann solver.

Notwithstanding the advent of the Gamma-ray Large Area Telescope, theoretical models predict that a significant fraction of the cosmic gamma-ray background (CGB), at the level of 20% of the currently measured value, will remain unresolved. The angular power spectrum of intensity fluctuations of the CGB contains information on its origin. We show that probing the latter from a few tens of arcmin to several degree scales, together with complementary GLAST observations of gamma-ray emission from galaxy clusters and the blazars luminosity function, can discriminate between a background that originates from unresolved blazars or cosmic rays accelerated at structure formation shocks.

We present a new numerical algorithm for the solution of coupled collisional and collisionless systems, based on the block structured adaptive mesh and time refinement strategy (AMR). We describe the issues associated with the discretization of the system equations and the synchronization of the numerical solution on the hierarchy of grid levels. We implement a code based on a higher order, conservative and directionally unsplit Godunov's method for hydrodynamics; a symmetric, time centered modified symplectic scheme for collisionless component; and a multilevel, multigrid relaxation algorithm for the elliptic equation coupling the two components. Numerical results that illustrate the accuracy of the code and the relative merit of various implemented schemes are also presented.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.

Hyperbolic conservation laws with stiff source terms appear in the study of a variety of physical systems. Early work showed that the use of formally second-order accurate semi-implicit methods could lead to a substantial loss of accuracy, due to inconsistencies between the flux calculation without sources and the limiting equilibrium behavior of the gas. In this paper we present an efficient second order accurate scheme to treat stiff source terms within the framework of higher order Godunov's methods. We employ Duhamel's formula to devise a modified predictor step which accounts for the effects of stiff source terms on the conservative fluxes and recovers the correct isothermal behavior in the limit of an infinite cooling/reaction rate. Source term effects on the conservative quantities are fully accounted for by means of a one-step, second order accurate semi-implicit corrector scheme based on the deferred correction method of Dutt et. al. We demostrate the accurate, stable and convergent results of the proposed method through a set of benchmark problems for a variety of stiffness conditions and source types.

Nearby sources of cosmic rays up to a ZeV(=10^21 eV) could be observed with a multi-messenger approach including secondary gamma-rays and neutrinos. If cosmic rays above ~10^18 eV are produced in magnetized environments such as galaxy clusters, the flux of secondary gamma-rays below ~1 TeV can be enhanced up to several orders of magnitudes compared to unmagnetized sources. A particular source of enhancement are synchrotron and cascade photons from e^+e^- pairs produced by protons from sources with relatively steep injection spectra proportional to E^-2.6. Such sources should be visible at the same time in ultra-high energy cosmic ray experiments and gamma-ray telescopes.

A 50 ksec Chandra observation of the unidentified TeV source in Cygnus reported by the HEGRA collaboration reveals no obvious diffuse X-ray counterpart. However, 240 Pointlike X-ray sources are detected within or nearby the extended TeV J2032+4130 source region, of which at least 36 are massive stars and 2 may be radio emitters. That the HEGRA source is a composite, having as counterpart the multiple point-like X-ray sources we observe, cannot be ruled out. Indeed, the distribution of point-like X-ray sources appears non-uniform and concentrated broadly within the extent of the TeV source region. We offer a hypothesis for the origin of the very high energy gamma-ray emission in Cyg OB2 based on the local acceleration of TeV range cosmic rays and the differential distribution of OB vs. less massive stars in this association.

Using the mosaic of nine XMM-Newton observations, we study the hydrodynamic state of the merging cluster of galaxies Abell 3266. The high quality of the spectroscopic data and large field of view of XMM-Netwon allow us to determine the thermodynamic conditions of the intracluster medium on scales of order of 50 kpc. A high quality entropy map reveals the presence of an extended region of low entropy gas, running from the primary cluster core toward the northeast along the nominal merger axis. The mass of the low entropy gas amounts to approximately 2e13 solar masses, which is comparable to the baryonic mass of the core of a rich cluster. We test the possibility that the origin of the observed low entropy gas is either related to the disruption a preexisting cooling core in Abell 3266 or to the stripping of gas from an infalling subcluster companion. We find that both the radial pressure and entropy profiles as well as the iron abundance of Abell 3266 do not resemble those in other known cooling core clusters (Abell 478). Thus we conclude that the low entropy region is subcluster gas in the process of being stripped off from its dark matter halo. In this scenario the subcluster would be falling onto the core of A3266 from the foreground. This would also help interpret the observed high velocity dispersion of the galaxies in the cluster center, provided that the mass of the subcluster is at most a tenth of the mass of the main cluster.

We combine a high-resolution hydro-simulation of the LambdaCDM cosmology with two radiative transfer schemes (for continuum and line radiation) to predict the properties, spectra and spatial distribution of fluorescent Ly-alpha emission at z~3. We focus on line radiation produced by recombinations in the dense intergalactic medium ionized by UV photons. In particular, we consider both a uniform background and the case where gas clouds are illuminated by a nearby quasar. We find that the emission from optically thick regions is substantially less than predicted from the widely used static, plane-parallel model. The effects induced by a realistic velocity field and by the complex geometric structure of the emitting regions are discussed in detail. We make predictions for the expected brightness and size distributions of the fluorescent sources.Our results account for recent null detections and can be used to plan new observational campaigns both in the field (to measure the intensity of the diffuse UV background) and in the proximity of bright quasars (to understand the origin of high colum-density absorbers).

We compare the propagation of iron and proton nuclei above 10^19 eV in a structured Universe with source and magnetic field distributions obtained from a large scale structure simulation and source densities about 10^(-5) Mpc^(-3). All relevant cosmic ray interactions are taken into account, including photo-disintegration and propagation of secondary products. Iron injection predicts spectral shapes different from proton injection which disagree with existing data below about 30 EeV. Injection of light nuclei or protons must therefore contribute at these energies. However, at higher energies, existing data are consistent with injection of pure iron with spectral indices between 2 and 2.4. This allows a significant recovery of the spectrum above roughly 100 EeV, especially in the case of large deflections. Significant auto-correlation and anisotropy, and considerable cosmic variance are also predicted in this energy range. The mean atomic mass A fluctuates considerably between different scenarios. At energies below 60 EeV, if the observed A > 35, magnetic fields must have a negligible effect on propagation. At the highest energies the observed flux will be dominated by only a few sources whose location may be determined by next generation experiments to within 10-20 degrees even if extra-galactic magnetic fields are important.

We discuss the influence of large scale cosmic magnetic fields on the propagation of hadronic cosmic rays above 10^19 eV based on large scale structure simulations. Our simulations suggest that rather substantial deflection up to several tens of degrees at 10^20 eV are possible for nucleon primaries. Further, spectra and composition of cosmic rays from individual sources can depend on magnetic fields surrounding these sources in intrinsically unpredictable ways. This is true even if deflection from such individual sources is small. We conclude that the influence of large scale cosmic magnetic fields on ultra-high energy cosmic ray propagation is currently hard to quantify. We discuss possible reasons for discrepant results of simulations by Dolag et al. which predict deflections of at most a few degrees for nucleons. We finally point out that even in these latter simulations a possible heavy component would in general suffer substantial deflection.

We study signatures of a structured universe in the multi-pole moments, auto-correlation function, and cluster statistics of ultra-high energy cosmic rays above 10^19 eV. We compare scenarios where the sources are distributed homogeneously or according to the baryon density distribution obtained from a cosmological large scale structure simulation. The influence of extragalactic magnetic fields is studied by comparing the case of negligible fields with fields expected to be produced along large scale shocks with a maximal strength consistent with observations. We confirm that strongly magnetized observers would predict considerable anisotropy on large scales, which is already in conflict with current data. In the best fit scenario only the sources are strongly magnetized, although deflection can still be considerable, of order 20 degrees up to 10^20 eV, and a pronounced GZK cutoff is predicted. We then discuss signatures for future large scale full-sky detectors such as the Pierre Auger and EUSO projects. Auto-correlations are sensitive to the source density only if magnetic fields do not significantly affect propagation. In contrast, for a weakly magnetized observer, degree scale auto-correlations below a certain level indicate magnetized discrete sources. It may be difficult even for next generation experiments to distinguish between structured and unstructured source distributions.

Spatially-resolved gas pressure maps of the Coma galaxy cluster are obtained from a mosaic of XMM-Newton observations in the scale range between a resolution of 20 kpc and an extent of 2.8 Mpc. A Fourier analysis of the data reveals the presence of a scale-invariant pressure fluctuation spectrum in the range between 40 and 90 kpc and is found to be well described by a projected Kolmogorov/Oboukhov-type turbulence spectrum. Deprojection and integration of the spectrum yields the lower limit of $\sim 10$ percent of the total intracluster medium pressure in turbulent form. The results also provide observational constraints on the viscosity of the gas.

We present an observational study of the Lx vs Lb sigma^2 relation for early-type galaxies in the Coma cluster based on the XMM-Newton survey data. Compared to a similar relation for a sample dominated by field early-type galaxies, the Coma cluster galaxies show a flatter slope. Our calculations show that adiabatic compression produces a flattening in the Lx vs Lb sigma^2 relation that is in remarkable agreement with the observed effect. Our scenario is further supported by the observed compactness of the X-ray emission of Coma galaxies.

Fluctuations in the redshifted 21 centimeter emission from neutral hydrogen probe the epoch of reionization. We examine the observability of this signal and the impact of extragalactic foreground radio sources. We use cosmological simulations to predict the angular correlation functions of intensity fluctuations due to unresolved radio galaxies, cluster radio halos and relics and free-free emission from the interstellar and intergalactic medium at the frequencies and angular scales relevant for the proposed 21cm tomography. In accord with previous findings, the brightness temperature fluctuations due to foreground sources are much larger than those from the primary 21cm signal at all scales. In particular, diffuse cluster radio emission, which has been previously neglected, provides the most significant foreground contamination. However, we show that the contribution to the angular fluctuations at scales \theta > 1' is dominated by the spatial clustering of bright foreground sources. This excess can be removed if sources above flux levels S > 0.1 mJy are detected and removed. Hence, efficient source removal may be sufficient to allow the detection of angular fluctuations in the 21cm emission free of extragalactic foregrounds at \theta > 1 arcmin. In addition, the removal of sources above S=0.1 mJy also reduces the foreground fluctuations to roughly the same level as the 21cm signal at scales \theta < 1 arcmin. This should allow the substraction of the foreground components in frequency space, making it possible to observe in detail the topology and history of reionization.

The contribution to the ultraviolet background (UVB) from thermal emission due to gas shock heated during cosmic structure formation is assessed with an updated version of Press-Schechter (Sheth & Tormen 1999) formalism. The calculation is consistent with empirical estimates based on the observed properties of galaxies and the observed cosmic star formation history. The bulk of the radiation turns out to be produced by objects in the mass range 10^{11- 13} M_solar, i.e. large galaxies and small groups. When compared to more conventional components (QSOs, stellar) it is found that near 1 Ry, thermal emission is comparable to stellar contributions and amounts to about 10 %, 20 % and 35 % of the total flux at redshifts of 3, 4.5 and higher, respectively; more importantly, near the ionization threshold for HeII, the thermal contribution is comparable to the QSO intensity already at redshift ~3 and dominates at redshifts above 4. In addition, thermal photons alone are enough to produce and sustain HeII reionization already at z~6. Finally, the observed Gunn-Peterson effect at high redshifts (3<z<6) constrains the escape fraction of ionizing photons from galaxies to less than a few percent.

I present a review of nonthermal processes in the large scale structure of the universe. After examining the properties of cosmological shock waves and their role as particle accelerators, I discuss the main observational facts, from radio to \gr and describe the processes that are thought be responsible for the observed nonthermal emissions. Finally, I emphasize the important role of $\gamma$-ray astronomy for the progress in the field. Non detections (upper limits) at these photon energies have already allowed us to draw important conclusions. Future observations will tell us more about the physical conditions in the intracluster medium, physics of shocks dissipation, aspects of CR acceleration.

We calculate the contribution to the ultraviolet background (UVB) from thermal emission from gas shock heated by cosmic structure formation. Our main calculation is based on an updated version of Press-Schechter theory. It is consistent with a more empirical estimate based on the observed properties of galaxies and the observed cosmic star formation history. Thermal UVB emission is characterized by a hard spectrum extending well beyond 4 Ry. The bulk of the radiation is produced by objects in the mass range 10^11-13 M_solar, i.e. large galaxies and small groups. We compute a composite UVB spectrum due to QSO, stellar and thermal components. The ratio of the UVB intensities at the H and He Lyman limits increases from 60 at z=2 to more than 300 at z=6. A comparison of the resulting photoionization rates to the observed Gunn-Peterson effect at high redshifts constrains the escape fraction of ionizing photons from galaxies to be less than a few percent. Near 1 Ry, thermal and stellar emission are comparable amounting to about 10 %, 20 % and 35 % of the total flux at redshifts of 3, 4.5 and higher, respectively. However, near the ionization threshold for HeII, the thermal contribution is much stronger. It is comparable to the QSO intensity already at redshift ~3 and dominates at redshifts above 4. Thermal photons alone are enough to produce and sustain HeII reionization already at z~6. We discuss the possible implications of our results for the thermal history of the intergalactic medium, in particular for HeII reionization.

The environmental influence of cluster media on its member galaxies, known as Butcher--Oemler effect, has recently been subject to revision due to numerous observations of strong morphological transformations occurring outside the cluster virial radii, caused by some unidentified gas removal processes. In this context we present new XMM-Newton observations of M86 group. The unique combination of high spatial and spectral resolution and large field of view of XMM-Newton allows an in-depth investigation of the processes involved in the spectacular disruption of this object. We identify a possible shock with Mach number of ~1.4 in the process of crushing the galaxy in the North-East direction. The latter is ascribed to the presence of a dense X-ray emitting filament, previously revealed in the RASS data. The shock is not associated with other previously identified features of M86 X-ray emission, such as the plume, the north-eastern arm and the southern extension, which are found to have low entropy, similar to the inner 2 kpc of M86. Finally, mere existence of the large scale gas halo around the M86 group, suggests that the disruptions of M86's X-ray halo may be caused by small-scale types of interactions such as galaxy-galaxy collisions.

We investigate the impact of a structured universe in the multi-pole moments, auto-correlation function, and cluster statistics of cosmic rays above 10^19 eV. We compare structured and uniform source distributions with and without magnetic fields obtained from a cosmological simulation. We find that current data marginally favor structured source distributions and magnetic fields reaching a few micro Gauss in galaxy clusters but below 0.1 micro Gauss in our local extragalactic neighborhood. A pronounced GZK cutoff is also predicted in this scenario. Future experiments will make the degree scale auto-correlation function a sensitive probe of micro Gauss fields surrounding the sources.

In this paper we compute the luminosity function of radio relics. In our calculation we include only those relics that are produced by the compression of former radio cocoons. This compression is provided by shocks that are generated in the process of structure formation. Starting from an analytical model for the luminosity evolution of ageing radio cocoons, the luminosity function of radio galaxies and the statistics of shocks as inferred from cosmological simulations, we are able to make the first estimates of the brightness distribution of radio relics. The computed luminosity function is consistent with current observations and predicts that more than $10^3$ radio relics should be discovered with the upcoming generation of low-frequency radio telescopes. Moreover, we predict that radio relics are predominantly found in low-pressure regions outside the cores of clusters.

The HEGRA Cherenkov telescope array group recently reported a steady and extended unidentified TeV gamma-ray source lying at the outskirts of Cygnus OB2. This is the most massive stellar association known in the Galaxy, estimated to contain ~2600 OB type members alone. It has been previously argued that the large scale shocks and turbulence induced by the multiple interacting supersonic winds from the many young stars in such associations may play a role in accelerating Galactic cosmic rays. Indeed, Cyg OB2 also coincides with the non-variable MeV-GeV range unidentified EGRET source, 3EG 2033+4118. We report on the near-simultaneous follow-up observations of the extended TeV source region with the CHANDRA X-ray Observatory and the Very Large Array (VLA) radio telescope obtained in order to explore this possibility. Analysis of the CO, HI, and IRAS 100 micron emissions shows that the TeV source region coincides with an outlying sub-group of powerful OB stars which have evacuated or destroyed much of the ambient atomic, molecular and dust material, and which may be related to the very high-energy emissions. An interesting SNR-like structure is also revealed near the TeV source region in the CO, HI and radio emission maps. Applying a numerical simulation which accurately tracks the radio to gamma-ray emission from primary hadrons as well as primary and secondary e+/-, we find that the broadband spectrum of the TeV source region favors a predominantly nucleonic - rather than electronic - origin of the high-energy flux, though deeper X-ray and radio observations are needed to confirm this. A very reasonable, ~0.1%, conversion efficiency of Cyg OB2's extreme stellar wind mechanical luminosity to nucleonic acceleration to ~PeV (10^15 eV) energies is sufficient to explain the multifrequency emissions.

We simulate propagation of cosmic ray nucleons above 10^{19} eV in scenarios where both the source distribution and magnetic fields within about 50 Mpc from us are obtained from an unconstrained large scale structure simulation. We find that consistency of predicted sky distributions with current data above 4 x 10^{19} eV requires magnetic fields of ~0.1 microGauss in our immediate environment, and a nearby source density of ~10^{-4}-10^{-3} Mpc^{-3}. Radio galaxies could provide the required sources, but only if both high and low-luminosity radio galaxies are very efficient cosmic ray accelerators. Moreover, at ~10^{19} eV an additional isotropic flux component, presumably of cosmological origin, should dominate over the local flux component by about a factor three in order to explain the observed isotropy. This argues against the scenario in which local astrophysical sources of cosmic rays above ~10^{19} eV reside in strongly magnetized (B~0.1 microGauss) and structured intergalactic medium. Finally we discuss how future large scale full-sky detectors such as the Pierre Auger project will allow to put much more stringent constraints on source and magnetic field distributions.

The origins of all extra-solar cosmic rays -- both the ultra high-energy and lower energy Galactic variety -- remain unclear. It has been argued that the large scale shocks and turbulence induced by the multiple interacting supersonic winds from the many hot stars in young stellar associations may play a role in accelerating Galactic cosmic rays -- with or without the associated multiple supernova remnants (eg. Cesarsky & Montmerle, 1983; Bykov 2001). In this context, the report by the HEGRA Cherenkov telescope group of a steady and extended unidentified TeV gamma-ray source lying at the outskirts of Cygnus OB2 is particularly significant. This is the most massive stellar association known in the Galaxy, estimated to contain ~2600 OB type members alone indeed, Cyg OB2 also coincides with the non-variable MeV-GeV range unidentified EGRET source, 3EG 2033+4118. We summarise here the near-simultaneous follow-up observations of the extended TeV source region with the CHANDRA X-ray Observatory and the Very Large Array (VLA) radio telescope. Applying a numerical simulation which accurately tracks the radio to gamma-ray emission from primary hadrons as well as primary and secondary e+/-, we find that the broadband spectrum of the TeV source region favors a predominantly nucleonic rather than electronic origin of the high-energy flux, though deeper X-ray and radio observations will help confirm this. A very reasonable, ~0.1%, conversion efficiency of Cyg OB2's extreme stellar wind mechanical luminosity to nucleonic acceleration to ~PeV (10^15 eV) energies is sufficient to explain the multifrequency emissions.

We investigate the spatial and spectral properties of non-thermal emission from clusters of galaxies at gamma-ray energies between 10 keV and 10 TeV due to inverse-Compton (IC) emission, pion-decay and non-thermal bremsstrahlung (NTB) from cosmic-ray(CR) ions and electrons accelerated at cosmic shock and secondary e+- from inelastic p-p collisions. We identify two main emission region, namely the core (also bright in thermal X-ray) and the outskirts region where accretion shocks occur. IC emission from shock accelerated CR electrons dominate the emission in the outer regions of galaxy clusters, provided that at least a fraction of a percent of the shock ram pressure is converted into CR electrons. A clear detection of this component and of its spatial distribution will allow us direct probing of cosmic accretion shocks. In the cluster core, gamma-ray emission above 100 MeV is dominated by pion-decay mechanism and, at lower energies, by IC emission from secondary e+-. However, IC emission from shock accelerated electrons projected onto the cluster core will not be negligible. We emphasize the importance of separating these emission components for a correct interpretation of the experimental data and outline a strategy for that purpose. Failure in addressing this issue will produce unsound estimates of the intra-cluster magnetic field strength and CR ion content. According to our estimate future space borne and ground based gamma-ray facilities should be able to measure the whole nonthermal spectrum both in the cluster core and at its outskirts. The importance of such measurements in advancing our understanding of non-thermal processes in the intra-cluster medium is discussed.

We address the role of gamma-ray astronomy in the investigation of nonthermal processes in the large scale structure of the universe. Based on EGRET upper limits on nearby galaxy clusters (GCs) we constrain the acceleration efficiency of CR electrons at intergalactic shocks to <=1 % than the shock ram pressure. That implies a contribution to the cosmic gamma-ray background from intergalactic shocks of order 25 % of the measured level. We model spatial and spectral properties of nonthermal gamma-ray emission due to shock accelerated cosmic-rays (CRs) in GCs and emphasize the importance of imaging capability of upcoming gamma-ray facilities for a correct interpretation of the observational results. GC observations at this photon energy will help us understand the origin of the radio emitting particles, the possible level of CR pressure and the strength of magnetic fields in intracluster environment and possibly will allow us detect the accretion shocks.

We investigate numerically the contribution to the cosmic gamma-ray background from cosmic-rays ions and electrons accelerated at intergalactic shocks associated with cosmological structure formation. We show that the kinetic energy of accretion flows in the low-red-shift inter-galactic medium is thermalized primarily through moderately strong shocks, which allow for an efficient conversion of shock ram pressure into cosmic-ray pressure. Cosmic-rays accelerated at these shocks produce a diffuse gamma-ray flux which is dominated by inverse Compton emission from electrons scattering off cosmic microwave background photons. Decay of neutral pi-mesons generated in p-p inelastic collisions of the ionic cosmic-ray component with the thermal gas contribute about 30% of the computed emission. Based on experimental upper limits on the photon flux above 100 MeV from nearby clusters we constrain the efficiency of conversion of shock ram pressure into relativistic CR electrons to < 1%. Thus, we find that cosmic-rays of cosmological origin can generate an overall significant fraction of order 20% and no more than 30% of the measured gamma-ray background.

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.

Gamma-ray astronomy will play a crucial role in the investigation of nonthermal processes in the large scale structure of the universe. Particularly, galaxy clusters (GC) observations at this photon energy will help us understand the origin of radio emitting high energy particles, the possible level of cosmic-ray (CR) pressure in intracluster environment, and the strength of intracluster magnetic fields. In addition here we point out the importance of these observations for a possible detection of cluster accretion shocks and for constraining their CR acceleration efficiency. We model spatial and spectral properties of gamma-ray emission due to shock accelerated CRs in GC and emphasize the importance of imaging capability of upcoming gamma-ray facilities for a correct interpretation of observational results.

We report the discovery of large-scale diffuse radio emission from what appears to be a large-scale filamentary network of galaxies in the region of cluster ZwCl 2341.1+0000, and stretching over an area of at least $6 h^{-1}_{50}$ Mpc in diameter. Multicolour CCD observations yield photometric redshifts indicating that a significant fraction of the optical galaxies in this region is at a redshift of z=0.3. This is supported by spectroscopic measurements of 4 galaxies in the SDSS survey at a mean z=0.27. We present VLA images at 20 cm (NVSS) and 90 cm wavelengths, showing the detailed radio structure of the filaments. Comparison with the VLA high resolution FIRST radio survey shows that the diffuse emission is not due to known individual point sources. The diffuse radio-emission has a spectral index $\alpha \lesssim -0.5$, and is most likely synchrotron emission from relativistic charged particles in an inter-galactic magnetic field. Furthermore, this optical/radio structure is detected in X-rays by the ROSAT all-sky survey. It has a 0.1--2.4 keV luminosity of about $10^{44}$ erg s$^{-1}$ and shows an extended highly non-relaxed morphology. These observations suggest that ZwCl 2341.1+0000 is possibly a proto-cluster of galaxies in which we are witnessing the process of structure formation. We show (both analytically and by numerical simulations) that the energetics of accretion shocks generated in forming large-scale structures are sufficient to produce enough high energy cosmic-ray (CR) electrons required to explain the observed radio emission, provided a magnetic field of strength of about 1 micro Gauss is present there.Thus it is the first evidence of cosmic-ray particle acceleration and magnetic fields occuring on a super-cluster scale. (Abridged)

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.

We present COSMOCR, a numerical code for the investigation of cosmic ray related studies in computational cosmology. The code follows the diffusive shock acceleration, the mechanical and radiative energy losses and the spatial transport of the supra-thermal particles in cosmic environment. Primary cosmic ray electrons and ions are injected at shocks according to the thermal leakage prescription. Secondary electrons are continuously injected as a results of p-p inelastic collisions of primary cosmic ray ions and thermal background nuclei. The code consists of a conservative, finite volume method with a power-law sub-grid model in momentum space. Two slightly different schemes are implemented depending on the stiffness of the cooling terms. Comparisons of numerical results with analytical solution for a number of tests of direct interest show remarkable performance of the present code.

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 present results from three-dimensional (3-D) numerical simulations investigating the magnetohydrodynamics of cloud-wind interactions. The initial cloud is spherical while the magnetic field is uniform and transverse to the cloud motion. A simplified analytical model that describes the magnetic energy evolution in front of the cloud is developed and compared with simulation results. In addition, it is found the interaction of the cloud with a magnetized interstellar medium (ISM) results in the formation of a highly structured magnetotail. The magnetic flux in the wake of the cloud organizes into flux ropes and a reconnection, current sheet is developed, as field lines of opposite polarity are brought close together near the symmetry axis. At the same time, magnetic pressure is strongly enhanced at the leading edge of the cloud from the stretching of the field lines that occurs there. This has an important dynamical effect on the subsequent evolution of the cloud, since some unstable modes tend to be strongly enhanced.

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.

We present results from the first three-dimensional numerical simulations of moderately supersonic cloud motion through a tenuous, magnetized medium. We show that the interaction of the cloud with a magnetic field perpendicular to its motion has a great dynamical impact on the development of instabilities at the cloud surface. Even for initially spherical clouds, magnetic field lines become trapped in surface deformations and undergo stretching. The consequent field amplification that occurs there and particularly its variation across the cloud face then dramatically enhance the growth rate of Rayleigh-Taylor unstable modes, hastening the cloud disruption.

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.

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

We extend previous studies of the physics of interstellar cloud collisions by beginning investigation of the role of magnetic fields through 2D magnetohydrodynamic (MHD) numerical simulations. We study head-on collisions between equal mass, mildly supersonic diffuse clouds. We include a moderate magnetic field and two limiting field geometries, with the field lines parallel (aligned) and perpendicular (transverse) to the colliding cloud motion. We explore both adiabatic and radiative cases, as well as symmetric and asymmetric ones. We also compute collisions between clouds evolved through prior motion in the intercloud medium and compare with unevolved cases. We find that: In the (i) aligned case, adiabatic collisions, like their HD counterparts, are very disruptive, independent of the cloud symmetry. However, when radiative processes are taken into account, partial coalescence takes place even in the asymmetric case, unlike the HD calculations. In the (ii) transverse case, collisions between initially adjacent unevolved clouds are almost unaffected by magnetic fields. However, the interaction with the magnetized intercloud gas during the pre-collision evolution produces a region of very high magnetic energy in front of the cloud. In collisions between evolved clouds with transverse field geometry, this region acts like a ``bumper'', preventing direct contact between the clouds, and eventually reverses their motion. The ``elasticity'', defined as the ratio of the final to the initial kinetic energy of each cloud, is about 0.5-0.6 in the cases we considered. This behavior is found both in adiabatic and radiative cases.

A description is given for preserving ${\bmsy\nabla}\cdot{\vec B}=0$ in a magnetohydrodynamic (MHD) code that employs the upwind, Total Variation Diminishing (TVD) scheme and the Strang-type operator splitting for multi-dimensionality. The method is based on the staggered mesh technique to constrain the transport of magnetic field: the magnetic field components are defined at grid interfaces with their advective fluxes on grid edges, while other quantities are defined at grid centers. The magnetic field at grid centers for the upwind step is calculated by interpolating the values from grid interfaces. The advective fluxes on grid edges for the magnetic field evolution are calculated from the upwind fluxes at grid interfaces. Then, the magnetic field can be maintained with ${\bmsy\nabla}\cdot{\vec B}=0$ exactly, if this is so initially, while the upwind scheme is used for the update of fluid quantities. The correctness of the code is demonstrated through tests comparing numerical solutions either with analytic solutions or with numerical solutions from the code using an explicit divergence-cleaning method. Also the robustness is shown through tests involving realistic astrophysical problems.

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

We present a study of the kinetic energy dissipation in interstellar cloud collisions. The main aim is to understand the dependence of the elasticity (defined as the ratio of the final to the initial kinetic energy of the clouds) on the velocity and mass ratio of the colliding clouds, magnetic field strength, and gas metallicity for head-on collisions. The problem has been studied both analytically and via numerical simulations. We have derived handy analytical relationships that well approximate the analogous numerical results. The main findings of this work are: (i) the kinetic energy dissipation in cloud collisions is minimum (i.e. the collision elasticity is maximum) for a cloud relative velocity $v_r \simeq 30 km s^{-1}$; (ii) the above minimum value is proportional $Z L_c^2$, where $Z$ is the metallicity and $L_c$ is the cloud size: the larger is $Z L_c^2$ the more dissipative (i.e. inelastic) the collision will be; (iii) in general, we find that the energy dissipation decreases when the magnetic field strength, and mass ratio of the clouds are increased and the metallicity is decreased, respectively. We briefly discuss the relevance of this study to the global structure of the interstellar medium and to galaxy formation and evolution.

We consider the stability of clouds surrounded by a hotter confining medium with respect to which they are in motion, against Kelvin-Helmholtz instabilities (KHI). In the presence of cooling, sound waves are damped by dissipation. Whenever cooling times are shorter than sound crossing times, as they are in the normal interstellar medium, this implies that the instability generated at the interface of the two media cannot propagate far from the interface itself. To study how this influences the overall stability, first we derive an analytic dispersion relation for cooling media, separated by a shear layer. The inclusion of dissipation does not heal the instability, but it is shown that only a small volume around the interface is affected, the perturbation decaying exponentially with distance from the surface; this is confirmed by numerical simulations. Numerical simulations of spherical clouds moving in a surrounding intercloud medium by which they are pressure confined show that these clouds develop a core/halo structure, with a turbulent halo, and a core in laminar flow nearly unscathed by the KHI. The related and previously reported ``champagne effect'', whereby clouds seem to explode from their top sides, is cured by the inclusion of radiative losses.