BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ben Moore

Ben Moore appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-phastro-ph.EPhep-phastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

47works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

47 published item(s)

preprint2016arXiv

Close-in planetesimal formation by pile-up of drifting pebbles

The consistency of planet formation models suffers from the disconnection between the regime of small and large bodies. This is primarily caused by so-called growth barriers: the direct growth of larger bodies is halted at centimetre-sized objects and particular conditions are required for the formation of larger, gravitationally bound planetesimals. We aim to connect models of dust evolution and planetesimal formation to identify regions of protoplanetary discs that are favourable for the formation of kilometre-sized bodies and the first planetary embryos. We combine semi-analytical models of viscous protoplanetary disc evolution, dust growth and drift including backreaction of the dust particles on the gas, and planetesimal formation via the streaming instability into one numerical code. We investigate how planetesimal formation is affected by the mass of the protoplanetary disc, its initial dust content, and the stickiness of dust aggregates. We find that the dust growth and drift leads to a global redistribution of solids. The pile-up of pebbles in the inner disc provides local conditions where the streaming instability is effective. Planetesimals form in an annulus with its inner edge lying between 0.3 AU and 1 AU and its width ranging from 0.3 AU to 3 AU. The resulting surface density of planetesimals follows a radial profile that is much steeper than the initial disc profile. These results support formation of terrestrial planets in the solar system from a narrow annulus of planetesimals, which reproduces their peculiar mass ratios.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Internal Alignments of Red Versus Blue Discs in Dark Matter Halos

