Search

Massive stars lead the ionization and mechanical feedback within young star-forming regions. The Large Magellanic Cloud (LMC) is an ideal galaxy for studying individual massive stars and quantifying their feedback contribution to the environment. We analyze eight exemplary targets in LMC N11 B from the Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) program, using novel spectra from HST (COS and STIS) in the UV, and from VLT (X-shooter) in the optical. We model the spectra of early to late O-type stars by using state-of-the-art PoWR atmosphere models. We determine the stellar and wind parameters (e.g., $T_\star$, $\log g$, $L_{\star}$, $\dot{M}$, $v_\infty$) of the analyzed objects, chemical abundances (C, N, O), ionizing and mechanical feedback ($Q_\mathrm{H}$, $Q_\mathrm{He{\small{I}}}$, $Q_\mathrm{He{\small{II}}}$, $L_\mathrm{mec}$) and X-rays. We report ages of $2-4.5$ Myr and masses of $30-60$ $M_\odot$ for the analyzed stars in N11 B, consistent with a scenario of sequential star formation. We note that the observed wind-momentum luminosity relation is consistent with theoretical predictions. We detect nitrogen enrichment in most of the stars, up to a factor of seven. However, we do not find a correlation between nitrogen enrichment and projected rotational velocity. Finally, based on their spectral type, we estimate the total ionizing photons injected from the O-type stars in N11 B into its environment. We report $\log$ ($\sum$ $Q_\mathrm{H}$)$=50.5$ ph s$^{-1}$, $\log$ ($\sum$ $Q_\mathrm{He{\small{I}}}$)$=49.6$ ph s$^{-1}$ and $\log$ ($\sum$ $Q_\mathrm{He{\small{II}}}$)$=44.4$ ph s$^{-1}$, consistent with the total ionizing budget in N11., Comment: Accepted in Astronomy and Astrophysics

We present a detailed study of the large-scale shock front in Stephan's Quintet, a byproduct of past and ongoing interactions. Using integral-field spectroscopy from the new William Herschel Telescope Enhanced Area Velocity Explorer (WEAVE), recent 144 MHz observations from the LOFAR Two-metre Sky Survey (LoTSS), and archival data from the Very Large Array and James Webb Space Telescope (JWST), we obtain new measurements of key shock properties and determine its impact on the system. Harnessing the WEAVE large integral field unit's (LIFU) field of view (90 $\times$ 78 arcsec$^{2}$), spectral resolution ($R\sim2500$) and continuous wavelength coverage across the optical band, we perform robust emission line modeling and dynamically locate the shock within the multi-phase intergalactic medium (IGM) with higher precision than previously possible. The shocking of the cold gas phase is hypersonic, and comparisons with shock models show that it can readily account for the observed emission line ratios. In contrast, we demonstrate that the shock is relatively weak in the hot plasma visible in X-rays (with Mach number of $\mathcal{M} \sim 2 - 4$), making it inefficient at producing the relativistic particles needed to explain the observed synchrotron emission. Instead, we propose that it has led to an adiabatic compression of the medium, which has increased the radio luminosity ten-fold. Comparison of the Balmer line-derived extinction map with the molecular gas and hot dust observed with JWST suggests that pre-existing dust may have survived the collision, allowing the condensation of H$_{2}$ - a key channel for dissipating the shock energy., Comment: 23 pages, 15 figures, accepted for publication in MNRAS

