The Milky Way as a High Redshift Galaxy: The Importance of Thick Disk Formation in Galaxies

We compare the star-formation history and dynamics of the Milky Way (MW) with the properties of distant disk galaxies. During the first ~4 Gyr of its evolution, the MW formed stars with a high star-formation intensity (SFI), Sigma_SFR~0.6 Msun/yr/kpc2 and as a result, generated outflows and high turbulence in its interstellar medium. This intense phase of star formation corresponds to the formation of the thick disk. The formation of the thick disk is a crucial phase which enables the MW to have formed approximately half of its total stellar mass by z~1 which is similar to "MW progenitor galaxies" selected by abundance matching. This agreement suggests that the formation of the thick disk may be a generic evolutionary phase in disk galaxies. Using a simple energy injection-kinetic energy relationship between the 1-D velocity dispersion and SFI, we can reproduce the average perpendicular dispersion in stellar velocities of the MW with age. This relationship, its inferred evolution, and required efficiency are consistent with observations of galaxies from z~0-3. The high turbulence generated by intense star formation naturally resulted in a thick disk, a chemically well-mixed ISM, and is the mechanism that links the evolution of MW to the observed characteristics of distant disk galaxies.Comment: 5 pages, 4 figures; accepted to ApJ Letter

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

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&

Long-term radial-velocity variations of the Sun as a star: The HARPS view

Stellar radial velocities play a fundamental role in the discovery of extrasolar planets and the measurement of their physical parameters as well as in the study of stellar physical properties. We investigate the impact of the solar activity on the radial velocity of the Sun using the HARPS spectrograph to obtain measurements that can be directly compared with those acquired in the extrasolar planet search programs. We use the Moon, the Galilean satellites, and several asteroids as reflectors to measure the radial velocity of the Sun as a star and correlate it with disc-integrated chromospheric and magnetic indexes of solar activity that are similar to stellar activity indexes. We discuss in detail the systematic effects that affect our measurements and the methods to account for them. We find that the radial velocity of the Sun as a star is positively correlated with the level of its chromospheric activity at about 95 percent significance level. The amplitude of the long-term variation measured in the 2006-2014 period is 4.98 \pm 1.44 m/s, in good agreement with model predictions. The standard deviation of the residuals obtained by subtracting a linear best fit is 2.82 m/s and is due to the rotation of the reflecting bodies and the intrinsic variability of the Sun on timescales shorter than the activity cycle. A correlation with a lower significance is detected between the radial velocity and the mean absolute value of the line-of-sight photospheric magnetic field flux density. Our results confirm similar correlations found in other late-type main-sequence stars and provide support to the predictions of radial velocity variations induced by stellar activity based on current models.Comment: 11 pages, 7 figures, 2 tables, 1 Appendix; accepted by Astronomy and Astrophysic

Determining the Mass of Kepler-78b With Nonparametric Gaussian Process Estimation

Kepler-78b is a transiting planet that is 1.2 times the radius of Earth and orbits a young, active K dwarf every 8 hours. The mass of Kepler-78b has been independently reported by two teams based on radial velocity measurements using the HIRES and HARPS-N spectrographs. Due to the active nature of the host star, a stellar activity model is required to distinguish and isolate the planetary signal in radial velocity data. Whereas previous studies tested parametric stellar activity models, we modeled this system using nonparametric Gaussian process (GP) regression. We produced a GP regression of relevant Kepler photometry. We then use the posterior parameter distribution for our photometric fit as a prior for our simultaneous GP + Keplerian orbit models of the radial velocity datasets. We tested three simple kernel functions for our GP regressions. Based on a Bayesian likelihood analysis, we selected a quasi-periodic kernel model with GP hyperparameters coupled between the two RV datasets, giving a Doppler amplitude of 1.86 $\pm$ 0.25 m s$^{-1}$ and supporting our belief that the correlated noise we are modeling is astrophysical. The corresponding mass of 1.87 $^{+0.27}_{-0.26}$ M$_{\oplus}$ is consistent with that measured in previous studies, and more robust due to our nonparametric signal estimation. Based on our mass and the radius measurement from transit photometry, Kepler-78b has a bulk density of 6.0$^{+1.9}_{-1.4}$ g cm$^{-3}$. We estimate that Kepler-78b is 32$\pm$26% iron using a two-component rock-iron model. This is consistent with an Earth-like composition, with uncertainty spanning Moon-like to Mercury-like compositions.Comment: 10 pages, 5 figures, accepted to ApJ 6/16/201

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

