Magnetic nulls and super-radial expansion in the solar corona

TL;DR

This paper provides the first observational evidence of coronal magnetic nulls in the solar corona using polarized light, introduces a new method to measure magnetic expansion, and highlights the diagnostic power of linear polarization for understanding solar wind and eruptions.

Contribution

It presents the first clear detection of coronal magnetic nulls in off-limb observations and introduces a novel technique to calculate expansion factors directly from polarization data.

Findings

Detection of magnetic nulls in the corona.

A new method to compute expansion factors without photospheric extrapolation.

Linear polarization as a diagnostic tool for coronal magnetic topology.

Abstract

Magnetic fields in the sun's outer atmosphere -- the corona -- control both solar-wind acceleration and the dynamics of solar eruptions. We present the first clear observational evidence of coronal magnetic nulls in off-limb linearly polarized observations of pseudostreamers, taken by the Coronal Multichannel Polarimeter (CoMP) telescope. These nulls represent regions where magnetic reconnection is likely to act as a catalyst for solar activity. CoMP linear-polarization observations also provide an independent, coronal proxy for magnetic expansion into the solar wind, a quantity often used to parameterize and predict the solar wind speed at Earth. We introduce a new method for explicitly calculating expansion factors from CoMP coronal linear-polarization observations, which does not require photospheric extrapolations. We conclude that linearly-polarized light is a powerful new diagnostic of critical coronal magnetic topologies and the expanding magnetic flux tubes that channel the solar wind.

Figures

a. Although pseudostreamers are three-dimensional magnetic structures, this two-dimensional idealization is a reasonable representation of the middle portion of a pseudostreamer that is oriented along an observer’s LOS. As discussed above and further demonstrated below, such an optimal orientation is in fact a necessary selection factor for obtaining clear linear polarization signatures. We will refer to magnetic nulls throughout this paper, however, we note that this effectively two-dimensional interpretation does not distinguish between a true, three-dimensional magnetic null point and topol

