Hydrogen termination following Cu deposition on Si(001)

Similar documents
A report (dated September 20, 2011) on. scientific research carried out under Grant: FA

Report Documentation Page

REGENERATION OF SPENT ADSORBENTS USING ADVANCED OXIDATION (PREPRINT)

System Reliability Simulation and Optimization by Component Reliability Allocation

Use of Wijsman's Theorem for the Ratio of Maximal Invariant Densities in Signal Detection Applications

P. Kestener and A. Arneodo. Laboratoire de Physique Ecole Normale Supérieure de Lyon 46, allée d Italie Lyon cedex 07, FRANCE

Periodic Magnetoresistance Oscillations in Side-Gated Quantum Dots

CRS Report for Congress

Crowd Behavior Modeling in COMBAT XXI

Modulation Instability of Spatially-Incoherent Light Beams and Pattern Formation in Incoherent Wave Systems

Attribution Concepts for Sub-meter Resolution Ground Physics Models

Quantitation and Ratio Determination of Uranium Isotopes in Water and Soil Using Inductively Coupled Plasma Mass Spectrometry (ICP-MS)

REPORT DOCUMENTATION PAGE. Theoretical Study on Nano-Catalyst Burn Rate. Yoshiyuki Kawazoe (Tohoku Univ) N/A AOARD UNIT APO AP

Analysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions. Dr. Bijan Khatib-Shahidi & Rob E.

Metrology Experiment for Engineering Students: Platinum Resistance Temperature Detector

Diagonal Representation of Certain Matrices

Report Documentation Page

VLBA IMAGING OF SOURCES AT 24 AND 43 GHZ

High Resolution Surface Characterization from Marine Radar Measurements

Babylonian resistor networks

REPORT DOCUMENTATION PAGE

PIPS 3.0. Pamela G. Posey NRL Code 7322 Stennis Space Center, MS Phone: Fax:

Closed-form and Numerical Reverberation and Propagation: Inclusion of Convergence Effects

Real-Time Environmental Information Network and Analysis System (REINAS)

DIRECTIONAL WAVE SPECTRA USING NORMAL SPREADING FUNCTION

Impact of Typhoons on the Western Pacific Ocean DRI: Numerical Modeling of Ocean Mixed Layer Turbulence and Entrainment at High Winds

Electrodeposition of nanomaterials. W. Schwarzacher H. H. Wills Physics Laboratory, University of Bristol

Optimizing Robotic Team Performance with Probabilistic Model Checking

INFRARED SPECTROSCOPY OF HYDROGEN CYANIDE IN SOLID PARAHYDROGEN (BRIEFING CHARTS)

FRACTAL CONCEPTS AND THE ANALYSIS OF ATMOSPHERIC PROCESSES

Exact Solution of a Constrained. Optimization Problem in Thermoelectric Cooling

Marginal Sea - Open Ocean Exchange

SW06 Shallow Water Acoustics Experiment Data Analysis

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER

Unusual Optical Properties of Aligned Carbon Nanotube Mats in Infrared Energy Region

Computer Simulation of Sand Ripple Growth and Migration.

Determining the Stratification of Exchange Flows in Sea Straits

Dynamics of Droplet-Droplet and Droplet-Film Collision. C. K. Law Princeton University

Development and Application of Acoustic Metamaterials with Locally Resonant Microstructures

Convection and Shear Flow in TC Development and Intensification

Volume 6 Water Surface Profiles

Ocean Acoustics Turbulence Study

INFRARED SPECTRAL MEASUREMENTS OF SHUTTLE ENGINE FIRINGS

DISTRIBUTION A: Distribution approved for public release.

Scattering of Internal Gravity Waves at Finite Topography

Report Documentation Page

The Bernoulli equation in a moving reference frame

MTIE AND TDEV ANALYSIS OF UNEVENLY SPACED TIME SERIES DATA AND ITS APPLICATION TO TELECOMMUNICATIONS SYNCHRONIZATION MEASUREMENTS

Thermo-Kinetic Model of Burning for Polymeric Materials

Heat Conduction by Molecular Dynamics Technique

