Next Article in Journal
2-Amino-N-Phenethylbenzamides for Irritable Bowel Syndrome Treatment
Previous Article in Journal
Preparation and Application of Co-Doped Zinc Oxide: A Review
Previous Article in Special Issue
Synthesis of Gold Clusters and Nanoparticles Using Cinnamon Extract—A Mechanism and Kinetics Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploration of Free Energy Surface of the Au10 Nanocluster at Finite Temperature

by
Francisco Eduardo Rojas-González
1,
César Castillo-Quevedo
2,
Peter Ludwig Rodríguez-Kessler
3,
José Oscar Carlos Jimenez-Halla
4,
Alejandro Vásquez-Espinal
5,
Rajagopal Dashinamoorthy Eithiraj
6,
Manuel Cortez-Valadez
7 and
José Luis Cabellos
8,*
1
Departamento de Física, Edificio 3F, Universidad de Sonora, Hermosillo 83000, Sonora, Mexico
2
Departamento de Fundamentos del Conocimiento, Centro Universitario del Norte, Universidad de Guadalajara, Carretera Federal No. 23, km. 191, Colotlán 46200, Jalisco, Mexico
3
Centro de Investigaciones en Óptica, A.C. (CIO) Lomas del Bosque 115, León 37150, Guanajuato, Mexico
4
Departamento de Química, División de Ciencias Exactas y Naturales, Universidad de Guanajuato, Noria Alta s/n, Guanajuato 36050, Guanajuato, Mexico
5
Química y Farmacia, Facultad de Ciencias de la Salud, Universidad Arturo Prat. Casilla 121, Iquique 1100000, Chile
6
Department of Physics, School of Advanced Sciences, Vellore Institute of Technology (VIT), Chennai 600 127, Tamil Nadu, India
7
CONAHCYT-Departamento de Investigación en Física, Universidad de Sonora, Apdo. Postal 5-88, Hermosillo 83190, Sonora, Mexico
8
Coordinación de Investigación y Desarrollo Tecnológico, Universidad Politécnica de Tapachula, Carretera Tapachula a Puerto Madero km. 24, Tapachula 30830, Chiapas, Mexico
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(14), 3374; https://doi.org/10.3390/molecules29143374
Submission received: 26 May 2024 / Revised: 8 July 2024 / Accepted: 15 July 2024 / Published: 18 July 2024
(This article belongs to the Special Issue Clusters—between Atoms and Nanoparticles)

Abstract

:
The first step in comprehending the properties of Au10 clusters is understanding the lowest energy structure at low and high temperatures. Functional materials operate at finite temperatures; however, energy computations employing density functional theory (DFT) methodology are typically carried out at zero temperature, leaving many properties unexplored. This study explored the potential and free energy surface of the neutral Au10 nanocluster at a finite temperature, employing a genetic algorithm coupled with DFT and nanothermodynamics. Furthermore, we computed the thermal population and infrared Boltzmann spectrum at a finite temperature and compared it with the validated experimental data. Moreover, we performed the chemical bonding analysis using the quantum theory of atoms in molecules (QTAIM) approach and the adaptive natural density partitioning method (AdNDP) to shed light on the bonding of Au atoms in the low-energy structures. In the calculations, we take into consideration the relativistic effects through the zero-order regular approximation (ZORA), the dispersion through Grimme’s dispersion with Becke–Johnson damping (D3BJ), and we employed nanothermodynamics to consider temperature contributions. Small Au clusters prefer the planar shape, and the transition from 2D to 3D could take place at atomic clusters consisting of ten atoms, which could be affected by temperature, relativistic effects, and dispersion. We analyzed the energetic ordering of structures calculated using DFT with ZORA and single-point energy calculation employing the DLPNO-CCSD(T) methodology. Our findings indicate that the planar lowest energy structure computed with DFT is not the lowest energy structure computed at the DLPN0-CCSD(T) level of theory. The computed thermal population indicates that the 2D elongated hexagon configuration strongly dominates at a temperature range of 50–800 K. Based on the thermal population, at a temperature of 100 K, the computed IR Boltzmann spectrum agrees with the experimental IR spectrum. The chemical bonding analysis on the lowest energy structure indicates that the cluster bond is due only to the electrons of the 6 s orbital, and the Au d orbitals do not participate in the bonding of this system.

1. Introduction

Nanoclusters of transition metals enhance catalytic activity due to their high surface-to-volume ratio and high surface energy [1,2]. Moreover, the transition metal atoms have incompletely filled d orbitals, allowing them to easily donate and accept electrons from other ions [3,4]. Gold is the most electronegative metallic atom [5], and when it forms bulk, gold is chemically inert [6,7,8,9]. However, when gold atoms make small clusters at nanoscale size, gold is not a noble, nonreactive material [10]. The Au clusters find applications in many fields of chemistry, physics, and nanomaterials [11]. Gold nanoclusters with diameters below 10 nm [6] are highly catalytic active due to quantum size effects [7,12], high oxidation states [13,14], and low coordinated atoms [15]., e.g., supported gold nanoclusters catalyze CO oxidation at low temperatures [8,16].
The lowest energy configurations of neutral Au clusters have been investigated from experimental and theoretical perspectives [11,17]. From a theoretical approach, there is a need to determine the cluster size at which the transition from 2D to 3D occurs [8,11,18], considering relativistic corrections, van der Waals forces, and temperature. This paper aims to provide insight in this direction.
In a study by Bravo-Pérez et al. [19], the neutral Aun (n = 3–7) clusters were investigated using density functional theory (DFT), and it was found that the lowest energy structure is planar. From an experimental standpoint, there is evidence that small Au clusters up to n = 7 are planar [20]. By using scalar-relativistic pseudopotentials for the valence electrons of gold and employing DFT, the lowest energy structure of Aun (n = 2–10) anionic and neutral clusters were located, revealing 3D structures [21]. Previous research has also looked into the lowest-energy structures of Aun (n = 2–20) clusters using a genetic algorithm coupled with DFT and empirical potential; the authors found that the low-energy structures are planar. Moreover, the 2D–3D transitions take place at Au atomic clusters with seven atoms [22]. For neutral Au clusters, previous high-level ab-initio calculations show that a structural transition from 2D planar to 3D structures occurs within the size range n = 7–10 [23,24,25,26,27]. Moreover, the transition could be between 12 and 15 atoms [28,29]. Assadollahzadeh et al. [30] reported a systematic search for the lowest energy structures of Aun (n = 2–20) clusters employing a genetic algorithm coupled to DFT with a relativistic pseudopotential. They found planar structures up to a cluster size of Au10. Sarosi et al. [31] reported a planar putative lowest energy structure employing DFT at the level of theory B3PW91/def2-TZVPD. Recently, Pham Vu Nhat et al. [11] suggested that the 2D–3D transitions for Au clusters occur at the size Au10 cluster. However, previous studies on the structure, stability, and bonding in Au10 clusters reported the putative global minimum as the three-dimensional structure at the MP2 level of theory [32]. The 2D–3D transition in neutral Aun clusters depends on the level of theory employed [31] and temperature. Furthermore, to establish a close-real value of energies, we suggest that energies must be computed at a single-point CCSD(T) level of theory [33,34].
The chemistry and physics of gold are strongly dominated by relativistic effects [35], e.g., the Au–Au bonding distance is shortened due to relativistic effects. Previous studies on the structure of small, neutral gold clusters employed the ZORA as a relativistic correction [10,36]. Romaniello et al. [37] included relativistic effects for calculating the linear response properties of metals, particularly the Au atom. Soto et al. [38] employed ZORA approximation to study the gold cluster cation Au7, and they found a putative lowest-energy puckered 3D structure. Due to relativistic effects, the preference for the planar shape of neutral gold clusters continues up to 11 atoms [27,30,39].
The van der Waals interactions (vdW) are critical for predicting the stability of molecules [40,41,42] and can change the energetic ordering of isomers [32,34,40,43,44]. Previous studies on small gold clusters in the size range of 12–14 atoms found that, considering van der Waals interactions in the energy calculations, the 2D–3D transition size is 13 atoms [45]. Besides relativistic effects and vdW interactions, the temperature is a critical factor in all molecular systems, mainly transition-metal nanoclusters exhibit catalytic activity at high temperatures. Gold nanoclusters, specifically, have recently received attention because of the discovery of their catalytic activity at low temperatures [10].
The experiments are performed at finite temperatures [10,46,47], and the materials also work at finite temperatures [32,34,44,47]. However, from the theoretical point of view, DFT calculations are typically performed at zero temperature and are valid at zero temperature and high vacuum conditions. However, using statistical mechanics, we can include the effects of temperature and pressure. In the mid-1960s, Mermin et al. [48] studied thermal properties in an inhomogeneous electron gas. Previous works have also addressed temperature in DFT [49,50]. More recently, DFT has been developed to account for finite temperatures [51,52,53]. Currently, temperature effects can be considered using two methods. The first is ab initio molecular dynamics, in which temperature is controlled by a thermostat [54,55]. Second, employing quantum statistical mechanics that requires the vibrational modes or phonon spectrum to compute the molecular partition function that contains all thermodynamic information [32,34,56,57,58]. A change in energetic separation among isomers distribution and dynamic structural rearrangements are the first effects of temperature [10,34,41,59]. Particularly in small gold clusters, the transitions from 2D to 3D structure depend on temperature. However, few studies have examined the gold isomers’ stability at finite temperatures [59]. Goldsmith et al. [41] studied neutral gas gold clusters at finite temperatures by replica-exchange ab initio molecular dynamics. A previous study on neutral Au12 predicts the dynamic coexistence of multiple planar and non-planar structures at room temperature [60]. Ghiringhelli et al. [10] studied the structure of small neutral gold clusters Au3, Au4, and Au7 and their infrared spectra. Previous work observed the IR spectra in small neutral gold clusters containing up to eight Au atoms and compared them with computed IR spectra [61,62]. In previous research on chemical bonding in Au clusters, Rodríguez et al. [63] used the QTAIM scheme to calculate the atomic properties and determine the chemical bonding at the Au-gold–thiol interface. Recently, Chebotaev et al. [64] studied the interactions between pterin and gold clusters using QTAIM theory. Additionally, Zubarev et al. [65] used the AdNDP method to study chemical bonding in Au20 cluster.
In this paper, we utilized a genetic algorithm coupled with DFT using the version 5.0 of the ORCA quantum chemistry program suite [66] to investigate the potential and free energy surfaces of the Au10 cluster. Our approach included local optimizations with the ZORA relativistic approximation and the dispersion D3BJ [67], considering thermal effects for temperatures ranging from 50 to 800 K based on quantum statistical mechanics. By employing a Boltzmann weights scheme [33,34,43,44], we computed the temperature-dependent Boltzmann–IR spectrum utilizing the calculated thermal population. Our study successfully validated that the computed Boltzmann–IR spectrum aligns with IR experimental data. Notably, we identified the 2D elongated hexagon configuration as the probable global minimum within the temperature range of 50–800 K, indicating that the transition from 2D to 3D does not occur in the Au10 cluster. Additionally, we rigorously examined the energetic ordering of low-energy structures computed at the DFT level against single-point energy calculations using the domain-based local pair natural orbital coupled cluster method with single-, double-, and perturbative triple-excitations (DLPNO-CCSD(T)) level of theory.
The remaining sections of this manuscript are organized: Section 2 presents an overview of the theory and computational details. In Section 3, we delve into the results covering various aspects such as the low-energy structures considering ZORA and DBJ3 dispersion, the chemical bonding with AdNDP and QTAIM, the energetic ordering of isomers computed at the DLPNO-CCSD(T) level of theory, the computed IR spectrum temperature dependent and compared with observed data, and the thermal population. These findings culminate in the conclusions presented in Section 4.

