A&A 412, 875-878 (2003)
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
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 (
is the continuum
intensity) of any spectral line are square integrable functions of
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.
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.