Role of Si adatoms in the Si 111 -Au 5 2 quasi-one-dimensional system

Similar documents
LOW-TEMPERATURE Si (111) HOMOEPITAXY AND DOPING MEDIATED BY A MONOLAYER OF Pb

In order to determine the energy level alignment of the interface between cobalt and

Studying Metal to Insulator Transitions in Solids using Synchrotron Radiation-based Spectroscopies.

Spin-resolved photoelectron spectroscopy

Direct observation of a Ga adlayer on a GaN(0001) surface by LEED Patterson inversion. Xu, SH; Wu, H; Dai, XQ; Lau, WP; Zheng, LX; Xie, MH; Tong, SY

Supplementary Figures

Photoelectron Interference Pattern (PEIP): A Two-particle Bragg-reflection Demonstration

SUPPLEMENTARY INFORMATION

Lecture 5. X-ray Photoemission Spectroscopy (XPS)

Energy Spectroscopy. Ex.: Fe/MgO

Energy Spectroscopy. Excitation by means of a probe

Organic semiconductor heterointerfaces containing bathocuproine

Name: (a) What core levels are responsible for the three photoelectron peaks in Fig. 1?

High Resolution Photoemission Study of the Spin-Dependent Band Structure of Permalloy and Ni

Electronic structure of Sn/Si(111) 3 3: Indications of a low-temperature phase

SUPPLEMENTARY INFORMATION

Electron Spectroscopy

Energy-level alignment at interfaces between metals and the organic semiconductor 4,4 -N,N -dicarbazolyl-biphenyl

Properties of Individual Nanoparticles

Spectroscopies for Unoccupied States = Electrons

Surface Sensitivity & Surface Specificity

Published in Physical Review B 63, issue , 1-9, 2001 which should be used for any reference to this work

Application of single crystalline tungsten for fabrication of high resolution STM probes with controlled structure 1

A Momentum Space View of the Surface Chemical Bond - Supplementary Information

Surface Characte i r i zat on LEED Photoemission Phot Linear optics

Supplementary Figure 1 PtLuSb RHEED and sample structure before and after capping layer

Photon Interaction. Spectroscopy

Quantum Condensed Matter Physics Lecture 12

disordered, ordered and coherent with the substrate, and ordered but incoherent with the substrate.

Photoemission Spectroscopy

SUPPLEMENTARY FIGURES

Spin- and angle-resolved photoemission spectroscopy study of the Au(1 1 1) Shockley surface state

SUPPLEMENTARY INFORMATION

M2 TP. Low-Energy Electron Diffraction (LEED)

Author(s) Okuyama, H; Aruga, T; Nishijima, M. Citation PHYSICAL REVIEW LETTERS (2003), 91(

Spectroscopy of Nanostructures. Angle-resolved Photoemission (ARPES, UPS)

Electronic and atomic structures of a Sn induced 3 3x3 3 superstructure on the Ag/Ge(111) 3x 3 surface

STRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY

One-dimensional Atomic and Electronic Structures of Submonolayer Potassium Deposited on Stepped Surface of Ni(755)

Optimizing Graphene Morphology on SiC(0001)

Reflection high energy electron diffraction and scanning tunneling microscopy study of InP(001) surface reconstructions

Concepts in Surface Physics

SUPPLEMENTARY INFORMATION

Solid Surfaces, Interfaces and Thin Films

Unoccupied surface states on Si 111 ) )-Ag

The electronic properties of Si(001)-Bi(2 n)

Energy level alignment and two-dimensional structure of pentacene on Au 111 surfaces

Transverse momentum of ionized atoms and diatomic molecules acquired in collisions with fast highly-charged heavy ion. V. Horvat and R. L.

Supporting Information. for. Angew. Chem. Int. Ed. Z Wiley-VCH 2003

SUPPLEMENTARY INFORMATION

MODERN TECHNIQUES OF SURFACE SCIENCE

Imaging Methods: Scanning Force Microscopy (SFM / AFM)

Energy band of manipulated atomic structures on an insulator substrate

Supplementary Figures

X- ray Photoelectron Spectroscopy and its application in phase- switching device study

Observation of Bulk Defects by Scanning Tunneling Microscopy and Spectroscopy: Arsenic Antisite Defects in GaAs

