Showing 1–50 of 50 results for author: Bellinger, E P

Formation of Sub-Chandrasekhar Mass Black Holes and Red Stragglers via Hawking Stars in Ultra-Faint Dwarf Galaxies

Authors: Andrew D. Santarelli, Matthew E. Caplan, Earl P. Bellinger

Abstract: Presently, primordial black holes (PBHs) in the asteroid-mass window from $10^{-16}$ M$_\odot$ to $10^{-10}$ M$_\odot$ are a popular dark matter candidate. If they exist, some stars would capture them upon formation, and they would slowly accrete the star over gigayears. Such Hawking stars -- stars with a central PBH -- provide a novel channel for the formation of both sub-Chandrasekhar mass black… ▽ More Presently, primordial black holes (PBHs) in the asteroid-mass window from $10^{-16}$ M$_\odot$ to $10^{-10}$ M$_\odot$ are a popular dark matter candidate. If they exist, some stars would capture them upon formation, and they would slowly accrete the star over gigayears. Such Hawking stars -- stars with a central PBH -- provide a novel channel for the formation of both sub-Chandrasekhar mass black holes as well as red straggler stars. Here we report on stellar evolution models that extend our previous work to Hawking stars with masses between 0.5 and 1.4 M$_\odot$. We explore three accretion schemes, and find that a wide range of PBHs in the asteroid-mass window can robustly accrete stars as small as 1 M$_{\odot}$ within the age of the Universe. This mechanism of producing sub-solar mass black holes is highly dependent on the assumed accretion physics and stellar metallicity. Lower-metallicity stars are generally accreted more rapidly, suggesting that it may be more likely for sub-Chandrasekhar mass Hawking stars formed in the early universe, such as those in ultra-faint dwarf (UFD) galaxies, to transmute their star into a sub-Chandrasekhar mass black hole within a Hubble time. We present a stellar population synthesis of a Draco II-like UFD galaxy containing Hawking stars and show that the number of red stragglers they produce can qualitatively match the observed population. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 11 pages, 5 figures, to be submitted to the ApJ

An Expanded Set of Los Alamos OPLIB Tables in MESA: Type-1 Rosseland-mean Opacities and Solar Models

Authors: Ebraheem Farag, Christopher J. Fontes, F. X. Timmes, Earl P. Bellinger, Joyce A. Guzik, Evan B. Bauer, Suzannah R. Wood, Katie Mussack, Peter Hakel, James Colgan, David P. Kilcrease, Manolo E. Sherrill, Tryston C. Raecke, Morgan T. Chidester

Abstract: We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range to $Z$\,=\,10… ▽ More We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range to $Z$\,=\,10$^{-6}$ to $Z$\,=\,0.2 which allows improved accuracy of opacities at low and high metallicity, increases the table density in the metallicity range $Z$\,=\,10$^{-4}$ to $Z$\,=\,0.1 to enhance the accuracy of opacities drawn from interpolations across neighboring metallicities, and adds entries for hydrogen mass fractions between $X$\,=\,0 and $X$\,=\,0.1 including $X$\,=\,$10^{-2}, 10^{-3}, 10^{-4}, 10^{-5}, 10^{-6}$ that can improve stellar models of hydrogen deficient stars. We implement these new OPLIB radiative opacity tables in \MESA, and find that calibrated solar models agree broadly with previously published helioseismic and solar neutrino results. We find differences between using the new 1194 OPLIB opacity tables and the 126 OPAL opacity tables range from $\approx$\,20--80\% across individual chemical mixtures, up to $\approx$\,8\% and $\approx$\,15\% at the bottom and top of the solar convection zone respectively, and $\approx$\,7\% in the solar core. We also find differences between standard solar models using different opacity table sources that are on par with altering the initial abundance mixture. We conclude that this new, open-access set of OPLIB opacity tables does not solve the solar modeling problem, and suggest the investigation of physical mechanisms other than the atomic radiative opacity. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 28 pages, 12 (13) figures. Accepted for Publication in ApJ

The K2 Asteroseismic KEYSTONE sample of Dwarf and Subgiant Solar-Like Oscillators. I: Data and Asteroseismic parameters

Authors: Mikkel N. Lund, Sarbani Basu, Allyson Bieryla, Luca Casagrande, Daniel Huber, Saskia Hekker, Lucas Viani, Guy R. Davies, Tiago L. Campante, William J. Chaplin, Aldo M. Serenelli, J. M. Joel Ong, Warrick H. Ball, Amalie Stokholm, Earl P. Bellinger, Michaël Bazot, Dennis Stello, David W. Latham, Timothy R. White, Maryum Sayeed, Víctor Aguirre Børsen-Koch, Ashley Chontos

Abstract: The KEYSTONE project aims to enhance our understanding of solar-like oscillators by delivering a catalogue of global asteroseismic parameters (${Δν}$ and ${ν_{\rm max}}$) for 173 stars, comprising mainly dwarfs and subgiants, observed by the K2 mission in its short-cadence mode during campaigns 6-19. We derive atmospheric parameters and luminosities using spectroscopic data from TRES, astrometric… ▽ More The KEYSTONE project aims to enhance our understanding of solar-like oscillators by delivering a catalogue of global asteroseismic parameters (${Δν}$ and ${ν_{\rm max}}$) for 173 stars, comprising mainly dwarfs and subgiants, observed by the K2 mission in its short-cadence mode during campaigns 6-19. We derive atmospheric parameters and luminosities using spectroscopic data from TRES, astrometric data from $\textit{Gaia}$, and the infrared flux method (IRFM) for a comprehensive stellar characterisation. Asteroseismic parameters are robustly extracted using three independent methods, complemented by an iterative refinement of the spectroscopic analyses using seismic ${\log g}$ values to enhance parameter accuracy. Our analysis identifies new detections of solar-like oscillations in 159 stars, providing an important complement to already published results from previous campaigns. The catalogue provides homogeneously derived atmospheric parameters and luminosities for the majority of the sample. Comparison between spectroscopic ${T_{\rm eff}}$ and those obtained from the IRFM demonstrates excellent agreement. The iterative approach to spectroscopic analysis significantly enhances the accuracy of the stellar properties derived. △ Less

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

Comments: 22 pages, 18 figures, accepted for publication in Astronomy & Astrophysics (A&A)

Bridging theory and observations in stellar pulsations: The impact of convection and metallicity on the instability strips of Classical and Type-II Cepheids

Authors: Mami Deka, Earl P. Bellinger, Shashi M. Kanbur, Sukanta Deb, Anupam Bhardwaj, Hugh Riley Randall, Selim Kalici, Susmita Das

Abstract: The effect of metallicity on the theoretical and empirical period-luminosity (PL) relations of Cepheid variables is not well understood and remains a highly debated issue. Here, we examine empirical colour-magnitude diagrams (CMDs) of Classical and Type-II Cepheids in the Magellanic Clouds and compare those with the theoretically predicted instability strip (IS) edges. We explore the effects of in… ▽ More The effect of metallicity on the theoretical and empirical period-luminosity (PL) relations of Cepheid variables is not well understood and remains a highly debated issue. Here, we examine empirical colour-magnitude diagrams (CMDs) of Classical and Type-II Cepheids in the Magellanic Clouds and compare those with the theoretically predicted instability strip (IS) edges. We explore the effects of incorporating turbulent flux, turbulent pressure, and radiative cooling into the convection theory on the predicted IS at various metallicities using MESA-RSP. We find that the edges become redder with the increasing complexity of convection physics incorporated in the fiducial convection sets, and are similarly shifted to the red with increasing metallicity. The inclusion of turbulent flux and pressure improves the agreement of the red edge of the IS, while their exclusion leads to better agreement with observations of the blue edge. About 90% of observed stars are found to fall within the predicted bluest and reddest edges across the considered variations of turbulent convection parameters. Furthermore, we identify and discuss discrepancies between theoretical and observed CMDs in the low effective temperature and high luminosity regions for stars with periods greater than ~ 20 days. These findings highlight the potential for calibrating the turbulent convection parameters in stellar pulsation models or the prediction of a new class of rare, long-period, 'red Cepheids', thereby improving our understanding of Cepheids and their role in cosmological studies. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS on 25-04-2024. 19 pages, 15 figures, 8 tables

Stupendously Large Black Hole Coalescence and Hubble Tensions

Authors: Earl Patrick Bellinger, Jakob Stegmann, Tom Wagg

Abstract: Local and distant measurements of the Hubble constant are in significant tension: local measurements of the Hubble constant appear to show a Universe that is significantly contracted when compared to distant measurements. From the point of view of an observer, a passing gravitational wave could cause the Universe to appear locally contracted and expanded in a quadrupolar pattern. The inspiral of a… ▽ More Local and distant measurements of the Hubble constant are in significant tension: local measurements of the Hubble constant appear to show a Universe that is significantly contracted when compared to distant measurements. From the point of view of an observer, a passing gravitational wave could cause the Universe to appear locally contracted and expanded in a quadrupolar pattern. The inspiral of a pair of stupendously large black holes (SLABs) with a chirp mass of $\sim10^{21}~\rm{M}_\odot$ may produce gravitational radiation with sufficiently large amplitude and wavelength to increase $H_0$ in one direction, and decrease it by the same amount in the other. These tensions would then oscillate with a period corresponding to half the orbital period of the binary. If such a gravitational wave were aligned with the plane of the Milky Way, most readily visible galaxies would appear closer than they actually are, thereby causing the apparent Hubble tension. Due to the long binary period, we would be in the same phase of the gravitational wave for the complete history of astronomical observation. The opposite tension would be visible in the orthogonal directions, thus giving the opportunity to falsify the existence of inspiraling SLAB binaries. As a corollary, the Hubble tension may place an upper limit on the maximum mass of inspiraling black holes in the Universe. △ Less

Submitted 28 March, 2024; originally announced March 2024.

Comments: 4 pages, 3 figures

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

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

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

Submitted 8 March, 2024; originally announced March 2024.

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

Asteroseismic Inversions for Internal Sound Speed Profiles of Main-sequence Stars with Radiative Cores

