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

James Wadsley

James Wadsley 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.SRastro-ph.EPastro-ph.HEastro-ph.IMhep-phnlin.CD

Catalog footprint

What is connected

70works
9topics
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

70 published item(s)

preprint2019arXiv

LRP2020: Signposts of planet formation in protoplanetary disks

Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Clumpy high-z galaxies as a testbed for feedback-regulated galaxy formation

We study the dependence of fragmentation in massive gas-rich galaxy disks at $z > 1$ on feedback model and hydrodynamical method, employing the GASOLINE2 SPH code and the lagrangian mesh-less code GIZMO in finite mass mode. We compare non-cosmological galaxy disk runs with standard blastwave supernovae (SN)feedback, which introduces delayed cooling in order to drive winds, and runs with the new superbubble SN feedback, which produces winds naturally by modelling the detailed physics of SN-driven bubbles and leads to efficient self-regulation of star formation. We find that, with blastwave feedback, massive star forming clumps form in comparable number and with very similar masses in GASOLINE2 and GIZMO. The typical masses are in the range $10^7-10^8 M_{\odot}$, lower than in most previous works, while giant clumps with masses above $10^9 M_{\odot}$ are exceedingly rare. With superbubble feedback, instead, massive bound star forming clumps do not form because galaxies never undergo a phase of violent disk instability. Only sporadic, unbound star forming overdensities lasting only a few tens of Myr can arise that are triggered by perturbations of massive satellite companions. We conclude that there is a severe tension between explaining massive star forming clumps observed at $z > 1$ primarily as the result of disk fragmentation driven by gravitational instability and the prevailing view of feedback-regulated galaxy formation. The link between disk stability and star formation efficiency should thus be regarded as a key testing ground for galaxy formation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Direct detection of precursors of gas giants formed by gravitational instability with the Atacama Large Millimetre/sub-millimetre Array

Phases of gravitational instability are expected in the early phases of disk evolution, when the disk mass is still a substantial fraction of the mass of the star. Disk fragmentation into sub-stellar objects could occur in the cold exterior part of the disk. Direct detection of massive gaseous clumps on their way to collapse into gas giant planets would offer an unprecedented test of the disk instability model. Here we use state-of-the-art 3D radiation-hydro simulations of disks undergoing fragmentation into massive gas giants, post-processed with the RADMC-3D ray-tracing code to produce dust continuum emission maps. These are then fed into the Common Astronomy Software Applications (CASA) ALMA simulator. The synthetic maps show that both overdense spiral arms and actual clumps at different stages of collapse can be detected with the Atacama Large Millimetre/sub-millimetre Array (ALMA) in the full configuration at the distance of the Ophiuchus star forming region (125 pc). The detection of clumps is particularly effective at shorter wavelengths (690 GHz) combining two resolutions with multi-scale clean. Furthermore, we show that a flux-based estimate of the mass of a protoplanetary clump can be from comparable to a factor of 3 higher than the gravitationally bound clump mass. The estimated mass depends on the assumed opacity, and on the gas temperature, which should be set using the input of radiation-hydro simulations. We conclude that ALMA has the capability to detect "smoking gun" systems that are a signpost of the disk instability model for gas giant planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effects of Galactic Disc Inclination and Resolution on Observed GMC Properties and Larson's Scaling Relations

With ALMA making it possible to resolve giant molecular clouds (GMCs) in other galaxies, it is becoming necessary to quantify the observational bias on measured GMC properties. Using a hydrodynamical simulation of a barred spiral galaxy, we compared the physical properties of GMCs formed in position-position-position space (PPP) to the observational position-position-velocity space (PPV). We assessed the effect of disc inclination: face-on (PPV_face) and edge-on (PPV_edge), and resolution: 1.5 pc versus 24 pc, on GMC properties and the further implications of using Larson's scaling relations for mass-radius and velocity dispersion-radius. The low-resolution PPV data are generated by simulating ALMA Cycle 3 observations using the CASA package. Results show that the median properties do not differ strongly between PPP and PPV_face under both resolutions, but PPV_edge clouds deviate from these two. The differences become magnified when switching to the lower, but more realistic resolution. The discrepancy can lead to opposite results for the virial parameter's measure of gravitational binding, and therefore the dynamical state of the clouds. The power-law indices for the two Larson's scaling relations decrease going from PPP, PPV_face to PPV_edge and decrease from high to low resolutions. We conclude that the relations are not entirely driven by the underlying physical origin and therefore have to be used with caution when considering the environmental dependence, dynamical state, and the extragalactic CO-to-H2 conversion factor of GMCs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

NIHAO project II: Halo shape, phase-space density and velocity distribution of dark matter in galaxy formation simulations

We use the NIHAO (Numerical Investigation of Hundred Astrophysical Objects) cosmological simulations to study the effects of galaxy formation on key properties of dark matter (DM) haloes. NIHAO consists of $\simeq 90$ high-resolution SPH simulations that include (metal-line) cooling, star formation, and feedback from massive stars and SuperNovae, and cover a wide stellar and halo mass range: $10^6 < M_* / M_{\odot} < 10^{11}$ ( $10^{9.5} < M_{\rm halo} / M_{\odot} < 10^{12.5}$). When compared to DM-only simulations, the NIHAO haloes have similar shapes at the virial radius, R_{\rm vir}, but are substantially rounder inside $\simeq 0.1R_{\rm vir}$. In NIHAO simulations $c/a$ increases with halo mass and integrated star formation efficiency, reaching $\sim 0.8$ at the Milky Way mass (compared to 0.5 in DM-only), providing a plausible solution to the long-standing conflict between observations and DM-only simulations. The radial profile of the phase-space $Q$ parameter ($ρ/σ^3$) is best fit with a single power law in DM-only simulations, but shows a flattening within $\simeq 0.1R_{\rm vir}$ for NIHAO for total masses $M>10^{11} M_{\odot}$. Finally, the global velocity distribution of DM is similar in both DM-only and NIHAO simulations, but in the solar neighborhood, NIHAO galaxies deviate substantially from Maxwellian. The distribution is more symmetric, roughly Gaussian, with a peak that shifts to higher velocities for Milky Way mass haloes. We provide the distribution parameters which can be used for predictions for direct DM detection experiments. Our results underline the ability of the galaxy formation processes to modify the properties of dark matter haloes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Supermassive black hole pairs in clumpy galaxies at high redshift: delayed binary formation and concurrent mass growth

Massive gas-rich galaxy discs at $z \sim 1-3$ host massive star-forming clumps with typical baryonic masses in the range $10^7-10^8$~M$_{\odot}$ which can affect the orbital decay and concurrent growth of supermassive black hole (BH) pairs. Using a set of high-resolution simulations of isolated clumpy galaxies hosting a pair of unequal-mass BHs, we study the interaction between massive clumps and a BH pair at kpc scales, during the early phase of the orbital decay. We find that both the interaction with massive clumps and the heating of the cold gas layer of the disc by BH feedback tend to delay significantly the orbital decay of the secondary, which in many cases is ejected and then hovers for a whole Gyr around a separation of 1--2 kpc. In the envelope, dynamical friction is weak and there is no contribution of disc torques: these lead to the fastest decay once the orbit of the secondary BH has circularised in the disc midplane. In runs with larger eccentricities the delay is stronger, although there are some exceptions. We also show that, even in discs with very sporadic transient clump formation, a strong spiral pattern affects the decay time-scale for BHs on eccentric orbits. We conclude that, contrary to previous belief, a gas-rich background is not necessarily conducive to a fast BH decay and binary formation, which prompts more extensive investigations aimed at calibrating event-rate forecasts for ongoing and future gravitational-wave searches, such as with Pulsar Timing Arrays and the future evolved Laser Interferometer Space Antenna.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Uncovering Mass Segregation with Galaxy Analogues in Dark Matter Simulations

We investigate mass segregation in group and cluster environments by identifying galaxy analogues in high-resolution dark matter simulations. Subhalos identified by the AHF and ROCKSTAR halo finders have similar mass functions, independent of resolution, but different radial distributions due to significantly different subhalo hierarchies. We propose a simple way to classify subhalos as galaxy analogues. The radial distributions of galaxy analogues agree well at large halo-centric radii for both AHF and ROCKSTAR but disagree near parent halo centres where the phase-space information used by ROCKSTAR is essential. We see clear mass segregation at small radii (within $0.5\,r_{vir}$) with average galaxy analogue mass decreasing with radius. Beyond the virial radius, we find a mild trend where the average galaxy analogue mass increases with radius. These mass segregation trends are strongest in small groups and dominated by the segregation of low mass analogues. The lack of mass segregation in massive galaxy analogues suggests that the observed trends are driven by the complex accretion histories of the parent halos rather than dynamical friction.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A lower fragmentation mass scale in high redshift galaxies and its implications on giant clumps: a systematic numerical study

