TY - JOUR A1 - Ramiaramanantsoa, Tahina A1 - Moffat, Anthony F. J. A1 - Harmon, Robert A1 - Ignace, R. A1 - St-Louis, Nicole A1 - Vanbeveren, Dany A1 - Shenar, Tomer A1 - Pablo, Herbert A1 - Richardson, Noel D. A1 - Howarth, Ian D. A1 - Stevens, Ian R. A1 - Piaulet, Caroline A1 - St-Jean, Lucas A1 - Eversberg, Thomas A1 - Pigulski, Andrzej A1 - Popowicz, Adam A1 - Kuschnig, Rainer A1 - Zoclonska, Elzbieta A1 - Buysschaert, Bram A1 - Handler, Gerald A1 - Weiss, Werner W. A1 - Wade, Gregg A. A1 - Rucinski, Slavek M. A1 - Zwintz, Konstanze A1 - Luckas, Paul A1 - Heathcote, Bernard A1 - Cacella, Paulo A1 - Powles, Jonathan A1 - Locke, Malcolm A1 - Bohlsen, Terry A1 - Chené, André-Nicolas A1 - Miszalski, Brent A1 - Waldron, Wayne L. A1 - Kotze, Marissa M. A1 - Kotze, Enrico J. A1 - Böhm, Torsten T1 - BRITE-Constellation high-precision time-dependent photometry of the early O-type supergiant zeta Puppis unveils the photospheric drivers of its small- and large-scale wind structures JF - Monthly notices of the Royal Astronomical Society N2 - From 5.5 months of dual-band optical photometric monitoring at the 1 mmag level, BRITE-Constellation has revealed two simultaneous types of variability in the O4I(n)fp star ζ Puppis: one single periodic non-sinusoidal component superimposed on a stochastic component. The monoperiodic component is the 1.78-d signal previously detected by Coriolis/Solar Mass Ejection Imager, but this time along with a prominent first harmonic. The shape of this signal changes over time, a behaviour that is incompatible with stellar oscillations but consistent with rotational modulation arising from evolving bright surface inhomogeneities. By means of a constrained non-linear light-curve inversion algorithm, we mapped the locations of the bright surface spots and traced their evolution. Our simultaneous ground-based multisite spectroscopic monitoring of the star unveiled cyclical modulation of its He ii λ4686 wind emission line with the 1.78-d rotation period, showing signatures of corotating interaction regions that turn out to be driven by the bright photospheric spots observed by BRITE. Traces of wind clumps are also observed in the He ii λ4686 line and are correlated with the amplitudes of the stochastic component of the light variations probed by BRITE at the photosphere, suggesting that the BRITE observations additionally unveiled the photospheric drivers of wind clumps in ζ Pup and that the clumping phenomenon starts at the very base of the wind. The origins of both the bright surface inhomogeneities and the stochastic light variations remain unknown, but a subsurface convective zone might play an important role in the generation of these two types of photospheric variability. KW - techniques: photometric KW - techniques: spectroscopic KW - stars: massive KW - stars: rotation KW - starspots KW - supergiants KW - stars: winds, outflows Y1 - 2017 U6 - https://doi.org/10.1093/mnras/stx2671 SN - 0035-8711 SN - 1365-2966 VL - 473 IS - 4 SP - 5532 EP - 5569 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Sander, Andreas Alexander Christoph A1 - Fürst, F. A1 - Kretschmar, P. A1 - Oskinova, Lida A1 - Todt, Helge Tobias A1 - Hainich, Rainer A1 - Shenar, Tomer A1 - Hamann, Wolf-Rainer T1 - Coupling hydrodynamics with comoving frame radiative transfer BT - Stellar wind stratification in the high-mass X-ray binary Vela X-1 JF - Astronomy and astrophysics : an international weekly journal N2 - Aims. To gain a realistic picture of the donor star in Vela X-1, we constructed a hydrodynamically consistent atmosphere model describing the wind stratification while properly reproducing the observed donor spectrum. To investigate how X-ray illumination affects the stellar wind, we calculated additional models for different X-ray luminosity regimes. Methods. We used the recently updated version of the Potsdam Wolf-Rayet code to consistently solve the hydrodynamic equation together with the statistical equations and the radiative transfer. Results. The wind flow in Vela X-1 is driven by ions from various elements, with Fe III and S III leading in the outer wind. The model-predicted mass-loss rate is in line with earlier empirical studies. The mass-loss rate is almost unaffected by the presence of the accreting NS in the wind. The terminal wind velocity is confirmed at u(infinity) approximate to 600 km s(-1). On the other hand, the wind velocity in the inner region where the NS is located is only approximate to 100 km s(-1), which is not expected on the basis of a standard beta-velocity law. In models with an enhanced level of X-rays, the velocity field in the outer wind can be altered. If the X-ray flux is too high, the acceleration breaks down because the ionization increases. Conclusions. Accounting for radiation hydrodynamics, our Vela X-1 donor atmosphere model reveals a low wind speed at the NS location, and it provides quantitative information on wind driving in this important HMXB. KW - stars: mass-loss KW - stars: winds, outflows KW - stars: early-type KW - stars: atmospheres KW - stars: massive KW - X-rays: binaries Y1 - 2018 U6 - https://doi.org/10.1051/0004-6361/201731575 SN - 1432-0746 VL - 610 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Ramachandran, Varsha A1 - Hamann, Wolf-Rainer A1 - Hainich, Rainer A1 - Oskinova, Lida A1 - Shenar, Tomer A1 - Sander, Andreas Alexander Christoph A1 - Todt, Helge Tobias A1 - Gallagher, John S. T1 - Stellar population of the superbubble N206 in the LMC II. Parameters of the OB and WR stars, and the total massive star feedback JF - Astronomy and astrophysics : an international weekly journal N2 - Context. Clusters or associations of early-type stars are often associated with a "superbubble" of hot gas. The formation of such superbubbles is caused by the feedback from massive stars. The complex N206 in the Large Magellanic Cloud (LMC) exhibits a superbubble and a rich massive star population. Aims. Our goal is to perform quantitative spectral analyses of all massive stars associated with the N206 superbubble in order to determine their stellar and wind parameters. We compare the superbubble energy budget to the stellar energy input and discuss the star formation history of the region. Results. We present the stellar and wind parameters of the OB stars and the two Wolf-Rayet (WR) binaries in the N206 complex. Twelve percent of the sample show Oe/Be type emission lines, although most of them appear to rotate far below critical. We found eight runaway stars based on their radial velocity. The wind-momentum luminosity relation of our OB sample is consistent with the expectations. The Hertzsprung-Russell diagram (HRD) of the OB stars reveals a large age spread (1-30 Myr), suggesting different episodes of star formation in the complex. The youngest stars are concentrated in the inner part of the complex, while the older OB stars are scattered over outer regions. We derived the present day mass function for the entire N206 complex as well as for the cluster NGC2018. The total ionizing photon flux produced by all massive stars in the N206 complex is Q(0) approximate to 5 x 10(50) s(-1), and the mechanical luminosity of their stellar winds amounts to L-mec = 1.7 x 10(38) erg s(-1). Three very massive Of stars are found to dominate the feedback among 164 OB stars in the sample. The two WR winds alone release about as much mechanical luminosity as the whole OB star sample. The cumulative mechanical feedback from all massive stellar winds is comparable to the combined mechanical energy of the supernova explosions that likely occurred in the complex. Accounting also for the WR wind and supernovae, the mechanical input over the last five Myr is approximate to 2.3 x 10(52) erg. Conclusions. The N206 complex in the LMC has undergone star formation episodes since more than 30 Myr ago. From the spectral analyses of its massive star population, we derive a current star formation rate of 2.2 x 10(-3) M-circle dot yr(-1). From the combined input of mechanical energy from all stellar winds, only a minor fraction is emitted in the form of X-rays. The corresponding input accumulated over a long time also exceeds the current energy content of the complex by more than a factor of five. The morphology of the complex suggests a leakage of hot gas from the superbubble. KW - stars: massive KW - Magellanic Clouds KW - stars: winds, outflows KW - Hertzsprung-Russell and C-M diagrams KW - techniques: spectroscopic KW - ISM: bubbles Y1 - 2018 U6 - https://doi.org/10.1051/0004-6361/201832816 SN - 1432-0746 VL - 615 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Kubatova, Brankica A1 - Szecsi, D. A1 - Sander, Andreas Alexander Christoph A1 - Kubat, Jiří A1 - Tramper, F. A1 - Krticka, Jiri A1 - Kehrig, C. A1 - Hamann, Wolf-Rainer A1 - Hainich, Rainer A1 - Shenar, Tomer T1 - Low-metallicity massive single stars with rotation BT - II. Predicting spectra and spectral classes of chemically homogeneously evolving stars JF - Astronomy and astrophysics : an international weekly journal N2 - Context. Metal-poor massive stars are assumed to be progenitors of certain supernovae, gamma-ray bursts, and compact object mergers that might contribute to the early epochs of the Universe with their strong ionizing radiation. However, this assumption remains mainly theoretical because individual spectroscopic observations of such objects have rarely been carried out below the metallicity of the Small Magellanic Cloud. Aims. Here we explore the predictions of the state-of-the-art theories of stellar evolution combined with those of stellar atmospheres about a certain type of metal-poor (0.02 Z(circle dot)) hot massive stars, the chemically homogeneously evolving stars that we call Transparent Wind Ultraviolet INtense (TWUIN) stars. Methods. We computed synthetic spectra corresponding to a broad range in masses (20 130 M-circle dot) and covering several evolutionary phases from the zero-age main-sequence up to the core helium-burning stage. We investigated the influence of mass loss and wind clumping on spectral appearance and classified the spectra according to the Morgan-Keenan (MK) system. Results. We find that TWUIN stars show almost no emission lines during most of their core hydrogen-burning lifetimes. Most metal lines are completely absent, including nitrogen. During their core helium-burning stage, lines switch to emission, and even some metal lines (oxygen and carbon, but still almost no nitrogen) are detected. Mass loss and clumping play a significant role in line formation in later evolutionary phases, particularly during core helium-burning. Most of our spectra are classified as an early-O type giant or supergiant, and we find Wolf-Rayet stars of type WO in the core helium-burning phase. Conclusions. An extremely hot, early-O type star observed in a low-metallicity galaxy could be the result of chemically homogeneous evolution and might therefore be the progenitor of a long-duration gamma-ray burst or a type Ic supernova. TWUIN stars may play an important role in reionizing the Universe because they are hot without showing prominent emission lines during most of their lifetime. KW - stars: massive KW - stars: winds, outflows KW - stars: rotation KW - galaxies: dwarf KW - radiative transfer Y1 - 2019 U6 - https://doi.org/10.1051/0004-6361/201834360 SN - 1432-0746 SN - 0004-6361 VL - 623 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Hamann, Wolf-Rainer A1 - Gräfener, G. A1 - Liermann, A. A1 - Hainich, Rainer A1 - Sander, Andreas Alexander Christoph A1 - Shenar, Tomer A1 - Ramachandran, Varsha A1 - Todt, Helge Tobias A1 - Oskinova, Lida T1 - The Galactic WN stars revisited BT - Impact of Gaia distances on fundamental stellar parameters JF - Astronomy and astrophysics : an international weekly journal N2 - Comprehensive spectral analyses of the Galactic Wolf-Rayet stars of the nitrogen sequence (i.e. the WN subclass) have been performed in a previous paper. However, the distances of these objects were poorly known. Distances have a direct impact on the "absolute" parameters, such as luminosities and mass-loss rates. The recent Gaia Data Release (DR2) of trigonometric parallaxes includes nearly all WN stars of our Galactic sample. In the present paper, we apply the new distances to the previously analyzed Galactic WN stars and rescale the results accordingly. On this basis, we present a revised catalog of 55 Galactic WN stars with their stellar and wind parameters. The correlations between mass-loss rate and luminosity show a large scatter, for the hydrogen-free WN stars as well as for those with detectable hydrogen. The slopes of the log L - log M correlations are shallower than found previously. The empirical Hertzsprung-Russell diagram (HRD) still shows the previously established dichotomy between the hydrogen-free early WN subtypes that are located on the hot side of the zero-age main sequence (ZAMS), and the late WN subtypes, which show hydrogen and reside mostly at cooler temperatures than the ZAMS (with few exceptions). However, with the new distances, the distribution of stellar luminosities became more continuous than obtained previously. The hydrogen-showing stars of late WN subtype are still found to be typically more luminous than the hydrogen-free early subtypes, but there is a range of luminosities where both subclasses overlap. The empirical HRD of the Galactic single WN stars is compared with recent evolutionary tracks. Neither these single-star evolutionary models nor binary scenarios can provide a fully satisfactory explanation for the parameters of these objects and their location in the HRD. KW - stars: mass-loss KW - stars: winds, outflows KW - stars: Wolf-Rayet KW - stars: atmospheres KW - stars: evolution KW - stars: distances Y1 - 2019 U6 - https://doi.org/10.1051/0004-6361/201834850 SN - 1432-0746 VL - 625 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Munoz, Melissa A1 - Moffat, Anthony F. J. A1 - Hill, Grant M. A1 - Shenar, Tomer A1 - Richardson, Noel D. A1 - Pablo, Herbert A1 - St-Louis, Nicole A1 - Ramiaramanantsoa, Tahina T1 - WR 148: identifying the companion of an extreme runaway massive binary JF - Monthly notices of the Royal Astronomical Society N2 - WR 148 (HD 197406) is an extreme runaway system considered to be a potential candidate for a short-period (4.3173 d) rare WR + compact object binary. Provided with new high-resolution, high signal-to-noise spectra from the Keck observatory, we determine the orbital parameters for both the primary WR and the secondary, yielding respective projected orbital velocity amplitudes of 88.1 ± 3.8 km s−1 and 79.2 ± 3.1 km s−1 and implying a mass ratio of 1.1 ± 0.1. We then apply the shift-and-add technique to disentangle the spectra and obtain spectra compatible with a WN7ha and an O4-6 star. Considering an orbital inclination of ∼67°, derived from previous polarimetry observations, the system's total mass would be a mere 2–3M⊙⁠, an unprecedented result for a putative massive binary system. However, a system comprising a 37M⊙ secondary (typical mass of an O5V star) and a 33M⊙ primary (given the mass ratio) would infer an inclination of ∼18°. We therefore reconsider the previous methods of deriving the orbital inclination based on time-dependent polarimetry and photometry. While the polarimetric results are inconclusive requiring better data, the photometric results favour low inclinations. Finally, we compute WR 148’s space velocity and retrace the runaway's trajectory back to the Galactic plane (GP). With an ejection velocity of 198 ± 27 km s−1 and a travel time of 4.7 ± 0.8 Myr to reach its current location, WR 148 was most likely ejected via dynamical interactions in a young cluster. KW - binaries: spectroscopic KW - stars: individual: WR 148 KW - stars: kinematics and dynamics KW - stars: mass-loss KW - stars: winds, outflows KW - stars: Wolf-Rayet Y1 - 2017 U6 - https://doi.org/10.1093/mnras/stw2283 SN - 0035-8711 SN - 1365-2966 VL - 467 SP - 3105 EP - 3121 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Sander, Andreas Alexander Christoph A1 - Hamann, Wolf-Rainer A1 - Todt, Helge Tobias A1 - Hainich, Rainer A1 - Shenar, Tomer T1 - Coupling hydrodynamics with comoving frame radiative transfer I. A unified approach for OB and WR stars JF - Astronomy and astrophysics : an international weekly journal N2 - Context. For more than two decades, stellar atmosphere codes have been used to derive the stellar and wind parameters of massive stars. Although they have become a powerful tool and sufficiently reproduce the observed spectral appearance, they can hardly be used for more than measuring parameters. One major obstacle is their inconsistency between the calculated radiation field and the wind stratification due to the usage of prescribed mass-loss rates and wind-velocity fields. Aims. We present the concepts for a new generation of hydrodynamically consistent non-local thermodynamical equilibrium (nonLTE) stellar atmosphere models that allow for detailed studies of radiation-driven stellar winds. As a first demonstration, this new kind of model is applied to a massive O star. Methods. Based on earlier works, the PoWR code has been extended with the option to consistently solve the hydrodynamic equation together with the statistical equations and the radiative transfer in order to obtain a hydrodynamically consistent atmosphere stratification. In these models, the whole velocity field is iteratively updated together with an adjustment of the mass-loss rate. Results. The concepts for obtaining hydrodynamically consistent models using a comoving-frame radiative transfer are outlined. To provide a useful benchmark, we present a demonstration model, which was motivated to describe the well-studied O4 supergiant zeta Pup. The obtained stellar and wind parameters are within the current range of literature values. Conclusions. For the first time, the PoWR code has been used to obtain a hydrodynamically consistent model for a massive O star. This has been achieved by a profound revision of earlier concepts used for Wolf-Rayet stars. The velocity field is shaped by various elements contributing to the radiative acceleration, especially in the outer wind. The results further indicate that for more dense winds deviations from a standard beta-law occur. KW - stars: mass-loss KW - stars: winds, outflows KW - stars: early-type KW - stars: atmospheres KW - stars: fundamental parameters KW - stars: massive Y1 - 2017 U6 - https://doi.org/10.1051/0004-6361/201730642 SN - 1432-0746 VL - 603 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Oskinova, Lida A1 - Huenemoerder, D. P. A1 - Hamann, Wolf-Rainer A1 - Shenar, Tomer A1 - Sander, Andreas Alexander Christoph A1 - Ignace, R. A1 - Todt, Helge Tobias A1 - Hainich, Rainer T1 - On the Binary Nature of Massive Blue Hypergiants: High-resolution X-Ray Spectroscopy Suggests That Cyg OB2 12 is a Colliding Wind Binary JF - The astrophysical journal : an international review of spectroscopy and astronomical physics N2 - The blue hypergiant Cyg OB2 12 (B3Ia(+)) is a representative member of the class of very massive stars in a poorly understood evolutionary stage. We obtained its high-resolution X-ray spectrum using the Chandra observatory. PoWR model atmospheres were calculated to provide realistic wind opacities and to establish the wind density structure. We find that collisional de-excitation is the dominant mechanism depopulating the metastable upper levels of the forbidden lines of the He-like ions Si XIV and Mg XII. Comparison between the model and observations reveals that X-ray emission is produced in a dense plasma, which could reside only at the photosphere or in a colliding wind zone between binary components. The observed X-ray spectra are well-fitted by thermal plasma models, with average temperatures in excess of 10 MK. The wind speed in Cyg OB2 12 is not high enough to power such high temperatures, but the collision of two winds in a binary system can be sufficient. We used archival data to investigate the X-ray properties of other blue hypergiants. In general, stars of this class are not detected as X-ray sources. We suggest that our new Chandra observations of Cyg OB2 12 can be best explained if Cyg OB2 12 is a colliding wind binary possessing a late O-type companion. This makes Cyg OB2 12 only the second binary system among the 16 known Galactic hypergiants. This low binary fraction indicates that the blue hypergiants are likely products of massive binary evolution during which they either accreted a significant amount of mass or already merged with their companions. KW - stars: individual (Cyg OB2 12) KW - stars: massive KW - stars: mass-loss KW - stars: winds, outflows KW - supergiants KW - X-rays: stars Y1 - 2017 U6 - https://doi.org/10.3847/1538-4357/aa7e79 SN - 0004-637X SN - 1538-4357 VL - 845 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Gimenez-Garcia, Ana A1 - Shenar, Tomer A1 - Torrejon, J. M. A1 - Oskinova, Lida A1 - Martinez-Nunez, S. A1 - Hamann, Wolf-Rainer A1 - Rodes-Roca, J. J. A1 - González-Galan, A. A1 - Alonso-Santiago, J. A1 - González-Fernández, C. A1 - Bernabeu, Guillermo A1 - Sander, Andreas Alexander Christoph T1 - Measuring the stellar wind parameters in IGR J17544-2619 and Vela X-1 constrains the accretion physics in supergiant fast X-ray transient and classical supergiant X-ray binaries JF - Siberian Mathematical Journal N2 - Aims. To close this gap, we perform a comparative analysis of the optical companion in two important systems: IGR J175442619 (SFXT) and Vela X-1 (SGXB). We analyze the spectra of each star in detail and derive their stellar and wind properties. As a next step, we compare the wind parameters, giving us an excellent chance of recognizing key differences between donor winds in SFXTs and SGXBs. Methods. We use archival infrared, optical and ultraviolet observations, and analyze them with the non-local thermodynamic equilibrium (NLTE) Potsdam Wolf-Rayet model atmosphere code. We derive the physical properties of the stars and their stellar winds, accounting for the influence of X-rays on the stellar winds. Results. We find that the stellar parameters derived from the analysis generally agree well with the spectral types of the two donors: O9I (IGR J17544-2619) and B0.5Iae (Vela X-1). The distance to the sources have been revised and also agree well with the estimations already available in the literature. In IGR J17544-2619 we are able to narrow the uncertainty to d = 3.0 +/- 0.2 kpc. From the stellar radius of the donor and its X-ray behavior, the eccentricity of IGR J17544-2619 is constrained to e < 0.25. The derived chemical abundances point to certain mixing during the lifetime of the donors. An important difference between the stellar winds of the two stars is their terminal velocities (v(infinity) = 1500 km s(-1) in IGR J17544-2619 and v(infinity) = 700 km s(-1) in Vela X-1), which have important consequences on the X-ray luminosity of these sources. Conclusions. The donors of IGR J17544-2619 and Vela X-1 have similar spectral types as well as similar parameters that physically characterize them and their spectra. In addition, the orbital parameters of the systems are similar too, with a nearly circular orbit and short orbital period. However, they show moderate differences in their stellar wind velocity and the spin period of their neutron star which has a strong impact on the X-ray luminosity of the sources. This specific combination of wind speed and pulsar spin favors an accretion regime with a persistently high luminosity in Vela X-1, while it favors an inhibiting accretion mechanism in IGR J17544-2619. Our study demonstrates that the relative wind velocity is critical in class determination for the HMXBs hosting a supergiant donor, given that it may shift the accretion mechanism from direct accretion to propeller regimes when combined with other parameters. KW - accretion, accretion disks KW - methods: observational KW - techniques: spectroscopic KW - stars: atmospheres KW - X-rays: binaries KW - stars: winds, outflows Y1 - 2016 U6 - https://doi.org/10.1051/0004-6361/201527551 SN - 1432-0746 VL - 591 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Richardson, Noel D. A1 - Shenar, Tomer A1 - Roy-Loubier, Olivier A1 - Schaefer, Gail A1 - Moffat, Anthony F. J. A1 - St-Louis, Nicole A1 - Gies, Douglas R. A1 - Farrington, Chris A1 - Hill, Grant M. A1 - Williams, Peredur M. A1 - Gordon, Kathryn A1 - Pablo, Herbert A1 - Ramiaramanantsoa, Tahina T1 - The CHARA Array resolves the long-period Wolf-Rayet binaries WR 137 and WR 138 JF - Monthly notices of the Royal Astronomical Society N2 - We report on interferometric observations with the CHARAArray of two classical Wolf-Rayet (WR) stars in suspected binary systems, namely WR 137 and WR 138. In both cases, we resolve the component stars to be separated by a few milliarcseconds. The data were collected in the H band, and provide a measure of the fractional flux for both stars in each system. We find that the WR star is the dominant H-band light source in both systems (fWR, 137 = 0.59 +/- 0.04; fWR, 138 = 0.67 +/- 0.01), which is confirmed through both comparisons with estimated fundamental parameters for WR stars and O dwarfs, as well as through spectral modelling of each system. Our spectral modelling also provides fundamental parameters for the stars and winds in these systems. The results on WR 138 provide evidence that it is a binary system which may have gone through a previous mass-transfer episode to create the WR star. The separation and position of the stars in the WR 137 system together with previous results from the IOTA interferometer provides evidence that the binary is seen nearly edgeon. The possible edge-on orbit of WR 137 aligns well with the dust production site imaged by the Hubble Space Telescope during a previous periastron passage, showing that the dust production may be concentrated in the orbital plane. KW - binaries: visual KW - stars: individual: WR 137 KW - stars: individual: WR 138 KW - stars: mass-loss KW - stars: winds, outflows KW - stars: Wolf-Rayet Y1 - 2016 U6 - https://doi.org/10.1093/mnras/stw1585 SN - 0035-8711 SN - 1365-2966 VL - 461 SP - 4115 EP - 4124 PB - Oxford Univ. Press CY - Oxford ER -