Source author record

B. K. Gibson

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

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

82 published item(s)

preprint2026arXiv

Optimizing the potential of KM3NeT in detecting core-collapse supernovae

Core-collapse supernovae mark the end of life of massive stars. However, despite their importance in astrophysics, their underlying mechanisms remain unclear. Neutrinos that emerge from the dense core of the star offer a promising way to study supernova dynamics. A strategy is presented to improve the potential of the KM3NeT neutrino telescope to detect core-collapse supernovae in our Galaxy or the Large Magellanic Cloud by further exploiting the properties of its optical modules equipped with multiple photomultipliers. A supernova burst is expected to produce a sudden hit rate increase in the KM3NeT detectors. New observables have been defined for individual optical modules that exploit the geometry and time distribution of the detected hits, enabling a better discrimination between signal and background signatures. In addition, a thorough investigation of the related systematic uncertainties is presented for the first time. When implemented, this new methodology allowed KM3NeT to probe 46% more Galactic core-collapse supernova candidates than with the previous trigger strategy, reaching the dense Galactic bulge. It is now expected that, once completed, KM3NeT will achieve full Galactic sensitivity to core-collapse supernovae.

preprint2020arXiv

The RAdial Velocity Experiment (RAVE): Parameterisation of RAVE spectra based on convolutional neural networks

In the context of large spectroscopic surveys of stars, data-driven methods are key in deducing physical parameters for millions of spectra in a short time. Convolutional neural networks (CNNs) enable us to connect observables (e.g. spectra, stellar magnitudes) to physical properties (atmospheric parameters, chemical abundances, or labels in general). We trained a CNN, adopting stellar atmospheric parameters and chemical abundances from APOGEE DR16 (resolution R=22500) data as training set labels. As input, we used parts of the intermediate-resolution RAVE DR6 spectra (R~7500) overlapping with the APOGEE DR16 data as well as broad-band ALL_WISE and 2MASS photometry, together with Gaia DR2 photometry and parallaxes. We derived precise atmospheric parameters Teff, log(g), and [M/H] along with the chemical abundances of [Fe/H], [alpha/M], [Mg/Fe], [Si/Fe], [Al/Fe], and [Ni/Fe] for 420165 RAVE spectra. The precision typically amounts to 60K in Teff, 0.06 in log(g) and 0.02-0.04 dex for individual chemical abundances. Incorporating photometry and astrometry as additional constraints substantially improves the results in terms of the accuracy and precision of the derived labels. We provide a catalogue of CNN-trained atmospheric parameters and abundances along with their uncertainties for 420165 stars in the RAVE survey. CNN-based methods provide a powerful way to combine spectroscopic, photometric, and astrometric data without the need to apply any priors in the form of stellar evolutionary models. The developed procedure can extend the scientific output of RAVE spectra beyond DR6 to ongoing and planned surveys such as Gaia RVS, 4MOST, and WEAVE. We call on the community to place a particular collective emphasis and on efforts to create unbiased training samples for such future spectroscopic surveys.

preprint2016arXiv

Galactic Archaeology and Minimum Spanning Trees

Chemical tagging of stellar debris from disrupted open clusters and associations underpins the science cases for next-generation multi-object spectroscopic surveys. As part of the Galactic Archaeology project TraCD (Tracking Cluster Debris), a preliminary attempt at reconstructing the birth clouds of now phase-mixed thin disk debris is undertaken using a parametric minimum spanning tree (MST) approach. Empirically-motivated chemical abundance pattern uncertainties (for a 10-dimensional chemistry-space) are applied to NBODY6-realised stellar associations dissolved into a background sea of field stars, all evolving in a Milky Way potential. We demonstrate that significant population reconstruction degeneracies appear when the abundance uncertainties approach 0.1 dex and the parameterised MST approach is employed; more sophisticated methodologies will be required to ameliorate these degeneracies.

preprint2016arXiv

Galactic Chemical Evolution

We analyze the evolution of oxygen abundance radial gradients resulting from our chemical evolution models calculated with different prescriptions for the star formation rate (SFR) and for the gas infall rate, in order to assess their respective roles in shaping gradients. We also compare with cosmological simulations and confront all with recent observational datasets, in particular with abundances inferred from planetary nebulae. We demonstrate the critical importance in isolating the specific radial range over which a gradient is measured, in order for their temporal evolution to be useful indicators of disk growth with redshift.

preprint2016arXiv

HI4PI: A full-sky HI survey based on EBHIS and GASS

Measurement of the Galactic neutral atomic hydrogen (HI) column density, NHI, and brightness temperatures, Tb, is of high scientific value for a broad range of astrophysical disciplines. In the past two decades, one of the most-used legacy HI datasets has been the Leiden/Argentine/Bonn Survey (LAB). We release the HI 4$π$ survey (HI4PI), an all-sky database of Galactic HI, which supersedes the LAB survey. The HI4PI survey is based on data from the recently completed first coverage of the Effelsberg-Bonn HI Survey (EBHIS) and from the third revision of the Galactic All-Sky Survey (GASS). EBHIS and GASS share similar angular resolution and match well in sensitivity. Combined, they are ideally suited to be a successor to LAB. The new HI4PI survey outperforms the LAB in angular resolution (16.2', FWHM) and sensitivity (RMS: 43 mK). Moreover, it has full spatial sampling and thus overcomes a major drawback of LAB, which severely undersamples the sky. We publish all-sky column density maps of the neutral atomic hydrogen in the Milky Way, along with full spectroscopic data, in several map projections including HEALPix.

preprint2016arXiv

Identification of Globular Cluster Stars in RAVE data II: Extended tidal debris around NGC 3201

We report the identification of extended tidal debris potentially associated with the globular cluster NGC 3201, using the RAVE catalogue. We find the debris stars are located at a distance range of 1-7 kpc based on the forthcoming RAVE distance estimates. The derived space velocities and integrals of motion show interesting connections to NGC 3201, modulo uncertainties in the proper motions. Three stars, which are among the 4 most likely candidates for NGC 3201 tidal debris, are separated by 80 degrees on the sky yet are well matched by the 12 Gyr, [Fe/H] = -1.5 isochrone appropriate for the cluster. This is the first time tidal debris around this cluster has been reported over such a large spatial extent, with implications for the cluster$'$s origin and dynamical evolution.

preprint2016arXiv

Impacts of a Flaring Star-forming Disc and Stellar Radial Mixing on the Vertical Metallicity Gradient

Using idealised N-body simulations of a Milky Way-sized disc galaxy, we qualitatively study how the metallicity distributions of the thin disc star particles are modified by the formation of the bar and spiral arm structures. The thin disc in our numerical experiments initially has a tight negative radial metallicity gradient and a constant vertical scale-height. We show that the radial mixing of stars drives a positive vertical metallicity gradient in the thin disc. On the other hand, if the initial thin disc is flared, with vertical scale-height increasing with galactocentric radius, the metal poor stars originally in the outer disc become dominant in regions above the disc plane at every radii. This process can drive a negative vertical metallicity gradient, which is consistent with the current observed trend. This model mimics a scenario where the star-forming thin disc was flared in the outer region at earlier epochs. Our numerical experiment with an initial flared disc predicts that the negative vertical metallicity gradient of the mono-age relatively young thin disc population should be steeper in the inner disc, and the radial metallicity gradient of the mono-age population should be shallower at greater heights above the disc plane. We also predict that the metallicity distribution function of mono-age young thin disc populations above the disc plane would be more positively skewed in the inner disc compared to the outer disc.

preprint2016arXiv

The evolution of the radial gradient of Oxygen abundance in spiral galaxies

The aim of this work is to present our new series of chemical evolution models computed for spiral and low mass galaxies of different total masses and star formation efficiencies. We analyze the results of models, in particular the evolution of the radial gradient of oxygen abundance. Furthermore, we study the role of the infall rate and of the star formation history on the variations of this radial gradient. The relations between the O/H radial gradient and other spiral galaxies characteristics as the size or the stellar mass are also shown. We find that the radial gradient is mainly a scale effect which basically does not change with the redshift (or time) if it is measured within the optical radius. Moreover, when it is measured as a function of a normalized radius, show a similar value for all galaxies masses, showing a correlation with a dispersion around an average value which is due to the differences star formation efficiencies, in agreement with the idea of an universal O/H radial gradient

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.

preprint2016arXiv

Weighing the Galactic disk using the Jeans equation: lessons from simulations

Using three-dimensional stellar kinematic data from simulated galaxies, we examine the efficacy of a Jeans equation analysis in reconstructing the total disk surface density, including the dark matter, at the "Solar" radius. Our simulation dataset includes galaxies formed in a cosmological context using state-of-the-art high resolution cosmological zoom simulations, and other idealised models. The cosmologically formed galaxies have been demonstrated to lie on many of the observed scaling relations for late-type spirals, and thus offer an interesting surrogate for real galaxies with the obvious advantage that all the kinematical data are known perfectly. We show that the vertical velocity dispersion is typically the dominant kinematic quantity in the analysis, and that the traditional method of using only the vertical force is reasonably effective at low heights above the disk plane. At higher heights the inclusion of the radial force becomes increasingly important. We also show that the method is sensitive to uncertainties in the measured disk parameters, particularly the scale lengths of the assumed double exponential density distribution, and the scale length of the radial velocity dispersion. In addition, we show that disk structure and low number statistics can lead to significant errors in the calculated surface densities. Finally we examine the implications of our results for previous studies of this sort, suggesting that more accurate measurements of the scale lengths may help reconcile conflicting estimates of the local dark matter density in the literature.

preprint2015arXiv

Carbon-rich presolar grains from massive stars. Subsolar 12C/13C and 14N/15N ratios and the mystery of 15N

Carbon-rich grains with isotopic anomalies compared to the Sun are found in primitive meteorites. They were made by stars, and carry the original stellar nucleosynthesis signature. Silicon carbide grains of Type X and C, and low-density graphites condensed in the ejecta of core-collapse supernovae. We present a new set of models for the explosive He shell and compare them with the grains showing 12C/13C and 14N/15N ratios lower than solar. In the stellar progenitor H was ingested into the He shell and not fully destroyed before the explosion. Different explosion energies and H concentrations are considered. If the SN shock hits the He-shell region with some H still present, the models can reproduce the C and N isotopic signatures in C-rich grains. Hot-CNO cycle isotopic signatures are obtained, including a large production of 13C and 15N. The short-lived radionuclides 22Na and 26Al are increased by orders of magnitude. The production of radiogenic 22Ne from the decay of 22Na in the He shell might solve the puzzle of the Ne-E(L) component in low-density graphite grains. This scenario is attractive for the SiC grains of type AB with 14N/15N ratios lower than solar, and provides an alternative solution for SiC grains originally classified as nova grains. Finally, this process may contribute to the production of 14N and 15N in the Galaxy, helping to produce the 14N/15N ratio in the solar system.

preprint2015arXiv

Identification of Globular Cluster Stars in RAVE data I: Application to Stellar Parameter Calibration

We present the identification of potential members of nearby Galactic globular clusters using radial velocities from the RAdial Velocity Experiment Data Release 4 (RAVE-DR4) survey database. Our identifications are based on three globular clusters -- NGC 3201, NGC 5139 ($ω$ Cen) and NGC 362 -- all of which are shown to have |RV|>100 km/s. The identification of globular cluster stars in RAVE DR4 data offers a unique opportunity to test the precision and accuracy of the stellar parameters determined with the currently available Stellar Parameter Pipelines (SPPs) used in the survey, as globular clusters are ideal testbeds for the validation of stellar atmospheric parameters, abundances, distances and ages. For both NGC 3201 and $ω$ Cen, there is compelling evidence for numerous members (> 10) in the RAVE database; in the case of NGC 362 the evidence is more ambiguous, and there may be significant foreground and/or background contamination in our kinematically-selected sample. A comparison of the RAVE-derived stellar parameters and abundances with published values for each cluster and with BASTI isochrones for ages and metallicities from the literature reveals overall good agreement, with the exception of the apparent underestimation of surface gravities for giants, in particular for the most metal-poor stars. Moreover, if the selected members are part of the main body of each cluster our results would also suggest that the distances from Binney et al. 2013, where only isochrones more metal-rich than -0.9 dex were used, are typically underestimated by ~ 40% with respect to the published distances for the clusters, while the distances from Zwitter et al. 2010 show stars ranging from 1 to ~ 6.5 kpc -- with indications of a trend toward higher distances at lower metallicities -- for the three clusters analysed in this study.

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.

preprint2015arXiv

The imprint of satellite accretion on the chemical and dynamical properties of disc galaxies

Aims: We study the effects of the cosmological assembly history on the chemical and dynamical properties of the discs of spiral galaxies as a function of radius. Methods: We make use of the simulated Milky-Way mass, fully-cosmological discs, from {\tt RaDES} (Ramses Disc Environment Study). We analyse their assembly history by examining the proximity of satellites to the galactic disc, instead of their merger trees, to better gauge which satellites impact the disc. We present stellar age and metallicity profiles, Age-Metallicity Relation (AMR), Age-Velocity dispersion Relation (AVR), and Stellar Age Distribution (SAD) in several radial bins for the simulated galaxies. Results: Assembly histories can be divided into three different stages: i) a merger dominated phase, when a large number of mergers with mass ratios of $\sim$1:1 take place (lasting $\sim$3.2$\pm$0.4 Gyr on average); ii) a quieter phase, when $\sim$1:10 mergers take place (lasting $\sim$4.4$\pm$2.0 Gyr) - these two phases are able to kinematically heat the disc and produce a disc that is chemically mixed over its entire radial extension; iii) a "secular" phase where the few mergers that take place have mass ratios below 1:100, and which do not affect the disc properties (lasting $\sim$5.5$\pm$2.0 Gyr). Phase ii ends with a final merger event (at time $t_\mathrm{jump}$) marking the onset of important radial differences in the AMR, AVR, and SAD. Conclusions: Inverted AMR trends in the outer parts of discs, for stars younger than $t_\mathrm{jump}$, are found as the combined effect of radial motions and star formation in satellites temporarily located in these outer parts. "U-shaped" stellar age profiles change to an old plateau ($\sim$10 Gyr) in the outer discs for the entire {\tt RaDES} sample. This shape is a consequence of inside-out growth of the disc, radial motions of disc stars ... [abridged]

