Advances in Physical Chemistry

Volume 2011, Article ID 541375, 14 pages

http://dx.doi.org/10.1155/2011/541375

## Invariant Graphical Method for Electron-Atom Scattering Coupled-Channel Equations

^{1}School of Physics, The University of Western Australia, WA 6009, Australia^{2}Institute of Theoretical Physics, Curtin University of Technology, WA 6102, Australia

Received 9 September 2010; Revised 28 December 2010; Accepted 13 January 2011

Academic Editor: Wybren Buma

Copyright © 2011 J. B. Wang and A. T. Stelbovics. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We present application examples of a graphical method for the efficient construction of potential matrix elements in quantum physics or quantum chemistry. The simplicity and power of this method are illustrated through several examples. In particular, a complete set of potential matrix elements for electron-lithium scattering are derived for the first time using this method, which removes the frozen core approximation adopted by previous studies. This method can be readily adapted to study other many-body quantum systems.

#### 1. Introduction

The fundamental problem encountered by a theoretician working in the field of quantum physics or quantum chemistry is to evaluate the so-called observables to certain precision. Here, is the wave function of the system under investigation, and is the operator representing the physical quantity one needs to evaluate. This seemingly easy task encounters many obstacles, one of which is the complexity of the physical structure of a composite system. For example, there are several electrons and nuclei interacting with each other in an atomic or molecular system. Even though the exact form of such interactions is well known, it is exceedingly difficult to obtain accurate solutions from the corresponding Schrödinger or Dirac equation.

However, most physical systems have symmetry properties that influence their dynamics and structures in a way independent of the detailed interactions. A good example is the motion of a particle in a central field. The symmetry of the interaction leads to the conservation of angular momentum that is independent of the detailed form of the interaction potential . Angular momentum theory deals with rotational symmetries of a composite system. In classical physics, angular momentum treatments are straightforward following the simple manipulation rules of vectors. In the quantum world, not only the magnitude of angular momenta is quantized, but also their spatial orientations.

Angular momentum calculations in quantum theory prove to be most laborious, especially when many-body systems or complicated operators are involved. Furthermore, the final expressions are often incomprehensible even for experienced researchers in the field. Unless one derives the equations himself by a painstakingly tedious procedure, one cannot really tell where the phases and the - symbols come from. It is also quickly apparent, once one begins to derive the complete reduced potential matrix-elements for systems beyond four electrons, that the expressions become rapidly intractable to algebraic derivation by conventional methods. The problems are simply due to the vast number of Clebsch-Gordan coefficients that comprise the expressions, each of which is summed over several angular momentum labels. In addition, there are many pages of practical formulas that one often needs in order to perform these calculations, such as those listed in the appendixes of Brink and Satchler [1].

In view of these difficulties, we started to search for a method that is transparent, free of the incomprehensible phases, makes minimal use of formulas, and leads easily to a correct final expression. We will demonstrate in this paper that the invariant graphical method initiated by Danos and Fano [2] meets all the above criteria. We will then apply this method to derive the reduced matrix elements of selected electron-atom interaction potentials, which are not available mainly due to the difficulties mentioned above. This method can be readily adapted to study other many-body systems such as two-photon double ionization of atoms [3, 4], two-nucleon knockout reactions [5], spin networks [6], and entangled spin-1/2 qubits [7, 8].

A central point to this method is the concept that most physical quantities are invariant, that is, they are independent of the coordinate system one chooses for their description. This is the case for energies, cross-sections, and transition probabilities. In other words, these quantities are scalar quantities and have zero angular momentum. Consequently, the tensors they may contain would couple to a total of zero angular momentum. We will demonstrate, in the next section, that angular momentum calculations become much simpler once we take this invariance into account. Atomic units are used throughout this paper.

#### 2. Elements of the Invariant Graphical Method

In the invariant graphical method, there is only one basic graph (Figure 1(a)) for angular momentum coupling and recoupling. The essential elements are the horizontal lines, each representing a tensor with a given angular momentum labeled above the line. Here, tensors are used to represent all quantities including wave functions, operators, amplitudes, and so forth. In this way, a scalar quantity is a tensor of rank zero. A vector is a tensor of rank one and so on. As an example, the state with angular momentum can be represented by a tensor of the rank . This graph provides the basic transformation associated with four angular momenta to rearrange them into a different coupling scheme. The transformation is represented by a recoupling box, mathematically associated with a square 9- symbol such that where is represented by . The square 9- symbol is related to the Wigner 9- coefficients by where .

This basic graph is the key building block of this method. All possible recoupling transformation can be performed based upon this graph. For example, the recoupling of three angular momenta can be derived using the graph illustrated in Figure 1(b) as In this method, there is no need to use 6- symbols, which carry extra phases and normalization constants that can lead to tedious book-keeping and are prone to errors. There is also no need to remember the many sum rules, which are required to eliminate unnecessary summation indexes in other conventional methods.

