The relic abundance of thermal dark matter particles is generally assumed to be inversely proportional to their annihilation rate, which is therefore constrained by the present matter density, <sigma v> ~ 10^{-26} Omega_{dm}^{-1} cm^3 sec^{-1}. Here we point out that much lower values of <sigma v> are possible for heavy dark matter candidates (m < 10 TeV) that couple to other particle species through the electroweak force. With heavy dark matter particles present the early universe may evolve according to the following scenario. After an early entry into matter-dominated phase, dark matter particles form self-gravitating microhalos. Collisional interaction between dark matter particles and the surrounding radiation field eventually leads to microhalos gravothermal collapse and annihilation of most dark matter particles. For sufficiently heavy dark matter candidates (m < 10 TeV) the universe can return to radiation-dominated phase before the nucleosynthesis and thereafter follow the "standard" scenario.

We reexamine the limits on charged dark matter particles. We show that if their mass and charge fall in the range 100(q_X/e)^2< m_X < 10^8(q_X/e) TeV, then magnetic fields prevent particles in the halo from entering the galactic disk, while those initially trapped inside are accelerated through the Fermi mechanism and ejected within about 0.1-1 Gyrs. Consequently, previous constraints on charged dark matter based on terrestrial non-observation are invalid within that range. Further, we find that charged massive particles may simultaneously solve several long-standing astrophysical problems, including the underabundance of dwarf galaxies, the shallow density profiles in the cores of the LSB galaxies, the absence of cooling flows in the cores of galaxy clusters, and several others.

We reexamine the impact of dark matter (DM) annihilation on the intergalactic medium, taking into account the clumping of DM particles. We find that energy injection from the annihilation of the thermal relic DM particles may significantly raise the gas temperature at high redshifts and leave a strong imprint on the cosmological 21-cm signal, provided the particle mass is below ~1 TeV. Further, we find that while the energy injection from DM annihilation could not alone complete the reionization of the Universe, it could make a significant contribution to the electron optical depth.

During the epoch of reionization, the formation of radiation sources is accompanied by the growth of a H- photodissociating flux. We estimate the impact of this flux on the formation of molecular hydrogen and cooling in the first galaxies, assuming different types of radiation sources (e.g. Pop II and Pop III stars, miniquasars). We find that H- photodissociation reduces the formation of H2 molecules by a factor of ~1+1000k_s*x/(f_esc*delta), where x is the mean ionized fraction in the IGM, f_esc is the fraction of ionizing photons that escape from their progenitor halos, delta is the local gas overdensity and k_s is an order unity constant which depends on the type of radiation source. By the time a significant fraction of the universe becomes ionized, H- photodissociation may significantly reduce the H2 abundance and, with it, the primordial star formation rate, delaying the progress of reionization.

During the epoch of reionization the 21cm signal is sensitive to the scattering rate of the ultraviolet photons, redshifting across the Lyman_alpha resonance. Here we calculate the photon scattering rate profile for a single ultraviolet source. After taking into account previously neglected natural broadening of the resonance line, we find that photons approach the resonance frequency and experience most scatterings at a significantly smaller distance from the source than naively expected r=(dnu/nu_0)(c/H), where dnu=nu-nu_0 is the initial frequency offset, and the discrepancy increases as the initial frequency offset decreases. As a consequence, the scattering rate P(r) drops much faster with increasing distance than the previously assumed 1/r^2 profile. Near the source (r<1Mpc comoving), the scattering rate of photons that redshift into the Ly_alpha resonance converges to P(r) \propto r^{-7/3}. The scattering rate of Ly_alpha photons produced by splitting of photons that redshift into a higher resonance (Ly_gamma, Ly_delta, etc.) is only weakly affected by the radiative transfer, while the sum of scattering rates of Ly_alpha photons produced from all higher resonances also converges to P(r) \propto r^{-7/3} near the source. At 15<z<35, on scales of ~0.01-20Mpc/h (comoving), the total scattering rate of Ly_alpha photons from all Lyman resonances is found to be higher by a factor of ~1+0.3[(1+z)/20]^{2/3} than obtained without full radiative transfer. Consequently, during the early stage of reionization, the differential brightness of 21cm signal against the cosmic microwave background is also boosted by a similar factor.