Authors: Lynn Buchele, Earl P. Bellinger, Saskia Hekker, Sarbani Basu, Warrick Ball, Jørgen Christensen-Dalsgaard

Abstract: The theoretical oscillation frequencies of even the best asteroseismic models of solar-like oscillators show significant differences from observed oscillation frequencies. Structure inversions seek to use these frequency differences to infer the underlying differences in stellar structure. While used extensively to study the Sun, structure inversion results for other stars have so far been limited… ▽ More The theoretical oscillation frequencies of even the best asteroseismic models of solar-like oscillators show significant differences from observed oscillation frequencies. Structure inversions seek to use these frequency differences to infer the underlying differences in stellar structure. While used extensively to study the Sun, structure inversion results for other stars have so far been limited. Applying sound-speed inversions to more stars allows us to probe stellar theory over a larger range of conditions, as well as look for overall patterns that may hint at deficits in our current understanding. To that end, we present structure inversion results for 12 main-sequence solar-type stars with masses between 1M$_\odot$ and 1.15M$_\odot$. Our inversions are able to infer differences in the isothermal sound speed in the innermost 30% by radius of our target stars. In half of our target stars, the structure of our best-fit model fully agrees with the observations. In the remainder, the inversions reveal significant differences between the sound-speed profile of the star and that of the model. We find five stars where the sound speed in the core of our stellar models is too low and one star showing the opposite behavior. For the two stars which our inversions reveal the most significant differences, we examine whether changing the microphysics of our models improves them and find that changes to nuclear reaction rates or core opacities can reduce, but do not fully resolve, the differences. △ Less

Submitted 27 December, 2023; v1 submitted 19 December, 2023; originally announced December 2023.

Comments: 24 pages, 15 figures, accepted for publication in The Astrophysical Journal

Journal ref: ApJ 961 (2024) 198

Is there a black hole in the center of the Sun?

Authors: Matthew E. Caplan, Earl P. Bellinger, Andrew D. Santarelli

