Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

Currently, HOPS and AIPS are the primary choices for the time-consuming process of (millimeter) Very Long Baseline Interferometry (VLBI) data calibration. However, for a full end-to-end pipeline, they either lack the ability to perform easily scriptable incremental calibration or do not provide full control over the workflow with the ability to manipulate and edit calibration solutions directly. The Common Astronomy Software Application (CASA) offers all these abilities, together with a secure development future and an intuitive Python interface, which is very attractive for young radio astronomers. Inspired by the recent addition of a global fringe-fitter, the capability to convert FITS-IDI files to measurement sets, and amplitude calibration routines based on ANTAB metadata, we have developed the the CASA-based Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD). The pipeline will be able to handle data from multiple arrays: EHT, GMVA, VLBA and the EVN in the first release. Polarization and phase-referencing calibration are supported and a spectral line mode will be added in the future. The large bandwidths of future radio observatories ask for a scalable reduction software. Within CASA, a message passing interface (MPI) implementation is used for parallelization, reducing the total time needed for processing. The most significant gain is obtained for the time-consuming fringe-fitting task where each scan be processed in parallel.

To assess the effectiveness of optical emission as a probe of spatial asymmetry in core-collapse supernovae (CCSNe), we apply the radiative transfer software, \supernu, to a unimodal CCSN model. The \snsph\ radiation-hydrodynamics software was used to simulate an asymmetric explosion of a 16 M$_{\odot}$ ZAMS binary star. The ejecta has 3.36 M$_{\odot}$ with 0.024 M$_{\odot}$ of radioactive $^{56}$Ni, with unipolar asymmetry along the z-axis. For 96 discrete angular views, we find the ratio between maximum and minimum peak total luminosities is $\sim$1.36. The brightest light curves emerge from views orthogonal to the z-axis. Multigroup spectra from UV to IR are obtained. We find a shift in wavelength with viewing angle in a near-IR Ca II emission feature, consistent with Ca being mostly in the unimode. We compare emission from the grey gamma-ray transfer in \supernu\ and from the detailed gamma-ray transfer code \maverick. Relative to the optical light curves, the brightness of the gamma-ray emission is more monotonic with respect to viewing angle. UBVRI broad-band light curves are also calculated. Parallel with the unimode, the U and B bands have excess luminosity at $\gtrsim 10$ days post-explosion, due to $^{56}$Ni on the unimode. We compare our CCSN model with SN 2002ap, which is thought to have a similar ejecta morphology.

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.

With an increasing number of superluminous supernovae (SLSNe) discovered the question of their origin remains open and causes heated debates in the supernova community. Currently, there are three proposed mechanisms for SLSNe: (1) pair-instability supernovae (PISN), (2) magnetar-driven supernovae, and (3) models in which the supernova ejecta interacts with a circumstellar material ejected before the explosion. Based on current observations of SLSNe, the PISN origin has been disfavoured for a number of reasons. Many PISN models provide overly broad light curves and too reddened spectra, because of massive ejecta and a high amount of nickel. In the current study we re-examine PISN properties using progenitor models computed with the GENEC code. We calculate supernova explosions with FLASH and light curve evolution with the radiation hydrodynamics code STELLA. We find that high-mass models (200 and 250 solar masses) at relatively high metallicity (Z=0.001) do not retain hydrogen in the outer layers and produce relatively fast evolving PISNe Type I and might be suitable to explain some SLSNe. We also investigate uncertainties in light curve modelling due to codes, opacities, the nickel-bubble effect and progenitor structure and composition.

