{'fraction_non_alphanumeric': 0.030623141500599648, 'fraction_numerical': 0.03808178218797746, 'mean_word_length': 4.546746856205576, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Large Late-time Asphericities in Three Type IIP Supernovae

Type II-plateau supernovae (SNe IIP) are the results of the explosions of red supergiants and are the most common subclass of core-collapse supernovae. Past observations have shown that the outer layers of the ejecta of SNe IIP are largely spherical, but the degree of asphericity increases toward the core. We present evidence for high degrees of asphericity in the inner cores of three recent SNe IIP (SNe 2006my, 2006ov, and 2007aa), as revealed by late-time optical spectropolarimetry. The three objects were all selected to have very low interstellar polarization (ISP), which minimizes the uncertainties in ISP removal and allows us to use the continuum polarization as a tracer of asphericity. The three objects have intrinsic continuum polarizations in the range of 0.83-1.56% in observations taken after the end of the photometric plateau, with the polarization dropping to almost zero at the wavelengths of strong emission lines. Our observations of SN 2007aa at earlier times, taken on the photometric plateau, show contrastingly smaller continuum polarizations (~0.1%). The late-time H-alpha and [O I] line profiles of SN 2006ov provide further evidence for asphericities in the inner ejecta. Such high core polarizations in very ordinary core-collapse supernovae provide further evidence that essentially all core-collapse supernova explosions are highly aspherical, even if the outer parts of the ejecta show only small deviations from spherical symmetry.


INTRODUCTION
Core-collapse supernovae (SNe) are the explosive deaths of massive stars. Indirect inferences that these energetic explosions are aspherical have come from such evidence as the high space velocities of the pulsars produced by the SNe (Lyne & Lorimer 1994) and the structure of Galactic supernova remnants (e.g., Fesen 2001). A direct observational test for the presence of asphericity in SNe was proposed by Shapiro & Sutherland (1982), who pointed out that the polarizing effects of electron scattering in a supernova atmosphere would not cancel if the supernova were aspherical and spatially unresolved by a distant observer.
The first published detections of supernova polarization were for SN 1987A (Barrett 1988;Cropper et al. 1988;Mendez et al. 1988). Subsequent observations of additional SNe led to the discovery that all core-collapse supernovae with data of sufficient quality exhibit polarization, as recently reviewed by Wang & Wheeler (2008). An important test of our understanding was provided by the late-time Hubble Space Telescope (HST ) images of SN 1987A which spatially resolve the ejecta and show that they are elongated along a position angle (P.A.) consistent with that inferred from the polarimetry (Wang et al. 2002).
Our study here focuses on SNe IIP, which represent the single most common outcome for the ends of the lives of massive stars. Approximately 50% of all corecollapse SNe are of this subtype (Li et al., in prep.; see also Smartt et al. 2009). Their progenitors are relatively low-mass red supergiants in the mass range of 8.5−16.5 M ⊙ ) that are either isolated or sufficiently separated from their binary companions so as to leave their evolution and outer envelopes unaffected (Heger et al. 2003;Li et al. 2006;Smartt et al. 2009). In many ways, their progenitors and explosions should be the simplest ones to understand and model.
Despite this, published spectropolarimetry of SNe IIP is surprisingly limited. Most of the objects of this type listed in the compilation of Wang & Wheeler (2008) as having available polarimetry are either unpublished or have only photometric polarimetry. Our focus here is on measuring continuum polarization and separating intrinsic continuum polarization from the effects of interstellar polarization (ISP); with only photometric data, this is extremely difficult, if not impossible.
The two SNe IIP that have been the most well-studied spectropolarimetrically are SNe 1999em and 2004dj. SN 1999em showed low polarization shortly after explosion, indicating that the outer ejecta were nearly spherical. The level of continuum polarization slowly increased with time in a manner consistent with a single axis of symmetry while simultaneously the polarization signatures in the lines grew Wang & Wheeler 2008). SN 2004dj showed no detectable polarization (after correction for ISP) in the continuum or lines while on the plateau. At the end of the plateau stage, the continuum polarization jumped up to 0.56%, revealing an aspherical core (Leonard et al. 2006). The polarization percentage slowly decreased thereafter in a manner attributed to geometric dilution of the scattering electrons as the ejecta expanded. The observations of SN 2004dj prompted the present study in an attempt to determine whether the results for SN 2004dj were common to SNe IIP.
Three recent SNe IIP in the Virgo Cluster have provided an opportunity to study this class of objects in more detail. SN 2006my was discovered by K. Itagaki in NGC 4651 on 2006 November 8.82 (UT dates are used throughout this paper; Nakano & Itagaki 2006a). The same observer also discovered SN 2006ov in M61 on 2006 November 24.86 (Nakano & Itagaki 2006b). Both objects were soon spectroscopically confirmed as SNe IIP a few months past explosion (Stanishev & Nielsen 2006;Blondin et al. 2006). SN 2007aa in NGC 4030 was discovered by T. Doi on 2007 February 18.308 and spectra obtained by N. Morrell on Feb. 19.24 showed that it was a SN IIP similar to SN 1999em at about 20 days after explosion (Doi 2007;Folatelli & Morrell 2007). Li et al. (2007) presented photometric and spectroscopic evidence that both SN 2006my and SN 2006ov were normal Type IIP SNe discovered near the end of the photometric plateau stage.
While Li et al. (2007) identified possible progenitors for SNe 2006my and 2006ov in pre-explosion HST images, subsequent work has rejected those identifications (Leonard et al. 2009;Smartt et al. 2009;Crockett et al. 2009). Instead, Smartt et al. (2009) and Crockett et al. (2009) give upper limits for the zero-age main-sequence masses of any progenitor in the HST images of 13, 10, and 12 M ⊙ for SNe 2006my, 2006ov, and 2007aa, respectively. Smartt et al. (2009) also quote values for the 56 Ni masses produced by SN 2006ov and SN 2006my of 0.003 ± 0.002 M ⊙ and 0.03 ± 0.015 M ⊙ , respectively, based on analyses of the late-time light curves. The 56 Ni mass inferred for SN 2006ov is extremely low, comparable to that of the underluminous SN 1999br (Hamuy 2003;Pastorello et al. 2004), despite SN 2006ov having a normal luminosity on the plateau.
In this paper, we present spectropolarimetry of the above three SNe IIP obtained after the end of the photometric plateau. In §2, we describe the observations and corrections for small amounts of ISP from the Galaxy. We present the evidence for very high intrinsic continuum polarizations in these objects in §3 and discuss the interpretation and implications of our results in §4.

