The hydrogen-rich envelopes accreted by white dwarf stars from their red dwarf companions lead to thermonuclear runaways observed as classical nova eruptions peaking at up to 1 Million solar luminosities. Virtually all nova progenitors are novalike binaries exhibiting high rates of mass transfer to their white dwarfs before and after an eruption. It is a puzzle that binaries indistinguishable from novalikes, but with much lower mass transfer rates, and resulting dwarf nova outbursts, co-exist at the same orbital periods. Nova shells surrounding several dwarf novae demonstrate that at least some novae become dwarf novae between successive nova eruptions, though the mechanisms and timescales governing mass transfer rate variations are poorly understood. Here we report simulations of the multiGyr evolution of novae which self-consistently model every eruption's thermonuclear runaway, mass and angular momentum losses, feedback due to irradiation and variable mass transfer, and orbital size and period changes. The simulations reproduce the observed wide range of mass transfer rates at a given orbital period, with large and cyclic changes in white dwarf-red dwarf binaries emerging on kyr to Myr timescales. They also demonstrate that deep hibernation, (complete stoppage of mass transfer for long periods), occurs only in short-period binaries; that initially very different binaries converge to become nearly identical systems; that while almost all prenovae should be novalike binaries, dwarf novae should also occasionally be observed to give rise to novae; and that the masses of white dwarfs decrease only slightly while their red dwarf companions are consumed.

The supersoft X-ray and optical transient ASASSN-16oh has been interpreted by Maccarone et al. (2019) as having being induced by an accretion event on a massive white dwarf, resembling a dwarf nova super-outburst. These authors argued that the supersoft X-ray spectrum had a different origin than in an atmosphere heated by shell nuclear burning, because no mass was ejected. We find instead that the event's timescale and other characteristics are typical of non-mass ejecting thermonuclear runaways, as already predicted by Shara et al. (1977) and the extensive grid of nova models by Yaron et al. (2005). We suggest that the low X-ray and bolometric luminosity in comparison to the predictions of the models of nuclear burning are due to an optically thick accretion disk, hiding most of the white dwarf surface. If this is the case, we calculated that the optical transient can be explained as a non-ejective thermonuclear event on a WD of $\simeq$1.1M$_\odot$ accreting at the rate of $\simeq3.5{-}5{\times}10^{-7}$M$_\odot$yr$^{-1}$. We make predictions that should prove whether the nature of the transient event was due to thermonuclear burning or to accretion; observational proof should be obtained in the next few years, because a new outburst should occur within $\simeq$10-15 years of the event.

Models have long predicted that the frequency-averaged masses of white dwarfs in Galactic classical novae are twice as large as those of field white dwarfs. Only a handful of dynamically well-determined nova white dwarf masses have been published, leaving the theoretical predictions poorly tested. The recurrence time distributions and mass accretion rate distributions of novae are even more poorly known. To address these deficiencies, we have combined our extensive simulations of nova eruptions with the Strope et al (2010) and Schaefer et al (2010) databases of outburst characteristics of Galactic classical and recurrent novae to determine the masses of 92 white dwarfs in novae. We find that the mean mass (frequency averaged mean mass) of 82 Galactic classical novae is 1.06 (1.13) Msun, while the mean mass of 10 recurrent novae is 1.31 Msun. These masses, and the observed nova outburst amplitude and decline time distributions allow us to determine the long-term mass accretion rate distribution of classical novae. Remarkably, that value is just 1.3 x 10^{-10} Msun/yr, which is an order of magnitude smaller than that of cataclysmic binaries in the decades before and after classical nova eruptions. This predicts that old novae become low mass transfer rate systems, and hence dwarf novae, for most of the time between nova eruptions. We determine the mass accretion rates of each of the 10 known Galactic RN, finding them to be in the range 10^{-7} - 10^{-8} $ Msun/yr. We are able to predict the recurrence time distribution of novae and compare it with the predictions of population synthesis models.

