The Second Workshop on Lineshape Code Comparison: Isolated Lines

: In this work, we brieﬂy summarize the theoretical aspects of isolated line broadening. We present and discuss test run comparisons from different participating lineshape codes for the 2s-2p transition for LiI, B III and NV.

In a plasma, free electrons and ions interact and perturb atomic states, resulting in line broadening and shift [1,2].Each free plasma particle coming appreciably close to the emitter will contribute to the broadening and, in principle, also, the shift.This contribution will depend not only on the proximity of the emitter and perturber, but also on the duration of the interaction, i.e., the time interval during which the interaction is appreciable.Isolated lines are lines where for the upper and lower levels, the closest perturbing level is energetically much further away than the inverse collision duration.For instance, for the 2s-2p lines considered, the 2s-level states are perturbed by the 2p states in the sense that the important broadening process is: |2s⟩ → (via a collision) |2p⟩ → (via a collision) |2s⟩ (1) In general, we consider the broadening process: For isolated lines, the relevant energy spacings, ω αα ′ , between the collisionally connected states, α and α ′ , is of paramount importance.These energy spacings effectively reduce the interaction time and the effective impact parameters.For example, the relevant plasma-related quantity is, for perturbative collisions and a long-range dipole interaction [3]: where Here, E(t) is the electric field at the emitter, due to the perturbing particle with impact parameter ρ and velocity (at an infinite distance from the emitter) v.We note the decreasing width with increasing temperature, e.g., lnT √ T [4]; the width is dominated by a T −1/2 decay [5] for straight line trajectories.The imaginary exponential has the following important effects compared to the collisionally degenerate (hydrogenic, ω αα ′ = 0) case: 1.Only a part of the collision duration is effective, i.e., the effective collision duration is shortened, not the entire collision.2. Since the impact parameter affects the collision duration (∝ ρ/v), the effective impact parameters are also smaller.3. The effective velocities increase (slow collisions are adiabatic).4. As a corollary to all the above, broadening is decreased.
For strong collisions, we also have (canceling) higher order terms in the Dyson expansion.However, it is still true that the collision duration is effectively shortened by the imaginary exponentials.
In all calculations, only the 2s and 2p levels were considered.The only exception was the fully quantum-mechanical code [6], which used 13 nonrelativistic configurations 1s 2 nl(n = 2, 3, 4, 5, l ≤ 4) for B III and NV.

Electron vs. Ion Broadening
For hydrogen, given a velocity and impact parameter, a straight line trajectory gives the same contributions for any single electron and singly-charged ion.
For a hydrogen-like emitter, electrons are attracted; hence, they come closer, but also are accelerated, hence spending less time in the vicinity of the emitter.Ions are repelled, hence staying further away, but also slowing down; hence, they stay longer in the vicinity of the emitter.The action is ∫ V (t)dt, with V (t) the emitter-perturber interaction; hence, the net effect is that the same velocity electron and singly charged ion contribute the same to broadening [3].Of course, ion perturbers win overall, due to their velocity distribution favoring smaller velocities, both for hydrogen and hydrogen-like emitters.
For isolated lines, the action is ∫ V (t)e ıω αα ′ t dt.Therefore, only times of the order 1/ω αα ′ contribute.For electrons, this is not too bad: electrons rely on coming close, and since they are accelerated, their effective collision duration is small anyway.For ions, however, the consequences are far greater: ions were kept away and relied on a longer duration interaction to match the electrons.Now that a longer duration is negated by the shrinking of the effective collision duration to ∝ 1/ω αα ′ , the ions are ineffective [3].
The bottom line is that, for isolated lines, ions normally cannot compete with electrons, with two exceptions: (A) when ω αα ′ is small compared to kT, in which case, the collision duration was very short anyway; (B) ions can compete in non-dipole interactions.In contrast to the dipole case, we can have the channel |α⟩ As a result, isolated lines allow a test of electron broadening alone.For the calculations submitted to the workshop and considered in the present paper, ion broadening was explicitly neglected for both widths and shifts.The electron densities used were n e = 10 17 e/cm 3 for LiI, n e = 10 18 e/cm 3 for B III and n e = 10 19 e/cm 3 for NV.However, the results are presented, normalized to a density of 10 17 e/cm 3 , that is, the B III results were divided by 10 and the NV results by 100.

