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

H. M. P. Couchman

H. M. P. Couchman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-phastro-ph.SR

Catalog footprint

What is connected

26works
4topics
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

26 published item(s)

preprint2016arXiv

The anatomy of a star-forming galaxy: Pressure-driven regulation of star formation in simulated galaxies

We explore the regulation of star formation in star-forming galaxies through a suite of high-resolution isolated galaxy simulations. We use the SPH code GASOLINE, including photoelectric heating and metal cooling, which produces a multi-phase interstellar medium. We show that representative star formation and feedback sub-grid models naturally lead to a weak, sub-linear dependence between the amount of star formation and changes to star formation parameters. We incorporate these sub-grid models into an equilibrium pressure-driven regulation framework. We show that the sub-linear scaling arises as a consequence of the non-linear relationship between scale height and the effective pressure generated by stellar feedback. Thus, simulated star-formation regulation is sensitive to how well vertical structure in the ISM is resolved. Full galaxy disks experience density waves which drive locally time-dependent star formation. We develop a simple time-dependent, pressure-driven model that reproduces the response extremely well.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Energetics of Cusp Destruction

We present a new analytic estimate for the energy required to create a constant density core within a dark matter halo. Our new estimate, based on more realistic assumptions, leads to a required energy that is orders of magnitude lower than is claimed in earlier work. We define a core size based on the logarithmic slope of the dark matter density profile so that it is insensitive to the functional form used to fit observed data. The energy required to form a core depends sensitively on the radial scale over which dark matter within the cusp is redistributed within the halo. Simulations indicate that within a region of comparable size to the active star forming regions of the central galaxy that inhabits the halo, dark matter particles have their orbits radially increased by a factor of 2--3 during core formation. Thus the inner properties of the dark matter halo, such as halo concentration, and final core size, set the energy requirements. As a result, the energy cost increases slowly with halo mass as M$_{\rm{h}}^{0.3-0.7}$ for core sizes $\lesssim1$ kpc. We use the expected star formation history for a given dark matter halo mass to predict dwarf galaxy core sizes. We find that supernovae alone would create well over 4 kpc cores in $10^{10}$ M$_{\odot}$ dwarf galaxies \emph{if} 100% of the energy were transferred to dark matter particle orbits. We can directly constrain the efficiency factor by studying galaxies with known stellar content and core size, such as Fornax. We find that the efficiency of coupling between stellar feedback and dark matter orbital energy need only be at the 1% level or less to explain Fornax's 1 kpc core.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Properties of Bound and Unbound Molecular Cloud Populations Formed in Galactic Disc Simulations

We explore the effect of galactic environment on properties of molecular clouds. Using clouds formed in a large-scale galactic disc simulation, we measure the observable properties from synthetic column density maps. We confirm that a significant fraction of unbound clouds forms naturally in a galactic disc environment and that a mixed population of bound and unbound clouds can match observed scaling relations and distributions for extragalactic molecular clouds. By dividing the clouds into inner and outer disc populations, we compare their distributions of properties and test whether there are statistically significant differences between them. We find that clouds in the outer disc have lower masses, sizes, and velocity dispersions as compared to those in the inner disc for reasonable choices of the inner/outer boundary. We attribute the differences to the strong impact of galactic shear on the disc stability at large galactocentric radii. In particular, our Toomre analysis of the disc shows a narrowing envelope of unstable masses as a function of radius, resulting in the formation of smaller, lower mass fragments in the outer disc. We also show that the star formation rate is affected by the environment of the parent cloud, and is particularly influenced by the underlying surface density profile of the gas throughout the disc. Our work highlights the strengths of using galaxy-scale simulations to understand the formation and evolution of cloud properties - and the star formation within them - in the context of their environment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Systematic Problems With Using Dark Matter Simulations to Model Stellar Halos

