MID-INFRARED POLARIMETRY: INSIGHTS INTO MAGNETIC FIELDS AND DUST GRAIN PROPERTIES IN YOUNG STELLAR OBJECTS AND PROTOPLANETARY DISKS
By HAN ZHANG
A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY UNIVERSITY OF FLORIDA
2018 ⃝c 2018 Han Zhang When I look at your heavens, the work of your fingers, the moon and the stars, which you have set in place, what is mankind that you are mindful of them, human beings that you care for them? – Psalm 8:3-4 ACKNOWLEDGMENTS Firstly, I would like to thank my advisor, Charles M. Telesco, for his patience, guidance, and generous financial supports for the past five years. It is a great privilege for me to work with Charlie and be part of the CanariCam science team.
I want to thank all my collaborators, without your contributions, I can not finish the work presented in this thesis. I am particularly thankful to Dan Li, Eric Pantin, Aigen Li, Christopher
M. Wright, Peter Barnes, and Naib´ıMari˜nas.Thank you for your insightful discussions and comments.
I thank my fellow graduate students here at University of Florida. My great classmates
Rebecca, Nolan, Chen and Pekki, I remember the days we worked on the homework and projects together. I want to thank Jingzhe, my first-year roommate at UF, and Xiao, my driving coach. I want to thank Emily, Wenli, Chutipong, Shuo, Nahathai, Hanna, Tahlia,
Yinan, Amanda, Rachel, Alan, Krittapas, Ben Wu, and Billy. Thank you all for your company and love. Go Gators! Finally, a special thanks to my parents for always believing in me and supporting me. I am grateful for all my church friends both in Gainesville and China. Thank you for your prayers, comfort and encouragement. We love because He first loved.
4 TABLE OF CONTENTS page
ACKNOWLEDGMENTS ...... 4
LIST OF TABLES ...... 7
LIST OF FIGURES ...... 8 ABSTRACT ...... 10
CHAPTER
1 INTRODUCTION ...... 12
1.1 Magnetic Fields in Star Formation ...... 12 1.2 Protoplaneatry Disks ...... 13 1.3 Polarimetry: An Important Technique to Study Magnetic Fields ...... 15 1.4 Challenges: Scattering, Radiative Grain Alignment ...... 17 1.5 Thesis Outline ...... 19
2 OBSERVATIONS AND DATA REDUCTION ...... 23
2.1 CanariCam ...... 23 2.2 Stokes Parameters ...... 23 2.3 Data Reduction ...... 24 2.3.1 Polarimetric Imaging Data Reduction ...... 25 2.3.2 Spectropolarimetry Data Reduction ...... 27
3 THE MID-INFRARED POLARIZATION OF THE HERBIG AE STAR WL 16: AN INTERSTELLAR ORIGIN? ...... 30
3.1 Introduction ...... 30 3.2 Observation ...... 32 3.3 Orign of the Polarization ...... 34 3.3.1 Extinction toward WL 16 ...... 34 3.3.2 Decomposition of Absorptive and Emissive Polarization ...... 40 3.3.3 Relationship of WL 16 and Regional Magnetic Fields ...... 43 3.3.4 Intrinsic Polarization ...... 44 3.4 Discussions ...... 45 3.4.1 PAHs in WL 16 ...... 45 3.4.2 Disk Morphology ...... 46 3.5 Summary ...... 48
4 DETECTION OF POLARIZED INFRARED EMISSION BY POLYCYCLIC AROMATIC HYDROCARBONS IN THE MWC 1080 NEBULA ...... 50
4.1 Introduction ...... 50 4.2 Observation and Data Reduction ...... 51
5 4.3 Results ...... 53 4.4 Discussion ...... 54 4.4.1 Numerical Calculations using SD09 ...... 55 4.4.2 Alignment with Magnetic Fields ...... 56 4.4.3 Relationship between Polarization Angles and the Ambient Magnetic Field ...... 58 4.4.4 Marginally Detected 10.3 µm Polarization Feature ...... 59 4.5 Summary ...... 60 5 MODELING POLARIZATION OF YOUNG STELLAR OBJECTS AND PROTOPLANETARY DISKS AT MID-IR ...... 64
5.1 Introduction ...... 64 5.2 Theoretical Understanding of Dust Polarization ...... 65 5.2.1 Polarization from Dichnoic Emission and Absorption ...... 66 5.2.2 Dust Scattering ...... 68 5.3 Model Description ...... 69 5.3.1 Magnetic Fields Setup ...... 69 5.3.2 Fiducial Model: Spherical Power-law Envelope ...... 70 5.3.3 Radiation Transfer with RADMC-3D ...... 70 5.3.4 Disk Model ...... 72 5.3.5 Results ...... 73 5.3.6 Example: AB Aur ...... 75 5.3.7 Discussion ...... 76 5.4 Summary ...... 77
6 UNDERSTANDING THE MAGNETIC FIELDS IN W51 IRS2 USING MID-IR POLARIMETRY 90 6.1 Introduction ...... 90 6.2 Observation and Data Reduction ...... 92 6.3 Discussion ...... 93 6.3.1 Polarization Components – The Aitken Method ...... 93 6.3.2 Polarization Results of W51 IRS2 ...... 94 6.3.3 Magnetic Field Structure ...... 96 6.3.4 Gas Emission from VLA ...... 97 6.3.5 Magnetically Driven Gas Flow? ...... 98 6.4 Summary ...... 99
7 CONCLUSIONS ...... 106
7.1 Future Directions ...... 108 APPENDIX: AITKEN’S METHOD ...... 109
REFERENCES ...... 111
BIOGRAPHICAL SKETCH ...... 119
6 LIST OF TABLES Table page
3-1 Basic Properties of WL 16 ...... 32
3-2 Observing Log ...... 34
3-3 Polarization Measurements of WL 16 ...... 43 3-4 Extinction and Polarization of WL 16 and Elias 29 ...... 44
4-1 Observing Log ...... 60
4-2 Different models and polarization at 11.3 µm ...... 61
5-1 Model Parameters ...... 78 6-1 Observation Log ...... 100
6-2 Polarization and Flux Measurements ...... 100
7 LIST OF FIGURES Figure page
1-1 The formation of a low mass star...... 20
1-2 Illustration of the structures and physical processes in a protoplanetary disk. .... 21
1-3 Illustration of absorptive and emissive polarization in magnetically dust alignment theory...... 22
2-1 The frame of the object with o and e ray images...... 28
2-2 CanariCam wave-plate efficiency and Instrumental Correction...... 29 3-1 Total intensity maps of WL 16 at 8.7, 10.3, and 12.5 µm...... 35
3-2 The 8.7- µm linear polarization map of WL 16 superimposed on (total intensity) contours...... 36
3-3 The low-resolution (R≈50) spectrum of the brightest central 1′′.6 (21 pixels) region of WL 16...... 37
3-4 Polarization and emission/absorption decomposition results of WL 16...... 38
3-5 Polarization and emission/absorption decomposition results of Elias 29...... 39
3-6 The JHK color-color diagram...... 41 3-7 Comparison of the polarization profiles of WL 16 (black) and Elias 29 (blue). .... 42
3-8 Intensity and charge state of PAHs at two sides of the disk...... 47
4-1 Intensity map (contours) of MWC 1080 system at 11.2 µm...... 61
4-2 Intensity and polarization spectra of NW nebula...... 62 4-3 Signal-to-noise (S/N) ratio of polarized intensity, Stokes u(U/I ), and Stokes q(Q/I ) of the NW nebula...... 63
5-1 Ratio of the absorptive efficiency perpendicular and parallel to the shorter axis of the spheroid Qabs,⊥/Qabs,∥ covering wavelengths λ from 1.0 to 1000.0 µm...... 79
5-2 Ratio of the absorptive efficiency along two directions Q⊥/Q∥ vs changing grain radius a...... 80
5-3 Emissive polarization profiles vs wavelength...... 81 5-4 The product of the degree of polarization at 90◦ due to single scattering (P) and the dust albedo (ω) for dust size distributions with different values of amax (the value of amin is fixed at 0.01 µm)...... 82
8 5-5 Simulated linear polarization maps in a poloidal B-field configuration for a spherical envelope at λ=10.0 µm...... 83
5-6 Simulated linear polarization map in a toroidal B-field configuration for a spherical envelop at λ=10.0 µm...... 84
5-7 Dust scattering induced polarization map at λ=10.0 µm...... 84
5-8 Linear polarization map of a poloidal shape B-field configuration...... 85
5-9 Linear polarization map of a toroidal shape B-field configuration...... 86
5-10 Linear polarization map of an hour-glass shape B-field configuration...... 86 5-11 Linear polarization map of Aitken model IVa B-field configuration...... 87
5-12 Linear polarization map of Aitken model V B-field configuration...... 87
5-13 Linear polarization map of Aitken model VI B-field configuration...... 88
5-14 Linear polarization map of a dipole shape B-field configuration...... 88
5-15 Observation and modeling of AB Aur...... 89 6-1 Composite image of W51 IRS2...... 101
6-2 Polarization maps of W51 IRS2 at 8.7, 10.3, and 12.5 µm filters from top to the bottom...... 102 6-3 Decomposition results of each source in W51 IRS2 region...... 103
6-4 Maps of magnetic field orientations...... 104
6-5 Magnetic fields overlaid with 14GHz continuum emission from VLA...... 105
A-1 Polarization profiles...... 110
9 Abstract of Dissertation Presented to the Graduate School of the University of Florida in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy MID-INFRARED POLARIMETRY: INSIGHTS INTO MAGNETIC FIELDS AND DUST GRAIN PROPERTIES IN YOUNG STELLAR OBJECTS AND PROTOPLANETARY DISKS
By
Han Zhang August 2018
Chair: Charles M. Telesco Major: Astronomy
Polarization observations in star formation regions have been made covering different wavelengths, different scales, and different stellar evolution stages, from embedded cores to disks. One of the important goals of polarimetry observations is to explore the role of magnetic
fields in star formation regions, which are believed to be important in regulating star formation.
However, recent theoretical work and observations challenge the ability of polarimetry as an important technique to trace the magnetic field morphology, as other non-magnetic mechanisms, e.g., dust scattering, can also produce significant polarization signals.
This dissertation work focuses on two fundamental questions concerning polarimetry observations in star formation regions, (i) polarization mechanisms at mid-infrared wavelengths, and (ii) interpretation of the inferred magnetic field structure. We address these two questions using simulations and observations from the mid-infrared facility camera CanariCam on Gran
Telescopio Canarias (GTC), a 10.4-meter telescope in La Palma, Spain. To explore the origins of the polarization: 1) we study WL16, a Herbig star with an extended disk, and we find that the contamination from foreground extinction accounts for our detected polarization; 2) we model dust-scattering-induced polarization in the disk, reproducing the observations of AB Aur, a Herbig Ae star, and thereby constrains the size of dust grains to be micron-sized; 3) we discover ∼2% polarization from interstellar polycyclic aromatic hydrocarbon molecules, in the nebulous structure surrounding MWC 1080.
10 To be able to infer the magnetic field structure in protoplanetary disk from polarimetry observations, we conduct simulations of polarized disks expected for a variety of magnetic field structures, e.g, poloidal or toroidal. In the case of W51 IRS2, a star cluster formation region, by comparing the inferred magnetic fields with the ionized gas flow, we find that the gas flow is likely to be driven by the magnetic fields.
11 CHAPTER 1 INTRODUCTION
With efforts and studies over centuries, we have gradually come to know our place in
the Universe, one with billions upon billions of stars. How was our Sun formed? How did the Earth come to be? Motivated by answering these questions, we look at an abundant sample of
star-planet systems at the early phases of their lifetime in the Universe.
1.1 Magnetic Fields in Star Formation
Star formation is expected to start from the dense regions in giant molecular clouds, by
gravitationally contraction (Shu et al., 1987; McKee & Ostriker, 2007). Most stars are born in groups or clusters, and very few are born alone (Lada & Lada, 2003). When the stable
dense core is sufficiently disturbed or the mass of the core exceeds the Jeans mass, it tends to
collapse to form a star (Lada, 2005). Magnetic fields, together with turbulence and gravity, are
key forces in controlling star formation (McKee & Ostriker, 2007).
The ratio of gravitational energy to magnetic energy (the mass-to-flux ratio), M/Mϕ √ (where Mϕ≡ϕ/2π G, ϕ is the magnetic flux), determines whether the cloud is magnetically
supercritical or subcritical. In a magnetically subcritical cloud (M when the core mass M reaches and exceeds the Mϕ, the core becomes supercritical and the collapse can happen. The ratio M/Mϕ may vary from place to place within the cloud complexes, and ambipolar diffusion may convert a subcritical cloud to a supercritical state. The relative importance of turbulent energy to magnetic energy, which parameterizes whether the gas is sub-Alfv´enicor super-Alfv´enic.It can be studied by analyzing the linear polarization dispersion in the molecular clouds (detailed discussion of the technique is in Section 1.3). The magnetic fields inferred from the linear polarimetry appear to be fairly uniform and tend to preserve their orientations at different scales (Li et al., 2014b), suggesting that magnetic fields are relatively stronger than the turbulence, in other words, the turbulence is sub-Alfv´enic.In 12 summary, magnetic fields are important in the evolution of molecular clouds, and possible observational probes have been discussed in Li et al. (2014b). In the canonical star formation theory (Shu et al., 1987), the original ordered magnetic fields lines in the cloud will be pinched to form an hourglass morphology along the major axis of the dense core (Galli & Shu, 1993). The strength of the pinch is expected to be larger in weak magnetic fields than strong magnetic fields (Crutcher, 2012). Girart et al. (2006) report a clear hour-glass shaped B-field morphology in a binary protostellar system NGC 1333 IRAS 4A. The hour-glass morphology is expected particularly for contracting clouds with flux-frozen magnetic fields (Galli & Shu, 1993). Simulations suggest that a toroidal magnetic field dominates during the accretion phase after the protostar forms (Safarzadeh et al., 2017; Cho & Lazarian, 2007; Tomisaka, 2011). Determining magnetic field morphology is therefore an important observational test of star formation theory (Shu et al., 1987; Crutcher, 2012). Both analytic arguments and ideal magnetohydrodynamic (MHD) simulations (Machida et al., 2011) imply that, in the gravitationally collapsed cloud, disk formation, once thought to be the natural consequence of the conservation of angular momentum, is suppressed in the presence of magnetic fields, if the magnetic fields are aligned with the infalling envelope’s rotation axis (Li et al., 2014c). Possible solutions include non-ideal MHD effects (ambipolar diffusion, Ohmic dissipation and Hall effect), turbulence, misalignment between the magnetic field and rotation axis in weakly magnetized cores (Li et al., 2013, 2014c). The last mechanism has been investigated observationally by Chapman et al. (2013) on the ∼10,000 AU scale, Hull et al. (2013) on the ∼1000AU scale and Davidson et al. (2011), but with different conclusions, perhaps because of the different observation scales. The role of magnetic fields in star formation and planetary formation is an on-going topic of study (Crutcher, 2012; Reissl et al., 2016a). 1.2 Protoplaneatry Disks Star formation starts with the collapse of dense regions in molecular clouds. An accretion disk is developed during the collapse phase, as a consequence of removal of high angular 13 momentum material. The existence of circumstellar disks has been confirmed by, for example, Hubble Space Telescope observations, where the disks are seen as shadows against the bright nebular background. Conventionally, Young Stellar Objects (YSOs) have been classified into four different stages, Class 0, Class I, Class II, and Class III, depending on the appearance of their spectral energy distribution (SED) (Fig.1-1)(Lada & Wilking, 1984; Andre et al., 1993). This classification scheme is an empirical classification based on the spectral index in the near-infrared (NIR) and mid-infrared (mid-IR). It has been found to quantitatively fit into the theoretical scenario of star formation (Adams et al., 1987) consisted with four evolutionary stages that young stars go through. During the first stage of star formation, the Class 0 YSOs are deeply embedded within the optically thick gas and dust, and Class 0 sources are almost invisible in the optical and NIR wavelengths but have significant sub-millimeter emissions (Andre et al., 1993). At this stage, a rotationally supported disk may not be formed, and the fattened inner envelope is called a pesudodisk, extending to a few thousand AU (Chapman et al., 2013). A rotational signature of a disk is clearly detected in the Class I YSO, and it is still embedded in an envelope of infalling material, which is usually associated with strong outflows or jets. The SED of a Class I YSO shows rising infrared continuum that peaks at mid-to-far infrared wavelengths (Lada, 1987). In the third stage (Class II; (Williams & Cieza, 2011)), as the accretion proceeds, the envelope is eventually lost and the disk and central star become optically visible. The individual star and disk systems evolve differently depending on the mass of the stars and environment. The disk mass now is only a few percent of the central stellar mass and can be called a protoplanetary disk. The life time of the disk is of order several million years (Ribas et al., 2015). Once the gas disk has dissipated, and the central star reaches the main-sequence, it is the Class III phase, in which the disk evolves to a debris disk that is optically thin and gas poor, with larger dust grains, planetesimals and/or planets. Protoplanetary disks undergo a variety of processes, including dust aggregation, dust settling, dynamical interaction with (sub)stellar companion and photo-evaporation (Williams & 14 Cieza, 2011). Thanks to the multi-wavelength high-resolution and high-sensitivity observations of circumstellar disks, we are able to see the subtle features in these disks, especially through the polarimetric imaging in NIR and (sub-)millimeter continuum and gas observations. Multi-wavelength observations of protoplanetray disks reveal emission from different regions of the disk, illustrated in Fig.1-2. Disks are usually optically thick in the NIR and mid-IR, and the emission is mainly from dust grains at the surface layer and near the surface layer, respectively (Mari˜naset al., 2006), while the (sub-)millimeter observations trace larger and cooler dust grains at the disk mid-plane. Intriguing asymmetric features in protoplanetary disks, such as spiral arms, gaps and warps, have been modeled as a result from the perturbation of the embedded planet(s) (Zhu, 2015; Dong et al., 2015a), e.g., MWC 758 (Benisty et al., 2015) and HD 142527 (Fukagawa et al., 2006). Observations by Atacama Large Millimeter/submillimeter Array (ALMA) are beginning to reveal the structure of protoplanetary disks at their midplanes. With the revolution starting with the fabulous ring structures seen in HL Tau (ALMA Partnership et al., 2015), a structure that can be reproduced by the simulation including perturbation by multiple planets in the disk. Observations of, e.g., TW Hydrae (Andrews et al., 2016) also show ring structures, and in Elias 2-27, spiral arms are detected (P´erezet al., 2016; Meru et al., 2017). Magnetic fields must play a crucial role in the evolution of protoplanetary disks and the formation of planets. Magnetic fields influence the transport of dust and gas and even the migration of planetesimals and planets (Bertrang et al., 2017). In magnetized protoplanetary disks, MHD turbulence driven by the magnetorotational instabiliy (MRI) provides the source of viscosity, determining accretion rates and lifetime of the disks (Wardle, 2007). However, the observational constraints on magnetic fields remain a difficult task, as we discussed in the following. 1.3 Polarimetry: An Important Technique to Study Magnetic Fields Even though magnetic fields are important, it is very difficult to directly quantify magnetic fields observationally (Robishaw, 2008). There are several observational techniques that are 15 available for observing magnetic fields. Zeeman splitting of spectral lines is used to measure the strengths and directions of magnetic fields along our line of sight. Spectral lines will be split in a region that is embedded in magnetic fields and the splitting is proportional to the strength of magnetic fields along the line of sight. However, the interstellar magnetic fields are usually very weak, which cannot produce strong enough magnetic momentum for the detection of line splitting. Faraday rotation (Akahori & Ryu, 2010), effects the light as it propagates through the interstellar medium, with the plane of linear polarized background radiation being rotated in the magnetic fields by an angle β, which is proportional to the strengths of magnetic fields along our line of sight, the electron density, and the wavelength of the radiation (Crutcher, 2012). The fingerprint of magnetic fields is also embedded in the polarization of the received interstellar radiation. Dust-induced polarization is widely used as a tracer of magnetic fields in star formation regions and is linked to the strengths and orientations of magnetic fields (Matthews et al., 2009). Optical interstellar polarization was first reported independently by Hall (1949) and Hiltner (1949), and recognized as due to the preferential dichroic extinction of the background starlight. The reason why interstellar dust grains can produce polarized light is because of their alignment. The details of the alignment mechanism for dust grains are still unclear. The classical dust alignment theory suggests that non-spherical dust grains tend to spin, perhaps as a result of radiative alignment torques from an anisotropic radiation field (Lazarian, 2007; Cho & Lazarian, 2007). In the radiative alignment theory, when light shines on a irregular dust grain it starts to rotate. The grains will precess around the local magnetic field. The process has been explained in a series of papers, e.g., Hoang & Lazarian (2014) and Cho & Lazarian (2007), is the most favored and widely accepted mechanism. Consequently, as shown in Fig. 1-3, in the presence of an external magnetic field, these spinning non-spherical grains tend to align with their longer axis perpendicular to the magnetic field lines. In this way, at optical and NIR, aligned elongated dust grains polarize the background starlight and the resultant polarization is along the direction of interstellar 16 magnetic fields (Left panel in Fig. 1-3), e.g., Goodman et al. (1990) and Cashman & Clemens (2014). At far-infrared and (sub)-millimeter, aligned dust grains emit polarized emission with the direction perpendicular to the magnetic fields (Right panel in Fig. 1-3). In the latter case, to infer the magnetic field directions, we need to rotated the detected polarization vectors by 90◦, e.g., Stephens et al. (2014). Mid-IR polarimetry is an incisive probe of B-field morphologies. At mid-IR wavelengths, the situation becomes more complicated, since the observed polarization can be a combination of dichroic absorption, emission, and/or scattering (Aitken et al., 2004). Scattering polarization is discussed in Section 1.4. Using template absorption and derived emission polarization profiles, Aitken et al. (2004) proposes a method to separate the relative contributions from dichroic absorption and emission (Appendix A). Separating these two components offers an opportunity to probe directly the line-of-sight magnetic fields structure, unique to mid-IR, and we call it ‘polarimetric tomograpy’. Barnes et al. (2015) has demonstrated the potential of this method by probing the line-of-sight magnetic fields structure of the massive star formation region K3-50. We also emphasize that beside being a tracer of magnetic field, the degree of polarization depends on the geometry and composition of dust grains, and can therefore provide information about dust properties (Cho & Lazarian, 2007) (see Chapter 5). 1.4 Challenges: Scattering, Radiative Grain Alignment The origin of the polarization signal is still not fully understood. Besides the classical dust grain alignment in magnetic fields, there have been two other mechanisms proposed that can produce polarized emission, but are irrelevant to the study of magnetic fields. These are dust scattering (Kataoka et al., 2015), and radiative grain alignment (Tazaki et al., 2017). Even though scattering was regarded as negligible at wavelengths longer than NIR, both observations and theoretical studies indicate that scattering can not be ignored if there are dust grains with sizes comparable with observation wavelengths. Polarization observations towards the circumstellar disks around HD 142527 (Kataoka et al., 2016b), HL Tau (Stephens 17 et al., 2017; Kataoka et al., 2017, 2016a) and IM Lup (Hull et al., 2018) by ALMA show evidence of polarization due to scattering instead of magnetical dust alignment. Kataoka et al. (2016b) finds a sharp change in the directions of polarization vectors from the inner disk to outer disk in HD 142527, which is hard to explain in the context of magnetically aligned dust, but is consistent with the dust scattering. Stephens et al. (2017) shows that the polarization morphologies of HL Tau at three ALMA bands (band 3, band 6 and band 7) are different, which indicates a change in the dominated polarization mechanism(s) with wavelengths. And note that scattering polarization is very dependent on the wavelengths. Dust scattering produces polarization vectors along the azimuthal direction of the disk for a face-on disk (Yang et al., 2016). However, scattering depends on the dust properties, the geometry of the disks, inclinations, and anisotropy of the radiation, the scattering morphology can be complex. In the case of IM Lup, polarization vectors are roughly along the minor axis, which is indicative of polarization produced by dust scattering (Hull et al., 2018). If the observed polarization patterns are indeed from scattering-induced polarization, instead of providing information about the orientations of magnetic fields, it provides diagnostic constraints on the size of dust grains in protoplanetary disk, and evidence for dust growth (Kataoka et al., 2016a). Tazaki et al. (2017) revisits the radiative torque alignment theory in protoplanetary disks and finds that larger dust grains are expected to be aligned by the radiation flux at mid-plane, which would correspond to the polarization vectors tracing the direction of radiation fields instead of the magnetic fields. For a face-on disk, the polarization pattern can be in the azimuthal direction if the dust is aligned by the radiation fields, even in the presence of magnetic fields. Kataoka et al. (2017) finds that polarization vectors of HL Tau are along azimuthal direction at ALMA band 3, which supports the radiative torque alignment. Of course, the discussion of the dust alignment is not limited to these two. For instance, if grains are mechanically aligned by outflows (Lazarian & Hoang, 2007), then the polarization is expected to be parallel to the B-field (Hull et al., 2013). 18 There are various physical effects that can produce polarized radiation at infrared and millimeter wavelengths, and interpreting the polarization data is complicated. Multi-wavelength polarimetric observations are essential and significant to improve our understanding of the polarization signal, such as the multi-band polarimetry observations of HL Tau by ALMA (Stephens et al., 2017). We expect future polarimetry observations to solve the problem with observations of more protoplanetary disks systems. 1.5 Thesis Outline My dissertation is organized as follows. In Chapter 2, we give a brief introduction to the mid-IR facility Canaricam on GTC, and then describe the data reduction. In Chapter 3, we show the polarimetric imaging and spectra of a well-resolved protoplanetary disk, WL 16. In Chapter 4, we report the detection of the PAH polarization in the nebula around the stellar cluster MWC 1080. Chapter 5 displays a library of maps of linear polarization of a protoplanetary disk at 10 µm. In Chapter 6, we show the polarimetry imaging of a protocluster W51 IRS2 and introduce the method of “polarimetric tomography” to infer the magnetic field structures in the region. We summarize this thesis work in Chaper 7. 19 Figure 1-1. The formation of a low mass star. The plot is adapted from Dauphas & Chaussidon (2011). Left: SED signatures of five different phases, from prestellar core to Class 0-III YSOs. Right: Illustration of different stages. From top to bottom: star formation starts from the collapse of a dense molecular core. A protostar is formed (Class 0 YSO). Rotation disk and the star is embedded in an envelope (Class I YSO). The envelope dissipates and the disk is visible (Class II YSO). The disk continues to evolve and becomes a debris disk (Class III YSO). 20 MIR Scattered Light (sub-)mm 1 2 3 4 a b d c Distance in AU 1 10 100 1 Turbulent Mixing (radial or vertical) 2 Vertical Settling 0.35 mm 3.0 mm ALMA 3 Radial Drift 10 µm VLTI/MATISSE 4 a) Sticking b) Bouncing 2 µm 10 µm EELT c) Fragmentation with mass transfer d) Fragmentation JWST/MIRI Figure 1-2. Illustration of the structures and physical processes in a protoplanetary disk from Testi et al. (2014). The right side shows the regions of the disk that can be probed at various wavelengths. Near-infrared, mid-IR and sub-mm emission are from dust grains at the very surface layer, below surface-layer, and mid-plane of the disk, respectively. The left side shows the dust transportation and collision mechanisms. 21 Figure 1-3. Illustration of absorptive and emissive polarization in magnetically dust alignment theory, adapted from Lazarian (2007). (a)Left panel – Polarization of starlight passing through a cloud of aligned dust grains. The direction of polarization is parallel to the plane of the sky direction of magnetic field. (b)Right panel– Polarization of radiation from a optically thin cloud of aligned dust grains. The direction of polarization is perpendicular to the plane of the sky direction of magnetic field. 22 CHAPTER 2 OBSERVATIONS AND DATA REDUCTION 2.1 CanariCam CanariCam is a multi-mode mid-IR (8-25 µm) facility camera developed at UF (Telesco et al., 2003), which started science operation in 2012 at the 10.4-meter Gran Telescopio Canarias (GTC), on La Palma, Spain. GTC is a partnership of Spain, Mexico, and the University of Florida, and located at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofs´ıcade Canarias, on the island of La Palma. CanariCam was mounted at the GTC Nasmyth-A platform. It has a 320×240 Raytheon array, corresponding to a pixel scale of 0′′.08 and a field of view of 25′′.6×19′′.2. It allows Nyquist sampling at 8 µm. The theoretical angular resolution of the camera at 10 µm is 0′′.25 (1.22λ/D), but the seeing at the GTC degrades this somewhat, often being 0′′.3–0′′.4, depending on observing conditions. CanariCam is equipped with a suite of narrow-band and broad-band 10 and 20 µm filters and provides the low and moderate-resolution (R=100 –1300) slit spectroscopy in the 10 and 20 µm regions. It permits a convenient switch between different modes, imaging, spectrospopy, polarimetry and coronagraphy. It uses the standard chop-nod techniques to account for the thermal background radiation and the offset in the thermal background by the telescope. In this thesis work, observations are mainly conducted using the polarimetry mode provided by CanariCam (both imaging and spectropolarimetry). We will discuss in detail the polarimetric observations and polarimetric data reduction with CanariCam in Section 2.3. 2.2 Stokes Parameters The stokes parameters quantify the polarization properties of the radiation field and can fully describe the characteristics of a partially polarized light. Single electromagnetic waves traveling in the z direction can be describe by i(kz−ωt+ϕx ) Ex (z, t) = E0x e (2–1) i(kz−ωt+ϕy ) Ey (z, t) = E0y e (2–2) 23 Stokes parameters are defined as: ⟨ ⟩ ⟨ ⟩ ∗ ∗ I = Ex Ex + Ey Ey (2–3) ⟨ ⟩ ⟨ ⟩ ∗ − ∗ Q = Ex Ex Ey Ey (2–4) ⟨ ⟩ ⟨ ⟩ ∗ ∗ U = Ex Ey + Ey Ex (2–5) ⟨ ⟩ ⟨ ⟩ ∗ − ∗ V = i( Ex Ey Ey Ex ) (2–6) Or can be written as ◦ ◦ ◦ ◦ I = I0 + I90 = I45 + I135 = IRHC + ILHC (2–7) ◦ − ◦ Q = I0 I90 (2–8) ◦ − ◦ U = I45 I135 (2–9) V = IRHC − ILHC (2–10) where I denotes the total intensity, Q and U describe the two-dimensional state of linear polarization, V describes circular polarization, RHC is right-hand circular light, and LHC is √ left-hand circular light. The degree of polarization is then: p = Q2 + U2 + V 2/I . When a unpolarized light is incident on a briefringent material, it will split into two polarized components, one with its vibration direction along the fast axis, called the ordinary ray, and the other in the direction of the slow axis, and called the extraordinary ray. The directions of the oscillations of o ray and e ray are orthogonal to each other, their intensities are often written as Iα and Iα+90. The sum of the o and e ray intensities equals to the total intensity, i.e., the Stokes parameter I. 2.3 Data Reduction For polarimetric observations, the incident beam propagates within CanariCam as follows (see CanariCam user manual). First the light passes through the focal plane mask. The focal plane mask is used to eliminate overlapping of the orthogonally polarized beams on the array. It blocks 50% of the field of view, and the field of view in the polarimetry mode is therefore 24 reduced to 26′′× 9′′.6 made up of several smaller pieces, each of which is 26′′× 3′′.2. Then the light passes through the half-wave-plate. The half-wave-plate is used to modulate the position angle of polarization by the Wollaston prism rather than rotating the instrument. A mechanical rotation of θ results in an optical rotation of 2θ, i.e., rotating the half-wave plate by 22.5◦ is equivalent to rotating the analyzer by 45◦. With the half wave plate at 0◦ and 22.5◦, both Stokes Q and U can be derived based on Eqs.7-10, but two independent confirmation at angles of 45◦ and 67.5◦ are recommended to correct for the sensitivity difference in the two channels of the polarimeter. The light beam finally passes through the Wollaston prism, which is made of the birefringent crystal calcite, where the light is split into the two orthogonally polarized beams, o and e rays. These two rays are displayed on the detector plane simultaneously, as shown in Fig.2-1. In the field of view, two perpendicular polarized images are present and the mask is used to prevent the overlapping of the images. I adapt the data reduction procedures presented in Li (2014) and Berry et al. (2005), and describe how we reduce the imaging and spectropolarimetry data in detail below. The data reduction package can be found at github. 2.3.1 Polarimetric Imaging Data Reduction 1. Read in the FITS file image A FITS file generated by CanariCam has multiple extensions. The first extension (i.e., Ext. 0) only contains a ASCII header, which stores the complete information, such as the object coordinates, observation time, filters, and the history of data. For the other extensions (Ext. 1, 2, 3 etc), each contains a 3-dimensional array of 320×240×2, where 320×240 is the size of the field of view in pixels and 2 corresponds to a nod pair, one is with source and the other is off-source. To get the emission from the target, we need to subtract the off-source frame from the frame with source. In the polarimetry mode, half-wave-plate is rotated to four orientations 0, 22.5, 45, and 67.5 degrees and that is stored in the parameter nsave. After reading in the FITS file, we obtain four images, corresponding to four wave plate orientations, and each image contains a 320×240 array. 2. Extraction of the o and e ray images The entire field of view of a polarimetry image is divided into six slots (i.e., three o and e ray pairs) by the polo-mask, shown in Fig. 2-1 (we use the data of W51 IRS2 as an illustration). The next step is to separate and extract o and e ray images. o ray and e ray images need to be identical and save into two separate arrays for further processing. The o and e ray images have a 35 pixels (2′′.8) separation because of the position of the 25 mask. We divide the image into two sub images, so o and e ray images are separated. Then we use the centroid function to find the center of o and e ray images and crop the image into sizes that we want. The cropped o and e ray images are displayed in the right panel in Fig. 2-1. 3. Image alignment Now we have eight images, corresponding to o and e ray images at four different wave-plate angles. The precise alignment between these different image is important to derive the accurate polarization measurements. To register the images, we used the cross-correlation algorithm, which outputs the offsets between images depending on the location of the peak in the cross-correlation matrix. We select a frame as the reference, by default the o ray image at wave-plate angle 0◦. A pair of o and e ray images before and after registration is plotted in the right panel in Fig. 2-1. 4. Stokes parameters The Stokes parameters I, Q, and U provide a full description of the linear polarization. The mathematical connection between the Stokes parameters and the intensities are described in the following equation. I0 − I90 = Ipcos(2θ) = Q, I0 + I90 = I (2–11) I45 − I135 = Ipsin(2θ) = U, I45 + I135 = I (2–12) Where Iα represents the intensity of the o ray image, and Iα+90 is the intensity of the e ray image. Normalized Stokes parameters q and u are defined as q=Q/I and u=U/I . By subtraction and addition between o and e ray images at different wave-plate angles, we are able to derive the normalized stokes parameter, q and u. 5. Correction for instrumental effects The tertiary mirror of the telescope introduces significant instrumental polarization (IP). To characterize the degree of instrumental polarization, we select a sample of intrinsically unpolarized stars from Cohen standard (Cohen et al. 1999)and measure their polarization, which should come totally from the instrument. We find that the amplitude of the instrumental polarization is between 0.6% and 0.7%. The position angle of the instrumental polarization is given by the empirical equation θIP =ZD-RMA-31.9, where ZD (zenith distance) and RMA (rotator mechanical angle) are recorded in the FITS header. u = uobs − pIP cos(2θIP ) (2–13) q = qobs − pIP sin(2θIP ) (2–14) We subtract the instrumental effects from the observed Stokes uobs and qobs . Wave-plate efficiency also needs to be corrected to account for the loss in transmission, and the polarimetric efficiency profile is plotted in the right panel in Fig.2-2. 6. Calculation of the polarization percentage and polarization angle 26 With instrumental effect corrected Stokes parameter u and q, the polarization degree and orientation with corresponding uncertainties are given √ √ 2 2 − 2 2 2 p = q + u σp, σp = (udu) + (qdq) /p θ = 0.5atan(u, q), σθ = σp/2p The uncertainties σq and σu associated with the normalized Stokes parameters are derived using a 3-sigma clipping algorithm (Robinson, 1987), and they are then propagated through the analysis to obtain the polarization uncertainty σp and polarization PA uncertainty σθ = σp/2. The obtained polarization position angle needs to be calibrated to the North. There are two approaches, one is to use the polarization calibration star, with known polarization position angles, given in Smith et al. (2000). The other way to use the empirical relationship displayed in the left panel in Fig. 2-2. We examine a large sample of calibrations stars and find a constant offset value between the observed polarization position angles and the true position angle, relative to the North. For the imaging flux calibration, we use the flux-calibrated standard stars in Cohen et al. (1999). Aperture photometry of the standard star is performed in ADUs and then we compute the ratio between a number of counts (ADUs) per second and the real flux of the star. We apply the ratio to the target images and derive the flux of the science object. The accuracy of photometry is ∼10% using this method. The final polarization parameters are usually presented graphically in the form of a polarization map. Each measurement is represented as a vector, in which the length of the vector is proportional to the degree of polarization, and the orientation of the vector is the polarization position angles. The vectors are usually binned to reduce the number of vectors and better show the results. 2.3.2 Spectropolarimetry Data Reduction The first steps of spectropolarimetry data reduction are similar to the polarimetric imaging reduction. After correctly stacking the o and e ray 2 dimensional spectra, which are rectangular regions containing all source flux, we only need to do extraction and calibration, following the procedures in Gonz´alez-Mart´ınet al. (2013). 1. Extract the one-dimensional spectra from the two-dimensional spectra We extracted one-dimensional spectra by integrating the central region that we are interested along the slit direction. The width of the extraction region depends on whether it is a point-like or extended sources. Usually we integrate across the whole region to get the maximum signal. 2. Wavelengths and flux calibration We use sky emission features in the reference stacked raw spectra to do the wavelength calibration. There are about 8 sky lines (refer to GTC/CanariCam website) with known 27 200 150 Pixels 100 50 0 0 50 100 150 200 250 300 Pixels Figure 2-1. Left: The field of view of CanariCam in the case of W51 IRS2 in polarimetry mode. It has 320×240 pixels. The o and e ray images have a 35 pixel separation with a mask. Right: Comparison of reference image and off-set images before and after the alignment. Each panel contains the reference image, offset image, and the cross-correlation. peak wavelengths in N bands. After identifying at least three lines, we fit a second-order (or higher) polynomial function for the wavelength calibration. For the flux calibration, we use the ratio between the observed associated standard and its corresponding flux calibrated template given by Cohen et al. (1999). The final flux-calibrated spectrum is obtained by dividing the spectrum (in counts; ADUs) by this ratio. 28 1.00 Error-weighted mean = 120.3 +/- 0.8 0.95 150 0.90 100 0.85 Efficency 0.80 50 0.75 Calibration offset (degree) 0.70 0 7 8 9 10 11 12 13 Wavelength (μm) 0 5 10 15 20 25 30 File # Figure 2-2. CanariCam wave-plate efficiency and instrumental correction. Left panel: The wave-plate efficiency shows the trend through the N-band as window. Polarization measurements should be corrected for that. For imaging polarization at Si2 (8.7 µm), Si4 (10.3 µm), and Si6 (12.5 µm), we use 0.9, 0.99, and 0.88 to calibrate the efficiency. Right panel: The constant polarization position angle calibration value in the sample. 29 CHAPTER 3 THE MID-INFRARED POLARIZATION OF THE HERBIG AE STAR WL 16: AN INTERSTELLAR ORIGIN? 3.1 Introduction Magnetic fields (B-fields) play an important role in almost all stages of star formation as discussed in the comprehensive star formation review by McKee & Ostriker (2007). However, there are still many uncertainties about how B-fields regulate the protoplanetary disk formation and evolution. For example, magnetically driven core-collapse models (Shu et al., 1987; Galli & Shu, 1993) predict an hourglass-shaped B-field geometry at early stages in the disk evolution, a scenario challenged by recent observations (e.g. Chapman et al., 2013; Davidson et al., 2011; Hull et al., 2013). The only way to address these issues is with high angular resolution observations of the B-field morphologies in a variety of young disks and their environments. Polarimetry is a potentially incisive observational probe of B-field morphology (Barnes et al., 2015; Crutcher, 2012; Matthews et al., 2009; Hull et al., 2014). Dust grains can polarize light by scattering, dichroic absorption, or dichroic emission, the latter two processes attributed to non-spherical dust grains with their long axes aligned perpendicular to the B-field lines, perhaps as a result of radiative torque (Lazarian, 2007; Hoang & Lazarian, 2014). Dichroic absorption by aligned non-spherical dust grains can partially polarize background starlight with the transmitted E-field direction parallel to the aligning B-field lines (Cashman & Clemens, 2014). Light scattering by spherical and non-spherical dust grains can produce high fractional polarization prominently in the optical and near-IR. At far-IR and sub-mm wavelengths, aligned non-spherical dust grains emit polarized light with the emitted E-field direction being perpendicular to the B-field lines. In the mid-infrared (mid-IR, 8–30 µm), the situation becomes more complicated, since the observed polarization can be a combination of dichroic absorption, emission, and/or scattering (Aitken, 2005). Nevertheless, mid-IR polarimetry has some distinct advantages compared to other wavelength regions: the predicted emissive polarization is much larger than that at longer wavelengths (Cho & Lazarian, 2007), and we can potentially disentangle both the absorptive and emissive polarization components 30 simultaneously along the line of sight and thereby infer the three-dimensional structure of the B-fields (Barnes et al., 2015; Smith et al., 2000). Herbig Ae/Be stars (1.5 (Herbig, 1960; Hillenbrand et al., 1992). WL 16 was discovered by Wilking & Lada (1983) and identified as a late-stage Herbig Ae star embedded in the ρ Ophiuchus molecular cloud L1688 (Ressler & Barsony (2003), herafter RB03). Key physical properties of WL 16 are given in Table 3-1. No outflow and only weak 1.3 mm emission are observed for this object (RB03). Because of its high extinction (AV ≈ 25–30 mag), WL 16 is undetectable in the optical but displays extended, resolved emission in the mid-IR. The most plausible interpretation of the extended emission is that it is a disk with a diameter of nearly 900 AU. Kinematic modeling of near-IR CO vibrational overtone emission arising from the innermost region suggests that there is indeed a flat Keplerian gas disk (Carr et al., 1993; Najita et al., 1996). Nevertheless, whether the entire structure is an unusually large protoplanetary disk or a smaller disk associated with a remnant of the collapsing envelope is still unclear. The very extended emission is generally ascribed to a population of very small grains (VSGs) and polycyclic aromatic hydrocarbons (PAHs). The latter population is indicated by the IR spectra of WL 16 rich in PAH emission features (DeVito & Hayward, 1998; Geers et al., 2007). Because of their small size (several dozens to hundreds of carbon atoms), PAHs are dynamically well coupled to the gas and less affected by dust settling. The PAHs observable in the mid-IR likely reside in the disk surface layers. PAHs can also be easily ionized by stellar UV radiation due to their low ionization potentials, and ionized PAHs can be the tracer of low density optically thin regions (Maaskant et al., 2014). In this paper, we present high-resolution (∼0′′.4) mid-IR polarimetric imaging and spectropolarimetry of WL 16. In addition, we find it illuminating to consider complementary observations of Elias 29 (our calibration star), which is a low-mass luminous Class I protostar (36 L⊙) and a neighbor of WL 16 in the ρ Ophiuchus molecular cloud (Boogert et al., 2002). 31 Table 3-1. Basic Properties of WL 16 Properties Values References Distance 125 pc RB03 Stellar mass 4 M⊙ RB03 Luminosity 250 L⊙ RB03 Disk inclinationa 62.2◦ 0.4◦ RB03 Disk PAb 60◦ 2◦ RB03 −6.8 −1 Accretion rate 10 M⊙ yr Natta et al. (2006) Disk mass ¡0.001 M⊙ Andrews & Williams (2007) Disk diameter 900 AU RB03 a 0◦for face-on b Position angle of the major axis of the disk The paper is organized as follows. In Section 3.2 we describe our observations. In Section 3.3, we discuss the origin of the observed polarization from polarimetric imaging and spectroscopy, with detailed consideration given to the contribution associated with the foreground extinction and B-field. In Section 3.4 we examine the spatial distribution of the ionized and neutral PAHs along the major axis of the disk and draw attention to morphological features that may provide clues to WL 16’s dynamical evolution. Finally, in Section 3.5, we summarize this work. 3.2 Observation CanariCam is the mid-IR (8–25 µm) multi-mode facility camera of the 10.4-meter Gran Telescopio Canarias (GTC) on La Palma, Spain (Telesco et al., 2003). It has a field of view of 26′′ × 19′′ with a pixel scale of 0′′.079, which provides Nyquist sampling at 8.7 µm. In the polarimetry mode, the actual field of view is reduced to 26′′ × 2.′′6 after a focal mask (to avoid overlapping between o and e beams) is used. We obtained polarimetric images of WL 16 in three filters near 10 µm on 6 August 2013 and a low-resolution (R≈50) polarimetric spectrum of WL 16 from 7.5 to 13 µm on 4 and 5 July 2015 (Table 3-2), as part of the CanariCam Science Team guaranteed time program (PI: C. M. Telesco) at the GTC. The imaging and spectroscopic observations of WL 16 were interlaced with observations of standard star HD 145897 (Cohen et al., 1999) for flux and point-spread-function (PSF) calibration, and the standard star Elias 29, selected from Smith et al. (2000) to calibrate the polarization position angle (PA). The achieved angular resolution (full-width at half maximum intensity) for the 32 polarimetric images was 0′′.4–0′′.6 (Table 3-2). CanariCam was rotated so that the polarimetry field mask and the detector array’s long axis were along the major axis of the disk at PA=60◦. The standard chop-nod technique was applied with a 15′′ chop throw in the North-South direction. For spectropolarimetry, we positioned the 1′′.04×2′′.08 slit with the slit’s longer axis oriented at 72◦ from the North to cover the brightest part of the disk. The chop throw was set at 8′′ in the North-South direction. The data were reduced using custom IDL software, as described in Berry & Gledhill (2014); Li (2014). We computed normalized Stokes parameters q and u, where q=Q/I and √ 2 2 − 2 u=U/I. The degree of polarization p = q + u σp, where the last term (the “debias” term) was introduced to remove the positive offset in the signal floor resulting from noise. The polarization PA is computed as θ = 0.5arctan(u/q). The uncertainties σq and σu associated with the normalized Stokes parameters were derived using a 3-sigma clipping algorithm (Robinson, 1987), and they were then propagated through the analysis to obtain the polarization uncertainty σp and polarization PA uncertainty σθ = σp/2. The total intensity (Stokes I ) images of WL 16 in the three passbands are presented in Fig. 3-1. The polarization image at 8.7 µm, where the highest signal-to-noise ratio is achieved, is shown in Fig. 3-2. Polarizations in the Si2 (8.7 µm) and Si4 (10.3 µm) filters measured with an aperture of 1′′ in radius centered on the star are given in Table 3-3. The Si6 (12.5 µm) data are too noisy to provide meaningful polarization information. Note that in the polarimetric images, the upper and lower edges of the WL 16 disk are truncated by the focal mask, and the polarization vectors near the edges are unreliable and thus not displayed in Fig. 3-2. The data presented in Fig. 3-2 have been smoothed by 5×5 ′′ ′′ pixel (0.4 × 0.4) binning. Vectors are only plotted where p/σp ⩾ 2.0 and p ⩽ 6%. For spectropolarimetry, we extracted one-dimensional spectra by integrating the central 21 pixels (1′′.6) along the slit. We fitted a second-order polynomial of identified skylines to calibrate the wavelength. We rebinned the raw polarization data (o and e ray spectra) into 0.1 µm wavelength (5 pixels) bins, i.e., downsampling, and then applied an additional 5-pixel (0.1 33 Table 3-2. Observing Log Imaging Filters ∆λ Date Integration Sensitivity FWHM (PSF) (µm) (µm) (UT) Time (s) (mJy/10 σ/1 h) (′′) Si2(8.7) 1.1 2013 Aug 6 946 8.3 0.50 Si4(10.3) 0.9 2013 Aug 6 873 10.8 0.60 Si6(12.5) 0.7 2013 Aug 6 952 20.7 0.40 Spectropolarimetry Source(s) Date Integration (UT) Time (s) WL 16 2015 Jul 4&5 2648 Elias 29 2015 Jul 4&5 530 µm) boxcar-smoothing to the data to further increase the signal-to-noise ratio resulting in the equivalent spectral resolution R≈25. We masked out the region between 9.4–10.0 µm, which was dominated by the atmospheric ozone features. For the total intensity spectrum presented in Fig. 3-3 , we smoothed the raw data with a 3-pixel (0.06 µm) boxcar and the resultant spectral resolution R≈50. Polarimetric spectra of WL 16 and Elias 29 are shown in Figs. 3-4 and 3-5. 3.3 Orign of the Polarization Our observations (Fig. 3-1) confirm that WL 16 is well-resolved and extended in the mid-IR, which was previously ascribed to the emission from PAHs and VSGs (RB03). Our spectrum (Stokes I ) of WL 16 indeed shows PAH emission features at 8.6, 11.2, 12.0, and 12.7 µm (Fig. 3-3) that dominate the spectrum. Most pertinent to our primary focus on B-fields in disks is our finding that the 8.7 µm polarization vectors in WL 16 are roughly uniform in both magnitude and orientation across most of the field of view. This rough uniformity implies that the observed polarization results from grain alignment in a correspondingly uniform B-field. 3.3.1 Extinction toward WL 16 The interstellar extinction toward WL 16 is relatively high, and our first task before characterizing the B-field giving rise to the polarization is to determine the relative roles of the polarization arising in the immediate disk environment and in the general intervening interstellar medium. The 2MASS stellar extinction map implies AV = 26.40.6 mag at the location of WL 16 (Lombardi et al., 2008). It is not trivial to deduce the interstellar extinction 34 Figure 3-1. Total intensity maps of WL 16 at 8.7, 10.3, and 12.5 µm from top to bottom. Contours are surface brightness and logarithmically spaced. Parts of the disk structure are truncated by the CanariCam mask. Some peculiar features of the disk, including a spiral-arm-like structure (green dots) and a dark lane in the SW disk, which may indicate a disk warp, are highlighted. 35 Figure 3-2. The 8.7- µm linear polarization map of WL 16 superimposed on (total intensity) contours. The polarization vectors are plotted at the center of each 5 × 5 binned ′′ ′′ pixels, corresponding to 0.4 × 0.4 (50 × 50 AU), and are only plotted where p/σp ⩾ 2.0 and p ⩽ 6%. The lengths of polarization vectors are scaled to the polarization percentage and orientations correspond to the polarization PA. The green vector at the center, which is oriented at 30◦from North, is the polarization PA observed in the optical (Goodman et al., 1990) and near-IR (Sato et al., 1988) bands. for young stars, since they can also exhibit optical and near-IR excess emissions from their accreting disks. To help disentangle the interstellar extinction from the properties intrinsic to WL 16, we consider the two-color diagram as a tool to estimate the foreground interstellar extinction. Meyer et al. (1997) combined the observed near-IR properties of T Tauri stars and accretion disk models, and found that the dereddened colors of T Tauri stars occupy a narrow band, called the locus of classical T Tauri stars (CTTS locus), in the two-color (JHKL) diagram. While recognizing that the CTTS locus is derived from T Tauri stars, and Herbig stars may have different properties, we nevertheless apply this method to WL 16 and Elias 29, taking their near-IR photometry (Table 3-4) from Cutri et al. (2003) and plot them on the JHK two-color diagram in Fig. 3-6. We decompose the total reddening into a component of interstellar extinction and a component of intrinsic reddening. We use the following relation to derive the corrected J-band interstellar extinction for objects with excess in the near-IR (Gorlova et al., 2010): 36 Figure 3-3. The low-resolution (R≈50) spectrum of the brightest central 1′′.6 (21 pixels) region of WL 16. The slit (white rectangle) is shown on the inset image of WL 16 at 8.7 µm. The raw data were smoothed with a boxcar of 3 pixels (0.06 µm) in width. PAH emission features are seen at 8.6, 11.2, 12.0, and 12.7 µm. The 8.6 µm feature originates from C-H in-plane bending. The 11.2, 12.0, and 12.7 µm features originate from C-H out-of-plane bending. A /E(J − H) A = J J,corr k −1 − E(H − K) E(J − H) {CTTS / } (J − H) − (J − H)0 − [(H − K) − (H − K)0] (3–1) kCTTS where AJ is the total extinction, E(J-H) and E(H-K) are the color excesses, H-K and J-H are the observed colors, (H-K)0 and (J-H)0 are the intrinsic colors depending on the spectral 37 Figure 3-4. Polarization and emission/absorption decomposition results of WL 16. Black dots with 1 σ error bars are measurements from spectropolarimetry. The atmospheric ozone band is masked out. The data are integrated across the central 1′′.6 (21 pixels) region. We downsample the original data over 0.1 µm wavelength (5 pixels) bins and then apply a 5-pixel boxcar smoothing to increase the signal-to-noise ratio. The resultant spectral resolution is about 25. Overlayed red points with 1 σ-error are measurements from polarimetric imaging at 8.7 and 10.3 µm. The dashed and the dotted lines are the best fitting absorptive and emissive components, and the solid lines are the combination. From upper left to the bottom right: a) Normalized Stokes parameter q (Q/I). (b) Stokes u (U/I). (c) Polarization degree p. (d) Position Angle (PA) θ. θ is approximately 33◦ regardless of the wavelength. The decomposition fitting indicates that it is absorptive-polarization dominant. 38 Figure 3-5. Polarization and emission/absorption decomposition results of Elias 29. Black dots with 1 σ error bars are measurements from spectropolarimetry. The atmospheric ozone band is masked out. The data are integrated across the central 1′′.6 (21 pixels) region. We apply a 5-pixel boxcar smoothing to increase the signal-to-noise ratio. The resultant spectral resolution is about 48. Overlayed red points with 1 σ-error are measurements from polarimetric imaging at 8.7,10.3, and 12.5 µm. The dashed and the dotted lines are the best fitting absorptive and emissive components, and the solid lines are the combination. From upper left to the bottom right: a) Normalized Stokes parameter q (Q/I). (b) Stokes u (U/I). (c) Polarization degree p. (d) Position Angle (PA) θ. θ is approximately 20◦ regardless of the wavelength. The decomposition fitting indicates that it is absorptive-polarization dominant. 39 type of stars (Kenyon & Hartmann, 1995), and kCTTS = 0.58 0.11 is the slope of CTTS locus for JHK. Adopting the near-IR extinction relations AJ /AK =2.840.46, E(J-H)/E(H- K)=1.770.13 (Wang & Jiang, 2014) and the Ophiuchus extinction law (Chapman et al., 2009) to convert near-IR extinction to visual and mid-IR extinction, we derive a value for the interstellar extinction (AV =28 mag) toward the central star. Thus, the interstellar material accounts for most of the extinction along our line of sight. We derive the extinction toward WL 16 at 10.3 µm, A10.3, to be 1.970.30 mag. The interstellar visual extinction of Elias 29 is about 43 mag and A10.3 is 3.00.37 mag using the same approach (Table 3-4). 3.3.2 Decomposition of Absorptive and Emissive Polarization We have shown that the extinction due to the foreground interstellar medium is significant, and therefore it may also account for much of the observed mid-IR polarization. However, the dust inside the disk can potentially emit polarized thermal emission that contributes to the net polarization as we see toward WL 16. To evaluation the contribution from emissive polarization, we adopted the method present in Aitken et al. (2004), which showed that, in the mid-IR, the emissive and absorptive polarization components can be identified, and separated from each other, using spectral polarimetry. Aitken’s method is based on their finding that the emissive and absorptive polarization profiles (pem(λ) and pabs(λ)) across the 10- µm silicate feature correlate with each other through the extinction, i.e., pem(λ) = pabs(λ)/τ(λ), where τ(λ) is the extinction (due to silicate) near 10 µm. In practice, the absorptive profile, pabs(λ), is taken as that of the Becklin-Neugebauer (BN) object, which contains absorptive polarization alone, and the extinction curve, τ(λ), is derived from the observations of the Trapezium region of Orion. Aitken’s method has been successfully used by a number of studies to separate emissive and absorptive polarization components from the observed total polarization when spectral polarimetry, or imaging polarimetry at multiple wavelengths, is available (e.g., Smith et al. 2000; Barnes et al. 2015). We apply this method to the polarization data of WL 16 and Elias 29, with the results shown in Figs. 3-4 and 3-5, respectively. Normalized Stokes parameters q and u, polarization 40 Figure 3-6. The JHK color-color diagram. The green dots are the loci of dwarfs (Kenyon & Hartmann, 1995), which represent the intrinsic color of stars of different spectral types. The black dashed line is the interstellar reddening vector (Kenyon et al., 1998), with plus symbols indicating AV =10 mag intervals along the vector. The magenta dot-dashed line is the CTTS locus (Meyer et al., 1997). The black dot and blue triangle with 1 σ errorbars are the near-IR photometry measurements of WL 16 and Elias 29 (Cutri et al., 2003). We decompose the total extinction to the sum of an interstellar extinction vector and an intrinsic reddening vector, as shown by the dotted lines. 41 Figure 3-7. Comparison of the polarization profiles of WL 16 (black) and Elias 29 (blue). Both stars are located in/behind the same molecular cloud. At 10.3 µm, where the polarization peaks, we find their fractional polarizations are proportional to their interstellar extinctions. fraction p, and polarization position angle for WL 16 and Elias 29 are shown in each plot. The position angles θ for both WL 16 and Elias 29 are almost constant with wavelength from 8.0 to 13.0 µm, indicating that a single mechanism probably produces the polarization (Efstathiou et al., 1997). Should there be multiple components, we would usually expect θ to change with wavelengths (i.e., unless the absorbing and emitting fields are at 90◦). In addition, the value of the polarization p peaks around 10.3 µm, which is expected from the absorptive polarization of astronomical-silicate feature in this spectral region. The decomposition indicates that an absorptive-dominant polarization profile provides an excellent fit for both 42 Table 3-3. Polarization Measurements of WL 16 λ Flux Density p θ ( µm) (Jy) (%) (◦) 8.7 5.32(0.50) 0.91(0.12) 27.09(3.71) 10.3 1.68(0.17) 2.05(0.48) 29.02(6.70) 12.5 3.30(0.30) .. .. Values in parentheses are 1 σ uncertainties of measurements. All position angles are calibrated East from North. stars, with negligible emissive components (less than 0.25% for both cases). This result strengthens our conclusion that the mid-IR polarization of WL 16 is primarily due to absorption by interstellar silicate dust. Since absorption by elongated interstellar particles produces polarization parallel to the projected B-field, the direction of projected foreground magnetic field is about 33◦4◦(measured from spectral polarimetry and averaged over the entire 10- µm band). The observed polarizations at 10.3 µm for both WL 16 and Elias 29 are proportional to their respective values of interstellar extinctions as shown in Fig. 3-7. From this ratio, we derive a value for the polarization efficiency, which is defined as the ratio of polarization percentage −1 to the extinction (Cashman & Clemens, 2014), (p10.3/A10.3) ≃ 1% mag in Table 3-4. The parameter can be useful for interpreting the polarization properties of other sources and understanding the dust alignment efficiency, B-field strength, and alignment mechanisms in dense molecular clouds. This value is within the 0–3% range derived in Smith et al. (2000). We note however that the observed polarization corresponds only to the B-field component projected on the sky (integrated along the line of sight), whereas the extinction depends on the total dust column density along the line of sight. Since there is no reason to assume that the projected B-field is the same everywhere, we recognize that the diagnostic power of this value of the ratio (p10.3/A10.3) beyond the immediate environment of WL 16 may be limited. 3.3.3 Relationship of WL 16 and Regional Magnetic Fields Goodman et al. (1990) probed the large-scale B-field morphology of the Ophiuchus dark cloud complex at optical wavelengths. The distribution of polarization of background starlight 43 Table 3-4. Extinction and Polarization of WL 16 and Elias 29 WL 16 Elias 29 Ja 14.164 (0.029) 16.788 (0.178) Ha 10.478 (0.023) 11.049 (0.044) K a 8.064 (0.016) 7.140 (0.021) A10.3 1.97 (0.30) mag 3.00 (0.37) mag p10.3 2.00 (0.24)% 3.21 (0.07)% −1 −1 p10.3/A10.3 1.02 (0.20) % mag 1.07(0.13) % mag a(Cutri et al., 2003), in units of magnitude Values in parentheses are 1 σ uncertainties of measurements. across this region is fairly smooth, and our polarization measurements agree well with this trend of the large-scale field (Fig. 3-2). In addition, the K-band polarimetry of WL 16 by Sato et al. (1988) and Beckford et al. (2008), indicates a fractional polarization p=4.870.23% with θ=27.01.4◦ and p=5.100.05% with θ=33.900.03◦, respectively. These values of θ are consistent with our measurements across 8–13 µm, thus strengthening our conclusion that our polarization measurements do not constrain the B-field inside the disk, but instead, trace the foreground B-field. 3.3.4 Intrinsic Polarization While the observed polarization of WL 16 is mainly from the foreground, and the spectropolarimetry decomposition suggests a very low emissive polarization value (less than 0.25%), can we still place a limit on the intrinsic polarization from the disk? The emissive polarization percentage detected in the case of AB Aur, an archetypal Herbig Ae star (Li et al., 2016), is 0.5%. Assuming the foreground polarization is 2% with a fixed position angle, a 0.5% emissive polarization component could result in, at most, 7◦deviation from the assumed foreground polarization orientation. The uncertainties of our measurements (Table 3-3 and Fig. 3-4) prevent us extracting from the total observed polarization any intrinsic emissive polarization component if it is less than 0.5%. Reiterating the limitation that we are only sensitive to the projected B-field, we place an upper limit of 0.5% on the intrinsic polarization in WL 16. The intrinsic emissive mid-IR polarization in WL 16 is therefore likely to be much lower than a few percent predicted by Cho & Lazarian (2007). Polarization 44 from protoplanetary disks depends on the dust properties, e.g., dust sizes, oblateness, and composition, and also the entanglement of magnetic fields with dust grains, e.g., dust alignment and the strength of magnetic fields (Hughes et al., 2013). In the case of WL 16, it may result from a combination of these factors. 3.4 Discussions 3.4.1 PAHs in WL 16 WL 16 is rich in PAH emission features, making it useful for studying PAH properties in the environments of Herbig stars (Geers et al., 2007; DeVito & Hayward, 1998). ISO/SWS and Spitzer/IRS spectra of WL 16 show emission features at 6.2, 7.7, 8.6, 11.2, 12.7, 16.4, and 17.0 µm identified with de-excitation via C–C or C–H vibrational transitions of the UV-excited PAHs (Leger & Puget, 1984; Draine & Li, 2001). The 3.3 µm PAH feature is not detected (Geers et al., 2007). Because a PAH molecule can be excited by a single UV photon, PAHs can trace the disk emission up to large distances from the star. The intensities of the C–C stretching and C–H in-plane bending modes, which fall in the 6–9 µm range, are generally stronger for ionized PAHs than PAH neutrals. We use the total intensity spectrum (Stokes I ) from spectropolarimetry observations to derive the ratio of the 8.6 to 11.2 µm surface brightness profiles, I8.6 and I11.2 (band-width 0.33 µm), in order to trace the charge state of PAHs along the disk major axis (Joblin et al., 1996). The PAH features are often blended with neighboring PAH features (e.g., the 8.6 µm feature is on the wing of 7.7 µm feature). To correct the blended PAH emission features and subtract the continuum, we adopt the package PAHFIT by Smith et al. (2007) to derive the power emitted in each PAH features. There is as much as a 10% brightness difference in the PAH features and continuum emission between the integrated NE and SW side disk flux. Figure 3-8a shows the derived spatial distribution of PAHs at 8.6 and 11.2 µm in the central 2′′.0 region along the major axis of the disk. Although the spatial behaviors of 8.6 and 11.2 µm are slightly different, the SW disk intensity is much lower than that of the NE disk. Moreover, I8.6/I11.2 reaches the maximum (about 2.0) at the central region and decreases farther from the star (Fig. 3-8b), which means there are more 45 ionized PAHs near the star, with neutral PAHs dominating the external region, as expected (Weingartner & Draine, 2001). Interestingly, we also notice (Fig. 3-8b) that, while the values of I8.6/I11.2 on the two sides of the disk are roughly comparable within 0′′.3 (38 AU) of the star, they are markedly ′′ different between the two sides of the disk beyond about 0.5 (63 AU), with I8.6/I11.2 being higher on the SW side than that of the NE side in these outer regions as shown in Fig. 8c. If the disk at this large radius is optical thin in the mid-IR, the asymmetry can be explained by the neutralization of PAHs through electron recombination: PAH cations recombine more effectively with electrons in the NE than they do in the SW (Li & Lunine, 2003). However, WL 16, like most young Herbig-star disks (e.g., AB Aur, Mari˜naset al. (2006)), probably has an optically thick disk, the mid-IR emission that we see most likely arises in a thin layer near the disk surface (the disk ‘atmosphere’) (Testi et al., 2014). Thus, the interpretation of the anti-correlation is not obvious. Complicated disk structures, such as an asymmetric puffed-up disk inner rim and disk warps, combined with the inclined viewing angle of disk, may provide explanations. To solve the puzzle, multi-wavelengths observations are necessary. We do not see any polarization features at PAH emission bands (e.g., 8.6 and 11.2 µm) in the case of WL 16 (Fig. 3-4). Analytical modeling does predict detectable PAH polarization in astrophysical environments (Sironi & Draine, 2009). Though the astrophysical conditions they considered may be different from the condition in protoplanetary disks, PAH polarization is too small (0.1–0.5%) to be distinguished from the contribution of linear dichroism by aligned foreground dust. 3.4.2 Disk Morphology The total intensity maps (Fig. 3-1) of WL 16 reveal intriguing asymmetric features in the disk, including: (1) an S-shaped, spiral-arm-like structure extending to both sides of the disk (green dots in Fig. 3-1); (2) a dark lane immediately outside the spiral-arm-like structure in the SW disk, at 1′′.3 (163 AU) from the star; (3) twisted surface brightness contours, i.e., position angles of the contours’ major axes are changing from 95◦ (at smaller radii) to 50◦ 46 3.0 2.0 NE SW 1.8 NE SW 2.5 11.2 /F 1.6 8.6 R=F 2.0 1.4 (b) 1.2 1.5 1.0 0.5 0.0 -0.5 -1.0 offset (arsec) 1.2 Flux density (Jy) 1.0 1.1 NE /R 0.5 SW 8.6 µm 1.0 R (a) 11.2 µm (c) 0.0 0.9 1.0 0.5 0.0 -0.5 -1.0 0.0 0.2 0.4 0.6 0.8 1.0 offset (arsec) offset (arsec) Figure 3-8. (a): The inner 2′′.0 spatial distribution of intensity at 8.6 and 11.2 µm along the major axis of the disk. There is a clear brightness asymmetry of SW (blue) and NE (red) sides of the disk at both wavelengths. (b): R=I8.6/I11.2 along the major axis of the disk as a tracer of the charge state of PAHs. The larger the value, the more ionized are the PAHs. The ratio is decreasing farther from the central star as the decreasing UV radiation. (c): The relative RSW and RNE . There are more ionized PAHs at the SW side of the disk than the NE. (at larger radii); and (4) an asymmetric brightness distribution, the NE side being significantly brighter than the SW side of the disk. Though especially obvious in the PAH bands, some of these features are also seen in the continuum. In general, they appear to be qualitatively consistent with the scenario of an warped inner disk with respect to the outer disk as the result of precession induced by an unseen planet or stellar companions (Dong et al., 2015b). In this scenario, the dark lane in the SW could be a shadow cast by the disk warp. The illumination of a warped inner disk can mimic spiral features (Quillen, 2006). Although the twisted contours could be a result of a stellar companion of WL 16, a search for such a companion turned out 47 to be unsuccessful (Barsony et al., 2003). At 40′′ away, Elias 29 is the nearest (in projection) star of WL 16, but there is no obvious evidence that these two objects are dynamically related. We expect high-resolution and high-sensitivity radio telescopes, such as ALMA and IRAM, could lead to our new understanding of the dynamic and structure of the object. 3.5 Summary We present mid-IR polarimetric imaging and spectral observations of WL 16 obtained with GTC/CanariCam. WL 16 has a well resolved, extended disk (diameter ∼ 900 AU) in the mid-IR with ∼ 2% polarization. Our main conclusions are as follows: 1. Spectropolarimetry of WL 16 firmly supports the hypothesis that the observed mid-IR polarization is dominated by absorptive polarization arising from aligned non-spherical dust grains in the foreground. Polarized emission from dust inside the disk is non-measurable with an upper limit 0.5%. Because, in the most widely accepted dust alignment mechanisms, the absorptive polarization is parallel to the direction of B-field, we interpret the observed polarization map as indicating that the B-field in the molecular cloud is fairly uniform with projected orientation of about 33◦ from North. The direction is consistent with the near-IR polarization at WL 16 and the large-scale optical ρ Ophiuchus star formation region polarimetry. Though our original goal was to probe the B-field inside the protoplanetary disk, our study shows the importance of characterizing the foreground polarization as well. 2. The maximum values of the polarization of WL 16 and the nearby-polarized standard Elias 29 are proportional to their interstellar extinction. Using this ratio, we are able to characterize the polarization efficiency of dust grains in the dense molecular cloud, −1 (p10.3/A10.3) ≃1% mag . Keeping in mind that the observed polarization is associated with the projected B-field component, the parameter may be useful for constraining the dust alignment efficiency and properties in this region and for interpreting the observed polarization to other sources. 3. WL 16 is rich in PAH emission features and we have detected the 8.6, 11.2, 12.0, and 12.7 µm features in its disk. We see an asymmetry in the ratio I8.6/I11.2 between the two sides of the disk, with the NE (SW) side being brighter (fainter) at both 8.6 and 11.2 µm but with a lower (higher) value of I8.6/I11.2. This anti-correlation may be explained by complicated disk structures, e.g., warps and asymmetric disk inner rim. 4. The total intensity maps of WL 16 reveal asymmetric features, such as the S-shaped spiral-arm-like structure, twisted contours, asymmetric brightness distributions, and a dark lane on the SW side of the disk. These may indicate a disk warp, with future observations, especially with ALMA and IRAM, likely to fully clarify this picture. 48 This chapter, with minor differences, was published in its entirety under the same title in Monthly Notices of the Royal Astronomical Society, Volume 465, Issue 3, p.2983-2990 49 CHAPTER 4 DETECTION OF POLARIZED INFRARED EMISSION BY POLYCYCLIC AROMATIC HYDROCARBONS IN THE MWC 1080 NEBULA 4.1 Introduction A distinctive series of infrared (IR) emission features generally attributed to the stretching and bending vibrational modes of planar polycyclic aromatic hydrocarbon (PAH) molecules are observed in most dusty astronomical objects at 3.3, 6.2, 7.7, 8.6, 11.3, and 12.7 µm (Leger & Puget, 1984; Allamandola et al., 1985). Leger (1988, hereafter L88) first noted that these IR emission features, if due to PAHs, are expected to be polarized as a result of anisotropic illumination by a source of ultraviolet (UV) photons (e.g., stars). UV absorption is favored when the molecular plane faces the illuminating source. If the spinning of IR-emitting PAHs do not deviate significantly from their initial orientations at UV absorption, the IR emission will preferentially come from those PAHs that are facing the illuminating source at the time of UV absorption. Therefore, PAHs will emit polarized light, with the polarization direction of out-of-plane modes (11.3, 12.7 µm) being along the source-molecule direction and that of in-plane modes (3.3, 6.2, 7.7, 8.6 µm) being perpendicular to it. Sironi & Draine (2009, hereafter SD09) revisited the L88 scheme by considering realistic rotational dynamics of PAHs as well as an arbitrary degree of internal alignment between the grain symmetry axis and its angular momentum. Using realistic estimates of rotational temperatures for a typical PAH molecule of 200 carbon atoms, SD09 derived a value for the polarization fraction of 0.53% for the 11.3 µm feature in the case of the Orion Bar. Sellgren et al. (1988) performed the first systematic search for the polarization of the 3.3 and 11.3 µm PAH features in a number of astronomical sources, obtaining upper limits of 1% and 3% for the 3.3 and 11.3 µm features, respectively, in the Orion Bar. To test the PAH identification of these IR features and gain insight into the alignment of PAHs, we searched for linearly polarized PAH emission in the nebula associated with MWC 1080, a stellar cluster located at a distance of 2.2 kpc (Wang et al., 2008). The primary 4 star, MWC 1080A, is classified as a B0e star and has a luminosity of L/L⊙ ≈ 10 , which, 50 together with its stellar companions, illuminates the surrounding gas and dust in the adjacent molecular cloud. The mid-IR image of MWC 1080 at 11.2 µm resembles a pinwheel, with opposing arms curving off to the northwest and southeast (Fig. 4-1). The extended mid-IR emission resembling spiral arms or wings around MWC 1080 actually traces the internal surfaces of a biconical cavity created by the outflow from MWC 1080A (Li et al., 2014a; Sakon et al., 2007). The brightest part of the nebula (green rectangle in Fig. 4-1) lies 0.03 pc in projection to the northwest (hereafter NW nebula) of MWC1080A and extends ∼0.1 pc from the northeast to the southwest. That may well provide an optimal, Orion-Bar-like viewing geometry, i.e., a long column density along the line of sight through the photodissociation region and at an angle between the line of sight and illumination direction (α≈90◦) that is almost ideal for observing maximum polarization (SD09). Anomalous microwave emission (AME) is often ascribed to the rotational emission from rapidly-rotating nanoparticles (Draine & Lazarian, 1998a,b). PAHs are often considered to be associated with the AME due to their abundances and small sizes (Kogut et al., 1996; Leitch et al., 1997), although there now appears to be some doubt about the hypothesis (Hensley et al., 2016). Nitrogen-substituted PAHs (i.e., nitrogen in place of carbon), with greater dipole moments, may also be important components of the carriers of the AME (Hudgins et al., 2005; Mattioda et al., 2008). Therefore, our observations have broader implications for determining the alignment and polarization of rapidly spinning PAHs, which bears on the quest for the cosmic microwave background (CMB) B-mode (Hoang et al., 2013). 4.2 Observation and Data Reduction CanariCam is the mid-IR (8-25 µm) multi-mode facility spectrometer and camera on the 10.4 m Gran Telescopio CANARIAS (GTC) in La Palma, Spain (Telesco et al., 2003). It employs a 320×240-pixel Raytheon detector array with a pixel scale of 0′′.079, which provides a field of view of 26′′×19′′with Nyquist sampling (two pixels per λ/D) of the diffraction-limited point-spread function at 8 µm. Polarimetry is accomplished through the use of a Wollaston prism and a half-wave plate rotated to angles of 0◦, 22.5◦, 45◦, and 67.5◦. A Wollaston prism 51 in the optical path divides incoming light into two beams (ordinary and extraordinary), which are recorded by the detector simultaneously. We obtained low-resolution (R≡λ/∆λ≈50) spectropolarimetry observations of the NW nebula of MWC 1080 on 2015 July 31, August 5, and August 7 spanning the wavelength range 7.5–13.0 µm. We made four separate measurements, as indicated in the observation log presented in Table 4-1. The spectroscopic observations of the NW nebula were interlaced with observations of the Cohen standard star HD 21330 (Cohen et al., 1999) for flux and point-spread-function (PSF) calibration, and the standard star AFGL 2591, selected from Smith et al. (2000) to calibrate the polarization position angle. The standard mid-IR chop-nod technique was applied with an 11′′ chop throw in the northwest-southeast direction. We positioned the 1′′.04×2′′.08 slit with the slit’s longer axis oriented at 45◦ from the North to intersect the brightest part of the nebula, enclosed by the green rectangle in Fig. 4-1. The 11.2 µm image of MWC 1080 is adopted from Li et al. (2014a), and the data were taken using Michelle, the facility mid-IR camera at Gemini North. The data were reduced using custom IDL software, as described in Berry & Gledhill (2014) and Li (2014). We extracted one-dimensional spectra by integrating the central 21 pixels (1′′.6) along the slit direction. Wavelength calibration was done using several sky lines identified in the raw images. We computed normalized Stokes parameters q (q=Q/I ) and u (u=U/I ) using both the difference and ratio methods, with each providing consistent results (Tinbergen, 2005). The data were calibrated for instrumental efficiency and polarization. The estimated instrumental polarization was 0.6% as measured with HD 21330 and was subtracted from the observations of the NW nebula in the Q − U plane. The degree of √ polarization p = q2 + u2 − σ2, where the last term (the ‘debias’ term) was introduced to remove a positive offset in the signal floor resulting from squared background noise. Debiasing may introduce negative values if the noise fluctuations are stronger than the signal. The polarization position angle was computed as θ = 0.5arctan(u/q). The uncertainties σq and σu associated with the normalized Stokes parameters were derived using a standard 3-sigma 52 clipping algorithm (Robinson, 1987), and were then propagated through the analysis to obtain the polarization uncertainty σp and position angle uncertainty σθ = σp/2p (Patat & Romaniello, 2006). All the polarization position angles were calibrated east from north. We masked out the region between 9.2–10.0 µm, which is dominated by the atmospheric ozone feature. The three-pixel (0.06 µm) boxcar-smoothed intensity spectrum of the nebula (Stokes I) is presented in Fig. 4-2a. To further increase the signal-to-noise ratio (S/N) of the polarization measurements, we rebinned the ordinary and extraordinary ray spectra into 0.12 µm wavelength (6 pixels) bins (downsampling) and then applied an additional three-pixel (0.36 µm) boxcar-smoothing to the data. That results in an equivalent spectral resolution of the polarization spectrum of R≈32 (Fig. 4-2b). Stokes u and q are plotted in Figs. 4-3b and c, respectively. To emphasize the statistical significance of the measurements, we plot the S/N of the polarized intensity in Fig. 4-3a. 4.3 Results We present in Fig. 4-2 mid-IR intensity and polarization spectra of the NW nebula covering the 8.0–13.0 µm wavelength range. The PAH emission features are clearly seen in Fig. 4-2a, including the in-plane C–H bending feature at 8.6 µm and the out-of-plane C–H bending features at 11.3 (solo-CH), 12.0 (duet-CH), and 12.7 µm (trio-CH) (Allamandola et al., 1989). The relative strengths of the PAH features depend on the size, structure, and charging of PAHs (Allamandola et al., 1999; Draine & Li, 2007), and the physical conditions of the region where they are found (Bakes & Tielens 1994; Weingartner & Draine 2001). By fitting the spectrum of MWC 1080 obtained with the Infrared Space Observatory (ISO), which exhibits a more complete set of PAH emission features but mixed emission from both stars and the nebula, Seok & Li (2017) determine that the best fit of PAHs in this environment are mostly neutral and large. Exposed to the energetic photons from a B0e star with an effective temperature of ∼30,000 K (Wang et al., 2008), small PAHs are probably unstable, and the only PAHs that survive are ones that are large and cata-condensed in structure. 53 The polarization spectrum shown together with 1σ error bars of the NW nebula is presented in Fig. 4-2b. Two spectral regions show significant polarization. One is in the range ◦ ◦ of 10.9–11.7 µm, with p11.3=1.90.2% and position angle of 77.2 3.2 , the latter indicated by the solid red lines superimposed on the intensity map of MWC 1080 in Fig. 4-1. There is good consistency among the four separate measurements listed in Table 4-1. The other significant polarization is at 10.0–10.7 µm, peaking around 10.3 µm, and with p10.3=5.41.6% and position angle of 46.7◦8.2◦. While the 10.3 µm polarization feature is seemingly higher than that at 11.3 µm, the statistical significance of the results needs to be considered. We note here that the detection of polarization near 10.3 µm is only marginally significant, being barely 3σ, whereas the 11.3 µm detection is robust, being around 9σ (Fig. 4-3a). The 3σ upper limits for the polarization at 8.6, 12.0 and 12.7 µm are 1.9%, 2.0%, and 2.8%, respectively (Fig. 4-2b). While, theoretically, the emission at longer wavelengths is expected to have higher polarization (SD09), the 8.6, 12.0 and 12.7 µm PAH emission features are much weaker than the 11.3 µm feature and the lack of detection of the polarization in these regions is not surprising. Examining Stokes u and q explicitly (Fig. 4-3) demonstrates that the 10.3 µm polarization is apparent mainly in Stokes u, while the 11.3 µm polarization is mainly in Stokes q. This difference implies a different origin for the 10.3 µm if real and 11.3 µm features, almost certainly indicating that they originate in different dust populations. 4.4 Discussion The polarization feature between 10.9 and 11.7 µm is well correlated in wavelength position with the 11.3 µm PAH emission feature (Fig. 4-3), and we conclude that the polarized emission is indeed due to PAH molecules. The peak value p11.3 ≃ 2% is consistent with the previous PAH polarization search by Sellgren et al. (1988) who established an upper limit of 3% on its polarization in the Orion Bar. The polarization position angle of the feature has an offset angle of ∼60◦from the position angle of the projected illumination direction from the star to the nebula (∼315◦), as shown in Fig.4-1. 54 4.4.1 Numerical Calculations using SD09 We adopt the up-to-date SD09 models to model the degree of polarization in the environment of a reflection nebula. SD09 calculated the polarization for randomly oriented PAHs. They considered the intramolecular vibrational-rotation energy transfer (IVRET) process, which allows the efficient energy exchange between rotation and vibration. SD09 model the polarization of PAH emission by introducing two parameters, γ0 = Trot/T0 and γir = Trot/Tir, where Trot is the rotational temperature determined by gas-grain interactions, photon absorption and IR emission, T0 is the vibrational grain temperature that determines internal alignment prior to UV absorption, and Tir is the temperature during IR emission. Alignment by an external magnetic field is ignored in their model, and the polarization degree of PAH emission features is essentially determined by two parameters, γir and γ0. Since the internal alignment temperature Tia ≡ Tir is fixed for the different PAHs, and T0 is determined by the radiation field, the polarization degree is determined by Trot. Consequently, the knowledge of Trot is critical for achieving realistic predictions for the polarization of PAH emission. Collisions with gas neutrals and ions, UV absorption and subsequent IR emission, and electric dipole emission contribute to the rotation of PAHs (see Draine & Lazarian 1998b; Hoang et al. 2010). We assume in a reflection nebula that the gas temperature Tgas=100 3 −3 3 K, gas density nH = 10 cm , and the radiation field parameter U = 10 . We consider the PAH geometry as in Draine & Lazarian (1998b) and a typical grain size of a=7.5 A˚ (200 carbon atoms) in which most of the PAH mass is concentrated as in the ISM (SD09). We then carry out simulations of PAH dynamics using the Langevin code (Hoang et al., 2010) and calculate the degree of polarization. However, the values of polarization obtained using the SD09 model with the consideration of different hydrogen ionization fractions xH, i.e., nH+ /nH (nH+ is the number density of ionized hydrogen), are less than 0.5% (Table 4.2), much smaller than our measured value of ∼2%. It appears that other mechanisms that can enhance the PAH alignment need to be considered. 55 4.4.2 Alignment with Magnetic Fields SD09 discussed two possibilities that can significantly enhance the polarization of PAH emission, including suprathermal rotation (i.e., the PAH rotational temperature Trot is much higher than the gas temperature Tgas) and perfect internal alignment (i.e., the molecule principal axis is aligned with the angular momentum during both UV absorption and IR emission). The former is unlikely, since there are no obvious physical processes that can spin-up nanoparticles to suprathermal rotation rates (Draine & Lazarian 1998b; Hoang et al. 2010). On the other hand, perfect internal alignment can only be achieved at very low dust temperatures (T0 of a few K) during the UV photon absorption, which requires very efficient energy exchange among the vibrational-rotational modes. This is also very unlikely due to quantum suppression that may occur in nanoparticles (Draine & Hensley 2016). Alternatively, we recognize that the enhanced polarization may arise from external alignment, i.e., the partial alignment of the angular momentum with an ambient magnetic field. Therefore, we have repeated the modeling incorporating the mechanism of resonance paramagnetic relaxation to align nanoparticles (Lazarian & Draine, 2000). Subject to an external magnetic field, protons in aromatic hydrocarbons in the laboratory are known to experience stronger shielding (or deshielding) effects than regular hydrocarbons, because the π-electrons are delocalized and are free to circulate. Astronomical magnetic fields can induce diamagnetic ring currents and polarizabilities in the π-electron clouds, resulting in coupling between the magnetic fields and two-dimensional PAHs, thereby forcing some alignment. In our models, the magnetic field strength is assumed to be B=100 µG(Crutcher, 2012). The calculation is described here briefly, with more details presented in Hoang (2017). ∥,⊥ We define uˆ-vˆ to be the plane of the sky (see Hoang 2017). We define Fu,v to be the in-plane (∥) and out-of-plane (⊥) emission by a PAH molecule with the electric field E in the ∥,⊥ ∥,⊥ uˆ-vˆ plane. Then Iu,v is the emission intensity from the PAH. The flux Fu,v depends on: 1) fLTE(θ, J), the probability distribution of the principal axis of the PAH plane being aligned with the grain angular momentum J (LTE stands for local thermal equilibrium); and 2) fJ (J), the 56 probability distribution of the angular momentum J being aligned with the direction of the magnetic field. The grain angular momentum J has spherical angles θ and ϕ. β is the nutation angle between J and the principal axis of the grain. The emission Iw, with w = (u, v), is obtained by integrating over the distribution functions ∫ ∫ π ∥,⊥ Iw (α) = fJ (J)dJ fLTE,0(θ0, J)dθ0 J∫ 0 π × ∥,⊥ fLTE,ir(θ, J)dθA⋆(β, θ0)Fw (β, ϕ, θ, α), (4–1) 0 where A⋆ is the cross-section of UV absorption as given in SD09. fLTE,0 and fLTE,ir describe the thermal fluctuations of the principal axis of PAHs before UV absorption and during IR emission. fLTE(θ, J) can be described by the Boltzmann distribution (Lazarian & Roberge, ∫ π 1997) with 0 fLTE(θ, J)sinθdθ = 1. To simplify the calculation and derive the maximum value of polarization, we assume that the magnetic field is parallel to the stellar incident radiation direction. We simulate the distribution of angular momentum from the Langevin equation assuming the ergodic system approximation and compute numerically the intensity of radiation using Equation (4–1). The resulting degree of polarization with the viewing angle α for in-plane and out-of-plane modes is ∥ ⊥ ∥ ⊥ I , (α) − I , (α) p∥,⊥( ) = u v , α ∥,⊥ ∥,⊥ (4–2) Iu (α) + Iv (α) We show the modeling results in Table 4.2 for the reflection nebula, which includes the ratio of the rotational temperature Trot and the gas temperature Tgas, the degree of alignment of the angular momentum with the magnetic field QJ, and the estimated polarization for the ◦ different hydrogen ionization fractions xH with a viewing angle α=90 . The polarization varies from 0.14% to 0.34% in model a, increasing to 0.87%–2.1% when the external alignment is taken into account as in model b. In Table 4.2, the degrees of polarization p computed for both models increase with the increasing hydrogen ionization fraction xH. The results suggest that the polarization of PAH emission is dominated by PAHs in regions with a higher fraction of hydrogen in ionic form. In regions with a higher xH (i.e., higher nH+ /nH), there will be more 57 electrons available to neutralize the PAH ions created by photoionization (Tielens, 2005). Indeed, as shown in Seok & Li (2017), the aromatic hydrocarbon emission features observed in MWC 1080 are best modeled in terms of a mixture of PAHs with ∼ 80% being neutral and ∼ 20% being ionized. As mentioned earlier, the delocalized π-electrons in neutral PAHs may play a crucial role in coupling neutral PAHs with the magnetic field. Partial alignment of PAHs with the magnetic field at a level of QJ ≃ 0.08–0.1 (the averaged degree of alignment of angular momentum J with B) is required to reproduce the observed ∼2% polarization fraction at 11.3 µm. This scenario also accounts for the observation that the polarization angle is offset from the illumination direction. When PAHs are aligned with the magnetic field, even though only partially, the polarization direction of the out-of-plane mode emission is expected to be along the magnetic field (Hoang 2017). 4.4.3 Relationship between Polarization Angles and the Ambient Magnetic Field Both L88 and SD09 predict that the polarization associated with emission arising from the out-of-plane vibration mode should be along the illumination direction, which contrasts with our observed polarization angle having a ∼60◦ offset. We explore one possible explanation, namely, PAH alignment by an external magnetic field (Table 4.2). If magnetic alignment is important, we expect the polarization direction to match that of the ambient magnetic field lines. Indeed, we find that the optical polarimetry observations of background stars within a few degrees on the sky from MWC1080 (serveral hundred parsecs in distance), indicates a fairly uniform optical polarization position angle of ∼80◦, which suggests a smooth interstellar magnetic field threads the whole region (Manset & Bastien 2001; Maheswar et al. 2002; Heiles 2000). The value of the position angle agrees well with our measured polarization position angle of 77.23.2◦ at 11.3 µm, which supports the hypothesis that PAHs are at least partially aligned with the ambient interstellar magnetic field threading the nebula and its neighborhood. Nevertheless, the emission and alignment of PAHs depend on local astrophysical conditions and the detailed properties of PAHs, especially their sizes. Nanoparticles with radii ≲10 A˚ are thought to be negligibly polarized with the greatest quantum suppression 58 of alignment (Draine & Hensley, 2016). Based on our results, it seems that other physical processes such as Faraday rotation braking that facilitate the alignment of nanoparticles need to be considered (Papoular, 2016), since it is evident that the starlight anisotropy scheme alone in L88 is not sufficient to explain the measured high level of polarization. It is also worth noting that, given their abundances and small sizes, emission by rapidly spinning PAHs is widely believed to be the origin of AME in the 10–60 GHz frequency range (Draine & Lazarian, 1998b,a). If true, the considerable alignment of PAHs as suggested by our detection, naturally produces polarized spinning dust emission for which the polarization level is proportional to the degree of alignment of the PAH angular momentum with the magnetic field at ∼GHz frequencies (Hoang et al., 2013). It implies that polarized emission from spinning Galactic-foreground PAHs can indeed constitute an obstacle to the detection of the CMB B-mode signal. 4.4.4 Marginally Detected 10.3 µm Polarization Feature As shown in Fig. 4-3a, we have a barely significant polarization detection at 10.3 µm (3σ). If real, the different behaviors of the 10.3 and 11.3 µm features in Stokes u and q (Fig. 4-3) suggest that they originate in different dust populations. Therefore, it does not affect our interpretation of the high S/N (9σ) polarization detection at 11.3 µm, our main focus of this work. We do note, however, that there is no distinct Stokes I spectral fingerprint coinciding with the 10.3 µm polarization. It is unlikely that the well-known silicate feature or one of its variants can account for this polarization, since the silicate absorptive polarization profiles are broad, spanning the entire 8–13 µm region (e.g., Henning & Stognienko 1993; Smith et al. 2000; Zhang et al. 2017) rather than relatively narrow and sharp as the feature we see here. Other possibilities, including nanoparticles with silicate (Tielens, 2013) or metallic Fe compositions, e.g., hygrogenated iron nanoparticles (Bilalbegovi´cet al., 2017), might be worth investigating (Hoang & Lazarian, 2016) if further observations confirm and gain insight into the feature. 59 Table 4-1. Observing Log a a UT Date Target RA DEC Integration Airmass p11.3 θ11.3 J2000 J2000 seconds % ◦ 2015 Jul 31 NW nebula 23 17 25.18 60 50 44.56 993 1.18–1.23 2.46(0.42) 80.2(4.8) 2015 Aug 5 NW nebula 23 17 25.18 60 50 44.56 993 1.18–1.20 1.48(0.44) 88.7(8.2) 2015 Aug 5 NW nebula 23 17 25.18 60 50 44.56 993 1.21–1.28 1.37(0.38) 72.0(7.7) 2015 Aug 7 NW bebula 23 17 25.18 60 50 44.56 993 1.20-1.27 1.95(0.35) 70.2(5.1) Table 4-2. Different models and polarization at 11.3 µm a b xH Trot/Tgas γir γ0 p(%) QJ p(%) 0.001 0.62 1.54 0.31 0.14 0.058 0.87 0.003 0.707 1.76 0.35 0.19 0.069 1.1 0.005 0.76 1.90 0.38 0.22 0.076 1.6 0.010 0.94 2.35 0.47 0.34 0.088 2.1 4.5 Summary We report the unambiguous detection of polarized PAH emission at 11.3 µm with a position angle of 77.23.2◦ and polarization degree of 1.90.2%, which confirms the PAH hypothesis that PAH molecules can indeed emit polarized light. The detection of polarization indicates that the alignment of PAHs is considerable. We find that the starlight anisotropy scheme alone is not sufficient to account for this polarization. The PAHs are at least partially aligned by the ambient magnetic field threading this young stellar region and its neighborhood, a conclusion strongly supported by the fact that the measured polarization angle is identical to the large-scale interstellar magnetic field spanning this region. This observation could have important consequences for the accurate estimate of Galactic foreground polarization, a consideration relevant to current goals to detect the CMB B-mode signal. We expect future polarimetry observations, e.g., with SOFIA/HAWC+ and GTC/CanariCam, covering the complete suite of PAH emission features (e.g., the 6.2 µm band dominated by PAH cations and 3.3 µm band by small PAHs) and various astrophysical environments, will deepen our understanding of the properties and alignment of PAHs, e.g., the effects of their sizes and charge states This chapter, with minor differences, was published in its entirety under the same title in The Astrophysical Journal, Volume 844, Issue 1, article id. 6, 7 pp. (2017). 60 N E 46'' DEC (J2000) DEC 41'' 60 ◦ 50'36'' 27.5s 26.5s 23h 17m 25.5s 24s RA (J2000) Figure 4-1. Intensity map (contours) of MWC 1080 system at 11.2 µm adopted from Li et al. (2014a). The slit (dashed rectangle) is positioned to enclose the brightest part of the reflection nebula in CanariCam/GTC spectropolarimetry observations. We integrate the polarization signal from the region enclosed by the dashed rectangle at 11.3 µm. The thicker segment, together with the two thinner segments, shows the derived polarization position angle with 1σ uncertainties of the 11.3 µm PAH emission feature. The position angle of the projected illumination direction from the star to the nebula is ∼315◦. 61 µ