preprint2015arXiv

The imprints of the Galactic bar on the thick disk with RAVE

We study the kinematics of a local sample of stars, located within a cylinder of 500 pc radius centered on the Sun, in the RAVE dataset. We find clear asymmetries in the $v_R$-$v_ϕ$ velocity distributions of thin and thick disk stars: here are more stars moving radially outwards for low azimuthal velocities and more radially inwards for high azimuthal velocities. Such asymmetries have been previously reported for the thin disk as being due to the Galactic bar, but this is the first time that the same type of structures are seen in the thick disk. Our findings imply that the velocities of thick disk stars should no longer be described by Schwarzschild's, multivariate Gaussian or purely axisymmetric distributions. Furthermore, the nature of previously reported substructures in the thick disk needs to be revisited as these could be associated with dynamical resonances rather than to accretion events. It is clear that dynamical models of the Galaxy must fit the 3D velocity distributions of the disks, rather than the projected 1D, if we are to understand the Galaxy fully.

preprint2014arXiv

APASS Landolt-Sloan BVgri photometry of RAVE stars. I. Data, effective temperatures and reddenings

We provide APASS photometry in the Landolt BV and Sloan g'r'i' bands for all the 425,743 stars included in the latest 4th RAVE Data Release. The internal accuracy of the APASS photometry of RAVE stars, expressed as error of the mean of data obtained and separately calibrated over a median of 4 distinct observing epochs and distributed between 2009 and 2013, is 0.013, 0.012, 0.012, 0.014 and 0.021 mag for B, V, g', r' and i' band, respectively. The equally high external accuracy of APASS photometry has been verified on secondary Landolt and Sloan photometric standard stars not involved in the APASS calibration process, and on a large body of literature data on field and cluster stars, confirming the absence of offsets and trends. Compared with the Carlsberg Meridian Catalog (CMC-15), APASS astrometry of RAVE stars is accurate to a median value of 0.098 arcsec. Brightness distribution functions for the RAVE stars have been derived in all bands. APASS photometry of RAVE stars, augmented by 2MASS JHK infrared data, has been chi2 fitted to a densely populated synthetic photometric library designed to widely explore in temperature, surface gravity, metallicity and reddening. Resulting Teff and E(B-V), computed over a range of options, are provided and discussed, and will be kept updated in response to future APASS and RAVE data releases. In the process it is found that the reddening caused by an homogeneous slab of dust, extending for 140 pc on either side of the Galactic plane and responsible for E(B-V,poles)=0.036 +/- 0.002 at the galactic poles, is a suitable approximation of the actual reddening encountered at Galactic latitudes |b|>=25 deg.

preprint2014arXiv

Characterizing the High-Velocity Stars of RAVE: The Discovery of a Metal-Rich Halo Star Born in the Galactic Disk

We aim to characterize high-velocity (HiVel) stars in the solar vicinity both chemically and kinematically using the fourth data release of the RAdial Velocity Experiment (RAVE). We used a sample of 57 HiVel stars with Galactic rest-frame velocities larger than 275 km s$^{-1}$. With 6D position and velocity information, we integrated the orbits of the HiVel stars and found that, on average, they reach out to 13 kpc from the Galactic plane and have relatively eccentric orbits consistent with the Galactic halo. Using the stellar parameters and [$α$/Fe] estimates from RAVE, we found the metallicity distribution of the HiVel stars peak at [M/H] = -1.2 dex and is chemically consistent with the inner halo. There are a few notable exceptions that include a hypervelocity star (HVS) candidate, an extremely high-velocity bound halo star, and one star that is kinematically consistent with the halo but chemically consistent with the disk. High-resolution spectra were obtained for the metal-rich HiVel star candidate and the second highest velocity star in the sample. Using these high-resolution data, we report the discovery of a metal-rich halo star that has likely been dynamically ejected into the halo from the Galactic thick disk. This discovery could aid in explaining the assembly of the most metal-rich component of the Galactic halo.

preprint2014arXiv

Chemical gradients in the Milky Way from the RAVE data. II. Giant stars

We provide new constraints on the chemo-dynamical models of the Milky Way by measuring the radial and vertical chemical gradients for the elements Mg, Al, Si, Ti, and Fe in the Galactic disc and the gradient variations as a function of the distance from the Galactic plane ($Z$). We selected a sample of giant stars from the RAVE database using the gravity criterium 1.7$<$log g$<$2.8. We created a RAVE mock sample with the Galaxia code based on the Besan\c con model and selected a corresponding mock sample to compare the model with the observed data. We measured the radial gradients and the vertical gradients as a function of the distance from the Galactic plane $Z$ to study their variation across the Galactic disc. The RAVE sample exhibits a negative radial gradient of $d[Fe/H]/dR=-0.054$ dex kpc$^{-1}$ close to the Galactic plane ($|Z|<0.4$ kpc) that becomes flatter for larger $|Z|$. Other elements follow the same trend although with some variations from element to element. The mock sample has radial gradients in fair agreement with the observed data. The variation of the gradients with $Z$ shows that the Fe radial gradient of the RAVE sample has little change in the range $|Z|\lesssim0.6$ kpc and then flattens. The iron vertical gradient of the RAVE sample is slightly negative close to the Galactic plane and steepens with $|Z|$. The mock sample exhibits an iron vertical gradient that is always steeper than the RAVE sample. The mock sample also shows an excess of metal-poor stars in the [Fe/H] distributions with respect to the observed data. These discrepancies can be reduced by decreasing the number of thick disc stars and increasing their average metallicity in the Besan\c con model.

preprint2014arXiv

Chemodynamics of a Simulated Disc Galaxy: Initial Mass Functions and SNIa Progenitors

We trace the formation and advection of several elements within a cosmological adaptive mesh refinement simulation of an L* galaxy. We use nine realisations of the same initial conditions with different stellar Initial Mass Functions (IMFs), mass limits for type-II and type-Ia supernovae (SNII, SNIa) and stellar lifetimes to constrain these sub-grid phenomena. Our code includes self-gravity, hydrodynamics, star formation, radiative cooling and feedback from multiple sources within a cosmological framework. Under our assumptions of nucleosynthesis we find that SNII with progenitor masses of up to 100 Msun are required to match low metallicity gas oxygen abundances. Tardy SNIa are necessary to reproduce the classical chemical evolution knee in [O/Fe]-[Fe/H]: more prompt SNIa delayed time distributions do not reproduce this feature. Within our framework of hydrodynamical mixing of metals and galaxy mergers we find that chemical evolution is sensitive to the shape of the IMF and that there exists a degeneracy with the mass range of SNII. We look at the abundance plane and present the properties of different regions of the plot, noting the distinct chemical properties of satellites and a series of nested discs that have greater velocity dispersions, are more alpha-rich and metal poor with age.

