Source author record

Brad K. Gibson

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

60works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

60 published item(s)

preprint2022arXiv

Radioactive nuclei in the early Solar System: analysis of the 15 isotopes produced by core-collapse supernovae

Short-lived radioactive isotopes (SLRs) with half-lives between 0.1 to 100 Myr can be used to probe the origin of the Solar System. In this work, we examine the core-collapse supernovae production of the 15 SLRs produced: $^{26}$Al, $^{36}$Cl, $^{41}$Ca, $^{53}$Mn, $^{60}$Fe, $^{92}$Nb, $^{97}$Tc, $^{98}$Tc, $^{107}$Pd, $^{126}$Sn, $^{129}$I, $^{135}$Cs, $^{146}$Sm, $^{182}$Hf, and $^{205}$Pb. We probe the impact of the uncertainties of the core-collapse explosion mechanism by examining a collection of 62 core-collapse models with initial masses of 15, 20, and 25M$_{\odot}$, explosion energies between 3.4$\times$10$^{50}$ and 1.8$\times$10$^{52}$ ergs and compact remnant masses between 1.5M$_{\odot}$and 4.89M$_{\odot}$. We identify the impact of both explosion energy and remnant mass on the final yields of the SLRs. Isotopes produced within the innermost regions of the star, such as $^{92}$Nb and $^{97}$Tc, are the most affected by the remnant mass, $^{92}$Nb varying by five orders of magnitude. Isotopes synthesised primarily in explosive C-burning and explosive He-burning, such as $^{60}$Fe, are most affected by explosion energies. $^{60}$Fe increases by two orders of magnitude from the lowest to the highest explosion energy in the 15M$_{\odot}$model. The final yield of each examined SLR is used to compare to literature models.

preprint2021arXiv

A Large-scale Approach to Modelling Molecular Biosignatures: The Diatomics

This work presents the first steps to modelling synthetic rovibrational spectra for all molecules of astrophysical interest using a new approach implemented in the Prometheus code. The goal is to create a new comprehensive source of first-principles molecular spectra, thus bridging the gap for missing data to help drive future high-resolution studies. Our primary application domain is for molecules identified as signatures of life in planetary atmospheres (biosignatures), but our approach is general and can be applied to other systems. In this work we evaluate the accuracy of our method by studying four diatomic molecules H$_2$, O$_2$, N$_2$ and CO, all of which have well-known spectra. Prometheus uses the Transition-Optimised Shifted Hermite (TOSH) theory to account for anharmonicity for the fundamental $ν=0 \rightarrow ν=1$ band, along with thermal profile modeling for the rotational transitions. To this end, we expand TOSH theory to enable the modeling of rotational constants. We show that our simple model achieves results that are a better approximation of the real spectra than those produced through a harmonic approach. We compare our results with high-resolution HITRAN and ExoMol spectral data. We find that modelling accuracy tends to diminish for rovibrational transition away from the band origin, thus highlighting the need for the theory to be further adapted.

preprint2020arXiv

Explaining the chemical trajectories of accreted and in-situ halo stars of the Milky Way

The Milky Way underwent its last significant merger ten billion years ago, when the Gaia-Enceladus-Sausage (GES) was accreted. Accreted GES stars and progenitor stars born prior to the merger make up the bulk of the inner halo. Even though these two main populations of halo stars have similar $durations$ of star formation prior to their merger, they differ in [$α$/Fe]-[Fe/H] space, with the GES population bending to lower [$α$/Fe] at a relatively low value of [Fe/H]. We use cosmological simulations of a 'Milky Way' to argue that the different tracks of the halo stars through the [$α$/Fe]-[Fe/H] plane are due to a difference in their star formation history and efficiency, with the lower mass GES having its low and constant star formation regulated by feedback whilst the higher mass main progenitor has a higher star formation rate prior to the merger. The lower star formation efficiency of GES leads to lower gas pollution levels, pushing [$α$/Fe]-[Fe/H] tracks to the left. In addition, the increasing star formation rate maintains a higher relative contribution of Type~II SNe to Type~Ia SNe for the main progenitor population that formed during the same time period, thus maintaining a relatively high [$α$/Fe]. Thus the different positions of the downturns in the [$α$/Fe]-[Fe/H] plane for the GES stars are not reflective of different star formation durations, but instead reflect different star formation efficiencies. We argue that cosmological simulations match a wide range of independent observations, breaking degeneracies that exist in simpler models.

preprint2020arXiv

The $R$-Process Alliance: Fourth Data Release from the Search for $r$-Process-Enhanced Stars in the Galactic Halo

This compilation is the fourth data release from the $R$-Process Alliance (RPA) search for $r$-process-enhanced stars, and the second release based on "snapshot" high-resolution ($R \sim 30,000$) spectra collected with the du Pont 2.5m Telescope. In this data release, we propose a new delineation between the $r$-I and $r$-II stellar classes at $\mathrm{[Eu/Fe]} = +0.7$, instead of the empirically chosen $\mathrm{[Eu/Fe]} = +1.0$ level previously in use, based on statistical tests of the complete set of RPA data released to date. We also statistically justify the minimum level of [Eu/Fe] for definition of the $r$-I stars, [Eu/Fe] $> +0.3$. Redefining the separation between $r$-I and $r$-II stars will aid in analysis of the possible progenitors of these two classes of stars and whether these signatures arise from separate astrophysical sources at all. Applying this redefinition to previous RPA data, the number of identified $r$-II and $r$-I stars changes to 51 and 121, respectively, from the initial set of data releases published thus far. In this data release, we identify 21 new $r$-II, 111 new $r$-I (plus three re-identified), and 7 new (plus one re-identified) limited-$r$ stars out of a total of 232 target stars, resulting in a total sample of 72 new $r$-II stars, 232 new $r$-I stars, and 42 new limited-$r$ stars identified by the RPA to date.

preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- I: Survey Description, Spectra and Radial Velocities

