This paper presents a study of the specific merger rate as a function of group membership, local environment, and redshift in a very large, $500h^{-1} Mpc$, cosmological N-body simulation, the \textit{Millennium Simulation}. The goal is to provide environmental diagnostics of major merger populations in order to test simulations against observations and provide further constraints on major merger driven galaxy evolution scenarios. A halo sample is defined using the maximum circular velocity, which is both well defined for subhalos and closely correlated with galaxy luminosity. Subhalos, including the precursors of major mergers, are severely tidally stripped. Major mergers between subhalos are therefore extremely rare. Tidal stripping also suppresses dynamical friction, resulting in long major merger time scales when the more massive halo does not host other subhalos. In contrast, when other subhalos are present major merger time scales are several times shorter. This enhancement is likely due to inelastic unbound collisions between subhalos. Following these results, we predict that major mergers in group environments are dominated by mergers involving the central galaxy, that the specific merger rate is suppressed in groups, and that the frequency of fainter companions is enhanced for mergers and their remnants. We also observe an `assembly bias' in the major merger rate in that mergers of galaxy-like halos are slightly suppressed in overdense environments while mergers of group-like halos are slightly enhanced. A dynamical explanation for this trend is advanced which calls on both tidal effects and interactions between bound halos beyond the virial radii of locally dynamically dominant halos.

I review some results on the galaxy/light-mass connection obtained by dissipationless simulations in combination with a simple, non-parametric model to connect halo circular velocity to the luminosity of the galaxy they would host. I focus on the galaxy-mass correlation and mass-to-light ratios obtained from galaxy up to cluster scales. The predictions of this simple scheme are shown to be in very good agreement with SDSS observations.

The last decade has seen applications of Adaptive Mesh Refinement (AMR) methods for a wide range of problems from space physics to cosmology. With the advent of these methods, in which space is discretized into a mesh of many individual cubic elements, the contemporary analog of the extensively studied line radiative transfer (RT) in a semi-infinite slab is that of RT in a cube. In this study we provide an approximate solution of the RT equation, as well as analytic expressions for the probability distribution functions (pdfs) of the properties of photons emerging from a cube, and compare them with the corresponding slab problem. These pdfs can be used to perform fast resonant-line RT in optically thick AMR cells where, otherwise, it could take unrealistically long times to transfer even a handful of photons.

We develop a Ly-alpha radiative transfer (RT) Monte Carlo code for cosmological simulations.High resolution,along with appropriately treated cooling can result in simulated environments with very high optical depths.Thus,solving the Ly-alpha RT problem in cosmological simulations can take an unrealistically long time.For this reason,we develop methods to speed up the Ly-alpha RT.With these accelerating methods,along with the parallelization of the code,we make the problem of Ly-alpha RT in the complex environments of cosmological simulations tractable.We test the RT code against simple Ly-alpha emitter models,and then we apply it to the brightest Ly-alpha emitter of a gasdynamics+N-body Adaptive Refinement Tree (ART) simulation at z~8.We find that recombination rather than cooling radiation Ly-alpha photons is the dominant contribution to the intrinsic Ly-alpha luminosity of the emitter,which is ~4.8x10e43 ergs/s.The size of the emitter is pretty small,making it unresolved for currently available instruments.Its spectrum before adding the Ly-alpha Gunn-Peterson absorption (GP) resembles that of static media,despite some net inward radial peculiar motion.This is because for such high optical depths as those in ART simulations,velocities of order some hundreds km/s are not important.We add the GP in two ways.First we assume no damping wing,corresponding to the situation where the emitter lies within the HII region of a very bright quasar,and second we allow for the damping wing.Including the damping wing leads to a maximum line brightness suppression by roughly a factor of ~62.The line fluxes,even though quite faint for current ground-based telescopes,should be within reach for JWST.

We investigate the properties of the pairwise velocity dispersion as a function of galaxy luminosity in the context of a halo model. We derive the distribution of velocities of pairs at a given separation taking into account both one-halo and two-halo contributions. We show that pairwise velocity distribution in real space is a complicated mixture of host-satellite, satellite-satellite and two-halo pairs. The peak value is reached at around 1$h^{-1}$Mpc and does not reflect the velocity dispersion of a typical halo hosting these galaxies, but is instead dominated by the satellite-satellite pairs in high mass clusters. This is true even for cross-correlations between bins separated in luminosity. As a consequence the velocity dispersion at a given separation can decrease with luminosity, even if the underlying typical halo host mass is increasing, in agreement with recent observations. We compare our findings to numerical simulations and find a good agreement. Numerical simulations also suggest a luminosity dependent velocity bias, which depends on the subhalo mass. We develop models of the auto- and cross-correlation function of luminosity subsamples of galaxies in the observable $r_\proj - \pi$ space and calculate the inferred velocity dispersion as a function of wave vector if dispersion model is fit to the redshift space power spectrum. We find that so derived pairwise velocity dispersion also exhibits a bump at $k\sim 1 h/{\rm Mpc}$.

Galaxy-galaxy lensing is a powerful probe of the relation between galaxies and dark matter halos, but its theoretical interpretation requires a careful modeling of various contributions, such as the contribution from central and satellite galaxies. For this purpose, a phenomenological approach based on the halo model has been developed, allowing for fast exploration of the parameter space of models. In this paper, we investigate the ability of the halo model to extract information from the g-g weak lensing signal by comparing it to high-resolution dissipationless simulations that resolve subhalos. We find that the halo model reliably determines parameters such as the host halo mass of central galaxies, the fraction of galaxies that are satellites, and their radial distribution inside larger halos. If there is a significant scatter present in the central galaxy host halo mass distribution, then the mean and median mass of that distribution can differ significantly from one another, and the halo model mass determination lies between the two. This result suggests that when analyzing the data, galaxy subsamples with a narrow central galaxy halo mass distribution, such as those based on stellar mass, should be chosen for a simpler interpretation of the results.

We analyze the mass accretion histories (MAHs) and density profiles of cluster-size halos in a flat LCDM cosmology.To fit the MAHs of all systems we need to generalize the MAH fit found in previous systematic studies of predominantly galactic halos.We find that the concentration of the density distribution is tightly correlated with the halo MAH and its formation redshift. During the early period of fast mass growth the concentration remains roughly constant and low (~3-4), while during the slow accretion phase it increases with decreasing redshift as 1/(1+z).We consider fits of three widely discussed analytic density profiles.We find that there is no unique best fit for all clusters.At the same time, if a cluster is best fit by a particular analytic profile at z=0, the same is usually true at earlier epochs out to z ~ 1-2.The local logarithmic slope of the density profiles at 3% of the virial radius ranges from -1.2 to -2.0.In addition, the logarithmic slope becomes shallower with decreasing radius without reaching an asymptotic value down to the smallest resolved scale (<1% of the virial radius).During the early MAH period of rapid mass growth the density profiles can be well described by a single power law rho(r) ~ r^-a with a ~ 1.5-2.The relatively shallow power law slopes result in low concentrations at these stages of evolution, as the scale radius where the density profiles reaches the slope of -2 is at large radii.This indicates that the inner power law like density distribution of halos is built up during the periods of rapid mass accretion and active merging, while the outer steeper profile is formed when the mass accretion slows down.

We present a new non-parametric method for determining mean 3D density and mass profiles from weak lensing measurements around stacked samples of galaxies or clusters, that is, from measurement of the galaxy-shear or cluster-shear correlation functions. Since the correlation function is statistically isotropic, this method evades problems, such as projection of large-scale structure along the line of sight or halo asphericity, that complicate attempts to infer masses from weak lensing measurements of individual objects. We demonstrate the utility of this method in measuring halo profiles, galaxy-mass and cluster-mass cross-correlation functions, and cluster virial masses. We test this method on an N-body simulation and show that it correctly and accurately recovers the 3D density and mass profiles of halos. We find no evidence of problems due to a mass sheet degeneracy in the simulation. Cross-correlation lensing provides a powerful method for calibrating the mass-observable relation for use in measurement of the cluster mass function in large surveys. It can also be used on large scales to measure and remove the halo bias and thereby provide a direct measurement of Omega_m * sigma_8.

I present results with respect to the gamma-ray and synchrotron fluxes due to neutralino annihilation in the Large Magellanic Cloud based on Tasitsiomi et al., Astrop. Phys. 21, 637 (2004).

We analyze the mass accretion histories (MAHs) and density profiles of cluster- size halos with virial masses of 0.6-2.5x10^14/h Msun in a flat LCDM cosmology. In agreement with previous studies,we find that the concentration of the density distribution is tightly correlated with the halo's MAH and its formation redshift.During the period of fast mass growth the concentration remains approximately constant and low c_v ~ 3-4,while during the slow accretion stages it increases with decreasing redshift as c_v ~ (1+z)^-1.We consider fits of three widely discussed analytic density profiles to the simulated clusters focusing on the most relaxed inner regions.We find that there is no unique best fit analytic profile for all the systems.If,however,a cluster is best fit by a particular analytic profile at z=0,the same is usually true at earlier epochs out to z ~ 1-2.The local logarithmic slope of the density profiles at 3% of the virial radius ranges from -1.2 to -2.0,a remarkable diversity for the relatively narrow mass range of our cluster sample.For all the studied clusters the logarithmic slope becomes shallower with decreasing radius without reaching an asymptotic value down to the smallest resolved scale < ~1% of the virial radius.We do not find a clear correlation of the inner slope with the formation redshift or the shape of the halo's MAH.We do find, however,that during the period of rapid mass growth the density profiles can be well described by a single power law rho(r) ~ r^-gamma with gamma ~ 1.5-2.The relatively shallow power law slopes result in low concentrations at these stages of evolution.This indicates that the inner power law like density distribution of halos is built up during the periods of rapid mass accretion and active merging,while outer steeper profile is formed when the mass accretion is slow.

We use high-resolution dissipationless simulations of the concordance flat LCDM model to make predictions for the galaxy--mass (GM) correlations and compare them to the recent SDSS weak lensing measurements.We use a simple observationally motivated scheme to assign luminosities and colors to the halos.This allows us to closely match the selection criteria used to define observational samples.The simulations reproduce the observed GM correlation function and its observed dependencies on luminosity and color.The luminosity dependence of the correlation function is primarily determined by the changing relative contribution of central and satellite galaxies at different luminosities. The color dependence of the GM correlations reflects the difference in the typical environments of blue and red galaxies. We also find agreement between the predicted and observed cross-bias, b_x=b/r,at all probed scales.The GM correlation coefficient, r, is close to unity on scales >1/h Mpc.The cross bias is thus expected to measure the actual bias of galaxy clustering on these scales.The aperture mass-to-light ratio is independent of galaxy color.The aperture mass scales approximately linearly with luminosity at L_r>10^{10}h^{-2} Lsun, while at lower luminosities the scaling is shallower: L_r^{0.5}. We show that most of the luminous galaxies (M_r<-21) are near the centers of their halos and their GM correlation function at r<100/h kpc can therefore be interpreted as the average dark matter density profile of these galaxies. We find that for galaxies in a given narrow luminosity range, there is a broad and possibly non-gaussian distribution of halo virial masses. Therefore, the average relation between mass and luminosity derived from the weak lensing analyses should be interpreted with caution.

We calculate the expected flux of gamma-ray and radio emission from the LMC due to neutralino annihilation. Using rotation curve data to probe the density profile and assuming a minimum disk, we describe the dark matter halo of the LMC using models predicted by N-body simulations. We consider a range of density profiles including the NFW profile, a modified NFW profile proposed by Hayashi et al.(2003) to account for the effects of tidal stripping, and an isothermal sphere with a core. We find that the gamma-ray flux expected from these models may be detectable by GLAST for a significant part of the neutralino parameter space. The prospects for existing and upcoming Atmospheric Cherenkov Telescopes are less optimistic, as unrealistically long exposures are required for detection. However, the effects of adiabatic compression due to the baryonic component may improve the chances for detection by ACTs. The maximum flux we predict is well below EGRET's measurements and thus EGRET does not constrain the parameter space. The expected synchrotron emission generally lies below the observed radio emission from the LMC in the frequency range of 19.7 to 8550 MHz. As long as <sigma v> <2x 10^-26 cm^3 s^-1 for a neutralino mass of 50 GeV, the observed radio emission is not primarily due to neutralinos and is consistent with the assumption that the main source is cosmic rays. We find that the predicted fluxes, obtained by integrating over the entire LMC, are not very strongly dependent on the inner slope of the halo profile, varying by less than an order of magnitude for the range of profiles we considered.

We discuss the detectability of dark matter clumps in the Milky Way halo due to neutralino annihilation. We then focus on a known ``clump'', the Large Magellanic Cloud.

High resolution N-body simulations have revealed the survival of considerable substructure within galactic halos. Assuming that the predicted dark matter clumps are composed of annihilating neutralinos, we examine their detectability via Atmospheric Cherenkov Telescopes (ACTs). Depending on their density profile, individual neutralino clumps should be observable via their gamma-ray continuum and line emissions. We find that the continuum signal is the most promising signal for detecting a neutralino clump, being significantly stronger than the monochromatic signals. Limits from the line detectability can help lift degeneracies in the supersymmetric (SUSY) parameter space. We show that by combining the observations of different mass clumps, ACTs can explore most of the SUSY parameter space. ACTs can play a complementary role to accelerator and gamma-ray satellite limits by exploring relatively large neutralino masses and less concentrated clumps. We develop a strategy for dark matter clump studies by future ACTs based on VERITAS specifications and encourage the development of techniques to identify primaries. This can reduce the background by an order of magnitude.

The Cold Dark Matter model describes successfully our Universe on large scales, as has been verified by a wide range of observations. Nevertheless, in the last years, and especially with the advances in numerical simulations, a number of apparent inconsistencies arose between what is observed and what the CDM model predicts on small scales. In this work, the current status of observations on galactic and subgalactic scales is reviewed. Furthermore, theory and observation are brought together in order to reveal the nature of the inconsistencies and, consequently, to reveal their severity. Lastly, the progress towards the resolution of each one of these conflicts is briefly reviewed.