Estimation of Vertical Distributions of Water Vapor and Aerosols from Spaceborne Observations of Scattered Sunlight

Hurricane Wave Topography and Directional Wave Spectra in Near Real-Time

SENSORS FOR MEASURING THE VOLUME SCATTERING FUNCTION OF OCEANIC WATERS

Sensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts

On Applying Point-Interval Logic to Criminal Forensics

The Navy Precision Optical Interferometer for SSA applications: an update

Predicting Tropical Cyclone Formation and Structure Change

Predictive Model for Archaeological Resources. Marine Corps Base Quantico, Virginia John Haynes Jesse Bellavance

Broadband matched-field source localization in the East China Sea*

Coastal Engineering Technical Note

High-Fidelity Computational Simulation of Nonlinear Fluid- Structure Interaction Problems

Dynamics of Boundary Currents and Marginal Seas

Mine Burial Studies with a Large Oscillating Water-Sediment Tunnel (LOWST)

The Mechanics of Bubble Growth and Rise in Sediments

Flocculation, Optics and Turbulence in the Community Sediment Transport Model System: Application of OASIS Results

NAVGEM Platform Support

Grant Number: N IP To compare obtained theoretical results with NPAL experimental data.

USMC Enlisted Endstrength Model

An Observational and Modeling Study of Air-Sea Fluxes at Very High Wind Speeds

Wavelet Spectral Finite Elements for Wave Propagation in Composite Plates

EFFECTS OF LOCAL METEOROLOGICAL VARIABILITY ON SURFACE AND SUBSURFACE SEISMIC-ACOUSTIC SIGNALS

Understanding Near-Surface and In-cloud Turbulent Fluxes in the Coastal Stratocumulus-topped Boundary Layers

Z-scan Measurement of Upconversion in Er:YAG

Analysis of Infrared Measurements of Microbreaking and Whitecaps

Topographic Effects on Stratified Flows

Internal Tide Generation in the Indonesian Seas

INTERACTION AND REMOTE SENSING OF SURFACE WAVES AND TURBULENCE

Nano and Biological Technology Panel: Quantum Information Science

Characterization of Caribbean Meso-Scale Eddies

Effects of Constrictions on Blast-Hazard Areas

Erik L. Swanberg 1 and Steven G. Hoffert 2. Veridian Systems 1, Autometric 2. Sponsored by Defense Threat Reduction Agency

Report Documentation Page

IMPROVED SYNTHESIS AND REACTION CHEMISTRY OF FN 3

G.H. Share, W.N. Johnson, J.D. Kurfess, R.J. Murphy. A. Connors. B.L. Dingus, B.E. Schaefer. W. Collmar, V. Schonfelder

Curating and Maintenance of the Sediment Library and Dredge Collection

Improvement of Mesoscale Numerical Weather Prediction For Coastal Regions of Complex Terrain FY2003

Super-Parameterization of Boundary Layer Roll Vortices in Tropical Cyclone Models

Vortex Rossby Waves and Hurricane Evolution in the Presence of Convection and Potential Vorticity and Hurricane Motion

A 1/10th Degree Global Ocean Simulation Using the Parallel Ocean Program

Modeling the Formation and Offshore Transport of Dense Water from High-Latitude Coastal Polynyas

Discovery of Planetary Systems With SIM

STUDY OF DETECTION LIMITS AND QUANTITATION ACCURACY USING 300 MHZ NMR

Shallow Water Fluctuations and Communications

Study of Electromagnetic Scattering From Material Object Doped Randomely WithThin Metallic Wires Using Finite Element Method

Coupled Ocean-Atmosphere Modeling of the Coastal Zone

Mass Transport by Second Mode Internal Solitary Waves

Surface Fluxes and Wind-Wave Interactions in Weak Wind Conditions

LAGRANGIAN MEASUREMENTS OF EDDY CHARACTERISTICS IN THE CALIFORNIA CURRENT

Swash Zone Dynamics: Modeling and Data Analysis

Analysis of Mixing and Dynamics Associated with the Dissolution of Hurricane-Induced Cold Wakes

Transcription:

