Next Article in Journal
Analytical Results for the Three-Body Radiative Attachment Rate Coefficient, with Application to the Positive Antihydrogen Ion H+
Previous Article in Journal
A Precis of Threshold Laws for Positron vs. Electron Impact Ionization of Atoms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Elastic Photon Scattering on Hydrogenic Atoms near Resonances

by
Dmitrii Samoilenko
1,2,3,*,
Andrey V. Volotka
1,2 and
Stephan Fritzsche
1,2,3
1
Helmholtz-Institut Jena, D-07743 Jena, Germany
2
GSI Helmholtzzentrum für Schwerionenforschung, D-64291 Darmstadt, Germany
3
Theoretisch-Physikalisches Institut, Friedrich-Schiller-Universität, D-07743 Jena, Germany
*
Author to whom correspondence should be addressed.
Atoms 2020, 8(2), 12; https://doi.org/10.3390/atoms8020012
Submission received: 29 February 2020 / Revised: 13 April 2020 / Accepted: 16 April 2020 / Published: 18 April 2020

Abstract

:
Scattering of light on relativistic heavy ion beams is widely used for characterizing and tuning the properties of both the light and the ion beam. Its elastic component—Rayleigh scattering—is investigated in this work for photon energies close to certain electronic transitions because of its potential usage in the Gamma Factory initiative at CERN. The angle-differential cross-section, as well as the degree of polarization of the scattered light are investigated for the cases of 1 s 2 p 1 / 2 and 1 s 2 p 3 / 2 resonance transitions in H-like lead ions. In order to gauge the validity and uncertainty of frequently used approximations, we compare different methods. In particular, rigorous quantum electrodynamics calculations are compared with the resonant electric-dipole approximation evaluated within the relativistic and nonrelativistic formalisms. For better understanding of the origin of the approximation, the commonly used theoretical approach is explained here in detail. We find that in most cases, the nonrelativistic resonant electric-dipole approximation fails to describe the properties of the scattered light. At the same time, its relativistic variant agrees with the rigorous treatment within a level of 10% to 20%. These findings are essential for the design of an experimental setup exploiting the scattering process, as well as for the determination of the scattered light properties.

1. Introduction

Absorption and emission of a photon by an atom are the fundamental processes in atomic physics, which are of utmost importance for the characterization of the atomic system. In particular, the laser-induced fluorescence spectroscopy of atoms and molecules is based on an electronic excitation and subsequent fluorescence, which defines the electronic structure of a sample [1]. The energy scale of the fluorescence spectroscopy is typically located between infrared and ultraviolet regions. In order to describe the electronic structure of atoms and ions with binding energies in this energy range, it is often enough to employ the nonrelativistic theory. However, as the interest of atomic physics spreads to a higher energy domain, experimental findings reveal that in some cases, even relativistic consideration might not be enough. One of the most common examples of such a situation is the consideration of core (deeply bound) electrons in heavy atoms or ions. A proper description of the spectral properties involving the strongly bound electrons requires the incorporation of quantum electrodynamics (QED) effects [2]. Both relativistic and QED considerations can be, however, very tedious and demand from the researcher a background in the corresponding fields and sometimes considerable computational resources. At the same time, the well-known and easy-to-use Bohr theory based method of the description of the electronic structure remains attractive in its simplicity. It can provide preliminary estimations even for modern applications. The necessity of more complex considerations is, of course, to be decided in every particular case individually.
Here, we consider the photon absorption and emission (fluorescence spectroscopy) of H-like lead ion Pb 81 + in its ground state with photon energy being close (but not exact) to 1 s 2 p 1 / 2 and 1 s 2 p 3 / 2 resonances. Ab initio, the problem can be formulated in terms of elastic photon scattering or Rayleigh scattering. First, the incident photon is absorbed by an ion, and the ion undergoes the transition from the initial state into an intermediate virtual (or real) state. After that, the ion decays to the final state and emits the scattered photon. Rayleigh scattering is an important and widely studied process in atomic physics. Being one of the basic second-order processes of QED, Rayleigh scattering is extensively studied theoretically [3,4,5,6,7,8,9,10,11], as well as experimentally [12,13,14]. Apart from the fundamental interest, Rayleigh scattering is actively utilized in various fields, including the study of the solid-state, complex molecules, or nano-objects [15,16,17,18], astrophysics [19,20], as well as in medical diagnostics [21].
In this work, we investigate the angular distribution of two observables: cross-section and polarization of the scattered light. We carry out the calculations with three different approaches and compare the obtained results. The simplest approximation we refer to as the nonrelativistic resonant electric dipole one. Electron binding energies and wave functions are obtained by solving the ordinary Schrodinger equation; as an intermediate state, only the nearest resonant level is taken into account; and the electric dipole operator describes the photon-ion interaction. The second approach is the relativistic resonant electric dipole approximation, wherein wave functions and binding energies are the solutions of the Dirac equation. However, we still assume only one (resonant) intermediate state and electric dipole transition operator. Finally, the third consideration we refer to as exact. Here, we add to the relativistic resonant electric dipole approximation two things: (i) sum over all possible intermediate states and (ii) include all terms in the multipole expansion of the transition operator. We believe that the obtained results can be helpful for many applications that exploit the near-resonant Rayleigh scattering. In particular, our example focuses on the Gamma Factory initiative [22]. This project is pursuing building at CERN (European Organization for Nuclear Research) a γ -ray source with photon energies several orders of magnitude higher than available currently anywhere in the world with beam intensities comparable with those produced by free-electron lasers. The idea is to store relativistic highly charged heavy ions at CERN’s Large Hadron Collider (LHC) and scatter optical laser light on these ions. If we choose a laser wavelength to match (keeping in mind the relativistic Doppler effect) a certain electronic transition in the ions, the scattering cross-section (and, consequently, the output intensity) grows drastically. With values of the Lorentz factor γ ≈3000, feasible at LHC, Doppler-enhanced optical frequencies correspond to 1 s 2 p transitions in elements from the lead region of the periodic table. Recently, there have been the first successful proof-of-principle experiments with Pb 81 + ions [23,24], produced with Al stripper foil [25]. For that reason, we have chosen lead among other high-Z elements for consideration.

