Next Article in Journal
Exposure to Ultrafine Particles in the Ferroalloy Industry Using a Logbook Method
Previous Article in Journal
Broadband Visible Nonlinear Absorption and Ultrafast Dynamics of the Ti3C2 Nanosheet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Kinetics and Stoichiometry of Metal Cation Reduction on Multi-Crystalline Silicon in a Dilute Hydrofluoric Acid Matrix

Department of Physical Chemistry, Faculty of Environment and Natural Sciences, Brandenburg University of Technology Cottbus-Senftenberg, Universitätsplatz 1, 01968 Senftenberg, Germany
*
Author to whom correspondence should be addressed.
Nanomaterials 2020, 10(12), 2545; https://doi.org/10.3390/nano10122545
Submission received: 9 November 2020 / Revised: 11 December 2020 / Accepted: 13 December 2020 / Published: 17 December 2020

Abstract

:
In this study, the process of metal cation reduction on multi-crystalline silicon in a dilute hydrofluoric acid (HF) matrix is described using Ag(I), Cu(II), Au(III) and Pt(IV). The experimental basis utilized batch tests with various solutions of different metal cation and HF concentrations and multi-crystalline silicon wafers. The metal deposition kinetics and the stoichiometry of metal deposition and silicon dissolution were calculated by means of consecutive sampling and analysis of the solutions. Several reaction mechanisms and reaction steps of the process were discussed by overlaying the results with theoretical considerations. It was deduced that the metal deposition was fastest if the holes formed during metal ion reduction could be transferred to the valence bands of the bulk and surface silicon with hydrogen termination. By contrast, the kinetics were lowest when the redox levels of the metal ion/metal half-cells were weak and the equilibrium potential of the H3O+/H2 half-cells was high. Further minima were identified at the thresholds where H3O+ reduction was inhibited, the valence transfer via valence band mechanism was limited by a Schottky barrier and the dissolution of oxidized silicon was restricted by the activity of the HF species F, HF2 and H2F3. The findings of the stoichiometric conditions provided further indications of the involvement of H3O+ and H2O as oxidizing agents in addition to metal ions, and the hydrogen of the surface silicon termination as a reducing agent in addition to the silicon. The H3O+ reduction is the predominant process in dilute metal ion solutions unless it is disabled due to the metal-dependent equilibrium potential of the H3O+/H2 half-cell and the energetic level of the valence bands of the silicon. As silicon is not oxidized up to the oxidation state +IV by the reduction of the metal ions and H3O+, water is suspected of acting as a secondary oxidant. The stoichiometric ratios increased up to a maximum with higher molalities of the metal ions, in the manner of a sigmoidal function. If, owing to the redox level of the metal half-cells and the energetic level of the valence band at the metal–silicon contact, the surface silicon can be oxidized, the hydrogen of the termination is the further reducing agent.

1. Introduction

The wet chemical structuring of silicon is of great importance for various applications in photonics, microsystems technology, microelectronics and sensor technology [1]. The process of metal-assisted etching to create special spatial structures is the topic of interest [2,3,4]. Despite very intensive investigations, not all details of the reactions occurring during metal-assisted etching have yet been clarified [5]. An essential partial aspect of the research question refers to the process of electroless precious metal cation reduction on silicon in a diluted hydrofluoric acid (HF) matrix.
At the beginning of this process, the metal nucleation starts mainly at crystal defects or kinks on the silicon surface [6,7,8]. Further metal growth is a highly dynamic process in which already formed metallic nuclei are oxidized again to metal cations in order to subsequently be deposited again as nuclei at energetically more favorable locations [9,10]. The diffusion of the metal cations to the silicon surface and further to the nuclei is the rate-limiting factor of metal deposition [11,12,13,14,15]. The nuclei expand to nanoclusters [16], islands or dendrites [17,18,19,20] on the silicon surface according to the Stranski–Krastanov and Volmer–Weber principles [21,22], and then form thin (nm to µm) [6,12,19,23,24], usually microporous layers [16,17] on the silicon surface. Adhesion forces between metal and silicon are usually low [6,25,26,27]. Depending on the metal, different shapes are formed on the silicon surface, such as pits (Cu) [18,28,29], pores (Ag, Au, Pt) [6,16,19,30], threads or nanowires (Ag, Au) [16,18,19,23,24,29,30,31,32,33,34]. The silicon removal occurs linearly in time, independent of the type of growth of the metals, the type of silicon doping, the amount of doping or the orientation of the silicon surface [10,19,35]. The metal deposition is also linear up to a certain thickness of the metal layer [36] and then follows first-order kinetics [12,36,37].
The reactions behind this process are reported in different ways in the literature. According to a prominent theory, there is a kinetically controlled transfer of two electronic holes hv+ from the metal cation/metal redox pair Mez+/Me to the valence band of the silicon (divalent mechanism, Equation (1)) [16,38,39,40], or in the case of copper deposition, an electron uptake of the metal ions via the conduction band mechanism of silicon [11]. Furthermore, a twofold electron transfer ec from the silicon conduction band to the hydrogen ions or other HF species contained in the solution results in the formation of Si4+ and molecular hydrogen (Equation (2)) [16,38,39,40].
2 z   Me z + 2 z   Me   + 2   h v +   and   Si   +   2   h v +   Si 2 +
Si 2 + Si 4 + +   2   e c   and   2   H + +   2   e c     H 2
According to another prominent thesis, there is a fourfold valence transfer between the metal cations and the silicon atoms via the silicon valence band (tetravalent mechanism, Equation (3)) without the formation of molecular hydrogen. It is characterized by a fast reaction with bulk silicon atoms [18,24,29,41] and rather slow kinetics with the silicon atoms of the hydrogen surface terminations [26,41] formed by the HF matrix.
4 z   Me z + 4 z   Me   + 4   h v +   and   Si   +   4   h v +   Si 4 +
It was shown in the study by daRosa et al. [42] that the stoichiometric ratios of Cu2+ reduction to metallic Cu and the associated oxidation and dissolution of silicon during the reaction of Cu2+ with silicon in HF matrix varied between 3:4 and 2:1 mol∙kg−1:mol∙kg−1. This was interpreted as an indication of a change between the divalent and tetravalent reaction mechanism depending on the concentration of Cu2+ and HF. In this context, Chartier et al. [43] formulated Equation (4) for the anodic reaction on silicon.
Si   +   6   HF   + n   h v + H 2 SiF 6 + n   H + + ( 4 n 2 ) H 2
The difference in the redox potentials [17,44,45], i.e., free enthalpy of the redox half-cells involved, is considered as the driving force of the reactions [11,17]. After the initial metal nucleation, the consecutive metal cation reduction occurs at the metal–silicon interface. The background is the Fermi energy alignment at the metal–silicon interface and a valence and conduction band bending of the silicon resulting from the difference in the work function for the metal and the silicon. Depending on the difference in work functions, this process causes the transfer of holes from the metal cation reduction at the level of Fermi energy through the higher work function element (usually the metal) to the valence band of the lower work function element (usually the silicon) [12,19,20,29,38,46]. Metals with particularly high work functions, such as platinum, can also develop a Schottky barrier behavior at the metal–silicon contact, which acts as a potential barrier and inhibits valence transfer kinetically [47]. The lower the Schottky barrier height [20,46,48] and the greater the metal valence band overlap with the silicon [29], the faster the holes resulting from metal reduction should be transferred to the silicon [20,23,24,33,46,48].
The roles of HF and its dissociation products are ambivalent in the process. Hydrogen ions hypothetically participate in silicon oxidation, as described in Equations (2) and (4) [16,38,39,40,43,49,50,51,52]. Undissociated HF (HFundiss.) [13,39,53,54] and/or the anionic dissociation products of the HF (F [33,36,55,56,57], HF2 [43,58,59]) are suspected of being responsible for the dissolution reaction of the oxidized silicon, resulting in the formation of the formal end product SiF62− [53], according to Equation (5).
Si 4 + + 6   F or   3   HF 2 or   6   HF undiss . SiF 6 2 ±   3..6   H +
The kinetics of the spatial growth of a metal on the silicon and the dissolution of the silicon should be faster with higher HF concentrations [13,50,51]. By contrast, the formation of a dense metal deposit on the silicon surface can restrict the access of HF species to the oxidized silicon, and thus, reduce the rate-limiting process of silicon dissolution [33].
The formation of hydrogen terminations on the silicon in the HF-matrix [60] has the consequence that only those metals can be deposited whose Mez+/Me half-cell redox level is higher than that of the hydrogen half-cell [11,56]. Furthermore, it is known that fluoride adsorption on the metal deposited can generate a strong dipole moment in the boundary layer, which increases the work function of the metal [61], thereby influencing the band bending mechanism at the metal–silicon contact, and in turn, the valence transfer between metal cation and silicon.
In conclusion, it can be stated that the process of electroless metal deposition is characterized by contradictory findings and theories. So far, a comparative study of the electroless deposition of different metals on silicon under identical conditions is missing.
The aim of this study was to determine and compare the kinetics of silver, copper, gold and platinum deposition on silicon experimentally, to determine the influences of different concentrations of the metal cations Ag+ and Cu2+; the anionic complexes AuCl4 and PtCl62−; and HF on the kinetics of the respective metal depositions and to draw conclusions regarding the enhancing or limiting factors of the processes.
The data obtained experimentally were further evaluated regarding the stoichiometric ratios of metal cation reduction to metal state, silicon oxidation and consecutive dissolution. Conclusions concerning general and metal-dependent reaction processes will be drawn in the summary of the results.

2. Materials and Methods

Common procedure: The experiments performed consisted of batch tests with solutions of different compositions of ultrapure water (18 MΩ resistance, PURELAB® flex, ELGA VEOLIA, Celle, Germany); HF (40%, w/w, p. A. quality, EMSURE®, Merck KGaA, Darmstadt, Germany); and Ag+, Cu2+, AuCl4 or PtCl62− in different concentrations. Multi-crystalline silicon wafers were added in all experiments in the same amount. The batches were stirred under argon atmosphere and sampled several times according to a fixed schedule. The dissolved Si and Ag, Cu, Au or Pt species were subsequently determined by inductively coupled plasma atomic emission spectroscopy (ICP OES).
Preparation of the solutions: A stock solution based on AgF was prepared for the silver batch tests. It was obtained by reacting AgNO3 (≥99%, ACS reagent, Sigma-Aldrich, St. Louis, MO, USA) with NaOH (25%, w/w, extra pure, Carl Roth, Karlsruhe, Germany) to form Ag2O and its subsequent conversion with 40% HF solution (w/w), resulting in a final concentration of 2.9 mol∙kg−1 AgF/2.9 mol∙kg−1 HF. A total of 27 solutions with a respective mass msol of 0.060 ± 0.001 kg were prepared from the dilution of this stock solution and the addition of 40% HF solution (w/w). Initial HF molalities of b HF (diss., t = 0) = 0.059 ± 0.001, 0.295 ± 0.004 and 1.472 ± 0.072 mol∙kg−1 were set for one-third each. The initial Ag+ molalities b Ag (diss., t = 0) were between 1.2 × 10−6 and 1.3 × 10−1 mol∙kg−1.
The Cu(II) solutions were based on the addition of CuSO4·5H2O (≥98%, ACS reagent, Sigma-Aldrich, St. Louis, MO, USA). In summary, 24 solutions were prepared with respective mass msol of 0.061 ± 0.001 kg. The initial Cu(II) molalities b Cu (diss., t = 0) were between 2.1 × 10−5 and 4.1 × 10−1 mol∙kg−1. The initial HF molalities were 0.060 ± 0.001 mol∙kg−1 in eight batches, 0.299 ± 0.004 mol∙kg−1 in nine batches and 1.499 ± 0.004 mol∙kg−1 in seven batches.
A commercially available HAuCl4 solution (40–44% Au, w/w, ≥98%, ACS reagent, Sigma-Aldrich, St. Louis, MO, USA) was used as the stock solution for the Au(III) batches. From this, nine solutions with initial AuCl4 molalities of b Au (diss., t = 0) = 3.3 × 10−7 to 9.3 × 10−2 mol∙kg−1 were prepared. The masses of the solutions (msol) were 0.060 ± 0.001 kg. The initial HF molalities of these solutions were 1.499 ± 0.006 mol∙kg−1. Furthermore, four identically concentrated HF solutions with AuCl4 molalities between 6.7 × 10−5 and 3.7 × 10−3 mol∙kg−1 were prepared with the addition of 2.2 × 10−4 mol∙kg−1 of para-Chlorobenzoic acid (≥98%, Alfa Aesar, Kandel, Germany).
The Pt(IV) etching solutions contained additions of H2PtCl6·6H2O (≥37.5% Pt base, ACS reagent, Sigma-Aldrich, St. Louis, MO, USA). Ten batches were prepared with initial masses msol of 0.060 ± 0.001 kg and initial PtCl62− molalities between b Pt (diss., t = 0) = 4.2 × 10−7 and 9.7 × 10−3 mol∙kg−1 at initial HF molalities of 1.493 ± 0.033 mol∙kg−1. Analogously to the Au(III) batches, four further HF solutions with initial PtCl62− molalities between 1.6 × 10−5 and 1.7 × 10−3 mol∙kg−1 with the addition of 2.2 × 10−4 mol∙kg−1 of para-Chlorobenzoic acid were prepared.
Silicon wafers: ten undoped multi-crystalline silicon wafers (15 mm edge length, 0.175 mm thickness, SCHOTT Solar, Mainz, Germany) pretreated with 1% HF solution were used for each experiment, each with a total mass of mSi = 0.92 ± 0.03 g and a total surface area of ASi = 46 cm2.
Experimental procedure: Each HF solution was transferred into a 150 mL low-density polyethylene vessel (VITLAB, Großostheim, Germany) and stirred using a polytetrafluoroethylene (PTFE)-coated magnetic stir bar (AlNiCo, 8 × 20 mm, VWR collection, Radnor, PA, USA) and a stirring plate (RT 10 power, IKA®, Staufen, Germany) at about 125 rpm at T = 295 ± 2 K. A GL-45 screw cap made of PTFE with three attached GL-14 ports was used to close the vessel (Bohlender, Grünsfeld, Germany). The HF pretreated silicon wafers were attached to one of the ports using a PTFE thread. The vessel was rinsed via the other two ports with an argon flow of 45 ± 1 mL via PTFE hose lines (using a Gas Flow Meter, Yokogawa, Ratingen Germany and Ar 6.0, Praxair, Danbury, CT, USA). The silicon wafers were then transferred into the etching solution by pulling the PTFE thread externally. The stirring of the etching solution and rinsing with argon was maintained continuously. The solutions were sampled seven times in a period of t = 3600 s (Ag, Cu, Au) or t = 14,400 s (Pt) via the free GL14 port for the time of t = 0, 300, 600, 1200, 1800, 2700 and 3600 s (Ag, Cu, Au) or t = 0, 1200, 2400, 4800, 7200, 10,800 and 14,400 s (Pt).
Sample analysis: All samples were diluted in 1% HF matrix (v/v, dilution of 40% HF solution). They were subsequently analyzed by ICP OES (iCAP Duo 6500, Thermo-Fisher, Waltham, MA, USA). The spectral intensities λ = 243.7, 328.0 and 338.2 nm for the Ag; λ = 224.7, 324.7 and 327.3 nm for Cu; λ = 201.2, 242.7 and 267.5 nm for Au; λ = 214.4, 265.9 and 306.4 nm for Pt; and λ = 212.4, 251.6 and 252.8 nm for Si were measured axially to the torch. The molalities of the elements (+/− confidence interval) b (mol∙kg−1) were determined using external calibrations based on dilution series of 1000 mg∙kg−1 standard stock solutions (Cu, Ag, Au, Pt and Si, Certipur®, Merck KGaA, Darmstadt, Germany).

3. Results and Discussion

3.1. Kinetics of the Metal Depositions at One Initial Metal Ion and HF Concentration Levels

