We report on nano-patterning with multiply charged argon ions that facilitates the self-assembly of epitaxial Cu3Si nanostructures. In particular, we show that the impact sites formed from the dissipation of the incident ion potential energy for Arq+ (q=1,4,8) modulate the growth density and growth rate for silicide nanostructures. The observed nanostructure densities were found to vary as q0.9 for strain-driven, shape transition-type growth, and the observed growth rates far exceeded those obtained under thermal conditions. Relating the growth density to an underlying sputter yield for SiO2, we find a dependence on the ion potential energy relatively similar to that observed by others for Iq+ ions incident on a thermally grown oxide.

Metal silicide films are an integral part of modern microelectronics and find use in multiple contexts, such as Schottky diodes and local interconnects in metal-oxide–semiconductor devices.1,2 At the nanometer scale, they have been synthesized using a variety of techniques, including silicidation, chemical vapor transport, and chemical vapor deposition. Metal silicides of nanometer dimensions naturally form geometric structures that have distinct shapes and sizes governed by both their nucleation conditions and by strain evolution during growth. One example of such structures is Co on Si which initially forms square shapes that transition to clusters with large aspect ratios. While this and similar systems grown in the low deposition limit exhibit transitional growth behaviors consistent with the model of Tersoff and Tromp, their nucleation (both location and density) and growth rate are dependent on intrinsic defect sites in the material system.3,4 Attempts to control the growth have included the preferential placement of silicides using methods such as laser-induced periodic surface structures.5 

The Cu3Si system, in particular, has motivated recent interest based on its possible use in energy storage applications, such as batteries and solar cells.6,7 For example, Xu et al. investigated the use of Cu3Si in battery anodes and showed that an excellent electrochemical performance was obtained for a Si–Cu3Si–Cu composite electrode structure.8 Also, it has been shown recently that Cu3Si nanowires exhibit superior field-emission properties and can serve as highly efficient anti-reflective layers and as an effective catalyst in reactions.9 

The motivation for the work we present here is an investigation of the role that void formation from multicharged ion impacts has on the formation and growth of Cu3Si nanostructures. Specifically, it is known that annealing is required to successfully form Cu3Si nanostructures on SiO2/Si.10–12 As we showed in our recent work, SEM images taken under various annealing conditions reveal a linear, temperature-dependent variation in the density of voids or bare Si regions at the SiO2 surface.4 The observed dependence is linked to the importance of voids in Cu3Si nanostructure formation, particulary within an optimal temperature range of 450–600 °C. Multicharged ions are known to create distinct, fluence-dependent features at insulating surfaces, such as Al2O3 and SiO2.13–18 Therefore, it is reasonable to consider whether a controlled fluence of multicharged ions could supplement or even replace the annealing steps required for void formation in Cu3Si nanostructure growth.

In Sec. II, we present the methods and materials relevant to this work. The data obtained for multicharged ions inserted into the processing steps for Cu3Si nanostructure formation are shown in Sec. III. These results are discussed in the context of the energy and charge dissipation necessary to neutralize an incident multicharged ion in Sec. IV followed by our conclusions in Sec. V.

Single crystal, 3-inch diameter, boron doped p-type silicon wafers of (111) and (100) orientations with resistivities of 1–10 Ω cm were employed as substrates for this work. Prior to the oxide growth, the wafers were cleaned following standard RCA clean, which involves cleaning of the substrates in a chemical bath (1:1:5 solution of NH4OH + H2O2+H2O) for five minutes under ultrasonic agitation. This is followed by etching of cleaned surface with 1% dilute HF for two minutes and triple rinse in de-ionized water for six minutes to remove any native oxide. Following this cleaning, an ultra-thin thermal oxide layer of approximately 1.5±0.05 nm was grown via dry oxidation of Si at 200 °C for 2 h in a hot furnace. Oxide thickness measurements were performed at five different locations on each wafer (center and four corners) using a Sopra model GES-5E spectroscopic ellipsometer. The oxide grown Si wafers were diced into 10 × 10 mm2 squares and mounted on Omicron style sample platen and loaded into the target chamber.

