We report on our observing campaign of the compact binary merger GW190814, detected by the Advanced LIGO and Advanced Virgo detectors on August 14th, 2019. This signal has the best localisation of any observed gravitational wave (GW) source, with a 90% probability area of 18.5 deg$^2$, and an estimated distance of ~ 240 Mpc. We obtained wide-field observations with the Deca-Degree Optical Transient Imager (DDOTI) covering 88% of the probability area down to a limiting magnitude of $w$ = 19.9 AB. Nearby galaxies within the high probability region were targeted with the Lowell Discovery Telescope (LDT), whereas promising candidate counterparts were characterized through multi-colour photometry with the Reionization and Transients InfraRed (RATIR) and spectroscopy with the Gran Telescopio de Canarias (GTC). We use our optical and near-infrared limits in conjunction with the upper limits obtained by the community to constrain the possible electromagnetic counterparts associated with the merger. A gamma-ray burst seen along its jet's axis is disfavoured by the multi-wavelength dataset, whereas the presence of a burst seen at larger viewing angles is not well constrained. Although our observations are not sensitive to a kilonova similar to AT2017gfo, we can rule out high-mass (> 0.1 M$_{\odot}$) fast-moving (mean velocity >= 0.3c) wind ejecta for a possible kilonova associated with this merger.

The breakout of a supernova blast wave from its progenitor star provides strong constraints on the star and its immediate surroundings. These surroundings are shaped by mass loss from the star and can include a wide variety of inhomogeneities. Here we present results of multi-dimensional radiation-hydrodynamics calculations of the interactions of the supernova blast wave with inhomogeneities in the immediate surroundings of a massive Wolf-Rayet star, calculating the effect these interactions have on the shock breakout signal from supernovae.

The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g. in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide "curtaining", which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent P Cygni features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances.

We explore the impact of pulsar electromagnetic dipole and fallback accretion emission on the luminosity of a suite of kilonova models. The pulsar models are varied over pulsar magnetic field strength, pulsar lifetime, ejecta mass, and elemental abundances; the fallback models are varied over fallback accretion rate and ejecta mass. For the abundances, we use Fe and Nd as representatives of the wind and dynamical ejecta, respectively. We simulate radiative transfer in the ejecta in either 1D spherical or 2D cylindrical spatial geometry. For the grid of 1D simulations, the mass fraction of Nd is 0, $10^{-4}$, or $10^{-3}$ and the rest is Fe. Our models that fit the bolometric luminosity of AT 2017gfo (the kilonova associated with the first neutron star merger discovered in gravitational waves, GW170817) do not simultaneously fit the B, V, and I time evolution. However, we find that the trends of the evolution in B and V magnitudes are better matched by the fallback model relative to the pulsar model, implying the time dependence of the remnant source influences the color evolution. Further exploration of the parameter space and model deficiencies is needed before we can describe AT 2017gfo with a remnant source.

The 2017 detection of the inspiral and merger of two neutron stars in gravitational waves and gamma rays was accompanied by a quickly-reddening transient. Such a transient was predicted to occur following a rapid neutron capture (r-process) nucleosynthesis event, which synthesizes neutron-rich, radioactive nuclei and can take place in both dynamical ejecta and in the wind driven off the accretion torus formed after a neutron star merger. We present the first three-dimensional general relativistic, full transport neutrino radiation magnetohydrodynamics (GRRMHD) simulations of the black hole-accretion disk-wind system produced by the GW170817 merger. We show that the small but non-negligible optical depths lead to neutrino transport globally coupling the disk electron fraction, which we capture by solving the transport equation with a Monte Carlo method. The resulting absorption drives up the electron fraction in a structured, continuous outflow, with electron fraction as high as $Y_e\sim 0.4$ in the extreme polar region. We show via nuclear reaction network and radiative transfer calculations that nucleosynthesis in the disk wind will produce a blue kilonova.