Introduction to X-ray Photoelectron Spectroscopy (XPS) XPS which makes use of the photoelectric effect, was developed in the mid-1960

E kin. E vac. χφξ. Mott- Hubbard bands. E u CBM. E g E F. E o VBM. Electronic structure of SiC(0001) surfaces studied by two-photon photoemission

Surface Defects on Natural MoS 2

X-Ray Photoelectron Spectroscopy (XPS) Prof. Paul K. Chu

SUPPLEMENTARY INFORMATION

Supplementary Figure S1: Number of Fermi surfaces. Electronic dispersion around Γ a = 0 and Γ b = π/a. In (a) the number of Fermi surfaces is even,

X-ray Standing-Wave Investigation of (1X2)Rb/ Cu(110)

Scaling during shadowing growth of isolated nanocolumns

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 29 Sep 1999

Atomic and electronic structure of the Si 001 -Rb chemisorption system at 0.5 and 1.0 monolayer coverage

Spin pumping in Ferromagnet-Topological Insulator-Ferromagnet Heterostructures Supplementary Information.

SUPPLEMENTARY INFORMATION

Atomic and electronic structures of Si(111)-( 3x 3)R30 o -Au and (6x6)-Au surfaces

ELECTRON CURRENT IMAGE DIFFRACTION FROM CRYSTAL SURFACES AT LOW ENERGIES

Electronic Properties of Ultimate Nanowires. F. J. Himpsel, S. C. Erwin, I. Barke,

Experimental methods in physics. Local probe microscopies I

Scanning Tunneling Microscopy Studies of the Ge(111) Surface

Bremsstrahlung Isochromat Spectroscopy (BIS)

Electronic Supplementary Information

Electronic Structure of Surfaces

Soft X-ray Physics DELNOR-WIGGINS PASS STATE PARK

High resolution STM imaging with oriented single crystalline tips

Sponsored by. Contract No. N K-0073: Modification P00006 DARPA Order 5674 NR

Surface Studies by Scanning Tunneling Microscopy

a b c Supplementary Figure S1

CO Adsorption on Co(0001)-Supported Pt Overlayers

Origin of Metallic States at Heterointerface between Band Insulators LaAlO 3 and SrTiO 3

Minimal Update of Solid State Physics

Supplementary Information. for. Controlled Scalable Synthesis of Uniform, High-Quality Monolayer and Fewlayer

DFT EXERCISES. FELIPE CERVANTES SODI January 2006

Supporting Online Material for

Intrinsic vacancy induced nanoscale wire structure in heteroepitaxial Ga 2 Se 3 /Si(001)

Study of interface asymmetry in InAs GaSb heterojunctions

Semiconductor Physics and Devices Chapter 3.

Photon Energy Dependence of Contrast in Photoelectron Emission Microscopy of Si Devices

C. D. Lee and R. M. Feenstra Dept. Physics, Carnegie Mellon University, Pittsburgh, PA 15213

Scanning Tunneling Microscopy. how does STM work? the quantum mechanical picture example of images how can we understand what we see?

Surface and bulk electronic structure of thin-film wurtzite GaN

5.8 Auger Electron Spectroscopy (AES)

Electron momentum spectroscopy of metals

Dynamics of step fluctuations on a chemically heterogeneous surface of AlÕSi )Ã)

Extreme band bending at MBE-grown InAs(0 0 1) surfaces induced by in situ sulphur passivation

Hydrogen termination following Cu deposition on Si(001)

Auger Electron Spectrometry. EMSE-515 F. Ernst

Transcription:

PHYSICAL REVIEW B VOLUME 55, NUMBER 23 15 JUNE 1997-I Role of Si adatoms in the Si 111 -Au 5 2 quasi-one-dimensional system I. G. Hill and A. B. McLean Department of Physics, Queen s University, Kingston, Ontario, Canada K7L 3N6 Received 10 December 1996 The unoccupied electronic states of the Si 111 -Au(5 2) system have been studied, using inverse photoemission. Striking similarities exist between the unoccupied bands of the Si 111 -Au(5 2) system and the Si 111 -(7 7) and Si 111 -Au( 3 3)R30 systems. We argue that these similarities suggest that the Si 111 -Au(5 2) system contains Si adatoms, and that the adatom-induced surface states dominate both the inverse photoemission and photoemission spectra in the vicinity of the Fermi level. If this suggestion is correct, it means that the bands crossing the Fermi level may not be derived from the overlap of Au wave functions. The Au atoms may actually inhibit the formation of a Si adatom-induced band in the direction perpendicular to the Au chains. S0163-1829 97 04023-X INTRODUCTION Recent photoemission studies of quasi-one-dimensional systems, in their metallic phase, have demonstrated that there is vanishing emission intensity in the vicinity of the Fermi level. 1 5 This is in contrast to the behavior of two- and threedimensional metals, where well-defined Fermi edges are observed. 1 Although the studies mentioned above have been performed on bulk quasi-one-dimensional systems, it is also possible to fabricate quasi-one-dimensional, chainlike systems on surfaces. One example is the Si 111 -Au(5 2) system 6,7 Fig. 1. Although the details of this surface structure are still being worked out, it appears that the interchain spacing is significantly larger than the interatomic spacing along the chains. 6 So, at least from a structural point of view, the ordering of the Au atoms, on the Si surface, is quasi-onedimensional. In agreement with bulk quasi-one-dimensional systems, photoemission studies have found evidence for a strong anisotropy in the band dispersion. However, in contrast, they have observed a well-defined Fermi edge. 8 This difference in electronic properties begs the question, do bulk and surface quasi-one-dimensional metals belong to the same physical class? To investigate these questions further, we have performed inverse photoemission studies of the Si 111 -Au(5 2) system. We have paid particular attention to the region in the vicinity of the Fermi level, in an attempt to detect the Fermi edge, and the surface band identified in the photoemission studies 8 which crosses the Fermi level. fixed. The overall energy resolution of this system was estimated to be 0.6 ev full width at half maximum from the Fermi level emission onset of sputtered polycrystalline Ta, using the method described by Dose to extract the instrument pass function. 11 The overall momentum resolution of the system is both energy and angle dependent and has been esti- EXPERIMENT The experiments were performed in a -metal ultrahigh vacuum chamber with a base pressure of 7 10 11 torr. The inverse photoemission apparatus consists of a low-energy electron gun of the Erdman-Zipf design, 9 and a bandpass photon detector that utilizes a CuBe mesh electron multiplier Johnston Laboratories MM1 and a CaF 2 entrance window as described by Schäfer et al. 10 In contrast to Schäfer et al. we evaporated KCl onto the first dynode, to act as the photoemissive material. Angular dispersion is achieved by rotating the electron gun while the sample and detector remain FIG. 1. The Si 111 -Au(5 2) structure proposed by Bauer includes chains of Au atoms, which are common to many other models. The figure is included to define the crystallographic directions referenced in the text. In addition to the structure shown, Bauer also suggests that there are extra gold atoms on top of the chain rows, which increase the coverage from 0.4 to 0.45 ML. The spacing between these extra atoms is 4 2n bulk lattice constants along the chains, where n is a positive integer. The two parallelograms are the oblique 5 2 and 1 1 unit cells. 0163-1829/97/55 23 /15664 5 /$10.00 55 15 664 1997 The American Physical Society

