Revealing magnetic field structure at the surfaces of protoplanetary disks via near-infrared circular polarization

TL;DR

This paper proposes a new method using near-infrared circular polarization to probe magnetic field structures at the surfaces of protoplanetary disks, providing insights into disk wind mechanisms and dust properties.

Contribution

It introduces a novel observational technique to diagnose magnetic field configurations at disk surfaces through NIR circular polarization analysis.

Findings

Circular polarization can indicate magnetic field deviations from purely toroidal configurations.

Azimuthal dependence of polarization shows a double peaked profile for wind-driving fields.

Net circular polarization may be detectable even in unresolved disks.

Abstract

Context. Magnetic fields play a fundamental role in the dynamical evolution of protoplanetary disks, in particular via magnetically induced disk winds. The magnetic field structure at the disk surface is crucial for driving the disk winds; however, it is still poorly understood observationally. Aims. We explore a new method to probe the magnetic field structure at the disk surface using near-infrared (NIR) circular polarization. Near-infrared circular polarization arises when unpolarized stellar light is scattered by magnetically aligned grains at the disk surface. In this study, we aim to clarify to what extent the observed circular polarization pattern can be used to diagnose the magnetic field structure. Methods. We first calculated light scattering properties of aligned spheroids, and the results were then used to create expected observational images of the degree of circular polarization at a NIR wavelength. Results. Magnetically aligned grains can produce circular polarization, particularly when the field configuration deviates from a purely toroidal field. We find that disk azimuthal dependence of the degree of circular polarization tends to exhibit a double peaked profile when the field structure is favorable for driving disk winds by centrifugal force. We also find that even if the disk is spatially unresolved, a net circular polarization can possibly be nonzero. We also show that the amplitude of circular polarization is strongly dependent on grain composition and axis ratio. Conclusions. Our results suggest that circular polarization observations would be useful to study the magnetic field structure and dust properties at the disk surface.

Equations

$$\pmatrix{I}_{s}\\ Q_{s}\\ U_{s}\\ V_{s}\propto\pmatrix{S}_{11}&S_{12}S_{13}S_{14}\\ S_{21}S_{22}S_{23}S_{24}\\ S_{31}S_{32}S_{33}S_{34}\\ S_{41}S_{42}S_{43}S_{44}\pmatrix{I}_{i}\\ Q_{i}\\ U_{i}\\ V_{i}.$$

$$S_{11}=\frac{1}{2}(|S_1{|}^2+|S_2{|}^2+|S_3{|}^2+|S_4{|}^2),$$

$$S_{41}=\text{Im}\{S_2^{\ast }{S}_4+S_3^{\ast }{S}_1\},$$

$$\frac{h}{r}=0.1617\cdot {\left(\frac{r}{100\text{leavevmode}\mathrm{au}}\right)}^{0.219}.$$

$${⟨P_{\text{circular}}⟩}_{\varphi }=\frac{1}{2\pi }{\int }_0^{2\pi }{P}_{\text{circular}}(\varphi )d\varphi ,$$

Full text

*Aims. *We explore a new method to probe the magnetic field structure at the disk surface using near-infrared (NIR) circular polarization. Near-infrared circular polarization arises when unpolarized stellar light is scattered by magnetically aligned grains at the disk surface. In this study, we aim to clarify to what extent the observed circular polarization pattern can be used to diagnose the magnetic field structure.

*Methods. *We first calculated light scattering properties of aligned spheroids, and the results were then used to create expected observational images of the degree of circular polarization at a NIR wavelength.

*Results. *Magnetically aligned grains can produce circular polarization, particularly when the field configuration deviates from a purely toroidal field. We find that disk azimuthal dependence of the degree of circular polarization tends to exhibit a double peaked profile when the field structure is favorable for driving disk winds by centrifugal force. We also find that even if the disk is spatially unresolved, a net circular polarization can possibly be nonzero. We also show that the amplitude of circular polarization is strongly dependent on grain composition and axis ratio.

*Conclusions. *Our results suggest that circular polarization observations would be useful to study the magnetic field structure and dust properties at the disk surface.

1 Introduction

Magnetic fields are a crucial factor in determining how protoplanetary disks form, evolve, and disperse. Recent studies suggest that magnetohydrodynamical (MHD) disk winds play an important role in disk evolution. The occurrence of MHD disk winds is strongly dependent on the magnetic field structure at the disk surface. A poloidal magnetic field structure with a large angle with respect to the disk normal would give rise to disk winds via centrifugal force (Blandford & Payne, 1982), while a more toroidal field would drive winds through a magnetic pressure gradient (e.g., Bai, 2017). Therefore, to study how MHD disk winds are launched, it is essential to understand the magnetic field geometry at the disk surface.

