Hindawi Publishing Corporation International Journal of Spectroscopy Volume 2010, Article ID 506346, 6 pages doi:10.1155/2010/506346

Research Article

Hydrogen Stark Broadened Brackett lines

C. Stehle and S. Fouquet

LERMA, UMR8112, Observatoire de Paris, CNRS et Universite Pierre et Marie Curie, 5 place Jules Janssen, 92195 Meudon, France Correspondence should be addressed to C. Stehlé, chantal.stehle@obspm.fr Received 30 April 2009; Accepted 2 July 2009 Academic Editor: Manuel A. González

Copyright © 2010 C. Stehle and S. Fouquet. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Stark-broadened lines of the hydrogen Brackett series are computed for the conditions of stellar atmospheres and circumstellar envelopes. The computation is performed within the Model Microfield Method, which includes the ion dynamic effects and makes the bridge between the impact limit at low density and the static limit at high density and in the line wings. The computation gives the area normalized line shape, from the line core up to the static line wings.

1. Introduction

Hydrogen is the most abundant element in the universe. Its broad lines give noticeable features in the spectra of stellar atmospheres [1-4]. These lines are very sensitive to the interaction between hydrogen radiating atoms and the surrounding charges, electrons, and ions (mostly protons), which is connected to the random electric field generated by these charges. The electric field has two components, with different time scales: the rapidly varying electronic field and the slowly varying ionic electric field. The net field induces a strong mixing of the atomic states with the same principal quantum number, from which the Stark broadening originates. Astrophysical applications need to know the full line shape, from line centre up to the line wings, for a wide range of plasma conditions (i.e., electron density Ne and temperature T).

Whereas hydrogen line shapes have been widely used in stellar physics for the determination of gravity and/or temperature in the visible and UV ranges, the recent instrumental developments on various telescopes, such as AMBER [5] and CRIRES [6] on VLTI and VLT, require the availability of precise hydrogen lines in the infrared (i.e., between 1 and 5 ^m). The latter are badly known [7].

This paper thus aims at providing a coherent description of the line shapes of Brackett a, ft, y, which connect the levels of principal quantum numbers n equal to 4, for the lower state of the transition, and n' equal to 5, 6, and 7, for

the upper state, and have central wavelengths of 4.05, 2.63, and 2.12 ^m. The chosen plasma conditions are relevant to stellar photospheres and circumstellar environments (electron densities between 1010 and 1019cm-3 and temperatures between 103 and 107K). The plasma charges are assumed to be electrons and protons. This is a standard approximation for this type of study, although some improvements may be possible by including the effect of ionization of He (which is, by number of atoms, only 10% as abundant as hydrogen) and traces of heavier elements.

We suppose that plasma collective effects are included in the Coulomb interactions between the hydrogen bound electron and the plasma charges by Debye screening. This requires that the Debye length Xd is larger than the mean distance r0 between the protons or that the parameter a = r0/Xd is smaller than unity. Our tables will thus be limited to the values of temperatures and densities satisfying the following condition:

Ne< 2.6 109 T3 (cm-3, K), (1)

which is fulfilled for standard stellar atmospheres.

