Search

The Villafranca project is studying Galactic stellar groups with OB stars combining information from $Gaia$ and ground-based surveys. We summarize the status of the project and we present its most important results. The Villafranca project has been used to produce a new astrometric calibration for $Gaia$ (E)DR3, which improves the previous one significantly for bright stars. We have discovered that dynamical interactions among massive stars at a very young age ($\sim$1 Ma or less) can play a significant interaction in the dynamical evolution of clusters. As a consequence, our current view of the massive-star IMF may be distorted and the number of free-floating neutron stars and black holes higher than previously considered., Comment: To appear in Highlights of Spanish Astrophysics XIII, P254, version with updated style file

Massive stars (those larger than 8 solar masses at formation) have radiative envelopes that cannot sustain a dynamo, the mechanism that produces magnetic fields in lower-mass stars. Despite this, approximately 7\% of massive stars have observed magnetic fields, the origin of which is debated. We used multi-epoch interferometric and spectroscopic observations to characterize HD 148937, a binary system of two massive stars. We found that only one star is magnetic and that it appears younger than its companion. The system properties and a surrounding bipolar nebula can be reproduced with a model in which two stars merged (in a previous triple system) to produce the magnetic massive star. Our results provide observational evidence that magnetic fields form in at least some massive stars through stellar mergers., Comment: Full paper with supplementary materials. 59 pages, 18 figures

Stellar rotation is of key importance for the formation process, evolution, and final fate of massive stars. In this paper we review results from the study of the spin rate properties of a sample of more than 400 Galactic O-type stars surveyed by the IACOB and OWN projects. By combining vsini, Teff, and logg estimates (resulting from a detailed quantitative spectroscopic analysis) with information about the spectroscopic binarity status for an important fraction of the stars in the sample, we provide a renewed overview about how the empirical distribution of projected rotational velocities in the O-star domain depends on mass, evolutionary and binary status. The obtained distributions are then compared with predictions of several state-of-the-art evolutionary models for single stars, as well as from population synthesis simulations including binary interaction, and used to provide hints about the initial velocity distribution of stars with masses in the range ~15-80 Msol., Comment: 24 pages, 19 figures, accepted for publication in Astronomy & Astrophysics

At the time of this meeting, the latest Gaia data release is EDR3, published on 3 December 2020, but the next one, DR3, will appear soon, on 13 June 2022. This contribution describes, on the one hand, Gaia EDR3 results on massive stars and young stellar clusters, placing special emphasis on how a correct treatment of the astrometric and photometric calibration yields results that are simultaneously precise and accurate. On the other hand, it gives a brief description of the exciting results we can expect from Gaia DR3., Comment: To appear in Massive Stars Near and Far, Proceedings IAU Symposium No. 361, 2022. N. St-Louis, J. S. Vink & J. Mackey, eds

Multiplicity is a ubiquitous characteristic of massive stars. Multiple systems offer us a unique observational constraint on the formation of high-mass systems. Herschel 36 A is a massive triple system composed of a close binary (Ab1-Ab2) and an outer component (Aa). We measured the orbital motion of the outer component of Herschel 36 A using infrared interferometry with the AMBER and PIONIER instruments of ESO's Very Large Telescope Interferometer. Our immediate aims are to constrain the masses of all components of this system and to determine if the outer orbit is co-planar with the inner one. Reported spectroscopic data for all three components of this system and our interferometric data allow us to derive full orbital solutions for the outer orbit Aa-Ab and the inner orbit Ab1-Ab2. For the first time, we derive the absolute masses of mAa = 22.3 +/- 1.7 M_sun, mAb1 = 20.5 +/- 1.5 M_sun and mAb2 = 12.5 +/- 0.9 M_sun. Despite not being able to resolve the close binary components, we infer the inclination of their orbit by imposing the same parallax as the outer orbit. Inclinations derived from the inner and outer orbits imply a modest difference of about 22 deg. between the orbital planes. We discuss this result and the formation of Herschel 36 A in the context of Core Accretion and Competitive Accretion models, which make different predictions regarding the statistic of the relative orbital inclinations., Comment: accepted to be published by MNRAS, 8 pages, 5 figures

We use Gaia EDR3 astrometry to propose that a dynamical interaction between the multiple system $\theta^1$ Ori C and $\theta^1$ Ori F ejected the latter as a walkaway star ~1100 years ago (without deceleration) or somewhat later (with a more likely deceleration included). It is unclear whether the final 3-D velocity of $\theta^1$ Ori F will be large enough to escape the Orion nebula cluster., Comment: Accepted in RNAAS

CONTEXT. Dynamical interactions in young stellar clusters can eject massive stars early in their lives and significantly alter their mass functions. If all of the most massive stars are lost, we are left with an orphan cluster. AIMS. We study the Bermuda cluster (Villafranca O-014 NW), the most significant young stellar group in the North America and Pelican nebulae, and the massive stars that may have been ejected from it to test if it has been orphaned. METHODS. We use Gaia EDR3 parallaxes and proper motions to search for walkaway/runaway stars in the vicinity of the North America and Pelican nebulae. The candidates are analyzed with spectroscopy and photometry to assess their nature and their trajectories are traced back in time to determine at what time they left the Bermuda cluster. RESULTS. We detect three ejection events (Bajamar, Toronto, and HD 201 795 events) that expelled 5, 2, and 2 systems, respectively, or 6, 3, and 3 stars if we count the individual components in spectroscopic/eclipsing binaries. The events took place 1.611$\pm$0.011 Ma, 1.496$\pm$0.044 Ma, and 1.905$\pm$0.037 Ma ago, respectively, but our analysis is marginally consistent with the first two being simultaneous. We detect bow shocks in WISE images associated with four of the ejected systems. Combining the three events, the Bermuda cluster has lost >200 M_Sol, including its three most massive stars, so it can be considered an orphan cluster. One consequence is that the PDMF of the cluster has been radically altered from its top-heavy initial value to one compatible with a Kroupa-like function. Another one is that the cluster is currently expanding with a dynamical time scale consistent with the cause being the ejection events. A scenario in which the Bermuda cluster was formed in a conveyor belt fashion over several hundreds of ka or even 1 Ma is consistent with all the observables. [ABRIDGED], Comment: Accepted for publication in A&A. New version after acceptance with some added content suggested by referee

