format: text

text:
# 𝐵=3 Tetrahedrally Symmetric Solitons in the Chiral Quark Soliton Model

## Abstract

In this paper, $`B=3`$ soliton solutions with tetrahedral symmetry are obtained numerically in the chiral quark soliton model using the rational map ansatz. The solution exhibits a triply degenerate bound spectrum of the quark orbits in the background of tetrahedrally symmetric pion field configuration. The corresponding baryon density is tetrahedral in shape. Our numerical technique is independent on the baryon number and its application to $`B4`$ is straightforward.

The Chiral Quark Soliton Model (CQSM) was developed in 1980’s as an effective theory of QCD interpolating between the Constituent Quark Model and Skyrme Model diakonov ; cqsm . In the large $`N_c`$ limit, these models are identical manohar .

The CQSM is derived from the instanton liquid model of the QCD vacuum and incorporates the non-perturbative feature of the low-energy QCD, spontaneous chiral symmetry breaking. The vacuum functional is defined by;

$`𝒵={\displaystyle 𝒟\pi 𝒟\psi 𝒟\psi ^{}\mathrm{exp}\left[id^4x\overline{\psi }\left(i\overline{)}mU^{\gamma _5}\right)\psi \right]}`$ (1)
where the SU(2) matrix

$`U^{\gamma _5}={\displaystyle \frac{1+\gamma _5}{2}}U+{\displaystyle \frac{1\gamma _5}{2}}U^{}\mathrm{with}U=\mathrm{exp}\left(i\stackrel{}{\tau }\stackrel{}{\pi }/f_\pi \right)`$
describes chiral fields, $`\psi `$ is quark fields and $`m`$ is the dynamical quark mass. $`f_\pi `$ is the pion decay constant and experimentally $`f_\pi 93\mathrm{M}\mathrm{e}\mathrm{V}`$. Since our concern is the tree-level pions and one-loop quarks according to the Hartree mean field approach, the kinetic term of the pion fields which gives a contribution to higher loops can be neglected. Due to the interaction between the valence quarks and the Dirac sea, soliton solutions appear as bound states of quarks in the background of self-consistent mean chiral field. $`N_c`$ valence quarks fill the each bound state to form a baryon. The baryon number is thus identified with the number of bound states filled by the valence quarks kahana .

For $`B=1`$ and $`2`$, the spherically symmetric soliton meissner ; reinhardt ; wakamatsu  and the axially symmetric soliton sawado  were found respectively. Upon quantization, the intermediate states of nucleon and deuteron between the Constituent Quark Model and Skyrme Model were obtained.

The vacuum functional in Eq.(1) can be integrated over the quark fields to obtain the effective action

$`S_{\mathrm{eff}}[U]`$ $`=`$ $`iN_c\mathrm{lndet}\left(i\overline{)}mU^{\gamma _5}\right)`$ (2)
$`=`$ $`{\displaystyle \frac{i}{2}}N_c\mathrm{Spln}D^{}D`$ (3)
where $`D=i\overline{)}mU^{\gamma _5}`$. This determinant is ultraviolet divergent and must be regularized. Using the proper-time regularization scheme, we can write

$`S_{\mathrm{eff}}^{\mathrm{reg}}[U]={\displaystyle \frac{i}{2}}N_c{\displaystyle _{1/\mathrm{\Lambda }^2}^{\mathrm{}}}{\displaystyle \frac{d\tau }{\tau }}\mathrm{Sp}\left(\mathrm{e}^{D^{}D\tau }\mathrm{e}^{D_0^{}D_0\tau }\right)={\displaystyle \frac{i}{2}}N_cT{\displaystyle _{\mathrm{}}^{\mathrm{}}}{\displaystyle \frac{d\omega }{2\pi }}{\displaystyle _{1/\mathrm{\Lambda }^2}^{\mathrm{}}}{\displaystyle \frac{d\tau }{\tau }}\mathrm{Sp}\left[\mathrm{e}^{\tau (H^2+\omega ^2)}\mathrm{e}^{\tau (H_0^2+\omega ^2)}\right]`$ (4)
where $`T`$ is the Euclidean time separation, $`\mathrm{\Lambda }`$ is a cut-off parameter evaluated by the condition that the derivative expansion of Eq.(2) reproduces the pion kinetic term with the correct coefficient $`i.e.`$

$`f_\pi ^2={\displaystyle \frac{N_cm^2}{4\pi ^2}}{\displaystyle _{1/\mathrm{\Lambda }^2}^{\mathrm{}}}{\displaystyle \frac{d\tau }{\tau }}\mathrm{e}^{\tau m^2},`$ (5)
and $`H`$ is the Dirac one-quark Hamiltonian defined by

$`H={\displaystyle \frac{\stackrel{}{\alpha }\stackrel{}{}}{i}}+\beta mU^{\gamma _5}.`$ (6)
$`D_0D(U=1)`$ and $`H_0H(U=1)`$ correspond to the vacuum sectors.

At $`T\mathrm{}`$, we have $`\mathrm{e}^{iS_{\mathrm{eff}}}\mathrm{e}^{iE_{sea}T}`$. Integrating over $`\omega `$ in (4) and constructing a complete set of eigenstates of $`H`$ with

$`H|\nu =ϵ_\nu |\nu ,H_0|\nu ^{(0)}=ϵ_\nu ^{(0)}|\nu ^{(0)},`$ (7)
one obtains the sea quark energy meissner2

$`E_{sea}[U]={\displaystyle \frac{1}{4\sqrt{\pi }}}N_c{\displaystyle _{1/\mathrm{\Lambda }^2}^{\mathrm{}}}{\displaystyle \frac{d\tau }{\tau ^{3/2}}}\left({\displaystyle \underset{\nu }{}}\mathrm{e}^{\tau ϵ_\nu ^2}{\displaystyle \underset{\nu }{}}\mathrm{e}^{\tau ϵ_\nu ^{(0)2}}\right).`$ (8)
In the Hartree picture, the baryon states are the quarks occupying all negative Dirac sea and valence levels. Hence, if we define the total soliton energy $`E_{total}`$, the valence quark energy should be added;

$`E_{total}[U]=N_c{\displaystyle \underset{i}{}}E_{val}^{(i)}[U]+E_{sea}[U].`$ (9)
where $`E_{val}^i`$ is the valence quark contribution to the $`i`$ th baryon.

The baryon density $`b_0`$ for the baryon number $`B`$ soliton is defined by the zeroth component of the baryon current reinhardt ;

$`b_0`$ $`=`$ $`{\displaystyle \frac{1}{N_cB}}\overline{\psi }\gamma _0\psi ={\displaystyle \frac{1}{N_cB}}\left[{\displaystyle \underset{\nu }{}}\left(n_\nu \theta (ϵ_\nu )+\mathrm{sign}(ϵ_\nu )𝒩(ϵ_\nu )\right)\nu |\stackrel{}{r}\stackrel{}{r}|\nu {\displaystyle \underset{\nu }{}}\mathrm{sign}(ϵ_\nu ^{(0)})𝒩(ϵ_\nu ^{(0)})\nu |\stackrel{}{r}^{(0)}\stackrel{}{r}|\nu ^{(0)}\right]`$ (10)
where

$`𝒩(ϵ_\nu )={\displaystyle \frac{1}{\sqrt{4\pi }}}\mathrm{\Gamma }({\displaystyle \frac{1}{2}},\left({\displaystyle \frac{ϵ_\nu }{\mathrm{\Lambda }}}\right)^2)`$
and $`n_\nu `$ is the valence quark occupation number.

For $`B=3`$, it is expected the solution to have a tetrahedral symmetry from the study of the Skyrme model braaten . Therefore, we shall impose the same symmetry on the chiral fields using the rational map ansatz. According to the ansatz, the chiral field can be expressed as manton

$`U(r,z)=\mathrm{exp}\left(if(r)\stackrel{}{n}_R\stackrel{}{\tau }\right)`$ (11)
where

$`\stackrel{}{n}_R={\displaystyle \frac{1}{1+|R(z)|^2}}(2\mathrm{R}\mathrm{e}[R(z)],2\mathrm{I}\mathrm{m}[R(z)],1|R(z)|^\mathit{2})`$
and $`R(z)`$ is a rational map.

Rational maps are maps from $`CP(1)`$ to $`CP(1)`$ (equivalently, from $`S^2`$ to $`S^2`$) classified by winding number. It was shown in manton  that $`B=N`$ skyrmions can be well-approximated by rational maps with winding number $`N`$. The rational map with winding number $`N`$ possesses $`(2N+1)`$ complex parameters whose values can be determined by imposing the symmetry of the skyrmion. Thus, $`B=3`$ rational map with tetrahedral symmetry takes the form;

