Researcher profile

R. H. D. Townsend

R. H. D. Townsend contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
35works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2020arXiv

Improved Asymptotic Expressions for the Eigenvalues of Laplace's Tidal Equations

Laplace's tidal equations govern the angular dependence of oscillations in stars when uniform rotation is treated within the so-called traditional approximation. Using a perturbation expansion approach, I derive improved expressions for the eigenvalue associated with these equations, valid in the asymptotic limit of large spin parameter $q$. These expressions have a relative accuracy of order $q^{-3}$ for gravito-inertial modes, and $q^{-1}$ for Rossby and Kelvin modes; the corresponding absolute accuracy is of order $q^{-1}$ for all three mode types. I validate my analysis against numerical calculations, and demonstrate how it can be applied to derive formulae for the periods and eigenfunctions of Rossby modes.

preprint2020arXiv

The Aarhus red giants challenge II. Stellar oscillations in the red giant branch phase

Context. The large quantity of high-quality asteroseismic data that obtained from space-based photometric missions and the accuracy of the resulting frequencies motivate a careful consideration of the accuracy of computed oscillation frequencies of stellar models, when applied as diagnostics of the model properties. Aims. Based on models of red-giant stars that have been independently calculated using different stellar evolution codes, we investigate the extent to which the differences in the model calculation affect the model oscillation frequencies. Methods. For each of the models, which cover four different masses and different evolution stages on the red-giant branch, we computed full sets of low-degree oscillation frequencies using a single pulsation code and, from these frequencies, typical asteroseismic diagnostics. In addition, we carried out preliminary analyses to relate differences in the oscillation properties to the corresponding model differences. Results. In general, the differences in asteroseismic properties between the different models greatly exceed the observational precision of these properties, in particular for the nonradial modes whose mixed acoustic and gravity-wave character makes them sensitive to the structure of the deep stellar interior. In some cases, identifying these differences led to improvements in the final models presented here and in Paper I; here we illustrate particular examples of this. Conclusions. Further improvements in stellar modelling are required in order fully to utilise the observational accuracy to probe intrinsic limitations in the modelling. However, our analysis of the frequency differences and their relation to stellar internal properties provides a striking illustration of the potential of the mixed modes of red-giant stars for the diagnostics of stellar interiors.

preprint2020arXiv

The Contour Method: a new approach to finding modes of non-adiabatic stellar pulsations

The contour method is a new approach to calculating the non-adiabatic pulsation frequencies of stars. These frequencies can be found by solving for the complex roots of a characteristic equation constructed from the linear non-adiabatic stellar pulsation equations. A complex-root solver requires an initial trial frequency for each non adiabatic root. A standard method for obtaining initial trial frequencies is to use a star's adiabatic pulsation frequencies, but this method can fail to converge to non-adiabatic roots, especially as the growth and/or damping rate of the pulsations becomes large. The contour method provides an alternative way for obtaining initial trial frequencies that robustly converges to non-adiabatic roots, even for stellar models with extremely non-adiabatic pulsations and thus large growth/damping rates. We describe the contour method implemented in the GYRE stellar pulsation code and use it to calculate the non-adiabatic pulsation frequencies of $10\,\rm{M_{\odot}}$ and $20\,\rm{M_{\odot}}$ $β$ Cephei star models, and of a $0.9\,\rm{M_{\odot}}$ extreme helium star model.

preprint2020arXiv

The effects of surface fossil magnetic fields on massive star evolution: II. Implementation of magnetic braking in MESA and implications for the evolution of surface rotation in OB stars

The time evolution of angular momentum and surface rotation of massive stars is strongly influenced by fossil magnetic fields via magnetic braking. We present a new module containing a simple, comprehensive implementation of such a field at the surface of a massive star within the Modules for Experiments in Stellar Astrophysics (MESA) software instrument. We test two limiting scenarios for magnetic braking: distributing the angular momentum loss throughout the star in the first case, and restricting the angular momentum loss to a surface reservoir in the second case. We perform a systematic investigation of the rotational evolution using a grid of OB star models with surface magnetic fields ($M_\star=5-60$ M$_\odot$, $Ω/Ω_{\rm crit} =0.2-1.0$, $B_{\rm p} =1-20$ kG). We then employ a representative grid of B-type star models ($M_\star=5, 10, 15$ M$_\odot$, $Ω/Ω_{\rm crit} =0.2 , 0.5, 0.8$, $B_{\rm p} = 1, 3 ,10, 30$ kG) to compare to the results of a recent self-consistent analysis of the sample of known magnetic B-type stars. We infer that magnetic massive stars arrive at the zero age main sequence with a range of rotation rates, rather than with one common value. In particular, some stars are required to have close-to-critical rotation at the ZAMS. However, magnetic braking yields surface rotation rates converging to a common low value, making it difficult to infer the initial rotation rates of evolved, slowly-rotating stars.

preprint2020arXiv