The samples were exposed to the focused multicharged ion beams produced in the CUEBIT laboratory which is discussed in detail in Ref. 19. The pressure inside the target chamber during the irradiations was on the order of 108 mbar and the ion fluence was in the range of 10111012 per cm2. The ion species used for the irradiations were Ar1+, Ar4+, and Ar8+ which have potential or neutralization energies of 15.8, 143.9, and 577.7 eV, respectively. The kinetic energy of all the ion beams was 1 keV. Beam profiles were measured using a (HRBIS-4000 from Beam Imaging Solutions) as well as a translatable Faraday cup. An example beam profile, along with relevant information regarding profiles and fluence calculations, is shown in Ref. 20.

After irradiation, the samples were unloaded from the target chamber and transferred to a deposition chamber. Thin Cu films were deposited on these irradiated samples at temperature of 600 °C using a McAllister eBeam Evaporator at a pressure of 5×106 Torr. The deposition was kept to a 0.1 mA flux current and a deposition time of 5 min. The deposition rate was estimated to be 1 monolayer/min (ML/min). The sequence of steps involved in the growth process of the resulting Cu3Si nanostructures is shown in Fig. 1. Following deposition, surface morphologies of the sample were investigated using a Hitachi S4800 scanning electron microscope. Orientation and chemical morphology of the nano-islands were established using energy dispersive X-ray (EDX) analysis. Figures 2(a) and 2(b) and Table I show the size distribution (length, width, area, density) of the nanostructures formed for SiO2/p-Si(100).

FIG. 1.

Growth process of MCI-assisted Cu3Si nanostructures and SEM image showing irradiated and non-irradiated growth regions.

FIG. 1.

Growth process of MCI-assisted Cu3Si nanostructures and SEM image showing irradiated and non-irradiated growth regions.

Close modal
FIG. 2.

(a) SEM images showing Cu3Si islands grown on SiO2/p-Si(100) as a function of fluence and charge state with preferred orientation as measured by EDX. (b) Dependence of length (closed) and width (open) circles for SiO2/p-Si(100) showing a critical size consistent with the model of Tersoff and Trump.3 

FIG. 2.

(a) SEM images showing Cu3Si islands grown on SiO2/p-Si(100) as a function of fluence and charge state with preferred orientation as measured by EDX. (b) Dependence of length (closed) and width (open) circles for SiO2/p-Si(100) showing a critical size consistent with the model of Tersoff and Trump.3 

Close modal
TABLE I.

Size distributions of nano-islands on SiO2/p-Si(100) and SiO2/p-Si(111).

Growth modeLength (min–max)Width (min–max)Density (min–max)Area (min–max)Deposition timeVoid nucleation
method: substrate(nm)(nm)(×106 cm−2)(×103 nm2)(min)method
RDE: p(100) 286–1402 88–281 2.12–25.5 44–123.7 1–5 Thermal annealing 
HCI: p(100) 388–2614 137–493 1.1–13.7 170.3–762.0 HCI irradiation 
RDE: p(111) 282–2877 160–400 0.87–7.3 44.9–980 1–5 Thermal annealing 
HCI: p(111) 444–4281 291–818 1.2–8.9 120.6–758.2 HCI irradiation 
Growth modeLength (min–max)Width (min–max)Density (min–max)Area (min–max)Deposition timeVoid nucleation
method: substrate(nm)(nm)(×106 cm−2)(×103 nm2)(min)method
RDE: p(100) 286–1402 88–281 2.12–25.5 44–123.7 1–5 Thermal annealing 
HCI: p(100) 388–2614 137–493 1.1–13.7 170.3–762.0 HCI irradiation 
RDE: p(111) 282–2877 160–400 0.87–7.3 44.9–980 1–5 Thermal annealing 
HCI: p(111) 444–4281 291–818 1.2–8.9 120.6–758.2 HCI irradiation 