2. Materials and Methods

2.1. Statement of the Problem

The process under consideration can be illustrated as follows:
| i + γ i | a | f + γ f ,
where the incident photon γ i with polarization vector ϵ i and momentum k i is absorbed by an ion, and the ion is brought from the initial state | i = | 1 s to the intermediate state | a , which in the resonant case is restricted to 2 p 1 / 2 or 2 p 3 / 2 states only. After that, the ion decays to the final state | f = | 1 s and emits the scattered photon γ f with momentum k f and polarization vector ϵ f . The geometry of the described process is shown in Figure 1. The ion is assumed to be located at the origin of the coordinate system. The z-axis in the Cartesian coordinates is oriented along k i , and the x-axis is chosen in the direction of ϵ i if the light is linearly polarized. The scattered photon propagates in the direction of k f . Polar angle ϑ is measured between the propagation directions of the initial and the scattered photons and referred to as the scattering angle. At the same time, the azimuthal angle φ is the angle between the x-axis and the projection of k f on the x y -plane.
The problem can be conveniently formulated in terms of the density matrix formalism. Density operators of the incident ( ρ ^ i ) and the scattered (final) ( ρ ^ f ) photons are connected by a well-known relation [26]:
ρ ^ f = M ^ ρ ^ i M ^ ,
with M ^ being a scattering operator that characterizes the scattering process. The explicit expressions for this operator are given in the following subsection. Since the incident and scattered photons have well defined wave vectors, we write the photon density matrices in the helicity representation as follows:
k λ | ρ ^ | k λ = I 2 1 + P 3 P 1 i P 2 P 1 + i P 2 1 P 3 ,
where I stands for the intensity of the light and the Stokes parameters P 1 , 2 and P 3 define the degree of linear and circular polarization, respectively. For the scattering cross-section calculations, we assume the normalized incident light, i.e., its intensity is equaled to one. Keeping this in mind, the angle-differential cross-section is given by the trace of the density matrix of the scattered photon [26]:
d σ d Ω = k f 1 | ρ ^ f | k f 1 + k f + 1 | ρ ^ f | k f + 1 .
In turn, the degree of linear polarization of the scattered light P l can be found as:
P l = P 1 2 + P 2 2 .
Now, let us turn to the details of the calculations.

2.2. General Theory of Transition Matrix Amplitude

We adopt the Furry picture for our calculations. That means that the electron-nucleus interaction is incorporated in the unperturbed Hamiltonian of the system and the ion-photon interaction is considered as perturbation. For the determination of the properties of the scattered light, we then evaluate the second-order scattering operator M ^ , the matrix elements of which are defined as (in units = m e = c = 1 ) [11,27]:
M ( m i , m f , λ i , λ f ) = α ω a f | α · ϵ λ f * e i k f · r | a a | α · ϵ λ i e i k i · r | i E i + ω E a + f | α · ϵ λ i e i k i · r | a a | α · ϵ λ f * e i k f · r | i E f ω E a ,
where α is the fine structure constant, α is the vector of Dirac matrices, ω stays for the incident and scattered photons energies, and E i , E f , and E a define the energies of the initial, final, and intermediate states, respectively. The quantum numbers m i and m f are the total angular momentum projections of the initial and final states, respectively. Summing over a means that we take into account the contributions of all possible intermediate states instead of just the dominant one. One may notice that if the photon energy ω is close to E a E i , there is a singularity in the first term in Equation (6). This singularity is exactly the case of the resonant intermediate state. In order to regularize the expression, we introduce the corresponding width Γ a in the denominator, i.e., we replace E a with E a i Γ a / 2 in the vicinity of the resonance [28,29]. This prescription corresponds to performing an infinite resummation of the radiative correction insertions for the resonant level [30,31].
Angular properties can be most conveniently described with the help of the multipole expansion of the transition operator [32]. Transition operators inside the matrix elements take the following form:
α · ϵ λ e i k · r = 2 π p J M i J [ J ] ( i λ ) p D M λ J ( φ k , ϑ k , 0 ) α · a J M ( p ) ( ω , r ) ,
where a J M ( p ) ( ω , r ) denotes the electric ( p = 1 ) and magnetic ( p = 0 ) multipole components of the electromagnetic field, J represents the photon angular momentum, and M stands for its projection. The electric and magnetic multipole field components are defined as:
a J M ( 0 ) ( ω , r ) = j J ( ω r ) Y J J M ( r ^ ) ,
a J M ( 1 ) ( ω , r ) = j J 1 ( ω r ) J + 1 2 J + 1 Y J J 1 M ( r ^ ) j J + 1 ( ω r ) J 2 J + 1 Y J J + 1 M ( r ^ ) .
In the above equations, r ^ is the r unit vector, j J ( ω r ) are spherical Bessel functions of order J, and Y J L M ( r ^ ) are vector spherical harmonics, defined by:
Y J L M ( r ^ ) = m μ L m 1 μ | J M Y L m ( r ^ ) χ μ ,
with Clebsch–Gordan coefficients L m 1 μ | J M , ordinary spherical harmonics Y L m ( r ^ ) , and complex spherical unit vectors χ μ expressed via Cartesian unit vectors as:
χ μ = 1 2 ( x ^ i y ^ ) for   μ = 1 , z ^ for   μ = 0 , 1 2 ( x ^ + i y ^ ) for   μ = 1 .
After the expansion (7) is applied, it is obvious that the angular distribution of the emitted radiation is fully defined by the Wigner rotation D-matrices D M λ J ( φ k , ϑ k , 0 ) . As we agreed to the z-axis being coincident with the propagation direction of the incident photon, all the angles included in the corresponding D-matrix are equal to zero and D M λ i J ( 0 , 0 , 0 ) = δ M λ i [33]. Therefore, the angular dependence of the scattering amplitude (6) is defined by a single rotation matrix D M λ f J ( φ , ϑ , 0 ) with φ and ϑ defined as shown in Figure 1. From now on, for brevity, we omit the arguments of the rotation matrix and always assume them to be ( φ , ϑ , 0 ) .
Now, after we write the transition operator in the multipole expansion (7), we can employ the angular momentum theory to deal with the matrix elements effectively. The sum over a in Equation (6) implies not only summation over different states (with different principal quantum numbers n and total angular momentum j n ), but also summation over all possible magnetic substates. These substates have different values of projection m n of the angular momentum j n . A common way to simplify calculations with such a summation is to make use of the Wigner–Eckart theorem:
n 2 j 2 m 2 | α · a J M ( p ) | n 1 j 1 m 1 = ( 1 ) 2 J j 1 m 1 J M | j 2 m 2 n 2 j 2 | | α · a J ( p ) | | n 1 j 1 2 j 2 + 1 .
Applying the Wigner–Eckart theorem (12) together with the multipole expansion (7) to Equation (6), we come to the following form of the scattering amplitude:
M ( m i , m f , λ i , λ f ) = 2 π α ω n j n m n J 1 J 2 M 1 M 2 p 1 p 2 ( i ) J 1 J 2 ( i λ i ) p 1 ( i λ f ) p 2 [ J 1 ] [ J 2 ] D M 2 λ f J 2 * δ λ i M 1 × j n m n J 2 M 2 | j f m f j i m i J 1 M 1 | j n m n [ j n ] 2 f | | α · a J 2 ( p 2 ) | | n j n n j n | | α · a J 1 ( p 1 ) | | i E i + ω E n j n + j n m n J 1 M 1 | j f m f j i m i J 2 M 2 | j n m n [ j f ] [ j i ] f | | α · a J 1 ( p 1 ) | | n j n n j n | | α · a J 2 ( p 2 ) | | i E f ω E n j n
with [ J ] = 2 J + 1 . All the angle dependence is expressed now by the Wigner D-matrix, which depends on the indices J 2 , M 2 , λ f . The function can be effectively dealt with using the angular momentum algebra apparatus.