We assume also that the proton/electron density is small enough to ensure that each line 4-n'(n' = 5,6 or 7) remains distinguishable from the subsequent line 4 - (n' + 1) of the Brackett series. Using the Inglis-Teller [8] criterium, the upper limit to the electron density Ne is thus

log(Ne,max) = 22 - 7.5log(n') (cm-3). (2)

This, respectively, gives log(Ne,max) = 16.75,16.16, and 15.66 for the n' (= 5, 6 and 7) values considered in this paper. However, in order to allow interpolation within the tables for astrophysical applications, the line shapes for higher values of Ne, up to a decade, have been computed.

Different methods can be used to generate high-quality spectral line shapes for hydrogen lines perturbed by protons and electrons: Molecular Dynamics for describing precisely the ion dynamics effects in the line center [9, 10], quantum theory for the electron contribution to the line wings [11], and short range H-H+ molecular interactions also for the line wings, leading to the apparition of quasimolecular satellites, which are observed in the atmospheres of white dwarfs [12]. However, they are limited either to part of the profile, to restricted plasmas conditions, or to simple lines like Lyman or Balmer lines. Thus, as for astrophysical purposes, the tabulations go from the line centre up to the line wings, and it is necessary to find a compromise between accuracy and description of the whole profile.

In this context, the tabulations of Vidal et al. [13], using Unified Theory for the electrons and static approximation for the ions, have been used for stellar atmospheres, despite the intrinsic lack of accuracy in the line center due to the neglect of ion dynamics effects. The tabulations of Stehle et al., for the Lyman, Balmer [14, 15], and Paschen lines [16], using Model Microfield Method, which was initially developed by Brissaud et Frisch [17, 18] brought an important improvement by taking into account the ion dynamics effects. They are now used for atmospheres and for the computation of radiative diffusion processes in the radiative stellar envelopes [19]. In the case of partly ionized plasmas, like for the atmospheres of cool stars, the contribution of neutral broadening by hydrogen has to be included in the line shape, especially in the line wings [20, 21]. This effect will be neglected in the following.

In this paper we will present the Stark broadened profiles of Brackett lines. They will be computed within the formalism of Model Microfield Method, hereafter denoted by MMM. We will neglect the fine structure effects, which play a role, at low density, in the core of the lines with low n quantum mumbers, like Lya or Ha [22]).

For Bra, the profiles are computed at densities log10(Ne(cm-3)) ranging from 10 to 18.5, by step of 0.5. For each density, the profiles are computed at temperatures equal to 1000, 2500, 5000, 10000, 19550, 39810, 79810, 158500, 316200, 63100, and 1259000 K, assuming that condition 1 is satisfied. For Br^ , the computation stops at log10(Ne(cm-3)) = 18 and for Bry at 17.5, as explained previously.

2. Method

The broadening of spectral lines results from the interactions between the radiating hydrogen atom and the free ions and protons. These two contributions can be described in terms of interaction potentials, with the corresponding electronic and ionic plasma microfields Fel and Fion. Neglecting quadrupolar and other contributions that play a role at high

densities [23], the dipolar potential of interaction between the bound electron and the microfields may be written as

F(t) = -d ■ (Fel(t)+ Fion(t)).

The spectral line profile I(w) (with area normalized to unity) is thus defined in the Liouville space [24, 25], spanned by the states | i, f » (which stands for | nj, lj, mf, nf, lf, mf ») as

I (w) =

v j ji~Re X di,f ■ d* f' nXi,f di,f ■ d*f ^f >> i f (4)

x <U(w)>

el,ion; if;i' f ''

where (U(w))dion is the Fourier transform of the evolution operator of radiating Hydrogen atom in the Liouville space, averaged over the realizations of the stochastic dynamic electronic and ionic microfields Fel and Fion. The term d;,f ■ dj',f is the product of dipole operator elements between initial low states (denoted by i, i') and final upper states (denoted by f, f'), of the hydrogen bound electron. As fine structure and inelastic effects are neglected, one has Ef - Ej = E'f - E' = hw0.

The two microfields are stochastic processes. It is thus possible to define two distribution functions P(F) [26, 27], respectively, associated to the slowly varying ionic and rapidly varying electronic microfields. In order to take into account the dynamic effects of these microfields, a model for the dynamical statistics of field fluctuations is necessary. In MMM, the microfield fluctuations are handled with a statistical process model, where the microfield (electronic or ionic) is assumed to be constant during a given time interval. The microfield then jumps instantaneously to another constant value for the next time interval. The jumping times are assumed to follow a Poisson law, with a field dependent frequency v(F). The jumping frequency v(F) is chosen to reproduce the true field autocorrelation function [17, 18, 28]. This method has been tested against asymptotic impact and quasistatic limits and has been proved to lead to very good results for hydrogen [29] and hydrogenic ion lines [28].

The method has been already described in Stehle and Hutcheon [15], and we refer the reader to this paper for the details. An important point is that it is possible to disentangle the contributions of ions and fast electrons by introducing a frequency dependent electronic relaxation operator yel(w), which is independent from ionic fields and thus may be computed separately. The Fourier transform of the evolution operator, (U(w))el, averaged over all the realisations of the electronic fields, may be written as

<U(w)>el = i [AwI + !>el(w)]"

where Aw = w - is the detuning from line center, and I is the identity operator in the Liouville space.

Thus, this electron damping is first computed to account for average effect of the electronic fields. Then, the static Fourier transform of the evolution operator, (U(w))el,ion,

B m £y

1012 1013

Ne (cm-3)

Figure 1: HWHM of Bra at 2500 K, in A, versus electron density Ne in cm-3 (black and squares: MMM Stark only; red and circles: analytical impact formula; dashed green: MMM Stark profile with Doppler convolution).

averaged over the realisations of the electronic and ionic fields, may be written as

№)>e,;0n = i P(Fion )dFion

