Self-Consistent Field Theory of Ionization Collisions

A self-consistent field (SCF) theory for treating scattering systems was formulated previously and is extended here to the ionization process, in which the continuum orbitals are made square integrable by an amputation procedure. The method is applied to the electron-hydrogen scattering system in the zero angular momentum coupling model, and the differential cross section is compared with the recent results obtained by several other approaches. It is shown that the amputated continuum functions provide an effective projection of the scattering equation. SCF continuum functions generated in the present formalism are used to analyze the effective charge approximation.


Introduction
For many years there has been a great deal of experimental and theoretical work devoted to the study of electron impact ionization processes and its inverse, the three body recombination.Such processes are of interest in many fields, such as astrophysics, plasma physics, short-wave laser development, fusion research and medical physics.The theoretical description of this process is difficult because the system involves at least two particles in the continua, which interact via long range Coulomb potentials.A self-consistent field (SCF) treatment is even more problematic, due to the non-square integrability of the continuum wave functions.Recently, much theoretical analyses focused on these difficulties have been made, both formally and numerically, and various models have been treated by accurate numerical or variational calculations [1]- [9].Despite recent theoretical advances [10] [11] the ionization problem still remains as one of the outstanding subjects of theoretical interest.
In view of the successful application of the Hartree-Fock approach to manyparticle bound state problems in the past, we have attempted a similar treatment for the scattering systems.The conventional theories which have been used in the past for the bound state problems [12] [13] [14] [15] are not quite applicable for scattering and ionization problems, mainly because there are one or more continuum wave functions involved in the description of the final states.Evidently, the continuum functions are not square-integrable i.e. not L2, and this fact complicates the mathematical as well as the computational aspect of the theory.In order to remove this fundamental difficulty, we have developed previously a new procedure [16] [17] that introduces an amputation of the scattering function at large distance, making it square-integrable.Note that, except for the phase shift, the tail of the continuum function carries minimal dynamic information.The use of the amputated functions thus allowed the development of a self-consistent procedure for the evaluation of the scattering amplitude in a systematic way.It is termed a generalized Hartree-Fock theory (GHF).The theory was tested for its validity and applicability [18] [19] as it was successfully applied to the positron-hydrogen and electron-hydrogen scattering systems, and further extended to the positron-helium and electron-helium scattering systems, where the exact target functions are not known.In fact both the scattering and target functions are simultaneously determined self-consistently by the theory.
In the present paper, we extend the GHF theory to the ionization problem, where at least two continuum functions are to be made square-integrable (L2) by the amputation procedure.As the GHF theory focuses on properly treating the continuum function, it is eminently suited for the ionization problem.As preliminary to a full treatment of many physically interesting cases, we study in this report a simple electron-hydrogen scattering in the zero angular momentum approximation.Although the angular momentum is uncoupled, this model problem still contains much of the difficulties of a three-body Coulomb problem, and it also serves as a non-trivial test of the new approach to the ionization problem.In Section 2, the GHF theory for the elastic scattering is summarized and extended to the ionization process.In Section 3, we present the details of ionization of hydrogen by electron impact.Our results are summarized in Section 4. In Section 5, we completely analyze the GHF wave functions and critically examine the effective charge approximation.

Generalized Hartee-Fock Approch
The GHF for the scattering systems is summarized below, and then extended to the ionization problem.The theory is developed closely paralleling the SCF approach for the bound state case, and the differences in the scattering case are emphasized.

