Showing 1–50 of 180 results for author: Hauschildt, P H

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

Authors: Atul Mohan, Sven Wedemeyer, Peter H. Hauschildt, Sneha Pandit, Maryam Saberi

Abstract: 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 differen… ▽ More 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. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: Accepted in A&A Letters

Journal ref: A&A 664, L9 (2022)

EMISSA -- Exploring Millimeter Indicators of Solar-Stellar Activity I. The Initial mm-cm Main Sequence Star Sample

Authors: Atul Mohan, Sven Wedemeyer, Sneha Pandit, Maryam Saberi, Peter H. Hauschildt

Abstract: Due to their wide wavelength coverage across the millimetre to centimetre (mm - cm) range and their increased sensitivity, modern interferometric arrays facilitate observations of the thermal and non-thermal emission from different stellar atmospheric layers. We study the spectral energy distribution ($S_{obs}(ν)$) of main sequence stars using archival mm - cm data with the aim to study their atmo… ▽ More Due to their wide wavelength coverage across the millimetre to centimetre (mm - cm) range and their increased sensitivity, modern interferometric arrays facilitate observations of the thermal and non-thermal emission from different stellar atmospheric layers. We study the spectral energy distribution ($S_{obs}(ν)$) of main sequence stars using archival mm - cm data with the aim to study their atmospheric stratification as a function of stellar type. The main-sequence stars with significant detection in mm bands were identified in the ALMA Science Archive. These data were complemented with spectral flux data in the Ultra violet to centimetre range as compiled from various catalogues and observatory archives. We compare the resultant $S_{obs}(ν)$ of each star with a photospheric emission model ($S_{mod}(ν)$) calculated with the PHOENIX code. The departures of $S_{obs}(ν)$ from the model are quantified in terms of a "spectral flux excess" parameter ($ΔS/S_{mod} (ν)$) for every star in the sample. The initial sample consists of 12 stars across a range of spectral type from A1 to M3.5 and the Sun-as-a-star as reference. The stars with $T_{eff} = 4000 - 7000\,K$ (F - M type) showed a systematically higher $ΔS/S_{mod}$ in the mm - cm range, with the values rising with decreasing $ν$. The steepness of this rise is higher for cooler stars, though the fully convective 3000 K star in the sample deviated from this trend. For the A-type stars, $ΔS/S_{mod} \sim 0$ within errors. The high $ΔS/S_{mod}$ in cool stars points to the presence of hotter upper atmospheric layers, i.e. a chromosphere and corona, like for the Sun. The mm - cm $ΔS/S_{mod}$ spectrum offers a way to estimate the efficiency of the heating mechanisms across various atmospheric layers and thereby to understand their structure and activity. We emphasise the need for more mm - cm data. △ Less

Submitted 9 November, 2021; v1 submitted 25 October, 2021; originally announced October 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 655, A113 (2021)

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

Authors: Peter H. Hauschildt, E. Baron

Abstract: 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 architectu… ▽ More 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. △ Less

Submitted 4 February, 2021; v1 submitted 2 February, 2021; originally announced February 2021.

Comments: 31 pages, 8 figures, to appear in Astronomy and Computing

Time series of optical spectra of Nova V659 Sct

Authors: Dennis Jack, Klaus-Peter Schröder, Philippe Eenens, Uwe Wolter, José Nicolás González-Pérez, Jürgen H. M. M. Schmitt, Peter H. Hauschildt

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

Submitted 24 June, 2020; originally announced June 2020.

Comments: Accepted for publication in Astronomical Notes, 11 pages, 8 figures

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

Authors: D. Hintz, B. Fuhrmeister, S. Czesla, J. H. M. M. Schmitt, A. Schweitzer, E. Nagel, E. N. Johnson, J. A. Caballero, M. Zechmeister, S. V. Jeffers, A. Reiners, I. Ribas, P. J. Amado, A. Quirrenbach, G. Anglada-Escudé, F. F. Bauer, V. J. S. Béjar, M. Cortés-Contreras, S. Dreizler, D. Galadí-Enríquez, E. W. Guenther, P. H. Hauschildt, A. Kaminski, M. Kürster, M. Lafarga , et al. (3 additional authors not shown)

Abstract: 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 believe… ▽ More 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. △ Less

Submitted 13 May, 2020; originally announced May 2020.

Comments: Accepted for publication in A&A, 13 pages, 9 figures, 4 tables

Journal ref: A&A 638, A115 (2020)

