Improving the cosmic distance ladder. Distance and structure of the Large Magellanic Cloud

The Large Magellanic Cloud (LMC) is widely considered as the first step of the cosmological distance ladder, since it contains many different distance indicators. An accurate determination of the distance to the LMC allows one to calibrate these distance indicators that are then used to measure the distance to far objects. The main goal of this thesis is to study the distance and structure of the LMC, as traced by different distance indicators. For these purposes three types of distance indicators were chosen: Classical Cepheids,``hot'' eclipsing binaries and RR Lyrae stars. These objects belong to different stellar populations tracing, in turn, different sub-structures of the LMC. The RR Lyrae stars (age >10 Gyr) are distributed smoothly and likely trace the halo of the LMC. Classical Cepheids are young objects (age 50-200 Myr), mainly located in the bar and spiral arm of the galaxy, while ``hot'' eclipsing binaries mainly trace the star forming regions of the LMC. Furthermore, we have chosen these distance indicators for our study, since the calibration of their zero-points is based on fundamental geometric methods. The ESA cornerstone mission Gaia, launched on 19 December 2013, will measure trigonometric parallaxes for one billion stars with an accuracy of 20 micro-arcsec at V=15 mag, and 200 micro-arcsec at V=20 mag, thus will allow us to calibrate the zero-points of Classical Cepheids, eclipsing binaries and RR Lyrae stars with an unprecedented precision.

The Evolution of Massive Young Stellar Objects in the Large Magellanic Cloud

This thesis presents an analysis of the largest catalog to date of infrared spectra of massive young stellar objects in the Large Magellanic Cloud. Evidenced by their very different spectral features, the luminous objects span a range of evolutionary states from those most embedded in their natal molecular material to those that have dissipated and ionized their surroundings to form compact HII regions and photodissociation regions. We quantify the contributions of the various spectral features using the statistical method of principal component analysis. Using this analysis, we classify the YSO spectra into several distinct groups based upon their dominant spectral features: silicate absorption (S Group), silicate absorption and fine-structure line emission (SE), polycyclic aromatic hydrocarbon (PAH) emission (P Group), PAH and fine-structure line emission (PE), and only fine-structure line emission (E). Based upon the relative numbers of sources in each category, we are able to estimate the amount of time massive YSOs spend in each evolutionary stage. We find that approximately 50% of the sources have ionic fine-structure lines, indicating that a compact HII region forms about half-way through the YSO lifetime probed in our study. Of the 277 YSOs we collected spectra for, 41 have ice absorption features, indicating they are surrounded by cold ice-bearing dust particles. We have decomposed the shape of the ice features to probe the composition and thermal history of the ice. We find that most the CO2 ice is embedded a polar ice matrix that has been thermally processed by the embedded YSO. The amount of thermal processing may be correlated with the luminosity of the YSO. Using the Australia Telescope Compact Array, we imaged the dense gas around a subsample of our sources in the HII complexes N44, N105, N113, and N159 using HCO+ and HCN as dense gas tracers. We find that the molecular material in star forming environments is highly clumpy, with clumps that range from subparsec to ~2 parsecs in size and with masses between 10^2 to 10^4 solar masses. We find that there are varying levels of star formation in the clumps, with the lower-mass clumps tending to be without massive YSOs. These YSO-less clumps could either represent an earlier stage of clump to the more massive YSO-bearing ones or clumps that will never form a massive star. Clumps with massive YSOs at their centers have masses larger than those with massive YSOs at their edges, and we suggest that the difference is evolutionary: edge YSO clumps are more advanced than those with YSOs at their centers. Clumps with YSOs at their edges may have had a significant fraction of their mass disrupted or destroyed by the forming massive star. We find that the strength of the silicate absorption seen in YSO IR spectra feature is well-correlated with the on-source HCO+ and HCN flux densities, such that the strength of the feature is indicative of the embeddedness of the YSO. We estimate that ~40% of the entire spectral sample has strong silicate absorption features, implying that the YSOs are embedded in circumstellar material for about 40% of the time probed in our study.

A non-LTE analysis of the spectra of two narrow lined main sequence stars in the SMC

An analysis of high-resolution VLT/UVES spectra of two B-type main sequence stars, NGC 346-11 and AV 304, in the Small Magellanic Cloud (SMC), has been undertaken, using the non-LTE tlusty model atmospheres to derive the stellar parameters and chemical compositions of each star. The chemical compositions of the two stars are in reasonable agreement. Moreover, our stellar analysis agrees well with earlier analyses of H II regions. The results derived here should be representative of the current base-line chemical composition of the SMC interstellar medium as derived from B-type stars.

