Next Article in Journal
Analytical Solution of the Hanle Effect in View of CLASP and Future Polarimetric Solar Studies
Previous Article in Journal
Critical Stability of the Negatively Charged Positronium-Like Ions with Yukawa Potentials and Varying Z
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transition Rates for 3s3p2 4P–3s3p4s 4Po Transitions in Al i

by
Charlotte Froese Fischer
1,* and
James F. Babb
2,*
1
Department of Computer Science, University of British Colombia, 2366 Main Mall, Vancouver, BC V6T1Z4, Canada
2
Institute for Theoretical Atomic, Molecular, and Optical Physics (ITAMP), Center for Astrophysics | Harvard & Smithsonian, 60 Garden St., MS 14, Cambridge, MA 02138, USA
*
Authors to whom correspondence should be addressed.
Submission received: 1 May 2019 / Revised: 23 May 2019 / Accepted: 26 May 2019 / Published: 4 June 2019
(This article belongs to the Section Atomic, Molecular and Nuclear Spectroscopy and Collisions)

Abstract

:
Fully relativistic calculations have been performed for two multiplets, 3 s 3 p 2 4 P and 3 s 3 p 4 s 4 P o , in Al i. Wave functions were obtained for all levels of these multiplets using the grasp programs. Reported are the E1 transitions rates for all transitions between levels of these multiplets. Transition energies and transition rates are compared with observed values and other theory. Our calculated transition rates are smaller by about 10% than observed rates, reducing a large discrepancy between earlier calculations and experiments.

1. Introduction

Atomic spectra are of vital importance as plasma diagnostics and reliable wavelengths and transition probabilities are essential for applications. Recently, Hermann et al. [1] deduced transition probabilities between the fine structure lines connecting the 3 s 3 p 2 4 P and 3 s 3 p 2 4 P o multiplets from emission coefficients measured in laser ablation of aluminum in argon. The transition probabilities obtained were roughly a factor of two larger than those listed in the semi-empirical calculations of Kurucz and Peytremann [2], and no other values were found in existing tabulations. The fine structure lines connecting the 3 s 3 p 2 4 P and 3 s 3 p 2 4 P o multiplets occur in the same general wavelength region (305–310 nm) and with comparable strength as the well-studied (see [3] and the references therein) fine structure lines connecting the 3 s 2 3 p 2 P and 3 s 2 3 d 2 D multiplets, making the discrepancy of concern for applications in the UV wavelengths.
Multiconfiguration Dirac–Hartree–Fock (MCDHF) and relativistic configuration interaction (RCI) calculations have been performed by Papoulia et al. [3] for 28 states in neutral Al. The configurations of interest were 3 s 2 n l for n = 3 , 4 , 5 with l = 0–4, as well as 3 s 3 p 2 and 3 s 2 6 l for l = 0 , 1 , 2 . Lifetimes and transition data for radiative electric dipole (E1) transitions were reported. There was a significant improvement in accuracy, in particular for the more complex system of neutral Al i, which may prove useful for astrophysical applications to Al abundance determinations in stars. Omitted were the levels of the 3 s 3 p 4 s 4 P o multiplet, which lies above the first 3 s 2 ionization limit [4].
This paper reports transition rates for all E1 transitions between the 3 s 3 p 2 4 P and 3 s 3 p 4 s 4 P o multiplets using the variational multiconfiguration Dirac–Hartree–Fock (MCDHF) method [5], as implemented in the grasp programs [6]. The accuracy of the results is based on the accuracy of the theoretically-predicted transition energies compared with available measurements, as well as the agreement between length and velocity rates.

2. Underlying Theory