WOCS 5379: Detailed Analysis of the Evolution of a Post-Mass-Transfer Blue Straggler

The blue-straggler binary WOCS 5379 is a member of the old (6-7 Gyr) open cluster NGC 188. WOCS 5379 comprises a blue straggler star with a white dwarf companion in a 120-day eccentric orbit. Combined with the orbital period, this helium white dwarf is evidence of previous mass transfer by a red giant. Detailed models of the system evolution from a progenitor main-sequence binary, including mass transfer, are made using the Modules for Experiments in Stellar Astrophysics (MESA). Both of the progenitor stars are evolved in the simulation. WOCS 5379 is well reproduced with a primary star of initial mass 1.19 $M_{\odot}$, whose core becomes the white dwarf. The secondary star initially is 1.01 $M_{\odot}$. 300 Myr ago, the secondary finished receiving mass from the donor, having moved beyond the NGC 188 turnoff as a 1.20 $M_{\odot}$ blue straggler. The successful model has a mass transfer efficiency of 22\%. This non-conservative mass transfer is key to expanding the orbit fast enough to permit stable mass transfer. Even so, the mass transfer begins with a short unstable phase, during which half of the accreted mass is transferred. With increasing mass, the secondary evolves from a radiative core to a convective core. The final blue straggler interior is remarkably similar to a 2.1 Gyr-old 1.21 $M_{\odot}$ main-sequence star at the same location in the HR diagram. The white dwarf effective temperature is also reproduced, but the modeled white dwarf mass of 0.33 $M_{\odot}$ is smaller than the measured mass of 0.42 $M_{\odot}$.

preprint2019arXiv

The Aarhus Red Giants Challenge I: Stellar structures in the red giant branch phase

(Abridged). We introduce the Aarhus Red Giants Challenge, a series of detailed comparisons between widely used stellar evolution and oscillation codes aiming at establishing the minimum level of uncertainties in properties of red giants arising solely from numerical implementations. Using 9 state-of-the-art stellar evolution codes, we defined a set of input physics and physical constants for our calculations and calibrated the convective efficiency to a specific point on the main sequence. We produced evolutionary tracks and stellar structure models at fixed radius along the red-giant branch for masses of 1.0 M$_\odot$, 1.5 M$_\odot$, 2.0 M$_\odot$, and 2.5 M$_\odot$, and compared the predicted stellar properties. Once models have been calibrated on the main sequence we find a residual spread in the predicted effective temperatures across all codes of ~20 K at solar radius and ~30-40 K in the RGB regardless of the considered stellar mass. The predicted ages show variations of 2-5% (increasing with stellar mass) which we track down to differences in the numerical implementation of energy generation. The luminosity of the RGB-bump shows a spread of about 10% for the considered codes, which translates into magnitude differences of ~0.1 mag in the optical V-band. We also compare the predicted [C/N] abundance ratio and found a spread of 0.1 dex or more for all considered masses. Our comparisons show that differences at the level of a few percent still remain in evolutionary calculations of red giants branch stars despite the use of the same input physics. These are mostly due to differences in the energy generation routines and interpolation across opacities, and call for further investigations on these matters in the context of using properties of red giants as benchmarks for astrophysical studies.

preprint2016arXiv

The Variable Line Width of Achernar

Spectroscopic observations of Achernar over the past decades, have shown the photospheric line width, as measured by the rotational parameter $v \sin i$, to vary in correlation with the emission activity. Here we present new observations, covering the most recent activity phase, and further archival data collected from the archives. The $v \sin i$ variation is confirmed. On the basis of the available data it cannot be decided with certainty whether the increased line width precedes the emission activity, i.e. is a signature of the ejection mechanism, or postdates is, which would make it a signature of re-accretion of some of the disk-material. However, the observed evidence leans towards the re-accretion hypothesis. Two further stars showing the effect of variable line width in correlation with emission activity, namely 66 Oph and $π$ Aqr, are presented as well.

preprint2015arXiv

Long term variability of Cygnus X-1: VII. Orbital variability of the focussed wind in Cyg X-1 / HDE 226868 system

Binary systems with an accreting compact object are a unique chance to investigate the strong, clumpy, line-driven winds of early type supergiants by using the compact object's X-rays to probe the wind structure. We analyze the two-component wind of HDE 226868, the O9.7Iab giant companion of the black hole Cyg X-1 using 4.77 Ms of RXTE observations of the system taken over the course of 16 years. Absorption changes strongly over the 5.6 d binary orbit, but also shows a large scatter at a given orbital phase, especially at superior conjunction. The orbital variability is most prominent when the black hole is in the hard X-ray state. Our data are poorer for the intermediate and soft state, but show signs for orbital variability of the absorption column in the intermediate state. We quantitatively compare the data in the hard state to a toy model of a focussed Castor-Abbott-Klein-wind: as it does not incorporate clumping, the model does not describe the observations well. A qualitative comparison to a simplified simulation of clumpy winds with spherical clumps shows good agreement in the distribution of the equivalent hydrogen column density for models with a porosity length on the order of the stellar radius at inferior conjunction; we conjecture that the deviations between data and model at superior conjunction could be either due to lack of a focussed wind component in the model or a more complicated clump structure.

