Next Article in Journal / Special Issue
Charlotte Froese Fischer—Her Work and Her Impact
Previous Article in Journal / Special Issue
The Belgian Repository of Fundamental Atomic Data and Stellar Spectra (BRASS)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coulomb (Velocity) Gauge Recommended in Multiconfiguration Calculations of Transition Data Involving Rydberg Series

1
Department of Materials Science and Applied Mathematics, Malmö University, SE-20506 Malmö, Sweden
2
Division of Mathematical Physics, Department of Physics, Lund University, SE-22100 Lund, Sweden
3
Institute of Theoretical Physics and Astronomy, Vilnius University, Saulėtekio av. 3, LT-10222 Vilnius, Lithuania
4
Chimie Quantique et Photophysique, Université libre de Bruxelles, B-1050 Brussels, Belgium
5
Hebei Key Lab of Optic-electronic Information Materials, College of Physics Science and Technology, Hebei University, Baoding 071002, China
*
Author to whom correspondence should be addressed.
Atoms 2019, 7(4), 106; https://doi.org/10.3390/atoms7040106
Submission received: 23 October 2019 / Revised: 16 November 2019 / Accepted: 21 November 2019 / Published: 26 November 2019

Abstract

:
Astronomical spectroscopy has recently expanded into the near-infrared (nIR) wavelength region, raising the demands on atomic transition data. The interpretation of the observed spectra largely relies on theoretical results, and progress towards the production of accurate theoretical data must continuously be made. Spectrum calculations that target multiple atomic states at the same time are by no means trivial. Further, numerous atomic systems involve Rydberg series, which are associated with additional difficulties. In this work, we demonstrate how the challenges in the computations of Rydberg series can be handled in large-scale multiconfiguration Dirac–Hartree–Fock (MCDHF) and relativistic configuration interaction (RCI) calculations. By paying special attention to the construction of the radial orbital basis that builds the atomic state functions, transition data that are weakly sensitive to the choice of gauge can be obtained. Additionally, we show that the Babushkin gauge should not always be considered as the preferred gauge, and that, in the computations of transition data involving Rydberg series, the Coulomb gauge could be more appropriate for the analysis of astrophysical spectra. To illustrate the above, results from computations of transitions involving Rydberg series in the astrophysically important C IV and C III ions are presented and analyzed.

1. Introduction

The starlight emitted at optical or shorter wavelengths is efficiently scattered by intervening interstellar and intergalactic dust particles. To observe stars deeper into the galactic center and go even beyond the Milky Way, astrophysical missions and spectrographs were recently designed to observe nIR radiation, which has higher transmission through dust clouds [1,2,3]. Accurate transition data from the IR part of the spectrum are thus required to interpret the spectra of distant astronomical objects observed, and to carry out chemical abundance studies.
The interest in the nIR region is relatively recent, and atomic data corresponding to wavelengths from 1 to 5 μ m are scarce. Due to the limited resources and the numerous possible transitions, laboratory measurements are insufficient to provide astrophysicists with complete sets of atomic transition data. Critically evaluated theoretical data are, therefore, necessary to complement experiments and to allow for accurate chemical abundance analyses of stars. In the long wavelength IR regime, lines of atoms are produced by transitions between states lying close in energy, which often correspond to transitions between highly excited states. The latter instance necessitates atomic structure calculations over a large portion of a spectrum. Extensive spectrum calculations that produced transition data in the nIR region were formerly carried out as part of the Opacity Project [4]. The latter non-relativistic calculations were based on the close-coupling approximation of the R-matrix theory.
Performing spectrum calculations, in which multiple atomic states are targeted at the same time, is generally not trivial. In multiconfiguration calculations, the correlation between the electrons is taken into account by expanding the targeted states in a number of symmetry adapted basis functions, which are built from products of spin-orbitals. To accurately predict the energies of all the targeted states, the shapes of the radial parts of the spin-orbitals must be such that they account for the L S -term dependencies; i.e., the way the electrons are coupled to form different terms from the same configuration [5]. Additionally, many studies involve states that are part of Rydberg series. Perturbers often enter the Rydberg series and the atomic state expansions must correctly predict their positions [6]. Computations of Rydberg series have to further describe states with electron distributions occupying different regions in space, extending far out from the atomic core. The above challenges require that special attention is paid to the optimization scheme of the wave functions; i.e., how the orbital basis is generated. The challenges in the computations of Rydberg series become more apparent when computing transition data.
The transition parameters (line strengths, oscillator strengths, and transition rates) are expressed in terms of reduced matrix elements of the transition operator. Different choices of gauge, Babushkin and Coulomb, for the transition operator lead to alternative expressions for the reduced matrix elements, and consequently, the transition parameters. Gauge invariance of the transition data is a straightforward matter for hydrogenic systems. Yet, the use of approximate wave functions results in different values for transition data expressed in different gauges. During the past years, recommendations for choosing the appropriate gauge became contradictory, suggesting further work in the field [7,8,9,10,11,12]. The Babushkin gauge (or length form) is sensitive to the outer part of the wave functions that governs the atomic transitions, and transition data expressed in this gauge are often considered to be more reliable than transition data expressed in the Coulomb gauge (or velocity form) [13]. It is, however, argued that provided reasonably accurate approximate wave functions, the Coulomb gauge (or velocity form) may give the best results when the transition energy is not very small [14]. Recent work suggests that the Coulomb gauge gives more accurate results and is the preferred gauge for transitions involving high Rydberg states [15].
In this paper, we present and analyze results from computations of Rydberg series in the C IV and C III ions. Although the latter are of astrophysical interest, the goal of the paper is not benchmarking transition data for these two ions against other theoretical methods, but instead assessing the relative reliability of the MCDHF/RCI results obtained with the two different gauges. Using the MCDHF method, we apply different computational strategies for optimizing the radial orbital basis used for constructing the wave functions and compare the results. For transitions involving low-lying states, the transition data are accurately computed in both the Babushkin and the Coulomb gauge, independently of how the radial orbitals are optimized. On the other hand, transitions involving high Rydberg states are problematic, and the Babushkin gauge does not provide trustworthy results when conventional optimization strategies are applied. However, by paying special attention to the construction of the radial orbital basis that builds the atomic state functions, transition data that are weakly sensitive to the choice of gauge are produced for all the computed transitions in the ions we study. The present article is an extended transcript of the poster presentation given on 24 June 2019 at the 13th International Colloquium on Atomic Spectra and Oscillator Strengths (ASOS2019) for Astrophysical and Laboratory Plasmas that took place at Fudan University in Shanghai, China (https://asos2019.fudan.edu.cn).

2. Theory

2.1. MultiConfiguration Calculations

Numerical representations of atomic state functions (ASFs), which are approximations to the exact wave functions, are obtained using the fully relativistic MCDHF method [16,17]. In the MCDHF method, the ASFs Ψ ( γ π J M J ) are expanded over N C S F antisymmetrized basis functions Φ ( γ ν π J M J ) , which are known as configuration state functions (CSFs), i.e.,
Ψ ( γ π J M J ) = ν = 1 N C S F c ν Φ ( γ ν π J M J ) .
In the expression above, J and M J are the angular momentum quantum numbers, π is the parity, and γ ν denotes other appropriate labeling of the CSF ν , such as orbital occupancy and angular coupling tree. The CSFs are coupled products of one-electron Dirac orbitals ψ n κ , m , which have the general form:
ψ n κ , m ( r ) = 1 r P n κ ( r ) χ κ , m ( θ , φ ) i Q n κ ( r ) χ κ , m ( θ , φ ) ,
where χ ± κ , m ( θ , φ ) are the two-component spin-angular functions and { P n κ ( r ) , Q n κ ( r ) } are, respectively, the radial functions of the large and small components, which are represented on a logarithmic grid. The selection of the CSFs to be included in the expansion (1) depends on the shell structure of the atom at hand and the computed properties, as explained in Section 3. The shape of the radial functions { P n κ ( r ) , Q n κ ( r ) } is determined by the effective field in which the considered electron moves, which is in turn established by the included CSFs [18].
The expansion coefficients c ν , together with the radial parts of the spin-orbitals, are obtained in a self-consistent field (SCF) procedure. The set of SCF equations to be iteratively solved results from applying the variational principle on a weighted energy functional of all the targeted atomic states according to the extended optimal level (EOL) scheme [19]. In fully relativistic calculations, the energy functional is estimated from the expectation value of the Dirac-Coulomb Hamiltonian [17]. The angular integrations needed for the construction of the energy functional are based on the second quantization formalism in the coupled tensorial form [20,21].
The MCDHF method is employed to generate an orbital basis. Given this basis, the final wave functions Ψ ( γ π J M J ) of the targeted states are determined in subsequent RCI calculations. In the RCI calculations, the spin-orbitals defining the basis are fixed and only the expansion coefficients c ν are evaluated by diagonalizing the Hamiltonian matrix. At this step, the expansions based on Equation (1) can be augmented to include CSFs that account for additional electron correlation effects. All MCDHF and RCI calculations were performed using the relativistic atomic structure package GRASP2018 [22].

2.2. Transition Parameters

Once the wave functions Ψ ( γ π J M J ) have been determined, transition parameters can be computed. In this work, we focus on the computation of transition rates (or probabilities) for electric dipole (E1) transitions. Electric dipole transitions are much stronger than electric quadrupole (E2) and magnetic multipole (Mk) transitions. For the transition rate A ( k ) of electric dipole ( k = 1 ) emission from an upper state γ π J M J to any of the 2 J + 1 states γ π J M J of lower energy, we have the following proportionality
A ( 1 ) ( γ π J , γ π J ) ( E γ π J E γ π J ) 3 S ( 1 ) ( γ π J , γ π J ) 2 J + 1 ,
where E γ π J E γ π J is the transition energy and S ( 1 ) ( γ π J , γ π J ) is the line strength given by
S ( 1 ) ( γ π J , γ π J ) = | Ψ ( γ π J ) | | O ( 1 ) | | Ψ ( γ π J ) | 2 .
The E1 transition rates are therefore expressed in terms of reduced matrix elements of the electric dipole transition operator O ( 1 ) . From Equation (1), it follows that
Ψ ( γ π J ) | | O ( 1 ) | | Ψ ( γ π J ) = k , l c k c l Φ ( γ k π J ) | | O ( 1 ) | | Φ ( γ l π J ) .
The choice of gauge parameter determines whether the electric dipole matrix elements are computed in the Babushkin or the Coulomb gauge, which in non-relativistic calculations correspond to the length and the velocity form, respectively [10]. The two forms are equivalent for hydrogenic wave functions, but they result in different values when approximate many-electron wave functions are used. As shown later, in Section 4, these values reveal a strong dependence on the generated orbital basis and the captured correlation effects. Although the present results arise from fully relativistic calculations, similar behavior is observed when non-relativistic multiconfiguration calculations are performed [15].
The explicit expressions of the electric dipole reduced matrix elements in the Babushkin and the Coulomb gauge are given in [10]. Taking for convenience the non-relativistic limit, the electric dipole reduced matrix elements are, in the length and the velocity form, respectively, given by
Ψ ( γ π J ) | | i = 1 N r i C ( 1 ) ( i ) | | Ψ ( γ π J )
and
1 E γ π J E γ π J Ψ ( γ π J ) | | i = 1 N ( 1 ) ( i ) | | Ψ ( γ π J ) ,
where the summation runs over the number N of electrons and C ( 1 ) is the renormalized spherical harmonic of rank 1 [23]. The reduced matrix elements of (6) and (7) involve, respectively, sums over radial transition integrals of the kind
0 P ( r ) r P ( r ) d r
and
0 P ( r ) d d r P ( r ) d r ,
weighted with the products of the expansion coefficients of the CSFs and the angular coefficients [20,21,23]. P ( r ) and P ( r ) are the radial functions of the large components of the Dirac one-electron spin-orbitals (2) that build the CSFs of the initial state γ π J and the final state γ π J , respectively. In the present work, the initial and final states belonging to different parities are built from a common orbital basis.
In the computation of the integrals (8), the emphasis is given to the tail of the radial orbitals, while in the integrals (9) the emphasis is given instead to the inner part of the radial orbitals. In the simple Hartree–Fock (HF) model, the approximate wave functions usually display a correct asymptotic behavior towards large r (see also Section 5), and since the former integrals are also computationally simpler, the transition rates are traditionally provided in the length form [7,8,9,14]. As discussed in Section 5, when multiconfiguration methods concurrently target multiple atomic states, all wave functions are not always well approximated at large r; and the velocity form, or correspondingly, the Coulomb gauge, may by contrast, give the best results.
The agreement between the transition rates A B and A C , respectively, evaluated in the Babushkin and the Coulomb gauge is used as an indicator of accuracy. This is particularly useful when laboratory measurements are not available for comparison. The uncertainty of the computed transition rates in the preferred gauge can be estimated as
d T = | A B A C | max ( A B , A C ) ,
which reflects the relative discrepancy between the Babushkin and the Coulomb gauge of the computed line strengths [24,25]. Although the accuracy indicators d T should be used in a statistical manner for a group of transitions with similar properties (see [6]), individual d T values can point out problematic transitions, which could further be analyzed.

3. Computational Methodology—Optimization of the Orbital Basis

The accuracy of multiconfiguration calculations relies on the CSF expansion of Equation (1). A first approximation of the ASFs is acquired by performing an MCDHF calculation on expansions that are built from the configurations that define what is known as the multi-reference (MR) [17]. The orbitals that take part in this initial calculation are called spectroscopic orbitals and are kept frozen in all subsequent calculations. The initial approximation of the ASFs is improved by augmenting the expansion with CSFs that interact with the ones that are generated by the MR configurations. Such CSFs are built from configurations that differ by either a single (S) or a double (D) electron substitution from the configurations in the MR [17,26]. Following the SD-MR scheme, the interacting configurations are obtained by allowing substitutions of electrons from the spectroscopic orbitals to an active set of correlation orbitals, which is systematically increased (each step introducing an additional correlation orbital layer) [27,28]. These configurations produce CSFs that can be classified, based on the nature of the substitutions, into CSFs that capture valence–valence (VV), core–valence (CV), and core–core (CC) electron correlation effects.
Building accurate wave functions requires a very large orbital basis. Even so, a large but incomplete orbital basis does not ensure that the wave functions give accurate properties other than energies. In the MCDHF calculations, the correlation orbitals are obtained by applying the variational principle on the weighted energy functional of all the targeted atomic states. Thus, the orbitals of the first correlation layers will overlap with the spectroscopic orbitals that account for the effects that minimize the energy the most [18,29]. The energetically dominant effects must first be saturated to obtain orbitals localized in other regions of space, which might describe effects that do not lower the energy much, but are important for, e.g., transition parameters. One must, therefore, carefully choose the orbital basis with respect to the computed properties [30].
Valence atomic transitions are governed by the outer part of the wave functions and this part must be properly described by the correlation orbitals to obtain reliable transition parameters. States that are part of Rydberg series encompass valence orbitals of increasing principal quantum number n. Spectrum calculations that involve Rydberg series need, thus, to describe states with electron distributions localized in different regions of space extending far out from the atomic core. Since the overlap between orbitals describing Rydberg states is in some cases minor, generating an optimal orbital basis is not straightforward [6]. This raises the need to explore different computational strategies.

3.1. C IV

In lithium-like carbon, the configurations being studied are 1 s 2 n l with n = 2 to 8 and l = 0 to 4 and 1 s 2 6 h . These configurations define the MR and correspond to 53 targeted atomic states of both even and odd parity, which are simultaneously optimized. For simple systems such as three-electron systems, the MCDHF calculations are conventionally performed using CSF expansions that are produced by SD-MR electron substitutions from all spectroscopic orbitals. In this manner, the CSFs capture all valence (V), CV, and CC correlation effects. The 1 s 1 s pair-correlation effect is energetically very important and the orbitals of the first correlation layers overlap with the 1 s core orbital accounting for this effect (see Table 1). After building six correlation layers by utilizing this conventional approach, we see that all correlation orbitals up to 14 s , 14 p , 14 d , 12 f , 12 g , 8 h , and 7 i are rather contracted in comparison with the outer Rydberg orbitals. As a consequence, the wave functions are not properly described for all states, and in particular, not for the higher Rydberg states considered.
For a more appropriate description of the wave functions, the correlation orbitals must occupy the space between the 1 s core orbital and the inner valence orbitals. This can be accomplished by imposing restrictions on the allowed substitutions for obtaining the orbital basis. Thus, the MCDHF calculations are alternatively performed using CSF expansions that are produced by SD-MR substitutions with the restriction of allowing maximum one hole in the 1 s core shell. In this case, the shape of the correlation orbitals is established by CSFs accounting for V and CV correlation effects. The resulting correlation orbitals are, as shown in Table 1, more extended, overlapping with orbitals of higher Rydberg states.
The final wave functions of the targeted states are determined in subsequent RCI calculations, where D substitutions from the 1 s core orbital and triple (T) substitutions from all the spectroscopic orbitals, are included. The number of CSFs in the final even and odd state expansions are, respectively, 1,077,872 and 1,287,706, distributed over the different J symmetries.

3.2. C III

In beryllium-like carbon, the configurations in question are 1 s 2 2 s n l with n = 2 to 7 and l = 0 to 4 and 1 s 2 2 p 2 , 1 s 2 2 p 3 s , 1 s 2 2 p 3 p , and 1 s 2 2 p 3 d . These configurations define the MR and correspond to 114 targeted atomic states of both even and odd parity, which are simultaneously optimized. Having introduced two correlation orbitals, 8 s and 8 p —specifically targeted to account for the L S -term dependence [31], i.e., the difference between the n s orbitals for 2 s n s 3 S and 2 s n s 1 S and the difference between the n p orbitals for 2 s n p 3 P and 2 s n p 1 P —the MCDHF calculations are conventionally performed using CSF expansions that are produced by SD-MR electron substitutions from all spectroscopic orbitals with the restriction that only one excitation is allowed from the 1 s 2 atomic core. In this manner, the CSFs capture VV and CV correlation effects. The 1 s n l pair-correlation effect is comparatively important, and the orbitals of the first correlation layers are spatially localized between the 1 s orbital and the 2 s and 2 p orbitals. As the CV correlation effects start to saturate, the correlation orbitals are gradually located further away from the 1 s 2 atomic core (see Table 2). The correlation orbitals up to 12 s , 12 p , 12 d , 12 f , 11 g , and 8 h are, however, still contracted in comparison with the outer Rydberg orbitals.
For a more appropriate description of the wave functions, the correlation orbitals must occupy the space of the valence orbitals. In the alternative approach, this is accomplished by allowing SD substitutions only from the outer valence orbitals accounting for VV correlation. The resulting correlation orbitals are, as shown in Table 2, more extended, overlapping with orbitals of higher Rydberg states.
The final wave functions of the targeted states are determined in subsequent RCI calculations, where SDT substitutions from all the spectroscopic orbitals are included, with the restriction that only one substitution is allowed from the 1 s 2 atomic core. The number of CSFs in the final even and odd state expansions are, respectively, 1,578,620 and 1,274,147, distributed over the different J symmetries.

4. Results

Excitation energies are produced, based on the conventional and alternative computational strategies that were described in Section 3, and are compared with the critically evaluated data from the National Institute of Standards and Technology’s (NIST’s) Atomic Spectra Database (ASD) [32]. In C IV, the computed excitation energies are in excellent agreement with the NIST’s recommended values. Both computational approaches give similar energies and the relative differences from the NIST values are less than 0.01 % . For the more complex system of C III, the computed excitation energies agree also well with the energies proposed by NIST. The relative differences between theoretical and critically compiled energies are, on average, of the order of 0.1 % and 0.02 % , when following the conventional and the alternative approach, respectively. The NIST database does not provide excitation energies for the 2 s 6 s 2 S , 2 s 7 s 2 S , and 2 s 7 p 3 P states, which are included in the computations.
Transition rates A are produced based on the two different computational strategies. In the present computations, the uncertainties in the predicted excitation energies of two states associated with a transition most often cancel out, and consequently, the majority of the evaluated transition energies are ultimately in perfect agreement with the NIST values. The uncertainties of the computed transition rates solely emerge from the disagreement of the computed line strengths in the Babushkin and the Coulomb gauge, which are then reflected in the d T values. When the conventional strategy is applied, most of the transition rates are predicted with uncertainties d T lower than 1 % and 5 % , in lithium-like and beryllium-like carbon, respectively. Yet, for transitions involving high Rydberg states, the uncertainties increase remarkably, especially for the more complex C III ion. The alternative strategy for optimizing the radial orbitals yields transition rates that are overall more accurate. The improvement in accuracy is significant for transitions that involve high Rydberg states.
The uncertainties d T of the transition rates computed with the conventional and alternative approaches are presented and analyzed for groups of transitions in the studied carbon ions. Each group is selected to include transitions between a fixed state and Rydberg states described by electron distributions that are gradually localized farther from the atomic core. Accordingly, Figure 1a,b illustrates the uncertainties d T for the 2 p 2 P 1 / 2 n s 2 S 1 / 2 and n p 2 P 1 / 2 8 s 2 S 1 / 2 groups of transitions in C IV. Similarly, Figure 2a,b illustrates the d T values for the 2 s 2 1 S 0 2 s n p 1 P 1 and 2 s n s 1 S 0 2 s 7 p 1 P 1 groups of transitions in C III.
Figure 1a demonstrates the uncertainties for the series of transitions between the low-lying 2 p 2 P 1 / 2 state and successively higher Rydberg n s 2 S 1 / 2 states. The uncertainty d T of the transition rates computed with the conventional approach grows almost exponentially with the increasing principal quantum number n. The transition rate for the 2 p 2 P 1 / 2 8 s 2 S 1 / 2 transition, which is the transition between the two states with the largest energy difference in the plot, eventually exhibits the highest uncertainty— 12.4 % . When the alternative approach is utilized instead, the uncertainties range between 0 % and 0.4 % for the respective transitions. The same trends are also observed in other groups of transitions in C IV, such as the 2 p 2 P n d 2 D , the 2 s 2 S n p 2 P , and so forth.
Having as a starting point the 2 p 2 P 1 / 2 8 s 2 S 1 / 2 transition, Figure 1b demonstrates the uncertainties for the series of transitions between the high Rydberg 8 s 2 S 1 / 2 state and successively higher Rydberg n p 2 P 1 / 2 states. As n increases, the transition energy gets smaller. The uncertainties of the transition rates computed with the conventional approach exhibit a nearly exponential decay with increasing n. Similarly to Figure 1a, when following the alternative strategy, the uncertainties in the transition rates are substantially reduced, ranging between 0.1 % and 0.4 % . Other groups of transitions in C IV, such as the n p 2 P 8 d 2 D series and the n s 2 S 8 p 2 P series, follow analogous trends.
Figure 2a displays the d T values for transitions between the low-lying 2 s 2 1 S 0 state and successively higher Rydberg 2 s n p 1 P 1 states. Looking at Figure 2a, when the conventional approach is applied the uncertainties d T increase sharply for n > 5 . For the 2 s 2 1 S 0 2 s 7 p 1 P 1 transition, i.e., the transition between the two states with the largest energy difference, the d T rises to 63 % . The latter is about five times larger than the highest estimated d T value in the figures above. Once more, when the radial orbitals are optimized using the alternative strategy, the uncertainties drop dramatically, ranging between 0 % and 1.4 % for the respective transitions. More groups of transitions in C III that reveal similar behavior are the 2 s 2 p 1 P 1 2 s n s 1 S 0 series and the 2 p 2 3 P 0 2 s n p 3 P 1 series.
Starting with the 2 s 2 1 S 0 2 s 7 p 1 P 1 transition, Figure 2b displays the d T values for transitions between the high Rydberg 2 s 7 p 1 P 1 state and successively higher Rydberg 2 s n s 1 S 0 states. Likewise, in Figure 1b, the increase in n corresponds to transitions between states that gradually come closer in energy. The uncertainties of the transition rates computed with the conventional approach reduce rapidly as n increases. Applying the alternative strategy results in much lower uncertainties, which extend between 0.5 % and 6.6 % . A similar trend is also observed in the 2 s n p 1 P 1 2 s 7 s 1 S 0 series of transitions in C III. Although the last two points in Figure 2b correspond to transitions between states lying close in energy, the uncertainties are comparatively high. Nevertheless, the alternative strategy still predicts the transition rates with lower uncertainties.
Altogether, for transitions between low-lying states, the transition rates are accurately predicted independently of whether the conventional or the alternative computational strategy is employed. Further, when both states involved in a transition in C IV are high Rydberg states, the transition rates are also predicted with high accuracy in both computations. On average, the same holds for transitions between high Rydberg states in C III. The line strengths of transitions between two states close in energy, and with the outer electrons occupying nearly the same region of space, are relatively large, and therefore, weakly affected by the optimization strategy of the radial orbitals. Quite the contrary, transitions between a low-lying state and a high Rydberg state are problematic in both carbon ions. The line strength of transitions between two states with large energy differences, and with the outer electrons occupying different parts of space, take smaller values, which are more sensitive to how the radial orbitals are optimized with regard to correlation.
To better understand the origins of the large d T values in transitions between low-lying states and high Rydberg states, the convergences of the individual transition rates A B and A C , computed in the Babushkin and the Coulomb gauges respectively, are studied with respect to the increasing active set of correlation orbitals. In connection with the figures above, this is done for the 2 p 2 P 1 / 2 8 s 2 S 1 / 2 transition in C IV (see Figure 3a) and 2 s 2 1 S 0 2 s 7 p 1 P 1 transition in C III (see Figure 3b); i.e., the transitions with the highest uncertainties d T . In Figure 3a,b, the convergences of the A B and A C values are illustrated for the two different computational approaches.
As seen in Figure 3a,b, when the computations are performed in the alternative manner, the transition rates given by A B and A C ultimately come really close in value. Considering the small final d T value, the agreement between the A B and A C values is expected. One observes that the transition rate given by A C is rather stable with respect to the increasing orbital set. The A C value varies by only 1 % and 6.2 % for each of the transitions displayed in the figures below. On the contrary, the A B value varies by 13.7 % and 27.2 % , respectively. In the 2 p 2 P 1 / 2 8 s 2 S 1 / 2 transition, it takes five correlation layers for the A B value to start converging, while in the 2 s 2 1 S 0 2 s 7 p 1 P 1 transition it takes three layers of correlation orbitals for the A B value to converge.
Looking at Figure 3a,b, when the conventional approach is applied, the individual A B and A C values do not converge, as the large final uncertainties d T reveal. The transition rate given by A C is, however, again stable and is also consistent with the A B and A C values provided by the alternative computational strategy. Throughout the optimization of the radial orbitals in the conventional manner, the A C value varies by only 0.9 % and 2.3 % for each of the transitions displayed in Figure 3a,b, respectively. Although it seems that the A B values will eventually approach the A C ones, this would require a very large orbital basis, which is beyond the reach of the available computational resources. One may deduce that when the conventional computational strategy is applied the transition rates A B , computed in the Babushkin gauge, are problematic and unreliable.

5. Discussion

Transition data, such as transition rates A, are expressed in terms of reduced matrix elements of the transition operator (see Equation (5)), which can be computed in different gauges. According to Equation (10), the uncertainty of the computed A values is assessed by the agreement of the transition rates computed in the different gauges. Computations of reduced matrix elements in different gauges, however, probe separate parts of the wave functions. Hence, the radial parts of the wave functions must be well approximated as a whole to obtain gauge invariant transition rates.
For transitions between low-lying states, both computational strategies yield reduced matrix elements of the transition operator that almost reach gauge invariance, and the transition rates are, overall, accurately predicted. There are enough correlation orbitals spatially localized between the core and the inner valence orbitals that make up the low-lying states. As a result, the inner parts of the wave functions are adequately approximated. Moreover, the spectroscopic outer valence orbitals, which make up the higher Rydberg states and are localized farther from the atomic core, improve the description of the outer parts of the wave functions for representing the low-lying states, ensuring that they have the correct asymptotic behavior. The radial parts of the latter wave functions are then effectively described at all r values, being insensitive to the choice of the optimization strategy with regard to correlation.
For transitions between a low-lying state and a high Rydberg state, the conventional computational strategy fails to produce accurate transition rates. The correlation orbitals are significantly contracted compared to the outer Rydberg orbitals. Further, there are no spectroscopic orbitals farther localized to correct for the fact that the asymptotic behavior of the tail of the wave functions that represent the higher Rydberg states is not well approximated. Thus, the Babushkin gauge that probes the outer part of the wave functions does not produce trustworthy results. The inner parts of the wave functions representing the higher Rydberg states are, however, adequately approximated, and as a result, the Coulomb gauge yields transition rates that are more reliable (see also, Figure 3a,b).
The alternative computational strategy generates correlation orbitals that are more extended, increasing the overlap with the spectroscopic orbitals that make up the higher Rydberg states. In this case, the correlation orbitals are properly localized to ably describe the asymptotic behavior of the outer part of the wave functions representing the higher Rydberg states. That being so, after the final MCDHF and RCI computations in the alternative manner, the reduced matrix elements of the transition operator are practically gauge invariant and the transition rates are also accurately predicted for the transitions between low-lying states and high Rydberg states (see also, Figure 3a,b).
The radial transition integrals (8) and (9) that take part in the computations of the reduced matrix elements of the transition operator have an upper integration bound that goes to infinity. In (8) and (9), P ( r ) and P ( r ) are the radial parts of the spectroscopic and correlation orbitals that are included in the computations. If we express the transition integrals as a function of the upper integration bound R, we get
0 R P ( r ) r P ( r ) d r
and
0 R P ( r ) d d r P ( r ) d r ,
respectively. We can keep R = for the spectroscopic orbitals so that they extend to their full values and only introduce a cut-off value for R in the transition integrals involving correlation orbitals. In this manner, the effect on the transition rate values, from correlation orbitals gradually localized farther from the origin, can be studied. In connection with Figure 3a, the effect that the shape of the correlation orbitals has on the computation of transition rates is, in Figure 4a, examined for the 2 p 2 P 1 / 2 8 s 2 S 1 / 2 transition in C IV. In Figure 4a, the transition rates are computed by employing the alternative computational strategy and both Babushkin and Coulomb gauges are displayed. One observes that the two gauges are affected differently by the outer parts of the correlation orbitals.
In Figure 4a, the transition rate computed in the Coulomb gauge is mainly influenced by the correlation orbitals that are localized close to the origin and in the vicinity of the atomic core. Correlation orbitals occupying regions with R > 1 have an insignificant effect on the Coulomb gauge. This explains the fact that the conventional computational strategy, which generates more contracted orbitals, still predicts with accuracy, the transition rates in the Coulomb gauge. Oppositely, the Babushkin gauge is hugely affected by the correlation orbitals occupying the region between R = 5 and R = 6 . Looking at Figure 4b, the 8 s radial orbital, which extends far out from the 1 s 2 atomic core, begins its asymptotic decay at about R = 5 and dies out at R 6 . Only when we have orbitals extending into this region, the asymptotic behavior of the wave function representing the 8 s 2 S 1 / 2 state is well described, and thus, the Babushkin gauge will yield accurate transition rates. As previously seen, the conventional computational strategy fails to do so.
A similar study was performed for a transition between two high Rydberg states. In Figure 5a, the effect of the shape of the correlation orbitals on the computed transition rates is examined for the 7 p 2 P 1 / 2 8 s 2 S 1 / 2 transition in C IV. As seen in Figure 5a, correlation has nearly the same effect on both gauges. Although the 7 p and 8 s orbitals extend far out from the atomic core (see Figure 5b), correlation orbitals occupying the large R region remain unimportant in the Babushkin gauge. For transitions between states close in energy, the line strengths take large values and the change in the transition rates due to correlation is very small. For this reason, the conventional computational strategy also yields accurate transition rates for transitions between high Rydberg states.
To clarify the fact that the asymptotic behavior of the wave function at large distances is (or is not) well approximated, depending on the alternative (or the conventional) optimization strategy, Brillouin’s theorem [33,34] can be put forward to emphasize the importance of the variational content of the wave functions. When being interested into the description of Rydberg states 1 s 2 n l 2 L ( L = l ) in the single-configuration non-relativistic HF approximation, node counting of the valence radial function, n c = n l 1 , provides a simple and efficient way to select the desired state in the self-consistent procedure [34]. Each separately optimized state implicitly contains all single-electron excitations n l n l of both lower and upper parts of the spectrum, including the continuum 1 s 2 ϵ l , with the associated interesting property
Φ HF ( 1 s 2 n l 2 L ) | H | Φ ( 1 s 2 n l 2 L ) = Φ HF ( 1 s 2 n l 2 L ) | H | Φ ( 1 s 2 ϵ l 2 L ) = 0 , ( n , ϵ )
where H is the scalar non-relativistic Hamiltonian that is used in the energy functional to derive the HF equations. The annihilation property of the ( M L , M S ) -independent interaction matrix elements between the reference HF CSF built with the optimized HF orbitals, i.e., Φ HF , and all single-electron excitation CSFs Φ , defines Brillouin’s theorem and explains the reasonable accuracy of the HF approximation through the richness of its variational content. The above discussion can be extended to the relativistic framework by considering n κ ( n , ϵ ) κ single-electron excitations and taking for H the relativistic Hamiltonian that is used for deriving the DHF equations [17].
In the present work, it is hard to define the variational content of the MCDHF approach due to the complexity of the energy functional, but one should keep in mind that the optimization strategy is based (i) on a layer-by-layer approach in which only the last layer is variational while the previous ones are kept frozen, and (ii) on the use of the EOL method targeting, simultaneously, a large number of states for a spectrum calculation. The resulting lack of variational freedom for the individual states can (partially) be counterbalanced by the inclusion of enough interacting states in the Hamiltonian matrix. Going back to the single-configuration approximation case mentioned above, any member of a Rydberg series can be described through configuration interaction involving Brillouin one-electron excitations with a resulting CI-expansion strictly equivalent to the single approximation HF wave function if the basis of single-excitation CSFs is large and rich enough. This equivalence has been exploited to solve convergence problems encountered in the MCHF study of Rydberg series in strontium [35] or to demonstrate the correspondence between different orbital optimization schemes for describing the discrete-continuum interactions in complex systems [36]. In the context of our work, one illustrates the inadequacy of the orbitals obtained in the conventional approach that are used to compensate the lack of variational freedom in the representation of the high-lying Rydberg members. On the contrary, the alternative strategy proposed produces orbitals that have a better localization for describing the single-electron excitations, which would have been implicitly included with a fully optimized MCDHF wave function targeting a single Rydberg state.

6. Summary and Conclusions

The computations of transition data in the systems of lithium-like and beryllium-like carbon are examples of spectrum calculations that involve Rydberg series. In this work, we showed that, independently of the optimization scheme of the radial orbitals, transition parameters corresponding to the lower part of the spectrum are computed with high accuracy. As astronomical spectroscopy raises the demand on atomic data, highly accurate transition parameters are, however, also required for transitions that involve high Rydberg states. We demonstrated how this can be achieved by paying special attention to the optimization scheme of the radial orbitals with respect to correlation. Finally, we showed that the Babushkin gauge should not, by default, be considered as the preferred gauge, and that, in the computations of Rydberg series, it might be required that the transition rates in the Coulomb gauge are used as a reference for the interpretation of astrophysical observations.

Author Contributions

All authors contributed jointly to conceptualization, investigation, validation, and writing—reviewing and editing. A.P. and P.J. cooperatively took part in the methodology and visualization. A.P., P.J., and J.E. conducted the formal analysis. A.P. performed the writing—original draft preparation and data curation.

Funding

This research was funded by the Swedish Research Council (VR) under contract 2015-04842 and 2016-04185. A.P. was financially supported by the Royal Physiographic Society of Lund (https://www.fysiografen.se/en/) to attend the ASOS 2019 conference in Shanghai. M.G. acknowledges support from the FWO-FNRS Excellence of Science Programme (grant number EOS-O022818F). S.S. is a FRIA grantee of the Fonds de la Recherche Scientifique—FNRS.

Acknowledgments

M.G., G.G., and P.J., were research associates at Vanderbilt University in 1981, 1994–97, and 1995–97, respectively. They are grateful to Charlotte for her endless training and encouragement in many atomic physics research projects over the years, for which her deep knowledge in computational and mathematical physics was crucial. Her outstanding contribution to theoretical atomic physics has been a boundless source of inspiration for all of us.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Käufl, H.-U.; Ballester, P.; Biereichel, P.; Delabre, B.; Donaldson, R.; Dorn, R.; Fedrigo, E.; Finger, G.; Fischer, G.; Franza, F.; et al. CRIRES: A high-resolution infrared spectrograph for ESO’s VLT. Ground-Based Instrum. Astron. 2004, 5492, 1218–1227. [Google Scholar]
  2. Young, E.T.; Becklin, E.; Marcum, P.M.; Roellig, T.L.; Buizer, J.M.D.; Herter, T.L.; Güsten, R.; Dunham, E.W.; Temi, P.; Andersson, B.-G.; et al. Early Science with SOFIA, the Stratospheric Observatory for Infrared Astronomy. Astrophys. J. Lett. 2012, 749, L17. [Google Scholar] [CrossRef]
  3. Dorn, R.J.; Anglada-Escude, G.; Baade, D.; Bristow, P.; Follert, R.; Gojak, D.; Grunhut, J.; Hatzes, A.; Heiter, U.; Hilker, M.; et al. CRIRES+: Exploring the Cold Universe at High Spectral Resolution. Messenger 2014, 156, 7–11. [Google Scholar]
  4. Atomic Databases from the Opacity Project and the Iron Project. Available online: http://cdsweb.u-strasbg.fr/topbase/testop/home.html (accessed on 13 November 2019).
  5. Galvez, F.J.; Buendia, E.; Sarsa, A. 1s22p3 and 1s22s23l, l = s, p, d, excited states of boron isoelectronic series from explicitly correlated wave functions. J. Chem. Phys. 2005, 123, 034302. [Google Scholar] [CrossRef]
  6. 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] [CrossRef]
  7. Crossley, R.J.S. The Calculation of Atomic Transition Probabilities. Adv. At. Mol. Phys. 1969, 5, 237–296. [Google Scholar]
  8. Starace, A.F. Length and Velocity Formulas in Approximate Oscillator-Strength Calculations. Phys. Rev. A 1971, 3, 1242. [Google Scholar] [CrossRef]
  9. Starace, A.F. Comment on “Length and Velocity Formulas in Approximate Oscillator-Strength Calculations”. Phys. Rev. A 1973, 8, 1141. [Google Scholar] [CrossRef]
  10. Grant, I.P. Gauge invariance and relativistic radiative transitions. J. Phys. B At. Mol. Opt. Phys. 1974, 7, 1458. [Google Scholar] [CrossRef]
  11. Grant, I.P.; Starace, A.F. Gauge invariance and radiative transition probabilities. J. Phys. B At. Mol. Opt. Phys. 1975, 8, 1999. [Google Scholar] [CrossRef]
  12. Crossley, R.J.S. Fifteen Years On—The Calculation of Atomic Transition Probabilities Revisited. Phys. Scr. 1984, T8, 117–128. [Google Scholar] [CrossRef]
  13. Hibbert, A. Oscillator strengths of transitions in the beryllium sequence. J. Phys. B At. Mol. Opt. Phys. 1974, 7, 1417. [Google Scholar] [CrossRef]
  14. Cowan, R.D. Electric dipole radiation. In The Theory of Atomic Structure and Spectra; University of California Press: London, UK, 1981; pp. 400–402. [Google Scholar]
  15. Pehlivan Rhodin, A.; Hartman, H.; Nilsson, H.; Jönsson, P. Experimental and theoretical oscillator strengths of Mg I for accurate abundance analysis. Astron. Astrophys. 2017, 598, A102. [Google Scholar] [CrossRef]
  16. Grant, I.P. Relativistic Quantum Theory of Atoms and Molecules; Springer: New York, NY, USA, 2007. [Google Scholar]
  17. Froese Fischer, C.; 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]
  18. Verdebout, S.; Jönsson, P.; Gaigalas, G.; Godefroid, M.; Froese Fischer, C. Exploring biorthonormal transformations of pair-correlation functions in atomic structure variational calculations. J. Phys. B At. Mol. Opt. Phys. 2010, 43, 074017. [Google Scholar] [CrossRef]
  19. Dyall, K.G.; Grant, I.P.; Johnson, C.T.; Parpia, F.A.; Plummer, E.P. GRASP: A general-purpose relativistic atomic structure program. Comput. Phys. Commun. 1989, 55, 425–456. [Google Scholar] [CrossRef]
  20. Gaigalas, G.; Rudzikas, Z.; Froese Fischer, C. An efficient approach for spin-angular integrations in atomic structure calculations. J. Phys. B At. Mol. Opt. Phys. 1997, 30, 3747–3771. [Google Scholar] [CrossRef]
  21. Gaigalas, G.; Fritzsche, S.; Grant, I.P. Program to calculate pure angular momentum coefficients in jj-coupling. Comput. Phys. Commun. 2001, 139, 263–278. [Google Scholar] [CrossRef]
  22. Froese Fischer, C.; 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]
  23. Froese Fischer, C.; Godefroid, M.R. Programs for computing LS and LSJ transitions from MCHF wave functions. Comput. Phys. Commun. 1991, 64, 501–519. [Google Scholar] [CrossRef]
  24. Froese Fischer, C. Evaluating the accuracy of theoretical transition data. Phys. Scr. 2009, T134, 014019. [Google Scholar] [CrossRef]
  25. Ekman, J.; Godefroid, M.; Hartman, H. Validation and Implementation of Uncertainty Estimates of Calculated Transition Rates. Atoms 2014, 2, 215–224. [Google Scholar] [CrossRef]
  26. Froese Fischer, C.; Brage, T.; Jönsson, P. Structure of ψ1. In Computational Atomic Structure—An MCHF Approach; CRC Press: Boca Raton, FL, USA, 1997; pp. 70–71. [Google Scholar]
  27. Olsen, J.; Roos, B.O.; Jørgensen, P.; Jensen, H.J.A. Determinant based configurationinteraction algorithms for complete and restricted configuration interaction spaces. J. Chem. Phys. 1988, 89, 2185. [Google Scholar] [CrossRef]
  28. Sturesson, L.; Jönsson, P.; Froese Fischer, C. JJGEN: A flexible program for generating lists of jj-coupled configuration state functions. Comput. Phys. Commun. 2007, 177, 539–550. [Google Scholar] [CrossRef]
  29. Godefroid, M.; Jönsson, P.; Froese Fischer, C. Atomic structure variational calculations in spectroscopy. Phys. Scr. 1998, T78, 33–46. [Google Scholar] [CrossRef]
  30. Froese Fischer, C.; Jönsson, P.; Godefroid, M. Some two-electron properties of sodium. Phys. Rev. A 1998, 57, 1753. [Google Scholar] [CrossRef]
  31. Froese Fischer, C.; Brage, T.; Jönsson, P. Term dependence. In Computational Atomic Structure—An MCHF Approach; CRC Press: Boca Raton, FL, USA, 1997; pp. 55–56. [Google Scholar]
  32. Kramida, A.; Ralchenko, Y.U.; Reader, J.; NIST ASD Team. NIST Atomic Spectra Database (ver. 5.6.1); National Institute of Standards and Technology: Gaithersburg, MD, USA, 2018. Available online: https://physics.nist.gov/asd (accessed on 24 September 2019).
  33. Godefroid, M.; Liévin, J.; Metz, J.-Y. Brillouin’s theorem for complex atomic configurations. J. Phys. B At. Mol. Opt. Phys. 1987, 20, 3283–3296. [Google Scholar] [CrossRef]
  34. Froese Fischer, C. The Hartree-Fock Method for Atoms. A Numerical Approach; John Wiley and Sons: New York, NY, USA, 1977. [Google Scholar]
  35. Vaeck, N.; Godefroid, M.; Hansen, J.E. Multiconfiguration Hartree-Fock calculations for singlet terms in neutral strontium. Phys. Rev. A 1988, 38, 2830–2845. [Google Scholar] [CrossRef]
  36. Cowan, R.D.; Hansen, J.E. Discrete-continuum interactions in Cl I and S I. JOSA 1981, 71, 60. [Google Scholar] [CrossRef]
Figure 1. (a) The uncertainty d T of the computed transition rates for transitions between the 2 p 2 P 1 / 2 state and Rydberg n s 2 S 1 / 2 states of increasing principal quantum number n in C IV. The black squares and magenta diamonds, respectively, correspond to the results from the conventional and the alternative strategies for optimizing the radial orbitals. (b) Same as the first panel, but for transitions between the 8 s 2 S 1 / 2 state and successive Rydberg n p 2 P 1 / 2 states in C IV.
Figure 1. (a) The uncertainty d T of the computed transition rates for transitions between the 2 p 2 P 1 / 2 state and Rydberg n s 2 S 1 / 2 states of increasing principal quantum number n in C IV. The black squares and magenta diamonds, respectively, correspond to the results from the conventional and the alternative strategies for optimizing the radial orbitals. (b) Same as the first panel, but for transitions between the 8 s 2 S 1 / 2 state and successive Rydberg n p 2 P 1 / 2 states in C IV.
Atoms 07 00106 g001
Figure 2. (a) The uncertainty d T of the computed transition rates for transitions between the 2 s 2 1 S 0 state and Rydberg 2 s n p 1 P 1 states of increasing principal quantum number n in C III. The black squares and magenta diamonds, respectively, correspond to the results from the conventional and alternative strategies for optimizing the radial orbitals. (b) Same as the first panel, but for transitions between the 2 s 7 p 1 P 1 state and successive Rydberg 2 s n s 1 S 0 states in C III.
Figure 2. (a) The uncertainty d T of the computed transition rates for transitions between the 2 s 2 1 S 0 state and Rydberg 2 s n p 1 P 1 states of increasing principal quantum number n in C III. The black squares and magenta diamonds, respectively, correspond to the results from the conventional and alternative strategies for optimizing the radial orbitals. (b) Same as the first panel, but for transitions between the 2 s 7 p 1 P 1 state and successive Rydberg 2 s n s 1 S 0 states in C III.
Atoms 07 00106 g002
Figure 3. (a) The transition rates A in the Babushkin (circles) and the Coulomb (triangles) gauges for the 2 p 2 P 1 / 2 8 s 2 S 1 / 2 transition in C IV, as a function of the increasing number of correlation layers. The transition rates computed in the conventional and the alternative manner are respectively shown in black and magenta. (b) Same as the first panel, but for the 2 s 2 1 S 0 2 s 7 p 1 P 1 transition in C III.
Figure 3. (a) The transition rates A in the Babushkin (circles) and the Coulomb (triangles) gauges for the 2 p 2 P 1 / 2 8 s 2 S 1 / 2 transition in C IV, as a function of the increasing number of correlation layers. The transition rates computed in the conventional and the alternative manner are respectively shown in black and magenta. (b) Same as the first panel, but for the 2 s 2 1 S 0 2 s 7 p 1 P 1 transition in C III.
Atoms 07 00106 g003
Figure 4. (a) The transition rates A in the Babushkin and the Coulomb gauges for the 2 p 2 P 1 / 2 8 s 2 S 1 / 2 transition in C IV, as a function of the square root of the upper integration bound R in the radial transition integrals (11) and (12) involving correlation orbitals. The radial transition integrals involving spectroscopic orbitals extend to their full values, so that R = 0 corresponds to transition rates computed from wave functions exclusively built from spectroscopic orbitals. The wave functions are produced by the alternative computational strategy. (b) The spectroscopic 2 p and 8 s radial orbitals in C IV as a function of r . The two orbitals occupy different regions in space and their overlap is minor. The 8 s orbital extends far out from the atomic core.
Figure 4. (a) The transition rates A in the Babushkin and the Coulomb gauges for the 2 p 2 P 1 / 2 8 s 2 S 1 / 2 transition in C IV, as a function of the square root of the upper integration bound R in the radial transition integrals (11) and (12) involving correlation orbitals. The radial transition integrals involving spectroscopic orbitals extend to their full values, so that R = 0 corresponds to transition rates computed from wave functions exclusively built from spectroscopic orbitals. The wave functions are produced by the alternative computational strategy. (b) The spectroscopic 2 p and 8 s radial orbitals in C IV as a function of r . The two orbitals occupy different regions in space and their overlap is minor. The 8 s orbital extends far out from the atomic core.
Atoms 07 00106 g004
Figure 5. (a) Same as Figure 4a, but for the 7 p 2 P 1 / 2 8 s 2 S 1 / 2 transition in C IV. (b) The spectroscopic 7 p and 8 s radial orbitals in C IV as functions of r . Both orbitals occupy nearly the same regions in space, overlapping to a great extent.
Figure 5. (a) Same as Figure 4a, but for the 7 p 2 P 1 / 2 8 s 2 S 1 / 2 transition in C IV. (b) The spectroscopic 7 p and 8 s radial orbitals in C IV as functions of r . Both orbitals occupy nearly the same regions in space, overlapping to a great extent.
Atoms 07 00106 g005
Table 1. The mean radii r (a.u.) of the spectroscopic and correlation orbitals that belong to the s and p symmetries in C IV. The correlation orbitals result from two different optimization schemes, the conventional and the alternative, and they occupy different regions in space.
Table 1. The mean radii r (a.u.) of the spectroscopic and correlation orbitals that belong to the s and p symmetries in C IV. The correlation orbitals result from two different optimization schemes, the conventional and the alternative, and they occupy different regions in space.
SpectroscopicCorrelationSpectroscopicCorrelation
ConventionalAlternative ConventionalAlternative
1 s 0.27 9 s 0.511.12 9 p 0.440.87
2 s 1.31 10 s 0.430.92 2 p 1.28 10 p 0.411.03
3 s 3.00 11 s 0.420.84 3 p 2.95 11 p 0.401.00
4 s 5.55 12 s 0.460.87 4 p 5.64 12 p 0.451.18
5 s 8.81 13 s 0.560.87 5 p 8.99 13 p 0.482.59
6 s 12.82 14 s 0.401.26 6 p 13.10 14 p 0.775.94
7 s 17.58 7 p 17.94
8 s 23.09 8 p 23.54
Table 2. Same as Table 1, but for radial orbitals in C III. The correlation orbitals 8 s and 8 p , which are introduced to account for the L S -term dependencies, are the same in both optimization schemes and fairly diffuse in comparison with the rest of the correlation orbitals.
Table 2. Same as Table 1, but for radial orbitals in C III. The correlation orbitals 8 s and 8 p , which are introduced to account for the L S -term dependencies, are the same in both optimization schemes and fairly diffuse in comparison with the rest of the correlation orbitals.
SpectroscopicCorrelationSpectroscopicCorrelation
ConventionalAlternative ConventionalAlternative
1 s 0.26 9 s 1.054.86 9 p 1.003.56
2 s 1.28 10 s 1.483.67 2 p 1.23 10 p 1.243.37
3 s 3.57 11 s 1.883.25 3 p 3.74 11 p 1.563.37
4 s 6.63 12 s 1.878.40 4 p 7.04 12 p 1.529.04
5 s 10.80 5 p 11.37
6 s 15.95 6 p 16.71
7 s 22.10 7 p 23.04
term corr. term corr.
8 s 8.22 8 p 5.55

Share and Cite

MDPI and ACS Style

Papoulia, A.; Ekman, J.; Gaigalas, G.; Godefroid, M.; Gustafsson, S.; Hartman, H.; Li, W.; Radžiūtė, L.; Rynkun, P.; Schiffmann, S.; et al. Coulomb (Velocity) Gauge Recommended in Multiconfiguration Calculations of Transition Data Involving Rydberg Series. Atoms 2019, 7, 106. https://doi.org/10.3390/atoms7040106

AMA Style

Papoulia A, Ekman J, Gaigalas G, Godefroid M, Gustafsson S, Hartman H, Li W, Radžiūtė L, Rynkun P, Schiffmann S, et al. Coulomb (Velocity) Gauge Recommended in Multiconfiguration Calculations of Transition Data Involving Rydberg Series. Atoms. 2019; 7(4):106. https://doi.org/10.3390/atoms7040106

Chicago/Turabian Style

Papoulia, Asimina, Jörgen Ekman, Gediminas Gaigalas, Michel Godefroid, Stefan Gustafsson, Henrik Hartman, Wenxian Li, Laima Radžiūtė, Pavel Rynkun, Sacha Schiffmann, and et al. 2019. "Coulomb (Velocity) Gauge Recommended in Multiconfiguration Calculations of Transition Data Involving Rydberg Series" Atoms 7, no. 4: 106. https://doi.org/10.3390/atoms7040106

APA Style

Papoulia, A., Ekman, J., Gaigalas, G., Godefroid, M., Gustafsson, S., Hartman, H., Li, W., Radžiūtė, L., Rynkun, P., Schiffmann, S., Wang, K., & Jönsson, P. (2019). Coulomb (Velocity) Gauge Recommended in Multiconfiguration Calculations of Transition Data Involving Rydberg Series. Atoms, 7(4), 106. https://doi.org/10.3390/atoms7040106

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