The merger of neutron star binaries is believed to eject a wide range of heavy elements into the universe. By observing the emission from this ejecta, scientists can probe the ejecta properties (mass, velocity and composition distributions). The emission (a.k.a. kilonova) is powered by the radioactive decay of the heavy isotopes produced in the merger and this emission is reprocessed by atomic opacities to optical and infra-red wavelengths. Understanding the ejecta properties requires calculating the dependence of this emission on these opacities. The strong lines in the optical and infra-red in lanthanide opacities have been shown to significantly alter the light-curves and spectra in these wavelength bands, arguing that the emission in these wavelengths can probe the composition of this ejecta. Here we study variations in the kilonova emission by varying individual lanthanide (and the actinide uranium) concentrations in the ejecta. The broad forest of lanthanide lines makes it difficult to determine the exact fraction of individual lanthanides. Nd is an exception. Its opacities above 1 micron are higher than other lanthanides and observations of kilonovae can potentially probe increased abundances of Nd. Similarly, at early times when the ejecta is still hot (first day), the U opacity is strong in the 0.2-1 micron wavelength range and kilonova observations may also be able to constrain these abundances.

The electromagnetic observations of GW170817 were able to dramatically increase our understanding of neutron star mergers beyond what we learned from gravitational waves alone. These observations provided insight on all aspects of the merger from the nature of the gamma-ray burst to the characteristics of the ejected material. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called kilonovae or macronovae. Characteristics of the ejecta include large velocity gradients, relative to supernovae, and the presence of heavy $r$-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound-bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-binned opacities that preserve the integral of the opacity over frequency. Advantages of this area-preserving approach over the traditional expansion-opacity formalism include the ability to pre-calculate opacity tables that are independent of the type of hydrodynamic expansion and that eliminate the computational expense of calculating opacities within radiation-transport simulations. Tabular opacities are generated for all 14 lanthanides as well as a representative actinide element, uranium. We demonstrate that spectral simulations produced with the line-binned opacities agree well with results produced with the more accurate continuous Monte Carlo Sobolev approach, as well as with the commonly used expansion-opacity formalism. Additional investigations illustrate the convergence of opacity with respect to the number of included lines, and elucidate sensitivities to different atomic physics approximations, such as fully and semi-relativistic approaches.

We investigate $^{60}$Fe in massive stars and core-collapse supernovae focussing on uncertainties that influence its production in 15, 20 and 25 $M_\odot$ stars at solar metallicity. We find that the $^{60}$Fe yield is a monotonic increasing function of the uncertain $^{59}$Fe$(n,\gamma)^{60}$Fe cross section and that a factor of 10 reduction in the reaction rate results in a factor 8-10 reduction in the $^{60}$Fe yield; while a factor of 10 increase in the rate increases the yield by a factor 4-7. We find that none of the 189 simulations we have performed are consistent with a core-collapse supernova triggering the formation of the Solar System, and that only models using $^{59}$Fe$(n,\gamma)^{60}$Fe cross section that is less than or equal to that from NON-SMOKER can reproduce the observed $^{60}$Fe/$^{26}$Al line flux ratio in the diffuse ISM. We examine the prospects of detecting old core-collapse supernova remnants (SNRs) in the Milky Way from their $\gamma$-ray emission from the decay of $^{60}$Fe, finding that the next generation of gamma-ray missions could be able to discover up to $\sim100$ such old SNRs as well as measure the $^{60}$Fe yields of a handful of known Galactic SNRs. We also predict the X-ray spectrum that is produced by atomic transitions in $^{60}$Co following its ionization by internal conversion and give theoretical X-ray line fluxes as a function of remnant age as well as the Doppler and fine-structure line broadening effects. The X-ray emission presents an interesting prospect for addressing the missing SNR problem with future X-ray missions.

About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides in this intermediate region. We put forward that the next frontier for AMO science is to explore the AMO complexity that describes most of the Cosmos.