In the MCDHF method [5], the wave function Ψ ( γ P J M J ) for a state labeled γ P J M J , where J and M J are the angular quantum numbers and P the parity, is expanded in antisymmetrized and coupled configuration state functions (CSFs):
Ψ ( γ P J M J ) = j = 1 M c j Φ ( γ j P J M J ) .
The labels { γ j } denote other appropriate information about the CSFs, such as orbital occupancy and the coupling scheme. The CSFs are built from products of one-electron orbitals, having the general form:
ψ n κ , m ( r ) = 1 r P n κ ( r ) χ κ , m ( θ , φ ) ı Q n κ ( r ) χ κ , m ( θ , φ ) ,
where χ ± κ , m ( θ , φ ) are two-component spin-angular functions. The expansion coefficients and the radial functions are determined iteratively. In the present work, the Dirac–Coulomb Hamiltonian H D C was used [5], which included a correction for the finite size of the nucleus.
The radial functions { P n κ ( r ) , Q n κ ( r ) } were determined numerically as solutions of differential equations,
w a V ( a ; r ) c d d r κ a r c d d r + κ a r V ( a ; r ) 2 c 2 P a ( r ) Q a ( r ) = b ϵ a b δ κ a κ b P b ( r ) Q b ( r ) ,
where V ( a ; r ) = V n u c ( r ) + Y ( a ; r ) + X ¯ ( a ; r ) is a potential consisting of nuclear, direct, and exchange contributions arising from both diagonal and off-diagonal Φ α | H D C | Φ β matrix elements [5].
For a given set of radial functions, expansion coefficients c = ( c 1 , , c M ) t were obtained as solutions to the configuration interaction (CI) problem,
H c = E c ,
where H is the CI matrix of dimension M × M with elements
H i j = Φ ( γ i P J M J ) | H | Φ ( γ j P J M J ) .
Once self-consistent solutions have been obtained—sometimes referred to as the relativistic MCDHF or RMCDHF phase—an RCI calculation was performed using an extended Hamiltonian that included the transverse photon (Breit) and QED corrections. Wave functions from the latter Hamiltonian were used to compute the E1 transitions rates.

3. Systematic Procedures

Systematic procedures were used in which the orbital set used for defining the wave function expansion increased systematically within a correlation model.
The states of Al consist of a neon-like ( 1 s 2 2 s 2 2 p 6 ) core and three valence electrons. Wave function expansions were obtained from single- and double- (SD) excitations from a multireference (MR) set that interacted significantly with the CSFs of interest. For the even multiplet, the MR set included the CSFs from 3 s 3 p 2 , 3 s 3 d 2 , and 3 p 2 3 d configurations and for the odd multiplet, CSFs from 3 s 3 p 4 s , 3 p 3 d 4 s , 3 p 2 4 p , and 3 s 3 d 4 p . Because the odd multiplet was above the first ionization limit, all 3 s 2 n l CSFs were removed. In this paper, the orbital sets that define the set of excitations were classified according to the largest n l of the orbital set when the latter were ordered globally by n (the principal quantum number) and within n by l (the orbital quantum number). Thus, an n = 3 orbital set includes all orbitals up to 3 s , 3 p , 3 d and all n = 5 f orbitals up to 5 s , 5 p , 5 d , 5 f ( 5 g not included).
Our first model was the valence correlation (VV) model, in which all excitations involved only valence electrons. n = 3 calculations were performed for an average energy functional of the lowest even parity J = 1 / 2 , 3 / 2 , 5 / 2 states. This calculation defined the core orbitals for all subsequent calculations. The n = 4 , 5 f , 6 f calculations each varied only the new orbitals. For the odd multiplet, the first calculation had orbitals up to 4 s , 4 p , and 3 d , but is still referred to as an n = 3 calculation in this paper, with remaining sets being the regular n = 4 , 5 f , 6 f orbital sets. Table 1 shows the convergence of the fine-structure of the two multiplets and their separation.
The results of the converged valence correlation calculations, when compared with observation, suggest that the energy structure is not significantly affected by the core-valence (CV) that accounts for the polarization for the core. To confirm this conclusion, calculations were performed in which SD excitations included a single excitation from the 2 p -shell along with a single excitation of a valence electron. Wave function expansions were considerably larger and convergence a bit slower. Table 2 shows the convergence of the energy structure. The fine structure of the odd multiplet increased slightly and was in somewhat better agreement with the data from observation. At the same time, the transition energy for 4 P 1 / 2 4 P 1 / 2 o for an n = 7 f calculation was not in as good agreement with observed as the n = 6 f valence correlation calculation reported in Table 1.

