# Interplay of dust alignment, grain growth, and magnetic fields in polarization: lessons from the emission-to-extinction ratio

L. Fanciullo, V. Guillet, Francois Boulanger, A. Jones

► To cite this version:

L. Fanciullo, V. Guillet, Francois Boulanger, A. Jones. Interplay of dust alignment, grain growth, and magnetic fields in polarization: lessons from the emission-to-extinction ratio. Astronomy & Astrophysics - A&A, 2017, 602, pp.A7. 10.1051/0004-6361/201630373. hal-02022165

HAL Id: hal-02022165

https://hal.science/hal-02022165v1

Submitted on 26 Jun 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Interplay of dust alignment, grain growth, and magnetic fields in polarization: lessons from the emission-to-extinction ratio

L. Fanciullo, V. Guillet, F. Boulanger, and A. P. Jones

Institut d'Astrophysique Spatiale (IAS), Bâtiment 121, Université Paris-Sud 11 and CNRS, 91405 Orsay, France
e-mail: lapo.fanciullo@ias.u-psud.fr

Received 29 December 2016 / Accepted 27 January 2017

## ABSTRACT

**Context.** Polarized extinction and emission from dust in the interstellar medium (ISM) are hard to interpret, as their dependence on dust optical properties, grain alignment, and magnetic field orientation is complex. This is particularly true in molecular clouds. The aforementioned phenomena are usually considered independently in polarization studies, while it is likely that they all contribute and their effects have yet to be disentangled.

**Aims.** The data available today are not yet used to their full potential. The combination of emission and extinction, in particular, provides information not available from either of them alone. We combine data from the scientific literature on polarized dust extinction with Planck data on polarized emission, and we use them to constrain the possible variations in dust and environmental conditions inside molecular clouds, and especially translucent lines of sight, taking the magnetic field orientation into account.

**Methods.** We focused on the dependence between $\lambda_{max}$ (the wavelength of maximum polarization in extinction) and other observables such as the extinction polarization, the emission polarization, and the ratio between the two. We set out to reproduce these correlations using Monte Carlo simulations in which we varied the relevant quantities in a dust model, which are grain alignment, size distribution, and magnetic field orientation, to mimic the diverse conditions that are expected inside molecular clouds.

**Results.** None of the quantities we chose can explain the observational data on their own: the best results are obtained when all quantities vary significantly across and within clouds. However, some of the data, most notably the stars with a low ratio of polarization in emission to polarization in extinction, are not reproduced by our simulation.

**Conclusions.** Our results suggest not only that dust evolution is necessary to explain polarization in molecular clouds, but that a simple change in size distribution is not sufficient to explain the data. Our results also point the way for future and more sophisticated models.

**Key words.** polarization – ISM: clouds – dust, extinction – evolution – submillimeter: ISM

## 1. Introduction

The light of stars often shows a degree of polarization correlated to interstellar extinction, up to a degree of a few percent per $A_V$ magnitude. This phenomenon has long been recognized as the effect of non-spherical cosmic dust grains that are aligned with interstellar magnetic field lines (Hall 1949; Hiltner 1949). Dust extinguishes starlight, and the component of the electric field parallel to the longer axis of a grain is more extinguished than the orthogonal field component. Furthermore, interstellar dust grains align their shorter axes with the interstellar magnetic field, so that they are not generally randomly oriented (e.g., Andersson et al. 2015, and references therein). The overall result is that the dusty magnetized interstellar medium (ISM) polarizes the starlight that was not originally polarized. The polarization fraction $p$ and the polarization angle $\psi$ of starlight therefore provide information on both interstellar dust and the Galactic magnetic field, or at least the component of the field that is parallel to the plane of the sky.

It is mainly the large grains that are aligned (e.g. Kim & Martin 1995) and, in typical ISM conditions, their thermal emission mainly falls in the far-infrared (FIR) and submillimeter (submm) range. This emission is also polarized, since emission is more efficient for the electric field component parallel to the longer axis, and it is an important complement to observations of polarized extinction in the optical and near-infrared

(NIR; e.g., Hildebrand 1988; Planck Collaboration Int. XIX 2015). It should be noted that since radiation polarized parallel to the longer axis of a grain is least intense in extinction and most intense in emission, we expect $\psi$ in the submm to be orthogonal to $\psi$ in the optical.

The main factors that determine the polarization fraction $p$ are the optical properties of the dust, the alignment efficiency, and the orientation of the magnetic field lines. Polarization can be expressed as (Lee & Draine 1985)

$$p = p_0 R \cos^2 \gamma,$$

where $p_0$ is the highest possible polarization given the dust properties, the parameter $R$, which is between 0 and 1, accounts for the effects of imperfect alignment<fnref n="1" />, and $\gamma$ is the angle between the magnetic field lines and the plane of the sky.

The wavelength dependence of polarization in extinction usually follows the so-called Serkowski curve (Serkowski et al. 1975),

$$p(\lambda) = p_{\max} \cdot \exp(-K \cdot \ln(\lambda/\lambda_{\max})^2),$$

<footnote n="1">¹ When grains are in the Rayleigh regime, as is the case for the thermal emission of large dust grains, the parameter $R$ can be calculated analytically, and it is called the Rayleigh reduction factor (Greenberg 1968, p. 328).</footnote>
where the polarization has a maximum $p_{max}$ at a wavelength $\lambda_{max}$, usually falling in the visible; the value of the parameter K, tied to the inverse of the FWHM, is usually around unity. Since the polarization efficiency of a grain peaks at ~2π times its size (Kim & Martin 1995), $\lambda_{max}$ traces the typical size of aligned grains: variations in $\lambda_{max}$ between lines of sight may indicate a change in grain size distribution, in the dependence of alignment on size, or both (e.g. Andersson & Potter 2007, hereafter AP07). The problem is further complicated because $\lambda_{max}$ also shows some dependence on the magnetic field angle $\gamma$ (Voshchinnikov et al. 2016).

Since polarization depends on many factors at once, interpreting it is a degenerate and difficult problem. This is especially true of the dense and complex environments that are molecular clouds, where magnetic field orientation, grain alignment, and dust properties are expected to change on small scales. Despite this, studies on dust polarization are often focused on constraining the grain alignment efficiency (e.g. AP07) or the structure of the magnetic field (e.g. Planck Collaboration Int. XX 2015) without accounting for other variables. One way of confronting the problem is to construct a cloud model that includes dust evolution, magnetic field structure and grain alignment, but such models are complex and very demanding computationally. The judicious combination of observational data can also provide interesting insights on dust physics, while requiring far lighter calculations.

One example of such labor-saving combinations is the complementary use of extinction data and polarized dust thermal emission. This last has the dimensions of an intensity and it is usually observed in the FIR and submm. All-sky surveys in the submm like Planck (Planck Collaboration Int. XIX 2015) are opening up new possibilities for this type of multiwavelength analysis. The idea that extinction and emission combined can be more informative than either of them alone is explored in Planck Collaboration Int. XXI (2015), which examines the $P_{353}/p_V$ ratio<fnref n="2" /> in the diffuse ISM, $P_{353}$ being the polarized intensity in emission in the Planck 353 GHz channel, and $p_V$ being the starlight polarization degree in the V band. This ratio is measured in MJy sr$^{-1}$, and since $P_{353}$ and $p_V$ have (at first approximation) the same dependence on alignment and $\gamma$, it should provide strong constraints on the properties of aligned dust grains. Among the new results made possible by Planck is the determination of average $P_{353}/p_V$ in the diffuse ISM: 5.4 MJy sr$^{-1}$, which is about 2.5 times higher than predicted by previous dust models (Planck Collaboration Int. XXI 2015).

