Next Article in Journal
Current Status of the Bioactive Properties of Royal Jelly: A Comprehensive Review with a Focus on Its Anticancer, Anti-Inflammatory, and Antioxidant Effects
Next Article in Special Issue
Synthesis, Empirical and Theoretical Investigations on New Histaminium Bis(Trioxonitrate) Compound
Previous Article in Journal
Improved Antioxidant and Mechanical Properties of Food Packaging Films Based on Chitosan/Deep Eutectic Solvent, Containing Açaí-Filled Microcapsules
Previous Article in Special Issue
Comprehensive Study of the Ammonium Sulfamate–Urea Binary System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation on Gold–Ligand Interaction for Complexes from Gold Leaching: A DFT Study

School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(3), 1508; https://doi.org/10.3390/molecules28031508
Submission received: 22 December 2022 / Revised: 27 January 2023 / Accepted: 1 February 2023 / Published: 3 February 2023
(This article belongs to the Special Issue Advances in Density Functional Theory (DFT) Calculation)

Abstract

:
Gold leaching is an important process to extract gold from ore. Conventional alkaline cyanide process and alternative nontoxic lixiviants including thiosulfate, thiourea, thiocyanate, and halogen have been widely investigated. However, density functional theory (DFT) study on the gold complexes Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 required for discovering and designing new highly efficient and environmentally friendly gold leaching reagents is lacking, which is expected to support constructive information for the discovery and designation of new high-efficiency and environmentally friendly gold leaching reagents. In this study, the structure information, electron-transferring properties, orbital interaction, and chemical bond composition for complexes Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 depending on charge decomposition analysis (CDA), natural bond orbital (NBO), natural resonance theory (NRT), electron localization function (ELF), and energy decomposition analysis (EDA) were performed based on DFT calculation. The results indicate that there is not only σ-donation from ligand to Au+, but also electron backdonation from Au+ to ligands, which strengthens the coordinate bond between them. Compared with Cl, ligands CN, S2O32−, SC(NH2)2, and SCN have very large covalent contribution to the coordinate bond with Au+, which explains the special stability of Au-CN and Au-S bonds. The degree of covalency and bond energy in Au–ligand bonding decreases from Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, to AuCl2, which interprets the stability of the five complexes: Au(CN)2 > Au(S2O3)23− > Au[SC(NH2)2]2+ > Au(SCN)2 > AuCl2.

1. Introduction

As one of the most important noble metals, gold is not only used as an ornament, a concentrated form of wealth, and money, it also has critical technological applications in the electronic, computer, and space fields [1]. To meet the use of gold in society, gold extraction has shown its importance. Au possesses a nontoxic nature, stability against metal leaching, and resistance to excessive oxidation by O2 [2]. However, gold becomes unstable and dissolves in the positive potential region of electrochemical processes, especially in solutions containing CN, Cl, and sulfur-containing species such as S2O32−, SC(NH2)2, and SCN, which make the extraction of gold from ore feasible [3,4,5].
The conventional alkaline cyanide process is the dominant method for gold recovery from ores due to its simplicity and high efficiency. The gold leaching process by cyanide is represented by Equation (1), where complex Au(CN)2 is formed. However, the cyanide process has a limitation due to its environmental issues. In the past decades, alternative nontoxic lixiviants for gold have drawn extensive attention, including thiosulfate (S2O32−), thiourea (SC(NH2)2), thiocyanate (SCN), and halogen [6,7,8,9].
4 Au + 8 CN + O 2 + 2 H 2 O = 4 Au CN 2 + 4 OH
Thiosulfate is the most promising candidate since it is environmentally friendly, inexpensive, and has good selectivity towards gold (Equation (2)). Normally, thiosulfate is used with an ammonia–copper solution, since copper sulfate can be a catalyst and ammonia can stabilize thiosulfate. Although it has environmental benefit, it shows a high reagent consumption [10,11,12,13,14]. Thiourea is an organosulfur compound that forms a soluble cationic complex with gold in a sulfuric acid solution. Fe3+ is added to facilitate the oxidation of gold (Equation (3)). However, the use of thiourea would cause a high reagent consumption and a surface passivation [1,2,4,7,15]. Compared to cyanide, thiocyanate is far less toxic and is safer for the environment. Compared to thiourea and thiosulfate, thiocyanate leaching systems require lower reagent consumption and cost, and have better stability than thiourea. Fe3+ is normally used as the oxidant (Equation (4)). However, the reaction occurs in an acidic environment, which would be corrosive to equipment [8,16,17,18]. Chlorine is also used as an alternative lixiviant for gold due to its high dissolution rate (Equation (5)), but chlorine is strongly corrosive, highly volatile, and hazardous [9,19,20,21]. The stability constants and standard reduction potentials for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 are listed in Table 1 [1,5,14]. The higher the stability constant and the lower the standard reduction potential, the more stable the complex is, so the stability of the five complexes is: Au(CN)2 > Au(S2O3)23− > Au[SC(NH2)2]2+ > Au(SCN)2 > AuCl2. Generating a stable gold–ligand complex is the prerequisite for gold leaching.
4 Au + 8 S 2 O 3 2 + O 2 + 2 H 2 O = 4 [ Au ( S 2 O 3 ) 2 ] 3 + 4 OH
Au + 2 SC ( NH 2 ) 2 + Fe 3 + = Au [ SC ( NH 2 ) 2 ] 2 + + Fe 2 +
Au + 2 SCN + Fe 3 + = Au ( SCN ) 2 + Fe 2 +
2 Au + 2 Cl + Cl 2 = 2 AuCl 2
In addition to the alternative lixiviants mentioned above, new highly efficient and environmentally friendly reagents have been discovered or designed, such as ionic liquid, dicyanamide, ferrocyanide, etc. [22,23,24]. However, those new gold leaching reagents are in the preliminary research phase. The bonding properties of Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 in the quantum chemistry perspective would support instructive information for the discovery and design of new highly efficient and environmentally friendly gold leaching reagents. Although there are many studies on gold leaching by thiosulfate (S2O32−), thiourea (SC(NH2)2), thiocyanate (SCN), and halogen, there is a lack of systematic wavefunction analysis of the gold complexes Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2.
Therefore, in this study, the structure information, electron transferring properties, orbital interaction, chemical bond composition, and bond energy for complexes Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 depending on charge decomposition analysis (CDA), natural bond orbital (NBO), natural resonance theory (NRT), electron localization function (ELF), and energy decomposition analysis (EDA) have been performed based on density functional theory (DFT) calculation, according to which the stability of the five complexes is also analyzed. The information is expected to support constructive instruction to discover the potential gold leaching reagent.