CONTEXT. This is the second paper of a series on Galactic OB groups that uses astrometric and photometric data from Gaia and spectral classifications from GOSSS and LiLiMaRlin. The previous paper was based on Gaia DR2, this is based on Gaia EDR3. AIMS. The two aims of this paper are to revise the results for the sample from paper I using Gaia EDR3 data and to expand the sample of analyzed stellar groups to 26, from Villafranca O-001 to Villafranca O-026. METHODS. We use GOSSS to select Galactic stellar groups with O stars and an updated version of the method in paper 0 of this series, combining Gaia EDR3 G+GBP+GRP photometry, positions, proper motions, and parallaxes to assign memberships and measure distances. We present 99 spectra from GOSSS and 32 from LiLiMaRlin for stars in the analyzed groups or in their foreground. RESULTS. We derive distances to the 26 stellar groups with unprecedented precision and accuracy, with total uncertainties <1% within 1 kpc and of ~3% around 3 kpc, values that are almost 4x better than for Gaia DR2. We provide homogeneous spectral types for 110 stars and correct a number of errors in the literature, especially for objects in the Orion nebula cluster. For each of the groups we discuss its membership and present possible runaway/walkaway stars. At least two groups, Villafranca O-O12 S and Villafranca O-014 NW, are orphan clusters in which the most massive stars have been ejected by dynamical interactions, leaving objects with a capped mass function. The existence of such clusters has important consequences for the study of the IMF, the distribution of SNe across the Galaxy, and the population and dynamics of isolated compact objects. We fit PMS isochrones to four clusters to derive ages of 2.0$\pm$0.5 Ma for the sigma Orionis cluster, 4$\pm$2 Ma for NGC 2264, 5.0$\pm$0.5 Ma for NGC 2362, and 8$\pm$2 Ma for the gamma Velorum cluster. [ABRIDGED], Comment: Accepted for publication in A&A. Original arXiv version was submitted prior to paper being accepted

The GALANTE optical photometric survey is observing the northern Galactic plane and some adjacent regions using seven narrow- and intermediate-filters, covering a total of 1618 square degrees. The survey has been designed with multiple exposure times and at least two different air masses per field to maximize its photometric dynamic range, comparable to that of Gaia, and ensure the accuracy of its photometric calibration. The goal is to reach at least 1% accuracy and precision in the seven bands for all stars brighter than AB magnitude 17 while detecting fainter stars with lower values of the signal-to-noise ratio.The main purposes of GALANTE are the identification and study of extinguished O+B+WR stars, the derivation of their extinction characteristics, and the cataloguing of F and G stars in the solar neighbourhood. Its data will be also used for a variety of other stellar studies and to generate a high-resolution continuum-free map of the H{\alpha} emission in the Galactic plane. We describe the techniques and the pipeline that are being used to process the data, including the basis of an innovative calibration system based on Gaia DR2 and 2MASS photometry., Comment: 18 pages, accepted for publication in MNRAS

[ABRIDGED] AIMS: The MONOS project is collecting information and studying O-type spectroscopic binaries with delta > -20 deg. In this 2nd paper, we tackle the study of the 35 single-line spectroscopic binary (SB1) systems identified in the previous paper of the series (arXiv:1904.11385) by analyzing our data and reviewing the literature orbits of such systems. METHODS: We have measured the radial velocities for the ~700 spectra in our database using two different methods: Gaussian fitting for several diagnostic lines per object and cross-correlation using synthetic spectra. We also explored the TESS database and analyzed the light curves for 31 of the systems. RESULTS: We have confirmed 21 SB1 systems, discarded the binary nature of 6 stars (9 Sge, HD 192 281, HDE 229 232 AB, 68 Cyg, HD 108 and \alpha Cam), and left 6 stars as inconclusive due to lack of data. The remaining two stars are 15 Mon Aa which has been classified as SB2, and Cyg OB2-22 C, for which we find evidence that it is most likely a triple system where the O star is orbiting an eclipsing SB1. We have also recalculated 20 new orbital solutions, including the first spectroscopic orbital solution for V747 Cep. For Cyg OB2-22 C we have obtained new ephemerides but no new orbit., Comment: Table 10 and Appendix C are available at the CDS

R 144 is the visually brightest WR star in the Large Magellanic Cloud (LMC). R 144 was reported to be a binary, making it potentially the most massive binary thus observed. We perform a comprehensive spectral, photometric, orbital, and polarimetric analysis of R 144. R 144 is an eccentric (e=0.51) 74.2-d binary comprising two relatively evolved (age~2 Myr), H-rich WR stars. The hotter primary (WN5/6h, T=50 kK) and the cooler secondary (WN6/7h,T=45kK) have nearly equal masses. The combination of low rotation and H-depletion observed in the system is well reproduced by contemporary evolution models that include boosted mass-loss at the upper-mass end. The systemic velocity of R 144 and its relative isolation suggest that it was ejected as a runaway from the neighbouring R 136 cluster. The optical light-curve shows a clear orbital modulation that can be well explained as a combination of two processes: excess emission stemming from wind-wind collisions and double wind eclipses. Our light-curve model implies an orbital inclination of i=60.4+-1.5deg, resulting in accurately constrained dynamical masses of 74+-4 and 69+-4 Msun. Assuming that both binary components are core H-burning, these masses are difficult to reconcile with the derived luminosities (logL1,2 = 6.44, 6.39 [Lsun]), which correspond to evolutionary masses of the order of 110 and 100Msun, respectively. Taken at face value, our results imply that both stars have high classical Eddington factors of Gamma_e = 0.78+-0.1. If the stars are on the main sequence, their derived radii (~25Rsun) suggest that they are only slightly inflated, even at this high Eddington factor. Alternatively, the stars could be core-He burning, strongly inflated from the regular size of classical Wolf-Rayet stars (~1Rsun), a scenario that could help resolve the observed mass discrepancy., Comment: Accepted to A&A on the April 06 2021. 17 pages + 5 page appendix. Version 28.10.2021: corrected typo in Eq. (2); numerator and denominator were inversed on RHS (only typo; no impact on results)

We cross-match the Alma catalog of OB stars with Gaia DR2 astrometry and photometry as a first step towards producing a clean sample of massive stars in the solar neighbourhood with a high degree of completeness. We analyze the resulting colour-absolute magnitude diagram to divide our sample into categories and compare extinction estimates from two sources, finding problems with both of them. The distances obtained with three different priors are found to have few differences among them, indicating that Gaia DR2 distances are robust. An analysis of the 3-D distribution of massive stars in the solar neighbourhood is presented. We show that a kinematically distinct structure we dub the Cepheus spur extends from the Orion-Cygnus spiral arm towards the Perseus arm and is located above the Galactic mid-plane, likely being related to the recently discovered Radcliffe wave. We propose that this corrugation pattern in the Galactic disk may be responsible for the recent enhanced star formation at its crests and troughs. We also discuss our plans to extend this work in the immediate future., Comment: Accepted for publication in MNRAS

