-
An Expanded Set of Los Alamos OPLIB Tables in MESA: Type-1 Rosseland-mean Opacities and Solar Models
Authors:
Ebraheem Farag,
Christopher J. Fontes,
F. X. Timmes,
Earl P. Bellinger,
Joyce A. Guzik,
Evan B. Bauer,
Suzannah R. Wood,
Katie Mussack,
Peter Hakel,
James Colgan,
David P. Kilcrease,
Manolo E. Sherrill,
Tryston C. Raecke,
Morgan T. Chidester
Abstract:
We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range to $Z$\,=\,10…
▽ More
We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range to $Z$\,=\,10$^{-6}$ to $Z$\,=\,0.2 which allows improved accuracy of opacities at low and high metallicity, increases the table density in the metallicity range $Z$\,=\,10$^{-4}$ to $Z$\,=\,0.1 to enhance the accuracy of opacities drawn from interpolations across neighboring metallicities, and adds entries for hydrogen mass fractions between $X$\,=\,0 and $X$\,=\,0.1 including $X$\,=\,$10^{-2}, 10^{-3}, 10^{-4}, 10^{-5}, 10^{-6}$ that can improve stellar models of hydrogen deficient stars. We implement these new OPLIB radiative opacity tables in \MESA, and find that calibrated solar models agree broadly with previously published helioseismic and solar neutrino results. We find differences between using the new 1194 OPLIB opacity tables and the 126 OPAL opacity tables range from $\approx$\,20--80\% across individual chemical mixtures, up to $\approx$\,8\% and $\approx$\,15\% at the bottom and top of the solar convection zone respectively, and $\approx$\,7\% in the solar core. We also find differences between standard solar models using different opacity table sources that are on par with altering the initial abundance mixture. We conclude that this new, open-access set of OPLIB opacity tables does not solve the solar modeling problem, and suggest the investigation of physical mechanisms other than the atomic radiative opacity.
△ Less
Submitted 4 June, 2024;
originally announced June 2024.
-
Photons from neutrinos: the gamma ray echo of a supernova neutrino burst
Authors:
Cecilia Lunardini,
Joshua Loeffler,
Mainak Mukhopadhyay,
Matthew J. Hurley,
Ebraheem Farag,
F. X. Timmes
Abstract:
When a star undergoes core collapse, a vast amount of energy is released in a ~10 s long burst of neutrinos of all species. Inverse beta decay in the star's hydrogen envelope causes an electromagnetic cascade which ultimately results in a flare of gamma rays - an "echo" of the neutrino burst - at the characteristic energy of 0.511 MeV. We study the phenomenology and detectability of this flare. It…
▽ More
When a star undergoes core collapse, a vast amount of energy is released in a ~10 s long burst of neutrinos of all species. Inverse beta decay in the star's hydrogen envelope causes an electromagnetic cascade which ultimately results in a flare of gamma rays - an "echo" of the neutrino burst - at the characteristic energy of 0.511 MeV. We study the phenomenology and detectability of this flare. Its luminosity curve is characterized by a fast, seconds-long, rise and an equally fast decline, with a minute- or hour-long plateau in between. For a near-Earth star (distance D<1 kpc) the echo will be observable at near future gamma ray telescopes with an effective area of 10^3 cm^2 or larger. Its observation will inform us on the envelope size and composition. In conjunction with the direct detection of the neutrino burst, it will also give information on the neutrino emission away from the line of sight and will enable tests of neutrino propagation effects between the stellar surface and Earth.
△ Less
Submitted 9 July, 2024; v1 submitted 20 December, 2023;
originally announced December 2023.
-
Stellar Neutrino Emission Across The Mass-Metallicity Plane
Authors:
Ebraheem Farag,
F. X. Timmes,
Morgan T. Chidester,
Samalka Anandagoda,
Dieter H. Hartmann
Abstract:
We explore neutrino emission from nonrotating, single star models across six initial metallicities and seventy initial masses from the zero-age main sequence to the final fate. Overall, across the mass spectrum, we find metal-poor stellar models tend to have denser, hotter and more massive cores with lower envelope opacities, larger surface luminosities, and larger effective temperatures than thei…
▽ More
We explore neutrino emission from nonrotating, single star models across six initial metallicities and seventy initial masses from the zero-age main sequence to the final fate. Overall, across the mass spectrum, we find metal-poor stellar models tend to have denser, hotter and more massive cores with lower envelope opacities, larger surface luminosities, and larger effective temperatures than their metal-rich counterparts. Across the mass-metallicity plane we identify the sequence (initial CNO $\rightarrow$ $^{14}$N $\rightarrow$ $^{22}$Ne $\rightarrow$ $^{25}$Mg $\rightarrow$ $^{26}$Al $\rightarrow$ $^{26}$Mg $\rightarrow$ $^{30}$P $\rightarrow$ $^{30}$Si) as making primary contributions to the neutrino luminosity at different phases of evolution. For the low-mass models we find neutrino emission from the nitrogen flash and thermal pulse phases of evolution depend strongly on the initial metallicity. For the high-mass models, neutrino emission at He-core ignition and He-shell burning depends strongly on the initial metallicity. Anti-neutrino emission during C, Ne, and O burning shows a strong metallicity dependence with $^{22}$Ne($α$,$n$)$^{25}$Mg providing much of the neutron excess available for inverse-$β$ decays. We integrate the stellar tracks over an initial mass function and time to investigate the neutrino emission from a simple stellar population. We find average neutrino emission from simple stellar populations to be 0.5--1.2 MeV electron neutrinos. Lower metallicity stellar populations produce slightly larger neutrino luminosities and average $β$ decay energies. This study can provide targets for neutrino detectors from individual stars and stellar populations. We provide convenient fitting formulae and open access to the photon and neutrino tracks for more sophisticated population synthesis models.
△ Less
Submitted 19 October, 2023;
originally announced October 2023.
-
Seismic Signatures of the $^{12}$C($α$, $γ$)$^{16}$O Reaction Rate in White Dwarf Models with Overshooting
Authors:
Morgan T. Chidester,
F. X. Timmes,
Ebraheem Farag
Abstract:
We consider the combined effects that overshooting and the $^{12}$C($α$, $γ$)$^{16}$O reaction rate have on variable white dwarf stellar models. We find that carbon-oxygen white dwarf models continue to yield pulsation signatures of the current experimental $^{12}$C($α$, $γ$)$^{16}$O reaction rate probability distribution function when overshooting is included in the evolution. These signatures ho…
▽ More
We consider the combined effects that overshooting and the $^{12}$C($α$, $γ$)$^{16}$O reaction rate have on variable white dwarf stellar models. We find that carbon-oxygen white dwarf models continue to yield pulsation signatures of the current experimental $^{12}$C($α$, $γ$)$^{16}$O reaction rate probability distribution function when overshooting is included in the evolution. These signatures hold because the resonating mantle region, encompassing $\simeq$\,0.2\,\Msun\ in a typical $\simeq$\,0.6\,\Msun\ white dwarf model, still undergoes radiative helium burning during the evolution to a white dwarf. Our specific models show two potential low-order adiabatic g-modes, $g_2$ and $g_6$, that signalize the $^{12}$C($α$, $γ$)$^{16}$O reaction rate probability distribution function. Both g-mode signatures induce average relative period shifts of $ΔP/P = 0.44 \%$ and $ΔP/P = 1.33\%$ for $g_2$ and $g_6$ respectively. We find that $g_6$ is a trapped mode, and the $g_2$ period signature is inversely proportional to the $^{12}$C($α$, $γ$)$^{16}$O reaction rate. The $g_6$ period signature generally separates the slower and faster reaction rates, and has a maximum relative period shift of $ΔP/P = 3.45\%$. We conclude that low-order g-mode periods from carbon-oxygen white dwarfs may still serve as viable probes for the $^{12}$C($α$, $γ$)$^{16}$O reaction rate probability distribution function when overshooting is included in the evolution.
△ Less
Submitted 8 July, 2023;
originally announced July 2023.
-
Resolving The Peak Of The Black Hole Mass Spectrum
Authors:
Ebraheem Farag,
Mathieu Renzo,
Robert Farmer,
Morgan T. Chidester,
F. X. Timmes
Abstract:
Gravitational wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernova (PISN). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections from merging BHs. We use \MESA\ to evolve single, no…
▽ More
Gravitational wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernova (PISN). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections from merging BHs. We use \MESA\ to evolve single, non-rotating, massive helium cores with a metallicity of $Z = 10^{-5}$ until they either collapse to form a BH or explode as a PISN without leaving a compact remnant. We calculate the boundaries of the lower BH mass gap for S-factors in the range S(300 keV) = (77,203) keV b, corresponding to the $\pm 3σ$ uncertainty in our high resolution tabulated $^{12}$C($α$,$γ$)$^{16}$O reaction rate probability distribution function. We extensively test the temporal and mass resolution to resolve the theoretical peak of the BH mass spectrum across the BH mass gap. We explore the convergence with respect to convective mixing and nuclear burning, finding that significant time resolution is needed to achieve convergence. We also test adopting a minimum diffusion coefficient to help lower resolution models reach convergence. We establish a new lower edge of the upper mass gap as M\textsubscript{lower} $\simeq$\,60$^{+32}_{-14}$\,\Msun\ from the $\pm 3σ$ uncertainty in the $^{12}\text{C}(α, γ) ^{16}\text{O}$ rate. We explore the effect of a larger 3-$α$ rate on the lower edge of the upper mass gap, finding M\textsubscript{lower} $\simeq$\,69$^{+34}_{-18}$\,\Msun. We compare our results with BHs reported in the Gravitational-Wave Transient Catalog.
△ Less
Submitted 20 August, 2022;
originally announced August 2022.
-
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.
-
On Trapped Modes In Variable White Dwarfs As Probes Of The $^{12}$C($α, γ$)$^{16}$O Reaction Rate
Authors:
Morgan T. Chidester,
Ebraheem Farag,
F. X. Timmes
Abstract:
We seek signatures of the current experimental $^{12}$C$(α,γ)^{16}$O reaction rate probability distribution function in the pulsation periods of carbon-oxygen white dwarf models. We find that adiabatic g-modes trapped by the interior carbon-rich layer offer potentially useful signatures of this reaction rate probability distribution function. Probing the carbon-rich region is relevant because it f…
▽ More
We seek signatures of the current experimental $^{12}$C$(α,γ)^{16}$O reaction rate probability distribution function in the pulsation periods of carbon-oxygen white dwarf models. We find that adiabatic g-modes trapped by the interior carbon-rich layer offer potentially useful signatures of this reaction rate probability distribution function. Probing the carbon-rich region is relevant because it forms during the evolution of low-mass stars under radiative helium burning conditions, mitigating the impact of convective mixing processes. We make direct quantitative connections between the pulsation periods of the identified trapped g-modes in variable WD models and the current experimental $^{12}$C$(α,γ)^{16}$O reaction rate probability distribution function. We find an average spread in relative period shifts of $ΔP/P \simeq \pm$ 2\% for the identified trapped g-modes over the $\pm$ 3$σ$ uncertainty in the $^{12}$C$(α,γ)^{16}$O reaction rate probability distribution function -- across the effective temperature range of observed DAV and DBV white dwarfs and for different white dwarf masses, helium shell masses, and hydrogen shell masses. The g-mode pulsation periods of observed white dwarfs are typically given to 6-7 significant figures of precision. This suggests that an astrophysical constraint on the $^{12}$C$(α,γ)^{16}$O reaction rate could, in principle, be extractable from the period spectrum of observed variable white dwarfs.
△ Less
Submitted 5 July, 2022;
originally announced July 2022.
-
Observing intermediate-mass black holes and the upper--stellar-mass gap with LIGO and Virgo
Authors:
Ajit Kumar Mehta,
Alessandra Buonanno,
Jonathan Gair,
M. Coleman Miller,
Ebraheem Farag,
R. J. deBoer,
M. Wiescher,
F. X. Timmes
Abstract:
Using ground-based gravitational-wave detectors, we probe the mass function of intermediate-mass black holes (IMBHs) wherein we also include BHs in the upper mass gap $\sim 60-130~M_\odot$. Employing the projected sensitivity of the upcoming LIGO and Virgo fourth observing (O4) run, we perform Bayesian analysis on quasi-circular non-precessing, spinning IMBH binaries (IMBHBs) with total masses…
▽ More
Using ground-based gravitational-wave detectors, we probe the mass function of intermediate-mass black holes (IMBHs) wherein we also include BHs in the upper mass gap $\sim 60-130~M_\odot$. Employing the projected sensitivity of the upcoming LIGO and Virgo fourth observing (O4) run, we perform Bayesian analysis on quasi-circular non-precessing, spinning IMBH binaries (IMBHBs) with total masses $50\mbox{--} 500\, M_\odot$, mass ratios 1.25, 4, and 10, and dimensionless spins up to 0.95, and estimate the precision with which the source-frame parameters can be measured. We find that, at $2σ$, the mass of the heavier component of IMBHBs can be constrained with an uncertainty of $\sim 10-40\%$ at a signal-to-noise ratio of $20$. Focusing on the stellar-mass gap with new tabulations of the $^{12}\text{C}(α, γ)^{16} \text{O}$ reaction rate and its uncertanties, we evolve massive helium core stars using \MESA\, to establish the lower and upper edge of the mass gap as $\simeq$\,59$^{+34}_{-13}$\,$M_{\odot}$ and $\simeq$\,139$^{+30}_{-14}$\,$M_{\odot}$ respectively, where the error bars give the mass range that follows from the $\pm 3σ$ uncertainty in the $^{12}\text{C}(α, γ) ^{16} \text{O}$ nuclear reaction rate. We find that high resolution of the tabulated reaction rate and fine temporal resolution are necessary to resolve the peak of the BH mass spectrum. We then study IMBHBs with components lying in the mass gap and show that the O4 run will be able to robustly identify most such systems. Finally, we re-analyse GW190521 with a state-of-the-art aligned-spin waveform model, finding that the primary mass lies in the mass gap with 90\% credibility.
△ Less
Submitted 19 October, 2021; v1 submitted 13 May, 2021;
originally announced May 2021.
-
On The Impact Of 22Ne On The Pulsation Periods Of Carbon-Oxygen White Dwarfs With Helium Dominated Atmospheres
Authors:
Morgan T. Chidester,
F. X. Timmes,
Josiah Schwab,
Richard H. D. Townsend,
Ebraheem Farag,
Anne Thoul,
C. E. Fields,
Evan B. Bauer,
Michael H. Montgomery
Abstract:
We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M$_\odot$ carbon-oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of $^{22}$Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6$-$7 significant figures of precision. Usually white dwarf models wi…
▽ More
We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M$_\odot$ carbon-oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of $^{22}$Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6$-$7 significant figures of precision. Usually white dwarf models without $^{22}$Ne are fit to the observed periods and other properties. The root-mean-square residuals to the $\simeq$ 150$-$400 s low-order g-mode periods are typically in the range of $σ_{rms}$ $\lesssim$ 0.3 s, for a fit precision of $σ_{rms}/ P$ $\lesssim$ 0.3 %. We find average relative period shifts of $ΔP/P$ $\simeq$ $\pm$ 0.5 % for the low-order dipole and quadrupole g-mode pulsations within the observed effective temperature window, with the range of $ΔP/P$ depending on the specific g-mode, abundance of $^{22}$Ne, effective temperature, and mass of the white dwarf model. This finding suggests a systematic offset may be present in the fitting process of specific white dwarfs when $^{22}$Ne is absent. As part of the fitting processes involves adjusting the composition profiles of a white dwarf model, our study on the impact of $^{22}$Ne can provide new inferences on the derived interior mass fraction profiles. We encourage routinely including $^{22}$Ne mass fraction profiles, informed by stellar evolution models, to future generations of white dwarf model fitting processes.
△ Less
Submitted 20 January, 2021;
originally announced January 2021.
-
On Stellar Evolution In A Neutrino Hertzsprung-Russell Diagram
Authors:
Ebraheem Farag,
F. X. Timmes,
Morgan Taylor,
Kelly M. Patton,
R. Farmer
Abstract:
We explore the evolution of a select grid of solar metallicity stellar models from their pre-main sequence phase to near their final fates in a neutrino Hertzsprung-Russell diagram, where the neutrino luminosity replaces the traditional photon luminosity. Using a calibrated \MESA\ solar model for the solar neutrino luminosity ($L_{ν,\odot}$ = 0.02398 $\cdot$ $L_{γ,\odot}$ = 9.1795 $\times$ 10…
▽ More
We explore the evolution of a select grid of solar metallicity stellar models from their pre-main sequence phase to near their final fates in a neutrino Hertzsprung-Russell diagram, where the neutrino luminosity replaces the traditional photon luminosity. Using a calibrated \MESA\ solar model for the solar neutrino luminosity ($L_{ν,\odot}$ = 0.02398 $\cdot$ $L_{γ,\odot}$ = 9.1795 $\times$ 10$^{31}$ erg s$^{-1}$) as a normalization, we identify $\simeq$ 0.3 MeV electron neutrino emission from helium burning during the helium flash (peak $L_ν / L_{ν,\odot} \simeq$ 10$^4$, flux $Φ_{ν, {\rm He \ flash}} \simeq$ 170 (10 pc/$d$)$^{2}$ cm$^{-2}$ s$^{-1}$ for a star located at a distance of $d$ parsec, timescale $\simeq$ 3 days) and the thermal pulse (peak $L_ν / L_{ν,\odot} \simeq$ 10$^9$, flux $Φ_{ν, {\rm TP}} \simeq$ 1.7$\times$10$^7$ (10 pc/$d$)$^{2}$ cm$^{-2}$ s$^{-1}$, timescale $\simeq$ 0.1 yr) phases of evolution in low mass stars as potential probes for stellar neutrino astronomy. We also delineate the contribution of neutrinos from nuclear reactions and thermal processes to the total neutrino loss along the stellar tracks in a neutrino Hertzsprung-Russell diagram. We find, broadly but with exceptions, that neutrinos from nuclear reactions dominate whenever hydrogen and helium burn, and that neutrinos from thermal processes dominate otherwise.
△ Less
Submitted 12 March, 2020;
originally announced March 2020.
-
Investigating Opacity Modifications and Reaction Rate Uncertainties to Resolve the Cepheid Mass Discrepancy
Authors:
Joyce A. Guzik,
Ebraheem Farag,
Jakub Ostrowski,
Nancy R. Evans,
Hilding Neilson,
Sofia Moschou,
Jeremy J. Drake
Abstract:
Cepheid masses derived from pulsations or binary dynamics are generally lower than those derived from stellar evolution models. Recent efforts have been dedicated to investigating the effects of abundances, mass loss, rotation, convection and overshooting prescriptions for modifying the evolution tracks to reduce or remove this Cepheid mass discrepancy. While these approaches are promising, either…
▽ More
Cepheid masses derived from pulsations or binary dynamics are generally lower than those derived from stellar evolution models. Recent efforts have been dedicated to investigating the effects of abundances, mass loss, rotation, convection and overshooting prescriptions for modifying the evolution tracks to reduce or remove this Cepheid mass discrepancy. While these approaches are promising, either alone or in combination, more work is required to distinguish between possible solutions. Here we investigate nuclear reaction rate and opacity modifications on Cepheid evolution using the MESA code. We discuss the effects of opacity increases at envelope temperatures of 200,000-400,000 K proposed to explain the pulsation properties of hybrid main-sequence beta Cep/Slowly Pulsating B (SPB) variables which will evolve into Cepheids. We make use of the RSP nonlinear radial pulsation modeling capability in MESA to calculate periods and radial velocity amplitudes of Galactic Cepheids V1334 Cyg, Polaris, and delta Cep.
△ Less
Submitted 12 April, 2020; v1 submitted 10 February, 2020;
originally announced February 2020.
-
Sound speed and oscillation frequencies for solar models evolved with Los Alamos ATOMIC opacities
Authors:
Joyce A. Guzik,
C. J. Fontes,
P. Walczak,
S. R. Wood,
K. Mussack,
E. Farag
Abstract:
Los Alamos National Laboratory has calculated a new generation of radiative opacities (OPLIB data using the ATOMIC code) for elements with atomic number Z=1-30 with improved physics input, updated atomic data, and finer temperature grid to replace the Los Alamos LEDCOP opacities released in the year 2000. We calculate the evolution of standard solar models including these new opacities, and compar…
▽ More
Los Alamos National Laboratory has calculated a new generation of radiative opacities (OPLIB data using the ATOMIC code) for elements with atomic number Z=1-30 with improved physics input, updated atomic data, and finer temperature grid to replace the Los Alamos LEDCOP opacities released in the year 2000. We calculate the evolution of standard solar models including these new opacities, and compare with models evolved using the Lawrence Livermore National Laboratory OPAL (Iglesias and Rogers 1996) opacities. We use the solar abundance mixture of Asplund et al. (2009). The new Los Alamos ATOMIC opacities have steeper opacity derivatives than those of OPAL for temperatures and densities of the solar interior radiative zone. We compare the calculated nonadiabatic solar oscillation frequencies and solar interior sound speed to observed frequencies and helioseismic inferences. The calculated sound-speed profiles are similar for models evolved using either the updated Iben evolution code (see \cite{Guzik2010}), or the MESA evolution code (Paxton et al., 2015). The LANL ATOMIC opacities partially mitigate the "solar abundance problem".
△ Less
Submitted 14 May, 2016;
originally announced May 2016.