55 15 665 ROLE OF Si ADATOMS IN THE Si 111 -Au(5 2)... mated to be between 0.1 and 0.2 Å 1 for the energies and angles in this study. A Princeton Research Instruments RVL 8-120 reverse view low-energy electron diffraction LEED system was used for sample alignment and surface characterization. The rotational symmetry of the diamond crystal structure about the 111 axis defines three equivalent 5 2 domains. However, it has been demonstrated that a single domain sample can be prepared using a vicinal surface. 12 The intentional vicinal miscut breaks the threefold degeneracy and favors the nucleation of a single domain, specifically that with Au chains parallel to the surface steps. The samples used in this experiment were cut from vicinal Si 111 wafers, miscut by 4 0.5 towards 112. All of the samples were p type, with resistivities of 5 cm. The overlayers were prepared in situ by direct current heating of the sample and deposition of Au from a Knudsen cell evaporator. All currents were applied in the 1 1 2 uphill direction to avoid step bunching on the surface. 13 A clean, well-ordered Si 111 surface was prepared by flashing the Si sample to 1100 C. The sample was heated to 650 C during the Au deposition. The quantity of Au deposited was calibrated using a quartz crystal rate meter. Coverages were measured by timed exposures to be between 0.40 and 0.45 ML, where 1 ML is defined to be the number of silicon atoms per unit area in the top layer of an ideally terminated Si 111 plane, 7.83 10 14 cm 2. The pressure during evaporation typically remained below 1 10 9 torr. LEED was used to determine the sample orientation and surface quality and to verify the presence of a single domain. Inverse photoemission spectra were recorded along the 1 10-11 0 azimuth, parallel to the Au chains, and along the 1 1 2-112 azimuth, perpendicular to the chains. These directions were chosen so that the present study would provide complementary information to that of the previous photoemission study. 8 Great care was exercised to ensure alignment of the 111 axis of the sample, which differs by 4 from the sample normal due to the vicinal miscut. It was found that the surface quality degraded as either a function of time, or by exposure to the electron beam. The quality of the LEED images was noticeably inferior after the collection of three inverse photoemission spectra. For this reason, the sample was annealed at 650 C for 2 min after the collection of every second inverse spectrum. The pressure during the collection of inverse photoemission data and annealing remained below 3 10 10 torr. The spectra presented are representative of the more than 15 samples studied. The shapes of the spectra and the peak positions were found to be reproducible from sample to sample. RESULTS The sample flashing process described in the previous section reliably produced a sharp 7 7 pattern, as shown in Fig. 2 a. The splitting of the 1 1 spots, visible at LEED energies around 40 ev, was taken as evidence of an evenly spaced step/terrace structure on the vicinal surface. After Au deposition, the LEED patterns consisted of a well-defined 5 1 pattern, with diffuse half-order streaks Fig. 2 b. The weakness of the half-order intensities has been attributed to a lack of registry between well-ordered neighboring chains. 6,14 FIG. 2. a The top panel shows a typical 7 7 LEED pattern prior to Au deposition. b The bottom panel is an example of a single-domain Si 111 -Au(5 2) LEED pattern. The sample is oriented such that the 112 direction is pointing upward. In both cases, the energy of the electrons is 80 ev. The Au chains are horizontal. Examples of the spectra collected perpendicular to the chains are included in Fig. 3. There are three distinct features, two of which overlap near the zone center. The feature seen at 2.7 ev above the Fermi level at normal incidence is believed to consist of two states: one that exhibits no detectable dispersion, and one that disperses away from the Fermi level reaching 3.5 ev at 15 off normal. The third feature occurs at 4.5 ev above the Fermi level at normal incidence, and disperses sharply upward, reaching 5.5 ev at 10. In addition to these three features, there is an additional, less obvious, feature that appears less than 1 ev above the Fermi level, although the exact energy position is difficult to determine as it only appears as a shoulder on the much stronger 2.7 ev peak. This shoulder is visible at virtually every angle along both azimuths studied. At no point is there any evidence of a state crossing the Fermi level. Possible origins of the observed features will be discussed in the following section. Spectra collected along the chains are shown in Fig. 4. The states seen in this direction are very similar to those seen perpendicular to the chains. Two states appear to overlap near the zone center, one with negligible dispersion and the other dispersing upwards with increasing angles. A third feature is visible at 4.5 ev at normal incidence, which disperses