Ten weeks of daily imaging of the giant elliptical galaxy M87 with the Hubble Space Telescope (HST) has yielded 41 nova light curves of unprecedented quality for extragalactic cataclysmic variables. We have recently used these light curves to demonstrate that the observational scatter in the so-called Maximum-Magnitude Rate of Decline (MMRD) relation for classical novae is so large as to render the nova-MMRD useless as a standard candle. Here we demonstrate that a modified Buscombe - de Vaucouleurs hypothesis, namely that novae with decline times t2 > 10 days converge to nearly the same absolute magnitude about two weeks after maximum light in a giant elliptical galaxy, is supported by our M87 nova data. For 13 novae with daily-sampled light curves, well determined times of maximum light in both the F606W and F814W filters, and decline times $t2 > 10 days we find that M87 novae display M(606W,15) = -6.37 +/- 0.46 and M(814W,15) = -6.11 +/- 0.43. If very fast novae with decline times t2 < 10 days are excluded, the distances to novae in elliptical galaxies with stellar binary populations similar to those of M87 should be determinable with 1 sigma accuracies of +/-20% with the above calibrations.

We describe Castalia, a proposed mission to rendezvous with a Main Belt Comet (MBC), 133P/Elst-Pizarro. MBCs are a recently discovered population of apparently icy bodies within the main asteroid belt between Mars and Jupiter, which may represent the remnants of the population which supplied the early Earth with water. Castalia will perform the first exploration of this population by characterising 133P in detail, solving the puzzle of the MBC's activity, and making the first in situ measurements of water in the asteroid belt. In many ways a successor to ESA's highly successful Rosetta mission, Castalia will allow direct comparison between very different classes of comet, including measuring critical isotope ratios, plasma and dust properties. It will also feature the first radar system to visit a minor body, mapping the ice in the interior. Castalia was proposed, in slightly different versions, to the ESA M4 and M5 calls within the Cosmic Vision programme. We describe the science motivation for the mission, the measurements required to achieve the scientific goals, and the proposed instrument payload and spacecraft to achieve these.