CONTEXT. Gaia EDR3 has produced parallaxes for 1.468x10^9 sources but there are calibration issues that require corrections to the published values and uncertainties. AIMS. We want to characterize the behavior of the uncertainties of the Gaia EDR3 parallaxes. We also aim to provide a procedure for the calculation of distances to stars and stellar clusters. METHODS. We reanalyze some of the data in the calibration papers for QSO and LMC parallaxes and combine those results with measurements for six bright GCs. We calculate the angular covariance of EDR3 parallaxes at small separations based on the LMC results and combine it with the results for larger angles using QSOs to obtain an analytical formula for the angular covariance over the whole sky. The results for the six GCs are used to validate the parallax bias correction as a function of magnitude, color, and ecliptic latitude and to determine the constant used to convert internal uncertainties to external ones. RESULTS. The angular covariance at zero separation is 106.2 muas^2, yielding a minimum uncertainty for EDR3 parallaxes of 10.3 muas for individual stars. That value can be only slightly reduced for GCs after considering the behavior of the angular covariance of the parallaxes for small separations. The Lindegren et al. parallax bias correction works quite well, except for the brighter magnitudes, suggesting improvements may be possible there. The value of k is 1.1-1.7 and depends on G. Stars with moderately large values of RUWE can still provide useful parallaxes albeit with larger values of k. We give accurate and precise Gaia EDR3 distances to the six GCs and for the specific case of 47 Tuc we are able to beat the angular covariance limit and derive a high-precision distance of 4.53+-0.06 kpc. Finally, a recipe for the derivation of distances to stars and stellar clusters using Gaia EDR3 parallaxes is given. [ABRIDGED], Comment: Version accepted in A&A

CONTEXT. Many massive stars have nearby companions whose presence hamper their characterization through spectroscopy. AIMS. We want to continue obtaining spatially resolved spectroscopy of close massive visual binaries to derive their spectral types. METHODS. We have used lucky spectroscopy to obtain many short long-slit spectroscopic exposures of 19 close visual binaries under good seeing conditions. We selected those with the best characteristics, extracted the spectra using multiple-profile fitting, and combined the results to derive spatially separated spectra. The results are analyzed in combination with data from lucky imaging, regular intermediate-resolution single-order spectroscopy, and \'echelle high-resolution spectroscopy. RESULTS. The new application of lucky spectroscopy has allowed us to [a] spatially disentangle for the first time two O stars (FN CMa B and 6 Cas B) with brighter BA supergiant companions; [b] determine that two B stars (alpha Sco B and HD 164492 B) with close and more massive companions are fast rotators; [c] extend the technique to cases with extreme magnitude differences, shorter separations, and fainter primary magnitudes down to B=11; [d] spatially disentangle the spectra of stars with companions as diverse as an A supergiant (6 Cas A), a WR star (WR 157), and an M supergiant (alpha Sco A); [e] discover the unexpected identity of some targets such as two previously unknown bright O stars (HD 51756 B and BD +60 544) and a new member of the rare OC category (HD 8768 A); and [f] identify and classify which of the components of four visual binaries is a double-lined spectroscopic binary and for another seven systems detect signs of spectroscopic binarity using high-spectral-resolution spectroscopy. We also present a determination of the limits of the technique. [ABRIDGED], Comment: Accepted for publication in A&A

Milky Way globular clusters (MW GCs) are difficult to identify at low Galactic latitudes because of high differential extinction and heavy star crowding. The new deep near-IR images and photometry from the VISTA Variables in the Via L\'actea Extended Survey (VVVX) allow us to chart previously unexplored regions. Our long term aim is to complete the census of MW GCs. The immediate goals are to estimate the astrophysical parameters, measuring their reddenings, extinctions, distances, total luminosities, proper motions, sizes, metallicities and ages. We use the near-IR VVVX survey database, in combination with Gaia DR2 optical photometry, and with the Two Micron All Sky Survey (2MASS) photometry. We report the detection of a heretofore unknown Galactic Globular Cluster at $RA =$ 14:09:00.0; $DEC=-$65:37:12 (J2000). We calculate a reddening of $E(J-K_s)=(0.3\pm 0.03)$ mag and an extinction of $A_{K_s}=(0.15\pm 0.01)$ mag for this new GC. Its distance modulus and corresponding distance were measured as $(m-M)=(15.93\pm0.03)$ mag and $D=(15.5\pm1.0)$ kpc, respectively. We estimate the metallicity and age by comparison with known GCs and by fitting PARSEC and Dartmouth isochrones, finding $[Fe/H]=(-0.70\pm0.2)$ dex and $t=(11.0\pm1.0)$ Gyr. The mean GC PMs from Gaia are $\mu_{\alpha^\ast}=(-4.68 \pm 0.47 )$ mas $yr^{-1}$ and $\mu_{\delta}=(-1.34 \pm 0.45)$ mas $yr^{-1}$. The total luminosity of our cluster is estimated to be $M_{Ks}=(-7.76\pm 0.5)$ mag. We have found a new low-luminosity, old and metal-rich globular cluster, situated in the far side of the Galactic disk, at $R_{G}=11.2$ kpc from the Galactic centre, and at $z=1.0$ kpc below the plane. Interestingly, the location, metallicity and age of this globular cluster are coincident with the Monoceros Ring (MRi) structure.

CONTEXT. The GOSSS spectral classifications and Gaia data have significantly improved our ability to measure distances and determine memberships of stellar groups with OB stars. AIMS. We have started a program to identify, measure distances, and determine the membership of Galactic stellar groups with OB stars. We start with the identification and distance determinations of groups with O stars. In this paper we concentrate on groups that contain stars with the earliest spectral subtypes. METHODS. We use GOSSS to select Galactic stellar groups with O2-O3.5 stars and a method that combines Gaia DR2 photometry, positions, proper motions, and parallaxes to assign robust memberships and measure distances. We also include the two clusters in that paper to generate our first list of 16 O-type Galactic stellar groups. RESULTS. We derive distances, determine the membership, and analyze the structure of sixteen Galactic stellar groups with O stars, Villafranca O-001 to O-016, including the 14 groups with the earliest-O-type optically-accessible stars known in the Milky Way. We compare our distance with previous results and establish that the best consistency is with VLBI parallaxes and the worst is with kinematic distances. Our results indicate that massive stars can form in relatively low-mass clusters or even in near-isolation, as is the case for the Bajamar star in the North America nebula. This lends support to the hierarchical scenario of star formation, where some stars are born in well-defined bound clusters but others are born in associations that are unbound from the beginning: groups of newborn stars come in many shapes and sizes. We propose that HD 64 568 and HD 64 315 AB could have been ejected simultaneously from Villafranca O-012 S. Our results are consistent with a difference of approx. 20 microas in the Gaia DR2 parallax zero point between bright and faint stars. (ABRIDGED), Comment: Accepted for publication in A&A

