We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$\alpha$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

SN 2018hti is a Type I superluminous supernova (SLSN~I) with an absolute $g$-band magnitude of $-22.2$ at maximum brightness, discovered in a metal-poor galaxy at a redshift of 0.0612. We present extensive photometric and spectroscopic observations of this supernova, covering the phases from $\sim -35$ days to more than +340 days from the $r$-band maximum. Combining our $BVgri$-band photometry with {\it Swift} UVOT optical/ultraviolet photometry, we calculated the peak luminosity as $\sim 3.5\times10^{44}$ erg s$^{-1}$. Modeling the observed light curve reveals that the luminosity evolution of SN 2018hti can be produced by an ejecta mass of 5.8 $M_\odot$ and a magnetar with a magnetic field of $B=1.8\times10^{13}$~G having an initial spin period of $P_0=1.8$ ms. Based on such a magnetar-powered scenario and a larger sample, a correlation between the spin of the magnetar and the kinetic energy of the ejecta can be inferred for most SLSNe~I, suggesting a self-consistent scenario. Like for other SLSNe~I, the host galaxy of SN 2018hti is found to be relatively faint ($M_{g} = -17.75$ mag) and of low metallicity ($Z=0.3~Z_\odot$), with a star-formation rate of 0.3 $M_\odot$ yr$^{-1}$. According to simulation results of single-star evolution, SN 2018hti could originate from a massive, metal-poor star with a zero-age main sequence (ZAMS) mass of 25--40 $M_\odot$, or from a less massive rotating star with $M_\mathrm{ZAMS} \approx 16$--25 $M_\odot$. For the case of a binary system, its progenitor could also be a star with $M_\mathrm{ZAMS} \gtrsim 25$ $M_\odot$.

We present the discovery and high-cadence follow-up observations of SN 2018ivc, an unusual Type II supernova that exploded in NGC 1068 (D=10.1 Mpc). The light curve of SN 2018ivc declines piecewise-linearly, changing slope frequently, with four clear slope changes in the first 30 days of evolution. This rapidly changing light curve indicates that interaction between the circumstellar material and ejecta plays a significant role in the evolution. Circumstellar interaction is further supported by a strong X-ray detection. The spectra are rapidly evolving and dominated by hydrogen, helium, and calcium emission lines. We identify a rare high-velocity emission-line feature blueshifted at ~7800 km/s (in Ha, Hb, Pb, Pg, HeI, CaII), which is visible from day 18 until at least day 78 and could be evidence of an asymmetric progenitor or explosion. From the overall similarity between SN 2018ivc and SN 1996al, the \Ha{} equivalent width of its parent HII region, and constraints from pre-explosion archival Hubble Space Telescope images, we find that the progenitor of SN 2018ivc could be as massive as 52 Msun but is more likely <12 Msun. SN 2018ivc demonstrates the importance of the early discovery and rapid follow-up observations of nearby supernovae to study the physics and progenitors of these cosmic explosions.

We present multi-wavelength observations of two gap transients followed by the Carnegie Supernova Project-II and supplemented with data obtained by a number of different programs. Here in the first of two papers, we focus on the intermediate luminosity red transient (ILRT) designated SNhunt120, while in a companion paper we examine the luminous red novae AT 2014ej. Our data set for SNhunt120 consists of an early optical discovery, estimated to be within 3 days after outburst, the subsequent optical and near-infrared broadband followup extending over a $\sim$2 month period, two visual- and two near-infrared wavelength spectra, and Spitzer Space Telescope observations extending from early ($+$28 d) to late ($+$1155 d) phases. SNhunt120 resembles other ILRTs such as NGC 300-2008-OT and SN 2008S, and like these other ILRTs, SNhunt120 exhibits prevalent mid-infrared emission at both early and late phases. From the comparison of SNhunt120 and other ILRTs to electron-capture supernova simulations, we find that the current models underestimate the explosion kinetic energy and thereby produce synthetic light curves that over-estimate the luminosity. Finally, examination of pre-outburst Hubble Space Telescope images yields no progenitor detection.

We present optical and near-infrared broadband photometry and optical spectra of AT 2014ej from the the Carnegie Supernova Project-II. These observations are complemented with data from the CHilean Automatic Supernova sEarch, the Public ESO Spectroscopic Survey of Transient Objects, and from the Backyard Observatory Supernova Search. Observational signatures of AT 2014ej reveal that it is similar to other members of the gap-transient subclass known as luminous red novae (LRNe), including the ubiquitous double hump light curve and spectral properties akin to the LRN SN 2017jfs. A medium-dispersion, visual-wavelength spectrum of AT 2014ej taken the Magellan Clay telescope, exhibits a P Cygni H$\alpha$ feature characterized by a blue velocity at zero intensity of $\approx 110$ km s$^{-1}$ and a P Cygni minimum velocity of $\approx70$ km s$^{-1}$, and which we attribute to emission from a circumstellar wind. Inspection of pre-outbust Hubble Space Telescope images yields no conclusive progenitor detection. In comparison with a sample of LRNe from the literature, AT 2014ej lies among the brighter end of the luminosity distribution. Comparison of the ultra-violet, optical, infrared (UVOIR) light curves of well-observed LRNe to common-envelope evolution models from the literature, indicates the models under predict the luminosity of the comparison sample at all phases and also produce inconsistent time-scales of the secondary peak. Future efforts to model LRNe should expand upon the current parameter space explored and therefore may consider more massive systems and a wider range of dynamical timescales.

We present a new method to photometrically delineate between various sub-types of type Ia supernovae (SNe Ia). Using the color-stretch parameters, $s_{BV}$ or $s_{gr}$, and the time of i-band primary maximum relative to the B-band or g-band maximum it is demonstrated that 2003fg-like, 1991bg-like, and 2002cx-like SNe Ia can readily be identified. In the cases of these extreme SNe Ia, their primary i-band maximum occurs after the time of the B or g band maxima. We suggest that the timing of the i-band maximum can reveal the physical state of the SN Ia explosion as it traces: i) the speed of the recombination front of iron group elements in the ejecta, ii) the temperature evolution and rate of adiabatic cooling in the ejecta and, iii) the presence of interaction with a stellar envelope. This photometric sub-typing can be used in conjunction with other SNe analysis, such as the Branch diagram, to examine the physics and diversity of SNe Ia. The results here can also be used to screen out non-Ia SNe from cosmological samples that do not have complete spectroscopic typing. Finally, as future surveys like LSST create large databases of light curves of many objects this photometric identification can be used to readily identify and study the rates and bulk properties of peculiar SNe Ia.

Supernova LSQ13abf was discovered soon after explosion by the La Silla-QUEST Survey and followed by the CSP II at optical and near-IR wavelengths. Our analysis indicates LSQ13abf was discovered within two days of explosion and its first 10 days of evolution reveal a B-band light curve with an abrupt drop in luminosity. Contemporaneously, the V-band light curve exhibits a rise towards a first peak and the r- and i-band light curves show no early peak. The early light-curve evolution of LSQ13abf is reminiscent of the post explosion cooling phase observed in the Type Ib SN 2008D, and the similarity between the two objects extends over weeks. Spectroscopically, LSQ13abf resembles SN 2008D with P Cygni He I features that strengthen over time. Spectral energy distributions are constructed from broad-band photometry, and by fitting black-body (BB) functions a UVOIR light curve is constructed, and the underlying BB-temperature and BB-radius profiles are estimated. Explosion parameters are estimated by simultaneously fitting an Arnett model to the UVOIR light curve and the velocity evolution derived from spectral features, and a post-shock breakout cooling model to the first two epochs of the bolometric evolution. This combined model suggests an explosion energy of 1.3x10$^{51}$ ergs, a relatively high ejecta mass of 5.94 M$_{\odot}$, a Ni mass of 0.16 M$_{\odot}$, and a progenitor-star radius of 28.0 R$_{\odot}$. The ejecta mass suggests the origins of LSQ13abf lie with a >25 M$_{\odot}$ ZAMS progenitor and its radius is three and nine times larger than values estimated from the same analysis applied to observations of SNe 2008D and 1999ex, respectively. Alternatively, comparison of hydrodynamical simulations of >20-25 M$_{\odot}$ ZAMS progenitors that evolve to pre-SN envelope masses around 10 M$_{\odot}$ and extended (~100 R$_{\odot}$) envelopes also match the observations of LSQ13abf.

Correct estimates of stellar extreme ultraviolet (EUV; 100 - 1170 \AA) flux are important for studying the photochemistry and stability of exoplanet atmospheres, as EUV radiation ionizes hydrogen and contributes to the heating, expansion, and potential escape of a planet's upper atmosphere. Contamination from interstellar hydrogen makes observing EUV emission from M stars particularly difficult, and impossible past 100 pc, and necessitates other means to predict the flux in this wavelength regime. We present EUV -- infrared (100 \AA - 5.5 $\mu$m) synthetic spectra computed with the PHOENIX atmospheric code of three early M dwarf planet hosts: GJ 832 (M1.5 V), GJ 176 (M2.5 V), and GJ 436 (M3.5 V). These one-dimensional semiempirical nonlocal thermodynamic equilibrium models include simple temperature prescriptions for the stellar chromosphere and transition region, from where ultraviolet (UV; 100 - 3008 \AA) fluxes originate. We guide our models with Hubble Space Telescope far- and near-UV spectra and discuss the ability to constrain these models using Galaxy Evolution Explorer UV photometry. Our models closely reproduce the observations and predict the unobservable EUV spectrum at a wavelength resolution of < 0.1 \AA. The temperature profiles that best reproduce the observations for all three stars are described by nearly the same set of parameters, suggesting that early M type stars may have similar thermal structures in their upper atmospheres. With an impending UV observation gap and the scarcity of observed EUV spectra for stars less luminous and more distant than the Sun, upper-atmosphere models such as these are important for providing realistic spectra across short wavelengths and for advancing our understanding of the effects of radiation on planets orbiting M stars.

We present $81$ near-infrared (NIR) spectra of $30$ Type II supernovae (SNe II) from the Carnegie Supernova Project-II (CSP-II), the largest such dataset published to date. We identify a number of NIR features and characterize their evolution over time. The NIR spectroscopic properties of SNe II fall into two distinct groups. This classification is first based on the strength of the He I $\lambda1.083\,\mu$m absorption during the plateau phase; SNe II are either significantly above (spectroscopically strong) or below $50$ angstroms (spectroscopically weak) in pseudo equivalent width. However between the two groups, other properties, such as the timing of CO formation and the presence of Sr II, are also observed. Most surprisingly, the distinct weak and strong NIR spectroscopic classes correspond to SNe II with slow and fast declining light curves, respectively. These two photometric groups match the modern nomenclature of SNe IIP and IIL. Including NIR spectra previously published, 18 out of 19 SNe II follow this slow declining-spectroscopically weak and fast declining-spectroscopically strong correspondence. This is in apparent contradiction to the recent findings in the optical that slow and fast decliners show a continuous distribution of properties. The weak SNe II show a high-velocity component of helium that may be caused by a thermal excitation from a reverse-shock created by the outer ejecta interacting with the red supergiant wind, but the origin of the observed dichotomy is not understood. Further studies are crucial in determining whether the apparent differences in the NIR are due to distinct physical processes or a gap in the current data set.

We present a comprehensive dataset of optical and near-infrared photometry and spectroscopy of type~Ia supernova (SN) 2016hnk, combined with integral field spectroscopy (IFS) of its host galaxy, MCG -01-06-070, and nearby environment. Properties of the SN local environment are characterized by means of single stellar population synthesis applied to IFS observations taken two years after the SN exploded. SN 2016hnk spectra are compared to other 1991bg-like SNe Ia, 2002es-like SNe Ia, and Ca-rich transients. In addition, abundance stratification modelling is used to identify the various spectral features in the early phase spectral sequence and the dataset is also compared to a modified non-LTE model previously produced for the sublumnious SN 1999by. SN 2016hnk is consistent with being a sub-luminous (M$_{\rm B}=-16.7$ mag, s$_{\rm BV}$=0.43$\pm$0.03), highly reddened object. IFS of its host galaxy reveals both a significant amount of dust at the SN location, as well as residual star formation and a high proportion of old stellar populations in the local environment compared to other locations in the galaxy, which favours an old progenitor for SN 2016hnk. Inspection of a nebular spectrum obtained one year after maximum contains two narrow emission lines attributed to the forbidden [Ca II] $\lambda\lambda$7291,7324 doublet with a Doppler shift of 700 km s$^{-1}$. Based on various observational diagnostics, we argue that the progenitor of SN 2016hnk was likely a near Chandrasekhar-mass ($M_{\rm Ch}$) carbon-oxygen white dwarf that produced 0.108 $M_\odot$ of $^{56}$Ni. Our modeling suggests that the narrow [Ca II] features observed in the nebular spectrum are associated with $^{48}$Ca from electron capture during the explosion, which is expected to occur only in white dwarfs that explode near or at the $M_{\rm Ch}$ limit.

Our recent work demonstrates a correlation between the high-velocity blue edge, $v_{edge}$, of the iron-peak Fe/Co/Ni $H$-band emission feature and the optical light curve shape of normal, transitional and sub-luminous type Ia Supernovae (SNe Ia). We explain this correlation in terms of SN Ia physics. $v_{edge}$ corresponds to the sharp transition between the complete and incomplete silicon burning regions in the ejecta. It measures the point in velocity space where the outer $^{56}$Ni mass fraction, $X_{\rm{Ni}}$, falls to the order of 0.03-0.10. For a given $^{56}$Ni mass, $M(^{56}Ni)$, $v_{edge}$ is sensitive to the specific kinetic energy $E_{\rm kin}$($M(^{56}Ni)/M_{WD}$) of the corresponding region. Combining $v_{edge}$ with light curve parameters (i.e., s$_{BV}$, $\Delta m_{15,s}$ in $B$ and $V$) allows us to distinguish between explosion scenarios. The correlation between $v_{edge}$ and light-curve shape is consistent with explosion models near the Chandrasekhar limit. However, the available sub-$M_{Ch}$ WD explosion model based on SN 1999by exhibits velocities which are too large to explain the observations. Finally, the sub-luminous SN 2015bo exhibits signatures of a dynamical merger of two WDs demonstrating diversity among explosion scenarios at the faint end of the SNe Ia population.

We present the $H$-band wavelength region of thirty post-maximum light near-infrared (NIR) spectra of fourteen transitional and sub-luminous type Ia supernovae (SNe Ia), extending from $+$5d to +20d relative to the epoch of $B$-band maximum. We introduce a new observable, the blue-edge velocity, $v_{edge}$, of the prominent Fe/Co/Ni-peak $H$-band emission feature which is quantitatively measured. The $v_{edge}$ parameter is found to slowly decrease over sub-type ranging from around $-$13,000km/s for transitional SNe~Ia, down to $-$5,000km/s for the sub-luminous SNe Ia. Furthermore, inspection of the +10$\pm$3d spectra indicates that $v_{edge}$ is correlated with the color-stretch parameter, s$_{BV}$, and hence with peak luminosity. These results follow the previous findings that brighter SNe Ia tend to have $^{56}$Ni located at higher velocities as compared to sub-luminous objects. As $v_{edge}$ is a model-independent parameter, we propose it can be used in combination with traditional observational diagnostics to provide a new avenue to robustly distinguish between leading SNe Ia explosion models.

A major scientific goal of JWST is to probe the epoch of re-ionization of the Universe at z above 6, and up to 20 and beyond. At these redshifts, galaxies are just beginning to form and the observable objects are early black holes, supernovae, and cosmic infrared background. The JWST has the necessary sensitivity to observe these targets individually, but a public deep and wide science enabling survey in the wavelength range from 2-5 $\mu$m is needed to discover these black holes and supernovae and to cover the area large enough for cosmic infrared background to be reliably studied. This enabling survey will also discover a large number of other transients and enable sciences such as supernova cosmology up to z $\sim$ 5, star formation history at high redshift through supernova explosions, faint stellar objects in the Milky Way, and galaxy evolution up to z approaching 10. The results of this survey will also serve as an invaluable target feeder for the upcoming era of ELT and SKA.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

The high energy radiation environment around M dwarf stars strongly impacts the characteristics of close-in exoplanet atmospheres, but these wavelengths are difficult to observe due to geocoronal and interstellar contamination. On account of these observational restrictions, a stellar atmosphere model may be used to compute the stellar extreme ultraviolet (EUV; 100 - 912 \AA) spectrum. We present a case study of the ultra-cool M8 dwarf star, TRAPPIST-1, which hosts seven transiting short-period terrestrial sized planets whose atmospheres will be probed by the James Webb Space Telescope. We construct semi-empirical non-LTE model spectra of TRAPPIST-1 that span EUV to infrared wavelengths (100 \AA\ - 2.5 $\mu$m) using the atmosphere code PHOENIX. These upper-atmosphere models contain prescriptions for the chromosphere and transition region and include newly added partial frequency redistribution capabilities. In the absence of broadband UV spectral observations, we constrain our models using HST Ly$\alpha$ observations from TRAPPIST-1 and GALEX FUV and NUV photometric detections from a set of old M8 stars ($>$1 Gyr). We find that calibrating the models using both data sets separately yield similar FUV and NUV fluxes, and EUV fluxes that range from (1.32 - 17.4) $\times$ 10$^{-14}$ ergs s$^{-1}$ cm$^{-2}$. The results from these models demonstrate that the EUV emission is very sensitive to the temperature structure in the transition region. Our lower activity models predict EUV fluxes similar to previously published estimates derived from semi-empirical scaling relationships, while the highest activity model predicts EUV fluxes a factor of ten higher. Results from this study support the idea that the TRAPPIST-1 habitable zone planets likely do not have much liquid water on their surfaces due to the elevated levels of high energy radiation emitted by the host star.

We present optical and ultraviolet observations of nearby type Ic supernova SN 2017ein as well as detailed analysis of its progenitor properties from both the early-time observations and the prediscovery Hubble Space Telescope (HST) images. The optical light curves started from within one day to $\sim$275 days after explosion, and optical spectra range from $\sim$2 days to $\sim$90 days after explosion. Compared to other normal SNe Ic like SN 2007gr and SN 2013ge, \mbox{SN 2017ein} seems to have more prominent C{\footnotesize II} absorption and higher expansion velocities in early phases, suggestive of relatively lower ejecta mass. The earliest photometry obtained for \mbox{SN 2017ein} show indications of shock cooling. The best-fit obtained by including a shock cooling component gives an estimate of the envelope mass as $\sim$0.02 M$_{\odot}$ and stellar radius as 8$\pm$4 R$_{\odot}$. Examining the pre-explosion images taken with the HST WFPC2, we find that the SN position coincides with a luminous and blue point-like source, with an extinction-corrected absolute magnitude of M$_V$$\sim$$-$8.2 mag and M$_I$$\sim$$-$7.7 mag.Comparisons of the observations to the theoretical models indicate that the counterpart source was either a single WR star or a binary with whose members had high initial masses, or a young compact star cluster. To further distinguish between different scenarios requires revisiting the site of the progenitor with HST after the SN fades away.

Shifting the focus of Type Ia supernova (SN Ia) cosmology to the near-infrared (NIR) is a promising way to significantly reduce the systematic errors, as the strategy minimizes our reliance on the empirical width-luminosity relation and uncertain dust laws. Observations in the NIR are also crucial for our understanding of the origins and evolution of these events, further improving their cosmological utility. Any future experiments in the rest-frame NIR will require knowledge of the SN Ia NIR spectroscopic diversity, which is currently based on a small sample of observed spectra. Along with the accompanying paper, Phillips et al. (2018), we introduce the Carnegie Supernova Project-II (CSP-II), to follow up nearby SNe Ia in both the optical and the NIR. In particular, this paper focuses on the CSP-II NIR spectroscopy program, describing the survey strategy, instrumental setups, data reduction, sample characteristics, and future analyses on the data set. In collaboration with the Harvard-Smithsonian Center for Astrophysics (CfA) Supernova Group, we obtained 661 NIR spectra of 157 SNe Ia. Within this sample, 451 NIR spectra of 90 SNe Ia have corresponding CSP-II follow-up light curves. Such a sample will allow detailed studies of the NIR spectroscopic properties of SNe Ia, providing a different perspective on the properties of the unburned material, radioactive and stable nickel produced, progenitor magnetic fields, and searches for possible signatures of companion stars.

Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically-confirmed type Ia supernova (SN Ia) observed in the $Kepler$ field. The $Kepler$ data revealed an excess emission in its early light curve, allowing to place interesting constraints on its progenitor system (Dimitriadis et al. 2018, Shappee et al. 2018b). Here, we present extensive optical, ultraviolet, and near-infrared photometry, as well as dense sampling of optical spectra, for this object. SN 2018oh is relatively normal in its photometric evolution, with a rise time of 18.3$\pm$0.3 days and $\Delta$m$_{15}(B)=0.96\pm$0.03 mag, but it seems to have bluer $B - V$ colors. We construct the "uvoir" bolometric light curve having peak luminosity as 1.49$\times$10$^{43}$erg s$^{-1}$, from which we derive a nickel mass as 0.55$\pm$0.04M$_{\odot}$ by fitting radiation diffusion models powered by centrally located $^{56}$Ni. Note that the moment when nickel-powered luminosity starts to emerge is +3.85 days after the first light in the Kepler data, suggesting other origins of the early-time emission, e.g., mixing of $^{56}$Ni to outer layers of the ejecta or interaction between the ejecta and nearby circumstellar material or a non-degenerate companion star. The spectral evolution of SN 2018oh is similar to that of a normal SN Ia, but is characterized by prominent and persistent carbon absorption features. The C II features can be detected from the early phases to about 3 weeks after the maximum light, representing the latest detection of carbon ever recorded in a SN Ia. This indicates that a considerable amount of unburned carbon exists in the ejecta of SN 2018oh and may mix into deeper layers.

