# Radiation hydrodynamics of SN 1987A: I. Global analysis of the light curve for the first 4 months

## 1 Introduction

SN 1987A has provided us with an excellent opportunity to test the theory of massive star evolution, nucleosynthesis, and supernova explosion. The broad band photometric observations, ranging from ultraviolet to far-infrared, and the resultant bolometric light curve enable us to probe the physical processes occurring in the interior of SN 1987A.

The light curve is so sensitive to the hydrodynamics that it is a useful tool to infer the progenitor’s radius, the distribution of elements, mass of the ejecta (in particular, the mass of the hydrogen-rich envelope, $`M_{\mathrm{env}}`$), and the explosion energy, $`E`$ (e.g., Nomoto et al. 1994 for a review). Although earlier theoretical models with flux-limited diffusion approximation are generally in good agreement with observations (see, e.g., Arnett et al. 1989a, and Hillebrandt & Höflich (1989) for a review), there are still some uncertainties. One of the most important complications is that the supernova atmosphere is scattering dominated so that the color temperature is much higher than the effective temperature; the spectrum is a superposition of spectra emerging from layers with different depths and temperatures (Imshennik & Utrobin 1977; Shigeyama et al. 1987; Pizzochero 1990; Höflich & Wheeler 1999 and references therein). This effect is crucial to include in order to constrain model parameters.

Previous works on modeling this effect have used supernova atmospheric codes which rely either on the temperature structure derived from equilibrium diffusion models (Höflich 1991), or on the time-dependent luminosity from one-group radiation-hydro models (Hauschildt & Ensman 1994). To produce light curves, the latter work also used the temperature structure prescribed by the one-group radiation-hydro results. In the present work we attempt to solve the problem of a scattering dominated supernova envelope doing all calculations time-dependently. We calculate the temperature structure self-consistently, and we do not make assumptions on the radiative equilibrium which is strongly violated during the shock breakout.

We analyze the light curve of SN 1987A with a multi-group radiation hydrodynamics code called stella (Blinnikov & Bartunov (1993); Blinnikov et al. (1998)). The calculated broad band ( $`UBV`$, IUE UV) and bolometric light curves are compared with observations for the first $`4`$ months after core collapse. From this we update the constraints on model parameters such as the explosion energy and the extent of mixing.

We discuss the uncertainties in our predictions due to NLTE effects (which are not included in our modeling) by comparing with NLTE atmospheric calculations done by other authors. We believe that those uncertainties can be bracketed for the first months of SN 1987A light by two extreme assumptions on the scattering of photons in spectral lines. It is found that the predictions for the first hours and days are fairly insensitive to the assumptions on spectral lines, because of the overwhelming dominance of electron scattering. So, for those epochs we produce reliable predictions for the soft X-ray/extreme UV flash of the supernova.

## 2 Radiation Hydrodynamics

Earlier modeling of the light curve of SN 1987A has adopted the following types of numerical approaches:

1. Equilibrium-diffusion radiation hydrodynamics, with a flux-limiter to ensure a smooth transition from diffusion to free-streaming regimes (Shigeyama et al. 1987; Arnett 1987, 1988; Grasberg et al. (1987); Woosley 1988; Utrobin (1993)). Here equilibrium means one-temperature approximation where the radiation and gas have the same temperature. Local thermodynamical equilibrium (LTE) is assumed, and black body spectra are used to obtain monochromatic broad-band light curves.
2. Non-equilibrium, one-energy group (gray) radiation hydrodynamics, where equilibrium between the gas and the radiation is no longer assumed (i.e., two-temperature transfer). LTE is assumed and a diluted black body approximation is used for the monochromatic light curves (Ensman & Burrows 1992; Mair et al. 1992).
3. Multi-energy group (non-gray) atmospheric codes combined with gray radiation hydrodynamics. Here the atmospheric structures, i.e., the distributions of temperature, density, and velocity at $`\tau _{\mathrm{sc}}100`$ are obtained from the gray hydro code described above (Höflich 1991; Hauschildt & Ensman 1994) and the emerging spectra are computed taking into account NLTE effects.
4. Full multi-energy group radiation hydrodynamics. This is the approach in our study, and the code we have used is called stella (Blinnikov & Bartunov 1993; Blinnikov et al. 1998); a detailed description of new technical features of the code is given in Sorokina, Blinnikov, & Nomoto (1999). LTE for ionization and atomic level populations is assumed.

stella solves the time-dependent equations for the angular moments of intensity averaged over fixed frequency bands, using up to 300 zones for the Lagrangean coordinate and up to 100 frequency bins (i.e., energy groups). This high number of frequency groups allows one to have a reasonably accurate representation of non-equilibrium continuum radiation. There is no need to ascribe any temperature to the radiation: the photon energy distribution can be quite arbitrary.

The coupling of multi-group radiative transfer with hydrodynamics means that we can obtain the color temperature in a self-consistent calculation, and that no additional estimates of thermalization depth as in the one-energy group model of Ensman & Burrows (1992) are needed. Variable Eddington factors are computed, which fully take into account scattering and redshifts for each frequency group in each mass zone. The gamma-ray transfer is calculated using a one-group approximation for the non-local deposition of the energy of radioactive nuclei. Here we follow Swartz, Sutherland, & Harkness (1995; see also Jeffery 1998), and we only use a purely absorptive opacity. This should be a good approximation. 

In the equation of state, LTE ionizations and recombinations are taken into account. The effect of line opacity is treated as an expansion opacity according to the prescription of Eastman & Pinto (1993; see also Blinnikov et al. 1998). Their approach is different from that of Shigeyama & Nomoto (1990) who used Rosseland mean opacities for scattering and absorption processes, where the line opacities were assumed to be 0.009 and 0.01 cm<sup>2</sup> g<sup>-1</sup> for helium and heavier elements, respectively (Los Alamos opacity library). Our opacities are also different from those in other equilibrium-diffusion or one-group radiation hydro models, since instead of one single energy-averaged opacity we need opacities for all our energy groups.

## 3 Models

We have studied two progenitor models in some detail: the evolutionary model of Nomoto & Hashimoto (1988) and Saio, Nomoto, & Kato (1988b) (see also Shigeyama et al. (1988); Shigeyama & Nomoto (1990); Yamaoka et al. (1991); Saio, Kato, & Nomoto 1988a), and the non-evolutionary model of Utrobin (1993). Here we concentrate on the model studied by Nomoto and co-workers. The results for Utrobin’s model (which gives one of the best fits to the bolometric light of SN 1987A among the equilibrium diffusion models) are presented elsewhere.