2.3. Exact Calculations

In order to come now to the observables, we rewrite Equation (2) in terms of the matrix elements of the scattering operator in the following way:
λ f | ρ ^ f | λ f = 1 2 j i + 1 λ i λ i m i m f M ( m i , m f , λ i , λ f ) M * ( m i , m f , λ i , λ f ) λ i | ρ ^ i | λ i .
Here, we assume that the ion is initially unpolarized, and its polarization in the final state remains unobserved. The next part of the derivation is quite tedious, so we will only mention the operations that need to be performed. First, the product of two transition matrix amplitudes M M * in Equation (14) results in a product of two rotation matrices, as follows from Equation (13). This product can be written in terms of a double sum of a single rotation matrix multiplied by two Clebsch–Gordan coefficients. Next, we exploit formulas for the sums of products of Clebsch–Gordan coefficients (to sum over m i , m f , M 1 , M 1 ) and leave only the factors able to reveal some features of the angular distribution. Moreover, for the photon energies much smaller than the electron rest-mass energy, the second term in the square brackets of Equation (13) is much smaller than the first one. For that reason, we write down explicitly only those terms of Equation (14), which include only terms originated from the first term. We band all the other terms together into a single symbolic notation and arrive at the following look of Equation (14):
λ f | ρ ^ f | λ f n j n n j n λ i λ i K M N J 1 J 2 J 1 J 2 p 1 p 2 p 1 p 2 J 2 j f j n j n K J 2 j n j i J 1 J 1 K j n ( i λ i ) p 1 ( i λ f ) p 2 ( i λ i ) p 1 ( i λ f ) p 2 × ( 1 ) j n K ( i ) J 1 J 2 J 1 + J 2 [ J 1 ] [ J 2 ] [ J 1 ] [ J 2 ] D M N K J 2 λ f J 2 λ f | K N J 1 λ i J 1 λ i | K M × f | | α · a J 2 ( p 2 ) | | n j n n j n | | α · a J 1 ( p 1 ) | | i E i + ω E n j n f | | α · a J 2 ( p 2 ) | | n j n n j n | | α · a J 1 ( p 1 ) | | i E f + ω E n j n λ i | ρ ^ i | λ i + O E i + ω E n j n ( n j n ) E f ω E n j n ( n j n ) .
In the above equation, K symbolizes angular momentum transferred to the system as a result of its interaction with two photons and indices M , N effectively play the roles of its projections.
We conclude this section with a short remark concerning the numerics used in this work. As has been mentioned above, the calculation procedure implies a sum over the complete Dirac spectrum. In the numerical calculations, the entire spectrum is typically substituted by a finite set of discrete pseudostates (see, e.g., [34,35]). These pseudostates are then treated in the same manner as the bound states. Setting up the pseudostates requires a set of basis functions. One of the most popular choices of such a basis is a set of splines (piecewise polynomials), as they provide many advantages when it comes to numerical calculations [36]. Within the class of piecewise polynomials, B-splines (for details, see [34]) are especially well suited for atomic physics numerical tasks [35,37], and that is why they are used here as well. In the works [11,38], the B-splines basis set was also employed for the calculation of the Rayleigh scattering.

2.4. Relativistic Resonant Electric-Dipole Approximation

