Attosecond Time Delay Trends across the Isoelectronic Noble Gas Sequence

: The analysis and measurement of Wigner time delays can provide detailed information about the electronic environment within and around atomic and molecular systems, with one the key differences being the lack of a long-range potential after a halogen ion undergoes photoionization. In this work, we use relativistic random-phase approximation to calculate the average Wigner delay from the highest occupied subshells of the atomic pairings (2p, 2s in Fluorine, Neon), (3p, 3s in Chlorine, Argon), (4p, 4s, 3d, in Bromine, Krypton), and (5p, 5s, 4d in Iodine, Xenon). The qualitative behaviors of the Wigner delays between the isoelectronic pairings were found to be similar in nature, with the only large differences occurring at photoelectron energies less than 20 eV and around Cooper minima. Interestingly, the relative shift in Wigner time delays between negatively charged halogens and noble gases decreases as atomic mass increases. All atomic pairings show large differences at low energies, with noble gas atoms showing large positive Wigner delays, while negatively charged halogen ions show negative delays. The implications for photoionization studies in halide-containing molecules is also discussed.


Introduction
The recent advancement of attosecond extreme ultraviolet infrared (XUV-IR) laser metrology over the past decade [1][2][3][4][5][6] has enabled access for observing ultrafast phenomena across a variety of atomic and molecular systems at the natural time scale of their electronic motion. One common experimental technique utilizes an XUV-IR pump-probe process [7,8], where an electron is first ionized through the absorption of an XUV photon and subsequently streaked by the few-cycle IR laser field, which imprints itself on the photoelectron's final energy and momentum. By varying the time delay between the XUV and IR pulses, it is possible to measure the photoionization time delay relative to a reference. Another common technique is reconstruction of attosecond beating by interference of two-photon transitions (RABBITT) [9,10], where the target is first ionized by an XUV attosecond pulse train of high-order harmonics, and the photoelectron can then either absorb or emit a secondary IR photon in the continuum. By adjusting the delay between the IR laser field and the high-order harmonics, it is possible to extract the time delay for a particular transition.
These techniques have been used to measure photoionization delays in atoms [11][12][13][14][15], molecules [16][17][18][19][20], and solid-state systems [21][22][23], thereby providing new information about their electronic structure. The total time delay τ, also known as the streaking delay, given by a pump-probe experiment is frequently separated into two components, τ = τ w + τ CLC . This convention of separating the total delay into two separate terms is also followed in traditional RABBITT experiments, with the only difference being that τ CLC is replaced by the continuum-continuum delay τ CC [10].The first contribution is the Wigner time delay [24,25] τ w , which describes the group delay of the ionized photoelectron