The present paper aims to extend the study of $P_{353}/p_V$ to denser environments, namely translucent lines of sight in molecular clouds, using an updated method and a dust model optimized for the high $P_{353}/p_V$ found by Planck (see Guillet et al. 2016). The paper is organized as follows: Sect. 2 introduces the observational data we used, mostly from translucent lines of sight, and the selection that had to be made so that the comparison of extinction and emission would be meaningful. Section 3 presents the dust model we used, which was especially constructed to reproduce the $P_{353}/p_V$ ratio and classic dust observables. The dust model was created for the diffuse ISM, while we study translucent clouds where dust evolution may be taking place: the modification that is needed for the model to fit the data suggest what the nature of dust evolution might be in these areas. The model results are compared to observations in Sect. 4, and the meaning of this comparison is discussed in Sect. 5. Finally, Sect. 6 contains the conclusions and perspectives.

<table><caption>**Table 1.** Molecular clouds with their approximate position on the sky in Galactic coordinates.</caption><thead><tr><th>Cloud</th><th><i>l</i></th><th><i>b</i></th><th><i>D</i> (pc)</th><th>References</th></tr></thead><tbody><tr><td>Chamaeleon</td><td>297°</td><td>-15°.5</td><td>120-150</td><td>1, 2, 9</td></tr><tr><td>Musca</td><td>301°</td><td>-8°.0</td><td>120-150</td><td>3, 9</td></tr><tr><td>Ophiuchus</td><td>354°</td><td>15°.0</td><td>120</td><td>1, 4, 5, 7</td></tr><tr><td>R CrA</td><td>0°</td><td>-19°.5</td><td>130</td><td>1, 8</td></tr><tr><td>Taurus</td><td>174°</td><td>-14°.0</td><td>140</td><td>1, 6, 7</td></tr><tr><td>"</td><td>168°.5</td><td>-16°.5</td><td>"</td><td>"</td></tr></tbody></table>

**Notes.** Since the Taurus cloud is elongated, stars were sampled from around two different centers. The clouds include most of the lines of sight in this study but not those from Anderson et al. (1996).

**References.** (1) Whittet et al. (1992); (2) Covino et al. (1997); (3) Arnal et al. (1993); (4) Vrba et al. (1993); (5) Martin et al. (1992); (6) Whittet et al. (2001); (7) Loinard (2013); (8) Neuhäuser & Forbrich (2008); (9) Corradi et al. (2004). References in bold provided only the distance estimate and not the polarimetric data.

## 2. Data

Our work combines measures of starlight polarization in the near-ultraviolet (NUV) to NIR range that we recovered from the literature, and Planck measurements of total and polarized dust emission at 353 GHz for the same lines of sight. We use a total of 132 objects, which are reduced to 70 after selection (Sect. 2.3).

### 2.1. Extinction in the NUV to NIR

Most of our data points are from AP07, who include data on several clouds from the literature (Chamaeleon, Musca, Ophiuchus, R Coronae Australis (RCrA), and Taurus, Covino et al. 1997; Arnal et al. 1993; Vrba et al. 1993; Whittet et al. 1992, 2001, see Table 1) in a polarimetric study of the Coalsack nebula. The AP07 study is particularly useful for this purpose because the authors did not employ the Serkowski fit results from the literature, but they used the photometric and polarimetric data therein to conduct their own fits, thus minimizing the systematic effects from different fitting procedures. To increase the statistics, we also included data from Martin et al. (1992), who provide more lines of sight in Ophiuchus, and Anderson et al. (1996), who provide lines of sight that are for the most part<fnref n="3" /> not associated with the aforementioned clouds.

From the literature we obtained the Serkowski polarization parameters $p_{max}$ and $\lambda_{max}$ for all stars, as well as the polarization angle $\psi_{ext}$. The values of K are not calculated consistently in the literature: AP07 often use $K = 1.15$ and only fit K as a free parameter if it constitutes a statistically significant improvement, while Anderson et al. (1996) impose that K be a linear function of $\lambda_{max}$. For this reason, we chose not to include K in our work. The AP07 data retrieved from the Vizier online database<fnref n="4" /> do not include the uncertainties on $\psi_{ext}$ for Ophiuchus, therefore we complemented the data with the original article (Vrba et al. 1993): we used the average on $\psi$ in the various bands as the value of $\psi_{ext}$ and their standard deviation as the uncertainty. We excluded stars with standard deviations greater than 7°, which we interpreted as stars where angles in different bands are not compatible. We also obtained the value of the V-band polarization $p_V$.

<footnote n="3">³ The star HD 147933 from Anderson et al. (1996) is associated with the Ophiuchus cloud, but this star was eliminated from our sample in the selection process (Sect. 2.3).</footnote>

<footnote n="4">⁴ http://vizier.u-strasbg.fr</footnote>

<footnote n="2">² This ratio is called $R_{p/p}$ in Planck Collaboration Int. XXI (2015).</footnote>
for most of the AP07 stars by adopting the data directly from
the references (Covino et al. 1997; Arnal et al. 1993; Vrba et al.
1993; Whittet et al. 1992, 2001). The data from Anderson et al.
(1996) are in the form of (polarized) spectra rather than multi-
band photometry; for their lines of sight we took the polarization
at $\lambda = 545$ nm as the value for $p_V$. We did not use $p_V$ for the
stars in Martin et al. (1992), who fit data from multiple sources
and thus provide non-unique values for each band. Finally, we
obtained the extinction parameters for most of the stars: $A_V$ and
$E(B-V)$ for the AP07 stars, $E(B-V)$ for the Martin et al. (1992)
and Anderson et al. (1996) stars.

## 2.2. Emission: Planck and IRAS submm maps

Our submm data consist of the *Planck* 353 GHz (850 µm) maps for the *I*, *Q*, and *U* Stokes parameters, from which the polarized intensity $P_{353}$ and angle $\psi_{353}$ were obtained. We did not use any other frequencies because their signal-to-noise ratio (S/N) was lower. For selection purposes we also used the all-sky submm dust opacity maps created by Marc-Antoine Miville-Deschênes using *Planck* and *IRAS* data (Planck Collaboration XI 2014).

We used the second *Planck* public data release<fnref n="5" />, which consists of HEALPix all-sky maps of ten quantities: the Stokes parameters *I*, *Q*, and *U*, the number of hits, the variances of the Stokes parameters *I*, *QQ*, and *UU*, and the covariances *IQ*, *IU*, and *QU*. Maps are in NESTED ordering and Galactic coordinates; they have a pixelization $N_{side} = 2048$ for a total of $12 \cdot 2048^2 = 50331648$ pixels with 1'.7 side lengths, so that the beam of the instrument (*FWHM* ~ 5') is well sampled. The maps are in units of $K_{CMB}$ and are converted into MJy/sr with a conversion factor of 287.45 at 353 GHz (Planck Collaboration IX 2014). To obtain the value of *I*, *Q*, and *U* at the position of each star and increase the S/N, we employed the same technique as Planck Collaboration Int. XXI (2015); we averaged the values for the Stokes parameters on a Gaussian point-spread function (PSF) centered on the star coordinates and with a FWHM of 5', bringing the effective resolution to ~7'. In the case of *Q* and *U*, since we are working on a flat map recovered from a spherical map, we had to account for the fact that the direction of north changes from pixel to pixel. We did this by rotating the doublet (*Q*, *U*) until it was on the equator in the local reference frame.