preprint2015arXiv

Revisiting the Rigidly Rotating Magnetosphere model for $σ$ Ori E - II. Magnetic Doppler imaging, arbitrary field RRM, and light variability

The initial success of the Rigidly Rotating Magnetosphere (RRM) model application to the B2Vp star sigma OriE by Townsend, Owocki & Groote (2005) triggered a renewed era of observational monitoring of this archetypal object. We utilize high-resolution spectropolarimetry and the magnetic Doppler imaging (MDI) technique to simultaneously determine the magnetic configuration, which is predominately dipolar, with a polar strength Bd = 7.3-7.8 kG and a smaller non-axisymmetric quadrupolar contribution, as well as the surface distribution of abundance of He, Fe, C, and Si. We describe a revised RRM model that now accepts an arbitrary surface magnetic field configuration, with the field topology from the MDI models used as input. The resulting synthetic Ha emission and broadband photometric observations generally agree with observations, however, several features are poorly fit. To explore the possibility of a photospheric contribution to the observed photometric variability, the MDI abundance maps were used to compute a synthetic photospheric light curve to determine the effect of the surface inhomogeneities. Including the computed photospheric brightness modulation fails to improve the agreement between the observed and computed photometry. We conclude that the discrepancies cannot be explained as an effect of inhomogeneous surface abundance. Analysis of the UV light variability shows good agreement between observed variability and computed light curves, supporting the accuracy of the photospheric light variation calculation. We thus conclude that significant additional physics is necessary for the RRM model to acceptably reproduce observations of not only sigma Ori E, but also other similar stars with significant stellar wind-magnetic field interactions.

preprint2015arXiv

Structural glitches near the cores of red giants revealed by oscillations in g-mode period spacings from stellar models

With recent advances in asteroseismology it is now possible to peer into the cores of red giants, potentially providing a way to study processes such as nuclear burning and mixing through their imprint as sharp structural variations -- glitches -- in the stellar cores. Here we show how such core glitches can affect the oscillations we observe in red giants. We derive an analytical expression describing the expected frequency pattern in the presence of a glitch. This formulation also accounts for the coupling between acoustic and gravity waves. From an extensive set of canonical stellar models we find glitch-induced variation in the period spacing and inertia of non-radial modes during several phases of red-giant evolution. Significant changes are seen in the appearance of mode amplitude and frequency patterns in asteroseismic diagrams such as the power spectrum and the échelle diagram. Interestingly, along the red-giant branch glitch-induced variation occurs only at the luminosity bump, potentially providing a direct seismic indicator of stars in that particular evolution stage. Similarly, we find the variation at only certain post-helium-ignition evolution stages, namely, in the early phases of helium-core burning and at the beginning of helium-shell burning signifying the asymptotic-giant-branch bump. Based on our results, we note that assuming stars to be glitch-free, while they are not, can result in an incorrect estimate of the period spacing. We further note that including diffusion and mixing beyond classical Schwarzschild, could affect the characteristics of the glitches, potentially providing a way to study these physical processes.

preprint2015arXiv

The magnetic field and spectral variability of the He-weak star HR 2949

We analyze a high resolution spectropolarimetric dataset collected for the He-weak B3p IV star HR 2949. The Zeeman effect is visible in the circularly polarized component of numerous spectral lines. The longitudinal magnetic field varies between approximately $-650$ and $+150$ G. The polar strength of the surface magnetic dipole is calculated to be 2.4$^{+0.3}_{-0.2}$ kG. The star has strong overabundances of Fe-peak elements, along with extremely strong overabundances of rare-earth elements; however, He, Al, and S are underabundant. This implies that HR 2949 is a chemically peculiar star. Variability is seen in all photospheric lines, likely due to abundance patches as seen in many Ap/Bp stars. Longitudinal magnetic field variations measured from different spectral lines yield different results, likely a consequence of uneven sampling of the photospheric magnetic field by the abundance patches. Analysis of photometric and spectroscopic data for both HR 2949 and its companion star, HR 2948, suggests a revision of HR 2949's fundamental parameters: in particular, it is somewhat larger, hotter, and more luminous than previously believed. There is no evidence of optical or ultraviolet emission originating in HR 2949's magnetosphere, despite its moderately strong magnetic field and relatively rapid rotation; however, when calculated using theoretical and empirical boundaries on the initial rotational velocity, the spindown age is compatible with the stellar age. With the extensive phase coverage presented here, HR 2949 will make an excellent subject for Zeeman Doppler Imaging.

preprint2013arXiv

MOST Observations of sigma Ori E: Challenging the Centrifugal Breakout Narrative