A giant exoplanet orbiting a very low-mass star challenges planet formation models

Authors: J. C. Morales, A. J. Mustill, I. Ribas, M. B. Davies, A. Reiners, F. F. Bauer, D. Kossakowski, E. Herrero, E. Rodríguez, M. J. López-González, C. Rodríguez-López, V. J. S. Béjar, L. González-Cuesta, R. Luque, E. Pallé, M. Perger, D. Baroch, A. Johansen, H. Klahr, C. Mordasini, G. Anglada-Escudé, J. A. Caballero, M. Cortés-Contreras, S. Dreizler, M. Lafarga , et al. (157 additional authors not shown)

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

Submitted 26 September, 2019; originally announced September 2019.

Comments: Manuscript author version. 41 pages, 11 figures

The CARMENES search for exoplanets around M dwarfs -- Photospheric parameters of target stars from high-resolution spectroscopy. II. Simultaneous multiwavelength range modeling of activity insensitive lines

Authors: V. M. Passegger, A. Schweitzer, D. Shulyak, E. Nagel, P. H. Hauschildt, A. Reiners, P. J. Amado, J. A. Caballero, M. Cortés-Contreras, A. J. Domínguez-Fernández, A. Quirrenbach, I. Ribas, M. Azzaro, G. Anglada-Escudé, F. F. Bauer, V. J. S. Béjar, S. Dreizler, E. W. Guenther, T. Henning, S. V. Jeffers, A. Kaminski, M. Kürster, M. Lafarga, E. L. Martín, D. Montes , et al. (3 additional authors not shown)

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

Submitted 1 July, 2019; originally announced July 2019.

Journal ref: A&A 627, A161 (2019)

The CARMENES search for exoplanets around M dwarfs. Two temperate Earth-mass planet candidates around Teegarden's Star

Authors: M. Zechmeister, S. Dreizler, I. Ribas, A. Reiners, J. A. Caballero, F. F. Bauer, V. J. S. Béjar, L. González-Cuesta, E. Herrero, S. Lalitha, M. J. López-González, R. Luque, J. C. Morales, E. Pallé, E. Rodríguez, C. Rodríguez López, L. Tal-Or, G. Anglada-Escudé, A. Quirrenbach, P. J. Amado, M. Abril, F. J. Aceituno, J. Aceituno, F. J. Alonso-Floriano, M. Ammler-von Eiff , et al. (160 additional authors not shown)

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

Submitted 13 September, 2019; v1 submitted 17 June, 2019; originally announced June 2019.

Comments: A&A 627, A49. 26 pages, 17 figures, 6 tables. Press release available at http://www.astro.physik.uni-goettingen.de/~zechmeister/teegarden/teegarden.html. v2: two authors and one reference added

Journal ref: A&A 627, A49 (2019)

The CARMENES search for exoplanets around M dwarfs. Chromospheric modeling of M2-3 V stars with PHOENIX

Authors: D. Hintz, B. Fuhrmeister, S. Czesla, J. H. M. M. Schmitt, E. N. Johnson, A. Schweitzer, J. A. Caballero, M. Zechmeister, S. V. Jeffers, A. Reiners, I. Ribas, P. J. Amado, A. Quirrenbach, G. Anglada-Escudé, F. F. Bauer, V. J. S. Béjar, M. Cortés-Contreras, S. Dreizler, D. Galadí-Enríquez, E. W. Guenther, P. H. Hauschildt, A. Kaminski, M. Kürster, M. Lafarga, M. López del Fresno , et al. (4 additional authors not shown)

Abstract: 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 framewo… ▽ More 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. △ Less

Submitted 11 February, 2019; originally announced February 2019.

Comments: 16 pages + 6 pages appendix, 13 figures, 2 + 4 tables

Journal ref: A&A 623, A136 (2019)

Predicting the Extreme Ultraviolet Radiation Environment of Exoplanets Around Low-Mass Stars: the TRAPPIST-1 System