Observing the magnetic field in protoplanetary disks is a notoriously difficult task. Polarized thermal emission from magnetically aligned grains has been the most successful technique in the interstellar medium (Lazarian, 2007; Andersson et al., 2015). However, this technique is not always reliable in protoplanetary disks, because larger grains that are responsible for thermal emission are not necessarily aligned with the magnetic field (Tazaki et al., 2017; Yang, 2021). Alternatively, molecular line polarization might contain information about both $B$-field orientation and strength, but requires exquisite sensitivity, which can only marginally be reached with current telescopes (Teague et al., 2021; Harrison et al., 2021; Stephens et al., 2020). Recently, Yang & Li (2022) explored a new technique, finding that linear polarization angles of near-infrared (NIR) disk scattered light carry information on the magnetic field, as long as micron-sized grains at the disk surfaces are magnetically aligned. One of the great advantages of this technique is its sensitivity to the disk-surface magnetic field, which is crucial for understanding MHD disk winds.

Although Yang & Li (2022) mainly focused on linearly polarized scattered light arising from magnetically aligned grains, we can also expect circular polarization (CP) of scattered light (e.g., Bandermann & Kemp, 1973; Fukushima et al., 2020). However, there have been no quantitative studies on how disk-scattered light is circularly polarized at near-IR and to what extent it can be used for diagnosing disk magnetic fields.

In this study, we aim to trace the magnetic field structure in protoplanetary disks using circular polarized scattered light. To achieve this goal, we will create synthetic maps showing the degree of circular polarization throughout a protoplanetary disk and study how and to what extent the $B$-field is reflected in the circular polarization pattern.

This paper is organized as follows. We introduce our method and the dust model in Sect. 2. In Sect. 3, we present the CP patterns. We discuss the implications for future observations in Sect. 4. Finally, we conclude in Sect. 5.

2 Methods and models

To calculate the CP of disk-scattered light in the NIR, we follow a similar method to Yang & Li (2022), where the authors calculated linearly polarized scattered light by aligned grains at the disk surface. In this study, we assume single scattering; the stellar radiation will be scattered by a grain at the disk surface only once before reaching the observer.

2.1 Circular polarization calculation

At the disk surface, dust particles scatter unpolarized stellar light to the observer. The scattered light is in general linearly and circularly polarized. The polarization of scattered light is described by the Stokes parameter $(I_{s},Q_{s},U_{s},V_{s})$, which can be calculated from the Stokes parameter of incident light ($I_{i},Q_{i},U_{i},V_{i}$) and the scattering matrix (Bohren & Huffman, 1983):

[TABLE]

The scattering matrix elements $S_{ij}$ describe how the polarization state of light changes upon scattering by a grain. Considering unpolarized incident light from the star, the degree of circular polarization in scattered light $V_{s}/I_{s}$ is equal to $S_{41}$/$S_{11}$. The scattering matrix elements are given by (Bohren & Huffman, 1983)

[TABLE]

where $S_{i}\,(i=1,2,3,4)$ represent the scattering amplitude. To calculate Eqs. (2) and (3), we assume the dipole approximation (Rayleigh limit), where a grain has to be small compared to the wavelength. We also assume oblate dust grains with an axis ratio of 1.5. With these assumptions, we can find analytical expressions for the scattering amplitude $S_{i}$ (Bohren & Huffman, 1983). In the Rayleigh limit, $S_{41}/S_{11}$ does not depend on the size parameter, which is determined by grain size and wavelength, but depends solely on dust composition (refractive index) and the axis ratio of the oblate grain. The validity of the dipole approximation is discussed in Sect. 4.2.

We assume a dust composition model proposed by Birnstiel et al. (2018), where each grain consists of a mixture of water ice (Warren & Brandt, 2008), carbonaceous material, troilite (Henning & Stognienko, 1996), and astronomical silicate (Draine, 2003). For carbonaceous material, we consider two possibilities: one is refractory organics (Henning & Stognienko, 1996), as with Birnstiel et al. (2018), and the other is amorphous carbon (Zubko et al., 1996). Although dust composition including refractory organics has been commonly used in the literature (Pollack et al., 1996; D’Alessio et al., 2001; Kataoka et al., 2014; Birnstiel et al., 2018), a composition including amorphous carbon has also been used to model disk observational quantities (e.g., DIANA opacity: Woitke et al., 2016). The detailed effect of dust composition on CP is discussed in Sect. 3.4. Using optool (Dominik et al., 2021), we derived an effective refractive index at $\lambda=1.6\,\rm{\mu m}$ using the Bruggeman mixing rule (Bruggeman, 1935). As a result, we find a refractive index of $m=1.5625+i0.0204$ and $m=2.0217+i0.4308$ for a mixture comprising refractory organics and amorphous carbon, respectively.

