HST/STIS Spectroscopy and Modeling of the Long Term Cooling of WZ Sagittae following the July 2001 Outburst

We present the last HST/STIS E140M FUV spectrum (1150-1725A) of the dwarf nova (DN) WZ Sge, obtained in July 2004, 3 years following the early superoutburst of July 2001. Single white dwarf (WD) synthetic spectral fits (log{g}=8.5) to the data indicate that the WD has a temperature T~15,000K, about ~1500K above its quiescent temperature and it is still showing the effect of the outburst. Taking into account temperature estimates of the earlier phase of the cooling, we model the cooling curve of WZ Sge, over a period of 3 years, using a stellar evolution code including accretion and the effects of compressional heating. Assuming that compressional heating alone is the source of the energy released during the cooling phase, we find that (1) the mass of the white dwarf must be quite large (~1.0 Msun); and (2) the mass accretion rate must have a time-averaged (over 52 days of outburst) value of the order of 1.E-8 Msun/yr or larger. The outburst mass accretion rate derived from these compressional heating models is larger than the rates estimated from optical observations and from a FUV spectral fit by up to one order of magnitude. This implies that during the cooling phase the energy released by the WD is not due to compressional heating alone.