@article{ParkinBroosTownsleyetal.2011,
  author    = {Parkin, E. R. and Broos, Patrick S. and Townsley, L. K. and Pittard, J. M. and Moffat, Anthony F. J. and Naze, Y. and Rauw, G. and Oskinova, Lida and Waldron, W. L.},
  title     = {X-RAY EMISSION FROM THE DOUBLE-BINARY OB-STAR SYSTEM QZ CAR (HD 93206)},
  series = {ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES},
  volume    = {194},
  journal   = {ASTROPHYSICAL JOURNAL SUPPLEMENT SERIES},
  number    = {1},
  publisher = {IOP PUBLISHING LTD},
  address   = {BRISTOL},
  issn      = {0067-0049},
  doi       = {10.1088/0067-0049/194/1/8},
  pages     = {11},
  year      = {2011},
  abstract  = {X-ray observations of the double-binary OB-star system QZ Car (HD 93206) obtained with the Chandra X-ray Observatory over a period of roughly 2 years are presented. The respective orbits of systems A (O9.7 I+b2 v, P-A = 21 days) and B (O8 III+o9 v, P-B = 6 days) are reasonably well sampled by the observations, allowing the origin of the X-ray emission to be examined in detail. The X-ray spectra can be well fitted by an attenuated three-temperature thermal plasma model, characterized by cool, moderate, and hot plasma components at kT similar or equal to 0.2, 0.7, and 2 keV, respectively, and a circumstellar absorption of similar or equal to 0.2 x 10(22) cm(-2). Although the hot plasma component could be indicating the presence of wind-wind collision shocks in the system, the model fluxes calculated from spectral fits, with an average value of similar or equal to 7x10(-13) erg s(-1) cm(-2), do not show a clear correlation with the orbits of the two constituent binaries. A semi-analytical model of QZ Car reveals that a stable momentum balance may not be established in either system A or B. Yet, despite this, system B is expected to produce an observed X-ray flux well in excess of the observations. If one considers the wind of the O8 III star to be disrupted by mass transfer, the model and observations are in far better agreement, which lends support to the previous suggestion of mass transfer in the O8 III+o9 v binary. We conclude that the X-ray emission from QZ Car can be reasonably well accounted for by a combination of contributions mainly from the single stars and the mutual wind-wind collision between systems A and B.},
  language  = {en}
}
@article{BurgemeisterGvaramadzeStringfellowetal.2013,
  author    = {Burgemeister, S. and Gvaramadze, Visily V. and Stringfellow, G. S. and Kniazev, Alexei Y. and Todt, Helge Tobias and Hamann, Wolf-Rainer},
  title     = {WR 120bb and WR 120bc: a pair of WN9h stars with possibly interacting circumstellar shells},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {429},
  journal   = {Monthly notices of the Royal Astronomical Society},
  number    = {4},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/sts588},
  pages     = {3305 -- 3315},
  year      = {2013},
  abstract  = {Two optically obscured Wolf-Rayet (WR) stars have been recently discovered by means of their infrared (IR) circumstellar shells, which show signatures of interaction with each other. Following the systematics of the WR star catalogues, these stars obtain the names WR 120bb and WR 120bc. In this paper, we present and analyse new near-IR, J-, H- and K-band spectra using the Potsdam Wolf-Rayet model atmosphere code. For that purpose, the atomic data base of the code has been extended in order to include all significant lines in the near-IR bands. The spectra of both stars are classified as WN9h. As their spectra are very similar the parameters that we obtained by the spectral analyses hardly differ. Despite their late spectral subtype, we found relatively high stellar temperatures of 63 kK. The wind composition is dominated by helium, while hydrogen is depleted to 25 per cent by mass. Because of their location in the Scutum-Centaurus Arm, WR 120bb and WR 120bc appear highly reddened, A(Ks) approximate to 2 mag. We adopt a common distance of 5.8 kpc to both stars, which complies with the typical absolute K-band magnitude for the WN9h subtype of -6.5 mag, is consistent with their observed extinction based on comparison with other massive stars in the region, and allows for the possibility that their shells are interacting with each other. This leads to luminosities of log(L/L-circle dot) = 5.66 and 5.54 for WR 120bb and WR 120bc, with large uncertainties due to the adopted distance. The values of the luminosities of WR 120bb and WR 120bc imply that the immediate precursors of both stars were red supergiants (RSG). This implies in turn that the circumstellar shells associated with WR 120bb and WR 120bc were formed by interaction between the WR wind and the dense material shed during the preceding RSG phase.},
  language  = {en}
}
@article{ShenarHainichTodtetal.2016,
  author    = {Shenar, Tomer and Hainich, Rainer and Todt, Helge Tobias and Sander, Andreas Alexander Christoph and Hamann, Wolf-Rainer and Moffat, Anthony F. J. and Eldridge, J. J. and Pablo, H. and Oskinova, Lida and Richardson, N. D.},
  title     = {Wolf-Rayet stars in the Small Magellanic Cloud II. Analysis of the binaries},
  series = {American mineralogist : an international journal of earth and planetary materials},
  volume    = {591},
  journal   = {American mineralogist : an international journal of earth and planetary materials},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201527916},
  pages     = {25},
  year      = {2016},
  abstract  = {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.},
  language  = {en}
}
@article{MeyerPetrovPohl2020,
  author    = {Meyer, Dominique M.-A. and Petrov, Mykola and Pohl, Martin},
  title     = {Wind nebulae and supernova remnants of very massive stars},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {493},
  journal   = {Monthly notices of the Royal Astronomical Society},
  number    = {3},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/staa554},
  pages     = {3548 -- 3564},
  year      = {2020},
  abstract  = {A very small fraction of (runaway) massive stars have masses exceeding 60-70 M-circle dot and are predicted to evolve as luminous blue variable and Wolf-Rayet stars before ending their lives as core-collapse supernovae. Our 2D axisymmetric hydrodynamical simulations explore how a fast wind (2000 km s(-1)) and high mass-loss rate (10(-5)M(circle dot) yr(-1)) can impact the morphology of the circumstellar medium. It is shaped as 100 pc-scale wind nebula that can be pierced by the driving star when it supersonically moves with velocity 20-40 km s(-1) through the interstellar medium (ISM) in the Galactic plane. The motion of such runaway stars displaces the position of the supernova explosion out of their bow shock nebula, imposing asymmetries to the eventual shock wave expansion and engendering Cygnus-loop-like supernova remnants. We conclude that the size (up to more than 200 pc) of the filamentary wind cavity in which the chemically enriched supernova ejecta expand, mixing efficiently the wind and ISM materials by at least 10 per cent in number density, can be used as a tracer of the runaway nature of the very massive progenitors of such 0.1Myr old remnants. Our results motivate further observational campaigns devoted to the bow shock of the very massive stars BD+43 degrees 3654 and to the close surroundings of the synchrotron-emitting Wolf-Rayet shell G2.4+1.4.},
  language  = {en}
}
@article{KurfuerstFeldmeierKrticka2018,
  author    = {Kurf{\"u}rst, P. and Feldmeier, Achim and Krticka, Jiri},
  title     = {Two-dimensional modeling of density and thermal structure of dense circumstellar outflowing disks},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {613},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201731300},
  pages     = {24},
  year      = {2018},
  abstract  = {Context. Evolution of massive stars is affected by a significant loss of mass either via (nearly) spherically symmetric stellar winds or by aspherical mass-loss mechanisms, namely the outflowing equatorial disks. However, the scenario that leads to the formation of a disk or rings of gas and dust around massive stars is still under debate. It is also unclear how various forming physical mechanisms of the circumstellar environment affect its shape and density, as well as its kinematic and thermal structure. Results. Our models show the geometric distribution and contribution of viscous heating that begins to dominate in the central part of the disk for mass-loss rates higher than (M) over dot greater than or similar to 10(-10) M-circle dot yr(-1). In the models of dense viscous disks with (M) over dot > 10(-8) M-circle dot yr(-1), the viscosity increases the central temperature up to several tens of thousands of Kelvins, however the temperature rapidly drops with radius and with distance from the disk midplane. The high mass-loss rates and high viscosity lead to instabilities with significant waves or bumps in density and temperature in the very inner disk region. Conclusions. The two-dimensional radial-vertical models of dense outflowing disks including the full Navier-Stokes viscosity terms show very high temperatures that are however limited to only the central disk cores inside the optically thick area, while near the edge of the optically thick region the temperature may be low enough for the existence of neutral hydrogen, for example.},
  language  = {en}
}
@article{GomezMoranOskinova2018,
  author    = {Gomez-Moran, Ada Nebot and Oskinova, Lida},
  title     = {The X-ray catalog of spectroscopically identified Galactic O stars Investigating the dependence of X-ray luminosity on stellar and wind parameters},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {620},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201833453},
  pages     = {18},
  year      = {2018},
  abstract  = {The X-ray emission of O-type stars was first discovered in the early days of the Einstein satellite. Since then many different surveys have confirmed that the ratio of X-ray to bolometric luminosity in O-type stars is roughly constant, but there is a paucity of studies that account for detailed information on spectral and wind properties of O-stars. Recently a significant sample of O stars within our Galaxy was spectroscopically identified and presented in the Galactic O-Star Spectroscopic Survey (GOSS). At the same time, a large high-fidelity catalog of X-ray sources detected by the XMM-Newton X-ray telescope was released. Here we present the X-ray catalog of O stars with known spectral types and investigate the dependence of their X-ray properties on spectral type as well as stellar and wind parameters. We find that, among the GOSS sample, 127 O-stars have a unique XMM-Newton source counterpart and a Gaia data release 2 (DR2) association. Terminal velocities are known for a subsample of 35 of these stars. We confirm that the X-ray luminosities of dwarf and giant O stars correlate with their bolometric luminosity. For the subsample of O stars with measure terminal velocities we find that the X-ray luminosities of dwarf and giant O stars also correlate with wind parameters. However, we find that these correlations break down for supergiant stars. Moreover, we show that supergiant stars are systematically harder in X-rays compared to giant and dwarf O-type stars. We find that the X-ray luminosity depends on spectral type, but seems to be independent of whether the stars are single or in a binary system. Finally, we show that the distribution of log(L-X/L-bol) in our sample stars is non-Gaussian, with the peak of the distribution at log(L-X/L-bol) approximate to -6.6.},
  language  = {en}
}
@article{SanderTodtHainichetal.2014,
  author    = {Sander, Andreas Alexander Christoph and Todt, Helge Tobias and Hainich, Rainer and Hamann, Wolf-Rainer},
  title     = {The Wolf-Rayet stars in M31 I. Analysis of the late-type WN stars},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {563},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {0004-6361},
  doi       = {10.1051/0004-6361/201323240},
  pages     = {22},
  year      = {2014},
  abstract  = {Context. Comprehensive studies of Wolf-Rayet stars were performed in the past for the Galactic and the LMC population. The results revealed significant differences, but also unexpected similarities between the WR populations of these different galaxies. Analyzing the WR stars in M 31 will extend our understanding of these objects in different galactic environments. Aims. The present study aims at the late-type WN stars in M 31. The stellar and wind parameters will tell about the formation of WR stars in other galaxies with different metallicity and star formation histories. The obtained parameters will provide constraints to the evolution of massive stars in the environment of M 31. Methods. We used the latest version of the Potsdam Wolf-Rayet model atmosphere code to analyze the stars via fitting optical spectra and photometric data. To account for the relatively low temperatures of the late WN I 0 and WN I I subtypes, our WN models have been extended into this temperature regime. Results. Stellar and atmospheric parameters are derived for all known late-type WN stars in M 31 with available spectra. All of these stars still have hydrogen in their outer envelopes, some of them up to 50\% by mass. The stars are located on the cool side of the zero age main sequence in the Hertzsprung-Russell diagram, while their luminosities range from 105 to 1064). It is remarkable that no star exceeds 106 L. Conclusions. If formed via single-star evolution, the late-type WN stars in M 31 stem from an initial mass range between 20 and 60 M-circle dot. From the very late-type WN9-11 stars, only one star is located in the S Doradus instability strip. We do not find any late-type WN stars with the high luminosities known in the Milky Way.},
  language  = {en}
}
@article{ShenarSablowskiHainichetal.2019,
  author    = {Shenar, Tomer and Sablowski, D. P. and Hainich, Rainer and Todt, Helge Tobias and Moffat, Anthony F. J. and Oskinova, Lida and Ramachandran, Varsha and Sana, Hugues and Sander, Andreas Alexander Christoph and Schnurr, O. and St-Louis, N. and Vanbeveren, D. and Gotberg, Y. and Hamann, Wolf-Rainer},
  title     = {The Wolf-Rayet binaries of the nitrogen sequence in the Large Magellanic Cloud Spectroscopy, orbital analysis, formation, and evolution},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {627},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {0004-6361},
  doi       = {10.1051/0004-6361/201935684},
  pages     = {68},
  year      = {2019},
  abstract  = {Context. Massive Wolf-Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core collapse. It is not known whether core He-burning WR stars (classical WR; cWR) form predominantly through wind stripping (w-WR) or binary stripping (b-WR). Whereas spectroscopy of WR binaries has so-far largely been avoided because of its complexity, our study focuses on the 44 WR binaries and binary candidates of the Large Magellanic Cloud (LMC; metallicity Z approximate to 0.5 Z(circle dot)), which were identified on the basis of radial velocity variations, composite spectra, or high X-ray luminosities. Aims. Relying on a diverse spectroscopic database, we aim to derive the physical and orbital parameters of our targets, confronting evolution models of evolved massive stars at subsolar metallicity and constraining the impact of binary interaction in forming these stars. Methods. Spectroscopy was performed using the Potsdam Wolf-Rayet (PoWR) code and cross-correlation techniques. Disentanglement was performed using the code Spectangular or the shift-and-add algorithm. Evolutionary status was interpreted using the Binary Population and Spectral Synthesis (BPASS) code, exploring binary interaction and chemically homogeneous evolution. Results. Among our sample, 28/44 objects show composite spectra and are analyzed as such. An additional five targets show periodically moving WR primaries but no detected companions (SB1); two (BAT99 99 and 112) are potential WR + compact-object candidates owing to their high X-ray luminosities. We cannot confirm the binary nature of the remaining 11 candidates. About two-thirds of the WN components in binaries are identified as cWR, and one-third as hydrogen-burning WR stars. We establish metallicity-dependent mass-loss recipes, which broadly agree with those recently derived for single WN stars, and in which so-called WN3/O3 stars are clear outliers. We estimate that 45 +/- 30\% of the cWR stars in our sample have interacted with a companion via mass transfer. However, only approximate to 12 +/- 7\% of the cWR stars in our sample naively appear to have formed purely owing to stripping via a companion (12\% b-WR). Assuming that apparently single WR stars truly formed as single stars, this comprises approximate to 4\% of the whole LMC WN population, which is about ten times less than expected. No obvious differences in the properties of single and binary WN stars, whose luminosities extend down to log L approximate to 5.2 [L-circle dot], are apparent. With the exception of a few systems (BAT99 19, 49, and 103), the equatorial rotational velocities of the OB-type companions are moderate (v(eq) less than or similar to 250 km s(-1)) and challenge standard formalisms of angular-momentum accretion. For most objects, chemically homogeneous evolution can be rejected for the secondary, but not for the WR progenitor. Conclusions. No obvious dichotomy in the locations of apparently single and binary WN stars on the Hertzsprung-Russell diagram is apparent. According to commonly used stellar evolution models (BPASS, Geneva), most apparently single WN stars could not have formed as single stars, implying that they were stripped by an undetected companion. Otherwise, it must follow that pre-WR mass-loss/mixing (e.g., during the red supergiant phase) are strongly underestimated in standard stellar evolution models.},
  language  = {en}
}
@article{AlmeidaSanaTayloretal.2017,
  author    = {Almeida, Leonardo A. and Sana, H. and Taylor, W. and Barb{\´a}, Rodolfo and Bonanos, Alceste Z. and Crowther, Paul and Damineli, Augusto and de Koter, A. and de Mink, Selma E. and Evans, C. J. and Gieles, Mark and Grin, Nathan J. and H{\´e}nault-Brunet, V. and Langer, Norbert and Lennon, D. and Lockwood, Sean and Ma{\´i}z Apell{\´a}niz, Jes{\´u}s and Moffat, A. F. J. and Neijssel, C. and Norman, C. and Ram{\´i}rez-Agudelo, O. H. and Richardson, N. D. and Schootemeijer, Abel and Shenar, Tomer and Soszyński, Igor and Tramper, Frank and Vink, J. S.},
  title     = {The tarantula massive binary monitoring},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {598},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201629844},
  pages     = {36},
  year      = {2017},
  abstract  = {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⊙.},
  language  = {en}
}
@article{ShenarRichardsonSablowskietal.2017,
  author    = {Shenar, Tomer and Richardson, N. D. and Sablowski, Daniel P. and Hainich, Rainer and Sana, H. and Moffat, A. F. J. and Todt, Helge Tobias and Hamann, Wolf-Rainer and Oskinova, Lida and Sander, Andreas Alexander Christoph and Tramper, Frank and Langer, Norbert and Bonanos, Alceste Z. and de Mink, Selma E. and Gr{\"a}fener, G. and Crowther, Paul and Vink, J. S. and Almeida, Leonardo A. and de Koter, A. and Barb{\´a}, Rodolfo and Herrero, A. and Ulaczyk, Krzysztof},
  title     = {The tarantula massive binary monitoring},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {598},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201629621},
  pages     = {16},
  year      = {2017},
  abstract  = {We present the first SB2 orbital solution and disentanglement of the massive Wolf-Rayet binary R145 (P = 159 d) located in the Large Magellanic Cloud. The primary was claimed to have a stellar mass greater than 300 M-circle dot, making it a candidate for being the most massive star known to date. While the primary is a known late-type, H-rich Wolf-Rayet star (WN6h), the secondary has so far not been unambiguously detected. Using moderate-resolution spectra, we are able to derive accurate radial velocities for both components. By performing simultaneous orbital and polarimetric analyses, we derive the complete set of orbital parameters, including the inclination. The spectra are disentangled and spectroscopically analyzed, and an analysis of the wind-wind collision zone is conducted. The disentangled spectra and our models are consistent with a WN6h type for the primary and suggest that the secondary is an O3.5 If*/WN7 type star. We derive a high eccentricity of e = 0 : 78 and minimum masses of M-1 sin(3) i approximate to M-2 sin(3) i = 13 +/- 2 M-circle dot, with q = M-2/M-1 = 1.01 +/- 0.07. An analysis of emission excess stemming from a wind-wind collision yields an inclination similar to that obtained from polarimetry (i = 39 +/- 6 degrees). Our analysis thus implies M-1 = 53(-20)(+40) and M2 = 54(-20)(+40) M-circle dot, excluding M-1 > 300 M-circle dot. A detailed comparison with evolution tracks calculated for single and binary stars together with the high eccentricity suggests that the components of the system underwent quasi-homogeneous evolution and avoided mass-transfer. This scenario would suggest current masses of approximate to 80 M-circle dot and initial masses of M-i,M-1 approximate to 10(5) and M-i,M-2 approximate to 90 M-circle dot, consistent with the upper limits of our derived orbital masses, and would imply an age of approximate to 2.2 Myr.},
  language  = {en}
}