We study the effect of sub-grid physics, galaxy mass, structural parameters and resolution on the fragmentation of gas-rich galaxy discs into massive star forming clumps. The initial conditions are set up with the aid of the ARGO cosmological hydrodynamical simulation. Blast-wave feedback does not suppress fragmentation, but reduces both the number of clumps and the duration of the unstable phase. Once formed, bound clumps cannot be destroyed by our feedback model. Widespread fragmentation is promoted by high gas fractions and low halo concentrations. Yet giant clumps $M > 10^8 M_{\odot}$ lasting several hundred Myr are rare and mainly produced by clump-clump mergers. They occur in massive discs with maximum rotational velocities $V_{max} > 250$ km/s at $z \sim 2$, at the high mass end of the observed galaxy population at those redshifts. The typical gaseous and stellar masses of clumps in all runs are in the range $\sim 10^7-10^8 M_{\odot}$ for galaxies with disc mass in the range $10^{10}-8\times 10^{10} M_{\odot}$. Clumps sizes are usually in the range $100-400$ pc, in agreement with recent clump observations in lensed high-z galaxies. \\ We argue that many of the giant clumps identified in observations are not due to in-situ fragmetation, or are the result of blending of smaller structures owing to insufficient resolution. Using an analytical model describing local collapse inside spiral arms, we can predict the characteristic gaseous masses of clumps at the onset of fragmentation ($\sim 3-5 \times 10^7 M_{\odot}$) quite accurately, while the conventional Toomre mass overestimates them. Due to their moderate masses, clumps which migrate to the centre have marginal effect on bulge growth.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

NIHAO project I: Reproducing the inefficiency of galaxy formation across cosmic time with a large sample of cosmological hydrodynamical simulations

We introduce project NIHAO (Numerical Investigation of a Hundred Astrophysical Objects), a set of 100 cosmological zoom-in hydrodynamical simulations performed using the gasoline code, with an improved implementation of the SPH algorithm. The haloes in our study range from dwarf to Milky Way masses, and represent an unbiased sampling of merger histories, concentrations and spin parameters. The particle masses and force softenings are chosen to resolve the mass profile to below 1% of the virial radius at all masses, ensuring that galaxy half-light radii are well resolved. Using the same treatment of star formation and stellar feedback for every object, the simulated galaxies reproduce the observed inefficiency of galaxy formation across cosmic time as expressed through the stellar mass vs halo mass relation, and the star formation rate vs stellar mass relation. We thus conclude that stellar feedback is the chief piece of physics required to limit the efficiency of star formation in galaxies less massive than the Milky Way.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Star Formation in Disk Galaxies. III. Does stellar feedback result in cloud death?

Stellar feedback, star formation and gravitational interactions are major controlling forces in the evolution of Giant Molecular Clouds (GMCs). To explore their relative roles, we examine the properties and evolution of GMCs forming in an isolated galactic disk simulation that includes both localised thermal feedback and photoelectric heating. The results are compared with the three previous simulations in this series which consists of a model with no star formation, star formation but no form of feedback and star formation with photoelectric heating in a set with steadily increasing physical effects. We find that the addition of localised thermal feedback greatly suppresses star formation but does not destroy the surrounding GMC, giving cloud properties closely resembling the run in which no stellar physics is included. The outflows from the feedback reduce the mass of the cloud but do not destroy it, allowing the cloud to survive its stellar children. This suggests that weak thermal feedback such as the lower bound expected for supernova may play a relatively minor role in the galactic structure of quiescent Milky Way-type galaxies, compared to gravitational interactions and disk shear.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The distribution of mass components in simulated disc galaxies

Using 22 hydrodynamical simulated galaxies in a LCDM cosmological context we recover not only the observed baryonic Tully-Fisher relation, but also the observed "mass discrepancy--acceleration" relation, which reflects the distribution of the main components of the galaxies throughout their disks. This implies that the simulations, which span the range 52 < V$_{\rm flat}$ < 222 km/s where V$_{\rm flat}$ is the circular velocity at the flat part of the rotation curve, and match galaxy scaling relations, are able to recover the observed relations between the distributions of stars, gas and dark matter over the radial range for which we have observational rotation curve data. Furthermore, we explicitly match the observed baryonic to halo mass relation for the first time with simulated galaxies. We discuss our results in the context of the baryon cycle that is inherent in these simulations, and with regards to the effect of baryonic processes on the distribution of dark matter.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The effect of stellar feedback on a Milky Way-like galaxy and its gaseous halo

We present the study of a set of N-body+SPH simulations of a Milky Way-like system produced by the radiative cooling of hot gas embedded in a dark matter halo. The galaxy and its gaseous halo evolve for 10 Gyr in isolation, which allows us to study how internal processes affect the evolution of the system. We show how the morphology, the kinematics and the evolution of the galaxy are affected by the input supernova feedback energy E$_{\rm SN}$, and we compare its properties with those of the Milky Way. Different values of E$_{\rm SN}$ do not significantly affect the star formation history of the system, but the disc of cold gas gets thicker and more turbulent as feedback increases. Our main result is that, for the highest value of E$_{\rm SN}$ considered, the galaxy shows a prominent layer of extra-planar cold (log(T)<4.3) gas extended up to a few kpc above the disc at column densities of $10^{19}$ cm$^{-2}$. The kinematics of this material is in agreement with that inferred for the HI halos of our Galaxy and NGC 891, although its mass is lower. Also, the location, the kinematics and the typical column densities of the hot (5.3<log(T)<5.7) gas are in good agreement with those determined from the O$_{\rm VI}$ absorption systems in the halo of the Milky Way and external galaxies. In contrast with the observations, however, gas at log(T)<5.3 is lacking in the circumgalactic region of our systems.

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
preprint2015arXiv

What is a GMC? Are Observers and Simulators Discussing the Same Star-forming Clouds?

As both simulations and observations reach the resolution of the star-forming molecular clouds, it becomes important to clarify if these two techniques are discussing the same objects in galaxies. We compare clouds formed in a high resolution galaxy simulation identified as continuous structures within a contour, in the simulator's position-position-position (PPP) co-ordinate space and the observer's position-position-velocity space (PPV). Results indicate that the properties of the cloud populations are similar in both methods and up to 70% of clouds have a single counterpart in the opposite data structure. Comparing individual clouds in a one-to-one match reveals a scatter in properties mostly within a factor of two. However, the small variations in mass, radius and velocity dispersion produce significant differences in derived quantities such as the virial parameter. This makes it difficult to determine if a structure is truely gravitationally bound. The three cloud types originally found in the simulation in Fujimoto et al. (2014) are identified in both data sets, with around 80% of the clouds retaining their type between identification methods. We also compared our results when using a peak decomposition method to identify clouds in both PPP and PPV space. The number of clouds increased with this technique, but the overall cloud properties remained similar. However, the more crowded environment lowered the ability to match clouds between techniques to 40%. The three cloud types also became harder to separate, especially in the PPV data set. The method used for cloud identification therefore plays a critical role in determining cloud properties, but both PPP and PPV can potentially identify the same structures.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dark MaGICC: the effect of Dark Energy on galaxy formation. Cosmology does matter

We present the Dark MaGICC project, which aims to investigate the effect of Dark Energy (DE) modeling on galaxy formation via hydrodynamical cosmological simulations. Dark MaGICC includes four dynamical Dark Energy scenarios with time varying equations of state, one with a self-interacting Ratra-Peebles model. In each scenario we simulate three galaxies with high resolution using smoothed particle hydrodynamics (SPH). The baryonic physics model is the same used in the Making Galaxies in a Cosmological Context (MaGICC) project, and we varied only the background cosmology. We find that the Dark Energy parameterization has a surprisingly important impact on galaxy evolution and on structural properties of galaxies at z=0, in striking contrast with predictions from pure Nbody simulations. The different background evolutions can (depending on the behavior of the DE equation of state) either enhance or quench star formation with respect to a LCDM model, at a level similar to the variation of the stellar feedback parameterization, with strong effects on the final galaxy rotation curves. While overall stellar feedback is still the driving force in shaping galaxies, we show that the effect of the Dark Energy parameterization plays a larger role than previously thought, especially at lower redshifts. For this reason, the influence of Dark Energy parametrization on galaxy formation must be taken into account, especially in the era of precision cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Direct formation of supermassive black holes in metal-enriched gas at the heart of high-redshift galaxy mergers

