Adiabatic Contraction Research Articles

On the universality of density profiles

We use the secondary infall model described in Del Popolo (2009), which takes into account the effect of dynamical friction, ordered and random angular momentum, baryons adiabatic contraction and dark matter baryons interplay, to study how in- ner slopes of relaxed LCDM dark matter (DM) halos with and without baryons (baryons+DM, and pure DM) depend on redshift and on halo mass. We apply the quoted method to structures on galactic scales and clusters of galaxies scales. We find that the inner logarithmic density slope, of dark matter halos with baryons has a significant dependence on halo mass and redshift with slopes ranging from 0 for dwarf galaxies to 0.4 for objects of M = 10^13M_solar and 0.94 for M = 10^15M_solar clusters of galaxies. Structures slopes increase with increasing redshift and this trend reduces going from galaxies to clusters. In the case of density profiles constituted just of dark matter the mass and redshift dependence of slope is very slight. In this last case, we used the Merrit et al. (2006) analysis who compared N-body density profiles with various parametric models finding systematic variation in profile shape with halo mass. This last analysis suggests that the galaxy-sized halos obtained with our model have a different shape parameter, i.e. a different mass distribution, than the cluster-sized halos, obtained with the same model. The results of the present paper argue against universality of density profiles constituted by dark matter and baryons and confirm claims of a systematic variation in profile shape with halo mass, for dark matter halos.

Testing adiabatic contraction with Sloan Digital Sky Survey elliptical galaxies

We study the profiles of 75 086 elliptical galaxies from the Sloan Digital Sky Survey (SDSS) at both large (50-500 kpc/h) and small (~3 kpc/h) scales. Weak lensing observations in the outskirts of the halo are combined with measurements of the stellar velocity dispersion in the interior regions of the galaxy for stacked galaxy samples. The weak lensing measurements are well characterized by a Navarro, Frenk and White (NFW) profile. The dynamical mass measurements exceed the extrapolated NFW profile even after the estimated stellar masses are subtracted, providing evidence for the modification of the dark matter profile by the baryons. This excess mass is quantitatively consistent with the predictions of the adiabatic contraction (AC) hypothesis. Our finding suggests that the effects of AC during galaxy formation are stable to subsequent bombardment from major and minor mergers. We explore several theoretical and observational systematics and conclude that they cannot account for the inferred mass excess. The most significant source of systematic error is in the IMF, which would have to increase the stellar mass estimates by a factor of two relative to masses from the Kroupa IMF to fully explain the mass excess without AC. Such an increase would create tension with results from SAURON (Cappellari et al. 2006). We demonstrate a connection between the level of contraction of the dark matter halo profile and scatter in the size-luminosity relation, which is a projection of the fundamental plane. Whether or not AC is the mechanism supplying the excess mass, models of galaxy formation and evolution must reconcile the observed halo masses from weak lensing with the comparatively large dynamical masses at the half light radii of the galaxies.

Galaxy-induced transformation of dark matter haloes

We use N-body/gasdynamical LambdaCDM cosmological simulations to examine the effect of the assembly of a central galaxy on the shape and mass profile of its dark halo. Two series of simulations are compared; one that follows only the evolution of the dark matter component and a second one where a baryonic component is added. These simulations include radiative cooling but neglect star formation and feedback, leading most baryons to collect at the halo center in a disk which is too small and too massive when compared with typical spiral. This unrealistic model allows us, nevertheless, to gauge the maximum effect that galaxies may have in transforming their dark halos. We find that the shape of the halo becomes more axisymmetric: halos are transformed from triaxial into essentially oblate systems, with well-aligned isopotential contours of roughly constant flattening (c/a ~ 0.85). Halos always contract as a result of galaxy assembly, but the effect is substantially less pronounced than predicted by the "adiabatic contraction" hypothesis. The reduced contraction helps to reconcile LambdaCDM halos with constraints on the dark matter content inside the solar circle and should alleviate the long-standing difficulty of matching simultaneously the scaling properties of galaxy disks and the luminosity function. The halo contraction is also less pronounced than found in earlier simulations, a disagreement that suggests that halo contraction is not solely a function of the initial and final distribution of baryons. Not only how much baryonic mass has been deposited at the center of a halo matters, but also the mode of its deposition. It might prove impossible to predict the halo response without a detailed understanding of a galaxy's assembly history. (Abriged)

Testing adiabatic contraction of dark matter in fossil group candidates

