## Journal of the Brazilian Chemical Society

##
*Print version* ISSN 0103-5053*On-line version* ISSN 1678-4790

### J. Braz. Chem. Soc. vol.17 no.1 São Paulo Jan./Feb. 2006

#### http://dx.doi.org/10.1590/S0103-50532006000100023

**ARTICLE**

**A theoretical study on the photoionization of the valence orbitals of phosphine**

**Edmar M. Nascimento ^{I}; Luiz E. Machado^{II, }^{*}; Evandro M.S. Ribeiro^{III}; Luiz M. Brescansin^{IV}; Mu-Tao Lee^{V}**

^{I}Instituto de Física, Universidade Federal da Bahia, 40210-340 Salvador-BA, Brazil

^{II}Departamento de Física, Universidade Federal de São Carlos, 13565-905 São Carlos-SP, Brazil

^{III}Instituto de Física de São Carlos, Universidade de São Paulo, CP 369 13560-970 São Carlos-SP, Brazil

^{IV}Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, 13083-970 Campinas-SP, Brazil

^{V}Departamento de Química, Universidade Federal de São Carlos, 13565-905 São Carlos-SP, Brazil

**ABSTRACT**

We report a theoretical study on the photoionization of phosphine in the static-exchange level and frozen core approximation, using the method of continued fractions. The main subject of the present study is to investigate in which extent the Hartree-Fock description of the target applied to molecular photoionization is valid. Also, the role played by multichannel coupling is analysed. Our study shows that single-channel Hartree-Fock calculations can provide reliable results except for photon energies near the photoionization threshold.

**Keywords:** photoionization cross section, photoionization asymmetry parameters, MCF, PH_{3}, phosphine

**RESUMO**

Apresentamos um trabalho teórico sobre a fotoionização da fosfina ao nível estático-troca e da aproximação de caroço congelado, usando o método de frações continuadas. O objetivo principal do presente estudo é investigar até que ponto a descrição Hartree-Fock do alvo é válida, quando aplicada em estudos de fotoionização molecular. Também analisamos o papel desempenhado pelos acoplamentos multicanais. Nosso estudo mostra que cálculos Hartree-Fock monocanais podem fornecer resultados confiáveis, exceto para energias do fóton próximas ao limiar da fotoionização.

**Introduction**

Phosphine is a colorless, flamable, and toxic gas. This compound has important applications in several fields. For instance, in semiconductor technology, phosphine is used in silicon processing as a source for phosphorous implantation.^{1} Also in agriculture, it has been widely used as a pesticide.^{2,3} Particularly, for the latter application, photodecomposition is a very important mechanism for pesticide degradation and therefore studies on photoionization of this chemical species can contribute to environmental protection. Despite that, both experimental and/or theoretical studies on photoionization of phosphine are very limited. To our knowledge, there is only one experimental measurement of photoabsorption and photoionization cross sections for this molecule reported by Zarate *et al*.^{4} Experimental measurements of asymmetry parameters for photoionization from the valence orbitals of this molecule have been reported by Cauletti *et al*.^{5} Theoretically, one of the pioneer studies on photoionization of PH_{3} was performed in 1998 by Stener and Decleva^{6} using a B-spline calculation based on the density functional theory (DFT) approach. Lately, these authors have improved their method by using a time-dependent density functional theory (TD-DFT) approach^{7,8} in which an explicit treatment of the unbound continuum wave function and a linear response of the system to the external weak electromagnetic field are considered. Also, an *ab initio* study using a combination of the K-matrix (KM) technique and the random-phase approximation (RPA) was performed by Carravetta and Cacelli.^{9} In this latter study, the KM-RPA method is used as a tool to solve a configuration interaction problem in the continuum through the resolution of the RPA equations, projected on an L^{2} basis set, by the K-matrix employed. Therefore these methods account, in some extent, for the electron correlation effects which may be important in molecular photoionization calculations. In addition, the KM-RPA can also describe the coupling of different ionization channels that can be opened in a range of a few eV.

