Synthetic Stellar Spectra Research Articles

Atomic data for stellar spectroscopy

AbstractHigh-precision spectroscopy of large stellar samples plays a crucial role for several topical issues in astrophysics, such as studying the chemical evolution of the Milky Way Galaxy. Data are accumulating from instruments that obtain high-quality spectra of stars in the ultraviolet, optical and infrared wavelength regions on a routine basis. The interpretation of these spectra is often based on synthetic stellar spectra, either calculated on the fly or taken from a spectral library. One of the most important ingredients of these spectra is a set of high-quality transition data for numerous species, in particular neutral and singly ionized atoms. We rely heavily on the continuous activities of laboratory astrophysics groups that produce and improve the relevant experimental and theoretical atomic data. As an example, we briefly describe the efforts done in the context of the Gaia-ESO Public Spectroscopic Survey to compile and assess the best available data in a standard way, providing a list of recommended lines for analysis of optical spectra of FGK stars. The line data, together with specialised analysis methods, allow different surveys to obtain abundances with typical precisions of ∼0.1 dex on an industrial scale for ∼10 chemical elements. Several elements with urgent need for better atomic data have been identified.

The Stagger-grid: Synthetic stellar spectra and broad-band photometry

An abstract is not available for this content so a preview has been provided. As you have access to this content, a full PDF is available via the ‘Save PDF’ action button.

Hot star extension to theHubbleSpace Telescope stellar spectral library

Context.Libraries of stellar spectra find many uses in astrophysics, from photometric calibration to stellar population synthesis.Aims.We present low resolution spectra of 40 stars from 0.2μm (ultraviolet) to 1.0μm (near infrared) with excellent fluxing. The stars include normal O-type stars, helium-burning stars, and post-asymptotic giant branch (PAGB) stars.Methods.Spectra were obtained with the Space Telescope Imaging Spectrograph (STIS) installed in theHubbleSpace Telescope (HST) using three low resolution gratings, G230LB, G430L, and G750L. Cosmic ray hits and fringing in the red were corrected. A correction for scattered light was applied, significant only for our coolest stars. Cross-correlation was used to bring the spectra to a common, final, zero velocity wavelength scale. Finally, synthetic stellar spectra were used to estimate line of sight dust extinction to each star, and a five-parameter dust extinction model was fit, or a one-parameter fit in the case of low extinction.Results.These spectra dovetail with the similar Next Generation Stellar Library (NGSL) spectra, extending the NGSL’s coverage of stellar parameters, and extending to helium burning stars and stars that do not fuse.Conclusions.The fitted dust extinction model showed considerable variation from star to star, indicating variations in dust properties for different lines of sight. Interstellar absorption lines are present in most stars, notably MgII.

FITspec: A New Algorithm for the Automated Fit of Synthetic Stellar Spectra for OB Stars

Abstract In this paper we describe the FITspec code, a data mining tool for the automatic fitting of synthetic stellar spectra. The program uses a database of 27,000 cmfgen models of stellar atmospheres arranged in a six-dimensional (6D) space, where each dimension corresponds to one model parameter. From these models a library of 2,835,000 synthetic spectra were generated covering the ultraviolet, optical, and infrared regions of the electromagnetic spectrum. Using FITspec we adjust the effective temperature and the surface gravity. From the 6D array we also get the luminosity, the metallicity, and three parameters for the stellar wind: the terminal velocity ( ), the β exponent of the velocity law, and the clumping filling factor (F cl). Finally, the projected rotational velocity ( ) can be obtained from the library of stellar spectra. Validation of the algorithm was performed by analyzing the spectra of a sample of eight O-type stars taken from the iacob spectroscopic survey of Northern Galactic OB stars. The spectral lines used for the adjustment of the analyzed stars are reproduced with good accuracy. In particular, the effective temperatures calculated with the FITspec are in good agreement with those derived from spectral type and other calibrations for the same stars. The stellar luminosities and projected rotational velocities are also in good agreement with previous quantitative spectroscopic analyses in the literature. An important advantage of FITspec over traditional codes is that the time required for spectral analyses is reduced from months to a few hours.

