Milky Way demographics with the VVV survey I. The 84-million star colour-magnitude diagram of the Galactic bulge

Context. The Milky Way (MW) bulge is a fundamental Galactic component for understanding the formation and evolution of galaxies, in particular our own. The ESO Public Survey VISTA Variables in the Via Lactea is a deep near-IR survey mapping the Galactic bulge and southern plane. Particularly for the bulge area, VVV is covering similar to 315 deg(2). Data taken during 2010 and 2011 covered the entire bulge area in the JHKs bands. Aims. We used VVV data for the whole bulge area as a single and homogeneous data set to build for the first time a single colour-magnitude diagram (CMD) for the entire Galactic bulge. Methods. Photometric data in the JHK(s) bands were combined to produce a single and huge data set containing 173 150 467 sources in the three bands, for the similar to 315 deg(2) covered by VVV in the bulge. Selecting only the data points flagged as stellar, the total number of sources is 84 095 284. Results. We built the largest colour-magnitude diagrams published up to date, containing 173.1+ million sources for all data points, and more than 84.0 million sources accounting for the stellar sources only. The CMD has a complex shape, mostly owing to the complexity of the stellar population and the effects of extinction and reddening towards the Galactic centre. The red clump (RC) giants are seen double in magnitude at b similar to -8 degrees-10 degrees, while in the inner part (b similar to -3 degrees) they appear to be spreading in colour, or even splitting into a secondary peak. Stellar population models show the predominance of main-sequence and giant stars. The analysis of the outermost bulge area reveals a well-defined sequence of late K and M dwarfs, seen at (J - K-s) similar to 0.7-0.9 mag and K-s greater than or similar to 14 mag. Conclusions. The interpretation of the CMD yields important information about the MW bulge, showing the fingerprint of its structure and content. We report a well-defined red dwarf sequence in the outermost bulge, which is important for the planetary transit searches of VVV. The double RC in magnitude seen in the outer bulge is the signature of the X-shaped MW bulge, while the spreading of the RC in colour, and even its splitting into a secondary peak, are caused by reddening effects. The region around the Galactic centre is harder to interpret because it is strongly affected by reddening and extinction.

The star formation rate in the inner Milky Way Galaxy

The present star formation rate (SFR) in the inner Galaxy is puzzling for the chemical evolution models (CEM). No static CEM is able to reproduce the peak of the SFR in the 4 kpc ring. The main reason is probably a shortage of gas, which could be due to the dynamical effects produced by the galactic bar, not considered by these models. We developed a CEM that includes radial gas flows in order to mimic the effects of the galactic bar in the first 5 kpc of the galactic disk. In this model, the star formation (SF) is a two-step process: first, the diffuse gas forms molecular clouds. Then, stars form from cloud-cloud collisions or by the interaction between massive stars and the molecular gas. The former is called spontaneous and the latter induced SF. The mass in the different phases of each region changes by the processes associated with the stellar formation and death by: the SF due to spontaneous fragmentation of gas in the halo; formation of gas clouds in the disk from the diffuse gas; induced SF in the disk due to the interaction between massive stars and gas clouds; and finally, the restitution of the diffuse gas associated to these process of cloud and star formation. In the halo, the star formation rate for the diffuse gas follows a Schmidt law with a power n = 1.5. In the disk, the stars form in two steps: first, molecular clouds are formed from the diffuse gas also following a Schmidt law with n=1.5 and a proportionality factor. Including a specific pattern of radial gas flows, the CEM is able to reproduce with success the peak in the SFR at 4 kpc (fig. 1).

The gaseous halo of the Milky Way

In the last decade, sensitive observations have revealed that disc galaxies are surrounded by multiphase gaseous halos produced by the circulation of gas from the discs to the environment and vice-versa. This Thesis is a study of the gaseous halo of the Milky Way carried out via the modelling of the HI emission and the available absorption-line data. We fitted simple kinematical models to the HI LAB Survey and found that the Galaxy has a massive (~3x10^8 Mo) HI halo extending a few kiloparsecs above the plane. This layer rotates more slowly than the disc and shows a global inflow motion, a kinematics similar to that observed in the HI halos of nearby galaxies. We built a dynamical model of the galactic fountain to reproduce the properties of this layer. In this model, fountain clouds are ejected from the disc by SN feedback and - as suggested by hydrodynamical simulations - triggers the cooling of coronal gas, which is entrained by the cloud wakes and accretes onto the disc when the clouds fall back. For a proper choice of the parameters, the model reproduces well the HI data and predicts an accretion of coronal gas onto the disc at a rate of 2 Mo/yr. We extended this model to the warm-hot component of the halo, showing that most of the ion absorption features observed towards background sources are consistent with being produced in the turbulent wakes that lag behind the fountain clouds. Specifically, the column densities, positions, and velocities of the absorbers are well reproduced by our model. Finally, we studied the gas content of galaxies extracted from a cosmological N-body+SPH simulation, and found that an HI halo with the forementioned properties is not observed, probably due ti the relatively low resolution of the simulations.