We have detected a broad interstellar absorption band centred close to 7700 \r{A} and with a FWHM of 176.6$\pm$3.9 \r{A}. This is the first such absorption band detected in the optical range and is significantly wider than the numerous diffuse interstellar bands (DIBs). It remained undiscovered until now because it is partially hidden behind the A telluric band produced by O$_2$. The band was discovered using STIS@HST spectra and later detected in a large sample of stars of diverse type (OB stars, BA supergiants, red giants) using further STIS and ground-based spectroscopy. The EW of the band is measured and compared with our extinction and K I $\lambda\lambda$7667.021,7701.093 measurements for the same sample. The carrier is ubiquitous in the diffuse and translucent Galactic ISM but is depleted in the environment around OB stars. In particular, it appears to be absent or nearly so in sightlines rich in molecular carbon. This behaviour is similar to that of the $\sigma$-type DIBs, which originate in the low/intermediate-density UV-exposed ISM but are depleted in the high-density UV-shielded molecular clouds. We also present an update on our previous work on the relationship between $E(4405-5495)$ and $R_{5495}$ and incorporate our results into a general model of the ISM., Comment: 17 pages, accepted for publication in MNRAS

We have started an ambitious program to determine if the full diversity of extinction laws is real or if some of it is due to calibration or methodological issues. Here we start by analyzing the information on NIR extinction in a 2MASS stellar sample with good quality photometry and very red colours. We calculate the extinction at 1 $\mu$m, $A_1$, and the power-law exponent, $\alpha$ ($A_\lambda = A_1 \lambda^{-\alpha}$), for the 2MASS stars located in the extinction trajectory in the $H-K$ vs. $J-H$ plane expected for red giants with $A_1 > 5$ mag. We test the validity of the assumption about the nature of those stars, whether a single or multiple values of $\alpha$ are needed, and the spatial variations of the results. Most ($\sim$83%) of those stars can indeed be explained by high-extinction red giants and the rest is composed of extinguished AGB stars (mostly O-rich), blended sources, and smaller numbers of other objects, a contaminant fraction that can be reduced with the help of Gaia} DR2 data. Galactic red giants experience a NIR extinction with $\alpha\sim 2.27$ and an uncertainty of a few hundredths of a magnitude. There is no significant spread in $\alpha$ even though our sample is widely distributed and has a broad range of extinctions. Differences with previous results are ascribed to the treatment of non-linear photometric effects and/or the contaminant correction. Future research should concentrate in finding the correct functional form for the NIR extinction law. In the appendix we detail the treatment of non-linear photometric effects in the 2MASS bands., Comment: 13 pages, accepted for publication in MNRAS

Lack of high-precision long-term continuous photometric data for large samples of stars has prevented the large-scale exploration of pulsational variability in the OB star regime. As a result, the candidates for in-depth asteroseismic modelling remained limited to a few tens of dwarfs. The TESS nominal space mission has surveyed the southern sky, yielding continuous data of at least 27 d for hundreds of OB stars. We aim to couple TESS data in the southern sky with spectroscopy to study the variability over mass and evolution. We focus mainly on the presence of coherent pulsation modes that may or may not be present in the theoretical instability domains and unravel all frequency behaviour in the amplitude spectra of the TESS data. We compose a sample of 98 OB-type stars observed by TESS in Sectors 1-13 and with available high-resolution spectroscopy gathered by the IACOB and OWN surveys. We present the short-cadence 2-min light curves of dozens of OB-type stars, that have one or more spectra in the IACOB or OWN database. Based on these light curves and their Lomb-Scargle periodograms we perform variability classification and frequency analysis, and place the stars in the spectroscopic Hertzsprung-Russell diagram to interpret the variability in an evolutionary context. We deduce diverse origins of the variability found in all of the 98 OB stars in the TESS data. Among these we find several new variable stars, including three hybrid pulsators, three eclipsing binaries, high frequency modes in a Be star, and potential heat-driven pulsations in two Oe stars. We identify stars for which future asteroseismic modelling is possible, provided mode identification is achieved. By comparing the position of the variables to theoretical instability strips we discuss the current shortcomings in non-adiabatic pulsation theory, and the distribution of pulsators in the upper Hertzsprung-Russell diagram., Comment: Accepted for publication in A&A. 18(+47) pages, 15(+100) figures, 3(+2) tables. Language edited and typos removed. Fig. 3, 14, 15 corrected for small plotting error (results unchanged)

The apparent lack of massive O-type stars near the zero-age main sequence (at ages < 2 Myr) is a topic widely discussed. Different explanations for this elusive detection have been proposed, but no firm conclusions have been reached yet. We reassess this empirical result benefiting from the high-quality spectroscopic observations of >400 Galactic O-type stars gathered by the IACOB and OWN surveys. We used temperatures and gravities from a iacob-gbat/fastwind spectroscopic analysis to locate our sample in the Kiel and spectroscopic HR diagrams. We evaluated the completeness of our sample of stars, observational biases using information from the Galactic O star catalog (GOSC), systematics of our methodology, and compare with other recent studies using smaller samples of Galactic O-type stars. We base our discussion on the spectroscopic HR diagram to avoid the use of uncertain distances. We performed a detailed study of the young cluster Trumpler-14 as an example of how Gaia cluster distances can help to construct the associated classical HR diagram. The apparent lack of massive O-type stars near the ZAMS with masses between 30 and 70 Msol persist even when spectroscopic results from a large, non-biased sample of stars are used. We do not find correlation between the dearth of stars and observational biases, limitations of our methodology, or the use of spectroscopic HR diagram instead of the classical one. Investigating the efficiency of mass accretion during the formation process we conclude that an adjustment of the accretion rate towards lower values could reconcile the hotter boundary of detected O-type stars and the theoretical birthline. Last, we discuss that the presence of a small sample of O2-O3.5 stars found closer to the ZAMS might be explained taking into account non-standard star evolution (e.g. binary interaction, mergers, or homogeneous evolution)., Comment: 20 pages, 15 figures, accepted for publication in Astronomy & Astrophysics

We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD~93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites \textit{Chandra} and \textit{NuSTAR} and the $\gamma$-ray satellite \textit{AGILE}. High-energy emission coincident with HD~93129A was detected in the X-ray band up to $\sim$18~keV, whereas in the $\gamma$-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, $f_\mathrm{NT,e} < 0.02$. In addition, we set a lower limit for the magnetic field in the wind-collision region as $B_\mathrm{WCR} > 0.3$~G. We also argue a putative interpretation of the emission from which we estimate $f_\mathrm{NT,e} \approx 0.006$ and $B_\mathrm{WCR} \approx 0.5$~G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries., Comment: 11 pages, 10 figures, accepted for publication in MNRAS