Hydrogen termination following Cu deposition on Si(001) L. A. Baker, A. R. Laracuente,* and L. J. Whitman Naval Research Laboratory, Washington, DC 20375-5342, USA Received 9 September 2004; published 8 April 2005 We describe the surface structures following submonolayer Cu deposition on Si 001 and subsequent hydrogen termination as characterized by scanning tunneling microscopy. Cu adsorption at 870 K results in a characteristic 2 8 island+vacancy structure, as previously reported. In addition, occasional structures are observed attributed to Cu in surface interstitial sites. After H termination, the dominant features of the island+vacancy structure remain, but the size and distribution of the structures are significantly altered. Based on the atomic-scale appearance of both the clean and H-terminated structures, we propose that within the 2 8 island+vacancy structure all surface atoms are Si, with all Cu subsurface, contrary to previous structural models. DOI: 10.1103/PhysRevB.71.153302 PACS number s : 68.37.Ef, 68.43.Fg, 68.35. p The deposition and surface chemistry of metals on silicon is of considerable importance both to current technology and emerging nanoscience. Technologically, the CuuSi 001 interface is of particular significance, with current integrated circuit technology incorporating Cu interconnects formed by techniques such as ion-assisted trench filling. 1 Although Cu is a key component of such interconnects, interdiffusion of Cu from the surface into the bulk of the silicon chip can adversely affect device performance. 2 Therefore, a further understanding of the Cu/ Si interface, Cu contamination in Si, and the potential effects of hydrogen ubiquitous in processing could be significant in combating surface-tobulk diffusion of Cu on Si and in improving the processing of valuable Cu-based interconnects. In addition to the importance of metal-si interactions in current device technologies, such interactions are fundamentally important in the design of nanostructured surfaces. Transition metals, such as Co, Ni, and Cu, are known to react readily with Si to form metal silicides. 3 In the case of Ni, this reactivity has been exploited to create templates for metal deposition with nanometer scale order. 4 In the case of Cu, the high diffusion 5 and reactivity of Cu in Si make this system an intriguing candidate for the construction of similar nanoscale structures. The surface reactivity and stability can be further modulated by chemically functionalizing the surface. For example, hydrogen termination of Si is particularly well understood, both for manipulating the surface chemistry and clarifying the surface structure. 6 The deposition and growth of Cu on Si has been studied previously by a variety of techniques. 7 11 Scanning tunneling microscopy STM has been especially useful in exploring the Cu/Si surface on different crystal faces. 12 17 On the clean Si 001 surface, STM studies by Ichinokawa and co-workers, 16 18 Lindsay and co-workers, 12 and Nogami and co-workers 14 have all shown that the deposition of Cu results in the formation of islands and vacancies. The location of Cu on these surfaces, however, has yet to be resolved. Herein we report an STM analysis of the initial growth of Cu on Si 001 and the effects of subsequent hydrogen termination on the Cu-induced surface structure. Hydrogen termination causes numerous changes in the surface structure, some of which strongly suggest that previous structures proposed for Cu/Si 001 are incorrect. To our knowledge, this work presents the first report of the effects of hydrogen following submonolayer metal depostion on Si 001. Experiments were performed in ultrahigh vacuum UHV base pressure 5 10 11 Torr using Sb-doped Si 001 wafers 0.01 0.02 cm oriented to within 0.1 of 001. Samples were initially degassed in UHV overnight. Clean surfaces were prepared by resistively heating for 30 s at 1500 K, cooling for at least 5 min, and then briefly heating again to 1500 K for 5 s. The heating power was then reduced until the sample temperature desired for deposition was reached. Sample temperatures were determined from the heating power as described previously. 19 Copper was evaporated from a pure Cu wire wrapped in a resistively heated W wire. Samples were held at 870 K during Cu deposition to reproduce the highest deposition temperature previously reported. 14,15 In order to mimic the exposure conditions to be used during H termination, after depositing Cu the sample was post annealed in UHV for 10 min at 615 K. A total coverage of 0.6 to 1.25 ML of Cu 1 ML=6.78 10 14 cm 2 was deposited at a rate of 0.6 ML min 1 as measured with a quartz crystal microbalance. It is important to note that our substrate deposition temperature is at the high end of those reported previously 720 873 K, 14 16,18 which increases Cu diffusion into the bulk, lowering the effective Cu coverages on our surfaces compared to those dosed at lower substrate temperatures. To terminate a surface with H, a sample was held at 615 K and exposed for 2 min to atomic H produced by a hot tungsten filament located 1 cm away from the sample and a background H 2 pressure of 1 10 5 Torr. After H exposure, each sample was annealed again in UHV for 10 min at 615 K. On clean surfaces etching is not observed under these conditions. 20 All STM images were collected at room temperature starting 1 h after preparing the surface. All results shown are filled-state, gray-scale images collected at 2 3 V with constant currents of 0.1 1 na. A typical STM image of Cu 1.2 ML deposited at 870 K on clean Si 001 is shown in Fig. 1, where a large number of islands bright protrusions and vacancies dark recesses are observed. At higher Cu coverage an overall c 8 8 reconstruction occurs. 14 The atomic-scale structure of the island+ vacancy structure observed at low coverage is shown in Fig. 1 b, where a dashed box highlights one 2 8 unit cell. The islands are always oriented normal to the direction of the dimer rows and typically consist of three 1098-0121/2005/71 15 /153302 4 /$23.00 153302-1

Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. 1. REPORT DATE 2004 2. REPORT TYPE 3. DATES COVERED 00-00-2004 to 00-00-2004 4. TITLE AND SUBTITLE Hydrogen termination following Cu deposition on Si(001) 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) 5d. PROJECT NUMBER 5e. TASK NUMBER 5f. WORK UNIT NUMBER 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) Naval Research Laboratory,4555 Overlook Avenue SW,Washington,DC,20375 8. PERFORMING ORGANIZATION REPORT NUMBER 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR S ACRONYM(S) 12. DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release; distribution unlimited 13. SUPPLEMENTARY NOTES 14. ABSTRACT 11. SPONSOR/MONITOR S REPORT NUMBER(S) 15. SUBJECT TERMS 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF ABSTRACT a. REPORT b. ABSTRACT c. THIS PAGE Same as Report (SAR) 18. NUMBER OF PAGES 4 19a. NAME OF RESPONSIBLE PERSON Standard Form 298 (Rev. 8-98) Prescribed by ANSI Std Z39-18