An application of deep learning in the analysis of stellar spectra

Spectroscopic surveys require fast and efficient analysis methods to maximize their scientific impact. Here we apply a deep neural network architecture to analyze both SDSS-III APOGEE DR13 and synthetic stellar spectra. When our convolutional neural network model (StarNet) is trained on APOGEE spectra, we show that the stellar parameters (temperature, gravity, and metallicity) are determined with similar precision and accuracy as the APOGEE pipeline. StarNet can also predict stellar parameters when trained on synthetic data, with excellent precision and accuracy for both APOGEE data and synthetic data, over a wide range of signal-to-noise ratios. In addition, the statistical uncertainties in the stellar parameter determinations are comparable to the differences between the APOGEE pipeline results and those determined independently from optical spectra. We compare StarNet to other data-driven methods; for example, StarNet and the Cannon 2 show similar behaviour when trained with the same datasets, however StarNet performs poorly on small training sets like those used by the original Cannon. The influence of the spectral features on the stellar parameters is examined via partial derivatives of the StarNet model results with respect to the input spectra. While StarNet was developed using the APOGEE observed spectra and corresponding ASSET synthetic data, we suggest that this technique is applicable to other wavelength ranges and other spectral surveys.

The 15 273 Å diffuse interstellar band in the dark cloud Barnard 68

High obscuration of background stars behind dark clouds precludes the detection of optical diffuse interstellar bands (DIBs) and hence our knowledge of DIB carriers in these environments. Taking advantage of the reduced obscuration of starlight in the near-infrared (NIR) we used one of the strongest NIR DIBs at 15 273 Å to probe the presence and properties of its carrier throughout the nearby interstellar dark cloud Barnard 68. We measured equivalent widths (EW) for different ranges of visual extinction AV, using VLT/KMOS H-band (1.46–1.85 μm) moderate-resolution (R ~ 4000) spectra of 43 stars situated behind the cloud. To do so, we fitted the data with synthetic stellar spectra from the APOGEE project and TAPAS synthetic telluric transmissions appropriate for the observing site and time period. The results show an increase of DIB EW with increasing AV. However, the rate of increase is much flatter than expected from the EW-AV quasi-proportionality established for this DIB in the Galactic diffuse interstellar medium. Based on a simplified inversion assuming sphericity, it is found that the volume density of the DIB carrier is 2.7 and 7.9 times lower than this expected average value in the external and central regions of the cloud, which have nH≃ 0.4 and 3.5 × 105 cm-3, respectively. Further measurements with multiplex NIR spectrographs should allow detailed modeling of such an edge effect of this DIB and other bands and help clarify its actual origin.

Grids of synthetic stellar spectra

Grids of stellar spectra are extensively used in astrophysics in a range of applications, from derivation of atmospheric parameters in observed stellar spectra, to archaeological studies in integrated light of galaxies. Theoretical stellar libraries have been increasingly used in recent years to overcome limitations of empirical libraries by, for example, exploring atmospheric parameter spaces not well represented in the latter. Nevertheless, there are still caveats and stellar libraries can only partially model observations with accuracy. Here we discuss the current need for synthetic libraries, in particular for applications in stellar population studies. Next we briefly present a non-extensive sample of libraries available in the literature, and we then highlight how better atomic and molecular opacities are much needed to overcome the deficiencies of current synthetic libraries.

The effect of ISM absorption on stellar activity measurements and its relevance for exoplanet studies

