Monochromatic Opacities Research Articles

Radiative acceleration calculation methods and abundance anomalies in Am stars

Atomic diffusion with radiative levitation is a major transport process to consider to explain abundance anomalies in Am stars. Radiative accelerations vary from one species to another, yielding different abundance anomalies at the stellar surface. Radiative accelerations can be computed using different methods: some evolution codes use an analytical approximation, while others calculate them from monochromatic opacities. We compared the abundance evolutions predicted using these various methods. Our models were computed with the Toulouse-Geneva evolution code, in which both an analytical approximation (the single-valued parameter method) and detailed calculations from Opacity Project (OP) atomic data are implemented for the calculation of radiative accelerations. The time evolutions of the surface abundances were computed using macroscopic transport processes that are able to reproduce observed Am star surface abundances in presence of atomic diffusion, namely an ad hoc turbulent model or a global mass loss. The radiative accelerations obtained with the various methods are globally in agreement for all the models below the helium convective zone, but can be much greater between the bottom of the hydrogen convective zone and that of helium. The time evolutions of the surface abundances mostly agree within the observational error, but the abundance of some elements can exceed this error for the least massive mass-loss model. The gain in computing time from using analytical approximations is significant compared to sequential calculations from monochromatic opacities for the turbulence models and for the least massive wind model; the gain is small otherwise. Test calculations of turbulence models with the tabulated OPAL opacities yield quite similar abundances as OP for most elements but in a much shorter time, meaning that determining Am star parameters can be done using a two-step method.

Overview of the contributions from all lanthanide elements to kilonova opacity in the temperature range from 25 000 to 40 000 K

Context. It is now well established that the neutron star (NS) merger is at the origin of the production of trans-iron heavy elements in the universe. These elements are therefore present in large quantities in the ejected matter, whose electromagnetic radiation, called kilonova, is characterized by a significant opacity due to the high density of spectral lines belonging to many heavy ions. Among these, the lanthanide ions play an essential role since, with their open 4f subshell, they have a considerable number of transitions that can absorb emitted light. The knowledge of the atomic structure and the radiative parameters of these ions as well as the determination of the corresponding opacities is therefore of paramount importance for the spectral analysis of kilonovae. Aims. The main goal of the present work is to determine the relative contributions of the different lanthanide elements to the opacity of the emission spectrum of a kilonova in its early phase, that is, a few hours after the NS merger, where the conditions are such that the temperature is between 25 000 and 40 000 K. At these temperatures, the lanthanide ions whose charge states are between V and VII are predominant. Methods. We used the pseudo-relativistic Hartree–Fock (HFR) method extensively to calculate the relevant atomic data (energy levels, wavelengths, and oscillator strengths) in La-Lu V-VII ions. The corresponding monochromatic opacities were estimated from the expansion formalism. Results. We calculated the spectroscopic parameters for a total of more than 800 million radiative transitions in all the ions considered. These data were used to estimate the expansion opacities and Planck mean opacities for all the lanthanide elements at early-phase kilonova conditions between 25 000 and 40 000 K, making it possible to deduce the respective contributions of each element as a function of temperature. Atomic calculations were also carried out with the fully relativistic Multiconfiguration Dirac-Hartree-Fock (MCDHF) method in the specific case of the Yb V ion, as the available experimental data had not yet been compared with the theoretical calculations in our previous studies on lanthanide ions.

Opacity for realistic 3D MHD simulations of cool stellar atmospheres

Context. Realistic three-dimensional time-dependent simulations of stellar near-surface convection employ the opacity binning method for the efficient and accurate computation of the radiative energy exchange. The method provides several orders of magnitude of speedup, but its implementation includes a number of free parameters. Aims. Our aim is to evaluate the accuracy of the opacity binning method as a function of the choice of these free parameters. Methods. The monochromatic opacities computed with the SYNSPEC code were used to construct opacity distribution function (ODF) that was then verified through detailed comparison with the results of the ATLAS code. The opacity binning method was implemented with the SYNSPEC opacities for four representative cool main-sequence stellar spectral types (F3V, G2V, K0V, and M2V). Results. The ODFs from SYNSPEC and ATLAS show consistent results for the opacity and bolometric radiative energy exchange rate Q in the case of the F-, G-, and K-type stars. Significant differences, coming mainly from the molecular line lists, are found for the M-type star. It is possible to optimise a small number of bins to reduce the deviation of the results coming from the opacity grouping with respect to the ODF for the F-, G-, and K-type stars. In the case of the M-type star, the inclusion of splitting in wavelength is needed in the grouping to get similar results, with a subsequent increase in computing time. In the limit of a large number of bins, the deviation for all the binning configurations tested saturates and the results do not converge to the ODF solution. Due to this saturation, the Q rate cannot be improved by increasing the number of bins to more than about 20 bins. The more effective strategy is to select the optimal location of fewer bins.