The Radial Velocity Experiment (RAVE) is a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A). The 6th and final data release (DR6 or FDR) is based on 518387 observations of 451783 unique stars. RAVE observations were taken between 12 April 2003 and 4 April 2013. Here we present the genesis, setup and data reduction of RAVE as well as wavelength-calibrated and flux-normalized spectra and error spectra for all observations in RAVE DR6. Furthermore, we present derived spectral classification and radial velocities for the RAVE targets, complemented by cross matches with Gaia DR2 and other relevant catalogs. A comparison between internal error estimates, variances derived from stars with more than one observing epoch and a comparison with radial velocities of Gaia DR2 reveals consistently that 68% of the objects have a velocity accuracy better than 1.4 km/s, while 95% of the objects have radial velocities better than 4.0 km/s. Stellar atmospheric parameters, abundances and distances are presented in subsequent publication. The data can be accessed via the RAVE Web (http://rave-survey.org) or the Vizier database.

preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- II: Stellar Atmospheric Parameters, Chemical Abundances and Distances

We present part 2 of the 6th and final Data Release (DR6 or FDR) of the Radial Velocity Experiment (RAVE), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A) and span the complete time frame from the start of RAVE observations on 12 April 2003 to their completion on 4 April 2013. In the second of two publications, we present the data products derived from 518387 observations of 451783 unique stars using a suite of advanced reduction pipelines focussing on stellar atmospheric parameters, in particular purely spectroscopically derived stellar atmospheric parameters (Teff, log(g), and the overall metallicity), enhanced stellar atmospheric parameters inferred via a Bayesian pipeline using Gaia DR2 astrometric priors, and asteroseismically calibrated stellar atmospheric parameters for giant stars based on asteroseismic observations for 699 K2 stars. In addition, we provide abundances of the elements Fe, Al, and Ni, as well as an overall [alpha/Fe] ratio obtained using a new pipeline based on the GAUGUIN optimization method that is able to deal with variable signal-to-noise ratios. The RAVE DR6 catalogs are cross matched with relevant astrometric and photometric catalogs, and are complemented by orbital parameters and effective temperatures based on the infrared flux method. The data can be accessed via the RAVE Web site (http://rave-survey.org) or the Vizier database.

preprint2016arXiv

Chemical separation of disc components using RAVE

We present evidence from the RAdial Velocity Experiment (RAVE) survey of chemically separated, kinematically distinct disc components in the solar neighbourhood. We apply probabilistic chemical selection criteria to separate our sample into $α$-low (`thin disc') and $α$-high (`thick disc') sequences. Using newly derived distances, which will be utilized in the upcoming RAVE DR5, we explore the kinematic trends as a function of metallicity for each of the disc components. For our $α$-low disc, we find a negative trend in the mean rotational velocity ($V_{\mathrmϕ}$) as a function of iron abundance ([Fe/H]). We measure a positive gradient $\partial V_{\mathrmϕ}$/$\partial$[Fe/H] for the $α$-high disc, consistent with results from high-resolution surveys. We also find differences between the $α$-low and $α$-high discs in all three components of velocity dispersion. We discuss the implications of an $α$-low, metal-rich population originating from the inner Galaxy, where the orbits of these stars have been significantly altered by radial mixing mechanisms in order to bring them into the solar neighbourhood. The probabilistic separation we propose can be extended to other data sets for which the accuracy in [$α$/Fe] is not sufficient to disentangle the chemical disc components a priori. For such datasets which will also have significant overlap with Gaia DR1, we can therefore make full use of the improved parallax and proper motion data as it becomes available to investigate kinematic trends in these chemical disc components.

preprint2016arXiv

The role of gas infall in the evolution of disc galaxies

Spiral galaxies are thought to acquire their gas through a protracted infall phase resulting in the inside-out growth of their associated discs. For field spirals, this infall occurs in the lower density environments of the cosmic web. The overall infall rate, as well as the galactocentric radius at which this infall is incorporated into the star-forming disc, plays a pivotal role in shaping the characteristics observed today. Indeed, characterising the functional form of this spatio-temporal infall in-situ is exceedingly difficult, and one is forced to constrain these forms using the present day state of galaxies with model or simulation predictions. We present the infall rates used as input to a grid of chemical evolution models spanning the mass spectrum of discs observed today. We provide a systematic comparison with alternate analytical infall schemes in the literature, including a first comparison with cosmological simulations. Identifying the degeneracies associated with the adopted infall rate prescriptions in galaxy models is an important step in the development of a consistent picture of disc galaxy formation and evolution.

preprint2015arXiv

Galactic Chemical Evolution: Stellar Yields and the Initial Mass Function

We present a set of 144 galactic chemical evolution models applied to a Milky Way analogue, computed using four sets of low and intermediate star nucleosynthetic yields, six massive star yield compilations, and six functional forms for the initial mass function. The integrated or true yields for each combination are derived. A comparison is made between a grid of multiphase chemical evolution models computed with these yield combinations and empirical data drawn from the Milky Way's disc, including the solar neighbourhood. By means of a chi2 methodology, applied to the results of these multiphase models, the best combination of stellar yields and initial mass function capable of reproducing these observations is identified.

preprint2015arXiv

New quasars behind the Magellanic Clouds. Spectroscopic confirmation of near-infrared selected candidates

Quasi--stellar objects (quasars) located behind nearby galaxies provide an excellent absolute reference system for astrometric studies, but they are difficult to identify because of fore- and background contamination. Deep wide--field, high angular resolution surveys spanning the entire area of nearby galaxies are needed to obtain a complete census of such quasars. We embarked on a program to expand the quasar reference system behind the Large and the Small Magellanic Clouds, the Magellanic Bridge, and the Magellanic Stream, connecting the Clouds with the Milky Way. Hundreds of quasar candidates were selected based on their near--infrared colors and variability properties from the ongoing public ESO VISTA Magellanic Clouds survey. A subset of 49 objects was followed up with optical spectroscopy. We confirmed the quasar nature of 37 objects (34 new identifications), four are low redshift objects, three are probably stars, and the remaining three lack prominent spectral features for a secure classification; bona fide quasars, judging from their broad absorption lines are located, as follows: 10 behind the LMC, 13 behind the SMC, and 14 behind the Bridge. The quasars span a redshift range from z~0.5 to z~4.1. Upon completion the VMC survey is expected to yield a total of ~1500 quasars with Y<19.32 mag, J<19.09 mag, and Ks<18.04 mag.

preprint2015arXiv

Origin of the Metallicity Distribution in the Thick Disc

Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends. Methods. We analyse five simulated spiral galaxies with masses comparable to the Milky Way. The simulations span a range of star formation and energy feedback strengths and prescriptions, particle- and grid-based hydrodynamical implementations, as well as initial conditions/assembly history. Results. Consistently, we find that the steeper, negative, radial metallicity gradients seen in the mid-plane flatten with increasing height away from the plane. In simulations with stronger (and/or more spatially-extended) feedback, the negative radial gradients invert, becoming positive for heights in excess of 1 kpc. Such behaviour is consistent with that inferred from recent observations. Our measurements of the vertical metallicity gradients show no clear correlation with galactocentric radius, and are in good agreement with those observed in the Milky Way's thick disc (locally). Conclusions. Simulations employing stronger/more extended feedback prescriptions possess radial and vertical metallicity and kinematic gradients more in line with recent observations. The inverted, positive, radial metallicity gradients seen in the simulated thick stellar discs originate from a population of younger, more metal-rich, stars formed in-situ, superimposed upon a background population of older migrators from the inner disc; the contrast provided by the former increases radially, due to the inside-out growth of the disc. A similar behaviour may be responsible for the same flattening seen in the radial gradients with scaleheight in the Milky Way.

preprint2015arXiv

Simple and accurate modelling of the gravitational potential produced by thick and thin exponential disks

We present accurate models of the gravitational potential produced by a radially exponential disk mass distribution. The models are produced by combining three separate Miyamoto-Nagai disks. Such models have been used previously to model the disk of the Milky Way, but here we extend this framework to allow its application to disks of any mass, scalelength, and a wide range of thickness from infinitely thin to near spherical (ellipticities from 0 to 0.9). The models have the advantage of simplicity of implementation, and we expect faster run speeds over a double exponential disk treatment. The potentials are fully analytical, and differentiable at all points. The mass distribution of our models deviates from the radial mass distribution of a pure exponential disk by <0.4% out to 4 disk scalelengths, and <1.9% out to 10 disk scalelengths. We tabulate fitting parameters which facilitate construction of exponential disks for any scalelength, and a wide range of disk thickness (a user-friendly, web-based interface is also available). Our recipe is well suited for numerical modelling of the tidal effects of a giant disk galaxy on star clusters or dwarf galaxies. We consider three worked examples; the Milky Way thin and thick disk, and a disky dwarf galaxy.

preprint2015arXiv

The Sensitivity of Harassment to Orbit: Mass Loss from Early-Type Dwarfs in Galaxy Clusters

We conduct a comprehensive numerical study of the orbital dependence of harassment on early-type dwarfs consisting of 168 different orbits within a realistic, Virgo-like cluster, varying in eccentricity and pericentre distance. We find harassment is only effective at stripping stars or truncating their stellar disks for orbits that enter deep into the cluster core. Comparing to the orbital distribution in cosmological simulations, we find that the majority of the orbits (more than three quarters) result in no stellar mass loss. We also study the effects on the radial profiles of the globular cluster systems of early-type dwarfs. We find these are significantly altered only if harassment is very strong. This suggests that perhaps most early-type dwarfs in clusters such as Virgo have not suffered any tidal stripping of stars or globular clusters due to harassment, as these components are safely embedded deep within their dark matter halo. We demonstrate that this result is actually consistent with an earlier study of harassment of dwarf galaxies, despite the apparent contradiction. Those few dwarf models that do suffer stellar stripping are found out to the virial radius of the cluster at redshift=0, which mixes them in with less strongly harassed galaxies. However when placed on phase-space diagrams, strongly harassed galaxies are found offset to lower velocities compared to weakly harassed galaxies. This remains true in a cosmological simulation, even when halos have a wide range of masses and concentrations. Thus phase-space diagrams may be a useful tool for determining the relative likelihood that galaxies have been strongly or weakly harassed.

preprint2015arXiv

Tracking Cluster Debris (TraCD)- I. Dissolution of clusters and searching the solar cradle

The capability to reconstruct dissolved stellar systems in dynamical and chemical space is a key factor in improving our understanding of the evolution of the Milky Way. Here we concentrate on the dynamical aspect and given that a significant portion of the stars in the Milky Way have been born in stellar associations or clusters that have lived a few Myr up to several Gyr, we further restrict our attention to the evolution of star clusters. We have carried out our simulations in two steps: (1) we create a simulation of dissolution and mixing processes which yields a close fit to the present-day Milky Way dynamics and (2) we have evolved three sets of stellar clusters with masses of 400, 1000 and 15000 M$_{\odot}$ to dissolution. The birth location of these sets was 4, 6, 8 and 10 kpc for the 400 and 1000 M$_{\odot}$ clusters and 4, 6, 8, 10 and 12 kpc for the 15000 M$_{\odot}$. We have focused our efforts on studying the state of the escapers from these clusters after 4.5 Gyr of evolution with particular attention to stars that reach the Solar annulus, i.e. 7.5 $\le$ R$_{\rm{gc}}$ $\le$ 8.5 kpc. We give results for Solar twins and siblings over a wide range of radii and cluster masses for two dissolution mechanisms. From kinematics alone, we conclude that the Sun was $\sim$50 per cent more likely to have been born near its current Galactocentric radius, rather than have migrated (radially) $\sim$2 kpc since birth. We conclude our analysis by calculating magnitudes and colours of our single stars for comparison with the samples that the $Gaia$, $Gaia$-ESO and GALAH-AAO surveys will obtain. In terms of reconstructing dissolved star clusters we find that on short time-scales we cannot rely on kinematic evolution alone and thus it will be necessary to extend our study to include information on chemical space.

preprint2014arXiv

Cloud angular momentum and effective viscosity in global SPH simulations with feedback

We examine simulations of isolated galaxies to analyse the effects of localised feedback on the formation and evolution of molecular clouds. Feedback contributes to turbulence and the destruction of clouds, leading to a population of clouds that is younger, less massive, and with more retrograde rotation. We investigate the evolution of clouds as they interact with each other and the diffuse ISM, and determine that the role of cloud interactions differs strongly with the presence of feedback: in models without feedback, scattering events dramatically increase the retrograde fraction, but in models with feedback, mergers between clouds may slightly increase the prograde fraction. We also produce an estimate of the viscous time-scale due to cloud-cloud collisions, which increases with increasing strength of feedback (~20 Gyr vs ~10 Gyr), but is still much smaller than previous estimates (~1000 Gyr); although collisions become more frequent with feedback, less energy is lost in each collision than in the models without feedback.

preprint2014arXiv

Pseudo Three-dimensional Maps of the Diffuse Interstellar Band at 862 nm

The diffuse interstellar bands (DIBs) are absorption lines observed in visual and near infrared spectra of stars. Understanding their origin in the interstellar medium (ISM) is one of the oldest problems in astronomical spectroscopy, as DIBs ahave been known since 1922. In a completely new approach to understanding DIBs, we combined information from nearly 500,000 stellar spectra obtained by the massive spectroscopic survey RAVE (Radial Velocity Experiment) to produce the first pseudo three-dimensional map of the strength of the DIB at 8620~Ångstroms covering the nearest 3 kiloparsecs from the Sun, and show that it follows our independently constructed spatial distribution of extinction by interstellar dust along the Galactic plane. Despite having a similar distribution in the Galactic plane, the DIB 8620 carrier has a significantly larger vertical scale height than the dust. Even if one DIB may not represent the general DIB population, our observations outline the future direction of DIB research.

preprint2014arXiv

Stellar origin of the 182Hf cosmochronometer and the presolar history of solar system matter

Among the short-lived radioactive nuclei inferred to be present in the early solar system via meteoritic analyses, there are several heavier than iron whose stellar origin has been poorly understood. In particular, the abundances inferred for 182Hf (half-life = 8.9 million years) and 129I (half-life = 15.7 million years) are in disagreement with each other if both nuclei are produced by the rapid neutron-capture process. Here, we demonstrate that contrary to previous assumption, the slow neutron-capture process in asymptotic giant branch stars produces 182Hf. This has allowed us to date the last rapid and slow neutron-capture events that contaminated the solar system material at roughly 100 million years and 30 million years, respectively, before the formation of the Sun.

preprint2013arXiv

Decoding the message from meteoritic stardust silicon carbide grains

Micron-sized stardust grains that originated in ancient stars are recovered from meteorites and analysed using high-resolution mass spectrometry. The most widely studied type of stardust is silicon carbide (SiC). Thousands of these grains have been analysed with high precision for their Si isotopic composition. Here we show that the distribution of the Si isotopic composition of the vast majority of stardust SiC grains carry the imprints of a spread in the age-metallicity distribution of their parent stars and of a power-law increase of the relative formation efficiency of SiC dust with the metallicity. This result offers a solution for the long-standing problem of silicon in stardust SiC grains, confirms the necessity of coupling chemistry and dynamics in simulations of the chemical evolution of our Galaxy, and constrains the modelling of dust condensation in stellar winds as function of the metallicity.

preprint2013arXiv

Numerical Simulations of Bubble-Induced Star Formation in Dwarf Irregular Galaxies with a Novel Stellar Feedback Scheme

To study the star formation and feedback mechanism, we simulate the evolution of an isolated dwarf irregular galaxy (dIrr) in a fixed dark matter halo, similar in size to WLM, using a new stellar feedback scheme. We use the new version of our original N-body/smoothed particle chemodynamics code, GCD+, which adopts improved hydrodynamics, metal diffusion between the gas particles and new modelling of star formation and stellar wind and supernovae (SNe) feedback. Comparing the simulations with and without stellar feedback effects, we demonstrate that the collisions of bubbles produced by strong feedback can induce star formation in a more widely spread area. We also demonstrate that the metallicity in star forming regions is kept low due to the mixing of the metal-rich bubbles and the metal-poor inter-stellar medium. Our simulations also suggest that the bubble-induced star formation leads to many counter-rotating stars. The bubble-induced star formation could be a dominant mechanism to maintain star formation in dIrrs, which is different from larger spiral galaxies where the non-axisymmetric structures, such as spiral arms, are a main driver of star formation.

preprint2012arXiv

Calibrating an updated SPH scheme within GCD+

We adapt a modern scheme of smoothed particle hydrodynamics (SPH) to our tree N-body/SPH galactic chemodynamics code GCD+. The applied scheme includes imple- mentations of the artificial viscosity switch and artificial thermal conductivity pro- posed by Morris & Monaghan (1997), Rosswog & Price (2007) and Price (2008), to model discontinuities and Kelvin-Helmholtz instabilities more accurately. We first present hydrodynamics test simulations and contrast the results to runs undertaken without artificial viscosity switch or thermal conduction. In addition, we also explore the different levels of smoothing by adopting larger or smaller smoothing lengths, i.e. a larger or smaller number of neighbour particles, Nnb. We demonstrate that the new version of GCD+ is capable of modelling Kelvin-Helmholtz instabilities to a simi- lar level as the mesh code, Athena. From the Gresho vortex, point-like explosion and self-similar collapse tests, we conclude that setting the smoothing length to keep the number of neighbour particles as high as Nnb~58 is preferable to adopting smaller smoothing lengths. We present our optimised parameter sets from the hydrodynamics tests.

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.

preprint2012arXiv

RAMSES-CH: A New Chemodynamical Code for Cosmological Simulations

We present a new chemodynamical code - Ramses-CH - for use in simulating the self-consistent evolution of chemical and hydrodynamical properties of galaxies within a fully cosmological framework. We build upon the adaptive mesh refinement code Ramses, which includes a treatment of self-gravity, hydrodynamics, star formation, radiative cooling, and supernovae feedback, to trace the dominant isotopes of C, N, O, Ne, Mg, Si, and Fe. We include the contribution of Type Ia and II supernovae, in addition to low- and intermediate-mass asymptotic giant branch stars, relaxing the instantaneous recycling approximation. The new chemical evolution modules are highly flexible and portable, lending themselves to ready exploration of variations in the underpining stellar and nuclear physics. We apply Ramses-CH to the cosmological simulation of a typical L* galaxy, demonstrating the successful recovery of the basic empirical constraints regarding, [α/Fe]-[Fe/H] and Type Ia/II supernovae rates.

preprint2012arXiv

RAMSES-CH: A New Chemodynamics Code

We present a new chemodynamical code based on the adaptive mesh refinement code RAMSES. The new code uses Eulerian hydrodynamics and N-body dynamics in a cosmological framework to trace the production and advection of several chemical species. It is the first such code to follow the self-consistent evolution of chemical elements in cosmological volumes while maintaining sub-kiloparsec resolution. The code will be used to simulate disk galaxies and explore the influence of chemical evolution models and star formation on galactic abundance ratios.

preprint2011arXiv

A Comparison of Galaxy Group Luminosity Functions from Semi-Analytic Models

Semi-analytic models (SAMs) are currently one of the primary tools with which to model statistically significant ensembles of galaxies. The underlying physical prescriptions inherent to each SAM are, in many cases, different from one another. Several SAMs have been applied to the dark matter merger trees extracted from the Millennium Run, including those associated with the well-known Munich and Durham lineages. We compare the predicted luminosity distributions of galaxy groups using four publicly available SAMs (De Lucia et al. 2006; Bower et al. 2006; Bertone et al. 2007; Font et al. 2008), in order to explore a galactic environment in which the models have not been explored to the same degree as they have in the field or in rich clusters. We identify a characteristic "wiggle" in the group galaxy luminosity function generated using the De Lucia et al. (2006) SAM, that is not present in the Durham-based models, consistent to some degree with observations. However, a comparison between conditional luminosity functions of groups between the models and observations of Yang et al. (2007) suggest that neither model is a particularly good match. The luminosity function wiggle is interpreted as the result of the two-mode AGN feedback implementation used in the Munich models, which itself results in flattened magnitude gap distribution. An associated analysis of the magnitude gap distribution between first- and second-ranked group galaxies shows that while the Durham models yield distributions with approximately equal luminosity first- and second-ranked galaxies, in agreement with observations, the De Lucia et al. models favours the scenario in which the second-ranked galaxy is approximately one magnitude fainter than the primary,especially when the dynamic range of the mock data is limited to 3 magnitudes.

preprint2011arXiv

Chemical signature of gas-rich disc-disc mergers at high redshift

We performed numerical simulations of mergers between gas-rich disc galaxies, which result in the formation of late-type galaxies. Stars formed during the merger end up in a thick disc that is partially supported by velocity dispersion and has high [alpha/Fe] ratios at all metallicities. Stars formed later end up in a thin, rotationally supported disc which has lower [alpha/Fe] ratios. While the structural and kinematical properties of the merger remnants depend strongly upon the orbital parameters of the mergers, we find a clear chemical signature of gas-rich mergers.

preprint2011arXiv

Disc Heating: Comparing the Milky Way with Cosmological Simulations

We present the analysis of a suite of simulations run with different particle-and grid-based cosmological hydrodynamical codes and compare them with observational data of the Milky Way. This is the first study to make comparisons of properties of galaxies simulated with particle and grid-based codes. Our analysis indicates that there is broad agreement between these different modelling techniques. We study the velocity dispersion - age relation for disc stars at z=0 and find that four of the simulations are more consistent with observations by Holmberg et al. (2008) in which the stellar disc appears to undergo continual/secular heating. Two other simulations are in better agreement with the Quillen & Garnett (2001) observations that suggest a "saturation" in the heating profile for young stars in the disc. None of the simulations have thin discs as old as that of the Milky Way. We also analyse the kinematics of disc stars at the time of their birth for different epochs in the galaxies' evolution and find that in some simulations old stars are born cold within the disc and are subsequently heated, while other simulations possess old stellar populations which are born relatively hot. The models which are in better agreement with observations of the Milky Way's stellar disc undergo significantly lower minor-merger/assembly activity after the last major merger - i.e. once the disc has formed. All of the simulations are significantly "hotter" than the Milky Way disc; on top of the effects of mergers, we find a "floor" in the dispersion that is related to the underlying treatment of the heating and cooling of the interstellar medium, and the low density threshold which such codes use for star formation. This finding has important implications for all studies of disc heating that use hydrodynamical codes.

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.

preprint2011arXiv

Metal-Poor Lithium-Rich Giants in the Radial Velocity Experiment Survey

We report the discovery of eight lithium-rich field giants found in a high resolution spectroscopic sample of over 700 metal-poor stars ([Fe/H]<-0.5) selected from the RAVE survey. The majority of the Li-rich giants in our sample are very metal-poor ([Fe/H]<-1.9), and have a Li abundance (in the form of 7Li), A(Li)=log(n(Li)/n(H))+12, between 2.30 and 3.63, well above the typical upper red giant branch limit, A(Li)<0.5, while two stars, with A(Li)~1.7-1.8, show similar lithium abundances to normal giants at the same gravity. We further included two metal-poor, Li-rich globular cluster giants in our sample, namely the previously discovered M3-IV101 and newly discovered (in this work) M68-A96. This comprises the largest sample of metal-poor Li-rich giants to date. We performed a detailed abundance analysis of all stars, finding that the majority our sample stars have elemental abundances similar to that of Li-normal halo giants. Although the evolutionary phase of each Li-rich giant cannot be definitively determined, the Li-rich phase is likely connected to extra mixing at the red giant branch bump or early asymptotic giant branch that triggers cool bottom processing in which the bottom of the outer convective envelope is connected to the H-burning shell in the star. The surface of a star becomes Li-enhanced as 7Be (which burns to 7Li) is transported to the stellar surface via the Cameron-Fowler mechanism. We discuss and discriminate among several models for the extra mixing that can cause Li-production, given the detailed abundances of the Li-rich giants in our sample.

preprint2011arXiv

Metallicity gradients of disc stars for a cosmologically simulated galaxy

We analyse for the first time the radial abundance gradients of the disc stars of a disc galaxy simulated with our three dimensional, fully cosmological chemodynamical galaxy evolution code GCD+. We study how [Fe/H], [N/O], [O/Fe], [Mg/Fe] and [Si/Fe] vary with galactocentric radius. For the young stars of the disc, we found a negative slope for [Fe/H] and [N/O] but a positive [O/Fe], [Mg/Fe] and [Si/Fe] slope with radius. By analysing the star formation rate (SFR) at different radii, we found that the simulated disc contains a greater fraction of young stars in the outer regions, while the old stars tend to be concentrated in the inner parts of the disc. This can explain the positive [alpha/Fe] gradient as well as the negative [N/O] gradient with radius. This radial trend is a natural outcome of an inside-out formation of the disc, regardless of its size and can thus explain the recently observed positive [alpha/Fe] gradients in the Milky Way disc open clusters.

preprint2011arXiv

Observational Properties of the Metal-Poor Thick Disk of the Milky Way Galaxy and Insights into Its Origins

We have undertaken the study of the elemental abundances and kinematic properties of a metal-poor sample of candidate thick-disk stars selected from the RAVE spectroscopic survey of bright stars to differentiate among the present scenarios of the formation of the thick disk. In this paper, we report on a sample of 214 red giant branch, 31 red clump/horizontal branch, and 74 main-sequence/sub-giant branch metal-poor stars, which serves to augment our previous sample of only giant stars. We find that the thick disk [alpha/Fe] ratios are enhanced, and have little variation (<0.1 dex), in agreement with our previous study. The augmented sample further allows, for the first time, investigation of the gradients in the metal-poor thick disk. For stars with [Fe/H] < -1.2, the thick disk shows very small gradients, <0.03 +/- 0.02 dex/kpc, in alpha-enhancement, while we find a +0.01 +/- 0.04 dex/kpc radial gradient and a -0.09 +/- 0.05 dex/kpc vertical gradient in iron abundance. In addition, we show that the peak of the distribution of orbital eccentricities for our sample agrees better with models in which the stars that comprise the thick disk were formed primarily in the Galaxy, with direct accretion of stars contributing little. Our results thus disfavor direct accretion of stars from dwarf galaxies into the thick disk as a major contributor to the thick disk population, but cannot discriminate between alternative models for the thick disk, such as those that invoke high-redshift (gas-rich) mergers, heating of a pre-existing thin stellar disk by a minor merger, or efficient radial migration of stars.

preprint2011arXiv

Star formation history of barred disc galaxies

We present the first results of a pilot study aimed at understanding the influence of bars on the evolution of galaxy discs through the study of their stellar content. We examine here the kinematics, star formation history, mass-weighted, luminosity-weighted, and single stellar population (SSP) equivalent ages and metallicities for four galaxies ranging from lenticulars to late-type spirals. The data employed extends to 2-3 disc scalelengths, with S/N(A)>50. Several techniques are explored to derive star formation histories and SSP-equivalent parameters, each of which are shown to be compatible. We demostrate that the age-metallicity degeneracy is highly reduced by using spectral fitting techniques --instead of indices-- to derive these parameters. We found that the majority of the stellar mass in our sample is composed of old (~10 Gyr) stars. This is true in the bulge and the disc region, even beyond two disc scalelengths. In the bulge region, we find that the young, dynamically cold, structures produced by the presence of the bar (e.g., nuclear discs or rings) are responsible for shaping the bulges' age and metallicity gradients. In the disc region, a larger fraction of young stars is present in the external parts of the disc compared with the inner disc. The disc growth is, therefore, compatible with a moderate inside-out formation scenario, where the luminosity weighted age changes from ~10 Gyrs in the centre, to ~4 Gyrs at two disc scalelengths, depending upon the galaxy. For two galaxies, we compare the metallicity and age gradients of the disc major axis with that of the bar, finding very important differences. In particular, the stellar population of the bar is more similar to the bulge than to the disc, indicating that, at least in those two galaxies, bars formed long ago and have survived to the present day. (abridged)

preprint2011arXiv

The Chemical Evolution of Globular Clusters - II. Metals and Fluorine

In the first paper in this series, we proposed a new framework in which to model the chemical evolution of globular clusters. This model, is predicated upon the assumption that clusters form within an interstellar medium enriched locally by the ejecta of a single Type Ia supernova and varying numbers of asymptotic giant branch stars, superimposed on an ambient medium pre-enriched by low-metallicity Type II supernovae. Paper I was concerned with the application of this model to the observed abundances of several reactive elements and so-called non-metals for three classical intermediate-metallicity clusters, with the hallmark of the work being the successful recovery of many of their well-known elemental and isotopic abundance anomalies. Here, we expand upon our initial analysis by (a) applying the model to a much broader range of metallicities (from the factor of three explored in Paper I, to now, a factor of ~50; i.e., essentially, the full range of Galactic globular cluster abundances, and (b) incorporating a broader suite of chemical species, including a number of iron-peak isotopes, heavier alpha-elements, and fluorine. While most empirical globular cluster abundance trends are reproduced, our model would suggest the need for a higher production of Ca, Si, and Cu in low-metallicity (or so-called "prompt") Type Ia supernovae than predicted in current stellar models in order to reproduce the observed trends in NGC 6752, and a factor of two reduction in carbon production from asymptotic giant branch stars to explain the observed trends between carbon and nitrogen. Observations of heavy-element isotopes produced primarily by Type Ia supernovae, including those of titanium, iron, and nickel, could support/refute unequivocally our proposed framework. Hydrodynamical simulations would be necessary to study its viability from a dynamical point of view.

preprint2010arXiv

The circum-galactic gas around cosmologically simulated disks

We analyze the physical properties and infall rates of the circum-galactic gas around disks obtained in multi-resolved, cosmological, AMR simulations. At intermediate and low redshifts, disks are embedded into an extended, hot, tenuous corona that contributes largely in fueling the disk with non-enriched gas whereas the accretion of enriched gas from tidal streams occurs throughout episodic events. We derive an infall rate close to the disk of the same value as the one of the star formation rate in the disk and its temporal evolution as a function of galacto-centric radius nicely shows that the growth of galactic disks proceeds according to an inside-out formation scenario.

preprint2010arXiv

The Dawning of the Stream of Aquarius in RAVE

We identify a new, nearby (0.5 < d < 10 kpc) stream in data from the RAdial Velocity Experiment (RAVE). As the majority of stars in the stream lie in the constellation of Aquarius we name it the Aquarius Stream. We identify 15 members of the stream lying between 30 < l < 75 and -70< b <-50, with heliocentric line-of-sight velocities V_los~-200 km/s. The members are outliers in the radial velocity distribution, and the overdensity is statistically significant when compared to mock samples created with both the Besançon Galaxy model and newly-developed code Galaxia. The metallicity distribution function and isochrone fit in the log g - T_eff plane suggest the stream consists of a 10 Gyr old population with [m/H]~-1.0. We explore relations to other streams and substructures, finding the stream cannot be identified with known structures: it is a new, nearby substructure in the Galaxy's halo. Using a simple dynamical model of a dissolving satellite galaxy we account for the localization of the stream. We find that the stream is dynamically young and therefore likely the debris of a recently disrupted dwarf galaxy or globular cluster. The Aquarius stream is thus a specimen of ongoing hierarchical Galaxy formation, rare for being right in the solar suburb.

preprint2009arXiv

Are Dry Mergers Dry, Moist, or Wet?

We present a spectral analysis of a sample of red-sequence galaxies identified by van Dokkum (2005) as dry merger remnants and ongoing dry mergers. Kinematics, stellar population absorption features, and ionisation from emission lines, are derived. We find that approximately half of the sample showing strong tidal features have younger stellar populations than a control sample at a given velocity dispersion. Conversely, galaxies with weak tidal tails and/or ongoing mergers -- with the exception of one galaxy -- do not show this young component. This seems to indicate that the young stellar populations observed in a significant fraction of ellipticals is the consequence of star formation triggered by mergers. This young component is consistent with a light "frosting" of young stars (<2% by mass) superimposed upon a dominant, old (~11 Gyr), stellar population. In terms of stellar populations, these mergers are, in fact, fairly dry. We found, however, that merger remnants with young stellar populations are supported by rotation, contrary to the expectations of a major dry merger. This suggests that the small amount of gas involved has been sufficient to produce a dynamically cold stellar component. Half of the galaxies with strong tidal distortion, however, are slow rotating and have stellar populations compatible with the control sample at a given velocity dispersion. Remarkably, none of the galaxies with velocity dispersions in excess of 250 km/s have a young stellar component, independent of the merger stage.

preprint2009arXiv

Chemodynamical analysis of bulge stars for simulated disc galaxies

We analyse the kinematics and chemistry of the bulge stars of two simulated disc galaxies using our chemodynamical galaxy evolution code GCD+. First we compare stars that are born inside the galaxy with those that are born outside the galaxy and are accreted into the centre of the galaxy. Stars that originate outside of the bulge are accreted into it early in its formation within 3 Gyrs so that these stars have high [alpha/Fe] as well as having a high total energy reflecting their accretion to the centre of the galaxy. Therefore, higher total energy is a good indicator for finding accreted stars. The bulges of the simulated galaxies formed through multiple mergers separated by about a Gyr. Since [alpha/Fe] is sensitive to the first few Gyrs of star formation history, stars that formed during mergers at different epochs show different [alpha/Fe]. We show that the [Mg/Fe] against star formation time relation can be very useful to identify a multiple merger bulge formation scenario, provided there is sufficiently good age information available. Our simulations also show that stars formed during one of the merger events retain a systematically prograde rotation at the final time. This demonstrates that the orbit of the ancient merger that helped to form the bulge could still remain in the kinematics of bulge stars.

preprint2009arXiv

Galactic chemical evolution in hierarchical formation models - I. Early-type galaxies in the local Universe

We study the metallicities and abundance ratios of early-type galaxies in cosmological semi-analytic models (SAMs) within the hierarchical galaxy formation paradigm. To achieve this we implemented a detailed galactic chemical evolution (GCE) model and can now predict abundances of individual elements for the galaxies in the semi-analytic simulations. This is the first time a SAM with feedback from Active Galactic Nuclei (AGN) has included a chemical evolution prescription that relaxes the instantaneous recycling approximation. We find that the new models are able to reproduce the observed mass-metallicity (M*-[Z/H]) relation and, for the first time in a SAM, we reproduce the observed positive slope of the mass-abundance ratio (M*-[$alpha$/Fe]) relation. Our results indicate that in order to simultaneously match these observations of early-type galaxies, the use of both a very mildly top-heavy IMF (i.e., with a slope of x=1.15 as opposed to a standard x=1.3), and a lower fraction of binaries that explode as Type Ia supernovae appears to be required. We also examine the rate of supernova explosions in the simulated galaxies. In early-type (non-star forming) galaxies, our predictions are also consistent with the observed SNe rates. However, in star-forming galaxies, a higher fraction of SN Ia binaries than in our preferred model is required to match the data. If, however, we deviate from the classical model and introduce a population of SNe Ia with very short delay times, our models simultaneously produce a good match to the observed metallicities, abundance ratios and SN rates.

preprint2009arXiv

Structure, Kinematics, and Chemical Enrichment Patterns after Major Gas-Rich Disc-Disc Mergers

We used an N-body smoothed particle hydrodynamics algorithm, with a detailed treatment of star formation, supernovae feedback, and chemical enrichment, to perform eight simulations of mergers between gas-rich disc galaxies. We vary the mass ratio of the progenitors, their rotation axes, and their orbital parameters and analyze the kinematic, structural, and chemical properties of the remnants. Six of these simulations result in the formation of a merger remnant with a disc morphology as a result of the large gas-fraction of the remnants. We show that stars formed during the merger (a sudden starburst occur in our simulation and last for 0.2-0.3 Gyr) and those formed after the merger have different kinematical and chemical properties. The first ones are located in thick disc or the halo. They are partially supported by velocity dispersion and have high [alpha/Fe] ratios even at metallicities as high as [Fe/H]=-0.5. The former ones -- the young component -- are located in a thin disc rotationally supported and have lower [alpha/Fe] ratios. The difference in the rotational support of both components results in the rotation of the thick disc lagging that of the thin disc by as much as a factor of two, as recently observed.We find that, while the kinematic and structural properties of the merger remnant depends strongly upon the orbital parameters of the mergers, there is a remarkable uniformity in the chemical properties of the mergers. This suggests that general conclusions about the chemical signature of gas-rich mergers can be drawn.

preprint2008arXiv

Hydrodynamical Adaptive Mesh Refinement Simulations of Disk Galaxies

To date, fully cosmological hydrodynamic disk simulations to redshift zero have only been undertaken with particle-based codes, such as GADGET, Gasoline, or GCD+. In light of the (supposed) limitations of traditional implementations of smoothed particle hydrodynamics (SPH), or at the very least, their respective idiosyncrasies, it is important to explore complementary approaches to the SPH paradigm to galaxy formation. We present the first high-resolution cosmological disk simulations to redshift zero using an adaptive mesh refinement (AMR)-based hydrodynamical code, in this case, RAMSES. We analyse the temporal and spatial evolution of the simulated stellar disks' vertical heating, velocity ellipsoids, stellar populations, vertical and radial abundance gradients (gas and stars), assembly/infall histories, warps/lopsideness, disk edges/truncations (gas and stars), ISM physics implementations, and compare and contrast these properties with our sample of cosmological SPH disks, generated with GCD+. These preliminary results are the first in our long-term Galactic Archaeology Simulation program.

preprint2008arXiv

Simulating the Mass-Metallicity Relation from z=1

We use 112 N-body/hydrodynamical simulations in the standard Cold Dark Matter universe, to follow the formation of galaxy-sized halos and investigate the chemical enrichment of both the stellar component and the interstellar medium of galaxies, with stellar masses larger than 1e9 Msun. The resulting chemical properties of the simulated galaxies are broadly consistent with the observations. The predicted relationship between the mean metallicity and the galaxy stellar mass for both the stellar and the gaseous components at z=0 are in agreement with the relationships observed locally. The predicted scatter about these relationships, which is traced to the differing merging histories amongst the simulated galaxies with similar final masses, is similar to that observed. The predicted correlations between the total mass and the stellar mass of galaxies in our simulated sample from the present epoch up to z=1 agree with observed ones. The stellar mass versus metallicity relation and its associated scatter are reproduced by the simulations as consequences of the increasing efficiency of the conversion of gas into stars with stellar mass, and the differing merging histories amongst the galaxies with similar masses. The old ages of simulated low mass galaxies at z=0, and the weak level of chemical evolution for massive galaxies suggest however that our modeling of the supernova feedback may be incomplete, or that other feedback processes have been neglected. (Abridged)

preprint2006arXiv

The Spatial Distribution of the Galactic First Stars II: SPH Approach

We use cosmological, chemo-dynamical, smoothed particle hydrodynamical simulations of Milky-Way-analogue galaxies to find the expected present-day distributions of both metal-free stars that formed from primordial gas and the oldest star populations. We find that metal-free stars continue to form until z~4 in halos that are chemically isolated and located far away from the biggest progenitor of the final system. As a result, if the Population III initial mass function allows stars with low enough mass to survive until z=0 (< 0.8 Msol), they would be distributed throughout the Galactic halo. On the other hand, the oldest stars form in halos that collapsed close to the highest density peak of the final system, and at z=0 they are located preferentially in the central region of the Galaxy, i.e., in the bulge. According to our models, these trends are not sensitive to the merger histories of the disk galaxies or the implementation of supernova feedback. Furthermore, these full hydrodynamics results are consistent with our N-body results in Paper I, and lend further weight to the conclusion that surveys of low-metallicity stars in the Galactic halo can be used to directly constrain the properties of primordial stars. In particular, they suggest that the current lack of detections of metal-free stars implies that their lifetimes were shorter than a Hubble time, placing constraints on the metal-free initial mass function.

preprint2005arXiv

Gas Rich Mergers in Disk Formation

In order to explain disk galaxy formation within the hierarchical structure formation, it seems that gas rich mergers must play an important role. We review here our previous studies which have shown the importance of mergers at high redshift being gas rich, in the formation of both the stellar halo and thick disk components of disk galaxies. Regulation of star formation in the building blocks of our galaxy is required to form a low mass low metallicity stellar halo. This regulation results in high redshift, gas rich mergers during which the thick disk forms. In these proceedings, we categorise stars from our simulated disk galaxy into thin and thick disk components by using the Toomre diagram.

preprint2005arXiv

Kinematics of the Galactic Halo from Horizontal Branch stars in the Hamburg/ESO Survey

Large samples of Field Horizontal Branch (FHB) stars make excellent tracers of the Galactic halo; by studying their kinematics, one can infer important physical properties of our Galaxy. Here we present the results of a medium-resolution spectroscopic survey of 530 FHB stars selected from the Hamburg/ESO survey. The stars have a mean distance of ~7 kpc and thus probe the inner parts of the Milky Way halo. We measure radial velocities from the spectra in order to test the model of Sommer-Larsen et al., who suggested that the velocity ellipsoid of the halo changes from radially-dominated orbits to tangentially-dominated orbits as one proceeds from the inner to the outer halo. We find that the present data are unable to discriminate between this model and a more simple isothermal ellipsoid; we suggest that additional observations towards the Galactic centre might help to differentiate them.

preprint2002arXiv

Cosmology on a Mesh

An adaptive multi grid approach to simulating the formation of structure from collisionless dark matter is described. MLAPM (Multi-Level Adaptive Particle Mesh) is one of the most efficient serial codes available on the cosmological 'market' today. As part of Swinburne University's role in the development of the Square Kilometer Array, we are implementing hydrodynamics, feedback, and radiative transfer within the MLAPM adaptive mesh, in order to simulate baryonic processes relevant to the interstellar and intergalactic media at high redshift. We will outline our progress to date in applying the existing MLAPM to a study of the decay of satellite galaxies within massive host potentials.

preprint2002arXiv

Simulating the Hot X-ray Emitting Gas in Elliptical Galaxies

We study the chemo-dynamical evolution of elliptical galaxies and their hot X-ray emitting gas using high-resolution cosmological simulations. Our Tree N-body/SPH code includes a self-consistent treatment of radiative cooling, star formation, supernovae feedback, and chemical enrichment. We present a series of LCDM cosmological simulations which trace the spatial and temporal evolution of heavy element abundance patterns in both the stellar and gas components of galaxies. X-ray spectra of the hot gas are constructed via the use of the vmekal plasma model, and analysed using XSPEC with the XMM EPN response function. Simulation end-products are quantitatively compared with the observational data in both the X-ray and optical regime. We find that radiative cooling is important to interpret the observed X-ray luminosity, temperature, and metallicity of the interstellar medium of elliptical galaxies. However, this cooled gas also leads to excessive star formation at low redshift, and therefore results in underlying galactic stellar populations which are too blue with respect to observations.

preprint2002arXiv

The White Dwarf Cooling Sequence of the Globular Cluster Messier 4

We present the white dwarf sequence of the globular cluster M4, based on a 123 orbit Hubble Space Telescope exposure, with limiting magnitude V = 30, I = 28. The white dwarf luminosity function rises sharply for I >25.5, consistent with the behaviour expected for a burst population. The white dwarfs of M4 extend to approximately 2.5 magnitudes fainter than the peak of the local Galactic disk white dwarf luminosity function. This demonstrates a clear and significant age difference between the Galactic disk and the halo globular cluster M4. Using the same standard white dwarf models (Hansen 1999) to fit each luminosity function yields ages of 7.3 +/- 1.5 Gyr for the disk and 12.7 +/- 0.7 Gyr for M4 (2-sigma statistical errors).

preprint1999arXiv

A Database of Cepheid Distance Moduli and TRGB, GCLF, PNLF and SBF Data Useful for Distance Determinations

We present a compilation of Cepheid distance moduli and data for four secondary distance indicators that employ stars in the old stellar populations: the planetary nebula luminosity function (PNLF), the globular cluster luminosity function (GCLF), the tip of the red giant branch (TRGB), and the surface brightness fluctuation (SBF) method. The database includes all data published as of July 15, 1999. The main strength of this compilation resides in all data being on a consistent and homogeneous system: all Cepheid distances are derived using the same calibration of the period-luminosity relation, the treatment of errors is consistent for all indicators, measurements which are not considered reliable are excluded. As such, the database is ideal for inter-comparing any of the distance indicators considered, or for deriving a Cepheid calibration to any secondary distance indicator. Specifically, the database includes: 1) Cepheid distances, extinctions and metallicities; 2) apparent magnitudes of the PNLF cutoff; 3) apparent magnitudes and colors of the turnover of the GCLF (both in the V- and B-bands); 4) apparent magnitudes of the TRGB (in the I-band) and V-I colors at and 0.5 magnitudes fainter than the TRGB; 5) apparent surface brightness fluctuation magnitudes I, K', K_short, and using the F814W filter with the HST/WFPC2. In addition, for every galaxy in the database we give reddening estimates from DIRBE/IRAS as well as HI maps, J2000 coordinates, Hubble and T-type morphological classification, apparent total magnitude in B, and systemic velocity. (Abridged)