Double-degenerate (DD) mergers of carbon-oxygen white dwarfs have recently emerged as a leading candidate for normal Type Ia supernovae (SNe Ia). However, many outstanding questions surround DD mergers, including the characteristics of their light curves and spectra. We have recently identified a spiral instability in the post-merger phase of DD mergers and demonstrated that this instability self-consistently leads to detonation in some cases. We call this the spiral merger SN Ia model. Here, we utilize the SuperNu radiative transfer software to calculate three-dimensional synthetic light curves and spectra of the spiral merger simulation with a system mass of 2.1 $M_\odot$ from Kashyap et al. Because of their large system masses, both violent and spiral merger light curves are slowly declining. The spiral merger resembles very slowly declining SNe Ia, including SN 2001ay, and provides a more natural explanation for its observed properties than other SN Ia explosion models. Previous synthetic light curves and spectra of violent DD mergers demonstrate a strong dependence on viewing angle, which is in conflict with observations. Here, we demonstrate that the light curves and spectra of the spiral merger are less sensitive to the viewing angle than violent mergers, in closer agreement with observation. We find that the spatial distribution of 56Ni and IMEs follows a characteristic hourglass shape. We discuss the implications of the asymmetric distribution of 56Ni for the early-time gamma-ray observations of 56Ni from SN 2014J. We suggest that DD mergers that agree with the light curves and spectra of normal SNe Ia will likely require a lower system mass.

In Type Ia Supernovae (\sneia), the relative abundances of chemical elements are affected by the neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures in the spectra near maximum light, spectral features may be used to constrain the composition of the progenitor. We calculate the nucleosynthetic yields for three \snia simulations, assuming single degenerate, Chandrasekhar mass progenitors, for a wide range of progenitor metallicities, and calculate synthetic light curves and spectra to explore correlations between progenitor metallicity and the strength of spectral features. We use two 2D simulations of the deflagration-detonation-transition scenario with different $^{56}$Ni yields and the W7 simulation to control for differences between explosion models and total yields. While the overall yields of intermediate mass elements (16 $<$ A $\leq$ 40) differ between the three cases, trends in the yields are similar. With increasing metallicity, $^{28}$Si yields remain nearly constant, $^{40}$Ca yields decline, and Ti and $^{54}$Fe yields increase. In the synthetic spectra, we identify two features at 30 days post explosion that appear to deepen with progenitor metallicity: a Ti feature around 4200\,\AA\ and a Fe feature around 5200\,\AA\@. In all three simulations, their pseudo equivalent widths show a systematic trend with progenitor metallicity. This suggests that these two features may allow differentiation among progenitor metallicities of observed \sneia and potentially help reduce the intrinsic Hubble scatter.

The recurrent symbiotic nova V745 Sco exploded on 2014 February 6 and was observed on February 22 and 23 by the Chandra X-ray Observatory Transmission Grating Spectrometers. By that time the supersoft source phase had already ended and Chandra spectra are consistent with emission from a hot, shock-heated circumstellar medium with temperatures exceeding 10^7K. X-ray line profiles are more sharply peaked than expected for a spherically-symmetric blast wave, with a full width at zero intensity of approximately 2400 km/s, a full width at half maximum of 1200 +/- 30 km/s and an average net blueshift of 165 +/- 10 km/s. The red wings of lines are increasingly absorbed toward longer wavelengths by material within the remnant. We conclude that the blast wave was sculpted by an aspherical circumstellar medium in which an equatorial density enhancement plays a role, as in earlier symbiotic nova explosions. Expansion of the dominant X-ray emitting material is aligned close to the plane of the sky and most consistent with an orbit seen close to face-on. Comparison of an analytical blast wave model with the X-ray spectra, Swift observations and near-infrared line widths indicates the explosion energy was approximately 10^43 erg, and confirms an ejected mass of approximately 10^-7 Msun. The total mass lost is an order of magnitude lower than the accreted mass required to have initiated the explosion, indicating the white dwarf is gaining mass and is a supernova Type 1a progenitor candidate.