Equation (15) is not restricted to any particular kind of transition or certain states. However, when the photon energy is very close to the transition energy into a particular intermediate state (resonant state), the dominant contribution in the sum over n j n in Equation (15) is given by this state. Moreover, the restriction to the electric dipole transitions only allows us to obtain simple analytical expressions as for the angle-differential cross-section, as well as for the polarization. Under these assumptions, we arrive at the relativistic resonant electric-dipole approximation. To obtain this approximation from Equation (15), we restrict first the sum over n , j n , n , j n into the resonant state only. In the second place, we neglect the terms denoted by O ( ) and imply the electric-dipole approximation J 1 = J 2 = J 1 = J 2 = 1 , as well as p 1 = p 2 = p 1 = p 2 = 1 . With all these changes, we obtain:
λ f | ρ ^ f | λ f rel . resonant E 1 λ i λ i K M N 1 j f j n j n K 1 j n j i 1 1 K j n λ i λ i λ f λ f ( 1 ) j n K D M N K × 1 λ f 1 λ f | K N 1 λ i 1 λ i | K M | f | | α · a 1 ( 1 ) | | n j n | 4 ( E f + ω E n j n ) 2 + Γ n j n 2 / 4 .
For the calculations of the cross-section and degree of polarization of the scattered light, we have to calculate the trace or non-diagonal matrix elements of the resulting density matrix; see Equations (4) and (5). For this, we have to specify the total angular momentum of the intermediate (resonant) state j n . Here, we consider p 1 / 2 and p 3 / 2 states. The first one, the p 1 / 2 state, results in the following angle-differential cross-section:
d σ rel , p 1 / 2 d Ω ( circ . pol . ) d σ rel , p 1 / 2 d Ω ( lin . pol . ) 1
and the linear polarization Stokes parameters:
P 1 rel , p 1 / 2 ( circ . pol . ) = P 1 rel , p 1 / 2 ( lin . pol . ) = 0 , P 2 rel , p 1 / 2 ( circ . pol . ) = P 2 rel , p 1 / 2 ( lin . pol . ) = 0
for circular (circ. pol.) and linear (lin. pol.) polarizations of the incident light. These results can be understood by the following argumentation. For the intermediate state with j n = 1 / 2 6j-symbols in Equation (16), ensure that K = 0 , 1 . Because of the symmetry of the Clebsch–Gordan coefficients, only even values of K will survive, which leaves us with K = 0 and D 00 0 = 1 . Since only the D-matrix defines the angular distribution, for this intermediate state, all scattering directions and linear polarization states are equally possible.
For the second one, the p 3 / 2 state, we start again with Equation (16) and see that K = 2 is now allowed. That means that a second-order Legendre polynomial is appearing in the angular distribution. Derivation results in the following angle-differential cross-section:
d σ rel , p 3 / 2 d Ω ( circ . pol . ) D 00 0 + 1 4 D 00 2 1 + 3 7 cos 2 ϑ
and the linear polarization Stokes parameters:
P 1 rel , p 3 / 2 ( circ . pol . ) = sin 2 ϑ 7 / 3 + cos 2 ϑ , P 2 rel , p 3 / 2 ( circ . pol . ) = 0
for the case of circular polarization of the incident light.
As for the linear polarization case, we obtain:
d σ rel , p 3 / 2 d Ω ( lin . pol . ) 1 3 5 cos 2 φ sin 2 ϑ
and:
P 1 rel , p 3 / 2 ( lin . pol . ) = sin 2 φ cos 2 φ cos 2 ϑ 5 / 3 cos 2 φ sin 2 ϑ , P 2 rel , p 3 / 2 ( lin . pol . ) = 2 sin φ cos φ cos ϑ 5 / 3 cos 2 φ sin 2 ϑ .

2.5. Nonrelativistic Resonant Electric Dipole Approximation

It is well established by now that the nonrelativistic approximation can be employed to obtain rough estimations of the scattering processes. However, in many cases, it fails to describe the experimental results, even qualitatively [39,40,41]. Since the nonrelativistic resonant approximation has been used in the original Gamma Factory proposal [22] (see also [42]), here, we are going to work out the details of this approximation and compare it with previously outlined approaches.
In nonrelativistic treatment 2 p 1 / 2 and 2 p 3 / 2 states are considered as one degenerate 2 p state. We can formally come to this type of consideration from Equation (16) if we characterize electron states with orbital quantum number l instead of total angular momentum quantum number j. For the chosen states, that means the following substitutions: j i , f l i , f = 0 and j n , n l n , n = 1 . After that, we deal with the cross-section and polarization in the same manner as before and obtain well-known expressions:
d σ nonrel d Ω ( circ . pol . ) 1 + cos 2 ϑ ,
P 1 nonrel ( circ . pol . ) = sin 2 ϑ 1 + cos 2 ϑ , P 2 nonrel ( circ . pol . ) = 0 ,
d σ nonrel d Ω ( lin . pol . ) sin 2 φ + cos 2 φ cos 2 ϑ ,
P 1 nonrel ( lin . pol . ) = sin 2 φ cos 2 φ cos 2 ϑ sin 2 φ + cos 2 φ cos 2 ϑ , P 2 nonrel ( lin . pol . ) = 2 sin φ cos φ cos ϑ sin 2 φ + cos 2 φ cos 2 ϑ .
These expressions can be also found in [6,43,44].

3. Results and Discussion

3.1. Ion Rest Frame