OBSERVATIONS
Unfiltered images of all three SNe in our sample were obtained as part of our regular monitoring of their host galaxies for the Lick Observatory Supernova Search using the 0.76-m Katzman Automatic Imaging Telescope (KAIT; Filippenko et al. 2001), although we artificially increased the frequency of observations from the default search strategy to better sample the light curves of these objects. Portions of the SN 2006ov and SN 2006my light curves have already been presented and the reductions described by Li et al. (2007). Briefly, a point-spread function fitting technique was used to subtract pre-SN template images from the SN observations and aperture photometry was performed on the subtracted images using DAOPHOT (Stetson 1987) in IRAF 3 . KAIT unfiltered observations approximate the standard R band (Li et al. 2003), so the SN observations were then calibrated to an average of USNO-B1 stars in the field (up to 20 per field), resulting in systematic uncertainties in the zero points for each field of 0.05 − 0.18 mag. However, only the shapes of the light curves are relevant for this study.
We obtained spectopolarimetry of the three SNe using both the Kast spectrograph on the Lick 3-m Shane telescope (Miller & Stone 1993) and the Low Resolution Imaging Spectrometer (LRIS) mounted on the Keck I 10-m telescope (Oke et al. 1995). The Kast data were obtained using the red camera only, with the 300/7500 grating tilted to cover a wavelength range of approximately 4400−9900Å in conjunction with a GG455 orderblocking filter. A small amount of second-order light is likely present at wavelengths greater than ∼ 9000Å, but none of our conclusions depend on data at those wavelengths. A 3 ′′ -wide slit gave a spectral resolution of ∼ 16Å.
The LRIS data were taken with both halves of the spectrograph. All observations used an identical setup with the D560 dichroic beam-splitter to split the spectrum near 5500 − 5600Å, along with the 600/5000 grism on the blue side and the 400/7500 grating on the red side. We used a 1. ′′ 5-wide slit to achieve spectral resolutions of ∼ 6Å and 9Å on the blue and red sides, respectively. The total spectral coverage was 3170 − 9240Å, but the polarization data plotted in this paper are typically truncated below ∼ 4000Å because the low blue flux in SNe IIP due to line blanketing, and decreasing detector sensitivity at short wavelengths, combine to make the polarimetry extremely noisy at the blue end. In addition, this spectral region suffers from background-subtraction uncertainties due to the underlying blue starlight from the host galaxy. Further details of the individual observations are given in Table 1.
Basic two-dimensional image processing and spectral extraction were accomplished using standard tasks in IRAF. Flux calibration and removal of atmospheric absorption bands (Wade & Horne 1988) was performed using our own IDL tasks (Matheson et al. 2000). The individual one-dimensional SN spectra were combined on a pixel-by-pixel basis, weighted using the error estimates from the optimal spectral extraction technique (Horne 1986), and the final blue-side and red-side spectra (for LRIS) were scaled and combined across the overlap region. In those cases where differential loss of blue light was apparent because the observed P.A. of the observations deviated from the parallactic angle (Filippenko 1982), the combined observation was warped to the shape of the observation obtained closest to the parallactic angle using a low-order polynomial fit to the ratio of the spectra.
We also extracted some of the emission from nearby H II regions along the slits of selected observations. We fit Gaussians to the nebular emission lines and obtained recession velocities of 920, 1660, and 1650 km s −1 for the H II regions near SN 2006my, SN 2006ov, and SN 2007aa, respectively. These values are within a typical value for a galactic rotation velocity (100 − 200 km s −1 ) of those listed for the host-galaxy nuclei in the NASA Extragalactic Database (NED), so we adopt these redshifts in this paper and have removed them from all plots. All references to wavelengths are in the rest frame of each object.

Spectropolarimetry
The spectropolarimetric observations were reduced following the standard method of Miller et al. (1988), as implemented by . The angle zero points of the waveplates were determined using polarimetric standards from Schmidt et al. (1992). Typically, two standards were observed per night to confirm the zero point. In addition, we noted that the zero points for each instrument were generally stable over time with small shifts. Low-polarization standards from the literature (Mathewson & Ford 1970;Clemens & Tapia 1990;Schmidt et al. 1992;Berdyugin et al. 1995) were observed on each night to have negligible polarizations, setting a limit on the instrumental polarization of P < 0.1%. This study hinges on our ability to accurately measure the continuum polarization of SNe. ISP in both our own Galaxy and the host galaxies of these SNe can mimic the presence of continuum polarization and contaminate our results. All three SNe selected for this study are at high Galactic latitudes (b > 59 • ) and hence low Galactic ISP is expected (Serkowski et al. 1975) based on the low columns (E(B − V ) < 0.03 mag) measured along these lines of sight (Schlegel et al. 1998). We can directly measure the Galactic contribution to the ISP by following the lead of Tran (1995a) and selecting stars of low intrinsic polarizations and sufficiently large distances (estimated via spectroscopic parallaxes) to place them out of the Galactic plane. Such probe stars should sample most of the Galactic contribution to the ISP along these lines of sight. The probe stars for SN 2006my were BD +17 • 2526 and BD +16 • 2409, those for SN 2006ov were BD +05 • 2618 and BD +04 • 2608, while BD +00 • 2876 and BD −00 • 2514 were selected for SN 2007aa.
The probe stars were observed immediately after the object on one observation date per SN with the slit kept at the same P.A. Each pair of stars gave consistent results (within 0.1% in each Stokes parameter), so the two stars were averaged to form our estimate of the Galactic ISP for each SN. As expected, the Galactic ISPs were very small, with V -band polarizations of 0.07%, 0.06%, and 0.16% measured for the SN 2006my, SN 2006ov, and SN 2007aa probes, respectively. For the sake of completeness, we subtracted low-order polynomial fits to the Stokes q and u of the combined probes from the Stokes parameters measured during each SN observation, although the effect is almost negligible for SN 2006my and SN 2006ov.
As a final check on our instrumental stability, we observed the same probe star, BD +05 • 2618, on all three occasions that we obtained spectropolarimetry of SN 2006ov. The variation in each Stokes parameter was < 0.04%, consistent with similar experiments we have performed in the past with repeated observations of bright stars on different nights. This sets an upper limit to any systematic offset between the Lick and Keck polarimeters as well as to any night-to-night instrumental instability. We note that Leonard et al. (2006) also found similar results using the Kast polarimeter. Therefore, we adopt a value of 0.04% for any systematic error contribution where appropriate in the analysis below. Li et al. (2007) for SN 2006my and SN 2006ov, and below for SN 2007aa, all three objects in our study exhibit standard light-curve shapes for SNe IIP, with the three phases of a relatively flat plateau followed by a rapid decline phase and then a steady decline at a rate consistent with the radioactive decay of 56 Co. This characteristic shape is well understood theoretically (Eastman et al. 1994;Kasen & Woosley 2009) and defines this subclass of SNe (Barbon et al. 1979;Doggett & Branch 1985). The supernova shock deposits thermal energy in the envelope which is released as the hydrogen in the envelope recombines during the plateau phase. When the photosphere recedes all the way back through the hydrogen layer, recombination ceases to be a source of energy and the SN light curve rapidly falls until the luminosity matches the instantaneous rate of energy input from the decay of 56 Co.