We present deep XMM observations and ESO WFI optical imaging of two X-ray-selected fossil group candidates, RXCJ0216.7-4749 and RXCJ2315.7-0222. Using the X-ray data, we derive total mass profiles under the hydrostatic equilibrium assumption. The central regions of RXCJ0216.7-4749 are found to be dominated by an X-ray bright AGN, and although we derive a mass profile, uncertainties are large and the constraints are significantly weakened due to the presence of the central source. The total mass profile of RXCJ2315.7-0222 is of high quality, being measured in fifteen bins from [0.075 - 0.75]R500 and containing three data points interior to 30 kpc, allowing comprehensive investigation of its properties. We probe several mass models based on the standard NFW profile or on the Sersic-like model recently suggested by high-resolution N-body simulations. We find that the addition of a stellar component due to the presence of the central galaxy is necessary for a good analytical model fit. In all mass profile models fitted, the mass concentration is not especially high compared to non-fossil systems. In addition, the modification of the dark matter halo by adiabatic contraction slightly improves the fit. However, our result depends critically on the choice of IMF used to convert galaxy luminosity to mass, which leads to a degeneracy between the central slope of the dark matter profile and the normalisation of the stellar component. While we argue on the basis of the range of M_*/L_R ratios that lower M_*/L_R ratios are preferred on physical grounds and that adiabatic contraction has thus operated in this system, better theoretical and observational convergence on this problem is needed to make further progess.

Dark matter response to galaxy formation

We have resimulated the six galaxy-sized haloes of the Aquarius Project including metal-dependent cooling, star formation and supernova feedback. This allows us to study not only how dark matter haloes respond to galaxy formation, but also how this response is affected by details of halo assembly history. In agreement with previous work, we find baryon condensation to lead to increased dark matter concentration. Dark matter density profiles differ substantially in shape from halo to halo when baryons are included, but in all cases the velocity dispersion decreases monotonically with radius. Some haloes show an approximately constant dark matter velocity anisotropy with $ \beta \approx 0.1-02$, while others retain the anisotropy structure of their baryon-free versions. Most of our haloes become approximately oblate in their inner regions, although a few retain the shape of their dissipationless counterparts. Pseudo-phase-space densities are described by a power law in radius of altered slope when baryons are included. The shape and concentration of the dark matter density profiles are not well reproduced by published adiabatic contraction models. The significant spread we find in the density and kinematic structure of our haloes appears related to differences in their formation histories. Such differences already affect the final structure in baryon-free simulations, but they are reinforced by the inclusion of baryons, and new features are produced. The details of galaxy formation need to be better understood before the inner dark matter structure of galaxies can be used to constrain cosmological models or the nature of dark matter.

The central dark matter content of early-type galaxies: scaling relations and connections with star formation histories

We examine correlations between the masses, sizes, and star formation histories for a large sample of low-redshift early-type galaxies, using a simple suite of dynamical and stellar populations models. We confirm an anti-correlation between size and stellar age, and survey for trends with the central content of dark matter (DM). An average relation between central DM density and galaxy size of <rho_DM> ~ Reff^-2 provides the first clear indication of cuspy DM haloes in these galaxies -- akin to standard LCDM haloes that have undergone adiabatic contraction. The DM density scales with galaxy mass as expected, deviating from suggestions of a universal halo profile for dwarf and late-type galaxies. We introduce a new fundamental constraint on galaxy formation by finding that the central DM fraction decreases with stellar age. This result is only partially explained by the size-age dependencies, and the residual trend is in the opposite direction to basic DM halo expectations. Therefore we suggest that there may be a connection between age and halo contraction, and that galaxies forming earlier had stronger baryonic feedback which expanded their haloes, or else lumpier baryonic accretion that avoided halo contraction. An alternative explanation is a lighter initial mass function for older stellar populations.

DISSIPATIONAL VERSUS DISSIPATIONLESS GALAXY FORMATION AND THE DARK MATTER CONTENT OF GALAXIES