In the evolutionary model, a star with the initial mass $`M_{\mathrm{ms}}`$ = 23 $`M_{}`$ and low metallicity $`Z`$ = 0.005 is evolved from the main-sequence and onwards, with the Schwarzschild criterion applied for convection. During the evolution from blue to red, the stellar mass decreases to 16.3 $`M_{}`$ and a helium core of $`M_{\mathrm{core}}`$ = 6.7 $`M_{}`$ is formed. During the red phase, 0.7 $`M_{}`$ of helium is mixed out into the hydrogen-rich envelope yielding $`M_{\mathrm{core}}`$ = 6.0 $`M_{}`$ and $`M_{\mathrm{env}}`$ = 10.3 $`M_{}`$. The dredge up of the helium enhances the surface helium abundance to $`Y_{\mathrm{surf}}`$ = 0.43, which is large enough to move the star in HR-diagram from the red to the location of Sk –69202 in the blue. The resultant luminosity and radius are $`L_0`$ = $`1.3\times 10^5`$ $`L_{}`$ and $`R_0`$ = 48.5 $`R_{}`$, respectively. Because the radius of the progenitor is not well constrained in the evolutionary model, we have also constructed a hydrostatic envelope model with a radius of $`R_0`$ = 40 $`R_{}`$ and of $`R_0`$ = 58 $`R_{}`$ which are fitted to the evolved 6 $`M_{}`$ He core. The 6 $`M_{}`$ He core was evolved through the iron core collapse (Nomoto & Hashimoto 1988). The resultant explosion and explosive nucleosynthesis were calculated as in Hashimoto, Nomoto, & Shigeyama (1989) and Thielemann, Hashimoto, & Nomoto (1990). We have assumed in all our models (i.e., irrespective of the explosion energy) that the mass cut is located at $`M_\mathrm{r}M_\mathrm{c}`$ = 1.6 $`M_{}`$, so that the mass of the ejecta is $`M_{\mathrm{ej}}=14.7`$ $`M_{}`$. Following Shigeyama & Nomoto (1990) we denote our standard model 14E1. Various suffixes have been added to distinguish between different models which have different explosion energies and initial radii, as well as other different physical properties.

For our multi-group computation it is inconvenient to use the initial evolutionary 14E1 model directly. Instead, the star was constructed in hydrostatic equilibrium in the same way as was done in Blinnikov et al. (1998) for SN 1993J: as we need a much finer zoning in the outer layers in our calculations than was used in Nomoto & Hashimoto (1988) and Saio et al. (1988b), a remap of the original model was done onto another grid. We assume that at the outer boundary (i.e., at $`m=M`$, $`M`$ being the total mass of the star) the material pressure vanishes, $`p=0`$, and that there is no radiation coming from the outside. The density structure found in this way is shown in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months. The density of the free blue wind is only $`10^{16}`$ g cm<sup>-3</sup> when scaled as $`r^2`$ from $`10^{17}`$ cm and inward to $`3\times 10^{12}`$ cm (Lundqvist 1999). This falls outside the plot in Figure 1.

As the shock wave propagates through the star, the interfaces where the composition changes suddenly from C+O dominated to helium dominated, and then further out to hydrogen-dominated, are strongly Rayleigh-Taylor unstable. This induces mixing of the material before the shock breakout at the surface (e.g., Bandiera (1984); Ebisuzaki, Shigeyama, & Nomoto (1989); Arnett et al. 1989b ; Benz & Thielemann (1990), Hachisu et al. 1990, 1991, Müller, Fryxell, & Arnett (1991); Basko (1994)). Particularly important for the light curve behavior is the mixing of hydrogen down to the central region. Mixing of <sup>56</sup>Ni into the hydrogen-rich envelope also affects the light curve, and is decisive for spectral (Utrobin, Chugai, & Andronova (1995)) and X-ray observations. The mixed abundance distribution in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months is inferred from the comparison of the X-ray and $`\gamma `$-ray light curves and spectra with observations (Kumagai et al. 1989). The outermost mixed <sup>56</sup>Ni corresponds to the velocity $`4000`$ /kms in the best fitting model of the light curve (see Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months below). The mixing is artificial, since 2D modeling failed to distribute <sup>56</sup>Ni further out than $`2500`$ /kms (Hachisu et al. . 1990, 1991). The density profile was only marginally changed by the mixing.

We have tested explosion energies in the range $`E=(0.71.5)\times 10^{51}`$ ergs and named the runs 14E0.7, 14E1, 14E1.3, …, etc., where ‘14’ denotes the ejecta mass in solar units, see Table 1. All models in the runs labeled without the suffix ‘U’ have the mixed composition shown in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months already prior to the model explosion. Each model was exploded by the deposition of heat energy in a layer of mass $`0.03`$ $`M_{}`$ outside of 1.6 $`M_{}`$. Since stella does not include nuclear burning, preservation of the same mixed composition in the ejecta is assured. The explosion energies given above refers to the asymptotic kinetic energy of ejecta. The heat energy injected for a simulated explosion is higher by $`0.7\times 10^{51}`$ ergs, i.e., by the gravitational binding energy for the presupernova model. A small fraction ($`7\times 10^{48}`$ ergs) goes to the photon energy which is radiated away. For other models, 14E1U, and 14E1.2U, the unmixed composition in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months was used, and in 14E1A, we treated the opacity as pure absorption for the same total extinction (which, however, is often dominated by scattering in reality). For the latter model the entry for “forced $`\chi _{\mathrm{abs}}`$” is “yes”, i.e., $`\chi _{\mathrm{abs}}`$ was artificially set equal to the total extinction. The suffix ‘H’ denotes the model having the same density structure as the standard one, 14E1, but the abundance of hydrogen in the outer layers is enhanced to the solar value. The suffix ‘R’ is added for the models with non-standard initial radius (see Table 1). Some other models from the Table 1 are discussed below.

## 4 Hydrodynamics and Shock Breakout

The shock wave arrives at the surface of the star at time $`t_{\mathrm{prop}}`$, which for different values of the initial radius $`R_0`$, the ejected mass $`M_{\mathrm{ej}}`$, and explosion energy $`E`$, could be approximated by:

$$t_{\mathrm{prop}}1.6\left(\frac{R_0}{50R_{}}\right)\left[\left(\frac{M_{\mathrm{ej}}}{10M_{}}\right)/\left(\frac{E}{1\times 10^{51}\mathrm{erg}}\right)\right]^{\frac{1}{2}}\mathrm{hours}.$$
(1)   
This expression is consistent with Shigeyama et al. (1987). We compare the estimate (1) with computations in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months. The change in the velocity profile in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months shows how the materials are accelerated near the shock breakout until they reach homologous expansion. The maximum velocity at the outer edge is $`32,500`$ km s<sup>-1</sup> for 14E1 (which is somewhat lower than for 14E1.3 in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months). Further acceleration is limited by the inefficiency of the radiative precursor of the shock (see a semi-analytic approach reviewed by Nadyozhin 1994; see also Imshennik & Nadyozhin 1988). We will discuss the influence of various assumptions on the maximum velocity and other details of shock breakout elsewhere. Here we just note that velocities of the order $`(34)\times 10^4`$ km s<sup>-1</sup> are in very good agreement with the value found from the absorption feature of Mg II $`\lambda `$2800 in the early IUE spectra (Pun et al. (1995)). This has important implications for our understanding of the density structure of the circumstellar medium of the supernova (Chevalier 1999; Lundqvist 1999). More specifically, it indicates that the density of the blue supergiant wind must have been very low, corresponding to a mass loss rate of only $`10^8M_{}\mathrm{yr}^1`$.

The density distribution in the outer part is well approximated by a power law $`r^{8.6}`$ as found in Shigeyama & Nomoto (1990), but the very outermost layers are much steeper (see Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months). In between there is a dense shell, which was also found in non-equilibrium radiation hydrodynamic modeling (Blinnikov & Nadyozhin (1991); Blinnikov, Nadyozhin, & Bartunov (1991); Ensman & Burrows (1992)), but missed in the equilibrium diffusion modeling. Note that the density in the central parts computed by stella is much smoother than shown in Figure 6 by Shigeyama & Nomoto (1990).

## 5 Early Light Curve

After the shock breakout, the early light curve up to $`t`$ 25 days is powered by the diffusive release of the internal energy of the radiation field that is established by the shock wave. The bolometric light curve reaches its maximum luminosity of $`L_{\mathrm{bol}}`$ (4–9) $`\times 10^{44}`$ ergs s<sup>-1</sup> (Table 2 and Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months) immediately after shock breakout, and then drops rapidly by many orders of magnitude in $`10`$ days.

The total energy radiated during the first two days amounts to $``$ $`10^{47}`$ ergs (Table 2), but most of the radiation is emitted in a soft X-ray/EUV burst, and was not observed. However, the burst had the important effect that it ionized the surrounding gas (e.g., Lundqvist & Fransson 1996; Sonneborn et al. 1997; Lundqvist 1999). The resultant ionization of the circumstellar material is commented on briefly in §7.1, and will be compared with the observations in greater detail in Lundqvist, Blinnikov, & Bartunov (1999a).

After the burst, the ejecta expand so rapidly that the interior temperature (both of the matter and radiation) decreases almost adiabatically as $`r^1`$. As a result, the bolometric luminosity decreases sharply to $`L_{\mathrm{bol}}`$ (2 – 3) $`\times 10^{41}`$ ergs s<sup>-1</sup> to form a minimum of the bolometric light curve. Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months demonstrates that the model with $`E=10^{51}`$ ergs (i.e, 14E1) gives the best agreement with the observed bolometric flux. Note that the agreement for 14E1 is much better than in Figure 7 of Shigeyama & Nomoto (1990)). (Note also the higher resolution in our figure than in the figure of Shigeyama & Nomoto (1990)). The light curves computed by the flux-limited diffusion in Shigeyama & Nomoto (1990) produce a short plateau, $``$ 20 days (see Figs. 16 – 19 in Shigeyama & Nomoto (1990)). This is neither seen in the observations, nor in our models. The light curves computed here by the full radiative transport are in better agreement with observations for the same models (cf. Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months) suggesting that the more accurate method of modeling gives results which are closer to reality. 

The luminosity at this phase is lower than for typical Type II-P supernovae by a factor of 10 – 20. This is due to the small initial radius which leads to a low luminosity because a much larger fraction of the radiation field energy is lost by $`P`$d$`V`$ work than in ordinary Type II-P supernovae. This is well-known from early modeling of low-luminosity Type II-P supernovae (Imshennik & Nadyozhin (1964); Chevalier 1976). For SN 1987A, the low luminosity was successfully demonstrated in the models of Shigeyama et al. (1987), Arnett (1987), Grasberg et al. (1987), Woosley (1988), Woosley, Pinto, & Eastman (1988) and Utrobin (1993).

In Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months, we show the changes in the color temperature, $`T_\mathrm{c}`$, of the best blackbody fit to the flux, along with the effective temperature, $`T_{\mathrm{eff}}`$, defined by the luminosity and the radius of last scattering $`R`$ through $`L=4\pi \sigma T_{\mathrm{eff}}^4R^2`$ (see Blinnikov et al. , 1998, for details of finding $`R`$ and from that $`T_{\mathrm{eff}}`$). The maximum value of $`T_\mathrm{c}`$ is $`1.2\times 10^6`$ K for model 14E1 (see Table 2), which is higher than the $`7.6\times 10^5`$ K in Höflich’s (1991) NLTE time-dependent calculation for 14E1.25, but similar to the temperatures found by Ensman & Burrows (1992). Our results are in very good agreement with the estimates of Imshennik & Nadyozhin (1988, 1989). We emphasize that our multi-group radiative transfer with hydrodynamics obtains this temperature in a self-consistent way, and no additional estimates of the thermalization depth (like in the one-group model of Ensman & Burrows 1992) are needed. The large difference between color and effective temperatures is due to scattering (e.g., Sobolev 1980; Kolesov & Sobolev 1982; Höflich 1991; Wagoner & Montes 1993); the average energy of the photons is higher than that corresponding to the value of $`T_{\mathrm{eff}}`$. The effect is a deficit of photons in the visual at maximum light compared to a model with forced absorption (Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months). This effect continues throughout the first day after shock breakout. In Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months we demonstrate the dominance of scattering in extinction for high temperatures.

Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months shows $`T_\mathrm{c}`$ for the same three runs as in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months, and one can see that $`T_\mathrm{c}`$ is rather insensitive to $`E`$ (except for the maximum $`T_\mathrm{c}`$ at shock breakout, as displayed in Table 2).