During the epoch of reionization a large number of photons were produced with frequencies below the hydrogen Lyman limit. After redshifting into the closest resonance, these photons underwent multiple scatterings with atoms. We examine the effect of these scatterings on the temperature of the neutral intergalactic medium (IGM). Continuum photons, emitted between the Ly_alpha and Ly_gamma frequencies, heat the gas after being redshifted into the H Ly_alpha or D Ly_beta resonance. By contrast, photons emitted between the Ly_gamma and Ly-limit frequencies, produce effective cooling of the gas. Prior to reionization, the equilibrium temperature of ~100 K for hydrogen and helium atoms is set by these two competing processes. At the same time, Ly_beta resonance photons thermally decouple deuterium from other species, raising its temperature as high as 10^4 K. Our results have important consequences for the cosmic 21-cm background and the entropy floor of the early IGM which can affect star formation and reionization.

Recent WMAP observations indicate that the primordial abundance of lithium should have been 2-4 times higher than the value measured in low-metallicity stars. We show that this discrepancy can be explained by the process of ambipolar diffusion in protostellar clouds, which leads to depletion of elements with low ionization potential, such as lithium and sodium. In high-metallicity stars, [Fe/H]>~-1, the depletion of lithium due to this process is found to be negligible, but for [Fe/H]<~-1.5 the lithium abundance falls by a factor of ~3, consistent with WMAP observations.

We derive the density profile of cold dark matter halos using a self-similar accretion model. We show that if the clumpiness of the infalling matter is taken into account, then the inner density slope, d(log{rho})/d(log{r}), is close to -1. Compared with the density profiles predicted by different numerical simulations, we find that outside ~0.1% of the virial radius, our solutions agree best with the fitting formula proposed by Navarro et al. (2004), d(ln{rho})/d(ln{r})=-2(r/R)^0.17, where R is a characteristic radius, inside which the density profile becomes shallower than isothermal.

At the beginning of the reionization epoch, radiation sources produce fluctuations in the redshifted 21-cm background. We show that different types of sources (such as miniquasars, Pop II and III stars, supernovae, etc.) produce distinct signatures in the 21-cm signal radial profiles and statistical fluctuations, through which they can be identified. Further, we show that the 21-cm signal from X-ray emitting sources is much easier to observe than was expected, due to a previously neglected pumping mechanism.

We present new analytic calculations of the coupling between ultraviolet resonance photons and the population of the hyperfine states in the light elements (H, D, He3+) which include several previously neglected physical processes. Among these are the backreaction of resonant scattering on the pumping radiation, the scattering of Ly_beta photons and the effect of local departure from pure Hubble flow. The application of the new treatment to the redshifted hydrogen 21 and deuterium 92 cm lines from the high-redshift universe results in an amplitude correction of up to an order of magnitude. We further show that the standard assumption that ultraviolet pumping drives the spin temperature towards the kinetic temperature does not hold for deuterium, whose spin temperature is generally negative.

We show that element diffusion can produce large fluctuations in the initial helium abundance of the star-forming clouds. Diffusion time-scale, which in stars is much larger than the Hubble time, can fall below 10^8 years in the neutral gas clouds dominated by collisionless dark matter or with dynamically important radiation or magnetic pressure. Helium diffusion may therefore explain the recent observations of globular clusters, which are inconsistent with initially homogeneous helium distribution.

Protons gain energy in short range collisions with heavier dark matter particles (DMPs) of comparable velocity dispersion. We examine the conditions under which the heating of baryons by scattering off DMPs can offset radiative cooling in the cores of galaxy clusters. Collisions with a constant cross section independent of the relative velocity of the colliding particles, cannot produce stable thermal balance. In this case, avoiding an unrealistic increase of the central temperatures yields the upper bound on the cross-section, \sigma_xp<10^-25 cm^2 (m_x/m_p), where m_x and m_p are the DMP and proton mass, respectively. A stable balance, however, can be achieved for a power law dependence on the relative velocity, V, of the form \sigma_xp \propto V^a with a<-3. An advantage of this heating mechanism is that it preserves the metal gradients observed in clusters.