We present novel 3D multi-scale SPH simulations of gas-rich galaxy mergers between the most massive galaxies at $z \sim 8 - 10$, designed to scrutinize the direct collapse formation scenario for massive black hole seeds proposed in \citet{mayer+10}. The simulations achieve a resolution of 0.1 pc, and include both metallicity-dependent optically-thin cooling and a model for thermal balance at high optical depth. We consider different formulations of the SPH hydrodynamical equations, including thermal and metal diffusion. When the two merging galaxy cores collide, gas infall produces a compact, optically thick nuclear disk with densities exceeding $10^{-10}$ g cm$^3$. The disk rapidly accretes higher angular momentum gas from its surroundings reaching $\sim 5$ pc and a mass of $\gtrsim 10^9$ $M_{\odot}$ in only a few $10^4$ yr. Outside $\gtrsim 2$ pc it fragments into massive clumps. Instead, supersonic turbulence prevents fragmentation in the inner parsec region, which remains warm ($\sim 3000-6000$ K) and develops strong non-axisymmetric modes that cause prominent radial gas inflows ($> 10^4$ $M_{\odot}$ yr$^{-1}$), forming an ultra-dense massive disky core. Angular momentum transport by non-axisymmetric modes should continue below our spatial resolution limit, quickly turning the disky core into a supermassive protostar which can collapse directly into a massive black hole of mass $10^8-10^9$ $M_{\odot}$ via the relativistic radial instability. Such a "cold direct collapse"' explains naturally the early emergence of high-z QSOs. Its telltale signature would be a burst of gravitational waves in the frequency range $10^{-4} - 10^{-1}$ Hz, possibly detectable by the planned eLISA interferometer.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evolving Molecular Cloud Structure and the Column Density Probability Distribution Function

The structure of molecular clouds can be characterized with the probability distribution function (PDF) of the mass surface density. In particular, the properties of the distribution can reveal the nature of the turbulence and star formation present inside the molecular cloud. In this paper, we explore how these structural characteristics evolve with time and also how they relate to various cloud properties as measured from a sample of synthetic column density maps of molecular clouds. We find that, as a cloud evolves, the peak of its column density PDF will shift to surface densities below the observational threshold for detection, resulting in an underlying lognormal distribution which has been effectively lost at late times. Our results explain why certain observations of actively star-forming, dynamically older clouds, such as the Orion molecular cloud, do not appear to have any evidence of a lognormal distribution in their column density PDFs. We also study the evolution of the slope and deviation point of the power-law tails for our sample of simulated clouds and show that both properties trend towards constant values, thus linking the column density structure of the molecular cloud to the surface density threshold for star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Faint dwarfs as a test of DM models: WDM vs. CDM

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $Λ$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

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
preprint2014arXiv

The effect of models of the interstellar media on the central mass distribution of galaxies

We compare the central mass distribution of galaxies simulated with three different models of the interstellar medium (ISM) with increasing complexity: primordial (H+He) cooling down to 10^4K, additional cooling via metal lines and to lower temperatures, and molecular hydrogen (H_2) with shielding of atomic and molecular hydrogen, in addition to metal line cooling. In order to analyze the effect of these models, we follow the evolution of four field galaxies with V_peak < 120 km/s to a redshift of zero using high-resolution Smoothed Particle Hydrodynamic simulations in a fully cosmological LCDM context. The spiral galaxies produced in simulations with either primordial cooling or H_2 physics have realistic, rising rotation curves. In contrast, the simulations with metal line cooling and otherwise similar feedback and star formation produced galaxies with the peaked rotation curves typical of most previous LCDM simulations of spiral galaxies. The less-massive bulges and non-peaked rotation curves in the galaxies simulated with primordial cooling or H_2 are linked to changes in the angular momentum distribution of the baryons. These galaxies had smaller amounts of low-angular momentum baryons because of increased gas loss from stellar feedback. When there is only primordial cooling, the star forming gas is hotter and the feedback-heated gas cools more slowly than when metal line cooling is included and so requires less energy to be expelled. When H_2 is included, the accompanying shielding produces large amounts of clumpy, cold gas where H_2 forms. Star formation in clumpy gas results in more concentrated supernova feedback and greater efficiency of mass loss. The higher feedback efficiency causes a decrease of low-angular momentum material. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Role of Cold Flows and Reservoirs in Galaxy Formation With Strong Feedback

We examine gas accretion and subsequent star formation in representative galaxies from the McMaster Unbiased Galaxy Simulations (Stinson et al. 2010). Accreted gas is bimodal with a natural temperature division at $10^5$ K, near the peak of the cooling curve. Cold-mode accretion dominates inflows at early times, creating a peak in total accretion at redshift z=2-4 and declining exponentially below z$\sim$2. Hot-mode accretion peaks near z=1-2 and declines gradually. Hot-mode exceeds cold-mode accretion at z$\sim$1.8 for all four galaxies rather than when the galaxy reaches a characteristic mass. Cold-mode accretion can fuel immediate star formation, while hot-mode accretion preferentially builds a large, hot gas reservoir in the halo. Late-time star formation relies on reservoir gas accreted 2-8 Gyr prior. Thus, the reservoir allows the star formation rate to surpass the current overall gas accretion rate. Stellar feedback cycles gas from the interstellar medium back into the hot reservoir. Stronger feedback results in more gas cycling, gas removal in a galactic outflow and less star formation overall, enabling simulations to match the observed star formation history. For lower mass galaxies in particular, strong feedback can delay the star formation peak to z=1-2 from the accretion peak at z=2-4.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Unbound Star-forming Molecular Clouds

We explore whether observed molecular clouds could include a substantial population of unbound clouds. Using simulations which include only turbulence and gravity, we are able to match observed relations and naturally reproduce the observed scatter in the cloud size-linewidth coefficient, at fixed surface density. We identify the source of this scatter as a spread in the intrinsic virial parameter. Thus these observational trends do not require that clouds exist in a state of dynamical equilibrium. We demonstrate that cloud virial parameters can be accurately determined observationally with an appropriate size estimator. All our simulated clouds eventually form collapsing cores, regardless of whether the cloud is bound overall. This supports the idea that molecular clouds do not have to be bound to form stars or to have observed properties like those of nearby low-mass clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Parameter-Space Study of Carbon-Oxygen White Dwarf Mergers

The merger of two carbon-oxygen white dwarfs can lead either to a spectacular transient, stable nuclear burning or a massive, rapidly rotating white dwarf. Simulations of mergers have shown that the outcome strongly depends on whether the white dwarfs are similar or dissimilar in mass. In the similar-mass case, both white dwarfs merge fully and the remnant is hot throughout, while in the dissimilar case, the more massive, denser white dwarf remains cold and essentially intact, with the disrupted lower mass one wrapped around it in a hot envelope and disk. In order to determine what constitutes "similar in mass" and more generally how the properties of the merger remnant depend on the input masses, we simulated unsynchronized carbon-oxygen white dwarf mergers for a large range of masses using smoothed-particle hydrodynamics. We find that the structure of the merger remnant varies smoothly as a function of the ratio of the central densities of the two white dwarfs. A density ratio of 0.6 approximately separates similar and dissimilar mass mergers. Confirming previous work, we find that the temperatures of most merger remnants are not high enough to immediately ignite carbon fusion. During subsequent viscous evolution, however, the interior will likely be compressed and heated as the disk accretes and the remnant spins down. We find from simple estimates that this evolution can lead to ignition for many remnants. For similar-mass mergers, this would likely occur under sufficiently degenerate conditions that a thermonuclear runaway would ensue.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Angular Momentum Acquisition in Galaxy Halos

We use high-resolution cosmological hydrodynamic simulations to study the angular momentum acquisition of gaseous halos around Milky Way sized galaxies. We find that cold mode accreted gas enters a galaxy halo with ~70% more specific angular momentum than dark matter averaged over cosmic time (though with a very large dispersion). In fact, we find that all matter has a higher spin parameter when measured at accretion than when averaged over the entire halo lifetime, and is well characterized by λ~0.1, at accretion. Combined with the fact that cold flow gas spends a relatively short time (1-2 dynamical times) in the halo before sinking to the center, this naturally explains why cold flow halo gas has a specific angular momentum much higher than that of the halo and often forms "cold flow disks". We demonstrate that the higher angular momentum of cold flow gas is related to the fact that it tends to be accreted along filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Giant Molecular Cloud Formation in Disk Galaxies: Characterizing Simulated versus Observed Cloud Catalogues

