Showing 1–50 of 63 results for author: Hénault-Brunet, V

Probing populations of dark stellar remnants in the globular clusters 47 Tuc and Terzan 5 using pulsar timing

Authors: Peter J. Smith, Vincent Hénault-Brunet, Nolan Dickson, Mark Gieles, Holger Baumgardt

Abstract: We present a new method to combine multimass equilibrium dynamical models and pulsar timing data to constrain the mass distribution and remnant populations of Milky Way globular clusters (GCs). We first apply this method to 47 Tuc, a cluster for which there exists an abundance of stellar kinematic data and which is also host to a large population of millisecond pulsars. We demonstrate that the pul… ▽ More We present a new method to combine multimass equilibrium dynamical models and pulsar timing data to constrain the mass distribution and remnant populations of Milky Way globular clusters (GCs). We first apply this method to 47 Tuc, a cluster for which there exists an abundance of stellar kinematic data and which is also host to a large population of millisecond pulsars. We demonstrate that the pulsar timing data allow us to place strong constraints on the overall mass distribution and remnant populations even without fitting on stellar kinematics. Our models favor a small population of stellar-mass BHs in this cluster (with a total mass of $412^{+77}_{-146} \mathrm{M_\odot}$), arguing against the need for a large ($ > 2000 \ \mathrm{M_\odot}$) central intermediate-mass black hole. We then apply the method to Terzan 5, a heavily obscured bulge cluster which hosts the largest population of millisecond pulsars of any Milky Way GC and for which the collection of conventional stellar kinematic data is very limited. We improve existing constraints on the mass distribution and structural parameters of this cluster and place stringent constraints on its black hole content, finding an upper limit on the mass in BHs of $\sim 4500 \ \mathrm{M_\odot}$. This method allows us to probe the central dynamics of GCs even in the absence of stellar kinematic data and can be easily applied to other GCs with pulsar timing data, for which datasets will continue to grow with the next generation of radio telescopes. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 32 pages, 19 figures. Submitted to ApJ, comments welcome

FORECASTOR -- I. Finding Optics Requirements and Exposure times for the Cosmological Advanced Survey Telescope for Optical and UV Research mission

Authors: Isaac Cheng, Tyrone E. Woods, Patrick Côté, Jennifer Glover, Dhananjhay Bansal, Melissa Amenouche, Madeline A. Marshall, Laurie Amen, John Hutchings, Laura Ferrarese, Kim A. Venn, Michael Balogh, Simon Blouin, Ryan Cloutier, Nolan Dickson, Sarah Gallagher, Martin Hellmich, Vincent Hénault-Brunet, Viraja Khatu, Cameron Lawlor-Forsyth, Cameron Morgan, Harvey Richer, Marcin Sawicki, Robert Sorba

Abstract: The Cosmological Advanced Survey Telescope for Optical and ultraviolet Research (CASTOR) is a proposed Canadian-led 1m-class space telescope that will carry out ultraviolet and blue-optical wide-field imaging, spectroscopy, and photometry. CASTOR will provide an essential bridge in the post-Hubble era, preventing a protracted UV-optical gap in space astronomy and enabling an enormous range of disc… ▽ More The Cosmological Advanced Survey Telescope for Optical and ultraviolet Research (CASTOR) is a proposed Canadian-led 1m-class space telescope that will carry out ultraviolet and blue-optical wide-field imaging, spectroscopy, and photometry. CASTOR will provide an essential bridge in the post-Hubble era, preventing a protracted UV-optical gap in space astronomy and enabling an enormous range of discovery opportunities from the solar system to the nature of the Cosmos, in conjunction with the other great wide-field observatories of the next decade (e.g., Euclid, Roman, Vera Rubin). FORECASTOR (Finding Optics Requirements and Exposure times for CASTOR) will supply a coordinated suite of mission-planning tools that will serve as the one-stop shop for proposal preparation, data reduction, and analysis for the CASTOR mission. We present the first of these tools: a pixel-based, user-friendly, extensible, multi-mission exposure time calculator (ETC) built in Python, including a modern browser-based graphical user interface that updates in real time. We then provide several illustrative examples of FORECASTOR's use that advance the design of planned legacy surveys for the CASTOR mission: a search for the most massive white dwarfs in the Magellanic Clouds; a study of the frequency of flaring activity in M stars, their distribution and impacts on habitability of exoplanets; mapping the proper motions of faint stars in the Milky Way; wide and deep galaxy surveys; and time-domain studies of active galactic nuclei. △ Less

Submitted 30 March, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: Updated references and acknowledgements to match published version. 24 pages, 16 figures, 3 tables, published in AJ

A closer look at the binary content of NGC 1850

Authors: S. Saracino, S. Kamann, N. Bastian, M. Gieles, T. Shenar, N. Reindl, J. Müller-Horn, C. Usher, S. Dreizler, V. Hénault-Brunet

Abstract: Studies of young clusters have shown that a large fraction of O-/early B-type stars are in binary systems, where the binary fraction increases with mass. These massive stars are present in clusters of a few Myrs, but gradually disappear for older clusters. The lack of detailed studies of intermediate-age clusters has meant that almost no information is available on the multiplicity properties of s… ▽ More Studies of young clusters have shown that a large fraction of O-/early B-type stars are in binary systems, where the binary fraction increases with mass. These massive stars are present in clusters of a few Myrs, but gradually disappear for older clusters. The lack of detailed studies of intermediate-age clusters has meant that almost no information is available on the multiplicity properties of stars with M $<$ 4 $M_{\odot}$. In this study we present the first characterization of the binary content of NGC 1850, a 100 Myr-old massive star cluster in the Large Magellanic Cloud, relying on a VLT/MUSE multi-epoch spectroscopic campaign. By sampling stars down to M = 2.5 $M_{\odot}$, we derive a close binary fraction of 24 $\pm$ 5 \% in NGC 1850, in good agreement with the multiplicity frequency predicted for stars of this mass range. We also find a trend with stellar mass (magnitude), with higher mass (brighter) stars having higher binary fractions. We modeled the radial velocity curves of individual binaries using The Joker and constrained the orbital properties of 27 systems, $\sim$17\% of all binaries with reliable radial velocities in NGC 1850. This study has brought to light a number of interesting objects, such as four binaries showing mass functions f(M) $>$ 1.25 $M_{\odot}$. One of these, star #47, has a peculiar spectrum, explainable with the presence of two disks in the system, around the visible star and the dark companion, which is a black hole candidate. These results confirm the importance and urgency of studying the binary content of clusters of any age. △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: 24 pages, 22 Figures, 1 Table; Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

The young massive SMC cluster NGC 330 seen by MUSE III. Stellar parameters and rotational velocities

Authors: J. Bodensteiner, H. Sana, P. L. Dufton, C. Wang, N. Langer, G. Banyard, L. Mahy, A. de Koter, S. E. de Mink, C. J. Evans, Y. G��tberg, V. Hénault-Brunet, L. R. Patrick, F. R. N. Schneider

Abstract: The origin of initial rotation rates of stars, and how a star's surface rotational velocity changes during the evolution, either by internal angular momentum transport or due to interactions with a binary companion, remain open questions in stellar astrophysics. Here, we aim to derive the physical parameters and study the distribution of (projected) rotational velocities of B-type stars in the 35… ▽ More The origin of initial rotation rates of stars, and how a star's surface rotational velocity changes during the evolution, either by internal angular momentum transport or due to interactions with a binary companion, remain open questions in stellar astrophysics. Here, we aim to derive the physical parameters and study the distribution of (projected) rotational velocities of B-type stars in the 35 Myr-old, massive cluster NGC 330 in the Small Magellanic Cloud. NGC 330 is in an age range where the number of post-interaction binaries is predicted to be high near the cluster turnoff (TO). We develop a simultaneous photometric and spectroscopic grid-fitting method adjusting atmosphere models on multi-band Hubble Space Telescope photometry and Multi Unit Spectroscopic Explorer spectroscopy. This allows us to homogeneously constrain the physical parameters of over 250 B and Be stars, brighter than mF814W = 18.8 mag. The rotational velocities of Be stars in NGC 330 are significantly higher than the ones of B stars. The rotational velocities vary as a function of the star's position in the color-magnitude diagram, qualitatively following predictions of binary population synthesis. A comparison to younger clusters shows that stars in NGC 330 rotate more rapidly on average. The rotational velocities of the 35 Myr old population in NGC 330 quantitatively agree with predictions for a stellar population that underwent significant binary interactions: the bulk of the B stars could be single stars or primaries in pre-interaction binaries. The rapidly spinning Be stars could be mass and angular momentum gainers in previous interactions, while those Be stars close to the TO may be spun-up single stars. The slowly rotating, apparently single stars above the TO could be merger products. The different vsini-characteristics of NGC 330 compared to younger populations can be understood in this framework. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 18 pages (incl. appendix), 15 figures, 3 tables, accepted for publication in A&A

Journal ref: A&A 680, A32 (2023)

Multimass modelling of Milky Way globular clusters -- II. present-day black hole populations

Authors: Nolan Dickson, Peter J. Smith, Vincent Hénault-Brunet, Mark Gieles, Holger Baumgardt

Abstract: Populations of stellar-mass black holes (BHs) in globular clusters (GCs) influence their dynamical evolution and have important implications on one of the main formation channels for gravitational wave sources. Inferring the size of these populations remains difficult, however. In this work, multimass models of 34 Milky Way GCs, first presented in Dickson et al., are used to explore the present-da… ▽ More Populations of stellar-mass black holes (BHs) in globular clusters (GCs) influence their dynamical evolution and have important implications on one of the main formation channels for gravitational wave sources. Inferring the size of these populations remains difficult, however. In this work, multimass models of 34 Milky Way GCs, first presented in Dickson et al., are used to explore the present-day BH populations. Direct constraints on both the total and visible mass components provided by several observables allow these models to accurately determine the distribution of the dark mass (including BHs) within clusters, as we demonstrate in a proof-of-concept fitting of the models to mock observations extracted from Monte Carlo cluster models. New constraints on the BH population retained to the present-day in each cluster are inferred from our models. We find that BH mass fractions ranging from 0 to 1 per cent of the total mass are typically required to explain the observations, except for Omega Cen, for which we infer a mass fraction above 5 per cent, in agreement with previous works. Relationships between the dark remnant populations and other cluster parameters are examined, demonstrating a clear anti-correlation between the amount of BHs and mass segregation between visible stars, as well as a correlation between remnant mass fractions and the dynamical age of clusters. Our inferred BH populations are in good agreement overall with other recent studies using different methodologies, but with notable discrepancies for individual clusters. △ Less

Submitted 12 February, 2024; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: 18 pages, 9 figures, MNRAS accepted

Multimass modelling of Milky Way globular clusters -- I. Implications on their stellar initial mass function above 1 M$_{\odot}$

Authors: Nolan Dickson, Vincent Hénault-Brunet, Holger Baumgardt, Mark Gieles, Peter Smith

Abstract: The distribution of stars and stellar remnants (white dwarfs, neutron stars, black holes) within globular clusters holds clues about their formation and long-term evolution, with important implications for their initial mass function (IMF) and the formation of black hole mergers. In this work, we present best-fitting multimass models for 37 Milky Way globular clusters, which were inferred from var… ▽ More The distribution of stars and stellar remnants (white dwarfs, neutron stars, black holes) within globular clusters holds clues about their formation and long-term evolution, with important implications for their initial mass function (IMF) and the formation of black hole mergers. In this work, we present best-fitting multimass models for 37 Milky Way globular clusters, which were inferred from various datasets, including proper motions from Gaia EDR3 and HST, line-of-sight velocities from ground-based spectroscopy and deep stellar mass functions from HST. We use metallicity dependent stellar evolution recipes to obtain present-day mass functions of stars and remnants from the IMF. By dynamically probing the present-day mass function of all objects in a cluster, including the mass distribution of remnants, these models allow us to explore in detail the stellar (initial) mass functions of a large sample of Milky Way GCs. We show that, while the low-mass mass function slopes are strongly dependent on the dynamical age of the clusters, the high-mass slope ($α_3; m > 1 M_\odot$) is not, indicating that the mass function in this regime has generally been less affected by dynamical mass loss. Examination of this high-mass mass function slope suggests an IMF in this mass regime consistent with a Salpeter IMF is required to reproduce the observations. This high-mass IMF is incompatible with a top-heavy IMF, as has been proposed recently. Finally, based on multimass model fits to our sample of Milky Way GCs, no significant correlation is found between the high-mass IMF slope and cluster metallicity. △ Less

