SCF discrete-variational Xα calculations of the NH4Cl Auger spectrum

SCF discrete-variational Xα calculations of the NH4Cl Auger spectrum

CHEMICAL PHYSICS I. IntrodIIctioIl \ugc~ spc’~-Iros~c~p~ JII\csII~.IIJ~I oI‘Il1~ 15 wideI>’ &x-Ironic used structure 3s J Illefhod for of ...

456KB Sizes 0 Downloads 10 Views

CHEMICAL

PHYSICS

I. IntrodIIctioIl

\ugc~ spc’~-Iros~c~p~ JII\csII~.IIJ~I

oI‘Il1~

15

wideI>’

&x-Ironic

used

structure

3s J Illefhod

for

of nlolccules,

.\llhdS Jld c~lIII~WJlIdS (111 SlldXCs 11 ]_ %\VCVeT, the ~w~prmiiwt 01 Auger bpcctn veryuitenneedsquite III\ uh’d

~YJh1~~11~111~~

111~Ihc1ds.

Some

cskuk!tions

c.mIed WI and they WCIC able to describe Iirx4~qm .ind lint structures of Auger spectra quanIII.III\C~! . In p.IIIImlJr. .Ib Initio calcuiarions of Augc~ cnCIpies JS differences of the initial- and fin,& ‘~I.IJL! Ik11.drnerpm glvc good agrecnient with e.\periIllen 13 61. I Io\seVer. Ihe clnpioywx~t of ab initio 1J1~1hld~ IS IcsI~Icrcd Iu ~~11 niolecules for computaIICIIIJI 15.~boJJ5. CbJJlg Iiie srniicInpIric~1 approxh of JL’IIIIISOI~17.51 me nny c.~kul~Ic Auger transitiori ~SI~;ICS I;lr qu~fc’ Iic0y mulcculrs. Unfortunately. I~IIS.IpprL).xh m]IIIrcs knowledge of the energy of the hQ hole mIc’~.Ic’tJorI and c.\perimental ionization p~re1ir1;115 (II’) \\liich ni~y not be known in every case. .\ c~lIllpJOIllISC be~\ve~n db initio Inethods and

II.I\LT l1w11

~~lllI-~Jll~lirlCJ~

.J~qmlJcheS

is represented

by

1983

nectcd wirb calculations of energies and intensities of the Auger transitions within the Xa model are yet to be answered. In this work the self-consistent discrete-variational Xa Inethod(DVh1) 1151 is applied to calculations of the KVV hueshape for solid NH,Cl_ A new formalim 1o determine the singlet-triplet splitting of the Auger transition energies is suggested.

2. Calculations

of the Auger energies and line&apes

The kinetic energy of an Auger electron is equal to the difference of total energies for an initial state with a vacancy on the inner shell and a final state with fwo vac3ncies on valence NOS.

(1)

the

~IicIllo~is ~vtliclJ e~nploy the local eschange potentials, c g _\a. l~~~cntly SOIW applications of such methods N err’ Icportrd. for example. LCXO Xa to O2 [9] and It> SdJd ~YNO; 1lo]. SW Sa IO CF4 11 I], Co [I?] .ml .~t~ms 1131. ,1nt1 lo La-,‘33 in a cluster model 11-i I. I IuWvcx. sonw inctllodologicsl questions con-

70

75 hlarch

LETTERS

where $‘i are occupation numbers of the ith spinorbitals with spin oi participating in the transition and X is the final electronic configuration. For systems whose ground state is a closed shell, and NH,Cl in particular, the electronic configuration of the final state with two vacancies on valence MOs 0 0092614/83/0000-0000/S

03.00 0 1983 North-Holland

may be either singlet or triplet.

Because in +&e aingle-

determinant spin-polarized approximation of the wavefunction, spin and orbital angular momentum

G?t33 _

=O,n~~=O);

=E(II~

=0)-E,,,,

,

then average energies of singlet (S) and triplet (T) multiplets according to this sum rule are equal to ET@“2 = 0 ,n$ =0)= Es@

Ea2_

,

-4%pY3

(2)

(3)

-

From eqs. (l)-(3) for Auger energies of S and T states we obtain the following expressions:

(4) K&3 = 2&$a

-K&3

-

,n3 =1)-E(r~l=1,tr2=0,n3=0)

+-.&+e.&

= -e&

(7)

with an accuracy up to the third terms of (6). Here e&are orbital energies detemrined for TS. Estimations of the Auger energies using (4) and (5) require in the common case two spin-polarized calculations for every transition with occupation numbers determined by the corresponding configuration of configuraholes. The calculation on a clfla~~P-a~~~a tion gives immediately the triplet component of the Auger transition K 53 = -ey