2.2 Disk model and magnetic field geometries

To model the disk-scattering surface at NIR wavelengths, we assumed an axisymmetric flared disk. We determined the disk flaring geometry motivated by NIR observations (Avenhaus et al., 2018). By fitting ring structures seen in scattered light images, Avenhaus et al. (2018) determined the heights of disk scattering surfaces. As a result, they obtained an empirical fit to the ratio of the height of the scattering surface $h$ and the midplane distance from the star $r$ as follows:

[TABLE]

At the disk surface, the oblate grains are aligned with the magnetic field. For the sake of simplicity, we assume that grains are perfectly aligned with their short axis parallel to the local magnetic field direction. We introduce two parameters ($\theta_{B},\phi_{B}$) to specify the local magnetic field direction, as done in Yang & Li (2022) (see Fig. 1). For example, when $\theta_{B}=0^{\circ}$, the magnetic field is vertical. When $\theta_{B}=90^{\circ}$ and $\phi_{B}=0^{\circ}$, the field is radial. When $\theta_{B}=90^{\circ}$ and $\phi_{B}=90^{\circ}$, the field is azimuthal (toroidal). By considering the range of [$0^{\circ},90^{\circ}$] for both $\theta_{B}$ and $\phi_{B}$, we can mimic various magnetic field structures at the disk surface.

Given the disk geometry, the orientation of aligned oblate grains, and the inclination angle, we can calculate scattering angles at each location of the disk surface (see e.g., Stolker et al., 2016, for the case of spherical grains). By applying these scattering angles to the CP calculation method described in Sect. 2.1, we can make a projected disk image of the degree of CP. We also calculated the linear polarization of scattered light and confirmed that our code can successfully reproduce the results presented in Yang & Li (2022).

3 Results

We calculated the degree of circular polarization produced by aligned dust grains at the disk surface for a set of assumed magnetic field structures and analyzed the imprint of the $B$-field on the CP pattern.

3.1 Integrated CP values

We show two synthetic maps of the degree of circular polarization projected at the surface of the disk model in Fig. 2. For both images, we assumed an inclination of $i=55^{\circ}$, and polar magnetic field angle $\theta_{B}=45^{\circ}$, but we assumed different values of $\phi_{B}$. The upper panel represents a disk with $\phi_{B}=0^{\circ}$, which corresponds to a purely poloidal field. The CP pattern that arises from this disk surface is mirror symmetric (except for the sign) with respect to the disk minor axis. In the lower panel, we added a toroidal component to the magnetic field by increasing $\phi_{B}$ to $60^{\circ}$. In this situation, the mirror symmetry is broken, which demonstrates the effect of a different magnetic field structure on the CP pattern.

To further analyze the disk CP patterns, we calculated the azimuthally averaged degree of circular polarization, given by

[TABLE]

where $\phi$ is the disk azimuthal angle. We examined the integrated CP for a $50\times 50$ grid of $\phi_{B}$ and $\theta_{B}$ values, assuming $i=55^{\circ}$ and the dust composition with amorphous carbon. Figure 3 shows the azimuthally averaged CP values for various magnetic field configurations. The majority of possible magnetic field angles produce an integrated negative degree of CP, although this could also be positive if we assume negative $\phi_{B}$ values. The averaged CP amplitude tends to become larger for larger $\theta_{B}$ at intermediate $\phi_{B}$ values. This suggests that even if the disk is spatially unresolved, CP signals can remain nonzero, particularly when the magnetic field makes a large angle with the $z$-axis.