2. Results and Discussion

2.1. Optimized Geometries of Gold Complexes

The optimized geometries of five gold complexes are shown in Figure 1, whose calculated geometry and structural parameters are listed in Table 2. Au(CN)2 and AuCl2 have a D∞h symmetry with a linear structure, while Au(S2O3)23−, Au[SC(NH2)2]2+, and Au(SCN)2 have a C2 symmetry where Au+ and the ligands are not in the same surface and there is a slight bend for S-Au-S whose angle is around 176~177°. All the bond lengths in the complexes are within the range of the M–X covalent single-bond lengths [25]. The bond lengths from the two levels of theory are consistent with each other within 0.005 Å or better for all the species. The bond lengths and angles of Au-C in Au(CN)2, Au-S in Au(SCN)2, and Au-Cl in AuCl2 at the level of energy-consistent quasi-relativistic pseudopotential basis set TPSSh/aug-cc-pVTZ are consistent with those calculated at the level of scalar relativistic basis set SR-ZORA PBE/TZ2P reported previously [26,27]. The results indicate that the bond length of Au-C in Au(CN)2 is the shortest among the four complexes, which is 1.987 Å, and that of Au-X in the other four complexes is similar, which is around 2.3 Å. As bond strength increases, the atoms in the bond are pulled more tightly together. Therefore, the bond energy of Au-C in Au(CN)2 is larger than that of Au-S and Au-Cl in other four complexes. The charge population analyses are listed in Table 3. All methods (Hirshfeld, ADCH, CM5, and CHELPG) show that the Au atom carries little charge, while the negative charge is mainly distributed on ligands. The atomic charge on Au is less than 1, indicating that there is electron transfer from ligands to Au atom.

2.2. Charge Decomposition Analysis (CDA)

Charge decomposition analysis (CDA) is a partitioning scheme to analyze donor–acceptor interactions. It was developed by Dapprich and Frenking as a quantitative expression of the well-known Dewar–Chatt–Duncanson DCD model, which describes the metal–ligand interactions in terms of donation and backdonation chemical mechanisms [32,33,34]. The CDA method has proven to be very useful for estimating the relative strength of charge donation and backdonation in a series of transition metal complexes [35,36].
Based on CDA, the number of electrons donated from ligands to Au+ and the number of electrons backdonated from Au+ to ligands are listed in Table 4. It indicates that the coordinate bonds of all five gold complexes are not a pure σ-donation, and there is an electron backdonation from Au+ to ligands, although the number of electrons donated from ligands to Au+ is larger than that backdonated from Au+ to ligands. The backdonation strengthened the bond between Au+ and the ligands [37]. The percentage contribution of the backdonation relative to the donation has the tendency to decrease from Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, to AuCl2, which corresponds to the stability of the five complexes. Especially in Au(CN)2, its donation and backdonation are both the largest when compared with the other complexes, which would explain its high stability among those complexes.
The isocontour surfaces of molecular orbitals (MOs) from HOMO to HOMO-8 for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 derived from TPSSh/aug-cc-pVTZ calculation are presented in Figure 2. Based on CDA, the main MOs contributing to σ-donations are highlighted by the red-dash rectangular frames shown in Figure 2.
The MOs that contribute to σ-donation and their orbital composition from Au+ and ligands are shown in Figure 3. It can be observed that the spherical s orbital in Au+ participates in the formation of the σ bond for all five complexes. For Au(CN)2, there are two orbitals mainly contributing to the σ-bonding between Au+ and (CN)22−, orbital 23 (HOMO, σg) and orbital 18 (HOMO-5, σu). Orbital 23 is composed of orbital 10 (6s) from Au+ and orbital 13 from (CN)22−, and orbital 18 consisting of orbital 11 (6px) from Au+ and orbital 14 from (CN)22−. For Au(S2O3)23−, the number of electrons donated from ligands to Au+ is the most in orbital 64 (HOMO-3, πu), which is assembled by orbital 10 from Au+ and orbital 55 from (S2O3)24−. For Au[SC(NH2)2]2+, in orbital 47 (HOMO-2, σg), [SC(NH2)2]2 transfers the most electrons to Au+ through the σ-bonding between orbital 10 from Au+ and orbital 39 from [SC(NH2)2]2. For Au(SCN)2, σ-bonding between orbital 10 from Au+ and orbital 27 from (SCN)22− forms the orbital 37 (HOMO-2, σg), which contributes the most for the electrons donated from (SCN)22− to Au+. For AuCl2, σ-bonds in orbital 25 (HOMO-2, σg) and orbital 20 (HOMO-7, σu) mainly contribute to the σ-donation between Au+ and Cl22–. Orbital 10 from Au+ and orbital 13 from Cl22− compose orbital 25, and orbital 12 from Au+ and orbital 18 from Cl22– constitute orbital 20.
The MOs that contribute to π-backdonation and their orbital composition from Au+ and ligands are displayed in Figure 4. It shows that d orbitals from Au+ participate in the formation of the π bond in all five complexes. For Au(CN)2, the π-backdonation is mainly generated in orbital 13, which is composed of orbital 6 (5d) in Au+ and orbital 25 in (CN)22−. For Au(S2O3)23−, orbital 49 is in charge of the π-backdonation from Au+ to ligands, which is mainly composed of orbital 8 (5d) in Au+ and orbital 62 in (S2O3)24−. For Au[SC(NH2)2]2+, the π-backdonation from Au+ to ligands occurs in orbital 40, which is the combination of orbital 7 (5d) in Au+ and orbital 49 in [SC(NH2)2]2. For Au(SCN)2, the π-backdonation takes place from orbital 6 in Au+ to orbital 38 in (SCN)22−, which forms orbital 32 in the complex. For AuCl2, orbital 9 (5d) in Au+ and orbital 23 in Cl22− forms bonding orbital 17 in the complex, which is in charge with the π-backdonation from metal to ligands.
In the complexes, σ-bonds in those orbitals are all linear or nearly linear. The Au 5d orbitals transform as σg + πg + δg in D∞h symmetry. Owing to the strong relativistic effects in gold, the sd z 2 hybridization is enhanced [38]. The Au 5d atomic orbitals would not be involved in net bonding with main-group elements if there were no sd hybridizations. The sd hybridization results in two hybrid orbitals, which are 180° relative to each other. The geometrical characters of the two sd hybrid orbitals lead to the most effective overlaps with the AO in ligands. The linear geometry results in the maximum overlap between the AOs and the reduction in repulsion between the two negatively charged ligands [27], which is beneficial to stabilize the complexes.
The CDA results give an important insight into the orbital interaction and electron transfer between Au+ and ligands. Based on the donation and backdonation analysis in Au-C (Au(CN)2), Au-S (Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2), and Au-Cl (AuCl2), the strength of σ-bonds for the complexes can be explained. From the detailed analysis of orbital composition, the composition and shapes of the σ-bonds and π-backdonation bonds are demonstrated.