We present results from three weeks' photometric monitoring of the magnetic helium-strong star sigma Ori E using the MOST microsatellite. The star's light curve is dominated by twice-per-rotation eclipse-like dimmings arising when magnetospheric clouds transit across and occult the stellar disk. However, no evidence is found for any abrupt centrifugal breakout of plasma from the magnetosphere, either in the residual flux or in the depths of the light minima. Motivated by this finding we compare the observationally inferred magnetospheric mass against that predicted by a breakout analysis. The large discrepancy between the values leads us to argue that centrifugal breakout does not play a significant role in establishing the magnetospheric mass budget of sigma Ori E.

preprint2013arXiv

Variable rotational line broadening in the Be star Achernar

The main theoretical problem for the formation of a Keplerian disk around Be stars is how to supply angular momentum from the star to the disk, even more so since Be stars probably rotate somewhat sub-critically. For instance, nonradial pulsation may transport angular momentum to the stellar surface until (part of) this excess supports the disk formation/replenishment. The nearby Be star Achernar is presently building a new disk and offers an excellent opportunity to observe this process from relatively close-up. Spectra from various sources and epochs are scrutinized to identify the salient stellar parameters characterizing the disk life cycle as defined by Hαemission. Variable strength of the non-radial pulsation is confirmed, but does not affect the further results. For the first time it is demonstrated that the photospheric line width does vary in a Be star, by as much as Δv sin i \lesssim 35kms^{-1}. However, contrary to assumptions in which a photospheric spin-up accumulates during the diskless phase and then is released into the disk as it is fed, the apparent photospheric spin-up is positively correlated with the appearance of Hαline emission: The photospheric line widths and circumstellar emission increase together, and the apparent stellar rotation declines to the value at quiescence after the Hαline emission becomes undetectable.

preprint2012arXiv

A Magnetic Confinement vs. Rotation Classification of Massive-Star Magnetospheres

Building on results from the Magnetism in Massive Stars (MiMeS) project, this paper shows how a two-parameter classification of massive-star magnetospheres in terms of the magnetic wind confinement (which sets the Alfvén radius RA) and stellar rotation (which sets the Kepler co-rotation radius RK) provides a useful organisation of both observational signatures and theoretical predictions. We compile the first comprehensive study of inferred and observed values for relevant stellar and magnetic parameters of 64 confirmed magnetic OB stars with Teff > 16 kK. Using these parameters, we locate the stars in the magnetic confinement-rotation diagram, a log-log plot of RK vs. RA. This diagram can be subdivided into regimes of centrifugal magnetospheres (CM), with RA > RK, vs. dynamical magnetospheres (DM), with RK > RA. We show how key observational diagnostics, like the presence and characteristics of Halpha emission, depend on a star's position within the diagram, as well as other parameters, especially the expected wind mass-loss rates. In particular, we identify two distinct populations of magnetic stars with Halpha emission: namely, slowly rotating O-type stars with narrow emission consistent with a DM, and more rapidly rotating B-type stars with broader emission associated with a CM. For O-type stars, the high mass-loss rates are sufficient to accumulate enough material for line emission even within the relatively short free-fall timescale associated with a DM: this high mass-loss rate also leads to a rapid magnetic spindown of the stellar rotation. For the B-type stars, the longer confinement of a CM is required to accumulate sufficient emitting material from their relatively weak winds, which also lead to much longer spindown timescales. [abbreviated]

preprint2012arXiv

Basic parameters and properties of the rapidly rotating magnetic helium strong B star HR7355

The spectral and magnetic properties and variability of the B2Vnp emission-line magnetosphere star HR7355 were analyzed. The object rotates at almost 90% of the critical value, meaning it is a magnetic star for which oblateness and gravity darkening effects cannot be ignored any longer. A detailed modeling of the photospheric parameters indicate that the star is significantly cooler than suggested by the B2 spectral type, with T_eff=17500K atypically cool for a star with a helium enriched surface. The spectroscopic variability of helium and metal lines due to the photospheric abundance pattern is far more complex than a largely dipolar, oblique magnetic field of about 11 to 12kG may suggest. Doppler imaging shows that globally the most He enriched areas coincide with the magnetic poles and metal enriched areas with the magnetic equator. While most of the stellar surface is helium enriched with respect to the solar value, some isolated patches are depleted. The stellar wind in the circumstellar environment is governed by the magnetic field, i.e. the stellar magnetosphere is rigidly corotating with the star. The magnetosphere of HR7355 is similar to the well known σOri E: the gas trapped in the magnetospheric clouds is fairly dense, and at the limit to being optically thick in the hydrogen emission. Apart from a different magnetic obliquity, HR7355 and the more recently identified HR5907 have virtually identical stellar and magnetic parameters.

preprint2012arXiv

Solving the Puzzle of the Massive Star System Theta 2 Orionis A