CONTEXT. Many O-type stars have nearby companions whose presence hamper their characterization through spectroscopy. AIMS. We want to obtain spatially resolved spectroscopy of close O-type visual binaries to derive their spectral types. METHODS. We use the Space Telescope Imaging Spectrograph (STIS) of the Hubble Space Telescope (HST) to obtain long-slit blue-violet spectroscopy of eight Galactic O-type stars with nearby visual companions and use spatial-profile fitting to extract the separate spectra. We also use the ground-based Galactic O-Star Spectroscopic Survey (GOSSS) to study more distant visual components. RESULTS. We spatially resolve seven of the eight systems, present spectra for their components, and obtain their spectral types. Those seven multiple systems are iota Ori Aa,Ab,B, 15 Mon Aa,Ab,C, tau CMa Aa,Ab,B,C,D,E, HD 206 267 Aa,Ab,C,D, HD 193 443 A,B, HD 16 429 Aa,Ab, and IU Aur A,B. This is the first time that spatially resolved spectroscopy of the close visual binaries of those systems is obtained. We establish the applicability of the technique as a function of the separation and magnitude difference of the binary., Comment: Accepted for publication in A&A

The process of massive star ($M\geq8~M_\odot$) formation is still poorly understood. Observations of massive young stellar objects (MYSOs) are challenging due to their rarity, short formation timescale, large distances, and high circumstellar extinction. Here, we present the results of a spectroscopic analysis of a population of MYSOs in the Large Magellanic Cloud (LMC). We took advantage of the spectral resolution and wavelength coverage of X-shooter (300-2500 nm), mounted on the European Southern Observatory Very Large Telescope, to detect characteristic spectral features in a dozen MYSO candidates near 30 Doradus, the largest starburst region in the Local Group hosting the most massive stars known. The X-shooter spectra are strongly contaminated by nebular emission. We used a scaling method to subtract the nebular contamination from our objects. We detect H$\alpha,\beta$, [O I] 630.0 nm, Ca II infrared triplet, [Fe II] 1643.5 nm, fluorescent Fe II 1687.8 nm, H$_2$ 2121.8 nm, Br$\gamma$, and CO bandhead emission in the spectra of multiple candidates. This leads to the spectroscopic confirmation of 10 candidates as bona fide MYSOs. We compare our observations with photometric observations from the literature and find all MYSOs to have a strong near-infrared excess. We compute lower limits to the brightness and luminosity of the MYSO candidates, confirming the near-infrared excess and the massive nature of the objects. No clear correlation is seen between the Br$\gamma$ luminosity and metallicity. Combining our sample with other LMC samples results in a combined detection rate of disk features such as fluorescent Fe II and CO bandheads which is consistent with the Galactic rate (40\%). Most of our MYSOs show outflow features., Comment: Accepted to A&A, 26 pages, 19 figures, 4 tables

Binary systems composed of a red supergiant and a B star are useful probes of stellar evolution. We have searched the literature to create a catalog of 108 Galactic systems of such type, which is presented here., Comment: Table with 108 objects included

Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. We aim to derive the atmospheric parameters of the 31 SB2s in the TMBM sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. 32 epochs of FLAMES/GIRAFFE spectra are analysed using spectral disentangling to construct the individual spectra of 62 components. We apply the CMFGEN atmosphere code to determine their stellar parameters and their He, C and N surface abundances. From these properties, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. Components filling their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are rejuvenated. Their locations in the N-vsini diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched and rapidly rotating non-enriched) that cannot be reproduced through single-star population synthesis. This sample is the largest sample of binaries to be studied in such a homogeneous way. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks., Comment: accepted to A&A Core of the paper 15 pages, Appendix: 41 pages

An improved high-resolution and deep A$_{Ks}$ foreground dust extinction map is presented for the Galactic disk area within $295^{\circ} \lesssim l \lesssim 350^{\circ}$, $-1.0^{\circ} \lesssim b \lesssim +1.0^{\circ}$. At some longitudes the map reaches up to $|b|\sim2.25^{\circ}$, for a total of $\sim$148 deg$^2$. The map was constructed via the Rayleigh-Jeans Color Excess (RJCE) technique based on deep near-infrared (NIR) and mid-infrared (MIR) photometry. The new extinction map features a maximum bin size of 1', and relies on NIR observations from the Two Micron All-Sky Survey (2MASS) and new data from ESO's Vista Variables in the V\'ia L\'actea (VVV) survey, in concert with MIR observations from the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE). The VVV photometry penetrates $\sim$4 magnitudes fainter than 2MASS, and provides enhanced sampling of the underlying stellar populations in this heavily obscured region. Consequently, the new results supersede existing RJCE maps tied solely to brighter photometry, revealing a systematic underestimation of extinction in prior work that was based on shallower data. The new high-resolution and large-scale extinction map presented here is readily available to the community through a web query interface., Comment: 11 pages, 9 figures, published in MNRAS

The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus., Comment: Updated to reflect the published version (MNRAS, 489, 241). For a short introductory video of the S-PLUS project, see https://youtu.be/yc5kHrHU9Jk - The S-PLUS Data Release 1 is available at http://datalab.noao.edu/splus

[ABRIDGED] AIMS. MONOS is collecting information and studying O-type spectroscopic binaries with delta > -20 deg. In this 1st paper we describe the sample and provide spectral classifications and additional information for objects with previous SB or EB orbits. In the future we will test the validity of previous solutions and calculate new SB orbits. METHODS.The spectra in this paper have 2 sources: GOSSS, which is obtaining blue-violet R~2500 spectroscopy of thousands of massive stars, and LiLiMaRlin, a library of libraries of high-resolution spectroscopy of massive stars obtained from 4 different surveys and additional data from our own observing programs and public archives. We also use lucky images from AstraLux. RESULTS. We present homogeneous spectral classifications for 92 O-type spectroscopic multiple systems and 10 optical companions. We discuss the visual multiplicity of each system with the support of AstraLux images and additional sources. For 11 O-type objects and for 6 B-type objects we present their first GOSSS spectral classifications. For 2 known EBs we detect SB2/SB1 lines for the first time, to which we add a third system already reported by us recently. For 2 previous SB1 systems we detect their SB2 nature for the first time and give their first separate spectral classifications, something we also do for a 3rd object just recently identified as a SB2. We also detect 9 new astrometric companions and provide updated information on several others. For sigma Ori AaAbB we provide spectral classifications for the 3 components with a single observation for the first time thanks to a lucky spectroscopy observation obtained close to the Aa,Ab periastron and for theta^1 Ori CaCb we add it to the class of Galactic Of?p stars, raising the number of its members to 6. Our sample of O-type spectroscopic binaries contains more triple/higher-order systems than double systems., Comment: 36 pages, accepted for publication in A&A