The present-day chemical composition of the SMC from UVES spectra of the sharp-lined, B-type dwarf AV 304

High-resolution spectroscopic VLT/UVES observations are presented for the B-type main-sequence star, AV 304, in the Small Magellanic Cloud (SMC). These spectra have been analysed using LTE model-atmosphere techniques, to derive stellar atmospheric parameters and chemical compositions. As AV 304 is located within the hydrogen burning main-sequence band, its chemical composition should reflect that of the SMC interstellar medium (ISM). A detailed line-by-line differential analysis has been undertaken relative to a Galactic comparison star. A general metal deficiency for the a-process elements O, Si & S of -0.43 +/- 0.05 dex is found for AV 304, with iron having a similar underabundance. Oxygen may be relatively over- abundant by similar to0.1 dex and carbon and aluminium underabundant by similar to0.2 dex. A large nitrogen underabundance (of -1.2 dex relative to hydrogen and -0.7 dex relative to iron) is found. This is interpreted in terms of the CNO bi-cycle having been suppressed in the SMC. Furthermore, the large nitrogen deficiency is in excellent agreement with that found for SMC H II regions. Indeed, this represents a first for stellar astrophysics - confirming the low base-line nitrogen composition of the SMC ISM (viz. 12+log(N/H) similar to 6.66 +/- 0.10 dex).

The VLT-FLAMES survey of massive stars: mass loss and rotation of early-type stars in the SMC

We have studied the optical spectra of a sample of 31 O- and early B-type stars in the Small Magellanic Cloud, 21 of which are associated with the young massive cluster NGC 346. Stellar parameters are determined using an automated fitting method (Mokiem et al. 2005, A&A, 441, 711), which combines the stellar atmosphere code FASTWIND (Puls et al. 2005, A&A, 435, 669) with the genetic algorithm based optimisation routine PIKAIA (Charbonneau 1995, ApJS, 101, 309). Comparison with predictions of stellar evolution that account for stellar rotation does not result in a unique age, though most stars are best represented by an age of 1-3 Myr. The automated method allows for a detailed determination of the projected rotational velocities. The present day v(r) sin i distribution of the 21 dwarf stars in our sample is consistent with an underlying rotational velocity (v(r)) distribution that can be characterised by a mean velocity of about 160-190 km s(-1) and an effective half width of 100-150 km s(-1). The vr distribution must include a small percentage of slowly rotating stars. If predictions of the time evolution of the equatorial velocity for massive stars within the environment of the SMC are correct (Maeder & Meynet 2001, A&A, 373, 555), the young age of the cluster implies that this underlying distribution is representative for the initial rotational velocity distribution. The location in the Hertzsprung-Russell diagram of the stars showing helium enrichment is in qualitative agreement with evolutionary tracks accounting for rotation, but not for those ignoring vr. The mass loss rates of the SMC objects having luminosities of log L-star/L-circle dot greater than or similar to 5.4 are in excellent agreement with predictions by Vink et al. (2001, A&A, 369, 574). However, for lower luminosity stars the winds are too weak to determine. M accurately from the optical spectrum. Three targets were classified as Vz stars, two of which are located close to the theoretical zero-age main sequence. Three lower luminosity targets that were not classified as Vz stars are also found to lie near the ZAMS. We argue that this is related to a temperature effect inhibiting cooler from displaying the spectral features required for the Vz luminosity class.

Multiple major outbursts from a restless luminous blue variable in NGC 3432

We present new photometric and spectroscopic observations of an unusual luminous blue variable (LBV) in NGC 3432, covering three major outbursts in 2008 October, 2009 April and 2009 November. Previously, this star experienced an outburst also in 2000 (known as SN 2000ch). During outbursts the star reached an absolute magnitude between -12.1 and -12.8. Its spectrum showed H, He I and Fe II lines with P-Cygni profiles during and soon after the eruptive phases, while only intermediate-width lines in pure emission (including He II lambda 4686) were visible during quiescence. The fast-evolving light curve soon after the outbursts, the quasi-modulated light curve, the peak magnitude and the overall spectral properties are consistent with multiple episodes of variability of an extremely active LBV. However, the widths of the spectral lines indicate unusually high wind velocities (1500-2800 km s-1), similar to those observed in Wolf-Rayet stars. Although modulated light curves are typical of LBVs during the S-Dor variability phase, the luminous maxima and the high frequency of outbursts are unexpected in S-Dor variables. Such extreme variability may be associated with repeated ejection episodes during a giant eruption of an LBV. Alternatively, it may be indicative of a high level of instability shortly preceding the core-collapse or due to interaction with a massive, binary companion. In this context, the variable in NGC 3432 shares some similarities with the famous stellar system HD 5980 in the Small Magellanic Cloud, which includes an erupting LBV and an early Wolf-Rayet star.

