Showing 1–50 of 65 results for author: Johnston, C

Binarity at LOw Metallicity (BLOeM): I. a spectroscopic VLT monitoring survey of massive stars in the SMC

Authors: T. Shenar, J. Bodensteiner, H. Sana, P. A. Crowther, D. J. Lennon, M. Abdul-Masih, L. A. Almeida, F. Backs, S. R. Berlanas, M. Bernini-Peron, J. M. Bestenlehner, D. M. Bowman, V. A. Bronner, N. Britavskiy, A. de Koter, S. E. de Mink, K. Deshmukh, C. J. Evans, M. Fabry, M. Gieles, A. Gilkis, G. González-Torà, G. Gräfener, Y. Götberg, C. Hawcroft , et al. (52 additional authors not shown)

Abstract: Surveys in the Milky Way and Large Magellanic Cloud revealed that the majority of massive stars will interact with companions during their lives. However, knowledge of the binary properties of massive stars at low metallicity, which approaches the conditions of the Early Universe, remains sparse. We present the Binarity at LOw Metallicity (BLOeM) campaign - an ESO large programme designed to obtai… ▽ More Surveys in the Milky Way and Large Magellanic Cloud revealed that the majority of massive stars will interact with companions during their lives. However, knowledge of the binary properties of massive stars at low metallicity, which approaches the conditions of the Early Universe, remains sparse. We present the Binarity at LOw Metallicity (BLOeM) campaign - an ESO large programme designed to obtain 25 epochs of spectroscopy for 929 massive stars in the SMC - the lowest metallicity conditions in which multiplicity is probed to date (Z = 0.2 Zsun). BLOeM will provide (i) the binary fraction, (ii) the orbital configurations of systems with periods P < 3 yr, (iii) dormant OB+BH binaries, and (iv) a legacy database of physical parameters of massive stars at low metallicity. The stars are observed with the LR02 setup of the giraffe instrument of the Very Large Telescope (3960-4570A, resolving power R=6200; typical signal-to-noise ratio S/N=70-100). This paper utilises the first 9 epochs obtained over a three-month time. We describe the survey and data reduction, perform a spectral classification of the stacked spectra, and construct a Hertzsprung-Russell diagram of the sample via spectral-type and photometric calibrations. The sample covers spectral types from O4 to F5, spanning the effective temperature and luminosity ranges 6.5<Teff/kK<45 and 3.7<log L/Lsun<6.1 and initial masses 8<Mini/Msun<80. It comprises 159 O-type stars, 324 early B-type (B0-3) dwarfs and giants (luminosity classes V-III), 309 early B-type supergiants (II-I), and 137 late-type supergiants. At least 75 stars are Oe/Be stars: 20 O-type and 55 B-type (13% and 10% of the respective samples). In addition, it includes four high-mass X-ray binaries, three stars resembling luminous blue variables, two bloated stripped-star candidates, two candidate magnetic stars, and 74 eclipsing binaries. △ Less

Submitted 19 July, 2024; originally announced July 2024.

Comments: Submitted to A&A. Comments welcome

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

The BlackGEM telescope array I: Overview

Authors: Paul J. Groot, S. Bloemen, P. Vreeswijk, J. van Roestel, P. G. Jonker, G. Nelemans, M. Klein-Wolt, R. Le Poole, D. Pieterse, M. Rodenhuis, W. Boland, M. Haverkorn, C. Aerts, R. Bakker, H. Balster, M. Bekema, E. Dijkstra, P. Dolron, E. Elswijk, A. van Elteren, A. Engels, M. Fokker, M. de Haan, F. Hahn, R. ter Horst , et al. (49 additional authors not shown)

Abstract: The main science aim of the BlackGEM array is to detect optical counterparts to gravitational wave mergers. Additionally, the array will perform a set of synoptic surveys to detect Local Universe transients and short time-scale variability in stars and binaries, as well as a six-filter all-sky survey down to ~22nd mag. The BlackGEM Phase-I array consists of three optical wide-field unit telescopes… ▽ More The main science aim of the BlackGEM array is to detect optical counterparts to gravitational wave mergers. Additionally, the array will perform a set of synoptic surveys to detect Local Universe transients and short time-scale variability in stars and binaries, as well as a six-filter all-sky survey down to ~22nd mag. The BlackGEM Phase-I array consists of three optical wide-field unit telescopes. Each unit uses an f/5.5 modified Dall-Kirkham (Harmer-Wynne) design with a triplet corrector lens, and a 65cm primary mirror, coupled with a 110Mpix CCD detector, that provides an instantaneous field-of-view of 2.7~square degrees, sampled at 0.564\arcsec/pixel. The total field-of-view for the array is 8.2 square degrees. Each telescope is equipped with a six-slot filter wheel containing an optimised Sloan set (BG-u, BG-g, BG-r, BG-i, BG-z) and a wider-band 440-720 nm (BG-q) filter. Each unit telescope is independent from the others. Cloud-based data processing is done in real time, and includes a transient-detection routine as well as a full-source optimal-photometry module. BlackGEM has been installed at the ESO La Silla observatory as of October 2019. After a prolonged COVID-19 hiatus, science operations started on April 1, 2023 and will run for five years. Aside from its core scientific program, BlackGEM will give rise to a multitude of additional science cases in multi-colour time-domain astronomy, to the benefit of a variety of topics in astrophysics, such as infant supernovae, luminous red novae, asteroseismology of post-main-sequence objects, (ultracompact) binary stars, and the relation between gravitational wave counterparts and other classes of transients △ Less

Submitted 30 May, 2024; v1 submitted 29 May, 2024; originally announced May 2024.

Comments: 14 pages, submitted to Astronomy & Astrophysics

Planet Hunters TESS V: a planetary system around a binary star, including a mini-Neptune in the habitable zone

Authors: Nora L. Eisner, Samuel K. Grunblatt, Oscar Barragán, Thea H. Faridani, Chris Lintott, Suzanne Aigrain, Cole Johnston, Ian R. Mason, Keivan G. Stassun, Megan Bedell, Andrew W. Boyle, David R. Ciardi, Catherine A. Clark, Guillaume Hebrard, David W. Hogg, Steve B. Howell, Baptiste Klein, Joe Llama, Joshua N. Winn, Lily L. Zhao, Joseph M. Akana Murphy, Corey Beard, Casey L. Brinkman, Ashley Chontos, Pia Cortes-Zuleta , et al. (39 additional authors not shown)

Abstract: We report on the discovery and validation of a transiting long-period mini-Neptune orbiting a bright (V = 9.0 mag) G dwarf (TOI 4633; R = 1.05 RSun, M = 1.10 MSun). The planet was identified in data from the Transiting Exoplanet Survey Satellite by citizen scientists taking part in the Planet Hunters TESS project. Modeling of the transit events yields an orbital period of 271.9445 +/- 0.0040 days… ▽ More We report on the discovery and validation of a transiting long-period mini-Neptune orbiting a bright (V = 9.0 mag) G dwarf (TOI 4633; R = 1.05 RSun, M = 1.10 MSun). The planet was identified in data from the Transiting Exoplanet Survey Satellite by citizen scientists taking part in the Planet Hunters TESS project. Modeling of the transit events yields an orbital period of 271.9445 +/- 0.0040 days and radius of 3.2 +/- 0.20 REarth. The Earth-like orbital period and an incident flux of 1.56 +/- 0.2 places it in the optimistic habitable zone around the star. Doppler spectroscopy of the system allowed us to place an upper mass limit on the transiting planet and revealed a non-transiting planet candidate in the system with a period of 34.15 +/- 0.15 days. Furthermore, the combination of archival data dating back to 1905 with new high angular resolution imaging revealed a stellar companion orbiting the primary star with an orbital period of around 230 years and an eccentricity of about 0.9. The long period of the transiting planet, combined with the high eccentricity and close approach of the companion star makes this a valuable system for testing the formation and stability of planets in binary systems. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 24 pages, 16 figures, 4 tables

Journal ref: Published in AJ, 2024

The Asteroseismic Imprints of Mass Transfer: A Case Study of a Binary Mass Gainer in the SPB Instability Strip

Authors: Tom Wagg, Cole Johnston, Earl P. Bellinger, Mathieu Renzo, Richard Townsend, Selma E. de Mink

Abstract: We present new simulations investigating the impact of mass transfer on the asteroseismic signals of slowly pulsating B stars. We use MESA to simulate the evolution of a binary star system and GYRE to compute the asteroseismic properties of the accretor star. We show that, compared to a single star of the same final mass, a star that has undergone accretion (of non-enriched material) has a signifi… ▽ More We present new simulations investigating the impact of mass transfer on the asteroseismic signals of slowly pulsating B stars. We use MESA to simulate the evolution of a binary star system and GYRE to compute the asteroseismic properties of the accretor star. We show that, compared to a single star of the same final mass, a star that has undergone accretion (of non-enriched material) has a significantly different internal structure, evident in both the hydrogen abundance profile and Brunt-Väisälä frequency profile. These differences result in significant changes in the observed period spacing patterns, implying that one may use this as a diagnostic to test whether a star's core has been rejuvenated as a result of accretion. We show that it is essential to consider the full multimodal posterior distributions when fitting stellar properties of mass-gainers to avoid drawing misleading conclusions. Even with these considerations, stellar ages will be significantly underestimated when assuming single star evolution for a mass-gainer. We find that future detectors with improved uncertainties would rule out single star models with the correct mass and central hydrogen fraction. Our proof of principle analysis demonstrates the need to further investigate the impact of binary interactions on stellar asteroseismic signals for a wide range of parameters, such as initial mass, amount of mass transferred and the age of the accretor star at the onset of mass transfer. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Comments: 10 pages (main text), 8 figures, submitted to A&A, comments welcome. Interactive plots available at https://www.tomwagg.com/html/interact/mass-gainer-asteroseismology.html

Journal ref: A&A 687, A222 (2024)

Advancing parabolic operators in thermodynamic MHD models II: Evaluating a Practical Time Step Limit for Unconditionally Stable Methods

Authors: Ronald M. Caplan, Craig D. Johnston, Lars K. S. Daldoff, Jon A. Linker