Penetrating Collisions
These are collisions where the perturbing particle penetrates the emitter wavefunction extent for the states involved.As already discussed above, since smaller impact parameters gain in importance, penetration is more important for isolated lines: For instance, for neutrals, the collision duration is ∝ ρ/v, so that the impact parameters, ρ, and velocities, v, that contribute significantly are such that ρω αα ′ ≤ v, e.g., small ρ.Hence, close collisions, which typically dominate for isolated lines, are significantly affected by penetration [7], which softens these collisions.This effect can be quite important and typically produces a factor of a two-width decrease [8].
Thus, if penetration is accounted for, perturbation theory is typically valid.This is convenient, since by softening the interaction, penetration makes those important close collisions perturbative; hence, much easier to solve.On the downside, there is no one-size fits all formula, i.e., instead of universal broadening functions, now we have per species and line broadening functions, although universal formulas that account for penetration have been given for hydrogen [10] and are possible for H-like ions.Furthermore, wavefunction information is needed, and oscillator strengths are no longer adequate [8].

Quantal Calculations
In contrast to semiclassical calculations, which assume the perturbing particles to be classical particles moving in well-defined trajectories, quantal calculations model the perturbing particles in terms of their wavefunctions and consider quantal effects, such as temporary electron capture and resonances.Only one fully quantal calculation participated in the workshop (DARC).
Fully quantal calculations have been performed by a number of different methods, e.g., Coulomb-Born (CB), distorted waves (DW), up to close coupling (CC) (R-matrix [11] and convergent close coupling (CCC) [12]).CB agreed very well with CC, which is understandable, due to the weakness of the interaction.Comparisons of quantal and semiclassical calculations have been done and summarized [8,13].Good agreement is found if the penetration is taken into account in semiclassical calculations.What is puzzling is a recent DW calculation by Elabidi et al. [14] citing the agreement with the experiments; these results are a factor of two larger than the aforementioned CC calculations.The present K-matrix CC calculations of DARC also agree with experiments and show the factor of two disagreements with the previous CC calculations.It is not clear why the DW [14] and DARC show such a large difference with R-matrix and CCC.

Work to be Done on Quantal Calculations
This is an open issue, and work is needed to pinpoint the origin of these differences.Comparing partial wave results and/or cross-sections may be helpful.We have not resolved the discrepancy with experiments.A code comparison may be needed to assess the relative importance of various factors in quantal calculations.

Work to be Done with Semiclassical Calculations
If penetration and non-perturbative aspects are accounted for, the remaining differences between the quantal and semiclassical approach are due to [8,13]: a.The demarcation line between what may and what may not be treated semiclassically.b.Back-reaction issues, i.e., in classical terminology, deviations from the (straight line or hyperbolic) trajectory, due to the energy transfer between the perturber and emitter; e.g., for collisional excitation/deexcitation, in our case, between the 2s and 2p levels, given that these energy transfers are comparable to the thermal velocities.
Answers to these questions may be given by Feynman path calculations.Such calculations have started appearing, but only for the simplest cases [15].

(Dynamic) Shifts
Shifts are caused by the finite collision duration.The same considerations apply as for widths: Ion (dipole) shifts are typically negligible; penetration is also an important factor.In contrast to widths, where one is looking for the guaranteed nonnegative deviation of Re(I − S a S * b ) of a quantity, possibly small, from unity, with shifts, one is looking for the deviation of a possibly small quantity, Im(I − S a S * b ) from zero.Furthermore, canceling opposite sign contributions arise [16], and the end result is that a small error in the calculation or deviation from thermal distribution functions can change the shift sign.This study was repeated [17] and confirmed the results for ∆n = 0, but found consistent red shifts for ∆n ̸ = 0.The recommendation of [17] is that ∆n ̸ = 0 shifts are reliable, while ∆n = 0 are not.

