Showing 1–50 of 68 results for author: Howarth, I D

Rapid polarization variations in the O4 supergiant $ζ$ Puppis

Authors: Jeremy Bailey, Ian D. Howarth, Daniel V. Cotton, Lucyna Kedziora-Chudczer, Ain De Horta, Sarah L. Martell, Colin Eldridge, Paul Luckas

Abstract: We present time-series linear-polarization observations of the bright O4 supergiant $ζ$ Puppis. The star is found to show polarization variation on timescales of around an hour and longer. Many of the observations were obtained contemporaneously with Transiting Exoplanet Survey Satellite (TESS) photometry. We find that the polarization varies on similar timescales to those seen in the TESS light-c… ▽ More We present time-series linear-polarization observations of the bright O4 supergiant $ζ$ Puppis. The star is found to show polarization variation on timescales of around an hour and longer. Many of the observations were obtained contemporaneously with Transiting Exoplanet Survey Satellite (TESS) photometry. We find that the polarization varies on similar timescales to those seen in the TESS light-curve. The previously reported 1.78-day photometric periodicity is seen in both the TESS and polarization data. The amplitude ratio of photometry to polarization is ~9 for the periodic component and the polarization variation is oriented along position angle ~70 deg-160 deg. Higher-frequency stochastic variability is also seen in both datasets with an amplitude ratio of ~19 and no preferred direction. We model the polarization expected for a rotating star with bright photospheric spots and find that models that fit the photometric variation produce too little polarization variation to explain the observations. We suggest that the variable polarization is more likely the result of scattering from the wind, with corotating interaction regions producing the periodic variation and a clumpy outflow producing the stochastic component. The H$α$ emission line strength was seen to increase by 10% in 2021 with subsequent observations showing a return to the pre-2018 level. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 19 pages, 16 figures, Accepted by MNRAS

Discovery of Dipolar Chromospheres in Two White Dwarfs

Authors: J. Farihi, J. J. Hermes, S. P. Littlefair, I. D. Howarth, N. Walters, S. G. Parsons

Abstract: This paper reports the ULTRACAM discovery of dipolar surface spots in two cool magnetic white dwarfs with Balmer emission lines, while a third system exhibits a single spot, similar to the prototype GD 356. The light curves are modeled with simple, circular, isothermal dark spots, yielding relatively large regions with minimum angular radii of 20 deg. For those stars with two light curve minima, t… ▽ More This paper reports the ULTRACAM discovery of dipolar surface spots in two cool magnetic white dwarfs with Balmer emission lines, while a third system exhibits a single spot, similar to the prototype GD 356. The light curves are modeled with simple, circular, isothermal dark spots, yielding relatively large regions with minimum angular radii of 20 deg. For those stars with two light curve minima, the dual spots are likely observed at high inclination (or colatitude), however, identical and antipodal spots cannot simultaneously reproduce both the distinct minima depths and the phases of the light curve maxima. The amplitudes of the multi-band photometric variability reported here are all several times larger than that observed in the prototype GD 356; nevertheless, all DAHe stars with available data appear to have light curve amplitudes that increase toward the blue in correlated ratios. This behavior is consistent with cool spots that produce higher contrasts at shorter wavelengths, with remarkably similar spectral properties given the diversity of magnetic field strengths and rotation rates. These findings support the interpretation that some magnetic white dwarfs generate intrinsic chromospheres as they cool, and that no external source is responsible for the observed temperature inversion. Spectroscopic time-series data for DAHe stars is paramount for further characterization, where it is important to obtain well-sampled data, and consider wavelength shifts, equivalent widths, and spectropolarimetry. △ Less

Submitted 18 July, 2023; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: 9 pages, 5 figures, 3 tables, accepted to MNRAS

A study of the rapid rotator $ζ$ Aql: differential surface rotation?

Authors: Ian D. Howarth, Jeremy Bailey, Daniel V. Cotton, Lucyna Kedziora-Chudczer

Abstract: We report new, extremely precise, photopolarimetry of the rapidly-rotating A0 main-sequence star $ζ$ Aql, covering the wavelength range $\sim$400--900nm, which reveals a rotationally-induced signal. We model the polarimetry, together with the flux distribution and line profiles, in the framework of Roche geometry with $ω$-model gravity darkening, to establish the stellar parameters. An additional… ▽ More We report new, extremely precise, photopolarimetry of the rapidly-rotating A0 main-sequence star $ζ$ Aql, covering the wavelength range $\sim$400--900nm, which reveals a rotationally-induced signal. We model the polarimetry, together with the flux distribution and line profiles, in the framework of Roche geometry with $ω$-model gravity darkening, to establish the stellar parameters. An additional constraint is provided by TESS photometry, which shows variability with a period, $P_{\rm phot}$, of 11.1 hr. Modelling based on solid-body surface rotation gives rotation periods, $P_{\rm rot}$, that are in only marginal agreement with this value. We compute new ESTER stellar-structure models to predict horizontal surface velocity fields, which depart from solid-body rotation at only the $\sim$2% level (consistent with a reasonably strong empirical upper limit on differential rotation derived from the line-profile analysis). These models bring the equatorial rotation period, $P_{\rm rot,e}$, into agreement with $P_{\rm phot}$, without requiring any 'fine tuning' (for the Gaia parallax). We confirm that surface abundances are significantly subsolar ($\mbox{[M/H]} \simeq -0.5$). The star's basic parameters are established with reasonably good precision: ${M = 2.53\pm0.16\,\mbox{M}_\odot}$, $\log{L/\mbox{L}_\odot} = 1.72\pm0.02$, $R_{\rm p} = 2.21\pm 0.02\,\mbox{R}_\odot$, $T_{\rm eff} = 9693 \pm 50~\mbox{K}$, $i = 85{^{+5}_{-7}}^\circ$, and $ω/ω_{\rm c} = 0.95\pm0.02$. Comparison with single-star, solar-abundance stellar-evolution models incorporating rotational effects shows excellent agreement (but somewhat poorer agreement for models at $\mbox{[M/H]} \simeq -0.4$). △ Less

Submitted 28 January, 2023; v1 submitted 12 January, 2023; originally announced January 2023.

Comments: MNRAS, accepted Minor typographical edits in v3

A study of the F-giant star $θ$ Scorpii A: a post-merger rapid rotator?

Authors: Fiona Lewis, Jeremy Bailey, Daniel V. Cotton, Ian D. Howarth, Lucyna Kedziora-Chudczer, Floor van Leeuwen

