We explore the appearance of light clusters at high densities of collapsing stellar cores. Special attention is paid to the unstable isotope H4, which was not included in previous studies. The importance of light clusters in the calculation of rates for neutrino matter interaction is discussed. The main conclusion is that thermodynamic quantities are only weakly sensitive to the chemical composition. The change in pressure and hence the direct change in collapse dynamics will be minor. But the change in neutrino heating and neutronization processes can be significant.

The general solution of hydrostatic equilibrium equations for a two-component fluid of ions and electrons without a local electroneutrality constraint is found in the framework of Newtonian gravity theory. In agreement with the Poincar\'e theorem on analyticity and in the context of Dyson's argument, the general solution is demonstrated to possess a fixed (essential) singularity in the gravitational constant $G$ at $ G = 0 $. The regular component of the general solution can be determined by perturbation theory in $G$ starting from a locally neutral solution. The non-perturbative component obtained using the method of Wentzel, Kramers and Brillouin is exponentially small in the inner layers of the star and grows rapidly in the outward direction. Near the surface of the star, both components are comparable in magnitude, and their non-linear interplay determines the properties of an electro- or ionosphere. The stellar charge varies within the limits of $- 0.1 $ to $150$ C per solar mass. The properties of electro- and ionospheres are exponentially sensitive to variations of the fluid densities in the central regions of the star. The general solutions of two exactly solvable stellar models without a local electroneutrality constraint are also presented.

The variational principle for stars with a phase transition has been investigated. The term outside the integral in the expression for the second variation of the total energy of a star is shown to be obtained by passage to the limit from the integration over the region of mixed states in the star. The form of the trial functions ensuring this passage has been found. All of the results have been generalized to the case where general relativity is applicable. The known criteria for the dynamical stability of a star when a new phase appears at its center are shown to follow automatically from the variational principle. Numerical calculations of hydrostatically equilibrium models for hybrid stars with a phase transition have been performed. The form of the trial functions for the second variation of the total energy of a star that describes almost exactly the stability boundaries of such stellar models is proposed.

The characteristics of the gravitational collapse of a supernova and the fluxes of active and sterile neutrinos produced during the formation of its protoneutron core have been calculated numerically. The relative yields of active and sterile neutrinos in core matter with different degrees of neutronization have been calculated for various input parameters and various initial conditions. A significant increase in the fraction of sterile neutrinos produced in superdense core matter at the resonant degree of neutronization has been confirmed. The contributions of sterile neutrinos to the collapse dynamics and the total flux of neutrinos produced during collapse have been shown to be relatively small. The total luminosity of sterile neutrinos is considerably lower than the luminosity of electron neutrinos, but their spectrum is considerably harder at high energies.

It is shown that the unusual thermodynamic properties of matter within the region of two-phase coexistence in hybrid stars result in a change of the standard condition for beginning of convection. In particular, the thermal flux transported by convection may be directed towards the stellar center. We discuss favorable circumstances leading to such an effect of "inverse convection" and its possible influence on the thermal evolution of hybrid stars.

The relative yields of active and sterile neutrinos in the matter with a high density and different degree of neutronization are calculated. A significant increase in the proportion of sterile neutrinos produced in superdense matter when degree of neutronization approaches the value of two is found. The results obtained can be used in the calculations of the neutrino fluxes for media with a high density and different neutronization degrees in astrophysical processes such as the formation of protoneutron core of a supernova.

An analytical study that explains the existence of a very small region on the mass radius diagram of hybrid stars where all of the lines representing the sequences of models with different constant values of the bag constant B intersect is presented. This circumstance is shown to be a consequence of the linear dependence of pressure on energy density in the quark cores of hybrid stars.

The nucleosynthesis yields from neutrino-induced spallation of 4He are calculated in case of penetration of some 4He into the carbon-oxygen-neon shell. The neutrino contribution to creation of the weak r-process component becomes noticeable in case 4He is dragged down at radii R < 10^9 cm.