2.3. Chemical Bonding Analyses

The natural bond orbital (NBO) and natural resonance theory (NRT) analyses were performed to illustrate the participation of atomic orbitals (AOs) in the bonding orbitals of the complexes, which are listed in Table 5. The NBO results demonstrate significant s-d hybridization in Au+ for all five complexes. Additionally, except for Au(CN)2, a small amount of orbital p (0.01) also participates in the hybridization. For complexes Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2, the percentage of Au+ is similar, which is around 24~25%, while that in AuCl2 is less, which is no more than 20%.
Natural resonance theory (NRT) is a theoretical method based on quantum chemical calculations to describe molecules with significant resonance structures in terms of classic valence-bond concepts [39]. The NRT approach allows the expansion of total wave function in terms of chemically intuitive resonance structures [40].
The NRT analysis indicates the covalent and ionic contribution to the Au-L bond. For complexes Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, and Au(SCN)2, they have similar covalent contribution, which is above 0.4, showing covalent character. While for AuCl2, its covalent bond order is 0.176, which shows more ionic character due to the large electronegativity of Cl22−. More covalent contribution indicates that the interaction between metal and ligands is stronger, namely, the bonds would be more stable. For the three complexes containing the Au-S bond, their covalent contribution slightly decreases from Au(S2O3)23−, Au[SC(NH2)2]2+, to Au(SCN)2, which corresponds to their stability.
Electron localization function (ELF) measures the excess of kinetic energy density due to the Pauli repulsion. In the region of space where the Pauli repulsion is strong, ELF is close to 1, whereas where the probability of finding the same-spin electrons close together is high, ELF tends to 0. Regions where the value of ELF is close to 1 correspond to well-localized electrons [41]. ELF is a chemically intuitive tool to estimate the probability of finding electron pairs in space, and an ELF analysis can show the difference among covalent bonds, ionic bonds, and lone pairs [39].
The calculated ELFs shown in Figure 5 reveal the distributions of the electron pairs in the real space of gold complexes. It is revealed from the ELF analyses that electron pairing densities in the Au-C bond in Au(CN)2 and Au–S bonds in Au(S2O3)23−, Au[SC(NH2)2]2+, and Au(SCN)2 are all around 0.3, while that in the Au-Cl bond in AuCl2 is ~0.2. Since a high value of ELF corresponds to well-localized electrons, the electrons in Au-C bonds (Au(CN)2) and Au-S bonds (Au(S2O3)23−, Au[SC(NH2)2]2+, and Au(SCN)2) are more localized than that in Au-Cl bonds (AuCl2). Localized electrons are associated with covalent bonds, so Au-C bonds in Au(CN)2 and Au-S bonds in Au(S2O3)23−, Au[SC(NH2)2]2+, and Au(SCN)2 show more covalent character than Au-Cl bonds in AuCl2, which is in accord with the NRT analyses results.
NBO analysis supports the information of orbital composition. The NRT and ELF results demonstrate the bond properties between metal and ligands, which gives an important insight into the special stability of Au-CN and Au-S bonds. Compared with Cl, ligands CN, S2O32−, SC(NH2)2, and SCN have a very large covalent contribution to the coordinate bond with Au+, which explains the small stability constant of AuCl2.

2.4. Energy Decomposition Analysis (EDA)

Bonding energy is an important criterion for the stability of a chemical bond. The decomposition of the bonding energy between Au and the ligand provides further insight into the nature of the Au–ligand bonding. According to energy decomposition analysis (EDA), the interaction energy between two fragments can be decomposed into the sum of electrostatic interactions, Pauli repulsion, and the orbital interaction energy [26]:
Δ E bonding = Δ E electrostatic + Δ E pauli + Δ E orbint = Δ E steric + Δ E orbint
The steric interaction energy is the sum of the electrostatic interaction energy and Pauli repulsion. The orbital interaction can be viewed as an electron transfer from one fragment to the other or as the orbital mixing between two fragments, which indicates the degree of covalency between the fragments. An ionic bond can be described as the interaction between two frozen densities plus intramolecular relaxations, which is ΔEsteric in EDA. Therefore, the orbital interaction energy in the total bonding energy reflects the degree of covalency in the chemical bond between two atoms or fragments.
The EDA results for one-ligand deposition are shown in Table 6. It displays that the steric interaction is positive, indicating that the interaction introduced by steric is repulsive. The value of steric repulsion is AuCl2 < Au(CN)2 < Au(SCN)2 < Au[SC(NH2)2]2+ < Au(S2O3)23−, which can be explained by the ligand size and the symmetry of the complex. The ligand size increases from Cl, CN, SCN, SC(NH2)2, to S2O34−. Additionally, for AuCl2 and Au(CN)2, all atoms are in the same surface and the complex is in a linear shape. However, for Au(SCN)2, Au[SC(NH2)2]2+,and Au(S2O3)23−, the metal atom and ligand atoms are not in the same surface and the S-Au-S bond is not linear, which also lead to the steric repulsion; the orbital interactions are all negative, which indicates the attractive interaction between Au+ and ligands. The value of orbital interaction increases from Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, to AuCl2 meaning the degree of covalency between the fragments decreases. Combined the steric interaction and orbital interaction, the total bonding energy is negative, revealing the attraction between L-Au and ligand. In addition, it also indicates that the degree of covalency is larger than that of ionic. The value of total interaction is in the order Au(CN)2 > Au(S2O3)23− > Au[SC(NH2)2]2+ > Au(SCN)2, except AuCl2 whose total interaction is a little larger than that of Au(SCN)2 because of its low steric interaction.
The EDA results for two-ligand deposition are illustrated in Table 7, which has similar tendency with that in Table 6. The steric interaction is positive, while the orbital interactions are all negative. Compared with the values in Table 6, it is found that the energy for two-ligands deposition is more than twice of that for one-ligand deposition. The value of orbital interaction also increases from Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, to AuCl2, meaning the degree of covalency between the fragments decreases. The value of total interaction is also in the order Au(CN)2 > Au(S2O3)23− > Au[SC(NH2)2]2+ > Au(SCN)2 > AuCl2.
The covalency in the coordinate bond and the total interaction value are related to the stability of the complex. Based on the results in Table 6 and Table 7, it can be concluded that the EDA results are consistent with the stability of the five complexes, which is Au(CN)2 > Au(S2O3)23− > Au[SC(NH2)2]2+ > Au(SCN)2 > AuCl2.

3. Computational Methods

All DFT calculations were performed in the Gaussian 16 [42] program by using hybrid meta-GGA functional TPSSh [43] with aug-cc-pVTZ basis set [44,45,46] for nonmetallic atoms and Stuttgart pseudopotentials ECP60MDF with corresponding augmented valence basis aug-cc-pVTZ-PP for metallic atom Au [47,48]. For comparison, all complexes were also optimized at high-quality B2PLYP/def2-QZVP level [49,50]. The polarizable continuum model (PCM) model [51] was used to incorporate solvation effects with water for all calculations. All of the geometries were confirmed to be minima by performing a vibrational frequency calculation and had no imaginary frequency.
The natural bond orbital (NBO) and natural resonance theory (NRT) calculations were implemented in the NBO 5.0 package [52]. The charge decomposition analysis (CDA) [33,53], electron localization function (ELF) [54], and energy decomposition analysis (EDA) calculations were performed in Multiwfn 3.8 (dev) software [55].

4. Conclusions

The charge decomposition analysis (CDA), natural bond orbital (NBO), natural resonance theory (NRT), electron localization function (ELF), and energy decomposition analysis (EDA) for complexes Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 were performed by DFT calculation, based on which their structure information, electron transferring properties, orbital interaction, chemical bond composition, and bond energy were comprehensively elaborated. This study aimed to support the bonding information to the potential gold leaching reagent. The following conclusions can be drawn:
(1)
Based on the CDA results, the orbital interaction and electron transferring between Au+ and ligands for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 are interpreted. There is not only σ-donation from ligand to Au+, but also electron backdonation from Au+ to ligands, which strengthens the coordinate bond between them. The percentage of π-backdonation in Au(CN)2 is the largest contributing to its high stability.
(2)
From the perspective of NRT and ELF, compared with Cl, ligands CN, S2O32−, SC(NH2)2, and SCN have very large covalent contribution to the coordinate bond with Au+, which explains the special stability of Au-CN and Au-S bonds.
(3)
The decomposition of the bonding energy between Au and the ligand provides further insight into the nature of the Au–ligand bonding. The EDA results show the increased value of orbital interaction from Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, to AuCl2, corresponding to the decreased degree of covalency in Au–ligand bonding, which interprets the stability of the five complexes: Au(CN)2 > Au(S2O3)23− > Au[SC(NH2)2]2+ > Au(SCN)2 > AuCl2.

Author Contributions

Conceptualization, methodology, software, writing—original draft preparation, and writing—review and editing, N.Z.; validation, investigation, and resources, J.K.; project administration and funding acquisition, C.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant number: 51974016).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, J.; Miller, J.D. A Review of Gold Leaching in Acid Thiourea Solutions. Miner. Process. Extr. Metall. Rev. 2006, 27, 177–214. [Google Scholar] [CrossRef]
  2. Li, L.; Pan, C.; Shan, J.; She, W.; You, X.; Ji, C.; Gao, Q. Pit-Induced Electrochemical Layer Dissolution and Wave Propagation on an Au(111) Surface in an Acidic Thiourea Solution. J. Phys. Chem. C 2020, 124, 19112–19118. [Google Scholar] [CrossRef]
  3. Zhang, Y.; Xu, B.; Cui, M.; Li, Q.; Liu, X.; Jiang, T.; Lyu, X. Thiosulfate leaching of gold catalyzed by hexaamminecobalt(III): Electrochemical behavior and mechanisms. Electrochim. Acta 2021, 399, 139393. [Google Scholar] [CrossRef]
  4. Li, W.-J.; Zhou, H.; Bai, A.-P.; Song, Y.-S.; Cai, L.-L.; Zheng, S.-L.; Zhang, Q.-D.; Cao, S. Electrochemical adsorption and passivation on gold surface in alkaline thiourea solutions. Rare Met. 2020, 39, 951–958. [Google Scholar] [CrossRef]
  5. Birich, A.; Stopic, S.; Friedrich, B. Kinetic Investigation and Dissolution Behavior of Cyanide Alternative Gold Leaching Reagents. Sci. Rep. 2019, 9, 7191. [Google Scholar] [CrossRef]
  6. Zhang, X.M.; Senanayake, G. A Review of Ammoniacal Thiosulfate Leaching of Gold: An Update Useful for Further Research in Non-cyanide Gold Lixiviants. Miner. Process. Extr. Metall. Rev. 2016, 37, 385–411. [Google Scholar] [CrossRef]
  7. Whitehead, J.; Zhang, J.; McCluskey, A.; Lawrance, G. Comparative leaching of a sulfidic gold ore in ionic liquid and aqueous acid with thiourea and halides using Fe(III) or HSO5−oxidant. Hydrometallurgy 2009, 98, 276–280. [Google Scholar] [CrossRef]
  8. Azizitorghabeh, A.; Wang, J.; Ramsay, J.A.; Ghahreman, A. A review of thiocyanate gold leaching—Chemistry, thermodynamics, kinetics and processing. Miner. Eng. 2021, 160, 106689. [Google Scholar] [CrossRef]
  9. Ilyas, S.; Srivastava, R.R.; Kim, H. Gold recovery from secondary waste of PCBs by electro-Cl2 leaching in brine solution and solvo-chemical separation with tri-butyl phosphate. J. Clean. Prod. 2021, 295, 126389. [Google Scholar] [CrossRef]
  10. Senanayake, G. Gold leaching by thiosulphate solutions: A critical review on copper(II)–thiosulphate–oxygen interactions. Miner. Eng. 2005, 18, 995–1009. [Google Scholar] [CrossRef]
  11. Grosse, A.C.; Dicinoski, G.W.; Shaw, M.J.; Haddad, P. Leaching and recovery of gold using ammoniacal thiosulfate leach liquors (a review). Hydrometallurgy 2003, 69, 1–21. [Google Scholar] [CrossRef]
  12. Sitando, O.; Senanayake, G.; Dai, X.; Nikoloski, A.; Breuer, P. A review of factors affecting gold leaching in non-ammoniacal thiosulfate solutions including degradation and in-situ generation of thiosulfate. Hydrometallurgy 2018, 178, 151–175. [Google Scholar] [CrossRef]
  13. Xu, B.; Kong, W.; Li, Q.; Yang, Y.; Jiang, T.; Liu, X. A Review of Thiosulfate Leaching of Gold: Focus on Thiosulfate Consumption and Gold Recovery from Pregnant Solution. Metals 2017, 7, 222. [Google Scholar] [CrossRef]
  14. Aylmore, M.G.; Muir, D.M. Thiosulfate leaching of gold—A review. Miner. Eng. 2001, 14, 135–174. [Google Scholar] [CrossRef]
  15. Yu, B.; Liu, Y.; Peng, X.; Hua, S.; Zhou, G.; Yan, K.; Liu, Y. Synthesis, characterization, and antitumor properties of Au(i)-thiourea complexes. Metallomics 2020, 12, 104–113. [Google Scholar] [CrossRef]
  16. Yang, X.; Moats, M.S.; Miller, J.D.; Wang, X.; Shi, X.; Xu, H. Thiourea–thiocyanate leaching system for gold. Hydrometallurgy 2011, 106, 58–63. [Google Scholar] [CrossRef]
  17. Kholmogorov, A.; Kononova, O.; Pashkov, G.; Kononov, Y. Thiocyanate solutions in gold technology. Hydrometallurgy 2002, 64, 43–48. [Google Scholar] [CrossRef]
  18. Azizitorghabeh, A.; Mahandra, H.; Ramsay, J.; Ghahreman, A. Gold leaching from an oxide ore using thiocyanate as a lixiviant: Process optimization and kinetics. ACS Omega 2021, 6, 17183–17193. [Google Scholar] [CrossRef]
  19. Nikoloski, A.; Stockton, B. Application of alternative lixiviants for secondary heap leaching of gold. In Proceedings of the 7th Mill Operators Conference, Kalgoorlie, Australia, 12–14 October 2000. [Google Scholar]
  20. Dönmez, B.; Sevim, F.; Çolak, S. A Study on Recovery of Gold from Decopperized Anode Slime. Chem. Eng. Technol. 2001, 24, 91–95. [Google Scholar] [CrossRef]
  21. Wang, S.; Li, L.; Wang, H.; Wu, G.-D. Extraction of platinum and gold from copper anode slimes by a process of chlorinating roasting followed by chlorinating leaching. J. Min. Met. Sect. B Met. 2020, 56, 193–202. [Google Scholar] [CrossRef]
  22. Rodríguez, M.; Ayala, L.; Robles, P.; Sepúlveda, R.; Torres, D.; Carrillo-Pedroza, F.R.; Jeldres, R.I.; Toro, N. Leaching Chalcopyrite with an Imidazolium-Based Ionic Liquid and Bromide. Metals 2020, 10, 183. [Google Scholar] [CrossRef]
  23. Li, G.-Z.; Kou, J.; Xing, Y.; Hu, Y.; Han, W.; Liu, Z.-Y.; Sun, C.-B. Gold-leaching performance and mechanism of sodium dicyanamide. Int. J. Miner. Met. Mater. 2021, 28, 1759–1768. [Google Scholar] [CrossRef]
  24. Liu, Z.; Kou, J.; Xing, Y.; Sun, C. Recovery of Gold from Ore with Potassium Ferrocyanide Solution under UV Light. Minerals 2021, 11, 387. [Google Scholar] [CrossRef]
  25. Pyykkö, P.; Atsumi, M. Molecular single-bond covalent radii for elements 1-118. Chemistry 2009, 15, 186–197. [Google Scholar] [CrossRef] [PubMed]
  26. Xiong, X.G.; Wang, Y.L.; Xu, C.Q.; Qiu, Y.H.; Wang, L.S.; Li, J. On the gold-ligand covalency in linear [AuX2](−) complexes. Dalton Trans. 2015, 44, 5535–5546. [Google Scholar] [CrossRef]
  27. Xu, C.Q.; Xiong, X.G.; Li, W.L.; Li, J. Periodicity and Covalency of [MX2] (M=Cu, Ag, Au, Rg;X=H, Cl, CN) Complexes. Eur. J. Inorg. Chem. 2016, 2016, 1395–1404. [Google Scholar] [CrossRef]
  28. Hirshfeld, F.L. Bonded-atom fragments for describing molecular charge densities. Theor. Chim. Acta 1977, 44, 129–138. [Google Scholar] [CrossRef]
  29. Lu, T.; Chen, F. Atomic dipole moment corrected Hirshfeld population method. J. Theor. Comput. Chem. 2012, 11, 163–183. [Google Scholar] [CrossRef]
  30. Marenich, A.V.; Jerome, S.V.; Cramer, C.J.; Truhlar, D.G. Charge Model 5: An Extension of Hirshfeld Population Analysis for the Accurate Description of Molecular Interactions in Gaseous and Condensed Phases. J. Chem. Theory Comput. 2012, 8, 527–541. [Google Scholar] [CrossRef]
  31. Breneman, C.M.; Wiberg, K. Determining atom-centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis. J. Comput. Chem. 1990, 11, 361–373. [Google Scholar] [CrossRef]
  32. Chatt, J.; Duncanson, L.A. 586. Olefin co-ordination compounds. Part III. Infrared spectra and structure: Attempted preparation of acetylene complexes. J. Chem. Soc. 1953, 2939–2947. [Google Scholar] [CrossRef]
  33. Dapprich, S.; Frenking, G. Investigation of donor-acceptor interactions: A charge decomposition analysis using fragment molecular orbitals. J. Phys. Chem. 1995, 99, 9352–9362. [Google Scholar] [CrossRef]
  34. Hernández, M.G.; Beste, A.; Frenking, G.; Illas, F. Charge decomposition analysis of the chemisorption bond. Chem. Phys. Lett. 2000, 320, 222–228. [Google Scholar] [CrossRef]
  35. Ehlers, A.W.; Dapprich, S.; Vyboishchikov, S.F.; Frenking, G. Structure and Bonding of the Transition-Metal Carbonyl Complexes M (CO) 5L (M= Cr, Mo, W) and M (CO) 3L (M= Ni, Pd, Pt; L= CO, SiO, CS, N2, NO+, CN, NC, HCCH, CCH2, CH2, CF2, H2). Organometallics 1996, 15, 105–117. [Google Scholar] [CrossRef]
  36. Szilagyi, R.K.; Frenking, G.J.O. Structure and Bonding of the Isoelectronic Hexacarbonyls [Hf(CO)6]2−, [Ta(CO)6], W(CO)6, [Re(CO)6]+, [Os(CO)6]2+, and [Ir(CO)6]3+: A Theoretical Study. Organometallics 1997, 16, 4807–4815. [Google Scholar] [CrossRef]
  37. Crabtree, R.H. The Organometallic Chemistry of the Transition Metals; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  38. Pyykko, P.J.C.R. Relativistic effects in structural chemistry. Chem. Rev. 1988, 88, 563–594. [Google Scholar] [CrossRef]
  39. Ning, C.G.; Xiong, X.G.; Wang, Y.L.; Li, J.; Wang, L.S. Probing the electronic structure and chemical bonding of the “staple” motifs of thiolate gold nanoparticles: Au(SCH3)2− and Au2(SCH3)3. Phys. Chem. Chem. Phys. 2012, 14, 9323–9329. [Google Scholar] [CrossRef]
  40. Glendening, E.D.; Weinhold, F. Natural resonance theory: I. General formalism. J. Comput. Chem. 1998, 19, 593–609. [Google Scholar] [CrossRef]
  41. Poater, J.; Duran, M.; Solà, M.; Silvi, B. Theoretical Evaluation of Electron Delocalization in Aromatic Molecules by Means of Atoms in Molecules (AIM) and Electron Localization Function (ELF) Topological Approaches. Chem. Rev. 2005, 105, 3911–3947. [Google Scholar] [CrossRef]
  42. Frisch, M.E.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.P.G.A.; Petersson, G.A.; Nakatsuji, H.J.R.A.; et al. Gaussian 16; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  43. Tao, J.; Perdew, J.P.; Staroverov, V.N.; Scuseria, G.E. Climbing the Density Functional Ladder: Nonempirical Meta–Generalized Gradient Approximation Designed for Molecules and Solids. Phys. Rev. Lett. 2003, 91, 146401. [Google Scholar] [CrossRef] [Green Version]
  44. Dunning, T.H., Jr. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. [Google Scholar] [CrossRef]
  45. Woon, D.E.; Dunning, T.H., Jr. Gaussian basis sets for use in correlated molecular calculations. III. The atoms aluminum through argon. J. Chem. Phys. 1993, 98, 1358–1371. [Google Scholar] [CrossRef]
  46. Wilson, A.K.; Woon, D.E.; Peterson, K.A.; Dunning, T.H., Jr. Gaussian basis sets for use in correlated molecular calculations. IX. The atoms gallium through krypton. J. Chem. Phys. 1999, 110, 7667–7676. [Google Scholar] [CrossRef]
  47. Figgen, D.; Rauhut, G.; Dolg, M.; Stoll, H. Energy-consistent pseudopotentials for group 11 and 12 atoms: Adjustment to multi-configuration Dirac–Hartree–Fock data. Chem. Phys. 2005, 311, 227–244. [Google Scholar] [CrossRef]
  48. Peterson, K.A.; Puzzarini, C. Systematically convergent basis sets for transition metals. II. Pseudopotential-based correlation consistent basis sets for the group 11 (Cu, Ag, Au) and 12 (Zn, Cd, Hg) elements. Theor. Chem. Acc. 2005, 114, 283–296. [Google Scholar] [CrossRef]
  49. Grimme, S. Semiempirical hybrid density functional with perturbative second-order correlation. J. Chem. Phys. 2006, 124, 034108. [Google Scholar] [CrossRef]
  50. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef]
  51. Cammi, R.; Tomasi, J. Remarks on the use of the apparent surface charges (ASC) methods in solvation problems: Iterative versus matrix-inversion procedures and the renormalization of the apparent charges. J. Comput. Chem. 1995, 16, 1449–1458. [Google Scholar] [CrossRef]
  52. Weinhold, F.; Glendening, E.D. NBO 5.0 Program Manual: Natural Bond Orbital Analysis Programs; Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin: Madison, WI, USA, 2001; p. 53706. [Google Scholar]
  53. Xiao, M.; Lu, T. Generalized charge decomposition analysis (GCDA) method. J. Adv. Phys. Chem. 2015, 4, 111–124. [Google Scholar] [CrossRef]
  54. Becke, A.D.; Edgecombe, K.E. A simple measure of electron localization in atomic and molecular systems. J. Chem. Phys. 1990, 92, 5397–5403. [Google Scholar] [CrossRef]
  55. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
Figure 1. Optimized structure of Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 at TPSSh/aug-cc-pVTZ level.
Figure 1. Optimized structure of Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 at TPSSh/aug-cc-pVTZ level.
Molecules 28 01508 g001
Figure 2. Molecular orbitals (MOs) and corresponding orbital number and type for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 (isosurface = 0.03 a.u.). The bonding orbitals are highlighted by red-dash rectangular frames.
Figure 2. Molecular orbitals (MOs) and corresponding orbital number and type for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 (isosurface = 0.03 a.u.). The bonding orbitals are highlighted by red-dash rectangular frames.
Molecules 28 01508 g002
Figure 3. MOs that contribute to σ-donation and their orbital composition from Au+ and ligands for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 (isosurface = 0.04 a.u.).
Figure 3. MOs that contribute to σ-donation and their orbital composition from Au+ and ligands for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 (isosurface = 0.04 a.u.).
Molecules 28 01508 g003
Figure 4. MOs that contribute to π-backdonation and their orbital composition from Au+ and ligands for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 (isosurface = 0.04 a.u.).
Figure 4. MOs that contribute to π-backdonation and their orbital composition from Au+ and ligands for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 (isosurface = 0.04 a.u.).
Molecules 28 01508 g004
Figure 5. The electron localization functions (ELFs) of Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 derived from TPSSh/aug-cc-pVTZ calculation.
Figure 5. The electron localization functions (ELFs) of Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 derived from TPSSh/aug-cc-pVTZ calculation.
Molecules 28 01508 g005
Table 1. Stability constants and standard reduction potentials for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 [1,5,14].
Table 1. Stability constants and standard reduction potentials for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 [1,5,14].
ComplexAu(CN)2Au(S2O3)23−Au[SC(NH2)2]2+Au(SCN)2AuCl2
logβi38.328.721.317.19.1
E0/V−0.570.170.350.661.11
Table 2. The geometry and structural parameters of Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 at TPSSh/aug-cc-pVTZ (a) and B2PLYP/def2-QZVP (b) levels. The lengths are given in Å and angles are given in degree.
Table 2. The geometry and structural parameters of Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2 at TPSSh/aug-cc-pVTZ (a) and B2PLYP/def2-QZVP (b) levels. The lengths are given in Å and angles are given in degree.
ComplexSymm.R(Au − X)(X = C, Cl, S)∠(X – Au − X)(X = C, Cl, S)
abab
Au(CN)2D∞h1.9871.987180.00180.00
Au(S2O3)23−C22.2992.297176.45176.86
Au[SC(NH2)2]2+C22.3042.299176.18176.15
Au(SCN)2C22.3102.305176.94177.16
AuCl2D∞h2.2912.287180.00180.00
Table 3. The atomic charges on Au and ligand atoms for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2.
Table 3. The atomic charges on Au and ligand atoms for Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2.
ComplexAtomCharge
Hirshfeld [28]ADCH [29]CM5 [30]CHELPG [31]
Au(CN)2Au0.0890.1000.2070.197
C−0.173−0.077−0.109−0.020
Au(S2O3)23−Au−0.0920.0310.0420.290
S−0.409−0.502−0.469−0.727
Au[SC(NH2)2]2+Au0.1350.1530.2710.427
S−0.178−0.224−0.247−0.438
Au(SCN)2Au0.0660.2670.1950.262
S−0.178−0.238−0.252−0.393
AuCl2Au−0.036−0.0040.0560.270
Cl−0.482−0.498−0.528−0.636
Table 4. Relative contributions to the coordinate bond and percentage contribution of the backdonation relative to the donation from CDA of Hartree–Fock wavefunctions at TPSSh/aug-cc-pVTZ level.
Table 4. Relative contributions to the coordinate bond and percentage contribution of the backdonation relative to the donation from CDA of Hartree–Fock wavefunctions at TPSSh/aug-cc-pVTZ level.
ComplexDonationBackdonationBackdonation/Donation (%)
Au(CN)20.80350.146418.22
Au(S2O3)23−0.52940.082915.66
Au[SC(NH2)2]2+0.44580.072216.20
Au(SCN)20.66730.093814.06
AuCl20.62260.04537.28
Table 5. Natural bond orbital (NBO) analysis and natural resonance theory (NRT) bond order (BO) of Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2.
Table 5. Natural bond orbital (NBO) analysis and natural resonance theory (NRT) bond order (BO) of Au(CN)2, Au(S2O3)23−, Au[SC(NH2)2]2+, Au(SCN)2, and AuCl2.
ComplexNBO AnalysisNRT BO
CovalentIonicCovalent/Total
Au(CN)225.37% (sd0.27)Au +
74.63% (sp0.87)C
0.2150.3120.408
Au(S2O3)23−24.01% (sp0.01d0.18)Au + 75.99% (sp5.82d0.04)S0.2150.2570.456
Au[SC(NH2)2]2+24.70% (sp0.01d0.18)Au + 75.30% (sp6.54d0.04)S0.2160.2650.449
Au(SCN)224.22% (sp0.01d0.18)Au + 75.78% (sp9.76d0.05)S0.1940.2530.434
AuCl219.41% (sp0.01d0.21)Au + 80.59% (sp6.61d0.02)Cl0.1760.8240.176
Table 6. Energy decomposition analyses for [XAuX] → AuX + X with the water model. The results are based on the TPSSh/aug-cc-pVTZ calculations. All energies are given in eV.
Table 6. Energy decomposition analyses for [XAuX] → AuX + X with the water model. The results are based on the TPSSh/aug-cc-pVTZ calculations. All energies are given in eV.
ComplexSteric
Interaction
Orbital
Interaction
TotalOrb./Total
Au(CN)21.32−4.09−2.771.477
Au(S2O3)23−1.49−3.64−2.141.696
Au[SC(NH2)2]2+1.46−3.43−1.971.740
Au(SCN)21.34−3.17−1.831.737
AuCl20.83−2.70−1.861.446
Table 7. Energy decomposition analyses for [XAuX] → X + Au+ X with the water model. The results are based on the TPSSh/aug-cc-pVTZ calculations. All energies are given in eV.
Table 7. Energy decomposition analyses for [XAuX] → X + Au+ X with the water model. The results are based on the TPSSh/aug-cc-pVTZ calculations. All energies are given in eV.
ComplexSteric
Interaction
Orbital
Interaction
TotalOrb./Total
Au(CN)23.41−9.37−5.971.571
Au(S2O3)23−2.68−7.84−5.161.520
Au[SC(NH2)2]2+3.17−7.43−4.261.744
Au(SCN)22.84−6.90−4.051.701
AuCl22.35−6.22−3.871.606
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, N.; Kou, J.; Sun, C. Investigation on Gold–Ligand Interaction for Complexes from Gold Leaching: A DFT Study. Molecules 2023, 28, 1508. https://doi.org/10.3390/molecules28031508

AMA Style

Zhang N, Kou J, Sun C. Investigation on Gold–Ligand Interaction for Complexes from Gold Leaching: A DFT Study. Molecules. 2023; 28(3):1508. https://doi.org/10.3390/molecules28031508

Chicago/Turabian Style

Zhang, Na, Jue Kou, and Chunbao Sun. 2023. "Investigation on Gold–Ligand Interaction for Complexes from Gold Leaching: A DFT Study" Molecules 28, no. 3: 1508. https://doi.org/10.3390/molecules28031508

APA Style

Zhang, N., Kou, J., & Sun, C. (2023). Investigation on Gold–Ligand Interaction for Complexes from Gold Leaching: A DFT Study. Molecules, 28(3), 1508. https://doi.org/10.3390/molecules28031508

Article Metrics

Back to TopTop