We perform numerical calculations for both initial circular and initial linear polarization. Moreover, to investigate the limits of the resonant approximation, we present the results for two different photon energies: at and off the resonance. In the “resonance” scenario, we set the photon energy to match the corresponding transition energy. In the “off-resonance” case, the photon energy is chosen to be 30 Γ 2 p 1 / 2 or 30 Γ 2 p 3 / 2 above the corresponding resonances 2 p 1 / 2 or 2 p 3 / 2 , respectively. For the case of the H-like lead ion, these widths take the following values Γ 2 p 1 / 2 = 19.4 eV and Γ 2 p 3 / 2 = 16.9 eV. The resonance and off-resonance scenarios lead to essentially different values of the cross-section. Since we are mostly interested in the form of the angular distribution, we normalize the cross-section results for all scenarios and compare the line shapes. Here, we note that the shapes obtained via the nonrelativistic and relativistic resonant electric-dipole approximations are independent of the photon energy and are given by Equations (23)–(26) and (17)–(22), respectively.
Figure 2 shows the cross-section and polarization angular distributions in the case when the incident light is circularly polarized. For the 2 p 1 / 2 intermediate resonance (left column), we see that for the “resonance” photon energy, the relativistic resonant electric-dipole approximation (black dashed line) almost coincides with the exact calculations (solid blue line). The results differ for less than 0.3% of the mean value, and this difference is not visible in the figure. We can also see that both cross-section and polarization distributions are flat, as expected form Equations (17) and (18). Since the nondipole photon-ion interaction has no effect on the | 2 p 1 / 2 | 1 s 1 / 2 transition, the 2 p 1 / 2 intermediate state is a perfect candidate for the demonstration of an effect beyond the resonant approximation. One can see that in the “off-resonance” case, the angle-differential cross-section deviates from the flat line by around 20%. At the same time, the degree of linear polarization reaches 0.25 for ϑ = 90 . The middle column in Figure 2 shows the results for the 2 p 3 / 2 resonant state. First of all, we notice that the shape of the cross-section distribution obtained by the exact calculations differs significantly from the relativistic resonant electric-dipole approximation, Equations (19) and (20). Comparing the cross-section distributions in the “resonance” and “off-resonance” cases, we see that even that far from the resonance (30 line widths), the effect of nondipole contributions is more pronounced. On the other hand, the polarization seems to be more sensitive to the contributions from nonresonant states. When we consider 2 p as the intermediate state, we imply nonrelativistic approximation. In this approximation, the line shapes of both cross-section and polarization distributions do not change with changing the photon energy and agree with Equations (23) and (24).
In Figure 3, we present the results for the linearly polarized incident light. When we switch from circular to linear polarization, the system loses circular symmetry in the x y -plane. That means that the results are different for different values of azimuthal observation angle φ . Here, we consider only the case of φ = 0 . The discussion of the results would be quite similar to the previous figure, so we only make several observations: (i) for the 2 p 1 / 2 resonant state, we see a similar picture as in Figure 2, except that the effect of the nonresonant states is even more pronounced; (ii) for the 2 p 3 / 2 resonant state, we note that the results in the relativistic resonant electric-dipole approximation follow Equation (21) and Equation (22); one can also see that both distributions, in this case, vary significantly more strongly than for the circular polarization; (iii) for the nonrelativistic case, the shapes of both lines are essentially different as compared to the 2 p 1 / 2 and 2 p 3 / 2 cases. In accordance with Equations (25) and (26), they are also different from the case of circular initial polarization.
The large detuning of the photon energy in the “off-resonance” scenario is a limiting case and considered here for studies of the importance of the nonresonant state’s contribution. For the smaller values of the detuning, the nonresonant state’s contribution is lower. Specifically, in all the cases considered above for photon energies in the range of 5 Γ 2 p 1 / 2 or 5 Γ 2 p 3 / 2 from the resonance, the results of the relativistic resonant electric-dipole approximation differ from exact calculations for less than 12%. Therefore, we expect that this relativistic approximation performs well enough for the design of experiments [22], in contrast to the presently employed nonrelativistic one.

3.2. Laboratory Reference Frame

Since the primary application of the results obtained in this work is expected to be the Gamma Factory initiative, we make two additional points. The first one is that the results presented above can be transformed into the laboratory reference frame for any given value of Lorentz factor γ . Assuming that the incident photon and the ion move towards each other (so, the ion also moves along the z-axis), this can be done by combining any angular distribution with the well-known relativistic aberration formula:
cos ϑ = cos ϑ β 1 β cos ϑ ,
where ϑ is the scattering angle in the laboratory reference frame counted from the z-axis and β is the ion velocity relative to the speed of light. The latter can be expressed via γ from:
γ = 1 1 β 2 .
In this head-on collision scenario, the azimuthal angle φ transforms trivially to φ , i.e., φ = φ . Corresponding calculations are presented in Figure 4 for the case of circular initial polarization and 2 p 3 / 2 taken as the intermediate state. One can see from the figure that with growing γ , the distribution is concentrated in the vicinity of ϑ = π . For a fixed high value of γ , one obtains a well-known narrow cone-like radiation distribution, which is characteristic for relativistically moving particles. Moreover, as we discussed in the previous subsection, the relativistic resonant electric-dipole approximation provides a reasonable estimation of the angle distributions. It means that for the design of experiments, one can transform to the laboratory frame the distributions given by Equations (17)–(22).
The second point is that, although we so far considered a single process (meaning that we have just one scattered photon), in the experimental case, one has a full beam of photons. Then, the value of the differential cross-section in a certain direction describes the number of photons emitted in that direction. Similarly, the value of P l ( ϑ ) describes the polarization state of the photon beam in the direction of ϑ . Because of the concentration of the radiation in the direction of ϑ = π , all the photons at different polarization states form a single beam. The polarization properties of this beam are still to be defined. The recipe for how to define them will be discussed elsewhere, but we point out the obvious importance of the consideration of the polarization distribution for studying this issue.

4. Conclusions