The chemical composition of two supergiants in the dwarf irregular galaxy WLM

The chemical composition of two stars in WLM has been determined from high-quality Ultraviolet-Visual Echelle Spectrograph (UVES) data obtained at the VLT-UT2. The model atmospheres analysis shows that they have the same metallicity, [Fe/H] = - 0.38 +/- 0.20 (+/- 0.29). Reliable magnesium abundances are determined from several lines of two ionization states in both stars resulting in [Mg/Fe] = - 0.24 +/- 0.16 (+/- 0.28). This result suggests that the [alpha(Mg)/Fe] ratio in WLM may be suppressed relative to solar abundances ( also supported by differential abundances relative to similar stars in NGC 6822 and the Small Magellanic Cloud [SMC]). The absolute Mg abundance, [Mg/H] = -0.62, is high relative to what is expected from the nebulae though, where two independent spectroscopic analyses of the H II regions in WLM yield [O/H] = - 0.89. Intriguingly, the oxygen abundance determined from the O I lambda6158 feature in one WLM star is [O/H] = - 0.21 +/- 0.10 (+/- 0.05), corresponding to 5 times higher than the nebular oxygen abundance. This is the first time that a significant difference between stellar and nebular oxygen abundances has been found, and currently, there is no simple explanation for this difference. The two stars are massive supergiants with distances that clearly place them in WLM. They are young ( less than or equal to 10 Myr) and should have a similar composition to the ISM. Additionally, differential abundances suggest that the O/Fe ratio in the WLM star is consistent with similar stars in NGC 6822 and the SMC, galaxies where the average stellar oxygen abundances are in excellent agreement with the nebular results. If the stellar abundances reflect the true composition of WLM, then this galaxy lies well above the metallicity-luminosity relationship for dwarf irregular galaxies. It also suggests that WLM is more chemically evolved than currently interpreted from its color-magnitude diagram. The similarities between the stars in WLM and NGC 6822 suggest that these two galaxies may have had similar star formation histories.

The VLT-FLAMES survey of massive stars: NGC 346-013 as a test case for massive close binary evolution

Context. NGC 346-013 is a peculiar double-lined eclipsing binary in the Small Magellanic Cloud (SMC) discovered by the VLT-FLAMES survey of massive stars.

Discovery of two distinct red clumps in NGC 419: a rare snapshot of a cluster at the onset of degeneracy

