The evolution of massive clumps in star forming regions

In this thesis two related arguments are investigated: - The first stages of the process of massive star formation, investigating the physical conditions and -properties of massive clumps in different evolutionary stages, and their CO depletion; - The influence that high-mass stars have on the nearby material and on the activity of star formation. I characterise the gas and dust temperature, mass and density of a sample of massive clumps, and analyse the variation of these properties from quiescent clumps, without any sign of active star formation, to clumps likely hosting a zero-age main sequence star. I briefly discuss CO depletion and recent observations of several molecular species, tracers of Hot Cores and/or shocked gas, of a subsample of these clumps. The issue of CO depletion is addressed in more detail in a larger sample consisting of the brightest sources in the ATLASGAL survey: using a radiative tranfer code I investigate how the depletion changes from dark clouds to more evolved objects, and compare its evolution to what happens in the low-mass regime. Finally, I derive the physical properties of the molecular gas in the photon-dominated region adjacent to the HII region G353.2+0.9 in the vicinity of Pismis 24, a young, massive cluster, containing some of the most massive and hottest stars known in our Galaxy. I derive the IMF of the cluster and study the star formation activity in its surroundings. Much of the data analysis is done with a Bayesian approach. Therefore, a separate chapter is dedicated to the concepts of Bayesian statistics.

A Stellar Census of the Tucana-Horologium Moving Group

We report the selection and spectroscopic confirmation of 129 new late-type (SpT = K3-M6) members of the Tucana-Horologium moving group, a nearby (d similar to 40 pc), young (tau similar to 40 Myr) population of comoving stars. We also report observations for 13 of the 17 known Tuc-Hor members in this spectral type range, and that 62 additional candidates are likely to be unassociated field stars; the confirmation frequency for new candidates is therefore 129/191 = 67%. We have used radial velocities, Ha emission, and Li-6708 absorption to distinguish between contaminants and bona fide members. Our expanded census of Tuc-Hor increases the known population by a factor of similar to 3 in total and by a factor of similar to 8 for members with SpT >= K3, but even so, the K-M dwarf population of Tuc-Hor is still markedly incomplete. Our expanded census allows for a much more detailed study of Tuc-Hor than was previously feasible. The spatial distribution of members appears to trace a two-dimensional sheet, with a broad distribution in X and Y, but a very narrow distribution (+/- 5 pc) in Z. The corresponding velocity distribution is very small, with a scatter of +/- 1.1 km s(-1) about the mean UVW velocity for stars spanning the entire 50 pc extent of Tuc-Hor. We also show that the isochronal age (tau similar to 20-30 Myr) and the lithium depletion boundary age (tau similar to 40 Myr) disagree, following the trend in other pre-main-sequence populations for isochrones to yield systematically younger ages. The H alpha emission line strength follows a trend of increasing equivalent width with later spectral type, as is seen for young clusters. We find that moving group members have been depleted of measurable lithium for spectral types of K7.0-M4.5. None of our targets have significant infrared excesses in the WISE W3 band, yielding an upper limit on warm debris disks of F < 0.7%. Finally, our purely kinematic and color-magnitude selection procedure allows us to test the efficiency and completeness for activity-based selection of young stars. We find that 60% of K-M dwarfs in Tuc-Hor do not have ROSAT counterparts and would have been omitted in X-ray-selected samples. In contrast, GALEX UV-selected samples using a previously suggested criterion for youth achieve completeness of 77% and purity of 78%, and we suggest new SpT-dependent selection criteria that will yield > 95% completeness for tau similar to 40 Myr populations with GALEX data available.

Is there a Population II analogy to the F star lithium dip?

Observers have found a small number of lithium-depleted halo stars in the temperature range of the Spite plateau. The current status of the mass-loss hypothesis for producing the observed lithium dip in Population (Pop) I stars is briefly discussed and extended to Pop II stars as a possible explanation for these halo objects. Based on detections of F-type main-sequence variables, mass loss is assumed to occur in a narrow temperature region corresponding to this “instability strip.” As Pop II main-sequence stars evolve to the blue, they enter this narrow temperature region, then move back through the lower temperature area of the Spite plateau. If 0.05 M⊙ (solar mass) or more have been lost, they will show lithium depletion. This hypothesis affects the lithium-to- beryllium abundance, the ratio of high- to low-lithium stars, and the luminosity function. Constraints on the mass-loss hypothesis due to these effects are discussed. Finally, mass loss in this temperature range would operate in stars near the turnoff of metal-poor globular clusters, resulting in apparent ages 2 to 3 Gyr (gigayears) older than they actually are.