With the value of the submillimeter *Q* and *U* for all the lines of sight, we calculated the polarized intensity in emission, $P_{353} = \sqrt{Q^2 + U^2}$, and the polarization angle $\psi_{353} = \frac{1}{2} \arctan (U, Q)$ using the HEALPix angle convention (where the relative signs of *Q* and *U* are inverted with respect with the IAU convention). Because it is a quadratic function of measures with finite noise, $P_{353}$ has a positive bias, and it was debiased with the conventional formula (Wardle & Kronberg 1974):
$$P_{\text{deb}} = \sqrt{P_{\text{bias}}^2 - \sigma_p^2}.$$ We did not apply the more recent debiasing method (e.g. Plaszczynski et al. 2014; Montier et al. 2015) because after the smoothing, the environments we studied have a high S/N, and therefore a low bias. There was no need to apply corrections for the cosmic microwave and infrared backgrounds, which are negligible at this wavelength and for our dataset.

For each star we also calculated the polarization angle dispersion function *S* (Planck Collaboration Int. XIX 2015), which is a tracer of disorder in polarization and therefore in the magnetic field orientation. The (*I*, *Q*, *U*) triplet was again smoothed

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Illustration of line-of-sight and beam effects when comparing extinction and emission, from Planck Collaboration Int. XXI (2015). The measured extinction is $E(B-V)$, the extinction obtained from conversion of submillimeter emission is $E(B-V)_S$.</figcaption></figure>

to increase S/N, using a Gaussian PSF with 5' FWHM and bringing the maps to a ~7' resolution. The maps thus obtained are oversampled (four pixels per beam), so we also degrade the pixelization of the *Q* and *U* maps to $N_{side} = 1024$ to approach the Nyquist criterion. The dispersion function *S* was then computed for the pixel containing the star, with a lag $\delta = 5'$.

## 2.3. Selection

Since this paper compares different phenomena (dust extinction and thermal emission) that are observed at different wavelengths (NUV to NIR vs. submm), we need to ensure that the comparison is meaningful (see e.g. Planck Collaboration Int. XXI 2015). This means that first, we must be sure that we are observing the same type of grain at the two wavelengths and that second, the two wavelengths must probe the same volumes of ISM.

The first condition is met, in first approximation, as a consequence of dust physics: only large grains contribute to polarization, and this is true both in extinction and in emission. Large grains are also the main contributor to submm emission as well as visual and IR extinction (small grains are important in the UV). It should be kept in mind, however, that polarization and overall extinction do not necessarily trace the same grains: a population of large grains that are spherical or unaligned would contribute to extinction and emission, but not to polarization.

The second condition is not trivial. Extinction measured on background stars, as is the case of the data described in Sect. 2.1, only probes the matter in front of the star itself; emission has no such limitation, especially in the submm, where the ISM is optically thin (see Fig. 1). In presence of a background to the star, the total intensity *I* measures systematically more dust than what is observed in the optical. The effect of background on polarization is more complex, and we detail it in Sect. 2.4. To compare extinction and emission, we therefore need to discard those lines of sight that have significant dust emission from behind the star. As shown in Planck Collaboration Int. XXI (2015), this selection can be based on three criteria:

- Galactic latitude: all stars close to the Galactic plane are very likely to have significant background; we therefore only keep stars with Galactic latitude $|b| \ge 2^\circ$. This forces us to exclude some well-studied clouds, such as the Coalsack nebula (AP07), which is located on the Galactic plane.

- Polarization angle: dust polarization in extinction should be orthogonal to that in emission. In lines of sight where the angles $\psi_{ext}$ and $\psi_{353}$ are not orthogonal, extinction and emission do not come from the same dust; we exclude these lines of

<footnote n="5">⁵ http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/</footnote>
sight with a tolerance of $3\sigma$ or $10^\circ$, whichever is smaller. For stars in Anderson et al. (1996), whose $\psi_{\text{ext}}$ are given without uncertainties, we assume $\sigma(\psi_{\text{ext}}) = 0$. Since Planck angle uncertainties are usually larger than V-band angle uncertainties, this should not make a large difference.

- Column density: the dust submm optical depth $\tau_{353}$ can be converted into an expected $A_V$ or $E(B-V)$ and compared to the actual extinction measured; lines of sight where the $\tau_{353}$-derived extinction shows an excess have significant background. For the conversion we used the empirical factor $E(B-V)/\tau_{353} = 1.49 \times 10^4$ that was obtained by Planck Collaboration XI (2014) for the diffuse ISM. In molecular clouds, however, $\tau_{353}/E(B-V)$ is known to increase by a factor ~2–3 compared to the diffuse ISM (e.g. Stepnik et al. 2003; Ysard et al. 2013; Planck Collaboration XI 2014). We therefore decided to relax this condition and we kept all lines of sight where the $\tau_{353}$-derived $E(B-V)$ is lower than three times the measured value.

While none of these selection procedures intrinsically exclude all of the contaminated sightlines, combining them and using them together with the selection made by AP07 gives more reliable results.

We improved the data quality by another selection. We only kept those stars for which we had an S/N greater than 3 in $P_{353}$ and greater than 5 in $\lambda_{\text{max}}$. For a few of the stars in AP07 the quality of the Serkowski fit was low and the Serkowski parameters were not an adequate representation of the polarization curve; we recovered the observational data from Covino et al. (1997), Arnal et al. (1993), Vrba et al. (1993), Whittet et al. (1992, 2001) and excluded those stars that, at a visual inspection, do not follow Serkowski. Finally, we excluded those stars that according to Anderson et al. (1996) are likely to have intrinsic polarization. The combined selection left us with the values of $\lambda_{\text{max}}$, $p_{\text{max}}$, $E(B-V)$, $I_{353}$, and $P_{353}$ for 70 lines of sight, 56 of which also have information on $A_V$ and $p_V$.

## 2.4. Line-of-sight and beam depolarization

The magnetic field in the ISM has a non-negligible disordered (or “meandering”) component that introduces a confounding variable called “depolarization”. When there is confusion between polarized sources with different orientation angles, it is possible for the orthogonal components of the polarizations to cancel each other out, so that the overall polarization that is observed may be lower than that of each source taken separately. Depolarization may occur if the interstellar magnetic field changes orientation along a line of sight (line-of-sight depolarization), or if an instrument has a finite observational beam and the magnetic field changes orientation on scales smaller than said beam (beam depolarization). In most polarization studies, the two effects are subsumed under the name of “beam depolarization” or simply “depolarization”. However, since the two types of depolarization have different effects on the extinction and the emission, we treat them separately here.

The line-of-sight depolarization, at first approximation, has the same effect on extinction and emission if extinction and emission probe the same ISM. Complications arise if there is significant emission from the background to the star: if the magnetic field orientation is very different in the foreground and in the background, depolarization in emission may be very different from that in extinction. This would give unreliable measurements of the ratio $P_{353}/p_V$, for example. The selection described

in Sect. 2.3, if effective, should ensure that the line-of-sight depolarization affects extinction and emission in the same way. We remark that having a uniform magnetic field orientation on the line of sight is not equivalent to having no background emission: since $P_{353}$ is additive, in this case we would observe an excess of polarization in emission as compared to extinction, and overestimate $P_{353}/p_V$.

The beam depolarization affects observations that have finite beam size. This is usually the case at FIR and submm wavelengths: the Planck beams measure $5'$ or more. Extinction observations on stars, on the other hand, are point-like and suffer no beam effects, so that beam depolarization affects only polarized emission. Unlike line-of-sight depolarization, beam effects are unaffected by our selection criteria. However, the amount of beam depolarization can be estimated from observational data, such as the function $S$ (Sect. 2.2), which measures field disorder.

