Showing 1–50 of 134 results for author: Reese, D

An improved asteroseismic age of the rapid rotator Altair from TESS data

Authors: Michel Rieutord, Daniel R. Reese, Joey S. G. Mombarg, Stéphane Charpinet

Abstract: Understanding the effects of rotation in stellar evolution is key to modelling early-type stars, half of which have equatorial velocities over 100 km/s. The nearby star Altair is an example of such fast-rotating stars, and furthermore, it has the privilege of being modelled by a detailed 2D concordance model that reproduces most of its observables. The aim of this paper is to include new asterosei… ▽ More Understanding the effects of rotation in stellar evolution is key to modelling early-type stars, half of which have equatorial velocities over 100 km/s. The nearby star Altair is an example of such fast-rotating stars, and furthermore, it has the privilege of being modelled by a detailed 2D concordance model that reproduces most of its observables. The aim of this paper is to include new asteroseismic frequencies to improve our knowledge of Altair, especially its age. We processed images of Altair obtained during July 2022 by the Transiting Exoplanet Survey Satellite using the halo photometry technique to obtain its light curve over this observation period. By analysing the light curve, we derived a set of 22 new frequencies in the oscillation spectrum of Altair and confirmed 12 previously known frequencies. Compared with model predictions, we could associate ten frequencies with ten axisymmetric modes. This identification is based on the modelled visibility of the modes. Moreover, nine of the modelled frequencies can be adjusted to simultaneously match their corresponding observed frequencies, once the core hydrogen mass fraction of the concordance model is set to $X_{\rm core}/X_{\rm ini}\simeq0.972$, with $X_{\rm ini}=0.739$. Using the combined results of a 1D MESA model computing the pre-main sequence and a 2D time-dependent ESTER model computing the main sequence, we find that this core hydrogen abundance sets the age of Altair to 88$\pm$10 Myrs, which is slightly younger than previous estimates. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 9 pages, 8 figures, accepted for publication in A&A

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.

Predicted asteroseismic detection yield for solar-like oscillating stars with PLATO

Authors: M. J. Goupil, C. Catala, R. Samadi, K. Belkacem, R. M. Ouazzani, D. R. Reese, T. Appourchaux, S. Mathur, J. Cabrera, A. Börner, C. Paproth, N. Moedas, K. Verma, Y. Lebreton, M. Deal, J. Ballot, W. J. Chaplin, J. Christensen-Dalsgaard, M. Cunha, A. F. Lanza, A. Miglio, T. Morel, A. Serenelli, B. Mosser, O. Creevey , et al. (4 additional authors not shown)

Abstract: We determine the expected yield of detections of solar-like oscillations for the PLATO ESA mission. We used a formulation from the literature to calculate the probability of detection and validated it with Kepler data. We then applied this approach to the PLATO P1 and P2 samples with the lowest noise level and the much larger P5 sample, which has a higher noise level. We used the information avail… ▽ More We determine the expected yield of detections of solar-like oscillations for the PLATO ESA mission. We used a formulation from the literature to calculate the probability of detection and validated it with Kepler data. We then applied this approach to the PLATO P1 and P2 samples with the lowest noise level and the much larger P5 sample, which has a higher noise level. We used the information available in in the PIC 1.1.0, including the current best estimate of the signal-to-noise ratio. We also derived relations to estimate the uncertainties of seismically inferred stellar mass, radius and age and applied those relations to the main sequence stars of the PLATO P1 and P2 samples with masses equal to or below 1.2 $\rm{M}_\odot$ for which we had obtained a positive seismic detection. We found that one can expect positive detections of solar-like oscillations for more than 15 000 FGK stars in one single field after a two-years run of observation. For main sequence stars with masses $\leq 1.2 \rm{M}_\odot$, we found that about 1131 stars satisfy the PLATO requirements for the uncertainties of the seismically inferred stellar masses, radii and ages in one single field after a two-year run of observation. The baseline observation programme of PLATO consists in observing two fields of similar size (in the Southern and Northern hemispheres) for two years each. The expected seismic yields of the mission are more 30000 FGK dwarfs and subgiants with positive detections of solar-like oscillations, enabling to achieve the mission stellar objectives. The PLATO mission should produce a sample of seismically extremely well characterized stars of quality equivalent to the Kepler Legacy sample but containing a number of stars $\sim$ 80 times larger if observing two PLATO fields for two years each. They will represent a goldmine which will make possible significant advances in stellar modelling. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: 25 pages, 26 figure. Accepted in Astronomy and Astrophysics

Asteroseismic modelling strategies in the PLATO era. II. Automation of seismic inversions and quality assessment procedure

Authors: Jérôme Bétrisey, Gaël Buldgen, Daniel R. Reese, Georges Meynet

Abstract: *Context*. In the framework of the PLATO mission, to be launched in late 2026, seismic inversion techniques will play a key role in the mission precision requirements of the stellar mass, radius, and age. It is therefore relevant to discuss the challenges of the automation of seismic inversions, which were originally developed for individual modelling.\\ *Aims*. We tested the performance of our ne… ▽ More *Context*. In the framework of the PLATO mission, to be launched in late 2026, seismic inversion techniques will play a key role in the mission precision requirements of the stellar mass, radius, and age. It is therefore relevant to discuss the challenges of the automation of seismic inversions, which were originally developed for individual modelling.\\ *Aims*. We tested the performance of our newly developed quality assessment procedure of seismic inversions, which was designed in the perspective of a pipeline implementation.\\ *Methods*. We applied our assessment procedure on a testing set composed of 26 reference models. We divided our testing set into two categories, calibrator targets whose inversion behaviour is well known from the literature and targets for which we assessed manually the quality of the inversion. We then compared the results of our assessment procedure with our expectations as a human modeller for three types of inversions, the mean density inversion, the acoustic radius inversion, and the central entropy inversion.\\ *Results*. We found that our quality assessment procedure performs as well as a human modeller. The mean density inversion and the acoustic radius inversion are suited for a large-scale application, but not the central entropy inversion, at least in its current form.\\ *Conclusions*. Our assessment procedure showed promising results for a pipeline implementation. It is based on by-products of the inversion and therefore requires few numerical resources to assess quickly the quality of an inversion result. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

PlatoSim: An end-to-end PLATO camera simulator for modelling high-precision space-based photometry

Authors: N. Jannsen, J. De Ridder, D. Seynaeve, S. Regibo, R. Huygen, P. Royer, C. Paproth, D. Grießbach, R. Samadi, D. R. Reese, M. Pertenais, E. Grolleau, R. Heller, S. M. Niemi, J. Cabrera, A. Börner, S. Aigrain, J. McCormac, P. Verhoeve, P. Astier, N. Kutrowski, B. Vandenbussche, A. Tkachenko, C. Aerts

Abstract: PLAnetary Transits and Oscillations of stars (PLATO) is the ESA M3 space mission dedicated to detect and characterise transiting exoplanets including information from the asteroseismic properties of their stellar hosts. The uninterrupted and high-precision photometry provided by space-borne instruments such as PLATO require long preparatory phases. An exhaustive list of tests are paramount to desi… ▽ More PLAnetary Transits and Oscillations of stars (PLATO) is the ESA M3 space mission dedicated to detect and characterise transiting exoplanets including information from the asteroseismic properties of their stellar hosts. The uninterrupted and high-precision photometry provided by space-borne instruments such as PLATO require long preparatory phases. An exhaustive list of tests are paramount to design a mission that meets the performance requirements, and as such, simulations are an indispensable tool in the mission preparation. To accommodate PLATO's need of versatile simulations prior to mission launch - that at the same time describe accurately the innovative but complex multi-telescope design - we here present the end-to-end PLATO simulator specifically developed for the purpose, namely PlatoSim. We show step-by-step the algorithms embedded into the software architecture of PlatoSim that allow the user to simulate photometric time series of CCD images and light curves in accordance to the expected observations of PLATO. In the context of the PLATO payload, a general formalism of modelling, end-to-end, incoming photons from the sky to the final measurement in digital units is discussed. We show the strong predictive power of PlatoSim through its diverse applicability and contribution to numerous working groups within the PLATO Mission Consortium. This involves the on-going mechanical integration and alignment, performance studies of the payload, the pipeline development and assessments of the scientific goals. PlatoSim is a state-of-the-art simulator that is able to produce the expected photometric observations of PLATO to a high level of accuracy. We demonstrate that PlatoSim is a key software tool for the PLATO mission in the preparatory phases until mission launch and prospectively beyond. △ Less

Submitted 8 May, 2024; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: 30 pages, 22 figures

Journal ref: A&A 681, A18 (2024)

Asteroseismic modelling strategies in the PLATO era I. Mean density inversions and direct treatment of the seismic information

Authors: Jérôme Bétrisey, Gaël Buldgen, Daniel R. Reese, Martin Farnir, Marc-Antoine Dupret, Saniya Khan, Marie-Jo Goupil, Patrick Eggenberger, Georges Meynet

Abstract: Asteroseismic modelling will be part of the pipeline of the PLATO mission and will play a key role in the mission precision requirements on stellar mass, radius and age. It is therefore crucial to compare how current modelling strategies perform, and discuss the limitations and remaining challenges for PLATO, such as the so-called surface effects, the choice of physical ingredients, and stellar ac… ▽ More Asteroseismic modelling will be part of the pipeline of the PLATO mission and will play a key role in the mission precision requirements on stellar mass, radius and age. It is therefore crucial to compare how current modelling strategies perform, and discuss the limitations and remaining challenges for PLATO, such as the so-called surface effects, the choice of physical ingredients, and stellar activity. In this context, we carried out a systematic study of the impact of surface effects on the estimation of stellar parameters. In this work, we demonstrated how combining a mean density inversion with a fit of frequencies separation ratios can efficiently damp the surface effects and achieve precise and accurate stellar parameters for ten Kepler LEGACY targets, well within the PLATO mission requirements. We applied and compared two modelling approaches, directly fitting the individual frequencies, or coupling a mean density inversion with a fit of the ratios, to six synthetic targets with a patched 3D atmosphere from Sonoi et al. (2015) and ten actual targets from the LEGACY sample. The fit of the individual frequencies is unsurprisingly very sensitive to surface effects and the stellar parameters tend to be biased, which constitutes a fundamental limit to both accuracy and precision. In contrast, coupling a mean density inversion and a fit of the ratios efficiently damps the surface effects, and allows us to get both precise and accurate stellar parameters. The average statistical precision of our selection of LEGACY targets with this second strategy is 1.9% for the mass, 0.7% for the radius, and 4.1% for the age, well within the PLATO requirements. Using the mean density in the constraints significantly improves the precision of the mass, radius and age determinations, on average by 20%, 33%, and 16%, respectively. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 676, A10 (2023)

RUBIS: a simple tool for calculating the centrifugal deformation of stars and planets

Authors: Pierre S. Houdayer, Daniel R. Reese

