On the evolutionary status of high-latitude
variable V534 Lyr
Based on the high resolution spectral monitoring conducted at the 6-m
BTA telescope, we study the optical spectrum of the high-latitude variable V534 Lyr.
Heliocentric radial velocities corresponding to the positions of all metal
absorption components, as well as the Na I D and H lines were measured
during all the observational dates. The analysis of the velocity field examining the
lines of various nature revealed a low-amplitude variability of based on the
lines with a high excitation potential, which are formed in deep layers of the
stellar atmosphere, and allowed to estimate the systemic velocity of km s
( km s). The distance estimate of kpc for
the star leads us to its absolute magnitude of , what corresponds
to the MK spectral classification. The previously undetected for this star spectral phenomenon
was revealed: at certain times a splitting of the profiles of low-excited absorptions
is observed, reaching – km s. A combination of the
parameters: reduced metallicity , increased nitrogen abundance
, large spatial velocity, high luminosity, a strong
variability of the emission-absorption profiles of H I lines, splitting of
metal absorptions at different times of observations and the variability of
the velocity field in the atmosphere allow us to classify V534 Lyr as a
pulsating star near the HB and belonging to the thick disk of our Galaxy.
Keywords: stars, evolution, post-AGB stars, pulsating stars, optical spectroscopy.
During the past two decades, the 6-meter telescope BTA hosts an investigation of the spectra of numerous supergiants with a large excess of IR-radiation, presumably belonging to the stage of the asymptotic giant branch (AGB) and post-AGB. The program also includes a study of some related stars with an unclear evolutionary status. At the post-AGB stage, far evolved intermediate-mass stars are observed, with initial masses at the main sequence . Before the AGB stage, these stars are cool (with effective temperature of – K) red supergiants. After they have exhausted helium, a degenerate carbon-oxygen core is formed in the core of AGB stars with an initial mass between 2–4 . This core is surrounded by alternately energetically active layers of He and H combustion. Most of the time, the hydrogen layer provides the energy release, but as the products of hydrogen burning join the helium layer, a brief ignition of helium occurs in the latter. This configuration of the internal structure of the star is unstable, the theory predicts a fairly effective mixing and dredge-up (due to the penetrating convection) of heavy metals synthesized in nuclear reactions, accompanying these processes of energy release into the atmosphere of the star (see the survey [Herwig ()] and references therein).
Our initial goal was to determine the fundamental parameters of investigated stars and to search for the chemical composition anomalies in their atmospheres, related with the synthesis of chemical elements during the previous evolutionary stages. As the program runs, additional research became crucial, aimed to search for the time variability of spectral features and the velocity fields in the atmospheres and envelopes of studied stars. This expansion of the task required conducting for each program object of multiple observations with high spectral resolution in a broad wavelength range. The main results we have obtained from the spectroscopic data at the 6-m telescope are summarized in a number of recent publications [Envelopes (); Enrich (); Conf ()]. Let us briefly note the most important of them. Firstly, the parameters and chemical composition of the atmospheres for several tens of stars with an excess in infrared flux were determined. Secondly, evolutionary variations in chemical composition (a large excess of the -process or hot-bottom-process elements) were found for seven stars [04296 (); 23304 (); Egg1 (); 20000 (); 20572 ()]. Based on the representative sample of stars, studied by us and the authors of [Winckel ()], we formed a subsample of post-AGB stars with the atmospheres rich in the -process carbon and heavy metals, and with the carbon-rich circumstellar envelopes. The analysis of properties of this subsample led us to the conclusion on the interrelation of the peculiarity of line profiles, manifested in the presence of the emission component of the Na I doublet D lines, in the character of molecular features, asymmetry, permitted or forbidden metal emissions, formed in the envelopes, and in the splitting of strong absorption profiles with a low lower-level excitation potential, with kinematic and chemical properties of the circumstellar envelope and with the type of its morphology [Envelopes ()].
The presence of the above-mentioned features is the main difference of the spectra of post-AGB stars from the spectra of massive supergiants. The splitting into the components of the profiles of strongest heavy-metal absorptions in the spectra of selected post-AGB supergiants we have detected allowed us to conclude that the process of formation of a structured circumstellar envelope is accompanied with an enrichment by the products of stellar nucleosynthesis [Enrich ()]. Variability of the observed profiles of the absorption-emission H line and metal lines, as well as the type change (absorption/emission) of the C molecule Swan bands registered in several objects are caused by the variations of the parameters or the structure of the circumstellar envelope. The H profile type (pure absorption, pure emission, P Cyg or inverse P Cyg-type, with two emission components in the wings) is not related to the chemical composition of the atmosphere of the central star. The main factors influencing the H profile type and its variability are the mass loss rate, the velocity of stellar wind, kinematics and optical thickness of the envelope.
We have discovered the peculiarities, previously unknown in the optical spectra of the AGB-candidate stars. In particular, based on the observations, also performed with the NES echelle spectrograph of the 6-m telescope, the spectral features and the velocity fields in the atmosphere and the envelope of the cool AGB-supergiant V1027 Cyg [V1027Cyg ()], the optical component of the IR source IRAS 20004+2955 were studied. For the first time, the spectrum of this star revealed the splitting of the cores of strong metal absorptions and their ions (Si II, Ni I, Ti I, Ti II, Sc II, Cr I, Fe I, Fe II, Ba II). The broad profile of these lines contains a stably located weak emission in the core, the position of which can be considered as the velocity of the center of mass of the system (for brevity, we shall further refer to it as the systemic velocity) km s. From the symmetric small and moderate-intensity absorptions, weak radial velocity variations were revealed with an amplitude of 5–6 km s, caused by pulsations. A red-wave H profile shift was observed, caused by the distortion of the line core. In the red spectral region, numerous weak lines of the CN molecule and the K I 7696 Å line, which has a P Cyg-type profile are identified. An agreement of radial velocity, measured from the symmetric absorptions of metals and the CN lines indicates formation of a CN spectrum in the atmosphere of the star. Numerous interstellar bands, DIBs, were identified, the position of which in the spectrum, km s, corresponds to the velocity of the interstellar medium in the Local Arm of the Galaxy.
Absolutely new results were obtained for a nameless faint star [RAFGL5081 ()], associated with a powerful source of IR radiation RAFGL 5081. Its optical spectrum has been studied for the first time based on the long-term spectral monitoring with high spectral resolution. The estimates of the spectral class of the star Sp G5–8 II and its effective temperature K were made. We detected a hard-to-explain spectral phenomenon: a splitting of the medium and low-intensity profiles of stationary absorptions. Stationarity of absorptions eliminates the possibility of explaining the double peaks owing to the spectral binarity of the star. Radial velocities for the wind components of the Na I D lines and H line profiles reach km s and km s respectively. These profiles contain narrow components, the number, the depths and positions of which vary with time. The time-variable multicomponent structure of the Na I D lines and H line profiles indicates the heterogeneity and instability of the circumstellar envelope of RAFGL 5081. According to the presence in the Na I (1) line profiles of the components with the velocity of km s, it is concluded that RAFGL 5081 is located behind the Perseus Arm, i.e. not closer than 2 kpc. Note that this object is associated with a reflective nebula GN 02.44.7. The problem of photometric monitoring for determining the variability parameters of the central star RAFGL 5081 is set.
In the study of sources with a large IR flux excess, we identified the stars outside the expected AGB and post-AGB stages. It turned out that among such objects, in addition to AGB- and post-AGB stars, there exist massive stars of extremely high luminosity with extended and structured circumstellar envelopes. The most famous object that has for many years been preferably considered as a star in the post-AGB stage is V1302 Aql, which has numerous peculiarities in the spectrum. For a long time, the evolutionary status of this supergiant, associated with a powerful IR source IRC+10420, was unclear. The set of observed properties of V1302 Aql allowed to consider it as a star at the post-AGB stage or as a very massive star, past the red supergiant stage. It is obvious that depending on the accepted status, and hence on the luminosity of the object too, the estimate of its distance from the observer may differ several times. One of the decisive arguments, confirming the status of a massive evolved star for V1302 Aql was obtained by us in the analysis of spectral data at the 6-m telescope when the chemical composition of its atmosphere was first determined and a significant excess of nitrogen was found [IRC1 ()]. Detection of a rapid growth of the effective temperature [IRC1 (); Oudm98 ()] allowed to assume that the star quickly evolves with the increase rate of about 120 K a year. The phenomenon was for us an incentive to continue the spectral monitoring of this mysterious object. In the study of the V1302 Aql spectra over the period of observations of 1992–2014 the conclusion was drawn that the yellow hypergiant entered the phase of deceleration (or the suspension) of growth of the effective temperature and approached at the H–R diagram the high-temperature boundary of the Yellow Void [IRC3 ()]. It is appropriate to mention here the star HD 179821, which has for several decades appeared in the papers as a post-AGB candidate, and only a detailed analysis of the high-quality data of a long-term spectral monitoring with the echelle spectrographs of the SAO RAS 6-m BTA telescope and the McDonald observatory 2.7-m telescope (USA) led to the conclusion that it belonged to a family of massive evolved stars [HD179821 ()].
The experience obtained in the analysis of the spectra of stars of different luminosity and mass allows us to hope for a success in the study of stars with an obscure evolutionary status. The object of the present study is the star V534 Lyr (HD 172324, HIP 91359), located at the galactic latitude . Our interest to this relatively bright () star is due to the fact that Trams et al. [Trams1991 ()] included V534 Lyr in the list of high-luminosity stars at high Galactic latitudes (HGL). These authors attributed the major part of investigated objects to the post-AGB candidates. However, due to the lack of data on the IR flux, they classified V534 Lyr only as an alleged candidate in post-AGB. An additional stimulus to the study of V534 Lyr for us is the emission in the H I lines, found in the early paper by Bonsack and Greenstein [Bonsack1956 ()] and confirmed later [Giridhar2001 ()]. This star was repeatedly studied using various photometric systems [Rufener1976 (); Fernie1983 (); Hauck1998 ()] and spectral methods, but up to now, none of the available publications has a definite conclusion on its evolutionary status (see Table 1, in which the basic papers examining V534 Lyr are listed [Veen (); Giridhar2001 (); Giridhar2005 ()]).
|IR Photometry||Lack of IR excess||Post-AGB star?||[Veen ()]|
|Optical spectroscopy,||, ,||Post-AGB star ?||[Giridhar2001 ()]|
|method of model atmospheres||chemical composition|
|Optical Spectroscopy,||Variability and||High Speed|
|method of model atmospheres||H I profiles, chemical composition||HGL star ?||[Giridhar2005 ()]|
The purpose of this paper is to determine the main parameters of V534 Lyr using the high-quality optical spectral data, to search for the possible variability in time of the spectrum and the pattern of radial velocities in the atmosphere of the star, and to refine its evolutionary status. Section 2 briefly describes the used observational data and details of its processing. The data on the spectral features obtained from the high-resolution spectra and the detected profile features are given in Section 3. We make an attempt here to apply our data for the determination of luminosity, remoteness and evolutionary status of the star, as well as present a method for determining the parameters of the stellar atmosphere and calculating the abundances of chemical elements. Section 4 discusses the results of our study. The main conclusions are summarized in the final Section 5.
2 Observations and processing of spectra
In this study we used seven high-resolution spectra of V534 Lyr (), obtained with the NES echelle spectrograph [nes2 ()] of the 6-m BTA telescope on arbitrary dates in 2010 and 2017, as well as the earliest of our spectra, which was obtained on July 8, 2000 with the PFES echelle spectrograph installed in the primary focus of the 6-m telescope with a lower resolution of [PFES1 (); PFES2 ()], which imposed some restrictions on its use. Extraction of one-dimensional spectra from two-dimensional echelle frames was executed using a modified version [Yushkin ()] of the ECHELLE context from the MIDAS software suite. The removal of traces of cosmic particles was done by the median averaging of two spectra obtained successively one after another. The wavelength calibration was carried out using the spectra of the hollow-cathode Th-Ar lamp. For the subsequent spectrophotometric and positional reduction of one-dimensional spectra we used a modified version of the DECH20t code [dech ()]. The control of instrumental matching of the spectra of the star and the hollow-cathode lamp is made from the telluric lines [O I], O and HO.
In more detail, the procedure for measuring the radial velocity based on the spectra obtained with the NES spectrograph, and sources of error are described in [VR ()]. The mean square measurement error for stars with narrow absorptions in the spectrum does not exceed 1.0 km s (one-line accuracy [VR ()]). The initial data concerning the obtained spectra are presented in the first three columns of Table 2.
|Date||UT||, nm||, K||, km s|
|Fe II||He I+S II+Si II||Fe II (high)||Fe II (low)|
|–the average value based on Fe II (42) lines,|
|–the average value based on Fe II lines (multiplets 48, 49, 74).|
3 Main results
3.1 Peculiarities of the V534 Lyr spectrum
A peculiar feature of the optical spectrum of V534 Lyr is a powerful emission of the H profile, which was noted in all the papers dealing with the spectroscopy of this star (see, for example, [Giridhar2001 (); Giridhar2005 ()] and references therein). The H line, which is present in seven of eight our spectra, is the most intense. The variations of its profile with a two-peak emission are significant. In the earliest spectrum (July 8, 2000) the entire emission profile is located above the continuum. However, in the subsequent spectra the core of the absorption component which is formed in the surface layers of the atmosphere falls below the level of the continuum. The position of the absorption component, as well as the intensities and positions of the emission components vary, but the intensity of the long-wave peak is always higher than that of the short-wave peak (Fig. 1). The shift of the blue border of the absorption component relative to the vertical dashed line, corresponding to increases in the spectrum of September 24, 2010: km s.
A drastic profile variation is observed in the 2017 spectrum. Instead of the two-peak profile we observe a P Cyg-type profile, the absorption component of which is shifted by km s. The H profile is registered in a smaller number of dates, but this is enough to notice how the stronger of the emission components passes from the red wing to the blue wing (Fig. 2). Similar variations of H and H profiles in the spectrum of V534 Lyr were found in the paper [Giridhar2005 ()]. As it can be seen in Fig. 2, in the spectrum obtained on June 8, 2017, there was a change in the H profile type to the P Cyg-type with the maximum shift of the absorption component. A few days later, on June 13, 2017, we can see the splitting of the H absorption component in the absence of H profile variation.
Note that the V534 Lyr spectrum is characterized by a strengthening of neutral helium He I lines, with an appreciable weakening of metal lines. This is well illustrated in the upper panel of Fig. 3, where a fragment of the spectrum of this star with the He I 4026 Å line is compared with an analogous spectral region of a massive supergiant Cyg. The MK-spectral class of Cyg, namely A2 Ia, is close to V534 Lyr, but the intensities of the He I 4026 Å line are markedly different.
3.2 Radial velocity pattern
Generalized results of radial velocity measurements for each time of observations are presented in Table 2. Columns (5)–(8) contain the values obtained by averaging the velocities for: the Fe II 6318, 6384 and 6385 Å emission, in the 2017 spectra the Fe II 6493, 7496, 7513 Å, emissions, the He I, S II absorption cores, and the doublet Si II (2) absorption core, the cores of the components of Fe II absorptions, respectively. The last column lists velocity from the Fe II lines with a low lower-level excitation potential. In the spectra for four dates, these lines are split into two components, hence two mean velocity values are given: based on the long-wave components, and below these values the mean velocity for the short-wave components is indicated. Let us stress that for all the times when it is present in the spectrum the splitting reaches large values: – km s. In the 2000 spectrum, due to its reduced spectral resolution the absorption component measurements are unreliable. For this reason, the data for this spectrum are not given in Fig. 4, which presents the dependences of radial velocity from the central residual intensity of the corresponding line, . One or two signs (in cases of the split absorptions) correspond to each line. The circles mark the high-excitation Fe II absorptions, the squares – the low-excitation Fe II absorptions, the rings – He I, S II, the horizontal lines – mainly the Fe II 6318, 6384, 6385 Å, emissions, and in the spectra of 2017 – the Fe II 6493, 7496, 7513 Å emissions as well. These Fe II emissions are formed in extended envelopes and, being stationary for some super- and hypergigants, are used to estimate their systemic velocities [Chentsov2004 ()]. According to our data, in the investigated star V534 Lyr the velocities determined based on the emissions at different dates of 2010 differ by approximately 5 km s, and different lines are shifted by various values, though all in one direction. In the 2013 spectrum, the position of these emissions is significantly different from the previous ones. At the same time, in each of the dates, the position of the emissions is close to those for the unsplit absorptions, and hence the velocities in column (5) are close to the velocities in column (6), measured from the upper parts of absorptions.
The half-width of the Mg II 4481 Å line in the V534 Lyr spectrum is Å. Using the dependences of the half-widths of this line on the rotation velocity of the star from [Kinman2000 ()], we get a low rotation velocity – km s.
Table 2 and Fig. 4 show temporal variations of the positions of all lines and profile shapes of some of them. The latter phenomenon is mainly associated with the upper layers of the atmosphere. This can be seen, for example, from the fact that the low-excitation Fe II absorptions in the spectra obtained on 6 April, 2010, June 1, 2010 and October 12, 2013 are split, and the high-excitation absorptions formed deeper remain single in these dates. Figure 5 demonstrates this for the Fe II 5363 and 5506 Å lines (the low-level excitation potentials 3.2 eV and 10.2 eV respectively), while Fig. 6 – for the Fe II 5169 Å line, the strongest of the of the low-excitation absorptions in the visible part of the spectrum (its low-level excitation potential is 2.9 eV), and Si II 6371 Å (the lower-level potential is 8.1 eV). In the bottom plot of Fig. 5 (July 30, 2010 and September 24, 2010) the profile of Fe II 5363 Å is only slightly asymmetric, and in the top plot (April 6, 2010) it is split into two components, close in depth and spaced from each other by 48 km s. The profile of Fe II 5169 Å (Fig. 6) from April 6, 2010 also consists of two components, though possessing different depths, separated by 56 km s. Figure 6 shows an intermediate case of profile transformation. The absorption of Si II 6371 Å (the lower-level excitation potential of 8.1 eV) in the results of April 6, 2010 observations is clearly not split, but it is noticeably broadened as compared with the data obtained on September 24, 2010 (in the spectrum of July 30, 2010 it is absent).
|Fe II 5506|
|Fe II 5363|
|Fe II 5169|
|Si II 6371|
Table 3 contains the parameters of the line profiles, presented in Figs. 5 and 6, their depths , widths in the wings (at ) and equivalent widths . In the spectrum of April 6, 2010 the absorptions are shallower and wider in the wings than in the data from July 30, 2010 and September 24, 2010, while these differences are relatively small in the high-excitation Fe II lines, as well as in the He I and S II lines, and are noticeably larger in the low-excitation Fe II lines. The depths of Fe II 5506 Å, determined from the spectra on the dates indicated are practically identical, while for Fe II 5363 Å and Fe II 5169 Å they differ by 2.5 and 2.6 times, respectively. The absorption growth of the Fe II 5363 Å line in the observations of April 6, 2010 was by September 24, 2010 compensated by the narrowing the way that its equivalent width has returned to the previous value. We believe that the broadening and splitting of the low-excitation Fe II lines are related to the features of the velocity field in the upper layers of the atmosphere, rather than to the appearance of the emission components in them. In the region of the V534 Lyr spectrum available to us, the latter is confidently observed only in the hydrogen lines. At certain moments, as noted above, the spectrum revealed weak Fe II 6318, 6384, 6385 Å emissions. In the 2017 spectra, the Fe II emissions in the longer-wavelength region of 7388, 7496, 7513 Å were observed.
The difference in velocities is small: the variations from date to date are from 11 km s for He I and S II to 24 km s for the high-excitation Fe II absorptions. From the data given in the tables and figures with line profiles, it can be concluded that the radial velocity gradient in the atmosphere of V534 Lyr was minimal, and it was the most stable on July 30, 2010, September 24, 2010 and on June 8, 2017: the differential shifts of lines on these dates are close to the measurement errors, the widths and anomalies of their profile shapes are minimal. Assuming that the average radial velocity on these dates were close to the velocity of the center of mass of the star, we accept km s ( km s) as the first approximation. This estimate is also close to the average velocity values in columns (5)–(7) of Table 2. Close values of the means and amplitudes are given in [Bonsack1956 ()] and [Giridhar2005 ()], but one must bear in mind that in these papers the velocities were obtained by averaging over all the measured lines, not analyzing their features.
3.3 Luminosity and distance of the star
The first estimate of the two-dimensional spectral class of V534 Lyr was made by V. Morgan: B9 Ib [Morgan1950 ()]. Bonsack and Greenstein [Bonsack1956 ()] changed its spectral class to A0 Iabe, based on the fact that the Balmer series in the spectrum extends to H, and the H and H lines have the emission components, while the amplitude of radial velocity variation is higher than the typical for the Ib supergiants. The latter argument reinforces the significant and variable differential shifts of the lines, in particular, the Balmer progress observed on October 12, 2013 (the growth of from H to H from km s from km s), and noted by the authors of [Giridhar2005 (); Bonsack1956 ()] and observed by us variability of intensities and positions of the emission features of the H I line profiles. As can be seen in Fig. 1, in the 2000 and 2010 spectra the H emission is two-peak, and the intensity of the red peak is always above the blue one. However, in June 2017, H has a P Cyg-type profile, which in a couple of months acquires the character of a two-peak emission. The variability of the H profile is similar: as illustrated in Fig. 2, its emission component passes from the red wing to the blue one. According to the results of observations from June 8, 2017, the H profile is of a P Cyg-type with a powerful and shifted into the short-wave region absorption, which gets split in the spectrum obtained a few days later, on June 13, 2017.
The spectral class A0 Iab corresponds to the distance to V534 Lyr of kpc [Bonsack1956 ()]. The interstellar Na I (1) and Ca II (1) line profiles in the stellar spectrum also imply its great remoteness. The presence in them of the components with km s (Fig. 7), taking into account the data of Brand and Blitz [BrandBlitz ()] points to kpc. A great remoteness of the star could be confirmed by its parallax mas, measured by the GAIA, however, in the DR1 catalog it is burdened with a large error of mas. A noticeable own motion in the Leeuven () catalog at kpc corresponds to the velocity 103 km s, directed to the plane of the Galaxy. In combination with km s this yields the spatial velocity of about 140 km s.
The presence of interstellar components in the Na I (1) and Ca II (1) lines forced us to undertake a search for possible DIBs. Luna et al. [Luna2008 ()] in their study of interstellar and circumstellar absorption for the sample of post-AGB stars, found two such features (5780 and 5797 Å) in the spectrum of V534 Lyr. However, their measurements of the DIB positions are widely scattered. In our spectra, there are no features that could be reliably identified with the known DIBs. The only feature that can be identified with the strongest of the known bands, 5780 Å is very weak in the spectrum of V534 Lyr – its depth is only around 0.015.
The obtained distance estimate of kpc for the high-latitude star with the apparent stellar magnitude of [Fernie1983 ()] leads us to the value of its absolute magnitude: , which, according to [Straizys ()], corresponds to the spectral class of V534 Lyr. Vickers et al. [Vickers15 ()], modeling the energy distribution in the spectrum stars, obtained the distance of kpc. This reduced distance estimate leads to a decrease in the absolute magnitude, , and worsens the agreement of spectral class and luminosity.
The luminosity of the star can be determined from the equivalent width of the oxygen IR triplet O I 7773 Å in the spectrum of V534 Lyr: Å. Using the calibration from [Giridhar2003 (); IRC2 ()], we get a very high luminosity, . However, this estimate is burdened by considerable errors, since both calibrations were obtained by the authors of [Giridhar2003 (); IRC2 ()] for the population I supergiants. The star we are investigating, V534 Lyr, obviously does not refer to massive supergiants. In addition, its luminosity estimate is also distorted by the specific chemical composition of the atmosphere, primarily, by an anomalous oxygen abundance.
3.4 Determination of model atmosphere parameters and chemical abundances
To determine the main parameters of the model: effective temperature and surface gravity , we used our standard method [ThetaLeo ()], which was successfully used in the studies of spectra of various types of stars. The parameters and were determined based on the requirement of the ionization balance, i.e. the equality of iron abundance, calculated from the Fe I and Fe II lines. Microturbulent velocity was also found using the standard method from the condition of absence of the dependence of iron abundances (Fe I, Fe II), determined from a set of lines, on their equivalent widths . To account for possible variations, the parameters and were found for each available time of observations. The obtained individual temperatures , given in column (4) of the Table 2, differ from each other only slightly. Their differences are within the error limits of determining this parameter; therefore, to calculate the abundances of chemical elements, we used the mean value of K. This estimate agrees with the previous estimate from [Giridhar2005 ()] perfectly well, which is an additional indication of the constancy of stellar temperature. Note that earlier in the paper by Arellano Ferro et al. [Giridhar2001 ()] the condition of solar helium abundance in the atmosphere was accepted for determining the temperature of V534 Lyr which led to a higher value K.
We controlled the identification of features in the spectrum of V1 534 Lyr involving the atlas [Deneb ()] with high-resolution spectra of A-supergiants. The content of chemical elements in the atmosphere of the studied star is calculated from the equivalent line widths in the LTE approximation. The calculation of the model atmosphere, the abundance of chemical elements and the calculation of synthetic spectra are performed using the WIDTH9 code based on the Kurucz model grid [kurucz ()], adapted for the OS Linux environment [synth ()]. The excitation potentials and oscillator strengths for all the lines, as well as the broadening constants are taken from the Vienna Atomic Line Database (VALD) [vald ()].
Taking into account the above-described peculiarities of the V534 Lyr spectrum, we calculated its chemical composition from the equivalent widths of the absorptions measured in the spectrum of June 8, 2017. In the choice of observational material, we were guided by the registered wavelength range and the absence of absorption splitting in the spectrum (Fig. 4 and Table 2). To increase the number of lines we partially involved the data on the equivalent widths of unsplit lines from the 2013 spectrum. The results of the calculation performed with the parameters K, , are presented in Table 4, where the number of lines used and the root-mean-square error of determining the element abundances are indicated. In the last column of the table the relative values of are given, for the calculations of which we used the information about the element abundances in the photosphere of the Sun according to Asplund et al. [sun ()].
The presence of an emission-absorption profile of the H line, in which the position of the emission components and their intensity ratio vary with time is a sign of passage of a shock wave in the stellar atmosphere. In addition, we see a low-amplitude variability of radial velocity from the lines with a high excitation potential, which are formed in deep layers of the stellar atmosphere. Therefore, we have registered the signs of pulsational instability of the star. In the spectra of V534 Lyr we found a previously unknown feature: splitting of profiles of selected metal absorptions at separate times of observations. This kind of metal line splitting indicates the pulsations of W Vir-type stars. A good example is the velocity field in the atmosphere of W Vir itself [Kovtyukh ()]. An additional argument confirming the presence of pulsations is given to us by the analysis of individual velocities based on the split lines in the spectrum of June 13, 2017, where the degree of splitting of the low-excitation absorptions varies depending on the depth of absorptions. This difference is illustrated by the data in the last line of Table 2 and the panel of dependences in Fig. 4, obtained on June 13, 2017.
The parameters of the model atmosphere and the chemical element abundances in the atmosphere of V534 Lyr we have determined differ little from the previously published [Giridhar2005 ()]. Here we will consider the features of the chemical composition only briefly. A detailed analysis of the chemical composition of this star is difficult for a number of reasons: a very limited set of obtained abundances of chemical elements; application of the standard model atmosphere in the case of a star whose atmosphere is unstable and probably subject to the influence of shock waves. A contribution of error due to the neglect of the effect of deviation from the local thermodynamic equilibrium is also possible. However, their impact on the metallicity estimate for a hot star is small [Cohen ()].
The main groups of elements whose relative abundances allow us to judge on the stage of evolution of the star and its belonging to this or that population of the Galaxy are as follows: the CNO-triad, the iron group (Cr, V, Fe, Ni), the -process light metals (Mg, Al, Si, P, S) and the -process heavy metals. The abundance of iron in the atmosphere of V534 Lyr is slightly lowered: , which, in combination with a high radial velocity km s points to the fact that the star belongs to the thick disk of the Galaxy. The abundance of metals of the iron group (Cr, V, Mn, Ni) is slightly different of the iron abundance: .
To determine the status of a far-evolved star, the principal in the element prevalence is the abundance of elements of the CNO group. We have reliably determined a high nitrogen abundance from two N I lines of low intensity. The nitrogen excess in the atmosphere of the supergiant can be the result of the first mixing, at which the removal of products of the CN-cycle during the burning of hydrogen in the core is in force. A large helium excess, is a result of a sequence of nuclear reactions in the stellar core and their subsequent mixing. The excess of oxygen in V534 Lyr is illustrated in the bottom panel of Fig. 3, where the fragments of the spectra of V534 Lyr and Cyg containing the O I 6155–6157 Å oxygen lines are mapped. We can also clearly see here a weakening of the Fe II ion lines in the spectrum of V534 Lyr in comparison with the spectrum of Cyg. The abundances of light metals (Mg, Al, Si, P, S), synthesized owing to the -process are strengthened in the atmosphere of V534 Lyr on the average by . The strengthening of light metals is consistent with the fact that the star pertains to the thick disk of the Galaxy, as evidenced by the results of the extensive sample of thick-disk stars with a close metallicity [Reddy2006 ()].
With respect to the observed set of properties, a hot star HD 105262, which is located in the Galaxy at a high latitude of and has the B9–A0 spectral class can be regarded as an analog of V534 Lyr. In the SIMBAD database this star is cataloged as a post-AGB supergiant. Earlier this star, having no IR excess, was considered as a representative of an earlier stage of evolution, namely, of the horizontal branch (hereafter – HB). Klochkova and Panchuk [KP1988 ()], using the photographic spectra determined the parameters of K, the surface gravity , reduced metallicity and a detailed chemical composition of the atmosphere. In their recent work, Giridhar et al. [Giridhar2010 ()] obtained the same parameters of the model, though the metallicity is much lower: . This difference can be explained by a higher spectral resolution of the spectra in [Giridhar2010 ()]. For the comparison we make between V534 Lyr and HD 105262, a large excess of nitrogen, found in the atmosphere of HD 105262 [Giridhar2010 ()] is important. Both stars have a low rotation velocity, km s. For HD 105262, rotation velocity is measured from the high-resolution spectra in [Martin ()]. Belonging to the HB stage can explain the absence in V534 Lyr and HD 105262 of the circumstellar dust and infrared excess caused by it.
The spectral features, similar to the peculiarities of the spectrum V534 Lyr are observed in the A-supergiant BD+481220. The optical high-resolution spectra based on which a significant variability of the H I and metal line profiles were obtained over several observational seasons at the 6-m BTA telescope combined with an echelle spectrograph [IBVS ()]. The analysis of this spectral material by the atmospheric model method with the parameters K, , has shown that the metallicity BD+481220 is close to solar: [IRAS05040 ()], and the chemical composition of its atmosphere differs little from that of V534 Lyr. As in the case of V534 Lyr, a large helium excess was detected, and an equally significant excess of oxygen dex. At that, the carbon excess is small, , and the ratio . The abundances of light metals are changed: at . But the most important – in the atmosphere of BD+481220 a large excess of lithium is detected, what indicates the transfer of this element, synthesized at the AGB stage, into the atmosphere. In addition, it was concluded on the probable efficiency of the mechanism of selective separation of chemical elements onto the dust particles of the envelope. A full set of the available data (luminosity , velocity km s, reduced metallicity and the features of the optical spectrum and chemical composition) confirms the status of the O-rich post-AGB star with the initial mass of . Thus, the star BD+481220, having a number of properties close to V534 Lyr is a supergiant at a more advanced stage, after AGB.
We therefore see that the observed properties of V534 Lyr, with the exception of its high luminosity and location outside the plane of the Galaxy do not give any grounds for classifying it as a post-AGB star. The main fact is that V534 Lyr does not have an IR excess or chemical composition anomalies expected for the post-AGB supergiant. Let us consider other possible variants of its status. The lack of a dust envelope at a sufficiently high luminosity allows to consider it as an object at the stage after the red giant branch, located at an evolutionary stage above the HB. The source of energy release in the stars at this stage, past the helium flash in the core is the burning of helium in the core and hydrogen in the layer. A large mass loss is not expected in these supergiants. To confirm this status of the star one can rely on a convenient diagram – (Fig. 1 in [Martin ()]). The set of fundamental parameters of V534 Lyr corresponds to the HB stars in the diagram. A low rotation velocity – km s is also consistent with the supposed status near the BHB: it is known that HB stars rotate slowly [Peterson (); Cohen ()]. The results we obtained here for V534 Lyr make up another confirmation of the conclusion [Envelopes (); Evolution ()] on the heterogeneity of the sample of PPN candidates.
A combination of the observed features of V534 Lyr: the presence of pulsations in deep layers of the atmosphere, reduced metallicity, the type and variability of the emission-absorption H profile allows us to assume that the star belongs to pulsating population II stars that are located in the instability band above the HB and evolve to the AGB. Depending on the mass and hence on the period of pulsations, it can be a BL Her or W Vir-type star. A direct indication of the status of the pulsating star is provided by the optical spectrum features: two-peak emissions and time-variable H and H line profiles (Figs. 1 and 2), the presence at some times of observations of the absorption splitting (Fig. 5) and the presence of a velocity gradient in the atmosphere, registered on June 13, 2017, which is clearly visible in Fig. 4.
It is useful to pay attention here to the study of the chemical composition of a sample of 19 variable population II stars [Maas ()], based on the high-resolution spectroscopy. To calculate the abundances of chemical elements, the authors used only the spectra in which there were no signs of splitting (or asymmetry) of absorptions or emissions in the Balmer line series. Maas et al. [Maas ()] made a conclusion about the fundamental difference of the chemical composition in BL Her or W Vir-type stars: in the atmospheres of the BL Her-type stars, the Na abundance is higher than that in the W Vir-type stars. This conclusion gives us grounds to more likely classify the investigated star V534 Lyr to the Virginids: we failed to find any of the Na I lines in the spectra of V534 Lyr, expected in the case of an excess of this element in the atmosphere.
Typical representatives of the W Vir-type stars are the globular cluster members. An example is the Virginid V1 (K 307), a member of the globular cluster M 12. The authors of [M12 ()] determined its main parameters ( K, ) and metallicity , which agrees with the metallicity of other members of M 12. The atmosphere of this Virginid revealed altered abundances of CNO-elements, while nitrogen is in excess: dex. Relative abundances of Na and the -process elements: Mg, Al, Si, Ca, and Ti are also strengthened in a various degree. The main feature of the spectrum of K 307 is the presence of complex absorption-emission profiles of the H and H lines (see Fig. 2 in [M12 ()]), similar to the profiles in the spectrum of V534 Lyr.
The basic kinematic and chemical parameters for an extensive sample of HB stars in the thick disk are considered by Kinman et al. [Kinman2009 ()]. Our assumption that V534 Lyr can be a thick disk star in the stage above HB (post-HB) is confirmed by a large spatial velocity of the studied star, the absence of a dusty circumstellar envelope, and a specific chemical composition of the atmosphere. A factor complicating the fixation of the stellar status is that stars of several evolutionary stages can be located above the HB: it can be young enough stars, evolving to RGB, as well as the stars more advanced in the stellar evolution after or yet before the AGB.
The answer to the question critical for more specific conclusions is on the possible binarity of V534 Lyr. The variations with time are probably caused by pulsations, however, a spectral SB 1-type binarity is not yet excluded. The found velocity variability based on the Fe II emission (2) can be considered an indication on the possible binarity of the star. Note here the results of Şahin and Jeffery [S2007 ()], who, having analyzed the variability of photometric data for V534 Lyr did not find an explicit period and assumed that the variations can be short-period. There are currently no reliable photometric and spectral data required for the study and determination of the variability parameters of brightness and spectrum. Therefore, it is extremely necessary to conduct a long-term and detailed monitoring of V534 Lyr for the final conclusions about the nature of its variability.
Based on the observations at the 6-m telescope with the NES echelle spectrograph ( = 60 000) we investigated the features of the optical spectrum of V534 Lyr, a high-latitude star with an uncertain status. For all the observation dates, heliocentric radial velocities correspond to the positions of all the components of metal absorption, as well as the Na I D and H lines were measured. The analysis of the velocity field from the lines of different nature revealed a small-amplitude variability from the lines with a high excitation potential, which are formed in deep layers of stellar atmosphere, and allowed us to estimate the systemic velocity km s ( km s). The estimated distance kpc for a high-latitude star leads to the absolute magnitude estimate of , which agrees with the spectral classification of V534 Lyr.
We discovered an as yet unknown for this star spectral phenomenon: splitting of profiles of the selected metal absorptions at separate times of observation. For all the moments when splitting is present in the spectrum, it reaches large values: – km s.
The spectral class of the star is close to A0 Ib, while its effective temperature is K. Metallicity, reliably determined by metal lines of the iron group (Cr, V, Mn, Ni) slightly differs from the abundance of this element: . The nitrogen and helium excess indicates an advanced evolutionary stage of the star. A low iron-group metal abundance in combination with a high radial velocity indicates that the star belongs to the thick disk of the Galaxy.
The whole of observed peculiarities of V534 Lyr: the presence of probable pulsations in deep layers of the atmosphere, splitting of metal absorption profiles with the low lower-level excitation potential registered at separate times of observations, reduced metallicity, type and variability of the emission-absorption H and H profile allows us to place the star amid the pulsating population II stars, which are located in the instability band near the HB. In general, we can conclude that there is a complete discrepancy between the properties of V534 Lyr we have registered and the post-AGB stage, where the star was referred to in the previously published papers.
The work was supported by the Russian Science Foundation (grant No. 14–50–00043, direction “Magnetometry of Stars” ). E.G.S. thanks the Russian Foundation for Basic Research for the financial support (project 16–02–00587 a). This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France, and NASA’s Astrophysics Data System.
- (1) 1. F. Herwig, ARA&A 43, 435 (2005).
- (2) 2. V.G. Klochkova, Astrophys. Bull. 69, 279 (2014).
- (3) 3. V.G. Klochkova and V.E. Panchuk, Astron. Rep. 60, 344 (2016).
- (4) 4. V.G. Klochkova, ASP Conf. Ser. 510, 121 (2017).
- (5) 5. V.G. Klochkova, R. Szczerba, V.E. Panchuk, and K. Volk, Astron and Astrophys.345, 905 (1999).
- (6) 6. V.G. Klochkova, R. Szczerba, and V.E. Panchuk, Astron. Let. 26, 88 (2000).
- (7) 7. V.G. Klochkova, R. Szczerba, and V.E. Panchuk, Astron. Let. 26, 439 (2000).
- (8) 8. V.G. Klochkova and T. Kipper, Baltic Astronomy 15, 395 (2006).
- (9) 9. V.G. Klochkova, E.L. Chentsov, and V.E. Panchuk, Astrophys. Bull. 63, 112 (2008).
- (10) 10. H. van Winckel and M. Reyniers, Astron and Astrophys. 354, 135 (2000).
- (11) 11. V.G. Klochkova, V.E. Panchuk, and N.S. Tavolzhanskaya, Astron. Let. 42, 815 (2016).
- (12) 12. V.G. Klochkova, E.L. Chentsov, V.E. Panchuk, et al., Astron. Rep. 61, 962 (2017)
- (13) 13. V.G. Klochkova, E.L. Chentsov, and V.E. Panchuk, MNRAS 292, 19 (1997).
- (14) 14. R.D. Oudmaijer, Astron and Astrophys. Suppl. 129, 541 (1998).
- (15) 15, V.G. Klochkova, E.L. Chentsov, A.S. Miroshnichenko, et al., MNRAS 459, 4183 (2016).
- (16) 16. T. Şahin, D.L. Lambert, V.G. Klochkova, and V.E. Panchuk, MNRAS 461, 4071 (2016).
- (17) 17. N.R. Trams, L.B.F.M. Waters, H.J.G.L.M. Lamers, et al., Astron and Astrophys. Suppl. 87, 361 (1991).
- (18) 18. W.K. Bonsack and J.L. Greenstein, PASP 68, 249 (1956).
- (19) 19. A. Arellano Ferro, S. Giridhar, and P. Mathias, Astron and Astrophys. 368, 250 (2001).
- (20) 20. F. Rufener, Astron and Astrophys. Suppl. 26, 275 (1976).
- (21) 21. J.D. Fernie, ApJS 52, 7 (1983).
- (22) 22. B. Hauck and M. Mermilliod, Astron and Astrophys. Suppl. 129, 431 (1998).
- (23) 23. W.E.C.J. van der Veen, L.B.F.M. Waters, N.R. Trams, and H.E. Matthews, Astron and Astrophys. 285, 551 (1994).
- (24) 24. S. Giridhar and A. Arellano Ferro, Astron and Astrophys. 443, 297 (2005).
- (25) 25. V.E. Panchuk, V.G. Klochkova, and M.V. Yushkin, Astron. Rep. 61, 820 (2017).
- (26) 26. V. E. Panchuk, I. D. Najdenov, V. G. Klochkova, et al., Bull. Spec. Astrophys.Observ. 44, 127 (1997).
- (27) 27. V. Panchuk, M. Yushkin, V. Klochkova, and M. Sachkov, Proc. SPIE, 9908, 99086Y (2016).
- (28) 28. M.V. Yushkin and V.G. Klochkova, Preprint No. 206, Spec. Astrophys. Obs. of RAS, Nizhnij Arkhyz, (2005).
- (29) 29. G.A. Galazutdinov, Preprint No. 92, Spec. Astrophys. Obs. of RAS, Nizhnij Arkhyz, (1996).
- (30) 30. V.G. Klochkova, V.E. Panchuk, M.V. Yushkin, and D.S. Nasonov, Astrophys. Bull. 63, 386 (2008).
- (31) 31. E.L. Chentsov and E.S. Gorda, Astron. Let. 30, 461 (2004).
- (32) 32. T. Kinman, F. Castelli, C. Cacciari, et al., Astron and Astrophys. 364, 102 (2000).
- (33) 33. W.W. Morgan and N.G. Roman, ApJ 112, 362 (1950).
- (34) 34. J. Brand and L. Blitz, Astron and Astrophys. 275, 67 (1993).
- (35) 35. F. van Leeuwen, Astron and Astrophys. 474, 653 (2007).
- (36) 36. R. Luna, N.L.J. Cox, M.A. Satorre, et al., Astron and Astrophys. 480, 133 (2008).
- (37) 37. V. Straizys and G. Kuriliene, Astrophys Space Sci, 80, 353 (1981).
- (38) 38. S.B. Vickers, D.J. Frew, Q.A. Parker, and I.S. Bojic̃ić, MNRAS 447, 1673 (2015).
- (39) 39. A. Arellano Ferro, S. Giridhar, and E. Rojo Arellano, Rev. Mexicana Astronom. Astrofis. 39, 3 (2003).
- (40) 40. V.G. Klochkova, M.V. Yushkin, E.L. Chentsov, and V.E. Panchuk, Astron. Rep. 46, 139 (2002).
- (41) 41. V.G. Klochkova, V.E. Panchuk, and V.V. Tsymbal, Bull. Spec. Astrophys.Observ. 19, 19 (1985).
- (42) 42. V.G. Klochkova, E.G. Sendzikas, and E.L. Chentsov, Astrophys. Bull. 70, 99 (2015).
- (43) 43. R. Kurucz, CD-ROM No. 13 (Smithsonian Astrophys. Obs., Cambridge, 1993).
- (44) 44. V. Tsymbal, ASP Conf. Ser. 108, 198 (1996).
- (45) 45. N.E. Piskunov, F. Kupka, T.A. Ryabchikova, et al., Astron and Astrophys. Suppl. 112, 525 (1995).
- (46) 46. M. Asplund, N. Grevesse, A.J. Sauval, and P. Scott, ARA&A 47, 481 (2009).
- (47) 47. V.V. Kovtyukh, G. Wallerstein, S.M. Andrievsky, et al., Astron and Astrophys. 526, A116 (2011).
- (48) 48. J.G. Cohen and J.K. McCarthy, AJ 113, 1353 (1997).
- (49) 49. B.E. Reddy, D.L. Lambert, and C. Allende Prieto, MNRAS 367, 1329 (2006).
- (50) 50. V.G. Klochkova and V.E. Panchuk, Astrofiz. Issled. (Izvestiya SAO) 26, 23 (1988).
- (51) 51. S. Giridhar, R. Molina, A.A. Ferro, and G. Selvakumar, MNRAS 406, 290 (2010).
- (52) 52. J.C. Martin, AJ 128, 2474 (2004).
- (53) 53. V.G. Klochkova, E.L. Chentsov, V.E. Panchuk, and M.V. Yushkin, IBVS, No. 5584 (2004).
- (54) 54. V.G. Klochkova, E.L. Chentsov, N.S. Tavolzhanskaya, and V.E. Panchuk, Astron. Rep. 51, 642 (2007).
- (55) 55. R.C. Peterson, ApJ 275, 737 (1983).
- (56) 56. V.G. Klochkova, Astrophys. Bull. 67, 385 (2012).
- (57) 57. T. Maas, S. Giridhar, and D.L. Lambert, ApJ 666, 378 (2007).
- (58) 58. V.G. Klochkova, V.E. Panchuk, N.S. Tavolganskaya, and V. V.Kovtyukh, Astron. Let. 29, 748 (2003).
- (59) 59. T.D. Kinman, H.L. Morrison, and W.R. Brown, AJ 137, 3198 (2009).
- (60) 60. T. Şahin and C.S. Jeffery, Astron. Nachr. 328, 848 (2007).
Translated by A. Zyazeva