## 2.5. Final observables

Our observational data that are to be compared to a dust model are plotted in Fig. 2. The top panel shows the normalized polarization in extinction, $p_V/\tau_V$, as a function of $\lambda_{\text{max}}$. The two quantities have a clear negative correlation; we also see that the values for polarization are very widely scattered, and their upper limit marks an envelope, as is typical of polarized observations, whose shape may be partly determined by line-of-sight depolarization (e.g. Planck Collaboration Int. XX 2015). A very similar behavior can be seen in the submm polarization fraction $P_{353}/I_{353}$ as a function of $\lambda_{\text{max}}$, shown in the central panel. The bottom panel of Fig. 2 shows a different observable: $P_{353}/p_{\text{max}}$. This quantity, like the ratio $P_{353}/p_V$ used in Planck Collaboration Int. XXI (2015), is meant to trace the optical properties of grains by normalizing out the effects of alignment and magnetic field orientation that affect both emission and extinction. We decided to use $p_{\text{max}}$ in the construction of this ratio, as opposed to $p_V$, to avoid introducing spurious correlations: many of our stars have high values of $\lambda_{\text{max}}$ (>0.6 µm), and $p_V$ and $\lambda_{\text{max}}$ are going to be negatively correlated in that range. The bottom panel of Fig. 2 shows that even if the scatter in $P_{353}/p_{\text{max}}$ is quite large, it is still small compared to the scatter observed in $p_V/\tau_V$ and $P_{353}/I_{353}$; the dependence of $P_{353}/p_{\text{max}}$ on $\lambda_{\text{max}}$ is also much less pronounced. This is consistent with our expectations that this quantity be nearly independent of alignment and magnetic field orientation.

# 3. Model: DUSTEM with polarization

The dust model we use should ideally have the three following characteristics: it should predict polarization in both extinction and emission, it should be compatible with the latest results from the Planck mission (especially the $P_{353}/p_{\text{max}}$ ratio, which is underpredicted by pre-Planck models), and it should allow to modify the dust properties to simulate dust evolution. Unfortunately, while models that allow detailed dust evolution exist, to our knowledge, they either do not predict polarization (e.g. Jones et al. 2013) or are calibrated on extinction alone and cannot be expected to reproduce the correct $P_{353}/p_{\text{max}}$ ratio (e.g. Hirashita & Voshchinnikov 2014). We decided to use a model optimized for fitting the latest Planck data instead, at the cost of a simplified treatment of dust evolution where only grain size is accounted for.

We adopted the dust model recently developed by Guillet et al. (2016), called “Model A”, which is a modified version of the model by Compiègne et al. (2011). The computation
<figure><img src="image_3.png" /><figcaption>Fig. 2. Observational data to reproduce. Different colors indicate different clouds, different symbols indicate different references (circles: AP07; triangles: Anderson et al. 1996, squares: Martin et al. 1992). Top: the “classical” $p_V/\tau$ vs. $\lambda_{max}$ plot. This panel does not contain the stars without $A_V$ and $p_V$ measurements. Center: comparison of $\lambda_{max}$ in extinction with the polarization $P_{353}/I_{353}$ in emission. Bottom: the ratio of polarization $P_{353}/p_{max}$ as a function of $\lambda_{max}$.</figcaption></figure>

is performed with the DustEM Fortran numerical tool<fnref n="6" /> and its IDL wrapper<fnref n="7" />. The model populations and size distributions are chosen to minimize the number of free parameters; the parameters themselves are calculated by fitting the observables typical

<figure><img src="image_4.png" /><figcaption>Fig. 3. Size distribution for the model we used (model A from Guillet et al. 2016).</figcaption></figure>

of the low-latitude diffuse ISM ($|b| < 30^{\circ}$): the extinction curve, the polarization in extinction up to 4 µm, and the spectral energy distribution (SED) in emission, polarized and non, updated with Planck results. The model therefore reproduces the average observations for the diffuse ISM, including ratios such as $I_{353}/A_V$ and $P_{353}/p_V$.

The model includes three grain types (see Table 2 and Fig. 3): a population of neutral polycyclic aromatic hydrocarbons (PAHs) with a lognormal size distribution, plus two populations of large grains (amorphous carbon and silicates) that are distributed as power laws: $dn/da = a^{\alpha}$. We are mainly interested in observables where the large grain contribution is dominant, therefore the model, unlike the model of Compègne et al. (2011), has no separate population for very small carbonaceous grains: the very small grains are included in the amorphous carbon population, which is why the power law for carbons is weighted toward small sizes. Large grains are prolate spheroids with an axial ratio of 3 (oblate grains of the same axial ratio cannot reproduce the high $P_{353}/I_{353}$ observed by Planck: Guillet et al. 2016). The neutral PAHs and the amorphous carbon grains have the same compositions as their counterparts in Compègne et al. (2011); the silicate grains have the same composition as in Weingartner & Draine (2001), with added porosity: 20% of their volume consists of vacuum inclusions. The porosity of the silicate grains is essential in increasing their $P_{353}/p_V$ ratio to the value observed by Planck.

In the model silicate grains are aligned according to the phenomenological alignment function provided by DustEM:

$$f(a) = \frac{1}{2} f_{\max} \left( 1 + \tanh \left( \frac{\ln(a/a_{\text{alig}})}{p_{\text{stiff}}} \right) \right), \quad (1)$$

where $f(a)$ is the fraction of grains aligned as a function of the equivalent radius<fnref n="8" /> $a$, $f_{\max}$ is the maximum alignment efficiency, $a_{\text{alig}}$ is the size threshold for grain alignment and the parameter $p_{\text{stiff}}$ denotes the width of the transition. This alignment function is designed to increase monotonically with size, since small grains are generally unaligned. This parametric function is not designed to test any particular alignment process, and while its shape resembles the typical result of the radiative torque model (Lazarian & Hoang 2007), it is also compatible with magnetic

<footnote n="6">⁶ https://www.ias.u-psud.fr/DUSTEM/</footnote>

<footnote n="7">⁷ http://dustemwrap.irap.omp.eu/</footnote>

<footnote n="8">⁸ For non-spherical grains, the equivalent radius is the radius of a sphere of corresponding volume.</footnote>
<table><caption>**Table 2.** Standard version of the dust model (model A from Guillet et al. 2016).</caption><thead><tr><th>Population</th><th>Mass (per H)</th><th>Depletion (ppm)</th><th>a<sub>max</sub> (μm)</th><th>α</th><th>a<sub>alig</sub> (μm)</th><th>p<sub>stiff</sub></th><th>f<sub>max</sub></th></tr></thead><tbody><tr><td>PAHs</td><td>7.10 × 10<sup>-4</sup></td><td>59</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Carbon BGs</td><td>1.32 × 10<sup>-3</sup></td><td>110</td><td>0.07</td><td>-4.14</td><td>—</td><td>—</td><td>—</td></tr><tr><td>Silicate BGs</td><td>6.52 × 10<sup>-3</sup></td><td>37.9</td><td>0.52</td><td>-3.32</td><td>0.108</td><td>0.27</td><td>1.00</td></tr></tbody></table>

Notes. “Depletion” in the case of silicates refers to Si, Mg and Fe. We used the Mathis et al. (1983) spectrum for the solar neighborhood as interstellar radiation field.

<figure><img src="image_5.png" /><figcaption>**Table 3.** Effective wavelengths used to interpolate the model polarization curve and produce synthetic *UBVRIJH* observations for the Serkowski fit.</figcaption></figure>