#### 3. Applications

Coupled-channel equations have been used extensively in the description of electron-atom scattering. Their application has been particularly successful in the so-called convergent close-coupling approach (CCC) implemented by Bray and Stelbovics [9] and subsequently extended to other atoms with one or two valence electrons outside an inert core [10, 11, 13]. To perform such calculations, one requires the reduced matrix elements of the coupled-channel potentials, which we denote by with representing the complete set of quantum numbers specifying a particular configuration of the system under study. For electron scattering from an -electron target, the coupled-channel reduced matrix elements have the following structure which includes both the direct and exchange terms. Here, is the permutation operator, is the system energy, subscript 0 denotes the incoming electron, and subscripts 1- denote the target electrons. Using the explicitly antisymmetrized target states, the reduced matrix elements can be simplified to Each of these terms will be derived explicitly for e-He and e-Li scattering in Sections 3.2 and 3.3.

By employing the invariant graphical method, we will derive a complete set of direct and exchange potential matrix elements for electron-atom scattering, which are beyond the frozen core approximation adopted by Wu et al. [12], Fursa and Bray [13], Bray et al. [10, 16], and Zhang et al. [14].

##### 3.1. Two-Electron System

To illustrate the use of the invariant graphical method, we first consider the direct Coulomb matrix elements for a two-electron system, that is, , which was also discussed in detail in Danos and Fano [2]. The complete coupling and recoupling graph is shown in Figure 2. As described earlier, each line in the figure represents a tensor with a given angular momentum indicated by the label above the line. The joining of two lines into one line represents coupling of their angular momenta. Recoupling into a different tensor set is indicated by a transformation box and needs to be carried out until each component of the graph is expressible in terms of single-particle matrix elements (e.g., the end boxes and in this case). The invariance of such a triple product provides a means of eliminating unnecessary intermediate indexes. For example, it demands that the coupling of and must give rise to a tensor with angular momentum instead of an arbitrary index. It also helps in identifying the most direct and economical recoupling scheme. We first expand the two-body Coulomb potential into its multipole components where . The corresponding final expression can then be directly read off from Figure 2 as where and is the regular spherical harmonic tensor. This definition is introduced to avoid unnecessary phases entering the coupling and recoupling scheme.

Each square 9- symbol in (7) corresponds to a recoupling box in Figure 2. The end-box represents an invariant matrix element containing three spherical harmonic tensors written as , which equals For numerical computation purposes, (7) can be treated as a final expression. Nevertheless, further simplification can be made considering the many zeros in the 9- symbols. This leads to the well-known result (see, e.g, [15])

##### 3.2. Electron-He Scattering (A Three-Electron System)

As a second more complicated example, we will derive the momentum-space direct and exchange matrix elements for a three-electron system, such as electron scattering from a helium atom. The general configuration of the system can be represented by , where , , and are, respectively, the linear, orbital angular and spin momentum of the incoming electron, and are the angular momenta of the target electrons, and are the spin momenta of the target electrons, and are the total orbit angular and spin momentum of the target atom, and and are the total orbit angular and spin momentum of the entire system. Note that the linear radial integral part of the direct and exchange elements is straightforward to work out and has the general form of where are the target atomic orbitals with quantum numbers and , and is the th partial wave of the projectile electron. The coupling and recoupling graph for the direct and exchange potential spin part is shown in Figure 3 and those for the orbital part are shown in Figures 4, 5, 6, 7, and 8. The corresponding final expressions can then be directly read off from the figures as Minor simplification of (11)–(15) gives rise to the same expressions presented in [11], where the conventional algebraic approach was used. However, we emphasize that the invariant graphical method is so much simpler with the derivation completed by drawing the diagrams.

##### 3.3. Electron-Li Scattering (A Four-Electron System)

Derivations for these momentum-space potentials have been previously given for e-H [9], for e-Li [11], for e-Na [16], and for e-He [16], all using the conventional algebraic approach. The frozen core approximation was adopted for e-Li and e-Na in these studies, which could be the main reason for the small discrepancies between the theoretical calculations and experiments [10]. With the enormous success of quantum scattering theories in describing scattering from one- and two-electron targets, one is naturally seeking to perform calculations with more complex systems and free of approximations. However, such an extension has been proven to be extremely tedious and in some cases even intractable. In the following, we derive a complete set of direct and exchange matrix elements for electron-Lithium scattering (previously not available) using the invariant graphical method. In this way, we are able to remove the frozen core approximation adopted by Bray et al. [10, 16] and Zhang et al. [14]. Again, the linear radial integral part of the direct and exchange elements is straightforward to work out and has the general form of The corresponding coupling and recoupling graphs are shown in Figures 9, 10, 11, 12, 13, 14, and 15, and the results are as the following:

#### 4. Conclusions

The invariant graphical method developed by Danos and Fano [2] is very powerful. The complex manipulation of angular momenta is reduced to the drawing of a compact graph, from which the final expressions can be readily read off. The procedure is transparent and simple to apply. It avoids writing out the intermediate expansions, yields immediately the selection rules for intermediate angular momenta, and helps in finding the most direct and economical intermediate recoupling. It also requires a minimal use of rules in comparison with pages of formulas involved in other conventional methods.

#### Acknowledgments

J. B. Wang would like to thank H. Y. Wu for cross-checking the phases in (11) using the conventional algebraic method. The authors also acknowledge support from The University of Western Australia, Murdoch University, and Curtin University of Technology.

#### References

- D. M. Brink and G. R. Satchler,
*Angular Momentum*, Oxford University Press, Oxford, UK, 1993. - M. Danos and U. Fano, “Graphical analysis of angular momentum for collision products,”
*Physics Report*, vol. 304, no. 4, pp. 155–227, 1998. View at Google Scholar - A. S. Kheifets, “Sequential two-photon double ionization of noble gas atoms,”
*Journal of Physics B*, vol. 40, no. 22, pp. F313–F318, 2007. View at Publisher · View at Google Scholar - A. S. Kheifets, “Photoelectron angular correlation pattern in sequential two-photon double ionization of neon,”
*Journal of Physics B*, vol. 42, no. 13, Article ID 134016, 2009. View at Publisher · View at Google Scholar - E. C. Simpson, J. A. Tostevin, D. Bazin, and A. Gade, “Longitudinal momentum distributions of the reaction residues following fast two-nucleon knockout reactions,”
*Physical Review C*, vol. 79, no. 6, Article ID 064621, 2009. View at Publisher · View at Google Scholar - V. Aquilanti, A. C. P. Bitencourt, C. Da, A. Marzuoli, and M. Ragni, “Combinatorics of angular momentum recoupling theory: spin networks, their asymptotics and applications,”
*Theoretical Chemistry Accounts*, vol. 123, no. 3-4, pp. 237–247, 2009. View at Publisher · View at Google Scholar - J. Suzuki, G. N. M. Tabia, and B. G. Englert, “Symmetric construction of reference-frame-free qudits,”
*Physical Review A*, vol. 78, no. 5, Article ID 052328, 2008. View at Publisher · View at Google Scholar - A. Maser, U. Schilling, T. Bastin, E. Solano, C. Thiel, and J. Von Zanthier, “Generation of total angular momentum eigenstates in remote qubits,”
*Physical Review A*, vol. 79, no. 3, Article ID 033833, 2009. View at Publisher · View at Google Scholar - I. Bray and A. T. Stelbovics, “Convergent close-coupling calculations of electron-hydrogen scattering,”
*Physical Review A*, vol. 46, no. 11, pp. 6995–7011, 1992. View at Publisher · View at Google Scholar - I. Bray, J. Beck, and C. Plottke, “Spin-resolved electron-impact ionization of lithium,”
*Journal of Physics B*, vol. 32, no. 17, pp. 4309–4320, 1999. View at Publisher · View at Google Scholar - D. V. Fursa and I. Bray, “Calculation of electron-helium scattering,”
*Physical Review A*, vol. 52, no. 2, pp. 1279–1297, 1995. View at Publisher · View at Google Scholar - H. Wu, I. Bray, D. V. Fursa, and A. T. Stelbovics, “Convergent close-coupling calculations of positron-helium scattering at intermediate to high energies,”
*Journal of Physics B*, vol. 37, no. 6, pp. 1165–1172, 2004. View at Google Scholar - D. V. Fursa and I. Bray, “Electron-impact ionization of the helium metastable 2 S state,”
*Journal of Physics B*, vol. 36, no. 8, pp. 1663–1671, 2003. View at Publisher · View at Google Scholar - X. Zhang, C. T. Whelan, and H. R. J. Walters, “Electron impact ionization of lithium-spin asymmetry of the triple differential cross section,”
*Journal of Physics B*, vol. 25, pp. L457–L462, 1992. View at Google Scholar - R. D. Cowan,
*The Theory of Atomic Structure and Spectra*, University of California Press, Berkeley, Calif, USA, 1981. - I. Bray, “Convergent close-coupling method for the calculation of electron scattering on hydrogenlike targets,”
*Physical Review A*, vol. 49, no. 2, pp. 1066–1082, 1994. View at Publisher · View at Google Scholar