Photoionization and Opacity

Opacity determines radiation transport through material media. In a plasma source, the primary contributors to atomic opacity are bound–bound line transitions and bound-free photoionization into the continuum. We review the theoretical methodology for state-of-the-art photoionization calculations based on the R-matrix method as employed in the Opacity Project, the Iron Project, and solution of the heretofore unsolved problem of plasma broadening of autoionizing resonances due to electron impact, Stark (electric microfields), Doppler (thermal), and core-excitations. R-matrix opacity calculations entail huge amount of atomic data and calculations of unprecedented complexity. It is shown that in high-energy-density (HED) plasmas, photoionization cross sections become 3-D energy–temperature–density-dependent owing to considerable attenuation of autoionizing resonance profiles. Hence, differential oscillator strengths and monochromatic opacities are redistributed in energy. Consequently, Rosseland and Planck mean opacities are affected significantly.

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

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.

Quantitative comparison of opacities calculated using the R-matrix and distorted-wave methods: Fe xvii

ABSTRACT We present here a detailed calculation of opacities for Fe xvii at the physical conditions corresponding to the base of the Solar convection zone. Many ingredients are involved in the calculation of opacities. We review the impact of each ingredient on the final monochromatic and mean opacities (Rosseland and Planck). The necessary atomic data were calculated with the R-matrix and the distorted-wave (DW) methods. We study the effect of broadening, of resolution, of the extent of configuration sets and of configuration interaction to understand the differences between several theoretical predictions as well as the existing large disagreement with measurements. New Dirac R-matrix calculations including all configurations up to the n = 4, 5, and 6 complexes have been performed as well as corresponding Breit–Pauli DW calculations. The DW calculations have been extended to include autoionizing initial levels. A quantitative contrast is made between comparable DW and R-matrix models. We have reached self-convergence with n = 6 R-matrix and DW calculations. Populations in autoionizing initial levels contribute significantly to the opacities and should not be neglected. The R-matrix and DW results are consistent under the similar treatment of resonance broadening. The comparison with the experiment shows a persistent difference in the continuum while the filling of the windows shows some improvement. This study defines our path to the next generation of opacities and opacity tables for stellar modelling.

Multiconfiguration Dirac–Hartree–Fock radiative parameters for emission lines in Ce ii–ivions and cerium opacity calculations for kilonovae

ABSTRACTLarge-scale calculations of atomic structures and radiative properties have been carried out for singly, doubly, and trebly ionized cerium. For this purpose, the purely relativistic multiconfiguration Dirac–Hartree–Fock (MCDHF) method was used, taking into account the effects of valence–valence and core–valence electronic correlations in detail. The results obtained were then used to calculate the expansion opacities characterizing the kilonovae observed as a result of neutron star mergers. Comparisons with previously published experimental and theoretical studies have shown that the results presented in this work are the most complete currently available, in terms of quantity and quality, concerning the atomic data and monochromatic opacities for Ce ii, Ce iii, and Ce iv ions.

Analysis of radiative opacities for optically thin and thick astrophysical plasmas

