Little is known about the origin and basic properties of magnetic fields in clusters of galaxies. High conductivity in magnetized interstellar plasma suggests that galactic magnetic fields are (at least partly) ejected into intracluster (IC) space by the same processes that enrich IC gas with metals. We explore the dispersal of galactic fields by hydrodynamical simulations with our new {\em Enzo-Galcon} code, which is capable of tracking a large number galaxies during cluster assembly, and modeling the processes that disperse their interstellar media. Doing so we are able to describe the evolution of the mean strength of the field and its profile across the cluster. With the known density profile of dispersed gas and an estimated range of coherence scales, we predict the spatial distribution of Faraday rotation measure and find it to be consistent with observational data.

We present our recently developed {\em galcon} approach to hydrodynamical cosmological simulations of galaxy clusters - a subgrid model added to the {\em Enzo} adaptive mesh refinement code - which is capable of tracking galaxies within the cluster potential and following the feedback of their main baryonic processes. Galcons are physically extended galactic constructs within which baryonic processes are modeled analytically. By identifying galaxy halos and initializing galcons at high redshift ($z \sim 3$, well before most clusters virialize), we are able to follow the evolution of star formation, galactic winds, and ram-pressure stripping of interstellar media, along with their associated mass, metals and energy feedback into intracluster (IC) gas, which are deposited through a well-resolved spherical interface layer. Our approach is fully described and all results from initial simulations with the enhanced {\em Enzo-Galcon} code are presented. With a galactic star formation rate derived from the observed cosmic star formation density, our galcon simulation better reproduces the observed properties of IC gas, including the density, temperature, metallicity, and entropy profiles. By following the impact of a large number of galaxies on IC gas we explicitly demonstrate the advantages of this approach in producing a lower stellar fraction, a larger gas core radius, an isothermal temperature profile in the central cluster region, and a flatter metallicity gradient than in a standard simulation.

Measurement sensitivity in the energetic gamma-ray region has improved considerably, and is about to increase further in the near future, motivating a detailed calculation of high-energy (>100 MeV) and very-high-energy (VHE: >100 GeV) gamma-ray emission from the nearby starburst galaxy NGC253. Adopting the convection-diffusion model for energetic electron and proton propagation, and accounting for all the relevant hadronic and leptonic processes, we determine the steady-state energy distributions of these particles by a detailed numerical treatment. The electron distribution is directly normalized by the measured synchrotron radio emission from the central starburst region; a commonly expected theoretical relation is then used to normalize the proton spectrum in this region. Doing so fully specifies the electron spectrum throughout the galactic disk, and with an assumed spatial profile of the magnetic field, the predicted radio emission from the full disk matches well the observed spectrum, confirming the validity of our treatment. The resulting radiative yields of both particles are calculated; the integrated HE and VHE fluxes from the entire disk are predicted to be f(>100 MeV)~2x10^-8 cm^-2 s^-1 and f(>100 GeV)~4x10^-12 cm^-2 s^-1, respectively. We discuss the feasibility of measuring emission at these levels with the space-borne Fermi and the ground-based Cherenkov telescopes.

We describe a subgrid model for including galaxies into hydrodynamical cosmological simulations of galaxy cluster evolution. Each galaxy construct- or galcon- is modeled as a physically extended object within which star formation, galactic winds, and ram pressure stripping of gas are modeled analytically. Galcons are initialized at high redshift (z~3) after galaxy dark matter halos have formed but before the cluster has virialized. Each galcon moves self-consistently within the evolving cluster potential and injects mass, metals, and energy into intracluster (IC) gas through a well-resolved spherical interface layer. We have implemented galcons into the Enzo adaptive mesh refinement code and carried out a simulation of cluster formation in a LambdaCDM universe. With our approach, we are able to economically follow the impact of a large number of galaxies on IC gas. We compare the results of the galcon simulation with a second, more standard simulation where star formation and feedback are treated using a popular heuristic prescription. One advantage of the galcon approach is explicit control over the star formation history of cluster galaxies. Using a galactic SFR derived from the cosmic star formation density, we find the galcon simulation produces a lower stellar fraction, a larger gas core radius, a more isothermal temperature profile, and a flatter metallicity gradient than the standard simulation, in better agreement with observations.

Spurred by the improved measurement sensitivity in the very-high-energy (VHE: >100 GeV) gamma-ray band, we assess the feasibility of detection of the nearby starburst galaxy M82. VHE emission is expected to be predominantly from the decay of neutral pions which are produced in energetic proton interactions with ambient protons. An estimate of VHE emission from this process is obtained by an approximate, semi-quantitative calculation, and also by a detailed numerical treatment based on a convection-diffusion model for energetic electron and proton propagation and energy losses. All relevant hadronic and leptonic processes are considered, gauged by the measured synchrotron radio emission from the inner disk region. We estimate an integrated flux f(>100 GeV) 2E-12 1/(cm^2 s), possibly detectable by the current northern-hemisphere imaging air Cherenkov telescopes, MAGIC and VERITAS, and a good candidate for detection with the upcoming MAGIC II telescope. We also estimate f(>100 MeV) E-8 1/(cm^2 s), a level of emission that can be detected by GLAST/LAT based on the projected sensitivity for a one-year observation.

A2163 was observed by the RXTE satellite for 530 ks during a 6 month period starting in August 2004. The cluster primary emission is from very hot intracluster gas with kT~15 keV, but this component does not by itself provide the best fitting model. A secondary emission component is quite clearly needed, and while this could also be thermal at a temperature significantly lower than kT~15 keV, the best fit (to the combined PCA and HEXTE datasets) is obtained with a power law secondary spectral component. The deduced parameters of the non-thermal (NT) emission imply a significant fractional flux amounting to ~25% of the integrated 3-50 keV emission. NT emission is expected given the intense level of radio emission, most prominently from a large extended (`halo') central region of the cluster. Interpreting the deduced NT emission as Compton scattering of the radio-emitting relativistic electrons by the CMB, we estimate the volume-averaged value of the magnetic field in the extended radio region to be B=0.4+/-0.2 microG.

There are some basic differences between the observed properties of galaxies and clusters and the predictions from current hydrodynamical simulations. These are particularly pronounced in the central regions of galaxies and clusters. The popular NFW (Navarro, Frenk, and White) profile, for example, predicts a density cusp at the center, a behavior that (unsurprisingly) has not been observed. While it is not fully clear what are the reasons for this discrepancy, it perhaps reflects (at least partly) insufficient spatial resolution of the simulations. In this paper we explore a purely phenomenological approach to determine dark matter density profiles that are more consistent with observational results. Specifically, we deduce the gas density distribution from measured X-ray brightness profiles, and substitute it in the hydrostatic equilibrium equation in order to derive the form of dark matter profiles. Given some basic theoretical requirements from a dark matter profile, we then consider a number of simple profiles that have the desired asymptotic form. We conclude that a dark matter density profile of the form 1/(1+r/r_a)^3 is most consistent with current observational results.