As shown by
Both SN 2006my and SN 2006ov were discovered several months after explosion near the ends of their plateaus.
In order to make temporal comparisons amongst our observations of these SNe when the dates of explosion are so poorly known, we found it convenient to define a common time to represent the end of the plateau. As discussed below, the timing relative to the end of the plateau is probably the important physical parameter for the polarization, not the time elapsed since explosion. Hamuy (2003) defined the end of the plateau phase, t p , as the time when the SN magnitude falls halfway between its value on the plateau and the value at the start of the radioactive decay tail. The plateau magnitude can be difficult to uniquely define in some cases due to the varied shapes of SN IIP light curves. The compilation of SN IIP photometry by Poznanski et al. (2009) shows that some objects have "plateau" phases that actually rise for many days before falling while others steadily, if slowly, decline from the beginning. The photometric band used to define the light curve can also have an effect as the shapes of SN IIP light curves are a function of wavelength. Hamuy (2003) used the V band, but we have only unfiltered photometry for our three objects.
Despite these uncertainties in the definition of the plateau magnitude, t p occurs during the rapid decline phase and hence the uncertainty in t p is at worst a few days, which is negligible for our purposes. We measure dates for t p of 2006 Dec. 15.0, 2006Dec. 25.5, and 2007tively, based on the KAIT unfiltered photometry. Below, we measure t p for SN 1999em in the R band (which most closely approximates our KAIT unfiltered photometry), and our value for t p differs by only 0.2 days from the value derived by Hamuy (2003) in V , so we can conclude that bandpass effects have negligible effect on our measured values for t p . Hereafter, dates listed as "day X" refer to epochs of observation X days after the end of the plateau as defined by t p . Li et al. (2007) found that SN 2006my and SN 2006ov were not highly reddened and had no detectable Na I D absorption lines from the interstellar media of their hosts. We confirm the lack of narrow Na I absorption in our spectra. Our highest signal-to-noise ratio (S/N) spectra of each object are plotted in Figure 1. SN 2007aa does have a weak dip near 5888Å that may be from absorption in NGC 4030. The strength of the feature in SN 2007aa is difficult to quantify at our low spectral resolution, especially with the curvature of the spectrum making the continuum hard to define. We estimate the equivalent width of the dip to be ∼ 0.08Å, with an uncertainty of about 50%. If two-thirds of the absorption is due to the D1 component of the doublet (reflecting the ratio of the oscillator strengths, as expected in the low optical depth limit), the relationship of Munari & Zwitter (1997) predicts a host-galaxy reddening of E(B − V ) ≈ 0.014 mag, which corresponds to an ISP of P ≈ 0.04% for typical lines of sight in the Galaxy, up to an empirical maximum of P ≈ 0.13% (Serkowski et al. 1975). SNe 2006my and 2006ov appear to have even less absorption from their host galaxies. Even if our estimates of the reddening are off by a factor of a few, the conclusion that we should expect our objects to have low host-galaxy ISP will be unchanged.

SN 2006ov
We obtained three epochs of spectropolarimetry on SN 2006ov, two during the steep fall off the photometric plateau and one a month later. Our highest S/N data were obtained from Keck on day 0 and were fortuitously exactly synchronized with the end of the plateau. These data are plotted in the q − u plane in Figure 2, with each black circle representing a 100Å bin of data. As described above, the S/N degrades to the blue, so we only plot data at wavelengths greater than 4500Å. The cloud of points is clearly elongated in the positive q, positive u direction by an amount in excess of the error bars, indicating substantial intrinsic supernova polarization.
We selected two broad wavelength intervals, 6800-7200 and 7820-8140Å, as representative of the line-free continuum (the exact choice is justified below) and plotted the integrated polarizations 4 as the blue diamond and SN 2007aa shows a tiny dip near 5888Å that may be Na I absorption in its host galaxy. The other two SNe show even less evidence for host-galaxy absorption, indicating that we should expect negligible ISP from the host galaxies of these objects.
solid red circle in Figure 2. The two continuum points are very consistent with each other and have an average polarization (weighted by the inverses of their variances) of P = 1.56 ± 0.03% at P.A.= 8. • 2 ± 0. • 6. Many of the 100Å bins of data are clustered around these continuum points and along the line connecting them to the origin. The P.A. of 8. • 2 therefore is significant in this object and probably represents the projection of the SN symmetry axis onto the plane of the sky. Since the Stokes parameters (q,u) are referenced to the arbitrary direction of North on the sky, it is useful to rotate our coordinate system to align one Stokes parameter with the symmetry axis of the SN. We call these new coordinates rotated Stokes parameters (RSP; Trammell et al. 1993a;Tran 1995b). If the SN were axisymmetric, q RSP would be an estimator of the total polarization and u RSP would show any deviations from axisymmetry. The polarization data after rotation are plotted in Figure 3 in 20Å bins. The most striking aspect of Figure 3 is the high level of polarization (P ≈ 1.5%) in the continuum, with many strong, sharp depolarizations down to almost zero in q RSP present at the lines. We integrated the polarimetry over the wavelength range 6540-6600Å to determine the polarization at the peak of Hα and measured (q RSP , u RSP ) of (0.08%, 0.00%), with error bars of 0.05% in each Stokes parameter. The similar depolarizations present for many other lines are convincing evidence that the host-galaxy contribution to the ISP must be very small ( 0.1%), as might be expected from the low reddening to SN 2006ov. To estimate the ISP, many past workers have made the assumption that the intrinsic polarization in supernovae at the peaks of strong lines is zero (Trammell et al. 1993b;Höflich et al. 1996;Tran et al. 1997), but that assumption is clearly wrong in some cases (Maund et al. 2007b), especially at early times when He I λ6678 is blended with Hα. Here, the continuum polarization we measure must be almost entirely intrinsic to SN 2006ov and the lines are in fact

