Quadrant polarization parameters for the scattered light of circumstellar disks

Analysis of debris disk models and observations of HR 4796A

ETH Zurich, Institute for Particle Physics and Astrophysics, Wolfgang-Pauli-Strasse 27, 8093 Zurich, Switzerland
e-mail: schmid@astro.phys.ethz.ch

Received: 23 January 2021
Accepted: 19 August 2021

Abstract

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₀₀₀, Q₀₉₀, Q₁₈₀, Q₂₇₀ for Stokes Q and U₀₄₅, U₁₃₅, U₂₂₅, U₃₁₅ 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_ϕ and U_xxx∕Q_ϕ or quadrant ratios like Q₀₀₀∕Q₁₈₀ depend only on the disk inclination i and the polarized scattering phase function f_ϕ(θ) 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_ϕ(θ) for the phase angle range θ 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_ϕ(θ, 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_ϕ(θ, g₁, g₂, 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.