Theory
We summarize the salient properties of the GHF procedure for the simple elastic scattering.First, for an N-electron bound system consisting of a valence electron and a neutral atomic core of N − 1 electrons, the core state configuration assumes in its lowest order approximation a simple anti-symmetrized product of single-particle orbitals, all of which are L2.Thus, in a slightly more convenient form for our purpose below, we have , 1, 2, , with , i j ij ϕ ϕ δ = for each orbital i ϕ , and where i − Ψ are the (N − 1) elec- tron clusters which may be assumed to be fully antisymmetrized within the clusters.A is the anti-symmetrization operator.The nonlinear set of coupled equations for the individual orbitals i ϕ is then obtained by the projection, as Evidently, all the integrals involved are finite.
For the scattering problem with one or more open channels, k ϕ represents the continuum orbitals.The total wave function then behaves asymptotically as where k − Ψ describes N − 1 electron target state without the k-th (continuum) particle, and S k and C k are the free sine and cosine waves asymptotically.In the case of single channel scattering, the reactance matrix K assumes a simple function tanδ, where δ is the scattering phase shift.The normalizability condition on the single-electron orbitals breaks down for the continuum function k ϕ , i.e.In Equation (2.4), E k is the scattering energy and T k is the kinetic energy operator for the k-th continuum particle.D k cuts off the tail of k ϕ , and as a result the AWF X's are now square-integrable.The crucial point to note here is that the amputation does not sacrifice any vital information, since the tail of a continuum function carries very little dynamical information; it mainly contains the phase shift and the wave number for that channel.Therefore, no essential physics is lost by the amputation, but we have gained the L2 property.Furthermore, we replace the strong asymptotic condition (2.3) and the exact k − Ψ by a weak asymptotic condition (WAC) where ( ) and where E -kt is the approximate (HF) cluster energy associated the approximate trial function kt − Ψ .This is an im- portant relaxation of the original boundary condition; Equation (2.5) not only sets a weak accuracy requirement on the cluster function kt − Ψ , which will eventually be improved with the addition of more configurations to the solution, but also allows a SCF determination of both the target and scattering orbitals.The generalized SCF equations for scattering are given by the projections, as where t GHF Ψ ≅ Ψ is of the form (2.3).Furthermore, for the j-th bound orbitals, 0, where -k-j denotes the N-2 particles that exclude both the k-th and j-th particles.This is the desired set of the generalized Hartree-Fock (GHF) equations.Equations (2.6) can be reduced to compact forms where k U and i U are the effective SCF potentials seen by the continuum and bound state functions, respectively.The Y's are the exchange terms involving integrals of ij V , and E dti are the effective energies for the bound orbitals φ i for electron i in the potential U i .Because of the non-L2 difficulty, U i with i ≠ k for the target orbitals have never been derived previously within the SCF approach.
Evidently, Equation (2.7b) is a well-defined equation for the bound state orbitals distorted by the continuum orbital.Also U are totally new result which displays for the first time the SCF potential for the target system during a collision.As a result, the SCF target functions become mildly energy dependent.
As usual, the multi-configuration mixing is realized by writing the total wave functions as ∑ , where the last sum is over the closed (or some open) channels in the case of single (or multiple) channel scattering.This generalization will be considered in future reports.
In the previous reports we showed that the GHF theory, with two key ansatze of WAC and AWF, is effective in treating the many-body scattering problem.
We have explicitly demonstrated [18] [19] that the WAC converges to the strong (exact) asymptotic condition as the total wave function is improved via multiconfiguration mixing.

Ionization by Electron Impact
The GHF for scattering summarized above can immediately be adapted to the ionization problem of interest here in a natural way, because the amputation can now be carried out on the two continuum functions.The final state wave function is written as where Coulomb potentials may not be completely screened for the residual ions.This point becomes clearer in the discussion of specified systems (See Section 4).

Ionization of Hydrogen by Electron Impact
We consider the non-relativistic electron scattering by atomic hydrogen scattering with total orbital angular momentum L = 0 and total angular spin S = 0, 1.In the center of mass frame the non-relativistic Schrödinger equation for the interaction of electrons with a hydrogen atom is where 1 r and 2 r are the coordinates of the atomic electron and the incident particle respectively, relative to the infinitely heavy proton considered to be at rest, E t is the total energy of the system, and In the present formalism, the generalized Hartree-Fock (GHF) approach is chosen in its lowest approximation to be of the form , for the singlet (+) and triplet (−) cases respectively.This form automatically defines our model, i.e.
12 1 1 r r > → .Note that its converse is not valid, i.e. the model does not necessarily lead to (3.2).In the following the explicit reference to l 1 = l 2 = 0 is omitted to simplify notation.The corresponding amputated wave functions are defined as where the operator where V i (r) is the electron-electron interaction potential and given by ( ) The non-local potential y i describes the electron exchange in which the 2 electrons interchange roles and is defined as The V i 's are the direct SCF potentials for the u i 's, which have never been possible to derive for the continuum electrons.
The transition amplitude T fi is given by f Ψ and Ψ represent the wave functions describing the whole system in its fully interacting final state and non-interacting initial state, respectively.The fully differential cross section is given by ( ) The initial state, which takes into account the incident and the bound electrons, will be written as a product of a plane wave and an exact wave function of the hydrogen atom: which eventually is taken in the zero angular momentum state.The final state