<table><thead><tr><th>Band</th><th>Effective λ (μm)</th></tr></thead><tbody><tr><td>U</td><td>0.36</td></tr><tr><td>B</td><td>0.44</td></tr><tr><td>V</td><td>0.55</td></tr><tr><td>R</td><td>0.65</td></tr><tr><td>I</td><td>0.80</td></tr><tr><td>J</td><td>1.25</td></tr><tr><td>H</td><td>1.60</td></tr></tbody></table>

Notes. Adapted from Whittet et al. (1992, 2001), Bessell & Murphy (2012).

<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Top: grain alignment function for different values of $a_{alig}$. Bottom: corresponding polarization curves in extinction, normalized to $N_H = 10^{21} \text{ cm}^{-2}$. Dots indicate the value of $a_{alig}$ (top) and that of $\lambda_{max}$ from free-$K$ fits (bottom).</figcaption></figure>

alignment for grains with superparamagnetic inclusions (Mathis 1986; Voshchinnikov et al. 2016). The variation of the alignment function for varying $a_{alig}$ and its effects on the polarization curve in extinction, are shown in Fig. 4.

In addition to causing polarization, grain alignment affects dust extinction and emission as well. The resulting correction is very small and it is generally ignored in dust models; nonetheless, DustEM provides the option of including alignment effects in extinction and emission. We here chose to ignore these effects to allow a more direct comparison to the results of Guillet et al. (2016), where minor alignment effects were ignored as well.

### 3.1. Model results: fitting

DustEM provides the full extinction curve and dust emission for the model, including polarization; most of the parameters introduced in Sect. 2.5 have to be obtained from a fit to the DustEM output.

We took the extinction curve interpolated at 550 nm to be $A_V$ of the model, and the thermal emission at 353 GHz, plus a color correction to account for the spectral response of the corresponding Planck band, as the value $I_{353}$ of the SED. The same interpolation, operated on the polarized extinction and emission, gives us the model prediction for $p_V$ and $P_{353}$. The $p_{max}$, $\lambda_{max}$, and $K$ of the model were calculated by interpolating the model polarized extinction at the effective central wavelengths of the *UBVRIJH* photometric bands (Table 3) and fitting a Serkowski function to the synthetic observations thus obtained, keeping $K$ as a free parameter.

We also fit the model SED with a modified blackbody: $I_{\lambda} = B_{\lambda}(T) \cdot \tau_0 \cdot (\lambda/\lambda_0)^{-\beta}$. The fit was performed on the model emission interpolated and color-corrected at the wavelengths of the Planck HFI bands (350, 550, 850, 1380, and 2100 µm) and the IRAS 100 µm band. Since emission at those wavelengths is dominated by large grains, integrating the modified blackbody over wavelength provides the radiance $R$, or emitted power, of the large grain populations (see Planck Collaboration XI 2014).

### 3.2. Model variations

The model so far described was developed to fit the average observables in the diffuse ISM. Inside molecular clouds, however, evolution alters the properties of dust considerably: the main alteration is grain growth through accretion and coagulation (Boulanger et al. 1990; Stepinik et al. 2003; Köhler et al. 2012). Translucent lines of sight, such as those studied in this paper, typically probe regions at the onset of such an evolution: for
instance, Stepnik et al. (2003) and Ysard et al. (2013) found that coagulation takes place where $A_V$ is greater than 2 or 3. This raises the question of whether our sample can be explained by a dust model designed for the diffuse ISM. To answer this question, we studied how the model output is affected by grain alignment efficiency, magnetic field orientation, and grain size. By comparing these results to the observations we can determine whether the data can be explained by the variation of alignment and field orientation alone, using the same dust as in the diffuse ISM, or if dust growth is necessary, and to what extent. The modifications to our model are purely phenomenological and are not based on simulations of grain growth, dust alignment or magnetic field structure; however, they are useful for estimating the variations that physical models will have to reproduce.

### 3.2.1. Variations in alignment function

Loss of grain alignment inside molecular clouds is sometimes invoked to explain the decrease in polarization efficiency observed at high $A_V$ (e.g. Andersson et al. 2015, and references therein). This weakening of polarization is also in qualitative agreement with some alignment theories, e.g., the radiative torque model, which predicts that only the largest grains are aligned inside molecular clouds.

The alignment efficiency in DustEM is a function of the three parameters $a_{alig}$, $p_{stiff}$, and $f_{max}$, as shown in Eq. (1). We simulated different alignment efficiencies by running the model with different values of $a_{alig}$, keeping the same $p_{stiff}$ for simplicity. We did not change $f_{max}$ (equal to 1 in our model) because it has the exact same effects on $p_V$ and $P_{353}$, and it could have no effect on $P_{353}/p_V$; the parameter $a_{alig}$, on the other hand, affects $p_V$ and $P_{353}$ similarly but not identically, since the polarization cross-section has a different size-dependent behavior in the visible (where scattering is dominant) and in the submm.

Figure 4 shows the change in $f(a)$ as a function of $a_{alig}$ and the effect that this has on the polarization curve; higher $a_{alig}$ corresponds to lower $p_{max}$ because fewer grains are aligned, and to higher $\lambda_{max}$ because the average aligned grain is larger.

### 3.2.2. Variations in magnetic field orientation

An ordered magnetic field forming an angle $\gamma$ with the plane of the sky introduces a factor $\cos^2 \gamma$ in the polarized intensity, as shown in Sect. 1. Dust polarization models would greatly benefit from measures of the angle $\gamma$; unfortunately, dust only traces the magnetic field component parallel to the plane of the sky, which means that this information is not usually available. The angle $\gamma$ is therefore another variable parameter in our model: we ran the model for $\gamma = 0^\circ, 30^\circ, 45^\circ$, and $60^\circ$.

The Galactic magnetic field is the sum of an ordered component and a disordered, or meandering, component: this latter component causes the phenomena known as line-of-sight and beam depolarization, as explained in Sect. 2.4. Our model does not include a disordered magnetic field component and therefore it cannot predict depolarization; however, the polarization angle dispersion $S$ (Sect. 2.2) can be used as a measure of field disorder (see e.g. Planck Collaboration Int. XX 2015, which compare the observed $S$ with MHD simulations). Using this we were able to assess some of the effects of the disordered magnetic field, as shown in Sect. 4.1.

### 3.2.3. Variations in grain size distribution

Gas accretion on grain surfaces (e.g. Jones et al. 2013) and formation of aggregates are known to increase grain sizes inside molecular clouds. This growth is supported by theoretical studies (Köhler et al. 2012; Hirashita & Voshchinnikov 2014) and it is consistent with observed phenomena such as the flattening of extinction curves in dense environments (e.g. Fitzpatrick 1999; Weingartner & Draine 2001) and the coreshine observed in the NIR (Pagani et al. 2010).

As we described, no dust models currently treat dust evolution realistically while reproducing the emission-to-extinction polarization ratio revealed by Planck. We therefore opted for a model that is compatible with Planck data (model A from Guillet et al. 2016) at the cost of a simplified treatment of dust evolution. Specifically, we focused on a single aspect affected by dust evolution: the size distribution of grains. In our model, the size distribution for large grains is a power law defined by three parameters: the minimum and maximum grain sizes $a_{min}$ and $a_{max}$, and the power-law index $\alpha$.