Theoretical Overview of the Relativistic Random-Phase Approximation (RRPA)
The same theoretical formalism was applied to negative halogen ions and neutral noble gas atoms in order to calculate the photoionization dipole transition matrix amplitudes and phases for an initial bound state. This overview follows the same outline as the previous RRPA photoionization studies of Kheifets and Deshmukh [50] and the original multi-channel paper of Johnson and Lin [49]. Atomic units ( = e = m e = 1) are used throughout, unless stated otherwise. For a time-dependent perturbation of the form v + e −iωt + v − e iωt , the probability amplitude for a transition from the ground state u i to an excited state w i± , stimulated by said perturbation, is given by For an electromagnetic interaction described in the Coulomb gauge, it is possible to rewrite the transition amplitude T as a function of the vector potential A, where the perturbations v ±\ are described by, v Therefore, a photon of frequency ω (or equivalently of wavevector k ) and polarization can be described by the vector potential in the Coulomb gauge as A =ˆ e −i k · r , which can then be expanded in terms of its multipole components, a (λ) JM . Note that an upper index of λ = 1, 0 corresponds to the electric and magnetic multipoles, respectively. In the single active electron approximation, the transition amplitude describe by Equation (2) can be reduced even further to T (λ) where J and M are the angular momentum quantum numbers describing the incoming photon. It is common to label the initial bound state (u i ) of the electron by the quantum numbers ljm and the final continuum state (w i ) by the numbers The spin of the electron is given by the spinor χ ν , where ν = ±1/2. Any final state can be uniquely described by the index − κ = ∓ − j + 1/2 , where − j = − l ± 1/2 is the total angular momentum of the outgoing electron in the continuum. The final state can also be written as a partial wave expansion, which is given explicitly in Equation (40) of [49]. Inserting this expansion into Equation (3) results in The ionized photoelectron's energy and momentum are represented by E and p, respectively, while Ω− κ − m (p) is described in terms of Clebsch-Gordan coefficients and spherical harmonics.
The six-term bracket to the right of the Ω− where R (λ) J − a, a is a radial integral listed in [49] and π − l , l, J − λ + 1 is simply the parity factor responsible for imposing the necessary selection rules for a given transition.
Despite its use for describing single-electron transitions, Equation (4) is also valid for any closed-shell atomic species, as the only change required in Equations (4) and (6) is that the single-electron reduced matrix element is modified slightly to include multi-electron RRPA effects in both the initial and final states. An explanation of this modification is given In this work, we restrict ourselves to electric dipole transitions (λ = 1 and J = 1) with linearly polarized light oriented along theẑ axis (M = 0). From these restrictions, Equation (4) simplifies to Here, we have followed the convention of ref. [51] by defining separate transition amplitudes for the spin-up and spin-down cases and by omitting the scaling factor of 2πi √ 3Ep ω. By choosing the linear polarization of the incoming photon to be purely alongẑ, it is possible to take advantage of the axial symmetry of the system. Therefore, we will introduce the shorthand Y lm ≡ Y lm (p) = Y lm (θ, φ = 0), where θ = 0 corresponds to a photoelectron emission parallel to the direction of the initial photon.
It useful to observe how Equation (10) directly reduces in the simple case of electric dipole transitions for the np 1/2 and np 3/2 states. The expressions for the nd 3/2 and nd 5/2 transition amplitudes can be found in ref. [50].
The reduced matrix elements D lj→ − l − j can be evaluated numerically for both its real and imaginary components, as doing so allows for the phase η and Wigner time delay to be calculated using the standard formulation of For a given subshell nl j , the angle-dependent time delay can be calculated as a weighted average over all possible transition amplitudes and spin states. For the purposes of this work, we will only consider the case of θ = 0, as it is commonly the most dominant direction of photoelectron emission. As with Equation (10), it is informative to see how Equation (14) simplifies into the simple weighted average of the spin-up and spin-down Wigner delays for the case of an np 1/2 state.
This process of averaging Wigner delays was performed for every subshell listed in Table 1 below. By taking the average Wigner delays τ np 1/2 and τ np 3/2 and weighting them by their respective differential cross-sections, the total average time delay τ np was also calculated. An analogous process was also used to compute τ ns and τ nd for every halogen and noble gas.
It should be noted that the accuracy of this averaging process depends not only on the values of the Wigner delays, but also with regard to the accuracy of the transition amplitudes themselves. Because the photoionization cross-section can be computed from the transition amplitudes T (1ν) 10 m nl j , it is possible to estimate the accuracy of the RRPA calculations by simply comparing the predicted cross-sections with those of experimental measurements. Reference [38] lists the predicted RRPA cross-sections for all of the noble gases being studied in this work, and found there to be a good agreement with experimental values. This result implies that the calculated RRPA transition amplitudes given by Equation (10) should also be quite accurate.

Results and Discussion
In this section, we present the calculated binding energies and average Wigner delays for each of the highest subshells of the halide ions and noble gas atoms. It should be noted that RRPA often produces autoionization resonances in the time delay spectra for any given noble gas. However, they are not the focus of this paper and, therefore, have been filtered out to prevent the obfuscation of more general time delay features. The locations of autoionization resonances in noble gases are well documented [38], but they generally occur at photon energies close to the binding energies of orbital subshells. In Argon, for example, the autoionization resonances produced by RRPA occur around photon energies of 35 eV, which corresponds to the ionization threshold of the 3s 1/2 subshell.