Principally important for the description the physical processes in the collapsing stellar cores topics are surveyed. They are: the neutrino heat conduction theory, equation of state under the conditions of nuclear statistical equilibrium and possible phase transitions in dense, subnuclear and nuclear environment.

Some features of the physics of radiation-dominated shock waves are discussed with emphasis on the peculiarities which are important for correct numerical modeling of shock breakouts in supernova. With account of those peculiarities, a number of models for different supernova types is constructed based on multigroup radiation transfer coupled to hydrodynamics. We describe the implementation of a new algorithm RADA, designed for modeling photon transfer at extremely-relativistic motions of matter, into our older code STELLA. The results of numerical simulations of light curves, and continuum spectra are presented. The influence of effects of photon scattering on electrons, of thermalization depth and of special relativity in transfer equation is considered. Some cases are demonstrated, when the appearance of hard X-ray emission is possible at the shock breakout. The necessary refinements in numerical algorithms for radiative transfer and hydrodynamics are pointed out. Prospects for using the results of numerical simulation to analyze and interpret available and future data from space observatories are discussed.

Phase transition from hadronic matter to quark-gluon matter is discussed for various regimes of temperature and baryon number density. For small and medium densities, the phase transition is accurately described in the framework of the Field Correlation Method, whereas at high density predictions are less certain and leave room for the phenomenological models. We study formation of multiquark states (MQS) at zero temperature and high density. Relevant MQS components of the nuclear matter can be described using a previously developed formalism of the quark compound bags (QCB). Partial-wave analysis of nucleon-nucleon scattering indicates the existence of 6QS which manifest themselves as poles of $P$-matrix. In the framework of the QCB model, we formulate a self-consistent system of coupled equations for the nucleon and 6QS propagators in nuclear matter and the G-matrix. The approach provides a link between high-density nuclear matter with the MQS components and the cumulative effect observed in reactions on the nuclei, which requires the admixture of MQS in the wave functions of nuclei kinematically. 6QS determine the natural scale of the density for a possible phase transition into the MQS phase of nuclear matter. Such a phase transition can lead to dynamic instability of newly born protoneutron stars and dramatically affect the dynamics of supernovae. Numerical simulations show that the phase transition may be a good remedy for the triggering supernova explosions in the spherically symmetric supernova models. A specific signature of the phase transition is an additional neutrino peak in the neutrino light curve. For a Galactic core-collapse supernova, such a peak could be resolved by the present neutrino detectors. The possibility of extracting the parameters of the phase of transition from observation of the neutrino signal is discussed also.

The main observational properties and resulting classification of supernovae (SNe) are briefly reviewed. Then we discuss the progress in modeling of two basic types of SNe - the thermonuclear and core-collapse ones, with special emphasis being placed on difficulties relating to a consistent description of thermonuclear flame propagation and the detachment of supernova envelope from the collapsing core (a nascent neutron star). The properties of the neutrino flux expected from the corecollapse SNe, and the lessons of SN1987A, exploded in the Large Magellanic Cloud, are considered as well.

The similarity theory of stellar structure is used to study the properties of very massive stars when one can neglect all the sources of opacity except for the Thomson scattering. The dimensionless internal structure of such stars is practically independent of the energy generation law. It is shown that the mass-luminosity relation can be approximated by an analytical expression which is virtually universal regarding the chemical composition and energy generation law. A detailed comparison with the Eddington standard model is given. The application of the obtained results to the observations of massive stars is briefly discussed.

Global properties of different types of supernovae are overviewed. A special emphasis is made on the issues of interest both for the supernova and cosmic ray physics such as supernova shock waves and neutrino signal.