preprint2014arXiv

Constraints on the Galactic Bar from the Hercules stream as traced with RAVE across the Galaxy

Non-axisymmetries in the Galactic potential (spiral arms and bar) induce kinematic groups such as the Hercules stream. Assuming that Hercules is caused by the effects of the Outer Lindblad Resonance of the Galactic bar, we model analytically its properties as a function of position in the Galaxy and its dependence on the bar's pattern speed and orientation. Using data from the RAVE survey we find that the azimuthal velocity of the Hercules structure decreases as a function of Galactocentric radius, in a manner consistent with our analytical model. This allows us to obtain new estimates of the parameters of the Milky Way's bar. The combined likelihood function of the bar's pattern speed and angle has its maximum for a pattern speed of $Ω_b=(1.89 \pm 0.08) \times Ω_0$ where $Ω_0$ is the local circular frequency. Assuming a Solar radius of 8.05 kpc and a local circular velocity of 238 km/s, this corresponds to $Ω_b=56\pm 2$ km/s/kpc. On the other hand, the bar's orientation $ϕ_b$ cannot be constrained with the available data. In fact, the likelihood function shows that a tight correlation exists between the pattern speed and the orientation, implying that a better description of our best fit results is given by the linear relation $Ω_b/Ω_0=1.91 + 0.0044[ϕ_b(°)-48]$, with standard deviation of 0.02. For example, for an angle of $ϕ_b=30°$ the pattern speed is $54.0\pm0.5$ km/s/kpc. These results are not very sensitive to the other Galactic parameters such as the circular velocity curve or the peculiar motion of the Sun, and are robust to biases in distance.

preprint2014arXiv

Kinematic modelling of the Milky Way using the RAVE and GCS stellar surveys

We investigate the kinematic parameters of the Milky Way disc using the RAVE and GCS stellar surveys. We do this by fitting a kinematic model to the data taking the selection function of the data into account. For stars in the GCS we use all phase-space coordinates, but for RAVE stars we use only $(l,b,v_{\rm los})$. Using MCMC technique, we investigate the full posterior distributions of the parameters given the data. We investigate the `age-velocity dispersion' relation for the three kinematic components ($σ_R,σ_ϕ,σ_z$), the radial dependence of the velocity dispersions, the Solar peculiar motion ($U_{\odot},V_{\odot}, W_{\odot} $), the circular speed $Θ_0$ at the Sun and the fall of mean azimuthal motion with height above the mid-plane. We confirm that the Besançon-style Gaussian model accurately fits the GCS data, but fails to match the details of the more spatially extended RAVE survey. In particular, the Shu distribution function (DF) handles non-circular orbits more accurately and provides a better fit to the kinematic data. The Gaussian distribution function not only fits the data poorly but systematically underestimates the fall of velocity dispersion with radius. We find that correlations exist between a number of parameters, which highlights the importance of doing joint fits. The large size of the RAVE survey, allows us to get precise values for most parameters. However, large systematic uncertainties remain, especially in $V_{\odot}$ and $Θ_0$. We find that, for an extended sample of stars, $Θ_0$ is underestimated by as much as $10\%$ if the vertical dependence of the mean azimuthal motion is neglected. Using a simple model for vertical dependence of kinematics, we find that it is possible to match the Sgr A* proper motion without any need for $V_{\odot}$ being larger than that estimated locally by surveys like GCS.

preprint2014arXiv

Spectroscopic Signatures of Extra-Tidal Stars Around the Globular Clusters NGC 6656 (M22), NGC 3201 and NGC 1851 from RAVE

Stellar population studies of globular clusters have suggested that the brightest clusters in the Galaxy might actually be the remnant nuclei of dwarf spheroidal galaxies. If the present Galactic globular clusters formed within larger stellar systems, they are likely surrounded by extra-tidal halos and/or tails made up of stars that were tidally stripped from their parent systems. The stellar surroundings around globular clusters are therefore one of the best places to look for the remnants of an ancient dwarf galaxy. Here an attempt is made to search for tidal debris around the supernovae enriched globular clusters M22 and NGC 1851 as well as the kinematically unique cluster NGC 3201. The stellar parameters from the Radial Velocity Experiment (RAVE) are used to identify stars with RAVE metallicities, radial velocities and elemental-abundances consistent with the abundance patterns and properties of the stars in M22, NGC 1851 and NGC 3201. The discovery of RAVE stars that may be associated with M22 and NGC 1851 are reported, some of which are at projected distances of ~10 degrees away from the core of these clusters. Numerous RAVE stars associated with NGC 3201 suggest that either the tidal radius of this cluster is underestimated, or that there are some unbound stars extending a few arc minutes from the edge of the cluster's radius. No further extra-tidal stars associated with NGC 3201 could be identified. The bright magnitudes of the RAVE stars make them easy targets for high resolution follow-up observations, allowing an eventual further chemical tagging to solidify (or exclude) stars outside the tidal radius of the cluster as tidal debris. In both our radial velocity histograms of the regions surrounding NGC 1851 and NGC 3201, a peak of stars at 230 km/s is seen, consistent with extended tidal debris from omega Centauri.

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.

preprint2014arXiv

The Rich Are Different: Evidence from the RAVE Survey for Stellar Radial Migration

Using the RAdial Velocity Experiment fourth data release (RAVE DR4), and a new metallicity calibration that will be also taken into account in the future RAVE DR5, we investigate the existence and the properties of super-solar metallicity stars ([M/H] > +0.1 dex) in the sample, and in particular in the Solar neighbourhood. We find that RAVE is rich in super-solar metallicity stars, and that the local metallicity distribution function declines remarkably slowly up to +0.4 dex. Our results show that the kinematics and height distributions of the super-solar metallicity stars are identical to those of the [M/H] < 0 thin-disc giants that we presume were locally manufactured. The eccentricities of the super-solar metallicity stars indicate that half of them are on a roughly circular orbit (e < 0.15), so under the assumption that the metallicity of the interstellar medium at a given radius never decreases with time, they must have increased their angular momenta by scattering at corotation resonances of spiral arms from regions far inside the Solar annulus. The likelihood that a star will migrate radially does not seem to decrease significantly with increasing amplitude of vertical oscillations within range of oscillation amplitudes encountered in the disc.

preprint2014arXiv

The Role of Feedback in Shaping the Structure of the Interstellar Medium

We present an analysis of the role of feedback in shaping the neutral hydrogen (HI) content of simulated disc galaxies. For our analysis, we have used two realisations of two separate Milky Way-like (~L*) discs - one employing a conservative feedback scheme (MUGS), the other significantly more energetic (MaGICC). To quantify the impact of these schemes, we generate zeroth moment (surface density) maps of the inferred HI distribution; construct power spectra associated with the underlying structure of the simulated cold ISM, in addition to their radial surface density and velocity dispersion profiles. Our results are compared with a parallel, self-consistent, analysis of empirical data from THINGS (The HI Nearby Galaxy Survey). Single power-law fits (P~k^gamma) to the power spectra of the stronger-feedback (MaGICC) runs (over spatial scales corresponding to 0.5 kpc to 20 kpc) result in slopes consistent with those seen in the THINGS sample (gamma = -2.5). The weaker-feedback (MUGS) runs exhibit shallower power law slopes (gamma = -1.2). The power spectra of the MaGICC simulations are more consistent though with a two-component fit, with a flatter distribution of power on larger scales (i.e., gamma = -1.4 for scales in excess of 2 kpc) and a steeper slope on scales below 1 kpc (gamma = -5), qualitatively consistent with empirical claims, as well as our earlier work on dwarf discs. The radial HI surface density profiles of the MaGICC discs show a clear exponential behaviour, while those of the MUGS suite are essentially flat; both behaviours are encountered in nature, although the THINGS sample is more consistent with our stronger (MaGICC) feedback runs.

preprint2014arXiv

Weighing the local dark matter with RAVE red clump stars

We determine the Galactic potential in the solar neigbourhood from RAVE observations. We select red clump stars for which accurate distances, radial velocities, and metallicities have been measured. Combined with data from the 2MASS and UCAC catalogues, we build a sample of 4600 red clump stars within a cylinder of 500 pc radius oriented in the direction of the South Galactic Pole, in the range of 200 pc to 2000 pc distances. We deduce the vertical force and the total mass density distribution up to 2 kpc away from the Galactic plane by fitting a distribution function depending explicitly on three isolating integrals of the motion in a separable potential locally representing the Galactic one with four free parameters. Because of the deep extension of our sample, we can determine nearly independently the dark matter mass density and the baryonic disc surface mass density. We find (i) at 1kpc Kz/(2piG) = 68.5 pm 1.0 Msun/pc2, and (ii) at 2 kpc Kz/(2piG) = 96.9 pm 2.2 Msun/pc2. Assuming the solar Galactic radius at R0 = 8.5 kpc, we deduce the local dark matter density rhoDM (z=0) = 0.0143 pm 0.0011Msun pc3 = 0.542 pm 0.042 Gev/cm3 and the baryonic surface mass density Sigma = 44.4 pm 4.1 Msun/pc2 . Our results are in agreement with previously published Kz determinations up to 1 kpc, while the extension to 2 kpc shows some evidence for an unexpectedly large amount of dark matter. A flattening of the dark halo of order 0.8 can produce such a high local density in combination with a circular velocity of 240 km/s . Another explanation, allowing for a lower circular velocity, could be the presence of a secondary dark component, a very thick disc resulting either from the deposit of dark matter from the accretion of multiple small dwarf galaxies, or from the presence of an effective phantom thick disc in the context of effective galactic-scale modifications of gravity.

preprint2013arXiv

A new stellar chemo-kinematic relation reveals the merger history of the Milky Way disc

The velocity dispersions of stars near the Sun are known to increase with stellar age, but age can be difficult to determine so a proxy like the abundance of alpha elements (e.g., Mg) with respect to iron, [alpha/Fe], is used. Here we report an unexpected behavior found in the velocity dispersion of a sample of giant stars from the RAdial Velocity Experiment (RAVE) survey with high quality chemical and kinematical information, in that it decreases strongly for stars with [Mg/Fe] > 0.4 dex (i.e., those that formed in the first Gyr of the Galaxy's life). These findings can be explained by perturbations from massive mergers in the early Universe, which have affected more strongly the outer parts of the disc, and the subsequent radial migration of stars with cooler kinematics from the inner disc. Similar reversed trends in velocity dispersion are also found for different metallicity subpopulations. Our results suggest that the Milky Way disc merger history can be recovered by relating the observed chemo-kinematic relations to the properties of past merger events.

preprint2013arXiv

A RAVE Investigation on Galactic open Clusters I. Radial velocities and metallicities

Context. Galactic open clusters (OCs) mainly belong to the young stellar population in the Milky Way disk, but are there groups and complexes of OCs that possibly define an additional level in hierarchical star formation? Current compilations are too incomplete to address this question, especially regarding radial velocities (RVs) and metallicities ($[M/H]$). Aims. Here we provide and discuss newly obtained RV and $[M/H]$ data, which will enable us to reinvestigate potential groupings of open clusters and associations. Methods.We extracted additional RVs and $[M/H]$ from the RAdial Velocity Experiment (RAVE) via a cross-match with the Catalogue of Stars in Open Cluster Areas (CSOCA). For the identified OCs in RAVE we derived RV and $[M/H]$ from a cleaned working sample and compared the results with previous findings. Results. Although our RAVE sample does not show the same accuracy as the entire survey, we were able to derive reliable RV for 110 Galactic open clusters. For 37 OCs we publish RV for the first time. Moreover, we determined $[M/H]$ for 81 open clusters, extending the number of OCs with $[M/H]$ by 69.

preprint2013arXiv

Chemical gradients in the Milky Way from the RAVE data. I. Dwarf stars

Aim: We aim at measuring the chemical gradients of the elements Mg, Al, Si, and Fe along the Galactic radius to provide new constraints on the chemical evolution models of the Galaxy and Galaxy models such as the Besancon model. Methods: We analysed three different samples selected from three independent datasets: a sample of 19,962 dwarf stars selected from the RAVE database, a sample of 10,616 dwarf stars selected from the Geneva-Copenhagen Survey (GCS) dataset, and a mock sample (equivalent to the RAVE sample) created by using the GALAXIA code, which is based on the Besancon model. We measured the chemical gradients as functions of the guiding radius (Rg) at different distances from the Galactic plane reached by the stars along their orbit (Zmax). Results: The chemical gradients of the RAVE and GCS samples are negative and show consistent trends, although they are not equal: at Zmax<0.4 kpc and 4.5<Rg(kpc)<9.5, the iron gradient for the RAVE sample is d[Fe/H]/dRg=-0.065 dex kpc^{-1}, whereas for the GCS sample it is d[Fe/H]/dRg=-0.043 dex kpc^{-1} with internal errors +-0.002 and +-0.004 dex kpc^{-1}, respectively. The gradients of the RAVE and GCS samples become flatter at larger Zmax. Conversely, the mock sample has a positive iron gradient of d[Fe/H]/dRg=+0.053+-0.003 dex kpc^{-1} at Zmax<0.4 kpc and remains positive at any Zmax. These positive and unrealistic values originate from the lack of correlation between metallicity and tangential velocity in the Besancon model. The discrepancies between the observational samples and the mock sample can be reduced by i) decreasing the density, ii) decreasing the vertical velocity, and iii) increasing the metallicity of the thick disc in the Besancon model.

preprint2013arXiv

Constraining Sub-Grid Physics with High-Redshift Spatially-Resolved Metallicity Distributions

We examine the role of energy feedback in shaping the distribution of metals within cosmological hydrodynamical simulations of L* disc galaxies. While negative abundance gradients today provide a boundary condition for galaxy evolution models, in support of inside-out disc growth, empirical evidence as to whether abundance gradients steepen or flatten with time remains highly contradictory. We made use of a suite of L* discs, realised with and without `enhanced' feedback. All the simulations were produced using the smoothed particle hydrodynamics code Gasoline, and their in situ gas-phase metallicity gradients traced from redshift z~2 to the present-day. Present-day age-metallicity relations and metallicity distribution functions were derived for each system. The `enhanced' feedback models, which have been shown to be in agreement with a broad range of empirical scaling relations, distribute energy and re-cycled ISM material over large scales and predict the existence of relatively `flat' and temporally invariant abundance gradients. Enhanced feedback schemes reduce significantly the scatter in the local stellar age-metallicity relation and, especially, the [O/Fe]-[Fe/H] relation. The local [O/Fe] distribution functions for our L* discs show clear bimodality, with peaks at [O/Fe]=-0.05 and +0.05 (for stars with [Fe/H]>-1), consistent with our earlier work on dwarf discs. Our results with `enhanced' feedback are inconsistent with our earlier generation of simulations realised with `conservative' feedback. We conclude that spatially-resolved metallicity distributions, particularly at high-redshift, offer a unique and under-utilised constraint on the uncertain nature of stellar feedback processes.