Abstract: We present RUBIS (Rotation code Using Barotropy conservation over Isopotential Surfaces), a fully Python-based centrifugal deformation program available at \url{https://github.com/pierrehoudayer/RUBIS}. The code has been designed to calculate the centrifugal deformation of stars and planets resulting from a given cylindrical rotation profile, starting from a spherically symmetric non-rotating mode… ▽ More We present RUBIS (Rotation code Using Barotropy conservation over Isopotential Surfaces), a fully Python-based centrifugal deformation program available at \url{https://github.com/pierrehoudayer/RUBIS}. The code has been designed to calculate the centrifugal deformation of stars and planets resulting from a given cylindrical rotation profile, starting from a spherically symmetric non-rotating model. The underlying assumption in RUBIS is that the relationship between density and pressure is preserved during the deformation process. This leads to many procedural simplifications. For instance, RUBIS only needs to solve Poisson' equation, either in spheroidal or spherical coordinates depending on whether the 1D model has discontinuities or not. We present the benefits of using RUBIS to deform polytropic models and more complex barotropic structures, thus providing, to a certain extent, insights into baroclinic models. The resulting structures can be used for a wide range of applications, including the seismic study of models. Finally, we illustrate how RUBIS is beneficial specifically in the analysis of Jupiter's gravitational moments, thanks to its ability to handle discontinuous models while retaining a high accuracy compared to current methods. △ Less

Submitted 25 April, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

Journal ref: A&A 675, A181 (2023)

Spectroscopic detection of Altair's non-radial pulsations

Authors: M. Rieutord, P. Petit, D. Reese, T. Böhm, A. Lopez Ariste, G. Mirouh, A. Domiciano de Souza

Abstract: Context: Rapid rotation is a common feature of early-type stars but which remains a challenge for the models. The understanding of its effect on stellar evolution is however imperative to interpret the observed properties of numerous stars. Aims: We wish to bring more observational constraints on the properties of fast rotating stars, especially on their oscillation modes. Methods: We focus on the… ▽ More Context: Rapid rotation is a common feature of early-type stars but which remains a challenge for the models. The understanding of its effect on stellar evolution is however imperative to interpret the observed properties of numerous stars. Aims: We wish to bring more observational constraints on the properties of fast rotating stars, especially on their oscillation modes. Methods: We focus on the nearby star Altair which is known as a very rapidly rotating star with an equatorial velocity estimated recently at 313 km/s. We observed this star with the high-resolution spectropolarimeter Neo-Narval during six nights, with one night of interruption, in September 2020. Results: We detect significant line profile variations on the mean line profile of the spectra. Their time-frequency analysis shows that these variations are induced by gravito-inertial waves propagating at Altair's surface with azimuthal wavenumbers of order $m=10-15$. With a preliminary computation of the eigenspectrum using the most recent concordance model of Altair we can give a first modelling of the observed waves. Conclusions: Altair was known as the brightest $δ$ Scuti star. We now see that it is the brightest hybrid oscillating star with excited gravito-inertial waves and acoustic waves. Clearly, more observations and more advanced models are needed to explain the observations in greater details △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 8 pages, 8 figures, submitted to Astronomy and Astrophysics

Journal ref: A&A 669, A99 (2023)

Inversions of stellar structure from asteroseismic data

Authors: Gaël Buldgen, Jérôme Bétrisey, Ian W. Roxburgh, Sergei V. Vorontsov, Daniel R. Reese

Abstract: The advent of space-based photometry missions in the early 21st century enabled the application to asteroseismic data of advanced inference techniques until then restricted to the field of helioseismology. The high quality of the observations, the discovery of mixed modes in evolved solar-like oscillators and the need for an improvement in the determination of stellar fundamental parameters such a… ▽ More The advent of space-based photometry missions in the early 21st century enabled the application to asteroseismic data of advanced inference techniques until then restricted to the field of helioseismology. The high quality of the observations, the discovery of mixed modes in evolved solar-like oscillators and the need for an improvement in the determination of stellar fundamental parameters such as mass, radius and age led to the development of sophisticated modelling tools, amongst which seismic inversions play a key role. In this review, we will discuss the existing inversion techniques for the internal structure of distant stars adapted from helio- to asteroseismology. We will present results obtained for various Kepler targets, their coupling to other existing modelling techniques as well as the limitations of seismic analyses and the perspectives for future developments of these approaches in the context of the current TESS and the future PLATO mission, as well as the exploitation of the mixed modes observed in post-main sequence solar-like oscillators, for which variational formulations might not provide sufficient accuracy. △ Less

Submitted 23 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Frontiers in Astronomy and Space Sciences, special topic "Challenges of Asteroseismology in the Era of Space Missions", editors: Javier Pascual Granado, Tiago Campante, Antonio García Hernández, Zhao Guo

Properties of the ionisation glitch II. Seismic signature of the structural perturbation

Authors: Pierre S. Houdayer, Daniel R. Reese, Marie-Jo Goupil

Abstract: In the present paper, we aim to constrain the properties of the ionisation region in a star from the oscillation frequency variation (a so-called glitch) caused by rapid structural variations in this very region. In particular, we seek to avoid the use of calibration based on stellar models thus providing a truly independent estimate of these properties. These include both the helium abundance and… ▽ More In the present paper, we aim to constrain the properties of the ionisation region in a star from the oscillation frequency variation (a so-called glitch) caused by rapid structural variations in this very region. In particular, we seek to avoid the use of calibration based on stellar models thus providing a truly independent estimate of these properties. These include both the helium abundance and other physical quantities that can have a significant impact on the oscillation frequencies such as the electronic degeneracy parameter or the extent of the ionisation region. Taking as a starting point our first paper, we applied structural perturbations of the ionisation zone to the wave equation for radial oscillations in an isentropic region. The resulting glitch model is thus able to exploit the information contained in the fast frequency oscillation caused by the helium ionisation but also in the slow trend accompanying that of hydrogen. This information can directly be expressed in terms of parameters related respectively to the helium abundance, electronic degeneracy and extent of the ionisation region. Using a Bayesian inference, we show that a substantial recovery of the properties at the origin of the glitch is possible. A degeneracy between the helium abundance and the electronic degeneracy is found to exist, which particularly affects the helium estimate. Extending the method to cases where the glitch is subject to contamination (e.g. surface effects), we noted the importance of the slow glitch trend associated with hydrogen ionisation. We propose using a Gaussian process to disentangle the frequency glitch from surface effects. △ Less

Submitted 17 May, 2022; v1 submitted 9 February, 2022; originally announced February 2022.

Journal ref: A&A 663, A60 (2022)

The SAPP pipeline for the determination of stellar abundances and atmospheric parameters of stars in the core program of the PLATO mission

Authors: Matthew Raymond Gent, Maria Bergemann, Aldo Serenelli, Luca Casagrande, Jeffrey M. Gerber, Ulrike Heiter, Mikhail Kovalev, Thierry Morel, Nicolas Nardetto, Vardan Adibekyan, Víctor Silva Aguirre, Martin Asplund, Kevin Belkacem, Carlos del Burgo, Lionel Bigot, Andrea Chiavassa, Luisa Fernanda Rodríguez Díaz, Marie-Jo Goupil, Jonay I. González Hernández, Denis Mourard, Thibault Merle, Szabolcs Mészáros, Douglas J. Marshall, Rhita-Maria Ouazzani, Bertrand Plez , et al. (3 additional authors not shown)

Abstract: We introduce the SAPP (Stellar Abundances and atmospheric Parameters Pipeline), the prototype of the code that will be used to determine parameters of stars observed within the core program of the PLATO space mission. The pipeline is based on the Bayesian inference and provides effective temperature, surface gravity, metallicity, chemical abundances, and luminosity. The code in its more general ve… ▽ More We introduce the SAPP (Stellar Abundances and atmospheric Parameters Pipeline), the prototype of the code that will be used to determine parameters of stars observed within the core program of the PLATO space mission. The pipeline is based on the Bayesian inference and provides effective temperature, surface gravity, metallicity, chemical abundances, and luminosity. The code in its more general version can have a much wider range of applications. It can also provide masses, ages, and radii of stars and can be used for stars of stellar types not targeted by the PLATO core program, such as red giants. We validate the code on a set of 27 benchmark stars that includes 19 FGK-type dwarfs, 6 GK-type sub-giants, and 2 red giants. Our results suggest that combining various observables is the optimal approach, as it allows to break degeneracies between different parameters and yields more accurate values of stellar parameters and more realistic uncertainties. For the PLATO core sample, we obtain a typical uncertainty of 27 ($\rm{syst.}$) $\pm$ 37 ($\rm{stat.}$) K for T$_{\rm{eff}}$, 0.00 $\pm$ 0.01 dex for log$g$, 0.02 $\pm$ 0.02 dex for metallicity [Fe/H], -0.01 $\pm$ 0.03 R$_\odot$ for radii, -0.01 $\pm$ 0.05 M$_\odot$ for stellar masses, and -0.14 $\pm$ 0.63 Gyrs for ages. We also show that the best results are obtained by combining the $ν_{max}$ scaling relation and stellar spectra. This resolves the notorious problem of degeneracies, which is particularly important for F-type stars. △ Less

Submitted 20 January, 2022; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: 25 pages, 42 figures, accepted for publication in A&A

Journal ref: A&A 658, A147 (2022)

Five young $δ$ Scuti stars in the Pleiades seen with Kepler/K2

Authors: Simon J. Murphy, Timothy R. Bedding, Timothy R. White, Yaguang Li, Daniel Hey, Daniel Reese, Meridith Joyce

Abstract: We perform mode identification for five $δ$ Scuti stars in the Pleiades star cluster, using custom light curves from K2 photometry. By creating échelle diagrams, we identify radial and dipole mode ridges, comprising a total of 28 radial and 16 dipole modes across the five stars. We also suggest possible identities for those modes that lie offset from the radial and dipole ridges. We calculate non-… ▽ More We perform mode identification for five $δ$ Scuti stars in the Pleiades star cluster, using custom light curves from K2 photometry. By creating échelle diagrams, we identify radial and dipole mode ridges, comprising a total of 28 radial and 16 dipole modes across the five stars. We also suggest possible identities for those modes that lie offset from the radial and dipole ridges. We calculate non-rotating stellar pulsation models to verify our mode identifications, finding good agreement within the age and metallicity constraints of the cluster. We also find that for all stars, the least dense models are preferred, reflecting the lower density of these oblate, rotating stars. Three of the five stars show rotationally-split multiplets. We conclude that the sample shows promise for asteroseismic rotation rates, masses, and ages with rotating models in the future. Our preliminary modelling also indicates some sensitivity to the helium abundance. △ Less

Submitted 24 January, 2022; v1 submitted 7 November, 2021; originally announced November 2021.

Comments: Submitted to MNRAS

PLATO Hare-and-Hounds exercise: Asteroseismic model fitting of main-sequence solar-like pulsators

Authors: M. S. Cunha, I. W. Roxburgh, V. Aguirre Børsen-Koch, W. H. Ball, S. Basu, W. J. Chaplin, M. -J. Goupil, B. Nsamba, J. Ong, D. R. Reese, K. Verma, K. Belkacem, T. Campante, J. Christensen-Dalsgaard, M. T. Clara, S. Deheuvels, M. J. P. F. G. Monteiro, A. Noll, R. M. Ouazzani, J. L. Rørsted, A. Stokholm, M. L. Winther

Abstract: Asteroseismology is a powerful tool to infer fundamental stellar properties. The use of these asteroseismic-inferred properties in a growing number of astrophysical contexts makes it vital to understand their accuracy. Consequently, we performed a hare-and-hounds exercise where the hares simulated data for 6 artificial main-sequence stars and the hounds inferred their properties based on different… ▽ More Asteroseismology is a powerful tool to infer fundamental stellar properties. The use of these asteroseismic-inferred properties in a growing number of astrophysical contexts makes it vital to understand their accuracy. Consequently, we performed a hare-and-hounds exercise where the hares simulated data for 6 artificial main-sequence stars and the hounds inferred their properties based on different inference procedures. To mimic a pipeline such as that planned for the PLATO mission, all hounds used the same model grid. Some stars were simulated using the physics adopted in the grid, others a different one. The maximum relative differences found (in absolute value) between the inferred and true values of the mass, radius, and age were 4.32 per cent, 1.33 per cent, and 11.25 per cent, respectively. The largest systematic differences in radius and age were found for a star simulated assuming gravitational settling, not accounted for in the model grid, with biases of -0.88 per cent (radius) and 8.66 per cent (age). For the mass, the most significant bias (-3.16 per cent) was found for a star with a helium enrichment ratio outside the grid range. Moreover, a ~7 per cent dispersion in age was found when adopting different prescriptions for the surface corrections or shifting the classical observations by $\pm 1σ$. The choice of the relative weight given to the classical and seismic constraints also impacted significantly the accuracy and precision of the results. Interestingly, only a few frequencies were required to achieve accurate results on the mass and radius. For the age the same was true when at least one $l=2$ mode was considered. △ Less

Submitted 7 October, 2021; originally announced October 2021.

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

Properties of the ionisation glitch: I. Modelling the ionisation region

Authors: Pierre S. Houdayer, Daniel R. Reese, Marie-Jo Goupil, Yveline Lebreton

Abstract: Determining the properties of solar-like oscillating stars can be subject to many biases. A particularly important example is the helium-mass degeneracy, where the uncertainties regarding the internal physics can cause a poor determination of both the mass and surface helium content. Accordingly, an independent helium estimate is needed to overcome this degeneracy. A promising way to obtain such a… ▽ More Determining the properties of solar-like oscillating stars can be subject to many biases. A particularly important example is the helium-mass degeneracy, where the uncertainties regarding the internal physics can cause a poor determination of both the mass and surface helium content. Accordingly, an independent helium estimate is needed to overcome this degeneracy. A promising way to obtain such an estimate is to exploit the so-called ionisation glitch, i.e. a deviation from the asymptotic oscillation frequency pattern caused by the rapid structural variation in the He ionisation zones. Although progressively becoming more sophisticated, the glitch-based approach faces problems inherent to its current modelling such as the need for calibration by realistic stellar models. This requires a physical model of the ionisation region explicitly involving the parameters of interest such as the surface helium abundance, $Y_s$. Through a thermodynamic treatment of the ionisation region, an analytical approximation for the first adiabatic exponent $Γ_1$ is presented. The induced stellar structure is found to depend on only three parameters including the surface helium abundance $Y_s$ and the electron degeneracy $ψ_\textrm{CZ}$ in the convective region. The model thus defined allows a wide variety of structures to be described and, in particular, is able to approximate a realistic model in the ionisation region. The modelling work conducted enables us to study the structural perturbations causing the glitch. More elaborate forms of perturbations than the ones usually assumed are found. It is also suggested that there might be a stronger dependence of the structure on both the electron degeneracy in the convection zone and on the position of the ionisation region rather than on the amount of helium itself. △ Less

Submitted 1 October, 2021; originally announced October 2021.

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

Oscillations of 2D ESTER models. I. The adiabatic case

Authors: D. R. Reese, G. M. Mirouh, F. Espinosa Lara, M. Rieutord, B. Putigny

Abstract: Recent numerical and theoretical considerations have shown that low-degree acoustic modes in rapidly rotating stars follow an asymptotic formula and recent observations of pulsations in rapidly rotating delta Scuti stars seem to match these expectations. However, a key question is whether strong gradients or discontinuities can adversely affect this pattern to the point of hindering its identifica… ▽ More Recent numerical and theoretical considerations have shown that low-degree acoustic modes in rapidly rotating stars follow an asymptotic formula and recent observations of pulsations in rapidly rotating delta Scuti stars seem to match these expectations. However, a key question is whether strong gradients or discontinuities can adversely affect this pattern to the point of hindering its identification. Other important questions are how rotational splittings are affected by the 2D rotation profiles expected from baroclinic effects and whether it is possible to probe the rotation profile using these splittings. Accordingly, we numerically calculate pulsation modes in continuous and discontinuous rapidly rotating models produced by the 2D ESTER (Evolution STEllaire en Rotation) code. This spectral multi-domain code self-consistently calculates the rotation profile based on baroclinic effects and allows us to introduce discontinuities without loss of numerical accuracy. Pulsations are calculated using an adiabatic version of the Two-dimensional Oscillation Program (TOP) code. The variational principle is used to confirm the high accuracy of the pulsation frequencies and to derive an integral formula that closely matches the generalised rotational splittings, except when modes are involved in avoided crossings. This potentially allows us to probe the the rotation profile using inverse theory. Acoustic glitch theory, applied along the island mode orbit deduced from ray dynamics, can correctly predict the periodicity of the glitch frequency pattern produced by a discontinuity or the Gamma1 dip related to the He II ionisation zone in some of the models. The asymptotic frequency pattern remains sufficiently well preserved to potentially allow its detection in observed stars. △ Less

Submitted 21 October, 2020; originally announced October 2020.

Comments: 21 pages, 20 figures, 6 tables, accepted for publication in A&A

Journal ref: A&A 645, A46 (2021)

SPInS, a pipeline for massive stellar parameter inference: A public Python tool to age-date, weigh, size up stars, and more

Authors: Yveline Lebreton, Daniel Roy Reese

Abstract: Stellar parameters are required in a variety of contexts, ranging from the characterisation of exoplanets to Galactic archaeology. Among them, the age of stars cannot be directly measured, while the mass and radius can be measured in some particular cases (binary systems, interferometry). Stellar ages, masses, and radii have to be inferred from stellar evolution models by appropriate techniques. W… ▽ More Stellar parameters are required in a variety of contexts, ranging from the characterisation of exoplanets to Galactic archaeology. Among them, the age of stars cannot be directly measured, while the mass and radius can be measured in some particular cases (binary systems, interferometry). Stellar ages, masses, and radii have to be inferred from stellar evolution models by appropriate techniques. We have designed a Python tool named SPInS. It takes a set of photometric, spectroscopic, interferometric, and/or asteroseismic observational constraints and, relying on a stellar model grid, provides the age, mass, and radius of a star, among others, as well as error bars and correlations. We make the tool available to the community via a dedicated website. SPInS uses a Bayesian approach to find the PDF of stellar parameters from a set of classical constraints. At the heart of the code is a MCMC solver coupled with interpolation within a pre-computed stellar model grid. Priors can be considered, such as the IMF or SFR. SPInS can characterise single stars or coeval stars, such as members of binary systems or of stellar clusters. We illustrate the capabilities of SPInS by studying stars that are spread over the Hertzsprung-Russell diagram. We then validate the tool by inferring the ages and masses of stars in several catalogues and by comparing them with literature results. We show that in addition to the age and mass, SPInS can efficiently provide derived quantities, such as the radius, surface gravity, and seismic indices. We demonstrate that SPInS can age-date and characterise coeval stars that share a common age and chemical composition. The SPInS tool will be very helpful in preparing and interpreting the results of large-scale surveys, such as the wealth of data expected or already provided by space missions, such as Gaia, Kepler, TESS, and PLATO. △ Less

Submitted 31 August, 2020; originally announced September 2020.

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

Journal ref: A&A 642, A88 (2020)

Very regular high-frequency pulsation modes in young intermediate-mass stars

Authors: Timothy R. Bedding, Simon J. Murphy, Daniel R. Hey, Daniel Huber, Tanda Li, Barry Smalley, Dennis Stello, Timothy R. White, Warrick H. Ball, William J. Chaplin, Isabel L. Colman, Jim Fuller, Eric Gaidos, Daniel R. Harbeck, J. J. Hermes, Daniel L. Holdsworth, Gang Li, Yaguang Li, Andrew W. Mann, Daniel R. Reese, Sanjay Sekaran, Jie Yu, Victoria Antoci, Christoph Bergmann, Timothy M. Brown , et al. (11 additional authors not shown)

Abstract: Asteroseismology is a powerful tool for probing the internal structures of stars by using their natural pulsation frequencies. It relies on identifying sequences of pulsation modes that can be compared with theoretical models, which has been done successfully for many classes of pulsators, including low-mass solar-type stars, red giants, high-mass stars and white dwarfs. However, a large group of… ▽ More Asteroseismology is a powerful tool for probing the internal structures of stars by using their natural pulsation frequencies. It relies on identifying sequences of pulsation modes that can be compared with theoretical models, which has been done successfully for many classes of pulsators, including low-mass solar-type stars, red giants, high-mass stars and white dwarfs. However, a large group of pulsating stars of intermediate mass--the so-called delta Scuti stars--have rich pulsation spectra for which systematic mode identification has not hitherto been possible. This arises because only a seemingly random subset of possible modes are excited, and because rapid rotation tends to spoil the regular patterns. Here we report the detection of remarkably regular sequences of high-frequency pulsation modes in 60 intermediate-mass main-sequence stars, allowing definitive mode identification. Some of these stars have space motions that indicate they are members of known associations of young stars, and modelling of their pulsation spectra confirms that these stars are indeed young. △ Less

Submitted 13 May, 2020; originally announced May 2020.

Comments: published in Nature https://www.nature.com/articles/s41586-020-2226-8

The MUSTANG-2 Galactic Plane Survey (MGPS90) pilot

Authors: Adam Ginsburg, L. D. Anderson, Simon Dicker, Charles Romero, Brian Svoboda, Mark Devlin, Roberto Galván-Madrid, Remy Indebetouw, Hauyu Baobab Liu, Brian Mason, Tony Mroczkowski, W. P. Armentrout, John Bally, Crystal Brogan, Natalie Butterfield, Todd R. Hunter, Erik D. Reese, Erik Rosolowsky, Craig Sarazin, Yancy Shirley, Jonathan Sievers, Sara Stanchfield

Abstract: We report the results of a pilot program for a Green Bank Telescope (GBT) MUSTANG Galactic Plane survey at 3 mm (90 GHz), MGPS90. The survey achieves a typical $1σ$ depth of $1-2$ mJy beam$^{-1}$ with a 9" beam. We describe the survey parameters, quality assessment process, cataloging, and comparison with other data sets. We have identified 709 sources over seven observed fields selecting some of… ▽ More We report the results of a pilot program for a Green Bank Telescope (GBT) MUSTANG Galactic Plane survey at 3 mm (90 GHz), MGPS90. The survey achieves a typical $1σ$ depth of $1-2$ mJy beam$^{-1}$ with a 9" beam. We describe the survey parameters, quality assessment process, cataloging, and comparison with other data sets. We have identified 709 sources over seven observed fields selecting some of the most prominent millimeter-bright regions between $0°< \ell < 50°$ (total area $\approx 7.5 °^2$). The majority of these sources have counterparts at other wavelengths. By applying flux selection criteria to these sources, we successfully recovered several known hypercompact HII (HCHII) regions, but did not confirm any new ones. We identify 126 sources that have mm-wavelength counterparts but do not have cm-wavelength counterparts and are therefore candidate HCHII regions; of these, 10 are morphologically compact and are strong candidates for new HCHII regions. Given the limited number of candidates in the extended area in this survey compared to the relatively large numbers seen in protoclusters W51 and W49, it appears that most HCHII regions exist within dense protoclusters. Comparing the counts of HCHII to ultracompact HII (UCHII) regions, we infer the HCHII region lifetime is 16-46% that of the UCHII region lifetime. We additionally separated the 3 mm emission into dust and free-free emission by comparing with archival 870 $μ$m and 20 cm data. In the selected pilot fields, most ($\gtrsim80$%) of the 3 mm emission comes from plasma, either through free-free or synchrotron emission. △ Less

Submitted 22 April, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

Comments: Accepted to ApJS. Data publicly released, with links in paper. Reposted to fix missing figure labels

A realistic two-dimensional model of Altair

Authors: K. Bouchaud, A. Domiciano de Souza, M. Rieutord, D. R. Reese, P. Kervella

Abstract: Fast rotation is responsible for important changes in the structure and evolution of stars. Optical long baseline interferometry now permits the study of its effects on the stellar surface, mainly gravity darkening and flattening. We aim to determine the fundamental parameters of the fast-rotating star Altair, in particular its evolutionary stage, mass, and differential rotation, using state-of-th… ▽ More Fast rotation is responsible for important changes in the structure and evolution of stars. Optical long baseline interferometry now permits the study of its effects on the stellar surface, mainly gravity darkening and flattening. We aim to determine the fundamental parameters of the fast-rotating star Altair, in particular its evolutionary stage, mass, and differential rotation, using state-of-the-art stellar interior and atmosphere models together with interferometric, spectroscopic, and asteroseismic observations. We use ESTER 2D stellar models to produce the relevant surface parameters needed to create intensity maps from atmosphere models. Interferometric and spectroscopic observables are computed from these intensity maps and several stellar parameters are then adjusted using the MCMC algorithm Emcee. We determined Altair's equatorial radius to be 2.008 +/- 0.006 Rsun, the position angle 301.1 +/- 0.3 degrees, the inclination 50.7 +/- 1.2 degrees, and the equatorial angular velocity 0.74 +/- 0.01 times the Keplerian angular velocity. This angular velocity leads to a flattening of 0.220 +/- 0.003. We also deduce from the spectroscopically derived vsini ~ 243 km/s, a true equatorial velocity of ~314 km/s corresponding to a rotation period of 7h46m (~3 c/d). The data also impose a strong correlation between mass, metallicity, hydrogen abundance, and core evolution. Thanks to asteroseismic data, we constrain the mass of Altair to 1.86 +/- 0.03 Msun and further deduce its metallicity Z = 0.019 and its core hydrogen mass fraction Xc = 0.71, assuming an initial solar hydrogen mass fraction X = 0.739. These values suggest that Altair is ~100 Myrs old. Finally, the 2D ESTER model also gives the internal differential rotation of Altair, showing that its core rotates approximately 50% faster than the envelope, while the surface differential rotation does not exceed 6%. △ Less

Submitted 6 December, 2019; originally announced December 2019.

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

Journal ref: A&A 633, A78 (2020)

Chemical mixing in low mass stars: I. Rotation against atomic diffusion including radiative accelerations

Authors: M. Deal, M-J. Goupil, J. P. Marques, D. Reese, Y. Lebreton

Abstract: When modelling stars with masses larger than 1.2Msun with no observed chemical peculiarity, atomic diffusion is often neglected because, on its own, it causes unrealistic surface abundances compared with those observed. The reality is that atomic diffusion is in competition with other transport processes. The purpose of this study is to quantify the opposite or conjugated effects of atomic diffusi… ▽ More When modelling stars with masses larger than 1.2Msun with no observed chemical peculiarity, atomic diffusion is often neglected because, on its own, it causes unrealistic surface abundances compared with those observed. The reality is that atomic diffusion is in competition with other transport processes. The purpose of this study is to quantify the opposite or conjugated effects of atomic diffusion and rotationally induced mixing in stellar models of low mass stars. Our second goal is to estimate the impact of neglecting both rotational mixing and atomic diffusion in stellar parameter inferences for stars with masses larger than 1.3Msun. Using the AIMS code, we infer the masses and ages of a set of representative artificial stars for which models were computed with the CESTAM evolution code, taking into account rotationally induced mixing and atomic diffusion, including radiative accelerations. We show that for masses lower than 1.3Msun, rotation dominates the transport of chemical elements, and strongly reduces the effect of atomic diffusion, with net surface abundance modifications similar to solar ones. At larger mass, atomic diffusion and rotation are competing equally. Above 1.44Msun, atomic diffusion dominates in stellar models with initial rotation smaller than 80km.s-1 producing a chemical peculiarity which is not observed in Kepler-legacy stars. This indicates that a transport process of chemical elements is missing. Importantly, neglecting rotation and atomic diffusion (including radiative accelerations) in the models, when inferring the parameters of F-type stars, may lead to errors of 5%, 2.5% and 25% respectively for stellar masses, radii and ages. Atomic diffusion (including radiative accelerations) and rotational mixing should be taken into account in stellar models in order to determine accurate stellar parameters. △ Less

Submitted 31 October, 2019; originally announced October 2019.

Comments: 15 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 633, A23 (2020)

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

An ALMA+ACA measurement of the shock in the Bullet Cluster

Authors: Luca Di Mascolo, Tony Mroczkowski, Eugene Churazov, Maxim Markevitch, Kaustuv Basu, Tracy E. Clarke, Mark Devlin, Brian S. Mason, Scott W. Randall, Erik D. Reese, Rashid Sunyaev, Daniel R. Wik

Abstract: The thermal Sunyaev-Zeldovich (SZ) effect presents a relatively new tool for characterizing galaxy cluster merger shocks, traditionally studied through X-ray observations. Widely regarded as the "textbook example" of a cluster merger bow shock, the western shock front in the Bullet Cluster (1E0657-56) represents the ideal test case for such an SZ study. We aim to reconstruct a parametric model for… ▽ More The thermal Sunyaev-Zeldovich (SZ) effect presents a relatively new tool for characterizing galaxy cluster merger shocks, traditionally studied through X-ray observations. Widely regarded as the "textbook example" of a cluster merger bow shock, the western shock front in the Bullet Cluster (1E0657-56) represents the ideal test case for such an SZ study. We aim to reconstruct a parametric model for the shock SZ signal by directly and jointly fitting deep, high-resolution interferometric data from the Atacama Large Millimeter/submillimeter Array (ALMA) and Atacama Compact Array (ACA) in Fourier space. The ALMA+ACA data are primarily sensitive to the electron pressure difference across the shock front. To estimate the shock Mach number $M$, this difference can be combined with the value for the upstream electron pressure derived from an independent Chandra X-ray analysis. In the case of instantaneous electron-ion temperature equilibration, we find $M=2.08^{+0.12}_{-0.12}$, in $\approx 2.4σ$ tension with the independent constraint from Chandra, $M_X=2.74\pm0.25$. The assumption of purely adiabatic electron temperature change across the shock leads to $M=2.53^{+0.33}_{-0.25}$, in better agreement with the X-ray estimate $M_X=2.57\pm0.23$ derived for the same heating scenario. We have demonstrated that interferometric observations of the SZ effect provide constraints on the properties of the shock in the Bullet Cluster that are highly complementary to X-ray observations. The combination of X-ray and SZ data yields a powerful probe of the shock properties, capable of measuring $M$ and addressing the question of electron-ion equilibration in cluster shocks. Our analysis is however limited by systematics related to the overall cluster geometry and the complexity of the post-shock gas distribution. To overcome these limitations, a joint analysis of SZ and X-ray data is needed. △ Less

Submitted 30 October, 2019; v1 submitted 17 July, 2019; originally announced July 2019.

Comments: 9 pages (including appendices), 5 figures, and 1 table; accepted for publication in A&A

Journal ref: A&A 628, A100 (2019)

Revisiting the pulsational characteristics of the exoplanet host star $β$ Pictoris

Authors: K. Zwintz, D. R. Reese, C. Neiner, A. Pigulski, R. Kuschnig, M. Muellner, S. Zieba, L. Abe, T. Guillot, G. Handler, M. Kenworthy, R. Stuik, A. F. J. Moffat, A. Popowicz, S. M. Rucinski, G. A. Wade, W. W. Weiss, J. I. Bailey III, S. Crawford, M. Ireland, R. Kuhn, B. Lomberg, E. E. Mamajek, S. N. Mellon, G. J. Talens

Abstract: Exoplanet properties crucially depend on their host stars' parameters. In case the exoplanet host star shows pulsations, asteroseismology can be used for an improved description of the stellar parameters. We aim to revisit the pulsational properties of beta Pic and identify its pulsation modes from normalised amplitudes in five different passbands. We also investigate the potential presence of a m… ▽ More Exoplanet properties crucially depend on their host stars' parameters. In case the exoplanet host star shows pulsations, asteroseismology can be used for an improved description of the stellar parameters. We aim to revisit the pulsational properties of beta Pic and identify its pulsation modes from normalised amplitudes in five different passbands. We also investigate the potential presence of a magnetic field. We conduct a frequency analysis using three seasons of BRITE-Constellation observations in the BRITE blue and red filters, the ~620-day long bRing light curve and the nearly 8-year long SMEI photometric time series. We calculate normalised amplitudes using all passbands including previously published values obtained from ASTEP observations. We investigate the magnetic properties of beta Pic using spectropolarimetric observations conducted with the HARPSpol instrument. Using 2D rotating models, we fit the normalised amplitudes and frequencies through Monte Carlo Markov Chains. We identify 15 pulsation frequencies in the range from 34 to 55c/d, where two display clear amplitude variability. We use the normalised amplitudes in up to five passbands to identify the modes as three l = 1, six l = 2 and six l = 3 modes. beta Pic is shown to be non-magnetic with an upper limit of the possible undetected dipolar field of 300G. Multiple fits to the frequencies and normalised amplitudes are obtained including one with a near equator-on inclination for beta Pic, which corresponds to our expectations based on the orbital inclination of beta Pic b and the orientation of the circumstellar disk. This solution leads to a rotation rate of 27% of the Keplerian break-up velocity, a radius of 1.497+-0.025Rsun, and a mass of 1.797+-0.035Msun. The ~2% errors in radius and mass do not account for uncertainties in the models and a potentially erroneous mode-identification. △ Less

Submitted 29 May, 2019; originally announced May 2019.

Comments: 26 pages, 19 figures, accepted in Astronomy & Astrophysics

Journal ref: A&A 627, A28 (2019)

The PLATO Solar-like Light-curve Simulator: A tool to generate realistic stellar light-curves with instrumental effects representative of the PLATO mission

Authors: R. Samadi, A. Deru, D. Reese, V. Marchiori, E. Grolleau, J. J. Green, M. Pertenais, Y. Lebreton, S. Deheuvels, B. Mosser, K. Belkacem, A. Borner, A. M. S. Smith

Abstract: The preparation of science objectives of the ESA's PLATO space mission will require the implementation of hare-and-hound exercises relying on the massive generation of representative simulated light-curves. We developed a light-curve simulator named the PLATO Solar-like Light-curve Simulator (PSLS) in order to generate light-curves representative of typical PLATO targets, i.e. showing simultaneous… ▽ More The preparation of science objectives of the ESA's PLATO space mission will require the implementation of hare-and-hound exercises relying on the massive generation of representative simulated light-curves. We developed a light-curve simulator named the PLATO Solar-like Light-curve Simulator (PSLS) in order to generate light-curves representative of typical PLATO targets, i.e. showing simultaneously solar-like oscillations, stellar granulation, and magnetic activity. At the same time, PSLS also aims at mimicking in a realistic way the random noise and the systematic errors representative of the PLATO multi-telescope concept. To quantify the instrumental systematic errors, we performed a series of simulations at pixel level that include various relevant sources of perturbations expected for PLATO. From the simulated pixels, we extract the photometry as planned on-board. The simulated light-curves are then corrected for instrumental effects using the instrument Point Spread Functions reconstructed on the basis of a microscanning technique that will be operated during the in-flight calibration phases. These corrected light-curves are then fitted by a parametric model, which we incorporated in PSLS. We show that the instrumental systematic errors dominate the signal only at frequencies below 20muHz and are found to mainly depend on stellar magnitude and on the detector charge transfer inefficiency. To illustrate how realistic our simulator is, we compared its predictions with observations made by Kepler on three typical targets and found a good qualitative agreement with the observations. PSLS reproduces the main properties of expected PLATO light-curves. Its speed of execution and its inclusion of relevant stellar signals as well as sources of noises representative of the PLATO cameras make it an indispensable tool for the scientific preparation of the PLATO mission. △ Less

Submitted 23 March, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 19 pages, 16 figures, Astronomy & Astrophysics (in press)

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

AIMS - A new tool for stellar parameter determinations using asteroseismic constraints

Authors: B. M. Rendle, G. Buldgen, A. Miglio, D. Reese, A. Noels, G. R. Davies, T. L. Campante, W. J. Chaplin, M. N. Lund, J. S. Kuszlewicz, L. J. A. Scott, R. Scuflaire, W. H. Ball, J. Smetana, B. Nsamba

Abstract: A key aspect in the determination of stellar properties is the comparison of observational constraints with predictions from stellar models. Asteroseismic Inference on a Massive Scale (AIMS) is an open source code that uses Bayesian statistics and a Markov Chain Monte Carlo approach to find a representative set of models that reproduce a given set of classical and asteroseismic constraints. These… ▽ More A key aspect in the determination of stellar properties is the comparison of observational constraints with predictions from stellar models. Asteroseismic Inference on a Massive Scale (AIMS) is an open source code that uses Bayesian statistics and a Markov Chain Monte Carlo approach to find a representative set of models that reproduce a given set of classical and asteroseismic constraints. These models are obtained by interpolation on a pre-calculated grid, thereby increasing computational efficiency. We test the accuracy of the different operational modes within AIMS for grids of stellar models computed with the Liège stellar evolution code (main sequence and red giants) and compare the results to those from another asteroseismic analysis pipeline, PARAM. Moreover, using artificial inputs generated from models within the grid (assuming the models to be correct), we focus on the impact on the precision of the code when considering different combinations of observational constraints (individual mode frequencies, period spacings, parallaxes, photospheric constraints,...). Our tests show the absolute limitations of precision on parameter inferences using synthetic data with AIMS, and the consistency of the code with expected parameter uncertainty distributions. Interpolation testing highlights the significance of the underlying physics to the analysis performance of AIMS and provides caution as to the upper limits in parameter step size. All tests demonstrate the flexibility and capability of AIMS as an analysis tool and its potential to perform accurate ensemble analysis with current and future asteroseismic data yields. △ Less

Submitted 9 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS. 17 pages, 17 figures

Imaging the Thermal and Kinematic Sunyaev-Zel'dovich Effect Signals in a Sample of Ten Massive Galaxy Clusters: Constraints on Internal Velocity Structures and Bulk Velocities

Authors: Jack Sayers, Alfredo Montaña, Tony Mroczkowski, Grant W. Wilson, Michael Zemcov, Adi Zitrin, Nathália Cibirka, Sunil Golwala, David Hughes, Daisuke Nagai, Erik D. Reese, David Sánchez, John Zuhone

Abstract: We have imaged the Sunyaev-Zel'dovich (SZ) effect signals at 140 and 270 GHz towards ten galaxy clusters with Bolocam and AzTEC/ASTE. We also used Planck data to constrain the signal at large angular scales, Herschel-SPIRE images to subtract the brightest galaxies that comprise the cosmic infrared background (CIB), Chandra imaging to map the electron temperature $T_e$ of the intra-cluster medium (… ▽ More We have imaged the Sunyaev-Zel'dovich (SZ) effect signals at 140 and 270 GHz towards ten galaxy clusters with Bolocam and AzTEC/ASTE. We also used Planck data to constrain the signal at large angular scales, Herschel-SPIRE images to subtract the brightest galaxies that comprise the cosmic infrared background (CIB), Chandra imaging to map the electron temperature $T_e$ of the intra-cluster medium (ICM), and HST imaging to derive models of each galaxy cluster's mass density. The galaxy clusters gravitationally lens the background CIB, which produced an on-average reduction in brightness towards the galaxy clusters' centers after the brightest galaxies were subtracted. We corrected for this deficit, which was between 5-25% of the 270 GHz SZ effect signal within $R_{2500}$. Using the SZ effect measurements, along with the X-ray constraint on $T_e$, we measured each galaxy cluster's average line of sight (LOS) velocity $v_z$ within $R_{2500}$, with a median per-cluster uncertainty of +-700 km/s. We found an ensemble-mean <$v_z$> of 430+-210 km/s, and an intrinsic cluster-to-cluster scatter $σ_{int}$ of 470+-340 km/s. We also obtained maps of $v_z$ over each galaxy cluster's face with an angular resolution of 70". All four galaxy clusters previously identified as having a merger oriented along the LOS showed an excess variance in these maps at a significance of 2-4$σ$, indicating an internal $v_z$ rms of $\gtrsim$1000 km/s. None of the six galaxy clusters previously identified as relaxed or plane of sky mergers showed any such excess variance. △ Less

Submitted 25 June, 2019; v1 submitted 17 December, 2018; originally announced December 2018.

Comments: Accepted for publication in ApJ

Model physics in low-mass solar-type stars: atomic diffusion and metallicity mixture

Authors: Benard Nsamba, Tiago L. Campante, Mário J. P. F. G. Monteiro, Margarida S. Cunha, Ben M. Rendle, Daniel R. Reese, Kuldeep Verma

Abstract: Using asteroseismic data from the Kepler satellite, we explore the systematic uncertainties arising from changes in the input physics used when constructing evolution models of solar-type stars. We assess the impact of including atomic diffusion and of varying the metallicity mixture on the determination of global stellar parameters (i.e., radius, mass, and age). We find significant systematic unc… ▽ More Using asteroseismic data from the Kepler satellite, we explore the systematic uncertainties arising from changes in the input physics used when constructing evolution models of solar-type stars. We assess the impact of including atomic diffusion and of varying the metallicity mixture on the determination of global stellar parameters (i.e., radius, mass, and age). We find significant systematic uncertainties on global stellar parameters when diffusion is included in stellar grids. Furthermore, we find the systematic uncertainties on the global stellar parameters to be comparable to the statistical uncertainties when a different metallicity mixture is employed in stellar grids. △ Less

Submitted 2 December, 2018; originally announced December 2018.

Comments: 3 pages, 3 figures, proceeding of PHOST (Physics of Oscillating Stars) - a conference in honour of Prof. H. Shibahashi, 2-7 Sept. 2018, Banyuls-sur-mer (France)

Mode Classification in Fast-Rotating Stars using a Convolutional Neural Network: Model-Based Regular Patterns in $δ$ Scuti Stars

Authors: Giovanni M. Mirouh, George C. Angelou, Daniel R. Reese, Guglielmo Costa

Abstract: Oscillation modes in fast-rotating stars can be split into several subclasses, each with their own properties. To date, seismology of these stars cannot rely on regular pattern analysis and scaling relations. However, recently there has been the promising discovery of large separations observed in spectra of fast-rotating $δ$ Scuti stars: they were attributed to the island-mode subclass, and linke… ▽ More Oscillation modes in fast-rotating stars can be split into several subclasses, each with their own properties. To date, seismology of these stars cannot rely on regular pattern analysis and scaling relations. However, recently there has been the promising discovery of large separations observed in spectra of fast-rotating $δ$ Scuti stars: they were attributed to the island-mode subclass, and linked to the stellar mean density through a scaling law. In this work, we investigate the relevance of this scaling relation by computing models of fast-rotating stars and their oscillation spectra. In order to sort the thousands of oscillation modes thus obtained, we train a convolutional neural network isolating the island modes with 96\% accuracy. Arguing that the observed large separation is systematically smaller than the asymptotic one, we retrieve the observational $Δν- \overlineρ$ scaling law. This relation will be used to drive forward modelling efforts, and is a first step towards mode identification and inversions for fast-rotating stars. △ Less

Submitted 14 November, 2018; originally announced November 2018.

Comments: 6 pages, 6 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society Letters

Mean density inversions for red giants and red clump stars

Authors: Gaël Buldgen, B. Rendle, T. Sonoi, G. R. Davies, A. Miglio, S. J. A. J. Salmon, D. R. Reese, D. Bossini, P. Eggenberger, A. Noels, R. Scuflaire

Abstract: Since the CoRoT and Kepler missions, the availability of high quality seismic spectra for red giants has made them the standard clocks and rulers for Galactic Archeology. With the expected excellent data from the TESS and PLATO missions, red giants will again play a key role in Galactic studies and stellar physics, thanks to the precise masses and radii determined by asteroseismology. The determin… ▽ More Since the CoRoT and Kepler missions, the availability of high quality seismic spectra for red giants has made them the standard clocks and rulers for Galactic Archeology. With the expected excellent data from the TESS and PLATO missions, red giants will again play a key role in Galactic studies and stellar physics, thanks to the precise masses and radii determined by asteroseismology. The determination of these quantities is often based on so-called scaling laws, which have been used extensively for main-sequence stars. We show how the SOLA inversion technique can provide robust determinations of the mean density of red giants within 1 per cent of the real value, using only radial oscillations. Combined with radii determinations from Gaia of around 2 per cent precision, this approach provides robust, less model-dependent masses with an error lower than 10 per cent. It will improve age determinations, helping to accurately dissect the Galactic structure and history. We present results on artificial data of standard models, models including an extended atmosphere from averaged 3D simulations and non-adiabatic frequency calculations to test surface effects, and on eclipsing binaries. We show that the inversions provide very robust mean density estimates, using at best seismic information. However, we also show that a distinction between red-giant branch and red-clump stars is required to determine a reliable estimate of the mean density. The stability of the inversion enables an implementation in automated pipelines, making it suitable for large samples of stars. △ Less

Submitted 25 August, 2018; originally announced August 2018.

Comments: Accepted for publication in MNRAS

Asteroseismic modelling of solar-type stars: internal systematics from input physics and surface correction methods

Authors: B. Nsamba, T. L. Campante, M. J. P. F. G. Monteiro, M. S. Cunha, B. M. Rendle, D. R. Reese, K. Verma

Abstract: Asteroseismic forward modelling techniques are being used to determine fundamental properties (e.g. mass, radius, and age) of solar-type stars. The need to take into account all possible sources of error is of paramount importance towards a robust determination of stellar properties. We present a study of 34 solar-type stars for which high signal-to-noise asteroseismic data is available from multi… ▽ More Asteroseismic forward modelling techniques are being used to determine fundamental properties (e.g. mass, radius, and age) of solar-type stars. The need to take into account all possible sources of error is of paramount importance towards a robust determination of stellar properties. We present a study of 34 solar-type stars for which high signal-to-noise asteroseismic data is available from multi-year Kepler photometry. We explore the internal systematics on the stellar properties, that is, associated with the uncertainty in the input physics used to construct the stellar models. In particular, we explore the systematics arising from: (i) the inclusion of the diffusion of helium and heavy elements; and (ii) the uncertainty in solar metallicity mixture. We also assess the systematics arising from (iii) different surface correction methods used in optimisation/fitting procedures. The systematics arising from comparing results of models with and without diffusion are found to be 0.5%, 0.8%, 2.1%, and 16% in mean density, radius, mass, and age, respectively. The internal systematics in age are significantly larger than the statistical uncertainties. We find the internal systematics resulting from the uncertainty in solar metallicity mixture to be 0.7% in mean density, 0.5% in radius, 1.4% in mass, and 6.7% in age. The surface correction method by Sonoi et al. and Ball & Gizon's two-term correction produce the lowest internal systematics among the different correction methods, namely, ~1%, ~1%, ~2%, and ~8% in mean density, radius, mass, and age, respectively. Stellar masses obtained using the surface correction methods by Kjeldsen et al. and Ball & Gizon's one-term correction are systematically higher than those obtained using frequency ratios. △ Less

Submitted 13 April, 2018; originally announced April 2018.

Comments: 13 pages, 10 figures, 4 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society Journal

Mode identification in rapidly rotating stars from BRITE data

Authors: Daniel R. Reese, Marc-Antoine Dupret, Michel Rieutord

Abstract: Apart from recent progress in Gamma Dor stars, identifying modes in rapidly rotating stars is a formidable challenge due to the lack of simple, easily identifiable frequency patterns. As a result, it is necessary to look to observational methods for identifying modes. Two popular techniques are spectroscopic mode identification based on line profile variations (LPVs) and photometric mode identific… ▽ More Apart from recent progress in Gamma Dor stars, identifying modes in rapidly rotating stars is a formidable challenge due to the lack of simple, easily identifiable frequency patterns. As a result, it is necessary to look to observational methods for identifying modes. Two popular techniques are spectroscopic mode identification based on line profile variations (LPVs) and photometric mode identification based on amplitude ratios and phase differences between multiple photometric bands. In this respect, the BRITE constellation is particularly interesting as it provides space-based multi-colour photometry. The present contribution describes the latest developments in obtaining theoretical predictions for amplitude ratios and phase differences for pulsation modes in rapidly rotating stars. These developments are based on full 2D non-adiabatic pulsation calculations, using models from the ESTER code, the only code to treat in a self-consistent way the thermal equilibrium of rapidly rotating stars. These predictions are then specifically applied to the BRITE photometric bands to explore the prospects of identifying modes based on BRITE observations. △ Less

Submitted 20 March, 2018; originally announced March 2018.

Comments: 8 pages, 3 figures, proceedings of the 3rd BRITE Science Workshop

Non-adiabatic oscillations of fast-rotating stars: the example of Rasalhague

Authors: Giovanni M. Mirouh, Daniel R. Reese, Michel Rieutord, Jérôme Ballot

Abstract: Early-type stars generally tend to be fast rotators. In these stars, mode identification is very challenging as the effects of rotation are not well known. We consider here the example of $α$ Ophiuchi, for which dozens of oscillation frequencies have been measured. We model the star using the two-dimensional structure code ESTER, and we compute both adiabatic and non-adiabatic oscillations using t… ▽ More Early-type stars generally tend to be fast rotators. In these stars, mode identification is very challenging as the effects of rotation are not well known. We consider here the example of $α$ Ophiuchi, for which dozens of oscillation frequencies have been measured. We model the star using the two-dimensional structure code ESTER, and we compute both adiabatic and non-adiabatic oscillations using the TOP code. Both calculations yield very complex spectra, and we used various diagnostic tools to try and identify the observed pulsations. While we have not reached a satisfactory mode-to-mode identification, this paper presents promising early results. △ Less

Submitted 16 November, 2017; originally announced November 2017.

Comments: 4 pages, 3 figures. SF2A 2017 proceedings

Constraining convective regions with asteroseismic linear structural inversions

Authors: G. Buldgen, D. R. Reese, M. A. Dupret

Abstract: Context. Convective regions in stellar models are always associated with uncertainties, for example due to extra-mixing or the possible inaccurate position of the transition from convective to radiative transport of energy. These have a strong impact on stellar models and their fundamental parameters. The most promising method to reduce these uncertainties is to use asteroseismology to derive diag… ▽ More Context. Convective regions in stellar models are always associated with uncertainties, for example due to extra-mixing or the possible inaccurate position of the transition from convective to radiative transport of energy. These have a strong impact on stellar models and their fundamental parameters. The most promising method to reduce these uncertainties is to use asteroseismology to derive diagnostics probing the structural characteristics of these regions. Aims. We wish to use custom-made integrated quantities to improve the capabilities of seismology to probe convective regions in stellar interiors. We hope to increase the number of indicators obtained with structural seismic inversions to provide additional constraints on stellar models and the fundamental parameters determined from theoretical modeling. Methods. First, we present new kernels associated with a proxy of the entropy in stellar interiors. We then show how these kernels can be used to build custom-made integrated quantities probing convective regions inside stellar models. We present two indicators suited to probe convective cores and envelopes, respectively, and test them on artificial data. Results. We show that it is possible to probe both convective cores and envelopes using appropriate indicators obtained with structural inversion techniques. These indicators provide direct constraints on a proxy of the entropy of the stellar plasma, sensitive to the characteristics of convective regions. These constraints can then be used to improve the modeling of solar-like stars by providing an additional degree of selection of models obtained from classical forward modeling approaches. We also show that in order to obtain very accurate indicators, we need l = 3 modes for the envelope but that the core-conditions indicator is more flexible in terms of the seismic data required for its use. △ Less

Submitted 14 November, 2017; originally announced November 2017.

Comments: Accepted for publication in Astronomy and Astrophysics

Asteroseismic Stellar Modelling with AIMS

Authors: Mikkel N. Lund, Daniel R. Reese

Abstract: The goal of AIMS (Asteroseismic Inference on a Massive Scale) is to estimate stellar parameters and credible intervals/error bars in a Bayesian manner from a set of asteroseismic frequency data and so-called classical constraints. To achieve reliable parameter estimates and computational efficiency, it searches through a grid of pre-computed models using an MCMC algorithm -- interpolation within t… ▽ More The goal of AIMS (Asteroseismic Inference on a Massive Scale) is to estimate stellar parameters and credible intervals/error bars in a Bayesian manner from a set of asteroseismic frequency data and so-called classical constraints. To achieve reliable parameter estimates and computational efficiency, it searches through a grid of pre-computed models using an MCMC algorithm -- interpolation within the grid of models is performed by first tessellating the grid using a Delaunay triangulation and then doing a linear barycentric interpolation on matching simplexes. Inputs for the modelling consist of individual frequencies from peak-bagging, which can be complemented with classical spectroscopic constraints. AIMS is mostly written in Python with a modular structure to facilitate contributions from the community. Only a few computationally intensive parts have been rewritten in Fortran in order to speed up calculations. △ Less

Submitted 6 November, 2017; originally announced November 2017.

Comments: 11 pages, 4 figures. Tutorial presented at the IVth Azores International Advanced School in Space Sciences on "Asteroseismology and Exoplanets: Listening to the Stars and Searching for New Worlds" (arXiv:1709.00645), which took place in Horta, Azores Islands, Portugal in July 2016

Report number: AsteroseismologyExoplanets/2017/08

Stellar Inversion Techniques

Authors: Daniel R. Reese

Abstract: Stellar seismic inversions have proved to be a powerful technique for probing the internal structure of stars, and paving the way for a better understanding of the underlying physics by revealing some of the shortcomings in current stellar models. In this lecture, we provide an introduction to this topic by explaining kernel-based inversion techniques. Specifically, we explain how various kernels… ▽ More Stellar seismic inversions have proved to be a powerful technique for probing the internal structure of stars, and paving the way for a better understanding of the underlying physics by revealing some of the shortcomings in current stellar models. In this lecture, we provide an introduction to this topic by explaining kernel-based inversion techniques. Specifically, we explain how various kernels are obtained from the pulsation equations, and describe inversion techniques such as the Regularised Least-Squares (RLS) and Optimally Localised Averages (OLA) methods. △ Less

Submitted 6 November, 2017; originally announced November 2017.

Comments: 20 pages, 8 figures. Lecture presented at the IVth Azores International Advanced School in Space Sciences on "Asteroseismology and Exoplanets: Listening to the Stars and Searching for New Worlds" (arXiv:1709.00645), which took place in Horta, Azores Islands, Portugal in July 2016

Report number: AsteroseismologyExoplanets/2017/04

Non-adiabatic pulsations in ESTER models

Authors: Daniel Roy Reese, Marc-Antoine Dupret, Michel Rieutord

Abstract: One of the greatest challenges in interpreting the pulsations of rapidly rotating stars is mode identification, i.e. correctly matching theoretical modes to observed pulsation frequencies. Indeed, the latest observations as well as current theoretical results show the complexity of pulsation spectra in such stars, and the lack of easily recognisable patterns. In the present contribution, the lates… ▽ More One of the greatest challenges in interpreting the pulsations of rapidly rotating stars is mode identification, i.e. correctly matching theoretical modes to observed pulsation frequencies. Indeed, the latest observations as well as current theoretical results show the complexity of pulsation spectra in such stars, and the lack of easily recognisable patterns. In the present contribution, the latest results on non-adiabatic effects in such pulsations are described, and we show how these come into play when identifying modes. These calculations fully take into account the effects of rapid rotation, including centrifugal distortion, and are based on models from the ESTER project, currently the only rapidly rotating models in which the energy conservation equation is satisfied, a prerequisite for calculating non-adiabatic effects. Non-adiabatic effects determine which modes are excited and play a key role in the near-surface pulsation-induced temperature variations which intervene in multi-colour amplitude ratios and phase differences, as well as line profile variations. △ Less

Submitted 20 October, 2017; originally announced October 2017.

Comments: Proceedings for the Joint TASC2 & KASC9 Workshop, Terceira, Azores, 2016

Determining the metallicity of the solar envelope using seismic inversion techniques

Authors: Gaël Buldgen, S. J. A. J. Salmon, A. Noels, R. Scuflaire, M. A. Dupret, D. R. Reese

Abstract: The solar metallicity issue is a long-lasting problem of astrophysics, impacting multi- ple fields and still subject to debate and uncertainties. While spectroscopy has mostly been used to determine the solar heavy elements abundance, helioseismologists at- tempted providing a seismic determination of the metallicity in the solar convective enveloppe. However, the puzzle remains since two independ… ▽ More The solar metallicity issue is a long-lasting problem of astrophysics, impacting multi- ple fields and still subject to debate and uncertainties. While spectroscopy has mostly been used to determine the solar heavy elements abundance, helioseismologists at- tempted providing a seismic determination of the metallicity in the solar convective enveloppe. However, the puzzle remains since two independent groups prodived two radically different values for this crucial astrophysical parameter. We aim at provid- ing an independent seismic measurement of the solar metallicity in the convective enveloppe. Our main goal is to help provide new information to break the current stalemate amongst seismic determinations of the solar heavy element abundance. We start by presenting the kernels, the inversion technique and the target function of the inversion we have developed. We then test our approach in multiple hare-and-hounds exercises to assess its reliability and accuracy. We then apply our technique to solar data using calibrated solar models and determine an interval of seismic measurements for the solar metallicity. We show that our inversion can indeed be used to estimate the solar metallicity thanks to our hare-and-hounds exercises. However, we also show that further dependencies in the physical ingredients of solar models lead to a low accuracy. Nevertheless, using various physical ingredients for our solar models, we determine metallicity values between 0.008 and 0.014. △ Less

Submitted 1 September, 2017; originally announced September 2017.

Comments: Accepted for publication in MNRAS

Inversions of the Ledoux discriminant: a closer look at the tachocline

Authors: Gaël Buldgen, S. J. A. J. Salmon, M. Godart, A. Noels, R. Scuflaire, M. A. Dupret, D. R. Reese, J. Colgan, C. J. Fontes, P. Eggenberger, P. Hakel, D. P. Kilcrease, O. Richard

Abstract: Modelling the base of the solar convective envelope is a tedious problem. Since the first rotation inversions, solar modellers are confronted with the fact that a region of very limited extent has an enormous physical impact on the Sun. Indeed, it is the transition region from differential to solid body rotation, the tachocline, which furthermore is influenced by turbulence and is also supposed to… ▽ More Modelling the base of the solar convective envelope is a tedious problem. Since the first rotation inversions, solar modellers are confronted with the fact that a region of very limited extent has an enormous physical impact on the Sun. Indeed, it is the transition region from differential to solid body rotation, the tachocline, which furthermore is influenced by turbulence and is also supposed to be the seat of the solar magnetic dynamo. Moreover, solar models show significant disagreement with the sound speed profile in this region. In this paper, we show how helioseismology can provide further constraints on this region by carrying out an inversion of the Ledoux discriminant. We compare these inversions for Standard Solar Models built using various opacity tables and chemical abundances and discuss the origins of the discrepancies between Solar Models and the Sun. △ Less

Submitted 1 September, 2017; originally announced September 2017.

Comments: Accepted for publication as a letter in MNRAS

Precise surface gravities of $δ$ Scuti stars from asteroseismology

Authors: A. García Hernández, J. C. Suárez, A. Moya, Mário J. P. F. G. Monteiro, Z. Guo, D. R. Reese, J. Pascual-Granado, S. Barceló Forteza, S. Martín-Ruiz, R. Garrido, J. Nieto

Abstract: The work reported here demonstrates that it is possible to accurately determine surface gravities of $δ$ Sct stars using the frequency content from high precision photometry and a measurement of the parallax. Using a sample of 10 eclipsing binary systems with a $δ$ Sct component and the unique $δ$ Sct star discovered with a transiting planet, WASP-33, we were able to refine the $Δν-\barρ$ relation… ▽ More The work reported here demonstrates that it is possible to accurately determine surface gravities of $δ$ Sct stars using the frequency content from high precision photometry and a measurement of the parallax. Using a sample of 10 eclipsing binary systems with a $δ$ Sct component and the unique $δ$ Sct star discovered with a transiting planet, WASP-33, we were able to refine the $Δν-\barρ$ relation. Using this relation and parallaxes, we obtained independent values for the masses and radii, allowing us to calculate the surface gravities without any constraints from spectroscopic or binary analysis. A remarkably good agreement was found between our results and those published, extracted from the analysis of the radial velocities and light curves of the systems. This reinforces the potential of $Δν$ as a valuable observable for $δ$ Sct stars and settles the degeneracy problem for the $\log g$ determination through spectroscopy. △ Less

Submitted 21 July, 2017; originally announced July 2017.

Comments: 5 pages, 2 figures, 1 table

Seismic inversion of the solar entropy: A case for improving the Standard Solar Model

Authors: G. Buldgen, S. J. A. J. Salmon, A. Noels, R. Scuflaire, D. R. Reese, M-A. Dupret, J. Colgan, C. J. Fontes, P. Eggenberger, P. Hakel, D. P. Kilcrease, S. Turck-Chièze

Abstract: The Sun is the most constrained and well-studied of all stars. As a consequence, the physical ingredients entering solar models are used as a reference to study all other stars observed in the Universe. However, our understanding of the solar structure is still imperfect, as illustrated by the current debate on the heavy element abundances in the Sun. We wish to provide additional information on t… ▽ More The Sun is the most constrained and well-studied of all stars. As a consequence, the physical ingredients entering solar models are used as a reference to study all other stars observed in the Universe. However, our understanding of the solar structure is still imperfect, as illustrated by the current debate on the heavy element abundances in the Sun. We wish to provide additional information on the solar structure by carrying out structural inversions of a new physical quantity, a proxy of the entropy of the solar plasma which properties are very sensitive to the temperature gradient below the convective zone. We use new structural kernels to carry out direct inversions of an entropy proxy of the solar plasma and compare the solar structure to various standard solar models built using various opacity tables and chemical abundances. We also link our results to classical tests commonly found in the literature. Our analysis allows us to probe more efficiently the uncertain regions of the solar models, just below the convective zone, paving the way for new in-depth analyses of the Sun taking into account additional physical uncertainties of solar models beyond the specific question of chemical abundances. △ Less

Submitted 17 July, 2017; originally announced July 2017.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 607, A58 (2017)

Comparing Jupiter interior structure models to Juno gravity measurements and the role of a dilute core

Authors: S. M. Wahl, W. B. Hubbard, B. Militzer, T. Guillot, Y. Miguel, N. Movshovitz, Y. Kaspi, R. Helled, D. Reese, E. Galanti, S. Levin, J. E. Connerney, S. J. Bolton

Abstract: The Juno spacecraft has measured Jupiter's low-order, even gravitational moments, $J_2$--$J_8$, to an unprecedented precision, providing important constraints on the density profile and core mass of the planet. Here we report on a selection of interior models based on ab initio computer simulations of hydrogen-helium mixtures. We demonstrate that a dilute core, expanded to a significant fraction o… ▽ More The Juno spacecraft has measured Jupiter's low-order, even gravitational moments, $J_2$--$J_8$, to an unprecedented precision, providing important constraints on the density profile and core mass of the planet. Here we report on a selection of interior models based on ab initio computer simulations of hydrogen-helium mixtures. We demonstrate that a dilute core, expanded to a significant fraction of the planet's radius, is helpful in reconciling the calculated $J_n$ with Juno's observations. Although model predictions are strongly affected by the chosen equation of state, the prediction of an enrichment of $Z$ in the deep, metallic envelope over that in the shallow, molecular envelope holds. We estimate Jupiter's core to contain an 7--25 Earth mass of heavy elements. We discuss the current difficulties in reconciling measured $J_n$ with the equations of state, and with theory for formation and evolution of the planet. △ Less

Submitted 6 July, 2017; originally announced July 2017.

Comments: 21 Pages, 3 figures, plus Supplementary Information

Journal ref: Geophysical Research Letters, 44, 4649-4659 (2017)

What CoRoT tells us about δ Scuti stars existence of a regular pattern and seismic indices to characterize stars

Authors: Eric Michel, Marc-Antoine Dupret, Daniel Reese, Rhita-Maria Ouazzani, Jonas Debosscher, Antonio García Hernández, Kevin Belkacem, Reza Samadi, Sébastien Salmon, Juan Carlos Suarez, Sebastia Barceló Forteza

Abstract: Inspired by the so appealing example of red giants, where going from a handful of stars to thousands revealed the structure of the eigenspectrum, we inspected a large homogeneous set of around 1860 δ Scuti stars observed with CoRoT. This unique data set reveals a common regular pattern which appears to be in agreement with island modes featured by theoretical non-perturbative treatments of fast ro… ▽ More Inspired by the so appealing example of red giants, where going from a handful of stars to thousands revealed the structure of the eigenspectrum, we inspected a large homogeneous set of around 1860 δ Scuti stars observed with CoRoT. This unique data set reveals a common regular pattern which appears to be in agreement with island modes featured by theoretical non-perturbative treatments of fast rotation. The comparison of these data with models and linear stability calculations suggests that spectra can be fruitfully characterized to first order by a few parameters which might play the role of seismic indices for δ Scuti stars, as {Δν} and {ν_{max}} do for red giants. The existence of this pattern offers an observational support for guiding further theoretical works on fast rotation. It also provides a framework for further investigation of the observational material collected by CoRoT and Kepler. Finally, it sketches out the perspective of using δ Scuti stars pulsations for ensemble asteroseismology. △ Less

Submitted 10 May, 2017; originally announced May 2017.

Comments: 7 pages, 9 figures, proc. of the Conference Seismology of the Sun and the Distant Stars 2016 Using Today s Successes to Prepare the Future, Joint TASC2 and KASC9 Workshop , SPACEINN and HELAS8 Conference, 11-15 July 2016, Angra do Heroísmo, Terceira-Açores, Portugal

Linear approximation of seismic inversions: new kernels and structural effects

Authors: Gaël Buldgen, Daniel Reese, Marc-Antoine Dupret

Abstract: Thanks to the space-based photometry missions CoRoT and Kepler, we now benefit from a wealth of seismic data for stars other than the sun. In the future, K2, Tess and Plato will provide further observations. The quality of this data may allow kernel-based linear structural inversion techniques to be used for stars other than the sun. To understand the limitations of this approach, we analyse the v… ▽ More Thanks to the space-based photometry missions CoRoT and Kepler, we now benefit from a wealth of seismic data for stars other than the sun. In the future, K2, Tess and Plato will provide further observations. The quality of this data may allow kernel-based linear structural inversion techniques to be used for stars other than the sun. To understand the limitations of this approach, we analyse the validity of the linear assumption used in these inversion techniques. We inspect various structural pairs and see how they are affected by structural changes. We show that uncertainties in radius strongly affect structural pairs of nondimensional variables, and that various other effects might come into play. Amongst these, the importance of micro-physics give the most striking example of how uncertainties in stellar models impact the verification of the linear relations. We also point out that including seismic constraints in the forward modelling fit helps with satisfying the linear relations. △ Less

Submitted 23 February, 2017; originally announced February 2017.

Comments: To appear in the proceedings of the Astrofluid 2016 Conference

Asteroseismic inversions in the Kepler era: application to the Kepler Legacy sample

Authors: Gaël Buldgen, Daniel Reese, Marc-Antoine Dupret

Abstract: In the past few years, the CoRoT and Kepler missions have carried out what is now called the space photometry revolution. This revolution is still ongoing thanks to K2 and will be continued by the Tess and Plato2.0 missions. However, the photometry revolution must also be followed by progress in stellar modelling, in order to lead to more precise and accurate determinations of fundamental stellar… ▽ More In the past few years, the CoRoT and Kepler missions have carried out what is now called the space photometry revolution. This revolution is still ongoing thanks to K2 and will be continued by the Tess and Plato2.0 missions. However, the photometry revolution must also be followed by progress in stellar modelling, in order to lead to more precise and accurate determinations of fundamental stellar parameters such as masses, radii and ages. In this context, the long-lasting problems related to mixing processes in stellar interior is the main obstacle to further improvements of stellar modelling. In this contribution, we will apply structural asteroseismic inversion techniques to targets from the Kepler Legacy sample and analyse how these can help us constrain the fundamental parameters and mixing processes in these stars. Our approach is based on previous studies using the SOLA inversion technique to determine integrated quantities such as the mean density, the acoustic radius, and core conditions indicators, and has already been successfully applied to the 16Cyg binary system. We will show how this technique can be applied to the Kepler Legacy sample and how new indicators can help us to further constrain the chemical composition profiles of stars as well as provide stringent constraints on stellar ages. △ Less

Submitted 24 February, 2017; v1 submitted 23 February, 2017; originally announced February 2017.

Comments: To appear in the proceedings of the Kasc 9 Tasc 2 workshop

Frequency regularities of acoustic modes and multi-colour mode identification in rapidly rotating stars

Authors: D. R. Reese, F. Lignières, J. Ballot, M. -A. Dupret, C. Barban, C. van 't Veer-Menneret, K. B. MacGregor

Abstract: Context: Mode identification has remained a major obstacle in the interpretation of pulsation spectra in rapidly rotating stars. Aims: We would like to test mode identification methods and seismic diagnostics in rapidly rotating stars, using oscillation spectra based on new theoretical predictions. Methods: We investigate the auto-correlation function and Fourier transform of theoretically cal… ▽ More Context: Mode identification has remained a major obstacle in the interpretation of pulsation spectra in rapidly rotating stars. Aims: We would like to test mode identification methods and seismic diagnostics in rapidly rotating stars, using oscillation spectra based on new theoretical predictions. Methods: We investigate the auto-correlation function and Fourier transform of theoretically calculated frequency spectra, in which modes are selected according to their visibilities. Given the difficulties in predicting intrinsic mode amplitudes, we experimented with various ad-hoc prescriptions for setting these, including using random values. Furthermore, we analyse the ratios between mode amplitudes observed in different photometric bands. Results: When non-random intrinsic mode amplitudes are used, our results show that it is possible to extract the large frequency separation or half its value, and sometimes twice the rotation rate, from the auto-correlation function. The Fourier transforms are mostly sensitive to the large frequency separation or half its value. When the intrinsic mode amplitudes include random factors, the results are far less favourable. We also find that amplitude ratios provide a good way of grouping together modes with similar characteristics. By analysing the frequencies of these groups, it is possible to constrain mode identification as well as determine the large frequency separation and the rotation rate. △ Less

Submitted 31 January, 2017; originally announced January 2017.

Comments: 15 pages, 12 figures, accepted in A&A

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

Standing on the shoulders of Dwarfs: the Kepler asteroseismic LEGACY sample II - radii, masses, and ages

Authors: Victor Silva Aguirre, Mikkel N. Lund, H. M. Antia, Warrick H. Ball, Sarbani Basu, Jørgen Christensen-Dalsgaard, Yveline Lebreton, Daniel R. Reese, Kuldeep Verma, Luca Casagrande, Anders B. Justesen, Jakob R. Mosumgaard, William J. Chaplin, Timothy R. Bedding, Guy R. Davies, Rasmus Handberg, Günter Houdek, Daniel Huber, Hans Kjeldsen, David W. Latham, Timothy R. White, Hugo R. Coelho, Andrea Miglio, Ben Rendle

Abstract: We use asteroseismic data from the Kepler satellite to determine fundamental stellar properties of the 66 main-sequence targets observed for at least one full year by the mission. We distributed tens of individual oscillation frequencies extracted from the time series of each star among seven modelling teams who applied different methods to determine radii, masses, and ages for all stars in the sa… ▽ More We use asteroseismic data from the Kepler satellite to determine fundamental stellar properties of the 66 main-sequence targets observed for at least one full year by the mission. We distributed tens of individual oscillation frequencies extracted from the time series of each star among seven modelling teams who applied different methods to determine radii, masses, and ages for all stars in the sample. Comparisons among the different results reveal a good level of agreement in all stellar properties, which is remarkable considering the variety of codes, input physics and analysis methods employed by the different teams. Average uncertainties are of the order of $\sim$2\% in radius, $\sim$4\% in mass, and $\sim$10\% in age, making this the best-characterised sample of main-sequence stars available to date. Our predicted initial abundances and mixing-length parameters are checked against inferences from chemical enrichment laws $ΔY / ΔZ$ and predictions from 3D atmospheric simulations. We test the accuracy of the determined stellar properties by comparing them to the Sun, angular diameter measurements, Gaia parallaxes, and binary evolution, finding excellent agreement in all cases and further confirming the robustness of asteroseismically-determined physical parameters of stars when individual frequencies of oscillation are available. Baptised as the Kepler dwarfs LEGACY sample, these stars are the solar-like oscillators with the best asteroseismic properties available for at least another decade. All data used in this analysis and the resulting stellar parameters are made publicly available for the community. △ Less

Submitted 2 December, 2016; v1 submitted 26 November, 2016; originally announced November 2016.

Comments: 23 pages, 14 figures, ApJ, accepted

Asteroseismic modelling of the Binary HD 176465

Authors: B. Nsamba, M. J. P. F. G. Monteiro, T. L. Campante, D. R. Reese, T. R. White, A. García Hernández, C. Jiang

Abstract: The detection and analysis of oscillations in binary star systems is critical in understanding stellar structure and evolution. This is partly because such systems have the same initial chemical composition and age. Solar-like oscillations have been detected by Kepler in both components of the asteroseismic binary HD 176465. We present an independent modelling of each star in this binary system. S… ▽ More The detection and analysis of oscillations in binary star systems is critical in understanding stellar structure and evolution. This is partly because such systems have the same initial chemical composition and age. Solar-like oscillations have been detected by Kepler in both components of the asteroseismic binary HD 176465. We present an independent modelling of each star in this binary system. Stellar models generated using MESA (Modules for Experiments in Stellar Astrophysics) were fitted to both the observed individual frequencies and complementary spectroscopic parameters. The individual theoretical oscillation frequencies for the corresponding stellar models were obtained using GYRE as the pulsation code. A Bayesian approach was applied to find the probability distribution functions of the stellar parameters using AIMS (Asteroseismic Inference on a Massive Scale) as the optimisation code. The ages of HD 176465 A and HD 176465 B were found to be 2.81 $\pm$ 0.48 Gyr and 2.52 $\pm$ 0.80 Gyr, respectively. These results are in agreement when compared to previous studies carried out using other asteroseismic modelling techniques and gyrochronology. △ Less

Submitted 17 November, 2016; originally announced November 2016.

Journal ref: EPJ Web of Conferences, Proceedings of the joint TASC2/KASC9 workshop - SPACEINN & HELAS8 conference 2016

Analysis of the linear approximation of seismic inversions for various structural pairs

Authors: G. Buldgen, D. R. Reese, M. A. Dupret

Abstract: Thanks to space-based photometry missions CoRoT and Kepler, we benefit from a wealth of seismic data for stars other than the sun. In the future, K2, Tess, and Plato will complement this data and provide observations in addition to those already at hand. The availability of this data leads to questions on how it is feasible to extend linear structural inversion techniques to stars other than the s… ▽ More Thanks to space-based photometry missions CoRoT and Kepler, we benefit from a wealth of seismic data for stars other than the sun. In the future, K2, Tess, and Plato will complement this data and provide observations in addition to those already at hand. The availability of this data leads to questions on how it is feasible to extend linear structural inversion techniques to stars other than the sun. Linked to this problem is the question of the validity of the linear assumption. In this study, we analyse its limitations with respect to changes of structural variables.We wish to provide a more extended theoretical background to structural linear inversions by doing a study of the validity of the linear assumption for various structural variables. First, we recall the origins of the linear assumption for structural stellar inversions and explain its importance for asteroseismic studies. We recall the impact of unknown structural quantities such as the mass and the radius of the star on structural inversion results. We explain how kernels for new structural variables can be derived using two methods, one suited to asteroseismic targets, the other to helioseismic targets. For this second method, we present a new structural pair, namely the (A, Y) structural kernels. The kernels are tested in various numerical experiments that enable us to evaluate the weaknesses of different pairs and the domains of validity of their respective linear regime. The numerical tests we carry out allow us to disentangle the impact of various uncertainties in stellar models on the verification of the linear integral relations. We show the importance of metallicity, the equation of state, extra-mixing, and inaccuracies in the microphysics. We discuss the limitations of the method of conjugated functions due to the lack of extremely precise determinations of masses or radii in the asteroseismology. △ Less

Submitted 5 October, 2016; originally announced October 2016.

Comments: 14 pages. Accepted for publication in Astronomy and Astrophysics

In-depth study of 16CygB using inversion techniques

Authors: G. Buldgen, S. J. A. S. Salmon, D. R. Reese, M-A. Dupret

Abstract: The 16Cyg binary system hosts the solar-like Kepler targets with the most stringent observational constraints. Moreover, this system is particularly interesting since both stars are very similar in mass but the A component is orbited by a red dwarf, whereas the B component is orbited by a Jovian planet and thus could have formed a more complex planetary system. In our previous study, we showed tha… ▽ More The 16Cyg binary system hosts the solar-like Kepler targets with the most stringent observational constraints. Moreover, this system is particularly interesting since both stars are very similar in mass but the A component is orbited by a red dwarf, whereas the B component is orbited by a Jovian planet and thus could have formed a more complex planetary system. In our previous study, we showed that seismic inversions of integrated quantities could be used to constrain microscopic diffusion in the A component. In this study, we analyse the B component in the light of a more regularised inversion. We wish to analyse independently the B component of the 16Cyg binary system using the inversion of an indicator dedicated to analyse core conditions, denoted tu. Using this independent determination, we wish to analyse any differences between both stars due to the potential influence of planetary formation on stellar structure and/or their respective evolution. First, we recall the observational constraints for 16CygB and the method we used to generate reference stellar models. The inversion results were then used to analyse the differences between the A and B components. The tu indicator for 16CygB shows a disagreement with models including microscopic diffusion and sharing the chemical composition previously derived for 16CygA. Small changes in chemical composition are insufficient to solve the problem but that extra mixing can account for the differences seen between both stars. We use a parametric approach to analyse the impact of extra mixing in the form of turbulent diffusion on the behaviour of the tu values. We conclude on the necessity of further investigations using models with a physically motivated implementation of extra mixing including additional constraints to further improve the accuracy with which the fundamental parameters of this system are determined. △ Less

Submitted 19 August, 2016; originally announced August 2016.

Comments: Accepted for publication in Astronomy & Astrophysics

SpaceInn hare-and-hounds exercise: Estimation of stellar properties using space-based asteroseismic data

Authors: D. R. Reese, W. J. Chaplin, G. R. Davies, A. Miglio, H. M. Antia, W. H. Ball, S. Basu, G. Buldgen, J. Christensen-Dalsgaard, H. R. Coelho, S. Hekker, G. Houdek, Y. Lebreton, A. Mazumdar, T. S. Metcalfe, V. Silva Aguirre, D. Stello, K. Verma

Abstract: Context: Detailed oscillation spectra comprising individual frequencies for numerous solar-type stars and red giants are or will become available. These data can lead to a precise characterisation of stars. Aims: Our goal is to test and compare different methods for obtaining stellar properties from oscillation frequencies and spectroscopic constraints, in order to evaluate their accuracy and th… ▽ More Context: Detailed oscillation spectra comprising individual frequencies for numerous solar-type stars and red giants are or will become available. These data can lead to a precise characterisation of stars. Aims: Our goal is to test and compare different methods for obtaining stellar properties from oscillation frequencies and spectroscopic constraints, in order to evaluate their accuracy and the reliability of the error bars. Methods: In the context of the SpaceInn network, we carried out a hare-and-hounds exercise in which one group produced "observed" oscillation spectra for 10 artificial solar-type stars, and various groups characterised these stars using either forward modelling or acoustic glitch signatures. Results: Results based on the forward modelling approach were accurate to 1.5 % (radius), 3.9 % (mass), 23 % (age), 1.5 % (surface gravity), and 1.8 % (mean density). For the two 1 Msun stellar targets, the accuracy on the age is better than 10 % thereby satisfying PLATO 2.0 requirements. The average accuracies for the acoustic radii of the base of the convection zone, the He II ionisation, and the Gamma_1 peak were 17 %, 2.4 %, and 1.9 %, respectively. Glitch fitting analysis seemed to be affected by aliasing problems for some of the targets. Conclusions: Forward modelling is the most accurate approach, but needs to be complemented by model-independent results from, e.g., glitch analysis. Furthermore, global optimisation algorithms provide more robust error bars. △ Less

Submitted 28 April, 2016; originally announced April 2016.

Comments: 29 pages, 13 figures, 20 tables, accepted in A&A

Journal ref: A&A 592, A14 (2016)