Dirac-Hartree-Fock (DHF) Orbital Subshell Ionization Calculations
RRPA requires the use of Dirac-Hartree-Fock orbitals in order to account for ab initio relativistic effects and to obtain accurate subshell ionization potential energies. It is important to note that DHF calculations are typically the most accurate for the highest occupied shells regardless of the atom being studied; however, they can reasonably predict the binding energy of lower-lying subshells as atomic mass increases. Table 1 confirms this trend, as the predicted value of the 4d orbitals in xenon are within~4 eV of experimental measurements. The binding energies of halide ions are not well known, yet it is possible to approximate their highest experimental binding energies with electron affinity measurements. The electron affinity values were found to closely match the calculated DHF energies, with the approximation being increasingly valid for the higher-mass ions of bromide and iodide. These trends and absolute energies were also found to agree with the calculated values of Saha et al. [38,41] and Lindroth and Dahlstrom [39]. The congruence between our calculated DHF binding energy for F − and the 2p 3/2 energy reported in [39], which utilized a non-relativistic HF theory with exchange, is of particular interest as it implies that despite not being necessary, relativistic effects do not negatively impact Wigner time delay calculations for lighter-mass atomic systems.  Table 1 also lists the number of coupled photoionization channels used for the subsequent calculations of the reduced matrix elements. For Neon and Fluorine, all nine possible channels (2p, 2s, 1s) were coupled. Argon and Chlorine used 14 channels (3p, 3s, 2p, 2s) and omitted the 1s channels. Krypton and Bromine used 20 channels (4p, 4s, 3d, 3p, 3s) and omitted the 2p, 2s, and 1s channels. Finally, Xenon and Iodine used 33 channels (5p, 5s, 4d, 4p, 4s, 3d, 3s) and omitted the core 2p, 2s, and 1s channels. The omitted channels could be ignored due to the fact that they are significantly farther away in energy from the other channels and do not have substantial impact at the photon energies of this study.

Individual Photoionization Channel Phases and Wigner Time Delays of Neon and F −
While a previous study compared the Wigner delays of individual ionization channels in Argon and Chlorine [41], to the best of our knowledge, no similar study has been performed for the lighter pair of Neon and Fluorine. The phase η of each channel was calculated directly by computing the argument of the reduced matrix element associated with each channel (i.e., via Equations (6) and (13)). By simply taking the energy derivative of the phase, the Wigner time delays could also be determined, as seen in Figure 1. Due to their low atomic mass and relatively small number of electrons, relativistic effects do not play a crucial role in time delays in Neon and Fluorine. This is reflected in the behavior of both the phase and Wigner time delays for any given channel, since the j and − j values of a particular l → l + 1 or l → l − 1 transition have a negligible impact (e.g., the 2p 1/2 → s 1/2 and 2p 3/2 → s 1/2 both produce the same η and Wigner delay). This fact also holds true for the 3p 1/2 and 3p 3/2 channels in both Chlorine and Argon [41], but begins to break down for the lower-lying 3s 1/2 → p 1/2 and 3s 1/2 → p 3/2 channels. While not the primary focus of this work, the study of time delays in individual channels is often useful for analyzing the results of the average orbital delays, since it is possible to determine where a given channel dominates in a particular energy region and instructive to see how the average delay results from Wigner delays of individual transitions. For the sake of brevity, the individual channel phases and time delays for Cl − , Ar, Br − , Kr, I − , and Xe are omitted.  (6) and (13)). By simply taking the energy derivative of the phase, the Wigner time delays could also be determined, as seen in Figure 1. Due to their low atomic mass and relatively small number of electrons, relativistic effects do not play a crucial role in time delays in Neon and Fluorine. This is reflected in the behavior of both the phase and Wigner time delays for any given channel, since the and ̅ values of a particular → + 1 → − 1 transition have a negligible impact (e.g., the 2 ⁄ → ⁄ and 2 ⁄ → ⁄ both produce the same and Wigner delay). This fact also holds true for the 3 ⁄ and 3 ⁄ channels in both Chlorine and Argon [41], but begins to break down for the lower-lying 3 ⁄ → ⁄ and 3 ⁄ → ⁄ channels. While not the primary focus of this work, the study of time delays in individual channels is often useful for analyzing the results of the average orbital delays, since it is possible to determine where a given channel dominates in a particular energy region and instructive to see how the average delay results from Wigner delays of individual transitions. For the sake of brevity, the individual channel phases and time delays for Cl , Ar, Br , Kr, I , and Xe are omitted.
By comparing Figures 1b and 2b, it is apparent that the average 2 delay in Ne is dominated by the 2 → channels across all photon energies, whereas the average 2 delay in F primarily corresponds to 2 → transitions at higher energies and 2 → channels below photon energies near 10 eV. This agrees with the fact that for typical photoionization in noble gases, the → + 1 photoionization channels are known to dominate regardless of energy, with the only exception being Cooper minima. However, in photodetachment, the → − 1 channels dominate near the threshold and the → + 1 channels only begin to dominate once photon energy increases [41,45]. A cursory comparison of the time delays between Neon and Fluorine shows a strong agreement for photon energies greater than 50 eV for any given 2 transition. The same can also be said of 2 → transitions at photon energies above 75 eV. There is, however, a sharp contrast between the time delays at lower energies. The individual labels are a descriptive shorthand to describe the identical behavior of multiple ionization channels. For instance, the label 2 → relates to the two ionization channels 2 ⁄ → ⁄ and 2 ⁄ → ⁄ , while the label 2 → corresponds to the following three channels: 2 ⁄ → ⁄ , 2 ⁄ → ⁄ , and 2 ⁄ → ⁄ . The behavior of 2 → is equivalent to the behavior of the 2 ⁄ → ⁄ and 2 ⁄ → ⁄ channels. In the case of Ne, a small autoionization resonance was removed near 48 eV for the 2 → and 2 → channels.

