Showing 1–50 of 155 results for author: Joyce, M

Barium stars as tracers of s-process nucleosynthesis in AGB stars III. Systematic deviations from the AGB models

Authors: B. Világos, B. Cseh, A. Yagüe López, M. Joyce, A. Karakas, G. Tagliente, M. Lugaro

Abstract: Barium (Ba) stars help to verify asymptotic giant branch (AGB) star nucleosynthesis models since they experienced pollution from an AGB binary companion and thus their spectra carry the signatures of the slow neutron capture process (s process). For 180 Ba stars, we searched for AGB stellar models that match the observed abundance patterns. We employed three machine learning algorithms as classifi… ▽ More Barium (Ba) stars help to verify asymptotic giant branch (AGB) star nucleosynthesis models since they experienced pollution from an AGB binary companion and thus their spectra carry the signatures of the slow neutron capture process (s process). For 180 Ba stars, we searched for AGB stellar models that match the observed abundance patterns. We employed three machine learning algorithms as classifiers: a Random Forest method, developed for this work, and the two classifiers used in our previous study. We studied the statistical behaviour of the s-process elements in the observational sample to investigate if the AGB models systematically under- or overpredict the abundances observed in the Ba stars and show the results in the form of violin plots of the residuals between spectroscopic abundances and model predictions. We find a significant trend in the residuals that implies an underproduction of the elements Nb, Mo, and Ru in the models relative to the observations. This may originate from a process (e.g. the intermediate neutron-capture process, i process) at the metallicity of the Ba stars not yet included in the AGB models. Correlations are found between the residuals of these elements, suggesting a common origin for the deviations. In addition, there is a weak metallicity dependence of their residuals. The s-process temperatures derived with the [Zr/Fe] - [Nb/Fe] thermometer have an unrealistic value for the majority of our stars. The most likely explanation is that at least a fraction of these elements are not produced in a steady-state s process, and instead may be due to processes not included in the AGB models. The mass distribution of the identified models confirms that our sample of Ba stars was polluted by low-mass AGB stars. Most of the matching AGB models require low accreted mass, but a few systems with high accreted mass are needed to explain the observations. (abridged) △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 20 pages, 20 figures

At least one in a dozen stars exhibits evidence of planetary ingestion

Authors: Fan Liu, Yuan-Sen Ting, David Yong, Bertram Bitsch, Amanda Karakas, Michael T. Murphy, Meridith Joyce, Aaron Dotter, Fei Dai

Abstract: Stellar chemical compositions can be altered by ingestion of planetary material and/or planet formation which removes refractory material from the proto-stellar disc. These "planet signatures" appear as correlations between elemental abundance differences and the dust condensation temperature. Detecting these planet signatures, however, is challenging due to unknown occurrence rates, small amplitu… ▽ More Stellar chemical compositions can be altered by ingestion of planetary material and/or planet formation which removes refractory material from the proto-stellar disc. These "planet signatures" appear as correlations between elemental abundance differences and the dust condensation temperature. Detecting these planet signatures, however, is challenging due to unknown occurrence rates, small amplitudes, and heterogeneous star samples with large differences in stellar ages, and therefore stars born together (i.e., co-natal) with identical compositions can facilitate such detections. While previous spectroscopic studies were limited to small number of binary stars, the Gaia satellite provides new opportunities for detecting stellar chemical signatures of planets among co-moving pairs of stars confirmed to be co-natal. Here we report high-precision chemical abundances for a homogeneous sample of 91 co-natal pairs of stars with a well-defined selection function and identify at least seven new instances of planetary ingestion, corresponding to an occurrence rate of 8%. An independent Bayesian indicator is deployed, which can effectively disentangle the planet signatures from other factors, such as random abundance variation and atomic diffusion. Our study provides new evidence of planet signatures and facilitates a deeper understanding of the star-planet-chemistry connection by providing new observational constraints on the mechanisms of planet engulfment, formation and evolution. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 29 pages, 11 figures. Author's submitted version before final edits. Published in Nature on March 21, 2024: https://www.nature.com/articles/s41586-024-07091-y

Stellar Evolution in Real Time II: R Hydrae and an Open-Source Grid of >3000 Seismic TP-AGB Models Computed with MESA

Authors: Meridith Joyce, László Molnár, Giulia Cinquegrana, Amanda Karakas, Jamie Tayar, Dóra Tarczay-Nehéz

Abstract: We present a comprehensive characterization of the evolved thermally pulsing asymptotic giant branch (TP-AGB) star R Hydrae, building on the techniques applied in Stellar Evolution in Real Time I (Molnár et al. 2019) to T Ursae Minoris. We compute over 3000 theoretical TP-AGB pulse spectra using MESA and GYRE and combine these with classical observational constraints and nearly 400 years of measur… ▽ More We present a comprehensive characterization of the evolved thermally pulsing asymptotic giant branch (TP-AGB) star R Hydrae, building on the techniques applied in Stellar Evolution in Real Time I (Molnár et al. 2019) to T Ursae Minoris. We compute over 3000 theoretical TP-AGB pulse spectra using MESA and GYRE and combine these with classical observational constraints and nearly 400 years of measurements of R Hya's period evolution to fit R Hya's evolutionary and asteroseismic features. Two hypotheses for the mode driving R Hya's period are considered. Solutions that identify this as the fundamental mode (FM) as well as the first overtone (O1) are consistent with observations. Using a variety of statistical tests, we find that R Hya is most likely driven by the FM and currently occupies the ``power down'' phase of an intermediate pulse (TP ~ 9-16). We predict that its pulsation period will continue to shorten for millennia. Using supplementary calculations from the Monash stellar evolution code, we also find that R Hya is likely to have undergone third dredge-up in its most recent pulse. The MESA+GYRE model grid used in this analysis includes exact solutions to the adiabatic equations of stellar oscillation for the first 10 radial-order pressure modes for every time step in every evolutionary track. The grid is fully open-source and packaged with a data visualization application. This is the first publicly available grid of TP-AGB models with seismology produced with MESA. △ Less

Submitted 19 June, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: accepted to ApJ March 2024. Final revisions complete. Github repository associated to this project: https://github.com/mjoyceGR/AGB_grid_visualizer All grids available on Zenodo: https://zenodo.org/records/11280179 There are FOUR Zenodo DOIs associated to this project, the first of which is above

RR Lyrae Stars Belonging to the Candidate Globular Cluster Patchick 99

Authors: Evan Butler, Andrea Kunder, Zdenek Prudil, Kevin R. Covey, Macy Ball, Carlos Campos, Kaylen Gollnick, Julio Olivares Carvajal, Joanne Hughes, Kathryn Devine, Christian I. Johnson, A. Katherina Vivas, Michael R. Rich, Meridith Joyce, Iulia T. Simon, Tommaso Marchetti, Andreas J. Koch-Hansen, William I. Clarkson, Rebekah Kuss

Abstract: Patchick 99 is a candidate globular cluster located in the direction of the Galactic bulge, with a proper motion almost identical to the field and extreme field star contamination. A recent analysis suggests it is a low-luminosity globular cluster with a population of RR Lyrae stars. We present new spectra of stars in and around Patchick 99, targeting specifically the 3 RR Lyrae stars associated w… ▽ More Patchick 99 is a candidate globular cluster located in the direction of the Galactic bulge, with a proper motion almost identical to the field and extreme field star contamination. A recent analysis suggests it is a low-luminosity globular cluster with a population of RR Lyrae stars. We present new spectra of stars in and around Patchick 99, targeting specifically the 3 RR Lyrae stars associated with the cluster as well as the other RR Lyrae stars in the field. A sample of 53 giant stars selected from proper motions and a position on CMD are also observed. The three RR Lyrae stars associated with the cluster have similar radial velocities and distances, and two of the targeted giants also have radial velocities in this velocity regime and [Fe/H] metallicities that are slightly more metal-poor than the field. Therefore, if Patchick 99 is a bonafide globular cluster, it would have a radial velocity of -92+/-10 km s-1, a distance of 6.7+/-0.4 kpc (as determined from the RR Lyrae stars), and an orbit that confines it to the inner bulge. △ Less

Submitted 25 January, 2024; v1 submitted 24 January, 2024; originally announced January 2024.

Comments: Accepted to The Astrophysical Journal Letters. Replaced due to a typo in the title

A theoretical framework for BL Her stars -- II. New period-luminosity relations in the Gaia passbands

Authors: Susmita Das, László Molnár, Shashi M. Kanbur, Meridith Joyce, Anupam Bhardwaj, Harinder P. Singh, Marcella Marconi, Vincenzo Ripepi, Radoslaw Smolec

Abstract: We present new theoretical period-luminosity (PL) and period-Wesenheit (PW) relations for a fine grid of convective BL Her, the shortest period T2Cs, models computed using MESA-RSP and compare our results with the empirical relations from Gaia DR3. We use the state-of-the-art 1D non-linear radial stellar pulsation tool MESA-RSP to compute models of BL Her stars over a wide range of input parameter… ▽ More We present new theoretical period-luminosity (PL) and period-Wesenheit (PW) relations for a fine grid of convective BL Her, the shortest period T2Cs, models computed using MESA-RSP and compare our results with the empirical relations from Gaia DR3. We use the state-of-the-art 1D non-linear radial stellar pulsation tool MESA-RSP to compute models of BL Her stars over a wide range of input parameters - metallicity (-2.0 dex $\leq$ [Fe/H] $\leq$ 0.0 dex), stellar mass (0.5M$_{\odot}$-0.8M$_{\odot}$), stellar luminosity (50L$_{\odot}$-300L$_{\odot}$) and effective temperature (full extent of the instability strip; in steps of 50K). The BL Her stars in the All Sky region exhibit statistically different PL slopes compared to the theoretical PL slopes computed using the four sets of convection parameters. We find the empirical PL and PW slopes from BL Her stars in the Magellanic Clouds to be statistically consistent with the theoretical relations computed using the different convection parameter sets in the Gaia passbands. There is negligible effect of metallicity on the PL relations in the individual Gaia passbands. However, there exists a small but significant negative coefficient of metallicity in the PWZ relations for the BL Her models using the four sets of convection parameters. This could be attributed to the increased sensitivity of bolometric corrections to metallicities at wavelengths shorter than the V band. Our BL Her models also suggest a dependence of the mass-luminosity relation on metallicity. We found the observed Fourier parameter space to be covered well by our models. Higher mass models (> 0.6M$_{\odot}$) may be needed to reliably model the observed light curves of BL Her stars in the All Sky region. We also found the theoretical light curve structures (especially the Fourier amplitude parameters) to be affected by the choice of convection parameters. △ Less

Submitted 7 February, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: 20 pages, 8 figures, accepted in Astronomy & Astrophysics

The Milky Way Bulge extra-tidal star survey: BH 261 (AL 3)

Authors: Andrea Kunder, Zdenek Prudil, Kevin Covey, Joanne Hughes, Meridith Joyce, Iulia T. Simion, Rebekah Kuss, Carlos Campos, Christian I. Johnson, Catherine A. Pilachowski, Kristen A. Larson, Andreas J. Koch-Hansen, Tommaso Marchetti, Michael R. Rich, Evan Butler, William I. Clarkson, Michael J. Rivet, Kathryn Devine, A. Katherina Vivas, Gabriel I. Perren, Mario Soto, Erika Silva