Most of the mass in silicates is in the large grains (see Fig. 3), meaning that the size distribution is most sensitive to $a_{max}$; furthermore, none of the observables we use are in the UV where the contribution of small grains is important. Therefore, we decided to vary the $a_{max}$ of silicates between 350 nm and 1 µm (the standard value is ~500 nm) while keeping $\alpha$ and $a_{min}$ fixed. Although we chose to fix $\alpha$ mainly as a matter of convenience, we note that this is consistent with the model of grain growth by Hirashita & Voshchinnikov (2014), where the slope of the size distribution does not change much during evolution, and the largest variation is in the upper size cutoff.

The case of carbon grains is different, as their distribution is weighted toward small sizes: the mass available for large grains is now also dependent on $a_{min}$ and on the amount of PAHs, so that a realistic model becomes a necessity. For this reason we only varied the size distribution of silicates while leaving that of carbon grains fixed. While this choice does not give realistic results for the variation of extinction and emission with size distribution, it still allows us to predict the dust polarization, which in our model depends on silicates alone.

### 3.2.4. Multiparameter study: Monte Carlo simulation

The phenomena described in Sects. 3.2.1 to 3.2.3 are all expected to occur in molecular clouds, so that variations of a single model parameter at a time are not realistic, even if studying their effect can be instructive. We decided to use a Monte Carlo simulation to explore the effects of simultaneous variations of many parameters. As explained in the previous section, we can only vary the size distribution of silicates, which gives realistic results for polarization but not for unpolarized observables; therefore our Monte Carlo results can only be compared to polarization observables. The model was run one thousand times, and the values of $a_{alig}$ and $a_{max}$ for silicates were uniformly distributed within the ranges

$$50 \text{ nm} < a_{alig} < 300 \text{ nm}$$

$$350 \text{ nm} < a_{max} < 1000 \text{ nm}.$$

The operation was repeated for four values of $\gamma$ ($0^\circ, 30^\circ, 45^\circ$, and $60^\circ$), bringing the simulations to a total of 4000.

We found that not all combinations of $a_{alig}$, $a_{max}$ and $\gamma$ in the ranges chosen are realistic: the synthetic observables obtained for some such combinations have values that are never observed. We set out to find a realistic range of parameters by imposing that
model results have a range as close as possible to that of actual observations (e.g., that $0.35 < \lambda_{max} < 0.8 \mu m$ and $0.5 < K < 1.5$ as per Voshchinnikov & Hirashita 2014). Restricting $a_{alig}$ and $a_{max}$ to the following ranges eliminates most of the unrealistic values for $\lambda_{max}$ and K:

$$ 
\begin{array}{l}
75 \text{ nm} < a_{\text{alig}} < 150 \text{ nm} \\
350 \text{ nm} < a_{\text{max}} < 800 \text{ nm},
\end{array}
$$

while the same four values for $\gamma$ are kept. This selection left us with 844 Monte Carlo iterations.

# 4. Results

## 4.1. Alignment efficiency and magnetic field orientation

The effects of dust alignment efficiency and magnetic field orientation are shown in Fig. 5, which compares the model results with observational data. Dots represent the observed values of $p_V/\tau_V$ (top), $P_{353}/I_{353}$ (middle), and $P_{353}/p_{max}$ (bottom) as a function of $\lambda_{max}$. Curves represent the model; within each curve, $a_{alig}$ varies between 75 and 150 nm, and the four curves correspond to the four values of $\gamma$, 0°, 30°, 45°, and 60°.

The combination of variable alignment and magnetic field orientation can reproduce most of the observations in the case of $p_V/\tau_V$ and $P_{353}/I_{353}$, both the general trends and the dispersion. The curve for $\gamma = 0°$ coincides roughly with the highly polarized, low-$\lambda_{max}$ lines of sight in our sample. For higher values of $\gamma$ the polarization decreases, but the dispersion in $\lambda_{max}$ caused by the variation of $a_{alig}$ increases, pushing the maximum $\lambda_{max}$ to higher values: as a result, the model predicts that weakly polarized lines of sight that can have either small or large $\lambda_{max}$, which is indeed the trend found in the observational data. The relation between $\gamma$ and $\lambda_{max}$ described by Voshchinnikov et al. (2016) is evident in Fig. 5: although $\lambda_{max}$ is mainly affected by the alignment size threshold $a_{alig}$, the model curves with a higher $\gamma$ are clearly shifted to higher values of $\lambda_{max}$. The figure, however, reveals something more: the strength of the $\gamma$-$\lambda_{max}$ relation itself increases with $a_{alig}$. The leftmost tips of the model curves, corresponding to $a_{alig} = 75$ nm, all have very similar $\lambda_{max}$. In contrast, the rightmost tips, corresponding to $a_{alig} = 150$ nm, show wide differences in $\lambda_{max}$, comparable to the differences that are due to $a_{alig}$ itself. It should be noted that our model assumes a uniform magnetic field, and therefore it does not include line-of-sight or beam depolarization. If these effects were important, it might mean that the ordered component of the field is closer to the plane of the sky than our model predicts.

Figure 5 also shows that alignment and magnetic field orientation have very little effect on the model results for $P_{353}/p_{max}$, as was indeed expected: the different curves are close to each other. In fact, the dispersion observed in the observed $P_{353}/p_{max}$ is much larger than predicted by the model, which suggests that variations in the dust optical properties occur in translucent lines of sight. Again, one possible confounding factor in this interpretation is the depolarization caused by meandering of the magnetic field: we will now attempt to assess its effects.

While our model assumes a uniform magnetic field, the disorder in the field lines can be estimated from the angle dispersion function S (Sect. 2.2). The top panel of Fig. 6 shows an anticorrelation between S and the polarization fraction in emission, a well-known effect that is commonly attributed to line-of-sight depolarization caused by meandering of the field (e.g. Planck Collaboration Int. XX 2015). Beam depolarization is also a possible cause, but we can see that beam effects appear negligible in our sample: the bottom panel of Fig. 6 shows

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>Fig. 5. Observational data compared with our variable-alignment model ($75 \le a_{alig} \le 150$ nm). The color and symbol scheme are the same as in Fig. 2. The black square marks the position of the standard model. The top panel does not contain the stars without $A_V$ and $p_V$ measurements.</figcaption></figure>

that there is no clear influence of field meandering on the polarization ratio $P_{353}/p_{max}$. Beam depolarization, if present, should affect emission but not extinction, introducing an anticorrelation between S and $P_{353}/p_{max}$. The fact that we see no such correlation suggests that beam depolarization is negligible in our sample: this supports the idea that we are probing dust in relatively homogeneous regions, which is an important assumption in
<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 6.** Correlation between the observables we studied and the magnetic field meandering as measured by *S*. Top: polarization fraction $P_{353}/I_{353}$ as a function of $\lambda_{max}$. Bottom: polarization ratio $P_{353}/p_{max}$ as a function of $\lambda_{max}$. The color and symbol scheme are the same as in Fig. 2.</figcaption></figure>

extinction/emission comparison. The effect of line-of-sight depolarization on the polarization fraction is unfortunately impossible to assess without more advanced modeling, but even if present, it should have little effect on $P_{353}/p_{max}$ following our selection (Sects. 2.3 and 2.4).

## 4.2. Grain size distribution

As explained in Sect. 3.2.3, we ignored variations of size distribution in carbonaceous grains, which affect $I_{353}$ and $A_V$, but not $P_{353}$ and $p_{max}$. Because of this, our results on the effects of grain growth on observations are unrealistic for $P_{353}/I_{353}$ and $p_V/\tau_V$, but realistic for $P_{353}/p_{max}$, and we show only this last observable.