$`R(z)={\displaystyle \frac{\sqrt{3}iz^21}{z(z^2\sqrt{3}i)}}`$ (12)
where the complex coordinate $`z`$ on $`CP(1)`$ is identified with the polar coordinates on $`S^2`$ by $`z=\mathrm{tan}(\theta /2)\mathrm{e}^{i\phi }`$ via stereographic projection. Substituting (12) into (11), one obtains the complete form of the chiral fields with tetrahedral symmetry and winding number 3.

Apparently, the chiral fields in (11) takes a spherically symmetric form. Therefore one can apply the numerical technique developed for $`B=1`$ to find $`B=3`$ with tetrahedral symmetry wakamatsu .

Demanding that the total energy in (9) be stationary with respect to variation of the profile function $`f(r)`$,

$`{\displaystyle \frac{\delta }{\delta f(r)}}E_{total}=0,`$
yeilds the field equation

$`S(r)\mathrm{sin}f(r)=P(r)\mathrm{cos}f(r)`$ (13)
where

$`S(r)=N_c{\displaystyle \underset{\nu }{}}\left(n_\nu \theta (ϵ_\nu )+\mathrm{sign}(ϵ_\nu )𝒩(ϵ_\nu )\right)\nu |\gamma ^0\delta (|x|r)|\nu `$ (14)
$`P(r)=N_c{\displaystyle \underset{\nu }{}}\left(n_\nu \theta (ϵ_\nu )+\mathrm{sign}(ϵ_\nu )𝒩(ϵ_\nu )\right)\nu |i\gamma ^0\gamma ^5\stackrel{}{n}_R\stackrel{}{\tau }\delta (|x|r)|\nu .`$ (15)
The procedure to obtain the self-consistent solution of Eq.(13) is that $`1)`$ solve the eigenequation in (7) under an assumed initial profile function $`f_0(r)`$, $`2)`$ use the resultant eigenfunctions and eigenvalues to calculate $`S(r)`$ and $`P(r)`$, $`3)`$ solve Eq.(13) to obtain a new profile function, $`4)`$ repeat $`1)3)`$ until the self-consistency is attained.

For convenience, we shall take $`f_0(r)=\pi \mathrm{e}^{r/X}`$. To solve Eq.(7), we construct the trial function using the Kahana-Ripka basis kahana ;

$`\mathrm{\Psi }^\pm =\underset{K_{max}\mathrm{}}{lim}{\displaystyle \underset{i=1}{\overset{4}{}}}{\displaystyle \underset{K=0}{\overset{K_{max}}{}}}{\displaystyle \underset{M=K}{\overset{K}{}}}\alpha _{KM}^{(i)\pm }\phi _{KM}^{(i)\pm }(r,\theta ,\varphi )`$ (16)
where $`\mathrm{\Psi }^+`$ and $`\mathrm{\Psi }^{}`$ stand for parity $`(1)^K`$ and $`(1)^{K+1}`$ respectively, $`\phi `$ is the Kahana-Ripka basis and $`K`$ is the grand spin operator which is a good quantum number in the case of $`B=1`$ hedgehog. The basis is discretized by imposing an appropriate boundary condition for the radial wavefunctions at the radius $`r_{max}`$ chosen to be sufficiently larger than the soliton size. And also, the basis is made finite by including only those states with the momentum $`k`$ as $`k<k_{max}`$. The results should be, however, independent on $`r_{max}`$ and $`k_{max}`$.

According to the Rayleigh-Ritz variational method bransden , the upper bound of the spectrum can be obtained from the secular equation for each parity;

$`\mathrm{det}\left(𝐀^\pm ϵ𝐁^\pm \right)=0`$ (17)
where

$`𝐀^\pm (K^{}M^{},KM)`$ $`=`$ $`{\displaystyle \underset{ij=1}{\overset{4}{}}}{\displaystyle d^3x\phi _{K^{}M^{}}^{(i)\pm }H\phi _{KM}^{(j)\pm }}`$
$`𝐁^\pm (K^{}M^{},KM)`$ $`=`$ $`{\displaystyle \underset{ij=1}{\overset{4}{}}}{\displaystyle d^3x\phi _{K^{}M^{}}^{(i)\pm }\phi _{KM}^{(j)\pm }}.`$
For $`K\mathrm{}`$, the spectrum $`ϵ`$ becomes exact. Eq.(17) can be solved numerically.

Since the chiral field in Eq.(11) is less symmetric than the $`B=1`$ hedgehog, the hamiltonian has no grand spin symmetry. As a result, the states with different grand spin couple strongly and level splittings within the $`K`$ blocks occur. In Table 1 we present the schematic picture of the matrix elements $`𝐀(K^{}M^{},KM)`$. Although the size of the matrix $`𝐀(K^{}M^{},KM)`$ becomes quite large, due to the symmetry of the chiral fields, the functional space can be rearranged to reduce the size. Consequently, the space is divided with four blocks for each parity.

Fig. 1 shows the spectrum of the quark orbits as a function of the soliton-size parameter $`X`$. The $`P^+`$ orbit diving into the negative energy region is triply degenerate. As discussed in kahana , baryon number of the soliton equals to the number of diving levels occupied by $`N_c`$ valence quarks. Thus putting $`N_c=3`$ valence quarks on each of the degenerate levels, one obtains the $`B=3`$ soliton solution.

Fig. 2 shows the corresponding baryon density. As can be seen, it is tetrahedral in shape. Therefore, we confirm that the lowest lying $`B=3`$ configuration is tetrahedrally symmetric. This result is consistent with the $`B=3`$ skyrmion obtained by Braaten et.al braaten .

Fig. 3 shows the self-consistent profile function. For the total energy of the solution we obtain $`E_{total}=3596`$ MeV which is almost comparable to three times of the $`B=1`$ mass, and root mean square radius is $`\sqrt{r^2}0.6`$ fm. Our soliton seems to be tight object. This is mainly due to the missing of higher components of $`K`$ in our calculation. Their contribution becomes significant near the surface of the soliton and hence inclusion of the higher components will improve the size of the soliton.

Finally, we would like to mention that our result verifies the validity of the rational map ansatz for the Chiral Quark Soliton Model. The numerical technique used here is quite general and its application to $`B4`$ will be straightforward.

Acknowledgements

We are grateful to N.S.Manton for encouraging us to work on this subject and useful comments. We also thank V.B.Kopeliovich for suggesting this topic.


text:
# Enhancement of spontaneous emission in a quantum well by resonant surface plasmon coupling

## Abstract

Using time-resolved photoluminescence measurements, the recombination rate in an In0.18Ga0.82N/GaN quantum well (QW) is shown to be greatly enhanced when spontaneous emission is resonantly coupled to a silver surface plasmon. The rate of enhanced spontaneous emission into the surface plasmon was as much as 92 times faster than normal QW spontaneous emission. A calculation, based on Fermi’s golden rule, reveals the enhancement is very sensitive to silver thickness and indicates even greater enhancements are possible for QWs placed closer to the surface metal coating.

The spontaneous emission (SE) decay constant $`\tau `$ for radiating dipoles at $`\stackrel{}{r}_e`$ is given by Fermi’s golden rule

$$\frac{1}{\tau }=\frac{2\pi }{\mathrm{}}|f|\stackrel{}{d}\stackrel{}{E}(\stackrel{}{r}_e)|i|^2\rho (\mathrm{}\omega ),$$
(1)
where $`\rho (\mathrm{}\omega )`$ is the photon density of states (DOS) and $`f|\stackrel{}{d}\stackrel{}{E}(\stackrel{}{r}_e)|i`$ is the dipole emission matrix element. As pointed out by Purcell, SE may be enhanced by altering the photon DOSPurcellPR69 . For example, the ratio of enhanced to free space emission (the Purcell factor $`F`$) has been measured as large as 5 in an atomic system by placing the radiating atoms in a high $`Q`$, low volume cavityGoldsteinSpontaneous ; DutraPRA53 . A Purcell factor of up to 6 has been observed from quantum well (QW) and quantum dot emitters in vertical cavity surface emitting laser structures, while an enhancement of 15 has been observed from quantum dots in a microdisk cavityWeisbuchJL85 ; GerardPRL81 . Photonic crystals and distributed Bragg gratings have also been used to enhance the SE rate by as much as a factor of 4.5VuckovicJQE36 ; BoroditskyJLT17 ; BabaJLT17 . Such enhanced SE rates, achieved by increasing the photonic DOS in a small cavity, permit lower threshold, higher modulation frequency lasers as well as more efficient light emitting diodes.