The Codes
Three types of codes participated in the workshop: DARC [6], a fully quantal (K-matrix) close-coupling code that, in principle, should be the benchmark, but uses a different physical model from the rest, full simulation codes (SimU [18] and Simulation [19]), which, in principle, may also account for perturber-perturber interactions, though for the cases presented in this comparison, perturbing electrons were modeled as non-interacting quasiparticles interacting with the emitter via Debye-shielded fields and semiclassical impact codes (SCP [20][21][22][23], a perturbative impact code and Starcode [8]).Starcode is a non-perturbative impact code, which solves the Schrödinger equation if it considers that the perturbative calculation has too large calculation uncertainties (and with the current tolerance, it invariably does).Although Starcode is a single code, it has a number of options, and in the present paper, we refer to two Starcode versions: one version, referred to as Starcode, which accounts for penetrating collisions, and one, called Starcode-NP, which does not.In practice, this is controlled by a single option in the input file.
The problem specification was that only dipole electron broadening should be considered.Except for DARC, which includes non-dipole interactions, all other codes have adhered, either by default or by selection of the appropriate options, to this specification.However, judging from Starcode calculations with monopoles and quadrupoles, the inclusion of monopole and quadrupole interactions makes too small a difference to explain the differences with DARC.Previous quantal close-coupling calculations have also produced widths that are a factor of two smaller than the ones quoted by DARC, as discussed above.
In addition, the problem specifications were to neglect fine structure, but only the simulation codes adhered to this specification.Among the codes, only SimU uses a Debye-shielded electron-emitter interaction, while all others use either a pure Coulomb interaction (DARC, Simulation) or a cutoff length of the order of the Debye length.In view of the remarks on short impact parameters being most important, these differences are not expected to be significant.For example, for SimU, the difference in the calculations using pure Coulomb or Debye-shielded fields is of the order of 0.1%.Small differences when using Coulomb or Debye-shielded fields have also been confirmed in tests using Simulation.
Atomic data are discussed in the Appendix.

Strong Collisions
The term "strong collisions" has been used in the literature to denote either "collisions that may not be treated by perturbation theory" or "collisions that may not be treated within the model used", e.g., collisions for which the semiclassical approximation is not valid or collisions for which the normally employed long-range approximation is not valid.Note that the two definitions may be in conflict.For example, penetrating collisions may be non-perturbative within the long-range approximation and perturbative if the long-range approximation is not used.Impact codes SCP and Starcode use the second definition and estimate the strong collision contribution based on unitarity.Starcode, which is non-perturbative, actually computes the non-perturbative contribution exactly within the impact approximation, but excludes from the computation the phase space for which its model is physically not valid.Simulation and SimU make no distinction between strong and weak collisions, as they are treated in exactly the same way.The simulation codes compute strong collisions (even if the model used is physically invalid due to quantal effects and/or the long-range dipole approximation), as well as weak ones.However, we consider simulation codes to employ the first definition when we wish to consider the differences between simulation codes and perturbative results.Hence, in principle, differences between the simulation codes and Starcode-NP should be due to either non-impact (overlapping) collisions or the phase space part that is excluded from Starcode and included in the simulation.For DARC, this distinction is rather meaningless.
Generally speaking, good agreement is expected between all codes for the non-strong (in either meaning) contributions, while things are not as clear for strong collisions in either meaning.