preprint1999arXiv

The Extragalactic Distance Scale Key Project XXVII. A Derivation of the Hubble Constant Using the Fundamental Plane and Dn-Sigma Relations in Leo I, Virgo, and Fornax

Using published photometry and spectroscopy, we construct the fundamental plane and D_n-Sigma relations in Leo I, Virgo and Fornax. The published Cepheid P-L relations to spirals in these clusters fixes the relation between angular size and metric distance for both the fundamental plane and D_n-Sigma relations. Using the locally calibrated fundamental plane, we infer distances to a sample of clusters with a mean redshift of cz \approx 6000 \kms, and derive a value of H_0=78+- 5+- 9 km/s/Mpc (random, systematic) for the local expansion rate. This value includes a correction for depth effects in the Cepheid distances to the nearby clusters, which decreased the deduced value of the expansion rate by 5% +- 5%. If one further adopts the metallicity correction to the Cepheid PL relation, as derived by the Key Project, the value of the Hubble constant would decrease by a further 6%+- 4%. These two sources of systematic error, when combined with a +- 6% error due to the uncertainty in the distance to the Large Magellanic Cloud, a +- 4% error due to uncertainties in the WFPC2 calibration, and several small sources of uncertainty in the fundamental plane analysis, combine to yield a total systematic uncertainty of +- 11%. We find that the values obtained using either the CMB, or a flow-field model, for the reference frame of the distant clusters, agree to within 1%. The Dn-Sigma relation also produces similar results, as expected from the correlated nature of the two scaling relations. A complete discussion of the sources of random and systematic error in this determination of the Hubble constant is also given, in order to facilitate comparison with the other secondary indicators being used by the Key Project.

