Quantifying stellar radial migration in an N-body simulation: blurring, churning, and the outer regions of galaxy discs

Radial stellar migration in galactic discs has received much attention in studies of galactic dynamics and chemical evolution, but remains a dynamical phenomenon that needs to be fully quantified. In this work, using a Tree-SPH simulation of an Sb-type disc galaxy, we quantify the effects of blurring (epicyclic excursions) and churning (change of guiding radius). We quantify migration (either blurring or churning) both in terms of flux (the number of migrators passing at a given radius), and by estimating the population of migrators at a given radius at the end of the simulation compared to non-migrators, but also by giving the distance over which the migration is effective at all radii. We confirm that the corotation of the bar is the main source of migrators by churning in a bar-dominated galaxy, its intensity being directly linked to the episode of a strong bar, in the first 1-3 Gyr of the simulation. We show that within the outer Lindblad resonance (OLR), migration is strongly dominated by churning, while blurring gains progressively more importance towards the outer disc and at later times. Most importantly, we show that the OLR limits the exchange of angular momentum, separating the disc in two distinct parts with minimal or null exchange, except in the transition zone, which is delimited by the position of the OLR at the epoch of the formation of the bar, and at the final epoch. We discuss the consequences of these findings for our understanding of the structure of the Milky Way disc. Because the Sun is situated slightly outside the OLR, we suggest that the solar vicinity may have experienced very limited churning from the inner disc.Comment: Accepted for publication in Astronomy and Astrophysics (acceptance date: 27/04/15), 24 pages, 24 figure

Hiding its age: the case for a younger bulge

The determination of the age of the bulge has led to two contradictory results. On the one side, the color-magnitude diagrams in different bulge fields seem to indicate a uniformly old ($>$10 Gyr) population. On the other side, individual ages derived from dwarfs observed through microlensing events seem to indicate a large spread, from $\sim$ 2 to $\sim$ 13 Gyr. Because the bulge is now recognised as being mainly a boxy peanut-shaped bar, it is suggested that disk stars are one of its main constituents, and therefore also stars with ages significantly younger than 10 Gyr. Other arguments as well point to the fact that the bulge cannot be exclusively old, and in particular cannot be a burst population, as it is usually expected if the bulge was the fossil remnant of a merger phase in the early Galaxy. In the present study, we show that given the range of metallicities observed in the bulge, a uniformly old population would be reflected into a significant spread in color at the turn-off which is not observed. Inversely, we demonstrate that the correlation between age and metallicity expected to hold for the inner disk would conspire to form a color-magnitude diagram with a remarkably small spread in color, thus mimicking the color-magnitude diagram of a uniformly old population. If stars younger than 10 Gyr are part of the bulge, as must be the case if the bulge has been mainly formed through dynamical instabilities in the disk, then a very small spread at the turn-off is expected, as seen in the observations.Comment: 11 pages, 11 figures. Accepted for publication in A&

The age structure of stellar populations in the solar vicinity. Clues of a two-phase formation history of the Milky Way disk

We analyze high quality abundances data of solar neighborhood stars and show that there are two distinct regimes of [alpha/Fe] versus age which we identify as the epochs of the thick and thin disk formation. A tight correlation between metallicity and [alpha/Fe] versus age is clearly identifiable on thick disk stars, implying that this population formed from a well mixed ISM, over a time scale of 4-5 Gyr. Thick disk stars vertical velocity dispersion correlate with age, with the youngest objects having as small scale heights as those of thin disk stars. A natural consequence of these two results is that a vertical metallicity gradient is expected in this population. We suggest that the thick disk set the initial conditions for the formation of the inner thin disk. This provides also an explanation of the apparent coincidence between the step in metallicity at 7-10 kpc in the thin disk and the confinment of the thick disk at about R<10 kpc. We suggest that the outer thin disk developped outside the influence of the thick disk, but also that the high alpha-enrichment of the outer regions may originate from a primordial pollution by the gas expelled from the thick disk. Local metal-poor thin disk stars, whose properties are best explained by an origin in the outer disk, are shown to be as old as the youngest thick disk (9-10 Gyr), implying that the outer thin disk started to form while the thick disk formation was still on-going in the inner Galaxy. We point out that, given the tight age-abundance relations in the thick disk, an inside-out process would give rise to a radial gradient in abundances in this population which is not observed. Finally, we argue that the data discussed here leave little room for radial migration, either to have contaminated the solar vicinity, or, to have redistributed stars in significant proportion across the solar annulus.Comment: Accepted in A&A, Revised version with new figures and extended discussio

When the Milky Way turned off the lights: APOGEE provides evidence of star formation quenching in our Galaxy