Calculational Uncertainties
With respect to calculational uncertainties, the codes use quite different ways of estimating them.
• DARC.Calculation uncertainties were estimated by increasing the orbital angular momentum of the colliding electron up to 25.They turned out to be ≪5% [6].• SCP.For simple spectra, an error bar of 20% was adopted, based on the average agreement with the experimental results for He I lines [24] and Si-like ions [25], namely 17% for widths and 20% for shifts.• SimU.This is a full simulation code that computes the line profile directly instead of first computing the autocorrelation function.It is a semiclassical, long-range dipole code that handles non-perturbative aspects, as well as non-impact effects rigorously.SimU calculates | ⃗ D(ω)| 2 , where ⃗ D(t) is the dipole operator and ⃗ D(ω) is its Fourier transform.The procedure is repeated N times (N ≫ 1) and averaged, which corresponds to an averaging over a statistically representative ensemble of radiators.The uncertainty due to a finite statistical sampling is ∝1/ √ N , and the proportionality coefficient is inferred by observing convergence over the course of simulations.For details, see [26].Another important factor affecting the calculational accuracy is the finite step of the frequency grid, equal to t −1 max , where t max is the time over which ⃗ D(t) is calculated in each run prior to averaging.Contrary to the first (statistical) factor, the latter always results in an overestimate of the line width by ∼t −1  max .For the shift values, this factor contributes ±∼t −1 max /2.The total uncertainties were within 10% for the cases considered.
• Simulation.The spectra is obtained through a Fourier transform of the autocorrelation function, which has been obtained by computer simulation and averaged over a large enough number of samples.The error is estimated using several sets of configurations to average the autocorrelation function and comparing the results obtained in each case.These fluctuations are about 3% using sets of 10,000 runs.• Starcode.Starcode has two sources of errors: numerical errors, e.g., from integrations or solutions of the Schrödinger equation, which are controlled through appropriate tolerances and are generally negligible, and errors arising from the model used, which are the ones quoted here.Starcode recognizes that there is a part of the phase space for which the model used (e.g., semiclassical trajectories and, in the version without penetration, trajectories that penetrate the wavefunction extent) is not reliable.However, the contribution of that part of the phase space may be bound by unitarity.Thus, one first defines v 0 = h mρmax .This is a velocity below which all collisions have a de Broglie wavelength larger than the screening length, ρ max , and, hence, are definitely non-classical.The corresponding phase space makes a contribution: = πnα with n the electron density, f (v) the velocity distribution, S a and S b the S-matrices for the upper and lower levels, respectively, w 0 = mv 2 0 2kT and α is an estimate for the average which is between zero and two, due to unitarity.In practice, this contribution is negligible in our examples.For v > v 0 , the semiclassical approach is not valid for ρ < bh mv , where b is a parameter, taken to be one in the present calculations, which specifies how much larger than the de Broglie wavelength an impact parameter must be in order to be treatable semiclassically.Hence, if penetration is allowed, this non-semiclassical phase space provides the dominant contribution: If penetration is not allowed, then the upper limit of the ρ integral is max(R, bh mv ) with R the relevant wavefunction extent, e.g., the maximum spatial extent of the 2s and 2p levels in our case.Hence, E 2 in the case with penetration is replaced by E 2N P + E 3 with: with v 1 = h mR and w 1 = 2kT and This ensures that, provided one trusts the choices of b and R, only truly semiclassical paths are computed, and a rigorous, unitarity-based error bar may be returned for the part of the phase space that is not treatable semiclassically.Of course, other sources of error, e.g., the back reaction, are not accounted for.
For the cases treated here with penetration not allowed, the wavefunction extent was dominant for all temperatures in the LiI calculations, for T = 15 and 50 eV for B III and T = 50 eV for NV.The de Broglie cutoff was dominant for the remaining cases.Hence, if the de Broglie cutoff is dominant, the Starcode and Starcode-NP calculations should have the same error bar, and Starcode-NP should have a larger width, due to the stronger interaction.If R is the dominant cutoff, then the error bar of Starcode-NP should be significantly larger and the computed width of Starcode might even exceed (with the error bar not accounted for) Starcode-NP, due to the larger phase space that is contributing.Of course, the situation will be reversed when the calculation uncertainties are accounted for.
Width results are quoted as (the width from the part of phase space that is computed reliably + 1 2 error bar) ± 1 2 error bar.Shifts are quoted as (the shift from the part of phase space that is computed reliably) ± error bar.All Starcode and Starcode-NP strong collision estimates assume a constant Re(I − S a S * b ) = 1 for that part of the phase space.

Results: Widths
Widths are shown in Table