= 0 , n$ = 0) = ET - 2(ET - Eaza3)

=2Gzp3

algebra E(rtl=0,n2=1

multiplet splitting cannot be treated uambiguously, we use in calculations of multiplet energies the sum rule [ 161. In this approach the energies of singlet and triplet states are averages of all states (with unknown weights in a common case) which arise from the given electronic configuration. if we define E o2 + = E(n;’

25 March 1983

CHEMICAL PHYSICS LETTERS

Volume 96, number 1

+e7Q +&J .

@I

The calculation-k a c1flav~/2%~/2fl configuration yields the energy Kt& = -rzy + &q + $ which is required by (5) to determine the singlet component of a transition_ It is pertinent to emphasize the model character of the TS computations_ The high multiplicity in the Auger energy calculations is due to the mathematical formalism [see (6) and (7)] and not to the physics of the Auger process. For closed-shell ground-state molecules the intensities of the Auger transitions are connected with Coulomb

(9

The spin index in K123 is omitted because the difference between hdoZa3 and K$a o3 does not exceed tenths of one eV according to our calculations. For valence orbitals this index is omitted due to equality ofEazp3 and E8ap3 in the case of closed-shell systems. Now we expand the total statistical Xo energy in a Taylor series in occupation numbers of MOs involved in the Auger process,

(6) and following the Slater transition state(TS) concept [ 171, select a state in the middle of the initial (11~ = 0, I’2 = 1, "3 = 1) and final (n, = 1, “2 = 0, it3 = 0) states. For the TS should be assigned occupation numbers nl = l/2, n2 = 112 and ns = 112. Using the expansion (6) with respect to TS one may deduce after some

where Q and b denote valence MOs involved in the Auger process and Q(E) is the wavefunction of an Auger electron. In the one-center approsimation which has been commonly used [2,5,6] Q(E) is a spherical wave centered at the atom with the primary vacancy (at N in our case)_ @I,,,(e) = RI(~) YI,,, >

(10)

where R](E) is the radial function of a continuum wavefunction and Ylm is the spherical harmonic_ The one-electron MOs are determined in DVhl by the usual LCAO expansion: pi = c

i

c&a-

J

(11) 72

Voluttxe 96. number 1

CHEhflCAL

where X~ are atomic basis functions. In the one-center approximation

where pi are atomic functions centered at an atom with the primary vacancy. The matrix elements in (117) Jlnve been calculated [ 181 using HFS numerical functions. Here we use numerical HF functions which are practically tile same as HFS ones. Notice that the Auger energies given in the tables are calculated within the doubfc-zeta basis of Clenlenti [20]. Auger energies obtairwd within tile HF basis slightly differ from tile latter ones mainly due to the distinction of calculdtcd energies of the N Is A0 in both bases. The N 1s 11’is 3 11.9 or 406.3 CV against the experimental v&c 403 eV [ 191 of a free atom. Only main Auger lines are calculated witllout taking into account llie processes of shake-up and shake-off types The half-~ idths of a gaussian decomposition of the t?rpcrimental spectrum of NH_+Cl are used wflen plotting the tfleoreticai spectrum.

3. Results and discussion The geometry

of the NJ-J_+Clmolecule

25 hfareh 1983

PHYSICS LFXI-ERS

is taken

partly from the X-ray data on solid NH,CJ 1211, i.e. the point group is C&, R(N-H) = 2.1 au, and from the ab initio calculations of Clementi [22], i-e_ R(H-Cl) = 3-4 au- The hydrogen atoms are assumed to form an undistorted tetrahedron. The cafculations are carried out as within HF using Clementi’s doublezeta f201 functions on N and CLand fs the H. In table 1 we give the IPs, Mulliken charges on atoms and experimental 1Ps of the constituent atoms N, CJ and H, table 2 presents the energies of Auger transitions calculated within the dot&e-zeta basis according to (5) and f7), and the integraf intensities determined by (I f) using the LCAO coefficients of the ground-state wavefunctions within the HF basis. Because the number of fine-structure lines for the main Auger transitions is quite large, seven peaks are selected for comparison with features of the expcrimental spectrum. The theoretical intensities are caJculated to be integral intensities of Auger Iines determining the peaks. In table 3 are given the calculated Auger energies and relative intensities of peaks 1-7 and the experimental v&es obtained by a Gauss decomposition of the experimental spectrum. The latter is recorded for polycrystalline N&Cl using the Mg I& We_ A more detailed description has been published elsewhere [23]. The spectrum after substraction of the back-

l‘&lc I The

wIct11~~d

1%. MuIlAen

ctuqe+

on

Itzoms of NH~CI. and rhe experimental

IPs of N. H and Cl atoms (in eV)

Experimental IFS a)

bfcJ double-zeta - ._.--------7.8

8.2

7-9

17.4 18.9 20.3 19.1

17.2 18.9 20.0

3VI

29.2

2SN

207.0 208.5 264.6 411.9 2832.9

Qll

Qh

---.. 3)

72

- -.. -

horn

ret 119).

I ----

QN QCI

-___-.__

I_~-

7.9

0.37 0.30 -0.80 -0.62

13.7 13.2 13.6

3PC1 2px 1SH

199.8 201-4 258.9 406.3 2806.5 0.28 0.70 -0.40 -0.60

2%2,3, 2sfl 1Sh’

kl

Cl

25.3 26.0 206 204 274 403 2827

CHEMICAL PHYSICS UZl-TERS

Volume 46. number 1

25 March 1983

Table 2 The calculated Auger energies (in ev) and integral intensities of the N KVV transitions in NH,Cl Confiiuration

Spin

Kinetic energy

integral intensity

Number of peak

0

336.8

0.590

1

0

347-O 349.6

0.026

1

2

1

345.6 350.4

0.730 I 0.112

Sa;”7ai'

0 1

351.2 352.2

0.528 0.074 I 0.170

6a;* 6a;’ Ze;’

0 0 1

359.7 358.8 360.4

0

1

358.3 361.7

0

361.9

1

362.3

0

362.7 362.9

0

3

1

0 0

1 6a;’ &I;’

0 1

73;’

sa;”

8ai’ 2e;’

8ai2

0

362.8 362.7 362-V

I.020 0.038

1

0.379

0.010

372.0 374.2

0.536 0.020

373.7

0.666

373.6

0.034 O-1 24 j

381.6

0.274

eV to match the first peaks of the

spectra. This shift is due to different phases of NH4CI in the calcnlations and the experiment. After the shift

a close resemblance of both curves is evident. Note that the more traditional approach in calculations for solids is a cluster model within which one should calculate the cl?ster NH4@-. However, as

I

5

0.018 J 0.584

368.1 370.3

ground is presented in fig. 1 together with the theoretical one. The experimental Auger energies are

displaced -9.5

4

6

7

has been noticed f24], a cluster simulating a solid should preferab3y consewe the stoichiometry of the solid. From this point of view the molecule NH&C1is the minimal stoichiometric cluster for solid NH,Cl. Most calculated Auger lines are located near experimental peak positions within limits of experimental errors. The exceptions are peaks 4 and 7 orig-&ating from the 2e- 2 and Sai_ 2 configurations.

agreement for 2e-’

The dis-

may be explained by the one-de73

Volume

96. number

CHEhIICAL

1

T.Me 3 The c\perimcnwl and calculated enersies(in intensirxs for the N KVV i\uger transitions ~-_ . -~- _-__-_ Number

EIhcor

EexP

&or a)

(C’V)

(eV)

(5)

lfheor IS G)

Iexp (%)

1 ,

336.8 34x5 351.2 35s.3 362.7 371.8 3Sl.h

336.S 317 352 362 363 373 379.5

7.7 11 10 ‘2 16 I9 3.7

5.6 13 5.6 41 14 18 2.8

9.2 11 10 20 23 22 4.1

3) fi;s arc hcor and IIS lllr: I \tb LCt$@ coetkicnts