Opacity enhancements for stellar interior conditions have been explored to explain observed pulsation frequencies and to extend the pulsation instability region for B-type main-sequence variable stars. For these stars, the pulsations are driven in the region of the opacity bump of Fe-group elements at $\sim$200,000 K in the stellar envelope. Here we explore effects of opacity enhancements for the somewhat cooler main-sequence A-type stars, in which $p$-mode pulsations are driven instead in the second helium ionization region at $\sim$50,000 K. We compare models using the new LANL OPLIB vs. LLNL OPAL opacities for the AGSS09 solar mixture. For models of 2 solar masses and effective temperature 7600 K, opacity enhancements have only a mild effect on pulsations, shifting mode frequencies and/or slightly changing kinetic-energy growth rates. Increased opacity near the bump at 200,000 K can induce convection that may alter composition gradients created by diffusive settling and radiative levitation. Opacity increases around the hydrogen and 1st He ionization region (13,000 K) can cause additional higher-frequency $p$ modes to be excited, raising the possibility that improved treatment of these layers may result in prediction of new modes that could be tested by observations. New or wider convective zones and higher convective velocities produced by opacity increases could also affect angular momentum transport during evolution. More work needs to be done to quantify the effects of opacity on the boundaries of the pulsation instability regions for A-type stars.

The effects of opacities in driving pulsations and extending the pulsation instability region for B-type stars (SPB and $\beta$ Cep variables) have recently been explored by, e.g., Walczak et al. 2015, 2017, Daszynska-Daszkiewicz et al. 2017; Moravveji 2016, and Cugier 2014. For these stars, the pulsations are driven by the opacity bump of Fe-group elements at around 200,000 K in the stellar envelope. For $\delta$ Scuti pulsations in A-type stars, pulsations are driven instead by the $\kappa$ effect in the second helium ionization region at around 50,000 K. Opacity modifications investigated for B-type stars may also have an effect on the structure of A-type stars and the instability region. Studies of the pulsation instability region for A-type stars show that other factors also influence the pulsations, for example, helium abundance and diffusive settling. Here we explore some of these effects, focusing on modeling of A-type stars. We compare results using the new LANL OPLIB vs. LLNL OPAL opacities for the AGSS09 solar abundance mixture. We comment on prospects for using BRITE observations of A-type stars to constrain opacities.

We report the discovery and monitoring of the near-infrared counterpart (AT2017gfo) of a binary neutron-star merger event detected as a gravitational wave source by Advanced LIGO/Virgo (GW170817) and as a short gamma-ray burst by Fermi/GBM and Integral/SPI-ACS (GRB170817A). The evolution of the transient light is consistent with predictions for the behaviour of a "kilonova/macronova", powered by the radioactive decay of massive neutron-rich nuclides created via r-process nucleosynthesis in the neutron-star ejecta. In particular, evidence for this scenario is found from broad features seen in Hubble Space Telescope infrared spectroscopy, similar to those predicted for lanthanide dominated ejecta, and the much slower evolution in the near-infrared Ks-band compared to the optical. This indicates that the late-time light is dominated by high-opacity lanthanide-rich ejecta, suggesting nucleosynthesis to the 3rd r-process peak (atomic masses A~195). This discovery confirms that neutron-star mergers produce kilo-/macronovae and that they are at least a major - if not the dominant - site of rapid neutron capture nucleosynthesis in the universe.

With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the binary neutron star merger GW170817. The bright, rapidly fading ultraviolet emission indicates a high mass ($\approx0.03$ solar masses) wind-driven outflow with moderate electron fraction ($Y_{e}\approx0.27$). Combined with the X-ray limits, we favor an observer viewing angle of $\approx 30^{\circ}$ away from the orbital rotation axis, which avoids both obscuration from the heaviest elements in the orbital plane and a direct view of any ultra-relativistic, highly collimated ejecta (a gamma-ray burst afterglow).