Pathways to quiescence: SHARDS view on the star formation histories of massive quiescent galaxies at 1.0 < z < 1.5

We present star formation histories (SFHs) for a sample of 104 massive (stellar mass M > 10^10 M_⊙) quiescent galaxies (MQGs) at z = 1.0–1.5 from the analysis of spectrophotometric data from the Survey for High-z Absorption Red and Dead Sources (SHARDS) and HST/WFC3 G102 and G141 surveys of the GOODS-North field, jointly with broad-band observations from ultraviolet (UV) to far-infrared (far-IR). The sample is constructed on the basis of rest-frame UVJ colours and specific star formation rates (sSFRs = SFR/Mass). The spectral energy distributions (SEDs) of each galaxy are compared to models assuming a delayed exponentially declining SFH. A Monte Carlo algorithm characterizes the degeneracies, which we are able to break taking advantage of the SHARDS data resolution, by measuring indices such as MgUV and D4000. The population of MQGs shows a duality in their properties. The sample is dominated (85 per cent) by galaxies with young mass-weighted ages, t_M t_M < 2 Gyr, short star formation time-scales, 〈τ〉 ∼ 60–200 Myr, and masses log(M/M_⊙) ∼ 10.5. There is an older population (15 per cent) with t_M t_M = 2–4 Gyr, longer star formation time-scales, 〈τ〉∼ 400 Myr, and larger masses, log(M/M_⊙) ∼ 10.7. The SFHs of our MQGs are consistent with the slope and the location of the main sequence of star-forming galaxies at z > 1.0, when our galaxies were 0.5–1.0 Gyr old. According to these SFHs, all the MQGs experienced a luminous infrared galaxy phase that lasts for ∼500 Myr, and half of them an ultraluminous infrared galaxy phase for ∼100 Myr. We find that the MQG population is almost assembled at z ∼ 1, and continues evolving passively with few additions to the population.

No surviving evolved companions of the progenitor of SN 1006

Type Ia supernovae are thought to occur when a white dwarf made of carbon and oxygen accretes sufficient mass to trigger a thermonuclear explosion(1). The accretion could be slow, from an unevolved (main-sequence) or evolved (subgiant or giant) star(2,3) (the single-degenerate channel), or rapid, as the primary star breaks up a smaller orbiting white dwarf(3,4) (the double-degenerate channel). A companion star will survive the explosion only in the single-degenerate channel(5). Both channels might contribute to the production of type Ia supernovae(6,7), but the relative proportions of their contributions remain a fundamental puzzle in astronomy. Previous searches for remnant companions have revealed one possible case for SN 1572 (refs 8, 9), although that has been questioned(10). More recently, observations have restricted surviving companions to be small, main-sequence stars(11-13), ruling out giant companions but still allowing the single-degenerate channel. Here we report the results of a search for surviving companions of the progenitor of SN 1006 (ref. 14). None of the stars within 4 arc minutes of the apparent site of the explosion is associated with the supernova remnant, and we can firmly exclude all giant and subgiant stars from being companions of the progenitor. In combination with previous results, our findings indicate that fewer than 20 per cent of type Ia supernovae occur through the single-degenerate channel.

Post T Tauri stars in the solar neighborhood: isolated or members of young associations and moving groups

Post T Tauri stars (PTTS) are late-type stars in the age range between 10 and 100 Myr filling the gap between T Tauri (TTs) and zero-age: main sequence phases. This period of evolution remains ambiguous and until now different studies of young stars have failed to find the numbers of PTTS that are expected. In the last years, some PTTS have been identified among the X-ray detected pre-main sequence stars in some star-forming regions. More recently, additional PTTS have been identified in young associations and moving groups (β Pic, TW Hya, Tucana/Horologium, and the AB Dor). However, many isolated PTTS still remain undiscovered. In this contribution, we compiled the PTTS previously identified in the literature, and identified new candidates using the information provided by the high resolution spectra obtained during our surveys of late-type stars possible members to young moving groups, FGK stars in the solar neighborhood, and RasTyc sample. To identify PTTS we applied an age-oriented definition using relative age indicators (Li abundance, chromospheric and coronal emission and the kinematics) as well as color-magnitude diagrams and pre-main sequence isochrones.