We present the results of a study of simulated Giant Molecular Clouds (GMCs) formed in a Milky Way-type galactic disk with a flat rotation curve. This simulation, which does not include star formation or feedback, produces clouds with masses ranging between 10^4 Msun and 10^7 Msun. We compare our simulated cloud population to two observational surveys; The Boston University- Five College Radio Astronomy Observatory Galactic Ring Survey and the BIMA All-Disk Survey of M33. An analysis of the global cloud properties as well as a comparison of Larson's scaling relations is carried out. We find that simulated cloud properties agree well with the observed cloud properties, with the closest agreement occurring between the clouds at comparable resolution in M33. Our clouds are highly filamentary - a property that derives both from their formation due to gravitational instability in the sheared galactic environment, as well as to cloud- cloud gravitational encounters. We also find that the rate at which potentially star forming gas accumulates within dense regions - wherein n_{thresh} > 10^4 cm^{-3} - is 3% per 10 Myr, in clouds of roughly 10^6 Msun. This suggests that star formation rates in observed clouds are related to the rates at which gas can be accumulated into dense subregions within GMCs via filamentary flows. The most internally well-resolved clouds are chosen for listing in a catalogue of simulated GMCs; the first of its kind. The catalogued clouds are available as an extracted data set from the global simulation.

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-WDM: the effects of warm dark matter in hydrodynamical simulations of disc galaxy formation

We study the effect of warm dark matter (WDM) on hydrodynamic simulations of galaxy formation as part of the Making Galaxies in a Cosmological Context (MaGICC) project. We simulate three different galaxies using three WDM candidates of 1, 2 and 5 keV and compare results with pure cold dark matter simulations. WDM slightly reduces star formation and produces less centrally concentrated stellar profiles. These effects are most evident for the 1 keV candidate but almost disappear for $m_{\mathrm{WDM}}>2$ keV. All simulations form similar stellar discs independent of WDM particle mass. In particular, the disc scale length does not change when WDM is considered. The reduced amount of star formation in the case of 1 keV particles is due to the effects of WDM on merging satellites which are on average less concentrated and less gas rich. The altered satellites cause a reduced starburst during mergers because they trigger weaker disc instabilities in the main galaxy. Nevertheless we show that disc galaxy evolution is much more sensitive to stellar feedback than it is to WDM candidate mass. Overall we find that WDM, especially when restricted to current observational constraints ($m_{\mathrm{WDM}}>2$ keV), has a minor impact on disc galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Overview of the Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) Project

MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The circumgalactic medium of massive galaxies at z~3: a test for stellar feedback, galactic outflows, and cold streams

We present new results on the kinematics, thermal and ionization state, and spatial distribution of metal-enriched gas in the circumgalactic medium (CGM) of massive galaxies at redshift 3, using the "Eris" suite of cosmological "zoom-in" simulations. The reference run adopts a blastwave scheme for supernova feedback that produces galactic outflows, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, and a model for the diffusion of metals and thermal energy. Synthetic spectra through the multiphase CGM produce interstellar absorption line strengths of Lya, CII, CIV, SiII, and SiIV as a function of galactocentric impact parameter (scaled to the virial radius) that are in broad agreement with those observed at high-redshift by Steidel et al. (2010). Only about one third of all the gas within R_vir is outflowing. The fraction of sightlines within one virial radius that intercept optically thick material is 27%, in agreement with recent observations by Rudie et al. (2012). Such optically thick absorption is shown to trace inflowing "cold" streams that penetrate deep inside the virial radius. The streams, enriched to metallicities above 0.01 solar, give origin to strong (log N > 13) CII absorption with a covering factor of 22% within R_vir and 10% within 2 R_vir. Galactic outflows do not cause any substantial suppression of the cold accretion mode. The central galaxy is surrounded by a large OVI halo, with a typical column density log N>14 and a near unity covering factor maintained all the way out to 150 kpc. This matches the trends recently observed in star-forming galaxies at low redshift by Tumlinson et al. (2011). Our zoom-in simulations of this single system appear to reproduce quantitatively the complex baryonic processes that determine the exchange of matter, energy, and metals between galaxies and their surroundings. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The dependence of dark matter profiles on the stellar to halo mass ratio: a prediction for cusps vs cores

We use 31 simulated galaxies from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter (DM) haloes. The sample covers a wide mass range: 9.4e9<Mhalo/Msun<7.8e11, hosting galaxies with stellar masses: 5.0e5<M*/Msun<8.3e10, i.e. from dwarf to L*. The galaxies are simulated with several baryonic prescriptions, including a range of stellar feedbacks. The main result is a clear dependence of the inner slope of the DM density profile, α in ρr^α, on the ratio between stellar-to-halo mass (M*/Mhalo). This relation is independent of the stellar feedback scheme, allowing a prediction for cusp vs core formation. When M*/Mhalo is low, ~0.01%, energy from stellar feedback is insufficient to significantly alter the inner DM density and the galaxy retains a cuspy profile. At higher M*/Mhalo, feedback drives the expansion of the DM and generates cored profiles. The flattest profiles form where M*/Mhalo~0.5%. Above this ratio, stars formed in the central regions deepen the gravitational potential enough to oppose this supernova-driven expansion process, resulting in smaller cores and cuspier profiles. Combining the dependence of α on M*/Mhalo with the abundance matching relation between M* and Mhalo provides a prediction for how α varies with M*. Further, using the Tully-Fisher relation allows a prediction for the dependence of the DM inner slope on the observed rotation velocity of galaxies. The most cored galaxies are expected to have Vrot~50km/s, with α decreasing for more massive disc galaxies: spirals with Vrot~150km/s have central slopes α<-0.8, approaching the NFW profile. This novel prediction for the dependence of α on disc galaxy mass can be tested using current observational data sets, and can be applied to theoretical modeling of mass profiles and populations of disc galaxies.

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
preprint2013arXiv

The Pressure of the Star Forming ISM in Cosmological Simulations

We examine the pressure of the star-forming interstellar medium (ISM) of Milky-Way sized disk galaxies using fully cosmological SPH+N-body, high resolution simulations. These simulations include explicit treatment of metal-line cooling in addition to dust and self-shielding, $\mathrm{H_{2}}$ based star formation. The 4 simulated halos have masses ranging from a few times $10^{10}$ to nearly $10^{12}$ solar masses. Using a kinematic decomposition of these galaxies into present-day bulge and disk components, we find that the typical pressure of the star-forming ISM in the present-day bulge is higher than that in the present-day disk by an order of magnitude. We also find that pressure of the star-forming ISM at high redshift is on average, higher than ISM pressures at low redshift. This explains the why the bulge forms at higher pressures: the disk assembles at lower redshift, when the ISM is lower pressure and the bulge forms at high redshift, when the ISM is at higher pressure. If ISM pressure and IMF variation are tied together as suggested in studies like \cite{Conroy2012}, these results could indicate a time-dependent IMF in Milky-Way like systems, as well as a different IMF in the bulge and the disk.

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

Baryons Matter: Why Luminous Satellite Galaxies Have Reduced Central Masses

Using high resolution cosmological hydrodynamical simulations of Milky Way-massed disk galaxies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright (-15 < M_V < -8) satellite galaxies. These simulations resolve high density regions, comparable to giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H_2 formation and destruction that ties star formation to the presence of shielded, molecular gas. Before infall, supernova feedback from the clumpy, bursty star formation captured by this physically motivated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the massive satellite progenitors (Mvir > 10^9 Msun, Mstar > 10^7 Msun) compared to DM-only simulations. The progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to both heating at reionization and gas loss from initial star forming events, preserving the steep inner density profile predicted by DM-only simulations. After infall, tidal stripping acts to further reduce the central densities of the luminous satellites, particularly those that enter with cored dark matter halos, increasing the discrepancy in the central masses predicted by baryon+DM and DM-only simulations. We show that DM-only simulations, which neglect the baryonic effects described in this work, produce denser satellites with larger central velocities. We provide a simple correction to the central DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations should be used with great caution when interpreting kinematic observations of the Milky Way's dwarf satellites.

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

Connecting the Dots: Analyzing Synthetic Observations of Star-Forming Clumps in Molecular Clouds