preprint1999arXiv

The HST Key Project on the Extragalactic Distance Scale XXV. A Recalibration of Cepheid Distances to Type Ia Supernovae and the Value of the Hubble Constant

Cepheid-based distances to seven Type Ia supernovae (SNe)-host galaxies have been derived using the standard HST Key Project on the Extragalactic Distance Scale pipeline. For the first time, this allows for a transparent comparison of data accumulated as part of three different HST projects, the Key Project, the Sandage et al. Type Ia SNe program, and the Tanvir et al. Leo I Group study. Re-analyzing the Tanvir et al. galaxy and six Sandage et al. galaxies we find a mean (weighted) offset in true distance moduli of 0.12+/-0.07 mag -- i.e., 6% in linear distance -- in the sense of reducing the distance scale, or increasing H0. Adopting the reddening-corrected Hubble relations of Suntzeff et al. (1999), tied to a zero point based upon SNe~1990N, 1981B, 1998bu, 1989B, 1972E and 1960F and the photometric calibration of Hill et al. (1998), leads to a Hubble constant of H0=68+/-2(random)+/-5(systematic) km/s/Mpc. Adopting the Kennicutt et al. (1998) Cepheid period-luminosity-metallicity dependency decreases the inferred H0 by 4%. The H0 result from Type Ia SNe is now in good agreement, to within their respective uncertainties, with that from the Tully-Fisher and surface brightness fluctuation relations.