Hα emission fluxes and lithium abundances of low-mass stars in the young open cluster IC 4665

As part of a long term effort to understand pre-main sequence Li burning, we have obtained high resolution spectroscopic observations of 14 late type stars (G0-M1) in the young open cluster IC 4665. Most of the stars have Hα filled-in and Li I absorption, as expected for their young age. From the equivalent widths of Hα emission excess (obtained using the spectral subtraction technique) and the the Li i λ6708 feature, we have derived Hα emission fluxes and photospheric Li abundances. The mean Li abundance of IC 4665 solar-type stars is log N(Li) = 3.1; the same as in other young clusters (α Per, Pleiades) and T Tauri stars. Our results support the conclusions from previous works that PMS Li depletion is very small for masses ∼ 1 M_⨀ . Among the IC 4665 late-G and early K-type stars, there is a spread in Li abundances of about one order of magnitude. The Li-poor IC 4665 members have low Hα excess and vsini≤10. Hence, the Li-activity-rotation connection which has been clearly established in the Pleiades also seems to hold in IC 4665. One M-type IC 4665 star that we have observed does not show Li, implying a very efficient Li depletion as observed in α Per stars of the same spectral type. The level of chromospheric activity and Li depletion among the low-mass stars of IC 4665 is similar to that in the Pleiades. In fact, we note that the Li abundance distributions in several young clusters (α Per, Pleiades, IC 2391, IC 4665) and in post T Tauri stars are strikingly similar. This result suggests that Hα emission and Li abundance not well correlated with age for low-mass stars between 20 and 100 Myr old. We argue that a finer age indicator, the ''LL-clock'', would be the luminosity at which the transition between efficient Li depletion and preservation takes place for fully convective objects. The LL-clock could allow in the near future to derive the relative ages of young open clusters, and clarify the study of PMS evolution of cool stars.

Stratigraphic and geochemical study of the organic-rich black shales in the Tarcău Nappe of the Moldavidian Domain (Carpathian Chain, Romania)

An integrated stratigraphic analysis has been made of the Tarcău Nappe (Moldavidian Domain, Eastern Romanian Carpathians), coupled with a geochemical study of organic-rich beds. Two Main Sequence Boundaries (Early Oligocene and near to the Oligocene–Aquitanian boundary, respectively) divide the sedimentary record into three depositional sequences. The sedimentation occurred in the central area of a basin supplied by different and opposite sources. The high amount of siliciclastics at the beginning of the Miocene marks the activation of the “foredeep stage”. The successions studied are younger than previously thought and they more accurately date the deformation of the different Miocene phases affecting the Moldavidian Basin. The intervals with black shales identified are related to two main separate anoxic episodes with an age not older than Late Rupelian and not before Late Chattian. The most important organic-rich beds correspond to the Lower Menilites, Bituminous Marls and Lower Dysodilic Shales Members (Interval 2). These constitute a good potential source rock for petroleum, with homogeneous Type II oil-prone organic matter, highly lipidic and thermally immature. The deposition of black shales has been interpreted as occurring within a deep, periodically isolated and tectonically controlled basin.

An accurate distance scale to the extremely massive cluster Stephenson 2

Discoveries during the last two years have revealed the existence of a vast region of star formation close to the base of the Scutum Arm, where at least five clusters of red supergiants have been found. In order to understand the nature of this region, we need to determine accurate distances to the clusters. We present here the first results of an ongoing program to derive fundamental parameters (such as age, distance, etc.) to the massive cluster Stephenson 2 studying for the first time its main sequence stars.

A VLT/FLAMES survey for massive binaries in Westerlund 1. III. The WC9d binary W239 and implications for massive stellar evolution