Abstract: There is probably not a black hole in the center of the sun. Despite this detail, our goal in this work to convince the reader that this question is interesting and that work studying stars with central black holes is well motivated. If primordial black holes exist then they may exist in sufficiently large numbers to explain the dark matter in the universe. While primordial black holes may form at… ▽ More There is probably not a black hole in the center of the sun. Despite this detail, our goal in this work to convince the reader that this question is interesting and that work studying stars with central black holes is well motivated. If primordial black holes exist then they may exist in sufficiently large numbers to explain the dark matter in the universe. While primordial black holes may form at almost any mass, the asteroid-mass window between $10^{-16} - 10^{-10} ~ \textrm{M}_\odot$ remains a viable dark matter candidate and these black holes could be captured by stars upon formation. Such a star, partially powered by accretion luminosity from a microscopic black hole in its core, has been called a `Hawking star.' Stellar evolution of Hawking stars is highly nontrivial and requires detailed stellar evolution models, which were developed in our recent work. We present here full evolutionary models of solar mass Hawking stars using two accretion schemes: one with a constant radiative efficiency, and one that is new in this work that uses an adaptive radiative efficiency to model the effects of photon trapping. △ Less

Submitted 19 January, 2024; v1 submitted 12 December, 2023; originally announced December 2023.

Comments: 12 pages, 4 figures, Invited article in Ap&SS

Solar evolution models with a central black hole

Authors: Earl P. Bellinger, Matt E. Caplan, Taeho Ryu, Deepika Bollimpalli, Warrick H. Ball, Florian Kühnel, R. Farmer, S. E. de Mink, Jørgen Christensen-Dalsgaard

Abstract: Hawking (1971) proposed that the Sun may harbor a primordial black hole whose accretion supplies some of the solar luminosity. Such an object would have formed within the first 1 s after the Big Bang with the mass of a moon or an asteroid. These light black holes are a candidate solution to the dark matter problem, and could grow to become stellar-mass black holes (BHs) if captured by stars. Here… ▽ More Hawking (1971) proposed that the Sun may harbor a primordial black hole whose accretion supplies some of the solar luminosity. Such an object would have formed within the first 1 s after the Big Bang with the mass of a moon or an asteroid. These light black holes are a candidate solution to the dark matter problem, and could grow to become stellar-mass black holes (BHs) if captured by stars. Here we compute the evolution of stars having such a BH at their center. We find that such objects can be surprisingly long-lived, with the lightest black holes having no influence over stellar evolution, while more massive ones consume the star over time to produce a range of observable consequences. Models of the Sun born about a BH whose mass has since grown to approximately $10^{-6}~\rm{M}_\odot$ are compatible with current observations. In this scenario, the Sun would first dim to half its current luminosity over a span of 100 Myr as the accretion starts to generate enough energy to quench nuclear reactions. The Sun would then expand into a fully-convective star, where it would shine luminously for potentially several Gyr with an enriched surface helium abundance, first as a sub-subgiant star, and later as a red straggler, before becoming a sub-solar-mass BH. We also present results for a range of stellar masses and metallicities. The unique internal structures of stars harboring BHs may make it possible for asteroseismology to discover them, should they exist. We conclude with a list of open problems and predictions. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 20 pages, 6 figures, accepted for publication in The Astrophysical Journal

Geometry of the LMC based on multi-phase analysis of multi-wavelength Cepheid light curves using OGLE-IV and Gaia DR3 data

Authors: Gautam Bhuyan, Sukanta Deb, Shashi M. Kanbur, Earl P. Bellinger, Mami Deka, Anupam Bhardwaj

Abstract: The period-luminosity (PL) relation of Cepheids in the Large Magellanic Cloud (LMC) plays a pivotal role in extra-galactic distance measurement and the determination of the Hubble constant $(H_{0})$. In this work, we probe the geometry of the LMC through a detailed study of multi-phase PL relations of these Cepheids, leveraging data from the OGLE-IV and Gaia DR3 databases. We analyse the light cur… ▽ More The period-luminosity (PL) relation of Cepheids in the Large Magellanic Cloud (LMC) plays a pivotal role in extra-galactic distance measurement and the determination of the Hubble constant $(H_{0})$. In this work, we probe the geometry of the LMC through a detailed study of multi-phase PL relations of these Cepheids, leveraging data from the OGLE-IV and Gaia DR3 databases. We analyse the light curves of a combined sample of $\sim$3300 fundamental (FU) and first overtone (FO) mode classical Cepheids. We obtain multi-phase data with $50$ phase points over a complete pulsation cycle from the OGLE $(V, I)$ and Gaia $(G,G_{\rm BP}, G_{\rm RP})$ photometric bands. We determine the distance modulus and reddening values of individual Cepheids by fitting a simultaneous reddening law to the apparent distance modulus values. We calculate the LMC viewing angle parameters: the inclination angle $(i)$ and position angle of line of nodes $(θ_{\rm lon})$ by fitting a plane of the form $z = f(x,y)$ to the three-dimensional distribution of Cepheids in Cartesian coordinates $(x,y,z)$. The values of LMC viewing angles from multi-phase PL relations are found to be: $i = 22\rlap{.}^{\circ}87 \pm 0\rlap{.}^{\circ}43 ~\textrm{(stat.)} \pm 0\rlap{.}^{\circ}53 ~\textrm{(syst.)}$, $θ_{\rm lon} = 154\rlap{.}^{\circ}76 \pm 1\rlap{.}^{\circ}16 ~\textrm{(stat.)} \pm 1\rlap{.}^{\circ}01 ~\textrm{(syst.)}$, respectively. The use of multi-phase PL relations in multiple bands results in lower uncertainties for the LMC viewing angle parameters as compared to those derived from the mean light PL relations. This shows that the use of multi-phase PL relations with multi-wavelength photometry significantly improves the precision of these measurements, allowing better constraints on the morphology and the structure of the LMC. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 16 pages, 11 figures

The Potential of Asteroseismology to Resolve the Blue Supergiant Problem

Authors: Earl Patrick Bellinger, Selma E. de Mink, Walter E. van Rossem, Stephen Justham

Abstract: Despite major progress in our understanding of massive stars, concerning discrepancies still remain between observations and theory. Most notable are the numerous stars observed beyond the theoretical main sequence, an evolutionary phase expected to be short-lived and hence sparsely populated. This is the "Blue Supergiant Problem." Stellar models with abnormal internal structures can provide long-… ▽ More Despite major progress in our understanding of massive stars, concerning discrepancies still remain between observations and theory. Most notable are the numerous stars observed beyond the theoretical main sequence, an evolutionary phase expected to be short-lived and hence sparsely populated. This is the "Blue Supergiant Problem." Stellar models with abnormal internal structures can provide long-lived solutions for this problem: core hydrogen-burning stars with oversized cores may explain the hotter ones, and core helium-burning stars with undersized cores may explain the cooler ones. Such stars may result from enhanced or suppressed mixing in single stars, or, more likely, as the products of binary interaction and stellar mergers. Here we investigate the potential of asteroseismology to uncover the nature of blue supergiants. We construct stellar models for the above scenarios and show that they predict $g$-mode period spacings that differ by an order of magnitude: $\sim$ 200~min versus $\sim$ 20~min for long-lived core H and He burning stars, respectively. For the classical scenario of H-shell-burning stars rapidly crossing the Hertzsprung gap, we furthermore predict changes of the order of $10^{-2}\,μ\rm{Hz\,yr}^{-1}$ in high-frequency modes; this effect would be in principle observable from $\sim 5$~yr of asteroseismic monitoring if these modes can be identified. This raises the possibility of revealing the internal structure of blue supergiants, and thus determining whether these stars are indeed binary merger products. These asteroseismic diagnostics may be measurable through long time-series observations from the ongoing TESS mission and upcoming PLATO mission, thereby laying a path toward resolving the blue supergiant problem. △ Less

Submitted 27 May, 2024; v1 submitted 31 October, 2023; originally announced November 2023.

Comments: Accepted for publication in the Astrophysical Journal Letters

Variability of Blue Supergiants in the LMC with TESS

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

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

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

Comments: Accepted to ApJ

On the feasibility of structure inversions for gravity-mode pulsators

Authors: Vincent Vanlaer, Conny Aerts, Earl P. Bellinger, Jørgen Christensen-Dalsgaard

Abstract: Gravity-mode asteroseismology has significantly improved our understanding of mixing in intermediate mass stars. However, theoretical pulsation periods of stellar models remain in tension with observations, and it is often unclear how the models of these stars should be further improved. Inversions provide a path forward by directly probing the internal structure of these stars from their pulsatio… ▽ More Gravity-mode asteroseismology has significantly improved our understanding of mixing in intermediate mass stars. However, theoretical pulsation periods of stellar models remain in tension with observations, and it is often unclear how the models of these stars should be further improved. Inversions provide a path forward by directly probing the internal structure of these stars from their pulsation periods, quantifying which parts of the model are in need of improvement. This method has been used for solar-like pulsators, but has not yet been applied to main-sequence gravity-mode pulsators. Our aim is to determine whether structure inversions for gravity-mode pulsators are feasible. We focus on the case of slowly rotating SPB stars. We computed and analyzed dipole mode kernels for three variables pairs: $(ρ,c), (N^2,c)$, and $(N^2,ρ)$. We assessed the potential of these kernels by predicting the oscillation frequencies of a model after perturbing its structure. We then tested two inversion methods, RLS and SOLA, using a model grid computed with MESA and GYRE. We find that changing the stellar structure affects the oscillation frequencies in a nonlinear way. The oscillation modes for which this nonlinear dependency is the strongest are in resonance with the near-core peak in the buoyancy frequency. The near core region of the star can be probed with SOLA, while RLS requires fine tuning to obtain accurate results. Both RLS and SOLA are strongly affected by the nonlinear dependencies on the structure differences, as these methods are based on a first-order approximation. These inversion methods need to be modified for meaningful applications of inversions to SPB stars. Our results show that inversions of gravity-mode pulsators are possible in principle, but that the typical inversion methods developed for solar-like oscillators are not applicable. [abridged] △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: Accepted to A&A, 21 pages, 33 figures

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

Predicting Physical Parameters of Cepheid and RR Lyrae variables in an Instant with Machine Learning

Authors: Anupam Bhardwaj, Earl P. Bellinger, Shashi M. Kanbur, Marcella Marconi

Abstract: We present a machine learning method to estimate the physical parameters of classical pulsating stars such as RR Lyrae and Cepheid variables based on an automated comparison of their theoretical and observed light curve parameters at multiple wavelengths. We train artificial neural networks (ANNs) on theoretical pulsation models to predict the fundamental parameters (mass, radius, luminosity, and… ▽ More We present a machine learning method to estimate the physical parameters of classical pulsating stars such as RR Lyrae and Cepheid variables based on an automated comparison of their theoretical and observed light curve parameters at multiple wavelengths. We train artificial neural networks (ANNs) on theoretical pulsation models to predict the fundamental parameters (mass, radius, luminosity, and effective temperature) of Cepheid and RR Lyrae stars based on their period and light curve parameters. The fundamental parameters of these stars can be estimated up to 60 percent more accurately when the light curve parameters are taken into consideration. This method was applied to the observations of hundreds of Cepheids and thousands of RR Lyrae in the Magellanic Clouds to produce catalogs of estimated masses, radii, luminosities, and other parameters of these stars. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: Proceedings of IAU GA Symposium - Machine Learning in Astronomy: Possibilities and Pitfalls, in Busan, South Korea, 2022

Improved asteroseismic inversions for red-giant surface rotation rates

Authors: F. Ahlborn, E. P. Bellinger, S. Hekker, S. Basu, D. Mokrytska

Abstract: Asteroseismic observations of internal stellar rotation have indicated a substantial lack of angular momentum transport in theoretical models of subgiant and red-giant stars. Accurate core and surface rotation rate measurements are therefore needed to constrain internal transport processes included in the models. We eliminate substantial systematic errors of asteroseismic surface rotation rates fo… ▽ More Asteroseismic observations of internal stellar rotation have indicated a substantial lack of angular momentum transport in theoretical models of subgiant and red-giant stars. Accurate core and surface rotation rate measurements are therefore needed to constrain internal transport processes included in the models. We eliminate substantial systematic errors of asteroseismic surface rotation rates found in previous studies. We propose a new objective function for the Optimally Localized Averages method of rotational inversions for red-giant stars, which results in more accurate envelope rotation rate estimates obtained from the same data. We use synthetic observations from stellar models across a range of evolutionary stages and masses to demonstrate the improvement. We find that our new inversion technique allows us to obtain estimates of the surface rotation rate that are independent of the core rotation. For a star at the base of the red-giant branch, we reduce the systematic error from about 20% to a value close to 0, assuming constant envelope rotation. We also show the equivalence between this method and the method of linearised rotational splittings. Our new rotational inversion method substantially reduces the systematic errors of red-giant surface rotation rates. In combination with independent measures of the surface rotation rate, this will allow better constraints to be set on the internal rotation profile. This will be a very important probe to further constrain the internal angular momentum transport along the lower part of the red-giant branch. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 17 pages, 22 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 668, A98 (2022)

Period-Colour and Amplitude-Colour relations for OGLE-$δ$ Scuti stars in the Galactic Bulge and LMC

Authors: Mami Deka, Shashi M. Kanbur, Sukanta Deb, Susmita Das, Kerdaris Kurbah, Earl P. Bellinger, Anupam Bhardwaj

Abstract: We present an analysis on the behaviour of the Galactic bulge and the Large Magellanic Cloud (LMC) $δ$ Scuti stars in terms of period-colour and amplitude-colour (PCAC) relations at maximum, mean and minimum light. The publicly available Optical Gravitational Lensing Experiment-IV (OGLE-IV) light curves for Galactic bulge and OGLE-III light curves for LMC $δ$ Scuti stars are exploited for the anal… ▽ More We present an analysis on the behaviour of the Galactic bulge and the Large Magellanic Cloud (LMC) $δ$ Scuti stars in terms of period-colour and amplitude-colour (PCAC) relations at maximum, mean and minimum light. The publicly available Optical Gravitational Lensing Experiment-IV (OGLE-IV) light curves for Galactic bulge and OGLE-III light curves for LMC $δ$ Scuti stars are exploited for the analysis. It has been found that the Galactic bulge $δ$ Scuti stars obey flat PC relations at maximum/mean/minimum light while the LMC $δ$ Scutis have sloped/sloped/flat PC relations at maximum/mean/minimum light. Both the Galactic bulge and the LMC $δ$ Scutis have sloped/flat/sloped AC relations at maximum/mean/minimum. These relations also show that Galactic $δ$ Scutis are hotter as compared to their LMC counterparts. The period-amplitude (PA) relations for $δ$ Scutis exhibit different behaviour in the Galactic bulge and the LMC. The LMC variables are found to have higher amplitudes at a given period. The amplitude of the Galactic bulge $δ$ Scuti shows a bimodal distribution which can be modelled using a two-component Gaussian Mixture Model: one component with a lower amplitude and another with a higher amplitude. The observed behaviour of the $δ$ Scuti PCAC relations can be explained using the theory of the interaction of hydrogen ionization front (HIF) and stellar photosphere as well as the PA diagram. We use MESA-RSP to calculate theoretical non-linear hydrodynamical pulsation models for $δ$ Scuti stars with input metallicities of $Z=0.02$ and $Z=0.008$ appropriate for the Galactic bulge and LMC, respectively. The observed PCAC relations and theoretical calculations support the HIF-photosphere interaction theory. △ Less

Submitted 3 September, 2022; v1 submitted 26 August, 2022; originally announced August 2022.

Comments: 12 pages, 8 figures, accepted in MNRAS

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

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

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

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

Comments: 50 pages, 29 figures; Accepted to ApJS

Algebraic and machine learning approach to hierarchical triple-star stability

Authors: Pavan Vynatheya, Adrian S. Hamers, Rosemary A. Mardling, Earl P. Bellinger

Abstract: We present two approaches to determine the dynamical stability of a hierarchical triple-star system. The first is an improvement on the Mardling-Aarseth stability formula from 2001, where we introduce a dependence on inner orbital eccentricity and improve the dependence on mutual orbital inclination. The second involves a machine learning approach, where we use a multilayer perceptron (MLP) to cla… ▽ More We present two approaches to determine the dynamical stability of a hierarchical triple-star system. The first is an improvement on the Mardling-Aarseth stability formula from 2001, where we introduce a dependence on inner orbital eccentricity and improve the dependence on mutual orbital inclination. The second involves a machine learning approach, where we use a multilayer perceptron (MLP) to classify triple-star systems as `stable' and `unstable'. To achieve this, we generate a large training data set of 10^6 hierarchical triples using the N-body code MSTAR. Both our approaches perform better than previous stability criteria, with the MLP model performing the best. The improved stability formula and the machine learning model have overall classification accuracies of 93 % and 95 % respectively. Our MLP model, which accurately predicts the stability of any hierarchical triple-star system within the parameter ranges studied with almost no computation required, is publicly available on Github in the form of an easy-to-use Python script. △ Less

Submitted 6 September, 2022; v1 submitted 7 July, 2022; originally announced July 2022.

Comments: Accepted for publication by MNRAS

Towards solar measurements of nuclear reaction rates

Authors: Earl Patrick Bellinger, Jørgen Christensen-Dalsgaard

Abstract: Nuclear reaction rates are a fundamental yet uncertain ingredient in stellar evolution models. The astrophysical S-factor pertaining to the initial reaction in the proton-proton chain is uncertain at the 1% level, which contributes a systematic but generally unpropagated error of similar order in the theoretical ages of stars. In this work, we study the prospect of improving the measurement of thi… ▽ More Nuclear reaction rates are a fundamental yet uncertain ingredient in stellar evolution models. The astrophysical S-factor pertaining to the initial reaction in the proton-proton chain is uncertain at the 1% level, which contributes a systematic but generally unpropagated error of similar order in the theoretical ages of stars. In this work, we study the prospect of improving the measurement of this and other reaction rates in the pp chain and CNO cycle using helioseismology and solar neutrinos. We show that when other aspects of the solar model are improved, then it shall be possible using current solar data to improve the precision of this measurement by nearly an order of magnitude, and hence the corresponding uncertainty on the ages of low-mass stars by a similar amount. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 8 pages, 11 figures, accepted for publication in MNRAS

