skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: An Expanded Set of Los Alamos OPLIB Tables in MESA: Type-1 Rosseland-mean Opacities and Solar Models
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 toZ= 10−6toZ= 0.2, which allows improved accuracy of opacities at low and high metallicity, increases the table density in the metallicity rangeZ= 10−4toZ= 0.1 to enhance the accuracy of opacities drawn from interpolations across neighboring metallicities, and adds entries for hydrogen mass fractions betweenX= 0 andX= 0.1 includingX= 10−2, 10−3, 10−4, 10−5, 10−6that can improve stellar models of hydrogen deficient stars. We implement these new OPLIB radiative opacity tables inMESAand 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 ≈20% to 80% across individual chemical mixtures, up to ≈8% and ≈15% at the bottom and top of the solar convection zone respectively, and ≈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.  more » « less
Award ID(s):
1927130
PAR ID:
10513869
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
968
Issue:
2
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 56
Size(s):
Article No. 56
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Observed scatter in the Lyαopacity of quasar sightlines atz< 6 has motivated measurements of the correlation between Lyαopacity and galaxy density, as models that predict this scatter make strong and sometimes opposite predictions for how they should be related. Our previous work associated two highly opaque Lyαtroughs atz∼ 5.7 with a deficit of Lyαemitting galaxies (LAEs). In this work, we survey two of the most highly transmissive lines of sight at this redshift toward thez= 6.02 quasar SDSS J1306+0356 and thez= 6.17 quasar PSO J359-06. We find that both fields are underdense in LAEs within 10h−1Mpc of the quasar sightline, somewhat less extensive than underdensities associated with Lyαtroughs. We combine our observations with three additional fields from the literature and find that while fields with extreme opacities are generally underdense, moderate opacities span a wider density range. The results at high opacities are consistent with models that invoke UV background fluctuations and/or late reionization to explain the observed scatter in intergalactic medium (IGM) Lyαopacities. There is tension at low opacities, however, as the models tend to associate lower IGM Lyαopacities with higher densities. Although the number of fields surveyed is still small, the low-opacity results may support a scenario in which the ionizing background in low-density regions increases more rapidly than some models suggest after becoming ionized. Elevated gas temperatures from recent reionization may also be making these regions more transparent. 
    more » « less
  2. Abstract The cold neutral medium (CNM) is where neutral atomic hydrogen (Hi) is converted into molecular clouds, so the structure and kinematics of the CNM are key drivers of galaxy evolution. Here we provide new constraints on the vertical distribution of the CNM using the recently developedkinematic_scaleheightsoftware package and a large catalog of sensitive Hiabsorption observations. We estimate the thickness of the CNM in the solar neighborhood to beσz∼ 50–90 pc, assuming a Gaussian vertical distribution. This is a factor of ∼2 smaller than typically assumed, indicating that the thickness of the CNM in the solar neighborhood is similar to that found in the inner Galaxy, consistent with recent simulation results. If we consider only structures with Hioptical depthsτ> 0.1 or column densitiesN(Hi) > 1019.5cm−2, which recent work suggests are thresholds for molecule formation, we findσz∼ 50 pc. Meanwhile, for structures withτ< 0.1 or column densitiesN(Hi) < 1019.5cm−2, we findσz∼ 120 pc. These thicknesses are similar to those derived for the thin- and thick-disk molecular cloud populations traced by CO emission, possibly suggesting that cold Hiand CO are well mixed. Approximately 20% of CNM structures are identified as outliers, with kinematics that are not well explained by Galactic rotation. We show that some of these CNM structures—perhaps representing intermediate-velocity clouds—are associated with the Local Bubble wall. We compare our results to recent observations and simulations, and we discuss their implications for the multiphase structure of the Milky Way’s interstellar medium. 
    more » « less
  3. Abstract The dust grain size distribution (GSD) likely varies significantly across star-forming environments in the Universe, but its impact on star formation remains unclear. This ambiguity arises because the GSD interacts nonlinearly with processes like heating, cooling, radiation, and chemistry, which have competing effects and varying environmental dependencies. Processes such as grain coagulation, expected to be efficient in dense star-forming regions, reduce the abundance of small grains and increase that of larger grains. Motivated by this, we investigate the effects of similar GSD variations on the thermochemistry and evolution of giant molecular clouds (GMCs) using magnetohydrodynamic simulations spanning a range of cloud masses and grain sizes, which explicitly incorporate the dynamics of dust grains within the full-physics framework of the STARFORGE project. We find that grain size variations significantly alter GMC thermochemistry: the leading-order effect is that larger grains, under fixed dust mass, GSD dynamic range, and dust-to-gas ratio, result in lower dust opacities. This reduced opacity permits interstellar radiation field and internal radiation photons to penetrate more deeply. This leads to rapid gas heating and inhibited star formation. Star formation efficiency is highly sensitive to grain size, with an order-of-magnitude reduction when grain size dynamic range increases from 10−3–0.1μm to 0.1–10μm. Additionally, warmer gas suppresses low-mass star formation, and decreased opacities result in a greater proportion of gas in diffuse ionized structures. 
    more » « less
  4. Abstract Photochemical hazes are suspected to obscure molecular features, such as water, from detection in the transmission spectra of exoplanets with atmospheric temperatures <800 K. The opacities of laboratory produced organic compounds (tholins) from Khare et al. have become a standard for modeling haze in exoplanet atmospheres. However, these tholins were grown in an oxygen-free, Titan-like environment that is very different from typical assumptions for exoplanets, where C/O ∼ 0.5. This work presents the 0.13–10 μ m complex refractive indices derived from laboratory transmission measurements of tholins grown in environments with different oxygen abundances. With the increasing uptake of oxygen, absorption increases across the entire wavelength range, and a scattering feature around 6 μ m shifts toward shorter wavelengths and becomes more peaked around 5.8 μ m, due to a C = O stretch resonance. Using GJ 1214 b as a test case, we examine the transmission spectra of a sub-Neptune planet with C/O ratios of solar, 1, and 1000 to evaluate the effective differences between our opacities and those of Khare. For an atmosphere with solar hydrogen and helium abundances, we find a difference of 200–1500 ppm, but for high-metallicity ( Z = 1000) environments, the difference may only be 20 ppm. The 1–2 μ m transmission data for GJ 1214 b rule out the Titan-like haze model, and are more consistent with C/O = 1 and C/O = solar haze models. This work demonstrates that using haze opacities that are more consistent with underlying assumptions about bulk atmospheric composition are important for building self-consistent models that appropriately constrain the atmospheric C/O ratio, even when molecular features are obscured. 
    more » « less
  5. Abstract Analysis of inclusions in primitive meteorites reveals that several short-lived radionuclides (SLRs) with half-lives of 0.1–100 Myr existed in the early solar system (ESS). We investigate the ESS origin of107Pd,135Cs, and182Hf, which are produced byslowneutron captures (thes-process) in asymptotic giant branch (AGB) stars. We modeled the Galactic abundances of these SLRs using theOMEGA+galactic chemical evolution (GCE) code and two sets of mass- and metallicity-dependent AGB nucleosynthesis yields (Monash and FRUITY). Depending on the ratio of the mean-lifeτof the SLR to the average length of time between the formations of AGB progenitorsγ, we calculate timescales relevant for the birth of the Sun. Ifτ/γ≳ 2, we predict self-consistent isolation times between 9 and 26 Myr by decaying the GCE predicted107Pd/108Pd,135Cs/133Cs, and182Hf/180Hf ratios to their respective ESS ratios. The predicted107Pd/182Hf ratio indicates that our GCE models are missing 9%–73% of107Pd and108Pd in the ESS. This missing component may have come from AGB stars of higher metallicity than those that contributed to the ESS in our GCE code. Ifτ/γ≲ 0.3, we calculate instead the time (TLE) from the last nucleosynthesis event that added the SLRs into the presolar matter to the formation of the oldest solids in the ESS. For the 2M,Z= 0.01 Monash model we find a self-consistent solution ofTLE= 25.5 Myr. 
    more » « less