Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.izmiran.rssi.ru/~obridko/papers/171.pdf Дата изменения: Tue Mar 22 14:24:43 2011 Дата индексирования: Sun Apr 10 01:59:37 2016 Кодировка: Поисковые слова: cosmic rays

FOURIER PARAMETERS AND MOMENTS POLARIZATION PROFILES OF MAGNETICALLY LINES. FOURIER VECTOR

OF ACTIVE

MAGNETOGRAPH

B. IOSHPA and V. OBRIDKO
IZMIRAN, Troitsk, Moscow Region, Russia

and I. KOZHEVATOV
NIRFI, Nizhny Novgorod, Russia

Abstract. A new method is proposed to determine all components of the solar magnetic fields using the cumulants of the profile of a magnetic sensitive line. The method is based on polarization measurements in a number of points of the line profile and subsequent calculation of the amplitudes and phases of its two first Fourier-harmonics. Key words: Solar magnetic field - Fourier-magnetograph

1.

Introduction

The methods employed to obtain vector magnetic fields on the Sun are usually based on measurements of the Stokes parameters in magnetoactive lines followed by interpretation with adequate radiative transfer modeling. As a matter of fact, one measures the polarized intensities either in selected regions of the profile (using Babcock-type magnetograph), or all over the profile (using Stokesmeter). The measurements can be carried out simultaneously or by scanning along the profile. The magnetic field parameters, calculated from the measured polarization, are very sensitive to variations of the line profile due to the physical conditions in the line forming layer (such as temperature, pressure, and inhomogeneities), and strongly depend on the adopted model atmosphere. The difficulties involved in the methods which employ a single wavelength can be avoided by using integral parameters of the profile (displacement of the gravity center, width, asymmetry etc.). Semel (1970) was one of the first to use relative displacements of the gravity centers of the a-components of magnetoactive lines as a measure of the longitudinal magnetic field. The present work is, in some sense, an extension of this idea to the vector magnetic field.
Solar Physics 164: 373-380, 1996. 9 1996 Kluwer Academic Publishers. Printed in Belgium.


374
2.

B. IOSHPA

ET AL.

Cumulants and their relations to magnetic fields

As integral parameters of the line we shall take the cumulants determined from the following relations:

mn = j-n "dnlnI(w) ]

Jw=0'

(1)

where I(w) = f I(.k)eJ~;~d)~ is the Fourier transform of the line depression profile I(,k), and j = x/Q--f. There exist relationships between the cumulants and the central and starting moments of I (Malakhov, 1978): 8s1 ~-~ OL1 - center of gravity of the line profile se2 = #2 = c~2 - c~ - square of the width of the line profile se3 -- P3 -- oz3 - 3oqo~2 + 2ol3 - asymmetry of the line profile where f·
]~n

(~ - ~1) n d~

fI()~)d()~)

'

are the central moments,

~:

f I(A) An dA f I(,~)d()~)

'

are the starting moments, and I(A) is the line depression profile as a function of the wavelength measured from the nominal value. Though the first cumulants coincide with the central moments, we shall rather use the system of cumulants owing to their additivity in the convolution procedure (Malakhov, 1978). As shown below, the knowledge of the two first cumulants of different states of polarization is sufficient to calculate all components of the vector magnetic field. Let us discuss one of the possible procedures. Let solar emission in a magnetoactive line with Stokes parameters Io,Vo,Qo and U0 be analyzed by a polarization analyzer consisting of 2 elements: a controlled birefringent plate and a linear polarizer. The operation of the analyzer will be described in the reference frame in which the axes coincide with those of the polarization ellipse, consequently U0 = 0 *. The angle between the transverse magnetic field component and the plate axis is denoted by fl, and the angle between the beam and the magnetic field vector by %
* In general, the polarization ellipse does not maintain a single orientation along the line profile and, consequently, such a single reference system cannot be found. However, some solutions of the radiative transfer equation do allow such a definition, in particular, those used at the end of this section.


FOURIER VECTOR MAGNETOGRAPH

3"75