The averaged CP percentage is zero in some situations: when the magnetic field is vertical ($\theta_{B}=0^{\circ}$), poloidal ($\phi_{B}=0^{\circ}$), or toroidal ($\theta_{B}=90^{\circ}$, $\phi_{B}=90^{\circ}$). This can be caused by either an averaging effect or no intrinsic CP production at all. To distinguish which of these two causes is responsible, Fig. 4 shows the nonaveraged maximum absolute value of CP at a disk surface ring of $h/r=0.2$. The nonzero maximum values for $\phi_{B}=0^{\circ}$ or $\theta_{B}=0^{\circ}$ reveal that oblates aligned with these magnetic fields produce a symmetric CP pattern, where one side of the disk has negative CP values and the other side has positive (see Fig. 2, upper panel). As a result, CP signals cancel out when they are averaged. When the field is purely toroidal ($\theta_{B}=90^{\circ};\phi_{B}=90^{\circ}$), the maximum CP is zero. Alignment with this magnetic field structure results in the symmetry axis of the oblate being perpendicular to the incident light, which results in no intrinsic circular polarization (Bandermann & Kemp, 1973). Because of this, we conclude that the toroidal magnetic field seems to be better probed by the linear polarization technique instead of the CP-based technique (see Sect. 4.1). In Fig. 4, we also find that there is another case where the intrinsic CP value becomes zero, that is, when $\theta_{B}=78^{\circ}$ and $\phi_{B}=0^{\circ}$. In this case, the incident light is parallel to the symmetry axis of the oblate grain. As a result, the $S_{41}$ component vanishes, and therefore no circular polarization is produced (Bandermann & Kemp, 1973).

3.2 Azimuthal profiles

To study the effect of the $B$-field on the CP pattern in more detail, we created azimuthal profiles. These profiles show the CP of a disk surface ring at $h/r=0.2$, measured counterclockwise from the near side of the disk. After evaluating all possible magnetic field angles, we distinguished three different shapes of azimuthal CP dependence: a dip followed by a peak (“regular”), a peak followed by a dip (“reversed”), and two dips (“double peaked”). We present examples for these three profile shapes in Fig. 5, where we assume the dust composition with amorphous carbon and an inclination of $i=55^{\circ}$. Because we find that the reversed profile occurs in only a relatively small part of the parameter space, we mainly focus on the regular and doubled-peaked profiles hereafter.

We find that the regular and double-peaked azimuthal CP profiles can place additional constraints on the magnetic field configuration, as shown in Figs. 3 and 4. In those plots, we overlaid the distribution of the three azimuthal profile shapes. In addition, we indicated the boundary of magnetic field structures at which MHD disk winds by centrifugal force might occur with a blue dashed line. In order to launch a wind along the magnetic field lines, the centrifugal force needs to exceed the stellar gravity. This requirement is met when $\sin\theta_{B}\geq 1/\sqrt{3\cos^{2}\phi_{B}+1}$. When $\phi_{B}=0$, we obtain $\theta_{B}=30^{\circ}$, which agrees with the result of Blandford & Payne (1982). Magnetic field configurations at the right side of the line have field lines that make a sufficiently large angle with the $z$-axis to give rise to MHD disk winds by centrifugal force. Almost all double-peaked CP profiles are situated inside the disk-wind boundary, which suggests that these types of profiles are produced by magnetic field lines that make a large angle with the $z$-axis.

3.3 Effect of disk inclination

In order to study the impact of disk inclination angle, we also compare CP at $i=35^{\circ}$ and $75^{\circ}$ in Fig. 6. The effect of disk inclination can be seen in the differences between the boundaries separating the regular and double-peaked profiles. For the case of $i=35^{\circ}$, the double-peaked profile occurs even if $\theta_{B}$ is as small as $40^{\circ}$. For the case of $i=75^{\circ}$, the double-peaked profiles only occur at even higher $\theta_{B}$ values than in a lower inclination case. Therefore, as long as the disk is moderately or highly inclined, the occurrence of the double-peaked profile indicates a magnetic field configuration with a high $\theta_{B}$ value that is favorable for launching winds by centrifugal force.

3.4 Effect of other dust-grain parameters

We examine how different dust models would affect our results by evaluating the CP patterns for another dust composition and other axis ratios. In the alternative dust composition, we replaced amorphous carbon with refractory organics (Birnstiel et al., 2018). We find that this less absorbing material results in less circular polarized light, with a maximum value of $\sim\mkern-4.0mu0.4\%$ compared to $\sim\mkern-4.0mu6\%$ for the mixture with amorphous carbon. These lower CP values can be physically explained as follows. Circular polarization is produced when there is a phase lag between dipole radiation along the major and minor axis of the grain. The amount of the phase lag is related to the absorption efficiency and becomes zero when grains are nonabsorbing (Bandermann & Kemp, 1973). Consequently, less absorbing material results in a smaller phase lag and therefore weaker CP production. Lower axis ratios also result in less CP because circular polarization arises from the difference in polarizability along each spheroid axis, which is smaller for more spherical particles. Fortunately, although both grain composition and axis ratio affect the absolute value of CP, the shape of the azimuthal profiles does not change. This demonstrates that the occurrence of double-peaked profiles is independent of the dust composition and grain-axis ratio, and is therefore a reliable probe of the magnetic field structure, as long as the Rayleigh approximation holds. Therefore, if we detect CP, its amplitude would be a useful diagnostic of grain properties.

