With our robotic 1.2 m TIGRE telescope, we were able to obtain eight optical spectra with intermediate resolution (R = 20,000) of the Nova V659 Sct during different phases of its outburst. We present a list of the lines found in the Nova spectra. The most common features are H I, O I, Na I, Fe II and Ca II. Studying the spectral evolution of the strong features we found that the absorption features move to higher expansion velocities before disappearing and the emission features show (different) asymmetries. Thanks to the intermediate spectral resolution we identified and analysed the interstellar medium absorption features present in the spectra. We detected atomic absorption features of Na I and Ca II. The sodium D lines show more complex substructures with three main absorption features at around a velocity of -10, 30 and 85 km/s. We identified several DIBs in the Nova V659 Sct spectra and determined their velocities and equivalent widths.

The He I infrared (IR) line at a vacuum wavelength of 10833 A is a diagnostic for the investigation of atmospheres of stars and planets orbiting them. For the first time, we study the behavior of the He I IR line in a set of chromospheric models for M-dwarf stars, whose much denser chromospheres may favor collisions for the level population over photoionization and recombination, which are believed to be dominant in solar-type stars. For this purpose, we use published PHOENIX models for stars of spectral types M2 V and M3 V and also compute new series of models with different levels of activity following an ansatz developed for the case of the Sun. We perform a detailed analysis of the behavior of the He I IR line within these models. We evaluate the line in relation to other chromospheric lines and also the influence of the extreme ultraviolet (EUV) radiation field. The analysis of the He I IR line strengths as a function of the respective EUV radiation field strengths suggests that the mechanism of photoionization and recombination is necessary to form the line for inactive models, while collisions start to play a role in our most active models. Moreover, the published model set, which is optimized in the ranges of the Na I D2, H$\alpha$, and the bluest Ca II IR triplet line, gives an adequate prediction of the He I IR line for most stars of the stellar sample. Because especially the most inactive stars with weak He I IR lines are fit worst by our models, it seems that our assumption of a 100% filling factor of a single inactive component no longer holds for these stars.

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.

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

Context. Teegarden's Star is the brightest and one of the nearest ultra-cool dwarfs in the solar neighbourhood. For its late spectral type (M7.0V), the star shows relatively little activity and is a prime target for near-infrared radial velocity surveys such as CARMENES. Aims. As part of the CARMENES search for exoplanets around M dwarfs, we obtained more than 200 radial-velocity measurements of Teegarden's Star and analysed them for planetary signals. Methods. We find periodic variability in the radial velocities of Teegarden's Star. We also studied photometric measurements to rule out stellar brightness variations mimicking planetary signals. Results. We find evidence for two planet candidates, each with $1.1M_\oplus$ minimum mass, orbiting at periods of 4.91 and 11.4 d, respectively. No evidence for planetary transits could be found in archival and follow-up photometry. Small photometric variability is suggestive of slow rotation and old age. Conclusions. The two planets are among the lowest-mass planets discovered so far, and they are the first Earth-mass planets around an ultra-cool dwarf for which the masses have been determined using radial velocities.

We present precise photospheric parameters of 282 M dwarfs determined from fitting the most recent version of PHOENIX models to high-resolution CARMENES spectra in the visible (0.52 - 0.96 $\mu$m) and near-infrared wavelength range (0.96 - 1.71 $\mu$m). With its aim to search for habitable planets around M dwarfs, several planets of different masses have been detected. The characterization of the target sample is important for the ability to derive and constrain the physical properties of any planetary systems that are detected. As a continuation of previous work in this context, we derived the fundamental stellar parameters effective temperature, surface gravity, and metallicity of the CARMENES M-dwarf targets from PHOENIX model fits using a $\chi^2$ method. We calculated updated PHOENIX stellar atmosphere models that include a new equation of state to especially account for spectral features of low-temperature stellar atmospheres as well as new atomic and molecular line lists. We show the importance of selecting magnetically insensitive lines for fitting to avoid effects of stellar activity in the line profiles. For the first time, we directly compare stellar parameters derived from multiwavelength range spectra, simultaneously observed for the same star. In comparison with literature values we show that fundamental parameters derived from visible spectra and visible and near-infrared spectra combined are in better agreement than those derived from the same spectra in the near-infrared alone.

Modeling the submillimeter to centimeter emission of stars is challenging due to a lack of sensitive observations at these long wavelengths. We launched an ongoing campaign to obtain new observations entitled Measuring the Emission of Stellar Atmospheres at Submillimeter/millimeter wavelengths (MESAS). Here we present ALMA, GBT, and VLA observations of Sirius A, the closest main-sequence A-type star, that span from 1.4 to 9.0 millimeters. These observations complement our previous millimeter data on Sirius A and are entirely consistent with the PHOENIX stellar atmosphere models constructed to explain them. We note that accurate models of long wavelength emission from stars are essential not only to understand fundamental stellar processes, but also to determine the presence of dusty debris in spatially unresolved observations of circumstellar disks.

Chromospheric modeling of observed differences in stellar activity lines is imperative to fully understand the upper atmospheres of late-type stars. We present one-dimensional parametrized chromosphere models computed with the atmosphere code PHOENIX using an underlying photosphere of 3500 K. The aim of this work is to model chromospheric lines of a sample of 50 M2-3 dwarfs observed in the framework of the CARMENES, the Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs, exoplanet survey. The spectral comparison between observed data and models is performed in the chromospheric lines of Na I D2, H$\alpha$, and the bluest Ca II infrared triplet line to obtain best-fit models for each star in the sample. We find that for inactive stars a single model with a VAL C-like temperature structure is sufficient to describe simultaneously all three lines adequately. Active stars are rather modeled by a combination of an inactive and an active model, also giving the filling factors of inactive and active regions. Moreover, the fitting of linear combinations on variable stars yields relationships between filling factors and activity states, indicating that more active phases are coupled to a larger portion of active regions on the surface of the star.

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.

Our understanding of stellar atmospheres and our ability to infer architectures of extrasolar planetary systems rely on understanding the emission of stars at submillimeter to centimeter wavelengths. In this chapter we describe how unconstrained stellar emission can interfere with the accurate characterization of circumstellar debris. The ngVLA is the only facility with the sensitivity that allows for the observations of a broad range of stellar spectral types in a feasible amount of time. The observations will enable the building and testing of accurate models of stellar emission, which in turn are required for evaluating both the occurrence and abundance of debris over the proposed wavelength range of the ngVLA

In the early stages of planet formation, small dust grains grow to become mm sized particles in debris disks around stars. These disks can in principle be characterized by their emission at submillimeter and millimeter wavelengths. Determining both the occurrence and abundance of debris in unresolved circumstellar disks of A-type main-sequence stars requires that the stellar photospheric emission be accurately modeled. To better constrain the photospheric emission for such systems, we present observations of Sirius A, an A-type star with no known debris, from the JCMT, SMA, and VLA at 0.45, 0.85, 0.88, 1.3, 6.7, and 9.0 mm. We use these observations to inform a PHOENIX model of Sirius A's atmosphere. We find the model provides a good match to these data and can be used as a template for the submm/mm emission of other early A-type stars where unresolved debris may be present. The observations are part of an ongoing observational campaign entitled Measuring the Emission of Stellar Atmospheres at Submm/mm wavelengths (MESAS)

We present the results of local, vertically stratified, radiation magnetohydrodynamic shearing box simulations of magnetorotational instability (MRI) turbulence for a (hydrogen poor) composition applicable to accretion disks in AM CVn type systems. Many of these accreting white dwarf systems are helium analogues of dwarf novae (DNe). We utilize frequency-integrated opacity and equation of state tables appropriate for this regime to accurately portray the relevant thermodynamics. We find bistability of thermal equilibria in the effective temperature, surface mass density plane typically associated with disk instabilities. Along this equilibrium curve (i.e. the S-curve) we find that the stress to thermal pressure ratio $\alpha$ varied with peak values of $\sim 0.15$ near the tip of the upper branch. Similar to DNe, we found enhancement of $\alpha$ near the tip of the upper branch caused by convection; this increase in $\alpha$ occurred despite our choice of zero net vertical magnetic flux. Two notable differences we find between DN and AM CVn accretion disk simulations are that AM CVn disks are capable of exhibiting persistent convection in outburst, and ideal MHD is valid throughout quiescence for AM CVns. In contrast, DNe simulations only show intermittent convection, and non-ideal MHD effects are likely important in quiescence. By combining our previous work with these new results, we also find that convective enhancement of the MRI is anticorrelated with mean molecular weight.

The CARMENES radial velocity (RV) survey is observing 324 M dwarfs to search for any orbiting planets. In this paper, we present the survey sample by publishing one CARMENES spectrum for each M dwarf. These spectra cover the wavelength range 520--1710nm at a resolution of at least $R > 80,000$, and we measure its RV, H$\alpha$ emission, and projected rotation velocity. We present an atlas of high-resolution M-dwarf spectra and compare the spectra to atmospheric models. To quantify the RV precision that can be achieved in low-mass stars over the CARMENES wavelength range, we analyze our empirical information on the RV precision from more than 6500 observations. We compare our high-resolution M-dwarf spectra to atmospheric models where we determine the spectroscopic RV information content, $Q$, and signal-to-noise ratio. We find that for all M-type dwarfs, the highest RV precision can be reached in the wavelength range 700--900nm. Observations at longer wavelengths are equally precise only at the very latest spectral types (M8 and M9). We demonstrate that in this spectroscopic range, the large amount of absorption features compensates for the intrinsic faintness of an M7 star. To reach an RV precision of 1ms$^{-1}$ in very low mass M dwarfs at longer wavelengths likely requires the use of a 10m class telescope. For spectral types M6 and earlier, the combination of a red visual and a near-infrared spectrograph is ideal to search for low-mass planets and to distinguish between planets and stellar variability. At a 4m class telescope, an instrument like CARMENES has the potential to push the RV precision well below the typical jitter level of 3-4ms$^{-1}$.

Context: The main goal of the CARMENES survey is to find Earth-mass planets around nearby M-dwarf stars. Seven M-dwarfs included in the CARMENES sample had been observed before with HIRES and HARPS and either were reported to have one short period planetary companion (GJ15A, GJ176, GJ436, GJ536 and GJ1148) or are multiple planetary systems (GJ581 and GJ876). Aims: We aim to report new precise optical radial velocity measurements for these planet hosts and test the overall capabilities of CARMENES. Methods: We combined our CARMENES precise Doppler measurements with those available from HIRES and HARPS and derived new orbital parameters for the systems. Bona-fide single planet systems are fitted with a Keplerian model. The multiple planet systems were analyzed using a self-consistent dynamical model and their best fit orbits were tested for long-term stability. Results: We confirm or provide supportive arguments for planets around all the investigated stars except for GJ15A, for which we find that the post-discovery HIRES data and our CARMENES data do not show a signal at 11.4 days. Although we cannot confirm the super-Earth planet GJ15Ab, we show evidence for a possible long-period ($P_{\rm c}$ = 7025$_{-629}^{+972}$ d) Saturn-mass ($m_{\rm c} \sin i$ = 51.8$_{-5.8}^{+5.5}M_\oplus$) planet around GJ15A. In addition, based on our CARMENES and HIRES data we discover a second planet around GJ1148, for which we estimate a period $P_{\rm c}$ = 532.6$_{-2.5}^{+4.1}$ d, eccentricity $e_{\rm c}$ = 0.34$_{-0.06}^{+0.05}$ and minimum mass $m_{\rm c} \sin i$ = 68.1$_{-2.2}^{+4.9}M_\oplus$. Conclusions: The CARMENES optical radial velocities have similar precision and overall scatter when compared to the Doppler measurements conducted with HARPS and HIRES. We conclude that CARMENES is an instrument that is up to the challenge of discovering rocky planets around low-mass stars.

We report on the first star discovered to host a planet detected by radial velocity (RV) observations obtained within the CARMENES survey for exoplanets around M dwarfs. HD 147379 ($V = 8.9$ mag, $M = 0.58 \pm 0.08$ M$_{\odot}$), a bright M0.0V star at a distance of 10.7 pc, is found to undergo periodic RV variations with a semi-amplitude of $K = 5.1\pm0.4$ m s$^{-1}$ and a period of $P = 86.54\pm0.06$ d. The RV signal is found in our CARMENES data, which were taken between 2016 and 2017, and is supported by HIRES/Keck observations that were obtained since 2000. The RV variations are interpreted as resulting from a planet of minimum mass $m_{\rm p}\sin{i} = 25 \pm 2$ M$_{\oplus}$, 1.5 times the mass of Neptune, with an orbital semi-major axis $a = 0.32$ au and low eccentricity ($e < 0.13$). HD 147379b is orbiting inside the temperate zone around the star, where water could exist in liquid form. The RV time-series and various spectroscopic indicators show additional hints of variations at an approximate period of 21.1d (and its first harmonic), which we attribute to the rotation period of the star.

We present results of our dense spectroscopic monitoring of Nova V5668 Sgr. Starting on March 19 in 2015, only a few days after discovery, we have obtained a series of spectra with the TIGRE telescope and its HEROS echelle spectrograph which offers a resolution of R = 20,000 and covers the optical wavelength range from 3800 to 8800 {\AA}. We performed a line identification of the discernible features for four spectra which are representative for the respective phases in the light curve evolution of that nova. By simultaneously analysing the variations in the visual light curve and the corresponding spectra of Nova V5668 Sgr, we found that during the declining phases of the nova the absorption features in all hydrogen and many other lines had shifted to higher expansion velocities of -2000 km s^-1. Conversely, during the rise towards the following maximum, these observed absorption features had returned to lower expansion velocities.We found that the absorption features of some Fe II lines displayed the same behaviour, but in addition disappeared for a few days during some declining phases. Features of several N I lines also disappeared while new N II lines appeared in emission for a few days during some of the declining phases of the light curve of Nova V5668 Sgr. The shape of the emission features is changing during the evolution and shows a clear double peak structure after the deep minimum.

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.

Brown dwarfs and giant planets orbiting close to a host star are subjected to significant irradiation that can modify the properties of their atmospheres. In order to test the atmospheric models that are used to describe these systems, it is necessary to obtain accurate observational estimates of their physical properties (masses, radii, temperatures, albedos). Interacting compact binary systems provide a natural laboratory for studying strongly irradiated sub-stellar objects. As the mass-losing secondary in these systems makes a critical, but poorly understood transition from the stellar to the sub-stellar regime, it is also strongly irradiated by the compact accretor. In fact, the internal and external energy fluxes are both expected to be comparable in these objects, providing access to an unexplored irradiation regime. However, the atmospheric properties of such donors have so far remained largely unknown. Here, we report the direct spectroscopic detection and characterisation of an irradiated sub-stellar donor in an accreting white dwarf binary system. Our near-infrared observations allow us to determine a model-independent mass estimate for the donor of $M_2=0.055\pm0.008M_{\odot}$ and an average spectral type of ${\rm L1}\pm{\rm1}$, supporting both theoretical predictions and model-dependent observational constraints. Our time-resolved data also allow us to estimate the average irradiation-induced temperature difference between the day and night sides on the sub-stellar donor, $\Delta{\rm T}\simeq57$~K, and the maximum difference between the hottest and coolest parts of its surface, of $\Delta{\rm T}_{max}\simeq200$~K. The observations are well described by a simple geometric reprocessing model with a bolometric (Bond) albedo of $A_B<0.54$ at 2-$\sigma$ confidence level, consistent with high reprocessing efficiency, but poor lateral heat redistribution in the donor's atmosphere.

We surveyed the 25 Ori association for direct-imaging companions. This association has an age of only few million years. Among other targets, we observed CVSO 30, which has recently been identified as the first T Tauri star found to host a transiting planet candidate. We report on photometric and spectroscopic high-contrast observations with the Very Large Telescope, the Keck telescopes, and the Calar Alto observatory. They reveal a directly imaged planet candidate close to the young M3 star CVSO 30. The JHK-band photometry of the newly identified candidate is at better than 1 sigma consistent with late-type giants, early-T and early-M dwarfs, and free-floating planets. Other hypotheses such as galaxies can be excluded at more than 3.5 sigma. A lucky imaging z' photometric detection limit z'= 20.5 mag excludes early-M dwarfs and results in less than 10 MJup for CVSO 30 c if bound. We present spectroscopic observations of the wide companion that imply that the only remaining explanation for the object is that it is the first very young (< 10 Myr) L-T-type planet bound to a star, meaning that it appears bluer than expected as a result of a decreasing cloud opacity at low effective temperatures. Only a planetary spectral model is consistent with the spectroscopy, and we deduce a best-fit mass of 4-5 Jupiter masses (total range 0.6-10.2 Jupiter masses). This means that CVSO 30 is the first system in which both a close-in and a wide planet candidate are found to have a common host star. The orbits of the two possible planets could not be more different: they have orbital periods of 10.76 hours and about 27000 years. The two orbits may have formed during a mutual catastrophic event of planet-planet scattering.