Pair Instability Supernovae have been suggested as candidates for some Super Luminous Supernovae, such as SN 2007bi, and as one of the dominant types of explosion occurring in the early Universe from massive, zero-metallicity Population III stars. The progenitors of such events can be rapidly rotating, therefore exhibiting different evolutionary properties due to the effects of rotationally-induced mixing and mass-loss. Proper identification of such events requires rigorous radiation hydrodynamics and radiative transfer calculations that capture not only the behavior of the light curve but also the spectral evolution of these events. We present radiation hydrodynamics and radiation transport calculations for 90-300 Msun rotating pair-instability supernovae covering both the shock break-out and late light curve phases. We also investigate cases of different initial metallicity and rotation rate to determine the impact of these parameters on the detailed spectral characteristics of these events. In agreement with recent results on non-rotating pair instability supernovae, we find that for a range of progenitor masses and rotation rates these events have intrinsically red colors in contradiction with observations of super-luminous supernovae. The spectroscopic properties of rotating pair instability supernovae are similar to those of non-rotating events with stripped hydrogen and helium envelopes. We find that the progenitor metallicity and rotation rate properties are erased after the explosion and cannot be identified in the resulting model spectra. It is the combined effects of pre-supernova mass-loss and the basic properties of the supernova ejecta such as mass, temperature and velocity that have the most direct impact in the model spectra of pair instability supernovae.

We present the results of an intensive multiwavelength campaign on nova LMC 2012. This nova evolved very rapidly in all observed wavelengths. The time to fall two magnitudes in the V band was only 2 days. In X-rays the super soft phase began 13$\pm$5 days after discovery and ended around day 50 after discovery. During the super soft phase, the \Swift/XRT and \Chandra\ spectra were consistent with the underlying white dwarf being very hot, $\sim$ 1 MK, and luminous, $\sim$ 10$^{38}$ erg s$^{-1}$. The UV, optical, and near-IR photometry showed a periodic variation after the initial and rapid fading had ended. Timing analysis revealed a consistent 19.24$\pm$0.03 hr period in all UV, optical, and near-IR bands with amplitudes of $\sim$ 0.3 magnitudes which we associate with the orbital period of the central binary. No periods were detected in the corresponding X-ray data sets. A moderately high inclination system, $i$ = 60$\pm$10$^{\arcdeg}$, was inferred from the early optical emission lines. The {\it HST}/STIS UV spectra were highly unusual with only the \ion{N}{5} (1240\AA) line present and superposed on a blue continuum. The lack of emission lines and the observed UV and optical continua from four epochs can be fit with a low mass ejection event, $\sim$ 10$^{-6}$ M$_{\odot}$, from a hot and massive white dwarf near the Chandrasekhar limit. The white dwarf, in turn, significantly illuminated its subgiant companion which provided the bulk of the observed UV/optical continuum emission at the later dates. The inferred extreme white dwarf characteristics and low mass ejection event favor nova LMC 2012 being a recurrent nova of the U Sco subclass.