The young O9.5 V spectroscopic binary Theta 2 Ori A shows moderately hard X-ray emission and relatively narrow X-ray lines, suggesting that it may be a Magnetically Confined Wind Shock (MCWS) source, similar to its more massive analogue Theta 1 Ori C. X-ray flares occurring near periastron led to the suggestion that the flares are produced via magnetic reconnection as magnetospheres on both components of the Theta 2 Ori A binary interact at closest approach. We use a series of high-resolution spectropolarimetric observations of Theta 2 Ori A to place an upper limit on the magnetic field strength of 135 G (95% credible region). Such a weak dipole field would not produce magnetic confinement, or a large magnetosphere. A sub-pixel analysis of the Chandra ACIS images of Theta 2 Ori A obtained during quiescence and flaring show that the hard, flaring X-rays are offset from the soft, quiescent emission by 0.4 arcsec. If the soft emission is associated with the A1/A2 spectroscopic binary, the offset and position angle of the hard, flaring source place it at the location of the intermediate-mass A3 companion, discovered via speckle interferometry. The spectropolarimetric and X-ray results taken together point to the A3 companion, not the massive A1/A2 binary, as the source of hard, flaring X-ray emission. We also discuss a similar analysis performed for the magnetic Bp star Sigma Ori E. We find a similar origin for its X-ray flaring.

preprint2011arXiv

HR 5907: Discovery of the most rapidly rotating magnetic B-type star by the MiMeS Collaboration

We report the discovery and analysis of a very strong magnetic field in the rapidly rotating early B-type star HR 5907, based on observations obtained as part of the Magnetism in Massive Stars (MiMeS) project. We infer a rotation period of 0.508276 +0.000015/-0.000012 d from photometric and Hα EW measurements, making this the shortest period, non-degenerate, magnetic massive star known to date. From the comparison of IUE UV and optical spectroscopy with LTE BRUCE/KYLIE models we find a solid-angle integrated, uniform black-body temperature of 17 000 \pm 1000 K, a projected rotational velocity of 290 \pm 10 km/s, an equatorial radius of 3.1 \pm 0.2 R_sun, a stellar mass of 5.5 \pm 0.5 M_sun, and an inclination angle of the rotation axis to our line-of-sight of 70 \pm 10\circ. Our measurements of the longitudinal magnetic field, which vary between -500 and -2000 G, phase coherently with the rotation period and imply a surface dipole field strength of \sim15.7 kG. On the other hand, from fits to mean Least-Squares Deconvolved Stokes V line profiles we infer a dipole field strength of \sim10.4 kG. This disagreement may result from a magnetic configuration more complex than our model, and/or from the non-uniform helium surface abundance distribution. In either case we obtain a magnetic obliquity nearly aligned with the rotation axis (β = 7+2/-1\circ). Our optical spectroscopy also shows weak variability in carbon, silicon and nitrogen lines. The emission variability in hydrogen Balmer and Paschen lines indicates the presence of a dense, highly structured magnetosphere, interpreted as a centrifugally supported, magnetically confined circumstellar disk.

preprint2011arXiv

Revisiting the Rigidly Rotating Magnetosphere model for sigma Ori E. I. Observations and Data Analysis

We have obtained 18 new high-resolution spectropolarimetric observations of the B2Vp star sigma Ori E with both the Narval and ESPaDOnS spectropolarimeters. The aim of these observations is to test, with modern data, the assumptions of the Rigidly Rotating Magnetosphere (RRM) model of Townsend & Owocki (2005), applied to the specific case of sigma Ori E by Townsend et al. (2005). This model includes a substantially offset dipole magnetic field configuration, and approximately reproduces previous observational variations in longitudinal field strength, photometric brightness, and Halpha emission. We analyze new spectroscopy, including H I, He I, C II, Si III and Fe III lines, confirming the diversity of variability in photospheric lines, as well as the double S-wave variation of circumstellar hydrogen. Using the multiline analysis method of Least-Squares Deconvolution (LSD), new, more precise longitudinal magnetic field measurements reveal a substantial variance between the shapes of the observed and RRM model time-varying field. The phase resolved Stokes V profiles of He I 5876 A and 6678 A lines are fit poorly by synthetic profiles computed from the magnetic topology assumed by Townsend et al. (2005). These results challenge the offset dipole field configuration assumed in the application of the RRM model to sigma Ori E, and indicate that future models of its magnetic field should also include complex, higher-order components.

preprint2010arXiv

Discovery of a strong magnetic field in the rapidly rotating B2Vn star HR 7355

We report on the detection of a strong, organized magnetic field in the helium-variable early B-type star HR 7355 using spectropolarimetric data obtained with ESPaDOnS on CFHT by the MiMeS large program. We also present results from new V-band differential photometry obtained with the CTIO 0.9m telescope. We investigate the longitudinal field, using a technique called Least-Squares Deconvolution (LSD), and the rotational period of HR 7355. These new observations strongly support the proposal that HR 7355 harbors a structured magnetosphere similar to that in the prototypical helium-strong star, sigma Ori E.