We analyze the behavior of the outer envelope in a massive star during and after the collapse of its iron core into a protoneutron star (PNS) in terms of the equations of one-dimensional spherically symmetric ideal hydrodynamics. The profiles obtained in the studies of the evolution of massive stars up to the final stages of their existence, immediately before a supernova explosion (Boyes et al. 1999), are used as the initial data for the distribution of thermodynamic quantities in the envelope.We use a complex equation of state for matter with allowances made for arbitrary electron degeneracy and relativity, the appearance of electron-positron pairs, the presence of radiation, and the possibility of iron nuclei dissociating into free nucleons and helium nuclei. We performed calculations with the help of a numerical scheme based on Godunov's method. These calculations allowed us to ascertain whether the emersion of the outer envelope in a massive star is possible through the following two mechanisms: first, the decrease in the gravitational mass of the central PNS through neutrino-signal emission and, second, the effect of hot nucleon bubbles, which are most likely formed in the PNS corona, on the envelope emersion. We show that the second mechanism is highly efficient in the range of acceptable masses of the nucleon bubbles ($\leq 0.01M_\odot$) simulated in our hydrodynamic calculations in a rough, spherically symmetric approximation.

The hydrodynamical modelling of Type II plateau supernova light curves predicts a correlation between three observable parameters (the plateau duration, the absolute magnitude and photospheric velocity at the middle of the plateau) on the one side and three physical parameters (the explosion energy E, the mass of the envelope expelled M, and the presupernova radius R) on the other side. The correlation is used, together with adopted EPM distances, to estimate E, M, and R for a dozen of well-observed SNe IIP. For this set of supernovae, the resulting value of E varies within a factor of 6 [0.5< E/(10^51 erg)<3], whereas the envelope mass remains within the limits 10< M/M_sun <30. The presupernova radius is typically (200-600)R_sun, but can exceed 1000R_sun for the brightest supernovae (e.g., SN 1992am). A new method of determining the distance of SNe IIP is proposed. It is based on the assumption of a correlation between the explosion energy E and the Ni56 mass required to power the post-plateau light curve tail through Co56 decay. The method is useful for SNe IIP with well-observed bolometric light curves both during the plateau and radioactive tail phases. The resulting distances and future improvements are discussed.

We simulate shock propagation in a wide range of core-collapsing presupernovae: from compact WR stars exploding as SNe Ib/c through very extended envelopes of the narrow-line SNe IIn. We find that the same physical phenomenon of radiating shocks can produce outbursts of X-ray radiation (with photon energy $3kT \sim 1$ keV) lasting only a second in SNe Ib/c, as well as a very high flux of visual light, lasting for months, in SNe IIn.

Numerical modeling of shock breakout in compact presupernovae of type Ib/c is done by two independent methods. Peak color temperatures reach 4 to 8 million K for explosion energies 1 to 9 foe with a typical local time-scale of 0.03 seconds. In the obverver's frame the burst is smeared over few seconds due to light travel time correction R/c.

The huge neutrino pulse that occurs during the collapse of a massive stellar core, is expected to contribute to the origination of a number of isotopes both of light chemical elements and heavy ones. It is shown that, in general, the heating of stellar matter due to the neutrino scattering off electrons and the heat released from the neutrino-helium breakup followed by the thermonuclear reactions should be taken into account. On the base of kinetic network, using all the important reactions up to Z=8, the main features and the time-dependent character of the neutrino- driven neutron flux are investigated. The time-dependent densities of free neutrons produced in helium breakup, Y_n(t), were used to calculate the r-process nucleosynthesis with another full kinetic network for 3200 nuclides. It was found that in the case of metal-deficient stars, Z < 0.01 Z(solar), the resulting density of free neutrons seems to be high enough to drive the r-process efficiently under favorable conditions. But it is impossible to obtain a sufficient amount of heavy nuclei in neutrino-induced r-process in a helium shell at radii R > R_cr \approx 10^9 cm. We speculate that to make the neutrino-induced r-process work efficiently in the shell, one has to invoke nonstandard presupernova models in which helium hopefully is closer to the collapsed core owing, for instance, to a large scale mixing or/and rotation and magnetic fields. Apart from this exotic possibility, the neutrino-induced nucleosynthesis in the helium shell is certainly not strong enough to explain the observed solar r-process abundances.