The resolution of the radiative transfer equation in radiation-hydrodynamic simulations of astrophysical plasmas require radiative opacities. In this work, an analysis of the monochromatic and multigroup opacities of an astrophysical plasma mixture has been carried out. The study has been made in ranges of electron temperatures and densities of 1−1000 eV and 1011−1020 cm−3, respectively, a wide range of plasma conditions that can be found in several astrophysical scenarios where local and non-local thermodynamic regimes are attained. Collisional-radiative calculations were performed to obtain the plasma level populations and the monochromatic opacities in that range of plasma conditions, covering both thermodynamic regimes. Since the astrophysical mixture includes chemical elements from hydrogen to iron, their contribution to the total opacity will depend on the plasma conditions and we have made a characterization of their contribution as a function of the electron density and temperature and also of the photon frequency. Multigroup and gray approaches are commonly used in the radiative transfer equation in radiation-hydrodynamic calculations. We have analyzed the influence of the number of the groups in the accuracy of the multigroup opacities and we have showed that, for a given plasma condition, the opacity of the multicomponent plasma in the gray approach may be considerably influenced by only some of the contributing elements, due to the influence of the weighting function in the mean, which can lead to great differences with respect to the monochromatic opacity. Finally, since there are situations in which the self-absorption of the plasma radiation becomes relevant due to the dimensions of the plasma, we have performed an analysis of the influence of the radiation trapping in the monochromatic and multigroup opacities in terms of the plasma conditions and the width of the plasma slab, assuming the plasma with planar geometry and we have also studied the departures of the local thermodynamic regime.

Computation of Atomic Astrophysical Opacities

The revision of the standard Los Alamos opacities in the 1980–1990s by a group from the Lawrence Livermore National Laboratory (OPAL) and the Opacity Project (OP) consortium was an early example of collaborative big-data science, leading to reliable data deliverables (atomic databases, monochromatic opacities, mean opacities, and radiative accelerations) widely used since then to solve a variety of important astrophysical problems. Nowadays the precision of the OPAL and OP opacities, and even of new tables (OPLIB) by Los Alamos, is a recurrent topic in a hot debate involving stringent comparisons between theory, laboratory experiments, and solar and stellar observations in sophisticated research fields: the standard solar model (SSM), helio and asteroseismology, non-LTE 3D hydrodynamic photospheric modeling, nuclear reaction rates, solar neutrino observations, computational atomic physics, and plasma experiments. In this context, an unexpected downward revision of the solar photospheric metal abundances in 2005 spoiled a very precise agreement between the helioseismic indicators (the radius of the convection zone boundary, the sound-speed profile, and helium surface abundance) and SSM benchmarks, which could be somehow reestablished with a substantial opacity increase. Recent laboratory measurements of the iron opacity in physical conditions similar to the boundary of the solar convection zone have indeed predicted significant increases (30–400%), although new systematic improvements and comparisons of the computed tables have not yet been able to reproduce them. We give an overview of this controversy, and within the OP approach, discuss some of the theoretical shortcomings that could be impairing a more complete and accurate opacity accounting.

Analysis of microscopic properties of radiative shock experiments performed at the Orion laser facility

In this work we have conducted a study on the radiative and spectroscopic properties of the radiative precursor and the post-shock region from experiments with radiative shocks in xenon performed at the Orion laser facility. The study is based on post-processing of radiation-hydrodynamics simulations of the experiment. In particular, we have analyzed the thermodynamic regime of the plasma, the charge state distributions, the monochromatic opacities and emissivities, and the specific intensities for plasma conditions of both regions. The study of the intensities is a useful tool to estimate ranges of electron temperatures present in the xenon plasma in these experiments and the analysis performed of the microscopic properties commented above helps to better understand the intensity spectra. Finally, a theoretical analysis of the possibility of the onset of isobaric thermal instabilities in the post-shock has been made, concluding that the instabilities obtained in the radiative-hydrodynamic simulations could be thermal ones due to strong radiative cooling.

Solar neutrino flux at keV energies

We calculate the solar neutrino and antineutrino flux in the keV energy range. The dominant thermal source processes are photo production (γ e→ e νν̄), bremsstrahlung (e+Ze→ Ze+e+νν̄), plasmon decay (γ→νν̄), and νν̄ emission in free-bound and bound-bound transitions of partially ionized elements heavier than hydrogen and helium. These latter processes dominate in the energy range of a few keV and thus carry information about the solar metallicity. To calculate their rate we use libraries of monochromatic photon radiative opacities in analogy to a previous calculation of solar axion emission. Our overall flux spectrum and many details differ significantly from previous works. While this low-energy flux is not measurable with present-day technology, it could become a significant background for future direct searches for keV-mass sterile neutrino dark matter.

Radiative properties of matter based on quantum statistical method