The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.

One of the necessary parameters needed for the interpretation of the light curves of transiting exoplanets or eclipsing binaries, as well as interferometric measurements of a star or microlensing events is how the intensity and polarization of a light change from the center to the limb. Scattering and absorption processes in stellar atmosphere affect both the center-to limb variation of intensity (CLVI) and polarization (CLVP). In this paper, we present a study of the CLVI and CLVP in continuum spectra considering different contributions of scattering and absorption opacity for different spectral type stars with spherical atmospheres. We solve the polarized radiative transfer equation in the presence of continuum scattering, considering spherical stellar model atmospheres. We developed two independent codes based on Feautrier and short characteristics methods to cross-check our results. We calculate the CLVI and CLVP in continuum for the Phoenix grid of spherical stellar model atmospheres for a range of $T_{eff} = 4000 - 7000 \rm K$, $\log g = 1.0 - 5.5$ and $\lambda = 4000 - 7000 \rm \AA$, which are tabulated and available at the CDS. For sub-giant and dwarf stars ($\log g = 3.0 - 4.5$), lower $\log g$ and lower $T_{eff}$ of a star lead to higher limb polarization of the star. For giant and supergiant stars ($\log g = 1.0 - 2.5$), the highest effective temperature yields the largest polarization. By decreasing of the $T_{eff}$ of a star down to $4500 - 5500 \rm K$ (depending on $\log g$) the limb polarization decreases and reaches a local minimum. It increases again down to $T_{eff}$ of $4000 \rm K$. For the most compact dwarf stars ($\log g = 5.0 - 5.5$) the limb polarization degree shows a maximum for models with $T_{eff}$ in the range $4200 - 4600 \rm K$ (depending on $\log g$) and decreases toward higher and lower temperatures.

We investigate the frequency of high carbon-to-oxygen (C/O $= 0.9$) M dwarf stars in the solar neighbourhood. Using synthetic spectra, we find that such M dwarfs would have weaker TiO bands relative to hydride features. Similar weakening has already been detected in M-subdwarf (sdM) stars. By comparing to existing spectroscopic surveys of nearby stars, we show that less than one percent of nearby stars have high carbon-to-oxygen ratios. This limit does not include stars with C/O$=0.9$, [m/H]$>0.3$, and [C/Fe]$>0.1$, which we predict to have low-resolution optical spectra similar to solar metallicity M dwarfs.

AA Dor is a close, totally eclipsing, post common-envelope binary with an sdOB-type primary and an extremely low-mass secondary, located close to the mass limit of stable central hydrogen burning. Within error limits, it may either be a brown dwarf or a late M-type dwarf. We aim to extract the secondary's contribution to the phase-dependent composite spectra. The spectrum and identified lines of the secondary decide on its nature. In January 2014, we measured the phase-dependent spectrum of AA Dor with XSHOOTER over one complete orbital period. Since the secondary's rotation is presumable synchronized with the orbital period, its surface strictly divides into a day and night side. Therefore, we may obtain the spectrum of its cool side during its transit and of its hot, irradiated side close to its occultation. We developed the Virtual Observatory (VO) tool TLISA to search for weak lines of a faint companion in a binary system. We identified 53 spectral lines of the secondary in the ultraviolet-blue, visual, and near-infrared XSHOOTER spectra that are strongest close to its occultation. We identified 57 (20 additional) lines in available UVES (Ultraviolet and Visual Echelle Spectrograph) spectra from 2001. The lines are mostly from C II-III and O II, typical for a low-mass star that is irradiated and heated by the primary. We verified the orbital period of P = 22597.033201 +/- 0.00007 s and determined the orbital velocity Ksec = 232.9 (+16.6 / -6.5) km/s of the secondary. The mass of the secondary is Msec = 0.081 (+0.018 / -0.010) Msun and, hence, it is not possible to reliably determine a brown dwarf or an M-type dwarf nature. Although we identified many emission lines of the secondary's irradiated surface, the resolution and signal-to-noise ratio of our UVES and XSHOOTER spectra are not good enough to extract a good spectrum of the secondary's nonirradiated hemisphere.

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.

We present the atmospheric structure and the fundamental parameters of three red supergiants, increasing the sample of RSGs observed by near-infrared spectro-interferometry. Additionally, we test possible mechanisms that may explain the large observed atmospheric extensions of RSGs. We carried out spectro-interferometric observations of 3 RSGs in the near-infrared K-band with the VLTI/AMBER instrument at medium spectral resolution. To comprehend the extended atmospheres, we compared our observational results to predictions by available hydrostatic PHOENIX, available 3-D convection, and new 1-D self-excited pulsation models of RSGs. Our near-infrared flux spectra are well reproduced by the PHOENIX model atmospheres. The continuum visibility values are consistent with a limb-darkened disk as predicted by the PHOENIX models, allowing us to determine the angular diameter and the fundamental parameters of our sources. Nonetheless, in the case of V602 Car and HD 95686, the PHOENIX model visibilities do not predict the large observed extensions of molecular layers, most remarkably in the CO bands. Likewise, the 3-D convection models and the 1-D pulsation models with typical parameters of RSGs lead to compact atmospheric structures as well, which are similar to the structure of the hydrostatic PHOENIX models. They can also not explain the observed decreases in the visibilities and thus the large atmospheric molecular extensions. The full sample of our RSGs indicates increasing observed atmospheric extensions with increasing luminosity and decreasing surface gravity, and no correlation with effective temperature or variability amplitude, which supports a scenario of radiative acceleration on Doppler-shifted molecular lines.

The main goal of this research is to determine the angular size and the atmospheric structures of cool giant stars and to compare them with hydrostatic stellar model atmospheres, to estimate the fundamental parameters, and to obtain a better understanding of the circumstellar environment. We conducted spectro-interferometric observations of epsilon Oct, beta Peg, NU Pav, and psi Peg in the near-infrared K band (2.13-2.47 microm), and gamma Hya (1.9-2.47 microm) with the VLTI/AMBER instrument at medium spectral resolution. To obtain the fundamental parameters, we compared our data with hydrostatic atmosphere models (PHOENIX). We estimated the Rosseland angular diameters of epsilon Oct, beta Peg, NU Pav, psi Peg, and gamma Hya. Together with distances and bolometric fluxes, we estimated radii, effective temperatures, and luminosities of our targets. In the beta Peg visibility, we observed a molecular layer of CO with a size similar to that modeled with PHOENIX. However, there is an additional slope in absorption starting around 2.3 microm. This slope is possibly due to a shell of water that is not modeled with PHOENIX. The visibility of psi Peg shows a low increase in the CO bands, compatible with the modeling of the PHOENIX model. The visibility data of the other sources show no increase in molecular bands. The spectra and visibilities predicted by the PHOENIX atmospheres agree with the spectra and the visibilities observed in our stars (except for beta Peg). This indicates that the opacity of the molecular bands is adequately included in the model, and the atmospheres of our targets have an extension similar to the modeled atmospheres. The atmosphere of beta Peg is more extended than that predicted by the model. The targets are located close to the red limits of the evolutionary tracks of the STAREVOL model.

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.

In 2010 a sub-stellar companion to the solar analog pre-main sequence star PZ Tel and member of the about 12 Myr old Beta Pic moving group was found by high-contrast direct imaging independently by two teams. In order to determine the basic parameters of this companion more precisely and independent of evolutionary models, hence age independent, we obtained follow-up spectroscopic observations of primary and companion. We use the Spectrograph for INtegral Field Observations in the Near Infrared (SINFONI) at the Very Large Telescope Unit 4/YEPUN of ESO's Paranal Observatory in H+K band and process the data using the spectral deconvolution technique. The resulting spectrum of the companion is then compared to a grid of Drift-Phoenix synthetic model spectra, a combination of a general-purpose model atmosphere code with a non-equilibrium, stationary cloud and dust model, using a chi^2 minimization analysis. We find a best fitting spectral type of G6.5 for PZ Tel A. The extracted spectrum of the sub-stellar companion, at a spatial position compatible with earlier orbit estimates, yields a temperature Teff= 2500 +138-115 K, a visual extinction A_V= 0.53 +0.84-0.53 mag, a surface gravity of log g= 3.50 +0.51-0.30 dex, and a metallicity at the edge of the grid of [M/H]= 0.30 -0.30 dex. We derive a luminosity of log(Lbol/Lsun)= -2.66 +0.06-0.08, a radius of R= 2.42 +0.28-0.34 R_Jup and a mass of M= 7.5 +16.9-4.3 M_Jup for the PZ Tel companion, being consistent with most earlier estimates using photometry alone. Combining our results with evolutionary models, we find a best fitting mass of about 21 Jupiter masses at an age corresponding to the recently determined lithium depletion age of 7 +4-2 Myr. Hence, the PZ Tel companion is most likely a wide brown dwarf companion in the 12 +8-4 Myr old Beta Pic moving group.

We present the atmospheric structure and the fundamental properties of the red supergiants (RSGs) AH Sco, UY Sct, and KW Sgr based on VLTI/AMBER observations. We carried out spectro-interferometric observations of AH Sco, UY Sct, and KW Sgr in the near-infrared K band with the VLTI/AMBER instrument, and compared the data to a new grid of hydrostatic PHOENIX model atmospheres. In our visibility data, we observe molecular layers of water and CO in extended atmospheres. For a uniform disk modeling, we observe size increases at the water band of 10% to 25% and at the CO bandheads of 20%-35% with respect to the near continuum bandpass. The PHOENIX atmosphere models predict the spectra and the continuum visibility values, but cannot reproduce the large extensions of the molecular layers. This indicates that the opacities of the molecular bands are included, but that the model atmospheres are too compact compared to the observations. The observed extended layers may be levitated by processes such as pulsation or convection, which are not included in the hydrostatic atmospheric models. Comparing the continuum visibility values to PHOENIX models, we estimate the Rosseland-mean photospheric angular diameters. Together with the distance and the spectro-photometry, we calculate radii and effective temperatures. The location of the targets in the HR-diagram is confirmed to be close to, and possibly slightly to the right of, the Hayashi limit of recent evolutionary tracks corresponding to masses between about 20 Msun and 40 Msun

We present a new library of high-resolution synthetic spectra based on the stellar atmosphere code PHOENIX that can be used for a wide range of applications of spectral analysis and stellar parameter synthesis. The spherical mode of PHOENIX was used to create model atmospheres and to derive detailed synthetic stellar spectra from them. We present a new self-consistent way of describing micro-turbulence for our model atmospheres. The synthetic spectra cover the wavelength range from 500AA to 50.000AA with resolutions of R=500.000 in the optical and near IR, R=100.000 in the IR and a step size of 0.1AA in the UV. The parameter space covers 2.300K<=Teff<=12.000K, 0.0<=log(g)<=+6.0, -4.0<=[Fe/H]<=+1.0, and -0.2<=[alpha/Fe]<=+1.2. The library is a work in progress and we expect to extend it up to Teff=25.000 K.

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.

We investigate the atmospheric structure and fundamental properties of the red supergiant VY CMa. We obtained near-infrared spectro-interferometric observations of VY CMa with spectral resolutions of 35 and 1500 using the AMBER instrument at the VLTI. The visibility data indicate the presence of molecular layers of water vapor and CO in the extended atmosphere with an asymmetric morphology. The uniform disk diameter in the water band around 2.0 mu is increased by \sim20% compared to the near-continuum bandpass at 2.20-2.25 mu and in the CO band at 2.3-2.5 mu it is increased by up to \sim50%. The closure phases indicate relatively small deviations from point symmetry close to the photospheric layer, and stronger deviations in the extended H2O and CO layers. Making use of the high spatial and spectral resolution, a near-continuum bandpass can be isolated from contamination by molecular and dusty layers, and the Rosseland-mean photospheric angular diameter is estimated to 11.3 +/- 0.3 mas based on a PHOENIX atmosphere model. Together with recent high-precision estimates of the distance and spectro-photometry, this estimate corresponds to a radius of 1420 +/- 120 Rsun and an effective temperature of 3490 +/- 90 K. VY CMa exhibits asymmetric, possibly clumpy, atmospheric layers of H2O and CO, which are not co-spatial, within a larger elongated dusty envelope. Our revised fundamental parameters put VY CMa close to the Hayashi limit of recent evolutionary tracks of initial mass 25 Msun with rotation or 32 Msun without rotation, shortly before evolving blueward in the HR-diagram.

We present a new library of synthetic spectra based on the stellar atmosphere code PHOENIX. It covers the wavelength range from 500{\AA} to 55000{\AA} with a resolution of R=500000 in the optical and near IR, R=100000 in the IR and {\Delta}{\lambda}=0.1{\AA} in the UV. The parameter space covers 2300K<=Teff<=8000K, 0.0<=log(g)<=6.0, -4.0<=[Fe/H]<=+1.0 and -0.3<=[{\alpha}/Fe]<=+0.8. The library is work-in-progress and going to be extended to at least Teff=25000K. We use a new self-consistent way of describing the microturbulence for our model atmospheres. The entire library of synthetic spectra will be available for download. Futhermore we present a method for fitting spectra, especially designed to work with the new 2nd generation VLT instrument MUSE. We show that we can determine stellar parameters (Teff, log(g), [Fe/H] and [{\alpha}/Fe]) and even single element abundances.

With a uniform VLT SINFONI data set of nine targets, we have developed an empirical grid of J,H,K spectra of the atmospheres of objects estimated to have very low substellar masses of \sim5-20 MJup and young ages of \sim1-50 Myr. Most of the targets are companions, objects which are especially valuable for comparison with atmosphere and evolutionary models, as they present rare cases in which the age is accurately known from the primary. Based on the sample youth, all objects are expected to have low surface gravity, and this study investigates the critical early phases of the evolution of substellar objects. The spectra are compared with grids of five different theoretical atmosphere models. This analysis represents the first systematic model comparison with infrared spectra of young brown dwarfs. The fits to the full JHK spectra of each object result in a range of best fit effective temperatures of +/-150-300K whether or not the full model grid or a subset restricted to lower log(g) values is used. This effective temperature range is significantly larger than the uncertainty typically assigned when using a single model grid. Fits to a single wavelength band can vary by up to 1000K using the different models. Since the overall shape of these spectra is governed more by the temperature than surface gravity, unconstrained model fits did not find matches with low surface gravity or a trend in log(g) with age. This suggests that empirical comparison with spectra of unambiguously young objects targets (such as these SINFONI data) may be the most reliable method to search for indications of low surface gravity and youth. For two targets, the SINFONI data are a second epoch and the data show no variations in morphology over time. The analysis of two other targets, AB Pic B and CT Cha B, suggests that these objects may have lower temperatures, and consequently lower masses, than previously estimated.

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.

The existence of predominantly cold non-baryonic dark matter is unambiguously demonstrated by several observations (e.g., structure formation, big bang nucleosynthesis, gravitational lensing, and rotational curves of spiral galaxies). A candidate well motivated by particle physics is a weakly interacting massive particle (WIMP). Self-annihilating WIMPs would affect the stellar evolution especially in the early universe. Stars powered by self-annihilating WIMP dark matter should possess different properties compared with standard stars. While a direct detection of such dark matter powered stars seems very challenging, their cumulative emission might leave an imprint in the diffuse metagalactic radiation fields, in particular in the mid-infrared part of the electromagnetic spectrum. In this work the possible contributions of dark matter powered stars (dark stars; DSs) to the extragalactic background light (EBL) are calculated. It is shown that existing data and limits of the EBL intensity can already be used to rule out some DS parameter sets.

We present a grid of atmospheric models and synthetic spectral energy distributions (SEDs) for late-type dwarfs and giants of solar and 1/3 solar metallicity with many opacity sources computed in self-consistent Non-Local Thermodynamic Equilibrium (NLTE), and compare them to the LTE grid of Short & Hauschildt (2010) (Paper I). We describe, for the first time, how the NLTE treatment affects the thermal equilibrium of the atmospheric structure (T(tau) relation) and the SED as a finely sampled function of Teff, log g, and [A/H] among solar metallicity and mildly metal poor red giants. We compare the computed SEDs to the library of observed spectrophotometry described in Paper I across the entire visible band, and in the blue and red regions of the spectrum separately. We find that for the giants of both metallicities, the NLTE models yield best fit Teff values that are ~30 to 90 K lower than those provided by LTE models, while providing greater consistency between \log g values, and, for Arcturus, Teff values, fitted separately to the blue and red spectral regions. There is marginal evidence that NLTE models give more consistent best fit Teff values between the red and blue bands for earlier spectral classes among the solar metallicity GK giants than they do for the later classes, but no model fits the blue band spectrum well for any class. For the two dwarf spectral classes that we are able to study, the effect of NLTE on derived parameters is less significant.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

We check how the change in surface conditions between the Sun and red giant branch stars changes the characteristic surface convection length scale to be used in models. We investigate the question in the case of the mixing length theory and of the phenomenology of full spectrum of turbulence. For the observational part, we rely on independent measurements of effective temperatures and interferometric radii of nearby red giants. We find that the local red giant branch cannot be explained taking into account the solar calibrated convective length scale.

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.

Context. Radial velocity monitoring of very cool dwarfs such as late M- and hot L-dwarfs has become a promising tool to search for rocky planets as well as to follow-up planetary candidates around dwarfs found by transit surveys. These stars are faint at optical wavelengths, as their spectral flux distribution peaks at near-infrared (NIR) wavelengths. For this reason, it is desirable to measure the radial velocities in this wavelength regime. However, in the NIR there are only very few medium- and high-resolution spectrographs available which are mounted at large telescopes. In the near future, high-resolution spectrographs for the NIR will be built, which will allow us to search for rocky planets around cool M-dwarfs and L-dwarfs from radial velocities monitoring. Methods. Stellar atmosphere synthetic models for an M- and an L-dwarf with temperatures of 2200 K and 1800 K, respectively, and a theoretical spectrum of the Earth's transmission in the spectral range from 0.9 to 2.5 \mum are used. We simulate a series of Doppler-shifted spectra observed with different resolving powers and signal-to-noise ratios, and for different rotational broadenings of the dwarf. For different combinations of the input parameters we recover the radial velocity by means of cross-correlation with a high signal-to-noise ratio template and determine the associate uncertainties. Results. The highest precision in radial velocity measurements for the cool M-dwarf is found in the Y band around 1.0 \mum, while for the L-dwarf it is determined in the J band around 1.25 \mum. We note that synthetic models may lack of some faint absorption features or underestimate their abundances. Conversely, some instrumental/calibration aspects that are not taken into account in our estimations would rise the uncertainties.

We present the discovery of a substellar companion on a wide orbit around the ~2.5 Msun star HIP 78530, which is a member of the 5 Myr-old Upper Scorpius association. We have obtained follow-up imaging over two years and show that the companion and primary share common proper motion. We have also obtained JHK spectroscopy of the companion and confirm its low surface gravity, in accordance with the young age of the system. A comparison with Drift-Phoenix synthetic spectra indicates an effective temperature of 2800+/-200 K and a comparison with template spectra of young and old dwarfs indicates a spectral type of M8+/-1. The mass of the companion is estimated to be 19-26 Mjup based on its bolometric luminosity and the predictions of evolutionary models. The angular separation of the companion is 4.5", which at the distance of the primary star, 156.7 pc, corresponds to a projected separation of ~710 AU. This companion features one of the lowest mass ratios (~0.009) of any known companion at separations greater than 100 AU.

The phenomenological models of convection use characteristic length scales they do not determine but that are chosen to fit solar or stellar observations. We investigate if changes of these length scales are required between the Sun and low mass stars on the red giant branch (RGB). The question is addressed jointly in the frameworks of the mixing length theory and of the full spectrum of turbulence model. For both models, the convective length scale is assumed to be a fixed fraction of the local pressure scale height. We use constraints coming from the observed effective temperatures and linear radii independently. We rely on a sample of 38 nearby giants and subgiants for which surface temperatures and luminosities are known accurately and the radii are determined through interferometry to better than 10%. For the few cases where the stellar masses were determined by asteroseismological measurements, we computed dedicated models. First we calibrate the solar models. Then, with the same physics, we compute RGB models for masses between 0.9 Mo and 2.5 Mo and metallicities ranging from $\rm [Fe/H]=-0.34$ to solar. The evolution is followed up to 1000 Lo. A special attention is given to the opacities and to the non grey atmosphere models used as boundary conditions for which the model of convection is the same as in the interior. We find that for both the mixing length theory and the full spectrum of turbulence model the characteristic solar length scale for convection has to be slightly reduced to fit the lower edge of the observed RGB. The corresponding models also better match the expected mass distribution on the RGB and are in better agreement to the seismic constraints. These results are robust whether effective temperatures determined spectroscopically or radii determined interferometrically are used.

We present a systematic study of the physical properties of late-M dwarfs based on high-quality dynamical mass measurements and near-infrared (NIR) spectra. We use astrometry from Keck NGS and LGS AO imaging to determine orbits for late-M binaries. We find that LP 349-25 (M7.5+M8) is a pair of young brown dwarfs (Mtot = 0.120 Msun) for which Lyon and Tucson evolutionary models jointly predict an age of 140+/-30 Myr. This is consistent with the age of the Pleiades, but at least LP 349-25A defies the empirical Pleiades lithium depletion boundary, implying that the system is in fact older and that evolutionary models underpredict the component luminosities. We find that LHS 1901AB (M6.5+M6.5) is a pair of very low-mass stars (Mtot = 0.194 Msun) with model-derived ages consistent with limits from its lack of activity (> 6 Gyr). Our improved orbit for Gl 569Bab (M8.5+M9) results in a higher mass for this binary (Mtot = 0.140 Msun) compared to previous work (0.125 Msun). We use these masses along with our published results for 2MASS J2206-2047AB (M8+M8) to test four sets of ultracool model atmospheres currently in use. Fitting these models to our NIR integrated-light spectra provides temperature estimates warmer by ~250 K than those derived independently from Dusty evolutionary models given the measured masses and luminosities. We propose that model atmospheres are more likely to be the source of this discrepancy, as it would be difficult to explain a uniform temperature offset over such a wide range of masses, ages, and activity levels in the context of evolutionary models. Our results contrast those of Konopacky et al. as we find an opposite and smaller mass discrepancy from what they report when we adopt their model-testing approach since our Teff estimates from fitting spectra are ~650 K higher than from their fitting of broadband photometry alone.

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 a grid of LTE atmospheric models and synthetic spectra that cover the spectral class range from mid-G to mid-K, and luminosity classes from V to III, that is dense in Teff sampling (Delta Teff=62.5 K), for stars of solar metallicity and moderately metal poor scaled solar abundance ([A/H]=0.0 and -0.5). All models have been computed with two choices of atomic line list: a) the "big" line lists of Kurucz (1992) that best reproduce the broad-band solar blue and near UV flux level, and b) the "small" lists of Kurucz & Peytremann (1975) that provide the best fit to the high resolution solar blue and near-UV spectrum. We compare our model SEDs to a sample of stars carefully selected from the large catalog of uniformly re-calibrated spectrophotometry of Burnashev (1985) with the goal of determining how the quality of fit varies with stellar parameters, especially in the historically troublesome blue and near-UV bands. We confirm that our models computed with the "big" line list recover the derived Teff values of the PHOENIX NextGen grid, but find that the models computed with the "small" line list provide greater internal self-consistency among different spectral bands, and closer agreement with the empirical Teff scale of Ramirez & Melendez (2005), but not to the interferometrically derived Teff values of Baines et al. (2010). We find no evidence that the near UV band discrepancy between models and observations for Arcturus (alpha Boo) reported by Short & Hauschildt (2003 and 2009) is pervasive, and that Arcturus may be peculiar in this regard.

We present results from the earliest observations of DEBRIS, a Herschel Key Programme to conduct a volume- and flux-limited survey for debris discs in A-type through M-type stars. PACS images (from chop/nod or scan-mode observations) at 100 and 160 micron are presented toward two A-type stars and one F-type star: beta Leo, beta UMa and eta Corvi. All three stars are known disc hosts. Herschel spatially resolves the dust emission around all three stars (marginally, in the case of beta UMa), providing new information about discs as close as 11 pc with sizes comparable to that of the Solar System. We have combined these data with existing flux density measurements of the discs to refine the SEDs and derive estimates of the fractional luminosities, temperatures and radii of the discs.

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.

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.

Substellar objects have extremely long life-spans. The cosmological consequence for older objects are low abundances of heavy elements, which results in a wide distribution of objects over metallicity, hence over age. Within their cool atmosphere, dust clouds become a dominant feature, affecting the opacity and the remaining gas phase abundance of heavy elements. We investigate the influence of the stellar metallicity on the dust formation in substellar atmospheres and on the dust cloud structure and its feedback on the atmosphere. We utilize numerical simulations in which we solve a set of moment equations in order to determine the quasi-static dust cloud structure (DRIFT). These equations model the nucleation, the kinetic growth of composite particles, their evaporation and the gravitational settling as a stationary dust formation process. Element conservation equations augment this system of equations including the element replenishment by convective overshooting. The integration with an atmosphere code (PHOENIX) allows to determine a consistent (T, p, v_conv)-structure, and, hence, also to calculate synthetic spectra. A grid of DRIFT-PHOENIX model atmospheres was calculated for a wide range of metallicity to allow for a systematic study of atmospheric cloud structures throughout the evolution of the universe. We find dust clouds in even the most metal-poor ([M/H]=-6.0) atmosphere of brown dwarfs. Only the most massive among the youngest brown dwarfs and giant gas planets can resist dust formation. For very low heavy element abundances, a temperature inversion develops which has a drastic impact on the dust cloud structure. We further show that the dust-to-gas ratio does not scale linearly with the object's [M/H] for a given effective temperature.

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.

Gas and dust disks are common objects in the universe and can be found around various objects, e.g. young stars, cataclysmic variables, active galactic nuclei, or white dwarfs. The light that we receive from disks provides us with clues about their composition, temperature, and density. In order to better understand the physical and chemical dynamics of these disks, self-consistent radiative transfer simulations are inevitable. Therefore, we have developed a 1+1D radiative transfer code as an extension to the well-established model atmosphere code \verb!PHOENIX!. We will show the potential of the application of our code to model the spectra of white dwarf debris disks.

We determine the effective temperature, surface gravity and projected rotational velocity of nine T dwarfs from the comparison of high-resolution near-infrared spectra and synthetic models, and estimate the mass and age of the objects from state-of-the-art models. We use the AMES-COND cloudless solar metallicity models provided by the PHOENIX code to match the spectra of the T dwarfs observed with the NIR high-resolution spectrograph NIRSPEC using ten echelle orders to cover part of the J band from 1.165 to 1.323 microns with R~20,000. The projected rotational velocity, Teff and logg of the objects are determined based on the minimum RMS of the differences between the modelled and observed relative fluxes. The modelled spectra reproduce quite well the observed features for most of the T dwarfs, with Teff in the range of 922-1009 K, and logg between 4.3 and 5 (cm s^-2). Our results support the assumption of a dust free atmosphere for T dwarfs later than T5, where dust grains form and then gravitationally sediment into the low atmosphere. The modelled spectra do not accurately mimic some individual very strong lines like the KI doublet at 1.2436 and 1.2525 microns. Our modelled spectra does not match well the observed spectra of the two T dwarfs with earlier spectral types, namely SDSSp J125453.90-012247.4 (T2) and 2MASS J05591914-1404488 (T4.5), which is likely due to the presence of condensate clouds that are not incorporated in the models used here. By comparing our results and their uncertainties to evolutionary models, we estimate masses in the interval ~5-75 MJ for T dwarfs later than T5, which are in good agreement with those found in the literature. We found apparent young ages that are typically between 0.1 and a few Gyr for the same T dwarfs, which is consistent with recent kinematical studies.

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 present a new code for the calculation of the 1D structure and synthetic spectra of accretion disks. The code is an extension of the general purpose stellar atmosphere code PHOENIX and is therefore capable of including extensive lists of atomic and molecular lines as well as dust in the calculations. We assume that the average viscosity can be represented by a critical Reynolds number in a geometrically thin disk and solve the structure and radiative transfer equations for a number of disk rings in the vertical direction. The combination of these rings provides the total disk structure and spectrum. Since the warm inner regions of protoplanetary disks show a rich molecular spectrum, they are well suited for a spectral analysis with our models. In this paper we test our code by comparing our models with high-resolution VLT CRIRES spectra of the T Tauri star GQ Lup.

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.

Red optical and near-infrared spectroscopy are presented for SDSS J125637.13-022452.4, one of only four L subdwarfs reported to date. These data confirm the low-temperature, metal-poor nature of this source, as indicated by prominent metal-hydride bands, alkali lines, and collision-induced H2 absorption. The optical and near-infrared spectra of SDSS J1256-0224 are similar to those of the sdL4 2MASS J16262034+3925190, and we derive a classification of sdL3.5 based on the preliminary scheme of Burgasser, Cruz, & Kirkpatrick. The kinematics of SDSS J1256-0224 are consistent with membership in the Galactic inner halo, with estimated $UVW$ space velocities indicating a slightly prograde, eccentric and inclined Galactic orbit (3.5 <~ R <~ 11 kpc; |Zmax| = 7.5 kpc). Comparison to synthetic spectra computed with the Phoenix code, including the recent implementation of kinetic condensate formation (Drift-Phoenix), indicate Teff ~ 2100-2500 K and [M/H] ~ -1.5 to -1.0 for logg ~ 5.0-5.5 (cgs), although there are clear discrepancies between model and observed spectra particularly in the red optical region. The stronger metal-oxide bands present in the Drift-Phoenix model spectra, a result of phase-non-equilibrium abundances of grain species, appears to contradict prior suggestions that grain formation is inhibited in metal-poor atmospheres; conclusive statements on the metallicity dependence of grain formation efficiency are as yet premature. In addition, an apparent shift in the temperature scale of L subdwarfs relative to L dwarfs may obviate the need for modified grain chemistry to explain some of the former's unique spectral characteristics.

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.

We investigate the ability of both LTE and Non-LTE models to fit the near UV band absolute flux distribution and individual spectral line profiles of three standard stars for which high quality spectrophotometry and high resolution spectroscopy are available: The Sun (G2 V), Arcturus (K2 III), and Procyon (F5 IV-V). We investigate 1) the effect of the choice of atomic line list on the ability of NLTE models to fit the near UV band flux level, 2) the amount of a hypothesized continuous thermal absorption extinction source required to allow NLTE models to fit the observations, and 3) the semi-empirical temperature structure required to fit the observations with NLTE models and standard continuous near UV extinction. We find that all models that are computed with high quality atomic line lists predict too much flux in the near UV band for Arcturus, but fit the warmer stars well. The variance among independent measurements of the solar irradiance in the near UV is sufficiently large that we cannot definitely conclude that models predict too much near UV flux, in contrast to other recent results. We surmise that the inadequacy of current atmospheric models of K giants in the near UV band is best addressed by hypothesizing that there is still missing continuous thermal extinction, and that the missing near UV extinction becomes more important with decreasing effective temperature for spectral classes later than early G, suggesting a molecular origin.