preprint1999arXiv

The HST Key Project on the Extragalactic Distance Scale XXVI. The Calibration of Population II Secondary Distance Indicators and the Value of the Hubble Constant

A Cepheid-based calibration is derived for four distance indicators that utilize stars in the old stellar populations: the tip of the red giant branch (TRGB), the planetary nebula luminosity function (PNLF), the globular cluster luminosity function (GCLF) and the surface brightness fluctuation method (SBF). The calibration is largely based on the Cepheid distances to 18 spiral galaxies within cz =1500 km/s obtained as part of the HST Key Project on the Extragalactic Distance Scale, but relies also on Cepheid distances from separate HST and ground-based efforts. The newly derived calibration of the SBF method is applied to obtain distances to four Abell clusters in the velocity range between 3800 and 5000 km/s, observed by Lauer et al. (1998) using the HST/WFPC2. Combined with cluster velocities corrected for a cosmological flow model, these distances imply a value of the Hubble constant of H0 = 69 +/- 4 (random) +/- 6 (systematic) km/s/Mpc. This result assumes that the Cepheid PL relation is independent of the metallicity of the variable stars; adopting a metallicity correction as in Kennicutt et al. (1998), would produce a (5 +/- 3)% decrease in H0. Finally, the newly derived calibration allows us to investigate systematics in the Cepheid, PNLF, SBF, GCLF and TRGB distance scales.