Here, we demonstrate the use of MCIs as a processing tool in the growth of epitaxial copper silicide on p-doped Si substrates. Multiple samples were exposed to beams of 1 keV Arq+ ions (q=1,4,8) at the CUEBIT facility19 as shown schematically in Fig. 1. Subsequent copper deposition at 600°C for 5 min at a deposition rate of 1 ML/min led to the formation of a high density of nano-islands of Cu3Si. As is evident from the SEM image in Fig. 1 the high density growth corresponds to the region irradiated by the MCIs. Traditionally, efficient epitaxial nanostructure growth on Si/SiO2 substrates requires the formation of voids through annealing of the samples at high temperatures (600–1000 °C) for 10–12 hours.4,21,22 In our case, this step was replaced by irradiation with MCIs and nanostructures were still formed.

Representative data obtained for a Si(100) substrate are shown in Figs. 2 and 3. From Fig. 2, one can see a variation in the shapes of the Cu3Si nanostructures that depends on the Arq+(1q8) charge state used and on the ion fluence. For low fluence values (12.5×1011cm2) and q4, the observed structures are squares reflecting the fourfold symmetry of the substrate. At a higher fluence (2.5×1012cm2) and q4, a change in the observed shape from square to nanowire is seen. No square-shaped nanostructures were observed for Ar8+ ions even at the lowest fluence, and the observed shapes were either nanowire or a combination of nanowire and rectangular islands. In the latter case, the rectangular islands are observed to be present on the top of nanowires and positioned approximately at the center. The orientations of the nanowires were found to be along [1¯10] and [110] directions without any preference.

FIG. 3.

Mean density of Cu3si nanostructures as a function of multicharged ion fluence on SiO2/p-Si(100) substrates.

FIG. 3.

Mean density of Cu3si nanostructures as a function of multicharged ion fluence on SiO2/p-Si(100) substrates.

Close modal

Figure 3 shows the density of the Cu3Si nano-islands as a function of fluence with a representative SEM image as the inset. The mean density values were obtained by processing multiple SEM images that spanned the irradiated region of each sample. The density scales linearly indicating that each additional ion-surface interaction (increased fluence) is leading to the formation of additional Cu3Si nanostructures (increased density) on the SiO2/p-Si(100) substrate. However, this is also q-dependent, with the higher charge state (+8) leading to significantly more nanostructure formation. Similar results were seen on the Si(111) substrate. The fluence-dependent area of the Cu3Si nanostructures is shown in Fig. 4 for Ar4+ and Ar8+ irradiations and deposition times of 5 min. These data are compared with the time-dependent area of the Cu3Si nanostructures formed through the deposition and thermal annealing process (no irradiation) over a 1–5 min deposition time range. It is evident that the areas of the ion-assisted Cu3Si nanostructures are much larger compared to those formed via the thermally annealed process.

FIG. 4.

Area of Cu3Si nanostructures as a function of fluence for Ar4+ and Ar8+. For comparison, areas obtained from Cu3Si structures formed under a thermal annealing process (no ion irradiation) are also shown.

FIG. 4.

Area of Cu3Si nanostructures as a function of fluence for Ar4+ and Ar8+. For comparison, areas obtained from Cu3Si structures formed under a thermal annealing process (no ion irradiation) are also shown.

Close modal

The results we have obtained for Arq+ irradiation of p-type Si indicate that the ions and their charge state can play a significant role in the formation of Cu3Si nanostructures. As we noted in Sec. I, voids in the SiO2 layer are required to promote the formation of the silicide islands. In the absence of a thermal annealing step, the ion beam irradiation appears to be a generator of these voids as evidenced by the clear delineation between irradiated and non-irradiated growth regions in the SEM image of Fig. 1. Moreover, the island density can be tied quantitatively to the specific ion irradiation fluence as shown in Fig. 3. For both Ar4+ and Ar8+ ions the linear relationship between fluence and feature density indicates that each impact generates a possible island nucleation site with some q-dependent probability. This is consistent with fluence-dependent features observed in multicharged ion measurements on other substrates such as Al2O3.13 