inre_er~lintcnsitirs wlculated using for the ground- and initial-state ~.~\ciuncrionc respccti\ely. The expcrinlental Au_ecrenerSIC<.IIL’zhlftcd h? -9.5 eV.

Icrnunxlt approximation used because we do not resolve rhe multlplet structure of this configuration and dccounf for the configurarion interactions. and for S.I~ ’ ir may be connected with a lattice potential inllucnrc. The theoretical inrensities of peaks l-7 are 111a good agreement with experimental ones if one uses t!lc l-C.40 coefficients of the ground-state wave-

function. Thr intensilies were calculated using the Xl0 LCAO coefficients of the initial state wavefunctiorl with a vdcmcy in the K shell of N as well. In the

H--

LE-ITERS

25 March 1983

latter case the contribution

eV) and peak in NbCl

of pc.A

; J 5 6 7

PHYSICS

of p AOs does increase to the influence of the core vacancy field, leading to a considerable relative Increase of amplitude of peak 4. The same trend is observed when calculating the intensities with the use of LCAO coefficients for TS of the 7ai’2e? ’ configuration. Recently [23 ] we have qualitatively assigned the fine structure for the N KVV Auger spectra of NHdX, X = Br and Cl, assuming the approximate equivalency of the NH; and CH4 electronic structures, accounting for their isoelectronic character. The high-energy peaks 7 and 8 were considered to be satellite ones for transitions from the highest valence MOs composed mainly of N AOs. According to results of the DVM calculations the appearance of peaks 7 and 8 (see fig_ 1) is due to a substantial delocaliiation of the N-H-Cl bond owing to hybridization of N 2s, 2p,, the bridge H Is and Cl 3s. 3p, AOs. It is interesting to note that the shape of the other peaks l-6 is similar to that of the carbon Auger lines in CH,. It allows us to connect the intensities of peaks 7 and 8 with the covalency strength of the Cl-NH, bond. due mainly

