arXiv:astro-ph/0104040v1 3 Apr 2001 olwrdsitglxe eg,Ueae l 1999a). al. et Umeda (e.g., red-shift high galaxies con- from red-shift thus “evolve” low Ia systems, to SNe binary not or progenitor whether the firming identify to sary ta.19;Rese l 98 eas fterams uni- almost their the of derive because groups luminosities, 1998) both maximum These al. luminosities. et maximum Riess form 1999; (Perlmutter promising al. parameters a et cosmological provide determining that been for indicators have tool distance Ia SNe good Recently, as stars. used of events explosive nous clyfo h ih uv hp LS fnab N Ia, SNe nearby of same (LCS) the shape that assume curve and light the from ically aueo N aadterpoeiossol etesame the be should progenitors their high- and between Ia SNe of nature e-hf ( red-shift rpittpstuigL using typeset Journal Preprint Astrophysical the in appear to ti ieyacpe htteepoigsa tefis itself star exploding the that accepted widely is It yeI uenve(N a r n ftems lumi- most the of one are Ia) (SNe supernovae Ia Type fcro n xgn ehv locaie h esnwyol Cr T only why I reason not. headings: Type do the Subject a systems clarified as three also re explode have other the We soon the in will but oxygen. WDs and and the limit carbon mass that of Chandrasekhar indicate the strongly near results to These Sgr). (V3890 yr asi lw ntewn,ec Dcnrti usata ato the o of rates part mass-increasing substantial net a obtained retain have can we WD outburst. Thus, each last wind, ends. the the burning in before blown phase is quiescent mass the during Sgr) V3890 1 V80Sr,idctn naeaems crto rate accretion mass average an indicating Sgr), (V3890 1 eue h neoems nteW.Ec neoems tteop declin the of at speed mass the envelope and Each WD. duration ∆ nova the be on the mass to determining envelope in the role reduces key a play u 1 but hroula uaa oe navr asv ht wr close lig dwarf visual white early massive very The a inst outbursts. on warping CrB i.e., ca model radiation-induced T runaway curve a of thermonuclear Light include peak a with also second together companion. We disk the giant accretion disk. reproduce o red the accretion consists and a the model star by and companion Our disk the of Sagittarii. accretion effects V3890 an g red and with a Scorpii, with (WD) those V745 are Ophiuchi, here studied RS objects four The parameters. . . z 2 35 − hoeia ih uvso orrcretnvei ubrtare outburst in novae recurrent four of curves light Theoretical N a nodrta hsmto ok,the works, method this that order In Ia. SNe ) × 1 M ± 10 TCB,0 CrB), (T z . 377 EURN OA SAPOEIO YTMO YEI SUPERNOVAE Ia TYPE OF SYSTEM PROGENITOR A AS NOVAE RECURRENT eateto at cec n srnm,Cleeo rsa Arts of College Astronomy, and Science Earth of Department WD 0 n low- and − . M 01 7 M ± A SOHUH UCAS—SSESWT E IN COMPANION GIANT RED A WITH SYSTEMS — SUBCLASS OPHIUCHI RS 1 = 1. T eateto srnm,Ki nvriy ioh,Kouhok Hiyoshi, University, Keio Astronomy, of Department ∼ M ⊙ E 0 tl mltajv 14/09/00 v. emulateapj style X . ⊙ introduction 01 yr . 3 L 37 z o 39 g.Otclytikwns hc lwfo h D uigth during WDs the from blow which winds, thick Optically Sgr. V3890 for × − . max 12 yboi tr:idvda R p,TCB 39 g,V4 Sco) V745 Sgr, V3890 CrB, T Oph, (RS individual stars: — symbiotic M L 1 ± iais ls oa,ctcymcvrals—sas supernovae stars: — variables cataclysmic novae, — close binaries: N a hrfr,i sneces- is it Therefore, Ia. SNe 10 max 1 r,R p) 0 Oph), RS yrs, (18 ⊙ × fSeI opeeyempir- completely Ia SNe of , 0 − . o SOhwt o ealct ( metallicity low with Oph RS for 01 10 LSrlto od nhigh in holds relation -LCS 6 M − M ⊙ 7 M ⊙ TCB,2 CrB), (T ⊙ o r,1 CrB, T for yr − 1 oapa nteAtohsclJournal Astrophysical the in appear to euok,Tko1380,Japan 153-8902, Tokyo Meguro-ku, R p) 0 Oph), (RS [email protected] . 9 × [email protected] × . 10 10 35 zm Hachisu Izumi aioKato Mariko − − ABSTRACT ± 6 7 M M 0 . . and ⊙ ⊙ 05 01 1 yr R p) 5 Oph), (RS × M − oe fotuss ahs ta.(00)hv elucidated have analyz- (2000a) light-curve By theoretical al. the et 1999). of Hachisu basis outbursts, a of on al. model data et new Munari these Matsumoto, 2000; ing 2000a; phase Hachisu mid-plateau phase & al. the rising et Kato, Hachisu through the (e.g., phase from eclipses decline including obtained final curve been light the has complete to a outburst 1999 time, the first 25, the of February For in 1999). outburst (Schmeer recorded 2000, sixth underwent Livio e.g., (see, etc 1992) summary). stars recent symbiotic novaefor Renzini recurrent 1985), & Truran 1984), Munari & (e.g., Sparks, Webbink Starrfield, 1984; con- (e.g., (e.g., been dwarfs Tutukov ever white & have double Iben merging objects include is of which system) types sidered, binary Several observed known. the the not thus However, Nomoto (and 1984). (e.g., star Yokoi binary companion & a Thielemann, in Nomoto, (WD) 1982; dwarf white accreting an 10 1 ⊙ Z n ftecniae,tercretnv Scorpii U nova recurrent the candidates, the of One M 5 r,V4 c) 1 Sco), V745 yrs, (52 − o SOhwt oa ealct ( metallicity solar with Oph RS for ˙ 0 = 7 acc M dSine,Uiest fTko Komaba, Tokyo, of University Sciences, nd ⊙ . ∼ 0) 1 004), atcmain .. ooa Borealis, Coronae T i.e., companion, iant -u ooaa2382,Japan 223-8521, Yokohama u-ku, yr lhuhalrepr fteenvelope the of part large a Although 0 × hdwn ffc ntecompanion the on effect shadowing a − . h D as WDs the f 4 oee ooti aiu physical various obtain to modeled 1 10 tcre r elrpoue with reproduced well are curves ht × urn oa aenwgonup grown now have novae current oteCadaehrms limit, mass Chandrasekhar the to V4 c) 0 Sco), (V745 . − 35 ia aiu sas estimated also is maximum tical 10 6 hw eodr maximum secondary a shows B uenv fteWsconsist WDs the if supernova a blt fteaceinds to disk accretion the of ability M ± − neoems fe hydrogen after mass envelope eymsiewiedwarf white massive very a f ⊙ 7 eas h idquickly wind the because e clto nldsreflection includes lculation 0 M . 01 V4 c) 3 Sco), (V745 . 1 ⊙ × M yr M 10 ˙ ⊙ − He . − 11 1 o 75So and Sco, V745 for 7 stars: — 8 r,TCrB), T yrs, (80 M ∼ × ⊙ 0 10 . yr 1 outbursts, e − × × − 7 Z I. . M 1 10 10 0 = 2 yrs, (28 ⊙ − − 7 6 yr . M M 02) − ⊙ ⊙ 1 2 Hachisu and Kato the outburst nature of U Sco much more accurately than wald 1995), V394 Coronae Austrinae with P = 0.757 the previous estimates (e.g., Starrfield, Sparks, & Shaviv days (Schaefer 1990), Large Magellanic Cloud Nova 1990 1988; Kato 1990). We summarize some critically impor- No.2, hereafter LMC RN, with P =unknown but P 1.3 tant values they have obtained. (1) The white dwarf mass days was suggested (Sekiguchi 1992), CI Aquilae∼ with in U Sco is 1.37 0.01M⊙. (2) The companion is a slightly P = 0.618 days (Mennickent & Honeycutt 1995), and ± evolved main-sequence star of 1.4—1.6 M⊙ that expands red giant (RG) companions, T Coronae Borealis with to fill the Roche lobe after a large part of the central hy- P = 227.6 days (Lines, Lines, & McFaul 1988), RS Ophi- drogen has been consumed and is now transferring mass uchi with P = 460 days, (Dobrzycka & Kenyon 1994), to the white dwarf in a thermally unstable manner. (3) recently revised to be P = 457 days (Fekel et al. 2000), The white dwarf is now growing in mass at an average rate V745 Scorpii with P =unknown, V3890 Sagittarii with −7 −1 of 1 10 M⊙ yr . Therefore, the white dwarf will P =unknown, where P is the orbital period of the binary ∼ × reach the critical mass (MIa = 1.378M⊙, in W7 model of in units of day. These are sometimes called T Pyx sub- Nomoto et al. 1984) for the SN Ia explosion in the quite class, U Sco subclass, and T CrB (or RS Oph) subclass, near future. They concluded that U Sco is a very strong respectively (e.g., Warner 1995). We have summarized candidate for the immediate progenitor of SNe Ia. Other these features in Table 1. recurrent novae are also expected to be a strong candidate Classical novae, sharing many aspects with the recur- for SNe Ia. rent novae, are also binary systems consisting of a white Thus, we have started our light curve analysis of recur- dwarf and a less massive main-sequence star that fills the rent novae in which the mass and the growth rate of white Roche lobe and transfers its mass to the white dwarf. The dwarfs are determined as accurately as possible and ex- white dwarf accumulates hydrogen-rich matter and ignites amine the possibility of an SN Ia explosion in the quite hydrogen atop the white dwarf when the hydrogen-rich en- near future. In this series of papers, details of our results velope on the white dwarf reaches a critical mass. A part on recurrent novae are presented, although a part of our of the hydrogen diffuses into the white dwarf before the results has already been published in several short papers ignition so that a very surface layer of the white dwarf (Hachisu & Kato 1999 for T CrB; Hachisu et al. 2000a,b is dredged up into the hydrogen-rich envelope and blown for U Sco; Hachisu & Kato 2000a for RS Oph; Hachisu & off in the outburst wind (e.g., Prialnik 1986; Kovetz & Kato 2000b for V394 CrA). Prialnik 1994). As a result, nova ejecta contain white Here, in the first paper of this series, we have examined dwarf matter, i.e., carbon and oxygen for carbon-oxygen four recurrent novae which have a red giant companion, white dwarfs or oxygen, neon, and magnesium for O-Ne- T CrB, RS Oph, V745 Sco, V3890 Sgr. Two of them (T Mg white dwarfs (see, e.g., Prialnik 1986; Kovetz & Pri- CrB and RS Oph) are reexamined with new observational alnik 1994; Warner 1995, and references therein). Thus, knowledge. We first summarize the nature of recurrent no- the white dwarf is eroded gradually after many cycles of vae and explain the physical reason why they become so nova outbursts, recurrence periods of which are usually massive as near the Chandrasekhar mass limit in 2. Then, longer than thousands of years. we show details of our numerical methods for calculating§ Recurrent novae show some characteristic differences: evolution of bloated white dwarf envelopes in 3 and for (1) heavy elements such as carbon, oxygen and neon are obtaining theoretical light curves in 4, because§ details of not enriched in ejecta but similar to the solar values (e.g., our numerical method have not been§ fully described yet. Livio & Truran 1992; Warner 1995, for summary), indi- New numerical results on determination of white dwarf cating that the white dwarf is not eroded; (2) very short masses for four recurrent novae are shown and discussed recurrence periods from a decade to a century theoreti- in 5 for T CrB, in 6 for RS Oph, in 7 for V745 Sco, and cally require very massive white dwarfs close to the Chan- in §8 for V3890 Sgr.§ Conclusions follow§ in 9. drasekhar mass limit (e.g., Fig. 9 of Nomoto 1982). These § § two arguments indicate that the mass of the white dwarfs 2. recurrent novae as a progenitor of type ia in recurrent novae increases toward the Chandrasekhar supernovae mass limit (e.g., Starrfield, Sparks, & Truran 1985; Star- rfield, Sparks, & Shaviv 1988; Livio & Truran 1992; Della Recurrent novae are characterized by nova eruptions Valle & Livio 1996). We expect that, if the white dwarf is with its recurrence time scale from a decade to a century, made of carbon and oxygen, the white dwarf ignites carbon its brightness amplitude larger than 6 magnitudes, and its at the center and explodes as a Type Ia supernova when −1 expansion velocity faster than 300 km s (e.g., Webbink it reaches the critical mass of MIa = 1.378M⊙ (Nomoto et al. 1987). They are considered to be binary systems in 1982; Nomoto et al. 1984). However, the evolutionary which the secondary companion star overflows the inner path to recurrent novae has not long been elucidated as critical Roche lobe and transfers matter onto the primary well as the reason why so massive white dwarfs can exist white dwarf (WD) component. When the hydrogen-rich in recurrent novae. envelope on the WD reaches a critical mass, hydrogen ig- nites to trigger a nova eruption (hydrogen shell-flash). We call this a thermonuclear runaway (TNR) event. 2.1. Accretion Wind Evolution vs. Common Envelope Recurrent novae are morphologically divided into three Evolution groups depending on its companion star: dwarf com- In this subsection, we briefly explain how white dwarfs panions, only one member, T Pyxidis with P = 0.076 in recurrent novae have grown its mass toward the Chan- days, (Schaefer 1990, 1992; Patterson et al. 1998); drasekhar mass limit. The reason why the standard binary slightly evolved main-sequence (MS) companions, U Scor- evolution theory (will be cited below) has failed to explain pii with P =1.23 days (Schaefer 1990; Schaefer & Ring- the evolutionary path to recurrent novae is mainly in the Recurrent Novae 3 theoretical prediction of a common envelope formation and binaries that finally produce an SN Ia. Figure 3 shows the ensuing spiral-in process, because these processes in- such two regions of the WD+MS systems (SSS channel) hibit the growth of white dwarfs in binary systems. As and the WD+RG systems (symbiotic channel). We have discussed in many previous papers on binary evolutions, added in this figure the final parameter regions just before it had been widely accepted that hydrogen-rich envelopes an SN Ia explosion occurs. Starting from the initial region on mass-accreting white dwarfs expand to a red giant size encircled by the thin solid lines, binary systems evolve and when its mass accretion rate, M˙ acc, exceeds a critical limit, can explode as an SN Ia in the regions enclosed by the thick −6 −1 M˙ cr, i.e., M˙ acc > M˙ cr 1 10 M⊙ yr (see, e.g., Fig. solid lines. Hachisu et al. (1999a,b) showed that some re- 9 of Nomoto 1982), and∼ easily× forms a common envelope current novae lie on the boarder of their parameter region (e.g., Nomoto, Nariai, & Sugimoto 1979). Once a common that can produce an SN Ia. Among six recurrent novae envelope is formed, two stars begin to spiral-in each other with known orbital periods, five fall in the final regions of due to viscous drag; thus, producing a double degenerate SN Ia progenitors just before the explosion. Only T Pyx system (e.g., Iben & Tutukov 1984; Webbink 1984). is far out of the regions of SNe Ia. However, Hachisu, Kato, & Nomoto (1996) found that white dwarfs begin to blow strong winds (v 1000 km 3. optically thick winds in decay phase of novae wind ∼ −1 ˙ & −6 −1 s and Mwind 1 10 M⊙ yr ) when the mass ac- In the thermonuclear runaway (TNR) model of nova × ˙ ˙ cretion rate exceeds the critical rate, i.e., Macc > Mcr, as explosions, WD envelopes expand greatly as large as ∼ illustrated in Figure 1; thus, preventing the binary from 100 R⊙ or more and then the photospheric radius gradu- collapsing. The binary does not spiral-in but keeps its ally shrinks to the original size of the white dwarfs (e.g., separation. Hydrogen burns steadily and therefore the RWD = 0.0032 R⊙ for MWD = 1.37 M⊙). Here, we call helium layer of the white dwarf can grow at a rate of this decay phase, i.e., the phase after the maximum ex- ˙ ˙ MHe Mcr. The other transferred matter is blown in pansion of the photosphere. In the decay phase of novae, the wind≈ (M˙ M˙ M˙ ), where M˙ is the wind 4 wind ≈ acc − cr wind the photospheric temperature increases from Tph 10 K mass loss rate. We call this accretion wind because it be- to 106 K with the bolometric luminosity being∼ almost gins to blow when the accretion rate exceeds the critical constant∼ (near the Eddington limit). Therefore, the main rate. emitting region moves from optical to soft X-ray through Therefore, as shown in Figure 2, we have to revise Figure UV. Correspondingly, the optical luminosity reaches its 9 of Nomoto (1982), which has been widely used in binary maximum at the maximum expansion of the photosphere evolution scenarios as basic processes of mass-accreting and then decays toward the level in quiescent phase. Since white dwarfs (e.g., Kahabka & van den Heuvel 1997). The the WD envelope reaches a steady-state in the decay phase most important difference of Figure 2 from old Nomoto’s of novae (e.g., Prialnik 1986; Kato & Hachisu 1994), we (1982) is that white dwarfs grow its mass in a much wider are able to approximate the development of the envelope parameter region. In our new diagram, the status of a structure with a unique sequence of steady-state wind so- mass-accreting white dwarf is specified by the following lutions having a different envelope mass (∆M) as shown three phases: (1) accretion wind phase (M˙ acc > M˙ cr); (2) by Kato & Hachisu (1994). steady shell-burning phase (M˙ std < M˙ acc < M˙ cr); and (3) Optically thick winds, blowing from WDs in the decay intermittent shell-flash phase (M˙ acc < M˙ std), where M˙ std phase of novae, play a key role in determining the nova means the lowest limit of mass accretion rate for steady duration, because the winds eject a large part of the enve- hydrogen shell burning. The new growing region of white lope mass and drastically reduce the fuel. The winds blow when the WD photosphere is larger than 0.05 0.1R⊙, dwarfs, M˙ acc > M˙ std, is much wider than old Nomoto’s i.e., when the photospheric temperature∼ is lower− than (1982) narrow strip, M˙ std < M˙ acc < M˙ cr. log Tph . 5.5, because the wind is driven by a strong peak at log T 5.2 in OPAL opacity (e.g., Iglesias & Rogers 2.2. Evolutionary Paths to Recurrent Novae and Type Ia ∼ Supernovae 1996), and log Tph 5.5 corresponds to the shoulder of the peak on the high∼ temperature side. This opacity peak Based on the new mechanism of Figures 1 and 2, strongly blocks photons and, as a reaction, matter is ef- Hachisu et al. (1999a,b) have found two paths to SNe Ia fectively accelerated to blow a wind deep inside the pho- (see also, Li & van den Heuvel 1997). In these two paths, tosphere (see, e.g., Figs. 2, 3, and 4 in Kato & Hachisu white dwarfs can accrete mass continuously from the com- 1994). −6 −1 panion and grow at a rate of M˙ cr 1 10 M⊙ yr . Assuming the solar abundance of heavy elements (Z = Therefore, candidates of SN Ia progenitors∼ × are systems at 0.02), we have calculated the steady-state envelope solu- the final stages on these paths; i.e., one is a supersoft X-ray tions and obtained sequences for WDs with various masses source (SSS) consisting of a white dwarf and a lobe-filling of MWD = 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.34, 1.35, more massive main-sequence star (WD+MS system, e.g., 1.36, 1.37, and 1.377M⊙. Here, we have adopted the up- Li & van den Heuvel 1997) and the other is a symbiotic dated OPAL opacity (Iglesias & Rogers 1996) with lower star consisting of a white dwarf and a mass-losing red gi- hydrogen contents of X =0.70, 0.50, 0.35, 0.15, 0.10, 0.08, ant (WD+RG system, e.g., Hachisu et al. 1996). Both 0.07, 0.06, 0.05, and 0.04. The numerical methods and var- the systems contribute to main parts of the SN Ia birth ious assumptions are essentially the same as those in Kato rate (Hachisu et al. 1999a,b) (see also Livio 2000, for a & Hachisu (1994) except the opacities. Some of them are recent review). plotted in Figure 4. Hachisu et al. (1999a,b) have followed many evolution- For RS Oph and LMC RN, we have also obtained the ary paths and obtained the initial parameter regions of envelope solutions for lower metallicities, i.e., Z = 0.001, 4 Hachisu and Kato

0.002, 0.004, and 0.01 with X = 0.70, 0.50, 0.35, 0.15, respectively. The efficiency of hydrogen shell-flashes (nova 0.10, 0.08, 0.07, 0.06, 0.05, and 0.04. Some of them are outbursts) is defined by also shown in Figure 5. ∆MHe ∆MHe ηH = , (4) 3.1. Evolution of White Dwarf Envelope ∆Mmax ≈ ∆Mwind + ∆MHe Hydrogen shell-flashes occur when the mass transfer rate where ∆Mmax is the envelope mass at the optical maxi- is below the lowest limit for the steady hydrogen shell- mum. −3 −1 ˙ ˙ Assuming a high constant rate of M˙ acc 10 M⊙ yr burning, i.e., Macc < Mstd. In this subsection, we describe ∼ such a rapid development of white dwarf envelopes. As in equation (1), we start the calculation. The envelope already mentioned above, optically thick winds blow in mass is quickly increased and the optical luminosity soon the decay phase of hydrogen shell-flashes. Since each wind reaches the observational maximum in a day or so. Then, −7 −1 we drop the rate to M˙ 1 10 M⊙ yr or zero and solution is a unique function of the envelope mass ∆M acc ∼ × for a given WD mass, we can follow the development of further follow the decay phase of nova. Therefore, our ris- the envelope structure by calculating the envelope mass, ing phase is not a realistic one but just a rough sketch. which is decreasing due to the wind mass loss at a rate Once the envelope has reached the optical maximum, our theoretical light curve in the decay phase represents a real- of M˙ wind(∆M) and hydrogen shell burning at a rate of istic one because the flow in the envelope is approximately M˙ (∆M), i.e., nuc in a steady-state as already mentioned above. d ∆M = M˙ acc M˙ wind M˙ nuc, (1) dt − − 3.2. White Dwarf Envelope in Steady Mass Accretion ˙ where Macc is the mass accretion rate of the WD. Directly ˙ integrating equation (1), we obtain the time-development If the mass accretion rate (Macc) does not vary, the enve- of the envelope, ∆M, and then the evolutionary changes of lope solution is obtained from d(∆M)/dt = 0 in equation physical quantities such as the photospheric temperature (1). The mass accretion rate is in a balance with the mass decreasing rate of the envelope, Tph(∆M), photospheric radius Rph(∆M), photospheric ˙ velocity vph(∆M), wind mass loss rate Mwind(∆M), and M˙ acc = M˙ wind + M˙ nuc. (5) nuclear burning rate M˙ nuc(∆M), all of which are a unique function of the envelope mass, ∆M. If the steady mass-accretion rate is given, we can directly We have shown, in Figures 4 and 5, the envelope mass (without integration) obtain the corresponding state of the ∆M, the photospheric radius R , the photospheric tem- envelope from Figures 4 and 5. The critical envelope-mass ph decreasing rate, i.e., at the cease point of wind, exactly co- perature Tph, and the photospheric velocity vph against the ˙ ˙ incides with the critical mass-accretion rate, i.e., at the be- mass decreasing rate of the envelope, Mwind + Mnuc. The ˙ ˙ ˙ −4 −3 ginning point of wind, Mcr (Macc)cr = (Mnuc)cr, which wind mass loss rate becomes as large as 10 —10 M⊙ ≡ −1 can be approximately expressed by yr at the maximum expansion of the photosphere Rph −6 −1 ∼ 100R⊙, and then decreases down to 1 10 M⊙ yr just −7 1.7 X MWD −1 × M˙ cr =5.3 10 − 0.40 M⊙ yr , (6) before the wind stops, i.e., when Rph 0.1R⊙. The pho- × X  M⊙ −  tospheric temperature increases from 1∼ 104 Kto 3 105 × × K. The photospheric velocity remains as high as 1000 for various WD masses (MWD) and hydrogen contents (X) km s−1 during the wind phase. ∼ with solar metallicity (Z =0.02). This relation is reduced When the envelope mass decreases to below the critical from our wind solutions for X =0.70, X =0.50, X =0.35, mass, the wind stops and the photosphere rapidly shrinks X = 0.15, and X = 0.10 (Z = 0.02) and has accuracy towards its original size before the outburst, for exam- within to 3—5%. ple, from Rph = 0.06R⊙ to 0.005R⊙ within a week for In the same steady mass-accretion mentioned above, we MWD = 1.37 M⊙. Shortly after the wind stops, steady- also obtain the lower limit of the mass accretion rate for state shell burning on the WD ends but hydrogen burning steady hydrogen shell-burning as itself is still continuing and supplying a part of the lu- −7 1.7 X MWD −1 minosity. In this phase, the envelope mass is decreased M˙ std =2.6 10 − 0.40 M⊙ yr only by nuclear burning. The rest of the luminosity comes × X  M⊙ −  from the thermal energy of the hydrogen envelope and the 0.5M˙ cr, (7) hot ash (helium) below the hydrogen layer. The thermal ≈ energy amounts to several times 1043 ergs for the case of for solar metallicity (Z = 0.02). This relation is reduced from our steady-state solutions for X = 0.70, X = 0.50, MWD = 1.37 M⊙, X = 0.05, and Z = 0.02 (for the 1999 outburst of U Sco). It can supply a bolometric luminosity X = 0.35, X = 0.15, and X = 0.10 (Z = 0.02) and has of 1038 ergs s−1 for about ten days or so. accuracy better than 10%. The mass lost by the wind, ∆Mwind, and the mass added to the helium layer of the WD, ∆MHe, are calculated from 4. theoretical light curves One of our binary models is graphically shown in Figure ∆Mwind = M˙ wind dt, (2) Z 6. A circular orbit is assumed. Our light curves include the and contributions of three components: photosphere of WD, red giant (RG) companion, and accretion disk (ACDK or ∆MHe = M˙ nuc dt, (3) Z DK). Recurrent Novae 5

4.1. White Dwarf Photosphere If the accretion disk around the WD blocks the light We assume a black-body photosphere of the WD enve- from the WD photosphere, it makes a shadow on the sur- face of the companion star. Such an effect is also included lope and calculate the V -magnitude with a window func- in our calculation as explained below. tion given by Allen (1973). The photospheric surface is divided into 32 pieces in the latitudinal angle (∆θ = π/32) 4.3. Accretion Disk Surface and into 64 pieces in the longitudinal angle (∆φ =2π/64) as illustrated in Figure 6. The contributions to the V light We have to include the luminosity coming from the ac- are summed up from each piece. We do not consider the cretion disk when it is irradiated by the WD photosphere. limb-darkening effect because it is as small as ∆V . 0.2 Here, we assume a circular accretion disk, the edge of mag. which is defined by We have assumed the luminosity of the WD given by ∗ Rdisk = αR1, (11) 1 GMWDM˙ acc where α is a numerical factor indicating a size of the accre- LWD = LWD,0 + , (8) tion disk, and R∗ the effective radius of the inner critical 2 RWD 1 Roche lobe given by Eggleton’s (1983) formula. where the first term is the intrinsic luminosity of the WD The surface of the accretion disk is strongly irradiated calculated by our wind and steady hydrogen shell burning and dragged outward by the strong wind from the WD solutions and the second term is the accretion luminosity during the nova outburst. Its thickness may be increasing (e.g., Starrfield, Sparks, & Shaviv 1988); RWD is the ra- outward by both the ablation and the Kelvin-Helmholtz dius of the WD (e.g., R =0.0032R⊙ for 1.37M⊙ WD). WD instability. We here approximate such a surface by The accretion luminosity is as large as 700 L⊙ for a mass ν −7 −1 ̟ accretion rate of M˙ 1 10 M⊙ yr onto a 1.37M⊙ acc ∼ × h = βRdisk , (12) WD. Therefore, the accretion luminosity itself is negligi- Rdisk  bly small during the steady hydrogen shell burning phase where h is the height of the surface from the equatorial 4 (LWD,0 & 10 L⊙) but may be important after the hydro- plane, ̟ the distance on the equatorial plane from the gen burning diminishes. In such a case, we recalculate the center of the WD, and β is a numerical factor showing the surface temperature of the white dwarf, i.e., degree of thickness. For comparison, we have calculated L the three cases of ν, i.e., ν = 9/8 for the case of stan- σT 4 = WD , (9) ph 4πR2 dard accretion disk model by Shakura & Sunyaev (1973), WD ν = 2 and ν = 3 for the case of flaring-up disk (Schandl, including the accretion luminosity. Meyer-Hofmeister, & Meyer 1997). The surface of the accretion disk is divided into 32 pieces 4.2. Companion’s Irradiated Photosphere logarithmically evenly in the radial direction and into 64 We have to include the contribution of the companion pieces evenly in the azimuthal angle as shown in Figure 6. star when it is strongly irradiated by the WD photosphere. The outer edge of the accretion disk is also divided into 64 The companion star is assumed to fill the inner critical pieces in the azimuthal direction and 2 (or 8) pieces in the Roche lobe as shown in Figure 6 for T CrB, but underfill vertical direction by rectangles. When the photosphere for RS Oph, V745 Sco, and V3890 Sgr as described later. of the WD becomes very small, e.g., Rdisk/Rph > 10, we Dividing the latitudinal angle into 32 pieces (∆θ = π/32) attribute the first 16 meshes to the outer region (from and the longitudinal angle into 64 pieces (∆φ =2π/64), we ̟ = Rdisk to ̟ = Rdisk/√10) to avoid coarse meshes in have also summed up the contribution from each area, but the outer part, and then 16 meshes to the inner region we neglect both the limb-darkening effect and the gravity- (from ̟ = Rdisk/√10 to ̟ = Rph), each region of which darkening effect of the companion star because the both ef- is divided logarithmically evenly. fects are negligibly small (∆V . 0.2 mag) for the strongly The surface of the accretion disk absorbs photons and irradiated companion. Here, we assume that 50% of the reemits in the same way as the companion does. Each area absorbed energy is reemitted from the hemisphere of the of the heated disk surface emits blackbody photons at its companion with a black-body spectrum at a local temper- local temperature, Tdisk. We assume that 50% of the ab- ature, that is, the efficiency of irradiation is ηRG = 0.5 sorbed energy is emitted from the surface while the other unless otherwise specified. is carried into interior of the accretion disk and eventually The original (non-irradiated) photospheric temperature brought into the WD, that is, the efficiency of irradiation of the companion star is obtained from the light curve is ηDK =0.5. fitting, for example, TRG,org = 3100 K for T CrB. The The non-irradiated temperature of the disk surface is irradiated surface temperature is roughly estimated by assumed to be determined by the viscous heating of the standard accretion disk model (Shakura & Sunyaev 1973). L σT 4 η WD cos θ + σT 4 , (10) Then, the original disk surface temperature is given by RG RG 4πr2 RG,org ≈ ˙ where r is the distance from the WD and cos θ is the incli- 4 3GMWDMacc σTdisk,org = . (13) nation angle of the surface. The irradiated temperature is 8π̟3 as high as TRG 5000 K for cos θ =0.5, MWD =1.37 M⊙ The viscous temperature is estimated to be Tdisk,org ∼38 −1 ∼ (LWD = 2 10 ergs s ), and r 150R⊙ (T CrB). 9,000 K at ̟ = 1R⊙ for MWD = 1.37M⊙ and M˙ acc = × ∼ −7 −1 Therefore, the effect of the irradiation becomes very im- 1 10 M⊙ yr . However, the viscous heating is neg- portant in the late phase of outbursts when the irradiated ligibly× small when the disk surface is strongly irradiated hemisphere faces towards the Earth. because of T 40,000—50,000 K at ̟ =1R⊙ for the disk ∼ 6 Hachisu and Kato luminosity of the WD, LWD = 50, 000L⊙, and the average we obtain an approximate normal unit vector of n to patch inclination of the surface, cos θ = 0.1, which is estimated ✷ABCD by from (e.g., Schandl et al. 1997) S n + S n + S n + S n n = 1 1 2 2 3 3 4 4 , (23) L 3GM M˙ S1 + S2 + S3 + S4 σT 4 η WD cos θ + WD acc . (14) disk ≈ DK 4π̟2 8π̟3 where each area of the triangles is corresponding to The outer rim of the accretion disk is not irradiated by the S1 = ABG, S2 = BCG, WD photosphere so that the temperature of the disk rim S = △CDG, S = △DAG, (24) 3 △ 4 △ is simply assumed, for example, to be 2000 K for T CrB. as shown in Figure 7, and the total area of this patch ✷ABCD is given by 4.4. Radiation-induced Warping S = S1 + S2 + S3 + S4. (25) Radiation-induced warping of accretion disks have been Then, patch i receives flux from patch j as much as suggested by Pringle (1996). Here, we have introduced 4 n r n r a warping surface calculated by Maloney, Begelman, & σTj Sj ( j ji)Si( i ij ) ∆L = · · , (26) Pringle (1996). In their formulation, when the tilt vector ij 2π r 4 | ji| of the accretion disk surface is defined by where r r r l(̟,t) = (sin β˜ cosγ, ˜ sin β˜ sinγ, ˜ cos β˜), (15) ij = G,j G,i, (27) − one of their stationary solutions, r = r r , (28) ji G,i − G,j W = βe˜ iγ˜, (16) as illustrated in Figure 8. We sum up all combinations of pair patches on the condition of can be analytically solved as (n r ) 0, and (n r ) 0, (29) sin(x) j · ji ≥ i · ij ≥ M(1, 2, 2ix)= eix , (17) x and obtain the total flux from all patches, i.e., σT 4S (n r )S (n r ) where x √̟. So, we have j j j · ji i i · ij ∝ Li = ξi,j,k 4 . (30) 2π rji x = c1π√̟, (18) Xj Xk | | and Here, we set ˜ sin(c1π√̟) 1, for not-blocked β = c2 , (19) ξi,j,k = (31) c1π√̟ n 0, for blocked γ˜ = c π√̟, (20) depending whether the line connecting the centers of grav- 1 ity is blocked by patch k as shown in Figure 8. where c1 and c2 are numerical factors indicating a degree The surface (patch i) absorbs photon flux and reemits of warping. Here, we adopt c1 = 0.5 and c2 = 1.0, un- by a black-body spectrum with a local temperature of Ti, less otherwise specified. Then, we have transformed− the which is determined by coordinates of the disk, i.e., equation (12), by a rotation σT 4S = η L + σT 4 S , (32) matrix of i i BB i i,org i where η is the efficiency of the irradiation, T the cos β˜ cosγ ˜ sinγ ˜ sin β˜ cosγ ˜ BB i,org − − unheated temperature of the surface. We adopt 50% effi- T =  cos β˜ sinγ ˜ cosγ ˜ sin β˜ sinγ ˜  . (21) − ciency, i.e., ηBB =0.5 unless otherwise specified. sin β˜ 0 cos β˜ Using the window functions, W (λ) and W (λ), for V -   V B Such an example is shown in Figure 6 for T CrB. As will and B-filters (Allen 1973), respectively, we have calcu- be discussed later, the condition of warping instability is lated V -, B-, and bolometric magnitudes of the patches, satisfied only in T CrB. i.e., ∞ ∆Ei(V) = (ni e)Si WV (λ)Bλ(Ti)dλ, (33) 4.5. Numerical Method · Z0 ∞ Total V light of the system is calculated from the contri- n e bution of each patch with a window (response) function of ∆Ei(B) = ( i )Si WB(λ)Bλ(Ti)dλ, (34) · Z0 the V -filter (e.g., Allen 1973). Each patch is almost quad- ∆E (bol) = (n e)S σT 4, (35) rangular but, strictly speaking, not quadrangular for the i i · i i Roche geometry of the companion star and for the warping respectively. Here, the unit vector e denotes the direction surface of the accretion disk, because the four vertices are from the patch to the Earth, Bλ(Ti) the Planck function not on a plane. Therefore, we define approximately the for the temperature of Ti. We have summed up the patches normal unit vector to each patch as follows: the center of having (ni e) > 0, i.e., gravity of each patch can be defined by · E(V )= ∆E (V ), for (n e) > 0 (36) 1 i i · r = (r + r + r + r ) , (22) Xi G 4 A B C D E(B)= ∆E (B), for (n e) > 0 (37) where G denotes the center of gravity of the patch, and i i · A, B, C, and D are the apexes of each quadrangle as illus- Xi trated in Figure 7. Defining n1 as the normal unit vector E(bol) = ∆Ei(bol), for (ni e) > 0. (38) to ABG, n2 to BCG, n3 to CDG, and n4 to DAG, · △ △ △ △ Xi Recurrent Novae 7

−1 We have finally calculated the V -magnitude, mV , B- s (cool component) indicated the mass ratio of q = magnitude, mB, and the bolometric magnitude, mbol, as M2/M1 1.4 (Kraft 1958; Paczynski 1965; Kenyon & ∼ ◦ 5 d Garcia 1986). Adopting the inclination angle of i = 68 , mV = ( log E(V )+33.582) + mV,0 + 5 log( ), (39) the resulting masses were M = M 1.9M⊙ and 2 − 10 1 hot ∼ M = M 2.6M⊙, thus placing the hot compo- 5 d 2 cool ∼ mB = ( log E(B)+33.582)+ mB,0 + 5 log( ), (40) nent above the Chandrasekhar limit. Therefore, Webbink 2 − 10 (1976) and Webbink et al. (1987) proposed an outburst 5 d mechanism of T CrB based on their main-sequence accre- m = ( log E(bol) + 33.582)+4.75+5 log( ), (41) bol 2 − 10 tor model: the outburst was caused by the transfer of a chunk of matter ejected by the giant and, then, by the where d is the distance to the object in units of pc, m V,0 subsequent shock dissipation of the kinetic energy. This is a constant determined to satisfy accretion event model was investigated further by 3D nu- mV = mbol, at T = 6500 K, (42) merical simulations (Cannizzo & Kenyon 1992; Ruffert, and mB,0 is also a constant determined to satisfy Cannizzo, & Kenyon 1993). However, Selvelli, Cassatella, & Gilmozzi (1992) have mB = mV , at T = 9600 K (A0), (43) been opposed to the main-sequence accretor model from so that we use their analysis of IUE data in quiescence, which indicates mV,0 =4.52, (44) the existence of a mass-accreting white dwarf (WD) as fol- lows: (1) the bulk of the luminosity is emitted in the UV and with little contribution to the optical; (2) the presence of m =4.91. (45) B,0 strong He II and N V emission lines suggests tempera- In addition, in order to compare our brightness with the tures of the order of 105 K; (3) the rotation broadening of recent CCD photometry, we also calculate the Cousins the high-excitation lines is larger than 1300 km s−1 with- Rc- and Ic-magnitudes with the window functions given out inclination correction. Moreover, the detection of X- by Bessell (1990) along the same way as those described rays (Cordova, Mason, & Nelson 1981) and the presence above for V - and B-magnitudes. Then, we adopt of flickering in the optical light curves (e.g., Ianna 1964; Zamanov & Bruch 1998) are also naturally explained in mR,0 =4.91, (46) terms of accretion onto a WD. and The mass accretion rate estimated by Selvelli et al. mI,0 =4.50. (47) −8 (1992) in quiescence is very high (M˙ acc 2.5 10 M⊙ −1 ∼ × The most time-consuming part of the calculation is pro- yr ) and is exactly required by the thermonuclear run- portional to cubic of the patch number N, that is, N 3. away (TNR) theory to produce a TNR event every 80 yr on Therefore, it is hard to increase the patch number∝ sub- a massive (& 1.3M⊙) WD. They also pointed out that the stantially more than the present maximum case, i.e., measurement of the emission-line radial velocities (Kraft 1958) was intrinsically difficult and uncertain. Thus, the N = 32 64+32 64+(32 64 2+8 64)= 8704, (48) × × × × × 1866 and 1946 outbursts can be interpreted in terms of a because it takes about a half day on the Compaq (old TNR event on a very massive WD (e.g., Starrfield et al. DEC) Alpha 21264 (500 MHz) machine but about two 1985). days on the Intel Pentium III (500 MHz) machine to com- pute one time sequence for the 1999 U Sco outburst, about 5.2. Very Massive White Dwarf in T CrB 3000 steps for the mV and mB light curve calculations (Hachisu et al. 2000a). Rapid decline rates of the light curves indicate a very massive WD close to the Chandrasekhar limit. Assuming 5. t coronae borealis the solar metallicity (Z =0.02) of the WD envelope, Kato (1995, 1999) calculated nova light curves for the WD mass 5.1. Nature of T CrB Outbursts of MWD =1.2, 1.3, 1.35 and 1.377 M⊙, and found that the T Coronae Borealis (T CrB) is one of the well observed light curve of the 1.377M⊙ model is in better agreement recurrent novae and characterized by a secondary maxi- with the observational light curve of T CrB than the other mum occurring 150 days after the primary peak. His- lower mass models. torically, T CrB∼ bursted twice, in 1866 and 1946, with the Recently, other observational supports of a massive WD light curves very similar each other (e.g., Pettit 1946a, b, in T CrB have been reported. Belczy´nski & Mikolajewska c, d). Large ellipsoidal variations of optical and infrared (1998) derived a permitted range of binary parameters, light curves in quiescence suggest that an M3 giant fills its M = 1.2 0.2M⊙ and q = M /M 0.6, i.e., WD ± RG WD ∼ Roche lobe (e.g., Bailey 1975; Lines et al. 1988; Yudin & MRG 0.6—0.8M⊙, from amplitude of the ellipsoidal Munari 1993; Leibowitz et al. 1997; Shahbaz et al. 1997; variability∼ and constraints from the orbital solution of Belczy´nski & Mikolajewska 1998). From the ellipsoidal M-giants. In Shahbaz et al. (1997), a massive WD variations, the orbital period and ephemeris were deter- of MWD = 1.3 2.5M⊙ is suggested from the infrared mined to be HJD 2,431,931.05 + 227.67E at the epoch of light curve fitting.− Moreover, the measurement of the spectroscopic conjunction with the M-giant in front (Lines radial velocities (Kraft 1958) has been revised includ- et al. 1988). ing the effect of gas-streams around the hot component There have been debates on the nature of the hot com- (Hric et al. 1998). Hric et al. (1998) estimated the ponent of the binary system. Radial velocity curves of masses, MWD =1.2 0.2 M⊙ for the hot component and −1 ± K 31 km s (hot component) and K 23 km M = 1.38 0.2 M⊙ for the M-giant, assuming the in- 1 ∼ 2 ∼ RG ± 8 Hachisu and Kato clination of i = 68◦. Combining these permitted ranges of Schandl et al. 1997). However, Selvelli et al. (1992) es- the WD mass in T CrB, we may conclude that a mass of timated the optically thick region of the accretion disk in the WD is between M =1.3 1.4M⊙, which is very con- quiescence is as small as R . 1 R⊙ because the disk lu- WD − disk sistent with the light curve analysis MWD 1.37 1.38M⊙ minosity contributes mostly to the satellite UV. Belczy´nski by Kato (1999). ∼ − & Mikolajewska (1998) pointed out that the optical con- tribution of the accretion disk is clearly visible in the U, B, 5.3. Secondary Maximum of Outbursts and V bands in 1981–85 (optically high state) while it is practically absent in 1990–94 (optically low state) as first The secondary maximum in outbursts is rarely observed shown by Payne-Gaposchkin & Wright (1946). Thus, Bel- in fast novae and in other recurrent novae. Selvelli et al. czy´nski & Mikolajewska derived the disk radius (hot spot) (1992) suggested a possibility of an irradiated stationary of R 0.1a (20 R⊙) in the optically high state. shell around the binary system, although the presence of disk Therefore,∼ we here adopt a small size of the optically such a shell is just speculation. Recently, Hachisu & Kato thick disk, i.e., R 1 R⊙ just before the outburst. A (1999) proposed another mechanism of the second peak, disk ∼ large accretion disk of 50 R⊙ requires that (1) the ir- i.e., radiation-induced tilting disk around a massive WD. radiation efficiency is too∼ small and that (2) the disk is The main results of their analysis are: (1) the first peak inclined from too early phase of the outburst. If we intro- is naturally reproduced by a fast developing photosphere duce a small size of the accretion disk, we do not need to of the WD envelope based on the TNR model of a very assume such unnatural assumptions. This accretion disk massive WD, M 1.35M⊙; (2) the second peak is not WD satisfies the warping condition (49) so that we adopt an fully reproduced by∼ an irradiated M-giant model as sim- analytic form proposed by Maloney et al. (1996). ply estimated by Webbink et al. (1987); (3) instead, the second peak can be well reproduced if an irradiated tilting 5.5. Light Curve of Early Phase accretion disk around the WD is introduced together with the partly irradiated M-giant companion. The very fast decay during the first 10 days can be well Such a radiation-induced, tilting instability of an accre- reproduced only by a WD photosphere of a very massive WD, M 1.35—1.37 M⊙, as already shown by Kato tion disk has been suggested by Pringle (1996) when a WD ∼ central star is as luminous as the Eddington limit. This (1999) and Hachisu & Kato (1999). We have found that radiation-induced instability sets in if the condition the decay rate of the 1.37 M⊙ light curve is in better agree- ment with the observational points than the others. The 1/2 M˙ acc Rdisk Lbol difference from Hachisu & Kato (1999) comes from the − − . − 10 7 M⊙ yr 1 10 R⊙  2 1038 ergs s 1  difference of the disk size which affects the light curve 10— 1/2 × −1/2 30 days after maximum. Hachisu & Kato (1999) adopted R M WD WD (49) a large accretion disk which in general contributes more × 0.003 R⊙  1.37 M⊙  to the luminosity. Therefore, they reduced the WD lumi- is satisfied (Southwell, Livio, & Pringle 1997), where Rdisk nosity to be consistent with the observational points by is the radius at the edge of the optically thick accretion occulting the WD photosphere with the accretion disk. In disk. Selvelli et al. (1992) estimated the accretion rate of the present model, on the other hand, the small disk can- −7 −1 not occult a large part of the WD photosphere during the T CrB as low as M˙ acc 0.25 10 M⊙ yr , which meets condition (49) above∼ even if× the optically thick radius is day 10—30. Thus, we need a more massive WD to gen- erate a faster decline. Therefore, we have determined the as small as Rdisk 1 R⊙. Therefore, it is likely that the ∼ WD mass of M = 1.37 0.01 M⊙ in T CrB, giving a radiation induced instability grows during the outburst in WD ± T CrB. The growth timescale of warping is estimated by much better fit with the observational points even in the Pringle (1996) and Livio & Pringle (1996) as day 0—10. Deutsch (1948) suggested from the observations of Ash- −1/2 Mdisk MWD Rdisk brook (1946), Wright (1946), and Pettit (1946a) that the τprec 30 − ≃ 10 8 M⊙  1.37 M⊙   R⊙  color index of T CrB was about zero February 12 (4 days −1 after maximum). Our model indicates the photospheric Lbol day , (50) color of (B V )0 = 0.2 (Tph 18, 500 K, Rph 18 R⊙), × 2 1038 ergs s−1  − − ∼ ∼ × which is consistent with the observational color because where we assume that the α-parameter in the Shakura- the reddening is about E(B V )=0.15. − Sunyaev (1973) standard accretion disk is α 0.1. Thus, The mass of the companion is estimated as follows: ∼ the growth timescale is short enough to excite warping of Adopting MWD =1.37 M⊙ and the radial velocity of the −1 the accretion disk. cool component K2 = 23.32 km s (Kenyon & Garcia 1986), i.e., 3 3 5.4. Revised Model of Warping Accretion Disk MWD sin i f(M)= 2 =0.30 0.01M⊙, (51) We have revised Hachisu & Kato (1999) model to be (MWD + MRG) ± consistent with the observations described below. The we have MRG = 1.30, 1.25, 1.16, and 1.0 M⊙, for the main difference is in the accretion disk, which plays an inclination angle of i = 70◦, 68◦, 65◦, and 60◦, respec- essential role in the mid/late phase of the V light curve. tively. Following Belczy´nski & Mikolajewska (1998), Hachisu & Kato (1999) assumed that the size of the ac- q M /M 0.6 and an inclination angle of i 60◦, ≡ RG WD ∼ ∼ cretion disk is as large as 0.7 times the Roche lobe size, we have assumed the red giant mass of MRG =1.0 M⊙ (i = ◦ i.e., R 50 R⊙. This is because the mass transfer is 60 ). So we show here only the results for M =1.0M⊙. disk ∼ RG not a wind-fed type but a Roche lobe overflow type (e.g., In this case, the separation is a = 209.2R⊙, the effective Recurrent Novae 9 radius of the inner critical Roche lobe for the WD compo- 5.6. Light Curve Model of the Second Peak ∗ nent is R1 = 85.0R⊙, the effective radius for the red giant ∗ Warping (or wobbling) instability of the accretion disk companion is R = R = 73.6R⊙. We have examined the 2 2 begins to grow after the strong wind stops 60 days af- other two cases of the red giant mass, i.e., MRG =1.25 M⊙ ∼ ◦ ◦ ter maximum (on HJD 2,431,920) and, as a result, the V (i = 68 ) and MRG = 1.16 M⊙ (i = 65 ), and obtained ◦ light attains another peak. We expect that the radiation similar light curves to that of M =1.0 M⊙ (i = 60 ). RG induced warping of the accretion disk grows in a time scale Figure 9 shows the best fitted model to the 1946 out- of a few tens of days as derived in equation (50) and further burst together with the observational points. Here, we that the mass transfer from the red giant to the WD in- assume the solar metallicity (Z = 0.02) and hydrogen creases because the red giant surface is strongly irradiated content of X = 0.70 for the WD envelope. We also as- and the Roche lobe overflow is enhanced. Such enhance- sume that the V magnitude reached its maximum on HJD ment of the mass transfer certainly increases the radius of 2,431,861. The dotted curve denotes the total V light the optically thick region of the accretion disk up to 10— from the WD photosphere and the non-irradiated red gi- ∼ 20 R⊙ as suggested by Belczy´nski & Mikolajewska (1998). ant at TRG,org = 3100 K. Ellipsoidal light variations can Thus, we conclude that the radius of the optically thick ac- be seen in the late stage (dotted). The effective tempera- cretion disk increased from α =0.007 (Rdisk 0.6 R⊙) to ture of the M-giant is estimated to be TRG,org = 3100 K ∼ α = 0.07 (Rdisk 6 R⊙) during the day 60—150, by from our light-curve fitting, which is a little bit lower than fitting with the observational∼ points as shown∼ in Figure 9. T = 3560 K by Belczy´nski & Mikolajewska (1998), but eff The second peaks of the 1866 and 1946 outbursts are roughly consistent with T = 3200 K by Shahbaz et al. eff very similar. This is explained as follows: (1) The first (1999). peak ends up when the wind stops about 60 days after If we include irradiation of the companion, the light maximum. This period ( 60 days), between the opti- curve changes to the dash-dotted curve. This irradiation cal peak and the termination∼ of the strong wind, is al- effect can be well seen both 10—40 days after the optical most uniquely determined by the WD mass. (2) The sec- maximum and 120—260 days after the optical maximum. ond peak begins to rise when the warping of the accretion The irradiation effect was once examined and discarded by disk grows, which is about 120 days after maximum. The Webbink et al. (1987) as a main mechanism of the sec- growth time is estimated to be about a few times 30 days ondary maximum, mainly because its maximum largely as in equation (50). This growth time depends on the mass deviates from the observational secondary maximum as (Mdisk) and the optically thick radius (Rdisk) of disk, which seen in the figure. are determined mainly by the mass transfer rate from the The solid curve includes further the irradiated accre- red giant companion. We can naturally expect the similar tion disk. Here, we have started the calculation by as- time scales between the 1866 and 1946 outbursts, because suming the disk size of α = 0.007 in equation (11), i.e., the irradiation of the red giant by the WD photosphere is R =0.6 R⊙. disk almost the same. (3) Thus, we are able to obtain quite The photospheric radius expands to Rph 100 R⊙ in similar light curves if the orbital phases of the 1866 and order to fit the peak luminosity of the V light.∼ We have 1946 outbursts happen to be the same (see the light curves derived the distance modulus of (m M) = 10.2. The − V of the 1866 and 1946 outbursts, Pettit 1946d). Acciden- light curve during the first 10 days of the 1946 outburst tally, the two orbital phases were almost similar to each can be well reproduced only with the fast-developing WD other. photosphere of M = 1.37 M⊙ as already shown by WD The precessing period of the warping disk is similar to Kato (1999). The envelope mass of the WD at the op- −6 the growth time of warping (e.g., Livio & Pringle 1996). tical maximum is ∆M = 3.3 10 M⊙ from our wind About 60% faster precessing angular velocity than the or- solutions. Then, the mass accretion× rate was estimated to bital motion is required from the phase relation between ˙ −8 −1 be Macc = 4 10 M⊙ yr during the quiescent phase the rising shoulder of the second peak near HJD 2,431,960, between 1866× and 1946. This value is roughly consistent −8 −1 a small dip near HJD 2,432,090, and then a small bump with the estimation of 2.5 10 M⊙ yr by Selvelli et × near HJD 2,432,120 as shown in Figure 9. These dips are al. (1992). caused by a large shadow on the companion cast by the After the WD photosphere shrinks to smaller than accretion disk. In other words, the accretion disk blocks R 0.6 R⊙, the accretion disk reappears but its contri- disk ∼ the light from the WD photosphere and makes its shadow bution to the V light is small as shown in Figure 9, because on the companion surface. Finally, we have determined radiation emits mainly in UV. We assume a flaring-up rim Pprec = 145 days by fitting. This period is similar to the of the accretion disk, i.e., β =0.30 in equation (12), dur- 164 days precessing period of SS 433. ing the strong wind phase for the 1946 outburst of T CrB, Webbink et al. (1987) summarized the observational since Hachisu et al. (2000a) suggested that the rim of the development of the 1946 outburst, showing the color in- accretion disk is flaring-up during the 1999 outburst of U dex at the second peak of B V 0.4 after correcting the Sco, Then, the V light is reduced 20—30 days after maxi- reddening of E(B V )=0.15.− They∼ further pointed out mum because the disk rim makes shadow on the red giant that the reflection− by the M-giant is about B V 1.35, component. The warping of the accretion disk does not therefore, rejected the possibility of reflection− by∼ the M- occur yet, because the strong wind suppresses the warp- giant. Instead, they suggested that the much bluer color ing. The spherically symmetric wind momentum is much index of 0.4 is typical of the accretion disks seen in out- larger than the radiation momentum as shown by Kato & bursting∼ accretion-powered symbiotic stars (e.g., Kenyon Hachisu (1994). & Webbink 1984). In our warping-disk model, the color index is ranging from B V = 0.15 (on HJD 2,432,000) − 10 Hachisu and Kato to B V = 0.35 (on HJD 2,432,050), which is reasonably that the mass of the WD in T CrB has been increased up consistent− with the observational color at the second peak to the present mass from a much lower initial mass (see because the reflection effect by the warping disk dominates e.g., the evolutionary path to SNe Ia proposed by Hachisu during the second peak rather than the irradiated M-giant. et al. 1999a). Moreover, it is theoretically predicted that Ever since the emission lines of O III, N III, He II,[O III], an upper limit of degenerate carbon-oxygen core is about [Ne III], and [Ne V] first appeared, they have been always 1.07 M⊙, i.e., MC+O,0 < 1.07 M⊙ (Umeda et al. 1999b). present until two years after the outburst, except in early Thus, we may conclude that T CrB will explode as a Type June, at the most, from May 15 to July 17, 1946 ( 105— Ia supernova in quite a near future if the hot component 135 days after maximum, Sanford 1949). In our∼ model- is a carbon-oxygen white dwarf. ing, the warping-up disk occulted both the WD and the innermost region of the disk during these period, being 5.8. Summary of T CrB Outburst consistent with the disappearance of the emission lines. The main features of the T CrB outbursts can be well understood with a thermonuclear runaway model on a very 5.7. The Distance, Envelope Mass of the White Dwarf, massive white dwarf having a warping accretion disk. We and Possibility of SN Ia Explosion have summarized the following important points from the The distance to T CrB is estimated to be d = 0.94 fitting of the modeled light curve with observation: kpc with the absorption of AV = 0.35 (Harrison, John- (1) The early phase ( 0—10 days) visual light curve of T son, & Spyromilio 1993) because the distance modulus CrB outbursts can be∼ well reproduced only by a bloated is (m M)0 = (m M)V AV = 10.2 0.35 = 9.85. white dwarf photosphere of MWD = 1.37 0.01 M⊙ for This is− consistent with− both Harrison− et al.’s− estimation of solar metallicity of Z =0.02. This is because± decline rates d = 1100 pc, Belczy´nski & Mikolajewska’s (1998) estima- of optical light curves depends sensitively on white dwarf tion of d = 960 150 pc, and Bailey’s (1981) estimation of mass especially for white dwarf masses near the Chan- d = 1180 pc. Harrison± et al. determined the infrared spec- drasekhar mass limit. Thus, we are able to specify the tral type of the cool companion as M3 III. Assuming the mass of white dwarf with a great accuracy (three digits). absolute magnitude of M = 5.3 (M3 III), they obtain (2) The ensuing early phase ( 10—40 days) visual light K − ∼ the distance modulus of (m M)0 = mK MK AK = curve cannot be reproduced by any combination of a white 4.79 ( 5.3) 0.04 = 10.05,− where they− used a− relation dwarf photosphere and a non-irradiated red giant. It re- − − − between the absorptions of AK and AV given by Rieke & quires an irradiation effect of the red giant component. Lebofsky (1985), i.e., AK /AV =0.112. This heated-up red giant surface is consistent with the The envelope mass of the WD at the first peak is very early appearance of M-type spectrum (TiO bands) −6 ∆M = 3.3 10 M⊙ from our wind solutions. There- shortly after the outburst, indicating the luminosity of M- fore, the mass× accretion rate can be estimated to be type giant much brighter than 10th magnitude (Deutsch −8 −1 M˙ acc = 4 10 M⊙ yr if the envelope mass had been 1948). accumulated× during the quiescent phase between 1866 and (3) Through the mid and late phase ( 40—300 days), the ∼ 1946. This value is roughly consistent with the estimation secondary maximum can be well modeled if we introduce −8 −1 of 2.5 10 M⊙ yr by Selvelli et al. (1992). The wind the radiation-induced warping instability of the accretion × ˙ −4 disk, which sets in after the wind stops 60 days after mass loss rate at the first peak is Mwind =3.7 10 M⊙ ∼ −1 × maximum. Introducing an analytic form of warping (Mal- yr while the nuclear burning rate is as small as M˙ nuc = −7 −1 oney et al. 1996), we obtain the size of the accretion disk 9.3 10 M⊙ yr . The wind mass loss rate is quickly de- × ˙ −6 −1 as large as Rdisk 6 R⊙ and the precessing period of the creasing to Mwind 1 10 M⊙ yr about 30 days after warping disk as fast∼ as 140 days to reproduce the late the first peak. Then,∼ 90%× of the envelope mass has been phase light curve. ∼ blown in the wind. The residual 10% is left on the WD as (4) The envelope mass of the bloated white dwarf is a part of the helium layer. Therefore, the net growth rate −6 −9 −1 3.3 10 M⊙, suggesting the mass accretion rate of of the WD is M˙ =4 10 M⊙ yr . × −7 −1 He × 0.4 10 M⊙ yr between the 1866 and 1946 outbursts. At the end of this section, we discuss whether or not (5)× The wind carries away about 90% of the envelope mass, T CrB will explode as a Type Ia supernova in the near so that the residual 10% can be accumulated to the helium future. The conditions for SN Ia explosions are summa- layer on the white dwarf. As a result, the white dwarf rized as follows (see, e.g., Nomoto & Kondo 1991): (1) −8 −1 grows at a rate of 0.4 10 M⊙ yr . the accretion rate of the C+O WD should be higher than (6) T CrB will certainly× explode as a Type Ia supernova the critical accretion rate under which helium detonation in quite a near future if the white dwarf consists of carbon occurs, i.e., and oxygen. −8 −1 M˙ acc > M˙ He−det 1 10 M⊙ yr , (52) ∼ × 6. rs ophiuchi for Population I stars; (2) the initial mass of the C+O WD at the beginning of the mass transfer is less massive than RS Oph is also one of the well-observed recurrent novae. It is characterized by a long orbital period of 460 days (Do- M < 1.2 M⊙. (53) C+O,0 brzycka & Kenyon 1994, but recently 457 days by Fekel et T CrB satisfies the first condition as has already been sug- al. 2000), a relatively short recurrence period of 10—20 −8 −1 ∼ gested to be M˙ acc 4 10 M⊙ yr (even if M˙ acc yrs compared with 80 yrs of T CrB. It has been suggested −8 −1∼ × ∼ 2.5 10 M⊙ yr by Selvelli et al. 1992). For the sec- that a companion (M-giant) star is underfilling the Roche ond× condition, we cannot estimate the initial WD mass of lobe and losing its mass by massive stellar winds (e.g., T CrB. However, evolutionary scenario strongly suggests Dobrzycka et al. 1996). RS Oph underwent five recorded Recurrent Novae 11 outbursts (in 1898, 1933, 1958, 1967, and 1985), with the 6.1. The Model of the 1985 Outburst light curves being very similar to each other (e.g., Rosino Assuming an extremely massive white dwarf of M = 1987). The latest (1985) outburst has been observed at WD 1.377 M⊙ just prior to the SN Ia explosion, we have cal- all wave lengths from radio to X-rays (see papers in Bode culated light curves for various metallicities of Z =0.001, 1987). 0.004, 0.01, and 0.02 as shown in Figure 10. The V light of Although there had been intense debates on the mecha- the first five days (HJD 2,446,092 – 2,446,096) can be re- nism of RS Oph outbursts (e.g., Livio, Truran, & Webbink produced only by the bloated WD photosphere as already 1986; Webbink et al. 1987), various observational as- shown in Hachisu & Kato (2000a). The case of Z =0.001 pects favor TNR models on a very massive WD (see, e.g., is too slow to be compatible with the observations. Thus, Anupama & Mikolajewska 1999, for a recent summary). we may conclude that the metallicity is Z & 0.004. Here, Rapid decline rates of the V light curves also indicate a we assume Z =0.004 for RS Oph as a lower limit of metal- very massive WD close to the Chandrasekhar limit. Kato licity. For the set of M =1.377 M⊙ and Z =0.004, we (1991) has first calculated RS Oph outburst light curves WD have obtained the envelope mass at the optical maximum, for the WD masses of 1.33, 1.35, 1.36 and 1.37 M⊙ and −6 ∆M = 2.2 10 M⊙, which indicates a mass accretion found that the light curve of the 1.36M⊙ model is in better × −7 −1 rate of 1.2 10 M⊙ yr during the quiescent phase be- agreement with the observational light curve of RS Oph tween the 1968× and the 1985 outbursts. than the other mass models. The irradiation effects of the accretion disk and the red Since Kato (1991) adopted an old opacity in her calcula- giant (RG) are necessary to reproduce the optical and UV tion, Hachisu & Kato (2000a) recalculated the light curves light curves of the outburst as already shown by Hachisu of RS Oph outbursts with the new opacity (e.g., Iglesias & Kato (2000a). Dobrzycka & Kenyon (1994) obtained & Rogers 1996), including irradiation effects of both an the radial velocities of the M-type feature and the A-type M-giant companion and an accretion disk around the WD. feature but concluded that the A-type absorption feature Their main results are summarized as: (1) the first 4 days does not associate with any of the components. Thus, RS light curve of the outbursts can be well reproduced with Oph is a single line binary with a TNR model based on a M = 1.35 0.01M⊙ and an WD ± 3 3 apparent distance modulus of (m M)V = 11.09. (2) MWD sin i − 2 0.1 M⊙, (54) To reproduce the V and UV light curves in the late stage (MWD + MRG) ≈ (t 4—100 days after the optical maximum), the contri- ∼ from the semiamplitude of the radial velocity KRG = butions of the irradiated M-giant and the irradiated ac- 12.8 2.0 km s−1 and the orbital period of P = 460 cretion disk must be introduced. (3) The distance to RS days.± They suggested a relatively low inclination angle Oph is estimated to be 0.6 kpc both from the V and UV ◦ of i 30—40 . If MWD = 1.377 M⊙, then we obtain light curve fittings with the absorption of A = 2.3. (4) ∼ ◦ ◦ V MRG = 0.5 M⊙ for i = 30 , MRG = 0.7 M⊙ for i = 33 , About 90% envelope mass of the WD has been blown off in ◦ ◦ MRG =0.9 M⊙ for i = 36 , MRG =1.0 M⊙ for i = 37 , or the wind but the residual 10% has been left on the WD. ◦ MRG = 1.2 M⊙ for i = 39 . We assume MRG = 1.0 M⊙ Thus, Hachisu & Kato (2000a) has concluded that RS and i = 37◦ in the present paper, although we have al- Oph will explode as an SN Ia in the near future if the hot ready confirmed that the other set of MRG =0.7 M⊙ and component is a carbon-oxygen white dwarf. i = 33◦ reproduces almost the same light curves. In Hachisu & Kato’s (2000a) calculations, they assumed For the cool component of the RS Oph system, as sug- the solar metallicity (Z = 0.02) of the M-giant enve- gested by Dobrzycka et al. (1996), the RG lies well within lope. However, it has recently been suggested that the the inner critical Roche lobe, i.e., cool component of RS Oph is a metal-poor giant. Scott et ∗ al. (1994) obtained the carbon abundance of [C/H]= 3 RRG = γR2 (γ < 1), (55) − ∗ for RS Oph 143 days after the 1985 outburst and of where R2 is the effective radius of the inner critical Roche [C/H]= 4 in 1992. Similar depletions of carbon, oxy- lobe for the RG component and γ is a numerical factor. gen, neon,− sulphur, and iron with respect to solar values Dobrzycka et al. (1996) adopted γ =0.4 for the distance are suggested by Contini, Orio, & Prialnik (1995). They of d =1.5 kpc to RS Oph. On the other hand, Hachisu & obtained one or two tenths of the solar values for C, O, Kato (2000a) estimated the distance of d =0.6 kpc to RS Ne, S, and Fe, in order to reproduce the line ratios based Oph much shorter than Dobrzycka et al.’s value, so that on their shock models of nova shells. Smith et al. (1996) we adopt γ = 0.25 for d = 0.6 kpc in this paper. Even estimated metal abundance of AG Dra as [Fe/H]= 1.3 if we use γ = 0.4, the light curve itself is not so differ- and suggested that all yellow symbiotic stars including− RS ent from that for γ = 0.25 as shown in Hachisu & Kato Oph are metal-poor. They attributed K6 I–II to the cool (2000a). Only the difference is in the irradiation effect of component of RS Oph. the M-giant, but the irradiation itself is relatively small If the envelope of the white dwarf is also metal-poor, the compared with the case of T CrB (γ = 1). mass determination of the white dwarf from the light curve Recent studies on the radii of M-giants in symbiotic stars indicate γ 0.5 for the normal S-type symbiotic stars fitting may not be correct because the optically thick wind ∼ is driven by iron lines and the decay timescale is much (e.g., M¨urset & Schmid 1999). The cool component of longer in low metallicity envelopes. Therefore, we must the RS Oph system is not a normal giant. It is suggested reexamine the light curves of the 1985 outburst by intro- to be metal-poor and sometimes classified to K-giant (e.g., ducing lower metallicities and study the possible range of Smith et al. 1996) rather than M-giant. Therefore, we the metallicity and the WD mass. may adopt the relatively small radius of γ =0.25. As for the orbit of the binary system, the ephemeris for the inferior conjunction of the red giant in front is 12 Hachisu and Kato

2,444,999.9+460 E (Dobrzycka & Kenyon 1994). The 6.2. Development of the 1985 Outburst × separation is a = 325.0R⊙, the effective radii of the inner In this subsection, we summarize the observational de- critical Roche lobes for the WD component and the RG ∗ ∗ velopments of the RS Oph 1985 outburst and discuss them component are R1 = 138.6R⊙ and R2 = 108.2R⊙. For ∗ in relation to our model calculation. γ =0.25, we have RRG =0.25R2 27R⊙. Rosino (1987) and Rosino & Iijima (1987) divided the We assume 50% efficiency of the∼ irradiation of the red outburst into several phases along with the development giant (RG) component (η = 0.5). The nonirradiated RG of the optical spectrum: [[Phase 0—2 days]] after the opti- photospheric temperature of the RG is a parameter which cal maximum, broad emission lines of H I, He I, and Fe II is determined to be fitted with observation. Here, we adopt are observed with narrow absorption component (40—60 T = 3400 K against γ =0.25 as shown in Figure 11. − RG km s 1) and broad P Cygni troughs ( 2700 to 3700 km We have also included the contribution of irradiated ac- −1 II − ˚ −II I cretion disk (ACDK) luminosity. Here, we also assume s ). [[Phase 2—14 days]], He 4686 A , [N ], [O ], [O III], [Ne III], [S II], and [Fe II] 6374 A˚ appear. [[Phase 50% efficiency of the ACDK irradiation (ηDK =0.5). The viscous heating is neglected because it is much smaller 14—29 days]], He II 4686 A˚ and N III 4640 A˚ become than that of the irradiation effects. The temperature of strong. [[Phase 29—51 days]], coronal lines appear and be- the unheated surface of the ACDK including the rim is come strong. [[Phase 51—72 days]], the coronal line phase assumed to be Tdisk = 2000 K. We have checked two other lasts. [[Phase 72—119 days]], the coronal lines disappear. II ˚ cases of Tdisk = 1000 and 0 K and have found no significant Most of high excitation lines, including He 4686 A, are differences in the light curves. The power of ν is assumed fading. to be ν = 9/8 from the standard disk model. We have Snijders (1987a) presented an overall development of checked the dependency of the light curves on the param- the UV light curve during the 1985 outburst of RS Oph as eter ν by changing from ν =9/8 to ν = 2 but cannot find already shown in Figures 11 and 12. any significant differences as far as the disk rim is small Since RS Oph exploded in a circumstellar envelope (β =0.01—0.05). that is formed by the slow wind of the red giant com- The late-phase V light curve (t 5—100 days after ponent (Wallerstein 1958), the rapidly expanding ejecta − maximum) indicates strong irradiation∼ of the ACDK. We ( 3, 000—4,000 km s 1, e.g., Rosino 1987; Snijders ∼ have examined the case of a lobe-filling RG companion 1987a) of the outburst is decelerated by the slowly expand- − (γ = 1), but the irradiation of the RG is too luminous ing (. 20 km s 1) circumstellar matter to form a strong to be compatible with the observational light curves. On shock as noticed by Pottasch (1967). The secular varia- the other hand, a smaller size of the RG such as γ =0.25 tion of coronal lines was interpreted by Gorbatskii (1972) ( 27R⊙) gives a reasonable fit with the observations as in terms of a decelerating shock wave (see also Wallerstein shown∼ in Figures 11 and 12. & Garnavich 1986, for the 1985 outburst). This shock- To fit the late-phase light curve, we have calculated the heated layer is also responsible for soft X-ray (Mason et two cases of the disk parameters, i.e., (α, β) = (0.1, 0.05) al. 1987) and radio (Hjellming et al. 1986) emission as and (0.008, 0.05), as shown in Figure 11. Here, two solid will be mentioned below. lines denote the cases of α = 0.1 (Rdisk = 14 R⊙) and In our numerical modeling, the photospheric temper- β = 0.05 and α = 0.008 (Rdisk = 1 R⊙) and β = 0.05, ature of the WD envelope has once decreased down to respectively. Something between these two can roughly Tph = 12, 200 K at the optical maximum (HJD 2,446,093). reproduce the light curve of the 1985 outburst but the Then, it increases gradually to [[Phase 0—2]] Tph = 12, 200 gradual transition from α =0.1 to α =0.008 is much bet- — 17,700 K for the WD photospheric radius of Rph = 45 ter for reproducing the observations as shown in Figure — 21 R⊙, [[Phase 2—14]] Tph = 17, 700 — 67,100 K 12. It is very likely that the surface of the accretion disk for Rph = 21 — 1.6 R⊙, [[Phase 14—29]] Tph = 67, 100 has been dragged by the strong wind and been gradually — 114,000 K for Rph = 1.6 — 0.56 R⊙, [[Phase 29— blown in the strong wind. As a result, the accretion disk 51]] Tph = 114, 000 — 181,000 K for Rph = 0.56 — has become smaller during the strong wind phase ( 70 0.23 R⊙, [[Phase 51—72]] Tph = 181, 000 — 302,000K for ∼ days after maximum). Thus, we adopt Rph =0.23 — 0.083 R⊙, [[Phase 72—119]] Tph = 302, 000 — 1,020,000 K for Rph =0.083 — 0.0037 R⊙. The accre- (t−t0)/72 α1 tion disk reappears 4 days after the optical maximum, i.e., α = α0 , for t0

[Ne III], [S II], and [Fe II] 6374 A˚ in [[Phase 2—14]] and yr−1 and this causes the decrease in the electron density 2 [[Phase 14—29]] is consistent with our model because the because of ρ M˙ wind/4πr vwind. (2) After the hot wind WD photospheric temperature increases to high enough to from the WD∼ stops, the cool red giant wind resumes to excite these elements. The appearance of the coronal lines expand slowly and the newly ionized red giant wind dom- are related with the appearance of X-rays heated by the inates the spectrum. blast shock as shown below (Mason et al. 1987) because RS Oph became a strong X-ray source with characteris- they need high temperature photons of T & 1 106 K as in tic temperature of a few times 106 K (Mason et al. 1987). × the solar corona. The disappearance of the coronal lines Mason et al. observed X-rays from 55 days after the optical is usually understood by the shock breakout at the day maximum through 251 days with the medium-energy pro- 70. However, the quick disappearance of the coronal portional counter array (ME) and the low-energy imaging ∼ lines in [[Phase 72—119]] might be caused by the extinc- telescope (LE) onboard EXOSAT. They made six separate tion of the hydrogen shell burning, as suggested by Shore observations, 55, 62, 74, 83, 93, and 251 days after the op- et al. (1996), because it is consistent with the extinction tical maximum. The source became very weak in the ME time of the steady hydrogen shell burning in our model. from day 83 and only upper limits were derived. How- Evans et al. (1988) presented infrared observation (J, ever, there was still some residual X-ray emission detected H, K, and L bands) of the RS Oph 1985 outburst. The even more than 250 days after maximum only with the decline rates of J, H, K, and L are almost the same as LE band. The authors interpreted the early X-ray emis- t2 30 days (t2 is the time taken to drop 2 mag from max- sion in terms of emissions by the pre-outburst M-giant’s ∼ imum). At each wave length it is noticeable that the IR wind that had been shock-heated by the outburst ejecta flux of RS Oph had returned to its pre-outburst value by (see also Bode & Kahn 1985; O’Brien & Kahn 1987; Itoh the day 85 (after the optical maximum, but Evans et al. & Hachisu 1990; O’Brien, Bode, & Kahn 1992). On the ∼ took the optical maximum on HJD 2,466,095.5 as a origin other hand, the late X-ray emissions have been interpreted of time) and that the decline continued beyond this level. in terms of the continued hydrogen shell burning on top However, the IR fluxes had reverted to their pre-outburst of the white dwarf. This interpretation is not consistent levels until the day 400. The post-outburst minimum with our model calculation because the steady hydrogen ∼ lasted during about 300 days, i.e., from the day 100 to shell-burning ended long before the day 251, i.e., at the the day 400. This post-outburst minimum is discussed∼ day 110. We will discuss this point again in the next ∼ below in relation to the super-soft X-ray detection on the subsection∼ in relation to the distance to RS Oph. day 251 and the distance to RS Oph. Evans et al. (1988) also estimated the IR flux from the 6.3. The Distance to RS Ophiuchi circumstellar material in the RS Oph system by removing the contribution from the M-giant companion. Assuming The distance to RS Oph is determined to be 0.57 kpc that the He II 1640 A˚ line and the IR continuum do in from the dereddened UV flux (Snijders 1987a) as shown fact arise in the same region, they fitted the IR continuum in Figures 11 and 12. We assume the accretion disk exists on HJD 2,466,119 by a nebular continuum having electron during the outburst. We cannot reproduce the UV light 5 temperature Te 1.6 10 K. In our model calculation, curve if the disk does not exist or its size is smaller than we have the photospheric∼ × temperature of the bloated WD α < 0.003, i.e., . 0.5R⊙, as already shown in Hachisu & 5 envelope Tph 1.1 10 K on the same day, being roughly Kato (2000a). The UV light curves depend hardly on the consistent with∼ Evans× et al.’s estimation, though Evans et disk parameters of α, β, or ν when α & 0.01, mainly be- al. attributed it to the cooled shocked layer. Evans et cause the UV light is coming from the innermost part of al. also estimated that the spectral type of RS Oph cool the accretion disk. The distance depends weakly on the component is K8 III ( 2 sub-divisions) from [2.35] [2.2] irradiation efficiency (ηDK), for example, d =0.70 kpc for color. ± − ηDK = 1.0 or d = 0.50 kpc for ηDK = 0.25. The red gi- Shore et al. (1996) also divided the outburst into two ant does not contribute to the UV because its irradiated temperature is as low as 4000 K. Thus, the distance phases along with the development of the UV spectrum: ∼ The first phase, [[Phase 0—60 days]] after the optical max- estimation is rather robust. We have calculated the absorption of A = 11.09 imum, is the broad line phase. The high velocity material V − dominated the emission. This is roughly consistent with 5 log(570/10) = 2.3 for d = 0.57 kpc from our light curve the strong wind phase of our numerical modeling, i.e., until fitting. This absorption of AV = 2.3 is consistent with HJD 2,446,170. In the second phase, i.e., [[Phase 60—130 the color excess of E(B V )=0.73 (Snijders 1987b) be- cause of A =3.1 E(B −V ). This distance of 0.6 kpc is days]], the newly ionized red giant wind dominated the V − ∼ spectrum. This period is also consistent with the steady not consistent with other estimations of d =1.6 kpc from hydrogen burning phase of our numerical modeling, i.e., the hydrogen column density (Hjellming et al. 1986), until HJD 2,446,202. d = 1290 pc (Harrison et al. 1993), and d = 1 kpc Shore et al. also discussed that the line ratio of Si III] (Sekiguchi et al. 1990), both of which are determined 1892A˚ toC III] 1910 A,˚ which measures the electron from the K-band luminosity. In the following, therefore, density (Nussbaumer & Stencel 1989), steadily decreased we discuss extensively other distance estimations support- from an initial value of 1.0 on HJD 2,446,106 to a mini- ing our short distance of 0.6 kpc to RS Oph. mum of 0.3 on HJD 2,446,164. Then, the ratio began to The VLBI map 77 days after the optical maximum (Tay- lor et al. 1989) shows an elongated structure with a long rise thereafter, reaching 0.6 by HJD 2,446,202. This trend ′′ ′′ can be explained by our wind model: (1) The wind mass axis of 0. 22 and a short axis of 0. 08. If we assume the −4 −7 loss rate is gradually decreasing from 10 to 10 M⊙ distance of d = 1.6 kpc to RS Oph, this fastest aver- ∼ ∼ age expansion velocity should be 4,000 km s−1. This is 14 Hachisu and Kato not consistent with the initial velocity of ejecta 3, 000— during the post-outburst minimum 100—400 days after 4, 000 km s−1 because the ejecta is strongly decelerated∼ by maximum pointed by Evans et al.∼ 1988). If our value −7 −1 the circumbinary matter that is formed by the red giant of M˙ acc = 1.2 10 M⊙ yr is adopted as the mass wind in quiescence. If the distance is d =0.6 kpc, on the accretion rate before× the 1985 outburst, these line fluxes other hand, the average expansion velocity is much slower before/after the outburst are roughly consistent with the − ( 1, 400 km s 1), which is consistent with the initial ve- ratio of the estimated mass accretion rates before and af- ∼ locity of the ejecta. Then, the slowest average expansion ter the outburst (1.2/0.2) if the line fluxes depend pro- − velocity is estimated to be 500 km s 1, which is also portionally on the gravitational energy release. The day ∼ −1 consistent with the optical line width (e.g., 490 km s 251 was just in the period of the post-outburst minimum ∼ 110 days after optical maximum, Snijders 1987a). (m 11.5) of the day 100—400, after which the visual V ∼ ∼ The outer edge of the circumbinary matter (CBM) that luminosity increased to mV 10.5, by about one magni- 15 is formed by the red giant wind is about RCBM . 1 10 tude (e.g., Fig. 1 of Evans et∼ al. 1988). −1 × cm, since an upper limit of 20 km s to the velocity of Bohigas et al. (1989) estimated the Hβ intensity pho- the red giant wind has been given by Gorbatskii (1972) toionized by X-ray photons from the central hot WD and based on spectra by Dufay et al. (1964). Here, we assume compared it with the observation 201 days after maximum. that the ejecta of the previous 1967 outburst swept up the Their theoretically estimated value of previously existing CBM far away. The soft X-ray rapidly I(Hβ)=1.24 10−10ergs cm−1s−1 fell down from 60—70 days after the optical maximum so × − R 2 d 2 that the shock had broken out of the circumbinary matter WD (60) on that day (Mason et al. 1987; Itoh & Hachisu 1990; ×  1 109 cm 1.6 kpc × Lloyd et al. 1993). If the outer edge of the circumbinary is 3.5 times larger than the observed value of I(Hβ) = matter corresponds to the elongated radius of the VLBI 3.55 10−11 ergs cm−1 s−1 when we adopt R = 109 cm map on the day 77 (0′′. 11 of half axis), we obtain the dis- × WD 21 and d = 1.6 kpc. However, if we use RWD =0.003 R⊙ and tance to RS Oph, d . R /(∆θ)=2 10 cm = 0.65 −11 CBM × d = 0.6 kpc, the estimated value is I(Hβ)=3.8 10 kpc. ergs cm−1 s−1, being consistent with the observed× inten- Very soft X-rays were observed 251 days after the optical sity. maximum in the 1985 outburst of RS Oph (Mason et al. X-rays were also observed at the quiescent phase (Orio 1987). Mason et al. estimated a blackbody temperature 1993), but the flux is too low to be compatible with the of 3.5 105 K and a total energy flux of × TNR model. One possible explanation is an absorption by d 2 the massive cool wind (e.g., Shore et al. 1996) from the L =1 1037 ergs s−1, (57) X,obs × 1.6 kpc RG component as suggested by Anupama & Mikolajewska (1999). On the other hand, soft X-rays were observable 21 −1 for a hydrogen column density of NH =3 10 cm , and just after the outburst, because the ejecta swept away the × concluded that these soft X-rays come from a WD pho- cool wind around the RS Oph system. The cool red giant tosphere with steady hydrogen shell burning. However, wind fills again its Roche lobe and can cover the hot com- −1 steady hydrogen shell burning has stopped 109 days af- ponent in a time scale of ∆t = a/v 300R⊙/10 km s = ter the optical maximum (HJD 2,446,202) in our model of 300 days, which is roughly consistent∼ with the recovery of X =0.70, Z =0.004, and MWD =1.377M⊙ (Figs. 11 and the quiescent V luminosity as shown in Figure 1 of Evans 12), that is, this soft X-ray detection is indicating accretion et al. (1988). luminosity rather than hydrogen shell burning luminosity. Hjellming et al. (1986) estimated the distance to 5 Using a blackbody temperature of TBB 3.5 10 K and RS Oph of d = 1.6 kpc from H I absorption line mea- ∼ × a WD radius of 0.003 R⊙ for 1.377 M⊙, we estimate a surements, using the hydrogen column density of NH = 21 −2 total luminosity of (2.4 0.6) 10 cm and the relation of NH/(Ts d) = 2 4 1.59±1019 cm× −2 K−1 kpc−1 together with T = 100K.Ifa LX,WD =4πR σT 120L⊙. (58) s WD BB ∼ large× part of this hydrogen column density stems not from Thus, a rather short distance of 0.4 kpc to RS Oph the Galactic absorption but from the local one belonging can be derived by equating equation (57) and (58), i.e., to RS Oph, the distance is overestimated. For example, LX,obs = LX,WD. Thus, our short distance of 0.6 kpc is in the direction of the recurrent nova U Sco, Kahabka et roughly consistent with the soft X-ray observation by Ma- al. (1999) reported a hydrogen column density of (3— son et al. (1987) about eight months after the optical 4) 1021 cm−2, which is much higher than the Galactic maximum. This low luminosity corresponds to the accre- hydrogen× column density of 1.4 1021 cm−2. If RS Oph tion luminosity of is in the same case as in U Sco,× the distance of 1.6 kpc is 2LX,WDRWD −7 −1 just an upper limit and we possibly derive a much shorter M˙ acc = 0.2 10 M⊙ yr . (59) GMWD ∼ × distance. This accretion rate is six times lower than our estimated −7 −1 6.4. The Ejected Envelope Mass and Possibility of SN Ia mass accretion rate of M˙ acc =1.2 10 M⊙ yr averaged between 1967 and 1985. × Explosion Dobrzycka et al. (1996) pointed out that the line The envelope mass at the optical maximum is estimated −6 fluxes in H I and He I before the 1985 outburst are to be ∆M = 2.2 10 M⊙, which indicates a mass ac- × −7 −1 about four times larger than the fluxes after the 1985 cretion rate of 1.2 10 M⊙ yr during the quiescent outburst, thus indicating a decrease in the mass accre- phase between the× 1968 and 1985 outbursts. In our cal- −6 tion rate after the outburst (see also the low luminosity culation, about 90% of the envelope mass (2.0 10 M⊙) × Recurrent Novae 15 has been blown off in the optically thick wind, and the condition. Thus, we may conclude that RS Oph is an −6 residual 10% (0.2 10 M⊙) has been left and added immediate progenitor of an SN Ia even if the donor is a to the helium layer× of the WD. The residual mass itself metal-poor star. depends on both the hydrogen content X and the WD mass. It is roughly ranging from 10% (X = 0.70) to 6.5. Absence of Secondary Maximum 15%(X = 0.50) including ambiguities of the WD mass Among the four recurrent novae with a red giant com- (MWD = 1.377 0.01M⊙). Therefore, the net mass in- panion, only T CrB shows a secondary maximum of ∆V ± −8 −1 creasing rate of the WD is M˙ =1.2 10 M⊙ yr . ∼ He × 2 mag brighter than the quiescent level 120—190 days after Our calculated ejecta (wind) mass seems to be inconsis- the optical maximum. We discuss in more detail on the tent with observation. Bohigas et al. (1989) estimated difference between T CrB and the others. The essential the (fully ionized hydrogen) shell mass from their obser- difference exists in the fact that the red giant component vation 201 days after maximum. The obtained value is in T CrB fills its inner critical Roche lobe but the others −5 2 M 1.5 10 (d/1 kpc) M⊙. Then, they obtained shell ∼ × do not (see the following sections for V745 Sco and V3890 the ejecta mass form their equation (8), i.e., Mejecta Sgr). 2 ∼ 3(Vs/Ve) Mshell, where Vs is the shell velocity after the In the T CrB system, the mass transfer continues from shock has broken out of the circumstellar envelope and the red giant component to the white dwarf even in the Ve is the velocity of the ejecta before it collides with the outburst (hot) wind phase. Moreover, the strong irra- −1 circumstellar matter. Assuming Vs = 200 km s and diation increases the mass transfer rate from the red gi- −1 Ve = 2000 km s , Bohigas et al. estimated Mejecta ant component; thus increasing the radius of the optically −7 2 ∼ 4.5 10 (d/1 kpc) M⊙. This ejecta mass is not consis- thick region of the accretion disk. The unstable warping × −6 tent with our numerical model (∆M =2.0 10 M⊙) wind × grows shortly after the hot wind stops because the mass if we use d =0.6 kpc. transfer rate satisfies the warping condition (49). This inconsistency comes from the inappropriate pair On the other hand, the RS Oph system cannot main- of velocities because the estimated mass depends sensi- tain the mass transfer because the cool red giant wind has tively on the ratio Vs/Ve. If we adopt an appropriate pair once been evacuated by the outburst (hot) wind/ejecta of velocities, our wind model is consistent with the ob- from the white dwarf. During the day 0—80 (hot wind servational shell mass. Starting from the set of d = 0.6 phase), the growth of warping was suppressed by the hot −6 kpc and Mejecta ∆Mwind 2.0 10 M⊙, we obtain wind because the hot wind (matter) momentum is much ≡−6 ∼ × Mshell 5.4 10 M⊙ and Vs/Ve 0.35. Adopting larger than the photon momentum as shown by Kato & ∼ × −1 ∼ Ve 4000 km s from the observation (e.g., Rosino & Hachisu (1994). During the day 80—120 (steady burning ∼ −1 Iijima 1987; Snijders 1987a,b) gives Vs 1400 km s . phase), the warping instability may grow to some extent. −1 ∼ This velocity of 1400 km s is consistent with the ra- However, the disk cannot contribute much to the V light, dio map obtained 77 days after maximum (Taylor et al. since the optically thick region of the accretion disk has 1989). The VLBI map shows an elongated structure of an- become as small as or smaller than 1R⊙. gular dimensions of 0′′. 22 (long axis) and 0′′. 08 (short ∼ ∼ ∼ The mass transfer to the white dwarf resumes 400 axis), which is consistent with the distance of 0.6 kpc and days after the optical maximum, that is, after the hot∼ wind the fastest expansion velocity of 1400 km s−1 after shock, ′′ − −1 ′′ stops ( 80 days) and then the cool red giant wind fills i.e., 0. 092 (d/1 kpc) 1 (V /1000 km s ) 0. 22. Thus, ∼ −1 · s ∼ again the binary system (∆t = a/v 300 R⊙/10 km s we have obtained a consistent set of physical quantities, = 300 days) as already discussed earlier.∼ At that time, −6 −6 Mshell 5.4 10 M⊙, Mejecta 2.0 10 M⊙, and the luminosity of the white dwarf has once decreased as ∼ ×−1 ∼ × Vs = 1400 km s for the distance of d =0.6 kpc. low as 100—200 L⊙ and too low (two orders of mag- Thus, we obtain the cool wind mass loss rate of 4 nitude∼ lower) to excite the warping instability as easily −7 −1 × 10 M⊙ yr (= Mshell/18 yr), indicating that one fourth understood from the condition for RS Oph (Southwell et of the cool wind has been captured by the WD. Obser- al. 1997), i.e., vationally, a very high mass loss rate of the cool wind, −5 −1 ˙ 1/2 10 M⊙ yr , was suggested by Shore et al. (1996) Macc Rdisk Lbol ∼ − − . − based on the UV data. However, Dobrzycka et al. (1996) 3 10 8 M⊙ yr 1  R⊙  2 1038 ergs s 1  × × − argued against such a high mass loss rate of the M-giant R 1/2 M 1/2 cool wind. If the wind mass loss rate approaches such WD WD (62). −5 −1 × 0.003 R⊙  1.377 M⊙  a high rate of 10 M⊙ yr , a very large 10µm excess should be observed. But, the giant of RS Oph has only In the cool wind accretion phase, our estimated accretion −7 −1 a modest 10—20µm excess. They estimated a few times rate of the WD in RS Oph is M˙ acc 1.2 10 M⊙ yr −7 −1 36 −1 ∼ × 10 M⊙ yr (e.g., Kenyon, Fernandez-Castro, & Sten- (therefore, Lbol 5 10 ergs s for the accretion lumi- cel 1988). This mass loss rate of the cool wind is very nosity), which also∼ does× not meet the above condition even consistent with the present estimation. if Rdisk = 10 R⊙, so that the radiation-induced instability Finally, we examine whether or not RS Oph will explode does not occur in RS Oph. This is the reason why RS Oph as a Type Ia supernova in quite a near future. One of the does not show a second peak. conditions for SN Ia explosions is changed as follows: The other two, V745 Sco and V3890 Sgr, are similar −8 −1 M˙ > M˙ − 4 10 M⊙ yr , (61) systems to RS Oph, so that a secondary maximum was acc He det ∼ × for low metallicity stars (Nomoto 1982; Nomoto & Kondo not observed, as will be shown in the following sections. 1991) instead of equation (52). Our estimated mass accre- −7 −1 6.6. Summary of RS Oph Outburst tion rate of M˙ = 1.2 10 M⊙ yr meets the above acc × 16 Hachisu and Kato

The main features of the RS Oph outbursts can be recurrent novae with a red giant secondary. well understood with a thermonuclear runaway model on a very massive white dwarf close to the Chandrasekhar 7.1. The Model of the 1989 Outburst mass limit. We have summarized the following important We have calculated V light curves for five cases of the points: WD mass (MWD = 1.3, 1.35, 1.36, 1.37, and 1.377M⊙) (1) It has been suggested that the M-giant in RS Oph by assuming the hydrogen content X =0.70, 0.50, or 0.35 is metal-poor. Assuming an upper limit of white dwarf with the solar metallicity (Z = 0.02). By fitting light mass, i.e., MWD =1.377 M⊙, we have obtain a lower limit curves in the first 25 days of the 1989 outburst, we have de- of metallicity, i.e., Z 0.004 from the light curve fitting of termined the WD mass of M =1.35 0.01M⊙ as shown ≥ WD the early phase ( 0—4 days). In other words, the early in Figure 13. Here, we have shifted the± visual magnitude ∼ visual light curve can be well reproduced only by a bloated obtained with orange filter or 2415 films (Liller 1989a,b) white dwarf photosphere of MWD = 1.377 0.01 M⊙ by 0.8 mag down (crosses, ) because they are sys- ± for a low metallicity of Z = 0.004. On the other hand, tematically∼ about 0.8 mag brighter× than the V -magnitude MWD =1.35 0.01 M⊙ is obtained for the solar metallic- obtained by Sekiguchi et al. (1990). The decline rate of ± ity of Z =0.02. the V light in the early phase depends very sensitively on (2) The ensuing phase ( 5—80 days) visual light curve the WD mass, but depends hardly on the hydrogen con- ∼ cannot be reproduced only by a white dwarf photosphere. tent X, the companion (its mass and radius), or the size It requires an irradiation effect of the accretion disk, the of the accretion disk, so that we can safely determine the size of which is gradually decreasing from 10 R⊙ to WD mass from the early light curve (see also Kato 1999). ∼ 1 R⊙. Our fitting indicates an apparent distance modulus of ∼ (3) The UV light curve can also be well reproduced with (m M)V = 17.0. Thus, the distance to V745 Sco is the same model as that of the visual luminosity. The early derived− to be d = 5.0 kpc if we adopt the absorption of phase UV flux mainly comes from the bloated white dwarf AV = 3.1 E(B V )=3.5, where E(B V )=1.1 (Har- photosphere while the late phase UV flux comes from the rison et al.× 1993;− Williams 1994). − irradiated disk surface because the main emitting region From our solutions, the envelope mass at the maximum −6 of the white dwarf photosphere has shift from UV to soft V light is estimated to be ∆M =4.6 10 M⊙, indicat- X-ray in the late stage. × −7 −1 ing a mass accretion rate of M˙ acc = 0.9 10 M⊙ yr (4) The optically thick wind stops on the day 80. The × ∼ between the 1937 and 1989 outbursts. About 95% of the WD photosphere quickly shrinks. The sharp drop in the envelope mass has been blown off by the wind and the −7 UV flux around the day 80 can be understood from this residual 5% (2.3 10 M⊙) has been left and added to quick shrink of the WD photosphere. the helium layer of× the WD. Therefore, we obtain the net (5) The distance to RS Oph is estimated to be 0.6 kpc −8 mass increasing rate of the WD as M˙ He =0.45 10 M⊙ from the fitting with the dereddened UV flux. This dis- yr−1. × tance is consistent with the V light curve fitting. (6) The envelope mass of the bloated white dwarf is 7.2. Irradiation Effect of Companion M-giant −6 2.2 10 M⊙, suggesting the mass accretion rate of × −7 −1 In the late phase (t 25—200 days after the optical 1.2 10 M⊙ yr . maximum), contributions∼ of the irradiated red giant (RG) (7)× The wind carries away about 90% of the envelope mass, and accretion disk (ACDK) to the V light may become so that the residual 10% can be accumulated to the helium dominant like as in T CrB and RS Oph. First, in this sub- layer of the white dwarf. As a result, the white dwarf grows −8 −1 section, we examine the irradiation effect of the M-giant, at a rate of 1.2 10 M⊙ yr . and then, in the next subsection, we also study the irradi- (8) RS Oph will× certainly explode as a Type Ia supernova ation of the ACDK. in quite a near future if the white dwarf consists of carbon In order to examine companion’s irradiation effect, we and oxygen. need to know the separation, mass ratio, inclination an- gle, and the radius of the red giant, but none of them is 7. v745 scorpii known. The lack of both the orbital information and the V745 Scorpii underwent the second recorded outburst visual luminosity observation in the late phase clouds the in 1989 (the first recorded outburst in 1937), and has been accurate determination of the outburst model because of recognized as a member of recurrent novae (e.g., Sekiguchi ample ambiguity in the system parameters. et al. 1990, and references therein). Its light curve is char- The orbit is assumed to be circular. It is observationally acterized by a very rapid rise (in about half a day) and suggested that orbital periods of S-type symbiotic stars are related to their spectral type by decline (t3 9 days, t3 is the time taken to drop 3 mag from maximum).∼ The general spectral evolution closely P (1.31)S 190 days, (63) ∼ × resembles that of RS Oph (e.g., Sekiguchi et al. 1990). where the power of S means the sub-spectral type of M- Although its orbital period, the most important system giants (S = 4 for M4 III, see, e.g., M¨urset & Schmid 1999; parameter, is still not known, the secondary cool compo- Harries & Howarth 2000). From this relation, we assume nent has been classified as M6—8 III (from optical spectra, P = 560 days for V745 Sco (M4 III, e.g., Harrison et al. Sekiguchi et al. 1990; Duerbeck & Seitter 1990; Williams 1993) and the ephemeris for inferior conjunction of the M- et al. 1991), M4 III (from the infrared spectra and colors, giant in front, φmin = HJD 2, 448, 700.0+ 560E. In what Harrison et al. 1993), or M4—6 III (from the absorption follows, we refer this ephemeris simply as φmin = 700 day. line indices, Anupama & Mikolajewska 1999). Therefore, It has also been suggested that M-giants in S-type sym- V745 Sco is a twin system to RS Oph, a member of the biotic stars lies well within the inner critical Roche lobe, Recurrent Novae 17 that is, γ 0.3 0.6 (e.g., M¨urset & Schmid 1999). surface of the ACDK including the rim is assumed to ∼ − Therefore, we assume roughly the ratio of be Tdisk = 1000 K. We have examined another case of R 1 Tdisk = 0 K, but found no significant differences in the V RG γ . (64) R∗ ≡ ∼ 2 light curves. The power of ν is assumed to be ν = 9/8 2 from the standard disk model. In addition, we also assume 50% efficiency of the irradia- We have finally obtain allowable sets of the disk param- tion (ηRG = 0.5). For the companion M-giant, Sekiguchi eters of α . 0.01 and β . 0.01 for the orbital inclination et al. (1990) suggested mV = 17 in quiescence from angle of i 80◦ as shown in Figure 13. Here, the thin the V magnitude just after the outburst. Anupama & solid line denotes∼ the case of α =0.01 and β =0.01, while Mikolajewska (1999) obtained m5500 = 17.7 in March the dash-dotted line corresponds to α = 0 (no accretion 1998. We adopt here mV = 17.7 in quiescence. disk). Both can reproduce the V light curve of the 1989 The light curves are calculated for six cases of the com- outburst. panion mass, i.e., MRG = 0.5, 0.6, 0.7, 0.8, 1.0, and We have examined six cases of the inclination angle of 1.2M⊙. Since we have obtained similar light curves for the orbit, i = 80◦, 70◦, 60◦, 50◦, 40◦, and 30◦ and shown all of these six masses, we show here only the results three of them in Figure 14. The V light curves are more for MRG = 1.0M⊙. For this case, the separation is luminous for a lower inclination angle of the orbit, since a a = 380.1R⊙, the effective radii of the inner critical Roche ∗ larger area of the ACDK can be seen from the Earth. This lobes are R1 = 154.0R⊙ for the WD component, and ∗ means that lower inclinations of the orbit contradicts the R2 = 134.3R⊙ for the RG component. If γ =0.5, then we ∗ existence of a large ACDK. To summarize, for any case, have R = R = 0.5R 67R⊙. Then, the nonirradi- 2 RG 2 ∼ the ACDK should not contribute so much to the V light. ated photospheric temperature of the RG is estimated to Figure 15 shows the dependency of the hydrogen content be Tph,RG = 2700 K from the light curve fitting in quies- for the case with an accretion disk. The V light curves sug- cence. gest that the hydrogen content is not clearly identified by We have finally adopted the best fit set of the inclina- ◦ the observation but it may be somewhere between X =0.5 tion angle i = 80 for the ephemeris of φmin = 700 day and 0.7 (solid curves). as shown in Figure 13. Dash-dotted curves in Figure 14 Thus, a very small size of the optically thick ACDK, i.e., show V light curves for various inclination angles (i = 30◦, ◦ ◦ Rdisk . 1 R⊙, in the late phase indicates no secondary 60 , 80 ) of the orbit for the ephemeris of φmin = 700 day. maximum of the light curve as discussed in RS Oph even Their light curves depends hardly on the inclination angle if the warping instability sets in. because the non-irradiated hemisphere of the M-giant faces towards the Earth. On the other hand, V light is much 7.4. The Distance, Envelope Mass of the White Dwarf, brighter (dotted curves) if the ephemeris of φmin = 500 day is adopted because the irradiated hemisphere of the and Possibility of SN Ia Explosion M-giant faces towards the Earth, which is not compatible The distance to V745 Sco is estimated to be d = 5.0 with the observation. kpc with the absorption of AV =3.5. This is roughly con- Figure 15 shows the dependency of the hydrogen con- sistent with Harrison et al.’s (1993) estimation of d =4.6 tent, X. The early-phase V light curve depends hardly on kpc. Harrison et al. attributed M4 III to the cool com- the hydrogen content, while the late-phase V light curve is ponent from the infrared CO lines. Assuming the abso- weakly identified by the hydrogen content of the WD en- lute magnitude of MK = 5.5 (M4 III), they obtain the velope. The optically thick, outburst wind stops 59, 69, 82 distance modulus of (m −M) = m M A = − 0 K − K − K days after maximum, and the steady hydrogen shell burn- 8.2 ( 5.5) 0.4= 13.3, where they adopted mK = 8.2 ing ends 71, 90, 123 days after maximum, for X = 0.35, from− Sekiguchi− − et al. (1990) and used a relation between 0.50, and 0.70, respectively. The V light curve indicates the absorptions of AK and AV given by Rieke & Lebofsky that the hydrogen content is somewhere between X =0.5 (1985), i.e., AK /AV =0.112. and X =0.7 (dash-dotted curves). Sekiguchi et al.’s (1990) K-band luminosity mK =8.21 We have also calculated light curves for a much longer (87 days after maximum) is about 0.2 mag brighter than orbital period (P = 1000 days, for M6 III, e.g., Anupama Harrison et al.’s (1993) K-band luminosity mK = 8.41 & Mikolajewska 1999), and have obtained similar results (two years after the 1989 outburst). This is consistent to those for P = 560 days. In this case, the binary param- with our models as seen in Figure 13, because the cool ∗ ∗ eters are a = 559.4R⊙, R1 = 226.7R⊙, R2 = 197.6R⊙, component had been strongly irradiated by the WD photo- R2 = 99R⊙. We have estimated TRG = 2500 K for the sphere always during Sekiguchi et al.’s observations. Using best fit model in the late phase. m =8.41, we obtain a distance modulus of (m M) = K − 0 mK MK AK =8.41 ( 5.5) 0.4= 13.5, i.e., d =5.0 7.3. Small Irradiation Effect of Accretion Disk kpc,− being− coincided with− our− distance− estimation. In this subsection, we show that the ACDK should not At the end of this section, we discuss the possibility contribute so much to the V light. Since the visual light of SN Ia explosion. We have already obtained the average −7 −1 curve is well reproduced by M-giant’s irradiation with no mass accretion rate of M˙ acc =0.9 10 M⊙ yr between accretion disk, one may imagine that the accretion disk the 1937 and 1989 outbursts from× the envelope mass of −6 does not contribute much. Here, we determine the upper ∆M = 4.6 10 M⊙. This meets condition (52) for SN limit of the ACDK size in relation to the reason why a Ia explosions× (Nomoto & Kondo 1991), For all the three secondary maximum does not occur in V745 Sco. cases of X =0.70, 0.50, and 0.35, about 95% of the enve- We here again assume 50% efficiency of the ACDK irra- lope mass has been blown off in the outburst wind and the −7 diation (η = 0.5). The temperature of the unheated residual 5% (2.3 10 M⊙) has been left and added to the DK × 18 Hachisu and Kato helium layer of the WD. Therefore, we obtain the net mass condition of MWD & 1.35M⊙ is required from the consis- −8 −1 increasing rate of the WD of M˙ He =0.45 10 M⊙ yr . tencies with the TNR model, i.e., the mass accretion rate × −7 −1 Thus, we may conclude that V745 Sco is also an immedi- is smaller than M˙ acc . 1 10 M⊙ yr , and the rapid ate progenitor of Type Ia supernova if the hot component decline of the V light curve× 20—40 days after maximum as is a carbon-oxygen white dwarf. shown in Figure 16. The V light curve of a 1.34M⊙ WD model lies above the observational points on the day 26 8. v3890 sagittarii and the upper limits on the day 33 and 37, being not con- sistent with the observation. On the other hand, V light V3890 Sagittarii underwent the second recorded out- curve of a 1.35M⊙ or 1.37M⊙ WD is less bright during the burst in 1990 (since the 1962 first recorded outburst), and day 3—40. The lack of the luminosity during these days has been recognized as a member of the recurrent novae can be supplied by the irradiation of a large size accretion (Jones 1990a; Liller 1990). Its light curve is character- disk. ized by a very rapid rise (within two days) and fast decline Figure 17 shows V light curves of fixed-size disk models (t3 17 days; see, e.g., Fig. 1 in Gonz´alez-Riestra 1992). for the WD masses of 1.35 M⊙ (solid curves) and 1.37 M⊙ The∼ general spectral evolution closely resembles that of (dashed curves). The disk size parameters are α =0.7, 0.1, RS Oph (e.g., Sekiguchi 1990). The cool component has and 0.01. It is clear that the fixed-size disk model cannot been classified from optical spectra as M4 III (Sekiguchi reproduce the V light curve in the later stage. This is 1990), M8 III (Williams et al. 1991), or M5 III (Anupama simply because the large disk contributes too much in the & Mikolajewska 1999). The infrared spectra and colors later phase. indicate a spectral type M5 III (Harrison et al. 1993). Here, we summarize our parameters: the orbital pe- Thus, V3890 Sgr is a similar system to RS Oph and V745 riod of P = 730 days from equation (63) and the spectral Sco, a member of the recurrent novae with an M-giant type of M5 III (Anupama & Mikolajewska 1999) for the secondary. cool component, the ephemeris for the M-giant in front is φ = HJD 2,448,050.0 +730 E, γ =0.5 from equation 8.1. The Model of the 1990 Outburst min × (64), 50% efficiency of the irradiation (ηRG = 0.5), the We have calculated the total V light of the WD pho- M-giant mass of MRG =1.0 M⊙. In this case, the separa- tosphere and the non-irradiated RG photosphere for the tion is a = 453.5R⊙, the effective radii of the inner critical ∗ WD mass of MWD = 1.2, 1.3, 1.34, 1.35, 1.36, 1.37, and Roche lobes are R1 = 183.8R⊙ for the WD component ∗ 1.377M⊙ by assuming the hydrogen content X =0.70 and and R2 = 160.2R⊙ for the RG component. Then, we have ∗ the solar metallicity of Z = 0.02. Four of them (1.2, 1.3, R = 0.5R 80R⊙. The nonirradiated photospheric RG 2 ∼ 1.35, and 1.37 M⊙) are plotted in Figure 16 (solid curves). temperature of the RG is estimated to be Tph,RG = 2600 Here, we assume that the optical maximum was reached K from fitting with the quiescent V magnitude of mV 16 on HJD 2,2448,009. This is because the irradiation of the (Anupama & Mikolajewska 1999). For the accretion∼ disk, M-giant cannot contribute to the early phase V light. we assume 50% efficiency of the irradiation (ηDK = 0.5), The observed V light curve in the first 20 days can be the unheated surface temperature of Tdisk = 1000 K. roughly reproduced by the model of 1.3M⊙ WD, but it deviates greatly in the ensuing phase (20—40 days after 8.3. Evaporating Accretion Disk maximum). The envelope mass of the 1.3M⊙ WD at the The optically thick strong wind blows during the out- optical maximum is too large to be compatible with the burst so that the disk surface is certainly blown off in −6 TNR model. The envelope mass of ∆M 5 10 M⊙ the wind due to drag and ablation. Moreover, the mass −7 −1 ∼ × suggests M˙ acc 2 10 M⊙ yr , which is just below the transfer from the M-giant is certainly inhibited during the lower limit of steady∼ × hydrogen shell burning for X = 0.7, strong wind phase because the cool red giant wind has i.e., M˙ std in equation (7) and also in Figure 2. been evacuated by the hot wind from the WD as discussed Only a bloated WD photosphere cannot reproduce the V in RS Oph. In such a situation, it is very likely that the light curve even in the early phase of outburst. Moreover, edge of the disk is gradually dissipating with time. We our simple picture of the WD photosphere never shows a here simulate such an effect by shrinking the disk size ex- sharp spike on the light curve as seen in Figure 16. These ponentially, i.e., (t−t0)/30 features strongly suggest that the irradiation of the ac- α1 cretion disk is important even in the early phase: if the α = α0 , (65) α0  optically thick region of the accretion disk is as large as where we adopt α0 = 0.7 and α1 = 0.01, t0 is the time or larger than 10 R⊙, the decay of the V light becomes ∼ of the optical maximum, and the time t is in units of day. much slower as seen in RS Ohp (see, e.g., Fig. 11). There- We set the minimum size of the ACDK at α = 0.001 in fore, we seek a reasonable set of system parameters based our calculations. on the assumption that the V3890 Sgr system consists of Thus, we have obtained the best fit model for the WD the WD photosphere, the irradiated M-giant, and the ir- mass of MWD = 1.35 0.01 M⊙ as shown in Figures 17 radiated accretion disk. and 18. The envelope± mass at the maximum V light is −6 estimated to be ∆M = 3.1 10 M⊙, indicating a mass 8.2. Large Size of Accretion Disk in Early Phase × −7 −1 transfer rate of M˙ acc = 1.1 10 M⊙ yr between the In order to obtain a consistent set of system parame- 1962 and 1990 outbursts. The× efficiency of the ACDK irra- ters, we have to adopt the WD mass more massive than diation has been examined for other two cases of ηDK =1.0 MWD & 1.35M⊙, and a size of the irradiated ACDK larger (dot-dashed line) and 0.25 (dotted line) only for X =0.70 than Rdisk & 10 R⊙ at the initial phase of outburst. The and MWD =1.35 M⊙ as shown in Figure 18. Recurrent Novae 19

We have examined the dependency of the V light curve (Harrison et al. 1993), because the distance modulus on the orbital inclination angle. The luminosity of the is (m M)0 = (m M)V AV = 14.7 1.6 = 13.1. ACDK depends largely on the the orbital inclination an- This is− roughly consistent− with− Harrison et− al.’s (1993) gle mainly because a larger area of the ACDK can be seen estimation of d = 5.2 kpc. Harrison et al. determined from the Earth for a smaller inclination angle. Therefore, the infrared spectral type of the cool companion as M5 the inclination angle should be as low as i . 30◦ to re- III. Assuming M = 5.5 (M5 III), they obtained the K − produce the V light curve during t 5—20 days after distance modulus of (m M)0 = mK MK AK = maximum. ∼ 8.24 ( 5.5) 0.17 =− 13.6 by using− Rieke &− Lebof- − − − We have also calculated V -light curves for other four sky (1985) relation of AK /AV = 0.112. It should be cases of the companion mass, i.e., MRG =0.6, 0.7, 0.8, noted here that there are a large divergence in the ob- and 1.2 M⊙, but obtained similar light curves to MRG = served color excesses: Gonz´alez-Riestra (1992) deter- 1.0 M⊙. We have checked the other case of Tdisk = 0 mined E(B V )=1.1 from UV lines; Williams (1994) K, but found no significant differences in the light curves. obtained E(−B V )=0.76 from H I and He II line ratios; The power of ν is fixed to be ν =9/8, but we cannot find Anupama & Sethi− (1994) obtained E(B V )=1.1 from any significant differences in the light curves even if we the He I line ratios, and E(B V )=1.2 from− the Hα/Hβ adopt ν = 2. Finally, we have study the effect of hydro- ratio. − gen content. Either the early- or late-phase V light curve Finally, we examine the possibility of Type Ia supernova depends hardly on the hydrogen content X as shown in explosion. The envelope mass at the maximum V light is −6 Figure 18. The optically thick, outburst wind stops 56, estimated to be ∆M =3.1 10 M⊙ from our solutions, × −7 68, and 80 days after maximum, and the steady hydrogen indicating a mass transfer rate of M˙ acc = 1.1 10 M⊙ shell burning ends 68, 89, and 122 days after maximum, yr−1 between the 1962 and 1990 outbursts.× Moreover, for X =0.35, 0.50, and 0.70, respectively. about 90% of the envelope mass has been blown off in −7 Anupama & Sethi (1994) obtained optical spectra 18— the outburst wind and the residual 10% (3 10 M⊙) has 20 days after maximum. They pointed out that these been left and added to the helium layer of the× WD. There- line features closely resembles the spectrum of RS Oph fore, we obtain the net mass increasing rate of M˙ He = −8 −1 60 days after maximum. Anupama & Sethi derived the 1.1 10 M⊙ yr . This meets condition (52) of SN Ia Zanstra temperature T∗ = 300, 000 K of the ionizing cen- explosion× (e.g., Nomoto & Kondo 1991), if the WD con- tral source from the helium and hydrogen line (Hβ and sists of carbon and oxygen. Thus, we may conclude that He II 4686 A)˚ ratios and its radius R∗ =0.06R⊙ from Hβ V3890 Sgr is an immediate progenitor of Type Ia super- line flux. For RS Oph, our model gives us Tph = 300, 000 nova if the hot component is a carbon-oxygen white dwarf. K and Rph =0.08 R⊙ about 70 days after maximum (see The white dwarf can grow to 1.378M⊙ even if mild helium the previous section), being roughly consistent with the shell flashes repeatedly occur after the helium layer grows Zanstra temperature derived by Anupama & Sethi. For to the critical mass at which helium ignites (see, e.g., Kato V3890 Sgr, however, our model suggests Tph = 83, 000 & Hachisu 1999). K and Rph =1.0 R⊙ on the day 18, being not consistent. This may indicate the importance of the contribution from 9. conclusions the shocked layer photons. Mukai et al. (1990) obtained optical spectrum on May We have developed a numerical method for calculating 14 (17 days after maximum). Comparing them with the light curves of recurrent nova outbursts. Our basic model earlier spectra observed by Wagner, Bertram, & Starrfield consists of a white dwarf with an accretion disk and a (1990) and Buckley et al. (1990), Mukai et al. pointed out red giant companion. The calculated light curves includes that V3890 Sgr had progressed to a higher ionization state contributions from the white dwarf photosphere, the irra- on that day. In our modeling of the light curve, the accre- diated accretion disk, and the irradiated red giant. We tion disk plays a key role in these days. It is very likely summarize our main results on the light curve analysis of that the conspicuous change in the spectrum 20 days the four recurrent novae with a red giant companion, T after maximum is caused by the evaporation of∼ the accre- CrB, RS Oph, V745 Sco, and V3890 Sgr. tion disk. As seen in Figure 17, the disk size decreases to . 1 R⊙ on the day 20. Then, the mean temperature of 9.1. Relevance to Type Ia Supernova Progenitors the disk surface increases∼ to & 50, 000 K. The same thing 1. happens in RS Oph around the day 60—70. The mass of white dwarf (WD) has been estimated, We thus suggest a very small size∼ of the optically thick from the decline of early-phase light curve, to be MWD = 1.37 0.01 M⊙ for T CrB, 1.35 0.01 M⊙ with solar metal- ACDK or no disk, that is, Rdisk . 0.1 R⊙, in the late ± ± licity (Z =0.02) or 1.377 0.01 M⊙ with low metallicity phase (after HJD 2,448,036, indicated by ”disk vanishes” ± (Z = 0.004) for RS Oph, 1.35 0.01 M⊙ for V745 Sco, in Figure 17). Therefore, no secondary maximum of the ± light curve is expected as well as in RS Oph and V745 Sco. and 1.35 0.01 M⊙ for V3890 Sgr. 2. The apparent± distance moduli are obtained, with the light curve fitting, to be (m M)V = 10.2 for T CrB, 8.4. The Distance, Envelope Mass of the White Dwarf, (m M) = 11.2 for RS Oph,− (m M) = 17.0 for and Possibility of SN Ia Explosion − V − V V745 Sco, and (m M)V = 14.7 for V3890 Sgr. There- Fitting indicates an apparent distance modulus of (m fore, the distances are− estimated to be d = 0.94 kpc with − M)V = 14.7. Thus, the distance to V3890 Sgr is de- AV = 0.35 for T CrB, d = 0.57 kpc with AV = 2.3 for rived to be d = 4.2 kpc if we adopt the absorption of RS Oph, d = 5.0 kpc with AV = 3.5 for V745 Sco, and A = 3.1 E(B V ) = 1.6, where E(B V ) = 0.5 d =4.2 kpc with A =1.6 for V3890 Sgr. V × − − V 20 Hachisu and Kato

3. Each envelope mass at the optical maximum is also cal- accretion disk large enough to show a secondary maximum −6 −6 culated to be ∆M 3 10 M⊙ (T CrB), 2 10 M⊙ after a tilting instability grows up. −6∼ × −6 × (RS Oph), 5 10 M⊙ (V745 Sco), 3 10 M⊙ (V3890 8. The other three systems, RS Oph, V745 Sco, and V3890 Sgr), which indicates× an average mass× accretion rate of Sgr, cannot maintain the mass transfer during the outburst −7 −1 −7 M˙ acc 0.4 10 M⊙ yr (80 yrs, T CrB), 1.2 10 M⊙ because the cool red giant wind has once been evacuated −1 ∼ × −7 −1 × yr (18 yrs, RS Oph), 0.9 10 M⊙ yr (52 yrs, V745 by the ejecta of the outburst. The light curves of RS Oph −7 −1 × Sco), 1.1 10 M⊙ yr (28 yrs, V3890 Sgr) during qui- and V3890 Sgr require a relatively large size accretion disk escent phase.× of Rdisk 10 R⊙ in the mid-phase but a small size of ∼ 4. Although a large part of the envelope mass is blown Rdisk . 1 R⊙ in the late-phase. The late-phase UV data in the wind, each white dwarf can retain a substantial of RS Oph can be well reproduced with a small accretion part of the envelope mass after hydrogen burning. Thus, disk of Rdisk 0.5R⊙. No or a very small accretion disk is ∼ we have obtained net mass-increasing rates of the WDs, also derived in V745 Sco. These accretion disks are quickly −7 −1 −7 −1 . ⊙ M˙ He 0.10 10 M⊙ yr (T CrB), 0.12 10 M⊙ yr dissipating to become as small as Rdisk 1 R , at least, ∼ × −7 −1 × −7 (RS Oph), 0.05 10 M⊙ yr (V745 Sco), 0.11 10 M⊙ in the late phase of the outbursts, because no mass is fed. yr−1 (V3890 Sgr).× × Such a small disk cannot contribute much to the V light. 5. These results strongly indicate that the white dwarfs Even if the tilting instability begins to grow after the wind in recurrent novae have now grown up to near the Chan- stops, it is therefore very unlikely that secondary maxima drasekhar mass limit and will soon explode as a Type Ia like in T CrB are observed in these three systems. 9. supernova if the white dwarf consists of carbon and oxy- The post-outburst minima of RS Oph outbursts, which gen. is 1 mag (in V ) below the quiescent level and continues for 300 days from 100 to 400 days after the optical 9.2. Appearance/Absence of Secondary Maximum maximum,∼ is understood from the evacuation of the cool 6. Among the four recurrent novae with a red giant com- red giant wind by the ejecta of nova. After the hot WD panion, only T CrB shows a secondary maximum 120—190 wind stops 80 days after the optical maximum, the cool days after the optical maximum. This secondary max- red giant wind expands to recover again in 300 days. ∼ imum can be well reproduced if an irradiated accretion Then, 400 days after maximum, mass accretion from ∼ disk of Rdisk 6 R⊙ around the white dwarf tilts due to the red giant wind to the white dwarf resumes. Supersoft a radiation-induced∼ warping instability. If we neglect the X-rays from the white dwarf, which were observed about accretion disk, we cannot fully reproduce the secondary 250 days after maximum, are very consistent with the end peak by an irradiated M-giant model. of hot WD wind and the evacuation of the cool red giant 7. The difference between T CrB and the others is in the wind. companion: it fills its inner critical Roche lobe only in T CrB but in the others does not. Therefore, in T CrB the mass transfer is mainly driven by the Roche lobe overflow This research has been supported in part by the while in the other systems matter is fed by cool red giant Grant-in-Aid for Scientific Research (08640321, 09640325, winds. In T CrB, the mass transfer continues even in the 11640226) of the Japanese Ministry of Education, Science, outburst wind phase, which maintains an optically thick Culture, and Sports.

REFERENCES

Albinson, J. S., Bode, M. F., Evans, A., & Goldsmith, M. J. 1985, Danskin, K. 1985, IAU Circ. No. 4031 IAU Circ. No. 4091 Della Valle, M., Livio, M. 1996, ApJ, 473, 240 Alcock, G. E. D. 1985, IAU Circ. No. 4031 Deutsch, A. J. 1948, PASP, 60, 120 Allen, C. W. 1973, Astrophysical Quantities (London: The Athlone Dobrzycka, D., & Kenyon, S. J. 1994, AJ, 108, 2259 Press), chap. 10 Dobrzycka, D., Kenyon, S. J., Proga, D., Miko lajewska, J., Wade, R. Anupama, G. C., & Miko lajewska, J. 1999, A&A, 344, 177 A. 1996, AJ, 111, 2090 Anupama, G. C., & Sethi, S. 1994, MNRAS, 269, 105 Ducoty, R. 1985a, IAU Circ. No. 4031 Ashbrook, J. 1946, Harvard Bull., 918, 7 Ducoty, R. 1985b, IAU Circ. No. 4036 Bailey, J. 1975, JBAS, 85, 217 Ducoty, R. 1985c, IAU Circ. No. 4066 Bailey, J. 1981, MNRAS, 197, 31 Duerbeck, H. W., & Schwarz, H. E. 1989, IAU Circ., No. 4825 Belczy´nski, K., & Miko lajewska, J. 1998, MNRAS, 296, 77 Duerbeck, H. W., Schwarz, H. E., & Augusteijn, T. 1989, IAU Circ., Bessell, M. S. 1990, PASP, 102, 1181 No. 4844 Bessell, M. S., & Brett, J. M. 1988, PASP, 100, 1134 Duerbeck, H. W., Seitter, W. C. 1990, in IAU Colloquium 122, The Boattini, A. 1990, IAU Circ., No. 5007 Physics of Classical Novae, eds. A. Cassatella & R. Viotti (Berlin: Bode, M. F. 1987, RS Ophiuchi (1985) and the Recurrent Novae Springer), 425 Phenomenon (VNU Science Press, Utrecht) Dufay, J., Bloch, M., Bertaud, Ch., & Dufay, M. 1964, Annales Bode, M. F., & Kahn, F. D., 1985, MNRAS, 217, 205 d’Astrophysique, 27, 555 Bohigas, J., Echevarr´ia, J., Diego, F., & Sarmiento, J. A. 1989, Eggleton, P. P. 1983, ApJ, 268, 368 MNRAS, 238, 1395 Ellerbe, J. 1985, IAU Circ. No. 4066 Bortle, J. 1985a, IAU Circ. No. 4048 Evans, A., Callus, C. M., Albinson, J. S., Whitelock, P. A., Glass, I. Bortle, J. 1985b, IAU Circ. No. 4049 S., Carter, B., & Roberts, G. 1988, MNRAS, 234, 755 Bortle, J. 1985c, IAU Circ. No. 4074 Fekel, F. C., Joyce, R. R., Hinkle, K. H., & Skrutskie, M. 2000, AJ, Buckley, D. A., Wargau, W. F., & Soltynski, M. G. 1990, IAU Circ., 119, 1375 No. 5019 Gonz´alez-Riestra, R. 1992, A&A, 265, 71 Cannizzo, J. K., & Kenyon, S. J. 1992, ApJ, 386, L17 Gorbatskii, V. G. 1972, Soviet Astr., 16, 32 Contini, M., Orio, M., & Prialnik, D. 1995, MNRAS, 275, 195 Hachisu, I., & Kato, M. 1999, ApJ, 517, L47 Cook, L. 1985a, IAU Circ. No. 4036 Hachisu, I., & Kato, M. 2000a, ApJ, 536, L93 Cook, L. 1985b, IAU Circ. No. 4048 Hachisu, I., & Kato, M. 2000b, ApJ, 540, 447 Cook, L. 1985c, IAU Circ. No. 4066 Hachisu, I., Kato, M., Kato, T., & Matsumoto, K. 2000a, ApJ, 528, Cordova, F. A., Mason, K. O., & Nelson, J. E. 1981, ApJ, 245, 609 L97 Recurrent Novae 21

Hachisu, I., Kato, M., Kato, T., Matsumoto, K., & Nomoto, K. Nussbaumer, H., & Stencel, R. E. 1989, Exploring the Universe with 2000b, ApJ, 534, L189 the IUE Satellite, ed. Y. Kondo (Dordrecht: Kluwer), 203 Hachisu, I., Kato, M., & Nomoto, K. 1996, ApJ, 470, L97 O’Brien, T. J., Bode, M. F., & Kahn, F. D. 1992, MNRAS, 255, 683 Hachisu, I., Kato, M., & Nomoto, K. 1999a, ApJ, 522, 487 O’Brien, T. J., & Kahn, F. D. 1987, MNRAS, 228, 277 Hachisu, I., Kato, M., Nomoto, K., & Umeda, H. 1999b, ApJ, 519, Orio, M. 1993, A&A, 274, L41 314 Overbeek, D. 1989, IAU Circ., No. 4035 Harries, T. J., & Howarth, I. D. 2000, A&A, 361, 139 Paczynski, B. 1965, Acta Astronomica, 15, 197 Harrison, T. E., Johnson, J. J., & Spyromilio, J. 1993, AJ, 105, 320 Padin, S., Davis, R. J., & Bode, M. F. 1985, Nature, 315, 306 Herbig, G. H., & Neubauer, F. J. 1946, PASP, 58, 196 Patterson, J. et al. 1998, PASP, 110, 380 Hjellming, R. M. 1989, IAU Circ., No. 4853 Payne-Gaposchkin, C., & Wright, F. W. 1946, ApJ, 104, 75 Hjellming, R. M., van Gorkom, J. H., Taylor, A. R., Seaquist, E. R., Pearce, A. 1990a, IAU Circ., No. 5002 Padin, S., Davis, R. J., & Bode, M. F. 1986, ApJ, 305, L71 Pearce, A. 1990b, IAU Circ., No. 5004 Honda, M. 1989, IAU Circ., No. 4829 Pearce, A. 1990c, IAU Circ., No. 5007 Hric, L., Petrik, K., Urban, Z., Niarchos, P., & Anupama, C. G. 1998, Pearce, A. 1990d, IAU Circ., No. 5015 A&A, 339, 449 Pearce, A. 1990e, IAU Circ., No. 5021 Hurst, G. M. 1985a, IAU Circ. No. 4031 Pearce, A. 1990f, IAU Circ., No. 5028 Hurst, G. M. 1985b, IAU Circ. No. 4048 Pereira, A. 1990a, IAU Circ., No. 5007 Ianna, P. A. 1964, ApJ, 139, 780 Pereira, A. 1990b, IAU Circ., No. 5015 Iben, I. Jr., & Tutukov, A. V. 1984, ApJS, 54, 335 Perlmutter, S. et al. 1999, ApJ, 1999, 517, 565 Iglesias, C. A., & Rogers, F. 1996, ApJ, 464, 943 Pettit, E. 1946a, PASP, 58, 153 Itoh, H., & Hachisu, I. 1990, ApJ, 358, 551 Pettit, E. 1946b, PASP, 58, 213 Jones, A. F. 1990a, IAU Circ., No. 5002 Pettit, E. 1946c, PASP, 58, 255 Jones, A. F. 1990b, IAU Circ., No. 5004 Pettit, E. 1946d, PASP, 58, 359 Jones, A. F. 1990c, IAU Circ., No. 5007 Pottasch, S. R. 1967, Bull. Astr. Insts. Neth., 19, 227 Kahabka, P., Hartmann, H. W., Parmar, A. N., & Negueruela, I. Prialnik, D. 1986, ApJ, 310, 222 1999, A&A, 374, L43 Pringle, J. E. 1996, MNRAS, 281, 357 Kahabka, P., & van den Heuvel, E. P. J. 1997, ARA&A, 35, 69 Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618 Kato, M. 1983, PASJ, 35, 507 Riess, A. G. et al. 1998, AJ, 116, 1009 Kato, M. 1990, ApJ, 355, 277 Riggs, J. 1985, IAU Circ. No. 4036 Kato, M. 1991, ApJ, 369, 471 Rosino, L. 1987, RS Ophiuchi (1985) and the Recurrent Novae Kato, M. 1994, A&A, 281, L49 Phenomenon, ed. M. F. Bode (VNU Science Press, Utrecht), 1 Kato, M. 1995, in Cataclysmic Variables, eds. A. Bianchini, M. Della Rosino, L., & Iijima, T., 1987, RS Ophiuchi (1985) and the Recurrent Valle, & M. Orio (Dordrecht: Kluwer), 243 Novae Phenomenon, ed. M. F. Bode (VNU Science Press, Utrecht), Kato, M. 1999, PASJ, 51, 525 27 Kato, M., & Hachisu, I., 1994, ApJ, 437, 802 Ruffert, M., Cannizzo, J. K., & Kenyon, S. J. 1993, ApJ, 419, 780 Kato, M., Hachisu, I., 1999, ApJ, 513, L41 Sanford, R. F. 1946, PASP, 58, 156 Kenyon, S. J., Fernandez-Castro, T., & Stencel, R. E. 1988, AJ, 95, Sanford, R. F. 1947, PASP, 59, 87 1817 Sanford, R. F. 1949, ApJ, 109, 81 Kenyon, S. J., & Garcia, M. R. 1986, AJ, 91, 125 Schaefer, B. 1990, ApJ, 355, L39 Kenyon, S., & Webbink, R. 1984, ApJ, 279, 252 Schaefer, B. E., Landold, A. U., Vogt, N., Buckley, D., Warner, B., King, R. 1985, IAU Circ. No. 4036 Walker, A. R., & Bond, H. B. 1992, ApJS, 81, 321 Kovetz, A., Prialnik, D. 1994, ApJ, 424, 319 Schaefer, B., & Ringwald, F. A. 1995, ApJ, 447, L45 Kraft, R. P. 1958, ApJ, 127, 625 Schandl, S., Meyer-Hofmeister, E., & Meyer, F. 1997, A&A, 318, 73 Li, X.-D., van den Heuvel, E. P. J. 1997, A&A, 322, L9 Schmeer, P. 1999, IAU Circ., No. 7113 Liller, W. 1989a, IAU Circ., No. 4820 Schmeer, P. 1990a, IAU Circ., No. 5015 Liller, W. 1989b, IAU Circ., No. 4826 Schmeer, P. 1990b, IAU Circ., No. 5021 Liller, W. 1990, IAU Circ., No. 5010 Schmeer, P. 1990c, IAU Circ., No. 5028 Leibowitz, E. M., Ofek, E. O., & Mattei, J. A. 1997, MNRAS, 287, Scovil, C. 1985, IAU Circ. No. 4031 634 Scott, A. D., Rawlings, J. M. C., Krautter, J., & Evans, A. 1994, Lines, H. C., Lines, R. D., & McFaul, T. G. 1988, AJ, 95, 1505 MNRAS, 268, 749 Livio, M. 2000, Type Ia Supernovae: Theory and Cosmology, Sekiguchi, K. 1990, IAU Circ., No. 5047 (Cambridge: Cambridge Univ. Press), in press (astro-ph/9903264) Sekiguchi, K., Whitelock, P. A., Feast, M. W., Barrett, P. E., Livio, M., & Pringle, J. E. 1996, ApJ, 465, L55 Caldwell, J. A. R., Carter, B. S., Catchpole, R. M., Laing, J. D., Livio, M., & Truran, J. W. 1992, ApJ, 389, 695 Laney, C. D., Marang, F., & van Wyk, F. 1990, MNRAS, 246, 78 Livio, M., Truran, J. W., & Webbink, R. F. 1986, ApJ, 308, 736 Sekiguchi, K. 1992, Variable Stars and Galaxies, ed. B. Warner (ASP Lloyd, H. M., Bode, M. F., O’Brien, T. J., & Kahn, F. D. 1993, Conference Series), Vol. 30, p. 345 MNRAS, 265, 457 Selvelli, P. L., Cassatella, A., & Gilmozzi, R. 1992, ApJ, 393, 289 Lowder, W. 1985, IAU Circ. No. 4031 Shahbaz, T., Somers, M., Yudin, B., & Naylor, T. 1997, MNRAS, Lubbock, S. 1985a, IAU Circ. No. 4036 288, 1027 Lubbock, S. 1985b, IAU Circ. No. 4048 Shahbaz, T., Hauschildt, P. H., Naylor, T., & Ringwald, F. 1999, Makino, F. 1989, IAU Circ., No. 4285 MNRAS, 306, 675 Mason, K. O., C´ordova, F. A., Bode, M. F., Barr, P. 1987, RS Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 Ophiuchi (1985) and the Recurrent Novae Phenomenon, ed. M. Shore, S. N., Kenyon, S. J., Starrfield, S., & Sonneborn, G. 1996, F. Bode (VNU Science Press, Utrecht), 167 ApJ, 456, 717 Maloney, P., Begelman, M. C., & Pringle, J. E. 1996, ApJ, 472, 582 Shore, S., Sonneborn, G., & Starrfield, S. 1989, IAU Circ., No. 4826 Matsumoto, K., Kato, T., & Hachisu, I. 2000, PASJ, submitted Smith, V. V., Cunha, K., Jorissen, A., & Boffin, H. M. J. 1996, A&A, McLaughlin, D. B., 1946, PASP, 58, 159 315, 179 McNaught, R. H. 1985a, IAU Circ. No. 4031 Snijders, M. A. J. 1987a, RS Ophiuchi (1985) and the Recurrent McNaught, R. H. 1985b, IAU Circ. No. 4036 Novae Phenomenon, ed. M. F. Bode (VNU Science Press, Utrecht), McNaught, R. H. 1985c, IAU Circ. No. 4048 p. 51 McNaught, R. H. 1985d, IAU Circ. No. 4066 Snijders, M. A. J. 1987b, Ap&SS, 130, 243 Medway, K. 1985, IAU Circ. No. 4048 Southwell, K. A., Livio, M., & Pringle, J. E. 1997, ApJ, 478, L29 Mennickent, R. E.; Honeycutt, R. K. 1995, IBVS, No. 4232 Spoelstra, T. A. T., Taylor, A. R., Pooley, G. G., Evans, A., & Morgan, W. W., & Deutsch, A. J. 1947, ApJ, 106, 362 Albinson, J. S. 1987, MNRAS, 224, 791 Mukai, K., Rosen, S. R., Supelli, K., & Allen, D. 1990, IAU Circ., Starrfield, S., Sparks, W. M., & Shaviv, G. 1988, ApJ, 325, L35 No. 5015 Starrfield, S., Sparks, W. M., & Truran, J. W. 1985, ApJ, 291, 136 Munari, U., Renzini, A. 1992, ApJ, 397, L87 Sventek, P. 1985, IAU Circ. No. 4066 Munari, U., Zwitter, T., Tomov, T., Bonifacio, P., Molaro, P., Taylor, A. R.; Davis, R. J.; Porcas, R. W.; Bode, M. F. 1989, Selvelli, P., Tomasella, L., Niedzielski, A., & Pearce, A. 1999, A&A, MNRAS, 237, 81 374, L39 Umeda, H., Nomoto, K., Kobayashi, C., Hachisu, I., & Kato, M. M¨urset, U., & Schmid, H. M. 1999, A&AS, 137, 473 1999a, ApJ, 522, L43 Nomoto, K. 1982, ApJ, 253, 798 Umeda, H., Nomoto, K., Yamaoka, H., Wanajo, S. 1999b, ApJ, 513, Nomoto, K., Nariai, K., & Sugimoto, D. 1979, PASJ, 31, 287 861 Nomoto, K., & Kondo, Y. 1991, ApJ, 367, L19 Verdenet, M. 1985, IAU Circ. No. 4066 Nomoto, K., Thielemann, F., & Yokoi, K. 1984, ApJ, 286, 644 Wallerstein, G. 1958, PASP, 70, 537 Wallerstein, G. 1961, PASP, 73, 153 22 Hachisu and Kato

Wallerstein, G., & Garnavich, P. M. 1985, PASP, 98, 875 Weier, D. 1985, IAU Circ. No. 4036 Wagner, R. M., Bertram, R., & Starrfield, S. 1990, IAU Circ., No. Williams, R. E. 1992, AJ, 104, 725 5006 Williams, R. E. 1994, ApJ, 426, 279 Wagner, R. M., Bertram, R., Starrfield, S., & Shore, S. 1989, IAU Williams, R. E., Hamuy, M., Phillips, M. M., Heathcote, S. R., Wells, Circ., No. 4822 L., & Navarrette, M. 1991, ApJ, 376, 721 Warner, B. 1995, Cataclysmic Variable Stars, (Cambridge: Williams, R. E., Phillips, M. M., & Hamuy, M. 1994, ApJS, 90, 297 Cambridge University Press), chaps. 4 and 5 Williams, R. E., Sparks, W. M., Gallagher, J. S., Ney, E. P., Webbink, R. F. 1976, Nature, 262, 271 Starrfield, S. G., & Truran, J. W. 1981, ApJ, 251, 221 Webbink, R. F. 1984, ApJ, 277, 355 Wright, F. W. 1946, Harvard Bull., 918, 9 Webbink, R. F., Livio, M., Truran, J. W., & Orio, M. 1987, ApJ, Yudin, B, Munari, U. 1993, A&A, 270, 165 314, 653 Zamanov, R. K., & Bruch, A. 1998, A&A, 338, 988 Recurrent Novae 23

Table 1 Recurrent novae

Star Subclass Porb Outburst Dates (day) RS Oph RS Oph 457 1898,1933,1958,1967,1985 T CrB 227.6 1866,1946 V745 Sco — 1937,1989 V3890 Sgr — 1962,1990 U Sco U Sco 1.23 1863,1906,1936,1979,1987,1999 V394 CrA 0.757 1949,1987 LMC RN — 1968,1990 CI Aql 0.618 1917,2000 T Pyx T Pyx 0.076 1890,1902,1920,1944,1966

Table 2 Physical properties of Recurrent novae outbursts

Star Porb MWD X ∆Mmax τrec M˙ acc ηH M˙ He E(B−V ) d Fate −1 −1 (day) (M⊙) (M⊙) (yr) (M⊙ yr ) (M⊙ yr ) (kpc) 6 8 8 T CrB 227.67 1.37 0.70 3 × 10− 80 4.1 × 10− 0.10 0.4 × 10− 0.11 0.9 SN Ia a 6 7 8 RS Oph 460 1.35 0.70 2 × 10− 18 1.2 × 10− 0.10 1.2 × 10− 0.73 0.6 SN Ia b 6 7 8 RS Oph 460 1.377 0.70 2 × 10− 18 1.2 × 10− 0.10 1.2 × 10− 0.73 0.6 SN Ia 6 8 8 V745 Sco ··· 1.35 0.70 5 × 10− 52 9.0 × 10− 0.05 0.5 × 10− 1.1 5 SNIa 6 7 8 V3890 Sgr ··· 1.35 0.70 3 × 10− 28 1.1 × 10− 0.10 1.1 × 10− 0.5 4 SNIa c 6 7 7 U Sco 1.23056 1.37 0.05 3 × 10− 12 2.5 × 10− 0.40 1.0 × 10− 0.20 6 SN Ia d 6 7 8 V394 CrA 0.7577 1.37 0.05 6 × 10− 39 1.5 × 10− 0.23 3.4 × 10− 0.26 4 SN Ia

a taken from Hachisu & Kato (2000a) of solar metallicity model for RS Oph b the present low metallicity (Z = 0.004) model for RS Oph c taken from Hachisu et al. (2000a) d taken from Hachisu & Kato (2000b) 24 Hachisu and Kato

Optically Thick Winds

Fig. 1.— Optically thick winds blow from mass-accreting white dwarfs (left figure) when the mass transfer rate from a lobe-filling companion (right figure) exceeds a critical rate, i.e., M˙ acc > M˙ cr. Here, we assume that the white dwarf accretes mass from equatorial region and, at the same time, blows winds from polar regions as illustrated in the figure. Recurrent Novae 25

Fig. 2.— Response of white dwarfs to mass accretion is illustrated in the white dwarf mass and the mass accretion rate plane, i.e., in MWD-M˙ acc plane. Strong optically thick winds blow above the line of M˙ acc > M˙ cr. The wind mass loss rate is M˙ wind ≈ M˙ acc − M˙ cr. Steady hydrogen shell burning with no optically thick winds occur between M˙ std ≤ M˙ acc ≤ M˙ cr. There is no steady-state burning below M˙ acc < M˙ std. Instead, intermittent shell flashes occur. The envelope mass at which hydrogen shell flash occurs is also shown (taken from Fig. 9 of Nomoto 1982). 26 Hachisu and Kato

Fig. 3.— Regions producing SNe Ia are plotted in the log P − Md (orbital period — donor mass) plane for the WD+MS system (left in the figure) and the WD+RG system (right in the figure). Here we assume an initial white dwarf mass of 1.0M⊙. The initial system inside the region encircled by a thin solid line (labelled by ”initial”) is increasing its white dwarf mass up to the critical mass (MIa = 1.378M⊙) for Type Ia supernova explosion, the regions of which are encircled by a thick solid line (labelled by ”final”). Currently known positions of each recurrent nova are indicated by a star mark (⋆) for U Sco (e.g., Hachisu et al. 2000a), a filled triangle for T CrB (e.g., Belczy´nski & Miko lajewska 1998), an open rectangle for T Pyx (e.g., Patterson et al. 1998), but by arrows for the other three recurrent novae, V394 CrA, CI Aql, and RS Oph, because the mass of the companion is not yet available explicitly. Two subclasses of the recurrent novae, U Sco subclass and RS Oph (or T CrB) subclass correspond to the helium-rich supersoft X-ray source channel (Hachisu et al. 1999b) and the symbiotic channel (Hachisu et al. 1999a) of Type Ia supernovae, respectively. Recurrent Novae 27

Fig. 4.— Wind mass loss rate (dashed lines) and total mass decreasing rate of the envelope (solid lines labeled by ∆M), i.e., wind mass loss rate plus nuclear burning rate, M˙ wind + M˙ nuc, are plotted against the envelope mass for WDs with masses of 0.8, 1.0, 1.2, 1.3 and 1.377M⊙. There exist only wind solutions above the breaks on the solid lines while only static (no wind) solutions exist below the breaks. Thus, an optically thick wind stops at this break. The photospheric radius, Rph, the photospheric temperature, Tph, and the photospheric velocity, vph, are also plotted against the total decreasing rate of the envelope mass in the figure. We have shown four different cases of hydrogen content of the envelope material, i.e., (a) X = 0.35, (b) 0.15, (c) 0.10, and (d) 0.05, from top to bottom panel, for the solar metallicity of Z = 0.02. Each sequence is plotted from the wind solution with log Tph = 4.0 to the static solution with M˙ nuc = M˙ std, i.e., the lower limit for steady hydrogen burning. 28 Hachisu and Kato

Fig. 5.— Same as Fig. 4, but for low metallicities of Z = 0.01 (upper two panels) and Z = 0.004 (lower two panels), i.e., (a) X = 0.35, Z = 0.01, (b) X = 0.10, Z = 0.01, (c) X = 0.35, Z = 0.004, and (d) X = 0.10, Z = 0.004. Recurrent Novae 29

Fig. 6.— Model configuration near the second peak of the recurrent nova T CrB. The cool component (right in the figure) is a red giant filling up its inner critical Roche lobe. The hemisphere is heated up by the hot component (1.37 M⊙ white dwarf, left in the figure). The photospheric radius of the hot component near the second peak is as small as ∼ 0.003 R⊙, about ∼ 0.00003 times the size of the cool component, but it is exaggerated in this figure to easily see it. We assume a warping accretion disk around the hot component, which is precessing at a period of about 140 days. The surface of the accretion disk is also heated up by the hot component. The radius of the accretion disk near the second peak is as small as ∼ 6 R⊙, i.e., about ∼ 0.07 times the inner critical Roche lobe size. The accretion disk is also exaggerated in this figure to easily see it. 30 Hachisu and Kato

D

n 4 n n 3 S 3 C S 4 n G 2 A n 1 S 2 S 1

B

Fig. 7.— Each patch is divided into four triangle parts as shown in the figure, because each patch is almost a rectangular, but strictly speaking, not a rectangular for the Roche geometry of the companion star and for the warping disk surface. So that the normal unit vector to each patch, n, is approximately calculated from the four triangle parts by equation (22). Recurrent Novae 31

C

D k

n G n j j k k n i r G i ji B j G i A

Fig. 8.— When patch k blocks the light from patch j to patch i, we do not sum up the contribution from j to i as shown in equations (29), (30), and (31). 32 Hachisu and Kato

Fig. 9.— Model light curves are plotted against time (HJD 2,430,000+) together with the observational points (open circles: Pettit 1946a,b,c,d). Dotted: the total V light of the white dwarf (WD) photosphere and the red giant (RG) photosphere without irradiation. Large arrows attached by ”min” indicate epochs at the spectroscopic conjunction with the M-giant in front. An ellipsoidal variation can be seen in the late phase. Dash-dotted: the total V light of the WD photosphere and the RG photosphere irradiated by the WD. The observed second peak cannot be reproduced only by the irradiation of the M-giant. Solid: the total V light of the WD photosphere, the RG photosphere with irradiation, and the accretion disk (ACDK) surface heated-up by the hot component. The parameters specifying the light curves are shown in the figure. Recurrent Novae 33

Fig. 10.— Model V light curves for various metallicities are plotted against time (HJD 2,446,000+) together with the observational points of RS Oph outbursts. Filled circles indicate observational points with the previous outbursts superposed (taken from Rosino 1987, see also Hachisu & Kato 2000a). Open circles correspond to the observational points of the 1985 outburst, which are taken from IAU Circulars (Scovil 1985; Lowder 1985; Alcock 1985; Hurst 1985a,b; Danskin 1985; Ducoty 1985a,b,c; McNaught 1985a,b,c,d; Lubbock 1985a,b; King 1985; Weier 1985; Riggs 1985; Cook 1985a,b,c; Medway 1985; Bortle 1985a,b,c; Ellerbe 1985; Verdenet 1985; Sventek 1985; Albinson et al. 1985), i.e., No.4031, 4036, 4048, 4049, 4066, 4074, and 4091. The WD mass is assumed to be the extreme case of MWD = 1.377 M⊙ just before the SN Ia explosion. The metallicity is attached to each line (for Z = 0.001, Z = 0.004, Z = 0.01, Z = 0.02). Each line denotes the V light of the WD photosphere only. The hydrogen content of the WD envelope is assumed to be X = 0.70 for all models. The apparent distance modulus is (m − M)V = 11.2. The lowest metallicity case of Z = 0.001 cannot reproduce the rapid decline of the V light curve of RS Oph outburst. Since the hydrogen content, X, hardly affect the decline rate of the V light curve (e.g., Kato 1999), we may conclude that Z & 0.004 in RS Oph outbursts. 34 Hachisu and Kato

Fig. 11.— Same as Fig. 10 but for the entire outburst phase for MWD = 1.377M⊙, X = 0.7, and Z = 0.004. Dashed line: Sum of the contributions from the white dwarf (WD) photosphere and the non-irradiated red giant (RG) photosphere with the surface temperature of Tph = 3400 K. The RG lies well within the inner critical Roche lobe, i.e., its radius is 0.25 times the Roche lobe size (γ = 0.25). Dashed- doted: Sum of the WD photosphere and the irradiated RG photosphere. Solid: the total V light of the WD photosphere, the irradiated RG photosphere and the accretion disk (ACDK) surface. Dotted: ultra-violet (UV) radiation (911–3250 A)˚ from the WD photosphere, the irradiated RG photosphere, and the irradiated ACDK surface. The UV data (large open circles) are taken from Snijders (1987a). Here both the disk shape and size are assumed to be constant during the outburst, i.e., α = 0.01 (Rdisk ∼ 14 R⊙) and β = 0.05 for the upper solid (and upper dotted) line, and α = 0.008 (Rdisk ∼ 1 R⊙) and β = 0.05 for the lower solid (and lower dotted) line, respectively. Two epochs with ”wind stops” and ”end of steady H burning” are indicated by arrows. The optically thick wind blows during the period from the first phase of the outburst (HJD 2,446,092) to 79 days after maximum (HJD 2,446,170). The steady hydrogen shell burning ends at 112 days after maximum (HJD 2,446,203). Recurrent Novae 35

Fig. 12.— Same as Fig. 11 but for a case of gradually evaporating accretion disk. Here, we assume that the radius of the accretion disk is gradually shrinking from ∼ 14R⊙ (α = 0.1) to ∼ 1R⊙ (α = 0.008) for X = 0.7, to ∼ 2R⊙ (α = 0.012) for X = 0.5, to ∼ 3R⊙ (α = 0.02) for X = 0.35, during the wind phase. Its time variation is exponential as given by equation (56). 36 Hachisu and Kato

Fig. 13.— Model light curves are plotted against time (HJD 2,447,000+) together with the observational points of the V745 Sco 1989 outburst. Here, the inclination angle of i = 80◦ is assumed. The other system parameters are shown in the figure. Filled circles indicate the observational points with V filter (taken from Sekiguchi et al. 1990). Crosses (taken from Liller 1989a,b) are shifted by 0.8 mag down because they are systematically 0.8 mag brighter than the V magnitude. Small allows (downwards) indicate an upper limit before the 1989 outburst (taken from Honda 1989). Open circles are visual magnitudes taken from Shore et al. (1989) and Overbeek (1989). Dashed line denotes the V light curve of the WD photosphere and the nonirradiated RG only. Dash-dotted line corresponds to the V light from the WD photosphere and the irradiated RG. Solid line represents the total V light of the WD photosphere, the irradiated RG, and the irradiated accretion disk. The photospheric radius of the WD is indicated above the V light curve. Two epochs are indicated by arrows with ”wind stops” and ”end of steady H-burning.” The optically thick wind stops at HJD 2,447,819. The steady hydrogen shell burning ends at HJD 2,447,860. Fitting with the observations indicates an apparent distance modulus of (m − M)V = 17.0. Recurrent Novae 37

Fig. 14.— Same as Fig. 13 but for different inclination angles (i = 30, 60, 80◦ from top to bottom). Dot-dashed lines denote the model light curves for the cases of no accretion disk while solid lines denote the cases with an accretion disk. The model light curves with an accretion disk depend largely on the inclination angle, while we do not find any significant differences for the cases of no accretion disk. Dotted lines indicate a different ephemeris of φmin = 500 day, which seems not to be compatible with the observation for all three inclination angles. 38 Hachisu and Kato

Fig. 15.— Same as Fig. 13 but for three different hydrogen content of X = 0.35, 0.50, and 0.70. We do not determine the hydrogen content from these observational points, but may conclude that it is somewhere between X = 0.5 and 0.7. Recurrent Novae 39

Fig. 16.— Model light curves are plotted against time (HJD 2,448,000+) together with the observational points of the V3890 Sgr 1990 outburst. Here, the inclination angle of i = 30◦ is assumed. The other system parameters are shown in the figure. Open circles indicate the observational points taken from IAU Circulars (Jones 1990a,b,c; Pearce 1990a,b,c,d,e,f; Pereira 1990a,b; Boattini 1990; Liller 1990; Schmeer 1990a,b,c) and Figure 1 of Gonz´alez-Riestra (1992), while small downward arrows show an upper limit. Solid lines denote the total V light of the WD photosphere and the nonirradiated M-giant secondary. Dot-dashed lines depict the total V light of the WD photosphere and the irradiated M-giant secondary with two different φmin, i.e., φmin = 50 day (lower) and φmin = 250 day (upper) only for MWD = 1.35 M⊙. Fitting with the observations suggests an apparent distance modulus of (m−M)V = 14.7. Here, we assume that the WD photosphere reached its maximum expansion on HJD 2,448,009. 40 Hachisu and Kato

Fig. 17.— Same as Fig. 16 but for two white dwarf masses, i.e., 1.35 M⊙ (solid lines) and 1.37 M⊙ (dashed lines). Both solid and dashed lines represent the total V light of the WD photosphere, the irradiated M-giant, and the irradiated accretion disk, the size of which is fixed during the outburst as defined by equations (11) and (12). The optically thick wind blows during the outburst, from the beginning of the outburst (HJD 2,448,008) to 80 days after maximum (HJD 2,448,089), and then the steady hydrogen shell-burning ends at 122 days after maximum (HJD 2,448,131) for the 1.35 M⊙ WD and the hydrogen content of X = 0.7, while 57 days and 78 days, respectively, for the 1.37 M⊙ WD and the hydrogen content of X = 0.7. The slow decline just after the optical maximum can be reproduced by the relatively large disk size around the 1.35 M⊙ WD but not by the 1.37 M⊙ WD. Recurrent Novae 41

Fig. 18.— Same as Fig. 17 but for the case of evaporating accretion disk with three different hydrogen content of X = 0.35, 0.50, and 0.70 (solid lines). We do not pose any constraints on the hydrogen content from these observational points, because of a lack of data points in the later phase. For comparison, two other cases with different irradiation efficiencies, ηDK = 1.0 (dot-dashed) and 0.25 (dotted), are added only for X = 0.70.