Submitted 24 April, 2023; v1 submitted 2 March, 2023; originally announced March 2023.

Comments: 21 pages, 10 figures, MNRAS accepted

Evidence for a bottom-light initial mass function in massive star clusters

Authors: Holger Baumgardt, Vincent Henault-Brunet, Nolan Dickson, Antonio Sollima

Abstract: We have determined stellar mass functions of 120 Milky Way globular clusters and massive LMC/SMC star clusters based on a comparison of archival Hubble Space Telescope photometry with a large grid of direct N-body simulations. We find a strong correlation of the global mass function slopes of star clusters with both their internal relaxation times as well as their lifetimes. Once dynamical effects… ▽ More We have determined stellar mass functions of 120 Milky Way globular clusters and massive LMC/SMC star clusters based on a comparison of archival Hubble Space Telescope photometry with a large grid of direct N-body simulations. We find a strong correlation of the global mass function slopes of star clusters with both their internal relaxation times as well as their lifetimes. Once dynamical effects are being accounted for, the mass functions of most star clusters are compatible with an initial mass function described by a broken power-law distribution $N(m) \sim m^α$ with break masses at 0.4 M$_\odot$ and 1.0 M$_\odot$ and mass function slopes of $α_{Low}=-0.3$ for stars with masses $m<0.4$ M$_\odot$, $α_{High}=-2.30$ for stars with $m>1.0$ M$_\odot$ and $α_{Med}=-1.65$ for intermediate-mass stars. Alternatively, a log-normal mass function with a characteristic mass $\log M_C = -0.36$ and width $σ_C=0.28$ for low-mass stars and a power-law mass function for stars with $m>1$ M$_\odot$ also fits our data. We do not find a significant environmental dependency of the initial mass function with either cluster mass, density, global velocity dispersion or metallicity. Our results lead to a larger fraction of high-mass stars in globular clusters compared to canonical Kroupa/Chabrier mass functions, increasing the efficiency of self-enrichment in clusters and helping to alleviate the mass budget problem of multiple stellar populations in globular clusters. By comparing our results with direct N-body simulations we finally find that only simulations in which most black holes are ejected by natal birth kicks correctly reproduce the observed correlations. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: 20 pages, 9 figures, accepted for publication in MNRAS, data is available from http://people.smp.uq.edu.au/HolgerBaumgardt/globular/index.html

Dynamics in the outskirts of four Milky Way globular clusters: it's the tides that dominate

Authors: Zhen Wan, Anthony D. Arnold, William H. Oliver, Geraint F. Lewis, Holger Baumgardt, Mark Gieles, Vincent Hénault-Brunet, Thomas de Boer, Eduardo Balbinot, Gary Da Costa, Dougal Mackey, Denis Erkal, Annette Ferguson, Pete Kuzma, Elena Pancino, Jorge Penarrubia, Nicoletta Sanna, Antonio Sollima, Roeland P. van der Marel, Laura L. Watkins

Abstract: We present the results of a spectroscopic survey of the outskirts of 4 globular clusters -- NGC 1261, NGC 4590, NGC 1904, and NGC 1851 -- covering targets within 1 degree from the cluster centres, with 2dF/AAOmega on the Anglo-Australian Telescope (AAT) and FLAMES on the Very Large Telescope (VLT). We extracted chemo-dynamical information for individual stars, from which we estimated the veloc… ▽ More We present the results of a spectroscopic survey of the outskirts of 4 globular clusters -- NGC 1261, NGC 4590, NGC 1904, and NGC 1851 -- covering targets within 1 degree from the cluster centres, with 2dF/AAOmega on the Anglo-Australian Telescope (AAT) and FLAMES on the Very Large Telescope (VLT). We extracted chemo-dynamical information for individual stars, from which we estimated the velocity dispersion profile and the rotation of each cluster. The observations are compared to direct $N$-body simulations and appropriate {\sc limepy}/{\sc spes} models for each cluster to interpret the results. In NGC 1851, the detected internal rotation agrees with existing literature, and NGC 1261 shows some rotation signal beyond the truncation radius, likely coming from the escaped stars. We find that the dispersion profiles for both the observations and the simulations for NGC 1261, NGC 1851, and NGC 1904 do not decrease as the {\sc limepy}/{\sc spes} models predict beyond the truncation radius, where the $N$-body simulations show that escaped stars dominate; the dispersion profile of NGC 4590 follows the predictions of the {\sc limepy}/{\sc spes} models, though the data do not effectively extend beyond the truncation radius. The increasing/flat dispersion profiles in the outskirts of NGC 1261, NGC 1851 and NGC 1904, are reproduced by the simulations. Hence, the increasing/flat dispersion profiles of the clusters in question can be explained by the tidal interaction with the Galaxy without introducing dark matter. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: 17 pages, 10 figures. Accepted by MNRAS. Comments are welcome

The VLT-FLAMES Tarantula Survey: Observational evidence for two distinct populations of massive runaway stars in 30 Doradus

Authors: H. Sana, O. H. Ramírez-Agudelo, V. Hénault-Brunet, L. Mahy, L. A. Almeida, A. de Koter, J. M. Bestenlehner, C. J. Evans, N. Langer, F. R. N. Schneider, P. A. Crowther, S. E. de Mink, A. Herrero, D. J. Lennon, M. Gieles, J. Maíz Apellániz, M. Renzo, E. Sabbi, J. Th. van Loon, J. S. Vink

Abstract: Two main scenarios have been proposed for origin of massive runaway stars -- dynamical ejection or release from a binary at the first core collapse -- but their relative contribution remains debated. Using two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active… ▽ More Two main scenarios have been proposed for origin of massive runaway stars -- dynamical ejection or release from a binary at the first core collapse -- but their relative contribution remains debated. Using two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active star-forming region in the Local group. We use RV measurements of the O-type star populations in 30 Doradus obtained by the VLT-FLAMES Tarantula Survey and the Tarantula Massive Binary Monitoring to identify single and binary O-type runaways. We discuss their rotational properties and qualitatively compare observations with expectations of ejection scenarios. We identify 23 single and one binary O-type runaway objects, most of them outside the main star-forming regions in 30 Doradus. We find an overabundance of rapid rotators (vsini > 200km/s) among the runaway population, providing an explanation of the overabundance of rapidly rotating stars in the 30 Doradus field. Considerations of the projected rotation rates and runaway line-of-sight (los) velocities reveal a conspicuous absence of rapidly rotating (vsini > 210k/ms), fast moving (v_{los} > 60km/s) runaways, and suggest the presence of two different populations of runaway stars: a population of rapidly-spinning but slowly moving runaways and a population of fast moving but slowly rotating ones. These are detected with a ratio close to 2:1 in our sample. We argue that slowly moving but rapidly spinning runaways result from binary ejections, while rapidly moving but slowly spinning runaways could result from dynamical ejections. Given that detection biases will more strongly impact the slow-moving population, our results suggest that the binary evolution scenario dominates the current massive runaway population in 30 Doradus. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: Accepted for publication in A&A Letters; 9 pages, 5 figures

Journal ref: A&A 668, L5 (2022)

An X-ray quiet black hole born with a negligible kick in a massive binary within the Large Magellanic Cloud

Authors: Tomer Shenar, Hugues Sana, Laurent Mahy, Kareem El-Badry, Pablo Marchant, Norbert Langer, Calum Hawcroft, Matthias Fabry, Koushik Sen, Leonardo A. Almeida, Michael Abdul-Masih, Julia Bodensteiner, Paul A. Crowther, Mark Gieles, Mariusz Gromadzki, Vincent Henault-Brunet, Artemio Herrero, Alex de Koter, Patryk Iwanek, Szymon Kozłowski, Daniel J. Lennon, Jesus Maız Apellaniz, Przemysław Mroz, Anthony F. J. Moffat, Annachiara Picco , et al. (13 additional authors not shown)

Abstract: Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-ty… ▽ More Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Comments: Accepted to Nature Astronomy, 64 pages, 15 figures, 2 tables; ESO press release: https://www.eso.org/public/news/eso2210/; Nat Asr paper URL: https://www.nature.com/articles/s41550-022-01730-y

A black hole detected in the young massive LMC cluster NGC 1850

Authors: S. Saracino, S. Kamann, M. G. Guarcello, C. Usher, N. Bastian, I. Cabrera-Ziri, M. Gieles, S. Dreizler, G. S. Da Costa, T. -O. Husser, V. Hénault-Brunet

Abstract: We report the detection of a black hole (NGC 1850 BH1) in the $\sim$100 Myr-old massive cluster NGC~1850 in the Large Magellanic Cloud. It is in a binary system with a main-sequence turn-off star (4.9 $\pm$ 0.4 M${_\odot}$), which is starting to fill its Roche Lobe and becoming distorted. Using 17 epochs of VLT/MUSE observations we detected radial velocity variations exceeding 300 km/s associated… ▽ More We report the detection of a black hole (NGC 1850 BH1) in the $\sim$100 Myr-old massive cluster NGC~1850 in the Large Magellanic Cloud. It is in a binary system with a main-sequence turn-off star (4.9 $\pm$ 0.4 M${_\odot}$), which is starting to fill its Roche Lobe and becoming distorted. Using 17 epochs of VLT/MUSE observations we detected radial velocity variations exceeding 300 km/s associated to the target star, linked to the ellipsoidal variations measured by OGLE-IV in the optical bands. Under the assumption of a semi-detached system, the simultaneous modelling of radial velocity and light curves constraints the orbital inclination of the binary to ($38 \pm 2$)$^{\circ}$, resulting in a true mass of the unseen companion of $11.1_{-2.4}^{+2.1}$ $M_{\odot}$. This represents the first direct dynamical detection of a black hole in a young massive cluster, opening up the possibility of studying the initial mass function and the early dynamical evolution of such compact objects in high-density environments. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 11 pages, 8 figures, 4 Tables. Accepted for publication by Monthly Notices of the Royal Astronomical Society

MUSE narrow field mode observations of the central kinematics of M15

Authors: Christopher Usher, Sebastian Kamann, Mark Gieles, Vincent Hénault-Brunet, Emanuele Dalessandro, Eduardo Balbinot, Antonio Sollima

Abstract: We present observations of the stellar kinematics of the centre of the core collapsed globular cluster M15 obtained with the MUSE integral field spectrograph on the VLT operating in narrow field mode. Thanks to the use of adaptive optics, we obtain a spatial resolution of 0.1arcsec and are able to reliably measure the radial velocities of 864 stars within 8 arcsec of the centre of M15 thus providi… ▽ More We present observations of the stellar kinematics of the centre of the core collapsed globular cluster M15 obtained with the MUSE integral field spectrograph on the VLT operating in narrow field mode. Thanks to the use of adaptive optics, we obtain a spatial resolution of 0.1arcsec and are able to reliably measure the radial velocities of 864 stars within 8 arcsec of the centre of M15 thus providing the largest sample of radial velocities ever obtained for the innermost regions of this system. Combined with previous observations of M15 using MUSE in wide field mode and literature data, we find that the central kinematics of M15 are complex with the rotation axis of the core of M15 offset from the rotation axis of the bulk of the cluster. While this complexity has been suggested by previous work, we confirm it at higher significance and in more detail. △ Less

Submitted 24 February, 2021; v1 submitted 23 February, 2021; originally announced February 2021.

