Dark matter halo profiles in scale-free cosmologies

We explore the dependence of the central logarithmic slope of dark matter halo density profiles α on the spectral index n of the linear matter power spectrum P ( k ) using cosmological N -body simulations of scale-free models [i.e. P ( k ) ∝ k ⁿ ]. These simulations are based on a set of clear, reproducible and physically motivated criteria that fix the appropriate starting and stopping times for runs, and allow one to compare haloes across models with different spectral indices and mass resolutions. For each of our simulations we identify samples of well-resolved haloes in dynamical equilibrium and we analyse their mass profiles. By parametrizing the mass profile using a 'generalized' Navarro, Frenk & White profile in which the central logarithmic slope α is allowed to vary while preserving the r ⁻³ asymptotic form at large radii, we obtain preferred central slopes for haloes in each of our models. There is a strong correlation between α and n , such that α becomes shallower as n becomes steeper. However, if we normalize our mass profiles by r ₋₂, the radius at which the logarithmic slope of the density profile is −2, we find that these differences are no longer present. This is apparent if we plot the maximum slope     as a function of r / r ₋₂– we find that the profiles are similar for haloes forming in different n models. This reflects the importance of concentration, and reveals that the concentrations of haloes forming in steep- n cosmologies tend to be smaller than those of haloes forming in shallow- n cosmologies. We conclude that there is no evidence for convergence to a unique central asymptotic slope, at least on the scales that we can resolve.     

Voids in a ΛCDM universe

We study the formation and evolution of voids in the dark matter distribution using various simulations of the popular Λ cold dark matter cosmogony. We identify voids by requiring them to be regions of space with a mean overdensity of −0.8 or less – roughly the equivalent of using a spherical overdensity group finder for haloes. Each of the simulations contains thousands of voids. The distribution of void sizes in the different simulations shows good agreement when differences in particle and grid resolution are accounted for. Voids very clearly correspond to minima in the smoothed initial density field. Apart from a very weak dependence on the mass resolution, the rescaled mass profiles of voids in the different simulations agree remarkably well. We find a universal void mass profile of the form  ρ(< r )/ρ( r _eff) ∝ exp[( r / r _eff)^α]  , where r _eff is the effective radius of a void and  α∼ 2  . The mass function of haloes in voids is steeper than that of haloes that populate denser regions. In addition, the abundances of void haloes seem to evolve somewhat more strongly between redshifts ∼1 and 0 than the global abundances of haloes.     

The distribution function of dark matter in massive haloes

We study the distribution function (DF) of dark matter particles in haloes of mass range  10¹⁴–10¹⁵ M_⊙  . In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard Λ cold dark matter (ΛCDM) model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L , a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form     . This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (β₀ and  β_∞  ) and the scale of transition between them ( L ₀). The energy part   f _E ( E )  is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase-space density   Q =ρ( r )/σ³( r )  .     

Hot stars in old stellar populations: a continuing need for intermediate ages

We present numerical investigations into the formation of massive stars from turbulent cores of density structure  ρ∝ r ^−1.5  . The results of five hydrodynamical simulations are described, following the collapse of the core, fragmentation and the formation of small clusters of protostars. We generate two different initial turbulent velocity fields corresponding to power-law spectra   P ∝ k ⁻⁴  and   P ∝ k ^−3.5  , and we apply two different initial core radii. Calculations are included for both completely isothermal collapse, and a non-isothermal equation of state above a critical density  (10⁻¹⁴ g cm⁻³)  . Our calculations reveal the preference of fragmentation over monolithic star formation in turbulent cores. Fragmentation was prevalent in all the isothermal cases. Although disc fragmentation was largely suppressed in the non-isothermal runs due to the small dynamic range between the initial density and the critical density, our results show that some fragmentation still persisted. This is inconsistent with previous suggestions that turbulent cores result in the formation of a single massive star. We conclude that turbulence cannot be measured as an isotropic pressure term.     

Dark matter subhaloes in numerical simulations