Recently, we have developed a computational code for molecular photoionization calculations based on the method of continued fractions (MCF).^{10} In our method, both the ground- and the continuum-state wave functions are described by the Hartree-Fock (HF) approximation. As a first application, the HF-MCF was applied to study the photoionization of the outer-valence orbitals of ammonia.^{10} The comparison of our calculated cross sections and asymmetry parameters with the existent theoretical and experimental results was very encouraging. In the present work, we apply the HF-MCF to study the photoionization of the outer-valence orbitals of phosphine at the static-exchange level and frozen core approximation. The comparison of our calculated data with the experimental results can provide information on which extent the HF description of the target applied to molecular photoionization studies is valid. Also, the comparison with the calculated data obtained using a more elaborate theory would provide an insight on the role played by the multichannel-coupling and electron-correlation effects, which can be helpful for improving our method.

**Theory and Calculations**

Details of the MCF have been given elsewhere^{10-14} so only the essential aspects of the theory will be discussed here. The photoelectron differential cross sections averaged over molecular orientations are given by:

where s^{(L,V)} is the total photoionization cross section, obtained with the length (*L*) or velocity (*V*) form of the dipole moment operator. For nonlinear molecules, s^{(L,V)} is given by:

The asymmetry parameter appearing in equation (1) can be written as:

The quantity *E* in equation (2) is the incident photon energy, whereas * *in equations (2) and (3) are the partial-wave components of the dynamical coefficients,

where

and

In equations (5) and (6) F* _{i}* is the target ground state wave function, the final state (incoming-wave normalized) wave function of the system (ion plus photoelectron), represents the unit vector in the direction of polarization of the radiation and is the photoelectron momentum.

In the present study, F* _{i}* is an one-determinant wave function calculated at the Hartree-Fock level. The final molecular state is described by a single electronic configuration in which the ionic orbitals are constrained to be identical to those of the initial ground state. In this approximation, the photoelectron orbital is a so-lution of the one-electron Schrödinger equation

where satisfies the appropriate boundary conditions and *V _{N-1}* is the static-exchange (SE) potential of the molecular ion with the long-range Coulomb potential of the ionic core removed.

^{10}

To proceed, equation (7) is rewritten in an integral form, the Lippmann-Schwinger equation, as follows:

and then is solved using the MCF. In equation (8) is the state vector corresponding to the principal-value of the photoelectron wave function, * *is the principal value coulombic Green's operator, _{} is a state vector corresponding to the principal value of the coulombic wave function and *U*^{(0)} = *2V _{N}*

_{-1}is the reduced potential operator. Accordingly, the reactance

**K**matrix is given by

The MCF projects the interaction potential on a finite set of functions generated at each step of an iterative procedure and requires no basis functions for the solution of the continuum wave function. So we define the *n*th-order weakened potential operator *U ^{(n)}* as

Therefore, the reactance **K** matrix results in the form of a continued fraction. The *n ^{th}*-order correction of the matrix is given by

where

and

The reactance matrix is given by

where

with

and

It is expected the operator *U ^{(n)}* to become weaker and weaker with increasing

*n*. As a result, the iterative procedure can be stopped after some number of steps, when the desired convergence is achieved. The converged

**K**matrix corresponds to the exact solution for a given potential

*U*

^{(0)}in equation (8).

**K**

_{1}can be obtained backwardly using equation (11) by setting

**K**

_{n+}_{1}= 0. Then, the

**D**matrix is obtained using equation (16). In order to constrain the continuum solution to be orthogonal to the occupied bound orbital (f

*), the short-range portion of the static-exchange potential is calculated using a generalized Phillips-Kleinman pseudopotential defined by*

_{i}^{15}

where *V _{SE}* is the usual static-exchange potential and

In equation (20) e is the photoelectron energy.

**Results and Discussion**

The ground state of neutral PH_{3} belongs to C_{3v} symmetry and has the electronic configuration

