Source author record

Peter H. Hauschildt

Peter H. Hauschildt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2022arXiv

EMISSA -- Exploring Millimetre Indicators of Solar-Stellar Activity II. Towards a robust indicator of stellar activity

An activity indicator, which can provide a robust quantitative mapping between the stellar activity and the physical properties of its atmosphere, is important in exploring the physics of activity across spectral types. But the common activity indicators show large variability in their values which makes defining a robust quantitative scale difficult. Millimetre (mm) wavelengths probe the different atmospheric layers within the stellar chromosphere providing a tomographic view of the atmospheric dynamics. The project aims to define a robust mm-based activity indicator for the cool main-sequence stars ($\mathrm{T_{eff}} \sim$ 5000 - 7000 K). We derive the mm-brightness temperature ($\mathrm{T_B(ν)}$) spectral indices ($\mathrm{α_{mm}}$) for cool stars including the Sun using archival data in the 30 - 1000 GHz range. The derived values for $\mathrm{α_{mm}}$ are explored as a function of various physical parameters and empirical power-law functions were derived. $\mathrm{α_{mm}}$ estimates were also compared with other activity indicators. Despite the estimation errors, $\mathrm{α_{mm}}$ values could well distinguish the cool stars, unlike common activity indicators. The low estimation errors on the derived trends of $\mathrm{α_{mm}}$ versus physical parameters suggest that $\mathrm{α_{mm}}$ could be a robust activity indicator. $\mathrm{α_{mm}}$, which is linked to chromospheric thermal stratification and activity in cool stars can well distinguish and physically characterise the stars more robustly than common activity indicators. We emphasise the need for multi-frequency data across the mm-band for stars, with a range of physical parameters and gathered at multiple epochs during activity cycles. This will help explore $\mathrm{α_{mm}}$ in a statistically robust manner and study the emergence of chromospheric heating on the main-sequence.

preprint2021arXiv

A 3D radiative transfer framework: XII. Many-core, vector and GPU methods

3D detailed radiative transfer is computationally taxing, since the solution of the radiative transfer equation involves traversing the six dimensional phase space of the 3D domain. With modern supercomputers the hardware available for wallclock speedup is rapidly changing, mostly in response to requirements to minimize the cost of electrical power. Given the variety of modern computing architectures, we aim to develop and adapt algorithms for different computing architectures to improve performance on a wide variety of platforms. We implemented the main time consuming kernels for solving 3D radiative transfer problems for vastly different computing architectures using MPI, OpenMP, OpenACC and vector algorithms. Adapted algorithms lead to massively improved speed for all architectures, making extremely large model calculations easily feasible. These calculations would have previously been considered impossible or prohibitively expensive. Efficient use of modern computing devices is entirely feasible, but unfortunately requires the implementation of specialized algorithms for them.

preprint2016arXiv

An irradiated brown-dwarf companion to an accreting white dwarf

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, $Δ{\rm T}\simeq57$~K, and the maximum difference between the hottest and coolest parts of its surface, of $Δ{\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-$σ$ confidence level, consistent with high reprocessing efficiency, but poor lateral heat redistribution in the donor's atmosphere.

preprint2015arXiv

M dwarfs and the fraction of high carbon-to-oxygen stars in the solar neighbourhood

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.

preprint2014arXiv

A 3D radiative transfer framework: XI. multi-level NLTE

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.

preprint2013arXiv

A new extensive library of PHOENIX stellar atmospheres and synthetic spectra

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.

preprint2012arXiv

A 3D radiative transfer framework: X. Arbitrary Velocity Fields in the Co-moving Frame

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.

preprint2012arXiv

Modeling the near-UV band of GK stars, Paper II: NLTE models

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.

preprint2012arXiv

Near-infrared light curves of type Ia supernovae

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

preprint2011arXiv

A 3D radiative transfer framework: XIII. OpenCL implementation

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.

preprint2011arXiv

Detecting Planets around Very Cool Dwarfs at Near Infrared Wavelengths with the Radial Velocity Technique

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.

preprint2010arXiv

Hydrogen Recombination with Multilevel atoms

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.

preprint2009arXiv

A 3D radiative transfer framework: IV. spherical & cylindrical coordinate systems

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

preprint2009arXiv

A 3D radiative transfer framework: V. Homologous Flows

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.

preprint2009arXiv

A 3D radiative transfer framework: VI. PHOENIX/3D example applications

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.

preprint2002arXiv

Detailed Spectroscopic Analysis of SN 1987A: The Distance to the LMC using the SEAM method

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.

preprint2002arXiv

Detectibility of Hydrogen Mixing in Type Ia Supernova pre-Maximum Spectra

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.

preprint2000arXiv

An Effective Temperature Scale for Late M and L Dwarfs, from Resonance Absorption Lines of CsI and RbI

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.

preprint2000arXiv

Preliminary Spectral Analysis of the Type II Supernova 1999em

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.

preprint2000arXiv

Spherically symmetric model atmospheres for low mass pre-Main Sequence stars with effective temperatures between 2000 and 6800 K

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.

preprint2000arXiv

TiO and H$_2$O absorption lines in cool stellar atmospheres

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 $δ$ and $ϕ$ bands of TiO best reproduce the observed (V-I) color indices.

preprint1997arXiv

Abundance analysis of the slow nova PW Vulpeculae 1984

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ä 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ä et al. (1991) used an ionization correction method to obtain their abundances and verified their results with a photoionization model (Andreä et al. 1994). Our analysis of their data shows that the absolute fluxes of the optical emission lines used by Andreä are underestimated by 15% leading to a factor of 2 increase in their derived abundances. We also find the photoionization model used by Andreä et al. (1994) predicts 2 times more carbon than the photoionization code we used even when fitting the same data with similar model parameters.

preprint1997arXiv

Brown Dwarfs from the Stellar Perspective

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.

preprint1997arXiv

NLTE effects of Ti~I in M dwarfs and giants

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.

preprint1996arXiv

Analysis of Keck High Resolution Spectra of VB10

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

preprint1996arXiv

NLTE Model Atmospheres for M Dwarfs and Giants

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.

preprint1996arXiv

Non-LTE model atmosphere analysis of the early ultraviolet spectra of nova OS Andromedae 1986

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, $ρ\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.

preprint1996arXiv

Non-LTE Treatment of Fe II in Astrophysical Plasmas

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.

preprint1996arXiv

The Effects of Fe II NLTE on Nova Atmospheres and Spectra

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.