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

J. Wadsley

J. 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-ph.SRastro-phastro-ph.EPastro-ph.HEphysics.flu-dyn

Catalog footprint

What is connected

38works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

38 published item(s)

preprint2021arXiv

Cosmic filaments delay quenching inside clusters

We investigate how large-scale cosmic filaments impact the quenching of galaxies within one virial radius of 324 simulated clusters from The Three Hundred project. We track cosmic filaments with the versatile, observation-friendly program DisPerSE and identify halos hosting galaxies with VELOCIRaptor. We find that cluster galaxies close to filaments tend to be more star-forming, bluer, and contain more cold gas than their counterparts further away from filaments. This effect is recovered at all stellar masses. This is in stark contrast with galaxies residing outside of clusters, where galaxies close to filaments show clear signs of density related pre-processing. We first show that the density contrast of filaments is reduced inside the intra-cluster medium. Moreover, examination of flows around and into cluster galaxies shows that the gas flows in intra-cluster filaments are colder and tend to stream along with galaxies in their midst, partially shielding them from strangulation by the hot, dense intra-cluster medium. This also preserves accretion onto satellites and limit ram pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Magnetorotational instability with smoothed particle hydrodynamics

We present a thorough numerical study on the MRI using the smoothed particle magnetohydrodynamics method (SPMHD) with the geometric density average force expression (GDSPH). We perform shearing box simulations with different initial setups and a wide range of resolution and dissipation parameters. We show, for the first time, that MRI with sustained turbulence can be simulated successfully with SPH, with results consistent with prior work with grid-based codes. In particular, for the stratified boxes, our simulations reproduce the characteristic butterfly diagram of the MRI dynamo with saturated turbulence for at least 100 orbits. On the contrary, traditional SPH simulations suffer from runaway growth and develop unphysically large azimuthal fields, similar to the results from a recent study with mesh-less methods. We investigated the dependency of MRI turbulence on the numerical Prandtl number in SPH, focusing on the unstratified, zero net-flux case. We found that turbulence can only be sustained with a Prandtl number larger than $\sim$2.5, similar to the critical values of physical Prandtl number found in grid-code simulations. However, unlike grid-based codes, the numerical Prandtl number in SPH increases with resolution, and for a fixed Prandtl number, the resulting magnetic energy and stresses are independent of resolution. Mean-field analyses were performed on all simulations, and the resulting transport coefficients indicate no $α$-effect in the unstratified cases, but an active $αΩ$ dynamo and a diamagnetic pumping effect in the stratified medium, which are generally in agreement with previous studies. There is no clear indication of a shear-current dynamo in our simulation, which is likely to be responsible for a weaker mean-field growth in the tall, unstratified, zero net-flux simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmological Galaxy Evolution with Superbubble Feedback II: The Limits of Supernovae

We explore when supernovae can (and cannot) regulate the star formation and bulge growth in galaxies based on a sample of 18 simulated galaxies. The simulations include key physics such as evaporation and conduction, neglected in prior work, and required to correctly model superbubbles resulting from stellar feedback. We show that for galaxies with virial masses $>10^{12}\;M_\odot$, supernovae alone cannot prevent excessive star formation. This failure occurs due to a shutdown of galactic winds, with wind mass loadings falling from $η\sim10$ to $η<1$. In more massive systems, this transfer of baryons to the circumgalactic medium falters earlier on and the galaxies diverge significantly from observed galaxy scaling relations and morphologies. The decreasing efficiency is simply due to a deepening potential well preventing gas escape. This implies that non-supernova feedback mechanisms must become dominant for galaxies with stellar masses greater than $\sim4\times10^{10}\;M_\odot$. The runaway growth of the central stellar bulge, strongly linked to black hole growth, suggests that feedback from active galactic nuclei is the probable mechanism. Below this mass, supernovae alone are able to produce a realistic stellar mass fraction, star formation history and disc morphology.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The MaGICC Baryon Cycle: The Enrichment History of Simulated Disc Galaxies