Cooling Delays from Iron Sedimentation and Iron Inner Cores in White Dwarfs

Authors: M. E. Caplan, I. F. Freeman, C. J. Horowitz, A. Cumming, E. P. Bellinger

Abstract: Do white dwarfs have inner cores made of iron? Neutron rich nuclei like $^{56}$Fe experience a net gravitational force and sediment toward the core. Using new phase diagrams and molecular dynamics simulations, we show that $^{56}$Fe should separate into mesoscopic Fe-rich crystallites due to its large charge relative to the background. At solar abundances, these crystallites rapidly precipitate an… ▽ More Do white dwarfs have inner cores made of iron? Neutron rich nuclei like $^{56}$Fe experience a net gravitational force and sediment toward the core. Using new phase diagrams and molecular dynamics simulations, we show that $^{56}$Fe should separate into mesoscopic Fe-rich crystallites due to its large charge relative to the background. At solar abundances, these crystallites rapidly precipitate and form an inner core of order 100 km and $10^{-3} M_\odot$ that may be detectable with asteroseismology. Associated cooling delays could be up to a Gyr for low mass white dwarfs but are only $\sim$0.1 Gyr for massive white dwarfs, so while this mechanism may contribute to the Q-branch the heating is insufficient to fully explain it. △ Less

Submitted 25 August, 2021; originally announced August 2021.

Comments: 8 pages, 5 figures, accepted for publication in ApJ Letters

Age-Dating Red Giant Stars Associated with Galactic Disk and Halo Substructures

Authors: Samuel K. Grunblatt, Joel C. Zinn, Adrian M. Price-Whelan, Ruth Angus, Nicholas Saunders, Marc Hon, Amalie Stokholm, Earl P. Bellinger, Sarah L. Martell, Benoit Mosser, Emily Cunningham, Jamie Tayar, Daniel Huber, Jakob Lysgaard Rørsted, Victor Silva Aguirre

Abstract: The vast majority of Milky Way stellar halo stars were likely accreted from a small number ($\lesssim$3) of relatively large dwarf galaxy accretion events. However, the timing of these events is poorly constrained, relying predominantly on indirect dynamical mixing arguments or imprecise age measurements of stars associated with debris structures. Here, we aim to infer robust stellar ages for star… ▽ More The vast majority of Milky Way stellar halo stars were likely accreted from a small number ($\lesssim$3) of relatively large dwarf galaxy accretion events. However, the timing of these events is poorly constrained, relying predominantly on indirect dynamical mixing arguments or imprecise age measurements of stars associated with debris structures. Here, we aim to infer robust stellar ages for stars associated with galactic substructures to more directly constrain the merger history of the Galaxy. By combining kinematic, asteroseismic, and spectroscopic data where available, we infer stellar ages for a sample of 10 red giant stars that were kinematically selected to be associated with the stellar halo, a subset of which are associated with the Gaia-Enceladus-Sausage halo substructure, and compare their ages to 3 red giant stars in the Galactic disk. Despite systematic differences in both absolute and relative ages determined by this work, age rankings of stars in this sample are robust. Passing the same observable inputs to multiple stellar age determination packages, we measure a weighted average age for the Gaia-Enceladus-Sausage stars in our sample of 8 $\pm$ 3 (stat.) $\pm$ 1 (sys.) Gyr. We also determine hierarchical ages for the populations of Gaia-Enceladus-Sausage, in situ halo and disk stars, finding a Gaia-Enceladus-Sausage population age of 8.0$^{+3.2}_{-2.3}$ Gyr. Although we cannot distinguish hierarchical population ages of halo or disk structures with our limited data and sample of stars, this framework should allow distinct characterization of Galactic substructures using larger stellar samples and additional data available in the near future. △ Less

Submitted 10 September, 2021; v1 submitted 21 May, 2021; originally announced May 2021.

Comments: 28 pages, 21 figures. Published by AAS Journals

Journal ref: ApJ 916 (2021) Issue 2, id.88, 19 pp

Asteroseismic Inference of the Central Structure in a Subgiant Star

Authors: Earl P. Bellinger, Sarbani Basu, Saskia Hekker, Jørgen Chrisensen-Dalsgaard, Warrick H. Ball

Abstract: Asteroseismic measurements enable inferences of the underlying stellar structure, such as the density and the speed of sound at various points within the interior of the star. This provides an opportunity to test stellar evolution theory by assessing whether the predicted structure of a star agrees with the measured structure. Thus far, this kind of inverse analysis has only been applied to the Su… ▽ More Asteroseismic measurements enable inferences of the underlying stellar structure, such as the density and the speed of sound at various points within the interior of the star. This provides an opportunity to test stellar evolution theory by assessing whether the predicted structure of a star agrees with the measured structure. Thus far, this kind of inverse analysis has only been applied to the Sun and three solar-like main-sequence stars. Here we extend the technique to stars on the subgiant branch, and apply it to one of the best-characterized subgiants of the Kepler mission, HR 7322. The observation of mixed oscillation modes in this star facilitates inferences of the conditions of its inert helium core, nuclear-burning hydrogen shell, and the deeper parts of its radiative envelope. We find that despite significant differences in the mode frequencies, the structure near to the center of this star does not differ significantly from the predicted structure. △ Less

Submitted 29 July, 2021; v1 submitted 10 May, 2021; originally announced May 2021.

Comments: 12 pages, 15 figures. Accepted for publication in the Astrophysical Journal

Inverse analysis of asteroseismic data: a review

Authors: Earl P. Bellinger, Sarbani Basu, Saskia Hekker

Abstract: Asteroseismology has emerged as the best way to characterize the global and internal properties of nearby stars. Often, this characterization is achieved by fitting stellar evolution models to asteroseismic observations. The star under investigation is then assumed to have the properties of the best-fitting model, such as its age. However, the models do not fit the observations perfectly. This is… ▽ More Asteroseismology has emerged as the best way to characterize the global and internal properties of nearby stars. Often, this characterization is achieved by fitting stellar evolution models to asteroseismic observations. The star under investigation is then assumed to have the properties of the best-fitting model, such as its age. However, the models do not fit the observations perfectly. This is due to incorrect or missing physics in stellar evolution calculations, resulting in predicted stellar structures that are discrepant with reality. Through an inverse analysis of the asteroseismic data, it is possible to go further than fitting stellar models, and instead infer details about the actual internal structure of the star at some locations in its interior. Comparing theoretical and observed stellar structures then enables the determination of the locations where the stellar models have discrepant structure, and illuminates a path for improvements to our understanding of stellar evolution. In this invited review, we describe the methods of asteroseismic inversions, and outline the progress that is being made towards measuring the interiors of stars. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: 12 pages, 1 figure. Invited review, Dynamics of the Sun and Stars

Journal ref: Astrophysics and Space Science Proceedings, vol. 57 (2020) 171-183

TESS observations of Cepheid stars: first light results

Authors: E. Plachy, A. Pál, A. Bódi, P. Szabó, L. Molnár, L. Szabados, J. M. Benkő, R. I. Anderson, E. P. Bellinger, A. Bhardwaj, M. Ebadi, K. Gazeas, F. -J. Hambsch, A. Hasanzadeh, M. I. Jurkovic, M. J. Kalaee, P. Kervella, K. Kolenberg, P. Mikołajczyk, N. Nardetto, J. M. Nemec, H. Netzel, C. -C. Ngeow, D. Ozuyar, J. Pascual-Granado , et al. (11 additional authors not shown)

Abstract: We present the first analysis of Cepheid stars observed by the TESS space mission in Sectors 1 to 5. Our sample consists of 25 pulsators: ten fundamental mode, three overtone and two double-mode classical Cepheids, plus three Type II and seven anomalous Cepheids. The targets were chosen from fields with different stellar densities, both from the Galactic field and from the Magellanic System. Three… ▽ More We present the first analysis of Cepheid stars observed by the TESS space mission in Sectors 1 to 5. Our sample consists of 25 pulsators: ten fundamental mode, three overtone and two double-mode classical Cepheids, plus three Type II and seven anomalous Cepheids. The targets were chosen from fields with different stellar densities, both from the Galactic field and from the Magellanic System. Three targets have 2-minute cadence light curves available by the TESS Science Processing Operations Center: for the rest, we prepared custom light curves from the full-frame images with our own differential photometric FITSH pipeline. Our main goal was to explore the potential and the limitations of TESS concerning the various subtypes of Cepheids. We detected many low amplitude features: weak modulation, period jitter, and timing variations due to light-time effect. We also report signs of non-radial modes and the first discovery of such a mode in an anomalous Cepheid, the overtone star XZ Cet, which we then confirmed with ground-based multicolor photometric measurements. We prepared a custom photometric solution to minimize saturation effects in the bright fundamental-mode classical Cepheid, $β$ Dor with the lightkurve software, and we revealed strong evidence of cycle-to-cycle variations in the star. In several cases, however, fluctuations in the pulsation could not be distinguished from instrumental effects, such as contamination from nearby sources which also varies between sectors. Finally, we discuss how precise light curve shapes will be crucial not only for classification purposes but also to determine physical properties of these stars. △ Less

Submitted 17 December, 2020; originally announced December 2020.

Comments: 30 pages, 18 figures, accepted for publication in ApJS

A study of the stellar photosphere-hydrogen ionisation front interaction in pulsating variables using period-color relations

Authors: Susmita Das, Shashi M. Kanbur, Earl P. Bellinger, Anupam Bhardwaj, Harinder P. Singh

