264 Acta Polytechnica CTU Proceedings 2(1): 264–268, 2015 264 doi: 10.14311/APP.2015.02.0264 On the 2011 Outburst of the Recurrent Nova T Pyxidis L. Izzo1, M. della Valle2, A. Ederoclite3, M. Henze4 1Dip. di Fisica, Sapienza University of Rome, P.le A. Moro 2, Rome; ICRANet, P.le della Repubblica 10, Pescara, Italy 2INAF - Osservatorio di Capodimonte, Salita Moiariello 16, 80131 Napoli, Italy 3Centro de Estudios de Fisica del Cosmos de Aragón, Plaza San Juan 1, Planta 2, Teruel, E44001, Spain 4European Space Astronomy Centre, P.O. Box 78, 28692 Villanueva de la Cañada, Madrid, Spain Corresponding author: luca.izzo@gmail.com Abstract We discuss the nebular phase emission during the 2011 outburst of the recurrent nova T Pyxidis and present preliminary results on the analysis of the line profiles. We also present some discussions about the binary system configurations and the X-ray emission, showing that the white dwarf mass should be larger than 0.8 M�. Keywords: spectroscopy - X-rays - individual: T Pyx. 1 Introduction The 2011 outburst of the recurrent nova T Pyx shed more light on some interesting characteristics of this peculiar nova. The optical evolution was not very dif- ferent from the last (1966) outburst, showing the 8 days ’shoulder’ before the final rise to the maximum lumi- nosity, which was reached ∼ 30 days from the initial outburst. The distance of T Pyx was determined from light echoes to be 4.8 ± 0.5 kpc [1], which is marginally consistent within the errors with the value of 3.5 ± 0.4, estimated from the application of the maximum mag- nitude versus rate of decline (MMRD) relation [2] to the light curve of the 1966 outburst [3]. The presence of a clumpy ring with an inclination of 30-40 degrees, embedded within a spherical shell, was also recently re- ported using Hubble Space Telescope observations dur- ing the nebular phase [1]. A variation of the orbital period after the 2011 outburst was detected, indicating an ejected mass in the outburst of 3 × 10−5 M� [4]. Also, the orbital period was found to have increased during the long quiescence period between the last two outbursts of T Pyx. Moreover, with respect to the 1966 epoch, when no advanced X-ray and radio detectors were available, we have now a good knowledge of the emission of T Pyx at these frequencies during its out- burst [5]. In the previous work [6], we have presented the early spectral evolution, up to two months after the initial outburst. We have shown the hybrid spectral transi- tion He/N - Fe II - He/N and initial ejecta velocities. In this work, we present the optical observations of the nebular phase, which was marked by the onset of the supersoft X-ray source (SSS) emission which happened during the seasonal gap [7]. The optical observations presented here were obtained with the high-resolution spectrograph SARG, mounted on the 3.5m telescope Nazionale Galileo (TNG) located at La Palma. Spectra with power resolution of R = 29000 were obtained but the complete sample, which includes also observations obtained at the Very Large Telescope, will be presented in Ederoclite et al. in preparation. Then we briefly discuss the X-ray emission and finally we provide an updated scheme of the binary system configuration in T Pyx, taking into account the recent results. 2 The Nebular Phase of T Pyx The nebular phase in novae is easily recognized by the absence of absorption lines and the presence of forbid- den emission lines, e.g. [OIII] and [NII], in the opti- cal spectra, similarly to planetary nebula cases. In- deed, these emission lines are due to photo-ionization of the ejecta by the central hot white dwarf (WD). At the same time, forbidden transitions originate mainly in rarified environments, so that any structure observed in corresponding line profiles is mainly due to the distri- bution of the material and not to self-absorption mech- anisms. The analysis of the [OIII] 4959-5007 and [NII] 5755 λλ profiles can consequently provide hints on the geometry of the ejecta. These line profiles are well described by a large flat top base (vexp,1 ≈ 1700 km/s) overlapped by a saddle- shaped and brighter emission, with vexp,2 ≈ 600 km/s, see Fig. 1 for the case of the [OIII] 5007 λ. The high resolution allows us to detect a castellated structure on 264 http://dx.doi.org/10.14311/APP.2015.02.0264 On the 2011 Outburst of the Recurrent Nova T Pyxidis the top of the profiles which is due to the presence of a clumpy ejecta structure. The saddle-shaped configura- tion suggests a non-spherical geometry of the ejecta, as generally discussed in [8] and further developed in sev- eral publications. In order to estimate the asphericity from forbidden lines we have used a method first pre- sented in [9]. An analytic expression for the line profiles I(l) from an anisotropic distribution of the nova ejecta was provided: I(l) = K ( c v0 )α+1 ×[ lα−22 − l α−2 1 α − 2 l2(1 − 3cos2i) + lα2 − lα1 α (1 − sin2i) ] , l ≤ l1;[ lα−22 − l α−2 α − 2 l2(1 − 3cos2i) + lα2 − lα α (1 − sin2i) ] , l ≥ l1; (1) where c is the speed of light, v0 is the velocity at the inner radius of the expanding shell (we have as- sumed the maximum velocity observed for the bright profile, v0 ≈ 600 km/s), α describes the variation of the shell’s velocity with the distance form the WD, and l = (λ − λ0)/λ0 = v/c is our variable. The quantity -0.006-0.004-0.002 0.000 0.002 0.004 0.006 0 5.´10-13 1.´10-12 1.5´10-12 displacement l from @OIIIDΛΛ 5007 fl ux He rg �c m 2 �s �A ng L Figure 1: The saddle-shaped profile in the [OIII] 5007 line and the best fit (red line) obtained with the method described in the text. l2 = v2/c corresponds to the zero intensity along the wing of the line profile. A grid-search algorithm was developed to find the best fit for the inclination an- gle i that the binary orbital plane forms with the line of sight. For all the three lines, we find a large angle i ≈ 60 ± 10 degree, see also Tab. 1, which contrasts some estimation of the i given before the 2011 outburst [14], and suggests a larger value for the binary inclina- tion angle, see also [1]. We note that this approach does not take into account the presence of a clumpy struc- ture. A more detailed approach is needed, and will be presented elsewhere. 3 The Configuration of the Binary System It is clear from the nebular spectra that the emission in T Pyx is well described by quasi-spherical ejecta, which give rise to the larger flat top profile, and an equato- rial ring of material, which is slower and gives rise to the saddle-shaped profiles. The origin of this equato- rial ring is unknown: it is possible that it formed dur- ing the interaction of the main ejecta with pre-existing gas, which is equatorially distributed around the binary system (an hypothesis presented in [13]); another pos- sibility can be that the early interaction of the ejecta with the secondary star could influence the subsequent expansion [10], but other mechanisms can also work. In all cases, it is necessary to shed light on the observed asphericity in T Pyx, and the first ingredient we need is the mass of the WD. Table 1: Fit results of the line profiles of [OIII] 4959, 5007 and [NII] 5755, using the approach developed in [9]. Ident. wavelength ( Å) i (degree) v1 (km/s) α K (erg/cm 2/s/Ang) 4959 59 ± 10 465 ± 6 -2.2 ± 0.3 (3.7 ± 0.1) × 10−15 5007 59 ± 10 477 ± 6 -2.3 ± 0.3 (1.09 ± 0.03) × 10−14 5755 60 ± 10 444 ± 6 -2.8 ± 0.3 (6.8 ± 0.5) × 10−16 265 L. Izzo et al. � � � � � � � � Figure 2: The location of the T Pyx parameters (green stars) with respect to the M 31 nova correlations presented in [20]: (1) the ton-toff correlation, for T Pyx measured from the time of maximum brightness; (2) the kT-toff relation, with kT given by the average value kT=55.7 eV; (3) the relation between ton and t2 and (4) the correlation between ton and the expansion velocity as measured from the Hα emission line. X-ray observations of the SSS phase can in princi- ple provide this physical quantity. The SSS emission is caused by stable hydrogen burning in the matter that was not ejected during the nova outburst. This burn- ing matter forms an envelope around the WD and can be observed once the ejected material becomes optically thin to soft X-rays. The duration of the SSS phase is in- versely related to the WD mass [11], since more massive WDs need to accrete less matter to trigger an outburst. Given its recurrent nature, T Pyx is expected to contain a massive WD. The SSS phase of T Pyx was first detected about 90 days after the outburst [7], when an increase in soft X-rays was observed by the Swift-XRT [12]. T Pyx re- mained in this state for ≈ 200 days, although around day 190 the luminosity started to decline. We do not en- ter here into a discussion on the details of this emission, which will be described in a forthcoming work. What is interesting in the context of this work is the relatively short duration of the SSS phase, which suggest a large WD mass. Observations obtained before 2011 reported a WD mass of 0.7 M� and a mass ratio of q = 0.2 [14], which implies a secondary star of 0.14 M�. The orbital period was observed with high accuracy to be Porb = 1.8295 hrs. With some assumption on the type of the compan- ion star, we can infer the mass of the secondary which overfills its Roche lobe: ∆M2 = ρ2 4 3 π(∆R2) 3, where ρ2 is the density and ∆R2 the dimensions of the overfilling envelope of the secondary. From the q value and the ab- sence of signatures of the secondary in optical spectra during the quiescent phase [3], we can assume that it is a low-mass M-type star. Its density can be well approxi- mated using the Eggleton formula ρ2 = 107Porb(hrs) −2 g cm−3 [15], while its radius can be derived from the Patterson et al. formula R2 = 0.62M 0.61 2 R� [16]. The Roche lobe size of the secondary is related to the q pa- rameter through L2 = 1.631q 1/3M 1/3 1 Porb(hrs) 2/3 1010 cm [17]. The result is that the envelope overfilling the sec- ondary Roche lobe is related to the WD mass. From the recent estimate of the mass ejecta after the 2011 266 On the 2011 Outburst of the Recurrent Nova T Pyxidis outburst (3 × 10−5M�, [4]), and assuming that every observed T Pyx outbursts ejected similar masses, we ob- tain a total mass ejected from T Pyx of about 2 × 10−4 M�. If this mass was accreted onto the WD from the secondary we can obtain a lower limit on the WD mass from the above formula. This results in MWD ≥ 0.8 M�. 4 Final Considerations Our last result requires some important remarks. First, it was proposed that the mass ejected in T Pyx is larger than the one accreted [3]. We do not consider this pos- sibility in this work, i.e., that all the accreted mass results to be the ejected one. We know that this is not exactly true, since the SSS phase is due to the burning of accreted hydrogen which is not ejected in the nova outburst. The mass of this burning envelope can be estimated assuming that the energy source for the SSS phase is the gravitational contraction [18], but this contribution can be neglected. Second, there is the possibility that a fraction of the matter overfilling the secondary’s Roche lobe would flow out from the ex- ternal Lagrangian point L2, leading to the formation of an asymmetric structure around the binary system. Such a mechanism would operate only if the primary Roche lobe has been filled by accreted gas, leading to the detection of transient heavy elements in absorption (THEA) during the early phases of the outburst [19], whose presence, in this case, is not certain. Moreover, in tight systems where the secondary is less massive than the WD (q < 1) also centrifugal forces can contribute to such an outflow. There is again observational evi- dence that T Pyx shares many properties common to other novae. As an example of this, we compare T Pyx to the unbiased sample of M31 novae [20], in particular with respect to the correlations between some optical and X-ray observables: the blackbody temperature of the SSS emitter, the turn-on and turn-off time of the SSS emission, the expansion velocity of the ejecta and the t2 determined from the R-band light curve [21], see Fig. 2. We see that T Pyx lies in the same region oc- cupied by classical novae observed in M31. All these results need to be figured in a unique picture for T Pyx and can shed more light on CV systems with a mass ratio of less than unity. Further results and conclusions will given in Izzo et al. (2015) in preparation. Acknowledgement We are very grateful to R. Williams, S. N. Shore, E. Mason, D. de Martino and R. Gilmozzi for precious discussions. M.H. acknowledges support from an ESA fellowship. References [1] Sokoloski, J. L., Crotts, A. P. S., Uthas, H., et al.: 2013, ApJ 770, 33. [2] Della Valle, M. & Livio, M.: 1998, ApJ 506, 818. [3] Selvelli, P., Cassatella, A., Gilmozzi, R., et al.: 2008, A&A 492, 787. [4] Patterson, J., Oksanen, A., Monard, B., et al.: 2013, arXiv/1303.0736. [5] Nelson, T., Chomiuk, L., Roy, N., et al.: 2012, ApJ, 785, 78. [6] Izzo, L., et al.: 2012, Mem. S.A.It. S. 83, 830. [7] Schwarz, G. J., Osborne, J. P., Page, K., et al.: 2011, ATEL 3754. [8] Payne-Gaposchkin, C. H.: 1957, The galactic novae. [9] Vainu Bappu, M. K. & Menzel, D. H.: 1954, ApJ 119, 508. doi:10.1086/145857 [10] Hutchings, J. B.: 1972, MNRAS 155, 357. doi:10.1093/mnras/155.3.357 [11] Starrfield, S., Truran, J. W., sparks, W. M., et al.: 1991, in Extreme Ultraviolet Astronomy, ed. R. Malina & S. Bowyer, NewYork, Pergamon, p. 168. [12] Gehrels, N., et al.: 2004, ApJ, 611, 1005. doi:10.1086/422091 [13] Williams, R.: 2012, AJ, 144, 98. [14] Uthas, H., Knigge, C. & Steeghs, D.: 2010, MNRAS, 409,237. doi:10.1111/j.1365-2966.2010.17046.x [15] Eggleton, P. P.: 1983, ApJ, 268, 368. doi:10.1086/160960 [16] Patterson, J., Kemp, J., Harvey, D. A., et al.: 2005, PASP, 117, 1204. doi:10.1086/447771 [17] Paczynski, B.: 1971, ARA&A, 9, 183. doi:10.1146/annurev.aa.09.090171.001151 [18] Krautter, J., Oegelman, H., Starrfield, S., et al.: 1996, ApJ, 456, 788. [19] Williams, R., Mason, E., Della Valle, M., et al.: 2008, ApJ, 685, 451. doi:10.1086/590056 [20] Henze, M., Pietsch, W., Haberl, F., et al.: 2014, A&A, 563, A2. [21] della Valle, M. & Livio, M.: 1995, ApJ, 452, 704. 267 http://dx.doi.org/10.1086/145857 http://dx.doi.org/10.1093/mnras/155.3.357 http://dx.doi.org/10.1086/422091 http://dx.doi.org/10.1111/j.1365-2966.2010.17046.x http://dx.doi.org/10.1086/160960 http://dx.doi.org/10.1086/447771 http://dx.doi.org/10.1146/annurev.aa.09.090171.001151 http://dx.doi.org/10.1086/590056 L. Izzo et al. DISCUSSION CHRISTIAN KNIGGE: The system parameter we estimated in Uthas et al. 2011 should be used with caution. They are based on a radial velocity study, and parameters resulting from such studies are known to be biased quite often. CHRISTIAN KNIGGE: You consider R2 and L2 differing by ≈ 20%. Are you really saying this could be the physical configuration, i.e. that the donor is over- flowing the Roche lobe by this amount ? LUCA IZZO: The results that I’ve presented comes out naturally from formulae and physical considerations based on fixed binary parameters, as MWD = 0.7 M�, q = 0.2, Porb = 1.8295 hrs, and considering the donor as a low mass M-type star. Following a work by Patter- son et al., we obtain the R2 value, while from Paczynski the value of L2. Their interpretation depends obviously on a binary parameter as q presented in the Uthas et al. work. Anyway, other works consider that value of q as the actual one working in T Pyx case. 268 Introduction The Nebular Phase of T Pyx The Configuration of the Binary System Final Considerations