Source author record

P. H. Hauschildt

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

26works
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

26 published item(s)

preprint2020arXiv

The CARMENES search for exoplanets around M dwarfs. The He I infrared triplet lines in PHOENIX models of M2-3 V stars

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$α$, 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.

preprint2016arXiv

Direct Imaging discovery of a second planet candidate around the possibly transiting planet host CVSO 30

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.

preprint2015arXiv

Center-to-limb variation of intensity and polarization in continuum spectra of FGK stars for spherical atmospheres

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 $λ= 4000 - 7000 \rm Å$, 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.

preprint2015arXiv

Identification of the feature that causes the I-band secondary maximum of a type Ia supernova

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

preprint2015arXiv

Search with UVES and XSHOOTER for signatures of the low-mass secondary in the post common-envelope binary AA Dor

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.

preprint2015arXiv

What causes the large extensions of red-supergiant atmospheres? Comparisons of interferometric observations with 1-D hydrostatic, 3-D convection, and 1-D pulsating model atmospheres

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.

preprint2014arXiv

First spectroscopic observations of the sub-stellar companion of the young debris disk star PZ Telescopii

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.

preprint2014arXiv

VLTI/AMBER observations of cold giant stars: atmospheric structures and fundamental parameters

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.

preprint2012arXiv

A 3D radiative transfer framework IX. Time dependence

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.

preprint2012arXiv

Dark matter powered stars: Constraints from the extragalactic background light

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.

preprint2012arXiv

Fundamental properties and atmospheric structure of the red supergiant VY CMa based on VLTI/AMBER spectro-interferometry

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.

preprint2011arXiv

EChO - Exoplanet Characterisation Observatory

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.

preprint2011arXiv

Theoretical light curves of type Ia supernovae

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.

preprint2009arXiv

3D Radiative Transfer with PHOENIX

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.

preprint2009arXiv

Dust in brown dwarfs and extra-solar planets II. Cloud formation for cosmologically evolving abundances

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.

preprint2009arXiv

On the Hydrogen Recombination Time in Type II Supernova Atmospheres

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.

preprint2009arXiv

Physical parameters of T dwarfs derived from high-resolution near-infrared spectra

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.

preprint2009arXiv

Radiative transfer in circumstellar disks - I. 1D models for GQ Lupi

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.

preprint2009arXiv

Time-dependent radiative transfer with PHOENIX

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.

preprint2008arXiv

Non-LTE modeling of the near UV band of late-type stars

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.

preprint2005arXiv

Evidence for a co-moving sub-stellar companion of GQ Lup

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.

preprint2005arXiv

Review on low-mass stars and brown dwarfs

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.

preprint2002arXiv

The Molecular Line Opacity of MgH in Cool Stellar Atmospheres

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.

preprint1998arXiv

Infrared Colors at the Stellar/Substellar Boundary

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

preprint1998arXiv

Numerical Solution of the Expanding Stellar Atmosphere Problem

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.