Showing 1–50 of 97 results for author: Townsend, R H D

The impact of radiative levitation on mode excitation of main-sequence B-type pulsators

Authors: R. Rehm, J. S. G. Mombarg, C. Aerts, M. Michielsen, S. Burssens, R. H. D. Townsend

Abstract: Numerical computations of stellar oscillations for models representative of B-type stars predict fewer modes to be excited than observations reveal from modern space-based photometric data. One shortcoming of state-of-the-art evolution models of B-type stars that may cause a lack of excited modes is the absence of microscopic diffusion in most such models. We investigate whether the inclusion of m… ▽ More Numerical computations of stellar oscillations for models representative of B-type stars predict fewer modes to be excited than observations reveal from modern space-based photometric data. One shortcoming of state-of-the-art evolution models of B-type stars that may cause a lack of excited modes is the absence of microscopic diffusion in most such models. We investigate whether the inclusion of microscopic diffusion in stellar models of B-type stars, notably radiative levitation experienced by isotopes, leads to extra mode driving by the opacity mechanism compared to the case of models that do not include microscopic diffusion. We consider the case of slowly to moderately rotating stars and use non-rotating equilibrium models, while we account for (uniform) rotation in the computations of the pulsation frequencies. We calculate 1D stellar models with and without microscopic diffusion and examine the effect of radiative levitation on mode excitation, for both low-radial order pressure and gravity modes and for high-radial order gravity modes. We find systematically more modes to be excited for the stellar models including microscopic diffusion compared to those without it, in agreement with observational findings of pulsating B-type dwarfs. Furthermore, the models with microscopic diffusion predict that excited modes occur earlier on in the evolution compared to modes without it. In order to maintain realistic surface abundances during the main sequence, we include macroscopic envelope mixing by internal gravity waves. While radiative levitation has so far largely been neglected in stellar evolution computations of B-type stars for computational convenience, it impacts mode excitation predictions for stellar models of such stars. We conclude that the process of radiative levitation is able to reduce the discrepancy between predicted and observed excited pulsation modes in B-type stars. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: Paper accepted for publication in Astronomy & Astrophysics (10 pages, 11 figures). Full abstract in PDF

The intermittently-resonant coevolution of migrating planets and their pulsating stars

Authors: Jared Bryan, Julien de Wit, Meng Sun, Zoe L. de Beurs, Richard H. D. Townsend

Abstract: Hot Jupiters are expected to form far from their host star and move toward close-in, circular orbits via a smooth, monotonic decay due to mild and constant tidal dissipation. Yet, three systems have recently been found exhibiting planet-induced stellar pulsations suggesting unexpectedly strong tidal interactions. Here we combine stellar evolution and tide models to show that dynamical tides raised… ▽ More Hot Jupiters are expected to form far from their host star and move toward close-in, circular orbits via a smooth, monotonic decay due to mild and constant tidal dissipation. Yet, three systems have recently been found exhibiting planet-induced stellar pulsations suggesting unexpectedly strong tidal interactions. Here we combine stellar evolution and tide models to show that dynamical tides raised by eccentric gas giants can give rise to chains of resonance locks with multiple modes, enriching the dynamics seen in single-mode resonance locking of circularized systems. These series of resonance locks yield orders-of-magnitude larger changes in eccentricity and harmonic pulsations relative to those expected from a single episode of resonance locking or nonresonant tidal interactions. Resonances become more frequent as a star evolves off the main sequence providing an alternative explanation to the origin of some stellar pulsators and yielding the concept of "dormant migrating giants". Evolution trajectories are characterized by competing episodes of inward/outward migration and spin-up/-down of the star which are sensitive to the system parameters, revealing a new challenge in modeling migration paths and in contextualizing the observed populations of giant exoplanets and stellar binaries. This sensitivity however offers a new window to constrain the stellar properties of planetary hosts via tidal asteroseismology. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Comments: Revised after positive initial review report

Modelling time-dependent convective penetration in 1D stellar evolution

Authors: Cole Johnston, Mathias Michielsen, Evan H. Anders, Mathieu Renzo, Matteo Cantiello, P. Marchant, Jared A. Goldberg, Richard H. D. Townsend, Gautham Sabhahit, Adam S. Jermyn

Abstract: 1D stellar evolution calculations produce uncertain predictions for quantities like the age, core mass, core compactness, and nucleo-synthetic yields; a key source of uncertainty is the modeling of interfaces between regions that are convectively stable and those that are not. Theoretical and numerical work has demonstrated that there should be numerous processes adjacent to the convective boundar… ▽ More 1D stellar evolution calculations produce uncertain predictions for quantities like the age, core mass, core compactness, and nucleo-synthetic yields; a key source of uncertainty is the modeling of interfaces between regions that are convectively stable and those that are not. Theoretical and numerical work has demonstrated that there should be numerous processes adjacent to the convective boundary that induce chemical and angular momentum transport, as well as modify the thermal structure of the star. One such process is called convective penetration, wherein vigorous convection extends beyond the nominal convective boundary and alters both the composition and thermal structure. In this work, we incorporate the process of convective penetration in stellar evolution calculations using the stellar evolution software instrument mesa. We implement convective penetration according to the description presented by Anders et al. (2022a) to calculate a grid of models from the pre main sequence to He core depletion. The extent of the convective penetration zone is self-consistently calculated at each time step without introducing new free parameters. We find both a substantial penetration zone in all models with a convective core and observable differences to global stellar properties such as the luminosity and radius. We preset how the predicted radial extent of the penetration zone scales with the total stellar mass, age and the metallicity of the star. We discuss our results in the context of existing numerical and observational studies. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: Accepted for publication in The Astrophysical Journal, 16 pages, 8 figures, accompanying repository (https://github.com/mathren/dbconvpen)

MESA-Web: A cloud resource for stellar evolution in astronomy curricula

Authors: Carl E. Fields, Richard H. D. Townsend, A. L. Dotter, Michael Zingale, F. X. Timmes

Abstract: We present MESA-Web, a cloud resource with an online interface to the Modules for Experiments in Stellar Astrophysics (MESA) software instrument. MESA-Web allows learners to evolve stellar models without the need to download and install MESA. Since being released in 2015, MESA-Web has delivered over 17,000 calculations to over 2,200 unique learners and currently performs about 11 jobs per day. MES… ▽ More We present MESA-Web, a cloud resource with an online interface to the Modules for Experiments in Stellar Astrophysics (MESA) software instrument. MESA-Web allows learners to evolve stellar models without the need to download and install MESA. Since being released in 2015, MESA-Web has delivered over 17,000 calculations to over 2,200 unique learners and currently performs about 11 jobs per day. MESA-Web can be used as an educational tool for stars in the classroom or for scientific investigations. We report on new capabilities of MESA-Web introduced since its 2015 release including learner-supplied nuclear reaction rates, custom stopping conditions, and an expanded selection of input parameters. To foster collaboration we have created a Zenodo MESA-Web community hub where instructors can openly share examples of using MESA-Web in the classroom. We discuss two examples in the current community hub. The first example is a lesson module on Red Giant Branch stars that includes a suite of exercises designed to fit a range of learners and a Jupyter workbook for additional analysis. The second example is lesson materials for an upper-level Astronomy majors course in Stars and Radiation that includes an assignment verifying some of the expected trends that are presented in a popular stellar physics textbook. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: Accepted for publication in the Astronomy Education Journal. We encourage users of MESA-Web to share their educational materials with the larger community using our Zenodo hub at https://zenodo.org/communities/mesa-web/

Corrected Weight Functions for Stellar Oscillation Eigenfrequencies

Authors: R. H. D. Townsend, S. D. Kawaler

Abstract: Kawaler et al. (1985) present a variational expression for the eigenfrequencies associated with stellar oscillations. We highlight and correct a typographical error in the weight functions appearing in these expressions, and validate the correction numerically. Kawaler et al. (1985) present a variational expression for the eigenfrequencies associated with stellar oscillations. We highlight and correct a typographical error in the weight functions appearing in these expressions, and validate the correction numerically. △ Less

Submitted 3 August, 2023; v1 submitted 1 August, 2023; originally announced August 2023.

Comments: 3 pages, 1 figure

Journal ref: RNAAS, 7, 166 (023)

Overview and Validation of the Asteroseismic Modeling Portal v2.0

Authors: Travis S. Metcalfe, Richard H. D. Townsend, Warrick H. Ball

Abstract: The launch of NASA's Kepler space telescope in 2009 revolutionized the quality and quantity of observational data available for asteroseismic analysis. While Kepler was able to detect solar-like oscillations in hundreds of main-sequence and subgiant stars, the Transiting Exoplanet Survey Satellite (TESS) is now making similar observations for thousands of the brightest stars in the sky. The Astero… ▽ More The launch of NASA's Kepler space telescope in 2009 revolutionized the quality and quantity of observational data available for asteroseismic analysis. While Kepler was able to detect solar-like oscillations in hundreds of main-sequence and subgiant stars, the Transiting Exoplanet Survey Satellite (TESS) is now making similar observations for thousands of the brightest stars in the sky. The Asteroseismic Modeling Portal (AMP) is an automated and objective stellar model-fitting pipeline for asteroseismic data, which was originally developed to use models from the Aarhus Stellar Evolution Code (ASTEC). We briefly summarize an updated version of the AMP pipeline that uses Modules for Experiments in Stellar Astrophysics (MESA), and we present initial modeling results for the Sun and several solar analogs to validate the precision and accuracy of the inferred stellar properties. △ Less

Submitted 30 July, 2023; originally announced July 2023.

Comments: 3 pages, 1 table, AAS Journals accepted. Software available at https://github.com/travismetcalfe/amp2

Journal ref: Res. Notes AAS 7, 164 (2023)

A calibration point for stellar evolution from massive star asteroseismology

Authors: Siemen Burssens, Dominic M. Bowman, Mathias Michielsen, Sergio Simón-Díaz, Conny Aerts, Vincent Vanlaer, Gareth Banyard, Nicolas Nardetto, Richard H. D. Townsend, Gerald Handler, Joey S. G. Mombarg, Roland Vanderspek, George Ricker

Abstract: Massive stars are progenitors of supernovae, neutron stars and black holes. During the hydrogen-core burning phase their convective cores are the prime drivers of their evolution, but inferences of core masses are subject to unconstrained boundary mixing processes. Moreover, uncalibrated transport mechanisms can lead to strong envelope mixing and differential radial rotation. Ascertaining the effi… ▽ More Massive stars are progenitors of supernovae, neutron stars and black holes. During the hydrogen-core burning phase their convective cores are the prime drivers of their evolution, but inferences of core masses are subject to unconstrained boundary mixing processes. Moreover, uncalibrated transport mechanisms can lead to strong envelope mixing and differential radial rotation. Ascertaining the efficiency of the transport mechanisms is challenging because of a lack of observational constraints. Here we deduce the convective core mass and robustly demonstrate non-rigid radial rotation in a supernova progenitor, the $12.0^{+1.5}_{-1.5}$ solar-mass hydrogen-burning star HD 192575, using asteroseismology, TESS photometry, high-resolution spectroscopy, and Gaia astrometry. We infer a convective core mass ($M_{\rm cc} = 2.9^{+0.5}_{-0.8}$ solar masses), and find the core to be rotating between 1.4 and 6.3 times faster than the stellar envelope depending on the location of the rotational shear layer. Our results deliver a robust inferred core mass of a massive star using asteroseismology from space-based photometry. HD 192575 is a unique anchor point for studying interior rotation and mixing processes, and thus also angular momentum transport mechanisms inside massive stars. △ Less

Submitted 23 June, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: 41 pages, 5 figures, 1 table. Version comment: updated erroneous affiliation

The photometric variability of massive stars due to gravity waves excited by core convection

Authors: Evan H. Anders, Daniel Lecoanet, Matteo Cantiello, Keaton J. Burns, Benjamin A. Hyatt, Emma Kaufman, Richard H. D. Townsend, Benjamin P. Brown, Geoffrey M. Vasil, Jeffrey S. Oishi, Adam S. Jermyn

Abstract: Massive stars die in catastrophic explosions, which seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion's character and the remnant mass depend on models of the star's evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometri… ▽ More Massive stars die in catastrophic explosions, which seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion's character and the remnant mass depend on models of the star's evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometric observations reveal a ubiquitous red noise signal on massive main sequence stars; a hypothesized source of this noise is gravity waves driven by core convection. We present the first 3D simulations of massive star convection extending from the star's center to near its surface, with realistic stellar luminosities. Using these simulations, we make the first prediction of photometric variability due to convectively-driven gravity waves at the surfaces of massive stars, and find that gravity waves produce photometric variability of a lower amplitude and lower characteristic frequency than the observed red noise. We infer that the photometric signal of gravity waves excited by core convection is below the noise limit of current observations, so the red noise must be generated by an alternative process. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: As accepted for publication in Nature Astronomy except for final editorial revisions. Supplemental materials available online at https://doi.org/10.5281/zenodo.7764997 . We have also sonified our results to make them more accessible, see https://github.com/evanhanders/gmode_variability_paper/blob/main/sound/gmode_sonification.pdf

Journal ref: Nature Astronomy, 2023, 10.1038/s41550-023-02040-7

Discrepant Approaches to Modeling Stellar Tides, and the Blurring of Pseudosynchronization

Authors: R. H. D. Townsend, M. Sun

Abstract: We examine the reasons for discrepancies between two alternative approaches to modeling small-amplitude tides in binary systems. The 'direct solution' (DS) approach solves the governing differential equations and boundary conditions directly, while the 'modal decomposition' (MD) approach relies on a normal-mode expansion. Applied to a model for the primary star in the heartbeat system KOI-54, the… ▽ More We examine the reasons for discrepancies between two alternative approaches to modeling small-amplitude tides in binary systems. The 'direct solution' (DS) approach solves the governing differential equations and boundary conditions directly, while the 'modal decomposition' (MD) approach relies on a normal-mode expansion. Applied to a model for the primary star in the heartbeat system KOI-54, the two approaches predict quite different behavior of the secular tidal torque. The MD approach exhibits the pseudosynchronization phenomenon, where the torque due to the equilibrium tide changes sign at a single, well-defined and theoretically predicted stellar rotation rate. The DS approach instead shows 'blurred' pseudosynchronization, where positive and negative torques intermingle over a range of rotation rates. We trace a major source of these differences to an incorrect damping coefficient in the profile functions describing the frequency dependence of the MD expansion coefficients. With this error corrected some differences between the approaches remain; however, both are in agreement that pseudosynchronization is blurred in the KOI-54 system. Our findings generalize to any type of star for which the tidal damping depends explicitly or implicitly on the forcing frequency. △ Less

Submitted 10 June, 2023; originally announced June 2023.

Comments: 8 pages, 6 figures, accepted by ApJ

Strong Variability in AzV 493, an Extreme Oe-Type Star in the SMC

Authors: M. S. Oey, N. Castro, M. Renzo, I. Vargas-Salazar, M. W. Suffak, M. Ratajczak, J. D. Monnier, M. K. Szymanski, G. D. Phillips, N. Calvet, A. Chiti, G. Micheva, K. C. Rasmussen, R. H. D. Townsend

Abstract: We present 18 years of OGLE photometry together with spectra obtained over 12 years, revealing that the early Oe star AzV 493 shows strong photometric (Delta I < 1.2 mag) and spectroscopic variability with a dominant, 14.6-year pattern and ~40-day oscillations. We estimate stellar parameters T_eff = 42000 K, log L/L_sun = 5.83 +/- 0.15, M/M_sun = 50 +/- 9, and vsini = 370 +/- 40 km/s. Direct spect… ▽ More We present 18 years of OGLE photometry together with spectra obtained over 12 years, revealing that the early Oe star AzV 493 shows strong photometric (Delta I < 1.2 mag) and spectroscopic variability with a dominant, 14.6-year pattern and ~40-day oscillations. We estimate stellar parameters T_eff = 42000 K, log L/L_sun = 5.83 +/- 0.15, M/M_sun = 50 +/- 9, and vsini = 370 +/- 40 km/s. Direct spectroscopic evidence shows episodes of both gas ejection and infall. There is no X-ray detection, and it is likely a runaway star. AzV 493 may have an unseen companion on a highly eccentric (e > 0.93) orbit. We propose that close interaction at periastron excites ejection of the decretion disk, whose variable emission-line spectrum suggests separate inner and outer components, with an optically thick outer component obscuring both the stellar photosphere and the emission-line spectrum of the inner disk at early phases in the photometric cycle. It is plausible that AzV 493's mass and rotation have been enhanced by binary interaction followed by the core-collapse supernova explosion of the companion, which now could be either a black hole or neutron star. This system in the Small Magellanic Cloud can potentially shed light on OBe decretion disk formation and evolution, massive binary evolution, and compact binary progenitors. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Comments: 24 pages, 13 figures, 1 table. Accepted by ApJ

gyre_tides: Modeling binary tides within the gyre stellar oscillation code

Authors: Meng Sun, R. H. D. Townsend, Zhao Guo

Abstract: We describe new functionality in the GYRE stellar oscillation code for modeling tides in binary systems. Using a multipolar expansion in space and a Fourier-series expansion in time, we decompose the tidal potential into a superposition of partial tidal potentials. The equations governing the small-amplitude response of a spherical star to an individual partial potential are the linear, non-radial… ▽ More We describe new functionality in the GYRE stellar oscillation code for modeling tides in binary systems. Using a multipolar expansion in space and a Fourier-series expansion in time, we decompose the tidal potential into a superposition of partial tidal potentials. The equations governing the small-amplitude response of a spherical star to an individual partial potential are the linear, non-radial, non-adiabatic oscillation equations with an extra inhomogeneous forcing term. We introduce a new executable, gyre_tides, that directly solves these equations within the GYRE numerical framework. Applying this to selected problems, we find general agreement with results in the published literature but also uncover some differences between our direct solution methodology and the modal decomposition approach adopted by many authors. In its present form gyre_tides can model equilibrium and dynamical tides of aligned binaries in which radiative diffusion dominates the tidal dissipation (typically, intermediate and high-mass stars on the main sequence). Milestones for future development include incorporation of other dissipation processes, spin-orbit misalignment, and the Coriolis force arising from rotation. △ Less

Submitted 18 January, 2023; v1 submitted 16 January, 2023; originally announced January 2023.

Comments: 15 pages, 8 figures, accepted for publication in ApJ; This new functionality, along with GYRE version 7.0, will be release in a few days from today (Jan 16th, 2023)

Modules for Experiments in Stellar Astrophysics (MESA): Time-Dependent Convection, Energy Conservation, Automatic Differentiation, and Infrastructure

Authors: Adam S. Jermyn, Evan B. Bauer, Josiah Schwab, R. Farmer, Warrick H. Ball, Earl P. Bellinger, Aaron Dotter, Meridith Joyce, Pablo Marchant, Joey S. G. Mombarg, William M. Wolf, Tin Long Sunny Wong, Giulia C. Cinquegrana, Eoin Farrell, R. Smolec, Anne Thoul, Matteo Cantiello, Falk Herwig, Odette Toloza, Lars Bildsten, Richard H. D. Townsend, F. X. Timmes

Abstract: We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite difference approximations. We significantly enhance the treatment of the growth and decay of convection in MES… ▽ More We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite difference approximations. We significantly enhance the treatment of the growth and decay of convection in MESA with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron degenerate ignition events. We strengthen MESA's implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in MESA we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator split nuclear burning mode. We close by discussing major updates to MESA's software infrastructure that enhance source code development and community engagement. △ Less

Submitted 30 December, 2022; v1 submitted 7 August, 2022; originally announced August 2022.

Comments: 50 pages, 29 figures; Accepted to ApJS

A New Window to Tidal Asteroseismology: Non-linearly Excited Stellar Eigenmodes and the Period Spacing Pattern in KOI-54

Authors: Zhao Guo, Gordon I. Ogilvie, Gang Li, Richard H. D. Townsend, Meng Sun

Abstract: We revisit the Tidally Excited Oscillations (TEOs) in the A-type, main-sequence, eccentric binary KOI-54, the prototype of heartbeat stars. Although the linear tidal response of the star is a series of orbital-harmonic frequencies which are not stellar eigenfrequencies, we show that the non-linearly excited non-orbital-harmonic TEOs are eigenmodes. By carefully choosing the modes which satisfy the… ▽ More We revisit the Tidally Excited Oscillations (TEOs) in the A-type, main-sequence, eccentric binary KOI-54, the prototype of heartbeat stars. Although the linear tidal response of the star is a series of orbital-harmonic frequencies which are not stellar eigenfrequencies, we show that the non-linearly excited non-orbital-harmonic TEOs are eigenmodes. By carefully choosing the modes which satisfy the mode-coupling selection rules, a period spacing ($ΔP$) pattern of quadrupole gravity modes ($ΔP \approx 2520-2535$ sec) can be discerned in the Fourier spectrum, with a detection significance level of $99.9\%$. The inferred period spacing value agrees remarkably well with the theoretical $l=2,m=0$ g modes from a stellar model with the measured mass, radius and effective temperature. We also find that the two largest-amplitude TEOs at $N=90, 91$ harmonics are very close to resonance with $l=2,m=0$ eigenmodes, and likely come from different stars. Previous works on tidal oscillations primarily focus on the modeling of TEO amplitudes and phases, the high sensitivity of TEO amplitude to the frequency detuning (tidal forcing frequency minus the closest stellar eigenfrequency) requires extremely dense grids of stellar models and prevents us from constraining the stellar physical parameters easily. This work, however, opens the window of real tidal asteroseismology by using the eigenfrequencies of the star inferred from the non-linear TEOs and possibly very-close-to-resonance linear TEOs. Our seismic modeling of these identified eigen g-modes shows that the best-matching stellar models have ($M \approx 2.20, 2.35 M_{\odot}$) and super-solar metallicity, in good agreement with previous measurements. △ Less

Submitted 11 September, 2022; v1 submitted 12 February, 2022; originally announced February 2022.

Comments: MNRAS Accepted, comments welcome

On The Impact Of 22Ne On The Pulsation Periods Of Carbon-Oxygen White Dwarfs With Helium Dominated Atmospheres

Authors: Morgan T. Chidester, F. X. Timmes, Josiah Schwab, Richard H. D. Townsend, Ebraheem Farag, Anne Thoul, C. E. Fields, Evan B. Bauer, Michael H. Montgomery

Abstract: We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M$_\odot$ carbon-oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of $^{22}$Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6$-$7 significant figures of precision. Usually white dwarf models wi… ▽ More We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M$_\odot$ carbon-oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of $^{22}$Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6$-$7 significant figures of precision. Usually white dwarf models without $^{22}$Ne are fit to the observed periods and other properties. The root-mean-square residuals to the $\simeq$ 150$-$400 s low-order g-mode periods are typically in the range of $σ_{rms}$ $\lesssim$ 0.3 s, for a fit precision of $σ_{rms}/ P$ $\lesssim$ 0.3 %. We find average relative period shifts of $ΔP/P$ $\simeq$ $\pm$ 0.5 % for the low-order dipole and quadrupole g-mode pulsations within the observed effective temperature window, with the range of $ΔP/P$ depending on the specific g-mode, abundance of $^{22}$Ne, effective temperature, and mass of the white dwarf model. This finding suggests a systematic offset may be present in the fitting process of specific white dwarfs when $^{22}$Ne is absent. As part of the fitting processes involves adjusting the composition profiles of a white dwarf model, our study on the impact of $^{22}$Ne can provide new inferences on the derived interior mass fraction profiles. We encourage routinely including $^{22}$Ne mass fraction profiles, informed by stellar evolution models, to future generations of white dwarf model fitting processes. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 23 pages, 19 figures total. Accepted to the ApJ. Presented results at 2021 AAS January meeting

TESS Asteroseismology of $α$ Mensae: Benchmark Ages for a G7 Dwarf and its M-dwarf Companion

Authors: Ashley Chontos, Daniel Huber, Travis A. Berger, Hans Kjeldsen, Aldo M. Serenelli, Victor Silva Aguirre, Warrick H. Ball, Sarbani Basu, Timothy R. Bedding, William J. Chaplin, Zachary R. Claytor, Enrico Corsaro, Rafael A. García, Steve B. Howell, Mia S. Lundkvist, Savita Mathur, Travis S. Metcalfe, Martin B. Nielsen, Jia Mian Joel Ong, Zeynep Çelik Orhan, Sibel Örtel, Maïssa Salama, Keivan G. Stassun, R. H. D. Townsend, Jennifer L. van Saders , et al. (5 additional authors not shown)

Abstract: Asteroseismology of bright stars has become increasingly important as a method to determine fundamental properties (in particular ages) of stars. The Kepler Space Telescope initiated a revolution by detecting oscillations in more than 500 main-sequence and subgiant stars. However, most Kepler stars are faint, and therefore have limited constraints from independent methods such as long-baseline int… ▽ More Asteroseismology of bright stars has become increasingly important as a method to determine fundamental properties (in particular ages) of stars. The Kepler Space Telescope initiated a revolution by detecting oscillations in more than 500 main-sequence and subgiant stars. However, most Kepler stars are faint, and therefore have limited constraints from independent methods such as long-baseline interferometry. Here, we present the discovery of solar-like oscillations in $α$ Men A, a naked-eye (V=5.1) G7 dwarf in TESS's Southern Continuous Viewing Zone. Using a combination of astrometry, spectroscopy, and asteroseismology, we precisely characterize the solar analog alpha Men A (Teff = 5569 +/- 62 K, R = 0.960 +/- 0.016 Rsun, M = 0.964 +/- 0.045 Msun). To characterize the fully convective M dwarf companion, we derive empirical relations to estimate mass, radius, and temperature given the absolute Gaia magnitude and metallicity, yielding M = 0.169 +/- 0.006, R = 0.19 +/- 0.01 and Teff = 3054 +/- 44 K. Our asteroseismic age of 6.2 +/- 1.4 (stat) +/- 0.6 (sys) Gyr for the primary places $α$ Men B within a small population of M dwarfs with precisely measured ages. We combined multiple ground-based spectroscopy surveys to reveal an activity cycle of 13.1 +/- 1.1 years, a period similar to that observed in the Sun. We used different gyrochronology models with the asteroseismic age to estimate a rotation period of ~30 days for the primary. Alpha Men A is now the closest (d=10pc) solar analog with a precise asteroseismic age from space-based photometry, making it a prime target for next-generation direct imaging missions searching for true Earth analogs. △ Less

Submitted 4 December, 2021; v1 submitted 19 December, 2020; originally announced December 2020.

Comments: Accepted to The Astrophysical Journal; 15 pages, 10 figures

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

Authors: Meng Sun, Robert D. Mathieu, Emily M. Leiner, R. H. D. Townsend

Abstract: 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 t… ▽ More 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}$. △ Less

Submitted 15 December, 2020; originally announced December 2020.

Comments: 17 pages, 16 figures, accepted for publication in ApJ

Stability and pulsation of the first dark stars

Authors: Tanja Rindler-Daller, Katherine Freese, Richard H. D. Townsend, Luca Visinelli

Abstract: The first bright objects to form in the Universe might not have been "ordinary" fusion-powered stars, but "Dark Stars" (DSs) powered by the annihilation of dark matter (DM) in the form of Weakly Interacting Massive Particles (WIMPs). If discovered, DSs can provide a unique laboratory to test DM models. DSs are born with a mass of order $M_\odot$ and may grow to a few million solar masses; in this… ▽ More The first bright objects to form in the Universe might not have been "ordinary" fusion-powered stars, but "Dark Stars" (DSs) powered by the annihilation of dark matter (DM) in the form of Weakly Interacting Massive Particles (WIMPs). If discovered, DSs can provide a unique laboratory to test DM models. DSs are born with a mass of order $M_\odot$ and may grow to a few million solar masses; in this work we investigate the properties of early DSs with masses up to $\sim \! 1000 \, M_\odot$, fueled by WIMPS weighing $100$ GeV. We improve the previous implementation of the DM energy source into the stellar evolution code MESA. We show that the growth of DSs is not limited by astrophysical effects: DSs up to $\sim \! 1000 \, M_\odot$ exhibit no dynamical instabilities; DSs are not subject to mass-loss driven by super-Eddington winds. We test the assumption of previous work that the injected energy per WIMP annihilation is constant throughout the star; relaxing this assumption does not change the properties of the DSs. Furthermore, we study DS pulsations, for the first time investigating non-adiabatic pulsation modes, using the linear pulsation code GYRE. We find that acoustic modes in DSs of masses smaller than $\sim \! 200 \, M_\odot$ are excited by the $κ-γ$ and $γ$ mechanism in layers where hydrogen or helium is (partially) ionized. Moreover, we show that the mass-loss rates potentially induced by pulsations are negligible compared to the accretion rates. △ Less

Submitted 12 February, 2021; v1 submitted 31 October, 2020; originally announced November 2020.

Comments: minor revisions; accepted by MNRAS; 15 pages, 6 Figures

Report number: NORDITA-2020-078; LCTP-20-21; UTTG-20-2020

Journal ref: Mon. Not. Roy. Astron. Soc. 503 3, 3677-3691 (2021)

How the breakout-limited mass in B-star centrifugal magnetospheres controls their circumstellar H-alpha emission

Authors: Stanley P. Owocki, Matt E. Shultz, Asif ud-Doula, Jon O. Sundqvist, Richard H. D. Townsend, Steven R. Cranmer

Abstract: Strongly magnetic B-type stars with moderately rapid rotation form `centrifugal magnetospheres' (CMs), from the magnetic trapping of stellar wind material in a region above the Kepler co-rotation radius. A longstanding question is whether the eventual loss of such trapped material occurs from gradual drift and/or diffusive leakage, or through sporadic `{\em centrifugal break out}' (CBO) events, wh… ▽ More Strongly magnetic B-type stars with moderately rapid rotation form `centrifugal magnetospheres' (CMs), from the magnetic trapping of stellar wind material in a region above the Kepler co-rotation radius. A longstanding question is whether the eventual loss of such trapped material occurs from gradual drift and/or diffusive leakage, or through sporadic `{\em centrifugal break out}' (CBO) events, wherein magnetic tension can no longer contain the built-up mass. We argue here that recent empirical results for Balmer-$α$ emission from such B-star CMs strongly favor the CBO mechanism. Most notably, the fact that the onset of such emission depends mainly on the field strength at the Kepler radius, and is largely {\em independent} of the stellar luminosity, strongly disfavors any drift/diffusion process, for which the net mass balance would depend on the luminosity-dependent wind feeding rate. In contrast, we show that in a CBO model the {\em maximum confined mass} in the magnetosphere is independent of this wind feeding rate, and has a dependence on field strength and Kepler radius that naturally explains the empirical scalings for the onset of H$α$ emission, its associated equivalent width, and even its line profile shapes. However, the general lack of observed Balmer emission in late-B and A-type stars could still be attributed to a residual level of diffusive or drift leakage that does not allow their much weaker winds to fill their CMs to the breakout level needed for such emission; alternatively this might result from a transition to a metal-ion wind that lacks the requisite Hydrogen. △ Less

Submitted 25 September, 2020; originally announced September 2020.

Comments: 13 pages, 11 figures, accepted for MNRAS

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

Authors: J. Goldstein, R. H. D. Townsend

Abstract: 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 f… ▽ More 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. △ Less

Submitted 21 July, 2020; v1 submitted 23 June, 2020; originally announced June 2020.

Comments: 12 pages, 14 figures, Accepted to ApJ

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

Authors: R. H. D. Townsend

Abstract: 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… ▽ More 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. △ Less

Submitted 21 July, 2020; v1 submitted 22 June, 2020; originally announced June 2020.

Comments: 10 pages, 2 figures, accepted by MNRAS

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

Authors: J. Christensen-Dalsgaard, V. Silva Aguirre, S. Cassisi, M. Miller Bertolami, A. Serenelli, D. Stello, A. Weiss, G. Angelou, C. Jiang, Y. Lebreton, F. Spada, E. P. Bellinger, S. Deheuvels, R. M. Ouazzani, A. Pietrinferni, J. R. Mosumgaard, R. H. D. Townsend, T. Battich, D. Bossini, T. Constantino, P. Eggenberger, S. Hekker, A. Mazumdar, A. Miglio, K. B. Nielsen , et al. (1 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 10 February, 2020; v1 submitted 7 February, 2020; originally announced February 2020.

Comments: Astronomy and Astrophysics, in the press. Paper resubmitted to arXiv to correct author list in metadata

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

Authors: Z. Keszthelyi, G. Meynet, M. E. Shultz, A. David-Uraz, A. ud-Doula, R. H. D. Townsend, G. A. Wade, C. Georgy, V. Petit, S. P. Owocki

Abstract: 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 magneti… ▽ More 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. △ Less

Submitted 21 January, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: accepted for publication in MNRAS, a zenodo record is available at: https://zenodo.org/record/3250412

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

Authors: V. Silva Aguirre, J. Christensen-Dalsgaard, S. Cassisi, M. Miller Bertolami, A. M. Serenelli, D. Stello, A. Weiss, G. Angelou, C. Jiang, Y. Lebreton, F. Spada, E. P. Bellinger, S. Deheuvels, R. M. Ouazzani, A. Pietrinferni, J. R. Mosumgaard, R. H. D. Townsend, T. Battich, D. Bossini, T. Constantino, P. Eggenberger, S. Hekker, A. Mazumdar, A. Miglio, K. B. Nielsen , et al. (1 additional authors not shown)

Abstract: (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 c… ▽ More (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. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 635, A164 (2020)

Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation

Authors: Bill Paxton, R. Smolec, Josiah Schwab, A. Gautschy, Lars Bildsten, Matteo Cantiello, Aaron Dotter, R. Farmer, Jared A. Goldberg, Adam S. Jermyn, S. M. Kanbur, Pablo Marchant, Anne Thoul, Richard H. D. Townsend, William M. Wolf, Michael Zhang, F. X. Timmes

Abstract: We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For exam… ▽ More We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development. △ Less

Submitted 16 May, 2019; v1 submitted 4 March, 2019; originally announced March 2019.

Comments: 57 pages, 57 figures; Accepted to ApJS

The Impact of White Dwarf Luminosity Profiles on Oscillation Frequencies

Authors: F. X. Timmes, Richard H. D. Townsend, Evan B. Bauer, Anne Thoul, C. E. Fields, William M. Wolf

Abstract: KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportion… ▽ More KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportional to the enclosed mass, $L_r \propto M_r$, independent of the effective temperature $T_{\rm eff}$. Evolutionary models of young white dwarfs with $T_{\rm eff} \gtrsim$ 25,000 K suggest neutrino emission gives rise to luminosity profiles with $L_r$ $\not\propto$ $M_r$. We explore this contrast by comparing the oscillation frequencies between two nearly identical white dwarf models: one with an enforced $L_r \propto M_r$ luminosity profile and the other with a luminosity profile determined by the star's previous evolution history. We find the low order g-mode frequencies differ by up to $\simeq$ 70 $μ$Hz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an $L_r \propto M_r$ profile may have uncorrected, effectively-random errors at the level of tens of $μ$Hz. A mean frequency difference of 30 $μ$Hz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of $\simeq$ 12% in white dwarf mass, $\simeq$ 9% in the radius, and $\simeq$ 3% in the central oxygen mass fraction. △ Less

Submitted 8 October, 2018; originally announced October 2018.

Comments: 7 pages, 5 figures; Accepted for publication in ApJ Letters

The Tayler Instability in the Anelastic Approximation

Authors: J. Goldstein, R. H. D. Townsend, E. G. Zweibel

Abstract: The Tayler instability (TI) is a non-axisymmetric linear instability of an axisymmetric toroidal magnetic field in magneto-hydrostatic equilibrium (MHSE). Spruit (1999, 2002) has proposed that in a differentially rotating radiative region of a star, the TI drives a dynamo which generates magnetic fields that can efficiently transport angular momentum; a parameterized version of this dynamo has bee… ▽ More The Tayler instability (TI) is a non-axisymmetric linear instability of an axisymmetric toroidal magnetic field in magneto-hydrostatic equilibrium (MHSE). Spruit (1999, 2002) has proposed that in a differentially rotating radiative region of a star, the TI drives a dynamo which generates magnetic fields that can efficiently transport angular momentum; a parameterized version of this dynamo has been implemented in stellar structure and evolution codes and shown to be important for determining interior spin. Numerical simulations, however, have yet to definitively demonstrate the operation of the dynamo. A criterion for the MHSE to develop the TI was derived using fully-compressible magneto-hydrodynamics, while numerical simulations of dynamical processes in stars frequently use an anelastic approximation. This motivates us to derive a new anelastic Tayler instability (anTI) criterion. We find that some MHSE configurations unstable in the fully-compressible case, become stable in the anelastic case. We find and characterize the unstable modes of a simple family of cylindrical MHSE configurations using numerical calculations, and discuss the implications for fully non-linear anelastic simulations. △ Less

Submitted 18 July, 2019; v1 submitted 27 August, 2018; originally announced August 2018.

Comments: 11 pages, 3 figures, 1 appendix

On the sensitivity of gravito-inertial modes to differential rotation in intermediate-mass main-sequence stars

Authors: T. Van Reeth, J. S. G. Mombarg, S. Mathis, A. Tkachenko, J. Fuller, D. M. Bowman, B. Buysschaert, C. Johnston, A. García Hernández, J. Goldstein, R. H. D. Townsend, C. Aerts

Abstract: Context. While rotation has a major impact on stellar structure and evolution, its effects are not well understood. Thanks to high- quality and long timebase photometric observations obtained with recent space missions, we are now able to study stellar rotation more precisely. Aims. We aim to constrain radial differential rotation profiles in gamma Doradus (gamma Dor) stars, and to develop new the… ▽ More Context. While rotation has a major impact on stellar structure and evolution, its effects are not well understood. Thanks to high- quality and long timebase photometric observations obtained with recent space missions, we are now able to study stellar rotation more precisely. Aims. We aim to constrain radial differential rotation profiles in gamma Doradus (gamma Dor) stars, and to develop new theoretical seismic diagnosis for such stars with rapid and potentially non-uniform rotation. Methods. We derive a new asymptotic description which accounts for the impact of weak differential near-core rotation on gravity- mode period spacings. The theoretical predictions are illustrated from pulsation computations with the code GYRE and compared with observations of gamma Dor stars. When possible, we also derive the surface rotation rates in these stars by detecting and analysing signatures of rotational modulation, and compute the core-to-surface rotation ratios. Results. Stellar rotation has to be strongly differential before its effects on period spacing patterns can be detected, unless multiple period spacing patterns can be compared. Six stars in our sample exhibit a single unexplained period spacing pattern of retrograde modes. We hypothesise that these are Yanai modes. Finally, we find signatures of rotational spot modulation in the photometric data of eight targets. Conclusions. If only one period spacing pattern is detected and analysed for a star, it is difficult to detect differential rotation. A rigidly rotating model will often provide the best solution. Differential rotation can only be detected when multiple period spacing patterns have been found for a single star or its surface rotation rate is known as well. This is the case for eight stars in our sample, revealing surface-to-core rotation ratios between 0.95 and 1.05. △ Less

Submitted 26 June, 2018; v1 submitted 10 June, 2018; originally announced June 2018.

Comments: 15 pages, 13 figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 618, A24 (2018)

Non-Radial Pulsations in Post-Outburst Novae

Authors: William M. Wolf, Richard H. D. Townsend, Lars Bildsten

Abstract: After an optical peak, a classical or recurrent nova settles into a brief (days to years) period of quasi-stable thermonuclear burning in a compact configuration nearly at the white dwarf (WD) radius. During this time, the underlying WD becomes visible as a strong emitter of supersoft X-rays. Observations during this phase have revealed oscillations in the X-ray emission with periods on the order… ▽ More After an optical peak, a classical or recurrent nova settles into a brief (days to years) period of quasi-stable thermonuclear burning in a compact configuration nearly at the white dwarf (WD) radius. During this time, the underlying WD becomes visible as a strong emitter of supersoft X-rays. Observations during this phase have revealed oscillations in the X-ray emission with periods on the order of tens of seconds. A proposed explanation for the source of these oscillations are internal gravity waves excited by nuclear reactions at the base of the hydrogen-burning layer. In this work, we present the first models exhibiting unstable surface $g$-modes with periods similar to oscillation periods found in galactic novae. However, when comparing mode periods of our models to the observed oscillations of several novae, we find that the modes which are excited have periods shorter than that observed. △ Less

Submitted 5 February, 2018; originally announced February 2018.

Comments: 13 pages, 7 figures. Accepted to ApJ

Angular momentum transport by heat-driven g-modes in slowly pulsating B stars

Authors: R. H. D. Townsend, J. Goldstein, E. G. Zweibel

Abstract: Motivated by recent interest in the phenomenon of waves transport in massive stars, we examine whether the heat-driven gravity (g) modes excited in slowly-pulsating B (SPB) stars can significantly modify the stars' internal rotation. We develop a formalism for the differential torque exerted by g modes, and implement this formalism using the GYRE oscillation code and the MESASTAR stellar evolution… ▽ More Motivated by recent interest in the phenomenon of waves transport in massive stars, we examine whether the heat-driven gravity (g) modes excited in slowly-pulsating B (SPB) stars can significantly modify the stars' internal rotation. We develop a formalism for the differential torque exerted by g modes, and implement this formalism using the GYRE oscillation code and the MESASTAR stellar evolution code. Focusing first on a $4.21$ $M_\odot$ model, we simulate 1,000 years of stellar evolution under the combined effects of the torque due to a single unstable prograde g mode (with an amplitude chosen on the basis of observational constraints), and diffusive angular momentum transport due to convection, overshooting, and rotational instabilities. We find that the g mode rapidly extracts angular momentum from the surface layers, depositing it deeper in the stellar interior. The angular momentum transport is so efficient that by the end of the simulation the initially non-rotating surface layers are spun in the retrograde direction to $\approx30\%$ of the critical rate. However, the additional inclusion of magnetic stresses in our simulations, almost completely inhibits this spin-up. Expanding our simulations to cover the whole instability strip, we show that the same general behavior is seen in all SPB stars. After providing some caveats to contextualize our results, we hypothesize that the observed slower surface rotation of SPB stars (as compared to other B-type stars) may be the direct consequence of the angular momentum transport that our simulations demonstrate. △ Less

Submitted 6 December, 2017; originally announced December 2017.

Comments: Accepted for publication in MNRAS (13 pages, 8 figures)

Stellar parameters of Be stars observed with X-shooter

Authors: A. Shokry, Th. Rivinius, A. Mehner, C. Martayan, W. Hummel, R. H. D. Townsend, A. Mérand, B. Mota, D. M. Faes, M. A. Hamdy, M. M. Beheary, K. A. K Gadallah, M. S. Abo-Elazm

Abstract: Aims. The X-shooter archive of several thousand telluric star spectra was skimmed for Be and Be-shell stars to derive the stellar fundamental parameters and statistical properties, in particular for the less investigated late type Be stars, and the extension of the Be phenomenon into early A stars. Methods. An adapted version of the BCD method is used, utilizing the Balmer discontinuity parameters… ▽ More Aims. The X-shooter archive of several thousand telluric star spectra was skimmed for Be and Be-shell stars to derive the stellar fundamental parameters and statistical properties, in particular for the less investigated late type Be stars, and the extension of the Be phenomenon into early A stars. Methods. An adapted version of the BCD method is used, utilizing the Balmer discontinuity parameters to determine effective temperature and surface gravity. This method is optimally suited for late B stars. The projected rotational velocity was obtained by profile fitting to the Mg ii lines of the targets, and the spectra were inspected visually for the presence of peculiar features such as the infrared Ca ii triplet or the presence of a double Balmer discontinuity. The Balmer line equivalent widths were measured, but due to uncertainties in determining the photospheric contribution are useful only in a subsample of Be stars for determining the pure emission contribution. Results. A total of 78 Be stars, mostly late type ones, were identified in the X-shooter telluric standard star archive, out of which 48 had not been reported before. The general trend of late type Be stars having more tenuous disks and being less variable than early type ones is confirmed. The relatively large number (48) of relatively bright (V > 8.5) additional Be stars casts some doubt on the statistics of late type Be stars; they are more common than currently thought: The Be/B star fraction may not strongly depend on spectral subtype. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: Accepted for publication in A&A

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

Modules for Experiments in Stellar Astrophysics (MESA): Convective Boundaries, Element Diffusion, and Massive Star Explosions

Authors: Bill Paxton, Josiah Schwab, Evan B. Bauer, Lars Bildsten, Sergei Blinnikov, Paul Duffell, R. Farmer, Jared A. Goldberg, Pablo Marchant, Elena Sorokina, Anne Thoul, Richard H. D. Townsend, F. X. Timmes

Abstract: We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white… ▽ More We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact. △ Less

Submitted 3 January, 2018; v1 submitted 23 October, 2017; originally announced October 2017.

Comments: 64 pages, 61 figures; Accepted to AAS Journals

Super-Eddington stellar winds: unifying radiative-enthalpy vs. flux-driven models

Authors: Stanley P. Owocki, Richard H. D. Townsend, Eliot Quataert

Abstract: We derive semi-analytic solutions for optically thick, super-Eddington stellar winds, induced by an assumed steady energy addition $Δ{\dot E}$ concentrated around a near-surface heating radius $R$ in a massive star of central luminosity $L_\ast$. We show that obtaining steady wind solutions requires both that the resulting total luminosity $L_o = L_\ast + Δ{\dot E}$ exceed the Eddington luminosity… ▽ More We derive semi-analytic solutions for optically thick, super-Eddington stellar winds, induced by an assumed steady energy addition $Δ{\dot E}$ concentrated around a near-surface heating radius $R$ in a massive star of central luminosity $L_\ast$. We show that obtaining steady wind solutions requires both that the resulting total luminosity $L_o = L_\ast + Δ{\dot E}$ exceed the Eddington luminosity, $Γ_o \equiv L_o/L_{Edd} > 1$, and that the induced mass loss rate be such that the "photon-tiring" parameter $m \equiv {\dot M} GM/R L_o \le 1-1/Γ_o$, ensuring the luminosity is sufficient to overcome the gravitational potential $GM/R$. Our analysis unifies previous super-Eddington wind models that either: (1) assumed a direct radiative flux-driving without accounting for the advection of radiative enthalpy that can become important in such an optically thick flow; or (2) assumed that such super-Eddington outflows are adiabatic, neglecting the effects of the diffusive radiative flux. We show that these distinct models become applicable in the asymptotic limits of small vs. large values of $m Γ_o $, respectively. By solving the coupled differential equations for radiative diffusion and wind momentum, we obtain general solutions that effectively bridge the behaviours of these limiting models. Two key scaling results are for the terminal wind speed to escape speed, which is found to vary as $v_\infty^2/v_{esc}^2 = Γ_o/(1+m Γ_o) -1$, and for the final observed luminosity $L_{ obs}$, which for all allowed steady-solutions with $m < 1 - 1/Γ_o$, exceeds the Eddington luminosity, $L_{obs} > L_{Edd}$. Our super-Eddington wind solutions have potential applicability for modeling phases of eruptive mass loss from massive stars, classical novae, and the remnants of stellar mergers. △ Less

Submitted 25 August, 2017; originally announced August 2017.

Comments: 13 pages, 11 figures, accepted for MNRAS

Magnetic massive stars as progenitors of "heavy" stellar-mass black holes

Authors: V. Petit, Z. Keszthelyi, R. MacInnis, D. H. Cohen, R. H. D. Townsend, G. A. Wade, S. L. Thomas, S. P. Owocki, J. Puls, J. A. ud-Doula

Abstract: The groundbreaking detection of gravitational waves produced by the inspiralling and coalescence of the black hole (BH) binary GW150914 confirms the existence of "heavy" stellar-mass BHs with masses >25 Msun. Initial modelling of the system by Abbott et al. (2016a) supposes that the formation of black holes with such large masses from the evolution of single massive stars is only feasible if the w… ▽ More The groundbreaking detection of gravitational waves produced by the inspiralling and coalescence of the black hole (BH) binary GW150914 confirms the existence of "heavy" stellar-mass BHs with masses >25 Msun. Initial modelling of the system by Abbott et al. (2016a) supposes that the formation of black holes with such large masses from the evolution of single massive stars is only feasible if the wind mass-loss rates of the progenitors were greatly reduced relative to the mass-loss rates of massive stars in the Galaxy, concluding that heavy BHs must form in low-metallicity (Z < 0.25-0.5 Zsun) environments. However, strong surface magnetic fields also provide a powerful mechanism for modifying mass loss and rotation of massive stars, independent of environmental metallicity (ud-Doula & Owocki 2002; ud-Doula et al. 2008). In this paper we explore the hypothesis that some heavy BHs, with masses >25 Msun such as those inferred to compose GW150914, could be the natural end-point of evolution of magnetic massive stars in a solar-metallicity environment. Using the MESA code, we developed a new grid of single, non-rotating, solar metallicity evolutionary models for initial ZAMS masses from 40-80 Msun that include, for the first time, the quenching of the mass loss due to a realistic dipolar surface magnetic field. The new models predict TAMS masses that are significantly greater than those from equivalent non-magnetic models, reducing the total mass lost by a strongly magnetized 80 Msun star during its main sequence evolution by 20 Msun. This corresponds approximately to the mass loss reduction expected from an environment with metallicity Z = 1/30 Zsun. △ Less

Submitted 27 November, 2016; originally announced November 2016.

Comments: 10 pages, 7 figures, 1 table. Accepted for publication in MNRAS

A JVLA survey of the high frequency radio emission of the massive magnetic B- and O-type stars

Authors: Sushma Kurapati, Poonam Chandra, Gregg Wade, David H. Cohen, Alexandre David-Uraz, Marc Gagne, Jason Grunhut, Mary E. Oksala, Veronique Petit, Matt Shultz, Jon Sundqvist, Richard H. D. Townsend, Asif ud-Doula

Abstract: We conducted a survey of seven magnetic O and eleven B-type stars with masses above $8M_{\odot}$ using the Very Large Array in the 1cm, 3cm and 13cm bands. The survey resulted in a detection of two O and two B-type stars. While the detected O-type stars - HD 37742 and HD 47129 - are in binary systems, the detected B-type stars, HD 156424 and ALS 9522, are not known to be in binaries. All four star… ▽ More We conducted a survey of seven magnetic O and eleven B-type stars with masses above $8M_{\odot}$ using the Very Large Array in the 1cm, 3cm and 13cm bands. The survey resulted in a detection of two O and two B-type stars. While the detected O-type stars - HD 37742 and HD 47129 - are in binary systems, the detected B-type stars, HD 156424 and ALS 9522, are not known to be in binaries. All four stars were detected at 3cm, whereas three were detected at 1cm and only one star was detected at 13cm. The detected B-type stars are significantly more radio luminous than the non-detected ones, which is not the case for O-type stars. The non-detections at 13cm are interpreted as due to thermal free-free absorption. Mass-loss rates were estimated using 3cm flux densities and were compared with theoretical mass-loss rates, which assume free-free emission. For HD 37742, the two values of the mass-loss rates were in good agreement, possibly suggesting that the radio emission for this star is mainly thermal. For the other three stars, the estimated mass-loss rates from radio observations were much higher than those expected from theory, suggesting either a possible contribution from non- thermal emission from the magnetic star or thermal or non-thermal emission due to interacting winds of the binary system, especially for HD 47129. All the detected stars are predicted to host centrifugal magnetospheres except HD 37742, which is likely to host a dynamical magnetosphere. This suggests that non-thermal radio emission is favoured in stars with centrifugal magnetospheres. △ Less

Submitted 1 November, 2016; originally announced November 2016.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Effect of a magnetic field on massive star winds I: mass-loss and velocity for a dipole field

Authors: Christopher Bard, Richard H. D. Townsend

Abstract: We generalize the Rigid-Field Hydrodynamic equations to accommodate arbitrary magnetic field topologies, resulting in a new Arbitrary Rigid-Field hydrodynamic (ARFHD) formalism. We undertake a critical point calculation of the steady-state ARFHD equations with a CAK-type radiative acceleration and determine the effects of a dipole magnetic field on the usual CAK mass-loss rate and velocity structu… ▽ More We generalize the Rigid-Field Hydrodynamic equations to accommodate arbitrary magnetic field topologies, resulting in a new Arbitrary Rigid-Field hydrodynamic (ARFHD) formalism. We undertake a critical point calculation of the steady-state ARFHD equations with a CAK-type radiative acceleration and determine the effects of a dipole magnetic field on the usual CAK mass-loss rate and velocity structure. Enforcing the proper optically-thin limit for the radiative line-acceleration is found to decrease both the mass-loss and wind acceleration, while rotation boosts both properties. We define optically-thin-correction and rotation parameters to quantify these effects on the global mass-loss rate and develop scaling laws for the surface mass-flux as a function of surface colatitude. These scaling laws are found to agree with previous laws derived from magnetohydrodynamic simulations of magnetospheres. The dipole magnetosphere velocity structure is found to differ from a global beta-velocity law, which contradicts a central assumption of the previously-developed XADM model of X-ray emission from magnetospheres. △ Less

Submitted 29 July, 2016; originally announced July 2016.

Comments: 18 pages, 12 figures. Accepted for publication in MNRAS

An `Analytic Dynamical Magnetosphere' formalism for X-ray and optical emission from slowly rotating magnetic massive stars

Authors: Stanley P. Owocki, Asif ud-Doula, Jon O. Sundqvist, Veronique Petit, David H. Cohen, Richard H. D. Townsend

Abstract: Slowly rotating magnetic massive stars develop "dynamical magnetospheres" (DM's), characterized by trapping of stellar wind outflow in closed magnetic loops, shock heating from collision of the upflow from opposite loop footpoints, and subsequent gravitational infall of radiatively cooled material. In 2D and 3D magnetohydrodynamic (MHD) simulations the interplay among these three components is spa… ▽ More Slowly rotating magnetic massive stars develop "dynamical magnetospheres" (DM's), characterized by trapping of stellar wind outflow in closed magnetic loops, shock heating from collision of the upflow from opposite loop footpoints, and subsequent gravitational infall of radiatively cooled material. In 2D and 3D magnetohydrodynamic (MHD) simulations the interplay among these three components is spatially complex and temporally variable, making it difficult to derive observational signatures and discern their overall scaling trends.Within a simplified, steady-state analysis based on overall conservation principles, we present here an "analytic dynamical magnetosphere" (ADM) model that provides explicit formulae for density, temperature and flow speed in each of these three components -- wind outflow, hot post-shock gas, and cooled inflow -- as a function of colatitude and radius within the closed (presumed dipole) field lines of the magnetosphere. We compare these scalings with time-averaged results from MHD simulations, and provide initial examples of application of this ADM model for deriving two key observational diagnostics, namely hydrogen H-alpha emission line profiles from the cooled infall, and X-ray emission from the hot post-shock gas. We conclude with a discussion of key issues and advantages in applying this ADM formalism toward derivation of a broader set of observational diagnostics and scaling trends for massive stars with such dynamical magnetospheres. △ Less

Submitted 28 July, 2016; originally announced July 2016.

Comments: 15 pages, 11 figures, accepted for MNRAS

Sub-Inertial Gravity Modes in the B8V Star KIC 7760680 Reveal Moderate Core Overshooting and Low Vertical Diffusive Mixing

Authors: Ehsan Moravveji, Richard H. D. Townsend, Conny Aerts, Stephane Mathis

Abstract: KIC 7760680 is so far the richest slowly pulsating B star, by exhibiting 36 consecutive dipole ($\ell=1$) gravity (g-) modes. The monotonically decreasing period spacing of the series, in addition to the local dips in the pattern confirm that KIC 7760680 is a moderate rotator, with clear mode trapping in chemically inhomogeneous layers. We employ the traditional approximation of rotation to incorp… ▽ More KIC 7760680 is so far the richest slowly pulsating B star, by exhibiting 36 consecutive dipole ($\ell=1$) gravity (g-) modes. The monotonically decreasing period spacing of the series, in addition to the local dips in the pattern confirm that KIC 7760680 is a moderate rotator, with clear mode trapping in chemically inhomogeneous layers. We employ the traditional approximation of rotation to incorporate rotational effects on g-mode frequencies. Our detailed forward asteroseismic modelling of this g-mode series reveals that KIC 7760680 is a moderately rotating B star with mass $\sim3.25$ M$_\odot$. By simultaneously matching the slope of the period spacing, and the number of modes in the observed frequency range, we deduce that the equatorial rotation frequency of KIC 7760680 is 0.4805 day$^{-1}$, which is 26\% of its Roche break up frequency. The relative deviation of the model frequencies and those observed is less than one percent. We succeed to tightly constrain the exponentially-decaying convective core overshooting parameter to $f_{\rm ov}\approx0.024\pm0.001$. This means that convective core overshooting can coexist with moderate rotation. Moreover, models with exponentially-decaying overshoot from the core outperform those with the classical step-function overshoot. The best value for extra diffusive mixing in the radiatively stable envelope is confined to $\log D_{\rm ext}\approx0.75\pm0.25$ (with $D_{\rm ext}$ in cm$^2$ sec$^{-1}$), which is notably smaller than theoretical predictions. △ Less

Submitted 10 April, 2016; originally announced April 2016.

Comments: 12 Figures, 2 Tables, all data publicly available for download; accepted for publication in Astrophysical Journal

The Variable Line Width of Achernar

Authors: Th. Rivinius, R. H. D. Townsend, D. Baade, A. C. Carciofi, N. Leister, S. Štefl

Abstract: 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 availa… ▽ More 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. △ Less

Submitted 10 February, 2016; originally announced February 2016.

Comments: To appear in conference proceedings "Bright Emissaries: Be Stars as Messengers of Star-Disk Physics", eds. A. Sigut and C. Jones (11 Aug 2014 - 13 Aug 2014, London, Ontario, Canada)

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)

Modules for Experiments in Stellar Astrophysics (MESA): Binaries, Pulsations, and Explosions

Authors: Bill Paxton, Pablo Marchant, Josiah Schwab, Evan B. Bauer, Lars Bildsten, Matteo Cantiello, Luc Dessart, R. Farmer, H. Hu, N. Langer, R. H. D. Townsend, Dean M. Townsley, F. X. Timmes

Abstract: We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of… ▽ More We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development. △ Less

Submitted 8 January, 2017; v1 submitted 9 June, 2015; originally announced June 2015.

Comments: Published in the ApJ Supplement Series; updated to account for published erratum

Journal ref: ApJS, 220, 15 (2015)

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

Authors: M. E. Oksala, O. Kochukhov, J. Krticka, R. H. D. Townsend, G. A. Wade, M. Prvak, Z. Mikulasek, J. Silvester, S. P. Owocki

Abstract: 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 d… ▽ More 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. △ Less

Submitted 18 May, 2015; originally announced May 2015.

Comments: 16 pages, 17 figures, accepted for publication in MNRAS

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

Authors: M. Shultz, Th. Rivinius, C. P. Folsom, G. A. Wade, R. H. D. Townsend, J. Sikora, J. Grunhut, O. Stahl, the MiMeS Collaboration