4. Results

The wave functions from RCI expansions, determined using the Dirac–Coulomb–Breit-QED Hamiltonian, were used to compute the E1 transition rates for all transitions between these two multiplets. Table 3 reports the transition energy Δ E (cm 1 ), the wavelength λ (nm) in a vacuum, A ( μ s 1 ), and g f , in the length form for calculations of Table 1. Furthermore included is an indicator of accuracy d T = ( A l A v ) / max ( A l , A v ) , where A l and A v are transition rates from length and velocity forms, respectively. The average discrepancy between the two forms was 1.5%, and in all cases, the velocity form had a larger value than the length form.
The effect of core-valence is shown in Table 4 where results are reported both for the n = 5 f and n = 7 f calculation. The former was included because of the remarkable agreement in length and velocity forms, yet the transition energy (as shown in Table 2) was not in as good agreement with the observed as before. Since the transition rate is proportional to ( Δ E ) 3 , correcting the transition rate for this factor would introduce a 4.0% reduction. Indeed, the n = 7 f transition rates were smaller with the transition energy more accurate, but length and velocity were not in as good agreement.
In summary, valence correlation predicted the best transition energy and the best agreement in length and velocity for accurate transition energy.

5. Summary and Conclusions

Table 5 compares the predicted transition rates (based on observed transition energies or wavelengths in a vacuum, rather than computed transition energies as in Table 4), with values derived from observations by Hermann et al. [1] and the values reported by Kurucz and Peytremann [2]. The latter used a semi-empirical approach in which Slater parameters were determined empirically from observed energy levels and transition probabilities calculated by the use of scaled Thomas–Fermi–Dirac wave functions. As seen in Table 5, the present predicted transition rates are about 10% smaller than observed values, whereas the Kurucz and Peytremann values are about a half those of the observed rates. Thus, the discrepancy between theory and experiment has been reduced significantly.

Author Contributions

The authors C.F.F. and J.F.B. contributed jointly to conceptualization, methodology, validation, formal analysis, writing—original draft preparation, and writing—review and editing.

Funding

This research was funded by Canada NSERC Discovery Grant 2017-03851 (CFF) and US NSF Grant No. PHY-1607396 (JFB). The APC was funded by MDPI.

Acknowledgments

The authors (CFF and JB) acknowledge support, respectively, from the Canada NSERC Discovery Grant 2017-03851 and from ITAMP, which is supported in part by Grant No. PHY-1607396 from the NSF to Harvard University and the Smithsonian Astrophysical Observatory.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Hermann, J.; Lorusso, A.; Perrone, A.; Strafella, F.; Dutouquet, C.; Torralba, B. Simulation of emission spectra from nonuniform reactive laser-induced plasmas. Phys. Rev. E 2015, 92, 053103. [Google Scholar] [CrossRef] [PubMed]
  2. Kurucz, R.L.; Peytremann, E. Research in Space Science, SAO Special Report No. 362; Smithsonian Institution; Astrophysical Observatory: Cambridge, MA, USA, 1975. [Google Scholar]
  3. Papoulia, A.; Ekman, J.; Jönsson, P. Extended transition rates and lifetimes in Al I and Al II from systematic multiconfiguration calculations. Astron. Astrophys. 2019, 621, A16. [Google Scholar]
  4. Eriksson, K.B.S.; Isberg, H.B.S. The spectrum of atomic aluminium, Al I. Ark. Fys. 1963, 23, 527. [Google Scholar]
  5. Fischer, C.F.; Godefroid, M.; Brage, T.; Jönsson, P.; Gaigalas, G. Advanced multiconfiguration methods for complex atoms: I. Energies and wave functions. J. Phys. B At. Mol. Opt. Phys. 2016, 49, 182004. [Google Scholar] [CrossRef] [Green Version]
  6. Fischer, C.F.; Gaigalas, G.; Jönsson, P.; Bieroń, J. GRASP2018-A Fortran 95 version of the General Relativistic Atomic Structure Package. Comput. Phys. Commun. 2019, 237, 184–187. [Google Scholar] [CrossRef]
  7. Kramida, A.; Ralchenko, Y.; Reader, J.; NIST ASD Team. NIST Atomic Spectra Database (ver. 5.5.6); National Institute of Standards and Technology: Gaithersburg, MD, USA, 2019. Available online: https://physics.nist.gov/asd (accessed on 15 May 2019).