We examine two extreme models for the build-up of the stellar component of luminous elliptical galaxies. In one case, we assume the build-up of stars is dissipational, with centrally accreted gas radiating away its orbital and thermal energy; the dark matter halo will undergo adiabatic contraction and the central dark matter density profile will steepen. For the second model, we assume the central galaxy is assembled by a series of dissipationless mergers of stellar clumps that have formed far from the nascent galaxy. In order to be accreted, these clumps lose their orbital energy to the dark matter halo via dynamical friction, thereby heating the central dark matter and smoothing the dark matter density cusp. The central dark matter density profiles differ drastically between these models. For the isolated elliptical galaxy, NGC 4494, the central dark matter densities follow the power-laws r^(-0.2) and r^(-1.7) for the dissipational and dissipationless models, respectively. By matching the dissipational and dissipationless models to observations of the stellar component of elliptical galaxies, we examine the relative contributions of dissipational and dissipationless mergers to the formation of elliptical galaxies and look for observational tests that will distinguish between these models. Comparisons to strong lensing brightest cluster galaxies yield median M*/L_B ratios of 2.1+/-0.8 and 5.2+/-1.7 at z=0.39 for the dissipational and dissipationless models, respectively. For NGC 4494, the best-fit dissipational and dissipationless models have M*/L_B=2.97 and 3.96. Comparisons to expected stellar mass-to-light ratios from passive evolution and population syntheses appear to rule out a purely dissipational formation mechanism for the central stellar regions of giant elliptical galaxies.

MAGIC GAMMA-RAY TELESCOPE OBSERVATION OF THE PERSEUS CLUSTER OF GALAXIES: IMPLICATIONS FOR COSMIC RAYS, DARK MATTER, AND NGC 1275

The Perseus galaxy cluster was observed by the MAGIC Cherenkov telescope for a total effective time of 24.4 hr during 2008 November and December. The resulting upper limits on the gamma-ray emission above 100 GeV are in the range of 4.6 to 7.5 x 10^{-12} cm^{-2} s^{-1} for spectral indices from -1.5 to -2.5, thereby constraining the emission produced by cosmic rays, dark matter annihilations, and the central radio galaxy NGC1275. Results are compatible with cosmological cluster simulations for the cosmic-ray-induced gamma-ray emission, constraining the average cosmic ray-to-thermal pressure to <4% for the cluster core region (<8% for the entire cluster). Using simplified assumptions adopted in earlier work (a power-law spectrum with an index of -2.1, constant cosmic ray-to-thermal pressure for the peripheral cluster regions while accounting for the adiabatic contraction during the cooling flow formation), we would limit the ratio of cosmic ray-to-thermal energy to E_CR/E_th<3%. The upper limit also translates into a level of gamma-ray emission from possible annihilations of the cluster dark matter (the dominant mass component) that is consistent with boost factors of ~10^4 for the typically expected dark matter annihilation-induced emission. Finally, the upper limits obtained for the gamma-ray emission of the central radio galaxy NGC1275 are consistent with the recent detection by the Fermi-LAT satellite. Due to the extremely large Doppler factors required for the jet, a one-zone synchrotron self-Compton model is implausible in this case. We reproduce the observed spectral energy density by using the structured jet (spine-layer) model which has previously been adopted to explain the high-energy emission of radio galaxies.

Effects of dark matter substructures on gravitational lensing: results from the Aquarius simulations

We use high-resolution Aquarius simulations of Milky Way-sized haloes in the LCDM cosmology to study the effects of dark matter substructures on gravitational lensing. Each halo is resolved with ~ 10^8 particles (at a mass resolution ~ 10^3-4 M_sun/h) within its virial radius. Subhaloes with masses larger than 10^5 M_sun/h are well resolved, an improvement of at least two orders of magnitude over previous lensing studies. We incorporate a baryonic component modelled as a Hernquist profile and account for the response of the dark matter via adiabatic contraction. We focus on the "anomalous" flux ratio problem, in particular on the violation of the cusp-caustic relation due to substructures. We find that subhaloes with masses less than ~ 10^8 M_sun/h play an important role in causing flux anomalies; such low mass subhaloes have been unresolved in previous studies. There is large scatter in the predicted flux ratios between different haloes and between different projections of the same halo. In some cases, the frequency of predicted anomalous flux ratios is comparable to that observed for the radio lenses, although in most cases it is not. The probability for the simulations to reproduce the observed violations of the cusp lenses is about 0.001. We therefore conclude that the amount of substructure in the central regions of the Aquarius haloes is insufficient to explain the observed frequency of violations of the cusp-caustic relation. These conclusions are based purely on our dark matter simulations which ignore the effect of baryons on subhalo survivability.

Density profiles of dark matter haloes on galactic and cluster scales