Past UV and optical observations of stars hosting hot Jupiters have shown that some of these stars present an anomalously low chromospheric activity, significantly below the basal level. For WASP-13, observations have shown that the apparent lack of activity is possibly caused by absorption from the intervening ISM. Inspired by this result, we study the effect of ISM absorption on activity measurements (S and logR'$_{\rm HK}$ indices) for main-sequence late-type stars. To this end, we employ synthetic stellar photospheric spectra combined with varying amounts of chromospheric emission and ISM absorption. We present the effect of ISM absorption on activity measurements by varying several instrumental, stellar, and ISM parameters. We find that for relative velocities between the stellar and ISM lines smaller than 30-40 km/s and for ISM CaII column densities logN$_{\rm CaII}$>12, the ISM absorption has a significant influence on activity measurements. Direct measurements and three dimensional maps of the Galactic ISM absorption indicate that an ISM CaII column density of logN$_{\rm CaII}$=12 is typically reached by a distance of about 100 pc along most sight lines. In particular, for a Sun-like star lying at a distance greater than 100 pc, we expect a depression (bias) in the logR'$_{\rm HK}$ value larger than 0.05-0.1 dex, about the same size as the typical measurement and calibration uncertainties on this parameter. This work shows that the bias introduced by ISM absorption must always be considered when measuring activity for stars lying beyond 100 pc. We also consider the effect of multiple ISM absorption components. We discuss the relevance of this result for exoplanet studies.

A study on missing lines in the synthetic solar spectrum near the Ca triplet

Synthetic stellar spectra are extensively used for many different applications in astronomy, from stellar studies (such as in the determination of atmospheric parameters of observed stellar spectra), to extragalactic studies (e.g. as one of the main ingredients of stellar population models). One of the main ingredients of synthetic spectral libraries are the atomic and molecular line lists, which contain the data required to model all the absorption lines that should appear in these spectra. Although currently available line lists contain millions of lines, a relatively small fraction of these lines have accurate derived or measured transition parameters. As a consequence, many of these lines contain errors in the electronic transition parameters that can reach up to 200%. Furthermore, even for the Sun, our closest and most studied star, state-of-the-art synthetic spectra does not reproduce all the observed lines, indicating transitions that are missing in the line lists of the computed synthetic spectra. Given the importance and wide range of applications of these models, improvement of their quality is urgently necessary. In this work we catalogued missing lines in the atomic and molecular line lists used for the calculation of the synthetic spectra in the region of GAIA, comparing a solar model computed via a recent line list with a high quality solar atlas available in the literature. After that, we attempted the calibration of their atomic parameters with the code ALLiCE; the calibrated line parameters are publicly available for use.

Measuring diffuse interstellar bands with cool stars

DIBs are ubiquitous in stellar spectra. Traditionally, they have been studied through their extraction from hot stars, because of their smooth continuum. In an era where there are several going-on or planned massive Galactic surveys using multi-object spectrographs, cool stars constitute an appealing set of targets. From the technical point of view, the extraction of DIBs in their spectra is more challenging due to the complexity of the continuum. In this contribution we will provide the community with an improved set of stellar lines in the spectral regions associated to the strong DIBs at l6196, l6269, l6284, and l6379. These lines will allow for the creation of better stellar synthetic spectra, reproducing the background emission and a more accurate extraction of the magnitudes associated with a given DIB. The Sun and Arcturus were used as representative examples of dwarf and giant stars, respectively. A high quality spectrum for each of them was modeled using TURBOSPECTRUM and the VALD stellar line list. The oscillator strength log(gf) and/or wavelength of specific lines were modified to create synthetic spectra where the residuals in both the Sun and Arcturus were minimized. The synthetic spectra based on the improved line lists reproduce the observed spectra for the Sun and Arcturus in the mentioned spectral ranges with greater accuracy. Residuals between the synthetic and observed spectra are always <10%, much better than with previously existing options. The new line list has been tested with some characteristic spectra, from a variety of stars, including both giant and dwarf stars, and under different degrees of extinction. As it happened with the Sun and Arcturus residuals in the fits used to extract the DIB information are smaller when using synthetic spectra made with the updated line lists. Tables with the updated parameters are provided to the community.

COSMIC REIONIZATION ON COMPUTERS. ULTRAVIOLET CONTINUUM SLOPES AND DUST OPACITIES IN HIGH REDSHIFT GALAXIES

ABSTRACT We compare the properties of stellar populations of model galaxies from the Cosmic Reionization On Computers (CROC) project with the exiting ultraviolet (UV) and IR data. Since CROC simulations do not follow cosmic dust directly, we adopt two variants of the dust-follows-metals ansatz to populate model galaxies with dust. Using the dust radiative transfer code Hyperion, we compute synthetic stellar spectra, UV continuum slopes, and IR fluxes for simulated galaxies. We find that the simulation results generally match observational measurements, but, perhaps, not in full detail. The differences seem to indicate that our adopted dust-follows-metals ansatzes are not fully sufficient. While the discrepancies with the exiting data are marginal, the future James Webb Space Telescope (JWST) data will be of much higher precision, rendering highly significant any tentative difference between theory and observations. It is, therefore, likely, that in order to fully utilize the precision of JWST observations, fully dynamical modeling of dust formation, evolution, and destruction may be required.

Stellar activity as noise in exoplanet detection – II. Application to M dwarfs

The ubiquity of M dwarf stars combined with their low masses and luminosities make them prime targets in the search for nearby, habitable exoplanets. We investigate the effects of starspot-induced radial velocity (RV) jitter on detection and characterization of planets orbiting M dwarfs. We create surface spot configurations with both random spot coverage and active regions. Synthetic stellar spectra are calculated from a given spot map, and RV measurements are obtained using cross-correlation technique. We add the RV signal of an orbiting planet to these jitter measurements, and reduce the data to "measure" the planetary parameters. We investigate the detectability of planets around M dwarfs of different activity levels, and the recovery of input planetary parameters. When studying the recovery of the planetary period we note that while our original orbital radius places the planet inside the HZ of its star, even at a filling factor of 2% a few of our measurements fall outside the "conservative Habitable Zone". Higher spot filling factors result in more and higher deviations. Our investigations suggest that caution should be used when characterizing planets discovered with the RV method around stars that are (or are potentially) active.

TheGaia-ESO Survey: Extracting diffuse interstellar bands from cool star spectra

We study how diffuse interstellar bands (DIBs) measured toward distance-distributed target stars can be used to locate dense interstellar (IS) clouds in the Galaxy and probe a line-of-sight (LOS) kinematical structure, a potential useful tool when gaseous absorption lines are saturated or not available in the spectral range. Cool target stars are numerous enough for this purpose. We have devised automated DIB fitting methods appropriate to cool star spectra and multiple IS components. The data is fitted with a combination of a synthetic stellar spectrum, a synthetic telluric transmission, and empirical DIB profiles. In parallel, stellar distances and extinctions are estimated self-consistently by means of a 2D Bayesian method, from spectroscopically-derived stellar parameters and photometric data. We have analyzed Gaia-ESO Survey (GES) and previously recorded spectra that probe between $\sim$ 2 and 10 kpc long LOS in five different regions of the Milky Way. Depending on the observed spectral intervals, we extracted one or more of the following DIBs: $\lambda\lambda$ 6283.8, 6613.6 and 8620.4. For each field, we compared the DIB strengths with the Bayesian distances and extinctions, and the DIB Doppler velocities with the HI emission spectra. For all fields, the DIB strength and the target extinction are well correlated. In case of targets widely distributed in distance, marked steps in DIBs and extinction radial distance profiles match with each other and broadly correspond to the expected locations of spiral arms. For all fields, the DIB velocity structure agrees with HI emission spectra and all detected DIBs correspond to strong NaI lines. This illustrates how DIBs can be used to locate the Galactic interstellar gas and to study its kinematics at the kpc scale.

Chemical abundances in the globular clusters NGC6229 and NGC6779

Long-slit medium-resolution spectra of the Galactic globular clusters (GCs) NGC6229 and NGC6779, obtained with the CARELEC spectrograph at the 1.93-m telescope of the Haute-Provence observatory, have been used to determine the age, helium abundance (Y), and metallicity [Fe/H] as well as the first estimate of the abundances of C, N, O, Mg, Ca, Ti, and Cr for these objects. We solved this task by comparing the observed spectra and the integrated synthetic spectra, calculated with the use of the stellar atmosphere models with the parameters preset for the stars from these clusters. The model mass estimates, T eff, and log g were derived by comparing the observed “color-magnitude” diagrams and the theoretical isochrones. The summing-up of the synthetic blanketed stellar spectra was conducted according to the Chabrier mass function. To test the accuracy of the results, we estimated the chemical abundances, [Fe/H], log t, and Y for the NGC5904 and NGC6254 clusters, which, according to the literature, are considered to be the closest analogues of the two GCs of our study. Using the medium-resolution spectra from the library of Schiavon et al., we obtained for these two clusters a satisfactory agreement with the reported estimates for all the parameters within the errors. We derived the following cluster parameters. NGC6229: [Fe/H] = −1.65 dex, t = 12.6 Gyr, Y = 0.26, [α/Fe] = 0.28 dex; NGC6779: [Fe/H] = −1.9 dex, t = 12.6 Gyr, Y = 0.23, [α/Fe] = 0.08 dex; NGC5904: [Fe/H] = −1.6 dex, t = 12.6 Gyr, Y = 0.30, [α/Fe] = 0.35 dex; NGC6254: [Fe/H] = −1.52 dex, t = 11.2 Gyr, Y = 0.30, [α/Fe] = 0.025 dex. The value [α/Fe] denotes the average of the Ca and Mg abundances.

Impact of micro-telluric lines on precise radial velocities and its correction

Context: In the near future, new instruments such as ESPRESSO will arrive, allowing us to reach a precision in radial-velocity measurements on the order of 10 cm/s. At this level of precision, several noise sources that until now have been outweighed by photon noise will start to contribute significantly to the error budget. The telluric lines that are not neglected by the masks for the radial velocity computation, here called micro-telluric lines, are one such noise source. Aims: In this work we investigate the impact of micro-telluric lines in the radial velocities calculations. We also investigate how to correct the effect of these atmospheric lines on radial velocities. Methods: The work presented here follows two parallel lines. First, we calculated the impact of the micro-telluric lines by multiplying a synthetic solar-like stellar spectrum by synthetic atmospheric spectra and evaluated the effect created by the presence of the telluric lines. Then, we divided HARPS spectra by synthetic atmospheric spectra to correct for its presence on real data and calculated the radial velocity on the corrected spectra. When doing so, one considers two atmospheric models for the synthetic atmospheric spectra: the LBLRTM and TAPAS. Results: We find that the micro-telluric lines can induce an impact on the radial velocities calculation that can already be close to the current precision achieved with HARPS, and so its effect should not be neglected, especially for future instruments such as ESPRESSO. Moreover, we find that the micro-telluric lines' impact depends on factors, such as the radial velocity of the star, airmass, relative humidity, and the barycentric Earth radial velocity projected along the line of sight at the time of the observation.

The effects of spectrograph slit modes on the accuracy of stellar radial velocity measurement and atmospheric parameter estimation

Spectrograph slit is conventionally used to enhance the spectral resolution and manage how much light can be allowed to enter spectrograph. The narrow slit provides a higher resolution but sacrifices efficiency of spectrograph and results in a low signal to noise ratio (S/N) spectra product. We take GuoShouJing telescope as an example and carry out a series of experiments to study how its 2/3 slit mode affects the precision of stellar radial velocity measurement and atmosphere parameters estimate. By transforming the resolution and adding a Gaussian White Noise to the extremely high quality spectra from the Sloan Digital Sky Survey, we generate synthetic stellar spectra of various brightness with different S/Ns. Comparing the measurements on these noise added spectra with the original high quality ones, we summarize the influences of the 2/3 slit mode on the measurement accuracy of stellar radial velocity and atmospheric parameters.

A new extensive library of PHOENIX stellar atmospheres and synthetic spectra

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

On the interpolation of model atmospheres and high-resolution synthetic stellar spectra

We present tests carried out on optical and infrared stellar spectra to evaluate the accuracy of different types of interpolation. Both model atmospheres and continuum normalized fluxes were interpolated. In the first case we used linear interpolation, and in the second linear, cubic spline, cubic-Bezier and quadratic-Bezier methods. We generated 400 ATLAS9 model atmospheres with random values of the atmospheric parameters for these tests, spanning between -2.5 and +0.5 in [Fe/H], from 4500 to 6250 K in effective temperature, and 1.5 to 4.5 dex in surface gravity. Synthesized spectra were created from these model atmospheres, and compared with spectra derived by interpolation. We found that the most accurate interpolation algorithm among those considered in flux space is cubic-Bezier, closely followed by quadratic-Bezier and cubic splines. Linear interpolation of model atmospheres results in errors about a factor of two larger than linear interpolation of fluxes, and about a factor of four larger than high order flux interpolations.

AN ONLINE CATALOG OF CATACLYSMIC VARIABLE SPECTRA FROM THEFAR-ULTRAVIOLET SPECTROSCOPIC EXPLORER

We present an online catalog containing spectra and supporting information for cataclysmic variables that have been observed with the Far Ultraviolet Spectroscopic Explorer (FUSE). For each object in the catalog we list some of the basic system parameters such as (RA,Dec), period, inclination, white dwarf mass, as well as information on the available FUSE spectra: data ID, observation date and time, and exposure time. In addition, we provide parameters needed for the analysis of the FUSE spectra such as the reddening E(B-V), distance, and state (high, low, intermediate) of the system at the time it was observed. For some of these spectra we have carried out model fits to the continuum with synthetic stellar and/or disk spectra using the codes TLUSTY and SYNSPEC. We provide the parameters obtained from these model fits; this includes the white dwarf temperature, gravity, projected rotational velocity and elemental abundances of C, Si, S and N, together with the disk mass accretion rate, the resulting inclination and model-derived distance (when unknown). For each object one or more figures are provided (as gif files) with line identification and model fit(s) when available. The FUSE spectra as well as the synthetic spectra are directly available for download as ascii tables. References are provided for each object as well as for the model fits. In this article we present 36 objects, and additional ones will be added to the online catalog in the future. In addition to cataclysmic variables, we also include a few related objects, such as a wind accreting white dwarf, a pre-cataclysmic variable and some symbiotics.

STELLAR LIFETIME AND ULTRAVIOLET PROPERTIES OF THE OLD METAL-RICH GALACTIC OPEN CLUSTER NGC 6791: A PATHWAY TO UNDERSTAND THE ULTRAVIOLET UPTURN OF ELLIPTICAL GALAXIES

The evolutionary properties of the old metal-rich Galactic open cluster NGC6791 are assessed, based on deep UB photometry and 2Mass JK data. For 4739 stars in the cluster, bolometric luminosity and effective temperature have been derived from theoretical (U-B) and (J-K) color fitting. The derived H-R diagram has been matched with the UVBLUE grid of synthetic stellar spectra to obtain the integrated SED of the system, together with a full set UV (Fanelli) and optical (Lick) narrow-band indices. The cluster appears to be a fairly good proxy of standard elliptical galaxies, although with significantly bluer infrared colors, a shallower 4000A Balmer break, and a lower Mg2 index. The confirmed presence of a dozen hot stars, along their EHB evolution, leads the cluster SED to consistently match the properties of the most active UV-upturn galaxies, with 1.7+/-0.4% of the total bolometric luminosity emitted shortward of 2500A. The cluster Helium abundance results Y=0.30 +/-0.04, while the Post-MS implied stellar lifetime from star number counts fairly agrees with the theoretical expectations from both the Padova and BASTI stellar tracks. A Post-MS fuel consumption of 0.43 +/- 0.01 M_sun is found for NGC6791 stars, in close agreement with the estimated mass of cluster He-rich white dwarfs. Such a tight figure may lead to suspect that a fraction of the cluster stellar population does actually not reach the minimum mass required to effectively ignite He in the stellar core.