Abstract: We report high-precision observations of the linear polarization of the F1$\,$III star $θ$ Scorpii. The polarization has a wavelength dependence of the form expected for a rapid rotator, but with an amplitude several times larger than seen in otherwise similar main-sequence stars. This confirms the expectation that lower-gravity stars should have stronger rotational-polarization signatures as a co… ▽ More We report high-precision observations of the linear polarization of the F1$\,$III star $θ$ Scorpii. The polarization has a wavelength dependence of the form expected for a rapid rotator, but with an amplitude several times larger than seen in otherwise similar main-sequence stars. This confirms the expectation that lower-gravity stars should have stronger rotational-polarization signatures as a consequence of the density dependence of the ratio of scattering to absorption opacities. By modelling the polarization, together with additional observational constraints (incorporating a revised analysis of Hipparcos astrometry, which clarifies the system's binary status), we determine a set of precise stellar parameters, including a rotation rate $ω\, (= Ω/Ω_{\rm c})\ge 0.94$, polar gravity $\log{g_p} = 2.091 ^{+0.042}_{-0.039}$ (dex cgs), mass $3.10 ^{+0.37}_{-0.32}$ solar masses, and luminosity $\log(L/Lsun) =3.149^{+0.041}_{-0.028}$. These values are incompatible with evolutionary models of single rotating stars, with the star rotating too rapidly for its evolutionary stage, and being undermassive for its luminosity. We conclude that $θ$ Sco A is most probably the product of a binary merger. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: 12 pages, 10 figures, Accepted by MNRAS

Properties of the Be-type stars in 30 Doradus

Authors: P. L. Dufton, D. J. Lennon, J. I. Villasenor, I. D. Howarth, C. J. Evans, S. E. de Mink, H. Sana, W. D. Taylor

Abstract: The evolutionary status of Be-type stars remains unclear, with both single-star and binary pathways having been proposed. Here, VFTS spectroscopy of 73 Be-type stars, in the spectral-type range, B0--B3, is analysed to estimate projected rotational velocities, radial velocities and stellar parameters. They are found to be rotating faster than the corresponding VFTS B-type sample but simulations imp… ▽ More The evolutionary status of Be-type stars remains unclear, with both single-star and binary pathways having been proposed. Here, VFTS spectroscopy of 73 Be-type stars, in the spectral-type range, B0--B3, is analysed to estimate projected rotational velocities, radial velocities and stellar parameters. They are found to be rotating faster than the corresponding VFTS B-type sample but simulations imply that their projected rotational velocities are inconsistent with them all rotating at near critical velocities. The de-convolution of the projected rotational velocities estimates leads to a mean rotational velocity estimate of 320-350 km/s, approximately 100 km/s larger than that for the corresponding B-type sample. There is a dearth of targets with rotational velocities less than 0.4 of the critical velocity, with a broad distribution reaching up to critical rotation. Our best estimate for the mean or median of the rotational velocitiy is 0.68 of the critical velocity. Rapidly-rotating B-type stars are more numerous than their Be-type counterparts, whilst the observed frequency of Be-type stars identified as binary systems is significantly lower than that for normal B-type stars, consistent with their respective radial-velocity dispersions. The semi-amplitudes for the Be-type binaries are also smaller. Similar results are found for a SMC Be-type sample centred on NGC346 with no significant differences being found between the two samples. These results are compared with the predictions of single and binary stellar evolutionary models for Be-type stars. Assuming that a single mechanism dominated the production of classical Be-type stars, our comparison would favour a binary evolutionary history. △ Less

Submitted 7 March, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS, Appendicesand Tables 2 and 4 published as supplementary material

Carbon-enhanced stars with short orbital and spin periods

Authors: L. J. Whitehouse, J. Farihi, I. D. Howarth, S. Mancino, N. Walters, A. Swan, T. G. Wilson, J. Guo

Abstract: Many characteristics of dwarf carbon stars are broadly consistent with a binary origin, including mass transfer from an evolved companion. While the population overall appears to have old-disc or halo kinematics, roughly 2$\,$per cent of these stars exhibit H$α$ emission, which in low-mass main-sequence stars is generally associated with rotation and relative youth. Its presence in an older popula… ▽ More Many characteristics of dwarf carbon stars are broadly consistent with a binary origin, including mass transfer from an evolved companion. While the population overall appears to have old-disc or halo kinematics, roughly 2$\,$per cent of these stars exhibit H$α$ emission, which in low-mass main-sequence stars is generally associated with rotation and relative youth. Its presence in an older population therefore suggests either irradiation or spin-up. This study presents time-series analyses of photometric and radial-velocity data for seven dwarf carbon stars with H$α$ emission. All are shown to have photometric periods in the range 0.2--5.2$\,$d, and orbital periods of similar length, consistent with tidal synchronisation. It is hypothesised that dwarf carbon stars with emission lines are the result of close-binary evolution, indicating that low-mass, metal-weak or metal-poor stars can accrete substantial material prior to entering a common-envelope phase. △ Less

Submitted 3 August, 2021; v1 submitted 30 April, 2021; originally announced May 2021.

Comments: 16 pages, 7 figures, Accepted to MNRAS

The Tarantula Massive Binary Monitoring V. R 144: a wind-eclipsing binary with a total mass > 140 Msun

Authors: T. Shenar, H. Sana, P. Marchant, B. Pablo, N. Richardson, A. F. J. Moffat, T. Van Reeth, R. H. Barba, D. M. Bowman, P. Broos, P. A. Crowther, J. S. Clark, A. de Koter, S. E. de Mink, K. Dsilva, G. Graefener, I. D. Howarth, N. Langer, L. Mahy, J. Maiz Apellaniz, A. M. T. Pollock, F. R. N. Schneider, L. Townsley, J. S. Vink

Abstract: R 144 is the visually brightest WR star in the Large Magellanic Cloud (LMC). R 144 was reported to be a binary, making it potentially the most massive binary thus observed. We perform a comprehensive spectral, photometric, orbital, and polarimetric analysis of R 144. R 144 is an eccentric (e=0.51) 74.2-d binary comprising two relatively evolved (age~2 Myr), H-rich WR stars. The hotter primary (W… ▽ More R 144 is the visually brightest WR star in the Large Magellanic Cloud (LMC). R 144 was reported to be a binary, making it potentially the most massive binary thus observed. We perform a comprehensive spectral, photometric, orbital, and polarimetric analysis of R 144. R 144 is an eccentric (e=0.51) 74.2-d binary comprising two relatively evolved (age~2 Myr), H-rich WR stars. The hotter primary (WN5/6h, T=50 kK) and the cooler secondary (WN6/7h,T=45kK) have nearly equal masses. The combination of low rotation and H-depletion observed in the system is well reproduced by contemporary evolution models that include boosted mass-loss at the upper-mass end. The systemic velocity of R 144 and its relative isolation suggest that it was ejected as a runaway from the neighbouring R 136 cluster. The optical light-curve shows a clear orbital modulation that can be well explained as a combination of two processes: excess emission stemming from wind-wind collisions and double wind eclipses. Our light-curve model implies an orbital inclination of i=60.4+-1.5deg, resulting in accurately constrained dynamical masses of 74+-4 and 69+-4 Msun. Assuming that both binary components are core H-burning, these masses are difficult to reconcile with the derived luminosities (logL1,2 = 6.44, 6.39 [Lsun]), which correspond to evolutionary masses of the order of 110 and 100Msun, respectively. Taken at face value, our results imply that both stars have high classical Eddington factors of Gamma_e = 0.78+-0.1. If the stars are on the main sequence, their derived radii (~25Rsun) suggest that they are only slightly inflated, even at this high Eddington factor. Alternatively, the stars could be core-He burning, strongly inflated from the regular size of classical Wolf-Rayet stars (~1Rsun), a scenario that could help resolve the observed mass discrepancy. △ Less

Submitted 28 October, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: Accepted to A&A on the April 06 2021. 17 pages + 5 page appendix. Version 28.10.2021: corrected typo in Eq. (2); numerator and denominator were inversed on RHS (only typo; no impact on results)

Journal ref: A&A 650, A147 (2021)

The rotation of alpha Oph investigated using polarimetry

Authors: Jeremy Bailey, Daniel V. Cotton, Ian D. Howarth, Fiona Lewis, Lucyna Kedziora-Chudczer

Abstract: Recently we have demonstrated that high-precision polarization observations can detect the polarization resulting from the rotational distortion of a rapidly rotating B-type star. Here we investigate the extension of this approach to an A-type star. Linear-polarization observations of $α$ Oph (A5IV) have been obtained over wavelengths from 400 to 750 nm. They show the wavelength dependence expecte… ▽ More Recently we have demonstrated that high-precision polarization observations can detect the polarization resulting from the rotational distortion of a rapidly rotating B-type star. Here we investigate the extension of this approach to an A-type star. Linear-polarization observations of $α$ Oph (A5IV) have been obtained over wavelengths from 400 to 750 nm. They show the wavelength dependence expected for a rapidly-rotating star combined with a contribution from interstellar polarization. We model the observations by fitting rotating-star polarization models and adding additional constraints including a measured $v_e \sin{i}$. However, we cannot fully separate the effects of rotation rate and inclination, leaving a range of possible solutions. We determine a rotation rate $ω= Ω/Ω_ c$ between 0.83 and 0.98 and an axial inclination i > 60 deg. The rotation-axis position angle is found to be 142 $\pm$ 4 deg, differing by 16 deg from a value obtained by interferometry. This might be due to precession of the rotation axis due to interaction with the binary companion. Other parameters resulting from the analysis include a polar temperature Tp = 8725 $\pm$ 175 K, polar gravity $\log{g_p} = 3.93 \pm 0.08$ (dex cgs), and polar radius $R_{\rm p} = 2.52 \pm 0.06$ Rsun. Comparison with rotating-star evolutionary models indicates that $α$ Oph is in the later half of its main-sequence evolution and must have had an initial $ω$ of 0.8 or greater. The interstellar polarization has a maximum value at a wavelength ($λ_{\rm max}$) of $440 \pm 110$ nm, consistent with values found for other nearby stars. △ Less

Submitted 17 March, 2020; originally announced March 2020.

Comments: 14 pages, 11 figures, 5 tables, Accepted in MNRAS

A search for strong magnetic fields in massive and very massive stars in the Magellanic Clouds

Authors: S. Bagnulo, G. A. Wade, Y. Naze, J. H. Grunhut, M. E. Shultz, D. J. Asher, P. A. Crowther, C. J. Evans, A. David-Uraz, I. D. Howarth, N. Morrell, M. S. Munoz, C. Neiner, J. Puls, M. K. Szymanski, J. S. Vink

Abstract: Despite their rarity, massive stars dominate the ecology of galaxies via their strong, radiatively-driven winds throughout their lives and as supernovae in their deaths. However, their evolution and subsequent impact on their environment can be significantly affected by the presence of a magnetic field. While recent studies indicate that about 7% of OB stars in the Milky Way host strong, stable, o… ▽ More Despite their rarity, massive stars dominate the ecology of galaxies via their strong, radiatively-driven winds throughout their lives and as supernovae in their deaths. However, their evolution and subsequent impact on their environment can be significantly affected by the presence of a magnetic field. While recent studies indicate that about 7% of OB stars in the Milky Way host strong, stable, organised (fossil) magnetic fields at their surfaces, little is known about the fields of very massive stars, nor the magnetic properties of stars outside our Galaxy. We aim to continue searching for strong magnetic fields in a diverse set of massive and very massive stars (VMS) in the Large and Small Magellanic Clouds (LMC/SMC), and we evaluate the overall capability of FORS2 to usefully search for and detect stellar magnetic fields in extra-galactic environments. We have obtained FORS2 spectropolarimetry of a sample of 41 stars, which principally consist of spectral types B, O, Of/WN, WNh, and classical WR stars in the LMC and SMC. Four of our targets are Of?p stars; one of them was just recently discovered. Each spectrum was analysed to infer the longitudinal magnetic field. No magnetic fields were formally detected in our study, although Bayesian statistical considerations suggest that the Of?p star SMC159-2 is magnetic with a dipolar field of the order of 2.4 to 4.4kG. In addition, our first constraints of magnetic fields in VMS provide interesting insights into the formation of the most massive stars in the Universe. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Comments: Accepted by A&A

Journal ref: A&A, 635, A163 (15p, 2020)

Two New Rapidly-Rotating ON Stars Found With LAMOST

Authors: Guang-Wei Li, Ian D. Howarth

Abstract: The ON stars are a rare subtype of O stars, of uncertain origin. We report two new, rapidly-rotating ON stars found in data acquired with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, LAMOST. LS I +61 28 is an ON8.5 Vn dwarf with a projected equatorial rotational velocity of $v_{\textrm e}\sin{i} \simeq 298$ km$\;$s$^{-1}$, while HDE 236672 is an ON9 IVn subgiant with… ▽ More The ON stars are a rare subtype of O stars, of uncertain origin. We report two new, rapidly-rotating ON stars found in data acquired with the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, LAMOST. LS I +61 28 is an ON8.5 Vn dwarf with a projected equatorial rotational velocity of $v_{\textrm e}\sin{i} \simeq 298$ km$\;$s$^{-1}$, while HDE 236672 is an ON9 IVn subgiant with $v_{\textrm e}\sin{i} \simeq 253$ km$\;$s$^{-1}$ The former is the first rapidly-rotating ON dwarf to be found, and the latter is only the third ON subgiant. The luminosity classes of non-supergiant ON stars appear to be influenced by the axial inclination angle $i$: the rapidly-rotating giants are close to equator-on, while ON dwarfs with lower $v_{\textrm e}\sin{i}$ values are viewed more nearly pole-on. Combining parallaxes and proper motions from Gaia DR2 with radial-velocity measurements, we investigate the kinematics of non-supergiant ON stars, and infer that the dynamics, rapid rotation, and surface-nitrogen characteristics may all be consequences of binary interaction. △ Less

Submitted 16 November, 2019; originally announced November 2019.

Comments: 8 pages, 4 figures, 2 tables

Journal ref: The Astrophysical Journal, 888:81 (6pp), 2020 January 10

The distance, rotation, and physical parameters of zeta Pup

Authors: Ian D. Howarth, Floor van Leeuwen

Abstract: We scrutinize the Hipparcos parallax for the bright O supergiant zeta Pup, and confirm that the implied distance of 332+/-11pc appears to be reliable. We then review the implications for the star's physical parameters, and the consequences for the interpretation of P(phot), the 1.78-d photometric period. The equatorial rotation period is <3.7d (with 95% confidence), ruling out a proposed ~5.1d val… ▽ More We scrutinize the Hipparcos parallax for the bright O supergiant zeta Pup, and confirm that the implied distance of 332+/-11pc appears to be reliable. We then review the implications for the star's physical parameters, and the consequences for the interpretation of P(phot), the 1.78-d photometric period. The equatorial rotation period is <3.7d (with 95% confidence), ruling out a proposed ~5.1d value. If the photometric period is the rotation period then i, the inclination of the rotation axis to the line of sight, is 33.2+/-1.8 degrees. The inferred mass, radius, and luminosity are securely established to be less than canonical values for the spectral type, and are not in agreement with single-star evolution models. The runaway status, rapid rotation, and anomalous physical properties are all indicative of an evolutionary history involving binary (or multiple-star) interaction. We perform simple starspot modelling to show that the low axial inclination required if P(rot) = 1.78d has testable spectroscopic consequences, which have not been identified in existing time series. If P(phot) is directly related to drivers of systematic, high-velocity stellar-wind variability (`discrete absorption components') in zeta Pup, antisolar differential rotation is required. Model line profiles calculated on that basis are at variance with observations. △ Less

Submitted 23 January, 2019; originally announced January 2019.

Comments: MNRAS, submitted

A remarkable change of the spectrum of the magnetic Of?p star HD 148937 reveals evidence of an eccentric, high-mass binary

Authors: G. A. Wade, J. V. Smoker, C. J. Evans, I. D. Howarth, R. Barba, N. L. J. Cox, N. Morrell, Y. Nazé, J. Cami, A. Farhang, N. R. Walborn, J. Arias, R. Gamen

Abstract: We report new spectroscopic observations of the magnetic Of?p star HD 148937 obtained since 2015 that differ qualitatively from its extensive historical record of weak, periodic spectral variations. This remarkable behaviour represents clear evidence for an unprecedented change in the character of variability of the star. In this paper we describe the new spectral properties and compare them to th… ▽ More We report new spectroscopic observations of the magnetic Of?p star HD 148937 obtained since 2015 that differ qualitatively from its extensive historical record of weak, periodic spectral variations. This remarkable behaviour represents clear evidence for an unprecedented change in the character of variability of the star. In this paper we describe the new spectral properties and compare them to the previous line profiles. Based on measurements of the radial velocities of the C iii/N iii} emission lines near 4640 Å and the C iv absorption lines near 5800 Å, we infer that HD 148937 is likely a high-mass, double-lined spectroscopic binary. Combining the spectroscopic orbit with an archival interferometric measurement of the apparent separation of the equal-brightness components, we tentatively conclude that HD 148937 consists of two O-type stars with masses of approximately 34 and $49~M_\odot$, orbiting in an eccentric ($e=0.75$), long-period ($P_{\rm orb}\sim 26$ y) orbit. We discuss the potential relationship of the binary system to the peculiar properties of HD 148937, and propose future observations to refine the orbital and stellar properties. △ Less

Submitted 27 November, 2018; originally announced November 2018.

Comments: 12 pages, accepted by MNRAS

Journal ref: MNRAS, 483, 2581-2591 (2019)

Are the O stars in WR+O binaries exceptionally rapid rotators?

Authors: Dominic Reeve, Ian D. Howarth

Abstract: We examine claims of strong gravity-darkening effects in the O-star components of WR+O binaries. We generate synthetic spectra for a wide range of parameters, and show that the line-width results are consistent with extensive measurements of O stars that are either single or are members of `normal' binaries. By contrast, the WR+O results are at the extremes of, or outside, the distributions of bot… ▽ More We examine claims of strong gravity-darkening effects in the O-star components of WR+O binaries. We generate synthetic spectra for a wide range of parameters, and show that the line-width results are consistent with extensive measurements of O stars that are either single or are members of `normal' binaries. By contrast, the WR+O results are at the extremes of, or outside, the distributions of both models and other observations. Remeasurement of the WR+O spectra shows that they can be reconciled with other results by judicious choice of pseudo-continuum normalization. With this interpretation, the supersynchronous rotation previously noted for the O-star components in the WR+O binaries with the longest orbital periods appears to be unexceptional. Our investigation is therefore consistent with the aphorism that if the title of a paper ends with a question mark, the answer is probably `no'. △ Less

Submitted 9 May, 2018; originally announced May 2018.

Comments: Accepted in MNRAS

Polarization due to rotational distortion in the bright star Regulus

Authors: Daniel V. Cotton, Jeremy Bailey, Ian D. Howarth, Kimberly Bott, Lucyna Kedziora-Chudczer, P. W. Lucas, J. H. Hough

Abstract: Polarization in stars was first predicted by Chandrasekhar [1] who calculated a substantial linear polarization at the stellar limb for a pure electron-scattering atmosphere. This polarization will average to zero when integrated over a spherical star but could be detected if the symmetry is broken, for example by the eclipse of a binary companion. Nearly 50 years ago, Harrington and Collins [2] m… ▽ More Polarization in stars was first predicted by Chandrasekhar [1] who calculated a substantial linear polarization at the stellar limb for a pure electron-scattering atmosphere. This polarization will average to zero when integrated over a spherical star but could be detected if the symmetry is broken, for example by the eclipse of a binary companion. Nearly 50 years ago, Harrington and Collins [2] modeled another way of breaking the symmetry and producing net polarization - the distortion of a rapidly rotating hot star. Here we report the first detection of this effect. Observations of the linear polarization of Regulus, with two different high-precision polarimeters, range from +42 parts-per-million (ppm) at a wavelength of 741 nm to -22 ppm at 395 nm. The reversal from red to blue is a distinctive feature of rotation-induced polarization. Using a new set of models for the polarization of rapidly rotating stars we find that Regulus is rotating at 96.5$\substack{+0.6-0.8}$% of its critical angular velocity for breakup, and has an inclination greater than 76.5 degrees. The rotation axis of the star is at a position angle of 79.5$\pm$0.7 degrees. The conclusions are independent of, but in good agreement with, the results of previously published interferometric observations of Regulus [3]. The accurate measurement of rotation in early-type stars is important for understanding their stellar environments [4], and course of their evolution [5]. △ Less

Submitted 18 April, 2018; originally announced April 2018.

Comments: This is the full published article including the Methods and Supplementary Information sections: 33 pages, 10 figures, 8 tables

Journal ref: Nature Astronomy 1 (2017) 690-696

A detailed X-ray investigation of zeta Puppis IV. Further characterization of the variability

Authors: Yael Naze, Tahina Ramiaramanantsoa, Ian R. Stevens, Ian D. Howarth, Anthony F. J. Moffat

Abstract: Previously, the X-ray emission of zeta Puppis was found to be variable with light curves harbouring "trends" with a typical timescale longer than the exposure length. The origin of these changes was proposed to be linked to large-scale structures in the wind, but further characterization of the variability at high energies was needed. Since then, a number of new X-ray observations have become avai… ▽ More Previously, the X-ray emission of zeta Puppis was found to be variable with light curves harbouring "trends" with a typical timescale longer than the exposure length. The origin of these changes was proposed to be linked to large-scale structures in the wind, but further characterization of the variability at high energies was needed. Since then, a number of new X-ray observations have become available. Furthermore, a cyclic behaviour with a 1.78d period was identified in long optical photometric runs, which is thought to be associated with the launching mechanism of large-scale wind structures. We analysed these new X-ray data, revisited the old data, and compared X-ray with optical data, including when simultaneous. We found that the behaviour in X-rays cannot be explained in terms of a perfect clock because the amplitude and shape of its variations change with time. For example, zeta Puppis was much more strongly variable between 2007 and 2011 than before and after this interval. Comparing the X-ray spectra of the star at maximum and minimum brightness yields no compelling difference beyond the overall flux change: the temperatures, absorptions, and line shapes seem to remain constant, well within errors. The only common feature between X-ray datasets is that the variation amplitudes appear maximum in the medium (0.6-1.2keV) energy band. Finally, no clear and coherent correlation can be found between simultaneous X-ray and optical data. Only a subgroup of observations may be combined coherently with the optical period of 1.78d, although the simultaneous optical behaviour is unknown. The currently available data do not reveal any obvious, permanent, and direct correlation between X-ray and optical variations. The origin of the X-ray variability therefore still needs to be ascertained, highlighting the need for long-term monitoring in multiwavelengths, i.e. X-ray, UV, and optical. △ Less

Submitted 8 December, 2017; v1 submitted 25 October, 2017; originally announced October 2017.

Comments: accepted for publication by A&A

Journal ref: A&A 609, A81 (2018)

BRITE-Constellation high-precision time-dependent photometry of the early-O-type supergiant $ζ$ Puppis unveils the photospheric drivers of its small- and large-scale wind structures

Authors: Tahina Ramiaramanantsoa, Anthony F. J. Moffat, Robert Harmon, Richard Ignace, Nicole St-Louis, Dany Vanbeveren, Tomer Shenar, Herbert Pablo, Noel D. Richardson, Ian D. Howarth, Ian R. Stevens, Caroline Piaulet, Lucas St-Jean, Thomas Eversberg, Andrzej Pigulski, Adam Popowicz, Rainer Kuschnig, Elżbieta Zocłońska, Bram Buysschaert, Gerald Handler, Werner W. Weiss, Gregg A. Wade, Slavek M. Rucinski, Konstanze Zwintz, Paul Luckas , et al. (11 additional authors not shown)

Abstract: 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/SMEI, but this time along with a prominent f… ▽ More 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/SMEI, 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 multi-site spectroscopic monitoring of the star unveiled cyclical modulation of its He II $\lambda4686$ wind emission line with the $1.78$-day rotation period, showing signatures of Corotating Interaction Regions (CIRs) 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 $\lambda4686$ 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. △ Less

Submitted 23 October, 2017; originally announced October 2017.

Comments: 44 pages, 28 figures, 6 Tables; To appear in the Monthly Notices of the Royal Astronomical Society (MNRAS)

Rapid rotators revisited: absolute dimensions of KOI-13

Authors: Ian D. Howarth, Giuseppe Morello

Abstract: We analyse Kepler light-curves of the exoplanet KOI-13b transiting its moderately rapidly rotating (gravity-darkened) parent star. A physical model, with minimal ad hoc free parameters, reproduces the time-averaged light-curve at the ca. 10 parts per million level. We demonstrate that this Roche-model solution allows the absolute dimensions of the system to be determined from the star's projected… ▽ More We analyse Kepler light-curves of the exoplanet KOI-13b transiting its moderately rapidly rotating (gravity-darkened) parent star. A physical model, with minimal ad hoc free parameters, reproduces the time-averaged light-curve at the ca. 10 parts per million level. We demonstrate that this Roche-model solution allows the absolute dimensions of the system to be determined from the star's projected equatorial rotation speed, v(e)sin(i), without any additional assumptions; we find a planetary radius 1.33+/-0.05 R(Jup), stellar polar radius 1.55+/-0.06 R(sun), combined mass M(*) + M(P) (\simeq M*) = 1.47 +/- 0.17 M(sun), and distance d \simeq 370+/-25 pc, where the errors are dominated by uncertainties in relative flux contribution of the visual-binary companion KOI-13B. The implied stellar rotation period is within ca. 5% of the non-orbital, 25.43-hr signal found in the Kepler photometry. We show that the model accurately reproduces independent tomographic observations, and yields an offset between orbital and stellar-rotation angular-momentum vectors of 60.25+/-0.05 degrees. △ Less

Submitted 20 May, 2017; originally announced May 2017.

Comments: Accepted in MNRAS

High-precision stellar limb-darkening in exoplanetary transits

Authors: Giuseppe Morello, Angelos Tsiaras, Ian D. Howarth, Derek Homeier

Abstract: Characterization of the atmospheres of transiting exoplanets relies on accurate measurements of the extent of the optically thick area of the planet at multiple wavelengths with a precision $\lesssim$100 parts per million (ppm). Next-generation instruments onboard the James Webb Space Telescope (JWST) are expected to achieve $\sim$10 ppm precision for several tens of targets. A similar precision c… ▽ More Characterization of the atmospheres of transiting exoplanets relies on accurate measurements of the extent of the optically thick area of the planet at multiple wavelengths with a precision $\lesssim$100 parts per million (ppm). Next-generation instruments onboard the James Webb Space Telescope (JWST) are expected to achieve $\sim$10 ppm precision for several tens of targets. A similar precision can be obtained in modelling only if other astrophysical effects, including the stellar limb-darkening, are accounted for properly. In this paper, we explore the limits on precision due to the mathematical formulas currently adopted to approximate the stellar limb-darkening, and to the use of limb-darkening coefficients obtained either from stellar-atmosphere models or empirically. We propose a new limb-darkening law with two coefficients, `power-2', which outperforms other two-coefficient laws adopted in the literature in most cases, and particularly for cool stars. Empirical limb-darkening based on two-coefficient formulas can be significantly biased, even if the light-curve residuals are nearly photon-noise limited. We demonstrate an optimal strategy to fitting for the four-coefficients limb-darkening in the visible, using prior information on the exoplanet orbital parameters to break some of the degeneracies that otherwise would prevent the convergence of the fit. Infrared observations taken with the James Webb Space Telescope (JWST) will provide accurate measurements of the exoplanet orbital parameters with unprecedented precision, which can be used as priors to improve the stellar limb-darkening characterization, and therefore the inferred exoplanet parameters, from observations in the visible, such as those taken with Kepler/K2, JWST, other past and future instruments. △ Less

Submitted 31 August, 2017; v1 submitted 26 April, 2017; originally announced April 2017.

Journal ref: 2017, AJ, 154, 111

The VLT-FLAMES Tarantula Survey XXVII. Physical parameters of B-type main-sequence binary systems in the Tarantula nebula

Authors: R. Garland, P. L. Dufton, C. J. Evans, P. A. Crowther, I. D. Howarth, A. de Koter, S. E. de Mink, N. J. Grin, N. Langer, D. J. Lennon, C. M. McEvoy, H. Sana, F. R. N. Schneider, S. Símon Díaz, W. D. Taylor, A. Thompson, J. S. Vink

Abstract: A spectroscopic analysis has been undertaken for the B-type multiple systems (excluding those with supergiant primaries) in the VLT-FLAMES Tarantula Survey (VFTS). Projected rotational velocities, $v$sin$i$, for the primaries have been estimated using a Fourier Transform technique and confirmed by fitting rotationally broadened profiles. A subset of 33 systems with $v$sin$i$$\leq$ 80 km/s have bee… ▽ More A spectroscopic analysis has been undertaken for the B-type multiple systems (excluding those with supergiant primaries) in the VLT-FLAMES Tarantula Survey (VFTS). Projected rotational velocities, $v$sin$i$, for the primaries have been estimated using a Fourier Transform technique and confirmed by fitting rotationally broadened profiles. A subset of 33 systems with $v$sin$i$$\leq$ 80 km/s have been analysed using a TLUSTY grid of model atmospheres to estimate stellar parameters and surface abundances for the primaries. The effects of a potential flux contribution from an unseen secondary have also been considered. For 20 targets it was possible to reliably estimate their effective temperatures (T$_{eff}$) but for the other 13 objects it was only possible to provide a constraint of 20,000$\leq$T$_{eff}$$\leq$26,000 K -- the other parameters estimated for these targets will be consequently less reliable. The estimated stellar properties are compared with evolutionary models and are generally consistent with their membership of 30 Doradus, while the nature of the secondaries of 3 SB2 system is discussed. A comparison with a sample of single stars with $v$sin$i$$\leq$ 80 km/s obtained from the VFTS and analysed with the same techniques implies that the atmospheric parameters and nitrogen abundances of the two samples are similar. However, the binary sample may have a lack of primaries with significant nitrogen enhancements, which would be consistent with them having low rotational velocities and having effectively evolved as single stars without significant rotational mixing. This result, which may be actually a consequence of the limitations of the pathfinder investigation presented in this paper, should be considered as a motivation for spectroscopic abundance analysis of large samples of binary stars, with high quality observational data. △ Less

Submitted 24 April, 2017; originally announced April 2017.

Comments: 14 pages, 5 figures

First constraints on the magnetic field strength in extra-Galactic stars: FORS2 observations of Of?p stars in the Magellanic Clouds

Authors: S. Bagnulo, Y. Naze, I. D. Howarth, N. Morrell, J. S. Vink, G. A. Wade, N. Walborn, M. Romaniello, R. Barba

Abstract: Massive O-type stars play a dominant role in our Universe, but many of their properties remain poorly constrained. In the last decade magnetic fields have been detected in all Galactic members of the distinctive Of?p class, opening the door to a better knowledge of all O-type stars. With the aim of extending the study of magnetic massive stars to nearby galaxies, to better understand the role of m… ▽ More Massive O-type stars play a dominant role in our Universe, but many of their properties remain poorly constrained. In the last decade magnetic fields have been detected in all Galactic members of the distinctive Of?p class, opening the door to a better knowledge of all O-type stars. With the aim of extending the study of magnetic massive stars to nearby galaxies, to better understand the role of metallicity in the formation of their magnetic fields and magnetospheres, and to broaden our knowledge of the role of magnetic fields in massive star evolution, we have carried out spectropolarimetry of five extra-Galactic Of?p stars, as well as a couple of dozen neighbouring stars. We have been able to measure magnetic fields with typical error bars from 0.2 to 1.0 kG, depending on the apparent magnitude and on weather conditions. No magnetic field has been firmly detected in any of our measurements, but we have been able to estimate upper limits to the field values of our target stars. One of our targets, 2dFS 936, exhibited an unexpected strengthening of emission lines. We confirm the unusual behaviour of BI 57, which exhibits a 787 d period with two photometric peaks and one spectroscopic maximum. The observed strengthening of the emission lines of 2dFS 936, and the lack of detection of a strong magnetic field in a star with such strong emission lines is at odd with expectations. Together with the unusual periodic behaviour of BI 57, it represents a challenge for the current models of Of?p stars. The limited precision that we obtained in our field measurements (in most cases as a consequence of poor weather) has led to field-strength upper limits that are substantially larger than those typically measured in Galactic magnetic O stars. Further higher precision observations and monitoring are clearly required. △ Less

Submitted 2 March, 2017; originally announced March 2017.

Comments: Accepted by A&A

Journal ref: A&A 601, A136 (2017)

The R136 star cluster dissected with Hubble Space Telescope/STIS. I. Far-ultraviolet spectroscopic census and the origin of HeII 1640 in young star clusters

Authors: Paul A. Crowther, S. M. Caballero-Nieves, K. A. Bostroem, J. Maiz Apellaniz, F. R. N. Schneider, N. R. Walborn, C. R. Angus, I. Brott, A. Bonanos, A. de Koter, S. E. de Mink, C. J. Evans, G. Grafener, A. Herrero, I. D. Howarth, N. Langer, D. J. Lennon, J. Puls, H. Sana, J. S. Vink

Abstract: We introduce a HST/STIS stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS/MAMA spectroscopy of R136 using 17 contiguous 52x0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90% of the 57 sources brighter than m_F555W = 16.0 mag within a radius of 0.5… ▽ More We introduce a HST/STIS stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS/MAMA spectroscopy of R136 using 17 contiguous 52x0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90% of the 57 sources brighter than m_F555W = 16.0 mag within a radius of 0.5 parsec of R136a1, plus 8 additional nearby sources including R136b (O4\,If/WN8). We measure wind velocities for 52 early-type stars from CIV 1548-51, including 16 O2-3 stars. For the first time we spectroscopically classify all Weigelt & Baier members of R136a, which comprise three WN5 stars (a1-a3), two O supergiants (a5-a6) and three early O dwarfs (a4, a7, a8). A complete Hertzsprung-Russell diagram for the most massive O stars in R136 is provided, from which we obtain a cluster age of 1.5+0.3_-0.7 Myr. In addition, we discuss the integrated ultraviolet spectrum of R136, and highlight the central role played by the most luminous stars in producing the prominent HeII 1640 emission line. This emission is totally dominated by very massive stars with initial masses above ~100 Msun. The presence of strong HeII 1640 emission in the integrated light of very young star clusters (e.g A1 in NGC 3125) favours an initial mass function extending well beyond a conventional upper limit of 100 Msun. We include montages of ultraviolet spectroscopy for LMC O stars in the Appendix. Future studies in this series will focus on optical STIS/CCD medium resolution observations. △ Less

Submitted 16 March, 2016; originally announced March 2016.

Comments: 20 pages plus four Appendices providing LMC UV O spectral templates, UV spectral atlas in R136, wind velocities of LMC O stars and photometry of additional R136 sources

Journal ref: MNRAS 458, 624-659 (2016)

A reappraisal of parameters for the putative planet PTFO 8-8695b and its potentially precessing parent star

Authors: Ian D. Howarth

Abstract: Published photometry of fading events in the PTFO 8-8695 system is modelled using improved treatments of stellar geometry, surface intensities, and, particularly, gravity darkening, with a view to testing the planetary-transit hypothesis. Variability in the morphology of fading events can be reproduced by adopting convective-envelope gravity darkening, but near-critical stellar rotation is require… ▽ More Published photometry of fading events in the PTFO 8-8695 system is modelled using improved treatments of stellar geometry, surface intensities, and, particularly, gravity darkening, with a view to testing the planetary-transit hypothesis. Variability in the morphology of fading events can be reproduced by adopting convective-envelope gravity darkening, but near-critical stellar rotation is required. This leads to inconsistencies with spectroscopic observations; the model also predicts substantial photometric variability associated with stellar precession, contrary to observations. Furthermore, the empirical ratio of orbital to rotational angular momenta is at odds with physically plausible values. An exoplanet transiting a precessing, gravity-darkened star may not be the correct explanation of periodic fading events in this system. △ Less

Submitted 27 January, 2016; originally announced January 2016.

The MiMeS Survey of Magnetism in Massive Stars: Introduction and overview

Authors: G. A. Wade, C. Neiner, E. Alecian, J. H. Grunhut, V. Petit, B. de Batz, D. A. Bohlender, D. H. Cohen, H. F. Henrichs, O. Kochukhov, J. D. Landstreet, N. Manset, F. Martins, S. Mathis, M. E. Oksala, S. P. Owocki, Th. Rivinius, M. E. Shultz, J. O. Sundqvist, R. H. D. Townsend, A. ud-Doula, J. -C. Bouret, J. Braithwaite, M. Briquet, A. C. Carciofi , et al. (25 additional authors not shown)

Abstract: The MiMeS project is a large-scale, high resolution, sensitive spectropolarimetric investigation of the magnetic properties of O and early B type stars. Initiated in 2008 and completed in 2013, the project was supported by 3 Large Program allocations, as well as various programs initiated by independent PIs and archival resources. Ultimately, over 4800 circularly polarized spectra of 560 O and B s… ▽ More The MiMeS project is a large-scale, high resolution, sensitive spectropolarimetric investigation of the magnetic properties of O and early B type stars. Initiated in 2008 and completed in 2013, the project was supported by 3 Large Program allocations, as well as various programs initiated by independent PIs and archival resources. Ultimately, over 4800 circularly polarized spectra of 560 O and B stars were collected with the instruments ESPaDOnS at the Canada-France-Hawaii Telescope, Narval at the Télescope Bernard Lyot, and HARPSpol at the European Southern Observatory La Silla 3.6m telescope, making MiMeS by far the largest systematic investigation of massive star magnetism ever undertaken. In this paper, the first in a series reporting the general results of the survey, we introduce the scientific motivation and goals, describe the sample of targets, review the instrumentation and observational techniques used, explain the exposure time calculation designed to provide sensitivity to surface dipole fields above approximately 100 G, discuss the polarimetric performance, stability and uncertainty of the instrumentation, and summarize the previous and forthcoming publications. △ Less

Submitted 26 November, 2015; originally announced November 2015.

Comments: Accepted by MNRAS

Journal ref: MNRAS, 456, 2-22 (2016)

Limb-darkening coefficients from line-blanketed non-LTE hot-star model atmospheres

Authors: D. C. Reeve, I. D. Howarth

Abstract: We present grids of limb-darkening coefficients computed from non-LTE, line-blanketed TLUSTY model atmospheres, covering effective-temperature and surface-gravity ranges of 15--55kK and 4.75 dex (cgs) down to the effective Eddington limit, at 1x, 1x, 0.5x (LMC), 0.2x (SMC), and 0.1x solar. Results are given for the Bessell UBVRIJKHL, Sloan ugriz, Stromgren ubvy, WFCAM ZYJHK, Hipparcos, Kepler, and… ▽ More We present grids of limb-darkening coefficients computed from non-LTE, line-blanketed TLUSTY model atmospheres, covering effective-temperature and surface-gravity ranges of 15--55kK and 4.75 dex (cgs) down to the effective Eddington limit, at 1x, 1x, 0.5x (LMC), 0.2x (SMC), and 0.1x solar. Results are given for the Bessell UBVRIJKHL, Sloan ugriz, Stromgren ubvy, WFCAM ZYJHK, Hipparcos, Kepler, and Tycho passbands, in each case characterized by several different limb-darkening `laws'. We examine the sensitivity of limb darkening to temperature, gravity, metallicity, microturbulent velocity, and wavelength, and make a comparison with LTE models. The dependence on metallicity is very weak, but limb darkening is a moderately strong function of log(g) in this temperature regime. △ Less

Submitted 6 November, 2015; originally announced November 2015.

Comments: Accepted for MNRAS

The VLT-FLAMES Tarantula Survey. XXIII: two massive double-lined binaries in 30~Doradus

Authors: Ian D. Howarth, P. L. Dufton, P. R. Dunstall, C. J. Evans, L. A. Almeida, A. Z. Bonanos, J. S. Clark, N. Langer, H. Sana, S. Simon-Diaz, I. Soszynski, W. D. Taylor

Abstract: We investigate the characteristics of two newly discovered short-period, double-lined, massive binary systems, VFTS 450 (O9.7$\;$II--Ib$\,$+$\,$O7::) and VFTS 652 (B1$\;$Ib$\,+\,$O9:$\;$III:). We perform model-atmosphere analyses to characterise the photospheric properties of both members of each binary (denoting the `primary' as the spectroscopically more conspicuous component). Radial velocities… ▽ More We investigate the characteristics of two newly discovered short-period, double-lined, massive binary systems, VFTS 450 (O9.7$\;$II--Ib$\,$+$\,$O7::) and VFTS 652 (B1$\;$Ib$\,+\,$O9:$\;$III:). We perform model-atmosphere analyses to characterise the photospheric properties of both members of each binary (denoting the `primary' as the spectroscopically more conspicuous component). Radial velocities and optical photometry are used to estimate the binary-system parameters. We estimate $T_{\rm eff}=27$ kK, $\log{(g)}=2.9$ (cgs) for the VFTS 450 primary spectrum (34kK, 3.6: for the secondary spectrum); and $T_{\rm eff} = 22$kK, $\log{(g)}=2.8$ for the VFTS 652 primary spectrum (35kK, 3.7: for the secondary spectrum). Both primaries show surface nitrogen enrichments (of more than 1 dex for VFTS 652), and probable moderate oxygen depletions relative to reference LMC abundances. We determine orbital periods of 6.89d and 8.59d for VFTS 450 and VFTS 652, respectively, and argue that the primaries must be close to filling their Roche lobes. Supposing this to be the case, we estimate component masses in the range $\sim$20--50M$_\odot$. The secondary spectra are associated with the more massive components, suggesting that both systems are high-mass analogues of classical Algol systems, undergoing case-A mass transfer. Difficulties in reconciling the spectroscopic analyses with the light-curves and with evolutionary considerations suggest that the secondary spectra are contaminated by (or arise in) accretion disks. △ Less

Submitted 24 August, 2015; originally announced August 2015.

Journal ref: A&A 582, A73 (2015)

Spectral Variations of Of?p Oblique Magnetic Rotator Candidates in the Magellanic Clouds

Authors: Nolan R. Walborn, Nidia I. Morrell, Yael Naze, Gregg A. Wade, Stefano Bagnulo, Rodolfo H. Barba, Jesus Maiz Apellaniz, Ian D. Howarth, Christopher J. Evans, Alfredo Sota

Abstract: Optical spectroscopic monitoring has been conducted of two O stars in the Small and one in the Large Magellanic Cloud, the spectral characteristics of which place them in the Of?p category, which has been established in the Galaxy to consist of oblique magnetic rotators. All of these Magellanic stars show systematic spectral variations typical of the Of?p class, further strengthening their magneti… ▽ More Optical spectroscopic monitoring has been conducted of two O stars in the Small and one in the Large Magellanic Cloud, the spectral characteristics of which place them in the Of?p category, which has been established in the Galaxy to consist of oblique magnetic rotators. All of these Magellanic stars show systematic spectral variations typical of the Of?p class, further strengthening their magnetic candidacy to the point of virtual certainty. The spectral variations are related to photometric variations derived from OGLE data by Naze et al. (2015) in a parallel study, which yields rotational periods for two of them. Now circular spectropolarimetry is required to measure their fields, and ultraviolet spectroscopy to further characterize their low-metallicity, magnetically confined winds, in support of hydrodynamical analyses. △ Less

Submitted 9 July, 2015; originally announced July 2015.

Comments: 18 pages, 6 figures, accepted for publication by AJ

Journal ref: AJ, 150, id 99 (7p, 2015)

The VLT-FLAMES Tarantula Survey XXII. Multiplicity properties of the B-type stars

Authors: P. R. Dunstall, P. L. Dufton, H. Sana, C. J. Evans, I. D. Howarth, S. Simón-Díaz, S. E. de Mink, N. Langer, J. Maíz Apellániz, W. D. Taylor

Abstract: We investigate the multiplicity properties of 408 B-type stars observed in the 30 Doradus region of the Large Magellanic Cloud with multi-epoch spectroscopy from the VLT-FLAMES Tarantula Survey (VFTS). We use a cross-correlation method to estimate relative radial velocities from the helium and metal absorption lines for each of our targets. Objects with significant radial-velocity variations (and… ▽ More We investigate the multiplicity properties of 408 B-type stars observed in the 30 Doradus region of the Large Magellanic Cloud with multi-epoch spectroscopy from the VLT-FLAMES Tarantula Survey (VFTS). We use a cross-correlation method to estimate relative radial velocities from the helium and metal absorption lines for each of our targets. Objects with significant radial-velocity variations (and with an amplitude larger than 16 km/s) are classified as spectroscopic binaries. We find an observed spectroscopic binary fraction (defined by periods of <10^3.5 d and mass ratios >0.1) for the B-type stars, f_B(obs) = 0.25 +/- 0.02, which appears constant across the field of view, except for the two older clusters (Hodge 301 and SL 639). These two clusters have significantly lower fractions of 0.08 +/- 0.08 and 0.10 +/- 0.09, respectively. Using synthetic populations and a model of our observed epochs and their potential biases, we constrain the intrinsic multiplicity properties of the dwarf and giant (i.e. relatively unevolved) B-type stars in 30 Dor. We obtain a present-day binary fraction f_B(true) = 0.58 +/- 0.11, with a flat period distribution. Within the uncertainties, the multiplicity properties of the B-type stars agree with those for the O stars in 30 Dor from the VFTS. △ Less

Submitted 26 May, 2015; originally announced May 2015.

Comments: Accepted by A&A

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

Revisiting Spitzer transit observations with Independent Component Analysis: new results for the GJ436 system

Authors: G. Morello, I. P. Waldmann, G. Tinetti, I. D. Howarth, G. Micela, F. Allard

Abstract: We analyzed four Spitzer/IRAC observations at 3.6 and 4.5 μm of the primary transit of the exoplanet GJ436b, by using blind source separation techniques. These observations are important to investigate the atmospheric composition of the planet GJ436b. Previous analyses claimed strong inter-epoch variations of the transit parameters due to stellar variability, casting doubts on the possibility to e… ▽ More We analyzed four Spitzer/IRAC observations at 3.6 and 4.5 μm of the primary transit of the exoplanet GJ436b, by using blind source separation techniques. These observations are important to investigate the atmospheric composition of the planet GJ436b. Previous analyses claimed strong inter-epoch variations of the transit parameters due to stellar variability, casting doubts on the possibility to extract conclusively an atmospheric signal; those analyses also reported discrepant results, hence the necessity of this reanalysis. The method we used has been proposed in Morello et al. (2014) to analyze 3.6 μm transit light-curves of the hot Jupiter HD189733b; it performes an Independent Component Analysis (ICA) on a set of pixel-light-curves, i.e. time series read by individual pixels, from the same photometric observation. Our method only assumes the independence of instrumental and astrophysical signals, and therefore guarantees a higher degree of objectivity compared to parametric detrending techniques published in the literature. The datasets we analyzed in this paper represent a more challenging test compared to the previous ones. Contrary to previous results reported in the literature, our results (1) do not support any detectable inter-epoch variations of orbital and stellar parameters, (2) are photometrically stable at the level 10e-4 in the IR, and (3) the transit depth measurements at the two wavelengths are consistent within 1σ. We also (4) detect a possible transit duration variation (TDV) of 80 s (2 σ significance level), that has not been pointed out in the literature, and (5) confirm no transit timing variations (TTVs) >30 s. △ Less

Submitted 23 January, 2015; originally announced January 2015.

Comments: accepted for publication in The Astrophysical Journal

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

Time-series photometry of the O4 I(n)fp star zeta Puppis

Authors: Ian D. Howarth, Ian R. Stevens

Abstract: We report a time-series analysis of the O4 I(n)fp star zeta Pup, based on optical photometry obtained with the SMEI instrument on the Coriolis satellite, 2003--2006. A single astrophysical signal is found, with P = (1.780938 \pm 0.000093) d and a mean semi-amplitude of (6.9 \pm 0.3) mmag. There is no evidence for persistent coherent signals with semi-amplitudes in excess of ca. 2~mmag on any of th… ▽ More We report a time-series analysis of the O4 I(n)fp star zeta Pup, based on optical photometry obtained with the SMEI instrument on the Coriolis satellite, 2003--2006. A single astrophysical signal is found, with P = (1.780938 \pm 0.000093) d and a mean semi-amplitude of (6.9 \pm 0.3) mmag. There is no evidence for persistent coherent signals with semi-amplitudes in excess of ca. 2~mmag on any of the timescales previously reported in the literature. In particular, there is no evidence for a signature of the proposed rotation period, ca. 5.1~days; zeta Pup is therefore probably not an oblique magnetic rotator. The 1.8-day signal varies in amplitude by a factor ca. 2 on timescales of 10--100d (and probably by more on longer timescales), and exhibits modest excursions in phase, but there is no evidence for systematic changes in period over the 1000-d span of our observations. Rotational modulation and stellar-wind variability appear to be unlikely candidates for the underlying mechanism; we suggest that the physical origin of the signal may be pulsation associated with low-l oscillatory convection modes. △ Less

Submitted 18 September, 2014; originally announced September 2014.

Comments: MNRAS, in press

A new look at Spitzer primary transit observations of the exoplanet HD189733b

Authors: Giuseppe Morello, Ingo P. Waldmann, Giovanna Tinetti, Giovanni Peres, Giuseppina Micela, Ian D. Howarth

Abstract: Blind source separation techniques are used to reanalyse two exoplanetary transit lightcurves of the exoplanet HD189733b recorded with the IR camera IRAC on board the Spitzer Space Telescope at 3.6$μ$m during the "cold" era. These observations, together with observations at other IR wavelengths, are crucial to characterise the atmosphere of the planet HD189733b. Previous analyses of the same datas… ▽ More Blind source separation techniques are used to reanalyse two exoplanetary transit lightcurves of the exoplanet HD189733b recorded with the IR camera IRAC on board the Spitzer Space Telescope at 3.6$μ$m during the "cold" era. These observations, together with observations at other IR wavelengths, are crucial to characterise the atmosphere of the planet HD189733b. Previous analyses of the same datasets reported discrepant results, hence the necessity of the reanalyses. The method we used here is based on the Independent Component Analysis (ICA) statistical technique, which ensures a high degree of objectivity. The use of ICA to detrend single photometric observations in a self-consistent way is novel in the literature. The advantage of our reanalyses over previous work is that we do not have to make any assumptions on the structure of the unknown instrumental systematics. Such "admission of ignorance" may result in larger error bars than reported in the literature, up to a factor $1.6$. This is a worthwhile trade-off for much higher objectivity, necessary for trustworthy claims. Our main results are (1) improved and robust values of orbital and stellar parameters, (2) new measurements of the transit depths at 3.6$μ$m, (3) consistency between the parameters estimated from the two observations, (4) repeatability of the measurement within the photometric level of $\sim 2 \times 10^{-4}$ in the IR, (5) no evidence of stellar variability at the same photometric level within 1 year. △ Less

Submitted 12 March, 2014; originally announced March 2014.

Comments: 43 pages, 18 figures

Journal ref: Morello et al. 2014, ApJ, 786, 22

Amplitude variability in satellite photometry of the non-radially pulsating O9.5V star zeta Oph

Authors: Ian D. Howarth, K. J. F. Goss, I. R. Stevens, W. J. Chaplin, Y. Elsworth

Abstract: We report a time-series analysis of satellite photometry of the non-radially pulsating Oe star zeta Oph, principally using data from SMEI obtained 2003--2008, but augmented with MOST and WIRE results. Amplitudes of the strongest photometric signals, at 5.18, 2.96, and 2.67/d, each vary independently over the 6-year monitoring period (from ca. 30 to <2 mmag at 5.18/d), on timescales of hundreds of… ▽ More We report a time-series analysis of satellite photometry of the non-radially pulsating Oe star zeta Oph, principally using data from SMEI obtained 2003--2008, but augmented with MOST and WIRE results. Amplitudes of the strongest photometric signals, at 5.18, 2.96, and 2.67/d, each vary independently over the 6-year monitoring period (from ca. 30 to <2 mmag at 5.18/d), on timescales of hundreds of days. Signals at 7.19/d and 5.18/d have persisted (or recurred) for around two decades. Supplementary spectroscopic observations show an H-alpha emission episode in 2006; this coincided with small increases in amplitudes of the three strongest photometric signals. △ Less

Submitted 8 March, 2014; v1 submitted 26 February, 2014; originally announced February 2014.

Comments: MNRAS, in press

Emission-line stars in the LMC: the Armagh survey, and a metacatalogue

Authors: Ian D. Howarth

Abstract: [Aims] Accurate astrometry is required to reliably cross-match 20th-century photographic catalogues against 21st-century digital surveys. The present work provides modern-era identifications and astrometry for the 801 emission-line objects "of stellar appearance" in the Armagh survey (the largest of its nature to date). [Methods] Targets have been individually identified in digital images using th… ▽ More [Aims] Accurate astrometry is required to reliably cross-match 20th-century photographic catalogues against 21st-century digital surveys. The present work provides modern-era identifications and astrometry for the 801 emission-line objects "of stellar appearance" in the Armagh survey (the largest of its nature to date). [Methods] Targets have been individually identified in digital images using the Armagh Atlas and, in most cases, unambiguously matched to entries in the UCAC astrometric catalogues. [Results] Astrometry with sub-arcsecond precision is now available for all the major photographic spectroscopic surveys of the LMC. The results are used to compile an annotated metacatalogue of 1675 individual, spectroscopically identified candidate H-alpha-emission stars, including detailed cross-matching between catalogues, and resolving many (though not all) identification ambiguities in individual primary sources. △ Less

Submitted 7 June, 2013; originally announced June 2013.

Phase-resolved ultraviolet spectroscopy of the magnetic Of?p star HD 191612

Authors: W. L. F. Marcolino, J. -C. Bouret, J. O. Sundqvist, N. R. Walborn, A. W. Fullerton, I. D. Howarth, G. A. Wade, A. ud-Doula

Abstract: We present for the first time phase-resolved UV spectroscopy of an Of?p star, namely, HD 191612. The observations were acquired with the Space Telescope Imaging Spectrograph (STIS) on-board the Hubble Space Telescope (HST). We report the variability observed in the main photospheric and wind features and compare the results with previous findings for the Of?p star HD 108. We show that UV line stre… ▽ More We present for the first time phase-resolved UV spectroscopy of an Of?p star, namely, HD 191612. The observations were acquired with the Space Telescope Imaging Spectrograph (STIS) on-board the Hubble Space Telescope (HST). We report the variability observed in the main photospheric and wind features and compare the results with previous findings for the Of?p star HD 108. We show that UV line strengths, H(alpha), and longitudinal magnetic field, vary coherently according to the rotational period (P = 537.6d), providing additional support for the magnetic oblique rotator scenario. The stellar and wind parameters of HD 191612 are obtained based on NLTE expanding atmosphere models. The peculiar wind line profile variations revealed by the new STIS data - not reproduced by 1D atmosphere models - are addressed through non-spherical MHD simulations coupled with radiative transfer. The basic aspects of the UV variability observed are explained and the structure of the dynamical magnetosphere of HD 191612 is discussed. △ Less

Submitted 19 February, 2013; originally announced February 2013.

Comments: Accepted for publication in MNRAS (9 pages, 6 figures, 3 tables)

The VLT-FLAMES Tarantula Survey X: Evidence for a bimodal distribution of rotational velocities for the single early B-type stars

Authors: P. L. Dufton, N. Langer, P. R. Dunstall, C. J. Evans, I. Brott, S. E. de Mink, I. D. Howarth, M. Kennedy, C. McEvoy, A. T. Potter, O. H. Ramírez-Agudelo, H. Sana, S. Simón-Díaz, W. Taylor, J. S. Vink

Abstract: Aims: Projected rotational velocities (\vsini) have been estimated for 334 targets in the VLT-FLAMES Tarantula survey that do not manifest significant radial velocity variations and are not supergiants. They have spectral types from approximately O9.5 to B3. The estimates have been analysed to infer the underlying rotational velocity distribution, which is critical for understanding the evolution… ▽ More Aims: Projected rotational velocities (\vsini) have been estimated for 334 targets in the VLT-FLAMES Tarantula survey that do not manifest significant radial velocity variations and are not supergiants. They have spectral types from approximately O9.5 to B3. The estimates have been analysed to infer the underlying rotational velocity distribution, which is critical for understanding the evolution of massive stars. Methods: Projected rotational velocities were deduced from the Fourier transforms of spectral lines, with upper limits also being obtained from profile fitting. For the narrower lined stars, metal and non-diffuse helium lines were adopted, and for the broader lined stars, both non-diffuse and diffuse helium lines; the estimates obtained using the different sets of lines are in good agreement. The uncertainty in the mean estimates is typically 4% for most targets. The iterative deconvolution procedure of Lucy has been used to deduce the probability density distribution of the rotational velocities. Results: Projected rotational velocities range up to approximately 450 \kms and show a bi-modal structure. This is also present in the inferred rotational velocity distribution with 25% of the sample having $0\leq$\ve$\leq$100\,\kms and the high velocity component having \ve$\sim 250$\,\kms. There is no evidence from the spatial and radial velocity distributions of the two components that they represent either field and cluster populations or different episodes of star formation. Be-type stars have also been identified. Conclusions: The bi-modal rotational velocity distribution in our sample resembles that found for late-B and early-A type stars. While magnetic braking appears to be a possible mechanism for producing the low-velocity component, we can not rule out alternative explanations. △ Less

Submitted 11 December, 2012; originally announced December 2012.

Comments: to be publisged in A&A

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

Improved astrometry for the Bohannan & Epps catalogue

Authors: Ian D. Howarth

Abstract: Aims: Accurate astrometry is required to reliably cross-match 20th-century catalogues against 21st-century surveys. The present work aims to provide such astrometry for the 625 entries of the Bohannan & Epps (BE74) catalogue of H$α$ emission-line stars. Methods: BE74 targets have been individually identified in digital images and, in most cases, unambiguously matched to entries in the UCAC4 astrom… ▽ More Aims: Accurate astrometry is required to reliably cross-match 20th-century catalogues against 21st-century surveys. The present work aims to provide such astrometry for the 625 entries of the Bohannan & Epps (BE74) catalogue of H$α$ emission-line stars. Methods: BE74 targets have been individually identified in digital images and, in most cases, unambiguously matched to entries in the UCAC4 astrometric catalogue. Results: Sub-arcsecond astrometry is now available for almost all BE74 stars. Several identification errors in the literature illustrate the perils of relying solely on positional coincidences using poorer-quality astrometry. △ Less

Submitted 5 October, 2012; originally announced October 2012.

Comments: 3 pages, 1 figure. Accepted in A&A

Two new Wolf-Rayet stars in the LMC

Authors: Ian D. Howarth, Nolan R. Walborn

Abstract: We report the discovery of two previously unknown WN3 stars in the Large Magellanic Cloud. Both are bright (15th magnitude), isolated, and located in regions covered in earlier surveys, although both are relatively weak-lined. We suggest that there may be $\mathcal{O}(10)$ remaining undiscovered WNE stars in the LMC. We report the discovery of two previously unknown WN3 stars in the Large Magellanic Cloud. Both are bright (15th magnitude), isolated, and located in regions covered in earlier surveys, although both are relatively weak-lined. We suggest that there may be $\mathcal{O}(10)$ remaining undiscovered WNE stars in the LMC. △ Less

Submitted 30 July, 2012; originally announced July 2012.

A dynamical magnetosphere model for periodic Halpha emission from the slowly rotating magnetic O star HD191612

Authors: Jon O. Sundqvist, Asif ud-Doula, Stanley P. Owocki, Richard H. D. Townsend, Ian D. Howarth, Gregg A. Wade, the MiMeS Collaboration

Abstract: The magnetic O-star HD191612 exhibits strongly variable, cyclic Balmer line emission on a 538-day period. We show here that its variable Halpha emission can be well reproduced by the rotational phase variation of synthetic spectra computed directly from full radiation magneto-hydrodynamical simulations of a magnetically confined wind. In slow rotators such as HD191612, wind material on closed magn… ▽ More The magnetic O-star HD191612 exhibits strongly variable, cyclic Balmer line emission on a 538-day period. We show here that its variable Halpha emission can be well reproduced by the rotational phase variation of synthetic spectra computed directly from full radiation magneto-hydrodynamical simulations of a magnetically confined wind. In slow rotators such as HD191612, wind material on closed magnetic field loops falls back to the star, but the transient suspension of material within the loops leads to a statistically overdense, low velocity region around the magnetic equator, causing the spectral variations. We contrast such "dynamical magnetospheres" (DMs) with the more steady-state "centrifugal magnetospheres" of stars with rapid rotation, and discuss the prospects of using this DM paradigm to explain periodic line emission from also other non-rapidly rotating magnetic massive stars. △ Less

Submitted 5 March, 2012; originally announced March 2012.

Comments: 5 pages, 5 figures, accepted for publication in MNRAS letters

HST/STIS spectroscopy of the magnetic Of?p star HD 108: the low state at ultraviolet wavelengths

Authors: W. L. F. Marcolino, J. -C. Bouret, N. R. Walborn, I. D. Howarth, Y. Naze, A. W. Fullerton, G. A. Wade, D. J. Hillier, A. Herrero

Abstract: We present the first ultraviolet spectrum of the peculiar, magnetic Of?p star HD 108 obtained in its spectroscopic low state. The new data, obtained with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope, reveal significant changes compared to IUE spectra obtained in the high state: N V 1240, Si IV 1400, and C IV 1550 present weaker P-Cygni profiles (less absorption) in… ▽ More We present the first ultraviolet spectrum of the peculiar, magnetic Of?p star HD 108 obtained in its spectroscopic low state. The new data, obtained with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope, reveal significant changes compared to IUE spectra obtained in the high state: N V 1240, Si IV 1400, and C IV 1550 present weaker P-Cygni profiles (less absorption) in the new data, while N IV 1718 absorption is deeper, without the clear wind signature evident in the high state. Such changes contrast with those found in other magnetic massive stars, where more absorption is observed in the resonance doublets when the sightline is close to the plane of the magnetic equator. The new data show also that the photospheric Fe IV forest, at 1600--1700 angstroms, has strengthened compared to previous observations. The ultraviolet variability is large compared to that found in typical, non-magnetic O stars, but moderate when compared to the high-/low-state changes reported in the optical spectrum of HD 108 over several decades. We use non-LTE expanding-atmosphere models to analyze the new STIS observations. Overall, the results are in accord with a scenario in which the optical variability is mainly produced by magnetically constrained gas, close to the photosphere. The relatively modest changes found in the main ultraviolet wind lines suggest that the stellar wind is not substantially variable on a global scale. Nonetheless, multidimensional radiative-transfer models may be needed to understand some of the phenomena observed. △ Less

Submitted 27 February, 2012; originally announced February 2012.

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

Journal ref: MNRAS, 422, 2314-2321 (2012)

On the nature of the galactic early-B hypergiants

Authors: J. S. Clark, F. Najarro, I. Negueruela, B. W. Ritchie, M. A. Urbaneja, I. D. Howarth

Abstract: Despite their importance to a number of astrophysical fields, the lifecycles of very massive stars are still poorly defined. In order to address this shortcoming, we present a detailed quantitative study of the physical properties of four early-B hypergiants (BHGs); Cyg OB2 #12, zeta Sco, HD190603 and BP Cru. These are combined with an analysis of their long-term spectroscopic and photometric beha… ▽ More Despite their importance to a number of astrophysical fields, the lifecycles of very massive stars are still poorly defined. In order to address this shortcoming, we present a detailed quantitative study of the physical properties of four early-B hypergiants (BHGs); Cyg OB2 #12, zeta Sco, HD190603 and BP Cru. These are combined with an analysis of their long-term spectroscopic and photometric behaviour in order to determine their evolutionary status. The long-term datasets revealed that they are remarkably stable over long periods (>40yr), with the possible exception of zeta Sco prior to the 20th century, in contrast to the typical excursions that characterise luminous blue variables (LBVs). Zeta Sco, HD190603 and BP Cru possess physical properties intermediate between B supergiants and LBVs; we therefore suggest that BHGs are the immediate descendants and progenitors (respectively) of such stars (for initial masses in the range ~30-60Msun). In contrast, while the wind properties of Cyg OB2 #12 are consistent with this hypothesis, the combination of extreme luminosity and spectroscopic mass (~110Msun) and comparatively low temperature means it cannot be accommodated in such a scheme. Likewise, despite its co-location with several LBVs above the Humphreys-Davidson (HD) limit, the lack of long term variability and its unevolved chemistry apparently excludes such an identification. Since such massive stars are not expected to evolve to such cool temperatures, the properties of Cyg OB2 #12 are difficult to understand under current evolutionary paradigms. [ABRIDGED] △ Less

Submitted 17 February, 2012; originally announced February 2012.

Comments: 36 pages, 19 figures (of which 17 pages are online supplemental material). Accepted for publication in Astronomy and Astrophysics

The VLT-FLAMES Tarantula Survey: The fastest rotating O-type star and shortest period LMC pulsar - remnants of a supernova disrupted binary?

Authors: P. L. Dufton, P. R. Dunstall, C. J. Evans, I. Brott, M. Cantiello, A. de Koter, S. E. de Mink, M. Fraser, V. Hénault-Brunet, I. D. Howarth, N. Langer, D. J. Lennon, N. Markova, H. Sana, W. D. Taylor

Abstract: We present a spectroscopic analysis of an extremely rapidly rotating late O-type star, VFTS102, observed during a spectroscopic survey of 30 Doradus. VFTS102 has a projected rotational velocity larger than 500\kms\ and probably as large as 600\kms; as such it would appear to be the most rapidly rotating massive star currently identified. Its radial velocity differs by 40\kms\ from the mean for 30… ▽ More We present a spectroscopic analysis of an extremely rapidly rotating late O-type star, VFTS102, observed during a spectroscopic survey of 30 Doradus. VFTS102 has a projected rotational velocity larger than 500\kms\ and probably as large as 600\kms; as such it would appear to be the most rapidly rotating massive star currently identified. Its radial velocity differs by 40\kms\ from the mean for 30 Doradus, suggesting that it is a runaway. VFTS102 lies 12 pcs from the X-ray pulsar PSR J0537-6910 in the tail of its X-ray diffuse emission. We suggest that these objects originated from a binary system with the rotational and radial velocities of VFTS102 resulting from mass transfer from the progenitor of PSR J0537-691 and the supernova explosion respectively. △ Less

Submitted 6 December, 2011; v1 submitted 1 November, 2011; originally announced November 2011.

Comments: Accepted by ApJL; v2: minor edits in proofs

Journal ref: 2011, ApJ, 743, L22

The spectral variability and magnetic field characteristics of the Of?p star HD 148937

Authors: G. A. Wade, J. Grunhut, G. Graefener, I. D. Howarth, F. Martins, V. Petit, J. S. Vink, S. Bagnulo, C. P. Folsom, Y. Nazé, N. R. Walborn, R. H. D. Townsend, C. J. Evans, the MiMeS Collaboration

Abstract: We report magnetic and spectroscopic observations and modeling of the Of?p star HD 148937 within the context of the MiMeS LP at the CFHT. Thirty-two high signal-to-noise ratio circularly polarised (Stokes V) spectra and 13 unpolarised (Stokes I) spectra of HD 148937 were acquired in 2009 and 2010. A definite detection of a Stokes V Zeeman signature is obtained in the grand mean of all observations… ▽ More We report magnetic and spectroscopic observations and modeling of the Of?p star HD 148937 within the context of the MiMeS LP at the CFHT. Thirty-two high signal-to-noise ratio circularly polarised (Stokes V) spectra and 13 unpolarised (Stokes I) spectra of HD 148937 were acquired in 2009 and 2010. A definite detection of a Stokes V Zeeman signature is obtained in the grand mean of all observations (in both LSD mean profiles and individual spectral lines). The longitudinal magnetic field inferred from the Stokes V LSD profiles is consistently negative, in contrast to the essentially zero field strength measured from the diagnostic null profiles. A period search of equivalent width measurements confirms the previously-reported 7.03 d variability period. The variation of equivalent widths is not strictly periodic: we present evidence for evolution of the amount or distribution of circumstellar plasma. Interpreting the 7.03 d period as the stellar rotational period within the context of the ORM, we have phased the equivalent widths and longitudinal field measurements. The longitudinal field measurements show a weak sinusoidal variation of constant sign, with extrema out of phase with the Hα variation by about 0.25 cycles. The inferred magnetic configuration confirms the suggestion of Nazé et al (2010), who proposed that the weaker variability of HD 148937 as compared to other members of this class is a consequence of the stellar geometry. Based on the derived magnetic properties and published wind characteristics, we find a wind magnetic confinement parameter η\ast \simeq 20 and rotation parameter W = 0.12, supporting a picture in which the Halpha emission and other line variability have their origin in an oblique, rigidly rotating magnetospheric structure resulting from a magnetically channeled wind. (Abridged.) △ Less

Submitted 24 September, 2011; v1 submitted 24 August, 2011; originally announced August 2011.

Comments: 13 pages, MNRAS. Version 2, small change to Fig. 10

Journal ref: MNRAS, 419, 2459-2471 (2012)

On stellar limb darkening and exoplanetary transits

Authors: Ian D. Howarth

Abstract: This paper examines how to compare stellar limb-darkening coefficients evaluated from model atmospheres with those derived from photometry. Different characterizations of a given model atmosphere can give quite different numerical results (even for a given limb-darkening `law'), while light-curve analyses yield limb-darkening coefficients that are dependent on system geometry, and that are not dir… ▽ More This paper examines how to compare stellar limb-darkening coefficients evaluated from model atmospheres with those derived from photometry. Different characterizations of a given model atmosphere can give quite different numerical results (even for a given limb-darkening `law'), while light-curve analyses yield limb-darkening coefficients that are dependent on system geometry, and that are not directly comparable to any model-atmosphere representation. These issues are examined in the context of exoplanetary transits, which offer significant advantages over traditional binary-star eclipsing systems in the study of stellar limb darkening. `Like for like' comparisons between light-curve analyses and new model-atmosphere results, mediated by synthetic photometry, are conducted for a small sample of stars. Agreement between the resulting synthetic-photometry/atmosphere-model (SPAM) limb-darkening coefficients and empirical values ranges from very good to quite poor, even though the targets investigated show only a small dispersion in fundamental stellar parameters. △ Less

Submitted 5 August, 2011; v1 submitted 23 June, 2011; originally announced June 2011.

Comments: Accepted, MNRAS, 4/8/2011

Confirmation of the magnetic oblique rotator model for the Of?p star HD 191612

Authors: G. A. Wade, I. D. Howarth, R. H. D. Townsend, J. H. Grunhut, M. Shultz, J. -C. Bouret, A. Fullerton, W. Marcolino, F. Martins, Y. Nazé, A. ud Doula, N. R. Walborn, J. -F. Donati, the MiMeS Collaboration

Abstract: This paper reports high-precision Stokes V spectra of HD 191612 acquired using the ESPaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope, in the context of the Magnetism in Massive stars (MiMeS) Project. Using measurements of the equivalent width of the Halpha line and radial velocities of various metallic lines, we have updated both the spectroscopic and orbital ephemerides of this st… ▽ More This paper reports high-precision Stokes V spectra of HD 191612 acquired using the ESPaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope, in the context of the Magnetism in Massive stars (MiMeS) Project. Using measurements of the equivalent width of the Halpha line and radial velocities of various metallic lines, we have updated both the spectroscopic and orbital ephemerides of this star. We confirm the presence of a strong magnetic field in the photosphere of HD 191612, and detect its variability. We establish that the longitudinal field varies in a manner consistent with the spectroscopic period of 537.6 d, in an approximately sinusoidal fashion. This demonstrates a firm connection between the magnetic field and the processes responsible for the line and continuum variability. Interpreting the variation of the longitudinal magnetic field within the context of the dipole oblique rotator model we obtain a best-fit surface magnetic field model with obliquity beta=67\pm 5 deg and polar strength Bd=2450\pm 400 G . The inferred magnetic field strength implies an equatorial wind magnetic confinement parameter eta*~50, supporting a picture in which the Halpha emission and photometric variability have their origin in an oblique, rigidly rotating magnetospheric structure resulting from a magnetically channeled wind. This interpretation is supported by our successful Monte Carlo radiative transfer modeling of the photometric variation, which assumes the enhanced plasma densities in the magnetic equatorial plane above the star implied by such a picture. Predictions of the continuum linear polarisation resulting from Thompson scattering from the magnetospheric material indicate that the Stokes Q and U variations are highly sensitive to the magnetospheric geometry, and that expected amplitudes are in the range of current instrumentation. (abridged) △ Less

Submitted 15 June, 2011; originally announced June 2011.

Comments: 10 pages, accepted by MNRAS

Journal ref: MNRAS, 416, 3160-3169 (2011)

A lower mass for the exoplanet WASP-21b

Authors: S. C. C. Barros, D. L. Pollacco, N. P. Gibson, I. D. Howarth, F. P. Keenan, E. K. Simpson, I. Skillen, I. A. Steele, .

Abstract: We present high precision transit observations of the exoplanet WASP-21b, obtained with the RISE instrument mounted on 2.0m Liverpool Telescope. A transit model is fitted, coupled with an MCMC routine to derive accurate system parameters. The two new high precision transits allow to estimate the stellar density directly from the light curve. Our analysis suggests that WASP-21 is evolving off the m… ▽ More We present high precision transit observations of the exoplanet WASP-21b, obtained with the RISE instrument mounted on 2.0m Liverpool Telescope. A transit model is fitted, coupled with an MCMC routine to derive accurate system parameters. The two new high precision transits allow to estimate the stellar density directly from the light curve. Our analysis suggests that WASP-21 is evolving off the main sequence which led to a previous overestimation of the stellar density. Using isochrone interpolation, we find a stellar mass of 0.86 \pm 0.04 Msun which is significantly lower than previously reported (1.01 \pm 0.03 Msun). Consequently, we find a lower planetary mass of $0.27 \pm 0.01 Mjup$. A lower inclination (87.4 \pm 0.3 degrees) is also found for the system than previously reported, resulting in a slightly larger stellar (R_* =1.10 \pm 0.03 Rsun) and planetary radius (R_p = 1.14 \pm 0.04 Rjup). The planet radius suggests a hydrogen/helium composition with no core which strengthens the correlation between planetary density and host star metallicity. A new ephemeris is determined for the system, i.e., t0 =2455084.51974 \pm 0.00020 (HJD) and P=4.3225060 \pm 0.0000031 days. We found no transit timing variations in WASP-21b. △ Less

Submitted 10 June, 2011; originally announced June 2011.

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

The VLT-FLAMES Tarantula Survey I: Introduction and observational overview

Authors: C. J. Evans, W. D. Taylor, V. Henault-Brunet, H. Sana, A. de Koter, S. Simon-Diaz, G. Carraro, T. Bagnoli, N. Bastian, J. M. Bestenlehner, A. Z. Bonanos, E. Bressert, I. Brott, M. A. Campbell, M. Cantiello, J. S. Clark, E. Costa, P. A. Crowther, S. E. de Mink, E. Doran, P. L. Dufton, P. R. Dunstall, K. Friedrich, M. Garcia, M. Gieles , et al. (17 additional authors not shown)

Abstract: The VLT-FLAMES Tarantula Survey (VFTS) is an ESO Large Programme that has obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Here we introduce our scientific motivations and give an overview of the survey targets, including optical and near-infrared photometry and comprehensive details of the data reduction. One of the… ▽ More The VLT-FLAMES Tarantula Survey (VFTS) is an ESO Large Programme that has obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Here we introduce our scientific motivations and give an overview of the survey targets, including optical and near-infrared photometry and comprehensive details of the data reduction. One of the principal objectives was to detect massive binary systems via variations in their radial velocities, thus shaping the multi-epoch observing strategy. Spectral classifications are given for the massive emission-line stars observed by the survey, including the discovery of a new Wolf-Rayet star (VFTS 682, classified as WN5h), 2' to the northeast of R136. To illustrate the diversity of objects encompassed by the survey, we investigate the spectral properties of sixteen targets identified by Gruendl & Chu from Spitzer photometry as candidate young stellar objects or stars with notable mid-infrared excesses. Detailed spectral classification and quantitative analysis of the O- and B-type stars in the VFTS sample, paying particular attention to the effects of rotational mixing and binarity, will be presented in a series of future articles to address fundamental questions in both stellar and cluster evolution. △ Less

Submitted 28 April, 2011; v1 submitted 28 March, 2011; originally announced March 2011.

Comments: Accepted by A&A, 52 pages (main body: 19 pages, supplementary tables: 33 pages), v3: two classifications updated to match a parallel paper

New limb-darkening coefficients and synthetic photometry for model-atmosphere grids at Galactic, LMC, and SMC abundances

Authors: Ian D. Howarth

Abstract: New grids of Atlas9 models have been calculated using revised convection parameters and updated opacity-distribution functions, for chemical compositions intended to be representative of solar, [M/H] = +0.3, +0.5, Large Magellanic Cloud (LMC), and Small Magellanic Cloud (SMC) abundances. The grids cover T(eff) = 3.5-50kK, from log(g) = 5.0 to the effective Eddington limit. Limb-darkening coefficie… ▽ More New grids of Atlas9 models have been calculated using revised convection parameters and updated opacity-distribution functions, for chemical compositions intended to be representative of solar, [M/H] = +0.3, +0.5, Large Magellanic Cloud (LMC), and Small Magellanic Cloud (SMC) abundances. The grids cover T(eff) = 3.5-50kK, from log(g) = 5.0 to the effective Eddington limit. Limb-darkening coefficients and synthetic photometry are presented in the UBVRIJHKLM, uvby, ugriz, WFCAM, Hipparcos/Tycho, and Kepler passbands for these models, and for Castelli's comparable `new-ODF' grids. Flux distributions are given for the new models. The sensitivity of limb-darkening coefficients to the adopted physics is illustrated. △ Less

Submitted 13 November, 2010; v1 submitted 11 November, 2010; originally announced November 2010.

Detection of a magnetic field on HD108: clues to extreme magnetic braking and the Of?p phenomenon

Authors: F. Martins, J. -F. Donati, W. L. F. Marcolino, J. -C. Bouret, G. A. Wade, C. Escolano, I. D. Howarth

Abstract: We report the detection of a magnetic field on the Of?p star HD108. Spectropolarimetric observations conducted in 2007, 2008 and 2009 respectively with NARVAL@TBL and ESPaDOnS@CFHT reveal a clear Zeeman signature in the average Stokes V profile, stable on timescales of days to months and slowly increasing in amplitude on timescales of years. We speculate that this timescale is the same as that on… ▽ More We report the detection of a magnetic field on the Of?p star HD108. Spectropolarimetric observations conducted in 2007, 2008 and 2009 respectively with NARVAL@TBL and ESPaDOnS@CFHT reveal a clear Zeeman signature in the average Stokes V profile, stable on timescales of days to months and slowly increasing in amplitude on timescales of years. We speculate that this timescale is the same as that on which Ha emission is varying and is equal to the rotation period of the star. The corresponding longitudinal magnetic field, measured during each of the three seasons, increases slowly from 100 to 150G, implying that the polar strength of the putatively-dipolar large-scale magnetic field of HD108 is at least 0.5kG and most likely of the order of 1-2 kG. The stellar and wind properties are derived through a quantitative spectroscopic analysis with the code CMFGEN. The effective temperature is difficult to constrain because of the unusually strong HeI4471 and HeI5876 lines. Values in the range 33000-37000 K are preferred. A mass loss rate of about 1e-7 Msun/yr (with a clumping factor f=0.01) and a wind terminal velocity of 2000 km/s are derived. The wind confinement parameter eta_star is larger than 100, implying that the wind of HD108 is magnetically confined. Stochastic short-term variability is observed in the wind-sensitive lines but not in the photospheric lines, excluding the presence of pulsations. Material infall in the confined wind is the most likely origin for lines formed in the inner wind. Wind-clumping also probably causes part of the Ha variability. The projected rotational velocity of HD108 is lower than 50 km/s, consistent with the spectroscopic and photometric variation timescales of a few decades. Overall, HD108 is very similar to the magnetic O star HD191612 except for an even slower rotation. △ Less

Submitted 11 May, 2010; originally announced May 2010.

Comments: 10 pages, 6 figures; MNRAS accepted

A massive runaway star from 30 Doradus

Authors: C. J. Evans, N. R. Walborn, P. A. Crowther, V. Henault-Brunet, D. Massa, W. D. Taylor, I. D. Howarth, H. Sana, D. J. Lennon, J. Th. van Loon

Abstract: We present the first ultraviolet (UV) and multi-epoch optical spectroscopy of 30 Dor 016, a massive O2-type star on the periphery of 30 Doradus in the Large Magellanic Cloud. The UV data were obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope as part of the Servicing Mission Observatory Verification program after Servicing Mission 4, and reveal #016 to have one of the fas… ▽ More We present the first ultraviolet (UV) and multi-epoch optical spectroscopy of 30 Dor 016, a massive O2-type star on the periphery of 30 Doradus in the Large Magellanic Cloud. The UV data were obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope as part of the Servicing Mission Observatory Verification program after Servicing Mission 4, and reveal #016 to have one of the fastest stellar winds known. From analysis of the CIV 1548-51 doublet we find a terminal velocity, v_infty=3450 +/- 50km/s. Optical spectroscopy is from the VLT-FLAMES Tarantula Survey, from which we rule out a massive companion (with 2d<P<1yr) to a confidence of 98%. The radial velocity of #016 is offset from the systemic value by -85km/s, suggesting that the star has traveled the 120pc from the core of 30 Doradus as a runaway, ejected via dynamical interactions. △ Less

Submitted 29 April, 2010; originally announced April 2010.

Comments: 6 pages, ApJL, in press