preprint2010arXiv

Discovery of a strong magnetic field in the rapidly rotating B2Vn star HR 7355

We report the detection of a strong, organized magnetic field in the helium-variable early B-type star HR 7355 using spectropolarimetric data obtained with ESPaDOnS on the 3.6-m Canada-France-Hawaii Telescope within the context of the Magnetism in Massive Stars (MiMeS) Large Program. HR 7355 is both the most rapidly rotating known main-sequence magnetic star and the most rapidly rotating helium-strong star, with $v \sin i$ = 300 $\pm$ 15 km s$^{-1}$ and a rotational period of 0.5214404 $\pm$ 0.0000006 days. We have modeled our eight longitudinal magnetic field measurements assuming an oblique dipole magnetic field. Constraining the inclination of the rotation axis to be between $38^{\circ}$ and $86^{\circ}$, we find the magnetic obliquity angle to be between $30^{\circ}$ and $85^{\circ}$, and the polar strength of the magnetic field at the stellar surface to be between 13-17 kG. The photometric light curve constructed from HIPPARCOS archival data and new CTIO measurements shows two minima separated by 0.5 in rotational phase and occurring 0.25 cycles before/after the magnetic extrema. This photometric behavior coupled with previously-reported variable emission of the H$α$ line (which we confirm) strongly supports the proposal that HR 7355 harbors a structured magnetosphere similar to that in the prototypical helium-strong star, $σ$ Ori E.

preprint2010arXiv

Discovery of Rotational Braking in the Magnetic Helium-Strong Star Sigma Orionis E

We present new U-band photometry of the magnetic Helium-strong star Sigma Ori E, obtained over 2004-2009 using the SMARTS 0.9-m telescope at Cerro Tololo Inter-American Observatory. When combined with historical measurements, these data constrain the evolution of the star's 1.19 d rotation period over the past three decades. We are able to rule out a constant period at the p_null = 0.05% level, and instead find that the data are well described (p_null = 99.3%) by a period increasing linearly at a rate of 77 ms per year. This corresponds to a characteristic spin-down time of 1.34 Myr, in good agreement with theoretical predictions based on magnetohydrodynamical simulations of angular momentum loss from magnetic massive stars. We therefore conclude that the observations are consistent with Sigma Ori E undergoing rotational braking due to its magnetized line-driven wind.

preprint2010arXiv

Discovery of the most rapidly-rotating, non-degenerate, magnetic massive star by the MiMeS collaboration

We discuss the recent detection of a strong, organized magnetic field in the bright, broad-line B2V star, HD 142184, using the ESPaDOnS spectropolarimeter on the CFHT as part of the Magnetism in Massive Stars (MiMeS) survey. We find a rotational period of 0.50833 days, making it the fastest-rotating, non-degenerate magnetic star ever detected. Like the previous rapid-rotation record holder HR 7355 (also discovered by MiMeS: Oksala et al. 2010, Rivinius et al. 2010), this star shows emission line variability that is diagnostic of a structured magnetosphere.

preprint2010arXiv

Magnetic field detection in the B2Vn star HR 7355

The B2Vn star HR 7355 is found to be a He-rich magnetic star. Spectropolarimetric data were obtained with FORS1 at UT2 on Paranal observatory to measure the disk-averaged longitudinal magnetic field at various phases of the presumed 0.52 d cycle. A variable magnetic field with strengths between B_z = -2200 and +3200G was found, with confidence limits of 100 to 130G. The field topology is that of an oblique dipole, while the star itself is seen about equator-on. In the intensity spectra the HeI-lines show the typical equivalent width variability of He-strong stars, usually attributed to surface abundance spots. The amplitudes of the equivalent width variability of the HeI lines are extraordinarily strong compared to other cases. These results not only put HR 7355 unambiguously among the early-type magnetic stars, but confirm its outstanding nature: With v sin i = 320 km/s the parameter space in which He-strong stars are known to exist has doubled in terms of rotational velocity.

preprint2010arXiv

Rotationally-modulated g-modes in the rapidly-rotating delta Scuti star Rasalhague (alpha Ophiuchi)

Despite a century of remarkable progress in understanding stellar interiors, we know surprisingly little about the inner workings of stars spinning near their critical limit. New interferometric imaging of these so-called ``rapid rotators'' combined with breakthroughs in asteroseismology promise to lift this veil and probe the strongly latitude-dependent photospheric characteristics and even reveal the internal angular momentum distribution of these luminous objects. Here, we report the first high precision photometry on the low-amplitude delta cuti variable star Rasalhague (alpha Oph, A5IV, 2.18 Msun, omega/omega_c~0.88) based on 30 continuous days of monitoring using the MOST satellite. We have identified 57+/-1 distinct pulsation modes above a stochastic granulation spectrum with a cutoff of ~26 cycles per day. Remarkably, we have also discovered that the fast rotation period of 14.5~hours modulates low-frequency modes (1-10 day periods) that we identify as a rich family of g-modes (|m| up to 7). The spacing of the g-modes is surprisingly linear considering Coriolis forces are expected to strongly distort the mode spectrum, suggesting we are seeing prograde ``equatorial Kelvin'' waves (modes l=m). We emphasize the unique aspects of Rasalhague motivating future detailed asteroseismic modeling -- a source with a precisely measured parallax distance, photospheric oblateness, latitude temperature structure, and whose low-mass companion provides an astrometric orbit for precise mass determinations.