In this paper, we investigate the extent to which observations of molecular clouds can correctly identify and measure star-forming clumps. We produced a synthetic column density map and a synthetic spectral-line data cube from the simulated collapse of a 5000 M$_{\odot}$ molecular cloud. By correlating the clumps found in the simulation to those found in the synthetic observations, clump masses derived from spectral-line data cubes were found to be quite close to the true physical properties of the clumps. We also find that the `observed' clump mass function derived from the column density map is shifted by a factor of ~ 3 higher than the true clump mass function, due to projection of low-density material along the line of sight. Alves et al. (2007) first proposed that a shift of a clump mass function to higher masses by a factor of 3 can be attributed to a star formation efficiency of 30 %. Our results indicate that this finding may instead be due to an overestimate of clump masses determined from column density observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Implementing Molecular Hydrogen in Hydrodynamic Simulations of Galaxy Formation

Motivated by the observed connection between molecular hydrogen (H2) and star formation, we present a method for tracking the non-equilibrium abundance and cooling processes of H2 and H2-based star formation in Smoothed Particle Hydrodynamic simulations. The local abundances of H2 are calculated by integrating over the hydrogen chemical network. This calculation includes the gas-phase and dust grain formation of H2, shielding of HI and H2, and photodissociation of H2 by Lyman-Werner radiation from nearby stellar populations. Because this model does not assume equilibrium abundances, it is particularly well suited for simulations that model low-metallicity environments, such as dwarf galaxies and the early Universe. We further introduce an explicit link between star formation and local H2 abundance. This link limits star formation to "star-forming regions," represented by areas with abundant H2. With this implementation, we determine the effect of H2 on star formation in a cosmological simulation of a dwarf galaxy. This simulation is the first cosmological simulation with non-equilibrium H2 abundances to be integrated to a redshift of zero or to include efficient SN feedback. We find that our simulations are consistent with the observed Tully-Fisher, global Kennicutt-Schmidt, and resolved Kennicutt-Schmidt relations. We find that the inclusion of shielding of both the atomic and molecular hydrogen and, to a lesser extent, the additional cooling from H2 at temperatures between 200 and 5000 K increases the amount of cold gas in the galaxies. The changes to the ISM result in an increased amount of cold, dense gas in the disk of the galaxy and the formation of a clumpier interstellar media (ISM). The explicit link between star formation and H2 and the clumpier ISM results in a bluer galaxy with a greater spatial distribution of star formation at a redshift of zero. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

MAGICC haloes: confronting simulations with observations of the circumgalactic medium at z=0

We explore the circumgalactic medium (CGM) of two simulated star-forming galaxies with luminosities L ~ 0.1 and 1 L* generated using the smooth particle hydrodynamic code GASOLINE. These simulations are part of the Making Galaxies In a Cosmological Context (MAGICC) program in which the stellar feedback is tuned to match the stellar mass-halo mass relationship. For comparison, each galaxy was also simulated using a 'lower feedback' (LF) model which has strength comparable to other implementations in the literature. The 'MAGICC feedback' (MF) model has a higher incidence of massive stars and an approximately two times higher energy input per supernova. Apart from the low-mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends to approximately the virial radius. A significant fraction of this gas has been heated in supernova explosions in the disc and subsequently ejected into the CGM where it is predicted to give rise to substantial O VI absorption. The simulations do not yet address the question of what happens to the O VI when the galaxies stop forming stars. Our models also predict a reservoir of cool H I clouds that show strong Lyαabsorption to several hundred kpc. Comparing these models to recent surveys with the Hubble Space Telescope, we find that only the MF models have sufficient O VI and H I gas in the CGM to reproduce the observed distributions. In separate analyses, these same MF models also show better agreement with other galaxy observables (e.g. rotation curves, surface brightness profiles and H I gas distribution). We infer that the CGM is the dominant reservoir of baryons for galaxy haloes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The origin of metals in the circum-galactic medium of massive galaxies at z=3

We present a detailed study of the metal-enriched circum-galactic medium of a massive galaxy at z=3 using "ErisMC", a new cosmological hydrodynamic "zoom-in" simulations of a disk galaxy with mass comparable to the Milky Way. The run adopts a blastwave scheme for supernova feedback that generates galactic outflows without explicit wind particles, a star formation recipe based on a high gas density threshold, and high temperature metal cooling. ErisMC's main progenitor at z=3 resembles a "Lyman break" galaxy of mass M_vir=2.4e11 M_sun, virial radius R_vir=48 kpc, and star formation rate 18 M_sun/yr, and its metal-enriched CGM extends to 200 (physical) kpc from its center. Approximately 41, 9, and 50 percent of all gas-phase metals at z=3 are locked in a hot (T> 3e5 K), warm (3e5 >T> 3e4 K), and cold (T< 3e4 K) medium, respectively. We identify three sources of heavy elements: 1) the main host, responsible for 60% of all the metals found within 3R_vir; 2) its satellite progenitors, responsible for 28% of all the metals within 3R_vir, and for only 5% of those beyond 3R_vir; and nearby dwarfs, which give origin to 12% of all the metals within 3R_vir and 95% of those beyond 3R_vir. Late (z<5) galactic "superwinds" account for only 9% of all the metals observed beyond 2R_vir, the bulk having been released at redshifts 5< z < 8 by early star formation and outflows. In the CGM, lower overdensities are typically enriched by `older', colder metals. Heavy elements are accreted onto Eris along filaments via low-metallicity cold inflows, and are ejected hot via galactic outflows at a few hundred km/s. The outflow mass-loading factor is of order unity for the main halo, but can exceed 10 for nearby dwarfs. We stress that our "zoom-in" simulation focuses on the CGM of a single massive system and cannot describe the enrichment history of the IGM as a whole. (abridged)

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

Observing the End of Cold Flow Accretion using Halo Absorption Systems

We use cosmological SPH simulations to study the cool, accreted gas in two Milky Way-size galaxies through cosmic time to z=0. We find that gas from mergers and cold flow accretion results in significant amounts of cool gas in galaxy halos. This cool circum-galactic component drops precipitously once the galaxies cross the critical mass to form stable shocks, Mvir = Msh ~ 10^12 Msun. Before reaching Msh, the galaxies experience cold mode accretion (T<10^5 K) and show moderately high covering fractions in accreted gas: f_c ~ 30-50% for R<50 co-moving kpc and N_HI>10^16 cm^-2. These values are considerably lower than observed covering fractions, suggesting that outflowing gas (not included here) is important in simulating galaxies with realistic gaseous halos. Within ~500 Myr of crossing the Msh threshold, each galaxy transitions to hot mode gas accretion, and f_c drops to ~5%. The sharp transition in covering fraction is primarily a function of halo mass, not redshift. This signature should be detectable in absorption system studies that target galaxies of varying host mass, and may provide a direct observational tracer of the transition from cold flow accretion to hot mode accretion in galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Orbiting Circum-galactic Gas as a Signature of Cosmological Accretion

We use cosmological SPH simulations to study the kinematic signatures of cool gas accretion onto a pair of well-resolved galaxy halos. Cold-flow streams and gas-rich mergers produce a circum-galactic component of cool gas that generally orbits with high angular momentum about the galaxy halo before falling in to build the disk. This signature of cosmological accretion should be observable using background-object absorption line studies as features that are offset from the galaxy's systemic velocity by ~100 km/s. Accreted gas typically co-rotates with the central disk in the form of a warped, extended cold flow disk, such that the observed velocity offset is in the same direction as galaxy rotation, appearing in sight lines that avoid the galactic poles. This prediction provides a means to observationally distinguish accreted gas from outflow gas: the accreted gas will show large one-sided velocity offsets in absorption line studies while radial/bi-conical outflows will not (except possibly in special polar projections). This rotation signature has already been seen in studies of intermediate redshift galaxy-absorber pairs; we suggest that these observations may be among the first to provide indirect observational evidence for cold accretion onto galactic halos. Cold mode halo gas typically has ~3-5 times more specific angular momentum than the dark matter. The associated cold mode disk configurations are likely related to extended HI/XUV disks seen around galaxies in the local universe. The fraction of galaxies with extended cold flow disks and associated offset absorption-line gas should decrease around bright galaxies at low redshift, as cold mode accretion dies out.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Probing the shape and history of the Milky Way halo with orbital spectral analysis