Planets and Stellar Activity: Hide and Seek in the CoRoT-7 system

Since the discovery of the transiting super-Earth CoRoT-7b, several investigations have yielded different results for the number and masses of planets present in the system, mainly owing to the star's high level of activity. We re-observed CoRoT-7 in January 2012 with both HARPS and CoRoT, so that we now have the benefit of simultaneous radial-velocity and photometric data. This allows us to use the off-transit variations in the star's light curve to estimate the radial-velocity variations induced by the suppression of convective blueshift and the flux blocked by starspots. To account for activity-related effects in the radial-velocities which do not have a photometric signature, we also include an additional activity term in the radial-velocity model, which we treat as a Gaussian process with the same covariance properties (and hence the same frequency structure) as the light curve. Our model was incorporated into a Monte Carlo Markov Chain in order to make a precise determination of the orbits of CoRoT-7b and CoRoT-7c. We measure the masses of planets b and c to be 4.73 +/- 0.95 Mearth and 13.56 +/- 1.08 Mearth, respectively. The density of CoRoT-7b is (6.61 +/- 1.72)(Rp/1.58 Rearth)^(-3) g.cm^(-3), which is compatible with a rocky composition. We search for evidence of an additional planet d, identified by previous authors with a period close to 9 days. We are not able to confirm the existence of a planet with this orbital period, which is close to the second harmonic of the stellar rotation at around 7.9 days. Using Bayesian model selection we find that a model with two planets plus activity-induced variations is most favoured.Comment: Accepted 2014 July 2. Received 2014 June 30; in original form 2013 May 30 (17 pages, 9 figures

Bars & boxy/peanut bulges in thin & thick discs: I. Morphology and line-of-sight velocities of a fiducial model

We explore trends in the morphology and line-of-sight (los) velocity of stellar populations in the inner regions of disc galaxies, using N-body simulations with both a thin (kinematically cold) and a thick (kinematically hot) disc which form a bar and boxy/peanut (b/p) bulge. The bar in the thin disc component is $\sim$50\% stronger than the thick disc bar and is more elongated, with an axis ratio almost half that of the thick disc bar. The thin disc b/p bulge has a pronounced X-shape, while the thick disc b/p is weaker with a rather boxy shape. This leads to the signature of the b/p bulge in the thick disc to be weaker and further away from the plane than in the thin disc. Regarding the kinematics, we find that the los velocity of thick disc stars in the outer parts of the b/p bulge can be \emph{larger} than that of thin disc stars, by up to 40\% and 20\% for side-on and Milky Way-like orientations of the bar respectively. This is due to the different orbits followed by thin and thick disc stars in the bar-b/p region, which are affected by the fact that: i) thin disc stars are trapped more efficiently in the bar - b/p instability and thus lose more angular momentum than their thick disc counterparts and ii) thick disc stars have large radial excursions and therefore stars from large radii with high angular momenta can be found in the bar region. We also find that the difference between the los velocities of the thin and thick disc in the b/p bulge ($\Delta v_{los}$) correlates with the initial difference between the radial velocity dispersions of the two discs ($\Delta \sigma$) . We therefore conclude that stars in the bar - b/p bulge will have considerably different morphologies and kinematics depending on the kinematic properties of the disc population they originate from.Comment: Accepted for publication in A&A. 15 pages (2 page appendix). 16 figure