As long as we work on the Rayleigh limit, the effects of axis ratio and dust composition on the CP amplitude are degenerate. To better constrain dust properties, a combined analysis with the degree of CP and disk-scattered flux would be useful. This is because the imaginary part of the refractive index is not only important for CP, but also for the single scattering albedo, which in turn strongly affects disk-scattered flux.

4 Discussion

4.1 Role of linear and circular polarization in probing disk magnetic fields

Our study, together with previous studies, demonstrates that polarization signals in scattered optical and NIR light contain valuable information about the $B$-field in protoplanetary disks. Yang & Li (2022) explored how linearly polarized scattered light can be used to diagnose the magnetic field. Here, we discuss how two types of polarization observations (linear and circular) can complement each other and advance our understanding of disk magnetic fields.

It would be very useful to combine the linear and circular polarization analysis because this allows us to exclude polarization signals that are produced via pathways other than magnetically aligned grains, such as multiple scattering (Canovas et al., 2015) or instrumental errors. Observing both linear and circular polarization provides an independent check as to whether or not the polarization signals are truly produced by aligned grains.

In addition, linear and circular polarization in disk-scattered NIR light are sensitive to different magnetic field components. Linear polarization angles seem more suitable to indicate the presence of the toroidal component, whereas circular polarization profiles better probe the poloidal component with a large angle to the $z$-axis. The former $B$-field configuration drives the wind by a magnetic pressure gradient (e.g., Bai, 2017), while the latter does this via centrifugal force (Blandford & Payne, 1982). Therefore, the combination of these could narrow down the possible magnetic field structures and eventually reveal a launching mechanism of MHD winds.

To study the linear and circular polarization signals as described in the techniques of our study and of Yang & Li (2022), the disk needs to be resolved. However, the azimuthally averaged CP values can be measured without resolving the disk. We find that an integrated nonzero degree of circular polarization already points towards the presence of aligned grains at the disk surface, which is an essential requirement to study the magnetic field structure using linear and circular polarization. The integrated CP value can reveal disks that meet this requirement, which helps to search for suitable targets.

4.2 Circular polarization from large aligned grains

In this study, we adopt the Rayleigh approximation, which is only applicable for grains smaller than the wavelength. At a wavelength of $\lambda=1.6\leavevmode\nobreak\ \mu$m, the Rayleigh approximation is applicable when the grain radius is smaller than $\lambda/2\pi\sim 0.25\leavevmode\nobreak\ \mu$m. Here, we discuss how the presence of larger grains would affect the CP and whether or not the Rayleigh limit is a valid assumption.

Fukushima et al. (2020) calculated circularly polarized scattered light by aligned grains in interstellar slab-like clouds. These authors found that increasing grain radius results in increasing the degree of CP. They also found that the CP pattern outside the Rayleigh limit would become less simple than in the Rayleigh limit (see their Fig. 4). Therefore, the presence of micron-sized grains at the disk surface would make much more complicated CP patterns than those presented in this study.

However, one caveat of Fukushima et al. (2020) is that the authors only considered compact micron-sized grains. At the disk surface, micron-sized dust particles at the disk surface would be aggregates (e.g., Dominik & Tielens, 1997; Blum & Wurm, 2008), which are clusters of small single particles (monomers) glued together by intermolecular forces, such as the van der Waals force. Although these studies suggest that aggregates may have sizes beyond the Rayleigh limit, the degree of polarization of scattered light by large aggregates rather resembles that of individual monomers (e.g., West, 1991; Kozasa et al., 1993; Volten et al., 2007; Tazaki & Dominik, 2022). Tazaki & Dominik (2022) estimated the upper limit on the monomer radius to be $0.4\,\rm{\mu m}$ for a set of protoplanetary disks based on polarimetric observations. Therefore, as long as such dust aggregates are the dominant population of the disk surface, and monomers in the aggregate are aligned, our argument would remain approximately valid. Another way to circumvent possible deviations from the Rayleigh solution would be to observe the disk at longer wavelengths, such as $2.2\,\mu$m ($K$-band).