Surveys in the Milky Way and Large Magellanic Cloud revealed that the majority of massive stars will interact with companions during their lives. However, knowledge of the binary properties of massive stars at low metallicity, which approaches the conditions of the Early Universe, remains sparse. We present the Binarity at LOw Metallicity (BLOeM) campaign - an ESO large programme designed to obtain 25 epochs of spectroscopy for 929 massive stars in the SMC - the lowest metallicity conditions in which multiplicity is probed to date (Z = 0.2 Zsun). BLOeM will provide (i) the binary fraction, (ii) the orbital configurations of systems with periods P < 3 yr, (iii) dormant OB+BH binaries, and (iv) a legacy database of physical parameters of massive stars at low metallicity. The stars are observed with the LR02 setup of the giraffe instrument of the Very Large Telescope (3960-4570A, resolving power R=6200; typical signal-to-noise ratio S/N=70-100). This paper utilises the first 9 epochs obtained over a three-month time. We describe the survey and data reduction, perform a spectral classification of the stacked spectra, and construct a Hertzsprung-Russell diagram of the sample via spectral-type and photometric calibrations. The sample covers spectral types from O4 to F5, spanning the effective temperature and luminosity ranges 6.5

CONTEXT: The spectral analysis of hot, massive stars is a fundamental astrophysical method to obtain their intrinsic properties and their feedback. Quantitative spectroscopy for hot, massive stars requires detailed numerical modeling of the atmosphere and an iterative treatment to obtain the best solution within a given framework. AIMS: We present an overview of different techniques for the quantitative spectroscopy of hot stars employed within the X-Shooting ULLYSES collaboration, from grid-based approaches to tailored fits. By performing a blind test, we gain an overview about the similarities and differences of the resulting parameters. Our study aims to provide an overview of the parameter spread caused by different approaches. METHODS: For three different stars from the sample (SMC O5 star AzV 377, LMC O7 star Sk -69 50, and LMC O9 star Sk -66 171), we employ different atmosphere codes (CMFGEN, Fastwind, PoWR) and strategies to determine their best-fitting model. For our analyses, UV and optical spectra are used to derive the properties with some methods relying purely on optical data for comparison. To determine the overall spectral energy distribution, we further employ additional photometry from the literature. RESULTS: Effective temperatures for each of three sample stars agree within 3 kK while the differences in log g can be up to 0.2 dex. Luminosity differences of up to 0.1 dex result from different reddening assumptions, which seem to be larger for the methods employing a genetic algorithm. All sample stars are nitrogen-enriched. CONCLUSIONS: We find a reasonable agreement between the different methods. Tailored fitting tends to be able to minimize discrepancies obtained with more course or automatized treatments. UV spectral data is essential for the determination of realistic wind parameters. For one target (Sk -69 50), we find clear indications of an evolved status., Comment: 19+15 pages, 21+4 figures, accepted version (A&A 689, A30) including language editing, condensed abstract

Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z. To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver., Comment: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

The Gaia-ESO Survey is a public spectroscopic survey that has targeted $\gtrsim10^5$ stars covering all major components of the Milky Way from the end of 2011 to 2018, delivering its public final release in May 2022. Unlike other spectroscopic surveys, Gaia-ESO is the only survey that observed stars across all spectral types with dedicated, specialised analyses: from O ($T_\mathrm{eff} \sim 30,000-52,000$~K) all the way to K-M ($\gtrsim$3,500~K). The physics throughout these stellar regimes varies significantly, which has previously prohibited any detailed comparisons between stars of significantly different type. In the final data release (internal data release 6) of the Gaia-ESO Survey, we provide the final database containing a large number of products such as radial velocities, stellar parameters and elemental abundances, rotational velocity, and also, e.g., activity and accretion indicators in young stars and membership probability in star clusters for more than 114,000 stars. The spectral analysis is coordinated by a number of Working Groups (WGs) within the Survey, which specialise in the various stellar samples. Common targets are analysed across WGs to allow for comparisons (and calibrations) amongst instrumental setups and spectral types. Here we describe the procedures employed to ensure all Survey results are placed on a common scale to arrive at a single set of recommended results for all Survey collaborators to use. We also present some general quality and consistency checks performed over all Survey results., Comment: A&A accepted, minor revision, 36 pages, 38 figures