In this work, we investigated the excitation and subsequent decay of the 2 p electron states of the H-like lead ion using the well-developed theoretical apparatus of scattering theory. The investigation included the angular distribution of the cross-section and the degree of linear polarization of the scattered light. We derived the exact expression for density matrix elements, which were necessary for the calculation of these two observables. The expressions were explicitly tailored for two resonant electric-dipole approximations: nonrelativistic and relativistic. We performed numerical calculations with the exact formula, as well as with both approximations. The results were compared for different intermediate (resonant) states 2 p 1 / 2 and 2 p 3 / 2 , as well as for different photon energies with respect to the resonance. The findings showed two main points. First, in most cases, the results of the nonrelativistic approach were substantially different from the exact calculations. Second, for moderate detuning of the photon energy from the resonance transition energy, the relativistic resonant electric-dipole approximation provided accuracy around 10% for all considered scenarios. The analytical expressions derived in this approximation were quite general and could be used for other types of H-like ions, as well. Therefore, Equations (17), (19), and (21) derived for the cross-section and Equations (18), (20), and (22) for polarization could be used in the design of near-resonance Rayleigh scattering experiments and the interpretation of the results. In particular, we expect these formulas to be favorable to theoretical investigations at Gamma Factory. Additionally, it was demonstrated that the derived analytical expressions in the relativistic resonant electric-dipole approximation could be easily transformed into the laboratory reference frame.

Author Contributions