Context. There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. However the binary properties of massive stars – frequency, mass ratio & orbital separation – are still poorly constrained. Aims. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate binary W239. Methods. To accomplish this we have utilised our spectroscopic data in conjunction with multi-year optical and near-IR photometric observations in order to search for binary signatures. Comparison of these data to synthetic non-LTE model atmosphere spectra were used to derive the fundamental properties of the WC9 primary. Results. We found W239 to have an orbital period of only ~5.05 days, making it one of the most compact WC binaries yet identified. Analysis of the long term near-IR lightcurve reveals a significant flare between 2004-6. We interpret this as evidence for a third massive stellar component in the system in a long period (>6 yr), eccentric orbit, with dust production occuring at periastron leading to the flare. The presence of a near-IR excess characteristic of hot (~1300 K) dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short (<1 yr) period binaries, although confirmation will require further observations. Non-LTE model atmosphere analysis of the spectrum reveals the physical properties of the WC9 component to be fully consistent with other Galactic examples. Conclusions. The simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides compelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period O+OB binaries will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss – likely via case A mass transfer or a contact configuration – or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar evolution up to and beyond supernova and the resultant production of relativistic remnants.

NGC 7419 as a template for red supergiant clusters

Context. The open cluster NGC 7419 is known to contain five red supergiants and a very high number of Be stars. However, there are conflicting reports about its age and distance that prevent a useful comparison with other clusters. Aims. We intend to obtain more accurate parameters for NGC 7419, using techniques different from those of previous authors, so that it may be used as a calibrator for more obscured clusters. Methods. We obtained Strömgren photometry of the open cluster NGC 7419, as well as classification spectroscopy of ~20 stars in the area. We then applied standard analysis and classification techniques. Results. We find a distance of 4 ± 0.4 kpc and an age of 14 ± 2 Myr for NGC 7419. The main-sequence turn-off is found at spectral type B1, in excellent agreement. We identify 179 B-type members, implying that there are more than 1200 M⊙ in B stars at present. Extrapolating this to lower masses indicates an initial cluster mass of between 7000 and 10 000 M⊙, depending on the shape of the initial mass function. We find a very high fraction (≈40%) of Be stars around the turn-off, but very few Be stars at lower masses. We also report for the first time a strong variability in the emission characteristics of Be stars. We verified that the parameters of the red supergiant members can be used to obtain accurate cluster parameters. Conclusions. NGC 7419 is sufficiently massive to serve as a testbed for theoretical predictions and as a template to compare more obscured clusters. The distribution of stars above the main-sequence turn-off is difficult to accommodate with current evolutionary tracks. Though the presence of five red supergiants is marginally consistent with theoretical expectations, the high number of Be stars and very low number of luminous evolved B stars hint at some unknown physical factor that is not considered in current synthesis models.

The circumstellar environment and evolutionary status of the supergiant B[e] star Wd1-9

Context. Historically, supergiant (sg)B[e] stars have been difficult to include in theoretical schemes for the evolution of massive OB stars. Aims. The location of Wd1-9 within the coeval starburst cluster Westerlund 1 means that it may be placed into a proper evolutionary context and we therefore aim to utilise a comprehensive multiwavelength dataset to determine its physical properties and consequently its relation to other sgB[e] stars and the global population of massive evolved stars within Wd1. Methods. Multi-epoch R- and I-band VLT/UVES and VLT/FORS2 spectra are used to constrain the properties of the circumstellar gas, while an ISO-SWS spectrum covering 2.45−45μm is used to investigate the distribution, geometry and composition of the dust via a semi-analytic irradiated disk model. Radio emission enables a long term mass-loss history to be determined, while X-ray observations reveal the physical nature of high energy processes within the system. Results. Wd1-9 exhibits the rich optical emission line spectrum that is characteristic of sgB[e] stars. Likewise its mid-IR spectrum resembles those of the LMC sgB[e] stars R66 and 126, revealing the presence of equatorially concentrated silicate dust, with a mass of ~10−4M⊙. Extreme historical and ongoing mass loss (≳ 10−4M⊙yr−1) is inferred from the radio observations. The X-ray properties of Wd1-9 imply the presence of high temperature plasma within the system and are directly comparable to a number of confirmed short-period colliding wind binaries within Wd1. Conclusions. The most complete explanation for the observational properties of Wd1-9 is that it is a massive interacting binary currently undergoing, or recently exited from, rapid Roche-lobe overflow, supporting the hypothesis that binarity mediates the formation of (a subset of) sgB[e] stars. The mass loss rate of Wd1-9 is consistent with such an assertion, while viable progenitor and descendent systems are present within Wd1 and comparable sgB[e] binaries have been identified in the Galaxy. Moreover, the rarity of sgB[e] stars - only two examples are identified from a census of ~ 68 young massive Galactic clusters and associations containing ~ 600 post-Main Sequence stars - is explicable given the rapidity (~ 104yr) expected for this phase of massive binary evolution.