The corresponding SCF-HF wavefunction was obtained using the basis set shown in Table 1. At the equilibrium geometry this basis set gives an SCF energy of -342.46998 a.u., which can be compared with the SCF energy of -342.36424 a.u. of Carravetta and Cacelli^{9} and with the RHF value of -342.4934 of Clark *et al*.^{16} Our calculated dipole moment was 0.541 D, to be compared with the calculated value of 0.624 D of Clark *et al*.^{16} and with the experimental value of 0.578 D of Burrus.^{17} The static potential was evaluated using a single-center expansion truncated at *l _{max}* = 20, and all matrix elements arising in the MCF were calculated with

*l*= 6, retaining all allowed |

_{max}*m*|

__<__l for a given

*l*. In all of our calculations convergence was achieved within five iteractions.

Firstly, we intend to compare the results of our calculations at the Hartree-Fock (HF) and single-channel (SC) level with the available experimental results as well as with other reported calculations at the same SC level. Figures 1(a) and 2(a) show our calculated results of s for the photoionization out of the two outermost orbitals 5a_{1} (IP = 10.6 eV) and 2e (IP = 13.6 eV), respectively. The results are presented in both dipole-length (DL) and dipole-velocity (DV) forms, along with theoretical DL results obtained by Stener and Decleva^{6} using the B-spline technique and with those obtained by Carravetta and Cacelli^{9} using the SC-RPA. The experimental results of Zarate *et al*^{4} are also included for comparison. The observed differences between our DL and DV results can be attributed to the neglect of the electronic correlation in the target wavefunction.^{9,18} Our HF-SC results do not reproduce the sharp fall-off of s near the threshold seen in both the experimental and SC-RPA data. Probably, this is also due to that neglect. Away from this region, it is observed that the experimental results for s lay between our DL and DV results, in such a way that for the channel our DL results agree better with the experimental data, whereas for the 2e^{-1} channel our DV results do so. The DL results of the B-spline DFT calculation of Stener and Decleva^{6} lay well below the experimental results for both channels.

Comparison of our HF-SC theoretical results for b with experimental results of Cauletti *et al* ^{5} and with the theoretical ^{6,9} data is shown in the Figures 1(b) and 2(b), for photoionization out of the 5a_{1} and 2e orbitals, respectively. Both our DL and DV results are in reasonable good agreement with the experimental results^{5} and with the SC-RPA data of Carravetta and Cacelli ^{9} in the entire energy range, whereas the B-spline DFT results of Stener and Decleva^{6} underestimate the experimental data in the low-energy region, for both channels. Also, the positions of the maximum and minimum structures are generally well reproduced by our results, whereas the DL results of Stener and Decleva^{6} show minima that are shifted to lower energies, mainly for the photoionization out of the 2e orbital.

In order to estimate the effects of multi-channel (MC) coupling, in Figures 3 and 4 we compare our DL results for the 5a_{1} and 2e channels, respectively, with the corresponding MC-RPA results of Carravetta and Cacelli^{9}, as well as with the experimental results of Zarate *et al.*^{4} and Cauletti *et al*.^{5} The calculated DL data using the TD-DFT approach^{7,8} are also shown in these figures for comparison. For both s and b quantities, our results for energies above 12 eV (for ) and 15 eV (for 2e^{-1}) reproduce fairly well the multi-channel results of Carravetta and Cacelli.^{9} On the other hand, at lower energies, the MC-RPA s results of Carravetta and Cacelli also show the correct near-threshold behavior, thus confirming that such a behavior can only be described by the inclusion of the electronic correlation effects in the target wavefunction.^{9} In contrast, the MC and the SC results of Carravetta and Cacelli^{9} for b are not substantially different, except for the occurence of some structures that are seen in the MC results which can be attributed to the inter-channel coupling. In addition, the deficiency observed in the Kohn-Sham DFT results, seen in Figures 1 and 2, is completely recovered by the TD-DFT approach^{7,8} in such a way that the TD-DFT results for both s and b agree much better with the experimental data^{4} and with the RPA results of Carravetta and Cacelli^{9}. The response of the system to the external time-varying field has been attributed by the authors^{7,8} as the physical origin of this improvement. This statement, however, may not be true since our HF-SC results, although not accounting for such a response, are already in good agreement with both the TD-DFT and the experimental results. Therefore, we believe that the explicit treatment of the unbound continuum wave function may be responsible for that improvement. Moreover, the structures seen in the TD-DFT results for both s and b, due to the occurrence of Feshbach one-electron autoionization resonances, suggest that the TD-DFT calculations are of multichannel nature.

