Application of New Los Alamos OPLIB Opacities in Solar Modeling Using the Mesa Code
The Sun provides the only local laboratory to test our understanding of stellar physics. The production, then, of valid dynamical models of the Sun is critical to how we think about other stars. Despite improvements in modeling methodologies and abundance measurements, solar models continue to reproduce errors in convection zone depth and surface helium abundance compared to helioseismic observations. Using the 1D stellar evolution code MESAstar, this work first evaluates the effect of changes to timestep resolution and mixing length theory on model outputs. A standard test case is then constructed with both GS98 and AGSS09 abundances using OP and OPAL opacities before implementing the new Los Alamos OPLIB opacities. It is shown that, when using GS98 abundances, OPLIB opacities improve agreement with helioseismic inference of the convection zone base radius, R czb , by 1.68-σ and 1.19-σ compared to OP and OPAL, respectively, while worsening agreement with the inferred envelope helium mass fraction, Y surf , by 0.63-σ and 0.98-σ for OP and OPAL.