Comparison of Average Wigner Delays for Halogen Ions and Noble Gases
We observe in Figure 2a,b that Neon exhibits the familiar time delay behavior of having a large positive delay near the threshold energies of the 2 and 2 subshells, which The individual labels are a descriptive shorthand to describe the identical behavior of multiple ionization channels. For instance, the label 2p → s relates to the two ionization channels (2p 1/2 → s 1/2 and 2p 3/2 → s 1/2 ), while the label 2p → d corresponds to the following three channels: (2p 1/2 → d 3/2 , 2p 3/2 → d 3/2 , and 2p 3/2 → d 5/2 ). The behavior of 2s → p is equivalent to the behavior of the (2s 1/2 → p 1/2 and 2s 1/2 → p 3/2 ) channels. In the case of Ne, a small autoionization resonance was removed near 48 eV for the 2p → d and 2p → s channels.
By comparing Figures 1b and 2b, it is apparent that the average 2p delay in Ne is dominated by the 2p → d channels across all photon energies, whereas the average 2p delay in F − primarily corresponds to 2p → d transitions at higher energies and 2p → s channels below photon energies near 10 eV. This agrees with the fact that for typical photoionization in noble gases, the l → l + 1 photoionization channels are known to dominate regardless of energy, with the only exception being Cooper minima. However, in photodetachment, the l → l − 1 channels dominate near the threshold and the l → l + 1 channels only begin to dominate once photon energy increases [41,45]. A cursory comparison of the time delays between Neon and Fluorine shows a strong agreement for photon energies greater than 50 eV for any given 2p transition. The same can also be said of 2s → p transitions at photon energies above 75 eV. There is, however, a sharp contrast between the time delays at lower energies. eV. This minimum matches that of [41] although it is not shown due to the scale of the figure. The average 2 and 3 Wigner delays of F and Cl , respectively, originally displayed oscillations at photoelectron energies below 15 eV. These oscillations were determined to be caused by small variations in the average phase data and were subsequently removed by taking a best fit of the phase data and then computing the energy derivative of that best fit.
For the higher-mass systems of Bromine, Krypton, Iodine, and Xenon illustrated in Figure 3, the average Wigner time delays for each orbital were also found to diverge at low photoelectron energies, as Br and I both display large negative time delays (again on the order of −100 as) in the low energy region. This appears to indicate a universal time delay trend between negative charged halogens and noble gas atoms at photoelectron energies below 20 eV. While it is possible to explain the negative delay times as the result of a negative energy derivative of the phase, a physical explanation is less obvious. If one interprets a positive delay time as the retardation of the photoelectron wave with  [41] although it is not shown due to the scale of the figure. The average 2s and 3s Wigner delays of F − and Cl − , respectively, originally displayed oscillations at photoelectron energies below 15 eV. These oscillations were determined to be caused by small variations in the average phase data and were subsequently removed by taking a best fit of the phase data and then computing the energy derivative of that best fit.