Figure 1a–d shows the results of the Ag, Cu, Au and Pt deposition experiments with initial metal ion molalities (b Me (diss., t = 0) with Me = Ag, Cu, Au or Pt) of ≈ 1 × 10−2 and 1.5 mol∙kg−1 HF. The four following diagrams depict the molalities of the dissolved Ag, Cu, Au or Pt species at the respective sampling time t = i (b Me (diss., t = i)) and the respective findings for dissolved Si (b Si (diss., t = i)).
The decrease of b Me (diss., t = i) can be attributed to the reduction of Ag(I), Cu(II), Au(III) or Pt(IV) to the metal state, forming a metal deposit on the silicon wafers. The increase of dissolved silicon in the solutions is based on the oxidation of the silicon and the consecutive dissolution of the oxidized silicon by HF species.
A sectioned exponential trend of the molalities of the species can be observed for all types of metal depositions. This is in accordance with the behavior described in the literature [12,36,37]. Metal deposition in the first phase starts phenotypically at the edges of the silicon wafers, i.e., around the defects created by the cutting of the wafers. It also occurred successively at the wafer surfaces (nucleation). In the second phase, a metal film grows on the entire silicon surface, which appears sponge-like in the silver deposition due to its dendritic growth but appears compact and flat in the other metals’ deposits.
The mathematical separation of the two phases is the intersection t1/2 of the two exponential functions determined, shown in Figure 1a–d. The different exponents of the exponential functions of the first and second phase k1 (s−1) and k2 (s−1), further described as metal deposition rates, were calculated according to Equation (6). The parameter b Me (diss., t1/2) mentioned denotes the molality of the dissolved metal ion at the time t1/2. The value t95% was chosen arbitrarily and corresponds to the time at which the initial dissolved metal ion molality was reduced by 95% and b Me (diss., t95%) is the corresponding value (=0.05 ∙ b Me (diss., t = 0)).
k 1 = ln ( b   M e   ( d i s s . ,   t 1 / 2 ) b   M e   ( d i s s . ,   t   =   0 ) ) · ( t 1 / 2 ) 1   and   k 2 = ln ( b   M e   ( d i s s . ,     t 95   % ) b   M e   ( d i s s . ,   t   =   0 ) ) · ( t 95 % ) 1
It was found that the metal deposition rates of the first phase k1 were about a factor of 2 to 3 below the rates of the second phase k2. This implies that the metal deposits started relatively slowly. The metal deposition accelerated after a mean metal deposition of about 2.0 × 10−6 mol per cm² of wafer surface. Subsequently, the amount of metal deposition per time unit decreased constantly, resulting in an exponential decay of the metal cation molality. The duration of the first deposition phase from t = 0 to t1/2 is less than 360 s for the examples of Ag, Cu and Au shown in Figure 1. It is about 10 times longer for platinum deposition. The second phase, defined as the period from t1/2 to t95%, is about 5 (Ag, Au, Pt) to 10 times (Cu) longer than the first phase.
The kinetics of the metal deposits over both phases k95% were calculated by weighting the metal deposition rates k1 and k2 according to the duration of the two phases according to Equation (7).
k 95 % = k 1 · t 1 / 2 t   = 0 t 95 % + k 2 · t 95 % t 1 / 2 t 95 %
Based on this approach, Ag deposition was found to be the fastest (k95% ≈ 2.0 × 10−3 s−1), followed by Au (k95% ≈ 1.5 × 10−3 s−1), Cu (k95% ≈ 0.5 × 10−3 s−1) and Pt (k95% ≈ 0.2 × 10−3 s−1). It was noted that the amount of the silicon dissolution rate k95% corresponded to the respective metal deposition rate.

3.2. Kinetics of the Metal Depositions at Different Metal Ion and HF Concentration Levels

In Figure 2, the results of Figure 1 were taken (marked by red dots) and supplemented with the results of further batch tests with different initial metal ion and HF molalities. All experiments were performed and evaluated as described in Section 3.1. The results of these experiments are presented as calculated metal deposition rates k95% and plotted against the initial metal ion molalities b Me (diss., t = 0), and the initial redox levels of the Mez+/Me half-cells E (Mez+/Me) (t = 0). The calculation of E (Mez+/Me) (t = 0) was based on the Nernst equation using the standard potentials of the Mez+/Me half-cells [62,63] and the Bromley equations [64] to calculate the activities.
The classification of these results into the overall context reveals that the kinetics of all metal depositions vary considerably. The ranking of the metal deposition kinetics based on the findings in Figure 1 has no general validity. In principle, it can be seen that there is a positive correlation with the initial Ag(I), Cu(II), Au(III) or Pt(IV) molality or the redox level of the Mez+/Me half-cell, starting from the lowest selected metal ion molality up to a certain inflection point. This section is denoted as k95%,1. The functional relationship for silver and copper deposition is logarithmic, and is potential for gold and platinum deposition, as can be derived in the respective equations in Figure 2.
The inflection points identified are located at b Ag (diss., t = 0) = 5.5 × 10−3 mol∙kg−1 or E (Ag+/Ag) = 0.65 V vs. the standard hydrogen electrode (SHE) and at b Cu (diss., t = 0) = 1.0 × 10−2 mol∙kg−1 or E (Cu2+/Cu) = 0.27 V vs. SHE. There is also a narrow drop in the kinetics in copper deposition at b Cu (diss., t = 0) = 2.5 × 10−4 mol∙kg−1 or E (Cu2+/Cu) = 0.23 V vs. SHE. The gold and platinum depositions reveal two inflection points. They are situated at b Au (diss., t = 0) = 3.1 × 10−5 and 1.5 × 10−4 mol∙kg−1 or E (AuCl4/Au) = 0.91 V vs. SHE and 0.92 V vs. SHE, and at b Pt (diss., t = 0) = 5.1 × 10−5 and 6.0 × 10−5 mol∙kg−1, which corresponds to a redox potential of E (PtCl62−/Pt) ≈ 0.68 V vs. SHE.
When the inflection point is exceeded in the silver deposition experimental series, a significant forcing of the kinetics occurs, which is correlated with the initial Ag+ molality or the redox level of the Ag+/Ag half-cell via a potential function (cf. k95%,2 in Figure 2a). The maximum silver deposition rate determined was k95%,max = 6.1 × 10−3 s−1 at b Ag (diss., t = 0) = 1.5 × 10−1 mol∙kg−1. The results of three experiments within this test series do not follow this trend (cf. gray data points in Figure 2a). For Cu deposition experiments, k95%,max is 5.3 × 10−4 s−1 (= result in Figure 1b); for Au, k95%,max is 1.1 × 10−2 s−1; and for Pt, k95%,max is 6.4 × 10−4 s−1. The metal deposition rate above the inflection point k95%,2 in the Cu test series correlates negatively, potentially, with the initial Cu2+ molality or Cu2+/Cu redox level (cf. gray marked points in Figure 2b). The metal deposition rate within the gold deposition series initially drops to a minimum after the first inflection point has been reached.
After exceeding the second inflection point at b Au (diss., t = 0) = 1.5 × 10−4 mol∙kg−1, the metal deposition rate k95%,2 is positively logarithmically linked with the initial AuCl4 molality (cf. in Figure 2c).
The kinetics of platinum deposition are similar to the kinetics of gold deposition up to the second inflection point but at a lower level of k95%,1. In contrast to gold deposition, k95%,2 for platinum deposition decreases negatively logarithmically with an increasing initial Pt(IV) molality (cf. in Figure 2d).

3.3. Discussion of Metal Deposition Kinetics