preprint2013arXiv

Diffuse interstellar band at 8620 Å in RAVE: A new method for detecting the diffuse interstellar band in spectra of cool stars

Diffuse interstellar bands are usually observed in spectra of hot stars, where interstellar lines are rarely blended with stellar ones. The need for hot stars is a strong limitation in the number of sightlines we can observe and the distribution of sightlines in the Galaxy, as hot stars are rare and concentrated in the Galactic plane. We are introducing a new method, where interstellar lines can be observed in spectra of cool stars in large spectroscopic surveys. The method is completely automated and does not require prior knowledge of the stellar parameters. If known, the stellar parameters only reduce the computational time and are not involved in the extraction of the interstellar spectrum. The main step in extracting interstellar lines is a construction of the stellar spectrum, which is in our method done by finding other observed spectra that lack interstellar features and are otherwise very similar to the spectrum in question. Such spectra are then combined into a single stellar spectrum template, which matches the stellar component in an observed spectrum. We demonstrate the performance of this new method on a sample of 482,430 spectra observed in RAVE survey. However, many spectra have to be combined (48 on average) in order to achieve a S/N ratio high enough to measure the DIB's profile, hence limiting the spatial information about the ISM. Only one strong interstellar line is included in the RAVE spectral range, a diffuse interstellar band at 8620 Å. We compare its equivalent width with extinction maps and with Bayesian reddening, calculated for individual stars, and provide a linear relation between the equivalent width and reddening. Separately from the introduced method, we calculate equivalent widths of the diffuse interstellar band in spectra of hot stars with known extinction and compare all three linear relations with each other and with relations from the literature.

preprint2013arXiv

In the thick of it: metal-poor disc stars in RAVE

By selecting in the RAVE-DR4 survey the stars located between 1 and 2 kpc above the Galactic plane, we question the consistency of the simplest three-component model (thin disc, thick disc, halo) for the Milky Way. We confirm that the metallicity and azimuthal velocity distribution functions of the thick disc are not Gaussian. In particular, we find that the thick disc has an extended metallicity tail going at least down to [M/H]=-2 dex, contributing roughly at 3% of the entire thick disc population and having a shorter scale-length compared to the canonical thick disc. The mean azimuthal velocity of these metal-poor stars allows us to estimate the correlation between the metallicity and the orbital velocity, which is an important constraint on the formation mechanisms of the Galactic thick disc. Given our simple approach, we find dVphi/d[M/H] ~ 50 km/s/dex, which is in a very good agreement with previous literature values. We complete the study with a brief discussion on the implications of the formation scenarios for the thick disc, and suggest that given the above mentioned characteristics, a thick disc mainly formed by radial migration mechanisms seems unlikely.

preprint2013arXiv

Kinematic groups beyond the Solar neighbourhood with RAVE

We analyse the kinematics of disc stars observed by the RAVE survey in and beyond the Solar neighbourhood.We detect significant overdensities in the velocity distributions using a technique based on the wavelet transform.We find that the main local kinematic groups are large scale features, surviving at least up to ~1 kpc from the Sun in the direction of anti-rotation, and also at ~700 pc below the Galactic plane.We also find that for regions located at different radii than the Sun, the known groups appear shifted in the velocity plane. For example, the Hercules group has a larger azimuthal velocity for regions inside the Solar circle and a lower value outside. We have also discovered a new group at (U, V) = (92,-22) km/s in the Solar neighbourhood and confirmed the significance of other previously found groups. Some of these trends detected for the first time are consistent with dynamical models of the effects of the bar and the spiral arms. More modelling is required to definitively characterise the non-axisymmetric components of our Galaxy using these groups.

preprint2013arXiv

Phase mixing due to the Galactic potential: steps in the position and velocity distributions of popped star clusters

As star clusters are expected to form with low star formation efficiencies, the gas in the cluster is expelled quickly and early in their development: the star cluster "pops." This leads to an unbound stellar system, evolving in the Galactic potential. Previous N-body simulations have demonstrated the existence of a stepped number density distribution of cluster stars after popping, both in vertical position and vertical velocity, with a passing resemblance to a Christmas tree. Using numerical and analytical methods, we investigate the source of this structure, which arises due to the phase mixing of the out-of-equilibrium stellar system, determined entirely by the background analytic potential. Considering only the vertical motions, we construct a theoretical model to describe the time evolution of the phase space distribution of stars in a Miyamoto-Nagai disk potential and a full Milky-Way type potential comprising bulge, halo and disk components, which is then compared with N-body simulations. Using our theoretical model, we investigate the possible observational signatures and the feasibility of detection.

preprint2013arXiv

The [alpha/Fe] ratios of very metal-poor stars within the IGIMF theory