We use cosmological Λ cold dark matter (CDM) numerical simulations to model the evolution of the substructure population in 16 dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10¹¹ to  10¹⁴ M_⊙  at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation  d n /d v = (1/8)( v _cmax/ v _cmax,host)⁻⁴  . Halo to halo variance in the VDF is a factor of roughly 2 to 4. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for   r ≳ 0.3 r _vir  , where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease towards the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias towards more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power-law cosmologies of   P ∝ kⁿ   , and demonstrate that a steeper spectral index, n , results in significantly less substructure.     

Deep galaxy counts, extragalactic background light and the stellar baryon budget

We assess the constraints imposed by the observed extragalactic background light (EBL) on the cosmic history of star formation and the stellar-mass density today. The logarithmic slope of the galaxy number–magnitude relation from the Southern Hubble Deep Field imaging survey is flatter than 0.4 in all seven UBVIJHK optical bandpasses, i.e. the light from resolved galaxies has converged from the UV to the near-IR. We find a lower limit to the surface brightness of the optical extragalactic sky of about 15 nW m⁻² sr⁻¹, comparable to the intensity of the far-IR background from COBE data. Assuming a Salpeter initial mass function with a lower cut-off consistent with observations of M subdwarf disc stars, we set a lower limit of Ω_g+s h ²>0.0013  I ₅₀ to the visible (processed gas + stars) mass density required to generate an EBL at a level of 50  I ₅₀ nW m⁻² sr⁻¹; our 'best-guess' value is Ω_g+s h ²≈0.0031  I ₅₀. Motivated by the recent microlensing results of the MACHO collaboration, we consider the possibility that massive dark haloes around spiral galaxies are composed of faint white dwarfs, and show that only a small fraction (≲5 per cent) of the nucleosynthetic baryons can be locked in the remnants of intermediate-mass stars forming at z _F≲5, as the bright early phases of such haloes would otherwise overproduce the observed EBL.     

Degeneracy of semi-analytic solution for axially symmetric lens equations

We derive a simple semi-analytical approximation for lens equations with an arbitrary radially symmetric mass density ρ( r ), when   r /ξ₀≪ 1  and ξ₀ is the scalelength of the density profile. At the strong lensing regime, which is mostly constrained by the inner part of the mass density profile, we assume ρ∝ r ^α.
A dark matter (DM) haloes (GNFW model) are parametrized through a shape parameter α, a concentration parameter c ₁ and the total mass M . We apply our semi-analytical model to show how the solutions of the axially symmetric lens equations are degenerated in respect to the parameters α and c ₁.
In the case of an asymmetric dual image lens system, similar effective degeneracy is produced when the geometry of the lens is relaxed. Because it is impossible to determine the exact location of the source image, a family of solutions is acquired when the mass of the lens object and location of the observed images are fixed.
Our results indicate that the amount of degeneration is only weakly affected by the asymmetry in the lensing geometry set-up, e.g. the observational effective degeneracy is very close to the true physical degeneracy of the Einstein ring solutions. Basically with high-enough values for the concentration parameter, the degeneracy spawns the whole range for the shape parameter  α=[−2.0, −1.0]  .     

The circumstellar medium of the peculiar supernova SN1997ab

We report the detection of the slow-moving wind into which the compact supernova remnant SN 1997ab is expanding. Echelle spectroscopy provides clear evidence for a well-resolved narrow (full width at zero intensity, FWZI ∼180 km s⁻¹) P Cygni profile, both in Hα and Hβ, superimposed on the broad emission lines of this compact supernova remnant. From theoretical arguments we know that the broad and strong emission lines imply a circumstellar density ( n  ≥ 10⁷ cm⁻³). This, together with our detection, implies a massive and slow stellar wind experienced by the progenitor star shortly prior to the explosion.     

Constraints on dark matter physics from dwarf galaxies through galaxy cluster haloes