Large surveys have shown that red galaxies are preferentially aligned with their halos while blue galaxies have a more isotropic distribution. Since halos generally align with their filaments this introduces a bias in the measurement of the cosmic shear from weak lensing. It is therefore vitally important to understand why this difference arises. We explore the stability of different disc orientations within triaxial halos. We show that, in the absence of gas, the disc orientation is most stable when its spin is along the minor axis of the halo. Instead when gas cools onto a disc it is able to form in almost arbitrary orientation, including off the main planes of the halo (but avoiding an orientation perpendicular to the halo's intermediate axis). Substructure helps gasless galaxies reach alignment with the halo faster, but have less effect on galaxies when gas is cooling onto the disc. Our results provide a novel and natural interpretation for why red, gas poor galaxies are preferentially aligned with their halo, while blue, star-forming, galaxies have nearly random orientations, without requiring a connection between galaxies' current star formation rate and their merger history.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The same with less: The cosmic web of warm versus cold dark matter dwarf galaxies

We explore fundamental properties of the distribution of low mass dark matter halos within the cosmic web using warm dark matter (WDM) and cold dark matter (CDM) cosmological simulations. Using self abundance-matched mock galaxy catalogs, we show that the distribution of dwarf galaxies in a WDM universe, wherein low mass halo formation is heavily suppressed, is nearly indistinguishable to that of a CDM universe whose low mass halos are not seen because galaxy formation is suppressed below some threshold halo mass. However, if the scatter between dwarf galaxy luminosity and halo properties is large enough, low mass CDM halos would sometimes host relatively bright galaxies thereby populating CDM voids with the occasional isolated galaxy and reducing the numbers of completely empty voids. Otherwise, without high mass to light scatter, all mock galaxy clustering statistics that we consider--the auto-correlation function, the numbers and radial profiles of satellites, the numbers of isolated galaxies, and the PDF of small voids--are nearly identical in CDM and WDM. WDM voids are neither larger nor emptier than CDM voids, when constructed from abundance-matched halo catalogs. It is thus a challenge to determine whether the CDM problem of the over-abundance of small halos with respect to the number density of observed dwarf galaxies has a cosmological solution or an astrophysical solution. However, some clues about the dark matter particle and the scatter between the properties of dwarf galaxies and their dark matter halo hosts might be found in the cosmic web of galaxies in future surveys of the local volume.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Globular Cluster Toward M87 with a Radial Velocity < -1000 km/s: The First Hypervelocity Cluster

We report the discovery of an object near M87 in the Virgo Cluster with an extraordinary blueshift of -1025 km/s, offset from the systemic velocity by >2300 km/s. Evaluation of photometric and spectroscopic data provides strong evidence that this object is a distant massive globular cluster, which we call HVGC-1 in analogy to Galactic hypervelocity stars. We consider but disfavor more exotic interpretations, such as a system of stars bound to a recoiling black hole. The odds of observing an outlier as extreme as HVGC-1 in a virialized distribution of intracluster objects are small; it appears more likely that the cluster was (or is being) ejected from Virgo following a three-body interaction. The nature of the interaction is unclear, and could involve either a subhalo or a binary supermassive black hole at the center of M87.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Systematic Look at the Effects of Radiative Feedback on Disc Galaxy Formation

Galaxy formation models and simulations rely on various feedback mechanisms to reproduce the observed baryonic scaling relations and galaxy morphologies. Although dwarf galaxy and giant elliptical properties can be explained using feedback from supernova and active galactic nuclei, Milky Way-sized galaxies still represent a challenge to current theories of galaxy formation. In this paper, we explore the possible role of feedback from stellar radiation in regulating the main properties of disk galaxies such as our own Milky Way. We have performed a suite of cosmological simulations of the same $\sim10^{12} {\rm M}_{\odot}$ halo selected based on its rather typical mass accretion history. We have implemented radiative feedback from young stars using a crude model of radiative transfer for ultraviolet (UV) and infrared (IR) radiation. However, the model is realistic enough such that the dust opacity plays a direct role in regulating the efficiency of our feedback mechanism. We have explored various models for the dust opacity, assuming different constant dust temperatures, as well as a varying dust temperature model. We find that while strong radiative feedback appears as a viable mechanism to regulate the stellar mass fraction in massive galaxies, it also prevents the formation of discs with reasonable morphologies. In models with strong stellar radiation feedback, stellar discs are systematically too thick while the gas disc morphology is completely destroyed due to the efficient mixing between the feedback-affected gas and its surroundings. At the resolution of our simulation suite, we find it impossible to preserve spiral disc morphology while at the same time expelling enough baryons to satisfy the abundance matching constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Brightest Cluster Galaxies in Cosmological Simulations with Adaptive Mesh Refinement: Successes and Failures

A large sample of cosmological hydrodynamical zoom-in simulations with Adaptive Mesh Refinement (AMR) is analysed to study the properties of simulated Brightest Cluster Galaxies (BCGs). Following the formation and evolution of BCGs requires modeling an entire galaxy cluster, because the BCG properties are largely influenced by the state of the gas in the cluster and by interactions and mergers with satellites. BCG evolution is also deeply influenced by the presence of gas heating sources such as Active Galactic Nuclei (AGNs) that prevent catastrophic cooling of large amounts of gas. We show that AGN feedback is one of the most important mechanisms in shaping the properties of BCGs at low redshift by analysing our statistical sample of simulations with and without AGN feedback. When AGN feedback is included BCG masses, sizes, star formation rates and kinematic properties are closer to those of the observed systems. Some small discrepancies are observed only for the most massive BCGs and in the fraction of star-forming BCGs, effects that might be due to physical processes that are not included in our model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Precision cosmology in muddy waters: Cosmological constraints and N-body codes

Future large-scale structure surveys of the Universe will aim to constrain the cosmological model and the true nature of dark energy with unprecedented accuracy. In order for these surveys to achieve their designed goals, they will require predictions for the nonlinear matter power spectrum to sub-percent accuracy. Through the use of a large ensemble of cosmological N-body simulations, we demonstrate that if we do not understand the uncertainties associated with simulating structure formation, i.e. knowledge of the `true' simulation parameters, and simply seek to marginalize over them, then the constraining power of such future surveys can be significantly reduced. However, for the parameters {n_s, h, Om_b, Om_m}, this effect can be largely mitigated by adding the information from a CMB experiment, like Planck. In contrast, for the amplitude of fluctuations sigma8 and the time-evolving equation of state of dark energy {w_0, w_a}, the mitigation is mild. On marginalizing over the simulation parameters, we find that the dark-energy figure of merit can be degraded by ~2. This is likely an optimistic assessment, since we do not take into account other important simulation parameters. A caveat is our assumption that the Hessian of the likelihood function does not vary significantly when moving from our adopted to the 'true' simulation parameter set. This paper therefore provides strong motivation for rigorous convergence testing of N-body codes to meet the future challenges of precision cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The biasing of baryons on the cluster mass function and cosmological parameter estimation

We study the effect of baryonic processes on the halo mass function in the galaxy cluster mass range using a catalogue of 153 high resolution cosmological hydrodynamical simulations performed with the AMR code ramses. We use the results of our simulations within a simple analytical model to gauge the effects of baryon physics on the halo mass function. Neglect of AGN feedback leads to a significant boost in the cluster mass function similar to that reported by other authors. However, including AGN feedback not only gives rise to systems that are similar to observed galaxy clusters, but they also reverse the global baryonic effects on the clusters. The resulting mass function is closer to the unmodified dark matter halo mass function but still contains a mass dependent bias at the 5-10% level. These effects bias measurements of the cosmological parameters, such as $σ_8$ and $Ω_m$. For current cluster surveys baryonic effects are within the noise for current survey volumes, but forthcoming and planned large SZ, X-ray and multi-wavelength surveys will be biased at the percent level by these processes. The predictions for the halo mass function including baryonic effects need to be carefully studied with larger and improved simulations. However, simulations of full cosmological boxes with the resolution we achieve and including AGN feedback are still computationally challenging.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Toward an accurate mass function for precision cosmology

Cosmological surveys aim to use the evolution of the abundance of galaxy clusters to accurately constrain the cosmological model. In the context of LCDM, we show that it is possible to achieve the required percent level accuracy in the halo mass function with gravity-only cosmological simulations, and we provide simulation start and run parameter guidelines for doing so. Some previous works have had sufficient statistical precision, but lacked robust verification of absolute accuracy. Convergence tests of the mass function with, for example, simulation start redshift can exhibit false convergence of the mass function due to counteracting errors, potentially misleading one to infer overly optimistic estimations of simulation accuracy. Percent level accuracy is possible if initial condition particle mapping uses second order Lagrangian Perturbation Theory, and if the start epoch is between 10 and 50 expansion factors before the epoch of halo formation of interest. The mass function for halos with fewer than ~1000 particles is highly sensitive to simulation parameters and start redshift, implying a practical minimum mass resolution limit due to mass discreteness. The narrow range in converged start redshift suggests that it is not presently possible for a single simulation to capture accurately the cluster mass function while also starting early enough to model accurately the numbers of reionisation era galaxies, whose baryon feedback processes may affect later cluster properties. Ultimately, to fully exploit current and future cosmological surveys will require accurate modeling of baryon physics and observable properties, a formidable challenge for which accurate gravity-only simulations are just an initial step.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

What's Up in the Milky Way? The Orientation of the Disc Relative to the Triaxial Halo

Models of the Sagittarius Stream have consistently found that the Milky Way disc is oriented such that its short axis is along the intermediate axis of the triaxial dark matter halo. We attempt to build models of disc galaxies in such an `intermediate-axis orientation'. We do this with three models. In the first two cases we simply rigidly grow a disc in a triaxial halo such that the disc ends up perpendicular to the global intermediate axis. We also attempt to coax a disc to form in an intermediate-axis orientation by producing a gas+dark matter triaxial system with gas angular momentum about the intermediate axis. In all cases we fail to produce systems which remain with stellar angular momentum aligned with the halo's intermediate axis, even when the disc's potential flattens the inner halo such that the disc is everywhere perpendicular to the halo's local minor axis. For one of these unstable simulations we show that the potential is even rounder than the models of the Milky Way potential in the region probed by the Sagittarius Stream. We conclude that the Milky Way's disc is very unlikely to be in an intermediate axis orientation. However we find that a disc can persist off one of the principal planes of the potential. We propose that the disc of the Milky Way must be tilted relative to the principal axes of the dark matter halo. Direct confirmation of this prediction would constitute a critical test of Modified Newtonian Dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cores in warm dark matter haloes: a Catch 22 problem

The free streaming of warm dark matter particles dampens the fluctuation spectrum, flattens the mass function of haloes and imprints a fine grained phase density limit for dark matter structures. The phase space density limit is expected to imprint a constant density core at the halo center on the contrary to what happens for cold dark matter. We explore these effects using high resolution simulations of structure formation in different warm dark matter scenarios. We find that the size of the core we obtain in simulated haloes is in good agreement with theoretical expectations based on Liouville's theorem. However, our simulations show that in order to create a significant core, (r_c~1 kpc), in a dwarf galaxy (M~1e10 Msun), a thermal candidate with a mass as low as 0.1 keV is required. This would fully prevent the formation of the dwarf galaxy in the first place. For candidates satisfying large scale structure constrains (m_wdm larger than 1-2 keV) the expected size of the core is of the order of 10 (20) pc for a dark matter halo with a mass of 1e10 (1e8) Msun. We conclude that "standard" warm dark matter is not viable solution for explaining the presence of cored density profiles in low mass galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evidence for inhomogeneous reionization in the local Universe from metal-poor globular cluster systems

Exploiting a fundamental characteristic of galaxy assembly in the ΛCDM paradigm, the observed spatial biasing and kinematics of metal-poor globular star clusters are used to constrain the local reionization epoch around individual galaxies. Selecting three galaxies located in different environments, the first attempt at constraining the environmental propagation of reionization in the local Universe is carried out. The joint constraint from the three galaxies (z_reion = 10.5^{+1.0}_{-0.9}) agrees remarkably well with the latest WMAP constraint on z_reion for a simple instantaneous reionization model. More importantly, the range of z_reion values found here are consistent with the global range of z_reion estimates from other observations. We furthermore find a 1.7σ indication that reionization completed in low-density environments before the intergalactic medium in high-density environments was reionized. This is consistent with certain theoretical models that predict that reionization was globally prolonged in duration, with neutral hydrogen pockets surviving in high-density environments, even after the surrounding regions were reionized. More generally, this work provides a useful constraint on the formation history of galaxy stellar halos.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

LEDA 074886: A remarkable rectangular-looking galaxy

We report the discovery of an interesting and rare, rectangular-shaped galaxy. At a distance of 21 Mpc, the dwarf galaxy LEDA 074886 has an absolute R-band magnitude of -17.3 mag. Adding to this galaxy's intrigue is the presence of an embedded, edge-on stellar disk (of extent 2R_{e,disk} = 12 arcsec = 1.2 kpc) for which Forbes et al. reported V_rot/sigma ~ 1.4. We speculate that this galaxy may be the remnant of two (nearly edge-one) merged disk galaxies in which the initial gas was driven inward and subsequently formed the inner disk, while the stars at larger radii effectively experienced a dissipationless merger event resulting in this `emerald cut galaxy' having very boxy isophotes with a_4/a = -0.05 to -0.08 from 3 to 5 kpc. This galaxy suggests that knowledge from simulations of both `wet' and `dry' galaxy mergers may need to be combined to properly understand the various paths that galaxy evolution can take, with a particular relevance to blue elliptical galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear Evolution of Cosmological Structures in Warm Dark Matter Models

The dark energy dominated warm dark matter (WDM) model is a promising alternative cosmological scenario. We explore large-scale structure formation in this paradigm. We do this in two different ways: with the halo model approach and with the help of an ensemble of high resolution N-body simulations. Combining these quasi-independent approaches, leads to a physical understanding of the important processes which shape the formation of structures. We take a detailed look at the halo mass function, the concentrations and the linear halo bias of WDM. In all cases we find interesting deviations with respect to CDM. In particular, the concentration-mass relation displays a turnover for group scale dark matter haloes, for the case of WDM particles with masses of the order ~0.25 keV. This may be interpreted as a hint for top-down structure formation on small scales. We implement our results into the halo model and find much better agreement with simulations. On small scales the WDM halo model now performs as well as its CDM counterpart.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the origin of elemental abundances in the terrestrial planets

The abundances of elements in the Earth and the terrestrial planets provide the initial conditions for life and clues as to the history and formation of the Solar System. We follow the pioneering work of Bond et al. (2010) and combine circumstellar disk models, chemical equilibrium calculations and dynamical simulations of planet formation to study the bulk composition of rocky planets. We use condensation sequence calculations to estimate the initial abundance of solids in the circumstellar disk with properties determined from time dependent theoretical models. We combine this with dynamical simulations of planetesimal growth that trace the solids during the planet formation process. We calculate the elemental abundances in the resulting planets and explore how these vary with the choice of disk model and the initial conditions within the Solar Nebula. Although certain characteristics of the terrestrial planets in the Solar System could be reproduced, none of our models could reproduce the abundance properties of all the planets. We found that the choice of the initial planetesimal disk mass and of the disk model has a significant effect on composition gradients. Disk models that give higher pressure and temperature result in larger variations in the bulk chemical compositions of the resulting planets due to inhomogeneities in the element abundance profiles and due to the different source regions of the planets in the dynamical simulations. We observed a trend that massive planets and planets with relatively small semi-major axes are more sensitive to these variations than smaller planets at larger radial distance. Only these large variations in the simulated chemical abundances can potentially explain the diverse bulk composition of the Solar System planets, whereas Mercury's bulk composition can not be reproduced in our approach.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The effects of baryon physics, black holes and AGN feedback on the mass distribution in clusters of galaxies

The spatial distribution of matter in clusters of galaxies is mainly determined by the dominant dark matter component, however, physical processes involving baryonic matter are able to modify it significantly. We analyse a set of 500 pc resolution cosmological simulations of a cluster of galaxies with mass comparable to Virgo, performed with the AMR code RAMSES. We compare the mass density profiles of the dark, stellar and gaseous matter components of the cluster that result from different assumptions for the subgrid baryonic physics and galaxy formation processes. First, the prediction of a gravity only N-body simulation is compared to that of a hydrodynamical simulation with standard galaxy formation recipes, then all results are compared to a hydrodynamical simulation which includes thermal AGN feedback from Super Massive Black Holes (SMBH). We find the usual effects of overcooling and adiabatic contraction in the run with standard galaxy formation physics, but very different results are found when implementing SMBHs and AGN feedback. Star formation is strongly quenched, producing lower stellar densities throughout the cluster, and much less cold gas is available for star formation at low redshifts. At redshift z = 0 we find a flat density core of radius 10 kpc in both of the dark and stellar matter density profiles. We specu- late on the possible formation mechanisms able to produce such cores and we conclude that they can be produced through the coupling of different processes: (I) dynamical friction from the decay of black hole orbits during galaxy mergers; (II) AGN driven gas outflows producing fluctuations of the gravitational potential causing the removal of collisionless matter from the central region of the cluster; (III) adiabatic expansion in response to the slow expulsion of gas from the central region of the cluster during the quiescent mode of AGN activity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The formation of the brightest cluster galaxies in cosmological simulations: the case for AGN feedback

We use 500 pc resolution cosmological simulations of a Virgo-like galaxy cluster to study the properties of the brightest cluster galaxy (BCG) that forms at the center of the halo. We compared two simulations; one incorporating only supernovae feedback and a second that also includes prescriptions for black hole growth and the resulting AGN feedback from gas accretion. As previous work has shown, with supernovae feedback alone we are unable to reproduce any of the observed properties of massive cluster ellipticals. The resulting BCG is rotating quickly, has a high Sersic index, a strong mass excess in the center and a total central density profile falling more steeply than isothermal. Furthermore, it is far too efficient at converting most of the available baryons into stars which is strongly constrained by abundance matching. With a treatment of black hole dynamics and AGN feedback the BCG properties are in good agreement with data: they rotate slowly, have a cored surface density profile, a flat or rising velocity dispersion profile and a low stellar mass fraction. The AGN provides a new mechanism to create cores in luminous elliptical galaxies; the core expands due to the combined effects of heating from dynamical friction of sinking massive black holes and AGN feedback that ejects gaseous material from the central regions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

How common are Earth-Moon planetary systems?

The Earth's comparatively massive moon, formed via a giant impact on the proto-Earth, has played an important role in the development of life on our planet, both in the history and strength of the ocean tides and in stabilizing the chaotic spin of our planet. Here we show that massive moons orbiting terrestrial planets are not rare. A large set of simulations by Morishima et al., 2010, where Earth-like planets in the habitable zone form, provides the raw simulation data for our study. We use limits on the collision parameters that may guarantee the formation of a circumplanetary disk after a protoplanet collision that could form a satellite and study the collision history and the long term evolution of the satellites qualitatively. In addition, we estimate and quantify the uncertainties in each step of our study. We find that giant impacts with the required energy and orbital parameters for producing a binary planetary system do occur with more than 1 in 12 terrestrial planets hosting a massive moon, with a low-end estimate of 1 in 45 and a high-end estimate of 1 in 4.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Microhalos and Dark Matter Detection

Cosmological structure formation predicts that our galactic halo contains an enormous hierarchy of substructures and streams, the remnants of the merging hierarchy that began with tiny Earth mass microhalos. If these structures persist until the present time, they could influence dramatically the detection signatures of weakly interacting elementary particle dark matter (WIMP). Using numerical simulations that follow the tidal disruption within the Galactic potential and heating from stellar encounters, we find that neither microhalos nor streams have significant impact on direct detection, implying that dark matter constraints derived using simple smooth halo models are relatively robust. We also find that many dense central cusps survive, yielding a small enhancement in the signal for indirect detection experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multi-Wavelength Properties of Barred Galaxies in the Local Universe: Environment and evolution across the Hubble sequence

We investigate possible environmental and morphological trends in the $z\sim0$ bar fraction using two carefully selected samples representative of a low-density environment (the isolated galaxies from the AMIGA sample) and of a dense environment (galaxies in the Virgo cluster). Galaxies span a stellar mass range from $10^8$ to $10^{12}$M$_{\odot}$ and are visually classified using both high-resolution NIR (H-band) imaging and optical \texttt{rgb} images. We find that the bar fraction in disk galaxies is independent of environment suggesting that bar formation may occur prior to the formation of galaxy clusters. The bar fraction in early type spirals ($Sa-Sb$) is $\sim$50%, which is twice as high as the late type spirals ($Sbc-Sm$). The higher bar fraction in early type spirals may be due to the fact that a significant fraction of their bulges are pseudo-bulges which form via the buckling instability of a bar. i.e. a large part of the Hubble sequence is due to secular processes which move disc galaxies from late to early types. There is a hint of a higher bar fraction with higher stellar masses which may be due to the susceptibility to bar instabilities as the baryon fractions increase in halos of larger masses. Overall, the $S0$ population has a lower bar fraction than the $Sa-Sb$ galaxies and their barred fraction drops significantly with decreasing stellar mass. This supports the notion that $S0s$ form via the transformation of disk galaxies that enter the cluster environment. The gravitational harassment thickens the stellar disks, wiping out spiral patterns and eventually erasing the bar - a process that is more effective at lower galaxy masses.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the Stability of Tidal Streams

We explore the stability of tidal streams to perturbations, motivated by recent claims that the clumpy structure of the stellar streams surrounding the globular cluster Palomar 5 are the result of gravitational instability. We calculate the Jeans length of tidal streams by treating them as a thin expanding cylinder of collisionless matter. We also find a general relation between the density and the velocity dispersion inside a stream, which is used to determine the longitudinal Jeans criterion. Our analytic results are checked by following the time evolution of the phase space density within streams using numerical simulations. We conclude that tidal streams within our galactic halo are stable on all length scales and over all timescales.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Core creation in galaxies and haloes via sinking massive objects

We perform a detailed investigation into the disruption of central cusps via the transfer of energy from sinking massive objects. Constant density inner regions form at the radius where the enclosed mass approximately matches the mass of the infalling body. We explore parameter space using numerical simulations and give an empirical relation for the size of the resulting core within structures that have different initial cusp slopes. We find that infalling bodies always stall at the edge of these newly formed cores, experiencing no dynamical friction over many dynamical times. As applications, we consider the resulting decrease in the dark matter annihilation flux due to centrally destroyed cusps, and we present a new theory for the formation of close binary nuclei -- the `stalled binary' model. We focus on one particularly interesting binary nucleus system, the dwarf spheroidal galaxy VCC 128 which is dark matter dominated at all radii. We show that its nuclei would rapidly coalesce within a few million years if it has a central dark matter cusp slope steeper than r^{-1}. However, if its initial dark matter cusp is slightly shallower than a log slope of -0.75 at ~0.1% of the virial radius, then the sinking nuclei naturally create a core equal to their observed separation and stall. This is close to the log slope measured a recent billion particle CDM halo simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets

We present results from a suite of N-body simulations that follow the accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5--4.0 AU and the collisional growth is followed until the completion of planetary accretion (> 100 Myr). All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas, and orbits of Jupiter and Saturn. We end up with one to five planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today's eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Impact of Dark Matter Microhalos on Signatures for Direct and Indirect Detection

Detecting dark matter as it streams through detectors on Earth relies on knowledge of its phase space density on a scale comparable to the size of our solar system. Numerical simulations predict that our Galactic halo contains an enormous hierarchy of substructures, streams and caustics, the remnants of the merging hierarchy that began with tiny Earth mass microhalos. If these bound or coherent structures persist until the present time, they could dramatically alter signatures for the detection of weakly interacting elementary particle dark matter (WIMP). Using numerical simulations that follow the coarse grained tidal disruption within the Galactic potential and fine grained heating from stellar encounters, we find that microhalos, streams and caustics have a negligible likelihood of impacting direct detection signatures implying that dark matter constraints derived using simple smooth halo models are relatively robust. We also find that many dense central cusps survive, yielding a small enhancement in the signal for indirect detection experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mass Distribution in Galaxy Clusters: the Role of AGN Feedback

We use 1 kpc resolution cosmological AMR simulations of a Virgo-like galaxy cluster to investigate the effect of feedback from supermassive black holes (SMBH) on the mass distribution of dark matter, gas and stars. We compared three different models: (i) a standard galaxy formation model featuring gas cooling, star formation and supernovae feedback, (ii) a "quenching" model for which star formation is artificially suppressed in massive halos and finally (iii) the recently proposed AGN feedback model of Booth & Schaye (2009). Without AGN feedback (even in the quenching case), our simulated cluster suffers from a strong overcooling problem, with a stellar mass fraction significantly above observed values in M87. The baryon distribution is highly concentrated, resulting in a strong adiabatic contraction (AC) of dark matter. With AGN feedback, on the contrary, the stellar mass in the bright central galaxy (BCG) lies below observational estimates and the overcooling problem disappears. The stellar mass of the BCG is seen to increase with increasing mass resolution, suggesting that our stellar masses converges to the correct value from below. The gas and total mass distributions are in striking agreement with observations. We also find a slight deficit (~10%) of baryons at the virial radius, due to the effect of AGN-driven shock waves pushing gas to Mpc scales and beyond. This baryon deficit results in a slight adiabatic expansion of the dark matter distribution, that can be explained quantitatively by AC theory.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Multi-Wavelength Properties of Barred Galaxies in the Local Universe. I: Virgo Cluster

We study in detail how the barred galaxy fraction varies as a function of luminosity, HI gas mass, morphology and color in the Virgo cluster in order to provide a well defined, statistically robust measurement of the bar fraction in the local universe spanning a wide range in luminosity (factor of ~100) and HI gas mass. We combine multiple public data-sets (UKIDSS near-infrared imaging, ALFALFA HI gas masses, GOLDMine photometry). After excluding highly inclined systems, we define three samples where galaxies are selected by their B-band luminosity, H-band luminosity, and HI gas mass. We visually assign bars using the high resolution H-band imaging from UKIDSS. When all morphologies are included, the barred fraction is ~17-24% while for morphologically selected discs, we find that the barred fraction in Virgo is ~29-34%: it does not depend strongly on how the sample is defined and does not show variations with luminosity or HI gas mass. The barred fraction depends most strongly on the morphological composition of the sample: when the disc populations are separated into lenticulars (S0--S0/a), early-type spirals (Sa--Sb), and late-type spirals (Sbc--Sm), we find that the early-type spirals have a higher barred fraction (~45-50%) compared to the lenticulars and late-type spirals (~22-36%). This difference may be due to the higher baryon fraction of early-type discs which makes them more susceptible to bar instabilities. We do not find any evidence of barred galaxies being preferentially blue.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

New tools for probing the phase space structure of dark matter halos

We summarize recent developments in the use of spectral methods for analyzing large numbers of orbits in N-body simulations to obtain insights into the global phase space structure of dark matter halos. The fundamental frequencies of oscillation of orbits can be used to understand the physical mechanism by which the shapes of dark matter halos evolve in response to the growth of central baryonic components. Halos change shape primarily because individual orbits change their shapes adiabatically in response to the growth of a baryonic component, with those at small radii become preferentially rounder. Chaotic scattering of orbits occurs only when the central point mass is very compact and is equally effective for centrophobic long-axis tube orbits as it is for centrophilic box orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reionization of the Local Group of Galaxies

We present the first detailed structure formation and radiative transfer simulations of the reionization history of our cosmic neighbourhood. To this end, we follow the formation of the Local Group of galaxies and nearby clusters by means of constrained simulations, which use the available observational constraints to construct a representation of those structures which reproduces their actual positions and properties at the present time. We find that the reionization history of the Local Group is strongly dependent on the assumed photon production efficiencies of the ionizing sources, which are still poorly constrained. If sources are relatively efficient, i.e. the process is 'photon-rich', the Local Group is primarily ionized externally by the nearby clusters. Alternatively, if the sources are inefficient, i.e. reionization is 'photon-poor' the Local Group evolves largely isolated and reionizes itself. The mode of reionization, external vs. internal, has important implications for the evolution of our neighbourhood, in terms of e.g. its satellite galaxy populations and primordial stellar populations. This therefore provides an important avenue for understanding the young universe by detailed studies of our nearby structures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Systematic uncertainties in the determination of the local dark matter density

A precise determination of the local dark matter density and an accurate control over the corresponding uncertainties are of paramount importance for Dark Matter (DM) searches. Using very recent high-resolution numerical simulations of a Milky Way like object, we study the systematic uncertainties that affect the determination of the local dark matter density based on dynamical measurements in the Galaxy. In particular, extracting from the simulation with baryons the orientation of the Galactic stellar disk with respect to the DM distribution, we study the DM density for an observer located at $\sim$8 kpc from the Galactic center {\it on the stellar disk}, $ρ_0$. This quantity is found to be always larger than the average density in a spherical shell of same radius $\barρ_0$, which is the quantity inferred from dynamical measurements in the Galaxy, and to vary in the range $ρ_0/\barρ_0=1.01-1.41$. This suggests that the actual dark matter density in the solar neighbourhood is on average 21\% larger than the value inferred from most dynamical measurements, and that the associated systematic errors are larger than the statistical errors recently discussed in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The formation of disc galaxies in a LCDM universe

We study the formation of disc galaxies in a fully cosmological framework using adaptive mesh refinement simulations. We perform an extensive parameter study of the main subgrid processes that control how gas is converted into stars and the coupled effect of supernovae feedback. We argue that previous attempts to form disc galaxies have been unsuccessful because of the universal adoption of strong feedback combined with high star formation efficiencies. Unless extreme amounts of energy are injected into the interstellar medium during supernovae events, these star formation parameters result in bulge dominated S0/Sa galaxies as star formation is too efficient at z~3. We show that a low efficiency of star-formation more closely models the subparsec physical processes, especially at high redshift. We highlight the successful formation of extended disc galaxies with scale lengths r_d=4-5 kpc, flat rotation curves and bulge to disc ratios of B/D~1/4. Not only do we resolve the formation of a Milky Way-like spiral galaxy, we also observe the secular evolution of the disc as it forms a pseudo-bulge. The disc properties agree well with observations and are compatible with the photometric and baryonic Tully-Fisher relations, the Kennicutt-Schmidt relation and the observed angular momentum content of spiral galaxies. We conclude that underlying small-scale star formation physics plays a larger role than previously considered in simulations of galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Disk formation and the origin of clumpy galaxies at high redshift

Observations of high redshift galaxies have revealed a multitude of large clumpy rapidly star-forming galaxies. Their formation scenario and their link to present day spirals is still unknown. In this Letter we perform adaptive mesh refinement simulations of disk formation in a cosmological context that are unrivalled in terms of mass and spatial resolution. We find that the so called "chain-galaxies" and "clump-clusters" are a natural outcome of early epochs of enhanced gas accretion from cold dense streams as well as tidally and ram-pressured stripped material from minor mergers and satellites. Through interaction with the hot halo gas, this freshly accreted cold gas settles into a large disk-like system, not necessarily aligned to an older stellar component, that undergoes fragmentation and subsequent star formation, forming large clumps in the mass range 10^7-10^9 M_sun. Galaxy formation is a complex process at this important epoch when most of the central baryons are being acquired through a range of different mechanisms - we highlight that a rapid mass loading epoch is required to fuel the fragmentation taking place in the massive arms in the outskirts of extended disks, an accretion mode that occurs naturally in the hierarchical assembly process at early epochs.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Effects of ram pressure on the gas distribution and star formation in the Large Magellanic Cloud

We use high resolution N-body/SPH simulations to study the hydrodynamical interaction between the Large Magellanic Cloud (LMC) and the hot halo of the Milky Way. We investigate whether ram-pressure acting on the satellite's ISM can explain the peculiarities observed in the HI distribution and the location of the recent star formation activity. Due to the present nearly edge-on orientation of the disk with respect to the orbital motion, compression at the leading edge can explain the high density region observed in HI at the south-east border. In the case of a face-on disk (according to Mastropietro et al. 2008 the LMC was moving almost face-on before the last perigalactic passage), ram-pressure directed perpendicularly to the disk produces a clumpy structure characterized by voids and high density filaments that resemble those observed by the Parkes HI survey. As a consequence of the very recent edge-on motion, the H-alpha emission is mainly concentrated on the eastern side where 30 Doradus and most of the supergiant shells are located, although some H-alpha complexes form a patchy distribution on the entire disk. In this scenario only the youngest stellar complexes show a progression in age along the leading border of the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The structure and evolution of cold dark matter halos

In the standard cosmological model a mysterious cold dark matter (CDM) component dominates the formation of structures. Numerical studies of the formation of CDM halos have produced several robust results that allow unique tests of the hierarchical clustering paradigm. Universal properties of halos, including their mass profiles and substructure properties are roughly consistent with observational data from the scales of dwarf galaxies to galaxy clusters. Resolving the fine grained structure of halos has enabled us to make predictions for ongoing and planned direct and indirect dark matter detection experiments. While simulations of pure CDM halos are now very accurate and in good agreement (recently claimed discrepancies are addressed in detail in this review), we are still unable to make robust, quantitative predictions about galaxy formation and about how the dark matter distribution changes in the process. Whilst discrepancies between observations and simulations have been the subject of much debate in the literature, galaxy formation and evolution needs to be understood in more detail in order to fully test the CDM paradigm. Whatever the true nature of the dark matter particle is, its clustering properties must not be too different from a cold neutralino like particle to maintain all the successes of the model in matching large scale structure data and the global properties of halos which are mostly in good agreement with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Quantifying the heart of darkness with GHALO - a multi-billion particle simulation of our galactic halo

We perform a series of simulations of a Galactic mass dark matter halo at different resolutions, our largest uses over three billion particles and has a mass resolution of 1000 M_sun. We quantify the structural properties of the inner dark matter distribution and study how they depend on numerical resolution. We can measure the density profile to a distance of 120 pc (0.05% of R_vir) where the logarithmic slope is -0.8 and -1.4 at (0.5% of R_vir). We propose a new two parameter fitting function that has a linearly varying logarithmic density gradient which fits the GHALO and VL2 density profiles extremely well. Convergence in the density profile and the halo shape scales as N^(-1/3), but the shape converges at a radius three times larger at which point the halo becomes more spherical due to numerical resolution. The six dimensional phase-space profile is dominated by the presence of the substructures and does not follow a power law, except in the smooth under-resolved inner few kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Concentration, Spin and Shape of Dark Matter Haloes: Scatter and the Dependence on Mass and Environment

We use a series of cosmological N-body simulations for a flat LCDM cosmology to investigate the properties of dark matter haloes in the mass range 3.0e9-3.0e13 Msun. These properties include the concentration parameter (c), the spin parameter (lambda) and the mean axis ratio (q). For the concentration-mass relation we find c~M^(-0.11) in agreement with the model proposed by Bullock et al. even if we find a lower normalization (15%). The results for lambda and q are in good agreement with previous studies, while c and lambda are anti-correlated. In an attempt to remove unrelaxed haloes, we use the offset parameter (xoff), defined as the distance between the most bound particle and the center of mass. Removing haloes with large xoff increases the c by ~10%, lowers the lambda by ~15%, and removes the most prolate haloes. In addition, it largely removes the anti-correlation between c and lambda though not entirely. We also investigate the effects of the large-scale environment. We find that more concentrated haloes live in denser environments. Note, however, that the trend is weak compared to the scatter. For the spin parameters we find no environment dependence, while there is a weak indication that the most spherical haloes reside in denser region. Finally, using a simple model for disk galaxy formation we show that haloes that host low surface brightness galaxies are expected to be hosted by a biased sub-set of haloes. Not only do these haloes have spin parameters that are larger than average, they also have c that are 15% percent lower than the average at a given halo mass. We discuss the implications of all these findings for the claimed disagreement between halo concentrations inferred from LSB rotation curves, and those expected for a LCDM cosmology. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Empirical models for Dark Matter Halos. I. Nonparametric Construction of Density Profiles and Comparison with Parametric Models

We use techniques from nonparametric function estimation theory to extract the density profiles, and their derivatives, from a set of N-body dark matter halos. We consider halos generated from LCDM simulations of gravitational clustering, as well as isolated, spherical collapses. The logarithmic density slopes gamma = d(log rho)/d(log r) of the LCDM halos are found to vary as power-laws in radius, reaching values of gamma ~ -1 at the innermost resolved radii (~0.01 r_virial). This behavior is significantly different from that of broken power-law models like the NFW profile, but similar to that of models like de Vaucouleurs'. Accordingly, we compare the N-body density profiles with various parametric models to find which provide the best fit. We consider an NFW-like model with arbitrary inner slope; Dehnen & McLaughlin's anisotropic model; Einasto's model (identical in functional form to Sersic's model but fit to the space density); and the density model of Prugniel & Simien that was designed to match the deprojected form of Sersic's R^{1/n} law. Overall, the best-fitting model to the LCDM halos is Einasto's, although the Prugniel-Simien and Dehnen-McLaughlin models also perform well. With regard to the spherical collapse halos, both the Prugniel-Simien and Einasto models describe the density profiles well, with an rms scatter some four times smaller than that obtained with either the NFW-like model or the 3-parameter Dehnen-McLaughlin model. Finally, we confirm recent claims of a systematic variation in profile shape with halo mass.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Empirical Models for Dark Matter Halos. II. Inner profile slopes, dynamical profiles, and rho/sigma^3

We have recently shown that both the Prugniel-Simien model and Sersic's function (hereafter referred to as the Einasto model when applied to internal density profiles) describe simulated dark matter halos better than an NFW-like model with equal number of parameters. Here we provide analytical expressions for the logarithmic slopes of these models, and compare them with data from real galaxies. Depending on the Einasto parameters of the dark matter halo, one can expect an extrapolated, inner (0.01--1 kpc), logarithmic profile slope ranging from -0.2 to -1.5, with a typical value at 0.1 kpc around -0.7. Application of this (better fitting) model therefore alleviates some of the past disagreement with observations on this issue. We additionally provide useful expressions for the concentration and assorted scale radii: r_s, r_{-2}, r_e, R_e, r_virial, and r_max -- the radius where the circular velocity profile has its maximum value. We also present the circular velocity profiles and the radial behavior of rho(r)/sigma(r)^3 for both the Einasto and Prugniel-Simien models. We find this representation of the phase-space density profile to be well approximated by a power-law with slope slightly shallower than -2 near r=r_{-2}.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Empirical Models for Dark Matter Halos. III. The Kormendy relation and the log(rho_e)-log(R_e) relation

We have recently shown that the 3-parameter density-profile model from Prugniel & Simien provides a better fit to simulated, galaxy- and cluster-sized, dark matter halos than an NFW-like model with arbitrary inner profile slope gamma (Paper I). By construction, the parameters of the Prugniel-Simien model equate to those of the Sersic R^{1/n} function fitted to the projected distribution. Using the Prugniel-Simien model, we are therefore able to show that the location of simulated (10^{12} M_sun) galaxy-sized dark matter halos in the <mu>_e-log(R_e) diagram coincides with that of brightest cluster galaxies, i.e., the dark matter halos appear consistent with the Kormendy relation defined by luminous elliptical galaxies. These objects are also seen to define the new, and equally strong, relation log(rho_e) = 0.5 - 2.5log(R_e), in which rho_e is the internal density at r=R_e. Simulated (10^{14.5} M_sun) cluster-sized dark matter halos and the gas component of real galaxy clusters follow the relation log(rho_e) = 2.5[1 - log(R_e)]. Given the shapes of the various density profiles, we are able to conclude that while dwarf elliptical galaxies and galaxy clusters can have dark matter halos with effective radii of comparable size to the effective radii of their baryonic component, luminous elliptical galaxies can not. For increasingly large elliptical galaxies, with increasingly large profile shapes `n', to be dark matter dominated at large radii requires dark matter halos with increasingly large effective radii compared to the effective radii of their stellar component.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

A universal density slope - velocity anisotropy relation

One can solve the Jeans equation analytically for equilibrated dark matter structures, once given two pieces of input from numerical simulations. These inputs are 1) a connection between phase-space density and radius, and 2) a connection between velocity anisotropy and density slope, the α-βrelation. The first (phase-space density v.s. radius) has been analysed through several different simulations, however the second (α-βrelation) has not been quantified yet. We perform a large set of numerical experiments in order to quantify the slope and zero-point of the α-βrelation. When combined with the assumption of phase-space being a power-law in radius this allows us to conclude that equilibrated dark matter structures indeed have zero central velocity anisotropy, central density slope of α_0 = -0.8, and outer anisotropy of approximately β_\infinity = 0.5.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Fate of Supermassive Black Holes and Evolution of the M(BH)-sigma Relation in Merging Galaxies: The Effect of Gaseous Dissipation