15 666 I. G. HILL AND A. B. McLEAN 55 FIG. 3. Inverse photoemission spectra collected with the projection of the electron wave vector onto the surface lying perpendicular to the Au chains. The positive numbers on the right are angles, in degrees, toward 112. Negative angles are toward 1 1 2. The numbers in parentheses are the components of the electron wave vector parallel to the surface in Å 1 calculated at E F 2.7 ev. The zone boundary is reached at a wave vector of 0.142 Å 1. The top spectrum shows the position of the Fermi level on sputtered polycrystalline tantalum in electrical contact with the sample. Tick marks are provided as a guide to the eye. FIG. 4. Inverse photoemission spectra collected with the projection of the electron wave vector onto the surface lying parallel to the Au chains. The positive numbers on the right are angles, in degrees, toward 11 0. Negative angles are toward 1 10. The numbers in parentheses are the components of the electron wave vector parallel to the surface in Å 1 calculated at E F 2.7 ev. The zone boundary is reached at a wave vector of 0.409 Å 1. The top spectrum shows the position of the Fermi level on sputtered polycrystalline tantalum in electrical contact with the sample. Tick marks are provided as a guide to the eye. sharply upward with increasing angles. The shoulder around 1 ev is again seen, and there is no evidence of states crossing the Fermi level. DISCUSSION It is informative to first compare the normal incidence spectra of the Si 111 -Au(5 2) and clean Si 111 -(7 7) surfaces as shown in Fig. 5. The state labeled U 1 appears at the same energy and, when the spectra are normalized to the incident electron flux, with the same intensity on both the clean and Au reconstructed surfaces. On the clean 7 7 surface, this state has been shown to be spatially localized on the Si adatom 15 and it is thought to be a surface state that is derived from the interaction of Si adatoms, which sit in the T 4 sites, with the dangling bond orbitals of the substrate. 16 19 Furthermore, the U 1 state on the 5 2 surface is strongly attenuated after approximately one day in ultrahigh vacuum, where the predominant residual gas is hydrogen. The other features remain unaffected. The sensitivity of U 1 to contamination is similar to that of the adatom state on the 7 7 surface. 17 The binding energy of the U 1 state and the time dependence of the state intensity implies the presence of Si adatoms in the 5 2 structure. If this suggestion is correct then the fact that the U 1 state has the same intensity in both systems suggests that the Si adatom density is comparable. On average, each 5 2 unit cell would require 2.4 Si adatoms to have the same number as the 7 7 structure 12 in the 7 7 cell. Although not all of the structures that have been proposed so far include Si adatoms, there is a model that does 7 and the Si adatom density is comparable to the density we extract from our spectra. Marks and Plass 7 propose a structure that requires one Si adatom per unit cell and suggests that additional adatoms may be present. The remaining features of Fig. 5 bear a strong resemblance to the features seen on Si 111 -Au( 3 3)R30. The states labeled A and B have been previously identified as bulk transitions. 20 The bulk transition B looks sharper on the Si 111 -Au( 3 3)R30 surface than on the clean 7 7 surface, as it does here. This sharpening has been previously attributed to a reduction in umklapp scattering by surface reciprocal lattice vectors compared to the clean surface. 20 The state labeled U 2 appears very similar to a state that