depolarized.
With that in mind, we can now justify our choice of continuum intervals plotted in Figure 2 and marked on Figure 3. Past work has noted the utility of measuring the continuum polarization in the red portion of the optical where few line features interfere with the measurement. Leonard et al. (2006) used a wide continuum region of 6800-8200Å in their study of SN 2004dj. Integrating the observed polarization over a continuum region that wide for SN 2006ov will clearly underestimate the true polarization as it will include a strong depolarization present at the K I/O I blend as well as a weaker one due to [Ca II] emission. We formally get P = 1.26±0.02% if we use the wide definition for the continuum, significantly less than the value of 1.56% quoted above. Therefore, we define two continuum regions to exclude the strong lines. The first one, 6800-7200Å, is designed to avoid the red wing of Hα at its blue end and also exclude the [Ca II] emission at the red end. The second region, 7820 − 8140Å, is bracketed by O I λ7774 and Na I λλ8183, 8194.
We also obtained SN 2006ov data on two additional epochs, days −5 and 27. We have plotted a comparison of q RSP on all three epochs in Figure 4. Due to the lower S/N of the additional data, we have rebinned to 50Å per pixel. Despite the significant noise, the overall pattern is the same: large continuum polarization with dips down to nearly zero present at the wavelengths of strong line features. The Lick data were taken only 5 days before the day 0 Keck data and do show some significant differences, indicating rapid evolution of the spectrum as the SN made the transition off the photometric plateau. For example, the depolarizations in the wavelength range 6000 − 6500Å are more blended together at day −5 than at day 0.
We calculated the polarization in both of our contin- -Spectropolarimetry of SN 2006ov on day 0, presented as 20Å bins in the rotated q RSP -u RSP coordinate system so that q RSP is an estimator of the polarization and u RSP is near zero. The total-flux spectrum of SN 2006ov is overpolotted in the top panel to guide the eye, with several major lines labeled. The continuum polarization between strong line features is very high (∼ 1.5%), but the lines result in depolarization. The two horizontal black bars mark the extent of our two relatively line-free regions (6800-7200 and 7820-8140Å), which we take to be representative of the continuum. The data points with error bars represent the Stokes parameter q RSP , which is an estimator of the total polarization. As a reference, the solid line in each figure is the total-flux spectrum from each epoch. Note the high continuum polarization, with depressions at the location of the strong lines. The high point near 8150Å in the day −5 polarization is spurious and due to relatively poor night-sky subtraction. The two horizontal black bars mark our two continuum regions.
uum regions for each of our three epochs. The results are plotted in Figure 5 and compared to the KAIT unfiltered light curve. The two continuum windows gave very consistent results, showing a significant increase in the weighted-average polarization in only five days, from 1.25±0.06% on day −5 to 1.56±0.03% on day 0, falling to 1.15±0.08% on day 27. SN 2004dj showed large variations in the polarization angle at the end of the plateau (Leonard et al. 2006), but the SN 2006ov continuum observations are consistent with a single P.A. to within 1 • . Unfortunately, SN 2006ov exploded while in solar conjunction and we were unable to obtain observations in the month between discovery and our day −5 data, so it was impossible to establish whether the high polarization apparent at late times represented a significant change from early times.

SN 2006my
We obtained a single epoch of spectropolarimetry of SN 2006my on day 38, 2007 Jan. 21. We analyzed the data in a similar manner to SN 2006ov above. First, we measured the weighted average polarization in our two continuum windows to be P = 0.97 ± 0.04% at P.A. = 174. • 1 ± 1. • 2. Then we rotated our Stokes parameters by 174. • 1 to align q RSP with this polarization angle. The results are plotted in Figure 6.
The data are of lower quality than for SN 2006ov but show broadly similar characteristics. High continuum polarization is seen in q RSP , with deep depolarizations present at the wavelengths of strong line features. Li et al. (2007) found that the extinction to SN 2006my was likely to be minimal, and hence low ISP should be expected. We determined the polarization at the peak of Hα using the same window as for SN 2006ov above and found (q RSP ,u RSP ) of (0.03%, 0.03%) with an un- certainty of 0.04% in each Stokes parameter. Just as for SN 2006ov, the data for SN 2006my are consistent with negligible ISP in the host galaxy. In this object, the continuum polarization is distinctly lower than in SN 2006ov, although it is unclear how much of the difference is due to the later epoch of the SN 2006my observations.