Abstract: The Milky Way Bulge extra-tidal star survey (MWBest) is a spectroscopic survey with the goal of identifying stripped globular cluster stars from inner Galaxy clusters. In this way, an indication of the fraction of metal-poor bulge stars that originated from globular clusters can be determined. We observed and analyzed stars in and around BH 261, an understudied globular cluster in the bulge. From… ▽ More The Milky Way Bulge extra-tidal star survey (MWBest) is a spectroscopic survey with the goal of identifying stripped globular cluster stars from inner Galaxy clusters. In this way, an indication of the fraction of metal-poor bulge stars that originated from globular clusters can be determined. We observed and analyzed stars in and around BH 261, an understudied globular cluster in the bulge. From seven giants within the tidal radius of the cluster, we measured an average heliocentric radial velocity of <RV> = -61 +- 2.6 km/s with a radial velocity dispersion of σ= 6.1 +- 1.9 km/s. The large velocity dispersion may have arisen from tidal heating in the cluster's orbit about the Galactic center, or because BH 261 has a high dynamical mass as well as a high mass-to-light ratio. From spectra of five giants, we measure an average metallicity of <[Fe/H]> = -1.1 +- 0.2 dex. We also spectroscopically confirm an RR Lyrae star in BH 261, which yields a distance to the cluster of 7.1 +- 0.4~kpc. Stars with 3D velocities and metallicities consistent with BH 261 reaching to ~0.5 degrees from the cluster are identified. A handful of these stars are also consistent with the spatial distribution of that potential debris from models focussing on the most recent disruption of the cluster. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: accepted for publication in The Astronomical Journal

The Blanco DECam Bulge Survey (BDBS) VIII: Chemo-kinematics in the southern Galactic bulge from 2.3 million red clump stars with Gaia DR3 proper motions

Authors: Tommaso Marchetti, Meridith Joyce, Christian Johnson, R. Michael Rich, William Clarkson, Andrea Kunder, Iulia T. Simion, Catherine A. Pilachowski

Abstract: The Blanco DECam Bulge Survey (BDBS) provides near-ultraviolet to near-infrared photometry for ~250 million unique stars. By combining BDBS photometry with the latest Gaia astrometry, we characterize the chemo-dynamics of red clump stars across the BDBS footprint, using an unprecedented sample size and sky coverage. We construct a sample of ~2.3 million red clump giants in the bulge with photometr… ▽ More The Blanco DECam Bulge Survey (BDBS) provides near-ultraviolet to near-infrared photometry for ~250 million unique stars. By combining BDBS photometry with the latest Gaia astrometry, we characterize the chemo-dynamics of red clump stars across the BDBS footprint, using an unprecedented sample size and sky coverage. We construct a sample of ~2.3 million red clump giants in the bulge with photometric metallicities, BDBS photometric distances, and proper motions. We study the kinematics of the red clump stars as a function of sky position and metallicity, by investigating proper motion rotation curves, velocity dispersions, and proper motion correlations across the southern Galactic bulge. We find that metal-poor red clump stars exhibit lower rotation amplitudes, at ~29 km s$^{-1}$ kpc^{-1}. The peak of the angular velocity is ~39 km s^{-1} kpc^{-1} for [Fe/H] ~ -0.2 dex, exhibiting declining rotation at higher [Fe/H]. The velocity dispersion is higher for metal-poor stars, while metal-rich stars show a steeper gradient with Galactic latitude, with a maximum dispersion at low latitudes along the bulge minor axis. Only metal-rich stars ([Fe/H] >~ -0.5 dex) show clear signatures of the bar in their kinematics, while the metal-poor population exhibits isotropic motions with an axisymmetric pattern around Galactic longitude l = 0. This work reports the largest sample of bulge stars with distance, metallicity, and astrometry and shows clear kinematic differences with metallicity. The global kinematics over the bulge agrees with earlier studies. However, we see striking changes with increasing metallicity and for the first time, see kinematic differences for stars with [Fe/H]>-0.5, suggesting that the bar itself may have kinematics that depends on metallicity. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 12 pages, Accepted for publication in A&A

C3PO: Towards a complete census of co-moving pairs of stars. I. High precision stellar parameters for 250 stars

Authors: David Yong, Fan Liu, Yuan-Sen Ting, Meridith Joyce, Bertram Bitsch, Fei Dai, Aaron Dotter, Amanda I. Karakas, Michael T. Murphy