In the present paper, we improve the "Extended Secondary Infall Model" (ESIM) of Williams et al. (2004) to obtain further insights on the cusp/core problem. The model takes into account the effect of ordered and random angular momentum, dynamical friction and baryon adiabatic contraction in order to obtain a secondary infall model more close to the collapse reality. The model is applied to structures on galactic scales (normal and dwarf spiral galaxies) and on cluster of galaxies scales. The results obtained suggest that angular momentum and dynamical friction are able, on galactic scales, to overcome the competing effect of adiabatic contraction eliminating the cusp. The NFW profile can be reobtained, in our model only if the system is constituted just by dark matter and the magnitude of angular momentum and dynamical friction are reduced with respect to the values predicted by the model itself. The rotation curves of four LSB galaxies from de Blok & Bosma (2002) are compared to the rotation curves obtained by the model in the present paper obtaining a good fit to the observational data. On scales smaller than $\simeq 10^{11} h^{-1} M_{\odot}$ the slope $\alpha \simeq 0$ and on cluster scales we observe a similar evolution of the dark matter density profile but in this case the density profile slope flattens to $\alpha \simeq 0.6$ for a cluster of $\simeq 10^{14} h^{-1} M_{\odot}$. The total mass profile, differently from that of dark matter, shows a central cusp well fitted by a NFW model.

Galaxy density profiles and shapes - II. Selection biases in strong lensing surveys

[Abridged] Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio, and image separation distribution. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sersic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population.

THE CUSP/CORE PROBLEM AND THE SECONDARY INFALL MODEL

We study the cusp/core problem using a secondary infall model (SIM) that takes into account the effect of ordered and random angular momentum, dynamical friction and baryons adiabatic contraction. The model is applied to structures on galactic scales (normal and dwarfs spiral galaxies) and on clusters of galaxies scales. Our analysis suggest that angular momentum and dynamical friction are able, on galactic scales, to overcome the competing effect of adiabatic contraction eliminating the cusp. The slope of density profile of inner haloes flattens with decreasing halo mass and the profile is well approximated by a Burkert's profile. In order to obtain the NFW profile, starting from the profiles obtained from our model, the magnitude of angular momentum and dynamical friction must be reduced with respect to the values predicted by the model itself. The rotation curves of four LSB galaxies from Gentile et al. (2004) are compared to the rotation curves obtained by the model in the present paper obtaining a good fit to the observational data. The time evolution of the density profile of a galaxy of $10^8-10^9 M_{\odot}$ shows that after a transient steepening, due to the adiabatic contraction, the density profile flattens to $\alpha \simeq 0$. On cluster scales we observe a similar evolution of the dark matter density profile but in this case the density profile slope flattens to $\alpha \simeq 0.6$ for a cluster of $\simeq 10^{14} M_{\odot}$. The total mass profile, differently from that of dark matter, shows a central cusp well fitted by a NFW model.

The impact of feedback on disc galaxy scaling relations

We use a disc galaxy evolution model to investigate the impact of mass outflows (a.k.a. feedback) on disc galaxy scaling relations, mass fractions and spin parameters. Our model follows the accretion, cooling, star formation and ejection of baryonic mass inside growing dark matter haloes, with cosmologically motivated angular momentum distributions, and dark matter halo structure. Models without feedback produce discs that are too small, too gas-poor and which rotate too fast. Feedback reduces the galaxy formation efficiency ∈ GF (defined as the fraction of the universally available baryons that end up as stars and cold gas in a given galaxy), resulting in larger discs with higher gas fractions and lower rotation velocities. Models with feedback can reproduce the zero-points of the scaling relations among rotation velocity, stellar mass and disc size, but only in the absence of adiabatic contraction. Our feedback mechanism is maximally efficient in expelling mass, but our successful models require 25 per cent of the supernova (SN) energy or 100 per cent of the SN momentum to drive an outflow. It remains to be seen whether such high efficiencies are realistic or not. Our energy- and momentum-driven wind models result in different slopes of various scaling relations. Energy-driven winds result in steeper slopes to the galaxy-mass―halo-mass and stellar-mass-halo-mass relations, a shallower slope to the galaxy-size-stellar-mass relation at z = 0 and a steeper slope to the cold gas metallicity-stellar-mass relation at z ≃ 2. Observations favour the energy-driven wind at stellar masses below M star ≤ 10 10.5 M ⊙ , but the momentum-driven wind model at high masses. The ratio between the specific angular momentum of the baryons to that of the halo (j gal /m gal ) is not unity in our models with inflow and outflow. Yet this is the standard assumption in models of disc formation. Above a halo mass of M vir ≃ 10 12 M ⊙ , cooling becomes increasingly inefficient, which results in (j gal /m gal ) decreasing with increasing halo mass. Below a halo mass of M vir ≃ 10 12 M ⊙ , feedback becomes increasingly efficient. Feedback preferentially ejects low-angular-momentum material because star formation is more efficient at smaller galactic radii and at higher redshifts. This results in (j gal /m gal ) increasing with decreasing halo mass. This effect helps to resolve the discrepancy between the high spin parameters observed for dwarf galaxies with the low spin parameters predicted from A cold dark matter.