Abstract: Period-color (PC) relations may be used to study the interaction of the stellar photosphere and the hydrogen ionization front (HIF). RR Lyraes (RRLs) and long period classical Cepheids (P > 10d) have been found to exhibit different PC behavior at minimum and maximum light which can be explained by the HIF-photosphere interaction based on their location on the HR diagram. In this work, we extend th… ▽ More Period-color (PC) relations may be used to study the interaction of the stellar photosphere and the hydrogen ionization front (HIF). RR Lyraes (RRLs) and long period classical Cepheids (P > 10d) have been found to exhibit different PC behavior at minimum and maximum light which can be explained by the HIF-photosphere interaction based on their location on the HR diagram. In this work, we extend the study to include type II Cepheids (T2Cs) with an aim to test the HIF-photosphere interaction theory across a broad spectrum of variable star types. We find W Vir stars and BL Her stars to have similar PC relations as those from long period and short period classical Cepheids, respectively. We also use MESA to compute RRL, BL Her and classical Cepheid models to study the theoretical HIF-photosphere distance and find the results to be fairly consistent with the HIF-photosphere interaction theory. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: 5 pages, 3 figures. Accepted for publication in ASP Conference Series, Vol. 529 as a conference proceeding for the RR Lyrae/Cepheid 2019 Conference "Frontiers of Classical Pulsators: Theory and Observations" held in Cloudcroft, New Mexico, USA, 13-18 October 2019

Asteroseismic Inference of Subgiant Evolutionary Parameters with Deep Learning

Authors: Marc Hon, Earl P. Bellinger, Saskia Hekker, Dennis Stello, James S. Kuszlewicz

Abstract: With the observations of an unprecedented number of oscillating subgiant stars expected from NASA's TESS mission, the asteroseismic characterization of subgiant stars will be a vital task for stellar population studies and for testing our theories of stellar evolution. To determine the fundamental properties of a large sample of subgiant stars efficiently, we developed a deep learning method that… ▽ More With the observations of an unprecedented number of oscillating subgiant stars expected from NASA's TESS mission, the asteroseismic characterization of subgiant stars will be a vital task for stellar population studies and for testing our theories of stellar evolution. To determine the fundamental properties of a large sample of subgiant stars efficiently, we developed a deep learning method that estimates distributions of fundamental parameters like age and mass over a wide range of input physics by learning from a grid of stellar models varied in eight physical parameters. We applied our method to four Kepler subgiant stars and compare our results with previously determined estimates. Our results show good agreement with previous estimates for three of them (KIC 11026764, KIC 10920273, KIC 11395018). With the ability to explore a vast range of stellar parameters, we determine that the remaining star, KIC 10005473, is likely to have an age 1 Gyr younger than its previously determined estimate. Our method also estimates the efficiency of overshooting, undershooting, and microscopic diffusion processes, from which we determined that the parameters governing such processes are generally poorly-constrained in subgiant models. We further demonstrate our method's utility for ensemble asteroseismology by characterizing a sample of 30 Kepler subgiant stars, where we find a majority of our age, mass, and radius estimates agree within uncertainties from more computationally expensive grid-based modelling techniques. △ Less

Submitted 1 October, 2020; v1 submitted 15 September, 2020; originally announced September 2020.

Comments: 18 pages, 15 figures. Published in MNRAS. Updated citations and amendments to Figure 2 reference in V2

Tess asteroseismology of the known planet host star $λ^2$ Fornacis

Authors: M. B. Nielsen, W. H. Ball, M. R. Standing, A. H. M. J. Triaud, D. Buzasi, L. Carboneau, K. G. Stassun, S. R. Kane, W. J. Chaplin, E. P. Bellinger, B. Mosser, I. W. Roxburgh, Z. Çelik Orhan, M. Yıldız, S. Örtel, M. Vrard, A. Mazumdar, P. Ranadive, M. Deal, G. R. Davies, T. L. Campante, R. A. García, S. Mathur, L. González-Cuesta, A. Serenelli

Abstract: The Transiting Exoplanet Survey Satellite (TESS) is observing bright known planet-host stars across almost the entire sky. These stars have been subject to extensive ground-based observations, providing a large number of radial velocity (RV) measurements. In this work we use the new TESS photometric observations to characterize the star $λ^2$ Fornacis, and following this to update the parameters o… ▽ More The Transiting Exoplanet Survey Satellite (TESS) is observing bright known planet-host stars across almost the entire sky. These stars have been subject to extensive ground-based observations, providing a large number of radial velocity (RV) measurements. In this work we use the new TESS photometric observations to characterize the star $λ^2$ Fornacis, and following this to update the parameters of the orbiting planet $λ^2$ For b. We measure the p-mode oscillation frequencies in $λ^2$ For, and in combination with non-seismic parameters estimate the stellar fundamental properties using stellar models. Using the revised stellar properties and a time series of archival RV data from the UCLES, HIRES and HARPS instruments spanning almost 20 years, we refit the orbit of $λ^2$ For b and search the RV residuals for remaining variability. We find that $λ^2$ For has a mass of $1.16\pm0.03$M$_\odot$ and a radius of $1.63\pm0.04$R$_\odot$, with an age of $6.3\pm0.9$Gyr. This and the updated RV measurements suggest a mass of $λ^2$ For b of $16.8^{+1.2}_{-1.3}$M$_\oplus$, which is $\sim5$M$_\oplus$ less than literature estimates. We also detect a periodicity at 33 days in the RV measurements, which is likely due to the rotation of the host star. While previous literature estimates of the properties of $λ^2$ are ambiguous, the asteroseismic measurements place the star firmly at the early stage of its subgiant evolutionary phase. Typically only short time series of photometric data are available from TESS, but by using asteroseismology it is still possible to provide tight constraints on the properties of bright stars that until now have only been observed from the ground. This prompts a reexamination of archival RV data from the past few decades to update the characteristics of the planet hosting systems observed by TESS for which asteroseismology is possible. △ Less

Submitted 1 July, 2020; originally announced July 2020.

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

Journal ref: A&A 641, A25 (2020)

Asteroseismic sensitivity to internal rotation along the red-giant branch

Authors: F. Ahlborn, E. P. Bellinger, S. Hekker, S. Basu, G. C. Angelou

Abstract: Transport of angular momentum in stellar interiors is currently not well understood. Asteroseismology can provide us with estimates of internal rotation of stars and thereby advances our understanding of angular momentum transport. We can measure core-rotation rates in red-giant stars and we can place upper bounds on surface-rotation rates using measurements of dipole ($l=1$) modes. Here, we aim t… ▽ More Transport of angular momentum in stellar interiors is currently not well understood. Asteroseismology can provide us with estimates of internal rotation of stars and thereby advances our understanding of angular momentum transport. We can measure core-rotation rates in red-giant stars and we can place upper bounds on surface-rotation rates using measurements of dipole ($l=1$) modes. Here, we aim to determine the theoretical sensitivity of modes of different spherical degree towards the surface rotation. Additionally, we aim to identify modes that can potentially add sensitivity at intermediate radii. We used asteroseismic rotational inversions to probe the internal stellar rotation profiles in red-giant models from the base of the red-giant branch up to the luminosity bump. We used the inversion method of multiplicative optimally localised averages (MOLA) to assess how well internal and surface rotation rates can be recovered from different mode sets and different synthetic rotation profiles. We confirm that dipole mixed modes are sufficient to set constraints on the average core-rotation rates in red giants. However, surface-rotation rates estimated with only dipole mixed modes are contaminated by the core rotation. We show that the sensitivity to the surface rotation decreases from the base of the red-giant branch until it reaches a minimum at 0.6-0.8$L_\text{bump}$ due to a glitch in the buoyancy frequency. Thereafter a narrow range of increased surface sensitivity just below the bump luminosity exists. Quadrupole and octopole modes have more sensitivity in the outer parts of the star. If observed, quadrupole and octopole modes enable us to distinguish between differential and solid body rotation in the convection zone. To obtain accurate estimates of rotation rates at intermediate radii, acoustic oscillation modes with a spherical degree of $l\approx10$ are needed. △ Less

Submitted 16 June, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: accepted for publication in Astronomy and Astrophysics, revised manuscript after language editing

Journal ref: A&A 639, A98 (2020)

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

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

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

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

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

Convective boundary mixing in low- and intermediate-mass stars I. Core properties from pressure-mode asteroseismology

Authors: George C Angelou, Earl P. Bellinger, Saskia Hekker, Yvonne Elsworth, Sarbani Basu, Achim Weiss

Abstract: Convective boundary mixing (CBM) is ubiquitous in stellar evolution. It is a necessary ingredient in the models in order to match observational constraints from clusters, binaries and single stars alike. We compute `effective overshoot' measures that reflect the extent of mixing and which can differ significantly from the input overshoot values set in the stellar evolution codes. We use constraint… ▽ More Convective boundary mixing (CBM) is ubiquitous in stellar evolution. It is a necessary ingredient in the models in order to match observational constraints from clusters, binaries and single stars alike. We compute `effective overshoot' measures that reflect the extent of mixing and which can differ significantly from the input overshoot values set in the stellar evolution codes. We use constraints from pressure modes to infer the CBM properties of Kepler and CoRoT main-sequence and subgiant oscillators, as well as in two radial velocity targets (Procyon A and $α$ Cen A). Collectively these targets allow us to identify how measurement precision, stellar spectral type, and overshoot implementation impact the asteroseismic solution. With these new measures we find that the `effective overshoot' for most stars is in line with physical expectations and calibrations from binaries and clusters. However, two F-stars in the CoRoT field (HD 49933 and HD 181906) still necessitate high overshoot in the models. Due to short mode lifetimes, mode identification can be difficult in these stars. We demonstrate that an incongruence between the radial and non-radial modes drives the asteroseismic solution to extreme structures with highly-efficient CBM as an inevitable outcome. Understanding the cause of seemingly anomalous physics for such stars is vital for inferring accurate stellar parameters from TESS data with comparable timeseries length. △ Less

Submitted 16 February, 2020; v1 submitted 6 February, 2020; originally announced February 2020.

Comments: 20 Pages, 9 Figures, 4 Tables and supplementary material. Accepted for publication in MNRAS

The stellar photosphere-hydrogen ionization front interaction in Classical Pulsators: a theoretical explanation for observed period-colour relations

Authors: Susmita Das, Shashi M. Kanbur, Earl P. Bellinger, Anupam Bhardwaj, Harinder P. Singh, Brett Meerdink, Nicholas Proietti, Anthony Chalmers, Ryan Jordan