One of the predictions of the standard cold dark matter model is that dark haloes have centrally divergent density profiles. An extensive body of rotation curve observations of dwarf and low surface brightness galaxies shows the dark haloes of those systems to be characterized by soft constant-density central cores. Several physical processes have been proposed to produce soft cores in dark haloes, each one with different scaling properties. With the aim of discriminating among them we have examined the rotation curves of dark-matter-dominated dwarf and low surface brightness galaxies and the inner mass profiles of two clusters of galaxies lacking a central cD galaxy and with evidence of soft cores in the centre. The core radii and central densities of these haloes scale in a well-defined manner with the depth of their potential wells, as measured through the maximum circular velocity. As a result of our analysis we identify self-interacting cold dark matter as a viable solution to the core problem, where a non-singular isothermal core is formed in the halo centre surrounded by a Navarro, Frenk & White profile in the outer parts. We show that this particular physical situation predicts core radii in agreement with observations. Furthermore, using the observed scalings, we derive an expression for the minimum cross-section ( σ ) which has an explicit dependence with the halo dispersion velocity ( v ). If m _x is the mass of the dark matter particle: σ m _x ≈4×10⁻²⁵ (100 km s⁻¹  v ⁻¹) cm² GeV⁻¹.     

A spectroscopic study of IRAS F10214 + 4724

The z  = 2.286  IRAS galaxy F10214 + 4724 remains one of the most luminous galaxies in the Universe, despite its gravitational lens magnification. We present optical and near-infrared spectra of F10214 + 4724, with clear evidence for three distinct components: lines of width ∼ 1000 km s⁻¹ from a Seyfert 2 nucleus; ≲ 200 km s⁻¹ lines which are likely to be associated with star formation; and a broad (∼ 4000 km s⁻¹) C  III ] 1909-Å emission line which is blueshifted by ∼ 1000 km s⁻¹ with respect to the Seyfert 2 lines. Our study of the Seyfert 2 component leads to several new results. (i) From the double-peaked structure in the Lyα line, and the lack of Lyβ, we argue that the Lyα photons have emerged through a neutral column of N _H ∼ 2.5 × 10²⁵ m⁻², possibly located within the AGN narrow-line region, as proposed for several high-redshift radio galaxies. (ii) The resonant O  VI 1032, 1036-Å doublet (previously identified as Lyβ) is in an optically thick (1:1) ratio. At face value this implies an extreme density ( n _e ∼ 10¹⁷ m⁻³) more typical of broad-line region clouds. However, we attribute this instead to the damping wings of Lyβ from the resonant absorption. (iii) A tentative detection of He  II 1086 suggests little extinction in the rest frame ultraviolet.     

Virial shocks in galactic haloes?

We investigate the conditions for the existence of an expanding virial shock in the gas falling within a spherical dark matter halo. The shock relies on pressure support by the shock-heated gas behind it. When the radiative cooling is efficient compared with the infall rate, the post-shock gas becomes unstable; it collapses inwards and cannot support the shock. We find for a monatomic gas that the shock is stable when the post-shock pressure and density obey     . When expressed in terms of the pre-shock gas properties at radius r it reads as  ρ r Λ( T )/ u ³ < 0.0126  , where ρ is the gas density, u is the infall velocity and Λ( T ) is the cooling function, with the post-shock temperature   T ∝ u ²  . This result is confirmed by hydrodynamical simulations, using an accurate spheri-symmetric Lagrangian code. When the stability analysis is applied in cosmology, we find that a virial shock does not develop in most haloes that form before   z ∼ 2  , and it never forms in haloes less massive than a few  10¹¹ M_⊙  . In such haloes, the infalling gas is not heated to the virial temperature until it hits the disc, thus avoiding the cooling-dominated quasi-static contraction phase. The direct collapse of the cold gas into the disc should have non-trivial effects on the star formation rate and on outflows. The soft X-ray produced by the shock-heated gas in the disc is expected to ionize the dense disc environment, and the subsequent recombination would result in a high flux of Lα emission. This may explain both the puzzling low flux of soft X-ray background and the Lα emitters observed at high redshift.     

Galactic halo cusp–core: tidal compression in mergers