The different findings for metal deposition kinetics can be derived by considering the state of the initial redox levels of the Mez+/Me half-cells, the valence transfer at the metal–silicon contact and the influence of the HF molality.
Silver: The inflection point of the silver deposition kinetics is an indication of the change of valence transfer between Ag+, the hydrogen-terminated Si surface and the bulk silicon, that can be explained by the simplified theoretical energy scheme shown in Figure 3a.
The illustration contains the range of the redox levels of the Ag+/Ag half-cells E (Ag+/Ag) shown in Figure 2a; the work functions of silver ΦAg = 4.37 eV ± 0.37 eV [65,66,67,68,69] and hydrogen-terminated silicon ΦSi–H = 4.36 eV ± 0.19 eV in vacuum [70]; the Fermi energy EF; and the energetic levels of valence and conduction bands of hydrogen-terminated silicon (Ev (Si–H) = −0.68 eV and Ec (Si–H) = 0.65 eV) [71]. The positions of the valence and conduction bands of bulk silicon (Ev (Sibulk) = −0.41 eV and Ec (Sibulk) = 0.71 eV) were calculated based on the difference of the work function of silicon ΦSi = 4.76 eV ± 0.34 eV [72], the electron affinity of silicon χ = 4.05 ± 0.03 eV [73] and the band gap of 1.12 eV [73].
The work functions of Ag ΦAg and hydrogen-terminated silicon ΦSi–H are, on average, almost identical under a vacuum. When Ag is deposited on Si and the Fermi energy EF at the silver–silicon contact is equalized, there is no significant band bending of the valence and conduction band of the silicon at the silver–silicon contact due to the nearly identical work functions (ΔE = ΦAgΦSi–H = 4.37 eV − 4.36 eV = 0.01 eV).
If it is hypothetically assumed that the difference in the work function of Ag and hydrogen-terminated silicon ΔE = 0.01 eV has an equal effect on the band bending of bulk silicon and hydrogen-terminated silicon, the energy level of the valence band of bulk silicon at the silver–silicon contact Ev (Sibulk/Ag–Si) is −0.40 eV (=Ev (Sibulk) + ΔE = −0.41 eV + 0.01 eV) and that of the valence band of hydrogen-terminated silicon at the silver–silicon contact Ev (Si–H/Ag-Si) is −0.67 eV (Ev (Si–H) + ΔE = −0.68 eV + 0.01 eV). The lowest redox level of the Ag+/Ag half-cell was 0.45 V vs. SHE (cf. Figure 2a). It would, therefore, be high enough to reach the energetic level of the valence band of the bulk silicon at the silver–silicon contact but not that of the valence band of the hydrogen-terminated silicon. Under this consideration, it would be possible to transfer up to four electronic holes (hv+) from Ag+ reduction to the bulk silicon via the valence band mechanism (Equation (8)).
4   Ag +   Ag + 4   h v +   and   Si bulk + 4   h v +     Si 4 +
If the redox level of the Ag+/Ag half-cell exceeds 0.67 V vs. SHE, the energy level of the valence band of the hydrogen-terminated silicon is theoretically reached. This would allow the hole transfer from the Ag+ reduction to the valence band of the hydrogen-terminated silicon. The inflection point of the silver deposition kinetics at E (Ag+/Ag) = 0.65 V vs. SHE marked in Figure 2a deviates only slightly from this theoretical threshold. Based on this consideration, it is assumed that the mean energetic level of the valence band of the hydrogen-terminated silicon at the silver–silicon contact is −0.65 eV. Exceeding this value, the valence transfer from Ag+ reduction to bulk silicon and hydrogen-terminated silicon is enabled. Apparently, it is associated with a significant increase in the silver deposition kinetics (k95%,1 vs. k95%,2).
Assuming the band bending is due to the differences in the work functions between silver and silicon with ΔE = −0.39 eV (=ΦAgΦSi = 4.37 eV − 4.76 eV), an energetic level of −0.80 eV for the valence band of the bulk silicon at the silver–silicon contact should result (=Ev (Sibulk) + ΔE = −0.41 eV + −0.39 eV) and −1.07 eV (Ev (Si–H) + ΔE = −0.68 eV −0.39 eV) for the valence band of the hydrogen-terminated silicon. With the redox level of the Ag+/Ag half-cells between 0.45 V vs. SHE and 0.75 V vs. SHE (Figure 2a), neither valence band is accessible for hole transfer. The process of silver deposition would occur only via the conduction band mechanism. Under the experimental conditions, this mechanism cannot be excluded as the thermally or light-induced excitation of electrons from the valence band into the conduction band of silicon is possible. However, there is no plausible explanation for the inflection point under this consideration. Furthermore, this approach seems to be unlikely, because the higher work function of silicon compared to silver should have resulted in a valence depletion zone on the silver side. Thus, silver would have been the reducing agent and silicon the oxidizing agent. This fundamentally contradicts the experimental findings.
Copper: The argumentation applied to the silver deposition process can also be adapted to copper deposition. If the band bending mechanism is considered based on the differences in the work function of copper (ΦCu = 4.57 eV ± 0.53 eV [65,67,68,69]) and silicon, a valence depletion zone would occur on the copper side due to the higher work function of silicon, at an average of 0.19 eV (=ΦSiΦCu = 4.76 eV − 4.57 eV). Consequently, copper would be the reducing agent and silicon the oxidizing agent. The redox strength of the Cu2+/Cu half-cell would have to exceed a level of at least 0.60 V vs. SHE to match the valence band level of the bulk silicon, which is impossible (E0 (Cu2+/Cu) = 0.34 V vs. SHE [62]).
If the formation of the copper–silicon contact is considered based on a Fermi energy equalization between copper and hydrogen-terminated silicon, the valence and conduction bands of the silicon at the copper–silicon contact are bent by ΔE = 0.21 eV (cf. Figure 3b). The energy level of the valence band of bulk silicon is then −0.20 eV (=Ev (Sibulk) + ΔE = −0.41 eV + 0.21 eV) and that of the hydrogen-terminated silicon is −0.47 eV (=Ev (Si–H) + ΔE = −0.68 eV + 0.21 eV) at the copper–silicon contact. Under the experimental conditions, the Cu2+/Cu half-cells with a redox level between E (Cu2+/Cu) = 0.20 V vs. SHE (≈b Cu (diss., t = 0) ≈ 2.0 × 10−5 mol∙kg−1) and 0.31 V vs. SHE (cf. Figure 2b) would have been strong enough that valence transfer via the valence band mechanism could occur. The valence transfer from the Cu2+/Cu half-cell to the valence band of the hydrogen-terminated silicon would be energetically impossible. Only the Cu+/Cu half-cell (E0 (Cu+/Cu) = 0.521 V vs. SHE [62]) reaches the energetic level of Ev (Si–H) = −0.47 eV.
Two issues are noticeable when examining the kinetics of copper deposition. Firstly, the kinetics at b Cu (diss., t = 0) = 2.5 × 10−4 mol∙kg−1 decrease in a narrow range at a level of k95% = 1.1 × 10−4 s−1 (cf. Figure 2b). This point matches the threshold value of the oxonium ion reaction discussed later in Section 3.6.
The inflection point at b Cu (diss., t = 0) = 1.0 × 10−2 mol∙kg−1 and the associated decrease in the copper deposition rate when it is exceeded are due to the limited silicon dissolution. As mentioned in Equation (5), the dissolution of the oxidized silicon is attributed to either the undissociated HF and/or its anionic dissociation products F and/or HF2 [4,33,36,39,43,53,54,55,56,57,58]. A limitation of HF species should lead to a slower transformation of the oxidized silicon into the soluble product SiF62− and the formation of a more or less thick layer of oxidized silicon, which acts as an insulator [74].
Such a correlation can be mathematically deduced from the activities of the anionic HF dissociation products F, HF2 and H2F3, the amount of silicon dissolution at t95% and the diminished metal deposition rates k95%.
Copper and silver deposition experiments were investigated at three initial levels of total HF concentrations: At b HF (diss., t = 0) ≈ 0.06 mol∙kg−1, 0.3 mol∙kg−1 and 1.5 mol∙kg−1. Based on the additions of CuSO4-HF or AgF/HF-HF selected and the findings of Messnaoui [75] regarding the molalities of the HF dissociation products in relation to the total HF concentration, the initial molalities of the anionic HF dissociation products F, HF2 and H2F3 were calculated for the individual solutions.
The cumulative initial molalities of these species for the silver batch tests were, depending on the AgF/HF and HF additions, between 7.3 × 10−3 and 3.0 × 10−2 mol∙kg−1 for b HF (diss., t = 0) ≈ 0.06 mol∙kg−1; between 2.5 × 10−2 and 1.4 × 10−1 mol∙kg−1 at b HF (diss., t = 0) ≈ 0.3 mol∙kg−1; and between 1.6 × 10−1 and 2.8 × 10−1 mol∙kg−1 at b HF (diss., t = 0) ≈ 1.5 mol∙kg−1. The Cu(II)/HF solutions were based on independent additions of Cu2+ and HF. Therefore, the cumulative initial molalities of F, HF2 and H2F3 were identical: 7.4 × 10−3 mol∙kg−1 for the batches with b HF (diss., t = 0) ≈ 0.06 mol∙kg−1, 2.6 × 10−2 mol∙kg−1 with b HF (diss., t = 0) ≈ 0.3 mol∙kg−1 and 1.6 × 10−1 mol∙kg−1 with b HF (diss., t = 0) ≈ 1.5 mol∙kg−1.
Based on the Bromley equations [64], the mean activity coefficients of the anionic HF species γ (F, HF2 and H2F3) were calculated. The γ (F, HF2 and H2F3) value for AgF/HF solutions was = 0.91 for the solution with the lowest ionic strength of I = 7.3 × 10−3 mol∙kg−1, and γ (F, HF2 and H2F3) was 0.66 for the solution with the highest ionic strength of I = 2.8 × 10−1 mol∙kg−1. In the Cu2+/HF solutions, γ (F, HF2 and H2F3) varied between 0.69 in the lowest concentrated solution (I = 7.5 × 10−3 mol∙kg−1) and 0.05 in the highest concentrated solution (I = 1.8 × 100 mol∙kg−1). The activities of the anionic HF species a (F, HF2 and H2F3) (t = 0) were calculated from the product of the initial molalities and the mean activity coefficients calculated. The results for the individual solutions are plotted against the initial dissolved molalities of Ag+ and Cu2+ in the left half of Figure 4a,b. The molalities of the dissolved silicon at t95% b Si (diss., t95%) are shown in the same diagrams. The quotient of both parameters in relation to the initial molalities of Ag+ and Cu2+ is given in the right half of Figure 4a,b. It is recognizable that more silicon was oxidized and converted into the dissolved form as a result of the increasing initial molality of Ag+ and Cu2+. If the ratio of a (F, HF2 and H2F3) (t = 0) and b Si (diss., t95%) was <6:1 mol∙kg−1:mol∙kg−1, there was a significant decrease of the metal deposition rates k95% compared to the experiments with a (F, HF2 and H2F3) (t = 0):b Si (diss., t95%) ratio of ≥6:1 mol∙kg−1:mol∙kg−1. This criterion found reflects the stoichiometric ratio of F and Si of the formal end product of the oxidized silicon dissolution SiF62−. This leads to the conclusion that the species F, HF2 and H2F3 must be responsible for the rapid dissolution of the oxidized silicon.
The data points with limited F, HF2 and H2F3 activity are marked in gray in Figure 4a,b and in Figure 2a,b for visualization. The comparison of the metal deposition rates k95% shown in Figure 2a,b and the a (F, HF2 and H2F3) (t = 0):b Si (diss., t95%) ratios in Figure 4a,b indicates that the smaller the a (F, HF2 and H2F3) (t = 0):b Si (diss., t95%) ratio was, the more the kinetics of metal deposition diminished. If the metal deposition rates k95% were set as a reference at a (F, HF2 and H2F3) (t = 0): b Si (diss., t95%) ratio of ≥ 6:1 mol∙kg−1:mol∙kg−1, the kinetics of metal deposition at an a (F, HF2 and H2F3) (t = 0):b Si (diss., t95%) ratio of 3:1 mol∙kg−1:mol∙kg−1 at nominally the same initial Ag+ or Cu2+ molality was reduced by about half. At an a (F, HF2 and H2F3) (t = 0):b Si (diss., t95%) ratio of 1:1 mol∙kg−1:mol∙kg−1, the kinetics of metal deposition was only about 1/10 of that without a calculated limitation of the activities of the anionic fluoride species.
The kinetic inhibition of metal deposition due to a restricted access of HF species to the silicon surface caused by the deposited metal film itself, postulated by Geyer et al. [33], does not seem to have been the decisive criterion for the experiments of this study. Otherwise, the metal deposition rates k95% of all metal depositions should have decreased in unison with the increasing initial metal ion molality (cf. Figure 2). Furthermore, it was also observed that the metal layers in all experiments had partially detached from the wafers during the deposition process due to the low adhesion forces between metal and silicon [6,25,26,27] and the continuous stirring of the solutions. As a result, the silicon surface was repeatedly in contact with the HF solution.
Gold: Two inflection points of kinetics have been identified in the gold deposition experiments (cf. Figure 2c). In contrast to the findings for copper, the inflection points were not due to a limitation of the anionic HF dissociation products. The a (F + HF2 + H2F3) (t = 0):b Si (diss., t95%) ratios were above the critical limit of 6:1 mol∙kg−1:mol∙kg−1, with the exception of the experiment with the highest initial Au(III) molality of b Au (diss., t = 0) = 9.3 × 10−2 mol∙kg−1. It is obvious that the kinetics of gold deposition are a mirror image of those of the silver deposition. The functional increase in the gold deposition rate is potential between b Au (diss., t = 0) = 3.0 × 10−7 mol∙kg−1 (E (AuCl4/Au) = 0.87 V vs. SHE) and 3.1 × 10−5 mol∙kg−1 (E (AuCl4/Au) = 0.91 V vs. SHE), analogously to the kinetics of silver deposition above the level of b Ag (diss., t = 0) = 6.0 × 10−3 mol∙kg−1 or E (Ag+/Ag) = 0.65 V vs. SHE (cf. Figure 2a,c). The redox levels of the Mez+/Me half-cells in both cases are strong enough that holes can be transferred from the metal side to the valence band of bulk silicon and hydrogen-terminated silicon (cf. Figure 2a,c and Figure 3a,c). Exceeding the first inflection point, the gold deposition kinetics drop to a minimum at b Au (diss., t = 0) = 1.5 × 10−5 mol∙kg−1. Above this point, the gold deposition rate increased up to a level which was determined for silver below the threshold of b Ag (diss., t = 0) = 6.0 × 10−3 mol∙kg−1 or E (Ag+/Ag) = 0.65 V vs. SHE.
The behavior of gold deposition rates indicates a partial inhibition of the valence transfer processes when the inflection point is exceeded. The amount of band bending at the gold–silicon contact in the gold deposition experiments is the decisive parameter for the kinetics (cf. Figure 3c). According to the theoretical values of the work functions of gold (ΦAu = 4.96 eV ± 0.51 eV [65,67,68,69] and hydrogen-terminated silicon in vacuum, the valence and conduction bands would have to be bent by an average of ΔE = +0.60 eV after Fermi energy equalization at the gold–silicon contact. This means a theoretical energy level of +0.19 eV at the gold–silicon contact (=Ev (Sibulk) + ΔE = −0.41 eV + 0.60 eV) for the valence band of bulk silicon. This provides a potential barrier for the valence transfer between gold and silicon along the Fermi energy level to the valence band of the bulk silicon. It is known as the Schottky barrier Esb. Thus, the valence transfer between AuCl4 and bulk silicon is inhibited, while the valence transfer via the valence band mechanism with the hydrogen-terminated silicon is accessible, since with the band bending, the energy level of the valence band of the hydrogen-terminated silicon at the gold–silicon contact reaches a theoretical energy level of −0.08 eV (=Ev (Si–H) + ΔE = −0.68 eV + 0.60 eV). According to the findings for the gold deposition rates (Figure 2c), the partial inhibition of the valence transfer claimed occurs above an initial AuCl4 molality of b Au (diss., t = 0) = 3.1 × 10−5 mol∙kg−1 and is assumed to be fully effective above an initial AuCl4 molality of b Au (diss., t = 0) = 1.5 × 10−4 mol∙kg−1.
The proposed dependence of the potential barrier on the initial AuCl4 molality means that the amount of band bending at the metal–silicon contact must correlate with b Au (diss., t = 0). Therefore, there is a requirement that the difference in work functions between gold and silicon is also positively correlated with the AuCl4 molality. This is possible because the effective work function is not only a material property. It is influenced by the chemical potential of the supernatant solution, and thus, in these experiments and others, by the AuCl4 molality. The effect of the adsorption of ions on the metal surface becomes relevant if a significant dipole moment is generated at the boundary layer of the metal surface, which impedes the valence transfer from the solution through the metal to the silicon resulting in an enhanced effective work function. Gossenberger et al. [60] have shown, for example, that a fluoride adsorption of the platinum surface causes a strong dipole moment on the metal surface, which increases the effective work function of platinum by about 2 eV when the metal is almost completely covered with fluoride. Based on the kinetics of gold deposition, a similar effect by the adsorption of AuCl4 and/or its reduction intermediates on the gold surface is assumed. Based on the hypothesis that valence transfer to the valence band of the hydrogen-terminated silicon was not affected by a potential barrier, the effective work function of the gold may have been increased by a maximum of ΔE = 0.68 eV (Ev (Si–H) + ΔE = −0.68 eV + 0.68 eV = 0 eV = EF) over the work function of hydrogen-terminated silicon (cf. Figure 3c).
Platinum: The kinetics of platinum deposition (Figure 2d) are quite similar to those of gold deposition up to the inflection points. However, the platinum deposition rates k95%,1 are lower than the analogous gold deposition rates relative to the initial metal ion molalities. This is due to the fact that the PtCl62−/Pt half-cell is weaker than the AuCl4/Au half-cell (cf. Figure 2c,d). The valence bands of the bulk and hydrogen-terminated silicon are hypothetically accessible for valence transfers until the inflection points. In contrast to gold deposition, the platinum deposition rates k95%,2 decrease logarithmically to the initial PtCl62− molality or the redox level of the PtCl62−/Pt half-cell after the second inflection point at b Pt (diss., t = 0) = 6.0 × 10−5 mol∙kg−1 has been exceeded. This would be plausible if Schottky barriers were constituted at the platinum–silicon interface above this level limiting the hole transfer to the valence band of bulk silicon and that of hydrogen-terminated silicon at higher PtCl62− molalities. The energy scheme at the platinum–silicon interface shown in Figure 3d indicates this assumption. The average work function of platinum of ΦPt = 5.53 ± 0.40 eV in vacuum [65,67,68,69] is higher by an average of ΔE = 1.17 eV compared to the average work function of the hydrogen-terminated silicon. The band bending results in an energy level of the valence band of the bulk silicon at the platinum–silicon contact of +0.76 eV (=Ev (Sibulk) + ΔE = −0.41 eV + 1.17 eV). The energy of the valence band of the hydrogen-terminated silicon has a value of +0.49 eV (=Ev (Si–H) + ΔE = −0.68 eV + 1.17 eV). Kuznetsov et al. [47] have experimentally determined a barrier height of up to 0.30 eV at a platinum–silicon contact in HF solution. The effect of the Schottky barrier on the kinetics of metal deposition for platinum is greatest compared to all other metals because platinum has the highest work function of all metals.

3.4. Stoichiometry of Metal Cation Reduction and Silicon Oxidation at Solitary Initial Metal Ion and HF Concentration Levels

In addition to the consideration of the kinetics, the stoichiometric ratio of metal deposition Δb Me and silicon dissolution Δb Si can be determined from the experimental data. The ratio of the decrease of dissolved metal ion molality and that of the increase of silicon molality from t = 0 to t95%b Me: Δb Si (diss., t95%)) was calculated for the four examples shown in Figure 1 according to Equation (9). This calculation approach provided the time-weighted average of the stoichiometric ratio.
Δ b   M e   Δ b   S i ( d i s s . ,   t 95 % ) = b   M e   ( d i s s . , t   =   0 ) b   M e   ( d i s s . , t 95 % ) b   S i   ( d i s s . , t 95 % ) b   S i   ( d i s s . , t   =   0 )
Different Δb Me: Δb Si (diss., t95%) ratios were obtained at roughly the same initial metal ion molalities of about 1 × 10−2 mol∙kg−1. The ratios were Δb Agb Si (diss., t95%) ≈ 3.7:1 mol∙kg−1:mol∙kg−1, Δb Cub Si (diss., t95%) ≈ 1.5:1 mol∙kg−1:mol∙kg−1, Δb Aub Si (diss., t95%) ≈ 1.0:1 mol∙kg−1:mol∙kg−1 and Δb Ptb Si (diss., t95%) ≈ 1.0:1 mol∙kg−1:mol∙kg−1. The valence balance derived from the reaction schemes formulated in Equations (1)–(3) implies that three valences would have been formally transferred in copper and gold deposition, between three and four valences in silver deposition and four valences would have been transferred in platinum deposition during metal ion reduction to the metal state for the oxidation of one silicon atom. Thus, theoretically, between zero and one valence would be missing to complete the balance regarding silicon oxidation up to the oxidation state of +IV in the final product SiF62−. Under this convention, at least one additional oxidant would have to be involved in the silicon oxidation to complete the missing part of the valence transfer. According to the current state of knowledge in the literature, this is attributed to the reduction of hydrogen ions to molecular hydrogen [16,38,39,40,43,49,50,51,52] according to the reaction equations formulated in Equations (2) and (4).
The changes in the molalities of the dissolved metal and silicon species of the experiments presented in Figure 1 were calculated with a different approach based on the results of the consecutive samples analogous to Equation (9). The time-dependent stoichiometric ratios calculated are shown in Figure 5. It is noticeable that the stoichiometric ratios are not constant in the period of t = 0 to t95% and beyond and not similar for the different metal deposits. The ranges of the stoichiometric Δb Meb Si ratios calculated indicate that between two (lower limit for copper deposition) and 4.8 valences (upper limit for platinum deposition) per dissolved silicon atom must have been formally exchanged. As can be deduced by the example of silver deposition above t95%, the stoichiometric ratio of Δb Agb Si ≈ 0:1 mol∙kg−1:mol∙kg−1 also reveals that no valence transfer by metal reduction can cause the silicon oxidation and consecutive silicon dissolution. In conclusion, there must be at least one additional oxidizing agent besides metal cations and at least one additional reducing agent besides silicon in order to achieve a complete valence balance.
Including the data of the other experiments, it can be stated that the metal-specific sequences of the stoichiometric Δb Meb Si ratios are, in most cases, similar to the examples shown in Figure 5.
Silver: The initial stoichiometric Δb Agb Si ratio (between t = 0 and t = 300 s) in the silver deposition experiments was mostly about 75% of the maximum Δb Agb Si ratio. Subsequently, the stoichiometric ratios increased up to the respective maximum. The higher the initial Ag+ molality, the later the maximum of the respective Δb Agb Si ratio was reached. It was achieved at t = 600 s in the example of Figure 5a. Subsequently, the Δb Agb Si ratios decreased mostly below the initial level. The Δb Agb Si ratio in a total of nine experiments dropped sharply after the molality of b Ag (diss., t = i) had decreased to a level below 3.3 × 10−4 mol∙kg−1 (cf. Figure 1a and Figure 5a, between t = 2400 s and t = 3600 s).
Copper: Copper deposition usually started with a Δb Cub Si ratio at about 50% of the maximum Δb Cub Si ratio. Contrary to the silver deposition, the higher the initial Cu2+ molality, the earlier the maximum of Δb Cub Si occurred. It was achieved at t = 1800 s in the example of Figure 5b. Afterwards, the Δb Cub Si ratios decreased again to about the level of the initial stoichiometric ratio. Analogously to the silver experiments, there seems to be a threshold value at b Cu (diss., t = i) ≈ 2.5 × 10−4 mol∙kg−1, below which the Δb Cub Si ratio is significantly reduced. This value corresponds to the finding of the kinetics at which the copper deposition rate k95% has reached a local minimum (cf. Figure 2b).
Gold: Gold depositions started with one exception, as shown in Figure 5c, with the highest Δb Aub Si ratio at t = 300 s. After that, the Δb Aub Si ratios potentially decreased in relation to the time t = i. The experiment with the highest initial AuCl4 concentration of 9.3 × 10−3 mol∙kg−1 shows a reversed temporal sequence of the Δb Aub Si ratios. Analogously to silver and copper deposition, when the Δb Aub Si ratio falls below a threshold of b Au (diss., t = i) ≈ 1.5 × 10−4 mol∙kg−1, there also appears to be a significant decrease in the Δb Aub Si ratio compared to the previous Δb Aub Si ratios (cf. t = 0–t = 2400 s vs. t = 2400 s–t = 3600 s in Figure 1c and Figure 5c). This value corresponds to the second kinetic inflection point of gold deposition marked in Figure 2c.
Platinum: The platinum deposition experiment shown in Figure 5d is the exception in this series. In this example, the Δb Ptb Si ratio started at about 60% of the maximum ratio, decreased at t = 2400 s to the minimum to about 50% of the maximum Δb Ptb Si ratio and then increased potentially to the maximum Δb Ptb Si ratio of 6:5 mol∙kg−1:mol∙kg−1 at t95%. The sequence of the Δb Ptb Si ratios looks mirror-inverted to the abscissa t for the platinum deposits with initial PtCl62− molalities of less than 1.0 × 10−4 mol∙kg−1. The initial Δb Ptb Si ratio started at 40% of the respective maximum ratio, reached its maxima at t = 2400 s and then fell with potential progression to a level below the initial Δb Ptb Si ratio. The stoichiometric Δb Ptb Si ratio decreased significantly in two experiments when the PtCl62− molality fell below a threshold of b Pt (diss., t = i) ≈ 6.0 × 10−4 mol∙kg−1. Similar to gold deposition, this value corresponds to the second kinetic inflection point marked in Figure 2d.