protrusions, each centered on a dimer row. Dimer vacancies, one of which is labeled DVVVVD in Fig. 1 b, are observed between and at the end of the islands. Previously proposed structural models have assumed that these vacancies consist of two missing dimers. 12,14 Our atomic resolution images clearly show that these vacancies actually include four missing dimers, invalidating all previous structural models. We believe this discrepancy is a consequence of image resolution; a dull tip can make islands seem wider and make vacancies seem smaller than actual. In addition, longer islands with more than three bright features are occasionally but rarely observed, forming 2 n features. In these 2 n islands, additional bright protrusions reside between the dimer rows not shown. Additional vacancies observed on the surface after Cu deposition are very similar to the ones observed on clean Si 001. An analysis of the observed vacancy type and frequency, along with a comparison of the dimer vacancies observed on clean Si 001 from Wang et al. 21 is presented in Table I. Interestingly, one specific pattern is observed with high frequency on our Cu/Si 001 surfaces, an example of which is labeled VDVVDV on the right side of Fig. 1 b. The pattern consists of a vacancy, a dimer, two vacancies, a dimer and a vacancy. After Cu deposition this pattern is observed 30 times more frequently than on the clean surface, indicating it is a Cu induced structure. Transition metal incorporation in the silicon surface often leads to distinctive dimer vacancy complex formations. A classic example is the case of Ni on Si 001, in which concentrations of as little as 5% create arrays of line vacancies. 22,23 We observe an unusual structure on this surface that does not appear to have been noted in the past. Two examples are highlighted with arrows in Fig. 1 b, and the general structure can be determined from Fig. 1 c. The new structure consists of two round protrusions centered in the middle of the dimer row, with each protrusion offset from the dimer position by half a dimer. Note that copper impurities in bulk Si reside either at substitutional or interstitial sites. It is reasonable to assume that substitutional Cu will produce features aligned with the Si lattice. In contrast, the unusual location of these protrusions with respect to the lattice suggests they are interstitial Cu atoms at the surface. The effects that H termination has on the Cu/Si 100 island+vacancy surface are shown in Fig. 2. There are several apparent differences produced by the hydrogen termination. First, the average island length increased. Whereas islands on the Cu/Si 001 surface preferentially consist of three bright maxima and a terminating vacancy on each end, after hydrogen termination longer islands are observed Fig. 3. Another difference is the position of the vacancies with respect to the islands the islands are no longer preferentially terminated by vacancies. On average, one terminating vacancy of every island appears to have diffused away from the island. The remaining vacancies show no recurring pattern, randomly occurring as on a clean surface. Figure 2 b shows two common island+vacancy structures observed after hydrogen termination. Qualitatively, the islands look strikingly similar to Si islands created during Si homoepitaxy, where islands composed of 5 9 consecutive dimers are often observed. On the atomic scale, all features are identical FIG. 1. a Filled-state STM image of 1.2 ML Cu deposited on Si 001 at 870 K. b The dashed box highlights a 2 8 unit cell. A four dimer-long terminating vacancy is labeled DVVVVD, and a frequently observed defect structure is labeled VDVVDV. The arrows point to features attributed to interstitial Cu atoms. c Highresolution image of two proposed interstitial Cu atoms. The lines highlight the location of adjacent dimers. to the symmetric dimers observed on monohydrideterminated Si 001, including the underlying structures within the vacancy islands. 153302-2