The atmospheres of substellar objects contain clouds of oxides, iron, silicates, and other refractory condensates. Water clouds are expected in the coolest objects. The opacity of these `dust' clouds strongly affects both the atmospheric temperature-pressure profile and the emergent flux. Thus any attempt to model the spectra of these atmospheres must incorporate a cloud model. However the diversity of cloud models in atmospheric simulations is large and it is not always clear how the underlying physics of the various models compare. Likewise the observational consequences of different modeling approaches can be masked by other model differences, making objective comparisons challenging. In order to clarify the current state of the modeling approaches, this paper compares five different cloud models in two sets of tests. Test case 1 tests the dust cloud models for a prescribed L, L--T, and T-dwarf atmospheric (temperature T, pressure p, convective velocity vconv)-structures. Test case 2 compares complete model atmosphere results for given (effective temperature Teff, surface gravity log g). All models agree on the global cloud structure but differ in opacity-relevant details like grain size, amount of dust, dust and gas-phase composition. Comparisons of synthetic photometric fluxes translate into an modelling uncertainty in apparent magnitudes for our L-dwarf (T-dwarf) test case of 0.25 < \Delta m < 0.875 (0.1 < \Delta m M 1.375) taking into account the 2MASS, the UKIRT WFCAM, the Spitzer IRAC, and VLT VISIR filters with UKIRT WFCAM being the most challenging for the models. (abr.)

In our ongoing search for close and faint companions around T Tauri stars, we found a very faint (Ks=14.9mag, Ks_0=14.4mag) object, just ~2.67" northwest of the Chamaeleon star-forming region member CT Cha corresponding to a projected separation of ~440AU at 165+/-30 pc. We show that CT Cha A and this faint object form a common proper motion pair from data of the VLT Adaptive Optics (AO) instrument NACO taken in February 2006 and March 2007 and that the companion is by >=4 sigma significance not a stationary background object. Our AO integral field spectroscopy with SINFONI in J, and H+K bands yields a temperature of 2600+/-250K for the companion and an optical extinction of A_V=5.2+/-0.8mag, when compared to spectra calculated from Drift-Phoenix model atmospheres. We demonstrate the validity of the model fits by comparison to several other well-known young sub-stellar objects. Relative flux calibration of the bands was achieved using photometry from the NACO imaging data. We conclude that the CT Cha companion is a very low-mass member of Chamaeleon and very likely a physical companion to CT Cha, as the probability for a by chance alignment is <=0.01. Due to a prominent Pa-Beta emission in the J-band, accretion is probably still ongoing onto the CT Cha companion. From temperature and luminosity (log(Lbol/Lsun)= -2.68+/-0.21), we derive a radius of R=2.20+0.81-0.60 R_Jup. We find a consistent mass of M=17+/-6 MJup for the CT Cha companion from both its luminosity and temperature when placed on evolutionary tracks. Hence, the CT Cha companion is most likely a wide brown dwarf companion or possibly even a planetary mass object.

We estimate cluster ages from lithium depletion in five pre-main-sequence groups found within 100 pc of the Sun: TW Hydrae Association, Eta Chamaeleontis Cluster, Beta Pictoris Moving Group, Tucanae-Horologium Association and AB Doradus Moving Group. We determine surface gravities, effective temperatures and lithium abundances for over 900 spectra through least squares fitting to model-atmosphere spectra. For each group, we compare the dependence of lithium abundance on temperature with isochrones from pre-main-sequence evolutionary tracks to obtain model dependent ages. We find that the Eta Chamaelontis Cluster and the TW Hydrae Association are the youngest, with ages of 12+/-6 Myr and 12+/-8 Myr, respectively, followed by the Beta Pictoris Moving Group at 21+/-9 Myr, the Tucanae-Horologium Association at 27+/-11 Myr, and the AB Doradus Moving Group at an age of at least 45 Myr (where we can only set a lower limit since the models -- unlike real stars -- do not show much lithium depletion beyond this age). Here, the ordering is robust, but the precise ages depend on our choice of both atmospheric and evolutionary models. As a result, while our ages are consistent with estimates based on Hertzsprung-Russell isochrone fitting and dynamical expansion, they are not yet more precise. Our observations do show that with improved models, much stronger constraints should be feasible: the intrinsic uncertainties, as measured from the scatter between measurements from different spectra of the same star, are very low: around 10 K in effective temperature, 0.05 dex in surface gravity, and 0.03 dex in lithium abundance.

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

We aim to understand cloud formation in substellar objects. We combined the non-equilibrium, stationary cloud model of Helling, Woitke & Thi (2008; seed formation, growth, evaporation, gravitational settling, element conservation) with the general-purpose model atmosphere code PHOENIX (radiative transfer, hydrostatic equilibrium, mixing length theory, chemical equilibrium) in order to consistently calculate cloud formation and radiative transfer with their feedback on convection and gas phase depletion. We calculate the complete 1D model atmosphere structure and the chemical details of the cloud layers. The DRIFT-PHOENIX models enable the first stellar atmosphere simulation that is based on the actual cloud formation process. The resulting (T,p) profiles differ considerably from the previous limiting PHOENIX cases DUSTY and COND. A tentative comparison with observations demonstrates that the determination of effective temperatures based on simple cloud models has to be applied with care. Based on our new models, we suggest a mean Teff=1800K for the L-dwarf twin-binary system DENIS J0205-1159 which is up to 500K hotter than suggested in the literature. We show transition spectra for gas-giant planets which form dust clouds in their atmospheres and evaluate photometric fluxes for a WASP-1 type system.

We aim to understand the correlation between cloud formation and alkali line formation in substellar atmospheres.We perform line profile calculations for Na I and K I based on the coupling of our kinetic model for the formation and composition of dust grains with 1D radiative transfer calculations in atmosphere models for brown dwarfs and giant gas planets. The Na I and K I line profiles sensibly depend on the way clouds are treated in substellar atmosphere simulations. The kinetic dust formation model results in the highest pseudo-continuum compared to the limiting cases.

A test case comparison is presented for different dust cloud model approaches applied in brown dwarfs and giant gas planets. We aim to achieve more transparency in evaluating the uncertainty inherent to theoretical modelling. We show in how far model results for characteristic dust quantities vary due to different assumptions. We also demonstrate differences in the spectral energy distributions resulting from our individual cloud modelling in 1D substellar atmosphere simulations

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 PHOENIX atmosphere models for metal-rich cool white dwarfs using improved line shapes for the Na I and K I resonance doublets. Profiles for collisional broadening due to H2 and He based on the adiabatic representation show strong deviations from Van der Waals interaction at short distances. Comparison with observed spectra that show extremely broadened Na I lines indicates that a He-rich atmospheric composition is required to explain the line strengths and spectral energy distributions. Our current synthetic spectra, using an expansion in powers of density to the third order optimised for brown dwarf atmosphere conditions, significantly underestimate the observed absorption in the far wings, even predicting smaller total line strength than a Lorentzian profile. This is due to the handling of multiple perturber interactions becoming inadequate for the extreme densities of the coolest white dwarfs. The density expansion would have to be extended at least to the 7th order for an accurate treatment of such conditions and might break down altogether in the densest objects. The results of a direct calculation of the unified profile should therefore be used for model atmospheres of cool metal-rich white dwarfs. Qualitative comparison of the full adiabatic profile to the spectrum of WD2356-209 indicates good agreement with the observed line shape. Observations of the coolest white dwarfs may therefore serve as a laboratory for testing the physics of the deeper atmospheres and interiors of brown dwarfs and giant planets.

We present a UVES/VLT high resolution atlas of three L dwarfs and one T dwarf system, spectral classes at which most of the objects are brown dwarfs. Our atlas covers the optical region from H$\alpha$ up to the near infrared at 1 $\mu$m. We present spectral details of ultra-cool atmospheres at very high resolution ($R \sim 33 000$) and compare the spectra to model calculations. Our comparison shows that molecular features from VO and CaH, and atomic features from Cs and Rb are reasonably well fit by current models. On the other hand, features due to TiO, CrH, and water, and atomic Na and K reveal large discrepancies between model calculations and our observations.

The phase of Super-Soft-Source (SSS) emission of the sixth recorded outburst of the recurrent nova RS Oph was observed twice with Chandra and once with XMM-Newton. The observations were taken on days 39.7, 54.0, and 66.9 after outburst. We confirm a 35-sec period on day 54.0 and found that it originates from the SSS emission and not from the shock. We discus the bound-free absorption by neutral elements in the line of sight, resonance absorption lines plus self-absorbed emission line components, collisionally excited emission lines from the shock, He-like intersystem lines, and spectral changes during an episode of high-amplitude variability. We find a decrease of the oxygen K-shell absorption edge that can be explained by photoionization of oxygen. The absorption component has average velocities of -1286+-267 km/s on day 39.7 and of -771+-65 km/s on day 66.9. The wavelengths of the emission line components are consistent with their rest wavelengths as confirmed by measurements of non-self absorbed He-like intersystem lines. We have evidence that these lines originate from the shock rather than the outer layers of the outflow and may be photoexcited in addition to collisional excitations. We found collisionally excited emission lines that are fading at wavelengths shorter than 15A that originate from the radiatively cooling shock. On day 39.5 we find a systematic blue shift of -526+-114 km/s from these lines. We found anomalous He-like f/i ratios which indicates either high densities or significant UV radiation near the plasma where the emission lines are formed. During the phase of strong variability the spectral hardness light curve overlies the total light curve when shifted by 1000sec. This can be explained by photoionization of neutral oxygen in the line of sight if the densities of order 10^{10}-10^{11} cm^{-3}.

We report optical photometry and optical through mid-infrared spectroscopy of the classical nova V1186 Sco. This slowly developing nova had an complex light curve with multiple secondary peaks similar to those seen in PW Vul. The time to decline 2 magnitudes, t$_2$, was 20 days but the erratic nature of the light curve makes determination of intrinsic properties based on the decline time (e.g., luminosity) problematic, and the often cited MMRD relationship of Della Valle and Livio (1995) fails to yield a plausible distance. Spectra covering 0.35 to 35 $\mu$m were obtained in two separate epochs during the first year of outburst. The first set of spectra, taken about 2 months after visible maximum, are typical of a CO-type nova with narrow line emission from \ion{H}{1}, \ion{Fe}{2}, \ion{O}{1} and \ion{He}{1}. Later data, obtained between 260 and 380 days after maximum, reveal an emerging nebular spectrum. \textit{Spitzer} spectra show weakening hydrogen recombination emission with the emergence of [\ion{Ne}{2}] (12.81 $\mu$m) as the strongest line. Strong emission from [\ion{Ne}{3}] (15.56 $\mu$m) is also detected. Photoionization models with low effective temperature sources and only marginal neon enhancement (Ne $\sim$ 1.3 Ne$_{\odot}$) are consistent with these IR fine-structure neon lines indicating that V1186 Sco did not occur on a ONeMg white dwarf. In contrast, the slow and erratic light curve evolution, spectral development, and photoionization analysis of the ejecta imply the outburst occurred on a low mass CO white dwarf. We note that this is the first time strong [\ion{Ne}{2}] lines have been detected so early in the outburst of a CO nova and suggests that the presence of mid-infrared neon lines is not directly indicative of a ONeMg nova event.

We provide a grid of PHOENIX spectra of red giant and supergiant stars, that extend through optical and near-IR wavelengths. For the first time, models are also provided with modified surface abundances of C, N and O, as a step towards accounting for the changes known to occur due to convective dredge-up (and to be enhanced in the case of rotation). The aims are (i) to assess how well current models reproduce observed spectra, (ii) to quantify the effects of the abundance changes on the spectra, and (iii) to determine how these changes affect estimates of fundamental stellar parameters. Observed giant star spectra can be fitted very well at solar metallicity down to about 3400K. Modified surface abundances are preferred in only a minority of cases for luminosity class II, possibly indicating mixing in excess of standard first dredge-up. Supergiant stars show a larger variety of near-IR spectra, and good fits are currently obtained for about one third of the observations only. Modified surface abundances help reproducing strong CN bands, but do not suffice to resolve all the difficulties. The effect of the abundance changes on the estimated Teff depends on the wavelength range of observation and can amount several 100K. Reasons for the remaining discrepancies are discussed.

We report on recent progress in the modelling of the near-IR spectra of young stellar populations, i.e. populations in which red supergiants (RSGs) are dominant. First, we discuss the determination of fundamental parameters of RSGs using fits to their near-IR spectra with new PHOENIX model spectra; RSG-specific surface abundances are accounted for and effects of the microturbulence parameter are explored. New population synthesis predictions are then described and, as an example, it is shown that the spectra of young star clusters in M82 can be reproduced very well from 0.5 to 2.4 micrometers. We warn of remaining uncertainties in cluster ages.

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.

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.

The first substellar companion of possibly planetary mass around a normal star, GQ Lup, has been directly imaged (Neuhaeuser et al., 2005). Besides the unknown formation history, the mass of such an object is a criterion to decide about its true nature. We aim to determine the physical properties of the GQ Lup companion - effective temperate and surface gravity, and thus its mass independently from evolutionary models. We use the adaptive optics integral-field spectrograph SINFONI at the VLT for near-infrared spectroscopy from 1.1 to 2.5 um with a resolution of R = 2500--4000. We compare these spectra with synthetic atmospheric models (GAIA v2.0 cond). From the complete set of spectra we find a consistent effective temperature of Teff = 2600 +/- 100 K and surface gravity of log{g} = 3.7 +/- 0.5. Combined with a slightly revised luminosity of log(L/L_(\sun)) = -2.25 +/- 0.24 for the companion, we determine a radius of R = 3.50 (+1.50/-1.03) Jupiterradii and thus a mass of ~ 25 Jupitermasses. Due to the large uncertainty of the surface gravity, the mass could range from 4 to 155 Jupitermasses. By comparing the paramaters of the companion of GQ Lup to the ones of 2MASS J05352184-0546085, published by Stassun et al. (2006), we conclude that the companion to GQ Lup A has a mass lower than 36 Jupitermasses.

We determine stellar parameters for the M dwarf GJ 436 that hosts a Neptune-mass planet. We employ primarily spectral modeling at low and high resolution, examining the agreement between model and observed optical spectra of five comparison stars of type, M0-M3. Modeling high resolution optical spectra suffers from uncertainties in TiO transitions, affecting the predicted strengths of both atomic and molecular lines in M dwarfs. The determination of Teff, gravity, and metallicity from optical spectra remains at ~10%. As molecules provide opacity both in lines and as an effective continuum, determing molecular transition parameters remains a challenge facing models such as the PHOENIX series, best verified with high resolution and spectrophotometric spectra. Our analysis of GJ 436 yields an effective temperature of Teff = 3350 +/- 300 K and a mass of 0.44 Msun. New Doppler measurements for GJ 436 with a precision of 3 m/s taken during 6 years improve the Keplerian model of the planet, giving a minimum mass, M sin i = 0.0713 Mjup = 22.6 Mearth, period, P = 2.6439 d, and e = 0.16 +/- 0.02. The noncircular orbit contrasts with the tidally circularized orbits of all close-in exoplanets, implying either ongoing pumping of eccentricity by a more distant companion, or a higher Q value for this low-mass planet. The velocities indeed reveal a long term trend, indicating a possible distant companion.

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.

We constructed hydrodynamical model atmospheres for mid M-type main-, as well as pre-main-sequence (PMS) objects. Despite the complex chemistry encountered in these cool atmospheres a reasonably accurate representation of the radiative transfer is possible, even in the context of time-dependent and three-dimensional models. The models provide detailed information about the morphology of M-type granulation and statistical properties of the convective surface flows. In particular, we determined the efficiency of the convective energy transport, and the efficiency of mixing by convective overshoot. The convective transport efficiency was expressed in terms of an equivalent mixing-length parameter alpha in the formulation of mixing-length theory (MLT) given by Mihalas (1978). Alpha amounts to values around 2 for matching the entropy of the deep, adiabatically stratified regions of the convective envelope, and lies between 2.5 and 3.0 for matching the thermal structure of the deep photosphere. For current spectral analysis of PMS objects this implies that MLT models based on alpha=2.0 overestimate the effective temperature by 100 K and surface gravities by 0.25 dex. The average thermal structure of the formally convectively stable layers is little affected by convective overshoot and wave heating, i.e., stays close to radiative equilibrium conditions. Our models suggest that the rate of mixing by convective overshoot declines exponentially with geometrical distance to the Schwarzschild stability boundary. It increases at given effective temperature with decreasing gravitational acceleration.

We describe a method for accurately determining M dwarf metallicities with spectral synthesis based on abundance analyses of visual binary stars. We obtained high resolution, high signal-to-noise spectra of each component of five visual binary pairs at McDonald Observatory. The spectral types of the components range from F7 to K3 V for the primaries and M0.5 to M3.5 V for the secondaries. We have determined the metallicities of the primaries differentially with respect to the Sun by fitting synthetic spectra to Fe \textsc{i} line profiles in the observed spectra. In the course of our analysis of the M dwarf secondaries we have made significant improvements to the PHOENIX cool-star model atmospheres and the spectrum analysis code MOOG. Our analysis yields a RMS deviation of 0.11 dex in metallicity values between the binary pairs. We estimate the uncertainties in the derived stellar parameters for the M dwarfs to be 48 K, 0.10 dex, 0.12 dex, 0.15 km s$^{-1}$, and 0.20 km s$^{-1}$ for {T}$_{eff}$, log \textsl{g}, [M/H], $\xi$, and $\eta$ respectively. Accurate stellar evolutionary models are needed to progress further in the analysis of cool-star spectra; the new model atmospheres warrant recalculation of the evolutionary models.

We present work on H alpha spectral line characteristics in PHOENIX stellar model atmospheres and their comparison to microlensing observations. We examine in detail the H alpha equivalent width (EW) and the line shape characteristics for effective temperatures of 4500K< Teff < 5600K where H alpha is a strong spectral feature. We find that H alpha EW in models calculated under the assumption of local thermodynamic equilibrium (LTE) is up to 15% smaller than in models without this assumption, non-LTE models (NLTE) and that line shapes vary significantly for the two model types. A comparison with available high quality microlensing data, capable of tracing H alpha absorption across the face of one G5III giant, shows that the LTE model that fits the EW best is about 100K hotter than and the best-fitting NLTE model has a similar Teff as predicted by the spectral type analysis of the observed star but agree within the uncertainties of the observationally derived temperature. Neither LTE nor NLTE models fit the line shape well. We suspect unmodelled chromospheric emission. Line shape diagnostics suggest lower gravities than derived for the star and are unacceptable low in the case of the LTE models. We show that EW alone is insufficient for comparison to stellar model atmospheres, but combined with a new shape parameter we define is promising. In stellar parameter ranges where the H alpha line is strong, a NLTE approach of modeling stellar atmospheres is not only beneficial but mandatory.

(Abridged) We investigate the effects of metallicity on the broad-band photometric colors of late-type giants, and make a comparison of synthetic colors with observed photometric properties of late-type giants over a wide range of effective temperatures (T_eff=3500-4800 K) and gravities (log g=0.0-2.5), at [M/H]=-1.0 and -2.0. The influence of metallicity on the synthetic photometric colors is generally small at effective temperatures above \~3800 K, but the effects grow larger at lower T_eff, due to the changing efficiency of molecule formation which reduces molecular opacities at lower [M/H]. To make a detailed comparison of the synthetic and observed photometric colors of late type giants in the T_eff--color and color--color planes, we derive a set of new T_eff--log g--color relations based on synthetic photometric colors, at [M/H]=-0.5, -1.0, -1.5, and -2.0. While differences between the new T_eff--color relations and those available from the literature are typically well within ~100 K, effective temperatures predicted by the scales based on synthetic colors tend to be slightly higher than those resulting from the T_eff--color relations based on observations, with the offsets up to ~100 K. This is clearly seen both at [M/H]=-1.0 and -2.0, especially in the T_eff--(B-V) and T_eff--(V-K) planes. The consistency between T_eff--log g--color scales based on synthetic colors calculated with different stellar atmosphere codes is very good, with typical differences being well within \Delta T_eff~70 K at [M/H]=-1.0 and \Delta T_eff~40 K at [M/H]=-2.0.

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.

We present atmospheric models of red giant stars of various metallicities, including extremely metal poor (XMP, [Fe/H]<-3.5) models, with many chemical species, including, significantly, the first two ionization stages of Strontium (Sr) and Barium (Ba), treated in Non-Local Thermodynamic Equilibrium (NLTE) with various degrees of realism. We conclude that 1) for all lines that are useful Sr and Ba abundance diagnostics the magnitude and sense of the computed NLTE effect on the predicted line strength is metallicity dependent, 2) the indirect NLTE effect of overlap between Ba and Sr transitions and transitions of other species that are also treated in NLTE non-negligibly enhances NLTE abundance corrections for some lines, 3) the indirect NLTE effect of NLTE opacity of other species on the equilibrium structure of the atmospheric model is not significant, 4) the computed NLTE line strengths differ negligibly if collisional b-b and b-f rates are an order of magnitude smaller or larger than those calculated with standard analytic formulae, and 5) the effect of NLTE upon the resonance line of Ba II at 4554.03 AA is independent of whether that line is treated with hyperfine splitting. As a result, the derivation of abundances of Ba and Sr for metal-poor red giant stars with LTE modeling that are in the literature should be treated with caution.

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.

To assess the current status in the theoretical modeling of the spectral properties of late-type giants, we provide a comparison of synthetic photometric colors of late-type giants (calculated with PHOENIX, MARCS and ATLAS model atmospheres) with observations, at [M/H]=0.0 and -2.0. Overall, there is a good agreement between observed and synthetic colors, and synthetic colors and published Teff-color relations, both at [M/H]=0.0 and -2.0. Deviations from the observed trends in Teff-color planes are generally within \pm 150K (or less) in the effective temperature range Teff=3500-4800K. Synthetic colors calculated with different stellar atmosphere models typically agree to ~100K, within a large range of effective temperatures and gravities. Some discrepancies are seen in the Teff-(B-V) plane below Teff~3800K at [M/H]=0.0, due to difficulties in reproducing the 'turn-off' to the bluer colors which is seen in the observed data at Teff~3600K. Note that at [M/H]=-2.0 effective temperatures given by the scale of Alonso et al. (1999) are generally lower than those resulting from other Teff-color relations based both on observed and synthetic colors.

We show that contrary to what is expected from 1D stationary model atmospheres, 3D hydrodynamical modeling predicts a considerable influence of convection on the spectral properties of late-type giants. This is due to the fact that convection overshoots into the formally stable outer atmospheric layers producing a notable granulation pattern in the 3D hydrodynamical models, which has a direct influence on the observable spectra and colors. Within the framework of standard 1D model atmospheres the average thermal stratification of the 3D hydro model can not be reproduced with any reasonable choice of the mixing length parameter and formulation of the turbulent pressure. The differences in individual photometric colors -- in terms of 3D versus 1D -- reach up to ~0.2 mag, or \Delta Teff~70K. We discuss the impact of full 3D hydrodynamical models on the interpretation of observable properties of late-type giants, briefly mentioning problems and challenges which need to be solved for bringing these models to a routine use within the astronomical community in 5-10 years from now.

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.

Using NACO on the VLT in the imaging mode we have detected an object at a distance of only 0.7 arcsec from GQ Lup. The object turns out to be co-moving. We have taken two K-band spectra with a resolution of lambda /Delta lambda=700. In here, we analyze the spectra in detail. We show that the shape of spectrum is not spoiled by differences in the Strehl ratio in the blue and in the red part, as well as differential refraction. We reanalyze the spectra and derive the spectral type of the companion using classical methods. We find that the object has a spectral type between M9V and L4V, which corresponds to a Teff between 1600 and 2500 K. Using GAIA-dusty models, we find that the spectral type derivation is robust against different log(g)-values. The Teff derived from the models is again in the range between 1800 and 2400 K. While the models reproduce nicely the general shape of the spectrum, the 12CO-lines in the spectrum have about half the depth as those in the model. We speculate that this difference might be caused by veiling, like in other objects of similar age, and spectral class. We also find that the absolute brightness of the companion matches that of other low-mass free-floating objects of similar age and spectral type. A comparison with the objects in USco observed by Mohanty et al. (2004) shows that the companion of GQ Lup has a lower mass than any of these, as it is of later spectral type, and younger. The same is as true, for the companion of AB Pic. To have a first estimate of the mass of the object we compare the derived Teff and luminosity with those calculated from evolutionary tracks. We also point out that future instruments, like NAHUAL, will finally allow us to derive the masses of such objects more precisely.

We present new synthetic broad-band photometric colors for late-type giants based on synthetic spectra calculated with the PHOENIX model atmosphere code. The grid covers effective temperatures T_eff=3000...5000K, gravities log g=-0.5...+3.5, and metallicities [M/H]=+0.5...-4.0. We show that individual broad-band photometric colors are strongly affected by model parameters such as molecular opacities, gravity, microturbulent velocity, and stellar mass. Our exploratory 3D modeling of a prototypical late-type giant shows that convection has a noticeable effect on the photometric colors too, as it alters significantly both the vertical and horizontal thermal structures in the outer atmosphere. The differences between colors calculated with full 3D hydrodynamical and 1D model atmospheres are significant (e.g., \Delta(V-K)~0.2 mag), translating into offsets in effective temperature of up to ~70K. For a sample of 74 late-type giants in the Solar neighborhood, with interferometric effective temperatures and broad-band photometry available in the literature, we compare observed colors with a new PHOENIX grid of synthetic photometric colors, as well as with photometric colors calculated with the MARCS and ATLAS model atmosphere codes. (abridged)

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.

In this review, we examine the successes and weaknesses of modern low-mass star and brown dwarf theory.(1) We first focus on the mechanical (equation of state) and thermal (atmosphere) properties and on the evolution. We then examine the current shortcomings of the theory and we discuss recent observational analysis which have suggested discrepancies between models and observations.(2) We then examine the stellar and brown dwarf IMF and suggest that a power-law above the average thermal Jeans mass (about 1 Msol) rolling over a lognormal form below this limit adequately reproduces the observations of field and young cluster stellar and brown dwarf distributions. This yields a reasonably accurate estimate of the stellar and brown dwarf Galactic census. Finally (3) we argue that the combination of turbulence driven fragmentation at large scale and gravity at small scales provides an appealing solution for the general star and brown dwarf formation mechanism. It also provides a physical ground for the aforementioned power-law + lognormal form for the IMF, whereas a series of different power laws lacks such a physical motivation. At last, we argue that the deuterium-burning limit as the distinction between stars and planets has no physical foundation in this modern star formation scheme. Opacity limited fragmentation extending down to a few (< 10) jupiter masses, due to shocks, anisotropy or magnetic fields, provides a much more robust limit, even though difficult to determine accurately. Therefore, the various "direct" detections of exoplanets claimed recently in the literature are most likely regular low-mass brown dwarfs and the direct detection of an extrasolar planet remains for now elusive.

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.

Recently the Spitzer Space Telescope observed the transiting extrasolar planets, TrES-1 and HD209458b. These observations have provided the first estimates of the day side thermal flux from two extrasolar planets orbiting Sun-like stars. In this paper, synthetic spectra from atmospheric models are compared to these observations. The day-night temperature difference is explored and phase-dependent flux densities are predicted for both planets. For HD209458b and TrES-1, models with significant day-to-night energy redistribution are required to reproduce the observations. However, the observational error bars are large and a range of models remains viable.

We present a comparison of our VLT/NACO K-band spectrum of the GQ Lupi companion with the new GAIA-dusty model atmosphere grid for T=2000 and 2900 K and log g from 0 to 4. Then, we discuss the mass estimate for GQ Lup companion.

We report a measurement of the shape of the source star in microlensing event MOA 2002-BLG-33. The lens for this event was a close binary whose centre-of-mass passed almost directly in front of the source star. At this time, the source star was closely bounded on all sides by a caustic of the lens. This allowed the oblateness of the source star to be constrained. We found that a/b = 1.02^{+0.04}_{-0.02} where a and b are its semi-major and semi-minor axes respectively. The angular resolution of this measurement is approximately 0.04 microarcsec. We also report HST images of the event that confirm a previous identification of the source star as an F8-G2 turn-off main-sequence star.

We present semi-empirical model chromospheres computed with the atmosphere code PHOENIX. The models are designed to fit the observed spectra of five mid- to late-type M dwarfs. Next to hydrogen lines from the Balmer series we used various metal lines, e. g. from Fe {\sc i}, for the comparison between data and models. Our computations show that an NLTE treatment of C, N, O impacts on the hydrogen line formation, while NLTE treatment of less abundant metals such as nickel influences the lines of the considered species itself. For our coolest models we investigated also the influence of dust on the chromospheres and found that dust increases the emission line flux. Moreover we present an (electronically published) emission line list for the spectral range of 3100 to 3900 and 4700 to 6800 \AA for a set of 21 M dwarfs and brown dwarfs. The line list includes the detection of the Na {\sc i} D lines in emission for a L3 dwarf.

We compare evolutionary models for close-in exoplanets coupling irradiation and evaporation due respectively to the thermal and high energy flux of the parent star with observations of recently discovered new transiting planets. The models provide an overall good agreement with observations, although at the very limit of the quoted error bars of OGLE-TR-10, depending on its age. Using the same general theory, we show that the three recently detected hot-Neptune planets (GJ436, $\rho$ Cancri, $\mu$ Ara) may originate from more massive gas giants which have undergone significant evaporation. We thus suggest that hot-Neptunes and hot-Jupiters may share the same origin and evolution history. Our scenario provides testable predictions in terms of the mass-radius relationships of these hot-Neptunes.

We present a companion of the \le 2 Myr young classical T Tauri star GQ Lup in the Lupus star forming region at 140 \pm 50 pc from imaging, astrometry, and spectroscopy. With direct K-band imaging using VLT/NACO, we detected an object 6 mag fainter than GQ Lup located 0.7 arc sec west of it. Compared to images obtained 2 to 5 years earlier with Subaru/CIAO and HST/PC, this object shares the proper motion of GQ Lup by 5 and 7 sigma, respectively, hence it is a co-moving companion. Its K-L' color is consistent with a spectral type early to mid L. Our NACO K-band spectrum yields spectral type M9-L4 with H2O and CO absorption, consistent with the new GAIA-Dusty template spectrum for log g = 2 to 3 and Teff = 2000 K with ~2 Rjup radius at ~140 pc, hence few Jupiter masses. Using the theoretical models from Wuchterl & Tscharnuter (2003), Burrows et al. (1997), and Baraffe et al. (2002), the mass lies between 1 and 42 Jupiter masses.

We present the results of a set of model atmospheres and synthetic spectra computed with the PHOENIX code. The models cover a range of effective temperatures 2700 K <= T_eff <= 10 000 K, gravities (-0.5 <= log(g)<= 5.5) and metalicities (-4.0 <= [Z/H] <= +0.5). In addition, variations of alpha elements are considered for each metalicity. The models are computed with a homogeneous set of input data in order to allow for direct relative comparison between the models. For example, all models use a mixing length of l/H_p = 2.0. We provide synthetic spectra with a resolution of 0.2 nm from the UV to the infrared for all models. We give a brief overview of the input physics and show illustrative results. All synthetic spectra are available via ftp: ftp.hs.uni-hamburg.de/pub/outgoing/phoenix/GAIA

Late type giants (RGB/AGB stars) will be important tracers of the Galactic morphology and evolution in the framework of Gaia, as they are intrinsically bright and thus can probe distant stellar populations or those obscured by interstellar extinction. A realistic representation of their atmospheres and spectra with stellar atmosphere models is thus of crucial importance, both for the design and optimization of Gaia instruments, as well as the interpretation of provided astrophysical data. Our analysis of synthetic photometric colors of late-type giants based on PHOENIX, MARCS and ATLAS model atmospheres indicates a general agreement between the current theoretical predictions and observations in the framework of stationary 1-D model atmospheres. Presently available models allow temperature determinations of RGB/AGB stars to an accuracy of ~\pm100 K. In an exploratory study we try to quantify possible residual systematic effects due to the approximations made in 1-D models using full 3-D hydrodynamical models. We find that differences in broad-band photometric colors calculated with 1-D and 3-D models are significant, translating to the offsets in effective temperature of up to ~70 K. Clearly, full 3-D hydrodynamical models will help to alleviate such ambiguities in current theoretical modeling. Additionally, they will allow to study new phenomena, to open qualitatively new windows for stellar astrophysics in the Gaia-era.

Previous computations of low temperature Rosseland and Planck mean opacities from Alexander & Ferguson (1994) are updated and expanded. The new computations include a more complete equation of state with more grain species and updated optical constants. Grains are now explicitly included in thermal equilibrium in the equation of state calculation, which allows for a much wider range of grain compositions to be accurately included than was previously the case. The inclusion of high temperature condensates such as Al$_2$O$_3$ and CaTiO$_3$ significantly affects the total opacity over a narrow range of temperatures before the appearance of the first silicate grains. The new opacity tables are tabulated for temperatures ranging from 30000 K to 500 K with gas densities from 10$^{-4}$ g cm$^{-3}$ to 10$^{-19}$ g cm$^{-3}$. Comparisons with previous Rosseland mean opacity calculations are discussed. At high temperatures, the agreement with OPAL and Opacity Project is quite good. Comparisons at lower temperatures are more divergent as a result of differences in molecular and grain physics included in different calculations. The computation of Planck mean opacities performed with the opacity sampling method are shown to require a very large number of opacity sampling wavelength points; previously published results obtained with fewer wavelength points are shown to be significantly in error. Methods for requesting or obtaining the new tables are provided.

We present the single stellar population (SSP) synthesis results of our new synthetic stellar atmosphere models library with a spectral sampling of 0.3 A, covering the wavelength range from 3000 A to 7000 A for a wide range of metallicities (twice solar, solar, half and 1/10 solar). The stellar library is composed of 1650 spectra computed with the latest improvements in stellar atmospheres. In particular it incorporates non-LTE line-blanketed models for hot (Teff >= 27500 K) and LTE line-blanketed models (Phoenix) for cool (3000 K =< Teff =< 4500 K) stars. Because of the high spectral resolution of this library, evolutionary synthesis models can be used to predict the strength of numerous weak absorption lines, and the evolution of the profiles of the strongest lines over a wide range of ages. The SSP results have been calculated for ages 1 Myr to 17 Gyr using the stellar evolutionary provided by the Geneva and Padova-tracks groups. For young stellar populations, our results have a very detailed coverage of high-temperature stars with similar results for Padova and Geneva isochrones. For intermediate and old stellar populations, our results, once degraded to a lower resolution, are similar to the ones obtained by other groups apart from limitations imposed by the stellar evolutionary physics. The limitations and advantages of our library for the analysis of integrated populations are described. The full set of the stellar library and the evolutionary models are available for retrieval at the websites http://www.iaa.csic.es/~rosa and http://www.iaa.csic.es/~mcs/Sed@, or on request from the first two authors.

We present a library of 1654 high-resolution stellar spectra, with a sampling of 0.3 A and covering the wavelength range from 3000 A ~ to 7000 A The library was computed with the latest improvements in stellar atmospheres, incorporating non-LTE line-blanketed models for hot, massive (Teff >= 27500 K) and line-blanketed models for cool (3000 K =< Teff =< 4500 K) stars. The total coverage of the grid is 3000 K =< Teff =< 55000 K and -0.5 =< log g =< 5.5, for four chemical abundance values: twice solar, solar, half solar and 1/10 solar. Evolutionary synthesis models using this library are presented in a companion paper. We tested the general behavior of the library by calculating and comparing equivalent widths of numerous H and HeI lines, and some of the commonly used metallic indices. We also compared the library with the empirical libraries STELIB and Indo-US. The full set of the synthetic stellar spectra is available from our websites (http://www.iaa.csic.es/~rosa and http://www.astro.iag.usp.br/~lucimara/library.html).

Two peculiar stars showing an apparent extremely broadened and strong NaI D absorption have been discovered in surveys for cool white dwarfs by Oppenheimer et al. (2001) and Harris et al. (SDSS, 2003). We discuss the nature of these objects using PHOENIX atmosphere models for metal-poor brown dwarfs/very low mass stars, and new white dwarf LTE and NLTE models for hydrogen- and helium-dominated atmospheres with metals. These include complete molecular formation in chemical equilibrium and a model for the alkali resonance line broadening based on the damping profiles of Allard et al. (2003), as well as new molecular line opacities for metal hydrides. First results of our calculations indicate good agreement with a hydrogen-dominated WD atmosphere with a Na abundance roughly consistent with a state of high accretion. We analyse deviations of the abundances of Na, K, Mg and Ca from the cosmic pattern and comment on implications of these results for standard accretion scenarios.

We use the PHOENIX code package to model the X-ray spectrum of Nova V4743 Sagittarii observed with the LETGS onboard the Chandra satellite on 19 March 2003. To analyze nova atmospheres and related systems with an underlying nuclear burning envelope at X-ray wavelengths, it was necessary to update the code with new microphysics. We demonstrate that the X-ray emission is dominated by thermal bremsstrahlung and that the hard X-rays are dominated by Fe and N absorption. Preliminary models are calculated assuming solar abundances. It is shown that the models can be used to determine element abundances in the nova ejecta by increasing the absorption in the shell.

We use the PHOENIX code package to model the X-ray spectrum of Nova V4743 Sagittarii observed with the LETGS onboard the Chandra satellite on March 2003. Our atmosphere models are 1D spherical, expanding, line blanketed, and in full NLTE. To analyze nova atmospheres and related systems with an underlying nuclear burning envelope at X-ray wavelengths, it was necessary to update the code with new microphysics, as discussed in this paper. We demonstrate that the X-ray emission is dominated by thermal bremsstrahlung and that the hard X-rays are dominated by Fe and N absorption. The best fit to the observation is provided at a temperature of T_eff = 5.8 x 10^5 K, with L_bol = 50 000 L_sun. The models are calculated for solar abundances. It is shown that the models can be used to determine abundances in the nova ejecta.

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 present LTE and NLTE atmospheric models of a star with solar parameters, and study the effect of treating many thousands of Iron group lines out of LTE on the computed atmospheric structure, overall absolute flux distribution, and the moderately high resolution spectrum in the visible and near UV bands. Our NLTE modeling includes the first two or three ionization stages of 20 chemical elements, up to and including much of the Fe-group, and includes about 20000 Fe I and II lines. We investigate separately the effects of treating the light metals and the Fe-group elements in NLTE. Our main conclusions are that 1) NLTE line blanketed models with direct multi-level NLTE for many actual transitions gives qualitatively similar results as the more approximate treatment of Anderson (1989) for both the Fe statistical equilibrium and the atmospheric temperature structure, 2) models with many Fe lines in NLTE have a temperature structure that agrees more closely with LTE semi-empirical models based on center-to-limb variation and a wide variety of spectra lines, whereas LTE models agree more with semi-empirical models based only on an LTE calculation of the Fe I excitation equilibrium, 3) the NLTE effects of Fe-group elements on the model structure and flux distribution are much more important than the NLTE effects of all the light metals combined, and serve to substantially increases the violet and near UV flux level as a result of NLTE Fe over-ionization. These results suggest that there may still be important UV opacity missing from the models.

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.

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.

Model atmospheres have been computed for M dwarfs that are strongly irradiated by nearby hot companions. A variety of primary and secondary spectral types are explored in addition to models specific to four known systems: GD 245, NN Ser, AA Dor, and UU Sge. This work demonstrates that a dramatic temperature inversion is possible on at least one hemisphere of an irradiated M dwarf and the emergent spectrum will be significantly different from an isolated M dwarf or a black body flux distribution. For the first time, synthetic spectra suitable for direct comparison to high-resolution observations of irradiated M dwarfs in non-mass transferring post-common envelope binaries are presented. The effects of departures from local thermodynamic equilibrium on the Balmer line profiles are also discussed.

We present a complete line list for the X 1Sigma+ electronic ground state of LiCl computed using fully quantum-mechanical techniques. This list includes transition energies and oscillator strengths in the spectral region 0.3-39,640.7 cm-1 for all allowed rovibrational transitions in absorption within the electronic ground state. The calculations were performed using an accurate hybrid potential constructed from a spectral inversion fit of experimental data and from recent multi-reference single- and double-excitation configuration interaction calculations. The line list was incorporated into the stellar atmosphere code PHOENIX to compute spectra for a range of young to old T dwarf models. The possibility of observing a signature of LiCl in absorption near 15.8 microns is addressed and the proposal to use this feature to estimate the total lithium elemental abundance for these cool objects is discussed.

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 include the effect of evaporation in our evolutionary calculations of close-in giant planets, based on a recent model for thermal evaporation taking into account the XUV flux of the parent star (Lammer et al. 2003). Our analysis leads to the existence of a critical mass for a given orbital distance $m_{\rm crit}(a)$ below which the evaporation timescale becomes shorter than the thermal timescale of the planet. For planets with initial masses below $m_{\rm crit}$, evaporation leads to a rapid expansion of the outer layers and of the total planetary radius, speeding up the evaporation process. Consequently, the planet does not survive as long as estimated by a simple application of mass loss rates without following consistently its evolution. We find out that the transit planet HD 209458b might be in such a dramatic phase, although with an extremely small probability. As a consequence, we predict that, after a certain time, only planets above a value $m_{\rm crit}(a)$ should be present at an orbital distance $a$ of a star. For planets with initial masses above $m_{\rm crit}$, evaporation does not affect the evolution of the radius with time.

We discuss high-resolution, time-resolved spectra of the caustic exit of the binary microlensing event OGLE 2002-BUL-69 obtained with UVES on the VLT. The source star is a G5III giant in the Galactic Bulge. During such events, the source star is highly magnified, and a strong differential magnification around the caustic resolves its surface. Using an appropriate model stellar atmosphere generated by the NextGEN code we obtained a model light curve for the caustic exit and compared it with a dense set of photometric observations obtained by the PLANET microlensing follow up network. We further compared predicted variations in the H alpha equivalent width with those measured from our spectra. While the model and observations agree in the gross features, there are discrepancies suggesting shortcomings in the model, particularly for the H alpha line core, where we have detected amplified emission from the stellar chromosphere as the source star's trailing limb exited the caustic. This achievement became possible by the provision of the OGLE-III Early Warning System, a network of small telescopes capable of nearly-continuous round-the-clock photometric monitoring, on-line data reduction, daily near-real-time modelling in order to predict caustic crossing parameters, and a fast and efficient response of a 8m-class telescope to a ``Target-Of-Opportunity'' observation request.

(abridged) We present an analysis of high resolution optical spectra for a sample of very young, mid- to late M, low-mass stellar and substellar objects in Upper Sco and Taurus. Effective temperatures (T_eff) and surface gravities are derived from a multi-feature spectral analysis using TiO, NaI and KI, through comparison with the latest synthetic spectra. In combination, they allow us to determine T_eff to within 50K and gravity to within 0.25 dex. Our high-resolution analysis does not require extinction estimates. Moroever, it yields T_eff and gravities independent of theoretical evolutionary models. We find that our gravities for most of the sample agree remarkably well with the isochrone predictions for the likely cluster ages. However, discrepancies appear in our coolest targets: these appear to have significantly lower gravity (by up to 0.75 dex) than our hotter objects, even though our entire sample covers a relatively narrow range in T_eff (about 300K). This drop in gravity is also implied by inter-comparisons of the data alone, without recourse to synthetic spectra. We consider, and argue against, dust opacity, cool stellar spots or metallicity differences leading to the observed spectral effects; a real decline in gravity is strongly indicated. Such gravity variations are contrary to the predictions of the evolutionary tracks, causing improbably low ages to be inferred from the tracks for our coolest targets. We venture that these results may arise from evolutionary model uncertainties related to accretion, deuterium-burning and/or convection effects. Finally, when combined with photometry and distances, our technique for deriving gravities and temperatures provides a way of obtaining masses and radii for substellar objects independent of evolutionary models, as presented in Paper II.

We review our present understanding of the physical properties of substellar objects, brown dwarfs and irradiated or non-irradiated gaseous exoplanets. This includes a description of their internal properties, mechanical structure and heat content, their atmospheric properties, thermal profile and emergent spectrum, and their evolution, in particular as irradiated companions of a close parent star. The general theory can be used to make predictions in term of detectability for the future observational projects. Special attention is devoted to the evolution of the two presently detected transit planets, HD209458B and OGLE-TR-56B. For this latter, we present a consistent evolution for its recently revised mass and show that we reproduce the observed radius within its error bars. We briefly discuss differences between brown dwarfs and gaseous planets, both in terms of mass function and formation process. We outline several arguments to show that the minimum mass for deuterium burning, recently adopted officially as the limit to distinguish the two types of objects, is unlikely to play any specific role in star formation, so that such a limit is of purely semantic nature and is not supported by a physical justification.

Extending the theory we derived recently for HD209458b to different cases of strongly irradiated gaseous exoplanets, we have calculated the consistent evolution of the new transiting planet, OGLE-TR-56b, for its recently revised mass determination. The theory is shown to successfully reproduce the observed radius, for the proper age of the system. We also examine the dissipation of kinetic energy at the planet's internal adiabat due to atmospheric winds, and place constraints on the efficiency of this process. We show that a fraction $\sim 0.1$ -- 0.5% of the incident flux transformed into thermal energy deposited at the adiabatic level can accommodate the observed radii of both OGLE-TR-56b and HD209458b. The present theory yields quantitative predictions on the evolution of the emergent spectrum and fundamental properties of hot-jupiters. The predictions for radius, luminosity, temperature as a function of the planet's mass and orbital distance can be used as benchmarks for future detections of transit planets.

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.

We report a measurement of limb darkening of a solar-like star in the very high magnification microlensing event MOA 2002-BLG-33. A 15 hour deviation from the light curve profile expected for a single lens was monitored intensively in V and I passbands by five telescopes spanning the globe. Our modelling of the light curve showed the lens to be a close binary system whose centre-of-mass passed almost directly in front of the source star. The source star was identified as an F8-G2 main sequence turn-off star. The measured stellar profiles agree with current stellar atmosphere theory to within ~4% in two passbands. The effective angular resolution of the measurements is <1 micro-arcsec. These are the first limb darkening measurements obtained by microlensing for a Solar-like star.

We report the discovery of a new eclipsing system less than one arcminute south of the pulsating DB white dwarf KUV 05134+2605. The object could be identified with the point source 2MASSJ0516288+260738 published by the Two Micron All Sky Survey. We present and discuss the first light curves as well as some additional colour and spectral information. The eclipse period of the system is 1.29d, and, assuming this to be identical to the orbital period, the best light curve solution yields a mass ratio of m2/m1=0.11, a radius ratio of r2/r1~1 and an inclination of 74 deg. The spectral anaylsis results in a Teff=4200K for the primary. On this basis, we suggest that the new system probably consists of a late K + Brown dwarf (which would imply a system considerably younger than ~0.01 Gyr to have r2/r1~1), and outline possible future observations.

V4743 Sgr (Nova Sgr 2002 No. 3) was discovered on 20 September 2002. We obtained a 5ks ACIS-S spectrum in November 2002 and found that the nova was faint in X-rays. We then obtained a 25ks CHANDRA LETGS observation on 19 March 2003. By this time, it had evolved into the Super Soft X-ray phase exhibiting a continuous spectrum with deep absorption features. The light curve from the observation showed large amplitude oscillations with a period of 1325s (22min) followed by a decline in total count rate after ~13ks of observations. The count rate dropped from ~40cts/s to practically zero within ~6ks and stayed low for the rest of the observation (~6ks). The spectral hardness ratio changed from maxima to minima in correlation with the oscillations, and then became significantly softer during the decay. Strong H-like and He-like lines of oxygen, nitrogen, and carbon were found in absorption during the bright phase, indicating temperatures between 1-2MK, but they were shifted in wavelength corresponding to a Doppler velocity of -2400km/s. The spectrum obtained after the decline in count rate showed emission lines of CVI, NVI, and NVII suggesting that we were seeing expanding gas ejected during the outburst, probably originating from CNO-cycled material. An XMM-Newton ToO observation, obtained on 4 April 2003 and a later LETGS observation from 18 July 2003 also showed oscillations, but with smaller amplitudes.

We present plane-parallel and spherical LTE and NLTE atmospheric models of a variety of stellar parameters of the red giant star Arcturus (alpha Boo, HD124897, HR5340) and study their ability to fit the measured absolute flux distribution. Our NLTE models include tens of thousands of the strongest lines in NLTE, and we investigate separately the effects of treating the light metals and the Fe group elements Fe and Ti in NLTE. We find that the NLTE effects of Fe group elements on the model structure and flux distribution are much more important than the NLTE effects of all the light metals combined, and serve to substantially increases the violet and near UV flux level as a result of NLTE Fe over-ionization. Both the LTE and NLTE models predict significantly more flux in the blue and UV bands than is observed. We find that within the moderately metal-poor metallicity range, the effect of NLTE on the overall UV flux level decreases with decreasing metallicity. These results suggest that there may still be important UV opacity missing from the models. We find that models of solar metallicity giants of similar spectral type to Arcturus fit well the observed flux distributions of those stars from the red to the near UV band. This suggests that the blue and near UV flux discrepancy is metallicity dependent, increasing with decreasing metallicity.

We obtain high-precision limb-darkening measurements in five bands (V, V_E, I_E, I, and H) for the K3 III (Teff=4200 K, [Fe/H]=+0.3, log(g)=2.3) source of the Galactic bulge microlensing event EROS BLG-2000-5. These measurements are inconsistent with the predictions of atmospheric models at >10 sigma. While the disagreement is present in all bands, it is most apparent in I, I_E and V_E, in part because the data are better and in part because the intrinsic disagreement is stronger. We find that when limb-darkening profiles are normalized to have unit total flux, the I-band models for a broad range of temperatures all cross each other at a common point. The solar profile also passes through this point. However, the profile as measured by microlensing does not. We conjecture that the models have incorporated some aspect of solar physics that is not shared by giant atmospheres.

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.

Recently, 59 low-luminosity object transits were reported from the Optical Gravitational Lensing Experiment (OGLE). Our follow-up low-resolution spectroscopy of 16 candidates provided two objects, OGLE-TR-3 and OGLE-TR-10, which have companions with radii compatible with those of gas-giant planets. Further high-resolution spectroscopy revealed a very low velocity variation (<500m/s) of the host star OGLE-TR-3 which may be caused by its unseen companion. An analysis of the radial velocity and light curve results in M<2.5 M_jup, R<1.6 R_jup, and an orbital separation of about 5 R_sol, which makes it the planet with the shortest period known. This allows to identify the low-luminosity companion of OGLE-TR-3 as a possible new gas-giant planet. If confirmed, this makes OGLE-TR-3 together with OGLE-TR-56 the first extrasolar planets detected via their transit light curves.

We present evolutionary models for cool brown dwarfs and extra-solar giant planets. The models reproduce the main trends of observed methane dwarfs in near-IR color-magnitude diagrams. We also present evolutionary models for irradiated planets, coupling for the first time irradiated atmosphere profiles and inner structures. We focus on HD 209458-like systems and show that irradiation effects can substantially affect the radius of sub-jovian mass giant planets. Irradiation effects, however, cannot alone explain the large observed radius of HD 209458b. Adopting assumptions which optimise irradiation effects and taking into account the extension of the outer atmospheric layers, we still find $\sim$ 20% discrepancy between observed and theoretical radii. An extra source of energy seems to be required to explain the observed value of the first transit planet.

We present a multiwavelength study of the ONeMg Galactic nova V382 Velorum 1999 using HST/STIS and FUSE ultraviolet spectra and comparisons with published groundbased optical spectra. We find a close match to the basic phenomenology of another well-studied ONeMg nova, V1974 Cygni (Nova Cygni 1992), in particular to the spectral development through the start of the nebular phase. Following an ``iron curtain'' phase, the nova proceeded through a stage of P Cygni line profiles on all important resonance lines, as in many ONeMg novae and unlike the CO class. Emergent emission lines displayed considerable structure, as seen in V1974 Cyg, indicating fragmentation of the ejecta at the earliest stages of the outburst. Analysis and modeling of our ultraviolet spectra suggest that 4 - 5 $\times$ 10$^{-4}$M$_{\odot}$ of material was ejected and that the distance to the nova is $\simeq 2.5$ kpc. Relative to solar values, we find the following abundances: He = 1.0, C = 0.6$\pm$0.3, N = 17$\pm$4, O = 3.4$\pm$0.3, Ne = 17$\pm$3, Mg = 2.6$\pm$0.1, Al = 21$\pm$2, and Si = 0.5$\pm$0.3. Finally, we briefly draw comparisons with Nova LMC 2000, another ONeMg nova, for which similar data were obtained with HST and FUSE.

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.

V1494 Aql (Nova Aql 1999 No. 2) was discovered on 2 December 1999. We obtained Chandra ACIS-I spectra on 15 April and 7 June 2000 which appear to show only emission lines. Our third observation, on 6 August, showed that its spectrum had evolved to that characteristic of a Super Soft X-ray Source. We then obtained Chandra LETG+HRC-S spectra on 28 September (8 ksec) and 1 October (17 ksec). We analyzed the X-ray light curve of our grating observations and found both a short time scale ``burst'' and oscillations. Neither of these phenomena have previously been seen in the light curve of a nova in outburst. The ``burst'' was a factor of 10 rise in X-ray counts near the middle of the second observation, and which lasted about 1000 sec; it exhibited at least two peaks, in addition to other structure. Our time series analysis of the combined 25 ksec observation shows a peak at 2500 s which is present in independent analyses of both the zeroth order image and the dispersed spectrum and is not present in similar analyses of grating data for HZ 43 and Sirius B. Further analyses of the V1494 Aql data find other periods present which implies that we are observing non-radial g+ modes from the pulsating, rekindled white dwarf.

We have updated our PHOENIX model atmospheres and theoretical spectra for ultracool dwarfs with new opacity data for methane based on line strength predictions with the STDS software. By extending the line list to rotational levels of J=40 we can significantly improve the shape of the near-IR absorption features of CH_4, and in addition find an enhanced blanketing effect, resulting in up to 50% more flux emerging in the J band than seen in previous models, which may thus contribute to the brightening in J and blue IR colors observed in T dwarfs.

We present the current status of PHOENIX model atmospheres for dwarfs of the spectral type T, typical for older field brown dwarfs and low-mass brown dwarfs. The results are based on new predictions of the CH_4 line opacities from theoretical calculations with the STDS software package, extrapolating to transitions from rotational levels up to J=40. While individual line positions and strengths are reproduced with moderate to fair accuracy, the cumulative band strength in the region of the IR methane bands is modelled much better thanks to the inclusion of large numbers of faint lines relevant at high temperatures.

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.

Based on detailed 2D and 3D numerical radiation-hydrodynamics (RHD) simulations of time-dependent compressible convection, we have studied the dynamics and thermal structure of the convective surface layers of a prototypical late-type M-dwarf (Teff~2800K log(g)=5.0, solar chemical composition). The RHD models predict stellar granulation qualitatively similar to the familiar solar pattern. Quantitatively, the granular cells show a convective turn-over time scale of ~100s, and a horizontal scale of 80km; the relative intensity contrast of the granular pattern amounts to 1.1%, and root-mean-square vertical velocities reach 240m/s at maximum. Deviations from radiative equilibrium in the higher, formally convectively stable atmospheric layers are found to be insignificant allowing a reliable modeling of the atmosphere with 1D standard model atmospheres. A mixing-length parameter of alpha=2.1 provides the best representation of the average thermal structure of the RHD model atmosphere while alternative values are found when fitting the asymptotic entropy encountered in deeper layers of the stellar envelope alpha=1.5, or when matching the vertical velocity field alpha=3.5. The close correspondence between RHD and standard model atmospheres implies that presently existing discrepancies between observed and predicted stellar colors in the M-dwarf regime cannot be traced back to an inadequate treatment of convection in the 1D standard models. The RHD models predict a modest extension of the convectively mixed region beyond the formal Schwarzschild stability boundary which provides hints for the distribution of dust grains in cooler (brown dwarf) atmospheres.

Understanding post-common-envelope evolution is important in the studies of close-binary systems. The majority of all interacting binaries with a compact object in their system (e.g. cataclysmic variables, X-ray binaries) are thought to have passed through a common-envelope (CE) phase. Current models of interacting binaries all assume, as a starting point, that there is no significant modification of the secondary star compared with a normal star of similar spectral type. The extent to which the CE significantly alters the composition of the secondary star has yet to be determined. We are studying the M-type secondary in the pre-CV system EG UMa, in order to determine its atmospheric parameters by comparison with synthetic spectra generated using PHOENIX. Absorption lines due to several elements have been used in investigating the stellar parameters of effective temperature, gravity and over elemental abundance. In addition, we are searching for anomalies due to s-process elements (such as Ba, Sr, Rb, Y), which may have been deposited in the atmosphere during the CE stage, and prove that CE evolution has occurred. EG UMa displays strong YO absorption bands, which are normally associated with giants and S stars. These bands may have been formed as a result of an Y-abundance enhancement introduced by the CE, which would be consistent with the observed possible Rb and Sr enrichment.

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.

Recently, 46 low-luminosity object transits were reported from the Optical Gravitational Lensing Experiment. Our follow-up spectroscopy of the 16 most promising candidates provides a spectral classification of the primary. Together with the radius ratio from the transit measurements, we derived the radii of the low-luminosity companions. This allows to examine the possible sub-stellar nature of these objects. Fourteen of them can be clearly identified as low-mass stars. Two objects, OGLE-TR-03 and OGLE-TR-10 have companions with radii of 0.15 R_sun which is very similar to the radius of the transiting planet HD209458B. The planetary nature of these two objects should therefore be confirmed by dynamical mass determinations.

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.

The opacity due to photodissociation of 24MgH is investigated in the atmospheres of cool stars. The lowest two electronic transitions A 2Pi -- X 2Sigma+ and B' 2Sigma+ -- X 2 Sigma+ are considered where the cross sections for the latter were published previously (Weck, Stancil, & Kirby 2002) while the former are presented in this work. Model atmospheres calculated with the PHOENIX code are used to investigate the effect of the photodissociation opacity on spectra of cool stars. The A -- X photodissociation cross sections are obtained using a combination of ab initio and experimentally derived potential curves and dipole transition moments. Partial cross sections have been evaluated over the accessible wavelength range 1770-4560 Angstrom for all rotational transitions from the vibrational levels v''=0-11. Assuming a Boltzmann distribution of the rovibrational levels of the X 2Sigma+ state, LTE photodissociation cross sections are presented for temperatures between 1000 and 5000 K. Shape resonances, arising from rotational predissociation of quasi-bound levels of the A 2Pi state near threshold, characterize the LTE photodissociation cross sections. A sum rule is proposed as a check on the accuracy of the photodissociation calculations.

A new, complete, theoretical rotational and vibrational line list for the A-X electronic transition in MgH is presented. The list includes transition energies and oscillator strengths for all possible allowed transitions and was computed using the best available theoretical potential energies and dipole transition moment function with the former adjusted to account for experimental data. The A-X line list, as well as new line lists for the B'-X and the X-X (pure rovibrational) transitions, were included in comprehensive stellar atmosphere models for M, L, and T dwarfs and solar-type stars. The resulting spectra, when compared to models lacking MgH, show that MgH provides significant opacity in the visible between 4400 and 5600 Angstrom. Further, comparison of the spectra obtained with the current line list to spectra obtained using the line list constructed by Kurucz (1993) show that the Kurucz list significantly overestimates the opacity due to MgH particularly for the bands near 5150 and 4800 Angstrom with the discrepancy increasing with decreasing effective temperature.

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.

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

We present a spectral analysis of a sample of late L dwarfs. We use our latest model atmospheres and synthetic spectra and optical and K band spectra to determine effective temperatures. We derive effective temperatures of 1400 K to 1700 K for L8 to L6 dwarfs. The analysis demonstrates that our recent models that rain out the formed dust completely are applicable to optical spectra of late L dwarfs and that more consistent models are needed for intermediate L dwarfs and for infrared spectra. We compare the results for the effective temperatures with the temperatures of the onset of Methane formation. Our models predict Methane absorption at 3.3micron to occur at about 400 K higher temperatures than Methane absorption at 2.2micron. This is consistent with our data and previous observations which show Methane absorption at 3.3micron but not at 2.2micron in late L dwarfs.

We present opacity sampling model atmospheres, synthetic spectra and colors for brown dwarfs and very low mass stars in two limiting case of dust grain formation: 1) inefficient gravitational settling i.e. the dust is distributed according to the chemical equilibrium predictions, 2) efficient gravitational settling i.e. the dust forms and depletes refractory elements from the gas, but their opacity does not affect the thermal structure. The models include the formation of over 600 gas phase species, and 1000 liquids and crystals, and the opacities of 30 different types of grains including corundum (Al$_2$O$_3$), the magnesium aluminum spinel MgAl$_2$O$_4$, iron, enstatite (MgSiO$_3$), forsterite (Mg$_2$SiO$_4$), amorphous carbon, SiC, and a number of calcium silicates. The models extend from the beginning of the grain formation regime well into the condensation regime of water ice ($\teff= 3000 - 100$ K) and encompasses the range of $\log g= 2.5 - 6.0$ at solar metallicity. We find that silicate dust grains can form abundantly in the outer atmospheric layers of red and brown dwarfs with spectral type later than M8. The greenhouse effects of dust opacities provide a natural explanation for the peculiarly red spectroscopic distribution of the latest M dwarfs and young brown dwarfs. The grainless (Cond) models on the other hand, correspond closely to methane brown dwarfs such as Gliese 229B. We also recover that the $\lambda$5891,5897\AA Na I D and $\lambda$7687,7701\AA K I resonance doublets plays a critical role in T dwarfs where their red wing define the pseudo-continuum from the $I$ to the $Z$ bandpass.

We present broad band spectra of a sample of 21 low luminosity sources in the Trapezium Cluster, with masses in the range 0.008 - 0.10 Msun (assuming an age of 1 Myr). These were selected for low extinction in most cases and are located west of the brighter nebulosity. The spectra are in the H bandpass (1.4-1.95 um) and K bandpass (1.9-2.5 um) also for most of the brighter sources, with a resolution of 50 nm. They were taken with UKIRT using the CGS4 spectrometer. Absorption by water vapour bands is detected in all the substellar candidates except one, which is a highly reddened object with strong H2 emission and an anomalously blue (I-J) colour, implying that it is a very young cluster member with circumstellar matter. The observation of prominent water vapour bands confirms the low Effective Temperatures implied by our (I-J) colour measurements in an earlier paper and would imply late M or L spectral types if these were older field dwarfs. However, the profiles of the H bandpass spectra are very different from those of field dwarfs with similar water absorption strength, demonstrating that they are not foreground or background objects. In addition, the CO absorption bands at 2.3 um and the NaI absorption feature at 2.21 um are very weak for such cool sources. All these features are quite well reproduced by the AMES-Dusty-1999 model atmospheres of Allard et al.(2000,2001), and arise from the low gravities predicted for the Trapezium sources. This represents a new proof of the substellar status of our sources, independent of the statistical arguments for low contamination, which are reexamined here. The very late spectral types of the planetary mass objects and very low mass brown dwarfs demonstrate that they are cluster members, since they are too luminous to be field dwarfs in the background.

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 modeled irradiated planets located near a dM5 and a G2 primary star. The impinging radiation field was explicitly included in the solution of the radiative transfer equation and in the computation of the atmospheric structure. We find that large errors in both the thermal and reflected flux will result from models which do not include the impinging radiation in a self-consistent manner. A cool (teff = 500K) and a hot (teff = 1000K) planet were modeled at various orbital separations from both the dM5 and the G2 primary. In all scenarios, we compared the effects of the irradiation in two extreme cases: one where dust clouds form and remain suspended in the atmosphere, and another where dust clouds form but completely settle out of the atmosphere. The atmospheric structure and emergent spectrum strongly depend on the presence or absence of dust clouds. We find that, in the absence of dust opacity, the impinging radiation significantly alters the innermost layers of an EGP atmosphere and that they are actually brighter in the optical than dusty planets. Our models also indicate that the planet-to-star brightness ratio in the optical will be less that 1 x 10^-5 for objects like tau bootis which is consistent with recently reported upper limit values.

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 present VLT FORS1 spectroscopy that temporally resolves the second caustic crossing of the Bulge K giant source of microlensing event EROS 2000-BLG-5, the first time this has been accomplished for several phases of a caustic transit. The ~1 angstrom H-alpha equivalent width of the source star increases slightly as the center of the star egresses the caustic and then plummets by 30% during the final limb crossing. These changes are not seen in contemporaneous spectra of control stars in the FORS1 slit, but are qualitatively consistent with expectations from stellar atmosphere models as the caustic differentially magnifies different portions of the stellar face of the target. Observations such as these in a variety of stellar lines are equivalent to atmospheric tomography and are expected to provide a direct test of stellar models.

We present analyses of high resolution and medium resolution spectra of early L dwarfs. We used our latest set of model atmospheres to reproduce and analyze the observed features. We can model the optical flux and the atomic line profiles with the best accuracy to date. The models used to reproduce the observations include dust condensation and dust opacities. Compared to previous studies using older models we find that our dust treatment is much improved. The derived parameters for the objects are well in the expected range for old very low mass objects. This is also supported by the absence of Li in most of the objects. For the objects showing Li we can be almost certain that those are brown dwarfs. However, a spectral analysis in general, and this one in particular can only very roughly determine mass and age.

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 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 have obtained 1.0-2.5um spectra at R~600 of 14 disk dwarfs with spectral types M6 to L7. For four of the dwarfs we have also obtained infrared spectra at R~3000 in narrow intervals. In addition, we present new L' photometry for four of the dwarfs in the sample, which allows improved determinations of their bolometric luminosities. We resolve the L-dwarf Denis-P J 0205-1159 into an identical pair of objects separated by 0.35". The spectra, with the published energy distribution for one other dwarf, are compared to synthetic spectra generated by upgraded model atmospheres. Good matches are found for 2200> Teff K>1900 (spectral types around M9 to L3), but discrepancies exist at Teff> 2300 K (M8) and for Teff<1800 K (L4-L7). At the higher temperatures the mismatches are due to incompleteness in the water vapor linelist. At the lower temperatures the disagreement is probably due to our treatment of dust: we assume a photospheric distribution in equilibrium with the gas phase. We derive effective temperatures for the sample from the comparison with synthetic spectra and also by comparing our observed total intrinsic luminosities to structural model calculations (which are mostly independent of the atmosphere but are dependent on the unknown masses and ages of the targets). The two derivations agree to ~200 K except for the faintest object in the sample where the discrepancy is larger. Agreement with other temperature determinations is also ~200 K, except for the L7 dwarf.

We have written a light curve synthesis code that makes direct use of model atmosphere specific intensities, in particular the NEXTGEN model atmosphere grid for cool giants (T_eff <= 6800 K and log(g) <= 3.5, Hauschildt et al. 1999). We point out that these models (computed using spherical geometry) predict a limb darkening behaviour that deviates significantly from a simple linear or two-parameter law (there is less intensity at the limb of the star). The presence of a significantly nonlinear limb darkening law has two main consequences. First, the ellipsoidal light curve computed for a tidally distorted giant using the NEXTGEN intensities is in general different from the light curve computed using the same geometry but with the black body approximation and a one- or two-parameter limb darkening law. In most cases the light curves computed with the NEXTGEN intensities have deeper minima than their black body counterparts. Thus the light curve solutions for binaries with a giant component obtained with models with near linear limb darkening (either black body or plane-parallel model atmosphere intensities) are biased. Observations over a wide wavelength range (i.e. both the optical and infrared) are particularly useful in discriminating between models with nearly linear limb darkening and the NEXTGEN models. Second, we show that rotational broadening kernels for Roche lobe filling (or nearly filling) giants can be significantly different from analytic kernels due to a combination of the nonspherical shape of the star and the radical departure from a simple limb darkening law. As a result, geometrical information inferred from V_rot*sin(i) measurements of cool giants in binary systems are likewise biased.

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.

We consider the depletion of primordial deuterium in the interior of substellar objects as a function of mass, age and absolute magnitude in several photometric passbands. We characterize potential spectroscopic signatures of deuterium in the lines of deuterated water HDO. These results will serve as a useful, independent diagnostic to characterize the mass and/or the age of young substellar objects, and to provide an independent age determination of very young clusters. These results can serve to identify objects at the deuterium-burning limit and to confront the theoretical prediction that D-burning is a necessary condition to form star-like objects.

We present a grid of spherically symmetric model atmospheres for young pre-MS stars. This grid spans the parameter range $2000\K \leq \Teff \leq 6800\K$ and $2.0 \leq \logg \leq 3.5$ for $M=0.1\Msun$, appropriate for low mass stars and brown dwarfs. A major improvement is the replacement of TiO and \water line lists with the newer line list calculated by the NASA-AMES group, for TiO (about 175 million lines of 5 isotopes) and for \water (about 350 million lines in 2 isotopes). We provide the model structures, spectra and broad-band colors in standard filters in electronic form.

We compare the structures of model atmospheres and synthetic spectra calculated using different line lists for TiO and water vapor. We discuss the effects of different line list combinations on the model strutures and spectra for both dwarf and giant stars. It is shown that recent improvements result in significantly improved spectra, in particular in the optical where TiO bands are important. The water vapor dominated near-IR region remains problematic as the current water line lists do not yet completely reproduce the shapes of the observed spectra. We find that the AMES TiO list provides more opacity in most bands and that the new, smaller oscillator strengths lead to systematically cooler temperatures for early type M dwarfs than previous models. These effects combine and will help to siginificantly improve the fits of models to observations in the optical as well as result in improved synthetic photometry of M stars. We show that the Davis 1986 \fel-values for the $\delta$ and $\phi$ bands of TiO best reproduce the observed (V-I) color indices.

We present evolutionary models for low mass stars and brown dwarfs ($m \le 1.2 \msol$) based on recent improvement of the theory: equation of state, atmosphere models, ... We concentrate on early evolutionary phases from the initial deuterium burning phase to the zero-age Main Sequence. Evolutionary models for young brown dwarfs are also presented. We discuss the uncertainties of the present models. We analyse the difficulties arising when comparing models with observations for very young objects, in particular concerning the problem of reddening.

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 present evolutionary calculations for very-low-mass stars and brown dwarfs based on synthetic spectra and non-grey atmosphere models which include dust formation and opacity, i.e. objects with $\te\simle 2800$ K. The interior of the most massive brown dwarfs is shown to develop a conductive core after $\sim 2$ Gyr which slows down their cooling. Comparison is made in optical and infrared color-magnitude diagrams with recent late-M and L-dwarf observations. The saturation in optical colors and the very red near-infrared colors of these objects are well explained by the onset of dust formation in the atmosphere. Comparison of the faintest presently observed L-dwarfs with these dusty evolutionary models suggests that dynamical processes such as turbulent diffusion and gravitational settling are taking place near the photosphere. As the effective temperature decreases below $\te\approx 1300-1400$ K, the colors of these objects move to very blue near-infrared colors, a consequence of the ongoing methane absorption in the infrared. We suggest the possibility ofa brown dwarf dearth in $J,H,K$ color-magnitude diagrams around this temperature.

We present Keck HIRES spectra of 6 late-M dwarfs and 11 L dwarfs. Our goal is to assign effective temperatures to the objects using detailed atmospheric models and fine analysis of the alkali resonance absorption lines of CsI and RbI. These yield mutually consistent results (+-150 K) when we use ``cleared-dust'' models, which account for the removal of refractory species from the molecular states but do not include dust opacities. We find a tendency for the RbI line to imply a slightly higher temperature, which we ascribe to an incomplete treatment of the overlying molecular opacities. The final effective temperatures we adopt are based on the CsI fits alone, though the RbI fits support the CsI temperature sequence. This work, in combination with results from the infrared, hints that dust in these atmospheres has settled out of the high atmosphere but is present in the deep photosphere. We also derive radial and rotational velocities for all the objects, finding that the previously discovered trend of rapid rotation for very low mass objects is quite pervasive. To improve on our analysis, there is a clear need for better molecular line lists and a more detailed understanding of dust formation and dynamics.

We have obtained infrared (1 to 2.5 micron) spectroscopy for 42 halo and disk dwarfs with spectral type M1 to M6.5. These data are compared to synthetic spectra generated by the latest model atmospheres of Allard & Hauschildt. Photospheric parameters metallicity, effective temperature and radius are determined for the sample. We find good agreement between observation and theory except for known problems due to incomplete molecular data for metal hydrides and water. The metal-poor M subdwarfs are well matched by the models as oxide opacity sources are less important in this case. The derived effective temperatures for the sample range from 3600K to 2600K; at these temperatures grain formation and extinction are not significant in the photosphere. The derived metallicities range from solar to one-tenth solar. The radii and effective temperatures derived agree well with recent models of low mass stars.

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.

The visual surface brightness of K/M giants and dwarfs with near-solar metallicity differ slightly in agreement with the gravity effects predicted by recent theoretical models. We show that M-dwarfs display also a metallicity dependence of the surface brightness in the infrared $K$-band in agreement with theory. Based on these results, we present improved Barnes-Evans type relations and estimate the radii of 60 single or presumed M and K-dwarfs

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 briefly summarize our present knowledge of the theory of low-mass stars and substellar objects and their contribution to the Galactic population.

We report photometry and spectroscopy of the evolution of Nova Sagittarii 1994 #1 (V4332 Sagittarii) during outburst. We compare the photometric and spectral evolution of this outburst to known classes of outbursts -- including classical novae and outbursts occurring on symbiotic stars -- and find this object does NOT conform to any known class of outburst. The closest match to the behavior of this unusual object is M31 RV, an extremely luminous and red variable object discovered in the bulge of M31 in 1988. However, the temporal behavior and maximum luminosity of the two events differ by several orders of magnitude, requiring substantial intrinsic variation if these two events are members the same type of outburst. Our model of the spectroscopic evolution of this outburst shows that the effective temperature cooled from 4400 K to 2300 K over the three month span of our observations. In combination with line diagnostics in our later spectra, including [OI] lambda 5577 and the dramatic increase in the Halpha to Hbeta ratio, we infer the existence of a cool, dense (N_e ~ 10^{8-9} cm^{-3}) envelope that is optically thick in the Hydrogen Balmer recombination lines (case C). We suggest that a nuclear event in a single star, in which a slow shock drove the photosphere outwards, can power the observed luminosity evolution and the emission spectrum.

We have obtained high resolution echelle spectroscopy of the recurrent nova T CrB. We find that the surface lithium abundance in T CrB is signifcantly enhanced compared to field M giants, where it is not detectable. We offer possible explanations for this in terms of either a delay in the onset of convection in the giant star, enhanced coronal activity due to star-spots or the enhancement of Li resulting from the nova explosion(s).

Based on consistent evolutionary and pulsation calculations, we analyse the effect of metallicity and of different convection treatments in the stellar models on period - magnitude, - color and - radius relationships. In order to perform an accurate comparison with observations, we have computed grids of atmosphere models and synthetic spectra for different metallicities, covering the range of effective temperatures and gravities relevant for Cepheids. The models are compared to recent observations of galactic and Magellanic Clouds Cepheids. Unprecedented level of agreement is found between models and observations. We show that within the range of metallicity for the Galaxy and the Magellanic Clouds, a change of slope in the period - luminosity (PL) relationship is predicted at low periods, due to the reduction of the blue loop during core He burning. The minimum mass undergoing a blue loop and consequently the critical period at which this change of slope occurs depend on the metallicity Z and on the convection treatment in the stellar models. However, besides this change of slope, we do not find any significant effect of metallicity on period - magnitude relationships from V to K bands, and on period - color relationships in IR colors. We only find a detectable effect of $Z$ on $(B-V)$ colors. These results are not affected by uncertainties inherent to current stellar models, mainly due to convection treatment.

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.

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.

We present new infrared photometry for 61 halo and disk stars around the stellar/substellar boundary. These data are combined with available optical photometry and astrometric data to produce color--color and absolute magnitude--color diagrams. The disk and halo sequences are compared to the predictions of the latest model atmospheres and structural models. We find good agreement between observation and theory except for known problems in the V and H passbands probably due to incomplete molecular data for TiO, metal hydrides and H$_2$O. The metal--poor M subdwarfs are well matched by the models as oxide opacity sources are less important in this case. The known extreme M subdwarfs have metallicities about one--hundredth solar, and the coolest subdwarfs have T$_{eff}\sim 3000$ K with masses $\sim $0.09M/M$_{\odot}$. The grainless models are not able to reproduce the flux distributions of disk objects with T$_{eff} <$ 2500 K, however a preliminary version of the NextGen--Dusty models which includes homogeneous formation and extinction by dust grains {\it is} able to match the colors of these very cool objects. The least luminous objects in this sample are GD165B, three DENIS objects --- DBD0205, DBD1058 and DBD1228 --- and Kelu-1. These have T$_{eff}\sim$ 2000 K and are at or below the stellar limit with masses $\leq$0.075M/M$_{\odot}$. Photometry alone cannot constrain these parameters further as the age is unknown, but published lithium detections for two of these objects (Kelu-1 and DBD1228) imply that they are young (aged about 1 Gyr) and substellar (mass $\leq$0.06M/M$_{\odot}$).

We present evolutionary models for low mass stars from 0.075 to 1 $\msol$ for solar-type metallicities [M/H]= 0 and -0.5. The calculations include the most recent interior physics and the latest generation of non-grey atmosphere models. We provide mass-age-color-magnitude relationships for both metallicities. The mass-M$_V$ and mass-M$_K$ relations are in excellent agreement with the empirical relations derived observationally. The theoretical color-magnitude diagrams are compared with the sequences of globular clusters (47 Tucanae) and open clusters (NGC2420 and NGC2477) observed with the Hubble Space Telescope. Comparison is also made with field star sequences in $M_V$-$(V-I)$, $M_K$-$(I-K)$ and $M_K$-$(J-K)$ diagrams. These comparisons show that the most recent improvements performed in low-mass star atmosphere models yield now reliable stellar models in the near-infrared. These models can be used for metallicity, mass, temperature and luminosity calibrations. Uncertainties still remain, however, in the optical spectral region below $T_{eff} \sim 3700K$, where predicted (V-I) colors are too blue by 0.5 mag for a given magnitude. The possible origins for such a discrepancy, most likely a missing source of opacity in the optical and the onset of grain formation are examined in detail.

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 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 determine the elemental abundances for the ejecta of the slow nova PW Vul 1984. Our technique uses a minimization of the emission line fits of a photoionization model to available ultraviolet, optical and infrared spectra. We find the following abundances (by number) with respect to solar: He/H = 1.0 $\pm$ 0.4, C/H = 7.0 $^{+7}_{-4}$, N/H = 85 $^{+59}_{-41}$ and O/H = 6 $^{+7}_{-2}$. In addition, there is weak evidence for solar Ne and Mg and twice solar Fe. Previous studies (Saizar et al. 1991 and Andre\"{a} et al. 1991, 1994) of PW Vul have yielded considerable differences in their derived elemental abundances for the ejecta. Our abundances fall in between the previous studies. To explain the discrepant abundances, we analyze in detail the data and methods used to obtain the previous results. The abundances of Saizar et al. (1991) are significantly smaller then our values because of the lower electron temperature used by Saizar et al. in deriving elemental abundances from ion abundances. Andre\"{a} et al. (1991) used an ionization correction method to obtain their abundances and verified their results with a photoionization model (Andre\"{a} et al. 1994). Our analysis of their data shows that the absolute fluxes of the optical emission lines used by Andre\"{a} are underestimated by 15% leading to a factor of 2 increase in their derived abundances. We also find the photoionization model used by Andre\"{a} et al. (1994) predicts 2 times more carbon than the photoionization code we used even when fitting the same data with similar model parameters.

In this paper, we review the current theory of very low mass stars model atmospheres including the coolest known M~dwarfs, M~subdwarfs, and brown dwarfs, i.e. T$_{eff} \leq 5,000 $K and $-2.0 \le [M/H] \le +0.0$. We discuss ongoing efforts to incorporate molecular and grain opacities in cool stellar spectra, as well as the latest progress in deriving the effective temperature scale of M~dwarfs. We also present the latest results of the models related to the search for brown dwarfs.

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 have performed evolutionary calculations of very-low-mass stars from 0.08 to 0.8 $\msol$ for different metallicites from [M/H]= -2.0 to -1.0 and we have tabulated the mechanical, thermal and photometric characteristics of these models. The calculations include the most recent interior physics and improved non-grey atmosphere models. The models reproduce the entire main sequences of the globular clusters observed with the Hubble Space Telescope over the afore-mentioned range of metallicity. Comparisons are made in the WFPC2 Flight system including the F555, F606 and F814 filters, and in the standard Johnson-Cousins system. We examine the effects of different physical parameters, mixing-length, $\alpha$-enriched elements, helium fraction, as well as the accuracy of the photometric transformations of the HST data into standard systems. We derive mass-effective temperature and mass-magnitude relationships and we compare the results with the ones obtained with different grey-like approximations. These latter are shown to yield inaccurate relations, in particular near the hydrogen-burning limit. We derive new hydrogen-burning minimum masses, and the corresponding absolute magnitudes, for the different metallicities. We predict color-magnitude diagrams in the infrared NICMOS filters, to be used for the next generation of the HST observations, providing mass-magnitudes relationships in these colors down to the brown-dwarf limit. We show that the expected signature of the stellar to substellar transition in color-magnitude diagrams is a severe blueshift in the infrared colors, due to the increasing collision-induced absorption of molecular hydrogen with increasing density and decreasing temperature.

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.

We use a preliminary version of our ``NextGen'' grid of cool star model atmospheres to compute synthetic line profiles which fit high resolution Keck spectra of the cool M~dwarf VB10 satisfactorily well. We show that the parameters derived from the Keck data are consistent with the parameters derived from lower resolution spectra with larger wavelength coverage. We discuss the treatment of van der Waals broadening in cool, molecular (mostly ${\rm H_2}$) dominated stellar atmospheres. The line profiles are dominated by van der Waals pressure broadening and are a sensitive indicator for the gravity and metallicity. Therefore, the high-resolution Keck spectra are useful for determining the parameters of M dwarfs. There is some ambiguity between the metallicity and gravity. For VB10, we find from the high-resolution spectra that $ 5.0 < \log(g) < 5.5$ and $0 < \left[\case{{\rm M}}{{\rm H}}\right] < +0.5$ for an adopted fixed effective temperature of 2700~K (Schweitzer, 1995), which is consistent with recent interior calculations (Baraffe etal, 1995).

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.

We present preliminary non-grey model atmospheres and interiors for cool brown dwarfs. The resulting synthetic spectra are compared to available spectroscopic and photometric observations of the coolest brown dwarf yet discovered, Gl229B (Nakajima \etal, 1995). Despite recognized shortcomings of current methane opacities, we find that the model spectra provide an excellent fit to most of the spectral features of the brown dwarf. We confirm the presence of methane absorption and the substellar nature of Gl229B. These preliminary models set an upper limit for the effective temperature of 1000~K. We also compute the evolution of brown dwarfs with solar composition and masses from 0.02 to 0.065 $\msol$. While uncertainties in the age of the system yield some undetermination for the mass of Gl229B, the most likely solution is $m\approx 0.04-0.055 \msol$. In any case, we can set an upper limit $m= 0.065\msol$ for a very unlikely age $t=10$ Gyr.

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.

We have calculated a grid of more than 700 model atmospheres valid for a wide range of parameters encompassing the coolest known M~dwarfs, M~subdwarfs and brown dwarf candidates: $1500\le \teff \le 4000\,$K, $3.5\le \log(g)\le 5.5$, and $-4.0\le { [M/H]}\le +0.5$. Our equation of state includes 105 molecules and up to 27 ionization stages of 39 elements. In the calculations of the base grid of model atmospheres presented here, we include over 300 molecular bands of 4 molecules (TiO, VO, CaH, FeH) in the JOLA approximation, the water opacity of Ludwig (1971), collision induced opacities, b-f and f-f atomic processes, as well as about 2 million spectral lines selected from a list with more than 42 million atomic and 24 million molecular (H$_2$, CH, NH, OH, MgH, SiH, C$_2$, CN, CO, SiO) lines. High-resolution synthetic spectra are obtained using an opacity sampling method. The model atmospheres and spectra are calculated with the generalized stellar atmosphere code PHOENIX, assuming LTE, plane-parallel geometry, energy (radiative plus convective) conservation, and hydrostatic equilibrium. The model spectra give close agreement with observations of M~dwarfs across a wide spectral range from the blue to the near-IR, with one notable exception: the fit to the water bands. We discuss several practical applications of our model grid, e.g., broadband colors derived from the synthetic spectra. In light of current efforts to identify genuine brown dwarfs, we also show how low-resolution spectra of cool dwarfs vary with surface gravity, and how the high-resolution line profile of the Li~I resonance doublet depends on the Li abundance.

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.

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.