The enigmatic type IIP SN 2016X demonstrates the unprecedented asphericity in the nebular H-alpha line profile, the absence of nebular [O I] emission, and the unusual occultation effect due to the internal dust. The hydrodynamic modelling of the bolometric light curve and expansion velocities suggests that the event is an outcome of the massive star explosion that ejected 28 Msun with the kinetic energy of 1.7x10^51 erg and 0.03 Msun of radioactive Ni-56. We recover the bipolar distribution of Ni-56 from the H-alpha profile via the simulation of the emissivity produced by non-spherical Ni-56 ejecta. The conspicuous effect of the dust absorption in the H-alpha profile rules out the occultation by the dusty sphere or dusty thick disk but turns out consistent with the thin dusty disk-like structure in the plane perpendicular to the bipolar axis. We speculate that the absence of the nebular [O I] emission might originate from the significant cooling of the oxygen-rich matter mediated by CO and SiO molecules.

With the same method as used previously, we investigate neutrino-driven explosions of a larger sample of blue supergiant models. The larger sample includes three new presupernova stars. The results are compared with light-curve observations of the peculiar type IIP SN 1987A. The explosions were modeled in 3D with the neutrino-hydrodynamics code PROMETHEUS-HOTB, and light-curve calculations were performed in spherical symmetry with the radiation-hydrodynamics code CRAB. Our results confirm the basic findings of the previous work: 3D neutrino-driven explosions with SN 1987A-like energies synthesize an amount of Ni-56 that is consistent with the radioactive tail of the light curve. Moreover, the models mix hydrogen inward to minimum velocities below 400 km/s as required by spectral observations. Hydrodynamic simulations with the new progenitor models, which possess smaller radii than the older ones, show much better agreement between calculated and observed light curves in the initial luminosity peak and during the first 20 days. A set of explosions with similar energies demonstrated that a high growth factor of Rayleigh-Taylor instabilities at the (C+O)/He composition interface combined with a weak interaction of fast Rayleigh-Taylor plumes, where the reverse shock occurs below the He/H interface, provides a sufficient condition for efficient outward mixing of Ni-56 into the hydrogen envelope. This condition is realized to the required extent only in one of the older stellar models, which yielded a maximum velocity of around 3000 km/s for the bulk of ejected Ni-56, but failed to reproduce the helium-core mass of 6 Msun inferred from the absolute luminosity of the presupernova star. We conclude that none of the single-star progenitor models proposed for SN 1987A to date satisfies all constraints set by observations. (Abridged)

Based on the work of Menon & Heger (2017), we present the bolometric light curvesand spectra of the explosions of blue supergiant progenitors from binary mergers. We study SN 1987A and two other peculiar Type IIP supernovae: SN 1998A and SN 2006V. The progenitor models were produced using the stellar evolution code Keplerand then exploded using the 1D radiation hydrodynamic code Crab. The explosions of binary merger models exhibit an overall better fit to the light curve of SN 1987A than previous single star models, because of their lower helium-core masses, larger envelope masses, and smaller radii. The merger model that best matches the observational constraints of the progenitor of SN 1987A and the light curve is a model with a radius of 37 solar radii, an ejecta mass of 20.6 solar masses, an explosion energy of 1.7 Bethe, a nickel mass of 0.073 solar masses, and a nickel mixing velocity of 3,000 km/s. This model also works for SN 1998A and is comparable with earlier estimates from semi-analytic models. In the case of SN 2006V, however, a model with a radius of 150 solar radii and ejecta mass of 19.1 solar masses matches the light curve. These parameters are significantly higher than predictions from semi-analytic models for the progenitor of this supernova.