We present the preliminary results of our code OPAQS (opacity calculation using quantum statistical model) that is based on the self consistent Hartree–Fock–Slater model for the average atom. The code is capable of performing robust calculations of average charge state, frequency-dependent and mean opacities. The accuracy of the atomic model is verified by comparing the calculations of average charge state with various published results. The monochromatic opacities for iron computed at different sets of temperatures and densities are compared with LEDCOP. The Rosseland and Planck opacities for iron and aluminum are validated with some state-of-the-art codes. The results are in good agreement with the published data.

Influence of atomic kinetics in the simulation of plasma microscopic properties and thermal instabilities for radiative bow shock experiments.

Numerical simulations of laboratory astrophysics experiments on plasma flows require plasma microscopic properties that are obtained by means of an atomic kinetic model. This fact implies a careful choice of the most suitable model for the experiment under analysis. Otherwise, the calculations could lead to inaccurate results and inappropriate conclusions. First, a study of the validity of the local thermodynamic equilibrium in the calculation of the average ionization, mean radiative properties, and cooling times of argon plasmas in a range of plasma conditions of interest in laboratory astrophysics experiments on radiative shocks is performed in this work. In the second part, we have made an analysis of the influence of the atomic kinetic model used to calculate plasma microscopic properties of experiments carried out on magpie on radiative bow shocks propagating in argon. The models considered were developed assuming both local and nonlocal thermodynamic equilibrium and, for the latter situation, we have considered in the kinetic model different effects such as external radiation field and plasma mixture. The microscopic properties studied were the average ionization, the charge state distributions, the monochromatic opacities and emissivities, the Planck mean opacity, and the radiative power loss. The microscopic study was made as a postprocess of a radiative-hydrodynamic simulation of the experiment. We have also performed a theoretical analysis of the influence of these atomic kinetic models in the criteria for the onset possibility of thermal instabilities due to radiative cooling in those experiments in which small structures were experimentally observed in the bow shock that could be due to this kind of instability.

SOLAR OPACITY CALCULATIONS USING THE SUPER-TRANSITION-ARRAY METHOD

ABSTRACT A new opacity model has been developed based on the Super-Transition-Array (STA) method for the calculation of monochromatic opacities of plasmas in local thermodynamic equilibrium. The atomic code, named STAR (STA-Revised), is described and used to calculate spectral opacities for a solar model implementing the recent AGSS09 composition. Calculations are carried out throughout the solar radiative zone. The relative contributions of different chemical elements and atomic processes to the total Rosseland mean opacity are analyzed in detail. Monochromatic opacities and charge-state distributions are compared with the widely used Opacity Project (OP) code, for several elements near the radiation–convection interface. STAR Rosseland opacities for the solar mixture show a very good agreement with OP and the OPAL opacity code throughout the radiation zone. Finally, an explicit STA calculation was performed of the full AGSS09 photospheric mixture, including all heavy metals. It was shown that, due to their extremely low abundance, and despite being very good photon absorbers, the heavy elements do not affect the Rosseland opacity.

A NEW GENERATION OF LOS ALAMOS OPACITY TABLES

ABSTRACT We present a new, publicly available set of Los Alamos OPLIB opacity tables for the elements hydrogen through zinc. Our tables are computed using the Los Alamos ATOMIC opacity and plasma modeling code, and make use of atomic structure calculations that use fine-structure detail for all the elements considered. Our equation of state model, known as ChemEOS, is based on the minimization of free energy in a chemical picture and appears to be a reasonable and robust approach to determining atomic state populations over a wide range of temperatures and densities. In this paper we discuss in detail the calculations that we have performed for the 30 elements considered, and present some comparisons of our monochromatic opacities with measurements and other opacity codes. We also use our new opacity tables in solar modeling calculations and compare and contrast such modeling with previous work.

Structure of a laser-driven radiative shock