In summary, we have shown that our HF-SC calculation can provide reliable results for both s and b in the energy range from a few eV above threshold to intermediate energies. However, the near-threshold fall-off of the s values can only be well described in the framework of approximations beyond the HF level. On the other hand, both s and b are only slightly influenced by a multi-channel treatment. Considering the numerical simplicity of our method, it can be considered an adequate tool for the treatment of valence-shell molecular photoionization, specially for heavier molecules. A more complete treatment for such systems would constitute a formidable computational task.

**References**

1. Adler, J.; Doering, H.R.; Grosse, H.J.; Gleisberg, F.; *Zfl. Mitt.* **1983**, *69*, 328. [ Links ]

2. Scudamore, K.A.; Goodship, G.; *Pestic. Sci.* **1986**, *17*, 385 [ Links ]

3. Leitao, J.; De Saint Blanquat, G.; Bailly, J.R.; *Appl. Environ. Microbiol.* **1987**, *53*, 2328. [ Links ]

4. Zarate, E.B.; Cooper, G.; Brion, C.E.; *Chem. Phys.* **1990**, *148*, 277. [ Links ]

5. Cauletti, C.; Adam, M.Y.; Piancastelli, M.N.; *J. Electron Spectrosc. Relat. Phenom.* **1989**, *48*, 379. [ Links ]

6. Stener, M.; Decleva, P.; *J. Electron Spectrosc. Relat. Phenom.* **1998**, *94*, 195. [ Links ]

7. Stener, M.; Decleva, P.; *J. Chem. Phys.* **2000**, *112*, 10871. [ Links ]

8. Fronzoni, G.; Stener, M.; Decleva, P.; *J. Chem. Phys.* **2003**, *118*, 10051. [ Links ]

9. Carravetta, V.; Cacelli, I.; *Chem. Phys.* **1999**, *243*, 77. [ Links ]

10. Nascimento, E.M.; Ribeiro, E.M.S.; Brescansin, L.M.; Lee, M.-T.; Machado, L.E.; *J. Phys. B* **2003**, *36*, 3621. [ Links ]

11. Horacek, J.; Sasakawa, T.; *Phys. Rev. A* **1983**, *28*, 2151. [ Links ]

12. Horacek, J.; Sasakawa, T.; *Phys. Rev. A* **1984**, *30*, 2274. [ Links ]

13. Lee, M.-T.; Iga, I.; Fujimoto, M.M.; Lara, O.; *J. Phys. B* **1995**, *28*, L299. [ Links ]

14. Ribeiro, E.M.S.; Machado, L.E.; Lee, M.-T.; Brescansin, L.M.; *Comp. Phys. Comm.* **2001**, *136*, 117. [ Links ]

15. Lucchese, R.R.; Raseev, G.; McKoy, V.; *Phys. Rev. A* **1982**, *25* 2572. [ Links ]

16. Clark, S.A.C.; Brion, C.E.; Davidson, E.R.; Boyle, C.; *Chem. Phys.* **1989**, *136*, 55. [ Links ]

17. Burrus, C.A.; *J. Chem. Phys* **1958**, *28*, 427. [ Links ]

18. Cacelli, I.; Moccia, R.; Montuoro, R.; *Phys. Rev. A* **2000**, *63*, 012512. [ Links ]

Received: March 24, 2005

Published on the web: January 13, 2006

FAPESP helped in meeting the publication costs of this article.

* e-mail: dlem@df.ufscar.br