The Aquila comparison project: the effects of feedback and numerical methods on simulations of galaxy formation
Cecilia ScannapiecoMarkus WadepuhlOwen H. ParryJulio F. NavarroAdrian JenkinsVolker SpringelRomain TeyssierE. M. CarlsonH. M. P. CouchmanRobert A. CrainClaudio Dalla VecchiaCarlos S. FrenkChiaki KobayashiPierluigi MonacoG. MuranteTakashi OkamotoThomas QuinnJoop SchayeG. S. StinsonTom TheunsJames WadsleySimon D. M. WhiteRory M. Woods
415
Citation
168
Reference
10
Related Paper
Citation Trend
Abstract:
We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge.Keywords:
Smoothed Particle Hydrodynamics
Structure formation
Abstract Nowadays, the scientific community has generally accepted the existence of dark matter. Astrophysicists investigate this assumption by astronomical observation. The galaxy rotation curve is one of the pieces of evidence for dark matter’s existence. The TLM18 telescope detects a 21 cm radio wave from the Milky Way to measure the rotation curve of the Milky Way. This article analyzes the frequency-power data to find out the rotation curve of the galaxy. The rotation curve is compared to existing theories, including the Kepler model, the exponential disk model, and the isothermal dark matter halo. Combining the exponential disk model with the isothermal dark matter halo, the calculation fits well with observation.
Cite
Citations (2)
I provide a model rotation curve for the Milky Way that matches the details of the terminal velocity curve normalized to the Galactocentric distance $R_0 = 8.122$ kpc obtained by the GRAVITY collaboration and the corresponding circular speed of the LSR $\Theta_0 = 233.3$ km/s. The model provides a numerical representation of the azimuthally averaged radial run of the gravitational potential of each mass component of the Galaxy (bulge-bar, stellar disk, gas disk, and dark matter) as represented by the rotation curve of each. It provides precise estimates of quantities like the stellar mass of the Galaxy ($6.16 \pm 0.31 \times 10^{10}\;\mathrm{M}_{\odot}$) and the local density of dark matter ($\rho_{DM}(R_0) = 6.76^{+0.08}_{-0.14} \times 10^{-3}\; \mathrm{M}_{\odot}\,\mathrm{pc}^{-3} = 0.257^{+0.003}_{-0.005}\; \mathrm{GeV}\,\mathrm{cm}^{-3}$). The dark matter density implied by the radial force is less than that found in many studies of the vertical force, perhaps indicating that the usual assumption of a spherical dark matter halo is no longer adequate.
Gravitational potential
Bar (unit)
Cite
Citations (1)
We predict the internal structure and dynamics of present-day disc galaxies using galaxy evolution models within a hierarchical formation scenario. The halo mass aggregation histories, for a flat cold dark matter model with cosmological constant, were generated and used to calculate the virialization of dark matter haloes. A diversity of halo density profiles were obtained, the most typical one being close to that suggested by Navarro, Frenk & White. We modelled the way in which discs in centrifugal equilibrium are built within the evolving dark haloes, using gas accretion rates proportional to the halo mass aggregation rates, and assuming detailed angular momentum conservation. We calculated the gravitational interactions between halo and disc — including the adiabatic contraction of the halo due to disc formation — and the hydrodynamics, star formation and evolution of the galaxy discs. We find that the slope and zero-point of the Tully-Fisher (TF) relation in the infrared bands may be explained as a direct consequence of the cosmological initial conditions. This relation is almost independent of the assumed disc mass fraction, when the disc component in the rotation curve decomposition is non-negligible. Thus, the power spectrum of fluctuations can be normalized at galaxy scales through the TF relation independently of the disc mass fraction assumed. The rms scatter of the model TF relation originates mainly from the scatter in the dark halo structure and, to a minor extension, from the dispersion of the primordial spin parameter λ. The scatter obtained from our models does not disagree with the observational estimates. Our models allow us to understand why the residuals of the TF relation do not correlate significantly with disc size or surface brightness. We can also explain why low and high surface brightness galaxies have the same TF relation; the key point is the dependence of the star formation efficiency on the disc surface density. The correlations between gas fraction and surface brightness, and between scalelength and Vmax obtained with our models agree with those observed. Discs formed within the growing haloes, where λ is assumed to be time independent, have nearly exponential surface density distributions. The shape of the rotation curves changes with disc surface brightness and is nearly flat for most cases. The rotation curve decompositions show a dominance of dark matter down to very small radii, in conflict with some observational inferences. The introduction of shallow cores in the dark halo attenuates this difficulty and produces haloes with slightly smaller rotation velocities. Other features of our galaxy models are not strongly influenced by the shallow core.
Structure formation
Cite
Citations (86)
Strict criteria are applied to the sample of spiral galaxies with measured rotation curves in order to select those objects for which the observed rotation curve is an accurate tracer of the radial force law. The resulting sub-sample of 10 galaxies is then considered in view of two suggested explanations for the discrepancy between the luminous mass and the conventional dynamical mass of galaxies: dark haloes and the modified Newtonian dynamics (MOND). This is done by means of least-squares fits to the rotation curves. Three-parameter dark-halo models (M/L for the visible disc, the core radius and the asymptotic circular velocity of the halo) work well in reproducing the observed rotation curves, and it is found that, for the higher luminosity galaxies, the visible matter dominates the mass distribution within the optically bright disc. However, in the low-luminosity gas-rich dwarfs the dark component is everywhere dominant. MOND, with one free parameter, (M/L for the visible disc) generally works well in predicting the form of the rotation curves, in some cases better than multi-parameter dark-halo fits. If the distance to the galaxy is also taken as a free parameter, then the MOND fits are as good as three parameter dark-halo models and, with one exception, the implied distances are consistent with the adopted distances within the probable uncertainty in the distance estimates. Restricting the number of parameter in dark-halo models by making use of the disc-halo coupling does not produce satisfactory fits to the rotation curves. The overall conclusion is that MOND is currently the best phenomenological description of the systematics of the discrepancy in galaxies.
Modified Newtonian dynamics
Cite
Citations (787)
Recent measurements of gas velocity in the outer parts of high redshift galaxies suggest that steeply falling rotation curves may be common, or even universal, in these galaxies, in contrast to the near universal flat, non-declining rotation curves in nearby galaxies. We investigate the implications of these postulated steeply falling rotation curves for the role of dark matter in galaxy formation. Using an established computer code, the collapse of dark matter and baryonic matter together, starting with a variety of initial conditions, is simulated for comparison with the observed rotation curves. As soon as a smooth stellar disc is formed in the baryonic matter, with properties similar to the observed high redshift galaxies, the computed rotation curves are, without exception, relatively flat to large radius in the gas disc. Only a simulation without a dark matter halo is able to reproduce the observed rotation curves. This would imply that, if the high redshift steeply falling rotation curves turn out to be common, then the standard scenario for galaxy formation for these galaxies, namely baryonic matter falling into the potential well of a massive dark matter halo, must be wrong, unless there is pressure support via velocity dispersion significantly higher than has so far been observed. It would also imply that for these galaxies the flat rotation curves at low redshift must be due to dark matter which has subsequently fallen into the galactic potential well, or there must be some other explanation for the contemporary flat rotation curves, other than dark matter.
Cuspy halo problem
Cold dark matter
Cite
Citations (1)
We show that the Tully-Fisher relation observed for spiral galaxies can be explained in the current scenario of galaxy formation without invoking subtle assumptions, provided that galactic-sized dark haloes have shallow, core-like central profiles, with a core radius proportional to halo circular velocity. In such a system, both disk and halo contribute significantly to the maximum rotation of the disk, and the interaction between the disk and halo components acts to reduce the scatter in the Tully-Fisher relation. With model parameters chosen in plausible ranges, the model can well accommodate the zero-point, slope, scatter of the observed Tully-Fisher relation, as well as the large range of disk sizes. The model predicts that LSB disks obey a Tully-Fisher relation similar to that of normal disks, if disk mass-to-light ratio is properly taken into account. The halo profile required by the Tully-Fisher relation is as shallow as that required by the rotation curves of faint disks, but much shallower than that predicted by conventional CDM models. Our results cannot be explained by some of the recent proposals for resolving the conflict between conventional CDM models and the rotation-curve shapes of faint galaxies. If dark matter self-interaction (either scattering or annihilation) is responsible for the shallow profile, the observed Tully-Fisher relation requires the interaction cross section \sigma_X to satisfy <\sigma_{X}|v|>/m_{X}~10^{-16} cm^3/s/GeV.
Tully–Fisher relation
Cite
Citations (58)
In this work, to investigate the dynamical nature of the kinematic asymmetry in the isolated gas-rich dwarf irregular galaxy W LM in the Local Group, we consider that the dark matter halo and the disk do not have the same center of mass (i.e., the disk lies off-center in the potential of the extended dark matter halo), which is one of the possible physical explanations for the kinematic lopsidedness. To do so, we generate a lopsided halo potential by considering two dark matter mass density models, ISO and Burkert, and we add up the contribution to the rotation curve of a perturbation term [Formula: see text] in the gravitational potential, which arises from the offset between the disk and the dark matter halo. We show that such an m = 1 perturbation improves the rotation curve modeling when compared to a non-perturbed potential and the shape of the HI gas rotation curves is fitted better in the approaching side if the perturbation term in the halo potential is taken into account for this galaxy dynamics. In fact, displacing the disk center by 0.1 kpc from the halo center is sufficient to provide such an improvement in modeling the rotation curve.
Gravitational potential
Cite
Citations (0)
We present the 21-cm rotation curve of the nearby galaxy M33 out to a galactocentric distance of 16 kpc (13 disc scalelengths). The rotation curve keeps rising out to the last measured point and implies a dark halo mass ≳5×1010 M⊙. The stellar and gaseous discs provide virtually equal contributions to the galaxy gravitational potential at large galactocentric radii, but no obvious correlation is found between the radial distribution of dark matter and the distribution of stars or gas. Results of the best fit to the mass distribution in M33 picture a dark halo which controls the gravitational potential from 3 kpc outward, with a matter density which decreases radially as R−1.3. The density profile is consistent with the theoretical predictions for structure formation in hierarchical clustering cold dark matter (CDM) models, and favours lower mass concentrations than those expected in the standard cosmogony.
Gravitational potential
Mass distribution
Cite
Citations (205)
I provide a model rotation curve for the Milky Way that matches the details of the terminal velocity curve normalized to the Galactocentric distance $R_0 = 8.122$ kpc obtained by the GRAVITY collaboration and the corresponding circular speed of the LSR $\Theta_0 = 233.3$ km/s. The model provides a numerical representation of the azimuthally averaged radial run of the gravitational potential of each mass component of the Galaxy (bulge-bar, stellar disk, gas disk, and dark matter) as represented by the rotation curve of each. It provides precise estimates of quantities like the stellar mass of the Galaxy ($6.16 \pm 0.31 \times 10^{10}\;\mathrm{M}_{\odot}$) and the local density of dark matter ($\rho_{DM}(R_0) = 6.76^{+0.08}_{-0.14} \times 10^{-3}\; \mathrm{M}_{\odot}\,\mathrm{pc}^{-3} = 0.257^{+0.003}_{-0.005}\; \mathrm{GeV}\,\mathrm{cm}^{-3}$). The dark matter density implied by the radial force is less than that found in many studies of the vertical force, perhaps indicating that the usual assumption of a spherical dark matter halo is no longer adequate.
Gravitational potential
Bar (unit)
Cite
Citations (2)
In this paper, the effect of halo substructures on galaxy rotation curves is investigated using a simple model of dark matter clustering. A dark matter halo density profile is developed based only on the scale-free nature of clustering that leads to a statistically self-similar distribution of the substructures at the galactic scale. A semi-analytical method is used to derive rotation curves for such a clumpy dark matter density profile. It is found that the halo substructures significantly affect the galaxy velocity field. Based on the fractal geometry of the halo, this self-consistent model predicts a Navarro–Frenk–White-like rotation curve and a scale-free power spectrum of the rotation velocity fluctuations.
Cuspy halo problem
Cite
Citations (1)