The aim of this paper is to quantify the amplitude of the predicted plateau in [alpha/Fe] ratios associated with the most metal-poor stars of a galaxy. We assume that the initial mass function in galaxies is steeper if the star formation rate (SFR) is low -- as per the integrated galactic initial mass function (IGIMF) theory. A variant of the theory, in which the IGIMF depends upon the metallicity of the parent galaxy, is also considered. The IGIMF theory predicts low [alpha/Fe] plateaus in dwarf galaxies, characterised by small SFRs. The [alpha/Fe] plateau is up to 0.7dex lower than the corresponding plateau of the Milky Way. For a universal IMF one should expect instead that the [alpha/Fe] plateau is the same for all the galaxies, irrespective of their masses or SFRs. Assuming a strong dependence of the IMF on the metallicity of the parent galaxy, dwarf galaxies can show values of the [alpha/Fe] plateau similar to those of the Milky Way, and almost independent on the SFR. The [Mg/Fe] ratios of the most metal-poor stars in dwarf galaxies satellites of the Milky Way can be reproduced either if we consider metallicity-dependent IMFs or if the early SFRs of these galaxies were larger than we presently think. Present and future observations of dwarf galaxies can help disentangle between these different IGIMF formulations.

preprint2013arXiv

The asymmetric drift, the local standard of rest, and implications from RAVE data

The determination of the LSR is still a matter of debate. The classical value of the tangential peculiar motion of the Sun with respect to the LSR was challenged in recent years, claiming a significantly larger value. We show that the RAdial Velocity Experiment (RAVE) sample of dwarf stars is an excellent data set to derive tighter boundary conditions to chemodynamical evolution models of the extended solar neighbourhood. We present an improved Jeans analysis, which allows a better interpretation of the measured kinematics of stellar populations in the Milky Way disc. We propose an improved version of the Strömberg relation with the radial scalelengths as the only unknown. Binning RAVE stars in metallicity reveals a bigger asymmetric drift (corresponding to a smaller radial scalelength) for more metal-rich populations. With the standard assumption of velocity-dispersion independent radial scalelengths in each metallicity bin, we redetermine the LSR. The new Strömberg equation yields a joint LSR value of V_\sun=3.06 \pm 0.68 km/s, which is even smaller than the classical value based on Hipparcos data. The corresponding radial scalelength increases from 1.6 kpc for the metal-rich bin to 2.9 kpc for the metal-poor bin, with a trend of an even larger scalelength for young metal-poor stars. When adopting the recent Schönrich value of V_\sun=12.24 km/s for the LSR, the new Strömberg equation yields much larger individual radial scalelengths of the RAVE subpopulations, which seem unphysical in part. The new Strömberg equation allows a cleaner interpretation of the kinematic data of disc stars in terms of radial scalelengths. Lifting the LSR value by a few km/s compared to the classical value results in strongly increased radial scalelengths with a trend of smaller values for larger velocity dispersions.

preprint2013arXiv

The environmental dependence of neutral hydrogen in the GIMIC simulations

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) cosmological hydrodynamic simulation at z=0 to study the distribution and environmental dependence of neutral hydrogen (HI) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly$α$ absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse HI in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a re-simulation of the average cosmic density GIMIC region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average HI radial profiles by stacking the individual profiles according to both mass and environment. We find high HI column densities at large impact parameters in group environments and markedly lower HI densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments.

preprint2013arXiv

The Influence of Ram Pressure on the Evolution of Tidal Dwarf Galaxies

The formation mechanism of tidal dwarf galaxies means they are expected to contain little or no dark matter. As such, they might be expected to be very sensitive to their environment. We investigate the impact of ram pressure on tidal dwarf galaxies in a parameter study, varying dwarf galaxy properties and ram pressures. We submit model tidal dwarf galaxies to wind-tunnel style tests using a toy ram pressure model. The effects of ram pressure are found to be substantial. If tidal dwarf galaxies have their gas stripped, they may be completely destroyed. Ram pressure drag causes acceleration of our dwarf galaxy models, and this further enhances stellar losses. The dragging can also cause stars to lie in a low surface brightness stellar stream that points in the opposite direction to the stripped gas, in a manner distinctive from tidal streams. We investigate the effects of ram pressure on surface density profiles, the dynamics of the stars, and discuss the consequences for dynamical mass measurements.

preprint2013arXiv

The RAdial Velocity Experiment (RAVE): Fourth data release

We present the stellar atmospheric parameters (effective temperature, surface gravity, overall metallicity), radial velocities, individual abundances and distances determined for 425 561 stars, which constitute the fourth public data release of the RAdial Velocity Experiment (RAVE). The stellar atmospheric parameters are computed using a new pipeline, based on the algorithms of MATISSE and DEGAS. The spectral degeneracies and the 2MASS photometric information are now better taken into consideration, improving the parameter determination compared to the previous RAVE data releases. The individual abundances for six elements (magnesium, aluminum, silicon, titanium, iron and nickel) are also given, based on a special-purpose pipeline which is also improved compared to that available for the RAVE DR3 and Chemical DR1 data releases. Together with photometric information and proper motions, these data can be retrieved from the RAVE collaboration website and the Vizier database.

preprint2013arXiv

The relation between chemical abundances and kinematics of the Galactic disc with RAVE

Aims: We study the relations between stellar kinematics and chemical abundances of a large sample of RAVE giants in search for selection criteria needed for disentangling different Galactic stellar populations. Methods: We select a sample of 2167 giant stars with signal-to-noise per spectral measurements above 75 from the RAVE chemical catalogue and follow the analysis performed by Gratton and colleagues on 150 subdwarf stars spectroscopically observed at high-resolution. We then use a larger sample of 9131 giants (with signal-to-noise above 60) to investigate the chemo-kinematical characteristics of our stars by grouping them into nine subsamples with common eccentricity ($e$) and maximum distance achieved above the Galactic plane ($Z_max$). Results: The RAVE kinematical and chemical data proved to be reliable by reproducing the results by Gratton et al. obtained with high-resolution spectroscopic data. Our analysis, based on the $e$-$Z_max$ plane combined with additional orbital parameters and chemical information, provides an alternative way of identifying different populations of stars. In addition to extracting canonical thick- and thin-disc samples, we find a group of stars in the Galactic plane ($Z_max<1$ kpc and 0.4 $< e < $0.6), which show homogeneous kinematics but differ in their chemical properties. We interpret this as a clear sign that some of these stars have experienced the effects of heating and/or radial migration, which have modified their original orbits. The accretion origin of such stars cannot be excluded.

preprint2013arXiv

The wobbly Galaxy: kinematics north and south with RAVE red clump giants

The RAVE survey, combined with proper motions and distance estimates, can be used to study in detail stellar kinematics in the extended solar neighbourhood (solar suburb). Using the red clump, we examine the mean velocity components in 3D between an R of 6 and 10 kpc and a Z of -2 to 2 kpc, concentrating on North-South differences. Simple parametric fits to the R, Z trends for VPHI and the velocity dispersions are presented. We confirm the recently discovered gradient in mean Galactocentric radial velocity, VR, finding that the gradient is more marked below the plane, with a Z gradient also present. The vertical velocity, VZ, also shows clear structure, with indications of a rarefaction-compression pattern, suggestive of wave-like behaviour. We perform a rigorous error analysis, tracing sources of both systematic and random errors. We confirm the North-South differences in VR and VZ along the line-of-sight, with the VR estimated independent of the proper motions. The complex three-dimensional structure of velocity space presents challenges for future modelling of the Galactic disk, with the Galactic bar, spiral arms and excitation of wave-like structures all probably playing a role.

preprint2012arXiv

Exploring the Morphology of RAVE Stellar Spectra

The RAdial Velocity Experiment (RAVE) is a medium resolution R~7500 spectroscopic survey of the Milky Way which already obtained over half a million stellar spectra. They present a randomly selected magnitude-limited sample, so it is important to use a reliable and automated classification scheme which identifies normal single stars and discovers different types of peculiar stars. To this end we present a morphological classification of 350,000 RAVE survey stellar spectra using locally linear embedding, a dimensionality reduction method which enables representing the complex spectral morphology in a low dimensional projected space while still preserving the properties of the local neighborhoods of spectra. We find that the majority of all spectra in the database ~90-95% belong to normal single stars, but there is also a significant population of several types of peculiars. Among them the most populated groups are those of various types of spectroscopic binary and chromospherically active stars. Both of them include several thousands of spectra. Particularly the latter group offers significant further investigation opportunities since activity of stars is a known proxy of stellar ages. Applying the same classification procedure to the sample of normal single stars alone shows that the shape of the projected manifold in two dimensional space correlates with stellar temperature, surface gravity and metallicity.

preprint2012arXiv

Gas and Metal Distributions within Simulated Disk Galaxies

We highlight two research strands related to our ongoing chemodynamical Galactic Archaeology efforts: (i) the spatio-temporal infall rate of gas onto the disk, drawing analogies with the infall behaviour imposed by classical galactic chemical evolution models of inside-out disk growth; (ii) the radial age gradient predicted by spectrophometric models of disk galaxies. In relation to (i), at low-redshift, we find that half of the infall onto the disk is gas associated with the corona, while half can be associated with cooler gas streams; we also find that gas enters the disk preferentially orthogonal to the system, rather than in-plane. In relation to (ii), we recover age gradient troughs/inflections consistent with those observed in nature, without recourse to radial migrations.

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

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

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

preprint2012arXiv

Metallicity Gradients in Simulated Disk Galaxies

The stellar metallicity and abundance ratio gradients from the fiducial late-type galaxy simulation of Stinson et al. (2010) are presented. Over 1-3 scalelengths, gradients are shown to flatten with time, consistent with empirical evidence at high- and low-redshifts. Kinematic effects, including radial migration, though, flatten these intrinsicly steep gradients such that by redshift z=0, the measured gradients of these (now) old stars are flatter than their young counterparts, in contradiction to what is observed locally. Conversely, the stellar [O/Fe] gradient is (to first order) robust against migration, remaining fairly flat for both young and old populations today.

preprint2012arXiv

Properties of simulated Milky Way-mass galaxies in loose group and field environments