Comments: 8 pages, 6 figures. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

The dynamics of the globular cluster NGC3201 out to the Jacobi radius

Authors: Zhen Wan, William Oliver, Holger Baumgardt, Geraint Lewis, Mark Gieles, Vincent Hénault-Brunet, Thomas de Boer, Eduardo Balbinot, Gary Da Costa, Dougal Mackey

Abstract: As part of a chemo-dynamical survey of five nearby globular clusters with 2dF/AAOmega on the AAT, we have obtained kinematic information for the globular cluster NGC3201. Our new observations confirm the presence of a significant velocity gradient across the cluster which can almost entirely be explained by the high proper motion of the cluster. After subtracting the contribution of this perspecti… ▽ More As part of a chemo-dynamical survey of five nearby globular clusters with 2dF/AAOmega on the AAT, we have obtained kinematic information for the globular cluster NGC3201. Our new observations confirm the presence of a significant velocity gradient across the cluster which can almost entirely be explained by the high proper motion of the cluster. After subtracting the contribution of this perspective rotation, we found a remaining rotation signal with an amplitude of $\sim1\ km/s$ around a different axis to what we expect from the tidal tails and the potential escapers, suggesting that this rotation is internal and can be a remnant of its formation process. At the outer part, we found a rotational signal that is likely a result from potential escapers. The proper motion dispersion at large radii reported by Bianchini et al. has previously been attributed to dark matter. Here we show that the LOS dispersion between 0.5-1 Jacobi radius is lower, yet above the predictions from an N-body model of NGC3201 that we ran for this study. Based on the simulation, we find that potential escapers cannot fully explain the observed velocity dispersion. We also estimate the effect on the velocity dispersion of different amounts of stellar-mass black holes and unbound stars from the tidal tails with varying escape rates and find that these effects cannot explain the difference between the LOS dispersion and the N-body model. Given the recent discovery of tidal tail stars at large distances from the cluster, a dark matter halo is an unlikely explanation. We show that the effect of binary stars, which is not included in the N-body model, is important and can explain part of the difference in dispersion. We speculate that the remaining difference must be the result of effects not included in the N-body model, such as initial cluster rotation, velocity anisotropy and Galactic substructure. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: 14 pages, 13 figures, accepted by MNRAS

Astronomy in a Low-Carbon Future

Authors: Christopher D. Matzner, Nicolas B. Cowan, René Doyon, Vincent Hénault-Brunet, David Lafrenère, Martine Lokken, Peter G. Martin, Sharon Morsink, Magdalen Normandeau, Nathalie Ouellette, Mubdi Rahman, Joel Roediger, James Taylor, Rob Thacker, Marten van Kerkwijk

Abstract: The global climate crisis poses new risks to humanity, and with them, new challenges to the practices of professional astronomy. Avoiding the more catastrophic consequences of global warming by more than 1.5 degrees requires an immediate reduction of greenhouse gas emissions. According to the 2018 United Nations Intergovernmental Panel report, this will necessitate a 45% reduction of emissions by… ▽ More The global climate crisis poses new risks to humanity, and with them, new challenges to the practices of professional astronomy. Avoiding the more catastrophic consequences of global warming by more than 1.5 degrees requires an immediate reduction of greenhouse gas emissions. According to the 2018 United Nations Intergovernmental Panel report, this will necessitate a 45% reduction of emissions by 2030 and net-zero emissions by 2050. Efforts are required at all levels, from the individual to the governmental, and every discipline must find ways to achieve these goals. This will be especially difficult for astronomy with its significant reliance on conference and research travel, among other impacts. However, our long-range planning exercises provide the means to coordinate our response on a variety of levels. We have the opportunity to lead by example, rising to the challenge rather than reacting to external constraints. We explore how astronomy can meet the challenge of a changing climate in clear and responsible ways, such as how we set expectations (for ourselves, our institutions, and our granting agencies) around scientific travel, the organization of conferences, and the design of our infrastructure. We also emphasize our role as reliable communicators of scientific information on a problem that is both human and planetary in scale. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: 16 pages. A White Paper prepared for the Canadian Long Range Plan 2020

On the black hole content and initial mass function of 47 Tuc

Authors: Vincent Hénault-Brunet, Mark Gieles, Jay Strader, Miklos Peuten, Eduardo Balbinot, Kaela E. K. Douglas

Abstract: The globular cluster (GC) 47 Tuc has recently been proposed to host an intermediate-mass black hole (IMBH) or a population of stellar-mass black holes (BHs). To shed light on its dark content, we present an application of self-consistent multimass models with a varying mass function and content of stellar remnants, which we fit to various observational constraints. Our best-fitting model successfu… ▽ More The globular cluster (GC) 47 Tuc has recently been proposed to host an intermediate-mass black hole (IMBH) or a population of stellar-mass black holes (BHs). To shed light on its dark content, we present an application of self-consistent multimass models with a varying mass function and content of stellar remnants, which we fit to various observational constraints. Our best-fitting model successfully matches the observables and correctly predicts the radial distribution of millisecond pulsars and their gravitational accelerations inferred from long-term timing observations. The data favours a population of BHs with a total mass of $430^{+386}_{-301}$ $M_{\odot}$, but the most likely model has very few BHs. Since our models do not include a central IMBH and accurately reproduce the observations, we conclude that there is currently no need to invoke the presence of an IMBH in 47 Tuc. The global present-day mass function inferred is significantly depleted in low-mass stars (power-law slope $α=-0.52^{+0.17}_{-0.16}$). Given the orbit and predicted mass-loss history of this massive GC, the dearth of low-mass stars is difficult to explain with a standard initial mass function (IMF) followed by long-term preferential escape of low-mass stars driven by two-body relaxation, and instead suggests that 47 Tuc may have formed with a bottom-light IMF. We discuss alternative evolutionary origins for the flat mass function and ways to reconcile this with the low BH retention fraction. Finally, by capturing the effect of dark remnants, our method offers a new way to probe the IMF in a GC above the current main-sequence turn-off mass, for which we find a slope of $-2.49\pm0.08$. △ Less

Submitted 25 October, 2019; v1 submitted 22 August, 2019; originally announced August 2019.

Comments: 16 pages, 7 figures, accepted to MNRAS after minor revision

The Detailed Science Case for the Maunakea Spectroscopic Explorer, 2019 edition

Authors: The MSE Science Team, Carine Babusiaux, Maria Bergemann, Adam Burgasser, Sara Ellison, Daryl Haggard, Daniel Huber, Manoj Kaplinghat, Ting Li, Jennifer Marshall, Sarah Martell, Alan McConnachie, Will Percival, Aaron Robotham, Yue Shen, Sivarani Thirupathi, Kim-Vy Tran, Christophe Yeche, David Yong, Vardan Adibekyan, Victor Silva Aguirre, George Angelou, Martin Asplund, Michael Balogh, Projjwal Banerjee , et al. (239 additional authors not shown)