During the first day, the visual luminosity increases because the intensity peak is rapidly shifted into the optical due to the decreasing photospheric temperature. In order for the optical flare-up of the supernova to be seen at 6 mag at $`t=3`$ h, the condition $`t_{\mathrm{prop}}<3`$ h (Eqn. ) should be satisfied, which requires a relatively large $`E/M_{\mathrm{env}}`$ and small $`R_0`$ (Shigeyama et al. 1987; Woosley et al. 1987). Also, the ejecta and the radiation field should have expanded rapidly so that the radiation temperature becomes lower and the radius of the photosphere becomes larger sufficiently fast. Therefore, the expansion velocity, and thus $`E/M_{\mathrm{env}}`$, should be larger than certain values for a given initial radius.

Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months shows the $`V`$ light curve for the model 14E1.21 with realistic scattering dominated opacity for the first hours of SN 1987A. Here, and for other models below, we have used a distance modulus of 18.5. (See Walker 1999, for a discussion on the uncertainty of this value.) We note that the $`V`$ curve exhibits an early local minimum which does not exist in equilibrium diffusion models (e.g., Woosley (1988); Arnett (1988); a small minimum was found also by Utrobin (1993)). Höflich (1991) ascribed this to a non-LTE effect, and used the location in time of the minimum to constrain $`E`$. However, we too recover the local minimum despite our LTE approach. The reason for the minimum is that when the bolometric flux continues to fall (Fig. 9) after shock breakout, the bolometric correction overcomes the effect of the falling bolometric flux. The visual luminosity then increases because the intensity peak is rapidly shifted into the optical due to the decreasing photospheric temperature. Regardless of the exact cause of the difference between our and Höflich’s results on the one hand, and equilibrium diffusion models on the other, the first maximum in $`V`$ in our models is just on the level of Jones’ limit (Wampler et al. (1987)) even for our 14E1.3 run, where the energy of explosion is $`9\times 10^{49}`$ ergs higher (see Table 1) than the model 14E1.25 used by Höflich (1991; see also Höflich & Wheeler 1999). The earlier appearance of the peak of the $`V`$-flux in our run is due to this higher energy.

Compared with observations, and especially the first $`V`$ observation by McNaught & Zoltowski (1987; later revised by West & McNaught 1992), the calculated $`V`$ curve in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months rises too slowly. In equilibrium diffusion modeling (e.g., Arnett (1988); Woosley (1988); Shigeyama & Nomoto (1990)) the $`V`$ flux rises much faster. The reason for this difference is that equilibrium diffusion models do not care whether the total extinction is absorptive or due to scattering. The spectra are assumed in those models to be blackbody with the temperature equal to $`T_{\mathrm{eff}}`$. This is not a good approximation, since in reality the thermalization of photons takes place well below the surface of the last scattering resulting in $`T_\mathrm{c}T_{\mathrm{eff}}`$. In our run 14E1A we have assumed purely absorptive extinction, i.e., we force even the electron scattering to act as true absorption so that the thermalization of the photons occurs at an optical depth of order unity. We thereby reproduce the results obtained in equilibrium diffusion models. As a matter of fact, the $`V`$ curve in 14E1A rises even faster than was observed. (See Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months which shows the $`UBV`$ light curves for 14E1 and 14E1A for the first 5 days.) There is of course no physical reason why electron scattering should be treated as absorptive. From this it is clear that that an accurate treatment of scattering is crucial to constrain $`E`$. We will discuss this now.

Höflich (1991) was the first to show that non-equilibrium effects influence the $`V`$ light curve of SN 1987A drastically; the flux in $`V`$ is $`2`$ magnitudes lower for the first hours than in the equilibrium gray atmosphere case. Höflich (1991) tried a higher explosion energy than in 14E1 of Shigeyama & Nomoto (1990) to compensate for this reduction of flux by adopting the model 14E1.25 computed in the flux-limited equilibrium diffusion approximation by Shigeyama & Nomoto (1990). He then found a marginal agreement with the first data points of McNaught & Zoltowski (1987). The same results have since been presented in Höflich & Wheeler (1999). Compared with our $`V`$ light curve of 14E1.21 and 14E1.3 (see Figs. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months and Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months) the rise of the $`V`$ curve in Höflich’s (1991) 14E1.25 model is significantly faster. To clarify this discrepancy, one should note that Höflich used the temperature structure of an equilibrium diffusion model, which can be appreciably different from that in our full transport models. We cannot say with certainty that this fully explains why our results and those of Höflich differ so much. It could be due to that the expansion opacity should be treated differently in the energy equation than is described in Höflich (1990; see Blinnikov 1996, 1997). 

There are, of course, uncertainties also in our models. In particular, the role of NLTE effects needs to be further examined, but we cannot envisage that they are able to explain the big difference between us and Höflich (1991) for the first hours after shock breakout. At this epoch the extinction is totally dominated by electron scattering and spectral lines are not so important, see Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months. We see virtually no difference in our results for the first day when we treat the lines as fully absorptive, or as totally scattering dominated. The NLTE effects set in later, and are very important after a couple of weeks (cf. Baron et al. (1996)). We are therefore confident that our results for this epoch are more accurate than Höflich’s (1991).

We note that even Höflich’s $`V`$-flux is lower than the observed after the new reductions by West & McNaught (1992). A possible explanation to the discrepancy between models and observations is that the stars, used to calibrate the early plates of the supernova by West & McNaught (1992), are too cool for an object with a color temperature of $`T_\mathrm{c}10^5`$K. (We find at $`t=0.128`$ day that $`T_\mathrm{c}=1.14\times 10^5,1.05\times 10^5,9.8\times 10^4`$ K for the runs 14E0.7, 14E1, and 14E1.3, respectively.) The temperature at $`t=0.128`$ days decreases with increasing explosion energy because of the earlier emergence of the shock and the faster adiabatic cooling. It should be emphasized that our models are much hotter than equilibrium diffusion models (e.g., Arnett (1988); Woosley (1988); Utrobin (1993)), and a comparison between observations and our results is therefore more sensitive to calibration errors than are equilibrium diffusion models. As we will see in §6.1 our models fit the early IUE observations well, and since these observations are less likely to have the same error, we cannot exclude calibration errors to be the cause of the mismatch in $`V`$. This would mean that the early true $`V`$ flux was lower than hitherto believed. We point out that we have changed various parameters in our models to try to make our $`V`$ flux increase faster and thereby fit the observations better. These experiments included enhancing the iron abundance, and varying the explosion energy and presupernova model within the limits allowed by the global light curve. However, none of these attempts reduced the discrepancy. (See, e.g., the results for two different initial radii in Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months.)