We analyze the effect of dissipation on the orbital evolution of supermassive black holes (SMBHs) using high-resolution self-consistent gasdynamical simulations of binary equal- and unequal-mass mergers of disk galaxies. The galaxy models are consistent with the LCDM paradigm of structure formation and the simulations include the effects of radiative cooling and star formation. We find that equal-mass mergers always lead to the formation of a close SMBH pair at the center of the remnant with separations limited solely by the adopted force resolution of ~ 100 pc. Instead, the final SMBH separation in unequal-mass mergers depends sensitively on how the central structure of the merging galaxies is modified by dissipation. In the absence of dissipation, the satellite galaxy can be entirely disrupted before the merger is completed leaving its SMBH wandering at a distance too far from the center of the remnant for the formation of a close pair. In contrast, we show that gas cooling facilitates the pairing process by increasing the resilience of the companion galaxy to tidal disruption. Moreover, we demonstrate that merging disk galaxies constructed to obey the M(BH)-sigma relation, move relative to it depending on whether they undergo a dissipational or collisionless merger, regardless of the mass ratio of the merging systems. Collisionless simulations reveal that remnants tend to move away from the mean relation highlighting the role of gas-poor mergers as a possible source of scatter. In dissipational mergers, the interplay between strong gas inflows associated with the formation of massive nuclear disks and the consumption of gas by star formation provides the necessary fuel to the SMBHs and allows the merger remnants to satisfy the relation.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Caustic Rings and Cold Dark Matter