Abstract: (Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the sc… ▽ More (Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more. △ Less

Submitted 9 April, 2019; originally announced April 2019.

Comments: 9 chapters, 301 pages, 100 figures. This version of the DSC is a comprehensive update of the original version, released in 2016, which can be downloaded at arXiv:1606.00043. A detailed summary of the design of MSE is available in the MSE Book 2018, available at arXiv:1810.08695

Globular cluster number density profiles using Gaia DR2

Authors: T. J. L. de Boer, M. Gieles, E. Balbinot, V. Henault-Brunet, A. Sollima, L. L. Watkins, I. Claydon

Abstract: Using data from Gaia DR2, we study the radial number density profiles of the Galactic globular cluster sample. Proper motions are used for accurate membership selection, especially crucial in the cluster outskirts. Due to the severe crowding in the centres, the Gaia data is supplemented by literature data from HST and surface brightness measurements, where available. This results in 81 clusters wi… ▽ More Using data from Gaia DR2, we study the radial number density profiles of the Galactic globular cluster sample. Proper motions are used for accurate membership selection, especially crucial in the cluster outskirts. Due to the severe crowding in the centres, the Gaia data is supplemented by literature data from HST and surface brightness measurements, where available. This results in 81 clusters with a complete density profile covering the full tidal radius (and beyond) for each cluster. We model the density profiles using a set of single-mass models ranging from King and Wilson models to generalised lowered isothermal limepy models and the recently introduced spes models, which allow for the inclusion of potential escapers. We find that both King and Wilson models are too simple to fully reproduce the density profiles, with King (Wilson) models on average underestimating(overestimating) the radial extent of the clusters. The truncation radii derived from the limepy models are similar to estimates for the Jacobi radii based on the cluster masses and their orbits. We show clear correlations between structural and environmental parameters, as a function of Galactocentric radius and integrated luminosity. Notably, the recovered fraction of potential escapers correlates with cluster pericentre radius, luminosity and cluster concentration. The ratio of half mass over Jacobi radius also correlates with both truncation parameter and PE fraction, showing the effect of Roche lobe filling. △ Less

Submitted 5 March, 2019; v1 submitted 23 January, 2019; originally announced January 2019.

Comments: 23 pages, 23 figures. Accepted in MNRAS

Linking the rotation of a cluster to the spins of its stars: The kinematics of NGC6791 and NGC6819 in 3D

Authors: Sebastian Kamann, Nate Bastian, Mark Gieles, Eduardo Balbinot, Vincent Hénault-Brunet

Abstract: The physics governing the formation of star clusters is still not entirely understood. One open question concerns the amount of angular momentum that newly formed clusters possess after emerging from their parent gas clouds. Recent results suggest an alignment of stellar spins and binary orbital spins in star clusters, which support a scenario in which clusters are born with net angular momentum c… ▽ More The physics governing the formation of star clusters is still not entirely understood. One open question concerns the amount of angular momentum that newly formed clusters possess after emerging from their parent gas clouds. Recent results suggest an alignment of stellar spins and binary orbital spins in star clusters, which support a scenario in which clusters are born with net angular momentum cascading down to stellar scales. In this paper, we combine Gaia data and published line of sight velocities to explore if NGC6791 and NGC6819, two of the clusters for which an alignment of stellar spins has been reported, rotate in the same plane as their stars. We find evidence for rotation in NGC6791 using both proper motions and line of sight velocities. Our estimate of the inclination angle is broadly consistent with the mean inclination that has been determined for its stars, but the uncertainties are still substantial. Our results identify NGC6791 as a promising follow-up candidate to investigate the link between cluster and stellar rotation. We find no evidence for rotation in NGC6819. △ Less

Submitted 20 November, 2018; originally announced November 2018.

Comments: 11 pages, 9 figures, accepted for publication in MNRAS

Gaia and HST astrometry of the very massive $\sim$150 $M_\odot$ candidate runaway star VFTS682

Authors: M. Renzo, S. E. de Mink, D. J. Lennon, I. Platais, R. P. van der Marel, E. Laplace, J. M. Bestenlehner, C. J. Evans, V. Hénault-Brunet, S. Justham, A. de Koter, N. Langer, F. Najarro, F. R. N. Schneider, J. S. Vink

Abstract: How very massive stars form is still an open question in astrophysics. VFTS682 is among the most massive stars known, with an inferred initial mass of $\sim$150 $M_\odot$ . It is located in 30 Doradus at a projected distance of 29 pc from the central cluster R136. Its apparent isolation led to two hypotheses: either it formed in relative isolation or it was ejected dynamically from the cluster. We… ▽ More How very massive stars form is still an open question in astrophysics. VFTS682 is among the most massive stars known, with an inferred initial mass of $\sim$150 $M_\odot$ . It is located in 30 Doradus at a projected distance of 29 pc from the central cluster R136. Its apparent isolation led to two hypotheses: either it formed in relative isolation or it was ejected dynamically from the cluster. We investigate the kinematics of VFTS682 as obtained by Gaia and Hubble Space Telescope astrometry. We derive a projected velocity relative to the cluster of $38 \pm 17 \mathrm{km \ s^{-1}}$ (1$σ$ confidence interval). Although the error bars are substantial, two independent measures suggest that VFTS682 is a runaway ejected from the central cluster. This hypothesis is further supported by a variety of circumstantial clues. The central cluster is known to harbor other stars more massive than 150 $M_\odot$ of similar spectral type and recent astrometric studies on VFTS16 and VFTS72 provide direct evidence that the cluster can eject some of its most massive members, in agreement with theoretical predictions. If future data confirm the runaway nature, this would make VFTS682 the most massive runaway star known to date. △ Less

Submitted 12 October, 2018; originally announced October 2018.

Comments: accepted for publication in MNRAS Letters, author's version with larger figures

Mass modelling globular clusters in the Gaia era: a method comparison using mock data from an $N$-body simulation of M4

Authors: Vincent Hénault-Brunet, Mark Gieles, Antonio Sollima, Laura L. Watkins, Alice Zocchi, Ian Claydon, Elena Pancino, Holger Baumgardt

Abstract: As we enter a golden age for studies of internal kinematics and dynamics of Galactic globular clusters (GCs), it is timely to assess the performance of modelling techniques in recovering the mass, mass profile, and other dynamical properties of GCs. Here, we compare different mass-modelling techniques (distribution-function (DF)-based models, Jeans models, and a grid of N-body models) by applying… ▽ More As we enter a golden age for studies of internal kinematics and dynamics of Galactic globular clusters (GCs), it is timely to assess the performance of modelling techniques in recovering the mass, mass profile, and other dynamical properties of GCs. Here, we compare different mass-modelling techniques (distribution-function (DF)-based models, Jeans models, and a grid of N-body models) by applying them to mock observations from a star-by-star N-body simulation of the GC M 4 by Heggie. The mocks mimic existing and anticipated data for GCs: surface brightness or number density profiles, local stellar mass functions, line-of-sight velocities, and Hubble Space Telescope- and Gaia-like proper motions. We discuss the successes and limitations of the methods. We find that multimass DF-based models, Jeans, and N-body models provide more accurate mass profiles compared to single-mass DF-based models. We highlight complications in fitting the kinematics in the outskirts due to energetically unbound stars associated with the cluster ("potential escapers", not captured by truncated DF models nor by N-body models of clusters in isolation), which can be avoided with DF-based models including potential escapers, or with Jeans models. We discuss ways to account for mass segregation. For example, three-component DF-based models with freedom in their mass function are a simple alternative to avoid the biases of single-mass models (which systematically underestimate the total mass, half-mass radius, and central density), while more realistic multimass DF-based models with freedom in the remnant content represent a promising avenue to infer the total mass and the mass function of remnants. △ Less

Submitted 22 November, 2018; v1 submitted 11 October, 2018; originally announced October 2018.

Comments: 27 pages, 12 figures, MNRAS accepted

A-type stars in the Canada-France Imaging Survey I. The stellar halo of the Milky Way traced to large radius by blue horizontal branch stars

Authors: Guillaume F. Thomas, Alan W. McConnachie, Rodrigo A. Ibata, Patrick Côté, Nicolas Martin, Else Starkenburg, Raymond Carlberg, Scott Chapman, Sébastien Fabbro, Benoit Famaey, Nicholas Fantin, Stephen Gwyn, Vincent Hénault-Brunet, Khyati Malhan, Julio Navarro, Annie C. Robin, Douglas Scott

Abstract: We present the stellar density profile of the outer halo of the Galaxy traced over a range of Galactocentric radii from $15< R_{GC} < 220$ kpc by blue horizontal branch (BHB) stars. These stars are identified photometrically using deep $u-$band imaging from the new Canada-France-Imaging-Survey (CFIS) that reaches 24.5 mag. This is combined with $griz$ bands from Pan-STARRS 1 and covers a total of… ▽ More We present the stellar density profile of the outer halo of the Galaxy traced over a range of Galactocentric radii from $15< R_{GC} < 220$ kpc by blue horizontal branch (BHB) stars. These stars are identified photometrically using deep $u-$band imaging from the new Canada-France-Imaging-Survey (CFIS) that reaches 24.5 mag. This is combined with $griz$ bands from Pan-STARRS 1 and covers a total of $\sim4000$ deg$^2$ of the northern sky. We present a new method to select BHB stars that has low contamination from blue stragglers and high completeness. We use this sample to measure and parameterize the three dimensional density profile of the outer stellar halo. We fit the profile using (i) a simple power-law with a constant flattening (ii) a flattening that varies as a function of Galactocentric radius (iii) a broken power law profile. We find that outer stellar halo traced by the BHB is well modelled by a broken power law with a constant flattening of $q=0.86 \pm 0.02$, with an inner slope of $γ=4.24 \pm 0.08$. This is much steeper than the preferred outer profile that has a slope of $β=3.21\pm 0.07$ after a break radius of $r_b=41.4^{+2.5}_{-2.4}$ kpc. The outer profile of the stellar halo trace by BHB stars is shallower than that recently measured using RR Lyrae, a surprising result given the broad similarity of the ages of these stellar populations. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Comments: 15 pages, 17 figures, accepted in MNRAS

Response to comment on "An excess of massive stars in the local 30 Doradus starburst"

Authors: F. R. N. Schneider, H. Sana, C. J. Evans, J. M. Bestenlehner, N. Castro, L. Fossati, G. Gräfener, N. Langer, O. H. Ramírez-Agudelo, C. Sabín-Sanjulián, S. Simón-Díaz, F. Tramper, P. A. Crowther, A. de Koter, S. E. de Mink, P. L. Dufton, M. Garcia, M. Gieles, V. Hénault-Brunet, A. Herrero, R. G. Izzard, V. Kalari, D. J. Lennon, J. Maíz Apellániz, N. Markova , et al. (7 additional authors not shown)

Abstract: Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradu… ▽ More Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradus. △ Less

Submitted 25 July, 2018; originally announced July 2018.

Comments: Authors' version of a response to a technical comment published in Science; 8 pages, 1 figure

The unexpected kinematics of multiple populations in NGC6362: do binaries play a role?

Authors: E. Dalessandro, A. Mucciarelli, M. Bellazzini, A. Sollima, E. Vesperini, J. Hong, V. Hénault-Brunet, F. R. Ferraro, R. Ibata, B. Lanzoni, D. Massari, M. Salaris

Abstract: We present a detailed analysis of the kinematic properties of the multiple populations (MPs) in the low-mass Galactic globular cluster NGC6362 based on a sample of about 500 member stars for which radial velocities (RVs), Fe and Na abundances have been homogeneously derived. At distances from the cluster center larger than about 0.5 $r_h$, we find that first (FG - Na-poor) and second generation (S… ▽ More We present a detailed analysis of the kinematic properties of the multiple populations (MPs) in the low-mass Galactic globular cluster NGC6362 based on a sample of about 500 member stars for which radial velocities (RVs), Fe and Na abundances have been homogeneously derived. At distances from the cluster center larger than about 0.5 $r_h$, we find that first (FG - Na-poor) and second generation (SG - Na-rich) stars show hints of different line-of-sight velocity dispersion profiles, with FG stars being dynamically hotter. This is the first time that differences in the velocity dispersion of MPs are detected by using only RVs. While kinematic differences between MPs in globular clusters are usually described in terms of anisotropy differences driven by the different radial distributions, this explanation seems hardly viable for NGC6362, where SG and FG stars are spatially mixed. We demonstrate that the observed difference in the velocity dispersion profiles can be accounted for by the effect of binary stars. In fact, thanks to our multi-epoch RV measurements we find that the binary fraction is significantly larger in the FG sample (f~14%) than in SG population (f<1%), and we show that such a difference can inflate the velocity dispersion of FG with respect to SG by the observed amount in the relevant radial range. Our results nicely match the predictions of state-of-the art $N$-body simulations of the co-evolution of MPs in globular clusters that include the effect of binaries. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Comments: 18 pages, 6 figures; Accepted for publication in ApJ

The effect of stellar-mass black holes on the central kinematics of omega Cen: a cautionary tale for IMBH interpretations

Authors: Alice Zocchi, Mark Gieles, Vincent Hénault-Brunet

Abstract: The search for intermediate-mass black holes (IMBHs) in the centre of globular clusters is often based on the observation of a central cusp in the surface brightness profile and a rise towards the centre in the velocity dispersion profiles. Similar signatures, however, could result from other effects, that need to be taken into account in order to determine the presence (or the absence) of an IMBH… ▽ More The search for intermediate-mass black holes (IMBHs) in the centre of globular clusters is often based on the observation of a central cusp in the surface brightness profile and a rise towards the centre in the velocity dispersion profiles. Similar signatures, however, could result from other effects, that need to be taken into account in order to determine the presence (or the absence) of an IMBH in these stellar systems. Following our previous exploration of the role of radial anisotropy in shaping these observational signatures, we analyse here the effects produced by the presence of a population of centrally concentrated stellar-mass black holes. We fit dynamical models to omega Cen data, and we show that models with ~5% of their mass in black holes (consistent with ~100% retention fraction after natal kicks) can reproduce the data. When simultaneously considering both radial anisotropy and mass segregation, the best-fit model includes a smaller population of remnants, and a less extreme degree of anisotropy with respect to the models that include only one of these features. These results underline that before conclusions about putative IMBHs can be made, the effects of stellar-mass black holes and radial anisotropy need to be properly accounted for. △ Less

Submitted 6 June, 2018; originally announced June 2018.

Comments: 13 pages, 5 figures. Accepted for publication in MNRAS

Concurrent formation of supermassive stars and globular clusters: implications for early self-enrichment

Authors: Mark Gieles, Corinne Charbonnel, Martin Krause, Vincent Henault-Brunet, Oscar Agertz, Henny Lamers, Nathan Bastian, Alessia Gualandris, Alice Zocchi, James Petts

Abstract: We present a model for the concurrent formation of globular clusters (GCs) and supermassive stars (SMSs, $>10^3\,{\rm M}_\odot$) to address the origin of the HeCNONaMgAl abundance anomalies in GCs. GCs form in converging gas flows and accumulate low-angular momentum gas, which accretes onto protostars. This leads to an adiabatic contraction of the cluster and an increase of the stellar collision r… ▽ More We present a model for the concurrent formation of globular clusters (GCs) and supermassive stars (SMSs, $>10^3\,{\rm M}_\odot$) to address the origin of the HeCNONaMgAl abundance anomalies in GCs. GCs form in converging gas flows and accumulate low-angular momentum gas, which accretes onto protostars. This leads to an adiabatic contraction of the cluster and an increase of the stellar collision rate. A SMS can form via runaway collisions if the cluster reaches sufficiently high density before two-body relaxation halts the contraction. This condition is met if the number of stars $\gtrsim10^6$ and the gas accretion rate $\gtrsim10^5\,{\rm M}_\odot$/Myr, reminiscent of GC formation in high gas-density environments, such as -- but not restricted to -- the early Universe. The strong SMS wind mixes with the inflowing pristine gas, such that the protostars accrete diluted hot-hydrogen burning yields of the SMS. Because of continuous rejuvenation, the amount of processed material liberated by the SMS can be an order of magnitude higher than its maximum mass. This `conveyor-belt' production of hot-hydrogen burning products provides a solution to the mass budget problem that plagues other scenarios. Additionally, the liberated material is mildly enriched in helium and relatively rich in other hot-hydrogen burning products, in agreement with abundances of GCs today. Finally, we find a super-linear scaling between the amount of processed material and cluster mass, providing an explanation for the observed increase of the fraction of processed material with GC mass. We discuss open questions of this new GC enrichment scenario and propose observational tests. △ Less

Submitted 24 April, 2018; v1 submitted 12 April, 2018; originally announced April 2018.

Comments: 19 pages, 5 figures, accepted to MNRAS (shortened abstract and included feedback from the community)

An excess of massive stars in the local 30 Doradus starburst

Authors: F. R. N. Schneider, H. Sana, C. J. Evans, J. M. Bestenlehner, N. Castro, L. Fossati, G. Gräfener, N. Langer, O. H. Ramírez-Agudelo, C. Sabín-Sanjulián, S. Simón-Díaz, F. Tramper, P. A. Crowther, A. de Koter, S. E. de Mink, P. L. Dufton, M. Garcia, M. Gieles, V. Hénault-Brunet, A. Herrero, R. G. Izzard, V. Kalari, D. J. Lennon, J. Maíz Apellániz, N. Markova , et al. (7 additional authors not shown)

Abstract: The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation s… ▽ More The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation started about $8\,\mathrm{Myr}$ ago and the star-formation rate seems to have declined in the last $1\,\mathrm{Myr}$. The IMF is densely sampled up to $200\,\mathrm{M}_\odot$ and contains $32\pm12\%$ more stars above $30\,\mathrm{M}_\odot$ than predicted by a standard Salpeter IMF. In the mass range $15-200\,\mathrm{M}_\odot$, the IMF power-law exponent is $1.90^{+0.37}_{-0.26}$, shallower than the Salpeter value of 2.35. △ Less

Submitted 9 January, 2018; originally announced January 2018.

Comments: Authors' version of paper published in Science at http://science.sciencemag.org/content/359/6371/69 ; 15 pages main text (3 figures), 47 pages supplementary materials (10 figures, 3 tables)

Journal ref: Science 359, 69 (2018)

Mass models of NGC 6624 without an intermediate-mass black hole

Authors: Mark Gieles, Eduardo Balbinot, Rashid Yaaqib, Vincent Henault-Brunet, Alice Zocchi, Miklos Peuten, Peter G. Jonker

Abstract: An intermediate-mass black hole (IMBH) was recently reported to reside in the centre of the Galactic globular cluster (GC) NGC 6624, based on timing observations of a millisecond pulsar (MSP) located near the cluster centre in projection. We present dynamical models with multiple mass components of NGC 6624 - without an IMBH - which successfully describe the surface brightness profile and proper m… ▽ More An intermediate-mass black hole (IMBH) was recently reported to reside in the centre of the Galactic globular cluster (GC) NGC 6624, based on timing observations of a millisecond pulsar (MSP) located near the cluster centre in projection. We present dynamical models with multiple mass components of NGC 6624 - without an IMBH - which successfully describe the surface brightness profile and proper motion kinematics from the Hubble Space Telescope (HST) and the stellar mass function at different distances from the cluster centre. The maximum line-of-sight acceleration at the position of the MSP accommodates the inferred acceleration of the MSP, as derived from its first period derivative. With discrete realizations of the models we show that the higher-order period derivatives - which were previously used to derive the IMBH mass - are due to passing stars and stellar remnants, as previously shown analytically in literature. We conclude that there is no need for an IMBH to explain the timing observations of this MSP. △ Less

Submitted 29 November, 2017; v1 submitted 20 September, 2017; originally announced September 2017.

Comments: 8 pages, 7 figures, MNRAS. Updated to match final journal style

Journal ref: Mon. Not. R. Astron. Soc. 473, 4832-4839 (2018)

Testing lowered isothermal models with direct N-body simulations of globular clusters - II: Multimass models

Authors: Miklos Peuten, Alice Zocchi, Mark Gieles, Vincent Hénault-Brunet

Abstract: Lowered isothermal models, such as the multimass Michie-King models, have been successful in describing observational data of globular clusters. In this study we assess whether such models are able to describe the phase space properties of evolutionary $N$-body models. We compare the multimass models as implemented in (Gieles \& Zocchi) to $N$-body models of star clusters with different retention… ▽ More Lowered isothermal models, such as the multimass Michie-King models, have been successful in describing observational data of globular clusters. In this study we assess whether such models are able to describe the phase space properties of evolutionary $N$-body models. We compare the multimass models as implemented in (Gieles \& Zocchi) to $N$-body models of star clusters with different retention fractions for the black holes and neutron stars evolving in a tidal field. We find that multimass models successfully reproduce the density and velocity dispersion profiles of the different mass components in all evolutionary phases and for different remnants retention. We further use these results to study the evolution of global model parameters. We find that over the lifetime of clusters, radial anisotropy gradually evolves from the low-mass to the high-mass components and we identify features in the properties of observable stars that are indicative of the presence of stellar-mass black holes. We find that the model velocity scale depends on mass as $m^{-δ}$, with $δ\simeq0.5$ for almost all models, but the dependence of central velocity dispersion on $m$ can be shallower, depending on the dark remnant content, and agrees well with that of the $N$-body models. The reported model parameters, and correlations amongst them, can be used as theoretical priors when fitting these types of mass models to observational data. △ Less

Submitted 18 July, 2017; v1 submitted 6 February, 2017; originally announced February 2017.

Comments: 28 pages, 22 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2017, Volume 470, Issue 3, p.2736-2761

Radial anisotropy in omega Cen limiting the room for an intermediate-mass black hole

Authors: Alice Zocchi, Mark Gieles, Vincent Hénault-Brunet

Abstract: Finding an intermediate-mass black hole (IMBH) in a globular cluster (or proving its absence) would provide valuable insights into our understanding of galaxy formation and evolution. However, it is challenging to identify a unique signature of an IMBH that cannot be accounted for by other processes. Observational claims of IMBH detection are indeed often based on analyses of the kinematics of sta… ▽ More Finding an intermediate-mass black hole (IMBH) in a globular cluster (or proving its absence) would provide valuable insights into our understanding of galaxy formation and evolution. However, it is challenging to identify a unique signature of an IMBH that cannot be accounted for by other processes. Observational claims of IMBH detection are indeed often based on analyses of the kinematics of stars in the cluster core, the most common signature being a rise in the velocity dispersion profile towards the centre of the system. Unfortunately, this IMBH signal is degenerate with the presence of radially-biased pressure anisotropy in the globular cluster. To explore the role of anisotropy in shaping the observational kinematics of clusters, we analyse the case of omega Cen by comparing the observed profiles to those calculated from the family of LIMEPY models, that account for the presence of anisotropy in the system in a physically motivated way. The best-fit radially anisotropic models reproduce the observational profiles well, and describe the central kinematics as derived from Hubble Space Telescope proper motions without the need for an IMBH. △ Less

Submitted 2 February, 2017; originally announced February 2017.

Comments: 14 pages, 7 figures. Accepted for publication in MNRAS

Differences in the rotational properties of multiple stellar populations in M 13: a faster rotation for the "extreme" chemical subpopulation

Authors: M. J. Cordero, V. Hénault-Brunet, C. A. Pilachowski, E. Balbinot, C. I. Johnson, A. L. Varri

Abstract: We use radial velocities from spectra of giants obtained with the WIYN telescope, coupled with existing chemical abundance measurements of Na and O for the same stars, to probe the presence of kinematic differences among the multiple populations of the globular cluster (GC) M13. To characterise the kinematics of various chemical subsamples, we introduce a method using Bayesian inference along with… ▽ More We use radial velocities from spectra of giants obtained with the WIYN telescope, coupled with existing chemical abundance measurements of Na and O for the same stars, to probe the presence of kinematic differences among the multiple populations of the globular cluster (GC) M13. To characterise the kinematics of various chemical subsamples, we introduce a method using Bayesian inference along with an MCMC algorithm to fit a six-parameter kinematic model (including rotation) to these subsamples. We find that the so-called "extreme" population (Na-enhanced and extremely O-depleted) exhibits faster rotation around the centre of the cluster than the other cluster stars, in particular when compared to the dominant "intermediate" population (moderately Na-enhanced and O-depleted). The most likely difference between the rotational amplitude of this extreme population and that of the intermediate population is found to be $\sim$4 km s$^{-1}$, with a 98.4% probability that the rotational amplitude of the extreme population is larger than that of the intermediate population. We argue that the observed difference in rotational amplitudes, obtained when splitting subsamples according to their chemistry, is not a product of the long-term dynamical evolution of the cluster, but more likely a surviving feature imprinted early in the formation history of this GC and its multiple populations. We also find an agreement (within uncertainties) in the inferred position angle of the rotation axis of the different subpopulations considered. We discuss the constraints that these results may place on various formation scenarios. △ Less

Submitted 28 October, 2016; originally announced October 2016.

Comments: Accepted for publication in MNRAS

The Tarantula Massive Binary Monitoring: I. Observational campaign and OB-type spectroscopic binaries

Authors: L. A. Almeida, H. Sana, W. Taylor, R. Barbá, A. Bonanos, P. Crowther, A. Damineli, A. de Koter, S. E. de Mink, C. J. Evans, M. Gieles, N. J. Grin, V. Hénault-Brunet, N. Langer, D. Lennon, S. Lockwood, J. Maíz Apellániz, A. F. J. Moffat, C. Neijssel, C. Norman, O. H. Ramírez-Agudelo, N. D. Richardson, A. Schootemeijer, T. Shenar, I. Soszyński , et al. (2 additional authors not shown)

Abstract: Massive binaries (MBs) play a crucial role in the Universe. Knowing the distributions of their orbital parameters (OPs) is important for a wide range of topics, from stellar feedback to binary evolution channels, from the distribution of supernova types to gravitational wave progenitors, yet, no direct measurements exist outside the Milky Way. The Tarantula Massive Binary Monitoring was designed t… ▽ More Massive binaries (MBs) play a crucial role in the Universe. Knowing the distributions of their orbital parameters (OPs) is important for a wide range of topics, from stellar feedback to binary evolution channels, from the distribution of supernova types to gravitational wave progenitors, yet, no direct measurements exist outside the Milky Way. The Tarantula Massive Binary Monitoring was designed to help fill this gap by obtaining multi-epoch radial velocity monitoring of 102 MBs in the 30 Dor. In this paper, we analyse 32 VLT/FLAMES observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single- and 31 double-lined spectroscopic binaries. Overall, the OPs and binary fraction are remarkably similar across the 30 Dor region and compared to existing Galactic samples (GSs). This indicates that within these domains environmental effects are of second order in shaping the properties of MBs. A small difference is found in the distribution of orbital periods (OrbPs), which is slightly flatter (in log space) in 30 Dor than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, OrbPs in 30 Dor can be as short as 1.1 d; somewhat shorter than seen in GSs. Equal mass binaries q>0.95 in 30 Dor are all found outside NGC 2070 the very young and massive cluster at 30 Dor's core. One outstanding exception however is the fact that earliest spectral types (O2-O7) tend to have shorter OrbPs than latter (O9.2-O9.7). Our results point to a relative universality of the incidence rate of MBs and their OPs in the metallicity range from solar ($Z_{\odot}$) to about $0.5Z_{\odot}$. This provides the first direct constraints on MB properties in massive star-forming galaxies at the Universes peak of star formation at redshifts z~1 to 2, which are estimated to have $Z~0.5Z_{\odot}$. △ Less

Submitted 13 March, 2017; v1 submitted 11 October, 2016; originally announced October 2016.

Comments: 36 pages, 19 figures, 10 tables, published in A\&A

Journal ref: A&A 598, A84 (2017)

A stellar-mass black hole population in the globular cluster NGC 6101?

Authors: Miklos Peuten, Alice Zocchi, Mark Gieles, Alessia Gualandris, Vincent Henault-Brunet

Abstract: Dalessandro et al. observed a similar distribution for blue straggler stars and main-sequence turn-off stars in the Galactic globular cluster NGC 6101, and interpreted this feature as an indication that this cluster is not mass-segregated. Using direct N-body simulations, we find that a significant amount of mass segregation is expected for a cluster with the mass, radius and age of NGC 6101. Ther… ▽ More Dalessandro et al. observed a similar distribution for blue straggler stars and main-sequence turn-off stars in the Galactic globular cluster NGC 6101, and interpreted this feature as an indication that this cluster is not mass-segregated. Using direct N-body simulations, we find that a significant amount of mass segregation is expected for a cluster with the mass, radius and age of NGC 6101. Therefore, the absence of mass segregation cannot be explained by the argument that the cluster is not yet dynamically evolved. By varying the retention fraction of stellar-mass black holes, we show that segregation is not observable in clusters with a high black hole retention fraction (>50% after supernova kicks and >50% after dynamical evolution). Yet all model clusters have the same amount of mass segregation in terms of the decline of the mean mass of stars and remnants with distance to the centre. We also discuss how kinematics can be used to further constrain the presence of a stellar-mass black hole population and distinguish it from the effect of an intermediate-mass black hole. Our results imply that the kick velocities of black holes are lower than those of neutron stars. The large retention fraction during its dynamical evolution can be explained if NGC 6101 formed with a large initial radius in a Milky Way satellite. △ Less

Submitted 6 September, 2016; originally announced September 2016.

Comments: 10 pages, 7 figures

Journal ref: MNRAS (November 01, 2016) 462 (3): 2333

Testing lowered isothermal models with direct N-body simulations of globular clusters

Authors: Alice Zocchi, Mark Gieles, Vincent Hénault-Brunet, Anna Lisa Varri

Abstract: Several self-consistent models have been proposed, aiming at describing the phase space distribution of stars in globular clusters. This study explores the ability of the recently proposed LIMEPY models (Gieles & Zocchi) to reproduce the dynamical properties of direct N-body models of a cluster in a tidal field, during its entire evolution. These dynamical models include prescriptions for the trun… ▽ More Several self-consistent models have been proposed, aiming at describing the phase space distribution of stars in globular clusters. This study explores the ability of the recently proposed LIMEPY models (Gieles & Zocchi) to reproduce the dynamical properties of direct N-body models of a cluster in a tidal field, during its entire evolution. These dynamical models include prescriptions for the truncation and the degree of radially-biased anisotropy contained in the system, allowing us to explore the interplay between the role of anisotropy and tides in various stages of the life of star clusters. We show that the amount of anisotropy in an initially tidally underfilling cluster increases in the pre-collapse phase, and then decreases with time, due to the effect of the external tidal field on its spatial truncation. This is reflected in the correspondent model parameters, and the best-fit models reproduce the main properties of the cluster at all stages of its evolution, except for the phases immediately preceding and following core collapse. We also notice that the best-fit LIMEPY models are significantly different from isotropic King models, especially in the first part of the evolution of the cluster. Our results put limits on the amount of radial anisotropy that can be expected for clusters evolving in a tidal field, which is important to understand other factors that could give rise to similar observational signatures, such as the presence of an intermediate-mass black hole. △ Less

Submitted 6 May, 2016; originally announced May 2016.

Comments: 20 pages, 11 figures. Accepted for publication in MNRAS

No evidence for younger stellar generations within the intermediate age massive clusters NGC 1783, NGC 1696 and NGC 411

Authors: I. Cabrera-Ziri, F. Niederhofer, N. Bastian, M. Rejkuba, E. Balbinot, W. E. Kerzendorf, S. S. Larsen, A. D. Mackey, E. Dalessandro, A. Mucciarelli, C. Charbonnel, M. Hilker, M. Gieles, V. Hénault-Brunet

Abstract: Recently, Li et al. (2016) claimed to have found evidence for multiple generations of stars in the intermediate age clusters NGC 1783, NGC 1696 and NGC 411 in the Large and Small Magellanic Clouds (LMC/SMC). Here we show that these young stellar populations are present in the field regions around these clusters and are not likely associated with the clusters themselves. Using the same datasets, we… ▽ More Recently, Li et al. (2016) claimed to have found evidence for multiple generations of stars in the intermediate age clusters NGC 1783, NGC 1696 and NGC 411 in the Large and Small Magellanic Clouds (LMC/SMC). Here we show that these young stellar populations are present in the field regions around these clusters and are not likely associated with the clusters themselves. Using the same datasets, we find that the background subtraction method adopted by the authors does not adequately remove contaminating stars in the small number Poisson limit. Hence, we conclude that their results do not provide evidence of young generations of stars within these clusters. △ Less

Submitted 20 April, 2016; originally announced April 2016.

Comments: 7 pages, 8 figures, accepted for publication in MNRAS

Chemistry and Kinematics of Red Supergiant Stars in the Young Massive Cluster NGC 2100

Authors: L. R. Patrick, C. J. Evans, B. Davies, R-P. Kudritzki, V. Hénault-Brunet, N. Bastian, E. Lapenna, M. Bergemann

Abstract: We have obtained K-band Multi-Object Spectrograph (KMOS) near-IR spectroscopy for 14 red supergiant stars (RSGs) in the young massive star cluster NGC 2100 in the Large Magellanic Cloud (LMC). Stellar parameters including metallicity are estimated using the J-band analysis technique, which has been rigorously tested in the Local Universe. We find an average metallicity for NGC 2100 of [Z]=$-$0.43… ▽ More We have obtained K-band Multi-Object Spectrograph (KMOS) near-IR spectroscopy for 14 red supergiant stars (RSGs) in the young massive star cluster NGC 2100 in the Large Magellanic Cloud (LMC). Stellar parameters including metallicity are estimated using the J-band analysis technique, which has been rigorously tested in the Local Universe. We find an average metallicity for NGC 2100 of [Z]=$-$0.43$\pm$0.10 dex, in good agreement with estimates from the literature for the LMC. Comparing our results in NGC 2100 with those for a Galactic cluster (at Solar-like metallicity) with a similar mass and age we find no significant difference in the location of RSGs in the Hertzsprung--Russell diagram. We combine the observed KMOS spectra to form a simulated integrated-light cluster spectrum and show that, by analysing this spectrum as a single RSG, the results are consistent with the average properties of the cluster. Radial velocities are estimated for the targets and the dynamical properties are estimated for the first time within this cluster. The data are consistent with a flat velocity dispersion profile, and with an upper limit of 3.9 km/s, at the 95% confidence level, for the velocity dispersion of the cluster. However, the intrinsic velocity dispersion is unresolved and could, therefore, be significantly smaller than the upper limit reported here. An upper limit on the dynamical mass of the cluster is derived as $M_{dyn}$ $\le$ $15.2\times10^{4}M_{\odot}$ assuming virial equilibrium. △ Less

Submitted 14 April, 2016; v1 submitted 8 February, 2016; originally announced February 2016.

Comments: Accepted in MNRAS

Multiple populations in globular clusters: constraints from kinematics and dynamics

Authors: Vincent Hénault-Brunet

Abstract: We discuss constraints on the formation of multiple populations in globular clusters (GCs) imposed by their present-day kinematics (velocity dispersion and anisotropy) and spatial distribution. We argue that the observational evidence collected so far in the outer parts of clusters is generally consistent with an enriched population forming more centrally concentrated compared to the primordial po… ▽ More We discuss constraints on the formation of multiple populations in globular clusters (GCs) imposed by their present-day kinematics (velocity dispersion and anisotropy) and spatial distribution. We argue that the observational evidence collected so far in the outer parts of clusters is generally consistent with an enriched population forming more centrally concentrated compared to the primordial population, in agreement with all the scenarios proposed to date (in some cases by design), but not sufficient to favour a particular scenario. We highlight that the differential rotation of subpopulations is a signature that may provide crucial new constraints and allow us to distinguish between various scenarios. Finally, we discuss the spatial distribution of subpopulations in the central regions of GCs and speculate that mass segregation between subpopulations may be due to a difference in their binary fraction. △ Less

Submitted 7 December, 2015; originally announced December 2015.

Comments: 6 pages, 3 figures, to appear in the proceedings of IAU symp. 316, "Formation, evolution, and survival of massive star clusters", eds. C. Charbonnel & A. Nota

The VLT-FLAMES Tarantula Survey XXI. Stellar spin rates of O-type spectroscopic binaries

Authors: O. H. Ramírez-Agudelo, H. Sana, S. E. de Mink, V. Hénault-Brunet, A. de Koter, N. Langer, F. Tramper, G. Gräfener, C. J. Evans, J. S. Vink, P. L. Dufton, W. D. Taylor

Abstract: The initial distribution of spin rates of massive stars is a fingerprint of their elusive formation process. It also sets a key initial condition for stellar evolution and is thus an important ingredient in stellar population synthesis. So far, most studies have focused on single stars. Most O stars are however found in multiple systems. By establishing the spin-rate distribution of a sizeable sam… ▽ More The initial distribution of spin rates of massive stars is a fingerprint of their elusive formation process. It also sets a key initial condition for stellar evolution and is thus an important ingredient in stellar population synthesis. So far, most studies have focused on single stars. Most O stars are however found in multiple systems. By establishing the spin-rate distribution of a sizeable sample of O-type spectroscopic binaries and by comparing the distributions of binary sub-populations with one another as well as with that of presumed single stars in the same region, we aim to constrain the initial spin distribution of O stars in binaries, and to identify signatures of the physical mechanisms that affect the evolution of the massive stars spin rates. We use ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS) to establish the projected equatorial rotational velocities (\vrot) for components of 114 spectroscopic binaries in 30 Doradus. The \vrot\ values are derived from the full-width at half-maximum (FWHM) of a set of spectral lines, using a FWHM vs. \vrot\ calibration that we derive based on previous line analysis methods applied to single O-type stars in the VFTS sample. The overall \vrot\ distribution of the primary stars resembles that of single O-type stars in the VFTS, featuring a low-velocity peak (at $\vrot < 200$ kms) and a shoulder at intermediate velocities ($200 < \vrot < 300$ kms). The distributions of binaries and single stars however differ in two ways. First, the main peak at $\vrot \sim$100 kms is broader and slightly shifted toward higher spin rates in the binary distribution compared to that of the presumed-single stars. Second, the \vrot distribution of primaries lacks a significant population of stars spinning faster than 300 kms while such a population is clearly present in the single star sample. △ Less

Submitted 8 July, 2015; originally announced July 2015.

Comments: 16 pages, 16 figures, paper accepted in Astronomy & Astrophysics

Journal ref: A&A 580, A92 (2015)

Biases in the determination of dynamical parameters of star clusters: today and in the Gaia era

Authors: A. Sollima, H. Baumgardt, A. Zocchi, E. Balbinot, M. Gieles, V. Henault-Brunet, A. L. Varri

Abstract: The structural and dynamical properties of star clusters are generally derived by means of the comparison between steady-state analytic models and the available observables. With the aim of studying the biases of this approach, we fitted different analytic models to simulated observations obtained from a suite of direct N-body simulations of star clusters in different stages of their evolution and… ▽ More The structural and dynamical properties of star clusters are generally derived by means of the comparison between steady-state analytic models and the available observables. With the aim of studying the biases of this approach, we fitted different analytic models to simulated observations obtained from a suite of direct N-body simulations of star clusters in different stages of their evolution and under different levels of tidal stress to derive mass, mass function and degree of anisotropy. We find that masses can be under/over-estimated up to 50% depending on the degree of relaxation reached by the cluster, the available range of observed masses and distances of radial velocity measures from the cluster center and the strength of the tidal field. The mass function slope appears to be better constrainable and less sensitive to model inadequacies unless strongly dynamically evolved clusters and a non-optimal location of the measured luminosity function are considered. The degree and the characteristics of the anisotropy developed in the N-body simulations are not adequately reproduced by popular analytic models and can be detected only if accurate proper motions are available. We show how to reduce the uncertainties in the mass, mass-function and anisotropy estimation and provide predictions for the improvements expected when Gaia proper motions will be available in the near future. △ Less

Submitted 11 May, 2015; originally announced May 2015.

Comments: 14 pages, 8 figures, accepted for publication by MNRAS

Multiple populations in globular clusters: the distinct kinematic imprints of different formation scenarios

Authors: Vincent Hénault-Brunet, Mark Gieles, Oscar Agertz, Justin I. Read

Abstract: Several scenarios have been proposed to explain the presence of multiple stellar populations in globular clusters. Many of them invoke multiple generations of stars to explain the observed chemical abundance anomalies, but it has also been suggested that self-enrichment could occur via accretion of ejecta from massive stars onto the circumstellar disc of low-mass pre-main sequence stars. These sce… ▽ More Several scenarios have been proposed to explain the presence of multiple stellar populations in globular clusters. Many of them invoke multiple generations of stars to explain the observed chemical abundance anomalies, but it has also been suggested that self-enrichment could occur via accretion of ejecta from massive stars onto the circumstellar disc of low-mass pre-main sequence stars. These scenarios imply different initial conditions for the kinematics of the various stellar populations. Given some net angular momentum initially, models for which a second generation forms from gas that collects in a cooling flow into the core of the cluster predict an initially larger rotational amplitude for the polluted stars compared to the pristine stars. This is opposite to what is expected from the accretion model, where the polluted stars are the ones crossing the core and are on preferentially radial (low-angular momentum) orbits, such that their rotational amplitude is lower. Here we present the results of a suite of $N$-body simulations with initial conditions chosen to capture the distinct kinematic properties of these pollution scenarios. We show that initial differences in the kinematics of polluted and pristine stars can survive to the present epoch in the outer parts of a large fraction of Galactic globular clusters. The differential rotation of pristine and polluted stars is identified as a unique kinematic signature that could allow us to distinguish between various scenarios, while other kinematic imprints are generally very similar from one scenario to the other. △ Less

Submitted 25 March, 2015; originally announced March 2015.

Comments: 22 pages, 16 figures + appendix. Accepted for publication in MNRAS

On the uniqueness of kinematical signatures of intermediate-mass black holes in globular clusters

Authors: Alice Zocchi, Mark Gieles, Vincent Hénault-Brunet

Abstract: Finding an intermediate-mass black hole (IMBH) in a globular cluster (GC), or proving its absence, is a crucial ingredient in our understanding of galaxy formation and evolution. The challenge is to identify a unique signature of an IMBH that cannot be accounted for by other processes. Observational claims of IMBH detection are often based on analyses of the kinematics of stars, such as a rise in… ▽ More Finding an intermediate-mass black hole (IMBH) in a globular cluster (GC), or proving its absence, is a crucial ingredient in our understanding of galaxy formation and evolution. The challenge is to identify a unique signature of an IMBH that cannot be accounted for by other processes. Observational claims of IMBH detection are often based on analyses of the kinematics of stars, such as a rise in the velocity dispersion profile towards the centre. In this contribution we discuss the degeneracy between this IMBH signal and pressure anisotropy in the GC. We show that that by considering anisotropic models it is possible to partially explain the innermost shape of the projected velocity dispersion profile, even though models that do not account for an IMBH do not exhibit a cusp in the centre. △ Less

Submitted 21 January, 2015; originally announced January 2015.

Comments: 4 pages, 2 figures. To be published in the Proceedings IAU Symposium No. 312, Star Clusters and Black Holes in Galaxies Across Cosmic Time

The VLT-FLAMES Tarantula Survey XVIII. Classifications and radial velocities of the B-type stars

Authors: C. J. Evans, M. B. Kennedy, P. L. Dufton, I. D. Howarth, N. R. Walborn, N. Markova, J. S. Clark, S. E. de Mink, A. de Koter, P. R. Dunstall, V. Hénault-Brunet, J. Maíz Apellániz, C. M. McEvoy, H. Sana, S. Simón-Díaz, W. D. Taylor, J. S. Vink

Abstract: We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of the radia… ▽ More We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of the radial velocities across the 30 Dor region, and use the results to identify candidate runaway stars. Excluding potential runaways and members of two older clusters in the survey region (SL 639 and Hodge 301), we determine a systemic velocity for 30 Dor of 271.6 +/- 12.2 km/s from 273 presumed single stars. Employing a 3-sigma criterion we identify nine candidate runaway stars (2.9% of the single stars with radial-velocity estimates). The projected rotational velocities of the candidate runaways appear to be significantly different to those of the full B-type sample, with a strong preference for either large (>345 km/s) or small (<65 km/s) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (-106.9 +/- 16.2 km/s), and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12+log(N/H) > ~8.5. Combined with a large rotational velocity (vsini = 345 +/- 22 km/s), this is suggestive of past binary interaction for this star. △ Less

Submitted 16 January, 2015; originally announced January 2015.

Comments: 39 pages (main body: 14 pages; ancillary material: 25 pages)

Journal ref: 2015, A&A, 574, A13

The VLT-FLAMES Tarantula Survey XIX. B-type Supergiants - Atmospheric parameters and nitrogen abundances to investigate the role of binarity and the width of the main sequence

Authors: C. M. McEvoy, P. L. Dufton, C. J. Evans, V. M. Kalari, N. Markova, S. Simón-Díaz, J. S. Vink, N. R. Walborn, P. A. Crowther, A. de Koter, S. E. de Mink, P. R. Dunstall, V. Hénault-Brunet, A. Herrero, N. Langer, D. J. Lennon, J. Maíz Apellániz, F. Najarro, J. Puls, H. Sana, F. R. N. Schneider, W. D. Taylor

Abstract: TLUSTY non-LTE model atmosphere calculations have been used to determine atmospheric parameters and nitrogen (N) abundances for 34 single and 18 binary B-type supergiants (BSGs). The effects of flux contribution from an unseen secondary were considered for the binary sample. We present the first systematic study of the incidence of binarity for a sample of BSGs across the theoretical terminal age… ▽ More TLUSTY non-LTE model atmosphere calculations have been used to determine atmospheric parameters and nitrogen (N) abundances for 34 single and 18 binary B-type supergiants (BSGs). The effects of flux contribution from an unseen secondary were considered for the binary sample. We present the first systematic study of the incidence of binarity for a sample of BSGs across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the BSGs it may be necessary to extend the TAMS to lower temperatures. This is consistent with the derived distribution of mass discrepancies, projected rotational velocities (vsini) and N abundances, provided that stars cooler than this temperature are post RSG objects. For the BSGs in the Tarantula and previous FLAMES surveys, most have small vsini. About 10% have larger vsini (>100 km/s) but surprisingly these show little or no N enhancement. All the cooler BSGs have low vsini of <70km/s and high N abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. A lack of cool binaries, possibly reflects the small sample size. Single star evolutionary models, which include rotation, can account for the N enhancement in both the single and binary samples. The detailed distribution of N abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. The first comparative study of single and binary BSGs has revealed that the main sequence may be significantly wider than previously assumed, extending to Teff=20000K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolution of BSGs. Full abstract in paper △ Less

Submitted 12 January, 2015; v1 submitted 8 December, 2014; originally announced December 2014.

Comments: 21 pages, 15 figures, 11 tables

Journal ref: A&A 575, A70 (2015)

The VLT-FLAMES Tarantula Survey XVII. Physical and wind properties of massive stars at the top of the main sequence

Authors: Joachim M. Bestenlehner, Götz Gräfener, Jorick S. Vink, F. Najarro, A. de Koter, H. Sana, C. J. Evans, P. A. Crowther, V. Hénault-Brunet, A. Herrero, N. Langer, F. R. N. Schneider, S. Simón-Díaz, W. D. Taylor, N. R. Walborn

Abstract: The evolution and fate of very massive stars (VMS) is tightly connected to their mass-loss properties. Their initial and final masses differ significantly as a result of mass loss. VMS have strong stellar winds and extremely high ionising fluxes, which are thought to be critical sources of both mechanical and radiative feedback in giant Hii regions. However, how VMS mass-loss properties change dur… ▽ More The evolution and fate of very massive stars (VMS) is tightly connected to their mass-loss properties. Their initial and final masses differ significantly as a result of mass loss. VMS have strong stellar winds and extremely high ionising fluxes, which are thought to be critical sources of both mechanical and radiative feedback in giant Hii regions. However, how VMS mass-loss properties change during stellar evolution is poorly understood. In the framework of the VLT-Flames Tarantula Survey (VFTS), we explore the mass-loss transition region from optically thin O to denser WNh star winds, thereby testing theoretical predictions. To this purpose we select 62 O, Of, Of/WN, and WNh stars, an unprecedented sample of stars with the highest masses and luminosities known. We perform a spectral analysis of optical VFTS as well as near-infrared VLT/SINFONI data using the non-LTE radiative transfer code CMFGEN to obtain stellar and wind parameters. For the first time, we observationally resolve the transition between optically thin O and optically thick WNh star winds. Our results suggest the existence of a kink between both mass-loss regimes, in agreement with recent MC simulations. For the optically thick regime, we confirm the steep dependence on the Eddington factor from previous theoretical and observational studies. The transition occurs on the MS near a luminosity of 10^6.1Lsun, or a mass of 80...90Msun. Above this limit, we find that - even when accounting for moderate wind clumping (with f = 0.1) - wind mass-loss rates are enhanced with respect to standard prescriptions currently adopted in stellar evolution calculations. We also show that this results in substantial helium surface enrichment. Based on our spectroscopic analyses, we are able to provide the most accurate ionising fluxes for VMS known to date, confirming the pivotal role of VMS in ionising and shaping their environments. △ Less

Submitted 7 July, 2014; originally announced July 2014.

Comments: Accepted for publication in A&A, 19 pages, 14 figures, 6 tables, (74 pages appendix, 68 figures, 4 tables)

Journal ref: A&A 570, A38 (2014)

Binary-corrected velocity dispersions from single- and multi-epoch radial velocities: massive stars in R136 as a test case

Authors: Michiel Cottaar, Vincent Hénault-Brunet

Abstract: Orbital motions from binary stars can broaden the observed line-of-sight velocity distribution of a stellar system, artificially inflating the measured line-of-sight velocity dispersion, which can in turn lead to erroneous conclusions about the dynamical state of the system. Cottaar et al. (2012b) proposed a maximum likelihood procedure to recover the intrinsic velocity dispersion of a resolved st… ▽ More Orbital motions from binary stars can broaden the observed line-of-sight velocity distribution of a stellar system, artificially inflating the measured line-of-sight velocity dispersion, which can in turn lead to erroneous conclusions about the dynamical state of the system. Cottaar et al. (2012b) proposed a maximum likelihood procedure to recover the intrinsic velocity dispersion of a resolved star cluster from a single epoch of radial velocity data of individual stars, which they achieved by simultaneously fitting the intrinsic velocity distribution of the single stars and the centres of mass of the binaries along with the velocity shifts caused by binary orbital motions. Assuming well-characterized binary properties, they showed that this procedure can accurately reproduce intrinsic velocity dispersions below 1 km s$^{-1}$ for solar-type stars. Here we investigate the systematic offsets induced in cases where the binary properties are uncertain, and we show how two epochs of radial velocity data with an appropriate baseline can help to mitigate these systematic effects. We first test the method above using Monte Carlo simulations, taking into account the large uncertainties in the binary properties of OB stars. We then apply it to radial velocity data in the young massive cluster R136, an example for which the intrinsic velocity dispersion of O-type stars is known from an intensive multi-epoch approach. For typical velocity dispersions of young massive clusters ($\gtrsim 4$ km s$^{-1}$) and with a single epoch of data, we demonstrate that the method can just about distinguish between a cluster in virial equilibrium and an unbound cluster. This is due to the higher spectroscopic binary fraction and more loosely constrained distributions of orbital parameters of OB stars compared to solar-type stars. By extending the maximum likelihood method to multi-epoch data, <abridged> △ Less

Submitted 26 November, 2013; v1 submitted 25 November, 2013; originally announced November 2013.

Comments: Accepted by A&A; minor corrections made on November 27

Studying the kinematics of the giant star-forming region 30 Doradus. I. The data

Authors: S. Torres-Flores, R. Barbá, J. Maíz Apellániz, M. Rubio, G. Bosch, V. Hénault-Brunet, C. J. Evans

Abstract: We present high-quality VLT-FLAMES optical spectroscopy of the nebular gas in the giant star-forming region 30 Doradus. In this paper, the first of a series, we introduce our observations and discuss the main kinematic features of 30 Dor, as revealed by the spectroscopy of the ionized gas in the region. The primary data set consists of regular grid of nebular observations, which we used to produce… ▽ More We present high-quality VLT-FLAMES optical spectroscopy of the nebular gas in the giant star-forming region 30 Doradus. In this paper, the first of a series, we introduce our observations and discuss the main kinematic features of 30 Dor, as revealed by the spectroscopy of the ionized gas in the region. The primary data set consists of regular grid of nebular observations, which we used to produce a spectroscopic datacube of 30 Dor, centered on the massive star cluster R136 and covering a field-of-view of 10'x10'. The main emission lines present in the datacube are from Halpha and [NII]6548,6584. The Halpha emission-line profile varies across the region from simple single-peaked emission to complex, multiple-component profiles, suggesting that different physical mechanisms are acting on the excited gas. To analyse the gas kinematics we fit Gaussian profiles to the observed Halpha features. Unexpectedly, the narrowest Halpha profile in our sample lies close to the supernova remnant 30 Dor B. We present maps of the velocity field and velocity dispersion across 30 Dor, finding five previously unclassified expanding structures. These maps highlight the kinematic richness of 30 Dor (e.g. supersonic motions), which will be analysed in future papers. △ Less

Submitted 30 April, 2013; originally announced May 2013.

Comments: 24 pages, 12 figures. Accepted for publication in Astronomy and Astrophysics

Long-term evolution of the neutron-star spin period of SXP 1062

Authors: R. Sturm, F. Haberl, L. M. Oskinova, M. P. E. Schurch, V. Hénault-Brunet, J. S. Gallagher III, A. Udalski

Abstract: The Be/X-ray binary SXP1062 is of especial interest owing to the large spin period of the neutron star, its large spin-down rate, and the association with a supernova remnant constraining its age. This makes the source an important probe for accretion physics. To investigate the long-term evolution of the spin period and associated spectral variations, we performed an XMM-Newton target-of-opportun… ▽ More The Be/X-ray binary SXP1062 is of especial interest owing to the large spin period of the neutron star, its large spin-down rate, and the association with a supernova remnant constraining its age. This makes the source an important probe for accretion physics. To investigate the long-term evolution of the spin period and associated spectral variations, we performed an XMM-Newton target-of-opportunity observation of SXP1062 during X-ray outburst. Spectral and timing analysis of the XMM-Newton data was compared with previous studies, as well as complementary Swift/XRT monitoring and optical spectroscopy with the SALT telescope were obtained. The spin period was measured to be P_s =(1071.01+-0.16) s on 2012 Oct 14. The X-ray spectrum is similar to that of previous observations. No convincing cyclotron absorption features, which could be indicative for a high magnetic field strength, are found. The high-resolution RGS spectra indicate the presence of emission lines, which may not completely be accounted for by the SNR emission. The comparison of multi-epoch optical spectra suggest an increasing size or density of the decretion disc around the Be star. SXP1062 showed a net spin-down with an average of (2.27+-0.44) s/yr over a baseline of 915 days. △ Less

Submitted 18 July, 2013; v1 submitted 22 April, 2013; originally announced April 2013.

Comments: 8 pages, 6 figures, accepted for publication in A&A

Discovery of X-ray emission from young suns in the Small Magellanic Cloud

Authors: L. M. Oskinova, W. Sun, C. J. Evans, V. Hénault-Brunet, Y. -H. Chu, J. S. Gallagher III, M. A. Guerrero, R. A. Gruendl, M. Güdel, S. Silich, Y. Chen, Y. Nazé, R. Hainich, J. Reyes-Iturbide

Abstract: We report the discovery of extended X-ray emission within the young star cluster NGC 602 in the Wing of the Small Magellanic Cloud (SMC) based on observations obtained with the Chandra X-ray Observatory. X-ray emission is detected from the cluster core area with the highest stellar density and from a dusty ridge surrounding the HII region. We use a census of massive stars in the cluster to demonst… ▽ More We report the discovery of extended X-ray emission within the young star cluster NGC 602 in the Wing of the Small Magellanic Cloud (SMC) based on observations obtained with the Chandra X-ray Observatory. X-ray emission is detected from the cluster core area with the highest stellar density and from a dusty ridge surrounding the HII region. We use a census of massive stars in the cluster to demonstrate that a cluster wind or wind-blown bubble is unlikely to provide a significant contribution to the X-ray emission detected from the central area of the cluster. We therefore suggest that X-ray emission at the cluster core originates from an ensemble of low- and solar-mass pre-main-sequence (PMS) stars, each of which would be too weak in X-rays to be detected individually. We attribute the X-ray emission from the dusty ridge to the embedded tight cluster of the new-born stars known in this area from infrared studies. Assuming that the levels of X-ray activity in young stars in the low-metallicity environment of NGC 602a are comparable to their Galactic counterparts, then the detected spatial distribution, spectral properties, and level of X-ray emission are largely consistent with those expected from low- and solar-mass PMS stars and young stellar objects (YSOs). This is the first discovery of X-ray emission attributable to PMS stars and YSOs in the SMC, which suggests that the accretion and dynamo processes in young, low-mass objects in the SMC resemble those in the Galaxy. △ Less

Submitted 15 January, 2013; originally announced January 2013.

Comments: Accepted for publication in ApJ

Journal ref: ApJ, 765, 73 (12p, 2013)

The slow X-ray pulsar SXP 1062 and associated supernova remnant in the Wing of the Small Magellanic Cloud

Authors: L. M. Oskinova, M. A. Guerrero, V. Henault-Brunet, W. Sun, Y. -H. Chu, C. Evans, J. S. Gallagher III, R. A. Gruendl, J. Reyes-Iturbide

Abstract: SXP 1062 is an exceptional case of a young neutron star in a wind-fed high-mass X-ray binary associated with a supernova remnant. A unique combination of measured spin period, its derivative, luminosity and young age makes this source a key probe for the physics of accretion and neutron star evolution. Theoretical models proposed to explain the properties of SXP 1062 shall be tested with new data. SXP 1062 is an exceptional case of a young neutron star in a wind-fed high-mass X-ray binary associated with a supernova remnant. A unique combination of measured spin period, its derivative, luminosity and young age makes this source a key probe for the physics of accretion and neutron star evolution. Theoretical models proposed to explain the properties of SXP 1062 shall be tested with new data. △ Less

Submitted 18 October, 2012; originally announced October 2012.

Comments: Proceedings of IAUS 291 "Neutron Stars and Pulsars: Challenges and Opportunities after 80 years", J. van Leeuwen (ed.)

The VLT-FLAMES Tarantula Survey. IX. The interstellar medium seen through Diffuse Interstellar Bands and neutral sodium

Authors: Jacco Th. van Loon, Mandy Bailey, Benjamin L. Tatton, Jesus Maiz Apellaniz, Paul A. Crowther, Alex de Koter, Christopher J. Evans, Vincent Henault-Brunet, Ian D. Howarth, Philipp Richter, Hugues Sana, Sergio Simon-Diaz, William Taylor, Nolan R. Walborn

Abstract: The Tarantula Nebula (30 Dor) is a spectacular star-forming region in the Large Magellanic Cloud, seen through gas in the Galactic Disc and Halo. Diffuse Interstellar Bands offer a unique probe of the diffuse, cool-warm gas in these regions. The aim is to use DIBs as diagnostics of the local interstellar conditions, whilst at the same time deriving properties of the yet-unknown carriers. Spectra o… ▽ More The Tarantula Nebula (30 Dor) is a spectacular star-forming region in the Large Magellanic Cloud, seen through gas in the Galactic Disc and Halo. Diffuse Interstellar Bands offer a unique probe of the diffuse, cool-warm gas in these regions. The aim is to use DIBs as diagnostics of the local interstellar conditions, whilst at the same time deriving properties of the yet-unknown carriers. Spectra of over 800 early-type stars from the VLT Flames Tarantula Survey (VFTS) were analysed. Maps were created, separately, for the Galactic and LMC absorption in the DIBs at 4428 and 6614 Ang and - in a smaller region near the central cluster R136 - neutral sodium (Na I D); we also measured the DIBs at 5780 and 5797 Ang. The maps show strong 4428 and 6614 Ang DIBs in the quiescent cloud complex to the south of 30 Dor but weak absorption in the harsher environments to the north (bubbles) and near the OB associations. The Na maps show at least five kinematic components in the LMC and a shell-like structure surrounding R136, and small-scale structure in the Milky Way. The strengths of the 4428, 5780, 5797 and 6614 Ang DIBs are correlated, also with Na absorption and visual extinction. The strong 4428 Ang DIB is present already at low Na column density but the 6614, 5780 and 5797 Ang DIBs start to be detectable at subsequently larger Na column densities. The relative strength of the 5780 and 5797 Ang DIBs clearly confirm the Tarantula Nebula and Galactic high-latitude gas to represent a harsh radiation environment. The resilience of the 4428 Ang DIB suggests its carrier is large, compact and neutral. Structure is detected in the distribution of cool-warm gas on scales between one and >100 pc in the LMC and as little as 0.01 pc in the Sun's vicinity. Stellar winds from the central cluster R136 have created an expanding shell; some infalling gas is also detected, reminiscent of a galactic "fountain". △ Less

Submitted 16 October, 2012; originally announced October 2012.

Comments: Accepted for publication in Astronomy and Astrophysics

The VLT-FLAMES Tarantula Survey VIII. Multiplicity properties of the O-type star population

Authors: H. Sana, A. de Koter, S. E. de Mink, P. R. Dunstall, C. J. Evans, V. Henault-Brunet, J. Maiz Apellaniz, O. H. Ramirez-Agudelo, W. D. Taylor, N. R. Walborn, J. S. Clark, P. A. Crowther, A. Herrero, M. Gieles, N. Langer, D. J. Lennon, J. S. Vink

Abstract: Aims. We analyze the multiplicity properties of the massive O-type star population. With 360 O-type stars, this is the largest homogeneous sample of massive stars analyzed to date. Methods. We use multi-epoch spectroscopy and variability analysis to identify spectroscopic binaries. We also use a Monte-Carlo method to correct for observational biases. Results. We observe a spectroscopic binary… ▽ More Aims. We analyze the multiplicity properties of the massive O-type star population. With 360 O-type stars, this is the largest homogeneous sample of massive stars analyzed to date. Methods. We use multi-epoch spectroscopy and variability analysis to identify spectroscopic binaries. We also use a Monte-Carlo method to correct for observational biases. Results. We observe a spectroscopic binary fraction of 0.35\pm0.03, which corresponds to the fraction of objects displaying statistically significant radial velocity variations with an amplitude of at least 20km/s. We compute the intrinsic binary fraction to be 0.51\pm0.04. We adopt power-laws to describe the intrinsic period and mass-ratio distributions: f_P ~ (log P)^π (with 0.15 < log P < 3.5) and f_q ~ q^κ with 0.1 < q=M_2/M_1 < 1.0. The power-law indexes that best reproduce the observed quantities are π= -0.45 +/- 0.30 and κ= -1.0\pm0.4. The obtained period distribution thus favours shorter period systems compared to an Oepik law. The mass ratio distribution is slightly skewed towards low mass ratio systems but remains incompatible with a random sampling of a classical mass function. The binary fraction seems mostly uniform across the field of view and independent of the spectral types and luminosity classes. The binary fraction in the outer region of the field of view (r > 7.8', i.e. approx117 pc) and among the O9.7 I/II objects are however significantly lower than expected from statistical fluctuations. Conclusions. Using simple evolutionary considerations, we estimate that over 50% of the current O star population in 30 Dor will exchange mass with its companion within a binary system. This shows that binary interaction is greatly affecting the evolution and fate of massive stars, and must be taken into account to correctly interpret unresolved populations of massive stars. △ Less

Submitted 20 September, 2012; originally announced September 2012.

Comments: Accepted for publication in A&A. 23 pages, 20 figures