3.5. Stoichiometry of Metal Cation Reduction and Silicon Oxidation at Different Initial Metal Ion and HF Concentration Levels

Figure 6 summarizes the Δb Meb Si ratios of all metal deposits in relation to the initial molality of the metal ion b Me (diss., t = 0) dissolved. The diagrams show the t95%-weighted stoichiometric Δb Meb Si ratios over the period from t = 0 to t95%b Meb Si (diss., t95%)), as shown in Figure 5 as points and the variations of the stoichiometric ratios within this period as bars. The results of Figure 5 have been colored red in Figure 6.
It can be recognized that a higher initial metal ion molality is linked with a larger stoichiometric Δb Meb Si (diss., t95%) ratio, but the Δb Meb Si (diss., t95%) ratios are limited at the lower and upper end of the range of initial metal ion molarities. The Δb Meb Si ratios related to the t95%-weighted findings are just over 0:1 mol∙kg−1:mol∙kg−1 at the initial metal ion molalities of less than 1.0 × 10−5 mol∙kg−1. The upper limits are metal-specific. The t95%-weighted Δb Agb Si ratio for the silver deposition is a maximum of 4.3:1 mol∙kg−1:mol∙kg−1. The analog maximum ratio for the copper deposition is 1.9:1 mol∙kg−1:mol∙kg−1. The maximum Δb Aub Si ratio for gold deposition is about 1.1:1 mol∙kg−1:mol∙kg−1, and the maximum ratio for platinum deposition is about 1.0:1 mol∙kg−1:mol∙kg−1, as shown in Figure 5d.
As discussed in Figure 5, the maxima and minima of the Δb Meb Si ratios between t = 0 and t95% differ significantly from the t95%-weighted ratios. Sigmoidal functional relationships were calculated for the maxima and minima in relation to the initial molalities of the metal ion and are shown in Figure 6. The functions f1 reflect the course of the maxima. The latter was calculated based on Equation (10). The parameter (Δb Meb Si)max included represents the absolute maximum of the stoichiometric ratios of Δb Meb Si. (Δb Agb Si)max is ≈ 4:1 mol∙kg−1:mol∙kg−1 up to b Ag (diss., t = 0) = 5.5 × 10−3 mol∙kg−1 for the silver deposition. Above this threshold, (Δb Agb Si)max is 5:1 mol∙kg−1:mol∙kg−1. There are also two maxima for the copper deposition. The first maximum of (Δb Cub Si)max is 2:1 mol∙kg−1:mol∙kg−1 up to a threshold of b Cu (diss., t = 0) = 8.2 × 10−3 mol∙kg−1. The second maximum of Δb Cub Si)max = 3:1 mol∙kg−1:mol∙kg−1 is located above. Regarding the platinum deposition, f1 converges to a maximum of 1:1 mol∙kg−1:mol∙kg−1 up to b Pt (diss., t = 0) = 5.1 × 10−5 mol∙kg−1 and above that to a maximum of 6:5 mol∙kg−1:mol∙kg−1. There is one maximum (Δb Aub Si)max = 4:3 mol∙kg−1:mol∙kg−1 for the gold deposition. With an approximation of b Me (diss., t = 0) to the value of 0 mol∙kg−1, f1 converges towards a Δb Meb Si ratio of 0:1 mol∙kg−1:mol∙kg−1. The functional increase between the lower and upper limits is described by the empirical term x∙b Me (diss., t = 0)0.01. The validity ranges for f1 regarding the initial metal ion molality b Me (diss., t = 0), the absolute maxima (Δb Meb Si)max and the empirical slope parameter x, can be found in Table 1.
f 1 or   f 2 = Δ b   M e Δ b   S i ( d i s s . , t 95 % ) = b   M e   ( d i s s . ,   t   =   0 )   b   M e   ( d i s s . ,   t   =   0 )   ( Δ b   M e Δ b   S i ) m a x   + x   · b   M e   ( d i s s . ,   t   =   0 ) 0.01
Analogously to the maximum function f1, the functions of the minimas f2 were calculated according to Equation (10). The levels of the Δb Meb Si minima ratios increase similarly to the maxima with increasing initial metal ion molalities b Me (diss., t = 0). The functions f2 converge to the maxima of the functions f1 for copper, gold and platinum deposition (the first maximum of 2:1 mol∙kg−1:mol∙kg−1 for copper deposition, and the second maximum of 6:5 mol∙kg−1:mol∙kg−1 for platinum deposition). The functions f1 and f2 differ in the slope parameter x. By deviation, the minimum function f2 of the silver deposition converges to a lower maximum ratio compared to the f1 function with (Δb Agb Si)max = 10:3 mol∙kg−1:mol∙kg−1. The parameterization of f2 is also shown in Table 1.
The exceeding of the level of b Me (diss., t = 0) between 6.5 × 10−4 mol∙kg−1 (Pt) and 2.4 × 10−3 mol∙kg−1 (Au) resulted in a significant deviation from the minimum functions f2 in the copper, gold and platinum deposition series (but not in the silver deposition). Instead of a convergence of the functions towards the Δb Meb Si maxima ratios of the f1 functions, an inverse sigmoidal course of the Δb Meb Si ratios corresponding to the function f3 according to Equation (11) occurred (parameterization in Table 1). This means that the minima of the Δb Meb Si ratios decreased with higher initial metal ion molalities b Me (diss., t = 0) towards an absolute minimum ratio (Δb Meb Si)min. These minimum stoichiometric ratios are 1:2 mol∙kg−1:mol∙kg−1 for copper and platinum deposition and 1:3 mol∙kg−1:mol∙kg−1 for gold deposition.
f 3 = Δ b   M e Δ b   S i ( d i s s . , t 95 % ) = ( Δ b   M e Δ b   S i ) m a x , f 2 b   M e   ( d i s s . ,   t   =   0 )   b   M e   ( d i s s . ,   t   =   0 )   ( Δ b   M e Δ b   S i ) m a x , f 2 ( Δ b   M e Δ b   S i ) m i n + x · b   M e   ( d i s s . ,   t   =   0 ) 0.01
Based on these results, certain analogous reaction mechanisms can be derived for all metal depositions. The various findings regarding the stoichiometric Δb Meb Si ratios are indications of the involvement of different oxidizing and reducing agents. In connection with the conclusions drawn from the kinetic data about the involvement of bulk and hydrogen-terminated silicon, four sections of the process of metal deposition on silicon can be interpreted, denoted as I, II, III and IV in Figure 6.

3.6. Interpretation of the Stoichiometric Findings

3.6.1. The Involvement of the Oxonium Ion (Section I)

The stoichiometric Δb Meb Si ratios are close to 0:1 mol∙kg−1:mol∙kg−1 for the initial metal ion molalities between 1.0 × 10−7 and ≈ 1.0 × 10−5 mol∙kg−1. Regarding the valence balance, this finding means that silicon oxidation up to the oxidation state of +IV cannot have been achieved by reducing the metal cations alone. It is assumed that the oxonium ions H3O+ resulting from the HF dissociation in water must have been involved in the oxidation process of the silicon up to a certain b Me (diss., t = 0) threshold. A comparison of the standard potential of the 2H3O+/H2 half-cell of 0 V vs. SHE [62] with the energy levels of the valence bands of the bulk and hydrogen-terminated silicon (Ev (Sibulk) = −0.41 eV and Ev (Si–H) = −0.68 eV), outlined in Figure 3, indicates that a valence exchange initiated via the valence band mechanism could never have occurred without the band bending mechanism at the metal–silicon contact. Taking into account the band bending mechanism and the standard potential of the 2H3O+/H2 half-cell of 0 V vs. SHE [62], the oxonium ion reduction would only be possible in the presence of the gold-silicon or platinum–silicon contact (cf. Figure 3c,d). The energy levels of the silicon valence bands at the silver–silicon and copper–silicon contacts are between 0.20 eV (Ev (Sibulk/Cu-Si)) and 0.67 eV (Ev (Si-H/Ag-Si)) above the energy level of the standard 2H3O+/H2 half-cell (cf. Figure 3a,b). From this point of view, the oxonium ion reduction has to be excluded for the silver and copper deposition experiments. However, as a result of the similar Δb Meb Si ratios in all metal deposits tested, the oxonium ion-based process is assumed for all kinds of metal depositions. The supposed paradox can be explained by the underpotential or overvoltage kinetic effect of molecular hydrogen formation on metal surfaces.
If a critical current density J of metal cations on the metal surface is undershot, a positive equilibrium potential is present (underpotential effect), and conversely, a negative equilibrium potential (overvoltage) of the 2H3O+/H2 half-cell on the metal surface is established. The amount of deviation from the standard potential of the 2H3O+/H2 half-cell depends on the metal and the current density. The correlation was determined experimentally by Hunt et al. [76] for Ag, Cu, Au and Pt electrodes. The current density J (A∙cm−1) for the experiments of this study was roughly calculated based on Equation (12). The numerator is the product of the initial metal ion molality b Me (diss, t = 0), the mass of the solution msol (=0.06 kg), the factor 0.95 (because of the meaning of t95%), the charge number z for Ag(I), Cu(II), Au(III) or Pt(IV) and the Faraday constant F (96485 A∙s∙mol−1); the denominator is the product of the metal surface area AMe, which was roughly equated with the silicon wafer surface area ASi (=46 cm²) and t95%.
J = b   M e   ( d i s s . ,   t   =   0 ) · m s o l · 0.95 · z · F A M e · t 95 %
The equilibrium potentials of the 2H3O+/H2 half-cells were derived from the findings of Hunt et al. [76] in comparison with the current densities J calculated. The averaged values of the equilibrium potentials given for bright metal electrodes and metallic colloids in the solutions in this publication served as the basis for the calculation. This approach is in accordance with the observations of the experiments in this study, whereas the metals were deposited on the surface and particles of the metal were present in the solutions due to partial metal detachment from the silicon surface. The equilibrium potentials of the 2H3O+/H2 half-cells determined are shown in Figure 7 in relation to the initial metal ion molalities b Me (diss., t = 0). The equilibrium potentials of the 2H3O+/H2 half-cells can also be found in the energy diagrams in Figure 3. In turn, the findings of Figure 3 regarding the energetic level of the valence bands of the bulk and hydrogen-terminated silicon and the Fermi energy EF at the metal–silicon contact were transferred to Figure 7. In addition, the initial metal ion molalities b Me (diss., t = 0) were drawn in Figure 7, and have the intersection with f1 at Δb Meb Si = 3:z mol∙kg−1:mol∙kg−1 and with f2 at Δb Meb Si = 1:z mol∙kg−1:mol∙kg−1, according to Figure 6. It is assumed that oxonium ion reduction occurs below and up to these threshold values. The range of the oxonium ion effect has been marked with I in Figure 6 and Figure 7.
According to the f1 and f2 functions shown in Figure 6a, the threshold for silver deposition is b Ag (diss., t = 0) = 3.3 × 10−4 mol∙kg−1.This is the value identified at which the Δb Agb Si ratios decreased significantly when it was undershot during silver deposition experiments (cf. Figure 1a and Figure 5a). Figure 7a shows that the equilibrium potentials of the 2H3O+/H2 half-cell below this threshold are above the energy level of the valence band of bulk silicon at the silver–silicon contact but below the energy level of the valence band of hydrogen-terminated silicon (hollow symbols). This means that the oxonium ion reaction with the bulk silicon can occur due to the underpotential effect of the 2H3O+/H2 half-cell on the silver surface. The reaction with the hydrogen-terminated silicon is not possible. If the threshold value mentioned is exceeded, the 2H3O+/H2 half-cell is too weak to oxidize the bulk silicon via the valence band mechanism (full symbols). The stoichiometric Δb Agb Si ratio at the threshold value is a maximum of 3:1 mol∙kg−1:mol∙kg−1 (f1) and a minimum of 1:1 mol∙kg−1:mol∙kg−1 (f2). According to the valence balance, a maximum of three valences were exchanged by the reduction of 3 Ag+ to 3 Ag with one atom of the bulk silicon and one valence was exchanged by the reduction of one oxonium ion to elemental H with the oxidation state ±0. At the minimum, three valences were transferred to one silicon atom of the bulk silicon by oxonium ion reduction and one valence by the reduction of Ag+.
The findings for copper deposition are similar (cf. Figure 7b). The threshold value of b Cu (diss., t = 0) derived from Figure 6b is 2.5 × 10−4 mol∙kg−1. This threshold value was also identified in the copper deposition experiments as the value below which the Δb Cub Si ratio has been significantly diminished. Furthermore, it is the point at which the copper deposition rate has reached an intermediate minimum (cf. Figure 2b). The maximum Δb Cub Si ratio calculated at b Cu (diss., t = 0) = 2.5 × 10−4 mol∙kg−1 was 3:2 mol∙kg−1:mol∙kg−1 (f1) and the minimum Δb Cub Si ratio was 1:2 mol∙kg−1:mol∙kg−1 (f2). Similarly to silver deposition, this means one up to three valences must have been transferred from the oxonium ion reduction forming H (±0), and three up to one valences must have been transferred from Cu2+ reduction to Cu state to one silicon atom of the bulk silicon. According to the higher energetic level of the valence band of the hydrogen-terminated silicon, no reaction with the hydrogen terminated silicon is possible.
The experimental threshold value of b Au (diss., t = 0) observed and calculated for gold deposition is 1.5 × 10−4 mol∙kg−1 (Figure 7c). Similarly to copper deposition, this point corresponds to the intermediate minimum of the gold deposition rate (Figure 2c). The Δb Aub Si intersection point with the f1 function is 1:1 mol∙kg−1:mol∙kg−1 and that of the f2 function is 1:3 mol∙kg−1:mol∙kg−1 (cf. Figure 6c). The valence balance is similar to the silver and copper deposition findings. In contrast to the silver and copper deposition, the reaction of the oxonium ions with the bulk and hydrogen-terminated silicon is possible, due to the relatively strong underpotential effect of the 2H3O+/H2 half-cell on the gold surface and the strong silicon valence band bending at the gold–silicon contact (cf. Figure 3c). Above the threshold value, the equilibrium potential of the 2H3O+/H2 half-cell decreases to <0 V vs. SHE and the oxonium ion-based reaction is no longer possible. Firstly, it is assumed that the hydrogen termination of the silicon is process limiting. It acts hypothetically like a standard 2H3O+/H2 half-cell, and thus limits the valence transfer to the silicon. Secondly, the Fermi energy level of the gold–silicon contact cannot be reached by the 2H3O+/H2 half-cell with a redox level of less than 0 V vs. SHE. This is the essential condition for the transfer of the electronic holes to the silicon via the valence band mechanism.
Basically, the same applies to platinum deposition (Figure 7d). The value observed experimentally and the threshold value for the oxonium ion reaction with silicon determined functionally is b Pt (diss., t = 0) = 6.0 × 10−5 mol∙kg−1 (cf. Figure 6d). This point corresponds to the second inflection point determined in the platinum deposition rate (cf. Figure 2d). At this point, the maximum Δb Ptb Si ratio is 3:4 mol∙kg−1:mol∙kg−1 and the minimum Δb Ptb Si ratio is 1:4 mol∙kg−1:mol∙kg−1. The valence balance for the platinum–silicon system at this threshold is identical to the silver, copper and gold deposition findings at their respective thresholds. Based on this knowledge, the reaction process within section I can be summarized according to the scheme in Equation (13).
( > 0 . . 3 z )   Me z + + ( < 4..1 ) H 3 O + + Si ( > 0 . . 3 z ) Me + ( < 4..1 ) · H + ( < 4..1 ) H 2 O + Si 4 +
Derived from the energy scheme in Figure 7d, the reaction of oxonium ions with the bulk silicon and the reaction with the hydrogen-terminated silicon are possible at the platinum–silicon contact. It is remarkable that this effect is almost exclusively due to the strong silicon valence band bending at the platinum–silicon contact (cf. Figure 3d). The underpotential effect of the 2H3O+/H2 half-cell at the platinum surface is the lowest of all metals. Nevertheless, the effect of the oxonium ions is limited precisely when the equilibrium potential of the 2H3O+/H2 half-cell of the solution falls below the level of 0 V vs. SHE.
The hypothetical limitation of the oxonium ion reduction to the threshold values mentioned was further verified experimentally. If the oxonium ions react with silicon at the metal–silicon contact, the reaction according to Equation (13) would have to produce hydroxyl radicals (∙OH) intermediately according to the scheme shown in Equation (14). The hydroxyl radicals are subsequently converted to H2O by the reaction with the intermediate hydrogen radicals.
Si + 4   H 3 O +     Si 4 + + 8   · H   +   4   · OH     Si 4 + + 2   H 2 + 4   H 2 O
The intermediate formation of hydroxyl radicals is formally impossible in other conceivable reactions of silicon, for example, with H2O or HF. If oxonium ions were the reaction partner, the interception of intermediately formed hydrogen radicals should result in the formation of free hydroxyl radicals. Free hydroxyl radicals are able to react with all other species, including other silicon atoms, due to the very high standard potential of E0 = 2.26 ± 0.48 V vs. SHE [77]. It is assumed that the addition of a hydrogen radical scavenger would cause an increase in the amount of oxidized and dissolved silicon if oxonium ions were be involved in the reaction process.
Such a test was performed for gold and platinum deposition with and without the addition of para-Chlorobenzoic acid (pCBA), which acts as a hydrogen radical scavenger in an acidic solution [78]. This method is less suitable for silver and copper deposition because the less soluble [79] AgCl and CuCl compounds might form partially.
Figure 8 shows the results of the experiments in terms of silicon dissolution after t = 3600 s for gold deposition (b Si (diss., t = 3600 s)) and after t = 14,400 s for platinum deposition (b Si (diss., t = 14,400 s)) in relation to the initial metal ion molalities b Me (diss., t = 0). The results for b Si (diss., t = 3600s) or b Si (diss., t = 14,400 s) without the addition of pCBA were taken from the results shown in Figure 2c,d and 6c,d.
The silicon dissolution was significantly increased in the presence of pCBA compared to tests without the addition of pCBA if the initial metal ion molalities were below the threshold values of b Au (diss., t = 0) = 1.5 × 10−4 mol∙kg−1 or b Pt (diss., t = 0) = 6.0 × 10−5 mol∙kg−1. Above these threshold values, the addition of pCBA had no effect compared with the tests without pCBA. These results confirm that oxonium ions are involved as an oxidizing agent below the threshold values. Above the threshold values, there must be another oxidizing agent.
Coincidence of the oxonium ion reduction thresholds with the significant decreases in copper, gold and platinum deposition rates is deducible (cf. Figure 2b–d). An explanation is the coupling of the reaction of the oxonium ion to the silicon with the silicon dissolution process. For this, the oxonium ion has to react as an ion pair with fluoride H3O+∙F, as has been proposed by Giguère and Turell [80].
The oxonium ion oxidizes the silicon, and the fluoride binds to the oxidized silicon, destabilizing the neighboring Si–Si bonds and promoting further oxidation of the silicon either by further oxonium ion–fluoride pairs or the metal cations, as shown in Figure 9. If the oxonium ion reduction is inhibited, the process of silicon dissolution slows down, and thus, the continuation of the metal cation reduction does too. Above the threshold values, the process of silicon dissolution apparently changes toward the previously postulated anionic HF species F, HF2 and H2F3.

