Yufei
Wang‡
a,
Roberto
Garcia-Carrillo‡
a and
Hang
Ren
*abc
aDepartment of Chemistry, The University of Texas at Austin, Austin TX 78712, USA. E-mail: hren@utexas.edu
bCenter for Electrochemistry, The University of Texas at Austin, Austin TX 78712, USA
cTexas Materials Institute, The University of Texas at Austin, Austin TX 78712, USA
First published on 16th December 2024
Electrochemical metal dissolution reaction is a fundamental process in various critical technologies, including metal anode batteries and nanofabrication. However, experimentally revealing the kinetics and dynamics of active sites of metal dissolution reactions is challenging. Herein, we investigate metal dissolution on near-perfect single-crystal surfaces of Ag within regions of a few hundred nanometers isolated by scanning electrochemical cell microscopy (SECCM). Potential oscillation is observed under constant current conditions for dissolution. The one-to-one correspondence between the dissolution charge and the geometry of the dissolution pit from colocalized imaging allows ambiguous correlation, which suggests that each oscillation cycle corresponds to the dissolution of one atomic layer. The oscillation behavior is further explained in a kinetic model, which reveals that the oscillation comes from the dynamic evolution of the number of different active sites as the dissolution progresses on each atomic layer. In addition to the fundamental interest, the ability to observe layer-by-layer dissolution in electrochemical measurements suggests a potential pathway for developing electrochemical atomic layer etching for fabricating structures and devices with atomic precision.
Electrochemical metal dissolution reaction is one of the elementary reactions in many important technical processes, such as corrosion,1–3 sensors,4–6 material synthesis,7,8 and battery technology.9–11 Kinetics and rate-determining steps of metal dissolution reaction can be studied electrochemically.12–15 However, in conventional electrochemical experiments, the electrode surface is not an ideal homogeneous surface. It inevitably exhibits diverse surface structures, such as grains with different crystal orientations and defects, including grain boundaries, dislocations, and point defects. These surface characteristics significantly affect thermodynamics, kinetics, and mechanism of metal dissolution reactions. Building the microscopic structure–activity relationship of metal dissolution reaction at different surface structures and sites is crucial in understanding the dissolution behavior at complex metal surfaces.
Silver is a model system for studying metal dissolution reactions. Despic and Bockris examined the kinetics of metal dissolution and highlighted the role of surface diffusion and adatom concentrations,16 which is also confirmed by Gerischer and coworkers.17,18 Hackerman studied this reaction on a polycrystalline silver electrode using electrochemical impedance spectroscopy (EIS) and concluded the reaction is controlled by adatom diffusion.19 In all previous examples, electrochemical kinetics were measured using ensemble approaches on surfaces with unknown distributions of surface structures, including different local densities of dislocations. Because the measured activity is an ensemble average over the entire electrode surface, it is difficult to develop a site-specific mechanism that can be directly verified by experiments.
To better define the relationship between surface structure and the activity of metal dissolution, a local electrochemical method is desirable because it allows the isolation of specific surface structures. One emerging technique that allows such isolation is SECCM, which was developed by the Unwin group.20 Several studies have investigated metal dissolution and its relationship with crystal orientations.21–25 The influence of different electrolytes on metal dissolution kinetics and mechanisms has also been explored.26–28 Mauzeroll and coworkers applied SECCM to study the electrochemical metal dissolution and passivation on different alloys and interpreted the results using a high-field model.26–28 Our group have previously shown that the heterogeneity of metal dissolution at different surface structures can be elucidated using SECCM.29 Combined with colocalized structural characterization using SEM, EBSD, and FIB-assisted TEM, we have explained the origin of heterogeneity in the metal dissolution reactions on different grain boundaries.30 In addition, the effect of local oxide layer thickness on the local breakdown potential was also obtained via combined SECCM and TOF-SIMS. At low-defect crystal planes, however, the atomistic mechanisms of metal dissolution remain elusive.31
In this report, we demonstrated the SECCM approach to study the kinetics and mechanism of Ag dissolution at the low-defect single crystal (100) and (111) surfaces. The low-defect single crystal surface was prepared using Budewski's electrodeposition method.32 In galvanostatic experiments, potential oscillation is observed during the anodic dissolution. A model based on site-specific kinetics of dissolution on kink vs. edge sites was developed, which explains how the potential oscillation stems from layer-by-layer dissolution, with each oscillation cycle corresponding to monolayer dissolution. The preparation of high-quality electrode surface coupled with shrinking of electrode area under investigation by SECCM provide a unique design that allows us to link new experimental data with atomistic dissolution process, allowing new insight into metal dissolution reactions at low-defect surfaces.
Ag microparticles were deposited using constant potential in 6.0 M AgNO3 (99.9+% metal basis, Thermo Scientific) and 100 mM HNO3 (Trace Grade, Spectrum Chemical). The electrodeposition procedures are described in ESI Section 1.†
The SECCM was performed in a home-built scanning probe system which was covered by a thermally insulated Faraday cage, similar to our previous report.29,30 The nanopipette was fabricated using a laser puller (Sutter, P-2000) and characterized using a scanning electron microscope. A silver wire (0.127 mm dia, 99.9% metals basis, Thermo Scientific) was used as a reference electrode for all SECCM experiments and microparticle deposition. The constant current Ag dissolution experiment is described in the ESI Section 2.† All SEM images were obtained on a FEI Quanta 650 ESEM system with the 10 keV accelerate voltage and 10 mm working distance.
A zero current was applied between the substrate and the Ag quasi-reference counter electrode (QRCE) inside the nanopipette when the nanopipette approached the microparticle. When the droplet was in contact with the substrate, the circuit was closed, resulting in an abrupt change in potential. This served as feedback to pause the pipette and prevent it from crashing onto the surface. A constant anodic current was then applied to induce Ag dissolution while a potential transient was recorded. After an electrochemical measurement of dissolution at one location was complete, the nanopipette was retracted and moved to a new location for a subsequent experiment.
When dissolution occurs at different crystallographic facets, the geometries of the dissolution pits vary. As shown in the SEM images in Fig. 2, the pits exhibit projected shapes of near-squares on (100) facets (Fig. 2a), and triangles on (111) facets (Fig. 2d). The 3-fold and 4-fold symmetries of the (111) and (100) facets are reflected on these projected geometries of the dissolution pits. The potential transients at 2 pA applied current in Fig. 2c and f show that at t = 0, the potential quickly drops initially, then oscillates with a peak-to-peak amplitude of a few mV. The average period of oscillation is on the orders of 100 to 200 ms, and varies slightly from cycle to cycle. Conveniently, this potential is measured vs. the equilibrium potential of Ag/Ag+ due to the use of Ag QRCE, and therefore represents the overpotential for Ag dissolution.
We hypothesize that the dissolution in this isolated low-defect region operates in a layer-by-layer dissolution mechanism, with each oscillation cycle corresponding to the dissolution of one atomic layer of Ag. This layer-by-layer mechanism would have several predictions: first, dissolution predominantly occurs along a direction parallel to a certain crystal plane. Because the surface atoms on the (111) plane have the highest coordination number, the preferred dissolution direction should be parallel to (111). Second, the charge associated with each oscillation corresponds to a monolayer of Ag. Third, the frequency of oscillation increases with the applied current. This is because increasing current means a faster dissolution rate, therefore reducing the time to dissolve a specific monolayer.
The first prediction on the geometry is consistent with the triangle cross-section of the dissolution pit left on Ag(100) and Ag(111). On (100) surfaces, the geometry of the dissolution pits is confirmed as square-pyramidal by AFM (Fig. 2b). On (111) surfaces, the dissolution pit assumes a triangular prism shape with a flat bottom face parallel to the (111) surface of the microparticle (Fig. 2e). Both geometries are consistent with the preferred dissolution along the (111) plane. The deviation of the cross-section from a perfect square in some samples is attributed to dissolution starting from two positions on the surface and fusing together during the simultaneous dissolution and the tilting of our probe with respect to the normal of the sample surface. In addition, an oscillation in potential was observed during the dissolution on both Ag (100) (Fig. 2c) and Ag (111) (Fig. 2f). More examples of potential transients during dissolution of (111) and (100) are shown in ESI Section 4† and at different times with 10 pA of applied current in ESI Section 5.†
To investigate the second prediction, the charge from dissolution at different applied currents is compared with that calculated from AFM topography as described in ESI Section 6.† The charge agrees with that calculated using the volume of the dissolution pit (Fig. S6†). Note that the colocalization between the AFM and the SECCM allows us to associate each local E–t curve for dissolution to the exact dissolution pit under the AFM, which greatly reduces the error in the analysis from the variation at different locations. Furthermore, the charge associated with a single potential oscillation cycle is calculated from the average oscillation cycle duration and the applied current. This charge obtained from electrochemical measurement is compared with the expected charge of the monolayer based on the geometry of the dissolution pit. As shown in ESI Section 7,† the monolayer charge for dissolution on the Ag(100) surface is calculated using the area of the four (111) faces of the square pyramidal dissolution pit, which is obtained by the AFM after the SECCM dissolution experiment (Fig. S6†). The surface area of the dissolution pit on Ag (100) continues to grow vs. charge, so the “monolayer” charge also increases as more Ag layers are dissolved (Fig. S7†). We observed that the monolayer charge measured electrochemically and derived from the AFM topography both increase with the applied current (the dissolution time is fixed), as shown in Fig. S8.† When the dissolution current is <6 pA, the expected charge calculated from AFM topography agrees with the electrochemically measured charge, supporting the mechanism of layer-by-layer dissolution. At higher applied currents (>6 pA), the agreement becomes worse, likely due to the simultaneous stripping of multiple monolayers or overlapping of the potential transients. We speculate that this effect is ultimately responsible for the absence of visible oscillations in bulk electrochemical measurements on polycrystalline Ag.15,16,34
To test the third prediction of the layer-by-layer mechanism, we measured the potential transient as a function of applied current. The potential oscillations become denser as the applied current increases on both Ag (100) and Ag (111), as shown in Fig. 3a and b, respectively. This trend qualitatively agrees with the prediction. Fast Fourier Transform (FFT) analysis of the E–t curves during dissolution shows no prevalent frequency, indicating that the potential oscillations are not periodic (Fig. S9†).
The duration of each oscillation (Δt) was further analyzed by obtaining its distribution following a peak-finding procedure described in ESI Section 8.† As shown in Fig. 4a and b, the distribution of Δt shifts to a shorter time as the applied current increases on both Ag (100) and Ag (111). On Ag (100), Δt decreases from 150 ± 6 to 49 ± 1 ms (mean ± standard deviation) when the current increases from 1 to 10 pA (Fig. 4c), which qualitatively consistent with the third prediction. We note that an increase of current by 10× does not lead to a 10× decrease in Δt. This is attributed to two reasons: (1) the area of active dissolving (111) planes is larger at a higher current as expected for the expansion of the pyramidal pit, which can also be observed in the SEM images in Fig. 3. The larger area of active (111) planes at a higher current means more charge to completely dissolve a monolayer, and therefore, the frequency will decrease. (2) The potential oscillation overlaps at a higher applied current, leading to an underestimation of Δt. On Ag (111) planes, the average oscillation duration decreases from 360 ± 30 to 30 ± 10 ms when the current decreases from 1 to 10 pA (Fig. 4d). We did not attempt to quantitatively compare the difference of Δt between Ag (100) and Ag (111) because different nanopipettes were used, and the variations in the pipette/droplet size lead to changes in the absolute charge associated with (111) monolayers exposed at the electrode/electrolyte interface. We also note that the size of the pits can be larger than the pipette opening, which is attributed to droplet wetting on Ag, and the variation of the pipette size.
We also showed that similar potential oscillation behavior is observed under N2 and O2 environments (ESI Section 9†), suggesting O2 is not involved in the oscillation. To further elucidate the dynamics of Ag dissolution on a low-defect surface, we simulate the E–t traces using a model that considers the site-specific kinetics and double-layer charging during galvanostatic electro-dissolution. The most reactive surface sites are the ones with fewer neighboring interactions. Therefore, kink sites (6 neighboring atoms) should dissolve faster than edge sites (7 neighboring atoms). The dynamics of formation and annihilation between kink sites and edge sites are shown in Fig. 5b–d, which are considered in the simulation. Typically, edge dissolution creates two new kink sites (Fig. 5b), while kink dissolution does not change the number of kink sites (Fig. 5c) unless two neighboring kink sites are annihilated (Fig. 5d). Other sites, such as adatoms or atoms within a perfect layer (coordination number of 9), are neglected for simplicity. The geometry of the dissolution pits and the number of atoms in dissolution in each layer are also considered in the model as illustrated in Fig. S14.† A description of the mathematical model is presented in ESI Section 11.†
Simulated potential transients at different applied currents are presented in Fig. 5e, which shows the oscillation behavior. From the model, the oscillation in the current comes from the dynamic oscillation in the number of kink and edge sites as the dissolution progresses in each layer (vide infra). The average duration of the peak increases as the current increases, consistent with the experiments. On Ag (100), the frequency of the potential oscillation decreases as more layers are dissolved as required by the expansion of the pyramidal dissolution pit (Fig. S7†). The frequency stabilizes once the dissolution pits become sufficiently large. This trend in the frequency of oscillation is also observed experimentally on Ag (100) at a low dissolution current (e.g., 1 pA in Fig. 3a). At higher applied currents, the pits quickly become large enough that the oscillation frequency appears constant.
We observe several differences between the experiment and simulation. First, the duration of each oscillation is more erratic in the experiment compared to the simulation. This is likely due to the neglect of dissolution at sites other than kink or edge. Second, the potential overall increases with the applied current in the simulation, and significant heterogeneity in the current oscillation is observed at higher current (see also Fig. S4†). This discrepancy is attributed to the simplified assumptions in our model: the simulation considers sublayer dissolution only occurring after top layer is completely dissolved. Simultaneous dissolution of both top layer and sub-layers becomes important at higher applied current, which causes overlap of current contributions from multiple layers and deviation from our model prediction. Some shape variations in the potential transients can also be attributed to an induction mechanism, where the potential rises slowly until the activation barrier for the next layer is overcomed, followed by a sharp decrease. The details of this possible mechanism are discussed in ESI Section S12.† A more sophisticated model considering this phenomenon is beyond the scope of this current paper and will be explored in future work.
Combining simulation and experimental results, the potential oscillation can be understood through the following mechanism. As a monolayer starts dissolving, the number of kink atoms (i.e., active sites) increases through edge dissolution, which lowers the potential required to sustain the same current. This is because the kinetics (e.g., exchange current) is proportional to the number of atoms at the active sites. As the dissolution continues the monolayer, kink annihilation decreases the number of active sites, causing the overpotential to increase. The potential increases until a layer of atoms is completely dissolved and then decreases again as the next layer begins to dissolve.
The insights into the metal dissolution reaction have broad applicability. The selective etching demonstrated in this work is directly relevant to nanofabrication techniques, including nanostructure sculpting for tuning plasmonic resonance,34,35 and controlled layer dissolution for lithography fabrication.36 Extending of the kinetic and mechanistic studies to 2D materials, such as transition metal dichalcogenides (e.g. MoS2), can establish the foundational knowledge required to precisely control 2D material layers through etching, which hold promise for next-generation electronics.37,38
This study can be extended to polycrystalline materials, where techniques such as Electron Backscatter Diffraction (EBSD) can map crystallographic orientations and correlate them with local dissolution behavior to reveal how different facets and grain boundaries influence etch rates. The site-specific insights from local electrochemistry can inform processes like atomic layer etching (ALE) in polycrystalline films, helping to address uniformity challenges in nanofabrication and semiconductor applications.
Footnotes |
† Electronic supplementary information (ESI) available: Experimental and simulation details; additional electrochemistry data; and additional simulation data. See DOI: https://doi.org/10.1039/d4sc05954a |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2025 |