d) and the PFSS forward-modeled azimuth angle (Fig. 3d) suggests that the PFSS extrapolation underestimates the super-radiality of magnetic expansion to the north of the pseudostreamer. To further investigate this, we create maps of an expansion-factor proxy from both the CoMP and the synthetic PFSS data (Fig. 4). Using a predictor-corrector Euler algorithm and bilinear interpolation, we integrate polarization field lines following the azimuth down to the lower CoMP occulter. A linear-polarization expansion factor ($LPF$) is then determined from the local polarization field-line deviation (see

b), may be calculated from measurements at the solar photosphere along with, e.g., PFSS model evaluations at higher altitudes. This is then a standard input into empirical models that predict solar-wind speed at the Earth (e.g., Arge & Pizzo ([2000](#bib.bib3))).

Equations

$$f(\vec{r},\Phi )$$

$$f(\vec{r},\mathrm{d}\Phi )$$

$$f^{(3D)}(\vec{r}(s))$$

$$f^{(2.5D,sph.)}(\vec{r}(s))$$

$$f^{(2.5D,cyl.)}(\vec{r}(s))$$

Full text

Sun: magnetic fields — Sun: corona — Sun: solar wind

††journal: ApJL

1 Background

Measuring the strength and direction of the solar coronal magnetic field is a fundamental requirement for understanding solar activity and evolution, but remains a significant challenge. When simplifying assumptions are used such as the current-free or potential limit to the magnetic field in the corona, observations of the magnetic field at the solar surface may be used to approximate the coronal field. Such photospheric extrapolations – e.g., the potential-field-source-surface (PFSS) model – provide a good first-order characterization of the global coronal and interplanetary magnetic field (Newkirk & Altschuler, 1969; Schatten et al., 1969). However, this model ignores currents, removing the component of the magnetic energy that drives solar eruptions. Another possibility is to use the full vector of the magnetic field at the photospheric boundary and solve a non-linear force-free field (NLFFF) problem (under the assumption that magnetic forces dominate in the highly-conductive corona). However, issues such as a lack of magnetic dominance at the photosphere lead to inconsistencies and non-uniqueness of NLFFF models (De Rosa et al., 2009). The observed photospheric boundary on its own is therefore not sufficient to model the coronal magnetic field, motivating the explicit incorporation of coronal observations into magnetic models.

The Coronal Multichannel Polarimeter (CoMP) (Tomczyk et al., 2008) obtains daily observations of coronal emission-line polarization in the solar atmosphere, providing unique constraints on coronal magnetic models. Polarimetry has been used to determine magnetic fields at the solar surface since Hale’s observations in 1908. More recently, routine measurements of magnetism in its lower atmosphere, or chromosphere, have become possible (Ariste, 2002; Harvey, 2006). At higher coronal altitudes, emission-line measurements of circular polarization – sensitive to line-of-sight (LOS) oriented magnetic-field strength – have been rare (Harvey, 1969; Lin & Casini, 2000; Lin et al., 2004; Tomczyk et al., 2008). Obtaining routine circular-polarization measurements at the level of one Gauss in 15 minutes generally requires a larger-aperture telescope than the 20-cm CoMP (Tomczyk et al., 2016). On the other hand, linear polarization, which is approximately 100 times brighter than circular polarization, is measured by CoMP on a daily basis. As we now demonstrate, linear-polarization observations are ideally suited for probing certain topological and expansion properties of the coronal magnetic field.

2 Linear-polarization diagnostics: magnitude and Azimuth

The Fe XIII coronal emission line at $1074.7~{}nm$ is linearly polarized when unpolarized light emitted from lower layers of the solar atmosphere undergoes scattering in the corona. However, magnetic fields in the corona partially depolarize it, a process known as the Hanle effect (Sahal-Bréchot, 1977; Querfeld & Smartt, 1984; Arnaud & Newkirk, 1987). In the corona, at this wavelength, the depolarization of the scattered light preserves information about the direction of the magnetic field, but not its strength. In particular, if the local magnetic vector has an angle $\vartheta_{B}$ relative to the solar radius vector ($\hat{r}$) that is less than a critical “van Vleck” angle (van Vleck, 1925), its projection into the plane of sky (POS) lies parallel to the observed linear-polarization vector. However, if $\vartheta_{B}$ is greater than this critical angle, then the direction of the linear-polarization vector (known as the azimuth) is rotated $90^{\circ}$. The critical van-Vleck angle, $\vartheta_{B}=54.74^{\circ}$, arises because of an atomic-alignment dependence of the form $3\cos^{2}\vartheta_{B}=1$; at this critical angle the fraction of linearly-polarized light (linear-polarization magnitude) is zero.

CoMP observations and forward modeling have demonstrated that, in certain magnetic topologies, the loci of van-Vleck angle crossings can be picked out of linear-polarization magnitude observations as coherent and elongated dark features. The corona is optically thin in the Fe XIII line, and observations are intensity-weighted along the LOS, so multiple bright structures of varying magnetic orientations can in general contribute to the linear polarization observed by CoMP. Thus, the presence of distinct van-Vleck loci implies a degree of isolation and/or extension along the LOS of the magnetic structure. Under such conditions, these loci have proved useful for identifying coronal topologies. For example, the linear-polarization signatures expected from forward modeling cylindrical vs. spheromak magnetic flux ropes are distinct (Rachmeler et al., 2013). CoMP observations have provided one example of a possible spheromak (Dove et al., 2011); far more examples have been found of coronal prominence cavities (Gibson, 2015) that, for many sizes and shapes, possess “lagomorphic” (rabbit-headed) signatures in linear-polarization magnitude, consistent with a cylindrical magnetic flux-rope topology (Ba̧k-Stȩślicka et al., 2013, 2014).

3 Coronal magnetic nulls from linear-polarization magnitude

Coronal pseudostreamers and associated magnetic nulls are likely locations for a restructuring of the magnetic field to be triggered, potentially leading to coronal mass ejections (CMEs) (Antiochos et al., 1999; Török et al., 2011; Lynch & Edmondson, 2013). Pseudostreamers consist of multipolar closed field surrounded by unipolar open magnetic field (Hundhausen, 1972; Zhao & Webb, 2003; Wang et al., 2007). They are distinguished from so-called “helmet” (or bipolar) streamers, which are surrounded by open field of opposite polarities. Magnetic null points are expected in pseudostreamer topologies, as shown in Fig. 1a. Although pseudostreamers are three-dimensional magnetic structures, this two-dimensional idealization is a reasonable representation of the middle portion of a pseudostreamer that is oriented along an observer’s LOS. As discussed above and further demonstrated below, such an optimal orientation is in fact a necessary selection factor for obtaining clear linear polarization signatures. We will refer to magnetic nulls throughout this paper, however, we note that this effectively two-dimensional interpretation does not distinguish between a true, three-dimensional magnetic null point and topologically related structures such as a superposition of a line of nulls, a magnetic separator, or a hyperbolic flux tube (see, e.g., Titov et al. (2002); Aulanier et al. (2005); Parnell et al. (2010); Titov et al. (2011)).

Forward models of coronal pseudostreamer magnetic configurations have predicted a characteristic signature of three linear-polarization lobes (bright regions outlined by van-Vleck dark features; see Fig. 1b), which come to a point at the magnetic null (Rachmeler et al., 2014). In that work, CoMP data were presented which were consistent with the predicted pseudostreamer configuration, but which did not extend high enough to include the magnetic null. We have now observed a pseudostreamer with CoMP data that unequivocally matches the forward-modeled prediction of the pseudostreamer topology, and explicitly contains the signature of the converging lobes indicative of magnetic null(s) at their top (Fig. 1c).

Figs. 2 and 3 present another example of a pseudostreamer, and directly compare CoMP data to PFSS model predictions. Synthetic data are generated using the ‘pfss’ and ‘FORWARD’ (Gibson et al., 2016) packages of SolarSoft IDL (Freeland & Handy, 1998). Plasma weighting along the LOS is achieved by defining density radial profiles consistent with hydrostatic equilibrium: closed-field regions are assigned profiles that match coronal streamer observations (Gibson et al., 1999), while open-field regions are given profiles to match coronal hole observations (Guhathakurta et al., 1999).

Note that, although more asymmetric than the pseudostreamer of Fig. 1, the lobes of this pseudostreamer are apparent in linear polarization for both CoMP (Fig. 2c) and forward-modeled PFSS (Fig. 3c). PFSS magnetic field plots show this topology even more clearly (Figs. 3a-b). However, a visual inspection of the CoMP observations indicate a higher altitude null than the PFSS model possesses. The location of this CoMP null coincides with the base of the sharp linear feature extending above the pseudostreamer dome and up into the pseudostreamer stalk, as seen by the Solar Dynamics Observatory Atmospheric Imaging Assembly (AIA) images (Figs. 2a-b). We note that this feature is somewhat to the left of the center of the dome, but also that the CoMP linear polarization feature is similarly asymmetric.. The PFSS model also does not fully capture the bulging nature of the northward lobe seen in the CoMP images. The AIA 193 Å image (Fig. 2b) shows the presence of a prominence cavity within this lobe, so, it is possible the expanded lobe is a consequence of currents associated with a flux rope that are missed by the (current-free) PFSS model.

4 Coronal nonradial magnetic expansion from linear-polarization Azimuth

Another important coronal property is the height-dependent expansion of open magnetic field surrounding closed-field structures such as helmet streamers and pseudostreamers. The degree of super-radial expansion of open-field regions is of particular interest because it has been shown to be inversely correlated to solar-wind speed (Wang & Sheeley, 1990). A magnetic expansion factor, $f$ (see Appendix, Equation A2; also Fig. 3b), may be calculated from measurements at the solar photosphere along with, e.g., PFSS model evaluations at higher altitudes. This is then a standard input into empirical models that predict solar-wind speed at the Earth (e.g., Arge & Pizzo (2000)).

CoMP linear-polarization observations present us with a new resource for probing magnetic flux-tube expansion as a function of height in the low corona, independent of model extrapolations. Figs. 2d) and 3d) illustrate that van-Vleck crossings manifest not only as dark features in linear-polarization magnitude, but also as discontinuities in linear-polarization azimuth, so that diverging vs. converging magnetic fields have distinct appearances. Thus, azimuths crossing the van-Vleck angle – for example within the closed-field loops that make up the pseudostreamer– result in sharp jumps from red to blue. In contrast, the diverging magnetic fields to the north – associated with magnetic flux tubes open to the solar wind – have a gradual transition from blue to red and no van-Vleck crossing. However, the sharp jump at the van-Vleck angle may be blurred if, for example, the locus of these angles do not superpose along the LOS; this appears to be the case for the southern van-Vleck crossings for both the CoMP observations and the synthetic PFSS data.

A visual comparison of the CoMP-observed azimuth angle (Fig. 2d) and the PFSS forward-modeled azimuth angle (Fig. 3d) suggests that the PFSS extrapolation underestimates the super-radiality of magnetic expansion to the north of the pseudostreamer. To further investigate this, we create maps of an expansion-factor proxy from both the CoMP and the synthetic PFSS data (Fig. 4). Using a predictor-corrector Euler algorithm and bilinear interpolation, we integrate polarization field lines following the azimuth down to the lower CoMP occulter. A linear-polarization expansion factor ($LPF$) is then determined from the local polarization field-line deviation (see Appendix). We find that the CoMP data in this region indicates a more super-radial expansion than the PFSS model predicts. We note that this technique relies on an interpretation of the azimuth as an effectively two-dimensional POS field which only holds true in regions that are not heavily impacted by LOS effects, and that have magnetic field orientation tilted less than the van-Vleck angle from radial. The clear identification and localization of the van-Vleck inversion to the bottom-left (see also Fig. 3) gives us confidence that both conditions are met for the diverging, open magnetic field shown in this figure and used to calculate the $LPF$.

5 Discussion

We have shown how CoMP linear-polarization observations may be used to quantify both magnetic topology (height of magnetic nulls) and super-radial expansion ($LPF$ as a function of height and latitude). The detection of magnetic nulls provides a new way to identify critical points in the coronal field that facilitate solar eruptions. The linear-polarization expansion factor proxy, $LPF$, provides a new means for validating PFSS (and other) model predictions.

It has been pointed out that pseudostreamer-sourced solar wind may be inconsistent with the empirical relationship between PFSS-model-derived expansion factor and wind speed (Riley et al., 2011). In particular, the authors observed solar wind associated with a pseudostreamer to be slow, in contrast to the fast wind predicted by low PFSS-derived expansion factors. Wang et al. (2012) have also pointed out that the null-point height affects magnetic expansion, and that nonmonotonically expanding flux tubes may result in an overestimate of solar wind speed. The CoMP observations we have presented demonstrate how discrepancies between observations and the empirically-predicted solar wind speed could at least partially arise from an underestimate of magnetic expansion and null-point height by the PFSS model due to intrinsic limitations of photospheric extrapolations. It is also possible that reconnections at the open-closed magnetic boundary rather than flux-tube expansion is the dominant physical process controlling the solar wind speed associated with pseudostreamers, as has been argued in other models (e.g., Antiochos et al. (2011); Del Zanna et al. (2011)). These models too require accurate determination of coronal magnetic topology, in particular magnetic null point location. In conclusion, coronal linear-polarization measurements represent a valuable new type of data for validating and developing next-generation coronal and solar-wind acceleration models.

Acknowledgements. NCAR is supported by the National Science Foundation. SEG, KD, GdT and ST acknowledge support from the Air Force Office of Space Research, FA9550-15-1-0030. KD acknowledges funding from the Computational and Information Systems Laboratory, the High Altitude Observatory, and the Advanced Study Program. AIA data are courtesy of NASA/SDO and the AIA, EVE, and HMI science teams. We acknowledge helpful discussions with Anna Malanushenko, Ed Deluca, Antonia Savcheva and Scott McIntosh.

Appendix A Computing a Linear-Polarization Expansion Factor ($LPF$)

A.1 Background

A magnetic expansion factor measures the degree to which the cross sections of open magnetic flux tubes expand (or contract) in the solar corona as compared with their cross sections at the solar photosphere. In particular, magnetic expansion factor can be defined as:

[TABLE]

where $\mathcal{S}(\vec{r},\Phi)$ is the cross-sectional area of a tube carrying magnetic flux, $\Phi$, at position, $\vec{r}$, along its axis, $R_{\odot}$ is the solar radius, and $R_{\odot}^{2}/r^{2}$ is a correction factor accounting for the natural expansion of surfaces as $r$ increases in spherical geometry (referred to as “radial expansion” in the text).

From Eq. A1 it is possible to define a completely local expansion factor which does not depend on the surface cross section, $\mathcal{S}$. For this we employ “elemental” flux tubes, which possess an infinitesimal cross section surrounding a magnetic field line. Let $\mathrm{d}\Phi=B\mathrm{d}\mathcal{S}$ be an elemental magnetic flux tube, where $B$ is the magnetic field strength along the flux tube axis and $\mathrm{d}\mathcal{S}$ is the elemental flux tube cross section. Then Eq. A1 applied to elemental magnetic flux tubes leads to

[TABLE]

where $f(\vec{r})$ is the local expression of the expansion factor.

When the magnetic field is known in the entire volume, one can compute the local expansion factor from Eq. A2 by simply integrating field lines and evaluating the magnetic field strength at each step along the field lines (Wang & Sheeley, 1990). When the magnetic field direction is known, but not the magnetic field strength, we must return to Eq. A1, which we will use to define and justify a linear-polarization expansion factor proxy, $LPF$, determined, e.g., from CoMP linear polarization measurements.

A.2 $LPF$ proxies

To evaluate the $LPF$, we build upon methodology developed to compute the squashing factor of magnetic flux tubes (Pariat & Démoulin, 2012; Tassev & Savcheva, 2017). The squashing factor is a measure of the distortion of magnetic flux tubes that provides a means of identifying quasi-separatrix layers (QSLs) in a 3D magnetic field (Demoulin et al., 1996; Titov et al., 2002), and has been extensively used in solar flare analysis to relate photospheric flare ribbons to the topology of the coronal magnetic field (Schmieder et al., 1997; Masson et al., 2009; Savcheva et al., 2012a; Dalmasse et al., 2015). QSLs are regions of strong gradients in connectivity of magnetic field lines, and are likely regions for current sheet/layers formation and magnetic reconnection (Aulanier et al., 2005; Janvier et al., 2013).

The method is based on integrating a field line and calculating the field line deviation along it, $\delta\vec{r}$, from the local variations of the magnetic field vector as obtained from its local interpolation (Longcope & Strauss, 1994; Titov et al., 2003; Tassev & Savcheva, 2017). Once $\delta\vec{r}$ has been integrated along a field line, we can compute $\delta\vec{r}_{\perp}$, the component of $\delta\vec{r}$ perpendicular to the field line. For a 3D flux tube of circular cross section, $\delta\vec{r}_{\perp}(s)$ effectively represents the diameter of the flux tube at position $\vec{r}_{s}$ along the flux tube axis. $\delta\vec{r}_{\perp}(s)$ can thus be used to compute an estimate of the local expansion factor using Eq. A1.

Because of the LOS integration effect, the azimuth computed from CoMP only provides an effectively two-dimensional magnetic field direction. Thus, the expansion factor computed from CoMP linear-polarization measurements can only be a proxy of the true, local expansion factor of the coronal magnetic field. Using Eq. A1, we define three, slightly different, proxies of expansion factor

[TABLE]

Eq. A3 assumes that the CoMP data is associated with 3D magnetic flux tubes of circular cross section, $\mathcal{S}(\vec{r}(s))=\pi(\|\delta\vec{r}_{\perp}(s)\|/2)^{2}$, and that axes are in the POS. Eq. A4 assumes axisymmetry around the North-South axis of the spherical Sun, in which case the cross section of a flux tube is $\mathcal{S}(\vec{r}(s))=\|\delta\vec{r}_{\perp}(s)\|2\pi r(s)\sin\left(\theta(s)\right)$. Finally, Eq. A5 assumes invariance in the direction of the line of sight, meaning that the Sun appears as a cylinder of infinite length, $L$, in the LOS direction; magnetic flux tubes thus have a cross section, $\mathcal{S}(\vec{r}(s))=\|\delta\vec{r}_{\perp}(s)\|L$, and the correction factor that must be applied to remove the natural expansion of surfaces as $r$ increases is now $R_{\odot}/r$. For the case presented in this paper, we have compared all three $LPF$ proxies and found that although they vary somewhat in magnitude, they are very similar in structure. The conclusions we make comparing the CoMP and PFSS $LPF$ measures are independent of which measure is used; in Figure 4 we show the 3D proxy defined by Eq. A3.