preprint2010arXiv

The (B0+?)+O6 system FN CMa: A case for tidal-pulsational interaction?

FN CMa is visually double with a separation of about 0.6arcsec. Sixty high-cadence VLT/UVES spectra permit the A and B components to be disentangled, as the relative contribution of each star to the total light entering the spectrograph fluctuates between exposures due to changes in seeing. Component A exhibits rapid line-profile variations, leading us to attribute the photometric variability seen by HIPPARCOS (with a derived P=0.08866d) to this component. From a total of 122 archival and new echelle spectra it is shown that component A is an SB1 binary with an orbital period of 117.55 days. The eccentricity of 0.6 may result in tidal modulation of the pulsation(s) of component Aa.

preprint2010arXiv

The rapid magnetic rotator HR7355

For early type magnetic stars slow, at most moderate rotational velocities have been considered an observational fact. The detection of a multi-kilogauss magnetic field in the B2Vpn star with P approx. 0.52d and v sin approx. 300km/s has brought down this narrative. We have obtained more than 100 high-resolution, high-S/N echelle spectra in 2009. These spectra provide the most detailed description of the variability of any He-strong star to date. The circumstellar environment is dominated by a rotationally locked magnetosphere out to several stellar radii, causing hydrogen emission. The photosphere is characterized by surface chemical abundance inhomogeneities, with much stronger amplitudes, at least for helium, than slower rotating stars like σOri E. The highly complex rotational line profile modulations of metal lines are probably a consequence the equatorial gravity darkening of HR 7355, and thus may offer an independent measurement of the von Zeipel parameter β.

preprint2009arXiv

Discovery of a magnetic field in the O9 sub-giant star HD 57682 by the MiMeS Collaboration

We report the detection of a strong, organised magnetic field in the O9IV star HD 57682, using spectropolarimetric observations obtained with ESPaDOnS at the 3.6-m Canada-France-Hawaii Telescope within the context of the Magnetism in Massive Stars (MiMeS) Large Program. From the fitting of our spectra using NLTE model atmospheres we determined that HD 57682 is a $17^{+19}_{-9}$ M$_{\odot}$ star with a radius of $7.0^{+2.4}_{-1.8}$ R$_\odot$, and a relatively low mass-loss rate of $1.4^{+3.1}_{-0.95}\times10^{-9}$ M$_{\odot}$ yr$^{-1}$. The photospheric absorption lines are narrow, and we use the Fourier transform technique to infer $v\sin i=15\pm3$ km s$^{-1}$. This $v\sin i$ implies a maximum rotational period of 31.5 d, a value qualitatively consistent with the observed variability of the optical absorption and emission lines, as well as the Stokes $V$ profiles and longitudinal field. Using a Bayesian analysis of the velocity-resolved Stokes $V$ profiles to infer the magnetic field characteristics, we tentatively derive a dipole field strength of $1680^{+134}_{-356}$ G. The derived field strength and wind characteristics imply a wind that is strongly confined by the magnetic field.

preprint2009arXiv

Gamma-ray variability from wind clumping in HMXBs with jets

In the subclass of high-mass X-ray binaries known as "microquasars", relativistic hadrons in the jets launched by the compact object can interact with cold protons from the star's radiatively driven wind, producing pions that then quickly decay into gamma rays. Since the resulting gamma-ray emissivity depends on the target density, the detection of rapid variability in microquasars with GLAST and the new generation of Cherenkov imaging arrays could be used to probe the clumped structure of the stellar wind. We show here that the fluctuation in gamma rays can be modeled using a "porosity length" formalism, usually applied to characterize clumping effects. In particular, for a porosity length defined by h=l/f, i.e. as the ratio of the characteristic size l of clumps to their volume filling factor f, we find that the relative fluctuation in gamma-ray emission in a binary with orbital separation a scales as sqrt(h/pi a) in the "thin-jet" limit, and is reduced by a factor 1/sqrt(1 + phi a/(2 l)) for a jet with a finite opening angle phi. For a thin jet and quite moderate porosity length h ~ 0.03 a, this implies a ca. 10 % variation in the gamma-ray emission. Moreover, the illumination of individual large clumps might result in isolated flares, as has been recently observed in some massive gamma-ray binaries.

preprint2006arXiv