## 6 Before and After the Peak of the Light Curve

After the minimum around day $``$ 10, the observed bolometric light curve showed an almost exponential increase up to day $``$ 60, and subsequently formed a plateau-like broad peak around day $``$ 100. After a relatively rapid drop, the luminosity then declined slowly between $`t`$ = 120 – 400 days at the rate of <sup>56</sup>Co-decay. The energy source responsible for the broad peak of the light curve, and the tail, is therefore without doubt the radioactive decay of <sup>56</sup>Ni$``$ <sup>56</sup>Co$``$ <sup>56</sup>Fe . The total mass of initial <sup>56</sup>Ni in our models is $`M_{\mathrm{Ni}}0.078M_{}`$, which is the same mass as in the models of Shigeyama & Nomoto (1990). The theoretical bolometric light curves for models 14E1, 14E1M, and 14E1U with a different extent of mixing are shown in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months. It is clearly seen that the shape of the modeled light curve is strongly dependent on the distribution of hydrogen and <sup>56</sup>Ni in the ejecta, i.e., the amount of mixing that has occurred.

For the model with standard mixing (14E1 with <sup>56</sup>Ni mixed out to $`v`$ 4000 km s<sup>-1</sup>, Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months), there is significant heating of the outer layers due to radioactivity. The heating becomes noticeable in the light curve already at $`10`$ days, and then forms a smooth increase in the optical light curve up to the peak, as observed.

For the unmixed model 14E1U (Figs. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months, Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months), the shape of its bolometric light curve starts to differ from that of 14E1 already at $``$ 10 days. For the unmixed case, the increase in the luminosity due to radioactive heating is delayed until $`t=35`$ days, which causes a dip to appear in the light curve around day 30, and makes the light curve in the subsequent phase rise faster than in the mixed model. These properties are clearly incompatible with the observations.

For the case with intermediate mixing (14E1M), <sup>56</sup>Ni is mixed out to $`v`$ 2500 km s<sup>-1</sup>, while hydrogen is mixed as in the mixed model 14E1 (Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months). Also in this model the appearance of radioactive heating occurs too late to be compatible with observations (see Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months). This suggests that mixing of <sup>56</sup>Ni out to $`v`$ 4000 km s<sup>-1</sup> is needed. This conclusion is independent of the radiation transfer scheme used (Nomoto, Shigeyama, & Hashimoto 1987; Woosley 1988; Shigeyama et al. 1988). The mixing out to large velocities is also supported by the redshifted feature at $`3900`$ km s<sup>-1</sup> observed by Haas et al. (1990; see also Utrobin, Chugai, & Andronova (1995)). We note, however, that Utrobin (1993) obtains a good bolometric light curve without <sup>56</sup>Ni mixing, and that Kozma & Fransson (1998a,b) do not need to mix nickel out to more than $``$ 2000 km s<sup>-1</sup> to model iron line profiles at late epochs. We postpone a more detailed discussion on mixing to a future paper.

It is also interesting to compare the results of our non-equilibrium radiative transfer modeling by stella with the models in Figures 16 – 19 of Shigeyama & Nomoto (1990). One striking difference, already noted above, is that the shape of the minimum near day 10 is reproduced much better with our non-equilibrium modeling. Another difference is that the shape of the light curve around maximum (at $`3`$ months) is smoother in models with realistic scattering opacity (see Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months below), than in models with forced absorption like 14E1A (see Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months), and hence also in equilibrium diffusion models since the latter two types of model are closely related. These two models also show the same postmaximum sharp decline (cf. Utrobin (1993)), which can be understood in terms of enhanced emission according to Kirchoff’s law.

We note that the rising part of the light curve is modeled better with enhanced absorption in spectral lines than with scattering lines (see Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months for the bolometric luminosity, and Figs. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months, Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months for $`UBV`$ colors). This can perhaps be explained by NLTE effects. In particular, the effect of fluorescence could be important (see Baron et al. (1996); Li, McCray, & Sunyaev (1993); Li & McCray (1996), and references therein). This cannot be studied directly by stella, but, as found by Baron et al. (1996), forced absorption in spectral lines can reproduce some properties of the fluorescence as fluorescence is a form of thermalization; spectra of LTE models with absorptive lines are very similar to full NLTE spectra, while LTE models with scattering lines are far from reality (see also Eastman 1997; Blinnikov et al. (1998)). This hints why the rising part is modeled well when we apply forced absorption in lines.

Another cause for the deviation of the rising part of the light curve from what was observed could be that the distribution of hydrogen is different from that in our model. This effect was investigated by Utrobin (1993). The hydrogen distribution affects the light curve because it determines how the hydrogen recombination front propagates into the ejecta (e.g., Nadyozhin 1994; Shigeyama & Nomoto 1990). Because electron scattering is the main source of opacity, the opacity decreases sharply in the outward direction at the recombination front. This causes the photosphere to become associated with the recombination front. When the ejecta pass through this front their temperature quickly decreases to $`5500`$ K. This is seen as a change in the temperature profile in Figures Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months and Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months. One can locate the photosphere in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months at the point where the radiation temperature starts to deviate from the material temperature appreciably. We note that the recombination front is much broader than that in Shigeyama & Nomoto (1990) because of the large contribution of line opacity to the total opacity in our calculations. The effects of non-equilibrium transport are also important for the width of the recombination front.

As seen in Figures Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months and Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months the photosphere propagates inward in mass, while the material expands outward. For a certain period, the recombination front is almost stationary in radius (Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months). If $`T_{\mathrm{eff}}`$ were constant, then this would result in an almost constant bolometric luminosity, i.e., in a perfect plateau of the light curve, where the duration of the plateau stage depends on how deep into the star hydrogen has been mixed. For a deeper mixing of hydrogen, the plateau lasts longer. When the photosphere enters into layers which lack sufficient amounts of hydrogen, the hydrogen recombination front disappears and the plateau phase is terminated. This is the case for typical Type II-P supernovae, where radioactive heating does not show up until the very end of the plateau stage (Eastman et al. 1994). SN 1987A is quite different in this respect, since radioactivity is important also at early epochs; the diffusion flux from the radioactive energy release starts to dominate over the diffusion flux from the recombination of hydrogen already at $`6`$ weeks after the explosion. This is clearly seen in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months (cf. Fig. 12 in Shigeyama & Nomoto (1990), and Fig. 11 in Utrobin (1993)).