55 15 667 ROLE OF Si ADATOMS IN THE Si 111 -Au(5 2)... FIG. 5. Comparison of clean Si 111 -(7 7) and Si 111 - Au(5 2) spectra at normal incidence. Feature U 1 is thought to arise from Si adatoms. Note that the normalized intensity of U 1 is the same for both the Si 111 -(7 7) and the Si 111 -Au(5 2) systems, indicating a comparable surface adatom density. arises from Si-Au bonding on Si 111 -Au( 3 3)R30. 20 If the U 1 band on the 5 2 surface does arise from Si adatoms, then this has important implications for the way we think about the electronic structure of this system. On the Si(7 7) surface both the unoccupied U 1 band and the occupied S 1 band have been shown to be adatom derived. Both ab initio density functional theory calculations 16 and studies of the surface using scanning tunneling microscopy spectroscopy 15 have made this association. Consequently, the occupied surface band that has been studied with photoemission 8 may not arise from the overlap of the Au wave functions. The Au atoms may actually inhibit the overlap of Si adatoms in the direction perpendicular to the Au chains, and this may have the effect of reducing the bandwidth of the adatom-induced bands. So the quasi-onedimensional behavior that has been observed may only indirectly be a result of the presence of the Au chains. Photoemission studies of the 7 7 system 21 have identified a relatively flat surface band that is believed to extend above the Fermi level. This metallic state is responsible for pinning the Fermi level 0.7 ev above the bulk valenceband maximum. 19 The dispersion of this band is not unlike the dispersion observed so far on 5 2. Furthermore, inverse photoemission studies of the 7 7 system have failed to see this metallic state emerge as a peak above the Fermi level, although metallic tails are observed. 17 This is completely analogous with the behavior we observe on the 5 2 system. Therefore, the picture we are developing here, that the near Fermi-level states are adatom derived, seems consistent with existing studies of the 7 7 system. 16,17,15 However, we still have one major problem to resolve. Why does the surface band on the 5 2 system possess a well-defined Fermi edge? Does this mean that the system is not truly one dimensional? The lack of dispersion of the occupied surface band perpendicular to the chains clearly indicates that the system is one dimensional. We believe the answer to this paradox will be resolved with high-resolution photoemission studies of the Fermi-level region. We note that an early study of the bulk quasi-one-dimensional material K 0.3 MoO 3 found evidence for a well-defined Fermi edge. 22 However, more recent studies, performed at considerably higher-energy resolution ( 20 mev compared with 150 mev conclusively demonstrated that the system had vanishing emission intensity at the Fermi level. 1 The photoemission study of the 5 2 system was performed with an energy resolution of 125 mev. 8 This is insufficient to unambiguously determine whether there is a well-defined Fermi edge or not. High-resolution photoemission studies of this interesting system are in order. CONCLUSION The unoccupied states of the Si 111 -Au(5 2) system have been studied with inverse photoemission. It was found that the inverse photoemission spectra possessed a feature with a binding energy, intensity, and sensitivity to contamination that are similar to the adatom-induced feature on the Si 111 -(7 7) surface. We have explored the possibility that this feature (U 1 ), and also the uppermost occupied surface band (S 1 ), arise from the overlap of the Si adatoms with dangling-bond orbitals of the substrate. In this scenario, the one-dimensional behavior that has been observed in a previous photoemission study, 8 would arise from the anisotropic distribution of adatoms on the surface and not from the chainlike distribution of Au atoms. In fact, it is possible that the Au chains may inhibit the formation of a Si adatom band perpendicular to the Au chains. We have also argued, using the historical precedent of the K 0.3 MoO 3 system, that it is too early to conclusively state that there is a well-defined Fermi edge for the Au(5 2) system. The photoemission studies have to be performed at significantly higher-energy resolution. Since the adatoms interact strongly with the substrate, 18 it will be interesting to see whether the valence band emission possesses one- or two-dimensional character. ACKNOWLEDGMENTS The authors gratefully acknowledge the technical assistance of J. Thompson in machining the inverse photoemission and sample manipulation apparatus. We would also like to thank H. Tran of the National Research Council of Canada for offering her time to dice our samples. This work was supported by the Natural Sciences and Engineering Research Council of Canada.

15 668 I. G. HILL AND A. B. McLEAN 55 1 B. Dardel et al., Phys. Rev. Lett. 67, 3144 1991. 2 M. Nakamura et al., Phys. Rev. B 49, 16 191 1994. 3 C. Coluzza et al., Phys. Rev. B 47, 6625 1993. 4 B. Dardel et al., Europhys. Lett. 24, 687 1993. 5 A. Sekiyama et al., Phys. Rev. B 51, 13 899 1995. 6 E. Bauer, Surf. Sci. Lett. 250, L379 1991. 7 L. D. Marks and R. Plass, Phys. Rev. Lett. 75, 2172 1995. 8 I. R. Collins et al., Surf. Sci. 325, 45 1995. 9 P. W. Erdman and E. C. Zipf, Rev. Sci. Instrum. 53, 225 1982. 10 I. Schäfer et al., Rev. Sci. Instrum. 58, 710 1987. 11 V. Dose, Prog. Surf. Sci. 13, 225 1983. 12 J. D. O Mahony et al., Surf. Sci. Lett. 277, L57 1992. 13 H. M. H. Tokumoto, K. Miki, and K. Kajimura, J. Vac. Sci. Technol. B 9, 699 1991. 14 H. Lipson and K. E. Singer, J. Phys. C 7, 12 1974. 15 R. J. Hamers, R. M. Tromp, and J. E. Demuth, Surf. Sci. 181, 346 1987. 16 K. D. Brommer, M. Galván, J. A. Dal Pino, and J. D. Joannopoulos, Surf. Sci. 314, 57 1994. 17 J. M. Nicholls and B. Reihl, Phys. Rev. B 36, 8071 1987. 18 J. E. Northrup, Phys. Rev. Lett. 57, 154 1986. 19 W. Mönch, Semiconductor Surfaces and Interfaces, 2nd ed. Springer, Heidelberg, 1995. 20 J. M. Nicholls, F. Salvan, and B. Reihl, Surf. Sci. 178, 10 1986. 21 P. Mȧrtensson et al., Phys. Rev. B 36, 3015 1987. 22 J. Y. Veuillen, R. C. Cinti, and E. A. K. Nemeh, Europhys. Lett. 3, 355 1987.