Figure 7 compares the observations and the model results for $350 \le a_{max} \le 10^3$ nm (solid black line). The effect of a variable $a_{max}$ on $\lambda_{max}$ is small compared to the effect of grain alignment; on the other hand, $a_{max}$ has a strong effect on $P_{353}/p_{max}$ and is a plausible contributor to the large dispersion observed in this quantity. For comparison, the picture also shows the model results for fixed $a_{max}$ and $\alpha$ varying of $\pm 0.5$ around the standard value, that is, $-3.82 \le \alpha \le -2.82$ (dashed line). The variations in $\alpha$ have a modest effect on $\lambda_{max}$ comparable to that of $a_{max}$;

<figure><img src="image_12.png" /><figcaption>**Fig. 7.** Observational polarization ratio $P_{353}/p_{max}$ compared with a variable size distribution model. The solid line shows the evolution of the model for $350 \le a_{max} \le 10^3$ nm (**left to right**). The dashed line shows the effect of variable power-law index: $-3.82 \le \alpha \le -2.82$ (**left to right**). The black square marks the position of the standard model. The color and symbol scheme are the same as in Fig. 2.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Fig. 8.** Observational polarization ratio $P_{353}/p_{max}$ compared with Monte Carlo iterations of the model (grey). The black square marks the position of the standard model. The color and symbol scheme are the same as in Fig. 2.</figcaption></figure>

however, unlike $a_{max}$, the parameter $\alpha$ has nearly no influence on the $P_{353}/p_{max}$ ratio.

## 4.3. Multiparameter Monte Carlo

We mentioned above that, realistically, all the model parameters we use will simultaneously vary inside a molecular cloud, and we chose to represent this with a Monte Carlo simulation (Sect. 3.2.4). As in the previous section, our modeling does not account for the size distribution of carbon grains, and we only show polarization results.

Figure 8 compares the Monte Carlo results to the observational data. The Monte Carlo simulation was run for $\gamma = 0^\circ, 30^\circ, 45^\circ$, and $60^\circ$; the other model parameters are uniformly distributed in the regions $75 \text{ nm} \le a_{\text{alig}} \le 150 \text{ nm}$ and $350 \le a_{max} \le 800 \text{ nm}$. The combined variation of dust alignment, field orientation, and grain size allows the model to reproduce most of the observations. However, some lines of sight, spanning the full
<figure><img src="image_14.png" /><figcaption>**Fig. 9.** *P*<sub>353</sub>/*p*<sub>max</sub> ratio as a function of the modified blackbody temperature (top panel) and as a function of the normalized radiance *R*/*A*<sub>V</sub> (bottom panel). The bottom panel does not include stars without *A*<sub>V</sub> and *p*<sub>V</sub> measurements. The grey curves show the model for varying intensity of the radiation field (0.1 ≤ *G*<sub>0</sub> ≤ 3). The grey square marks the position of the standard model. The color and symbol scheme for the observations are the same as in Fig. 2.</figcaption></figure>

range of observed $\lambda_{\max}$, have a lower $P_{353}/p_{\max}$ than the model can reproduce. This may be the result of variations in the dust polarization properties that are due to factors other than size (e.g., grain shape, porosity, or chemical composition), which are not considered in our model.

**5. Discussion**

We have shown that the grain alignment and the orientation of the magnetic field both have important effects on our observables: a combination of variable $a_{\text{alig}}$ and variable $\gamma$ reproduces the general trends and the dispersions for $p_{\text{V}}/\tau_{\text{V}}$ and $P_{353}/I_{353}$ as a function of $\lambda_{\max}$; in fact, we find the familiar envelope in the distribution of points (Fig. 5). There are some stars, however, mainly in the Musca cloud, with a higher polarization than the model can reproduce. This is because the model was made to reproduce a $p_{\text{V}}/\tau_{\text{V}}$ of ~3%, the usually accepted value for the diffuse ISM (e.g. Serkowski et al. 1975; Andersson et al. 2015), while the outlier stars in Fig. 5 have $p_{\text{V}}/\tau_{\text{V}} \sim 4\%$. We should point out, however, that the measures in question have relatively large error bars and most stars are within ~1$\sigma$ above the envelope.

Figure 8 also shows that a combination of variable alignment, magnetic field orientation, and grain growth reproduces the general trend and most of the data scatter in the $P_{353}/p_{\max}$ vs. $\lambda_{\max}$ relation. Again, some data are beyond the model range; in this case it is the stars with low $P_{353}/p_{\max}$. These lines of sight may belong to regions with different dust properties: while our model reproduces the high $P_{353}/p_{\max}$ ratio from Planck Collaboration Int. XXI (2015), most models give lower values. Another possibility is that the low $P_{353}/p_{\max}$ is an effect of beam depolarization, which would lower the observed value of $P_{353}$ without affecting $p_{\max}$; however, we showed in Sect. 4.1 that beam depolarization does not seem important in our sample. A low $P_{353}/p_{\max}$ could also be due to low dust heating, because $P_{353}$ is proportional to dust emission: our model results are calculated for the typical interstellar radiation field intensity<fnref n="9" /> in the diffuse ISM, $G_0 = 1$, but we would expect the radiation field to be less intense in clouds. However, in the long-wavelength range of the Planck function, which is the case of Planck observations of dust, emission tends to depend linearly on dust temperature, meaning that the effect of heating on dust SED may be weak.

A proper estimation of heating effects would be no trivial task, as there is no univocal relation between the observed $A_{\text{V}}$ and the extinction actually experienced by dust (see AP07 for a discussion), but a preliminary analysis is shown in Fig. 9. The gray curves in the figure represent the model results for $0.1 \le G_0 \le 3$. The top panel shows no visible correlation between $P_{353}/p_{\max}$ and dust temperature, suggesting that $G_0$ only has modest effects (it should be kept in mind, though, that temperature is also influenced by dust evolution: e.g., Stepnik et al. 2003; Planck Collaboration Int. XVII 2014). The same conclusion seems supported by the bottom panel: this shows that $P_{353}/p_{\max}$ is not correlated with the normalized radiance $R/A_{\text{V}}$. Since the radiance $R$ is the bolometric emission of large dust grains (Sect. 3.1), $R/A_{\text{V}}$ is a measure of the power emitted (and therefore absorbed) by large grains, averaged on the line of sight<fnref n="10" />.

**6. Conclusions and perspectives**

Dust polarization depends on the efficiency of grain alignment, the orientation and meandering of the magnetic field, and the optical properties of the dust itself. Many studies have attempted to explain polarization in molecular clouds in terms of one of these factors (see, e.g., AP07 for a study focused on alignment and Planck Collaboration Int. XX (2015) for one focused on magnetic field meandering); nonetheless, it is likely that all factors are at play at once.

In this paper we used dust model A from Guillet et al. (2016) and varied its alignment efficiency, grain size, and magnetic field orientation to reproduce the diverse conditions that one may find in molecular clouds. The results were compared to extinction and emission data from both bibliographic sources and the Planck survey, with particular attention to the polarization observables $p_{\text{V}}/\tau_{\text{V}}$, $P_{353}/I_{353}$, and $P_{353}/p_{\max}$ as a function of $\lambda_{\max}$.

We find that none of the model parameters employed can explain the full set of observations on its own. Monte Carlo

<footnote n="9">⁹ We used the default interstellar radiation field in DustEM, which is the Mathis et al. (1983) SED for the solar neighborhood multiplied by the dimensionless intensity parameter $G_0$.</footnote>