We test the validity of comparing simulated field disk galaxies with the empirical properties of systems situated within environments more comparable to loose groups, including the Milky Way's Local Group. Cosmological simulations of Milky Way-mass galaxies have been realised in two different environment samples: in the field and in environments with similar properties to the Local Group. Apart from the environments of the galaxies, the samples are kept as homogeneous as possible with equivalent ranges in last major merger time, halo mass and halo spin. Comparison of these two samples allow for systematic differences in the simulations to be identified. Metallicity gradients, disk scale lengths, colours, magnitudes and age-velocity dispersion relations are studied for each galaxy in the suite and the strength of the link between these and environment of the galaxies is studied. The bulge-to-disk ratio of the galaxies show that these galaxies are less spheroid dominated than many other simulated galaxies in literature with the majority of both samples being disk dominated. We find that secular evolution and mergers dominate the spread of morphologies and metallicity gradients with no visible differences between the two environment samples. In contrast with this consistency in the two samples there is tentative evidence for a systematic difference in the velocity dispersion-age relations of galaxies in the different environments. Loose group galaxies appear to have more discrete steps in their velocity dispersion-age relations. We conclude that at the current resolution of cosmological galaxy simulations field environment galaxies are sufficiently similar to those in loose groups to be acceptable proxies for comparison with the Milky Way provided that a similar assembly history is considered.

preprint2012arXiv

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

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

preprint2012arXiv

The Distribution of Young Stars and Metals in Simulated Cosmological Disk Galaxies

We examine the distribution of young stars associated with the spiral arms of a simulated L* cosmological disk galaxy. We find age patterns orthogonal to the arms which are not inconsistent with the predictions of classical density wave theory, a view further supported by recent observations of face-on Grand Design spirals such as M51. The distribution of metals within a simulated ~0.1L* disk is presented, reinforcing the link between star formation, the age-metallicity relation, and the metallicity distribution function.

preprint2012arXiv

The Halo Shape and Evolution of Polar Disc Galaxies

We examine the properties and evolution of a simulated polar disc galaxy. This galaxy is comprised of two orthogonal discs, one of which contains old stars (old stellar disc), and the other, containing both younger stars and the cold gas (polar disc) of the galaxy. By exploring the shape of the inner region of the dark matter halo, we are able to confirm that the halo shape is a oblate ellipsoid flattened in the direction of the polar disc. We also note that there is a twist in the shape profile, where the innermost 3 kpc of the halo flattens in the direction perpendicular to the old disc, and then aligns with the polar disc out until the virial radius. This result is then compared to the halo shape inferred from the circular velocities of the two discs. We also use the temporal information of the simulation to track the system's evolution, and identify the processes which give rise to this unusual galaxy type. We confirm the proposal that the polar disc galaxy is the result of the last major merger, where the angular moment of the interaction is orthogonal to the angle of the infalling gas. This merger is followed by the resumption of coherent gas infall. We emphasise that the disc is rapidly restored after the major merger and that after this event the galaxy begins to tilt. A significant proportion of the infalling gas comes from filaments. This infalling gas from the filament gives the gas its angular momentum, and, in the case of the polar disc galaxy, the direction of the gas filament does not change before or after the last major merger.

preprint2012arXiv

The properties of the local spiral arms from RAVE data: two-dimensional density wave approach

Using the RAVE survey, we recently brought to light a gradient in the mean galactocentric radial velocity of stars in the extended solar neighbourhood. This gradient likely originates from non-axisymmetric perturbations of the potential, among which a perturbation by spiral arms is a possible explanation. Here, we apply the traditional density wave theory and analytically model the radial component of the two-dimensional velocity field. Provided that the radial velocity gradient is caused by relatively long-lived spiral arms that can affect stars substantially above the plane, this analytic model provides new independent estimates for the parameters of the Milky Way spiral structure. Our analysis favours a two-armed perturbation with the Sun close to the inner ultra-harmonic 4:1 resonance, with a pattern speed Ω_p=18.6^{+0.3}_{-0.2} km/s/kpc and a small amplitude A=0.55 \pm 0.02% of the background potential (14% of the background density). This model can serve as a basis for numerical simulations in three dimensions, additionally including a possible influence of the galactic bar and/or other non-axisymmetric modes.

preprint2012arXiv

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

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

preprint2012arXiv

The VMC Survey - VI. Quasars behind the Magellanic system

The number and spatial distribution of confirmed quasi-stellar objects (QSOs) behind the Magellanic system is limited. This undermines their use as astrometric reference objects for different types of studies. We have searched for criteria to identify candidate QSOs using observations from the VISTA survey of the Magellanic Clouds system (VMC) that provides photometry in the YJKs bands and 12 epochs in the Ks band. The (Y-J) versus (J-Ks) diagram has been used to distinguish QSO candidates from Milky Way stars and stars of the Magellanic Clouds. Then, the slope of variation in the Ks band has been used to identify a sample of high confidence candidates. These criteria were developed based on the properties of 117 known QSOs presently observed by the VMC survey. VMC YJKs magnitudes and Ks light-curves of known QSOs behind the Magellanic system are presented. About 75% of them show a slope of variation in Ks>10^-4 mag/day and the shape of the light-curve is in general irregular and without any clear periodicity. The number of QSO candidates found in tiles including the South Ecliptic Pole and the 30 Doradus regions is 22 and 26, respectively, with a ~20% contamination by young stellar objects, planetary nebulae, stars and normal galaxies. By extrapolating the number of QSO candidates to the entire VMC survey area we expect to find about 1200 QSOs behind the LMC, 400 behind the SMC, 200 behind the Bridge and 30 behind the Stream areas, but not all will be suitable for astrometry. Further, the Ks band light-curves can help support investigations of the mechanism responsible for the variations.

preprint2012arXiv

The VMC Survey. V. First results for Classical Cepheids

The VISTA Magellanic Cloud (VMC, PI M.R. Cioni) survey is collecting deep Ks-band time-series photometry of the pulsating variable stars hosted by the system formed by the two Magellanic Clouds (MCs) and the "bridge" connecting them. In this paper we present the first results for Classical Cepheids, from the VMC observations of two fields in the Large Magellanic Cloud (LMC). The VMC Ks-band light curves of the Cepheids are well sampled (12-epochs) and of excellent precision. We were able to measure for the first time the Ks magnitude of the faintest Classical Cepheids in the LMC (Ks\sim17.5 mag), which are mostly pulsating in the First Overtone (FO) mode, and to obtain FO Period-Luminosity (PL), Period-Wesenheit (PW), and Period-Luminosity-Color (PLC) relations, spanning the full period range from 0.25 to 6 day. Since the longest period Cepheid in our dataset has a variability period of 23 day, we have complemented our sample with literature data for brighter F Cepheids. On this basis we have built a PL relation in the Ks band that, for the first time, includes short period pulsators, and spans the whole range from 1.6 to 100 days in period. We also provide the first ever empirical PW and PLC relations using the (V-Ks) color and time-series Ks photometry. The very small dispersion (\sim0.07 mag) of these relations makes them very well suited to study the three-dimensional (3D) geometry of the Magellanic system. The use of "direct" (parallax- and Baade-Wesselink- based) distance measurements to both Galactic and LMC Cepheids, allowed us to calibrate the zero points of the PL, PW, and PLC relations obtained in this paper, and in turn to estimate an absolute distance modulus of (m-M)0=18.46\pm0.03 for the LMC. This result is in agreement with most of the latest literature determinations based on Classical Cepheids.

preprint2012arXiv

Thick disk kinematics from RAVE and the solar motion

Radial velocity surveys such as the Radial Velocity Experiment (RAVE) provide us with measurements of hundreds of thousands of nearby stars most of which belong to the Galactic thin, thick disk or halo. Ideally, to study the Galactic disks (both thin and thick) one should make use of the multi-dimensional phase-space and the whole pattern of chemical abundances of their stellar populations. In this paper, with the aid of the RAVE Survey, we study the thin and thick disks of the Milky Way, focusing on the latter. We present a technique to disentangle the stellar content of the two disks based on the kinematics and other stellar parameters such as the surface gravity of the stars. Using the Padova Galaxy Model, we checked the ability of our method to correctly isolate the thick disk component from the Galaxy mixture of stellar populations. We introduce selection criteria in order to clean the observed radial velocities from the Galactic differential rotation and to take into account the partial sky coverage of RAVE. We developed a numerical technique to statistically disentangle thin and thick disks from their mixture. We deduce the components of the solar motion relative to the Local Standard of Rest (LSR) in the radial and vertical direction, the rotational lag of the thick disk component relative to the LSR, and the square root of the absolute value of the velocity dispersion tensor for the thick disk alone. The analysis of the thin disk is presented in another paper. We find good agreement with previous independent parameter determinations. In our analysis we used photometrically determined distances. In the Appendix we show that similar values can be found for the thick disk alone as derived in the main sections of our paper even without the knowledge of photometric distances.

preprint2012arXiv

Thin disc, Thick Disc and Halo in a Simulated Galaxy

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

preprint2012arXiv

Thin disk kinematics from RAVE and the solar motion

Aims. We study the Milky Way thin disk with the Radial Velocity Experiment (RAVE) survey. We consider the thin and thick disks as different Galactic components and present a technique to statistically disentangle the two populations. Then we focus our attention on the thin disk component. Methods. We disentangle the thin disk component from amixture of the thin and thick disks using a data set providing radial velocities, proper motions, and photometrically determined distances. Results. We present the trend of the velocity dispersions in the thin disk component of the Milky Way (MW) in the radial direction above and below the Galactic plane using data from the RAdial Velocity Experiment (RAVE). The selected sample is a limited subsample from the entire RAVE catalogue, roughly mapping up to 500 pc above and below the Galactic plane, a few degrees in azimuthal direction and covering a radial extension of 2.0 kpc around the solar position. The solar motion relative to the local standard of rest is also re-determined with the isolated thin disk component. Major results are the trend of the velocity mean and dispersion in the radial and vertical direction. In addition the azimuthal components of the solar motion relative to the local standard of rest and the velocity dispersion are discussed.

preprint2011arXiv

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

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

preprint2011arXiv

An HI Survey of Six Local Group Analogs. II. HI properties of group galaxies

We have conducted an HI 21 cm emission-line survey of six loose groups of galaxies chosen to be analogs to the Local Group. The survey was conducted using the Parkes Multibeam instrument and the Australia Telescope Compact Array (ATCA) over a ~1 Mpc^2 area and covering the full depth of each group, with a M(HI) sensitivity of ~7x10^5 M(sun). Our survey detected 110 sources, 61 of which are associated with the six groups. All of these sources were confirmed with ATCA observations or were previously cataloged by HIPASS. The sources all have optical counterparts and properties consistent with dwarf irregular or late-type spiral galaxies. We present here the HI properties of the groups and their galaxies. We derive an HI mass function for the groups that is consistent with being flatter than the equivalent field HIMF. We also derive a circular velocity distribution function, tracing the luminous dark matter halos in the groups, that is consistent with those of the Local Group and HIPASS galaxies, both of which are shallower than that of clusters or predictions from CDM models of galaxy formation.

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