MY Camelopardalis, a very massive merger progenitor

Context. The early-type binary MY Cam belongs to the young open cluster Alicante 1, embedded in Cam OB3. Aims. MY Cam consists of two early-O type main-sequence stars and shows a photometric modulation suggesting an orbital period slightly above one day. We intend to confirm this orbital period and derive orbital and stellar parameters. Methods. Timing analysis of a very exhaustive (4607 points) light curve indicates a period of 1.1754514 ± 0.0000015 d. High-resolution spectra and the cross-correlation technique implemented in the todcor program were used to derive radial velocities and obtain the corresponding radial velocity curves for MY Cam. Modelling with the stellar atmosphere code fastwind was used to obtain stellar parameters and create templates for cross-correlation. Stellar and orbital parameters were derived using the Wilson-Devinney code, such that a complete solution to the binary system could be described. Results. The determined masses of the primary and secondary stars in MY Cam are 37.7 ± 1.6 and 31.6 ± 1.4M⊙, respectively. The corresponding temperatures, derived from the model atmosphere fit, are 42 000 and 39 000 K, with the more massive component being hotter. Both stars are overfilling their Roche lobes, sharing a common envelope. Conclusions. MY Cam contains the most massive dwarf O-type stars found so far in an eclipsing binary. Both components are still on the main sequence, and probably not far from the zero-age main sequence. The system is a likely merger progenitor, owing to its very short period.

On the formation and evolution of the first Be star in a black hole binary MWC 656

We find that the formation of MWC 656 (the first Be binary containing a black hole) involves a common envelope phase and a supernova explosion. This result supports the idea that a rapidly rotating Be star can emerge out of a common envelope phase, which is very intriguing because this evolutionary stage is thought to be too fast to lead to significant accretion and spin up of the B star. We predict ∼10–100 of B-BH binaries to currently reside in the Galactic disc, among which around 1/3 contain a Be star, but there is only a small chance to observe a system with parameters resembling MWC 656. If MWC 656 is representative of intrinsic Galactic Be-BH binary population, it may indicate that standard evolutionary theory needs to be revised. This would pose another evolutionary problem in understanding black hole (BH) binaries, with BH X-ray novae formation issue being the prime example. Future evolution of MWC 656 with an ∼5 M⊙ BH and with an ∼13 M⊙ main-sequence companion on an ∼60 d orbit may lead to the formation of a coalescing BH–NS (neutron star) system. The estimated Advanced LIGO/Virgo detection rate of such systems is up to ∼0.2 yr−1. This empirical estimate is a lower limit as it is obtained with only one particular evolutionary scenario, the MWC 656 binary. This is only a third such estimate available (after Cyg X-1 and Cyg X-3), and it lends additional support to the existence of so far undetected BH–NS binaries.

GU Monocerotis: A high-mass eclipsing overcontact binary in the young open cluster Dolidze 25

Context. The eclipsing binary GU Mon is located in the star-forming cluster Dolidze 25, which has the lowest metallicity measured in a Milky Way young cluster. Aims. GU Mon has been identified as a short-period eclipsing binary with two early B-type components. We set out to derive its orbital and stellar parameters. Methods. We present a comprehensive analysis, including B and V light curves and 11 high-resolution spectra, to verify the orbital period and determine parameters. We used the stellar atmosphere code FASTWIND to obtain stellar parameters and create templates for cross-correlation. We obtained a model to fit the light and radial-velocity curves using the Wilson-Devinney code iteratively and simultaneously. Results. The two components of GU Mon are identical stars of spectral type B1 V with the same mass and temperature. The light curves are typical of an EW-type binary. The spectroscopic and photometric analyses agree on a period of 0.896640 ± 0.000007 d. We determine a mass of 9.0 ± 0.6 M⊙ for each component and for temperatures of 28 000 ± 2000 K. Both values are consistent with the spectral type. The two stars are overfilling their respective Roche lobes, sharing a common envelope and, therefore the orbit is synchronised and circularised. Conclusions. The GU Mon system has a fill-out factor above 0.8, containing two dwarf B-type stars on the main sequence. The two stars are in a very advanced stage of interaction, with their extreme physical similarity likely due to the common envelope. The expected evolution of such a system very probably leads to a merger while still on the main sequence.