Black holes in the Milky Way Galaxy

Extremely strong observational evidence has recently been found for the presence of black holes orbiting a few relatively normal stars in our Milky Way Galaxy and also at the centers of some galaxies. The former generally have masses of 4–16 times the mass of the sun, whereas the latter are “supermassive black holes” with millions to billions of solar masses. The evidence for a supermassive black hole in the center of our galaxy is especially strong.

Circumstellar emission in Be/X-ray binaries of the Magellanic Clouds and the Milky Way

Aims. We study the optical and near-infrared colour excesses produced by circumstellar emission in a sample of Be/X-ray binaries. Our main goals are exploring whether previously published relations, valid for isolated Be stars, are applicable to Be/X-ray binaries and computing the distance to these systems after correcting for the effects of the circumstellar contamination. Methods. Simultaneous UBVRI photometry and spectra in the 3500−7000 Å spectral range were obtained for 11 optical counterparts to Be/X-ray binaries in the LMC, 5 in the SMC and 12 in the Milky Way. As a measure of the amount of circumstellar emission we used the Hα equivalent width corrected for photospheric absorption. Results. We find a linear relationship between the strength of the Hα emission line and the component of E(B − V) originating from the circumstellar disk. This relationship is valid for stars with emission lines weaker than EW ≈ −15   Å. Beyond this point, the circumstellar contribution to E(B − V) saturates at a value ≈0.17   mag. A similar relationship is found for the (V − I) near infrared colour excess, albeit with a steeper slope and saturation level. The circumstellar excess in (B − V) is found to be about five times higher for Be/X-ray binaries than for isolated Be stars with the same equivalent width EW(Hα), implying significant differences in the physical properties of their circumstellar envelopes. The distance to Be/X-ray binaries (with non-shell Be star companions) can only be correctly estimated by taking into account the excess emission in the V band produced by free-free and free-bound transitions in the circumstellar envelope. We provide a simple method to determine the distances that includes this effect.

VdBH 222: a starburst cluster in the inner Milky Way

Context. It has been suggested that the compact open cluster VdBH 222 is a young massive distant object. Aims. We set out to characterise VdBH 222 using a comprehensive set of multi-wavelength observations. Methods. We obtained multi-band optical (UBVR) and near-infrared (JHKS) photometry of the cluster field, as well as multi-object and long-slit optical spectroscopy for a large sample of stars in the field. We applied classical photometric analysis, as well as more sophisticated methods using the CHORIZOS code, to determine the reddening to the cluster. We then plotted dereddened HR diagrams and determined cluster parameters via isochrone fitting. Results. We have identified a large population of luminous supergiants confirmed as cluster members via radial velocity measurements. We find nine red supergiants (plus one other candidate) and two yellow supergiants. We also identify a large population of OB stars. Ten of them are bright enough to be blue supergiants. The cluster lies behind ≈7.5 mag of extinction for the preferred value of RV = 2.9. Isochrone fitting allows for a narrow range of ages between 12 and 16 Ma. The cluster radial velocity is compatible with distances of ~6 and ~10 kpc. The shorter distance is inconsistent with the age range and Galactic structure. The longer distance implies an age ≈ 12 Ma and a location not far from the position where some Galactic models place the far end of the Galactic bar. Conclusions. VdBH 222 is a young massive cluster with a likely mass >20 000 M⊙. Its population of massive evolved stars is comparable to that of large associations, such as Per OB1. Its location in the inner Galaxy, presumably close to the end of the Galactic bar, adds to the increasing evidence for vigorous star formation in the inner regions of the Milky Way.

The Gaia-ESO Survey: radial metallicity gradients and age-metallicity relation of stars in the Milky Way disk