The SE rate can also be modified when semiconductor or dye emitters are coupled to a surface plasmon (SP) of a metallic film HeckerAPL75 ; BabaJJAP35 ; GontijoPRB60 ; BrillianteJMS79 . A single QW can experience strong quantum electrodynamic coupling to a SP mode if placed within the SP fringing field penetration depth. An electron-hole pair in the QW recombines and emits a photon into a SP mode instead of into free space. The degree of SE rate modification for a given wavelength depends on the SP DOS at that wavelength. The strongest enhancement occurs near the asymptotic limit of the SP dispersion branch, the SP ”resonance” energy $`E_{sp}`$, where the SP DOS is very high. Non-resonant, SP-mediated SE enhancements as large as 6 have been observed from GaAs QWs near thin Ag filmsHeckerAPL75 . Even greater enhancements are possible for wide bandgap semiconductors whose emission wavelength is coincident with $`E_{sp}`$. In this report, time-resolved photoluminescence (TRPL) measurements of a partially silver-coated InGaN/GaN QW directly demonstrate the SP-mediated resonant enhancement of the SE rate for the first time in a semiconductor QW.

An InGaN/GaN QW was used in these experiments, grown by metal-organic chemical vapor deposition (MOCVD) on sapphire substrateKellerJCG195 . Over a 1.5 $`\mu `$m Si-doped GaN buffer layer was grown a 28 nm In0.04Ga0.96N reference layer, a 6 nm GaN layer, and the 3 nm In0.18Ga0.82N QW as shown in Fig. 1. Above the QW was a 12 nm Si-doped GaN cap layer, placing the QW within the fringing field depth of the SP. A layer of silver, $``$ 8 nm thick, was deposited by electron beam evaporation on one-half of the sample surface. The other half was left bare to facilitate direct comparison of the silvered and unsilvered results.

The bulk plasmon energy of silver is 3.76 eV, but the SP energy of Ag is lowered by the GaN dielectric constantEhrenreichPR128 . The SP dispersion relation is derived, using Maxwell’s equations, from the known dielectric properties of Ag and GaNBrunnerJAP82 . Considering a silver film of thickness $`t`$ and permittivity $`\epsilon _2`$, sandwiched between GaN and air with permittivities $`\epsilon _1`$ and $`\epsilon _3`$ respectively, the boundary condition gives the SP dispersion relation

$$(\frac{\gamma _1}{\epsilon _1}+\frac{\gamma _2}{\epsilon _2})(\frac{\gamma _3}{\epsilon _3}+\frac{\gamma _2}{\epsilon _2})(\frac{\gamma _1}{\epsilon _1}\frac{\gamma _2}{\epsilon _2})(\frac{\gamma _3}{\epsilon _3}\frac{\gamma _2}{\epsilon _2})e^{2\gamma _2t}=0,$$
(2)
where $`\gamma _i=k^2\epsilon _i\omega ^2/c^2`$ , $`i=1,2,3`$, and $`k=2\pi /\lambda `$ is the wavevector. The SP dispersion contains tangential and normal mode branches (Fig. 1), indicating the dominant direction of current flow in the silver film. For silver films with $`t`$ 8 nm, the tangential SP branch asymptotically approaches $`E_{sp}=`$ 2.85 eV ($`\lambda _{sp}=`$ 436 nm), the SP ”resonance” energy. Because the photon DOS is proportional to $`dk/d\omega `$, the SP DOS and SE enhancement will be greatest at the SP resonance.

The resonant enhancement is measured by comparing the luminescence decay rate from the photoexcited QW on the silvered and unsilvered sides. Room temperature TRPL measurements were performed using a 100 MHz Kerr-lens mode-locked, frequency-doubled Ti:Sapphire(Ti:S) laser with average incident pump power of 10 mW ($``$ 13 $`\mu `$J/cm2). The pump excitation energy (3.14 eV or 395 nm) was chosen to be below the bandgap of the InGaN reference layer and GaN layers so electron-hole pairs were generated only in the QW. The luminescence signal was dispersed in a grating spectrometer (600 gr/mm) and measured simultaneously across three wavelength bands (2.55-2.68, 2.71-2.87, and 2.79-2.96 eV) using a Hamamatsu streak camera with a resolution of 15 ps. The three adjacent 25 nm ($``$ 150 meV) wide wavelength bands spanned the entire continuous wave photoluminescence (cw PL) emitted from the QW. Features narrower than 25 nm were resolved by sequentially comparing adjacent 5 or 10 nm data windows offset from each other in 1 - 4 nm steps. Of course, all TRPL traces represent the sum effect of components with wavelength dependent behavior. However, narrowing the bandwidth of the windows further did not significantly improve the ability to measure wavelength dependence because of the reduced signal-to-noise, especially on the silvered side.

An example of a TRPL trace is presented in Fig. 2a, comparing the temporal decay of the unsilvered and silvered QW PL for a 25 nm-wide wavelength band. The luminescence decay constants from the silvered and unsilvered sides were independently derived from exponential fits to the data (5 or 10 nm windows) and then compared under identical pump and detector parameters. The uncertainty in fitted rate constants for a 10 nm window is 2 - 5 % at wavelengths near the peak cw PL emission (2.75 eV) and rises to as much as 10 - 16 % at longer and shorter wavelength extremes where the emission is weaker.

On the unsilvered side, the QW exhibited a long single exponential decay whose decay constant $`\tau _0`$ was slowest ($`\tau _0=`$ 28 ns) at wavelengths near the peak PL emission (2.76 eV) (Fig. 2b). This very long decay constant, and the correspondingly high PL intensity, indicate a high internal quantum efficiency and insignificant non-radiative (nr) processes compared to radiative (r) recombination ($`1/\tau _0=1/\tau _{nr}+1/\tau _r1/\tau _r`$) MonemarSPIE3624 ; ChichibuCh5 . Away from the peak PL emission wavelength, the recombination rate accelerates ($`\tau _0=`$ 4 - 5 ns) at the longest and shortest wavelengths measured. The wavelength dependence of QW emission is well understood from earlier TRPL studies of InGaN QWs ChichibuCh5  but differs markedly from the SE rate of a dipole radiator (dipole moment $`d`$) in a dielectric medium (index of refraction $`n`$)

$$\frac{1}{\tau _r(\omega )}=\frac{4nd^2\omega ^3}{3\mathrm{}c^3},$$
(3)
especially at low frequency where shallow level traps, impurities, and quantum dot-like structures in the QW contribute to recombination. A phenomenological estimate is a Lorentzian

$$\frac{1}{\tau _0(\omega )}\frac{1}{\tau _r(\omega )}=\frac{1}{\tau _{r_0}}\frac{\mathrm{\Delta }\omega ^2}{(\omega \omega _c)^2+\mathrm{\Delta }\omega ^2},$$
(4)
whose peak at ($`\omega _c`$) and linewidth $`\mathrm{\Delta }\omega `$ also roughly coincide with that of the measured cw PL (Fig. 2b).

By contrast, the weaker PL intensity through the silver-coated surface exhibits a bi-exponential decay. The faster decay component, with decay constant $`\tau _1`$, corresponds to enhanced recombination mediated by the SP mode ($`1/\tau _1=1/\tau _{nr}+1/\tau _r+1/\tau _{sp}1/\tau _{sp}`$). As expected, this faster component is strongest near $`E_{sp}`$ (2.85 eV), but it is evident at energies 200meV lower than $`E_{sp}`$ and decays with an almost wavelength-independent constant of $`\tau _1\tau _{sp}`$ 235 ps. Evidently, the SP resonance is fairly broad. The slower TRPL relaxation component has a decay constant ($`\tau _2`$) which decreases with decreasing wavelength from 5 ns at 2.94 eV to 9 ns at 2.61 eV. This component is the remnant radiative relaxation component of the QW ($`\tau _r`$). The frequency dependence is approximately $`\omega ^3`$ as in (3), and the role of other nonradiative processes is not evident.

To summarize the wavelength dependence of the SP enhancement, Fig. 3 plots the Purcell factor ($`F_{sp}=1+\tau _0/\tau _1`$) derived from the measured TRPL decay constants. The data demonstrate a sudden rise at higher energies, a peak enhancement at 2.8 eV, and weaker enhancement at lower energies. The maximum values of $`F_{sp}`$ were 36 (at 2.83 eV) and 92 (at 2.79 eV) at sample positions $`P_1`$ and $`P_2`$, respectively. By comparing the relative cw PL intensities of the InGaN reference layers, it is evident that the silver film is thicker at $`P_2`$ than at $`P_1`$. The differing $`F_{sp}`$ values at $`P_1`$ and $`P_2`$ consequently suggests a sensitivity of $`F_{sp}`$ to Ag thickness.

Fermi’s golden rule (1) provides insight into the frequency dependence of the Purcell factor and reveals the sensitivity of $`F_{sp}`$ to the silver thickness $`t`$ and Ag-QW separation $`a`$GontijoPRB60 . First, the electric field of the SP mode at the QW must be calculated and used to derive the dipole matrix element. The SP electric field varies only in the $`z`$ direction, so the normalization of $`E(z)`$ to a half quantum of zero-point fluctuation in the dispersive medium becomes

$$\alpha ^2=\frac{S}{A}=\frac{\mathrm{}\omega /2}{\frac{A}{8\pi }_{\mathrm{}}^{\mathrm{}}𝑑z\frac{(\omega \epsilon (\omega ,z))}{\omega }|E(z)|^2},$$
(5)
where $`E(z)`$ is the unnormalized electric field at a distance $`z`$ from the Ag-GaN interface, $`|\alpha E(a)|^2`$ is the normalized electric field at QW depth $`a`$, $`A`$ is the quantization area, and $`\epsilon (\omega ,z`$) is the dielectric function of the GaN, Ag, or air. The enhanced recombination rate ($`1/\tau _{sp}`$) can then be estimated in the QW under the influence of the local electric field from the tangential SP mode

$`{\displaystyle \frac{1}{\tau _{sp}(\omega )}}={\displaystyle \frac{2\pi }{\mathrm{}}}({\displaystyle \frac{1}{3}}d^2|\alpha E(a)|^2)2\pi k{\displaystyle \frac{A}{(2\pi )^2}}{\displaystyle \frac{dk}{d(\mathrm{}\omega )}}`$
$`={\displaystyle \frac{S}{3\mathrm{}^2}}|dE(a)|^2k{\displaystyle \frac{dk}{d\omega }},`$ (6)
where the factor $`1/3`$ comes from spatial averaging of polarization. Inserting measured values for $`a`$ (12 nm) and $`t`$ (8 nm) while using the $`\tau _1`$ data to fit for $`d`$ (24 nm), the calculation confirms that $`\tau _{sp}`$ is relatively independent of frequency from 2.6 eV to $`E_{sp}`$.

The frequency dependence of $`F_{sp}`$ derives from (6) and the frequency dependence of the unsilvered QW recombination rate $`\tau _0`$ (Fig. 2b). For the parameters in this experiment, the frequency dependence of $`F_{sp}`$ below $`E_{sp}`$ derives primarily from $`\tau _0`$. Unfortunately, a predictive theory of the $`\tau _0`$ is not available, so an accurate estimate of $`F_{sp}`$ is not possible. However, if (4) and (6) are used as approximations for $`\tau _0`$ and $`\tau _1`$, respectively, then $`F_{sp}(\omega )1+\frac{\tau _0}{\tau _1}`$ may be plotted (Fig. 3) using the parameters for this sample. The predicted frequency peak and width of $`F_{sp}`$ agree well with the measured peaks and widths of the TRPL data. The peak value of $`F_{sp}`$ is a sensitive function of Ag thickness through $`\tau _{sp}`$. A reduction of $`t`$ from 8 to 6 nm halves the peak value of $`F_{sp}`$, suggesting that the differing values of $`F_{sp}`$ observed at $`P_1`$ and $`P_2`$ arise from small variations of the Ag thickness.

Recently, $`F_{sp}`$ was measured in a similar sample by means of cw PL using a HeCd excitation source ($`E_{ex}=`$ 3.82 eV, $`\lambda _{ex}=`$ 325 nm) GontijoPRB60 . A ratio of the PL emission for the unsilvered to silvered sides was measured as a function of wavelength, and a peak enhancement factor of 55 was estimated from the data. Those measurements were repeated here, and although the measured enhancement factors were similar to and consistent with the TRPL data, the energy of maximum enhancement in both cw PL measurements was blue shifted from $`E_{sp}`$. When time-integrated PL is measured using the Ti:S laser instead, this blue shift disappears (inset, Fig. 3). Note that the HeCd laser excites all layers of the sample while the Ti:S laser only excites the QW. Therefore, any cw PL ratio measured using a HeCd laser must be understood to represent $`F_{sp}`$ convolved with other excitation-dependent effects (particularly in the cap layer), while the TRPL data measures only $`F_{sp}`$. Furthermore, the role of nonradiative recombination was ignored in that work, leading to an overly simiplistic estimate of $`F_{sp}`$.

Using (4) and (6) to estimate $`F_{sp}`$, even larger enhancements are predicted over narrower frequency bands. For a prominent resonance in $`F_{sp}`$, Fig. 4 indicates that it is necessary for the Ag film to possess a thickness $``$ 6 nm for a QW 12 nm from the surface. In support of the earlier deduction that the Ag film is slightly thicker at $`P_2`$ than $`P_1`$, the value of $`F_{sp}`$ is predicted to increase with increasing film thickness and asymptotically approaches 422 for $`t`$ 26 nm ($`a=`$12 nm).

The strength and frequency dependence of $`F_{sp}`$ are even more sensitive functions of the Ag-QW separation, especially for small $`a`$. Resonant enhancements of more than $`10^4`$ are predicted for QWs only 4 nm below the surface. These predicted enhancements may actually be conservative because the synergistic ”back action” coupling between dipole emitters and SP field, which increases as $`a`$ decreases, is not included in this calculation. A more comprehensive calculation of this enhancement factor, including the necessary radiation reaction effects, is beyond the scope of this paper. Nevertheless, the enhancement will likely remain broad because inhomogeneous broadening ($`\mathrm{}\mathrm{\Delta }\omega _{inh}`$ 100 meV) will probably limit $`\frac{\omega _c}{\mathrm{\Delta }\omega }<`$ 30.

###### Acknowledgements.

The authors would like to thank S. Kellar, U. Mishra, and S. DenBaars of UCSB for providing the QW sample, and I. Gontijo, M. Boroditsky, and G. Khitrova for useful discussions. This work was supported by the Army Research Office and the National Research Council.


text:
# Supernova 1998bw - The final phases Based on observations collected at the European Southern Observatory, La Silla and Paranal, Chile and on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555.

## 1 Introduction

Supernova (SN) 1998bw was discovered in April 1998 in a spiral arm of the face-on SB galaxy ESO~184-G82 at a redshift of $`z=0.0085`$ (Galama et al. 1998a). This exceptional type Ic supernova was born famous because of its positional and temporal coincidence with the gamma-ray burst GRB~980425 (Galama et al. 1998b). This association was further supported by the very early onset of the bright radio luminosity (Kulkarni et al. 1998). Also in the optical, the supernova was unusually luminous and showed very high expansion velocities (Patat et al. 2001).

The earliest photometric evolution of SN~1998bw was presented in Galama et al. (1998b). McKenzie & Schaefer (1999) monitored the early exponential decay, and the late phases were added by Sollerman et al. (2000). Patat et al. (2001) combined these data to make a rough bolometric light curve, which has been the basis of several model investigations (Nomoto et al. 2001; Nakamura et al. 2001a). The aim of this paper is to improve and extend this light curve to the final phases of SN~1998bw.

Modelling of the early light curve and spectra suggested an extremely energetic explosion \[$`(25)\times 10^{52}`$ erg\] of a massive star, composed mainly of carbon and oxygen (Iwamoto et al. 1998; Woosley et al. 1999). As this kinetic energy is more than ten times that of a canonical core-collapse supernova, the term hypernova was suggested for this event. The mass of $`{}_{}{}^{56}\mathrm{Ni}_{}^{}`$ needed to power the early light curve in these early models, 0.5–0.7 $`\mathrm{M}_{}`$, is much larger than the $`0.1`$ $`\mathrm{M}_{}`$ typical for “normal” core-collapse SNe (e.g., Sollerman 2002).

This large mass of $`{}_{}{}^{56}\mathrm{Ni}_{}^{}`$ was disputed by Höflich et al. (1999), who suggested that the early light curve could be reproduced by a normal SN Ic, given the right viewing angle and degree of asymmetry. The modelling of Sollerman et al. (2000) indicated, however, that at least 0.3 $`\mathrm{M}_{}`$ of $`{}_{}{}^{56}\mathrm{Ni}_{}^{}`$ was required to power the late emission of the supernova. This conclusion was later supported by Nakamura et al. (2001a), see also Mazzali et al. (2001).