Horizontal branch (HB) morphology is a complex multiple-parameter problem. Besides the metallicity, two other leading parameters are the mass loss rate (MLR) and the initial He abundance of the HB progenitors. Using the STAREV stellar evolution code, we produce a wide array of Extreme Horizontal Branch (EHB) stars and also examine their post-HB evolution. EHB stars are produced in our calculations by the so called `delayed (late) hot core flash' scenario. The MLR is increased on the red giant branch (RGB) to the extent that, prior to reaching core flash conditions, only a very thin H-rich envelope remains and helium ignition takes place at hotter positions on the HRD. We perform an extensive, self-consistent parameter study, covering populations I and II (Z=0.0001-0.03), for both normal initial helium abundances and He-enriched models (up to Y=0.40). For completeness of the study and in order to point to complete trends, we chose NOT to cut out several combinations (or results of) that may extend beyond realistic limits. We present results and complete evolutionary tracks for the covered parameter space, showing in particular that: a) Increased He abundance ON ITS OWN -- without having a significant-enough MLR on the RGB -- DOES NOT lead to the production of EHB stars; however, b) The bluest (hottest) HB positions do result from the COMBINED effect of He-enhancement and increased MLR; c) The general trend is that the effective temperature on the HB increases with decreasing metallicity, but there is an indication for a halt, or even a reversal of this trend, as Z further decreases below $10^{-3}$; d) EHB stars can serve as major contributors to the UV flux emanating from their host system. Thus, the present comprehensive study both complements and lends support to previous, more restricted studies of the HB phase, and adds results for unexplored regions of the parameter space.

We examine the possibility that icy super-Earth mass planets, formed over long time scales (0.1--1~Gyr) at large distances ($\sim$ 200--1000~AU) from their host stars, will develop massive H-rich atmospheres. Within the interior of these planets, high pressure converts CH$_4$ into ethane, butane, or diamond and releases H$_2$. Using simplified models which capture the basic physics of the internal structure, we show that the physical properties of the atmosphere depend on the outflux of H$_2$ from the mantle. When this outflux is $\lesssim 10^{10}$ [molec cm$^{-2}$ s$^{-1}$], the outgassed atmosphere has base pressure $\lesssim$ 1 bar. Larger outflows result in a substantial atmosphere where the base pressure may approach $10^3 - 10^4$ bar. For any pressure, the mean density of these planets, 2.4--3 [g cm$^{-3}$], is much larger than the mean density of Uranus and Neptune, 1.3--1.6 [g cm$^{-3}$]. Thus, observations can distinguish between a Planet Nine with a primordial H/He-rich atmosphere accreted from the protosolar nebula and one with an atmosphere outgassed from the core.

The extensive grid of numerical simulations of nova eruptions of Yaron et al.(2005) first predicted that some classical novae might deviate significantly from the Maximum Magnitude - Rate of Decline (MMRD) relation, which purports to characterise novae as standard candles. Kasliwal et al. (2011) have announced the observational detection of an apparently new class of faint, fast classical novae in the Andromeda galaxy. These objects deviate strongly from the MMRD relationship, exactly as predicted by Yaron et al. (2005). Shara et al. (2016) recently reported the first detections of faint, fast novae in M87. These previously overlooked objects are as common in the giant elliptical galaxy M87 as they are in the giant spiral M31; they comprise about 40% of all classical nova eruptions and greatly increase the observational scatter in the MMRD relation. We use the extensive grid of nova simulations of Yaron et al. (2005) to identify the underlying causes of the existence of faint, fast novae. These are systems which have accreted, and can thus eject, only very low mass envelopes, of order 10^-7 - 10^-8 Msun, on massive white dwarfs. Such binaries include, but are not limited to, the recurrent novae. These same models predict the existence of ultrafast novae which display decline times t2 as short as five hours. We outline a strategy for their future detection.

We develop a long-term 1-D evolution model for icy satellites that couples multiple processes: water migration and differentiation, geochemical reactions and silicate phase transitions, compaction by self-gravity, and ablation. The model further considers the following energy sources and sinks: tidal heating, radiogenic heating, geochemical energy released by serpentinization or absorbed by mineral dehydration, gravitational energy and insolation, and heat transport by conduction, convection, and advection. We apply the model to Enceladus, by guessing the initial conditions that would render a structure compatible with present-day observations, assuming the initial structure to have been homogeneous. Assuming the satellite has been losing water continually along its evolution, we postulate that it was formed as a more massive, more icy and more porous satellite, and gradually transformed into its present day state due to sustained long-term tidal heating. We consider several initial compositions and evolution scenarios and follow the evolution for the age of the Solar System, testing the present day model results against the available observational constraints. Our model shows the present configuration to be differentiated into a pure icy mantle, several tens of km thick, overlying a rocky core, composed of dehydrated rock at the center and hydrated rock in the outer part. For Enceladus, it predicts a higher rock/ice mass ratio than previously assumed and a thinner ice mantle, compatible with recent estimates based on gravity field measurements. Although, obviously, the model cannot be used to explain local phenomena, it sheds light on the internal structure invoked in explanations of localized features and activities.

Can a white dwarf, accreting hydrogen-rich matter from a non-degenerate companion star, ever exceed the Chandrasekhar mass and explode as a type Ia supernova? We explore the range of accretion rates that allow a white dwarf (WD) to secularly grow in mass, and derive limits on the accretion rate and on the initial mass that will allow it to reach $1.4M_\odot$ --- the Chandrasekhar mass. We follow the evolution through a long series of hydrogen flashes, during which a thick helium shell accumulates. This determines the effective helium mass accretion rate for long-term, self-consistent evolutionary runs with helium flashes. We find that net mass accumulation always occurs despite helium flashes. Although the amount of mass lost during the first few helium shell flashes is a significant fraction of that accumulated prior to the flash, that fraction decreases with repeated helium shell flashes. Eventually no mass is ejected at all during subsequent flashes. This unexpected result occurs because of continual heating of the WD interior by the helium shell flashes near its surface. The effect of heating is to lower the electron degeneracy throughout the WD, and especially in the outer layers. This key result yields helium burning that is quasi-steady state, instead of explosive. We thus find a remarkably large parameter space within which long-term, self-consistent simulations show that a WD can grow in mass and reach the Chandrasekhar limit, despite its helium flashes.

T Pyxidis is the only recurrent nova surrounded by knots of material ejected in previous outbursts. Following the eruption that began on 2011 April 14.29, we obtained seven epochs (from 4 to 383 days after eruption) of Hubble Space Telescope narrowband Ha images of T Pyx . The flash of radiation from the nova event had no effect on the ejecta until at least 55 days after the eruption began. Photoionization of hydrogen located north and south of the central star was seen 132 days after the beginning of the eruption. That hydrogen recombined in the following 51 days, allowing us to determine a hydrogen atom density of at least 7e5 cm^-3 - at least an order of magnitude denser than the previously detected, unresolved [NII] knots surrounding T Pyx. Material to the northwest and southeast was photoionized between 132 and 183 days after the eruption began. 99 days later that hydrogen had recombined. Both then (282 days after outburst) and 101 days later, we detected almost no trace of hydrogen emission around T Pyx. There is a large reservoir of previously unseen, cold diffuse hydrogen overlapping the previously detected, [NII] - emitting knots of T Pyx ejecta. The mass of this newly detected hydrogen is probably an order of magnitude larger than that of the [NII] knots. We also determine that there is no significant reservoir of undetected ejecta from the outer boundaries of the previously detected ejecta out to about twice that distance, near the plane of the sky. The lack of distant ejecta is consistent with the Schaefer et al (2010) scenario for T Pyx, in which the star underwent its first eruption within five years of 1866 after many millennia of quiescence, followed by the six observed recurrent nova eruptions since 1890. This lack of distant ejecta is not consistent with scenarios in which T Pyx has been erupting continuously as a recurrent nova for many centuries or millennia.

A definitive determination of the progenitors of type Ia supernovae (SNIa) has been a conundrum for decades. The single degenerate scenario $-$ a white dwarf (WD) in a semi-detached binary system accreting mass from its secondary $-$ is a plausible path; however, no simulation to date has shown that such an outcome is possible. In this study, we allowed a WD with a near Chandrasekhar mass of $1.4M_\odot$ to evolve over tens of thousands of nova cycles, accumulating mass secularly while undergoing periodic nova eruptions. We present the mass accretion limits within which a SNIa can possibly occur. The results showed, for each parameter combination within the permitted limits, tens of thousands of virtually identical nova cycles where the accreted mass exceeded the ejected mass, i.e. the WD grew slowly but steadily in mass. Finally, the WD became unstable, the maximal temperature rose by nearly two orders of magnitude, heavy element production was enhanced by orders of magnitude and the nuclear and neutrino luminosities became enormous. We also found that this mechanism leading to WD collapse is robust, with WDs in the range $1.0 - 1.38M_\odot$, and an accretion rate of $5*10^{-7}M_\odot/yr$, all growing steadily in mass. These simulations of the onset of a SNIa event make observationally testable predictions about the light curves of pre-SN stars, and about the chemistry of SNIa ejecta.

The Palomar Transient Factory (PTF) is a synoptic survey designed to explore the transient and variable sky in a wide variety of cadences. We use PTF observations of fields that were observed multiple times (>=10) per night, for several nights, to find asteroids, construct their lightcurves and measure their rotation periods. Here we describe the pipeline we use to achieve these goals and present the results from the first four (overlapping) PTF fields analyzed as part of this program. These fields, which cover an area of 21 deg^2, were observed on four nights with a cadence of ~20 min. Our pipeline was able to detect 624 asteroids, of which 145 (~20%) were previously unknown. We present high quality rotation periods for 88 main-belt asteroids and possible period or lower limit on the period for an additional 85 asteroids. For the remaining 451 asteroids, we present lower limits on their photometric amplitudes. Three of the asteroids have lightcurves that are characteristic of binary asteroids. We estimate that implementing our search for all existing high-cadence PTF data will provide rotation periods for about 10,000 asteroids mainly in the magnitude range ~14 to ~20.

Durda et al. (2004), using numerical models, suggested that binary asteroids with large separation, called Escaping Ejecta Binaries (EEBs), can be created by fragments ejected from a disruptive impact event. It is thought that six binary asteroids recently discovered might be EEBs because of the high separation between their components (~100 > a/Rp > ~20). However, the rotation periods of four out of the six objects measured by our group and others and presented here show that these suspected EEBs have fast rotation rates of 2.5 to 4 hours. Because of the small size of the components of these binary asteroids, linked with this fast spinning, we conclude that the rotational-fission mechanism, which is a result of the thermal YORP effect, is the most likely formation scenario. Moreover, scaling the YORP effect for these objects shows that its timescale is shorter than the estimated ages of the three relevant Hirayama families hosting these binary asteroids. Therefore, only the largest (D~19 km) suspected asteroid, (317) Roxane, could be, in fact, the only known EEB. In addition, our results confirm the triple nature of (3749) Balam by measuring mutual events on its lightcurve that match the orbital period of a nearby satellite in addition to its distant companion. Measurements of (1509) Esclangona at different apparitions show a unique shape of the lightcurve that might be explained by color variations.

Extremely luminous, red eruptive variables like RV in M31 are being suggested as exemplars of a new class of astrophysical object. Our greatly extended series of nova simulations shows that classical nova models can produce very red, luminous eruptions. In a poorly studied corner of 3-D nova parameter space (very cold, low-mass white dwarfs, accreting at very low rates) we find bona fide classical novae that are very luminous and red because they eject very slowly moving, massive envelopes. A crucial prediction of these nova models - in contrast to the predictions of merging star ("mergeburst") models - is that a hot remnant, the underlying white dwarf, will emerge after the massive ejected envelope has expanded enough to become optically thin. This blue remnant must fade on a timescale of decades - much faster than a "mergeburst", which must fade on timescales of millennia or longer. Furthermore, the cooling nova white dwarf and its expanding ejecta must become redder in the years after eruption, while a contracting mergeburst must become hotter and bluer. We predict that red novae will always brighten to L~1000 Lsun for about 1 year before rising to maximum luminosity at L~10^6 - 10^7 Lsun. The maximum luminosity attainable by a nova is likely to be L~10^7 Lsun, corresponding to M-12. In an accompanying paper we describe a fading, luminous blue candidate for the remnant of M31-RV; it is observed with HST to be compatible only with the nova model.

M31-RV was an extraordinarily luminous (~10^6 Lsun) eruptive variable, displaying very cool temperatures (roughly 1000 Kelvins) as it faded. The photometric behavior of M31-RV (and several other very red novae, i.e. luminous eruptive red variables) has led to several models of this apparently new class of astrophysical object. One of the most detailed models is that of "mergebursts": hypothetical mergers of close binary stars. These are predicted to rival or exceed the brightest classical novae in luminosity, but to be much cooler and redder than classical novae, and to become slowly hotter and bluer as they age. This prediction suggests two stringent and definitive tests of the mergeburst hypothesis. First, there should always be a cool red remnant, and NOT a hot blue remnant at the site of such an outburst. Second, the inflated envelope of a mergeburst event should be slowly contracting, hence it must display a slowly rising effective temperature. We have located a luminous, UV-bright object within 0.4 arcsec (1.5 sigma of the astrometric position) of M31-RV in archival WFPC2 images taken 10 years after the outburst: it resembles an old nova. Twenty years after the outburst, the object remains much too hot to be a mergeburst. Its behavior remains consistent with that of theoretical nova models which erupt on a low mass white dwarf. Future Hubble UV and visible images could determine if the M31-RV analogs (in M85 and in M99) are also behaving like old novae.

The total masses ejected during classical nova eruptions are needed to answer two questions with broad astrophysical implications: Can accreting white dwarfs be pushed towards the Chandrasekhar mass limit to yield type Ia supernovae? Are Ultra-luminous red variables a new kind of astrophysical phenomenon, or merely extreme classical novae? We review the methods used to determine nova ejecta masses. Except for the unique case of BT Mon (nova 1939), all nova ejecta mass determinations depend on untested assumptions and multi-parameter modeling. The remarkably simple assumption of equipartition between kinetic and radiated energy (E_kin and E_rad, respectively) in nova ejecta has been invoked as a way around this conundrum for the ultra-luminous red variable in M31. The deduced mass is far larger than that produced by any classical nova model. Our nova eruption simulations show that radiation and kinetic energy in nova ejecta are very far from being in energy equipartition, with variations of four orders of magnitude in the ratio E_kin/E_rad being commonplace. The assumption of equipartition must not be used to deduce nova ejecta masses; any such "determinations" can be overestimates by a factor of up to 10,000. We data-mined our extensive series of nova simulations to search for correlations that could yield nova ejecta masses. Remarkably, the mass ejected during a nova eruption is dependent only on (and is directly proportional to) E_rad. If we measure the distance to an erupting nova and its bolometric light curve then E_rad and hence the mass ejected can be directly measured.

We present results of thermal evolution calculations for objects originating in the Kuiper belt and transferring inwards, to the region of the outer planets. Kuiper belt objects (KBOs) are considered to be part of a reservoir that supplies the flux of small icy bodies, mainly Centaurs and Jupiter-family comets, to regions interior to the orbit of Neptune. We study the internal thermal evolution, for yr, of three typical KBOs and use the end state of the simulation as initial conditions for evolutionary calculations of two typical Centaurs. Some evolutionary trends can be identified for the KBOs, depending on key physical parameters, such as size and composition. The subsequent evolution in the Centaur region results in both specific features for each modeled object (mainly surface and sub-surface composition) and common characteristics of thermally evolved Centaurs.

We present results of visible wavelengths spectroscopic measurements (0.45 to 0.72 microns) of two binary asteroids, obtained with the 1-m telescope at the Wise Observatory on January 2008. The asteroids (90) Antiope and (1509) Esclangona were observed to search for spectroscopic variations correlated with their rotation while presenting different regions of their surface to the viewer. Simultaneous photometric observations were performed with the Wise Observatory's 0.46-m telescope, to investigate the rotational phase behavior and possible eclipse events. (90) Antiope displayed an eclipse event during our observations. We could not measure any slope change of the spectroscopic albedo within the error range of 3%, except for a steady decrease in the total light flux while the eclipse took place. We conclude that the surface compositions of the two components do not differ dramatically, implying a common origin and history. (1509) Esclangona did not show an eclipse, but rather a unique lightcurve with three peaks and a wide and flat minimum, repeating with a period of 3.2524 hours. Careful measurements of the spectral albedo slopes reveal a color variation of 7 to 10 percent on the surface of (1509) Esclangona, which correlates with a specific region in the photometric lightcurve. This result suggests that the different features on the lightcurve are at least partially produced by color variations and could perhaps be explained by the existence of an exposed fresh surface on (1509) Esclangona.

We present a new stellar evolution code and a set of results, demonstrating its capability at calculating full evolutionary tracks for a wide range of masses and metallicities. The code is fast and efficient, and is capable of following through all evolutionary phases, without interruption or human intervention. It is meant to be used also in the context of modeling the evolution of dense stellar systems, for performing live calculations for both normal star models and merger-products. The code is based on a fully implicit, adaptive-grid numerical scheme that solves simultaneously for structure, mesh and chemical composition. Full details are given for the treatment of convection, equation of state, opacity, nuclear reactions and mass loss. Results of evolutionary calculations are shown for a solar model that matches the characteristics of the present sun to an accuracy of better than 1%; a 1 Msun model for a wide range of metallicities; a series of models of stellar populations I and II, for the mass range 0.25 to 64 Msun, followed from pre-main-sequence to a cool white dwarf or core collapse. An initial final-mass relationship is derived and compared with previous studies. Finally, we briefly address the evolution of non-canonical configurations, merger-products of low-mass main-sequence parents.

A 3-D numerical model of comet nuclei is presented. An implicit numerical scheme was developed for the thermal evolution of a spherical nucleus composed of a mixture of ice and dust. The model was tested against analytical solutions, simplified numerical solutions, and 1-D thermal evolution codes. The 3-D code was applied to comet 67P/Churyumov-Gerasimenko; surface temperature maps and the internal thermal structure was obtained as function of depth, longitude and hour angle. The effect of the spin axis tilt on the surface temperature distribution was studied in detail. It was found that for small tilt angles, relatively low temperatures may prevail on near-pole areas, despite lateral heat conduction. A high-resolution run for a comet model of 67P/Churyumov-Gerasimenko with low tilt angle, allowing for crystallization of amorphous ice, showed that the amorphous/crystalline ice boundary varies significantly with depth as a function of cometary latitude.

Multicycle nova evolution models have been calculated over the past twenty years, the number being limited by numerical constraints. Here we present a long-term evolution code that enables a continuous calculation through an unlimited number of nova cycles for an unlimited evolution time, even up to (or exceeding) a Hubble time. Starting with two sets of the three independent nova parameters -- the white dwarf mass, the temperature of its isothermal core, and the rate of mass transfer on to it -- we have followed the evolution of two models, with initial masses of 1 and 0.65 solar masses, accretion rates (constant throughout each calculation) of 1e-11 and 1e-9 solar-masses/yr, and relatively high initial temperatures (as they are likely to be at the onset of the outburst phase), through over 1000 and over 3000 cycles, respectively. The results show that although on the short-term consecutive outbursts are almost identical, on the long-term scale the characteristics change. This is mainly due to the changing core temperature, which decreases very similarly to that of a cooling white dwarf for a time, but at a slower rate thereafter. As the white dwarf's mass continually decreases, since both models lose more mass than they accrete, the central pressure decreases accordingly. The outbursts on the massive white dwarf change gradually from fast to moderately fast, and the other characteristics (velocity, abundance ratios, isotopic ratios) change, too. Very slowly, a steady state is reached, where all characteristics, both in quiescence and in outburst, remain almost constant. For the less massive white dwarf accreting at a high rate, outbursts are similar throughout the evolution.

We use a quasi 3-D thermal evolution model for a spherical comet nucleus, which takes into account the diurnal and latitudinal variation of the solar flux, but neglects lateral heat conduction. We model the thermal evolution and activity of Comet 9P/Tempel 1, in anticipation of the Deep Impact mission encounter with the comet. We also investigate the possible outcome of a projectile impact, assuming that all the energy is absorbed as thermal energy. An interesting result of this investigation, is that the estimated amount of dust ejected due to the impact is equivalent to 2--2.6 days of activity, during "quiet" conditions, at perihelion. We show that production rates of volatiles that are released in the interior of the nucleus depend strongly on the porous structure, in particular on the surface to volume ratio of the pores. We develop a more accurate model for calculating this parameter, based on a distribution of pore sizes, rather than a single, average pore size.

This paper is a sequel to an earlier paper devoted to multiple, multicycle nova evolution models (Prialnik & Kovetz 1995, first paper of the series), which showed that the different characteristics of nova outbursts can be reproduced by varying the values of three basic and independent parameters: the white dwarf mass-M_{WD}, the temperature of its isothermal core-T_{WD} and the mass transfer rate-Mdot. Apart from being the largest computational classical novae parameter-space survey, we show here that the parameter space is constrained by several analytical considerations and find its limiting surfaces. Consequently, we extend the grid of multicycle nova evolution models presented in Paper I almost to its limits, adding multicycle nova outburst calculations for a considerable number of new parameter combinations. In particular, the extended parameter pace that produces nova eruptions includes low mass transfer rates down to 5e-13 Msun/yr, and more models for low T_{WD}. Resulting characteristics of these runs are added to the former parameter combination results, to provide a full grid spanning the entire parameter space for Carbon-Oxygen white dwarfs. The full grid covers the entire range of observed nova characteristics, even those of peculiar objects, which have not been numerically reproduced until now. Most remarkably, runs for very low Mdot lead to very high values for some characteristics, such as outburst amplitude A>~20, high super-Eddington luminosities at maximum, heavy element abundance of the ejecta Z_{ej}~0.63 and high ejected masses m_{ej}~7e-4 Msun.

We analyze UV observations with FAUST of four sky fields in the general direction of the fourth Galactic quadrant, in which we detect 777 UV sources. This is ~50% more than detected originally by Bowyer et al (1995). We discuss the source detection process and the identification of UV sources with optical counterparts. For the first time in this project we use ground-based objective-prism information for two of the fields, to select the best-matching optical objects with which to identify the UV sources. Using this, and correlations with existing catalogs, we present reliable identifications for \~75% of the sources. Most of the remaining sources have assigned optical counterparts, but lacking additional information we offer only plausible identification. We discuss the types of objects found, and compare the observed population with predictions of our UV Galaxy model.