The phase space coordinates of individual halo stars obtained by Galactic surveys enable the computation of their full 3-dimensional orbits. Spectral analysis of halo orbits can be used to construct "frequency maps" which provide a compact representation of the 6-dimensional phase space distribution function. Frequency maps identify important major orbit families, and the orbital abundances reflect the shape and orientation of the dark matter halo relative to the disk. We apply spectral analysis to halo orbits in a series of controlled simulations of disk galaxies. Although the shape of the simulated halo varies with radius, frequency maps of local samples of halo orbits confined to the inner halo contain most of the information about the global shape of the halo and its major orbit families. Quiescent or adiabatic disk formation results in significant trapping of halo orbits in resonant orbit families (i.e. orbits with commensurable frequencies). If a good estimate of the Galactic potential in the inner halo (within ~50kpc) is available, the appearance of strong, stable resonances in frequency maps of halo orbits will allow us to determine the degree of resonant trapping induced by the disk potential. The locations and strengths of these resonant families are determined both by the global shape of the halo and its distribution function. Identification of such resonances in the Milky Way's stellar halo would therefore provide evidence of an extended period of adiabatic disk growth. If the Galactic potential is not known exactly, a measure of the diffusion rate of large sample of 10^4 halo orbits can help distinguish between the true potential and an incorrect potential. The orbital spectral analysis methods described in this paper provide a strong complementarity to existing methods for constraining the potential of the Milky Way halo and its stellar distribution function (ABRIDGED).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Properties of Carbon-Oxygen White Dwarf Merger Remnants

Recent studies have shown that for suitable initial conditions both super- and sub-Chandrasekhar mass carbon-oxygen white dwarf mergers produce explosions similar to observed SNe Ia. The question remains, however, how much fine tuning is necessary to produce these conditions. We performed a large set of SPH merger simulations, sweeping the possible parameter space. We find trends for merger remnant properties, and discuss how our results affect the viability of our recently proposed sub-Chandrasekhar merger channel for SNe Ia.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The First Massive Black Hole Seeds and Their Hosts

We investigate the formation of the first massive black holes in high redshift galaxies, with the goal of providing insights to which galaxies do or do not host massive black holes. We adopt a novel approach to forming seed black holes in galaxy halos in cosmological SPH+N-body simulations. The formation of massive black hole seeds is dictated directly by the local gas density, temperature, and metallicity, and motivated by physical models of massive black hole formation. We explore seed black hole populations as a function of halo mass and redshift, and examine how varying the efficiency of massive black hole seed formation affects the relationship between black holes and their hosts. Seed black holes tend to form in halos with mass between 10^7 and 10^9 Msun, and the formation rate is suppressed around z = 5 due to the diffusion of metals throughout the intergalactic medium. We find that the time of massive black hole formation and the occupation fraction of black holes are a function of the host halo mass. By z = 5, halos with mass M_halo > 3 x 10^9 Msun host massive black holes regardless of the efficiency of seed formation, while the occupation fraction for smaller halos increases with black hole formation efficiency. Our simulations explain why massive black holes are found in some bulgeless and dwarf galaxies, but we also predict that their occurrence becomes rarer and rarer in low-mass systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The fragmentation of protostellar discs: the Hill criterion for spiral arms

We present a new framework to explain the link between cooling and fragmentation in gravitationally unstable protostellar discs. This framework consists of a simple model for the formation of spiral arms, as well as a criterion, based on the Hill radius, to determine if a spiral arm will fragment. This detailed model of fragmentation is based on the results of numerical simulations of marginally stable protostellar discs, including those found in the literature, as well as our new suite of 3-D radiation hydrodynamics simulations of an irradiated, optically-thick protostellar disc surrounding an A star. Our set of simulations probes the transition to fragmentation through a scaling of the physical opacity. This model allows us to directly calculate the critical cooling time of Gammie (2001), with results that are consistent with those found from numerical experiment. We demonstrate how this model can be used to predict fragmentation in irradiated protostellar discs. These numerical simulations, as well as the model that they motivate, provide strong support for the hypothesis that gravitational instability is responsible for creating systems with giant planets on wide orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Consequences of cosmic microwave background-regulated star formation

It has been hypothesized that the cosmic microwave background (CMB) provides a temperature floor for collapsing protostars that can regulate the process of star formation and result in a top-heavy initial mass function at high metallicity and high redshift. We examine whether this hypothesis has any testable observational consequences. First we determine, using a set of hydrodynamic galaxy formation simulations, that the CMB temperature floor would have influenced the majority of stars formed at redshifts between z=3 and 6, and probably even to higher redshift. Five signatures of CMB-regulated star formation are: (1) a higher supernova rate than currently predicted at high redshift; (2) a systematic discrepancy between direct and indirect measurements of the high redshift star formation rate; (3) a lack of surviving globular clusters that formed at high metallicity and high redshift; (4) a more rapid rise in the metallicity of cosmic gas than is predicted by current simulations; and (5) an enhancement in the abundances of alpha elements such as O and Mg at metallicities -2 < [Fe/H] < -0.5. Observations are not presently able to either confirm or rule out the presence of these signatures. However, if correct, the top-heavy IMF of high-redshift high-metallicity globular clusters could provide an explanation for the observed bimodality of their metallicity distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Interpreting the Evolution of the Size - Luminosity Relation for Disk Galaxies from Redshift 1 to the Present

A sample of very high resolution cosmological disk galaxy simulations is used to investigate the evolution of galaxy disk sizes back to redshift 1 within the Lambda CDM cosmology. Artificial images in the rest frame B band are generated, allowing for a measurement of disk scale lengths using surface brightness profiles as observations would, and avoiding any assumption that light must follow mass as previous models have assumed. We demonstrate that these simulated disks are an excellent match to the observed magnitude - size relation for both local disks, and for disks at z=1 in the magnitude/mass range of overlap. We disentangle the evolution seen in the population as a whole from the evolution of individual disk galaxies. In agreement with observations, our simulated disks undergo roughly 1.5 magnitudes/arcsec^2 of surface brightness dimming since z=1. We find evidence that evolution in the magnitude - size plane varies by mass, such that galaxies with M* > 10^9 M_sun undergo more evolution in size than luminosity, while dwarf galaxies tend to evolve potentially more in luminosity. The disks grow in such a way as to stay on roughly the same stellar mass - size relation with time. Finally, due to an evolving stellar mass - SFR relation, a galaxy at a given stellar mass (or size) at z=1 will reside in a more massive halo and have a higher SFR, and thus a higher luminosity, than a counterpart of the same stellar mass at z=0.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Misaligned Angular Momentum in Hydrodynamic Cosmological Simulations: Warps, Outer Discs, and Thick Discs

We present a detailed analysis of a disc galaxy forming in a high-resolution fully cosmological simulation to investigate the nature of the outer regions of discs and their relevance for the disc formation process. Specifically, we focus on the phenomenon of misaligned disc components and find that the outer disc warp is a consequence of the misalignment between the inner disc and the surrounding hot gaseous halo. As the infalling cold gas sinks toward the centre of the galaxy, it is strongly torqued by the hot gas halo. By the time the fresh gas reaches the central disc-forming region its angular momentum is completely aligned with the spin of the hot gas halo. If the spin of the hot gas halo, in turn, is not aligned with that of the inner disc, a misaligned outer disc forms comprised of newly accreted material. The inner and outer components are misaligned with each other because they respond differently to infalling substructure and accretion. The warped disc feeds the main gas disc due to viscous angular momentum losses, but small amounts of star formation in the warp itself form a low-metallicity thick disc. We show that observations of resolved stellar populations in warped galaxies in the local universe could provide evidence for the presence of these processes and therefore indirectly reveal ongoing gas accretion and the existence of hot gas halos.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Star Formation and Feedback in Smoothed Particle Hydrodynamic Simulations II: Resolution Effects

We examine the effect of mass and force resolution on a specific star formation (SF) recipe using a set of N-body/Smooth Particle Hydrodynamic simulations of isolated galaxies. Our simulations span halo masses from 10^9 to 10^13 solar masses, more than four orders of magnitude in mass resolution, and two orders of magnitude in the gravitational softening length, epsilon, representing the force resolution. We examine the total global star formation rate, the star formation history, and the quantity of stellar feedback and compare the disk structure of the galaxies. Based on our analysis, we recommend using at least 10^4 particles each for the dark matter and gas component and a force resolution of epsilon approximately equal to 10^-3 R_vir when studying global SF and feedback. When the spatial distribution of stars is important, the number of gas and dark matter particles must be increased to at least 10^5 of each. Low mass resolution simulations with fixed softening lengths show particularly weak stellar disks due to two-body heating. While decreasing spatial resolution in low mass resolution simulations limits two-body effects, density and potential gradients cannot be sustained. Regardless of the softening, low-mass resolution simulations contain fewer high density regions where SF may occur. Galaxies of approximately 10^10 solar masses display unique sensitivity to both mass and force resolution. This mass of galaxy has a shallow potential and is on the verge of forming a disk. The combination of these factors give this galaxy the potential for strong gas outflows driven by supernova feedback and make it particularly sensitive to any changes to the simulation parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Dual Origin of Stellar Halos II: Chemical Abundances as Tracers of Formation History