We explore the well-observed type IIP SN 2013ej with peculiar luminosity evolution. It is found that the hydrodynamic model cannot reproduce in detail the bolometric luminosity at both the plateau and the radioactive tail. Yet the ejecta mass of 23-26 Msun and the kinetic energy of (1.2-1.4)x10^51 erg are determined rather confidently. We suggest that the controversy revealed in hydrodynamic simulations stems from the strong asphericity of the Ni-56 ejecta. An analysis of the asymmetric nebular H-alpha line and of the peculiar radioactive tail made it possible to recover parameters of the asymmetric bipolar Ni-56 ejecta with the heavier jet residing in the rear hemisphere. The inferred Ni-56 mass is 0.039 Msun, twice as large compared to a straightforward estimate from the bolometric luminosity at the early radioactive tail. The bulk of ejected Ni-56 has velocities in the range of 4000-6500 km/s. The linear polarization predicted by the model with the asymmetric ionization produced by bipolar Ni-56 ejecta is consistent with the observational value.

Type II-plateau supernovae (SNe IIP) are the most numerous subclass of core-collapse SNe originating from massive stars. In the framework of the neutrino-driven explosion mechanism, we study the SN outburst properties for a red supergiant progenitor model and compare the corresponding light curves with observations of the ordinary Type IIP SN 1999em. Three-dimensional (3D) simulations of (parametrically triggered) neutrino-driven explosions are performed with the (explicit, finite-volume, Eulerian, multifluid hydrodynamics) code PROMETHEUS, using a presupernova model of a 15 Msun star as initial data. At approaching homologous expansion, the hydrodynamical and composition variables of the 3D models are mapped to a spherically symmetric configuration, and the simulations are continued with the (implicit, Lagrangian radiation-hydrodynamics) code CRAB to follow the blast-wave evolution during the SN outburst. Our 3D neutrino-driven explosion model with an explosion energy of about 0.5x10^51 erg produces Ni-56 in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail of SN 1999em. The considered presupernova model, 3D explosion simulations, and light-curve calculations can explain the basic observational features of SN 1999em, except for those connected to the presupernova structure of the outer stellar layers. Our 3D simulations show that the distribution of Ni-rich matter in velocity space is asymmetric with a strong dipole component that is consistent with the observations of SN 1999em. The monotonic luminosity decline from the plateau to the radioactive tail in ordinary SNe IIP is a manifestation of the intense turbulent mixing at the He/H composition interface.

The importance of allowance for the time-dependent effect in the kinetics at the photospheric phase during a supernova explosion has been confirmed by several independent research groups. The time-dependent effect provides a higher degree of hydrogen ionization in comparison with the steady-state solutions and strengthens the H$\alpha$ line in the resulting simulated spectrum, with the intensity of the effect increasing with time. However, some researchers argue that the time-dependent ionization effect is unimportant. Its allowance leads to an insignificant strengthening of H$\alpha$ in their modeling only in the first days after explosion. We have demonstrated the importance of the time-dependent effect with the models of SN 1999em as an example using the new original LEVELS software package. The role of a number of factors that can weaken the time-dependent effect has been checked. We have confirmed that the intensity of the effect is affected by the abundance of metal admixtures in the envelope, while the addition of extra levels to themodel hydrogen atom weakens the time-dependent effect to a lesser degree and never removes it completely.

We present the results of 3D-hydrodynamical simulations of accretion flow in the eclipsing dwarf nova V1239 Her in quiescence. The model includes the optical star filling its Roche lobe, a gas stream emanating from the inner Lagrangian point of the binary system, and the accretion disc structure. A cold hydrogen gas stream is initially emitted towards a point-like gravitational centre. A stationary accretion disc is formed in about 15 orbital periods after the beginning of accretion. The model takes into account partial ionization of hydrogen and uses realistic cooling function for hydrogen. The light curve of the system is calculated as the volume emission of optically thin layers along the line of sight up to the optical depth \tau =2/3 calculated using Planck-averaged opacities. The calculated eclipse light curves show good agreement with observations, with the changing shape of pre-eclipse and post-eclipse light curves being explained entirely due to ~ 50% variations in the mass accretion rate through the gas stream.