3.6.2. The Divalent and Tetravalent Reaction Mechanism (Section II)

Within the section II, the divalent and tetravalent valence exchange between the metal cations and the bulk silicon is assumed, as has been proposed in the reaction schemes of Equations (1) and (3). There is a slight overlap of section I with section II when comparing the threshold values of b Me (diss., t = 0).
Section II has its formal starting point where two electronic holes are transferred from the metal cation reduction to silicon. According to the charge numbers of the metal cations considered, these points correspond with a stoichiometric Δb Agb Si ratio of 2:1 mol∙kg−1:mol∙kg−1, a Δb Cub Si ratio of 1:1 mol∙kg−1:mol∙kg−1, a Δb Aub Si ratio of 2:3 mol∙kg−1:mol∙kg−1 and a Δb Ptb Si ratio of 1:2 mol∙kg−1:mol∙kg−1. The intersections with the functions f1 shown in Figure 6 are b Ag (diss., t = 0) = 1.3 × 10−4 mol∙kg−1, b Cu (diss., t = 0) = 8.1 × 10−5 mol∙kg−1, b Au (diss., t = 0) = 3.1 × 10−5 mol∙kg−1 and b Pt (diss., t = 0) = 5.1 × 10−5 mol∙kg−1. The gold and platinum deposition thresholds correspond to the first inflection points of the deposition rate k95%, shown in Figure 2c,d.
The tetravalent valence exchange occurs where the stoichiometric Δb Meb Si ratios converge to Δb Agb Si = 4:1 mol∙kg−1:mol∙kg−1, Δb Cub Si = 2:1 mol∙kg−1:mol∙kg−1, Δb Aub Si = 4:3 mol∙kg−1:mol∙kg−1 and Δb Ptb Si = 1:1 mol∙kg−1:mol∙kg−1. These values correspond to the first maxima of the functions f1 (cf. Figure 6 and Table 1). The f1 functions in silver and copper deposition reach these maxima at initial metal ion molalities of b Me (diss, t = 0) ≈ 1.0 × 10−2 mol∙kg−1. However, in gold deposition, section II ends at b Au (diss, t = 0) = 1.5 × 10−4 mol∙kg−1 when the postulated Schottky barrier to the valence band of the bulk silicon becomes effective (cf. Figure 2c and Figure 3c). At this point, the maximum stoichiometric Δb Aub Si ratio is 1:1 mol∙kg−1:mol∙kg−1 (cf. function f1 in Figure 6c). Section II is much narrower in platinum deposition for the same reason. The reaction with the bulk silicon ends at b Pt (diss, t = 0) = 6.0 × 10−5 mol∙kg−1 with a maximum stoichiometric Δb Ptb Si ratio of 3:4 mol∙kg−1:mol∙kg−1 (cf. Figure 2d and Figure 6d).
The main difference in section I is that another oxidant has to fill the gap in the valence balance when fewer than four valences are exchanged between metal cation and silicon. The participation of the oxonium ions is excluded for the reaction process in section II, as discussed in section I. It is assumed that either H2O or a HF species has to be the secondary oxidant. In terms of the concentration, H2O with a molality between b = 54 and 55.5 mol∙kg−1 is more favored compared to b HF (diss., t = 0) = 0.06 mol∙kg−1 up to 1.5 mol∙kg−1. Another indicator for H2O is that there is no significant dependence of metal deposition kinetics on the HF molality in silver and copper deposition experiments without the delay in silicon dissolution. The experiments with the delay of the silicon dissolution also indicate that the metal deposition was limited, probably because of the intermediate formation of silanol groups or oxides which can only be formed by reactions with H2O (cf. Figure 2a,b).
The formulation of the reaction of H2O with silicon under the formation of SiO2, as found in the literature [18,28,38,41,43,81,82], does not appear to be suitable for metal deposition. The redox potential of this reaction E0 = −0.84 V vs. SHE [28,43] is too weak for the valence transfer via the valence band mechanism at the metal–silicon contact. It is more plausible that the silicon is first oxidized to the unstable species Si2+ [83] by the metal cation reduction and that only secondarily does H2O change Si2+ to Si4+ according to Equation (15), as has been previously postulated by Memming and Schwandt [53].
Si 2 + + 2   H 2 O     Si 4 + + 2 · H   +   2   OH     Si 4 + + H 2 + 2   OH
After the formation of Si4+, intermediate silanol groups or oxides are converted to the final product SiF62− by the reaction with the anionic HF dissociation products F, HF2 and H2F3 proposed. The entire reaction process in section II can be summarized with Equation (16).
( 2 z . . 4 z )   Me z + + ( 2..0 )   H 2 O + 6   F ( and   HF 2 and   H 2 F 3 ) + Si bulk ( 2 z . . 4 z ) Me + ( 2..0 ) · H + ( 2..0 ) OH + SiF 6 2

3.6.3. The Change from Bulk-Si Oxidation to Si–Hx Oxidation (Section III)

If a certain threshold value of b Me (diss., t = 0) is exceeded, the Mez+/Me half-cells are able to oxidize the hydrogen-terminated silicon atoms via the valence band mechanism. As discussed for the kinetic evaluation of the experiments, the redox potentials of the Mez+/Me half-cell required depends on the energy level of the valence band of the hydrogen-terminated silicon at the metal–silicon contact. This parameter, in turn, depends on the difference in work functions between the metal and the hydrogen-terminated silicon ΔE, which is further influenced by the dipole moment at the metal-solution interface (cf. Figure 3).
Silver: Regarding silver deposition, this process is enabled at a redox level of the Ag+/Ag half-cell of 0.65 V vs. SHE (cf. Figure 2a and Figure 3a). This corresponds to an initial Ag+ molality of b Ag (diss., t = 0) = 5.5 × 10−3 mol∙kg−1 (cf. Figure 2a and Figure 6a). At this point, the maximum stoichiometric Δb Agb Si ratio is ≈ 4:1 mol∙kg−1:mol∙kg−1 (cf. f1 in Figure 6a). This finding can be explained by an equal oxidation of the Si–H, Si=H2 and Si≡H3 groups according to the reaction scheme in Figure 10a. The maximum stoichiometric Δb Agb Si ratio at higher initial Ag+ molalities converges towards a ratio of 5:1 mol∙kg−1:mol∙kg−1 (cf. f1 in Figure 6a).
This ratio can result if the Si=H2 and Si≡H3 groups but not the Si–H groups were involved in the oxidation process according to the reaction scheme in Figure 10b. This finding relates to the experiments for which a low activity of anionic fluoride species was postulated for the silicon dissolution process (cf. Figure 4a and Figure 6a). This leads to the conclusion that the Si–H terminations could not be dissolved sufficiently quickly after the initial oxidation and remained, for example, as silanol groups. These groups would be passivated against further oxidation by Ag+, so that only the Si=H2 and Si≡H3 groups participated in the further reaction course.
Copper: The oxidation of the hydrogen-terminated silicon by the Cu2+/Cu half-cell is theoretically impossible, because the redox potential of this half-cell is at least 0.1 eV below the level of the valence band of the hydrogen-terminated silicon at the copper–silicon contact (cf. Figure 2b and Figure 3b). Therefore, only the oxidation of the bulk silicon, according to the reaction scheme of Figure 10c, would be plausible (section II). However, the Δb Cub Si ratio in the experiments with initial Cu2+ molalities of b Cu (diss, t = 0) ≥ 8.2 × 10−3 mol∙kg−1 exceeded the stoichiometric ratio of 2:1 mol∙kg−1:mol∙kg−1 and converged to a maximum ratio of 3:1 mol∙kg−1:mol∙kg−1 (cf. f1 in Figure 6b). This result is possible under two conditions. The first requirement is that an accumulation of Cu(I) species must have occurred during the experiments. The second requirement is that only the Si≡H3 groups have been oxidatively attacked. The redox potential of the Cu+/Cu half-cell is higher than that of the Cu2+/Cu half-cell [62]. With a Cu+ molality of 8.2 × 10−3 mol∙kg−1, the Cu+/Cu half-cell reaches a level of 0.40 V vs. SHE. Hypothetically, this is the energetic level from which the valence transfer via the valence band of the hydrogen-terminated silicon is enabled. Based on the energy scheme shown in Figure 3b, it would be possible if the amount of band bending ΔE at the copper–silicon contact is 0.28 eV and not 0.21 eV, resulting in an energy level of the valence band of the hydrogen-terminated silicon at the copper–silicon contact of −0.40 eV. The sufficient accumulation of the intermediate Cu+ would be conceivable if the reduction from Cu2+ to metal state had been limited. Such a limitation was proven for the respective experiments due to delayed silicon dissolution (cf. Figure 4b). It is further assumed that especially the silicon atoms of the oxidized Si–H and Si=H2 groups were not dissolved sufficiently fast and Si−OH and Si=(OH)2 groups were formed. As a result, only the Si≡H3 groups remained freely accessible. The proposed reaction scheme based on these considerations is shown in Figure 10d.
Gold: Above an initial AuCl4 molality of b Au (diss., t = 0) = 1.0 × 10−2 mol∙kg−1, the maximum Δb Aub Si ratio in gold deposition converged towards a ratio of 4:3 mol∙kg−1:mol∙kg−1. From this finding, the reaction scheme shown in Figure 10e can be concluded, according to which the Si–H, Si=H2 and Si≡H3 groups are oxidized stoichiometrically uniformly. The process did not start at this ratio in the experiment with the highest initial AuCl4 molality of b Au (diss., t = 0) = 9.3 × 10−2 mol∙kg−1, but at Δb Aub Si ≈ 1:3 mol∙kg−1:mol∙kg−1, and increased subsequently to the maximum ratio of 4:3 mol∙kg−1:mol∙kg−1.
This deviating finding can be explained by the fact that the silicon dissolution was limited under the conditions of this experiment. The a (F + HF2 + H2F3) (t = 0):b Si (diss., t95%) ratio was below the critical level of 6:1 mol∙kg−1:mol∙kg−1 at ≈ 1:1 mol∙kg−1:mol∙kg−1 proposed. It is assumed that hydrogen-terminated silicon groups were not sufficiently available at the beginning of the experiment due to the retarded silicon dissolution and silanol group formation. Consequently, the valence transfer shifted to the Si–Si backbonds of the near-surface silicon atoms, according to the reaction scheme in Figure 10f. It is supposed that the activity of the anionic HF species increased with decreasing AuCl4 molality. As a result, the silicon dissolution and re-hydrogen termination process became faster, so that the reaction could switch to the reaction scheme of Figure 10e.
Platinum: The maximum Δb Ptb Si ratio of 6:5 mol∙kg−1:mol∙kg−1 shown in Figure 6d indicates that the Si–H groups were not oxidized, but the Si=H2 and Si≡H3 groups were, according to the reaction scheme in Figure 10g. Unlike the other metal depositions, this is not caused by delayed silicon dissolution. It is assumed that the valence transfer to the Si–H groups is inhibited by a Schottky barrier. According to the calculated findings of Papaconstantopoulos and Economou [84], the energy levels of the valence bands of the silicon with different hydrogen terminations vary. The less hydrogen bound to the silicon, the lower the energy level of the valence band of the silicon. Due to the high difference in work functions between platinum and hydrogen terminated silicon (cf. Figure 3d), the valence band of the silicon of the Si–H termination at the platinum–silicon contact was probably bent to just above the level of the Fermi energy, so that the Schottky barrier was constituted to this group but not to the Si=H2 and Si≡H3 groups.
Analogously to gold deposition, a Δb Ptb Si ratio of ≈ 3:5 mol∙kg−1:mol∙kg−1 was determined for the platinum deposition experiment with the highest initial PtCl62− molality of b Pt (diss., t = 0) = 9.7 × 10−3 mol∙kg−1 at the beginning of the experiment, which later increased to the maximum ratio of 6:5 mol∙kg−1:mol∙kg−1 (cf. Figure 5d). This reaction process also indicates a temporarily slower silicon dissolution process. The a (F + HF2 + H2F3) (t = 0):b Si (diss., t95%) ratio of 9:1 mol∙kg−1:mol∙kg−1 was just above the critical ratio of 6:1 mol∙kg−1:mol∙kg−1 proposed. In contrast to the gold deposition process sketched in Figure 10f, only a part of the oxidized silicon on the surface was temporarily not dissolved fast enough in the platinum deposition experiment. It is assumed, therefore, that some of the valences resulting from the PtCl62− reduction were transferred to the hydrogen-terminated silicon and the others to the Si-Si backbonds of the surface-near silicon. A mixed reaction can be interpreted from the initially proven Δb Ptb Si ratio of about 3:5 mol∙kg−1:mol∙kg−1, according to the scheme of Figure 10h, which must have changed in the course of the experiment to the reaction scheme postulated in Figure 10g.

3.6.4. The Decrease in Δb Me:Δb Si Ratios toward a Minimum (Section IV)