This paper describes the characterization of the GALANTE photometric system, a seven intermediate- and narrow-band filter system with a wavelength coverage from 3000 $\r{A}$ to 9000 $\r{A}$ . We describe the photometric system presenting the full sensitivity curve as a product of the filter sensitivity, CCD, telescope mirror, and atmospheric transmission curves, as well as some first- and second-order moments of this sensitivity function. The GALANTE photometric system is composed of four filters from the J-PLUS photometric system, a twelve broad-to-narrow filter system, and three exclusive filters, specifically designed to measure the physical parameters of stars such as effective temperature $T_{\rm eff}$, $\log(g)$, metallicity, colour excess $E(4405-5495)$, and extinction type $R_{5495}$. Two libraries, the Next Generation Spectral Library (NGSL) and the one presented in Ma\'iz Apell\'aniz & Weiler (2018), have been used to determine the transformation equations between the Sloan Digital Sky Survey ($\textit{SDSS}$) $\textit{ugriz}$ photometry and the GALANTE photometric system. We will use this transformation to calibrate the zero points of GALANTE images. To this end, a preliminary photometric calibration of GALANTE has been made based on two different $\textit{griz}$ libraries ($\textit{SDSS}$ DR12 and ATLAS All-Sky Stellar Reference Catalog, hereinafter $\textit{RefCat2}$). A comparison between both zero points is performed leading us to the choice of $\textit{RefCat2}$ as the base catalogue for this calibration, and applied to a field in the Cyg OB2 association., Comment: Accepted in MNRAS

We report the first confirmed detection of the galaxy cluster VVV-J144321-611754 at very low latitudes (l = 315.836$^{\circ}$, b = -1.650$^{\circ}$) located in the tile d015 of the VISTA Variables in the V\'ia L\'actea (VVV) survey. We defined the region of 30$\times$ 30 $arcmin^2$ centered in the brightest galaxy finding 25 galaxies. For these objects, extinction-corrected median colors of (H - K$_{s}$) = 0.34 $\pm$ 0.05 mag, (J - H) = 0.57 $\pm$ 0.08 mag and (J - K$_{s}$) = 0.87 $\pm$ 0.06 mag, and R$_{1/2}$ = 1.59 $\pm$ 0.16 $arcsec$; C = 3.01 $\pm$ 0.08; and Sersic index, n = 4.63 $\pm$ 0.39 were estimated. They were visually confirmed showing characteristics of early-type galaxies in the near-IR images. An automatic clustering analysis performed in the whole tile found that the concentration of galaxies VVV-J144321-611754 is a real, compact concentration of early-type galaxies. Assuming a typical galaxy cluster with low X-ray luminosity, the photometric redshift of the brightest galaxy is $z = $ 0.196 $\pm$ 0.025. Follow-up near-IR spectroscopy with FLAMINGOS-2 at the Gemini-South telescope revealed that the two brighter cluster galaxies have typical spectra of early-type galaxies and the estimated redshift for the brightest galaxy VVV-J144321.06-611753.9 is $z =$ 0.234$\pm$0.022 and for VVV-J144319.02-611746.1 is $z =$ 0.232$\pm$0.019. Finally, these galaxies clearly follow the cluster Red Sequence in the rest-frame near-IR color--magnitude diagram with the slope similar to galaxy cluster at redshift of 0.2. These results are consistent with the presence of a bona fide galaxy cluster beyond the Milky Way disk., Comment: Manuscript accepted for publication in The Astrophysical Journal

We present the physical characterization of FSR 1758, a new large, massive object very recently discovered in the Galactic Bulge. The combination of optical data from the 2nd Gaia Data Release (GDR2) and the DECam Plane Survey (DECaPS), and near-IR data from the VISTA Variables in the V\'{\i}a L\'actea Extended Survey (VVVX) led to a clean sample of likely members. Based on this integrated dataset, position, distance, reddening, size, metallicity, absolute magnitude, and proper motion of this object are measured. We estimate the following parameters: $\alpha=17:31:12$, $\delta=-39:48:30$ (J2000), $D=11.5 \pm 1.0$ kpc, $E(J-Ks)=0.20 \pm 0.03$ mag, $R_c=10$ pc, $R_t=150$ pc, $[Fe/H]=-1.5 \pm 0.3$ dex, $M_i < -8.6 \pm 1.0$, $\mu_{\alpha} = -2.85$ mas yr$^{-1}$, and $\mu_{\delta} = 2.55$ mas yr$^{-1}$. The nature of this object is discussed. If FRS 1758 is a genuine globular cluster, it is one of the largest in the Milky Way, with a size comparable or even larger than that of $\omega$ Cen, being also an extreme outlier in the size vs. Galactocentric distance diagram. The presence of a concentration of long-period RR Lyrae variable stars and blue horizontal branch stars suggests that it is a typical metal-poor globular cluster of Oosterhoff type II. Further exploration of a larger surrounding field reveals common proper motion stars, suggesting either tidal debris or that FRS\,1758 is actually the central part of a larger extended structure such as a new dwarf galaxy, tentatively named as Scorpius. In either case, this object is remarkable, and its discovery graphically illustrates the possibility to find other large objects hidden in the Galactic Bulge using future surveys., Comment: Accepted Astrophysical Journal Letters

We report new spectroscopic observations of the magnetic Of?p star HD 148937 obtained since 2015 that differ qualitatively from its extensive historical record of weak, periodic spectral variations. This remarkable behaviour represents clear evidence for an unprecedented change in the character of variability of the star. In this paper we describe the new spectral properties and compare them to the previous line profiles. Based on measurements of the radial velocities of the C iii/N iii} emission lines near 4640 \AA\ and the C iv absorption lines near 5800 \AA, we infer that HD 148937 is likely a high-mass, double-lined spectroscopic binary. Combining the spectroscopic orbit with an archival interferometric measurement of the apparent separation of the equal-brightness components, we tentatively conclude that HD 148937 consists of two O-type stars with masses of approximately 34 and $49~M_\odot$, orbiting in an eccentric ($e=0.75$), long-period ($P_{\rm orb}\sim 26$ y) orbit. We discuss the potential relationship of the binary system to the peculiar properties of HD 148937, and propose future observations to refine the orbital and stellar properties., Comment: 12 pages, accepted by MNRAS