Abstract: Period-colour and amplitude-colour (PCAC) relations can be used to probe both the hydrodynamics of outer envelope structure and evolutionary status of Cepheids and RR Lyraes. In this work, we incorporate the PCAC relations for RR Lyraes, BL Her, W Vir and classical Cepheids in a single unifying theory that involves the interaction of the hydrogen ionization front (HIF) and stellar photosphere and… ▽ More Period-colour and amplitude-colour (PCAC) relations can be used to probe both the hydrodynamics of outer envelope structure and evolutionary status of Cepheids and RR Lyraes. In this work, we incorporate the PCAC relations for RR Lyraes, BL Her, W Vir and classical Cepheids in a single unifying theory that involves the interaction of the hydrogen ionization front (HIF) and stellar photosphere and the theory of stellar evolution. PC relations for RR Lyraes and classical Cepheids using OGLE-IV data are found to be consistent with this theory: RR Lyraes have shallow/sloped relations at minimum/maximum light whilst long-period ($P>10$ days) Cepheids exhibit sloped/flat PC relations at minimum/maximum light. The differences in the PC relations for Cepheids and RR Lyraes can be explained based on the relative location of the HIF and stellar photosphere which changes depending on their position on the HR diagram. We also extend our analysis of PCAC relations for type II Cepheids in the Galactic bulge, LMC and SMC using OGLE-IV data. We find that BL Her stars have sloped PC relations at maximum and minimum light similar to short-period ($P<10$ days) classical Cepheids. W Vir stars exhibit sloped/flat PC relation at minimum/maximum light similar to long-period classical Cepheids. We also compute state-of-the-art 1D radiation hydrodynamic models of RR Lyraes, BL Her and classical Cepheids using the radial stellar pulsation code in MESA to further test these ideas theoretically and find that the models are generally consistent with this picture. We are thus able to explain PC relations at maximum and minimum light across a broad spectrum of variable star types. △ Less

Submitted 19 January, 2020; originally announced January 2020.

Comments: 20 pages, 14 figures. Accepted for publication in MNRAS

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

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

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

Submitted 10 December, 2019; originally announced December 2019.

Comments: accepted for publication in Astronomy and Astrophysics

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

A seismic scaling relation for stellar age II: The red giant branch

Authors: Earl Patrick Bellinger

Abstract: Owing to their simplicity and ease of application, seismic scaling relations are widely used to determine the properties of stars exhibiting solar-like oscillations, such as solar twins and red giants. So far, no seismic scaling relations for determining the ages of red giant stars have been developed. Such relations would be desirable for galactic archaeology, which uses stellar ages to map the h… ▽ More Owing to their simplicity and ease of application, seismic scaling relations are widely used to determine the properties of stars exhibiting solar-like oscillations, such as solar twins and red giants. So far, no seismic scaling relations for determining the ages of red giant stars have been developed. Such relations would be desirable for galactic archaeology, which uses stellar ages to map the history of the Milky Way. The ages of red giants must instead be estimated with reference to grids of theoretical stellar models, which can be computationally intensive. Here I present an exhaustive search for scaling age relations involving different combinations of observable quantities. The candidate scaling relations are calibrated and tested using more than 1,000 red giant stars whose ages were obtained via grid-based modeling. I report multiple high-quality scaling relations for red giant branch stars, the best of which are shown to be approximately as accurate as grid-based modeling with typical uncertainties of 15%. Additionally, I present new scaling mass and radius relations for red giants as well. △ Less

Submitted 9 December, 2019; v1 submitted 26 November, 2019; originally announced November 2019.

Comments: 5 pages, 4 figures. Accepted for publication in MNRAS Letters. Source code provided

When a Period Is Not a Full Stop: Light Curve Structure Reveals Fundamental Parameters of Cepheid and RR Lyrae Stars

Authors: Earl P. Bellinger, Shashi M. Kanbur, Anupam Bhardwaj, Marcella Marconi

Abstract: The period of pulsation and the structure of the light curve for Cepheid and RR Lyrae variables depend on the fundamental parameters of the star: mass, radius, luminosity, and effective temperature. Here we train artificial neural networks on theoretical pulsation models to predict the fundamental parameters of these stars based on their period and light curve structure. We find significant improv… ▽ More The period of pulsation and the structure of the light curve for Cepheid and RR Lyrae variables depend on the fundamental parameters of the star: mass, radius, luminosity, and effective temperature. Here we train artificial neural networks on theoretical pulsation models to predict the fundamental parameters of these stars based on their period and light curve structure. We find significant improvements to estimates of these parameters made using light curve structure and period over estimates made using only the period. Given that the models are able to reproduce most observables, we find that the fundamental parameters of these stars can be estimated up to 60% more accurately when light curve structure is taken into consideration. We quantify which aspects of light curve structure are most important in determining fundamental parameters, and find for example that the second Fourier amplitude component of RR Lyrae light curves is even more important than period in determining the effective temperature of the star. We apply this analysis to observations of hundreds Cepheids in the Large Magellanic Cloud and thousands of RR Lyrae in the Magellanic Clouds and Galactic bulge to produce catalogs of estimated masses, radii, luminosities, and other parameters of these stars. As an example application, we estimate Wesenheit indices and use those to derive distance moduli to the Magellanic Clouds of $μ_{\text{LMC},\text{CEP}} = 18.688 \pm 0.093$, $μ_{\text{LMC},\text{RRL}} = 18.52 \pm 0.14$, and $μ_{\text{SMC},\text{RRL}} = 18.88 \pm 0.17$ mag. △ Less

Submitted 25 November, 2019; originally announced November 2019.

Comments: Accepted for publication in MNRAS. Source code available at https://github.com/earlbellinger/Cepheid-neural-network

Testing stellar evolution with asteroseismic inversions of a main sequence star harboring a small convective core

Authors: Earl P. Bellinger, Sarbani Basu, Saskia Hekker, Jørgen Christensen-Dalsgaard

Abstract: The goal of stellar evolution theory is to predict the structure of stars throughout their lifetimes. Usually, these predictions can be assessed only indirectly, for example by comparing predicted and observed effective temperatures and luminosities. Thanks now to asteroseismology, which can reveal the internal structure of stars, it becomes possible to compare the predictions from stellar evoluti… ▽ More The goal of stellar evolution theory is to predict the structure of stars throughout their lifetimes. Usually, these predictions can be assessed only indirectly, for example by comparing predicted and observed effective temperatures and luminosities. Thanks now to asteroseismology, which can reveal the internal structure of stars, it becomes possible to compare the predictions from stellar evolution theory to actual stellar structures. In this work, we present an inverse analysis of the oscillation data from the solar-type star KIC 6225718, which was observed by the \emph{Kepler} space observatory during its nominal mission. As its mass is about 20% greater than solar, this star is predicted to transport energy by convection in its nuclear-burning core. We find significant differences between the predicted and actual structure of the star in the radiative interior near to the convective core. In particular, the predicted sound speed is higher than observed in the deep interior of the star, and too low at a fractional radius of 0.25 and beyond. The cause of these discrepancies is unknown, and is not remedied by known physics in the form of convective overshooting or elemental diffusion. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: 9 pages, 6 figures, accepted for publication in The Astrophysical Journal

Asteroseismic constraints on the cosmic-time variation of the gravitational constant from an ancient main-sequence star

Authors: Earl Patrick Bellinger, Jørgen Christensen-Dalsgaard

Abstract: We investigate the variation of the gravitational constant $G$ over the history of the Universe by modeling the effects on the evolution and asteroseismology of the low-mass star KIC 7970740, which is one of the oldest (~11 Gyr) and best-observed solar-like oscillators in the Galaxy. From these data we find $\dot{G}/G = (1.2 \pm 2.6) \times 10^{-12}~\text{yr}^{-1}$, that is, no evidence for any va… ▽ More We investigate the variation of the gravitational constant $G$ over the history of the Universe by modeling the effects on the evolution and asteroseismology of the low-mass star KIC 7970740, which is one of the oldest (~11 Gyr) and best-observed solar-like oscillators in the Galaxy. From these data we find $\dot{G}/G = (1.2 \pm 2.6) \times 10^{-12}~\text{yr}^{-1}$, that is, no evidence for any variation in $G$. We also find a Bayesian asteroseismic estimate of the age of the Universe as well as astrophysical S-factors for five nuclear reactions obtained through a 12-dimensional stellar evolution Markov chain Monte Carlo simulation. △ Less

Submitted 11 November, 2019; v1 submitted 13 September, 2019; originally announced September 2019.

Comments: 5 pages, 3 figures, accepted for publication in The Astrophysical Journal Letters

The Detailed Science Case for the Maunakea Spectroscopic Explorer, 2019 edition

Authors: The MSE Science Team, Carine Babusiaux, Maria Bergemann, Adam Burgasser, Sara Ellison, Daryl Haggard, Daniel Huber, Manoj Kaplinghat, Ting Li, Jennifer Marshall, Sarah Martell, Alan McConnachie, Will Percival, Aaron Robotham, Yue Shen, Sivarani Thirupathi, Kim-Vy Tran, Christophe Yeche, David Yong, Vardan Adibekyan, Victor Silva Aguirre, George Angelou, Martin Asplund, Michael Balogh, Projjwal Banerjee , et al. (239 additional authors not shown)

