We describe a new code to simulate the stellar evolution of a close interacting binary system. It is then used to calculate the evolution of a classical nova system composed of a 1.25 Msun Main-Sequence (MS) star and a 1.0 Msun white dwarf (WD) companion. The system begins as a well separated non-interacting binary system. Initially, the two stars evolve independently of each other. However, Roche lobe overflow begins as the MS star expands on its way to become a Red Giant. We follow the mass accreted onto the WD and the ensuing nuclear runaways for several thousand flashes. The main finding is that the Roche-Lobe mass transfer rate is modulated by oscillations in the MS star, with a period that is somewhat shorter than the thermal time scale of the star. This periodically modulates the rate of thermonuclear flashes on the WD, between once every 12000 yrs, such that the WD can cool, to once every 300 yrs, such that it cannot. The system is further complicated by the secular drift in the secondary modulation. Such secondary modulation could explain systems like T Pyxidis. Last, we find that the overall process of mass gain by the WD has an efficiency of roughly 9%, thus requiring a donor with an initial mass of larger than about 5 Msun MS for an initial 1 Msun WD, if the WD is to reach the Chandrasekhar mass.

We study C/O white dwarfs with masses of 1.0 to 1.4 Msun accreting solar-composition material at very high accretion rates. We address the secular changes in the WDs, and in particular, the question whether accretion and the thermonuclear runaways result is net accretion or erosion. The present calculation is unique in that it follows a large number of cycles, thus revealing the secular evolution of the WD system. We find that counter to previous studies, accretion does not give rise to steady state burning. Instead, it produces cyclic thermonuclear runaways of two types. During most of the evolution, many small cycles of hydrogen ignition and burning build a helium layer over the surface of the white dwarf. This He layer gradually thickens and progressively becomes more degenerate. Once a sufficient amount of He has accumulated, several very large helium burning flashes take place and expel the accreted envelope, leaving no net mass accumulation. The results imply that such a system will not undergo an accretion induced collapse, nor will it lead to a SN Type Ia, unless a major new physical process is found.

We generalize the problem of the semi-gray model to cases in which a non-negligible fraction of the stellar radiation falls on the long-wavelength range, and/or that the planetary long-wavelength emission penetrates into the transparent short wavelength domain of the absorption. Second, applying the most general assumptions and independently of any particular properties of an absorber, we show that the greenhouse effect saturates and any Earth-like planet has a maximal temperature which depends on the type of and distance to its main-sequence star, its albedo and the primary atmospheric components which determine the cutoff frequency below which the atmosphere is optically thick. For example, a hypothetical convection-less planet similar to Venus, that is optically thin in the visible, could have at most a surface temperature of 1200-1300K irrespective of the nature of the greenhouse gas. We show that two primary mechanisms are responsible for the saturation of the runaway greenhouse effect, depending on the value of the wavelength above which the atmosphere becomes optically thick. Unless this wavelength is small and resides in the optical region, saturation is achieved by radiating the thermal flux of the planet through the short wavelength tail of the thermal distribution. This has the observational implication, the radiation from such a planet should be skewed towards the NIR. Otherwise, saturation takes place by radiating through windows in the FIR.