Table 1. Convergence of the energy level structure for a valence correlation calculation is compared with observed data [4,7]. All results are in cm 1 .
Table 1. Convergence of the energy level structure for a valence correlation calculation is compared with observed data [4,7]. All results are in cm 1 .
n = 3 n = 4 n = 5 f n = 6 f Observed
3 s 3 p 2 4 P Fine structure
3 s 3 p 2 4 P 1 / 2  00000.00
3 s 3 p 2 4 P 3 / 2  4646454546.55
3 s 3 p 2 4 P 5 / 2  122121120120122.37
3 s 3 p 4 s 4 P o Fine structure
3 s 3 p 4 s 4 P 1 / 2 o  00000.00
3 s 3 p 4 s 4 P 3 / 2 o  5758565556.10
3 s 3 p 4 s 4 P 5 / 2 o  156157153150152.08
Multiplet separation
3 s 3 p 4 s 4 P 1 / 2 o 3 s 3 p 2 4 P 1 / 2 33,21333,02932,72432,69632,671.05
Table 2. Convergence of the energy level structure from a core-valence plus valence correlation calculation is compared with observed data [4,7]. All results are in cm 1 .
Table 2. Convergence of the energy level structure from a core-valence plus valence correlation calculation is compared with observed data [4,7]. All results are in cm 1 .
n = 5 f n = 6 f n = 7 f Observed
3 s 3 p 2 P Fine structure
3 s 3 p 2 4 P 1 / 2  0000.00
3 s 3 p 2 4 P 3 / 2  48454546.55
3 s 3 p 2 4 P 5 / 2  128120120122.37
3 s 3 p 4 s 4 P o Fine structure
3 s 3 p 4 s 4 P 1 / 2 o  0000.00
3 s 3 p 4 s 4 P 3 / 2 o  60565656.10
3 s 3 p 4 s 4 P 5 / 2 o  167155153152.08
Multiplet separation
3 s 3 p 2 4 P 1 / 2 3 s 3 p 4 s 4 P 1 / 2 o 33,142.5832,934.2132,838.3332,671.05
Table 3. Ab initio electric dipole (E1) transition data for the 3 s 3 p 2 4 P to 3 s 3 p 4 s 4 P o transition computed in the length form from valence correlation results.
Table 3. Ab initio electric dipole (E1) transition data for the 3 s 3 p 2 4 P to 3 s 3 p 4 s 4 P o transition computed in the length form from valence correlation results.
UpperLower Δ E (cm 1 ) λ (nm)A ( μ s 1 ) gf dT
3 s 3 p 4 s 4 P 1 / 2 o 3 s 3 p 2 4 P 1 / 2 32,696305.8422029.930.08400.016
3 s 3 p 4 s 4 P 1 / 2 o 3 s 3 p 2 4 P 3 / 2 32,650306.27027149.010.41910.017
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 1 / 2 32,752305.3242175.110.41990.014
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 3 / 2 32,706305.7507423.890.13390.015
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 5 / 2 32,631306.4498280.430.45290.017
3 s 3 p 4 s 4 P 5 / 2 o 3 s 3 p 2 4 P 3 / 2 32,801304.8663454.370.45460.014
3 s 3 p 4 s 4 P 5 / 2 o 3 s 3 p 2 4 P 5 / 2 32,726305.56128125.941.05770.015
Table 4. Ab initio electric dipole (E1) transition data for the 3 s 3 p 2 4 P to 3 s 3 p 4 s 4 P o transition computed in the length form from valence and core-valence results.
Table 4. Ab initio electric dipole (E1) transition data for the 3 s 3 p 2 4 P to 3 s 3 p 4 s 4 P o transition computed in the length form from valence and core-valence results.
UpperLower Δ E (cm 1 ) λ (nm)A ( μ s 1 ) gf dT
n = 5 f
3 s 3 p 4 s 4 P 1 / 2 o 3 s 3 p 2 4 P 1 / 2 33,142301.72729.760.08120.004
3 s 3 p 4 s 4 P 1 / 2 o 3 s 3 p 2 4 P 3 / 2 33,093302.170148.190.40570.006
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 1 / 2 33,203301.17374.720.40640.001
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 3 / 2 33,154301.61523.880.13030.003
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 5 / 2 33,075302.34179.910.43800.008
3 s 3 p 4 s 4 P 5 / 2 o 3 s 3 p 2 4 P 3 / 2 33,261300.65253.830.43770.001
3 s 3 p 4 s 4 P 5 / 2 o 3 s 3 p 2 4 P 5 / 2 33,181301.373125.121.02230.003
n = 7 f
3 s 3 p 4 s 4 P 1 / 2 o 3 s 3 p 2 4 P 1 / 2 32,838304.52228.510.07930.051
3 s 3 p 4 s 4 P 1 / 2 o 3 s 3 p 2 4 P 3 / 2 32,792304.947141.940.39580.053
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 1 / 2 32,894304.00371.570.39660.049
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 3 / 2 32,848304.42622.850.12700.050
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 5 / 2 32,773305.12076.550.42740.055
3 s 3 p 4 s 4 P 5 / 2 o 3 s 3 p 2 4 P 3 / 2 32,946303.52351.610.42770.047
3 s 3 p 4 s 4 P 5 / 2 o 3 s 3 p 2 4 P 5 / 2 32,871304.213119.860.99790.050
Table 5. Comparison of the transition rates computed from valence correlation calculations (present) and observed wavelengths (in a vacuum) from NIST [4,7] with observed rates from Hermann et al. [1] and values reported by Kurucz and Peytremann [2]
Table 5. Comparison of the transition rates computed from valence correlation calculations (present) and observed wavelengths (in a vacuum) from NIST [4,7] with observed rates from Hermann et al. [1] and values reported by Kurucz and Peytremann [2]
UpperLower λ (nm)A ( μ s 1 )
NISTPresentPresentHermannKurucz
3 s 3 p 4 s 4 P 1 / 2 o 3 s 3 p 2 4 P 1 / 2 305.9924305.842229.8934.18.3
3 s 3 p 4 s 4 P 1 / 2 o 3 s 3 p 2 4 P 3 / 2 306.4290306.2703148.80160.89.2
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 1 / 2 305.4679305.324274.9978.44.9
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 3 / 2 305.9029305.750723.8528.14.2
3 s 3 p 4 s 4 P 3 / 2 o 3 s 3 p 2 4 P 5 / 2 306.6144306.449880.3290.47.7
3 s 3 p 4 s 4 P 5 / 2 o 3 s 3 p 2 4 P 3 / 2 305.0073304.866354.2859.32.1
3 s 3 p 4 s 4 P 5 / 2 o 3 s 3 p 2 4 P 5 / 2 305.7144305.5613125.75140.75.0

Share and Cite

MDPI and ACS Style

Froese Fischer, C.; Babb, J.F. Transition Rates for 3s3p2 4P–3s3p4s 4Po Transitions in Al i. Atoms 2019, 7, 54. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms7020054

AMA Style

Froese Fischer C, Babb JF. Transition Rates for 3s3p2 4P–3s3p4s 4Po Transitions in Al i. Atoms. 2019; 7(2):54. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms7020054

Chicago/Turabian Style

Froese Fischer, Charlotte, and James F. Babb. 2019. "Transition Rates for 3s3p2 4P–3s3p4s 4Po Transitions in Al i" Atoms 7, no. 2: 54. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms7020054

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop