Day 640 infrared line and continuum measurements: Dust formation in SN 1987A

We have measured day 640-645 line and continuum spectra of (Ni II) 6.6 micrometer (Ne II) 12.8 micrometer (line emission was not detected), and (Fe II) 17.9 and 26.0 micrometer from SN 1987A. The high velocity feature at v(sub HVF) approximately 3900 km/sec found in both of our day 410 (Fe II) spectra is again detected in the day 640 (Ni II) spectrum, although the signal-to-noise of the day 640 (Fe II) spectra is insufficient to show this feature. The continuum fluxes provide clear evidence for the formation of dust between day 410 and day 640 and are best fitted by a graybody spectrum with a temperature of 342 +/- 17 K at day 640 and a surface area corresponding to a minimum dust velocity v(sub dust) = 1910 +/- 170 km/sec. Optically thin dust emissivity laws proportional to lambda(exp -1) or lambda(exp -2) are inconsistent with the data. Either the dust grains are large (radius a much greater than 4 micrometer and radiate like individual blackbodies, or else they are located in clumps optically thick in the 6-26 micrometer range. The (Ni II) 6.6 micrometer line flux yields a minimum Ni(+) mass of 5.8 +/- 1.6 x 10(exp -4) solar mass and a Ni/Fe abundance ratio of 0.06 +/- 0.02, equal to the solar value. The ratio of the two (Fe II) line profiles implies a gas temperature 2600 +/- 700 K, a drop of 1800 +/- 800 K from our day 410 measurement. The (Fe II) 26.0 micrometer line flux has decreased by a factor of 2 and the day 640 (Ni II) profile is blueshifted by -440 +/- 270 km/sec, relative to observations before day 500. We show that the decrease in the (Fe II) flux and the blueshift are not produced by a decrease in electron scattering optical depth, electron density, or temperature, but rather are probably due to obscuration by the same dust which produces the infrared continuum. This supports the interpretation that the dust spectrum is produced by optically thick clumps. We discuss possible explanations for the discrepancy between the mass of Fe(+) detected and the total iron mass required to power the light curve. The decrease in the (Fe II) fluxes relative to the decrease required to account for the blueshifts of optical lines from non-iron-group elements and the similarity between v(sub dust) and the Ni(+) expansion velocity imply a spatial association between the dust clumps and the iron-group elements. In addition, the larger blueshift observed for the near and far-infrared, heavy metal transitions relative to non-iron-group lines suggests that the iron-group elements are somewhat segregated from lighter elements such as the Mg(sup 0) and O(sup 0) responsible for shorter wavelength lines. We speculate that FeS may be an important constituent of the dust. A comparison of our line profiles with radiative transfer models shows that while power law and exponential density distributions yield reasonable fits to the data, polytrope distributions provided significantly worse agreement. The best fits require a substantial fraction of the iron to be undetectable, and are consistent with maximum expansion velocities of v(sub max) approximately 3000 km/sec.