Abstract: (Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the sc… ▽ More (Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more. △ Less

Submitted 9 April, 2019; originally announced April 2019.

Comments: 9 chapters, 301 pages, 100 figures. This version of the DSC is a comprehensive update of the original version, released in 2016, which can be downloaded at arXiv:1606.00043. A detailed summary of the design of MSE is available in the MSE Book 2018, available at arXiv:1810.08695

Stellar Astrophysics and Exoplanet Science with the Maunakea Spectroscopic Explorer (MSE)

Authors: Maria Bergemann, Daniel Huber, Vardan Adibekyan, George Angelou, Daniela Barría, Timothy C. Beers, Paul G. Beck, Earl P. Bellinger, Joachim M. Bestenlehner, Bertram Bitsch, Adam Burgasser, Derek Buzasi, Santi Cassisi, Márcio Catelan, Ana Escorza, Scott W. Fleming, Boris T. Gänsicke, Davide Gandolfi, Rafael A. García, Mark Gieles, Amanda Karakas, Yveline Lebreton, Nicolas Lodieu, Carl Melis, Thibault Merle , et al. (48 additional authors not shown)

Abstract: The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fib… ▽ More The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars. △ Less

Submitted 7 March, 2019; originally announced March 2019.

Comments: 31 pages, 11 figures; To appear as a chapter for the Detailed Science Case of the Maunakea Spectroscopic Explorer

A seismic scaling relation for stellar age

Authors: Earl Patrick Bellinger

Abstract: A simple solar scaling relation for estimating the ages of main-sequence stars from asteroseismic and spectroscopic data is developed. New seismic scaling relations for estimating mass and radius are presented as well, including a purely seismic radius scaling relation (i.e., no dependence on temperature). The relations show substantial improvement over the classical scaling relations and perform… ▽ More A simple solar scaling relation for estimating the ages of main-sequence stars from asteroseismic and spectroscopic data is developed. New seismic scaling relations for estimating mass and radius are presented as well, including a purely seismic radius scaling relation (i.e., no dependence on temperature). The relations show substantial improvement over the classical scaling relations and perform similarly well to grid-based modeling. △ Less

Submitted 6 May, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 10 pages, including source code. Accepted for publication in MNRAS

Stellar ages, masses and radii from asteroseismic modeling are robust to systematic errors in spectroscopy

Authors: Earl P. Bellinger, Saskia Hekker, George C. Angelou, Amalie Stokholm, Sarbani Basu

Abstract: The search for twins of the Sun and Earth relies on accurate characterization of stellar and exoplanetary parameters: i.e., ages, masses, and radii. In the modern era of asteroseismology, parameters of solar-like stars are derived by fitting theoretical models to observational data, which include measurements of their oscillation frequencies, metallicity [Fe/H], and effective temperature Teff. Com… ▽ More The search for twins of the Sun and Earth relies on accurate characterization of stellar and exoplanetary parameters: i.e., ages, masses, and radii. In the modern era of asteroseismology, parameters of solar-like stars are derived by fitting theoretical models to observational data, which include measurements of their oscillation frequencies, metallicity [Fe/H], and effective temperature Teff. Combining this information with transit data furthermore yields the corresponding parameters for their exoplanets. While [Fe/H] and Teff are commonly stated to a precision of ~0.1 dex and ~100 K, the impact of errors in their measurement has not been studied in practice within the context of the parameters derived from them. Here we use the Stellar Parameters in an Instant (SPI) pipeline to estimate the parameters of nearly 100 stars observed by Kepler and Gaia, many of which are confirmed planet hosts. We adjust the reported spectroscopic measurements of these stars by introducing faux systematic errors and artificially increasing the reported uncertainties, and quantify the differences in the resulting parameters. We find that a systematic error of 0.1 dex in [Fe/H] translates to differences of only 4%, 2%, and 1% on average in the resulting stellar ages, masses, and radii, which are well within their uncertainties (~11%, 3.5%, 1.4%) as derived by SPI. We also find that increasing the uncertainty of [Fe/H] measurements by 0.1 dex increases the uncertainties by only 0.01 Gyr, 0.02 M_sun, and 0.01 R_sun, which are again well below their reported uncertainties (0.5 Gyr, 0.04 M_sun, 0.02 R_sun). The results for Teff at 100 K are similar. Stellar parameters from SPI are unchanged within uncertainties by errors of up to 0.14 dex or 175 K, and are even more robust to errors in Teff than the seismic scaling relations. Consequently, the parameters for their exoplanets are robust as well. △ Less

Submitted 5 February, 2020; v1 submitted 17 December, 2018; originally announced December 2018.

Comments: Accepted for publication in A&A

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

Time-series analyses of Cepheid and RR Lyrae variables in the wide-field variability surveys

Authors: A. Bhardwaj, S. M. Kanbur, M. Marconi, S. Das, E. P. Bellinger, H. P. Singh, M. Rejkuba, C. -C. Ngeow

Abstract: We discuss time-series analyses of classical Cepheid and RR Lyrae variables in the Galaxy and the Magellanic Clouds at multiple wavelengths. We adopt the Fourier decomposition method to quantify the structural changes in the light curves of Cepheid and RR Lyrae variables. A quantitative comparison of Cepheid Fourier parameters suggests that the canonical mass-luminosity models that lie towards the… ▽ More We discuss time-series analyses of classical Cepheid and RR Lyrae variables in the Galaxy and the Magellanic Clouds at multiple wavelengths. We adopt the Fourier decomposition method to quantify the structural changes in the light curves of Cepheid and RR Lyrae variables. A quantitative comparison of Cepheid Fourier parameters suggests that the canonical mass-luminosity models that lie towards the red-edge of the instability strip show a greater offset with respect to observations for short-period Cepheids. RR Lyrae models are consistent with observations in most period bins. We use ensemble light curve analysis to predict the physical parameters of observed Cepheid and RR Lyrae variables using machine learning methods. Our preliminary results suggest that the posterior distributions of mass, luminosity, temperature and radius for Cepheids and RR Lyraes can be well-constrained for a given metal abundance, provided a smoother grid of models is adopted in various input physical parameters. △ Less

Submitted 16 October, 2018; originally announced October 2018.

Comments: 4 pages, 2 figures, Proceedings for IAUS347, 'Early Science with ELTs', IAU-GA, Vienna, August 2018

Asteroseismology of KIC 8263801:Is it a member of NGC 6866 and a red clump star?

Authors: Yanke Tang, Sarbani Basu, Guy R. Davies, Earl P. Bellinger, Ning Gai

Abstract: We present an asteroseismic analysis of the Kepler light curve of KIC 8263801, a red-giant star in the open cluster NGC 6866 that has previously been reported to be a helium-burning red-clump star. We extracted the frequencies of the radial and quadrupole modes from its frequency power spectrum and determined its properties using a grid of evolutionary models constructed with MESA. The oscillation… ▽ More We present an asteroseismic analysis of the Kepler light curve of KIC 8263801, a red-giant star in the open cluster NGC 6866 that has previously been reported to be a helium-burning red-clump star. We extracted the frequencies of the radial and quadrupole modes from its frequency power spectrum and determined its properties using a grid of evolutionary models constructed with MESA. The oscillation frequencies were calculated using the GYRE code and the surface term was corrected using the Ball & Gizon(2014) prescription. We find that the star has a mass of $M/M_{\odot} = 1.793\pm 0.072$, age $t=1.48\pm 0.21$ Gyr and radius $R/R_{\odot} = 10.53\pm 0.28$. By analyzing the internal structure of the best-fitting model, we infer the evolutionary status of the star KIC 8263801 as being on the ascending part of the red giant branch, and not on the red clump. This result is verified using a purely asteroseismic diagnostic, the $ε_{c}-Δν_{c}$ diagram which can distinguish red giant branch stars from red clump stars. Finally, by comparing its age with NGC 6866 ($t=0.65 \pm 0.1$ Gyr) we conclude that KIC 8263801 is not a member of this open cluster. △ Less

Submitted 1 September, 2018; originally announced September 2018.

Inverse Problems in Asteroseismology

Authors: Earl Patrick Bellinger

Abstract: Asteroseismology allows us to probe the internal structure of stars through their global modes of oscillation. Thanks to missions such as the NASA Kepler space observatory, we now have high-quality asteroseismic data for nearly 100 solar-type stars. In this thesis, new techniques to measure the ages, masses, and radii of stars are presented, as well as a way to infer their internal structure. Asteroseismology allows us to probe the internal structure of stars through their global modes of oscillation. Thanks to missions such as the NASA Kepler space observatory, we now have high-quality asteroseismic data for nearly 100 solar-type stars. In this thesis, new techniques to measure the ages, masses, and radii of stars are presented, as well as a way to infer their internal structure. △ Less

Submitted 20 August, 2018; originally announced August 2018.

Comments: 256 pages, 68 figures, Ph.D. dissertation in Computer Science at the Georg-August-University School of Science (GAUSS) Göttingen undertaken in the context of the International Max Planck Research School for Solar System Science under the supervision of Dr. ir. Saskia Hekker (Group Leader, SAGE Group) and Prof. Dr. Sarbani Basu (Chair, Department of Astronomy, Yale University)

Model-independent measurement of internal stellar structure in 16 Cygni A and B

Authors: Earl P. Bellinger, Sarbani Basu, Saskia Hekker, Warrick H. Ball

Abstract: We present a method for measuring internal stellar structure based on asteroseismology that we call "inversions for agreement." The method accounts for imprecise estimates of stellar mass and radius as well as the relatively limited oscillation mode sets that are available for distant stars. By construction, the results of the method are independent of stellar models. We apply this method to measu… ▽ More We present a method for measuring internal stellar structure based on asteroseismology that we call "inversions for agreement." The method accounts for imprecise estimates of stellar mass and radius as well as the relatively limited oscillation mode sets that are available for distant stars. By construction, the results of the method are independent of stellar models. We apply this method to measure the isothermal sound speeds in the cores of the solar-type stars 16 Cyg A and B using asteroseismic data obtained from Kepler observations. We compare the asteroseismic structure that we deduce against best-fitting evolutionary models and find that the sound speeds in the cores of these stars exceed those of the models. △ Less

Submitted 7 December, 2017; v1 submitted 31 October, 2017; originally announced October 2017.

Comments: Accepted for publication in The Astrophysical Journal

Mitigating the mass dependence in the $Δν$ scaling relation of red-giant stars

Authors: Elisabeth Guggenberger, Saskia Hekker, George C. Angelou, Sarbani Basu, Earl P. Bellinger

Abstract: The masses and radii of solar-like oscillators can be estimated through the asteroseismic scaling relations. These relations provide a direct link between observables, i.e. effective temperature and characteristics of the oscillation spectra, and stellar properties, i.e. mean density and surface gravity (thus mass and radius). These scaling relations are commonly used to characterize large samples… ▽ More The masses and radii of solar-like oscillators can be estimated through the asteroseismic scaling relations. These relations provide a direct link between observables, i.e. effective temperature and characteristics of the oscillation spectra, and stellar properties, i.e. mean density and surface gravity (thus mass and radius). These scaling relations are commonly used to characterize large samples of stars. Usually, the Sun is used as a reference from which the structure is scaled. However, for stars that do not have a similar structure as the Sun, using the Sun as a reference introduces systematic errors as large as 10\% in mass and 5\% in radius. Several alternatives for the reference of the scaling relation involving the large frequency separation (typical frequency difference between modes of the same degree and consecutive radial order) have been suggested in the literature. In a previous paper, we presented a reference function with a dependence on both effective temperature and metallicity. The accuracy of predicted masses and radii improved considerably when using reference values calculated from our reference function. However, the residuals indicated that stars on the red-giant branch possess a mass dependence that was not accounted for. Here, we present a reference function for the scaling relation involving the large frequency separation that includes the mass dependence. This new reference function improves the derived masses and radii significantly by removing the systematic differences and mitigates the trend with $ν_{\rm max}$ (frequency of maximum oscillation power) that exists when using the solar value as a reference. △ Less

Submitted 19 May, 2017; originally announced May 2017.

Comments: 12 pages, 7 figures, accepted for publication in MNRAS

Stellar Parameters in an Instant with Machine Learning: Application to Kepler LEGACY Targets

Authors: Earl P. Bellinger, George C. Angelou, Saskia Hekker, Sarbani Basu, Warrick H. Ball, Elisabeth Guggenberger

Abstract: With the advent of dedicated photometric space missions, the ability to rapidly process huge catalogues of stars has become paramount. Bellinger and Angelou et al. (2016) recently introduced a new method based on machine learning for inferring the stellar parameters of main-sequence stars exhibiting solar-like oscillations. The method makes precise predictions that are consistent with other method… ▽ More With the advent of dedicated photometric space missions, the ability to rapidly process huge catalogues of stars has become paramount. Bellinger and Angelou et al. (2016) recently introduced a new method based on machine learning for inferring the stellar parameters of main-sequence stars exhibiting solar-like oscillations. The method makes precise predictions that are consistent with other methods, but with the advantages of being able to explore many more parameters while costing practically no time. Here we apply the method to 52 so-called "LEGACY" main-sequence stars observed by the Kepler space mission. For each star, we present estimates and uncertainties of mass, age, radius, luminosity, core hydrogen abundance, surface helium abundance, surface gravity, initial helium abundance, and initial metallicity as well as estimates of their evolutionary model parameters of mixing length, overshooting coefficient, and diffusion multiplication factor. We obtain median uncertainties in stellar age, mass, and radius of 14.8%, 3.6%, and 1.7%, respectively. The source code for all analyses and for all figures appearing in this manuscript can be found electronically at: https://github.com/earlbellinger/asteroseismology △ Less

Submitted 18 May, 2017; originally announced May 2017.

Comments: 4 pages, 3 figures, 2 tables, to appear in the proceedings of the joint TASC2/KASC9/SPACEINN/HELAS8 conference "Seismology of the Sun and the Distant Stars 2016"

On The Statistical Properties of the Lower Main Sequence

Authors: George C. Angelou, Earl P. Bellinger, Saskia Hekker, Sarbani Basu

Abstract: Astronomy is in an era where all-sky surveys are mapping the Galaxy. The plethora of photometric, spectroscopic, asteroseismic and astrometric data allows us to characterise the comprising stars in detail. Here we quantify to what extent precise stellar observations reveal information about the properties of a star, including properties that are unobserved, or even unobservable. We analyse the dia… ▽ More Astronomy is in an era where all-sky surveys are mapping the Galaxy. The plethora of photometric, spectroscopic, asteroseismic and astrometric data allows us to characterise the comprising stars in detail. Here we quantify to what extent precise stellar observations reveal information about the properties of a star, including properties that are unobserved, or even unobservable. We analyse the diagnostic potential of classical and asteroseismic observations for inferring stellar parameters such as age, mass and radius from evolutionary tracks of solar-like oscillators on the lower main sequence. We perform rank correlation tests in order to determine the capacity of each observable quantity to probe structural components of stars and infer their evolutionary histories. We also analyse the principal components of classic and asteroseismic observables to highlight the degree of redundancy present in the measured quantities and demonstrate the extent to which information of the model parameters can be extracted. We perform multiple regression using combinations of observable quantities in a grid of evolutionary simulations and appraise the predictive utility of each combination in determining the properties of stars. We identify the combinations that are useful and provide limits to where each type of observable quantity can reveal information about a star. We investigate the accuracy with which targets in the upcoming TESS and PLATO missions can be characterized. We demonstrate that the combination of observations from GAIA and PLATO will allow us to tightly constrain stellar masses, ages and radii with machine learning for the purposes of Galactic and planetary studies. △ Less

Submitted 29 March, 2017; originally announced March 2017.

Comments: 37 Pages, 13 Figures, 18 Tables, accepted for publication in ApJ

Fundamental Parameters of Main-Sequence Stars in an Instant with Machine Learning

Authors: Earl P. Bellinger, George C. Angelou, Saskia Hekker, Sarbani Basu, Warrick Ball, Elisabeth Guggenberger

Abstract: Owing to the remarkable photometric precision of space observatories like Kepler, stellar and planetary systems beyond our own are now being characterized en masse for the first time. These characterizations are pivotal for endeavors such as searching for Earth-like planets and solar twins, understanding the mechanisms that govern stellar evolution, and tracing the dynamics of our Galaxy. The volu… ▽ More Owing to the remarkable photometric precision of space observatories like Kepler, stellar and planetary systems beyond our own are now being characterized en masse for the first time. These characterizations are pivotal for endeavors such as searching for Earth-like planets and solar twins, understanding the mechanisms that govern stellar evolution, and tracing the dynamics of our Galaxy. The volume of data that is becoming available, however, brings with it the need to process this information accurately and rapidly. While existing methods can constrain fundamental stellar parameters such as ages, masses, and radii from these observations, they require substantial computational efforts to do so. We develop a method based on machine learning for rapidly estimating fundamental parameters of main-sequence solar-like stars from classical and asteroseismic observations. We first demonstrate this method on a hare-and-hound exercise and then apply it to the Sun, 16 Cyg A & B, and 34 planet-hosting candidates that have been observed by the Kepler spacecraft. We find that our estimates and their associated uncertainties are comparable to the results of other methods, but with the additional benefit of being able to explore many more stellar parameters while using much less computation time. We furthermore use this method to present evidence for an empirical diffusion-mass relation. Our method is open source and freely available for the community to use. The source code for all analyses and for all figures appearing in this manuscript can be found electronically at https://github.com/earlbellinger/asteroseismology △ Less

Submitted 6 July, 2016; originally announced July 2016.

Comments: 26 pages, 18 figures, accepted for publication in ApJ

Measuring amplitudes of harmonics and combination frequencies in variable stars

Authors: Earl P. Bellinger, Daniel Wysocki, Shashi M. Kanbur

Abstract: Discoveries of RR Lyrae and Cepheid variable stars with multiple modes of pulsation have increased tremendously in recent years. The Fourier spectra of these stars can be quite complicated due to the large number of combination frequencies that can exist between their modes. As a result, light-curve fits to these stars often suffer from undesirable ringing effects that arise from noisy observation… ▽ More Discoveries of RR Lyrae and Cepheid variable stars with multiple modes of pulsation have increased tremendously in recent years. The Fourier spectra of these stars can be quite complicated due to the large number of combination frequencies that can exist between their modes. As a result, light-curve fits to these stars often suffer from undesirable ringing effects that arise from noisy observations and poor phase coverage. These non-physical overfitting artifacts also occur when fitting the harmonics of single-mode stars as well. Here we present a new method for fitting light curves that is much more robust against these effects. We prove that the amplitude measurement problem is very difficult (NP-hard) and provide a heuristic algorithm for solving it quickly and accurately. △ Less

Submitted 30 November, 2015; originally announced December 2015.

Comments: 4 pages, 3 figures, RR Lyrae 2015, to appear in the No. 105 issue of the Communications from the Konkoly Observatory of the Hungarian Academy of Sciences

Journal ref: Communications from the Konkoly Observatory, 2016, Vol. 105, p. 101-104

Period-Luminosity Relations for Cepheid Variables: From Mid-Infrared to Multi-Phase

Authors: Chow-Choong Ngeow, Shashi M. Kanbur, Earl P. Bellinger, Marcella Marconi, Ilaria Musella, Michele Cignoni, Ya-Hong Lin

Abstract: This paper discusses two aspects of current research on the Cepheid period-luminosity (P-L) relation: the derivation of mid-infrared (MIR) P-L relations and the investigation of multi-phase P-L relations. The MIR P-L relations for Cepheids are important in the James Webb Space Telescope era for the distance scale issue, as the relations have potential to derive the Hubble constant within ~ 2% accu… ▽ More This paper discusses two aspects of current research on the Cepheid period-luminosity (P-L) relation: the derivation of mid-infrared (MIR) P-L relations and the investigation of multi-phase P-L relations. The MIR P-L relations for Cepheids are important in the James Webb Space Telescope era for the distance scale issue, as the relations have potential to derive the Hubble constant within ~ 2% accuracy - a critical constraint in precision cosmology. Consequently, we have derived the MIR P-L relations for Cepheids in the Large and Small Magellanic Clouds, using archival data from Spitzer Space Telescope. We also compared currently empirical P-L relations for Cepheids in the Magellanic Clouds to the synthetic MIR P-L relations derived from pulsational models. For the study of multi-phase P-L relations, we present convincing evidence that the Cepheid P-L relations in the Magellanic Clouds are highly dynamic quantities that vary significantly when considered as a function of pulsational phase. We found that there is a difference in P-L relations as a function of phase between the Cepheids in each of the Clouds; the most likely cause for this is the metallicity difference between the two galaxies. We also investigated the dispersion of the multi-phase P-L relations, and found that the minimum dispersion do not differ significantly from the mean light P-L dispersion. △ Less

Submitted 2 February, 2012; v1 submitted 1 February, 2012; originally announced February 2012.

Comments: 10 pages, 7 figures and 1 table; Proceeding for presentation at "The Fundamental Cosmic Distance Scale: State of the Art and the Gaia Perspective", Naples, May 2011. Ap&SS accepted