We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unusual two-component shape, where the flux rises with a steep linear gradient for the first few days, followed by a quadratic rise as seen for typical SNe Ia. This "flux excess" relative to canonical SN Ia behavior is confirmed in our $i$-band light curve, and furthermore, SN 2018oh is especially blue during the early epochs. The flux excess peaks 2.14$\pm0.04$ days after explosion, has a FWHM of 3.12$\pm0.04$ days, a blackbody temperature of $T=17,500^{+11,500}_{-9,000}$ K, a peak luminosity of $4.3\pm0.2\times10^{37}\,{\rm erg\,s^{-1}}$, and a total integrated energy of $1.27\pm0.01\times10^{43}\,{\rm erg}$. We compare SN 2018oh to several models that may provide additional heating at early times, including collision with a companion and a shallow concentration of radioactive nickel. While all of these models generally reproduce the early K2 light curve shape, we slightly favor a companion interaction, at a distance of $\sim$$2\times10^{12}\,{\rm cm}$ based on our early color measurements, although the exact distance depends on the uncertain viewing angle. Additional confirmation of a companion interaction in future modeling and observations of SN 2018oh would provide strong support for a single-degenerate progenitor system.

On 2018 Feb. 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B_{max}=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and pre-discovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System (ATLAS). The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN~Ia light curve, and it unambiguously shows a ~4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a non-degenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing, published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe~Ia light curves. Finally, we use Swift X-ray non-detections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant density CSM these non-detections constrain rho<4.5 * 10^5 cm^-3 at a radius of 4 *10^15 cm from the progenitor star. Assuming a wind-like environment, we place mass-loss limits of Mdot< 8 * 10^-6 M_sun yr^-1 for v_w=100 km s^-1, ruling out some symbiotic progenitor systems.

The Carnegie Supernova Project-II (CSP-II) was an NSF-funded, four-year program to obtain optical and near-infrared observations of a "Cosmology" sample of $\sim100$ Type Ia supernovae located in the smooth Hubble flow ($0.03 \lesssim z \lesssim 0.10$). Light curves were also obtained of a "Physics" sample composed of 90 nearby Type Ia supernovae at $z \leq 0.04$ selected for near-infrared spectroscopic time-series observations. The primary emphasis of the CSP-II is to use the combination of optical and near-infrared photometry to achieve a distance precision of better than 5%. In this paper, details of the supernova sample, the observational strategy, and the characteristics of the photometric data are provided. In a companion paper, the near-infrared spectroscopy component of the project is presented.

The type Ia supernova (SN) 2012fr displayed an unusual combination of its Si II {\lambda}{\lambda}5972, 6355 features. This includes the ratio of their pseudo equivalent widths, placing it at the border of the Shallow Silicon (SS) and Core Normal (CN) spectral subtype in the Branch diagram, while the Si II {\lambda}6355 expansion velocities places it as a High-Velocity (HV) object in the Wang et al. spectral type that most interestingly evolves slowly, placing it in the Low Velocity Gradient (LVG) typing of Benetti et al. Only 5% of SNe Ia are HV and located in the SS+CN portion of the Branch diagram and less than 10% of SNe Ia are both HV and LVG. These features point towards SN 2012fr being quite unusual, similar in many ways to the peculiar SN 2000cx. We modeled the spectral evolution of SN 2012fr to see if we could gain some insight into its evolutionary behavior. We use the parameterized radiative transfer code SYNOW to probe the abundance stratification of SN 2012fr at pre-maximum, maximum, and post-maximum light epochs. We also use a grid of W7 models in the radiative transfer code PHOENIX to probe the effect of different density structures on the formation of the Si II {\lambda}6355 absorption feature at post-maximum epochs. We find that the unusual features observed in SN 2012fr are likely due to a shell-like density enhancement in the outer ejecta. We comment on possible reasons for atypical Ca II absorption features, and suggest that they are related to the Si II features.

We present an early-phase $g$-band light curve and visual-wavelength spectra of the normal Type Ia supernova (SN) 2013gy. The light curve is constructed by determining the appropriate S-corrections to transform KAIT natural-system $B$- and $V$-band photometry and Carnegie Supernova Project natural-system $g$-band photometry to the Pan-STARRS1 $g$-band natural photometric system. A Markov Chain Monte Carlo calculation provides a best-fit single power-law function to the first ten epochs of photometry described by an exponent of $2.16^{+0.06}_{-0.06}$ and a time of first light of MJD 56629.4$^{+0.1}_{-0.1}$, which is $1.93^{+0.12}_{-0.13}$ days (i.e., $<48$~hr) before the discovery date (2013 December 4.84 UT) and $-19.10^{+0.12}_{-0.13}$ days before the time of $B$-band maximum (MJD 56648.5$\pm0.1$). The estimate of the time of first light is consistent with the explosion time inferred from the evolution of the Si II $\lambda$6355 Doppler velocity. Furthermore, discovery photometry and previous nondetection limits enable us to constrain the companion radius down to $R_c \leq 4\,R_{\odot}$. In addition to our early-time constraints, we use a deep +235 day nebular-phase spectrum from Magellan/IMACS to place a stripped H-mass limit of $< 0.018\,M_{\odot}$. Combined, these limits effectively rule out H-rich nondegenerate companions.

We examine the early phase intrinsic $(B-V)_{0}$ color evolution of a dozen Type~Ia supernovae discovered within three days of the inferred time of first light ($t_{first}$) and have $(B-V)_0$ color information beginning within 5 days of $t_{first}$. The sample indicates there are two distinct early populations. The first is a population exhibiting blue colors that slowlybevolve, and the second population exhibits red colors and evolves more rapidly. We find that the early-blue events are all 1991T/1999aa-like with more luminous slower declining light curves than those exhibiting early-red colors. Placing the first sample on the Branch diagram (i.e., ratio of \ion{Si}{2} $\lambda\lambda$5972, 6355 pseudo-Equivalent widths) indicates all blue objects are of the Branch Shallow Silicon (SS) spectral type, while all early-red events except for the 2000cx-like SN~2012fr are of the Branch Core-Normal (CN) or CooL (CL) type. A number of potential processes contributing to the early emission are explored, and we find that, in general, the viewing-angle dependance inherent in the companion collision model is inconsistent with all SS objects with early-time observations being blue and exhibiting an excess. We caution that great care must be taken when interpreting early-phase light curves as there may be a variety of physical processes that are possibly at play and significant theoretical work remains to be done.

We present detailed ultraviolet, optical and near-infrared light curves of the Type Ia supernova (SN) 2012fr, which exploded in the Fornax cluster member NGC 1365. These precise high-cadence light curves provide a dense coverage of the flux evolution from $-$12 to $+$140 days with respect to the epoch of $B$-band maximum (\tmax). Supplementary imaging at the earliest epochs reveals an initial slow, nearly linear rise in luminosity with a duration of $\sim$2.5 days, followed by a faster rising phase that is well reproduced by an explosion model with a moderate amount of $^{56}$Ni mixing in the ejecta. From an analysis of the light curves, we conclude: $(i)$ explosion occurred $< 22$ hours before the first detection of the supernova, $(ii)$ the rise time to peak bolometric ($\lambda > 1800 $\AA) luminosity was $16.5 \pm 0.6$ days, $(iii)$ the supernova suffered little or no host-galaxy dust reddening, $(iv)$ the peak luminosity in both the optical and near-infrared was consistent with the bright end of normal Type Ia diversity, and $(v)$ $0.60 \pm 0.15 M_{\odot}$ of $^{56}$Ni was synthesized in the explosion. Despite its normal luminosity, SN 2012fr displayed unusually prevalent high-velocity \ion{Ca}{2} and \ion{Si}{2} absorption features, and a nearly constant photospheric velocity of the \ion{Si}{2} $\lambda$6355 line at $\sim$12,000 \kms\ beginning $\sim$5 days before \tmax. Other peculiarities in the early phase photometry and the spectral evolution are highlighted. SN 2012fr also adds to a growing number of Type Ia supernovae hosted by galaxies with direct Cepheid distance measurements.

JWST was conceived and built to answer one of the most fundamental questions that humans can address empirically: "How did the Universe make its first stars?". Our First Lights At REionization (FLARE) project transforms the quest for the epoch of reionization from the static to the time domain. It targets the complementary question: "What happened to those first stars?". It will be answered by observations of the most luminous events: supernovae and accretion on to black holes formed by direct collapse from the primordial gas clouds. These transients provide direct constraints on star-formation rates (SFRs) and the truly initial Initial Mass Function (IMF), and they may identify possible stellar seeds of supermassive black holes (SMBHs). Furthermore, our knowledge of the physics of these events at ultra-low metallicity will be much expanded. JWST's unique capabilities will detect these most luminous and earliest cosmic messengers easily in fairly shallow observations. However, these events are very rare at the dawn of cosmic structure formation and so require large area coverage. Time domain astronomy can be advanced to an unprecedented depth by means of a shallow field of JWST reaching 27 mag AB in 2 and 4.4 microns over a field as large as 0.1 square degree visited multiple times each year. Such a survey may set strong constraints or detect massive Pop III SNe at redshifts beyond 10, pinpointing the redshift of the first stars, or at least their death. Based on our current knowledge of superluminous supernovae (SLSNe), such a survey will find one or more SLSNe at redshifts above 6 in five years and possibly several direct collapse black holes. Although JWST is not designed as a wide field survey telescope, we show that such a wide field survey is possible with JWST and is critical in addressing several of its key scientific goals.

We present an analysis of ultraviolet (UV) to near-infrared observations of the fast-declining Type Ia supernovae (SNe Ia) 2007on and 2011iv, hosted by the Fornax cluster member NGC 1404. The B-band light curves of SN 2007on and SN 2011iv are characterised by dm_15(B) decline-rate values of 1.96 mag and 1.77 mag, respectively. Although they have similar decline rates, their peak B- and H-band magnitudes differ by ~0.60 mag and ~0.35 mag, respectively. After correcting for the luminosity vs. decline rate and the luminosity vs. colour relations, the peak B-band and H-band light curves provide distances that differ by ~14% and ~9%, respectively. These findings serve as a cautionary tale for the use of transitional SNe Ia located in early-type hosts in the quest to measure cosmological parameters. Interestingly, even though SN 2011iv is brighter and bluer at early times, by three weeks past maximum and extending over several months, its B-V colour is 0.12 mag redder than that of SN 2007on. To reconcile this unusual behaviour, we turn to guidance from a suite of spherical one-dimensional Chandrasekhar-mass delayed-detonation explosion models. In this context, 56Ni production depends on both the so-called transition density and the central density of the progenitor white dwarf. To first order, the transition density drives the luminosity-width relation, while the central density is an important second-order parameter. Within this context, the differences in the B-V color evolution along the Lira regime suggests the progenitor of SN~2011iv had a higher central density than SN~2007on.

Stripped-envelope (SE) supernovae (SNe) include H-poor (Type IIb), H-free (Type Ib) and He-free (Type Ic) events thought to be associated with the deaths of massive stars. The exact nature of their progenitors is a matter of debate. Here we present the analysis of the light curves of 34 SE SNe published by the Carnegie Supernova Project (CSP-I), which are unparalleled in terms of photometric accuracy and wavelength range. Light-curve parameters are estimated through the fits of an analytical function and trends are searched for among the resulting fit parameters. We found a tentative correlation between the peak absolute $B$-band magnitude and $\Delta m_{15}(B)$, as well as a correlation between the late-time linear slope and $\Delta m_{15}$. Making use of the full set of optical and near-IR photometry, combined with robust host-galaxy extinction corrections, bolometric light curves are constructed and compared to both analytic and hydrodynamical models. From the hydrodynamical models we obtained ejecta masses of $1.1-6.2$ $M_{\odot}$, $^{56}$Ni masses of $0.03-0.35$ $M_{\odot}$, and explosion energies (excluding two SNe Ic-BL) of $0.25-3.0\times10^{51}$ erg. Our analysis indicates that adopting $\kappa = 0.07$ cm$^{2}$ g$^{-1}$ as the mean opacity serves to be a suitable assumption when comparing Arnett-model results to those obtained from hydrodynamical calculations. We also find that adopting He I and O I line velocities to infer the expansion velocity in He-rich and He-poor SNe, respectively, provides ejecta masses relatively similar to those obtained by using the Fe II line velocities. The inferred ejecta masses are compatible with intermediate mass ($M_{ZAMS} \leq 20$ $M_{\odot}$) progenitor stars in binary systems for the majority of SE SNe. Furthermore, the majority of our SNe is affected by significant mixing of $^{56}$Ni, particularly in the case of SNe Ic.

Recently Graphics Processing Units (GPUs) have been used to speed up very CPU-intensive gravitational microlensing simulations. In this work, we use the Xeon Phi coprocessor to accelerate such simulations and compare its performance on a microlensing code with that of NVIDIA's GPUs. For the selected set of parameters evaluated in our experiment, we find that the speedup by Intel's Knights Corner coprocessor is comparable to that by NVIDIA's Fermi family of GPUs with compute capability 2.0, but less significant than GPUs with higher compute capabilities such as the Kepler. However, the very recently released second generation Xeon Phi, Knights Landing, is about 5.8 times faster than the Knights Corner, and about 2.9 times faster than the Kepler GPU used in our simulations. We conclude that the Xeon Phi is a very promising alternative to GPUs for modern high performance microlensing simulations.

We examine the bolometric lightcurves of five Type II-P supernovae (SNe 1998A, 2000cb, 2006V, 2006au and 2009E) which are thought to originate from blue supergiant progenitors using a new python package named SuperBoL. With this code, we calculate SNe lightcurves using three different techniques common in the literature: the quasi-bolometric method, which integrates the observed photometry, the direct integration method, which additionally corrects for unobserved flux in the UV and IR, and the bolometric correction method, which uses correlations between observed colors and V-band bolometric corrections. We present here the lightcurves calculated by SuperBoL along with previously published lightcurves, as well as peak luminosities and Ni-56 yields. We find that the direct integration and bolometric correction lightcurves largely agree with previously published lightcurves, but with what we believe to be more robust error calculations, with $0.2 \leq \delta L_{bol}/L_{bol} \leq 0.5$. Peak luminosities and Ni-56 masses are similarly comparable to previous work. SN 2000cb remains an unusual member of this sub-group, owing to the faster rise and flatter plateau than the other supernovae in the sample. Initial comparisons with the NLTE atmosphere code PHOENIX show that the direct integration technique re-produces the luminosity of a model supernova spectrum to ~5% when given synthetic photometry of the spectrum as input. Our code is publicly available. The ability to produce bolometric lightcurves from observed sets of broad-band light curves should be helpful in the interpretation of other types of supernovae, particularly those that are not well characterized, such as extremely luminous supernovae and faint fast objects.

We present optical spectra of the nearby Type Ia supernova SN 2011fe at 100, 205, 311, 349, and 578 days post-maximum light, as well as an ultraviolet spectrum obtained with Hubble Space Telescope at 360 days post-maximum light. We compare these observations with synthetic spectra produced with the radiative transfer code PHOENIX. The day +100 spectrum can be well fit with models which neglect collisional and radiative data for forbidden lines. Curiously, including this data and recomputing the fit yields a quite similar spectrum, but with different combinations of lines forming some of the stronger features. At day +205 and later epochs, forbidden lines dominate much of the optical spectrum formation; however, our results indicate that recombination, not collisional excitation, is the most influential physical process driving spectrum formation at these late times. Consequently, our synthetic optical and UV spectra at all epochs presented here are formed almost exclusively through recombination-driven fluorescence. Furthermore, our models suggest that the ultraviolet spectrum even as late as day +360 is optically thick and consists of permitted lines from several iron-peak species. These results indicate that the transition to the "nebular" phase in Type Ia supernovae is complex and highly wavelength-dependent.

We analyze a time series of optical spectra of SN~2014J from almost two weeks prior to maximum to nearly four months after maximum. We perform our analysis using the SYNOW code, which is well suited to track the distribution of the ions with velocity in the ejecta. We show that almost all of the spectral features during the entire epoch can be identified with permitted transitions of the common ions found in normal SNe Ia in agreement with previous studies. We show that 2014J is a relatively normal SN Ia. At early times the spectral features are dominated by Si II, S II, Mg II, and Ca II. These ions persist to maximum light with the appearance of Na I and Mg I. At later times iron-group elements also appear, as expected in the stratified abundance model of the formation of normal type Ia SNe. We do not find significant spectroscopic evidence for oxygen, until 100 days after maximum light. The +100 day identification of oxygen is tentative, and would imply significant mixing of unburned or only slight processed elements down to a velocity of 6,000 km/s. Our results are in relatively good agreement with other analyses in the IR. We briefly compare SN 2011fe to SN 2014J and conclude that the differences could be due to different central densities at ignition or differences in the C/O ratio of the progenitors.

We present the intensive spectroscopic follow up of the type Ia supernova (SN Ia) 2014J in the starburst galaxy M82. Twenty-seven optical spectra have been acquired from January 22nd to September 1st 2014 with the Isaac Newton (INT) and William Herschel (WHT) Telescopes. After correcting the observations for the recession velocity of M82 and for Milky Way and host galaxy extinction, we measured expansion velocities from spectral line blueshifts and pseudo-equivalent width of the strongest features in the spectra, which gives an idea on how elements are distributed within the ejecta. We position SN 2014J in the Benetti (2005), Branch et al. (2006) and Wang et al. (2009) diagrams. These diagrams are based on properties of the Si II features and provide dynamical and chemical information about the SN ejecta. The nearby SN 2011fe, which showed little evidence for reddening in its host galaxy, is shown as a reference for comparisons. SN 2014J is a border-line object between the Core-normal (CN) and Broad-line (BL) groups, which corresponds to an intermediate position between Low Velocity Gradient (LVG) and High Velocity Gradient (HVG) objects. SN 2014J follows the R(Si II)-\Delta m15 correlation, which confirms its classification as a relatively normal SN Ia. Our description of the SN Ia in terms of the evolution of the pseudo-equivalent width of various ions as well as the position in the various diagrams put this specific SN Ia into the overall sample of SN Ia.

We use observed UV through near IR spectra to examine whether SN 2011fe can be understood in the framework of Branch-normal SNe Ia and to examine its individual peculiarities. As a benchmark, we use a delayed-detonation model with a progenitor metallicity of Z_solar/20. We study the sensitivity of features to variations in progenitor metallicity, the outer density profile, and the distribution of radioactive nickel. The effect of metallicity variations in the progenitor have a relatively small effect on the synthetic spectra. We also find that the abundance stratification of SN 2011fe resembles closely that of a delayed detonation model with a transition density that has been fit to other Branch-normal Type Ia supernovae. At early times, the model photosphere is formed in material with velocities that are too high, indicating that the photosphere recedes too slowly or that SN 2011fe has a lower specific energy in the outer ~0.1 M_sun than does the model. We discuss several explanations for the discrepancies. Finally, we examine variations in both the spectral energy distribution and in the colors due to variations in the progenitor metallicity, which suggests that colors are only weak indicators for the progenitor metallicity, in the particular explosion model that we have studied. We do find that the flux in the U band is significantly higher at maximum light in the solar metallicity model than in the lower metallicity model and the lower metallicity model much better matches the observed spectrum.

Active galactic nuclei (AGNs) and quasars are important astrophysical objects to understand. Recently, microlensing observations have constrained the size of the quasar X-ray emission region to be of the order of 10 gravitational radii of the central supermassive black hole. For distances within a few gravitational radii, light paths are strongly bent by the strong gravity field of the central black hole. If the central black hole has nonzero angular momentum (spin), a photon's polarization plane will be rotated by the gravitational Faraday effect. The observed X-ray flux and polarization will then be influenced significantly by the strong gravity field near the source. Consequently, linear gravitational lensing theory is inadequate for such extreme circumstances. We present simple algorithms computing strong lensing effects of Kerr black holes, including effects on polarization. Our algorithms are realized in a program "KERTAP" in two versions: MATLAB and Python. The key ingredients of KERTAP are: a graphic user interface, a {\it backward} ray-tracing algorithm, a polarization propagator dealing with gravitational Faraday rotation, and algorithms computing observables such as flux magnification and polarization angles. Our algorithms can be easily realized in other programming languages such as FORTRAN, C, and C++. The MATLAB version of KERTAP is parallelized using the MATLAB Parallel Computing Toolbox and the Distributed Computing Server. The Python code was sped up using Cython and supports full implementation of MPI using 'mpi4py' package. As an example, we investigate the inclination angle dependence of the observed polarization and the strong lensing magnification of AGN X-ray emission. We conclude that it is possible to perform complex numerical-relativity-related computations using interpreted languages such as MATLAB and Python.

The effect of metallicity on the observed light of Type Ia supernovae (SNe Ia) could lead to systematic errors as the absolute magnitudes of local and distant SNe Ia are compared to measure luminosity distances and determine cosmological parameters. The UV light may be especially sensitive to metallicity, though different modeling methods disagree as to the magnitude, wavelength dependence, and even the sign of the effect. The outer density structure, ^56 Ni, and to a lesser degree asphericity, also impact the UV. We compute synthetic photometry of various metallicity-dependent models and compare to UV/optical photometry from the Swift Ultra-Violet/Optical Telescope. We find that the scatter in the mid-UV to near-UV colors is larger than predicted by changes in metallicity alone and is not consistent with reddening. We demonstrate that a recently employed method to determine relative abundances using UV spectra can be done using UVOT photometry, but we warn that accurate results require an accurate model of the cause of the variations. The abundance of UV photometry now available should provide constraints on models that typically rely on UV spectroscopy for constraining metallicity, density, and other parameters. Nevertheless, UV spectroscopy for a variety of SN explosions is still needed to guide the creation of accurate models. A better understanding of the influences affecting the UV is important for using SNe Ia as cosmological probes, as the UV light may test whether SNe Ia are significantly affected by evolutionary effects.

We obtained a time series of spectra covering the secondary maximum in the I-band of the bright Type Ia supernova 2014J in M82 with the TIGRE telescope. Comparing the observations with theoretical models calculated with the time dependent extension of the PHOENIX code, we identify the feature that causes the secondary maximum in the I-band light curve. Fe II 3d6(3D)4s-3d6(5D)4p and similar high excitation transitions produce a blended feature at 7500 {\AA}, which causes the rise of the light curve towards the secondary maximum. The series of observed spectra of SN 2014J and archival data of SN 2011fe confirm this conclusion. We further studied the plateau phase of the Rband light curve of SN 2014J and searched for features which contribute to the flux. The theoretical models do not clearly indicate a new feature that may cause the Rband plateau phase. However, Co II features in the range of 6500 - 7000 {\AA} and the Fe II feature of the I-band are clearly seen in the theoretical spectra, but do not appear to provide all of the flux necessary for the R-band plateau.