Using cosmological galaxy formation simulations from the MaGICC project, spanning more than three magnitudes in stellar mass (~10^7-3x10^{10} Msun), we trace the baryonic cycle of infalling gas from the virial radius through to its participation in the star formation process. An emphasis is placed upon the temporal history of chemical enrichment during its passage through the corona and CGM. We derive the distributions of time between gas crossing the virial radius and being accreted to the star forming region (which allows mixing within the corona), as well as the time between gas being accreted to the star forming region and then forming stars (which allows mixing within the disc). Significant numbers of stars are formed from gas that cycles back through the hot halo after first accreting to the star forming region. Gas entering high mass galaxies is pre-enriched in low mass proto-galaxies prior to entering the virial radius of the central progenitor, with only small amounts of primordial gas accreted, even at high redshift (z~5). After entering the virial radius, significant further enrichment occurs prior to the accretion of the gas to the star forming region, with gas that is feeding the star forming region surpassing 0.1Z by z=0. Mixing with halo gas, itself enriched via galactic fountains, is thus crucial in determining the metallicity at which gas is accreted to the disc. The lowest mass simulation (Mvir~2x10^{10}Msun, with M*~10^7Msun), by contrast, accretes primordial gas through the virial radius and onto the disc at all times. Much like the classical analytical solutions to the `G-dwarf problem', overproduction of low-metallicity stars is ameliorated by the inefficiency of star formation. Finally, gas outflow/metal removal rates from star forming regions as a function of galactic mass are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

All about baryons: revisiting SIDM predictions at small halo masses

We use cosmological hydrodynamic simulations to consistently compare the assembly of dwarf galaxies in both $Λ$ dominated, Cold (CDM) and Self--Interacting (SIDM) dark matter models. The SIDM model adopts a constant cross section of 2 $cm^{2}/g$, a relatively large value to maximize its effects. These are the first SIDM simulations that are combined with a description of stellar feedback that naturally drives potential fluctuations able to create dark matter cores. Remarkably, SIDM fails to significantly lower the central dark matter density at halo peak velocities V$_{max}$ $<$ 30 Km/s. This is due to the fact that the central regions of very low--mass field halos have relatively low central velocity dispersion and densities, leading to time scales for SIDM collisions greater than a Hubble time. CDM halos with V$_{max}$ $<$ 30 km/s have inefficient star formation, and hence weak supernova feedback. At a fixed 2 cm2/g SIDM cross section, the DM content of very low mass CDM and SIDM halos differs by no more than a factor of two within 100-200pc. At larger halo masses ($\sim$ 10$^{10}$ solar masses), the introduction of baryonic processes creates field dwarf galaxies with dark matter cores and central DM$+$baryon distributions that are effectively indistinguishable between CDM and SIDM. Both models are in broad agreement with observed Local Group field galaxies across the range of masses explored. To significantly differentiate SIDM from CDM at the scale of faint dwarf galaxies, a velocity dependent cross section that rapidly increases to values larger than 2 $cm^{2}/g$ for halos with V$_{max}$ < 25-30 Km/s needs to be introduced.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmological Galaxy Evolution with Superbubble Feedback I: Realistic Galaxies with Moderate Feedback

We present the first cosmological galaxy evolved using the modern smoothed particle hydrodynamics (SPH) code GASOLINE2 with superbubble feedback. We show that superbubble-driven galactic outflows powered by Type II supernovae alone can produce $\rm{L^*}$ galaxies with flat rotation curves with circular velocities $\sim 200\; \rm{km/s}$, low bulge-to-disc ratios, and stellar mass fractions that match observed values from high redshift to the present. These features are made possible by the high mass loadings generated by the evaporative growth of superbubbles. Outflows are driven extremely effectively at high redshift, expelling gas at early times and preventing overproduction of stars before $z=2$. Centrally concentrated gas in previous simulations has often lead to unrealistically high bulge to total ratios and strongly peaked rotation curves. We show that supernova-powered superbubbles alone can produce galaxies that agree well with observed properties without the need for additional feedback mechanisms or increased feedback energy. We present additional results arising from properly modelled hot feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

NIHAO III: The constant disc gas mass conspiracy

We show that the cool gas masses of galactic discs reach a steady state that lasts many Gyr after their last major merger in cosmological hydrodynamic simulations. The mass of disc gas, M$_{\rm gas}$, depends upon a galaxy halo's spin and virial mass, but not upon stellar feedback. Halos with low spin have high star formation efficiency and lower disc gas mass. Similarly, lower stellar feedback leads to more star formation so the gas mass ends up nearly the same irregardless of stellar feedback strength. Even considering spin, the M$_{\rm gas}$ relation with halo mass, M$_{200}$ only shows a factor of 3 scatter. The M$_{\rm gas}$--M$_{200}$ relation show a break at M$_{200}$=$2\times10^{11}$ M$_\odot$ that corresponds to an observed break in the M$_{\rm gas}$--M$_\star$ relation. The constant disc mass stems from a shared halo gas density profile in all the simulated galaxies. In their outer regions, the profiles are isothermal. Where the profile rises above $n=10^{-3}$ cm$^{-3}$, the gas readily cools and the profile steepens. Inside the disc, rotation supports gas with a flatter density profile except where supernova explosions disrupt the disc. Energy injection from stellar feedback also provides pressure support to the halo gas to prevent runaway cooling flows. The resulting constant gas mass makes simpler models for galaxy formation possible, either using a "bathtub" model for star formation rates or when modeling chemical evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The history of stellar metallicity in a simulated disc galaxy

We explore the chemical distribution of stars in a simulated galaxy. Using simulations of the same initial conditions but with two different feedback schemes (MUGS and MaGICC), we examine the features of the age-metallicity relation (AMR), and the three-dimensional age-metallicity-[O/Fe] distribution, both for the galaxy as a whole and decomposed into disc, bulge, halo, and satellites. The MUGS simulation, which uses traditional supernova feedback, is replete with chemical substructure. This sub- structure is absent from the MaGICC simulation, which includes early feedback from stellar winds, a modified IMF and more efficient feedback. The reduced amount of substructure is due to the almost complete lack of satellites in MaGICC. We identify a significant separation between the bulge and disc AMRs, where the bulge is considerably more metal-rich with a smaller spread in metallicity at any given time than the disc. Our results suggest, however, that identifying the substructure in observations will require exquisite age resolution, on the order of 0.25 Gyr. Certain satellites show exotic features in the AMR, even forming a 'sawtooth' shape of increasing metallicity followed by sharp declines which correspond to pericentric passages. This fact, along with the large spread in stellar age at a given metallicity, compromises the use of metallicity as an age indicator, although alpha abundance provides a more robust clock at early times. This may also impact algorithms that are used to reconstruct star formation histories from resolved stellar populations, which frequently assume a monotonically-increasing AMR.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Superbubble Feedback Model for Galaxy Simulations

We present a new stellar feedback model that reproduces superbubbles. Superbubbles from clustered young stars evolve quite differently to individual supernovae and are substantially more efficient at generating gas motions. The essential new components of the model are thermal conduction, sub-grid evaporation and a sub-grid multi-phase treatment for cases where the simulation mass resolution is insufficient to model the early stages of the superbubble. The multi-phase stage is short compared to superbubble lifetimes. Thermal conduction physically regulates the hot gas mass without requiring a free parameter. Accurately following the hot component naturally avoids overcooling. Prior approaches tend to heat too much mass, leaving the hot ISM below $10^6$ K and susceptible to rapid cooling unless ad-hoc fixes were used. The hot phase also allows feedback energy to correctly accumulate from multiple, clustered sources, including stellar winds and supernovae. We employ high-resolution simulations of a single star cluster to show the model is insensitive to numerical resolution, unresolved ISM structure and suppression of conduction by magnetic fields. We also simulate a Milky Way analog and a dwarf galaxy. Both galaxies show regulated star formation and produce strong outflows.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Orbital Decay of Supermassive Black Hole Binaries in Clumpy Multiphase Merger Remnants

We simulate an equal-mass merger of two Milky Way-size galaxy discs with moderate gas fractions at parsec-scale resolution including a new model for radiative cooling and heating in a multi-phase medium, as well as star formation and feedback from supernovae. The two discs initially have a $2.6\times10^6\mathrm{~M_{\odot}}$ supermassive black hole (SMBH) embedded in their centers. As the merger completes and the two galactic cores merge, the SMBHs form a a pair with a separation of a few hundred pc that gradually decays. Due to the stochastic nature of the system immediately following the merger, the orbital plane of the binary is significantly perturbed. Furthermore, owing to the strong starburst the gas from the central region is completely evacuated, requiring $\sim10$~Myr for a nuclear disc to rebuild. Most importantly, the clumpy nature of the interstellar medium has a major impact on the the dynamical evolution of the SMBH pair, which undergo gravitational encounters with massive gas clouds and stochastic torquing by both clouds and spiral modes in the disk. These effects combine to greatly delay the decay of the two SMBHs to separations of a few parsecs by nearly two orders of magnitude, $\sim 10^8$ yr, compared to previous work. In mergers of more gas-rich, clumpier galaxies at high redshift stochastic torques will be even more pronounced and potentially lead to stronger modulation of the orbital decay. This suggests that SMBH pairs at separations of several tens of parsecs should be relatively common at any redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating disc galaxy bulges that are consistent with observed scaling relations

We present a detailed comparison between the photometric properties of the bulges of two simulated galaxies and those of a uniform sample of observed galaxies. This analysis shows that the simulated galaxies have bulges with realistic surface brightnesses for their sizes and magnitude. These two field disc galaxies have rotational velocities ~ 100 km/s and were integrated to a redshift of zero in a fully cosmological Lambda cold dark matter context as part of high-resolution smoothed particle hydrodynamic simulations. We performed bulge-disc decompositions of the galaxies using artificial observations, in order to conduct a fair comparison to observations. We also dynamically decomposed the galaxies and compared the star formation histories of the bulges to those of the entire galaxies. These star formation histories showed that the bulges were primarily formed before z = 1 and during periods of rapid star formation. Both galaxies have large amounts of early star formation, which is likely related to the relatively high bulge-to-disc ratios also measured for them. Unlike almost all previous cosmological simulations, the realistically concentrated bulges of these galaxies do not lead to unphysically high rotational velocities, causing them to naturally lie along the observed Tully-Fisher relation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Galaxy Formation with local photoionisation feedback I. Methods

We present a first study of the effect of local photoionising radiation on gas cooling in smoothed particle hydrodynamics simulations of galaxy formation. We explore the combined effect of ionising radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm used for gravity, so it is computationally efficient and well resolved. The method foregoes calculating absorption and scattering in favour of a constant escape fraction for young stars to keep the calculation efficient enough to simulate the entire evolution of a galaxy in a cosmological context to the present day. This allows us to quantify the effect of the local photoionisation feedback through the whole history of a galaxy`s formation. The simulation of a Milky Way like galaxy using the local photoionisation model forms ~ 40 % less stars than a simulation that only includes a standard uniform background UV field. The local photoionisation model decreases star formation by increasing the cooling time of the gas in the halo and increasing the equilibrium temperature of dense gas in the disc. Coupling the local radiation field to gas cooling from the halo provides a preventive feedback mechanism which keeps the central disc light and produces slowly rising rotation curves without resorting to extreme feedback mechanisms. These preliminary results indicate that the effect of local photoionising sources is significant and should not be ignored in models of galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cuspy No More: How Outflows Affect the Central Dark Matter and Baryon Distribution in Lambda CDM Galaxies