The Gaia-ESO survey sample of massive OB stars in the Carina Nebula consists of 234 stars. The addition of brighter sources from the Galactic O-Star Spectroscopic Survey and additional sources from the literature allows us to create the most complete census of massive OB stars done so far in the region. It contains a total of 316 stars, being 18 of them in the background and four in the foreground. Of the 294 stellar systems in Car OB1, 74 are of O type, 214 are of non-supergiant B type and 6 are of WR or non-O supergiant (II to Ia) spectral class. We identify 20 spectroscopic binary systems with an O-star primary, of which 6 are reported for the first time, and another 18 with a B-star primary, of which 13 are new detections. The average observed double-lined binary fraction of O-type stars in the surveyed region is 0.35, which represents a lower limit. We find a good correlation between the spectroscopic n-qualifier and the projected rotational velocity of the stars. The fraction of candidate runaways among the stars with and without the n-qualifier is 4.4% and 2.4%, respectively, although non resolved double-lined binaries can be contaminating the fast rotators sample., Comment: Accepted for publication in A&A, 62 pages, 8 figures

MEGARA is the optical integral field and multi-object spectrograph at the Gran Telescopio Canarias. We have created MEGASTAR, an empirical library of stellar spectra obtained using MEGARA at high resolution $R=20\,000$ (FWHM), available in two wavelength ranges: one centered in H${\alpha}$, from 6420 to 6790\,\AA\ and the other centered in the \ion{Ca}{ii} triplet, from 8370 to 8885\,\AA\ (\mbox{HR-R} and \mbox{HR-I} VPH-grating configurations). In this work, we use MEGASTAR spectra, combination of these two short wavelength intervals, to estimate the stellar parameters namely effective temperature, surface gravity and metallicity (and their associated errors) for a sample of 351 MEGASTAR members with spectral types earlier than B2. We have applied a $\chi^2$ technique by comparing MEGASTAR data to theoretical stellar models. For those stars with stellar parameters derived in the literature, we have obtained a good agreement between those published parameters and ours. Besides the stellar parameters, we also provide several products like the rectified spectra, radial velocities and stellar indices for this sample of stars. In a near future, we will use MEGASTAR spectra and their derived stellar parameters to compute stellar population evolutionary synthesis models, which will contribute to a better interpretation of star clusters and galaxies spectra obtained with MEGARA., Comment: 20 pages, 14 figures, Monthly Notices of the Royal Astronomical Society (in press)

The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper (arXiv:2206.02901) introduces the survey results. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. The Gaia-ESO Survey obtained 202,000 spectra of 115,000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022., Comment: 38 pages. A&A in press

In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars using FLAMES on the ESO VLT (both Giraffe and UVES spectrographs), covering all the Milky Way populations, with a special focus on open star clusters. This article provides an overview of the survey implementation (observations, data quality, analysis and its success, data products, and releases), of the open cluster survey, of the science results and potential, and of the survey legacy. A companion article (Gilmore et al.) reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. The GES has determined homogeneous good-quality radial velocities and stellar parameters for a large fraction of its more than 110,000 unique target stars. Elemental abundances were derived for up to 31 elements for targets observed with UVES. Lithium abundances are delivered for about 1/3 of the sample. The analysis and homogenisation strategies have proven to be successful; several science topics have been addressed by the Gaia-ESO consortium and the community, with many highlight results achieved. The final catalogue has been released through the ESO archive at the end of May 2022, including the complete set of advanced data products. In addition to these results, the Gaia-ESO Survey will leave a very important legacy, for several aspects and for many years to come., Comment: Accepted for publication in Astronomy and Astrophysics. 30 pages, 30 figures, 4 tables