Authors: Sarah Peacock, Travis Barman, Evgenya L. Shkolnik, Peter H. Hauschildt, E. Baron

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 present a… ▽ More 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 present a case study of the ultra-cool M8 dwarf star, TRAPPIST-1, which hosts seven transiting short-period terrestrial sized planets whose atmospheres will be probed by the James Webb Space Telescope. We construct semi-empirical non-LTE model spectra of TRAPPIST-1 that span EUV to infrared wavelengths (100 Å - 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 HST Ly$α$ observations from TRAPPIST-1 and GALEX FUV and NUV photometric detections from a set of old M8 stars ($>$1 Gyr). We find that calibrating the models using both data sets separately yield similar FUV and NUV fluxes, and EUV fluxes that range from (1.32 - 17.4) $\times$ 10$^{-14}$ ergs s$^{-1}$ cm$^{-2}$. The results from these models demonstrate that the EUV emission is very sensitive to the temperature structure in the transition region. Our lower activity models predict EUV fluxes similar to previously published estimates derived from semi-empirical scaling relationships, while the highest activity model predicts EUV fluxes a factor of ten higher. Results from this study support the idea that the TRAPPIST-1 habitable zone planets likely do not have much liquid water on their surfaces due to the elevated levels of high energy radiation emitted by the host star. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: 13 pages, 5 figures, accepted to ApJ

Convection Enhances Magnetic Turbulence in AM CVn Accretion Disks

Authors: Matthew S. B. Coleman, Omer Blaes, Shigenobu Hirose, Peter H. Hauschildt

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

Submitted 12 March, 2018; originally announced March 2018.

Comments: 15 pages, 10 figures, accepted for publication in ApJ

The CARMENES search for exoplanets around M dwarfs - HD 147379b: A nearby Neptune in the temperate zone of an early-M dwarf

Authors: A. Reiners, I. Ribas, M. Zechmeister, J. A. Caballero, T. Trifonov, S. Dreizler, J. C. Morales, L. Tal-Or, M. Lafarga, A. Quirrenbach, P. J. Amado, A. Kaminski, S. V. Jeffers, J. Aceituno, V. J. S. Béjar, J. Guàrdia, E. W. Guenther, H. -J. Hagen, D. Montes, V. M. Passegger, W. Seifert, A. Schweitzer, M. Cortés-Contreras, M. Abril, F. J. Alonso-Floriano , et al. (147 additional authors not shown)

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

Submitted 15 December, 2017; originally announced December 2017.

Comments: accepted for publication as A&A Letter

The CARMENES search for exoplanets around M dwarfs: High-resolution optical and near-infrared spectroscopy of 324 survey stars

Authors: A. Reiners, M. Zechmeister, J. A. Caballero, I. Ribas, J. C. Morales, S. V. Jeffers, P. Schöfer, L. Tal-Or, A. Quirrenbach, P. J. Amado, A. Kaminski, W. Seifert, M. Abril, J. Aceituno, F. J. Alonso-Floriano, M. Ammler-von Eiff, R. Antona, G. Anglada-Escudé, H. Anwand-Heerwart, B. Arroyo-Torres, M. Azzaro, D. Baroch, D. Barrado, F. F. Bauer, S. Becerril , et al. (148 additional authors not shown)

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

Submitted 9 February, 2018; v1 submitted 17 November, 2017; originally announced November 2017.

Comments: accepted for publication in A&A, 13 pages plus 40 pages spectral atlas, first 10 atlas pages are reduced in quality to fit arXiv size limit; one CARMENES spectrum for each of the 324 stars is published in electronic format at http://carmenes.cab.inta-csic.es/

Journal ref: A&A 612, A49 (2018)

The CARMENES search for exoplanets around M dwarfs. First visual-channel radial-velocity measurements and orbital parameter updates of seven M-dwarf planetary systems

Authors: T. Trifonov, M. Kürster, M. Zechmeister, L. Tal-Or, J. A. Caballero, A. Quirrenbach, P. J. Amado, I. Ribas, A. Reiners, S. Reffert, S. Dreizler, A. P. Hatzes, A. Kaminski, R. Launhardt, Th. Henning, D. Montes, V. J. S. Béjar, R. Mundt, A. Pavlov, J. H. M. M. Schmitt, W. Seifert, J. C. Morales, G. Nowak, S. V. Jeffers, C. Rodríguez-López , et al. (144 additional authors not shown)

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

Submitted 29 January, 2018; v1 submitted 4 October, 2017; originally announced October 2017.

Comments: Accepted for publication in A&A, 24 pages, 16 figures, 14 tables

Journal ref: A&A 609, A117 (2018)

Study of the variability of Nova V5668 Sgr, based on high resolution spectroscopic monitoring

Authors: D. Jack, J. de J. Robles Pérez, I. De Gennaro Aquino, K. -P. Schröder, U. Wolter, P. Eenens, J. H. M. M. Schmitt, M. Mittag, A. Hempelmann, J. N. González-Pérez, G. Rauw, P. H. Hauschildt

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

Submitted 3 February, 2017; originally announced February 2017.

Comments: 13 pages, 11 figures, accepted for publication in AN

Optical and ultraviolet spectroscopic analysis of SN 2011fe at late times

Authors: Brian Friesen, E. Baron, Jerod T. Parrent, R. C. Thomas, David Branch, Peter Nugent, Peter H. Hauschildt, Ryan J. Foley, Darryl E. Wright, Yen-Chen Pan, Alexei V. Filippenko, Kelsey I. Clubb, Jeffrey M. Silverman, Keiichi Maeda, Isaac Shivvers, Patrick L. Kelly, Daniel P. Cohen, Armin Rest, Daniel Kasen

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

Submitted 16 July, 2016; originally announced July 2016.

Comments: 22 pages, 21 figuress, 1 table, submitted to MNRAS

An irradiated brown-dwarf companion to an accreting white dwarf

Authors: Juan V. Hernández Santisteban, Christian Knigge, Stuart P. Littlefair, Rene P. Breton, Vikram S. Dhillon, Boris T. Gänsicke, Thomas R. Marsh, Magaretha L. Pretorius, John Southworth, Peter H. Hauschildt

Abstract: 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 p… ▽ More 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. △ Less

Submitted 23 May, 2016; originally announced May 2016.

Comments: 19 pages, 4 figures, published in Nature 19 May

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

Authors: T. O. B. Schmidt, R. Neuhäuser, C. Briceño, N. Vogt, St. Raetz, A. Seifahrt, C. Ginski, M. Mugrauer, S. Buder, C. Adam, P. H. Hauschildt, S. Witte, Ch. Helling, J. H. M. M. Schmitt

Abstract: 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 C… ▽ More 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. △ Less

Submitted 17 May, 2016; originally announced May 2016.

Comments: 15 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 593, A75 (2016)

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

Authors: N. M. Kostogryz, I. Milic, S. V. Berdyugina, P. H. Hauschildt

Abstract: 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… ▽ More 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. △ Less

Submitted 23 November, 2015; originally announced November 2015.

Comments: 13 pages, 12 figures, submitted to A&A

Journal ref: A&A 586, A87 (2016)

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

Authors: John E. Gizis, Zachary Marks, Peter H. Hauschildt

Abstract: 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… ▽ More 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. △ Less

Submitted 23 October, 2015; originally announced October 2015.

Comments: MNRAS, accepted

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

Authors: D. Hoyer, T. Rauch, K. Werner, P. H. Hauschildt, J. W. Kruk

Abstract: 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… ▽ More 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. △ Less

Submitted 29 April, 2015; originally announced April 2015.

Comments: 17 pages, 19 figures

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

Authors: D. Jack, E. Baron, P. H. Hauschildt

Abstract: 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 ex… ▽ More 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. △ Less

Submitted 10 March, 2015; originally announced March 2015.

Comments: 7 pages, 8 figures, accepted for publication in MNRAS

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

Authors: B. Arroyo-Torres, M. Wittkowski, A. Chiavassa, M. Scholz, B. Freytag, J. M. Marcaide, P. H. Hauschildt, P. R. Wood, F. J. Abellan

Abstract: 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 inst… ▽ More 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. △ Less

Submitted 7 January, 2015; originally announced January 2015.

Comments: Accepted for publication in A&A

Journal ref: A&A 575, A50 (2015)

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

Authors: B. Arroyo-Torres, I. Martí-Vidal, J. M. Marcaide, M. Wittkowski, J. C. Guirado, P. H. Hauschildt, A. Quirrenbach, J. Fabregat

Abstract: 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-… ▽ More 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. △ Less

Submitted 29 April, 2014; originally announced April 2014.

Comments: Accepted for publication in A&A

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

Authors: Peter H. Hauschildt, E. Baron

Abstract: 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… ▽ More 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. △ Less

Submitted 16 April, 2014; originally announced April 2014.

Comments: 20 pages, 14 figures, A&A, in press

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

Authors: T. O. B. Schmidt, M. Mugrauer, R. Neuhäuser, N. Vogt, S. Witte, P. H. Hauschildt, Ch. Helling, A. Seifahrt

Abstract: 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… ▽ More 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. △ Less

Submitted 10 April, 2014; originally announced April 2014.

Comments: 7 pages, 6 figures, accepted for publication in A&A

The atmospheric structure and fundamental parameters of the red supergiants AH Sco, UY Sct and KW Sgr

Authors: B. Arroyo-Torres, M. Wittkowski, J. M. Marcaide, P. H Hauschildt

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

Submitted 11 June, 2013; v1 submitted 27 May, 2013; originally announced May 2013.

Comments: Accepted for publication in A&A

Journal ref: Arroyo-Torres, B., Wittkowski, M., Marcaide, J. M., & Hauschildt, P. H. 2013, A&A, 554, A76

A new extensive library of PHOENIX stellar atmospheres and synthetic spectra

Authors: Tim-Oliver Husser, Sebastian Wende - von Berg, Stefan Dreizler, Derek Homeier, Ansgar Reiners, Travis Barman, Peter H. Hauschildt

Abstract: 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-turbu… ▽ More 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. △ Less

Submitted 24 April, 2013; v1 submitted 22 March, 2013; originally announced March 2013.

Comments: accepted for publication in A&A, 9 pages, 8 figures

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

Authors: E. Baron, Peter H. Hauschildt, Bin Chen, Sebastian Knop

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 7 November, 2012; v1 submitted 24 October, 2012; originally announced October 2012.

Comments: 9 pages, 9 figures, A&A, in press, new version matches copy edited version, definition restored

A 3D radiative transfer framework IX. Time dependence

Authors: D. Jack, P. H. Hauschildt, E. Baron

Abstract: 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 d… ▽ More 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. △ Less

Submitted 25 September, 2012; originally announced September 2012.

Comments: A&A in press, 7 pages, 14 figures

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

Authors: M. Wittkowski, P. H. Hauschildt, B. Arroyo Torres, J. M. Marcaide

Abstract: 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 un… ▽ More 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. △ Less

Submitted 23 March, 2012; originally announced March 2012.

Comments: 5 pages, 5 figures, accepted for publication in Astronomy and Astrophysics (A&A) as a Letter

A New Extensive Library of Synthetic Stellar Spectra from PHOENIX Atmospheres and its Application to Fitting VLT MUSE Spectra

Authors: T. -O. Husser, S. Kamann, S. Dreizler, Peter. H. Hauschildt

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

Submitted 8 March, 2012; originally announced March 2012.

Comments: 8 pages, 4 figures, International Workshop on Stellar Spectral Libraries 2011

Journal ref: Astronomical Society of India Conference Series (ASICS), P. Prugniel and H. P. Singh (Editors), 2012

Near-infrared light curves of type Ia supernovae

Authors: Dennis Jack, Peter H. Hauschildt, E. Baron

Abstract: 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. Includ… ▽ More 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. △ Less

Submitted 9 January, 2012; originally announced January 2012.

Comments: 5 pages, 12 figures, A&A in press

Dark matter powered stars: Constraints from the extragalactic background light

Authors: A. Maurer, M. Raue, T. Kneiske, D. Elsässer, P. H. Hauschildt, D. Horns

Abstract: 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… ▽ More 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. △ Less

Submitted 5 January, 2012; originally announced January 2012.

Comments: Accepted for publication in ApJ; 7 pages, 5 figures

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

Authors: C. Ian Short, Eamonn A. Campbell, Heather Pickup, Peter H. Hauschildt

Abstract: 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 the… ▽ More 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. △ Less

Submitted 4 January, 2012; originally announced January 2012.

Comments: Submitted to The Astrophysical Journal. Observed spectrophotometric library, and grids of NLTE and LTE) synthetic spectra for GK stars available at http://www.ap.smu.ca/~ishort/PHOENIX

