Xuefeng Ren*a and
Maozhong An*b
aSchool of Food and Environment, Dalian University of Technology, Panjin, 124221, China. E-mail: rxfhit@163.com; renxuefeng@dlut.edu.cn
bState Key Laboratory of Urban Water Resource and Environment, School of Chemistry and Chemical Engineering, Harbin Institute of Technology, Harbin, 150001, China. E-mail: mzan@hit.edu.cn
First published on 12th January 2018
Quantum chemical calculations based on density functional theory (DFT) were employed to determine an appropriate gold source for gold electroplating and to ascertain the stable structures of gold-complexes in cyanide-free electrolyte. Based on the charge distribution of 5,5-dimethylhydantoin (DMH) and the bonding energy of gold complexes, Au3+ is the appropriate gold source for DMH-based gold electroplating electrolyte to get greater cathodic polarization and [Au(DMH)4]− with 2N(4)–Au coordination structure is the most stable form of gold ion in the electrolyte. The influence of DMH, used as the complexing agent, on electrochemical behaviors was investigated using cathodic polarization, cyclic voltammetry, and chronoamperometry measurements. With DMH as the complexing agent, the cathodic polarization of gold electrodeposition was significantly enhanced. DMH concentration was determined as 0.30 mol L−1 based on the investigation of the influence of the DMH concentration on cathodic polarization and gold electrodeposit micromorphology. The kinetic features based on cyclic voltammogram measurements revealed that the electrodeposition was an irreversible process under diffusion control with 0.30 mol L−1 DMH as the complexing agent. The ion and electron transfers were obviously inhibited by DMH. The gold electrodeposition process displayed progressive nucleation according to the Scharifker and Hills nucleation model with various applied potentials. The growth rate of the crystal nucleus was reduced by DMH and promoted by a negative shift of Eap.
A possible alternate to cyanide is DMH, a low cost and commercially available derivative of hydantoin with good solubility and stability in aqueous solutions in a large temperature range. DMH is a promising environmental friendly complexing agent, which can coordinate with many metal ions.20–22 The DMH-based electrolyte is simple, stable, environmentally friendly, and compatible with various chemicals and substrates. One Au3+ ion and four DMH− ions can form a coordination complex [Au(DMH)4]− in the electrolyte.23
It is important to get an insight into the bonding interaction of gold ions and DMH− ions in the electrolyte as well as study the mechanisms of action of DMH on the gold electroplating process. These insights will provide important guidance to perfect the DMH-based gold electroplating electrolyte. Quantum chemical calculation and molecular modeling are fast emerging areas of research for the modeling of small chemical and biological systems in order to understand and predict their behaviors at the molecular level.24–28 Quantum chemical calculations based on DFT have become an effective technique for the calculation of many natural systems in pharmacology, chemistry, and biology.29–32 The basic foundation of DFT states that the external potential is a functional of the group-state density, and the density (an observable in 3D space) is used to describe the complicated physics behind the interactions between electrons.33–36 DFT has made an unparalleled impact on the application of quantum mechanics to interesting and challenging problems in chemistry for the past several decades.37–39 In coordinate systems, bonding interactions between metal ions and complexing agents can be predicted and investigated by quantum chemical calculations.40–43
In this work, the charge distributions of DMH as well as the structures of gold(I) and gold(III) complexes were studied by quantum chemical calculations to determine the appropriate gold source for the gold electroplating. The role of DMH used as the complexing agent in cyanide-free gold electroplating electrolyte was investigated by cathodic polarization, cyclic voltammetry, and chronoamperometry. The morphology of the gold electrodeposits was evaluated by scanning electron microscope (SEM).
All quantum chemical calculations were carried out using the B3LYP functional method. In the geometry optimization and single point energy calculation, the 6-311++G** basis set was used for the C, N, O, and H atoms of the studied systems, except Au ions, for which we used the LANL2DZ ECP basis set. All calculations in these systems under investigation were performed using the Gaussian 09 program package at 318 K with water as the solvent in the IEFPCM theoretical model. The bonding energies Ebonding between Au+ or Au3+ and DMH− were calculated using eqn (1) and (2), where EAu+, EAu3+, EDMH, EH2O, and EH3O+ were the energy of Au+, Au3+, DMH, H2O, and H3O+, respectively, and Ecomplexes was the energy of the gold complex.
Ebonding = EAu+ + 2EDMH + 2EH2O − Ecomplexes − 2EH3O+ | (1) |
Ebonding = EAu3+ + 4EDMH + 4EH2O − Ecomplexes − 4EH3O+ | (2) |
All electrochemical measurements were performed in a typical three-electrode cell using an electrochemical workstation at 318 K. Platinum foil with 1 cm2 surface area and a saturated calomel electrode (SCE) were employed as the counter electrode (CE) and the reference electrode (RE), respectively. A platinum rotating disk electrode (Pt-RDE) with a diameter of 5 mm (with 0.19625 cm2 surface area) was employed as the working electrode (WE) for the cathodic polarization measurements, and a glassy carbon electrode (GCE) with a diameter of 3 mm (with 0.07065 cm2 surface area) was used for cyclic voltammograms and chronoamperometry measurements.
Gold electrodeposits were deposited on a copper sheet coated with a nickel deposit. Gold electrodeposits were all prepared under galvanostatic conditions (8 mA cm−2) with mild agitation at 318 K for 5 min, either in the absence or presence of different DMH concentrations. The surface morphologies and microstructure of the gold electrodeposits were studied by field emission scanning electron microscopy (FE-SEM, Hitachi SU8000).
In this work, the charge distributions of DMH as well as the structures of gold(I) and gold(III) complexes were studied by quantum chemical calculations to determine the appropriate gold source for the gold electroplating. The molecular structure and charge distribution of DMH (a heterocyclic structure organic molecule) with atom symbols is displayed in Fig. 1.
As shown in Fig. 1, 2N (−0.307) and 4N (−0.218) of DMH were the most likely atoms to form coordinate bonds with gold ions. Three gold(I) complexes and 12 gold(III) complexes were possible structures of gold complexes in the electrolyte according to the charge distribution of DMH. The stability of the gold(I) and gold(III) complexes in the gold electroplating electrolyte was studied by quantum chemical calculations. The bonding energies of different complexes under study were calculated at the B3LYP/GENECP level. Fig. 2 and 3 present the optimized structures of the gold(I) and gold(III) complexes with DMH, respectively, based on their electronic properties.
Fig. 2 Optimized structures of the gold(I) complexes with DMH (the yellow balls represent Au; the symbols of other atoms are shown in Fig. 1), (a) 2N(2)–Au, (b) 2N(1)–Au–4N(1), (c) 4N(2)–Au. |
Fig. 3 Optimized structures of the gold(III) complexes with DMH (the yellow ball represents Au, the symbols of the other atoms are shown in Fig. 1), (a) 2N(4)–Au, (b) 4N(4)–Au, (c) 2N(3)–Au–4N(1), (d) 2N(3)–Au–4N(1), (e) 2N(3)–Au–4N(1), (f) 2N(3)–Au–4N(1), (g) 2N(2)–Au–4N(2), (h) 2N(2)–Au–4N(2), (i) 2N(1)–Au–4N(3), (j) 2N(1)–Au–4N(3), (k) 2N(1)–Au–4N(3), (l) 2N(1)–Au–4N(3). |
As shown in Fig. 2 and 3, these three gold(I) complexes and 12 gold(III) complexes were stable structures in the DMH-based electrolyte. The key geometric parameters and bonding energies of these studied complexes are summarized in Tables 1–3.
EAu3+ | EAu+ | EDMH | EH2O | EH3O+ |
---|---|---|---|---|
−352640.16 | −355088.37 | −1195927.71 | −200762.74 | −201783.32 |
Number | Structure | Ecomplexes | Ebonding |
---|---|---|---|
a | 2N(2)–Au | −2744863.54 | −39.08 |
b | 2N(1)–Au–4N(1) | −2744847.35 | −55.27 |
c | 4N(2)–Au | −2744830.70 | −71.92 |
Number | Structures | Ecomplexes | Ebonding |
---|---|---|---|
a | 2N(4)–Au | −5133369.61 | −1100.95 |
b | 4N(4)–Au | −5133267.69 | −999.03 |
c | 2N(3)–Au–4N(1) | −5133340.45 | −1071.79 |
d | 2N(3)–Au–4N(1) | −5133340.23 | −1071.57 |
e | 2N(3)–Au–4N(1) | −5133339.95 | −1071.29 |
f | 2N(3)–Au–4N(1) | −5133339.60 | −1070.94 |
g | 2N(2)–Au–4N(2) | −5133310.78 | −1042.12 |
h | 2N(2)–Au–4N(2) | −5133305.50 | −1036.83 |
i | 2N(1)–Au–4N(3) | −5133283.83 | −1015.17 |
j | 2N(1)–Au–4N(3) | −5133269.90 | −1001.24 |
k | 2N(1)–Au–4N(3) | −5133269.89 | −1001.23 |
l | 2N(1)–Au–4N(3) | −5133287.51 | −1018.85 |
In conclusion from Tables 1–3, the large and similar bonding energies (Ebonding) indicate that these three possible structures of gold(I) complexes and 12 possible structures of gold(III) complexes were all stable structures in the gold electroplating electrolyte. With total bonding energy of −71.92 kJ mol−1, 4N(2)–Au of [Au(DMH)2]− is the most stable structure of the gold(I) complexes in the DMH-based gold electroplating electrolyte. In contrast, the structure 2N(4)–Au with bonding energy of −1100.95 kJ mol−1 was the most possible and stable structure of gold(III) complexes in the DMH-based gold electroplating electrolyte. The higher bonding energy of gold(III) complexes than that of gold(I) complexes indicate that gold(III) complexes were much more stable in the DMH-based gold electroplating electrolyte, so Au3+ is the appropriate gold source for the DMH-based gold electroplating electrolyte to get greater cathodic polarization. Furthermore, with the higher bonding energy, the structure 2N(4)–Au of gold(III) complexes possessed high stability in the gold electroplating electrolyte. The conclusions of the quantum chemical calculations proved that the most stable form of gold ion was [Au(DMH)4]− with 2N(4)–Au coordination structure in the gold electroplating electrolyte.
The results of the quantum chemical calculations confirmed that Au3+ is the appropriate gold source for the DMH-based gold electroplating electrolyte and [Au(DMH)4]− with 2N(4)–Au coordination structure is the most stable form of gold(III)-complexes in the DMH-based gold electroplating electrolyte. The influences of gold ions and DMH on the cyanide-free gold electrodeposition process were further studied to perfect the cyanide-free gold electroplating electrolyte.
The SEM images displayed in Fig. 5 confirmed that the gold electrodeposit obtained from the electrolyte in the absence of DMH was much rougher with bigger particles than that from the electrolytes in the presence of various DMH concentrations. Additionally, the performance of the gold electrodeposits can be improved with the increase of the DMH concentrations, as displayed in Fig. 5(a)–(e). The size of the crystal particles decreased with increasing DMH concentration, however, it only decreased slightly when the DMH concentration was more than 0.30 mol L−1. In summary, ensuring all Au3+ in the solution forms a stable gold(III) complex and considering the influence of DMH on the cathodic polarization, the optimum DMH concentration was determined as 0.30 mol L−1 in the DMH-based cyanide-free gold electroplating electrolyte for subsequent investigations.
Fig. 6 Cyclic voltammograms with different scan rates on a GCE of the gold electroplating electrolytes: (a) without DMH and (b) with 0.30 mol L−1 DMH. |
As displayed in Fig. 6, the reduction peaks of the gold electroplating process in the absence and presence of DMH all shifted to more negative potentials with the increase of the scan rate, indicating that the gold electrodeposition is an irreversible process both with and without DMH as the complexing agent. At the initial scan along the negative direction, the cyclic voltammograms of the gold electrodeposition in the presence of 0.30 mol L−1 DMH all increase at more negative potentials than that in the absence of DMH. Furthermore, the cathodic density peaks increase with the increase of the scan rate in the electrolytes both with and without DMH. The cyclic voltammograms of the electrolytes with DMH as the complexing agent revealed a much lower cathodic density peak than that of the electrolytes in the absence of DMH at the same scan rate.
The linear relationship between the peak current density (jp) and the square root of the scan rate (v1/2) is the characteristic of a diffusion-controlled reaction, which can be expressed as eqn (3):
jp = 0.4958(nF)3/2(αDv)1/2(RT)−1/2c | (3) |
Fig. 7 The peak current density (jp) of cyclic voltammograms with various scan rates versus v1/2 of the gold electroplating electrolytes (a) without DMH and (b) with 0.30 mol L−1 DMH. |
The Delahay equations,48 as shown in eqn (4) and (5), were used to analyze the kinetic features of the gold electrodeposition process based on the cyclic voltammograms.
Ep − Ep/2 = −1.857RT/(αnF) | (4) |
(5) |
Fig. 8 The peak potential (Ep) of cyclic voltammograms with various scan rates versus log v of the gold electroplating electrolytes (a) without DMH and (b) with 0.30 mol L−1 DMH. |
v/mV s−1 | DMH: 0 mol L−1 | DMH: 0.30 mol L−1 | ||
---|---|---|---|---|
Ep − Ep/2/V | α | Ep − Ep/2/V | α | |
10 | −0.022 | 0.771 | −0.041 | 0.414 |
20 | −0.024 | 0.707 | −0.033 | 0.514 |
40 | −0.025 | 0.678 | −0.041 | 0.414 |
60 | −0.027 | 0.628 | −0.051 | 0.333 |
80 | −0.024 | 0.707 | −0.047 | 0.361 |
100 | −0.023 | 0.737 | −0.053 | 0.320 |
Average | 0.705 | 0.393 |
As displayed in Table 4, the transfer coefficients (α) are calculated from the cyclic voltammetry curves according to eqn (4); the α of the gold electroplating in the electrolyte without DMH is 0.705, bigger than that of the electrolyte with DMH as the complexing agent, which is 0.393. The almost doubling of α demonstrates that the electron transfer of the gold electrodeposition process was obviously influenced by the DMH in the electrolyte as the complexing agent for Au3+. Following eqn (3) and the fitted jp vs. v1/2 plots, the diffusion coefficient (D) of ions in the electrolyte were calculated as 8.69 × 10−6 cm2 s−1 and 4.19 × 10−6 cm2 s−1,49 respectively, indicating that D was reduced with the addition of DMH in the electrolyte. It indicated that with DMH as the complexing agent for Au3+ in the gold electroplating electrolyte, the ion transfer was inhibited. In conclusion of the kinetic features based on the cyclic voltammogram measurements, the electrodeposition processes are both irreversible under diffusion control in the studied gold electroplating electrolyte, both with and without 0.30 mol L−1 DMH as the complexing agent. The ion and electron transfer of the gold electrodeposition process were obviously inhibited by DMH in the electrolyte as the complexing agent for Au3+.
(6) |
(7) |
Fig. 10 Current transient curves on the GCE of different gold electroplating electrolytes with various DMH concentrations under different Eap: (a) −0.70 V, (b) −0.75 V, (c) −0.80 V, and (d) −0.85 V. |
As shown in Fig. 9(a) and 10, potentiostatic current transient curves measured in the electrolytes without DMH or with various DMH concentrations; an increase in jm and a decrease in tm are significantly observed when the applied potential (Eap) was set to more negative values. The existence of a current density peak in the current transient curves indicates that the charge transfer reaction controls the reaction kinetics at the early stages of the gold electrodeposition process, showing that there is a nucleation stage preceding the growth of the gold deposit. All these current transient curves can reach a steady state of current density at about 10 s of measurement, demonstrating that after the nucleation stage the reaction is under diffusion control. Moreover, as shown in Fig. 10, when the DMH concentrations increased, a decrease in jm and an increase in tm occurred at the same Eap, and the jm of DMH with 0.10 and 0.20 mol L−1 concentrations are obviously higher than that of other concentrations.
As in the non-dimensional plots from electrolyte without DMH displayed in Fig. 9(b), when the Eap was set to more negative values, the non-dimensional plots of experimental data transferred from following the model of progressive nucleation to instantaneous nucleation. This indicates that the nucleation model of the gold electrodeposition in the absence of DMH complexing agent is mostly influenced by Eap. In contrast with the plots from electrolytes with DMH shown in Fig. 11, the nucleation model of the gold electrodeposition process with various DMH concentrations under different Eap is progressive nucleation. The gold nucleation process is slightly influenced by Eap when the DMH concentrations are 0.10 and 0.20 mol L−1, while it is hardly influenced by Eap when the DMH concentrations are 0.30, 0.45, and 0.60 mol L−1.
During the metal electrodeposition process, based on the instantaneous and progressive nucleation models, parameters of diffusion coefficient, nucleation rate constant, and nucleation density can be acquired using jm and tm with the equations displayed in Table 5, where A is the nucleation rate constant (cm−2 s), N and N∞ are the nuclei number on the electrode (cm−2), n is the number of electrons involved in the reaction (n = 3 in this study), F is the Faraday constant, D is the diffusion coefficient (cm2 s−1), c is the concentration of electroactive species in the electrolyte (c = 0.025 mol L−1), M is the molar mass of gold (M = 196.97 g mol−1), and ρ is the density of gold (ρ = 19.32 g cm−3).
Based on the study of current transient curves, in the electrolyte without DMH, the nucleation process is mostly influenced by Eap, rather than following the instantaneous or progressive nucleation model. The parameters of the electrolyte without DMH cannot be obtained using this method. The nucleation model for the gold electrodeposition process with various DMH concentrations under different Eap is progressive nucleation. The parameters of the nucleation rate of gold on the electrode (AN∞), the density of the saturated gold nucleation number (Ns), and the epitaxial growth rate of the gold crystal (Kv) can be studied with eqn (8)–(10).
(8) |
(9) |
(10) |
With the study of current transient curves in electrolytes with various DMH concentrations, following Table 5 and eqn (8)–(10), the AN∞, Ns, and Kv can be obtained, as listed in Table 6.
E (V) | DMH (mol L−1) | jm (mA cm−2) | tm (s) | AN∞ (×107 cm−2 s) | Ns (×107 cm−2) | Kv (×10−8 mol s−1 cm−2) |
---|---|---|---|---|---|---|
−0.70 | 0.10 | −13.994 | 0.856 | 4.89 | 2.43 | 4.83 |
−0.70 | 0.20 | −11.111 | 1.614 | 1.16 | 1.08 | 3.84 |
−0.70 | 0.30 | −9.036 | 2.034 | 0.87 | 1.03 | 3.12 |
−0.70 | 0.45 | −8.020 | 2.350 | 0.72 | 0.98 | 2.77 |
−0.70 | 0.60 | −7.592 | 2.389 | 0.76 | 1.06 | 2.62 |
−0.75 | 0.10 | −19.84 | 0.509 | 11.56 | 3.41 | 6.86 |
−0.75 | 0.20 | −17.891 | 0.587 | 9.27 | 3.16 | 6.18 |
−0.75 | 0.30 | −12.579 | 1.106 | 2.80 | 1.80 | 4.35 |
−0.75 | 0.45 | −10.749 | 1.422 | 1.81 | 1.49 | 3.71 |
−0.75 | 0.60 | −10.494 | 1.395 | 2.01 | 1.62 | 3.63 |
−0.80 | 0.10 | −26.582 | 0.265 | 45.66 | 7.01 | 9.18 |
−0.80 | 0.20 | −23.694 | 0.323 | 31.73 | 5.94 | 8.19 |
−0.80 | 0.30 | −14.522 | 0.806 | 5.44 | 2.54 | 5.02 |
−0.80 | 0.45 | −12.600 | 0.930 | 4.70 | 2.53 | 4.35 |
−0.80 | 0.60 | −13.141 | 0.898 | 4.80 | 2.50 | 4.54 |
−0.85 | 0.10 | −32.512 | 0.176 | 104.17 | 10.63 | 11.23 |
−0.85 | 0.20 | −27.813 | 0.207 | 87.50 | 10.50 | 9.61 |
−0.85 | 0.30 | −17.325 | 0.581 | 10.20 | 3.43 | 5.98 |
−0.85 | 0.45 | −14.692 | 0.710 | 7.77 | 3.20 | 5.08 |
−0.85 | 0.60 | −14.282 | 0.742 | 7.21 | 3.10 | 4.93 |
As shown in Table 6, at the same DMH concentration, AN∞, Ns, and Kv all increased with the negative shift of Eap, indicating that the growth rate of the crystal nucleus can be promoted by the negative shift of Eap. Additionally, at the same Eap, AN∞, Ns, and Kv all decreased with the increase of the DMH concentration, indicating the inhibition influence of DMH on nucleation. However, when the DMH concentration reached 0.30 mol L−1, the downtrend of AN∞, Ns, and Kv is diminished with the increase of DMH concentration, indicating that with the DMH concentration higher than 0.30 mol L−1, the growth rate of crystal nucleus was hardly influenced by the DMH concentration.
The cathodic polarization of the gold electrodeposition was significantly enhanced with DMH as the complexing agent; the 500 mV negative shift of the onset potential demonstrates a significant kinetic change compared with the electrolyte without DMH. The sharp increase in the potential of the cathodic polarization curves and the micromorphology of the gold electrodeposits was improved with the increase of the DMH concentration. Considering the influence of the DMH concentration on the cathodic polarization and the gold electrodeposit micromorphology, the optimum DMH concentration was determined as 0.30 mol L−1.
In conclusion of the kinetic features based on cyclic voltammograms measurements, the electrodeposition processes are both irreversible under diffusion control in the studied gold electroplating electrolyte, both with and without 0.30 mol L−1 DMH as the complexing agent. The transfer coefficients of the gold electroplating electrolytes in the absence and presence of DMH were 0.705 and 0.393, respectively. The diffusion coefficients were calculated as 8.69 × 10−6 cm2 s−1 and 4.19 × 10−6 cm2 s−1 without and with DMH, respectively. The ion and electron transfer of gold electrodeposition process were obviously inhibited by DMH in the electrolyte as the complexing agent for Au3+. The gold electrodeposition process displayed progressive nucleation with DMH as the complexing agent according to the SH nucleation model with various applied potentials. The growth rate of the crystal nucleus was reduced by DMH and promoted by the negative shift of Eap, while when the DMH concentration was higher than 0.30 mol L−1, the growth rate of the crystal nucleus was hardly influenced by the DMH concentration.
This journal is © The Royal Society of Chemistry 2018 |