In this paper, we present high-resolution HST/COS observations of the extreme magnetic O star NGC 1624-2. These represent the first ultraviolet spectra of this archetypal object. We examine the variability of its wind-sensitive resonance lines, comparing it to that of other known magnetic O stars. In particular, the observed variations in the profiles of the CIV and SiIV doublets between low state and high state are the largest observed in any magnetic O-type star, consistent with the expected properties of NGC 1624-2's magnetosphere. We also observe a redshifted absorption component in the low state, a feature not seen in most stars. We present preliminary modelling efforts based on the Analytic Dynamical Magnetosphere (ADM) formalism, demonstrating the necessity of using non-spherically symmetric models to determine wind/magnetospheric properties of magnetic O stars., Comment: 13 pages, 6 figures, accepted for publication by MNRAS (23 November 2018)

GALANTE is an optical photometric survey with seven intermediate/narrow filters that has been covering the Galactic Plane since 2016 using the Javalambre T80 and Cerro Tololo T80S telescopes. The P.I.s of the northern part (GALANTE NORTE) are Emilio J. Alfaro and Jes\'us Ma\'iz Apell\'aniz. and the P.I. of the southern part (GALANTE SUR) is Rodolfo H. Barb\'a. The detector has a continuous 1.4 degr x 1.4 degr field of view with a sampling of 0.55"/pixel and the seven filters are optimized to detect obscured early-type stars. The survey includes long, intermediate, short, and ultrashort exposure times to reach a dynamical range close to 20 magnitudes, something never achieved for such an optical project before. The characteristics of GALANTE allow for a new type of calibration scheme using external Gaia, Tycho-2, and 2MASS data that has already led to a reanalysis of the sensitivity of the Gaia G filter. We describe the project and present some early results. GALANTE will identify the majority of the early-type massive stars within several kpc of the Sun and measure their amount and type of extinction. It will also map the Halpha nebular emission, identify emission-line stars, and do other studies of low- and intermediate-mass stars., Comment: To appear in Highlights on Spanish Astrophysics X, Proceedings of the XIII Scientific Meeting of the Spanish Astronomical Society held on July 16-20, 2018, in Salamanca, Spain

We present results on two different Gaia-related topics. First, we describe our efforts to calibrate the three Gaia photometric passbands G, G_BP, and G_RP. We have built a new spectrophotometric HST/STIS library and used it to derive new sensitivity curves and zero points for the three bands, including recipes on how to correct some cases. Second, we present our results on Galactic runaway stars using Gaia DR1 proper motions: we detect 76 runaway stars, 17 (possibly 19) of them not previously identified as such., Comment: To appear in Highlights on Spanish Astrophysics X, Proceedings of the XIII Scientific Meeting of the Spanish Astronomical Society held on July 16-20, 2018, in Salamanca, Spain

