Atmosphere Code PHOENIX Research Articles

The CARMENES search for exoplanets around M dwarfs

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

Predicting the Extreme Ultraviolet Radiation Environment of Exoplanets around Low-mass Stars: The TRAPPIST-1 System

Abstract 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 Å) spectrum. We construct semiempirical nonlocal thermodynamic equilibrium model spectra of the ultracool M8 star TRAPPIST-1 that span EUV to infrared wavelengths (100 Å−2.5 μ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 Hubble Space Telescope Lyman α observations from TRAPPIST-1 and Galaxy Evolution Explorer UV photometric detections from a set of old M8 stars (>1 Gyr). We find that calibrating the models using both data sets separately yield similar far-ultraviolet and NUV fluxes, and EUV fluxes that range from (1.32–17.4) × 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 semiempirical scaling relationships, while the highest activity model predicts EUV fluxes a factor of 10 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.

Bolometric Light Curves of Peculiar Type II-P Supernovae

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

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.

RADIATION ENERGY BALANCE METHOD FOR CALCULATING THE TIME EVOLUTION OF TYPE Ia SUPERNOVAE DURING THE POST-EXPLOSION PHASE

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

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

Stellar to Substellar Model Atmospheres

AbstractThe spectral transition from Very Low Mass stars (VLMs) to brown dwarfs (BDs) and planetary mass objects (Planemos) requires model atmospheres that can treat line, molecule, and dust-cloud formation with completeness and accuracy. One of the essential problems is the determination of the surface velocity field throughout the main sequence down to the BD and planemo mass regimes. We present local 2D and 3D radiation hydrodynamic simulations using the CO5BOLD code with binned Phoenix gas opacities, forsterite dust formation (and opacities) and rotation. The resulting velocity field vs depth and Teff has been used in the general purpose model atmosphere code Phoenix, adapted in static 1D spherical symmetry for these cool atmospheres. The result is a better understanding of the spectral transition from the stellar to substellar regimes. However, problems remain in reproducing the colors of the dustiest brown dwarfs. The global properties of rotation can change the averaged spectral properties of these objects. Our project for the period 2011-2015 is therefore to develop scaled down global 3D simulations of convection, cloud formation and rotation thanks to funding by the Agence Nationale de la Recherche in France.

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.

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.

Multi-wavelength observations of a giant flare on CN Leonis

Aims. In M dwarfs, optical emission lines and continua are sensitive to changing chromospheric conditions, e.g., during flares. To study flare conditions for an observed spectrum, a comparison to synthesised spectra from model atmospheres is needed. Methods. Using the stellar atmosphere code PHOENIX, we computed a set of 41 1D NLTE parameterised chromospheric models including the photosphere and parts of the transition region. By comparison of a linear combination of the synthesised spectra and a quiescent (observed) chromosphere to observed UVES/VLT spectra of a giant flare of the M 5.5 dwarf CN Leo (Gl406), we find the best-fitting flare model chromosphere. Results. Our model spectra give a fairly good overall description of the observed continua and emission lines. In the best-fitting model, the temperature minimum is deep in the atmosphere resulting in high electron pressure for the chromospheric flaring area. The inferred chromospheric filling factor of the flare is about 3 percent, which declines during the flare. The photospheric flare filling factor is about 0.3 percent.

On the hydrogen recombination time in Type II supernova atmospheres

Non-local thermodynamic equilibrium radiative transfer calculations of differentially expanding supernovae atmospheres are computationally intensive and are almost universally performed in a time-independent snapshot mode, where both the radiative transfer problem and the rate equations are solved assuming the steady-state approximation. The validity of the steady-state approximation in the rate equations has recently been questioned for Type II supernova (SN II) atmospheres after maximum light on to the plateau. 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 two principal quantum numbers, the classical value of the recombination time is recovered in the case of a multilevel hydrogen atom. We also find that the recombination time at most optical depths is smaller in the case of a multilevel 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 affect the spectra. The Hα 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α profile is roughly independent of whether the steady-state or time-dependent approach is used.

Structure and spectra of irradiated secondaries in close binaries

Context. The standard stellar model atmosphere ignores the influence of external radiation. This assumption, while sufficient for most stars, fails for many short-period binaries. Aims. In setting up combined model atmospheres for close binaries, we want to constrain the parameters of both components, especially in the case of a hot primary component strongly influencing its cool secondary companion. This situation can be found after common envelope evolution (CEE). The status of both components today allows one to retrace the CEE itself. Methods. We used our stellar atmosphere code PHOENIX, which includes the effect of irradiation in its radiation transport equation, to investigate the close binary star UU Sge. We combined our calculated spectra of both components, weighted by their visible size, and adjusted the input parameters until reasonable agreement with observations is reached. Results. We derive a range of 80 000−85 000 K for the effective temperature of the primary (Teff, p) and give a rough estimate for the primary’s abundances, particularly the nitrogen enrichment. The heated day-side of the secondary has an apparent “effective” or equilibrium temperature of 24 000−26 000 K, nearly independent of its intrinsic luminosity. It shows an enhancement in nitrogen and carbon. Conclusions. The evolution of the primary and secondary stars were strongly influenced by the other’s presence. Radiation from the primary on the secondary’s day-side is still an important factor in understanding the secondary’s atmospheric structure.

Methane line opacities in very cool stellar objects

<i>Aims. <i/> We investigate the effects of different line data for methane <sup>12<sup/>CH<sub>4<sub/> on the structures of model atmospheres and low resolution synthetic spectra for ultra-cool substellar objects. <i>Methods. <i/> For each set of methane line data we compare the resulting model atmospheres and spectra computed with the general purpose model atmosphere code PHOENIX. <i>Results. <i/> The new HGW methane line data compares well to the HITRAN2004 data. We find the the HITRAN2004 methane lines are in some bands stronger than the HGW lines, resulting in deeper absorption bands in the synthetic spectra. <i>Conclusions. <i/>

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.

Spectral synthesis of circumstellar disks – application to white dwarf debris disks

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 PHOENIX. We will show the potential of the application of our code to model the spectra of white dwarf debris disks.

Radiative transfer in circumstellar disks

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.

Consistent Simulations of Substellar Atmospheres and Nonequilibrium Dust Cloud Formation

We aim to understand cloud formation in substellar objects. We combined our nonequilibrium, stationary cloud model DRIFT (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 = 1800 K for the L dwarf twin-binary system DENIS J0205–1159, which is up to 500 K 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.

Influence of chemi-ionization and chemi-recombination processes on hydrogen line shapes in M dwarfs

Aims. We study the influence of chemi-ionization in H ∗ (n) + H(1s) collisions and inverse chemi-recombination processes on the population of higher levels and consequently on profiles of hydrogen lines in the atmospheres of late type (M) stars. Methods. Modeling, using general stellar atmosphere code PHOENIX to reveal the importance of the inclusion of such processes. Results. We demonstrate for the first time observationally detectable effects of these chemi-processes on stellar spectra. Conclusions. It is very important to include chemi-ionization and chemi-recombination processes in modeling of atmospheres of late type stars, especially if one wants to use line profiles for diagnostics of stellar photospheres and lower chromospheres.

Analyzing SN 2003Z with PHOENIX

Aims. 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. Methods. Synthetic spectra are calculated with our stellar atmosphere codePHOENIX. It solves the special relativistic equation of radiative transfer, including large NLTE-calculations and l ine 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 A was covered. The spectra were taken on Feb. 4, 5, 9, and 11, 2003. Results. 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.

PHOENIX model chromospheres of mid- to late-type M dwarfs

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 , 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 A for a set of 21 M dwarfs and brown dwarfs. The line list includes the detection of the Na  D lines in emission for a L3 dwarf.