The limits of available computing power have forced models for the structure of stellar halos to adopt one or both of the following simplifying assumptions: (1) stellar mass can be "painted" onto dark matter particles in progenitor satellites; (2) pure dark matter simulations that do not form a luminous galaxy can be used. We estimate the magnitude of the systematic errors introduced by these assumptions using a controlled set of stellar halo models where we independently vary whether we look at star particles or painted dark matter particles, and whether we use a simulation in which a baryonic disk galaxy forms or a matching pure dark matter simulation that does not form a baryonic disk. We find that the "painting" simplification reduces the halo concentration and internal structure, predominantly because painted dark matter particles have different kinematics than star particles even when both are buried deep in the potential well of the satellite. The simplification of using pure dark matter simulations reduces the concentration further, but increases the internal structure, and results in a more prolate stellar halo. These differences can be a factor of 1.5-7 in concentration (as measured by the half-mass radius) and 2-7 in internal density structure. Given this level of systematic uncertainty, one should be wary of overinterpreting differences between observations and the current generation of stellar halo models based on dark matter only simulations when such differences are less than an order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Celestial Buffet: multiple populations and globular cluster formation in dwarf galaxies

We present a framework that explains the commonly observed variation in light element abundances in globular clusters. If globular clusters form in the centres of dwarf galaxies, they will be pumped onto larger orbits as star formation progresses. The potential well will only retain the moderate velocity AGB ejecta, the expected source of enrichment, but not supernova ejecta. There is no need to increase the initial cluster mass, a requirement of self-enrichment scenarios, as all the stars within the dwarf can contribute. As the clusters move through the dwarf centre they sweep up a mix of AGB ejecta and in-falling pristine gas to form a second generation of stars. The specific mix will vary in time and is thus able to explain the spread in second generation abundances observed in different clusters. The globular clusters will survive to the present day or be stripped as part of the hierarchical merging process of larger galaxies. We illustrate how this process may operate using a high-resolution simulation of a dwarf galaxy at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Halo orbits in cosmological disk galaxies: tracers of formation history

We analyze the orbits of stars and dark matter particles in the halo of a disk galaxy formed in a cosmological hydrodynamical simulation. The halo is oblate within the inner ~20 kpc and triaxial beyond this radius. About 43% of orbits are short axis tubes - the rest belong to orbit families that characterize triaxial potentials (boxes, long-axis tubes and chaotic orbits), but their shapes are close to axisymmetric. We find no evidence that the self-consistent distribution function of the nearly oblate inner halo is comprised primarily of axisymmetric short-axis tube orbits. Orbits of all families, and both types of particles are highly eccentric with mean eccentricity >0.6. We find that randomly selected samples of halo stars show no substructure in "integrals of motion" space. However individual accretion events can be clearly identified in plots of metallicity versus formation time. Dynamically young tidal debris is found primarily on a single type of orbit. However, stars associated with older satellites become chaotically mixed during the formation process (possibly due to scattering by the central bulge and disk, and baryonic processes), and appear on all four types of orbits. We find that the tidal debris in cosmological hydrodynamical simulations experiences significantly more chaotic evolution than in collisionless simulations, making it much harder to identify individual progenitors using phase space coordinates alone. However by combining information on stellar ages and chemical abundances with the orbital properties of halo stars in the underlying self-consistent potential, the identification of progenitors is likely to be possible.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

MaGICC Thick Disk I: Comparing a Simulated Disk Formed with Stellar Feedback to the Milky Way

We analyse the structure and chemical enrichment of a Milky Way-like galaxy with a stellar mass of 2 10^{10} M_sun, formed in a cosmological hydrodynamical simulation. It is disk-dominated with a flat rotation curve, and has a disk scale length similar to the Milky Way's, but a velocity dispersion that is ~50% higher. Examining stars in narrow [Fe/H] and [α/Fe] abundance ranges, we find remarkable qualitative agreement between this simulation and observations: a) The old stars lie in a thickened distribution with a short scale length, while the young stars form a thinner disk, with scale lengths decreasing, as [Fe/H] increases. b) Consequently, there is a distinct outward metallicity gradient. c) Mono-abundance populations exist with a continuous distribution of scale heights (from thin to thick). However, the simulated galaxy has a distinct and substantive very thick disk (h_z~1.5 kpc), not seen in the Milky Way. The broad agreement between simulations and observations allows us to test the validity of observational proxies used in the literature: we find in the simulation that mono-abundance populations are good proxies for single age populations (<1 Gyr) for most abundances.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The MaGICC volume: reproducing statistical properties of high redshift galaxies

