VI/98A Stark broadening of H lines (Stehle+, 1999) ================================================================================ Extensive tabulations of Stark broadened hydrogen line profiles Stehle C., Huucheon R. =1999A&AS..140...93S ================================================================================ ADC_Keywords: Atomic physics Keywords: atomic data - atomic processes - lines profiles Description: New tables of Stark broadened hydrogen Lyman, Balmer and Paschen lines are presented for a wide range of plasma conditions, allowing the complete spectrum of these series to be modelled up to the Inglis Teller limit. The formalism, based on the Model Microfield Method (MMM) for both the electronic and ionic broadenings , gives an accurate description of the line profile from the centre to the wings. The profiles are calculated for a pure hydrogen plasma. All the transitions from n_lower = 1 (Lyman), 2(Balmer), 3(Paschen) to n_upper= (2 to 30) Lyman, (3 to 30) Balmer, (4 to 30) Paschen are reported. The temperature range is 2500 to 1.259E+06 K (10 temperatures) The electronic density range is 10e10 to 3.162E+19 (20 densities) A summary of the parameter values is reported in file "params.dat". The line profiles are tabulated in directories: ly/lyNN for Lyman lines, with NN = 02 .. 30 ba/baNN for Balmer lines, with NN = 03 .. 30 pa/paNN for Paschen lines, with NN = 04 .. 30 Each TRANSITION directory (for example ly/ly02, for Lyman{alpha}) contains at maximum 20 files, denoted by profilX.dat, ie, profil1.dat, ..., profil20.dat, which correspond to different increasing densities. Due to the plasma lowering of the ionization limit (Inglis Teller limit), the number of reported density calculations decreases when NN increases. The number of density files profilX.dat is reported for each transition directory in the file nraie.dat. File Summary: -------------------------------------------------------------------------------- FileName Lrecl Records Explanations -------------------------------------------------------------------------------- ReadMe 80 . This file params.dat 85 7922 Summary of physical parameters read_gen.f 100 796 Interpolation Fortran90 program to compute the line profiles at a given temperature,density ly/* . 29 Lyman lines (lower number = 1) ba/* . 28 Balmer lines (lower number = 2) pa/* . 27 Paschen lines (lower number = 3) -------------------------------------------------------------------------------- Byte-by-byte Description of file: params.dat -------------------------------------------------------------------------------- Bytes Format Units Label Explanations -------------------------------------------------------------------------------- 1- 2 I2 --- lower [1/3] Lower state (1=Lyman, 2=Balmer, 3=Paschen) 4- 5 I2 --- upper [2/30] Upper state (upper > lower) 7- 14 F8.2 0.1nm lambda0 Wavelength of the line center 16- 24 E9.3 --- Kalpha Asymptotic value of iStark*(alpha^2.5^) 26- 34 E9.3 cm-3 Ne Electron density 36- 44 E9.3 30kV/m E0 *Dlambda/Dalpha = E0, normal Holtsmark field strength 46- 54 E9.3 K T Temperature 56- 64 E9.3 --- R0/Debye *Ratio of the mean interelectronic distance to the electronic Debye length 66- 85 A20 --- FileName Name of the file detailing the profile -------------------------------------------------------------------------------- Note on E0: the normal Holstmark electric field strength is: E0 = e * (4*pi/3) * Ne^(2/3)^ Its value is (in ues-cgs electric field units uesVolts/cm): E0 = 1.25E-9 * Ne^(2/3)^ in CGS-ues units (Ne expressed in cm^-3^) Note on R0/Debye: the Debye length is Debye = sqrt(eps0*k*T/Ne)/e (eps0 = dielectric constant) = 6.92 sqrt(T/Ne) cm (CGS units, Ne in cm^-3^) = 69.2 sqrt(T/Ne) m (SI units, Ne in m^-3^) The ratio R0/Debye is: 0.0898 (Ne^(1/6)^)/(T^(1/2)^) -------------------------------------------------------------------------------- Description of files: */*/profil*.dat -------------------------------------------------------------------------------- Each of these files contains a maximum of 10 profiles (with and without Doppler convolution) for 10 plasma temperatures. T1, T2, ... T10. The header lines are: line 1 : n_lower , n_upper, wavelength of the transition (Angstroms) line 2 : electronic density in cm^-3^ line 3 : D{lambda}/D{alpha} = E0, normal Holtsmark field strength (see the "Note on E0" section above) line 4 : Asymptotic value of iStark*({alpha}^2.5^) line 5 : blank line 6 : Temperature in K ( T1,T2,T3, then T4,T5,T6,then T7,T8,T9, then T10) line 7 : R0/Debye Ratio (see the "Note on R0/Debye" section above) line 8 : HHWM of the line (without Doppler broadening) line 9 : " " " " (with " ) -------------------------------------------------------------------------------- The 7 columns contain: column 1 : Dalpha = Detuning from central wavelength, expressed in units of {alpha} = lambda0 / E0 where lambda0 is the unperturbed wavelength (in {AA} = 0.1nm) and E0 is the normal Holstmark electric field in uesVolts/cm column 2 : iDoppler = Intensity with Doppler (convolved) for T1 in units of 1/{alpha} column 3 : iStark = Intensity with Stark only (in brackets) for T1 in units of 1/{alpha} column 4 : iDoppler (convolved) for T2 column 5 : iStark (in brackets) for T2 column 6 : iDoppler (convolved) for T3 column 7 : iStark (in brackets) for T3 -------------------------------------------------------------------------------- Description of program: read_gen.f This program should be run in the directory related to a specific transition (e.g. in directory ly/ly02 for Lyman{alpha}). The input to this program consists in: 1. Temperature (K) Example: 10000 2. Electronic density (cm^-3^) Example: 1.e14 The output of this program consists in four columns with: Column 1: wavelength distance from line centre ({AA} = 0.1nm) Column 2: frequency distance from line centre (in Hz) Column 3: iDoppler (intensity convolved with Doppler) in s (seconds) Column 3: iStarck (intensity without convolution) in s (seconds) Acknowledgements: Chantal Stehle History: This catalog supersedes the previous version (VI/82) * 18-Sep-2000: a corrected version of the following files was introduced: ly/ly02/profil17.dat ly/ly03/profil17.dat ba/ba04/profil16.dat read_gen.f (one line changed) * 28-Aug-2001: a few corrections in the reading program (read_gen.f) * 23-Nov-2007: paschen parameters (files pa/*) corrected (from Chantal Stehle). Original files kept as pa.ori * 14-Jan-2008: improved version of the interpolation program. References: Brissaud A., Frisch U., 1971, J. Quant. Spectrosc. Radiat. Transfer 11, 1767 Frisch U., Brissaud A., 1971, J. Quant. Spectrosc. Radiat. Transfer 11, 1753 Vidal C.R., Cooper J., Smith E.W., 1973ApJS...25...37V ================================================================================ (End) Francois Ochsenbein, Patricia Bauer [CDS] 14-Jun-1999