<footnote n="10">¹⁰ Stars 10 and 24 of Ophiuchus (according to the nomenclature of Vrba et al. 1993) are not shown in the bottom panel of Fig. 9, because their large $R/A_{\text{V}}$ place them outside the plot area. The ($R/A_{\text{V}}$, $P_{353}/p_{\max}$) values for these outliers are 2.15 and 6.45 for star 10 and 3.24 and 6.98 for star 24 in the units shown on the figure axes.</footnote>
simulations show that most of the data can be reproduced by letting $a_{\text{alig}}$ vary between 75 and 150 nm, $a_{\text{max}}$ vary between 350 nm and 800 nm, and $\gamma$ vary between 0° and 60°. This means that any studies of polarization in molecular clouds need to take into account all these aspects to explain the full range of the data. In particular, the ratio $P_{353}/p_{\text{max}}$ is very useful in reducing the contributions of alignment and magnetic field orientation, and highlighting variations in dust properties, and especially size distribution. Within the context of our model, the variations observed in $P_{353}/p_{\text{max}}$ can be partly explained by varying the maximum grain size $a_{\text{max}}$, while the power-law index $\alpha$ has little effect.

Nonetheless, some of the observations fall beyond the range of the model results, most notably some lines of sight with very low values of $P_{353}/p_{\text{max}}$ that are found over the full range of observed $\lambda_{\text{max}}$. This probably indicates variations in dust properties other than size distribution, such as shape, structure, or chemical composition. Dust evolution that is not related to size may also influence our estimate of magnetic field orientation: for instance, lines of sight with a low $P_{353}/I_{353}$, which our model attributes to magnetic field lines nearly orthogonal to the plane of the sky (large $\gamma$), may be instead explained by dust with a low polarization cross-section. The lines of sight with low $P_{353}/p_{\text{max}}$ can in principle be explained without a need for dust evolution: a dim radiation field due to extinction, or the beam depolarization due to a disordered magnetic field, would lower the value of polarization in emission without affecting extinction. However the analysis we conducted in the previous sections suggests that neither dust heating nor beam depolarization have much effect on our sample.

It is possible that the width of the magnetic field orientation range found by our model (0° to 60°) is overestimated. In addition to the aforementioned influence of dust properties, field meandering, which is absent in our model, could introduce line-of-sight depolarization, the effects of which are degenerate with increasing angle $\gamma$ of (the ordered component of) the magnetic field. Including field meandering in an ISM model would make for a very interesting follow-up study, but a polarized radiation transfer code is needed for this. It would also be useful to independently determine $\gamma$ in future research, for example, using MHD simulations or where available, measures of the line-of-sight magnetic field such as the Zeeman polarization.

Not all potentially interesting cases were considered in this paper: this is a first application of this technique to the study of polarization in molecular clouds. The full implications of this technique will become clearer with more detailed modeling and more observational data. A continuation of this work would benefit from extending the dataset to NIR extinction and polarization. Observations in the NIR can probe denser lines of sight than those in the visible; furthermore, increasing the wavelength of observation means that we are closer to the Rayleigh limit, so that the observables are better trackers of the overall mass of aligned grains and less dependent on the details of alignment and size distribution. This makes NIR observations an interesting complement to observations in the optical. Finally, different

types of observations could improve our constraints on the model variables: maps of molecular lines and elemental depletion could be useful in constraining grain growth processes such as accretion and coagulation.

**Acknowledgements.** We would like to thank S. N. Shore for his stimulating discussion and insightful remarks on a variety of subjects, from ISM physics to data analysis to writing, during the editing of this article. We are also indebted to N. V. Voshchinnikov for his helpful comments.

## References

Anderson, C. M., Weitenbeck, A. J., Code, A. D., et al. 1996, AJ, 112, 2726

Andersson, B.-G., & Potter, S. B. 2007, ApJ, 665, 369

Andersson, B.-G., Lazarian, A., & Vaillancourt, J. E. 2015, ARA&A, 53, 501

Arnal, E. M., Morras, R., & Rizzo, J. R. 1993, MNRAS, 265, 1

Bessell, M., & Murphy, S. 2012, PASP, 124, 140

Boulanger, F., Falgarone, E., Puget, J. L., & Helou, G. 1990, ApJ, 364, 136

Compiègne, M., Verstraete, L., Jones, A., et al. 2011, A&A, 525, A103

Corradi, W. J. B., Franco, G. A. P., & Knude, J. 2004, MNRAS, 347, 1065

Covino, E., Palazzi, E., Penprase, B. E., Schwarz, H. E., & Terranegra, L. 1997, A&AS, 122, 95

Fitzpatrick, E. L. 1999, PASP, 111, 63

Greenberg, J. M. 1968, in Nebulae and Interstellar Matter (the University of Chicago Press), 221

Guillet, V., Fanciullo, L., Verstraete, L., et al. 2016, A&A, submitted

Hall, J. S. 1949, Science, 109, 166

Hildebrand, R. H. 1988, QJRAS, 29, 327

Hiltner, W. A. 1949, Science, 109, 165

Hirashita, H., & Voshchinnikov, N. V. 2014, MNRAS, 437, 1636

Jones, A. P., Fanciullo, L., Köhler, M., et al. 2013, A&A, 558, A62

Kim, S.-H., & Martin, P. G. 1995, ApJ, 444, 293

Köhler, M., Stepnik, B., Jones, A. P., et al. 2012, A&A, 548, A61

Lazarian, A., & Hoang, T. 2007, MNRAS, 378, 910

Lee, H. M., & Draine, B. T. 1985, ApJ, 290, 211

Loinard, L. 2013, in Advancing the Physics of Cosmic Distances, ed. R. de Grijs, IAU Symp., 289, 36

Martin, P. G., Adamson, A. J., Whittet, D. C. B., et al. 1992, ApJ, 392, 691

Mathis, J. S. 1986, ApJ, 308, 281

Mathis, J. S., Mezger, P. G., & Panagia, N. 1983, A&A, 128, 212

Montier, L., Plaszczyński, S., Levrier, F., et al. 2015, A&A, 574, A136

Neuhäuser, R., & Forbrich, J. 2008, in Handbook of Star Forming Regions, Vol. II (ASP Monographs; San Francisco, CA: ASP), 735

Pagani, L., Steinacker, J., Bacmann, A., Stutz, A., & Henning, T. 2010, Science, 329, 1622

Planck Collaboration IX. 2014, A&A, 571, A9

Planck Collaboration XI. 2014, A&A, 571, A11

Planck Collaboration Int. XVII. 2014, A&A, 566, A55

Planck Collaboration Int. XIX. 2015, A&A, 576, A104

Planck Collaboration Int. XX. 2015, A&A, 576, A105

Planck Collaboration Int. XXI. 2015, A&A, 576, A106

Plaszczyński, S., Montier, L., Levrier, F., & Tristram, M. 2014, MNRAS, 439, 4048

Serkowski, K., Mathewson, D. S., & Ford, V. L. 1975, ApJ, 196, 261

Stepnik, B., Abergel, A., Bernard, J.-P., et al. 2003, A&A, 398, 551

Voshchinnikov, N. V., & Hirashita, H. 2014, MNRAS, 445, 301

Voshchinnikov, N. V., Il'in, V. B., & Das, H. K. 2016, MNRAS, 462, 2343

Vrba, F. J., Coyne, G. V., & Tapia, S. 1993, AJ, 105, 1010

Wardle, J. F. C., & Kronberg, P. P. 1974, ApJ, 194, 249

Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296

Whittet, D. C. B., Martin, P. G., Hough, J. H., et al. 1992, ApJ, 386, 562

Whittet, D. C. B., Gerakines, P. A., Hough, J. H., & Shenoy, S. S. 2001, ApJ, 547, 872

Ysard, N., Abergel, A., Ristorcelli, I., et al. 2013, A&A, 559, A133