Comparison of Average Wigner Delays for Halogen Ions and Noble Gases
We observe in Figure 2a,b that Neon exhibits the familiar time delay behavior of having a large positive delay near the threshold energies of the 2s and 2p subshells, which slowly vanishes as photon energy increases. Fluorine, by contrast, exhibits a strong negative delay on the order of −100 as near the 2p threshold. This difference can be explained by comparing the differences in the calculated phases for both Fluorine and Neon illustrated in Figure 1a, as each of Neon's ionization channels displays a dramatic increase near their respective threshold energies, while Fluorine's 2s → p and 2p → s channels tend to decrease more gradually over a longer energy range. For noble gas atoms, the Wigner time delay of any given orbital will trend towards positive infinity at energies near the threshold due to the drastic increase in the Coulomb phase, which is known to dominate [40,41], and the individual channels of Ne confirm this trend. The Coulomb phase for negatively charged atoms is essentially nonexistent due to the lack of a long-range Coulomb potential and, therefore, the Wigner delays corresponding to photodetachment do not possess the same behavior of trending towards infinity at low energies. In the case of photodetachment, however, it is still possible for a short-range potential to play a role at energies extremely close to the threshold. For example, in Figure 1a, the phase of the 2p → d ionization channels in Fluorine rapidly decreases in an energy region of ∼ 1 eV near threshold and then slowly increases with photon energy. This also explains the sharp increase in the average 2p Wigner delays in Fluorine at low photoelectron energies below ∼ 10 eV and small positive delay times in the higher photoelectron energy region.
Just as with the individual channel analysis in Fluorine and Neon, we find that the average 2s and 2p Wigner delays also diverge at lower photoelectron energies and converge towards zero in the large photoelectron kinetic energy limit. The physical picture underlying this vanishing time delay is quite clear, as a photoelectron with high kinetic energy will spend less time near the perturbative effects of the remaining ion and instead behave more similarly to a free electron wave. Conversely, at low kinetic energies, the photoelectron will spend more time near the remaining ion and be more susceptible to collective electron effects. We also observe the same general feature of diverging Wigner delays between the 3s and 3p states in Chlorine and Argon at low energy (see Figure 2c,d), with the only difference being the introduction of Cooper minima. A Cooper minimum occurs when the transition matrix element changes sign and undergoes a phase shift of π. Typically, this happens when the initial state radial wavefunction contains at least one node and overlaps with the continuum wavefunctions. This is the process responsible for the well-known 3p Cooper minima illustrated for Argon and Chlorine in Figure 2d. However, a different mechanism is responsible for the observed minimum in the average 3s Wigner delay of Argon and Chlorine. Instead of being the result of a radial node in the initial wavefunction, the behavior of the Cooper minimum in the 3s channel is caused by a π shift in the phases of the 3s → p channels, which occurs due to the result of significant interchannel coupling with the 3p ionization channels [38,56,57]. Without the effect of interchannel coupling, this minimum does not appear in the 3s ionization channels, which is why the 3s minimum can be deemed an "induced Cooper minimum" as it is still the result of a π shift but its origin differs from that of the 3p Cooper minimum. However, the location of this induced Cooper minimum in Chlorine appears to be shifted ∼ 10 eV higher than the induced Cooper minimum in Argon. A similar shift occurs for the Cooper minimum [58] in the 3p spectra, with Chlorine again displaying a higher shift of ∼ 6 eV. Because Figure 2 plots both time delays with respect to the kinetic energy of the ionized electron, any difference in the location of the two Cooper minima must be the result of properties not related to the threshold energies of the atomic systems. If the shift was solely caused by a difference in binding energies, the two Wigner delay spectra would overlap.
For the higher-mass systems of Bromine, Krypton, Iodine, and Xenon illustrated in Figure 3, the average Wigner time delays for each orbital were also found to diverge at low photoelectron energies, as Br − and I − both display large negative time delays (again on the order of −100 as) in the low energy region. This appears to indicate a universal time delay trend between negative charged halogens and noble gas atoms at photoelectron energies below 20 eV. While it is possible to explain the negative delay times as the result of a negative energy derivative of the phase, a physical explanation is less obvious. If one interprets a positive delay time as the retardation of the photoelectron wave with respect to a free electron wave of the same kinetic energy, then a negative Wigner delay can be interpreted as an acceleration in the outgoing electron wave packet. In the case of the induced Cooper minimum in the average 3s state for Argon at low energy, the "dip" in the delay time can be thought to be the result of induced oscillations in the outer 3p subshells that screen the 3s electrons [57]. This only can occur for systems where the ns and np states are strongly coupled, which explains the absence of a similar feature in the 2s time delay spectra of Neon due to its 2s and 2p interchannel coupling being much weaker. and states are strongly coupled, which explains the absence of a similar feature in the 2 time delay spectra of Neon due to its 2 and 2 interchannel coupling being much weaker.
The Cooper minima in the average 5 Wigner delay for Xenon have also been explained to be the result of interchannel couplings, with the first Cooper minimum near 35 eV being the result of couplings with the 5 ⁄ and 5 ⁄ states and the second minimum at ~120 eV being the result of couplings with the 4 ⁄ and 4 ⁄ states [59,60]. In the case of I , we also find the same 5 Cooper minima at energies approximately equal to those of Xenon. A similar equivalence was observed between the average 5 Wigner delays of Iodine and Xenon, as well as for the 4 , 4 , and 3 delays in Br and Kr. However, the minimum in the 4 time delay spectra for I was found to lie ~11 eV higher than for Xenon. Comparing the relative shifts between Cooper minima for the average Wigner delays of each halogen and noble gas pairing, we find a qualitative trend where the difference decreases as atomic mass increases. The origin of this trend is still not entirely clear, although it could be related to the relative difference in the electronegativities and polarizabilities of the halogen and noble gas atoms, as relative shifts have been observed in the photodetachment cross-sections for different isotopes of Chlorine [61]. A previous theoretical study [62] also reported shifts in the respective 3 and 3 cross-sections for Ar and Cl and other high-Z isoelectronic species that were positively charged (i.e., Sc , Mn ) and concluded that interchannel coupling and initial-state correlation effects can impact the location of Cooper minima, although they do not directly account for the difference in binding energies between Ar and Cl . Despite not focusing on the relative shift in photoionization cross-sections, our work regarding the 3 Wigner delays between Ar The Cooper minima in the average 5s Wigner delay for Xenon have also been explained to be the result of interchannel couplings, with the first Cooper minimum near 35 eV being the result of couplings with the 5p 1/2 and 5p 3/2 states and the second minimum at ∼ 120 eV being the result of couplings with the 4d 3/2 and 4d 5/2 states [59,60]. In the case of I − , we also find the same 5s Cooper minima at energies approximately equal to those of Xenon. A similar equivalence was observed between the average 5p Wigner delays of Iodine and Xenon, as well as for the 4s, 4p, and 3d delays in Br − and Kr. However, the minimum in the 4d time delay spectra for I − was found to lie ∼ 11 eV higher than for Xenon.
Comparing the relative shifts between Cooper minima for the average Wigner delays of each halogen and noble gas pairing, we find a qualitative trend where the difference decreases as atomic mass increases. The origin of this trend is still not entirely clear, although it could be related to the relative difference in the electronegativities and polarizabilities of the halogen and noble gas atoms, as relative shifts have been observed in the photodetachment cross-sections for different isotopes of Chlorine [61]. A previous theoretical study [62] also reported shifts in the respective 3s and 3p cross-sections for Ar and Cl − and other high-Z isoelectronic species that were positively charged (i.e., Sc 3+ , Mn 7+ ) and concluded that interchannel coupling and initial-state correlation effects can impact the location of Cooper minima, although they do not directly account for the difference in binding energies between Ar and Cl − . Despite not focusing on the relative shift in photoionization crosssections, our work regarding the 3s Wigner delays between Ar and Cl − seems to agree with the conclusions of [62], since the location and appearance of an "induced Cooper minimum" in Figure 2c was found to be the result of significant coupling to the 3p ionization channels in Argon. However, more work must be carried out to determine the origin for the relative shift in the Cooper minima between isoelectronic systems.