We study the relationship between age, metallicity, and α-enhancement of FGK stars in the Galactic disk. The results are based upon the analysis of high-resolution UVES spectra from the Gaia-ESO large stellar survey. We explore the limitations of the observed dataset, i.e. the accuracy of stellar parameters and the selection effects that are caused by the photometric target preselection. We find that the colour and magnitude cuts in the survey suppress old metal-rich stars and young metal-poor stars. This suppression may be as high as 97% in some regions of the age-metallicity relationship. The dataset consists of 144 stars with a wide range of ages from 0.5 Gyr to 13.5 Gyr, Galactocentric distances from 6 kpcto 9.5 kpc, and vertical distances from the plane 0 < |Z| < 1.5 kpc. On this basis, we find that i) the observed age-metallicity relation is nearly flat in the range of ages between 0 Gyr and 8 Gyr; ii) at ages older than 9 Gyr, we see a decrease in [Fe/H] and a clear absence of metal-rich stars; this cannot be explained by the survey selection functions; iii) there is a significant scatter of [Fe/H] at any age; and iv) [Mg/Fe] increases with age, but the dispersion of [Mg/Fe] at ages >9 Gyr is not as small as advocated by some other studies. In agreement with earlier work, we find that radial abundance gradients change as a function of vertical distance from the plane. The [Mg/Fe] gradient steepens and becomes negative. In addition, we show that the inner disk is not only more α-rich compared to the outer disk, but also older, as traced independently by the ages and Mg abundances of stars.

An emission-line survey of the Milky Way /

Includes bibliographical references

Gleanings in Buddha-fields : and, the romance of the Milky Way /

Mode of access: Internet.

A photographic atlas of selected regions of the Milky Way,

Mode of access: Internet.

The wonders of the heavens, being a popular view of astronomy, including a full illustration of the mechanism of the heavens; embracing the sun, moon, and stars, with descriptions of the planets, comets, fixed stars, double stars, the constellations, the galaxy, or milky-way, the zodiacal light, aurora borealis, or northern lights, meteors, clouds, falling stars, aërolites, & c. Illustrated by numerous maps and engravings.

In double columns.

The large- and small-scale Ca IIK structure of the Milky Way from observations of Galactic and Magellanic sightlines

Aims. The large and small-scale (pc) structure of the Galactic interstellar medium can be investigated by utilising spectra of early-type stellar probes of known distances in the same region of the sky. This paper determines the variation in line strength of Ca ii at 3933.661 Å as a function of probe separation for a large sample of stars, including a number of sightlines in the Magellanic Clouds. 

Methods. FLAMES-GIRAFFE data taken with the Very Large Telescope towards early-type stars in 3 Galactic and 4 Magellanic open clusters in Ca ii are used to obtain the velocity, equivalent width, column density, and line width of interstellar Galactic calcium for a total of 657 stars, of which 443 are Magellanic Cloud sightlines. In each cluster there are between 43 and 111 stars observed. Additionally, FEROS and UVES Ca ii K and Na i D spectra of 21 Galactic and 154 Magellanic early-type stars are presented and combined with data from the literature to study the calcium column density - parallax relationship. 

Results. For the four Magellanic clusters studied with FLAMES, the strength of the Galactic interstellar Ca ii K equivalent width on transverse scales from ∼0.05-9 pc is found to vary by factors of ∼1.8-3.0, corresponding to column density variations of ∼0.3-0.5 dex in the optically-thin approximation. Using FLAMES, FEROS, and UVES archive spectra, the minimum and maximum reduced equivalent widths for Milky Way gas are found to lie in the range ∼35-125 mÅ and ∼30-160 mÅ for Ca ii K and Na i D, respectively. The range is consistent with a previously published simple model of the interstellar medium consisting of spherical cloudlets of filling factor ∼0.3, although other geometries are not ruled out. Finally, the derived functional form for parallax (π) and Ca ii column density (N_CaII) is found to be π(mas) = 1 / (2.39 × 10^-13 × N_CaII (cm^-2) + 0.11). Our derived parallax is ∼25 per cent lower than predicted by Megier et al. (2009, A&A, 507, 833) at a distance of ∼100 pc and ∼15 percent lower at a distance of ∼200 pc, reflecting inhomogeneity in the Ca ii distribution in the different sightlines studied.

Analytic modelling of the galactic fountain in the Milky Way potential