The MMM expression is more complex [17, 18] than this one, which corresponds to the usual Unified Theory with static ionic fields.

In order to reduce the dimensions of the Liouville space, (16 X 25 = 400 states for Bra, i.e.), we use the formalism of the reduced Liouville space, which takes advantage of the invariance of the different operators, like d ■ d , by angular average over all the orientations of the electric fields and of the fact that the dipole tensor d in (4) is of rank 1 (see e.g. [28, 29]). However, the number of reduced states \ ni, lf, nf, lf » (with \li - lf \ < 1) remains important (i.e., 10 reduced Liouville states of rank 1 for Bra). Thus, we use another approximation, already called "isotropic" approximation in [15], which uses the diagonal form yes°(w), with equal diagonal matrix elements, instead of yel(w). This scalar tensor is deduced from the pure electronic profile I( w)el by the following relation:

I (w)el = -1 Im[Aa + /yéWr1,

Si,fdi,f ■ tfj, Yel; if;i' f (w)

Si,f di,f ■ dtf

The Stark profile is thus obtained after computing the average over electronic field values, which gives yes°(w),

1012 1013

Ne (cm-3)

Figure 2: Same as Figure 1 for Brß.

and then the average over the ionic fields, both using the formalism of Model Microfield Method. The third step is then the Doppler convolution, which gives the final line shape.

3. Results

The three Brackett a, ft, y line shapes have been studied in a wide range of stellar conditions, but we report hereafter only a selection of results, relative to the line center and the line wings.

3.1. Line Center. We present in Figures 1, 2, and 3 the variations of the half-width (HWHM) of Bra, Brft, and Bry lines for different values of the electron density. This quantity is the detuning from the line center, at the point where the profile reaches half the maximum value of the line profile (which is not necessarily at the line center, as will be discussed below). The figures show the half width of the MMM profile with and without Doppler effect and also the value of the impact half-width in its own validity range.

Hence, impact limit has been proved to be valid, for both electron and ion contributions, in the line center and at low densities. Moreover, it has been proved that the value of the impact width is analytical and that the corresponding profile is Lorentzian in the line center [30]. The validity condition of the impact limit is that the half-width value is smaller than the ion plasma frequency. Let us take the example of Bra. At 1012 cm-3 and 2500 K the ion plasma frequency is equal to 2.6 109rd ■ s-1, whereas the ion impact contribution to the HWHM is equal to 1.8 109 rd ■ s-1. At 1014cm-3, they are respectively, equal to 2.6 1010 and 3.9 1010. The impact limit should thus be reached gradually as the density decreases below 1013cm-3.

Figure 1, relative to Bra, shows indeed that the half-width of the Stark profile (black) tends to converge towards the impact analytical limit (red) at these low densities.

Ne (cm-3) Figure 3: Same as Figure 1 for Bry.

However, the convolution with the Doppler profile increases the half-width value. As a consequence, the half width is dominated by the Doppler broadening at low densities. Similar behaviour occurs for the other lines, as may be on Figures 2 and 3 for Br3 and Bry.

At higher densities, the ion dynamic effects become smaller in the line center (they are negligible in the line wings, as will be seen below), and the line shape departs from the Lorentzian shape. We found that, depending on the temperature conditions, the Brackett 3 line may present a small dip in the line center at moderate densities of 1015 cm-3 and for the largest temperatures as seen in Figure 4. This effect was well known in earlier tabulations for the Lyman^ and H3 lines. The dependence versus the temperature is a consequence of the the electron impact broadening, which varies in T-1/2 and which may fill (at low T values) or not (at large T values) the central dip.

3.2. Line Shapes. Another typical behaviour of hydrogen lines is the convergence towards the (Holtsmark) static limit in the line wings, which scales as |Aw|-5'2 when the line shape is expressed in angular frequency units (which is the appropriate unit for the line shape computations). However, traditionally, the line intensity was expressed in units of Aa = AX/Fo , where Fo (esu) = 1.25 109(Ne)2/3 is the normal Holts-mark field. In these units, the Holtsmark limit is given by

„A ^ Ka ( Xq/Fq + Aa \ 1 (Aa) =

for |Aa|

where X0 is the central wavelength, and Ka a constant, which depends on the transition,

Ka = 1'512 for Bra, 2'401 for Br^,2'926 for Bry' (10)

< 0.06

i Brß

j.ea-a—1

\ \-cj X

0 1 2 3 4 5 6

Figure 4: MMM Br3 line at 1015 cm-3 versus Aa = AX/F0 without Doppler convolution for different temperatures (red : 1000 K; blue with empty circles : 2500 K; violet with empty squares: 10000 K; green: 19950 K).

' \ 1 Bra

i i i

Figure 5:MMMBra lineat Ne = 1012cm-3, with (red) and without (blue) Doppler broadening, and the asymptotic Holstmark limit (black) in |Aa|-5/2 for 11 different temperatures in K (1000, 2500, 5000, 104, 2 104, 3.98 104, 7.94 104, 1.58 105, 3.16 105, 6.31 105, 1.26 106). The corresponding values of the profiles in the line center are, respectively, 0.28, 0.30, 0.34, 0.46, 0.56, 0.69, 0.88, 1.23, 1.45, 1.87, without Doppler effect (blue), and, 0.102,0.071,0.053, 0.039,0.028, 0.020, 0.014, 0.010, 0.0072, 0.0051, 0.0036 including Doppler effect (red).

This variation in Aa of 9 introduces, at large detunings, a "trivial" asymmetry between I (Aa) and I (-Aa) (which does not exist between I (Aw) and I (-Aw)). Figures 5, 6 and 7 show the profiles I (Aa) of Bra, Br3, and Bry lines at 1012 cm-3 and various temperatures. The profiles, including Doppler effects are reported in red color, and the pure Stark profiles in blue.

-I io-6

1 \i Brß

Figure 6: Same as Figure 5 but for Brft, the values of the profiles in the line center are, for the 11 different temperatures of Figure 5, respectively, 0.124, 0.129, 0.136, 0.149, 0.168, 0.195, 0.233, 0.285, 0.355, 0.451, 0.579 without Doppler effect (blue), and 0.099, 0.082, 0.067, 0.053, 0.040, 0.030, 0.0215, 0.015, 0.011, 0.0078, 0.0056 including Doppler effect (red).

10-4 -

< 10-6

10-8 -

10-10 -

Figure 7: Same as Figure 5 but for Bry, the values of the profiles in the line center are, f for the 11 different temperatures of Figure 5, respectively, 0.107, 0.107, 0.109, 0.113, 0.121, 0.132, 0.150, 0.175, 0.211, 0.260, 0.327 without Doppler effect (blue), and 0.091, 0.080, 0.069, 0.057, 0.0450, 0.034, 0.025, 0.018, 0.013, 0.0094, 0.0067 including Doppler effect (red).

The profiles are area normalized, that is,

I (Aa) X d(Aa) = 1.

As a consequence, broad profiles have small values of I(Aa = 0). For the density considered in Figures 5-7, the Doppler profile dominates the Stark profile in the line

center and becomes indistinguishable, as expected, in the line wings, where they follow the asymptotic limit of (9). At higher densities (not reported here), the Doppler width is smaller than the Stark width, and Doppler convolution is no longer necessary.

4. Conclusions

This study shows that Stark broadened infrared Brackett lines of hydrogen follow the same trends as the lines of lower series which are more known theoretically and experimentally. This study will allow missing absorption in the spectra of stellar atmospheres in the infrared due to the lack of data to be filled in a next future. Dedicated experiments and comparisons with other methods, for instance, FFM [31], would be helpful to test these theoretical results, which will be also constrained by observational data. The corresponding tables will be accessible at http://amrel.obspm.fr/stark-h.

Acknowledgments

The authors would like to thank Nicolas Moreau for his helpful technical support, and Pauline Beghin, Edrice Bouteldja, Aurelie de Paz and Cecile Turc, students at University of Pierre et Marie Curie, Paris, for their contribution in checking the data. The support of French INSU program PNPS is acknowledged.

References

[1] P. S. Barklem, "Hydrogen lines," Physica Scripta, vol. 133, Article ID 014023, 5 pages, 2008.

[2] I. Hubeny and T. Lanz, "Non-LTE line-blanketed model atmospheres of hot stars. I. Hybrid complete linearization/accelerated lambda iteration method," Astrophysical Journal, vol. 439, no. 2, pp. 875-904, 1995.

[3] N. F. Allard and D. Koester, "Theoretical profiles of Lyman-alpha satellites and application to synthetic spectra of DA white dwarfs," Astronomy & Astrophysics, vol. 258, no. 2, pp. 464-468, 1992.

[4] M. F. Nieva and N. Przybilla, "Hydrogen and helium line formation in OB dwarfs and giants: a hybrid non-LTE approach," Astronomy & Astrophysics, vol. 467, no. 1, pp. 295309, 2007.

[5] R. G. Petrov, F. Malbet, G. Weigelt, et al., "AMBER, the near-infrared spectro-interferometric three-telescope VLTI instrument," Astronomy & Astrophysics, vol. 464, no. 1, pp. 112, 2007.

[6] H. U. Kaufl, P. Amico, P. Ballester, et al., "CRIRES: commissioning and first science results," in Ground-based and Airborne Instrumentation for Astronomy II, vol. 7014 of Proceedings ofSPIE, Marseille, France, June 2008.

[7] N. Przybilla and K. Butler, "Non-LTE line formation for hydrogen revisited," Astrophysical Journal, vol. 609, no. 2, pp. 1181-1191,2004.

[8] D. R. Teller and E. Inglis, "Ionic depression of series limits in one-electron spectra," Astrophysical Journal, vol. 90, no. 3, pp. 439-448, 1968.

[9] M. A. Gigosos and V. Cardenoso, "New plasma diagnosis tables of hydrogen Stark broadening including ion dynamics," Journal of Physics B, vol. 29, no. 20, pp. 4795-4838, 1996.

R. Stamm, E. W. Smith, and B. Talin, "Study of hydrogen Stark profiles by means of computer simulation," Physical Review A, vol. 30, no. 4, pp. 2039-2046, 1984.

C. Stehle, "Stark broadening of the hydrogen Lyman-a line from the center to the near line wings for low-density plasmas," Physical Review A, vol. 34, no. 5, pp. 4153-4162, 1986.

N. F. Allard, D. Koester, N. Feautrier, and A. Spielfiedel, "Freefree quasi-molecular absorption and satellites in Lyman-alpha due to collisions with H and H+," Astronomy & Astrophysics, vol. 108, pp. 417-431, 1994.

C. R. Vidal, J. Cooper, and E. W. Smith, "Hydrogen Stark-broadening tables," The Astrophysical Journal Supplement, vol. 25, no. 214, pp. 37-136, 1973.

C. Stehle, "Stark broadening of hydrogen Lyman and Balmer in the conditions of stellar envelopes," Astronomy & Astrophysics, vol. 104, pp. 509-527, 1994.

C. Stehle and R. Hutcheon, "Extensive tabulations of Stark broadened hydrogen line profiles," Astronomy & Astrophysics, vol. 140, pp. 93-97, 1999.

C. Stehle, "Paschen lines of hydrogen and He+ ion," Physica Scripta, no. T65, pp. 183-187, 1996.

U. Frisch and A. Brissaud, "Theory of Stark broadening—

I soluble scalar model as a test," Journal of Quantitative Spectroscopy & Radiative Transfer, vol. 11, no. 12, pp. 17531766, 1971.

A. Brissaud and U. Frisch, "Theory of Stark broadening—

II exact line profile with model microfield," Journal of Quantitative Spectroscopy & Radiative Transfer, vol. 11, no. 12, pp. 1767-1783, 1971.

F. Clausset, C. Stehle, and M. C. Artru, "Parametrized Stark profiles of hydrogen lines in hot stellar plasmas," Astronomy & Astrophysics, vol. 287, no. 2, pp. 666-675, 1994. P. S. Barklem, H. C. Stempels, C. Allende Prieto, O. P. Kochukhov, N. Piskunov, and B. J. O'Mara, "Detailed analysis of Balmer lines in cool dwarf stars," Astronomy & Astrophysics, vol. 385, pp. 951-967, 2002.

P. S. Barklem, N. Piskunov, and B. J. O'Mara, "Self-broadening in Balmer line wing formation in stellar atmospheres," Astronomy & Astrophysics, vol. 363, pp. 1091-1105, 2000. C. Stehle and N. Feautrier, "Stark broadening of the hydrogen Ha line at low densities: fine structure and spontaneous emission effects," Journal of Physics B, vol. 18, no. 7, pp. 12971306, 1985.

C. Stehle, D. Gilles, and A. V. Demura, "Asymmetry of Stark profiles," The European Physical Journal D, vol. 12, no. 2, pp. 355-367, 2000.

U. Fano, "Description of states in quantum mechanics by density operator techniques," Reviews of Modern Physics, vol. 29, no. 1,pp. 74-93, 1957.

A. Ben-Reuven, "Symmetry considerations in pressure-broadening theory," Physical Review, vol. 141, no. 1, pp. 34-40, 1966.

C. F. Hooper Jr., "Electric microfield distributions in plasmas," Physical Review, vol. 149, pp. 77-91, 1966. C. F. Hooper Jr., "Low-frequency component electric microfield distributions in plasmas," Physical Review, vol. 165, pp. 215-222, 1968.

C. Stehle, "Stark profiles of He+," Astronomy & Astrophysics, vol. 292, no. 2, pp. 699-712, 1994.

C. Stehle, A. Mazure, G. Nollez, and N. Feautrier, "Stark broadening of hydrogen lines—new results for the Balmer lines and astrophysical consequences," Astronomy & Astrophysics, vol. 127, no. 2, pp. 263-266, 1983.

[30] C. Stehle, "Analytical expressions for the lines of one-electron ions," Astronomy & Astrophysics, vol. 305, no. 2, pp. 677-686, 1996.

[31] B. Talin, A. Calisti, L. Godbert, R. Stamm, R. W. Lee, and L. Klein, "Frequency-fluctuation model for line-shape calculations in plasma spectroscopy," Physical Review A, vol. 51, no. 3, pp. 1918-1928, 1995.

Copyright of International Journal of Spectroscopy is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.