TABLE I. Notation, pattern, and observation of dimer vacancy defects from Ref. 21, and this study. The fraction n is the ratio of the number of defects N to the total number of observed dimers. a Does not include vacancies associated with islands. Previous reports of Cu deposition on Si 001 are somewhat inconsistent. Ichinokawa and co-workers proposed a Cu/Si 001 phase diagram with mixed 2 1, 2 2, 6 2, and c 12 10 phases for 0.1 20 ML Cu coverage and 670 870 K substrate temperatures. In the interpretation of their STM results, the bright islands observed at low coverage 0.1 ML were attributed to Si ejected from a dimer row by the deposition of Cu, which they proposed to be bound under the island. In contrast, STM studies by both Nogami and co-workers and Lindsay and co-workers reported only a c 8 8 phase. Lindsay and co-workers attributed the presence of a c 8 8 phase on Si 001 to Cu contamination, as identified by secondary ion mass spectrometry. No structural model was ever proposed for the actual position of the Cu contaminants on the reconstructed surface. Nogami and coworkers, who deposited submonolayer 0.05 0.20 ML Cu on Si 001 at 720 870 K, proposed a structural model with Cu residing at the island maxima and a two dimer vacancy, in contradiction to the ejected Si model of Ichinokawa and co-workers. We observed a c 8 8 surface similar to the one reported previously by Nogami and Lindsay, but believe that Nogami s model is incorrect. Copper normally occupies an interstitial site when present in bulk Si; however, theoretical simulations have shown there is a significant energetic advantage 2.71 ev for interstitial Cu + to react with Si vacancies and form substitutional Cu. 24 Therefore, it is probable that Cu initially goes in the surface interstitially, and then reacts through surface defects to occupy substitutional sites, ejecting Si atoms in the process which then become the observed homoepitaxial-like islands. This mechanisms is further supported by reports that Si deposition onto defect-rich surfaces results in similar island structures. 25 Ab initio calculations by Estreicher show that the Cu- Si bond is less stable than the Si-H bond by approximately 2 ev. 26 Therefore, H passivation should have a significant effect on any Cu present at the surface. If the islands were made out of Cu, as suggested by Nogami, one would expect FIG. 2. Filled-state STM images of 1.2 ML Cu deposited on Si 001 at 870 K, followed by H termination at 615 K. b The atomic-scale structure of a typical island. 153302-3