Abstract: Unconditionally stable time stepping schemes are useful and often practically necessary for advancing parabolic operators in multi-scale systems. However, serious accuracy problems may emerge when taking time steps that far exceed the explicit stability limits. In our previous work, we compared the accuracy and performance of advancing parabolic operators in a thermodynamic MHD model using an impl… ▽ More Unconditionally stable time stepping schemes are useful and often practically necessary for advancing parabolic operators in multi-scale systems. However, serious accuracy problems may emerge when taking time steps that far exceed the explicit stability limits. In our previous work, we compared the accuracy and performance of advancing parabolic operators in a thermodynamic MHD model using an implicit method and an explicit super time-stepping (STS) method. We found that while the STS method outperformed the implicit one with overall good results, it was not able to damp oscillatory behavior in the solution efficiently, hindering its practical use. In this follow-up work, we evaluate an easy-to-implement method for selecting a practical time step limit (PTL) for unconditionally stable schemes. This time step is used to `cycle' the operator-split thermal conduction and viscosity parabolic operators. We test the new time step with both an implicit and STS scheme for accuracy, performance, and scaling. We find that, for our test cases here, the PTL dramatically improves the STS solution, matching or improving the solution of the original implicit scheme, while retaining most of its performance and scaling advantages. The PTL shows promise to allow more accurate use of unconditionally stable schemes for parabolic operators and reliable use of STS methods. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 12 pages, 6 figures. ASTRONUM 23

Automated eccentricity measurement from raw eclipsing binary light curves with intrinsic variability

Authors: Luc W. IJspeert, Andrew Tkachenko, Cole Johnston, Andrej Prša, Mark A. Wells, Conny Aerts

Abstract: Abstract abridged. Eclipsing binary systems provide the opportunity to measure the fundamental parameters of their component stars in a stellar-model-independent way. This makes them ideal candidates for testing and calibrating theories of stellar structure and (tidal) evolution. Even without spectroscopic follow-up there is often enough information in their photometric time series to warrant anal… ▽ More Abstract abridged. Eclipsing binary systems provide the opportunity to measure the fundamental parameters of their component stars in a stellar-model-independent way. This makes them ideal candidates for testing and calibrating theories of stellar structure and (tidal) evolution. Even without spectroscopic follow-up there is often enough information in their photometric time series to warrant analysis, especially if there is an added value present in the form of intrinsic variability, such as pulsations. Our goal is to implement and validate a framework for the homogeneous analysis of large numbers of eclipsing binary light curves, such as the numerous high-duty-cycle observations from space missions like TESS. The aim of this framework is to be quick and simple to run and to limit the user's time investment when obtaining, amongst other parameters, orbital eccentricities. We developed a new and fully automated methodology for the analysis of eclipsing binary light curves with or without additional intrinsic variability. Our method includes a fast iterative pre-whitening procedure. Orbital and stellar parameters are measured under the assumption of spherical stars of uniform brightness. We tested our methodology in two settings: a set of synthetic light curves with known input and the catalogue of Kepler eclipsing binaries. The synthetic tests show that we can reliably recover the frequencies and amplitudes of the sinusoids included in the signal as well as the input binary parameters. Recovery of the tangential component of eccentricity is the most accurate and precise. Kepler results confirm a robust determination of orbital periods, with 80.5% of periods matching the catalogued ones. We present the eccentricities for this analysis and show that they broadly follow the theoretically expected pattern as a function of the orbital period. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 17 pages, 17 figures, 1 table, 5 appendices (43 additional figures, 1 additional table). Accepted for publication in A&A

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)

FINKER: Frequency Identification through Nonparametric KErnel Regression in astronomical time series

Authors: F. Stoppa, C. Johnston, E. Cator, G. Nelemans, P. J. Groot

Abstract: Optimal frequency identification in astronomical datasets is crucial for variable star studies, exoplanet detection, and asteroseismology. Traditional period-finding methods often rely on specific parametric assumptions, employ binning procedures, or overlook the regression nature of the problem, limiting their applicability and precision. We aim to introduce a universal, nonparametric kernel regr… ▽ More Optimal frequency identification in astronomical datasets is crucial for variable star studies, exoplanet detection, and asteroseismology. Traditional period-finding methods often rely on specific parametric assumptions, employ binning procedures, or overlook the regression nature of the problem, limiting their applicability and precision. We aim to introduce a universal, nonparametric kernel regression method for optimal frequency determination that is generalizable, efficient, and robust across various astronomical data types. FINKER uses nonparametric kernel regression on folded datasets at different frequencies, selecting the optimal frequency by minimizing squared residuals. This technique inherently incorporates a weighting system that accounts for measurement uncertainties and facilitates multiband data analysis. We evaluate our method's performance across a range of frequencies pertinent to diverse data types and compare it with an established period-finding algorithm, conditional entropy. The method demonstrates superior performance in accuracy and robustness compared to existing algorithms, requiring fewer observations to identify significant frequencies reliably. It exhibits resilience against noise and adapts well to datasets with varying complexity. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Journal ref: A&A 686, A158 (2024)

Variability of Blue Supergiants in the LMC with TESS

Authors: Linhao Ma, Cole Johnston, Earl P. Bellinger, Selma E. de Mink

Abstract: The blue supergiant (BSG) problem, namely, the overabundance of BSGs inconsistent with classical stellar evolution theory, remains an open question in stellar astrophysics. Several theoretical explanations have been proposed, which may be tested by their predictions for the characteristic time variability. In this work, we analyze the light curves of a sample of 20 BSGs obtained from the Transitin… ▽ More The blue supergiant (BSG) problem, namely, the overabundance of BSGs inconsistent with classical stellar evolution theory, remains an open question in stellar astrophysics. Several theoretical explanations have been proposed, which may be tested by their predictions for the characteristic time variability. In this work, we analyze the light curves of a sample of 20 BSGs obtained from the Transiting Exoplanet Survey Satellite (TESS) mission. We report a characteristic signal in the low-frequency ($f\lesssim2\;\mathrm{day}^{-1}$) range for all our targets. The amplitude spectrum has a peak frequency at $\sim0.2\;\mathrm{day}^{-1}$, and we are able to fit it by a modified Lorentzian profile. The signal itself shows strong stochasticity across different TESS sectors, suggesting its driving mechanism happens on short ($\lesssim\mathrm{months}$) timescales. Our signals resemble those obtained for a limited sample of hotter OB stars and yellow supergiants, suggesting their possible common origins. We discuss three possible physical explanations: stellar winds launched by rotation, convection motions that reach the stellar surface, and waves from the deep stellar interior. The peak frequency of the signal favors processes related to the convective zone caused by the iron opacity peak, and the shape of the spectra might be explained by the propagation of high-order, damped gravity waves excited from that zone. We discuss the uncertainties and limitations of all these scenarios. △ Less

Submitted 28 April, 2024; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: Accepted to ApJ

Heating and cooling in stellar coronae: coronal rain on a young Sun

Authors: Simon Daley-Yates, Moira M. Jardine, Craig D. Johnston

Abstract: Recent observations of rapidly-rotating cool dwarfs have revealed H$α$ line asymmetries indicative of clumps of cool, dense plasma in the stars' coronae. These clumps may be either long-lived (persisting for more than one stellar rotation) or dynamic. The fastest dynamic features show velocities greater than the escape speed, suggesting that they may be centrifugally ejected from the star, contrib… ▽ More Recent observations of rapidly-rotating cool dwarfs have revealed H$α$ line asymmetries indicative of clumps of cool, dense plasma in the stars' coronae. These clumps may be either long-lived (persisting for more than one stellar rotation) or dynamic. The fastest dynamic features show velocities greater than the escape speed, suggesting that they may be centrifugally ejected from the star, contributing to the stellar angular momentum loss. Many however show lower velocities, similar to coronal rain observed on the Sun. We present 2.5D magnetohydrodynamic simulations of the formation and dynamics of these condensations in a rapidly rotating ($P_{\rm rot}~=~ 1 \ \mathrm{day}$) young Sun. Formation is triggered by excess surface heating. This pushes the system out of thermal equilibrium and triggers a thermal instability. The resulting condensations fall back towards the surface. They exhibit quasi-periodic behaviour, with periods longer than typical periods for solar coronal rain. We find line-of-sight velocities for these clumps in the range $50 \ \mathrm{km} \ \mathrm{s}^{-1}$ (blue shifted) to $250 \ \mathrm{km} \ \mathrm{s}^{-1}$ (red shifted). These are typical of those inferred from stellar H$α$ line asymmetries, but the inferred clump masses of $3.6\times 10^{14}\ \mathrm{g}$ are significantly smaller. We find that a maximum of $\simeq~3\%$ of the coronal mass is cool clumps. We conclude that coronal rain may be common in solar like stars, but may appear on much larger scales in rapid rotators. △ Less

Submitted 8 September, 2023; originally announced September 2023.

Comments: 11 pages, 5 figures

The Thickness of Electric Current Sheets and Implications for Coronal Heating

Authors: James A. Klimchuk, James E. Leake, Lars K. S. Daldorff, Craig D. Johnston

Abstract: The thickness of current sheets is extremely important, especially as it relates to the onset of fast magnetic reconnection. Onset determines how much magnetic free energy can build up in a field before it is explosively released. This has implications for many phenomena on the Sun and throughout the universe, including the heating of the solar corona. Significant effort has been devoted to the qu… ▽ More The thickness of current sheets is extremely important, especially as it relates to the onset of fast magnetic reconnection. Onset determines how much magnetic free energy can build up in a field before it is explosively released. This has implications for many phenomena on the Sun and throughout the universe, including the heating of the solar corona. Significant effort has been devoted to the question of whether equilibrium current sheets in realistic geometries have finite or zero thickness. Using a simple force balance analysis, we show why current sheets without a guide field (2D) and with a guide field that is invariant in the guide field direction (2.5D) cannot be in equilibrium if they have both finite thickness and finite length. We then estimate the conditions under which the tension of a curved line-tied guide field can facilitate equilibrium in 3D sheets that are finite in all dimensions. Finally, we argue that some quasi-statically evolving current sheets undergoing slow stressing (e.g., when the coronal magnetic field is subjected to photospheric boundary driving) may reach a critical shear, at which point they lose equilibrium, spontaneously collapse, and reconnect. The critical shear is generally consistent with the heating requirements of solar active regions. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 15 pages, 5 figures, accepted by Frontiers in Physics

Identifying and characterising the population of hot sub-luminous stars with multi-colour MeerLICHT data

Authors: P. Ranaivomanana, C. Johnston, P. J. Groot, C. Aerts, R. Lees, L. IJspeert, S. Bloemen, M. Klein-Wolt, P. Woudt, E. Kording, R. Le Poole, D. Pieterse

Abstract: Colour-magnitude diagrams reveal a population of blue (hot) sub-luminous objects with respect to the main sequence. These hot sub-luminous stars are the result of evolutionary processes that require stars to expel their obscuring, hydrogen-rich envelopes to reveal the hot helium core. As such, these objects offer a direct window into the hearts of stars that are otherwise inaccessible to direct ob… ▽ More Colour-magnitude diagrams reveal a population of blue (hot) sub-luminous objects with respect to the main sequence. These hot sub-luminous stars are the result of evolutionary processes that require stars to expel their obscuring, hydrogen-rich envelopes to reveal the hot helium core. As such, these objects offer a direct window into the hearts of stars that are otherwise inaccessible to direct observation. We showcase MeerLICHT's capabilities of detecting faint hot subdwarfs and identifying the dominant frequency in the photometric variability of these compact hot stars, in comparison to their $Gaia$ DR3 data. We hunt for oscillations, which will be an essential ingredient for accurately probing stellar interiors in future asteroseismology. Comparative MeerLICHT and $Gaia$ colour-magnitude diagrams are presented as a way to select hot subdwarfs from our sample. A dedicated frequency determination technique is developed and applied to the selected candidates to determine their dominant variability using time-series data from MeerLICHT and $Gaia$ DR3. We explore the power of both datasets in determining the dominant frequency. Using the $g-i$ colour, MeerLICHT offers a colour-magnitude diagram that is comparable in quality to that of $Gaia$ DR3. The MeerLICHT colour-colour diagrams allow for the study of different stellar populations. The frequency analysis of MeerLICHT and $Gaia$ DR3 data demonstrates the superiority of our MeerLICHT multi-colour photometry in estimating the dominant frequency compared to the sparse $Gaia$ DR3 data. MeerLICHT's multi-band photometry leads to the discovery of high-frequency faint subdwarfs. Our MeerLICHT results are a proof-of-concept of the capacity of the BlackGEM instrument currently in the commissioning stage at ESO's La Silla Observatory in Chile. △ Less

Submitted 14 February, 2023; originally announced February 2023.

Comments: Accepted from A&A journal. Publication is still in process

MSC Class: 62M10; 65T50 (primary)

Journal ref: A&A 672, A69 (2023)

Tidally perturbed g-mode pulsations in a sample of close eclipsing binaries

Authors: T. Van Reeth, C. Johnston, J. Southworth, J. Fuller, D. M. Bowman, L. Poniatowski, J. Van Beeck

Abstract: Context. Thanks to the high-precision photometry from space missions such as Kepler and TESS, tidal perturbations and tilting of pulsations have been detected in more than a dozen binary systems. However, only two of these were g-mode pulsators. Aims. We aim to detect tidally perturbed g modes in additional binary systems and characterise them observationally. Methods. We perform a custom data red… ▽ More Context. Thanks to the high-precision photometry from space missions such as Kepler and TESS, tidal perturbations and tilting of pulsations have been detected in more than a dozen binary systems. However, only two of these were g-mode pulsators. Aims. We aim to detect tidally perturbed g modes in additional binary systems and characterise them observationally. Methods. We perform a custom data reduction of the available Kepler and TESS photometry of a well-studied sample of 35 binary systems with gamma Doradus pulsators. For each target, we model the binary signal using a sum of 100 sine waves, with frequencies at orbital harmonics, and measure significant pulsation frequencies by iteratively prewhitening the residual light curve. Pulsations are labelled as tidally perturbed g modes if they are part of both period-spacing patterns and orbital-frequency-spaced multiplets. After visual inspection and confirmation, the properties of these targets and g modes are characterised. Results. We detect tidally perturbed g-mode pulsations for five short-period binaries that are circularised and (almost) synchronously rotating: KIC3228863, KIC3341457, KIC4947528, KIC9108579, and KIC12785282. Tidally perturbed g modes that occur within the same star and have the same mode identification (k,m), are found to have near-identical relative amplitude and phase modulations, which are within their respective 1-sigma uncertainties also identical for the Kepler and TESS photometric passbands. By contrast, pulsations with different mode identification (k,m) are found to exhibit different modulations. Moreover, the observed amplitude and phase modulations are correlated, indicating that the binary tides primarily distort g-mode amplitudes on the stellar surface. The phase modulations are then primarily a geometric effect of the integration of the stellar flux over the visible stellar surface. (abbreviated) △ Less

Submitted 20 January, 2023; originally announced January 2023.

Comments: 21 pages, 33 figures, 2 tables. Accepted for publication in Astronomy and Astrophysics on 17 January 2023

Journal ref: A&A 671, A121 (2023)

Tidal perturbations and eclipse mapping in the pulsations in the hierarchical triple system U~Gru

Authors: Cole Johnston, Andrew Tkachenko, Timothy Van Reeth, Dominic M. Bowman, Kresimir Pavlovski, Hugues Sana, Sanjay Sekaran

Abstract: Context. Unambiguous examples of the influence of tides on self-excited, free stellar pulsations have recently been observationally detected in space-based photometric data. Aims. We aim to investigate U Gru and contextualise it within the growing class of tidally influenced pulsators. Initial analysis of U Gru revealed frequencies spaced by the orbital frequency that are difficult to explain by c… ▽ More Context. Unambiguous examples of the influence of tides on self-excited, free stellar pulsations have recently been observationally detected in space-based photometric data. Aims. We aim to investigate U Gru and contextualise it within the growing class of tidally influenced pulsators. Initial analysis of U Gru revealed frequencies spaced by the orbital frequency that are difficult to explain by currently proposed tidal mechanisms. Methods. We re-investigate the TESS photometry of U Gru alongside new uves spectroscopy. We analyse the uves spectroscopy with least-squares deconvolution and spectral disentangling techniques, and perform an atmospheric analysis. We remove the binary signature from the light curve using an effective model in order to investigate the pulsation signal in the residuals. We track the amplitudes and phases of the residual pulsations as a function of the orbital period to reveal their tidal influence. Results. We establish that U Gru is likely a hierarchical triple system. We identify a single p mode oscillation to exhibit amplitude and phase variation over the binary orbit. We propose a toy model to demonstrate that the series of frequencies separated by the orbital frequency can be reproduced by eclipse mapping. We find no evidence of modulation to the other independent oscillation modes. Conclusions. We demonstrate that U Gru hosts at least one tidally perturbed pulsation. Additionally we argue that eclipse mapping of the dominant, tidally perturbed mode can produce the series of frequencies separated by the observed orbital frequency. Our simulations show that the effects of eclipse mapping are mode dependent, and are not expected to produce an observable signature for all pulsation modes in an eclipse binary. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: Accepted for publication in Astronomy & Astrophysics. 13 pages, 10 figures, 3 tables, 1 appendix

Journal ref: A&A 670, A167 (2023)

Planet Hunters TESS IV: A massive, compact hierarchical triple star system TIC 470710327

Authors: Nora L. Eisner, Cole Johnston, Silvia Toonen, Abigail J. Frost, Soetkin Janssens, Chris J. Lintott, Suzanne Aigrain, Hugues Sana, Michael Abdul-Masih, Karla Z. Arellano-Córdova, Paul G. Beck, Emma Bordier, Emily Canon, Ana Escorza, Mattias Fabry, Lars Hermansson, Steve Howell, Grant Miller, Shreeya Sheyte, Safaa Alhassan, Elisabeth M. L. Baeten, Frank Barnet, Stewart. J. Bean, Mikael Bernau, David M. Bundy , et al. (15 additional authors not shown)

Abstract: We report the discovery and analysis of a massive, compact, hierarchical triple system (TIC 470710327) initially identified by citizen scientists in data obtained by NASA's Transiting Exoplanet Survey Satellite (TESS). Spectroscopic follow-up observations obtained with the HERMES spectrograph, combined with eclipse timing variations (ETVs), confirm that the system is comprised of three OB stars, w… ▽ More We report the discovery and analysis of a massive, compact, hierarchical triple system (TIC 470710327) initially identified by citizen scientists in data obtained by NASA's Transiting Exoplanet Survey Satellite (TESS). Spectroscopic follow-up observations obtained with the HERMES spectrograph, combined with eclipse timing variations (ETVs), confirm that the system is comprised of three OB stars, with a compact 1.10 d eclipsing binary and a non-eclipsing tertiary on a 52.04 d orbit. Dynamical modelling of the system (from radial velocity and ETVs) reveal a rare configuration wherein the tertiary star (O9.5-B0.5V; 14-17 M$_{\odot}$) is more massive than the combined mass of the inner binary (10.9-13.2 M$_{\odot}$). Given the high mass of the tertiary, we predict that this system will undergo multiple phases of mass transfer in the future, and likely end up as a double neutron star gravitational wave progenitor or an exotic Thorne-Zytkow object. Further observational characterisation of this system promises constraints on both formation scenarios of massive stars as well as their exotic evolutionary end-products. △ Less

Submitted 14 February, 2022; originally announced February 2022.

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

TESS Eclipsing Binary Stars. I. Short cadence observations of 4584 eclipsing binaries in Sectors 1-26

Authors: Andrej Prsa, Angela Kochoska, Kyle E. Conroy, Nora Eisner, Daniel R. Hey, Luc IJspeert, Ethan Kruse, Scott W. Fleming, Cole Johnston, Martti H. Kristiansen, Daryll LaCourse, Danielle Mortensen, Joshua Pepper, Keivan G. Stassun, Guillermo Torres, Michael Abdul-Masih, Joheen Chakraborty, Robert Gagliano, Zhao Guo, Kelly Hambleton, Kyeongsoo Hong, Thomas Jacobs, David Jones, Veselin Kostov, Jae Woo Lee , et al. (22 additional authors not shown)

Abstract: In this paper we present a catalog of 4584 eclipsing binaries observed during the first two years (26 sectors) of the TESS survey. We discuss selection criteria for eclipsing binary candidates, detection of hither-to unknown eclipsing systems, determination of the ephemerides, the validation and triage process, and the derivation of heuristic estimates for the ephemerides. Instead of keeping to th… ▽ More In this paper we present a catalog of 4584 eclipsing binaries observed during the first two years (26 sectors) of the TESS survey. We discuss selection criteria for eclipsing binary candidates, detection of hither-to unknown eclipsing systems, determination of the ephemerides, the validation and triage process, and the derivation of heuristic estimates for the ephemerides. Instead of keeping to the widely used discrete classes, we propose a binary star morphology classification based on a dimensionality reduction algorithm. Finally, we present statistical properties of the sample, we qualitatively estimate completeness, and discuss the results. The work presented here is organized and performed within the TESS Eclipsing Binary Working Group, an open group of professional and citizen scientists; we conclude by describing ongoing work and future goals for the group. The catalog is available from http://tessEBs.villanova.edu and from MAST. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 32 pages, 21 figures, accepted to ApJ Supplement Series; comments welcome

An all-sky sample of intermediate- to high-mass OBA-type eclipsing binaries observed by TESS

Authors: Luc W. IJspeert, Andrew Tkachenko, Cole Johnston, Stefano Garcia, Joris De Ridder, Timothy Van Reeth, Conny Aerts

Abstract: Context. Intermediate- to high-mass stars are the least numerous types of stars and they are less well understood than their more numerous low-mass counterparts in terms of their internal physical processes. Modelling the photometric variability of a large sample of main-sequence intermediate- to high-mass stars in eclipsing binary systems will help to improve the models for such stars. Aims. Our… ▽ More Context. Intermediate- to high-mass stars are the least numerous types of stars and they are less well understood than their more numerous low-mass counterparts in terms of their internal physical processes. Modelling the photometric variability of a large sample of main-sequence intermediate- to high-mass stars in eclipsing binary systems will help to improve the models for such stars. Aims. Our goal is to compose a homogeneously compiled sample of main-sequence intermediate- to high-mass OBA-type dwarfs in eclipsing binary systems from TESS photometry. We search for binaries with and without pulsations and determine their approximate ephemerides. Methods. Our selection starts from a catalogue of dwarfs with colours corresponding to those of OBA-type dwarfs in the TESS Input Catalog. We develop a new automated method aimed at detecting eclipsing binaries in the presence of strong pulsational and/or rotational signal relative to the eclipse depths and apply it to publicly available 30-min cadence TESS light curves. Results. Using targets with TESS magnitudes below 15 and cuts in the 2MASS magnitude bands of $J - H < 0.045$ and $J - K < 0.06$ as most stringent criteria, we arrive at a total of 189 981 intermediate- to high-mass candidates, 91193 of which have light curves from at least one of two data reduction pipelines. The eclipsing binary detection and subsequent manual check for false positives resulted in 3155 unique OBA-type eclipsing binary candidates. Conclusions. Our sample of eclipsing binary stars in the intermediate- to high-mass regime allows for future binary (and asteroseismic) modelling with the aim to better understand the internal physical processes in this hot part of the main sequence. △ Less

Submitted 21 July, 2021; originally announced July 2021.

Comments: 12 pages, 9 figures, 2 tables, 3 appendices (8 additional figures, 1 additional table). Accepted for publication in A&A

Journal ref: A&A 652, A120 (2021)

One size does not fit all: Evidence for a range of mixing efficiencies in stellar evolution calculations

Authors: Cole Johnston

Abstract: Context: Internal chemical mixing in intermediate- and high-mass stars represents an immense uncertainty in stellar evolution models.In addition to extending the main-sequence lifetime, chemical mixing also appreciably increases the mass of the stellar core. Several studies have made attempts to calibrate the efficiency of different convective boundary mixing mechanisms, with sometimes seemingly c… ▽ More Context: Internal chemical mixing in intermediate- and high-mass stars represents an immense uncertainty in stellar evolution models.In addition to extending the main-sequence lifetime, chemical mixing also appreciably increases the mass of the stellar core. Several studies have made attempts to calibrate the efficiency of different convective boundary mixing mechanisms, with sometimes seemingly conflicting results. Aims: We aim to demonstrate that stellar models regularly under-predict the masses of convective stellar cores. Methods: We gather convective core mass and fractional core hydrogen content inferences from numerous independent binary and asteroseismic studies, and compare them to stellar evolution models computed with the MESA stellar evolution code. Results: We demonstrate that core mass inferences from the literature are ubiquitously more massive than predicted by stellar evolution models without or with little convective boundary mixing. Conclusions: Independent of the form of internal mixing, stellar models require an efficient mixing mechanism that produces more massive cores throughout the main sequence to reproduce high-precision observations. This has implications for the post-main sequence evolution of all stars which have a well developed convective core on the main sequence. △ Less

Submitted 21 July, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 5 pages, 2 figures, 1 appendix table. Accepted for publication in A&A. References in table corrected

Journal ref: A&A 655, A29 (2021)

TESS Data for Asteroseismology (T'DA) Stellar Variability Classification Pipeline: Set-Up and Application to the Kepler Q9 Data

Authors: Jeroen Audenaert, James S. Kuszlewicz, Rasmus Handberg, Andrew Tkachenko, David J. Armstrong, Marc Hon, Refilwe Kgoadi, Mikkel N. Lund, Keaton J. Bell, Lisa Bugnet, Dominic M. Bowman, Cole Johnston, Rafael A. García, Dennis Stello, László Molnár, Emese Plachy, Derek Buzasi, Conny Aerts, the T'DA collaboration

Abstract: The NASA Transiting Exoplanet Survey Satellite (TESS) is observing tens of millions of stars with time spans ranging from $\sim$ 27 days to about 1 year of continuous observations. This vast amount of data contains a wealth of information for variability, exoplanet, and stellar astrophysics studies but requires a number of processing steps before it can be fully utilized. In order to efficiently p… ▽ More The NASA Transiting Exoplanet Survey Satellite (TESS) is observing tens of millions of stars with time spans ranging from $\sim$ 27 days to about 1 year of continuous observations. This vast amount of data contains a wealth of information for variability, exoplanet, and stellar astrophysics studies but requires a number of processing steps before it can be fully utilized. In order to efficiently process all the TESS data and make it available to the wider scientific community, the TESS Data for Asteroseismology working group, as part of the TESS Asteroseismic Science Consortium, has created an automated open-source processing pipeline to produce light curves corrected for systematics from the short- and long-cadence raw photometry data and to classify these according to stellar variability type. We will process all stars down to a TESS magnitude of 15. This paper is the next in a series detailing how the pipeline works. Here, we present our methodology for the automatic variability classification of TESS photometry using an ensemble of supervised learners that are combined into a metaclassifier. We successfully validate our method using a carefully constructed labelled sample of Kepler Q9 light curves with a 27.4 days time span mimicking single-sector TESS observations, on which we obtain an overall accuracy of 94.9%. We demonstrate that our methodology can successfully classify stars outside of our labeled sample by applying it to all $\sim$ 167,000 stars observed in Q9 of the Kepler space mission. △ Less

Submitted 13 July, 2021; originally announced July 2021.

Comments: 35 pages, 17 figures, 6 tables, Accepted for publication in The Astronomical Journal

A fast multi-dimensional magnetohydrodynamic formulation of the transition region adaptive conduction (TRAC) method

Authors: C. D. Johnston, A. W. Hood, I. De Moortel, P. Pagano, T. A. Howson

Abstract: We have demonstrated that the Transition Region Adaptive Conduction (TRAC) method permits fast and accurate numerical solutions of the field-aligned hydrodynamic equations, successfully removing the influence of numerical resolution on the coronal density response to impulsive heating. This is achieved by adjusting the parallel thermal conductivity, radiative loss, and heating rates to broaden the… ▽ More We have demonstrated that the Transition Region Adaptive Conduction (TRAC) method permits fast and accurate numerical solutions of the field-aligned hydrodynamic equations, successfully removing the influence of numerical resolution on the coronal density response to impulsive heating. This is achieved by adjusting the parallel thermal conductivity, radiative loss, and heating rates to broaden the transition region (TR), below a global cutoff temperature, so that the steep gradients are spatially resolved even when using coarse numerical grids. Implementing the original 1D formulation of TRAC in multi-dimensional magnetohydrodynamic (MHD) models would require tracing a large number of magnetic field lines at every time step in order to prescribe a global cutoff temperature to each field line. In this paper, we present a highly efficient formulation of the TRAC method for use in multi-dimensional MHD simulations, which does not rely on field line tracing. In the TR, adaptive local cutoff temperatures are used instead of global cutoff temperatures to broaden any unresolved parts of the atmosphere. These local cutoff temperatures are calculated using only local grid cell quantities, enabling the MHD extension of TRAC to efficiently account for the magnetic field evolution, without tracing field lines. Consistent with analytical predictions, we show that this approach successfully preserves the properties of the original TRAC method. In particular, the total radiative losses and heating remain conserved under the MHD formulation. Results from 2D MHD simulations of impulsive heating in unsheared and sheared arcades of coronal loops are also presented. These simulations benchmark the MHD TRAC method against a series of 1D models and demonstrate the versatility and robustness of the method in multi-dimensional magnetic fields. We show, for the first time, that pressure differences, ... △ Less

Submitted 21 June, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Comments: 18 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 654, A2 (2021)

Multi-Stranded Coronal Loops: Quantifying Strand Number and Heating Frequency from Simulated Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA) Observations

Authors: Thomas Williams, Robert W. Walsh, Stephane Regnier, Craig D. Johnston

Abstract: Coronal loops form the basic building blocks of the magnetically closed solar corona yet much is still to be determined concerning their possible fine-scale structuring and the rate of heat deposition within them. Using an improved multi-stranded loop model to better approximate the numerically challenging transition region, this paper examines synthetic NASA Solar Dynamics Observatory's (SDO) Atm… ▽ More Coronal loops form the basic building blocks of the magnetically closed solar corona yet much is still to be determined concerning their possible fine-scale structuring and the rate of heat deposition within them. Using an improved multi-stranded loop model to better approximate the numerically challenging transition region, this paper examines synthetic NASA Solar Dynamics Observatory's (SDO) Atmospheric Imaging Assembly (AIA) emission simulated in response to a series of prescribed spatially and temporally random, impulsive and localised heating events across numerous sub-loop elements with a strong weighting towards the base of the structure; the nanoflare heating scenario. The total number of strands and nanoflare repetition times are varied systematically in such a way that the total energy content remains approximately constant across all the cases analysed. Repeated time lag detection during an emission time series provides a good approximation for the nanoflare repetition time for low-frequency heating. Furthermore, using a combination of AIA 171/193 and 193/211 channel ratios in combination with spectroscopic determination of the standard deviation of the loop apex temperature over several hours alongside simulations from the outlined multi-stranded loop model, it is demonstrated that both the imposed heating rate and number of strands can be realised. △ Less

Submitted 26 May, 2021; originally announced May 2021.

A comparison of the dynamical and model-derived parameters of the pulsating eclipsing binary KIC9850387

Authors: S. Sekaran, A. Tkachenko, C. Johnston, C. Aerts

Abstract: We aim to unravel the interior mixing profile of the pulsating eclipsing binary KIC9850387 by comparing its dynamical parameters and the parameters derived through a combination of evolutionary and asteroseismic modelling. We created a grid of stellar evolutionary models using the stellar evolutionary code \textsc{mesa} and performed an isochrone-cloud (isocloud) based evolutionary modelling of th… ▽ More We aim to unravel the interior mixing profile of the pulsating eclipsing binary KIC9850387 by comparing its dynamical parameters and the parameters derived through a combination of evolutionary and asteroseismic modelling. We created a grid of stellar evolutionary models using the stellar evolutionary code \textsc{mesa} and performed an isochrone-cloud (isocloud) based evolutionary modelling of the system. We then generated a grid of pulsational models using the stellar pulsation code \textsc{gyre} based on the age constraints from the evolutionary modelling. Finally, we performed asteroseismic modelling of the observed $\ell=1$ and $\ell=2$ period-spacing patterns, utilising different combinations of observational constraints, merit functions, and asteroseismic observables to obtain strong constraints on the interior properties of the primary star. Through a combination of asteroseismic modelling and dynamical constraints, we found that the system comprises two main-sequence components at an age of $1.2\pm0.1$ Gyr. We also observed that asteroseismic modelling provided stronger constraints on the interior properties than evolutionary modelling. Overall, we found high levels of interior mixing and posited that this is a result of intrinsic non-tidal mixing mechanisms due to a similar observed behaviour in single stars. We investigated the high-frequency regime of KIC9850387 and found evidence of the surface effect as well as rotational splitting in the form of a prograde-retrograde dipole g$_{1}$ mode doublet with a missing zonal mode. We find that the dynamical parameters and the parameters extracted from the asteroseismic modelling of period-spacing patterns are only barely compliant, reinforcing the need for homogeneous analyses of samples of pulsating eclipsing binaries that aim to calibrating interior mixing profiles. △ Less

Submitted 22 February, 2021; originally announced February 2021.

Comments: 21 pages, 18 figures, accepted for publication in A&A

Journal ref: A&A 648, A91 (2021)

Characterization of the variability in the O+B eclipsing binary HD 165246

Authors: Cole Johnston, Nicolas Aimar, Michael Abdul-Masih, Dominic. M. Bowman, Tim R. White, Calum Hawcroft, Hugues Sana, Sanjay Sekaran, Karan Dsilva, Andrew Tkachenko, Conny Aerts

Abstract: O-stars are known to experience a wide range of variability mechanisms originating at both their surface and their near-core regions. Characterization and understanding of this variability and its potential causes are integral for evolutionary calculations. We use a new extensive high-resolution spectroscopic data set to characterize the variability observed in both the spectroscopic and space-bas… ▽ More O-stars are known to experience a wide range of variability mechanisms originating at both their surface and their near-core regions. Characterization and understanding of this variability and its potential causes are integral for evolutionary calculations. We use a new extensive high-resolution spectroscopic data set to characterize the variability observed in both the spectroscopic and space-based photometric observations of the O+B eclipsing binary HD~165246. We present an updated atmospheric and binary solution for the primary component, involving a high level of microturbulence ($13_{-1.3}^{+1.0}\,$km\,s$^{-1}$) and a mass of $M_1=23.7_{-1.4}^{+1.1}$~M$_{\odot}$, placing it in a sparsely explored region of the Hertzsprung-Russell diagram. Furthermore, we deduce a rotational frequency of $0.690\pm 0.003\,$d$^{-1}$ from the combined photometric and line-profile variability, implying that the primary rotates at 40\% of its critical Keplerian rotation rate. We discuss the potential explanations for the overall variability observed in this massive binary, and discuss its evolutionary context. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society, 14 pages, 10 figures, 6 tables, 1 appendix (4 additional figures)

Discovery and Characterization of a Rare Magnetic Hybrid $β$ Cephei Slowly Pulsating B-type Star in an Eclipsing Binary in the Young Open Cluster NGC 6193

Authors: Keivan G. Stassun, Guillermo Torres, Cole Johnston, Daniel J. Stevens, Dax L. Feliz, Marina Kounkel, Luke G. Bouma

Abstract: As many as 10\% of OB-type stars have global magnetic fields, which is surprising given their internal structure is radiative near the surface. A direct probe of internal structure is pulsations, and some OB-type stars exhibit pressure modes ($β$ Cep pulsators) or gravity modes (slowly pulsating B-type stars; SPBs); a few rare cases of hybrid $β$ Cep/SPBs occupy a narrow instability strip in the H… ▽ More As many as 10\% of OB-type stars have global magnetic fields, which is surprising given their internal structure is radiative near the surface. A direct probe of internal structure is pulsations, and some OB-type stars exhibit pressure modes ($β$ Cep pulsators) or gravity modes (slowly pulsating B-type stars; SPBs); a few rare cases of hybrid $β$ Cep/SPBs occupy a narrow instability strip in the H-R diagram. The most precise fundamental properties of stars are obtained from eclipsing binaries (EBs), and those in clusters with known ages and metallicities provide the most stringent constraints on theory. Here we report the discovery that HD 149834 in the $\sim$5 Myr cluster NGC 6193 is an EB comprising a hybrid $β$ Cep/SPB pulsator and a highly irradiated low-mass companion. We determine the masses, radii, and temperatures of both stars; the $\sim$9.7 M$_\odot$ primary resides in the instability strip where hybrid pulsations are theoretically predicted. The presence of both SPB and $β$ Cep pulsations indicates that the system has a near-solar metallicity, and is in the second half of the main-sequence lifetime. The radius of the $\sim$1.2 M$_\odot$ companion is consistent with theoretical pre-main-sequence isochrones at 5 Myr, but its temperature is much higher than expected, perhaps due to irradiation by the primary. The radius of the primary is larger than expected, unless its metallicity is super-solar. Finally, the light curve shows residual modulation consistent with the rotation of the primary, and Chandra observations reveal a flare, both of which suggest the presence of starspots and thus magnetism on the primary. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: Accepted for publication in The Astrophysical Journal, 15 pages, 14 figures

Planet Hunters TESS II: Findings from the first two years of TESS

Authors: Nora L. Eisner, Oscar Barragán, Chris Lintott, Suzanne Aigrain, Belinda Nicholson, Tabetha S. Boyajian, Steve B. Howell, Cole Johnston, Ben Lakeland, Grant Miller, Adam McMaster, Hannu Parviainen, Emily J. Safron, Megan E. Schwamb, Laura Trouille, Sophia Vaughan, Norbert Zicher, Campbell Allen, Sarah Allen, Mark Bouslog, Cliff Johnson, Molly N. Simon, Zach Wolfenbarger, Elisabeth M. L. Baeten, David M. Bundy , et al. (1 additional authors not shown)

Abstract: We present the results from the first two years of the Planet Hunters TESS citizen science project, which identifies planet candidates in the TESS data by engaging members of the general public. Over 22,000 citizen scientists from around the world visually inspected the first 26 Sectors of TESS data in order to help identify transit-like signals. We use a clustering algorithm to combine these clas… ▽ More We present the results from the first two years of the Planet Hunters TESS citizen science project, which identifies planet candidates in the TESS data by engaging members of the general public. Over 22,000 citizen scientists from around the world visually inspected the first 26 Sectors of TESS data in order to help identify transit-like signals. We use a clustering algorithm to combine these classifications into a ranked list of events for each sector, the top 500 of which are then visually vetted by the science team. We assess the detection efficiency of this methodology by comparing our results to the list of TESS Objects of Interest (TOIs) and show that we recover 85 % of the TOIs with radii greater than 4 Earth radii and 51 % of those with radii between 3 and 4 Earth radii. Additionally, we present our 90 most promising planet candidates that had not previously been identified by other teams, 73 of which exhibit only a single transit event in the TESS light curve, and outline our efforts to follow these candidates up using ground-based observatories. Finally, we present noteworthy stellar systems that were identified through the Planet Hunters TESS project. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: Accepted for publication in MNRAS (22 pages, 12 figures, 3 tables)

V772 Cas: an ellipsoidal HgMn star in an eclipsing binary

Authors: O. Kochukhov, C. Johnston, J. Labadie-Bartz, S. Shetye, T. A. Ryabchikova, A. Tkachenko, M. E. Shultz

Abstract: The late B-type star V772 Cas (HD 10260) was previously suspected to be a rare example of a magnetic chemically peculiar star in an eclipsing binary system. Photometric observations of this star obtained by the TESS satellite show clear eclipses with a period of 5.0137 d accompanied by a significant out-of-eclipse variation with the same period. High-resolution spectroscopy reveals V772 Cas to be… ▽ More The late B-type star V772 Cas (HD 10260) was previously suspected to be a rare example of a magnetic chemically peculiar star in an eclipsing binary system. Photometric observations of this star obtained by the TESS satellite show clear eclipses with a period of 5.0137 d accompanied by a significant out-of-eclipse variation with the same period. High-resolution spectroscopy reveals V772 Cas to be an SB1 system, with the primary component rotating about a factor two slower than the orbital period and showing chemical peculiarities typical of non-magnetic HgMn chemically peculiar stars. This is only the third eclipsing HgMn star known and, owing to its brightness, is one of the very few eclipsing binaries with chemically peculiar components accessible to detailed follow-up studies. Taking advantage of the photometric and spectroscopic observations available for V772 Cas, we performed modelling of this system with the PHOEBE code. This analysis provided fundamental parameters of the components and demonstrated that the out-of-eclipse brightness variation is explained by the ellipsoidal shape of the evolved, asynchronously rotating primary. This is the first HgMn star for which such variability has been definitively identified. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: 14 pages, 13 figures; accepted for publication in MNRAS

Tango of celestial dancers: A sample of detached eclipsing binary systems containing g-mode pulsating components. A case study of KIC9850387

Authors: S. Sekaran, A. Tkachenko, M. Abdul-Masih, A. Prsa, C. Johnston, D. Huber, S. J. Murphy, G. Banyard, A. W. Howard, H. Isaacson, D. M. Bowman, C. Aerts

Abstract: Eclipsing binary systems with components that pulsate in gravity modes ($g$ modes) allow for simultaneous and independent constraints of the chemical mixing profiles of stars. We aim to assemble a sample of $g$-mode pulsators in detached eclipsing binaries with the purpose of finding good candidates for future evolutionary and asteroseismic modelling. In addition, we present a case study of the ec… ▽ More Eclipsing binary systems with components that pulsate in gravity modes ($g$ modes) allow for simultaneous and independent constraints of the chemical mixing profiles of stars. We aim to assemble a sample of $g$-mode pulsators in detached eclipsing binaries with the purpose of finding good candidates for future evolutionary and asteroseismic modelling. In addition, we present a case study of the eclipsing binary KIC9850387, identified as our most promising candidate. We selected all of the detached eclipsing binaries in the Kepler eclipsing binary catalogue and performed a visual inspection to determine the presence and density of $g$ modes, and the presence of $g$-mode period-spacing patterns in their frequency spectra. We then characterised our sample based on their $g$-mode pulsational parameters and binary and atmospheric parameters. A spectroscopic follow-up of KIC9850387 was then performed. We then performed spectral disentangling followed by atmospheric modelling and abundance analysis for the primary star. We utilised an iterative approach to simultaneously optimise the pulsational and eclipse models, and subsequently performed an analysis of the pressure- ($p$-) and $g$-mode pulsational frequencies. We compiled a sample of 93 Kepler eclipsing binary stars with $g$-mode pulsating components and identified clear $g$-mode period-spacing patterns in the frequency spectra of seven of these systems. We found that the $g$-mode pulsational parameters and the binary and atmospheric parameters of our sample are weakly correlated at best. We find that the eclipsing binary KIC9850387 is a double-lined spectroscopic binary in a near-circular orbit with an intermediate-mass primary, and a solar-like secondary. We find $\ell=1$ and $\ell=2$ period-spacing patterns in the frequency spectrum of KIC9850387 spanning more than ten radial orders each. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 27 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 643, A162 (2020)

The TESS light curve of AI Phoenicis

Authors: P. F. L. Maxted, Patrick Gaulme, D. Graczyk, K. G. Hełminiak, C. Johnston, Jerome A. Orosz, Andrej Prša, John Southworth, Guillermo Torres, Guy R. Davies, Warrick Ball, William J Chaplin, .

Abstract: Accurate masses and radii for normal stars derived from observations of detached eclipsing binary stars are of fundamental importance for testing stellar models and may be useful for calibrating free parameters in these model if the masses and radii are sufficiently precise and accurate. We aim to measure precise masses and radii for the stars in the bright eclipsing binary AI Phe, and to quantify… ▽ More Accurate masses and radii for normal stars derived from observations of detached eclipsing binary stars are of fundamental importance for testing stellar models and may be useful for calibrating free parameters in these model if the masses and radii are sufficiently precise and accurate. We aim to measure precise masses and radii for the stars in the bright eclipsing binary AI Phe, and to quantify the level of systematic error in these estimates. We use several different methods to model the TESS light curve of AI Phe combined with spectroscopic orbits from multiple sources to estimate precisely the stellar masses and radii together with robust error estimates. We find that the agreement between different methods for the light curve analysis is very good but some methods underestimate the errors on the model parameters. The semi-amplitudes of the spectroscopic orbits derived from spectra obtained with modern echelle spectrographs are consistent to within 0.1%. The masses of the stars in AI Phe are $M_1 = 1.1938 \pm 0.0008 M_{\odot}$ and $M_2 = 1.2438 \pm 0.0008M_{\odot}$, and the radii are $R_1 = 1.8050 \pm 0.0022 R_{\odot}$ and $R_2 = 2.9332 \pm 0.0023 R_{\odot}$. We conclude that it is possible to measure accurate masses and radii for stars in bright eclipsing binary stars to a precision of 0.2% or better using photometry from TESS and spectroscopy obtained with modern echelle spectrographs. We provide recommendations for publishing masses and radii of eclipsing binary stars at this level of precision. △ Less

Submitted 8 June, 2020; v1 submitted 20 March, 2020; originally announced March 2020.

Comments: Accepted for publication in MNRAS. 12 pages, 4 tables and 8 figures

The mass discrepancy in intermediate- and high-mass eclipsing binaries: the need for higher convective core masses

Authors: A. Tkachenko, K. Pavlovski, C. Johnston, M. G. Pedersen, M. Michielsen, D. M. Bowman, J. Southworth, V. Tsymbal, C. Aerts

Abstract: (Abridged) Eclipsing, spectroscopic double-lined binary star systems (SB2) are excellent laboratories for calibrating theories of stellar interior structure and evolution. We aim to investigate the mass discrepancy in binary stars. We study the effect of near-core mixing on the mass of the convective core of the stars and interpret the results in the context of the mass discrepancy. Two scenarios… ▽ More (Abridged) Eclipsing, spectroscopic double-lined binary star systems (SB2) are excellent laboratories for calibrating theories of stellar interior structure and evolution. We aim to investigate the mass discrepancy in binary stars. We study the effect of near-core mixing on the mass of the convective core of the stars and interpret the results in the context of the mass discrepancy. Two scenarios are considered, where individual stellar components of a binary system are treated independent of each other and where they are forced to have the same age and initial chemical composition. We find that the mass discrepancy is present in our sample and that it is anti-correlated with the surface gravity of the star. No correlations are found with other fundamental and atmospheric parameters, including the stellar mass. The mass discrepancy can be partially accounted for by increasing the amount of near-core mixing in stellar evolution models. We also find that ignoring the microturbulent velocity and turbulent pressure in stellar atmosphere models of hot evolved stars results in overestimation of their effective temperature by up to 8%. Together with enhanced near-core mixing, this can almost entirely account for the 30% mass discrepancy found for the evolved primary component of V380 Cyg. We find a strong link between the mass discrepancy and the convective core mass. The mass discrepancy can be solved by considering the combined effect of extra near-core boundary mixing and consistent treatment in the spectrum analysis of hot evolved stars. Our binary modelling results in convective core masses between 17 and 35% of the stellar mass, in excellent agreement with results from gravity-mode asteroseismology of single stars. This implies larger helium core masses near the end of the main sequence than anticipated so far. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: 22 pages (17 pages of main text + 5 pages of appendix), 14 figures. Accepted for publication in A&A

Chromospheric evaporation and phase mixing of Alfvén waves in coronal loops

Authors: H. J. Van Damme, I. De Moortel, P. Pagano, C. D. Johnston

Abstract: Phase mixing of Alfvén waves has been studied extensively as a possible coronal heating mechanism but without the full thermodynamic consequences considered self-consistently. It has been argued that in some cases, the thermodynamic feedback of the heating could substantially affect the transverse density gradient and even inhibit the phase mixing process. In this paper, we use MHD simulations wit… ▽ More Phase mixing of Alfvén waves has been studied extensively as a possible coronal heating mechanism but without the full thermodynamic consequences considered self-consistently. It has been argued that in some cases, the thermodynamic feedback of the heating could substantially affect the transverse density gradient and even inhibit the phase mixing process. In this paper, we use MHD simulations with the appropriate thermodynamical terms included to quantify the evaporation following heating by phase mixing of Alfvén waves in a coronal loop and the effect of this evaporation on the transverse density profile. The numerical simulations were performed using the Lare2D code. We set up a 2D loop model consisting of a field-aligned thermodynamic equilibrium and a cross-field (background) heating profile. A continuous, sinusoidal, high-frequency Alfvén wave driver was implemented. As the Alfvén waves propagate along the field, they undergo phase mixing due to the cross-field density gradient in the coronal part of the loop. We investigated the presence of field-aligned flows, heating from the dissipation of the phase-mixed Alfvén waves, and the subsequent evaporation from the lower atmosphere. We find that phase mixing of Alfvén waves leads to modest heating in the shell regions of the loop and evaporation of chromospheric material into the corona with upflows of the order of only 5-20 m/s. Although the evaporation leads to a mass increase in the shell regions of the loop, the effect on the density gradient and, hence, on the phase mixing process, is insignificant. This paper self-consistently investigates the effect of chromospheric evaporation on the cross-field density gradient and the phase mixing process in a coronal loop. We found that the effects in our particular setup (small amplitude, high frequency waves) are too small to significantly change the density gradient. △ Less

Submitted 26 February, 2020; originally announced February 2020.

Comments: 11 pages, 12 figures, accepted for publication in A&A

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

Modelling the solar transition region using an adaptive conduction method

Authors: C. D. Johnston, P. J. Cargill, A. W. Hood, I. De Moortel, S. J. Bradshaw, A. C. Vaseekar

Abstract: Modelling the solar Transition Region with the use of an Adaptive Conduction (TRAC) method permits fast and accurate numerical solutions of the field-aligned hydrodynamic equations, capturing the enthalpy exchange between the corona and transition region, when the corona undergoes impulsive heating. The TRAC method eliminates the need for highly resolved numerical grids in the transition region an… ▽ More Modelling the solar Transition Region with the use of an Adaptive Conduction (TRAC) method permits fast and accurate numerical solutions of the field-aligned hydrodynamic equations, capturing the enthalpy exchange between the corona and transition region, when the corona undergoes impulsive heating. The TRAC method eliminates the need for highly resolved numerical grids in the transition region and the commensurate very short time steps that are required for numerical stability. When employed with coarse spatial resolutions, typically achieved in multi-dimensional magnetohydrodynamic codes, the errors at peak density are less than 5% and the computation time is three orders of magnitude faster than fully resolved field-aligned models. This paper presents further examples that demonstrate the versatility and robustness of the method over a range of heating events, including impulsive and quasi-steady footpoint heating. A detailed analytical assessment of the TRAC method is also presented, showing that the approach works through all phases of an impulsive heating event because (i) the total radiative losses and (ii) the total heating when integrated over the transition region are both preserved at all temperatures under the broadening modifications of the method. The results from the numerical simulations complement this conclusion. △ Less

Submitted 28 February, 2020; v1 submitted 5 February, 2020; originally announced February 2020.

Comments: 20 pages, 9 figures, 1 table, accepted for publication in A&A

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

MOBSTER -- III. HD 62658: a magnetic Bp star in an eclipsing binary with a non-magnetic 'identical twin'

Authors: M. E. Shultz, C. Johnston, J. Labadie-Bartz, V. Petit, A. David-Uraz, O. Kochukhov, G. A. Wade, J. Pepper, K. G. Stassun, J. E. Rodriguez, M. B. Lund, D. J. James

Abstract: HD 62658 (B9p V) is a little-studied chemically peculiar star. Light curves obtained by the Kilodegree Extremely Little Telescope (KELT) and Transiting Exoplanet Survey Satellite (TESS) show clear eclipses with a period of about 4.75 d, as well as out-of-eclipse brightness modulation with the same 4.75 d period, consistent with synchronized rotational modulation of surface chemical spots. High-res… ▽ More HD 62658 (B9p V) is a little-studied chemically peculiar star. Light curves obtained by the Kilodegree Extremely Little Telescope (KELT) and Transiting Exoplanet Survey Satellite (TESS) show clear eclipses with a period of about 4.75 d, as well as out-of-eclipse brightness modulation with the same 4.75 d period, consistent with synchronized rotational modulation of surface chemical spots. High-resolution ESPaDOnS circular spectropolarimetry shows a clear Zeeman signature in the line profile of the primary; there is no indication of a magnetic field in the secondary. PHOEBE modelling of the light curve and radial velocities indicates that the two components have almost identical masses of about 3 M$_\odot$. The primary's longitudinal magnetic field $\langle B_z \rangle$ varies between about $+100$ and $-250$ G, suggesting a surface magnetic dipole strength $B_{\rm d} = 850$~G. Bayesian analysis of the Stokes $V$ profiles indicates $B_{\rm d} = 650$~G for the primary and $B_{\rm d} < 110$ G for the secondary. The primary's line profiles are highly variable, consistent with the hypothesis that the out-of-eclipse brightness modulation is a consequence of rotational modulation of that star's chemical spots. We also detect a residual signal in the light curve after removal of the orbital and rotational modulations, which might be pulsational in origin; this could be consistent with the weak line profile variability of the secondary. This system represents an excellent opportunity to examine the consequences of magnetic fields for stellar structure via comparison of two stars that are essentially identical with the exception that one is magnetic. The existence of such a system furthermore suggests that purely environmental explanations for the origin of fossil magnetic fields are incomplete. △ Less

Submitted 7 October, 2019; originally announced October 2019.

Comments: 14 pages, 10 figures, accepted for publication in MNRAS

Isochrone-cloud fitting of the extended Main-Sequence Turn-Off of young clusters

Authors: C. Johnston, C. Aerts, M. G. Pedersen, N. Bastian

Abstract: Extended main-sequence turn-offs (eMSTO) are a commonly observed property of young clusters. A global theoretical interpretation for the eMSTOs is still lacking, but stellar rotation is considered a necessary ingredient to explain the eMSTO. We aim to assess the importance of core-boundary and envelope mixing in stellar interiors for the interpretation of eMSTOs in terms of one coeval population.… ▽ More Extended main-sequence turn-offs (eMSTO) are a commonly observed property of young clusters. A global theoretical interpretation for the eMSTOs is still lacking, but stellar rotation is considered a necessary ingredient to explain the eMSTO. We aim to assess the importance of core-boundary and envelope mixing in stellar interiors for the interpretation of eMSTOs in terms of one coeval population. We construct isochrone-clouds based on interior mixing profiles of stars with a convective core calibrated from asteroseismology of isolated galactic field stars. We fit these isochrone-clouds to the measured eMSTO to estimate the age and core mass of the stars in the two young clusters NGC 1850 and NGC 884, assuming one coeval population and fixing the metallicity to the one measured from spectroscopy. We assess the correlations between the interior mixing properties of the cluster members and their rotational and pulsation properties. We find that stellar models based on asteroseismically-calibrated interior mixing profiles lead to enhanced core masses of eMSTO stars and can explain a good fraction of the observed eMSTOs of the two considered clusters in terms of one coeval population of stars, with similar ages to those in the literature, given the large uncertainties. The rotational and pulsation properties of the stars in NGC 884 are not sufficiently well known to perform asteroseismic modelling, as it is achieved for field stars from space photometry. The stars in NGC 884 for which we have vsini and a few pulsation frequencies show no correlation between these properties and the core masses of the stars that set the cluster age. Future cluster space asteroseismology may allow to interpret the values of the core masses in terms of the physical processes that cause them, based on the modelling of the interior mixing profiles for the individual member stars with suitable identified modes. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: 11 pages, 8 figures, 1 table, NGC 1850, NGC 884, accepted for publication in A&A

Journal ref: A&A 632, A74 (2019)

The first view of $δ$ Scuti and $γ$ Doradus stars with the TESS mission

Authors: V. Antoci, M. S. Cunha, D. M. Bowman, S. J. Murphy, D. W. Kurtz, T. R. Bedding, C. C. Borre, S. Christophe, J. Daszyńska-Daszkiewicz, L. Fox-Machado, A. García Hernández, H. Ghasemi, R. Handberg, H. Hansen, A. Hasanzadeh, G. Houdek, C. Johnston, A. B. Justesen, F. Kahraman Alicavus, K. Kotysz, D. Latham, J. M. Matthews, J. Mønster, E. Niemczura, E. Paunzen , et al. (41 additional authors not shown)

Abstract: We present the first asteroseismic results for $δ$ Scuti and $γ$ Doradus stars observed in Sectors 1 and 2 of the TESS mission. We utilise the 2-min cadence TESS data for a sample of 117 stars to classify their behaviour regarding variability and place them in the Hertzsprung-Russell diagram using Gaia DR2 data. Included within our sample are the eponymous members of two pulsator classes, $γ$ Dora… ▽ More We present the first asteroseismic results for $δ$ Scuti and $γ$ Doradus stars observed in Sectors 1 and 2 of the TESS mission. We utilise the 2-min cadence TESS data for a sample of 117 stars to classify their behaviour regarding variability and place them in the Hertzsprung-Russell diagram using Gaia DR2 data. Included within our sample are the eponymous members of two pulsator classes, $γ$ Doradus and SX Phoenicis. Our sample of pulsating intermediate-mass stars observed by TESS also allows us to confront theoretical models of pulsation driving in the classical instability strip for the first time and show that mixing processes in the outer envelope play an important role. We derive an empirical estimate of 74% for the relative amplitude suppression factor as a result of the redder TESS passband compared to the Kepler mission using a pulsating eclipsing binary system. Furthermore, our sample contains many high-frequency pulsators, allowing us to probe the frequency variability of hot young $δ$ Scuti stars, which were lacking in the Kepler mission data set, and identify promising targets for future asteroseismic modelling. The TESS data also allow us to refine the stellar parameters of SX Phoenicis, which is believed to be a blue straggler. △ Less

Submitted 26 September, 2019; originally announced September 2019.

Comments: 33 pages (incl. appendix), accepted for publications in MNRAS

Discovery of tidally-perturbed pulsations in the eclipsing binary U Gru: a pioneering system for tidal asteroseismology

Authors: Dominic M. Bowman, Cole Johnston, Andrew Tkachenko, David E. Mkrtichian, Khemsinan Gunsriwiwat, Conny Aerts

Abstract: The interior physics of stars is currently not well constrained for early-type stars. This is particularly pertinent for multiple systems as binary interaction becomes more prevalent for more massive stars, which strongly affects their evolution. High-precision photometry from the Transiting Exoplanet Survey Satellite (TESS) mission offers the opportunity to remedy the dearth of observations of pu… ▽ More The interior physics of stars is currently not well constrained for early-type stars. This is particularly pertinent for multiple systems as binary interaction becomes more prevalent for more massive stars, which strongly affects their evolution. High-precision photometry from the Transiting Exoplanet Survey Satellite (TESS) mission offers the opportunity to remedy the dearth of observations of pulsating stars that show evidence of binary interaction, specifically pulsating mass-accreting components of semi-detached Algol-type eclipsing binary (oEA) systems. We present the TESS light curve of the circular eclipsing binary system U Gru (TIC 147201138), which shows evidence for free heat-driven pressure modes and a series of tidally-perturbed pressure modes. We highlight the asteroseismic potential of studying pulsating stars in binary systems, and demonstrate how tidal asteroseismology can be applied to infer the influence of binary interaction on stellar structure. △ Less

Submitted 22 August, 2019; originally announced August 2019.

Comments: 9 pages, 4 figures, accepted (in press) for publication in ApJ Letters

Modelling of the B-type binaries CW Cep and U Oph: A critical view on dynamical masses, core boundary mixing, and core mass

Authors: C. Johnston, K. Pavlovski, A. Tkachenko

Abstract: Context: Intermediate-Mass stars are often overlooked as they are not supernova progenitors but still host convective cores and complex atmospheres which require computationally expensive treatment. Due to this, there is a general lack of such stars modelled by state of the art stellar structure and evolution codes. Aims: This paper aims to use high-quality spectroscopy to update the dynamically o… ▽ More Context: Intermediate-Mass stars are often overlooked as they are not supernova progenitors but still host convective cores and complex atmospheres which require computationally expensive treatment. Due to this, there is a general lack of such stars modelled by state of the art stellar structure and evolution codes. Aims: This paper aims to use high-quality spectroscopy to update the dynamically obtained stellar parameters and produce a new evolutionary assessment of the bright B0.5+B0.5 and B5V+B5V binary systems CW Cep and U Oph. Methods: We use new spectroscopy obtained with the Hermes spectrograph to revisit the photometric binary solution of the two systems. The updated mass ratio and effective temperatures are incorporated to obtain new dynamical masses for the primary and secondary. With these, we perform isochrone-cloud based evolutionary modelling to investigate the core properties of these stars. Results: We report the first abundances for CW Cep and U Oph as well as report an updated dynamical solution for both systems. We find that we cannot uniquely constrain the amount of core boundary mixing in any of the stars we consider. Instead, we report their core masses and compare our results to previous studies. Conclusions: We find that the per-cent level precision on fundamental stellar quantities are accompanied with core mass estimates to between ~ 5-15%. We find that differences in analysis techniques can lead to substantially different evolutionary mode △ Less

Submitted 28 May, 2019; originally announced May 2019.

Comments: 15 pages, 7 figures, two appendices with 4 figures each. Accepted for publication in Astronomy & Astrophysics

Low-frequency gravity waves in blue supergiants revealed by high-precision space photometry

Authors: Dominic M. Bowman, Siemen Burssens, May G. Pedersen, Cole Johnston, Conny Aerts, Bram Buysschaert, Mathias Michielsen, Andrew Tkachenko, Tamara M. Rogers, Philipp V. F. Edelmann, Rathish P. Ratnasingam, Sergio Simón-Díaz, Norberto Castro, Ehsan Moravveji, Benjamin J. S. Pope, Timothy R. White, Peter De Cat

Abstract: Almost all massive stars explode as supernovae and form a black hole or neutron star. The remnant mass and the impact of the chemical yield on subsequent star formation and galactic evolution strongly depend on the internal physics of the progenitor star, which is currently not well understood. The theoretical uncertainties of stellar interiors accumulate with stellar age, which is particularly pe… ▽ More Almost all massive stars explode as supernovae and form a black hole or neutron star. The remnant mass and the impact of the chemical yield on subsequent star formation and galactic evolution strongly depend on the internal physics of the progenitor star, which is currently not well understood. The theoretical uncertainties of stellar interiors accumulate with stellar age, which is particularly pertinent for the blue supergiant phase. Stellar oscillations represent a unique method of probing stellar interiors, yet inference for blue supergiants is hampered by a dearth of observed pulsation modes. Here we report the detection of diverse variability in blue supergiants using the K2 and TESS space missions. The discovery of pulsation modes or an entire spectrum of low-frequency gravity waves in these stars allow us to map the evolution of hot massive stars towards the ends of their lives. Future asteroseismic modelling will provide constraints on ages, core masses, interior mixing, rotation and angular momentum transport. The discovery of variability in blue supergiants is a step towards a data-driven empirical calibration of theoretical evolution models for the most massive stars in the Universe. △ Less

Submitted 6 May, 2019; originally announced May 2019.

Comments: Publication date: 6 May 2019 at 16:00h (GMT)

Journal ref: Nature Astronomy, 2019

The effects of numerical resolution, heating timescales and background heating on thermal non-equilibrium in coronal loops

Authors: C. D. Johnston, P. J. Cargill, P. Antolin, A. W. Hood, I. De Moortel, S. J. Bradshaw

Abstract: Thermal non-equilibrium (TNE) is believed to be a potentially important process in understanding some properties of the magnetically closed solar corona. Through one-dimensional hydrodynamic models, this paper addresses the importance of the numerical spatial resolution, footpoint heating timescales and background heating on TNE. Inadequate transition region (TR) resolution can lead to significant… ▽ More Thermal non-equilibrium (TNE) is believed to be a potentially important process in understanding some properties of the magnetically closed solar corona. Through one-dimensional hydrodynamic models, this paper addresses the importance of the numerical spatial resolution, footpoint heating timescales and background heating on TNE. Inadequate transition region (TR) resolution can lead to significant discrepancies in TNE cycle behaviour, with TNE being suppressed in under-resolved loops. A convergence on the periodicity and plasma properties associated with TNE required spatial resolutions of less than 2 km for a loop of length 180 Mm. These numerical problems can be resolved using an approximate method that models the TR as a discontinuity using a jump condition, as proposed by Johnston et al. (2017a,b). The resolution requirements (and so computational cost) are greatly reduced while retaining good agreement with fully resolved results. Using this approximate method we (i) identify different regimes for the response of coronal loops to time-dependent footpoint heating including one where TNE does not arise and (ii) demonstrate that TNE in a loop with footpoint heating is suppressed unless the background heating is sufficiently small. The implications for the generality of TNE are discussed. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: 15 pages, 10 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 625, A149 (2019)

Barium & related stars and their white-dwarf companions II. Main-sequence and subgiant stars

Authors: A. Escorza, D. Karinkuzhi, A. Jorissen, L. Siess, H. Van Winckel, D. Pourbaix, C. Johnston, B. Miszalski, G-M. Oomen, M. Abdul-Masih, H. M. J. Boffin, P. North, R. Manick, S. Shetye, J. Mikołajewska

Abstract: Barium (Ba) dwarfs and CH subgiants are the less-evolved analogues of Ba and CH giants. They are F- to G-type main-sequence stars polluted with heavy elements by a binary companion when the latter was on the Asymptotic Giant Branch (AGB). This companion is now a white dwarf that in most cases cannot be directly detected. We present a large systematic study of 60 objects classified as Ba dwarfs or… ▽ More Barium (Ba) dwarfs and CH subgiants are the less-evolved analogues of Ba and CH giants. They are F- to G-type main-sequence stars polluted with heavy elements by a binary companion when the latter was on the Asymptotic Giant Branch (AGB). This companion is now a white dwarf that in most cases cannot be directly detected. We present a large systematic study of 60 objects classified as Ba dwarfs or CH subgiants. Combining radial-velocity measurements from HERMES and SALT high-resolution spectra with radial-velocity data from CORAVEL and CORALIE, we determine the orbital parameters of 27 systems. We also derive their masses by comparing their location in the Hertzsprung-Russell diagram with evolutionary models. We confirm that Ba dwarfs and CH subgiants are not at different evolutionary stages and have similar metallicities, despite their different names. Additionally, Ba giants appear significantly more massive than their main-sequence analogues. This is likely due to observational biases against the detection of hotter main-sequence post-mass-transfer objects. Combining our spectroscopic orbits with the Hipparcos astrometric data, we derive the orbital inclinations and the mass of the WD companion for four systems. Since this cannot be done for all systems in our sample yet (but should be with upcoming Gaia data releases), we also analyse the mass-function distribution of our binaries. We can model this distribution with very narrow mass distributions for the two components and random orbital orientation on the sky. Finally, based on BINSTAR evolutionary models, we suggest that the orbital evolution of low-mass Ba systems can be affected by a second phase of interaction along the Red Giant Branch of the Ba star, impacting on the eccentricities and periods of the giants. △ Less

Submitted 8 April, 2019; originally announced April 2019.

Comments: Accepted for publication in A&A on the 5th of April, 2019

Journal ref: A&A 626, A128 (2019)

Astro2020 Science White Paper: gravity-wave asteroseismology of intermediate- and high-mass stars

Authors: Andrew Tkachenko, Conny Aerts, Dominic M. Bowman, Timothy Van Reeth, Joris De Ridder, Cole Campbell Johnston, May Gade Pedersen, Siemen Burssens, Mathias Michielsen, Joey Mombarg, Sanjay Sekaran, Robin Bjorklund, Tamara Rogers, Philipp Edelmann, Rathish Previn Ratnasingam, Konstanze Zwintz, Juna Kollmeier, Jennifer Johnson, Hans-Walter Rix, Jamie Tayar

Abstract: The evolution of a star is driven by the physical processes in its interior making the theory of stellar structure and evolution the most crucial ingredient for not only stellar evolution studies, but any field of astronomy which relies on the yields along stellar evolution. High-precision time-series photometric data assembled by recent space missions revealed that current models of stellar struc… ▽ More The evolution of a star is driven by the physical processes in its interior making the theory of stellar structure and evolution the most crucial ingredient for not only stellar evolution studies, but any field of astronomy which relies on the yields along stellar evolution. High-precision time-series photometric data assembled by recent space missions revealed that current models of stellar structure and evolution show major shortcomings already in the two earliest nuclear burning phases, impacting all subsequent phases prior to the formation of the end-of-life remnant. This white paper focuses specifically on the transport of chemical elements and of angular momentum in the stellar structure and evolution models of stars born with convective core, as revealed by their gravity-mode oscillations. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: 8 pages. White paper submitted to the Astro2020 Decadal Survey

Two's a crowd? Characterising the effect of photometric contamination on the extraction of the global asteroseismic parameter $ν_{\text{max}}$ in red-giant binaries

Authors: Sanjay Sekaran, Cole Johnston, Andrew Tkachenko, Paul G. Beck, Andrej Prsa, Kelly M. Hambleton

Abstract: Theoretical scaling relations for solar-like oscillators and red giants are widely used to estimate fundamental stellar parameters. The accuracy and precision of these relations have often been questioned in the literature, with studies often utilising binarity for model-independent validation. However, it has not been tested if the photometric effects of binarity introduce a systematic effect on… ▽ More Theoretical scaling relations for solar-like oscillators and red giants are widely used to estimate fundamental stellar parameters. The accuracy and precision of these relations have often been questioned in the literature, with studies often utilising binarity for model-independent validation. However, it has not been tested if the photometric effects of binarity introduce a systematic effect on the extraction of the seismic properties of the pulsating component(s). In this paper, we present an estimation of the impact of a contaminating photometric signal with a distinct background profile on the global asteroseismic parameter $ν_{\text{max}}$ through the analysis of synthetic red-giant binary light curves. We generate the pulsational and granulation parameters for single red giants with different masses, radii and effective temperatures from theoretical scaling relations and use them to simulate single red-giant light curves with the characteristics of Kepler long-cadence photometric data. These are subsequently blended together according to their light ratio to generate binary red-giant light curves of various configurations. We then perform a differential analysis to characterise the systematic effects of binarity on the extraction of $ν_{\text{max}}$. We find that the $ν_{\text{max}}$ extraction for red-giant power spectra featuring overlapping power excesses is unreliable if unconstrained priors are used. Outside of this scenario, we obtain results that are nearly identical to single-star case. We conclude that i) the photometric effects of binarity on the extraction of $ν_{\text{max}}$ are largely negligible as long as the power excesses of the individual components do not overlap, and that ii) there is minimal advantage to using more than two super-Lorentzian components to model the granulation signal of a binary red-giant. △ Less

Submitted 22 March, 2019; v1 submitted 21 March, 2019; originally announced March 2019.

Comments: 17 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 624, A140 (2019)

A Fast and Accurate Method to Capture the Solar Corona/Transition Region Enthalpy Exchange

Authors: C. D. Johnston, S. J. Bradshaw

Abstract: The brightness of the emission from coronal loops in the solar atmosphere is strongly dependent on the temperature and density of the confined plasma. After a release of energy, these loops undergo a heating and upflow phase, followed by a cooling and downflow cycle. Throughout, there are significant variations in the properties of the coronal plasma. In particular, the increased coronal temperatu… ▽ More The brightness of the emission from coronal loops in the solar atmosphere is strongly dependent on the temperature and density of the confined plasma. After a release of energy, these loops undergo a heating and upflow phase, followed by a cooling and downflow cycle. Throughout, there are significant variations in the properties of the coronal plasma. In particular, the increased coronal temperature leads to an excess downward heat flux that the transition region is unable to radiate. This generates an enthalpy flux from the transition region to the corona, increasing the coronal density. The enthalpy exchange is highly sensitive to the transition region resolution in numerical simulations. With a numerically under-resolved transition region, major errors occur in simulating the coronal density evolution and, thus, the predicted loop emission. This Letter presents a new method that addresses the difficulty of obtaining the correct interaction between the corona and corona/chromosphere interface. In the transition region, an adaptive thermal conduction approach is used that broadens any unresolved parts of the atmosphere. We show that this approach, referred to as TRAC, successfully removes the influence of numerical resolution on the coronal density response to heating while maintaining high levels of agreement with fully resolved models. When employed with coarse spatial resolutions, typically achieved in multi-dimensional MHD codes, the peak density errors are less than 3% and the computation time is three orders of magnitude faster than fully resolved field-aligned models. The advantages of using TRAC in field-aligned hydrodynamic and multi-dimensional magnetohydrodynamic simulations are discussed. △ Less

Submitted 4 March, 2019; originally announced March 2019.

Comments: 8 pages, 3 figures, accepted for publication in ApJL

Asteroseismology of massive stars with the TESS mission: the runaway Beta Cep pulsator PHL 346 = HN Aqr

Authors: Gerald Handler, Andrzej Pigulski, Jadwiga Daszyńska-Daszkiewicz, Andreas Irrgang, David Kilkenny, Zhao Guo, Norbert Przybilla, Filiz Kahraman Aliçavuş, Thomas Kallinger, Javier Pascual-Granado, Ewa Niemczura, Tomasz Różański, Sowgata Chowdhury, Derek L. Buzasi, Giovanni M. Mirouh, Dominic M. Bowman, Cole Johnston, May G. Pedersen, Sergio Simón-Diaz, Ehsan Moravveji, Kosmas Gazeas, Peter De Cat, Roland K. Vanderspek, George R. Ricker

Abstract: We report an analysis of the first known Beta Cep pulsator observed by the TESS mission, the runaway star PHL 346 = HN Aqr. The star, previously known as a singly-periodic pulsator, has at least 34 oscillation modes excited, 12 of those in the g-mode domain and 22 p modes. Analysis of archival data implies that the amplitude and frequency of the dominant mode and the stellar radial velocity were v… ▽ More We report an analysis of the first known Beta Cep pulsator observed by the TESS mission, the runaway star PHL 346 = HN Aqr. The star, previously known as a singly-periodic pulsator, has at least 34 oscillation modes excited, 12 of those in the g-mode domain and 22 p modes. Analysis of archival data implies that the amplitude and frequency of the dominant mode and the stellar radial velocity were variable over time. A binary nature would be inconsistent with the inferred ejection velocity from the Galactic disc of 420 km/s, which is too large to be survivable by a runaway binary system. A kinematic analysis of the star results in an age constraint (23 +- 1 Myr) that can be imposed on asteroseismic modelling and that can be used to remove degeneracies in the modelling process. Our attempts to match the excitation of the observed frequency spectrum resulted in pulsation models that were too young. Hence, asteroseismic studies of runaway pulsators can become vital not only in tracing the evolutionary history of such objects, but to understand the interior structure of massive stars in general. TESS is now opening up these stars for detailed asteroseismic investigation. △ Less

Submitted 21 February, 2019; originally announced February 2019.

Comments: accepted for ApJL

Asteroseismic masses, ages, and core properties of $γ$ Doradus stars using gravito-inertial dipole modes and spectroscopy

Authors: Joey S. G. Mombarg, Timothy Van Reeth, May G. Pedersen, Geert Molenberghs, Dominic M. Bowman, Cole Johnston, Andrew Tkachenko, Conny Aerts

Abstract: The asteroseismic modelling of period spacing patterns from gravito-inertial modes in stars with a convective core is a high-dimensional problem. We utilise the measured period spacing pattern of prograde dipole gravity modes (acquiring $Π_0$), in combination with the effective temperature ($T_{\rm eff}$) and surface gravity ($\log g$) derived from spectroscopy, to estimate the fundamental stellar… ▽ More The asteroseismic modelling of period spacing patterns from gravito-inertial modes in stars with a convective core is a high-dimensional problem. We utilise the measured period spacing pattern of prograde dipole gravity modes (acquiring $Π_0$), in combination with the effective temperature ($T_{\rm eff}$) and surface gravity ($\log g$) derived from spectroscopy, to estimate the fundamental stellar parameters and core properties of 37 $γ~$Doradus ($γ~$Dor) stars whose rotation frequency has been derived from $\textit{Kepler}$ photometry. We make use of two 6D grids of stellar models, one with step core overshooting and one with exponential core overshooting, to evaluate correlations between the three observables $Π_0$, $T_{\rm eff}$, and $\log g$ and the mass, age, core overshooting, metallicity, initial hydrogen mass fraction and envelope mixing. We provide multivariate linear model recipes relating the stellar parameters to be estimated to the three observables ($Π_0$, $T_{\rm eff}$, $\log g$). We estimate the (core) mass, age, core overshooting and metallicity of $γ~$Dor stars from an ensemble analysis and achieve relative uncertainties of $\sim\!10$ per cent for the parameters. The asteroseismic age determination allows us to conclude that efficient angular momentum transport occurs already early on during the main sequence. We find that the nine stars with observed Rossby modes occur across almost the entire main-sequence phase, except close to core-hydrogen exhaustion. Future improvements of our work will come from the inclusion of more types of detected modes per star, larger samples, and modelling of individual mode frequencies. △ Less

Submitted 20 February, 2019; v1 submitted 18 February, 2019; originally announced February 2019.

Comments: Accepted for publication in MNRAS

Combined asteroseismology, spectroscopy, and astrometry of the CoRoT B2V target HD 170580

Authors: C. Aerts, M. G. Pedersen, E. Vermeyen, L. Hendriks, C. Johnston, A. Tkachenko, P. I. Pápics, J. Debosscher, M. Briquet, A. Thoul, M. Rainer, E. Poretti

Abstract: Space asteroseismology reveals that stellar structure and evolution models of intermediate- and high-mass stars are in need of improvement in terms of angular momentum and chemical element transport. We aim to probe the interior structure of a hot massive star in the core-hydrogen burning phase of its evolution. We analyse CoRoT space photometry, Gaia DR2 space astrometry, and high-resolution high… ▽ More Space asteroseismology reveals that stellar structure and evolution models of intermediate- and high-mass stars are in need of improvement in terms of angular momentum and chemical element transport. We aim to probe the interior structure of a hot massive star in the core-hydrogen burning phase of its evolution. We analyse CoRoT space photometry, Gaia DR2 space astrometry, and high-resolution high signal-to-noise HERMES and HARPS time-series spectroscopy of the slowly rotating B2V star HD 170580. From the time-series spectroscopy we derive $vsini=4\pm 2\,$km\,s$^{-1}$, where the uncertainty results from the complex pulsational line-profile variability. We detect 42 frequencies with amplitude above 5$σ$. Among these we identify 5 rotationally split triplets and 1 quintuplet. Asteroseismic modelling based on CoRoT, Gaia DR2 and spectroscopic data leads to a star of 8M$_\odot$ near core-hydrogen exhaustion and an extended overshoot zone. The detected low-order pressure-mode frequencies cannot be fit within the uncertainties of the CoRoT data by models without atomic diffusion. Irrespective of this limitation, the low-order gravity modes reveal HD 170580 to be a slow rotator with an average rotation period between 73 and 98 d and a hint of small differential rotation. Future Gaia DR3 data taking into account the multiplicity of the star, along with long-term TESS photometry would allow to put better observational constraints on the asteroseismic models of this blue evolved massive star. Improved modelling with atomic diffusion, including radiative levitation, is needed to achieve compliance with the low helium surface abundance of the star. This poses immense computational challenges but is required to derive the interior rotation and mixing profiles of this star. △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: Manuscript accepted for publication in Astronomy & Astrophysics; 13 pages, 8 figures, 2 tables

Journal ref: A&A 624, A75 (2019)

Diverse Variability of O and B Stars Revealed from 2-minute Cadence Light Curves in Sectors 1 and 2 of the TESS Mission: Selection of an Asteroseismic Sample

Authors: May G. Pedersen, Sowgata Chowdhury, Cole Johnston, Dominic Bowman, Conny Aerts, Gerald Handler, Peter De Cat, Coralie Neiner, Alexandre David-Uraz, Derek Buzasi, Andrew Tkachenko, Sergio Simon-Diaz, Ehsan Moravveji, James Sikora, Giovanni M. Mirouh, Catherine C. Lovekin, Matteo Cantiello, Jadwiga Daszynska-Daszkiewicz, Andrzej Pigulski, Roland K. Vanderspek, George R. Ricker

Abstract: Uncertainties in stellar structure and evolution theory are largest for stars undergoing core convection on the main sequence. A powerful way to calibrate the free parameters used in the theory of stellar interiors is asteroseismology, which provides direct measurements of angular momentum and element transport. We report the detection and classification of new variable O and B stars using high-pr… ▽ More Uncertainties in stellar structure and evolution theory are largest for stars undergoing core convection on the main sequence. A powerful way to calibrate the free parameters used in the theory of stellar interiors is asteroseismology, which provides direct measurements of angular momentum and element transport. We report the detection and classification of new variable O and B stars using high-precision short-cadence (2-min) photometric observations assembled by the Transiting Exoplanet Survey Satellite (TESS). In our sample of 154 O and B stars, we detect a high percentage (90%) of variability. Among these we find 23 multiperiodic pulsators, 6 eclipsing binaries, 21 rotational variables, and 25 stars with stochastic low-frequency variability. Several additional variables overlap between these categories. Our study of O and B stars not only demonstrates the high data quality achieved by TESS for optimal studies of the variability of the most massive stars in the Universe, but also represents the first step towards the selection and composition of a large sample of O and B pulsators with high potential for joint asteroseismic and spectroscopic modeling of their interior structure with unprecedented precision. △ Less

Submitted 24 January, 2019; v1 submitted 22 January, 2019; originally announced January 2019.

Comments: 16 pages, 4 figures, 1 Table, Accepted for publication in ApJL

Photometric detection of internal gravity waves in early-type stars observed by CoRoT

Authors: D. M. Bowman, C. Aerts, C. Johnston, M. G. Pedersen, T. M. Rogers, P. V. F. Edelmann, S. Simón-Díaz, T. Van Reeth, B. Buysschaert, A. Tkachenko, S. A. Triana

Abstract: Early-type stars are predicted to excite an entire spectrum of internal gravity waves (IGWs) at the interface of their convective cores and radiative envelopes. Numerical simulations of IGWs predict stochastic low-frequency variability in photometric observations, yet the detection of IGWs in early-type stars has been limited by a dearth of high-quality photometric time series. We present observat… ▽ More Early-type stars are predicted to excite an entire spectrum of internal gravity waves (IGWs) at the interface of their convective cores and radiative envelopes. Numerical simulations of IGWs predict stochastic low-frequency variability in photometric observations, yet the detection of IGWs in early-type stars has been limited by a dearth of high-quality photometric time series. We present observational evidence of stochastic low-frequency variability in the CoRoT photometry of a sample of O, B, A and F stars. The presence of this stochastic low-frequency variability in stars across the upper main-sequence cannot be universally explained as granulation or stellar winds, but its morphology is found to be consistent with predictions from IGW simulations. △ Less

Submitted 30 November, 2018; originally announced November 2018.

Comments: 6 pages, 3 figures, Proceedings of the PHOST "Physics of Oscillating Stars" conference (2-7 Sept. 2018, Banyuls-sur-mer, France), Edited by J. Ballot, S. Vauclair, & G. Vauclair

On the photometric detection of Internal Gravity Waves in upper main-sequence stars I. Methodology and application to CoRoT targets

Authors: D. M. Bowman, C. Aerts, C. Johnston, M. G. Pedersen, T. M. Rogers, P. V. F. Edelmann, S. Simón-Díaz, T. Van Reeth, B. Buysschaert, A. Tkachenko, S. A. Triana

Abstract: Context. Main sequence stars with a convective core are predicted to stochastically excite Internal Gravity Waves (IGWs), which effectively transport angular momentum throughout the stellar interior and explain the observed near-uniform interior rotation rates of intermediate-mass stars. However, there are few detections of IGWs, and fewer still made using photometry, with more detections needed t… ▽ More Context. Main sequence stars with a convective core are predicted to stochastically excite Internal Gravity Waves (IGWs), which effectively transport angular momentum throughout the stellar interior and explain the observed near-uniform interior rotation rates of intermediate-mass stars. However, there are few detections of IGWs, and fewer still made using photometry, with more detections needed to constrain numerical simulations. Aims. We aim to formalise the detection and characterisation of IGWs in photometric observations of stars born with convective cores (M > 1.5 M$_{\odot}$) and parameterise the low-frequency power excess caused by IGWs. Methods. Using the most recent CoRoT light curves for a sample of O, B, A and F stars, we parameterise the morphology of the flux contribution of IGWs in Fourier space using an MCMC numerical scheme within a Bayesian framework. We compare this to predictions from IGW numerical simulations and investigate how the observed morphology changes as a function of stellar parameters. Results. We demonstrate that a common morphology for the low-frequency power excess is observed in early-type stars observed by CoRoT. Our study shows that a background frequency-dependent source of astrophysical signal is common, which we interpret as IGWs. We provide constraints on the amplitudes of IGWs and the shape of their detected frequency spectrum across a range of mass, which is the first ensemble study of stochastic variability in such a diverse sample of stars. Conclusions. The evidence of a low-frequency power excess across a wide mass range supports the interpretation of IGWs in photometry of O, B, A and F stars. We also discuss the prospects of observing hundreds of massive stars with the Transiting Exoplanet Survey Satellite (TESS) in the near future. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Comments: Accepted for publication in A&A, 41 pages including 11 figures and 32 appendix figures

Journal ref: A&A 621, A135 (2019)

Binary Asteroseismic Modelling: isochrone-cloud methodology and application to Kepler gravity-mode pulsators

Authors: C. Johnston, A. Tkachenko, C. Aerts, G. Molenberghs, D. M. Bowman, M. G. Pedersen, B. Buysschaert, P. I. Papics

Abstract: The simultaneous presence of variability due to both pulsations and binarity is no rare phenomenon. Unfortunately, the complexities of dealing with even one of these sources of variability individually means that the other signal is often treated as a nuisance and discarded. However, both types of variability offer means to probe fundamental stellar properties in robust ways through asteroseismic… ▽ More The simultaneous presence of variability due to both pulsations and binarity is no rare phenomenon. Unfortunately, the complexities of dealing with even one of these sources of variability individually means that the other signal is often treated as a nuisance and discarded. However, both types of variability offer means to probe fundamental stellar properties in robust ways through asteroseismic and binary modelling. We present an efficient methodology that includes both binary and asteroseismic information to estimate fundamental stellar properties based on a grid-based modelling approach. We report parameters for three gravity mode pulsating {\it Kepler} binaries , such as mass, radius, age, as well the mass of the convective core and location of the overshoot region. We discuss the presence of parameter degeneracies and the way our methodology deals with them. We provide asteroseismically calibrated isochrone-clouds to the community; these are a generalisation of isochrones when allowing for different values of the core overshooting in the two components of the binary. △ Less

Submitted 1 October, 2018; originally announced October 2018.

Comments: 16 pages, 15 figures, 6 tables. Accepted for publication in MNRAS