2. Results and Discussion

2.1. The Lowest Energy Structures

The final optimized geometries depicted in Figure 1 are reported at PBE0-D3BJ [67,68] and ZORA [10,36] SARC-ZORA-TZVP [69] SARC/J [70] and employing AutoAux generation procedure [71] level of theory. All calculations were performed using the version 5.0 of the ORCA quantum chemistry program suite [72].
Figure 1 shows the 12 most energetically stable structures of Au10. The relative Gibbs free energies are computed at 100 K, 1 atm, and lying within the energy range of 0.0–8.6 kcal/mol. The harmonic vibrational analysis of these structures showed no imaginary frequencies, corroborating that they are energy minima. The putative global minimum, 1, can be described as a planar elongated hexagon formed by triangular units with symmetry C2v. The thermal population is 88.7%, and the average Au–Au bond length is 2.7 Å. This value is slightly longer than the experimental one, 2.4715 Å, obtained from high-resolution rotational spectroscopy in the Au dimer [17,73] This system, 1, has been reported as the putative global minimum by different authors [31].
The second and third isomers lie around 0.57 kcal/mol above 1 in Gibbs free energy. These structures resemble a trigonal prism capped by three and one atoms at both ends. Their corresponding thermal populations are 4.93% and 4.76%, respectively.
The next series of isomers, 4, 5, and 6 are located at 0.99–1.06 kcal/mol above the putative global minimum. Their corresponding thermal population are 0.59, 0.49, and 0.41%. This form is quite similar to that found in low-lying Sn-type clusters [74]. This isomer is similar with what was previously reported as low-lying energy structures [11].
The two following isomers, 7 and 8, lie at 2.43 and 2.47 kcal/mol higher in energy than 1, and they are distorted capped trigonal prisms with a zero thermal population. Nhat et al. [11] used the PNO-LCCSD(T)-F12/aug-cc-pVTZ-PP level of theory, plus ZPE energy correction based on PBE/cc-pVTZ-PP-optimized geometries. They found the capped trigonal prism at 1.8 kcal/mol above the trigonal edge-capped prism. Structures 9 and 11 (3.51 and 7.27 kcal/mol above 1) have been reported as quasi-planar at 300 K [41]. The former is located at 3.3 kcal/mol at the PNO-LCCSD(T)-F12/aug-cc-pVTZ-PP level of theory [11]. The distorted structure, 10, at 5.89 kcal/mol, is a folded flake structure that has not been reported. The symmetric structure, 12, at 8.59 kcal/mol, is a small distorted cage.
To explore spin effects in determining the lowest energy structure of Au10, we computed the triplet state. All neutral Au10 optimized with spin multiplicities singlets were taken as initial structures for optimization without any symmetry restriction being performed considering triplet states. Our results suggest that the lowest-lying structure in the triplet state is 30 kcal/mol above the singlet putative global minimum. The second low-energy triplet is at 33 kcal/mol. Thus, higher spin states (up to triplet) are not energetically favored.

2.2. The Chemical Bonding Analysis on Au10 Cluster

2.2.1. AdNDP Analysis

To shed light in the bonding of the Au10 cluster, we carry out the chemical bonding analysis on the 2D elongated hexagon structure depicted in Figure 1 (1) and also on the 3D capped trigonal prism displayed in Figure 1 (6). The results of the AdNDP analysis for the lowest energy structure, displayed in Figure 1 (1), are shown in Figure 2. This analysis reveals the presence of five lone pairs (LPs) on each of the Au atoms with ON values very close to the optimum value of 2 |e|, which indicates that the Au d orbitals do not participate in the bonding of this system. The bonding pattern is completed by five delocalized sigma bond elements as follows: two 6c–2e bonds involving the hexagonal ends, two 7c–2e bonds involving the hexagonal ends plus the Au atom in the center of each hexagon, and, finally, a 10c–2e bond delocalized over the entire structure.
The analysis f or the 3D structure reveals the five LPs on each Au atom, while the remaining bonding pattern consists of three 3c–2e bonds, each of them involving the Au atom that is capping the vertical axis of the prism and the two atoms of the same axis, and two delocalized sigma bonds: one of 6c–2e and the other of 7c–2e, whose main contribution comes from the lower and upper triangular faces of the prism, respectively. Figure 3 displays the results of the AdNDP analysis for the 3D low-energy structure of the Au10 system.

2.2.2. QTAIM Analysis

The QTAIM parameters computed at the BCP provide information about the nature of the interaction between the two Au atoms [64].
The total number of critical points of each type we found for the lowest energy structure of the Au10 cluster is 10 NCP, 19 BCP, and 10 ring critical points. Figure 4a displays the BCP points, Figure 4b displays the index of BCP, and Figure 4c displays the BP paths. We found 39 CP points, which fulfill the Poincaré–Hopf relationship. According to our computed CP (Figure 4a), 19 BCPs indicate chemical bonds among Au–Au atoms. Additionally, Figure 4c shows the BP indicative of chemical bonding. To shed light on the chemical bonding Au–Au, in each BCP, we evaluated these QTAIM parameters: the electronic density (ρ), the Laplacian of electron density (∇2ρ), and the energy density H(r), that are displayed in Table 1. Analyzing the data of Table 1, the Au–Au average bond distance between Au–Au atoms is 2.70 Å, which agrees with the bond distance reported for gold nanoclusters, which are in the range of 2.6–2.84 Å [63]. From Table 1, the Laplacian of electron density is positive for all BCPs; the positive value of the Laplacian at the BCP points indicates a depletion of the density [63]. Thus, according to the values displayed in Table 1, the Au–Au bond is non-covalent, which agrees with previous work [63]. Similar weak non-covalent Pd–Pd interactions have been reported in the metallic (Pd–Pd) and heterobimetallic (Pd–Ir) clusters [75]. The positive Laplacian of electron density and negative energy density values indicate a closed-shell Au–Au interaction with a weak or partially covalent character.

2.3. Energetics

To shed light on the proper energy hierarchy of the Au10 isomers and compare them with DFT energies. Single-point energies of each isomer were computed employing the method of DLPNO-CCSD(T) compared with those obtained at the DFT level of theory. Previous studies pointed out that, employing different methods to compute energies yields different results [34], particularly in small gold clusters [28].
Table 2, in the first row, displays the number of isomers following the label cluster given in Figure 1. The second row of Table 2 displays the relative DFT Gibbs free energy computed at PBE0-D3BJ and ZORA SARC-ZORA-TZVP SARC/J level of theory and at 100 K. The third row of Table 2 displays single-point energy computed at the DLPNO-CCSD(T) level of theory plus ZPE energy. This follows the energetic ordering of the second row of Table 2 at the DFT level of theory. The second row of Table 2 indicates that the putative global minimum is the planar elongated hexagon structure depicted in Figure 1 (1), according to single-point energy calculation at the DLPNO-CCSD(T) level of theory, as displayed in the third row of Table 2. This isomer lies at 4.26 kcal/mol above the putative lowest energy structure. This structure is composed of a trigonal prism as the main central structure that is edge-capped by four atoms, as depicted in Figure 1 (6). Moreover, it is labeled as isomer number six in Table 2. In summary, Table 2 indicates that the lowest energy structure is the 2D elongated hexagon structure computed at the DFT level of theory. In contrast, at the DLPNO-CCSD(T) + ZPE level of theory, the 3D capped trigonal prism is the lowest energy structure. According to DLPNO-CCSD(T) SP energies, the second isomer, labeled number 4 in Table 2, is located at 0.03 kcal/mol above the putative global minimum. The next isomer, labeled 6 in Table 2, is the putative global minimum at the DLPNO-CCSD(T) level of theory. The computed thermal population for those three isomers indicates that they contribute 0.59, 0.49, and 0.41%, respectively. The following isomers, labeled from 7 to 12 in Table 2, exhibit a continuous increase in energy without exchanges and in agreement with DFT energy ordering.

3. Theoretical Methods and Computational Details

3.1. Lowest Energy Structure Search