Results
We first choose the GHF wave function of the scattering system, as in the case of Hartree-Fock, in its lowest approximation: . We carried out calculations at total energies E t = 2.0 Ry and E t = 3.0 Ry, since the previous data are available for these energies for comparison [10] [20].Our results for the triplet  1 and compared with the exterior complex scaling [20] (ECS) method.The ECS method seems to be the most accurate procedure available for computing ionization cross section for the model.It uses a two-dimensional grid to solve for the outgoing scattered wave function without explicitly imposing the asymptotic boundary condition for the three charged particles.The coordinates are then scaled by a complex phase factor, beyond a certain radius where the tail of the Coulomb potential is ignored.As a result, the scattered wave function decays like a bound state wave function, which makes the asymptotic boundary condition simple to satisfy. Figure 1 contains additional data obtained by the Convergent Close Coupling (CCC) method [21].Agreement between our results and all the others is very good for the triplet case.This is especially surprising because the GHF approximation is considered here only in its lowest approximation.Apparently, the amputated wave functions carry much of the essential dynamical information contained in the scattering functions.Hence, we have been able to impose self consistency as a result of the amputation procedure which carefully gets rid of the long range tail represented by the asymptotic part of the scattering functions.Another feature of the GHF approach is that, by definition, our ionization differential cross sections are symmetric about the point

E ε =
. This property is absent in the close-coupling approach.In Figure 2 we display the interaction potential V GHF between the two continuum electrons.This is the first time the SCF potential is displayed in the ionization problem.It is then compared to the pure Coulombic potential V C . Figure 3 shows that the amputated wave function X which carries all the dynamics is of short range.It is presumably due to the character of the present model.Our results for the singlet case are reported in Table 2. Figure 4 shows that the singlet case is more problematic, as the existing theories are not all in agreement.Since the GHF is quite distinct from all the other approaches and based on the potentially powerful SCF approach, further analyses should help clarify the situation.

Analysis of the GHF Wave Functions
From the result of the GHF treatment given in the previous section, it is clear that the GHF functions in the single configuration approximation are already in nearly exact form for all values of E ε for S = 1 and for