Spectra of Type Ia Supernovae (SN Ia) show both continuum-like and absorption-like line features. The pseudo equivalent width (pEW) and Doppler shift of absorption-line-like features, such as Si II 6150A, Ca II HK 3750A, and the Ca II IR-triplet 8150A quantities that are often associated with the optical depth and the velocity of a shell of interacting material, are important tools for interpreting and classifying the SN Ia spectra. In this paper, we examine the nature of spectral features in SN Ia spectra using W7 model spectra and a technique that we call "knock-out" spectra. We show that the P-Cyg profiles of Si II 6150A and many other features are largely emission dominated, rather than absorption dominated, and that the concepts of "absorption line" and "continuum" are therefore not adequate to describe the nature of spectral features in SNe Ia. Apparent absorption features, like Si II 6150A, are frequently just coincidental troughs between two (or more) uncorrelated emission features. In this situation, the pEW measured between these emission peaks is little related to the true strength of the presumed absorption feature. Furthermore, using the same knock-out technique, we demonstrate how spectral features overlap each other at different times after the explosion. This overlap distorts individual line profiles and affects measured absorption-line velocities. With synthetic SN Ia ejecta stratifications that are tuned to match specific observations, the method presented in this paper can in principle be used to quantify the effects of specific line opacities on measured line velocity, line strength, and line blending, and improve the interpretation and informative value of observed SN Ia spectra.

Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) are methods used to stochastically solve the radiative transport and diffusion equations, respectively. These methods combine into a hybrid transport-diffusion method we refer to as IMC-DDMC. We explore a multigroup IMC-DDMC scheme that, in DDMC, combines frequency groups with sufficient optical thickness. We term this procedure "opacity regrouping". Opacity regrouping has previously been applied to IMC-DDMC calculations for problems in which the dependence of the opacity on frequency is monotonic. We generalize opacity regrouping to non-contiguous groups and implement this in \supernu, a code designed to do radiation transport in high-velocity outflows with non-monotonic opacities. We find that regrouping of non-contiguous opacity groups generally improves the speed of IMC-DDMC radiation transport. We present an asymptotic analysis that informs the nature of the Doppler shift in DDMC groups and summarize the derivation of the Gentile-Fleck factor for modified IMC-DDMC. We test \supernu\ using numerical experiments including a quasi-manufactured analytic solution, a simple ten-group problem, and the W7 problem for Type Ia supernovae. We find that the opacity regrouping is necessary to make our IMC-DDMC implementation feasible for the W7 problem and possibly Type Ia supernova simulations in general. We compare the bolometric light curves and spectra produced by the \supernu\ and \phoenix\ radiation transport codes for the W7 problem. The overall shape of the bolometric light curves are in good agreement, as are the spectra and their evolution with time. However, for the numerical specifications we considered, we find that the peak luminosity of the light curve calculated using \supernu\ is $\sim$10% less than that calculated using \phoenix.

We present a systematic study of the pure deflagration model of Type Ia supernovae using three-dimensional, high-resolution, full-star hydrodynamical simulations, nucleosynthetic yields calculated using Lagrangian tracer particles, and light curves calculated using radiation transport. We evaluate the simulations by comparing their predicted light curves with many observed SNe Ia using the SALT2 data-driven model and find that the simulations may correspond to under-luminous SNe Iax. We explore the effects of the initial conditions on our results by varying the number of randomly selected ignition points from 63 to 3500, and the radius of the centered sphere they are confined in from 128 to 384 km. We find that the rate of nuclear burning depends on the number of ignition points at early times, the density of ignition points at intermediate times, and the radius of the confining sphere at late times. The results depend primarily on the number of ignition points, but we do not expect this to be the case in general. The simulations with few ignition points release more nuclear energy $E_{\mathrm{nuc}}$, have larger kinetic energies $E_{\mathrm{K}}$, and produce more $^{56}$Ni than those with many ignition points, and differ in the distribution of $^{56}$Ni, Si, and C/O in the ejecta. For these reasons, the simulations with few ignition points exhibit higher peak B-band absolute magnitudes $M_\mathrm{B}$ and light curves that rise and decline more quickly; their $M_\mathrm{B}$ and light curves resemble those of under-luminous SNe Iax, while those for simulations with many ignition points are not.

We explore Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) for radiation transport in high-velocity outflows with structured opacity. The IMC method is a stochastic computational technique for nonlinear radiation transport. IMC is partially implicit in time and may suffer in efficiency when tracking Monte Carlo particles through optically thick materials. DDMC accelerates IMC in diffusive domains. Abdikamalov extended IMC and DDMC to multigroup, velocity-dependent transport with the intent of modeling neutrino dynamics in core-collapse supernovae. Densmore has also formulated a multifrequency extension to the originally grey DDMC method. We rigorously formulate IMC and DDMC over a high-velocity Lagrangian grid for possible application to photon transport in the post-explosion phase of Type Ia supernovae. This formulation includes an analysis that yields an additional factor in the standard IMC-to-DDMC spatial interface condition. To our knowledge the new boundary condition is distinct from others presented in prior DDMC literature. The method is suitable for a variety of opacity distributions and may be applied to semi-relativistic radiation transport in simple fluids and geometries. Additionally, we test the code, called SuperNu, using an analytic solution having static material, as well as with a manufactured solution for moving material with structured opacities. Finally, we demonstrate with a simple source and 10 group logarithmic wavelength grid that IMC-DDMC performs better than pure IMC in terms of accuracy and speed when there are large disparities between the magnitudes of opacities in adjacent groups. We also present and test our implementation of the new boundary condition.