A long-standing paradigm in astrophysics is that collisions- or mergers- of two neutron stars (NSs) form highly relativistic and collimated outflows (jets) powering gamma-ray bursts (GRBs) of short (< 2 s) duration. However, the observational support for this model is only indirect. A hitherto outstanding prediction is that gravitational wave (GW) events from such mergers should be associated with GRBs, and that a majority of these GRBs should be off-axis, that is, they should point away from the Earth. Here we report the discovery of the X-ray counterpart associated with the GW event GW170817. While the electromagnetic counterpart at optical and infrared frequencies is dominated by the radioactive glow from freshly synthesized r-process material in the merger ejecta, known as kilonova, observations at X-ray and, later, radio frequencies exhibit the behavior of a short GRB viewed off-axis. Our detection of X-ray emission at a location coincident with the kilonova transient provides the missing observational link between short GRBs and GWs from NS mergers, and gives independent confirmation of the collimated nature of the GRB emission.

Modelling the base of the solar convective envelope is a tedious problem. Since the first rotation inversions, solar modellers are confronted with the fact that a region of very limited extent has an enormous physical impact on the Sun. Indeed, it is the transition region from differential to solid body rotation, the tachocline, which furthermore is influenced by turbulence and is also supposed to be the seat of the solar magnetic dynamo. Moreover, solar models show significant disagreement with the sound speed profile in this region. In this paper, we show how helioseismology can provide further constraints on this region by carrying out an inversion of the Ledoux discriminant. We compare these inversions for Standard Solar Models built using various opacity tables and chemical abundances and discuss the origins of the discrepancies between Solar Models and the Sun.

The Sun is the most constrained and well-studied of all stars. As a consequence, the physical ingredients entering solar models are used as a reference to study all other stars observed in the Universe. However, our understanding of the solar structure is still imperfect, as illustrated by the current debate on the heavy element abundances in the Sun. We wish to provide additional information on the solar structure by carrying out structural inversions of a new physical quantity, a proxy of the entropy of the solar plasma which properties are very sensitive to the temperature gradient below the convective zone. We use new structural kernels to carry out direct inversions of an entropy proxy of the solar plasma and compare the solar structure to various standard solar models built using various opacity tables and chemical abundances. We also link our results to classical tests commonly found in the literature. Our analysis allows us to probe more efficiently the uncertain regions of the solar models, just below the convective zone, paving the way for new in-depth analyses of the Sun taking into account additional physical uncertainties of solar models beyond the specific question of chemical abundances.

The electromagnetic transients accompanying compact binary mergers ($\gamma$-ray bursts, afterglows and 'macronovae') are crucial to pinpoint the sky location of gravitational wave sources. Macronovae are caused by the radioactivity from freshly synthesised heavy elements, e.g. from dynamic ejecta and various types of winds. We study macronova signatures by using multi-dimensional radiative transfer calculations. We employ the radiative transfer code SuperNu and state-of-the-art LTE opacities for a few representative elements from the wind and dynamical ejecta (Cr, Pd, Se, Te, Br, Zr, Sm, Ce, Nd, U) to calculate synthetic light curves and spectra for a range of ejecta morphologies. The radioactive power of the resulting macronova is calculated with the detailed input of decay products. We assess the detection prospects for our most complex models, based on the portion of viewing angles that are sufficiently bright, at different cosmological redshifts ($z$). The brighter emission from the wind is unobscured by the lanthanides (or actinides) in some of the models, permitting non-zero detection probabilities for redshifts up to $z=0.07$. We also find the nuclear mass model and the resulting radioactive heating rate are crucial for the detectability. While for the most pessimistic heating rate (from the FRDM model) no reasonable increase in the ejecta mass or velocity, or wind mass or velocity, can possibly make the light curves agree with the observed nIR excess after GRB130603B, a more optimistic heating rate (from the Duflo-Zuker model) leads to good agreement. We conclude that future reliable macronova observations would constrain nuclear heating rates, and consequently help constrain nuclear mass models.