The hierarchical cold dark matter (CDM) model for structure formation is a well defined and testable model. Direct detection is the best technique for confirming the model yet predictions for the energy and density distribution of particles on earth remain inadequate. Axially symmetric collapse of collisionless dark matter can leave observable caustic rings in phase space and this model is frequently used to make experimental predictions (Sikivie 1999). Such cold collapses inevitably suffer from radial orbit instabilities that produce unrealistic bar-like halos. Moreover, this model bears no relation to the hierarchical growth of CDM galactic halos which form via a complicated sequence of mergers and violent relaxation. This process destroys any symmetry and phase wraps existing caustics on a scale comparable to the first objects to collapse. Since axions can cluster on microscopic scales and free streaming of neutralinos only erases structure smaller than \sim 100 (GeV/m_{CDM}) A.U., the dynamical effects of caustics in the Galactic halo are expected to be negligible.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The Metamorphosis of Tidally Stirred Dwarf Galaxies

We present results from high-resolution N-Body/SPH simulations of rotationally supported dwarf irregular galaxies moving on bound orbits in the massive dark matter halo of the Milky Way.The dwarf models span a range in disk surface density and the masses and sizes of their dark halos are consistent with the predictions of cold dark matter cosmogonies. We show that the strong tidal field of the Milky Way determines severe mass loss in their halos and disks and induces bar and bending instabilities that transform low surface brightness dwarfs (LSBs) into dwarf spheroidals (dSphs) and high surface brightness dwarfs (HSBs) into dwarf ellipticals (dEs) in less than 10 Gyr. The final central velocity dispersions of the remnants are in the range 8-30 km/s and their final $v/σ$ falls to values $< 0.5$, matching well the kinematics of early-type dwarfs. The transformation requires the orbital time of the dwarf to be $\simlt 3-4$ Gyr, which implies a halo as massive and extended as predicted by hierarchical models of galaxy formation to explain the origin of even the farthest dSph satellites of the Milky Way, Leo I and Leo II. Only dwarfs with central dark matter densities as high as those of Draco and Ursa Minor can survive for 10 Gyr in the proximity of the Milky Way: this is naturally achieved within hierarchical models, where the densest objects should have small orbital times due to their early formation epochs. Part of the gas is stripped and part is funneled to the center due to the bar, generating one strong burst of star formation in HSBs and smaller, multiple bursts in LSBs. Extended low-surface brightness stellar and gaseous streams originate from LSBs and, when projected along the line of sight, can lead to overestimate the mass-to-light ratio of the bound remnant by a factor $\simlt 2$,

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Collisional versus Collisionless Dark Matter

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 $σ_{SIDM}/σ_{CDM}=0.85$ and are more circular than CDM, $β_{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.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Tidal stirring and the origin of dwarf spheroidals in the Local Group

N-Body/SPH simulations are used to study the evolution of dwarf irregular galaxies (dIrrs) entering the dark matter halo of the Milky Way or M31 on plunging orbits. We propose a new dynamical mechanism driving the evolution of gas rich, rotationally supported dIrrs, mostly found at the outskirts of the Local Group (LG), into gas free, pressure supported dwarf spheroidals (dSphs) or dwarf ellipticals (dEs), observed to cluster around the two giant spirals. The initial model galaxies are exponential disks embedded in massive dark matter halos and reproduce nearby dIrrs. Repeated tidal shocks at the pericenter of their orbit partially strip their halo and disk and trigger dynamical instabilities that dramatically reshape their stellar component. After only 2-3 orbits low surface brightness (LSB) dIrrs are transformed into dSphs, while high surface brightness (HSB) dIrrs evolve into dEs. This evolutionary mechanism naturally leads to the morphology-density relation observed for LG dwarfs. Dwarfs surrounded by very dense dark matter halos, like the archetypical dIrr GR8, are turned into Draco or Ursa Minor, the faintest and most dark matter dominated among LG dSphs. If disks include a gaseous component, this is both tidally stripped and consumed in periodic bursts of star formation. The resulting star formation histories are in good qualitative agreement with those derived using HST color-magnitude diagrams for local dSphs.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Dark Matter Substructure in Galactic Halos

We use numerical simulations to examine the substructure within galactic and cluster mass halos that form within a hierarchical universe. Clusters are easily reproduced with a steep mass spectrum of thousands of substructure clumps that closely matches observations. However, the survival of dark matter substructure also occurs on galactic scales, leading to the remarkable result that galaxy halos appear as scaled versions of galaxy clusters. The model predicts that the virialised extent of the Milky Way's halo should contain about 500 satellites with circular velocities larger than Draco and Ursa-Minor i.e. bound masses > 10^8Mo and tidally limited sizes > kpc. The substructure clumps are on orbits that take a large fraction of them through the stellar disk leading to significant resonant and impulsive heating. Their abundance and singular density profiles has important implications for the existence of old thin disks, cold stellar streams, gravitational lensing and indirect/direct detection experiments.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

The Clustering Properties of IRAS Galaxies

We investigate the clustering of galaxies in the QDOT redshift survey of IRAS galaxies. We find that the redshift space two-point correlation function is well approximated by a power-law of slope -1.11+\-0.09, with clustering length 3.87+\-0.32 Mpc/h out to pair separations of 25 Mpc\h. On scales larger than 40 Mpc/h, the correlation function is consistent with zero. The r.m.s. fluctuation in the count of QDOT galaxies above the Poisson level in spheres of radius 8 Mpc/h is σ_8^{IRAS}=0.58+\-0.14, showing that fluctuations in the distribution of IRAS galaxies on these scales are smaller than those of optical galaxies by a factor of about 0.65, with an uncertainty of 25%. We find no detectable difference between the correlation functions measured in redshift space and in real space, leading to a $2σ$ limit of b_{IRAS}/Ω^{0.6} >1.05, where b_{IRAS} is the bias factor for IRAS galaxies and Ωis the cosmological density parameter. The QDOT autocorrelation function calculated in concentric shells increases significantly with shell radius. This difference is more likely due to sampling fluctuations than to an increase of the clustering strength with galaxy luminosity, but the two effects are difficult to disentangle; our data allow at most an increase of \sim 20% in clustering strength for each decade in luminosity. For separations greater than \sim 3 Mpc/h, the cross-correlation function of Abell clusters (with Richness R\ge 1) and QDOT galaxies is well approximated by a power-law of slope -1.81+\-0.10, with clustering length 10.10+\-0.45 Mpc/h, and no significant signal beyond \gsim 50 Mpc\h. This cross-correlation depends only weakly on cluster richness.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

The Origin of the Magellanic Stream

SHORTENED ABSTRACT: We present numerical investigations designed to critically test models of the origin of the Magellanic Stream. The most developed model is the tidal model which fails to reproduce several of its characteristic properties. We suggest an alternative model for the origin of the Stream which can explain all of its observed features and dynamics, as well as provide a strong constraint on the distribution of gas within the halo of the Milky Way. We propose that the Stream consists of material which was ram-pressure stripped from the Magellanic System during its last passage through an extended ionized disk of the Galaxy. This collision took place some 500 million years ago at a galacto-centric distance of about 65 kpc, and swept $\sim 20$\% of the least bound HI into the Stream. The gas with the lowest column density lost the most orbital angular momentum, and is presently at the tip of the Stream, having fallen to a distance of $\sim 20$ kpc from the Milky Way attaining a negative velocity of 200 \kms. To prevent the stripped material from leading the Magellanic Clouds and attaining too large an infall velocity, we postulate the existence of an extended dilute halo of diffuse ionized gas surrounding the Milky Way. If the halo gas is at the virial temperature of the potential well of the Milky Way, its thermal emission would contribute $\sim$ 40\% of the observed diffuse background radiation in the 0.5-1.0 keV (M) band, consistent with recent ROSAT measurements as well as pulsar dispersion measures. Ram pressure stripping

0 citations0 reviewsNo code linkedOpen access