The question of the nature of the progenitor of Type Ia supernovae (SNe Ia) is important both for our detailed understanding of stellar evolution and for their use as cosmological probes of the dark energy. Much of the basic features of SNe Ia can be understood directly from the nuclear physics, a fact which Gerry would have appreciated. We present an overview of the current observational and theoretical situation and show that it not incompatible with most SNe Ia being the results of thermonuclear explosions near the Chandrasekhar mass.

We present UV through NIR broad-band photometry, and optical and NIR spectroscopy of Type Iax supernova 2012Z. The data set consists of both early and late-time observations, including the first late phase NIR spectrum obtained for a spectroscopically classified SN Iax. Simple model calculations of its bolometric light curve suggest SN 2012Z produced ~0.3 M_sun of (56)Ni, ejected about a Chandrasekhar mass of material, and had an explosion energy of ~10^51 erg, making it one of the brightest and most energetic SN Iax yet observed. The late phase NIR spectrum of SN 2012Z is found to broadly resemble similar epoch spectra of normal SNe Ia; however, like other SNe Iax, corresponding visual-wavelength spectra differ substantially compared to all supernova types. Constraints from the distribution of IMEs, e.g. silicon and magnesium, indicate that the outer ejecta did not experience significant mixing during or after burning, and the late phase NIR line profiles suggests most of the (56)Ni is produced during high density burning. The various observational properties of SN 2012Z are found to be consistent with the theoretical expectations of a Chandrasekhar mass white dwarf progenitor that experiences a pulsational delayed detonation, which produced several tenths of a solar mass of (56)Ni during the deflagration burning phase and little (or no) (56)Ni during the detonation phase. Within this scenario only a moderate amount of Rayleigh-Taylor mixing occurs both during the deflagration and fallback phase of the pulsation, and the layered structure of the IMEs is a product of the subsequent denotation phase. The fact that the SNe Iax population does not follow a tight brightness-decline relation similar to SNe Ia can then be understood in the framework of variable amounts of mixing during pulsational rebound and variable amounts of (56)Ni production during the early subsonic phase of expansion.

We present near-infrared synthetic spectra of a delayed-detonation hydrodynamical model and compare them to observed spectra of four normal type Ia supernovae ranging from day +56.5 to day +85. This is the epoch during which supernovae are believed to be undergoing the transition from the photospheric phase, where spectra are characterized by line scattering above an optically thick photosphere, to the nebular phase, where spectra consist of optically thin emission from forbidden lines. We find that most spectral features in the near-infrared can be accounted for by permitted lines of Fe II and Co II. In addition, we find that [Ni II] fits the emission feature near 1.98 {\mu}m, suggesting that a substantial mass of 58Ni exists near the center of the ejecta in these objects, arising from nuclear burning at high density. A tentative identification of Mn II at 1.15 {\mu}m may support this conclusion as well.

We present the first J-band spectrum of Mrk 231, which reveals a large \ion{He}{1}*$\lambda 10830$ broad absorption line with a profile similar to that of the well-known \ion{Na}{1} broad absorption line. Combining this spectrum with optical and UV spectra from the literature, we show that the unusual reddening noted by \citet{veilleux13} is explained by a reddening curve like those previously used to explain low values of total-to-selective extinction in SNe Ia. The nuclear starburst may be the origin and location of the dust. Spatially-resolved emission in the broad absorption line trough suggests nearly full coverage of the continuum emission region. The broad absorption lines reveal higher velocities in the \ion{He}{1}* lines (produced in the quasar-photoionized \ion{H}{2} region) compared with the \ion{Na}{1} and \ion{Ca}{2} lines (produced in the corresponding partially-ionized zone). {\it Cloudy} simulations show that a density increase is required between the \ion{H}{2} and partially-ionized zones to produce ionic column densities consistent with the optical and IR absorption line measurements and limits, and that the absorber lies $\sim 100\rm \, pc$ from the central engine. These results suggest that the \ion{He}{1}* lines are produced in an ordinary quasar BAL wind that impacts upon, compresses, and accelerates the nuclear starburst's dusty effluent (feedback in action), and the \ion{Ca}{2} and \ion{Na}{1} lines are produced in this dusty accelerated gas. This unusual circumstance explains the rarity of \ion{Na}{1} absorption lines; without the compression along our line of sight, Mrk~231 would appear as an ordinary FeLoBAL.

Multi-level non-local thermodynamic equilibrium (NLTE) radiation transfer calculations have become standard throughout the stellar atmospheres community and are applied to all types of stars as well as dynamical systems such as novae and supernovae. Even today spherically symmetric 1D calculations with full physics are computationally intensive. We show that full NLTE calculations can be done with fully 3 dimensional (3D) radiative transfer. With modern computational techniques and current massive parallel computational resources, full detailed solution of the multi-level NLTE problem coupled to the solution of the radiative transfer scattering problem can be solved without sacrificing the micro physics description. We extend the use of a rate operator developed to solve the coupled NLTE problem in spherically symmetric 1D systems. In order to spread memory among processors we have implemented the NLTE/3D module with a hierarchical domain decomposition method that distributes the NLTE levels, radiative rates, and rate operator data over a group of processes so that each process only holds the data for a fraction of the voxels. Each process in a group holds all the relevant data to participate in the solution of the 3DRT problem so that the 3DRT solution is parallelized within a domain decomposition group. We solve a spherically symmetric system in 3D spherical coordinates in order to directly compare our well-tested 1D code to the 3D case. We compare three levels of tests: a) a simple H+He test calculation, b) H+He+CNO+Mg, c) H+He+Fe. The last test is computationally large and shows that realistic astrophysical problems are solvable now, but they do require significant computational resources. With presently available computational resources it is possible to solve the full 3D multi-level problem with the same detailed micro-physics as included in 1D modeling.

A comprehensive set of optical and near-infrared (NIR) photometry and spectroscopy is presented for the faint and fast 2008ha-like supernova (SN) 2010ae. Contingent on the adopted value of host extinction, SN 2010ae reached a peak brightness of -13.8 > M_V > -15.3 mag, while modeling of the UVOIR light curve suggests it produced 0.003--0.007 M_sun of (56)Ni, ejected 0.30--0.60 M_sun of material, and had an explosion energy of 0.04--0.30x10^51 erg. The values of these explosion parameters are similar to the peculiar SN 2008ha --for which we also present previously unpublished early phase optical and NIR light curves-- and places these two transients at the faint end of the 2002cx-like SN population. Detailed inspection of the post maximum NIR spectroscopic sequence indicates the presence of a multitude of spectral features, which are identified through SYNAPPS modeling to be mainly attributed to Co II. Comparison with a collection of published and unpublished NIR spectra of other 2002cx-like SNe, reveals that a Co II footprint is ubiquitous to this subclass of transients, providing a link to Type Ia SNe. A visual-wavelength spectrum of SN 2010ae obtained at +252 days past maximum shows a striking resemblance to a similar epoch spectrum of SN 2002cx. However, subtle differences in the strength and ratio of calcium emission features, as well as diversity among similar epoch spectra of other 2002cx-like SNe indicates a range of physical conditions of the ejecta, highlighting the heterogeneous nature of this peculiar class of transients.

We present predictions for hydrogen and helium emission line luminosities from circumstellar matter around Type Ia supernovae (SNe Ia) using time dependent photoionization modeling. ESO/VLT optical echelle spectra of the SN Ia 2000cx were taken before and up to 70 days after maximum. We detect no hydrogen and helium lines, and place an upper limit on the mass loss rate for the putative wind of less than 1.3EE{-5} solar masses per year, assuming a speed of 10 km/s and solar abundances for the wind. In a helium-enriched case, the best line to constrain the mass loss would be He I 10,830 A. We confirm the details of interstellar Na I and Ca II absorption towards SN 2000cx as discussed by Patat et al., but also find evidence for 6613.56 A Diffuse Interstellar Band (DIB) absorption in the Milky Way. We discuss measurements of the X-ray emission from the interaction between the supernova ejecta and the wind and we re-evaluate observations of SN 1992A obtained 16 days after maximum by Schlegel & Petre. We find an upper limit of 1.3EE{-5} solar masses per year. These results, together with the previous observational work on the normal SNe Ia 1994D and 2001el, disfavour a symbiotic star in the upper mass loss rate regime from being the likely progenitor scenario for these SNe. To constrain hydrogen in late time spectra, we present ESO/VLT and ESO/NTT optical and infrared observations of SNe Ia 1998bu and 2000cx 251-388 days after maximum. We see no hydrogen line emission in SNe 1998bu and 2000cx at these epochs, and we argue from modeling that the mass of such hydrogen-rich gas must be less than 0.03 solar masses for both supernovae. Comparing similar upper limits with recent models of Pan et al., it seems hydrogen-rich donors with a separation of less than 5 times the radius of the donor may be ruled out for the five SNe Ia 1998bu, 2000cx, 2001el, 2005am and 2005cf.

Radiative transfer studies of Type Ia supernovae (SNe Ia) hold the promise of constraining both the time-dependent density profile of the SN ejecta and its stratification by element abundance which, in turn, may discriminate between different explosion mechanisms and progenitor classes. Here we present a detailed analysis of Hubble Space Telescope ultraviolet (UV) and ground-based optical spectra and light curves of the SN Ia SN 2010jn (PTF10ygu). SN 2010jn was discovered by the Palomar Transient Factory (PTF) 15 days before maximum light, allowing us to secure a time-series of four UV spectra at epochs from -11 to +5 days relative to B-band maximum. The photospheric UV spectra are excellent diagnostics of the iron-group abundances in the outer layers of the ejecta, particularly those at very early times. Using the method of 'Abundance Tomography' we have derived iron-group abundances in SN 2010jn with a precision better than in any previously studied SN Ia. Optimum fits to the data can be obtained if burned material is present even at high velocities, including significant mass fractions of iron-group elements. This is consistent with the slow decline rate (or high 'stretch') of the light curve of SN 2010jn, and consistent with the results of delayed-detonation models. Early-phase UV spectra and detailed time-dependent series of further SNe Ia offer a promising probe of the nature of the SN Ia mechanism.

Recent quasar microlensing observations have constrained the sizes of X-ray emission regions to be within about 10 gravitational radii of the central supermassive black hole. Therefore, the X-ray emission from lensed quasars is first strongly lensed by the black hole before it is lensed by the foreground galaxy and star fields. We present a scheme that combines the initial strong lensing of a Kerr black hole with standard linearized microlensing by intervening stars. We find that X-ray microlensed light curves incorporating Kerr strong gravity can differ significantly from standard curves. The amplitude of the fluctuations in the light curves can increase or decrease by ~0.65-0.75 mag by including Kerr strong gravity. Larger inclination angles give larger amplitude fluctuations in the microlensing light curves. Consequently, current X-ray microlensing observations might have under or overestimated the sizes of the X-ray emission regions. We estimate this bias using a simple metric based on the amplitude of magnitude fluctuations. The half light radius of the X-ray emission region can be underestimated up to ~50% or overestimated up to ~20%. Underestimates are found in most situations we have investigated. The only exception is for a disk with large spin, radially flat emission profile, and observed nearly face on, where an overestimate is found. Thus, more accurate microlensing size constraints should be obtainable by including Kerr lensing. The caustic crossing time can differ by months after including Kerr strong gravity. A simultaneous monitoring of gravitational lensed quasars in both X-ray and optical bands with densely sampled X-ray light curves might reveal this feature. We conclude that it should be possible to constrain important parameters such as inclination angles and black hole spins from combined Kerr and microlensing effects.

Recent quasar microlensing observations have constrained the X-ray emission sizes of quasars to be about 10 gravitational radii, one order of magnitude smaller than the optical emission sizes. Using a new ray-tracing code for the Kerr space-time, we find that the observed X-ray flux is strongly influenced by the gravity field of the central black hole, even for observers at moderate inclination angles. We calculate inclination-dependent flux profiles of active galactic nuclei in the optical and X-ray bands by combining the Kerr lensing and projection effects for future references. We further study the dependence of the X-ray-to-optical flux ratio on the inclination angle caused by differential lensing distortion of the X-ray and optical emission, assuming several corona geometries. The strong lensing X-ray-to-optical magnification ratio can change by a factor of ~10 for normal quasars in some cases, and another factor of ~10 for broad absorption line quasars (BALs) and obscured quasars. Comparing our results with the observed distributions in normal and broad absorption line quasars, we find that the inclination angle dependence of the magnification ratios can change the X-ray-to-optical flux ratio distributions significantly. In particular, the mean value of the spectrum slope parameter $\alpha_{ox},$ $0.3838\log F_{2 keV}/F_{2500 {\AA}}$, can differ by ~0.1-0.2 between normal and broad absorption line quasars, depending on corona geometries, suggesting larger intrinsic absorptions in BALs.

3-D astrophysical atmospheres will have random velocity fields. We seek to combine the methods we have developed for solving the 1-D problem with arbitrary flows to those that we have developed for solving the fully 3-D relativistic radiative transfer problem in the case of monotonic flows. The methods developed in the case of 3-D atmospheres with monotonic flows, solving the fully relativistic problem along curves defined by an affine parameter, are very flexible and can be extended to the case of arbitrary velocity fields in 3-D. Simultaneously, the techniques we developed for treating the 1-D problem with arbitrary velocity fields are easily adapted to the 3-D problem. The algorithm we present allows the solution of 3-D radiative transfer problems that include arbitrary wavelength couplings. We use a quasi-analytic formal solution of the radiative transfer equation that significantly improves the overall computation speed. We show that the approximate lambda operator developed in previous work gives good convergence, even neglecting wavelength coupling. Ng acceleration also gives good results. We present tests that are of similar resolution to what has been presented using Monte-Carlo techniques, thus our methods will be applicable to problems outside of our test setup. Additional domain decomposition parallelization strategies will be explored in future work.

Context. Time-dependent, 3D radiation transfer calculations are important for the modeling of a variety of objects, from supernovae and novae to simulations of stellar variability and activity. Furthermore, time-dependent calculations can be used to obtain a 3D radiative equilibrium model structure via relaxation in time. Aims. We extend our 3D radiative transfer framework to include direct time dependence of the radiation field; i.e., the $\partial I/ \partial t$ terms are fully considered in the solution of radiative transfer problems. Methods. We build on the framework that we have described in previous papers in this series and develop a subvoxel method for the $\partial I/\partial t$ terms. Results. We test the implementation by comparing the 3D results to our well tested 1D time dependent radiative transfer code in spherical symmetry. A simple 3D test model is also presented. Conclusions. The 3D time dependent radiative transfer method is now included in our 3D RT framework and in PHOENIX/3D.

In supernova spectroscopy relatively little attention has been given to the properties of optically thick spectral lines in epochs following the photosphere's recession. Most treatments and analyses of post-photospheric optical spectra of supernovae assume that forbidden-line emission comprises most if not all spectral features. However, evidence exists which suggests that some spectra exhibit line profiles formed via optically thick resonance-scattering even months or years after the supernova explosion. To explore this possibility we present a geometrical approach to supernova spectrum formation based on the "Elementary Supernova" model, wherein we investigate the characteristics of resonance-scattering in optically thick lines while replacing the photosphere with a transparent central core emitting non-blackbody continuum radiation, akin to the optical continuum provided by decaying 56Co formed during the explosion. We develop the mathematical framework necessary for solving the radiative transfer equation under these conditions, and calculate spectra for both isolated and blended lines. Our comparisons with analogous results from the Elementary Supernova code SYNOW reveal several marked differences in line formation. Most notably, resonance lines in these conditions form P Cygni-like profiles, but the emission peaks and absorption troughs shift redward and blueward, respectively, from the line's rest wavelength by a significant amount, despite the spherically symmetric distribution of the line optical depth in the ejecta. These properties and others that we find in this work could lead to misidentification of lines or misattribution of properties of line-forming material at post-photospheric times in supernova optical spectra.

The nearby Type Ia supernova SN 2011fe in M101 (cz=241 km s^-1) provides a unique opportunity to study the early evolution of a "normal" Type Ia supernova, its compositional structure, and its elusive progenitor system. We present 18 high signal-to-noise spectra of SN 2011fe during its first month beginning 1.2 days post-explosion and with an average cadence of 1.8 days. This gives a clear picture of how various line-forming species are distributed within the outer layers of the ejecta, including that of unburned material (C+O). We follow the evolution of C II absorption features until they diminish near maximum light, showing overlapping regions of burned and unburned material between ejection velocities of 10,000 and 16,000 km s^-1. This supports the notion that incomplete burning, in addition to progenitor scenarios, is a relevant source of spectroscopic diversity among SNe Ia. The observed evolution of the highly Doppler-shifted O I 7774 absorption features detected within five days post-explosion indicate the presence of O I with expansion velocities from 11,500 to 21,000 km s^-1. The fact that some O I is present above C II suggests that SN 2011fe may have had an appreciable amount of unburned oxygen within the outer layers of the ejecta.

We present a study of the peculiar Type Ia supernova 2001ay (SN 2001ay). The defining features of its peculiarity are: high velocity, broad lines, and a fast rising light curve, combined with the slowest known rate of decline. It is one magnitude dimmer than would be predicted from its observed value of Delta-m15, and shows broad spectral features. We base our analysis on detailed calculations for the explosion, light curves, and spectra. We demonstrate that consistency is key for both validating the models and probing the underlying physics. We show that this SN can be understood within the physics underlying the Delta-m15 relation, and in the framework of pulsating delayed detonation models originating from a Chandrasekhar mass, white dwarf, but with a progenitor core composed of 80% carbon. We suggest a possible scenario for stellar evolution which leads to such a progenitor. We show that the unusual light curve decline can be understood with the same physics as has been used to understand the Delta-m15 relation for normal SNe Ia. The decline relation can be explained by a combination of the temperature dependence of the opacity and excess or deficit of the peak luminosity, alpha, measured relative to the instantaneous rate of radiative decay energy generation. What differentiates SN 2001ay from normal SNe Ia is a higher explosion energy which leads to a shift of the Ni56 distribution towards higher velocity and alpha < 1. This result is responsible for the fast rise and slow decline. We define a class of SN 2001ay-like SNe Ia, which will show an anti-Phillips relation.

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.

We present photometry and spectroscopy of the type IIP supernova 2009bw in UGC 2890 from few days after the outburst to 241 days. The light curve of SN 2009bw during the photospheric phase is similar to that of normal SNe IIP but with brighter peak and plateau (Mmax R = -17.82 mag, Mplateau R = -17.37 mag). The luminosity drop from the photospheric to the nebular phase is one of the fastest ever observed, ~2.2 mag in about 13 days. The radioactive tail of the bolometric light curve indicates that the amount of ejected 56 Ni is \approx 0.022 M\odot. The photospheric spectra reveal high velocity lines of H{\alpha} and H{\beta} until about 105 days after the shock breakout, suggesting a possible early interaction between the SN ejecta and pre-existent circumstellar material, and the presence of CNO elements. By modeling the bolometric light curve, ejecta expansion velocity and photospheric temperature, we estimate a total ejected mass of 8-12M\odot, a kinetic energy of ~0.3 foe and an initial radius of ~ 3.6 - 7 \times 10^13 cm.

We compare and analyze a time series of spectral observations obtained during the first 30 days of evolution of SN 2007od with the non-LTE code PHOENIX. Despite some spectroscopic particularities in the Balmer features, this supernova appears to be a normal Type II, and the fits proposed are generally in good agreement with the observations. As a starting point we have carried out an analysis with the parameterized synthetic spectrum code SYNOW to confirm line identifications and to highlight differences between the results of the two codes. The analysis computed using PHOENIX suggests the presence of a high velocity feature in H{\beta} and an H{\alpha} profile reproduced with a density profile steeper than that of the other elements. We also show a detailed analysis of the ions velocities of the 6 synthetic spectra. The distance is estimated for each epoch with the Spectral-fitting Expanding Atmosphere Method (SEAM). Consistent results are found using all the spectra which give the explosion date of JD 2454403 (29 October, 2007) and a distance modulus \mu = 32.2 \pm 0.3.

Aims. With our time-dependent model atmosphere code PHOENIX, our goal is to simulate light curves and spectra of hydrodynamical models of all types of supernovae. In this work, we simulate near-infrared light curves of SNe Ia and confirm the cause of the secondary maximum. Methods. We apply a simple energy solver to compute the evolution of an SN Ia envelope during the free expansion phase. Included in the solver are energy changes due to expansion, the energy deposition of {\gamma}-rays and interaction of radiation with the material. Results. We computed theoretical light curves of several SN Ia hydrodynamical models in the I, J, H, and K bands and compared them to the observed SN Ia light curves of SN 1999ee and SN 2002bo. By changing a line scattering parameter in time, we obtained quite reasonable fits to the observed near-infrared light curves. This is a strong hint that detailed NLTE effects in IR lines have to be modeled, which will be a future focus of our work. Conclusions. We found that IR line scattering is very important for the near-infrared SN Ia light curve modeling. In addition, the recombination of Fe III to Fe II and of Co III to Co II is responsible for the secondary maximum in the near-infrared bands. For future work the consideration of NLTE for all lines (including the IR subordinate lines) will be crucial.