DARK MATTER DENSITIES DURING THE FORMATION OF THE FIRST STARS AND IN DARK STARS

The first stars in the universe form inside ∼106 M☉ dark matter (DM) halos whose initial density profiles are laid down by gravitational collapse in hierarchical structure-formation scenarios. During the formation of the first stars in the universe, the baryonic infall compresses the DM further. The resultant DM density is presented here, using an algorithm originally developed by Young to calculate changes to the profile as a result of adiabatic infall in a spherical halo model; the Young prescription takes into account the noncircular motions of halo particles. The density profiles obtained in this way are found to be within a factor of 2 of those obtained using the simple adiabatic contraction prescription of Blumenthal and colleagues. Our results hold regardless of the nature of the DM or its interactions and rely merely on gravity. If the DM consists of weakly interacting massive particles, which are their own antiparticles, their densities are high enough that their annihilation in the first protostars can indeed provide an important heat source and prevent the collapse all the way to fusion. In short, a "Dark Star" phase of stellar evolution, powered by DM annihilation, may indeed describe the first stars in the universe.

DARK MATTER ANNIHILATION AND PRIMORDIAL STAR FORMATION

We investigate the effects of weakly-interacting massive particle (WIMP) dark matter annihilation on the formation of Population III.1 stars, which are theorized to form from the collapse of gas cores at the centers of dark matter minihalos. We consider the relative importance of cooling due to baryonic radiative processes and heating due to WIMP annihilation. We analyze the dark matter and gas profiles of several halos formed in cosmological-scale numerical simulations. The heating rate depends sensitively on the dark matter density profile, which we approximate with a power law rho_chi ~ r^{-alpha_chi}, in the numerically unresolved inner regions of the halo. If we assume a self-similar structure so that alpha_chi ~= 1.5 as measured on the resolved scales ~1pc, then for a fiducial WIMP mass of 100GeV, the heating rate is typically much smaller (<10^{-3}) than the cooling rate for densities up to n_H=10^{17}cm^{-3}. In one case, where alpha_chi=1.65, the heating rate becomes similar to the cooling rate by a density of n_H=10^{15}cm^{-3}. The dark matter density profile is expected to steepen in the central baryon-dominated region <~1pc due to adiabatic contraction, and we observe this effect (though with relatively low resolution) in our numerical models. From these we estimate alpha_chi~=2.0. The heating now dominates cooling above n_H~=10^{14}cm^{-3}, in agreement with the previous study of Spolyar, Freese & Gondolo. We expect this leads to the formation of an equilibrium structure with a baryonic and dark matter density distribution exhibiting a flattened central core. Examining such equilibria, we find total luminosities due to WIMP annihilation are relatively constant and ~10^3 L_sun, set by the radiative luminosity of the baryonic core. We discuss the implications for Pop III.1 star formation... (abridged)

Effects of galaxy-halo alignment and adiabatic contraction on gravitational lens statistics

We study the strong gravitational lens statistics of triaxial cold dark matter (CDM) halos occupied by central early-type galaxies. We calculate the image separation distribution for double, cusp and quad configurations. The ratios of image multiplicities at large separations are consistent with the triaxial NFW model, and at small separations are consistent with the singular isothermal ellipsoid (SIE) model. At all separations, the total lensing probability is enhanced by adiabatic contraction. If no adiabatic contraction is assumed, naked cusp configurations become dominant at approximately 2.5'', which is inconsistent with the data. We also show that at small-to-moderate separations, the image multiplicities depend sensitively on the alignment of the shapes of the luminous and dark matter projected density profiles. In constrast to other properties that affect these ratios, the degree of alignment does not have a significant effect on the total lensing probability. These correlations may therefore be constrained by comparing the theoretical image separation distribution to a sufficiently large lens sample from future wide and deep sky surveys such as Pan-Starrs, LSST and JDEM. Understanding the correlations in the shapes of galaxies and their dark matter halo is important for future weak lensing surveys.

Dark matter annihilation effects on the first stars

