A&A 390, 621-626 (2002)
DOI: 10.1051/0004-6361:20020732

Models of infrared spectra of Sakurai's Object (V4334 Sgr) in 1997[*]

Y. V. Pavlenko1,2 - T. R. Geballe 3


1 - Main Astronomical Observatory of Academy of Sciences of Ukraine, Golosiiv woods, Kyiv-127, 03680 Ukraine
2 - Isaac Newton Institute of Chile, Kiev Branch
3 - Gemini Observatory, 670 N. A'ohoku Place, Hilo, HI 96720 USA

Received 5 March 2002 / Accepted 17 May 2002

Abstract
Theoretical spectral energy distributions computed for a grid of hydrogen-deficient and carbon-rich model atmospheres have been compared with the observed infrared (1-2.5 $\mu $m) spectra of V4334 Sgr (Sakurai's Object) on 1997 April 21 and July 13. The comparison yields an effective temperature of $T_{\rm eff}$ = 5500 $\pm$ 200 K for the April date and $T_{\rm eff}$ = 5250 $\pm$ 200 K for July. The observed spectra are well fitted by Asplund et al. (1999) abundances, except that the carbon abundance is higher by 0.3 dex. Hot dust produces significant excess continuum at the long wavelength ends of the 1997 spectra.

Key words: stars: individual: V4334 Sgr (Sakurai's Object) - stars: AGB and post-AGB - stars: atmospheres -
stars: fundamental parameters


1 Introduction

V4334 Sgr (Sakurai's Object), the "novalike object in Sagittarius'' discovered by Y. Sakurai on February 20, 1996 (Nakano et al. 1996) is a very rare example of extremely fast evolution of a star during a very late final helium-burning event (Duerbeck & Benetti 1996). During the first few months after discovery, Sakurai's Object increased in visual brightness to ${V} \sim 12^{\rm m}$. In 1997 it increased further to ${V} \sim 11^{\rm m}$. In March 1997 the first evidence of dust formation was seen (Kimeswenger et al. 1997; Kamath & Ashok 1999; Kerber et al. 2000). In early 1998 the optical brightness of Sakurai's Object decreased (dimming first reported by Liller et al. 1998), but then recovered. However, during the second half of 1998 an avalanche-like growth of the dusty envelope occurred, causing a rapid decrease in optical brightness and the complete visual disappearance of the star in 1999. At present essentially only thermal emission by dust can be observed (Geballe et al. 2002). Our view of the born again star has been completely obscurred by the dust it has produced.

Abundance analyses by Asplund et al. (1997, 1999) and Kipper & Klochkova (1997) have found peculiarities similar to those of R CrB-like stars. Asplund et al. (1999) estimate that the logarithmic abundances of hydrogen, helium and carbon in atmosphere of Sakurai's Object in October 1996 were -2.42, -0.02 and -1.62, respectively[*], with hydrogen only the third most abundant element by number. All of the above studies are based on optical spectra obtained in 1996. At that time the spectrum of Sakurai's Object resembled that of an F-supergiant; molecular bands were absent or very weak. Cooling of the photosphere of Sakurai's Object resulted in its optical spectrum during 1997 and 1998 resembling those of C-giants with very strong bands of CN and C2 (Pavlenko et al. 2000). Modeling of some of these optical spectra have allowed estimates of the changes in $T_{\rm eff}$ and EB-V to be made during this period of rapid evolution of the optical spectrum (Pavlenko et al. 2000; Pavlenko & Duerbeck 2001).

Modeling of near infrared (1-2.5 $\mu $m) spectra of Sakurai's Object is of interest for several reasons. In addition to providing comparisons with results obtained from the optical spectrum and tests of the reliability of molecular and atomic data, it allows accurate determination of the effective temperature and sensitive tests for emission by hot dust. Use of the 1-2.5 $\mu $m region for modeling is especially important after 1996, when the bulk of the photospheric flux shifted from the optical into this waveband.

In this paper we present and compare model 1-2.5 $\mu $m spectra with those of Sakurai's Object obtained during 1997, on UT April 21 and July 13 at the United Kingdom Infrared Telescope (UKIRT). The observed spectra together with observational details were presented by Eyres et al. (1998) and the July spectrum is also shown in Geballe et al. (2002). The resolutions of these spectra as presented here are 1.4 nm (0.0014 $\mu $m) at 1.02-1.35 $\mu $m and 2.8 nm (0.0028 $\mu $m) at 1.42-2.52 $\mu $m. Narrow spectral features in the 1.82-1.95 $\mu $m portions of these spectra are due to incomplete removal of strong telluric lines.

2 Modeling procedure

Grids of plane-parallel model atmospheres in LTE, with no energy divergence were computed by the SAM12 program (Pavlenko 2002). This program is a modification of ATLAS12 (Kurucz 1999). Opacities due to C I bound-free absorption, of importance in atmospheres of hydrogen-deficient, carbon-rich stars over a wide (0.1-8 $\mu $m) wavelength region (see Pavlenko 1999, 2002; Asplund et al. 2000), were computed using the OPACITY PROJECT cross sections database (Seaton et al. 1992). The opacity of C- also was taken into account Myerscough & McDowell 1996; see Pavlenko 1999 for more details).

An opacity sampling approach (Sneden et al. 1976) was used to account for atomic and molecular line absorption. The source of the atomic line information was the VALD database (Kupka et al. 1999). Lists of diatomic molecular lines of 12CN, 13CN, 12C2, 13C2, 12C13C, 12CO, 13CO, SiH, and MgH were taken from Kurucz (1993). We adopted Voigt profiles for every absorption line; damping constants were computed following Unsold (1955). Microturbulent velocities of 3-6  km s-1 were adopted.

Two grids of model atmospheres with different abundances were computed. One grid used the chemical composition of Sakurai's Object determined by Asplund et al. (1999) for October 1996. The other is the same except that the abundance of carbon is increased by 0.6 dex. This "carbon-rich'' case is of interest because the carbon abundance of Asplund et al. (1997) is larger by 0.6 dex than that which was obtained by them later from the analysis of high resolution spectra (see Asplund et al. 1997, 1999, 2000. This "carbon problem'' appears to arise more from the analysis of C I lines than from C II or C2 lines (Asplund, private communication). For both grids the isotopic ratio 12C/13C = 5 was adopted (Asplund et al. 1997).

To determine molecular densities, a system of equations of chemical equilibrium was solved for a mixture of $\sim$70 atoms, ions and molecules, including the most abundant diatomic molecules containing carbon. We used the approach developed by Kurucz in ATLAS12 (Kurucz 1999) in which ratios of densities of atoms and nx, ...,nz and molecules nxy...z follow the equations:

 
$\displaystyle n_x*...*n_z/n_{x...z}=
\exp(-D_{0}/T_{\rm ev}+b
-c* (T+d*(T-e*(T+f*T)))+ 3/2*(m-k-1)*ln T),$     (1)

where D0 and $T_{\rm ev}$ are the dissociation potential and the temperature in eV (for neutral diatomic molecules m=1, k=0). The adopted molecular constants for the three most important molecules are given in Table 1.

 
Table 1: Dissociation equilibrium constants of CO, C2 and CN.
Molecule D0 b c*1.E2 d*1.E6 e*1.E10 f*1.E15
C2 6.116 48.75 .2192 .4149 .4121 1.550
CN 7.700 47.45 .1332 .1989 .1778 .6323
CO 11.105 49.45 .1651 .3103 .3100 1.168


Synthetic spectra[*] were computed using the WITA6 program (Pavlenko 1999) for the same grid of opacities, abundances, isotopic ratios and microturbulent velocites that were adopted for the model atmospheres. WITA6 computes spectra and spectral energy distributions (SEDs) taking into account "line by line'' absorption by atomic and molecular transitions. For Sakurai's Object spectra computed with wavelength step 0.05 nm were convolved with gaussian profiles with full widths at half maximum of 1.4 and 2.8 nm over the appropriate wavelength intervals. Computed and observed spectra were normalized at 1.7 $\mu $m for comparison; see Pavlenko et al. (2000) for more details.

3 Results

3.1 Principal spectral features

In Fig. 1 the principal features formed by the molecular species C2, CO, and CN are displayed as separate spectra. Atomic features are also shown, as these are also present in Sakurai's Object (see Eyres et al. 1998; Geballe et al. 2002). As in the optical spectrum (Pavlenko et al. 2000), absorption of only a few molecular species accounts for the main features in the IR spectrum. Only the most abundant isotopic species of each molecule is shown. Of the less abundant isotopic species, only bands of 13CO have been detected in the infrared (Eyres et al. 1998).

3.2 Dependences on T $\mathsfsl{_{eff}}$, log g, V $\mathsfsl{_t}$ and log N(H)

The model spectra of Sakurai's Object display a strong dependence on $T_{\rm eff}$  (Fig. 2).[*] In general, the dependence of the IR SED on $T_{\rm eff}$ is determined mainly by the variations of the molecular densities with temperature. The band strengths of CN, CO and C2 all increase as $T_{\rm eff}$ decreases. Changes in the continuum fluxes are much smaller. Similar effects are seen in model optical spectra (Pavlenko & Yakovina 2000). However, there the molecular bands are numerous, whereas in the infrared only the few strongest vibration-rotation bands of CN, C2, and CO are prominent.

  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{MS2440f1.eps}
\end{figure} Figure 1: Model spectra of species that produced the strongest absoption features in the 1.0-2.5 $\mu $m spectrum of Sakurai's Object spectra during 1997, computed for $T_{\rm eff}$ /log g = 5500/0.0 model atmosphere with Asplund et al. (1999) abundances for October 1996. The model spectrum due to atomic species alone, labelled VALD (see text) is also shown. Spectra are artificially shifted on the y-axis.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{MS2440f2.eps}
\end{figure} Figure 2: Dependence of the model IR spectrum on $T_{\rm eff}$. The model spectra use Asplund et al. (1999) abundances for October 1996. The observed spectrum of Sakurai's Object on July 13, 1997 is shifted on the y-axis.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{MS2440f3.eps}
\end{figure} Figure 3: Dependence of the model spectrum on log g.
Open with DEXTER

As can be seen in Fig. 3, the dependence of the spectrum on log g is generally rather weak. However, there are differences in the responses of different spectral regions. The strong molecular bands show rather weak dependence on log g, whereas the fluxes at 1.25-1.35, 1.60-1.75, 1.9-2.2 microns show more noticeable changes.

Previous abundance analyses of the spectra of Sakurai's Object and related R CrB stars indicate microturbulent velocities $V_{\rm t}$ in the range 5-8 km s-1 (cf. Asplund et al. 2000). The value of $V_{\rm t}$ affects the spectral distribution, as is shown in Fig. 4. The effect of $V_{\rm t}$ on the IR spectra of Sakurai's Object is larger at the heads of molecular bands than elsewhere, because the heads are formed by closely packed molecular lines whose overall absorption is sensitive to $V_{\rm t}$.

The main sources of line opacity in the model atmospheres approximating Sakurai's Object in 1997 are molecular (Pavlenko et al. 2000). Thus it is not surprising that the optical spectra which match Sakurai's Object respond weakly to changes in the hydrogen abundance. This is in contrast to the behavior of models corresponding to the star a year earlier (Asplund et al. 1997). Similarly, the model IR spectra of Sakurai's Object for $T_{\rm eff}$ = 5000-6000 K depend weakly on log N(H) (Fig. 5). The magnitude of the change in the spectrum when log N(H) is changed from -2.42 (the Asplund et al. 1999 value for October 1996) to -0.97 (i.e, a change of 1.5 dex) is comparable (in a qualitative sense) to lowering log g from 1 to 0 (Fig. 3).

3.3 Changes between 1997 April 21 and July 13

Fits to the spectra of Sakurai's Object on April 21 and July 13 are shown in Figs. 6 and 7. The long wavelength portion of the H band is of special interest for the "carbon problem'', because the strongest absorption bands of the C2 molecule, the Ballick-Ramsey bands, occur just longward of 1.768 $\mu $m. In the comparatively hot atmosphere of Sakurai's Object log N(C) > log N(O) (Asplund et al. 1997, 1999) and the abundance of C2 depends mainly on the elemental abundance of carbon. Therefore, these bands may provide the most accurate determination of log N(C). The fits imply that the carbon abundance is in the range log $N\rm (C) = -1.3 \pm 0.2$. The most likely value is 0.3 dex higher than that found by Asplund et al. (1999). The accuracy of the determination of log N(C) is limited mainly by the quality of the molecular line list.

The effective temperatures that best fit the 1.0-2.0 $\mu $m spectra in 1997 April and July are $5500\pm 200$ K and $5250 \pm 200$ K, respectively, indicating that the cooling evidenced by the dramatic spectral changes seen between 1996 and 1997 (e.g., Geballe et al. 2002) continued in 1997. Our estimated uncertainties in the above temperatures are rather large, despite the comparatively good fits to the observed spectra, because of questions concerning abundances, non-sphericity effects, and dynamical phenomena, and because of contamination of the spectra by dust emission (see below).

  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{MS2440f4.eps}
\end{figure} Figure 4: Dependence of the model spectrum on $V_{\rm t}$.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{MS2440f5.eps}
\end{figure} Figure 5: Dependence of the model spectrum of Sakurai's Object on log N(H).
Open with DEXTER

3.4 Hot dust

Emission by dust is evident in the 1997 spectra by the mismatch between the synthetic and observed spectra longward of 2.0 $\mu $m in Figs. 6 and 7. The difference between the observed and synthetic spectra is greater in the July spectrum, attesting to an increase in the amount of dust. The thermal emission from the dusty envelope overlaps the region of first overtone bands of 12CO and 13CO at $\lambda~>~$2.3 $\mu $m. Usually these bands are used for determination of carbon abundances and isotopic ratios (cf. Lazarro et al. 1991). The reduced equivalent widths of the CO bands in July 1997 cannot be reasonably attributed to a large decrease in the oxygen abundance, because (1) this is unlikely to have occurred in three months and because the continuum shortward of the CO bands also shows an excess. We note that in fitting spectra, the most frequent situation is that the computed spectra have excess flux due to the deficit of known or hypothetical opacities. To fit the observed spectra, opacities in the model would need to be decreased at $\lambda >$$\mu $m, an unrealistic possibility.

4 Discussion

Our analysis of the 1997 infrared spectra of Sakurai's Object strongly implies that dust was already present at that time. We note that Duerbeck (2002) did not find evidence for dust in the optical spectra from 1997. On the other hand, the fits by Pavlenko & Duerbeck (2001) to the observed SEDs at optical wavelengths indicate that EB-V had increased by 0.6 (from 0.7 to 1.3) from April 1997 to August 1998. However, the August 1998 data were best fit using $T_{\rm eff}$  = 5250 $\pm$ 200 K, the same value as for July 1997 in this paper. Between June 1997 and August 1998 there were some variations in the photospheric radiaton, probably caused by mass losses events, evolution of the dusty envelope, and dynamical processes in the photosphere - envelope system (see light curve of Sakurai's Object in Duerbeck 2002). Nevertheless, $T_{\rm eff}$ apparently remained nearly constant during this period.

In 1997, the year of maximum optical brightness of Sakurai's Object, the luminosity was still dominated by optical radiation. At that time "quasi-periodic fluctuations of increasing cycle length and amplitude were superimposed on the general brightness evolution'' (Duerbeck 2002). In general, the effective temperature during such fluctuations does not need to follow changes of luminosity. In fact, it can be anti-correlated, since an increased radius can more than compensate for a lower $T_{\rm eff}$ . On the other hand, a change of radius can change the thermodynamical properties in the radiating region (i.e. photosphere). That may explain the similarity of $T_{\rm eff}$ obtained in this paper for July 1997 and that found in August 1998 by Pavlenko & Duerbeck (2001). The decreased optical brightness of Sakurai's Object in 1998 was mainly caused by development of the dust envelope (Kimeswenger 1999).

One question arises - were the optical and 2 $\mu $m SED's being affected by the same dust in 1997-1998? The answer is probably yes. As mentioned earlier, the effective temperature remained constant during this time and thus cannot be the cause of the large change in EB-V. This suggests that the cause of the increase in EB-V was newly formed dust. The new dust would be expected to have been close to the star and thus quite hot. Indeed the full 1-5 $\mu $m spectrum from 1998 (e.g., Geballe et al. 2002) shows that the excess peaked close to 3 $\mu $m, indicating a mean dust temperature close to 1000 K at that time. The dust must have been hotter (and closer) in 1997; this is supported by the data from 1997 (Eyres et al. 1998; Geballe et al. 2002) which show that the continuum flux density decreased monotonically with wavelength in the observed wavelength range, 1-4 $\mu $m, at that time. Comparison of the 1997 and 1998 SEDs also show that much less dust was present in 1997. Thus we confirm that the first appearance of dust occured in 1997 and that the amount of dust increased through summer 1998, prior to its becoming totally dominant in the latter part of 1998 and since then.

  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{MS2440f6.eps}\par\includegraphics[width=8.8cm,clip]{MS2440f6a.eps}
\end{figure} Figure 6: Top: fits to observed spectrum of of Sakurai's Object on 1997 April 21. Bottom: details of the fits at 1.6-2.0 $\mu $m; much of the structure at 1.82-1.95 $\mu $m in the observed spectrum is due to incomplete removal of telluric absorption features. Synthetic spectra were computed for a microturbulent velocity of 6 km s-1.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=9cm,clip]{MS2440f7new.eps}\par\includegraphics[width=8.8cm,clip]{MS2440f7a.eps}
\end{figure} Figure 7: Top: fits to observed spectrum of of Sakurai's Object on 1997 July 13. Bottom: details of the fit at 1.6-2.0 $\mu $m. Synthetic spectra were computed for a microturbulent velocity of 6 km s-1.
Open with DEXTER

Acknowledgements

We thank the staff of the Joint Astronomy Centre for assistance in obtaining the spectra and the VALD database team for its helpful assistance. Partial financial support for YVP was provided by a Small Research Grant from the American Astronomical Society. TRG's research is supported by the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., on behalf of the international Gemini partnership of Argentina, Australia, Brazil, Canada, Chile, the United Kingdom, and the United States of America. We thank the referee, M. Asplund for several helpful suggestions.

References

 


Copyright ESO 2002