Abstract: We conduct a line-by-line differential analysis of a sample of 125 co-moving pairs of stars (dwarfs and subgiants near solar metallicity). We obtain high precision stellar parameters with average uncertainties in effective temperature, surface gravity and metallicity of 16.5 K, 0.033 dex and 0.014 dex, respectively. We classify the co-moving pairs of stars into two groups, chemically homogeneous (… ▽ More We conduct a line-by-line differential analysis of a sample of 125 co-moving pairs of stars (dwarfs and subgiants near solar metallicity). We obtain high precision stellar parameters with average uncertainties in effective temperature, surface gravity and metallicity of 16.5 K, 0.033 dex and 0.014 dex, respectively. We classify the co-moving pairs of stars into two groups, chemically homogeneous (conatal; |Delta[Fe/H]| $\le$ 0.04 dex) and inhomogeneous (non-conatal), and examine the fraction of chemically homogeneous pairs as a function of separation and effective temperature. The four main conclusions from this study are: (1) A spatial separation of \ds = 10$^6$ AU is an approximate boundary between homogeneous and inhomogeneous pairs of stars, and we restrict our conclusions to only consider the 91 pairs with \ds $\le$ 10$^6$ AU; (2) There is no trend between velocity separation and the fraction of chemically homogeneous pairs in the range \dv $\le$ 4 \kms; (3) We confirm that the fraction of chemically inhomogeneous pairs increases with increasing \teff\ and the trend matches a toy model of that expected from planet ingestion; (4) Atomic diffusion is not the main cause of the chemical inhomogeneity. A major outcome from this study is a sample of 56 bright co-moving pairs of stars with chemical abundance differences $\leq$ 0.02 dex (5\%) which is a level of chemical homogeneity comparable to that of the Hyades open cluster. These important objects can be used, in conjunction with star clusters and the \gaia\ ``benchmark'' stars, to calibrate stellar abundances from large-scale spectroscopic surveys. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: MNRAS in press (see source file for full versions of long tables)

Detailed asteroseismic modeling of RR Lyrae stars with non-radial modes

Authors: H. Netzel, L. Molnar, M. Joyce

Abstract: Photometric observations from the last decade have revealed additional low-amplitude periodicities in many classical pulsators that are likely due to pulsations in non-radial modes. One group of multi-mode RR Lyrae stars, the so-called 0.61 stars, is particularly interesting. In these stars, the radial first overtone is accompanied by additional signals with period ratios around 0.61. The most pro… ▽ More Photometric observations from the last decade have revealed additional low-amplitude periodicities in many classical pulsators that are likely due to pulsations in non-radial modes. One group of multi-mode RR Lyrae stars, the so-called 0.61 stars, is particularly interesting. In these stars, the radial first overtone is accompanied by additional signals with period ratios around 0.61. The most promising explanation for these signals is pulsation in non-radial modes of degrees 8 and 9. If the theory behind the additional signals in the 0.61 stars is substantiated, it would allow us to use non-radial modes to study classical pulsators. We aim to perform asteroseismic modeling of selected 0.61 stars with independently determined physical parameters to test whether this assumption behind the modeling leads to correct results. Namely, we test whether the additional signals are indeed due to non-radial modes of the proposed moderate degrees. We selected a number of and RR Lyrae stars that are also 0.61 stars and have good observational constraints on their other physical parameters. We assume that the nature of those modes is correctly explained with non-radial modes of degrees 8 or 9. Using this assumption and observational constraints on physical parameters, we performed asteroseismic modeling to test whether the observed periods and period ratios can be reproduced. For the majority of selected targets, we obtained a good match between observed and calculated periods and period ratios. For a few targets however, the results obtained are ambiguous and not straightforward to interpret. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 10 pages, 9 figures, 3 tables, accepted for publication in MNRAS

Bridging the Gap between Intermediate and Massive Stars II: $M_\text{mas}$ for the most metal-rich stars and implications for Fe CCSNe rates

Authors: Giulia C. Cinquegrana, Meridith Joyce, Amanda I. Karakas

Abstract: The minimum initial mass required for a star to explode as an Fe core collapse supernova, typically denoted $M_\text{mas}$, is an important quantity in stellar evolution because it defines the border between intermediate mass and massive stellar evolutionary paths. The precise value of $M_\text{mas}$ carries implications for models of galactic chemical evolution and the calculation of star formati… ▽ More The minimum initial mass required for a star to explode as an Fe core collapse supernova, typically denoted $M_\text{mas}$, is an important quantity in stellar evolution because it defines the border between intermediate mass and massive stellar evolutionary paths. The precise value of $M_\text{mas}$ carries implications for models of galactic chemical evolution and the calculation of star formation rates. Despite the fact that stars with super solar metallicities are commonplace within spiral and some giant elliptical galaxies, there are currently no studies of this mass threshold in super metal-rich models with $Z>0.05$. Here, we study the minimum mass necessary for a star to undergo an Fe core collapse supernova when its initial metal content falls in the range $2.5\times 10^{-3} \leq Z \leq 0.10$. Although an increase in initial $Z$ corresponds to an increase in the Fe ignition threshold for $Z \approx 1\times 10^{-3}$ to $Z\approx0.04$, we find that there is a steady reversal in trend that occurs for $Z > 0.05$. Our super metal-rich models thus undergo Fe core collapse at lower initial masses than those required at solar metallicity. Our results indicate that metallicity--dependent curves extending to $Z=0.10$ for the minimum Fe ignition mass should be utilised in galactic chemical evolution simulations to accurately model supernovae rates as a function of metallicity, particularly for simulations of metal-rich spiral and elliptical galaxies. △ Less

Submitted 11 August, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS (August 2023)

To Grow Old and Peculiar: A Survey of Anomalous Variable Stars in M80 and Age Determination using K2 and Gaia

Authors: László Molnár, Emese Plachy, Attila Bódi, András Pál, Meridith Joyce, Csilla Kalup, Christian I. Johnson, Zoltán Dencs, Szabolcs Mészáros, Henryka Netzel, Karen Kinemuchi, Juna A. Kollmeier, Jose Luis Prieto, Aliz Derekas

Abstract: The globular cluster Messier 80 was monitored by the Kepler space telescope for 80 days during the K2 mission. Continuous, high-precision photometry of such an old, compact cluster allows us to study its variable star population in unprecedented detail. We extract light curves for 27 variable stars using differential-image photometry. A search for new variables in the images led to the discovery o… ▽ More The globular cluster Messier 80 was monitored by the Kepler space telescope for 80 days during the K2 mission. Continuous, high-precision photometry of such an old, compact cluster allows us to study its variable star population in unprecedented detail. We extract light curves for 27 variable stars using differential-image photometry. A search for new variables in the images led to the discovery of two new variable stars: an RR Lyrae and a variable red giant star, respectively. Analysis of the RR Lyrae population reveals multiple RRc stars with additional modes and/or peculiar modulation cycles. We newly classify star V28 as a spotted extreme horizontal branch variable. Despite their faintness, we clearly detect the three SX Phe stars but we did not find new pulsation modes beyond the known ones in them. Spectra taken with the VLT and Magellan Clay telescopes, as well as absolute color-magnitude diagrams of the cluster based on Gaia and Pan-STARRS observations confirm the classification of the peculiar modulated variables as bona-fide RRc stars. We propose that they highlight a subgroup of overtone stars that may have been overlooked before. We fit MESA isochrones to the CMDs to estimate the age and metallicity of the cluster. We confirm that M80 is old and metal-poor, but show that isochrone fitting to old populations comes with numerous uncertainties. △ Less

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

Comments: 23 pages, 18 figures, accepted for publication in A&A. v2: minor edits and updates to the text

Journal ref: A&A 678, A104 (2023)

Convergence of halo statistics: code comparison between Rockstar and CompaSO using scale-free simulations

Authors: Sara Maleubre, Daniel J. Eisenstein, Lehman H. Garrison, Michael Joyce

Abstract: In this study, we perform a halo-finder code comparison between Rockstar and CompaSO. Based on our previous analysis aiming at quantifying resolution of $N$-body simulations by exploiting large (up to $N=4096^3$) simulations of scale-free cosmologies run using Abacus, we focus on convergence of the HMF, 2PCF and mean radial pairwise velocities of halo centres selected with the aforementioned two a… ▽ More In this study, we perform a halo-finder code comparison between Rockstar and CompaSO. Based on our previous analysis aiming at quantifying resolution of $N$-body simulations by exploiting large (up to $N=4096^3$) simulations of scale-free cosmologies run using Abacus, we focus on convergence of the HMF, 2PCF and mean radial pairwise velocities of halo centres selected with the aforementioned two algorithms. We establish convergence, for both Rockstar and CompaSO, of mass functions at the $1\%$ precision level and of the mean pairwise velocities (and also 2PCF) at the $2\%$ level. At small scales and small masses, we find that Rockstar exhibits greater self-similarity, and we also highlight the role played by the merger-tree post-processing of CompaSO halos on their convergence. Finally, we give resolution limits expressed as a minimum particle number per halo in a form that can be directly extrapolated to LCDM. △ Less

Submitted 4 December, 2023; v1 submitted 1 August, 2023; originally announced August 2023.

Comments: 13 pages, accepted in MNRAS. arXiv admin note: substantial text overlap with arXiv:2211.07607

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 3, January 2024, Pages 5603 to 5615

A grid of 200,000 models of young $δ$ Scuti stars using MESA and GYRE

Authors: Simon J. Murphy, Timothy R. Bedding, Anuj Gautam, Meridith Joyce

Abstract: The rapidly increasing number of delta Scuti stars with regular patterns among their pulsation frequencies necessitates modelling tools to better understand the observations. Further, with a dozen identified modes per star, there is potential to make meaningful inferences on stellar structure using these young $δ$ Sct stars. We compute and describe a grid of $>$200,000 stellar models from the earl… ▽ More The rapidly increasing number of delta Scuti stars with regular patterns among their pulsation frequencies necessitates modelling tools to better understand the observations. Further, with a dozen identified modes per star, there is potential to make meaningful inferences on stellar structure using these young $δ$ Sct stars. We compute and describe a grid of $>$200,000 stellar models from the early pre-main-sequence to roughly one-third of the main-sequence lifetime, and calculate their pulsation frequencies. From these, we also calculate asteroseismic parameters and explore how those parameters change with mass, age, and metal mass fraction. We show that the large frequency separation, $Δν$, is insensitive to mass at the zero-age main sequence. In the frequency regime observed, the $Δν$ we measure (from modes with $n\sim5$--9) differs from the solar scaling relation by $\sim$13%. We find that the lowest radial order is often poorly modelled, perhaps indicating that the lower-order pressure modes contain further untapped potential for revealing the physics of the stellar interior. We also show that different nuclear reaction networks available in MESA can affect the pulsation frequencies of young $δ$ Sct stars by as much as 5%. We apply the grid to five newly modelled stars, including two pre-main-sequence stars each with 15+ modes identified, and we make the grid available as a community resource. △ Less

Submitted 12 September, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: 17 pages including appendices. Submitted to MNRAS

A multi-wavelength analysis of BL Her stars: Models versus Observations

Authors: S. Das, L. Molnár, S. M. Kanbur, M. Joyce, A. Bhardwaj, H. P. Singh, M. Marconi, V. Ripepi, R. Smolec

Abstract: We present new theoretical period--luminosity (PL) and period--radius (PR) relations at multiple wavelengths (Johnson--Cousins--Glass and {\sl Gaia} passbands) for a fine grid of BL~Herculis models computed using {\sc mesa-rsp}. The non-linear models were computed for periods typical of BL~Her stars, i.e. $1\leq P ({\rm days}) \leq4$, covering a wide range of input parameters: metallicity ($-$2.0… ▽ More We present new theoretical period--luminosity (PL) and period--radius (PR) relations at multiple wavelengths (Johnson--Cousins--Glass and {\sl Gaia} passbands) for a fine grid of BL~Herculis models computed using {\sc mesa-rsp}. The non-linear models were computed for periods typical of BL~Her stars, i.e. $1\leq P ({\rm days}) \leq4$, covering a wide range of input parameters: metallicity ($-$2.0 dex $\leq$ [Fe/H] $\leq$ 0.0 dex), stellar mass (0.5--0.8 M$_{\odot}$), luminosity (50--300 L$_{\odot}$) and effective temperature (full extent of the instability strip; in steps of 50K). We investigate the impact of four sets of convection parameters on multi-wavelength properties. Most empirical relations match well with theoretical relations from the BL~Her models computed using the four sets of convection parameters. No significant metallicity effects are seen in the PR relations. Another important result from our grid of BL~Her models is that it supports combining PL relations of RR Lyrae and Type~II Cepheids together as an alternative to classical Cepheids for the extragalactic distance scale calibration. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 8 pages, 3 figures, Proceedings for IAUS376 'At the cross-roads of astrophysics and cosmology: Period-luminosity relations in the 2020s', Budapest, April 2023

Asteroseismic sounding of bulge globular clusters with the Roman Space Telescope

Authors: László Molnár, Csilla Kalup, Meridith Joyce

Abstract: Globular clusters are relics of the early Universe and they hold clues to many aspects of stellar and galactic evolution. We propose to point the Roman Space Telescope at one or more clusters either as a part of or as an extension of the Galactic Bulge Time Domain Survey. This would provide a unique opportunity to apply the powerful toolkit of asteroseismology to a globular cluster, an observation… ▽ More Globular clusters are relics of the early Universe and they hold clues to many aspects of stellar and galactic evolution. We propose to point the Roman Space Telescope at one or more clusters either as a part of or as an extension of the Galactic Bulge Time Domain Survey. This would provide a unique opportunity to apply the powerful toolkit of asteroseismology to a globular cluster, an observation that is largely out of reach for any other time-domain photometric missions. In this white paper we present the possible targets in the vicinity of the notional survey fields. Potential science cases include precise determination of stellar parameters throughout the cluster, accurate estimation of the integrated mass loss for metal-poor and metal-rich clusters, asteroseismic analysis and mass estimation for RR Lyrae stars, and determination of the seismic ages of clusters. We provide comparisons with other photometric missions and recommendations for maximizing the scientific return from a dedicated globular cluster observing run. △ Less

Submitted 26 June, 2023; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: White paper, 12 pages, 2 figures. Submitted to the Roman Core Community Surveys White Paper Call. v2: minor updates

Comment on the feasibility of carbon burning in Betelgeuse: a response to "The evolutionary stage of Betelgeuse inferred from its pulsation periods," arXiv:2306.00287

Authors: László Molnár, Meridith Joyce, Shing-Chi Leung

Abstract: The recent pre-print by Saio et al. 2023 argues that the supergiant Betelgeuse is already undergoing carbon burning, based on the assumption that all of its light variations are caused by radial pulsations. However, the angular diameter measurements of the star are in conflict with the stellar radius required by their models, as we show in this note. We discuss the feasibility that the Great Dimmi… ▽ More The recent pre-print by Saio et al. 2023 argues that the supergiant Betelgeuse is already undergoing carbon burning, based on the assumption that all of its light variations are caused by radial pulsations. However, the angular diameter measurements of the star are in conflict with the stellar radius required by their models, as we show in this note. We discuss the feasibility that the Great Dimming was caused by constructive mode interference using long-term brightness measurements and comment on differences in modeling frameworks adopted in Saio et al. 2023 vs Joyce et al. 2020. △ Less

Submitted 9 June, 2023; originally announced June 2023.

Comments: This note is a response to the pre-print arXiv:2306.00287. We present it in the spirit of community discourse. 3 pages, 1 figure, accepted to RNAAS. Feedback welcome

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

Testing growth rate dependence in cosmological perturbation theory using scale-free models

Authors: Azrul Pohan, Michael Joyce, David Benhaiem, Francesco Sylos Labini

Abstract: We generalize previously derived analytic results for the one-loop power spectrum (PS) in scale-free models (with linear PS $P(k) \propto k^n$) to a broader class of such models in which part of the matterlike component driving the Einstein de Sitter expansion does not cluster. These models can be conveniently parametrized by $α$, the constant logarithmic linear growth rate of fluctuations (with… ▽ More We generalize previously derived analytic results for the one-loop power spectrum (PS) in scale-free models (with linear PS $P(k) \propto k^n$) to a broader class of such models in which part of the matterlike component driving the Einstein de Sitter expansion does not cluster. These models can be conveniently parametrized by $α$, the constant logarithmic linear growth rate of fluctuations (with $α=1$ in the usual case). For $-3< n<-1$, where the one-loop PS is both infrared and ultraviolet convergent and thus explicitly self-similar, it is characterized conveniently by a single numerical coefficient $c(n, α)$. We compare the analytical predictions for $c(n=-2, α)$ with results from a suite of $N$-body simulations with $α\in [0.25, 1]$ performed with an appropriately modified version of the GADGET code. Although the simulations are of small ($256^3$) boxes, the constraint of self-similarity allows the identification of the converged PS at a level of accuracy sufficient to test the analytical predictions for the $α$ dependence of the evolved PS. Good agreement for the predicted dependence on $α$ of the PS is found. To treat the UV sensitivity of results which grows as one approaches $n =-1$, we derive exact results incorporating a regularization $k_c$ and obtain expressions for $c(n, α, k_c/k)$. Assuming that this regularization is compatible with self-similarity allows us to infer a predicted functional form of the PS equivalent to that derived in effective field theory (EFT). The coefficient of the leading EFT correction at one loop has a strong dependence on $α$, with a change in sign at $α\approx 0.16$, providing a potentially stringent test of EFT. △ Less

Submitted 28 July, 2023; v1 submitted 24 April, 2023; originally announced April 2023.

Comments: revised version published in Phys. Rev. D

A Review of the Mixing Length Theory of Convection in 1D Stellar Modeling

Authors: Meridith Joyce, Jamie Tayar

Abstract: We review the application of the one-dimensional Mixing Length Theory (MLT) model of convection in stellar interiors and low-mass stellar evolution. We summarize the history of MLT, present a derivation of MLT in the context of the 1D stellar structure equations, and discuss the physical regimes in which MLT is relevant. We review attempts to improve and extend the formalism, including to higher d… ▽ More We review the application of the one-dimensional Mixing Length Theory (MLT) model of convection in stellar interiors and low-mass stellar evolution. We summarize the history of MLT, present a derivation of MLT in the context of the 1D stellar structure equations, and discuss the physical regimes in which MLT is relevant. We review attempts to improve and extend the formalism, including to higher dimensions. We discuss the interactions of MLT with other modeling physics and demonstrate the impact of introducing variations in the convective mixing length, αMLT, on stellar tracks and isochrones. We summarize the process of performing a solar calibration of αMLT and the state-of-the-art on calibrations to non-solar targets. We discuss the scientific implications of changing the mixing length, using recent analyses as demonstration. We review the most prominent successes of MLT and remaining challenges, and we conclude by speculating on the future of this treatment of convection. △ Less

Submitted 31 May, 2023; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: invited review for MDPI Galaxies Special Issue "The Structure and Evolution of Stars" (resubmitted) Generally minor revisions in response to referee reports and community feedback. Improvements to Figure 1 and derivation in Section 6

Constraining accuracy of the pairwise velocities in $N$-body simulations using scale-free models

Authors: Sara Maleubre, Daniel J. Eisenstein, Lehman H. Garrison, Michael Joyce

Abstract: We present a continuation of an analysis that aims to quantify resolution of $N$-body simulations by exploiting large (up to $N=4096^3$) simulations of scale-free cosmologies run using Abacus. Here we focus on radial pairwise velocities of the matter field, both by direct estimation and through the cumulative-2PCF (using the pair conservation equation). We find that convergence at the $1\%$ level… ▽ More We present a continuation of an analysis that aims to quantify resolution of $N$-body simulations by exploiting large (up to $N=4096^3$) simulations of scale-free cosmologies run using Abacus. Here we focus on radial pairwise velocities of the matter field, both by direct estimation and through the cumulative-2PCF (using the pair conservation equation). We find that convergence at the $1\%$ level of the mean relative pairwise velocity can be demonstrated over a range of scales, evolving from a few times the grid spacing at early times to slightly below this scale at late times. We show the analysis of two different box sizes as well as from averaging results from the smaller boxes, and compare the power of the two aforementioned estimators in constraining accuracy at each scale. Down to scales of order of the smoothing parameter, convergence is obtained at $\sim5\%$ precision, and shows a behaviour indicating asymptotic stable clustering. We also infer for LCDM simulations conservative estimates on the evolution of the lower cut-off to resolution (at $1\%$ and $5\%$ precision) as a function of redshift. △ Less

Submitted 2 August, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 14 pages, 9 figures, 1 table. Accepted for publication in MNRAS. v2 does NOT overlap with arXiv:2308.00438

JWST Imaging of Earendel, the Extremely Magnified Star at Redshift $z=6.2$

Authors: Brian Welch, Dan Coe, Erik Zackrisson, S. E. de Mink, Swara Ravindranath, Jay Anderson, Gabriel Brammer, Larry Bradley, Jinmi Yoon, Patrick Kelly, Jose M. Diego, Rogier Windhorst, Adi Zitrin, Paola Dimauro, Yolanda Jimenez-Teja, Abdurro'uf, Mario Nonino, Ana Acebron, Felipe Andrade-Santos, Roberto J. Avila, Matthew B. Bayliss, Alex Benitez, Tom Broadhurst, Rachana Bhatawdekar, Marusa Bradac , et al. (38 additional authors not shown)

Abstract: The gravitationally lensed star WHL0137-LS, nicknamed Earendel, was identified with a photometric redshift $z_{phot} = 6.2 \pm 0.1$ based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera (NIRCam) images of Earendel in 8 filters spanning 0.8--5.0$μ$m. In these higher resolution images, Earendel remains a single unresolved point… ▽ More The gravitationally lensed star WHL0137-LS, nicknamed Earendel, was identified with a photometric redshift $z_{phot} = 6.2 \pm 0.1$ based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera (NIRCam) images of Earendel in 8 filters spanning 0.8--5.0$μ$m. In these higher resolution images, Earendel remains a single unresolved point source on the lensing critical curve, increasing the lower limit on the lensing magnification to $μ> 4000$ and restricting the source plane radius further to $r < 0.02$ pc, or $\sim 4000$ AU. These new observations strengthen the conclusion that Earendel is best explained by an individual star or multiple star system, and support the previous photometric redshift estimate. Fitting grids of stellar spectra to our photometry yields a stellar temperature of $T_{\mathrm{eff}} \simeq 13000$--16000 K assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from $\log(L) = 5.8$--6.6 $L_{\odot}$, which is in the range where one expects luminous blue variable stars. Follow-up observations, including JWST NIRSpec scheduled for late 2022, are needed to further unravel the nature of this object, which presents a unique opportunity to study massive stars in the first billion years of the universe. △ Less

Submitted 9 November, 2022; v1 submitted 18 August, 2022; originally announced August 2022.

Comments: Accepted to ApJL. Data products, lens models, and analysis code will be available online at https://cosmic-spring.github.io

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

Bridging the Gap between Intermediate and Massive Stars I: Validation of MESA against the State-of-the-Art Monash Stellar Evolution Program for a 2$M_{\odot}$ AGB Star

Authors: Giulia C. Cinquegrana, Meridith Joyce, Amanda I. Karakas

Abstract: One--dimensional stellar structure and evolution programs are built using different physical prescriptions and algorithms, which means there can be variations between models' predictions even when using identical input physics. This leads to questions about whether such deviations are physical or numerical; code validation studies are important and necessary tools for studying these questions. We… ▽ More One--dimensional stellar structure and evolution programs are built using different physical prescriptions and algorithms, which means there can be variations between models' predictions even when using identical input physics. This leads to questions about whether such deviations are physical or numerical; code validation studies are important and necessary tools for studying these questions. We provide the first direct comparison between the Monash stellar evolution program and MESA for a $2M_{\odot}$ model evolved from the zero-age main sequence to the tip of the thermally pulsing asymptotic giant branch. We compare the internal structure of the two models at six critical evolutionary points and find that they are in excellent agreement in characteristics like central temperature, central density and the temperature at the base of the convective envelope during the thermally pulsing asymptotic giant branch. The hydrogen-exhausted core mass between the models differs by less than 4.2% throughout the entire evolution, the final values vary only by 1.5%. Surface quantities such as luminosity and radius vary by less than 0.2% prior to the asymptotic giant branch. During thermal pulses, the difference extends to 3.4%, largely due to uncertainties in mixing and the treatment of atmospheric boundary conditions. Given that the veteran Monash code is closed source, the present work provides the first fully open-source computational analog. This increases accessibility to precision modeling on the asymptotic giant branch and lays the groundwork for higher-mass calculations that are performed with MESA but preserve the standards of the Monash code during the AGB. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: Submitted, ApJ

Journal ref: The Astrophysical Journal 939.1(2022):50

Spinning up the Surface: Evidence for Planetary Engulfment or Unexpected Angular Momentum Transport?

Authors: Jamie Tayar, Facundo D. Moyano, Melinda Soares-Furtado, Ana Escorza, Meridith Joyce, Sarah L. Martell, Rafael A. García, Sylvain N. Breton, Stéphane Mathis, Savita Mathur, Vincent Delsanti, Sven Kiefer, Sabine Reffert, Dominic M. Bowman, Timothy Van Reeth, Shreeya Shetye, Charlotte Gehan, Samuel K. Grunblatt

Abstract: In this paper, we report the potential detection of a nonmonotonic radial rotation profile in a low-mass lower-luminosity giant star. For most low- and intermediate-mass stars, the rotation on the main sequence seems to be close to rigid. As these stars evolve into giants, the core contracts and the envelope expands, which should suggest a radial rotation profile with a fast core and a slower enve… ▽ More In this paper, we report the potential detection of a nonmonotonic radial rotation profile in a low-mass lower-luminosity giant star. For most low- and intermediate-mass stars, the rotation on the main sequence seems to be close to rigid. As these stars evolve into giants, the core contracts and the envelope expands, which should suggest a radial rotation profile with a fast core and a slower envelope and surface. KIC 9267654, however, seems to show a surface rotation rate that is faster than its bulk envelope rotation rate, in conflict with this simple angular momentum conservation argument. We improve the spectroscopic surface constraint, show that the pulsation frequencies are consistent with the previously published core and envelope rotation rates, and demonstrate that the star does not show strong chemical peculiarities. We discuss the evidence against any tidally interacting stellar companion. Finally, we discuss the possible origin of this unusual rotation profile, including the potential ingestion of a giant planet or unusual angular momentum transport by tidal inertial waves triggered by a close substellar companion, and encourage further observational and theoretical efforts. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: 18 pages, 9 figures, submitted to AAS Journals

A prescription for the asteroseismic surface correction

Authors: Yaguang Li, Timothy R. Bedding, Dennis Stello, Daniel Huber, Marc Hon, Meridith Joyce, Tanda Li, Jean Perkins, Timothy R. White, Joel C. Zinn, Andrew W. Howard, Howard Isaacson, Daniel R. Hey, Hans Kjeldsen

Abstract: In asteroseismology, the surface effect refers to a disparity between the observed and the modelled frequencies in stars with solar-like oscillations. It originates from improper modelling of the surface layers. Correcting the surface effect usually requires using functions with free parameters, which are conventionally fitted to the observed frequencies. On the basis that the correction should va… ▽ More In asteroseismology, the surface effect refers to a disparity between the observed and the modelled frequencies in stars with solar-like oscillations. It originates from improper modelling of the surface layers. Correcting the surface effect usually requires using functions with free parameters, which are conventionally fitted to the observed frequencies. On the basis that the correction should vary smoothly across the H--R diagram, we parameterize it as a simple function of surface gravity, effective temperature, and metallicity. We determine this function by fitting a wide range of stars. The absolute amount of the surface correction decreases with luminosity, but the ratio between it and $ν_{\rm max}$ increases, suggesting the surface effect is more important for red giants than dwarfs. Applying the prescription can eliminate unrealistic surface correction, which improves parameter estimations with stellar modelling. Using two open clusters, we found a reduction of scatter in the model-derived ages for each star in the same cluster. As an important application, we provide a new revision for the $Δν$ scaling relation that, for the first time, accounts for the surface correction. The values of the correction factor, $f_{Δν}$, are up to 2\% smaller than those determined without the surface effect considered, suggesting decreases of up to 4\% in radii and up to 8\% in masses when using the asteroseismic scaling relations. This revision brings the asteroseismic properties into an agreement with those determined from eclipsing binaries. The new correction factor and the stellar models with the corrected frequencies are available at {https://www.github.com/parallelpro/surface}. △ Less

Submitted 15 May, 2023; v1 submitted 1 August, 2022; originally announced August 2022.

Comments: 13 pages, 11 figures. Accepted by MNRAS

Cosmological perturbation theory using generalized Einstein de Sitter cosmologies

Authors: Michael Joyce, Azrul Pohan

Abstract: The separable analytical solution in standard perturbation theory for an Einstein de Sitter (EdS) universe can be generalized to the wider class of such cosmologies (``generalized EdS'', or gEdS) in which a fraction of the pressure-less fluid does not cluster. We derive the corresponding kernels in both Eulerian perturbation theory (EPT) and Lagrangian perturbation theory, generalizing the canonic… ▽ More The separable analytical solution in standard perturbation theory for an Einstein de Sitter (EdS) universe can be generalized to the wider class of such cosmologies (``generalized EdS'', or gEdS) in which a fraction of the pressure-less fluid does not cluster. We derive the corresponding kernels in both Eulerian perturbation theory (EPT) and Lagrangian perturbation theory, generalizing the canonical EdS expressions to a one-parameter family where the parameter can be taken to be the exponent $α$ of the growing mode linear amplification $D(a) \propto a^α$. For the power spectrum (PS) at one loop in EPT, the contribution additional to standard EdS is given, for each of the `13' and `22' terms, as a function of two infra-red safe integrals. In the second part of the paper we show that the calculation of cosmology-dependent corrections in perturbation theory in standard (e.g. LCDM-like) models can be simplified, and their magnitude and parameter dependence better understood, by relating them to our analytic results for gEdS models. At second order the time dependent kernels are equivalent to the analytic kernels of the gEdS model with $α$ replaced by a single redshift dependent effective growth rate $α_2(z)$. At third order the time evolution can be conveniently parametrized in terms of two additional such effective growth rates. For the PS calculated at one loop order, the correction to the PS relative to the EdS limit can be expressed in terms of just $α_2(z)$, one additional effective growth rate function and the four infra-red safe integrals of the gEdS limit. This is much simplified compared to expressions in the literature that use six or eight red-shift dependent functions and are not explicitly infra-red safe. Using the analytic gEdS expression for the PS with $α=α_2(z)$ gives a good approximation (to $\sim 25 \%$) for the exact result. △ Less

Submitted 17 April, 2023; v1 submitted 10 July, 2022; originally announced July 2022.

Comments: revised version accepted in Phys. Rev. D., additional detailed comparison with previous works, new summary table, additional references

Blanco DECam Bulge Survey (BDBS) IV: Metallicity Distributions and Bulge Structure from 2.6 Million Red Clump Stars

Authors: Christian I. Johnson, R. Michael Rich, Iulia T. Simion, Michael D. Young, William I. Clarkson, Catherine A. Pilachowski, Scott Michael, Tommaso Marchetti, Mario Soto, Andrea Kunder, Andreas J. Koch-Hansen, A. Katherina Vivas, Meridith Joyce, Juntai Shen, Alexis Osmond

Abstract: We present photometric metallicity measurements for a sample of 2.6 million bulge red clump stars extracted from the Blanco DECam Bulge Survey (BDBS). Similar to previous studies, we find that the bulge exhibits a strong vertical metallicity gradient, and that at least two peaks in the metallicity distribution functions appear at b < -5. We can discern a metal-poor ([Fe/H] ~ -0.3) and metal-rich (… ▽ More We present photometric metallicity measurements for a sample of 2.6 million bulge red clump stars extracted from the Blanco DECam Bulge Survey (BDBS). Similar to previous studies, we find that the bulge exhibits a strong vertical metallicity gradient, and that at least two peaks in the metallicity distribution functions appear at b < -5. We can discern a metal-poor ([Fe/H] ~ -0.3) and metal-rich ([Fe/H] ~ +0.2) abundance distribution that each show clear systematic trends with latitude, and may be best understood by changes in the bulge's star formation/enrichment processes. Both groups exhibit asymmetric tails, and as a result we argue that the proximity of a star to either peak in [Fe/H] space is not necessarily an affirmation of group membership. The metal-poor peak shifts to lower [Fe/H] values at larger distances from the plane while the metal-rich tail truncates. Close to the plane, the metal-rich tail appears broader along the minor axis than in off-axis fields. We also posit that the bulge has two metal-poor populations -- one that belongs to the metal-poor tail of the low latitude and predominantly metal-rich group, and another belonging to the metal-poor group that dominates in the outer bulge. We detect the X-shape structure in fields with |Z| > 0.7 kpc and for stars with [Fe/H] > -0.5. Stars with [Fe/H] < -0.5 may form a spheroidal or "thick bar" distribution while those with [Fe/H] > -0.1 are strongly concentrated near the plane. △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: 26 pages, 22 figures, accepted for publication in MNRAS; the full data table is very long so only a stub table has been provided here; the full electronic table will be provided through MNRAS upon publication, but early access to the full table will be granted upon request to the authors

Blanco DECam Bulge Survey (BDBS) V. Cleaning the foreground populations from Galactic bulge colour-magnitude diagrams using Gaia EDR3

Authors: Tommaso Marchetti, Christian I. Johnson, Meridith Joyce, R. Michael Rich, Iulia Simion, Michael D. Young, William I. Clarkson, Catherine A. Pilachowski, Scott Michael, Andrea Kunder, Andreas J. Koch-Hansen

Abstract: The Blanco DECam Bulge Survey (BDBS) has imaged more than 200 square degrees of the southern Galactic bulge, providing photometry in the ugrizy filters for $\sim 250$ million unique stars. The presence of a strong foreground disk population, along with complex reddening and extreme image crowding, has made it difficult to constrain the presence of young and intermediate age stars in the bulge popu… ▽ More The Blanco DECam Bulge Survey (BDBS) has imaged more than 200 square degrees of the southern Galactic bulge, providing photometry in the ugrizy filters for $\sim 250$ million unique stars. The presence of a strong foreground disk population, along with complex reddening and extreme image crowding, has made it difficult to constrain the presence of young and intermediate age stars in the bulge population. We employed an accurate cross-match of BDBS with the latest data release (EDR3) from the Gaia mission, matching more than 140 million sources with BDBS photometry and Gaia EDR3 photometry and astrometry. We relied on Gaia EDR3 astrometry, without any photometric selection, to produce clean BDBS bulge colour-magnitude diagrams (CMDs). Gaia parallaxes were used to filter out bright foreground sources, and a Gaussian mixture model fit to Galactic proper motions could identify stars kinematically consistent with bulge membership. We applied this method to 127 different bulge fields of $1$ deg$^2$ each, with $|\ell| \leq 9.5^\circ$ and $-9.5^\circ \leq b \leq -2.5^\circ$. The astrometric cleaning procedure removes the majority of blue stars in each field, especially near the Galactic plane, where the ratio of blue to red stars is $\lesssim 10\%$, increasing to values $\sim 20\%$ at higher Galactic latitudes. We rule out the presence of a widespread population of stars younger than 2 Gyr. The vast majority of blue stars brighter than the turnoff belong to the foreground population, according to their measured astrometry. We introduce the distance between the observed red giant branch bump and the red clump as a simple age proxy for the dominant population in the field, and we confirm the picture of a predominantly old bulge. Further work is needed to apply the method to estimate ages to fields at higher latitudes, and to model the complex morphology of the Galactic bulge. △ Less

Submitted 26 June, 2022; originally announced June 2022.

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

Journal ref: A&A 664, A124 (2022)

Gender Disparity in Publishing Six Months after the KITP Workshop "Probes of Transport in Stars"

Authors: Meridith Joyce, Jamie Tayar, Daniel Lecoanet

Abstract: Conferences and workshops shape scientific discourse. The Kavli Institute for Theoretical Physics (KITP) hosts long-term workshops to stimulate scientific collaboration that would not otherwise have taken place. One goal of KITP programs is to increase diversity in the next generation of scientists. In this analysis, we examine gender trends in authorship of papers generated as a result of the KIT… ▽ More Conferences and workshops shape scientific discourse. The Kavli Institute for Theoretical Physics (KITP) hosts long-term workshops to stimulate scientific collaboration that would not otherwise have taken place. One goal of KITP programs is to increase diversity in the next generation of scientists. In this analysis, we examine gender trends in authorship of papers generated as a result of the KITP program \textit{Probes of Transport in Stars}, which ran from October 11th, 2021 to December 17th, 2021. While 38\% of workshop participants were women, only 19\% of publications produced between December 1st, 2021 and June 3rd, 2022 had female first-authors. Further, of these early publications, 61\% had all-male author lists. Among publications resulting from the KITP program, the proportions of both male first-author papers and papers with all-male author lists are higher than predicted by models that take into account the gender distribution of the KITP participants. These results motivate more thorough investigations of collaboration networks at scientific conferences and workshops. Importantly, they also suggest that programs, conferences, and workshops of any kind need to take steps beyond those implemented in this KITP program to enable more diverse collaborations and address gender disparities in science. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: submitted to PASP. Comments and discussion welcome

Is Thermohaline Mixing the Full Story? Evidence for Separate Mixing Events near the Red Giant Branch Bump

Authors: Jamie Tayar, Meridith Joyce

Abstract: The abundances of mixing--sensitive elements including lithium, [C/N], and 12C/13C are known to change near the red giant branch bump. The explanation most often offered for these alterations is double diffusive thermohaline mixing in the stellar interior. In this analysis, we investigate the ability of thermohaline mixing to explain the observed timing of these chemical depletion events. Recent o… ▽ More The abundances of mixing--sensitive elements including lithium, [C/N], and 12C/13C are known to change near the red giant branch bump. The explanation most often offered for these alterations is double diffusive thermohaline mixing in the stellar interior. In this analysis, we investigate the ability of thermohaline mixing to explain the observed timing of these chemical depletion events. Recent observational measurements of lithium and [C/N] show that the abundance of lithium decreases before the abundance of [C/N], whereas numerical simulations of the propagation of the thermohaline mixing region computed with MESA show that the synthetic abundances drop simultaneously. We therefore conclude that thermohaline mixing alone cannot explain the distinct events of lithium depletion and [C/N] depletion, as the simultaneity predicted by simulations is not consistent with the observation of separate drops. We thus invite more sophisticated theoretical explanations for the observed temporal separation of these chemical depletion episodes as well as more extensive observational explorations across a range of masses and metallicities. △ Less

Submitted 19 May, 2022; originally announced May 2022.

Comments: In revision at AAS Journals. 9 pages, 6 figures

The Ages of Galactic Bulge Stars with Realistic Uncertainties

Authors: Meridith Joyce, Christian I. Johnson, Tommaso Marchetti, R. Michael Rich, Iulia Simion, John Bourke

Abstract: Using modern isochrones with customized physics and carefully considered statistical techniques, we recompute the age distribution for a sample of 91 micro-lensed dwarfs in the Galactic bulge presented by Bensby et al. (2017) and do not produce an age distribution consistent with their results. In particular, our analysis finds that only 15 of 91 stars have ages younger than 7 Gyr, compared to the… ▽ More Using modern isochrones with customized physics and carefully considered statistical techniques, we recompute the age distribution for a sample of 91 micro-lensed dwarfs in the Galactic bulge presented by Bensby et al. (2017) and do not produce an age distribution consistent with their results. In particular, our analysis finds that only 15 of 91 stars have ages younger than 7 Gyr, compared to their finding of 42 young stars in the same sample. While we do not find a constituency of very young stars, our results do suggest the presence of an $\sim8$ Gyr population at the highest metallicities, thus contributing to long-standing debate about the age--metallicity distribution of the Galactic bulge. We supplement this with attempts at independent age determinations from two sources of photometry, BDBS and \textit{Gaia}, but find that the imprecision of photometric measurements prevents reliable age and age uncertainty determinations. Lastly, we present age uncertainties derived using a first-order consideration of global modeling uncertainties in addition to standard observational uncertainties. The theoretical uncertainties are based on the known variance of free parameters in the 1D stellar evolution models used to generate isochrones, and when included, result in age uncertainties of $2$--$5$ Gyr for this spectroscopically well-constrained sample. These error bars, which are roughly twice as large as typical literature values, constitute realistic lower limits on the true age uncertainties. △ Less

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

Comments: accepted to ApJ; revisions complete

Stellar and black hole assembly in z<0.3 infrared-luminous mergers: intermittent starbursts vs. super-Eddington accretion

Authors: Duncan Farrah, Andreas Efstathiou, Jose Afonso, Jeronimo Bernard-Salas, Joe Cairns, David L Clements, Kevin Croker, Evanthia Hatziminaoglou, Maya Joyce, Mark Lacy, Vianney Lebouteiller, Alix Lieblich, Carol Lonsdale, Seb Oliver, Chris Pearson, Sara Petty, Lura K Pitchford, Dimitra Rigopoulou, Michael Rowan-Robinson, Jack Runburg, Henrik Spoon, Aprajita Verma, Lingyu Wang

Abstract: We study stellar and black hole mass assembly in a sample of 42 infrared-luminous galaxy mergers at z<0.3 by combining results from radiative transfer modelling with archival measures of molecular gas and black hole mass. The ratios of stellar mass, molecular gas mass, and black hole mass to each other are consistent with those of massive gas-rich galaxies at z<0.3. The advanced mergers may show i… ▽ More We study stellar and black hole mass assembly in a sample of 42 infrared-luminous galaxy mergers at z<0.3 by combining results from radiative transfer modelling with archival measures of molecular gas and black hole mass. The ratios of stellar mass, molecular gas mass, and black hole mass to each other are consistent with those of massive gas-rich galaxies at z<0.3. The advanced mergers may show increased black hole mass to stellar mass ratios, consistent with the transition from AGN to ellipticals and implying substantial black hole mass growth over the course of the merger. Star formation rates are enhanced relative to the local main sequence, by factors of ~100 in the starburst and ~1.8 in the host. The starburst star formation rates appear distinct to star formation in the main sequence at all redshifts up to at least z~5. Starbursts may prefer late-stage mergers, but are observed at any merger stage. We do not find evidence that the starbursts in these low-redshift systems substantially increase the total stellar mass, with a soft upper limit on the stellar mass increase from starburst activity of about a factor of two. In contrast, 12 objects show evidence for super-Eddington accretion, associated with late-stage mergers, suggesting that many AGN in infrared-luminous mergers go through a super-Eddington phase. The super-Eddington phase may increase black hole mass by up to an order of magnitude at an accretion efficiency of 42+/-33% over a period of 44+/-22Myr. Our results imply that super-Eddington accretion is an important black hole growth channel in infrared-luminous galaxies at all redshifts. △ Less

Submitted 29 April, 2022; originally announced May 2022.

Comments: MNRAS Accepted

Solar calibration of $α_{\rm MLT}$ using ÆSOPUS opacities in MESA

Authors: Giulia C Cinquegrana, Meridith Joyce

Abstract: The simplistic but ubiquitous Mixing Length Theory (MLT) formalism is used to model convective energy transport within 1D stellar evolution calculations. The formalism relies on the free parameter $α_{\rm MLT}$, which must be independently calibrated within each stellar evolution program and for any given set of physical assumptions. We present a solar calibration of $α_{\text{MLT}}$ appropriate f… ▽ More The simplistic but ubiquitous Mixing Length Theory (MLT) formalism is used to model convective energy transport within 1D stellar evolution calculations. The formalism relies on the free parameter $α_{\rm MLT}$, which must be independently calibrated within each stellar evolution program and for any given set of physical assumptions. We present a solar calibration of $α_{\text{MLT}}$ appropriate for use with the AESOPUS opacities, which have recently been made available for use with the MESA stellar evolution software. We report a calibrated value of $α_{\rm MLT}=1.931$ and demonstrate the impact of using an appropriately calibrated value in simulations of a $3 M_{\odot}$ asymptotic giant branch star. △ Less

Submitted 18 April, 2022; originally announced April 2022.

Comments: 4 pages, 1 figure. Accepted for publication in RNAAS

Journal ref: Research Notes of the AAS 6.4 (2022): 77

Characterizing Observed Extra Mixing Trends in Red Giants using the Reduced Density Ratio from Thermohaline Models

Authors: Adrian E. Fraser, Meridith Joyce, Evan H. Anders, Jamie Tayar, Matteo Cantiello

Abstract: Observations show an almost ubiquitous presence of extra mixing in low-mass upper giant branch stars. The most commonly invoked explanation for this is thermohaline mixing. One-dimensional stellar evolution models include various prescriptions for thermohaline mixing, but the use of observational data directly to discriminate between thermohaline prescriptions has thus far been limited. Here, we p… ▽ More Observations show an almost ubiquitous presence of extra mixing in low-mass upper giant branch stars. The most commonly invoked explanation for this is thermohaline mixing. One-dimensional stellar evolution models include various prescriptions for thermohaline mixing, but the use of observational data directly to discriminate between thermohaline prescriptions has thus far been limited. Here, we propose a new framework to facilitate direct comparison: Using carbon-to-nitrogen measurements from the SDSS-IV APOGEE survey as a probe of mixing and a fluid parameter known as the reduced density ratio from one-dimensional stellar evolution programs, we compare the observed amount of extra mixing on the upper giant branch to predicted trends from three-dimensional fluid dynamics simulations. Using this method, we are able to empirically constrain how mixing efficiency should vary with the reduced density ratio. We find the observed amount of extra mixing is strongly correlated with the reduced density ratio and that trends between reduced density ratio and fundamental stellar parameters are robust across choices for modeling prescription. We show that stars with available mixing data tend to have relatively low density ratios, which should inform the regimes selected for future simulation efforts. Finally, we show that there is increased mixing at low reduced density ratios, which is consistent with current hydrodynamical models of thermohaline mixing. The introduction of this framework sets a new standard for theoretical modeling efforts, as validation for not only the amount of extra mixing, but trends between the degree of extra mixing and fundamental stellar parameters is now possible. △ Less

Submitted 11 October, 2022; v1 submitted 18 April, 2022; originally announced April 2022.

Comments: 22 pages, 7 figures, accepted for publication in ApJ. MESA histories available on Zenodo at https://doi.org/10.5281/zenodo.7152746 . MESA inlists available on GitHub at https://github.com/afraser3/Empirical-Magnetic-Thermohaline

Discovery of post-mass-transfer helium-burning red giants using asteroseismology

Authors: Yaguang Li, Timothy R. Bedding, Simon J. Murphy, Dennis Stello, Yifan Chen, Daniel Huber, Meridith Joyce, Dion Marks, Xianfei Zhang, Shaolan Bi, Isabel L. Colman, Michael R. Hayden, Daniel R. Hey, Gang Li, Benjamin T. Montet, Sanjib Sharma, Yaqian Wu

Abstract: A star expands to become a red giant when it has fused all the hydrogen in its core into helium. If the star is in a binary system, its envelope can overflow onto its companion or be ejected into space, leaving a hot core and potentially forming a subdwarf-B star. However, most red giants that have partially transferred envelopes in this way remain cool on the surface and are almost indistinguisha… ▽ More A star expands to become a red giant when it has fused all the hydrogen in its core into helium. If the star is in a binary system, its envelope can overflow onto its companion or be ejected into space, leaving a hot core and potentially forming a subdwarf-B star. However, most red giants that have partially transferred envelopes in this way remain cool on the surface and are almost indistinguishable from those that have not. Among $\sim$7000 helium-burning red giants observed by NASA's Kepler mission, we use asteroseismology to identify two classes of stars that must have undergone dramatic mass loss, presumably due to stripping in binary interactions. The first class comprises about 7 underluminous stars with smaller helium-burning cores than their single-star counterparts. Theoretical models show that these small cores imply the stars had much larger masses when ascending the red giant branch. The second class consists of 32 red giants with masses down to 0.5 M$_\odot$, whose implied ages would exceed the age of the universe had no mass loss occurred. The numbers are consistent with binary statistics, and our results open up new possibilities to study the evolution of post-mass-transfer binary systems. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: a Letter to Nature Astronomy

The X-ray Disk/Wind Degeneracy in AGN

Authors: M. L. Parker, G. A. Matzeu, J. H. Matthews, M. J. Middleton, T. Dauser, J. Jiang, A. M. Joyce

Abstract: Relativistic Fe K emission lines from accretion disks and from disk winds encode key information about black holes, and their accretion and feedback mechanisms. We show that these two processes can in principle produce indistinguishable line profiles, such that they cannot be disentangled spectrally. We argue that it is likely that in many cases both processes contribute to the net line profile, a… ▽ More Relativistic Fe K emission lines from accretion disks and from disk winds encode key information about black holes, and their accretion and feedback mechanisms. We show that these two processes can in principle produce indistinguishable line profiles, such that they cannot be disentangled spectrally. We argue that it is likely that in many cases both processes contribute to the net line profile, and their relative contributions cannot be constrained purely by Fe K spectroscopy. In almost all studies of Fe K emission to date, a single process (either disk reflection or wind Compton scattering) is assumed to dominate the total line profile. We demonstrate that fitting a single process emission model (pure reflection or pure wind) to a hybrid line profile results in large systematic biases in the estimates of key parameters, such as mass outflow rate and spin. We discuss various strategies to mitigate this effect, such as including high energy data covering the Compton hump, and the implications for future X-ray missions. △ Less

Submitted 28 March, 2022; originally announced March 2022.

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

Schwarzschild and Ledoux are equivalent on evolutionary timescales

Authors: Evan H. Anders, Adam S. Jermyn, Daniel Lecoanet, Adrian E. Fraser, Imogen G. Cresswell, Meridith Joyce, J. R. Fuentes

Abstract: Stellar evolution models calculate convective boundaries using either the Schwarzschild or Ledoux criterion, but confusion remains regarding which criterion to use. Here we present a 3D hydrodynamical simulation of a convection zone and adjacent radiative zone, including both thermal and compositional buoyancy forces. As expected, regions which are unstable according to the Ledoux criterion are co… ▽ More Stellar evolution models calculate convective boundaries using either the Schwarzschild or Ledoux criterion, but confusion remains regarding which criterion to use. Here we present a 3D hydrodynamical simulation of a convection zone and adjacent radiative zone, including both thermal and compositional buoyancy forces. As expected, regions which are unstable according to the Ledoux criterion are convective. Initially, the radiative zone adjacent to the convection zone is Schwarzschild-unstable but Ledoux-stable due to a composition gradient. Over many convective overturn timescales the convection zone grows via entrainment. The convection zone saturates at the size originally predicted by the Schwarzschild criterion, although in this final state the Schwarzschild and Ledoux criteria agree. Therefore, the Schwarzschild criterion should be used to determine the size of stellar convection zones, except possibly during short-lived evolutionary stages in which entrainment persists. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: Manuscript accepted in ApJ Letters. Supplemental materials are in a Zenodo repository (https://zenodo.org/record/6344868). Code is in a Github repository (https://github.com/evanhanders/schwarzschild_or_ledoux)

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

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

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

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

Comments: Submitted to MNRAS

The most metal-rich asymptotic giant branch stars

Authors: Amanda Karakas, Giulia Cinquegrana, Meridith Joyce

Abstract: We present new stellar evolutionary sequences of very metal-rich stars evolved with the Monash Stellar Structure code and with MESA. The Monash models include masses of $1-8M_{\odot}$ with metallicities $Z=0.04$ to $Z=0.1$ and are evolved from the main sequence to the thermally-pulsing asymptotic giant branch (AGB). These are the first $Z=0.1$ AGB models in the literature. The MESA models include… ▽ More We present new stellar evolutionary sequences of very metal-rich stars evolved with the Monash Stellar Structure code and with MESA. The Monash models include masses of $1-8M_{\odot}$ with metallicities $Z=0.04$ to $Z=0.1$ and are evolved from the main sequence to the thermally-pulsing asymptotic giant branch (AGB). These are the first $Z=0.1$ AGB models in the literature. The MESA models include intermediate-mass models with $Z=0.06$ to $Z=0.09$ evolved to the onset of the thermally-pulsing phase. Third dredge-up only occurs in intermediate-mass models $Z \le 0.08$. Hot bottom burning (HBB) shows a weaker dependence on metallicity, with the minimum mass increasing from 4.5$M_{\odot}$ for $Z=0.014$ to $\approx 5.5 M_{\odot}$ for Z = 0.04, $6M_{\odot}$ for $ 0.05 \le Z \le 0.07$ and above 6.5$ M_{\odot}$ for $Z\ge 0.08$. The behaviour of the $Z=0.1$ models is unusual; most do not experience He-shell instabilities owing to rapid mass-loss on the early part of the AGB. Turning off mass-loss produces He-shell instabilities, however thermal pulses are weak and result in no third dredge-up. The minimum mass for carbon ignition is reduced from 8$M_{\odot}$ for $Z=0.04$ to 7$M_{\odot}$ for $Z=0.1$, which implies a reduction in the minimum mass for core-collapse supernovae. MESA models of similarly high metallicity ($Z=0.06 - 0.09$) show the same lowering of the minimum mass for carbon ignition: carbon burning is detected in a $6 M_{\odot}$ model at the highest metallicity ($Z=0.09$) and in all $7 M_{\odot}$ models with $Z \ge 0.06$. This demonstrates robustness of the lowered carbon burning threshold across codes. △ Less

Submitted 1 November, 2021; originally announced November 2021.

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

TIC 308396022: a $δ$ Scuti--$γ$ Doradus hybrid with large-amplitude radial fundamental mode and regular g-mode period spacing

Authors: Tao-Zhi Yang, Zhao-Yu Zuo, Gang Li, Timothy R Bedding, Simon J Murphy, Meridith Joyce

Abstract: We analyse the pulsating behaviour of TIC 308396022 observed by the TESS mission. The star is a high-amplitude $δ$ Sct star (HADS) that shows a very rich amplitude spectrum using the 3-yr light curve. Among these frequencies, the strongest peak of $f_{1}= 13.20362567(12) \rm{d^{-1}}$ is identified as the radial fundamental mode, and we also find the first and second overtones ($f_2$ and $f_3$). In… ▽ More We analyse the pulsating behaviour of TIC 308396022 observed by the TESS mission. The star is a high-amplitude $δ$ Sct star (HADS) that shows a very rich amplitude spectrum using the 3-yr light curve. Among these frequencies, the strongest peak of $f_{1}= 13.20362567(12) \rm{d^{-1}}$ is identified as the radial fundamental mode, and we also find the first and second overtones ($f_2$ and $f_3$). In the low frequency range (< 2.5 $\rm{d^{-1}}$), 22 peaks are identified to be gravity modes, which show a regular period spacing of about 2460\,s and have the angular degree $l = 1$. The period spacing pattern does not show a significant downward trend, suggesting the star rotates slowly. We note that this is a $δ$ Sct--$γ$ Dor hybrid star containing a high-amplitude radial fundamental mode and a regular g-mode period spacing pattern. With O-C analysis, we find the star shows a significant time delay, implying that the star has a companion which is likely to be a white dwarf. The history of possible mass transfer provides a great opportunity to test the current theories of binary evolution, mass transfer, and pulsation. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: Accepted by A&A

First results on RR Lyrae stars with the TESS space telescope: untangling the connections between mode content, colors and distances

Authors: László Molnár, Attila Bódi, András Pál, Anupam Bhardwaj, Franz-Josef Hambsch, József M. Benkő, Aliz Derekas, Mohammad Ebadi, Meridith Joyce, Amir Hasanzadeh, Katrien Kolenberg, Michael B. Lund, James M. Nemec, Henryka Netzel, Chow-Choong Ngeow, Joshua Pepper, Emese Plachy, Zdeněk Prudil, Robert J. Siverd, Marek Skarka, Radosław Smolec, Ádám Sódor, Salma Sylla, Pál Szabó, Róbert Szabó , et al. (3 additional authors not shown)

Abstract: The TESS space telescope is collecting continuous, high-precision optical photometry of stars throughout the sky, including thousands of RR Lyrae stars. In this paper, we present results for an initial sample of 118 nearby RR Lyrae stars observed in TESS Sectors 1 and 2. We use differential-image photometry to generate light curves and analyse their mode content and modulation properties. We combi… ▽ More The TESS space telescope is collecting continuous, high-precision optical photometry of stars throughout the sky, including thousands of RR Lyrae stars. In this paper, we present results for an initial sample of 118 nearby RR Lyrae stars observed in TESS Sectors 1 and 2. We use differential-image photometry to generate light curves and analyse their mode content and modulation properties. We combine accurate light curve parameters from TESS with parallax and color information from the Gaia mission to create a comprehensive classification scheme. We build a clean sample, preserving RR Lyrae stars with unusual light curve shapes, while separating other types of pulsating stars. We find that a large fraction of RR Lyrae stars exhibit various low-amplitude modes, but the distribution of those modes is markedly different from those of the bulge stars. This suggests that differences in physical parameters have an observable effect on the excitation of extra modes, potentially offering a way to uncover the origins of these signals. However, mode identification is hindered by uncertainties when identifying the true pulsation frequencies of the extra modes. We compare mode amplitude ratios in classical double-mode stars to stars with extra modes at low amplitudes and find that they separate into two distinct groups. Finally, we find a high percentage of modulated stars among the fundamental-mode pulsators, but also find that at least 28% of them do not exhibit modulation, confirming that a significant fraction of stars lack the Blazhko effect. △ Less

Submitted 15 September, 2021; originally announced September 2021.

Comments: 35 pages, 26 figures, and a 19-page appendix. Accepted for publication in ApJS. Light curve data temporarily available: https://konkoly.hu/staff/lmolnar/tess_rrl_firstlight_all_lcs.txt

Accuracy of power spectra in dissipationless cosmological simulations

Authors: Sara Maleubre, Daniel Eisenstein, Lehman H. Garrison, Michael Joyce

Abstract: We exploit a suite of large \emph{N}-body simulations (up to N=$4096^3$) performed with \Abacus, of scale-free models with a range of spectral indices $n$, to better understand and quantify convergence of the matter power spectrum. Using self-similarity to identify converged regions, we show that the maximal wavenumber resolved at a given level of accuracy increases monotonically as a function of… ▽ More We exploit a suite of large \emph{N}-body simulations (up to N=$4096^3$) performed with \Abacus, of scale-free models with a range of spectral indices $n$, to better understand and quantify convergence of the matter power spectrum. Using self-similarity to identify converged regions, we show that the maximal wavenumber resolved at a given level of accuracy increases monotonically as a function of time. At the 1\% level it starts at early times from a fraction of $k_Λ$, the Nyquist wavenumber of the initial grid, and reaches at most, if the force softening is sufficiently small, $\sim 2-3 k_Λ$ at the very latest times we evolve to. At the $5\%$ level, accuracy extends up to wavenumbers of order $5k_Λ$ at late times. Expressed as a suitable function of the scale-factor, accuracy shows a very simple $n$-dependence, allowing a extrapolation to place conservative bounds on the accuracy of \emph{N}-body simulations of non-scale free models like LCDM. We note that deviations due to discretization in the converged range are not well modelled by shot noise, and subtracting it in fact degrades accuracy. Quantitatively our findings are broadly in line with the conservative assumptions about resolution adopted by recent studies using large cosmological simulations (e.g. Euclid Flagship) aiming to constrain the mildly non-linear regime. On the other hand, we remark that conclusions about small scale clustering (e.g. concerning the validity of stable clustering) obtained using PS data at wavenumbers larger than a few $k_Λ$ may need revision in light of our convergence analysis. △ Less

Submitted 29 March, 2022; v1 submitted 9 September, 2021; originally announced September 2021.

Comments: 14 pages, 10 figures, version to match published version. Additional subsection in Shot noise and discreteness effects

Journal ref: MNRAS (2022), 512, 1829-1842

The K2 Galactic Archaeology Program Data Release 3: Age-abundance patterns in C1-C8 and C10-C18

Authors: Joel C. Zinn, Dennis Stello, Yvonne Elsworth, Rafael A. García, Thomas Kallinger, Savita Mathur, Benoît Mosser, Marc Hon, Lisa Bugnet, Caitlin Jones, Claudia Reyes, Sanjib Sharma, Ralph Schönrich, Jack T. Warfield, Rodrigo Luger, Andrew Vanderburg, Chiaki Kobayashi, Marc H. Pinsonneault, Jennifer A. Johnson, Daniel Huber, Sven Buder, Meridith Joyce, Joss Bland-Hawthorn, Luca Casagrande, Geraint F. Lewis , et al. (6 additional authors not shown)

Abstract: We present the third and final data release of the K2 Galactic Archaeology Program (K2 GAP) for Campaigns C1-C8 and C10-C18. We provide asteroseismic radius and mass coefficients, $κ_R$ and $κ_M$, for $\sim 19,000$ red giant stars, which translate directly to radius and mass given a temperature. As such, K2 GAP DR3 represents the largest asteroseismic sample in the literature to date. K2 GAP DR3 s… ▽ More We present the third and final data release of the K2 Galactic Archaeology Program (K2 GAP) for Campaigns C1-C8 and C10-C18. We provide asteroseismic radius and mass coefficients, $κ_R$ and $κ_M$, for $\sim 19,000$ red giant stars, which translate directly to radius and mass given a temperature. As such, K2 GAP DR3 represents the largest asteroseismic sample in the literature to date. K2 GAP DR3 stellar parameters are calibrated to be on an absolute parallactic scale based on Gaia DR2, with red giant branch and red clump evolutionary state classifications provided via a machine-learning approach. Combining these stellar parameters with GALAH DR3 spectroscopy, we determine asteroseismic ages with precisions of $\sim 20-30\%$ and compare age-abundance relations to Galactic chemical evolution models among both low- and high-$α$ populations for $α$, light, iron-peak, and neutron-capture elements. We confirm recent indications in the literature of both increased Ba production at late Galactic times, as well as significant contribution to r-process enrichment from prompt sources associated with, e.g., core-collapse supernovae. With an eye toward other Galactic archaeology applications, we characterize K2 GAP DR3 uncertainties and completeness using injection tests, suggesting K2 GAP DR3 is largely unbiased in mass/age and with uncertainties of $2.9\%\,(\rm{stat.})\,\pm0.1\%\,(\rm{syst.})$ & $6.7\%\,(\rm{stat.})\,\pm0.3\%\,(\rm{syst.})$ in $κ_R$ & $κ_M$ for red giant branch stars and $4.7\%\,(\rm{stat.})\,\pm0.3\%\,(\rm{syst.})$ & $11\%\,(\rm{stat.})\,\pm0.9\%\,(\rm{syst.})$ for red clump stars. We also identify percent-level asteroseismic systematics, which are likely related to the time baseline of the underlying data, and which therefore should be considered in TESS asteroseismic analysis. △ Less

Submitted 10 March, 2022; v1 submitted 11 August, 2021; originally announced August 2021.

Comments: Published in ApJ; machine-readable tables available through the open access online journal article

Journal ref: ApJ 926 191 (2022)

Gravitational collapse from cold uniform asymmetric initial conditions

Authors: Francesco Sylos Labini, Michael Joyce

Abstract: Using controlled numerical N-body experiments, we show how, in the collapse dynamics of an initially cold and uniform distribution of particles with a generic asymmetric shape, finite $N$ fluctuations and perturbations induced by the anisotropic gravitational field compete to determine the physical properties of the asymptotic quasi-stationary state. When finite $N$ fluctuations dominate the dynam… ▽ More Using controlled numerical N-body experiments, we show how, in the collapse dynamics of an initially cold and uniform distribution of particles with a generic asymmetric shape, finite $N$ fluctuations and perturbations induced by the anisotropic gravitational field compete to determine the physical properties of the asymptotic quasi-stationary state. When finite $N$ fluctuations dominate the dynamics, the particle energy distribution changes greatly and the final density profile {decays outside its core} as $r^{-4}$ with an $N$-dependent amplitude. On the other hand, in the limit where the anisotropic perturbations dominate, the collapse is softer and the density profile shows a decay as $r^{-3}$, as is typical of halos in cosmological simulations. However, even in this limit, convergence with $N$ of the macroscopic properties of the virialized system, such as the particle energy distributions, the bound mass, and the density profile, is very slow and not clearly established, including for our largest simulations (with $N \sim 10^6$). Our results illustrate the challenges of accurately simulating the first collapsing structures in standard-type cosmological models △ Less

Submitted 4 June, 2021; originally announced June 2021.

Comments: 8 pages, 6 figures. Accepted for publication in Astronomy and Astrophysics

Revised Best Estimates for the Age and Mass of the Methuselah Star HD 140283 using MESA and Interferometry and Implications for 1D Convection

Authors: J. Tang, M. Joyce

Abstract: In light of recently revised observational measurements of the radius and spectroscopic parameters of the extremely old and metal-poor Gaia benchmark star HD 140283 -- also known as the Methuselah star due to prior suggestions that its age is in tension with the age of the Universe -- we present new, best estimates for the star's mass and age from stellar modeling. These are derived using 1D stell… ▽ More In light of recently revised observational measurements of the radius and spectroscopic parameters of the extremely old and metal-poor Gaia benchmark star HD 140283 -- also known as the Methuselah star due to prior suggestions that its age is in tension with the age of the Universe -- we present new, best estimates for the star's mass and age from stellar modeling. These are derived using 1D stellar evolutionary tracks computed with MESA and the most up-to-date measurements from CHARA interferometry. Excluding modeling variance from the uncertainties, we report a mass of $0.809 \pm 0.001 M_{\odot}$ and an age of $12.01 \pm 0.05$ Gyr ($1 σ$). When dominant sources of modeling uncertainty are taken into account, we report $0.81 \pm 0.05 M_{\odot}$ and $12 \pm 0.5$ Gyr, respectively. These results are consistent with recent literature, and the best-fitting age is not in conflict with the currently accepted age of the universe ($13.5$ Gyr; arXiv:1303.5089 [astro-ph.CO]). △ Less

Submitted 20 May, 2021; originally announced May 2021.

Journal ref: Res. Notes AAS 5(2021) 117

Test beam characterization of sensor prototypes for the CMS Barrel MIP Timing Detector

Authors: R. Abbott, A. Abreu, F. Addesa, M. Alhusseini, T. Anderson, Y. Andreev, A. Apresyan, R. Arcidiacono, M. Arenton, E. Auffray, D. Bastos, L. A. T. Bauerdick, R. Bellan, M. Bellato, A. Benaglia, M. Benettoni, R. Bertoni, M. Besancon, S. Bharthuar, A. Bornheim, E. Brücken, J. N. Butler, C. Campagnari, M. Campana, R. Carlin , et al. (174 additional authors not shown)

Abstract: The MIP Timing Detector will provide additional timing capabilities for detection of minimum ionizing particles (MIPs) at CMS during the High Luminosity LHC era, improving event reconstruction and pileup rejection. The central portion of the detector, the Barrel Timing Layer (BTL), will be instrumented with LYSO:Ce crystals and Silicon Photomultipliers (SiPMs) providing a time resolution of about… ▽ More The MIP Timing Detector will provide additional timing capabilities for detection of minimum ionizing particles (MIPs) at CMS during the High Luminosity LHC era, improving event reconstruction and pileup rejection. The central portion of the detector, the Barrel Timing Layer (BTL), will be instrumented with LYSO:Ce crystals and Silicon Photomultipliers (SiPMs) providing a time resolution of about 30 ps at the beginning of operation, and degrading to 50-60 ps at the end of the detector lifetime as a result of radiation damage. In this work, we present the results obtained using a 120 GeV proton beam at the Fermilab Test Beam Facility to measure the time resolution of unirradiated sensors. A proof-of-concept of the sensor layout proposed for the barrel region of the MTD, consisting of elongated crystal bars with dimensions of about 3 x 3 x 57 mm$^3$ and with double-ended SiPM readout, is demonstrated. This design provides a robust time measurement independent of the impact point of the MIP along the crystal bar. We tested LYSO:Ce bars of different thickness (2, 3, 4 mm) with a geometry close to the reference design and coupled to SiPMs manufactured by Hamamatsu and Fondazione Bruno Kessler. The various aspects influencing the timing performance such as the crystal thickness, properties of the SiPMs (e.g. photon detection efficiency), and impact angle of the MIP are studied. A time resolution of about 28 ps is measured for MIPs crossing a 3 mm thick crystal bar, corresponding to an MPV energy deposition of 2.6 MeV, and of 22 ps for the 4.2 MeV MPV energy deposition expected in the BTL, matching the detector performance target for unirradiated devices. △ Less

Submitted 16 July, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Journal ref: Journal of Instrumentation, Volume 16, July 2021

Good and Proper: Self-similarity of N-body Simulations with Proper Force Softening

Authors: Lehman H. Garrison, Michael Joyce, Daniel J. Eisenstein

Abstract: Analysis of self-similarity in scale-free $N$-body simulations reveals the spatial and temporal scales for which statistics measured in cosmological simulations are converged to the physical continuum limit. We examine how the range of scales in which the two-point correlation function is converged depends on the force softening length and whether it is held constant in comoving or proper coordina… ▽ More Analysis of self-similarity in scale-free $N$-body simulations reveals the spatial and temporal scales for which statistics measured in cosmological simulations are converged to the physical continuum limit. We examine how the range of scales in which the two-point correlation function is converged depends on the force softening length and whether it is held constant in comoving or proper coordinates. We find that a proper softening that reaches roughly 1/30th of the inter-particle spacing by the end of the simulation resolves the same spatial and temporal scales as a comoving softening of the same length while using a third fewer time steps, for a range of scale factors typical to $Λ$CDM simulations. We additionally infer an inherent resolution limit, set by the particle mass and scaling as $a^{-1/2}$, beyond which reducing the softening does not improve the resolution. We postulate a mapping of these results with spectral index $n=-2$ to $Λ$CDM simulations. △ Less

Submitted 15 April, 2021; v1 submitted 17 February, 2021; originally announced February 2021.

Comments: 14 pages, 12 figures. Minor changes to match accepted version. Comments welcome!

Characterisation of 92 Southern TESS Candidate Planet Hosts and a New Photometric [Fe/H] Relation for Cool Dwarfs

Authors: Adam D. Rains, Maruša Žerjal, Michael J. Ireland, Thomas Nordlander, Michael S. Bessell, Luca Casagrande, Christopher A. Onken, Meridith Joyce, Jens Kammerer, Harrison Abbot

Abstract: We present the results of a medium resolution optical spectroscopic survey of 92 cool ($3,000 \lesssim T_{\rm eff} \lesssim 4,500\,$K) southern TESS candidate planet hosts, and describe our spectral fitting methodology used to recover stellar parameters. We quantify model deficiencies at predicting optical fluxes, and while our technique works well for $T_{\rm eff}$, further improvements are neede… ▽ More We present the results of a medium resolution optical spectroscopic survey of 92 cool ($3,000 \lesssim T_{\rm eff} \lesssim 4,500\,$K) southern TESS candidate planet hosts, and describe our spectral fitting methodology used to recover stellar parameters. We quantify model deficiencies at predicting optical fluxes, and while our technique works well for $T_{\rm eff}$, further improvements are needed for [Fe/H]. To this end, we developed an updated photometric [Fe/H] calibration for isolated main sequence stars built upon a calibration sample of 69 cool dwarfs in binary systems, precise to $\pm0.19\,$dex, from super-solar to metal poor, over $1.51 < {\rm Gaia}~(B_P-R_P) < 3.3$. Our fitted $T_{\rm eff}$ and $R_\star$ have median precisions of 0.8% and 1.7%, respectively and are consistent with our sample of standard stars. We use these to model the transit light curves and determine exoplanet radii for 100 candidate planets to 3.5% precision and see evidence that the planet-radius gap is also present for cool dwarfs. Our results are consistent with the sample of confirmed TESS planets, with this survey representing one of the largest uniform analyses of cool TESS candidate planet hosts to date. △ Less

Submitted 14 July, 2021; v1 submitted 16 February, 2021; originally announced February 2021.

Comments: Accepted version. 33 Pages, 12 Figures, 9 Tables

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 504, Issue 4, July 2021, pp.5788-5805

The relationship between photometric and spectroscopic oscillation amplitudes from 3D stellar atmosphere simulations

Authors: Yixiao Zhou, Thomas Nordlander, Luca Casagrande, Meridith Joyce, Yaguang Li, Anish M. Amarsi, Henrique Reggiani, Martin Asplund

Abstract: We establish a quantitative relationship between photometric and spectroscopic detections of solar-like oscillations using ab initio, three-dimensional (3D), hydrodynamical numerical simulations of stellar atmospheres. We present a theoretical derivation as proof of concept for our method. We perform realistic spectral line formation calculations to quantify the ratio between luminosity and radial… ▽ More We establish a quantitative relationship between photometric and spectroscopic detections of solar-like oscillations using ab initio, three-dimensional (3D), hydrodynamical numerical simulations of stellar atmospheres. We present a theoretical derivation as proof of concept for our method. We perform realistic spectral line formation calculations to quantify the ratio between luminosity and radial velocity amplitude for two case studies: the Sun and the red giant $ε$ Tau. Luminosity amplitudes are computed based on the bolometric flux predicted by 3D simulations with granulation background modelled the same way as asteroseismic observations. Radial velocity amplitudes are determined from the wavelength shift of synthesized spectral lines with methods closely resembling those used in BiSON and SONG observations. Consequently, the theoretical luminosity to radial velocity amplitude ratios are directly comparable with corresponding observations. For the Sun, we predict theoretical ratios of 21.0 and 23.7 ppm/[m/s] from BiSON and SONG respectively, in good agreement with observations 19.1 and 21.6 ppm/[m/s]. For $ε$ Tau, we predict K2 and SONG ratios of 48.4 ppm/[m/s], again in good agreement with observations 42.2 ppm/[m/s], and much improved over the result from conventional empirical scaling relations which gives 23.2 ppm/[m/s]. This study thus opens the path towards a quantitative understanding of solar-like oscillations, via detailed modelling of 3D stellar atmospheres. △ Less

Submitted 3 February, 2021; originally announced February 2021.

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

The GALAH Survey: Chemical Clocks

Authors: Michael R. Hayden, Sanjib Sharma, Joss Bland-Hawthorn, Lorenzo Spina, Sven Buder, Martin Asplund, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Janez Kos, Geraint F. Lewis, Jane Lin, Karin Lind, Sarah L. Martell, Katharine J. Schlesinger, Jeffrey D. Simpson, Daniel B. Zucker, Tomaz Zwitter, Boquan Chen, Klemen Cotar, Diane Feuillet, Jonti Horner, Meridith Joyce, Thomas Nordlander , et al. (5 additional authors not shown)

Abstract: Previous studies have found that the elemental abundances of a star correlate directly with its age and metallicity. Using this knowledge, we derive ages for a sample of 250,000 stars taken from GALAH DR3 using only their overall metallicity and chemical abundances. Stellar ages are estimated via the machine learning algorithm $XGBoost$, using main sequence turnoff stars with precise ages as our i… ▽ More Previous studies have found that the elemental abundances of a star correlate directly with its age and metallicity. Using this knowledge, we derive ages for a sample of 250,000 stars taken from GALAH DR3 using only their overall metallicity and chemical abundances. Stellar ages are estimated via the machine learning algorithm $XGBoost$, using main sequence turnoff stars with precise ages as our input training set. We find that the stellar ages for the bulk of the GALAH DR3 sample are accurate to 1-2 Gyr using this method. With these ages, we replicate many recent results on the age-kinematic trends of the nearby disk, including the age-velocity dispersion relationship of the solar neighborhood and the larger global velocity dispersion relations of the disk found using $Gaia$ and GALAH. The fact that chemical abundances alone can be used to determine a reliable age for a star have profound implications for the future study of the Galaxy as well as upcoming spectroscopic surveys. These results show that the chemical abundance variation at a given birth radius is quite small, and imply that strong chemical tagging of stars directly to birth clusters may prove difficult with our current elemental abundance precision. Our results highlight the need of spectroscopic surveys to deliver precision abundances for as many nucleosynthetic production sites as possible in order to estimate reliable ages for stars directly from their chemical abundances. Applying the methods outlined in this paper opens a new door into studies of the kinematic structure and evolution of the disk, as ages may potentially be estimated for a large fraction of stars in existing spectroscopic surveys. This would yield a sample of millions of stars with reliable age determinations, and allow precise constraints to be put on various kinematic processes in the disk, such as the efficiency and timescales of radial migration. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: 13 pages, 15 figures, submitted to MNRAS

A precise asteroseismic age and metallicity for HD 139614: a pre-main-sequence star with a protoplanetary disc in Upper-Centaurus Lupus

Authors: Simon J. Murphy, Meridith Joyce, Timothy R. Bedding, Timothy R. White, Mihkel Kama

Abstract: HD\,139614 is known to be a $\sim$14-Myr-old, possibly pre-main-sequence star in the Sco-Cen OB association in the Upper Centaurus-Lupus subgroup, with a slightly warped circumstellar disc containing ring structures hinting at one or more planets. The star's chemical abundance pattern is metal-deficient except for volatile elements, which places it in the $λ$ Boo class and suggests it has recently… ▽ More HD\,139614 is known to be a $\sim$14-Myr-old, possibly pre-main-sequence star in the Sco-Cen OB association in the Upper Centaurus-Lupus subgroup, with a slightly warped circumstellar disc containing ring structures hinting at one or more planets. The star's chemical abundance pattern is metal-deficient except for volatile elements, which places it in the $λ$ Boo class and suggests it has recently accreted gas-rich but dust-poor material. We identify seven dipole and four radial pulsation modes among its $δ$ Sct pulsations using the TESS light curve and an échelle diagram. Precision modelling with the MESA stellar evolution and GYRE stellar oscillation programs confirms it is on the pre-main sequence. Asteroseismic, grid-based modelling suggests an age of $10.75\pm0.77$ Myr, a mass of $1.52\pm0.02$ M$_{\odot}$, and a global metal abundance of $Z=0.0100\pm0.0010$. This represents the first asteroseismic determination of the bulk metallicity of a $λ$ Boo star. The precise age and metallicity offer a benchmark for age estimates in Upper Centaurus--Lupus, and for understanding disc retention and planet formation around intermediate-mass stars. △ Less

Submitted 19 January, 2021; v1 submitted 23 November, 2020; originally announced November 2020.

Comments: Figures make use of colour. Accepted in MNRAS