Acknowledgement The authors wish to express their sincere gratitude to Drs. Ju.M. Shul’ga and A.F. Shestackov for discussion and critical remarks, and to the referee for valuable comments_

References

h--H-_CI

[ 1 ] R. Weissmann and K. hl’uller, Surface Sci. Rep. 105

.-__

307

_. *

327

347

367

387 Et,”

ieVl-

1 lg. 1. 1 he c.dsul~Icd and cxpcrmxnt.d N KVV Auger specII~ of jY14Cl. The positions of theoretical peaks 1-7 are tn.~rkxl x\lth \cr~x.d bxs. Their amplitudes are proportional to Ihc pc.k mrcnsirws. The pe_k to the left is due to admixturc 01 thr 0 KCUcwitarlon ofC1 2p levels.

73

(1981) 251. [ 2 ] H. Siegbahn. L. Asplund and P. Krlfve. Chem. Phi-s. Letters 35 (1975) 330. [3 J 1~. Agren, S. Svcnsson and U-1. Wahlgren. Chem. Phys. Letters 35 (1975) 336. [?I hl.T- Okland, K. Faegri Jr. and R. Mannc, Chem. Phys. Letters40 (1976) 185. 151 K. Faegri Jr. and H.P. Kelly. Phys. Rev_ A19 (1979) 1649. [6] H. Agren. J. Chem. Phys. 75 (1981) 1267. [7] D.R. Jennison, Chem. Phys. Letters 69 (1980) 435. [S] D-R. Jennison, J. Vat. Sci. Technol. 17 (1980) 172. (91 B-1. Dunlap. P.A. hlills and D.E. Ramaker, J_ Chem. Phys. 75 (1981) 300. [IO] F-L. Hutson. D.E. Ramalter, B-1. Dunlap, J-D. Ganjei and J-S. hfurday. J. Chcm. Phys. 76 (1982) 2181.

Volume 96, number 1

CHEMICAL PHYSICS LETTERS

111 J hf. Barber, I.D. Qark and A. ~~~h~ffe, Chem. Phys. Letters 48 (1977) 593. 1121 P.C. Desmukh and R.G. Hayes, Chem. Phys, Letters 88 (1982) 384. 1131 K-D. Sen, 3. Chem. Phys. 73 (1980) 4704. fl4] LA. Topal, privtrte communicatioa [ ISI G-L. Gutsev and A.A. Levi& Chem. Phys. 51 (1980) 459. [i6] T. Ziegler, A. Rauk and E.J. Baerends, Theoret. Chim. Acta 43 (1977) 261. iI71 3-C Slates, The self-consistent field for molecules and solids @lit, Moscow, 1978) ch. 2_ fl8l EJ. McGuire_ Phyr Rev_ 185 (1969) l_

25 March 1983

[ 19 I K. Siegbahn, C. NordIin#& k Fahhnan, R Nordberg, K. Hamrin, J. Hedman,G. Johanssdn,T. Bergmark, S.-E. Karlsson, I. Lindgren and B. Lindberg, ESCA, Atomic, molecular and solid state structures studied by means of electron spectroscopy (MI, Moscow, 1971) table43. I201 E. Clementi. R Matcha and A. Veillard, J. Chem. Phys. 47 (1967) 1865.

1211 Structure Reports 23,295, (221 E. Clementi, J. Chem. Phys. 46 (1967) 3851. (23 1 G-M. hlichailov and Ju.G. Borod’ko, Surface 6 (1982) 85 [in Russian]_ 1241 C Umrigu and D-E. Ellis, Phys. Rev. B21 (1980) 652.