Quadrant polarization parameters for the scattered light of circumstellar disks Analysis of debris disk models and observations of HR 4796A

Context. Modern imaging polarimetry provides spatially resolved observations for many circumstellar disks and quantitative results for the measured polarization which can be used for comparisons with model calculations and for systematic studies of disk samples. Aims. This paper introduces the quadrant polarization parameters Q(000), Q(090), Q(180), Q(270) for Stokes Q and U-045, U-135, U-225, U-315 for Stokes U for circumstellar disks and describes their use for the polarimetric characterization of the dust in debris disks. Methods. We define the quadrant polarization parameters Q(xxx) and U-xxx and illustrate their properties with measurements of the debris disk around HR 4796A from Milli et al. (2019, A&A, 626, A54).. We calculate quadrant parameters for simple models of rotationally symmetric and optically thin debris disks and the results provide diagnostic diagrams for the determination of the scattering asymmetry of the dust. This method is tested with data for HR 4796A and compared with detailed scattering phase curve extractions in the literature. Results. The parameters Q(xxx) and U-xxx are ideal for a well-defined and simple characterization of the azimuthal dependence of the polarized light from circumstellar disk because they are based on the natural Stokes Q and U quadrant pattern produced by circumstellar scattering. For optically thin and rotationally symmetric debris disks the quadrant parameters normalized to the integrated azimuthal polarization Q(xxx)/Q(phi) and U-xxx/Q(phi) or quadrant ratios like Q(000)/Q(180) depend only on the disk inclination i and the polarized scattering phase function f(phi)(theta) of the dust, and they do not depend on the radial distribution of the scattering emissivity. Because the disk inclination i is usually well known for resolved observations, we can derive the shape of f(phi)(theta) for the phase angle range theta sampled by the polarization quadrants. This finding also applies to models with vertical extensions as observed for debris disks. Diagnostic diagrams are calculated for all normalized quadrant parameters and several quadrant ratios for the determination of the asymmetry parameter g of the polarized Henyey-Greenstein scattering phase function f(phi)(theta, g). We apply these diagrams to the measurement of HR 4796A, and find that a phase function with only one parameter does not reproduce the data well. We find a better solution with a three-parameter phase function f(phi)(theta, g(1), g(2), w), but it is also noted that the well-observed complex disk of HR 4796A cannot be described in full detail with the simple quadrant polarization parameters. Conclusions. The described quadrant polarization parameters are very useful for quantifying the azimuthal dependence of the scattering polarization of spatially resolved circumstellar disks illuminated by the central star. They provide a simple test of the deviations of the disk geometry from axisymmetry and can be used to constrain the scattering phase function for optically thin disks without detailed model fitting of disk images. The parameters are easy to derive from observations and model calculations and are therefore well suited to systematic studies of the dust scattering in circumstellar disks.

Automated summary

This study introduces quadrant polarization parameters for circumstellar disks, providing a simple characterization of azimuthal dependence of polarized light from these disks. These parameters can be used for systematic studies of dust scattering without detailed model fitting.