Aims. Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundred objects from the southern hemisphere (UVES-POP). We aim to extend the sample, offering a finer coverage of effective temperatures and surface gravity with a uniform collection of spectra obtained in the northern hemisphere. Methods. Between 2010 and 2020, we acquired several thousand echelle spectra of bright stars with the Mercator-HERMES spectrograph located in the Roque de Los Muchachos Observatory in La Palma, whose pipeline offers high-quality data reduction products. We have also developed methods to correct for the instrumental response in order to approach the true shape of the spectral continuum. Additionally, we have devised a normalisation process to provide a homogeneous normalisation of the full spectral range for most of the objects. Results. We present a new spectral library consisting of 3256 spectra covering 2043 stars. It combines high signal-to-noise and high spectral resolution over the entire range of effective temperatures and luminosity classes. The spectra are presented in four versions: raw, corrected from the instrumental response, with and without correction from the atmospheric molecular absorption, and normalised (including the telluric correction).
Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundred objects from the southern hemisphere (UVES-POP) . We aim to extend the sample, offering a finer coverage of effective temperatures and surface gravity with a uniform collection of spectra obtained in the northern hemisphere. Between 2010 and 2020, we acquired several thousand echelle spectra of bright stars with the Mercator-HERMES spectrograph located in the Roque de Los Muchachos Observatory in La Palma, whose pipeline offers high-quality data reduction products. We have also developed methods to correct for the instrumental response in order to approach the true shape of the spectral continuum. Additionally, we have devised a normalisation process to provide a homogeneous normalisation of the full spectral range for most of the objects. We present a new spectral library consisting of 3256 spectra covering 2043 stars. It combines high signal-to-noise and high spectral resolution over the entire range of effective temperatures and luminosity classes. The spectra are presented in four versions: raw, corrected from the instrumental response, with and without correction from the atmospheric molecular absorption, and normalised (including the telluric correction).
Context. Protoplanetary disks show large diversity regarding their morphology and dust composition. With mid-infrared interferometry the thermal emission of disks can be spatially resolved, and the distribution and properties of the dust within can be studied. Aims. Our aim is to perform a statistical analysis on a large sample of 82 disks around low- and intermediate-mass young stars, based on mid-infrared interferometric observations. We intend to study the distribution of disk sizes, variability, and the silicate dust mineralogy. Methods. Archival mid-infrared interferometric data from the MIDI instrument on the Very Large Telescope Interferometer are homogeneously reduced and calibrated. Geometric disk models are used to fit the observations to get spatial information about the disks. An automatic spectral decomposition pipeline is applied to analyze the shape of the silicate feature. Results. We present the resulting data products in the form of an atlas, containing N band correlated and total spectra, visibilities, and differential phases. The majority of our data can be well fitted with a continuous disk model, except for a few objects, where a gapped model gives a better match. From the mid-infrared size-luminosity relation we find that disks around T Tauri stars are generally colder and more extended with respect to the stellar luminosity than disks around Herbig Ae stars. We find that in the innermost part of the disks (r less than or similar to 1 au) the silicate feature is generally weaker than in the outer parts, suggesting that in the inner parts the dust is substantially more processed. We analyze stellar multiplicity and find that in two systems (AB Aur and HD 72106) data suggest a new companion or asymmetric inner disk structure. We make predictions for the observability of our objects with the upcoming Multi-AperTure mid-Infrared Spectro Scopic Experiment (MATISSE) instrument, supporting the practical preparations of future MATISSE observations of T Tauri stars.
Context. DG Tau is a low-mass pre-main sequence star, whose strongly accreting protoplanetary disk exhibits a so-far enigmatic behavior: its mid-infrared thermal emission is strongly time-variable, even turning the 10 mu m silicate feature from emission to absorption temporarily.Aims. We look for the reason for the spectral variability at high spatial resolution and at multiple epochs.Methods. Infrared interferometry can spatially resolve the thermal emission of the circumstellar disk, also giving information about dust processing. We study the temporal variability of the mid-infrared interferometric signal, observed with the VLTI/MIDI instrument at six epochs between 2011 and 2014. We fit a geometric disk model to the observed interferometric signal to obtain spatial information about the disk. We also model the mid-infrared spectra by template fitting to characterize the profile and time dependence of the silicate emission. We use physically motivated radiative transfer modeling to interpret the mid-infrared interferometric spectra.Results. The inner disk (r < 1-3 au) spectra exhibit a 10 mu m absorption feature related to amorphous silicate grains. The outer disk (r > 1-3 au) spectra show a crystalline silicate feature in emission, similar to the spectra of comet Hale-Bopp. The striking difference between the inner and outer disk spectral feature is highly unusual among T Tauri stars. The mid-infrared variability is dominated by the outer disk. The strength of the silicate feature changed by more than a factor of two. Between 2011 and 2014 the half-light radius of the mid-infrared-emitting region decreased from 1.15 to 0.7 au.Conclusions. For the origin of the absorption we discuss four possible explanations: a cold obscuring envelope, an accretion heated inner disk, a temperature inversion on the disk surface and a misaligned inner geometry. The silicate emission in the outer disk can be explained by dusty material high above the disk plane, whose mass can change with time, possibly due to turbulence in the disk.
We present scattered light images of the TW Hya disk performed with the Spectro-Polarimetric High-contrast Exoplanet REsearch instrument in Polarimetric Differential Imaging mode at 0.63, 0.79, 1.24, and 1.62 mu m. We also present H2/H3-band angular differential imaging (ADI) observations. Three distinct radial depressions in the polarized intensity distribution are seen, around approximate to 85, approximate to 21, and less than or similar to 6 au.(21) The overall intensity distribution has a high degree of azimuthal symmetry; the disk is somewhat brighter than average toward the south and darker toward the north-west. The ADI observations yielded no signifiant detection of point sources in the disk. Our observations have a linear spatial resolution of 1-2 au, similar to that of recent ALMA dust continuum observations. The submicron- sized dust grains that dominate the light scattering in the disk surface are strongly coupled to the gas. We created a radiative transfer disk model with self-consistent temperature and vertical structure iteration and including grain size-dependent dust settling. This method may provide independent constraints on the gas distribution at higher spatial resolution than is feasible with ALMA gas line observations. We find that the gas surface density in the "gaps" is reduced by approximate to 50% to approximate to 80% relative to an unperturbed model. Should embedded planets be responsible for carving the gaps then their masses are at most a few 10M(circle plus). The observed gaps are wider, with shallower flanks, than expected for planet-disk interaction with such low-mass planets. If forming planetary bodies have undergone collapse and are in the "detached phase," then they may be directly observable with future facilities such as the Mid-Infrared E-ELT Imager and Spectrograph at the E-ELT.
We present a mid-IR interferometric survey of the circumstellar environment of a specific class of post-Asymptotic Giant Branch (post-AGB) binaries. For this class the presence of a compact dusty disk has been postulated on the basis of various spatially unresolved measurements. Our interferometric survey was performed with the MIDI instrument on the VLTI. In total 19 different systems were observed using variable baseline configurations. Combining all the visibilities at a single wavelength at 10.7 micron, we fitted two parametric models to the data: a uniform disk (UD) and a ring model mimicking a temperature gradient. We compared our observables of the whole sample, with synthetic data computed from a grid of radiative transfer models of passively irradiated disks in hydrostatic equilibrium. These models are computed with a Monte Carlo code that has been widely applied to describe the structure of protoplanetary disks around young stellar objects (YSO). The spatially resolved observations show that the majority of our targets cluster closely together in the distance-independent size-colour diagram, and have extremely compact N-band emission regions. The typical uniform disk diameter of the N-band emission region is about 40 mass which corresponds to a typical brightness temperature of 400-600 K. The resolved objects display very similar characteristics in the interferometric observables and in the spectral energy distributions. Therefore, the physical properties of the disks around our targets must be similar. The grid of protoplanetary disk models covers very well the observed objects. Much like for young stars, the spatially resolved N-band emission region is determined by the hot inner rim of the disk. Continued comparisons between post-AGB and protoplanetary disks will help to understand grain growth and disk evolution processes,
Context. Many post-asymptotic giant branch (post-AGB) stars in binary systems have an infrared (IR) excess arising from a dusty circumbinary disk. The disk formation, current structure, and further evolution are, however, poorly understood. Aims. We aim to constrain the structure of the circumstellar material around the post-AGB binary and RV Tauri pulsator AC Her. We want to constrain the spatial distribution of the amorphous and of the crystalline dust. Methods. We present very high-quality mid-IR interferometric data that were obtained with the MIDI/VLTI instrument. We analyze the MIDI visibilities and differential phases in combination with the full spectral energy distribution, using the MCMax radiative transfer code, to find a good structure model of AC Her’s circumbinary disk. We include a grain size distribution and midplane settling of dust self-consistently in our models. The spatial distribution of crystalline forsterite in the disk is investigated with the mid-IR features, the 69 μm band and the 11.3 μm signatures in the interferometric data. Results. All the data are well fitted by our best model. The inclination and position angle of the disk are precisely determined at i = 50 ± 8° and PA = 305 ± 10°. We firmly establish that the inner disk radius is about an order of magnitude larger than the dust sublimation radius. The best-fit dust grain size distribution shows that significant grain growth has occurred, with a significant amount of mm-sized grains now being settled to the midplane of the disk. A large total dust mass ≥10-3 M⊙ is needed to fit the mm fluxes. By assuming αturb = 0.01, a good fit is obtained with a small grain size power law index of 3.25, combined with a small gas/dust ratio ≤10. The resulting gas mass is compatible with recent estimates employing direct gas diagnostics. The spatial distribution of the forsterite is different from the amorphous dust, as more warm forsterite is needed in the surface layers of the inner disk. Conclusions. The disk in the AC Her system is in a very evolved state, as shown by its small gas/dust ratio and large inner hole. Mid-IR interferometry offers unique constraints, complementary to mid-IR features, for studying the mineralogy in disks. A better uv coverage is needed to constrain in detail the distribution of the crystalline forsterite in the disk of AC Her, but we find strong similarities with the protoplanetary disk HD 100546.
The disks around Herbig Ae/Be stars are commonly divided into group I and group II based on their far-infrared spectral energy distribution, and the common interpretation for that is flared and flat disks. Recent observations suggest that many flaring disks have gaps, whereas flat disks are thought to be gapless. The different groups of objects can be expected to have different structural signatures in high-angular-resolution data. Over the past 10 years, the MIDI instrument on the Very Large Telescope Interferometer has collected observations of several tens of protoplanetary disks. We model the large set of observations with simple geometric models. A population of radiative-transfer models is synthesized for interpreting the mid-infrared signatures. Objects with similar luminosities show very different disk sizes in the mid-infrared. Restricting to the young objects of intermediate mass, we confirm that most group I disks are in agreement with being transitional. We find that several group II objects have mid-infrared sizes and colors overlapping with sources classified as group I, transition disks. This suggests that these sources have gaps, which has been demonstrated for a subset of them. This may point to an intermediate population between gapless and transition disks. Flat disks with gaps are most likely descendants of flat disks without gaps. Gaps, potentially related to the formation of massive bodies, may therefore even develop in disks in a far stage of grain growth and settling. The evolutionary implications of this new population could be twofold. Either gapped flat disks form a separate population of evolved disks, or some of them may further evolve into flaring disks with large gaps. The latter transformation may be governed by the interaction with a massive planet, carving a large gap and dynamically exciting the grain population in the disk.
Gamma Doradus and delta Scuti pulsators cover the transition region between low mass and massive main-sequence stars, and as such, are critical for testing stellar models. When they reside in binary systems, we can combine two independent methods to derive critical information, such as precise fundamental parameters to aid asteroseismic modelling. In the Kepler light curve of KIC10080943, clear signatures of gravity- and pressure-mode pulsations have been found. Ground-based spectroscopy revealed this target to be a double-lined binary system. We present the analysis of four years of Kepler photometry and high-resolution spectroscopy to derive observational constraints with which to evaluate theoretical predictions of the stellar structure and evolution for intermediate-mass stars. We used the method of spectral disentangling to determine atmospheric parameters for both components and derive the orbital elements. With PHOEBE, we modelled the ellipsoidal variation and reflection signal of the binary in the light curve and used classical Fourier techniques to analyse the pulsation modes. We show that the eccentric binary system KIC10080943 contains two hybrid pulsators with masses $M_1=2.0\pm0.1~M_\odot$ and $M_2=1.9\pm0.1~M_\odot$, with radii $R_1=2.9\pm0.1~R_\odot$ and $R_2=2.1\pm0.2~R_\odot$. We detect rotational splitting in the g modes and p modes for both stars and use them to determine a first rough estimate of the core-to-surface rotation rates for the two components, which will be improved by future detailed seismic modelling.
We present our approach of visibility modeling of disks around low-mass (< 2M circle dot) young stellar objects (YSOs). We compiled an atlas based on mid-infrared interferometric ob servations from the MIDI instrument at the VLTI. We use three different models to fit the data. These models allow us to determine overall sizes (and the extent of the inner gaps) of the modeled circumstellar disks.
Binaries with circumbinary disks are commonly found among optically bright post-AGB stars. Although clearly linked to binary interaction processes, the formation, evolution and fate of these disks are still badly understood. Due to their compactness, interferometric techniques are required to resolve them. Here, we discuss our high-quality multiwavelength interferometric data of two prototypical yet very different post-AGB binaries, AC and 89 Herculis, as well as the modeling thereof with radiative transfer models. A detailed account of the data and models of both objects is published in three separate papers elsewhere; here we focus on comparing the modeling results for the two objects. In particular we discuss the successes and limitations of the models which were developed for protoplanetary disks around young stars. We conclude that multiwavelength high-angular-resolution observations and radiative transfer disk models are indispensible to understand these complex interacting objects and their place in the grand scheme of the (binary) evolution of low and intermediate mass stars.
Abstract We carried out an extensive photometric and spectroscopic investigation of the SPB binary, HD 25558 (see Fig. 1 for the time and geographic distribution of the observations). The ~2000 spectra obtained at 13 observatories during 5 observing seasons, the ground-based multi-colour light curves and the photometric data from the MOST satellite revealed that this object is a double-lined spectroscopic binary with a very long orbital period of about 9 years. We determined the physical parameters of the components, and have found that both lie within the SPB instability strip. Accordingly, both components show line-profile variations consistent with stellar pulsations. Altogether, 11 independent frequencies and one harmonic frequency were identified in the data. The observational data do not allow the inference of a reliable orbital solution, thus, disentangling cannot be performed on the spectra. Since the lines of the two components are never completely separated, the analysis is very complicated. Nevertheless, pixel-by-pixel variability analysis of the cross-correlated line profiles was successful, and we were able to attribute all the frequencies to the primary or secondary component. Spectroscopic and photometric mode-identification was also performed for several of these frequencies of both binary components. The spectroscopic mode-identification results suggest that the inclination and rotation of the two components are rather different. While the primary is a slow rotator with ~6 d rotation period, seen at ~60° inclination, the secondary rotates fast with ~1.2 d rotation period, and is seen at ~20° inclination. Our spectropolarimetric measurements revealed that the secondary component has a magnetic field with at least a few hundred Gauss strength, while no magnetic field was detected in the primary. The detailed analysis and results of this study will be published elsewhere.
The presence of disks and outflows is widespread among post-AGB binaries. In the first paper of this series, a surprisingly large fraction of optical light was found to be resolved in the 89 Her post-AGB system. The data showed this flux to arise from close to the central binary. Scattering off the inner rim of the circumbinary disk, or in a dusty outflow were suggested as two possible origins. With detailed dust radiative transfer models of the disk we aim to discriminate between these two configurations. By including Herschel/SPIRE photometry, we extend the SED such that it now fully covers UV to sub-mm wavelengths. The MCMax radiative transfer code is used to create a large grid of disk models. Our models include a self-consistent treatment of dust settling as well as of scattering. A Si-rich composition with two additional opacity sources, metallic Fe or amorphous C, are tested. The SED is fit together with mid-IR (MIDI) visibilities as well as the optical and near-IR visibilities of Paper I, to constrain the structure of the disk and in particular of its inner rim. The near-IR visibility data require a smooth inner rim, here obtained with a two-power-law parameterization of the radial surface density distribution. A model can be found that fits all the IR photometric and interferometric data well, with either of the two continuum opacity sources. Our best-fit passive models are characterized by a significant amount of mm-sized grains, which are settled to the midplane of the disk. Not a single disk model fits our data at optical wavelengths though, the reason being the opposing constraints imposed by the optical and near-IR interferometric data. A geometry in which a passive, dusty, and puffed-up circumbinary disk is present, can reproduce all the IR but not the optical observations of 89 Her. Another dusty, outflow or halo, component therefore needs to be added to the system.
Context. The presence of stable disks around post-asymptotic giant branch (post-AGB) binaries is a widespread phenomenon. Also, the presence of (molecular) outflows is now commonly inferred in these systems.
To enable optical long baseline interferometry toward faint objects, long integrations are necessary despite atmospheric turbulence. Fringe trackers are needed to stabilize the fringes and thus increase the fringe visibility and phase signal-to-noise ratio (S/N), with efficient controllers robust to instrumental vibrations and to subsequent path fluctuations and flux drop-outs. We report on simulations, analysis, and comparison of the performances of a classical integrator controller and of a Kalman controller, both optimized to track fringes under realistic observing conditions for different source magnitudes, disturbance conditions, and sampling frequencies. The key parameters of our simulations (instrument photometric performance, detection noise, turbulence, and vibrations statistics) are based on typical observing conditions at the Very Large Telescope observatory and on the design of the GRAVITY instrument, a 4-telescope single-mode long-baseline interferometer in the near-infrared, next in line to be installed at VLT Interferometer. We find that both controller performances follow a two-regime law with the star magnitude, a constant disturbance limited regime, and a diverging regime limited by both the detector and photon-noise. Moreover, we find that the Kalman controller is optimal in the high and medium S/N regime owing to its predictive commands based on an accurate disturbance model. In the low S/N regime, the model is not accurate enough to be more robust than an integrator controller. Identifying the disturbances from high S/N measurements improves the Kalman performances in the case of strong optical path difference disturbances.
For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner 3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.
Abstract For over a decade, the structure of the inner “hole” in the transition disk around TW Hydrae has been a subject of debate. To probe the innermost regions of the protoplanetary disk, observations at the highest possible spatial resolution are required. We present new interferometric data of TW Hya from near-infrared to millimeter wavelengths. We confront existing models of the disk structure with the complete data set and develop a new, detailed radiative-transfer model. This model is characterized by: 1) a spatial separation of the largest grains from the small disk grains; and 2) a smooth inner rim structure, rather than a sharp disk edge.
Context. The implementation of fringe tracking for optical interferometers is inevitable when optimal exploitation of the instrumental capacities is desired. Fringe tracking allows continuous fringe observation, considerably increasing the sensitivity of the interferometric system. In addition to the correction of atmospheric path-length differences, a decent control algorithm should correct for disturbances introduced by instrumental vibrations, and deal with other errors propagating in the optical trains.Aims. In an effort to improve upon existing fringe-tracking control, especially with respect to vibrations, we attempt to construct control schemes based on Kalman filters. Kalman filtering is an optimal data processing algorithm for tracking and correcting a system on which observations are performed. As a direct application, control schemes are designed for GRAVITY, a future four-telescope near-infrared beam combiner for the Very Large Telescope Interferometer (VLTI).Methods. We base our study on recent work in adaptive-optics control. The technique is to describe perturbations of fringe phases in terms of an a priori model. The model allows us to optimize the tracking of fringes, in that it is adapted to the prevailing perturbations. Since the model is of a parametric nature, a parameter identification needs to be included. Different possibilities exist to generalize to the four-telescope fringe tracking that is useful for GRAVITY.Results. On the basis of a two-telescope Kalman-filtering control algorithm, a set of two properly working control algorithms for four-telescope fringe tracking is constructed. The control schemes are designed to take into account flux problems and low-signal baselines. First simulations of the fringe-tracking process indicate that the defined schemes meet the requirements for GRAVITY and allow us to distinguish in performance. In a future paper, we will compare the performances of classical fringe tracking to our Kalman-filter control.Conclusions. Kalman-filter based control schemes will likely become the next standard for fringe tracking, providing the possibility to maximize the tracking performance. The results of the present study are currently being incorporated into the final design of GRAVITY.
We report on the discovery of an infrared emission band in the Spitzer spectrum of the S-type AGB star NP Aurigae that is caused by TiO molecules in the circumstellar environment. We modeled the observed emission to derive the temperature of the TiO molecules (approximate to 600K), an upper limit on the column density (approximate to 10(17.25) cm(-2)) and a lower limit on the spatial extent of the layer that contains these molecules. (approximate to 4.6 R-*). This is the first time that this TiO emission band is observed. A search for similar emission features in the sample of S-type stars yielded two additional candidates. However, owing to the additional dust emission, the identification is less stringent. By comparing the stellar characteristics of NP Aur to those of the other stars in our sample, we find that all stars with TiO emission show large-amplitude pulsations, s-process enrichment, and a low C/O ratio. These characteristics might be necessary requirements for a star to show TiO in emission, but they are not sufficient.