A Method for a Pseudo-Local Measurement of the Galactic Magnetic Field
aa r X i v : . [ a s t r o - ph . GA ] J a n Draft version January 5, 2021
Typeset using L A TEX
RNAAS style in AASTeX62
A Method for a Pseudo-Local Measurement of the Galactic Magnetic Field
Steven R. Spangler Department of Physics and Astronomy, University of Iowa
ABSTRACTMuch of the information about the magnetic field in the Milky Way and other galaxies comes frommeasurements which are path integrals, such as Faraday rotation and the polarization of synchrotronradiation of cosmic ray electrons. The measurement made at the radio telescope results from thecontributions of volume elements along a long line of sight. The inferred magnetic field is thereforesome sort of average over a long line segment. A magnetic field measurement at a given spatial locationis of much more physical significance. In this paper, we point out that HII regions fortuitously offersuch a “point” measurement, albeit of one component of the magnetic field, and averaged over thesightline through the HII region. However, the line of sight (LOS) through an HII region is muchsmaller (e.g. 30 - 50 pc) than one through the entire Galactic disk, and thus constitutes a “pseudo-local” measurement. We use published HII region Faraday rotation measurements to provide a newconstraint on the magnitude of magnetohydrodynamic (MHD) turbulence in the Galaxy, as well asto raise intriguing speculations about the modification of the Galactic field during the star formationprocess.
Keywords: interstellar medium — interstellar magnetic fields — interstellar plasma1. INTRODUCTIONIt has been known for decades that the interstellar medium is threaded by a magnetic field. This is not surprising,since almost all phases of the interstellar medium (ISM) have sufficient level of ionization to satisfy the conditionsfor a plasma parameter Λ e ≫ δbB where δb is the rms value ofthe turbulent component and B is the magnitude of the large scale component, range from a few tens of percent toa number in excess of unity. For example, Ferri`ere (2011) and Ferri`ere (2020) review our knowledge of the Galactic Corresponding author: Steven [email protected] field, and pay attention to the relative magnitudes of the regular and turbulent fields, citing results of δbB ≥ µG , that of anisotropic fluctuations 2 µG , and that of the regular component (ordered over Galactic scales) is 2 µG . Itshould be emphasized that these estimates are dependent on assumptions on the outer scale of the turbulence, i.e. thescale which distinguishes a Galactic-scale “deterministic” component from a turbulent component. Nonetheless, thesereviews make clear that there is substantial observational evidence for a turbulent component to the Galactic magneticfield that is comparable in magnitude, or even larger than an ordered, large scale field. In this paper, we will considerthe dimensionless turbulent amplitude δbB to be a free parameter, subject to the observational constraints cited above.Knowledge of the magnitude and spectrum of the Galactic MHD turbulence is important to Galactic astrophysics, andastrophysics in general, since this turbulence contributes to processes such as heating of the ISM through dissipation ofturbulence, confinement and propagation of the cosmic rays, turbulent transport of heavy elements, and modificationor regulation of the star formation process.One of the difficulties in determining the relative roles of a large scale and a turbulent component of the Galactic fieldis the fact that some of the main diagnostic techniques such as Faraday rotation, polarization of Galactic synchrotronradiation, and dust polarization consist of path integrals over long lines of sight through the Galaxy. What this meansis that we get some sort of average of the Galactic field over a long line of sight (LOS) through the Galaxy, rather thanthe value at a set of points along that line of sight. Zeeman splitting of OH and methanol maser transitions are notsubject to this limitation, since the region of emission is very localized (Crutcher and Kemball 2019). However, asdiscussed in Crutcher and Kemball (2019) regions of maser emission have very high densities compared to the generalISM, and are products of the star formation process. We will return to this point in Section 4.In this paper, we discuss a type of observation which utilizes the Faraday rotation technique to obtain an estimate ofthe Galactic field (strictly speaking one component of the field) in a localized region, rather than an entire line of sightthrough the Galaxy. The localized region is comprised of the volume of an HII region. In this sense, the techniqueprovides a “local” measurement. In what follows, we first summarize the basis of the technique of Faraday rotation,then discuss recent observational results which allow it to be used in a “local” measurement of the Galactic field.Faraday rotation is a fundamental plasma process which consists of rotation in the plane of polarization of radiowaves due to the presence of a magnetized plasma (Nicholson 1983). It is used as a diagnostic technique in laboratoryplasmas as well as in radio astronomy. In the present context, it is applied to radio waves emitted by extragalacticradio sources and which propagate through the plasma of the interstellar medium (ISM). The change in the positionangle relative to what would be measured in the absence of the plasma, ∆ χ , is given by∆ χ = (cid:20)(cid:18) e πm e c (cid:19) Z L n ~B · ~ds (cid:21) λ (1)In Equation (1) e, m e , and c are the customary fundamental physical constants, n is the electron density, ~B is thevector magnetic field, both functions of position along the line of sight, and ~ds is an increment in the line of sight fromthe source to the radio telescope. Faraday rotation is a path integral measurement rather than a point measurement,so there are contributions to plasma throughout the ISM along the line of sight. The quantity in square brackets inEquation (1) is termed the rotation measure (RM), and is usually quoted in SI units of radians/m . For most linesof sight through the ISM, the magnetic field measured is not only an average over a long line through the ISM, it isalso an average weighted by the plasma density. Information on many lines of sight around the sky is used to obtainestimates of the large scale structure of the Galactic field, as well as the contribution of the turbulent component.The observational impetus for this paper is a set of studies of the HII regions Rosette Nebula and W4 by A. Costaand coworkers (Savage et al. 2013; Costa et al. 2016; Costa and Spangler 2018). Costa and co-workers carried outFaraday rotation measurements with the Very Large Array (VLA) of background radio sources viewed through thesetwo HII regions. The main results of these investigations were: (1) rotation measures (RM) for sources viewed throughthe HII regions were typically much larger than for nearby sources whose lines of sight probed the Galaxy, but not theHII regions, (2) although there were large variations in the RM from source-to-source when the LOS passed throughthe HII region, there was a well-defined average RM (magnitude and sign) for the HII region, particularly the Rosette,(3) the sign of the HII region-induced RM (which dominated the total RM) was consistent with the sign of the large Table 1.
HII Region Average Magnetic Field ComponentsHII region l b < B k > µG Ref.Rosette 206.5 -2.1 +2.73 Costa et al. (2016)W4 135 +1.0 -2.29 Costa and Spangler (2018)Sh2-264 195 -12 +2.2 Harvey-Smith et al. (2011)Sivan 3 144.5 +14 -2.5 Harvey-Smith et al. (2011)Sh2-171 118.0 +6 -2.3 Harvey-Smith et al. (2011)Sh2-220 161.3 -12.5 -6.3 Harvey-Smith et al. (2011) scale Galactic field in that part of the sky. The RM due to the Rosette Nebula (Galactic longitude = 206.5 degrees)was dominantly positive, while that for W4 (Galactic longitude = 135 degrees) was dominantly negative.The important point is that the magnetic field in the “RM anomaly” associated with the HII region is confined tothe volume of the HII region rather than being an average over a line of sight of many kiloparsecs in the Galacticplane. The radii of both HII regions are approximately 20-25 parsecs, so the magnetic field is “local” in this sense.As discussed in the text earlier in this section, in view of the apparent large amplitude, turbulent component of theGalactic field, a reasonable a-priori expectation would be that the magnetic field at a point in the Galactic plane, andaveraged over a volume of radius 20 - 25 parsecs, would have poor correlation with the ordered, large scale Galacticfield. The observations of the Rosette Nebula and W4 do not bear out this expectation. In this paper, we explore theconsequences of this result for our inference about the Galactic magnetic field.2. ADDITIONAL OBSERVATIONAL RESULTS ON FARADAY ROTATION THROUGH HII REGIONSTwo measurements is, of course, very poor statistics, even by the base standards of astronomy. There is a need foradditional measurements of the sort presented by Costa and collaborators. Such measurements were provided by theobservations of Harvey-Smith et al. (2011), who carried out similar observations, for similar goals, of somewhat largerHII regions like the λ Orionis HII region. In collecting HII regions from the sample of Harvey-Smith et al. (2011),we restricted attention to objects in the outer two quadrants of the Galactic plane, i.e. those with 90 ◦ ≤ l ≤ ◦ .Harvey-Smith et al. (2011) reported the the average value of the line-of-sight component of the magnetic field < B k > for each HII region, and this quantity was also available from the analyses of Savage et al. (2013); Costa et al. (2016);Costa and Spangler (2018).In the studies of Costa and coworkers, the mean LOS component of the field, < B k > , was not explicitly measured;analyses of the data concentrated on fitting plasma shell models to the sets of RM values for each HII region. Inaddition to other parameters, these models were dependent on the magnitude and orientation of the magnetic fieldoutside the HII region (see Section 5.1 of Costa et al. (2016)). A limiting type of model was one in which there was noamplification of the external magnetic field in the HII region (referred to by Costa as the “Harvey-Smith model”). Inthis case, the model fits give exactly the same quantity as the parameter < B k > in Harvey-Smith et al. (2011). Wetook values from Table 4 of Costa et al. (2016) and Table 6 of Costa and Spangler (2018) to obtain values of < B k > ,which should be directly comparable to this quantity in Harvey-Smith et al. (2011). Once again, < B k > has a signand a magnitude. The values of < B k > from Harvey-Smith et al. (2011), Costa et al. (2016), and Costa and Spangler(2018) are plotted in Figure 1, and given in Table 1.The clear result from Figure 1 is that for all 6 HII regions, the “local” magnetic field defined by the average over thevolume of the HII region, has the same polarity as the Galactic-scale, mean magnetic field. In none of these 6 casesare fluctuations on scales of a few tens of parsecs capable of making the “local” polarity opposite to the average overseveral kiloparsecs.3. STATISTICS OF A LOCAL, LINE OF SIGHT MEASUREMENT OF ~B IN A TURBULENT MEDIUMIn this section, we discuss the statistics of the LOS component of the magnetic field at a given point, when themagnetic field consists of the vector sum of a large scale, ordered magnetic field ~B , and a superposed, spatiallyrandom turbulent field ~δb . For ease of following the discussion, the Appendix gives a glossary of mathematical symbolsand variables used in the following. ●● ●● ○○ ○○ ○○ ○○
120 140 160 180 200 220 l - - - b HII Region Magnetic Fields
Figure 1.
Measurements of < B k > from Faraday rotation measurements through 6 HII regions in the 2nd and 3rd Galacticquadrants. Solid (blue) symbols indicate < B k >> < B k >< p | < B k > | . The locationsof the 6 sources used are plotted in Galactic coordinates. In the simplest case of a purely azimuthal field, the projection of thelarge scale Galactic field on the LOS would change sign at l = 180 ◦ . In all cases, the polarity of < B k > is the same as that forthe large scale Galactic field in that quadrant. Geometry of the Line of Sight and the Galactic Magnetic Field
The geometry of the situation is shown in Figure 2. The z axis is defined by the path from the radio source throughthe HII region and to the observer, with the positive direction towards the observer. The x axis is in the planeof the sky, and points toward decreasing galactic longitude. We also simplify the situation by adopting the goodapproximation that the x-z plane coincides with the Galactic plane. Finally, the y axis is perpendicular to the Galacticplane and completes a right-handed coordinate system.The observational results of Harvey-Smith et al. (2011); Savage et al. (2013); Costa et al. (2016); Costa and Spangler(2018) show that there is a well-defined mean field throughout the HII regions, with a magnitude equal to or greaterthan the general Galactic field. Given the known presence of plasma turbulence in the interstellar medium (ISM), weinterpret this average HII region field to be the vector sum of the large scale galactic field ~B and a turbulent field ~δb .As indicated in Figure 2, the angle between the line-of-sight (LOS) and the large scale field is Θ. We also assume, forthis analysis, that the large-scale field is in the Galactic plane.Given the configuration illustrated in Figure 2, we define a second coordinate system which is convenient for describ-ing the turbulence. The k direction is in the direction of the mean field, the p direction is perpendicular to the meanfield, and in the Galactic plane, and the ⊥ direction completes the right-hand coordinate system and is perpendicularto the Galactic plane. These three directions are defined by the unit vectors ˆ e k , ˆ e p , and ˆ e ⊥ . Model for the Magnetic Field in the Vicinity of an HII Region
Given this latter coordinate system, we express the magnetic field (mean field plus turbulence) by ~B = ( B + δb k )ˆ e k + δb p ˆ e p + δb ⊥ ˆ e ⊥ (2)In analogy with the Alfv`enic turbulence in the solar wind, we expect < ( δb p ) > = < ( δb ⊥ ) > and < ( δb p ) > ≫ < ( δb k ) > .The model we adopt is that δb k , δb p , and δb ⊥ are random, 0-mean quantities with an amplitude dominated byfluctuations on the outer scale l, l ≥ R , where R is the radius of the HII region. In the cases of the Rosette Nebula(Costa et al. 2016) and W4 (Costa and Spangler 2018), the values of 2 R are 38 and 50 parsecs, respectively. Thesevalues are comparable to, or slightly less than the outer scale of the 2D component of interstellar turbulence claimedby Minter and Spangler (1996). Figure 2.
Cartoon illustrating geometry of observer, polarized extragalactic radio source, and HII region. A Galactic magneticfield is illustrated, consisting of a large scale, ordered component ~B inclined at an angle Θ with respect to the line of sight,and a superposed turbulent component ~δb . The plane of the paper is the Galactic plane, the z axis is defined by a ray fromthe radio source to the observer. The coordinates x and y complete a right-handed coordinate system. In this illustration, onlyone component of the turbulent field is shown, δb p , which is the component perpendicular to the large scale field and in theGalactic plane. To lend realism, the trace shown represents 6 hours of turbulent solar wind magnetic fluctuations measured bythe WIND spacecraft on September 25, 2020. We adopt a very simple model for the average observed RM to a set of extragalactic radio sources viewed throughan HII region, RM obs = RM bck + RM neb (3)where RM bck is the background RM due to the plasma in the Galaxy which would be observed in the absence ofthe nebula, and RM neb is the average RM for lines of sight through the HII region. RM neb is given by the simpleexpression RM neb = 2 CnRAB z (4)where C is the set of fundamental constants in curved brackets in Equation (1), n is the mean electron density inthe ionized portion of the HII region (determined by independent measurements such as intensity of thermal radioemission in Savage et al. (2013); Costa et al. (2016); Costa and Spangler (2018)), R again is the radius of the HIIregion, and B z is the z component of the interstellar field (mean field plus turbulence) at the location of the HIIregion. The variable A ≥ A = 1), while Costa et al. (2016);Costa and Spangler (2018) considered physically-based models in which A ≤
4. The case A = 1 was considered aswell. B z is obtained from Equation (2) as B z = ~B · ˆ e z (5)Substituting Equation (2) into this expression and evaluating the dot products of unit vectors gives B z = ( B + δb k ) cos Θ − δb p sin Θ (6)The δb ⊥ component of the turbulent field makes no contribution because, for lines of sight in or close to the Galacticplane, it is perpendicular to the LOS. The negative sign in the second term arises because positive δb p is pointing awayfrom the observer (Figure 2).Equation (6) shows the important, though obvious fact that the polarity of B z (i.e. the same as ~B or opposite)depends on both the amplitude of the fluctuations relative to the mean field as well as the angle Θ between the largescale field and the LOS.This expression can be rewritten in terms of dimensionless turbulence amplitudes (turbulent “modulation indices”)as B z = B (cos Θ + x cos Θ − y sin Θ) (7) x ≡ δb k B (8) y ≡ δb p B (9)where x, y are dimensionless, 0 mean random variables with assumed properties √ < x >, p < y > ≤
1. With allof this, the observed RM is then RM obs = RM bck + 2 CnRAB (cos Θ + x cos Θ − y sin Θ) (10) RM obs = RM bck + RM n (cos Θ + x cos Θ − y sin Θ) (11)with an obvious definition of the net nebular RM n . It is worth emphasizing that RM n is a property of the nebulaitself, independent of the ISM or the geometry of the line of sight.The question considered in this paper then is whether the second term has the same sign as the first, and if not, if itis large enough to cause the observed RM to have the opposite sign as RM bck . Equation (11) shows that this dependson two factors.1. The first is the relative magnitudes of RM bck and RM n . Savage et al. (2013); Costa et al. (2016); Costa and Spangler(2018) found that for both the Rosette Nebula and W4, the largest absolute magnitudes of the RMs throughthe nebula were of the order of 1000-1500 rad/m , whereas the background Galactic RMs in those parts of thesky had absolute magnitudes of the order of 100 - 150 rad/m . This indicates that the magnitude of RM n couldbe as high as 10 times that of RM bck .2. The statistics of the quantity ξ = cos Θ + x cos Θ − y sin Θ (12)are obviously crucial. For the observed RM to be of the opposite polarity to RM bck , it is necessary that ξ < ξ < − ∆ where ∆ is defined as the ratio of background to nebular RM.For the remainder of this section, we turn to a discussion of the statistics of ξ . Statistics of the Random Variable ξ We begin by making a further simplification and ignoring the x term in ξ . The justification for this is the assumptionstated above that for Alfv`enic turbulence, the fluctuations in the parallel component (to the mean field) are smallcompared to the transverse components. In addition, both the mean field and the parallel component contributionsto the stochastic variable ξ are proportional to cos Θ. We therefore proceed to consider the statistics of the variable ξ = cos Θ − y sin Θ (13)Only y is a random variable, so we model it as a Gaussian-distributed, random variable with zero mean, with aprobability distribution function p ( y ) = 1 √ πσ exp( − y / σ ) (14)The quantity σ is the dimensionless amplitude of the fluctuations of the p component of the turbulence, and is a crucialparameter in our analysis.Given this expression we can generate the probability distribution function (pdf) of the stochastic variable ξp ( ξ ) = 1 √ πσ (cid:18) − (cos Θ − ξ ) σ sin Θ (cid:19) (15)For economy of notation, let a ≡ cos Θ, b ≡ sin Θ, with the obvious identity b = √ − a , so that p ( ξ ) = 1 √ πσb exp (cid:18) − ( ξ − a ) σ b (cid:19) (16) p ( ξ ) = 1 √ πd exp (cid:18) − ( ξ − a ) d (cid:19) (17)with d ≡ σb . The distribution function of ξ is also Gaussian, with a nonzero mean value < ξ > = a = cos Θ.A sample plot of p ( ξ ) is shown in Figure 2. This figure shows that for sufficiently large-amplitude turbulence( σ = 0 .
70 in this case), and sufficiently large angles of the mean field with respect to the line of sight (Θ = 70 ◦ ) theprobability of the local field being reversed with respect to the mean field is substantial. - - ξ p ( ξ ) Figure 3.
The probability distribution function of the variable ξ for the case in which the dimensionless amplitude σ of themagnetic fluctuations in the p component of the magnetic field is σ = 0 .
70 and the angle Θ between the line of sight and themean magnetic field is 70 degrees. In this case, the probability of ξ being negative is substantial. Expression for P − , the Probability of a Polarity Reversal Given the expression in Equation (16) or (17) for the probability distribution function of the parameter ξ , we havethe quantity of interest, the total probability that the line-of-sight component of the magnetic field at a given locationin the ISM will be reversed from that of the mean field, We note this quantity by P − , and the obvious expression forit is P − = Z − ∆ −∞ dξp ( ξ ) (18)By inspection, it can be seen that this integral is the same as P − = Z − (∆+ a ) −∞ dyp G ( y, d ) = Z ∞ ∆+ a dyp G ( y, d ) (19)where p G ( y, d ) is a zero mean, normalized Gaussian pdf with standard deviation d .If we make a final change of variables, y → t ≡ y √ d , then Equation (19) becomes P − = 1 √ π Z ∞ X dte − t (20) X ≡ ∆ + a √ d = (∆ + cos Θ) √ σ sin Θ (21)Equation (20) is very close to the standard form of the Error Function (Beckmann 1967) erf ( X ) ≡ √ π Z X dte − t (22)Comparison of Equations (20) and (22) allows the final, compact formula for the total probability of polarity reversal, P − = 12 [1 − erf ( X )] = 12 erf c ( X ) (23)where erf c ( X ) is the complement of the error function.Equation (23) is plotted in Figure 4. As the argument X →
0, the probability of reversal P − → .
5. This makes X P ( X ) Reversal Probability Function
Figure 4.
The probability that the line of sight component of the magnetic field at the location of an HII region is opposite tothat of the global, Galactic field. The probability P − is a function of the parameter X = (∆+cos Θ) √ σ sin Θ . complete sense, since small argument X corresponds to strong turbulence and/or a mean field nearly perpendicularto the line of sight. In this case, the local magnetic field is dominated by the turbulent component δb p , which has a50 percent change of being positive and and a 50 percent chance of being negative.The case of X → ∞ corresponds to weak turbulence, and/or a mean field increasingly directed towards the observer.In this case, the probability of a reversal in the line-of-sight component of the magnetic field approaches zero.The probability of a reversal in the line of sight component is determined by the parameter X , and it is of interestto calculate this quantity for some observationally-relevant cases. Calculations Relevant to Observations of HII Regions
In this section, we briefly apply these formulas to the observations presented in Section 2. A total of 6 HII regionsdistributed over a range of Galactic longitudes precludes a convincing and sophisticated statistical analysis. However,the results presented in Figure 1 suggest that theoretical values of P − ≥ .
20 would begin to have an uncomfortableconfrontation with the observations. On the other hand, if characteristics of the HII regions, amplitude of ISMturbulence and direction of the line of sight indicate a value of P − ∼ .
01, our calculations are completely consistentwith the observations discussed in Section 2. The formula given in Eq. (21) shows that the parameter X depends onΘ, σ , and ∆. In what follows, we adopt a value of ∆ = 0 .
1, which is in accord with the observations of the RosetteNebula and W4 discussed above.Savage et al. (2013); Costa et al. (2016) in their study of the Rosette nebula, found that reasonable fits of simplemagnetized shell models to the set of RM measurements gave a value of Θ ≃ ◦ if the interstellar magnetic field isamplified in the shell by shock compression, and Θ ≃ ◦ if there were no change in the magnitude of the magneticfield in the shell, with respect to the general ISM (see Table 5 of Costa et al. 2016). Savage et al. (2013) noted thatan independent estimate, based on models of the Galactic magnetic field, would give Θ ≃ ◦ at the location of theRosette. Costa and Spangler (2018) used a value of Θ = 55 ◦ for W4, if the Galactic magnetic field is approximatedas azimuthal. Thus values of 50 ◦ ≤ Θ ≤ ◦ are relevant to the results shown in Figure 1.The amplitude of the turbulent magnetic field on spatial scales of tens of parsecs is not well known, but results fromthe literature discussed in Section 1 indicate that it is of order the mean, large scale field. Thus our parameter σ ifprobably of the order of unity, or a large fraction thereof. A value of σ = 0 . ≃ √ would have the magnitude offluctuating field (quadratic sum of δb p and δb ⊥ components) equal to the magnitude of the large scale field.A value of σ = 0 .
70 and Θ = 70 ◦ gives X = 0 .
48; reference to Figure 4 shows that the probability of reversal P − is about 24 %. In such a case it is mildly surprising that all 6 HII regions discussed in Section 2 produce “Faradayrotation anomalies” with the same polarity as the large scale field. On the other hand, a slightly smaller value of σ = 0 .
50 and Θ = 50 ◦ , also consistent with the Costa et al. (2016) results for the Rosette Nebula and W4, would have X = 1 .
37. In this case, Figure 4 shows that the probability of a polarity reversal for a single source is less than 5 %,and the results of Figure 1 are unremarkable.This discussion shows that future investigations on more HII regions could provide more stringent limits on thelarge scale turbulent fluctuations in the Galactic magnetic field. The formulas presented in this section provide amathematical vocabulary for discussing current as well as future observations. For example, it is clear that theprobability of a reversal becomes much greater for Θ ∼ π . HII regions in the close vicinity of the Galactic anticenterwould provide particularly sensitive diagnostics for the amplitude of large scale turbulence.4. CONCLUSIONS FROM ZEEMAN EFFECT MEASUREMENTS OF STAR FORMATION REGION MASERSThe analysis in Section 3 assumes that HII regions occur in the presence of an unmodified Galactic magnetic field. Inthe simplified stellar bubble models presented in Savage et al. (2013); Costa et al. (2016); Costa and Spangler (2018)the expanding bubble modifies the Galactic field, but the RM enhancement due to the bubble is proportional to theGalactic field (mean field plus turbulent contribution at the location of the HII region) in the absence of the HII region.In the model considered by Harvey-Smith et al. (2011) the HII region consists of an ionization of the interstellar gaswith no modification of the Galactic field. Harvey-Smith et al. (2011) describe this as the HII region “lighting up” theGalactic field that exists at the location of the HII region.These viewpoints may be physically naive. HII regions are produced by hot, massive stars that are products ofthe star formation process. Star formation regions, in turn, result from contraction of gas from the general ISM, itsconversion to molecular form, and an associated increase in gas density. It therefore seems likely that the process offorming dense clumps in molecular clouds from which stars form results in major modification of the Galactic field,perhaps to the extent of totally randomizing the vector field from its value in the absence of the star formation region.In this light, the observational results presented in Section 2 are even more remarkable; the preservation of the polarityof the large scale field occurs in spite of both the presence of natural turbulence in the Galactic field, as well as thepossibly major variations to the field induced by the star formation process.The HII region results presented in Section 2 are not the only “local” measurements of the magnetic field in the ISM.Hydroxyl, methanol, and water masers in star formation regions display the Zeeman effect, permitting measurementsof the magnitude and sign of the LOS component of the magnetic field in the maser region (Crutcher and Kemball2019). These masers occur in regions of very high density compared to the ISM far from star formation regions, andthus probe gas that has been substantially modified in density, and presumably magnetic field, by the star formationprocess (Crutcher and Kemball 2019).In spite of this, observations of the polarity of the magnetic field in maser regions show results similar to thosepresented in Section 2. Fish et al. (2003) examined Zeeman effect measurements for more than 50 massive starformation regions. Although they did not find a general Galactic-wide magnetic field direction, they did find that “. . . in the solar neighborhood the magnetic field outside the solar circle is oriented clockwise as viewed from the northGalactic pole . . . ”. This is the same sense as shown in Figure 1 on the basis of the HII region results. A subsequentstudy by Green et al. (2012) yielded similar results. Green et al. (2012) reported measurements of the magnetic field(sign and magnitude of the line-of-sight component of the magnetic field) for 14 Zeeman pairs in 6 high mass starformation regions in the Carina-Sagittarius arm (4th Galactic quadrant). Their results show a dominant magneticpolarity in this spiral arm (See Figure 7 of Green et al. (2012) . The significance of these observational results was0summarized in Han (2017), who noted that they implied “turbulence in molecular clouds and violent star-formingregions seem not to alter the mean magnetic field direction, although the field strength is much enhanced from a few µ G to a few mG ”.5. SPECULATIONS FROM THE STUDY OF HYDRODYNAMICAL TURBULENCEThe results for the HII regions presented above, and illustrated in Figure 1 indicate, at the very least, potentiallysignificant limits on the magnitude of turbulent magnetic field fluctuations in the Galactic field. The results on starformation region OH masers from Fish et al. (2003) and Green et al. (2012) indicate, and the HII regions resultsmay indicate, a far stronger and mysterious physical process: a memory of the large scale polarity of the Galactic fieldas it is modified during the contraction and increase in density associated with star formation.The obvious question which then arises is the meaning of “the persistence of memory” in the large scale, Galacticmagnetic field. The problem arises because of the assumed isotropic nature of turbulent fluctuations on scales muchsmaller than the outer scale of the turbulence. Such fluctuations, deep in the inertial subrange or close to the dissipationrange, should have no tendency to respect large scale asymmetries in the Galactic plasma. It is worthwhile to look forinsight from the study of hydrodynamical turbulence. There are obvious physical differences between the turbulencein the diffuse plasmas of the interstellar medium and the turbulence in gases and liquids in laboratory and industrialsettings. However, there are important similarities, too, and both the quality of the diagnostics and the shear amountof effort expended in the study of hydrodynamic turbulence exceeds that which is possible for astrophysical turbulence.A review of our understanding of “passive scalars” in hydrodynamic turbulence is presented in Warhaft (2000).Passive scalars are quantities which are convected by fluid turbulence, without affecting the flow field. Warhaft (2000)considered as a prime example temperature fluctuations in hydrodynamic turbulence. An astrophysical example of apassive scalar could be the plasma density fluctuations responsible for interstellar radio scintillations. Warhaft (2000)presents an illuminating discussion of the differences between the engineering and physics views of turbulence, andoffers an admonition that is certainly relevant for Galactic astronomers as well: “ The physicists, with their quest foruniversality, have concentrated on inertial and dissipation scales, the hope being, if the Reynolds (Re) and Peclet (Pe)numbers are high enough, that the small-scale behavior will be independent of the large scales. We address both thelarge and small scales in this review and argue that this division is artificial and dangerous”. Warhaft (2000) thengoes to present laboratory and theoretical evidence that small scale passive scalar fluctuations (primarily temperaturefluctuations) are anisotropic with an anisotropy that reflects large scale asymmetry. In the cases discussed by Warhaft(2000) this large scale asymmetry was most typically a transverse temperature gradient imposed on a turbulent shearflow. Warhaft (2000) also shows that the small scale passive scalar fluctuations possess a skewness and kurtosis thatexceeds that present in the fundamental turbulent velocity field.It is not clear to me if the clear and insightful presentation of Warhaft (2000) holds a key to understanding themysterious observation described in this paper, or if any such comparison would be a hallucinogenic extrapolation oflaboratory results far outside their range of applicability. While Warhaft’s review does describe a documented persis-tence of large scale anisotropy on small turbulent scales, as appears to be indicated in our Galactic ISM observations,it must be emphasized that the discussion of Warhaft (2000) is concerned with passive scalars in the turbulence.In the Galactic ISM context, that would be plasma density, temperature, ionization fraction, and metallicity. Themagnetic field, on the other hand, should share status with fluid velocity as a “primary excitation” of the turbulence,if the ISM turbulence is Alfv´enic. Nonetheless, it would be worthwhile to consider further if the extensive, laboratorybased studies of turbulence from an engineering perspective can illuminate the nature of plasma turbulence in theinterstellar medium. 6. SUMMARY AND CONCLUSIONS1. In the case of the Galactic-anticenter-direction HII regions Rosette Nebula and W4, the mean magnetic field inthe HII region has the same polarity as the general Galactic field in that part of space, despite the presence ofstrong MHD turbulence which would be expected to randomize the polarity. An additional 4 HII regions studiedby Harvey-Smith et al (2011) share this property.2. We derive formulas for the probability of the “point” magnetic field having the same polarity as the general field,which depends on properties of the turbulent magnetic fluctuations as well as the details of the line of sight. Forreasonable estimates, this probability could be non-negligible, being as high as 25 % for plausible characteristics.1As a result, the observation of 6 anticenter-direction HII regions, all having the same magnetic polarity of thegeneral field, becomes mildly curious.3. The results on HII regions support similar results on the polarity of star formation OH masers presented byFish et al. (2003) and Green et al. (2012). In those studies, the polarity of the magnetic field in maser regionsis correlated among star formation regions that are spatially separated, and the polarity is apparently determinedby that of the large scale Galactic field. The maser region results not only limit the properties of turbulentfluctuations in the Galactic field, but also restrict the field randomization that presumably occurs in the assemblyof star formation cores during the process of star formation.4. We point out that studies of turbulence in laboratory settings also show a persistence of large scale asymmetriesin small and inertial scale fluctuations, which could be termed “the persistence of memory” in the Galacticmagnetic field. These results are discussed in Warhaft (2000), and refer to passive scalars in hydrodynamicturbulence. Although it is unclear if the results discussed in Warhaft (2000) are applicable to Galactic MHDturbulence, it would appear to be a worthwhile path to investigate.2 REFERENCES
Beck, R. 2016, A&A Rv, 24, 4Beckmann, P. 1967,
Probability in CommunicationEngineering , Harcourt, Brace & World, p36Costa, A.H., Spangler, S.R., Sink, J.R. et al. 2016, ApJ,821, 92Costa, A.H. and Spangler, S.R. 2018, ApJ, 865, 65Crutcher, R.M. and Kemball, A.J. 2019, Frontiers inAstronomy and Space Science, 6, 66Ferri`ere, K. 2011, Mem. Soc. Astron. Italiana, 82, 824Ferri`ere, K. 2020, Plasma Physics and Controlled Fusion,62, 014014Fish, V.L., Reid, M.J., Argon, A.L., and Menten, K.M.2003, ApJ, 596, 328 Green, J.A., McClure-Griffiths, N.M., Caswell, J.L. et al.2012, MNRAS, 425, 2530Han, J.L. 2017, ARA&A, 55, 111Harvey-Smith, L., Madsen, G.J., and Gaensler, B.M. 2011,ApJ, 736, 83Minter, A.H. and Spangler, S.R. 1996, ApJ, 458, 194Nicholson, D.R.. 1983,
Introduction to Plasma Theory ,John Wiley & Sons, p3Savage, A.H., Spangler, S.R., and Fischer, P.D. 2013, ApJ,765, 42van Eck, C.L., Brown, J.C., Stil, J.M. et al. 2011, ApJ, 728,97Warhaft, Z. 2000, Annual Review of Fluid Mechanics, 32,203 A. GLOSSARY OF MATHEMATICAL VARIABLES
Variable Definition∆ χ change in polarization position angle due to Faraday rotation (Eq.1) ~B total magnetic field in ISM at location of HII region (Eq.2) ~B large scale, ordered component of Galactic magnetic field ~δb turbulent component of Galactic magnetic fieldΘ angle between LOS and large scale magnetic field ~B at location of HII region (Fig. 2)(ˆ e k , ˆ e p , ˆ e ⊥ ) unit vectors defining coordinate system aligned with large scale field ~B ( δb k , δb p , δb ⊥ ) components of the turbulent magnetic field in a coordinate system aligned with large scale field ~B RM obs observed RM for a background source observed through an HII region (Eq.3) RM bck RM due to ISM along LOS, but excluding HII region (Eq.3) RM neb RM due to plasma in HII region (Eq.3) C set of atomic constants appearing in expression for Faraday rotation (Eq. 1) n electron density in plasma causing Faraday rotation (Eq. 1) R radius of HII region B z LOS component of Galactic magnetic field at location of HII region A factor ( ≥
1) determining whether Galactic field is amplified within HII region. Model dependent (Eq. 4)ˆ e z unit vector in direction of LOS x dimensionless amplitude of turbulent magnetic fluctuations in the direction of the large scale field (Eq. 8) y dimensionless amplitude of turbulent magnetic fluctuations ⊥ to large scale field (Eq. 9) ξ derived variable dependent on Θ, x , and y (Eq. 12)∆ threshold value of ξ such that sign of RM reversed for ξ < − ∆ σ rms value of Gaussian-distributed y (Eq. 14) a secondary variable a ≡ cos Θ (Eq. 16) b secondary variable b ≡ sin Θ (Eq. 16) d secondary variable d ≡ σb (Eq. 17) p ( ξ ) normalized, Gaussian pdf of ξ (Eq. 15-17) p G ( ξ ) shifted, zero-mean version of p ( ξ ) (Eq. 19) P − probability that turbulent fluctuations cause reversal of LOS magnetic polarity (Eq. 18) t secondary variable t ≡ y √ d (Eq. 17) X single parameter determining P −−