We present a cosmological hydrodynamical simulation of a representative volume of the Universe, as part of the Making Galaxies in a Cosmological Context (MaGICC) project. MaGICC uses a thermal implementation for supernova and early stellar feedback. This work tests the feedback model at lower resolution across a range of galaxy masses, morphologies and merger histories. The simulated sample compares well with observations of high redshift galaxies ($z \ge 2$) including the stellar mass - halo mass ($M_\star - M_h$ ) relation, the Galaxy Stellar Mass Function (GSMF) at low masses ($M_\star \lt 5 \times 10^{10} M_\odot$ ) and the number density evolution of low mass galaxies. The poor match of $M_\star - M_h$ and the GSMF at high masses ($M_\star \ge 5 \times 10^{10} M_\odot$ ) indicates supernova feedback is insufficient to limit star formation in these haloes. At $z = 0$, our model produces too many stars in massive galaxies and slightly underpredicts the stellar mass around $L_\star$ mass galaxy. Altogether our results suggest that early stellar feedback, in conjunction with supernovae feedback, plays a major role in regulating the properties of low mass galaxies at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Building the Stellar Halo Through Feedback in Dwarf Galaxies

We present a new model for the formation of stellar halos in dwarf galaxies. We demonstrate that the stars and star clusters that form naturally in the inner regions of dwarfs are expected to migrate from the gas rich, star forming centre to join the stellar spheroid. For dwarf galaxies, this process could be the dominant source of halo stars. The effect is caused by stellar feedback-driven bulk motions of dense gas which, by causing potential fluctuations in the inner regions of the halo, couple to all collisionless components. This effect has been demonstrated to generate cores in otherwise cuspy cold dark matter profiles and is particularly effective in dwarf galaxy haloes. It can build a stellar spheroid with larger ages and lower metallicities at greater radii without requiring an outside-in formation model. Globular cluster-type star clusters can be created in the galactic ISM and then migrate to the spheroid on 100\thinspace Myr timescales. Once outside the inner regions they are less susceptible to tidal disruption and are thus long lived; clusters on wider orbits may be easily unbound from the dwarf to join the halo of a larger galaxy during a merger. A simulated dwarf galaxy ($\text{M}_{vir}\simeq10^{9}\text{M}_{\odot}$ at $z=5$) is used to examine this gravitational coupling to dark matter and stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Making Galaxies in a Cosmological Context: The Need for Early Stellar Feedback

We introduce the Making Galaxies in a Cosmological Context (MaGICC) program of smoothed particle hydrodynamics (SPH) simulations. We describe a parameter study of galaxy formation simulations of an L* galaxy that uses early stellar feedback combined with supernova feedback to match the stellar mass--halo mass relationship. While supernova feedback alone can reduce star formation enough to match the stellar mass--halo mass relationship, the galaxy forms too many stars before z=2 to match the evolution seen using abundance matching. Our early stellar feedback is purely thermal and thus operates like a UV ionization source as well as providing some additional pressure from the radiation of massive, young stars. The early feedback heats gas to >10^6 K before cooling to 10^4 K. The pressure from this hot gas creates a more extended disk and prevents more star formation prior to z=1 than supernovae feedback alone. The resulting disk galaxy has a flat rotation curve, an exponential surface brightness profile, and matches a wide range of disk scaling relationships. The disk forms from the inside-out with an increasing exponential scale length as the galaxy evolves. Overall, early stellar feedback helps to simulate galaxies that match observational results at low and high redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Metallicity Gradients in Disks: Do Galaxies Form Inside-Out?

We examine radial and vertical metallicity gradients using a suite of disk galaxy simulations, supplemented with two classic chemical evolution approaches. We determine the rate of change of gradient and reconcile differences between extant models and observations within the `inside-out' disk growth paradigm. A sample of 25 disks is used, consisting of 19 from our RaDES (Ramses Disk Environment Study) sample, realised with the adaptive mesh refinement code RAMSES. Four disks are selected from the MUGS (McMaster Unbiased Galaxy Simulations) sample, generated with the smoothed particle hydrodynamics (SPH) code GASOLINE, alongside disks from Rahimi et al. (GCD+) and Kobayashi & Nakasato (GRAPE-SPH). Two chemical evolution models of inside-out disk growth were employed to contrast the temporal evolution of their radial gradients with those of the simulations. We find that systematic differences exist between the predicted evolution of radial abundance gradients in the RaDES and chemical evolution models, compared with the MUGS sample; specifically, the MUGS simulations are systematically steeper at high-redshift, and present much more rapid evolution in their gradients. We find that the majority of the models predict radial gradients today which are consistent with those observed in late-type disks, but they evolve to this self-similarity in different fashions, despite each adhering to classical `inside-out' growth. We find that radial dependence of the efficiency with which stars form as a function of time drives the differences seen in the gradients; systematic differences in the sub-grid physics between the various codes are responsible for setting these gradients. Recent, albeit limited, data at redshift z=1.5 are consistent with the steeper gradients seen in our SPH sample, suggesting a modest revision of the classical chemical evolution models may be required.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Aquila comparison Project: The Effects of Feedback and Numerical Methods on Simulations of Galaxy Formation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Distribution of Metals in Cosmological Hydrodynamical Simulations of Dwarf Disk Galaxies

