Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userModeling Galaxy-Mass Correlations in Dissipationless Simulations
48
0
0.0
(
0
)
Authors
Argyro Tasitsiomi
Ask ChatGPT about the research
No Arabic abstract
We use high-resolution dissipationless simulations of the concordance flat LCDM model to make predictions for the galaxy--mass (GM) correlations and compare them to the recent SDSS weak lensing measurements.We use a simple observationally motivated scheme to assign luminosities and colors to the halos.This allows us to closely match the selection criteria used to define observational samples.The simulations reproduce the observed GM correlation function and its observed dependencies on luminosity and color.The luminosity dependence of the correlation function is primarily determined by the changing relative contribution of central and satellite galaxies at different luminosities. The color dependence of the GM correlations reflects the difference in the typical environments of blue and red galaxies. We also find agreement between the predicted and observed cross-bias, b_x=b/r,at all probed scales.The GM correlation coefficient, r, is close to unity on scales >1/h Mpc.The cross bias is thus expected to measure the actual bias of galaxy clustering on these scales.The aperture mass-to-light ratio is independent of galaxy color.The aperture mass scales approximately linearly with luminosity at L_r>10^{10}h^{-2} Lsun, while at lower luminosities the scaling is shallower: L_r^{0.5}. We show that most of the luminous galaxies (M_r<-21) are near the centers of their halos and their GM correlation function at r<100/h kpc can therefore be interpreted as the average dark matter density profile of these galaxies. We find that for galaxies in a given narrow luminosity range, there is a broad and possibly non-gaussian distribution of halo virial masses. Therefore, the average relation between mass and luminosity derived from the weak lensing analyses should be interpreted with caution.
rate research
Read More
Galaxy-galaxy lensing is a powerful probe of the relation between galaxies and dark matter halos, but its theoretical interpretation requires a careful modeling of various contributions, such as the contribution from central and satellite galaxies. For this purpose, a phenomenological approach based on the halo model has been developed, allowing for fast exploration of the parameter space of models. In this paper, we investigate the ability of the halo model to extract information from the g-g weak lensing signal by comparing it to high-resolution dissipationless simulations that resolve subhalos. We find that the halo model reliably determines parameters such as the host halo mass of central galaxies, the fraction of galaxies that are satellites, and their radial distribution inside larger halos. If there is a significant scatter present in the central galaxy host halo mass distribution, then the mean and median mass of that distribution can differ significantly from one another, and the halo model mass determination lies between the two. This result suggests that when analyzing the data, galaxy subsamples with a narrow central galaxy halo mass distribution, such as those based on stellar mass, should be chosen for a simpler interpretation of the results.
We present a method of including galaxy formation in dissipationless N-body simulations. Galaxies that form during the evolution are identified at several epochs and replaced by single massive soft particles. This allows one to produce two-component models containing galaxies and a background dark matter distribution. We applied this technique to obtain two sets of models: one for field galaxies and one for galaxy clusters. We tested the method for the standard CDM scenario for structure formation in the universe. A direct comparison of the simulated galaxy distribution to the observed one sets the amplitude of the initial density fluctuation spectrum, and thus the present time in the simulations. The rates of formation and merging compare very well to simulations that include hydrodynamics, and are compatible with observations. We also discuss the cluster luminosity function.
We exploit a suite of large N-body simulations (up to N=$4096^3$) performed with Abacus, of scale-free models with a range of spectral indices $n$, to better understand and quantify convergence of the matter power spectrum in dark matter only cosmological N-body simulations. Using self-similarity to identify converged regions, we show that the maximal wavenumber resolved at a given level of accuracy increases monotonically as a function of time. At the $1%$ level it starts at early times from a fraction of $k_Lambda$, the Nyquist wavenumber of the initial grid, and reaches at most, if the force softening is sufficiently small, $sim 2 k_Lambda$ at the very latest times we evolve to. At the $5%$ level accuracy extends up to slightly larger wavenumbers, of order $5k_Lambda$ at late times. Expressed as a suitable function of the scale-factor, accuracy shows a very simple $n$-dependence, allowing a straightforward extrapolation to place conservative bounds on the accuracy of N-body simulations of non-scale free models like LCDM. Quantitatively our findings are broadly in line with the conservative assumptions about resolution adopted by recent studies using large cosmological simulations (e.g. Euclid Flagship) aiming to constrain the mildly non-linear regime. On the other hand, we note that studies of the matter power spectrum in the literature have often used data at larger wavenumbers, where convergence to the physical result is poor. Even qualitative conclusions about clustering at small scales, e.g concerning the validity of the stable clustering approximation, may need revision in light of our results.
Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, energy conserving Sedov-Taylor, and momentum conserving snowplow$-$ with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, $f_{mathrm{hot}}$, as $(1 - f_{mathrm{hot}})^{-4/5}$. We also include winds from young massive stars and AGB stars, Stromgren sphere gas heating by massive stars, and a gas cooling limiting mechanism driven by radiative recombination of dense HII regions. We present initial tests for isolated Milky-Way like systems simulated with the GADGET based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a two order of magnitude variation in particle mass in the range (1.3$-$130)$times 10^4$ solar masses.
Dissipationless collapses in Modified Newtonian Dynamics (MOND) are studied by using a new particle-mesh N-body code based on our numerical MOND potential solver. We found that low surface-density end-products have shallower inner density profile, flatter radial velocity-dispersion profile, and more radially anisotropic orbital distribution than high surface-density end-products. The projected density profiles of the final virialized systems are well described by Sersic profiles with index m~4, down to m~2 for a deep-MOND collapse. Consistently with observations of elliptical galaxies, the MOND end-products, if interpreted in the context of Newtonian gravity, would appear to have little or no dark matter within the effective radius. However, we found impossible (under the assumption of constant mass-to-light ratio) to simultaneously place the resulting systems on the observed Kormendy, Faber-Jackson and Fundamental Plane relations of elliptical galaxies. Finally, the simulations provide strong evidence that phase mixing is less effective in MOND than in Newtonian gravity.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