LiLiMaRlin is a library of libraries of massive-star high-resolution optical spectra built by collecting data from [a] our spectroscopic surveys (OWN, IACOB. NoMaDS, and CAF\'E-BEANS) and programs and [b] searches in public archives. The current version has 18 077 spectra of 1665 stars obtained with seven different telescopes (HET 9.2 m, NOT 2.56 m, CAHA 2.2 m, MPG/ESO 2.2 m, OHP 1.93 m, Mercator 1.2 m, and Stella 1.2 m). All the spectra have been filtered to eliminate misidentifications and bad-quality ones, uniformly reprocessed, and placed on a common format. We present applications of this library of libraries to the analysis of spectroscopic binaries (OWN and MONOS, see poster by E. Trigueros P\'aez at this meeting) and the study of the interstellar medium (CollDIBs). We discuss our plans for the future., Comment: To appear in Highlights on Spanish Astrophysics X, Proceedings of the XIII Scientific Meeting of the Spanish Astronomical Society held on July 16-20, 2018, in Salamanca, Spain

Context. Most massive stars are in binary or multiple systems. Several massive stars have been detected as doublelined spectroscopic binaries and among these, the OWN Survey has detected a non-negligible number whose components show very different spectral line broadening (i.e., projected rotational velocities). This fact raises a discussion about the contributing processes, such as angular-momentum transfer and tidal forces. Aims. We seek to constrain the physical and evolutionary status of one of such systems, the O+O binary HD 93343. Methods. We analyzed a series of high-resolution multiepoch optical spectra to determine the orbital parameters, projected rotational velocities, and evolutionary status of the system. Results. HD 93343 is a binary system comprised of two O7.5 Vz stars that each have minimum masses of approximately 22 Mo in a wide and eccentric orbit (e = 0.398$\pm$0.004; P=50.432$\pm$0.001 d). Both stars have very similar stellar parameters, and hence ages. As expected from the qualitative appearance of the combined spectrum of the system, however, these stars have very different projected rotational velocities (~65 and ~325 km/s, respectively). Conclusions. The orbits and stellar parameters obtained for both components seem to indicate that their youth and relative separation is enough to discard the effects of mass transfer and tidal friction. Thus, non-synchronization should be intrinsic to their formation., Comment: 9 pages, 7 figures. Accepted to Astronomy & Astrophysics

CONTEXT.The first Gaia Data Release (DR1) significantly improved the previously available proper motions for the majority of the Tycho-2 stars. AIMS. We want to detect runaway stars using Gaia DR1 proper motions and compare our results with previous searches. METHODS. Runaway O stars and BA supergiants are detected using a 2-D proper-motion method. The sample is selected using Simbad, spectra from our GOSSS project, literature spectral types, and photometry processed using CHORIZOS. RESULTS. We detect 76 runaway stars, 17 (possibly 19) of them with no prior identification as such, with an estimated detection rate of approximately one half of the real runaway fraction. An age effect appears to be present, with objects of spectral subtype B1 and later having travelled for longer distances than runaways of earlier subtypes. We also tentatively propose that the fraction of runaways is lower among BA supergiants that among O stars but further studies using future Gaia data releases are needed to confirm this. The frequency of fast rotators is high among runaway O stars, which indicates that a significant fraction of them (and possibly a majority) is produced in supernova explosions., Comment: Accepted for publication in A&A

CONTEXT: Many massive stars have nearby companions whose presence hamper their characterization through spectroscopy. AIMS: We want to obtain spatially resolved spectroscopy of close massive visual binaries to derive their spectral types. METHODS: We obtain a large number of short long-slit spectroscopic exposures of five close binaries under good seeing conditions, select those with the best characteristics, extract the spectra using multiple-profile fitting, and combine the results to derive spatially separated spectra. RESULTS: We demonstrate the usefulness of Lucky Spectroscopy by presenting the spatially resolved spectra of the components of each system, in two cases with separations of only ~0.3". Those are delta Ori Aa+Ab (resolved in the optical for the first time) and sigma Ori AaAb+B (first time ever resolved). We also spatially resolve 15 Mon AaAb+B, zeta Ori AaAb+B (both previously resolved with GOSSS, the Galactic O-Star Spectroscopic Survey), and eta Ori AaAb+B, a system with two spectroscopic B+B binaries and a fifth visual component. The systems have in common that they are composed of an inner pair of slow rotators orbited by one or more fast rotators, a characteristic that could have consequences for the theories of massive star formation., Comment: Accepted for publication in A&A, 7 pages

[ABRIDGED] Context. O stars are excellent tracers of the intervening ISM because of their high luminosity, blue intrinsic SED, and relatively featureless spectra. We are currently conducting GOSSS, which is generating a large sample of O stars with accurate spectral types within several kpc of the Sun. Aims. To obtain a global picture of the properties of dust extinction in the solar neighborhood based on optical-NIR photometry of O stars with accurate spectral types. Methods. We have processed a photometric set with the CHORIZOS code to measure the amount and type of extinction towards 562 O-type stellar systems. We have tested three different families of extinction laws and analyzed our results with the help of additional archival data. Results. The Ma\'iz Apell\'aniz et al. (2014) family of extinction laws provides a better description of Galactic dust that either the Cardelli et al. (1989) or Fitzpatrick (1999) families, so it should be preferentially used. In many cases O stars and late-type stars experience similar amounts of extinction at similar distances but some O stars are located close to the molecular clouds left over from their births and have larger extinctions than the average for nearby late-type populations. In qualitative terms, O stars experience a more diverse extinction than late-type stars, as some are affected by the small-grain-size, low-R_5495 effect of molecular clouds and others by the large-grain-size, high-R_5495 effect of H II regions. Late-type stars experience a narrower range of grain sizes or R_5495, as their extinction is predominantly caused by the average, diffuse ISM. We propose that the reason for the existence of large-grain-size, high-R_5495 regions in the ISM in the form of H II regions and hot-gas bubbles is the selective destruction of small dust grains by EUV photons and possibly by thermal sputtering by atoms or ions., Comment: 33 pages, accepted for publication in A&A. This second version includes language editing, other minor changes to the text, and the solution to formatting issues in the main table

The IACOB and OWN surveys are two ambitious complementary observational projects which have made available a large multi-epoch spectroscopic database of optical high resolution spectra of Galactic massive O-type stars. As a first step in the study of the full sample of (more than 350) O stars surveyed by the IACOB/OWN projects, we have performed the quantitative spectroscopic analysis of a subsample of 128 stars included in the modern grid of O-type standards for spectral classification. We use semi-automatized tools to determine the set of spectroscopic parameters that can be obtained from the optical spectrum of O-type stars. We also benefit from the multi-epoch character of the surveys to perform a spectroscopic variability study of the sample, accounting for spectroscopic binarity and variability of the main wind diagnostic lines. We provide a general overview of the stellar and wind parameters of this reference sample, and updated recipes for the SpT\,--\,Teff/log g calibrations for Galactic O-type stars. We evaluate our semi-automatized analysis strategy with $\sim$40 stars from the literature, and find a good agreement. The agreement between the synthetic spectra associated with fastwind best fitting models and the observed spectra is good for most targets, but 46 stars present a particular behavior of the wind diagnostic lines that cannot be reproduced by our grid of spherically symmetric unclumped models. These are potential targets of interest for more detailed investigations of clumpy winds and/or the existence of additional circumstellar components. Last, our variability study has led to the detection of signatures of spectroscopic binarity in 27\% of the stars and small amplitude radial velocity variations in the photospheric lines of another 30\%. Additionally, 31\% of the investigated stars show variability in the wind diagnostic lines., Comment: 20 pages, 18 figures, accepted for publication in Astronomy & Astrophysics

Characterization of the dynamics of massive star systems and the astrophysical properties of the interacting components are a prerequisite for understanding their formation and evolution. Optical interferometry at milliarcsecond resolution is a key observing technique for resolving high-mass multiple compact systems. Here we report on VLTI/GRAVITY, Magellan/FIRE, and MPG2.2m/FEROS observations of the late-O/early-B type system HD 93206 A, which is a member of the massive cluster Collinder 228 in the Carina nebula complex. With a total mass of about 90 M_sun, it is one of the most compact massive-quadruple systems known. In addition to measuring the separation and position angle of the outer binary Aa - Ac, we observe BrG and HeI variability in phase with the orbital motion of the two inner binaries. From the differential phases analysis, we conclude that the BrG emission arises from the interaction regions within the components of the individual binaries, which is consistent with previous models for the X-ray emission of the system based on wind-wind interaction. With an average 3-sigma deviation of ~15 deg, we establish an upper limit of p ~ 0.157 mas (0.35 AU) for the size of the BrG line-emitting region. Future interferometric observations with GRAVITY using the 8m UTs will allow us to constrain the line-emitting regions down to angular sizes of 20 uas (0.05 AU at the distance of the Carina nebula)., Comment: 17pages, 12 figures Accepted for publication in ApJ

Massive O-type stars play a dominant role in our Universe, but many of their properties remain poorly constrained. In the last decade magnetic fields have been detected in all Galactic members of the distinctive Of?p class, opening the door to a better knowledge of all O-type stars. With the aim of extending the study of magnetic massive stars to nearby galaxies, to better understand the role of metallicity in the formation of their magnetic fields and magnetospheres, and to broaden our knowledge of the role of magnetic fields in massive star evolution, we have carried out spectropolarimetry of five extra-Galactic Of?p stars, as well as a couple of dozen neighbouring stars. We have been able to measure magnetic fields with typical error bars from 0.2 to 1.0 kG, depending on the apparent magnitude and on weather conditions. No magnetic field has been firmly detected in any of our measurements, but we have been able to estimate upper limits to the field values of our target stars. One of our targets, 2dFS 936, exhibited an unexpected strengthening of emission lines. We confirm the unusual behaviour of BI 57, which exhibits a 787 d period with two photometric peaks and one spectroscopic maximum. The observed strengthening of the emission lines of 2dFS 936, and the lack of detection of a strong magnetic field in a star with such strong emission lines is at odd with expectations. Together with the unusual periodic behaviour of BI 57, it represents a challenge for the current models of Of?p stars. The limited precision that we obtained in our field measurements (in most cases as a consequence of poor weather) has led to field-strength upper limits that are substantially larger than those typically measured in Galactic magnetic O stars. Further higher precision observations and monitoring are clearly required., Comment: Accepted by A&A