1. Introduction
The Compton scattering process is the scattering between a bound electron (
e) in an atomic or molecular system and an incident photon (
) in the electromagnetic field
It is one of the most important and mysterious electromagnetic processes in physics. This process was first noticed by A. H. Compton in 1923 [
1,
2], through which the quantum nature of X-rays was revealed successfully. Since Compton scattering was discovered in the 1920s, it has been carefully studied and extensively investigated for almost a century. The study of atomic Compton scattering could give us information on the interaction between light and matter, and it can also provide opportunities to reveal the underling nature of electric structures, electron correlations, electron momentum distributions, and other aspects of physics [
3,
4,
5]. Because of wide applications, Compton scattering is always an attractive topic in atomic physics, molecular physics, condensed matter physics, nuclear physics and material science [
5,
6,
7,
8,
9] (
Figure 1).
Since the 1920s, great efforts have been made to develop theoretical methods on ab initio calculations for the Compton scattering process. The simplest approach is the free electron approximation (FEA), which was first developed by O. Klein and Y. Nishina in 1929 [
10]. In the FEA approach, the electrons in atomic or molecular systems are treated as free electrons. The atomic binding, electron shielding, electron correlations, electron motions around atomic nucleus, and other many-body effects are all omitted for simplicity. In FEA, the angular distribution for Compton scattering process is given by the Klein-Nishina formula [
10,
11]
Furthermore, in FEA approach, the final state photon energy
after the scattering process is totally determined by the scattering angle
in Compton scattering process
Here,
is known as the Compton energy. Because of computational simplicity and clear physical meaning, the FEA approach is most widely used in many branches of science. over the past few decades, it has become a standard treatment for atomic Compton scattering in many textbooks. However, the FEA approach is much too simple, very little information on electronic structures and property of materials can be acquired. From Equation (
1), it is evident that the angular distribution
is independent of the electron momentum distributions in target materials.
Soon after the FEA approach was formulated by O. Klein and Y. Nishina in 1929, other approaches considering bound structures in atoms and molecules emerged. It was J. W. M. DuMond first realized that the atomic effects and electron motions would give inevitable effects on the scattering process. Then he introduced an approach to effectively include atomic binding effects as well as electron motion around atomic nuclei in the calculation of Compton scattering [
12,
13,
14]. His approach was later known as the impulse approximation (IA). The DuMond’s approach is the nonrelativistic impulse approximation (NRIA), and the same result was re-derived by P. Eisenberger et al. in the 1970s [
15,
16,
17,
18]. The relativistic impulse approximation (RIA) was developed until the 1970s–1980s by R. Ribberfors et al. [
18,
19,
20,
21,
22,
23]. In the IA formulation, the doubly differential cross-section (DDCS) for Compton scattering process can be factorized into two parts as follows:
where
is the solid angle for final state outgoing photon. In Equation (
3),
is a factor dependent on kinematical and dynamical properties of atomic Compton scattering, and irrelevant to the electronic structure of target materials. The factor
, which is known as the Compton profile, is related to the momentum distributions of bound electrons in the atomic or molecular systems [
24]. The Compton profile can reflect electronic structures, physical and chemical properties of target materials. It can be determined from theoretical calculations and experimental measurements. In actual ab initio calculations, the Compton profile
can be obtained using the nonrelativistic Hartree–Fock theory (HF), the relativistic Dirac–Hartree–Fock theory (DHF), and the density functional theory (DFT). It is worth noting that in the IA approach, all the many-body effects in atomic and molecular systems in Compton scattering processes are incorporated into Compton profiles
.
From the IA approach, the angular distribution of Compton scattering process can be obtained by integration of Equation (
3), and the results is also a simple correction for FEA:
Here,
is the angular distribution calculated in FEA approach as shown in Equation (
1), and the correction factor
is called the scattering function. It should be mentioned that apart from IA, there are other approaches which can give the same results in Equation (
4). In other words, scattering function
can also be calculated by other methods, such as the nonrelativistic Waller-Hartree theory [
25,
26]. This kind of approach, in which the angular distribution of Compton scattering is given by Equation (
4), is called the incoherent scattering function or incoherent scattering factor (ISF) approach. Same as Compton profile
, the scattering function
can also reflect the electronic properties of target materials. On the one hand,
and
can be measured from Compton scattering experiments with high precision. On the other hand, they can also be predicted by theoretical ab initio calculations in atomic, molecular, and condensed matter physics [
24,
27]. The Compton profile
and scattering function
provide a bridge between Compton scattering and interdisciplinary studies in many branches of science. They can offer opportunities to learn the electronic structure and properties of target materials from Compton scattering process. Recently, there are many studies in which electronic structures and properties (include electron interactions, electron correlations, electron momentum distributions, band structures and fermi surfaces) are investigated through Compton profile and scattering function [
6,
7,
8,
9].
In recent decades, several methods beyond IA and ISF were emerged and formulated [
28,
29,
30,
31,
32,
33,
34]. Among these approaches, the most successful one is based on the perturbation expansion of many-body electromagnetic interactions in atomic or molecular systems. In this approach, the DDCS of Compton scattering process can be calculated by the scattering matrix of a many-body theory
This approach is called the Scattering Matrix (SM) approach. In the SM, scattering matrix of Compton scattering process
can be calculated through the many-body interaction Hamiltonian
, which gives the many-body interaction between photon and bound electrons in atomic or molecular systems. The SM approach can take into account the factor of atomic bindings and electron interactions in atomic or molecular systems as much as possible. SM could handle with these many-body effects precisely in the dynamical process of Compton scattering. Recently, the SM formulation has revealed many nontrivial properties of Compton scattering, and it has attracted lots of interests in interdisciplinary studies [
29,
32].
Furthermore, there are several experiments which give supporting to the SM predictions [
35,
36,
37,
38,
39,
40,
41]. Recently, Max Kircher et al. conducted a kinematically complete Compton scattering experiment using X-rays produced from accelerators with energy about 2.1 keV. By measuring the angular distribution of scattered photons, the experimental observations present large deviations with the FEA results. However, when SM approach is employed and ab initio numerical calculations are carried out, the experimental data are consistent with theoretical predictions [
41]. These results indicate that the SM approach is becoming a promising tool to duel with Compton scattering with bound electrons, and it may have great impacts in this area in the near future.
This paper is organized as follows.
Section 1 gives an introduction on Compton scattering and the development of theoretical methods on ab initio calculations of Compton scattering process.
Section 2 gives a description of FEA,
Section 3 discuss the IA approach, and
Section 4 is devoted to the ISF. The most advanced SM approach is presented in
Section 5.
Section 6 presents the comparisons between theoretical calculations and experimental measurements. Databases and applications for Compton scattering are briefly introduced in
Section 7. Summaries are presented in
Section 8.
2. Free Electron Approximation
In the free electron approximation (FEA), the electrons in Compton scattering are treated as free electrons, all the atomic binding effects and many-body interactions of bound electrons in atomic or molecular systems are neglected in the scattering process. Furthermore, it is also assumed that the bound electrons are at rest before the Compton scattering.
The FEA approach for Compton scattering process was given by O. Klein and Y. Nishina in 1929 [
10]. In this formulation, the scattered photon energy
in Compton scattering process is totally determined by its scattering angle
via equation
Here,
is the mass of electron, and
is known as the Compton energy. When
, the energy of scattered photon
reaches its minimum, meanwhile the energy transfer
arrives at its maximum. They can be expressed in the following
They correspond to the Compton edge of the Compton scattering spectrum of
[
42,
43].
In FEA formulation, the angular distribution of Compton scattering process is given by the Klein-Nishina formula [
10,
11]
and the corresponding DDCS can be expressed as
Here,
is the classical radius of electron,
and
are polarization vectors for incoming and outgoing photons,
and
are the energy and solid angle of the outgoing scattered photon, respectively. In the FEA results, due to the Dirac delta function
in Equation (
10), the spectrum of DDCS is an isolated line located at Compton energy
, which is known as “Compton line”, as illustrated in
Figure 2.
In the nonrelativistic limit, where the incoming photon energy
, the angular distribution of Compton scattering reduced to
In this case, the Compton energy becomes
where
K is the momentum transfer in the Compton scattering process. Furthermore, in the elastic scattering limit
(namely the scattering angle
), the Equation (
11) further reduce to the Thomson formula
This is the differential cross-section of elastic scattering process between photon and free electron [
11,
24,
27,
44].
It should be mentioned that the Klein-Nishina formula in Equation (
9) in FEA works well only in the cases that atomic binding energies are negligible, and the atomic electrons are approximately free. When the incident photon energy
and energy transfer
T are extremely high, the FEA result is appropriate. However, in many cases, the incident photon energy is comparable to the X-ray characteristic energies, in which the atomic bindings and electron many-body interactions cannot be neglected. Then the FEA formulation becomes invalid, and it fails to fit the experimental data [
28].
3. Impulse Approximation
In the impulse approximation (IA) method, the atomic binding effects are effectively considered, and the electron pre-collision motions around atomic nuclei are also included. The basic starting point of IA approach can be shown in the following way. Because of atomic binding and electrons many-body interactions, the bound electrons in atomic and molecular systems have a momentum distribution
when moving around atomic or molecular nuclei. In principle, the momentum distribution for bound electrons is determined by ground state wavefunctions in atomic or molecular systems. In the Compton scattering process, suppose electron in momentum eigenstate
is scattered with incoming photon
very rapidly, like an impulse acting on the electron. This scattering process is too quick to be disturbed by other electrons. In this way, electron momentum eigenstate
scattered with incident photon
independently as free electrons. Many-body interactions and interference terms between electrons with different momentum eigenstates (
and
) in the dynamical process of Compton scattering are omitted for simplicity. Therefore, the DDCS of the Compton scattering process is achieved through a summation of the scattering probability for all possible momentum eigenstates
In this expression,
is the DDCS of Compton scattering between photon
and electron in momentum eigenstate
, which can be calculated by FEA. The
and
are the energies of electron before and after Compton scattering process. The Equation (
14) can be viewed as the basic starting point of IA approach, and different versions of IA treatments are achieved by applying different numerical schemes in calculating Equation (
14) [
42].
The aforementioned atomic binding effects and electron pre-collision motion make the IA results significantly different from FEA results. In the FEA results, there is a one-to-one correspondence between the final state photon energy
and the scattering angle
(see Equation (
6)). In the IA approach, due to the electron motion around atomic nuclei, the energy of the scattered photon
cannot be totally determined by its scattering angle
as in FEA. For the same scattering angle
, the outgoing photon energy
has a continuous distribution. The maximum cross-section is located at the Compton energy
, which is given by Equation (
6), forming a “Compton peak” in the spectrum of DDCS. In
Figure 2, the DDCS of Compton scattering predicted in FEA and IA approaches are illustrated schematically. Sometimes, it is more enlightening to interpret the emergence of “Compton peak” as the Doppler broadening effect of the “Compton line” (located at Compton energy
) [
45,
46]. This Doppler effect is due to the bound electrons’ motion around atomic nuclei.
The nonrelativistic impulse approximation (NRIA) was first developed by DuMond in the 1930s [
12,
13,
14]. In the 1970s, Eisenberger and Platzman re-derived the NRIA formulation based on the more accurate scattering matrix (SM) approach. It is shown that NRIA can be viewed as the leading order approximation of SM approach [
15,
16,
17,
18]. In NRIA, the DDCS of Compton scattering is given by [
14,
15]
Here,
K is the modulus of the momentum transfer vector
in Compton scattering process, and
is given by
In Equation (
15), the factor
relies on the dynamical and kinematical properties of Compton scattering process in the nonrelativistic limit. The
is a factor coming from the many-body effects in atomic or molecular systems, which is known as the “Compton profile” [
24]
where
is the electron momentum density of the atomic or molecular ground-states. The Compton profile is closely related to the electronic properties of the atomic or molecular systems, and it has been widely studied in atomic, molecular and condensed matter physics [
6,
7,
9]. Furthermore, for most of the atomic systems, the momentum distribution is spherical symmetric, then the Compton profile reduces to
In these cases, the Compton profile is bell-shaped and axisymmetric around the axis. However, many of the molecular or condensed matter systems do not have this spherical symmetric property.
The relativistic impulse approximation (RIA) was developed by R. Ribberfors et al. in 1975–1985 [
18,
19,
20,
21,
22,
23]. In this formulation, the DDCS of Compton scattering is given by [
21,
22,
47]
The same as in NRIA,
is a factor depends on kinematical and dynamical properties of Compton scattering. The factor relevant to atomic or molecular structure is the Compton profile
defined in Equation (
17). In Equation (
19),
K is the modulus of the momentum transfer vector
and
is the projection of the electron’s initial momentum on the momentum transfer direction
Moreover, the function
is defined to be
with
and
defined as
To summarize, combining Equations (
15) and (
19), it is clear that the DDCS of Compton scattering process calculated in IA approach factorizes into two parts
where
is a factor dependent on kinematical and dynamical properties of Compton scattering process, and the Compton profile
is relevant to the electron momentum distribution of the atomic or molecular ground-states. It is worth noting that in the IA approach, all the many-body effects in atomic and molecular systems in Compton scattering process are incorporated into Compton profiles
. The Compton profile can reflect electronic structures and properties of target materials. Therefore, it provides us an opportunity to learn the electronic structures, electron interactions and properties of materials from Compton scattering processes. In actual ab initio calculations, the Compton profile
can be calculated by the nonrelativistic Hartree–Fock theory (HF), the relativistic Dirac–Hartree–Fock theory (DHF), and the density functional theory (DFT) [
24,
48,
49,
50,
51,
52,
53]. Interestingly, the FEA result can be simply reduced from IA result by replacing the Compton profile with the Dirac delta distribution
. Considering the atomic binding effects and electron pre-collision motions, the IA formulation could overcome the shortcomings in the FEA method. It is a practical approach to calculate the Compton scattering of bound electrons with X-rays and gamma-rays. Later research shows that the IA approach could serve as a good approximation in energy region near the “Compton peak” [
29,
42] (see
Figure 2 and
Figure 3).
Although the IA formulation effectively takes into consideration the atomic binding effects and electron pre-collision motions around atomic nuclei, it still has limitations in dealing with Compton scattering. In IA formulation, all the many-body effects in Compton scattering process are incorporated into Compton profile
, which is an observable only related to the atomic or molecular ground-states. The many-body effects coming from the ionized states and the interference terms among different momentum eigenstates are still insufficient in the dynamical process of Compton scattering. In the past few years, several approaches beyond the IA had already been investigated [
28,
29,
30,
31,
32,
33,
34]. This research, which mainly employs the low-energy theorem (LET) and scattering matrix (SM) approaches, has revealed many nontrivial properties of Compton scattering and has attracted lots of interests in interdisciplinary studies. Through comparing IA with these more advanced approaches, it is clear that the validity region for IA approach is just near the Compton peak region. Furthermore, in the validity region of IA, the momentum transfer
K in Compton scattering is much larger than the average momentum
for bound electrons (namely
) [
28,
32,
54,
55]. We will specialize in the SM approach in
Section 5.
4. Incoherent Scattering Factor/Incoherent Scattering Function
In this section, we give an introduction of the incoherent scattering function/incoherent scattering factor (ISF) approach. First, we can use the IA result as an example to illustrate the ISF approach. Starting from the IA result of DDCS in Equations (
15) and (
19), the angular distribution of Compton scattering can be calculated by an integration of final state photon energy
over the allowed range
After this integration, the angular distribution
can be reduced to the product of Klein-Nishina result in Equation (
9) and a correction factor [
21,
22].
The correction factor
is called the scattering function or scattering factor. Apart from IA, there are other ways to give the same results in Equation (
27) [
25,
26]. This kind of approach, in which the angular distribution of Compton scattering is given by Equation (
27) and
, is called the ISF approach.
We can use the following way to understand the physical meaning of scattering function
. If the atomic or molecular system has N electrons, and each electron scattered with photon independently as free electrons. Then the angular distribution of Compton scattering for whose system becomes
Therefore, the scatting function
, which is defined to be the ratio between the contribution of total system and that of single electron
can be view as the number of activated electrons in the Compton scattering process.
Here, we would present the scatting function
calculated within the RIA approach, the scatting function calculated using NRIA can be obtained in a similar way. To obtain scatting function
from RIA, we need to substitute the DDCS in Equation (
19) into Equation (
26) and evaluate the integral. After using some approximations, the factorization results of angular distribution in Equation (
27) can be obtained. Finally, the scattering function
can be expressed as [
21,
47,
56]:
The Heaviside step function
guarantees that only the activated electrons are included. In Compton scattering process, the electron becomes activated when the transferred energy
is larger than the binding energy
for
i-th subshell. In Equation (
30),
is the number of electrons in
i-th subshell,
denotes the maximum value of
for the
i-th subshell electron
and function
is defined to be an integral for Compton profile
Here, the function
is the single-electron Compton profile for
i-th subshell. It can be expressed as
The
is the electron momentum distribution for
i-th subshell, which can be calculated by ground state momentum wavefunctions. Then the total Compton profile for atomic or molecular system is given by
The quantity
in Equation (
30) is the scattering function in RIA formulation. In principle, it is a two-variable function depending on initial photon energy
and scattering angle
. However, in the nonrelativistic limit, these two variables are related to each other and they cannot be fully separated, which makes the scattering function
further reduce to a single-variable function [
27,
56]. The interdependence of
and
is realized via a new variable
x used in Hubbell’s work [
27]
Actually, in this expression, x is proportional to the momentum transfer in the elastic scattering between photon and electron, which is the Rayleigh scattering process, at scattering angle .
It should be mentioned that although the above ISF result on angular distribution displayed in Equation (
27) is derived from IA approach, the same result can be derived from other theories. In other words, there are alternative methods to calculate the scattering function
in Compton scattering process. For instance, in the nonrelativistic Waller–Hatree theory, the scattering function is given by [
25]
where
N is the number of electrons in atomic or molecular systems, and
is form factor
Here,
is the single-electron wavefunction for
i-th electron, and
is momentum transfer vector. Later, in a work presented by J. H. Hubbell et al., a more accurate expression for scattering function was used to include ionized state effects [
27]
where the summation is over all electrons in atomic or molecular systems. In this equation,
is the ground state for atomic or molecular systems,
and
are position vectors for
i-th and
j-th electrons correspondingly, and
is the total form factor for atoms or molecules defined by
The above results predicted by ISF approach are valid when variable
x and momentum transfer
K are large, which has been confirmed by experiments [
36,
39,
57]. However, there are limitations on the ISF approach. First, the ISF approach based on Equation (
27) can only be used to calculate the angular distribution for Compton scattering process. The information coming from the more differential quantities, such as DDCS
, are lost in the integration. Secondly, in the ab initio calculations of scattering function
, contributions from Compton scattering and Raman scattering cannot be efficiently distinguished [
29,
40]. In numerical tabulations, contributions of these two processes are summed over to give a total result. In the high energy region, it is lucky the Compton scattering is dominant compared to the Raman scattering.
The same as Compton profile
discussed in
Section 3, the scattering function
could also reflects the properties of target materials. Therefore, scattering function also creates a bridge between Compton scattering and interdisciplinary studies. It offers opportunities to learn the electronic structure and properties of target materials from Compton scattering process [
27,
57].
5. Scattering Matrix Approach
In this section, we specialize in the scattering matrix (SM) approach. In this approach, atomic binding effects and many-body effects in atomic or molecular systems are taken into account more precisely. It not only provides us a more complete understanding of the main features of Compton scattering with bound electrons, but also helps us recognize the validity regions for other methods (such as FEA and IA). With these advantages, the SM approach has attracted large interests in interdisciplinary studies [
32,
34].
In the SM approach, the DDCS of Compton scattering process can be calculated by the scattering matrix of a many-body theory
Therefore, it is known as the scattering matrix (SM) approach. In this approach, the scattering matrix of Compton scattering process
can be calculated through the many-body interaction Hamiltonian
for atomic or molecular systems. The interaction Hamiltonian
is determined by electromagnetic interactions between bound electrons and radiation photon fields, and it is given by the Quantum Electrodynamics (QED) for atomic and molecular systems. There are two categories of SM approaches: the nonrelativistic and relativistic theories. In the nonrelativistic theories, the interaction Hamiltonian
is expressed as
Here,
is quantized electromagnetic vector potential which describe radiation field of incoming photon acting on atomic or molecular system. The
is the momentum for
i-th electron, and the summation is over all electrons. In the relativistic theories, the interaction Hamiltonian is written by
where
is the conventional Dirac matrices for
i-th electrons.
The earlier works on SM approach were carried out based on the nonrelativistic Hamiltonian
, and they were restricted to
term (only the first term in the square bracket of Equation (
42) was included) [
15,
58,
59,
60]. The contributions from the second term
were accomplished by M. Gavrila et al. in the 1970s [
61,
62,
63,
64]. The full relativistic treatment was first attempted by I. B. Whittingham [
65,
66] and then developed by P. M. Bergstrom, T. Surić, R. H. Pratt et al. in the 1990s [
28,
29,
54,
67,
68]. In these works, the initial and final states (
and
) in atomic or molecular systems were calculated within the independent particle model (IPM), in which the single-electron states are solved through the unperturbed Hamiltonian
Here,
and
are the momentum and position of the
i-th electron. In the IPM method, the potential for
i-th electron
is given by the mean field of atomic nuclear potential and many-body interactions from other bound electrons. The IPM treatment may not be a perfect method in dueling with electron non-local exchange and correlations, and new treatments to tame electron non-local exchange in the framework of SM is still in development [
40]. With the initial and final state wavefunctions
,
, and interaction Hamiltonian
, the DDCS of Compton scattering process can be achieved from Equation (
41). There are a lot of techniques to calculate the matrix elements
in Equation (
41), which is beyond the scope of this work. In particular, in the nonrelativistic theories, the scattering matrix element
reduces to the Kramers–Heisenberg–Waller (KHW) matrix element [
69,
70,
71].
The spectrum of DDCS in Compton scattering obtained from SM approach has three main features. First, when final state photon energy approach to zero, namely in the
limit, an infrared rise emerges in the spectrum [
32]. This infrared rise is an example of infrared divergence behavior, which is a common feature of QED [
11,
72]. The same infrared rise behavior is also predicted in the low-energy theorem (LET) method [
73,
74] (Infrared behavior arises in the
limit, and the energy transfer
. In these cases, the energy of outgoing photons is tiny, and it makes these outgoing photons very difficult to observe experimentally. The Compton scattering process is happened as if the incident photon energy
is totally absorbed by atomic or molecular systems without producing notable influences, similar to the photoionization process/photoelectric effect). Secondly, some resonant peaks appear near the characteristic X-rays energies, which are the transition energies between different atomic or molecular states (When the resonant peaks emerge at X-ray characteristic energies, the electron undergoes a deexcitation process to lower energy intermediate states before it is ionized in Compton scattering process. The resonant peak energies correspond to the transition energies in these deexcitation subprocesses). This resonant behavior mostly come from Compton scattering with L and M shell electrons [
29]. At last, the “Compton peak” is reproduced in the vicinity of Compton energy
. In the SM results, the center of Compton peak
is slightly different from the Compton energy
, and the difference
leads to the “Compton defect” or “asymmetry of Compton profile” [
32,
40,
75] (There are other ways to define “asymmetry of Compton profile”. For instance, in reference [
75], the “asymmetry of Compton profile” is defined as
. Here,
K is the modulus of momentum transfer vector
in Compton scattering process, and
is the effective Compton profile extracted from SM results. Actually, from a detailed analysis of kinematical behavior in Compton scattering process, the following conclusions can be obtained. If the summit of Compton peak in the SM result is located at Compton energy
exactly, namely
, the extracted Compton profile
would reach its maximum value at
. In this case,
is axisymmetric around
for spherical symmetric atoms, the same as conventional Compton profile in Equation (
18), and there is no “asymmetry”. However, when the summit of Compton peak is shifted from Compton energy
,
would not produce the maximal value at
. In this case,
is not axisymmetric around
for spherical symmetric atoms, so there is asymmetry in
between positive and negative
values). When the modulus of momentum transfer
K in Compton scattering is sufficiently large such that inequalities
and
are satisfied, the “Compton defect” becomes extremely small and it can be neglected [
28,
29,
32,
54,
68,
75]. Here,
is the average momentum in atomic or molecular systems,
a is parameter defined to be
, and
is the fine-structure constant. In these cases, the RIA result does not present notable deviations from the SM result in the Compton peak region. The spectrum of DDCS predicted by SM result is illustrated in
Figure 3. To summarize, in the SM results, DDCS of Compton scattering exhibit infrared rise, the resonant peak, and the Compton peak. The three categories of peaks arise in different energy ranges. When we discuss the Compton peak, the final photon energy
should near Compton energy
, and the incident photon energy
. When the resonant peak is observed, final photon energy
is just near the X-ray characteristic energy of atoms or molecules, and initial photon energy
is usually much larger than this characteristic energy. To see the infrared rise, final photon energy
should goes to zero (namely
), while the incident energy
is not necessary to be very small.
It should be noted that the SM approach is still in development now [
32,
33,
34,
37,
38,
40,
76,
77,
78]. In some studies, new treatments and techniques are pursued to handle electron non-local exchange and correlation by methods beyond IPM [
37,
38,
76]. Other studies are devoted to more complex scattering process, for instance, Hopersky et al. investigated the Compton scattering and Rayleigh scattering of two X-ray photons [
77,
78]. Apart from the theoretical explorations, there are several experiments which have provided evidence to confirm the SM approach [
35,
36,
37,
38,
39,
40,
41]. Recently, Max Kircher et al. conducted a kinematically complete Compton scattering experiment using X-rays produced from accelerators with energy about 2.1 keV. By measuring the angular distribution of the scattered photon, the experimental observations present large deviations with the FEA results, but the experimental data are consistent with theoretical predictions from SM approach [
41]. This observation indicates that the SM approach is becoming a promising tool to duel with Compton scattering with bound electrons. Furthermore, the resonant peaks in Compton spectrum near the characteristic X-rays energies predicted by SM approach have also been confirmed by experiments [
35]. However, despite lots of attempts, the infrared rise behavior in Compton spectrum predicted by SM approach (which is also predicted by LET approach) has not been confirmed in experiments yet.
Although there is still inadequacy in duel with some many-body effects (such as the electron non-local exchange), the SM is the most advanced and accurate approach in the ab initio calculations of Compton scattering with bound electrons over recent years. First, in SM approach, the initial state , final state , and the dynamical process of Compton scattering are all treated by many-body QED theory of atomic or molecular systems. It is a fully quantum many-body approach, not just making simple corrections to the FEA results, as with RIA and ISF approaches. Second, the SM approach can take the many-body effects in atomic or molecular systems into account as much as possible. In this approach, atomic bindings, electron motions around atomic nuclei and electron many-body interactions are all considered in the starting points of theoretical treatments. Third, SM results can reflect all the main features in Compton scattering process: the infrared rise at low energy, the resonant peak at X-ray characteristic energy, and the broaden “Compton peak” near the Compton energy . With the aforementioned superiority, SM can make more accurate predictions in the entire region of the spectrum. Furthermore, many other approaches, such as FEA and IA, can be derived from the SM approach by making appropriate and simplified approximations. Since it was developed, SM have attracted lots of interests in atomic and molecular physics, and it may have great impacts in these areas in the near future.
6. Comparisons between Theoretical Calculations and Experimental Measurements
In the section, to have a better understanding of characteristics and limitations of the approaches described above, the comparisons between theoretical calculations and experimental measurements are provided. Limited by the scope of this work, only a few representative examples are presented. For more examples, the readers can resort to references [
28,
36,
57,
79,
80,
81,
82,
83,
84,
85].
For the angular distribution of Compton scattering process, the theoretical and experimental results of Fe and Cu are presented in
Figure 4. In this figure, the incident photon energy is
keV. We choose Fe and Cu elements as representative examples of elemental metal and ferromagnetic metal, respectively. Results of other elements exhibit a similar behavior, and they are not displayed in this figure. The readers could find more examples in references [
79,
80,
81,
82]. The experimental measurements from references [
80,
81] and the theoretical predictions of FEA and ISF approaches [
81,
82] are plotted in this figure. In reference [
82], the scattering function
is calculated using Equation (
30) in RIA framework [
81,
82] (Please note that in
Section 4, we have mentioned that the scattering function
can be obtained from RIA, Waller-Hartree theory, and other methods. In particular, within the RIA framework, the scattering function
is calculated through the integration of Compton profile
, see Equations (
30)–(
33).). The results in
Figure 4 indicate that ISF is a better approach than FEA in calculating the angular distribution for Compton scattering process. The ISF results successfully reproduce experimental measurements at all angles, while the FEA results can bring about large discrepancies in the small angle regions. In this region, both the energy transfer and momentum transfer in Compton scattering process are very small, so that atomic electrons cannot be viewed as free electrons anymore. In the low-energy and low-momentum transfer region (correspond to small scattering angle
), only a small percentage of bound electrons in atomic or molecular systems are activated in the Compton scattering process. Therefore, the FEA results, in which all electrons are treated as free and activated, would tremendously overestimate the angular distribution of Compton scattering process.
For the DDCS of Compton scattering process, the theoretical and experimental results are presented in
Figure 5. This figure gives the results of Cu, Sn and Pb elements. It is worth noting that for the DDCS, the FEA result becomes singular at the Compton energy
, while it becomes zero at other energies. This is due to the Dirac delta function
in Equation (
10). For this reason, it is not that valuable to plot the FEA result, what we need to focus are the IA and SM predictions. In this figure, subfigures (a), (b) and (c) correspond to the following cases: (a). Compton scattering for Pb atom at
keV and
; (b). Compton scattering for K-shell electrons of Sn atom at
keV and
; (c). Compton scattering for K-shell electrons of Cu atom at
keV and
. The subfigures (a) and (b) display the Compton peak region, while the subfigure (c) presents the spectrum outside the Compton peak. In the subfigure (a), the SM result is calculated by Bergstrom et al. [
54], the RIA result is given by Qiao et al. [
42], and experimental measurements are given by Rullhusen and Schumacher [
83]. In the subfigures (b) and (c), the theoretical SM results are given through combined works of Bergstrom et al. and Gavrila [
28,
63], the RIA result is given by Bergstrom et al. [
28], and the experimental measurements are given by Basavaraju et al. [
84] and Manninen et al. [
85], respectively. In subfigure (c), the LET results are calculated by Bergstrom et al. [
28]. From these comparisons, some conclusions can be drawn. First, for the DDCS of Compton scattering, the FEA approach becomes deficient and inconvenient, because of the singular behavior in the spectrum. Second, the SM and RIA results are similar in Compton peak region
. Both RIA and SM results are consistent with experimental observations in the Compton peak region, when some discrepancies are included. Third, the SM result is largely different from the RIA result outside the Compton peak, due to the infrared rise mentioned in
Section 5 (and possible resonant peaks near X-ray characteristic energies, which are not emerged in subfigure (c) (In the SM results, resonant peak behavior often appears in the cases of L and M shell electrons [
28,
29]. In
Figure 5, subfigure (c) corresponds to the Compton scattering with K-shell electrons, so the resonant peaks do not emerge.)). In regions far from the Compton peak, especially the infrared region where final photon energy
is very small, more experimental data are required to test the SM results.
Since RIA result is reliable and consistent with experimental observations in the Compton peak region
, we can safely use Compton profile
in IA approach to tackle Compton scattering in peak region. In recent decades, much research has emerged to study Compton profiles using the theoretical IA approach combined with experimental measurements near Compton peak [
5,
7,
9]. In experiments, the Compton profiles for atomic or molecular systems can be exacted from experimental data via equation
Here,
is the experimental measured DDCS of Compton scattering. Experimental studies for Compton profile have attracted huge interest in recent years [
5,
53,
86]. Furthermore, there are other quantities similar to the conventional Compton profile
discussed above. For example, if the incident photon beams are polarized, the differential cross-section is connected with the magnetic Compton profile
[
87,
88,
89,
90,
91], where
and
are the spin polarized electron momentum densities. Other kinds of Compton profiles, such as the directional Compton profile obtained by setting
z axis along different crystallographic axes, are also widely studied in recent years [
5,
92,
93]. Limited to the scope of the present work, we only focus ourselves on the conventional Compton profile
defined in Equation (
17). Other kinds of Compton profiles are not discussed in detail. More information on various kinds of Compton profiles can be found in references [
5,
87].
7. Database and Applications
In recent decades, Compton scattering had been extensively applied into many branches of science, including atomic [
29], molecular, condensed matter [
5,
6,
7,
8,
9], astrophysical [
94], nuclear and elementary particle physics [
47,
95,
96]. A lot of experimental and theoretical investigations concerning X-rays and gamma-rays cannot be carried on without the help of Compton scattering [
36,
79,
97]. As discussed in
Section 3 and
Section 4, Compton scattering is a powerful tool to study momentum distribution for bound electron in atomic, molecular and condensed matter systems [
3,
4,
5], both for theoretical and experimental studies. With the help of Compton profile and Compton scattering experiments, electron correlations, Fermi surfaces and band structures in materials can be investigated [
6,
7,
8,
9]. The Compton profiles are also closely related to the positron annihilation angular correlation spectra [
98]. Furthermore, the development of modern gamma-ray spectrometer and imaging devices also benefits a lot from the Compton scattering [
99,
100,
101,
102].
In recent years, many databases on Compton scattering have already been built up. The most common quantities in tabulations and databases are Compton profile
, incoherent scattering function
, and total cross-section
(In the theoretical calculations on scattering function
and total cross-section
, contributions from Compton scattering and Raman scattering cannot be fully distinguished [
29]. In numerical tabulations, the contributions from these two processes are summed over to give a total result. In the high energy region, the Compton scattering process is dominant compared with the Raman scattering.). For engineering or industrial applications, the data of differential cross sections
and
in Compton scattering can be easily achieved from the tabulation of Compton profile and scattering function (using Equations (
3) and (
4)). The most widely used database on Compton profile
over past years was given by F. Biggs et al. in 1975, in which a complete study on atomic Compton profile for elements (
) was presented [
24]. In Biggs’s work, the nonrelativistic Hartree–Fock theory was used to calculate Compton profiles for light elements (
) and the relativistic Dirac–Hartree–Fock theory was used to calculate Compton profiles for heavy elements (
). On the scattering function
, J. H. Hubbell et al. provided extensive and widely available tabulations for elements (
) based on Equation (
39) in 1975, with ground state wavefunctions calculated by several methods [
27]. Later, Kahane gave refined calculations using RIA approach and Dirac–Hartree–Fock ground state wavefunctions for all elements (
) in 1998 [
56]. For the total cross-section
, J. H. Hubbell et al. also gave tabulations for elements (
) in their early work in 1975 [
27]. The up-to-date tabulations on total cross-section
are provided by EPDL and NIST databases, which combine the theoretical and experimental data [
103,
104].