General Remarks
1.In general, hyperbolic path calculations produced broader lines than straight line paths.This is expected, because due to the energy spacings (imaginary exponentials), only times closest to the time of the closest approach (t = 0) are important, but for these times, the electron is closer to the emitter if it is moving in a hyperbolic rather than a straight line path.Hence, the interaction is stronger, which results in larger widths.2. Non-perturbative collisions (impact or not) are treated correctly in simulations and by strong collision or error estimates in perturbative impact codes.This is not an issue for non-perturbative impact codes.Non-perturbative collisions are significant if penetration is not accounted for.3. Furthermore, in all cases, Starcode (with penetration) produces the lowest widths.This, however, is understood and expected, as Starcode with the option in question uses a different physical model that accounts for penetration, which, in turn, softens the interaction and results in smaller widths.4. DARC consistently produces large widths, in many cases, more than twice the results of Starcode, which also accounts for penetrating collisions.This is not understood and at variance with previous fully quantum mechanical close-coupling calculations, as discussed above.5.When investigating the differences between impact codes and simulations, particularly at low temperatures, one should keep in mind that impact codes treat all particles as impacting, while with simulation codes, the slow particles that are not impacting are treated correctly, i.e., differences could arise from non-impact effects for some particles.This would, in principle, mean lower widths for the simulations and should be less of an issue, at least for non-perturbative calculations at higher T. If that is correct, agreement should improve with non-perturbative impact calculations at higher T. As shown below, this is indeed the case for B III and NV, while the opposite is true for LiI.Note, however, that quantitatively non-impact effects are not expected to be an issue, as discussed below.6.Another issue is the temperature dependence, where we find some codes and cases with an increasing width with increasing temperature or a width that varies very little with temperature.Intuitively, one might expect that higher temperatures mean weaker collisions and, hence, smaller widths.We know, however, that an increase in temperature sometimes actually results in larger widths, e.g., for initially strong collisions, such as non-impact ions [27].Summing up, a non-decreasing temperature behavior, as observed in the simulation codes, may be the result of either non-impact effects or a significant or dominant strong collision term.The first possibility seems unlikely from simple estimates, discussed below, while if non-impact effects were negligible and the second were the case, then, since Starcode-NP is also a non-perturbative impact code, the differences from Starcode-NP would be due to the non-semiclassical phase space, as well as the emitter relevant wavefunction extent.This would have to be significant enough to change the widths vs. the temperature behavior.To check this conjecture, we list in Table 2 the Starcode-NP width contributions from the phase space part where Starcode-NP is presumably valid (labeled "non-strong" in Table 2, even though it includes non-perturbative collisions), as well as the "strong" collision contribution from the non-semiclassical and wavefunction extent phase space, where Starcode-NP (or any other workshop code, except DARC) is not valid.The actual width quoted is as described above, e.g., 0.83 + 0.23/2 ± 0.23/2 for the lowest LiI temperature.Therefore, we conclude that the differences between Starcode-NP and the simulation codes, as well as the different temperature trends exhibited are within the calculation uncertainties.

LiI Results
Figure 1 displays the results from five codes.Interestingly, only Starcode, in both versions, displays a decreasing width with temperature behavior (Within the semiclassical model, threshold issues do not arise, as we have virtual excitations and deexcitations, and a higher T means weaker interactions, due to their shorter duration.)All three other calculations produce an increase in the width as T rises from five to 15 eV.Simulation varies very little with temperature.With SimU, that trend continues out to T = 50 eV, whereas for SCP and Simulation, it is reversed.For SCP, this may be attributed to the symmetrization, which respects the quantum properties of the S-matrix rather than the method of computing the strong collision term, which is less than half the width already at 5 eV and decreases with T. For simulations, a possible explanation would be a non-impact behavior for electrons at low T.However, simple estimates do not support this conjecture.For example, at T = 1 eV, n = 10 17 e/cm 3 and L − α, electrons are well impacted [28], and it is hard to see why they would be non-impacted at 5 eV for Li with a longer lifetime (inverse HWHM), due to the smaller matrix elements and imaginary exponentials.Hence, the most likely origin is [27], the increasing width vs. temperature for strong collisions.As already discussed, however, within the calculational uncertainties, Starcode-NP could also be made to show an increasing width with temperature.