The reaction processes according to section III are followed by a successive decrease in the stoichiometric Δb Meb Si ratios to the level of the f2 function for silver and to the level of the f3 functions for the other metal deposits (cf. Figure 6). The Δb Agb Si ratio converges to a minimum of 10:3 mol∙kg−1:mol∙kg−1 (cf. f2 in Figure 6a, Table 1). The minima of the other Δb Meb Si ratios decrease the higher the initial metal ion molality was and tended towards a minimum ratio of Δb Cub Si = 1:2 mol∙kg−1:mol∙kg−1, Δb Aub Si = 1:3 mol∙kg−1:mol∙kg−1 and Δb Ptb Si = 1:2 mol∙kg−1:mol∙kg−1 (cf. f3 in Figure 6b–d, Table 1).
This finding can be interpreted as follows: The initial reactions of the metal cations according to section III lead to the oxidization of the upper Si–H, Si=H2 and Si≡H3 groups and their consecutive dissolution (cf. Figure 10). Through the involvement of H2O and the anionic HF species, the silicon underneath the previously removed silicon will be re-terminated with hydrogen. The process step of metal cation reduction and silicon oxidation then starts again. Through the multiple sequence of the process, the metal successively penetrates the surface of the silicon. The consecutive reaction by H2O is becoming increasingly important due to the successive decreasing molality of the metal cations. The change of the reaction course is metal-dependent as the different minimum ratios of Δb Meb Si indicate.
Silver: The minimum Δb Agb Si ratio of 10:3 mol∙kg−1:mol∙kg−1 in silver deposition leads to the conclusion that by reducing 10 Ag+, formally one Si–H and two Si=H2 groups have to be oxidized according to the reaction scheme in Figure 11a. This means that the Si≡H3 groups attacked initially (Figure 10a) would not be further oxidized with decreasing Ag+ molality. The shift towards a selective oxidation of the different hydrogen-terminated silicon groups can be explained by the fact that the energetic level of the Ag+/Ag half-cell decreases with the successive consumption of Ag+. If the Ag+ molality falls below a certain level, the energetic level of the valence band of the silicon of the Si≡H3 group can obviously no longer be reached. Based on the experimental findings, the level of Ev (Si≡H3) should be about −0.70 eV (cf. Figure 3a). The energy level corresponds to that of the Ag+/Ag half-cell at an Ag+ activity of 4.5 × 10−2 mol∙kg−1.
Copper: The determined minimum of the ratio Δb Cub Si of 1:2 mol∙kg−1:mol∙kg−1 indicates a change from the oxidation of the Si≡H3 groups via Cu+ reduction (cf. Figure 10d) to the oxidation of the Si–H groups formed after the removal of the Si≡H3 groups (cf. Figure 11b).
Gold: The trend of gold deposition from a maximum Δb Aub Si ratio of 4:3 mol∙kg−1:mol∙kg−1 toward a minimum Δb Aub Si ratio of 1:3 mol∙kg−1:mol∙kg−1 suggests that after the initial attack on all hydrogen-terminated silicon groups, according to the scheme of Figure 10e, the reaction has shifted successively to the Si–H and Si=H2 groups (Δb Aub Si = 1:1 mol∙kg−1:mol∙kg−1), and furthermore, exclusively to the Si–H groups (Δb Aub Si = 2:3 mol∙kg−1:mol∙kg−1) with decreasing AuCl4 molality. After the AuCl4 molality has decreased to the level of 1.5 × 10−4 mol∙kg−1, the bulk silicon would also have been accessible to oxidation again. The proven minimum Δb Aub Si ratio of 1:3 mol∙kg−1:mol∙kg−1 could then result from the reaction with the Si–H termination and the underlying bulk silicon, according to the reaction scheme outlined in Figure 11c. The reaction process assumed is associated with a high formation of molecular hydrogen compared to the other metal deposits. According to the optical impression, the gas formation during gold deposition was strongest compared to the other metal depositions.
Platinum: The process of platinum deposition is only partially similar to gold deposition. The oxidation of the Si=H2 and Si≡H3 groups is preferred in section III (cf. Figure 10g). As the platinum deposition penetrates the silicon, the Δb Ptb Si ratio decreases to a minimum of 1:2 mol∙kg−1:mol∙kg−1. This leads to the conclusion that the reaction path changes towards the Si–H groups with decreasing PtCl62− molality, as shown in Figure 11d. Only twofold valence transfer formally occurs in the reduction of Pt(IV) via the intermediate Pt(II) state to the Pt state. This has the consequence that no molecular hydrogen can be produced whose formation requires onefold valence transfers. This conclusion corresponds to the observation of the experimental course of platinum deposition, after which no gas formation was observable.

4. Conclusions

Based on batch experiments, the reductions of Ag+, Cu2+, AuCl4 and PtCl62− on multi-crystalline silicon in a HF matrix was investigated regarding kinetics and stoichiometry.
Metal deposition started slowly and almost linearly in time in all experiments. After exceeding a metal deposition of about 2.0 × 10−6 mol per cm2; silicon surface, the subsequent metal deposition resulted in an exponential decay. The kinetics of both phases differed by a factor of 2 to 3. This result corresponds to findings described in the literature [12,36,37]. The reason was the change from the initial metal nucleation to the three-dimensional growth of the metal layer on the silicon.
The kinetics of metal deposition correlate positively with the molalities of Ag+, Cu2+, AuCl4 and PtCl62−. The decisive criteria were the overlap of the energetic level of the redox potential of the Mez+/Me half-cells involved with the valence band of the bulk and hydrogen-terminated silicon at the metal–silicon contact and a possibly existing Schottky barrier behavior, as has been postulated previously in [20,21,22,23,24,33,46,48].
The metal deposition kinetics for silver increased logarithmically with higher initial Ag+ molality in the energetic overlap range of the Ag+/Ag half-cell with the valence band of the bulk silicon. The energetic level of the valence band of the hydrogen-terminated silicon with a redox level of E (Ag+/Ag) = 0.65 V vs. SHE (b Ag (diss., t = 0) = 5.5 × 10−3 mol∙kg−1) was presumably reached, and the kinetics of silver deposition correlated potentially with the initial Ag+ molality. The maximum metal deposition rate was found at the initial Ag+ molality of b Ag (diss., t = 0) of 1.5 × 10−1 mol∙kg−1.
The kinetics of gold deposition were mirror-inverted compared to the findings for silver deposition. The kinetics of gold deposition correlated potentially with the initial AuCl4 molality up to an initial AuCl4 molality of b Au (diss., t = 0) of 3.1 × 10−5 mol∙kg−1 and reached the maximum at this threshold. Above it, the kinetics of gold deposition decreased to a minimum at b Au (diss., t = 0) of 1.5 × 10−4 mol∙kg−1 and increased above logarithmically. It is suspected that the valence band of the bulk silicon at the threshold is bent above the energetic level of the Fermi energy. Thus, a Schottky barrier inhibited the reaction path with the bulk silicon. The valence transfer shifted above the threshold to the hydrogen-terminated silicon.
Based on the findings of silver and gold deposition kinetics, it can be deduced that the kinetics of the reactions with bulk and hydrogen-terminated silicon are similar. The proposed faster reaction with the silicon atoms of the bulk silicon compared to the hydrogen-terminated silicon atoms [18,24,26,29,41] cannot be confirmed according to the results of this study.
A partial Schottky barrier behavior could be derived for platinum deposition similar to gold deposition if an initial PtCl62− molality of b Pt (diss., t = 0) = 5.1 × 10−5 mol∙kg−1 was exceeded. The Schottky barrier affected the reaction with the bulk silicon and also the reaction with the Si–Hx groups. The kinetic inhibition of platinum deposition was more pronounced compared to gold deposition. This can be attributed to the stronger band bending effect and Schottky barrier height due to the comparatively greater work function of platinum [65,67,68,69].
The investigation of the kinetics of copper deposition revealed a kinetic inhibition of copper deposition in the case of limited silicon dissolution. The decrease in the kinetics of copper deposition occurred when the activities of the HF dissociation products F, HF2 and H2F3 had decreased below a ratio of 6:1 mol∙kg−1:mol∙kg−1 in relation to the molality of silicon dissolution. The criterion of limitation found reflects the stoichiometric ratio of the final product SiF62−. Another relationship between other HF species, as insinuated in the literature in different variants [13,33,36,39,43,53,54,55,56,57,58,59], could not be verified.
The stoichiometric ratios of metal deposition determined in relation to silicon dissolution Δb Meb Si indicate a complex reaction process. The stoichiometric Δb Meb Si ratios increased with rising metal cation molalities within a sigmoidal functional relationship. At the lowest level of b Me (diss., t = 0), the Δb Meb Si ratios tended to ratios shortly above 0:1 mol∙kg−1:mol∙kg−1, indicating the involvement of at least one further oxidizing agent.
The assumption of a metal cation reduction with coexisting hydrogen ion reduction via the silicon conduction band mechanism [16,38,39,40] was not confirmed. The results and calculations indicate an oxonium ion reduction with valence transfer via the valence band of silicon at the metal–silicon contact. In the case of a sufficiently low molality of the metal cation, an underpotential effect of the 2H3O+/H2 half-cell is formed on the metal surface. If the resulting equilibrium potential of the 2H3O+/H2 half-cell reaches at least the energy level of the valence band of the bulk silicon, and at most, the level of the Fermi energy at the metal–silicon contact, the valence transfer from the oxonium ion to the silicon can occur via the valence band mechanism. Therefore, the molality of b Ag (diss., t = 0) = 3.3 × 10−4 mol∙kg−1, b Cu (diss., t = 0) = 2.5 × 10−4 mol∙kg−1, b Au (diss., t = 0) = 1.5 × 10−4 mol∙kg−1 or b Pt (diss., t = 0) = 6.0 × 10−5 mol∙kg−1 must not be exceeded. Otherwise, the process of oxonium ion reduction at the metal–silicon contact via the valence band mechanism is impossible.
Above these thresholds, the oxidation of the silicon will be caused by metal cation and consecutive water reduction. The valence exchange between metal cation and silicon is divalent at low molalities. According to the charge numbers of the Ag(I), Cu(II), Au(III) and Pt(IV) species, the Δb Meb Si ratios are Δb Agb Si = 2:1 mol∙kg−1:mol∙kg−1, Δb Cub Si = 1:1 mol∙kg−1:mol∙kg−1, Δb Aub Si = 2:3 mol∙kg−1:mol∙kg−1 and Δb Ptb Si = 1:2 mol∙kg−1:mol∙kg−1. The missing valences for the oxidation of the silicon up to the oxidation state of +IV are produced by the reduction of water. At higher molalities, the proportion of tetravalent valence transfer increases and the proportion of water in the reaction process decreases. The functional relationship between the initial molalities of Ag+, Cu2+, AuCl4 and PtCl62− and stoichiometric Δb Meb Si ratios is sigmoidal. Above metal cation molalities of ≈ 1.0 × 10−2 mol∙kg−1, the valence exchange is tetravalent. At this level, the Δb Meb Si ratios reach their maxima. With the exception of gold deposition, the stoichiometric ratios are not twice as high as in the divalent reaction range. The maximum Δb Agb Si ratio was 5:1 mol∙kg−1:mol∙kg−1. The maximum Δb Cub Si ratio amounted to 3:1 mol∙kg−1:mol∙kg−1, and the maximum Δb Ptb Si ratio reached a level of 6:5 mol∙kg−1:mol∙kg−1. The valence balance implies the involvement of the hydrogen of the Si–Hx groups as a further reducing agent, because the oxidation of silicon above an oxidation state of +IV is not plausible. The criterion that the redox level of the Mez+/Me half-cells has to be at least equal to the energetic levels of the valence bands of the Si–Hx species is achieved for the initial molalities of b Ag (diss., t = 0) ≥ 5.5 × 10−3 mol∙kg−1. Much lower initial metal ion molalities are required in the gold and platinum deposition process, mainly due to the strong silicon valence band bending at the metal–silicon contact. If the threshold values of b Au (diss., t = 0) = 3.1 × 10−5 mol∙kg−1 and b Pt (diss., t = 0) = 5.1 × 10−5 mol∙kg−1 are exceeded, the valence transfer occurs exclusively to the valence band of the hydrogen-terminated silicon because the valence transfer to the bulk silicon atoms is inhibited by the constitution of a Schottky barrier. In copper deposition, the energetic level of the valence band of Si–Hx groups at the copper–silicon contact of –0.40 eV can be achieved by the Cu+/Cu half-cell at b Cu (diss., t = 0) = 8.2 × 10−3 mol∙kg−1. This process became significant in the experiments with limited silicon dissolution. The valence balance which can be derived from the maximum stoichiometric ratio of Δb Cub Si = 3:1 mol∙kg−1:mol∙kg−1 indicates that the reduction of Cu+ to Cu must have proceeded mainly with the oxidation of the Si≡H3 groups (→ 6 h+). The Δb Agb Si ratio of 5:1 mol∙kg−1:mol∙kg−1 for silver suggests an equal oxidation of the Si=H2 (→ 4 h+) and Si≡H3 (→ 6 h+) groups. In both cases, the reaction path was forced by the limited dissolution of previously oxidized bulk silicon and oxidized silicon of the former Si–H and partially Si=H2 groups induced by the low activity of F, HF2 and H2F3. The Δb Aub Si ratio of 4:3 mol∙kg−1:mol∙kg−1 for gold is due to an equal oxidation of all Si–Hx groups (Si–H (→ 2 h+), Si=H2 (→ 4 h+) and Si≡H3 (→ 6 h+)). The Δb Ptb Si ratio of 6:5 mol∙kg−1:mol∙kg−1 can be explained by the assumption that the oxidation of the Si–H groups is inhibited by the constitution of a Schottky barrier, and thus, only the valence exchange with the Si=H2 and Si≡H3 groups can occur.

Author Contributions