We present an analysis of high precision V light curves (LC) for 18 local Type Ia Supernovae, SNe Ia, as obtained with the same telescope and setup at the Las Campanas Observatory (LCO). This homogeneity provides an intrinsic accuracy a few hundreds of a magnitude both with respect to individual LCs and between different objects. Based on the Single Degenerate Scenario, SD, we identify patterns which have been predicted by model calculations as signatures of the progenitor and accretion rate which change the explosion energy and the amount of electron capture, respectively. Using these templates as principle components and the overdetermined system of SN pairs, we reconstruct the properties of progenitors and progenitor systems. All LCO SNe Ia follow the brightness decline relation but 2001ay. After subtraction of the two components, the remaining scatter is reduced to 0.01-0.03m. Type SNe Ia seem to originate from progenitors with Main Sequence masses of 3Mo with the exception of two subluminous SNe Ia with < 2Mo. The component analysis indicates a wide range of accretion rates in the progenitor systems closing the gap to accretion induced collapses (AIC). SN1991t-like objects show differences in $dm15$ but no tracers of our secondary parameters. This may point to a different origin such as DD-Scenario or the Pulsating Delayed Detonations. SN2001ay does not follow the decline relation. It can be understood in the framework of C-rich WDs, and this group may produce an anti-Phillips relation. We suggest that this may be a result of a common envelope phase and mixing during central He burning as in SN1987A.

We discuss an implementation of our 3D radiative transfer (3DRT) framework with the OpenCL paradigm for general GPU computing. We implement the kernel for solving the 3DRT problem in Cartesian coordinates with periodic boundary conditions in the horizontal $(x,y)$ plane, including the construction of the nearest neighbor $\Lstar$ and the operator splitting step. We present the results of a small and a large test case and compare the timing of the 3DRT calculations for serial CPUs and various GPUs. The latest available GPUs can lead to significant speedups for both small and large grids compared to serial (single core) computations.

Aims. We present the first theoretical SN Ia light curves calculated with the time-dependent version of the general purpose model atmosphere code PHOENIX. Our goal is to produce light curves and spectra of hydro models of all types of supernovae. Methods. We extend our model atmosphere code PHOENIX to calculate type Ia supernovae light curves. A simple solver was imple- mented which keeps track of energy conservation in the atmosphere during the free expansion phase. Results. The correct operation of the new additions to PHOENIX were verified in test calculations. Furthermore, we calculated theo- retical light curves and compared them to the observed SN Ia light curves of SN 1999ee and SN 2002bo. We obtained LTE as well as NLTE model light curves. Conclusions. We have verified the correct operation of our extension into the time domain. We have calculated the first SN Ia model light curves using PHOENIX in both LTE and NLTE. For future work the infrared model light curves need to be further investigated.

Spectroscopic analyses of Type Ia supernovae have shown there exist four spectroscopic groups---cools, broad line, shallow silicon, and core normal---defined by the widths of the Si II features at 5972 Angstroms and 6355 Angstroms. 1991bg-likes are classified as "cools". Cools are dim, undergo a rapid decline in luminosity, and produce significantly less nickel than normal Type Ia supernovae. They also have an unusually deep and wide trough in their spectra around 4200 Angstroms and a relatively strong Si II absorption attributed to the line at 5972 Angstroms. We examine the spectra of supernova (SN) 1991bg and the cools SN 1997cn, SN 1999by, and SN 2005bl using the highly parameterized synthetic spectrum code SYNOW, and find general agreement with similar spectroscopic studies. Our analysis reveals that this group of supernovae is fairly homogeneous, with many of the blue spectral features well fit by Fe II. The nature of the spectroscopic commonalities and the variations in the class are discussed. Finally, we examine intermediates such as SN 2004eo and discuss the spectroscopic subgroup distribution of Type Ia supernovae.

Theoretical arguments and observations suggest that the atmospheres of Brown Dwarfs and planets are very dynamic on chemical and on physical time scales. The modelling of such substellar atmospheres has, hence, been much more demanding than initially anticipated. This Splinter (http://star-www.st-and.ac.uk/~ch80/CS16/MultiDSplinter_CS16.html) has combined new developments in atmosphere modelling, with novel observational techniques, and new challenges arising from planetary and space weather observations.

A solution of the radiative-transfer problem in 3D with arbitrary velocity fields in the Eulerian frame is presented. The method is implemented in our 3D radiative transfer framework and used in the PHOENIX/3D code. It is tested by comparison to our well- tested 1D co-moving frame radiative transfer code, where the treatment of a monotonic velocity field is implemented in the Lagrangian frame. The Eulerian formulation does not need much additional memory and is useable on state-of-the-art computers, even large-scale applications with 1000's of wavelength points are feasible.

We present synthetic spectral fits of the typical Type Ib SN 1999dn and the Hydrogen Rich Ib SN 2000H using the generalized non-local thermodynamic equilibrium stellar atmospheres code \phx. We fit model spectra to five epochs of SN 1999dn ranging from ten days pre-maximum light to 17 days post-maximum light and the two earliest epochs of SN 2000H available, maximum light and six days post-maximum. Our goal is to investigate the possibility of hydrogen in Type Ib Supernovae (SNe Ib), specifically a feature around 6200\AA\ which has previously been attributed to high velocity H-alpha. In earlier work on SN 1999dn we found the most plausible alternative to H-alpha to be a blend of Si II and Fe II lines which can be adjusted to fit by increasing the metallicity. Our models are simple; they assume a powerlaw density profile with radius, homologous expansion, and solar compositions. The helium core is produced by burning 4H --> He in order to conserve nucleon number. For models with hydrogen the outer skin of the model consists of a shell of solar composition. The hydrogen mass of the standard solar composition shell is M_H less than about 0.001 times the mass of the sun in SN 1999dn and M_H less than about 0.2 times the mass of the sun for SN 2000H. Our models fit the observed spectra reasonably well, successfully reproducing most features including the characteristic He I absorptions. The hydrogen feature in SN 1999dn is clear, but much more pronounced in SN 2000H. We discuss a possible evolutionary scenario that accounts for the dichotomy in the hydrogen shell mass between these two supernovae.

Hydrogen recombination is one of the most important atomic processes in many astrophysical objects such as Type II supernova (SN~II) atmospheres, the high redshift universe during the cosmological recombination era, and H II regions in the interstellar medium. Accurate predictions of the ionization fraction can be quite different from those given by a simple solution if one takes into account many angular momentum sub-states, non-resonant processes, and calculates the rates of all atomic processes from the solution of the radiative transfer equation instead of using a Planck function under the assumption of thermal equilibrium. We use the general purpose model atmosphere code PHOENIX 1D to compare how the fundamental probabilities such as the photo-ionization probability, the escape probability, and the collisional de-excitation probability are affected by the presence of other metals in the environment, multiple angular momentum sub-states, and non-resonant processes. Our comparisons are based on a model of SN 1999em, a SNe Type II, 20 days after its explosion.

High-quality observations of $B$ and $V$ light curves obtained at Las Campanas Observatory for local Type Ia Supernovae (SNe Ia) show clear evidence that SNe Ia with the same brightness decline or stretch may have systematic and independent deviations at times < 5 days before and at times > 30 days after maximum light. This suggests the existence of two independent secondary parameters which control the shape of SN Ia light curves in addition to the brightness decline relation. stretch. The differences are consistent in morphology of the time dependence and size with predictions by models within the delayed detonation scenario. The secondary parameters may reflect two independent physical effects caused by variations in the progenitor and accretion rates, and link the LC variations in shape with the intrinsic, absolute brightness.

We demonstrate the application of our 3D radiative transfer framework in the model atmosphere code PHOENIX/3D for a number of spectrum synthesis calculations for very different conditions. The 3DRT framework discussed in the previous papers of this series was added to our general-purpose model atmosphere code PHOENIX/1D and an extended 3D version PHOENIX/3D was created. The \phxT code is parallelized via the MPI library using a hierarchical domain decomposition and displays very good strong scaling. We present the results of several test cases for widely different atmosphere conditions and compare the 3D calculations with equivalent 1D models to assess the internal accuracy of the 3D modeling. In addition, we show the results for a number of parameterized 3D structures. With presently available computational resources it is possible to solve the full 3D radiative transfer (including scattering) problem with the same micro-physics as included in 1D modeling.

NLTE radiative transfer calculations of differentially expanding supernovae atmospheres are computationally intensive and are almost universally performed in time-independent snapshot mode. The validity of the steady-state approximation in the rate equations has recently been questioned. We calculate the effective recombination time of hydrogen in SN II using our general purpose model atmosphere code PHOENIX. While we find that the recombination time for the conditions of SNe II at early times is increased over the classical value for the case of a simple hydrogen model atom with energy levels corresponding to just the first 2 principle quantum numbers, the classical value of the recombination time is recovered in the case of a multi-level hydrogen atom. We also find that the recombination time at most optical depths is smaller in the case of a multi-level atom than for a simple two-level hydrogen atom. We find that time dependence in the rate equations is important in the early epochs of a supernova's lifetime. The changes due to the time dependent rate equation (at constant input luminosity) are manifested in physical parameters such as the level populations which directly affects the spectra. The H-alpha profile is affected by the time dependent rate equations at early times. At later times time dependence does not significantly modify the level populations and therefore the H-alpha profile is roughly independent of whether the steady-state or time-dependent approach is used.

A solution of the radiative-transfer problem in arbitrary velocity fields introduced in a previous paper, has limitations in its applicability. For large-scale applications, the methods described also require large memory sets that are commonly not available to state-of-the-art computing hardware. In this work, we modify the algorithm to allow the computation of large-scale problems. We reduce the memory footprint via a domain decomposition. By introducing iterative Gauss-Seidel type solvers, we improve the speed of the overall computation. Because of the domain decomposition, the new algorithm requires the use of parallel-computing systems. The algorithm that we present permits large-scale solutions of radiative-transfer problems that include arbitrary wavelength couplings. In addition, we discover a quasi-analytic formal solution of the radiative transfer that significantly improves the overall computation speed. More importantly, this method ensures that our algorithm can be applied to multi-dimensional Lagrangian radiative-transfer calculations. In multi-dimensional atmospheres, velocity fields are in general chaotic ensuring that the inclusion of arbitrary wavelength couplings are mandatory.

Using the methods of general relativity Lindquist derived the radiative transfer equation that is correct to all orders in v/c. Mihalas developed a method of solution for the important case of monotonic velocity fields with spherically symmetry. We have developed the generalized atmosphere code PHOENIX, which in 1-D has used the framework of Mihalas to solve the radiative transfer equation (RTE) in 1-D moving flows. We describe our recent work including 3-D radiation transfer in PHOENIX and particularly including moving flows exactly using a novel affine method. We briefly discuss quantitative spectroscopy in supernovae.

Aims. We present first results and tests of a time-dependent extension to the general purpose model atmosphere code PHOENIX. We aim to produce light curves and spectra of hydro models for all types of supernovae. Methods. We extend our model atmosphere code PHOENIX to solve time-dependent non-grey, NLTE, radiative transfer in a special relativistic framework. A simple hydrodynamics solver was implemented to keep track of the energy conservation of the atmosphere during free expansion. Results. The correct operation of the new additions to PHOENIX were verified in test calculations. Conclusions. We have shown the correct operation of our extension to time-dependent radiative transfer and will be able to calculate supernova light curves and spectra in future work.

Observations and theoretical calculations have shown the importance of non-spherically symmetric structures in supernovae. Thus, the interpretation of observed supernova spectra requires the ability to solve the transfer equation in 3-D moving atmospheres. We present an implementation of the solution of the radiative transfer equation in 3-D homologously expanding atmospheres in spherical coordinates. The implementation is exact to all orders in v/c. We use a new affine method that makes use of the fact that photons travel on straight lines. We compare our results in 3-D for spherically symmetric test problems with high velocity fields and find excellent agreement. Our well-tested 1-D results are based on methods where the momentum directions vary along the characteristic (co-moving momentum directions). Thus, we are able to verify both the analytic framework and its numerical implementation. Additionally, we have been able to test the parallelization over characteristics. Using 512^2 momentum angles we ran the code on 16,384 Opteron processors and achieved excellent scaling. It is now possible to calculate synthetic spectra from realistic 3D hydro simulations. This should open an era of progress in hydro modeling, similar to that that occurred in the 1980s when 1-D models were confronted with synthetic spectra.

We extend our framework for 3D radiative transfer calculations with a non-local operator splitting methods along (full) characteristics to spherical and cylindrical coordinate systems. These coordinate systems are better suited to a number of physical problems than Cartesian coordinates. The scattering problem for line transfer is solved via means of an operator splitting (OS) technique. The formal solution is based on a full characteristics method. The approximate $\Lambda$ operator is constructed considering nearest neighbors exactly. The code is parallelized over both wavelength and solid angle using the MPI library. We present the results of several test cases with different values of the thermalization parameter for the different coordinate systems. The results are directly compared to 1D plane parallel tests. The 3D results agree very well with the well-tested 1D calculations.

In the next decade Type Ia supernovae (SNe Ia) will be used to test theories predicting changes in the Dark Energy equation of state with time. Ultimately this requires a dedicated space mission like JDEM. SNe Ia are mature cosmological probes --- their limitations are well characterized, and a path to improvement is clear. Dominant systematic errors include photometric calibration, selection effects, reddening, and population-dependent differences. Building on past lessons, well-controlled new surveys are poised to make strides in these areas: the Palomar Transient Factory, Skymapper, La Silla QUEST, Pan-STARRS, the Dark Energy Survey, LSST, and JDEM. They will obviate historical calibrations and selection biases, and allow comparisons via large subsamples. Some systematics follow from our ignorance of SN Ia progenitors, which there is hope of determining with SN Ia rate studies from 0<z<4. Aside from cosmology, SNe Ia regulate galactic and cluster chemical evolution, inform stellar evolution, and are laboratories for extreme physics. Essential probes of SNe Ia in these contexts include spectroscopy from the UV to the IR, X-ray cluster and SN remnant observations, spectropolarimetry, and advanced theoretical studies. While there are an abundance of discovery facilities planned, there is a deficit of follow-up resources. Living in the systematics era demands deep understanding rather than larger statistics. NOAO ReSTAR initiative to build 2-4m telescopes would provide necessary follow-up capability. Finally, to fully exploit LSST, well-matched wide-field spectroscopic capabilities are desirable.

- Constraining the cosmological parameters and understanding Dark Energy have tremendous implications for the nature of the Universe and its physical laws. - The pervasive limit of systematic uncertainties reached by cosmography based on Cepheids and Type Ia supernovae (SNe Ia) warrants a search for complementary approaches. - Type II SNe have been shown to offer such a path. Their distances can be well constrained by luminosity-based or geometric methods. Competing, complementary, and concerted efforts are underway, to explore and exploit those objects that are extremely well matched to next generation facilities. Spectroscopic follow-up will be enabled by space- based and 20-40 meter class telescopes. - Some systematic uncertainties of Type II SNe, such as reddening by dust and metallicity effects, are bound to be different from those of SNe Ia. Their stellar progenitors are known, promising better leverage on cosmic evolution. In addition, their rate - which closely tracks the ongoing star formation rate - is expected to rise significantly with look- back time, ensuring an adequate supply of distant examples. - These data will competitively constrain the dark energy equation of state, allow the determination of the Hubble constant to 5%, and promote our understanding of the processes involved in the last dramatic phases of massive stellar evolution.

We present synthetic spectra around maximum for the type II supernova SN 2003Z, which was first detected on January 29.7 2003. Comparison with observed spectra aim at the determination of physical parameters for SN 2003Z. Synthetic spectra are calculated with our stellar atmosphere code PHOENIX. It solves the special relativistic equation of radiative transfer, including large NLTE-calculations and line blanketing by design, in 1-dimensional spherical symmetry. The observed spectra were obtained at the 3.5 meter telescope at Calar Alto. The TWIN instrument was used so that a spectral range from about 3600 to 7500 Angstroem was covered. The spectra were taken on Feb. 4, 5, 9, and 11, 2003. The physical parameters of the models give the luminosities, a range of possible velocity profiles for the SN, an estimate of the colour excess, and the observed metalicity.

We present a new formal solution of the Lagrangian equation of radiative transfer that is useful in solving the equation of radiative transfer in the presence of arbitrary velocity fields. Normally a term due to the inclusion of the wavelength derivative in the Lagrangian equation of radiative transfer is associated with a generalised opacity. In non-monotonic velocity fields, this generalised opacity may become negative. To ensure that the opacity remains positive, this term of the derivative is included in the formal solution of the radiative transfer problem. The new definition of the generalised opacity allows for a new solution of the equation of radiative transfer in the presence of velocity fields. It is especially useful for arbitrary velocity fields, where it effectively prevents the occurrences of negative generalised opacities and still allows the explicit construction of the Lambda-operator of the system needed for an accelerated Lambda-iteration. We performed test calculations, where the results of old, established solutions were compared with the new solution. The relative deviations never exceeded 1% and so the new solution is indeed suitable for use in radiative-transfer modelling. Non-monotonic velocity fields along photon paths frequently occur in three-dimensional hydrodynamical models of astrophysical atmospheres. Therefore, the formal solution will be of use for multidimensional radiative transfer and has immediate applications in the modelling of pulsating stars and astrophysical shock fronts.

A comparative study of optical spectra of Type Ia supernovae (SNe Ia) is extended, in the light of new data. The discussion is framed in terms of the four groups defined in previous papers of this series: core normal (CN); broad line (BL); cool (CL); and shallow silicon (SS). Emerging features of the SN Ia spectroscopic diversity include evidence (1) that extreme CL SN 1991bg-likes are not a physically distinct subgroup and (2) for the existence of a substantial number of SN 1999aa-like SSs that are very similar to each other and distinguishable from CN even as late as three weeks after maximum light. SN 1999aa-likes may be relatively numerous, yet not a physically distinct subgroup. The efficacy of quantitative spectroscopic subclassification of SNe Ia based on the equivalent widths of absorption features near 5750 A and 6100 A near maximum light is discussed. The absolute magnitude dispersion of a small sample of CNs is no larger than the characteristic absolute magnitude uncertainty.

The first paper in this series explored the effects of altering the chemical mixture of the stellar population on an element by element basis on stellar evolutionary tracks and isochrones to the end of the red giant branch. This paper extends the discussion by incorporating the fully consistent synthetic stellar spectra with those isochrone models in predicting integrated colors, Lick indices, and synthetic spectra. Older populations display element ratio effects in their spectra at higher amplitude than younger populations. In addition, spectral effects in the photospheres of stars tend to dominate over effects from isochrone temperatures and lifetimes, but, further, the isochrone-based effects that are present tend to fall along the age-metallicity degeneracy vector, while the direct stellar spectral effects usually show considerable orthogonality.

We present a general method to solve radiative transfer problems including scattering in the continuum as well as in lines in 3D configurations with periodic boundary conditions. he scattering problem for line transfer is solved via means of an operator splitting (OS) technique. The formal solution is based on a full characteristics method. The approximate $\Lambda$ operator is constructed considering nearest neighbors exactly. The code is parallelized over both wavelength and solid angle using the MPI library. We present the results of several test cases with different values of the thermalization parameter and two choices for the temperature structure. The results are directly compared to 1D plane parallel tests. The 3D results agree very well with the well-tested 1D calculations.

We present the first results from our GALEX program designed to obtain ultraviolet (UV) spectroscopy of nearby core-collapse supernovae (SNe). Our first target, SN 2005ay in the nearby galaxy NGC 3938, is a typical member of the II-P SN subclass. Our spectra show remarkable similarity to those of the prototypical type II-P event SN 1999em, and resemble also Swift observations of the recent type II-P event SN 2005cs. Taken together, the observations of these three events trace the UV spectral evolution of SNe II-P during the first month after explosion, as required in order to interpret optical observations of high-redshift SNe II-P, and to derive cross-filter K-corrections. While still highly preliminary, the apparent UV homogeneity of SNe II-P bodes well for the use of these events as cosmological probes at high redshift.

We use the radiative transfer code PHOENIX to study the line formation of the wavelength region 5000-7000 Angstroms. This is the region where the SNe Ia defining Si II feature occurs. This region is important since the ratio of the two nearby silicon lines has been shown to correlate with the absolute blue magnitude. We use a grid of LTE synthetic spectral models to investigate the formation of line features in the spectra of SNe Ia. By isolating the main contributors to the spectral formation we show that the ions that drive the spectral ratio are Fe III, Fe II, Si II, and S II. While the first two strongly dominate the flux transfer, the latter two form in the same physical region inside of the supernova. We also show that the naive blackbody that one would derive from a fit to the observed spectrum is far different than the true underlying continuum.

The ever-expanding depth and quality of photometric and spectroscopic observations of stellar populations increase the need for theoretical models in regions of age-composition parameter space that are largely unexplored at present. Stellar evolution models that employ the most advanced physics and cover a wide range of compositions are needed to extract the most information from current observations of both resolved and unresolved stellar populations. The Dartmouth Stellar Evolution Database is a collection of stellar evolution tracks and isochrones that spans a range of [Fe/H] from -2.5 to +0.5, [alpha/Fe] from -0.2 to +0.8 (for [Fe/H] <=0) or +0.2 (for [Fe/H] >0), and initial He mass fractions from Y=0.245 to 0.40. Stellar evolution tracks were computed for masses between 0.1 and 4 Msun, allowing isochrones to be generated for ages as young as 250 Myr. For the range in masses where the core He flash occurs, separate He-burning tracks were computed starting from the zero age horizontal branch. The tracks and isochrones have been transformed to the observational plane in a variety of photometric systems including standard UBV(RI)c, Stromgren uvby, SDSS ugriz, 2MASS JHKs, and HST ACS-WFC and WFPC2. The Dartmouth Stellar Evolution Database is accessible through a website at http://stellar.dartmouth.edu/~models/ where all tracks, isochrones, and additional files can be downloaded.

We present detailed radiative transfer spectral synthesis models for the Iron Low Ionization Broad Absorption Line (FeLoBAL) active galactic nuclei (AGN) FIRST J121442.3+280329 and ISO J005645.1-273816. Detailed NLTE spectral synthesis with a spherically symmetric outflow reproduces the observed spectra very well across a large wavelength range. While exact spherical symmetry is probably not required, our model fits are of high quality and thus very large covering fractions are strongly implied by our results. We constrain the kinetic energy and mass in the ejecta and discuss their implications on the accretion rate. Our results support the idea that FeLoBALs may be an evolutionary stage in the development of more ``ordinary'' QSOs.

A comparative study of optical spectra of Type Ia supernovae (SNe Ia) obtained near 1 week, 3 weeks, and 3 months after maximum light is presented. Most members of the four groups that were defined on the basis of maximum light spectra in Paper II (core normal, broad line, cool, and shallow silicon) develop highly homogeneous postmaximum spectra, although there are interesting exceptions. Comparisons with SYNOW synthetic spectra show that most of the spectral features can be accounted for in a plausible way. The fits show that 3 months after maximum light, when SN Ia spectra are often said to be in the nebular phase and to consist of forbidden emission lines, the spectra actually remain dominated by resonance scattering features of permitted lines, primarily those of Fe II. Even in SN 1991bg, which is said to have made a very early transition to the nebular phase, there is no need to appeal to forbidden lines at 3 weeks postmaximum, and at 3 months postmaximum the only clear identification of a forbidden line is [Ca II] 7291, 7324. Recent studies of SN Ia rates indicate that most of the SNe Ia that have ever occurred have been "prompt" SNe Ia, produced by young (100,000,000 yr) stellar populations, while most of the SNe Ia that occur at low redshift today are "tardy", produced by an older (several Gyrs) population. We suggest that the shallow silicon SNe Ia tend to be the prompt ones.

We calculate detailed NLTE synthetic spectra of a Pulsating Reverse Detonation (PRD) model, a novel explosion mechanism for Type Ia supernovae. While the hydro models are calculated in 3-D, the spectra use an angle averaged hydro model and thus some of the 3-D details are lost, but the overall average should be a good representation of the average observed spectra. We study the model at 3 epochs: maximum light, seven days prior to maximum light, and 5 days after maximum light. At maximum the defining Si II feature is prominent, but there is also a prominent C II feature, not usually observed in normal SNe Ia near maximum. We compare to the early spectrum of SN 2006D which did show a prominent C II feature, but the fit to the observations is not compelling. Finally we compare to the post-maximum UV+optical spectrum of SN 1992A. With the broad spectral coverage it is clear that the iron-peak elements on the outside of the model push too much flux to the red and thus the particular PRD realizations studied would be intrinsically far redder than observed SNe Ia. We briefly discuss variations that could improve future PRD models.

We present extensive ugrizYHJK photometry and optical spectroscopy of SN 2005gj obtained by the SDSS-II and CSP Supernova Projects, which give excellent coverage during the first 150 days after the time of explosion. These data show that SN 2005gj is the second clear case, after SN 2002ic, of a thermonuclear explosion in a dense circumstellar environment. Both the presence of singly and doubly ionized iron-peak elements (FeIII and weak SII, SiII) near maximum light as well as the spectral evolution show that SN 2002ic-like events are Type Ia explosions. Independent evidence comes from the exponential decay in luminosity of SN 2005gj, pointing to an exponential density distribution of the ejecta. The interaction of the supernova ejecta with the dense circumstellar medium is stronger than in SN 2002ic: (1) the supernova lines are weaker; (2) the Balmer emission lines are more luminous; and (3) the bolometric luminosity is higher close to maximum light. The velocity evolution of the Halpha components suggest that the CSM around SN 2005gj is clumpy and it has a flatter density distribution compared with the steady wind solution, in agreement with SN 2002ic. An early X-ray observation with Chandra gives an upper-limit on the mass loss rate from the companion of < 2x10^{-4} Msun/yr.

The ACS Survey of Galactic Globular Clusters, an HST Treasury Project, will deliver high quality, homogeneous photometry of 65 globular clusters. This paper introduces a new collection of stellar evolution tracks and isochrones suitable for analyzing the ACS Survey data. Stellar evolution models were computed at [Fe/H]= -2.5, -2.0, -1.5, -1.0, -0.5, and 0; [alpha/Fe]= -0.2, 0, 0.2, 0.4, 0.6, and 0.8; and three initial He abundances for masses from 0.1 to 1.8 Msun and ages from 2 to 15 Gyr. Each isochrone spans a wide range in luminosity from Mv~14 up to the tip of the red giant branch. These are complemented by a set of He-burning tracks that extend from the zero age horizontal branch to the onset of thermal pulsations on the asymptotic giant branch. In addition, a set of computer programs are provided that make it possible to interpolate the isochrones in [Fe/H], generate luminosity functions from the isochrones, and create synthetic horizontal branch models. The tracks and isochrones have been converted to the observational plane with two different color-Teff transformations, one synthetic and one semi-empirical, in ground-based B, V, and I, and F606W and F814W for both ACS-WFC and WFPC2 systems. All models and programs presented in this paper are available from http://stellar.dartmouth.edu/~models/

Integrated light from distant galaxies is often compared to stellar population models via the equivalent widths of spectral features--spectral indices--whose strengths rely on the abundances of one or more elements. Such comparisons hinge not only on the overall metal abundance but also on relative abundances. Studies have examined the influence of individual elements on synthetic spectra but little has been done to address similar issues in the stellar evolution models that underlie most stellar population models. Stellar evolution models will primarily be influenced by changes in opacities. In order to explore this issue in detail, twelve sets of stellar evolution tracks and isochrones have been created at constant heavy element mass fraction Z that self-consistently account for varying heavy element mixtures. These sets include scaled-solar, alpha-enhanced, and individual cases where the elements C, N, O, Ne, Mg, Si, S, Ca, Ti, and Fe have been enhanced above their scaled-solar values. The variations that arise between scaled-solar and the other cases are examined with respect to the H-R diagram and main sequence lifetimes.

A comparative study of spectra of 21 Type Ia supernovae (SNe Ia) obtained about one week before maximum light, and 8 spectra obtained 11 or more days before maximum, is presented. To a large extent the premaximum spectra exhibit the defining characteristics of the four groups defined in Paper II (core-normal, broad-line, cool, and shallow-silicon). Comparisons with SYNOW synthetic spectra show that all strong features and most weak ones can be accounted for in a plausible way. The issues of detached high-velocity features, the possible ubiquity of carbon clumps, the maximum detectable ejecta velocities, and the possibility of blueshifted emission-line peaks are discussed.

We present spectral fits to five epochs of the typical Type Ib supernova 1999dn using the generalized, non-LTE, stellar atmospheres code PHOENIX. Our goal is threefold: to determine basic physical properties of the supernova ejecta, such as velocity, temperature, and density gradients; to reproduce He I absorption lines by invoking non-thermal excitation; and, to investigate possible spectral signatures of hydrogen, especially a feature around 6200 Angstrom, which has been attributed to high velocity $H_\alpha$. Our models assume an atmosphere with uniform composition devoid of any hydrogen. Our model spectra fit the observed spectra well, successfully reproducing most of the features, including the prominent He I absorptions. The most plausible alternative to $H_\alpha$ as the source of the 6200 Angstrom feature is a blend of Fe II and Si II lines, which can be made stronger to fit the observed feature better by increasing the metallicity of the ejecta. High-metallicity models fit well at early epochs, but not as well as solar-metallicity models after maximum light. While this blend of metal lines is a reasonable explanation of the source of the 6200 Angstrom feature, it is still important to investigate hydrogen as the source; therefore, a second paper will present models that include a thin shell of hydrogen around the main composition structure.

Two quantitative tests DIFF1 and DIFF2 for measuring goodness-of-fit between two locally-normalized supernova spectra are presented. Locally-normalized spectra are obtained by dividing a spectrum by the same spectrum smoothed over a wavelength interval relatively large compared to line features, but relatively small compared to continuum features. DIFF1 essentially measures the mean relative difference between the line patterns of locally-normalized spectra and DIFF2 is DIFF1 minimized by a relative logarithmic wavelength shift between the spectra: the shift is effectively an artificial relative Doppler shift. Both DIFF1 and DIFF2 measure the physical similarity of line formation, and thus of supernovae. DIFF1 puts more weight on overall physical similarity of the supernovae than DIFF2 because the DIFF2 shift compensates somewhat for some physical distinction in the supernovae. Both tests are useful in ordering supernovae into empirical groupings for further analysis. We present some examples of locally-normalized spectra for Type IIb supernova SN 1993J with some analysis of these spectra. The UV parts of two of the SN 1993J spectra are HST spectra that have not been published before. We also give an example of fitted locally-normalized spectra and, as an example of the utility of DIFF1 and DIFF2, some preliminary statistical results for hydrogen-deficient core-collapse (HDCC) supernova spectra. This paper makes use of and refers to material to found at the first author's online supernova spectrum database SUSPEND (SUpernovae Spectra PENDing further analysis: see http://www.nhn.ou.edu/~jeffery/astro/sne/spectra/spectra.html)

Higher resolution telescopes as well as 3D numerical simulations will require the development of detailed 3D radiative transfer calculations. Building upon our previous work we extend our method to include both continuum and line transfer. We present a general method to calculate radiative transfer including scattering in the continuum as well as in lines in 3D static atmospheres. The scattering problem for line transfer is solved via means of an operator splitting (OS) technique. The formal solution is based on a long-characteristics method. The approximate $\Lambda$ operator is constructed considering nearest neighbors {\em exactly}. The code is parallelized over both wavelength and solid angle using the MPI library. We present the results of several test cases with different values of the thermalization parameter and two choices for the temperature structure. The results are directly compared to 1D spherical tests. With our current grid setup the interior resolution is much lower in 3D than in 1D, nevertheless the 3D results agree very well with the well-tested 1D calculations. We show that with relatively simple parallelization that the code scales to very large number of processors which is mandatory for practical applications. Advances in modern computers will make realistic 3D radiative transfer calculations possible in the near future. Our current code scales to very large numbers of processors, but requires larger memory per processor at high spatial resolution.

Type Ia SNe (SNe Ia) are thought to come from carbon-oxygen white dwarfs that accrete mass from binary companions until they approach the Chandrasekhar limit, ignite carbon, and undergo complete thermonuclear disruption. A survey of the observed types of binaries that contain white dwarfs is presented. We propose that certain systems that seem most promising as SN Ia progenitors should be more intensively observed and modeled, to determine whether the white dwarfs in these systems will be able to reach the Chandrasekhar limit. In view of the number of promising single-degenerate systems and the dearth of promising double-degenerate systems, we suspect that single-degenerates produce most or perhaps all SNe Ia, while double-degenerates produce some or perhaps none.

We present detailed NLTE spectral synthesis models of the Type II supernova 2005cs, which occurred in M51 and for which the explosion time is well determined. We show that previous estimates for the reddening were significantly too high and briefly discuss how this will effect the inferred progenitor mass. We also show that standard CNO-burning enhanced abundances require far too large an oxygen depletion, although there is evidence for a single optical N II line and the sodium abundance shows clear evidence for enhancement over solar both as expected from CNO processing. Finally we calculate a distance using the SEAM method. Given the broad range of distances to M51 in the literature, the determination of a distance using Cepheid variables would be quite valuable.

We derive the radiative transfer equation for arbitrary stationary relativistic flows in stationary spacetimes, i.e. for steady-state transfer problems. We show how the standard characteristics method of solution developed by Mihalas and used throughout the radiative transfer community can be adapted to multi-dimensional applications with isotropic sources. Because the characteristics always coincide with geodesics and can always be specified by constants, direct integration of the characteristics derived from the transfer equation as commonly done in 1-D applications is not required. The characteristics are known for a specified metric from the geodesics. We give details in both flat and static spherically symmetric spacetimes. This work has direct application in 3-dimensional simulations of supernovae, gamma-ray bursts, and active galactic nuclei, as well as in modeling neutron star atmospheres.

We are developing a brand new stellar population models with flexible chemistry (isochrones plus stellar colors and spectra) in order to set a new standard of completeness and excellence. Here we present preliminary results to assess the effects of stellar evolution models and stellar model atmosphere to the well-known Lick indices at constant heavy element mass fraction Z that self-consistently account for varying heavy element mixtures. We have enhanced chemical elements one by one. Our ultimate goal is to demonstrate 10% absolute mean ages for a sample of local galaxies derived from an integrated light spectrum.

Synthetic spectra generated with the parameterized supernova synthetic-spectrum code SYNOW are compared to spectra of the unusual Type Ib supernova 2005bf. We confirm the discovery by Folatelli et al. (2006) that very early spectra (about 30 days before maximum light) contain both photospheric-velocity (8000 km/s) features of He I, Ca II, and Fe II, and detached high-velocity (14,000 km/s) features of H-alpha, Ca II, and Fe II. An early spectrum of SN 2005bf is an almost perfect match to a near-maximum-light spectrum of the Type Ib SN 1999ex. Although these two spectra were at very different times with respect to maximum light (20 days before maximum for SN 2005bf and five days after for SN 1999ex), they were for similar times after explosion - about 20 days for SN 2005bf and 24 days for SN 1999ex. The almost perfect match clinches the previously suggested identification of H-alpha in SN 1999ex and supports the proposition that many if not all Type Ib supernovae eject a small amount of hydrogen. The earliest available spectrum of SN 2005bf resembles a near-maximum-light spectrum of the Type Ic SN 1994I. These two spectra also were at different times with respect to maximum light (32 days before maximum for SN 2005bf and four days before for SN 1994I) but at similar times after explosion - about eight days for SN 2005bf and 10 days for SN 1994I. The resemblance motivates us to consider a reinterpretation of the spectra of Type Ic supernovae, involving coexisting photospheric-velocity and high-velocity features. The implications of our results for the geometry of the SN 2005bf ejecta, which has been suggested to be grossly asymmetric, are briefly discussed.

We present a general method to calculate radiative transfer including scattering in the continuum as well as in lines in spherically symmetric systems that are influenced by the effects of general relativity (GR). We utilize a comoving wavelength ansatz that allows to resolve spectral lines throughout the atmosphere. The used numerical solution is an operator splitting (OS) technique that uses a characteristic formal solution. The bending of photon paths and the wavelength shifts due to the effects of GR are fully taken into account, as is the treatment of image generation in a curved spacetime. We describe the algorithm we use and demonstrate the effects of GR on the radiative transport of a two level atom line in a neutron star like atmosphere for various combinations of continuous and line scattering coefficients. In addition, we present grey continuum models and discuss the effects of different scattering albedos on the emergent spectra and the determination of effective temperatures and radii of neutron star atmospheres.

Using models of the SN~IIP 2005cs, we show that detailed spectral analysis can be used to determine reddening and abundances.

We present four spectra of the Type Ia supernova (SN Ia) 2006D extending from -7 to +13 days with respect to B-band maximum. The spectra include the strongest signature of unburned material at photospheric velocities observed in a SN Ia to date. The earliest spectrum exhibits C II absorption features below 14,000 km/s, including a distinctive C II \lambda 6580 absorption feature. The carbon signatures dissipate as the SN approaches peak brightness. In addition to discussing implications of photospheric-velocity carbon for white dwarf explosion models, we outline some factors that may influence the frequency of its detection before and around peak brightness. Two effects are explored in this regard, including depopulation of the C II optical levels by non-LTE effects, and line-of-sight effects resulting from a clumpy distribution of unburned material with low volume-filling factor.

A previously unpublished ultralow-dispersion spectrum of Sanduleak -69 202, the stellar progenitor of SN 1987A, is presented and the uncertain presupernova evolution of Sanduleak -69 202 is discussed.

Howell et al. have reported the discovery of SN Ia SN 2003fg (SNLS-03D3bb) and conclude that SN 2003fg is very likely a super-Chandrasekhar-mass SN Ia perhaps with a mass of order 2 solar masses. Their work is the first strong evidence that has been presented for a super-Chandrasekhar SN Ia. We have performed an analysis of the SN 2003fg data using the Yoon & Langer binding energy formula for a rotating super-Chandrasekhar-mass white dwarf (also used by Howell et al.) included in a simple model of SNe Ia (which we call the SSC model for Simple Super-Chandrasekhar model for SNe Ia) which assumes spherically symmetric ejecta and relies on the approximations of an exponential density profile for SN Ia ejecta and of a sharp boundary of the SN Ia iron-peak-element core. Our results support the conclusion of Howell et al.: SN 2003fg is very probably super-Chandrasekhar and probably has mass of order 2 solar masses.

JEDI (Joint Efficient Dark-energy Investigation) is a candidate implementation of the NASA-DOE Joint Dark Energy Mission (JDEM). JEDI will probe dark energy in three independent methods: (1) type Ia supernovae, (2) baryon acoustic oscillations, and (3) weak gravitational lensing. In an accompanying paper, an overall summary of the JEDI mission is given. In this paper, we present further details of the supernova component of JEDI. To derive model-independent constraints on dark energy, it is important to precisely measure the cosmic expansion history, H(z), in continuous redshift bins from z \~ 0-2 (the redshift range in which dark energy is important). SNe Ia at z > 1 are not readily accessible from the ground because the bulk of their light has shifted into the near-infrared where the sky background is overwhelming; hence a space mission is required to probe dark energy using SNe. Because of its unique near-infrared wavelength coverage (0.8-4.2 microns), JEDI has the advantage of observing SNe Ia in the rest frame J band for the entire redshift range of 0 < z < 2, where they are less affected by dust, and appear to be nearly perfect standard candles. During the first year of JEDI operations, spectra and light curves will be obtained for ~4,000 SNe Ia at z < 2. The resulting constraints on dark energy are discussed, with special emphasis on the improved precision afforded by the rest frame near-infrared data.

By definition, a Type Ic supernova (SN Ic) does not have conspicuous lines of hydrogen or helium in its optical spectrum. SNe Ic usually are modelled in terms of the gravitational collapse of bare carbon-oxygen cores. We consider the possibility that the spectra of ordinary (SN 1994I-like) SNe Ic have been misinterpreted, and that SNe Ic eject hydrogen. An absorption feature usually attributed to a blend of Si II 6355 and C II 6580 may be produced by H-alpha. If SN 1994I-like SNe Ic eject hydrogen, the possibility that hypernova (SN 1998bw-like) SNe Ic, some of which are associated with gamma-ray bursts, also eject hydrogen should be considered. The implications of hydrogen for SN Ic progenitors and explosion models are briefly discussed.

The spectroscopic properties of a selected optical photospheric spectra of core collapse supernovae (CCSNe) are investigated.Special attention is devoted to traces of hydrogen at early phases. The generated spectra are found to match the observed ones reasonably well, including a list of only 23 candidate ions. Guided by SN Ib 1990I, the observed trough near 6300\AA is attributed to H$\alpha$ in almost all Type Ib events, although in some objects it becomes too weak to be discernible, especially at later phases. Alternative line identifications are discussed. Differences in the way hydrogen manifests its presence within CCSNe are highlighted. In Type Ib SNe, the H$\alpha$ contrast velocity (i.e. line velocity minus the photospheric velocity) seems to increase with time at early epochs, reaching values as high as 8000 km s$^{-1}$ around 15-20 days after maximum and then remains almost constant. The derived photospheric velocities, indicate a lower velocity for Type II SNe 1987A and 1999em as compared to SN Ic 1994I and SN IIb 1993J, while Type Ib events display a somewhat larger variation. The scatter, around day 20, is measured to be $\sim$5000 km s$^{-1}$. Following two simple approaches, rough estimates of ejecta and hydrogen masses are given. A mass of hydrogen of approximately 0.02 $M_\odot$ is obtained for SN 1990I, while SNe 1983N and 2000H ejected $\sim$0.008 $M_\odot$ and $\sim$0.08 $M_\odot$ of hydrogen, respectively. SN 1993J has a higher hydrogen mass, $\sim 0.7$ $M_\odot$ with a large uncertainty. A low mass and thin hydrogen layer with very high ejection velocities above the helium shell, is thus the most likely scenario for Type Ib SNe. Some interesting and curious issues relating to oxygen lines suggest future investigations.

We have performed detailed NLTE spectral synthesis modeling of 2 types of 1-D hydro models: the very highly parameterized deflagration model W7, and two delayed detonation models. We find that overall both models do about equally well at fitting well observed SNe Ia near to maximum light. However, the Si II 6150 feature of W7 is systematically too fast, whereas for the delayed detonation models it is also somewhat too fast, but significantly better than that of W7. We find that a parameterized mixed model does the best job of reproducing the Si II 6150 line near maximum light and we study the differences in the models that lead to better fits to normal SNe Ia. We discuss what is required of a hydro model to fit the spectra of observed SNe Ia near maximum light.

A short biography of Hans Bethe.

We describe a highly flexible framework to solve 3D radiation transfer problems in scattering dominated environments based on a long characteristics piece-wise parabolic formal solution and an operator splitting method. We find that the linear systems are efficiently solved with iterative solvers such as Gauss-Seidel and Jordan techniques. We use a sphere-in-a-box test model to compare the 3D results to 1D solutions in order to assess the accuracy of the method. We have implemented the method for static media, however, it can be used to solve problems in the Eulerian-frame for media with low velocity fields.

The absolute visual magnitudes of three Type IIb, 11 Type Ib and 13 Type Ic supernovae (collectively known as stripped-envelope supernovae) are studied by collecting data on the apparent magnitude, distance, and interstellar extinction of each event. Weighted and unweighted mean absolute magnitudes of the combined sample as well as various subsets of the sample are reported. The limited sample size and the considerable uncertainties, especially those associated with extinction in the host galaxies, prevent firm conclusions regarding differences between the absolute magnitudes of supernovae of Type Ib and Ic, and regarding the existence of separate groups of overluminous and normal-luminosity stripped-envelope supernovae. The spectroscopic characteristics of the events of the sample are considered. Three of the four overluminous events are known to have had unusual spectra. Most but not all of the normal luminosity events had typical spectra. Light curves of stripped-envelope supernovae are collected and compared. Because SN 1994I in M51 was very well observed it often is regarded as the prototypical Type Ic supernova, but it has the fastest light curve in the sample. Light curves are modeled by means of a simple analytical technique that, combined with a constraint on E/M from spectroscopy, yields internally consistent values of ejected mass, kinetic energy, and nickel mass.

A comparative study of near-maximum-light optical spectra of 24 Type Ia supernovae (SNe Ia) is presented. The spectra are quantified in two ways, and assigned to four groups. Seven "core-normal" SNe Ia have very similar spectra, except for strong high-velocity CaII absorption in SN 2001el. Seven SNe Ia are assigned to a "broad-line" group, the most extreme of which is SN 1984A. Five SNe Ia, including SN 1991bg, are assigned to a "cool" group. Five SNe Ia, including SN 1991T, are assigned to a "shallow-silicon" group. Comparisons with SYNOW synthetic spectra provide a basis for discussion of line identifications, and an internally consistent quantification of the maximum-light spectroscopic diversity among SNe Ia. The extent to which SN Ia maximum-light spectra appear to have a continuous distribution of properties, rather than consisting of discrete subtypes, is discussed.

Type Ia supernovae have played a crucial role in the discovery of the dark energy, via the measurement of their light curves and the determination of the peak brightness via fitting templates to the observed lightcurve shape. Two spectroscopic indicators are also known to be well correlated with peak luminosity. Since the spectroscopic luminosity indicators are obtained directly from observed spectra, they will have different systematic errors than do measurements using photometry. Additionally, these spectroscopic indicators may be useful for studies of effects of evolution or age of the SNe Ia progenitor population. We present several new variants of such spectroscopic indicators which are easy to automate and which minimize the effects of noise. We show that these spectroscopic indicators can be measured by proposed JDEM missions such as SNAP and JEDI.

We present early time high-resolution (VLT/UVES) and late time low-resolution (VLT/FORS) optical spectra of the normal type Ia supernova, SN 2001el. The high-resolution spectra were obtained at -9 and -2 days to allow the detection of narrow hydrogen and/or helium emission lines from the circumstellar medium of the SN. No such lines were detected, and we therefore use photoionisation models to derive upper limits of 9x10^-6 Msun/yr and 5x10^-5 Msun/yr for the mass loss rate from the progenitor system assuming velocities of 10 km/s and 50 km/s, respectively, for a wind extending to outside at least a few x 10^15 cm away from the SN explosion site. These limits exclude a symbiotic star in the upper mass loss rate regime from being the progenitor of SN 2001el. The low resolution spectrum was obtained in the nebular phase of the SN, 400 days after the maximum light, to search for any hydrogen rich gas originating from the SN progenitor system. However, we see no signs of Balmer lines in our spectrum. Therefore, we model the late time spectra to derive an upper limit of ~0.03 Msun for solar abundance material present at velocities lower than 1000 km/s within the SN explosion site. According to simulations of Marietta et al. (2000) this is less than the expected mass lost by a subgiant, red giant or main sequence secondary star at a small binary separation as a result of the SN explosion. Finally, we discuss the origin of high velocity Ca II lines. We see both the CaII IR triplet and the H&K lines in the -9 days spectrum at a very high velocity of up to 34000 km/s. The spectrum also shows a flat-bottomed Si II `6150 A' feature similar to the one previously observed in SN 1990N at -14 days. We compare these spectral features to those observed in SNe 1984A and 1990N at even higher velocities.

Detailed quantitative spectroscopy of Type Ia supernovae (SNe~Ia) provides crucial information needed to minimize systematic effects in both ongoing SNe Ia observational programs such as the Nearby Supernova Factory, ESSENCE, and the SuperNova Legacy Survey (SNLS) and in proposed JDEM missions such as SNAP, JEDI, and DESTINY. Quantitative spectroscopy is mandatory to quantify and understand the observational strategy of comparing ``like versus like''. It allows us to explore evolutionary effects, from variations in progenitor metallicity to variations in progenitor age, to variations in dust with cosmological epoch. It also allows us to interpret and quantify the effects of asphericity, as well as different amounts of mixing in the thermonuclear explosion.

We present high resolution VLT UVES observations of the active K dwarf GJ 117. 6^Li enhancement has been shown for energetic solar events, one chromospherically active binary, and several dwarf halo stars. Our analysis reveals the detection of 6^Li on this source with 6^Li/7^Li = 0.030+/-0.010. We found no significant contribution from other lines, including Ti I, in the Li profile of GJ 117 and a template star of similar spectral type and metallicity. We discuss the possibility for 6^Li production by spallation and find it to be consistent with the activity levels of the object.

We present analysis from simultaneous FUSE, ASCA, and EUVE observations, as well as a reanalysis of archival HST spectra, from the extreme Narrow-line Seyfert 1 Galaxy RE 1034+39 (KUG 1031+398). RE 1034+39 has an unusually hard spectral energy distribution (SED) that peaks in the soft X-rays. Its emission lines are unusual in that they can all be modelled as a Lorentzian centered at the rest wavelength with only a small range in velocity widths. In order to investigate whether the unusual SED influences the emission line ratios and equivalent widths, we present three complementary types of photoionization analysis. The FUSE spectrum was particularly important because it includes the high-ionization line OVI. First, we use the photoionization code Cloudy and the SED developed from the coordinated observations to confirm that the emission lines are consistent with observed hard SED. The best model parameters were an ionization parameter log(U) ~ -2 and a hydrogen number density log(n_H)=9.75 [cm^-2]. Second, we present a Locally Optimally-emitting Cloud model. This model produced enhanced OVI as observed, but also yielded far too strong MgII. Third, we develop a series of semi-empirical SEDs, run Cloudy models, and compare the results with the measured values using a figure of merit (FOM). The FOM minimum indicates similar SED and gas properties as were inferred from the one-zone model using the RE 1034+39 continuum. Furthermore, the FOM increases sharply toward softer continua, indicating that a hard SED is required by the data in the context of a one-zone model.

JEDI will probe dark energy in three independent ways by measuring the expansion history of the universe: (1) using type Ia supernovae as cosmological standard candles over a range of distances, (2) using baryon oscillations as a cosmological standard ruler over a range of cosmic epochs, (3) mapping the weak gravitational lensing distortion by foreground galaxies of the images of background galaxies at different distances. JEDI will unravel the nature of dark energy with accuracy and precision. JEDI is a 2m-class space telescope with the unique ability of simultaneous wide-field imaging (0.8-4.2 micron in five bands) and multi-slit spectroscopy (0.8-3.2 micron) with a field of view of 1 square degree. What makes JEDI efficient is its ability to simultaneously obtain high signal-to-noise ratio, moderate resolution slit spectra for all supernovae and ~ 5000 galaxies in its wide field of view, and to combine imaging and spectroscopy so that the appropriate balance of time is devoted to each. JEDI will measure the cosmic expansion history H(z) as a free function to < 2% accuracy in redshift bins of 0.2-0.3. Assuming a flat universe and \sigma(\Omega_m)=0.01 (0.03), JEDI could measure a departure from a vanilla \Lambda CDM model (w_0=-1, w'=0) with \sigma(w_0)=0.013 (0.031) and \sigma(w')=0.046 (0.063). JEDI will obtain the well-sampled lightcurves in Z, J, H, K, L bands and spectra of ~ 14,000 type Ia supernovae with redshifts ranging from 0 to 1.7; the redshifts of ~ 10-100 million galaxies to H ~ 23 and z ~ 4 over 1000-10,000 square degrees; and measurements of the shapes of galaxies over 1000-10,000 square degrees in Z,J,H,K,L for > 10^9 galaxies to H \~ 25.

We present spectroscopic and photometric observations of the peculiar Type II supernova (SN) 1998A. The light curves and spectra closely resemble those of SN 1987A, suggesting that the SN 1998A progenitor exploded when it was a compact blue supergiant. However, the comparison with SN 1987A also highlights some important differences: SN 1998A is more luminous and the spectra show bluer continua and larger expansion velocities at all epochs. These observational properties indicate that the explosion of SN 1998A is more energetic than SN 1987A and more typical of SNe II. Comparing the observational data to simulations, we deduce that the progenitor of SN 1998A was a massive star (~ 25 Mo) with a small pre-supernova radius (< 6 x 10^{12} cm). The Ba II lines, unusually strong in SN 1987A and some faint II--P events, are almost normal in the case of SN 1998A, indicating that the temperature plays a key role in determining their strength.

As the first step in a comprehensive, comparative, direct analysis of the spectra of Type Ia supernovae (SNe Ia), we use the parameterized supernova synthetic-spectrum code, SYNOW, to interpret 26 spectra of the well-observed SN 1994D. Our results are consistent with the traditional view that the composition structure (element abundance fractions versus ejection velocity) is radially stratified. We find that resonance-scattering features due to permitted lines of Ca II, Na I, and Fe II persist to more than 100 days after explosion. The fitting parameters for SN 1994D, together with those to be determined for other SNe Ia, will provide an internally consistent quantification of the spectroscopic diversity among SNe Ia, and shed light on how the various manifestations of observational diversity are related to their physical causes.

Due to their intrinsic brightness, supernovae make excellent cosmological probes. We describe the SEAM method for obtaining distances to Type IIP supernovae (SNe IIP) and present a distance to SN 1999em for which a Cepheid distance exists. Our models give results consistent with the Cepheid distance, even though we have not attempted to tune the underlying hydrodynamical model, we have simply chosen the best fits. This is in contradistinction to the expanding photosphere method (EPM) which yields a distance to SN 1999em that is 50% smaller than the Cepheid distance. We emphasize the differences between SEAM and EPM. We show that the dilution factors used in the EPM analysis were systematically too small at later epochs. We also show that the EPM blackbody assumption is suspect. Since SNe IIP are visible to redshifts as high as z less than about 6, with the JWST, SEAM may be a valuable probe of the early universe.

We briefly describe the current version of the PHOENIX code. We then present some results on the modeling of Type II supernovae and show that fits to observations can be obtained, when account is taken for spherically symmetric, line-blanketed, expanding atmospheres. We describe the SEAM method of obtaining distances to supernovae and briefly discuss its future prospects.

Ground-based and Hubble Space Telescope observations are presented for SN 1993J and SN 1998S. SN 1998S shows strong, relatively narrow circumstellar emission lines of N III-V and C III-IV, as well as broad lines from the ejecta. Both the broad ultraviolet and optical lines in SN 1998S indicate an expansion velocity of ~7,000 km/s. The broad emission components of Ly-alpha and Mg II are strongly asymmetrical after day 72 past the explosion, and differ in shape from H-alpha. Different models based on dust extinction from dust in the ejecta or shock region, in combination with H-alpha from a circumstellar torus, are discussed. It is concluded, however, that the double-peaked line profiles are more likely to arise as a result of optical depth effects in the narrow, cool, dense shell behind the reverse shock, than in a torus-like region. The ultraviolet lines of SN 1993J are broad, with a box-like shape, coming from the ejecta and a cool dense shell. The shapes of the lines are well fitted with a shell with inner velocity ~7,000 km/s and outer velocity ~10,000 km/s. For both SN 1993J and SN 1998S a strong nitrogen enrichment is found, with N/C~12.4 in SN 1993J and N/C~6.0 in SN 1998S. From a compilation of all supernovae with determined CNO ratios, we discuss the implications of these observations for the structure of the progenitors of Type II supernovae.

Recently, with the advances in computational speed and availability there has been a growth in the number and resolution of fully 3-D hydrodynamical simulations. However, all of these simulations are purely hydrodynamical and there has been little attempt to include the effects of radiative transfer except in a purely phenomenological manner because the computational cost is too large even for modern supercomputers. While there has been an effort to develop 3-D Monte Carlo radiative transfer codes, most of these have been for static atmospheres or have employed the Sobolev approximation, which limits their applicability to studying purely geometric effects such as macroscopic mixing. Also the computational requirements of Monte Carlo methods are such that it is difficult to couple with 3-D hydrodynamics. Here, we present an algorithm for calculating 1-D spherical radiative transfer in the presence of non-monotonic velocity fields in the co-moving frame. Non-monotonic velocity flows will occur in convective, and Raleigh--Taylor unstable flows, in flows with multiple shocks, and in pulsationally unstable stars such as Mira and Cepheids. This is a first step to developing fully 3-D radiative transfer than can be coupled with hydrodynamics. We present the computational method and the results of some test calculations.

In spite of the apparent lack of Si II and S II features in its spectra, SN 2002cx was classified as a peculiar Type Ia supernova (SN Ia) on the basis of its overall photometric and spectroscopic behavior. Spectra obtained near maximum light contained Fe III features, as in SN 1991T-like events, but the blueshifts of the Fe III absorptions were exceptionally low. The luminosity also was low. We use the supernova synthetic--spectrum code SYNOW to study line identifications in SN 2002cx. We find that the maximum-light spectra appear to contain weak features of Si II, S II, Si III, and Ca II, which strengthens the connection with SN 1991T-like events. We show that later spectra, obtained 12, 25, and 56 days after maximum, consist of P-Cygni resonance-scattering features due to permitted Fe II and Co II lines. SN 2002cx had been thought to have made the transition from a permitted-line to a forbidden-line spectrum between 25 and 56 days. Owing to the low expansion velocities the postmaximum spectral features are narrower and easier to identify than they are in other SNe Ia. SN 2002cx will lead to improved line identifications in other SNe Ia and clarify when the transition from a permitted-line to a forbidden-line spectrum occurs. In the context of current SN Ia explosion models, we suggest that the properties of SN 2002cx may be consistent with 3D deflagration models, which are not favored for normal SNe Ia.

Brown dwarfs bridge the gap between the stellar and planetary mass regimes. Evolving from conditions very similar to the lowest-mass stars, the atmospheres of older brown dwarfs closely resemble those expected in close-in extrasolar giant planets, and with cooler BDs still being discovered, more and more approach the properties of gas giants at wider separation. Interpreting the spectra of BDs is therefore a crucial step towards understanding and predicting the spectral and thermal properties of EGPs. Essential properties of substellar atmospheres are massive molecular line-blanketing and the condensation of species with decreasing Teff, changing the chemical equilibrium composition and causing absorption from dust grains. More complex details involve the distribution of dust clouds over the surface giving rise to temporal variability, and possible deviations from chemical equilibrium conditions. In the case of close-in EGPs and some BDs in binary systems, the effect of irradiation from the primary significantly affects the spectral properties and thermal evolution.

We find that spectra of certain Iron Low Ionization Broad Absorption Line (FeLOBAL) QSOs, which are characterized by low--ionization emission and blue shifted absorption lines, can be well matched with the spectral synthesis code SYNOW. SYNOW is a resonance scattering code and assumes that line emission comes from a single line forming region. This interpretation is novel as traditionally line emission and absorption in BALQSOs are thought to come from two different regions. We extend this analysis by using the detailed PHOENIX code to model the spectra. We present a SYNOW fit and a preliminary model result from PHOENIX.

The Type Ia SN 2000cx exhibited multiple peculiarities, including a lopsided B-band light-curve peak that does not conform to current methods for using shapes of light curves to standardize SN Ia luminosities. We use the parameterized supernova synthetic-spectrum code SYNOW to study line identifications in the photospheric-phase spectra of SN 2000cx. Previous work established the presence of Ca II infrared-triplet features forming above velocity about 20,000 km/s, much higher than the photospheric velocity of about 10,000 km/s. We find Ti II features forming at the same high velocity. High-velocity line formation is partly responsible for the photometric peculiarities of SN 2000cx: for example, B-band flux blocking by Ti II absorption features that decreases with time causes the B light curve to rise more rapidly and decline more slowly than it otherwise would. SN 2000cx contains an absorption feature near 4530 A that may be H-beta, forming at the same high velocity. The lack of conspicuous H-alpha and P-alpha signatures does not necessarily invalidate the H-beta identification if the high-velocity line formation is confined to a clump that partly covers the photosphere and the H-alpha and P-alpha source functions are elevated relative to that of resonance scattering. The H-beta identification is tentative. If it is correct, the high-velocity matter must have come from a nondegenerate companion star.

We describe two separate wavelength discretization schemes that can be used in the numerical solution of the comoving frame radiative transfer equation. We present an improved second order discretization scheme and show that it leads to significantly less numerical diffusion than previous scheme. We also show that due to the nature of the second order term in some extreme cases it can become numerically unstable. We stabilize the scheme by introducing a mixed discretization scheme and present the results from several test calculations.

High-velocity features in Type Ia supernova spectra provide a way to probe the outer layers of these explosions. The maximum-light spectra of the unique Type Ia supernova 2000cx exhibit interesting Ca II features with high-velocity components. The infrared triplet absorption is quadruply notched, while the H&K absorption is wide and flat. Stimulated by a three-dimensional interpretation of similar Ca II features in another Type Ia supernova (SN 2001el, Kasen et al. 2003), we present alternative spherically symmetric and three-dimensional ejecta models to fit the high-velocity (v > 16,000 km/s) Ca II features of SN 2000cx. We also present simple estimates of the high-velocity ejecta mass for a few trial compositions and discuss their implications for explosion modelling.

In this paper we present spectroscopic and photometric observations for four core collapse supernovae (SNe), namely SNe 1994N, 1999br, 1999eu and 2001dc. Together with SN 1997D, we show that they form a group of exceptionally low-luminosity events. These SNe have narrow spectral lines (indicating low expansion velocities) and low luminosities at every phase (significantly lower than those of typical core-collapse supernovae). The very low luminosity during the ^{56}Co radioactive decay tail indicates that the mass of ^{56}Ni ejected during the explosion is much smaller (M(Ni) ~ 2-8 x 10^{-3} Mo) than the average (M(Ni) ~ 6-10 x 10^{-2} Mo). Two supernovae of this group (SN 1999br and SN 2001dc) were discovered very close to the explosion epoch, allowing us to determine the lengths of their plateaux (~ 100 days) as well as establishiing the explosion epochs of the other, less-completely observed SNe. It is likely that this group of SNe represent the extreme low-luminosity tail of a single continuous distribution of SN II-P events. Their kinetic energy is also exceptionally low. Although an origin from low mass progenitors has also been proposed for low-luminosity core-collapse SNe, recent work provides evidence in favour of the high mass progenitor scenario. The incidence of these low--luminosity SNe could be as high as 4-5$ $ of all type II SNe.

Since April 2000, a program using the ESO/VLT/UVES to search for early circumstellar signatures from nearby supernovae has been conducted. Until now, two Type Ia supernovae have been observed, SNe 2000cx and 2001el. Here we report on preliminary results for SNe 2001el, and we discuss how the observations can be used together with detailed modeling to derive an upper limit on the putative wind from the progenitor system. For a hydrogen-rich wind with velocity 10 km/s, the mass loss rate for the progenitor system of SN 2001el is less than 1O^{-5} solar masses per year.

We present 29 optical spectra of the normal Type Ia SN 1998aq, ranging from 9 days before to 241 days after the time of maximum brightness. This spectroscopic data set, together with photometric data presented elsewhere, makes SN 1998aq one of the best observed Type Ia supernova at optical wavelengths. We use the parameterized supernova synthetic-spectrum code SYNOW to study line identifications in the early photospheric-phase spectra. The results include evidence for lines of singly ionized carbon, at ejection velocities as low as 11,000 km/sec. Implications for explosion models are discussed.

Motivated by recent 3-D calculations of the explosion of Type Ia supernova via a pure deflagration we calculate the observed spectra at 15-25 days past maximum light of a parameterized model which has a considerable fraction of unburnt C+O in the central regions. Rather than attempting a self consistent 3-D calculation, which is beyond the scope of current computer codes, we modify the composition structure of the 1-D deflagration model W7. In our exploratory parameterized calculations, we find that a central concentration of C+O is not ruled out by observations for the epochs we study. We briefly examine whether nebular phase spectra could be incompatible with observations.

SN 2001el is the first normal Type Ia supernova to show a strong, intrinsic polarization signal. In addition, during the epochs prior to maximum light, the CaII IR triplet absorption is seen distinctly and separately at both normal photospheric velocities and at very high velocities. The high-velocity triplet absorption is highly polarized, with a different polarization angle than the rest of the spectrum. The unique observation allows us to construct a relatively detailed picture of the layered geometrical structure of the supernova ejecta: in our interpretation, the ejecta layers near the photosphere (v \approx 10,000 km/s) obey a near axial symmetry, while a detached, high-velocity structure (v \approx 18,000-25,000 km/s) with high CaII line opacity deviates from the photospheric axisymmetry. By partially obscuring the underlying photosphere, the high-velocity structure causes a more incomplete cancellation of the polarization of the photospheric light, and so gives rise to the polarization peak and rotated polarization angle of the high-velocity IR triplet feature. In an effort to constrain the ejecta geometry, we develop a technique for calculating 3-D synthetic polarization spectra and use it to generate polarization profiles for several parameterized configurations. In particular, we examine the case where the inner ejecta layers are ellipsoidal and the outer, high-velocity structure is one of four possibilities: a spherical shell, an ellipsoidal shell, a clumped shell, or a toroid. The synthetic spectra rule out the spherical shell model, disfavor a toroid, and find a best fit with the clumped shell. We show further that different geometries can be more clearly discriminated if observations are obtained from several different lines of sight.

We present the results of a large grid of synthetic spectra and compare them to early spectroscopic observations of SN 1993W. This supernova was discovered close to its explosion date and at a recession velocity of 5400 km/s is located in the Hubble flow. We focus here on two early spectra that were obtained approximately 5 and 9 days after explosion. We parameterize the outer supernova envelope as a power-law density profile in homologous expansion. In order to extract information on the value of the parameters a large number of models was required. We show that very early spectra combined with detailed models can provide constraints on the value of the power law index, the ratio of hydrogen to helium in the surface of the progenitor, the progenitor metallicity and the amount of radioactive nickel mixed into the outer envelope of the supernova. The spectral fits reproduce the observed spectra exceedingly well. The spectral results combined with the early photometry predict that the explosion date was 4.7 \pm 0.7 days before the first spectrum was obtained. The ability to obtain the metallicity from early spectra make SN IIP attractive probes of chemical evolution in the universe and by showing that we have the ability to pin down the parameters of the progenitor and mixing during the supernova explosion, it is likely to make SN IIP useful cosmological distance indicators which are at the same time complementary to SNe Ia.

We briefly describe the current version of the PHOENIX code. We then present some illustrative results from the modeling of Type Ia and Type II supernovae, hot stars, and irradiated giant planets. Good fits to observations can be obtained, when account is taken for spherically symmetric, line-blanketed, static or expanding atmospheres.

The usual interpetation of the spectrum of a BAL QSO is that a broad-band continuum from the central engine plus broad emission lines from a surrounding region (the BELR) emerges from near the center of the QSO, and broad absorption lines are superimposed in a separate outlying region (the BALR). For FBQS 1214+2803, designated as an FeLoBAL QSO because its spectrum contains numerous absorption features from excited states of Fe II, we explore an alternative interpretation based on resonance scattering. In this model line emission and absorption occur in the same line-forming region (the LFR). A resonance-scattering synthetic spectrum computed with the supernova synthetic-spectrum code SYNOW fits the spectrum of FBQS 1214+2802 rather well, so the resonance-scattering model merits further study. Some implications of the model and its possible applicability to other QSOs are briefly discussed.

The presence of a small amount of hydrogen is expected in most single degenerate scenarios for producing a Type Ia supernova (SN Ia). While hydrogen may be detected in very early high resolution optical spectra, in early radio spectra, and in X-ray spectra, here we examine the possibility of detecting hydrogen in early low resolution spectra such as those that will be obtained by proposed large scale searches for nearby SNe Ia. We find that definitive detections will require both very early spectra (less than 5 days after explosion) and perhaps slightly higher amounts of hydrogen than are currently predicted to be mixed into the outer layers of SNe Ia. Thus, the non-detection of hydrogen so far does not in and of itself rule out any current progenitor models. Nevertheless, very early spectra of SNe Ia will provide significant clues to the amount of hydrogen present and hence to the nature of the SN Ia progenitor system. Spectral coverage in both the optical and IR will be required to definitively identify hydrogen in low resolution spectra.

We describe Brute, a code for generating synthetic spectra from 3D supernova models. It is based on a Monte Carlo implementation of the Sobolev approximation. The current version is highly parameterized, but in the future it will be extended for full NLTE spectrum modelling of supernovae.

We describe the method we have used to parallelize our spherically symmetric special relativistic short characteristics general radiative transfer code PHOENIX. We describe some possible parallelization strategies and show why they would be inefficient. We discuss the multiple parallelization strategy techniques that we have adopted. We briefly discuss generalizing these strategies to full 3-D (spatial) radiation transfer codes.

Photometric and spectroscopic observations of the peculiar Type Ib Supernova 1991D are presented. SN 1991D was exceptionally bright for a Type Ib supernova. The He I lines were rather weak and the velocity at the photosphere as a function of time was unusually low. Comparison of the observed and synthetic spectra indicates that either hydrogen was ejected with a minimum velocity of 12,000 km/s or the spectrum contained features caused by lines of Ne I. Light curve modelling suggests that the progenitor probably had a very large radius (~10^14 cm) and that a considerable amount of 56Ni was synthesized during the explosion (~0.7 Msolar). We suggest a progenitor model of SN 1991D that involves a binary system.

Supernova 1987A remains the most well-studied supernova to date. Observations produced excellent broad-band photometric and spectroscopic coverage over a wide wavelength range at all epochs. We model the observed spectra from Day 1 to Day 81 using a hydrodynamical model. We show that good agreement can be obtained at times up to about 60 days, if we allow for extended nickel mixing. Later than about 60 days the observed Balmer lines become stronger than our models can reproduce. We show that this is likely due to a more complicated distribution of gamma-rays than we allow for in our spherically symmetric calculations. We present synthetic light curves in UBVRIJHK and a synthetic bolometric light curve. Using this broad baseline of detailed spectroscopic models we find a distance modulus mu = 18.5 +/- 0.2 using the SEAM method of determining distances to supernovae. We find that the explosion time agrees with that of the neutrino burst and is constrained at 68% confidence to within +/- 0.9 days. We argue that the weak Balmer lines of our detailed model calculations casts doubt on the accuracy of the purely photometric EPM method. We also suggest that Type IIP supernovae will be most useful as distance indicators at early times due to a variety of effects.

We report and interpret HST/STIS long-slit observations of the optical and ultraviolet (1150 - 10270 Angstrom) emission-line spectra of the rapidly brightening Spot 1 on the equatorial ring of SN 1987A between 1997 September and 1999 October (days 3869 -- 4606 after outburst). The emission is caused by radiative shocks created where the supernova blast wave strikes dense gas protruding inward from the equatorial ring. We measure and tabulate line identifications, fluxes and, in some cases, line widths and shifts. We compute flux correction factors to account for substantial interstellar line absorption of several emission lines. Nebular analysis shows that optical emission lines come from a region of cool (T_e ~ 10^4 K) and dense (n_e ~ 10^6 cm^-3) gas in the compressed photoionized layer behind the radiative shock. The observed line widths indicate that only shocks with shock velocities V_s < 250 km/s have become radiative, while line ratios indicate that much of the emission must have come from yet slower (V_s < 135 k/ms) shocks. We are able to fit the UV fluxes with an idealized radiative shock model consisting of two shocks (V_s = 135 and 250 km/s). The observed UV flux increase with time can be explained by the increase in shock surface areas as the blast wave overtakes more of the protrusion. The observed flux ratios of optical to highly-ionized UV lines are greater by a factor of ~ 2 -- 3 than predictions from the radiative shock models and we discuss the possible causes. We also present models for the observed H-alpha line widths and profiles, which suggests that a chaotic flow exists in the photoionized regions of these shocks. We discuss what can be learned with future observations of all the spots present on the equatorial ring.

The recent announcement that sodium absorption has been observed in the atmosphere of HD209458b, the only EGP observed to transit its parent star, is the first direct detection of an EGP atmosphere. We explore the possibility that neutral sodium is {\em not} in local thermodynamic equilibrium (LTE) in the outer atmosphere of irradiated EGPs and that the sodium concentration may be underestimated by models that make the LTE assumption. Our results indicate that it may not be necessary to invoke excessive photoionization, low metallicity, or even high altitude clouds to explain the observations.

A stellar wind module has been developed for the PHOENIX stellar atmosphere code for the purpose of computing non-LTE, line-blanketed, expanding atmospheric structures and detailed synthetic spectra of hot luminous stars with winds. We apply the code to observations of Deneb, for which we report the first positive detections of mm and cm emission (obtained using the SCUBA and the VLA), as well a strong upper limit on the 850 micron flux (using the HHT). The slope of the radio spectrum shows that the stellar wind is partially ionized. We report a uniform-disk angular diameter measurement, 2.40 +/- 0.06 mas, from the Navy Prototype Optical Interferometer (NPOI). The measured bolometric flux and corrected NPOI angular diameter yield an effective temperature of 8600 +/- 500 K. Least-squares comparisons of synthetic spectral energy distributions from 1220 A to 3.6 cm with the observations provide estimates for the effective temperature and the mass-loss rate of 8400 +/- 100 K and 8 +/- 3 E-7 M_sun/yr, respectively. This range of mass-loss rates is consistent with that derived from high dispersion UV spectra when non-LTE metal-line blanketing is considered. We are unable achieve a reasonable fit to a typical Halpha P-Cygni profile with any model parameters over a reasonable range. This is troubling because the \ha profile is the observational basis for Wind Momentum-Luminosity Relationship.

We present 0.9-2.5um spectroscopy with R~800 and 1.12-1.22um spectroscopy with R~5800 for the M dwarfs Gl 229A and LHS 102A, and for the L dwarf LHS 102B. We also report IZJHKL' photometry for both components of the LHS 102 system, and L' photometry for Gl 229A. The data are combined with previously published spectroscopy and photometry to produce flux distributions for each component of the kinematically old disk M/L-dwarf binary system LHS 102 and the kinematically young disk M/T-dwarf binary system Gliese 229. The data are analyzed using synthetic spectra generated by the latest "AMES-dusty" and "AMES-cond" models by Allard & Hauschildt. Although the models are not able to reproduce the overall slope of the infrared flux distribution of the L dwarf, most likely due to the treatment of dust in the photosphere, the data for the M dwarfs and the T dwarf are well matched. We find that the Gl 229 system is metal-poor despite having kinematics of the young disk, and that the LHS 102 system has solar metallicity. The observed luminosities and derived temperatures and gravities are consistent with evolutionary model predictions if the Gl 229 system is very young (age ~30 Myr) with masses (A,B) of (0.38,>0.007)M(sun), and the LHS 102 system is older, aged 1-10 Gyr with masses (A,B) of (0.19,0.07)M(sun).

The Asiago Supernova Catalog is used to carry out a comparative study of supernova absolute-magnitude distributions. An overview of the absolute magnitudes of the supernovae in the current observational sample is presented, and the evidence for subluminous and overluminous events is examined. The fraction of supernovae that are underluminous (M_B > -15) appears to be higher (perhaps much higher) than one fifth but it remains very uncertain. The fraction that are overluminous (M_B < -20) is lower (probably much lower) than 0.01. The absolute-magnitude distributions for each supernova type, restricted to events within 1 Gpc, are compared. Although these distributions are affected by observational bias in favor of the more luminous events, they are useful for comparative studies. We find mean absolute blue magnitudes (for H_0=60) of -19.46 for normal Type Ia supernovae (SNe Ia), -18.04 for SNe Ibc, -17.61 and -20.26 for normal and bright SNe Ibc considered separately, -18.03 for SNe II-L, -17.56 and -19.27 for normal and bright SNe II-L considered separately, -17.00 for SNe II-P, and -19.15 for SNe IIn.

Supernova flux and polarization spectra bring vital information on the geometry, physical conditions, and composition structure of the ejected matter. For some supernovae the circumstellar matter is also probed by the observed spectra. Some of this information can be inferred directly from the observed line profiles and fluxes, but because of the Doppler broadening and severe line blending, interpretation often involves the use of synthetic spectra. The emphasis in this Chapter is on recent results obtained with the help of synthetic spectra.

Observations of high-redshift Type Ia supernovae (SNe Ia) have provided strong evidence that the dark energy is real, and making further accurate observations of high-redshift SNe Ia is the most promising way to probe the nature of the dark energy. We discuss one of the concerns about such a project - that of coping with SN Ia evolution. We emphasize that SN Ia evolution differs in an important respect from the kind of evolution that has foiled some past projects in observational cosmology, and we outline empirical strategies that will take it into account. The supporting role of physical models of SNe Ia also is discussed. Our conclusion is that systematic errors due to SN Ia evolution will be small.

We show that the shape of P-Cygni line profiles of photospheric phase supernova can be analytically inverted to extract both the optical depth and source function of the line -- i.e. all the physical content of the model for the case when the Sobolev approximation is valid. Under various simplifying assumptions, we derive formulae that give S(r) and tau(r) in terms of derivatives of the line flux with respect to wavelength. The transition region between the minimum and maximum of the line profile turns out to give especially interesting information on the optical depth near the photosphere. The formulae give insights into the relationship between line shape and physical quantities that may be useful in interpreting observed spectra and detailed numerical calculations.

We investigate the prospects for constraining the maximum scale of clumping in composition that is consistent with observed Type Ia supernova flux spectra. Synthetic spectra generated without purely spherical composition symmetry indicate that gross asymmetries make prominent changes to absorption features. Motivated by this, we consider the case of a single unblended line forming in an atmosphere with perturbations of different scales and spatial distributions. Perturbations of about 1% of the area of the photodisk simply weaken the absorption feature by the same amount independent of the line of sight. Conversely, perturbations of about 10% of the area of the photodisk introduce variation in the absorption depth which does depend on the line of sight. Thus, 1% photodisk area perturbations may be consistent with observed profile homogeneity but 10% photodisk area perturbations can not. Based on this, we suggest that the absence of significant variation in the depths of Si II 6355 absorption features in normal Type Ia spectra near maximum light indicates that any composition perturbations in these events are quite small. This also constrains future three-dimensional explosion models to produce ejecta profiles with only small scale inhomogeneities.

We describe two parallel algorithms for line opacity calculations based on a local file and on a global file approach. The performance and scalability of both approaches is discussed for different test cases and very different parallel computing systems. The results show that a global file approach is more efficient on high-performance parallel supercomputers with dedicated parallel I/O subsystem whereas the local file approach is very useful on farms of workstations, e.g., cheap PC clusters.

We have fit the normal, well observed, Type Ia Supernova (SN Ia) SN 1994D with non-LTE spectra of the deflagration model W7. We find that well before maximum luminosity W7 fits the optical spectra of SN 1994D. After maximum brightness the quality of the fits weakens as the spectrum forms in a core rich in iron peak elements. We show the basic structure of W7 is likely to be representative of the typical SN Ia. We have shown that like W7, the typical SN Ia has a layer of unburned C+O composition at v > 15000 \kmps, followed by layers of C-burned and O-burned material with a density structure similar to W7. We present UVOIR (UBVRIJKH) synthetic photometry and colors and compare with observation. We have computed the distance to the host galaxy, NGC 4526, obtaining a distance modulus of \mu = 30.8 \pm 0.3. We discuss further application of this direct measurement of SNe Ia distances. We also discuss some simple modifications to W7 that could improve the quality of the fits to the observations.

Supernova 1987A remains the most well-observed and well-studied supernova to date. Observations produced excellent broad-band photometric and spectroscopic coverage over a wide wavelength range at all epochs. Here, we focus on the very early spectroscopic observations. Only recently have numerical models been of sufficient detail to accurately explain the observed spectra. In SN 1987A, good agreement has been found between observed and synthetic spectra for day one, but by day four, the predicted Balmer lines become much weaker than the observed lines. We present the results of work based on a radiation-hydrodynamic model by Blinnikov and collaborators. Synthetic non-LTE spectra generated from this model by the general radiation transfer code PHOENIX strongly support the theory that significant mixing of nickel into the outer envelope is required to maintain strong Balmer lines. Preliminary results suggest a lower limit to the average nickel mass of 1.0 \times 10^{-5} solar masses is required above 5000 \kmps by day four. PHOENIX models thus have the potential to be a sensitive probe for nickel mixing in the outer layers of a supernova.

We have calculated fast direct spectral model fits to two early-time spectra of the Type-II plateau SN 1999em, using the SYNOW synthetic spectrum code. The first is an extremely early blue optical spectrum and the second a combined HST and optical spectrum obtained one week later. Spectroscopically this supernova appears to be a normal Type II and these fits are in excellent agreement with the observed spectra. Our direct analysis suggests the presence of enhanced nitrogen. We have further studied these spectra with the full NLTE general model atmosphere code PHOENIX. While we do not find confirmation for enhanced nitrogen (nor do we rule it out), we do require enhanced helium. An even more intriguing possible line identification is complicated Balmer and He I lines, which we show falls naturally out of the detailed calculations with a shallow density gradient. We also show that very early spectra such as those presented here combined with sophisticated spectral modeling allows an independent estimate of the total reddening to the supernova, since when the spectrum is very blue, dereddening leads to changes in the blue flux that cannot be reproduced by altering the ``temperature'' of the emitted radiation. These results are extremely encouraging since they imply that detailed modeling of early spectra can shed light on both the abundances and total extinction of SNe II, the latter improving their utility and reliability as distance indicators.

We present spectral analysis of early observations of the Type IIn supernova 1998S using the general non-local thermodynamic equilibrium atmosphere code \tt PHOENIX}. We model both the underlying supernova spectrum and the overlying circumstellar interaction region and produce spectra in good agreement with observations. The early spectra are well fit by lines produced primarily in the circumstellar region itself, and later spectra are due primarily to the supernova ejecta. Intermediate spectra are affected by both regions. A mass-loss rate of order $\dot M \sim 0.0001-0.001$\msol yr$^{-1}$ is inferred for a wind speed of 100-1000 \kmps. We discuss how future self-consistent models will better clarify the underlying progenitor structure.

We present a grid of nova models that have an extremely large number of species treated in non-LTE, and apply it to the analysis of an extensive time series of ultraviolet spectroscopic data for Nova Cygni 1992. We use ultraviolet colors to derive the time development of the effective temperature of the expanding atmosphere during the fireball phase and the first ten days of the optically thick wind phase. We find that the nova has a pure optically thick wind spectrum until about 10 days after the explosion. During this interval, we find that synthetic spectra based on our derived temperature sequence agree very well with the observed spectra. We find that a sequence of hydrogen deficient models provides an equally good fit providing the model effective temperature is shifted upwards by ~1000 K. We find that high resolution UV spectra of the optically thick wind phase are fit moderately well by the models. We find that a high resolution spectrum of the fireball phase is better fit by a model with a steep density gradient, similar to that of a supernova, than by a nova model.

SN 1984A shows unusually large expansion velocities in lines from freshly synthesized material, relative to typical Type Ia Supernovae (SNe Ia). SN 1984A is an example of a group of SNe Ia which have very large blue-shifts of the P-Cygni features, but otherwise normal spectra. We have modeled several early spectra of SN 1984A with the multi-purpose NLTE model atmosphere and spectrum synthesis code, PHOENIX. We have used as input two delayed detonation models: CS15DD3 (Iwamoto et al. 1999) and DD21c (Hoeflich, Wheeler & Thielemann 1998). These models show line expansion velocities which are larger than that for a typical deflagration model like W7, which we have previously shown to fit the spectra of normal SNe Ia quite well. We find these delayed detonation models to be reasonable approximations to large absorption feature blue-shift SNe Ia, like SN 1984A. Higher densities of newly synthesized intermediate mass elements at higher velocities, v > 15,000 km/s, are found in delayed detonation models than in deflagration models. We find that this increase in density at high velocities is responsible for the larger blue-shifts in the synthetic spectra. We show that the variations in line width in observed SNe Ia are likely due to density variations in the outer, high-velocity layers of their atmospheres.

A comparison of the ratio of the depths of two absorption features in the spectra of TypeIa supernovae (SNe Ia) near the time of maximum brightness with the blueshift of the deep red Si II absorption feature 10 days after maximum shows that the spectroscopic diversity of SNe Ia is multi-dimensional. There is a substantial range of blueshifts at a given value of the depth ratio. We also find that the spectra of a sample of SNe Ia obtained a week before maximum brightness can be arranged in a ``blueshift sequence'' that mimics the time evolution of the pre-maximum-light spectra of an individual SN Ia, the well observed SN 1994D. Within the context of current SN Ia explosion models, we suggest that some of the SNe Ia in our sample were delayed-detonations while others were plain deflagrations.

We present a technique to treat systems with very many levels, like molecules, in non-LTE. This method is based on a superlevel formalism coupled with rate operator splitting. Superlevels consist of many individual levels that are assumed to be in LTE relative to each other. The usage of superlevels reduces the dimensionality of the rate equations dramatically and, thereby, makes the problem computationally more easily treatable. Our superlevel formalism retains maximum accuracy by using direct opacity sampling (dOS) when calculating the radiative transitions and the opacities. We developed this method in order to treat molecules in cool dwarf model calculations in non-LTE. Cool dwarfs have low electron densities and a radiation field that is far from a black body radiation field, both properties may invalidate the conditions for the common LTE approximation. Therefore, the most important opacity sources, the molecules, need to be treated in non-LTE. As a case study we applied our method to carbon monoxide. We find that our method gives accurate results since the conditions for the superlevel method are very well met for molecules. Due to very high collisional cross sections with hydrogen, and the high densities of H_2 the population of CO itself shows no significant deviation from LTE.

We have calculated a grid of photospheric phase atmospheres of Type Ia supernovae (SNe Ia) with metallicities from ten times to one thirtieth the solar metallicity in the C+O layer of the deflagration model, W7. We have modeled the spectra using the multi-purpose NLTE model-atmosphere and spectrum-synthesis code, PHOENIX. We show models for the epochs 7, 10, 15, 20, and 35 days after explosion. When compared to observed spectra obtained at the approximately corresponding epochs these synthetic spectra fit reasonably well. The spectra show variation in the overall level of the UV continuum with lower fluxes for models with higher metallicity in the unburned C+O layer. This is consistent with the classical surface cooling and line blocking effect due to metals in the outer layers of C+O. The UV features also move consistently to the blue with higher metallicity, demonstrating that they are forming at shallower and faster layers in the atmosphere. The potentially most useful effect is the blueward movement of the Si II feature at 6150 Angstrom with increasing C+O layer metallicity. We also demonstrate the more complex effects of metallicity variations by modifying the 54Fe content of the incomplete burning zone in W7 at maximum light. We briefly address some shortcomings of the W7 Finally, we identify that the split in the Ca H+K feature produced in W7 and observed in some SNe Ia is due to a blending effect of Ca II and Si II and does not necessarily represent a complex abundance or ionization effect in Ca II. amodel when compared to observations.

We present the extension of our NextGen model atmosphere grid to the regime of giant stars. The input physics of the models presented here is nearly identical to the NextGen dwarf atmosphere models, however spherical geometry is used self-consistently in the model calculations (including the radiative transfer). We re-visit the discussion of the effects of spherical geometry on the structure of the atmospheres and the emitted spectra and discuss the results of NLTE calculations for a few selected models.

We present detailed non-local thermodynamic equilibrium (NLTE) synthetic spectra for comparison with a time series of observed optical spectra of the Type Ic supernova 1994I which occurred in M51. With the exceptions of Si I and S I, we treat the important species in the formation of the spectrum in full NLTE. We present results for both a hydrodynamic model that has been fit to the light curve and for an illustrative custom crafted model that is more massive. Both models give reasonable fits to the overall observed spectra; however, neither is able to reproduce all the observed features. Some conspicuous observed features are absent and some predicted features are unobserved. No model that we have explored is able to satisfactorily reproduce the observed infrared feature near 1 micron on April 15, 1994 (+7d), which has been attributed to the triplet He I lambda 10830 transition. The low-mass hydrodynamic model produces an infrared feature with a blend of He I, C I, O I, and Si I--II lines, but it predicts a strong unobserved absorption feature near 6100 Angstroms due to Fe III, and the observed feature just blueward of 6000 Angstroms most likely due to Na D is not reproduced. The more massive model does a better job of reproducing the observed infrared lineshape, but also predicts the unobserved feature near 6100 Angstroms. The early-time spectrum of the low-mass model is far too blue; thus, a more massive model may be slightly favored. Since the predicted infrared feature is produced by a blend of so many elements and there is no overwhelming evidence for other helium features such as lambda 5876, it may be premature to conclude that SNe Ic unambiguously contain helium. Thus, we conclude that pure C+O cores are still viable progenitors for SNe Ic.

Synthetic spectra generated with the parameterized supernova synthetic-spectrum code SYNOW are compared to observed photospheric-phase spectra of the Type Ic supernova 1994I. The observed optical spectra can be well matched by synthetic spectra that are based on the assumption of spherical symmetry. We consider the identification of the infrared absorption feature observed near 10,250 \AA, which previously has been attributed to He I $\lambda10830$ and regarded as strong evidence that SN 1994I ejected some helium. We have difficulty accounting for the infrared absorption with He I alone. It could be a blend of He I and C I lines. Alternatively, we find that it can be fit by Si I lines without compromising the fit in the optical region. In synthetic spectra that match the observed spectra, from 4 days before to 26 days after the time of maximum brightness, the adopted velocity at the photosphere decreases from 17,500 to 7000 \kms. Simple estimates of the kinetic energy carried by the ejected mass give values that are near the canonical supernova energy of $10^{51}$ ergs. The velocities and kinetic energies that we find for SN 1994I in this way are much lower than those that we find elsewhere for the peculiar Type Ic SNe 1997ef and 1998bw, which therefore appear to have been hyper-energetic.

Synthetic spectra generated with the parameterized supernova synthetic-spectrum code SYNOW are compared to spectra of the Type Ia SN 1994D that were obtained before the time of maximum brightness. Evidence is found for the presence of two-component Fe II and Ca II features, forming in high velocity ($\ge 20,000$ \kms) and lower velocity ($\le 16,000$ \kms) matter. Possible interpretations of these spectral splits, and implications for using early--time spectra of SNe Ia to probe the metallicity of the progenitor white dwarf and the nature of the nuclear burning front in the outer layers of the explosion, are discussed.

We have used our PHOENIX multi-purpose model atmosphere code to calculate atmospheric models that represent novae in the optically thick wind phases of their outburst. We have improved the treatment of NLTE effects by expanding the number of elements that are included in the calculations from 15 to 19, and the number of ionization stages from 36 to 87. The code can now treat a total of 10713 levels and 102646 lines in NLTE. Aluminum, P, K, and Ni are included for the first time in the NLTE treatment and most elements now have at least the lowest six ionization stages included in the NLTE calculation. We have investigated the effects of expanded NLTE treatment on the chemical concentration of astrophysically significant species in the atmosphere, the equilibrium structure of the atmosphere, and the emergent flux distribution. Although we have found general qualitative agreement with previous, more limited NLTE models, the expanded NLTE treatment leads to significantly different values for the size of many of the NLTE deviations. In particular, for the hottest model presented here (effective temperature of 35000 K), for which NLTE effects are largest, we find that the expanded NLTE treatment reduces the NLTE effects for these important variables: neutral Hydrogen concentration, pressure structure, and emergent far UV flux. Moreover, we find that the addition of new NLTE species may greatly affect the concentration of species that were already treated in NLTE, so that, generally, all species that contribute significantly to the electron reservoir or to the total opacity, or whose line spectrum overlaps or interlocks with that of a species of interest, must be treated in NLTE to insure an accurate result for any particular species.

We calculate a grid of spherically symmetric OB stellar atmospheres at low metallicities, including both non-local thermodynamic equilibrium (NLTE) and metal line blanketing effects. This is done to assess the uncertainties in helium abundance determination by nebular codes due to input stellar atmosphere models. The more sophisticated stellar atmosphere models we use can differ from LTE models by as much as 40 percent in the ratio of He to H-ionizing photons.

A systematic survey of ions that could be responsible for features in the optical spectra of supernovae is carried out. Six different compositions that could be encountered in supernovae are considered. For each composition, the LTE optical depth of one of the strongest optical lines of each ion is plotted against temperature. For each ion that can realistically be considered as a candidate to produce identifiable features in supernova spectra, a sample synthetic spectrum is displayed. The optical depth plots and the synthetic spectra can provide guidance to studies of line identifications in the optical spectra of all types of supernovae during their photospheric phases.

In this paper we discuss numerical methods and algorithms for the solution of NLTE stellar atmosphere problems involving expanding atmospheres, e.g., found in novae, supernovae and stellar winds. We show how a scheme of nested iterations can be used to reduce the high dimension of the problem to a number of problems with smaller dimensions. As examples of these sub-problems, we discuss the numerical solution of the radiative transfer equation for relativistically expanding media with spherical symmetry, the solution of the multi-level non-LTE statistical equilibrium problem for extremely large model atoms, and our temperature correction procedure. Although modern iteration schemes are very efficient, parallel algorithms are essential in making large scale calculations feasible, therefore we discuss some parallelization schemes that we have developed.

We present our NextGen Model Atmosphere grid for low mass stars for effective temperatures larger than $3000\K$. These LTE models are calculated with the same basic model assumptions and input physics as the VLMS part of the NextGen grid so that the complete grid can be used, e.g., for consistent stellar evolution calculations and for internally consistent analysis of cool star spectra. This grid is also the starting point for a large grid of detailed NLTE model atmospheres for dwarfs and giants (Hauschildt et al, in preparation). The models were calculated from $3000\K$ to $10000\K$ (in steps of $200\K$) for $3.5 \le \logg \le 5.5$ (in steps of 0.5) and metallicities of $-4.0 \le \mh \le 0.0$. We discuss the results of the model calculations and compare our results to the Kurucz 1994 grid. Some comparisons to standard stars like Vega and the Sun are presented and compared with detailed NLTE calculations.

A parameterized supernova synthetic-spectrum code is used to study line identifications in the photospheric-phase spectra of the peculiar Type Ia SN 1991T, and to extract some constraints on the composition structure of the ejected matter. The inferred composition structure is not like that of any hydrodynamical model for Type Ia supernovae. Evidence that SN 1991T was overluminous for an SN Ia is presented, and it is suggested that this peculiar event probably was a substantially super-Chandrasekhar explosion that resulted from the merger of two white dwarfs.

The observed map of 1.809 MeV gamma-rays from radioactive Al-26 (Oberlack et al, 1996) shows clear evidence of Galactic plane origin with an uneven distribution. We have simulated the map using a Monte Carlo technique together with simple assumptions about the spatial distributions and yields of Al-26 sources (clustered core-collapse supernovae and Wolf Rayet stars; low- and high-mass AGB stars; and novae). Although observed structures (e.g., tangents to spiral arms, bars, and known star-forming regions) are not included in the model, our simulated gamma-ray distribution bears resemblance to the observed distribution. The major difference is that the model distribution has a strong smooth background along the Galactic plane from distant sources in the disk of the Galaxy. We suggest that the smooth background is to be expected, and probably has been suppressed by background subtraction in the observed map. We have also found an upper limit of 1 Msun to the contribution of flux from low-yield, smoothly distributed sources (low-mass AGB stars and novae).

Until we can more carefully assess the validity of the physical assumptions made we withdraw this paper.

We use a spherical non-LTE fully line blanketed model atmosphere to fit the full multi-wavelength spectrum, including the extreme ultraviolet (EUV) continuum observed by the {\it Extreme Ultraviolet Explorer}, of the B2 II star \epscma. The available spectrophotometry of \epscma\ from 350 \AA\ to 25 \micron\ is best fit with model parameters \Teff = 21750\,K, \Logg = 3.5, and an angular diameter of 0.77 mas. The close agreement between the model and the measured EUV flux from \epscma\ is a result of the higher temperatures at the formation depths of the \ion{H}{1} and \ion{He}{1} Lyman continua compared to other models. The realistic model treatment of early B giants with spherical geometry and NLTE metal line blanketing results in the prediction of significantly larger EUV fluxes compared with plane-parallel models. We find that our metal line blanketed spherical models show significantly warmer temperature structures, 1-3 kK at the formation depth of the Lyman continua, and predict stronger EUV fluxes, up to a factor of 5 in the \ion{H}{1} Lyman continuum, compared with plane-parallel atmospheres that have identical model parameters. In contrast, we find spherical and plane-parallel models that do not include metal line blanketing are nearly identical. Our \Teff = 21000 K, \Logg = 3.2, spherical NLTE model predicts more than twice as many hydrogen ionizing photons and over 200 times more neutral helium ionizing photons than a standard hydrostatic plane-parallel LTE model with the same stellar parameters.

We describe an important addition to the parallel implementation of our generalized NLTE stellar atmosphere and radiative transfer computer program PHOENIX. In a previous paper in this series we described data and task parallel algorithms we have developed for radiative transfer, spectral line opacity, and NLTE opacity and rate calculations. These algorithms divided the work spatially or by spectral lines, that is distributing the radial zones, individual spectral lines, or characteristic rays among different processors and employ, in addition task parallelism for logically independent functions (such as atomic and molecular line opacities). For finite, monotonic velocity fields, the radiative transfer equation is an initial value problem in wavelength, and hence each wavelength point depends upon the previous one. However, for sophisticated NLTE models of both static and moving atmospheres needed to accurately describe, e.g., novae and supernovae, the number of wavelength points is very large (200,000--300,000) and hence parallelization over wavelength can lead both to considerable speedup in calculation time and the ability to make use of the aggregate memory available on massively parallel supercomputers. Here, we describe an implementation of a pipelined design for the wavelength parallelization of PHOENIX, where the necessary data from the processor working on a previous wavelength point is sent to the processor working on the succeeding wavelength point as soon as it is known. Our implementation uses a MIMD design based on a relatively small number of standard MPI library calls and is fully portable between serial and parallel computers.

We examine two scenarios for formation of the planetary nebula K648: a prompt scenario where the planetary nebula is ejected and formed immediately after a helium shell flash and a delayed scenario where a third dredge up occurs and the envelope is ejected during the following interpulse phase. We present models of both scenarios and find that each can produce K648-like systems. We suggest that the prompt scenario is more favorable but cannot rule out the delayed scenario.

We present new detailed NLTE calculations for model atmospheres of novae during outburst. This fully self-consistent NLTE treatment for a number of model atoms includes 3922 NLTE levels and 47061 NLTE primary transitions. We discuss the implication of departures from LTE for the strengths of the lines in nova spectra. The new results show that our large set of NLTE lines constitute the majority of the total line blanketing opacity in nova atmospheres. Although we include LTE background lines, their effect are small on the model structures and on the synthetic spectra. We demonstrate that the assumption of LTE leads to incorrect synthetic spectra and that NLTE calculations are required for reliably modeling nova spectra. In addition, we show that detailed NLTE treatment for a number of ionization stages of iron changes the results of previous calculations and improve the fit to observed nova spectra.

We present detailed NLTE Ti~I calculations in model atmospheres of cool dwarf and giant stars. A fully self-consistent NLTE treatment for a Ti~I model atom with 395 levels and 5279 primary bound-bound transitions is included, and we discuss the implication of departures from LTE in this atom for the strengths of Ti~I lines and TiO molecular bands in cool star spectra. We show that in the atmospheric parameter range investigated, LTE is a poor approximation to Ti~I line formation, as expected from the low collisional rates in cool stars. The secondary effects of Ti~I overionization on the TiO number density and the TiO molecular opacities, however, are found to be negligible in the molecular line forming region for the relatively small parameter range studied in this paper.

We describe the parallel implementation of our generalized stellar atmosphere and NLTE radiative transfer computer program PHOENIX. We discuss the parallel algorithms we have developed for radiative transfer, spectral line opacity, and NLTE opacity and rate calculations. Our implementation uses a MIMD design based on a relatively small number of MPI library calls. We report the results of test calculations on a number of different parallel computers and discuss the results of scalability tests.

We present detailed NLTE synthetic spectra of hydrodynamic SNe Ia models. We make no assumptions about the form of the spectrum at the inner boundary. We calculate both Chandrasekhar-mass deflagration models and sub-Chandrasekhar ``helium detonators.'' Gamma-ray deposition is handled in a simple, accurate manner. We have parameterized the storage of energy that arises from the time dependent deposition of radioactive decay energy in a reasonable manner, that spans the expected range. We find that the Chandrasekhar-mass deflagration model W7 of Nomoto etal shows good agreement with the observed spectra of SN 1992A and SN 1994D, particularly in the UV, where our models are expected to be most accurate. The sub-Chandrasekhar models do not reproduce the UV deficit observed in normal SNe Ia. They do bear some resemblance to sub-luminous SNe Ia, but the shape of the spectra (i.e. the colors) are opposite to that of the observed ones and the intermediate mass element lines such as Si II, and Ca II are extremely weak, which seems to be a generic difficulty of the models. Although the sub-Chandrasekhar models have a significant helium abundance (unlike Chandrasekhar-mass models), helium lines are not prominent in the spectra near maximum light and thus do not act as a spectral signature for the progenitor.

We have analyzed the early optically thick ultraviolet spectra of Nova OS And 1986 using a grid of spherically symmetric, non-LTE, line-blanketed, expanding model atmospheres and synthetic spectra with the following set of parameters: $5,000\le$ T$_{model}$ $\le 60,000$K, solar abundances, $\rho \propto r^{-3}$, $\v_{max} = 2000\kms$, $L=6 \times 10^{4}\Lsun$, and a statistical or microturbulent velocity of 50 $\kms$. We used the synthetic spectra to estimate the model parameters corresponding to the observed {\it IUE} spectra. The fits to the observations were then iteratively improved by changing the parameters of the model atmospheres, in particular T$_{model}$ and the abundances, to arrive at the best fits to the optically thick pseudo-continuum and the features found in the {\it IUE} spectra. The {\it IUE} spectra show two different optically thick subphases. The earliest spectra, taken a few days after maximum optical light, show a pseudo-continuum created by overlapping absorption lines. The later observations, taken approximately 3 weeks after maximum light, show the simultaneous presence of allowed, semi-forbidden, and forbidden lines in the observed spectra. Analysis of these phases indicate that OS And 86 had solar metallicities except for Mg which showed evidence of being underabundant by as much as a factor of 10. We determine a distance of 5.1 kpc to OS And 86 and derive a peak bolometric luminosity of $\sim$ 5 $\times$ 10$^4$ L$_{\odot}$. The computed nova parameters provide insights into the physics of the early outburst and explain the spectra seen by {\it IUE}. Lastly, we find evidence in the later observations for large non-LTE effects of Fe{\sc ii} which, when included, lead to much better agreement with the observations.

Supernovae, with their diversity of compositions, velocities, envelope masses, and interactions are good testing grounds for probing the importance of NLTE in expanding atmospheres. In addition to treating H, He, Li I, O I, Ne I, Na I, and Mg II in NLTE, we use a very large model atom of Fe II to test the importance of NLTE processes in both SNe Ia and II. Since the total number of potential line transitions that one has to include is enormous ($\approx 40$ million), approximations and simplifications are required to treat the problem accurately and in finite computer time. With our large Fe II model atom (617 levels, 13,675 primary NLTE line transitions) we are able to test several assumptions for treating the background opacity that are needed to obtain correct UV line blanketing which determines the shape of near-maximum light supernova spectra. We find that, due to interactions within the multiplets, treating the background lines as pure scattering (thermalization parameter $\epsilon = 0$) is a poor approximation and that an overall mean value of $\epsilon \sim 0.05 - 0.10 $ is a far better approximation. This is true even in SNe Ia, where the continuum absorption optical depth at 5000 \AA\ ($\equiv \tau_{\rm std}$) is $<< 1$. We also demonstrate that a detailed treatment of NLTE effects is required to properly determine the ionization states of both abundant and trace elements.

An interesting way to calibrate the absolute magnitudes of remote Type Ia supernovae (SNe Ia) that are well out in the Hubble flow, and thus determine the value of the Hubble constant, H_0, has been introduced by van den Bergh. His approach relies on calculations of the peak absolute magnitudes and broad--band colors for SN Ia explosion models. It does not require any corrections for extinction by interstellar dust, and no SNe Ia are excluded on grounds of peculiarity. Within the last few years distances have been determined to the parent galaxies of six SNe Ia by means of Cepheid variables. Cepheid--based distances also have become available for three other SNe Ia if one is willing to use the distance to a galaxy in the same group in lieu of the distance to the parent galaxy itself. Here we determine the value of H_0 in a way that is analogous to that of van den Bergh, but now using Cepheid--based distances instead of calculated light curves. We obtain H_0 = 55 km/s/Mpc. This value, with Lambda=0 and Omega=1, corresponds to a cosmic expansion time of 12 Gyr, which is consistent with several recent determinations of the ages of globular clusters.

The atmospheres of M stars are dominated by a small number of very strong molecular compounds (H$_2$O, TiO, H$_2$, CO, VO). Most of the hydrogen is locked in molecular H$_2$, most of the carbon in CO; and H$_2$O, TiO and VO opacities define a pseudo-continuum covering the entire flux distribution of these stars. The optical ``continuum'' is due to TiO vibrational bands which are often used as temperature indicators for these stars. These may be the depth of the bands relative to the troughs in between them; or the depth of the VO bands; or of the atomic lines relative to the local ``continuum''; or even the strength of the infrared water bands; all of these depend on the strength of the TiO bands and the amount of flux-redistribution to longer wavelengths exerted by them. Departures from LTE of the Ti I atom, and thus the concentration of the important TiO molecule, could, therefore, have severe and measurable consequences on the atmospheric structure and spectra of these stars. In this paper we discuss NLTE effects of Ti I in fully self-consistent models for a few representative M/Brown dwarf and M giant model atmospheres and spectra.

We describe our implementation of an extremely detailed model atom of singly ionized iron for NLTE computations in static and moving astrophysical plasmas. Our model atom includes 617 levels, 13675 primary permitted transitions and up to 1.2 million secondary transitions. Our approach guarantees that the total iron opacity is included at the correct wavelength with reasonable memory and CPU requirements. We find that the lines saturate the wavelength space, such that special wavelength points inserted along the detailed profile functions may be replaced with a statistical sampling method. We describe the results of various test calculations for novae and supernovae.

The atmospheres of novae at early times in their outbursts are very extended, expanding shells with low densities. Models of these atmospheres show that NLTE effects are very important and must be included in realistic calculations. We have, therefore, been improving our atmospheric studies by increasing the number of ions treated in NLTE. One of the most important ions is Fe II which has a complex structure and numerous lines in the observable spectrum. In this paper we investigate NLTE effects for Fe II for a wide variety of parameters. We use a detailed Fe II model atom with 617 level and 13675 primary lines, treated using a rate-operator formalism. We show that the radiative transfer equation in nova atmospheres {\em must} be treated with sophisticated numerical methods and that simple approximations, such as the Sobolev method, {\em cannot} be used because of the large number of overlapping lines in the co-moving frame. Our results show that the formation of the Fe II lines is strongly affected by NLTE effects. For low effective temperatures, $\Teff < 20,000\,$K, the optical Fe II lines are most influenced by NLTE effects, while for higher $\Teff$ the UV lines of Fe II are very strongly affected by NLTE. The departure coefficients are such that Fe II tends to be overionized in NLTE when compared to LTE. Therefore, Fe II-NLTE must be included with sophisticated radiative transfer in nova atmosphere models in order to reliably analyze observed nova spectra.

Statistical connections between the properties of Type Ia supernovae (SNe Ia) and the B--V colors of their parent galaxies are established. Compared to SNe Ia in blue galaxies [B-V$\la$0.75], SNe Ia in redder galaxies have (1) a wider dispersion in the blueshifts of their Si II $\lambda$6355 absorption features, ten days after maximum light; (2) more rapidly declining light curves; and (3) lower luminosities. Even when the spectroscopically peculiar, very subluminous SNe Ia such as SN 1991bg are disregarded, SNe Ia in red galaxies are less luminous than those in blue galaxies by about 0.3 magnitudes. When SNe Ia that are thought to have been significantly extinguished by dust in their parent galaxies are disregarded, those in blue galaxies have observational absolute--magnitude dispersions of only $\sigma_{obs}(M_B)$=0.20 and $\sigma_{obs}(M_V)$=0.17, which implies that their intrinsic absolute--magnitude dispersions are very small. We use six SNe Ia whose absolute magnitudes have been calibrated by means of Cepheids, which also indicate that the intrinsic absolute--magnitude dispersions of SNe Ia in blue galaxies are very small, to calibrate SNe Ia in blue galaxies and obtain $\rm H_0=57\pm4\ km\ s^{-1}\ Mpc^{-1}$. This value is in excellent agreement with that obtained by Saha et al. (1995b), in spite of the fact that they do not take into account any dependence of SN Ia absolute magnitude on the nature of the parent galaxy. Some implications of the statistical connections between SNe Ia and the colors of their parent galaxies, for identifying the progenitor binary systems of SNe Ia and for using high--redshift SNe Ia to measure q$_0$, are briefly discussed.

Modeling the atmospheres of SNe~Ia requires the solution of the NLTE radiative transfer equation. We discuss the formulation of the radiative transfer equation in the co-moving frame. For characteristic velocities larger than about 2000 km/s, the effects of advection on the synthetic spectra are non-negligible, and hence should be included in model calculations. We show that the time-independent or quasi-static approximation is adequate for SNe~Ia near maximum light, as well as for most other astrophysical problems; e.g., hot stars, novae, and other types of supernovae. We examine the use of the Sobolev approximation in modeling moving atmospheres and find that the number of overlapping lines in the co-moving frame make the approximation suspect in models that predict both lines and continua. We briefly discuss the form of the Rosseland mean opacity in the co-moving frame, and present a formula that is easy to implement in radiation hydrodynamics calculations.

We present optical spectra of the Type Ic supernova 1994I in M51 and preliminary non-LTE analysis of the spectra. Our models are not inconsistent with the explosions of C+O cores of massive stars. While we find no direct evidence for helium in the optical spectra, our models cannot rule out small amounts of helium. More than 0.1~\msol\ of helium seems unlikely.

In this Letter we present evidence for a spectral sequence among Type Ia supernovae (SNe Ia). The sequence is based on the systematic variation of several features seen in the near-maximum light spectrum. This sequence is analogous to the recently noted photometric sequence among SNe Ia which shows a relationship between the peak brightness of a SN Ia and the shape of its light curve. In addition to the observational evidence we present a partial theoretical explanation for the sequence. This has been achieved by producing a series of non-LTE synthetic spectra in which only the effective temperature is varied. The synthetic sequence nicely reproduces most of the differences seen in the observed one and presumably corresponds to the amount of 56Ni produced in the explosion.