The slope of the lines from Fig. 3 can be used to calculate a feature density per ion which is plotted in Fig. 5 as a function of the ion charge state. Fitting this dependence to a power law for our p-Si(100) substrate, i.e., density qγ, we find γ=0.95±0.04 or a nearly linear dependence on charge state. For the p-Si(111) substrate we find a slightly lower value of γ=0.90±0.04. According to the Coulomb explosion model proposed by Tona et al.,23 the relation between the sputter yield of positive ions (Yp) and potential energy of a multicharged ion (Ep) was determined to be YpEpΓ, where Γ was found to be 0.6 for the sputter yields of O+ ions due to Iq+ ion impacts on native SiO2/Si(111) targets. For our Arq+ ion charge states, q=1, 4 and 8 correspond to reneutralization or potential energies of 15.8, 143.9, and 577.7 eV, respectively. Assuming that our observed Cu3Si nanostructures densities are formed due to voids created by the multicharged ion irradiation at the surface, which in turn depend on the number of bonds broken (or equivalently, the sputter yield), we can assume that the densities of the nanostructures are directly related to the sputter yield for the SiO2 substrate. Hence treating the densities to be proportional to the underlying sputter yield, and following the same relation of YpEpΓ as in Tona et al., we find Γ=0.58±0.18 and 0.53±0.06 for Arq+ on SiO2/Si(100) and SiO2/Si(111), respectively. This is remarkably similar to the Tona et al. result,23 and the observed differences could be attributed to the differing oxides used in the two experiments (native vs thermally grown). Although we can see that our data are consistent with the Coulomb explosion model, it is worth noting that charge state enhanced preferential sputtering24 and defect mediated sputtering25 have been observed to play a role in other studies.

FIG. 5.

Normalized density per ion of Cu3Si nanostructures formed on p-Si(100) and p-Si(111) substrates irradiated by Arq+ ions.

FIG. 5.

Normalized density per ion of Cu3Si nanostructures formed on p-Si(100) and p-Si(111) substrates irradiated by Arq+ ions.

Close modal

Regarding the morphology of the nanostructures formed, we observe that there is a transition from squares to wires for q=4 and 8. Such a transition is consistent with the Tersoff and Trump model3—an energetics based strain driven shape transition model. According to this model, the nanostructure formation minimizes the strain energy of an island per unit volume. The model predicts a critical size below which a compact shape (L=W) is observed, but when exceeded, rectangular nanoshapes (L>W) are observed. Figure 2(b) shows the critical size at which the shape transitions from a square to a wire for SiO2/p-Si(100).

The formation of Cu3Si nanostructures on p-type Si depends on void formation in the silicon substrate. These voids are traditionally created by thermal annealing at high temperatures (450600°C) for a period of 12 h. In this work, we have shown that the annealing step can be replaced by controlled multicharged ion irradiation, as demonstrated by the comparison of the density of Cu3Si nanoislands for irradiated and non-irradiated samples. Specifically, silicon oxide samples were irradiated with multicharged ArQ+ ion beams with charge state Q (Q = 1,4,8) and corresponding neutralization (potential) energies Ep (Ep=15.8, 143.9, and 577.7 eV). The irradiated samples were subsequently exposed to a Cu thin film deposition process. The surface and chemical morphology of the resulting Cu3Si nanostructures were measured using SEM and energy dispersive x-ray analysis, respectively, and a significantly higher density of nanoislands was seen in the samples that were irradiated with multicharged ions. Further, the calculated feature density per ion showed a near linear dependence on charge state (slope = 0.90±0.04). Assuming that nanostructure formation depends on void formation, and void formation in turn directly depends on sputtering, sputter yields may be considered to be directly dependent on the densities of nanostructure formation. Based on a Coulomb explosion model, a power law relating the sputter yields to the neutralization energy of multicharged ions was proposed by Tona et al.23 By applying the proposed power law to the density of nanostructure formation, we obtain a sputter yield of 0.58±0.18 and 0.53±0.06 due to ArQ+ irradiation for oxides on Si(100) and Si(111) substrates, respectively, which is in close agreement to the sputter yield of 0.6 reported for O+ ions due to IQ+ irradiation on oxides grown on Si(111) substrates.

This research was funded by the Defense Advanced Research Projects Agency (Army Research Office Agreement No. W911NF-13-1-0042).

The authors have no conflicts to disclose.