Conceptualization, S.S. and J.A.; methodology, S.S.; validation, J.A.; formal analysis, S.S.; investigation, S.S.; resources, S.S.; data curation, S.S.; writing—original draft preparation, S.S.; writing—review and editing, J.A.; visualization, S.S.; supervision, J.A.; project administration, J.A.; funding acquisition, S.S. and J.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Deutsche Forschungsgemeinschaft, grant number 314964421. The APC was funded by Brandenburg University of Technology Cottbus-Senftenberg.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Li, J.; Van Toan, N.; Wang, Z.; Ono, T. Metal-assisted-chemical-etching of silicon nanowires for templating 3D graphene growth towards energy storage in microsystems. J. Micromech. Microeng. 2019, 29, 055007. [Google Scholar] [CrossRef]
  2. Nichkalo, S.; Druzhinin, A.; Yerokhov, V.; Ostapiv, O. Fabrication and Characterization of High-Performance Anti-reflecting Nanotextured Si Surfaces for Solar Cells. In Nanooptics, Nanophotonics, Nanostructures, and Their Applications; Fesenko, O., Yatsenko, L., Eds.; NANO 2017, Springer Proceedings in Physics, 210; Springer: Cham, Switzerland, 2017. [Google Scholar] [CrossRef]
  3. Stafiniak, A.; Prazmowska, J.; Macherzynski, W.; Paszkiewicz, R. Nanostructuring of Si substrates by a metal assisted chemical etching and dewetting process. RSC Adv. 2018, 8, 31224–31230. [Google Scholar] [CrossRef] [Green Version]
  4. Zheng, H.; Han, M.; Zheng, P.; Zheng, L.; Qin, H.; Deng, L. Porous silicon templates prepared by Cu-assisted chemical etching. Mat. Lett. 2014, 118, 146–149. [Google Scholar] [CrossRef]
  5. Kuntyi, O.; Zozulya, G.I.; Shepida, M.V.; Nichkalo, S.I. Deposition of nanostructured metals on the surface of silicon by galvanic replacement: A mini-review. Voprosy Khimii Khimicheskoi Tekhnologii 2019, 3, 74–82. [Google Scholar]
  6. Nagahara, L.A.; Ohmori, T.; Hashimoto, K.; Fujishima, A. The influence of hydrofluoric-acid concentration on electroless copper deposition onto silicon. J. Electroanal. Chem. 1992, 333, 363–369. [Google Scholar] [CrossRef]
  7. Homma, T.; Wade, C.P.; Chidsey, C.E.D. Nucleation of trace copper on the H-Si (111) surface in aqueous fluoride solutions. J. Phys. Chem. B 1998, 102, 7919–7923. [Google Scholar] [CrossRef]
  8. Hsu, C.-C.; Shen, F.-Y.; Huang, F.-S. Gold nanowires fabricated by immersion plating. Nanotechnology 2008, 19, 195302. [Google Scholar] [CrossRef]
  9. Milazzo, R.G.; Mio, A.M.; D’Arrigo, G.; Grimaldi, M.G.; Spinella, C.; Rimini, E. Coalescence of silver clusters by immersion in diluted HF solution. J. Chem. Phys. 2015, 143, 024306. [Google Scholar] [CrossRef]
  10. Smith, Z.R.; Smith, R.L.; Collins, S.D. Mechanism of nanowire formation in metal assisted chemical etching. Electrochim. Acta 2013, 92, 139–147. [Google Scholar] [CrossRef]
  11. Norga, G.J.; Platero, M.; Black, K.A.; Reddy, A.J.; Michel, J.; Kimerling, L.C. Mechanism of copper deposition on silicon from dilute hydrofluoric acid solution. J. Electrochem. Soc. 1997, 144, 2801–2810. [Google Scholar] [CrossRef]
  12. Lee, M.K.; Wang, H.D.; Wang, J.J. A Cu seed layer for Cu deposition on silicon. Solid State Electron. 1997, 41, 695–702. [Google Scholar] [CrossRef]
  13. Srinivasan, R.; Suni, I.I. Kinetic analysis of Au deposition from aqueous HF onto Si (111) by surface second harmonic generation. J. Electrochem. Soc. 1999, 146, 570–573. [Google Scholar] [CrossRef]
  14. Gorji, S.G.; Razak, K.A.; Cheong, K.Y. Gold nanoparticles deposited on linker-free silicon substrate and embedded in aluminum Schottky contact. J. Colloid Interface Sci. 2013, 408, 220–228. [Google Scholar] [CrossRef] [PubMed]
  15. Chyan, O.M.R.; Chen, J.-J.; Chien, H.Y.; Sees, J.; Hall, L. Copper deposition on HF etched silicon surfaces: Morphological and kinetic studies. J. Electrochem. Soc. 1996, 143, 92–96. [Google Scholar] [CrossRef]
  16. Brahiti, N.; Bouanik, S.-A.; Hadjersi, T. Metal-assisted electroless etching of silicon in aqueous NH4HF2 solution. Appl. Surf. Sci. 2012, 258, 5628–5637. [Google Scholar] [CrossRef]
  17. Scudiero, L.; Fasasi, A.; Griffiths, P.R. Characterization of a controlled electroless deposition of copper thin film on germanium and silicon surfaces. Appl. Surf. Sci. 2011, 257, 4422–4427. [Google Scholar] [CrossRef]
  18. Peng, K.Q.; Hu, J.J.; Yan, Y.J.; Wu, Y.; Fang, H.; Xu, Y.; Lee, S.T.; Zhu, J. Fabrication of single-crystalline silicon nanowires by scratching a silicon surface with catalytic metal particles. Adv. Funct. Mater. 2006, 16, 387–394. [Google Scholar] [CrossRef]
  19. Huang, Z.; Geyer, N.; Werner, P.; de Boor, J.; Gösele, U. Metal-assisted chemical etching of silicon: A review. Adv. Mater. 2011, 23, 285–308. [Google Scholar] [CrossRef]
  20. Li, S.; Ma, W.; Zhou, Y.; Chen, X.; Xiao, Y.; Ma, M.; Zhu, W.; Wie, F. Fabrication of porous silicon nanowires by MACE method in HF/H2O2/AgNO3 system at room temperature. Nanoscale Res. Lett. 2014, 9, 196–203. [Google Scholar] [CrossRef] [Green Version]
  21. Milazzo, R.G.; Mio, A.M.; D’Arrigo, G.; Smecca, E.; Alberti, A.; Fisichella, G.; Giannazzo, F.; Spinella, C.; Rimini, E. Influence of hydrofluoric acid treatment on electroless deposition of Au clusters. Beilstein J. Nanotechnol. 2017, 8, 183–189. [Google Scholar] [CrossRef] [Green Version]
  22. Warren, S.; Reitzle, A.; Kazimirov, A.; Ziegler, J.C.; Bunk, O.; Cao, L.X.; Renner, F.U.; Kolb, D.M.; Bedzyk, M.J.; Zegenhagen, J. A structure study of the electroless deposition of Au on Si (111): H. Surf. Sci. 2002, 496, 287–298. [Google Scholar] [CrossRef]
  23. Yae, S.; Nasu, N.; Matsumoto, K.; Hagihara, T.; Fukumuro, N.; Matsuda, H. Nucleation behavior in electroless displacement deposition of metals on silicon from hydrofluoric acid solutions. Electrochim. Acta 2007, 53, 35–41. [Google Scholar] [CrossRef]
  24. Han, H.; Huang, Z.; Lee, W. Metal-assisted chemical etching of silicon and nanotechnology applications. Nano Today 2014, 9, 271–304. [Google Scholar] [CrossRef] [Green Version]
  25. Nagahara, L.A.; Ohmori, T.; Hashimoto, K.; Fujishima, A. Effects of HF solution in the electroless deposition process on silicon surfaces. J. Vac. Sci. Technol. 1993, 11, 763–767. [Google Scholar] [CrossRef]
  26. dos Santos Filho, S.G.; Pasa, A.A.; Hasenack, C.M. A mechanism for electroless Cu plating onto Si. Microelectron. Eng. 1997, 33, 149–155. [Google Scholar] [CrossRef]
  27. Magagnin, L.; Maboudian, R.; Carraro, C. Gold deposition by galvanic displacement on semiconductor surfaces: Effect of substrate on adhesion. J. Phys. Chem. B 2002, 106, 401–407. [Google Scholar] [CrossRef]
  28. Morinaga, H.; Suyama, M.; Ohmi, T. Mechanism of metallic particle growth and metal-induced pitting on Si wafer surface in wet chemical processing. J. Electroanal. Chem. 1994, 141, 2834–2841. [Google Scholar] [CrossRef]
  29. Peng, K.Q.; Fang, H.; Hu, J.; Wu, Y.; Zhu, J.Z.; Yan, Y.; Lee, S.T. Metal-particle-induced, highly localized site-specific etching of Si and formation of single-crystalline Si nanowires in aqueous fluoride solution. Chem. Eur. J. 2006, 12, 7942–7947. [Google Scholar] [CrossRef]
  30. Scheeler, S.P.; Ullrich, S.; Kudera, S.; Pacholski, C. Fabrication of porous silicon by metal-assisted etching using highly ordered gold nanoparticle arrays. Nanoscale Res. Lett. 2012, 7, 450–456. [Google Scholar] [CrossRef] [Green Version]
  31. Lajvardi, M.; Eshghi, H.; Izadifard, M.; Ghazi, M.E.; Goodarzi, A. Effects of silver and gold catalytic activities on the structural and optical properties of silicon nanowires. Phys. E 2016, 75, 136–143. [Google Scholar] [CrossRef]
  32. Dawood, M.K.; Tripathy, S.; Dolmanan, S.B.; Ng, T.H.; Tan, H.; Lam, J. Influence of catalytic gold and silver metal nanoparticles on structural, optical, and vibrational properties of silicon nanowires synthesized by metal-assisted chemical etching. J. Appl. Phys. 2012, 112, 073509. [Google Scholar] [CrossRef]
  33. Geyer, N.; Fuhrmann, B.; Huang, Z.; de Boor, J.; Leipner, H.S.; Werner, P. Model for the mass transport during metal-assisted chemical etching with contiguous metal films as catalysts. J. Phys. Chem. C 2012, 116, 13446–13451. [Google Scholar] [CrossRef]
  34. Zhang, X.; Zhu, J.; Huang, X.; Qian, Q.; He, Y.; Chi, L.; Wang, Y. Controllable and facile fabrication of gold nanostructures for selective metal-assisted etching of silicon. Small 2014, 10, 2451–2458. [Google Scholar] [CrossRef] [PubMed]
  35. Ozdemir, B.; Kulakci, M.; Turan, R.; Unalan, H.E. Effect of electroless etching parameters on the growth and reflection properties of silicon nanowires. Nanotechnology 2011, 22, 155606. [Google Scholar] [CrossRef] [PubMed]
  36. daRosa, C.P.; Iglesia, E.; Maboudian, R. Dynamics of copper deposition onto silicon by galvanic displacement. J. Electrochem. Soc. 2008, 155, D244–D250. [Google Scholar] [CrossRef] [Green Version]
  37. Torcheux, L.; Mayeux, A.; Chemla, M. Electrochemical coupling effects on the corrosion of silicon samples in HF solutions. J. Electrochem. Soc. 1995, 142, 2037. [Google Scholar] [CrossRef]
  38. Jeske, M.; Schultze, J.W.; Thönissen, M.; Münder, H. Electrodeposition of metals into porous silicon. Thin Solid Film. 1995, 255, 63–66. [Google Scholar] [CrossRef]
  39. Xia, X.H.; Ashruf, C.M.A.; French, P.J.; Kelly, J.J. Galvanic cell formation in silicon/metal contacts: The effect on silicon surface morphology. Chem. Mater. 2000, 12, 1671–1678. [Google Scholar] [CrossRef]
  40. Kolasinski, K.W.; Barclay, W.B. The stoichiometry of electroless silicon etching in solutions of V2O5 and HF. Angew. Chem. Int. Ed. 2013, 52, 6731–6734. [Google Scholar] [CrossRef]
  41. Tsuboi, T.; Sakka, T.; Ogata, Y.H. Metal deposition into a porous silicon layer by immersion plating: Influence of halogen ions. J. Appl. Phys. 1998, 83, 4501–4506. [Google Scholar] [CrossRef]
  42. daRosa, C.P.; Maboudian, R.; Iglesia, E. Copper deposition onto silicon by galvanic displacement: Effect of silicon dissolution rate. J. Electrochem. Soc. 2008, 155, E70–E78. [Google Scholar] [CrossRef]
  43. Chartier, C.; Bastide, S.; Lévy-Clément, C. Metal-assisted chemical etching of silicon in HF–H2O2. Electrochim. Acta 2008, 53, 5509–5516. [Google Scholar] [CrossRef]
  44. Cerruti, M.; Doerk, G.; Hernandez, G.; Carraro, C.; Maboudian, R. Galvanic deposition of Pt clusters on silicon: Effect of HF concentration and application as catalyst for silicon nanowire growth. Langmuir 2010, 26, 432–437. [Google Scholar] [CrossRef] [PubMed]
  45. Ma, C.; Ye, W.; Shi, X.; Chang, Y.; Chen, Y.; Wang, C. Comparative study of electroless copper deposition based on the seed layers of Pd, PtPd and AuPd. Appl. Surf. Sci. 2009, 255, 3713–3718. [Google Scholar] [CrossRef]
  46. Li, X.; Xiao, Y.; Yan, C.; Zhou, K.; Schweizer, S.L.; Sprafke, A.; Lee, J.-H.; Wehrspohn, R.B. Influence of the mobility of Pt nanoparticles on the anisotropic etching properties of silicon. Solid State Lett. 2013, 2, 22–24. [Google Scholar] [CrossRef]
  47. Kuznetsov, G.V.; Skryshevsky, V.A.; Vdovenkova, T.A.; Tsyganova, A.I.; Gorostiza, P.; Sanz, F. Platinum electroless deposition on silicon from hydrogen fluoride solutions. Electrical properties. J. Electrochem. Soc. 2001, 148, C528–C532. [Google Scholar] [CrossRef]
  48. Geyer, N.; Wollschläger, N.; Fuhrmann, B.; Tonkikh, A.; Berger, A.; Werner, P.; Jungmann, M.; Krause-Rehberg, R.; Leipner, H.S. Influence of the doping level on the porosity of silicon nanowires prepared by metal-assisted chemical etching. Nanotechnology 2015, 26, 245301–245307. [Google Scholar] [CrossRef]
  49. Li, G.; Jiao, J.; Seraphin, S.; Raghavan, S.; Jeon, J.S. Masking effect of copper during anisotropic etching of silicon in buffered hydrofluoric acid solutions. J. Appl. Phys. 1999, 85, 1857–1863. [Google Scholar] [CrossRef]
  50. Lombardi, I.; Marchionna, S.; Zangari, G.; Pizzini, S. Effect of Pt particle size and distribution on photoelectrochemical hydrogen evolution by p-Si photocathodes. Langmuir 2007, 23, 12413–12420. [Google Scholar] [CrossRef]
  51. Cichoszewski, J.; Reuter, M.; Schwerdt, F.; Werner, J.H. Role of catalyst concentration on metal assisted chemical etching of silicon. Electrochim. Acta 2013, 109, 333–339. [Google Scholar] [CrossRef]
  52. Yae, S.; Abe, M.; Fukumuro, N.; Matsuda, H. Palladium enhanced etching of n-type silicon in hydrofluoric acid solution. Electrochemistry 2008, 76, 144–146. [Google Scholar] [CrossRef] [Green Version]
  53. Memming, R.; Schwandt, G. Anodic dissolution of silicon in hydrofluoric acid solutions. Surf. Sci. 1966, 4, 109–124. [Google Scholar] [CrossRef]
  54. Bandarenka, H.; Redko, S.; Nenzi, P.; Balucan, M.; Bondarenko, V. Optimization of chemical displacement deposition of copper on porous silicon. J. Nanosci. Nanotechnol. 2012, 12, 8725–8731. [Google Scholar] [CrossRef] [PubMed]
  55. Smith, R.L.; Collins, S.D. Porous silicon formation mechanisms. J. Appl. Phys. 1992, 71, R1. [Google Scholar] [CrossRef]
  56. Peng, K.Q.; Yan, Y.; Gao, S.; Zhu, J. Dendrite-assisted growth of silicon nanowires in electroless metal deposition. Adv. Funct. Mater. 2003, 13, 127–132. [Google Scholar] [CrossRef]
  57. Liu, F.M.; Green, M. Efficient SERS substrates made by electroless silver deposition into patterned silicon structures. J. Mater. Chem. 2004, 14, 1526–1532. [Google Scholar] [CrossRef]
  58. Knotter, M. Etching mechanism of vitreous silicon dioxide in HF-based solutions. J. Am. Chem. Soc. 2000, 122, 4345–4351. [Google Scholar] [CrossRef]
  59. Abramova, E.N.; Khort, A.M.; Tsygankov, V.N.; Yakovenko, A.G.; Shvets, V.I. The role of the etchant ion in the formation and growth of pores in silicon during its etching in hydrofluoric acid solutions. Dokl. Chem. 2016, 467, 61–63. [Google Scholar] [CrossRef]
  60. Beckmann, K.H. Investigation of the chemical properties of stain films on silicon by means of infrared spectroscopy. Surf. Sci. 1965, 3, 314–332. [Google Scholar] [CrossRef]
  61. Gossenberger, F.; Roman, T.; Forster-Tonigold, K.; Groß, A. Change of the work function of platinum electrodes induced by halide adsorption. Beilstein J. Nanotechnol. 2014, 5, 152–161. [Google Scholar] [CrossRef]
  62. Vanýsek, P. Electrochemical series. In CRC Handbook of Chemistry and Physics, 93th ed.; Haynes, W.M., Ed.; Taylor & Francis Ltd.: Milton Park, UK, 2012; pp. 2012–2013. [Google Scholar]
  63. Ginstrup, O. The redox system platinum(0)/platinum(II)/platinum(IV) with chloro and bromo ligands. Acta Chem. Scand. 1972, 26, 1527–1541. [Google Scholar] [CrossRef] [Green Version]
  64. Bromley, L.A. Thermodynamic properties of strong electrolytes in aqueous solutions. Aiche J. 1973, 19, 313–320. [Google Scholar] [CrossRef]
  65. Trasatti, S. Work function, electronegativity, and electrochemical behaviour of metals. II. Potentials of zero charge and “electrochemical” work functions. J. Electroanal. Chem. 1971, 33, 351–378. [Google Scholar] [CrossRef]
  66. Chelvayohann, M.; Mee, C.H.B. Work function measurements on (110), (100) and (111) surfaces of silver. J. Phys. C Solid State Phys. 1982, 15, 2305–2312. [Google Scholar] [CrossRef]
  67. Hölzl, J.; Schulte, F.K. Work functions of metals. In Solid Surface Physics; Höhler, G., Ed.; Springer: Berlin, Germany, 1979. [Google Scholar]
  68. Rivière, J.C. Work function: Measurements and results. In Solid Surface Science; Green, M., Ed.; Decker: New York, NY, USA, 1969. [Google Scholar]
  69. Michaelson, H.B. The work function of the elements and its periodicity. J. Appl. Phys. 1977, 48, 4729. [Google Scholar] [CrossRef] [Green Version]
  70. Bertocchi, M.; Amato, M.; Marri, I.; Ossicini, S. Tuning the work function of Si (100) surface by halogen absorption: A DFT study. Phys. Status Solidi C 2017, 14, 1700193. [Google Scholar] [CrossRef]
  71. Zvyagin, I.P.; Kurova, I.A.; Nal’gieva, M.A.; Ormont, N.N. Hopping ε2 conductivity of boron-doped a-Si: H films annealed in hydrogen at high temperature. Semiconductors 2006, 40, 108–112. [Google Scholar] [CrossRef]
  72. Novikov, A. Experimental measurement of work function in doped silicon surfaces. Solid State Electron. 2010, 54, 8–13. [Google Scholar] [CrossRef]
  73. Fujimura, N.; Ohta, A.; Makihara, K.; Miyazaki, S. Evaluation of valence band top and electron affinity of SiO2 and Si-based semiconductors using X-ray photoelectron spectroscopy. Jpn. J. Appl. Phys. 2016, 55, 08PC06. [Google Scholar] [CrossRef]
  74. Reed, M.L.; Fedder, G.K. Micro mechanical systems 2—Photolithographic microfabrication. In Handbook of Sensors and Actuators; Fukuda, T., Menz, W., Eds.; Elsevier Science: Amsterdam, The Netherlands, 1998. [Google Scholar]
  75. Messnaoui, B. Representation of excess properties and liquid composition of aqueous solutions of the HF+water system. J. Solut. Chem. 2008, 37, 715–726. [Google Scholar] [CrossRef]
  76. Hunt, H.; Chittum, J.F.; Ritchev, H.W. Overvoltage. In Seventy-Third Generation Meeting; The Electrochemical Society: Savannah, GA, USA, 1938; pp. 299–314. [Google Scholar]
  77. Buxton, G.V.; Greenstock, C.L.; Helman, W.P.; Ross, A.B. Critical review of rate constants for reactions of hydrated electrons, hydrogen atoms and hydroxyl radicals (∙OH/O) in aqueous solution. J. Phys. Chem. Ref. Data 1988, 17, 513. [Google Scholar] [CrossRef] [Green Version]
  78. Zona, R.; Solar, S.; Getoff, N.; Sehested, K.; Holcman, J. Reactivity of H atoms and hydrated electrons with chlorobenzoic acids. Radiat. Phys. Chem. 2008, 77, 162–168. [Google Scholar] [CrossRef]
  79. Atkins, P.W.; Paula, J. Physikalische Chemie; WILEY-VCH Verlag GmbH & Co. KGaA: Weinheim, Germany, 2006. [Google Scholar]
  80. Giguère, P.A.; Turrell, S. The nature of hydrofluoric acid. A spectroscopic study of the proton-transfer complex H3O+∙F. J. Am. Chem. Soc. 1980, 102, 5473–5477. [Google Scholar] [CrossRef]
  81. Es, F.; Kulakcı, M.; Turan, R. An alternative metal-assisted etching route for texturing silicon wafers for solar cell applications. IEEE J. Photovolt. 2016, 6, 440–446. [Google Scholar] [CrossRef]
  82. Harraz, F.A.; Sakka, T.; Ogata, Y.H. Immersion plating of nickel onto a porous silicon layer from fluoride solutions. Phys. Status Solidi A 2003, 197, 51–56. [Google Scholar] [CrossRef]
  83. Memming, R. Semiconductor Electrochemistry; WILEY-VCH Verlag GmbH & Co. KGaA: Weinheim, Germany, 2001. [Google Scholar]
  84. Papaconstantopoulos, D.A.; Economou, E.N. Calculation of the electronic properties of hydrogenated silicon. Phys. Rev. B 1981, 24, 7233–7246. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Dissolved (a) Ag, (b) Cu, (c) Au and (d) Pt species b Me (diss., t = i) and (ad) Si species b Si (diss., t = i) at t = i; the period of the metal deposition phases 1 (t = 0 until t1/2) and 2 (t1/2 until t95%); the metal deposition rates k1 and k2; the t95%-weighted deposition rate k95%; and the stoichiometric ratio of the change in molality of the dissolved metal ion Δb Me in relation to the change in molality of the dissolved silicon species Δb Si in the period of t = 0 until t95%b Meb Si (diss., t95%)).