EChO - Exoplanet Characterisation Observatory

Authors: G. Tinetti, J. P. Beaulieu, T. Henning, M. Meyer, G. Micela, I. Ribas, D. Stam, M. Swain, O. Krause, M. Ollivier, E. Pace, B. Swinyard, A. Aylward, R. van Boekel, A. Coradini, T. Encrenaz, I. Snellen, M. R. Zapatero-Osorio, J. Bouwman, J. Y-K. Cho, V. Coudé du Foresto, T. Guillot, M. Lopez-Morales, I. Mueller-Wodarg, E. Palle , et al. (109 additional authors not shown)

Abstract: 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-w… ▽ More 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. △ Less

Submitted 12 December, 2011; originally announced December 2011.

Comments: Accepted for publication in Experimental Astronomy, 23 pages, 15 figures

A 3D radiative transfer framework: XIII. OpenCL implementation

Authors: Peter H. Hauschildt, E. Baron

Abstract: 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 a… ▽ More 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. △ Less

Submitted 29 July, 2011; v1 submitted 28 July, 2011; originally announced July 2011.

Comments: A&A, in press

Theoretical light curves of type Ia supernovae

Authors: D. Jack, P. H. Hauschildt, E. Baron

Abstract: 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… ▽ More 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. △ Less

Submitted 18 May, 2011; v1 submitted 17 May, 2011; originally announced May 2011.