Conclusions
In this work, we performed RRPA Wigner time delay calculations for the halogen ion and noble gas pairings of F − , Ne , Cl − , Ar , (Br − , Kr), and I − , Xe . The individual photoionization channels were then averaged to obtain the average delay times for the highest occupied s, p, and d states. The Wigner delays were plotted with respect to photoelectron energy in order to account for energy shifts due to differences in binding energies. For photoelectron energies below 20 eV, negatively charged halogen ions exhibit large negative Wigner delays that sharply increase. This qualitative difference is due to the absence of a large Coulomb phase, which is known to dominate the Wigner delay behavior of noble gases at low energies. Overall, the qualitative time delay behaviors of halogens and noble gases were found to be similar, with each pairing displaying the same general features and Cooper minima. The location of the 3p Cooper minimum in Cl − was found to be shifted ∼6 eV higher than Ar. A similar shift of ∼10 eV was noticed in the 3s Wigner delays between Ar and Cl − . As atomic mass increased, this relative shift between Cooper minima was found to decrease, except in the case of the 4d Cooper minima of I − , which displayed an ∼11 eV shift above that of Xe. The physical process underlying this relative shift is still not clear, indicating the need for more detailed analysis of negatively charged ions and noble gas atoms in the regions of Cooper minima. Molecular photoionization studies are becoming of greater interest to the attosecond community, with Iodine-containing molecules having already been favored due to iodine's 4d orbitals. Our findings indicate that the average p and s Wigner delays of Br − and I − are the most similar to the average d orbital Wigner delays of their noble gas counterparts Kr and Xe. This implies that Bromine and Iodine are the best halogens for studying time delay phenomena in molecular systems.