The search for the lowest energy structure and the low-energy structures is a complicated task due to several factors: (a) the exploration of the potential energy surface should be systematic and unbiased [77]; (b) the number of possible combinations of atomic arrangements grows exponentially with the number of atoms leading to a combinatorial explosion problem [44,78]; (c) computation of the total energy requires quantum mechanical DLPNO-CCSD(T) level of theory to achieve a realistic value of energy; (d) sampling a large region of the configuration space [79]. Several methods and theoretical studies have been developed to search for the lowest energy structures. The design and use of algorithms, like simulated annealing [80,81,82,83,84,85], kick method [86,87], gradient embedded genetic algorithm (GEGA) [88,89,90], basin hopping [91,92], among others [93], helped to explore the potential energy surface. Previous work employed genetic algorithms [94,95,96,97,98] and kick methodology [97,99,100,101,102,103,104,105,106] coupled with density functional theory with the aim of exploring the potential energy surface. In this study, our computational procedure employs a genetic algorithm coupled to the Orca code and implemented in the GALGOSON code [34,43,107]. The methodology employs a three-step search strategy, the optimization in this first stage was at the PBE0 [68] and LANL2DZ [108] level of theory. In previous studies, PBE0 was employed to compute the density of states, which yielded good agreement with the observed spectrum for anion Au19 [109]. The LANL2DZ basis set is employed for transition metals due to its low computational cost [34,110], and as the second step, structures lying into ten kcal/mol found in the previous step were optimized employing the PBE0 functional and utilizing the SARC-ZORA-TZVP basis set [69], in conjunction with the auxiliary basis SARC/J [70] and using the AutoAux generation procedure [71]. The dispersion is taken into account by employing the atom-pairwise dispersion correction with the Becke–Johnson damping scheme (D3BJ) [67]. The relativistic effects are considered through SARC-ZORA-TZVP. Additionally, we ensure that each isomer is a proper local low-energy structure through the vibrational mode of each isomer, confirming that the lowest vibrational mode of each isomer is positive. As the third step, the single-point energy calculations at the DLPNO-CCSD(T) level were performed using the version 5.0 of the ORCA quantum chemistry program suite with TightPNO settings [72].

3.2. Thermochemistry Properties

The partition function describes all statistical thermodynamic properties of an ensemble of atoms at thermodynamic equilibrium. Here, we computed the thermodynamic properties employing the partition function Q displayed in Equation (1).
Q T = i g i e Δ E i / k B T
In Equation (1), g i is the degeneracy factor, k B is the Boltzmann constant, T is the absolute temperature, and Δ E i is the total electronic energy of a cluster [111,112]. Due to the coupling of the internal modes, the computation of Equation (1) could be complicated [43]. The starting point for computing the thermochemical data is the partition function given in Equation (1), and it can be decomposed as contributions from Δ E i rotational, translation, vibration, and electronic. Therefore, the total partition function of a molecule in its ground state is partition functions, so we can re-write Equation (1) as a product of partial partition functions [43,111,113] as given in Equation (2).
q = q t r a n s q r o t q v i b q e l e c
All Equations displayed in Table 3 are computed at a standard pressure of 1 atm and finite temperature, employing ideal gas statistical mechanics. The equations are equivalent to those given in Ref. [111] and any standard text of thermodynamics [112,113]. In this study, the q t r a n s is computed at a finite temperature and is employed to compute the translational entropy. The rotational contribution depends on moments of inertia and rotational symmetry number [33]; the so-called rigid rotor harmonic oscillator approximation is employed to separate the rotational and vibrational modes [43]. Previous works computed the difference in the rotational entropy computed with and without the symmetries [33]. To address the failure of the harmonic oscillator approximation at low frequencies, entropic contributions to the free energies are calculated by default using Grimme’s Quasi-RRHO approach [57]. The electronic partition function, q e l e c , is given by q e l e c = ω0. The contribution to the total energy of a molecule at a finite temperature is a sum of electronic, translational, vibrational, and rotational energies. We employ equations 3–6 to compute the internal energy (U), enthalpy (H), and Gibbs energy (G) of the Au10 cluster at finite temperature. The rotational contribution to the entropy is calculated using the expressions given by Herzberg [114].
U 0 = ε 0 + Z P E
U T = U 0 + E R o t + E T r a n s + E V i b + E e l e c t
H = U T + n R T
Δ G = Δ H Δ S T
In the Equations above, Z P E is the zero-point energy correction, ε 0 is the electronic energy, and E R o t + E T r a n s + E V i b + E e l e c t are the contributions to energy due to translation, rotation, electronic, and vibration as a function of temperature, respectively. In Equation (5), R is the ideal gas constant, n is the amount of substance, and T is the absolute temperature.

3.3. Thermal Population

In a molecule, the measured properties represent the statistical averages of the ensemble of geometric conformations or isomers available to the cluster [115] in an ensemble at thermal equilibrium. We computed the thermal population or probabilities of occurrence [34,56,79,111,116,117,118] employing Equation (7).
P T = e β Δ G K e β Δ G K
where β = 1 / k B T , k B is the Boltzmann constant, T is the temperature, and Δ G k is the Gibbs free energy of the kth isomer. Equation (7) is restricted so that the sum of all the thermal population at fixed temperature T , P i T , must be equal to 1, according to Equation (8).
i = 1 n P T = 1
In this study, I R B o l t is the Boltzmann-weighted spectrum at temperature T which is given by Equation (9).
I R B o l t = i = 1 n I R i x P i T
where n is the total number of clusters in the ensemble, and I R i is the IR of the ith isomer at temperature T = 0, and P i T is the probability of the ith isomer given by Equation (7).

3.4. The Chemical Bonding

The chemical bonding analysis was performed employing the adaptive natural density partitioning method [119], which can be seen as an extension of the natural bond orbital (NBO) [120,121,122] analysis in which both Lewis bonding (localized 1c–2e lone-pairs and 2c–2e bonds) and delocalized bonding elements are recovered from the analysis of the first-order reduced density matrix in the basis of the natural atomic orbitals (NAOs). The procedure consists of obtaining the transformation matrix between the basis set of atomic orbitals, the basis set of NAOs, and the first-order-reduced density matrix in the basis of these NAOs. Then, in the AdNDP analysis, a block diagonalization of this density matrix is performed, obtaining the corresponding eigenvectors (bonding elements) and their eigenvalues (occupation numbers, ON), which are selected if the ON is close enough to the threshold of 2.00 |e|. The transformation and first-order reduced density matrices were obtained at the PBE0/LANL2DZ level using the NBO7 program [123] interfaced with the Gaussian package. The AdNDP analysis was carried out with the freely MultiWFN program v3.8 [124].
Additionally, we used the quantum theory of atoms in molecules (QTAIM) [125,126] for determining the chemical bonding at the lowest energy structure of Au10 cluster employing the freely MultiWFN program v3.8 [124]. The electron density topology is based on the critical point (CP) concept. The CP is defined as the point where the first derivative of the electron density vanishes. The bond critical point (BCP) is found between two nuclear critical points (NCPs) and indicates the presence of a chemical bond.

4. Thermal Population

In this section, based on the computed free-energy differences among Au10 isomers, we calculated the thermal population [34,56,79,111,116,117,118] or Boltzmann distribution at finite temperature and as standard pressure, as given in Equation (7). The observed molecular properties are statistical averages over a collection of molecules. Consequently, the statistical averages are crucial for the experimental measurement of thermodynamic quantities [127]. Distinguishing isomers based on their energy differences is essential for determining the solid–solid transition point within the thermal population. The thermal population is related to several solid–solid transition points if free energy differences between the putative lowest energy structure and the first low-energy isomer are small. Conversely, if those free energy differences are considerably greater, the thermal population has no solid–solid transition points. Additionally, free-energy differences are related to the fundamental chemical quantities such as binding constants, solubilities, partition coefficients, and adsorption coefficients [128,129,130]. Figure 5a, displays the thermal population for each particular Au10 isomer for temperatures ranging from 50 to 800 K. We computed that the thermal population at the PBE0 Def2-TZVP def2/J RIJCOSX level of theory takes into account the D3BJ dispersion correction method. In Figure 5a, the thermal population of the lowest energy planar elongated hexagon configuration is depicted by a black solid line, and at temperatures ranging from 50 to 150 K, the thermal population is constant. At 150 K, it began to decay slowly exponentially up to room temperature, where it achieved 90% probability. Above room temperature, the thermal population decays almost linearly, achieving 50% at 800 K, as displayed in the upper panel of Figure 5a. The analysis of these results led to three interesting conclusions in temperatures ranging from 50 to 800 K. (1) There are no solid–solid transformation points, (2) the planar elongated hexagon configuration strongly dominates, and (3) all molecular properties of the Au10 cluster are attributed to planar-elongated hexagon configuration. The thermal population of the isomer, an edge-capped trigonal prism displayed in Figure 1 (4), is depicted by a red solid line in Figure 5a. At a temperature of 150 K, it starts to increase, achieving 5% probability at room temperature, and at a temperature of 800 K, it achieves almost 15% of probability. This low-energy structure does not contribute to the molecular properties. Figure 5b, displays the thermal population of Au10 computed in the temperature range of 50 to 800 K and at PBE0 ZORA SARC-ZORA-TZVP UseSym SARC/J AUTOAUX level of theory, taking into account the D3BJ dispersion correction method and relativistic effects through ZORA correction. The planar elongated hexagon structure strongly dominates in the temperatures range of 50–800 K. At 50 K, the probability of a planar elongated hexagon configuration is 98% and decays exponentially fast to achieve 58% probability at room temperature; above room temperature, the probability decays exponentially until it achieves 30% at 800 K. The solid–solid transformation point (Tss1) shown in Figure 5b, is located at 350 K, and occurs between isomers two and four located higher in energy above the putative global minimum, and other Tss points with lower probabilities less than 10%, may occurs around room temperature point. Notice that, in Figure 5a,b, the planar elongated hexagon configuration strongly dominates the thermal population. Moreover, relativistic effects do not significantly change the thermal population. At 50 K, the probability achieves almost 100%. Therefore, based on our calculations, the molecular properties of Au10 are mainly due to the putative global minimum energy structure, planar elongated hexagon configuration, in a temperature range of 50–800 K.

5. Infrared Spectra (IR) of Au10 Cluster