The Chemical and Dynamical Evolution of Isolated Dwarf Galaxies

Using a suite of simulations (Governato et al. 2010) which successfully produce bulgeless (dwarf) disk galaxies, we provide an analysis of their associated cold interstellar media (ISM) and stellar chemical abundance patterns. A preliminary comparison with observations is undertaken, in order to assess whether the properties of the cold gas and chemistry of the stellar components are recovered successfully. To this end, we have extracted the radial and vertical gas density profiles, neutral hydrogen velocity dispersion, and the power spectrum of structure within the ISM. We complement this analysis of the cold gas with a brief examination of the simulations' metallicity distribution functions and the distribution of alpha-elements-to-iron.

preprint2011arXiv

The Cold Gas Content of Bulgeless Dwarf Galaxies

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

preprint2011arXiv

The RAdial Velocity Experiment (RAVE): Third Data Release

We present the third data release of the RAdial Velocity Experiment (RAVE) which is the first milestone of the RAVE project, releasing the full pilot survey. The catalog contains 83,072 radial velocity measurements for 77,461 stars in the southern celestial hemisphere, as well as stellar parameters for 39,833 stars. This paper describes the content of the new release, the new processing pipeline, as well as an updated calibration for the metallicity based upon the observation of additional standard stars. Spectra will be made available in a future release. The data release can be accessed via the RAVE webpage: http://www.rave-survey.org.

preprint2011arXiv

The RAVE catalogue of stellar elemental abundances: first data release

We present chemical elemental abundances for $36,561$ stars observed by the RAdial Velocity Experiment (RAVE), an ambitious spectroscopic survey of our Galaxy at Galactic latitudes $|$b$|>25^{\circ}$ and with magnitudes in the range 9$<I_{DENIS}<$13. RAVE spectra cover the Ca-triplet region at 8410--8795Å with resolving power R$\sim$7500. This first data release of the RAVE chemical catalogue is complementary to the third RAVE data release of radial velocities and stellar parameters, and it contains chemical abundances for the elements Mg, Al, Si, Ca, Ti, Fe and Ni, with a mean error of $\sim$0.2 dex, as judged from accuracy tests performed on synthetic and real spectra. Abundances are estimated through a dedicated processing pipeline in which the curve of growth of individual lines is obtained from a library of absorption-line equivalent widths to construct a model spectrum that is then matched to the observed spectrum via a $χ^2$-minimization technique. We plan to extend this pipeline to include estimates for other elements, such as oxygen and sulfur, in future data releases.

preprint2010arXiv

A search for new members of the beta Pic, Tuc-Hor and epsilon Cha moving groups in the RAVE database

We report on the discovery of new members of nearby young moving groups, exploiting the full power of combining the RAVE survey with several stellar age diagnostic methods and follow-up high-resolution optical spectroscopy. The results include the identification of one new and five likely members of the beta Pictoris moving group, ranging from spectral types F9 to M4 with the majority being M dwarfs, one K7 likely member of the epsilon Cha group and two stars in the Tuc-Hor association. Based on the positive identifications we foreshadow a great potential of the RAVE database in progressing toward a full census of young moving groups in the solar neighbourhood.

preprint2010arXiv

Detection of a radial velocity gradient in the extended local disc with RAVE

Using a sample of 213,713 stars from the Radial Velocity Experiment (RAVE) survey, limited to a distance of 2 kpc from the Sun and to |z|<1 kpc, we report the detection of a velocity gradient of disc stars in the fourth quadrant, directed radially from the Galactic centre. In the direction of the Galactic centre, we apply a simple method independent of stellar proper motions and of Galactic parameters to assess the existence of this gradient in the RAVE data. This velocity gradient corresponds to |K+C| < 3 km/s/kpc, where K and C are the Oort constants measuring the local divergence and radial shear of the velocity field, respectively. In order to illustrate the effect, assuming a zero radial velocity of the Local Standard of Rest we then reconstruct the two-dimensional Galactocentric velocity maps using two different sets of proper motions and photometric distances based either on isochrone fitting or on K-band magnitudes, and considering two sets of values for the Galactocentric radius of the Sun and local circular speed. Further observational confirmation of our finding with line-of-sight velocities of stars at low latitudes, together with further modelling, should help constrain the non-axisymmetric components of the Galactic potential, including the bar, the spiral arms and possibly the ellipticity of the dark halo.

preprint2010arXiv

Origins of the Thick Disk as Traced by the Alpha-Elements of Metal-Poor Giant Stars Selected from RAVE

Theories of thick disk formation can be differentiated by measurements of stellar elemental abundances. We have undertaken a study of metal-poor stars selected from the RAVE spectroscopic survey of bright stars to establish whether or not there is a significant population of metal-poor thick-disk stars ([Fe/H] <~ -1.0) and to measure their elemental abundances. In this paper, we present abundances of four alpha-elements (Mg, Si, Ca, Ti) and iron for a subsample of 212 RGB and 31 RC/HB stars from this study. We find that the [alpha/Fe] ratios are enhanced implying that enrichment proceeded by purely core-collapse supernovae. This requires that star formation in each star forming region had a short duration. The relative lack of scatter in the [alpha/Fe] ratios implies good mixing in the ISM prior to star formation. In addition, the ratios resemble that of the halo, indicating that the halo and thick disk share a similar massive star IMF. We conclude that the alpha enhancement of the metal-poor thick disk implies that direct accretion of stars from dwarf galaxies similar to surviving dwarf galaxies today did not play a major role in the formation of the thick disk.

preprint2010arXiv

Testing formation mechanisms of the Milky Way's thick disc with RAVE

We study the eccentricity distribution of a thick disc sample of stars observed in the Radial Velocity Experiment (RAVE) and compare it to that expected in four simulations of thick disc formation in the literature (accretion of satellites, heating of a primordial thin disc during a merger, radial migration, and gas-rich mergers), as compiled by Sales et al. (2009). We find that the distribution of our sample is peaked at low eccentricities and falls off smoothly and rather steeply to high eccentricities. This distribution is fairly robust to changes in distances, thin disc contamination, and the particular thick disc sample used. Our results are inconsistent with what is expected for the pure accretion simulation, since we find that the dynamics of local thick disc stars implies that the majority must have formed "in situ". Of the remaining models explored, the eccentricity distribution of our stars appears to be most consistent with the gas-rich merger case.

preprint2010arXiv

The RAVE Survey: Rich in Very Metal-Poor Stars