Le galassie a spirale, come la Via Lattea, sono caratterizzate dalla presenza di gas freddo e formazione stellare e vengono perciò chiamate star-forming. Per creare nuove stelle è necessaria una sufficiente riserva di gas, la cui disponibilità governa l’evoluzione della galassia stessa. Finora, non è stato individuato con certezza un meccanismo che possa alimentare la formazione di nuove stelle nelle galassie star-forming. Una delle possibili sorgenti di tale gas è l’alone galattico caldo (corona galattica) il cui raffreddamento e successivo accrescimento possono essere stimolati dal processo di fontana galattica. L’esplosione di supernovae porta nubi di gas freddo in orbita al di sopra del disco stellare; queste nubi raggiungono altezze dell’ordine del kiloparsec, interagendo con la corona di gas caldo. Il moto delle nubi all’interno di un mezzo meno denso comporta l’instaurarsi dell’instabilità di Kelvin-Helmholtz, che ’strappa’ gas dalle nubi e causa la condensazione di materia coronale. Quest’ultima viene quindi accresciuta e, ricadendo sul disco, trasferisce nuovo materiale alla galassia e ne alimenta la formazione stellare. Lo scopo di questa tesi è derivare un modello analitico di fontana galattica che consenta di ottenere una formulazione analitica per il tempo orbitale, cioè il tempo richiesto alle nubi per ricadere sul disco galattico. Infatti, più u tempo le nubi impiegano per attraversare il materiale coronale caldo e ricadere sul disco, più materiale viene accresciuto durante l’orbita. Conoscendo i tempi orbitali sarebbe possibile calcolare il tasso di accrescimento legato al fenomeno di fontana e studiarne l’andamento con il raggio del disco. Questo modello potrebbe rivelarsi utile per lo studio dell’impatto della fontana nell’evoluzione globale del disco galattico e della chimica dell’intera galassia.

Mock HI observations of Milky Way-like simulated galaxies

L'obiettivo di questo lavoro di tesi consiste nel descrivere sia il processo necessario per la creazione di osservazioni sintetiche di galassie simulate simili alla Via Lattea nella riga di emissione a 21 cm dell'idrogeno neutro (HI), sia il lavoro di analisi fondamentale che serve a confrontare in modo efficace l'output generato con delle osservazioni di galassie reali. Come prima cosa è descritta la teoria quantistica che sta alla base dell'emissione a 21 cm di HI, illustrando l'importanza di tale riga di emissione nell'ambito dell'astronomia e come si possano ottenere informazioni fondamentali sulle sorgenti di questa radiazione a partire dai dati osservativi. Il lavoro poi si focalizza sull'utilizzo del software MARTINI per la creazione di osservazioni sintetiche della linea a 21 cm per una galassia simulata con proprietà simili alla Via Lattea generata utilizzando il modello numerico SMUGGLE. Infine, si passa ad una breve descrizione dell'analisi dei dati sintetici creati, e al loro confronto con dei dati provenienti da osservazioni reali di galassie con proprietà simili, per ottenere una valutazione qualitativa della bontà del modello SMUGGLE impiegato nella simulazione numerica.

Large change in the predicted number of small halos due to a small amplitude oscillating inflaton potential

A smooth inflaton potential is generally assumed when calculating the primordial power spectrum, implicitly assuming that a very small oscillation in the inflaton potential creates a negligible change in the predicted halo mass function. We show that this is not true. We find that a small oscillating perturbation in the inflaton potential in the slow-roll regime can alter significantly the predicted number of small halos. A class of models derived from supergravity theories gives rise to inflaton potentials with a large number of steps and many trans-Planckian effects may generate oscillations in the primordial power spectrum. The potentials we study are the simple quadratic (chaotic inflation) potential with superimposed small oscillations for small field values. Without leaving the slow-roll regime, we find that for a wide choice of parameters, the predicted number of halos change appreciably. For the oscillations beginning in the 10(7)-10(8) M(circle dot) range, for example, we find that only a 5% change in the amplitude of the chaotic potential causes a 50% suppression of the number of halos for masses between 10(7)-10(8) M(circle dot) and an increase in the number of halos for masses <10(6) M(circle dot) by factors similar to 15-50. We suggest that this might be a solution to the problem of the lack of observed dwarf galaxies in the range 10(7)-10(8) M(circle dot). This might also be a solution to the reionization problem where a very large number of Population III stars in low mass halos are required.