IR spectra are usually used to identify functional groups and bond chemical information; however, from the experimental point of view. Identifying which IR bands correspond to vibrational molecular modes can be complex and demands ab initio calculations. In these computations, temperature is often not considered, and discrepancies between experimental and theoretical IR spectra can arise from finite temperature and anharmonic effects. In this study, each spectrum of each isomer was computed using DFT, as it is implemented in the Gaussian 09 code. In Figure 6a–g, present the IR spectra for the presumed global minimum and the six lowest energy configurations, found within the relative energy range of 0–2.49 kcal/mol. The Figure 6h compares the computed Boltzmann spectrum of the Au10 cluster with the observed IR spectrum [11,41]. The Boltzmann weighted spectrum at 100 K is depicted in a black solid line. The red solid line represents the experimental infrared spectrum of Au10 taken from references [11,41]. The experimental spectrum in Figure 6h, shows two peaks located at 150 cm−1 and 180 cm−1, respectively. The spectrum is almost flat in the range of 50–150 cm−1. Interestingly, the computed IR Boltzmann spectrum is in good agreement with the experimental studies of Goldsmith et al. [11,41]. The peak of the IR Boltzmann spectrum, depicted with a black solid line located at 150 cm−1, shows good agreement with the experimental IR spectrum, as shown in Figure 6h.
However, the intensities of those peaks are similar; the Boltzmann spectrum is narrower than the experimental one. The vibrational mode of the peak is located at 150 cm−1, consisting of stretching and oscillating the two central atoms along the larger axis of the cluster and within the plane. Whereas the peak located at 180 cm−1 consists of the stretching and oscillation of the two central atoms along the shorter axis of the cluster and within the plane. Interestingly, 88.7% of the Boltzmann IR spectrum, depicted in the black solid line and displayed in Figure 6h, comes from the IR spectrum of the putative global minimum. The other 11.3% of the Boltzmann IR spectrum is composed of contributions from the six low-energy structures.

6. Conclusions

In summary, our study systematically explored the potential and free-energy surfaces of the Au10 cluster. We achieved this using an efficient genetic algorithm written in Python 2.7 coupled with density functional theory as implemented in the version 5.0 of the ORCA quantum chemistry program suite. The study highlights the findings of elucidate the lowest energy structure of the Au10 cluster, considering the relativistic corrections, dispersion, and temperature. We consider relativistic effects through the ZORA approximation, the dispersion through the D3BJ Grimme methodology, and the effects of temperature are included through nanothermodynamics.
The computed thermal population indicates that the planar elongated hexagon configuration strongly dominates in the temperature range from 0 to 800 K, regardless of whether ZORA is considered or not. Overall, relativistic corrections and temperature effects do not result in a transition from the 2D planar to a 3D Au10 structure. We also compared the energetic ordering of the low-energy structures computed at the DFT level of theory versus the single-point energies employing the DLPNO-CCSD(T) level of theory, which indicated an interchange in the energetic ordering of the isomers, determining a 3D structure as the lowest energy structure.
We computed the IR spectrum at 100 K employing the Boltzmann methodology. This agrees with experimental studies conducted by B. R. Goldsmith et al. The study highlights that the 2D elongated hexagon configuration is the putative global minimum in the temperature range from 0 to 800 K; consequently, the transition from 2D to 3D does not occur at the Au10 cluster.
The 2D elongated hexagon structure’s AdNDP analysis indicates that the bonding of the cluster is solely due to the electrons of the 6 s orbital, with no participation of the Au-d orbitals. The QTAIM analysis indicates that the positive Laplacian of electron density and negative energy density values (on BCP) indicate a closed-shell Au–Au interaction with a weak covalent character. For the capped trigonal prism-shaped isomer, the AdNDP analysis for this structure reveals the five LPs on each Au atom, while the remaining bonding pattern consists of three 3c–2e bonds. A forthcoming project involves computing the optical spectra of Au10 clusters using the Boltzmann scheme.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29143374/s1, The xyz atomic coordinates of the optimized Au10 cluster.

Author Contributions

F.E.R.-G. methodology, software, validation, performed the calculations, and drafted and wrote the manuscript; C.C.-Q. Methodology, software, validation, performed the calculations, and drafted the manuscript; J.L.C. conception of the study, software development, design, and validation, performed calculations and analysis, drafted the manuscript; performed the data analyses, investigation, resources, and revised and wrote the manuscript; J.O.C.J.-H. performed the calculations, analysis, revised, and wrote the manuscript; P.L.R.-K. performed the calculations, analysis, revised and wrote the manuscript; A.V.-E. performed the calculations, analysis, revised, and wrote the manuscript; R.D.E. performed the calculations, analysis, revised, and wrote the manuscript. M.C.-V. performed the calculations, analysis, revised, and wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received partial external funding. A.V.-E. thanks the financial support of the National Agency for Research and Development (ANID) through Project FONDECYT, Regular 1221019. F.E.R.-G. thanks CONAHCYT for the scholarship 492050.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Acknowledgments

We are grateful to. Ing. Jessica Vidal and L.C.C. Daniel Mendoza for granting us access to their clusters and computational support. Computational resources for this work were partially provided by ACARUS through the High-Performance Computing Area of the University of Sonora. We are grateful to the PAKAL supercomputer lab at the Polytechnic University of Tapachula for supporting this work.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

DFTDensity Functional Theory
DLPNO-CCSD(T)Domain-based Local Pair Natural Orbital Coupled-Cluster Theory
ZPEZero-Point Energy
IRVibrational Infrared Spectrum
BOFABoltzmann-Optics-Full-Ader code (Module Nanothermodynamics)
GALGOSONGlobal Genetic Algorithm

References

  1. Narayanan, R.; El-Sayed, M.A. Catalysis with Transition Metal Nanoparticles in Colloidal Solution Nanoparticle Shape Dependence and Stability. J. Phys. Chem. B 2005, 109, 12663–12676. [Google Scholar] [CrossRef] [PubMed]
  2. Khan, Y.; Sadia, H.; Ali Shah, S.Z.; Khan, M.N.; Shah, A.A.; Ullah, N.; Ullah, M.F.; Bibi, H.; Bafakeeh, O.T.; Khedher, N.B.; et al. Classification, Synthetic, and Characterization Approaches to Nanoparticles, and Their Applications in Various Fields of Nanotechnology: A Review. Catalysts 2022, 12, 1386. [Google Scholar] [CrossRef]
  3. Dutta, S.; Pati, S.K. Novel design of single transition metal atoms anchored on C6 N6 nanosheet for electrochemical and photochemical N2 reduction to Ammonia. Catal. Today 2023, 424, 113804. [Google Scholar] [CrossRef]
  4. Ma, X.; Hu, J.; Zheng, M.; Li, D.; Lv, H.; He, H.; Huang, C. N2 reduction using single transition-metal atom supported on defective WS2 monolayer as promising catalysts: A DFT study. Appl. Surf. Sci. 2019, 489, 684–692. [Google Scholar] [CrossRef]
  5. Kuhn, M.; Sammynaiken, R.; Sham, T. Electronic structure of Au-Ta alloys, An X-ray spectroscopy study. Phys. B Condens. Matter 1998, 252, 114–126. [Google Scholar] [CrossRef]
  6. Haruta, M. When Gold Is Not Noble: Catalysis by Nanoparticles. Chem. Rec. 2003, 3, 75–87. [Google Scholar] [CrossRef] [PubMed]
  7. Sun, K.; Kohyama, M.; Tanaka, S.; Takeda, S. Understanding of the activity difference between nanogold and bulk gold by relativistic effects. J. Energy Chem. 2015, 24, 485–489. [Google Scholar] [CrossRef]
  8. Yadav, J.; Saini, S. Atop adsorption of oxygen on small sized gold clusters Analysis of size and site reactivity from restructuring perspective. Comput. Theor. Chem. 2020, 1191, 113014. [Google Scholar] [CrossRef]
  9. Hammer, B.; Norskov, J.K. Why gold is the noblest of all the metals. Nature 1995, 376, 238–240. [Google Scholar] [CrossRef]
  10. Ghiringhelli, L.M.; Gruene, P.; Lyon, J.T.; Rayner, D.M.; Meijer, G.; Fielicke, A.; Scheffler, M. Not so loosely bound rare gas atoms: Finite-temperature vibrational fingerprints of neutral gold-cluster complexes. New J. Phys. 2013, 15, 083003. [Google Scholar] [CrossRef]
  11. Nhat, P.V.; Si, N.T.; Hang, N.T.N.; Nguyen, M.T. The lowest-energy structure of the gold cluster Au10: Planar vs. nonplanar? Phys. Chem. Chem. Phys. 2022, 24, 42–47. [Google Scholar] [CrossRef] [PubMed]
  12. Valden, M.; Lai, X.; Goodman, D.W. Onset of Catalytic Activity of Gold Clusters on Titania with the Appearance of Nonmetallic Properties. Science 1998, 281, 1647–1650. [Google Scholar] [CrossRef] [PubMed]
  13. Lopez, N.; Norskov, J.K. Catalytic CO Oxidation by a Gold Nanoparticle, A Density Functional Theory Study. J. Am. Chem. Soc. 2002, 124, 11262–11263. [Google Scholar] [CrossRef]
  14. Yang, G.; Wang, Y.; Peng, F.; Bergara, A.; Ma, Y. Gold as a 6p-Element in Dense Lithium Aurides. J. Am. Chem. Soc. 2016, 138, 4046–4052. [Google Scholar] [CrossRef] [PubMed]
  15. Boronat, M.; Leyva-Pérez, A.; Corma, A. Theoretical and Experimental Insights into the Origin of the Catalytic Activity of Subnanometric Gold Clusters Attempts to Predict Reactivity with Clusters and Nanoparticles of Gold. Acc. Chem. Res. 2014, 47, 834–844. [Google Scholar] [CrossRef] [PubMed]
  16. Schwank, J. Catalytic gold. Gold Bull. 1983, 16, 103–110. [Google Scholar] [CrossRef]
  17. Xiao, L.; Tollberg, B.; Hu, X.; Wang, L. Structural study of gold clusters. J. Chem. Phys. 2006, 124, 114309. [Google Scholar] [CrossRef] [PubMed]
  18. Liao, M.S.; Watts, J.D.; Huang, M.J. Theoretical Comparative Study of Oxygen Adsorption on Neutral and Anionic Agn and Aun Clusters (n = 2–25). J. Phys. Chem. C 2014, 118, 21911–21927. [Google Scholar] [CrossRef]
  19. Bravo-Pérez, G.; Garzón, I.; Novaro, O. Ab initio study of small gold clusters. J. Mol. Struct. THEOCHEM 1999, 493, 225–231. [Google Scholar] [CrossRef]
  20. Gilb, S.; Weis, P.; Furche, F.; Ahlrichs, R.; Kappes, M.M. Structures of small gold cluster cations Ion mobility measurements versus density functional calculations. J. Chem. Phys. 2002, 116, 4094–4101. [Google Scholar] [CrossRef]
  21. Häkkinen, H.; Landman, U. Gold clusters (AuN, 2<~ N <~10) and their anions. Phys. Rev. B 2000, 62, R2287–R2290. [Google Scholar] [CrossRef]
  22. Wang, J.; Wang, G.; Zhao, J. Density-functional study of Aun (n = 2−20) clusters: Lowest-energy structures and electronic properties. Phys. Rev. B 2002, 66, 035418. [Google Scholar] [CrossRef]
  23. Häkkinen, H.; Moseler, M.; Landman, U. Bonding in Cu, Ag, and Au Clusters: Relativistic Effects, Trends, and Surprises. Phys. Rev. Lett. 2002, 89, 033401. [Google Scholar] [CrossRef] [PubMed]
  24. Olson, R.M.; Varganov, S.; Gordon, M.S.; Metiu, H.; Chretien, S.; Piecuch, P.; Kowalski, K.; Kucharski, S.A.; Musial, M. Where Does the Planar-to-Nonplanar Turnover Occur in Small Gold Clusters. J. Am. Chem. Soc. 2005, 127, 1049–1052. [Google Scholar] [CrossRef] [PubMed]
  25. Han, Y.K. Structure of Au8: Planar or nonplanar? J. Chem. Phys. 2006, 124, 024316. [Google Scholar] [CrossRef] [PubMed]
  26. Bulusu, S.; Zeng, X.C. Structures and relative stability of neutral gold clusters, Aun, n = 15–19. J. Chem. Phys. 2006, 125, 154303. [Google Scholar] [CrossRef] [PubMed]
  27. Johansson, M.P.; Warnke, I.; Le, A.; Furche, F. At What Size Do Neutral Gold Clusters Turn Three-Dimensional. J. Phys. Chem. C 2014, 118, 29370–29377. [Google Scholar] [CrossRef]
  28. Nhat, P.V.; Si, N.T.; Anh, N.N.K.; Duong, L.V.; Nguyen, M.T. The Au12 Gold Cluster, Preference for a Non-Planar Structure. Symmetry 2022, 14, 1665. [Google Scholar] [CrossRef]
  29. Li, X.B.; Wang, H.Y.; Yang, X.D.; Zhu, Z.H.; Tang, Y.J. Size dependence of the structures and energetic and electronic properties of gold clusters. J. Chem. Phys. 2007, 126, 084505. [Google Scholar] [CrossRef] [PubMed]
  30. Assadollahzadeh, B.; Schwerdtfeger, P. A systematic search for minimum structures of small gold clusters, Aun n = 2–20. J. Chem. Phys. 2009, 131, 064306. [Google Scholar] [CrossRef]
  31. Sárosi, M.B.; Petrar, P.M.; King, R.B. Density functional study of bare gold clusters, the ten-vertex neutral system. J. Mol. Model. 2013, 19, 4585–4590. [Google Scholar] [CrossRef]
  32. David, J.; Guerra, D.; Restrepo, A. Structure, stability and bonding in the Au10 clusters. Chem. Phys. Lett. 2012, 539–540, 64–69. [Google Scholar] [CrossRef]
  33. Buelna-García, C.E.; Robles-Chaparro, E.; Parra-Arellano, T.; Quiroz-Castillo, J.M.; del Castillo-Castro, T.; Martínez-Guajardo, G.; Castillo-Quevedo, C.; de León-Flores, A.; Anzueto-Sánchez, G.; Martin-del Campo-Solis, M.F.; et al. Theoretical Prediction of Structures, Vibrational Circular Dichroism, and Infrared Spectra of Chiral Be4B8 Cluster at Different Temperatures. Molecules 2021, 26, 3953. [Google Scholar] [CrossRef] [PubMed]
  34. Castillo-Quevedo, C.; Buelna-Garcia, C.E.; Paredes-Sotelo, E.; Robles-Chaparro, E.; Zamora-Gonzalez, E.; Martin-del Campo-Solis, M.F.; Quiroz-Castillo, J.M.; del Castillo-Castro, T.; Martínez-Guajardo, G.; de Leon-Flores, A.; et al. Effects of Temperature on Enantiomerization Energy and Distribution of Isomers in the Chiral Cu13 Cluster. Molecules 2021, 26, 5710. [Google Scholar] [CrossRef] [PubMed]
  35. Schwerdtfeger, P. Relativistic effects in properties of gold. Heteroat. Chem. 2002, 13, 578–584. [Google Scholar] [CrossRef]
  36. van Lenthe, E.; Baerends, E.J.; Snijders, J.G. Relativistic total energy using regular approximations. J. Chem. Phys. 1994, 101, 9783–9792. [Google Scholar] [CrossRef]
  37. Romaniello, P.; de Boeij, P.L. The role of relativity in the optical response of gold within the time-dependent current-density-functional theory. J. Chem. Phys. 2005, 122, 164303. [Google Scholar] [CrossRef] [PubMed]
  38. Soto, J.R.; Molina, B.; Castro, J.J. Nonadiabatic structure instability of planar hexagonal gold cluster cation Au7 and its spectral signature. RSC Adv. 2014, 4, 8157–8164. [Google Scholar] [CrossRef]
  39. Pyykko, P. Relativistic effects in structural chemistry. Chem. Rev. 1988, 88, 563–594. [Google Scholar] [CrossRef]
  40. Luna-Valenzuela, A.; Cabellos, J.L.; Alonso, J.A.; Posada-Amarillas, A. Effects of van der Waals interactions on the structure and stability of Cu8-xPdx (x = 0, 4, 8) cluster isomers. Mater. Today Commun. 2021, 26, 102024. [Google Scholar] [CrossRef]
  41. Goldsmith, B.R.; Florian, J.; Liu, J.X.; Gruene, P.; Lyon, J.T.; Rayner, D.M.; Fielicke, A.; Scheffler, M.; Ghiringhelli, L.M. Two-to-three dimensional transition in neutral gold clusters: The crucial role of van der Waals interactions and temperature. Phys. Rev. Mater. 2019, 3, 016002. [Google Scholar] [CrossRef]
  42. Rydberg, H.; Dion, M.; Jacobson, N.; Schröder, E.; Hyldgaard, P.; Simak, S.I.; Langreth, D.C.; Lundqvist, B.I. Van der Waals Density Functional for Layered Structures. Phys. Rev. Lett. 2003, 91, 126402. [Google Scholar] [CrossRef] [PubMed]
  43. Buelna-Garcia, C.E.; Cabellos, J.L.; Quiroz-Castillo, J.M.; Martinez-Guajardo, G.; Castillo-Quevedo, C.; de Leon-Flores, A.; Anzueto-Sanchez, G.; Martin-del Campo-Solis, M.F. Exploration of Free Energy Surface and Thermal Effects on Relative Population and Infrared Spectrum of the Be6B11 Fluxional Cluster. Materials 2021, 14, 0112. [Google Scholar] [CrossRef] [PubMed]
  44. Buelna-García, C.E.; Castillo-Quevedo, C.; Quiroz-Castillo, J.M.; Paredes-Sotelo, E.; Cortez-Valadez, M.; Martin-del Campo-Solis, M.F.; López-Luke, T.; Utrilla-Vázquez, M.; Mendoza-Wilson, A.M.; Rodríguez-Kessler, P.L.; et al. Relative Populations and IR Spectra of Cu38 Cluster at Finite Temperature Based on DFT and Statistical Thermodynamics Calculations. Front. Chem. 2022, 10, 841964. [Google Scholar] [CrossRef] [PubMed]
  45. Kinaci, A.; Narayanan, B.; Sen, F.G.; Davis, M.J.; Gray, S.K.; Sankaranarayanan, S.K.R.S.; Chan, M.K.Y. Unraveling the Planar-Globular Transition in Gold Nanoclusters through Evolutionary Search. Sci. Rep. 2016, 6, 34974. [Google Scholar] [CrossRef]
  46. Guan, P.W.; Houchins, G.; Viswanathan, V. Uncertainty quantification of DFT-predicted finite temperature thermodynamic properties within the Debye model. J. Chem. Phys. 2019, 151, 244702. [Google Scholar] [CrossRef] [PubMed]
  47. Sutton, C.; Levchenko, S.V. First-Principles Atomistic Thermodynamics and Configurational Entropy. Front. Chem. 2020, 8, 757. [Google Scholar] [CrossRef] [PubMed]
  48. Mermin, N.D. Thermal Properties of the Inhomogeneous Electron Gas. Phys. Rev. 1965, 137, A1441–A1443. [Google Scholar] [CrossRef]
  49. Stoitsov, M.V.; Petkov, I.Z. Density functional theory at finite temperatures. Ann. Phys. 1988, 184, 121–147. [Google Scholar] [CrossRef]
  50. Pribram-Jones, A.; Pittalis, S.; Gross, E.K.U.; Burke, K. Thermal Density Functional Theory in Context. In Frontiers and Challenges in Warm Dense Matter; Graziani, F., Desjarlais, M.P., Redmer, R., Trickey, S.B., Eds.; Springer International Publishing: Cham, Switzerland, 2014; pp. 25–60. [Google Scholar]
  51. Gázquez, J.L.; Franco-Pérez, M.; Ayers, P.W.; Vela, A. Temperature-dependent approach to chemical reactivity concepts in density functional theory. Int. J. Quantum Chem. 2019, 119, e25797. [Google Scholar] [CrossRef]
  52. Gonis, A.; Däne, M. Extension of the Kohn-Sham formulation of density functional theory to finite temperature. J. Phys. Chem. Solids 2018, 116, 86–99. [Google Scholar] [CrossRef]
  53. Eschrig, H. T < 0 ensemble-state density functional theory via Legendre transform. Phys. Rev. B 2010, 82, 205120. [Google Scholar] [CrossRef]
  54. Jalife, S.; Liu, L.; Pan, S.; Cabellos, J.L.; Osorio, E.; Lu, C.; Heine, T.; Donald, K.J.; Merino, G. Dynamical behavior of boron clusters. Nanoscale 2016, 8, 17639–17644. [Google Scholar] [CrossRef]
  55. Martínez-Guajardo, G.; Luis Cabellos, J.; Díaz-Celaya, A.; Pan, S.; Islas, R.; Chattaraj, P.K.; Heine, T.; Merino, G. Dynamical behavior of Borospherene: A Nanobubble. Sci. Rep. 2015, 5, 11287. [Google Scholar] [CrossRef] [PubMed]
  56. Pracht, P.; Grimme, S. Calculation of absolute molecular entropies and heat capacities made simple. Chem. Sci. 2021, 12, 6551–6568. [Google Scholar] [CrossRef]
  57. Grimme, S. Supramolecular Binding Thermodynamics by Dispersion-Corrected Density Functional Theory. Chem.—A Eur. J. 2012, 18, 9955–9964. [Google Scholar] [CrossRef] [PubMed]
  58. Zaby, P.; Ingenmey, J.; Kirchner, B.; Grimme, S.; Ehlert, S. Calculation of improved enthalpy and entropy of vaporization by a modified partition function in quantum cluster equilibrium theory. J. Chem. Phys. 2021, 155, 104101. [Google Scholar] [CrossRef] [PubMed]
  59. Beret, E.C.; Ghiringhelli, L.M.; Scheffler, M. Free gold clusters: Beyond the static, monostructure description. Faraday Discuss. 2011, 152, 153–167. [Google Scholar] [CrossRef] [PubMed]
  60. Santarossa, G.; Vargas, A.; Iannuzzi, M.; Baiker, A. Free energy surface of two- and three-dimensional transitions of Au12 nanoclusters obtained by ab initio metadynamics. Phys. Rev. B 2010, 81, 174205. [Google Scholar] [CrossRef]
  61. Gruene, P.; Rayner, D.M.; Redlich, B.; van der Meer, A.F.G.; Lyon, J.T.; Meijer, G.; Fielicke, A. Structures of Neutral Au7, Au19, and Au20 Clusters in the Gas Phase. Science 2008, 321, 674–676. [Google Scholar] [CrossRef]
  62. Gruene, P.; Butschke, B.; Lyon, J.T.; Rayner, D.M.; Fielicke, A. Far-IR Spectra of Small Neutral Gold Clusters in the Gas Phase. Z. Für Phys. Chem. 2014, 228, 337–350. [Google Scholar] [CrossRef]
  63. Rodríguez, J.I.; Uribe, E.A.; Baltazar-Méndez, M.I.; Autschbach, J.; Castillo-Alvarado, F.; Gutiérrez-González, I. Size evolution relativistic DFT-QTAIM study on the gold cluster complexes Au4-S-CnH2n-S-Au4 (n = 2–5). Chem. Phys. Lett. 2016, 660, 287–294. [Google Scholar] [CrossRef]
  64. Chebotaev, P.P.; Plavskii, V.Y.; Kononov, A.I.; Buglak, A.A. Pterin interactions with gold clusters, A theoretical study. Dye. Pigment. 2023, 216, 111323. [Google Scholar] [CrossRef]
  65. Zubarev, D.Y.; Boldyrev, A.I. Deciphering Chemical Bonding in Golden Cages. J. Phys. Chem. A 2009, 113, 866–868. [Google Scholar] [CrossRef] [PubMed]
  66. Neese, F. Software update: The ORCA program system-Version 5.0. WIREs Comput. Mol. Sci. 2022, 12, e1606. [Google Scholar] [CrossRef]
  67. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef]
  68. Adamo, C.; Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys. 1999, 110, 6158–6170. [Google Scholar] [CrossRef]
  69. Pantazis, D.A.; Chen, X.Y.; Landis, C.R.; Neese, F. All-Electron Scalar Relativistic Basis Sets for Third-Row Transition Metal Atoms. J. Chem. Theory Comput. 2008, 4, 908–919. [Google Scholar] [CrossRef] [PubMed]
  70. Pantazis, D.A.; Neese, F. All-Electron Scalar Relativistic Basis Sets for the Lanthanides. J. Chem. Theory Comput. 2009, 5, 2229–2238. [Google Scholar] [CrossRef] [PubMed]
  71. Stoychev, G.L.; Auer, A.A.; Neese, F. Automatic Generation of Auxiliary Basis Sets. J. Chem. Theory Comput. 2017, 13, 554–562. [Google Scholar] [CrossRef]
  72. Liakos, D.G.; Guo, Y.; Neese, F. Comprehensive Benchmark Results for the Domain Based Local Pair Natural Orbital Coupled Cluster Method (DLPNO-CCSD(T)) for Closed- and Open-Shell Systems. J. Phys. Chem. A 2020, 124, 90–100. [Google Scholar] [CrossRef]
  73. Simard, B.; Hackett, P. High resolution study of the (0, 0) and (1, 1) bands of the A0u+-X0g+ system of Au2. J. Mol. Spectrosc. 1990, 142, 310–318. [Google Scholar] [CrossRef]
  74. Schafer, S.; Assadollahzadeh, B.; Mehring, M.; Schwerdtfeger, P.; Schäfer, R. Structure and Electric Properties of SnN Clusters (N = 6–20) from Combined Electric Deflection Experimentsand Quantum Theoretical Studies. J. Phys. Chem. A 2008, 112, 12312–12319. [Google Scholar] [CrossRef] [PubMed]
  75. 126 Lepetit, C.; Fau, P.; Fajerwerg, K.; Kahn, M.L.; Silvi, B. Topological analysis of the metal-metal bond: A tutorial review. Coord. Chem. Rev. 2017, 345, 150–181. [Google Scholar] [CrossRef]
  76. Grimme, S.; Muck-Lichtenfeld, C.; Erker, G.; Kehr, G.; Wang, H.; Beckers, H.; Willner, H. When Do Interacting Atoms Form a Chemical Bond? Spectroscopic Measurements and Theoretical Analyses of Dideuteriophenanthrene. Angew. Chem. Int. Ed. 2009, 48, 2592–2595. [Google Scholar] [CrossRef] [PubMed]
  77. Baletto, F.; Ferrando, R. Structural properties of nanoclusters: Energetic, thermodynamic, and kinetic effects. Rev. Mod. Phys. 2005, 77, 371–423. [Google Scholar] [CrossRef]
  78. Schuster, P. Taming combinatorial explosion. Proc. Natl. Acad. Sci. USA 2000, 97, 7678–7680. [Google Scholar] [CrossRef] [PubMed]
  79. Li, Z.H.; Jasper, A.W.; Truhlar, D.G. Structures, Rugged Energetic Landscapes, and Nanothermodynamics of Aln (2 ≤ n ≤ 65) Particles. J. Am. Chem. Soc. 2007, 129, 14899–14910. [Google Scholar] [CrossRef] [PubMed]
  80. Kirkpatrick, S.; Gelatt, C.D.; Vecchi, M.P. Optimization by Simulated Annealing. Science 1983, 220, 671–680. [Google Scholar] [CrossRef] [PubMed]
  81. Metropolis, N.; Rosenbluth, A.W.; Rosenbluth, M.N.; Teller, A.H.; Teller, E. Equation of State Calculations by Fast Computing Machines. J. Chem. Phys. 1953, 21, 1087–1092. [Google Scholar] [CrossRef]
  82. Xiang, Y.; Gong, X.G. Efficiency of generalized simulated annealing. Phys. Rev. E 2000, 62, 4473–4476. [Google Scholar] [CrossRef]
  83. Xiang, Y.; Gubian, S.; Suomela, B.; Hoeng, J. Generalized Simulated Annealing for Global Optimization: The GenSA Package for R. R J. 2013, 5, 13–29. [Google Scholar] [CrossRef]
  84. Vlachos, D.; Schmidt, L.; Aris, R. Comparison of small metal clusters: Ni, Pd, Pt, Cu, Ag, Au. Z. Phys. D At. Mol. Cl. 1993, 26, 156–158. [Google Scholar] [CrossRef]
  85. Granville, V.; Krivanek, M.; Rasson, J.P. Simulated annealing: A proof of convergence. IEEE Trans. Pattern Anal. Mach. Intell. 1994, 16, 652–656. [Google Scholar] [CrossRef]
  86. Saunders, M. Stochastic search for isomers on a quantum mechanical surface. J. Comput. Chem. 2004, 25, 621–626. [Google Scholar] [CrossRef] [PubMed]
  87. Saunders, M. Stochastic exploration of molecular mechanics energy surfaces. Hunt. Glob. Minimum. J. Am. Chem. Soc. 1987, 109, 3150–3152. [Google Scholar] [CrossRef]
  88. Alexandrova, A.N.; Boldyrev, A.I. Search for the Li0,+1,−1, (n = 5–7) Lowest-Energy Structures Using the ab-initio Gradient Embedded Genetic Algorithm (GEGA). Elucidation of the Chemical Bonding in the Lithium Clusters. J. Chem. Theory Comput. 2005, 1, 566–580. [Google Scholar] [CrossRef]
  89. Alexandrova, A.N.; Boldyrev, A.I.; Fu, Y.J.; Yang, X.; Wang, X.B.; Wang, L.S. Structure of the Na x Cl x + 1 (x = 1–4) clusters via ab-initio genetic algorithm and photoelectron spectroscopy. J. Chem. Phys. 2004, 121, 5709–5719. [Google Scholar] [CrossRef]
  90. Alexandrova, A.N. H·(H2O)n Clusters: Microsolvation of the Hydrogen Atom via Molecular ab-initio Gradient Embedded Genetic Algorithm (GEGA). J. Phys. Chem. A 2010, 114, 12591–12599. [Google Scholar] [CrossRef]
  91. Harding, D.; Mackenzie, S.R.; Walsh, T.R. Structural Isomers and Reactivity for Rh6 and Rh6+. J. Phys. Chem. B 2006, 110, 18272–18277. [Google Scholar] [CrossRef]
  92. Wales, D.J.; Doye, J.P.K. Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms. J. Phys. Chem. A 1997, 101, 5111–5116. [Google Scholar] [CrossRef]
  93. Balan, J.L.; Morato-Marquez, J.A.; Torres-Torres, J.G.; Cabellos, J.L.; Ortiz-Chi, F. Structural and optical properties of the Agn (n = 3–12); tyrosine complexes: A density functional theory study. R. Soc. Open Sci. 2023, 10, 230908. [Google Scholar] [CrossRef] [PubMed]
  94. Guo, J.C.; Feng, L.Y.; Wang, Y.J.; Jalife, S.; Vásquez-Espinal, A.; Cabellos, J.L.; Pan, S.; Merino, G.; Zhai, H.J. Coaxial Triple-Layered versus Helical Be6B11 Clusters: Dual Structural Fluxionality and Multifold Aromaticity. Angew. Chem. Int. Ed. 2017, 56, 10174–10177. [Google Scholar] [CrossRef] [PubMed]
  95. Dong, X.; Jalife, S.; Vásquez-Espinal, A.; Ravell, E.; Pan, S.; Cabellos, J.L.; Liang, W.Y.; Cui, Z.h.; Merino, G. Li2B12 and Li3B12: Prediction of the Smallest Tubular and Cage-like Boron Structures. Angew. Chem. Int. Ed. 2018, 57, 4627–4631. [Google Scholar] [CrossRef] [PubMed]
  96. Mondal, S.; Cabellos, J.L.; Pan, S.; Osorio, E.; Torres-Vega, J.J.; Tiznado, W.; Restrepo, A.; Merino, G. 10-π-Electron arenes a la carte: Structure and bonding of the [E-(Cn Hn)-E]n−6 (E = Ca, Sr, Ba; n = 6–8) complexes. Phys. Chem. Chem. Phys. 2016, 18, 11909–11918. [Google Scholar] [CrossRef]
  97. Ravell, E.; Jalife, S.; Barroso, J.; Orozco-Ic, M.; Hernandez-Juarez, G.; Ortiz-Chi, F.; Pan, S.; Cabellos, J.L.; Merino, G. Structure and Bonding in CE5(E = Al–Tl) Clusters: Planar Tetracoordinate Carbon versus Pentacoordinate Carbon. Chem.—Asian J. 2018, 13, 1467–1473. [Google Scholar] [CrossRef] [PubMed]
  98. Grande-Aztatzi, R.; Martínez-Alanis, P.R.; Cabellos, J.L.; Osorio, E.; Martínez, A.; Merino, G. Structural evolution of small gold clusters doped by one and two boron atoms. J. Comput. Chem. 2014, 35, 2288–2296. [Google Scholar] [CrossRef] [PubMed]
  99. Pan, S.; Moreno, D.; Cabellos, J.L.; Romero, J.; Reyes, A.; Merino, G.; Chattaraj, P.K. In Quest of Strong Be-Ng Bonds among the Neutral Ng-Be Complexes. J. Phys. Chem. A 2014, 118, 487–494. [Google Scholar] [CrossRef]
  100. Cui, Z.H.; Ding, Y.H.; Cabellos, J.L.; Osorio, E.; Islas, R.; Restrepo, A.; Merino, G. Planar tetracoordinate carbons with a double bond in CAl3E clusters. Phys. Chem. Chem. Phys. 2015, 17, 8769–8775. [Google Scholar] [CrossRef] [PubMed]
  101. Vargas-Caamal, A.; Pan, S.; Ortiz-Chi, F.; Cabellos, J.L.; Boto, R.A.; Contreras-Garcia, J.; Restrepo, A.; Chattaraj, P.K.; Merino, G. How strong are the metallocene–metallocene interactions? Cases of ferrocene, ruthenocene, and osmocene. Phys. Chem. Chem. Phys. 2016, 18, 550–556. [Google Scholar] [CrossRef]
  102. Vargas-Caamal, A.; Cabellos, J.L.; Ortiz-Chi, F.; Rzepa, H.S.; Restrepo, A.; Merino, G. How Many Water Molecules Does it Take to Dissociate HCl? Chem.—A Eur. J. 2016, 22, 2812–2818. [Google Scholar] [CrossRef]
  103. Cui, Z.h.; Vassilev-Galindo, V.; Luis Cabellos, J.; Osorio, E.; Orozco, M.; Pan, S.; Ding, Y.h.; Merino, G. Planar pentacoordinate carbon atoms embedded in a metallocene framework. Chem. Commun. 2017, 53, 138–141. [Google Scholar] [CrossRef] [PubMed]
  104. Vargas-Caamal, A.; Ortiz-Chi, F.; Moreno, D.; Restrepo, A.; Merino, G.; Cabellos, J.L. The rich and complex potential energy surface of the ethanol dimer. Theor. Chem. Acc. 2015, 134, 16. [Google Scholar] [CrossRef]
  105. Flórez, E.; Acelas, N.; Ibargüen, C.; Mondal, S.; Cabellos, J.L.; Merino, G.; Restrepo, A. Microsolvation of NO3: Structural exploration and bonding analysis. RSC Adv. 2016, 6, 71913–71923. [Google Scholar] [CrossRef]
  106. Castro, A.C.; Osorio, E.; Cabellos, J.L.; Cerpa, E.; Matito, E.; Solà, M.; Swart, M.; Merino, G. Exploring the Potential Energy Surface of E2 P4 ClustersEGroup13 Element the Quest for Inverse Carbon-Free Sandwiches. Chem.—A Eur. J. 2014, 20, 4583–4590. [Google Scholar] [CrossRef] [PubMed]
  107. Buelna-Garcia, C.E.; Castillo-Quevedo, C.; Paredes-Sotelo, E.; Martinez-Guajardo, G.; Cabellos, J.L. Boltzmann Populations of the Fluxional Be6B11 and Chiral Be4B8 Clusters at Finite Temperatures Computed by DFT and Statistical Thermodynamics. In Density Functional Theory; Glossman-Mitnik, D., Ed.; IntechOpen: Rijeka, Croatia, 2021; Chapter 4. [Google Scholar] [CrossRef]
  108. Hay, P.J.; Wadt, W.R. Ab-initio effective core potentials for molecular calculations. Potentials for K to Au including the outermost core orbitals. J. Chem. Phys. 1985, 82, 299–310. [Google Scholar] [CrossRef]
  109. Shao, N.; Huang, W.; Gao, Y.; Wang, L.M.; Li, X.; Wang, L.S.; Zeng, X.C. Probing the Structural Evolution of Medium-Sized Gold Clusters Aun (n = 27–35). J. Am. Chem. Soc. 2010, 132, 6596–6605. [Google Scholar] [CrossRef]
  110. Hay, P.J.; Wadt, W.R. Ab-initio effective core potentials for molecular calculations. Potentials for the transition metal atoms Sc to Hg. J. Chem. Phys. 1985, 82, 270–283. [Google Scholar] [CrossRef]
  111. Dzib, E.; Cabellos, J.L.; Ortíz-Chi, F.; Pan, S.; Galano, A.; Merino, G. Eyringpy: A program for computing rate constants in the gas phase and in solution. Int. J. Quantum Chem. 2019, 119, e25686. [Google Scholar] [CrossRef]
  112. McQuarrie, D.A.M. Statistical Mechanics; Chemistry Series; Harper & Row: New York, NY, USA, 1975. [Google Scholar]
  113. Hill, T. An Introduction to Statistical Thermodynamics; Addison-Wesley Series in Chemistry; Dover Publications: New York, NY, USA, 1986. [Google Scholar]
  114. Herzberg, G. Infrared and Raman Spectra of Polyatomic Molecules. In Molecular Spectra and Molecular Structure; Van Nostrand: New York, NY, USA, 1945. [Google Scholar]
  115. Teague, S.J. Implications of protein flexibility for drug discovery. Nat. Rev. Drug Discov. 2003, 2, 527–541. [Google Scholar] [CrossRef]
  116. Mendoza-Wilson, A.M.; Balandrán-Quintana, R.R.; Cabellos, J.L. Thermochemical behavior of sorghum procyanidin trimers with C4-C8 and C4-C6 interflavan bonds in the reaction with superoxide anion radical and H2O2-forming NADH-oxidase flavoenzyme. Comput. Theor. Chem. 2020, 1186, 112912. [Google Scholar] [CrossRef]
  117. Schebarchov, D.; Baletto, F.; Wales, D.J. Structure, thermodynamics, and rearrangement mechanisms in gold clusters-insights from the energy landscapes framework. Nanoscale 2018, 10, 2004–2016. [Google Scholar] [CrossRef] [PubMed]
  118. Grigoryan, V.G.; Springborg, M. Temperature and isomeric effects in nanoclusters. Phys. Chem. Chem. Phys. 2019, 21, 5646–5654. [Google Scholar] [CrossRef] [PubMed]
  119. Zubarev, D.Y.; Boldyrev, A.I. ”Developing paradigms of chemical bonding: Adaptive natural density partitioning. Phys. Chem. Chem. Phys. 2008, 10, 5207–5217. [Google Scholar] [CrossRef] [PubMed]
  120. Foster, J.P.; Weinhold, F. Natural hybrid orbitals. J. Am. Chem. Soc. 1980, 102, 7211–7218. [Google Scholar] [CrossRef]
  121. Reed, A.E.; Weinhold, F. Natural bond orbital analysis of near-Hartree-Fock water dimer. J. Chem. Phys. 1983, 78, 4066–4073. [Google Scholar] [CrossRef]
  122. Reed, A.E.; Weinhold, F. Natural localized molecular orbitals. J. Chem. Phys. 1985, 83, 1736–1740. [Google Scholar] [CrossRef]
  123. Glendening, E.D.; Badenhoop, J.K.; Reed, A.E.; Carpenter, J.E.; Bohmann, J.A.; Morales, C.M.; Karafiloglou, P.; Landis, C.R.; Weinhold, F. NBO 7.0. Theoretical Chemistry Institute, University of Wisconsin, Madison 2018. Available online: https://nbo.chem.wisc.edu/biblio_css.html (accessed on 15 March 2024).
  124. Lu, T.; Chen, F. Multiwfn, A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  125. Bader, R.F.W. A quantum theory of molecular structure and its applications. Chem. Rev. 1991, 91, 893–928. [Google Scholar] [CrossRef]
  126. Kumar, P.S.V.; Raghavendra, V.; Subramanian, V. Bader’s Theory of Atoms in Molecules (AIM) and its Applications to Chemical Bonding. J. Chem. Sci. 2016, 128, 1527–1536. [Google Scholar] [CrossRef]
  127. Palma, C.A.; Björk, J.; Klappenberger, F.; Arras, E.; Kühne, D.; Stafström, S.; Barth, J.V. Visualization and thermodynamic encoding of single-molecule partition function projections. Nat. Commun. 2015, 6, 6210. [Google Scholar] [CrossRef]
  128. van Gunsteren, W.; Daura, X.; Mark, A. Computation of Free Energy. Helv. Chim. Acta 2002, 85, 3113–3129. [Google Scholar] [CrossRef]
  129. Hansen, N.; van Gunsteren, W.F. Practical Aspects of Free-Energy Calculations: A Review. J. Chem. Theory Comput. 2014, 10, 2632–2647. [Google Scholar] [CrossRef] [PubMed]
  130. Christ, C.D.; van Gunsteren, W.F. Enveloping distribution sampling: A method to calculate free energy differences from a single simulation. J. Chem. Phys. 2007, 126, 184110. [Google Scholar] [CrossRef] [PubMed]
  131. van Lenthe, E.; Snijders, J.G.; Baerends, E.J. The zero-order regular approximation for relativistic effects, The effect of spin-orbit coupling in closed shell molecules. J. Chem. Phys. 1996, 105, 6505–6516. [Google Scholar] [CrossRef]
Figure 1. The optimized geometries of the neutral Au10 cluster at PBE0 ZORA SARC-ZORA-TZVP SARC/J D3BJ level of theory. The most important low-energy isomers appear in two orientations. The relative Gibbs free energies (in kcal/mol) are shown in parentheses, the point group symmetry is displayed in square brackets, while the thermal population at 100 K is given in red parentheses. The yellow spheres stand for the Au atoms (atomic coordinates available in Supplementary Materials).
Figure 1. The optimized geometries of the neutral Au10 cluster at PBE0 ZORA SARC-ZORA-TZVP SARC/J D3BJ level of theory. The most important low-energy isomers appear in two orientations. The relative Gibbs free energies (in kcal/mol) are shown in parentheses, the point group symmetry is displayed in square brackets, while the thermal population at 100 K is given in red parentheses. The yellow spheres stand for the Au atoms (atomic coordinates available in Supplementary Materials).
Molecules 29 03374 g001
Figure 2. (Color online) Results of the AdNDP analysis for the lowest energy structure of the Au10 system. Each green and white isosurface represents the N-center orbital at each vertex of the cluster, respectively.
Figure 2. (Color online) Results of the AdNDP analysis for the lowest energy structure of the Au10 system. Each green and white isosurface represents the N-center orbital at each vertex of the cluster, respectively.
Molecules 29 03374 g002
Figure 3. Shows the results of the AdNDP analysis for the 3D capped trigonal prism-shaped isomer displayed in Figure 1 (6). Each green and white isosurface represents the N-center orbital at each vertex of the cluster, respectively.
Figure 3. Shows the results of the AdNDP analysis for the 3D capped trigonal prism-shaped isomer displayed in Figure 1 (6). Each green and white isosurface represents the N-center orbital at each vertex of the cluster, respectively.
Molecules 29 03374 g003
Figure 4. (Color online) Molecular graphs for the lowest energy structure Au10 cluster. (a) The BCPs are depicted in red-yellow color small spheres; (b) Index of the BCPs (numbers); (c) Bond paths (BPs) are the lines of maximum density that connects a pair of nuclear critical points. The BPs are used as an indicator of chemical bonding [76] between a pair of Au atoms. The atoms are depicted as red-brown spheres.
Figure 4. (Color online) Molecular graphs for the lowest energy structure Au10 cluster. (a) The BCPs are depicted in red-yellow color small spheres; (b) Index of the BCPs (numbers); (c) Bond paths (BPs) are the lines of maximum density that connects a pair of nuclear critical points. The BPs are used as an indicator of chemical bonding [76] between a pair of Au atoms. The atoms are depicted as red-brown spheres.
Molecules 29 03374 g004
Figure 5. (Color online) The thermal population for all Au10 isomers at temperatures ranging from 50 to 800 K. For an easy comparison, in panel (a), we show the thermal population without relativistic effects. In panel (b), we show the thermal population where relativistic effects are taken through zero-order regular approximation (ZORA) [131].
Figure 5. (Color online) The thermal population for all Au10 isomers at temperatures ranging from 50 to 800 K. For an easy comparison, in panel (a), we show the thermal population without relativistic effects. In panel (b), we show the thermal population where relativistic effects are taken through zero-order regular approximation (ZORA) [131].
Molecules 29 03374 g005
Figure 6. (Color online) At 0 K, the IR spectra of Au10 isomers were computed in the frequency range of 50–250 cm−1. Panel (ag) display the IR spectra of each low–energy isomer computed with Version 5.0 of the ORCA quantum chemistry program suite. Panel (h) compares the IR experimental, and the Boltzmann spectra computed at 100 K and depicted as a black solid line. The red solid line represents the experimental IR spectrum from references [11,41].
Figure 6. (Color online) At 0 K, the IR spectra of Au10 isomers were computed in the frequency range of 50–250 cm−1. Panel (ag) display the IR spectra of each low–energy isomer computed with Version 5.0 of the ORCA quantum chemistry program suite. Panel (h) compares the IR experimental, and the Boltzmann spectra computed at 100 K and depicted as a black solid line. The red solid line represents the experimental IR spectrum from references [11,41].
Molecules 29 03374 g006
Table 1. (1) Indexes of BCP as they are displayed in Figure 4b. (2) Electronic density (ρ). (3) Laplacian of electron density (∇2ρ). (4) Energy density H(r). (5) Bond distances.
Table 1. (1) Indexes of BCP as they are displayed in Figure 4b. (2) Electronic density (ρ). (3) Laplacian of electron density (∇2ρ). (4) Energy density H(r). (5) Bond distances.
Indexes BCP ρ2ρH(r)Bond Distance
140.57480.1145−0.18102.6892
190.72270.1291−0.24362.6770
230.57480.1145−0.18102.6892
240.53580.1100−0.14782.7173
200.71300.8501−0.22872.7304
150.53580.1100−0.14782.7173
Table 2. The first row of the table displays the number of isomers following the ordering in Figure 1. The second row shows the relative DFT Gibbs free energy computed at the PBE0-D3BJ and ZORA SARC-ZORA-TZVP SARC/J level of theory, and at 100 K. The third row shows the relative single point energy (SPE) computed employing DLPNO-CCSD(T) level of theory, with the TightPNO setting and ZPE energy correction based on DFT.
Table 2. The first row of the table displays the number of isomers following the ordering in Figure 1. The second row shows the relative DFT Gibbs free energy computed at the PBE0-D3BJ and ZORA SARC-ZORA-TZVP SARC/J level of theory, and at 100 K. The third row shows the relative single point energy (SPE) computed employing DLPNO-CCSD(T) level of theory, with the TightPNO setting and ZPE energy correction based on DFT.
Energy DFT versus DLPNO-CCSD(T) SPE Energy
No. of Isomer
(Figure 1)
DFT Energy
(kcal/mol)
DLPNO-CCSD(T) SPE
Energy (kcal/mol)
10.04.26
20.578.44
30.589.34
40.990.03
51.030.53
61.060.0
72.433.18
82.473.73
93.513.92
105.894.56
117.277.11
128.5915.08
Table 3. The contribution of electronic, translational, vibrational, and rotational partition functions to Equation (1).
Table 3. The contribution of electronic, translational, vibrational, and rotational partition functions to Equation (1).
ContributionPartition Function
Translational q trans = 2 π m k B T h 2 3 2 k B T P
Rotational linear q rot l = 1 σ r T Θ r
Rotational nonlinear q r o t n l = π Θ A Θ B Θ C 1 / 2 T 3 / 2 σ r , Θ i = 2 2 I i k B , i = A , B , C
Vibrational q vib = i = 1 n v e Θ vib , i / 2 T 1 e Θ vib , i / T , Θ vib , i = h ν i k B
Electronic q e = ω 0
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