SN 2007aa
Unlike the other two objects, SN 2007aa was discovered shortly after explosion, allowing us the opportunity to obtain data at early times and have a comparison for the late-time data.
Our first epoch of spectropolarimetry was obtained on day −50, well before the end of the plateau, and only 28 days after discovery. As discussed above, Folatelli & Morrell (2007) described a spectrum taken one day after discovery as being similar to that of SN 1999em at 20 days after explosion, implying an epoch for our first observation of ∼ 47 days after explosion. We used the SuperNova IDentification code of Blondin & Tonry (2007) to cross-correlate our day −50 spectrum with a library of supernova templates, holding the redshift fixed. The top three individual SNe with spectral matches (ignoring spectra on multiple dates per object) were SNe 2004et, 2004dj, and 2006bp at 56, 56, and 58 days (respectively) after their estimated explosion dates (Sahu et al. 2006;Vinkó et al. 2006;Quimby et al. 2007).
These estimates are all consistent with SN 2007aa having a normal plateau duration of ∼ 100 days, with our first set of data being taken near the middle. Our second epoch, at day −22, was taken late in the plateau stage, but before the beginning of the steep decline phase, while our third epoch, on day 3, was taken during the steep decline and just before the beginning of the radioactivedecay tail.
The two plateau epochs of spectropolarimetry are plotted versus wavelength in Figure 7 and in the q − u plane in Figure 8. We have applied a rotation of 172 • to the data in Figure 7  To examine this in more detail, we have plotted the data points in the Hα line from day −22 in Figure 8 as the solid black circles. The data points are 50Å bins in wavelength from 6400 to 6750Å, representing a velocity range of −7500 to 8400 km s −1 relative to Hα. The radius of each circle is proportional to wavelength, so the smaller points are in the blue wing of the profile and the larger circles are from the red wing. The Hα points show a spread in both q and u in a manner inconsistent with a simple axisymmetric geometry. In axisymmetry, all the points would lie along a line in q − u space. Instead, the points show a q − u "loop." As one moves across the Hα profile from red to blue, one moves in a counterclockwise direction around the loop.
Such loops were first identified in supernova spectropolarimetry by Cropper et al. (1988) in the Hα and Ca II near-infrared (NIR) triplet lines of SN 1987A. The interpretation of the loops in SN 1987A has remained obscure for over 20 years despite the high quality of the data (Wang & Wheeler 2008). Subsequent investigations have found q − u loops in systems as different from SN 1987A as SNe IIn (Hoffman et al. 2008), stripped-envelope SNe (Maund et al. 2007a,b), and even SNe Ia (Kasen et al. 2003;Wang et al. 2003a;. Kasen et al. (2003) found that loops were a generic consequence of systems containing multiple components with misaligned axes of symmetry and constructed a few examples with very different geometries. The relevance of these models to SN II atmospheres during the plateau stage is unclear as both the photosphere and Hα-forming region are within the hydrogen envelope of the progenitor. Chugai (1992) has a more promising model, specifically designed for SN 1987A, which interprets the Hα loops as evidence of a small number of individual clumps data at wavelengths greater than 8200Å is unexplained and may not be real, although none of the standard or probe stars observed that night showed such an effect. In any case, none of the results presented in this paper rely on data from that portion of the spectrum. The horizontal black bar in each q RSP panel marks our continuum region.
of 56 Ni creating pockets of high ionization and Hα excitation. At different radial velocities, we see the effects of different clumps causing the net polarization angle to rotate with wavelength.
We defer additional analysis of the line polarization features in SN 2007aa to future work on the polarization of SNe II at early times because our focus here is on the continuum polarization and its temporal evolution. For now, we will note the following points. First, the long axis of the Hα loop on day −22 is mostly aligned with the P.A. of the late-time continuum polarization, possibly indicating a common symmetry axis. Second, the total polarization excursion in the line is ∼ 1% at an epoch where the continuum polarization is likely to be very much lower, as discussed below.
We will now discuss the continuum polarization in detail. We found our low-S/N late-time polarization data in the final epoch to be unreliable near 8100Å, so we have focused only on the bluer of our two continuum polarization windows, 6800-7200Å. Figure 8 has our continuum polarization measurements plotted for each of our three epochs of observation. As was clear from Figure 7, the observed continuum polarization on the plateau was very low. We measure formal polarizations of 0.11 ± 0.04% and 0.08 ± 0.05% on days −50 and −22, respectively. By contrast, the polarization after the end of the plateau jumped to 0.83 ± 0.12% on day 3.
The late-time spectropolarimetry data of SN 2007aa from day 3 are plotted in Figure 9. The data are extremely noisy due to the faintness of the object at this time, but when binned up to 100Å per pixel the basic appearance of the data is qualitatively similar to that shown previously for SN 2006ov and SN 2006my. The continuum polarization is near 1% in q RSP , but there are broad dips down to near zero polarization at the locations of strong lines such as Hα, with u RSP centered near zero. The polarization we measure at the peak of Hα is (q RSP , u RSP ) = (0.02%, 0.02%), with 0.1% error bars in each Stokes parameter, consistent with zero, as shown in Figure 8.
A few of the dips in polarization at the lines appear to reach negative values of q RSP , especially near 7650Å and 8400Å, but these are only marginally statistically significant. The wavelength bins that are the most negative are those with the least flux from the SN (due to absorptions from the telluric A band and the Ca II NIR triplet P-Cygni feature) and hence the lowest S/N, magnifying any systematic uncertainties in the sky subtraction. An integration over 8500−8800Å, representing the emission component of the Ca II NIR triplet, gives (q RSP , u RSP ) = (−0.09 ± 0.17%, −0.11 ± 0.17%), which like the polarization at the peak of Hα is completely consistent with zero.
If we interpret the late-time polarization of SN 2007aa as being similar to that of SN 2006ov and SN 2006my, then the contribution from ISP should be low in this object as well, given the low polarization in the lines. Further evidence for this can be seen in Figure 8, where we have plotted the peaks of the Hα lines for all three epochs (integrated over the same 6540−6600Å interval). The Hα peaks in the first and last epoch are near zero. The Hα polarization in the second epoch is higher, but still low (0.3%). That is also the epoch showing clear evidence of strong polarization throughout the line, so perhaps a non-zero polarization at the peak is to be expected, particularly if an aspherical excitation distribution is present (Chugai 1992). Small amounts of ISP have a similar effect to changing the origin in the q − u plane, so there is still a possibility of a pathological ISP canceling out some intrinsic supernova polarization. However, taken together with the low continuum polarization measured on the plateau, the simplest interpretation is that SN 2007aa, like SN 2006ov and SN 2006my, suffers minimal ISP ( 0.1%).
We have taken the continuum polarization measurements and plotted them versus time in Figure 10, along with the KAIT unfiltered light curve. The polarization shows a strong increase from the plateau stage to the late-time observations taken during the rapid decline phase. Even if we are mistaken about the ISP being negligible, the observation that there was a sharp change in the polarization associated with the end of the plateau is a robust result. This is strongly reminiscent of the jump in polarization seen in SN 2004dj (Leonard et al. 2006). In SN 2007aa, however, there was evidence of asphericity in the form of line polarization on the plateau, while SN 2004dj did not even show any line polarization features before the sudden increase in continuum polarization.  and are 50Å bins in wavelength from 6400 to 6750Å, representing a velocity range of −7500 to 8400 km s −1 relative to Hα. The points form a loop in the q − u plane. The radius of each circle is proportional to wavelength, so moving across the line profile from blue to red corresponds to moving around the loop in a counterclockwise direction. Also plotted are continuum polarization points (integrations over 6800-7200Å) and a set of points representing the peak of the Hα emission line (integrations over 6540-6600Å). Note the low polarizations in the Hα peaks and in the continuum for the first two epochs, which together imply minimal ISP in this object. The day 3 continuum polarization represents a significant jump in polarization. The dashed lines represent the rotated q RSP −u RSP coordinate system, designed to have the q RSP axis go through the day 3 continuum point.

DISCUSSION
The spectropolarimetry data presented above for the three objects after the end of the plateau are all qualitatively similar with each other and with SN 2004dj (Leonard et al. 2006). The continuum polarization for each object is high, but flat and constant with wavelength between the strong line features in the red. This is a signature of the wavelength-independent nature of Thomson scattering. The lines produce large depolarizations and many overlapping lines significantly decrease the net polarization in the blue. The minimal polarizations at the peaks of Hα and lack of significant interstellar Na I absorptions are consistent with negligible ISP in all three objects, and therefore the large continuum polarizations are intrinsic to the SNe.
Classically, the expected effect of lines forming in front of an aspherical electron-scattering photosphere is to produce "inverse P-Cygni" line polarization features (Jeffery 1989). Line emission near zero velocity produces a dip in polarization near the line center, and selective blocking of forward-scattered light produces a polarization peak in the blueshifted absorption trough. Higher polarizations can also be found in the red wing of the line as the opaque photosphere obstructs the observer from seeing unpolarized line emission from the far side of the ejecta. This simple picture provides an explanation for the line polarization features in photospheric-phase spectropolarimetry of SN 1987A and SN 1999em (Jeffery 1991a, and possibly our early data for SN 2007aa as well. However, we do not see emission peaks blueward of Hα or the Ca II NIR triplet in our highest S/N late-time data (Fig. 3). Instead, the polarization is simply low relative to the continuum throughout the line profile. There are a few contributing effects.
The first is simple flux dilution. The continuum photosphere and associated polarizing electrons are located relatively deep in the supernova ejecta and intrinsically unpolarized line emission from the optically thin outer layers of the ejecta travels unimpeded to the observer and dilutes the observed polarization from the continuum value. As an example, the [Ca II] λ7300 doublet in the day 0 SN 2006ov data represents emission from the low-density parts of the outer ejecta. The line emission strengthens over time as the ejecta thin out and the supernova begins the transition to the nebular phase. This line has a flux peak a little less than a factor of ∼2 above the continuum (the continuum flux level is hard to define at this epoch because of many overlapping P-Cygni features), and the polarization dips from 1.3 − 1.4% in the nearby continuum pixels down to ∼ 0.8% in the line.
The second reason is that photons scattering in lines are generally depolarized (Höflich et al. 1996). While the outer layers of the ejecta are optically thin in the continuum, it seems likely that the optical depth in the strong lines is still significant. The depolarizations over the full line profile from lines that are strong in the cool outer ejecta, such as Na I, Ca II, and Hα, are likely manifestations of this effect. Absorption and reemission in the line causes the photons to lose the directional information previously imprinted upon them by electron scattering. The net effect of a large number of overlapping lines in the blue is to result in a decrease in polarization at those wavelengths ). In addition, the line-formation region is very geometrically extended compared to the electron-scattering photosphere, so the selective blocking of forward-scattered light is not an important effect.
A few of the strong depolarizations do not correspond to obvious strong features in the total-flux spectrum. An example is the depolarization visible near 5700Å in Figure 3, which was also present in SN 2004dj (Leonard et al. 2006). It is unclear which weak spectral feature is responsible for that depolarization and why the many other weak spectral features do not also result in similarly strong depolarizations.
The time evolution of the polarization is of significant interest. SNe IIP generally do not exhibit strong polarization shortly after explosion (Wang & Wheeler 2008), and yet SN 2004dj had a strong jump in continuum polarization near the end of its plateau (Leonard et al. 2006). The main motivation for this study was to determine whether SN 2004dj was unique, or if other SNe IIP exhibit similar polarization behavior. We have shown timeseries plots for the two objects in our study with multiepoch polarimetry in Figures 5 and 10, but SN 2007aa was the only object with data taken on the plateau to establish the low early-time polarization.
To make further progress, we searched for SN IIP polarimetry in the literature. As discussed in §1, there are surprisingly few published datasets of SN IIP spectropolarimetry in the refereed literature. Since we are interested in the intrinsic continuum polarization, we must also exclude objects with large ISP corrections because estimates of ISP are uncertain and require making assumptions. SN 1999em Wang & Wheeler 2008) and SN 2004dj (Leonard et al. 2006) are the only objects from the published literature to meet our criteria. Of the three objects in the study by , SNe 1997ds and 1999gi both had large and uncertain ISP contamination, and SN 1998A was a SN 1987A analog of unclear relevance to our SNe IIP. Leonard & Filippenko (2005) also plotted continuum polarization versus time for several SNe IIP, but three of the five objects in their sample had ISP corrections of 1% or more.
We collected photometry for the two objects from the literature and determined the date of the end of the plateau for each. The combined set of polarization versus time since t p is shown in Figure 11. We must caution the reader that the data shown in this Figure are necessarily very heterogeneous.
We also reanalyzed the published polarization data in order to make our comparisons as uniform as possible. The published SN 1999em data generally do not extend sufficiently far to the red to include our favored continuum regions. We calculated the continuum polarization using as much of our 6800-7200Å window as was available in each observation from the dataset of . The 1999 Nov. 5 and 2000 Apr. 5/9 data included the full wavelength range, while for the 1999 Dec. 8 and 17 data we integrated from 6800Å to the red wavelength end of the data. The 1999 Dec. 17 continuum polarization is the only one of the four values to differ from the synthesized V -band polarizations quoted by  by more than the error bars, possibly due to the limited wavelength range (the data from that night only extend redward to 6850Å). The data have not been corrected for ISP because the effect is believed to be small (∼ 0.1%; ).
As discussed above, applying the very wide continuum region used by Leonard et al. (2006) in their study of SN 2004dj to our SN 2006ov data would reduce the overall polarization measured on day 0 from 1.56% to 1.26%. To correct for this effect, we obtained the SN 2004dj data from D. Leonard and calculated the polarization in our two continuum windows in an identical manner to the other objects in this paper. Therefore, the peak continuum polarization plotted in Figure 11 is now 0.60%, ∼10% higher than the published value of 0.56% given by Leonard et al. (2006).
Despite these caveats, we believe the basic trend shown in Figure 11 is sound. At early times, SNe IIP are seen to exhibit low polarization (Wang & Wheeler 2008). The large, extended, slowly rotating hydrogen envelopes of their progenitors are unlikely to be substantially aspherical prior to explosion. Even an aspherical initial explosion will become more spherical when propagating outward through the remainder of the star (Chevalier & Soker 1989;Kifonidis et al. 2006). At the end of the plateau, the thick hydrogen envelope becomes mostly optically thin, allowing us to see deep into the ejecta and observe the aspherical core of the explosion. In addition, the observed polarization for a given geometry is expected to be maximized when the optical depth to electron scattering is nearly unity (e.g., Höflich 1991). To be fair, the dramatic jump in polarization was only directly seen in SN 2004dj and SN 2007aa. However, SN 2006ov does show a consistent increase in polarization of 0.3% in both of our continuum regions over the course of just five days, indicating very rapid changes in the continuum polarization.
It is instructive to compare to the supernova with the most densely sampled polarimetric data, SN 1987A. Jeffery (1991b) presented ISP-corrected broad-band polarimetry from a number of sources (Barrett 1988;Mendez et al. 1988). The BV RI polarizations do show a common trend of rising polarization at late times, with a possible spike in the B and V polarization percentages near the beginning of the radioactive-decay tail. While the SNe IIP in Figure 11 show declining polarization after the start of the radioactive tail, the SN 1987A data continue rising up to ∼ 1% at 200-300 days after explosion. The Type IIb SN 2001ig (Maund et al. 2007b) also showed a distinct jump in polarization associated with the hydrogen layer becoming optically thin and the photosphere receding into the He core.
Prior to the publication of the SN 2004dj polarization results, Chugai et al. (2005) noticed unusual line profiles present in that object in the nebular phase. The Hα, Hβ, [O I], and [Ca II] lines all exhibited a double-peaked structure with the main peak initially blueshifted by about 1600 km s −1 and moving redward with time, along with a similar peak on the red side of the profile. Chugai et al. (2005) modeled these line profiles in terms of a bipolar 56 Ni component protruding from a more spherical "cocoon" that produced an ionization and excitation asymmetry. Chugai (2006) found that minimal adjustments to the initial model could roughly match the polarization results at the second nebular epoch, although his models did not predict the decline in polarization after the beginning of the radioactive-decay tail seen in the data or the different polarization angle seen in the first nebular dataset relative to the later data.
While asphericities in the 56 Ni distribution are likely related to the polarization we observe in the three objects in our sample, none of them exhibits late-time line profiles similar to those of SN 2004dj. The latetime Hα profiles of our objects are shown in Figure 12. SNe 2006my and 2007aa have mostly rounded Hα profiles, with SN 2006my showing a small notch or dip near 6575Å. The early-time spectra of SN 2006ov do show a pair of double notches, near 6535Å and 6590Å in the day −5 spectrum, corresponding to velocities of −1250 and 1200 km s −1 with respect to Hα. The notches disappeared in just five days and the line profiles after the end of the plateau were normal, by contrast with SN 2004dj, whose Hα line-profile asymmetries remained prominent for at least 330 days after explosion (Chugai et al. 2005). The blue and red notches in the Hα profile of SN 2006ov before the beginning of the radioactive-decay tail are reminiscent of the so-called "Bochum event" in SN 1987A (Hanuschik & Dachs 1987), although at significantly lower velocities (∼ 1200 km s −1 versus ∼ 4000 km s −1 in SN 1987A). Similar notches have been seen in other SNe IIP near the end of the plateau (e.g., at v ≈ 2500 km s −1 in SN 1999em; Leonard et al. 2002a) and do likely indicate asphericities or clumps in the 56 Ni distribution (Chugai 1992;Utrobin, Chugai, & Andronova 1995).
Given the evidence for a potentially clumpy 56 Ni distribution, it is worth considering whether clumps of 56 Ni or other deformations of the core-envelope boundary could be responsible for the continuum polarization that we measure, rather than an aspherical photosphere in the core. Leonard et al. (2006) hypothesized that the change in the polarization angle seen from the first to the second epoch after the end of the plateau in SN 2004dj was due to some non-axisymmetric portion of the 56 Ni distribution beginning to be uncovered before the whole core as the photosphere receded in the ejecta. Chugai (2006) proposed an alternative explanation, which he termed a "spotty photosphere," invoking an inhomogeneous distribution of the elements at the hydrogen-helium boundary due to fluid instabilities at that interface during the explosion (Müller et al. 1991). These inhomogeneities would result in brightness and opacity variations that would be manifested as a nonzero net polarization that suddenly appears when the overlying layers become optically thin and the photosphere passes through the coreenvelope boundary, even if the core of the SN were essentially spherical on large scales.
Our data for these three objects  5 SN 1999em SN 2004dj SN 2006my SN 2006ov SN 2007aa Fig. 11.-Compilation of SN IIP intrinsic continuum polarization measurements after correction for ISP. The SN 1999em data (purple squares) are from  and the SN 2004dj data (cyan circles) are from Leonard et al. (2006), reanalyzed in a manner similar to the current objects (see text for details). The SNe 2006ov (blue stars), 2006my (orange diamond), and 2007aa (red triangles) data points are from this work. The abscissa represents time in days relative to the end of the plateau, tp, for each object, as defined in the text. We used R-band photometry from Leonard et al. (2006) and Vinkó et al. (2006) to derive tp for SN 2004dj of 2004 Oct. 12.4, and R-band photometry from Leonard et al. (2002a) and Hamuy et al. (2001) were used to derive tp for SN 1999em of 2000 Feb 23.3. A systematic error contribution of 0.04% has been added in quadrature to all points to account for potential night-to-night instrumental variations, except for the SN 2004dj data, which already included a systematic error contribution. The black solid line is the well-sampled R-band light curve of SN 1999em to guide the eye. The data are consistent with a sharp increase in the continuum polarization in the last 10-20 days before the end of the plateau. epochs from day −5 through day 27. While the continuum polarization for SN 2007aa at early times is sufficiently small that the polarization angle is hard to reliably measure (a small amount of ISP could result in dramatic changes in the inferred angle), the long axis of the q − u loop seen in the Hα line on day −22 (Fig. 8) is approximately aligned with the late-time continuum polarization. The combination of low continuum polarization and significant line polarization (as seen in day −22 in SN 2007aa) may indicate that the Hα excitation is more sensitive to the underlying aspherical 56 Ni distribution than the continuum and hence the alignment with the late-time polarization is not accidental. In addition, the continuum polarization in SN 1999em remained at a constant P.A. over the course of the first 161 days after explosion ). The polarization angle variation of SN 2004dj appears to be an outlier compared to these other objects, or else the first observation after the end of the plateau was very fortunately timed to catch the polarization angle while it was rapidly varying.
The simplest explanation for the constancy of the polarization angles is that these SNe IIP have highly aspherical cores surrounded by hydrogen envelopes that are more spherical. Thus, we conclude that our late-time continuum polarization measurements are most easily understood as good tracers of the underlying asphericity of the photosphere in the core, and not of a "spotty photosphere" or other boundary effect that only produces a transient polarization signal when the photosphere enters the core. Observations taken well before the photosphere recedes through the boundary between the hydrogen-rich envelope and the helium-rich core already show signs of the asphericity that is present at late times.
Our maximum continuum polarization of 1.56% in SN 2006ov on day 0 implies a minimum axis ratio of 1.45:1 in the context of the aspherical electron-scattering atmospheres of Höflich (1991), if the object were viewed in the equatorial plane, and possibly higher if the symmetry axis were less inclined to the line of sight. Despite the lower continuum polarization observed in SN 2004dj, Leonard et al. (2006) derived a similar axis ratio for the core of that object after accounting for the preferred inclination of Chugai et al. (2005).
SN 2006ov also shows some evidence for inner asphericity in the late-time [O I] line profiles. The profiles of the λλ6300, 6364 doublet in our latest spectrum (day 51) are shown in Figure 13. Although the spectrum has not yet become fully nebular and the [O I] doublet does not yet completely dominate its portion of the spectrum, a double-peaked structure is apparent. The peaks in λ6300 are at velocities near 0 and 850 km s −1 . Confirmation  Li et al. (2007), while the day 51 spectrum was acquired on 2007 Feb. 14.655 with LRIS using a setup similar to that for the other spectra. Two notches (marked with arrows) are seen in the earliest SN 2006ov spectra near 6535 and 6590Å. At later times, none of the objects shows significant asymmetries in the Hα line profile. that this represents the oxygen distribution and not contamination by some other line can be seen from the profile of λ6364, which shows a peak near zero velocity as well as one to the red; however, the redder peak in λ6364 is closer to 500 km s −1 .
Double-peaked [O I] profiles have been seen in stripped-envelope supernovae and attributed to a toroidal oxygen distribution (Mazzali et al. 2005;Maeda et al. 2008;Modjaz et al. 2008), but the velocity splittings of the peaks in those objects are significantly higher. Alternative geometries are also possible (Taubenberger et al. 2009;Milisavljevic et al. 2010). Here, a more likely explanation is that we are seeing large-scale clumps or blobs in the inner oxygen distribution, especially given that the main peak is at zero velocity. Such oxygen clumps have been seen in stripped-envelope supernovae (Filippenko & Sargent 1989;Matheson et al. 2000;Silverman et al. 2009), SN 1987A (Stathakis et al. 1991, and the SN IIP 2002hh (Pozzo et al. 2006). The clumps are usually attributed to Rayleigh-Taylor instabilities created at the interface of the oxygen-rich and helium-rich layers in the core during the explosion (Müller et al. 1991), although other largescale fluid motions during the explosion are possible.

CONCLUSIONS
We have presented late-time spectropolarimetry of three SNe IIP. All three objects had large continuum polarizations after the end of the photometric plateau (0.8 P 1.6%), while SN 2007aa had contrastingly smaller continuum polarization on the plateau (P ≈ 0.1%). These results confirm that the polarization evolution of SN 2004dj (Leonard et al. 2006) was not unique and that highly aspherical cores exist in SNe IIP. This follows the trend identified by Wang et al. (2001) that the degree of polarization in core-collapse SNe generally  increases with time.
It is important to note that an intrinsic continuum polarization of 1.56% is quite high 5 . Although there are extreme cases such as the SN Ic 1997X ), very few core-collapse SNe exhibit intrinsic continuum polarizations higher than this. Most stripped-envelope SNe have significantly lower continuum polarization (after correction for ISP), including even extreme examples such as the well-studied broad-lined SN Ic 2002ap (Kawabata et al. 2002;Leonard et al. 2002b;Wang et al. 2003b). The line polarization features in some objects are stronger, however. SNe IIP are the single most common form of core-collapse supernova and represent the explosions of the most common massive stars, those at the low-mass end of the range of supernova progenitors. Our observations of these objects demonstrate that large degrees of asphericity are generic ingredients in even ordinary core-collapse explosions and are not just limited to various unusual or extreme objects.
This observational view parallels recent developments in the theoretical modeling of core-collapse SNe. Modern simulations of supernova explosions have revealed a rich phenomenology of multi-dimensional effects (Blondin et al. 2003;Burrows et al. 2006). Extreme models invoking jets (Khokhlov et al. 1999) have had some success explaining the polarization evolution and other tracers of asphericity in core-collapse SNe Couch et al. 2009), although the mechanism by which jets would be launched in ordinary SNe IIP remains obscure, as jet formation is usually thought to require atypically high angular momentum in the core (e.g., Burrows et al. 2007b). Moreover, the current generation of neutrino-powered corecollapse simulations (e.g., Blondin & Mezzacappa 2007;Burrows et al. 2007a;Marek & Janka 2009) all generate complex structures in the explosion without needing to invoke jets. In particular, these simulations commonly create aspherical shocks dominated by low-order deformations that produce unipolar or bipolar outflows that seem promising for generating the core asphericities we infer from our polarization measurements. It remains to be seen whether these features persist in the upcoming more-detailed three-dimensional models.
Unfortunately, most explosion models to date only simulate the central core of the star and stop when an outward-going shock is launched (or has failed to launch). However, the shock continues to evolve as it traverses the star and our observations can only probe the final state of the core, which makes detailed comparisons of the observed asphericities to those seen in the simulations difficult. Only a few studies (Kifonidis et al. 2006;Hammer et al. 2009) have followed the aspherical shock wave as it propagates outward through the star. Future work in this direction, along with radiative-transfer calculations, is necessary to see if the existing models can reproduce the observations.
It is interesting to note that SN 2006ov showed the largest continuum polarizations in our sample while it ejected the smallest amount of nickel and possibly had the lowest-mass progenitor. As mentioned above, Smartt et al. (2009)   ), yet it was a magnitude fainter than SN 2007aa at the beginning of the radioactive decay tail. The progenitor mass limit for SN 2006ov is also very low, M < 10 M ⊙ Crockett et al. 2009). In addition, the expansion velocities of SN 2006ov are lower than those of the other SNe IIP (Figure 12), possibly indicating a low explosion energy. A connection between the low nickel mass, low progenitor mass, or low explosion energy and the high continuum polarization would be exciting if confirmed in other objects.
We also encourage the enlargement of the published SN IIP spectropolarimetry sample, particularly of ob-jects lacking significant ISP. While observations at very late times are likely rare in the extant unpublished datasets, it would be fruitful to expand the number of objects plotted in Figure 11. Only then will we be able to look for relationships between the core asphericity and other supernova properties. Trends with progenitor mass, explosion energy, and 56 Ni production will all be interesting to examine and may provide constraints on the theoretical uncertainties in the core-collapse supernova mechanism.
We thank Douglas C. Leonard for supplying us with his SNe 1999em and 2004dj data and Ryan J. Foley, Mohan Ganeshalingam, Matthew Moore, and Thea N. Steele for their assistance with some of the observations. Most of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration; it was made possible by the generous financial support of the W. M. Keck Foundation. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Mauna Kea has always had within the indigenous Hawaiian community; we are most fortunate to have the opportunity to conduct observations from this mountain. We also would like to thank the expert assistance of the Keck and Lick staffs in making these observations possible. A.V.F.'s supernova group at U.C. Berkeley has been supported by NSF grants AST-0607485 and AST-0908886, as well as by the TABASGO Foundation. KAIT and its ongoing operation were made possible by donations from Sun Microsystems, Inc., the Hewlett-Packard Company, AutoScope Corporation, Lick Observatory, the NSF, the University of California, the Sylvia & Jim Katzman Foundation, and the TABASGO Foundation.