Quenching, the cessation of star formation, is one of the most significant events in the life cycle of galaxies. We show here the first evidence that the Milky Way experienced a generalised quenching of its star formation at the end of its thick disk formation $\sim$9 Gyr ago. Elemental abundances of stars studied as part of the APOGEE survey reveal indeed that in less than $\sim$2 Gyr the star formation rate in our Galaxy dropped by an order-of-magnitude. Because of the tight correlation between age and alpha abundance, this event reflects in the dearth of stars along the inner disk sequence in the [Fe/H]-[$\alpha$/Fe] plane. Before this phase, which lasted about 1.5 Gyr, the Milky Way was actively forming stars. Afterwards, the star formation resumed at a much lower level to form the thin disk. These events are very well matched by the latest observation of MW-type progenitors at high redshifts. In late type galaxies, quenching is believed to be related to a long and secular exhaustion of gas. In our Galaxy, it occurred on a much shorter time scale, while the chemical continuity before and after the quenching indicates that it was not due to the exhaustion of the gas. While quenching is generally associated with spheroids, our results show that it also occurs in galaxies like the Milky Way, possibly when they are undergoing a morphological transition from thick to thin disks. Given the demographics of late type galaxies in the local universe, in which classical bulges are rare, we suggest further that this may hold true generally in galaxies with mass lower than or approximately $M^*$, where quenching could be directly a consequence of thick disk formation. We emphasize that the quenching phase in the Milky Way could be contemporaneous with, and related to, the formation of the bar. We sketch a scenario on how a strong bar may inhibit star formation.Comment: 17 pages, 8 figures. Published versio

On the kinematic detection of accreted streams in the Gaia era: a cautionary tale

The $\Lambda$CDM cosmological scenario predicts that our Galaxy should contain hundreds of stellar streams at the solar vicinity, fossil relics of the merging history of the Milky Way and more generally of the hierarchical growth of galaxies. Because of the mixing time scales in the inner Galaxy, it has been claimed that these streams should be difficult to detect in configuration space but can still be identifiable in kinematic-related spaces like the energy/angular momenta spaces, E-Lz and Lperp-Lz, or spaces of orbital/velocity parameters. By means of high-resolution, dissipationless N-body simulations, containing between 25$\times10^6$ and 35$\times10^6$ particles, we model the accretion of a series of up to four 1:10 mass ratio satellites then up to eight 1:100 satellites and we search systematically for the signature of these accretions in these spaces. In all spaces considered (1) each satellite gives origin to several independent overdensities; (2) overdensities of multiple satellites overlap; (3) satellites of different masses can produce similar substructures; (4) the overlap between the in-situ and the accreted population is considerable everywhere; (5) in-situ stars also form substructures in response to the satellite(s) accretion. These points are valid even if the search is restricted to kinematically-selected halo stars only. As we are now entering the 'Gaia era', our results warn that an extreme caution must be employed before interpreting overdensities in any of those spaces as evidence of relics of accreted satellites. Reconstructing the accretion history of our Galaxy will require a substantial amount of accurate spectroscopic data, that, complemented by the kinematic information, will possibly allow us to (chemically) identify accreted streams and measure their orbital properties. (abridged)Comment: Accepted on A&A. A high-resolution version of the paper is available at http://aramis.obspm.fr/~paola/ELZ/Elz.pd

Composting paper and grass clippings with anaerobically treated palm oil mill effluent

Purpose The purpose of this study is to investigate the composting performance of anaerobically treated palm oil mill effluent (AnPOME) mixed with paper and grass clippings. Methods Composting was conducted using a laboratory scale system for 40 days. Several parameters were determined: temperature, mass reduction, pH, electrical conductivity, colour, zeta potential, phytotoxicity and final compost nutrients. Results The moisture content and compost mass were reduced by 24 and 18 %, respectively. Both final compost pH value and electrical conductivity were found to increase in value. Colour (measured as PtCo) was not suitable as a maturity indicator. The negative zeta potential values decreased from −12.25 to −21.80 mV. The phytotoxicity of the compost mixture was found to decrease in value during the process and the final nutrient value of the compost indicates its suitability as a soil conditioner. Conclusions From this study, we conclude that the addition of paper and grass clippings can be a potential substrate to be composted with anaerobically treated palm oil mill effluent (AnPOME). The final compost produced is suitable for soil conditioner

The Galactic disk mass-budget : II. Brown dwarf mass-function and density

In this paper, we extend the calculations conducted previously in the stellar regime to determine the brown dwarf IMF in the Galactic disk. We perform Monte Carlo calculations taking into account the brown dwarf formation rate, spatial distribution and binary fraction. Comparison with existing surveys seems to exclude a power-law MF as steep as the one determined in the stellar regime below 1 \msol and tends to favor a more flatish behaviour. Comparison with methane-dwarf detections tends to favor an eventually decreasing form like the lognormal or the more general exponential distributions determined in the previous paper. We calculate predicting brown dwarf counts in near-infrared color diagrams and brown dwarf discovery functions. These calculations yield the presently most accurate determination of the brown dwarf census in the Galactic disk. The brown dwarf number density is comparable to the stellar one, $n_{BD}\simeq n_\star\simeq 0.1$ pc$^{-3}$. The corresponding brown dwarf mass density, however, represents only about 10% of the stellar contribution, i.e. \rho_{BD}\simle 5.0\times 10^{-3} \mvol. Adding up the local stellar density determined previously yields the density of star-like objects, stars and brown dwarfs, in the solar neighborhood \rho_\odot \approx 5.0\times 10^{-2} \mvol.Comment: 39 pages, Latex file, uses aasms4.sty, to be published in ApJ, corrected version with correct figure