preprint1999arXiv

The HST Key Project on the Extragalactic Distance Scale. XXII. The Discovery of Cepheids in NGC 1326-A

We report on the detection of Cepheids and the first distance measurement to the spiral galaxy NGC 1326-A, a member of the Fornax cluster of galaxies. We have employed data obtained with the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope. Over a 49 day interval, a total of twelve V-band (F555W) and eight I-band (F814W) epochs of observation were obtained. Two photometric reduction packages, ALLFRAME and DoPHOT, have been employed to obtain photometry measures from the three Wide Field CCDs. Variability analysis yields a total of 17 Cepheids in common with both photometry datasets, with periods ranging between 10 and 50 days. Of these 14 Cepheids with high-quality lightcurves are used to fit the V and I period-luminosity relations and derive apparent distance moduli, assuming a Large Magellanic Cloud distance modulus (m-M) (LMC) = 18.50 +- 0.10 mag and color excess E(B-V) = 0.10 mag. Assuming A(V)/E(V-I) = 2.45, the DoPHOT data yield a true distance modulus to NGC 1326-A of (m-M)_0 = 31.36 +- 0.17 (random) +- 0.13 (systematic) mag, corresponding to a distance of 18.7 \pm 1.5 (random) \pm 1.2 (systematic) Mpc. The derived distance to NGC 1326-A is in good agreement with the distance derived previously to NGC 1365, another spiral galaxy member of the Fornax cluster. However the distances to both galaxies are significantly lower than to NGC 1425, a third Cepheid calibrator in the outer parts of the cluster.