B III Results
Figure 2 displays the widths from five codes.General trends are confirmed, with DARC giving the largest widths.SCP and Starcode-NP exhibit fairly good agreement, but the disagreement of SimU with Starcode-NP and the agreement of SimU with Starcode at low T (which is not expected, as SimU does not account for penetrating collisions, which make the difference between Starcode and Starcode-NP) and Starcode-NP at high T (which is expected) is puzzling.As shown in Table 2, the non-strong Starcode-NP contribution alone is significantly larger than SimU's total width, and the strong contribution will only increase the width of Starcode-NP and, hence, the difference.Figure 1.A comparison of LiI 2s-2p transition widths as computed by SCP (solid black), SimU (dotted red), Simulation (dashed green) and Starcode with (dash-dotted blue) and without (dash-double dotted orange) penetration accounted for.A comparison of hyperbolic trajectory B III 2s-2p transition widths as computed by SCP (solid black), SimU (dotted red), DARC (dashed green) and Starcode with (dash-dotted blue) and without (dot-double dash orange) penetration accounted for.Figure 3 compares straight and hyperbolic trajectory results where applicable.Overall, the comparison demonstrates the difference, particularly for SCP, that arises from including vs. neglecting the attraction by the emitter.One would have expected a closer agreement between SimU and Simulation, which use the same physical model.For SimU and Simulation, we see a width increase with increasing T, as already discussed.This increase is reversed for Simulation, but not for SimU.At high T, the calculations converge, as expected.

NV Results
Figure 4 displays the widths from five codes.SimU and Starcode-NP agree quite well, which is expected, at least for high temperatures, where non-impact effects should play no role.SCP and DARC also agree quite well with each other; however, their disagreement with other codes at the highest T value is not understood, since close coupling effects are expected to be negligible at this high T.All codes show the expected decrease for the widths as a function of T, as collisions weaken with increasing emitter charge Z [25].
Figure 5 compares straight and hyperbolic trajectory results, where applicable.Overall, the comparison demonstrates the significant difference that arises from including vs. neglecting the attraction by the emitter.Due to the higher Z, it is expected that these differences would be more pronounced than the differences found for B III.The close agreement of SimU-straight and Simulation-straight is expected, as both codes use the same physical model.

Results: Shifts
As discussed, shift calculations within the stated case definitions are probably inaccurate, but the comparisons may well illustrate trends and code differences.The results are summarized in Table 3. Shifts were normalized to a density of 10 17 e/cm 3 .In general, we expect that if shifts are perturbative, they should also decrease with increasing T, due to the weakening of the interaction.For the same reason, we expect a convergence of the straight line and hyperbolic trajectory results as T → ∞.No such statement of decreasing shifts with increasing T can be made with certainty in the non-perturbative case.2. Starcode, in particular (with or without penetration), often shows a very weak increase of shift with T, even in the version with penetration accounted for, which results in weak, perturbative collisions.The reason is that Starcode explicitly excludes the non-semiclassical phase space, e.g., impact parameters shorter than the de Broglie wavelength.(In addition, Starcode, without accounting for penetration, also excludes impact parameters inside the wavefunction extent, because, again, these may not be treated within a long-range approximation).When T decreases, this phase space, which may not be computed within the code's semiclassical framework, increases.Hence, the phase space that gives a shift that may presumably be reliably computed shrinks, while the phase space that may not be treated grows.For both shifts and widths, Starcode binds the contribution of the non-semiclassical phase space by unitarity and returns it as an error bar.Therefore, while the quoted shift may be fairly insensitive to T, the error bar is not.These results are to be interpreted as an uncertainty in shift calculations inherent in semiclassical calculations, i.e., the part of the phase space that may not be treated semiclassically can make a substantial shift contribution rather than a temperature-insensitive shift.3. Straight line SimU and Simulation agreement is excellent for the charged cases (B III and NV), as expected and, surprisingly, much poorer (though still good) for LiI. 4. SCP consistently gives significantly smaller shifts for ion lines (B III and NV).Furthermore, straight-line SCP calculations give significantly larger shifts than hyperbolic path SCP calculations.This is counterintuitive, since as discussed, hyperbolic paths result in a stronger effective interaction close to the emitter and, hence, action.Furthermore, the straight-line SCP shifts are monotonically decreasing as a function of T, in contrast to the hyperbolic trajectory results.This behavior follows the semiclassical b-functions [2], which result in lower shifts for hyperbolic compared to straight-line paths and, in addition, a decreasing b-function with increasing ion charge.