The relation between the duration of the plateau phase and the depth of the hydrogen layer can be given more quantitatively. Suppose that hydrogen is mixed down to a shell where the expansion velocity of the hydrogen-rich layer is $`v_\mathrm{H}`$. Then this velocity is related to the observed quantities for the plateau phase as

$$v_\mathrm{H}1300\left(\frac{L_{\mathrm{pl}}}{8.5\times 10^{41}\text{ ergs s}\text{-1}}\right)^{1/2}\left(\frac{t_{\mathrm{pl}}}{100\text{ d}}\right)^1\text{ km s}\text{-1}$$
(2)   
where $`t_{\mathrm{pl}}`$ is the time at the end of the plateau, $`L_{\mathrm{pl}}=4\pi \sigma T_{\mathrm{eff}}^4R^2`$ the luminosity at $`tt_{\mathrm{pl}}`$, and $`T_{\mathrm{eff}}`$ 5500 K because of the association of the photosphere with the hydrogen recombination front. The time $`t`$ is for the freely expanding ejecta, $`t=R/v`$. For example, if we take $`L_{\mathrm{pl}}=2\times 10^{41}`$ ergs s<sup>-1</sup>, and $`t_{\mathrm{pl}}=30`$ days, as in the unmixed model 14E1U (Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months), we find from Eqn.(2) $`v_\mathrm{H}2000`$ km s<sup>-1</sup> which is in good agreement with Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months.

In the standard model, 14E1, hydrogen is mixed down to $`M_\mathrm{r}2M_{}`$, i.e., only $`0.4M_{}`$ outside $`M_\mathrm{c}`$. The expansion velocity at that radius is only $``$ 1000 km s<sup>-1</sup>. If we adopt $`t_{\mathrm{pl}}100`$ days, Eqn. 2 gives $`v_\mathrm{H}`$ 1300 km s<sup>-1</sup> as is required from observations. Equation (2) thus appears to give a reasonable estimate of $`v_\mathrm{H}`$ when we put $`t_{\mathrm{pl}}`$ equal to the time when the second maximum ends. However, we caution that this result should not be overinterpreted. Equation (2) gives the photospheric velocity as long as $`T_{\mathrm{eff}}`$ is 5500 K, which is roughly the case for SN 1987A also for the second maximum. But the formula also assumes that the temperature is governed by the presence of hydrogen. This is not a unique statement, as the complicated thermal balance may settle around this temperature also without the presence of hydrogen. Even if hydrogen is mixed far into the core, the end of the second maximum is likely to give only limited information about the minimum velocity of hydrogen.

### 6.1 Broadband fluxes

Ideally, both spectra and $`UBV`$ colors should be obtained by full NLTE modeling (see arguments put forward by Eastman 1997). For SN 1987A spectral NLTE modeling has been made by several groups (see, e.g., Schmutz et al. (1990); Höflich (1990); Takeda (1991); Mazzali, Lucy, & Butler (1992); Duschinger et al. (1995); Mazzali & Chugai (1995); Höflich & Wheeler (1999) and references therein), while NLTE modeling of colors are given in the literature only for the first days after the explosion (Eastman & Kirshner (1989); Hauschildt & Ensman (1994)). Our approach cannot add to this NLTE modeling, but as discussed by, e.g., Höflich (1995, see also references therein), LTE modeling is useful even for such rapidly expanding objects as Type Ia SNe, and should therefore also give some insight to the physical conditions of SN 1987A. Inspired by this, we compare our LTE results with observations in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months. The $`B`$ and $`V`$ light curves are in surprisingly good agreement with observations. Likewise, for the $`U`$ band the agreement is satisfactory for the first days, though at later epochs the modeled flux is too high. Nevertheless, the shape of the light curve qualitatively reproduces the observations.

There is also a qualitative agreement between our calculated fluxes and the observed fluxes in IUE bands (Pun et al. (1995)) (see Figs. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 monthsRadiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months). In particular, the agreement is good for the first 10 days for the SWP and the two LWP bands, while the modeled flux overshoots early for the LWP 3000 – 3300 Å band, like it does for the modeled $`U`$ flux. Later, the modeled flux in the bands with the shortest wavelengths undershoots, while for the range 2500 – 3000 Å, the agreement is still fairly good. The modeled flux in the band with the longest wavelengths continues to be high for the first 100 days (Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months).

The disagreement between modeled and observed UV fluxes after $`100`$ days is not surprising, because the LTE modeling at that epoch becomes quite unrealistic when there is no longer a true photosphere. For earlier epochs, there can be several causes for deviations. The early drop of the modeled fluxes in the shortest IUE bands could signal that the photosphere is somewhat hotter than in the 14E1 model; the flux here falls into the Wien part of the spectrum and it is exponentially sensitive to the temperature. It is harder to explain why the observed UV flux for $`\lambda >3000`$ Å and the flux in the $`U`$ band fall below what we predict. Within the LTE approach this could mean that the expansion opacity is not complete in the Eastman-Pinto (1993) approximation. The line list of Eastman & Pinto includes $`10^5`$ lines, but they and Baron et al. (1996) have pointed out that it is also necessary to include millions of weaker lines (see also Höflich 1995), mostly of iron. Should the Eastman-Pinto list be sufficient, then one might think of enhancing the abundance of iron group elements in the outer layers of the supernova to increase the opacity. Our experiments with such an enhancement show that we can bring the $`U`$ flux in much better agreement with observations. There could also be other causes for our too strong UV flux, but it seems reasonable to assume that line opacities are somehow involved. This is highlighted by the good agreement between modeled and IUE fluxes in the 2500 – 3000 Å range where there are relatively few spectral lines. This could perhaps indicate that the envelope is contaminated with heavy elements which could give a higher opacity in the UV.

In this context we note that Pun et al. (1995) do not cite Wagoner, Perez, & Vasu (1991) correctly when they state that “the expansion opacities in the wavelength region 1000–4000 Å increase by a factor of more than 100 as the temperature of the atmosphere drops from 12,000 to 5000 K”. We point out that it is not the expansion opacity in Wagoner et al. (1991; see also Eastman & Pinto (1993)) which increases by this number. The correct statement is that the ratio of expansion opacity to electron scattering increases by a factor of $`100`$. It does so because the electron scattering opacity drops drastically due to recombination as the temperature is lowered.

We have not tried to include millions of weak lines in our models to check whether this can bring the UV observations and our LTE models in better agreement. If such an experiment would fail, and the metal content of the envelope is not unusually high (cf. above), then one has to consider NLTE effects already for a few days after the explosion. The main NLTE effect here could be the excitation of hydrogen from its second principal level, $`n=2`$, perhaps creating an optically thick Balmer continuum. From NLTE atmospheric calculations (see, e.g., Schmutz et al. (1990); Takeda (1991); Duschinger et al. (1995)) we know that early overpopulation of $`n=2`$ is present in SN 1987A, though it is not sufficient to explain the observed absorption in the Balmer range, unless there is a direct excitation due to radioactivity. While it is certainly important to include nonthermal excitation in late spectra hydrogen (e.g., Xu et al. (1992); Kozma & Fransson 1998b), the same effect operating at early times could be one more hint of extremely efficient outward mixing of radioactive material into the hydrogen-rich envelope. The possibility that hydrogen could be excited as a result of circumstellar interaction seems much more unlikely because of the very low density inferred for the circumstellar gas (§4; Chevalier 1999; Lundqvist 1999).

### 6.2 Parameters of the Photosphere

We present in Figures Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 monthsRadiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months the comparison of our numerical results with the “photospheric” parameters found by observers. We put this in quotes since what is given by observers are not really the parameters of the photosphere, but the best blackbody fit temperature $`T_{\mathrm{obs}}`$ (which is higher than $`T_{\mathrm{eff}}`$), and the radius, $`R_{\mathrm{obs}}`$, found from $`L=4\pi \sigma T_{\mathrm{obs}}^4R_{\mathrm{obs}}^2`$. We emphasize that $`R_{\mathrm{obs}}`$ is substantially smaller than the radius of the true photosphere, $`R_{\mathrm{ph}}`$, especially at the earliest stages.

## 7 Dependence on Model Parameters

### 7.1 Explosion Energy

In the above discussion, the light curve has mainly been used to probe the internal abundance distribution. Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months shows how the bolometric light curve depends on $`E`$. For a given ejecta mass and using the standard mixing, we can find constraints on $`E`$ from the light curve (Shigeyama et al. 1987, 1988; Nomoto et al. 1987, 1994; Woosley 1988; Woosley et al. 1988; Arnett & Fu 1989; Imshennik & Popov (1992)). First, the luminosity near the minimum around day 10 (which corresponds to the early short plateau in the $`V`$ curve) is almost proportional to $`E`$ (Litvinova & Nadezhin 1990, Popov 1993), thus providing an important constraint on $`E`$. Second, the time of the peak, $`t_{\mathrm{peak}}3`$ months, depends on $`E`$. For larger $`E`$, i.e., faster expansion, the rise starts earlier because of earlier appearance of heating due to radioactivity; the decline after the peak is earlier because of a larger velocity of the ejecta which causes the diffusion time-scale to be shorter. The analytical treatment of this epoch is given in detail by Imshennik & Popov (1992). We caution that it is not just $`E`$ that determines $`t_{\mathrm{peak}}`$ (for a given hydrogen distribution), but the combination $`E/M_{\mathrm{env}}`$. $`t_{\mathrm{peak}}`$ is therefore mainly determined by $`E/M_{\mathrm{env}}`$ and the hydrogen distribution.

Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months shows that both the first and the second parts of the light curve are reproduced well by 14E1 (see also Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months). 14E1.3 is too bright near the minimum, while 14E0.7 is too dim near the minimum and evolves too slowly. Compared with the flux-limited diffusion model (Shigeyama & Nomoto 1990) 14E1 evolves similarly, but in much better agreement with observations both near the minimum and near the peak when a more accurate radiative transfer scheme is used.

For the progenitor we have used in Figure Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months (i.e., the model with mixing and $`M_{\mathrm{env}}=10.3`$ $`M_{}`$), the best explosion energy in order to fit the observations is close to $`1.1\times 10^{51}`$ ergs. Considering the uncertainties of the progenitor model in terms of mixing and envelope mass, we obtain best fits to the bolometric light curve for the explosion energies in the range $`(0.851.35)\times 10^{51}`$ ergs. As noted in Figure 9, the “true” bolometric is likely to be intermediate to the SAAO and CTIO/ESO results used in our fits (see Suntzeff & Bouchet 1990). This allows for a $`\pm 10\%`$ span in luminosity. However, a larger error in luminosity ($`15\%`$, Lundqvist et al. 1999b) is due to the still prevailing uncertainty in distance modulus to the supernova (e.g., Walker 1999), which gives a combined error of approximately $`\pm 30\%`$ in absolute luminosity. The explosion energy should therefore be within the range $`(0.81.4)\times 10^{51}`$ ergs.

We note that $`E`$ is very similar to in the analytical diffusion models of Arnett & Fu (1989) and Imshennik & Popov (1992). (The initial theory of diffusion was developed by Arnett  for Type Ia supernovae.) All these models use an eigenvalue formulation of the problem, which is not strictly correct (Blinnikov & Popov 1993). However, the more correct, though more complicated, moving-boundary formulation produces results which agree rather well with the results of Imshennik & Popov (1992; see Popov 1995). We therefore support the use of the results of Imshennik & Popov (1992) to make reliable estimates of the supernova parameters from the near maximum light.

This rather limited range of energies we find is important for calculations of the ionization of the circumstellar gas (e.g., Lundqvist & Fransson 1996; Lundqvist 1999). Until now these calculations have been based on the models 500full1 and 500full2 of Ensman & Burrows (1992). Qualitatively, the 500full1 model is rather similar to the models in our preferred energy range. We will discuss this in detail in Lundqvist et al. (1999a), but from the analysis in Lundqvist & Fransson (1996) we note immediately that the ionization of the outer rings is particularly sensitive to the spectrum of the burst.

### 7.2 Radius of the Progenitor

The above analysis has been based on a progenitor model with $`R_0=`$ 48.5 $`R_{}`$. The uncertainty in the luminosity of the progenitor may imply that $`R_0`$ could be uncertain by $``$ 20 % (Woosley (1988); Saio et al. 1988b ). Thus we have calculated the light curve models 14E1.26R, 14E1.34R, and 14E1.45R for $`R_0=`$ 40$`R_{}`$, and one model, 14E1.4R6, for $`R_0=`$ 58$`R_{}`$ (Table 1 and Figs. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months, Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months). The early light curve for the first 2 days is not so different from the case of the standard $`R_0=`$ 48.5 $`R_{}`$, but $`t_{\mathrm{prop}}`$ scales as $`R_0/E^{1/2}`$.

However, the light curve near the minimum (plateau in $`V`$) is dimmer for $`R_0=`$ 40$`R_{}`$ because the luminosity at that epoch is approximately proportional to $`R_0`$. As a result, for $`R_0=`$ 40$`R_{}`$, we need an explosion energy larger than $`1.2\times 10^{51}`$ ergs to put the light curve in agreement with the observations near the minimum (Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months). But then the peak is too early. A larger radius, $`R_0=`$ 58$`R_{}`$, shows an opposite trend. In this case, one needs a lower energy for the brightness at the minimum, but then the peak is too late.

### 7.3 Hydrogen mass

We have also made some runs for models where the hydrogen abundance in the outer layers was artificially raised to 0.7 (as is also assumed in, e.g., the non-evolutionary models by Utrobin (1993)). For example, the model 14E1.25H (Table 1) has the same composition of metals as 14E1, but H is enhanced at the expense of He. So, the total mass of H is here $`M_\mathrm{H}7M_{}`$, while it was $`M_\mathrm{H}5.5M_{}`$ in the standard runs (14E1, 14E1.3 etc.). The hydrogen rich model 14E1.25H evolves similarly to the standard runs 14E1.3 and 14E1.21 (Fig. Radiation hydrodynamics of SN 1987A:  I. Global analysis of the light curve for the first 4 months). This explains the success of such models, as demonstrated by Utrobin (1993), but it is hard to justify a solar H abundance for SN 1987A from the evolutionary point of view, as well as in context of the abundances in the inner circumstellar ring (e.g., Lundqvist & Fransson 1996).

## 8 Conclusions

In this paper we have described an extensive set of full radiation hydrodynamics calculations aimed to improving the modeling of the first few months of the light curve of SN 1987A. We have shown that the improved models can reproduce the light curve and suggest that proper handling of the radiation transfer is indeed decisive for the success of model fits.

The full multi-group radiation hydrodynamic modeling is more reliable mainly because the effects of scattering are treated self-consistently. Our findings are that:

1. the color temperatures and broad band photometry (and full continuum spectra) are predicted correctly for the first days of the supernova.
2. the shape of the modeled light curve, especially near the bolometric minimum at $`10`$ days, and during the broad peak at $`100`$ days is much better in the multi-group approach than in the equilibrium diffusion one.
3. the density profiles of the supernova at various epochs are smoother.

Our code assumes LTE, but we have bracketed NLTE effects by extreme assumptions of the treatment of spectral lines (as either being scattering dominated, or fully absorptive). We find that the emission at shock break-out is not sensitive to those assumptions, which gives us confidence in our results. This is supported by a very good agreement with IUE observations for the first days. For the first 100 days, best agreement is obtained when NLTE effects are mimicked by treating the line opacity as absorptive, following the prescription of Baron et al. (1996).

Looking at our results in greater detail, we find that the color temperature around shock breakout exceeds $`10^6`$ K, which is higher than those obtained using a more approximate approach (e.g., Shigeyama & Nomoto 1990), but not so different from the model 500full1 of Ensman & Burrows (1992). The large difference between color and effective temperatures persists for the first hours. This implies that the rise in the $`V`$ luminosity is slower than in equilibrium diffusion models (which was first noticed by Höflich 1991). Because of the overwhelming dominance of electron scattering during the first hours, the rise in the $`V`$ band is too slow in all of our models. This raises the question on the calibration of the first photometry data of SN 1987A (McNaught & Zoltowski 1987; West & McNaught 1992).

It is indicated by our light curve modeling that mixing of <sup>56</sup>Ni up to $`v30004000`$ km s<sup>-1</sup> could be needed, but our analysis is unlikely to supersede those of spectral analysis (e.g., Utrobin, Chugai, & Andronova (1995), Kozma & Fransson 1998a,b), modeling of early X-ray emission (e.g., Kumagai et al. 1989), or direct observations of infrared IR lines (Erickson et al. 1988).

We have improved on the constraining of $`E`$. The earlier flux-limited diffusion calculations (Shigeyama & Nomoto 1990) provided a constraint on $`E`$ from both the pre-peak light curve and the plateau-like maximum light, concluding $`E=(1.1\pm 0.4)\times 10^{51}`$ ergs. We find from our more detailed analysis that the best agreement with the observations is obtained for $`E=(1.1\pm 0.3)\times 10^{51}`$ ergs. To arrive at this result we have assumed that the most likely range of $`M_{\mathrm{env}}`$ is $`M_{\mathrm{env}}`$ = 7 - 10 $`M_{}`$ (Saio et al. 1988b), in order for models of the presupernova evolution to be consistent with the enhancement of N/C and N/O in the circumstellar matter (Lundqvist & Fransson 1996). Knowing the energy with this accuracy, as well as having a detailed spectroscopic evolution from our models, we can constrain the ionization of the circumstellar gas much better than before. This will be discussed in Lundqvist et al. (1999a).

We are grateful to Ron Eastman, Stan Woosley, Wolfgang Hillebrandt, Vladimir Imshennik, Dmitriy Nadyozhin, Nikolai Chugai, Cecilia Kozma, Alexandra Kozyreva, Jason Pun and Victor Utrobin for discussions, and also Ron Eastman, Vlad Popolitov and Elena Sorokina for using their computer codes as part of ours. Preliminary results of this work were reported at the “CTIO/ESO/LCO Workshop on SN 1987A: Ten years after” in the contributions by Lundqvist & Sonneborn (1999) and Nomoto, Blinnikov, & Iwamoto (1999). We would like to thank the participants of this workshop for stimulating discussions. This work was supported in Russia by a grant from The International Science & Technology Center 97-370, and by the Russian Basic Research Foundation grants RBRF 96-02-19756 and RBRF 96-02-17604. We are also grateful to grants from The Royal Swedish Academy of Sciences and The Wenner-Gren Center Foundation for Scientific Research. P.L. receives further support from The Swedish Natural Science Research Council and The Swedish Board of Space Research. The project has in part also been supported by the grant-in-Aid for Scientific Research (05242102, 06233101) and COE research (07CE2002) of the Ministry of the Education, Science, Culture, and Sports in Japan.