Fully cosmological, high resolution N-Body + SPH simulations are used to investigate the chemical abundance trends of stars in simulated stellar halos as a function of their origin. These simulations employ a physically motivated supernova feedback recipe, as well as metal enrichment, metal cooling and metal diffusion. As presented in an earlier paper, the simulated galaxies in this study are surrounded by stellar halos whose inner regions contain both stars accreted from satellite galaxies and stars formed in situ in the central regions of the main galaxies and later displaced by mergers into their inner halos. The abundance patterns ([Fe/H] and [O/Fe]) of halo stars located within 10 kpc of a solar-like observer are analyzed. We find that for galaxies which have not experienced a recent major merger, in situ stars at the high [Fe/H] end of the metallicity distribution function are more [alpha/Fe]-rich than accreted stars at similar [Fe/H]. This dichotomy in the [O/Fe] of halo stars at a given [Fe/H] results from the different potential wells within which in situ and accreted halo stars form. These results qualitatively match recent observations of local Milky Way halo stars. It may thus be possible for observers to uncover the relative contribution of different physical processes to the formation of stellar halos by observing such trends in the halo populations of the Milky Way, and other local L* galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Enrichment of the Intergalactic Medium with Adiabatic Feedback I: Metal Cooling and Metal Diffusion

A study of the IGM metal enrichment using a series of SPH simulations is presented, employing metal cooling and turbulent diffusion of metals and thermal energy. An adiabatic feedback mechanism was adopted where gas cooling was prevented to generate galactic winds without explicit wind particles. The simulations produced a cosmic star formation history (SFH) that is broadly consistent with observations until z $\sim$ 0.5, and a steady evolution of the universal neutral hydrogen fraction ($Ω_{\rm H I}$). At z=0, about 40% of the baryons are in the warm-hot intergalactic medium (WHIM), but most metals (80%-90%) are locked in stars. At higher z the proportion of metals in the IGM is higher due to more efficient loss from galaxies. The IGM metals primarily reside in the WHIM throughout cosmic history. The metallicity evolution of the gas inside galaxies is broadly consistent with observations, but the diffuse IGM is under enriched at z $\sim$ 2.5. Galactic winds most efficiently enrich the IGM for halos in the intermediate mass range $10^{10}$M$_{\sun}$ - $10^{11}$ M$_{\sun}$. At the low mass end gas is prevented from accreting onto halos and has very low metallicities. At the high mass end, the fraction of halo baryons escaped as winds declines along with the decline of stellar mass fraction of the galaxies. This is likely because of the decrease in star formation activity and in wind escape efficiency. Metals enhance cooling which allows WHIM gas to cool onto galaxies and increases star formation. Metal diffusion allows winds to mix prior to escape, decreasing the IGM metal content in favour of gas within galactic halos and star forming gas. Diffusion significantly increases the amount of gas with low metallicities and changes the density-metallicity relation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Genesis of the Milky Way's Thick Disk via Stellar Migration

The separation of the Milky Way disk into a thin and thick component is supported by differences in the spatial, kinematic and metallicity distributions of their stars. These differences have led to the view that the thick disk formed early via a cataclysmic event and constitutes fossil evidence of the hierarchical growth of the Milky Way. We show here, using N-body simulations, how a double-exponential vertical structure, with stellar populations displaying similar dichotomies can arise purely through internal evolution. In this picture, stars migrate radially, while retaining nearly circular orbits, as described by Sellwood & Binney (2002). As stars move outwards they populate a thickened component. Such stars found at the present time in the solar neighborhood formed early in the disk's history at smaller radii where stars are more metal-poor and alpha-enhanced, leading to the properties observed for thick disk stars. Classifying stars as members of the thin or thick disk by either velocity or metallicity leads to an apparent separation in the other property as observed. This scenario is supported by the SDSS observation that stars in the transition region do not show any correlation between rotational velocity and metallicity. The good qualitative agreement between our simulation and observations in the Milky Way hint that the thick disk may be a ubiquitous galaxy feature generated by stellar migration. Nonetheless, we cannot exclude that some fraction of the thick disk is a fossil of a past more violent history, nor can this scenario explain thick disks in all galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Wandering Black Holes in Bright Disk Galaxy Halos

We perform SPH+N-body cosmological simulations of massive disk galaxies, including a formalism for black hole seed formation and growth, and find that satellite galaxies containing supermassive black hole seeds are often stripped as they merge with the primary galaxy. These events naturally create a population of "wandering" black holes that are the remnants of stripped satellite cores; galaxies like the Milky Way may host 5 -- 15 of these objects within their halos. The satellites that harbor black hole seeds are comparable to Local Group dwarf galaxies such as the Small and Large Magellanic Clouds; these galaxies are promising candidates to host nearby intermediate mass black holes. Provided that these wandering black holes retain a gaseous accretion disk from their host dwarf galaxy, they give a physical explanation for the origin and observed properties of some recently discovered off-nuclear ultraluminous X-ray sources such as HLX-1.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Numerical Study of Brown Dwarf Formation via Encounters of Protostellar Disks

The formation of brown dwarfs (BDs) due to the fragmentation of proto-stellar disks undergoing pairwise encounters was investigated. High resolution allowed the use of realistic initial disk models where both the vertical structure and the local Jeans mass were resolved. The results show that objects with masses ranging from giant planets to low mass stars can form during such encounters from initially stable disks. The parameter space of initial spin-orbit orientations and the azimuthal angles for each disk was explored. An upper limit on the initial Toomre Q value of ~2 was found for fragmentation to occur. Depending on the initial configuration, shocks, tidal-tail structures and mass inflows were responsible for the condensation of disk gas. Retrograde disks were generally more likely to fragment. When the interaction timescale was significantly shorter than the disks' dynamical timescales, the proto-stellar disks tended to be truncated without forming objects. The newly-formed objects had masses ranging from 0.9 to 127 Jupiter masses, with the majority in the BD regime. They often resided in star-BD multiples and in some cases also formed hierarchical orbiting systems. Most of them had large angular momenta and highly flattened, disk-like shapes. The objects had radii ranging from 0.1 to 10 AU. The disk gas was assumed to be locally isothermal, appropriate for the short cooling times in extended proto-stellar disks, but not for condensed objects. An additional case with explicit cooling that reduced to zero for optically thick gas was simulated to test the extremes of cooling effectiveness and it was still possible to form objects in this case. Detailed radiative transfer is expected to lengthen the internal evolution timescale for these objects, but not to alter our basic results.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

At the heart of the matter: the origin of bulgeless dwarf galaxies and Dark Matter cores

For almost two decades the properties of "dwarf" galaxies have challenged the Cold Dark Matter (CDM) paradigm of galaxy formation. Most observed dwarf galaxies consists of a rotating stellar disc embedded in a massive DM halo with a near constant-density core. Yet, models based on the CDM scenario invariably form galaxies with dense spheroidal stellar "bulges" and steep central DM profiles, as low angular momentum baryons and DM sink to the center of galaxies through accretion and repeated mergers. Processes that decrease the central density of CDM halos have been identified, but have not yet reconciled theory with observations of present day dwarfs. This failure is potentially catastrophic for the CDM model, possibly requiring a different DM particle candidate. This Letter presents new hydrodynamical simulations in a Lambda$CDM framework where analogues of dwarf galaxies, bulgeless and with a shallow central DM profile, are formed. This is achieved by resolving the inhomogeneous interstellar medium, resulting in strong outflows from supernovae explosions which remove low angular momentum gas. This inhibits the formation of bulges and decreases the dark-matter density to less than half within the central kiloparsec. Realistic dwarf galaxies are thus shown to be a natural outcome of galaxy formation in the CDM scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Feedback and the Formation of Dwarf Galaxy Stellar Halos

Stellar population studies show that low mass galaxies in all environments exhibit stellar halos that are older and more spherically distributed than the main body of the galaxy. In some cases, there is a significant intermediate age component that extends beyond the young disk. We examine a suite of Smoothed Particle Hydrodynamic (SPH) simulations and find that elevated early star formation activity combined with supernova feedback can produce an extended stellar distribution that resembles these halos for model galaxies ranging from $v_{200}$ = 15 km s$^{-1}$ to 35 km s$^{-1}$, without the need for accretion of subhalos.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Pairing of Supermassive Black Holes in unequal-mass galaxy mergers