The polarization analyzer has six states: 1) G (88 P (/3 § 45 ~ - delay between the ordinary and extraordinary beams produced by the plate is ~, angle between polarizer axis and the fast axis of the plate is 45 ~ The arguments of G and P describe state of the retarder and the polarizer, respectively. 2) G (-88 P (fl + 45 ~ - delay is -}, angles are the same. 3) G (0, fl); P (/3 § 45 ~ - delay is zero, angles are the same. 4) G (0,/3); P (/3) - delay is zero, polarizer axis turned by -45 ~ 5) G (0,/3); P (/3- 45 ~ - delay is zero, polarizer axis turned by -90 ~ 6) G (0,/~); P (/3-90 ~ - delay is zero, polarizer axis turned by -135 ~ It is easy to show that radiation intensities at the analyzer output for each of its 6 states can be described as follows:

IO) = 0.5Io - 0.5Vo f (2) = 0.5Io + 0.5Vo

i(3) = 0.5Io - 0.5 sin(2/3)Qo

(2)
I (4) = 0.5Io + 0.Scos(2/~)Qo i(5) = 0.5Io + 0.5sin(2/3)Qo
I (6) = o.5Io -

0.5 cos(2 )Qo.

Now, let us calculate the relationships between the original Stokes parameters and the cumulants for six states of the analyzer. For the sake of simplicity, assume that I0 and Q0 are strictly even functions of ,~ and ~ is strictly odd function of ~. By substituting I(~) -- i(6) in the expressions for moments we obtain:

=

=

-A(H,

cos

(3a)


376

B. IOSHPA ET AL.

~e~ = ~e~ = B(H, 3') - A2( H, 7)" cos2 3' 1) 2) ~e~3) =
~e(4)2 -

B - Csin 2/3 1-Dsin2/3
B + C cos 2/3

-- l+Dcos2/3 ~e~5) _ B + C sin 2fl l+Dsin2/3
-

(3b)

@)
where

_

B - C cos 2r

-

1 - D cos 2/3'

A(H,7) = f Vo),d)~. 1 f Iod~ cos 7'
B(H, 3`) - f I~ f lod), '

(4)

(5)

C(H, 3`) - f Q~ : Iod~ ' D(H, 7) - f Qod,~ f Iod~, '

(6) (7)

and H is the magnetic field strength. However, the signals under consideration depend on properties of the atmosphere (pressure, temperature, velocity and magnetic fields, and the variation of these parameters along the line-ofsight). In order to simplify the computation, we adopt the same approximation as Unno (1956). Under this approximation, A, B, C, and D, depend on constant values of H (the field strength), 7 (the inclination angle), ~0 (the line-to-continuum opacity ratio), g (the Lande factor), A),D (the Doppler halfwidth of the line) and the two coefficients which describe the source function. Figure 1 shows calculations based on this approach. The parameters have been chosen to correspond to those of the photospheric Fe I 5250 A line. The ratio r]0 of opacities at the line center to that in the continuum was taken equal to 10. The field strength is given in units of the Doppler halfwidth (H = 1 corresponds approximately to 1000 G). A is expressed in Doppler halfwidths, B and C in square of the Doppler halfwidths.


FOURIER VECTOR MAGNETOGRAPH

377

3.0

9

12.0 -

8.0

"

21 io., ~ I
H 6.0
I'

4.0

~176 ....... f.b ...... Y.b ...... :3.0
H

0.30
C(H,7)
I / /

,,

D(H,7.).-i

........

4-.0
,
I

,"/"
/f/
~"

0.15
,
II r

,'"
,..

H
Fig. 1. line FeI 1000 G. line for

...... X

H

3

Dependence of A, B, C, and D on the field strength H. The simulations refer to the 5250 · observed in the photosphere, so that H = 1 corresponds approximately to The solid line is used for "7 = 0, the dashed line for 3' = 7r/2, and the dash-dotted h' = 77/4.

3. Calculating cumulants by Fourier-analysis. Fourier vector magnetograph

Below we describe the method for calculating the cumulants by Fourier analysis of the line profile. For this purpose, we use expressions connecting the amplitudes and phases of the Fourier harmonics of the line profile with its cumulants (Didkovsky, Kozhevatov and Stepanyan, 1986):


378

B. IOSHPA ET AL.

(zxa)

ael

_

:

1,al + ..., (8)

r.,
s=

:

i(a)ea,