We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass -- halo mass relationship. Each realisation employed identical initial conditions and assembly histories, but differed from one another in their underlying baryonic physics prescriptions, including (a) the efficiency with which each supernova energy couples to the ISM, (b) the impact of feedback associated with massive star radiation pressure, (c) the role of the minimum shut-off time for radiative cooling of Type II SNe remnants, (d) the treatment of metal diffusion, and (e) varying the IMF. Our analysis focusses on the resulting stellar metallicity distribution functions (MDFs) in each simulated (analogous) `solar neighbourhood' and central `bulge' region. We compare the simulated MDFs' skewness, kurtosis, and dispersion (inter-quartile, inter-decile, inter-centile, and inter-tenth-percentile regions) with that of the empirical solar neighbourhood MDF and Local Group dwarfs. We find that the MDFs of the simulated discs are more negatively skewed, with higher kurtosis, than those observed locally. We can trace this difference to the simulations' tight and correlated age-metallicity relations (compared with that of the Milky Way), suggesting that these relations within `dwarf' discs might be steeper than in L* discs and/or the degree of stellar orbital re-distribution and migration inferred locally has not been captured in their entirety, at the resolution of our simulations. The important role of metal diffusion in ameliorating the over-production of extremely metal-poor stars is highlighted.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The luminosity function of diverse satellite galaxy systems

The high-resolution, SPH galaxies of the McMaster Unbiased Galaxy Survey (MUGS) are used to examine the satellite systems of sixteen model host galaxies. Each galaxy has a different mass, angular momentum and merger history that yield a rich set of satellite luminosity functions. With new observations of distant satellite systems, we can compare these luminosity functions to satellite systems beyond the Local Group. We find that the luminosity functions of our simulations compare well to observations when the luminosity functions are scaled according to host mass. We use the recently-found relationship between dwarf satellites and host mass in distant satellite systems (Trentham & Tully 2009) to normalize a theoretical, complete luminosity function for the Milky Way (Koposov et al. 2008). The luminosity function of satellites, expressed as a function of the host mass, is given by dN/dM_V= 3.5 x M_host^(0.91) x 10^(0.1 x M_V-10.2), where host mass is given in M_sun. The mass of a host galaxy can be used to predict the number of dwarf satellites and even when considering spiral and elliptical hosts separately this relation holds.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The stellar metallicity distribution of disc galaxies and bulges in cosmological simulations

By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The effects of the star formation density threshold are visible in the star formation histories, which show a modulation in their behaviour driven by the threshold. The derived MDFs show median metallicities lower by 0.2-0.3 dex than the MDF observed locally in the disc and in the Galactic bulge. Possible reasons for this apparent discrepancy include the use of low stellar yields and/or centrally-concentrated star formation. The dispersions are larger than the one of the observed MDF; this could be due to simulated discs being kinematically hotter relative to the Milky Way. The fraction of low metallicity stars is largely overestimated, visible from the more negatively skewed MDF with respect to the observational sample. For our fiducial Milky Way analog, we study the metallicity distribution of the stars born "in situ" relative to those formed via accretion (from disrupted satellites), and demonstrate that this low-metallicity tail to the MDF is populated primarily by accreted stars. Enhanced supernova and stellar radiation energy feedback to the surrounding interstellar media of these pre-disrupted satellites is suggested as an important regulator of the MDF skewness.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Halo expansion in cosmological hydro simulations: towards a baryonic solution of the cusp/core problem in massive spirals