2MASS J20395358+4222505 is an obscured early B supergiant near the massive OB star association Cyg OB2. Despite its bright infrared magnitude (K$_{s}$=5.82) it has remained largely ignored because of its dim optical magnitude (B=16.63, V=13.68). In a previous paper we classified it as a highly reddened, potentially extremely luminous, early B-type supergiant. We obtained its spectrum in the U, B and R spectral bands during commissioning observations with the instrument MEGARA@GTC. It displays a particularly strong H$\alpha$ emission for its spectral type, B1 Ia. The star seems to be in an intermediate phase between super- and hypergiant, a group that it will probably join in the near (astronomical) future. We observe a radial velocity difference between individual observations and determine the stellar parameters, obtaining T$_{eff}$ = 24000 K, logg$_{c}$= 2.88 $\pm$ 0.15. The rotational velocity found is large for a B-supergiant, vsini= 110 $\pm$ 25 km s$^{-1}$. The abundance pattern is consistent with solar, with a mild C underabundance (based on a single line). Assuming that J20395358+4222505 is at the distance of Cyg OB2 we derive the radius from infrared photometry, finding R= 41.2 $\pm$ 4.0 R$_{\odot}$, log(L/L$_{\odot}$)= 5.71 $\pm$ 0.04 and a spectroscopic mass of 46.5 $\pm$ 15.0 M$_{\odot}$. The clumped mass-loss rate (clumping factor 10) is very high for the spectral type, $\dot{M}$ = 2.4x10$^{-6}$ M$_{\odot}$ a$^{-1}$. The high rotational velocity and mass-loss rate place the star at the hot side of the bi-stability jump. Together with the nearly solar CNO abundance pattern, they may also point to evolution in a binary system, J20395358+4222505 being the initial secondary., Comment: Accepted for publication in MNRAS, 13 pages, 9 figures