We examine the process of element segregation by gravity in giant elliptical galaxies and X-ray clusters. We solve the full set of flow equations developed by Burgers for a multi-component fluid, under the assumption that magnetic fields slow the diffusion by a constant factor F_B. Compared to the previous calculations that neglected the residual heat flow terms, we find that diffusion is faster by ~20%. In clusters we find that the diffusion changes the local abundance typically by factors of 1+0.3(T/10^8 K)^1.5/F_B and 1+0.15(T/10^8 K)^1.5/F_B for helium and heavy elements, respectively, where T is the gas temperature. In elliptical galaxies, the corresponding factors are 1+0.2(T/10^7 K)^1.5/F_B and 1+0.1(T/10^7 K)^1.5/F_B, respectively. If the suppression factor F_B is modest, diffusion could significantly affect observational properties of hot X-ray clusters and cD galaxies. In particular, diffusion steepens the baryon distribution, increases the total X-ray luminosity, and changes the spectrum and evolution of stars that form out of the helium-rich gas. Detection of these diffusion signatures would allow to gauge the significance of the magnetic fields, which also inhibit thermal heat conduction as a mechanism for eliminating cooling flows in the same environments.

Magnetic fields are generated in ionized objects rotating with respect to a radiation background. Based on conservation of canonical ion momentum, it has been previously suggested that even if the radiation intensity is unrestricted the maximum field strength generated this way is ~(10^{-4}\Omega) Gauss, where \Omega is the initial angular velocity of the object. We show that this limit is valid only for an object made of fully ionized and optically thin plasma. The limit can be relaxed by orders of magnitude in the presence of a high neutral fraction or if the plasma is coupled to other forms of matter such as stars or compact clouds.

We revisit the process of gravitational sedimentation of helium and heavy elements in the intra-cluster medium. We find that helium applies an inward drag force on heavy elements, boosting their sedimentation speed to nearly half its own. This speed is almost independent of the mass and the electric charge of heavy elements. In the absence of small-scale magnetic fields, helium sedimentation can increase the He/H abundance ratio in the cores of hot clusters by three orders of magnitude. It also steepens the baryonic density profile yielding a higher X-ray luminosity, which offers an explanation of the observed luminosity-temperature relation. If the primordial He/H ratio is assumed, then the gas density inferred from the observed X-ray emissivity might be underestimated by 30% in the cores of clusters and overestimated by 7% in the outer regions. The dark matter density on the other hand might be overestimated by a factor of 8/3 in the cores and underestimated by 18% in outer regions.

We derive similarity solutions for the expansion of negative initial density perturbations \delta M/M \propto r^-s (s>0) in an Einstein-de Sitter universe filled with collisional baryonic gas and collisionless matter. The solutions are obtained for planar, cylindrical, and spherical perturbations. For steep perturbations, central cavities surrounded by over-dense shells expanding as t^{2(s+1)/3s} develop in both matter components. We also consider the case when baryonic shells are driven by internal pressure. Before redshift \sim 8-10 baryonic shells cool by inverse Compton scattering with the cosmic background radiation and can form bound objects of mass \sim 10^11\Omega_b^{-1/2} M_{\odot} (1+z)^{-3/2}. Hot shocks of radius R cause a Compton y-distortion of order y=10^{-5}\Omega_b (R/10 Mpc)^2 f_F, where f_F is the volume filling factor of these shells.

Similarity solutions are found for the adiabatic collapse of density perturbations $\delta M/M \propto r^{-s}$ $(s>0)$ in a flat universe containing collisional gas only. The solutions are obtained for planar, cylindrical, and spherical perturbations with zero initial pressure. For adiabatic index $\gamma\ge 4/3$, a shock develops at a fixed fraction of the current turnaround distance. Near the center of a spherical perturbations with $\gamma>4/3$ and $s>1/2$, the gas is in quasi-hydrostatic equilibrium (pressure supported) and has an asymptotic power law density profile, $\rho\sim r^{-3s/(s+1)}$, independent of $\gamma$. For $s\le 1/2$, the profile depends on $\gamma$, the pressure is finite, the temperature decreases inward, and gravity dominates pressure causing a continuous inward flow. Although for $1/2<s<2$ the temperature decreases at the center, the gas is pressure supported. The pressure is finite in cylindrical perturbations for $s\le 2(\gamma-1)/(3\gamma-4)$, and in planar perturbations for any $s>0$. We also derive the asymptotic behaviour of the gas variables near the center in a universe dominated by collisionless matter. In such a universe, the gas in a spherical perturbation with $s<2$ cannot be pressure supported and the temperature approaches a constant near the center. The solutions and the asymptotic behaviour are relevant for modelling the gas distribution in galaxy clusters and pancake-like superclusters, and determining the structure of haloes of self-interacting dark matter with large interaction cross section.