Gravitational wave observations need accompanying electromagnetic signals to accurately determine the sky positions of the sources. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called macronovae. Characteristics of the ejecta include large velocity gradients and the presence of heavy $r$-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound-bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-smeared opacities that preserve the integral of the opacity over frequency, which is an alternative approach to the expansion-opacity formalism. The use of area-preserving line profiles produces frequency-dependent opacities that are one to two orders of magnitude greater than those obtained with the use of expansion opacities in the Sobolev approximation. Opacities are calculated for four $r$-process elements (cerium, neodymium, samarium and uranium) using fully and semi-relativistic methods, as well as different amounts of configuration interaction in the atomic structure, in order to test the sensitivity of the emission to the underlying atomic physics. We determine the effect of these new opacities on simulated spectra and broad-band light curves by applying a multi-dimensional ray-trace method to the ejecta predicted from 3-dimensional merger calculations. With our substantially larger opacities the simulations yield slightly lower luminosities that peak in the mid-IR. These results suggest that those radioactively powered transients that are related to the very heaviest $r$-process material are more difficult to observe than previously believed.

The analysis of the BRITE oscillation spectrum of the main sequence early B-type star $\nu$ Eridani is presented. Only models with the modified mean opacity profile can account for the observed frequency ranges as well as for the values of some individual frequencies. The number of the $\kappa$ modified seismic models is constrained by the nonadiabatic parameter $f$, which is very sensitive to the opacity changes in the subphotospheric layers where the pulsations are driven. We present an example of the model that satisfies all the above conditions. It seems that the OPLIB opacities are preferred over those from the OPAL and OP projects. Moreover, we discuss additional consequences of the opacity modification, namely, an enhancement of the efficiency of convection in the Z-bump as well as an occurrence of close radial modes which is a kind of avoided-crossing phenomenon common for nonradial modes in standard main sequence models.

The hydrogen and helium accreted by X-ray bursting neutron stars is periodically consumed in runaway thermonuclear reactions that cause the entire surface to glow brightly in X-rays for a few seconds. With models of the emission, the mass and radius of the neutron star can be inferred from the observations. By simultaneously probing neutron star masses and radii, X-ray bursts are one of the strongest diagnostics of the nature of matter at extremely high densities. Accurate determinations of these parameters are difficult, however, due to the highly non-ideal nature of the atmospheres where X-ray bursts occur. Observations from X-ray telescopes such as RXTE and NuStar can potentially place strong constraints on nuclear matter once uncertainties in atmosphere models have been reduced. Here we discuss current progress on modeling atmospheres of X-ray bursting neutron stars and some of the challenges still to be overcome.

Recent R-matrix calculations claim to produce a significant enhancement in the opacity of Fe XVII due to atomic core excitations [S. N. Nahar & A.K. Pradhan, Phys. Rev. Letters 116, 235003 (2016), arXiv:1606.02731] and assert that this enhancement is consistent with recent measurements of higher-than-predicted iron opacities [J. E. Bailey et al., Nature 517, 56 (2015)]. This comment shows that the standard opacity models which have already been directly compared with experimental data produce photon absorption cross-sections for Fe XVII that are effectively equivalent to (and in fact larger than) the new R-matrix opacities. Thus, the new R-matrix results cannot be expected to significantly impact the existing discrepancies between theory and experiment because they produce neither a "large enhancement" nor account for "missing continuum plasma opacity" relative to standard models.

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".

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 (EOS) 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.

We have entered the era of explosive transient astronomy, in which upcoming real-time surveys like the Large Synoptic Survey Telescope (LSST), the Palomar Transient Factory (PTF) and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) will detect supernovae in unprecedented numbers. Future telescopes such as the James Webb Space Telescope may discover supernovae from the earliest stars in the universe and reveal their masses. The observational signatures of these astrophysical transients are the key to unveiling their central engines, the environments in which they occur, and to what precision they will pinpoint cosmic acceleration and the nature of dark energy. We present a new method for modeling supernova light curves and spectra with the radiation hydrodynamics code RAGE coupled with detailed monochromatic opacities in the SPECTRUM code. We include a suite of tests that demonstrate how the improved physics is indispensable to modeling shock breakout and light curves.

The merger of two white dwarfs (a.k.a. double degenerate merger) has often been cited as a potential progenitor of type Ia supernovae. Here we combine population synthesis, merger and explosion models with radiation-hydrodynamics light-curve models to study the implications of such a progenitor scenario on the observed type Ia supernova population. Our standard model, assuming double degenerate mergers do produce thermonuclear explosions, produces supernova light-curves that are broader than the observed type Ia sample. In addition, we discuss how the shock breakout and spectral features of these double degenerate progenitors will differ from the canonical bare Chandrasekhar-massed explosion models. We conclude with a discussion of how one might reconcile these differences with current observations.

Astronomers have proposed a number of mechanisms to produce supernova explosions. Although many of these mechanisms are now not considered primary engines behind supernovae, they do produce transients that will be observed by upcoming ground-based surveys and NASA satellites. Here we present the first radiation-hydrodynamics calculations of the spectra and light curves from three of these "failed" supernovae: supernovae with considerable fallback, accretion induced collapse of white dwarfs, and energetic helium flashes (also known as type .Ia supernovae).

New, more detailed calculations of the emission spectra of the He-like K ALPHA complex of calcium, iron and nickel have been carried out using data from both distorted-wave and R-matrix calculations. The value of the GD ratio (an extended definition of the G ratio that accounts for the effect of resolved and unresolved satellite lines) is significantly enhanced at temperatures below the temperature of He-like maximum abundance. Furthermore it is shown that satellite lines are important contributors to the GD ratio such that GD/G>1 at temperatures well above the temperature of maximum abundance. These new calculations demonstrate, with an improved treatment of the KLn (n>=3) satellite lines, that K ALPHA satellite lines need to be included in models of He like spectra even at relatively high temperatures. The excellent agreement between spectra and line ratios calculated from R-matrix and distorted-wave data also confirms the validity of models based on distorted-wave data for highly charged systems, provided the effect of resonances are taken into account as independent processes.

The General Spectral Modeling (GSM) code employs the quasi-static approximation, a standard, low-density methodology that assumes the ionization balance is separable from a determination of the excited-state populations that give rise to the spectra. GSM also allows for some states to be treated only as contributions to effective rates. While these two approximations are known to be valid at low densities, this work investigates using such methods to model high-density, non-LTE emission spectra and determines at what point the approximations break down by comparing to spectra produced by the LANL code ATOMIC which makes no such approximations. As both approximations are used by other astrophysical and low-density modeling codes, the results should be of broad interest. He-like K$\alpha$ emission spectra are presented for Ni, Fe, and Ca, in order to gauge the effect of both approximations employed in GSM. This work confirms that at and above the temperature of maximum abundance of the He-like ionization stage, the range of validity for both approximations is sufficient for modeling the low- and moderate-density regimes one typically finds in astrophysical and magnetically confined fusion plasmas. However, a breakdown does occur for high densities; we obtain quantitative limits that are significantly higher than previous works. This work demonstrates that, while the range of validity for both approximations is sufficient to predict the density-dependent quenching of the z line, the approximations break down at higher densities. Thus these approximations should be used with greater care when modeling high-density plasmas such as those found in inertial confinement fusion and electromagnetic pinch devices.

It is shown that above the temperature of maximum abundance, recombination rates into the excited states of He-like ions that are calculated using earlier, more approximate methods differ markedly from rates obtained from recent distorted-wave and R-Matrix calculations (unified recombination rate coefficients) for Ca, Fe, and Ni. The present rates lead to G ratios that are greatly lower than those resulting from the more approximate rates in previous works, by up to a factor of six at high electron temperatures. Excellent agreement between the distorted-wave and the R-Matrix rates, as well as excellent agreement in the G ratios calculated from them, provides support for the accuracy of these new values which have a broad applicability to the modelling and interpreting of X-ray spectra from a variety of astrophysical and laboratory sources.