preprint1998arXiv

The HST Key Project on the Extragalactic Distance Scale XII. The Discovery of Cepheids and a New Distance to NGC 2541

We report the detection of Cepheids and a new distance to the spiral galaxy NGC 2541, based on data obtained with the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope (HST). A total of 25 exposures (divided into 13 epochs) are obtained using the F555W filter (transformed to Johnson V), and nine exposures (divided into five epochs) using the F814W filter (transformed to Cousins I). Photometric reduction of the data is performed using two independent packages, DoPHOT and DAOPHOT II/ALLFRAME, which give very good agreement in the measured magnitudes. A total of 34 bona fide Cepheids, with periods ranging from 12 to over 60 days, are identified based on both sets of photometry. By fitting V and I period-luminosity relations, apparent distance moduli are derived assuming a Large Magellanic Cloud distance modulus and mean color excess of 18.50 +/- 0.10 mag and E(B-V) = 0.10 mag respectively. Adopting A(V)/E(V-I)=2.45, we obtain a true distance modulus to NGC 2541 of 30.47 +/- 0.11 (random) +/- 0.12 (systematic) mag (D = 12.4 +/- 0.6 (random) +/- 0.7 (systematic) Mpc), and a total (Galactic plus internal) mean color excess E(B-V) = 0.08 +/- 0.05 (internal error) mag.

preprint1998arXiv

The HST Key Project on the Extragalactic Distance Scale XVII. The Cepheid Distance to NGC 4725