E ε <
for the S = 0 state.For these regions of small E ε , no configuration mixing seems necessary.This is rather surprising, because the ionization process is presumably a more complicated process as compared to the other elastic and inelastic collisions.Although it was shown [16] formally that a "proper" projection of the scattering equation can yield the "exact" amplitude, it is in general difficult to construct such a projection.Apparently, the projections employed in the GHF in terms of the amputated scattering functions seem to be very effective.
Since we have nearly "exact" wave functions for low E ε , it is of interest to examine the wave functions generated by other approximation procedures, such as the plane waves, pure Coulomb waves, and the effective charge Coulomb E. Zerrad  waves.In particular, the effective charge approximation can be critically analyzed.In the conventional approach to ionization, the integrals for the "post" form of the ionization amplitude do not converge unless the so-called Peterkop condition is met.That is, for the two effective charges Z 1 and Z 2 , it is required that Thus, in the limit of k 1 → 0, we have Z 1 → 1. Table 3 and Table 4 illustrate the wave functions in the various approximations, as they are compared with the GHF functions.The effective charges used in the tables are adjusted to simulate Table 3. Wave functions generated for the S = 0 case with Z 1 = 0.1, Z 2 = 0.2 at 0.4 E ε = : B = plane waves, P = pure Coulomb waves, Z = effective charge Coulomb waves, and G:GHF.The "prior" form of the amplitude does not have the convergence difficulty and thus such a condition is not necessary.The pure Coulombic functions with Z = 1 also totally fail.
Another surprising result of the GHF treatment is that the k 1 component (<k 2 ) E. Zerrad of the GHF wave function turned out to be nearly identical to the plane wave, in spite of complicated scattering equations they satisfy, including the full exchange terms.This is contrary to the conventional physical picture of more distortion for lower energy scattering.Apparently, several components in the equation work to cancel each other's effects.It also strongly suggests that much of the earlier calculations all agree at small ε because the low energy component of the final wave functions are more or less that of free particle.Further, the higher energy component with k 2 is only partially distorted by the coulomb potential, with the effective charge of roughly Z 2 ≈ 0.4.
The triplet case is also of some interest, as all the calculations more or less agree with each other for all E ε .The reason for this may be that the scattering functions for the triplet state have a node at r 1 = r 2 , which places a stringent constraint on the wave functions.However, the magnitude of the wave functions at the peaks near r 2 ≈ r 1 /2 can still be quite different among the various theories.In Figures 5-8, we plotted the wave functions in the various approximations, at several sample points of r 1 (>r 2 ).It clearly shows that the various wave functions are all similar in the triplet case, but not for the singlet.
Eventually, configuration mixing is needed to improve the wave functions in the singlet case with 0.2 E ε > . Such mixing was carried out previously in the GHF context for the electron-hydrogen scattering system, by adding new channels and solving a set of coupled equations.Alternatively, a variational treatment may also be possible, as in the bound state case.
A part of our study in progress includes an extension of additional channels in the total wave function to test the convergence of the theory.We adopt a vartiational procedure to improve on the ionization differential cross section.This

Conclusion
The present study shows that the GHF approach is a viable method to treat complex collision systems.The effectiveness of the theory may be attributed to the crucial physics information carried by the amputated wave function.As emphasized in [16], it is possible in principle to obtain accurate solutions to the collision problem by adopting a "proper projection" of the scattering equations, although it is in general difficult to find the proper projection.However, the GHF seems to provide a partial answer to the projection problem, in terms of the amputated functions X.For application of GHF to more realistic and complex systems, it is important to better understand the role played by the X's.The multiconfiguration extension of the GHF is in progress to clarify the situation.
where D k and D p may be different, not only in the constant kinetic energies, but also possibly due to the presence of long-range interaction.In some cases, the E. Zerrad DOI: 10.4236/jmp.2017.890911525 Journal of Modern Physics to be evaluated in the SCF approximation.The continuum wave functions are normalized to a δ function in momentum space, as in(2.3).As in the Hartree-Fock case, the resulting system of mutually coupled SCF equations is solved by iteration, subject to the self consistency requirement.First one guesses a set of scattering wave functions.Then one computes the SCF potentials and the inhomogeneous terms.The set of Equation (3.4) is then solved.These new scattering wave functions are amputated and used to compute the new SCF potentials and the new inhomogeneous terms, which in turn are used to obtain a third set of scattering wave functions.The cycle continues until self consistency is reached within some preset tolerance.The numerical method we adopt for the scattering functions is the use of the distorted Green's functions.In this method, we transform the set of differential equations that the scattering functions satisfy, along with their corresponding boundary conditions, into integral equations.The final step in the numerical procedure is to calculate the proper integral expressions for the scattering amplitudes T fi and the corresponding ionization cross section.

Figure 1 .
Figure 1.The single-differential ionization cross sections (in units of

Figure 2 .
Figure 2. Comparison of the Coulombic potential (cross +) and the GHF potential (cross x) for the continuum electron for the S = 1 case at 0.5 E ε = .

Figure 3 .
Figure 3. Plot of the continuum wave function u obtained by the GHF procedure and its corresponding amputated wave function X, for the S = 1 case at 0.5 E ε =

Figure 4 .
Figure 4.The single-differential ionization cross sections (in units of

Table 2
Ry and 3 Ry are given as a function of the energy fraction E εfor the S = 0 case.The Exterior Complex Scaling is denoted by ECS and the Generalized Hartree-Fock approach is denoted by GHF.v.

Table 4 .
Wave functions generated for the S = 0 case with Z 1 = 0.05, Z 2 = 0.4 at