We examine the evolution of the inner dark matter (DM) and baryonic density profile of a new sample of simulated field galaxies using fully cosmological, Lambda CDM, high resolution SPH + N-Body simulations. These simulations include explicit H2 and metal cooling, star formation (SF) and supernovae (SNe) driven gas outflows. Starting at high redshift, rapid, repeated gas outflows following bursty SF transfer energy to the DM component and significantly flatten the originally `cuspy' central DM mass profile of galaxies with present day stellar masses in the 10^4.5 -- 10^9.8 Msolar range. At z=0, the central slope of the DM density profile of our galaxies (measured between 0.3 and 0.7 kpc from their centre) is well fitted by rhoDM propto r^alpha with alpha \simeq -0.5 + 0.35 log_10(Mstar/10^8Msolar) where Mstar is the stellar mass of the galaxy and 4 < log_10 Mstar < 9.4. These values imply DM profiles flatter than those obtained in DM--only simulations and in close agreement with those inferred in galaxies from the THINGS and LITTLE THINGS survey. Only in very small halos, where by z=0 star formation has converted less than ~ 0.03% of the original baryon abundance into stars, outflows do not flatten the original cuspy DM profile out to radii resolved by our simulations. The mass (DM and baryonic) measured within the inner 500 pc of each simulated galaxy remains nearly constant over four orders of magnitudes in stellar mass for Mstar 10^9 Msolar. This finding is consistent with estimates for faint Local Group dwarfs and field galaxies. These results address one of the outstanding problems faced by the CDM model, namely the strong discrepancy between the original predictions of cuspy DM profiles and the shallower central DM distribution observed in galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

In situ formation of SgrA* stars via disk fragmentation: parent cloud properties and thermodynamics

The formation of the massive young stars surrounding SgrA* is still an open question. In this paper, we simulate the infall of a turbulent molecular cloud towards the Galactic Center (GC). We adopt two different cloud masses (4.3x10^4 and 1.3x10^5 solar masses). We run five simulations: the gas is assumed to be isothermal in four runs, whereas radiative cooling is included in the fifth run. In all the simulations, the molecular cloud is tidally disrupted, spirals towards the GC, and forms a small, dense and eccentric disk around SgrA*. With high resolution simulations, we follow the fragmentation of the gaseous disk. Star candidates form in a ring at ~0.1-0.4 pc from the super-massive black hole (SMBH) and have moderately eccentric orbits (~0.2-0.4), in good agreement with the observations. The mass function of star candidates is top-heavy only if the local gas temperature is high (>~100 K) during the star formation and if the parent cloud is sufficiently massive (>~10^5 solar masses). Thus, this study indicates that the infall of a massive molecular cloud is a viable scenario for the formation of massive stars around SgrA*, provided that the gas temperature is kept sufficiently high (>~100 K).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

MaGICC Disks: Matching Observed Galaxy Relationships Over a Wide Stellar Mass Range

We use the same physical model to simulate four galaxies that match the relation between stellar and total mass, over a mass range that includes the vast majority of disc galaxies. The resultant galaxies, part of the Making Galaxies in a Cosmological Context (MaGICC) program, also match observed relations between luminosity, rotation velocity, size, colour, star formation rate, HI mass, baryonic mass, and metallicity. Radiation from massive stars and supernova energy regulate star formation and drive outflows, balancing the complex interplay between cooling gas, star formation, large scale outflows, and recycling of gas in a manner which correctly scales with the mass of the galaxy. Outflows also play a key role in simulating galaxies with exponential surface brightness profiles, flat rotation curves and dark matter cores. Our study implies that large scale outflows are the primary driver of the dependence of disc galaxy properties on mass. We show that the amount of outflows invoked in our model is required to meet the constraints provided by observations of OVI absorption lines in the circum-galactic-media of local galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radial Migration in Disk Galaxies I: Transient Spiral Structure and Dynamics

We seek to understand the origin of radial migration in spiral galaxies by analyzing in detail the structure and evolution of an idealized, isolated galactic disk. To understand the redistribution of stars, we characterize the time-evolution of properties of spirals that spontaneously form in the disk. Our models unambiguously show that in such disks, single spirals are unlikely, but that a number of transient patterns may coexist in the disk. However, we also show that while spirals are transient in amplitude, at any given time the disk favors patterns of certain pattern speeds. Using several runs with different numerical parameters we show that the properties of spirals that occur spontaneously in the disk do not sensitively depend on resolution. The existence of multiple transient patterns has large implications for the orbits of stars in the disk, and we therefore examine the resonant scattering mechanisms that profoundly alter angular momenta of individual stars. We confirm that the corotation scattering mechanism described by Sellwood & Binney (2002) is responsible for the largest angular momentum changes in our simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Reproducing the Stellar Mass/Halo Mass Relation in Simulated LCDM Galaxies: Theory vs Observational Estimates

We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses ranging from a few times 10^8 to nearly 10^12 solar masses. At z=0 our simulated galaxies have a baryon content and morphology typical of field galaxies. Over a stellar mass range of 2.2 x 10^3 to 4.5 x 10^10 solar masses, we find extremely good agreement between the SHM ratio in simulations and the present-day predictions from the statistical Abundance Matching Technique presented in Moster et al. (2012). This improvement over past simulations is due to a number systematic factors, each decreasing the SHM ratios: 1) gas outflows that reduce the overall SF efficiency but allow for the formation of a cold gas component 2) estimating the stellar masses of simulated galaxies using artificial observations and photometric techniques similar to those used in observations and 3) accounting for a systematic, up to 30 percent overestimate in total halo masses in DM-only simulations, due to the neglect of baryon loss over cosmic times. Our analysis suggests that stellar mass estimates based on photometric magnitudes can underestimate the contribution of old stellar populations to the total stellar mass, leading to stellar mass errors of up to 50 percent for individual galaxies. These results highlight the importance of using proper techniques to compare simulations with observations and reduce the perceived tension between the star formation efficiency in galaxy formation models and in real galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The luminosity function of diverse satellite galaxy systems

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thin disc, Thick Disc and Halo in a Simulated Galaxy

Within a cosmological hydrodynamical simulation, we form a disc galaxy with sub- components which can be assigned to a thin stellar disc, thick disk, and a low mass stellar halo via a chemical decomposition. The thin and thick disc populations so selected are distinct in their ages, kinematics, and metallicities. Thin disc stars are young (<6.6 Gyr), possess low velocity dispersion (σU,V,W = 41, 31, 25 km/s), high [Fe/H], and low [O/Fe]. The thick disc stars are old (6.6<age<9.8 Gyrs), lag the thin disc by \sim21 km/s, possess higher velocity dispersion (σU,V,W = 49, 44, 35 km/s), relatively low [Fe/H] and high [O/Fe]. The halo component comprises less than 4% of stars in the "solar annulus" of the simulation, has low metallicity, a velocity ellipsoid defined by (σU,V,W = 62, 46, 45 km/s) and is formed primarily in-situ during an early merger epoch. Gas-rich mergers during this epoch play a major role in fuelling the formation of the old disc stars (the thick disc). This is consistent with studies which show that cold accretion is the main source of a disc galaxy's baryons. Our simulation initially forms a relatively short (scalelength \sim1.7 kpc at z=1) and kinematically hot disc, primarily from gas accreted during the galaxy's merger epoch. Far from being a competing formation scenario, migration is crucial for reconciling the short, hot, discs which form at high redshift in ΛCDM, with the properties of the thick disc at z=0. The thick disc, as defined by its abundances maintains its relatively short scale-length at z = 0 (2.31 kpc) compared with the total disc scale-length of 2.73 kpc. The inside-out nature of disc growth is imprinted the evolution of abundances such that the metal poor α-young population has a larger scale-length (4.07 kpc) than the more chemically evolved metal rich α-young population (2.74 kpc).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Advanced Morphological Galaxy Classification: A Comparison of Real and Simulated Galaxies

Encoded within the morphological structure of galaxies are clues related to their formation and evolutionary history. Recent advances pertaining to the statistics of galaxy morphology include sophisticated measures of concentration (C), asymmetry (A), and clumpiness (S). In this study, these three parameters (CAS) have been applied to a suite of simulated galaxies and compared with observational results inferred from a sample of nearby galaxies. The simulations span a range of late-type systems, with masses between ~1e10 Msun and ~1e12 Msun, and employ star formation density thresholds between 0.1 cm^-3 and 100 cm^-3. We have found that the simulated galaxies possess comparable concentrations to their real counterparts. However, the results of the CAS analysis revealed that the simulated galaxies are generally more asymmetric, and that the range of clumpiness values extends beyond the range of those observed. Strong correlations were obtained between the three CAS parameters and colour (B-V), consistent with observed galaxies. Furthermore, the simulated galaxies possess strong links between their CAS parameters and Hubble type, mostly in-line with their real counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Growing Massive Black Hole Pairs in Minor Mergers of Disk Galaxies

We perform a suite of high-resolution smoothed particle hydrodynamics simulations to investigate the orbital decay and mass evolution of massive black hole (MBH) pairs down to scales of ~30 pc during minor mergers of disk galaxies. Our simulation set includes star formation and accretion onto the MBHs, as well as feedback from both processes. We consider 1:10 merger events starting at z~3, with MBH masses in the sensitivity window of the Laser Interferometer Space Antenna, and we follow the coupling between the merger dynamics and the evolution of the MBH mass ratio until the satellite galaxy is tidally disrupted. While the more massive MBH accretes in most cases as if the galaxy were in isolation, the satellite MBH may undergo distinct episodes of enhanced accretion, owing to strong tidal torques acting on its host galaxy and to orbital circularization inside the disk of the primary galaxy. As a consequence, the initial 1:10 mass ratio of the MBHs changes by the time the satellite is disrupted. Depending on the initial fraction of cold gas in the galactic disks and the geometry of the encounter, the mass ratios of the MBH pairs at the time of satellite disruption can stay unchanged or become as large as 1:2. Remarkably, the efficiency of MBH orbital decay correlates with the final mass ratio of the pair itself: MBH pairs that increase significantly their mass ratio are also expected to inspiral more promptly down to nuclear-scale separations. These findings indicate that the mass ratios of MBH pairs in galactic nuclei do not necessarily trace the mass ratios of their merging host galaxies, but are determined by the complex interplay between gas accretion and merger dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hierarchical formation of bulgeless galaxies II: Redistribution of angular momentum via galactic fountains

Within a fully cosmological hydrodynamical simulation, we form a galaxy which rotates at 140 km/s, and is characterised by two loose spiral arms and a bar, indicative of a Hubble Type SBc/d galaxy. We show that our simulated galaxy has no classical bulge, with a pure disc profile at z=1, well after the major merging activity has ended. A long-lived bar subsequently forms, resulting in the formation of a secularly-formed "pseudo" bulge, with the final bulge-to-total light ratio B/T=0.21. We show that the majority of gas which loses angular momentum and falls to the central region of the galaxy during the merging epoch is blown back into the hot halo, with much of it returning later to form stars in the disc. We propose that this mechanism of redistribution of angular momentum via a galactic fountain, when coupled with the results from our previous study which showed why gas outflows are biased to have low angular momentum, can solve the angular momentum/bulgeless disc problem of the cold dark matter paradigm.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hierarchical formation of bulgeless galaxies: Why outflows have low angular momentum

Using high resolution, fully cosmological smoothed particle hydro-dynamical simulations of dwarf galaxies in a Lambda cold dark matter Universe, we show how baryons attain a final angular momentum distribution which allows pure disc galaxies to form. Blowing out substantial amounts of gas through supernovae and stellar winds, which is well supported observationally, is a key ingredient in forming bulgeless discs. We outline why galactic outflows preferentially remove low angular momentum material, and show that this is a natural result when structure forms in a cold dark matter cosmology. The driving factors are a) the mean angular momentum of accreted material increases with time, b) lower potentials at early times, c) the existence of an extended reservoir of high angular momentum gas which is not within star forming regions, meaning that only gas from the inner region (low angular momentum gas) is expelled and d) the tendency for outflows to follow the path of least resistance which is perpendicular to the disc. We also show that outflows are enhanced during mergers, thus expelling much of the gas which has lost its angular momentum during these events, and preventing the formation of "classical", merger driven bulges in low mass systems. Stars formed prior to such mergers form a diffuse, extended stellar halo component.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mechanisms of Baryon Loss for Dark Satellites in Cosmological SPH Simulations

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Cold Gas Content of Bulgeless Dwarf Galaxies

We present an analysis of the neutral hydrogen (HI) properties of a fully cosmological hydrodynamical dwarf galaxy, run with varying simulation parameters. As reported by Governato et al. (2010), the high resolution, high star formation density threshold version of this galaxy is the first simulation to result in the successful reproduction of a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to compare in detail the HI distribution and kinematics of this simulated bulgeless disk with what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas density profiles, velocity dispersion (e.g., velocity ellipsoid, turbulence), and the power spectrum of structure within the cold interstellar medium from the simulations. The highest resolution dwarf, when using a high density star formation threshold comparable to densities of giant molecular clouds, possesses bulk characteristics consistent with those observed in nature, though the cold gas is not as radially extended as that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The lines-of-sight velocity dispersion radial profiles have values that are in good agreement with observed dwarf galaxies, but due to the fact that only the streaming velocities of particles are tracked, a correction to include the thermal velocities can lead to profiles that are quite flat. The ISM power spectra of the simulations appear to possess more power on smaller spatial scales than that of the SMC. We conclude that unavoidable limitations remain due to the unresolved physics of star formation and feedback within pc-scale molecular clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Properties of Prestellar Discs in Isolated and Multiple Prestellar Systems

We present high-resolution 3D smoothed particle hydrodynamics simulations of the formation and evolution of protostellar discs in a turbulent molecular cloud. Using a piecewise polytropic equation of state, we perform two sets of simulations. In both cases we find that isolated systems undergo a fundamentally different evolution than members of binary or multiple systems. When formed, isolated systems must accrete mass and increase their specific angular momentum, leading to the formation of massive, extended discs, which undergo strong gravitational instabilities and are susceptible to disc fragmentation. Fragments with initial masses of 5.5 M_jup, 7.4 M_jup and 12 M_jup are produced in our simulations. In binaries and small clusters, we observe that due to competition for material from the parent core, members do not accrete significant amounts of high specific angular momentum gas relative to isolated systems. We find that discs in multiple systems are strongly self-gravitating but that they are stable against fragmentation due to disc truncation and mass profile steeping by tides, accretion of high specific angular momentum gas by other members, and angular momentum being redirected into members' orbits. In general, we expect disc fragmentation to be less likely in clusters and to be more a feature of isolated systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cosmological Galaxy Formation Simulations Using SPH

We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$\times10^{11}$ M$_\odot$ to 2$\times10^{12}$ M$_\odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Measuring the Clump Mass Function in the Age of SCUBA2, Herschel, and ALMA

We use simulated images of star-forming regions to explore the effects of various image acquisition techniques on the derived clump mass function. In particular, we focus on the effects of finite image angular resolution, the presence of noise, and spatial filtering. We find that, even when the image has been so heavily degraded with added noise and lowered angular resolution that the clumps it contains clearly no longer correspond to pre-stellar cores, still the clump mass function is typically consistent with the stellar initial mass function within their mutual uncertainties. We explain this result by suggesting that noise, source blending, and spatial filtering all randomly perturb the clump masses, biasing the mass function toward a lognormal form whose high-mass end mimics a Salpeter power law. We argue that this is a consequence of the central limit theorem and that it strongly limits our ability to accurately measure the true mass function of the clumps. We support this conclusion by showing that the characteristic mass scale of the clump mass function, represented by the ``break mass'', scales as a simple function of the angular resolution of the image from which the clump mass function is derived. This strongly constrains our ability to use the clump mass function to derive a star formation efficiency. We discuss the potential and limitations of the current and next generation of instruments for measuring the clump mass function.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Angular Momentum Transport in Extended Galactic Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Forming a Large Disc Galaxy from a z<1 Major Merger

Using high resolution SPH simulations in a fully cosmological Lambda CDM context we study the formation of a bright disk dominated galaxy that originates from a "wet" major merger at z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly following the last major merger and the present day properties of the final remnant are typical of early type spiral galaxies, with an i band B/D ~0.65, a disk scale length of 7.2 kpc, g-r = 0.5 mag, an HI line width (W_{20}/2) of 238 km/sec and total magnitude i = -22.4. The key ingredients for the formation of a dominant stellar disk component after a major merger are: i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, ii) supernova feedback that is able to partially suppress star formation during mergers and iii) relative fading of the spheroidal component. The gas fraction of the progenitors' disks does not exceed 25% at z<3, emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of disks and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disk galaxies by the high likelihood of interactions and mergers for galaxy sized halos at relatively low z.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Role of Cold Flows in the Assembly of Galaxy Disks

We use high resolution cosmological hydrodynamical simulations to demonstrate that cold flow gas accretion, particularly along filaments, modifies the standard picture of gas accretion and cooling onto galaxy disks. In the standard picture, all gas is initially heated to the virial temperature of the galaxy as it enters the virial radius. Low mass galaxies are instead dominated by accretion of gas that stays well below the virial temperature, and even when a hot halo is able to develop in more massive galaxies there exist dense filaments that penetrate inside of the virial radius and deliver cold gas to the central galaxy. For galaxies up to ~L*, this cold accretion gas is responsible for the star formation in the disk at all times to the present. Even for galaxies at higher masses, cold flows dominate the growth of the disk at early times. Within this modified picture, galaxies are able to accrete a large mass of cold gas, with lower initial gas temperatures leading to shorter cooling times to reach the disk. Although star formation in the disk is mitigated by supernovae feedback, the short cooling times allow for the growth of stellar disks at higher redshifts than predicted by the standard model.

0 citations0 reviewsNo code linkedOpen access