The distance to NGC 4725 has been derived from Cepheid variables, as part of the Hubble Space Telescope Key Project on the Extragalactic Distance Scale. Thirteen F555W (V) and four F814W (I) epochs of cosmic-ray-split Wide Field and Planetary Camera 2 observations were obtained. Twenty Cepheids were discovered, with periods ranging from 12 to 49 days. Adopting a Large Magellanic Cloud distance modulus and extinction of 18.50+/-0.10 mag and E(V-I)=0.13 mag, respectively, a true reddening-corrected distance modulus (based on an analysis employing the ALLFRAME software package) of 30.50 +/- 0.16 (random) +/- 0.17 (systematic) mag was determined for NGC 4725. The corresponding of distance of 12.6 +/- 1.0 (random) +/- 1.0 (systematic) Mpc is in excellent agreement with that found with an independent analysis based upon the DoPHOT photometry package. With a foreground reddening of only E(V-I)=0.02, the inferred intrinsic reddening of this field in NGC 4725, E(V-I)=0.19, makes it one of the most highly-reddened, encountered by the HST Key Project, to date.

preprint1998arXiv

The Hubble Space Telescope Extragalactic Distance Scale Key Project. X. The Cepheid Distance to NGC 7331

The distance to NGC 7331 has been derived from Cepheid variables observed with HST/WFPC2, as part of the Extragalactic Distance Scale Key Project. Multi-epoch exposures in F555W (V) and F814W (I), with photometry derived independently from DoPHOT and DAOPHOT/ALLFRAME programs, were used to detect a total of 13 reliable Cepheids, with periods between 11 and 42 days. The relative distance moduli between NGC 7331 and the LMC, imply an extinction to NGC 7331 of A_V = 0.47+-0.15 mag, and an extinction-corrected distance modulus to NGC 7331 of 30.89+-0.14(random) mag, equivalent to a distance of 15.1 Mpc. There are additional systematic uncertainties in the distance modulus of +-0.12 mag due to the calibration of the Cepheid Period-Luminosity relation, and a systematic offset of +0.05+-0.04 mag if we applied the metallicity correction inferred from the M101 results of Kennicutt et al 1998.

preprint1997arXiv

Galactic Winds and the Photo-chemical Evolution of Elliptical Galaxies: The Classic Model Revisited

We consider the simultaneous chemical, photometric, and gaseous thermal energy evolution of elliptical galaxies. The evolution of chemical abundances in the intracluster medium (ICM) is set by the differing timescales for gas ejection, via supernovae (SNe)-driven winds, from dwarf, normal, and giant ellipticals, and is monitored concurrently. Emphasis is placed upon the influence of, and sensitivity to, the underlying stellar initial mass function (IMF), star formation efficiency, supernovae Type Ia rates, supernovae remnant (SNR) dynamics, and the most recent advances in stellar nucleosynthesis. Unlike many previous studies, we adhere to a wide range of optical (e.g. colour-metallicity-luminosity relationship) and x-ray (e.g. recent ASCA ICM abundance measurements) observational constraints. IMFs biased toward high mass stars, at least during the early phases of star formation, are implicated in order to satisfy all the observational constraints.

preprint1997arXiv

Infall Models of Elliptical Galaxies: Further Evidence for a Top-Heavy Initial Mass Function

Chemical and photometric models of elliptical galaxies with infall of primordial gas, and subsequent ejection of processed material via galactic winds, are described. Ensuring that these models are consistent with the present-day colour-luminosity relation and the measured intracluster medium (ICM) abundances, we demonstrate that the initial mass function (IMF) must be significantly flatter (i.e., x<=0.80) than the canonical Salpeter slope (i.e., x=1.35). Such a ``top-heavy'' IMF is in agreement with the earlier conclusions based upon closed-box models for elliptical galaxy evolution. On the other hand, the top-heavy IMF, in conjunction with these semi-analytic infall models, predicts an ICM gas mass which exceeds that observed by up to a factor three, in contrast with the canonical closed-box models. Time and position-dependent IMF formalisms may prove to be a fruitful avenue for future research, but those presently available in the literature are shown to be inconsistent with several important observational constraints.

preprint1997arXiv

Supernovae Types Ia/II and Intracluster Medium Enrichment

We re-examine the respective roles played by supernovae (SNe) Types Ia and II in enriching the intracluster medium (ICM) of galaxy clusters, in light of the recent downward shift of the ASCA abundance ratios of alpha-elements to iron favoured by Ishimaru & Arimoto (1997, PASJ, 49, 1). Because of this shift, Ishimaru & Arimoto conclude that >50% of the ICM iron must have originated from within Type Ia SNe progenitors. A point not appreciated in their study, nor in most previous analyses, is the crucial dependence of such a conclusion upon the adopted massive star physics. Employing several alternative Type II SN yield compilations, we demonstrate how uncertainties in the treatment of convection and mass-loss can radically alter our perception of the relative importance of Type Ia and II SNe as ICM polluters. If mass-loss of the form favoured by Maeder (1992, A&A, 264, 105) or convection of the form favoured by Arnett (1996, Supernovae and Nucleosynthesis) is assumed, the effect upon the oxygen yields would lead us to conclude that Type Ia SNe play no part in polluting the ICM, in contradiction with Ishimaru & Arimoto. Apparent dichotomies still exist (e.g. the mean ICM neon-to-iron ratio implies a 100% Type II Fe origin, while the mean sulphur ratio indicates a 100% Type Ia origin) that cannot be reconciled with the currently available yield tables.

preprint1995arXiv

CHEMICAL ABUNDANCES IN CLUSTERS OF GALAXIES

We study the origin of iron and alpha-elements (O, Mg, Si) in clusters of galaxies. In particular, we discuss the [O/Fe] ratio and the iron mass-to-luminosity ratio in the intracluster medium (ICM) and their link to the chemical and dynamical evolution of elliptical and lenticular galaxies. We adopt a detailed model of galactic evolution incorporating the development of supernovae- driven galactic winds which pollute the ICM with enriched ejecta. We demonstrate \it quantitatively \rm the crucial dependence upon the assumed stellar initial mass function in determining the evolution of the mass and abundances ratios of heavy elements in typical model ICMs. We show that completely opposite behaviours of [alpha/Fe] ratios (\ie positive versus negative ratios) can be obtained by varying the initial mass function without altering the classic assumptions regarding type Ia supernovae progenitors or their nucleosynthesis. Our results indicate that models incorporating somewhat flatter-than-Salpeter initial mass functions (ie x approx 1, as opposed to x=1.35) are preferred, provided the intracluster medium iron mass-to-luminosity ratio, preliminary [alpha/Fe]>0 ASCA results, and present-day type Ia supernovae rates, are to be matched. A simple Virgo cluster simulation which adheres to these constraints shows that approx 70% of the measured ICM iron mass has its origin in type II supernovae, with the remainder being synthesized in type Ia systems.

preprint1995arXiv

The Case Against Bimodal Star Formation in Elliptical Galaxies

We consider the present-day photometric and chemical properties of elliptical galaxies, adopting the bimodal star formation scenario of Elbaz, Arnaud & Vangioni-Flam (1995). These models utilise an initial mass function (IMF) biased heavily toward massive stars during the early phases of galactic evolution, leading to early Type II supernovae-driven galactic winds. A subsequent lengthy, milder, star formation phase with a normal IMF ensues, responsible supposedly for the stellar population observed today. Based upon chemical evolution arguments alone, this scenario has been invoked to explain the observed metal mass, and their abundance ratios, in the intracluster medium of galaxy clusters. Building upon the recent compilations of metallicity-dependent isochrones for simple stellar populations, we have constructed a coupled photometric and chemical evolution package for composite stellar populations in order to quantify the effects of such a model upon the photo-chemical properties of the resultant elliptical galaxies. We demonstrate that these predicted properties are incompatible with those observed at the current epoch.

preprint1994arXiv

Stellar Mass Loss-Driven Wind Models of Elliptical Galaxies?

Recent claims in the literature (Bressan, Chiosi & Fagotto 1994) that the epoch describing the onset of galactic winds in spheroidal star systems has been severely overestimated in the past, due to the neglect of energy deposited in the interstellar medium from stellar winds, is evaluated in the light of a more conservative approach to modeling the input of thermal energy to the system from massive stars undergoing mass loss. Applying the most recent models of stellar kinetic energy deposition, coupled with reasonable assumptions regarding the efficiency of thermalisation, it is shown that contrary to the aforementioned study, the influence of stellar winds in driving the galactic winds in ellipticals of mass M_G>1e9 M_Sun is most likely negligible.