4.3 Do grains align with magnetic fields?

In the above sections, we assume that grains at the disk surface perfectly align with magnetic fields. Here, we discuss to what extent this assumption is reasonable.

Tazaki et al. (2017) and Yang & Li (2022) suggest that grain alignment with the magnetic field, also known as $B$-RAT, is feasible at the disk surfaces when grains contain a sufficient amount of superparamagnetic inclusions. However, $B$-RAT is not the only known mechanism to align grains. Dust grains at the disk surface can also become aligned with the stellar radiation, which is called $k$-RAT (Lazarian & Hoang, 2007). $k$-RAT alignment predicts that the grain’s minor axis becomes parallel to the incoming radiation direction. In this case, we expect no circular polarization as long as the grain polarization property is approximated by the Rayleigh limit. This is because no CP is produced when the particle symmetry axis is either in the scattering plane ($k$-RAT for oblate grains) or perpendicular to the incident light direction ($k$-RAT for prolate grains) (Bandermann & Kemp, 1973).

Therefore, if CP observations end up with a nondetection, this could be attributed to $k$-RAT alignment, although other reasons still exist, such as very low alignment efficiency, $B$-RAT alignment with a symmetric $B$-field structure, and so on. Conversely, if we detect a CP signal, this is likely arising from $B$-RAT alignment. In this way, CP observations will shed light on grain-alignment processes at the disk surface.

It is also worth mentioning that CP can be produced by multiple scattering from nonaligned or spherical grains. However, such a CP signal will efficiently vanish when it is integrated over the disk as long as we can assume an axisymmetric disk structure. In contrast, as we discuss in Sect. 3.1, an integrated CP signal can be nonzero if it is due to aligned grains. Therefore, if we detect a net nonzero integrated CP from an axisymmetric disk, it is likely attributable to grains aligned with the $B$-field rather than the multiple-scattering origin.

4.4 Future prospects

The next step in this study would be to compare our results with direct observations of the degree of circular polarization in protoplanetary disks. No direct circular polarization measurements of protoplanetary disks have been published yet, but promising instruments and techniques are emerging. The Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument on the Subaru telescope will allow us to obtain circular polarization images of resolved protoplanetary disks at NIR wavelengths. In addition, although not originally designed for circular polarization measurements, VLT/SPHERE could also be used to observe CP, as proposed by van Holstein et al. (2020).

One caveat of our model is that we assume single scattering to simplify the disk radiative transfer processes. However, protoplanetary disks are optically thick, and in general, multiple scattering plays a role, which can affect the degree of polarization (e.g., Ma & Schmid, 2022). In addition, aligned grains also produce nonazimuthal linear polarization (Yang & Li, 2022) which can be converted into circular polarization via multiple scattering. Further clarification of these effects is left for our future study.

5 Summary

In this study, we explored a new method to study the magnetic field structure at the surface of protoplanetary disks using NIR circular polarization of scattered light. Assuming a disk model and various magnetic field geometries, we created synthetic maps of the degree of circular polarization produced by scattering oblate spheroids at the disk surface. By analyzing the azimuthal dependence and the integrated value, we studied to what extent the configuration of the $B$-field is reflected in the CP patterns. We also investigated the effects of disk inclination and dust grain parameters on our results. We examined what we can learn from the CP disk images about the launching mechanism of MHD disk winds, and how this method ties in with linear polarization techniques. Our primary findings are as follows.

• •

It is shown that magnetically aligned grains at the disk surface can produce circular polarization signals. Also, circular polarization signals remain nonzero even if the disk is not spatially resolved.

• •

It is found that a disk exhibiting a CP pattern with a double-peaked azimuthal profile likely has a poloidal magnetic field with a large angle with respect to the $z$-axis. In contrast, circular polarization becomes zero for a toroidal configuration.

• •

The degree of CP strongly depends on dust composition. A more absorbing composition, such as when grains are comprised of amorphous carbon, produces a higher degree of CP. The composition does not change the CP azimuthal profiles, as long as the Rayleigh approximation holds.

• •

We find that circular polarization is complementary to the linear polarization technique to probe the disk $B$-field because the combined analysis allows us to break degeneracy, and the polarization signals trace different magnetic field components.

Acknowledgements.

We thank an anonymous referee for useful comments. In addition, we thank C. Dominik and J. Hashimoto for a fruitful discussion. R.T. acknowledges the JSPS overseas research fellowship.