We explain in simple terms how the build-up of dark haloes by merging compact satellites, as in the cold dark matter (CDM) cosmology, inevitably leads to an inner cusp of density profile  ρ∝ r ^−α  with  α≳ 1  , as seen in cosmological N -body simulations. A flatter halo core with  α < 1  exerts on the satellites tidal compression in all directions, which prevents the deposit of stripped satellite material in the core region. This makes the satellite orbits decay from the radius where  α∼ 1  to the halo centre with no local tidal mass transfer, and thus causes a rapid steepening of the inner profile to  α > 1  . These tidal effects, the resultant steepening of the profile to a cusp, and the stability of this cusp to tandem mergers with compact satellites are demonstrated using N -body simulations. The transition at  α∼ 1  is then addressed using toy models in the limiting cases of impulse and adiabatic approximations and using tidal radii for satellites on radial and circular orbits. In an associated paper, we address the subsequent slow convergence from either side to an asymptotic stable cusp with  α≳ 1  . Our analysis thus implies that an inner cusp is enforced when small haloes are typically more compact than larger haloes, as in the CDM scenario, such that enough satellite material makes it intact into the inner halo and is deposited there. We conclude that a necessary condition for maintaining a flat core, as indicated by observations, is that the inner regions of the CDM satellite haloes be puffed up by about 50 per cent such that when they merge into a larger halo they would be disrupted outside the halo core. This puffing up could be due to baryonic feedback processes in small haloes, which may be stimulated by the tidal compression in the halo cores.     

On the prospect of inferring the halo structure and the masses of dark objects through parallax microlensing

We study the proposed use of parallax microlensing in the direction of the Large Magellanic Cloud (LMC) to separate the effects of the mass function of dark massive halo objects (MHOs or 'machos') on the one hand, and their spatial distribution and kinematics on the other. This disentanglement is supposed to allow a much better determination of the two than could be achieved entirely on the basis of the durations of events. We restrict our treatment to the same class of power-law spherical models for the halo of MHOs studied in a previous paper by Marković 38 Sommer-Larsen, and assume that one can eliminate microlensing events caused by massive objects outside the halo (e.g., the LMC halo). Whereas the duration-based error in the average MHO mass, μ¯ ≡  M ¯/M, exceeds (at N  = 100 events) μ¯ by a factor of 2 or more, parallax microlensing remarkably brings it down to 15–20 per cent of μ¯, regardless of the shape of the mass function. In addition, the slope α of the mass function, d n /dμ ∝ μ^α, can be inferred relatively accurately (σ_α < 0.4) for a broader range, −3 < α < 0. The improvement in the inference of the halo structure is also significant: the index γ of the density profile ( ρ ∼  R ^−γ) can be obtained with the error σ_γ < 0.4. While in a typical situation the errors for the parameters specifying the velocity dispersion profile are of about the same magnitude as the parameters themselves, virtually all the uncertainty is 'concentrated' in linear combinations of the parameters that may have little influence on the profile, thus allowing its reasonably accurate inference.     

Survival and disruption of subsystems during a cold collapse

Cold collapse of a cluster composed of small identical clumps, each of which is in virial equilibrium, is considered. Since the clumps have no relative motion with respect to each other initially, the cluster collapses under its own gravity. At the first collapse of the cluster, most of the clumps are destroyed, but some survive. In order to find the condition for the clumps to survive, we made a systematic study in two-parameter space: the number of the clumps N _c and the size of the clump r _v . We obtained the condition N _c ≫ 1 and n _k  ≥ 1, where n _k is related to r _v and the initial radius of the cluster R _ini through the relation R _ini/ r _v  = 2 N ^{( n k +5)/6}_c. A simple analytical argument supports the numerical result. This n _k corresponds to the index of the power spectrum of the density fluctuation in the cosmological hierarchical clustering, and thus our result may suggest that in the systems smaller than 2/Ω h ²)Mpc, the first violent collapse is strong enough to sweep away all the substructures that exist before the collapse.     