We study the effects of WIMP dark matter (DM) on the collapse and evolution of the first stars in the Universe. Using a stellar evolution code, we follow the pre-Main Sequence (MS) phase of a grid of metal-free stars with masses in the range 5-600 solar mass forming in the centre of a 1e6 solar mass halo at redhisft z=20. DM particles of the parent halo are accreted in the proto-stellar interior by adiabatic contraction and scattering/capture processes, reaching central densities of order 1e12 GeV/cm3 at radii of the order of 10 AU. Energy release from annihilation reactions can effectively counteract the gravitational collapse, in agreement with results from other groups. We find this stalling phase (known as "dark" star) is transients and lasts from 2.1e3 yr (M=600 solar mass) to 1.8e4 yr (M=9 solar mass). Later in the evolution, DM scattering/capture rate becomes high enough that energy deposition from annihilations significantly alters the pre-MS evolution of the star in a way that depends on DM (i) velocity dispersion, (ii) density, (iii) elastic scattering cross section with baryons. For our fiducial set of parameters (10 km/s, 1e11 GeV/cm3, 1e-38 cm2) we find that the evolution of stars of mass lower than 40 solar masses "freezes" on the HR diagram before reaching the ZAMS. Stars with bigger masses manage to ignite nuclear reactions; however, DM "burning" prolonges their lifetimes by a factor 2 (5) for a 600 (40) solar mass star.

A Cosmologically Motivated Description of the Dark Matter Halo Profile for the Low Surface Brightness Galaxy, Malin 1

ABSTRACTIn this paper we derive a possible mass profile for the low surface brightness galaxy, Malin 1, based upon previously published space-based and ground-based photometric properties and kinematics. We use properties of the bulge, normal disk, outer extended disk, and H i mass as inputs into mass profile models. We find that the dark matter halo model of Malin 1 is best described by a halo profile that has undergone adiabatic contraction, inconsistent with the findings for most disk galaxies to date, yet consistent with rotation curve studies of M31. More importantly, we find that Malin 1 is baryon dominated in its central regions out to a radius of (in the bulge region). Low surface brightness galaxies are often referred to as being dark matter dominated at all radii. If this is the case, then Malin 1 would seem to have characteristics similar to those of normal barred disk galaxies, as suggested by other recent work. We also find that Malin 1 falls on the rotation curve shear versus spiral arm pitch angle relation for normal galaxies, although more LSB galaxies need to be studied to determine if this is typical.

Formation of Disk Galaxies in Computer Simulations

The formation of disk galaxies is one of the most outstanding problems in modern astrophysics and cosmology. We review the progress made by numerical simulations carried out on large parallel supercomputers. These simulations model the formation of disk galaxies within the current structure formation paradigm in which the Universe is dominated by a cold dark matter component and a cosmological constant. We discuss how computer simulations have been an essential tool in advancing the field further over the last decade or so. Recent progress stems from a combination of increased resolution and improved treatment of the astrophysical processes modeled in the simulations, such as the phenomenological description of the interstellar medium and of the process of star formation. We argue that high mass and spatial resolution is a necessary condition in order to obtain large disks comparable with observed spiral galaxies avoiding spurious dissipation of angular momentum. A realistic model of the star formation history. gas-to-stars ratio and the morphology of the stellar and gaseous component is instead controlled by the phenomenological description of the non-gravitational energy budget in the galaxy. This includes the energy injection by supernovae explosions as well as by accreting supermassive black holes at scales below the resolution. We continue by showing that simulations of gas collapse within cold dark matter halos including a phenomenological description of supernovae blast-waves allow to obtain stellar disks with nearly exponential surface density profiles as those observed in real disk galaxies, counteracting the tendency of gas collapsing in such halos to form cuspy baryonic profiles. However, the ab-initio formation of a realistic rotationally supported disk galaxy with a pure exponential disk in a fully cosmological simulation is still an open problem. We argue that the suppression of bulge formation is related to the physics of galaxy formation during the merger of the most massive protogalactic lumps at high redshift, where the reionization of the Universe likely plays a key role. A sufficiently high resolution during this early phase of galaxy formation is also crucial to avoid artificial angular momentum loss and spurious bulge formation. Finally, we discuss the role of mergers in disk formation, adiabatic halo contraction during the assembly of the disk, cold flows, thermal instability and other aspects of galaxy formation, focusing on their relevance to the puzzling origin of bulgeless galaxies.

The impact of the dark matter-gas interaction on the collapse behavior of spherical symmetric systems

Context. If the gas in the evolving cosmic halos is dissipating energy (cooling), then the dark matter (DM) halo also undergoes a compactification because of the variation in the gravitational potential . This is well-known as adiabatic contraction (AC).