Radiative shocks are ubiquitous in stellar environments and are characterized by high temperature plasma emitting a considerable fraction of their energy as radiation. The physical structure of these shocks is complex and experimental benchmarks are needed to provide a deeper understanding of the physics at play. In addition, experiments provide unique data for testing radiation hydrodynamics codes which, in turn, are used to model astrophysical phenomena.Radiative shocks have been studied on various high-energy laser facilities for more than a decade, highlighting the importance of radiation on the plasma dynamics. Particularly the PALS facility has focused in producing radiative shocks with typical velocities of ∼50–60 km s−1 in xenon at a fraction of a bar. In addition PALS has the unique capability of producing the most powerful XUV laser available today (21.2 nm (58.4 eV), 0.15 ns), opening the door to new diagnostics of dense plasmas. Here we present results of XUV imaging of the precursor and post-shock structure of radiative shocks generated in xenon in this facility, together with time-and-space resolved measurements of the XUV self-emission using fast diode. The experimental results are interpreted with the help of 2D ARWEN radiative hydrodynamics simulations and state-of-the art monochromatic opacities.

Opacities and spectra of hydrogen atmospheres of moderately magnetized neutron stars

There is observational evidence that central compact objects (CCOs) in supernova remnants have moderately strong magnetic fields $B\sim10^{11}$ G. Meanwhile, available models of partially ionized hydrogen atmospheres of neutron stars with strong magnetic fields are restricted to $B\gtrsim10^{12}$ G. We extend the equation of state and radiative opacities, presented in previous papers for $10^{12}\mbox{ G}\lesssim B \lesssim 10^{15}$ G, to weaker fields. An equation of state and radiative opacities for a partially ionized hydrogen plasma are obtained at magnetic fields $B$, temperatures $T$, and densities $\rho$ typical for atmospheres of CCOs and other isolated neutron stars with moderately strong magnetic fields. The first- and second-order thermodynamic functions, monochromatic radiative opacities, and Rosseland mean opacities are calculated and tabulated, taking account of partial ionization, for $3\times10^{10}\mbox{ G}\lesssim B\lesssim 10^{12}$ G, $10^5$ K $\lesssim T\lesssim 10^7$ K, and a wide range of densities. Atmosphere models and spectra are calculated to verify the applicability of the results and to determine the range of magnetic fields and effective temperatures where the incomplete ionization of the hydrogen plasma is important.

UTILITARIAN OPACITY MODEL FOR AGGREGATE PARTICLES IN PROTOPLANETARY NEBULAE AND EXOPLANET ATMOSPHERES

As small solid grains grow into larger ones in protoplanetary nebulae, or in the cloudy atmospheres of exoplanets, they generally form porous aggregates rather than solid spheres. A number of previous studies have used highly sophisticated schemes to calculate opacity models for irregular, porous particles with size much smaller than a wavelength. However, mere growth itself can affect the opacity of the medium in far more significant ways than the detailed compositional and/or structural differences between grain constituents once aggregate particle sizes exceed the relevant wavelengths. This physics is not new; our goal here is to provide a model that provides physical insight and is simple to use in the increasing number of protoplanetary nebula evolution, and exoplanet atmosphere, models appearing in recent years, yet quantitatively captures the main radiative properties of mixtures of particles of arbitrary size, porosity, and composition. The model is a simple combination of effective medium theory with small-particle closed-form expressions, combined with suitably chosen transitions to geometric optics behavior. Calculations of wavelength-dependent emission and Rosseland mean opacity are shown and compared with Mie theory. The model's fidelity is very good in all comparisons we have made, except in cases involving pure metal particles or monochromatic opacities for solid particles with size comparable to the wavelength.

Solar axion flux from the axion-electron coupling

In non-hadronic axion models, where axions couple to electrons at tree level, the solar axion flux is completely dominated by the ABC reactions (Atomic recombination and deexcitation, Bremsstrahlung and Compton). In this paper the ABC flux is computed from available libraries of monochromatic photon radiative opacities (OP, LEDCOP and OPAS) by exploiting the relations between axion and photon emission cross sections. These results turn to be ∼ 30% larger than previous estimates due to atomic recombination (free-bound electron transitions) and deexcitation (bound-bound), which where not previously taken into account.

Ab initio calculation of the non-relativistic free–free Gaunt factor incorporating plasma screening

We present calculations of Gaunt factors for free–free absorption over a wide range of temperatures and densities. The calculations employ a partial wave expansion approach, which is able to account for plasma screening within the calculation of the free–free Gaunt factor. Much of the existing Gaunt factor data pertains to hydrogenic systems, and plasma screening is often incorporated in opacity calculations using approximate methods. The use of a more accurate method allows us to determine the accuracy of such approximations in calculations of the free–free monochromatic and mean opacities.