Colour-magnitude diagrams (CMDs) of the Small Magellanic Cloud (SMC) star cluster NGC 419, derived from Hubble Space Telescope (HST)/Advanced Camera for Surveys (ACS) data, reveal a well-delineated secondary clump located below the classical compact red clump typical of intermediate-age populations. We demonstrate that this feature belongs to the cluster itself, rather than to the underlying SMC field. Then, we use synthetic CMDs to show that it corresponds very well to the secondary clump predicted to appear as a result of He-ignition in stars just massive enough to avoid e(-)-degeneracy settling in their H-exhausted cores. The main red clump instead is made of the slightly less massive stars which passed through e(-) degeneracy and ignited He at the tip of the red giant branch. In other words, NGC 419 is the rare snapshot of a cluster while undergoing the fast transition from classical to degenerate H-exhausted cores. At this particular moment of a cluster`s life, the colour distance between the main-sequence turn-off and the red clump(s) depends sensitively on the amount of convective core overshooting, Lambda(c). By coupling measurements of this colour separation with fits to the red clump morphology, we are able to estimate simultaneously the cluster mean age (1.35(-0.04)(+0.11) Gyr) and overshooting efficiency (Lambda(c) = 0.47(-0.04)(+0.14)). Therefore, clusters like NGC 419 may constitute important marks in the age scale of intermediate-age populations. After eye inspection of other CMDs derived from HST/ACS data, we suggest that the same secondary clump may also be present in the Large Magellanic Cloud clusters NGC 1751, 1783, 1806, 1846, 1852 and 1917.

VELOCITY FIELD OF COMPRESSIBLE MAGNETOHYDRODYNAMIC TURBULENCE: WAVELET DECOMPOSITION AND MODE SCALINGS

We study compressible magnetohydrodynamic turbulence, which holds the key to many astrophysical processes, including star formation and cosmic-ray propagation. To account for the variations of the magnetic field in the strongly turbulent fluid, we use wavelet decomposition of the turbulent velocity field into Alfven, slow, and fast modes, which presents an extension of the Cho & Lazarian decomposition approach based on Fourier transforms. The wavelets allow us to follow the variations of the local direction of the magnetic field and therefore improve the quality of the decomposition compared to the Fourier transforms, which are done in the mean field reference frame. For each resulting component, we calculate the spectra and two-point statistics such as longitudinal and transverse structure functions as well as higher order intermittency statistics. In addition, we perform a Helmholtz-Hodge decomposition of the velocity field into incompressible and compressible parts and analyze these components. We find that the turbulence intermittency is different for different components, and we show that the intermittency statistics depend on whether the phenomenon was studied in the global reference frame related to the mean magnetic field or in the frame defined by the local magnetic field. The dependencies of the measures we obtained are different for different components of the velocity; for instance, we show that while the Alfven mode intermittency changes marginally with the Mach number, the intermittency of the fast mode is substantially affected by the change.

Reconnection studies under different types of turbulence driving

We study a model of fast magnetic reconnection in the presence of weak turbulence proposed by Lazarian and Vishniac (1999) using three-dimensional direct numerical simulations. The model has been already successfully tested in Kowal et al. (2009) confirming the dependencies of the reconnection speed V-rec on the turbulence injection power P-inj and the injection scale l(inj) expressed by a constraint V-rec similar to P(inj)(1/2)l(inj)(3/4)and no observed dependency on Ohmic resistivity. In Kowal et al. (2009), in order to drive turbulence, we injected velocity fluctuations in Fourier space with frequencies concentrated around k(inj) = 1/l(inj), as described in Alvelius (1999). In this paper, we extend our previous studies by comparing fast magnetic reconnection under different mechanisms of turbulence injection by introducing a new way of turbulence driving. The new method injects velocity or magnetic eddies with a specified amplitude and scale in random locations directly in real space. We provide exact relations between the eddy parameters and turbulent power and injection scale. We performed simulations with new forcing in order to study turbulent power and injection scale dependencies. The results show no discrepancy between models with two different methods of turbulence driving exposing the same scalings in both cases. This is in agreement with the Lazarian and Vishniac (1999) predictions. In addition, we performed a series of models with varying viscosity nu. Although Lazarian and Vishniac (1999) do not provide any prediction for this dependence, we report a weak relation between the reconnection speed with viscosity, V-rec similar to nu(-1/4).

DYNAMICAL EVOLUTION OF VISCOUS DISKS AROUND Be STARS. I. PHOTOMETRY

Be stars possess gaseous circumstellar disks that modify in many ways the spectrum of the central B star. Furthermore, they exhibit variability at several timescales and for a large number of observables. Putting the pieces together of this dynamical behavior is not an easy task and requires a detailed understanding of the physical processes that control the temporal evolution of the observables. There is an increasing body of evidence that suggests that Be disks are well described by standard alpha-disk theory. This paper is the first of a series that aims at studying the possibility of inferring several disk and stellar parameters through the follow-up of various observables. Here we study the temporal evolution of the disk density for different dynamical scenarios, including the disk buildup as a result of a long and steady mass injection from the star, the disk dissipation that occurs after mass injection is turned off, as well as scenarios in which active periods are followed by periods of quiescence. For those scenarios, we investigate the temporal evolution of continuum photometric observables using a three-dimensional non-LTE radiative transfer code. We show that light curves for different wavelengths are specific of a mass loss history, inclination angle, and alpha viscosity parameter. The diagnostic potential of those light curves is also discussed.

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.

The VLT-FLAMES survey of massive stars: atmospheric parameters and rotational velocity distributions for B-type stars in the Magellanic Clouds

Aims.We aim to provide the atmospheric parameters and rotational velocities for a large sample of O- and early B-type stars, analysed in a homogeneous and consistent manner, for use in constraining theoretical models. Methods: Atmospheric parameters, stellar masses, and rotational velocities have been estimated for approximately 250 early B-type stars in the Large (LMC) and Small (SMC) Magellanic Clouds from high-resolution VLT-FLAMES data using the non-LTE TLUSTY model atmosphere code. This data set has been supplemented with our previous analyses of some 50 O-type stars (Mokiem et al. 2006, 2007) and 100 narrow-lined early B-type stars (Hunter et al. 2006; Trundle et al. 2007) from the same survey, providing a sample of ~400 early-type objects. Results: Comparison of the rotational velocities with evolutionary tracks suggests that the end of core hydrogen burning occurs later than currently predicted and we argue for an extension of the evolutionary tracks. We also show that the large number of the luminous blue supergiants observed in the fields are unlikely to have directly evolved from main-sequence massive O-type stars as neither their low rotational velocities nor their position on the H-R diagram are predicted. We suggest that blue loops or mass-transfer binary systems may populate the blue supergiant regime. By comparing the rotational velocity distributions of the Magellanic Cloud stars to a similar Galactic sample, we find that (at 3s confidence level) massive stars (above 8 M?) in the SMC rotate faster than those in the solar neighbourhood. However there appears to be no significant difference between the rotational velocity distributions in the Galaxy and the LMC. We find that the v sin i distributions in the SMC and LMC can modelled with an intrinsic rotational velocity distribution that is a Gaussian peaking at 175 km s-1 (SMC) and 100 km s-1 (LMC) with a 1/e half width of 150 km s-1. We find that in NGC 346 in the SMC, the 10-25 M? main-sequence stars appear to rotate faster than their higher mass counterparts. It is not expected that O-type stars spin down significantly through angular momentum loss via stellar winds at SMC metallicity, hence this could be a reflection of mass dependent birth spin rates. Recently Yoon et al. (2006) have determined rates of GRBs by modelling rapidly rotating massive star progenitors. Our measured rotational velocity distribution for the 10-25 M? stars is peaked at slightly higher velocities than they assume, supporting the idea that GRBs could come from rapid rotators with initial masses as low as 14 M? at low metallicities.

Herschel key program HERITAGE:A far-infrared source catalog for the Magellanic Clouds

Observations from the HERschel Inventory of the Agents of Galaxy Evolution (HERITAGE ) have been used to identify dusty populations of sources in the Large and Small Magellanic Clouds (LMC and SMC). We conducted the study using the HERITAGE catalogs of point sources available from the Herschel Science Center from both the Photodetector Array Camera and Spectrometer (PACS; 100 and 160 μm) and Spectral and Photometric Imaging Receiver (SPIRE; 250, 350, and 500 μm) cameras. These catalogs are matched to each other to create a Herschel band-merged catalog and then further matched to archival Spitzer IRAC and MIPS catalogs from the Spitzer Surveying the Agents of Galaxy Evolution (SAGE) and SAGE-SMC surveys to create single mid- to far-infrared (far-IR) point source catalogs that span the wavelength range from 3.6 to 500 μm. There are 35,322 unique sources in the LMC and 7503 in the SMC. To be bright in the FIR, a source must be very dusty, and so the sources in the HERITAGE catalogs represent the dustiest populations of sources. The brightest HERITAGE sources are dominated by young stellar objects (YSOs), and the dimmest by background galaxies. We identify the sources most likely to be background galaxies by first considering their morphology (distant galaxies are point-like at the resolution of Herschel) and then comparing the flux distribution to that of the Herschel Astrophysical Terahertz Large Area Survey (ATLAS ) survey of galaxies. We find a total of 9745 background galaxy candidates in the LMC HERITAGE images and 5111 in the SMC images, in agreement with the number predicted by extrapolating from the ATLAS flux distribution. The majority of the Magellanic Cloud-residing sources are either very young, embedded forming stars or dusty clumps of the interstellar medium. Using the presence of 24 μm emission as a tracer of star formation, we identify 3518 YSO candidates in the LMC and 663 in the SMC. There are far fewer far-IR bright YSOs in the SMC than the LMC due to both the SMC's smaller size and its lower dust content. The YSO candidate lists may be contaminated at low flux levels by background galaxies, and so we differentiate between sources with a high ("probable") and moderate ("possible ") likelihood of being a YSO. There are 2493/425 probable YSO candidates in the LMC/SMC. Approximately 73% of the Herschel YSO candidates are newly identified in the LMC, and 35% in the SMC. We further identify a small population of dusty objects in the late stages of stellar evolution including extreme and post-asymptotic giant branch, planetary nebulae, and supernova remnants. These populations are identified by matching the HERITAGE catalogs to lists of previously identified objects in the literature. Approximately half of the LMC sources and one quarter of the SMC sources are too faint to obtain accurate ample FIR photometry and are unclassified.