Rojas-González, F.E.; Castillo-Quevedo, C.; Rodríguez-Kessler, P.L.; Jimenez-Halla, J.O.C.; Vásquez-Espinal, A.; Eithiraj, R.D.; Cortez-Valadez, M.; Cabellos, J.L. Exploration of Free Energy Surface of the Au10 Nanocluster at Finite Temperature. Molecules 2024, 29, 3374. https://doi.org/10.3390/molecules29143374

AMA Style

Rojas-González FE, Castillo-Quevedo C, Rodríguez-Kessler PL, Jimenez-Halla JOC, Vásquez-Espinal A, Eithiraj RD, Cortez-Valadez M, Cabellos JL. Exploration of Free Energy Surface of the Au10 Nanocluster at Finite Temperature. Molecules. 2024; 29(14):3374. https://doi.org/10.3390/molecules29143374

Chicago/Turabian Style

Rojas-González, Francisco Eduardo, César Castillo-Quevedo, Peter Ludwig Rodríguez-Kessler, José Oscar Carlos Jimenez-Halla, Alejandro Vásquez-Espinal, Rajagopal Dashinamoorthy Eithiraj, Manuel Cortez-Valadez, and José Luis Cabellos. 2024. "Exploration of Free Energy Surface of the Au10 Nanocluster at Finite Temperature" Molecules 29, no. 14: 3374. https://doi.org/10.3390/molecules29143374

APA Style

Rojas-González, F. E., Castillo-Quevedo, C., Rodríguez-Kessler, P. L., Jimenez-Halla, J. O. C., Vásquez-Espinal, A., Eithiraj, R. D., Cortez-Valadez, M., & Cabellos, J. L. (2024). Exploration of Free Energy Surface of the Au10 Nanocluster at Finite Temperature. Molecules, 29(14), 3374. https://doi.org/10.3390/molecules29143374

Article Metrics

Back to TopTop