A&A 412, 875-878 (2003)
DOI: 10.1051/0004-6361:20034429
Research Note
J. C. del Toro Iniesta1 -
A. López Ariste2,
1 - Instituto de Astrofísica de Andalucía (CSIC),
Apdo. de Correos 3004, 18080 Granada, Spain
2 -
High Altitude Observatory, NCAR, PO Box 3000, Boulder, CO
80307, USA
Received 2 October 2003 / Accepted 21 October 2003
Abstract
A family of well-known orthonormal functions, the set of
Hermite functions, is proposed as a suitable basis for expanding the
Stokes profiles of any spectral line. An expansion series thus
provides different degrees of approximation to the Stokes spectrum,
depending on the number of basis elements used (or on the number of
coefficients). Hence, an usually large number of wavelength samples,
may be substituted by a few such coefficients, thus reducing
considerably the size of data files and the analysis of observable
information. Moreover, since the set of Hermite functions is an
universal basis, it promises to help in modern inversion techniques
of the radiative transfer equation that infer the solar physical
quantities from previously compiled look-up tables or artificial
neural networks. These features appear to be particularly important
in modern solar applications producing huge amounts of
spectropolarimetric data and on near-future, on-line
applications aboard spacecrafts.
Key words: line: profiles - magnetic fields - Sun: magnetic fields - stars: magnetic fields - radiative transfer
Most of what we currently know about the Sun's magnetic fields is
provided to us by the Stokes profiles of spectral lines (Lites et al. 1994).
Spectropolarimetry appears to be the most useful among the current
means devised to observe the solar and stellar atmospheres (del Toro Iniesta 2003a) and, certainly, inversion techniques
(ITs) of the radiative transfer equation reveal themselves as the
more suitable techniques devoted to infer the atmospheric physical
quantities from the observables. Recent advances on instrumentation
are challenging our abilities to deal with the huge amounts of data
produced and some improvements or modifications of our currently
available ITs will be welcome as well as any speed increase in
numerically solving the radiative transfer equation, a procedure in
which most ITs lay somehow. The main spectropolarimetric observable
is the set of four Stokes parameters as functions of the wavelength:
the Stokes profiles. Every inference technique pretends to reproduce
the observable by means of some physical model of the stellar
atmosphere and some known physics of transport of
radiation. Certainly, if we are able to isolate repetitive features
in the profiles and represent them as a linear superposition of these
basic features, the analysis simplify and the computation time
may reduce significantly as compared to other, more classical,
inversion procedures (for reviews, see Socas-Navarro
2001; del Toro Iniesta 2003b; Bellot Rubio 2003). An expansion series in terms of
basis profiles is the fundamental of the Principal Component
Analysis (PCA, Rees et al. 2000). This method has proved to
produce reliable results (Socas-Navarro et al. 2001) and its coefficients to have relations with the solar
atmospheric parameters (Skumanich & López Ariste 2002)
but, from a practice point of view, it shows the main drawback of
having a non-unique set of basic profiles: this set must be calculated
once for every observation or one must acknowledge the risk of using
basic profiles that do not reproduce every particular observation.
Here we propose a new means for Stokes profile data management by
introducing a well-known universal basis for the vector space to
which any Stokes profile belongs:
,
the space of square
integrable functions. (More specifically, all three Stokes Q,
U, and V, and Stokes I in line depression.) By means of this
basis, several wavelength samples can be substituted by a few
expansion coefficients, thus reducing the dimensionality of the
inversion problem.
![]() |
Figure 1: First six Hermite functions. |
Open with DEXTER |
The structure of the paper is as follows: we introduce the set of Stokes basis functions in Sect. 2. In Sect. 3, a numerical experiment with a data base of more than 40 000 synthetic Stokes profiles is carried out in order to illustrate the various degrees of approximation depending on the number of basis profiles used. A fit to a real (i.e. observed) Stokes profile set is also presented in that section. Finally, the results are discussed in Sect. 4.
Considered as functions of the wavelength, all three Stokes Q, U,
and V profiles and Stokes I in line depression (
,
where
is the continuum
intensity) of any spectral line are square integrable functions of
wavelength:
.
This
property stems from the fact that all of them are continuous and drop
asymptotically to zero. As is well known,
is a Hilbert
space where a scalar product,
![]() |
Figure 2: Average root-mean-square deviations between Milne-Eddington profiles and their successive expansions in terms of Hermite functions. The number of terms is in the abscissae. The error bars represent the standard deviation among the 46 844 expansions. |
Open with DEXTER |
A glance to Fig. 1 readily illustrates our interest
in this set of functions:
has the shape of
,
looks like a typical Stokes V profile due to the Zeeman effect, and
is similar to
Stokes Q or U. Certainly, these three first elements of the
orthonormal basis have definite symmetry properties and would
represent the Zeeman Stokes spectrum of a line formed in the absence
of velocity gradients (see, e.g., Auer & Heasley 1978; Landi Degl'Innocenti & Landi Degl'Innocenti 1981). The
asymmetries present in real Stokes profiles of such
lines, formed in the presence of gradients, can be taken into
account through the contribution of higher order functions. Since the
set of Hermite functions is a basis of the
space, every
Stokes profile can be represented exactly by an infinite linear
combination of hn functions. The strong similarities between the
first Hermite functions and the actually observed profiles suggests
that just a few such functions are needed for an expansion of the
Stokes parameters as functions of the wavelength. If this is the
case, we have the means for substituting the (usually) many
wavelength samples of the observations by a few coefficients with the
consequent reduction in dimensionality. Note that an expansion of the
profiles as the one we propose is independent of any radiative
transfer physics and, therefore, the results can later be processed
with any other analysis technique that infers the atmospheric
physical quantities.
![]() |
Figure 3: Observed Stokes profiles (circles), and Hermite expansions up to n=10 (dotted lines) and up to n=19 (solid lines), of the Fe I line at 1558.8257 nm as formed in a sunspot penumbra. |
Open with DEXTER |
The calculation of the expansion coefficients is very simple and
based on the scalar product (1). In fact, if a generic
Stokes profile is denoted by
and is to be expanded to n Hermite functions,
Not only symmetric Stokes profiles like ME profiles can be
represented with a few Hermite functions. As an example, we show in
Fig. 3 two fits, with 11 (n=10; dotted line) and
20 terms (n=19; solid line), to real Stokes profiles observed in
the penumbra of a sunspot where del Toro Iniesta et al. (2001) discovered supersonic downflows. The rms differences
between observations and eleven-term fits are
,
,
,
and
for I, Q, U,
and V, respectively, clearly smaller than or equal to the noise
levels. The differences with the twenty-term fits are even slightly
smaller. This observation of the Fe I line at 1558.8257 nm has
been selected as an extreme example of asymmetry.
We have shown the properties of the Hermite family of functions as an orthonormal set of Stokes profiles. Every Stokes profiles, symmetric or not, can be fitted with a few Hermite functions, hence providing an excellent means of reducing the dimensionality of data representation. As compared to other Stokes expansions like that performed within the PCA analysis, the advantages of the Hermite Stokes basis profiles is twofold. On the one hand, the basis is universal and does not depend on the particular observation; we know it a priori and can even use this knowledge in on-line applications. On the other hand, the orthonormal basis is analytical and thus noise-free. These features allow us to consider the Hermite expansion of the profiles as an alternative representation of the observations and enable its use in whatever inversion technique of the radiative transfer equation, either current or future.
We can simply consider Hermite expansion as a tool for reducing the data dimensionality and use it, for example, to circumvent limited telemetry problems from spacecrafts. Or we can use the analytic approximations as to build the look-up tables needed in PCA-like inversion techniques, or even to train neural network machines (Carroll & Staude 2001; Socas-Navarro 2003). In both PCA and neural network cases, an acceleration of the procedures is guaranteed since we have to deal with much less observables: the set of Hermite coefficients. But the finite Hermite representation can even be used as an input to "regular'' inversion techniques like ME or SIR. Therefore, we are facing an extremely useful means in practice.
Besides these practical features, the mathematical properties of Hermite functions may in principle be used with much more theoretical purposes. In fact, the possibilities must be explored that these Hermite functions help in solving the radiative transfer equation in a more efficient and accurate way than those currently used.
Acknowledgements
This work was supported by the Spanish Programa Nacional de Astronomía y Astrofísica through Project AyA2001-1177 and Programa Nacional del Espacio through Project ESP2002-04256-C04-01, partly using European FEDER funds. J. Graham (HAO) kindly provided us with the set of ME profiles from ASP inversions used in Sect. 3.