On the distribution of haloes, galaxies and mass

The stochasticity in the distribution of dark haloes in the cosmic density field is reflected in the distribution function   P _V ( N _h| δ _m)  , which gives the probability of finding N _h haloes in a volume V with mass density contrast δ _m. We study the properties of this function using high-resolution N -body simulations, and find that   P _V ( N _h| δ _m)  is significantly non-Poisson. The ratio between the variance and the mean goes from ∼1 (Poisson) at  1+ δ _m≪1  to <1 (sub-Poisson) at  1+ δ _m∼1  to >1 (super-Poisson) at  1+ δ _m≫1  . The mean bias relation is found to be well described by halo bias models based on the Press–Schechter formalism. The sub-Poisson variance can be explained as a result of halo exclusion, while the super-Poisson variance at high δ _m may be explained as a result of halo clustering. A simple phenomenological model is proposed to describe the behaviour of the variance as a function of δ _m. Galaxy distribution in the cosmic density field predicted by semi-analytic models of galaxy formation shows similar stochastic behaviour. We discuss the implications of the stochasticity in halo bias to the modelling of higher order moments of dark haloes and of galaxies.     

Understanding the halo-mass and galaxy-mass cross-correlation functions

We use the Millennium Simulation (MS) to measure the cross-correlation between halo centres and mass (or equivalently the average density profiles of dark haloes) in a Lambda cold dark matter (ΛCDM) cosmology. We present results for radii in the range  10  h ⁻¹ kpc < r < 30  h ⁻¹ Mpc  and for halo masses in the range  4 × 10¹⁰ < M ₂₀₀ < 4 × 10¹⁴  h ⁻¹ M_⊙  . Both at   z = 0  and at   z = 0.76  these cross-correlations are surprisingly well fitted if the inner region is approximated by a density profile of NFW or Einasto form, the outer region by a biased version of the linear mass autocorrelation function, and the maximum of the two is adopted where they are comparable. We use a simulation of galaxy formation within the MS to explore how these results are reflected in cross-correlations between galaxies and mass. These are directly observable through galaxy–galaxy lensing. Here also we find that simple models can represent the simulation results remarkably well, typically to ≲10 per cent. Such models can be used to extend our results to other redshifts, to cosmologies with other parameters, and to other assumptions about how galaxies populate dark haloes. Our galaxy formation simulation already reproduces current galaxy–galaxy lensing data quite well. The characteristic features predicted in the galaxy–galaxy lensing signal should provide a strong test of the ΛCDM cosmology as well as a route to understanding how galaxies form within it.     

Global velocity field and bubbles in the blue compact dwarf galaxy Mrk 86

We have studied the velocity field of the blue compact dwarf galaxy Mrk 86 (NGC 2537) using data provided by 14 long-slit optical spectra obtained in 10 different orientations and positions. This kinematical information is complemented with narrow-band ([O  iii ]5007 Å and H α ) and broad-band ( B , V , Gunn r and K ) imaging. The analysis of the galaxy global velocity field suggests that the ionized gas could be distributed in a rotating inclined disc, with projected central angular velocity of Ω=34 km s⁻¹ kpc⁻¹. The comparison between the stellar, H  i and modelled dark matter density profile indicates that the total mass within its optical radius is dominated by the stellar component. Peculiarities observed in its velocity field can be explained by irregularities in the ionized gas distribution or local motions induced by star formation.
Kinematical evidences for two expanding bubbles, Mrk 86–B and Mrk 86–C, are given. They show expanding velocities of 34 and 17 km s⁻¹, H α luminosities of 3×10³⁸ and 1.7×10³⁹ erg s⁻¹, and physical radii of 374 and 120 pc, respectively. The change in the [S  ii ]/H α , [N  ii ]/H α , [O  ii ]/[O  iii ] and [O  iii ]/H β line ratios with the distance to the bubble precursor suggests a diminution in the ionization parameter and, in the case of Mrk 86–B, an enhancement of the shock-excited gas emission. The optical–near-infrared colours of the bubble precursors are characteristic of low‐metallicity star‐forming regions (∼0.2 Z_⊙) with burst strengths of about 1 per cent in mass.     