Centrifugal Breakout of Magnetically Confined Line-Driven Stellar Winds

We present 2D MHD simulations of the radiatively driven outflow from a rotating hot star with a dipole magnetic field aligned with the star's rotation axis. We focus primarily on a model with moderately rapid rotation (half the critical value), and also a large magnetic confinement parameter, $η_{\ast} \equiv B_{\ast}^2 R_{\ast}^{2} / \dot{M} V_{\infty} = 600$. The magnetic field channels and torques the wind outflow into an equatorial, rigidly rotating disk extending from near the Kepler corotation radius outwards. Even with fine-tuning at lower magnetic confinement, none of the MHD models produce a stable Keplerian disk. Instead, material below the Kepler radius falls back on to the stellar surface, while the strong centrifugal force on material beyond the corotation escape radius stretches the magnetic loops outwards, leading to episodic breakout of mass when the field reconnects. The associated dissipation of magnetic energy heats material to temperatures of nearly $10^{8}$K, high enough to emit hard (several keV) X-rays. Such \emph{centrifugal mass ejection} represents a novel mechanism for driving magnetic reconnection, and seems a very promising basis for modeling X-ray flares recently observed in rotating magnetic Bp stars like $σ$ Ori E.

preprint2005arXiv

Influence of the Coriolis force on the instability of slowly pulsating B stars

This paper explores the effect of rotation on the kappa-mechanism instability of slowly pulsating B stars. A new nonadiabatic code, that adopts the so-called `traditional approximation&#39; to treat the Coriolis force, is used to investigate the influence exerted by rotation over the stability of stellar models covering the mass range 2.5 M_sun <= M_* <= 13.0 M_sun. The principal finding is that, for all modes considered apart from the prograde sectoral class, rotation shifts the kappa-mechanism instability toward higher luminosities and effective temperatures; these shifts are accompanied by broadenings in the extent of instability strips. Such behaviour is traced to the shortening of mode periods under the action of the Coriolis force. Instability strips associated with prograde sectoral modes behave rather differently, being shifted to marginally lower luminosities and effective temperatures under the influence of rotation. The implications of these results are discussed in the context of the observational scarcity of pulsation in B-type stars having significant rotation; various scenarios are explored to explain the apparent dichotomy between theory and observations. Furthermore, the possible significance of the findings to Be stars is briefly examined.

preprint2005arXiv

On the evidence for discs around Blue Straggler stars

Recent observations of blue stragglers by De Marco et al. (2004) have revealed continuum deficits on the blue side of the Balmer discontinuity, leading these authors to infer the presence of discs around the stars. This intriguing possibility may throw light on aspects of the mechanisms responsible for at least some of these objects; current theories of blue straggler formation invoke stellar collisions or interacting binaries, both of which appear capable of forming a circumstellar disc. However, by synthesizing photospheric spectra for models of rotating blue stragglers, we demonstrate that the Balmer jump enhancements can be wholly attributed to the influence of oblateness and gravity darkening on the formation of the continuum. Therefore, we are led to conclude that the observations of De Marco et al. can be ascribed a more prosaic explanation, that of rapid stellar rotation arising from the merger/interaction formation process.

preprint2005arXiv

The Rigidly Rotating Magnetosphere of Sigma Ori E

We attempt to characterize the observed variability of the magnetic helium-strong star sigma Ori E in terms of a recently developed rigidly rotating magnetosphere model. This model predicts the accumulation of circumstellar plasma in two co-rotating clouds, situated in magnetohydrostatic equilibrium at the intersection between magnetic and rotational equators. We find that the model can reproduce well the periodic modulations observed in the star&#39;s light curve, H alpha emission-line profile, and longitudinal field strength, confirming that it furnishes an essentially correct, quantitative description of the star&#39;s magnetically controlled circumstellar environment.

preprint2004arXiv

A Rigidly Rotating Magnetosphere Model for Circumstellar Emission from Magnetic OB Stars

We present a semi-analytical approach for modeling circumstellar emission from rotating hot stars with a strong dipole magnetic field tilted at an arbitrary angle to the rotation axis. By assuming the rigid-field limit in which material driven (e.g., in a wind outflow) from the star is forced to remain in strict rigid-body co-rotation, we are able to solve for the effective centrifugal-plus-gravitational potential along each field line, and thereby identify the location of potential minima where material is prone to accumulate. Applying basic scalings for the surface mass flux of a radiatively driven stellar wind, we calculate the circumstellar density distribution that obtains once ejected plasma settles into hydrostatic stratification along field lines. The resulting accumulation surface resembles a rigidly rotating, warped disk, tilted such that its average surface normal lies between the rotation and magnetic axes. Using a simple model of the plasma emissivity, we calculate time-resolved synthetic line spectra for the disk. Initial comparisons show an encouraging level of correspondence with the observed rotational phase variations of Balmer-line emission profiles from magnetic Bp stars like sigma Ori E.