We propose a robust, quantitative method to compare the synthetic light curves of a Type Ia Supernova (SNIa) explosion model with a large set of observed SNeIa, and derive a figure of merit for the explosion model's agreement with observations. The synthetic light curves are fit with the data-driven model SALT2 which returns values for stretch, color, and magnitude at peak brightness, as well as a goodness-of-fit parameter. Each fit is performed multiple times with different choices of filter bands and epoch range in order to quantify the systematic uncertainty on the fitted parameters. We use a parametric population model for the distribution of observed SNIa parameters from large surveys, and extend it to represent red, dim, and bright outliers found in a low-redshift SNIa data set. We discuss the potential uncertainties of this population model and find it to be reliable given the current uncertainties on cosmological parameters. Using our population model, we assign each set of fitted parameters a likelihood of being observed in nature, and a figure of merit based on this likelihood. We define a second figure of merit based on the quality of the light curve fit, and combine the two measures into an overall figure of merit for each explosion model. We compute figures of merit for a variety of 1D, 2D and 3D explosion models and show that our evaluation method allows meaningful inferences across a wide range of light curve quality and fitted parameters.

Type Ia supernovae (SNe Ia) originate from the thermonuclear explosions of carbon-oxygen (C-O) white dwarfs (WDs). The single-degenerate scenario is a well-explored model of SNe Ia where unstable thermonuclear burning initiates in an accreting, Chandrasekhar-mass WD and forms an advancing flame. By several proposed physical processes the rising, burning material triggers a detonation, which subsequently consumes and unbinds the WD. However, if a detonation is not triggered and the deflagration is too weak to unbind the star, a completely different scenario unfolds. We explore the failure of the Gravitationally-Confined Detonation (GCD) mechanism of SNe Ia, and demonstrate through 2D and 3D simulations the properties of failed-detonation SNe. We show that failed-detonation SNe expel a few 0.1 solar masses of burned and partially-burned material and that a fraction of the material falls back onto the WD, polluting the remnant WD with intermediate-mass and iron-group elements, that likely segregate to the core forming an WD whose core is iron rich. The remaining material is asymmetrically ejected at velocities comparable to the escape velocity from the WD, and in response, the WD is kicked to velocities of a few hundred km/s. These kicks may unbind the binary and eject a runaway/hyper-velocity WD. Although the energy and ejected mass of the failed-detonation SN are a fraction of typical thermonuclear SNe, they are likely to appear as sub-luminous low-velocity SNe Ia. Such failed detonations might therefore explain or are related to the observed branch of peculiar SNe Ia, such as the family of low-velocity sub-luminous SNe (SN 2002cx/SN 2008ha-like SNe).

X-ray grating spectra provide the confirmation of continued mass loss from novae in the super-soft source (SSS) phase of the outburst. In this work expanding nova atmosphere models are developed and used to study the effect of mass loss on the SSS spectra. The very high temperatures combined with high expansion velocities and large radial extension make nova in the SSS phase very interesting but also difficult objects to model. The radiation transport code PHOENIX was applied to SSS novae before, but careful analysis of the old results has revealed a number of problems which lead to new methods and improvements to the code: 1) an improved NLTE module (a new opacity formalism, rate matrix solver, global iteration scheme, and temperature correction method); 2) a new hybrid hydrostatic-dynamic nova atmosphere setup; 3) the models are treated in pure NLTE (no LTE approximation for any opacity). With the new framework a modest amount of models (limited by computation time) are calculated. These show: 1) systematic behaviour for various atmospheric conditions, 2) the effect of expansion on the model spectrum is significant, and 3) the spectra are sensitive to the details of the atmospheric structure. The models are compared to the ten well-exposed grating spectra presently available: 5x V4743 Sgr, 3x RS Oph, and 2x V2491 Cyg. Although the models are on a coarse grid they do match the observations surprisingly well. Also, hydrostatic models are computed. The reproduction of the data is clearly inferior to the expanding models and, more importantly, their interpretation with hydrostatic models leads to conclusions opposite to those from expanding models. The models enable the derivation of accurate constraints on the physical conditions deep in the nova atmosphere that are revealed only in the SSS phase.