Recently Block published an astro-ph{{http://arxiv.org/abs/1106.3928 (2011).}} insinuating that Lemaitre discovery paper of the Expanding Universe was censored prior to its translation into English and publication in the Monthly Notices of the Royal Astronomical Society. Consequently, Lemaitre's credit for the discovery of the velocity-distance correlation was not recognized. We examine here the chain of events leading to the discovery of the 'Hubble law'. Our summary: (a) Lemaitre found a theoretical linear correlation between velocity and distance. (b) Lemaitre assumed the existence of a linear relation between velocity and distance and calculated the coefficient. (c) Hubble took the data plotted it and demonstrated that a linear relation represents the observed data. (d) Hubble never believed in Lemaitre's solution, namely in an expanding universe. Consequently, Hubble never cited Lemaitre. We conclude that the charge that Lemaitre's paper was censored or ignored let alone plagiarized by Hubble, is not founded, and explain why Lemaitre's earlier theoretical discovery and derived 'Hubble constant' was not cited or recognized, by Hubble as well as by many other leading researchers.

Radiation from accretion discs in cataclysmic variable stars (CVs) provides fundamental information about the properties of these close binary systems and about the physics of accretion in general. The detailed diagnostics of accretion disc structure can be achieved by including in its description all the relevant heating and cooling physical mechanism, in particular the convective energy transport that, although dominant at temperatures less than about 10 000 K, is usually not taken into account when calculating spectra of accretion discs. We constructed a radiative transfer code coupled with a code determining the disc's hydrostatic vertical structure. We have obtained for the first time model spectra of cold, convective accretion discs. As expected, these spectra are mostly flat in the optical wavelengths with no contribution from the UV, which in quiescence must be emitted by the white dwarf. The disc structures obtained with our radiative-transfer code compare well with the solutions of equations used to describe the dwarf-nova outburst cycle according to the thermal-viscous disc instability model thus allowing the two to be combined. Our code allows calculating the spectral evolution of dwarf nova stars through their whole outburst cycle, providing a new tool for testing models of accretion discs in cataclysmic variables. We show that convection plays an important role in determining the vertical disc structure and substantially affects emitted spectra when, as often the case, it is effective at optical depths tau ~ 1. The emergent spectrum is independent of the parameters of the convection model.(Abstract shortened)

This paper discusses new perspectives and approaches to the problem of disk dynamics where, in this study, we focus on the effects of viscous instabilities influenced by boundary effects. The Boussinesq approximation of the viscous large shearing box equations is analyzed in which the azimuthal length scale of the disturbance is much larger than the radial and vertical scales. We examine the stability of a non-axisymmetric potential vorticity mode, i.e. a PV-anomaly. in a configuration in which buoyant convection and the strato-rotational instability do not to operate. We consider a series of boundary conditions which show the PV-anomaly to be unstable both on a finite and semi-infinite radial domains. We find these conditions leading to an instability which is the disk analog of Tollmien-Schlichting waves. When the viscosity is weak, evidence of the instability is most pronounced by the emergence of a vortex sheet at the critical layer located away from the boundary where the instability is generated. For some boundary conditions a necessary criterion for the onset of instability for vertical wavelengths that are a sizable fraction of the layer's thickness and when the viscosity is small is that the appropriate Froude number of the flow be greater than one. This instability persists if more realistic boundary conditions are applied, although the criterion on the Froude number is more complicated. The unstable waves studied here share qualitative features to the instability seen in rotating Blasius boundary layers. The implications of these results are discussed. An overall new strategy for exploring and interpreting disk instability mechanisms is also suggested.

We present spectra of accretion discs around white-dwarfs calculated with an improved and updated version of the Shaviv & Wehrse (1991) model. The new version includes line opacities and convective energy transport and can be used to calculate spectra of hot discs in bright systems (nova--like variables or dwarf novae in outburst) as well as spectra of cold accretion discs in quiescent dwarf novae.

The purpose of this paper is to explore how accretion discs manifest the phenomenon of transient growth on a global scale. We investigate analytically the time response of a thin accretion disc to particular axisymmetric perturbations. To facilitate an analytical treatment we replace the energy equation with a general polytropic assumption. The asymptotic expansion of Klu\'zniak & Kita (2000), which extended the method of Regev (1983) to a full steady polytropic disc (with $n=3/2$), is further developed and implemented for both the steady (for any polytropic index) and time-dependent problems. The spatial form and temporal behaviour of selected dynamical disturbances are studied in detail. We identify the perturbation space which leads to transient growth and provide analytical solutions which manifest this expected transient growth behaviour. Three terms (physical causes) responsible for the appearance of transient growth are identified. Two depend explicitly on the viscosity while the third one is relevant also for inviscid discs. The main conclusion we draw is that the phenomenon of transient growth exists in discs on a global scale.

The linear stability of accretion disks is revisited. The governing equations are expanded asymptotically and solved to first order in the expansion parameter $\epsilon$ defined by the ratio of the disk's vertical thickness to its radial extent. Algebraically growing solutions are found for global perturbations on the radial accretion flow of thin inviscid compressible Keplerian disks. The algebraic temporal behavior is exhibited in the vertical velocities and the thermodynamic variables and has the form $t\sin\Omega_0 t$ locally in the disk where $\Omega_0$ is the Keplerian rotation rate. The physical implications and relations to the Solberg-Hoiland stability criteria are discussed.

The exact ionization state of Be7 in the solar core is crucial for the precise prediction of the solar B8 neutrino flux. We therefore examine the effect of pressure ionization on the ionization state of Be7 and all elements with 12 >= Z >= 4. We show that under the conditions prevailing in the solar core, one has to consider the effect of the nearest neighbor on the electronic structure of a given ion. To this goal, we first solve the Schroedinger and then the Kohn-Sham equations for an ion immersed in a dense plasma under conditions for which the mean interparticle distance is smaller than the Debye radius. The question of which boundary conditions should be imposed on the wave function is discussed, examined and found to be crucial. Contrary to previous estimates showing that Beryllium is partially ionized, we find that it is fully ionized. As a consequence, the predicted rate of the Be7 + e- reaction is reduced by 20-30%, depending on the exact solar model. Since Be7 is a trace element, its total production is controlled by the unchanged He4+He3 reaction rate, and its destruction is determined by the rate of electron capture. As the latter rate decreases when the Beryllium is fully ionized (relative to the case of partially ionized Be), the estimate for the abundance of Be7 increases and with it the B8 neutrino flux. The increase in phi_nu(B8) is by about 20-30%. The neutrino flux due to Be7 electron capture remains effectively unchanged because the change in the rate is compensated for by a change in the abundance. Hence the prediction for the ratio of phi_nu(B8) / phi_nu(Be7) changes as well.

In a recent paper, Bahcall et al. (2000) list various new approaches to the problem of screening of nuclear reactions in stellar plasma and assert that they are all wrong or irrelevant. Except for two, all approaches mentioned by Bahcall et al. assume the mean field approximation. The two exceptions are Carraro et al. (1988) and Shaviv & Shaviv (2000a). While Carraro et al. (1988) paper is discussed shortly and refuted by Bahcall et al. (2000) the Shaviv & Shaviv (2000a) paper is not discussed and refuted only by association. However, the association is totally unfounded because Shaviv & Shaviv (2000a) have shown that kinetic equations must be used to solve the screening problem and that the mean field approximation is inadequate for this problem. They also showed that the Carraro et al. (1988) approach is erroneous. Therefore we summarize here the method of $S^{2}$ and their main result. We contrast the kinetic equations method with the mean field approximation and expose the different assumptions and omissions in each method.

The $^8$B solar neutrino flux as measured by Super-Kamiokande is consistent with the $^{37}$Ar production rate in $^{37}$Cl at Homestake. GALLEX and SAGE, continue to observe $^{71}$Ge production rates in $^{71}$Ga that are consistent with the minimal signal expected from the solar luminosity. The observed $^8$B solar neutrino flux is in good agreement with that predicted by the standard solar model of Dar and Shaviv with nuclear reaction rates that are supported by recent measurements of nuclear fusion cross sections at low energies. The measurements of Super-Kamiokande, SAGE and GALLEX suggest that the expected the pep, $^7$Be and NO solar neutrino fluxes are strongly suppressed. This can be explained by neutrino oscillations and the Mikheyev-Smirnov-Wolfenstein effect. Since neither a flavor change, nor a terrestrial variation, nor a spectral distortion of the $^8$B solar neutrino flux has been observed yet, the solar neutrino problem does not provide conclusive evidence for neutrino properties beyond the standard electroweak model. The deviations of the experimental results from those predicted by the standard solar models may reflect the approximate nature of of solar models and of our knowledge of nuclear reaction rates, radiation transport and particle diffusion in dense stellar plasmas. Only future observations of spectral distortions, or terrestrial modulation or flavor change of solar neutrinos in solar neutrino experiments, such as Super-Kamiokande, SNO, Borexino and HELLAZ will be able to establish that neutrino properties beyond the minimal standard electroweak model are responsible for the solar neutrino problem.

We use the Shaviv and Wehrse (1991) code to model the vertical structure and the emission properties of quiescent dwarf nova discs. We find that in the case of HT Cas the quiescent disc must be optically thin, in contradiction with the requirements of the standard disc instability model. We find a viscosity parameter alpha > 1. Although this is much less than values (~ 100) obtained in isothermal slab models it is not consistent with the accretion disc model assumptions.

An improved standard solar model has been used to calculate the fluxes of standard solar neutrinos. It includes premain sequence evolution, element diffusion, partial ionization effects, and all the possible nuclear reactions between the main elements. It uses updated values for the initial solar element abundances, the solar age, the solar luminosity, the nuclear reaction rates and the radiative opacities. Neither nuclear equilibrium, nor complete ionization are assumed. The calculated solar neutrino fluxes are compared with published results from the four solar neutrino experiments. The calculated $^8$B solar neutrino flux is consistent, within the theoretical and experimental uncertainties, with the solar neutrino observations at Homestake and Kamiokande. The observations suggest that the $^7$Be solar neutrino flux is much smaller than that predicted. However, conclusive evidence for the suppression of the $^7$Be solar neutrino flux will require experiments like BOREXINO and HELLAZ. If the $^7$Be solar neutrino flux is suppressed, it still can be due either to standard physics and astrophysics or neutrino properties beyond the standard electroweak model. Only future neutrino experiments, such as SNO, Superkamiokande, BOREXINO and HELLAZ, will be able to show that the solar neutrino problem is a consequence of neutrino properties beyond the standard electroweak model.

In a recent paper Dar and Shaviv (1994) presented results of calculations with an improved standard solar model (SSM) that suggest a standard physics solution to the solar neutrino problem. In a subsequent publication with a similar title Bahcall et al. (1994) have claimed to refute the results obtained by Dar and Shaviv. Here we show that the criticism and conclusions of Bahcall et al. are unjustified. In particular, their attempts to reproduce the results of Dar and Shaviv failed because they did not use the same nuclear cross sections and did not include in their calculations major physics improvements which were included in the improved SSM code used by Dar and Shaviv.

An imroved standard solar model, with more accurate input parameters and more accurate treatment of plasma physics effects, predicts a solar neutrino flux which is consistent within experimental uncertainties with the solar neutrino observations after 1986.