A new look at the kinematics of the bulge from an N-body model

(Abridged) By using an N-body simulation of a bulge that was formed via a bar instability mechanism, we analyse the imprints of the initial (i.e. before bar formation) location of stars on the bulge kinematics, in particular on the heliocentric radial velocity distribution of bulge stars. Four different latitudes were considered: $b=-4^\circ$, $-6^\circ$, $-8^\circ$, and $-10^\circ$, along the bulge minor axis as well as outside it, at $l=\pm5^\circ$ and $l=\pm10^\circ$. The bulge X-shaped structure comprises stars that formed in the disk at different locations. Stars formed in the outer disk, beyond the end of the bar, which are part of the boxy peanut-bulge structure may show peaks in the velocity distributions at positive and negative heliocentric radial velocities with high absolute values that can be larger than 100 $\rm km$ $\rm s^{-1}$, depending on the observed direction. In some cases the structure of the velocity field is more complex and several peaks are observed. Stars formed in the inner disk, the most numerous, contribute predominantly to the X-shaped structure and present different kinematic characteristics. Our results may enable us to interpret the cold high-velocity peak observed in the APOGEE commissioning data, as well as the excess of high-velocity stars in the near and far arms of the X-shaped structure at $l$=$0^\circ$ and $b$=$-6^\circ$. When compared with real data, the kinematic picture becomes more complex due to the possible presence in the observed samples of classical bulge and/or thick disk stars. Overall, our results point to the existence of complex patterns and structures in the bulge velocity fields, which are generated by the bar. This suggests that caution should be used when interpreting the bulge kinematics: the presence of substructures, peaks and clumps in the velocity fields is not necessarily a sign of past accretion events.Comment: 21 pages, 18 figures. Accepted for publication in A&

The contrasting climate response to tropical and extratropical energy perturbations

This is the final version of the article. Available from Springer Verlag via the DOI in this record.The link between cross-equatorial energy transport, the double-intertropical convergence zone (DI) problem and biases in tropical and extratropical albedo and energy budgets in climate models have been investigated in multiple studies, though DI biases persist in many models. Here, a coupled climate model, HadGEM2-ES, is used to investigate the response to idealised energy perturbations in the tropics and extratropics, in both the northern and southern hemispheres, through the imposition of stratospheric aerosols that reflect incoming radiation. The impact on the tropical climate of high and low latitude forcing strongly contrasts, with large changes in tropical precipitation and modulation of the DI bias when the tropics are cooled as precipitation moves away from the cooled hemisphere. These responses are muted when the extratropics are cooled, as the meridional energy transport anomalies that are excited by these energy budget anomalies are partitioned between the atmosphere and ocean. The results here highlight the persistence of the DI bias in HadGEM2-ES, indicating why little progress has been made in rectifying these problems through many generations of climate models. A highly linear relationship between cross-equatorial atmospheric energy transport, tropical precipitation asymmetry and tropical sea surface temperature biases is also demonstrated, giving some suggestion as to where improvements in these large scale, persistent biases may be achieved.MKH, MC and JMH were supported by the Natural Environment Research Council/Department for International Development via the Future Climates for Africa (FCFA) funded project ’Improving Model Processes for African Climate’ (IMPALA, NE/M017265/1). JMH and AJ were supported by the Joint UK BEIS/Defra Met Office Hadley Centre Climate Programme (GA01101)

Climatic impacts of stratospheric geoengineering with sulfate, black carbon and titania injection

In this paper, we examine the potential climatic effects of geoengineering by sulfate, black carbon and titania injection against a baseline RCP8.5 scenario. We use the HadGEM2-CCS model to simulate scenarios in which the top-of-the-atmosphere radiative imbalance due to rising greenhouse gas concentrations is offset by sufficient aerosol injection throughout the 2020–2100 period. We find that the global-mean temperature is effectively maintained at historical levels for the entirety of the period for all 3 aerosol-injection scenarios, though there are a wide range of side-effects which are discussed in detail. The most prominent conclusion is that although the BC injection rate necessary to produce an equivalent global mean temperature-response is much lower, the severity of stratospheric temperature changes (> +70 °C) and precipitation impacts effectively exclude BC from being a viable option for geoengineering. Additionally, while it has been suggested that titania would be an effective particle because of its high scattering efficiency, it also efficiently absorbs solar ultraviolet radiation producing a significant stratospheric warming (> +20 °C). As injection rates for titania are close to those for sulfate, there appears little benefit of using titania when compared to injection of sulfur dioxide, which has the added benefit of being well modelled through extensive research that has been carried out on naturally occurring explosive volcanic eruptions.The authors would like to thank Valentina Aquila for supplying AVHRR and SAGE data, and to Peter Cox, Angus Ferraro, David Keith and Alan Robock for helpful discussions. A. C. Jones was supported by a Met Office/NERC CASE (ref. 580 009 183) PhD studentship; J. M. Haywood and A. Jones were supported by the Joint UK DECC/Defra Met Office Hadley Centre Climate Programme (GA01101)