E. S. Srinadhu: Conceptualization (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Writing – original draft (equal); Writing – review & editing (equal). D. D. Kulkarni: Investigation (equal); Writing – review & editing (equal). D. A. Field: Investigation (equal). J. E. Harriss: Data curation (equal); Investigation (equal). C. E. Sosolik: Conceptualization (equal); Funding acquisition (equal); Investigation (equal); Project administration (equal); Supervision (equal); Writing – original draft (equal); Writing – review & editing (equal).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
F.
Deng
,
R. A.
Johnson
,
P. M.
Asbeck
,
S. S.
Lau
,
W. B.
Dubbelday
,
T.
Hsiao
, and
J.
Woo
, “
Salicidation process using NiSi and its device application
,”
J. Appl. Phys.
81
,
8047
(
1997
).
2.
H.
Iwai
,
T.
Ohguro
, and
S.
Ohmi
, “
NiSi salicide technology for scaled CMOS
,”
Microelectron. Eng.
60
,
157
(
2002
).
3.
J.
Tersoff
and
R. M.
Tromp
, “
Shape transition in growth of strained islands: Spontaneous formation of quantum wires
,”
Phys. Rev. Lett.
70
,
2782
(
1993
).
4.
E. S.
Srinadhu
,
J. E.
Harriss
, and
C. E.
Sosolik
, “
Shape transitions of Cu3Si islands grown on Si(111) and Si(100)
,”
Appl. Surf. Sci.
465
,
201
(
2019
).
5.
P.
Nurnberger
,
H.
Reinhardt
,
D.
Rhinow
,
R.
Riedel
,
S.
Werner
, and
N.
Hampp
, “
Controlled growth of periodically aligned copper-silicide nanocrystal arrays on silicon directed by laser-induced periodic surface structures (LIPSS)
,”
Appl. Surf. Sci.
420
,
70
(
2017
).
6.
Y. T.
Bie
,
J. L.
Yu
,
J.
Yang
,
W.
Lu
,
Y. N.
Nuli
, and
J. L.
Wang
, “
Porous microspherical silicon composite anode material for lithium ion battery
,”
Electrochim. Acta
178
,
65
(
2015
).
7.
J. B.
Zhou
,
N.
Lin
,
Y.
Han
,
J.
Zhou
,
Y. C.
Zhu
,
J.
Du
, and
Y. T.
Qian
, “
Cu3Si@Si core-shell nanoparticles synthesized using a solid-state reaction and their performance as anode materials for lithium ion batteries
,”
Nanoscale
7
,
15075
(
2015
).
8.
K.
Xu
,
Y.
He
,
L.
Ben
,
H.
Li
, and
X.
Huang
, “
Enhanced electrochemical performance of Si–Cu–Ti thin films by surface covered with Cu3Si nanowires
,”
J. Power Sources
281
,
455
(
2015
).
9.
F. W.
Yuan
,
C. Y.
Wang
,
G. A.
Li
,
S. H.
Chang
,
L. W.
Chu
,
L. J.
Chen
, and
H. Y.
Tuan
, “
Solution-phase synthesis of single-crystal Cu3Si nanowire arrays on diverse substrates with dual functions as high-performance field emitters and efficient anti-reflective layers
,”
Nanoscale
5
,
9875
(
2013
).
10.
Y.
Cechal
,
J.
Polcak
,
M.
Kolibal
,
P.
Babor
, and
T.
Sikola
, “
Formation of copper islands on a native SiO2 surface at elevated temperatures
,”
Appl. Surf. Sci.
256
,
3636
(
2010
).
11.
N.
Benouattas
,
A.
Mosser
,
D.
Raiser
,
J.
Faerber
, and
A.
Bouabellou
, “
Behaviour of copper atoms in annealed Cu/SiOx/Si systems
,”
Appl. Surf. Sci.
153
,
79
(
2000
).
12.
S.
Li
,
H.
Cai
,
C. L.
Gan
,
J.
Guo
,
Z.
Dong
, and
J.
Ma
, “
Controlled synthesis of copper-silicide nanostructures
,”
Cryst. Growth Des.
10
,
2983
(
2010
).
13.
R.
Lake
,
J.
Pomeroy
, and
C.
Sosolik
, “
Energy dissipation of highly charged ions on Al oxide films
,”
J. Phys.: Condens. Matter
22
,
084008
(
2010
).
14.
J. M.
Pomeroy
,
R. E.
Lake
, and
C. E.
Sosolik
, “
Highly charged ion interactions with thin insulating films
,”
Nucl. Instrum. Methods Phys. Res. Sect. B
269
,
1238
(
2011
).
15.
R. E.
Lake
,
J. M.
Pomeroy
,
H.
Grube
, and
C. E.
Sosolik
, “
Charge state dependent energy deposition by ion impact
,”
Phys. Rev. Lett.
107
,
063202
(
2011
).
16.
P.
Karmakar
,
S.
Bhattacharjee
,
V.
Naik
,
A.
Sinha
, and
A.
Chakrabarti
, “
Coulomb explosion sputtering of selectively oxidized Si
,”
J. Phys.: Condens. Matter
22
,
175005
(
2010
).
17.
W.
Arnoldbik
,
N.
Tomozeiu
, and
F.
Habraken
, “
Electronic sputtering of thin SiO2 films by meV heavy ions
,”
Nucl. Instrum. Methods Phys. Res. Sect. B
203
,
151
(
2003
).
18.
M.
Sporn
,
G.
Libiseller
,
T.
Neidhart
,
M.
Schmid
,
F.
Aumayr
,
H.
Winter
, and
P.
Varga
, “
Potential sputtering of clean SiO2 by slow highly charged ions
,”
Phys. Rev. Lett.
79
,
945
(
1997
).
19.
R.
Shyam
,
D. D.
Kulkarni
,
D. A.
Field
,
E. S.
Srinadhu
,
D. B.
Cutshall
,
W. R.
Harrell
,
J. E.
Harriss
, and
C. E.
Sosolik
, “
First multicharged ion irradiation results from the CUEBIT facility at Clemson University
,”
AIP Conf. Proc.
1640
,
129
(
2015
).
20.
D. D.
Kulkarni
,
R. E.
Shyam
,
D. B.
Cutshall
,
D. A.
Field
,
J. E.
Harriss
,
W. R.
Harrell
, and
C. E.
Sosolik
, “
Tracking subsurface ion radiation damage with metal–oxide–semiconductor device encapsulation
,”
J. Mater. Res.
30
,
1413
1421
(
2015
).
21.
L. J.
Chen
,
J. J.
Chu
,
W.
Lur
,
H. F.
Hsu
, and
T. C.
Lee
, “Recent developments in the epitaxial growth of transition metal silicides on silicon,” in Heterostructures on Silicon: One Step Further with Silicon, NATO ASI Series (Series E: Applied Sciences), edited by Y. Nissim and E. Rosencher (Springer, Dordrecht, 1989), Vol. 160, pp. 231–238.
22.
H.
Hibino
,
M.
Uematsu
, and
Y.
Watanabe
, “
Void growth during thermal decomposition of silicon oxide layers studied by low-energy electron microscopy
,”
J. Appl. Phys.
100
,
113519
(
2006
).
23.
M.
Tona
,
S.
Takahashi
,
K.
Nagata
,
N.
Yoshiyasu
,
C.
Yamada
,
N.
Nakamura
,
S.
Ohtani
, and
M.
Sakurai
, “
Coulomb explosion potential sputtering induced by slow highly charged ion impact
,”
Appl. Phys. Lett.
87
,
224102
(
2005
).
24.
T.
Mizutani
, “
Preferential sputtering of oxygen from SiO2 by low-energy ion beam and neutral beam bombardment
,”
Jpn. J. Appl. Phys.
30
,
L628
(
1991
).
25.
M.
Sporn
,
G.
Libiseller
,
T.
Neidhart
,
M.
Schmid
,
F.
Aumayr
,
H.
Winter
,
P.
Varga
,
M.
Grether
,
D.
Niemann
, and
N.
Stolterfoht
, “
Potential sputtering of clean SiO2 by slow highly charged ions
,”
Phys. Rev. Lett.
79
,
945
(
1997
).