X-ray grating observations have revealed great detail in the spectra of Novae in the Super Soft Source (SSS) phase. Notable features in the SSS spectra are blue-shifted absorption lines, P-Cygni line profiles, and the absence of strong ionization edges, all of which are indicators of an expanding atmosphere. We present, and make publicly available, a set of 672 wind-type (WT) synthetic spectra, obtained from the expanding NLTE SSS models introduced in Van Rossum 2010 with the PHOENIX stellar atmosphere code. The set presented in this paper is limited to solar abundances with the aim to focus on the basic model parameters and their effect on the spectra, providing the basis upon which abundance effects can be studied using a much bigger non-solar set in the next paper in this series. We fit the WT spectra to the five grating spectra taken in the SSS phase of nova V4743 Sgr 2003 as an example application of the WT models. Within the limits of solar abundances we demonstrate that the following parameters are constrained by the data (in order of decreasing accuracy): column density N_H, bolometric luminosity L_bol, effective temperature T_eff, white dwarf radius R, wind asymptotic velocity v_inf, and the mass-loss rate M_dot. The models are also sensitive to the assumed white dwarf mass M_wd but the effect on the spectra can largely be compensated by the other model parameters. The WT spectra with solar abundances fit the data better than abundance optimized hydro-static models.

A new method is presented for calculating the time evolution of spherically symmetric Type Ia Supernova in the post-explosion phase, enabling light curves and spectra to be simulated in a physically self-consistent way. The commonly exploited radiative equilibrium, that is in essence a /gas energy balance/ condition, is unsuitable for this purpose for important physical and numerical reasons. Firstly, the RE depends on the heating and cooling rates of the gas by the radiation field, two quantities that almost completely cancel and are very difficult to calculate accurately. Secondly, the internal energy of the gas is only a tiny fraction of the total energy in the system (the vast majority of the energy resides in the radiation field), so that the vast majority of the energy is neglected in solving for the energy balance. The method presented in this paper, based on the /radiation energy balance/, addresses the bulk of the energy, does not depend on the heating/cooling rates, guarantees an accurate run of the bolometric luminosity over time while bringing the gas temperatures into consistence with the radiation field. We have implemented the method in the stellar atmosphere code PHOENIX and applied it to the classical W7 model. The results illustrate the importance of each of the four physical contributions to the energy balance as a function of time. The simulated spectra and light curves for W7 show good resemblance to the observations, which demonstrates what can be done using PHOENIX with the REB method.

Super Soft Source (SSS) spectra are powered by nuclear burning on the surface of a white dwarf. The released energy causes a radiatively-driven wind that leads to a radially extended atmosphere around the white dwarf. Significant blue shifts in photospheric absorption lines are found in the spectra of novae during their SSS phase, being an evidence of continued mass loss in this phase. We present spherically symmetric PHOENIX models that account for the expansion of the ejecta. A comparison to a plane parallel, hydrostatic atmosphere model demonstrates that the mass loss can have a significant impact on the model spectra. The dynamic model yields less pronounced absorption edges, and harder X-ray spectra are the result. Therefore, lower effective temperatures are needed to explain the observed spectra. Although both types of models are yet to be fine-tuned in order to accurately determine best fit parameters, the implications on the chemical abundances are going in opposite directions. With the expanding models the requirement for strong depletion of the crucial elements that cause these edges is now avoidable.