Investigation, D.S.; methodology, A.V.V.; supervision, S.F.; writing, original draft, D.S.; writing, review and editing, A.V.V. and S.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Demtröder, W. Laser Spectroscopy: Basic Concepts and Instrumentation; Springer: Heidelberg, Germany, 2003. [Google Scholar]
  2. Indelicato, P. QED tests with highly charged ions. J. Phys. B At. Mol. Opt. Phys. 2019, 52, 232001. [Google Scholar] [CrossRef] [Green Version]
  3. Franz, W. Rayleighsche Streuung harter Strahlung an schweren Atomen. Z. Phys. 1936, 98, 314. [Google Scholar] [CrossRef]
  4. Brown, G.E.; Peierls, R.E.; Woodward, J.B. The coherent scattering of γ-rays by K electrons in heavy atoms—I. Method. Proc. R. Soc. Lond. A 1954, 227, 51. [Google Scholar]
  5. Johnson, W.R.; Feiock, F.D. Rayleigh Scattering and the Electromagnetic Susceptibility of Atoms. Phys. Rev. 1968, 168, 22. [Google Scholar] [CrossRef]
  6. Kane, P.P.; Kissel, L.; Pratt, R.H.; Roy, S.C. Elastic scattering of γ-rays and X-rays by atoms. Phys. Rep. 1986, 140, 75. [Google Scholar] [CrossRef]
  7. Pratt, R.H.; Bergstrom, P.M., Jr.; Kissel, L. New relativistic S-matrix results for scattering —Beyond the usual anomalous factors/beyond impulse approximation. In Resonant Anomalous X-ray Scattering: Theory and Applications; Materlik, G., Sparks, C.J., Fischer, K., Eds.; Elsevier: Amsterdam, The Netherlands, 1994; pp. 9–33. [Google Scholar]
  8. Kissel, L.; Zhou, B.; Roy, S.C.; Sen Gupta, S.K.; Pratt, R.H. The validity of form-factor, modified-form-factor and anomalous-scattering-factor approximations in elastic scattering calculations. Acta Cryst. 1995, A51, 271. [Google Scholar] [CrossRef] [Green Version]
  9. Carney, J.P.J.; Pratt, R.H.; Manakov, N.L.; Meremianin, A.V. Dependence of photon-atom scattering on energy resolution and target angular momentum. Phys. Rev. A 2000, 61, 042704. [Google Scholar] [CrossRef]
  10. Surzhykov, A.; Yerokhin, V.A.; Stöhlker, T.; Fritzsche, S. Rayleigh X-ray scattering from many-electron atoms and ions. J. Phys. B At. Mol. Opt. Phys. 2015, 48, 189501. [Google Scholar] [CrossRef] [Green Version]
  11. Volotka, A.V.; Yerokhin, V.A.; Surzhykov, A.; Stöhlker, T.; Fritzsche, S. Many-electron effects on X-ray Rayleigh scattering by highly charged He-like ions. Phys. Rev. A 2016, 93, 023418. [Google Scholar] [CrossRef] [Green Version]
  12. Jung, M.; Dunford, R.W.; Gemmell, D.S.; Kanter, E.P.; Krässig, B.; LeBrun, T.W.; Southworth, S.H.; Young, L.; Carney, J.P.J.; LaJohn, L.; et al. Manifestations of Nonlocal Exchange, Correlation, and Dynamic Effects in X-ray Scattering. Phys. Rev. Lett. 1998, 81, 8. [Google Scholar] [CrossRef] [Green Version]
  13. Young, L.; Dunford, R.W.; Kanter, E.P.; Krässig, B.; Southworth, S.H.; Bonham, R.A.; Lykos, P.; Morong, C.; Timm, A.; Carney, J.P.J.; et al. Corrections to the usual X-ray scattering factors in rare gases: Experiment and theory. Phys. Rev. A 2001, 63, 052718. [Google Scholar] [CrossRef]
  14. Blumenhagen, K.-H.; Fritzsche, S.; Gassner, T.; Gumberidze, A.; Märtin, R.; Schell, N.; Seipt, D.; Spillmann, U.; Surzhykov, A.; Trotsenko, S.; et al. Polarization transfer in Rayleigh scattering of hard X-rays. New J.Phys. 2016, 18, 103034. [Google Scholar] [CrossRef]
  15. Sfeir, M.Y.; Wang, F.; Huang, L.; Chuang, C.-C.; Hone, J.; O’Brien, S.P.; Heinz, T.F.; Brus, L.E. Probing electronic transitions in individual carbon nanotubes by Rayleigh scattering. Science 2004, 306, 1540. [Google Scholar] [CrossRef] [Green Version]
  16. Kampel, N.S.; Griesmaier, A.; Hornbak Steenstrup, M.P.; Kaminski, F.; Polzik, E.S.; Müller, J.H. Effect of Light Assisted Collisions on Matter Wave Coherence in Superradiant Bose-Einstein Condensates. Phys. Rev. Lett. 2012, 108, 090401. [Google Scholar] [CrossRef] [Green Version]
  17. Kulik, L.V.; Ovchinnikov, K.; Zhuravlev, A.S.; Bisti, V.E.; Kukushkin, I.V.; Schmult, S.; Dietsche, W. Resonant Rayleigh scattering as a probe of spin polarization in a two-dimensional electron system. Phys. Rev. B 2012, 85, 113403. [Google Scholar] [CrossRef]
  18. Wu, W.; Yue, J.; Li, D.; Lin, X.; Zhu, F.; Yin, X.; Zhu, J.; Dai, X.; Liu, P.; Wei, Y.; et al. Interface dipole enhancement effect and enhanced Rayleigh scattering. Nano Res. 2015, 8, 303. [Google Scholar] [CrossRef]
  19. Maeda, K.; Terada, Y.; Kasen, D.; Röpke, F.K.; Bamba, A.; Diehl, R.; Nomoto, K.; Kromer, M.; Seitenzahl, I.R.; Yamaguchi, H.; et al. Prospect of Studying Hard X- and Gamma-Rays from Type Ia Supernovae. Astrophys. J. 2012, 760, 54. [Google Scholar] [CrossRef] [Green Version]
  20. The, L.-S.; Burrows, A. Expectations for the Hard X-ray Continuum and Gamma-Ray Line Fluxes from the Type Ia Supernova SN 2014J in M82. Astrophys. J. 2014, 786, 141. [Google Scholar] [CrossRef] [Green Version]
  21. Zhang, X.-D.; Wu, D.; Shen, X.; Chen, J.; Sun, Y.-M.; Liu, P.-X.; Liang, X.-J. Size-dependent radiosensitization of PEG-coated gold nanoparticles for cancer radiation therapy. Biomaterials 2012, 33, 6408. [Google Scholar] [CrossRef] [Green Version]
  22. Krasny, M.W. The Gamma Factory proposal for CERN. arXiv 2015, arXiv:1511.07794. [Google Scholar]
  23. Płaczek, W.; Abramov, A.; Alden, S.E.; Alemany Fernandez, R.; Antsiferov, P.S.; Apyan, A.; Bartosik, H.; Bessonov, E.G.; Biancacci, N.; Bieroń, J.; et al. Gamma Factory at CERN— Novel Research Tools Made of Light. Acta Phys. Pol. B 2019, 50, 1191. [Google Scholar] [CrossRef]
  24. Budker, D.; Crespo López-Urrutia, J.R.; Derevianko, A.; Flambaum, V.V.; Krasny, M.W.; Petrenko, A.; Pustelny, S.; Surzhykov, A.; Yerokhin, V.A.; Zolotorev, M. Atomic physics studies at the Gamma Factory at CERN. arXiv 2020, arXiv:2003.03855. [Google Scholar]
  25. Kröger, F.M.; Weber, G.; Shevelko, V.P.; Hirlander, S.; Alemany-Fernandez, R.; Cornelis, K.; Goddard, B.; Velotti, F.M.; Krasny, M.W.; Stöhlker, T. Charge State Tailoring of Highly Relativistic Heavy ion Beams—Comparison of Theory and Experiment; Helmholtz Institute Jena Annual Report; Helmholtz-Institut: Jena, Germany, 2018; p. 62. [Google Scholar]
  26. Blum, K. Density Matrix Theory and Applications; Springer: Heidelberg, Germany, 2012. [Google Scholar]
  27. Akhiezer, A.I.; Berestetskii, V.B. Quantum Electrodynamics; Wiley: New York, USA, 1965. [Google Scholar]
  28. Wu, Z.W.; Volotka, A.V.; Surzhykov, A.; Dong, C.Z.; Fritzsche, S. Level sequence and splitting identification of closely spaced energy levels by angle-resolved analysis of fluorescence light. Phys. Rev. A 2016, 93, 063413. [Google Scholar] [CrossRef] [Green Version]
  29. Wu, Z.W.; Volotka, A.V.; Surzhykov, A.; Fritzsche, S. Angle-resolved X-ray spectroscopic scheme to determine overlapping hyperfine splittings in highly charged heliumlike ions. Phys. Rev. A 2017, 96, 012503. [Google Scholar] [CrossRef] [Green Version]
  30. Shabaev, V.M. Two-time Green’s function method in quantum electrodynamics of high-Z few-electron atoms. Phys. Rep. 2002, 356, 119. [Google Scholar] [CrossRef] [Green Version]
  31. Andreev, O.Y.; Labzowsky, L.N.; Plunien, G.; Solovyev, D.A. QED theory of the spectral line profile and its applications to atoms and ions. Phys. Rep. 2008, 455, 135. [Google Scholar] [CrossRef]
  32. Rose, M.E. Elementary Theory of Angular Momentum; Wiley: New York, USA, 1957. [Google Scholar]
  33. Varshalovich, D.; Moskalev, A.; Khersonskii, V. Quantum Theory of Angular Momentum; World Scientific Publishing: Singapore, 1988. [Google Scholar]
  34. De Boor, C. A Practical Guide to Splines; Bloch, A., Epstein, C.L., Goriely, A., Greengard, L., Eds.; Springer: New York, NY, USA, 1978. [Google Scholar]
  35. Sapirstein, J. The use of basis splines in theoretical atomic physics. J. Phys. B At. Mol. Opt. Phys. 1996, 29, 5213. [Google Scholar] [CrossRef]
  36. Bachau, H.; Cormier, E.; Decleva, P.; Hansen, J.E.; Martin, F. Applications of B-splines in atomic and molecular physics. Rep. Prog. Phys. 2001, 64, 1815. [Google Scholar] [CrossRef]
  37. Froese, F.C. B-Splines in Variational Atomic Structure Calculations; Arimondo, E., Berman, P.R., Lin, C.C., Eds.; Academic Press: New York, NY, USA, 2008; pp. 235–291. [Google Scholar]
  38. Safari, L.; Amaro, P.; Fritzche, S.; Santos, J.P.; Fratini, F. Relativistic total cross-section and angular distribution for Rayleigh scattering by atomic hydrogen. Phys. Rev. A 2012, 85, 043406. [Google Scholar] [CrossRef] [Green Version]
  39. Chen, M.H.; Scofield, J.H. Relativistic effects on angular distribution and polarization of dielectronic satellite lines of hydrogenlike ions. Phys. Rev. A 1995, 52, 2057. [Google Scholar] [CrossRef]
  40. Surzhykov, A.; Fritzsche, S.; Gumberidze, A.; Stöhlker, T. Lyman-α1 Decay in Hydrogenlike Ions: Interference between the E1 and M2 Transition Amplitudes. Phys. Rev. Lett. 2002, 88, 153001. [Google Scholar] [CrossRef]
  41. Surzhykov, A.; Fritzsche, S.; Stöhlker, T. Polarization of the Lyman–α1 Line Following the Radiative Recombination of Bare, High-Z Ions. Hyperfine Interact. 2003, 146, 35. [Google Scholar] [CrossRef]
  42. Bessonov, E.G. Light sources based on relativistic ion beams. Nucl. Instrum. Meth. B 2013, 309, 92. [Google Scholar] [CrossRef] [Green Version]
  43. Surzhykov, A.; Yerokhin, V.A.; Jahrsetz, T.; Amaro, P.; Stöhlker, T.; Fritzsche, S. Polarization correlations in the elastic Rayleigh scattering of photons by hydrogenlike ions. Phys. Rev. A 2013, 88, 062515. [Google Scholar] [CrossRef]
  44. Manakov, N.L.; Meremianin, A.V.; Maquet, A.; Carney, J.P.J. Photon-polarization effects and their angular dependence in relativistic two-photon bound-bound transitions. J. Phys. B At. Mol. Opt. Phys. 2000, 33, 4425. [Google Scholar] [CrossRef]