In this paper we focus on the very late light curve of SN~1998bw. In Sect. 2 we present new VLT observations and the technique used to achieve the revised optical light curves. Here we also present some late near-IR photometry. In Sect. 3 we present the very late *Hubble Space Telescope* (*HST*) observations of SN~1998bw. The results are presented in Sect. 4. After a discussion in Sect. 5, we briefly summarize our conclusions in Sect. 6.

## 2 Ground Based Observations

### 2.1 Late-Time Optical Photometry

The photometry discussed here is based on the late observations presented in Sollerman et al. (2000; hereafter S00). S00 obtained data at eight epochs between 32 and 503 days past explosion, using the ESO 3.6m telescope on La Silla and the Very Large Telescope (VLT) on Paranal. Here we assume the date of explosion to be 980425.9, the moment of detection of GRB~980425 (Soffitta et al. 1998).

To the S00 data set we have added the late time observations presented in Patat et al. (2001; hereafter P01). These data were primarily obtained using the ESO 3.6m telescope and are added in particular to include the $`B`$-band photometry. We have also added previously unpublished photometry obtained at the VLT. The details of the ground-based optical observations are shown in Table 1.

As discussed in S00 and P01, the photometry of SN~1998bw is fairly complicated. This is because the supernova is superposed on an H ii region, and sits in a rather complex region of the host galaxy. This was strikingly revealed by the *HST*/STIS image published by Fynbo et al. (2000), where several objects can be seen within $`1\mathrm{}`$ of the position of the supernova. This is here shown in detail in Fig. 1.

This means that even sophisticated point-spread function (PSF) fitting photometry may not be adequate to accurately measure the SN brightness at late phases. S00 used DAOPhot to estimate the errors due to background subtraction, while P01 used the SNOOPY software. Although these methods were consistent, the very complex background of SN~1998bw as observed with *HST* called for a revision of the light curve. We decided to try another method to establish the late time light curve of this important supernova.

We therefore obtained deep $`B`$, $`V`$, $`R`$, and $`I`$ imaging of the region containing the supernova using FORS1 on the VLT. These observations were obtained in service mode under excellent atmospheric conditions in early October and in early November 2000, some 900 days after the supernova explosion. In each filter band 6–9 frames were coadded to achieve a deep combined image. The $`B`$ and $`V`$ frames have a total exposure time of 90 minutes each, while the $`R`$ and $`I`$ frames have 120 minutes. All combined frames have good image quality, with a FWHM less than one arcsecond.

In these images, no sign of the supernova could be detected, i.e., the background on the position of the supernova appeared smooth. These frames were therefore used as template frames to subtract the background from all the previous images (see e.g., Filippenko et al. 1986 for an early account of this method). The detailed procedure used is outlined in Schmidt et al. (1998), and we have used the same software as used by the High-$`z`$ Supernova Search Team for galaxy subtraction, magnitude measurements and error estimates.

After the usual bias subtraction and flat fielding, all images were aligned with the template image in the same filter band. The image with the best seeing was then convolved to match the seeing of the other image. Thereafter the image was scaled to the same intensity as the template, and the template was subtracted from the image in a region around the supernova. The program DoPhot (Schechter et al. 1993) was used to measure the magnitudes of the supernova and of the local standard stars in the subtracted image.

A PSF was constructed from stars in the field using DAOPhot (Stetson 1987). It was used to add ten artificial stars with the same flux as the supernova to user-specified locations on the original frame. The magnitudes of these stars were then retrieved using the same technique as for the supernova. This method was used to assess the accuracy of our measurements, and the standard deviations in the differences in the retrieved magnitudes of the artificial stars were adopted as the errors on the supernova magnitudes in Table 2.

In all cases, the success of the subtraction has been manually checked. Indeed, very clean results were achieved, where only the supernova remains on the subtracted frame. This is due to the fact that many of these observations have been obtained with the same telescope/instrument/CCD setup. Similar seeing conditions are also important, and in particular a high quality template frame is crucial.

In Fig. 2 we show an example of the region of the supernova in our template frame, and in our frame from day 503. This was the last epoch used in S00 and P01. Also shown is the subtraction, which clearly shows the supernova.

The instrumental magnitudes were finally converted using local standards measured by Galama et al. (1998b). This worked well for the $`B`$-band, where four standard stars showed a standard deviation of $`0.03`$ magnitudes. Most of this scatter is in fact due to the neglect of color terms, which, however, do not affect the supernova as we chose a primary standard with similar color as the SN itself. For the other bands, all Galama standards were saturated in the long exposures required at these late phases. We then established secondary, faint local standards. The magnitudes of four of these local standards are presented in Table 3, together with their angular offsets from the supernova. The standard deviations for these stars measured over all the epochs were never higher than 0.03 magnitudes. Adding this to the uncertainty reported by Galama et al. for their standards, we consider the accuracy of the secondary standards to be about 0.06 magnitudes.

### 2.2 Late Upper Limits

At our last epoch of VLT observations prior to the template images, i.e., at 706 days past explosion, the supernova is not detected in the subtracted images. We have used these difference images to estimate upper limits on the supernova magnitude in the following way. We measured the fluxes in a number of apertures located next to the supernova. From these measurements we estimated the standard deviations in a 5 pixel radius aperture, corresponding to the image FWHM. This estimate was then used to acquire $`3\sigma `$ upper limits within the aperture. These limits are 25.1, 25.1, 25.1, and 24.7 in $`B`$, $`V`$, $`R`$, and $`I`$, respectively. This procedure is more conservative than a simple measurement of the statistics in the region of the supernova. We also verified that a $`R=25.1`$ supernova would indeed have been detected by adding such a star to the original frame, and recovering the magnitude in the subtracted frame. This would thus provide conservative upper limits on the flux of the supernova at this epoch. The main assumption is that the supernova was indeed absent in the template frame. This is further elaborated below.

### 2.3 The Near-Infrared

Here, we present the last infrared (IR) observations of SN~1998bw of which we are aware.

SN~1998bw was observed using ISAAC on the VLT on 1999 May 1, 370 days after the explosion. The observations were obtained in Autojittermode, offsetting the telescope within a $`60\mathrm{}`$ wide box. Five frames were obtained in the $`J`$-band (1.11–1.39 $`\mu `$m), each one composed of a stack of five 30 second exposures. The total exposure time in the $`J`$-band was thus 750 seconds.

In the $`H`$-band (1.50–1.80 $`\mu `$m) five frames composed of five 20 seconds exposures yielded a total exposure time of 500 seconds. In addition, an infrared spectrum was obtained at this epoch, but the signal-to-noise ratio of this spectrum is low, and we will not discuss it further. Instead, we will use the two acquisition images to determine the IR contribution at these late phases.

The reductions were done using IRAF111Image Reduction and Analysis Facility (IRAF), a software system distributed by the National Optical Astronomy Observatories (NOAO).. From each image we subtracted a mean sky, constructed from combining all the other frames with suitable clipping to remove objects. Flat field corrections were applied using dome flat-fields. The images were then stacked and combined. The FWHM in the combined frames were $`0\stackrel{}{.}4`$ and $`0\stackrel{}{.}55`$ for $`J`$ and $`H`$, respectively.

For the calibration of the photometry, we used the zero points determined for this specific night at the VLT, from observations of photometric standards. To check the stability of the zero point, we also re-measured three local standard stars from the three earlier epochs obtained with NTT/SOFI (P01). The standard deviation in the standard stars on the three SOFI nights were 0.04 and 0.02 magnitudes in $`J`$ and $`H`$, respectively. The magnitudes agreed with the ISAAC measurements to better than 0.05 magnitudes, with a scatter of 0.02 mag.

In the combined frames, the SN magnitude were measured with PSF photometry using DAOPhot within IRAF. The measured magnitudes are $`J=20.16`$ and $`H=19.54`$. The aperture correction was obtained from the three local standards, with a negligible scatter. The DAOPhot measurement error was 0.03 ($`J`$) and 0.06 ($`H`$) magnitudes. The total error budget for the photometry is therefore certainly better than 0.1 magnitudes in both bands. These results should be treated with some caution, however, as we know from the optical data that PSF fitting may not be the best technique for this supernova.

## 3 The Last Phases with HST

### 3.1 STIS Observations

In order to follow the supernova to very late times, and to get a better understanding of the supernova environment, we have obtained detailed imaging at four epochs using the Space Telescope Imaging Spectrograph (STIS) onboard the *HST*. A log of our HST observations is given in Table 4.

