Normalized to: Jenkins, A.
[1]
oai:arXiv.org:2006.16249 [pdf] - 2124970
Primordial black holes from cusp collapse on cosmic strings
Submitted: 2020-06-29
Primordial black holes (PBHs) are of fundamental interest in cosmology and
astrophysics, and have received much attention as a dark matter candidate and
as a potential source of gravitational waves. One possible PBH formation
mechanism is the gravitational collapse of cosmic strings. Thus far, the
entirety of the literature on PBH production from cosmic strings has focused on
the collapse of (quasi)circular cosmic string loops, which make up only a tiny
fraction of the cosmic loop population. We demonstrate here a novel PBH
formation mechanism: the collapse of a small segment of cosmic string in the
neighbourhood of a cusp. Using the hoop conjecture, we show that collapse is
inevitable whenever a cusp appears on a macroscopically-large loop, forming a
PBH whose mass is smaller than that of the loop by a factor of the
dimensionless string tension $G\mu$. Since cusps are generic features of cosmic
string loops, and do not rely on finely-tuned loop configurations like circular
collapse, this implies that cosmic strings produce PBHs in far greater numbers
than has previously been recognised. The resulting PBHs are highly spinning and
boosted to ultrarelativistic velocities; they populate a unique region of the
BH mass-spin parameter space, and are therefore a ``smoking gun'' observational
signature of cosmic strings. We derive stringent new constraints on $G\mu$ from
the evaporation of cusp-collapse PBHs, and update existing constraints on
$G\mu$ from gravitational-wave searches.
[3]
oai:arXiv.org:2005.14495 [pdf] - 2104881
A high-resolution cosmological simulation of a strong gravitational lens
Submitted: 2020-05-29, last modified: 2020-06-01
We present a cosmological hydrodynamical simulation of a 10^13 Msun galaxy
group and its environment (out to 10 times the virial radius) carried out using
the EAGLE model of galaxy formation. Exploiting a novel technique to increase
the resolution of the dark matter calculation independently of that of the gas,
the simulation resolves dark matter haloes and subhaloes of mass 5x10^6 Msun .
It is therefore useful for studying the abundance and properties of the haloes
and subhaloes targeted in strong lensing tests of the cold dark matter model.
We estimate the halo and subhalo mass functions and discuss how they are
affected both by the inclusion of baryons in the simulation and by the
environment. We find that the halo and subhalo mass functions have lower
amplitude in the hydrodynamical simulation than in its dark matter only
counterpart. This reflects the reduced growth of haloes in the hydrodynamical
simulation due to the early loss of gas by reionisation and galactic winds and,
additionally, in the case of subhaloes, disruption by enhanced tidal effects
within the host halo due to the presence of a massive central galaxy. The
distribution of haloes is highly anisotropic reflecting the filamentary
character of mass accretion onto the cluster. As a result, there is significant
variation in the number of structures with viewing direction. The median number
of structures near the centre of the halo, when viewed in projection, is
reduced by a factor of two when baryons are included.
[4]
oai:arXiv.org:2005.12926 [pdf] - 2102978
The origin of X-ray coronae around simulated disc galaxies
Submitted: 2020-05-26
The existence of hot, accreted gaseous coronae around massive galaxies is a
long-standing central prediction of galaxy formation models in the $\Lambda$CDM
cosmology. While observations now confirm that extraplanar hot gas is present
around late-type galaxies, the origin of the gas is uncertain with suggestions
that galactic feedback could be the dominant source of energy powering the
emission. We investigate the origin and X-ray properties of the hot gas that
surrounds galaxies of halo mass, $(10^{11}-10^{14}) \mathrm{M}_\odot$, in the
cosmological hydrodynamical EAGLE simulations. We find that the central X-ray
emission, $\leq 0.10 R_{\mathrm{vir}}$, of halos of mass $\leq 10^{13}
\mathrm{M}_\odot$ originates from gas heated by supernovae (SNe). However,
beyond this region, a quasi-hydrostatic, accreted atmosphere dominates the
X-ray emission in halos of mass $\geq 10^{12} \mathrm{M}_\odot$. We predict
that a dependence on halo mass of the hot gas to dark matter mass fraction can
significantly change the slope of the $L_{\mathrm{X}}-M_{\mathrm{vir}}$
relation (which is typically assumed to be $4/3$ for clusters) and we derive
the scaling law appropriate to this case. As the gas fraction in halos
increases with halo mass, we find a steeper slope for the
$L_{\mathrm{X}}-M_{\mathrm{vir}}$ in lower mass halos, $\leq 10^{14}
\mathrm{M}_\odot$. This varying gas fraction is driven by active galactic
nuclei (AGN) feedback. We also identify the physical origin of the so-called
"missing feedback" problem, the apparently low X-ray luminosities observed from
high star-forming, low-mass galaxies. This is explained by the ejection of
SNe-heated gas from the central regions of the halo.
[6]
oai:arXiv.org:1909.11627 [pdf] - 2097205
Projection effects on the observed angular spectrum of the astrophysical
stochastic gravitational wave background
Submitted: 2019-09-25, last modified: 2020-04-23
The detection and characterization of the Stochastic Gravitational Wave
Background (SGWB) is one of the main goals of Gravitational Wave (GW)
experiments. The observed SGWB will be the combination of GWs from cosmological
(as predicted by many models describing the physics of the early Universe) and
astrophysical origins, which will arise from the superposition of GWs from
unresolved sources whose signal is too faint to be detected. Therefore, it is
important to have a proper modeling of the astrophysical SGWB (ASGWB) in order
to disentangle the two signals; moreover, this will provide additional
information on astrophysical properties of compact objects. Applying the Cosmic
Rulers formalism, we compute the observed ASGWB angular power spectrum, hence
using gauge invariant quantities, accounting for all effects intervening
between the source and the observer. These are the so-called projection
effects, which include Kaiser, Doppler and gravitational potentials effect. Our
results show that these projection effects are the most important at the
largest scales, and they contribute to up to tens of percent of the angular
power spectrum amplitude, with the Kaiser term being the largest at all scales.
While the exact impact of these results will depend on instrumental and
astrophysical details, a precise theoretical modeling of the ASGWB will
necessarily need to include all these projection effects.
[8]
oai:arXiv.org:2003.04321 [pdf] - 2061816
Setting the Stage: Structures from Gaussian Random Fields
Submitted: 2020-03-09
We study structure formation in a set of cosmological simulations to uncover
the scales in the initial density field that gave rise to the formation of
present-day structures. Our simulations share a common primordial power
spectrum (here Lambda-CDM), but the introduction of hierarchical variations of
the phase information allows us to systematically study the scales that
determine the formation of structure at later times. We consider the variance
in z=0 statistics such as the matter power spectrum and halo mass function. We
also define a criterion for the existence of individual haloes across
simulations, and determine what scales in the initial density field contain
sufficient information for the non-linear formation of unique haloes. We study
how the characteristics of individual haloes such as the mass and
concentration, as well as the position and velocity, are affected by variations
on different scales, and give scaling relations for haloes of different mass.
Finally, we use the example of a cluster-mass halo to show how our hierarchical
parametrisation of the initial density field can be used to create variants of
particular objects. With properties such as mass, concentration, kinematics and
substructure of haloes set on distinct and well-determined scales, and its
unique ability to introduce variations localised in real space, our method is a
powerful tool to study structure formation in cosmological simulations.
[12]
oai:arXiv.org:1902.07719 [pdf] - 2038203
Shot noise in the astrophysical gravitational-wave background
Submitted: 2019-02-20, last modified: 2020-01-27
We calculate the noise induced in the anisotropies of the astrophysical
gravitational-wave background by finite sampling of both the galaxy
distribution and the compact binary coalescence event rate. This shot noise
leads to a scale-invariant bias term in the angular power spectrum $C_\ell$,
for which we derive a simple analytical expression. We find that this bias
dominates over the true cosmological power spectrum in any reasonable observing
scenario, and that only with very long observing times and removal of a large
number of foreground sources can the true power spectrum be recovered.
[13]
oai:arXiv.org:1907.06642 [pdf] - 2038234
Estimating the angular power spectrum of the gravitational-wave
background in the presence of shot noise
Submitted: 2019-07-15, last modified: 2020-01-27
There has been much recent interest in studying anisotropies in the
astrophysical gravitational-wave (GW) background, as these could provide us
with interesting new information about galaxy clustering and large-scale
structure. However, this information is obscured by shot noise, caused by the
finite number of GW sources that contribute to the background at any given
time. We develop a new method for estimating the angular spectrum of
anisotropies, based on the principle of combining statistically-independent
data segments. We show that this gives an unbiased estimate of the true,
astrophysical spectrum, removing the offset due to shot noise power, and that
in the limit of many data segments, it is the most efficient (i.e.
lowest-variance) estimator possible.
[19]
oai:arXiv.org:1911.09720 [pdf] - 2002513
Universality in the structure of dark matter haloes over twenty orders
of magnitude in halo mass
Submitted: 2019-11-21
Dark matter haloes are the basic units of all cosmic structure. They grew by
gravitational amplification of weak initial density fluctuations that are still
visible on large scales in the cosmic microwave background radiation. Galaxies
formed within relatively massive haloes as gas cooled and condensed at their
centres, but many hypotheses for the nature of dark matter imply that the halo
population should extend to masses many orders of magnitude below those where
galaxies can form. Here, we use a novel, multi-zoom technique to create the
first consistent simulation of the formation of present-day haloes over the
full mass range populated when dark matter is aWeakly Interacting Massive
Particle (WIMP) of mass ~100 GeV. The simulation has a dynamic range of 30
orders of magnitude in mass, resolving the internal structure of hundreds of
Earth-mass haloes just as well as that of hundreds of rich galaxy clusters.
Remarkably, halo density profiles are universal over the entire mass range and
are well described by simple two-parameter fitting formulae. Halo mass and
concentration are tightly related in a way which depends on cosmology and on
the nature of the dark matter. At fixed mass, concentration is independent of
local environment for haloes less massive than those of typical galaxies. These
results are important for predicting annihilation radiation signals from dark
matter, since these should be dominated by contributions from the smallest
structures.
[20]
oai:arXiv.org:1811.12907 [pdf] - 1973647
GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers
Observed by LIGO and Virgo during the First and Second Observing Runs
The LIGO Scientific Collaboration;
the Virgo Collaboration;
Abbott, B. P.;
Abbott, R.;
Abbott, T. D.;
Abraham, S.;
Acernese, F.;
Ackley, K.;
Adams, C.;
Adhikari, R. X.;
Adya, V. B.;
Affeldt, C.;
Agathos, M.;
Agatsuma, K.;
Aggarwal, N.;
Aguiar, O. D.;
Aiello, L.;
Ain, A.;
Ajith, P.;
Allen, G.;
Allocca, A.;
Aloy, M. A.;
Altin, P. A.;
Amato, A.;
Ananyeva, A.;
Anderson, S. B.;
Anderson, W. G.;
Angelova, S. V.;
Antier, S.;
Appert, S.;
Arai, K.;
Araya, M. C.;
Areeda, J. S.;
Arène, M.;
Arnaud, N.;
Arun, K. G.;
Ascenzi, S.;
Ashton, G.;
Aston, S. M.;
Astone, P.;
Aubin, F.;
Aufmuth, P.;
AultONeal, K.;
Austin, C.;
Avendano, V.;
Avila-Alvarez, A.;
Babak, S.;
Bacon, P.;
Badaracco, F.;
Bader, M. K. M.;
Bae, S.;
Baker, P. T.;
Baldaccini, F.;
Ballardin, G.;
Ballmer, S. W.;
Banagiri, S.;
Barayoga, J. C.;
Barclay, S. E.;
Barish, B. C.;
Barker, D.;
Barkett, K.;
Barnum, S.;
Barone, F.;
Barr, B.;
Barsotti, L.;
Barsuglia, M.;
Barta, D.;
Bartlett, J.;
Bartos, I.;
Bassiri, R.;
Basti, A.;
Bawaj, M.;
Bayley, J. C.;
Bazzan, M.;
Bécsy, B.;
Bejger, M.;
Belahcene, I.;
Bell, A. S.;
Beniwal, D.;
Berger, B. K.;
Bergmann, G.;
Bernuzzi, S.;
Bero, J. J.;
Berry, C. P. L.;
Bersanetti, D.;
Bertolini, A.;
Betzwieser, J.;
Bhandare, R.;
Bidler, J.;
Bilenko, I. A.;
Bilgili, S. A.;
Billingsley, G.;
Birch, J.;
Birney, R.;
Birnholtz, O.;
Biscans, S.;
Biscoveanu, S.;
Bisht, A.;
Bitossi, M.;
Bizouard, M. A.;
Blackburn, J. K.;
Blackman, J.;
Blair, C. D.;
Blair, D. G.;
Blair, R. M.;
Bloemen, S.;
Bode, N.;
Boer, M.;
Boetzel, Y.;
Bogaert, G.;
Bondu, F.;
Bonilla, E.;
Bonnand, R.;
Booker, P.;
Boom, B. A.;
Booth, C. D.;
Bork, R.;
Boschi, V.;
Bose, S.;
Bossie, K.;
Bossilkov, V.;
Bosveld, J.;
Bouffanais, Y.;
Bozzi, A.;
Bradaschia, C.;
Brady, P. R.;
Bramley, A.;
Branchesi, M.;
Brau, J. E.;
Briant, T.;
Briggs, J. H.;
Brighenti, F.;
Brillet, A.;
Brinkmann, M.;
Brisson, V.;
Brockill, P.;
Brooks, A. F.;
Brown, D. D.;
Brunett, S.;
Buikema, A.;
Bulik, T.;
Bulten, H. J.;
Buonanno, A.;
Buskulic, D.;
Rosell, M. J. Bustamante;
Buy, C.;
Byer, R. L.;
Cabero, M.;
Cadonati, L.;
Cagnoli, G.;
Cahillane, C.;
Bustillo, J. Calderón;
Callister, T. A.;
Calloni, E.;
Camp, J. B.;
Campbell, W. A.;
Canepa, M.;
Cannon, K. C.;
Cao, H.;
Cao, J.;
Capocasa, E.;
Carbognani, F.;
Caride, S.;
Carney, M. F.;
Carullo, G.;
Diaz, J. Casanueva;
Casentini, C.;
Caudill, S.;
Cavaglià, M.;
Cavalier, F.;
Cavalieri, R.;
Cella, G.;
Cerdá-Durán, P.;
Cerretani, G.;
Cesarini, E.;
Chaibi, O.;
Chakravarti, K.;
Chamberlin, S. J.;
Chan, M.;
Chao, S.;
Charlton, P.;
Chase, E. A.;
Chassande-Mottin, E.;
Chatterjee, D.;
Chaturvedi, M.;
Chatziioannou, K.;
Cheeseboro, B. D.;
Chen, H. Y.;
Chen, X.;
Chen, Y.;
Cheng, H. -P.;
Cheong, C. K.;
Chia, H. Y.;
Chincarini, A.;
Chiummo, A.;
Cho, G.;
Cho, H. S.;
Cho, M.;
Christensen, N.;
Chu, Q.;
Chua, S.;
Chung, K. W.;
Chung, S.;
Ciani, G.;
Ciobanu, A. A.;
Ciolfi, R.;
Cipriano, F.;
Cirone, A.;
Clara, F.;
Clark, J. A.;
Clearwater, P.;
Cleva, F.;
Cocchieri, C.;
Coccia, E.;
Cohadon, P. -F.;
Cohen, D.;
Colgan, R.;
Colleoni, M.;
Collette, C. G.;
Collins, C.;
Cominsky, L. R.;
Constancio, M.;
Conti, L.;
Cooper, S. J.;
Corban, P.;
Corbitt, T. R.;
Cordero-Carrión, I.;
Corley, K. R.;
Cornish, N.;
Corsi, A.;
Cortese, S.;
Costa, C. A.;
Cotesta, R.;
Coughlin, M. W.;
Coughlin, S. B.;
Coulon, J. -P.;
Countryman, S. T.;
Couvares, P.;
Covas, P. B.;
Cowan, E. E.;
Coward, D. M.;
Cowart, M. J.;
Coyne, D. C.;
Coyne, R.;
Creighton, J. D. E.;
Creighton, T. D.;
Cripe, J.;
Croquette, M.;
Crowder, S. G.;
Cullen, T. J.;
Cumming, A.;
Cunningham, L.;
Cuoco, E.;
Canton, T. Dal;
Dálya, G.;
Danilishin, S. L.;
D'Antonio, S.;
Danzmann, K.;
Dasgupta, A.;
Costa, C. F. Da Silva;
Datrier, L. E. H.;
Dattilo, V.;
Dave, I.;
Davier, M.;
Davis, D.;
Daw, E. J.;
DeBra, D.;
Deenadayalan, M.;
Degallaix, J.;
De Laurentis, M.;
Deléglise, S.;
Del Pozzo, W.;
DeMarchi, L. M.;
Demos, N.;
Dent, T.;
De Pietri, R.;
Derby, J.;
De Rosa, R.;
De Rossi, C.;
DeSalvo, R.;
de Varona, O.;
Dhurandhar, S.;
Díaz, M. C.;
Dietrich, T.;
Di Fiore, L.;
Di Giovanni, M.;
Di Girolamo, T.;
Di Lieto, A.;
Ding, B.;
Di Pace, S.;
Di Palma, I.;
Di Renzo, F.;
Dmitriev, A.;
Doctor, Z.;
Donovan, F.;
Dooley, K. L.;
Doravari, S.;
Dorrington, I.;
Downes, T. P.;
Drago, M.;
Driggers, J. C.;
Du, Z.;
Ducoin, J. -G.;
Dupej, P.;
Dwyer, S. E.;
Easter, P. J.;
Edo, T. B.;
Edwards, M. C.;
Effler, A.;
Ehrens, P.;
Eichholz, J.;
Eikenberry, S. S.;
Eisenmann, M.;
Eisenstein, R. A.;
Essick, R. C.;
Estelles, H.;
Estevez, D.;
Etienne, Z. B.;
Etzel, T.;
Evans, M.;
Evans, T. M.;
Fafone, V.;
Fair, H.;
Fairhurst, S.;
Fan, X.;
Farinon, S.;
Farr, B.;
Farr, W. M.;
Fauchon-Jones, E. J.;
Favata, M.;
Fays, M.;
Fazio, M.;
Fee, C.;
Feicht, J.;
Fejer, M. M.;
Feng, F.;
Fernandez-Galiana, A.;
Ferrante, I.;
Ferreira, E. C.;
Ferreira, T. A.;
Ferrini, F.;
Fidecaro, F.;
Fiori, I.;
Fiorucci, D.;
Fishbach, M.;
Fisher, R. P.;
Fishner, J. M.;
Fitz-Axen, M.;
Flaminio, R.;
Fletcher, M.;
Flynn, E.;
Fong, H.;
Font, J. A.;
Forsyth, P. W. F.;
Fournier, J. -D.;
Frasca, S.;
Frasconi, F.;
Frei, Z.;
Freise, A.;
Frey, R.;
Frey, V.;
Fritschel, P.;
Frolov, V. V.;
Fulda, P.;
Fyffe, M.;
Gabbard, H. A.;
Gadre, B. U.;
Gaebel, S. M.;
Gair, J. R.;
Gammaitoni, L.;
Ganija, M. R.;
Gaonkar, S. G.;
Garcia, A.;
García-Quirós, C.;
Garufi, F.;
Gateley, B.;
Gaudio, S.;
Gaur, G.;
Gayathri, V.;
Gemme, G.;
Genin, E.;
Gennai, A.;
George, D.;
George, J.;
Gergely, L.;
Germain, V.;
Ghonge, S.;
Ghosh, Abhirup;
Ghosh, Archisman;
Ghosh, S.;
Giacomazzo, B.;
Giaime, J. A.;
Giardina, K. D.;
Giazotto, A.;
Gill, K.;
Giordano, G.;
Glover, L.;
Godwin, P.;
Goetz, E.;
Goetz, R.;
Goncharov, B.;
González, G.;
Castro, J. M. Gonzalez;
Gopakumar, A.;
Gorodetsky, M. L.;
Gossan, S. E.;
Gosselin, M.;
Gouaty, R.;
Grado, A.;
Graef, C.;
Granata, M.;
Grant, A.;
Gras, S.;
Grassia, P.;
Gray, C.;
Gray, R.;
Greco, G.;
Green, A. C.;
Green, R.;
Gretarsson, E. M.;
Groot, P.;
Grote, H.;
Grunewald, S.;
Gruning, P.;
Guidi, G. M.;
Gulati, H. K.;
Guo, Y.;
Gupta, A.;
Gupta, M. K.;
Gustafson, E. K.;
Gustafson, R.;
Haegel, L.;
Halim, O.;
Hall, B. R.;
Hall, E. D.;
Hamilton, E. Z.;
Hammond, G.;
Haney, M.;
Hanke, M. M.;
Hanks, J.;
Hanna, C.;
Hannam, M. D.;
Hannuksela, O. A.;
Hanson, J.;
Hardwick, T.;
Haris, K.;
Harms, J.;
Harry, G. M.;
Harry, I. W.;
Haster, C. -J.;
Haughian, K.;
Hayes, F. J.;
Healy, J.;
Heidmann, A.;
Heintze, M. C.;
Heitmann, H.;
Hello, P.;
Hemming, G.;
Hendry, M.;
Heng, I. S.;
Hennig, J.;
Heptonstall, A. W.;
Vivanco, Francisco Hernandez;
Heurs, M.;
Hild, S.;
Hinderer, T.;
Hoak, D.;
Hochheim, S.;
Hofman, D.;
Holgado, A. M.;
Holland, N. A.;
Holt, K.;
Holz, D. E.;
Hopkins, P.;
Horst, C.;
Hough, J.;
Howell, E. J.;
Hoy, C. G.;
Hreibi, A.;
Huang, Y.;
Huerta, E. A.;
Huet, D.;
Hughey, B.;
Hulko, M.;
Husa, S.;
Huttner, S. H.;
Huynh-Dinh, T.;
Idzkowski, B.;
Iess, A.;
Ingram, C.;
Inta, R.;
Intini, G.;
Irwin, B.;
Isa, H. N.;
Isac, J. -M.;
Isi, M.;
Iyer, B. R.;
Izumi, K.;
Jacqmin, T.;
Jadhav, S. J.;
Jani, K.;
Janthalur, N. N.;
Jaranowski, P.;
Jenkins, A. C.;
Jiang, J.;
Johnson, D. S.;
Johnson-McDaniel, N. K.;
Jones, A. W.;
Jones, D. I.;
Jones, R.;
Jonker, R. J. G.;
Ju, L.;
Junker, J.;
Kalaghatgi, C. V.;
Kalogera, V.;
Kamai, B.;
Kandhasamy, S.;
Kang, G.;
Kanner, J. B.;
Kapadia, S. J.;
Karki, S.;
Karvinen, K. S.;
Kashyap, R.;
Kasprzack, M.;
Katsanevas, S.;
Katsavounidis, E.;
Katzman, W.;
Kaufer, S.;
Kawabe, K.;
Keerthana, N. V.;
Kéfélian, F.;
Keitel, D.;
Kennedy, R.;
Key, J. S.;
Khalili, F. Y.;
Khan, H.;
Khan, I.;
Khan, S.;
Khan, Z.;
Khazanov, E. A.;
Khursheed, M.;
Kijbunchoo, N.;
Kim, Chunglee;
Kim, J. C.;
Kim, K.;
Kim, W.;
Kim, W. S.;
Kim, Y. -M.;
Kimball, C.;
King, E. J.;
King, P. J.;
Kinley-Hanlon, M.;
Kirchhoff, R.;
Kissel, J. S.;
Kleybolte, L.;
Klika, J. H.;
Klimenko, S.;
Knowles, T. D.;
Koch, P.;
Koehlenbeck, S. M.;
Koekoek, G.;
Koley, S.;
Kondrashov, V.;
Kontos, A.;
Koper, N.;
Korobko, M.;
Korth, W. Z.;
Kowalska, I.;
Kozak, D. B.;
Kringel, V.;
Krishnendu, N.;
Królak, A.;
Kuehn, G.;
Kumar, A.;
Kumar, P.;
Kumar, R.;
Kumar, S.;
Kuo, L.;
Kutynia, A.;
Kwang, S.;
Lackey, B. D.;
Lai, K. H.;
Lam, T. L.;
Landry, M.;
Lane, B. B.;
Lang, R. N.;
Lange, J.;
Lantz, B.;
Lanza, R. K.;
Lartaux-Vollard, A.;
Lasky, P. D.;
Laxen, M.;
Lazzarini, A.;
Lazzaro, C.;
Leaci, P.;
Leavey, S.;
Lecoeuche, Y. K.;
Lee, C. H.;
Lee, H. K.;
Lee, H. M.;
Lee, H. W.;
Lee, J.;
Lee, K.;
Lehmann, J.;
Lenon, A.;
Leroy, N.;
Letendre, N.;
Levin, Y.;
Li, J.;
Li, K. J. L.;
Li, T. G. F.;
Li, X.;
Lin, F.;
Linde, F.;
Linker, S. D.;
Littenberg, T. B.;
Liu, J.;
Liu, X.;
Lo, R. K. L.;
Lockerbie, N. A.;
London, L. T.;
Longo, A.;
Lorenzini, M.;
Loriette, V.;
Lormand, M.;
Losurdo, G.;
Lough, J. D.;
Lousto, C. O.;
Lovelace, G.;
Lower, M. E.;
Lück, H.;
Lumaca, D.;
Lundgren, A. P.;
Lynch, R.;
Ma, Y.;
Macas, R.;
Macfoy, S.;
MacInnis, M.;
Macleod, D. M.;
Macquet, A.;
Magaña-Sandoval, F.;
Zertuche, L. Magaña;
Magee, R. M.;
Majorana, E.;
Maksimovic, I.;
Malik, A.;
Man, N.;
Mandic, V.;
Mangano, V.;
Mansell, G. L.;
Manske, M.;
Mantovani, M.;
Marchesoni, F.;
Marion, F.;
Márka, S.;
Márka, Z.;
Markakis, C.;
Markosyan, A. S.;
Markowitz, A.;
Maros, E.;
Marquina, A.;
Marsat, S.;
Martelli, F.;
Martin, I. W.;
Martin, R. M.;
Martynov, D. V.;
Mason, K.;
Massera, E.;
Masserot, A.;
Massinger, T. J.;
Masso-Reid, M.;
Mastrogiovanni, S.;
Matas, A.;
Matichard, F.;
Matone, L.;
Mavalvala, N.;
Mazumder, N.;
McCann, J. J.;
McCarthy, R.;
McClelland, D. E.;
McCormick, S.;
McCuller, L.;
McGuire, S. C.;
McIver, J.;
McManus, D. J.;
McRae, T.;
McWilliams, S. T.;
Meacher, D.;
Meadors, G. D.;
Mehmet, M.;
Mehta, A. K.;
Meidam, J.;
Melatos, A.;
Mendell, G.;
Mercer, R. A.;
Mereni, L.;
Merilh, E. L.;
Merzougui, M.;
Meshkov, S.;
Messenger, C.;
Messick, C.;
Metzdorff, R.;
Meyers, P. M.;
Miao, H.;
Michel, C.;
Middleton, H.;
Mikhailov, E. E.;
Milano, L.;
Miller, A. L.;
Miller, A.;
Millhouse, M.;
Mills, J. C.;
Milovich-Goff, M. C.;
Minazzoli, O.;
Minenkov, Y.;
Mishkin, A.;
Mishra, C.;
Mistry, T.;
Mitra, S.;
Mitrofanov, V. P.;
Mitselmakher, G.;
Mittleman, R.;
Mo, G.;
Moffa, D.;
Mogushi, K.;
Mohapatra, S. R. P.;
Montani, M.;
Moore, C. J.;
Moraru, D.;
Moreno, G.;
Morisaki, S.;
Mours, B.;
Mow-Lowry, C. M.;
Mukherjee, Arunava;
Mukherjee, D.;
Mukherjee, S.;
Mukund, N.;
Mullavey, A.;
Munch, J.;
Muñiz, E. A.;
Muratore, M.;
Murray, P. G.;
Nagar, A.;
Nardecchia, I.;
Naticchioni, L.;
Nayak, R. K.;
Neilson, J.;
Nelemans, G.;
Nelson, T. J. N.;
Nery, M.;
Neunzert, A.;
Ng, K. Y.;
Ng, S.;
Nguyen, P.;
Nichols, D.;
Nielsen, A. B.;
Nissanke, S.;
Nitz, A.;
Nocera, F.;
North, C.;
Nuttall, L. K.;
Obergaulinger, M.;
Oberling, J.;
O'Brien, B. D.;
O'Dea, G. D.;
Ogin, G. H.;
Oh, J. J.;
Oh, S. H.;
Ohme, F.;
Ohta, H.;
Okada, M. A.;
Oliver, M.;
Oppermann, P.;
Oram, Richard J.;
O'Reilly, B.;
Ormiston, R. G.;
Ortega, L. F.;
O'Shaughnessy, R.;
Ossokine, S.;
Ottaway, D. J.;
Overmier, H.;
Owen, B. J.;
Pace, A. E.;
Pagano, G.;
Page, M. A.;
Pai, A.;
Pai, S. A.;
Palamos, J. R.;
Palashov, O.;
Palomba, C.;
Pal-Singh, A.;
Pan, Huang-Wei;
Pang, B.;
Pang, P. T. H.;
Pankow, C.;
Pannarale, F.;
Pant, B. C.;
Paoletti, F.;
Paoli, A.;
Papa, M. A.;
Parida, A.;
Parker, W.;
Pascucci, D.;
Pasqualetti, A.;
Passaquieti, R.;
Passuello, D.;
Patil, M.;
Patricelli, B.;
Pearlstone, B. L.;
Pedersen, C.;
Pedraza, M.;
Pedurand, R.;
Pele, A.;
Penn, S.;
Perego, A.;
Perez, C. J.;
Perreca, A.;
Pfeiffer, H. P.;
Phelps, M.;
Phukon, K. S.;
Piccinni, O. J.;
Pichot, M.;
Piergiovanni, F.;
Pillant, G.;
Pinard, L.;
Pirello, M.;
Pitkin, M.;
Poggiani, R.;
Pong, D. Y. T.;
Ponrathnam, S.;
Popolizio, P.;
Porter, E. K.;
Powell, J.;
Prajapati, A. K.;
Prasad, J.;
Prasai, K.;
Prasanna, R.;
Pratten, G.;
Prestegard, T.;
Privitera, S.;
Prodi, G. A.;
Prokhorov, L. G.;
Puncken, O.;
Punturo, M.;
Puppo, P.;
Pürrer, M.;
Qi, H.;
Quetschke, V.;
Quinonez, P. J.;
Quintero, E. A.;
Quitzow-James, R.;
Raab, F. J.;
Radkins, H.;
Radulescu, N.;
Raffai, P.;
Raja, S.;
Rajan, C.;
Rajbhandari, B.;
Rakhmanov, M.;
Ramirez, K. E.;
Ramos-Buades, A.;
Rana, Javed;
Rao, K.;
Rapagnani, P.;
Raymond, V.;
Razzano, M.;
Read, J.;
Regimbau, T.;
Rei, L.;
Reid, S.;
Reitze, D. H.;
Ren, W.;
Ricci, F.;
Richardson, C. J.;
Richardson, J. W.;
Ricker, P. M.;
Riemenschneider, G. M.;
Riles, K.;
Rizzo, M.;
Robertson, N. A.;
Robie, R.;
Robinet, F.;
Rocchi, A.;
Rolland, L.;
Rollins, J. G.;
Roma, V. J.;
Romanelli, M.;
Romano, R.;
Romel, C. L.;
Romie, J. H.;
Rose, K.;
Rosińska, D.;
Rosofsky, S. G.;
Ross, M. P.;
Rowan, S.;
Rüdiger, A.;
Ruggi, P.;
Rutins, G.;
Ryan, K.;
Sachdev, S.;
Sadecki, T.;
Sakellariadou, M.;
Salafia, O.;
Salconi, L.;
Saleem, M.;
Salemi, F.;
Samajdar, A.;
Sammut, L.;
Sanchez, E. J.;
Sanchez, L. E.;
Sanchis-Gual, N.;
Sandberg, V.;
Sanders, J. R.;
Santiago, K. A.;
Sarin, N.;
Sassolas, B.;
Sathyaprakash, B. S.;
Saulson, P. R.;
Sauter, O.;
Savage, R. L.;
Schale, P.;
Scheel, M.;
Scheuer, J.;
Schmidt, P.;
Schnabel, R.;
Schofield, R. M. S.;
Schönbeck, A.;
Schreiber, E.;
Schulte, B. W.;
Schutz, B. F.;
Schwalbe, S. G.;
Scott, J.;
Scott, S. M.;
Seidel, E.;
Sellers, D.;
Sengupta, A. S.;
Sennett, N.;
Sentenac, D.;
Sequino, V.;
Sergeev, A.;
Setyawati, Y.;
Shaddock, D. A.;
Shaffer, T.;
Shahriar, M. S.;
Shaner, M. B.;
Shao, L.;
Sharma, P.;
Shawhan, P.;
Shen, H.;
Shink, R.;
Shoemaker, D. H.;
Shoemaker, D. M.;
ShyamSundar, S.;
Siellez, K.;
Sieniawska, M.;
Sigg, D.;
Silva, A. D.;
Singer, L. P.;
Singh, N.;
Singhal, A.;
Sintes, A. M.;
Sitmukhambetov, S.;
Skliris, V.;
Slagmolen, B. J. J.;
Slaven-Blair, T. J.;
Smith, J. R.;
Smith, R. J. E.;
Somala, S.;
Son, E. J.;
Sorazu, B.;
Sorrentino, F.;
Souradeep, T.;
Sowell, E.;
Spencer, A. P.;
Srivastava, A. K.;
Srivastava, V.;
Staats, K.;
Stachie, C.;
Standke, M.;
Steer, D. A.;
Steinke, M.;
Steinlechner, J.;
Steinlechner, S.;
Steinmeyer, D.;
Stevenson, S. P.;
Stocks, D.;
Stone, R.;
Stops, D. J.;
Strain, K. A.;
Stratta, G.;
Strigin, S. E.;
Strunk, A.;
Sturani, R.;
Stuver, A. L.;
Sudhir, V.;
Summerscales, T. Z.;
Sun, L.;
Sunil, S.;
Suresh, J.;
Sutton, P. J.;
Swinkels, B. L.;
Szczepańczyk, M. J.;
Tacca, M.;
Tait, S. C.;
Talbot, C.;
Talukder, D.;
Tanner, D. B.;
Tápai, M.;
Taracchini, A.;
Tasson, J. D.;
Taylor, R.;
Thies, F.;
Thomas, M.;
Thomas, P.;
Thondapu, S. R.;
Thorne, K. A.;
Thrane, E.;
Tiwari, Shubhanshu;
Tiwari, Srishti;
Tiwari, V.;
Toland, K.;
Tonelli, M.;
Tornasi, Z.;
Torres-Forné, A.;
Torrie, C. I.;
Töyrä, D.;
Travasso, F.;
Traylor, G.;
Tringali, M. C.;
Trovato, A.;
Trozzo, L.;
Trudeau, R.;
Tsang, K. W.;
Tse, M.;
Tso, R.;
Tsukada, L.;
Tsuna, D.;
Tuyenbayev, D.;
Ueno, K.;
Ugolini, D.;
Unnikrishnan, C. S.;
Urban, A. L.;
Usman, S. A.;
Vahlbruch, H.;
Vajente, G.;
Valdes, G.;
van Bakel, N.;
van Beuzekom, M.;
Brand, J. F. J. van den;
Broeck, C. Van Den;
Vander-Hyde, D. C.;
van Heijningen, J. V.;
van der Schaaf, L.;
van Veggel, A. A.;
Vardaro, M.;
Varma, V.;
Vass, S.;
Vasúth, M.;
Vecchio, A.;
Vedovato, G.;
Veitch, J.;
Veitch, P. J.;
Venkateswara, K.;
Venugopalan, G.;
Verkindt, D.;
Vetrano, F.;
Viceré, A.;
Viets, A. D.;
Vine, D. J.;
Vinet, J. -Y.;
Vitale, S.;
Vo, T.;
Vocca, H.;
Vorvick, C.;
Vyatchanin, S. P.;
Wade, A. R.;
Wade, L. E.;
Wade, M.;
Walet, R.;
Walker, M.;
Wallace, L.;
Walsh, S.;
Wang, G.;
Wang, H.;
Wang, J. Z.;
Wang, W. H.;
Wang, Y. F.;
Ward, R. L.;
Warden, Z. A.;
Warner, J.;
Was, M.;
Watchi, J.;
Weaver, B.;
Wei, L. -W.;
Weinert, M.;
Weinstein, A. J.;
Weiss, R.;
Wellmann, F.;
Wen, L.;
Wessel, E. K.;
Weßels, P.;
Westhouse, J. W.;
Wette, K.;
Whelan, J. T.;
White, L. V.;
Whiting, B. F.;
Whittle, C.;
Wilken, D. M.;
Williams, D.;
Williamson, A. R.;
Willis, J. L.;
Willke, B.;
Wimmer, M. H.;
Winkler, W.;
Wipf, C. C.;
Wittel, H.;
Woan, G.;
Woehler, J.;
Wofford, J. K.;
Worden, J.;
Wright, J. L.;
Wu, D. S.;
Wysocki, D. M.;
Xiao, L.;
Yamamoto, H.;
Yancey, C. C.;
Yang, L.;
Yap, M. J.;
Yazback, M.;
Yeeles, D. W.;
Yu, Hang;
Yu, Haocun;
Yuen, S. H. R.;
Yvert, M.;
Zadrożny, A. K.;
Zanolin, M.;
Zappa, F.;
Zelenova, T.;
Zendri, J. -P.;
Zevin, M.;
Zhang, J.;
Zhang, L.;
Zhang, T.;
Zhao, C.;
Zhou, M.;
Zhou, Z.;
Zhu, X. J.;
Zimmerman, A. B.;
Zlochower, Y.;
Zucker, M. E.;
Zweizig, J.
Submitted: 2018-11-30, last modified: 2019-10-04
We present the results from three gravitational-wave searches for coalescing
compact binaries with component masses above 1$\mathrm{M}_\odot$ during the
first and second observing runs of the Advanced gravitational-wave detector
network. During the first observing run (O1), from September $12^\mathrm{th}$,
2015 to January $19^\mathrm{th}$, 2016, gravitational waves from three binary
black hole mergers were detected. The second observing run (O2), which ran from
November $30^\mathrm{th}$, 2016 to August $25^\mathrm{th}$, 2017, saw the first
detection of gravitational waves from a binary neutron star inspiral, in
addition to the observation of gravitational waves from a total of seven binary
black hole mergers, four of which we report here for the first time: GW170729,
GW170809, GW170818 and GW170823. For all significant gravitational-wave events,
we provide estimates of the source properties. The detected binary black holes
have total masses between $18.6_{-0.7}^{+3.2}\mathrm{M}_\odot$, and
$84.4_{-11.1}^{+15.8} \mathrm{M}_\odot$, and range in distance between
$320_{-110}^{+120}$ Mpc and $2840_{-1360}^{+1400}$ Mpc. No neutron star - black
hole mergers were detected. In addition to highly significant
gravitational-wave events, we also provide a list of marginal event candidates
with an estimated false alarm rate less than 1 per 30 days. From these results
over the first two observing runs, which include approximately one
gravitational-wave detection per 15 days of data searched, we infer merger
rates at the 90% confidence intervals of $110\, -\, 3840$
$\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$ for binary neutron stars and $9.7\, -\,
101$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$ for binary black holes assuming fixed
population distributions, and determine a neutron star - black hole merger rate
90% upper limit of $610$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$.
[21]
oai:arXiv.org:1909.00819 [pdf] - 2082792
Probing the gravitational wave background from cosmic strings with LISA
Auclair, Pierre;
Blanco-Pillado, Jose J.;
Figueroa, Daniel G.;
Jenkins, Alexander C.;
Lewicki, Marek;
Sakellariadou, Mairi;
Sanidas, Sotiris;
Sousa, Lara;
Steer, Daniele A.;
Wachter, Jeremy M.;
Kuroyanagi, Sachiko
Submitted: 2019-09-02
Cosmic string networks offer one of the best prospects for detection of
cosmological gravitational waves (GWs). The combined incoherent GW emission of
a large number of string loops leads to a stochastic GW background (SGWB),
which encodes the properties of the string network. In this paper we analyze
the ability of the Laser Interferometer Space Antenna (LISA) to measure this
background, considering leading models of the string networks. We find that
LISA will be able to probe cosmic strings with tensions $G\mu \gtrsim
\mathcal{O}(10^{-17})$, improving by about $6$ orders of magnitude current
pulsar timing arrays (PTA) constraints, and potentially $3$ orders of magnitude
with respect to expected constraints from next generation PTA observatories. We
include in our analysis possible modifications of the SGWB spectrum due to
different hypotheses regarding cosmic history and the underlying physics of the
string network. These include possible modifications in the SGWB spectrum due
to changes in the number of relativistic degrees of freedom in the early
Universe, the presence of a non-standard equation of state before the onset of
radiation domination, or changes to the network dynamics due to a string
inter-commutation probability less than unity. In the event of a detection,
LISA's frequency band is well-positioned to probe such cosmic events. Our
results constitute a thorough exploration of the cosmic string science that
will be accessible to LISA.
[24]
oai:arXiv.org:1811.12437 [pdf] - 2025350
Subhalo destruction in the Apostle and Auriga simulations
Richings, Jack;
Frenk, Carlos;
Jenkins, Adrian;
Robertson, Andrew;
Fattahi, Azadeh;
Grand, Robert J. J.;
Navarro, Julio;
Pakmor, Rudiger;
Gomez, Facundo A.;
Marinacci, Federico
Submitted: 2018-11-29, last modified: 2019-07-05
N-body simulations make unambiguous predictions for the abundance of
substructures within dark matter halos. However, the inclusion of baryons in
the simulations changes the picture because processes associated with the
presence of a large galaxy in the halo can destroy subhalos and substantially
alter the mass function and velocity distribution of subhalos. We compare the
effect of galaxy formation on subhalo populations in two state-of-the-art sets
of hydrodynamical CDM simulations of Milky Way mass halos, APOSTLE and AURIGA.
We introduce a new method for tracking the orbits of subhalos between
simulation snapshots that gives accurate results down to a few kiloparsecs from
the centre of the halo. Relative to a dark matter-only simulation, the
abundance of subhalos in APOSTLE is reduced by 50% near the centre and by 10%
within r200. In AURIGA the corresponding numbers are 80% and 40%. The velocity
distributions of subhalos are also affected by the presence of the galaxy, much
more so in AURIGA than in APOSTLE . The differences on subhalo properties in
the two simulations can be traced back to the mass of the central galaxies,
which in AURIGA are typically twice as massive as those in APOSTLE . We show
that some of the results from previous studies are inaccurate due to systematic
errors in the modelling of subhalo orbits near the centre of halos.
[26]
oai:arXiv.org:1808.08276 [pdf] - 1897843
Galaxy formation in the Planck Millennium: the atomic hydrogen content
of dark matter halos
Baugh, C. M.;
Gonzalez-Perez, Violeta;
Lagos, Claudia D. P.;
Lacey, Cedric G.;
Helly, John;
Jenkins, Adrian;
Frenk, Carlos S.;
Benson, Andrew;
Bower, Richard;
Cole, Shaun
Submitted: 2018-08-24, last modified: 2019-06-11
We present recalibrations of the GALFORM semi-analytical model of galaxy
formation in a new N-body simulation with the Planck cosmology. The Planck
Millennium simulation uses more than 128 billion particles to resolve the
matter distribution in a cube of $800$ Mpc on a side, which contains more than
77 million dark matter haloes with mass greater than $2.12 \times 10^{9} h^{-1}
{\rm M_{\odot}}$ at the present day. Only minor changes to a very small number
of model parameters are required in the recalibration. We present predictions
for the atomic hydrogen content (HI) of dark matter halos, which is a key input
into the calculation of the HI intensity mapping signal expected from the
large-scale structure of the Universe. We find that the HI mass $-$ halo mass
relation displays a clear break at the halo mass above which AGN heating
suppresses gas cooling, $\approx 3 \times 10^{11} h^{-1} M_{\rm \odot}$. Below
this halo mass, the HI content of haloes is dominated by the central galaxy;
above this mass it is the combined HI content of satellites that prevails. We
find that the HI mass - halo mass relation changes little with redshift up to
$z=3$. The bias of HI sources shows a scale dependence that gets more
pronounced with increasing redshift.
[28]
oai:arXiv.org:1810.03635 [pdf] - 1892184
No cores in dark matter-dominated dwarf galaxies with bursty star
formation histories
Bose, Sownak;
Frenk, Carlos S.;
Jenkins, Adrian;
Fattahi, Azadeh;
Gomez, Facundo A.;
Grand, Robert J. J.;
Marinacci, Federico;
Navarro, Julio F.;
Oman, Kyle A.;
Pakmor, Ruediger;
Schaye, Joop;
Simpson, Christine M.;
Springel, Volker
Submitted: 2018-10-08, last modified: 2019-05-30
Measurements of the rotation curves of dwarf galaxies are often interpreted
as requiring a constant density core at the centre, at odds with the "cuspy"
inner profiles predicted by $N$-body simulations of cold dark matter (CDM)
haloes. It has been suggested that this conflict could be resolved by
fluctuations in the inner gravitational potential caused by the periodic
removal of gas following bursts of star formation. Earlier work has suggested
that core formation requires a bursty and extended star formation history
(SFH). Here we investigate the structure of CDM haloes of dwarf galaxies
($M_{{\rm DM}} \sim 10^9-5\times10^{10}\,{\rm M}_\odot$) formed in the APOSTLE
('A Project of Simulating the Local Environment') and AURIGA cosmological
hydrodynamic simulations. Our simulations have comparable or better resolution
than others that make cores ($M_{{\rm gas}} \sim 10^4\,{\rm M}_\odot$,
gravitational softening $\sim 150$ pc). Yet, we do not find evidence of core
formation at {\it any} mass or any correlation between the inner slope of the
DM density profile and temporal variations in the SFH. APOSTLE and AURIGA
dwarfs display a similar diversity in their cumulative SFHs to available data
for Local Group dwarfs. Dwarfs in both simulations are DM-dominated on all
resolved scales at all times, likely limiting the ability of gas outflows to
alter significantly the central density profiles of their haloes. We conclude
that recurrent bursts of star formation are not sufficient to cause the
formation of cores, and that other conditions must also be met for baryons to
be able to modify the central DM cusp.
[33]
oai:arXiv.org:1810.13435 [pdf] - 1854106
Anisotropies in the astrophysical gravitational-wave background: The
impact of black hole distributions
Submitted: 2018-10-31, last modified: 2019-03-22
We use population inference to explore the impact that uncertainties in the
distribution of binary black holes (BBH) have on the astrophysical
gravitational-wave background (AGWB). Our results show that the AGWB monopole
is sensitive to the nature of the BBH population (particularly the local merger
rate), while the anisotropic $C_\ell$ spectrum is only modified to within a few
percent, at a level which is insignificant compared to other sources of
uncertainty (such as cosmic variance). This is very promising news for future
observational studies of the AGWB, as it shows that (i) the monopole can be
used as a new probe of the population of compact objects throughout cosmic
history, complementary to direct observations by LIGO and Virgo and (ii) we are
able to make surprisingly robust predictions for the $C_\ell$ spectrum, even
with only very approximate knowledge of the black hole population. As a result,
the AGWB anisotropies have enormous potential as a new probe of the large-scale
structure of the Universe, and of late-Universe cosmology in general.
[36]
oai:arXiv.org:1901.01078 [pdf] - 1809685
Response to Cusin et al's comment on arXiv:1810.13435
Submitted: 2019-01-04
We offer a brief response to the criticisms put forward by Cusin et al in
arXiv:1811.03582 about our work arXiv:1810.13435 and arXiv:1806.01718,
emphasising that none of these criticisms are relevant to our main results.
[40]
oai:arXiv.org:1809.06275 [pdf] - 1791025
Can we detect quantum gravity with compact binary inspirals?
Submitted: 2018-09-17, last modified: 2018-10-09
Treating general relativity as an effective field theory, we compute the
leading-order quantum corrections to the orbits and gravitational-wave emission
of astrophysical compact binaries. These corrections are independent of the
(unknown) nature of quantum gravity at high energies, and generate a phase
shift and amplitude increase in the observed gravitational-wave signal.
Unfortunately (but unsurprisingly), these corrections are undetectably small,
even in the most optimistic observational scenarios.
[41]
oai:arXiv.org:1809.09625 [pdf] - 1989039
The Milky Way's total satellite population and constraining the mass of
the warm dark matter particle
Submitted: 2018-09-25
The Milky Way's (MW) satellite population is a powerful probe of warm dark
matter (WDM) models as the abundance of small substructures is very sensitive
to the properties of the WDM particle. However, only a partial census of the
MW's complement of satellite galaxies exists because surveys of the MW's close
environs are incomplete both in depth and in sky coverage. We present a new
Bayesian analysis that combines the sample of satellites recently discovered by
the Dark Energy Survey (DES) with those found in the Sloan Digital Sky Survey
(SDSS) to estimate the total satellite galaxy luminosity function down to
$M_{\rm V}=0$. We find that there should be at least $124^{+40}_{-27}$ ($68\%$
CL, statistical error) satellites as bright or brighter than $M_{\rm V}=0$
within 300 kpc of the Sun, with only a weak dependence on MW halo mass. When it
comes online the Large Synoptic Survey Telescope should detect approximately
half of this population. We also show that WDM models infer the same number of
satellites as in $\Lambda$CDM, which will allow us to rule out those models
that produce insufficient substructure to be viable.
[42]
oai:arXiv.org:1806.01718 [pdf] - 1747943
Anisotropies in the astrophysical gravitational-wave background:
Predictions for the detection of compact binaries by LIGO and Virgo
Submitted: 2018-06-05, last modified: 2018-09-10
We develop a detailed anisotropic model for the astrophysical
gravitational-wave background, including binary mergers of two stellar-mass
black holes, two neutron stars, or one of each, which are expected to be the
strongest contributions in the LIGO-Virgo frequency band. The angular spectrum
of the anisotropies, quantified by the $C_\ell$ components, is calculated using
two complementary approaches: (i) a simple, closed-form analytical expression,
and (ii) a detailed numerical study using an all-sky mock light cone galaxy
catalogue from the Millennium simulation. The two approaches are in excellent
agreement at large angular scales, and differ by a factor of order unity at
smaller scales. These anisotropies are considerably larger in amplitude than
e.g. those in the temperature of the cosmic microwave background, confirming
that it is important to model these anisotropies, and indicating that this is a
promising avenue for future theoretical and observational work.
[43]
oai:arXiv.org:1802.06046 [pdf] - 1751752
Anisotropies in the stochastic gravitational-wave background: Formalism
and the cosmic string case
Submitted: 2018-02-16, last modified: 2018-09-10
We develop a powerful analytical formalism for calculating the energy density
of the stochastic gravitational wave background, including a full description
of its anisotropies. This is completely general, and can be applied to any
astrophysical or cosmological source. As an example, we apply these tools to
the case of a network of Nambu-Goto cosmic strings. We find that the angular
spectrum of the anisotropies is relatively insensitive to the choice of model
for the string network, but very sensitive to the value of the string tension
$G\mu$.
[44]
oai:arXiv.org:1708.04247 [pdf] - 1716943
The total satellite population of the Milky Way
Submitted: 2017-08-14, last modified: 2018-07-17
The total number and luminosity function of the population of dwarf galaxies
of the Milky Way (MW) provide important constraints on the nature of the dark
matter and on the astrophysics of galaxy formation at low masses. However, only
a partial census of this population exists because of the flux limits and
restricted sky coverage of existing Galactic surveys. We combine the sample of
satellites recently discovered by the Dark Energy Survey (DES) with the
satellites found in Sloan Digital Sky Survey (SDSS) Data Release 9 (together
these surveys cover nearly half the sky) to estimate the total luminosity
function of satellites down to $M_{\rm V}=0$. We apply a new Bayesian inference
method in which we assume that the radial distribution of satellites
independently of absolute magnitude follows that of subhaloes selected
according to their peak maximum circular velocity. We find that there should be
at least $124^{+40}_{-27}$ (68 per cent CL, statistical error) satellites
brighter than $M_{\rm V}=0$ within $300$ kpc of the Sun. As a result of our use
of new data and better simulations, and a more robust statistical method, we
infer a much smaller population of satellites than reported in previous studies
using earlier SDSS data only; we also address an underestimation of the
uncertainties in earlier work by accounting for stochastic effects. We find
that the inferred number of faint satellites depends only weakly on the assumed
mass of the MW halo and we provide scaling relations to extend our results to
different assumed halo masses and outer radii. We predict that half of our
estimated total satellite population of the MW should be detected by the Large
Synoptic Survey Telescope. The code implementing our estimation method is
available online.
[45]
oai:arXiv.org:1703.10907 [pdf] - 1582275
The Cluster-EAGLE project: global properties of simulated clusters with
resolved galaxies
Barnes, David J.;
Kay, Scott T.;
Bahe, Yannick M.;
Vecchia, Claudio Dalla;
McCarthy, Ian G.;
Schaye, Joop;
Bower, Richard G.;
Jenkins, Adrian;
Thomas, Peter A.;
Schaller, Matthieu;
Crain, Robert A.;
Theuns, Tom;
White, Simon D. M.
Submitted: 2017-03-31, last modified: 2017-07-06
We introduce the Cluster-EAGLE (C-EAGLE) simulation project, a set of
cosmological hydrodynamical zoom simulations of the formation of $30$ galaxy
clusters in the mass range $10^{14}<M_{200}/\mathrm{M}_{\odot}<10^{15.4}$ that
incorporates the Hydrangea sample of Bah\'e et al. (2017). The simulations
adopt the state-of-the-art EAGLE galaxy formation model, with a gas particle
mass of $1.8\times10^{6}\,\mathrm{M}_{\odot}$ and physical softening length of
$0.7\,\mathrm{kpc}$. In this paper, we introduce the sample and present the
low-redshift global properties of the clusters. We calculate the X-ray
properties in a manner consistent with observational techniques, demonstrating
the bias and scatter introduced by using estimated masses. We find the total
stellar content and black hole masses of the clusters to be in good agreement
with the observed relations. However, the clusters are too gas rich, suggesting
that the AGN feedback model is not efficient enough at expelling gas from the
high-redshift progenitors of the clusters. The X-ray properties, such as the
spectroscopic temperature and the soft-band luminosity, and the
Sunyaev-Zel'dovich properties are in reasonable agreement with the observed
relations. However, the clusters have too high central temperatures and
larger-than-observed entropy cores, which is likely driven by the AGN feedback
after the cluster core has formed. The total metal content and its distribution
throughout the ICM are a good match to the observations.
[46]
oai:arXiv.org:1603.04443 [pdf] - 1570104
Knowing the unknowns: uncertainties in simple estimators of galactic
dynamical masses
Campbell, David J. R.;
Frenk, Carlos S.;
Jenkins, Adrian;
Eke, Vincent R.;
Navarro, Julio F.;
Sawala, Till;
Schaller, Matthieu;
Fattahi, Azadeh;
Oman, Kyle A.;
Theuns, Tom
Submitted: 2016-03-14, last modified: 2017-04-21
The observed stellar kinematics of dispersion-supported galaxies are often
used to measure dynamical masses. Recently, several analytical relationships
between the stellar line-of-sight velocity dispersion, the projected (2D) or
deprojected (3D) half-light radius, and the total mass enclosed within the
half-light radius, relying on the spherical Jeans equation, have been proposed.
Here, we make use of the APOSTLE cosmological hydrodynamical simulations of the
Local Group to test the validity and accuracy of such mass estimators for both
dispersion and rotation-supported galaxies, for field and satellite galaxies,
and for galaxies of varying masses, shapes, and velocity dispersion
anisotropies. We find that the mass estimators of Walker et al. and Wolf et al.
are able to recover the masses of dispersion-dominated systems with little
systematic bias, but with a 1-sigma scatter of 25 and 23 percent, respectively.
The error on the estimated mass is dominated by the impact of the 3D shape of
the stellar mass distribution, which is difficult to constrain observationally.
This intrinsic scatter becomes the dominant source of uncertainty in the masses
estimated for galaxies like the dwarf spheroidal (dSph) satellites of the Milky
Way, where the observational errors in their sizes and velocity dispersions are
small. Such scatter may also affect the inner density slopes of dSphs derived
from multiple stellar populations, relaxing the significance with which
Navarro-Frenk-White profiles may be excluded, depending on the degree to which
the relevant properties of the different stellar populations are correlated.
Finally, we derive a new optimal mass estimator that removes the residual
biases and achieves a statistically significant reduction in the scatter to 20
percent overall for dispersion-dominated galaxies, allowing more precise and
accurate mass estimates.
[47]
oai:arXiv.org:1703.10610 [pdf] - 1582257
The Hydrangea simulations: galaxy formation in and around massive
clusters
Bahé, Yannick M.;
Barnes, David J.;
Vecchia, Claudio Dalla;
Kay, Scott T.;
White, Simon D. M.;
McCarthy, Ian G.;
Schaye, Joop;
Bower, Richard G.;
Crain, Robert A.;
Theuns, Tom;
Jenkins, Adrian;
McGee, Sean L.;
Schaller, Matthieu;
Thomas, Peter A.;
Trayford, James W.
Submitted: 2017-03-30
We introduce the Hydrangea simulations, a suite of 24 cosmological
hydrodynamic zoom-in simulations of massive galaxy clusters (M_200c =
10^14-10^15 M_Sun) with baryon particle masses of ~10^6 M_Sun. Designed to
study the impact of the cluster environment on galaxy formation, they are a key
part of the `Cluster-EAGLE' project (Barnes et al. 2017). They use a galaxy
formation model developed for the EAGLE project, which has been shown to yield
both realistic field galaxies and hot gas fractions of galaxy groups consistent
with observations. The total stellar mass content of the simulated clusters
agrees with observations, but central cluster galaxies are too massive, by up
to 0.6 dex. Passive satellite fractions are higher than in the field, and at
stellar masses Mstar > 10^10 M_Sun this environmental effect is quantitatively
consistent with observations. The predicted satellite stellar mass function
matches data from local cluster surveys. Normalized to total mass, there are
fewer low-mass (Mstar < 10^10 M_Sun) galaxies within the virial radius of
clusters than in the field, primarily due to star formation quenching.
Conversely, the simulations predict an overabundance of massive galaxies in
clusters compared to the field that persists to their far outskirts (>
5r_200c). This is caused by a significantly increased stellar mass fraction of
(sub-)haloes in the cluster environment, by up to ~0.3 dex even well beyond
r_200c. Haloes near clusters are also more concentrated than equally massive
field haloes, but these two effects are largely uncorrelated.
[48]
oai:arXiv.org:1610.01159 [pdf] - 1531869
The Auriga Project: the properties and formation mechanisms of disc
galaxies across cosmic time
Submitted: 2016-10-04, last modified: 2017-01-11
We introduce a suite of thirty cosmological magneto-hydrodynamical zoom
simulations of the formation of galaxies in isolated Milky Way mass dark
haloes. These were carried out with the moving mesh code \textlcsc{AREPO},
together with a comprehensive model for galaxy formation physics, including AGN
feedback and magnetic fields, which produces realistic galaxy populations in
large cosmological simulations. We demonstrate that our simulations reproduce a
wide range of present-day observables, in particular, two component disc
dominated galaxies with appropriate stellar masses, sizes, rotation curves,
star formation rates and metallicities. We investigate the driving mechanisms
that set present-day disc sizes/scale lengths, and find that they are related
to the angular momentum of halo material. We show that the largest discs are
produced by quiescent mergers that inspiral into the galaxy and deposit high
angular momentum material into the pre-existing disc, simultaneously increasing
the spin of dark matter and gas in the halo. More violent mergers and strong
AGN feedback play roles in limiting disc size by destroying pre-existing discs
and by suppressing gas accretion onto the outer disc, respectively. The most
important factor that leads to compact discs, however, is simply a low angular
momentum for the halo. In these cases, AGN feedback plays an important role in
limiting central star formation and the formation of a massive bulge.
[49]
oai:arXiv.org:1607.04569 [pdf] - 1531044
The redshift evolution of massive galaxy clusters in the MACSIS
simulations
Submitted: 2016-07-15, last modified: 2016-11-15
We present the MAssive ClusterS and Intercluster Structures (MACSIS) project,
a suite of 390 clusters simulated with baryonic physics that yields realistic
massive galaxy clusters capable of matching a wide range of observed
properties. MACSIS extends the recent BAHAMAS simulation to higher masses,
enabling robust predictions for the redshift evolution of cluster properties
and an assessment of the effect of selecting only the hottest systems. We study
the observable-mass scaling relations and the X-ray luminosity-temperature
relation over the complete observed cluster mass range. As expected, we find
the slope of these scaling relations and the evolution of their normalization
with redshift departs significantly from the self-similar predictions. However,
for a sample of hot clusters with core-excised temperatures
$k_{\rm{B}}T\geq5\,\rm{keV}$ the normalization and slope of the observable-mass
relations and their evolution are significantly closer to self-similar. The
exception is the temperature-mass relation, for which the increased importance
of non-thermal pressure support and biased X-ray temperatures leads to a
greater departure from self-similarity in the hottest systems. As a
consequence, these also affect the slope and evolution of the normalization in
the luminosity-temperature relation. The median hot gas profiles show good
agreement with observational data at $z=0$ and $z=1$, with their evolution
again departing significantly from the self-similar prediction. However,
selecting a hot sample of clusters yields profiles that evolve significantly
closer to the self-similar prediction. In conclusion, our results show that
understanding the selection function is vital for robust calibration of cluster
properties with mass and redshift.
[50]
oai:arXiv.org:1604.07409 [pdf] - 1528049
Substructure and galaxy formation in the Copernicus Complexio warm dark
matter simulations
Submitted: 2016-04-25, last modified: 2016-10-26
We use the Copernicus Complexio (COCO) high resolution $N$-body simulations
to investigate differences in the properties of small-scale structures in the
standard cold dark matter (CDM) model and in a model with a cutoff in the
initial power spectrum of density fluctuations consistent with both a thermally
produced warm dark matter (WDM) particle or a sterile neutrino with mass 7 keV
and leptogenesis parameter $L_6=8.7$. The latter corresponds to the "coldest"
model with this sterile neutrino mass compatible with the identification of the
recently detected 3.5 keV X-ray line as resulting from particle decay. CDM and
WDM predict very different number densities of subhaloes with mass $\leq
10^9\,h^{-1}\,M_\odot$ although they predict similar, nearly universal,
normalised subhalo radial density distributions. Haloes and subhaloes in both
models have cuspy NFW profiles, but WDM subhaloes below the cutoff scale in the
power spectrum (corresponding to maximum circular velocities
$V_{\mathrm{max}}^{z=0} \leq50~\mathrm{kms}^{-1}$) are less concentrated than
their CDM counterparts. We make predictions for observable properties using the
GALFORM semi-analytic model of galaxy formation. Both models predict Milky Way
satellite luminosity functions consistent with observations, although the WDM
model predicts fewer very faint satellites. This model, however, predicts
slightly more UV bright galaxies at redshift $z>7$ than CDM, but both are
consistent with observations. Gravitational lensing offers the best prospect of
distinguishing between the models.
[51]
oai:arXiv.org:1510.01320 [pdf] - 1462935
The EAGLE simulations of galaxy formation: public release of halo and
galaxy catalogues
McAlpine, Stuart;
Helly, John C.;
Schaller, Matthieu;
Trayford, James W.;
Qu, Yan;
Furlong, Michelle;
Bower, Richard G.;
Crain, Robert A.;
Schaye, Joop;
Theuns, Tom;
Vecchia, Claudio Dalla;
Frenk, Carlos S.;
McCarthy, Ian G.;
Jenkins, Adrian;
Rosas-Guevara, Yetli;
White, Simon D. M.;
Baes, Maarten;
Camps, Peter;
Lemson, Gerard
Submitted: 2015-10-05, last modified: 2016-08-25
We present the public data release of halo and galaxy catalogues extracted
from the EAGLE suite of cosmological hydrodynamical simulations of galaxy
formation. These simulations were performed with an enhanced version of the
GADGET code that includes a modified hydrodynamics solver, time-step limiter
and subgrid treatments of baryonic physics, such as stellar mass loss,
element-by-element radiative cooling, star formation and feedback from star
formation and black hole accretion. The simulation suite includes runs
performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with
numerical resolution chosen to marginally resolve the Jeans mass of the gas at
the star formation threshold. The free parameters of the subgrid models for
feedback are calibrated to the redshift z=0 galaxy stellar mass function,
galaxy sizes and black hole mass - stellar mass relation. The simulations have
been shown to match a wide range of observations for present-day and
higher-redshift galaxies. The raw particle data have been used to link galaxies
across redshifts by creating merger trees. The indexing of the tree produces a
simple way to connect a galaxy at one redshift to its progenitors at higher
redshift and to identify its descendants at lower redshift. In this paper we
present a relational database which we are making available for general use. A
large number of properties of haloes and galaxies and their merger trees are
stored in the database, including stellar masses, star formation rates,
metallicities, photometric measurements and mock gri images. Complex queries
can be created to explore the evolution of more than 10^5 galaxies, examples of
which are provided in appendix. (abridged)
[52]
oai:arXiv.org:1512.02219 [pdf] - 1384681
Vertical disc heating in Milky Way-sized galaxies in a cosmological
context
Submitted: 2015-12-07, last modified: 2016-04-02
Vertically extended, high velocity dispersion stellar distributions appear to
be a ubiquitous feature of disc galaxies, and both internal and external
mechanisms have been proposed to be the major driver of their formation.
However, it is unclear to what extent each mechanism can generate such a
distribution, which is likely to depend on the assembly history of the galaxy.
To this end, we perform 16 high resolution cosmological-zoom simulations of
Milky Way-sized galaxies using the state-of-the-art cosmological
magneto-hydrodynamical code \textlcsc{AREPO}, and analyse the evolution of the
vertical kinematics of the stellar disc in connection with various heating
mechanisms. We find that the bar is the dominant heating mechanism in most
cases, whereas spiral arms, radial migration, and adiabatic heating from
mid-plane density growth are all sub-dominant. The strongest source, though
less prevalent than bars, originates from external perturbations from
satellites/sub-halos of masses log$_{10} (M/\rm M_{\odot}) \gtrsim 10$.
However, in many simulations the orbits of newborn star particles become cooler
with time, such that they dominate the shape of the age-velocity dispersion
relation and overall vertical disc structure unless a strong external
perturbation takes place.
[53]
oai:arXiv.org:1505.06436 [pdf] - 1358783
The Copernicus Complexio: a high-resolution view of the small-scale
Universe
Submitted: 2015-05-24, last modified: 2016-02-11
We introduce Copernicus Complexio (COCO), a high-resolution cosmological
N-body simulation of structure formation in the $\Lambda{\rm CDM}{}$ model.
COCO follows an approximately spherical region of radius $\sim 17.4h^{-1}\,{\rm
Mpc}$ embedded in a much larger periodic cube that is followed at lower
resolution. The high resolution volume has a particle mass of
$1.135\times10^5h^{-1}{\rm M}_{\odot}$ (60 times higher than the Millennium-II
simulation). COCO gives the dark matter halo mass function over eight orders of
magnitude in halo mass; it forms $\sim 60$ haloes of galactic size, each
resolved with about 10 million particles. We confirm the power-law character of
the subhalo mass function, $\bar{N}(>\mu)\propto\mu^{-s}$, down to a reduced
subhalo mass $M_{sub}/M_{200}\equiv\mu=10^{-6}$, with a best-fit power-law
index, $s=0.94$, for hosts of mass $\langle M_{200}\rangle=10^{12}h^{-1}{\rm
M}_{\odot}$. The concentration-mass relation of COCO haloes deviates from a
single power law for masses $M_{200}<\textrm{a few}\times 10^{8}h^{-1}{\rm
M}_{\odot}$, where it flattens, in agreement with results by Sanchez-Conde et
al. The host mass invariance of the reduced maximum circular velocity function
of subhaloes, $\nu\equiv V_{max}/V_{200}$, hinted at in previous simulations,
is clearly demonstrated over five orders of magnitude in host mass. Similarly,
we find that the average, normalised radial distribution of subhaloes is
approximately universal (i.e. independent of subhalo mass), as previously
suggested by the Aquarius simulations of individual haloes. Finally, we find
that at fixed physical subhalo size, subhaloes in lower mass hosts typically
have lower central densities than those in higher mass hosts.
[54]
oai:arXiv.org:1411.0311 [pdf] - 1314304
Decaying dark matter: the case for a deep X-ray observation of Draco
Submitted: 2014-11-02, last modified: 2015-11-20
Recent studies of M31, the Galactic centre, and galaxy clusters have made
tentative detections of an X-ray line at ~3.5 keV that could be produced by
decaying dark matter. We use high resolution simulations of the Aquarius
project to predict the likely amplitude of the X-ray decay flux observed in the
GC relative to that observed in M31, and also of the GC relative to other parts
of the Milky Way halo and to dwarf spheroidal galaxies. We show that the
reported detections from M31 and Andromeda are compatible with each other, and
with upper limits arising from high galactic latitude observations, and imply a
decay time {\tau} ~10^28 seconds. We argue that this interpretation can be
tested with deep observations of dwarf spheroidal galaxies: in 95 per cent of
our mock observations, a 1.3 Msec pointed observation of Draco with XMM-Newton
will enable us to discover or rule out at the 3{\sigma} level an X-ray feature
from dark matter decay at 3.5 keV, for decay times {\tau} < 0.8x10^28 sec.
[55]
oai:arXiv.org:1511.01098 [pdf] - 1358941
The APOSTLE simulations: solutions to the Local Group's cosmic puzzles
Sawala, Till;
Frenk, Carlos S.;
Fattahi, Azadeh;
Navarro, Julio F.;
Bower, Richard G.;
Crain, Robert A.;
Vecchia, Claudio Dalla;
Furlong, Michelle;
Helly, John. C.;
Jenkins, Adrian;
Oman, Kyle A.;
Schaller, Matthieu;
Schaye, Joop;
Theuns, Tom;
Trayford, James;
White, Simon D. M.
Submitted: 2015-11-03
The Local Group of galaxies offer some of the most discriminating tests of
models of cosmic structure formation. For example, observations of the Milky
Way (MW) and Andromeda satellite populations appear to be in disagreement with
N-body simulations of the "Lambda Cold Dark Matter" ({\Lambda}CDM) model: there
are far fewer satellite galaxies than substructures in cold dark matter halos
(the "missing satellites" problem); dwarf galaxies seem to avoid the most
massive substructures (the "too-big-to-fail" problem); and the brightest
satellites appear to orbit their host galaxies on a thin plane (the "planes of
satellites" problem). Here we present results from APOSTLE (A Project Of
Simulating The Local Environment), a suite of cosmological hydrodynamic
simulations of twelve volumes selected to match the kinematics of the Local
Group (LG) members. Applying the Eagle code to the LG environment, we find that
our simulations match the observed abundance of LG galaxies, including the
satellite galaxies of the MW and Andromeda. Due to changes to the structure of
halos and the evolution in the LG environment, the simulations reproduce the
observed relation between stellar mass and velocity dispersion of individual
dwarf spheroidal galaxies without necessitating the formation of cores in their
dark matter profiles. Satellite systems form with a range of spatial
anisotropies, including one similar to that of the MW, confirming that such a
configuration is not unexpected in {\Lambda}CDM. Finally, based on the observed
velocity dispersion, size, and stellar mass, we provide new estimates of the
maximum circular velocity for the halos of nine MW dwarf spheroidals.
[56]
oai:arXiv.org:1507.03643 [pdf] - 1351450
The APOSTLE project: Local Group kinematic mass constraints and
simulation candidate selection
Fattahi, Azadeh;
Navarro, Julio F.;
Sawala, Till;
Frenk, Carlos S.;
Oman, Kyle A.;
Crain, Robert A.;
Furlong, Michelle;
Schaller, Matthieu;
Schaye, Joop;
Theuns, Tom;
Jenkins, Adrian
Submitted: 2015-07-13
We use a large sample of isolated dark matter halo pairs drawn from
cosmological N-body simulations to identify candidate systems whose kinematics
match that of the Local Group of Galaxies (LG). We find, in agreement with the
"timing argument" and earlier work, that the separation and approach velocity
of the Milky Way (MW) and Andromeda (M31) galaxies favour a total mass for the
pair of ~ 5*10^12 M_sun. A mass this large, however, is difficult to reconcile
with the small relative tangential velocity of the pair, as well as with the
small deceleration from the Hubble flow observed for the most distant LG
members. Halo pairs that match these three criteria have average masses a
factor of ~2 times smaller than suggested by the timing argument, but with
large dispersion, spanning more than a decade in mass. Guided by these results,
we have selected 12 halo pairs with total mass in the range 1.6-3.6 *10^12
M_sun for the APOSTLE project (A Project Of Simulations of The Local
Environment), a suite of resimulations at various numerical resolution levels
(reaching up to ~10^4 M_sun per gas particle) that use the hydrodynamical code
and subgrid physics developed for the EAGLE project. These simulations
reproduce, by construction, the main kinematics of the MW-M31 pair, and produce
satellite populations whose overall number, luminosities, and kinematics are in
good agreement with observations of the MW and M31 companions. These
diagnostics are sensitive to the total mass assumed for the MW-M31 pair;
indeed, the LG satellite population would be quite difficult to reproduce for
pair masses as high as indicated by the timing argument. The APOSTLE candidate
systems thus provide an excellent testbed to confront directly many of the
predictions of the Lambda-CDM cosmology with observations of our local
Universe.
[57]
oai:arXiv.org:1507.01998 [pdf] - 1311969
The COpernicus COmplexio: Statistical Properties of Warm Dark Matter
Haloes
Submitted: 2015-07-07, last modified: 2015-07-09
The recent detection of a 3.5 keV X-ray line from the centres of galaxies and
clusters by Bulbul et al. (2014a) and Boyarsky et al. (2014a) has been
interpreted as emission from the decay of 7 keV sterile neutrinos which could
make up the (warm) dark matter (WDM). As part of the COpernicus COmplexio
(COCO) programme, we investigate the properties of dark matter haloes formed in
a high-resolution cosmological $N$-body simulation from initial conditions
similar to those expected in a universe in which the dark matter consists of 7
keV sterile neutrinos. This simulation and its cold dark matter (CDM)
counterpart have $\sim13.4$bn particles, each of mass $\sim 10^5\, h^{-1}
M_\odot$, providing detailed information about halo structure and evolution
down to dwarf galaxy mass scales. Non-linear structure formation on small
scales ($M_{200}\, \leq\, 2 \times 10^9\,h^{-1}\,M_\odot$) begins slightly
later in COCO-Warm than in COCO-Cold. The halo mass function at the present day
in the WDM model begins to drop below its CDM counterpart at a mass $\sim 2
\times 10^{9}\,h^{-1}\,M_\odot$ and declines very rapidly towards lower masses
so that there are five times fewer haloes of mass $M_{200}=
10^{8}\,h^{-1}\,M_\odot$ in COCO-Warm than in COCO-Cold. Halo concentrations on
dwarf galaxy scales are correspondingly smaller in COCO-Warm, and we provide a
simple functional form that describes its evolution with redshift. The shapes
of haloes are similar in the two cases, but the smallest haloes in COCO-Warm
rotate slightly more slowly than their CDM counterparts.
[58]
oai:arXiv.org:1407.5627 [pdf] - 1241728
Surface photometry of BCGs and intracluster stars in Lambda-CDM
Submitted: 2014-07-21, last modified: 2015-07-05
We simulate the phase-space distribution of stellar mass in 9 massive
Lambda-CDM galaxy clusters by applying the semi-analytic particle tagging
method of Cooper et al. to the Phoenix suite of high-resolution N-body
simulations (M200 = 7.5 to 33 x 10^14 Msol). The resulting surface brightness
(SB) profiles of brightest cluster galaxies (BCGs) match well to observations.
On average, stars formed in galaxies accreted by the BCG account for ~90 per
cent of its total mass (the remainder is formed in situ). In circular
BCG-centred apertures, the superposition of multiple debris clouds (each ~10
per cent of the total BCG mass) from different progenitors can result in an
extensive outer diffuse component, qualitatively similar to a 'cD envelope'.
These clouds typically originate from tidal stripping at z < 1 and comprise
both streams and the extended envelopes of other massive galaxies in the
cluster. Stars at very low SB contribute a significant fraction of the total
cluster stellar mass budget: in the central 1 Mpc^2 of a z ~ 0.15 cluster
imaged at SDSS-like resolution, our fiducial model predicts 80-95 per cent of
stellar mass below a SB of mu_V = 26.5 mag arcsec^2 is associated with accreted
stars in the envelope of the BCG. The ratio of BCG stellar mass (including this
diffuse component) to total cluster stellar mass is ~30 per cent.
[59]
oai:arXiv.org:1409.8617 [pdf] - 1241307
Baryon effects on the internal structure of LCDM halos in the EAGLE
simulations
Schaller, Matthieu;
Frenk, Carlos S.;
Bower, Richard G.;
Theuns, Tom;
Jenkins, Adrian;
Schaye, Joop;
Crain, Robert A.;
Furlong, Michelle;
Vecchia, Claudio Dalla;
McCarthy, I. G.
Submitted: 2014-09-30, last modified: 2015-05-25
We investigate the internal structure and density profiles of halos of mass
$10^{10}-10^{14}~M_\odot$ in the Evolution and Assembly of Galaxies and their
Environment (EAGLE) simulations. These follow the formation of galaxies in a
$\Lambda$CDM Universe and include a treatment of the baryon physics thought to
be relevant. The EAGLE simulations reproduce the observed present-day galaxy
stellar mass function, as well as many other properties of the galaxy
population as a function of time. We find significant differences between the
masses of halos in the EAGLE simulations and in simulations that follow only
the dark matter component. Nevertheless, halos are well described by the
Navarro-Frenk-White (NFW) density profile at radii larger than ~5% of the
virial radius but, closer to the centre, the presence of stars can produce
cuspier profiles. Central enhancements in the total mass profile are most
important in halos of mass $10^{12}-10^{13}M_\odot$, where the stellar fraction
peaks. Over the radial range where they are well resolved, the resulting galaxy
rotation curves are in very good agreement with observational data for galaxies
with stellar mass $M_*<5\times10^{10}M_\odot$. We present an empirical fitting
function that describes the total mass profiles and show that its parameters
are strongly correlated with halo mass.
[60]
oai:arXiv.org:1410.3485 [pdf] - 1280754
Evolution of galaxy stellar masses and star formation rates in the EAGLE
simulations
Furlong, M.;
Bower, R. G.;
Theuns, T.;
Schaye, J.;
Crain, R. A.;
Schaller, M.;
Vecchia, C. Dalla;
Frenk, C. S.;
McCarthy, I. G.;
Helly, J.;
Jenkins, A.;
Rosas-Guevara, Y. M.
Submitted: 2014-10-13, last modified: 2015-04-22
We investigate the evolution of galaxy masses and star formation rates in the
Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations.
These comprise a suite of hydrodynamical simulations in a $\Lambda$CDM
cosmogony with subgrid models for radiative cooling, star formation, stellar
mass loss, and feedback from stars and accreting black holes. The subgrid
feedback was calibrated to reproduce the observed present-day galaxy stellar
mass function and galaxy sizes. Here we demonstrate that the simulations
reproduce the observed growth of the stellar mass density to within 20 per
cent. The simulation also tracks the observed evolution of the galaxy stellar
mass function out to redshift z = 7, with differences comparable to the
plausible uncertainties in the interpretation of the data. Just as with
observed galaxies, the specific star formation rates of simulated galaxies are
bimodal, with distinct star forming and passive sequences. The specific star
formation rates of star forming galaxies are typically 0.2 to 0.4 dex lower
than observed, but the evolution of the rates track the observations closely.
The unprecedented level of agreement between simulation and data makes EAGLE a
powerful resource to understand the physical processes that govern galaxy
formation.
[61]
oai:arXiv.org:1501.01311 [pdf] - 1223890
The EAGLE simulations of galaxy formation: calibration of subgrid
physics and model variations
Crain, Robert A.;
Schaye, Joop;
Bower, Richard G.;
Furlong, Michelle;
Schaller, Matthieu;
Theuns, Tom;
Vecchia, Claudio Dalla;
Frenk, Carlos S.;
McCarthy, Ian G.;
Helly, John C.;
Jenkins, Adrian;
Rosas-Guevara, Yetli M.;
White, Simon D. M.;
Trayford, James W.
Submitted: 2015-01-06, last modified: 2015-04-03
We present results from thirteen cosmological simulations that explore the
parameter space of the "Evolution and Assembly of GaLaxies and their
Environments" (EAGLE) simulation project. Four of the simulations follow the
evolution of a periodic cube L = 50 cMpc on a side, and each employs a
different subgrid model of the energetic feedback associated with star
formation. The relevant parameters were adjusted so that the simulations each
reproduce the observed galaxy stellar mass function at z = 0.1. Three of the
simulations fail to form disc galaxies as extended as observed, and we show
analytically that this is a consequence of numerical radiative losses that
reduce the efficiency of stellar feedback in high-density gas. Such losses are
greatly reduced in the fourth simulation - the EAGLE reference model - by
injecting more energy in higher density gas. This model produces galaxies with
the observed size distribution, and also reproduces many galaxy scaling
relations. In the remaining nine simulations, a single parameter or process of
the reference model was varied at a time. We find that the properties of
galaxies with stellar mass <~ M* (the "knee" of the galaxy stellar mass
function) are largely governed by feedback associated with star formation,
while those of more massive galaxies are also controlled by feedback from
accretion onto their central black holes. Both processes must be efficient in
order to reproduce the observed galaxy population. In general, simulations that
have been calibrated to reproduce the low-redshift galaxy stellar mass function
will still not form realistic galaxies, but the additional requirement that
galaxy sizes be acceptable leads to agreement with a large range of
observables.
[62]
oai:arXiv.org:1412.2748 [pdf] - 905335
Local Group galaxies emerge from the dark
Sawala, Till;
Frenk, Carlos S.;
Fattahi, Azadeh;
Navarro, Julio F.;
Bower, Richard G.;
Crain, Robert A.;
Vecchia, Claudio Dalla;
Furlong, Michelle;
Helly, John. C.;
Jenkins, Adrian;
Oman, Kyle A.;
Schaller, Matthieu;
Schaye, Joop;
Theuns, Tom;
Trayford, James;
White, Simon D. M.
Submitted: 2014-12-08
The "Lambda Cold Dark Matter" (LCDM) model of cosmic structure formation is
eminently falsifiable: once its parameters are fixed on large scales, it
becomes testable in the nearby Universe. Observations within our Local Group of
galaxies, including the satellite populations of the Milky Way and Andromeda,
appear to contradict LCDM predictions: there are far fewer satellite galaxies
than dark matter halos (the "missing satellites" problem), galaxies seem to
avoid the largest substructures (the "too big to fail" problem), and the
brightest satellites appear to orbit their host galaxies on a thin plane (the
"planes of satellites" problem). We present results from the first hydrodynamic
simulations of the Local Group that match the observed abundance of galaxies.
We find that when baryonic and dark matter are followed simultaneously in the
context of a realistic galaxy formation model, all three "problems" are
resolved within the LCDM paradigm.
[63]
oai:arXiv.org:1407.7040 [pdf] - 1215926
The EAGLE project: Simulating the evolution and assembly of galaxies and
their environments
Schaye, Joop;
Crain, Robert A.;
Bower, Richard G.;
Furlong, Michelle;
Schaller, Matthieu;
Theuns, Tom;
Vecchia, Claudio Dalla;
Frenk, Carlos S.;
McCarthy, I. G.;
Helly, John C.;
Jenkins, Adrian;
Rosas-Guevara, Y. M.;
White, Simon D. M.;
Baes, Maarten;
Booth, C. M.;
Camps, Peter;
Navarro, Julio F.;
Qu, Yan;
Rahmati, Alireza;
Sawala, Till;
Thomas, Peter A.;
Trayford, James
Submitted: 2014-07-25, last modified: 2014-10-02
We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical
simulations that follow the formation of galaxies and black holes in
representative volumes. We discuss the limitations of such simulations in light
of their finite resolution and poorly constrained subgrid physics, and how
these affect their predictive power. One major improvement is our treatment of
feedback from massive stars and AGN in which thermal energy is injected into
the gas without the need to turn off cooling or hydrodynamical forces, allowing
winds to develop without predetermined speed or mass loading factors. Because
the feedback efficiencies cannot be predicted from first principles, we
calibrate them to the z~0 galaxy stellar mass function and the amplitude of the
galaxy-central black hole mass relation, also taking galaxy sizes into account.
The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the
full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement
close to that attained by semi-analytic models, and unprecedented for
hydrodynamical simulations. We compare our results to a representative set of
low-redshift observables not considered in the calibration, and find good
agreement with the observed galaxy specific star formation rates, passive
fractions, Tully-Fisher relation, total stellar luminosities of galaxy
clusters, and column density distributions of intergalactic CIV and OVI. While
the mass-metallicity relations for gas and stars are consistent with
observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at
lower masses. The gas fractions and temperatures are too high for clusters of
galaxies, but for groups these discrepancies can be resolved by adopting a
higher heating temperature in the subgrid prescription for AGN feedback. EAGLE
constitutes a valuable new resource for studies of galaxy formation.
[64]
oai:arXiv.org:1405.5216 [pdf] - 1209644
Dwarf galaxies in CDM and SIDM with baryons: observational probes of the
nature of dark matter
Submitted: 2014-05-20, last modified: 2014-09-20
We present the first cosmological simulations of dwarf galaxies, which
include dark matter self-interactions and baryons. We study two dwarf galaxies
within cold dark matter, and four different elastic self-interacting scenarios
with constant and velocity-dependent cross sections, motivated by a new force
in the hidden dark matter sector. Our highest resolution simulation has a
baryonic mass resolution of $1.8\times 10^2\,{\rm M}_\odot$ and a gravitational
softening length of $34\,{\rm pc}$ at $z=0$. In this first study we focus on
the regime of mostly isolated dwarf galaxies with halo masses
$\sim10^{10}\,{\rm M}_\odot$ where dark matter dynamically dominates even at
sub-kpc scales. We find that while the global properties of galaxies of this
scale are minimally affected by allowed self-interactions, their internal
structures change significantly if the cross section is large enough within the
inner sub-kpc region. In these dark-matter-dominated systems, self-scattering
ties the shape of the stellar distribution to that of the dark matter
distribution. In particular, we find that the stellar core radius is closely
related to the dark matter core radius generated by self-interactions. Dark
matter collisions lead to dwarf galaxies with larger stellar cores and smaller
stellar central densities compared to the cold dark matter case. The central
metallicity within $1\,{\rm kpc}$ is also larger by up to $\sim 15\%$ in the
former case. We conclude that the mass distribution, and characteristics of the
central stars in dwarf galaxies can potentially be used to probe the
self-interacting nature of dark matter.
[65]
oai:arXiv.org:1409.0723 [pdf] - 1216637
Cosmological simulations of galaxy clusters with feedback from active
galactic nuclei: profiles and scaling relations
Submitted: 2014-09-02
We present results from a new set of 30 cosmological simulations of galaxy
clusters, including the effects of radiative cooling, star formation, supernova
feedback, black hole growth and AGN feedback. We first demonstrate that our AGN
model is capable of reproducing the observed cluster pressure profile at
redshift, z~0, once the AGN heating temperature of the targeted particles is
made to scale with the final virial temperature of the halo. This allows the
ejected gas to reach larger radii in higher-mass clusters than would be
possible had a fixed heating temperature been used. Such a model also
successfully reduces the star formation rate in brightest cluster galaxies and
broadly reproduces a number of other observational properties at low redshift,
including baryon, gas and star fractions; entropy profiles outside the core;
and the X-ray luminosity-mass relation. Our results are consistent with the
notion that the excess entropy is generated via selective removal of the
densest material through radiative cooling; supernova and AGN feedback largely
serve as regulation mechanisms, moving heated gas out of galaxies and away from
cluster cores. However, our simulations fail to address a number of serious
issues; for example, they are incapable of reproducing the shape and diversity
of the observed entropy profiles within the core region. We also show that the
stellar and black hole masses are sensitive to numerical resolution,
particularly the gravitational softening length; a smaller value leads to more
efficient black hole growth at early times and a smaller central galaxy.
[66]
oai:arXiv.org:1406.6362 [pdf] - 1351404
The chosen few: the low mass halos that host faint galaxies
Sawala, Till;
Frenk, Carlos S.;
Fattahi, Azadeh;
Navarro, Julio F.;
Theuns, Tom;
Bower, Richard G.;
Crain, Robert A.;
Furlong, Michelle;
Jenkins, Adrian;
Schaller, Matthieu;
Schaye, Joop
Submitted: 2014-06-24
Since reionization prevents star formation in most halos below 3 x 10^9 solar
masses, dwarf galaxies only populate a fraction of existing dark matter halos.
We use hydrodynamic cosmological simulations of the Local Group to study the
discriminating factors for galaxy formation in the early Universe and connect
them to the present-day properties of galaxies and halos. A combination of
selection effects related to reionization, and the subsequent evolution of
halos in different environments, introduces strong biases between the
population of halos that host dwarf galaxies, and the total halo population.
Halos that host galaxies formed earlier and are more concentrated. In addition,
halos more affected by tidal stripping are more likely to host a galaxy for a
given mass or maximum circular velocity, vmax, today. Consequently, satellite
halos are populated more frequently than field halos, and satellite halos of
10^8 - 10^9 solar masses or vmax of 12 - 20 km/s, similar to the Local Group
dwarf spheroidals, have experienced a greater than average reduction in both
mass and vmax after infall. They are on closer, more radial orbits with higher
infall velocities and earlier infall times. Together, these effects make dwarf
galaxies highly biased tracers of the underlying dark matter distribution.
[67]
oai:arXiv.org:1308.1399 [pdf] - 824986
The properties of warm dark matter haloes
Submitted: 2013-08-06, last modified: 2014-05-19
Well-motivated elementary particle candidates for the dark matter, such as
the sterile neutrino, behave as warm dark matter (WDM).For particle masses of
order a keV, free streaming produces a cutoff in the linear fluctuation power
spectrum at a scale corresponding to dwarf galaxies. We investigate the
abundance and structure of WDM haloes and subhaloes on these scales using high
resolution cosmological N-body simulations of galactic haloes of mass similar
to the Milky Way's. On scales larger than the free-streaming cutoff, the
initial conditions have the same power spectrum and phases as one of the cold
dark matter (CDM) haloes previously simulated by Springel et al as part of the
Virgo consortium Aquarius project. We have simulated four haloes with WDM
particle masses in the range 1.4-2.3keV and, for one case, we have carried out
further simulations at varying resolution. N-body simulations in which the
power spectrum cutoff is resolved are known to undergo artificial fragmentation
in filaments producing spurious clumps which, for small masses (<10^7Msun in
our case) outnumber genuine haloes. We have developed a robust algorithm to
identify these spurious objects and remove them from our halo catalogues. We
find that the WDM subhalo mass function is suppressed by well over an order
magnitude relative to the CDM case for masses <10^9Msun. Requiring that there
should be at least as many subhaloes as there are observed satellites in the
Milky Way leads to a conservative lower limit to the (thermal equivalent) WDM
particle mass of ~1.5\rmn{keV}. WDM haloes and subhaloes have cuspy density
distributions that are well described by NFW or Einasto profiles. Their central
densities are lower for lower WDM particle masses and none of the models we
have considered suffer from the "too big to fail" problem recently highlighted
by Boylan-Kolchin et al.
[68]
oai:arXiv.org:1404.3724 [pdf] - 1208978
Bent by baryons: the low mass galaxy-halo relation
Sawala, Till;
Frenk, Carlos S.;
Fattahi, Azadeh;
Navarro, Julio F.;
Bower, Richard G.;
Crain, Robert A.;
Vecchia, Claudio Dalla;
Furlong, Michelle;
Jenkins, Adrian;
McCarthy, Ian G.;
Qu, Yan;
Schaller, Matthieu;
Schaye, Joop;
Theuns, Tom
Submitted: 2014-04-14
The relation between galaxies and dark matter halos is of vital importance
for evaluating theoretical predictions of structure formation and galaxy
formation physics. We show that the widely used method of abundance matching
based on dark matter only simulations fails at the low mass end because two of
its underlying assumptions are broken: only a small fraction of low mass (below
10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower
rate due to the effect of baryons. In this regime, reliance on dark matter only
simulations for abundance matching is neither accurate nor self-consistent. We
find that the reported discrepancy between observational estimates of the halo
masses of dwarf galaxies and the values predicted by abundance matching does
not point to a failure of LCDM, but simply to a failure to account for baryonic
effects. Our results also imply that the Local Group contains only a few
hundred observable galaxies in contrast with the thousands of faint dwarfs that
abundance matching would suggest. We show how relations derived from abundance
matching can be corrected, so that they can be used self-consistently to
calibrate models of galaxy formation.
[69]
oai:arXiv.org:1203.5339 [pdf] - 1117525
The journey of QSO haloes from z=6 to the present
Submitted: 2012-03-23, last modified: 2013-09-25
We apply a recently developed scaling technique to the Millennium-XXL, one of
the largest cosmological N-body simulations carried out to date 3x10^11
particles within a cube of volume ~70Gpc^3). This allows us to investigate the
cosmological parameter dependence of the mass and evolution of haloes in the
extreme high-mass tail of the z=6 distribution. We assume these objects to be
likely hosts for the population of rare but ultraluminous high-redshift quasars
discovered by the Sloan Digital Sky Survey. Haloes with a similar abundance to
these quasars have a median mass of 9x10^12 Msun in the currently preferred
cosmology, but do not evolve into equally extreme objects at z=0. Rather, their
descendants span the full range conventionally assigned to present-day
clusters, 6x10^13 to 2.5x10^15 Msun for this same cosmology. The masses both at
z=6 and at z=0 shift up or down by factors exceeding two if cosmological
parameters are pushed to the boundaries of the range discussed in published
interpretations of data from the WMAP satellite. The main factor determining
the future growth of a high-mass z=6 halo is the mean overdensity of its
environment on scales of 7 to 14 Mpc, and descendant masses can be predicted 6
to 8 times more accurately if this density is known than if it is not. All
these features are not unique to extreme high-z haloes, but are generic to
hierarchical growth. Finally, we find that extreme haloes at z=6 typically
acquired about half of their total mass in the preceding 100 Myr, implying very
large recent accretion rates which may be related to the large black hole
masses and high luminosities of the SDSS quasars.
[70]
oai:arXiv.org:1203.3216 [pdf] - 1117317
Scaling relations for galaxy clusters in the Millennium-XXL simulation
Submitted: 2012-03-14, last modified: 2013-09-25
We present a very large high-resolution cosmological N-body simulation, the
Millennium-XXL or MXXL, which uses 303 billion particles to represent the
formation of dark matter structures throughout a 4.1Gpc box in a LambdaCDM
cosmology. We create sky maps and identify large samples of galaxy clusters
using surrogates for four different observables: richness estimated from galaxy
surveys, X-ray luminosity, integrated Sunyaev-Zeldovich signal, and lensing
mass. The unprecedented combination of volume and resolution allows us to
explore in detail how these observables scale with each other and with cluster
mass. The scatter correlates between different mass-observable relations
because of common sensitivities to the internal structure, orientation and
environment of clusters, as well as to line-of-sight superposition of
uncorrelated structure. We show that this can account for the apparent
discrepancies uncovered recently between the mean thermal SZ signals measured
for optically and X-ray selected clusters by stacking data from the Planck
satellite. Related systematics can also affect inferences from extreme clusters
detected at high redshift. Our results illustrate that cosmological conclusions
from galaxy cluster surveys depend critically on proper modelling, not only of
the relevant physics, but also of the full distribution of the observables and
of the selection biases induced by cluster identification procedures.
[71]
oai:arXiv.org:1306.5968 [pdf] - 1172284
A new way of setting the phases for cosmological multi-scale Gaussian
initial conditions
Submitted: 2013-06-25
We describe how to define an extremely large discrete realisation of a
Gaussian white noise field that has a hierarchical structure and the property
that the value of any part of the field can be computed quickly. Tiny
subregions of such a field can be used to set the phase information for
Gaussian initial conditions for individual cosmological simulations of
structure formation. This approach has several attractive features: (i) the
hierarchical structure based on an octree is particular well suited for
generating follow up resimulation or zoom initial conditions; (ii) the phases
are defined for all relevant physical scales in advance so that resimulation
initial conditions are, by construction, consistent both with their parent
simulation and with each other; (iii) the field can easily be made public by
releasing a code to compute it - once public, phase information can be shared
or published by specifying a spatial location within the realisation. In this
paper we describe the principles behind creating such realisations. We define
an example called PANPHASIA and in a companion paper on arXiv (Jenkins &
Booth), make public a code to compute it. With fifty octree levels Panphasia
spans a factor of more than 10^15 in linear scale - a range that significantly
exceeds the ratio of the current Hubble radius to the putative CDM
free-streaming scale. We show how to modify a code used for making cosmological
and resimulation initial conditions so that it can take the phase information
from Panphasia and, using this code, we demonstrate that it is possible to make
good quality resimulation initial conditions. We define a convention for
publishing phase information from Panphasia and publish the initial phases for
several of the Virgo Consortium's most recent cosmological simulations
including the 303 billion particle MXXL simulation.
[72]
oai:arXiv.org:1306.5771 [pdf] - 684409
Panphasia: a user guide
Submitted: 2013-06-24
We make a very large realisation of a Gaussian white noise field, called
PANPHASIA, public by releasing software that computes this field. Panphasia is
designed specifically for setting up Gaussian initial conditions for
cosmological simulations and resimulations of structure formation. We make
available both software to compute the field itself and codes to illustrate
applications including a modified version of a public serial initial conditions
generator. We document the software and present the results of a few basic
tests of the field. The properties and method of construction of Panphasia are
described in full in a companion paper Jenkins 2013.
[73]
oai:arXiv.org:1206.6495 [pdf] - 656088
The abundance of (not just) dark matter haloes
Submitted: 2012-06-27, last modified: 2013-04-24
We study the effect of baryons on the abundance of structures and
substructures in a Lambda-CDM cosmology, using a pair of high resolution
cosmological simulations from the GIMIC project. Both simulations use identical
initial conditions, but while one contains only dark matter, the other also
includes baryons. We find that gas pressure, reionisation, supernova feedback,
stripping, and truncated accretion systematically reduce the total mass and the
abundance of structures below ~10^12 solar masses compared to the pure dark
matter simulation. Taking this into account and adopting an appropriate
detection threshold lowers the abundance of observed galaxies with maximum
circular velocities below 100 km/s, significantly reducing the reported
discrepancy between Lambda-CDM and the measured HI velocity function of the
ALFALFA survey. We also show that the stellar-to-total mass ratios of galaxies
with stellar masses of ~10^5 - 10^7 solar masses inferred from abundance
matching of the (sub)halo mass function to the observed galaxy mass function
increase by a factor of ~2. In addition, we find that an important fraction of
low-mass subhaloes are completely devoid of stars. Accounting for the presence
of dark subhaloes below 10^10 solar masses further reduces the abundance of
observable objects, and leads to an additional increase in the inferred
stellar-to-total mass ratio by factors of 2 - 10 for galaxies in haloes of 10^9
- 10^10 solar masses. This largely reconciles the abundance matching results
with the kinematics of individual dwarf galaxies in Lambda-CDM. We propose
approximate corrections to the masses of objects derived from pure dark matter
calculations to account for baryonic effects.
[74]
oai:arXiv.org:1209.3018 [pdf] - 1151396
Dark Matter Halo Merger Histories Beyond Cold Dark Matter: I - Methods
and Application to Warm Dark Matter
Submitted: 2012-09-13, last modified: 2012-10-07
We describe a methodology to accurately compute halo mass functions,
progenitor mass functions, merger rates and merger trees in non-cold dark
matter universes using a self-consistent treatment of the generalized extended
Press-Schechter formalism. Our approach permits rapid exploration of the
subhalo population of galactic halos in dark matter models with a variety of
different particle properties or universes with rolling, truncated, or more
complicated power spectra. We make detailed comparisons of analytically derived
mass functions and merger histories with recent warm dark matter cosmological
N-body simulations, and find excellent agreement. We show that, once the
accretion of smoothly distributed matter is accounted for, coarse-grained
statistics such as the mass accretion history of halos can be almost
indistinguishable between cold and warm dark matter cases. However, the halo
mass function and progenitor mass functions differ significantly, with the warm
dark matter cases being strongly suppressed below the free-streaming scale of
the dark matter. We demonstrate the importance of using the correct solution
for the excursion set barrier first-crossing distribution in warm dark matter -
if the solution for a flat barrier is used instead the truncation of the halo
mass function is much slower, leading to an overestimate of the number of low
mass halos.
[75]
oai:arXiv.org:1208.5040 [pdf] - 554406
Gyroscopic Inflation
Submitted: 2012-08-24
We propose a new framework for multi-field inflation in which a nearly
constant potential energy is maintained during inflation before decreasing
rapidly, in a manner analogous to a classical top spinning upright for a long
time before falling down. We provide the simplest realization of such dynamics
as a well-controlled, weakly-coupled effective field theory with a global shift
symmetry. Nonperturbative quantum gravitational effects, which break global
symmetries, are suppressed in this model. Primordial gravitational waves may be
within experimental reach.
[76]
oai:arXiv.org:1201.1940 [pdf] - 558130
The Phoenix Project: the Dark Side of Rich Galaxy Clusters
Submitted: 2012-01-09, last modified: 2012-06-25
[abridged] We introduce the Phoenix Project, a set of $\Lambda$CDM
simulations of the dark matter component of nine rich galaxy clusters. Each
cluster is simulated at least at two different numerical resolutions. For eight
of them, the highest resolution corresponds to $\sim 130$ million particles
within the virial radius, while for one this number is over one billion. We
study the structure and substructure of these systems and contrast them with
six galaxy-sized dark matter haloes from the Aquarius Project, simulated at
comparable resolution. This comparison highlights the approximate mass
invariance of CDM halo structure and substructure. We find little difference in
the spherically-averaged mass, pseudo-phase-space density, and velocity
anisotropy profiles of Aquarius and Phoenix haloes. When scaled to the virial
properties of the host halo, the abundance and radial distribution of subhaloes
are also very similar, despite the fact that Aquarius and Phoenix haloes differ
by roughly three decades in virial mass. The most notable difference is that
cluster haloes have been assembled more recently and are thus significantly
less relaxed than galaxy haloes, which leads to decreased regularity, increased
halo-to-halo scatter and sizable deviations from the mean trends. This
accentuates the effects of the strong asphericity of individual clusters on
surface density profiles, which may vary by up to a factor of three at a given
radius, depending on projection. The high apparent concentration reported for
some strong-lensing clusters might very well reflect these effects. A more
recent assembly also explains why substructure in some Phoenix haloes is
slightly more abundant than in Aquarius, especially in the inner regions.
Resolved subhaloes nevertheless contribute only $11 \pm 3%$ of the virial mass
in Phoenix clusters. .
[77]
oai:arXiv.org:1112.0315 [pdf] - 1092092
The Aquila comparison Project: The Effects of Feedback and Numerical
Methods on Simulations of Galaxy Formation
Scannapieco, C.;
Wadepuhl, M.;
Parry, O. H.;
Navarro, J. F.;
Jenkins, A.;
Springel, V.;
Teyssier, R.;
Carlson, E.;
Couchman, H. M. P.;
Crain, R. A.;
Vecchia, C. Dalla;
Frenk, C. S.;
Kobayashi, C.;
Monaco, P.;
Murante, G.;
Okamoto, T.;
Quinn, T.;
Schaye, J.;
Stinson, G. S.;
Theuns, T.;
Wadsley, J.;
White, S. D. M.;
Woods, R.
Submitted: 2011-12-01, last modified: 2012-04-19
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.
[78]
oai:arXiv.org:1104.2929 [pdf] - 966686
The haloes of bright satellite galaxies in a warm dark matter universe
Lovell, Mark;
Eke, Vincent;
Frenk, Carlos;
Gao, Liang;
Jenkins, Adrian;
Theuns, Tom;
Wang, Jie;
White, Simon;
Boyarsky, Alexey;
Ruchayskiy, Oleg
Submitted: 2011-04-14, last modified: 2012-01-16
High resolution N-body simulations of galactic cold dark matter haloes
indicate that we should expect to find a few satellite galaxies around the
Milky Way whose haloes have a maximum circular velocity in excess of 40 kms.
Yet, with the exception of the Magellanic Clouds and the Sagittarius dwarf,
which likely reside in subhaloes with significantly larger velocities than
this, the bright satellites of the Milky Way all appear to reside in subhaloes
with maximum circular velocities below 40 kms. As recently highlighted by
Boylan-Kolchin et al., this discrepancy implies that the majority of the most
massive subhaloes within a cold dark matter galactic halo are much too
concentrated to be consistent with the kinematic data for the bright Milky Way
satellites. Here we show that no such discrepancy exists if haloes are made of
warm, rather than cold dark matter because these haloes are less concentrated
on account of their typically later formation epochs. Warm dark matter is one
of several possible explanations for the observed kinematics of the satellites.
[79]
oai:arXiv.org:1107.1916 [pdf] - 457340
Where will supersymmetric dark matter first be seen?
Submitted: 2011-07-10, last modified: 2011-09-16
If the dark matter consists of supersymmetric particles, Gamma-ray
observatories such as the Large Area Telescope aboard the Fermi satellite may
detect annihilation radiation from the haloes of galaxies and galaxy clusters.
Much recent effort has been devoted to searching for this signal around the
Milky Way's dwarf satellites. Using a new suite of high-resolution simulations
of galaxy cluster haloes (the Phoenix Project), together with the Aquarius
simulations of Milky-Way-like galaxy haloes, we show that higher
signal-to-noise and equally clean signals are, in fact, predicted to come from
nearby rich galaxy clusters. Most of the cluster emission is produced by small
subhaloes with masses less than that of the Sun. The large range of mass scales
covered by our two sets of simulations allows us to deduce a physically
motivated extrapolation to these small (and unresolved) masses. Since tidal
effects destroy subhaloes in the dense inner regions of haloes, most cluster
emission is then predicted to come from large radii, implying that the nearest
and brightest systems should be much more extended than Fermi's angular
resolution limit. The most promising targets for detection are clusters such as
Coma and Fornax, but detection algorithms must be tuned to the predicted
profile of the emission if they are to maximize the chance of finding this weak
signal.
[80]
oai:arXiv.org:1010.6197 [pdf] - 1041622
A Halo Expansion (HEX) Technique for Approximating Simulated Dark Matter
Haloes
Submitted: 2010-10-29, last modified: 2011-07-08
We apply a basis function expansion method to create a time-evolving
density/potential approximation of the late growth of simulated N-body dark
matter haloes. We demonstrate how the potential of a halo from the Aquarius
Project can be accurately represented by a small number of basis functions, and
show that the halo expansion (HEX) method provides a way to replay simulations.
We explore the level of accuracy of the technique as well as some of its
limitations. We find that the number of terms included in the expansion must be
large enough to resolve the large-scale distribution and shape of the halo but,
beyond this, additional terms result in little further improvement. Particle
and subhalo orbits can be integrated in this realistic, time-varying halo
potential approximation, at much lower cost than the original simulation, with
high fidelity for many individual orbits, and a good match to the distributions
of orbital energy and angular momentum. Statistically, the evolution of
structural subhalo properties, such as mass, half-mass radius and
characteristic circular velocity, are very well reproduced in the halo
expansion approximation over several gigayears. We demonstrate an application
of the technique by following the evolution of an orbiting subhalo at much
higher resolution than can be achieved in the original simulation. Our method
represents a significant improvement over commonly used techniques based on
static analytical descriptions of the halo potential.
[81]
oai:arXiv.org:1008.0484 [pdf] - 1034065
The Link between Galactic Satellite Orbits and Subhalo Accretion
Submitted: 2010-08-03, last modified: 2011-04-18
We calculate the orbital angular momentum of dark matter subhaloes in the
Aquarius simulations of cold dark matter galactic haloes. We calculate the
orientation of their angular momentum relative to that of the spin vector of
their host halo and find a variety of different configurations. All six
Aquarius haloes contain statistically significant populations of subhalo orbits
that are aligned with the main halo spin. All haloes posses a population of
subhaloes that rotates in the same direction as the main halo and three of them
possess, in addition, a population that rotates in the opposite direction.
These configurations arise from the filamentary accretion of subhaloes.
Quasi-planar distributions of coherently rotating satellites, such as those
inferred in the Milky Way and other galaxies, arise naturally in simulations of
a $\Lambda$CDM universe.
[82]
oai:arXiv.org:1102.0002 [pdf] - 958301
The Density and Pseudo-Phase-Space Density Profiles of CDM halos
Submitted: 2011-01-31
Cosmological N-body simulations indicate that the spherically-averaged
density profiles of cold dark matter halos are accurately described by Einasto
profiles, where the logarithmic slope is a power-law of adjustable exponent,
\gamma =dln\rho /dlnr ~ r^\alpha $. The pseudo-phase-space density (PPSD)
profiles of CDM halos also show remarkable regularity, and are well
approximated by simple power laws, Q(r)=\rho /\sigma ^3 ~ r^-\chi . We show
that this is expected from dynamical equilibrium considerations, since Jeans'
equations predict that the pseudo-phase-space density profiles of Einasto halos
should resemble power laws over a wide range of radii. For the values of \alpha
typical of CDM halos, the inner Q profiles of equilibrium halos deviate
significantly from a power law only very close to the center, and simulations
of extremely high-resolution would be needed to detect such deviations
unambiguously. We use an ensemble of halos drawn from the Millennium-II
simulation to study which of these two alternatives describe best the mass
profile of CDM halos. Our analysis indicates that at the resolution of the best
available simulations, both Einasto and power-law PPSD profiles (with
adjustable exponents \alpha and \chi, respectively) provide equally acceptable
fits to the simulations. A full account of the structure of CDM halos requires
understanding how the shape parameters that characterize departures from
self-similarity, like \alpha or \chi, are determined by evolutionary history,
environment or initial conditions.
[83]
oai:arXiv.org:1002.4539 [pdf] - 1025360
The evolution of galaxy cluster X-ray scaling relations
Submitted: 2010-02-24, last modified: 2010-12-15
We use numerical simulations to investigate, for the first time, the joint
effect of feedback from supernovae (SNe) and active galactic nuclei (AGN) on
the evolution of galaxy cluster X-ray scaling relations. Our simulations are
drawn from the Millennium Gas Project and are some of the largest
hydrodynamical N-body simulations ever carried out. Feedback is implemented
using a hybrid scheme, where the energy input into intracluster gas by SNe and
AGN is taken from a semi-analytic model of galaxy formation. This ensures that
the source of feedback is a population of galaxies that closely resembles that
found in the real universe. We show that our feedback model is capable of
reproducing observed local X-ray scaling laws, at least for non-cool core
clusters, but that almost identical results can be obtained with a simplistic
preheating model. However, we demonstrate that the two models predict opposing
evolutionary behaviour. We have examined whether the evolution predicted by our
feedback model is compatible with observations of high-redshift clusters.
Broadly speaking, we find that the data seems to favour the feedback model for
z<0.5, and the preheating model at higher redshift. However, a statistically
meaningful comparison with observations is impossible, because the large
samples of high-redshift clusters currently available are prone to strong
selection biases. As the observational picture becomes clearer in the near
future, it should be possible to place tight constraints on the evolution of
the scaling laws, providing us with an invaluable probe of the physical
processes operating in galaxy clusters.
[84]
oai:arXiv.org:1003.0671 [pdf] - 1025490
What is the (Dark) Matter with Dwarf Galaxies?
Submitted: 2010-03-03, last modified: 2010-12-07
We present cosmological hydrodynamical simulations of the formation of dwarf
galaxies in a representative sample of haloes extracted from the Millennium-II
Simulation. Our six haloes have a z = 0 mass of ~10^10 solar masses and show
different mass assembly histories which are reflected in different star
formation histories. We find final stellar masses in the range 5 x 10^7 - 10^8
solar masses, consistent with other published simulations of galaxy formation
in similar mass haloes. Our final objects have structures and stellar
populations consistent with dwarf elliptical and dwarf irregular galaxies.
However, in a Lambda CDM universe, 10^10 solar mass haloes must typically
contain galaxies with much lower stellar mass than our simulated objects if
they are to match observed galaxy abundances. The dwarf galaxies formed in our
own and all other current hydrodynamical simulations are more than an order of
magnitude more luminous than expected for haloes of this mass. We discuss the
significance and possible implications of this result.
[85]
oai:arXiv.org:1008.5114 [pdf] - 1034563
Assembly History and Structure of Galactic Cold Dark Matter Halos
Submitted: 2010-08-30, last modified: 2010-09-02
We use the Aquarius simulation series to study the imprint of assembly
history on the structure of Galaxy-mass cold dark matter halos. Our results
confirm earlier work regarding the influence of mergers on the mass density
profile and the inside-out growth of halos. The inner regions that contain the
visible galaxies are stable since early times and are significantly affected
only by major mergers. Particles accreted diffusely or in minor mergers are
found predominantly in the outskirts of halos. Our analysis reveals trends that
run counter to current perceptions of hierarchical halo assembly. For example,
major mergers (i.e. those with progenitor mass ratios greater than 1:10)
contribute little to the total mass growth of a halo, on average less than 20
per cent for our six Aquarius halos. The bulk is contributed roughly equally by
minor mergers and by "diffuse" material which is not resolved into individual
objects. This is consistent with modeling based on excursion-set theory which
suggests that about half of this diffuse material should not be part of a halo
of any scale. Interestingly, the simulations themselves suggest that a
significantly fraction is not truly diffuse, since it was ejected from earlier
halos by mergers prior to their joining the main system. The Aquarius
simulations resolve halos to much lower mass scales than are expected to retain
gas or form stars. These results thus confirm that most of the baryons from
which visible galaxies form are accreted diffusely, rather than through
mergers, and they suggest that only relatively rare major mergers will affect
galaxy structure at later times.
[86]
oai:arXiv.org:1006.2882 [pdf] - 295220
The statistics of the subhalo abundance of dark matter haloes
Submitted: 2010-06-15, last modified: 2010-08-31
We study the population statistics of the surviving subhaloes of LCDM dark
matter haloes using a set of very high resolution N-body simulations. These
include both simulations of representative regions of the Universe and
ultra-high resolution resimulations of individual dark matter haloes. We find
that more massive haloes tend to have a larger mass fraction in subhaloes. For
example, cluster size haloes typically have 7.5 percent of their mass within
R200 in substructures of fractional mass larger than 1e-5, which is 25 percent
higher than galactic haloes. There is, however, a large variance in the subhalo
mass fraction from halo to halo, whereas the subhalo abundance shows much
higher regularity. For dark matter haloes of fixed mass, the subhalo abundance
decreases by 30 percent between redshift 2 and 0. The subhalo abundance
function correlates with the host halo concentration parameter and formation
redshift. However, the intrinsic scatter is not significantly reduced for
narrow ranges of concentration parameter or formation redshift, showing that
they are not the dominant parameters that determine the subhalo abundance in a
halo.
[87]
oai:arXiv.org:0911.4484 [pdf] - 204072
There's no place like home? Statistics of Milky Way-mass dark matter
halos
Submitted: 2009-11-24, last modified: 2010-05-18
We present an analysis of the distribution of structural properties for Milky
Way-mass halos in the Millennium-II Simulation (MS-II). This simulation of
structure formation within the standard LCDM cosmology contains thousands of
Milky Way-mass halos and has sufficient resolution to properly resolve many
subhalos per host. It thus provides a major improvement in the statistical
power available to explore the distribution of internal structure for halos of
this mass. In addition, the MS-II contains lower resolution versions of the
Aquarius Project halos, allowing us to compare our results to simulations of
six halos at a much higher resolution. We study the distributions of mass
assembly histories, of subhalo mass functions and accretion times, and of
merger and stripping histories for subhalos capable of impacting disks at the
centers of halos. We show that subhalo abundances are not well-described by
Poisson statistics at low mass, but rather are dominated by intrinsic scatter.
Using the masses of subhalos at infall and the abundance-matching assumption,
there is less than a 10% chance that a Milky Way halo with M_vir =10^12 M_sun
will host two galaxies as bright as the Magellanic Clouds. This probability
rises to ~25% for a halo with M_vir=2.5 x 10^12 M_sun. The statistics relevant
for disk heating are very sensitive to the mass range that is considered
relevant. Mergers with infall mass : redshift zero virial mass greater than
1:30 could well impact a central galactic disk and are a near inevitability
since z=2, whereas only half of all halos have had a merger with infall mass :
redshift zero virial mass greater than 1:10 over this same period.
[88]
oai:arXiv.org:0910.3211 [pdf] - 268061
Galactic stellar haloes in the CDM model
Cooper, A. P.;
Cole, S.;
Frenk, C. S.;
White, S. D. M.;
Helly, J.;
Benson, A. J.;
De Lucia, G.;
Helmi, A.;
Jenkins, A.;
Navarro, J. F.;
Springel, V.;
Wang, J.
Submitted: 2009-10-19, last modified: 2010-04-14
We present six simulations of galactic stellar haloes formed by the tidal
disruption of accreted dwarf galaxies in a fully cosmological setting. Our
model is based on the Aquarius project, a suite of high resolution N-body
simulations of individual dark matter haloes. We tag subsets of particles in
these simulations with stellar populations predicted by the Galform
semi-analytic model. Our method self-consistently tracks the dynamical
evolution and disruption of satellites from high redshift. The luminosity
function and structural properties of surviving satellites, which agree well
with observations, suggest that this technique is appropriate. We find that
accreted stellar haloes are assembled between 1<z<7 from less than 5
significant progenitors. These progenitors are old, metal-rich satellites with
stellar masses similar to the brightest Milky Way dwarf spheroidals. In
contrast to previous stellar halo simulations, we find that several of these
major contributors survive as self-bound systems to the present day. Both the
number of these significant progenitors and their infall times are inherently
stochastic. This results in great diversity among our stellar haloes, which
amplifies small differences between the formation histories of their dark halo
hosts. The masses and density/surface-brightness profiles of the stellar haloes
are consistent with expectations from the Milky Way and M31. Each halo has a
complex structure, consisting of well-mixed components, tidal streams, shells
and other subcomponents. This structure is not adequately described by smooth
models. We find one example of an accreted thick disk. Contrasts in age and
metallicity between halo stars and those in surviving satellites are in broad
agreement with recent observations. [Abridged]
[89]
oai:arXiv.org:0909.0265 [pdf] - 951061
The properties of satellite galaxies in simulations of galaxy formation
Submitted: 2009-09-01, last modified: 2010-04-11
We investigate the properties of satellite galaxies in cosmological
N-body/SPH simulations of galaxy formation in Milky Way-sized haloes. Because
of their shallow potential wells, satellite galaxies are very sensitive to
heating processes which affect their gas content. Their properties can
therefore be used to constrain the nature of feedback processes that regulate
galaxy formation. In our simulations, we assume that all the energy produced by
supernovae is used as kinetic energy to drive galactic winds. Several of our
simulations produce bright, disc-dominated galaxies. We find that wind models
in which the wind speed, v_w, is proportional to local velocity dispersion of
dark matter, sigma, (and thus the wind mass-loading, eta_w \propto sigma^{-2})
have episodic star formation histories, reproduce the observed satellite
luminosity function reasonably well (down to M_v=-7) and match the
luminosity-metallicity relation observed in the Local Group satellites. By
contrast, models that assume a constant wind speed overproduce faint satellites
and predict an incorrect luminosity-metallicity relation. Our simulations
therefore suggest that the feedback processes that operate on the scale of
satellite galaxies should generate galactic outflows whose intensity varies
inversely with the depth of the potential.
[90]
oai:arXiv.org:1003.0974 [pdf] - 1025533
Full-sky map of the ISW and Rees-Sciama effect from Gpc simulations
Submitted: 2010-03-04, last modified: 2010-03-16
We present a new method for constructing maps of the secondary temperature
fluctuations imprinted on the cosmic microwave background (CMB) radiation by
photons propagating through the evolving cosmic gravitational potential. Large
cosmological N-body simulations are used to calculate the complete non-linear
evolution of the peculiar gravitational potential. Tracing light rays back
through the past lightcone of a chosen observer accurately captures the
temperature perturbations generated by linear (the integrated Sachs-Wolfe or
ISW effect) and non-linear (the Rees-Sciama or RS effect) evolution. These
effects give rise to three kinds of non-linear features in the temperature
maps. (a) In overdense regions, converging flows of matter induce cold spots of
order 100 Mpc in extent which can dominate over the ISW effect at high
redshift, and are surrounded by hot rings. (b) In underdense regions, the RS
effect enhances ISW cold spots which can be surrounded by weak hot rings. (c)
Transverse motions of large lumps of matter produce characteristic dipole
features, consisting of adjacent hot and cold spots separated by a few tens of
Megaparsecs. These non-linear features are not easily detectable; they modulate
the ISW sky maps at about the 10 percent level. The RS effect causes the
angular power spectrum to deviate from linear theory at l~50 and generates
non-Gaussianity, skewing the one-point distribution function to negative
values. Cold spots of similar angular size, but much smaller amplitude than the
CMB cold spot reported by Cruz et al. are produced. Joint analysis of our maps
and the corresponding galaxy distribution may enable techniques to be developed
to detect these non-linear, non-Gaussian features. Our maps are available at
http://astro.dur.ac.uk/~cai/ISW
[91]
oai:arXiv.org:0906.2785 [pdf] - 1434123
The angular momentum of cold dark matter haloes with and without baryons
Submitted: 2009-06-16, last modified: 2010-01-18
We investigate the magnitude and internal alignment of the angular momentum
of cold dark matter haloes in simulations with and without baryons. We analyse
the cumulative angular momentum profiles of hundreds of thousands of haloes in
the Millennium simulation and in a smaller, but higher resolution, simulation,
in total spanning 5 orders of magnitude in mass. For haloes of a given mass,
the median specific angular momentum increases as j(<r) proportional to r. The
direction of the vector varies considerably with radius: the median angle
between the inner (< 0.25 Rvir) and total (< Rvir) angular momentum vectors is
about 25degr. To investigate how baryons affect halo spin, we use another high
resolution simulation, which includes gas cooling, star formation and feedback.
This simulation produces a sample of galaxies with a realistic distribution of
disc-to-total ratios. The formation of the galaxy spins up the dark matter
within 0.1 Rvir such that the specific halo angular momentum increases by
approx 50% in the median. The dark matter angular momentum becomes better
aligned, but there remains a broad distribution of (mis-)alignments between the
halo and the central galaxy, with a median angle between their angular momenta
of ~ 30degr. Galaxies have a range of orientations relative to the shape of the
halo: half of them have their minor axes misaligned by more than 45degr.
Finally, we align a sample of haloes according to the orientation of their
galaxies and stack the projected mass distributions. Although the individual
haloes are aspherical, galaxy--halo misalignments produce a stacked mass
distribution that is indistinguishable from circular. If the misalignment found
in our simulations is realistic, it will be extremely difficult for weak
lensing studies to measure the shape of dark haloes using this technique.
[92]
oai:arXiv.org:1001.2310 [pdf] - 1024597
Secondary Infall and the Pseudo-Phase-Space Density Profiles of Cold
Dark Matter Halos
Submitted: 2010-01-13
We use N-body simulations to investigate the radial dependence of the density
and velocity dispersion in cold dark matter (CDM) halos. In particular, we
explore how closely Q rho/sigma^3, a surrogate measure of the phase-space
density, follows a power-law in radius. Our study extends earlier work by
considering, in addition to spherically-averaged profiles, local Q-estimates
for individual particles, Q_i; profiles based on the ellipsoidal radius
dictated by the triaxial structure of the halo, Q_i(r'); and by carefully
removing substructures in order to focus on the profile of the smooth halo,
Q^s. The resulting Q_i^s(r') profiles follow closely a power law near the
center, but show a clear upturn from this trend near the virial radius,
r_{200}. The location and magnitude of the deviations are in excellent
agreement with the predictions from Bertschinger's spherical secondary-infall
similarity solution. In this model, Q \propto r^{-1.875} in the inner,
virialized regions, but departures from a power-law occur near r_{200} because
of the proximity of this radius to the location of the first shell crossing -
the shock radius in the case of a collisional fluid. Particles there have not
yet fully virialized, and so Q departs from the inner power-law profile. Our
results imply that the power-law nature of $Q$ profiles only applies to the
inner regions and cannot be used to predict accurately the structure of CDM
halos beyond their characteristic scale radius.
[93]
oai:arXiv.org:0909.1593 [pdf] - 341093
The earliest stars and their relics in the Milky Way
Submitted: 2009-09-09, last modified: 2010-01-13
We have implemented a simple model to identify the likely sites of the first
stars and galaxies in the high-resolution simulations of the formation of
galactic dark matter halos of the Aquarius Project. The first star in a galaxy
like the Milky Way formed around redshift z=35; by z=10, the young galaxy
contained up to 30000 dark matter haloes capable of forming stars by molecular
hydrogen cooling. These minihaloes were strongly clustered and feedback may
have severely limited the actual number of Population III stars that formed. By
the present day, the remnants of the first stars would be strongly concentrated
to the centre of the main halo. If a second generation of long-lived stars
formed near the first (the first star relics), we would expect to find half of
them within 30kpc/h of the galactic centre and a significant fraction in
satellites where they may be singled out by their anomalous metallicity
patterns. The first halo in which gas could cool by atomic hydrogen line
radiation formed at z=25; by z=10, the number of such first galaxies had
increased to ~300. Feedback might have decreased the number of first galaxies
at the time they undergo significant star formation, but not the number that
survive to the present because near neighbours merge. Half of all the "first
galaxies" that form before z=10 merge with the main halo before z~3 and most
lose a significant fraction of their mass. However, today there should still be
more than 20 remnants orbiting within the central ~30kpc/h of the Milky Way.
These satellites have circular velocities of a few kilometers per second or
more, comparable to those of known Milky Way dwarfs. They are a promising
hunting ground for the remnants of the earliest epoch of star formation.
[94]
oai:arXiv.org:0910.0258 [pdf] - 1017799
Second-order Lagrangian perturbation theory initial conditions for
resimulations
Submitted: 2009-10-01, last modified: 2010-01-11
We describe and test a new method for creating initial conditions for
cosmological N-body dark matter simulations based on second-order Lagrangian
perturbation theory (2lpt). The method can be applied to multi-mass particle
distributions making it suitable for creating resimulation, or `zoom' initial
conditions. By testing against an analytic solution we show that the method
works well for a spherically symmetric perturbation with radial features
ranging over more than three orders of magnitude in linear scale and eleven
orders of magnitude in particle mass. We apply the method and repeat
resimulations of the rapid formation of a high mass halo at redshift fifty and
the formation of a Milky-Way mass dark matter halo at redshift zero. In both
cases we find that many properties of the final halo show a much smaller
sensitivity to the start redshift with the 2lpt initial conditions, than
simulations started from Zel'dovich initial conditions. For spherical overdense
regions structure formation is erroneously delayed in simulations starting from
Zel'dovich initial conditions, and we demonstrate for the case of the formation
the high redshift halo that this delay can be accounted for using the spherical
collapse model. In addition to being more accurate, 2lpt initial conditions
allow simulations to start later, saving computer time.
[95]
oai:arXiv.org:0903.3041 [pdf] - 177667
Resolving Cosmic Structure Formation with the Millennium-II Simulation
Submitted: 2009-03-18, last modified: 2009-09-18
We present the Millennium-II Simulation (MS-II), a very large N-body
simulation of dark matter evolution in the concordance LCDM cosmology. The
MS-II assumes the same cosmological parameters and uses the same particle
number and output data structure as the original Millennium Simulation (MS),
but was carried out in a periodic cube one-fifth the size (100 Mpc/h) with 5
times better spatial resolution (a Plummer equivalent softening of 1.0 kpc/h)
and with 125 times better mass resolution (a particle mass of 6.9 \times 10^6
Msun/h). By comparing results at MS and MS-II resolution, we demonstrate
excellent convergence in dark matter statistics such as the halo mass function,
the subhalo abundance distribution, the mass dependence of halo formation
times, the linear and nonlinear autocorrelations and power spectra, and halo
assembly bias. Together, the two simulations provide precise results for such
statistics over an unprecedented range of scales, from halos similar to those
hosting Local Group dwarf spheroidal galaxies to halos corresponding to the
richest galaxy clusters. The "Milky Way" halos of the Aquarius Project were
selected from a lower resolution version of the MS-II and were then resimulated
at much higher resolution. As a result, they are present in the MS-II along
with thousands of other similar mass halos. A comparison of their assembly
histories in the MS-II and in resimulations of 1000 times better resolution
shows detailed agreement over a factor of 100 in mass growth. We publicly
release halo catalogs and assembly trees for the MS-II in the same format
within the same archive as those already released for the MS.
[96]
oai:arXiv.org:0906.4350 [pdf] - 378157
Galaxies-Intergalactic Medium Interaction Calculation --I. Galaxy
formation as a function of large-scale environment
Crain, Robert A.;
Theuns, Tom;
Vecchia, Claudio Dalla;
Eke, Vincent R.;
Frenk, Carlos S.;
Jenkins, Adrian;
Kay, Scott T.;
Peacock, John A.;
Pearce, Frazer R.;
Schaye, Joop;
Springel, Volker;
Thomas, Peter A.;
White, Simon D. M.;
Wiersma, Robert P. C.
Submitted: 2009-06-23
[Abridged] We present the first results of hydrodynamical simulations that
follow the formation of galaxies to z=0 in spherical regions of radius ~20
Mpc/h drawn from the Millennium Simulation. The regions have overdensities that
deviate by (-2, -1, 0, +1, +2)sigma from the cosmic mean, where sigma is the
rms mass fluctuation on a scale of ~20Mpc/h at z=1.5. The simulations have mass
resolution of up to 10^6 Msun/h, cover the entire range of large-scale
environments and allow extrapolation of statistics to the entire 500 (Mpc/h)^3
Millennium volume. They include gas cooling, photoheating from an ionising
background, SNe feedback and winds, but no AGN. We find that the specific SFR
density at z <~ 10 varies systematically from region to region by up to an
order of magnitude, but the global value, averaged over all volumes, reproduces
observational data. Massive, compact galaxies, similar to those observed in the
GOODS fields, form in the overdense regions as early as z=6, but do not appear
in the underdense regions until z~3. These environmental variations are not
caused by a dependence of the star formation properties on environment, but
rather by a strong variation of the halo mass function from one environment to
another, with more massive haloes forming preferentially in the denser regions.
At all epochs, stars form most efficiently in haloes of circular velocity ~ 250
km/s. However, the star formation history exhibits a form of "downsizing" (even
in the absence of AGN): the stars comprising massive galaxies at z=0 have
mostly formed by z=1-2, whilst those comprising smaller galaxies typically form
at later times. However, additional feedback is required to limit star
formation in massive galaxies at late times.
[97]
oai:arXiv.org:0810.1522 [pdf] - 1000976
The Diversity and Similarity of Simulated Cold Dark Matter Halos
Submitted: 2008-10-08, last modified: 2009-06-11
We study the mass, velocity dispersion, and anisotropy profiles of
$\Lambda$CDM halos using a suite of N-body simulations of unprecedented
numerical resolution (the {\it Aquarius Project}). Our analysis confirms a
number of results claimed by earlier work, and clarifies a few issues where
conflicting claims may be found in the recent literature. The
spherically-averaged density profile becomes progressively shallower inwards
and, at the innermost resolved radius, the logarithmic slope is $\gamma \equiv
-$d$\ln\rho/$d$\ln r \simlt 1$. Asymptotic inner slopes as steep as the
recently claimed $\rho \propto r^{-1.2}$ are clearly ruled out. The radial
dependence of $\gamma$ is well approximated by a power-law, $\gamma \propto
r^{\alpha}$ (the Einasto profile). The shape parameter, $\alpha$, varies
slightly but significantly from halo to halo, implying that the mass profiles
of $\Lambda$CDM halos are not strictly universal: different halos cannot, in
general, be rescaled to look identical. Departures from similarity are also
seen in velocity dispersion profiles and correlate with those in density
profiles so as to preserve a power-law form for the spherically averaged
pseudo-phase-space density, $\rho/\sigma^3\propto r^{-1.875}$. Our conclusions
are reliable down to radii below 0.4% of the virial radius, providing
well-defined predictions for halo structure when baryonic effects are
neglected, and thus an instructive theoretical template against which the
modifications induced by the baryonic components of real galaxies can be
judged.
[98]
oai:arXiv.org:0903.4559 [pdf] - 1001772
Effects of Dark Matter Substructures on Gravitational Lensing: Results
from the Aquarius Simulations
Xu, D. D.;
Mao, S.;
Wang, J.;
Springel, V.;
Gao, L.;
White, S. D. M.;
Frenk, C. S.;
Jenkins, A.;
Li, G.;
Navarro, J. F.
Submitted: 2009-03-26, last modified: 2009-06-10
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.
[99]
oai:arXiv.org:0905.1696 [pdf] - 1002183
How common is the Milky Way - satellite system alignment?
Submitted: 2009-05-11
The highly flattened distribution of satellite galaxies in the Milky Way
presents a number of puzzles. Firstly, its polar alignment stands out from the
planar alignments commonly found in other galaxies. Secondly, recent proper
motion measurements reveal that the orbital angular momentum of at least 3, and
possibly as many as 8, of the Milky Ways satellites point (within 30 degrees)
along the axis of their flattened configuration, suggesting some form of
coherent motion. In this paper we use a high resolution cosmological simulation
to investigate whether this pattern conflicts with the expectations of the cold
dark matter model of structure formation. We find that this seemingly unlikely
set up occurs often: approximately 35% of the time we find systems in which the
angular momentum of 3 individual satellites point along, or close to, the short
axis of the satellite distribution. In addition, in 30% of the systems we find
that the net angular momentum of the 6 best aligned satellites lies within 35
degrees of the short axis of the satellite distribution, as observed for the
Milky Way.
[100]
oai:arXiv.org:0809.4488 [pdf] - 1000953
Towards accurate modelling of the ISW effect, the non-linear
contribution
Submitted: 2008-09-25, last modified: 2009-03-13
In a universe with a cosmological constant, the large-scale gravitational
potential varies in time and this is, in principle, observable. Using an N-body
simulation of a $\Lambda$CDM universe, we show that linear theory is not
sufficiently accurate to predict the power spectrum of the time derivative,
$\dot{\Phi}$, needed to compute the imprint of large-scale structure on the
cosmic microwave background (CMB). The linear part of the $\dot{\Phi}$ power
spectrum (the integrated Sachs-Wolfe effect or ISW) drops quickly as the
relative importance of $\Omega_{\Lambda}$ diminishes at high redshift, while
the non-linear part (the Rees-Sciama effect or RS) evolves more slowly with
redshift. Therefore, the deviation of the total power spectrum from linear
theory occurs at larger scales at higher redshifts. The deviation occurs at
$k\sim 0.1 $ $h$ Mpc$^{-1}$ at $z=0$. The cross-correlation power spectrum of
the density $\delta$ with $\dot{\Phi}$ behaves differently to the power
spectrum of $\dot{\Phi}$. Firstly, the deviation from linear theory occurs at
smaller scales ($k\sim 1 $ $h$ Mpc$^{-1}$ at $z=0$). Secondly, the correlation
becomes negative when the non-linear effect dominates. For the
cross-correlation power spectrum of galaxy samples with the CMB, the non-linear
effect becomes significant at $l\sim 500$ and rapidly makes the cross power
spectrum negative. For high redshift samples, the cross-correlation is expected
to be suppressed by $5-10%$ on arcminute scales. The RS effect makes a
negligible contribution to the large-scale ISW cross-correlation measurement.
However, on arc-minute scales it will contaminate the expected
cross-correlation signal induced by the Sunyaev-Zel'dovich effect.
[101]
oai:arXiv.org:0810.2300 [pdf] - 17392
Mock galaxy redshift catalogues from simulations: implications for
Pan-STARRS1
Submitted: 2008-10-13, last modified: 2008-12-22
We describe a method for constructing mock galaxy catalogues which are well
suited for use in conjunction with large photometric surveys. We use the
semi-analytic galaxy formation model of Bower et al. implemented in the
Millennium simulation. We apply our method to the specific case of the surveys
soon to commence with PS1, the first of 4 telescopes planned for the Pan-STARRS
system. PS1 has 5 photometric bands (grizy), and will carry out an all-sky 3pi
survey and a medium deep survey (MDS) over 84 sq.deg. We calculate the expected
magnitude limits for extended sources in the two surveys. We find that, after 3
years, the 3pi survey will have detected over 10^8 galaxies in all 5 bands, 10
million of which will lie at redshift z>0.9, while the MDS will have detected
over 10^7 galaxies with 0.5 million lying at z>2. These numbers at least double
if detection in the shallowest band, y is not required. We then evaluate the
accuracy of photometric redshifts estimated using an off-the-shelf photo-z
code. With the grizy bands alone it is possible to achieve an accuracy in the
3pi survey of Delta z/(1+z)~0.06 for 0.25<z<0.8, which could be reduced by
about 15% using near infrared photometry from the UKIDDS survey, but would
increase by about 25% for the deeper sample without the y band photometry. For
the MDS an accuracy of Delta z/(1+z)~0.05 is achievable for 0.02<z<1.5 using
grizy. A dramatic improvement in accuracy is possible by selecting only red
galaxies. In this case, Delta z/(1+z)~0.02-0.04 is achievable for ~100 million
galaxies at 0.4<z<1.1 in the 3pi survey and for 30 million galaxies in the MDS
at 0.4<z<2. We investigate the effect of using photo-z in the estimate of the
baryonic acoustic oscillation scale. We find that PS1 will achieve a similar
accuracy in this estimate as a spectroscopic survey of 20 million galaxies.
[102]
oai:arXiv.org:0804.0004 [pdf] - 11363
The clustering of the first galaxy halos
Submitted: 2008-04-01, last modified: 2008-12-10
We explore the clustering properties of high redshift dark matter halos,
focusing on halos massive enough to host early generations of stars or galaxies
at redshift 10 and greater. Halos are extracted from an array of dark matter
simulations able to resolve down to the "mini-halo" mass scale at redshifts as
high as 30, thus encompassing the expected full mass range of halos capable of
hosting luminous objects and sources of reionization. Halo clustering on
large-scales agrees with the Sheth, Mo & Tormen halo bias relation within all
our simulations, greatly extending the regime where large-scale clustering is
confirmed to be "universal" at the 10-20% level (which means, for example, that
3sigma halos of cluster mass at z=0 have the same large-scale bias with respect
to the mass distribution as 3sigma halos of galaxy mass at z=10). However, on
small-scales, the clustering of our massive halos (> ~10^9 Msun/h) at these
high redshifts is stronger than expected from comparisons with small-scale halo
clustering extrapolated from lower redshifts. This implies "non-universality"
in the scale-dependence of halo clustering, at least for the commonly used
parameterizations of the scale-dependence of bias that we consider. We provide
a fit for the scale-dependence of bias in our results. This study provides a
basis for using extraordinarily high redshift galaxies (redshift ~10) as a
probe of cosmology and galaxy formation at its earliest stages. We show also
that mass and halo kinematics are strongly affected by finite simulation
volumes. This suggests the potential for adverse affects on gas dynamics in
hydrodynamic simulations of limited volumes, such as is typical in simulations
of the formation of the "first stars", though further study is warranted.
[103]
oai:arXiv.org:0812.0362 [pdf] - 19059
Phase-space structure in the local dark matter distribution and its
signature in direct detection experiments
Submitted: 2008-12-02
We study predictions for dark matter phase-space structure near the Sun based
on high-resolution simulations of six galaxy halos taken from the Aquarius
Project. The local DM density distribution is predicted to be remarkably
smooth; the density at the Sun differs from the mean over a best-fit
ellipsoidal equidensity contour by less than 15% at the 99.9% confidence level.
The local velocity distribution is also very smooth, but it differs
systematically from a (multivariate) Gaussian distribution. This is not due to
the presence of individual clumps or streams, but to broad features in the
velocity modulus and energy distributions that are stable both in space and
time and reflect the detailed assembly history of each halo. These features
have a significant impact on the signals predicted for WIMP and axion searches.
For example, WIMP recoil rates can deviate by ~10% from those expected from the
best-fit multivariate Gaussian models. The axion spectra in our simulations
typically peak at lower frequencies than in the case of multivariate Gaussian
velocity distributions. Also in this case, the spectra show significant
imprints of the formation of the halo. This implies that once direct DM
detection has become routine, features in the detector signal will allow us to
study the dark matter assembly history of the Milky Way. A new field, "dark
matter astronomy", will then emerge.
[104]
oai:arXiv.org:0801.1127 [pdf] - 8758
The Unorthodox Orbits of Substructure Halos
Submitted: 2008-01-08, last modified: 2008-11-18
(Abridged) We use cosmological N-body simulations to study the properties of
substructure halos in galaxy-sized dark matter halos. We extend prior work on
the subject by considering the whole population of subhalos physically
associated with the main system. These are defined as subhalos that have at
some time in the past been within the virial radius of the halo's main
progenitor and that have survived as self-bound entities to $z=0$. We find that
this population extends beyond {\it three times} the virial radius, and
contains objects on extreme orbits. We trace the origin of these unorthodox
orbits to the tidal dissociation of bound groups of subhalos, which results in
the ejection of some subhalos along tidal streams. Ejected subhalos are
primarily low-mass systems, leading to mass-dependent biases in their spatial
distribution and kinematics: the lower the subhalo mass at accretion time, the
less centrally concentrated and kinematically hotter their descendant
population. The bias is strongest amongst the most massive subhalos, but
disappears at the low-mass end. Our findings imply that subhalos identified
within the virial radius represent an incomplete census of the substructure
physically related to a halo: only about {\it one half} of all associated
subhalos are found today within the virial radius of a halo. These results may
explain the age dependence of the clustering of low-mass halos, and has
implications for (i) the interpretation of the structural parameters and
assembly histories of halos neighboring massive systems; (ii) the existence of
low-mass dynamical outliers in the Local Group; and (iii) the presence of
evidence for evolutionary effects well outside the traditional virial boundary
of a galaxy system.
[105]
oai:arXiv.org:0809.0894 [pdf] - 15995
A blueprint for detecting supersymmetric dark matter in the Galactic
halo
Submitted: 2008-09-05
Dark matter is the dominant form of matter in the universe, but its nature is
unknown. It is plausibly an elementary particle, perhaps the lightest
supersymmetric partner of known particle species. In this case, annihilation of
dark matter in the halo of the Milky Way should produce gamma-rays at a level
which may soon be observable. Previous work has argued that the annihilation
signal will be dominated by emission from very small clumps (perhaps smaller
even than the Earth) which would be most easily detected where they cluster
together in the dark matter halos of dwarf satellite galaxies. Here we show,
using the largest ever simulation of the formation of a galactic halo, that
such small-scale structure will, in fact, have a negligible impact on dark
matter detectability. Rather, the dominant and likely most easily detectable
signal will be produced by diffuse dark matter in the main halo of the Milky
Way. If the main halo is strongly detected, then small dark matter clumps
should also be visible, but may well contain no stars, thereby confirming a key
prediction of the Cold Dark Matter (CDM) model.
[106]
oai:arXiv.org:0809.0898 [pdf] - 142165
The Aquarius Project: the subhalos of galactic halos
Submitted: 2008-09-05
We have performed the largest ever particle simulation of a Milky Way-sized
dark matter halo, and present the most comprehensive convergence study for an
individual dark matter halo carried out thus far. We have also simulated a
sample of 6 ultra-highly resolved Milky-way sized halos, allowing us to
estimate the halo-to-halo scatter in substructure statistics. In our largest
simulation, we resolve nearly 300,000 gravitationally bound subhalos within the
virialized region of the halo. Simulations of the same object differing in mass
resolution by factors up to 1800 accurately reproduce the largest subhalos with
the same mass, maximum circular velocity and position, and yield good
convergence for the abundance and internal properties of dark matter
substructures. We detect up to four generations of subhalos within subhalos,
but contrary to recent claims, we find less substructure in subhalos than in
the main halo when regions of equal mean overdensity are compared. The overall
substructure mass fraction is much lower in subhalos than in the main halo.
Extrapolating the main halo's subhalo mass spectrum down to an Earth mass, we
predict the mass fraction in substructure to be well below 3% within 100 kpc,
and to be below 0.1% within the Solar Circle. The inner density profiles of
subhalos show no sign of converging to a fixed asymptotic slope and are well
fit by gently curving profiles of Einasto form. The mean concentrations of
isolated halos are accurately described by the fitting formula of Neto et al.
down to maximum circular velocities of 1.5 km/s, an extrapolation over some 5
orders of magnitude in mass. However, at equal maximum circular velocity,
subhalos are more concentrated than field halos, with a characteristic density
that is typically ~2.6 times larger and increases towards the halo centre.
[107]
oai:arXiv.org:0711.0746 [pdf] - 6731
The redshift dependence of the structure of massive LCDM halos
Submitted: 2007-11-06, last modified: 2008-04-02
(Abridged) We use two very large cosmological simulations to study how the
density profiles of relaxed LCDM dark halos depend on redshift and on halo
mass. We confirm that these profiles deviate slightly but systematically from
the NFW form and are better approximated by the empirical formula, $d\log
\rho/d\log r \propto r^{\alpha}$, first used by Einasto to fit star counts in
the Milky Way. The best-fit value of the additional shape parameter, $\alpha$,
increases gradually with mass, from alpha ~ 0.16 for present-day galaxy halos
to alpha~ 0.3 for the rarest and most massive clusters. Halo concentrations
depend only weakly on mass at z=0, and this dependence weakens further at
earlier times. At z ~3 the average concentration of relaxed halos does not vary
appreciably over the mass range accessible to our simulations (M >=3e11Msun).
Furthermore, in our biggest simulation, the average concentration of the most
massive}, relaxed halos is constant at c200 ~3.5 to 4 for 0<=z<= 3. These
results have important implications for currently planned surveys of distant
clusters.
[108]
oai:arXiv.org:0709.3470 [pdf] - 314893
Constraints on sigma_8 from galaxy clustering in N-body simulations and
semi-analytic models
Submitted: 2007-09-21
We generate mock galaxy catalogues for a grid of different cosmologies, using
rescaled N-body simulations in tandem with a semi-analytic model run using
consistent parameters. Because we predict the galaxy bias, rather than fitting
it as a nuisance parameter, we obtain an almost pure constraint on sigma_8 by
comparing the projected two-point correlation function we obtain to that from
the SDSS. A systematic error arises because different semi-analytic modelling
assumptions allow us to fit the r-band luminosity function equally well.
Combining our estimate of the error from this source with the statistical
error, we find sigma_8=0.97 +/- 0.06. We obtain consistent results if we use
galaxy samples with a different magnitude threshold, or if we select galaxies
by b_J-band rather than r-band luminosity and compare to data from the 2dFGRS.
Our estimate for sigma_8 is higher than that obtained for other analyses of
galaxy data alone, and we attempt to find the source of this difference. We
note that in any case, galaxy clustering data provide a very stringent
constraint on galaxy formation models.
[109]
oai:arXiv.org:astro-ph/0702241 [pdf] - 89295
Virial Scaling of Massive Dark Matter Halos: Why Clusters Prefer a High
Normalization Cosmology
Evrard, A. E.;
Bialek, J.;
Busha, M.;
White, M.;
Habib, S.;
Heitmann, K.;
Warren, M.;
Rasia, E.;
Tormen, G.;
Moscardini, L.;
Power, C.;
Jenkins, A. R.;
Gao, L.;
Frenk, C. S.;
Springel, V.;
White, S. D. M.;
Diemand, J.
Submitted: 2007-02-08, last modified: 2007-08-29
We present a precise estimate of the bulk virial scaling relation of halos
formed via hierarchical clustering in an ensemble of simulated cold dark matter
cosmologies. The result is insensitive to cosmological parameters, the presence
of a trace, dissipationless gas component, and numerical resolution down to a
limit of ~1000 particles. The dark matter velocity dispersion scales with total
mass as log(sigma_{DM}(M,z)) = log(1082.9 +- 4.0 \kms) + (0.3361 +- 0.0026)
log(h(z)M_{200}/10^{15} Msun), with h(z) the dimensionless Hubble parameter. At
fixed mass, the velocity dispersion likelihood is nearly log-normal, with
scatter sigma_{ln sigma} = 0.0426 +- 0.015, except for a tail to higher
dispersions containing 10% of the population that are merger transients. We
combine this relation with the halo mass function in LCDM models, and show that
a low normalization condition, S_8 = sigma_8(Omega_m/0.3)^{0.35} \= 0.69,
favored by recent WMAP and SDSS analysis requires that galaxy and gas specific
energies in rich clusters be 50% larger than that of the underlying dark
matter. Such large energetic biases are in conflict with the current generation
of direct simulations of cluster formation. A higher normalization, S_8 = 0.80,
alleviates this tension and implies that the hot gas fraction within r_{500} is
(0.71 +- 0.09) h_{70}^{-3/2} Omega_b/\Omega_m, a value consistent with recent
Sunyaev-Zel'dovich observations.
[110]
oai:arXiv.org:astro-ph/0509323 [pdf] - 270440
The effects of ellipticity and substructure on estimates of cluster
density profiles based on lensing and kinematics
Submitted: 2005-09-13, last modified: 2007-07-18
We address the question of how well the density profile of galaxy clusters
can be determined by combining strong lensing and velocity dispersion data. We
use cosmological dark matter simulations of clusters to test the reliability of
the method, producing mock catalogues of tangential and radial gravitational
arcs and simulating the radial velocity dispersion profile of the cluster
brightest central galaxy. The density profiles of the simulated clusters
closely follow the NFW form, but we find that the recovered values of the inner
slope are systematically underestimated, by about 0.4 in the mean, if the lens
is assumed to be axially symmetric. However, if the ellipticity and orientation
of the iso-contours of the cluster lensing potential are taken into account,
then the inner slopes can be recovered quite accurately for a significant
subset of the clusters whose central surface density profiles appear the most
regular. These have lensing potentials with ellipticities in the range
0.15-0.4. Further simulations projecting one cluster along many random
lines-of-sight show that, even for lower ellipticities, the central slopes are
underestimated by $\sim 10-35%$. These simulations closely mimic past
observations (see e.g. Sand et al., 2004), suggesting that existing estimates
of the central slopes may be biased towards low values. For the remaining
clusters, where the lensing potential is strongly perturbed by active merging
or by substructure, the correct determination of the inner slope requires a
more accurate model for the lens. When the halo profile is modelled by a
generalised NFW profile, we find that the inferred scale radius and
characteristic density, unlike the inner slope, are generally poorly
constrained, since there is a strong degeneracy between these two parameters.
[111]
oai:arXiv.org:astro-ph/0608190 [pdf] - 84067
Formation of z ~ 6 quasars from hierarchical galaxy mergers
Li, Yuexing;
Hernquist, Lars;
Robertson, Brant;
Cox, Thomas J.;
Hopkins, Philip F.;
Springel, Volker;
Gao, Liang;
Di Matteo, Tiziana;
Zentner, Andrew R.;
Jenkins, Adrian;
Yoshida, Naoki
Submitted: 2006-08-09, last modified: 2007-06-26
The discovery of luminous quasars at redshift z ~ 6 indicates the presence of
supermassive black holes (SMBHs) of mass ~10^9 Msun when the Universe was less
than one billion years old. This finding presents several challenges for
theoretical models. Here, we present the first multi-scale simulations that,
together with a self-regulated model for the SMBH growth, produce a luminous
quasar at z ~ 6.5 in the LCDM paradigm. We follow the hierarchical assembly
history of the most massive halo in a ~ 3 Gpc^3 volume, and find that this halo
of ~ 8x 10^{12} Msun forming at z ~ 6.5 after several major mergers is able to
reproduce a number of observed properties of SDSS J1148+5251, the most distant
quasar detected at z =6.42 (Fan et al. 2003). Moreover, the SMBHs grow through
gas accretion below the Eddington limit in a self-regulated manner owing to
feedback. We find that the progenitors experience significant star formation
(up to 10^4 Msun/yr) preceding the major quasar phase such that the stellar
mass of the quasar host reaches 10^{12} Msun at z ~ 6.5, consistent with
observations of significant metal enrichment in SDSS J1148+5251. Our results
provide a viable formation mechanism for z ~ 6 quasars in the standard LCDM
cosmology, and demonstrate a common, merger-driven origin for the rarest
quasars and the fundamental SMBH-host correlation in a hierarchical
Universe.(Abridged)
[112]
oai:arXiv.org:0706.2919 [pdf] - 398516
The statistics of LCDM Halo Concentrations
Submitted: 2007-06-20
We use the Millennium Simulation (MS) to study the statistics of LCDM halo
concentrations at z = 0. Our results confirm that the average halo
concentration declines monotonically with mass; a power-law fits well the
concentration-mass relation for over 3 decades in mass, up to the most massive
objects to form in a LCDM universe (~ 10^15 h^-1 Msol). This is in clear
disagreement with the predictions of the model proposed by Bullock et al. for
these rare objects, and agrees better with the original predictions of Navarro,
Frenk, & White. The large volume surveyed, together with the unprecedented
numerical resolution of the MS, allow us to estimate with confidence the
distribution of concentrations and, consequently, the abundance of systems with
unusual properties. About one in a hundred cluster haloes (M200 >~ 3x10^14 h^-1
Msol) have concentrations exceeding c200 = 7.5, a result that may be used to
interpret the likelihood of unusually strong massive gravitational lenses, such
as Abell 1689, in the LCDM cosmogony. A similar fraction (1 in 100) of
galaxy-sized haloes (M200 ~ 10^12 h^-1 Msol) have c200 < 4.5, an important
constraint on models that attempt to reconcile the rotation curves of low
surface-brightness galaxies by appealing to haloes of unexpectedly low
concentration. We find that halo concentrations are independent of spin once
haloes manifestly out of equilibrium are removed from the sample. Compared to
their relaxed brethren, the concentrations of out-of-equilibrium haloes tend to
be lower and to have more scatter, while their spins tend to be higher...
(continue)
[113]
oai:arXiv.org:astro-ph/0610174 [pdf] - 85598
The first generation of stars in LCDM cosmology
Submitted: 2006-10-05, last modified: 2007-04-04
We have performed a large set of high-resolution cosmological simulations
using smoothed particle hydrodynamics to study the formation of the first
luminous objects in the LCDM cosmology. We follow the collapse of primordial
gas clouds in eight early structures and document the scatter in the properties
of the first star-forming clouds. Our first objects span formation redshifts
from z~10 to z~50 and cover an order of magnitude in halo mas. We find that the
physical properties of the central star-forming clouds are very similar in all
of the simulated objects despite significant differences in formation redshift
and environment. The physical properties of the clouds have little correlation
with spin, mass, or assembly history of the host halo. The collapse of
protostellar objects at higher redshifts progresses much more rapidly due to
the higher densities, which accelerates the formation of molecular hydrogen,
enhances initial cooling and shortens the dynamical timescales. The mass of the
star-forming clouds cover a broad range, from a few hundred to a few thousand
solar masses, and exhibit various morphologies: some of have disk-like
structures nearly rotational supported; others form flattened spheroids; still
others form bars. All of them develop a single proto-stellar `seed' which does
not fragment into multiple objects up to the moment that the central gas
becomes optically thick to H2 cooling lines. At this time, the instantaneous
mass accretion rate onto the centre varies significantly from object to object,
with disk-like structuress have the smallest mass accretion rates. The
formation epoch and properties of the star-forming clouds are sensitive to the
values of cosmological parameters.
[114]
oai:arXiv.org:astro-ph/0608607 [pdf] - 84484
The spin and shape of dark matter haloes in the Millennium simulation of
a LambdaCDM universe
Submitted: 2006-08-28, last modified: 2006-12-20
We investigate the spins and shapes of over a million dark matter haloes
identified at z=0 in the Millennium simulation. Our sample spans halo masses
ranging from dwarf galaxies to rich galaxy clusters. The very large dynamic
range of this simulation enables the distribution of spins and shapes and their
variation with halo mass and environment to be characterised with unprecedented
precision. We compare results for haloes identified using different algorithms,
including a novel one based on the branches of the halo merger trees. We
investigate (and remove) biases in the estimate of angular momentum introduced
by both the algorithms themselves and by numerical effects. We find that for
this many objects, the traditional lognormal function is no longer an adequate
description of the distribution of the dimensionless spin parameter lambda, and
we provide a different function that gives a better fit. The variation of spin
with halo mass is weak but detectable, although the trend depends strongly on
the halo definition used. The haloes exhibit a range of shapes, with a
preference for prolateness over oblateness. More massive haloes tend to be less
spherical and more prolate. We find that the more spherical haloes have spin in
the median, and those closest to spherical have a spin independent of mass. The
most massive have a spin independent of shape. The majority of haloes have
their angular momentum vector aligned with their minor axis and perpendicular
to their major axis. We find that higher tend to be more clustered, with a
stronger effect for more massive haloes. (abridged)
[115]
oai:arXiv.org:astro-ph/0607150 [pdf] - 83353
The halo mass function from the dark ages through the present day
Submitted: 2006-07-07, last modified: 2006-11-20
We use an array of high-resolution N-body simulations to determine the mass
function of dark matter haloes at redshifts 10-30. We develop a new method for
compensating for the effects of finite simulation volume that allows us to find
an approximation to the true ``global'' mass function. By simulating a wide
range of volumes at different mass resolution, we calculate the abundance of
haloes of mass 10^{5-12} Msun/h. This enables us to predict accurately the
abundance of the haloes that host the sources that reionize the universe. In
particular, we focus on the small mass haloes (>~10^{5.5-6} Msun/h) likely to
harbour population III stars where gas cools by molecular hydrogen emission,
early galaxies in which baryons cool by atomic hydrogen emission at a virial
temperature of ~10^4K (10^{7.5-8} Msun/h), and massive galaxies that may be
observable at redshift ~10. When we combine our data with simulations that
include high mass halos at low redshift, we find that the best fit to the halo
mass function depends not only on linear overdensity, as is commonly assumed in
analytic models, but also upon the slope of the linear power spectrum at the
scale of the halo mass. The Press-Schechter model gives a poor fit to the halo
mass function in the simulations at all epochs; the Sheth-Tormen model gives a
better match, but still overpredicts the abundance of rare objects at all times
by up to 50%. Finally, we consider the consequences of the recently released
WMAP 3-year cosmological parameters. These lead to much less structure at high
redshift, reducing the number of z=10 ``mini-haloes'' by more than a factor of
two and the number of z=30 galaxy hosts by more than four orders of magnitude.
Code to generate our best-fit halo mass function may be downloaded from
http://icc.dur.ac.uk/Research/PublicDownloads/genmf_readme.html
[116]
oai:arXiv.org:astro-ph/0610602 [pdf] - 86026
The baryon fraction of LambdaCDM haloes
Submitted: 2006-10-19
We investigate the baryon fraction in dark matter haloes formed in
non-radiative gas-dynamical simulations of the LambdaCDM cosmogony. By
combining a realisation of the Millennium Simulation (Springel et al.) with a
simulation of a smaller volume focussing on dwarf haloes, our study spans five
decades in halo mass, from 10^10 Msun/h to 10^15 Msun/h. We find that the
baryon fraction within the halo virial radius is typically 90% of the cosmic
mean, with an rms scatter of 6%, independently of redshift and of halo mass
down to the smallest resolved haloes. Our results show that, contrary to the
proposal of Mo et al. (2005), pre-virialisation gravitational heating is unable
to prevent the collapse of gas within galactic and proto-galactic haloes, and
confirm the need for non-gravitational feedback in order to reduce the
efficiency of gas cooling and star formation in dwarf galaxy haloes.
Simulations including a simple photoheating model (where a gas temperature
floor of T_{floor} = 2x10^4 K is imposed from z=11) confirm earlier suggestions
that photoheating can only prevent the collapse of baryons in systems with
virial temperatures T_{200} < ~2.2 T_{floor} ~ 4.4x10^4 K (corresponding to a
virial mass of M_{200} ~ 10^10 Msun/h and a circular velocity of V_{200} ~ 35
km/s). Photoheating may thus help regulate the formation of dwarf spheroidals
and other galaxies at the extreme faint-end of the luminosity function, but it
cannot, on its own, reconcile the abundance of sub-L* galaxies with the vast
number of dwarf haloes expected in the LambdaCDM cosmogony. The lack of
evolution or mass dependence seen in the baryon fraction augurs well for X-ray
cluster studies that assume a universal and non-evolving baryon fraction to
place constraints on cosmological parameters.
[117]
oai:arXiv.org:astro-ph/0607237 [pdf] - 83440
Satellite Systems around Galaxies in Hydrodynamic Simulations
Submitted: 2006-07-11
We investigate the properties of satellite galaxies formed in N-body/SPH
simulations of galaxy formation in the $\Lambda$CDM cosmology. The simulations
include the main physical effects thought to be important in galaxy formation
and, in several cases, produce realistic spiral discs. In total, a sample of 9
galaxies of luminosity comparable to the Milky Way was obtained. At magnitudes
brighter than the resolution limit, $M_V=-12$, the luminosity function of the
satellite galaxies in the simulations is in excellent agreement with data for
the Local Group. The radial number density profile of the model satellites, as
well as their gas fractions also match observations very well. In agreement
with previous N-body studies, we find that the satellites tend to be
distributed in highly flattened configurations whose major axis is aligned with
the major axis of the (generally triaxial) dark halo. In 2 out of 3 systems
with sufficiently large satellite populations, the satellite system is nearly
perpendicular to the plane of the galactic disc, a configuration analogous to
that observed in the Milk Way. The discs themselves are perpendicular to the
minor axis of their host halos in the inner parts, and the correlation between
the orientation of the galaxy and the shape of the halo persists even out to
the virial radius. However, in one case the disc's minor axis ends up, at the
virial radius, perpendicular to the minor axis of the halo. The angular momenta
of the galaxies and their host halo tend to be well aligned.
[118]
oai:arXiv.org:astro-ph/0508046 [pdf] - 74935
The many lives of active galactic nuclei: cooling flows, black holes and
the luminosities and colours of galaxies
Croton, Darren J.;
Springel, Volker;
White, Simon D. M.;
De Lucia, G.;
Frenk, C. S.;
Gao, L.;
Jenkins, A.;
Kauffmann, G.;
Navarro, J. F.;
Yoshida, N.
Submitted: 2005-08-02, last modified: 2006-06-14
We simulate the growth of galaxies and their central supermassive black holes
by implementing a suite of semi-analytic models on the output of the Millennium
Run, a very large simulation of the concordance LCDM cosmogony. Our procedures
follow the detailed assembly history of each object and are able to track the
evolution of all galaxies more massive than the Small Magellanic Cloud
throughout a volume comparable to that of large modern redshift surveys. In
this first paper we supplement previous treatments of the growth and activity
of central black holes with a new model for `radio' feedback from those AGN
that lie at the centre of a quasistatic X-ray emitting atmosphere in a galaxy
group or cluster. We show that for energetically and observationally plausible
parameters such a model can simultaneously explain: (i) the low observed mass
drop-out rate in cooling flows; (ii) the exponential cut-off at the bright end
of the galaxy luminosity function; and (iii) the fact that the most massive
galaxies tend to be bulge-dominated systems in clusters and to contain
systematically older stars than lower mass galaxies. This success occurs
because static hot atmospheres form only in the most massive structures, and
radio feedback (in contrast, for example, to supernova or starburst feedback)
can suppress further cooling and star formation without itself requiring star
formation. We discuss possible physical models which might explain the
accretion rate scalings required for our phenomenological `radio mode' model to
be successful.
[119]
oai:arXiv.org:astro-ph/0510488 [pdf] - 419228
A marked correlation function analysis of halo formation times in the
Millennium Simulation
Submitted: 2005-10-17, last modified: 2006-03-29
We study the environmental dependence of the formation epoch of dark matter
haloes in the Millennium Simulation: a ten billion particle N-body simulation
of standard LambdaCDM cosmology. A sensitive test of this dependence - the
marked correlation function - reveals highly significant evidence that haloes
of a given mass form earlier in denser regions. We define a marked
cross-correlation function, which helps quantify how this effect depends upon
the choice of the halo population used to define the environment. The mean halo
formation redshift as a function of the local overdensity in dark matter is
also well determined, and we see an especially clear dependence for
galaxy-sized haloes. This contradicts one of the basic predictions of the
excursion set model of structure formation, even though we see that this theory
predicts other features of the distribution of halo formation epochs rather
well. It also invalidates an assumption usually employed in the popular halo,
or HOD, models of galaxy clustering, namely that the distribution of halo
properties is a function of halo mass but not of halo environment.
[120]
oai:arXiv.org:astro-ph/0602065 [pdf] - 79641
Erratum - The many lives of AGN: cooling flows, black holes and the
luminosities and colours of galaxies
Croton, Darren J.;
Springel, Volker;
White, Simon D. M.;
De Lucia, G.;
Frenk, C. S.;
Gao, L.;
Jenkins, A.;
Kauffmann, G.;
Navarro, J. F.;
Yoshida, N.
Submitted: 2006-02-03
In Figure 6 we inadvertently labeled the proxy circular velocity as the
virial velocity of the dark matter halo instead of what is actually plotted,
the maximum circular velocity of the dark matter halo. The maximum halo
circular velocity is a much better estimate of the disk V_c than is V_vir. This
confusion influenced the discussion of the Tully-Fisher relation in our paper.
In fact, Figure 6 demonstrates that it is possible to simultaneously reproduce
both the local Tully-Fisher relation and luminosity function using
semi-analytic techniques applied to the standard LCDM cosmology, thus
contradicting previous studies of this issue and our own discussion in Section
3.6.
[121]
oai:arXiv.org:astro-ph/0504038 [pdf] - 72135
The First Generation of Star-Forming Haloes
Submitted: 2005-04-04, last modified: 2005-12-14
We model gas cooling in high-resolution N-body simulations in order to
investigate the formation of the first generation of stars. We follow a region
of a LCDM universe especially selected to contain a rich cluster by the present
day. The properties of the dark haloes that form in these sub-solar mass
resolution simulations are presented in a companion paper by Gao et al. The
first gas clouds able to cool by molecular hydrogen line emission collapse at
extremely high redshift, z ~ 47, when the mass of the dark halo is 2.4 x 10^5
Msun/h. By z ~ 30, a substantial population of haloes are capable of undergoing
molecular hydrogen cooling although their ability to form stars is dependent on
the efficiency of feedback processes such as dissociating Lyman-Werner
radiation. The mass of the main halo grows extremely rapidly and, by z ~ 36,
its virial temperature has reached 10^4K, at which point gas cooling becomes
dominated by more effective atomic line processes. By z ~ 30, a small ``group''
of such potential galaxies will have formed unless prevented from doing so by
feedback processes. By this redshift, massive (>~ 100 Msun) population III
stars are able to ionise gas well beyond their own host halo and neighbouring
HII regions can percolate to form an ionized superbubble. Such patches would be
too widely separated to contribute significantly to reionisation at this time.
The large number density of early cooling haloes in the pre-reionised universe
raises the exciting prospect that this ultra-early generation of stars may be
observable as gamma-ray bursts or supernovae.
[122]
oai:arXiv.org:astro-ph/0503676 [pdf] - 72058
The effects of feedback on the morphology of galaxy discs
Submitted: 2005-03-30, last modified: 2005-08-16
We have performed hydrodynamic simulations of galaxy formation in a LCDM
universe. We have followed galaxy formation in a dark matter halo, chosen to
have a relatively quiet recent merger history, using different models for star
formation and feedback. In all cases, we have adopted a multi-phase description
of the interstellar medium and modelled star formation in quiescent and burst
modes. We have explored two triggers for starbursts: strong shocks and high gas
density, allowing for the possibility that stars in the burst may form with a
top-heavy initial mass function. We find that the final morphology of the
galaxy is extremely sensitive to the modelling of star formation and feedback.
Starting from identical initial conditions, galaxies spanning the wide range of
Hubble types, with B-band disc-to-total luminosity ratios ranging from 0.25 to
0.9, can form in the same dark matter halo. Models in which starbursts are
induced by high gas density (qualitatively similar to models in which feedback
is produced by AGN) generate energetic winds and result in galaxies with an
early-type morphology. Models in which the starbursts are induced by strong
shocks lead to extended discs. In this case, the feedback associated with the
bursts suppresses the collapse of baryons in small halos, helping to create a
reservoir of hot gas that is available for cooling after $z = 1$, following the
bulk of the dynamical activity that builds up the halo. This gas then cools to
form a large, young stellar disc.
[123]
oai:arXiv.org:astro-ph/0503003 [pdf] - 71385
Early structure in LCDM
Submitted: 2005-03-01, last modified: 2005-08-12
We use a novel technique to simulate the growth of the most massive
progenitor of a protocluster region from redshift z~80, when its mass was about
10 Msun until the present day. Our nested sequence of N-body resimulations
allows us to study in detail the structure both of the dark matter object
itself and of its environment. Our effective resolution is optimal at redshifts
of 49, 29, 12, 5 and 0 when the dominant object has mass 1.2e5, 5e7, 2e10, 3e12
and 8e14 Msun/hrespectively, and contains ~1e6 simulation particles within its
virial radius. Extended Press-Schechter theory correctly predicts both this
rapid growth and the substantial overabundance of massive haloes we find at
early times in regions surrounding the dominant object. Although the
large-scale structure in these regions differs dramatically from a scaled
version of its present-day counterpart, the internal structure of the dominant
object is remarkably similar. Molecular hydrogen cooling could start as early
as z~49 in this object, while cooling by atomic hydrogen becomes effective at
z~39. If the first stars formed in haloes with virial temperature ~2000, their
comoving abundance by z=49 would be similar to that of dwarf galaxies today,
while their comoving correlation length would be ~2.5Mpc/h.
[124]
oai:arXiv.org:astro-ph/0501174 [pdf] - 70286
The 2dF Galaxy Redshift Survey: Power-spectrum analysis of the final
dataset and cosmological implications
Cole, S.;
Percival, W. J.;
Peacock, J. A.;
Norberg, P.;
Baugh, C. M.;
Frenk, C. S.;
Baldry, I.;
Bland-Hawthorn, J.;
Bridges, T.;
Cannon, R.;
Colless, M.;
Collins, C.;
Couch, W.;
Cross, N. J. G.;
Dalton, G.;
Eke, V. R.;
De Propris, R.;
Driver, S. P.;
Efstathiou, G.;
Ellis, R. S.;
Glazebrook, K.;
Jackson, C.;
Jenkins, A.;
Lahav, O.;
Lewis, I.;
Lumsden, S.;
Maddox, S.;
Madgwick, D.;
Peterson, B. A.;
Sutherland, W.;
Taylor, K.
Submitted: 2005-01-11, last modified: 2005-08-05
We present a power spectrum analysis of the final 2dF Galaxy Redshift Survey,
employing a direct Fourier method. The sample used comprises 221,414 galaxies
with measured redshifts. We investigate in detail the modelling of the sample
selection. A new angular mask is derived, based on revisions to the photometric
calibration. The redshift selection function is determined by dividing the
survey according to rest-frame colour, and deducing a self-consistent treatment
of k-corrections and evolution for each population. The covariance matrix for
the power-spectrum estimates is determined using two different approaches to
the construction of mock surveys which are used to demonstrate that the input
cosmological model can be correctly recovered. We are confident that the 2dFGRS
power spectrum can be used to infer the matter content of the universe. On
large scales, our estimated power spectrum shows evidence for the `baryon
oscillations' that are predicted in CDM models. Fitting to a CDM model,
assuming a primordial $n_{s}=1$ spectrum, $h=0.72$ and negligible neutrino
mass, the preferred parameters are $\Omega_{M} h = 0.168 \pm 0.016$ and a
baryon fraction $\Omega_{b} /\Omega_{M} = 0.185\pm0.046$ (1$\sigma$ errors).
The value of $\Omega_{M} h$ is $1\sigma$ lower than the $0.20 \pm 0.03$ in our
2001 analysis of the partially complete 2dFGRS. This shift is largely due to
the signal from the newly-sampled regions of space, rather than the refinements
in the treatment of observational selection. This analysis therefore implies a
density significantly below the standard $\Omega_{M} =0.3$: in combination with
CMB data from WMAP, we infer $\Omega_{M} =0.231\pm 0.021$. (Abridged.)
[125]
oai:arXiv.org:astro-ph/0504456 [pdf] - 72553
Constraints on the dark energy equation of state from the imprint of
baryons on the power spectrum of clusters
Submitted: 2005-04-21, last modified: 2005-06-13
Acoustic oscillations in the baryon-photon fluid leave a signature in the
matter power spectrum. The overall shape of the spectrum and the wavelength of
the oscillations depend upon the sound horizon scale at recombination. Using
the $\Lambda$ cold dark matter Hubble Volume simulation, we show that the
imprint of baryons is visible in the power spectrum of cluster-mass dark matter
haloes, in spite of significant differences between the halo power spectrum and
the prediction of linear perturbation theory. A measurement of the sound
horizon scale can constrain the dark energy equation of state. We show that a
survey of clusters at intermediate redshift ($ z\sim1 $), like the
Sunyaev-Zeldovich survey proposed by the South Pole Telescope or a red sequence
photometric survey with VISTA, could potentially constrain the sound horizon
scale to an accuracy of $\sim 2%$, in turn fixing the ratio of the pressure of
the dark energy to its density ($w$) to better than $\sim 10%$. Our approach
does not require knowledge of the cluster mass, unlike those that depend upon
the abundance of clusters.
[126]
oai:arXiv.org:astro-ph/0504097 [pdf] - 72194
Simulating the joint evolution of quasars, galaxies and their
large-scale distribution
Springel, Volker;
White, Simon D. M.;
Jenkins, Adrian;
Frenk, Carlos S.;
Yoshida, Naoki;
Gao, Liang;
Navarro, Julio;
Thacker, Robert;
Croton, Darren;
Helly, John;
Peacock, John A.;
Cole, Shaun;
Thomas, Peter;
Couchman, Hugh;
Evrard, August;
Colberg, Joerg;
Pearce, Frazer
Submitted: 2005-04-05, last modified: 2005-04-06
The cold dark matter model has become the leading theoretical paradigm for
the formation of structure in the Universe. Together with the theory of cosmic
inflation, this model makes a clear prediction for the initial conditions for
structure formation and predicts that structures grow hierarchically through
gravitational instability. Testing this model requires that the precise
measurements delivered by galaxy surveys can be compared to robust and equally
precise theoretical calculations. Here we present a novel framework for the
quantitative physical interpretation of such surveys. This combines the largest
simulation of the growth of dark matter structure ever carried out with new
techniques for following the formation and evolution of the visible components.
We show that baryon-induced features in the initial conditions of the Universe
are reflected in distorted form in the low-redshift galaxy distribution, an
effect that can be used to constrain the nature of dark energy with next
generation surveys.
[127]
oai:arXiv.org:hep-th/0501223 [pdf] - 141555
Spontaneous Lorentz Violation and the Long-Range Gravitational
Preferred-Frame Effect
Submitted: 2005-01-28, last modified: 2005-03-22
Lorentz-violating operators involving Standard Model fields are tightly
constrained by experimental data. However, bounds are more model-independent
for Lorentz violation appearing in purely gravitational couplings. The
spontaneous breaking of Lorentz invariance by the vacuum expectation value of a
vector field selects a universal rest frame. This affects the propagation of
the graviton, leading to a modification of Newton's law of gravity. We compute
the size of the long-range preferred-frame effect in terms of the coefficients
of the two-derivative operators in the low-energy effective theory that
involves only the graviton and the Goldstone bosons.
[128]
oai:arXiv.org:astro-ph/0503400 [pdf] - 71782
The Distribution of Satellite Galaxies: The Great Pancake
Submitted: 2005-03-18
The 11 known satellite galaxies within 250 kpc of the Milky Way lie close to
a great circle on the sky. We use high resolution N-body simulations of
galactic dark matter halos to test if this remarkable property can be
understood within the context of the cold dark matter cosmology. We construct
halo merger trees from the simulations and use a semianalytic model to follow
the formation of satellite galaxies. We find that in all 6 of our simulations,
the 11 brightest satellites are indeed distributed along thin, disk-like
structures analogous to that traced by the Milky Way's satellites. This is in
sharp contrast to the overall distributions of dark matter in the halo and of
subhalos within it which, although triaxial, are not highly aspherical. We find
that the spatial distribution of satellites is significantly different from
that of the most massive subhalos but is similar to that of the subset of
subhalos that had the most massive progenitors at earlier times. The elongated
disk-like structure delineated by the satellites has its long axis aligned with
the major axis of the dark matter halo. We interpret our results as reflecting
the preferential infall of satellites along the spines of a few filaments of
the cosmic web.
[129]
oai:arXiv.org:astro-ph/0502515 [pdf] - 71306
A Universal Density Profile for Dark and Luminous Matter?
Submitted: 2005-02-24
We explore similarities in the luminosity distribution of early type galaxies
and the mass profiles of LCDM halos. The spatial structure of these systems may
be accurately described by a simple law where the logarithmic slope of the
projected density is a power law of radius; the Sersic law. We show that this
law provides a significantly better fit than a three-parameter generalization
of the NFW profile and derive the best-fitting Sersic parameters for a set of
high-resolution LCDM halos spanning a wide range in mass. The mean Sersic n
values are 3.0 for dwarf- and galaxy-sized halos and 2.4 for cluster-sized
halos, similar to the values that characterize luminous elliptical galaxies. We
discuss possible reasons why the same law should describe dark and luminous
systems that span a range of over seven decades in mass.
[130]
oai:arXiv.org:hep-th/0407174 [pdf] - 117403
Anthropic Distribution for Cosmological Constant and Primordial Density
Perturbations
Submitted: 2004-07-21, last modified: 2004-09-29
The anthropic principle has been proposed as an explanation for the observed
value of the cosmological constant. Here we revisit this proposal by allowing
for variation between universes in the amplitude of the scale-invariant
primordial cosmological density perturbations. We derive a priori probability
distributions for this amplitude from toy inflationary models in which the
parameter of the inflaton potential is smoothly distributed over possible
universes. We find that for such probability distributions, the likelihood that
we live in a typical, anthropically-allowed universe is generally quite small.
[131]
oai:arXiv.org:astro-ph/0407058 [pdf] - 65900
Cosmological simulations of the intracluster medium
Submitted: 2004-07-02, last modified: 2004-09-10
We investigate the properties of the intracluster medium (ICM) that forms
within N-body/hydrodynamical simulations of galaxy clusters in a \LambdaCDM
cosmology. When radiative cooling and a simple model for galactic feedback are
included, our clusters have X-ray luminosities and temperatures in good
agreement with observed systems, demonstrating the required excess entropy in
their cores. More generally, cooling and feedback increases the entropy of the
ICM everywhere, albeit without significantly affecting the slope of the profile
(S prop. to r) at large radii. The temperature of the ICM is only modestly
increased by these processes, with projected temperature profiles being in
reasonable agreement with the observations. Star/galaxy formation is still too
efficient in our simulations, however, and so our gas mass fractions are around
60 per cent of the observed value at r_{2500}. Finally, we examine the
reliability of using the hydrostatic equilibrium equation to estimate cluster
masses and find that it underpredicts the true mass of our clusters by up to 20
per cent, due to incomplete thermalisation of the gas. Feedback reduces this
discrepancy, however, with estimates being accurate to within 10 per cent out
to r_{500}.
[132]
oai:arXiv.org:astro-ph/0404589 [pdf] - 64497
The subhalo populations of LCDM dark haloes
Submitted: 2004-04-29, last modified: 2004-09-02
We use a large cosmological simulation and a variety of high resolution
resimulations of individual cluster and galaxy haloes to study the systematics
of subhalo populations over ranges of 1000 in halo mass and 1000 in the ratio
of subhalo to parent halo mass. The subhalo populations of different haloes are
not scaled copies of each other, but vary systematically with halo properties.
On average, the amount of substructure increases with halo mass. At fixed mass,
it decreases with halo concentration and with halo formation redshift. These
trends are comparable in size to the scatter in subhalo abundance between
similar haloes. Averaged over all haloes of given mass, the abundance of low
mass subhaloes per unit parent halo mass is independendent of parent mass.It is
very similar to the abundance per unit mass of low mass haloes in the universe
as a whole, once differing boundary definitions for subhaloes and haloes are
accounted for. The radial distribution of subhaloes within their parent haloes
is substantially less centrally concentrated than that of the dark matter. It
varies at most weakly with the mass(or concentration) of the parent halo and
not at all with subhalo mass. It does depend on the criteria used to define the
subhalo population considered. About 90% of present-day subhaloes were accreted
after z=1 and about 70% after z=0.5. Only about 8% the total mass of all haloes
accreted at z=1 survives as bound subhaloes at z=0. For haloes accreted at z=2,
the survival mass fraction is just 2%. Subhaloes seen near the centre of their
parent typically were accreted earlier and retain less of their original mass
than those seen near the edge. These strong systematics mean that comparison
with galaxies in real clusters is only possible if the formation of the
luminous component is modelled appropriately.
[133]
oai:arXiv.org:astro-ph/0408132 [pdf] - 66623
Disk Galaxy Rotation Curves in Triaxial CDM Halos
Hayashi, E.;
Navarro, J. F.;
Jenkins, A.;
Frenk, C. S.;
Power, Chris;
White, S. D. M.;
Springel, V.;
Stadel, J.;
Quinn, T.;
Wadsley, J.
Submitted: 2004-08-06
We use N-body hydrodynamical simulations to study the structure of disks in
triaxial potentials resembling CDM halos. Our analysis focuses on the accuracy
of the dark mass distribution inferred from rotation curves derived from
simulated long-slit spectra. We consider a massless disk embedded in a halo
with axis ratios of 0.5:0.6:1.0 and with its rotation axis aligned with the
minor axis of the halo. Closed orbits for the gaseous particles deviate from
coplanar circular symmetry, resulting in a variety of long-slit rotation curve
shapes, depending on the orientation of the disk relative to the line of sight.
Rotation curves may thus differ significantly from the spherically-averaged
circular velocity profile of the dark matter halo. "Solid-body" rotation
curves--typically interpreted as a signature of a constant density core in the
dark matter distribution--are obtained about 25% of the time for random
orientations although the dark matter follows the cuspy density profile
proposed by Navarro, Frenk & White (NFW). We conclude that the discrepancies
reported between the shape of the rotation curve of low surface brightness
galaxies and the structure of CDM halos may be resolved once the complex
effects of halo triaxiality on the dynamics of the gas component is properly
taken into account.
[134]
oai:arXiv.org:astro-ph/0406043 [pdf] - 65200
Gradient instability for w < -1
Submitted: 2004-06-01, last modified: 2004-07-07
We show that in single scalar field models of the dark energy with equations
of state satisfying $w \equiv p / \rho < -1$, the effective Lagrangian for
fluctuations about the homogeneous background has a wrong sign spatial kinetic
term. In most cases, spatial gradients are ruled out by microwave background
observations. The sign of $w+1$ is not connected to the sign of the time
derivative kinetic term in the effective Lagrangian.
[135]
oai:arXiv.org:astro-ph/0312499 [pdf] - 61712
Early Formation and Late Merging of the Giant Galaxies
Submitted: 2003-12-18, last modified: 2004-06-19
The most luminous galaxies in the present Universe are found at the centers
of the most massive dark matter haloes, rich galaxy clusters. In the LCDM
cosmology, such massive halo cores are present at redshift z=6 with a comoving
number density (as a function of mass interior to ~10 kpc) that is comparable
to today's value. The identity of the matter in these central regions is,
however, predicted to change as major mergers bring together stars and dark
matter from initially well separated sub-units. We use N-body simulations to
investigate how these mergers push pre-existing matter outwards in the dominant
galaxy while preserving the inner density profile of collisionless matter. It
appears that the central regions of large galaxies end up dominated by stars
formed in a number of dense cores, well before the last major mergers. The
density profile of collisionless matter (stars and dark matter combined) in
these central regions appears to be stable and to have attractor-like behavior
under merging. This suggests that the baryon loading associated with
dissipative contraction and star formation may be erased as subsequent mergers
drive the mass distribution back to a universal profile. Such suppression of
the effects of baryon loading, along with the early assembly of mass
concentrations, may help resolve some apparent challenges to the CDM model for
structure formation.
[136]
oai:arXiv.org:astro-ph/0405010 [pdf] - 64530
Galaxies and subhaloes in LCDM galaxy clusters
Submitted: 2004-05-03, last modified: 2004-06-04
We combine 10 high resolution resimulations of cluster--sized dark haloes
with semi--analytic galaxy formation modelling in order to compare the number
density and velocity dispersion profiles of cluster galaxies with those of dark
matter substructures (subhaloes). While the radial distribution of galaxies
follows closely that of the dark matter, the distribution of dark matter
subhaloes is much less centrally concentrated. The velocity dispersion profiles
of galaxies are also very similar to those of the dark matter, while those for
subhaloes are biased high, particularly in the inner regions of clusters. We
explain how these differences, already clearly visible in earlier work, are a
consequence of the formation of galaxies at the centres of dark matter haloes.
Galaxies and subhaloes represent different populations and are not directly
comparable. Evolution produces a complex and strongly position--dependent
relation between galaxies and the subhaloes in which they reside. This relation
can be properly modelled only by appropriate physical representation of the
galaxy formation process.
[137]
oai:arXiv.org:astro-ph/0311231 [pdf] - 60811
The Inner Structure of LambdaCDM Halos III: Universality and Asymptotic
Slopes
Submitted: 2003-11-11, last modified: 2004-01-06
We investigate the mass profile of LambdaCDM halos using a suite of numerical
simulations spanning five decades in halo mass, from dwarf galaxies to rich
galaxy clusters. Our analysis confirms the proposal of Navarro, Frenk & White
(NFW) that the shape of LambdaCDM halo mass profiles differs strongly from a
power law and depends little on mass. The logarithmic slope of the
spherically-averaged density profile, as measured by beta=-dln(rho)/dln(r),
decreases monotonically towards the center and becomes shallower than
isothermal (beta<2) inside a characteristic radius, r_{-2}. Although the
fitting formula proposed by NFW provides a reasonably good approximation to the
density and circular velocity profiles of individual halos, systematic
deviations from the best NFW fits are also noticeable. Inside r_{-2}, the
profile of simulated halos gets shallower with radius more gradually than
predicted and, as a result, NFW fits tend to underestimate the dark matter
density in these regions. This discrepancy has been interpreted as indicating a
steeply divergent cusp, but our results suggest a different interpretation. We
use the density and enclosed mass at our innermost resolved radii to place
strong constraints on beta_{0}: density cusps as steep as r^{-1.5} are
inconsistent with most of our simulations, although beta_{0}=1 is still
consistent with our data. Our density profiles show no sign of converging to a
well-defined asymptotic inner power law. We propose a simple formula that
reproduces the radial dependence of the slope better than the NFW profile, and
so may minimize errors when extrapolating our results inward to radii not yet
reliably probed by numerical simulations.
[138]
oai:arXiv.org:astro-ph/0310576 [pdf] - 60239
The Inner Structure of LCDM Halos II: Halo Mass Profiles and LSB
Rotation Curves
Submitted: 2003-10-20
We use a set of high-resolution cosmological N-body simulations to
investigate the inner mass profile of galaxy-sized cold dark matter (CDM)
halos. These simulations extend the thorough numerical convergence study
presented in Paper I of this series (Power et al. 2003), and demonstrate that
the mass profile of CDM halos can be robustly estimated beyond a minimum
converged radius of order r_conv ~ 1 kpc/h in our highest resolution runs. The
density profiles of simulated halos become progressively shallow from the
virial radius inwards, and show no sign of approaching a well-defined power-law
behaviour near the centre. At r_conv, the logarithmic slope of the density
profile is steeper than the asymptotic \rho \propto r^-1 expected from the
formula proposed by Navarro, Frenk, and White (1996), but significantly
shallower than the steeply divergent \rho \propto r^-1.5 cusp proposed by Moore
et al. (1999). We perform a direct comparison of the spherically-averaged dark
matter circular velocity (V_c) profiles with rotation curves of low surface
brightness (LSB) galaxies from the samples of de Blok et al. (2001), de Blok
and Bosma (2002), and Swaters et al. (2003). Most (about two-thirds) LSB
galaxies in this dataset are roughly consistent with CDM halo V_c profiles.
However, about one third of LSBs in these samples feature a sharp transition
between the rising and flat part of the rotation curve that is not seen in the
V_c profiles of CDM halos. This discrepancy has been interpreted as excluding
the presence of cusps, but we argue that it might simply reflect the difference
between circular velocity and gas rotation speed likely to arise in gaseous
disks embedded within realistic, triaxial CDM halos.
[139]
oai:arXiv.org:astro-ph/0305286 [pdf] - 56759
Cosmic Structure Growth and Dark Energy
Submitted: 2003-05-15, last modified: 2003-08-20
Dark energy has a dramatic effect on the dynamics of the universe, causing
the recently discovered acceleration of the expansion. The dynamics are also
central to the behavior of the growth of large scale structure, offering the
possibility that observations of structure formation provide a sensitive probe
of the cosmology and dark energy characteristics. In particular, dark energy
with a time varying equation of state can have an influence on structure
formation stretching back well into the matter dominated epoch. We analyze this
impact, first calculating the linear perturbation results, including those for
weak gravitational lensing. These dynamical models possess definite observable
differences from constant equation of state models. Then we present a large
scale numerical simulation of structure formation, including the largest volume
to date involving a time varying equation of state. We find the halo mass
function is well described by the Jenkins et al. mass function formula. We also
show how to interpret modifications of the Friedmann equation in terms of a
time variable equation of state. The results presented here provide steps
toward realistic computation of the effect of dark energy in cosmological
probes involving large scale structure, such as cluster counts,
Sunyaev-Zel'dovich effect, or weak gravitational lensing.
[140]
oai:arXiv.org:astro-ph/0306568 [pdf] - 57649
Momentum transfer across shear flows in Smoothed Particle Hydrodynamic
simulations of galaxy formation
Submitted: 2003-06-26
We investigate the evolution of angular momentum in SPH simulations of galaxy
formation, paying particular attention to artificial numerical effects. We find
that a cold gas disc forming in an ambient hot gas halo receives a strong
hydrodynamic torque from the hot gas. By splitting the hydrodynamic force into
artificial viscosity and pressure gradients, we find that the angular momentum
transport is caused not by the artificial viscosity but by the pressure
gradients. Using simple test simulations of shear flows, we conclude that the
pressure gradient-based viscosity can be divided into two components: one due
to the noisiness of SPH and the other to ram pressure. The former is
problematic even with very high resolution because increasing resolution does
not reduce the noisiness. On the other hand, the ram pressure effect appears
only when a cold gas disc or sheet does not contain enough particles. In such a
case, holes form in the disc or sheet, and then ram pressure from intra-hole
hot gas, causes significant deceleration. In simulations of galactic disc
formation, star formation usually decreases the number of cold gas particles,
and hole formation leads to the fragmentation of the disc. To circumvent these
problem, we modify the SPH algorithm, decoupling the cold from the hot gas
phases, i.e. inhibiting the hydrodynamic interaction between cold and hot
particles. This, a crude modelling of a multi-phase fluid in SPH cosmological
simulations, leads to the formation of smooth extended cold gas discs and to
better numerical convergence. The decoupling is applicable in so far as the
self-gravitating gas disc with negligible external pressure is a good
approximation for a cold gas disc. (abridged)
[141]
oai:arXiv.org:astro-ph/0207664 [pdf] - 354083
Stable clustering, the halo model and nonlinear cosmological power
spectra
Smith, R. E.;
Peacock, J. A.;
Jenkins, A.;
White, S. D. M.;
Frenk, C. S.;
Pearce, F. R.;
Thomas, P. A.;
Efstathiou, G.;
Couchmann, H. M. P.;
Consortium, The Virgo
Submitted: 2002-07-31, last modified: 2003-06-11
We present the results of a large library of cosmological N-body simulations,
using power-law initial spectra. The nonlinear evolution of the matter power
spectra is compared with the predictions of existing analytic scaling formulae
based on the work of Hamilton et al. The scaling approach has assumed that
highly nonlinear structures obey `stable clustering' and are frozen in proper
coordinates. Our results show that, when transformed under the self-similarity
scaling, the scale-free spectra define a nonlinear locus that is clearly
shallower than would be required under stable clustering. Furthermore, the
small-scale nonlinear power increases as both the power-spectrum index n and
the density parameter Omega decrease, and this evolution is not well accounted
for by the previous scaling formulae. This breakdown of stable clustering can
be understood as resulting from the modification of dark-matter haloes by
continuing mergers. These effects are naturally included in the analytic `halo
model' for nonlinear structure; using this approach we are able to fit both our
scale-free results and also our previous CDM data. This approach is more
accurate than the commonly-used Peacock--Dodds formula and should be applicable
to more general power spectra. Code to evaluate nonlinear power spectra using
this method is available from
http://as1.chem.nottingham.ac.uk/~res/software.html Following publication, we
will make the power-law simulation data available through the Virgo website
http://www.mpa-garching.mpg.de/Virgo
[142]
oai:arXiv.org:astro-ph/0212357 [pdf] - 358218
The Halo Occupation Distribution and the Physics of Galaxy Formation
Berlind, Andreas A.;
Weinberg, David H.;
Benson, Andrew J.;
Baugh, Carlton M.;
Cole, Shaun;
Davé, Romeel;
Frenk, Carlos S.;
Jenkins, Adrian;
Katz, Neal;
Lacey, Cedric G.
Submitted: 2002-12-16, last modified: 2003-04-30
The halo occupation distribution (HOD) describes the bias between galaxies
and dark matter by specifying (a) the probability P(N|M) that a halo of virial
mass M contains N galaxies of a particular class and (b) the relative
distributions of galaxies and dark matter within halos. We calculate predicted
HODs for a Lambda-CDM cosmological model using an SPH hydrodynamic simulation
and a semi-analytic (SA) galaxy formation model. Although the two methods
predict different galaxy mass functions, their HOD predictions agree remarkably
well. For mass-selected samples, the mean occupation <N(M)> exhibits a sharp
cutoff at low halo masses, a slowly rising plateau for <N>~1-2, and a more
steeply rising high occupancy regime. At low <N>, the mean pair and triple
counts are well below Poisson expectations, with important consequences for
small scale behavior of 2- and 3-point correlation functions. The HOD depends
strongly on galaxy age, with high mass halos populated mainly by old galaxies
and low mass halos by young galaxies. The SPH simulation supports several
simplifying assumptions about HOD bias: the most massive galaxy in a halo
usually lies close to the center and moves near the halo's mean velocity;
satellite galaxies have the same radial profile and velocity dispersion as the
dark matter; and the mean occupation at fixed halo mass is independent of the
halo's larger scale environment. By applying the SPH and SA HODs to a large
volume N-body simulation, we show that both methods predict slight, observable
departures from a power-law galaxy correlation function. The predicted HODs are
closely tied to the underlying galaxy formation physics, they offer useful
guidance to theoretical models of galaxy clustering, and they will be tested
empirically by ongoing analyses of galaxy redshift surveys. (Shortened)
[143]
oai:arXiv.org:astro-ph/0201544 [pdf] - 47488
The Inner Structure of LambdaCDM Halos I: A Numerical Convergence Study
Submitted: 2002-01-31, last modified: 2003-01-14
We present a comprehensive set of convergence tests which explore the role of
various numerical parameters on the equilibrium structure of a simulated dark
matter halo. We report results obtained with two independent, state-of-the-art,
multi-stepping, parallel N-body codes: PKDGRAV and GADGET. We find that
convergent mass profiles can be obtained for suitable choices of the
gravitational softening, timestep, force accuracy, initial redshift, and
particle number. For softenings chosen so that particle discreteness effects
are negligible, convergence in the circular velocity is obtained at radii where
the following conditions are satisfied: (i) the timestep is much shorter than
the local orbital timescale; (ii) accelerations do not exceed a characteristic
acceleration imprinted by the gravitational softening; and (iii) enough
particles are enclosed so that the collisional relaxation timescale is longer
than the age of the universe. The most stringent requirement for convergence is
typically that imposed on the particle number by the collisional relaxation
criterion, which implies that in order to estimate accurate circular velocities
at radii where the density contrast may reach $\sim 10^6$, the region must
enclose of order 3000 particles (or more than a few times $10^6$ within the
virial radius). Applying these criteria to a galaxy-sized $\Lambda$CDM halo, we
find that the spherically-averaged density profile becomes progressively
shallower from the virial radius inwards, reaching a logarithmic slope
shallower than -1.2 at the innermost resolved point, $r \sim 0.005 r_{200}$,
with little evidence for convergence to a power-law behaviour in the inner
regions.
[144]
oai:arXiv.org:astro-ph/0110246 [pdf] - 45303
Galaxy Clusters in Hubble Volume Simulations: Cosmological Constraints
from Sky Survey Populations
Evrard, A. E.;
MacFarland, T. J.;
Couchman, H. M. P.;
Colberg, J. M.;
Yoshida, N.;
White, S. D. M.;
Jenkins, A. R.;
Frenk, C. S.;
Pearce, F. R.;
Peacock, J. A.;
Thomas, P. A.
Submitted: 2001-10-10, last modified: 2002-03-11
We use giga-particle N-body simulations to study galaxy cluster populations
in Hubble Volumes of LCDM (Omega_m=0.3, Omega_Lambda=0.7) and tCDM (Omega_m=1)
world models. Mapping past light-cones of locations in the computational space,
we create mock sky surveys of dark matter structure to z~1.4 over 10,000 sq deg
and to z~0.5 over two full spheres. Calibrating the Jenkins mass function at
z=0 with samples of ~1.5 million clusters, we show that the fit describes the
sky survey counts to <~20% acccuracy over all redshifts for systems larger than
poor groups (M>5e13 Msun/h). Fitting the observed local temperature function
determines the ratio beta of specific thermal energies in dark matter and
intracluster gas. We derive a scaling with power spectrum normalization beta
\propto sigma8^{5/3}, and measure a 4% error on sigma8 arising from cosmic
variance in temperature-limited cluster samples. Considering distant clusters,
the LCDM model matches EMSS and RDCS X-ray-selected survey observations under
economical assumptions for intracluster gas evolution. Using transformations of
mass-limited cluster samples that mimic sigma8 variation, we explore SZ search
expectations for a 10 sq deg survey complete above 10^{14} Msun/h. Cluster
counts are shown to be extremely sensitive to sigma8 uncertainty while redshift
statistics, such as the sample median, are much more stable. For LCDM, the
characteristic temperature at fixed sky surface density is a weak function of
redshift, implying an abundance of hot clusters at z>1. Assuming constant beta,
four kT>8 keV clusters lie at z>2 and 40 kT>5 keV clusters lie at z>3 on the
whole sky. Detection of Coma-sized clusters at z>1 violate LCDM at 95%
confidence if their surface density exceeds 0.003 per sq deg, or 120 on the
whole sky.
[145]
oai:arXiv.org:astro-ph/0106462 [pdf] - 43258
Including star formation and supernova feedback within cosmological
simulations of galaxy formation
Submitted: 2001-06-26, last modified: 2001-10-23
We investigate phenomenological models of star formation and supernova
feedback in N-body/SPH simulations of galaxy formation. First, we compare
different prescriptions in the literature for turning cold gas into stars
neglecting feedback effects. We find that most prescriptions give broadly
similar results: the ratio of cold gas to stars in the final galaxies is
primarily controlled by the range of gas densities where star formation is
allowed to proceed efficiently. In the absence of feedback, the fraction of gas
that cools is much too high resulting, for example, in a K-band luminosity
function that is much brighter than observed. This problem is ameliorated by
including a feedback model which either imparts radial kinetic perturbations to
galactic gas or directly reheats such material and prevents it from cooling for
a certain period of time. In both these models, a significant fraction of cold
gas is heated and expelled from haloes with an efficiency that varies inversely
with with halo circular velocity. Increasing the resolution of a simulation
allows a wider dynamic range in mass to be followed, but the average properties
of the resolved galaxy population remain largely unaffected. However, as the
resolution is increased, more and more gas is reheated by small galaxies; our
results suggest that convergence requires the full mass range of galaxies to be
resolved.
[146]
oai:arXiv.org:astro-ph/0011212 [pdf] - 39219
Simulations of Deep Pencil-Beam Redshift Surveys
Yoshida, N.;
Colberg, J.;
White, S. D. M.;
Evrard, A. E.;
MacFarland, T. J.;
Couchman, H. M. P.;
Jenkins, A.;
Frenk, C. S.;
Pearce, F. R.;
Efstathiou, G.;
Peacock, J. A.;
Thomas, P. A.
Submitted: 2000-11-10, last modified: 2001-03-13
We create mock pencil-beam redshift surveys from very large cosmological
$N$-body simulations of two Cold Dark Matter cosmogonies, an Einstein-de Sitter
model ($\tau$CDM) and a flat model with $\Omega_0 =0.3$ and a cosmological
constant ($\Lambda$CDM). We use these to assess the significance of the
apparent periodicity discovered by Broadhurst et al. (1990). Simulation
particles are tagged as `galaxies' so as to reproduce observed present-day
correlations. They are then identified along the past light-cones of
hypothetical observers to create mock catalogues with the geometry and the
distance distribution of the Broadhurst et al. data. We produce 1936 (2625)
quasi-independent catalogues from our $\tau$CDM ($\Lambda$CDM) simulation. A
couple of large clumps in a catalogue can produce a high peak at low
wavenumbers in the corresponding one-dimensional power spectrum, without any
apparent large-scale periodicity in the original redshift histogram. Although
the simulated redshift histograms frequently display regularly spaced clumps,
the spacing of these clumps varies between catalogues and there is no
`preferred' period over our many realisations. We find only a 0.72 (0.49) per
cent chance that the highest peak in the power spectrum of a $\tau$CDM
($\Lambda$CDM) catalogue has a peak-to-noise ratio higher than that in the
Broadhurst et al. data. None of the simulated catalogues with such high peaks
shows coherently spaced clumps with a significance as high as that of the real
data. We conclude that in CDM universes, the kind of regularity observed by
Broadhurst et al. has a priori probability well below $10^{-3}$.
[147]
oai:arXiv.org:astro-ph/0010587 [pdf] - 38934
Simulations of galaxy formation in a cosmological volume
Submitted: 2000-10-30
We present results of large N-body-hydrodynamic simulations of galaxy
formation. Our simulations follow the formation of galaxies in cubic volumes of
side 100Mpc, in two versions of the cold dark matter (CDM) cosmogony: the
standard, Omega=1 SCDM model and the flat, Omega=0.3 LCDM model. Over 2000
galaxies form in each of these simulations. We examine the rate at which gas
cools and condenses into dark matter halos. This roughly tracks the cosmic star
formation rate inferred from observations at various redshifts. Galaxies in the
simulations form gradually over time in the hierarchical fashion characteristic
of the CDM cosmogony. In the LCDM model, substantial galaxies first appear at
z=5 and the population builds up rapidly until z=1 after which the rate of
galaxy formation declines as cold gas is consumed and the cooling time of hot
gas increases. In both cosmologies, the present-day K-band luminosity function
of the simulated galaxies resembles observations. The galaxy autocorrelation
functions differ significantly from those of the dark matter. At the present
epoch there is little bias in either model between galaxies and dark matter on
large scales, but a significant anti-bias on scales of 1Mpc/h and a positive
bias on scales of 100kpc/h. The galaxy correlation function evolves little with
redshift in the range z=0-3, and depends on the luminosity of the galaxy
sample. The projected pairwise velocity dispersion of the galaxies is much
lower than that of the dark matter on scales less than 2Mpc/h. Applying a
virial mass estimator to the largest galaxy clusters recovers the cluster
virial masses in an unbiased way. (Abridged)
[148]
oai:arXiv.org:astro-ph/0005260 [pdf] - 36032
Mass function of dark matter halos
Submitted: 2000-05-12, last modified: 2000-09-25
We combine data from a number of N-body simulations to predict the abundance
of dark halos in Cold Dark Matter universes over more than 4 orders of
magnitude in mass. A comparison of different simulations suggests that the
dominant uncertainty in our results is systematic and is smaller than 10--30%
at all masses, depending on the halo definition used. In particular, our
``Hubble Volume'' simulations of \tcdm and \lcdm cosmologies allow the
abundance of massive clusters to be predicted with uncertainties well below
those expected in all currently planned observational surveys. We show that for
a range of CDM cosmologies and for a suitable halo definition, the simulated
mass function is almost independent of epoch, of cosmological parameters, and
of initial power spectrum when expressed in appropriate variables. This
universality is of exactly the kind predicted by the familiar Press-Schechter
model, although this model predicts a mass function shape which differs from
our numerical results, overestimating the abundance of ``typical'' halos and
underestimating that of massive systems.
[149]
oai:arXiv.org:astro-ph/9912220 [pdf] - 109883
A Comparison of Semi-Analytic and Smoothed Particle Hydrodynamics Galaxy
Formation
Submitted: 1999-12-10, last modified: 2000-08-25
We compare the statistical properties of galaxies found in two different
models of hierarchical galaxy formation: the semi-analytic model of Cole et al.
and the smoothed particle hydrodynamics (SPH) simulations of Pearce et al.
Using a `stripped-down' version of the semi-analytic model which mimics the
resolution of the SPH simulations and excludes physical processes not included
in them, we find that the two models produce an ensemble of galaxies with
remarkably similar properties, although there are some differences in the gas
cooling rates and in the number of galaxies that populate halos of different
mass. The full semi-analytic model, which has effectively no resolution limit
and includes a treatment of star formation and supernovae feedback, produces
somewhat different (but readily understandable) results. Agreement is
particularly good for the present-day global fractions of hot gas, cold dense
(i.e. galactic) gas and uncollapsed gas, for which the SPH and stripped-down
semi-analytic calculations differ by at most 25%. In the most massive halos,
the stripped-down semi-analytic model predicts, on the whole, up to 50% less
gas in galaxies than is seen in the SPH simulations. The two techniques
apportion this cold gas somewhat differently amongst galaxies in a given halo.
This difference can be tracked down to the greater cooling rate in massive
halos in the SPH simulation compared to the semi-analytic model. (abridged)
[150]
oai:arXiv.org:astro-ph/0007362 [pdf] - 37200
Public Release of N-body simulation and related data by the Virgo
consortium
Frenk, C. S.;
Colberg, J. M.;
Couchman, H. M. P.;
Efstathiou, G.;
Evrard, A. E.;
Jenkins, A.;
MacFarland, T. J.;
Moore, B.;
Peacock, J. A.;
Pearce, F. R.;
Thomas, P. A.;
White, S. D. M.;
Yoshida, N.
Submitted: 2000-07-25
We are making available on the WWW a selection of the archived data from
N-body simulations carried out by the Virgo consortium and related groups. This
currently includes: (i) time-slice, lightcone and cluster data from the two
$10^9$-particle Hubble volume simulations described by Evrard 1998; (ii)
time-slice data from simulations of 4 different cold dark matter cosmological
models with $256^3$ particles analysed by Jenkins et al 1998; (iii) Dark halo
catalogs, merger trees and galaxy catalogs from the GIF project described by
Kauffmann et al 1999. Basic software is supplied to read the data. The data can
be accessed from: http://www.mpa-garching.mpg.de/Virgo/data_download.html
[151]
oai:arXiv.org:astro-ph/0007348 [pdf] - 37186
A simluated TCDM cosmology cluster catalogue: the NFW profile and the
temperature-mass scaling relations
Submitted: 2000-07-24
We have extracted over 400 clusters, covering more than 2 decades in mass,
from three simulations of the TCDM cosmology. This represents the largest,
uniform catalogue of simulated clusters ever produced. The clusters exhibit a
wide variety of density-profiles. Only a minority are well-fit in their outer
regions by the widely used density profile of Navarro, Frenk & White (1977).
Others have steeper outer density profiles, show sharp breaks in their density
profiles, or have significant substructure. If we force a fit to the NFW
profile, then the best-fit concentrations decline with increasing mass, but
this is driven primarily by an increase in substructure as one moves to higher
masses. The measured temperature-mass relations for properties measured within
a sphere enclosing a fixed overdensity all follow the self-similar form,
T\propto M^{2/3}, however the normalisation is lower than in observed clusters.
The temperature-mass relations for properties measured within a fixed physical
radius are significantly steeper then this. Both can be accurately predicted
using the NFW model.
[152]
oai:arXiv.org:astro-ph/0005259 [pdf] - 36031
Clustering of Galaxy Clusters in CDM Universes
Colberg, J. M.;
White, S. D. M.;
Yoshida, N.;
MacFarland, T. J.;
Jenkins, A.;
Frenk, C. S.;
Pearce, F. R.;
Evrard, A. E.;
Couchman, H. M. P.;
Efstathiou, G.;
Peacock, J. A.;
Thomas, P. A.
Submitted: 2000-05-11, last modified: 2000-07-11
We use very large cosmological N--body simulations to obtain accurate
predictions for the two-point correlations and power spectra of mass-limited
samples of galaxy clusters. We consider two currently popular cold dark matter
(CDM) cosmogonies, a critical density model ($\tau$CDM) and a flat low density
model with a cosmological constant ($\Lambda$CDM). Our simulations each use
$10^9$ particles to follow the mass distribution within cubes of side
$2h^{-1}$Gpc ($\tau$CDM) and $3h^{-1}$Gpc ($\Lambda$CDM) with a force
resolution better than $10^{-4}$ of the cube side. We investigate how the
predicted cluster correlations increase for samples of increasing mass and
decreasing abundance. Very similar behaviour is found in the two cases. The
correlation length increases from $r_0=12$ -- 13$h^{-1}$Mpc for samples with
mean separation $d_{\rm c}=30h^{-1}$Mpc to $r_0=22$-- 27$h^{-1}$Mpc for samples
with $d_{\rm c}=100h^{-1}$Mpc. The lower value here corresponds to $\tau$CDM
and the upper to $\Lambda$CDM. The power spectra of these cluster samples are
accurately parallel to those of the mass over more than a decade in scale. Both
correlation lengths and power spectrum biases can be predicted to better than
10% using the simple model of Sheth, Mo & Tormen (2000). This prediction
requires only the linear mass power spectrum and has no adjustable parameters.
We compare our predictions with published results for the APM cluster sample.
The observed variation of correlation length with richness agrees well with the
models, particularly for $\Lambda$CDM. The observed power spectrum (for a
cluster sample of mean separation $d_{\rm c}=31h^{-1}$Mpc) lies significantly
above the predictions of both models.
[153]
oai:arXiv.org:astro-ph/9908107 [pdf] - 107779
Parameter Tests Within Cosmological Simulations of Galaxy Formation
Submitted: 1999-08-10, last modified: 2000-02-29
Numerical simulations of galaxy formation require a number of parameters.
Some of these are intrinsic to the numerical integration scheme (eg the
timestep), while others describe the physical model (eg the gas metallicity).
In this paper, we present results of a systematic exploration of the effects of
varying a subset of these parameters on simulations of galaxy formation. We use
N-body and ``Smoothed Particle Hydrodynamics'' techniques to follow the
evolution of cold dark matter and gas in a small volume. We compare a fiducial
model to 24 different simulations, in which one parameter at a time is varied,
focussing on properties such as the relative fraction of hot and cold gas, and
the abundance and masses of galaxies. We find that for reasonable choices of
numerical values, many parameters have relatively little effect on the
galaxies, with the notable exception of the parameters that control the
resolution of the simulation and the efficiency with which gas cools.
[154]
oai:arXiv.org:astro-ph/0002308 [pdf] - 275813
Collisional versus Collisionless Dark Matter
Submitted: 2000-02-16, last modified: 2000-02-17
We compare the structure and substructure of dark matter halos in model
universes dominated by collisional, strongly self interacting dark matter
(SIDM) and collisionless, weakly interacting dark matter (CDM). While SIDM
virialised halos are more nearly spherical than CDM halos, they can be
rotationally flattened by as much as 20% in their inner regions. Substructure
halos suffer ram-pressure truncation and drag which are more rapid and severe
than their gravitational counterparts tidal stripping and dynamical friction.
Lensing constraints on the size of galactic halos in clusters are a factor of
two smaller than predicted by gravitational stripping, and the recent detection
of tidal streams of stars escaping from the satellite galaxy Carina suggests
that its tidal radius is close to its optical radius of a few hundred parsecs
--- an order of magnitude smaller than predicted by CDM models but consistent
with SIDM. The orbits of SIDM satellites suffer significant velocity bias
$\sigma_{SIDM}/\sigma_{CDM}=0.85$ and are more circular than CDM,
$\beta_{SIDM}} \approx 0.5$, in agreement with the inferred orbits of the
Galaxy's satellites. In the limit of a short mean free path, SIDM halos have
singular isothermal density profiles, thus in its simplest incarnation SIDM is
inconsistent with galactic rotation curves.
[155]
oai:arXiv.org:astro-ph/9912238 [pdf] - 109901
Experimental Cosmic Statistics II: Distribution
Submitted: 1999-12-10
Colombi et al. 1999 (paper I) investigated the counts-in-cells statistics and
their respective errors in the $\tau$CDM Virgo Hubble Volume simulation. This
extremely large $N$-body experiment also allows a numerical investigation of
the {\em cosmic distribution function}, $\Upsilon(\tA)$ itself for the first
time. For a statistic $A$, $\Upsilon(\tA)$ is the probability density of
measuring the value $\tA$ in a finite galaxy catalog. $\Upsilon$ was evaluated
for the distribution of counts-in-cells, $P_N$, the factorial moments, $F_k$,
and the cumulants, $\xiav$ and $S_N$'s, using the same subsamples as paper I.
While paper I concentrated on the first two moments of $\Upsilon$, i.e. the
mean, the cosmic error and the cross-correlations, here the function $\Upsilon$
is studied in its full generality, including a preliminary analysis of joint
distributions $\Upsilon(\tA,\tB)$. The most significant, and reassuring result
for the analyses of future galaxy data is that the cosmic distribution function
is nearly Gaussian provided its variance is small. A good practical criterion
for the relative cosmic error is that $\Delta A/A \la 0.2$. This means that for
accurate measurements, the theory of the cosmic errors, presented by Szapudi &
Colombi (1996) and Szapudi, Colombi & Bernardeau (1999), and confirmed
empirically by paper I, is sufficient for a full statistical description and
thus for a maximum likelihood rating of models. As the cosmic error increases,
the cosmic distribution function $\Upsilon$ becomes increasingly skewed and is
well described by a generalization of the lognormal distribution. The cosmic
skewness is introduced as an additional free parameter. (...more in paper...)
[156]
oai:arXiv.org:astro-ph/9912236 [pdf] - 109899
Experimental Cosmic Statistics I: Variance
Submitted: 1999-12-10
Counts-in-cells are measured in the $\tau$CDM Virgo Hubble Volume simulation.
This large N-body experiment has 10^9 particles in a cubic box of size 2000
h^{-1} Mpc. The unprecedented combination of size and resolution allows for the
first time a realistic numerical analysis of the cosmic errors and cosmic
correlations of statistics related to counts-in-cells measurements, such as the
probability distribution function P_N itself, its factorial moments F_k and the
related cumulants $\xiav$ and S_N's. These statistics are extracted from the
whole simulation cube, as well as from 4096 sub-cubes of size 125 h^{-1}Mpc,
each representing a virtual random realization of the local universe. The
measurements and their scatter over the sub-volumes are compared to the
theoretical predictions of Colombi, Bouchet & Schaeffer (1995) for P_0, and of
Szapudi & Colombi (1996, SC) and Szapudi, Colombi & Bernardeau (1999a, SCB) for
the factorial moments and the cumulants. The general behavior of experimental
variance and cross-correlations as functions of scale and order is well
described by theoretical predictions, with a few percent accuracy in the weakly
non-linear regime for the cosmic error on factorial moments. (... more in paper
>...)
[157]
oai:arXiv.org:astro-ph/9907189 [pdf] - 107414
Inhomogeneous Reionization Regulated by Radiative and Stellar Feedbacks
Submitted: 1999-07-14
We study the inhomogeneous reionization in a critical density CDM universe
due to stellar sources, including Population III objects. The spatial
distribution of the sources is obtained from high resolution numerical N-body
simulations. We calculate the source properties taking into account a
self-consistent treatment of both radiative (ie ionizing and H2
-photodissociating photons) and stellar (ie SN explosions) feedbacks regulated
by massive stars. This allows us to describe the topology of the ionized and
dissociated regions at various cosmic epochs and derive the evolution of H, He,
and H2 filling factors, soft UV background, cosmic star formation rate and the
final fate of ionizing objects. The main results are: (i) galaxies reionize the
IGM by z~10 (with some uncertainty related to the gas clumping factor), whereas
H2 is completely dissociated already by z~25; (ii) reionization is mostly due
to the relatively massive objects which collapse via H line cooling, while
objects whose formation relies on H2 cooling alone are insufficient to this
aim; (iii) the diffuse soft UV background is the major source of radiative
feedback effects for z<15; at higher z direct flux from neighboring objects
dominates; (iv) the match of the calculated cosmic star formation history with
the one observed at lower redshifts suggests that the conversion efficiency of
baryons into stars is ~1%; (v) we find that a very large population of dark
objects which failed to form stars is present by z~8. We discuss and compare
our results with similar previous studies.
[158]
oai:arXiv.org:astro-ph/9906160 [pdf] - 106881
The Santa Barbara cluster comparison project: a comparison of
cosmological hydrodynamics solutions
Frenk, C. S.;
White, S. D. M.;
Bode, P.;
Bond, J. R.;
Bryan, G. L.;
Cen, R.;
Couchman, H. M. P.;
Evrard, A. E.;
Gnedin, N.;
Jenkins, A.;
Khokhlov, A. M.;
Klypin, A.;
Navarro, J. F.;
Norman, M. L.;
Ostriker, J. P.;
Owen, J. M.;
Pearce, F. R.;
Pen, U. -L.;
Steinmetz, M.;
Thomas, P. A.;
Villumsen, J. V.;
Wadsley, J. W.;
Warren, M. S.;
Xu, G.;
Yepes, G.
Submitted: 1999-06-09
We have simulated the formation of an X-ray cluster in a cold dark matter
universe using 12 different codes. The codes span the range of numerical
techniques and implementations currently in use, including SPH and grid methods
with fixed, deformable or multilevel meshes. The goal of this comparison is to
assess the reliability of cosmological gas dynamical simulations of clusters in
the simplest astrophysically relevant case, that in which the gas is assumed to
be non-radiative. We compare images of the cluster at different epochs, global
properties such as mass, temperature and X-ray luminosity, and radial profiles
of various dynamical and thermodynamical quantities. On the whole, the
agreement among the various simulations is gratifying although a number of
discrepancies exist. Agreement is best for properties of the dark matter and
worst for the total X-ray luminosity. Even in this case, simulations that
adequately resolve the core radius of the gas distribution predict total X-ray
luminosities that agree to within a factor of two. Other quantities are
reproduced to much higher accuracy. For example, the temperature and gas mass
fraction within the virial radius agree to about 10%, and the ratio of specific
kinetic to thermal energies of the gas agree to about 5%. Various factors
contribute to the spread in calculated cluster properties, including
differences in the internal timing of the simulations. Based on the overall
consistency of results, we discuss a number of general properties of the
cluster we have modelled.
[159]
oai:arXiv.org:astro-ph/9906032 [pdf] - 106753
Cosmological galaxy formation
Pearce, F. R.;
Frenk, C. S.;
Jenkins, A.;
Colberg, J. M.;
Thomas, P. A.;
Couchman, H. M. P.;
White, S. D. M.;
Efstathiou, G.;
Peacock, J. A.;
Nelson, A. H.
Submitted: 1999-06-02
We discuss the first results from two successful simulations of galaxy
formation within a cosmological volume. With over 2000 large objects forming in
each we have sufficient numbers to reliably produce both galaxy correlation and
luminosity functions. We find that the observed galaxy counts are well fitted
by these models and that the galaxies display an almost un-evolving correlation
function back to a redshift of 3 which closely resembles the featureless
observed form and amplitude.
[160]
oai:arXiv.org:astro-ph/9906039 [pdf] - 106760
Galaxy clustering determined from numerical cosmological simulations
Jenkins, Adrian;
Frenk, C. S.;
Pearce, F. R.;
Thomas, P. A.;
Colberg, J. M.;
White, S. D. M.;
Couchman, H. M. P.;
Peacock, J. A.;
Efstathiou, G.;
Nelson, A. H.
Submitted: 1999-06-02
We have simulated the growth of structure in two 100 Mpc boxes for LCDM and
SCDM universes. These N-body/SPH simulations include a gaseous component which
is able to cool radiatively. A fraction of the gas cools into cold dense
objects which we identify as galaxies. In this article we give a preliminary
analysis of the clustering behaviour of these galaxies concentrating on the
LCDM model. We find a galaxy correlation function which is very close to a
power law and which evolves relatively little with redshift. The pairwise
dispersions of the galaxies are significantly lower than the dark matter. The
LCDM model gives a surprisingly good match to the observational determinations
of the galaxy correlation function and pairwise dispersions.
[161]
oai:arXiv.org:astro-ph/9905160 [pdf] - 106486
A simulation of galaxy formation and clustering
Submitted: 1999-05-13
We discuss early results from the first large N-body/hydrodynamical
simulation to resolve the formation of galaxies in a cold dark matter universe.
The simulation follows the formation of galaxies by gas cooling within dark
halos of mass a few times $10^{11}\Msun$ and above, in a flat universe with a
positive cosmological constant. Over 2200 galaxies form in our simulated volume
of $(100 \Mpc)^3$. Assigning luminosities to the model galaxies using a
spectral population synthesis model results in a K-band luminosity function in
excellent agreement with observations. The two-point correlation function of
galaxies in the simulation evolves very little since $z=3$ and has a shape
close to a power-law over four orders of magnitude in amplitude. At the present
day, the galaxy correlation function in the simulation is antibiased relative
to the mass on small scales and unbiased on large scales. It provides a
reasonable match to observations.
[162]
oai:arXiv.org:astro-ph/9808257 [pdf] - 102601
Galaxy Clusters in the Hubble Volume Simulations
Colberg, Joerg M.;
White, Simon D. M.;
MacFarland, Thomas J.;
Jenkins, Adrian;
Frenk, Carlos S.;
Pearce, F. R.;
Evrard, August E.;
Couchman, H. M. P.;
Efstathiou, George;
Peacock, John A.;
Thomas, Peter A.
Submitted: 1998-08-24
We report on analyses of cluster samples obtained from the Hubble Volume
Simulations. These simulations, an $\Omega=1$ model named $\tau$CDM and a flat
low $\Omega$ model with a cosmological constant ($\Lambda$CDM), comprise the
largest computational efforts to date in numerical cosmology. We investigate
the presence of massive galaxy clusters at $z\approx 0.8$. The $\tau$CDM model
fails to form clusters at such a redshift. However, due to the small number of
observed clusters around $z\approx 0.8$ and the uncertainties in the
determinations of their masses, this conclusion still is somewhat preliminary.
We produce cluster catalogs at $z=0$ for both cosmologies and investigate their
two--point correlation function $\xi$. We show that the relationship between
the mean density of subsamples of clusters, expressed via their mean separation
$d_{\rm c}$, and the correlation length $r_0$, defined through $\xi(r_0) = 1$,
is not linear but turns over gently for large $d_{\rm c}$. An analytic
prediction by Mo & White (1996) overpredicts $r_0$. The results from the
analysis of the APM cluster data by Croft et al. (1997) are nicely matched by
the $\Lambda$CDM model.
[163]
oai:arXiv.org:astro-ph/9805078 [pdf] - 101299
Peculiar Velocities of Galaxy Clusters
Submitted: 1998-05-06
We investigate the peculiar velocities predicted for galaxy clusters by
theories in the cold dark matter family. A widely used hypothesis identifies
rich clusters with high peaks of a suitably smoothed version of the linear
density fluctuation field. Their peculiar velocities are then obtained by
extrapolating the similarly smoothed linear peculiar velocities at the
positions of these peaks. We test these ideas using large high resolution
N-body simulations carried out within the Virgo supercomputing consortium. We
find that at early times the barycentre of the material which ends up in a rich
cluster is generally very close to a high peak of the initial density field.
Furthermore the mean peculiar velocity of this material agrees well with the
linear value at the peak. The late-time growth of peculiar velocities is,
however, systematically underestimated by linear theory. At the time clusters
are identified we find their rms peculiar velocity to be about 40% larger than
predicted. Nonlinear effects are particularly important in superclusters. These
systematics must be borne in mind when using cluster peculiar velocities to
estimate the parameter combination $\sigma_8\Omega^{0.6}$.
[164]
oai:arXiv.org:astro-ph/9709010 [pdf] - 554701
Evolution of structure in cold dark matter universes
Consortium, The Virgo;
:;
Jenkins, A.;
Frenk, C. S.;
Pearce, F. R.;
Thomas, P. A.;
Colberg, J. M.;
White, S. D. M.;
Couchman, H. M. P.;
Peacock, J. A.;
Efstathiou, G.;
Nelson, A. H.
Submitted: 1997-09-02, last modified: 1998-01-09
We present an analysis of the clustering evolution of dark matter in four
cold dark matter (CDM) cosmologies. We use a suite of high resolution,
17-million particle, N-body simulations which sample volumes large enough to
give clustering statistics with unprecedented accuracy. We investigate both a
flat and an open model with Omega_0=0.3, and two models with Omega=1, one with
the standard CDM power spectrum and the other with the same power spectrum as
the Omega_0=0.3 models. The amplitude of primordial fluctuations is set so that
the models reproduce the observed abundance of rich galaxy clusters by the
present day. The mass 2-point correlation function and power spectrum of all
the simulations differ significantly from those of the observed galaxy
distribution, in both shape and amplitude. Thus, for any of these models to
provide an acceptable representation of reality, the distribution of galaxies
must be biased relative to the mass in a non-trivial, scale-dependent, fashion.
In the Omega=1 models the required bias is always greater than unity, but in
the Omega_0=0.3 models an "antibias" is required on scales smaller than \sim
5\hmpc. The mass correlation functions in the simulations are well fit by
recently published analytic models. The velocity fields are remarkably similar
in all the models, whether they be characterised as bulk flows, 1-particle or
pairwise velocity dispersions. This similarity is a direct consequence of our
adopted normalisation. The small-scale pairwise velocity dispersion of the dark
matter is somewhat larger than recent determinations from galaxy redshift
surveys, but the bulk-flows predicted by our models are broadly in agreement
with most available data.
[165]
oai:arXiv.org:astro-ph/9711040 [pdf] - 99173
Linking Cluster Formation to Large Scale Structure
Submitted: 1997-11-05
We use two high resolution CDM simulations to show that (i) when clusters of
galaxies form the infall pattern of matter is not random but shows clear
features which are correlated in time; (ii) in addition, the infall patterns
are correlated with the cluster's surrounding Large Scale Structure; (iii)
Large Scale Structure shows a mix of both filaments and sheets; (iv) the amount
of mass in filaments is slightly larger for a low $\Omega$ model.
[166]
oai:arXiv.org:astro-ph/9710368 [pdf] - 99127
Genus Statistics of the Virgo N-body simulations and the 1.2-Jy Redshift
Survey
Springel, V.;
White, S. D. M.;
Colberg, J. M.;
Couchman, H. M. P.;
Efstathiou, G. P.;
Frenk, C. S.;
Jenkins, A. R.;
Pearce, F. R.;
Nelson, A. H.;
Peacock, J. A.;
Thomas, P. A.
Submitted: 1997-10-31
We study the topology of the Virgo N-body simulations and compare it to the
1.2-Jy redshift survey of IRAS galaxies by means of the genus statistic. Four
high-resolution simulations of variants of the CDM cosmology are considered: a
flat standard model (SCDM), a variant of it with more large-scale power (tCDM),
and two low density universes, one open (OCDM) and one flat (LCDM). The fully
sampled N-body simulations are examined down to strongly nonlinear scales, both
with spatially fixed smoothing, and with an adaptive smoothing technique. While
the tCDM, LCDM, and OCDM simulations have very similar genus statistics in the
regime accessible to fixed smoothing, they can be separated with adaptive
smoothing at small mass scales. In order to compare the N-body models with the
1.2-Jy survey, we extract large ensembles of mock catalogues from the
simulations. These mock surveys are used to test for systematic effects in the
genus analysis and to establish the distribution of errors of the genus curve.
We find that a simple multivariate analysis of the genus measurements is
compromised both by non-Gaussian distributed errors and by noise that dominates
the covariance matrix. We therefore introduce a principal components analysis
of the genus curve. With a likelihood ratio test we find that the 1.2-Jy data
favours the LCDM, tCDM and OCDM models compared to SCDM. When genus
measurements for different smoothing scales are combined, the SCDM model can be
excluded at a 99% confidence level, while the other three models fit the 1.2-Jy
data well. (abridged)
[167]
oai:arXiv.org:astro-ph/9709229 [pdf] - 98682
On Arcs and Omega
Submitted: 1997-09-24
The gravitational lens effect of galaxy clusters can produce large arcs from
source galaxies in their background. Typical source redshifts of ~ 1 require
clusters at z ~ 0.3 for arcs to form efficiently. Given the cluster abundance
at the present epoch, the fewer clusters exist at z ~ 0.3 the higher Omega_0
is, because the formation epoch of galaxy clusters strongly depends on Omega_0.
In addition, at fixed Omega_0, clusters are less concentrated, and hence less
efficient lenses, when the cosmological constant is positive, Omega_Lambda > 0.
Numerical cluster simulations show that the expected number of arcs on the sky
is indeed a sensitive function of Omega_0 and Omega_Lambda. The numerical
results are compatible with the statistics of observed arcs only in a universe
with low matter density, Omega_0 ~ 0.3, and zero cosmological constant. Other
models fail by one or two orders of magnitude, rendering arc statistics a
sensitive probe for cosmological parameters.
[168]
oai:arXiv.org:astro-ph/9707167 [pdf] - 97991
Arc statistics with realistic cluster potentials. IV. Clusters in
different cosmologies
Submitted: 1997-07-15
We use numerical simulations of galaxy clusters in different cosmologies to
study their ability to form large arcs. The cosmological models are: Standard
CDM (SCDM; Omega_0=1, Omega_Lambda=0); tauCDM with reduced small-scale power
(parameters as SCDM, but with a smaller shape parameter of the power spectrum);
open CDM (OCDM; Omega_0=0.3, Omega_Lambda=0); and spatially flat, low-density
CDM (LambdaCDM; Omega_0=0.3, Omega_Lambda=0.7). All models are normalised to
the local number density of rich clusters. Simulating gravitational lensing by
these clusters, we compute optical depths for the formation of large arcs. For
large arcs with length-to-width ratio >= 10, the optical depth is largest for
OCDM. Relative to OCDM, the optical depth is lower by about an order of
magnitude for LambdaCDM, and by about two orders of magnitude for S/tauCDM.
These differences originate from the different epochs of cluster formation
across the cosmological models, and from the non-linearity of the strong
lensing effect. We conclude that only the OCDM model can reproduce the observed
arc abundance well, while the other models fail to do so by orders of
magnitude.
[169]
oai:arXiv.org:astro-ph/9707018 [pdf] - 554691
The structure of galaxy clusters in different cosmologies
Consortium, Virgo;
:;
Thomas, Peter A.;
Colberg, J. M.;
Couchman, H. M. P.;
Efstathiou, G. P.;
Frenk, C. S.;
Jenkins, A. R.;
Nelson, A. H.;
Hutchings, R. M.;
Peacock, J. A.;
Pearce, F. R.;
White, S. D. M.
Submitted: 1997-07-02
We investigate the internal structure of clusters of galaxies in
high-resolution N-body simulations of 4 different cosmologies. There is a
higher proportion of disordered clusters in critical-density than in
low-density universes, although the structure of relaxed clusters is very
similar in each. Crude measures of substructure, such as the shift in the
position of the centre-of-mass as the density threshold is varied, can
distinguish the two in a sample of just 20 or so clusters; it is harder to
differentiate between clusters in open and flat models with the same density
parameter. Most clusters are in a quasi-steady state within the virial radius
and are well-described by the density profile of Navarro, Frenk & White (1995).
[170]
oai:arXiv.org:astro-ph/9702086 [pdf] - 96605
The Virgo Consortium: The evolution and formation of galaxy clusters
Colberg, J. M.;
White, S. D. M.;
Jenkins, A.;
Pearce, F. R.;
Frenk, C. S.;
Thomas, P. A.;
Hutchings, R.;
Couchman, H. M. P.;
Peacock, J. A.;
Efstathiou, G. P.;
Nelson, A. H.
Submitted: 1997-02-10
We report on work done by the Virgo consortium, an international
collaboration set up in order to study the formation and evolution of Large
Scale Structure using N-body simulations on the latest generation of parallel
supercomputers. We show results of 256^3 particle simulations of the formation
of clusters in four Dark Matter models with different cosmological parameters.
Normalizing the models such that one obtains the correct abundance of rich
clusters yields an interesting result: The peculiar velocities of the clusters
are almost independent of \Omega, and depend only weakly on \Gamma, the shape
parameter of the power spectrum. Thus, it is nearly impossible to distinguish
between high and low \Omega models on the basis of the peculiar velocities.
[171]
oai:arXiv.org:astro-ph/9610206 [pdf] - 554677
The Virgo consortium: simulations of dark matter and galaxy clustering
Consortium, The Virgo;
:;
Jenkins, A.;
Frenk, C. S.;
Pearce, F. R.;
Thomas, P. A.;
Hutchings, R.;
Colberg, J. M.;
White, S. D. M.;
Couchman, H. M. P.;
Peacock, J. A.;
Efstathiou, G. P.;
Nelson, A. H.
Submitted: 1996-10-25
We report on work in progress by the Virgo consortium, a collaboration set up
to carry out large simulations of the formation of galaxies and large-scale
structure exploiting the latest generation of parallel supercomputers. We show
results of $256^3$ particle N-body simulations of the clustering evolution of
dark matter in four cold dark matter models with different cosmological
parameters. The high resolution and large volume of these simulations allows us
to determine reliably the mass autocorrelation function for pair separations in
the range $40\hkpc$ to $20\hmpc$. Comparison of these with the observed galaxy
correlation function shows that for any of these models to be viable, the
distribution of galaxies must be biased relative to the distribution of mass in
a non-trivial, scale-dependent fashion. In particular, low $\Omega_0$ models
require the galaxies to be more ``weakly'' clustered than the mass at small and
intermediate pair separations. Simulations which include the evolution of gas
show that cold gas knots form with approximately the abundance expected on
theoretical grounds, although a few excessively massive objects grow near the
centres of rich clusters. The locations where these cold gas knots form are, in
general, biased relative to the distribution of mass in a scale-dependent way.
Some of these biases have the required sign but they are, for the most part,
weaker than is necessary for agreement with observations. The antibias present
in our low $\Omega_0$ N-body/SPH simulation appears to be related to the
merging and disruption of galaxies in rich clusters.