Spectrophotometry of SN 1996al carried out throughout 15 years is presented. The early photometry suggests that SN 1996al is a Linear type-II supernova, with an absolute peak of Mv ~ -18.2 mag. Early spectra present broad, asymmetric Balmer emissions, with super-imposed narrow lines with P-Cygni profile, and He I features with asymmetric, broad emission components. The analysis of the line profiles shows that the H and He broad components form in the same region of the ejecta. By day +142, the Halpha profile dramatically changes: the narrow P-Cygni profile disappears, and the Halpha is fitted by three emission components, that will be detected over the remaining 15 yrs of the SN monitoring campaign. Instead, the He I emissions become progressively narrower and symmetric. A sudden increase in flux of all He I lines is observed between 300 and 600 days. Models show that the supernova luminosity is sustained by the interaction of low mass (~1.15 Msun) ejecta, expelled in a low kinetic energy (~ 1.6 x 10^50 erg) explosion, with highly asymmetric circumstellar medium. The detection of Halpha emission in pre-explosion archive images suggests that the progenitor was most likely a massive star (~25 Msun ZAMS) that had lost a large fraction of its hydrogen envelope before explosion, and was hence embedded in a H-rich cocoon. The low-mass ejecta and modest kinetic energy of the explosion are explained with massive fallback of material into the compact remnant, a 7-8 Msun black hole.

The well-studied type IIP SN 1987A, produced by the explosion of a blue supergiant (BSG) star, is a touchstone for massive-star evolution, simulations of neutrino-driven explosions, and modeling of light curves and spectra. In the framework of the neutrino-driven mechanism, we study the dependence of explosion properties on the structure of four different BSGs and compare the corresponding light curves with observations of SN 1987A. We perform 3D simulations with the PROMETHEUS code until about one day and map the results to the 1D code CRAB for the light curve calculations. All of our 3D models with explosion energies compatible with SN 1987A produce 56Ni in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail. One of the progenitors yields maximum velocities of ~3000 km/s for the bulk of ejected 56Ni, consistent with observations. In all of our models inward mixing of hydrogen during the 3D evolution leads to minimum H-velocities below 100 km/s, in good agreement with spectral observations. The considered BSG models, 3D explosion simulations, and light-curve calculations can thus explain basic observational features of SN 1987A. However, all progenitors have too large pre-SN radii to reproduce the narrow initial luminosity peak, and the structure of their outer layers is not suitable to match the observed light curve during the first 30-40 days. Only one stellar model has a structure of the He core and the He/H composition interface that enables sufficient outward mixing of 56Ni and inward mixing of hydrogen to produce a good match of the dome-like shape of the observed light-curve maximum. But this model falls short of the He-core mass of 6 Msun inferred from the absolute luminosity of the pre-SN star. The lack of an adequate pre-SN model for SN 1987A is a pressing challenge for the theory of massive-star evolution. (Abridged)

The explosion energy and the ejecta mass of a type IIP supernova (SN IIP) derived from hydrodynamic simulations are principal parameters of the explosion theory. However, the number of SNe IIP studied by hydrodynamic modeling is small. Moreover, some doubts exist in regard to the reliability of derived SN IIP parameters. The well-observed type IIP SN 2012A will be studied via hydrodynamic modeling. Their early spectra will be checked for a presence of the ejecta clumpiness. Other observational effects of clumpiness will be explored. Supernova parameters are determined by means of the standard hydrodynamic modeling. The early hydrogen Halpha and Hbeta lines are used for the clumpiness diagnostics. The modified hydrodynamic code is employed to study the clumpiness effect in the light curve and expansion kinematics. We found that SN 20012A is the result of the explosion of a red supergiant with the radius of 715 Rsun. The explosion energy is 5.25x10^50 erg, the ejecta mass is 13.1 Msun, and the total Ni-56 mass is 0.012 Msun. The estimated mass of a progenitor, a main-sequence star, is 15 Msun. The Halpha and Hbeta lines in early spectra indicate that outer ejecta are clumpy. Hydrodynamic simulations show that the clumpiness modifies the early light curve and increases the maximum velocity of the outer layers. The pre-SN 2012A was a normal red supergiant with the progenitor mass of about 15 Msun. The outer layers of ejecta indicate the clumpy structure. The clumpiness of the external layers can increase the maximum expansion velocity.