The observations on 2000 June 25, 2000 Nov. 21, and 2001 Aug. 28 were taken as part of the Supernova Intensive Study programme (P.I. Kirshner). The data from 2000 June 11 was taken as part of the Survey of Host Galaxies of Gamma-Ray Bursts (P.I. Holland, Holland et al. 2000). The total exposure times in the 50CCD (clear, CL) aperture were 1240, 5672, 4000, and 5674 seconds on the four epochs, respectively. The 50CCD aperture has a central wavelength of 5850 Å and a width of 4410 Å.

In addition to the CL frames we also obtained some images in the F20X50LP (long pass, LP) filter at two epochs. The total exposure times were 1185 and 1319 seconds, respectively. The F28x50LP aperture has a central wavelength of 7320 Å and a width of 2720 Å.

The 2000 June 11 images were taken using a four-point STIS-CCD-BOX dithering pattern with shifts of 2.5 pixels ($`=0\stackrel{}{.}127`$) between exposures. The other images were dithered using the same pattern with shifts of 5 pixels ($`=0\stackrel{}{.}254`$).

The data were preprocessed through the standard STIS pipeline and combined using the dither (v2.0) software (Fruchter & Hook 2002) as implemented in IRAF. We selected a final output scale of $`0\stackrel{}{.}0254`$/pixel, corresponding to exactly half the original size, and set the “pixfrac” parameter to 0.6. Since there were only two images taken using the LP filter on 2000 Nov. 21 we had to increase “pixfrac” to 1.0 for these data in order to fill all pixels in the output frame. This resulted in a reduction in the image quality of this data compared to the drizzled 2000 June 11 LP image. The six drizzled images were finally reregistered onto a common system.

### 3.2 Photometry

The supernova is located in a complex region of the host galaxy (Fig. 1). There are several objects within a radius of $`1\mathrm{}`$ from the supernova position, and SN~1998bw itself is located on top of an extended, filamentary structure. In order to perform photometry of the supernova without contamination from the underlying structure we again used the template subtraction method to subtract the 2001 Aug. 28 CL image from each of the other CL images. An example subtraction is given in Fig. 1.

We then constructed PSFs for each subtracted CL image using DAOPhot II/AllStar (Stetson 1987; Stetson & Harris 1988) and fit these to the image of the supernova on each subtracted frame. We computed aperture corrections for each frame using several isolated stars. The PSF magnitudes were corrected to an aperture with a radius of $`1\stackrel{}{.}108`$, which Table 14.3 of the STIS Instrument Handbook suggests contains 100% of the flux from a point source. These corrections were about 0.1 magnitudes at all epochs. The instrumental magnitudes were converted to the AB system using the zero points of Gardner et al. (2000).

In order to understand any systematic errors we performed a series of artificial star tests on each CL image. We added a set of stars with known magnitudes comparable to that of the supernova to locations on each frame where the underlying structure was similar to the underlying structure of SN~1998bw. We then subtracted the template image and performed photometry on these artificial stars in exactly the same way as for the supernova. The resulting median differences between the input and recovered magnitudes of the artificial stars were always below 0.2 magnitudes. We regard this as a conservative error on the difference photometry.

No template image was taken using the LP filter. Therefore, no reliable magnitudes could be obtained in this band. For the 2000 June 11 image we estimated the supernova magnitude by pure PSF fitting on the drizzled image. The estimated magnitude is LP$`{}_{\mathrm{PSF}}{}^{}=25.02\pm 0.07`$. Note that there was an error in the LP aperture corrections used by Fynbo et al. (2000). That paper also contains some color information on the surrounding objects, although a correction of -0.51 mag should be added to all those LP magnitudes. The PSF magnitude for the supernova will clearly also include some of the light from the underlying background, which will result in overestimating the flux from the supernova. As mentioned above, the 2000 Nov. 21 LP frame is less useful. Visual inspection of that frame does show the supernova, just as in the CL frame from the same epoch, but we were not able to perform a reliable PSF estimate of its magnitude, since the remaining PSF residuals were always large. Therefore, only the CL magnitudes will be discussed below.

## 4 Results

### 4.1 The Optical Light Curve from the Ground

The results of the ground based photometry are presented in Table 2, and shown in Fig. 3.

The light curves presented in Fig. 3 are all very smooth, with the exponential decay continuing all the way to the very last data points, some 540 days past explosion. The $`R`$ and $`I`$ curves, for example, are perfectly well fit by a linear decline of 1.5–1.6 magnitudes per 100 days. This is in contrast to the results published in S00 and P01, where the very late light curve appeared to flatten out. We interpret the difference as due to the complex background of SN~1998bw, and regard the new magnitudes achieved with template subtraction as more accurate. Having seen the very complex environment of the supernova in *HST* detail this may not be surprising. In fact, the total magnitude of the complex region surrounding the supernova amounts to $`m_\mathrm{V}21`$. In measuring supernovae some 2 magnitudes fainter than the background, even small errors in the PSF fitting will amount to substantial errors in the supernova magnitudes, an effect that will mimic a flattening of the late light curve. Core-collapse supernovae are known to explode in regions of star formation, and it is not unlikely that photometry of other supernovae at very late phases is hampered by similar effects (see Clocchiatti et al. 2001 for a similar discussion).

To investigate the total light emitted should try to build a bolometric light curve. We begin by constructing the $`L_{BVRI}`$ light curve, containing the optical emission from 4400 to 7900 Å. We have done this in the following way.

All magnitudes were converted to a monochromatic flux using the conversions in Bessell (1979). To add the missing $`B`$-magnitudes, a linear fit to these magnitudes was used, where the day 141 magnitude from McKenzie & Schaefer (1999) was included to constrain the early evolution. The light curve slope in the $`B`$-band is then the same as in the other bands. The fluxes were corrected for an extinction of $`E(BV)`$=0.06 mag (Schlegel et al. 1998; P01), and then the monochromatic fluxes were simply integrated from 4400 to 7900 Å. Finally we adopted a distance of 35 Mpc to convert to luminosity ($`z=0.0085`$, $`H_0=73`$ km s-1 Mpc-1). To the resulting $`L_{BVRI}`$ points we will add a constant IR contribution in Sect. 4.4 to arrive at the light curve presented in Fig. 4.

The slope of the late time decline is well fit by a linear decay in absolute magnitudes by 1.55 magnitudes per 100 days from 140 to 538 days past explosion. This is certainly faster than the decay rate of the radioactive $`{}_{}{}^{56}\mathrm{Co}_{}^{}`$ powering the supernova at this phase, again showing that most of the gamma-rays are slipping out of the ejecta at late phases (S00). We discuss this further in Sect. 5.

### 4.2 HST Observations

The photometry from the HST observations is presented in Table 5. The CL magnitudes are the AB magnitudes measured by DAOPHOT on the template subtracted frames and with aperture corrections applied (CLSub). The errors are simply the formal errors from DAOPhot added to the standard deviations in the aperture corrections. This may be a bit optimistic, and from the artificial star test presented in Sect. 3.2 we regard an additional systematic error of 0.2 magnitudes as a more reliable estimate of the errors. In Table 5 we have also included the magnitudes measured on the drizzled frames using DAOPhot PSF-fitting (CLPSF). These clearly overestimate the flux of the supernova by including also light from the underlying structure, emphasizing the need of template subtraction even at the resolution achieved by *HST*.

The very first result to note is, of course, the clear detection of the supernova in the subtracted frames (Fig. 1). Fynbo et al. (2000) used VLT imaging (S00) to establish the position of the supernova in their STIS images. Based on this, they suggested that SN~1998bw was identical to the object positioned close to the middle of this messy region. Our analysis has demonstrated that this object is indeed fading. This strongly reinforces the supernova identification and allows us to follow the supernova to very late phases.

Secondly, from the magnitudes in Table 5, and from Fig. 3, we note that the supernova light curve seems to be levelling out at the very late phases. The reason for this will be discussed in Sect. 5.