LiI
Figure 6 plots the shifts computed by five codes vs. electron temperature T. Except for Starcode (explained above), in both versions, the shift is decreasing with T, which is expected in view of the weakening of the interaction.The differences between the other three codes are generally better than 20%, although it is not clear why SimU and Simulation, which essentially use the same model and, hence, would, in principle, be expected to give the same results, differ by that amount.shows a small increase between five and 15 eV.SCP is significantly lower for small T, which might be attributable to the use of perturbative impact theory by SCP, as it reflects the behavior of the semiclassical b-function.This is consistent with the convergence seen for higher T.

NV
Figure 9 plots the shifts computed by five codes vs. electron temperature T. DARC produces significantly larger shifts at low T, with all codes converging for large T. The large difference with Starcode seems to indicate that higher order partial waves (except the s-wave, which is viewed as an error bar in Starcode) differ substantially from the semiclassical results.SCP is lowest and exhibits insensitivity to T. SimU and Starcode without penetration also agree quite well, except at the lowest T. The lowest temperature point discrepancy between SimU and Starcode-NP could be due to the fact that Starcode does not compute the short impact parameter (ρ < h mv ) contribution, which is not semiclassical, but adds an estimate to the shift, as discussed above.The calculations agree within their respective calculation uncertainties, although one might have expected a better agreement.However, attributing the difference to the phase space not included in the Starcode-NP calculation (but added as an error bar) does not seem to hold: For example, as already discussed, Starcode-NP calculations treat collisions with impact parameters larger than the de Broglie wavelength and use a unitarity-based estimate for smaller impact parameters.If instead, Starcode-NP treats all collisions with impact parameters larger than 0.1 times the de Broglie wavelength and only uses unitarity estimates for impact parameters smaller than 0.1 times the de Broglie wavelength, the shift actually decreases and is in the range of 1 = 1.5 cm −1 at the lowest temperature, compared to a value of 3 cm −1 for SimU.The three straight line results are in fairly good agreement, with SimU and Simulation being in excellent agreement.The hyperbolic trajectory and straight line SimU results are interesting in that the hyperbolic trajectory shift is larger at the lowest T (expected, as discussed in view of the stronger effective interaction for hyperbolic rather than straight line trajectories) and slightly smaller for the other two T values; however, that difference is within the calculation uncertainties, and the two calculations converge, as expected, for high T.
Figure 2.A comparison of hyperbolic trajectory B III 2s-2p transition widths as computed by SCP (solid black), SimU (dotted red), DARC (dashed green) and Starcode with (dash-dotted blue) and without (dot-double dash orange) penetration accounted for.

Figure 7
Figure7plots the shifts computed by five codes vs. electron temperature T. The general trend (shifts decreasing with T) is observed, with the exception of the non-penetrative Starcode, which gives identical (within errors) results for the 15 and 50 eV cases, which has been discussed above, and SCP, which

Figure 7 .
Figure 7.A comparison of B III 2s-2p transition shifts as computed by SCP (solid black), SimU (dotted red), DARC (dashed green) and Starcode with (dash-dotted blue) and without (dot-double dash orange) penetration accounted for.

Figure 8
Figure 8 plots the shifts computed by five codes vs. electron temperature T. The straight line SimU and Simulation agreement is very good, as discussed above.The SimU hyperbolic trajectory shifts are

Figure 10
Figure 10 plots the shifts computed by five codes vs. electron temperature T.

Table 1 .
1.All widths have been normalized to a density of 10 17 e/cm 3 .2s-2p Transition width (FWHM) comparisons.All widths are in cm −1 and normalized to a density of 10 17 e/cm 3 .

Table 2 .
2s-2p Transition width (FWHM) strong collision importance in Starcode-NP.All widths are in cm −1 and normalized to a density of 10 17 e/cm 3 .

Table 3 .
2s-2p Transition shift comparisons of all shifts are in cm −1 and normalized to a density of 10 17 e/cm 3 .