Researcher profile

G. Stinson

G. Stinson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2016arXiv

The different baryonic Tully-Fisher relations at low masses

We compare the Baryonic Tully-Fisher relation (BTFR) of simulations and observations of galaxies ranging from dwarfs to spirals, using various measures of rotational velocity Vrot. We explore the BTFR when measuring Vrot at the flat part of the rotation curve, Vflat, at the extent of HI gas, Vlast, and using 20% (W20) and 50% (W50) of the width of HI line profiles. We also compare with the maximum circular velocity of the parent halo, Vmax, within dark matter only simulations. The different BTFRs increasingly diverge as galaxy mass decreases. Using Vlast one obtains a power law over four orders of magnitude in baryonic mass, with slope similar to the observed BTFR. Measuring Vflat gives similar results as Vlast when galaxies with rising rotation curves are excluded. However, higher rotation velocities would be found for low mass galaxies if the cold gas extended far enough for Vrot to reach a maximum. W20 gives a similar slope as Vlast but with slightly lower values of Vrot for low mass galaxies, although this may depend on the extent of the gas in your galaxy sample. W50 bends away from these other relations toward low velocities at low masses. By contrast, Vmax bends toward high velocities for low mass galaxies, as cold gas does not extend out to the radius at which halos reach Vmax. Our study highlights the need for careful comparisons between observations and models: one needs to be consistent about the particular method of measuring Vrot, and precise about the radius at which velocities are measured.

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.

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.

preprint2014arXiv

The Main Sequence and the Fundamental Metallicity Relation in MaGICC Galaxies: Evolution and Scatter

Using cosmological galaxy simulations from the MaGICC project, we study the evolution of the stellar masses, star formation rates and gas phase abundances of star forming galaxies. We derive the stellar masses and star formation rates using observational relations based on spectral energy distributions by applying the new radiative transfer code GRASIL-3D to our simulated galaxies. The simulations match well the evolution of the stellar mass-halo mass relation, have a star forming main sequence that maintains a constant slope out to redshift z $\sim$ 2, and populate projections of the stellar mass - star formation - metallicity plane, similar to observed star forming disc galaxies. We discuss small differences between these projections in observational data and in simulations, and the possible causes for the discrepancies. The light-weighted stellar masses are in good agreement with the simulation values, the differences between the two varying between 0.06 dex and 0.20 dex. We also find a good agreement between the star formation rate tracer and the true (time-averaged) simulation star formation rates. Regardless if we use mass- or light-weighted quantities, our simulations indicate that bursty star formation cycles can account for the scatter in the star forming main sequence.

preprint2012arXiv

A two-phase scenario for bulge assembly in LCDM cosmologies

We analyze and compare the bulges of a sample of L* spiral galaxies in hydrodynamical simulations in a cosmological context, using two different codes, P-DEVA and GASOLINE. The codes regulate star formation in very different ways, with P-DEVA simulations inputing low star formation efficiency under the assumption that feedback occurs on subgrid scales, while the GASOLINE simulations have feedback which drives large scale outflows. In all cases, the marked knee-shape in mass aggregation tracks, corresponding to the transition from an early phase of rapid mass assembly to a later slower one, separates the properties of two populations within the simulated bulges. The bulges analyzed show an important early starburst resulting from the collapse-like fast phase of mass assembly, followed by a second phase with lower star formation, driven by a variety of processes such as disk instabilities and/or mergers. Classifying bulge stellar particles identified at z=0 into old and young according to these two phases, we found bulge stellar sub-populations with distinct kinematics, shapes, stellar ages and metal contents. The young components are more oblate, generally smaller, more rotationally supported, with higher metallicity and less alpha-element enhanced than the old ones. These results are consistent with the current observational status of bulges, and provide an explanation for some apparently paradoxical observations, such as bulge rejuvenation and metal-content gradients observed. Our results suggest that bulges of L* galaxies will generically have two bulge populations which can be likened to classical and pseudo-bulges, with differences being in the relative proportions of the two, which may vary due to galaxy mass and specific mass accretion and merger histories.

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.

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.

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.

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.

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.

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.

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.

preprint2010arXiv

High precision spectra at large redshift for dynamical DE cosmologies

The next generation mass probes will investigate DE nature by measuring non-linear power spectra at various z, and comparing them with high precision simulations. Producing a complete set of them, taking into account baryon physics and for any DE state equation w(z), would really be numerically expensive. Regularities reducing such duty are essential. This paper presents further n-body tests of a relation we found, linking models with DE state parameter w(z) to const.-w models, and also tests the relation in hydro simulations.