The circumstellar material around SN IIn 1997eg: another detection of a very narrow P Cygni profile

We report the detection of a very narrow P Cygni profile on top of the broad emission H α and H β lines of the Type IIn Supernova 1997eg. A similar feature has been detected in SN 1997ab, SN 1998S and SN 1995G . The detection of the narrow P Cygni profile indicates the existence of a dense circumstellar material (CSM), into which the ejecta of the supernova is expanding. From the analysis of the spectra of SN 1997eg we deduce (i) that such CSM is very dense  ( n ≳5×10⁷ cm^-3)  , (ii) that it has a low expanding velocity of about 160 km s⁻¹. The origin of such dense CSM can be either a very dense progenitor wind  ( M˙ ∼10^-2 M_⊙ yr^-1)  or a circumstellar shell product of the progenitor wind expanding into a high-pressure environment.     

Faint Lyman-break galaxies as a crucial test for galaxy formation models

It has recently been shown that galaxy formation models within the Λ cold dark matter cosmology predict that, compared to the observed population, small galaxies (with stellar masses  <10¹¹ M_⊙  ) form too early, are too passive since   z ∼ 3  and host too old stellar populations at   z = 0  . We then expect an overproduction of small galaxies at   z ≳ 4  that should be visible as an excess of faint Lyman-break galaxies. To check whether this excess is present, we use the morgana galaxy formation model and grasil spectrophotometric  +  radiative transfer code to generate mock catalogues of deep fields observed with Hubble Space Telescope Advanced Camera for Surveys. We add observational noise and the effect of Lyman α emission, and perform colour–colour selections to identify Lyman-break galaxies. The resulting mock candidates have plausible properties that closely resemble those of observed galaxies. We are able to reproduce the evolution of the bright tail of the luminosity function of Lyman-break galaxies (with a possible underestimate of the number of the brightest i -dropouts), but uncertainties and degeneracies in dust absorption parameters do not allow to give strong constraints to the model. Besides, our model shows a clear excess with respect to observations of faint Lyman-break galaxies, especially of   z ₈₅₀∼ 27 V   -dropouts at   z ∼ 5  . We quantify the properties of these 'excess' galaxies and discuss the implications: these galaxies are hosted in dark matter haloes with circular velocities in excess of 100 km s⁻¹, and their suppression may require a deep rethinking of stellar feedback processes taking place in galaxy formation.     

Analytical approach to subhalo population in dark matter haloes

In the standard model of cosmic structure formation, dark matter haloes form by gravitational instability. The process is hierarchical: smaller systems collapse earlier, and later merge to form larger haloes. The galaxy clusters, hosted by the largest dark matter haloes, are at the top of this hierarchy and representing the largest as well as the last structures formed in the Universe, while the smaller and first haloes are those Earth-sized dark subhaloes that have been both predicted by theoretical considerations and found in numerical simulations, though there do not exist any observational hints of their existence. The probability that a halo of mass m at redshift z will be part of a larger halo of mass M at the present time can be described in the frame of the extended Press & Schecter theory making use of the progenitor (conditional) mass function. Using the progenitor mass function, we calculate analytically, at redshift zero, the distribution of subhaloes in mass, formation epoch and rarity of the peak of the density field at the formation epoch. That is done for a Milky Way size system, assuming both a spherical and an ellipsoidal collapse model. Our calculation assumes that small progenitors do not lose mass due to dynamical processes after entering the parent halo, and that they do not interact with other subhaloes. For a Λ cold dark matter power spectrum, we obtain a subhalo mass function  d n /d m   proportional to   m ^−α  with a model-independent  α∼ 2  . Assuming that the dark matter is a weakly interacting massive particle, the inferred distributions are used to test the feasibility of an indirect detection in the γ-ray energy band of such a population of subhaloes with a Gamma-ray Large Area Space Telescope like satellite.