Figure 1. Dissolved (a) Ag, (b) Cu, (c) Au and (d) Pt species b Me (diss., t = i) and (ad) Si species b Si (diss., t = i) at t = i; the period of the metal deposition phases 1 (t = 0 until t1/2) and 2 (t1/2 until t95%); the metal deposition rates k1 and k2; the t95%-weighted deposition rate k95%; and the stoichiometric ratio of the change in molality of the dissolved metal ion Δb Me in relation to the change in molality of the dissolved silicon species Δb Si in the period of t = 0 until t95%b Meb Si (diss., t95%)).
Nanomaterials 10 02545 g001
Figure 2. Metal deposition rates k95% (k95%,1 and k95%,2) of the (a) Ag, (b) Cu, (c) Au and (d) Pt depositions in relation to the initial metal ion molality (b Me (t = 0), Me = Ag, Cu, Au, Pt) or to the redox potential of the Mez+/Me half-cells E (Mez+/Me). Red dots = examples from Figure 1, grey dots = findings with delayed silicon dissolution.
Figure 2. Metal deposition rates k95% (k95%,1 and k95%,2) of the (a) Ag, (b) Cu, (c) Au and (d) Pt depositions in relation to the initial metal ion molality (b Me (t = 0), Me = Ag, Cu, Au, Pt) or to the redox potential of the Mez+/Me half-cells E (Mez+/Me). Red dots = examples from Figure 1, grey dots = findings with delayed silicon dissolution.
Nanomaterials 10 02545 g002
Figure 3. Work functions Φ of the metals (a) Ag, (b) Cu, (c) Au and (d) Pt and that of hydrogen-terminated Si in vacuum [65,66,67,68,69,70]; the difference in work functions ΔE; the Fermi energy level EF and that of the valence and conduction band of hydrogen-terminated silicon [71] (Ev (Si–H) or Ec (Si–H)) and of bulk silicon [72,73] (EV (Sibulk) or Ec (Sibulk)); the energy levels of the valence bands at the metal–silicon contact (Ev (Sibulk/Me-Si) or Ev (Si–H/Me-Si)); the Schottky barrier Esb; the redox levels of the Mez+/Me half-cells E (Mez+/Me); and the equilibrium potentials of the 2H3O+/H2 half-cells at the metal surface Equation Pot. (2H3O+/H2).
Figure 3. Work functions Φ of the metals (a) Ag, (b) Cu, (c) Au and (d) Pt and that of hydrogen-terminated Si in vacuum [65,66,67,68,69,70]; the difference in work functions ΔE; the Fermi energy level EF and that of the valence and conduction band of hydrogen-terminated silicon [71] (Ev (Si–H) or Ec (Si–H)) and of bulk silicon [72,73] (EV (Sibulk) or Ec (Sibulk)); the energy levels of the valence bands at the metal–silicon contact (Ev (Sibulk/Me-Si) or Ev (Si–H/Me-Si)); the Schottky barrier Esb; the redox levels of the Mez+/Me half-cells E (Mez+/Me); and the equilibrium potentials of the 2H3O+/H2 half-cells at the metal surface Equation Pot. (2H3O+/H2).
Nanomaterials 10 02545 g003
Figure 4. Left: Sum of the initial activities of the anionic HF dissociation products F, HF2 and H2F3 (a (F + HF2 + H2F3) (t = 0)) and molality of the dissolved silicon at t95% (b Si (diss., t95%)) in relation to the initial metal ion molality: (a) b Ag (diss., t = 0) or (b) b Cu (diss., t = 0). Right: Quotient of a (F + HF2 + H2F3) (t = 0) and b Si (diss., t95%) in relation to b Ag (diss., t = 0) or b Cu (diss., t = 0). (Red dots = examples from Figure 1, grey dots = findings with delayed silicon dissolution).
Figure 4. Left: Sum of the initial activities of the anionic HF dissociation products F, HF2 and H2F3 (a (F + HF2 + H2F3) (t = 0)) and molality of the dissolved silicon at t95% (b Si (diss., t95%)) in relation to the initial metal ion molality: (a) b Ag (diss., t = 0) or (b) b Cu (diss., t = 0). Right: Quotient of a (F + HF2 + H2F3) (t = 0) and b Si (diss., t95%) in relation to b Ag (diss., t = 0) or b Cu (diss., t = 0). (Red dots = examples from Figure 1, grey dots = findings with delayed silicon dissolution).
Nanomaterials 10 02545 g004
Figure 5. Stoichiometric ratios between the metal deposition Δb Me ((a) Ag, (b) Cu, (c) Au and (d) Pt) and silicon dissolution Δb Si as a function of time tb Meb Si (diss., t = i)) based on the experiments of Figure 1 with an indication of the t95%-weighted stoichiometric ratios (red dots) and the corresponding minima and maxima (bars) denoted as (Δb Meb Si (diss., t95%)).
Figure 5. Stoichiometric ratios between the metal deposition Δb Me ((a) Ag, (b) Cu, (c) Au and (d) Pt) and silicon dissolution Δb Si as a function of time tb Meb Si (diss., t = i)) based on the experiments of Figure 1 with an indication of the t95%-weighted stoichiometric ratios (red dots) and the corresponding minima and maxima (bars) denoted as (Δb Meb Si (diss., t95%)).
Nanomaterials 10 02545 g005
Figure 6. Stoichiometric ratios between the metal deposition Δb Me ((a) Ag, (b) Cu, (c) Au and (d) Pt) and silicon dissolution Δb Si until t95%b Meb Si (diss., t95%)), in relation to the initial metal ion molality (b Me (diss., t = 0)), are indicated as points for the t95%-weighted values and as bars for the corresponding minima and maxima with delimitation by sigmoid functions f1, f2 and f3 and classification of the underlying processes in sections I, II, III and IV (red dots = examples from Figure 1, grey dots = findings with limited silicon dissolution).
Figure 6. Stoichiometric ratios between the metal deposition Δb Me ((a) Ag, (b) Cu, (c) Au and (d) Pt) and silicon dissolution Δb Si until t95%b Meb Si (diss., t95%)), in relation to the initial metal ion molality (b Me (diss., t = 0)), are indicated as points for the t95%-weighted values and as bars for the corresponding minima and maxima with delimitation by sigmoid functions f1, f2 and f3 and classification of the underlying processes in sections I, II, III and IV (red dots = examples from Figure 1, grey dots = findings with limited silicon dissolution).
Nanomaterials 10 02545 g006
Figure 7. Equilibrium potentials of the 2H3O+/H2 half-cells at the (a) Ag, (b) Cu, (c) Au or (d) Pt surfaces (bright metal state/colloidal metal solution [76]) in relation to the initial metal ion molalities (b Me (diss., t = 0)), the level of the Fermi energy EF and levels of the valence bands of the bulk silicon (Ev (Sibulk/Me-Si) and that of the hydrogen-terminated silicon (Ev (Si–H/Me-Si)) at the metal–silicon contact. (Hollow symbols = oxonium ion reduction to molecular hydrogen possible, full symbols = no reaction).
Figure 7. Equilibrium potentials of the 2H3O+/H2 half-cells at the (a) Ag, (b) Cu, (c) Au or (d) Pt surfaces (bright metal state/colloidal metal solution [76]) in relation to the initial metal ion molalities (b Me (diss., t = 0)), the level of the Fermi energy EF and levels of the valence bands of the bulk silicon (Ev (Sibulk/Me-Si) and that of the hydrogen-terminated silicon (Ev (Si–H/Me-Si)) at the metal–silicon contact. (Hollow symbols = oxonium ion reduction to molecular hydrogen possible, full symbols = no reaction).
Nanomaterials 10 02545 g007
Figure 8. Molality of dissolved silicon in (a) gold deposition after t = 3600 s (b Si (diss., t = 3600 s)) and (b) platinum deposition after t = 14,400 s (b Si (diss., t = 14,400 s) in relation to the initial metal ion molality (b Au (diss., t = 0) and b Pt (diss., t = 0), respectively) with and without the presence of para-chlorobenzoic acid (pCBA).
Figure 8. Molality of dissolved silicon in (a) gold deposition after t = 3600 s (b Si (diss., t = 3600 s)) and (b) platinum deposition after t = 14,400 s (b Si (diss., t = 14,400 s) in relation to the initial metal ion molality (b Au (diss., t = 0) and b Pt (diss., t = 0), respectively) with and without the presence of para-chlorobenzoic acid (pCBA).
Nanomaterials 10 02545 g008
Figure 9. Proposed scheme of the reaction of the ion pair H3O+∙F with silicon.
Figure 9. Proposed scheme of the reaction of the ion pair H3O+∙F with silicon.
Nanomaterials 10 02545 g009
Figure 10. Reaction schemes of (a,b) silver, (c,d) copper, (e,f) gold and (g,h) platinum deposition on hydrogen-terminated silicon within section III.
Figure 10. Reaction schemes of (a,b) silver, (c,d) copper, (e,f) gold and (g,h) platinum deposition on hydrogen-terminated silicon within section III.
Nanomaterials 10 02545 g010
Figure 11. Reaction schemes of (a) silver, (b) copper, (c) gold and (d) platinum deposition on hydrogen-terminated silicon within section IV.
Figure 11. Reaction schemes of (a) silver, (b) copper, (c) gold and (d) platinum deposition on hydrogen-terminated silicon within section IV.
Nanomaterials 10 02545 g011
Table 1. Parameterization of the sigmoidal functions f1, f2 and f3 according to Equations (10) and (11).
Table 1. Parameterization of the sigmoidal functions f1, f2 and f3 according to Equations (10) and (11).
ParameterUnitAgCuAuPt
f1b Me (diss., t = 0) valid rangemol7.1 × 10−75.5 × 10−32.1 × 10−58.2 × 10−33.3 × 10−74.2 × 10−76.0 × 10−5
kg5.5 × 10−31.3 × 10−18.2 × 10−34.1 × 10−19.3 × 10−25.1 × 10−59.7 × 10−3
(Δb Meb Si)maxmol∙kg−14:15:12:13:14:31:16:5
mol∙kg−1
xmol0.993.0 × 10−53.3 × 10−44.5 × 10−51.5 × 10−34.1 × 10−55.1 × 10−53.5 × 10−5
kg0.99
f2b Me (diss., t = 0) valid rangemol7.1 × 10−7–1.3 × 10−12.1 × 10−5–2.3 × 10−33.3 × 10−74.2 × 10−7–6.5 × 10−4
kg2.4 × 10−3
(Δb Meb Si)maxmol∙kg−110:32:14:36:5
mol∙kg−1
xmol0.992.5 × 10−44.0 × 10−43.7 × 10−42.1 × 10−4
kg0.99
f3b Me (diss., t = 0) valid rangemol---2.3 × 10−3–4.1 × 10−12.4 × 10−36.5 × 10−4–9.7 × 10−3
kg9.3 × 10−2
(Δb Meb Si)minmol∙kg−1---1:21:31:2
mol∙kg−1
xmol0.99---3.5 × 10−39.0 × 10−31.2 × 10−3
kg0.99
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Schönekerl, S.; Acker, J. The Kinetics and Stoichiometry of Metal Cation Reduction on Multi-Crystalline Silicon in a Dilute Hydrofluoric Acid Matrix. Nanomaterials 2020, 10, 2545. https://doi.org/10.3390/nano10122545

AMA Style

Schönekerl S, Acker J. The Kinetics and Stoichiometry of Metal Cation Reduction on Multi-Crystalline Silicon in a Dilute Hydrofluoric Acid Matrix. Nanomaterials. 2020; 10(12):2545. https://doi.org/10.3390/nano10122545

Chicago/Turabian Style

Schönekerl, Stefan, and Jörg Acker. 2020. "The Kinetics and Stoichiometry of Metal Cation Reduction on Multi-Crystalline Silicon in a Dilute Hydrofluoric Acid Matrix" Nanomaterials 10, no. 12: 2545. https://doi.org/10.3390/nano10122545

APA Style

Schönekerl, S., & Acker, J. (2020). The Kinetics and Stoichiometry of Metal Cation Reduction on Multi-Crystalline Silicon in a Dilute Hydrofluoric Acid Matrix. Nanomaterials, 10(12), 2545. https://doi.org/10.3390/nano10122545

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop