Search

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., Comment: Paper accepted for publication in Astronomy & Astrophysics (10 pages, 11 figures). Full abstract in PDF

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., Comment: 3 pages, 1 figure

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., Comment: 8 pages, 6 figures, accepted by ApJ

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., Comment: 24 pages, 13 figures, 1 table. Accepted by ApJ

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., Comment: 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)

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 $\alpha$ 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 $\alpha$ 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., Comment: Accepted to The Astrophysical Journal; 15 pages, 10 figures

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}$., Comment: 17 pages, 16 figures, accepted for publication in ApJ

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}}$ $\beta$ Cephei star models, and of a $0.9\,\rm{M_{\odot}}$ extreme helium star model., Comment: 12 pages, 14 figures, Accepted to ApJ

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., Comment: 10 pages, 2 figures, accepted by MNRAS

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., Comment: Astronomy and Astrophysics, in the press. Paper resubmitted to arXiv to correct author list in metadata

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$, $\Omega/\Omega_{\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$, $\Omega/\Omega_{\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., Comment: accepted for publication in MNRAS, a zenodo record is available at: https://zenodo.org/record/3250412

(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., Comment: accepted for publication in Astronomy and Astrophysics

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., Comment: 11 pages, 3 figures, 1 appendix

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., Comment: 15 pages, 13 figures. Accepted for publication in Astronomy and Astrophysics

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., Comment: Accepted for publication in MNRAS (13 pages, 8 figures)

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., Comment: Accepted for publication in A&A

Recent spectropolarimetric surveys of bright, hot stars have found that ~10% of OB-type stars contain strong (mostly dipolar) surface magnetic fields (~kG). The prominent paradigm describing the interaction between the stellar winds and the surface magnetic field is the magnetically confined wind shock (MCWS) model. In this model, the stellar wind plasma is forced to move along the closed field loops of the magnetic field, colliding at the magnetic equator, and creating a shock. As the shocked material cools radiatively it will emit X-rays. Therefore, X-ray spectroscopy is a key tool in detecting and characterizing the hot wind material confined by the magnetic fields of these stars. Some B-type stars are found to have very short rotational periods. The effects of the rapid rotation on the X-ray production within the magnetosphere have yet to be explored in detail. The added centrifugal force due to rapid rotation is predicted to cause faster wind outflows along the field lines, leading to higher shock temperatures and harder X-rays. However, this is not observed in all rapidly rotating magnetic B-type stars. In order to address this from a theoretical point of view, we use the X-ray Analytical Dynamical Magnetosphere (XADM) model, originally developed for slow rotators, with an implementation of new rapid rotational physics. Using X-ray spectroscopy from ESA's XMM-Newton space telescope, we observed 5 rapidly rotating B-type stars to add to the previous list of observations. Comparing the observed X-ray luminosity and hardness ratio to that predicted by the XADM allows us to determine the role the added centrifugal force plays in the magnetospheric X-ray emission of these stars., Comment: IAUS Conference Proceedings

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., Comment: 10 pages, 7 figures, 1 table. Accepted for publication in MNRAS

Magnetic B-type stars exhibit photometric variability due to diverse causes, and consequently on a variety of timescales. In this paper we describe interpretation of BRITE photometry and related ground-based observations of 4 magnetic B-type systems: $\epsilon$ Lupi, $\tau$ Sco, a Cen and $\epsilon$ CMa., Comment: 7 pages, proceedings of the 3rd BRITE Science Meeting

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 $\pi$ Aqr, are presented as well., Comment: 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 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\'elescope 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., Comment: Accepted by MNRAS

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., Comment: Published in the ApJ Supplement Series; updated to account for published erratum

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., Comment: 16 pages, 17 figures, accepted for publication in MNRAS

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., Comment: 22 pages, 21 figures, published in MNRAS

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 \'echelle 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., Comment: Accepted for publication in ApJ

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., Comment: proposed for acceptance in A&A, 11 pages, 11 figures (two in appendix)

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\'en 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., Comment: 18 pages, accepted by MNRAS Replaced 28 March 2015 to include corrected figure 10 (see MNRAS erratum to this effect)

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\alpha 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 \Delta 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\alpha 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\alpha line emission becomes undetectable., Comment: Accepted by A&A as letter

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., Comment: 14 pages, 7 figures, accepted by MNRAS