Figure 1. (Color online) The adopted geometry for the scattering process. The z-axis is chosen along the direction of the momentum of the incident photon k i , while the x-axis is fixed by its polarization vector ϵ i . The polar (scattering) ϑ and azimuthal φ angles uniquely define the direction of the scattered photon k f .
Figure 1. (Color online) The adopted geometry for the scattering process. The z-axis is chosen along the direction of the momentum of the incident photon k i , while the x-axis is fixed by its polarization vector ϵ i . The polar (scattering) ϑ and azimuthal φ angles uniquely define the direction of the scattered photon k f .
Atoms 08 00012 g001
Figure 2. Cross-section and polarization angular distribution for different resonant states: 2 p 1 / 2 —left column, 2 p 3 / 2 —middle column, and nonrelativistic 2 p state—right column. Different approaches are compared: dashed black line—relativistic resonant electric dipole approximation, dashed-dotted red line—exact calculations for the “off-resonance” scenario (see the text), solid blue line—exact calculations at the corresponding resonance. Incident light is assumed to be circularly polarized.
Figure 2. Cross-section and polarization angular distribution for different resonant states: 2 p 1 / 2 —left column, 2 p 3 / 2 —middle column, and nonrelativistic 2 p state—right column. Different approaches are compared: dashed black line—relativistic resonant electric dipole approximation, dashed-dotted red line—exact calculations for the “off-resonance” scenario (see the text), solid blue line—exact calculations at the corresponding resonance. Incident light is assumed to be circularly polarized.
Atoms 08 00012 g002
Figure 3. The same as Figure 2, but for the linearly polarized incident light and for the azimuthal observation angle φ = 0 .
Figure 3. The same as Figure 2, but for the linearly polarized incident light and for the azimuthal observation angle φ = 0 .
Atoms 08 00012 g003
Figure 4. Angular distribution of the scattered light in the laboratory reference frame at different values of the Lorentz factor. The color illustrates the values of the cross-section, normalized to one for a fixed value of γ .
Figure 4. Angular distribution of the scattered light in the laboratory reference frame at different values of the Lorentz factor. The color illustrates the values of the cross-section, normalized to one for a fixed value of γ .
Atoms 08 00012 g004

Share and Cite

MDPI and ACS Style

Samoilenko, D.; Volotka, A.V.; Fritzsche, S. Elastic Photon Scattering on Hydrogenic Atoms near Resonances. Atoms 2020, 8, 12. https://doi.org/10.3390/atoms8020012

AMA Style

Samoilenko D, Volotka AV, Fritzsche S. Elastic Photon Scattering on Hydrogenic Atoms near Resonances. Atoms. 2020; 8(2):12. https://doi.org/10.3390/atoms8020012

Chicago/Turabian Style

Samoilenko, Dmitrii, Andrey V. Volotka, and Stephan Fritzsche. 2020. "Elastic Photon Scattering on Hydrogenic Atoms near Resonances" Atoms 8, no. 2: 12. https://doi.org/10.3390/atoms8020012

APA Style

Samoilenko, D., Volotka, A. V., & Fritzsche, S. (2020). Elastic Photon Scattering on Hydrogenic Atoms near Resonances. Atoms, 8(2), 12. https://doi.org/10.3390/atoms8020012

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