In order to test if the supernova was still present in the 2001 Aug. 28 CL image we used DaoPhot II/AllStar to fit a PSF to the light at the location of SN~1998bw. Since there was no template to subtract from this image, the PSF fit will be strongly biased by the underlying light at the location of the supernova. The aperture-corrected PSF magnitude is CL$`{}_{\mathrm{PSF}}{}^{}=27.4\pm 0.1`$. However, an examination of the residuals after subtracting the PSF for the supernova suggests that the flux from the supernova is contaminated by light from the underlying structure and by the object $`0\stackrel{}{.}056`$ south of the supernova. To estimate an upper limit we again performed an artificial star test by injecting stars of different magnitudes onto regions of similar background as the supernova position. For input magnitudes fainter than $`2930`$, we found that the recovered PSF magnitudes were always close to 27.4, as measured for the supernova position. Using the standard deviation of the recovered magnitudes as a one sigma error, we find the $`3\sigma `$ upper limit of about 28.5 for the supernova. This can be regarded as an upper limit on the emission from the supernova in the CL band 1221 days past explosion. This suggests that the supernova continued to fade well into 2001 at a rate similar to that observed between June and November 2000.

### 4.3 Combining the datasets

In measuring the ground based photometry, we made the implicit assumption that the supernova had completely disappeared in our last VLT template images. From the *HST* results we now know that the supernova was still there at this phase, although very faint. In fact, at a magnitude of $`26.7`$ on day 907 the supernova will not affect any of the ground based measurements by more than 0.05 magnitudes. We will make no correction for this effect.

However, our late VLT upper limit may indeed be significantly affected, by up to 0.25 magnitudes. We have taken this into account in the combined upper limit on the luminosity presented in Fig. 4.

Better spatial resolution clearly allows more accurate subtractions and thus photometry. It is clear, however, that even with *HST* there may be unresolved substructures contributing to the PSF. At 35 Mpc each drizzled pixel corresponds to 4.3 parsecs on the sky. In Table 5 we have included the magnitudes obtained with PSF-fitting on the drizzled frames, and it is clear that such an approach overestimates the flux in the supernova. The effect increases at later phases and would again make the light curve flatten. This effect should, however, be correctly accounted for using the template subtraction. We therefore believe that the flattening of the late light curve as observed by *HST* is a robust result.

There is again the issue of contamination from remaining supernova flux in the *HST* template frame. Our estimate of the upper limit on the supernova CL magnitude in the 2001 August frame is $`M_{\mathrm{CL}}^{\mathrm{AB}}28.5`$. If the supernova was really at this magnitude in the template frame, the supernova would be brighter than estimated in June and November 2000. The effect could amount to 0.1 magnitudes in June and some 0.2–0.3 magnitudes in November. This is an uncertainty we can not overcome without further deep and very late imaging. Here we only note that such a contamination will not dramatically affect the late light curve slope. It would only flatten the late light curve between 778 and 940 days from a slope of 0.6 magnitudes per 100 days to 0.5 magnitudes per 100 days.

It is also non-trivial to compare directly the photometry obtained in the very broad CL filter with the ground based data. Established conversions as those presented in Rejkuba et al. (2000) are based on stars, and the late supernova spectrum differs from any stellar spectrum. Simply integrating a flat spectrum over the FWHM is a crude approach for such a wide filter. Nakamura et al. (2001a) assumed that ($`M_{\mathrm{AB}}=`$) $`M_\mathrm{V}=M_{\mathrm{bol}}`$ to put the photometry of Fynbo et al. (2000) onto their “bolometric” light curve. Such an estimate differs by a factor 1.8 from the simple integration.

The correct approach would of course be to use the supernova spectra obtained at the same epochs for the corrections, but no such spectra are available. The latest spectrum from S00 is from 504 days past explosion. However, this spectrum is rather noisy and may be contaminated by the nearby objects revealed by HST.

We have taken the following approach to compare the different observations. From the broad band VLT photometry we do get a gross spectral energy distribution (SED) for the epochs up to 540 days past explosion. This distribution is seen to evolve fairly slowly with time and should not be influenced by the nearby objects. At all epochs, most of the energy ($`\nu F_\nu `$) emerges in the $`R`$-band. From the late-time spectroscopy in S00 we know that this is mainly due to the strong \[O i\] $`\lambda \lambda `$ 6300,6364 lines. We will assume that the same gross spectral distribution is appropriate at the late HST phases, with a fairly small contribution from the $`U`$-band (P01) and a continuation of the SED into the IR as observed on day 370.

We then used SYNPHOT to scale this SED to give the measured count rates for the CL filter function. The scaled SED can then be integrated, after a correction for extinction (Fitzpatrick 1999), and converted to luminosity as above. We noted that the count rate in the CL filter does not depend much on the assumed SED outside the $`L_{BVRI}`$ range. Assuming no flux outside 4000–8000 Å decreased the measured counts by only $`14`$%, as both the SED and the CL filter curve peak in the middle of this interval. To avoid assumptions on the non-observed SED regimes we therefore integrated the SED only over the 4000–8000 Å range, to make it directly comparable to the ground based data.

### 4.4 The OIR Light Curve

An uncertainty in the bolometric light curve of P01 was the contribution of the IR emission. They assumed that the fraction of the $`UVOIR`$ emission escaping in the IR was the same at late times as it was on day 65, the last day covered by their IR observations. Then the IR ($`JHK`$) contained $`35`$% of the supernova emission, and this constant fraction was added to the optical light curve at all later phases (P01).

It is true that for the well observed SN~1987A, the fraction of the energy emerging from the supernova in the $`V`$, $`R`$, and $`I`$ bands ($`L_{VRI}`$) was virtually constant ($`40`$%) up to 400 days after the plateau phase (Schmidt et al. 1993). However, the physics of a rapidly expanding, hydrogen-free ejecta will differ from the case of SN~1987A. As the radioactive heating decreases and the ejecta expand, the gas will cool down and a larger fraction of the emission may therefore be pushed into the IR. The assumption of constant IR contribution at late phases should therefore be checked.

To estimate the contribution in the IR we followed the same procedure as outlined above for $`L_{BVRI}`$. The optical emission was interpolated on day 370 and integrated from $`B`$ to $`I`$. The total emission was then calculated by including also the $`J`$ and $`H`$-bands in the integral, where the conversions to flux from Wilson et al. (1972) were used. The result is that at this epoch 42% of the emission emerges in the near-IR. The same exercise for the data of day 65 gives a fraction of 31% (comparing $`L_{BVRI}`$ to $`L_{BVRIJH}`$). At all epochs, most of the energy emerges in the $`R`$-band.

It is thus clear that the importance of the near-IR has increased with time, although the effect is not dramatic.

To complete our OIR (optical-infrared) light curve we have simply added the IR fraction observed at day 370 to all the dates on the light curve presented in Fig. 3. This assumes that the IR contribution does not evolve at the last *HST* phases.

## 5 Discussion

The final OIR light curve is admittedly and necessarily based on a number of assumptions, but we are confident that the basic features of the light curve, the steep decline followed by a late flattening, are real. In this section we will discuss a few possible interpretations for such an evolution.

We first note the steep exponential decline seen in the ground based observations. This fast decline appears to have started already about 65 days past the explosion (McKenzie & Schaefer 1999; S00), and the light curve continues to fall significantly steeper than the decay rate of $`{}_{}{}^{56}\mathrm{Co}_{}^{}`$ up to more than 500 days past explosion. There is thus no sign for a “positron phase” in which the fully deposited kinetic energy from the positrons would dominate the light curve. Two different scenarios could account for this.

One is that the optical depth for gamma-rays, albeit decreasing, may always be large enough to dominate over the positrons. An example is shown below (Fig. 5). Alternatively, the positrons may indeed dominate at late phases, but some of the positrons are able to escape the ejecta before being thermalized.

The second characteristic to note in the very late-time light curve of SN~1998bw is the flattening of the light curve at the last phases, as observed by *HST*. The late-time light curves of core-collapse supernovae can be powered by several mechanisms (see e.g., Sollerman et al. 2001), as will be further discussed below.

### 5.1 Radioactive Powering?

The contributions from different radioactive elements to late supernova light curves were discussed by Woosley et al. (1989). Only in SN~1987A has a radioactive decay other than ($`{}_{}{}^{56}\mathrm{Ni}_{}^{}`$$``$)$`{}_{}{}^{56}\mathrm{Co}_{}^{}`$$``$$`{}_{}{}^{56}\mathrm{Fe}_{}^{}`$ been unambiguously observed to power the late light curve. Also for this supernova, a flattening of the light curve was observed after about 800 days. The reason for this was actually two-fold. First, at these late phases the radioactive decay of the more long-lived nucleus $`{}_{}{}^{57}\mathrm{Co}_{}^{}`$ became important. Secondly, the long recombination scale allowed energy stored at earlier epochs to contribute at later phases (Clayton et al. 1992; Fransson & Kozma 1993). This is known as the “freeze-out” effect.

For SN~1998bw, with its very high explosion energy, the discussion must also take into account that the optical depth to gamma-rays decrease very rapidly. The high explosion energy also opens the possibility for different nucleosynthesis (Nakamura et al. 2001b).

The decreasing optical depth to gamma-rays is obvious from the first part of the light curve, as mentioned above. A very simplistic model to account for this in terms of radioactive decay is presented in Fig. 5. In this model the flux from the decay of $`{}_{}{}^{56}\mathrm{Co}_{}^{}`$ evolves as $`\mathrm{e}^{t/111.3}\times \left(10.965\mathrm{e}^\tau \right)`$, where the optical depth, $`\tau =(t_1/t)^2`$, decreases due to the homologous expansion. Here, 111.3 days is the decay time of $`{}_{}{}^{56}\mathrm{Co}_{}^{}`$ and $`t_1`$ sets the time when the optical depth to gamma-rays is unity. Furthermore, 3.5% of the energy in these decays is in the form of the kinetic energy of the positrons, which are assumed to be fully trapped. This model can reasonably well mimic the observed light curve from day 64 to day 538 for a value of $`t_1=105`$ days (Fig. 5). This means that the trapping of gamma-rays is almost complete at day 64, but that the positron contribution is not dominating until after $`600`$ days.

At the *HST* phases the light curve flattens out and can no longer be explained in terms of $`{}_{}{}^{56}\mathrm{Co}_{}^{}`$. An interesting possibility is then a contribution from $`{}_{}{}^{57}\mathrm{Co}_{}^{}`$, in particular since the decay rate in the *HST* detections seems to agree well with the decay time of $`{}_{}{}^{57}\mathrm{Co}_{}^{}`$$``$$`{}_{}{}^{57}\mathrm{Fe}_{}^{}`$.

However, the decay of $`{}_{}{}^{57}\mathrm{Co}_{}^{}`$ has no positron channel. A fairly large amount of $`{}_{}{}^{57}\mathrm{Ni}_{}^{}`$ (which quickly decays to $`{}_{}{}^{57}\mathrm{Co}_{}^{}`$) must therefore be synthesized in the explosion, in order to make a significant contribution in this model.

In Fig. 5 we have included a ratio of $`{}_{}{}^{57}\mathrm{Ni}_{}^{}`$/$`{}_{}{}^{56}\mathrm{Ni}_{}^{}`$ which is three times larger than the ratio observed in SN~1987A. The optical depth to gamma rays is also enhanced by a factor of 2.4 as compared to $`{}_{}{}^{56}\mathrm{Ni}_{}^{}`$, to allow for the smaller gamma-ray energies involved (Woosley et al. 1989). In addition, $`{}_{}{}^{44}\mathrm{Ti}_{}^{}`$ has also been added with a three times larger ratio. This nucleus contributes at much later phases in SN~1987A (Kozma 2000; Lundqvist et al. 2001), but can become important at relatively earlier phases in SN~1998bw as it has a significant positron channel.

Although these abundance enhancements are quite arbitrarily chosen to match the observations, we emphasize that the nucleosynthesis yields may indeed be different for very energetic and possibly asymmetric explosions. Nakamura et al. (2001b) and Maeda et al. (2001) clearly show that enhanced amounts of 44Ti may be expelled in such circumstances. We therefore note the possibility that the late light curve of SN~1998bw can be explained in terms of radioactivity, if the abundances of the radioactive isotopes are enhanced compared to the case of SN~1987A.

Although the above model is clearly too simplistic to give accurate results, it can be used to estimate the amount of $`{}_{}{}^{56}\mathrm{Ni}_{}^{}`$ needed to power the supernova at a distance of 35 Mpc. With the radioactive luminosities properly included (Kozma & Fransson 1998) we have used a nickel mass of $`0.3\mathrm{M}_{}`$ to power the light curve in Fig. 5.

This is close to the lower limit obtained by S00. This is indeed a lower limit, as the simple model forces the gamma-rays to be almost fully trapped at day 65, at the start of the rapid decline. The models presented in S00, with a complete calculation of the gamma-deposition for real explosion models, trapped 6–10% of the gamma-ray energy at 400 days, while the toy model presented here traps $`7`$% at this epoch. Although detailed calculations are needed to accurately determine the nickel mass, it is clear that the estimates will be lower than in S00 simply because some of the light attributed to the supernova in that analysis was due to the complex background.

The above estimate is, however, not very far from the 0.4 $`\mathrm{M}_{}`$ required to power the peak of the light curve in the refined models of Nakamura et al. (2001a), especially when the different distance estimates are taken into account.

We note that a scenario where the steep light curve decline is instead due to positron escape will clearly require much higher nickel masses to account for the observed luminosity. In such a scenario the late time contribution of $`{}_{}{}^{57}\mathrm{Co}_{}^{}`$ is less likely, because with no positron channel an unrealistically high abundance of this isotope would be required.

### 5.2 Other Powering Mechanisms?

It must be clarified that apart from radioactive heating, several other mechanisms could power the late light curve. None of the usual suspects can really be ruled out based on the sparse observational constraints we have at the latest phases.

Interaction with circumstellar material (CSM) is not uncommon in core-collapse supernovae (e.g., Leibundgut 1994). In the context of SN~1998bw, Tan et al. (2001) proposed an association of GRB~980425 based on a model with a high pre-supernova mass loss of $``$few $`\times 10^4`$ $`\mathrm{M}_{}`$ yr-1. Weiler et al. (2001) also argued for a CSM from radio monitoring of the supernova. However, the very fast wind expected from a WR progenitor can give a rather modest density even with a high mass loss rate, and the interaction need not give rise to optical emission. We note that as long as the supernova was spectroscopically monitored, no spectral signatures from CSM interaction were ever detected in SN~1998bw (S00). Although absence of evidence is not evidence of absence, this may constrain the most massive progenitor models.

If SN~1998bw formed a black hole which accretes matter, this could in principle show up in the late light curve (Fynbo et al. 2000). We note that the luminosity at late phases is close to the Eddington luminosity, which would require very high efficiency. Moreover, an accretion scenario should provide a power-law decay rate (Balberg et al. 2000) , which is not observed. Black hole powering is mainly expected to be seen for supernovae with very little radioactive material, because the radioactive heating will otherwise outshine the accretion luminosity. This is quite the opposite of SN~1998bw.

Another scenario for late-time emission is a light echo, as observed in SN~1998bu (Cappellaro et al. 2001). Again, we have little information available to exclude such a possibility. We can only note that the supernova appears fairly red in the day 778 data (LP versus CL), while a light echo should reflect an early blue phase.

Finally, the late light curve might also be boosted by a freeze-out effect. As mentioned above, this was observed in SN~1987A and contributed to the light curve tail powered by $`{}_{}{}^{57}\mathrm{Co}_{}^{}`$, as illustrated in Fig. 1 of Fransson & Kozma (1993). Most of the freeze-out in SN~1987A took place in the extended hydrogen envelope, which is absent in SN~1998bw. On the other hand, the faster expansion of SN 1998bw favors a freeze-out scenario. If freeze-out was indeed important at these phases, smaller yields of the long-lived nuclei would be required to power the light curve. Detailed time-dependent modeling is needed to quantify this effect.

## 6 Summary

We have presented new and updated photometry of the famous GRB supernova SN~1998bw. Using the template subtraction technique we have revised the late-time light curve and shown that the fast decline rate continues to more than 500 days past explosion. This highlights the observational problem of accurate photometry in complex regions, where core-collapse supernovae are known to reside. We note that no sign of a fully trapped $`{}_{}{}^{56}\mathrm{Co}_{}^{}`$ positron phase is observed. A simple radioactive model with $`0.3`$ $`\mathrm{M}_{}`$ of ejected $`{}_{}{}^{56}\mathrm{Ni}_{}^{}`$ can fit the data, and can be regarded as a lower limit to the amount of ejected nickel in SN~1998bw. Our very late HST observations allowed us to follow the supernova to $`1000`$ days, and to detect a late flattening of the light curve. One of the possible scenarios to explain this is the contribution of radioactive isotopes with longer life times. This could require different nucleosynthesis for the very energetic SN~1998bw than for SN~1987A.

###### Acknowledgements.

We thank Paolo Mazzali and Keiichi Maeda for interesting discussions, and all members of the SINS team for support. This research was supported by NASA through grant GO-2563.001 to the SINS group from STSCI, which is operated by AURA, Inc., under NASA contract NAS5-26555. A.V.F. acknowledges the support of the Guggenheim Foundation Fellowship. S.T.H. acknowledge support from the NASA LTSA grant NAGS–9364.