A clear prediction of the Cold Dark Matter model is the existence of cuspy dark matter halo density profiles on all mass scales. This is not in agreement with the observed rotation curves of spiral galaxies, challenging on small scales the otherwise successful CDM paradigm. In this work we employ high resolution cosmological hydro-dynamical simulations to study the effects of dissipative processes on the inner distribution of dark matter in Milky-Way like objects (M~1e12 Msun). Our simulations include supernova feedback, and the effects of the radiation pressure of massive stars before they explode as supernovae. The increased stellar feedback results in the expansion of the dark matter halo instead of contraction with respect to N-body simulations. Baryons are able to erase the dark matter cuspy distribution creating a flat, cored, dark matter density profile in the central several kpc of a massive Milky-Way like halo. The profile is well fit by a Burkert profile, with fitting parameters consistent with the observations. In addition, we obtain flat rotation curves as well as extended, exponential stellar disk profiles. While the stellar disk we obtain is still partially too thick to resemble the MW thin disk, this pilot study shows that there is enough energy available in the baryonic component to alter the dark matter distribution even in massive disc galaxies, providing a possible solution to the long standing problem of cusps vs. cores.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mechanisms of Baryon Loss for Dark Satellites in Cosmological SPH Simulations

We present a study of satellites in orbit around a high-resolution, smoothed particle hydrodynamics (SPH) galaxy simulated in a cosmological context. The simulated galaxy is approximately the same mass as the Milky Way. The cumulative number of luminous satellites at z = 0 is similar to the observed system of satellites orbiting the Milky Way although an analysis of the satellite mass function reveals an order of magnitude more dark satellites than luminous. Some of the dark subhalos are more massive than some of the luminous subhalos at z = 0. What separates luminous and dark subhalos is not their mass at z = 0, but the maximum mass the subhalos ever achieve. We study the effect of four mass loss mechanisms on the subhalos: ultraviolet (UV) ionising radiation, ram pressure stripping, tidal stripping, and stellar feedback, and compare the impact of each of these four mechanisms on the satellites. In the lowest mass subhalos, UV is responsible for the majority of the baryonic mass loss. Ram pressure stripping removes whatever mass remains from the low mass satellites. More massive subhalos have deeper potential wells and retain more mass during reionisation. However, as satellites pass near the centre of the main halo, tidal forces cause significant mass loss from satellites of all masses. Satellites that are tidally stripped from the outside can account for the luminous satellites that are lower mass than some of the dark satellites. Stellar feedback has the greatest impact on medium mass satellites that had formed stars, but lost all their gas by z = 0. Our results demonstrate that the missing satellite problem is not an intractable issue with the cold dark matter cosmology, but is rather a manifestation of baryonic processes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Estimating small angular scale CMB anisotropy with high resolution N-body simulations: weak lensing

We estimate the impact of weak lensing by strongly nonlinear cosmological structures on the cosmic microwave background. Accurate calculation of large $\ell$ multipoles requires N-body simulations and ray-tracing schemes with both high spatial and temporal resolution. To this end we have developed a new code that combines a gravitational Adaptive Particle-Particle, Particle-Mesh (AP3M) solver with a weak lensing evaluation routine. The lensing deviations are evaluated while structure evolves during the simulation so that all evolution steps--rather than just a few outputs--are used in the lensing computations. The new code also includes a ray-tracing procedure that avoids periodicity effects in a universe that is modeled as a 3-D torus in the standard way. Results from our new simulations are compared with previous ones based on Particle-Mesh simulations. We also systematically investigate the impact of box volume, resolution, and ray-tracing directions on the variance of the computed power spectra. We find that a box size of $512 h^{-1}$ Mpc is sufficient to provide a robust estimate of the weak lensing angular power spectrum in the $\ell$-interval (2,000--7,000). For a reaslistic cosmological model the power $[\ell(\ell+1)C_{\ell}/2π]^{1/2}$ takes on values of a few $μK$ in this interval, which suggests that a future detection is feasible and may explain the excess power at high $\ell$ in the BIMA and CBI observations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Angular Momentum Transport in Extended Galactic Disks

We demonstrate a significant difference in the angular momentum transport properties of galactic disks between regions in which the interstellar medium is single phase or two phase. Our study is motivated by observations of HI in extended galactic disks which indicate velocity dispersions of nonthermal origin, suggesting that turbulence in the gas may be contributing significantly to the observed dispersion. To address this, we have implemented a shearing-box framework within the FLASH code. The new code was used to perform local simulations of galactic disks that incorporate differential rotation, self-gravity, vertical stratification, hydrodynamics and cooling. These simulations explore plausible mechanisms for driving turbulent motions via the thermal and self-gravitational instabilities coupling to differential rotation. Where a two-phase medium develops, gravitational angular momentum transporting stresses are much greater, creating a possible mechanism for transferring energy from galactic rotation to turbulence. In simulations where the disk conditions do not trigger the formation of a two-phase medium, it is found that perturbations to the flow damp without leading to a sustained mechanism for driving turbulence. The differing angular momentum transport properties of the single- and two-phase regimes of the disk suggest that a significant, dynamically motivated division can be drawn between the two, even when this division occurs far outside the star formation cutoff in a galactic disk.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Predictions of Quasar Clustering: Redshift, Luminosity and Selection Dependence

We show that current clustering observations of quasars and luminous AGN can be explained by a merger model augmented by feedback from outflows. Using numerical simulations large enough to study clustering out to 25 comoving h^{-1} Mpc, we calculate correlation functions, biases, and correlation lengths as a function of AGN redshift and optical and X-ray luminosity. At optical wavelengths, our results match a wide range of current observations and generate predictions for future data sets. We reproduce the weak luminosity dependence of clustering over the currently well-measured range, and predict a much stronger dependence at higher luminosities. The increase in the amplitude of binary quasar clustering observed in the Sloan Digital Sky Survey (SDSS) is also reproduced and is predicted to occur at higher redshift, an effect that is due to the one halo term in the correlation function. On the other hand, our results do not match the rapid evolution of the correlation length observed in the SDSS at z\simeq 3, a discrepancy that is at least partially due to differences in the scales probed by our simulation versus this survey. In fact, we show that changing the distances sampled from our simulations can produce changes as large as 40% in the fitted correlation lengths. Finally, in the X-ray, our simulations produce correlation lengths similar to that observed in the Chandra Deep Field (CDF) North, but not the significantly larger correlation length observed in the CDF South.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Implementing Feedback in Simulations of Galaxy Formation: A Survey of Methods

We present a detailed investigation of different approaches to modeling feedback in simulations of galaxy formation. Gas-dynamic forces are evaluated using Smoothed Particle Hydrodynamics (SPH) while star formation and supernova feedback are included using a three parameter model which determines the star formation rate normalization, feedback energy and lifetime of feedback regions. The star formation rate is calculated using a Lagrangian Schmidt Law for all gas particles which satisfy temperature, density and convergent flow criteria. Feedback is incorporated as thermal heating of the ISM. We compare the effects of distributing this energy over the smoothing scale or depositing it on a single particle. Radiative losses are prevented from heated particles by adjusting the density used in radiative cooling. We test the models on the formation of galaxies from cosmological initial conditions and also on isolated Milky Way and dwarf galaxies. Extremely violent feedback is necessary to produce a gas disk with angular momentum remotely close to that of observed disk galaxies. This is a result of the extreme central concentration of the dark halos in the sCDM model, and the pervasiveness of the core-halo angular momentum transport mechanism. We emphasize that the disks formed in hierarchical simulations are partially a numerical artifact produced by the minimum mass scale of the simulation acting as a highly efficient `support' mechanism. Disk formation is strongly affected by the treatment of dense regions in SPH, which along with the difficulty of representing the hierarchical formation process, means that realistic simulations of galaxy formation require far higher resolution than currently used.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Evolution of structure in cold dark matter universes

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.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

The structure of galaxy clusters in different cosmologies

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).

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

The Virgo consortium: simulations of dark matter and galaxy clustering

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 $Ω_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 $Ω_0$ N-body/SPH simulation appears to be related to the merging and disruption of galaxies in rich clusters.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Mergers of Systems Containing Gas

Several simple mergers between model galaxy clusters containing a mixture of gas and dark matter are examined, testing the coupling of the gas to the underlying collisionless material. The gas is shocked, irreversibly dissipating the energy fed into it by the collisionless component and forms a resolved constant-density core. For the dark matter, however, admixture of phase space vacuum is not very efficient and a constant-density core is not produced. In the final state the central gas has little residual kinetic energy, indicating that streaming motions do not help to support the gas.

0 citations0 reviewsNo code linkedOpen access