FIG. 3. Distribution of island lengths relative to dimer row width, 0.768 nm on Si 001 following Cu deposition filled bars and after subsequent hydrogen termination dashed bars. hydrogen termination to displace the Cu from the surface and therefore reduce or eliminate the islands. Although STM cannot unambiguously determine the position of Cu in the reconstruction, in light of the images obtained following hydrogen termination, where only H-terminated Si-like features are observed, we propose that all Cu is subsurface after H termination. We believe the results to date are best explained by a model where each Cu atom initially enters the Si lattice interstitially, then ejects a Si atom from the lattice to occupy a substitutional site either underneath a Si island or an adjacent vacancy. After H passivation, which changes the surface energies, surface Si atoms no longer bond to the Cu, driving the Cu into the bulk and allowing the islands and vacancies to reconstruct into the structures observed. Our essential structural results can be seen in composite in Fig. 4. Clean Si 001 exposed to Cu has a distinct island +vacancy structure which becomes a c 8 8 reconstruction at high coverage. We propose that Cu enters the surface at interstitial sites and relaxes into substitutional sites ejecting Si atoms in the process. Most of the near-surface Cu is underneath the Si islands and within the vacancies, with the FIG. 4. Composite of the surface structures observed in this study. a Clean Si 001, b Cu-dosed Si 001, c hydrogenterminated Si 001, d hydrogen-terminated, Cu-dosed Si 001. remaining residual in interstitial sites. When the Cu/Si 001 surface is hydrogen terminated, the island+ vacancy structure is retained but the size and distribution of the structures are altered. This change is in contrast to monohydrideterminated Si 001, where surfaces with extremely low defect densities are routinely achieved. Based on the atomicscale appearance of the H-terminated Cu/Si 001, we propose that all surface atoms are H-terminated Si, and that all Cu is actually subsurface. Note that such a surface structure Si 001 :H with a relatively narrow distribution of defects could make an interesting template for nucleating nanostructures. We are presently investigating how these defects affect the subsequent nucleation and growth of metal monolayers. This work was funded by the Office of Naval Research and NRL/NRC L.A.B.. *Electronic address: laracuente@nrl.navy.mil 1 O. R. Monteiro, J. Vac. Sci. Technol. B 17, 1094 1999. 2 A. A. Istratov and E. R. Weber, J. Electrochem. Soc. 149, G21 2002. 3 E. R. Weber, Appl. Phys. A: Solids Surf. 30, 1 1983. 4 J.-Z. Wang et al., Phys. Rev. B 65, 235303 2002. 5 A. A. Istratov et al., Phys. Rev. Lett. 81, 1243 1998. 6 K. Oura et al., Surf. Sci. Rep. 35, 1 1999. 7 S. Guillet et al., Solid State Commun. 277, L52 1992. 8 S. Guillet et al., Appl. Surf. Sci. 65/66, 742 1993. 9 P. N. Hai et al., Phys. Rev. B 56, 4620 1997. 10 M. Hanbucken and G. Le Lay, Surf. Sci. 168, 122 1986. 11 P. Mathiez et al., Surf. Sci. 168, 158 1986. 12 P. W. Murray et al., Phys. Rev. B 54, 13 468 1996. 13 T. Yasue and T. Koshikawa, Surf. Sci. 377 379, 923 1997. 14 B. Z. Liu et al., Surf. Sci. 453, 137 2000. 15 H. Itoh et al., Surf. Rev. Lett. 5, 747 1998. 16 T. Ichinokawa et al., J. Vac. Sci. Technol. B 12, 2070 1994. 17 T. Ikeda et al., Surf. Sci. 336, 76 1995. 18 T. Ichinokawa et al., Surf. Sci. 241, 416 1991. 19 A. Laracuente and L. J. Whitman, Surf. Sci. 476, L247 2001. 20 A. R. Laracuente and L. J. Whitman, Surf. Sci. 545, 70 2003. 21 J. Wang et al., Phys. Rev. B 47, 10 497 1993. 22 H. Niehus et al., J. Microsc. 152, 735 1988. 23 J.-Y. Koo et al., Phys. Rev. B 52, 17 269 1995. 24 S. K. Estreicher, Phys. Rev. B 60, 5375 1999. 25 H. W. Liu et al., Nanotechnology 13, 729 2002. 26 S. K. Estreicher, Physica B 273 274, 424 1999. 153302-4