Comments: 12 pages, 28 figures, corrected reference

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

Authors: Florian Rodler, Carlos del Burgo, Soeren Witte, Christiane Helling, Peter H. Hauschildt, Eduardo L. Martin, Carlos Alvarez

Abstract: 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 th… ▽ More 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. △ Less

Submitted 11 May, 2011; originally announced May 2011.

Comments: 8 pages; A&A, 2011

A 3D radiative transfer framework: VII. Arbitrary velocity fields in the Eulerian frame

Authors: A. M. Seelmann, P. H. Hauschildt, E. Baron

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

Submitted 20 July, 2010; originally announced July 2010.

Hydrogen Recombination with Multilevel atoms

Authors: Soma De, E. Baron, Peter H. Hauschildt

Abstract: 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 i… ▽ More 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. △ Less

Submitted 4 May, 2010; originally announced May 2010.

Comments: 29 pages, 12 figures, MNRAS, in press

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

Authors: Peter H. Hauschildt, E. Baron

Abstract: 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 paralle… ▽ More 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. △ Less

Submitted 20 November, 2009; v1 submitted 17 November, 2009; originally announced November 2009.

Comments: fixed PDF generation problem. Also available at tp://ftp.hs.uni-hamburg.de/pub/outgoing/phoenix/preprints/3DRT_paper6.pdf A&A, in press (Ref: AA/2009/13064)