where A,k = 2rr/w is the period of the Fourier harmonic, and .4 and 9 are its amplitude and phase, respectively. The two first eumulants are enough to determine the magnetic field under our simplifying assumptions (see Section 2). However, the determination of the first three cumulants might allow a more elaborate modeling (for instance, taking into account the asymmetry of the line profiles). Because of this future possibility, and since the determination of two or three harmonics requires a similar computational effort, here we describe how to compute r r and r To determine the three first cumulants, we need to know the amplitudes and phases of the two first harmonics of the Fourier-transform of the line profile (the domain of definition of the line profile is considered as the period of the first harmonic). As seen from Equation (8), the phases of the long-period harmonics (those with small 1/A~) are proportional to the shift of the center of gravity (ml). In order to determine a~t and ae3, we have to know the phases of two harmonics (with small 1/A)~1 and with somewhat larger 1/A~2). A similar procedure is applied to determine the even cumulants from the amplitudes. The period of the first harmonic is large enough, so that the first term dominates the series for the phases. The term A)~-3 and higher order terms can be neglected, and consequently, 51 is determined from the phase of the first harmonic. The second harmonic has a smaller period; in this case the series for the phases has two major terms while the series for the amplitudes has just one term. Since the first cumulant and S are known, the expressions for the phase and amplitude of the second harmonic are used to derive the third and the second cumulants. According to the Kotelnikov - Shannon criterion (see Bell, 1972), the characteristics of the Fourier harmonics can be determined by making records at discrete points spaced by no more than 1/3 of the harmonic period. Therefore, the intensities measured at six points along the profile are enough to determine the amplitude and phase of the two first harmonics. In this case one condition should be observed: before discretization, higher-frequency harmonics must be eliminated from the line profile to avoid aliasing (see, e.g., Bracewell, 1978). The smoothing can be done by selecting the appropriate resolution of the instrument. The procedure of determining magnetic fields by the proposed Fouriermethod comprises the following stages:


FOURIER VECTOR MAGNETOGRAPH

379

1. Measure intensities at six discrete equidistant points along the line profile at six positions of the polarization analyzer. 2. Calculate the phases and amplitudes of the two first harmonics of the line profile. 3. Use these to find two or three first cumulants of the line profile. 4. Determine the magnetic field parameters from Equations (3)-(7). These expressions yield at least 6 independent equations for 3 unknown values, H, 7, and/3. Consequently, one expects no unique solution and, therefore, some kind of least squares best fit has to be sought. This procedure has not been developed yet. Note that parameters of the Fourier-harmonics can be obtained by a different method than proposed in items 1 and 2 above, for example, directly by using Michelson interferometers (Bell, 1972; Didkovsky, Kozhevatov and Stepanyan, 1986).

4.

Conclusions

We believe that the proposed method for determining magnetic field parameters from the cumulants of the line profile offers a number of advantages, compared with the traditional Babcock method, based on the comparison of polarized light intensities at two symmetric points in the line profile. 1. The relation between the cumulants and the magnetic field value does not show saturation or back pass at large field intensities. 2. The magnetic field is determined from integral properties of the line profile, so that the result will not be affected strongly, by the asymmetries of the profile. 3. The use of integral properties will considerably reduce the effect of spatial irregularities (see Rees and Semel, 1979). 4. When determining the magnetic field with a Babcock-type filter magnetograph, the presence of non-compensated line-of-sight velocity may result in significant errors. Our method, in which the effects of the magnetic field and the line of-sight velocity are separated, does not have this disadvantage, and to a certain extent, it resembles the photographic method. (For detailed description of the photographic method see, for example, Bray and Loughead, 1964) 5. Our method requires measurements at a limited number of points along the line profile.

Acknowledgements
We thank Jorge Sgnchez Almeida for valuable comments on the manuscript.


380

B. IOSHPA ET AL. References

Bell, J.: 1972, Introductory Fourier Transform Spectroscopy, Academic Press, New York Bracewell, R.N.: 1978, The Fourier Transform and its Application, McGraw Hill, New York Bray, R and Loughead, R.: 1964, Sunspots, Chapman and Hall Ltd., London Didkovsky, L.V., Kozhevatov, I.E., and Stepanyan, N.N.: 1986, Izvestia Krimskoy Astrophysicheskoy Observatorii (SU) 74, 142 Malakhov, A.N.: 1978, Cumulant analysis of Random non-Gaussian Processes and their Transformation, Sov. Radio, Moscow Rees, D.E., and Semel, M.D.: 1979, Astron. Astrophys. 74, 5 Semel, M.D.: 1970, Astron. Astrophys. 51 330 Unno, W.: 1956, Publ.Astron.Soc.Japan 8, 108