We examine the pairing process of supermassive black holes (SMBHs) down to scales of 20-100 pc using a set of N-body/SPH simulations of binary mergers of disk galaxies with mass ratios of 1:4 and 1:10. Our numerical experiments are designed to represent merger events occurring at various cosmic epochs. The initial conditions of the encounters are consistent with the LambdaCDM paradigm of structure formation, and the simulations include the effects of radiative cooling, star formation, and supernovae feedback. We find that the pairing of SMBHs depends sensitively on the amount of baryonic mass preserved in the center of the companion galaxies during the last phases of the merger. In particular, due to the combination of gasdynamics and star formation, we find that a pair of SMBHs can form in 1:10 minor mergers provided that galaxies are relatively gas-rich (gas fractions of 30% of the disk mass) and that the mergers occur at relatively high redshift (z~3), when dynamical friction timescales are shorter. Since 1:10 mergers are most common events during the assembly of galaxies, and mergers are more frequent at high redshift when galaxies are also more gas-rich, our results have positive implications for future gravitational wave experiments such as the Laser Interferometer Space Antenna.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Redistributing hot gas around galaxies: do cool clouds signal a solution to the overcooling problem?

We present a pair of high-resolution smoothed particle hydrodynamics (SPH) simulations that explore the evolution and cooling behavior of hot gas around Milky-Way size galaxies. The simulations contain the same total baryonic mass and are identical other than their initial gas density distributions. The first is initialised with a low entropy hot gas halo that traces the cuspy profile of the dark matter, and the second is initialised with a high-entropy hot halo with a cored density profile as might be expected in models with pre-heating feedback. Galaxy formation proceeds in dramatically different fashion depending on the initial setup. While the low-entropy halo cools rapidly, primarily from the central region, the high-entropy halo is quasi-stable for ~4 Gyr and eventually cools via the fragmentation and infall of clouds from ~100 kpc distances. The low-entropy halo's X-ray surface brightness is ~100 times brighter than current limits and the resultant disc galaxy contains more than half of the system's baryons. The high-entropy halo has an X-ray brightness that is in line with observations, an extended distribution of pressure-confined clouds reminiscent of observed populations, and a final disc galaxy that has half the mass and ~50% more specific angular momentum than the disc formed in the low-entropy simulation. The final high-entropy system retains the majority of its baryons in a low-density hot halo. The hot halo harbours a trace population of cool, mostly ionised, pressure-confined clouds that contain ~10% of the halo's baryons after 10 Gyr of cooling. The covering fraction for HI and MgII absorption clouds in the high-entropy halo is ~0.4 and ~0.6, respectively, although most of the mass that fuels disc growth is ionised, and hence would be under counted in HI surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The nature of HI absorbers in GRB afterglows: clues from hydrodynamic simulations

In recent work, we have shown that it is possible to link quantitatively many aspects of damped Lyman alpha (DLA) absorbers in the spectra of quasars to high resolution simulations of galaxy formation. Using runs from the same series of hydrodynamic numerical studies, we consider the expected properties of Lyman alpha absorbers seen in the spectra of high redshift (z>2) gamma ray burst afterglows (GRB-DLAs). If GRBs are associated with the death of massive stars, their afterglows provide insights into otherwise unprobed regions of protogalactic objects, but detailed physical interpretations are currently embryonic. We find that median impact parameters (measured from the potential minimum) are approximately 1 kpc for GRBs compared with 4 kpc for QSO-DLAs. However, an equally important difference is that GRB-DLAs are predominantly associated with halos of mass 10^10<M_vir/M_sol<10^12, an order of magnitude larger than the hosts of QSO-DLAs. Accordingly, there are differences in the stellar properties of hosts. Our simulations accurately predict the form of the GRB-DLA HI column density distribution, producing quantitative agreement for N_HI>10^19 cm^-2, but they somewhat underpredict the incidence of low column densities N_HI<10^19 cm^-2. Line-of-sight neutral gas metallicities predicted by our simulations (10^-2 < Z/Z_sol < 1) are consistent with the modest observational constraints. Because of large internal dispersions in gas metallicities, this agreement is not significantly compromised by imposing a cut-off on the metallicity of stars able to launch GRBs (Z_star<Z_sol/3), confounding claims that the observed metallicity of GRB-DLAs poses a challenge to current GRB models. (Abridged.)

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Beyond Inside-Out Growth: Formation and Evolution of Disk Outskirts

We have performed a high mass and force resolution simulation of an idealized galaxy forming from dissipational collapse of gas embedded in a spherical dark matter halo. The simulation includes star formation and effects of stellar feedback. In our simulation a stellar disk forms with a surface density profile consisting of an inner exponential breaking to a steeper outer exponential. The break forms early on and persists throughout the evolution, moving outward as more gas is able to cool and add mass to the disk. The parameters of the break are in excellent agreement with observations. The break corresponds with a rapid drop in the star formation rate associated with a drop in the cooled gas surface density, but the outer exponential is populated by stars that were scattered outward on nearly circular orbits from the inner disk by spiral arms. The resulting profile and its associated break are therefore a consequence of the interplay between a radial star formation cutoff and redistribution of stellar mass by secular processes. A consequence of such evolution is a sharp change in the radial mean stellar age profile at the break radius.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Early Gas Stripping as the Origin of the Darkest Galaxies in the Universe

The known galaxies most dominated by dark matter (Draco, Ursa Minor and Andromeda IX) are satellites of the Milky Way and the Andromeda galaxies. They are members of a class of faint galaxies, devoid of gas, known as dwarf spheroidals, and have by far the highest ratio of dark to luminous matter. None of the models proposed to unravel their origin can simultaneously explain their exceptional dark matter content and their proximity to a much larger galaxy. Here we report simulations showing that the progenitors of these galaxies were probably gas-dominated dwarf galaxies that became satellites of a larger galaxy earlier than the other dwarf spheroidals. We find that a combination of tidal shocks and ram pressure swept away the entire gas content of such progenitors about ten billion years ago because heating by the cosmic ultraviolet background kept the gas loosely bound: a tiny stellar component embedded in a relatively massive dark halo survived until today. All luminous galaxies should be surrounded by a few extremely dark-matter-dominated dwarf spheroidal satellites, and these should have the shortest orbital periods among dwarf spheroidals because they were accreted early.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Multi-scale simulations of merging galaxies with supermassive black holes

We present the results of the first multi-scale N-Body+SPH simulations of merging galaxies containing central supermassive black holes (SMBHs) and having a spatial resolution of only a few parsecs. Strong gas inflows associated with equal-mass mergers produce non-axisymmetric nuclear disks with masses of order $10^9 M_{\odot}$, resolved by about $10^6$ SPH particles. Such disks have sizes of several hundred parsecs but most of their mass is concentrated within less than $50$ pc. We find that a close SMBH pair forms after the merger. The separation of the two SMBHs then shrinks further owing to dynamical friction against the predominantly gaseous background. The orbits of the SMBHs decay down to the minimum resolvable scale in a few million years for an ambient gas temperature and density typical of a region undergoing a starburst. These results suggest the initial conditions necessary for the eventual coalescence of the two holes arise naturally from the merging of two equal-mass galaxies whose structure and orbits are consistent with the predictions of the $Λ$CDM model. Our findings have important implications for planned gravitational wave detection experiments such as {\it LISA}.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Proto-Brown Dwarf Disks as Products of Protostellar Disk Encounters

The formation of brown dwarfs via encounters between proto-stars has been confirmed with high-resolution numerical simulations with a restricted treatment of the thermal conditions. The new results indicate that young brown dwarfs (BDs) formed this way are disk-like and often reside in multiple systems. The newly-formed proto-BDs disks are up to 18 AU in size and spin rapidly making small-scale bipolar outflows, fragmentation and the possible formation of planetary companions likely as have recently been observed for BDs. The object masses range from 2 to 73 Jupiter masses, distributed in a manner consistent with the observed sub-stellar initial mass function. The simulations usually form multiple BDs on eccentric orbits about a star. One such system was hierarchical, a BD binary in orbit around a star, which may explain recently observed hierarchical systems. One third of the BDs were unbound after a few thousand years and interactions among orbiting BDs may eject more or add to the number of binaries. Improvements over prior work include resolution down to a Jupiter mass, self-consistent models of the vertical structure of the initial disks and careful attention to avoid artificial fragmentation.

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
preprint2003arXiv

Gasoline: An adaptable implementation of TreeSPH

The key algorithms and features of the Gasoline code for parallel hydrodynamics with self-gravity are described. Gasoline is an extension of the efficient Pkdgrav parallel N-body code using smoothed particle hydrodynamics. Accuracy measurements, performance analysis and tests of the code are presented. Recent successful Gasoline applications are summarized. These cover a diverse set of areas in astrophysics including galaxy clusters, galaxy formation and gas-giant planets. Future directions for gasdynamical simulations in astrophysics and code development strategies for tackling cutting edge problems are discussed.

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

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