On the Hydrogen Recombination Time in Type II Supernova Atmospheres

Authors: Soma De, E. Baron, P. H. Hauschildt

Abstract: 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… ▽ More 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. △ Less

Submitted 3 October, 2009; originally announced October 2009.

Comments: Abstract shortened for arXiv. 39 pages, 16 figures, MNRAS in press

Comoving-frame radiative transfer in arbitrary velocity fields -- II. Large scale applications

Authors: S. Knop, P. H. Hauschildt, E. Baron

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

Submitted 23 September, 2009; originally announced September 2009.

3D Radiative Transfer with PHOENIX

Authors: E. Baron, Bin Chen, P. H. Hauschildt

Abstract: 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)… ▽ More 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. △ Less

Submitted 28 August, 2009; originally announced August 2009.

Comments: 13 pages, 9 figures, to appear in Recent Directions in Astrophysical Quantitative Spectroscopy and Radiation Hydrodynamics, Ed. I. Hubeny, American Institute of Physics (2009)

Journal ref: AIP Conf.Proc.1171:148-160,2009

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

Authors: S. Witte, Ch. Helling, P. H. Hauschildt

Abstract: 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… ▽ More 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. △ Less

Submitted 25 August, 2009; originally announced August 2009.

Comments: 14 pages, 12 figures

Time-dependent radiative transfer with PHOENIX

Authors: D. Jack, P. H. Hauschildt, E. Baron

Abstract: 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 im… ▽ More 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. △ Less

Submitted 9 July, 2009; originally announced July 2009.

Comments: 7 pages, 12 figures

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

Authors: C. del Burgo, E. L. Martin, M. R. Zapatero Osorio, P. H. Hauschildt

Abstract: 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 wi… ▽ More 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. △ Less

Submitted 25 March, 2009; originally announced March 2009.

Comments: 27 pages, 17 figures, A&A, in press

A 3D radiative transfer framework: V. Homologous Flows

Authors: E. Baron, Peter H. Hauschildt, Bin Chen

Abstract: 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 coordi… ▽ More 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. △ Less

Submitted 13 March, 2009; originally announced March 2009.

Comments: 6 pages, 3 figures, Astronomy & Astrophysics, in press

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

Authors: S. D. Huegelmeyer, S. Dreizler, P. H. Hauschildt, A. Seifahrt, D. Homeier, T. Barman

Abstract: 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 ge… ▽ More 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. △ Less

Submitted 13 March, 2009; originally announced March 2009.

Comments: 9 pages, 9 figures, accepted for publication in Astronomy & Astrophysics