The Gaia-ESO Survey (GES) is a large public spectroscopic survey that has collected, over a period of 6 years, spectra of ~ 10^5 stars. This survey provides not only the reduced spectra, but also the stellar parameters and abundances resulting from the analysis of the spectra. The GES dataflow is organised in 19 working groups. Working group 13 (WG13) is responsible for the spectral analysis of the hottest stars (O, B and A type, with a formal cut-off of Teff > 7000 K) that were observed as part of GES. We present the procedures and techniques that have been applied to the reduced spectra, in order to determine the stellar parameters and abundances of these stars. The procedure used is similar to that of other working groups in GES. A number of groups (called `Nodes') each independently analyse the spectra, using their state-of-the-art techniques and codes. Specific for the analysis in WG13 is the large temperature range that is covered (Teff = 7000 - 50,000 K), requiring the use of different analysis codes. Most Nodes can therefore only handle part of the data. Quality checks are applied to the results of these Nodes by comparing them to benchmark stars, and by comparing them one to another. For each star the Node values are then homogenised into a single result: the recommended parameters and abundances. Eight Nodes each analysed (part of) the data. In total 17,693 spectra of 6462 stars were analysed, most of them in 37 open star clusters. The homogenisation led to stellar parameters for 5584 stars. Abundances were determined for a more limited number of stars. Elements studied are He, C, N, O, Ne, Mg, Al, Si and Sc. Abundances for at least one of those elements were determined for 292 stars. The hot-star data analysed here, as well as the Gaia-ESO Survey data in general, will be of considerable use in future studies of stellar evolution and open clusters., Comment: 26 pages, 14 figures, accepted for publication in Astronomy & Astrophysics; language-edited version; two appendices merged

MEGARA is an optical integral field and multiobject fibre based spectrograph for the 10.4m Gran Telescopio CANARIAS that offers medium to high spectral resolutions (FWHM) of R $\simeq$ 6000, 12000, 20000. Commissioned at the telescope in 2017, it started operation as a common-user instrument in 2018. We are creating an instrument-oriented empirical spectral library from MEGARA-GTC stars observations, MEGASTAR, crucial for the correct interpretation of MEGARA data. This piece of work describes the content of the first release of MEGASTAR, formed by the spectra of 414 stars observed with R $\simeq$ 20000 in the spectral intervals from 6420 to 6790 $\r{A}$ and from 8370 to 8885 $\r{A}$, and obtained with a continuum average signal to noise ratio around 260. We describe the release sample, the observations, the data reduction procedure and the MEGASTAR database. Additionally, we include in Appendix A, an atlas with the complete set of 838 spectra of this first release of the MEGASTAR catalogue., Comment: 15 pages, 12 figures and 5 tables. Accepted for publication in MNRAS

Cygnus OB2 provides a unique insight into the high-mass stellar content in one of the largest groups of young massive stars in our Galaxy. Although several studies of its massive population have been carried out over the last decades, an extensive spectroscopic study of the whole known O-star population in the association is still lacking. In this work, we created the most complete spectroscopic census of O stars carried out so far in Cygnus OB2 using already existing and new spectroscopy. We present the spectra for 78 O-type stars, from which we identify new binary systems, obtain the distribution of rotational velocities, and determine the main stellar parameters for all the stars in the region that have not been detected as double-line spectroscopic binaries. We also derive radii, luminosities, and masses for those stars with reliable Gaia astrometry, in addition to creating the Hertzsprung-Russell Diagram to interpret the evolutionary status of the association. This work has shown the improvement reached when using accurate spectroscopic parameters and astrometry for the interpretation of the evolutionary status of a population, revealing, in the case of Cygnus OB2, at least two star-forming bursts at $\sim$3 and $\sim$5 Myr. We find an apparent deficit of very fast rotators in the distribution of rotational velocities. The inspection of the dynamical distribution of the sample has allowed us to identify nine O stars with peculiar proper motions and discuss a possible dynamical ejection scenario or past supernova explosions in the region., Comment: Accepted for publication in A&A, 25 pages, 11 figures

MEGARA (Multi Espectr{\'o}grafo en GTC de Alta Resoluci{\'o}n para Astronom{\'\i}a) is an optical (3650~--~9750\AA), fibre-fed, medium-high spectral resolution (R = 6000, 12000, 20000) instrument for the GTC 10.4m telescope, commissioned in the summer of 2017, and currently in operation. The scientific exploitation of MEGARA demands a stellar-spectra library to interpret galaxy data and to estimate the contribution of the stellar populations. This paper introduces the MEGARA-GTC spectral library, detailing the rationale behind the catalogue building. We present the spectra of 97 stars (21 individual stars and 56 members of the globular cluster M15, being both sub-samples taken during the commissioning runs; and 20 stars from our on-going GTC Open-Time program). The spectra have R~=~20000 in the HR-R and HR-I setups, centred at 6563 and 8633~\AA\ respectively. We describe the procedures to reduce and analyse the data. Then, we determine the best-fitting theoretical models to each spectrum through a $\chi^{2}$ minimisation technique to derive the stellar physical parameters and discuss the results. We have also measured some absorption lines and indices. Finally, this article introduces our project to complete the library and the database to make the spectra available to the community., Comment: 27 paginas, 18 figuras, MNRAS, in press

For the first time, we have explored the spatial substructure of the Cygnus OB2 association using parallaxes from the recent second Gaia data release. We find significant line-of-sight substructure within the association, which we quantify using a parameterised model that reproduces the observed parallax distribution. This inference approach is necessary due to the non-linearity of the parallax-distance transformation and the asymmetry of the resulting probability distribution. By using a Markov Chain Monte Carlo ensemble sampler and an unbinned maximum likelihood test we identify two different stellar groups superposed on the association. We find the main Cygnus OB2 group at ~1760 pc, further away than recent estimates have envisaged, and a foreground group at ~1350 pc. We also calculate individual membership probabilities and identify outliers as possible non-members of the association., Comment: 6 pages, 3 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

Cygnus OB2 is a rich OB association in the Galaxy which has experienced intense star formation in the last 20-25 Myr. Its stellar population shows a correlation between age and Galactic longitude. Exploring the chemical composition of its stellar content we will be able to check the degree of homogeneity of the natal molecular cloud and possible effects of self-enrichment processes. Our aim is to determine silicon and oxygen abundances for a sample of eight early-type slow rotators in Cygnus OB2 in order to check possible inhomogeneities across the whole association and whether there exists a correlation of chemical composition with Galactic longitude. We have performed a spectroscopic analysis of a sample of late O and early B stars with low rotational velocity, which have been chosen so as to cover the whole association area. We have carried out an analysis based on equivalent widths of metal lines, the wings of the H Balmer lines and FASTWIND stellar atmosphere models to determine their stellar fundamental parameters as well as the silicon and oxygen surface abundances. We derive a rather homogeneous distribution of silicon and oxygen abundances across the region, with average values of 12+log(Si/H)=7.53$\pm$0.08 dex and 12+log(O/H)=8.65$\pm$0.12 dex. We find a homogeneous chemical composition in Cygnus OB2 with no clear evidence for significant chemical self-enrichment, despite indications of strong stellar winds and possible supernovae during the history of the region. Comparison with different scenarios of chemical enrichment by stellar winds and supernovae point to star forming efficiencies not significantly above 10%. The degree of homogeneity that we find is consistent with the observed Milky Way oxygen gradient based on HII regions. We also find that the oxygen scatter within Cygnus OB2 is at least of the same order than among HII regions at similar Galactocentric distance., Comment: 15 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

The Cygnus complex is one of the most powerful star forming regions at a close distance from the Sun (~1.4 kpc). Its richest OB association Cygnus OB2 is known to harbor many tens of O-type stars and hundreds of B-type stars, providing a large homogeneous population of OB stars that can be analyzed. Many studies of its massive population have been developed in the last decades, although the total number of OB stars is still incomplete. Our aim is to increase the sample of O and B members of Cygnus OB2 and its surroundings by spectroscopically classifying 61 candidates as possible OB-type members of Cygnus OB2. We have obtained new blue intermediate-resolution spectra suitable for spectral classification of the 61 candidates in Cygnus OB2 and surroundings. We thus performed a spectral classification of the sample using He I-II and metal lines rates, as well as the Marxist Ghost Buster (MGB) software for O-type stars and the IACOB standards catalog for B-type stars. Out of the 61 candidates, we have classified 42 stars as new massive OB-type stars, earlier than B3, in Cygnus OB2 and surroundings, including 11 O-type stars. The other candidates are discarded as they display later spectral types inconsistent with membership in the association. However, the magnitude cutoff and dust extinction introduce an incompleteness. Many O and early B stars at B > 16 mag are still undiscovered in the region. Finally, we have studied the age and extinction distribution of our sample within the region, placing them in the Hertzsprung-Russell Diagram using different stellar models in order to assess age uncertainties. Massive star formation in Cygnus OB2 seems to have proceeded from lower to higher Galactic longitudes, regardless of the details of the models used. The correlation between age and Galactic longitude previously found in the region is now confirmed., Comment: 20 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

The Gaia-ESO Survey (GES) is obtaining high quality spectra of $\sim$10$^{5}$ stars in our Galaxy, providing an homogeneous and unique overview of all the main components of the Milky Way, its formation history and the evolution of young, mature and ancient Galactic populations. Our group is in charge of the early-type massive stars that define the youngest population in the survey. In this contribution, we present the results of the quantitative spectroscopic analysis of O-type stars in the Carina Nebula within the Gaia-ESO Survey. For this aim, we have used FASTWIND and CMFGEN stellar atmosphere codes, providing stellar parameters for the current sample (GES data release iDR4)., Comment: 4 pages, 3 figures, to be published in Highlights on Spanish Astrophysics IX, Proceedings of the XII Scientific Meeting of the Spanish Astronomical Society held on July 18-22, 2016, in Bilbao, Spain. F. Figueras, A. S\'anchez-Lavega, S. P\'erez-Hoyos, A. Alonso, S. Arribas, C. Hern\'andez Monteagudo (eds.)