Very metal-poor stars are of obvious importance for many problems in chemical evolution, star formation, and galaxy evolution. Finding complete samples of such stars which are also bright enough to allow high-precision individual analyses is of considerable interest. We demonstrate here that stars with iron abundances [Fe/H] < -2 dex, and down to below -4 dex, can be efficiently identified within the Radial Velocity Experiment (RAVE) survey of bright stars, without requiring additional confirmatory observations. We determine a calibration of the equivalent width of the Calcium triplet lines measured from the RAVE spectra onto true [Fe/H], using high spectral resolution data for a subset of the stars. These RAVE iron abundances are accurate enough to obviate the need for confirmatory higher-resolution spectroscopy. Our initial study has identified 631 stars with [Fe/H] <= -2, from a RAVE database containing approximately 200,000 stars. This RAVE-based sample is complete for stars with [Fe/H] < -2.5, allowing statistical sample analysis. We identify three stars with [Fe/H] <= -4. Of these, one was already known to be `ultra metal-poor', one is a known carbon-enhanced metal-poor star, but we obtain [Fe/H]= -4.0, rather than the published [Fe/H]=-3.3, and derive [C/Fe] = +0.9, and [N/Fe] = +3.2, and the third is at the limit of our S/N. RAVE observations are on-going and should prove to be a rich source of bright, easily studied, very metal-poor stars.

preprint2009arXiv

GASS: The Parkes Galactic All-Sky Survey. I. Survey Description, Goals, and Initial Data Release

The Parkes Galactic All-Sky Survey (GASS) is a survey of Galactic atomic hydrogen (HI) emission in the Southern sky covering declinations $δ\leq 1^{\circ}$ using the Parkes Radio Telescope. The survey covers $2π$ steradians with an effective angular resolution of ~16', at a velocity resolution of 1.0 km/s, and with an rms brightness temperature noise of 57 mK. GASS is the most sensitive, highest angular resolution survey of Galactic HI emission ever made in the Southern sky. In this paper we outline the survey goals, describe the observations and data analysis, and present the first-stage data release. The data product is a single cube at full resolution, not corrected for stray radiation. Spectra from the survey and other data products are publicly available online.

preprint2009arXiv

RAVE spectroscopy of luminous blue variables in the Large Magellanic Cloud

CONTEXT: The RAVE spectroscopic survey for galactic structure and evolution obtains 8400-8800 Ang spectra at 7500 resolving power at the UK Schmidt Telescope using the 6dF multi-fiber positioner. More than 300,000 9<I<12 and |b|>25 deg southern stars have been observed to date. AIMS: This paper presents the first intrinsic examination of stellar spectra from the RAVE survey, aimed at evaluating their diagnostic potential for peculiar stars and at contributing to the general understanding of Luminous Blue Variables (LBVs). METHODS: We used the multi-epoch spectra for all seven LBVs observed, between 2005 and 2008, in the Large Magellanic Cloud (LMC) by the RAVE survey. RESULTS: We demonstrate that RAVE spectra possess significant diagnostic potential when applied to peculiar stars and, in particular, LBVs. The behaviour of the radial velocities for both emission and absorption lines, and the spectral changes between outburst and quiescence states are described and found to agree with evidence gathered at more conventional wavelengths. The wind outflow signatures and their variability are investigated, with multi-components detected in S Doradus. Photoionisation modelling of the rich emission line spectrum of R 127 shows evidence of a massive detached ionised shell that was ejected during the 1982-2000 outburst. Surface inhomogeneities in the nuclear-processed material, brought to the surface by heavy mass loss, could have been observed in S Doradus, even if alternative explanations are possible. We also detect the transition from quiescence to outburst state in R 71. Finally, our spectrum of R 84 offers one of the clearest views of its cool companion.

preprint2008arXiv

The Magellanic Clouds as a template for the study of stellar populations and galaxy interactions

The Magellanic System represents one of the best places to study the formation and evolution of galaxies. Photometric surveys of various depths, areas and wavelengths have had a significant impact on our understanding of the system; however, a complete picture is still lacking. VMC (the VISTA near-infrared YJKs survey of the Magellanic System) will provide new data to derive the spatially resolved star formation history and to construct a three-dimensional map of the system. These data combined with those from other ongoing and planned surveys will give us an absolutely unique view of the system opening up the doors to truly new science!

preprint2007arXiv

An HI survey of six Local Group analogs: I. Survey description and the search for high-velocity clouds

We have conducted an HI 21 cm emission-line survey using the Parkes 20cm multibeam instrument and the Australia Telescope Compact Array (ATCA) of six loose groups of galaxies chosen to be analogs to the Local Group. The goal of this survey is to make a census of the HI-rich galaxies and high-velocity clouds (HVCs) within these groups and compare these populations with those in the Local Group. The Parkes observations covered the entire volume of each group with a rms M(HI) sensitivity of 4-10x10^5 M(sun) per 3.3 km/s channel. All potential sources detected in the Parkes data were confirmed with ATCA observations at ~2' resolution and the same M(sun) sensitivity. All the confirmed sources have associated stellar counterparts; no starless HI clouds--HVC analogs--were found in the six groups. In this paper, we present a description of the survey parameters, its sensitivity and completeness. Using the population of compact HVCs (CHVCs) around the Milky Way as a template coupled with the detailed knowledge of our survey parameters, we infer that our non-detection of CHVC analogs implies that, if similar populations exist in the six groups studied, the CHVCs must be clustered within 90 kpc of group galaxies, with average M(HI) < 4x10^5 M(sun) at the 95% confidence level. The corollary is that the same must apply to Milky Way CHVCs. This is consistent with our previous results from a smaller sample of groups, and in accordance with recent observational and theoretical constraints from other authors. These results confirm that there is very little neutral matter around galaxies, and that any substantial reservoir of baryons must be in other phases.

preprint2005arXiv

Stellar abundance gradients in galactic disks. I. Method and spectral line gradients

We describe the technique of absorption line imaging of galaxy disks using the Taurus Tunable Filter on the Anglo-Australian Telescope and demonstrate its sensitivity to the behaviour of spectral features associated with Mg and Fe. Radial profiles of Mg2 and Fe5270 line-strengths are presented for a sample of eight face-on spiral galaxies spanning a range of Hubble types. Signatures of phenomena including merger-induced star formation, HII rings and galactic bars are also reported. This study demonstrates the capacity of tunable filters to measure Mg and Fe line-strengths across the face of spiral galaxies, which can ultimately reveal clues about the star formation history and chemical evolution.

preprint2005arXiv

Studying Galaxy Formation in Loose Galaxy Groups

We present the results of our HI survey of six loose groups of galaxies analogous to the Local Group. The survey was conducted using the Parkes telescope and the Australia Telescope Compact Array to produce a census of all the gas-rich galaxies and analogs to the high-velocity clouds (HVCs) within these groups down to M(HI) < 10^7 M(sun) as a test of models of galaxy formation. We present the HI mass function and halo mass function of the loose groups and show that they are consistent with those of the Local Group. We discuss the possible role of HVCs in solving the ``missing satellite'' problem and discuss the implications of our observations for models of galaxy formation.

preprint2004arXiv

Globular Cluster and Galaxy Formation: M31, the Milky Way and Implications for Globular Cluster Systems of Spiral Galaxies

The globular cluster (GC) systems of the Milky Way and of our neighboring spiral galaxy, M31, comprise 2 distinct entities, differing in 3 respects. 1. M31 has young GCs, ages from ~100 Myr to 5 Gyr old, as well as old globular clusters. No such young GCs are known in the Milky Way. 2. We confirm that the oldest M31 GCs have much higher nitrogen abundances than do Galactic GCs at equivalent metallicities. 3. Morrison et al. found M31 has a subcomponent of GCs that follow closely the disk rotation curve of M31. Such a GC system in our own Galaxy has yet to be found. These data are interpreted in terms of the hierarchical-clustering-merging (HCM) paradigm for galaxy formation. We infer that M31 has absorbed more of its dwarf systems than has the Milky Way. This inference has 3 implications: 1. All spiral galaxies likely differ in their GC properties, depending on how many companions each galaxy has, and when the parent galaxy absorbs them. The the Milky Way ties down one end of this spectrum, as almost all of its GCs were absorbed 10-12 Gyr ago. 2. It suggests that young GCs are preferentially formed in the dwarf companions of parent galaxies, and then absorbed by the parent galaxy during mergers. 3. Young GCs seen in tidally-interacting galaxies might come from dwarf companions of these galaxies, rather than be made a-new in the tidal interaction. There is no ready explanation for the marked difference in nitrogen abundance for old M31 GCs relative to the oldest Galactic GCs. The predictions made by Li & Burstein regarding the origin of nitrogen abundance in globular clusters are consistent with what is found for the old M31 GCs compared to that for the two 5 Gyr-old M31 GCs.

preprint2003arXiv

The 1000 Brightest HIPASS Galaxies: HI Mass Function and Omega_HI

We present a new accurate measurement of the HI mass function of galaxies from the HIPASS Bright Galaxy Catalog, a sample of 1000 galaxies with the highest HI peak flux densities in the southern hemisphere (Koribalski et al. 2003). This sample spans nearly four orders of magnitude in HI mass (from log M_HI/M_sun=6.8 to 10.6, H0=75) and is the largest sample of HI selected galaxies to date. We develop a bivariate maximum likelihood technique to measure the space density of galaxies, and show that this is a robust method, insensitive to the effects of large scale structure. The resulting HI mass function can be fitted satisfactorily with a Schechter function with faint-end slope alpha=-1.30. This slope is found to be dependent on morphological type, with later type galaxies giving steeper slopes. We extensively test various effects that potentially bias the determination of the HI mass function, including peculiar motions of galaxies, large scale structure, selection bias, and inclination effects, and quantify these biases. The large sample of galaxies enables an accurate measurement of the cosmological mass density of neutral gas: Omega_HI=(3.8 +/- 0.6) x 10^{-4}. Low surface brightness galaxies contribute only 15% to this value, consistent with previous findings.

preprint2000arXiv

Final Results from the Hubble Space Telescope Key Project to Measure the Hubble Constant

We present here the final results of the Hubble Space Telescope Key Project to measure the Hubble constant. We summarize our method, the results and the uncertainties, tabulate our revised distances, and give the implications of these results for cosmology. The analysis presented here benefits from a number of recent improvements and refinements, including (1) a larger LMC Cepheid sample to define the fiducial period-luminosity (PL) relations, (2) a more recent HST Wide Field and Planetary Camera 2 (WFPC2) photometric calibration, (3) a correction for Cepheid metallicity, and (4) a correction for incompleteness bias in the observed Cepheid PL samples. New, revised distances are given for the 18 spiral galaxies for which Cepheids have been discovered as part of the Key Project, as well as for 13 additional galaxies with published Cepheid data. The new calibration results in a Cepheid distance to NGC 4258 in better agreement with the maser distance to this galaxy. Based on these revised Cepheid distances, we find values (in km/sec/Mpc) of H0 = 71 +/- 2 (random) +/- 6 (systematic) (type Ia supernovae), 71 +/- 2 +/- 7 (Tully-Fisher relation), 70 +/- 5 +/- 6 (surface brightness fluctuations), 72 +/- 9 +/- 7 (type II supernovae), and 82 +/- 6 +/- 9 (fundamental plane). We combine these results for the different methods with 3 different weighting schemes, and find good agreement and consistency with H0 = 72 +/- 8. Finally, we compare these results with other, global methods for measuring the Hubble constant.

preprint1998arXiv

The disruption of nearby galaxies by the Milky Way

Interactions between galaxies are common and are an important factor in determining their physical properties such as position along the Hubble sequence and star-formation rate. There are many possible galaxy interaction mechanisms, including merging, ram-pressure stripping, gas compression, gravitational interaction and cluster tides. The relative importance of these mechanisms is often not clear, as their strength depends on poorly known parameters such as the density, extent and nature of the massive dark halos that surround galaxies. A nearby example of a galaxy interaction where the mechanism is controversial is that between our own Galaxy and two of its neighbours -- the Large and Small Magellanic Clouds. Here we present the first results of a new HI survey which provides a spectacular view of this interaction. In addition to the previously known Magellanic Stream, which trails 100 degrees behind the Clouds, the new data reveal a counter-stream which lies in the opposite direction and leads the motion of the Clouds. This result supports the gravitational model in which leading and trailing streams are tidally torn from the body of the Magellanic Clouds.

preprint1998arXiv

The HST Key Project on the Extragalactic Distance Scale XIV. The Cepheids in NGC 1365

We report the detection of Cepheid variable stars in the barred spiral galaxy NGC 1365, located in the Fornax cluster, using the Hubble Space Telescope Wide Field and Planetary Camera 2. Twelve V (F555W) and four I (F814W) epochs of observation were obtained. The two photometry packages, ALLFRAME and DoPHOT, were separately used to obtain profile-fitting photometry of all the stars in the HST field. The search for Cepheid variable stars resulted in a sample of 52 variables, with periods between 14 and 60 days, in common with both datasets. ALLFRAME photometry and light curves of the Cepheids are presented. A subset of 34 Cepheids were selected on the basis of period, light curve shape, similar ALLFRAME and DoPHOT periods, color, and relative crowding, to fit the Cepheid period-luminosity relations in V and I for both ALLFRAME and DoPHOT. The measured distance modulus to NGC 1365 from the ALLFRAME photometry is 31.31 +/- 0.20 (random) +/- 0.18 (systematic) mag, corresponding to a distance of 18.3 +/- 1.7 (random) +/- 1.6 (systematic) Mpc. The reddening is measured to be E(V-I) = 0.16 +/- 0.08 mag. These values are in excellent agreement with those obtained using the DoPHOT photometry, namely a distance modulus of 31.26 +/- 0.10 mag, and a reddening of 0.15 +/- 0.10 mag (internal errors only).