Refine
Document Type
- Article (41)
- Doctoral Thesis (1)
Language
- English (42) (remove)
Is part of the Bibliography
- yes (42)
Keywords
- stars: massive (23)
- stars: Wolf-Rayet (16)
- stars: mass-loss (16)
- stars: winds, outflows (14)
- stars: atmospheres (13)
- stars: early-type (10)
- Magellanic Clouds (8)
- binaries: close (8)
- stars: evolution (8)
- binaries: spectroscopic (6)
Institute
- Institut für Physik und Astronomie (42) (remove)
We report the first-ever discovery of a Wolf-Rayet (WR) star in the Large Magellanic Cloud via detection of a circular shell with the Spitzer Space Telescope. Follow-up observations with Gemini-South resolved the central star of the shell into two components separated from each other by a parts per thousand 2 arcsec (or a parts per thousand 0.5 pc in projection). One of these components turns out to be a WN3 star with H and He lines both in emission and absorption (we named it BAT99 3a using the numbering system based on extending the Breysacher et al. catalogue). Spectroscopy of the second component showed that it is a B0 V star. Subsequent spectroscopic observations of BAT99 3a with the du Pont 2.5-m telescope and the Southern African Large Telescope revealed that it is a close, eccentric binary system, and that the absorption lines are associated with an O companion star. We analysed the spectrum of the binary system using the non-LTE Potsdam WR (powr) code, confirming that the WR component is a very hot (a parts per thousand 90 kK) WN star. For this star, we derived a luminosity of log L/ L-aS (TM) = 5.45 and a mass-loss rate of 10(- 5.8) M-aS (TM) yr(- 1), and found that the stellar wind composition is dominated by helium with 20 per cent of hydrogen. Spectroscopy of the shell revealed an He iii region centred on BAT99 3a and having the same angular radius (a parts per thousand 15 arcsec) as the shell. We thereby add a new example to a rare class of high-excitation nebulae photoionized by WR stars. Analysis of the nebular spectrum showed that the shell is composed of unprocessed material, implying that the shell was swept-up from the local interstellar medium. We discuss the physical relationship between the newly identified massive stars and their possible membership of a previously unrecognized star cluster.
Context. Massive Wolf-Rayet stars are recognized today to be in a very common, but short, evolutionary phase of massive stars. While our understanding of Wolf-Rayet stars has increased dramatically over the past decades, it remains unclear whether rapid rotators are among them. There are various indications that rapidly rotating Wolf-Rayet stars should exist. Unfortunately, due to their expanding atmospheres, rotational velocities of Wolf-Rayet stars are very difficult to measure. However, recently observed spectra of several Wolf-Rayet stars reveal peculiarly broad and round emission lines. Could these spectra imply rapid rotation?
Aims. In this work, we model the effects of rotation on the atmospheres of Wolf-Rayet stars. We further investigate whether the peculiar spectra of five Wolf-Rayet stars may be explained with the help of stellar rotation, infer appropriate rotation parameters, and discuss the implications of our results.
Methods. We make use of the Potsdam Wolf-Rayet (PoWR) non-LTE model atmosphere code. Since the observed spectra of WolfRayet stars are mainly formed in their expanding atmospheres, rotation must be accounted for with a 3D integration scheme of the formal integral. For this purpose, we assume a rotational velocity field consisting of an inner co-rotating domain and an outer domain, where the angular momentum is conserved.
Results. We find that rotation can reproduce the unique spectra analyzed here. However, the inferred rotational velocities at the stellar surface are large (similar to 200 km s(-1)), and the inferred co-rotation radii (similar to 10R.) suggest the existence of very strong photospheric magnetic fields (similar to 20 kG).
We present an overview of four deep phase-constrained Chandra HETGS X-ray observations of delta Ori A. Delta Ori A is actually a triple system that includes the nearest massive eclipsing spectroscopic binary, delta Ori Aa, the only such object that can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, delta Ori Aa2, has a much lower X-ray luminosity than the brighter primary (delta Ori Aa1), delta Ori Aa provides a unique system with which to test the spatial distribution of the X-ray emitting gas around delta Ori Aa1 via occultation by the photosphere of, and wind cavity around, the X-ray dark secondary. Here we discuss the X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time of nearly 500 ks and covering nearly the entire binary orbit. The companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities obtained simultaneously with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectra. We find that the X-ray emission is dominated by embedded wind shock emission from star Aa1, with little contribution from the tertiary star Ab or the shocked gas produced by the collision of the wind of Aa1 against the surface of Aa2. We find a similar temperature distribution to previous X-ray spectrum analyses. We also show that the line half-widths are about 0.3-0.5 times the terminal velocity of the wind of star Aa1. We find a strong anti-correlation between line widths and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the strong and weak lines of Fe XVII and Ne X are inconsistent with model predictions, which may be an effect of resonance scattering.
The distribution of angular momentum in massive stars is a critical component of their evolution, yet not much is known on the rotation velocities of Wolf-Rayet stars. There are various indications that rapidly rotating Wolf-Rayet stars should exist. Unfortunately, due to their expanding atmospheres, rotational velocities of Wolf-Rayet stars are very difficult to measure. In this work, we model the effects of rotation on the atmospheres of Wolf-Rayet stars by implementing a 3D integration scheme in the PoWR code. We further investigate whether the peculiar spectra of five Wolf-Rayet stars may imply rapid rotation, infer the corresponding rotation parameters, and discuss the implications of our results. We find that rotation helps to reproduce the unique spectra analyzed here. However, if rotation is indeed involved, the inferred rotational velocities at the stellar surface are large (∼ 200 km/s), and the implied co-rotation radii (∼ 10R∗) suggest the existence of very strong photospheric magnetic fields (∼ 20 kG).
In the last decades, stellar atmosphere codes have become a key tool in understanding massive stars, including precise calculations of stellar and wind parameters, such as temperature, massloss rate, and terminal wind velocity. Nevertheless, for these models the hydrodynamic equation is not solved in the wind. Motivated by the results of the CAK theory, the models typically use a beta velocity law, which however turns out not to be adequate for stars with very strong winds, and treat the mass-loss rate as a free parameter. In a new branch of the Potsdam Wolf-Rayet model atmosphere (PoWR) code, we solve the hydrodynamic equation consistently throughout the stellar atmosphere. The PoWR code performs the calculation of the radiative force without approximations (e.g. Sobolev). We show the impact of hydrodynamically consistent modelling on OB and WR stars in comparison to conventional models and discuss the obtained velocity fields and their impact on the observed spectral lines.
Macroclumping in WR 136
(2015)
Macroclumping proved to resolve the discordance between different mass-loss rate diagnostics for O-type stars, in particular between Hα and the P v resonance lines. In this paper, we report first results from a corresponding investigation for WR stars. We apply our detailed 3-D Monte Carlo (MC) line formation code to the P v resonance doublet and show, for the Galactic WNL star WR136, that macroclumping is require to bring this line in accordance with the mass-loss rate derived from the emission-line spectrum.
Context. Wolf-Rayet (WR) stars have a severe impact on their environments owing to their strong ionizing radiation fields and powerful stellar winds. Since these winds are considered to be driven by radiation pressure, it is theoretically expected that the degree of the wind mass-loss depends on the initial metallicity of WR stars.
Aims. Following our comprehensive studies of WR stars in the Milky Way, M31, and the LMC, we derive stellar parameters and mass-loss rates for all seven putatively single WN stars known in the SMC. Based on these data, we discuss the impact of a low-metallicity environment on the mass loss and evolution of WR stars.
Methods. The quantitative analysis of the WN stars is performed with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The physical properties of our program stars are obtained from fitting synthetic spectra to multi-band observations.
Results. In all SMC WN stars, a considerable surface hydrogen abundance is detectable. The majority of these objects have stellar temperatures exceeding 75 kK, while their luminosities range from 10(5.5) to 10(6.1) L-circle dot. The WN stars in the SMC exhibit on average lower mass-loss rates and weaker winds than their counterparts in the Milky Way, M31, and the LMC.
Conclusions. By comparing the mass-loss rates derived for WN stars in different Local Group galaxies, we conclude that a clear dependence of the wind mass-loss on the initial metallicity is evident, supporting the current paradigm that WR winds are driven by radiation. A metallicity effect on the evolution of massive stars is obvious from the HRD positions of the SMC WN stars at high temperatures and high luminosities. Standard evolution tracks are not able to reproduce these parameters and the observed surface hydrogen abundances. Homogeneous evolution might provide a better explanation for their evolutionary past.
Context. Spectroscopic analysis remains the most common method to derive masses of massive stars, the most fundamental stellar parameter. While binary orbits and stellar pulsations can provide much sharper constraints on the stellar mass, these methods are only rarely applicable to massive stars. Unfortunately, spectroscopic masses of massive stars heavily depend on the detailed physics of model atmospheres.
Aims. We demonstrate the impact of a consistent treatment of the radiative pressure on inferred gravities and spectroscopic masses of massive stars. Specifically, we investigate the contribution of line and continuum transitions to the photospheric radiative pressure. We further explore the effect of model parameters, e.g., abundances, on the deduced spectroscopic mass. Lastly, we compare our results with the plane-parallel TLUSTY code, commonly used for the analysis of massive stars with photospheric spectra.
Methods. We calculate a small set of O-star models with the Potsdam Wolf-Rayet (PoWR) code using different approaches for the quasi-hydrostatic part. These models allow us to quantify the effect of accounting for the radiative pressure consistently. We further use PoWR models to show how the Doppler widths of line profiles and abundances of elements such as iron affect the radiative pressure, and, as a consequence, the derived spectroscopic masses.
Results. Our study implies that errors on the order of a factor of two in the inferred spectroscopic mass are to be expected when neglecting the contribution of line and continuum transitions to the radiative acceleration in the photosphere. Usage of implausible microturbulent velocities, or the neglect of important opacity sources such as Fe, may result in errors of approximately 50% in the spectroscopic mass. A comparison with TLUSTY model atmospheres reveals a very good agreement with PoWR at the limit of low mass-loss rates.
We present new grids of Potsdam Wolf-Rayet (PoWR) model atmospheres for Wolf-Rayet stars of the nitrogen sequence (WN stars). The models have been calculated with the latest version of the PoWR stellar atmosphere code for spherical stellar winds. The WN model atmospheres include the non-LTE solutions of the statistical equations for complex model atoms, as well as the radiative transfer equation in the co-moving frame. Iron-line blanketing is treated with the help of the superlevel approach, while wind inhomogeneities are taken into account via optically thin clumps. Three of our model grids are appropriate for Galactic metallicity. The hydrogen mass fraction of these grids is 50%, 20%, and 0%, thus also covering the hydrogen-rich late-type WR stars that have been discovered in recent years. Three grids are adequate for LMC WN stars and have hydrogen fractions of 40%, 20%, and 0%. Recently, additional grids with SMC metallicity and with 60%, 40%, 20%, and 0% hydrogen have been added. We provide contour plots of the equivalent widths of spectral lines that are usually used for classification and diagnostics.
We report on both high-precision photometry from the Microvariability and Oscillations of Stars (MOST) space telescope and ground-based spectroscopy of the triple system delta Ori A, consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7 days, and a more distant tertiary (O9 IV P > 400 years). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for three weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy, we have a well-constrained radial velocity (RV) curve of Aa1. While we are unable to recover RV variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries that show evidence for tidally induced pulsations.
We present an analysis of a 78 ks Chandra high-energy transmission gratings observation of the B0I star QV Nor, the massive donor of the wind-accreting pulsar 4U1538-52. The neutron star (NS) orbits its companion in a very close orbit (r < 1.4R(*), in units of the stellar radii), thereby allowing probing of the innermost wind regions. The flux of the Fe K alpha line during eclipse reduces to only similar to 30% of the flux measured out of eclipse. This indicates that the majority of Fe fluorescence must be produced in regions close to the NS, at distances smaller than 1R(*) from its surface. The fact that the flux of the continuum decreases to only similar to 3% during eclipse allows for a high contrast of the Fe Ka line fluorescence during eclipse. The line is not resolved and centered at lambda = 1.9368(-0.0018)(+0.0032) angstrom. From the inferred plasma speed limit of v < c Delta lambda/lambda < 800 km s(-1) and range of ionization parameters of log xi =[-1, 2], together with the stellar density profile, we constrain the location of the cold, dense material in the stellar wind of QV Nor using simple geometrical considerations. We then use the Fe K alpha line fluorescence as a tracer of wind clumps and determine that these clumps in the stellar wind of QV Nor (B0I) must already be present at radii r < 1.25R(*), close to the photosphere of the star.
Eclipsing systems of massive stars allow one to explore the properties of their components in great detail. We perform a multi-wavelength, non-LTE analysis of the three components of the massive multiple system delta Ori A, focusing on the fundamental stellar properties, stellar winds, and X-ray characteristics of the system. The primary's distance-independent parameters turn out to be characteristic for its spectral type (O9.5 II), but usage of the Hipparcos parallax yields surprisingly low values for the mass, radius, and luminosity. Consistent values follow only if delta Ori lies at about twice the Hipparcos distance, in the vicinity of the sigma-Orionis cluster. The primary and tertiary dominate the spectrum and leave the secondary only marginally detectable. We estimate the V-band magnitude difference between primary and secondary to be Delta V approximate to 2.(m)8. The inferred parameters suggest that the secondary is an early B-type dwarf (approximate to B1 V), while the tertiary is an early B-type subgiant (approximate to B0 IV). We find evidence for rapid turbulent velocities (similar to 200 km s(-1)) and wind inhomogeneities, partially optically thick, in the primary's wind. The bulk of the X-ray emission likely emerges from the primary's stellar wind (logL(X)/L-Bol approximate to -6.85), initiating close to the stellar surface at R-0 similar to 1.1 R-*. Accounting for clumping, the mass-loss rate of the primary is found to be log (M) over dot approximate to -6.4 (M-circle dot yr(-1))., which agrees with hydrodynamic predictions, and provides a consistent picture along the X-ray, UV, optical, and radio spectral domains.
With a deep Chandra/HETGS exposure of WR 6, we have resolved emission lines whose profiles show that the X-rays originate from a uniformly expanding spherical wind of high X-ray-continuum optical depth. The presence of strong helium-like forbidden lines places the source of X-ray emission at tens to hundreds of stellar radii from the photosphere. Variability was present in X-rays and simultaneous optical photometry, but neither were correlated with the known period of the system or with each other. An enhanced abundance of sodium revealed nuclear-processed material, a quantity related to the evolutionary state of the star. The characterization of the extent and nature of the hot plasma in WR 6 will help to pave the way to a more fundamental theoretical understanding of the winds and evolution of massive stars.
We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution spectral data set of the delta Ori Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of approximate to 479 ks and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range of 5-25 is is confirmed, with a maximum amplitude of about +/- 15% within a single approximate to 125 ks observation. Periods of 4.76 and 2.04 days are found in the total X-ray flux, as well as an apparent overall increase in the flux level throughout the nine-day observational campaign. Using 40 ks contiguous spectra derived from the original observations, we investigate the variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at phi = 0.0 when the secondary delta Ori Aa2 is at the inferior conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission line width variability to the presence of a wind cavity created by a wind-wind collision, which is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary wind, thus producing phase-locked X-ray variability.
During the summer of 2013, a 4-month spectroscopic campaign took place to observe the variabilities in three Wolf-Rayet stars. The spectroscopic data have been analysed for WR 134 (WN6b), to better understand its behaviour and long-term periodicity, which we interpret as arising from corotating interaction regions (CIRs) in the wind. By analysing the variability of the He ii lambda 5411 emission line, the previously identified period was refined to P = 2.255 +/- 0.008 (s.d.) d. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 +/- 6 d, or similar to 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grey-scale difference images with theoretical grey-scales previously calculated from models including CIRs in an optically thin stellar wind, we find that two CIRs were likely present. A separation in longitude of Delta I center dot a parts per thousand integral 90A degrees was determined between the two CIRs and we suggest that the different maximum velocities that they reach indicate that they emerge from different latitudes. We have also been able to detect observational signatures of the CIRs in other spectral lines (C iv lambda lambda 5802,5812 and He i lambda 5876). Furthermore, a DAC was found to be present simultaneously with the CIR signatures detected in the He i lambda 5876 emission line which is consistent with the proposed geometry of the large-scale structures in the wind. Small-scale structures also show a presence in the wind, simultaneously with the larger scale structures, showing that they do in fact co-exist.
We report on interferometric observations with the CHARAArray of two classical Wolf-Rayet (WR) stars in suspected binary systems, namely WR 137 and WR 138. In both cases, we resolve the component stars to be separated by a few milliarcseconds. The data were collected in the H band, and provide a measure of the fractional flux for both stars in each system. We find that the WR star is the dominant H-band light source in both systems (fWR, 137 = 0.59 +/- 0.04; fWR, 138 = 0.67 +/- 0.01), which is confirmed through both comparisons with estimated fundamental parameters for WR stars and O dwarfs, as well as through spectral modelling of each system. Our spectral modelling also provides fundamental parameters for the stars and winds in these systems. The results on WR 138 provide evidence that it is a binary system which may have gone through a previous mass-transfer episode to create the WR star. The separation and position of the stars in the WR 137 system together with previous results from the IOTA interferometer provides evidence that the binary is seen nearly edgeon. The possible edge-on orbit of WR 137 aligns well with the dust production site imaged by the Hubble Space Telescope during a previous periastron passage, showing that the dust production may be concentrated in the orbital plane.
Aims. To close this gap, we perform a comparative analysis of the optical companion in two important systems: IGR J175442619 (SFXT) and Vela X-1 (SGXB). We analyze the spectra of each star in detail and derive their stellar and wind properties. As a next step, we compare the wind parameters, giving us an excellent chance of recognizing key differences between donor winds in SFXTs and SGXBs. Methods. We use archival infrared, optical and ultraviolet observations, and analyze them with the non-local thermodynamic equilibrium (NLTE) Potsdam Wolf-Rayet model atmosphere code. We derive the physical properties of the stars and their stellar winds, accounting for the influence of X-rays on the stellar winds. Results. We find that the stellar parameters derived from the analysis generally agree well with the spectral types of the two donors: O9I (IGR J17544-2619) and B0.5Iae (Vela X-1). The distance to the sources have been revised and also agree well with the estimations already available in the literature. In IGR J17544-2619 we are able to narrow the uncertainty to d = 3.0 +/- 0.2 kpc. From the stellar radius of the donor and its X-ray behavior, the eccentricity of IGR J17544-2619 is constrained to e < 0.25. The derived chemical abundances point to certain mixing during the lifetime of the donors. An important difference between the stellar winds of the two stars is their terminal velocities (v(infinity) = 1500 km s(-1) in IGR J17544-2619 and v(infinity) = 700 km s(-1) in Vela X-1), which have important consequences on the X-ray luminosity of these sources. Conclusions. The donors of IGR J17544-2619 and Vela X-1 have similar spectral types as well as similar parameters that physically characterize them and their spectra. In addition, the orbital parameters of the systems are similar too, with a nearly circular orbit and short orbital period. However, they show moderate differences in their stellar wind velocity and the spin period of their neutron star which has a strong impact on the X-ray luminosity of the sources. This specific combination of wind speed and pulsar spin favors an accretion regime with a persistently high luminosity in Vela X-1, while it favors an inhibiting accretion mechanism in IGR J17544-2619. Our study demonstrates that the relative wind velocity is critical in class determination for the HMXBs hosting a supergiant donor, given that it may shift the accretion mechanism from direct accretion to propeller regimes when combined with other parameters.
Context. Massive Wolf-Rayet (WR) stars are evolved massive stars (M-i greater than or similar to 20 M-circle dot) characterized by strong mass-loss. Hypothetically, they can form either as single stars or as mass donors in close binaries. About 40% of all known WR stars are confirmed binaries, raising the question as to the impact of binarity on the WR population. Studying WR binaries is crucial in this context, and furthermore enable one to reliably derive the elusive masses of their components, making them indispensable for the study of massive stars. Aims. By performing a spectral analysis of all multiple WR systems in the Small Magellanic Cloud (SMC), we obtain the full set of stellar parameters for each individual component. Mass-luminosity relations are tested, and the importance of the binary evolution channel is assessed. Methods. The spectral analysis is performed with the PotsdamWolf-Rayet (PoWR) model atmosphere code by superimposing model spectra that correspond to each component. Evolutionary channels are constrained using the Binary Population and Spectral Synthesis (BPASS) evolution tool. Results. Significant hydrogen mass fractions (0.1 < X-H < 0.4) are detected in all WN components. A comparison with mass-luminosity relations and evolutionary tracks implies that the majority of the WR stars in our sample are not chemically homogeneous. The WR component in the binary AB6 is found to be very luminous (log L approximate to 6.3 [L-circle dot]) given its orbital mass (approximate to 10 M-circle dot), presumably because of observational contamination by a third component. Evolutionary paths derived for our objects suggest that Roche lobe overflow had occurred in most systems, affecting their evolution. However, the implied initial masses (greater than or similar to 60 M-circle dot) are large enough for the primaries to have entered the WR phase, regardless of binary interaction. Conclusions. Together with the results for the putatively single SMC WR stars, our study suggests that the binary evolution channel does not dominate the formation of WR stars at SMC metallicity.
Via their powerful radiation, stellar winds, and supernova explosions, massive stars (Mini & 8 M☉) bear a tremendous impact on galactic evolution. It became clear in recent decades that the majority of massive stars reside in binary systems. This thesis sets as a goal to quantify the impact of binarity (i.e., the presence of a companion star) on massive stars. For this purpose, massive binary systems in the Local Group, including OB-type binaries, high mass X-ray binaries (HMXBs), and Wolf-Rayet (WR) binaries, were investigated by means of spectral, orbital, and evolutionary analyses.
The spectral analyses were performed with the non-local thermodynamic equillibrium (non-LTE) Potsdam Wolf-Rayet (PoWR) model atmosphere code. Thanks to critical updates in the calculation of the hydrostatic layers, the code became a state-of-the-art tool applicable for all types of hot massive stars (Chapter 2). The eclipsing OB-type triple system δ Ori served as an intriguing test-case for the new version of the PoWR code, and provided key insights regarding the formation of X-rays in massive stars (Chapter 3). We further analyzed two prototypical HMXBs, Vela X-1 and IGR J17544-2619, and obtained fundamental conclusions regarding the dichotomy of two basic classes of HMXBs (Chapter 4). We performed an exhaustive analysis of the binary R 145 in the Large Magellanic Cloud (LMC), which was claimed to host the most massive stars known. We were able to disentangle the spectrum of the system, and performed an orbital, polarimetric, and spectral analysis, as well as an analysis of the wind-wind collision region. The true masses of the binary components turned out to be significantly lower than suggested, impacting our understanding of the initial mass function and stellar evolution at low metallicity (Chapter 5). Finally, all known WR binaries in the Small Magellanic Cloud (SMC) were analyzed. Although it was theoretical predicted that virtually all WR stars in the SMC should be formed via mass-transfer in binaries, we find that binarity was not important for the formation of the known WR stars in the SMC, implying a strong discrepancy between theory and observations (Chapter 6).
Context:
Massive binaries play a crucial role in the Universe. Knowing the distributions of their orbital parameters is important for a wide range of topics from stellar feedback to binary evolution channels and from the distribution of supernova types to gravitational wave progenitors, yet no direct measurements exist outside the Milky Way.
Aims:
The Tarantula Massive Binary Monitoring project was designed to help fill this gap by obtaining multi-epoch radial velocity (RV) monitoring of 102 massive binaries in the 30 Doradus region.
Methods:
In this paper we analyze 32 FLAMES/GIRAFFE observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single-lined (SB1) and 31 double-lined (SB2) spectroscopic binaries.
Results:
Overall, the binary fraction and orbital properties across the 30 Doradus region are found to be similar to existing Galactic samples. This indicates that within these domains environmental effects are of second order in shaping the properties of massive binary systems. A small difference is found in the distribution of orbital periods, which is slightly flatter (in log space) in 30 Doradus than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, orbital periods in 30 Doradus can be as short as 1.1 d, somewhat shorter than seen in Galactic samples. Equal mass binaries (q> 0.95) in 30 Doradus are all found outside NGC 2070, the central association that surrounds R136a, the very young and massive cluster at 30 Doradus’s core. Most of the differences, albeit small, are compatible with expectations from binary evolution. One outstanding exception, however, is the fact that earlier spectral types (O2–O7) tend to have shorter orbital periods than later spectral types (O9.2–O9.7).
Conclusions:
Our results point to a relative universality of the incidence rate of massive binaries and their orbital properties in the metallicity range from solar (Z⊙) to about half solar. This provides the first direct constraints on massive binary properties in massive star-forming galaxies at the Universe’s peak of star formation at redshifts z ~ 1 to 2 which are estimated to have Z ~ 0.5 Z⊙.