Nanoscale Research Letters manuscript No. (will be inserted by the editor) Phase transition on the Si(001) clean surface prepared in UHV MBE chamber A study by high resolution STM and in situ RHEED 1 L. V. Arapkina · V. A. Yuryev · K. V. Chizh · V. M. Shevlyuga · 1 M. S. Storojevyh · L. A. Krylova 0 2 n a J Received: date/Accepted: date 8 1 Abstract The Si(001) surface deoxidizedby short an- Keywords Silicon · Surface reconstruction· Scanning ] nealing at T ∼925◦C in the ultrahigh vacuum molecu- tunnelling microscopy · Reflected high energy electron i c larbeamepitaxychamberhasbeenin situinvestigated diffraction · Clean surface preparation s - byhighresolutionscanningtunnellingmicroscopy(STM) l r andreflectedhighenergyelectrondiffraction(RHEED). mt RHEED patterns corresponding to (2×1) and (4×4) 1 Introduction structureswereobservedduringsampletreatment.The . at (4×4)reconstructionaroseatT >600◦Cafteranneal- Investigationsof cleansilicon surfacespreparedin con- m ing. The reconstruction was observed to be reversible: ditions of actual technological chambers are of great - the (4×4)structure turnedinto the (2×1)one atT ? interest due to the industrial requirements to operate d 600◦C, the (4×4) structure appeared again at recur- on nanometer and subnanometer scale when designing n o ring cooling. The c(8×8) reconstruction was revealed future nanoelectronicdevices [1]. In the nearestfuture, c by STM at room temperature on the same samples. A the sizes of structural elements of such devices will be [ fractionofthesurfaceareacoveredbythec(8×8)struc- close to the dimensions of structure features of Si(001) 2 ture decreasedas the sample cooling rate was reduced. surface, at least of its high-order reconstructions such v The (2×1) structure was observed on the surface free as c(8×8). Most of researches of the Si(001) surface 9 of the c(8× 8) one. The c(8 ×8) structure has been have thus far been carried out in specially refined con- 0 evidenced to manifest itself as the (4×4) one in the ditions which allowed one to study the most common 9 3 RHEED patterns. A model of the c(8 × 8) structure types of the surface reconstructions such as (2 × 1), . formationhasbeenbuiltonthebasisofthe STMdata. c(4×4),c(4×2)orc(8×8)[2,3,4,5,6,7,8,9,10,11,12,13, 9 0 Originofthehigh-orderstructureontheSi(001)surface 14]. Unfortunately, no or very few papers have thus far 0 and its connection with the epinucleation phenomenon been devoted to investigations of the Si surface which 1 are discussed. is formed as a result of the wafer cleaning and deoxi- : v dation directly in the device manufacturing equipment i X [14].ButanyonewhodealswithSi-basednanostructure r The research was supported by the Science and Innovations engineeringandthedevelopmentofsuchnanostructure a Agency of RFunder the State Contract No. 02.513.11.3130 and formationcyclescompatiblewithsomestandarddevice theEducationAgencyofRFundertheStateContractNo.Π2367. manufacturingprocessesmeetsthechallengingproblem LarisaV.Arapkina of obtaining the clean Si surface within the imposed A. M. Prokhorov General Physics Institute of RAS, 38 Vavilov technological restrictions which is one of the key ele- Street,Moscow,119991, Russia ments of the entire structure formationcycle [1,15,16]. Tel.:+7-499-5038318 The case is that the ambient in technological ves- E-mail:[email protected] sels such as molecular beam epitaxy (MBE) chambers VladimirA.Yuryev is usually not as pure as in specially refined ones de- Tel.:+7-499-5038144 Fax:+7-499-1350356 signed for surface studies. There are many sources of E-mail:[email protected] surfacecontaminantsintheprocesschambersincluding 2 materialsofwaferheatersorevaporatorsofelementsas atoms deposition on the Si(001)-(2×1) surface [7,10]; well as foreignsubstances used for epitaxy anddoping. similar structures were found to arise due to various In addition, due to technological reasons the tempera- treatments and low-temperature annealing of the orig- turetreatmentsapplicablefordevicefabricationfollow- inal Si(001)-(2×1) surface without deposition of any ing the standardprocessessuchasCMOSoftencannot foreign atoms [4,5,6]. Data of the STM studies of the be as aggressive as those used for surface preparation Si(001)-c(8×8) surface were presented in Refs. [5,10]. in the basic experiments. Moreover, the commercially ItmaybesupposedontheanalogywiththeSi(001)- available technological equipment sometimes does not c(4×4)reconstruction[12,31,32,33,34,35]thatthepres- enablethewishfulannealingofSiwafersatthetemper- ence of impurity atoms on the surface as well as in the ature of ∼ 1200◦C even if the early device formation subsurface regions is not the only reason of formation stage allows one to heat the wafer to such a high tem- of reconstructions different from the (2×1) one, and perature. Nevertheless, the technologist should always the conditions of thermal treatments should be taken be convinced that the entirely deoxidized and atomi- intoaccount.Theresultsofexplorationofeffectofsuch cally clean Si surface is reliably and reproducibly ob- factor as the rate of sample cooling from the anneal- tained. ing temperature to the room one on the process of the A detailed knowledge of the Si surface structure c(8×8) reconstruction formation are reported in the whichisformedintheaboveconditions—itsreconstruc- present article. It is shown by means of RHEED that tion,defectiveness,finestructuralpeculiarities,etc.—is thediffractionpatternscorrespondingtothe(2×1)sur- of great importance too because this structure may af- face structure reversibly turn into those corresponding fect the properties of nanostructured layers deposited to the c(8×8) one depending on the sample tempera- on it. For instance, the Si surface structure may affect ture,andapointofthisphasetransitionisdetermined. themagnitudeandthedistributionofthesurfacestress BasedontheSTMdataamodelofthec(8×8)structure of the Ge wetting layer on nanometer scale when the formation is brought forward. Ge/Si structure is grown, which in turn affect the Ge nanoclusternucleationandeventuallythe propertiesof 2 Methods and equipment quantum dot arraysformedon the surface [1,16,17,18, 19,20,21,22,23,24,25,26,27,28,29,30]. The experiments were made using an integrated ultra- Thus, it is evident from the above that the control- high-vacuum (UHV) system [27] based on the Riber lable formation of the clean Si(001) surface with the EVA32molecularbeamepitaxychamberequippedwith prescribed parameters required for technologicalcycles the Staib Instruments RH20 diffractometer of reflected ofnanofabricationcompatiblewiththestandarddevice high energy electrons and coupled through a transfer manufacturingprocessesshouldbeconsideredasanim- line with the GPI 300 UHV scanning tunnelling mi- portant goal, and this article presents a step to it. croscope [36,37,38]. This instrument enables the STM In the present paper, we report the results of in- study of samples at any stage of Si surface preparation vestigation of the Si(001) surface treated following the and MBE growth. The samples can be serially moved standardprocedureofSiwaferpreparationfortheMBE into the STM chamber for the analysis and back into growth of the SiGe/Si(001) or Ge/Si(001) heterostruc- the MBE vessel for further treatments as many times tures. A structure arising on the Si(001) surface as a as required never leaving the UHV ambient. RHEED result of short high-temperature annealing for SiO2 re- experiments can be carried out in situ, i.e. directly in movalisexplored.Itiswellknownthatsuchexperimen- the MBE chamber during the process. tal treatments favor the formation of nonequilibrium Samples for STM were 8×8 mm2 squares cut from structuresonthesurface.Themoststudiedofthemare the specially treated commercial B-doped CZ Si(100) presentlythe(2×1)andc(4×4)ones.Thisworkexper- wafers (p-type, ρ = 12 Ωcm). RHEED measurements ′′ imentally investigates by means of scanning tunneling were carried out at the STM samples and similar 2 ′′ microscopy (STM) and reflected high energy electron wafers;the2 sampleswereinvestigatedonlybymeans diffraction (RHEED) the formation and atomic struc- ofRHEED.Afterchemicaltreatmentfollowingthestan- ture of the less studied high-order c(8×8) (or c(8×n) dard procedure described elsewhere [1,39] (which in- [14,15,16]) reconstruction. Observations of this recon- cluded washing in ethanol, etching in the mixture of struction have already been reported in the literature HNO3 and HF and rinsing in the deionized water),the [4,5,6,10]butthereisnoclearcomprehensionofcauses samples were placed in the holders. The STM samples ofitsformationasthestructureslookinglikethec(8×8) weremountedonthemolybdenumSTMholdersandin- one appear after different treatments: The c(8×8) re- flexiblyclampedwiththetantalumfasteners.TheSTM construction was observed to be a result of the coper holders were placed in the holders for MBE made of 3 (a) Fig.1 Adiagramofsamplecoolingafterthethermaltreatment at925◦CmeasuredbyIRpyrometer;coolingratesareasfollows: ∼0.17◦C/sor“slowcooling”oftheSTMsamples(1);∼0.4◦C/s or“quenching” oftheSTMsamples(2)and2′′ wafers(3). ′′ molybdenumwithtantaluminserts.The2 waferswere inserted directly into the standard molybdenum MBE holdersand didnot havesohardfastening as the STM samples. (b) Thereuponthe sampleswereloadedintotheairlock andtransferredintothepreliminaryannealingchamber Fig. 2 STMimageoftheSi(001)surfacewiththeresidualsili- where outgassed at ∼ 600◦C and ∼ 5×10−9 Torr for conoxide(−1.5V,150pA), annealingat∼925◦C for∼2min. (a),theimageisinverted:darkareascorrespondtotheoxide,the about 6 hours. After that the samples were moved for lighterareasrepresentthedeoxidizedsurface;STMimageofthe finaltreatmentanddecompositionoftheoxidefilminto cleanSi(001)surface(+1.9V,70pA)withtheFouriertransform the MBE chamber evacuated down to ∼ 10−11Torr. pattern shown inthe insert, annealing at ∼925◦C for ∼3min. (b)[14]. There were two stages of annealing in the process of sample heating—at ∼ 600◦C for ∼ 5min and at ∼ 800◦C for ∼ 3min [1,14,27]. The final annealing was ent elevated temperatures in the process of the sample carried out at ∼ 925◦C.1 Then the temperature was treatment and at room temperature after cooling. The rapidly lowered to ∼ 850◦C. The rates of the further STM samples were always explored by RHEED before coolingdowntothe roomtemperaturewere∼0.4◦C/s moving into the STM chamber. (referredto as the “quenching”mode ofboth the STM The STM tips were ex situ made of the tungsten samplesand2′′ wafers)or∼0.17◦C/s(calledthe“slow wire and cleaned by ion bombardment [40] in a special cooling” mode of only the STM samples) (Fig. 1). The UHV chamber connected to the STM chamber. The pressure in the MBE chamber increased to ∼2×10−9 STM images were obtained in the constant tunnelling Torr during the process. current mode at room temperature. The STM tip was In both chambers, the samples were heated from zero-biased while the sample was positively or nega- the rear side by radiators of tantalum. The tempera- tively biased when scanned in empty or filled states ture was monitored with the IMPAC IS12-Si pyrome- imaging mode. terwhichmeasuredtheSisampletemperaturethrough The STM images were processed afterwords using chamber windows. The atmosphere composition in the the WSxM software [41]. MBE chamber was monitored using the SRS RGA-200 residual gas analyser before and during the process. After cooling, the STM samples were moved into 3 Experimental findings the STM chamber in which the pressure did not ex- ceed 1×10−10 Torr. RHEED patterns were obtained Fig.2 demonstratesSTMimagesofthe Si(001)surface for all samples directly in the MBE chamber at differ- afterannealingat∼925◦Cofdifferentduration.Fig.2a depicts the early phase of the oxide film removal; the 1 The samples were heated over 920◦C about a half of the annealing duration is 2 min. A part of the surface is final annealing time; a diagram of the thermal processing and someadditionaldetailscanbefoundinRef.[27]. still oxidized: the dark areas in the image correspond 4 5 (a) (b) (c) Fig. 3 Empty (a) and filled (b) state images of the same region on the Si(001) surface (+1.7 V, 100 pA and −2.0 V, 100 pA). Positionsofextremesoflinescanprofiles(c)matchexactlyfortheempty(1)andfilled(2)statedistributionsalongthecorresponding linesintheimages(a)and(b). to the surface coated with the oxide film. The lighter areas correspond to the purified surface. A structure of ordered “rectangles” (the grey features) is observed on the deoxidizes surface. After longer annealing (for 3 min.) and quenching (Fig. 1), the surface is entirely purified of the oxide (Fig. 2b). It consists of terraces separated by the SB or SA monoatomic steps with the heightof ∼1.4 ˚A [3]. Eachterraceis composedofrows running along [110] or [110] directions. The surface re- constructionisdifferentfromthe(2×1)one.Theinsert (a) of Fig. 2b demonstrates the Fourier transform of this image which corresponds to the c(8×8) structure [5]: ReflexesoftheFouriertransformcorrespondtothedis- tance ∼31 ˚A in both [110]and [110] directions. So the revealed structure have a periodicity of ∼ 31 ˚A that correspondsto8 translationsa onthe surfacelattice of Si(001)alongthe<110>directions(a=3.83˚Aisaunit translation length). Rows consisting of structurally ar- rangedrectangularblocksareclearlyseenintheempty ◦ (b) state STM image (Fig. 2b). They turn by 90 on the neighbouring terraces. Fig. 4 Reflected high energy electron diffraction patterns ob- served in the [010] (a) and [110] (b) azimuths; electron energy Fig. 3 demonstrates the empty and filled state im- was9.8keVand9.3keV,respectively. ages of the same surface region. Each block consists of two lines separated by a gap. This fine structure of tances on the surface corresponding to the reflex posi- blocks is clearly seen in the both pictures (a) and (b) tions in the diffraction pattern were calculated accord- butitsimagesaredifferentindifferentscanningmodes. ingtoRef.[42].Thederivedsurfacestructureis(4×4). Acharacteristicpropertymostclearlyseeninthe filled One sample showed the RHEED patterns correspond- state mode (Fig. 3b) is the presence of the brightness ingtothe(2×1)structure[42]afterthesametreatment maxima on both sides of the lines inside the blocks. though. These peculiar features are described below in more Temperature dependences of the RHEED patterns detail. Fig. 3c shows the profiles of the images taken in the [110] azimuth were investigated during sample alongthe white lines.Extremepositions ofbothcurves heating andcooling.It was found that the reflexes cor- are well fitted. Relative heights of the features outside responding to 2a were distinctly seen in the RHEED andinsidetheblockscanbeestimatedfromtheprofiles. patterns during annealing at ∼ 925◦C after 2 minutes Fig. 4 demonstratestypicalRHEEDpatterns taken of treatment. The reflexes corresponding to 4a started at room temperature from the STM sample annealed to appear during sample quenching and became defi- for 3 min. with further quenching. Characteristic dis- nitely visible at the temperature of ∼ 600◦C; a weak 5 riod of 8a remains in such rows.Two adjacent terraces aredesignatedinFig.5abyfigures‘1’and‘3’.Arowof “rectangles”markedas‘2’issituatedontheterrace‘3’; ithasthesameheightastheterrace‘1’.Thefilledstate image, which is magnified in comparison with the for- mer one, is given in Fig. 5b. A part of the surface free of the “rectangles” is occupied by the (2 × 1) recon- struction. Images of the dimer rows with the resolved Si atoms are marked as ‘B’ in Fig. 5b. The “rectan- gles” are also seen in the image (they are marked as ‘A’) as well as single defects: dimerized Si atoms (‘C’) andchaoticallylocatedonthesurfaceaccumulationsof (a) severaldimers.Mostofthesedimersareorientedparal- leltodimersofthelowersurfaceandlocatedstrictlyon the dimer row. Note that influence of the cooling rate onthesurfacestructurewasobservedbytheauthorsof Ref.[6]:whenthesamplecoolingratewasdecreasedthe surface reconstruction turned from c(8×8) to c(4×2) whichwasconsideredasthederivativereconstructionof the (2×1) one transformedbecause of dimer buckling. (b) Fig. 6 presents the STM images obtained for the Fig.5 STMimagesofthecleanSi(001)surfacepreparedinthe samples cooled in the quenching mode but containing slowcoolingmode:(a)thesurfacemainlycoveredbythe(2×1) areas free of “rectangles”. The images (a) and (b) of structure, +2.0V,100pA,‘1’and‘3’areterraces, theheightof the same place on the surface were obtained serially therow‘2’coincideswiththeheightoftheterrace‘1’;amagnified image taken with atomic resolution (b), −1.5 V, 150 pA, ‘A’ is one by one. We managed to image the surface struc- the“rectangle”,‘B’marksthedimerrowscomposingthe(2×1) ture between the areas occupied by the “rectangle” structure(separateatomsareseen),‘C’showsstructuraldefects, rows, but only in the filled state mode (see the insert i.e.thedimersof theuppermost layeroriented alongthe dimers ofthelower(2×1)rows(b). at Fig. 6b). Like in Fig. 5b this structure is seen to be formed by parallel dimer rows going 2a apart. The direction of these rows is perpendicular to the direc- (4×4) signal started to arise at ∼ 525◦C if the sam- tion of the rows of “rectangles”. The height difference ple was cooled slowly (Fig. 1). At the repeated heating of the rows of “rectangles” and the (2×1) rows is 1 fromroomtemperature to 925◦C,the (4×4)structure monoatomic step (∼1.4˚A). We did not succeed to ob- disappearedat∼600◦Cgivingplacetothe(2×1)one. tainagoodenoughimageofthesesubjacentdimerrows The (4×4) structure appearedagainat ∼600◦ during in the empty state mode. It should be noted also that recurring cooling. positions of the “rectangles” are always strictly fixed ′′ relative to the dimer rows of the lower layer: they oc- The RHEEDpatterns obtainedfrom2 samples al- wayscorrespondedtothe(2×1)reconstruction.Diffrac- cupy exactly three subjacent dimer rows. It also may be seen in the STM images presented in Refs. [5,10]. tion patterns for the STM sample which was not hard fastenedtotheholdercorrespondedtothe(4×4)struc- tureafterquenching(STMmeasurementswerenotmade 3.1 Fine structure of the observed reconstruction for this sample). Effectsofannealingdurationandcoolingrateonthe Let us consider the observed structure in detail. cleansurfacestructurewerestudiedbySTM.Itwases- A purified sample surface consists of monoatomic tablished that increase of annealing duration to 6 min. steps. Following the nomenclature by Chadi [3], they did not cause any changes of the surface structure. On aredesignatedasSA andSB inFig.2b.Eachterraceis the contrary, decrease of the sample cooling rate dras- composed by rows running along the [110] or [110] di- tically changes the structure of the surface. The STM rections.Eachrowconsistsofrectangularblocks(“rect- images of the sample surface for the slow cooling mode angles”). They may be regarded as surface structural (Fig. 1) are presented in Fig. 5. The difference of this units as they are present on the surface after thermal surface from that of the quenched samples (Fig. 2b) is treatment in any mode, irrespective of a degree of sur- that only a few rows of “rectangles” are observed on face coverage by them. Reflexes of the Fourier trans- it. The order of the “rectangle” positions with the pe- form of the picture shown in Fig. 2b correspondto the 6 (a) (a) (b) Fig.7 STMemptystateimagesoftheSi(001)surface;ac(8×8) (b) unitcellismarkedbythewhiteboxinimage(a)(+1.9V,50pA), Fig.6 Empty(a)andfilled(b)stateimagesofthesameregion distancesbetweentherowsmarkedby‘A’and‘B’equal3aand4a on the Si(001) surface (+2.0 V, 100 pA and −2.0 V, 100 pA); (thatcorrespondstoc(8×6)andc(8×8)structures,respectively), aninsertat (b) shows theimage ofthe (2×1)surfaceobtained twolong“rectangles”anddivacanciesarisingintheadjacentrows betweentherowsof“rectangles”. aremarked by ‘L’and ‘V’,respectively; a row wedging between two rows (‘W’) and lost blocks (‘LB’) are seen in (b) (+1.6 V, 100pA). distances ∼ 31 and ∼ 15˚A in both [110] and [110] di- rections. Hence the structure revealed in the long shot seems to have a periodicity of ∼31˚A that corresponds point defects are observedin this image. In addition, a to 8 translations a on the surface lattice of Si(001). It row wedging in between two rows and separating them resembles the Si(001)-c(8×8) surface [5]. Reflexes cor- by an additional distance a is seen in the centre of the responding to the distance of ∼ 15˚A (4a) arise due to upper side of the picture (‘W’). The total distance be- theperiodicityalongtherows.STMimagesobtainedat tween the wedged off rows becomes 5a. higher magnificationsgive anevidence that the surface Henceitmaybeconcludedthatthe orderandsome appears to be disordered, though. periodicitytakeplaceonlyalongthe rows—disordering Fig. 7 shows the magnified images of the investi- ofthe c(8×8)structureacrossthe rowsisrevealed(we gatedsurface.Therowsoftheblocksareseentobesitu- often refer to this structure as c(8×n)). atedatvaryingdistancesfromoneanother(hereinafter, The block length can possess two values: ∼ 15˚A thedistancesaremeasuredbetweencorrespondingmax- (4a)and∼23˚A(6a).Distancesbetweenequivalentpo- ima of features). A unit c(8×8) cell is marked with a sitionsoftheadjacentshortblocksintherowsare8a.If square box in Fig. 7a. The distances between the ad- the long block appears in a row, a divacancy is formed jacent rows of the rectangles are 4a in such structures in the adjacent row to restore the checkerboard order (‘B’ in Fig. 7a). The adjacent rows designated as ‘A’ of blocks. Fig. 7a illustrates this peculiarity. The long are 3a apart (c(8×6)). blockismarkedas‘L’,thedivacancyarisenintheadja- A structure with the rows going at 4a apart is pre- centrow is lettered by ‘V’. In addition, the long blocks sented in Fig. 7b. The lost blocks (‘LB’) that resemble werefoundtohaveonemorepeculiarity.Theyhaveex- 7 ima) of zigzag chains (Fig. 8b). The distances between the maxima are ∼4˚A along the rows. The presented STM data are interpreted by us as a structure composed by Si ad-dimers and divacancies. 4 Discussion 4.1 Structural model The above data allow us to bring forward a model of theobservedSi(001)surfacereconstruction.Themodel isbasedonthefollowingassumptions:(i)theoutermost (a) surfacelayerisformedbyad-dimers;(ii)theunderlying layer has a structure of (2×1); (iii) every rectangular block consists of ad-dimers and divacancies a number of which controls the block length. Fig. 9a shows a schematic drawing of the c(8×8) structure (a unit cell is outlined). This structure is a basic one for the model brought forward. The elemen- tary structural unit is a short rectangle. These blocks form raised rows running vertically (shown by empty circles). Smaller shaded circles show horizontal dimer rows of the lower terrace. The rest black circles show bulk atoms. Each “rectangle” consists of two dimer pairs separated with a dimer vacancy. The structures (b) on the Si(001) surface composed of close ad-dimers are Fig.8 Emptystate(a)andfilledstate(b)imagesoftheSi(001)- believed to be stable [6,13] or at least metastable [43]. c(8×n) surface (+1.7 V, 150 pA, and −2.2V, 120 pA). Corre- Inourmodel,apositionofthe“rectangles”isgoverned sponding schematic drawings of the surfacestructure aresuper- by the location of the dimer rows of the (2×1) struc- imposedonbothpictures.Thelightercircleisthehigherthecor- responding atom is situated inthe surfacestructure. The dimer ture of the underlying layer. The rows of blocks are al- buckling is observed in the filled state image, which is reflected ways normal to the dimer rows in the underlying layer in the drawing by larger open circles representing higher atoms toforma correctepiorientation[43].Everyrectangular ofthetiltedSidimersoftheuppermostlayerofthestructure. block is bounded by the dimer rows of the underly- ing layer from both short sides. Short sides of blocks form non-rebonded SB steps [3] with the underlying tra maxima in their central regions. The maxima are substrate(see Fig.5bandthree verticallyrunning (the not so pronounced as the main ones but nevertheless very left) rows of “rectangles” in Fig. 7a). they are quite recognizable in the pictures (Fig. 7a). Fig.9bdemonstratesthesamemodelforthecaseof Fig.8presentsmagnifiedSTMimagesofthe blocks thelongrectangle.Thisblockisformedduetothepres- (“shortrectangles”).Theimagesobtainedintheempty- enceofanadditionaldimerinthemiddleoftherectan- state (Fig. 8a) and filled-state (Fig. 8b) modes are dif- gle.Thestructureconsistingofonedimerismetastable ferent. In the empty-state mode, short blocks look like [6,13], so this type of blocks cannot be dominating in two lines separated by ∼ 8˚A (the distance is mea- thestructure.Eachlongblockisboundedonbothshort sured between brightness maxima in each line). It is a sides by the dimer rows of the underlying terrace, too. maximum measured value which can lessen depending The presence of the long rectangle results in the for- on scanning parameters. Along the rows, each block is mation of a dimer-vacancy defect in the adjacent row; formed by two parts. The distance between the bright- this is shown in Fig. 9b—the long block is drawn in nessmaximaofthesepartsis∼11.5˚A(orsomegreater the middle row, the dimer vacancy is present in the depending on scanning parameters). In the filled-state lastleft row.Accordingto our STM datathe surface is mode, the block division into two structurally identi- disorderedinthedirectionperpendiculartotherowsof cal parts remains. Depending on scanning conditions, theblocks.Thedistancesbetweentheneighboringrows each part looks like either bright coupled dashes and may be less than those in the c(8×8) structure.Hence blobs (Figs. 3b and 6b) or two links (brightness max- the structure presented in this paper may be classified 8 (a) (b) (c) Fig. 9 Aschematic drawingof thec(8×n)structure: c(8×8) withthe shortblocks (a), aunitcell isoutlined; the samestructure withthelongblock(b); c(8×6)structure(c). as follows: in the Fourier transform of the filled state image, reflexes corresponding to the distance of 8a are absentinthe[110]and[110]directions,whereasthere- flexescorrespondingto4aand2aarepresent(itshould be noticed that the image itself resembles that of the (4×4) reconstructed surface). If an empty state image is not available, it might be concluded that the (4×4) structure is arranged on the surface. An explanation of this observation is simple. Main contribution to the STMimage is made by ad-dimers situated onthe sides Fig. 10 TheSi(001)-c(8×8)surfacereciprocallattice. ofthe“rectangles”,i.e.ontopsoftheunderlyingdimer rows.Accordingtocalculationsmade,e.g.,inRefs.[44, as c(8 × n) one. Fig. 9c demonstrates an example of 45]dimers locatedin sucha wayarecloserto the STM such structure—the c(8×6) one. tip andlook in the images brighterthan those situated In Fig. 8, the presented structure is superimposed in the troughs. Hence, it may be concluded that the on STM images of the surface. The filled state image RHEED (4×4) pattern results from electron diffrac- (Fig. 8b) reveals dimer buckling in the blocks which tiononthe extreme dimers ofthe “rectangles”forming is often observed in this mode at some values of sam- the c(8×8) surface structure. ple bias and tunnelling current. Upper atoms of tilted The latter statement is illustrated by the STM 3-D dimers are shown by larger open circles. empty-statetopographshowninFig.11c.The extreme dimers located on the sides of the rectangular blocks are seen to be somewhat higher than the other ones 4.2 Comparison of STM and RHEED data of the dimer pairs; they form a superfine relief which turned out to be sufficient to backscatter fast electrons Now d discrepancy of results obtained by STM and incident on the surface at grazing angles. RHEED within the proposedmodel. Fig. 10 presents a sketchofthereciprocallatticeofc(8×8).TheRHEED patterns obtained in the [110] azimuth correspond to 4.3 Origin thec(8×8)structure;thepatternsobservedinthe[010] azimuthdonot(Fig.4).Thereasonofthisdiscrepancy The Si(001)-c(8×8) structure have formerly been ob- may be understood from the STM filled state image served and described in a number of publications [4,5, which corresponds to the electron density distribution 6,7,10]. Conditions of its formation were different: we of electrons paired in covalent bond of a Si–Si dimer. shallexplaintheobservecoperatomsweredepositedon Fig. 11 compares STM images of the same region on silicon (2×1) surface to form the c(8×8) reconstruc- one terrace obtained in the empty-state (a) and filled- tion [10], although it is known that Cu atoms are not state (b) modes; inserts show their Fourier transforms, absorbedontheSi(001)cleansurfaceifthesampletem- ◦ the differences in which for the two STM modes are peratureisgreaterthan600 C,andonthecontraryCu 9 (a) (b) (c) Fig. 11 STMimagesofthesameareaonthesurfaceobtained intheemptystate(a)andfrilledstate (b)modes (+1.96V,120pA and −1.96 V, 100 pA); for the convenience of comparison, ‘D’ indicates the same vacancy defect; corresponding Fourier transforms areshownintheinserts.A3-DSTMemptystatemicrograph(+2.0V,200pA)oftheSi(001)-c(8×8)surfaceisshownin(c). desorptionfromthesurfacetakesplace [7,10];fastcool- the defect surface reconstruction. However, the follow- ing from the annealing temperature of ∼ 1100◦C was ing arguments urge us to doubt about the Cu-based applied [4,5]; samples treated in advance by ion bom- model: (i) undetectable trace amounts of Cu were sug- bardmentwereannealedandrapidlycooled[6].There- gested in Ref. [5], the presence or absence of which is sultantsurfaces weremainly exploredby STMand low unprovable;(ii) evenifthe suggestionistrue,ourSTM energyelectrondiffraction(LEED).STMinvestigations imagesgiveanevidenceofadifferentamountofdimers yieldedalikeresults—abasicunit ofthe reconstruction intherectangularblocks,so,itisunclearwhyCuatoms wasa“rectangle”,butthestructureofthe“rectangles” formdifferent stable configurationsonsimilar surfaces; revealed by different authors was different. In general, and (iii) it is hard to explain why Cu atoms cyclically an origin of the Si(001)-c(8× 8) structure is unclear compose and decompose the rectangular blocks during now. the cyclical thermal treatments of the samples. It ap- plies equally to any other impurity or contamination. STMimagesmostresemblingourdatawerereported inRef.[5].Inthatpaper,thec(8×8)structurewasob- Now we consider a different interpretation of our served in samples without special treatment by coper: data.Asmentionedabove,literaturesuggeststwocauses thesamplesweresubjectedtoannealingatthetemper- of c(8×8) appearance.The first is an impact of impu- atureof∼1050◦Cfortheoxidefilmremoval.Formation rity atoms adsorbed on the surface even at trace con- of the c(8×8) reconstructionwas explained in that ar- centrations. The second is a thermal cycle of the oxide ticle by the presence of a trace amount of Cu atoms filmdecompositionandsamplecooling.Thefirstmodel the concentrationof which was beyond the Auger elec- seems to be hardly applicable for explanation of the tronspectroscopydetectionthreshold.TheSTMempty reported experimental results. According to our data, state images of the samples were similar to those pre- therearenoimpuritiesadsorbeddirectlyonthestudied sented in the current paper. A very important differ- surface:RHEEDpatternscorrespondtoacleanSi(001) ence is observed in the filled state images—we observe surface reconstructed in (2×1) or, at lower tempera- absolutely different configuration of dimers within the tures,(4×4)configuration.Cycliccontaminantdesorp- “rectangles”. Nevertheless, the presence of Cu cannot tionathightemperatures(?600◦C)andadsorptionon be completelyexcluded.Someamountofthe Cuatoms samplecoolingisunbelievable.Consecutivesegregation may come on the surface from the construction mate- and desegregation of an undetectable impurity in sub- rials of the MBE chamber (although there is a circum- surface layers also does not seem verisimilar. stance that to some extent contradicts this viewpoint: The second explanation looks more attractive. It Cuatomswerenotdetectedintheresidualatmosphere was found in Ref. [46] as a result of the STM studies of the MBE chamber within the sensitivity limit of the thattheSi(001)surfacesubjectedtothethermaltreat- SRS RGA-200 mass spectrometer) or even from the Si ment at ∼820◦C which was used for decomposition of wafer. Cu is known to be a poorly controllable impu- the thin (∼1 nm thick) SiO2 films obtained by chemi- rity and its concentrationin the subsurface layersof Si cal oxidation contained a high density of vacancy-type wafers which were not subjected to the gettering pro- defectsandtheiragglomeratesaswellasindividualad- cess may reach 1015 cm−3. This amount of Cu may dimers. So, the initial bricks for the considered surface appear to be sufficient to give rise to the formation of structure are abundant after the SiO2 decay. 10 (a) (b) (c) Fig. 12 Schematic representation ofthe surfacestress fields interactions duringformationofthe c(8×8)structure: (a) orderingof the“rectangles” withintherows;(b)orderingoftherowsrelativetoeachother;(c)theorderedc(8×8)structure. Literaturepresentsawideexperimentalmaterialon sists of abundant ad-atoms. On cooling the ad-dimers a different reconstructionof the Si(001)surface—c(4× have to migrate along the surface and be build in the 4)—which also arise at the temperatures of ? 600◦C. lattice. A number of competing sinks may exist on the Forexample,areviewofarticlesdescribingdifferentex- surface (steps, vacancies, etc.), but high cooling rate perimental investigations can be found in Refs. [12,31, may impede ad-atom annihilation slowing their migra- 32,33,34,35]. Based on the generalized data, an infer- tion to sinks and in such way creating supersaturation encecanbemadethatthec(4×4)structureformsinthe and favoring 2-D islanding, and freezing a high-order ◦ intervalfrom600to700 C.Mostlikely,atthesetemper- reconstruction. atures an appreciable migration of Si ad-atoms starts The followingscenariomaybe proposedtodescribe on surface. The structure is free of impurities. It irre- the c(8 × 8) structure formation. A large number of versibly transits to the (2×1) one at the temperature ad-dimers remains on the surface during the sample ◦ greater than 720 C. Ref. [47] demonstrates formation annealing after the oxide film removal. They form the of the Si(001)-(2×8) structure, also without impurity uppermost layer of the structure. The underlying layer atoms. In analogy with the above literature data, for- is (2×1) reconstructed.Ad-dimers are mobile and can mation of the c(8×8) reconstructionmay be expected form different complexes (islands). Calculations show asaresultoflow-temperatureannealingand/orfurther that the most energetically favorable island configura- quenching.Thestandardannealingtemperatureforob- tions are single dimer on a row in non-epitaxial orien- taining (2×1) structure is known to be in the interval tation[43,45,48,49](Fig.5b),complexesoftwodimers ◦ from 1200 to 1250 C. At these temperatures in UHV (pairs of dimers) in epi-orientation (metastable [43]) ambient not only oxide film removal from the surface and two dimers on a row in non-epitaxial orientation takesplace,butalsosiliconevaporationandcarbondes- separatedbya divacancy,andtripple-dimerepi-islands orptiongoeson.Unfortunately,wehavenotgotatech- consideredascriticalepinuclei[43].Thesemobiledimers nicalopportunitytocarryoutsuchahigh-temperature and complexes migrate in the stress field of the (2×1) ◦ annealing in our instrument. Treatment at 925 C that structure. The sinks for ad-dimers are (A) steps, (B) we apply likely does not result in substantial evapora- vacancydefects ofthe underlying (2×1) reconstructed tion of Si atoms from the surface, and C atoms, if any, layer, and (C) “fastening” them to the (2×1) surface may diffuse into subsurface layers. As a result, a great asa c(8×8)structure.The mainsinks athightemper- amount of ad-dimers arise on the surface, like it hap- aturesareAandB.Asthesampleiscooled,theCsink pens in the process described in Ref. [46]. Formation becomes dominating. Ad-dimers on the Si(001)-(2×1) processes of the (2×1) and c(8×8) structures are dif- surface are known to tend to form dimer rows [50]. In ferent.(2×1)ariseduringthehigh-temperatureanneal- this case such rows are formed by metastable dimer ingandad-atomsoftheuppermostlayerdonotneedto pairsgatheredinthe“rectangles”.The“rectangles”are migrate and be embedded into the lattice to form this orderedwithaperiodof8translationsintherows.The reconstruction.Onthecontrary,c(8×8)appearsduring ordering is likely controlled by the (2×1) structure of sample cooling, at rather low temperatures, and at the the underlying layer and interaction of the stress fields moment of a prior annealing the uppermost layer con- arising around each “rectangle”. Effect of the under-