Abstract: 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 overa… ▽ More 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. △ Less

Submitted 13 April, 2015; originally announced April 2015.

Comments: 22 pages, 21 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 2015 449 (4): 3945-3965

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

Authors: M. S. Cunha, D. Stello, P. P. Avelino, J. Christensen-Dalsgaard, R. H. D. Townsend

Abstract: 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 descr… ▽ More 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. △ Less

Submitted 31 March, 2015; originally announced March 2015.

Comments: Accepted for publication in ApJ

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

Authors: V. Grinberg, M. A. Leutenegger, N. Hell, K. Pottschmidt, M. Böck, J. A. García, M. Hanke, M. A. Nowak, J. O. Sundqvist, R. H. D. Townsend, J. Wilms

Abstract: 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.… ▽ More 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. △ Less

Submitted 25 February, 2015; originally announced February 2015.

Comments: proposed for acceptance in A&A, 11 pages, 11 figures (two in appendix)

Journal ref: A&A 576, A117 (2015)

Rotation, spectral variability, magnetic geometry and magnetosphere of the Of?p star CPD -28 2561

Authors: G. A. Wade, R. H. Barbá, J. Grunhut, F. Martins, V. Petit, J. O. Sundqvist, R. H. D. Townsend, N. R. Walborn, E. Alecian, E. J. Alfaro, J. Ma\' iz Apellániz, J. I. Arias, R. Gamen, N. Morrell, Y. Nazé, A. ud-Doula, the MiMeS Collaboration

Abstract: We report magnetic and spectroscopic observations and modeling of the Of?p star CPD -28 2561. Using more than 75 new spectra, we have measured the equivalent width variations and examined the dynamic spectra of photospheric and wind-sensitive spectral lines. A period search results in an unambiguous 73.41 d variability period. High resolution spectropolarimetric data analyzed using Least-Squares D… ▽ More We report magnetic and spectroscopic observations and modeling of the Of?p star CPD -28 2561. Using more than 75 new spectra, we have measured the equivalent width variations and examined the dynamic spectra of photospheric and wind-sensitive spectral lines. A period search results in an unambiguous 73.41 d variability period. High resolution spectropolarimetric data analyzed using Least-Squares Deconvolution yield a Zeeman signature detected in the mean Stokes V profile corresponding to phase 0.5 of the spectral ephemeris. Interpreting the 73.41 d period as the stellar rotational period, we have phased the equivalent widths and inferred longitudinal field measurements. The phased magnetic data exhibit a weak sinusoidal variation, with maximum of about 565 G at phase 0.5, and a minimum of about -335 G at phase 0.0, with extrema approximately in phase with the (double-wave) Halpha equivalent width variation. Modeling of the Halpha equivalent width variation assuming a quasi-3D magnetospheric model produces a unique solution for the ambiguous couplet of inclination and magnetic obliquity angles: (i, beta) or (beta, i)=(35 deg,90 deg). Adopting either geometry, the longitudinal field variation yields a dipole polar intensity Bd=2.6\pm 0.9~kG, consistent with that obtained from direct modelling of the Stokes V profiles. We derive a wind magnetic confinement parameter eta*\simeq 100, leading to an Alfvén radius RA\simeq 3-5~R*, and a Kepler radius RK\simeq 20~R*. This supports a physical scenario in which the Halpha emission and other line variability have their origin in an oblique, co-rotating 'dynamical magnetosphere' structure resulting from a magnetically channeled wind. Nevertheless, the details of the formation of spectral lines and their variability within this framework remain generally poorly understood. △ Less

Submitted 9 April, 2015; v1 submitted 28 November, 2014; originally announced November 2014.

Comments: 18 pages, accepted by MNRAS Replaced 28 March 2015 to include corrected figure 10 (see MNRAS erratum to this effect)

Journal ref: MNRAS, 447, 2551-2567 (2015)

Variable rotational line broadening in the Be star Achernar

Authors: Th. Rivinius, D. Baade, R. H. D. Townsend, A. C. Carciofi, S. Štefl

Abstract: 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… ▽ More 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. △ Less

Submitted 27 September, 2013; originally announced September 2013.

Comments: Accepted by A&A as letter

GYRE: An open-source stellar oscillation code based on a new Magnus Multiple Shooting Scheme

Authors: R. H D. Townsend, S. Teitler

Abstract: We present a new oscillation code, GYRE, which solves the stellar pulsation equations (both adiabatic and non-adiabatic) using a novel Magnus Multiple Shooting numerical scheme devised to overcome certain weaknesses of the usual relaxation and shooting schemes appearing in the literature. The code is accurate (up to 6th order in the number of grid points), robust, efficiently makes use of multiple… ▽ More We present a new oscillation code, GYRE, which solves the stellar pulsation equations (both adiabatic and non-adiabatic) using a novel Magnus Multiple Shooting numerical scheme devised to overcome certain weaknesses of the usual relaxation and shooting schemes appearing in the literature. The code is accurate (up to 6th order in the number of grid points), robust, efficiently makes use of multiple processor cores and/or nodes, and is freely available in source form for use and distribution. We verify the code against analytic solutions and results from other oscillation codes, in all cases finding good agreement. Then, we use the code to explore how the asteroseismic observables of a 1.5 M_sun star change as it evolves through the red-giant bump. △ Less

Submitted 13 August, 2013; originally announced August 2013.

Comments: 14 pages, 7 figures, accepted by MNRAS

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

Authors: R. H. D. Townsend, Th. Rivinius, J. F. Rowe, A. F. J. Moffat, J. M. Matthews, D. Bohlender, C. Neiner, J. H. Telting, D. B. Guenther, T. Kallinger, R. Kuschnig, S. M. Rucinski, D. Sasselov, W. W. Weiss

Abstract: 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, eit… ▽ More 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. △ Less

Submitted 8 April, 2013; originally announced April 2013.

Comments: 6 pages, 2 figures, accepted by the Astrophysical Journal

Polarimetric observations of $σ$ Orionis E

Authors: Alex C. Carciofi, Daniel M. Faes, Richard H. D. Townsend, Jon E Bjorkman

Abstract: Some massive stars possess strong magnetic fields that confine plasma in the circumstellar environment. These \textit{magnetospheres} have been studied spectroscopically, photometrically and, more recently, interferometrically. Here we report on the first firm detection of a magnetosphere in continuum linear polarization, as a result of monitoring of $σ$\,Ori\,E at the Pico dos Dias Observatory. {… ▽ More Some massive stars possess strong magnetic fields that confine plasma in the circumstellar environment. These \textit{magnetospheres} have been studied spectroscopically, photometrically and, more recently, interferometrically. Here we report on the first firm detection of a magnetosphere in continuum linear polarization, as a result of monitoring of $σ$\,Ori\,E at the Pico dos Dias Observatory. {The non-zero intrinsic polarization indicates an asymmetric structure, whose minor elongation axis is oriented $150\fdg0$ east of the celestial north.} A modulation of the polarization was observed, with a period of half of the rotation period, which supports the theoretical prediction of the presence of two diametrally opposed, co-rotating blobs of gas. A phase lag of -0.085 was detected between the polarization minimum and the primary minimum of the light curve, suggestive of a complex shape of the plasma clouds. We present a preliminary analysis of the data with the Rigidly Rotating Magnetosphere model, which could not reproduce simultaneously the photometric and polarimetric data. A toy model comprising two spherical co-rotating blobs {joined by a thin disk} proved more successful in reproducing the polarization modulation. {With this model we were able to determine that the total scattering mass of the thin disk is similar to the mass of the blobs ($2M_{\rm b}/M_{\rm d}=1.2$) and that the blobs are rotating counterclockwise on the plane of the sky.} This result shows that polarimetry can provide a diagnostic of the geometry of clouds, which will serve as an important constraint for {improving} the Rigidly Rotating Magnetosphere model. △ Less

Submitted 19 February, 2013; originally announced February 2013.

Comments: 5 pages, 4 figures, accepted for The Astrophysical Journal Letters

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

Authors: V. Petit, S. P. Owocki, G. A. Wade, D. H. Cohen, J. O. Sundqvist, M. Gagné, J. Maíz Apellániz, M. E. Oksala, D. A. Bohlender, Th. Rivinius, H. F. Henrichs, E. Alecian, R. H. D. Townsend, A. ud-Doula, the MiMeS Collaboration

Abstract: 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 compi… ▽ More 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] △ Less

Submitted 1 November, 2012; originally announced November 2012.

Comments: Accepted for publication in MNRAS 30 pages, 10 figures, 6 tables