We determine the oxygen density in the central zone of nine type IIP supernovae (SN~IIP) at the nebular stage using oxygen doublet [O I] 6300, 6364 \AA. Combined with two available estimates these data indicate that oxygen densities on day 300 are distributed in rather narrow range $(2.3\pm1)\times10^9$ cm$^{-3}$. The result does not depend on the distance, extinction, or model assumptions. We demonstrate that the found density distribution suggests that the explosion energy of SN~IIP increases with the stellar mass.

We study disparity between H\alpha\ and H\beta\ in early spectra of the type IIP supernova SN 2008in. The point is that these lines cannot be described simultaneously in a spherically-symmetric model with the smooth density distribution. It is shown that an assumption of a clumpy structure of external layers of the envelope resolves the problem. We obtain estimates of the velocity at the inner border of the inhomogeneous zone (\approx6100 km s^{-1}), the filing factor of inhomogeneities (\leq0.5), and the mass of the inhomogeneous layers (~ 0.03 M_{\odot}). The amplitude of flux fluctuations in the early spectrum of H\alpha\ (\Delta F/F ~ 10^{-2}) imposes a constraint on the size of inhomogeneities (\leq 200 km s^{-1}). A detection of fluctuations in the early H\alpha\ of type IIP supernovae might become an observational test of the inhomogeneous structure of their envelopes. We propose also the indirect test of the clumpy structure of external layers: the study of properties of the initial radiation outburst due to the shock breakout. The inhomogeneous structure of external layers of type IIP supernovae could be an outcome of density perturbations and density inversion in outer convective layers of presupernova red supergiant.

The explosion energy and the ejecta mass of a type IIP supernova make up the basis for the theory of explosion mechanism. So far, these parameters have only been determined for seven events. Type IIP supernova 2008in is another well-observed event for which a detailed hydrodynamic modeling can be used to derive the supernova parameters. Hydrodynamic modeling was employed to describe the bolometric light curve and the expansion velocities at the photosphere level. A time-dependent model for hydrogen ionization and excitation was applied to model the Halpha and Hbeta line profiles. We found an ejecta mass of 13.6 Msun, an explosion energy of 5.05x10^50 erg, a presupernova radius of 570 Rsun, and a radioactive Ni-56 mass of 0.015 Msun. The estimated progenitor mass is 15.5 Msun. We uncovered a problem of the Halpha and Hbeta description at the early phase, which cannot be resolved within a spherically symmetric model. The presupernova of SN 2008in was a normal red supergiant with the minimum mass of the progenitor among eight type IIP supernovae explored by means of the hydrodynamic modeling. The problem of the absence of type IIP supernovae with the progenitor masses <15 Msun in this sample remains open.

Among type IIP supernovae there are a few events that resemble the well-studied supernova 1987A produced by the blue supergiant in the Large Magellanic Cloud. We study a peculiar supernova 2000cb and compare it with the supernova 1987A. We carried out hydrodynamic simulations of the supernova in an extended parameter space to describe its light curve and spectroscopic data. The hydrogen H-alpha and H-beta lines are modeled using a time-dependent approach. We constructed the hydrodynamic model by fitting the photometric and spectroscopic observations. We infer a presupernova radius of 35 Rsun, an ejecta mass of 22.3 Msun, an explosion energy of 4.4x10^{51} erg, and a radioactive Ni-56 mass of 0.083 Msun. The estimated progenitor mass on the main sequence lies in the range of 24-28 Msun. The early H-alpha profile on day 7 is consistent with the density distribution found from hydrodynamic modeling, while the H-alpha line on day 40 indicates an extended Ni-56 mixing up to a velocity of 8400 km/s. We emphasize that the dome-like light curves of both supernova 2000cb and supernova 1987A are entirely powered by radioactive decay. This is unlike normal type IIP supernovae, the plateau of which is dominated by the internal energy deposited after the shock wave propagation through the presupernova. We find signatures of the explosion asymmetry in the photospheric and nebular spectra. The explosion energy of supernova 2000cb is higher by a factor of three compared to supernova 1987A, which poses a serious problem for explosion mechanisms of type IIP supernovae.

Unusually bright type IIP supernova (SN) 2009kf is studied employing the hydrodynamic modelling. We derived optimal values of the ejecta mass of 28.1 Msun, explosion energy of 2.2x10^{52} erg, and presupernova radius of 2x10^3 Rsun assuming that Ni-56 mass is equal to the upper limit of 0.4 Msun. We analyzed effects of the uncertainties in the extinction and Ni-56 mass and concluded that both the ejecta mass and explosion energy cannot be significantly reduced compared with the optimal values. The huge explosion energy of SN 2009kf indicates that the explosion is caused by the same mechanism which operates in energetic SNe Ibc (hypernovae), i.e., via a rapid disk accretion onto black hole. The ejecta mass combined with the black hole mass and the mass lost by stellar wind yields the progenitor mass of about 36 Msun. We propose a scenario in which massive binary evolution might result in the SN 2009kf event.

Previous studies of type IIP supernovae have inferred that progenitor masses recovered from hydrodynamic models are higher than 15 Msun. To verify the progenitor mass of this supernova category, we attempt a parameter determination of the well-observed luminous type IIP supernova 2004et. We model the bolometric light curve and the photospheric velocities of SN 2004et by means of hydrodynamic simulations in an extended parameter space. From hydrodynamic simulations and observational data, we infer a presupernova radius of 1500 Rsun, an ejecta mass of 24.5 Msun, an explosion energy of 2.3x10^{51} erg, and a radioactive Ni-56 mass of 0.068 Msun. The estimated progenitor mass on the main sequence is in the range of 25-29 Msun. In addition, we find clear signatures of the explosion asymmetry in the nebular spectra of SN 2004et. The measured progenitor mass of SN 2004et is significantly higher than the progenitor mass suggested by the pre-explosion images. We speculate that the mass inferred from hydrodynamic modeling is overestimated and crucial missing factors are multi-dimensional effects.

The progenitor mass of type IIP supernova can be determined from either hydrodynamic modeling of the event or pre-explosion observations. To compare these approaches, we determine parameters of the sub-luminous supernova 2005cs and estimate its progenitor mass. We compute the hydrodynamic models of the supernova to describe its light curves and expansion velocity data. We estimate a presupernova mass of 17.3 Msun, an explosion energy of 4.1x10^{50} erg, a presupernova radius of 600 Rsun, and a radioactive Ni-56 mass of 0.0082 Msun. The derived progenitor mass of SN 2005cs is 18.2 Msun, which is in-between those of low-luminosity and normal type IIP supernovae. The obtained progenitor mass of SN 2005cs is higher than derived from pre-explosion images. The masses of four type IIP supernovae estimated by means of hydrodynamic modeling are systematically higher than the average progenitor mass for the 9-25 Msun mass range. This result, if confirmed for a larger sample, would imply that a serious revision of the present-day view on the progenitors of type IIP supernovae is required.

We study a possibility of detection of circumstellar absorption lines of Na I D$_{1,2}$ and Ca II H,K in spectra of type IIP supernovae at the photospheric epoch. The modelling shows that the circumstellar lines of Na I doublet will not be seen in type IIP supernovae for moderate wind density, e.g., characteristic of SN 1999em, whereas rather pronounced Ca II lines with P Cygni profile should be detectable. A similar model is used to describe Na I and Ca II circumstellar lines seen in SN 1998S, type IIL with a dense wind. We show that line intensities in this supernova are reproduced, if one assumes an ultraviolet excess, which is caused primarily by the comptonization of supernova radiation in the shock wave.

The origin of low-luminosity Type IIP supernovae is unclear: they have been proposed to originate either from massive (about 25 Msun) or low-mass (about 9 Msun) stars. We wish to determine parameters of the low-luminosity Type IIP supernova 2003Z, to estimate a mass-loss rate of the presupernova, and to recover a progenitor mass. We compute the hydrodynamic models of the supernova to describe the light curves and the observed expansion velocities. The wind density of the presupernova is estimated using a thin shell model for the interaction with circumstellar matter. We estimate an ejecta mass of 14 Msun, an explosion energy of 2.45x10^50 erg, a presupernova radius of 229 Rsun, and a radioactive Ni-56 amount of 0.0063 Msun. The upper limit of the wind density parameter in the presupernova vicinity is 10^13 g/cm, and the mass lost at the red/yellow supergiant stage is less than 0.6 Msun assuming the constant mass-loss rate. The estimated progenitor mass is in the range of 14.4-17.4 Msun. The presupernova of SN 2003Z was probably a yellow supergiant at the time of the explosion. The progenitor mass of SN 2003Z is lower than those of SN 1987A and SN 1999em, normal Type IIP supernovae, but higher than the lower limit of stars undergoing a core collapse. We propose an observational test based on the circumstellar interaction to discriminate between the massive (about 25 Msun) and moderate-mass (about 16 Msun) scenarios.

We propose new diagnostics for circumstellar interaction in Type IIP supernovae by the detection of high velocity (HV) absorption features in Halpha and He I 10830 A lines during the photospheric stage. To demonstrate the method, we compute the ionization and excitation of H and He in supernova ejecta taking into account time-dependent effects and X-ray irradiation. We find that the interaction with a typical red supergiant wind should result in the enhanced excitation of the outer layers of unshocked ejecta and the emergence of corresponding HV absorption, i.e. a depression in the blue absorption wing of Halpha and a pronounced absorption of He I 10830 A at a radial velocity of about -10,000 km/s. We identify HV absorption in Halpha and He I 10830 A lines of SN 1999em and in Halpha of SN 2004dj as being due to this effect. The derived mass loss rate is close to 10^{-6} Msun/yr for both supernovae, assuming a wind velocity 10 km/s. We argue that, in addition to the HV absorption formed in the unshocked ejecta, spectra of SN 2004dj and SN 1999em show a HV notch feature that is formed in the cool dense shell (CDS) modified by the Rayleigh-Taylor instability. The CDS results from both shock breakout and radiative cooling of gas that has passed through the reverse shock wave. The notch becomes dominant in the HV absorption during the late photospheric phase, ~60 d. The wind density deduced from the velocity of the CDS is consistent with the wind density found from the HV absorption produced by unshocked ejecta.

There is still no consensus about progenitor masses of Type IIP supernovae. We study a normal Type IIP SN 1999em in detail and compare it to a peculiar Type IIP SN 1987A. We computed the hydrodynamic and time-dependent atmosphere models interpreting simultaneously both the photometric and spectroscopic observations. The bolometric light curve of SN 1999em and the spectral evolution of its H-alpha line are consistent with a presupernova radius of 500 Rsun, an ejecta mass of 19.0 Msun, an explosion energy of 1.3x10^51 erg, and a radioactive 56Ni mass of 0.036 Msun. A mutual mixing of hydrogen-rich and helium-rich matter in the inner layers of the ejecta guarantees a good fit of the calculated light curve to that observed. Based on the hydrodynamic models in the vicinity of the optimal model, we derive the approximate relationships between the basic physical and observed parameters. We find that the hydrogen recombination in the atmosphere of a normal Type IIP SN 1999em, as well as most likely other Type IIP supernovae at the photospheric epoch, is essentially a time-dependent phenomenon. It is also shown that in normal Type IIP supernovae the homologous expansion of the ejecta in its atmosphere takes place starting from nearly the third day after the supernova explosion. A comparison of SN 1999em with SN 1987A reveals two very important results for supernova theory. First, the comparability of the helium core masses and the explosion energies implies a unique explosion mechanism for these core collapse supernovae. Second, the optimal model for SN 1999em is characterized by a weaker 56Ni mixing up to 660 km/s compared to a moderate 56Ni mixing up to 3000 km/s in SN 1987A, hydrogen being mixed deeply downward to 650 km/s.

We study a time-dependent hydrogen ionization in the atmosphere of SN 1987A during the first month after the explosion. The model includes kinetics of hydrogen ionization and excitation, molecular hydrogen kinetics, and a time-dependent energy balance. The primary strong effect of the time-dependent ionization is the enhanced hydrogen ionization compared to the steady-state model. The time-dependent ionization provides a sufficient population of excited hydrogen levels to account for the observed H-alpha without invoking the external Ni-56. We find that the Ba II 6142 A line in SN 1987A can be reproduced for the LMC barium abundance. This resolves the long-standing problem of the unacceptably high barium overabundance in SN 1987A. The key missing factor that should be blamed for the "barium problem" is the time-dependent ionization. The modelling of the H-alpha profile on day 4.64 indicates the ratio of the kinetic energy to the ejected mass about 0.83 10^{50} erg/Msun.

We have studied the influence of the presupernova structure and the degree of Ni-56 mixing on the bolometric light curve of SN 1987A in terms of radiation hydrodynamics in the one-group approximation by abandoning LTE and by taking into account nonthermal ionization and the contribution of spectral lines to opacity. Our study shows that moderate Ni-56 mixing at velocities of < 2500 km/s can explain the observed light curve if the density of the outer layers of the presupernova exceeds the value obtained in the evolutionary model of a single nonrotating star severalfold. Abandoning LTE and allowing for nonthermal ionization when solving the equation of state and calculating the mean opacities and the thermal emission coefficient leads to a significant difference between the gas temperature and the radiation temperature in the optically thin layers of the supernova envelope. We demonstrate the fundamental role of the contribution of spectral lines to the opacity in an expanding envelope and of the accurate description of radiative transfer in reproducing the observed shape of the bolometric light curve. We have found that disregarding the contribution of spectral lines to the opacity introduces an error of 20% into the explosion energy, and that a similar error is possible when determining the mass of the ejected matter. The resonant scattering of radiation in numerous lines accelerates the outer layers to velocities of 36 000 km/s; this additional acceleration affects the outer layers with a mass of 10^{-6} Msun. Proper calculations of the supernova luminosity require that not only the delay effects, but also the limb-darkening effects be taken into account.

We have developed a time-dependent model of ionization, excitation and energy balance of SN 1987A atmosphere at the photospheric epoch to study early behavior of hydrogen and helium lines. The ionization freeze-out effects play a key role in producing both the strong H-alpha during nearly the first month and the He I 5876 A scattering line on day 1.76. Using an extended reaction network between hydrogen molecules and their ions demonstrates that ion-molecular processes are likely responsible for the blue peak in H-alpha profile at the Bochum event epoch.

We have shown that in the atmosphere of SN 1987A level populations of hydrogen are mostly controlled by an ionization freeze-out and ion-molecular processes up to 40 days. The ionization freeze-out effects are important for normal SNe II-P as well. The time dependent effects and the mutual neutralization between H^- and H^+ may result in a non-monotonic radial dependence of the Sobolev optical depth for H-alpha and then in a blue emission satellite of H-alpha observed in SN 1987A spectra at Bochum event phase. The relative abundance of molecular hydrogen about 0.0001-0.001 is high enough for the photospheric phase of SN 1987A. We emphasize that the poor knowledge of the far UV radiation field in the envelope of SN 1987A still forbids a truly reliable estimate of the Ba abundance but well-constructed and tested models may improve this situation.

We analyzed the spectra and light curve of the peculiar type II-P supernova 1997D to recover ejecta parameters. The optimal hydrodynamical model of SN 1997D, which meets observational constraints at the photospheric epoch, suggests a low explosion energy of about $10^{50}$ erg, ejecta mass around $6 M_{\odot}$, and presupernova radius near $85 R_{\odot}$. We confirm the previous result by Turatto et al. (\cite{tmy98}) that the ejecta contain a very low amount of radioactive $^{56}$Ni ($\approx 0.002 M_{\odot}$). Modelling the nebular spectrum supports the hydrodynamical model and permits us to estimate the mass of freshly synthesized oxygen (0.02--0.07 $M_{\odot}$). Combined with the basic results of stellar evolution theory the obtained parameters of SN 1997D imply that the progenitor was a star from the 8--12 $M_{\odot}$ mass range at the main sequence. The fact that at least some progenitors from this mass range give rise to core-collapse supernovae with a low kinetic energy ($\approx 10^{50}$ erg) and low amount of radioactive $^{56}$Ni ($\approx 0.002 M_{\odot}$) has no precedent and imposes important constraints on the explosion mechanism. We speculate that the galactic supernovae 1054 and 1181 could be attributed to SN 1997D-like events.