A&A 475, 349-357 (2007)
DOI: 10.1051/0004-6361:20078046
A. I. Shapiro^{1} - S. V. Berdyugina^{1,2} - D. M. Fluri^{1} - J. O. Stenflo^{1,3}
1 - Institute of Astronomy, ETH Zurich, 8092 Zurich, Switzerland
2 - Tuorla Observatory, University of Turku, 21500 Piikki
,
Finland
3 - Faculty of Mathematics & Science, University of Zurich, 8057 Zurich, Switzerland
Received 8 June 2007 / Accepted 4 September 2007
Abstract
Context. Weak entangled magnetic fields with mixed polarity occupy the main part of the quiet Sun. The Zeeman effect diagnostics fails to measure such fields because of cancellation in circular polarization. However, the Hanle effect diagnostics, accessible through the second solar spectrum, provides us with a very sensitive tool for studying the distribution of weak magnetic fields on the Sun.
Aims. Molecular lines are very strong and even dominate in some regions of the second solar spectrum. The CN
system is one of the richest and most promising systems for molecular diagnostics and well suited for the application of the differential Hanle effect method. The aim is to interpret observations of the CN
system using the Hanle effect and to obtain an estimation of the magnetic field strength.
Methods. We assume that the CN molecular layer is situated above the region where the continuum radiation is formed and employ the single-scattering approximation. Together with the Hanle effect theory this provides us with a model that can diagnose turbulent magnetic fields.
Results. We have succeeded in fitting modeled CN lines in several regions of the second solar spectrum to observations and obtained a magnetic field strength in the range from 10-30 G in the upper solar photosphere depending on the considered lines.
Key words: line: formation - Sun: magnetic fields - molecular processes - polarization - radiative transfer - scattering
These hidden magnetic fields can be accessed with the Hanle effect (Stenflo 1982), which modifies the properties of coherent scattering processes. The Hanle effect can thus only be employed through the spectrum formed by the coherent scattering, which is called the "second solar spectrum'' (Ivanov 1991) due to its richness and significance. The Hanle effect usually leads to a reduction and rotation of line polarization as compared to the non-magnetic case. The net-influence of rotation cancels out in a turbulent magnetic field, but the depolarization remains visible for most geometries even in a turbulent magnetic field. Moreover, the Hanle effect is extremely sensitive to the weak field. Therefore, the second solar spectrum provides us with a very sensitive tool for studying the distribution of weak magnetic fields on the Sun and gives the possibility to measure spatially unresolved mixed polarity magnetic fields.
Molecular lines are particularly valuable for magnetic field diagnostics thanks to their broad range of magnetic sensitivities within narrow spectral regions. This allows us to employ the differential Hanle effect (Stenflo et al. 1998), which greatly reduces the model dependence of deduced magnetic field strengths, and has lead to the first unambiguous detection of molecular Hanle effect in the second solar spectrum (Fluri & Berdyugina 2004). Recent studies of the molecular Hanle effect have yielded field strengths of the order of 10 G for the quiet solar photosphere, assuming a single-valued, isotropic, turbulent magnetic field (Berdyugina & Fluri 2004; Trujillo Bueno et al. 2004; Faurobert & Arnaud 2003), although larger field strengths are possible depending on the assumed collision rate, in particular when interpreting MgH lines (Bommier et al. 2006; Asensio Ramos & Trujillo Bueno 2005).
Diagnostics with molecular Hanle effect has concentrated so far only on two molecules, and MgH. It is very important to extend the technique to other molecules because in general different molecules form at different heights due to their large temperature sensitivity and can thus sample other layers of the solar atmosphere. Combined with atomic Hanle effect, which is well established since the work by Faurobert-Scholl (1992), only such a multitude of lines will allow us to constrain the 3-dimensional structure of the unresolved solar magnetic field.
In the present paper we apply the Hanle effect to a new molecule, namely the CN violet system, for which we have recently developed a complete theory of the Hanle effect taking into account the Paschen-Back effect and interference between doublet components (Shapiro et al. 2007). After describing our radiative transfer method we fit our model to observations in three spectral regions of the second solar spectrum, which enables us to determine the magnetic field strength in the quiet solar photosphere.
In this section we introduce our method for solving the radiative transfer problem. It is based on the modeling strategy developed by Faurobert & Arnaud (2003), who have calculated scattering polarization at the line centers of and MgH transitions. Here, we expand their method to obtain wavelength dependent spectra not restricted to line centers.
The radiative transfer model is optimized for relatively weak lines whose scattering polarization is calculated within the framework of a plane-parallel atmosphere using the single-scattering approximation. The CN lines are assumed to form in a homogeneous atmospheric layer situated above the continuum formation region. We fully account for line blending and for depolarization of the continuum polarization by line absorption and scattering. In this respect, the current model is more complex than the one employed by Berdyugina & Fluri (2004) for very weak lines. The CN bands that are modeled in this paper include stronger lines and require a consistent treatment of continuum depolarization, although we increase at the same time the model dependence of the magnetic field diagnostics, because we rely on a proper computation of the continuum polarization. Nonetheless, the model remains simple and ideally suited for differential Hanle effect diagnostics.
First we introduce the necessary notation and the source function. Then we discuss in detail the two-layer approximation, the free model parameters, and the fitting procedure for comparing the model with observations.
The total, monochromatic optical depth is given by
Now we can write the radiative transfer equation in the usual form
In our further discussion we employ the single-scattering approximation, since the CN lines for which we will apply this theory are optically thin. Then, the scattering part of the line source function takes the form
Within the weak-field approximation and under the assumption of an isotropic turbulent magnetic field the phase matrix
is parametrized as
With our choice of the phase matrix we implicitly neglect interferences between different fine structure components and perturbations by the Paschen-Back effect. The latter becomes relevant in the CN violet system only for field strengths exceeding 100 G (Shapiro et al. 2007), while the fields detected in the present study are weaker than 30 G.
Let us rewrite Eq. (4) in the form
The emergent Stokes vector is given by the formal solution of the radiative transfer equation
The "continuum'' layer corresponds to the lower photosphere of the average quiet Sun, defines the lower boundary condition of the CN layer, and is thus crucial for the formation of the CN lines. In particular, the continuum layer determines the anisotropy of the continuum radiation, without which scattering polarization would cancel out, and it sets the degree of continuum polarization, which is subsequently depolarized by the line opacity.
In our modeling strategy we take the continuum intensity and polarization from previous work, rather than recalculating these values. The continuum intensity is obtained from the analytical expressions for limb-darkening derived from observations of the quiet Sun (Neckel 1996). We assume the radiation field to be azimuthally symmetric. The values of the continuum polarization correspond to the calculations by Stenflo (2005) for . The continuum polarization for arbitrary -values are computed by scaling the value at with the analytical expressions introduced by Fluri & Stenflo (1999), which describe the -dependence of the continuum polarization and are based on theoretical modeling of the quiet solar photosphere.
The CN layer is described by the source function and the radiative transfer equation given in Sect. 2.1. In Eqs. (5) and (12) we insert the Stokes vector , which is incident into the CN layer at the lower boundary and given by the frequency and angle dependent continuum intensity and polarization formed in the continuum layer. The same equations require also the total optical thickness of the CN layer. Assuming , since the continuum is formed below the CN layer, the total optical thickness of the CN layer can be calculated with the expression
where the frequency independent coefficient is a free parameter of the model. It depends on the CN density and the geometrical thickness of the CN layer.
Within the two-layer approximation the total source function, Eq. (10), is further simplified since the opacity in the CN layer is only due to CN molecules, i.e. . We assume that the ratio of the scattering and absorption probabilities is the same for all molecular lines. This is basically equivalent to the assumption that the inelastic collision rate is independent of the involved transitions because the life time of the exited CN state is almost independent of its quantum numbers (with the exception of some states with very small J numbers). Furthermore we neglect in our model the elastic collisions. Equation (10) is then reduced to the form
Here the first term with the summation represents the scattered diffuse part of the radiation field, restricted within the CN layer to line contributions, while the second term corresponds to the thermal part.
The emergent Stokes vector is calculated employing the formal solution of the radiative transfer equation. To improve the accuracy of our results we divide the CN layer into k sublayers with equal optical thickness. The source function is assumed to be constant within each sublayer and evaluated at the center of each layer, giving us k different values . Introducing the above assumptions into Eq. (12), we find for the emergent radiation field at a given frequency
The first term describes the radiation transmitted through the CN layer while the summation in the second term adds up the contributions from the different CN sublayers to the emergent spectrum. Note that the depth dependence of the source function within the CN layer is introduced by the exponential factor in Eq. (5), which describes absorption of the incident radiation and thus the height variation of the radiation field even within the isothermal slab. The scattering coefficient, on the other hand, depends only on the probability of the collisional de-excitation of the upper CN states, which under our assumptions remains constant for the whole molecular layer.
The subdivision of the molecular layer, expressed by Eq. (15), is a further sophistication of the model introduced by Faurobert & Arnaud (2003), who have assumed one single value of the source function within the molecular layer. The sublayers take care of the depth dependence of the source function within the CN layer, which results from the optical depth in the exponent of Eq. (5). Note that this does not contradict to our assumption of an isothermal CN layer with constant density. Our test calculations have shown that the emergent Stokes vector becomes independent of the number of sublayers for k greater than five.
The model described above contains several free parameters, which have to be determined by fitting the theoretical spectrum to observations. These are the temperature of the isothermal CN layer , the scattering probability , the optical depth normalization coefficient , and the magnetic field strength B.
The temperature affects the emergent spectrum in two ways. First of all, it enters via thermal term of the source function (cf. Eq. (14)). This strongly influences the intensity, but not the other Stokes parameters, because of the strong temperature dependence of the Planck function and because the contribution from the thermal part is relatively large in the case of CN lines. Secondly, the line opacities depend on temperature via the Boltzmann factor (cf. Eq. 2), which affects the whole Stokes vector, although in a much weaker way than the thermal term.
The scattering probability defines the balance between the thermal and nonthermal parts of the total source function and influences the whole Stokes vector. The contribution to the Stokes vector from the CN layer linearly depends on .
The optical depth normalization coefficient describes the relative contributions to the Stokes vector from the CN layer and the continuum layer. It effectively determines the optical thickness of the CN layer, thus scaling both the molecular line depths in the intensity spectrum and the degree of line polarization.
Our main goal is the determination of the magnetic field strength B assuming a turbulent, isotropic field distribution within the CN layer. The magnetic field enters via Eq. (8) into the line source function. Since we consider weak fields and neglect Zeeman broadening, the magnetic field strength has virtually no effect on Stokes I and only causes differential Hanle effect in the second solar spectrum (see Sect. 3.1).
Figure 1: Differential Hanle effect in synthetic Stokes profiles of R-branch lines in the CN violet (0,0) band. For selected angular momenta N of the lower state we show the complete doublet, which causes the double peak structure for large N. Three values for the turbulent magnetic field strength were chosen, giving identical Stokes I profiles but large modifications in the Stokes Q profiles: 0 G (upper curve), 10 G (intermediate curve), and 100 G (lower curve). In this and all other figures we assume an elastic collision rate equal to zero. | |
Open with DEXTER |
Usually, the influence of B on the spectrum can be easily distinguished from the effects of and . The turbulent magnetic field causes a depolarization of the CN lines. The same can in principle also be obtained with a smaller scattering probability . Even an increased parameter can mimic the same behavior because it leads to stronger continuum depolarization. The crucial point, however, is that and affect all molecular lines with comparable strength, while the magnetic field leads to differential depolarization, because the Hanle effect efficiency strongly depends on the J quantum number (Fig. 1). Thus, the relative strength of CN lines in the second solar spectrum mainly depends on the magnetic field, but not so much on the other free parameters. Therefore, it is possible to diagnose the magnetic field strength when fitting several lines with different Hanle effect sensitivities.
We simultaneously fit the modeled polarization spectrum to observed
data and the modeled intensity spectrum to
observations. During the fitting procedure we vary the magnetic field strength B, the scattering probability
,
the effective temperature
and the optical depth normalization coefficient
.
The most suitable values for these free parameters are determined with a least-square method. More precisely, we minimize the quantity
The value in Eq. (16) is a measure of the discrepancy between observations and the theoretical model, and describes the accuracy of the fit. Wavelength points that are heavily affected by atomic blends are excluded from the summation and the fitting procedure since we neglect atomic line opacity in our model.
In this section we apply the previously developed theory to the interpretation of observations of the CN system, the so-called violet system, in several spectral regions, with the goal to determine the turbulent magnetic field strength in the quiet Sun. From the discussion in Sect. 2.3 it is clear that suitable wavelength regions should contain molecular lines with different Hanle effect sensitivities. Therefore, we first investigate the differential Hanle effect in the CN violet system. This will guide us to choose three ideal spectral windows in different molecular bands, which are then analyzed and discussed.
Figure 1 illustrates the different Hanle effect sensitivities of synthetic line profiles for the example of the CN violet (0,0) band (other CN violet bands give the same qualitative results). The relative degree of line polarization in the non-magnetic case is defined by three factors: the oscillator strength (which increases with the J number), the scattering polarizability W_{2} (which decreases with J and asymptotically approaches a value of 0.1), and the Boltzmann factor (which also decreases with J). However, the final polarization obtained within the lines not only depends on the intrinsic line polarization. Rather it results from a combination of different contributions, namely the continuum polarization, continuum depolarization due to the line opacity, and line polarization, where only the latter is subject to the Hanle effect. It is clearly visible in Fig. 1 that stronger lines lead to additional continuum depolarization (cf. Fluri & Stenflo 2003). Furthermore, we can identify the expected diverse behavior of the lines in a magnetic field. While the intensity profiles remain unaffected, the line polarization is depolarized with increasing field strength, reflecting the J dependence of . A magnetic field of 10 G already leads to a nearly saturated Hanle effect in the lines with and 5, while the polarization in the and 80 lines remains almost unchanged compared to the non-magnetic case. The latter lines require a magnetic field strength in the order of 100 G for a considerable depolarization.
Figure 2: The dependence of the Hanle depolarization factor , given in Eq. (7), on the magnetic field strength and the N quantum number for the CN violet (0,0) band. The three curves correspond to the depolarization factors equal to 0.9 (solid), 0.5 (dashed), and 0.3 (dotted), the last of which is already close to the asymptotic saturation limit of 0.2. | |
Open with DEXTER |
The magnetic field necessary for sizable depolarization by the Hanle effect is shown in Fig. 2 for different N quantum numbers, again for the CN violet (0,0) band. The different curves correspond to three values of and identify the magnetic field range in which a specific line is sensitive to the Hanle effect. For example, lines with react to the Hanle effect mainly in the field range from about 10 G to 60 G. For stronger fields these lines are maximally depolarized. On the other hand, their Hanle effect depolarization is negligible for weaker fields, which makes them ideal reference lines below 10 G when compared to lines with much smaller N.
The ideal choice of spectral windows thus depends on the actual magnetic field strength. Based on previous results we selected narrow wavelength regions containing lines with N greater than 40 and others with N around 10. An additional selection criterion is a small number of atomic or molecular blends since they would significantly alter the intensity and polarization spectra.
Figure 3: -contours of the fitting procedure in the first wavelength region ( 3865.9-3866.8 Å). The panels show 2-dimensional cuts in the parameter space through the minimum . Numbers on the contours indicate the values of as defined in Eq. (16). | |
Open with DEXTER |
We compare our model calculations with the data from the second solar spectrum atlas compiled by Gandorfer (2005). All spectra were recorded with the UV version of the Zurich Imaging Polarimeter ZIMPOL II (Gandorfer et al. 2004) at the 1.5 m McMath-Pierce facility of the National Solar Observatory on Kitt Peak (Arizona, USA). The spectrograph slit was set parallel to the solar limb at the distance corresponding to the cosine of the heliocentric angle . Due to the telescope guiding and seeing, the value of was not completely stable, resulting in a small uncertainty of the exact limb distance. Because the scattering polarization is very sensitive to small variations in very close to the solar limb, behaves like a semi-free parameter for our analysis within strict limits defined by the error bars.
One of the main problems for our fitting procedure is the uncertainty of the zero point of the polarization scale in the observations. To solve this problem, we selected several wavelengths as continuum reference points, which are far away from very strong lines and least affected by blends. Based on the theoretical continuum polarization (Stenflo 2005), we shifted the polarization scale of the observations appropriately. The same wavelength points were used to renormalize the data.
We have chosen three wavelength regions for our calculations. These regions were not observed simultaneously and correspond, in principle, to different areas of the solar surface. Therefore measured parameters can vary from region to region.
The line parameters for all three regions are given in Tables 2-4 (available only in electronic form), which provide also weaker CN lines considered in the fit.
The first wavelength region 3865.9-3866.8 Å contains four CN violet system doublets with N=6,7 and 52 from the (1,1) vibrational band and with N=69 from the (0,0) vibrational band. Because of the relatively large Landé factors , the lines with N=6 and 7 are very sensitive to the magnetic field and reach the saturated Hanle effect already at 5-10 G, while the lines with N=52 and 69 are much less affected by the magnetic field (see Fig. 2). Therefore, this region is ideal for applying the differential Hanle effect technique, in particular since it contains lines with small and large N quantum numbers from the same vibrational band, as the relative strength of the lines from different bands can be influenced by other factors. Furthermore, all lines can be measured simultaneously in a single observation because they are spread over less than 1 Å.
Figure 3 shows the accuracy of the fitted spectrum illustrated by -contours within the parameter space. All four free parameters are well constrained by the data, and their best values are listed in Table 1. In particular we find a turbulent magnetic field strength of G. The magnetic field B and the scattering probability exhibit a somewhat linear dependence (Fig. 3, leftmost panel), which is easily understood. A greater scattering probability leads to additional line polarization that has to be compensated with a stronger Hanle effect and thus a larger magnetic field. Nonetheless, both B and remain restricted to well defined limits because of the differential sensitivity to the Hanle effect of the involved CN lines and the roughly line independent scattering probability. Similarly, the parameters and also show a tendency towards a linear dependence, because a greater effective temperature reduces the line depth in the intensity spectrum that can be balanced by additional line opacity.
Figure 4: Fits to the observations (thin solid lines) of Stokes and in the first wavelength region ( Å). The assumed magnetic field strengths in the calculated spectra are 0 G (no depolarization, thick dashed line), 27 G (best fit, thick solid line), and 45 G (too much depolarization, dotted line). The labels at the bottom of the lower panel indicate the position of the lines in the CN doublets with the corresponding N quantum number of the lower states. The labels at the top of the lower panel indicate the wavelength points which were chosen for the fitting procedure. | |
Open with DEXTER |
Figure 5: -contours as in Fig. 4 but for the second wavelength region close to the (1,1) bandhead. We obtain a minimum value . | |
Open with DEXTER |
Figure 6: Fits to the observations (thin solid lines) of Stokes and Q/I_{c} in the second wavelength region close to the (1,1) bandhead. The assumed magnetic field strengths in the theoretical spectra are 0 G (no depolarization, dashed line), 25 G (best fit, thick solid line) and 45 G (too much depolarization, dotted line). Labels at the bottom of the lower panel indicate the N quantum numbers of the lower states and the line positions within each doublet (close to the bandhead N is omitted due to overlap). The labels at the top of the lower panel indicate the wavelength points which were selected for the fitting procedure. | |
Open with DEXTER |
The modeled spectra are compared with observations in Fig. 4. For the calculated curves we have assumed the best values of the parameters , , and as listed in Table 1 while varying the magnetic field strength. The differential Hanle effect is clearly visible and constrains the magnetic field strength to well defined limits. It has to be stronger than 15 G to saturate the highly sensitive N=7 doublet. The upper limit for the field strength is defined by the only modest observed depolarizations of the lines with N=52 and 69. The observed doubling of the N=52 doublet probably appears due to a strong Ti I blend neglected in our model. We obtain the best fit for the magnetic field strength of 27 G.
Molecular bandheads have a high potential for successful differential Hanle effect diagnostics. They contain a mixture of lines with different J-numbers and correspondingly different magnetic field sensitivities. On the other hand, all these lines belong to the same rotational branch and vibrational band so that many of the relevant line formation parameters nearly coincide. Here we present the results of our fitting procedure in the (1,1) bandhead of the CN violet system.
The fitting procedure is illustrated in Fig. 5 by the -contours. Due to numerous neglected blends from atoms and molecules other than CN the quality of the fit is worse than in the first wavelength region (Sect. 3.3). Furthermore, the noise level in this region is also a bit higher. Nonetheless, we determine a magnetic field strength of 255 G in good agreement with the first region.
Figure 7: -contours as in Fig. 4 but for the third wavelength region near the (0,0) bandhead. We obtain a minimum value . | |
Open with DEXTER |
Figure 8: Fits to the observations (thin solid lines) of the CN violet system in the third wavelength region close to the (0,0) bandhead. The magnetic field strengths in the modeled spectra were set to 0 G (no depolarization, dashed line), 11 G (best fit, thick solid line), and 25 G (too much depolarization, dotted line). Labels at the bottom of the lower panel indicate the CN doublets with the N quantum numbers of the lower states. The labels at the top of the lower panel indicate the wavelength points chosen for the fitting procedure. | |
Open with DEXTER |
The modeled spectra are compared to observations in Fig. 6. The bandhead of the (1,1) band is clearly visible near 3871.5 Å. The selected wavelength region also contains a few lines from the (0,0) band. Contrary to the first region, lines with different Hanle effect sensitivities can not be easily distinguished as all lines are blended with each other. Even when neglecting the magnetic field it is possible to fit the overall average degree of polarization. However, the exact shape of the spectrum can only be reproduced when accounting for a magnetic field strength within well defined limits.
In this wavelength region we had to adopt a cosine of the heliocentric angle , which still lies within the error limits of the observed data. In fact, we failed to obtain good fits in this region with , because the continuum polarization in the incident radiation field was too large.
The CN lines near the (0,0) bandhead are very strong and heavily blended. This results in a significant optical thickness larger than 1 at the (0,0) bandhead so that the single scattering approximation loses its validity. Therefore, we selected a region which is 4 Å away from the bandhead.
This third wavelength region contains a mixture of lines with and from the P branch. The main features relevant for diagnosing the magnetic field in this spectral window are groups of neighboring lines containing a doublet with and one with , which are particularly well seen in the intensity spectrum. The doublets within one group can however be distinguished due to a consistent wavelength shift and different polarization characteristics. The lines with have a larger effective polarizability (compared to W_{2} close to the asymptotic limit of 0.1 for lines with ), which in principle causes a larger line polarization in a non-magnetic medium. The Hanle effect introduces a competing contribution to the relative line polarization because the lines with small N are much more sensitive to the magnetic field. From Fig. 8 we see that the consistent modeling leads to smaller polarization and the expected stronger Hanle effect in the lines with small N. Therefore, we can indeed infer the magnetic field strength from the spectrum.
The results of the fitting procedure are presented in Figs. 7 and 8. Again all the free parameters are well constrained and we reach a reasonable quality of the fit, except in three narrow regions due to missing blend lines in our model at 3879.0 Å, 3879.85 Å, and 3880.2 Å, which is clearly seen in the intensity spectrum. The best fit is obtained for a magnetic field strength of G, which is significantly smaller than in the first two analyzed spectral windows.
Table 1: Least squares fits to the magnetic field strength B, the optical depth normalization coefficient , the scattering probability , the temperature of the isothermal CN layer and cosine of the heliocentic angle in three considered regions of spectrum
The above comparison of the modeled spectra to the observed data was performed independently in the three selected wavelength regions. The obtained sets of free parameters for each region are collected in Table 1. They agree quite well in spite of the simple modeling, but there exist some differences that should be addressed.
The magnetic field strength found in the third region lies significantly below the value determined for the first two spectral regions. This is however not necessarily a contradiction, because the analyzed data were not observed simultaneously and sampled in principle different spatial positions on the Sun.
Furthermore, the line strength and accordingly the optical thickness changes from the first region to the third one. This influences the formation heights of the CN lines, which leads to different physical conditions in the corresponding atmospheric layers. Lines in the third region are stronger than in the first two and therefore form higher in the atmosphere where the temperature is smaller. A comparison of the contribution functions calculated for the typical lines in each spectral window shows that the difference in temperatures is around 100 K, consistent with the difference of 160 K obtained with our least-squares method (see Table 1). Our results suggest that the magnetic field strength may decrease with height, but a detailed interpretation of these small differences lies outside the scope of the current model and will be clarified by future work.
In addition, the continuum normalization and polarization in the third region remains an issue due to broad wings of nearby strong atomic lines including the H8 line from the Balmer series. The intensity normalization has little effect on the inferred magnetic field strength because it affects both the theoretical and observed signal in the same way. However, the additional opacity from those nearby atomic lines would probably reduce the continuum polarization and slightly increase the magnetic field strength by a few Gauss.
Overall, the physical parameters determined with the fit lie in a reasonable range. The small scattering probability is consistent with expectations and with a relatively small influence of non-LTE effects. The total optical thickness of the CN layer is close to but still smaller than 1 (maximal optical thicknesses in the considered regions are 0.26, 0.5 and 0.32 for the first, second and third regions correspondingly) in accordance with our assumptions. Finally, the temperature of the CN layer around 4500-4600 K that is required for the best fits, corresponds well to the temperature in the upper solar photosphere where the CN lines are formed.
An interesting extension of our method is the interpretation of the center-to-limb variations of the continuum and line radiation field, as it was done for Sr I 4607 Å line by Stenflo et al. (1997) and for the and molecular lines by Faurobert & Arnaud (2003) . Such observation for the CN violet system are not available yet, but could be very useful for constraining the atmospheric models.
Acknowledgements
This work was supported by SNF grants 20002-103696 and PE002-104552. S.B. acknowledges the EURYI Award from the ESF.
Table 2: Line parameters used for the fitting procedure in the first spectral region (3865.9-3866.8 Å). The columns list the wavelength (), the line identification (ID), the polarizability (W_{2}), the lifetime of the upper level (), the turbulent magnetic field strength at which the line polarization is reduced to 30% of its non-magnetic value neglecting elastic collision (B_{30}), the vibrational number which coincides for the upper and lower states (v), and the oscillator strength (f_{vJS}).
Table 3: Line parameters used for the fitting procedure in the second spectral region (in the (1,1) bandhead), with the same parameters as in Table 2.
Table 4: Line parameters used for the fitting procedure in the third spectral region (close to the (0,0) bandhead), with the same parameters as in Table 2.