Fernando
Buendía
a,
Jorge A.
Vargas
a,
Marcela R.
Beltrán
*a,
Jack B. A.
Davis
b and
Roy L.
Johnston
b
aInstituto de Investigaciones en Materiales, Universidad Nacional Autónoma de México, Circ. ext. s/n Apdo. Postal 70-360, C.P. 04510, México D.F., Mexico. E-mail: mbeltran@unam.mx; Fax: +52 (55)56161251; Tel: +52 (55)56224624
bSchool of Chemistry, University of Birmingham, Birmingham, West Midlands B15 2TT, UK
First published on 13th July 2016
A comparative theoretical study has been performed of the gas phase and deposited AumRhn (4 ≤ m + n ≤ 6) clusters. The combined use of a genetic algorithm and Density Functional Theory (DFT) calculations allows us to explore the potential energy surface and, therefore, find efficiently and automatically the global minimum configuration for each composition. Our results show interesting effects on the geometries of the clusters on deposition. This occurs because the rhodium atoms (electronically) prefer to be in contact with the MgO surface, sometimes promoting planar clusters to become three-dimensional when deposited, and three-dimensional clusters in the gas phase to become two-dimensional. Together with the change in geometries, the magnetic moment is reduced from the gas phase, as the electrons rearrange themselves when the cluster interacts with the substrate.
In the present work we study binary gold rhodium clusters. Despite the fact that there have been very few studies,27,28 this alloy is found to possess very exciting properties. The presence of high magnetic moments and a very high reactivity has been recently reported,26,29 in comparison to pure gold clusters.30–33 From both experimental and theoretical points of view, there are advantages of studying gas-phase clusters, because of their interesting properties such as size selectivity, ionic species, and adsorption. The interplay between the electronic and geometric structures can be clearly established. In fact, gas phase experiments have already had a marked impact on the understanding of cluster structures.2–9 The deposition of clusters on a surface is, nonetheless, an important advancement in the development of future technologies based on more stable subnanoscale devices.
This study is based on the fact that many experimental and theoretical studies of small particles show important changes in both their electronic as well as their geometric structures once they are deposited.34–40 In this work, we use the well-studied and electronically stable MgO(100) surface.41–48 Its stability offers the opportunity to study the effects of the cluster geometries once deposited on the surface. Keeping the system as small as possible is very important for reducing computing time, but it is also relevant to the increased interest in catalysis by small size-selected clusters.49 The introduction of a surface increases the computational complexity of the optimization process, due to the increase in the number of atoms involved.
Additionally, when the gas-phase minima are deposited on the surface and then locally minimized (on the surface) the lowest energy surface-supported clusters are not necessarily found because of the presence of energy barriers to cluster rearrangement. Here, we use the Birmingham Parallel Genetic Algorithm (BPGA) which combined with the Vienna ab initio Simulation Package (VASP) can accurately predict the structure of clusters in the gas phase.50–53 The method can also be applied to finding the true lowest energy geometries for surface-supported clusters, optimised by the genetic algorithm in presence of the surface.52,54
In the next section, we briefly discuss the details of both the DFT and the BPGA methods used in this work. In Section 3, we present and discuss a comparative study between the structural, electronic and magnetic properties for AumRhn (4 ≤ m + n ≤ 6) clusters in the gas phase versus supported (where the global minimum search is performed over the MgO surface). Finally, we present our conclusions.
Within the BPGA, the application of crossover and mutation operators to the given geometries and their local relaxation is managed by independently working subprocesses synchronized with a global database (pool). The initial population is generated through the minimization (with VASP) of random structures. Subsequently, the crossover and mutation operations begin. For the crossover, a pair of clusters are taken from the pool through “roulette-wheel” selection.57 Offspring are then produced through single-point, weighted crossover, carried out according to the Deaven and Ho cut and splice method.65 For the mutation, there are two possibilities, one is moving 20% of the cluster atoms, and the other (possible only in bimetallic clusters) is swapping two atoms of different elements generating a homotop. In our case, 90% of the new candidates were obtained via crossover, while 10% were determined by applying the move mutation operator to an older structure from the pool.
The BPGA can work in two modes to find the global minimum structure: in the gas phase or by taking into account the interaction with a surface. We used both for this comparative study. It is important to differentiate between the two methods to simulate the clusters on a surface. This can be done by placing onto the surface the lowest-energy clusters previously obtained in the gas phase and allowing the combined system to relax. What BPGA does, however, is to globally optimise the cluster directly on the surface, which predicts the most stable cluster (of given size and composition) that can be grown on the substrate or the result of allowing a gas phase deposited cluster to anneal on the surface. This yields structural transformations of the deposited cluster, otherwise not found using the first method mentioned above. This has been corroborated in an initial stage of this study.
The initial random pool in the gas phase consists of 10 configurations, while for the supported clusters the pool has 15 initial configurations. The calculations have been performed within a cubic supercell leaving at least 10 Å of space between periodically repeated clusters to avoid cluster–cluster interactions. In the second part, the clusters have been placed 1.5 Å above a 6 × 6 × 2 slab of an ideal MgO(100) surface leaving 14.7 Å of vacuum between slabs in the z-direction. We consider only two MgO layers after comparing both the total energies and the geometrical reconfigurations varying the number of layers and observing small differences (<2%) which are not important for the purpose of this study.
The binding energies have been calculated from the following expression,
![]() | (1) |
The energies, Eads, of the surface-supported global minima were calculated as follows,
Eads = Eslab+AumRhn − (Eslab + EAumRhn), | (2) |
![]() | ||
Fig. 1 Geometry structures for (a) gas phase atomic neutral Au4−nRhn (n = 0–4) clusters and their corresponding side (b) and top (c) views when supported on MgO. |
![]() | ||
Fig. 2 Geometry structures for (a) gas-phase atomic neutral Au5−nRhn (n = 0–5) clusters and their corresponding side (b) and top (c) views when supported on MgO. |
![]() | ||
Fig. 3 Geometry structures for (a) gas-phase atomic neutral Au6−nRhn (n = 0–6) clusters and their corresponding side (b) and top (c) views when supported on MgO. |
For the 4-atom systems, we obtained a consistent decrease of the magnetic moment, both with the number of gold atoms and also in the presence of the surface. The binding energies increase with the number of rhodium atoms, due to the greater strength of Rh–Rh binding as evidenced by the higher cohesive energy of bulk Rh (5.75 eV per atom) compared to Au (3.81 eV per atom).66 The most stable cluster geometries on the substrate are often different from those optimised in the gas phase, because there is a tendency to increase the number of Rh atoms in contact with the surface as Rh–O interactions are stronger than Au–O. The clusters tend to change in geometry as they land on the surface depending on the availability of rhodium atoms to bind to the surface.
Minimal structural changes are observed for Rh4 between the gas phase (1a) and supported (1b) clusters, with only small differences in the bond lengths. The adsorption energy is −2.16 eV (Fig. 4a) and the magnetic moment of the cluster is quenched (see Fig. 5a). AuRh3 has a planar rhombic structure in the gas phase that turns into a tetrahedral structure when it is supported. Three rhodium atoms bond to the surface, giving rise to the highest adsorption energy (−2.69 eV) of all the supported 4-atom clusters. Its magnetic moment is also reduced by 2 μB when compared with the gas phase. The most drastic geometrical changes and a decrease of 2 μB in its magnetic moment happen for Au2Rh2. This phenomenon is due to the stronger interaction of the surface with rhodium rather than with gold atoms. Au3Rh has the same structure whether in the gas phase or supported, similar to what happens to pure Rh4, and has the lowest adsorption energy (−1.75 eV), but the multiplicity remains the same for the gas phase and the supported cluster. Finally Au4 has a 2D rhombohedral structure and one Au–Au bond breaks in order for stronger bonds to be made with the superficial (or uppermost) oxygen atoms, yielding an adsorption energy of −1.90 eV. The two systems (gas phase and supported) are both singlets.
![]() | ||
Fig. 4 Adsorption energies on eV for (a) Au4−nRhn (n = 0–4), (b) Au5−nRhn (n = 0–5) and (c) Au6−nRhn (n = 0–6) clusters. |
![]() | ||
Fig. 5 Total cluster magnetic moment (μB) for (a) Au4−nRhn (n = 0–4), (b) Au5−nRhn (n = 0–5) and (c) Au6−nRhn (n = 0–6) clusters. |
As a reference, we have considered the global minima in the gas phase and then relaxed them on the MgO surface. The geometries show little change, other than small changes in bond lengths. However, the magnetic moments are greatly affected by the interaction with the surface, diminishing in all cases in a similar way as in the clusters generated on the surface by means of the BPGA. The differences in energy between the surface deposited and BPGA surface optimised structures are the largest for AuRh3 (0.41 eV) and Au2Rh2 (0.33 eV), which also display the largest structural differences to the gas phase structures.
In the 5-atom systems, there are many differences between the geometrical structures of the gas phase and supported clusters. In the gas phase, the clusters are three dimensional (3D) apart from Au4Rh, which is planar. But when they are supported, the lowest energy configuration is 2D, with the exception of Rh5, which retains the 3D gas phase structure. Incidentally this cluster has the highest adsorption energy to the surface (−3.38 eV, see Fig. 4b). It also maintains its high multiplicity (M = 5). The ground state of AuRh4 exhibits structural changes when optimised on the surface; it transforms from an edge-bridged tetrahedron to an edge-bridged square (Au forms the bridge). Its adsorption energy is the lowest for these alloys; this may be due to the geometric transformation (see Fig. 4b). The magnetic moment of the gas-phase minimum is reduced from 7 to 5 μB when it is supported. Au2Rh3 and Au3Rh2 also present conformational changes when supported. The gas-phase clusters have 3D structures with a magnetic moment of 5 μB, for both clusters. However, the lowest minima for the supported clusters are W-like 2D structures, standing on their edge perpendicular to the surface. Their magnetic moment is the lowest (1 μB) of all the 5-atom clusters. The observed tendency for the planar structures standing perpendicular to the MgO surface is due to the “metal-on-top” effect, commonly known for Au clusters.52 Au4Rh and Au5 are very similar, both in their geometrical and electronic structures in the gas phase. They have 2D structures and only have one unpaired electron. Their adsorption energies are −2.11 and −1.52 eV, respectively (see Fig. 4b). Again if we only deposit the gas-phase minima on the MgO(100) surface, they do not undergo any significant geometry change, though their magnetic moment decreases.
In the 6-atom systems, we observe a very interesting behaviour both in the gas phase and when supported. Their geometry is dominated by the number of gold atoms, being 2D for high gold concentrations. The transition to 3D occurs at Au2Rh4 in the gas phase and at Au3Rh3 when supported on MgO. Their adsorption energies decrease with the number of gold atoms (see Fig. 4c). The magnetic moment for the mixed clusters generated on the surface is constant (2 μB see Fig. 5c).
The ground state structures for Rh6, both in the gas phase and when deposited, are trigonal prisms, they also have 6 unpaired electrons (see Fig. 5c), and the adsorption energy (−3.40 eV), corresponding to adsorption via one of the square faces of the trigonal prism, is the highest among all the systems studied here. The AuRh5 cluster has the same structure in the gas phase and when supported, with 4 and 2 unpaired electrons, respectively, and an adsorption energy of −2.86 eV. Au2Rh4 is slightly changed when supported, with an adsorption energy of −2.51 eV. In this case a Rh–Rh bond is broken when supported to allow two rhodium atoms to form strong bonds with the substrate and to stabilise the structure via the metal-on-top effect. The magnetic moment is significantly quenched relative to the free cluster (from 6 to 2 μB as can be seen in Fig. 5c).
Au3Rh3 is unique in the sense that its structure changes from planar in the gas phase to a 3D structure when deposited. This is possibly due to the huge interaction between the oxygen anions in the surface and the rhodium atoms. This cluster has the lowest adsorption energy of all the 6-atom systems with a value of −1.96 eV (as can be seen in Fig. 4c), and its multiplicity changes from a septet in the gas phase to a triplet for the supported cluster. In the gas phase, Au4Rh2 has a triangular structure with the rhodium atoms in central positions and has a multiplicity of five. The ground state of the supported cluster has a similar structure and lies perpendicular to the surface; the rhodium atoms change their position to the basis of the triangle to be in contact with the surface (see Fig. 3b; due to a combination of the metal-on-top effect and the stronger Rh–O bonding), and the arrangement of the rhodium atoms is different (i.e. a different homotop is selected), with 2 Rh atoms lying along the lower edge of the triangle, enabling 2 strong Rh–O bonds to be formed. A similar rearrangement has been observed for Au2Ir4/MgO,52 and the magnetic moment decreases from 4 to 2 μB. The adsorption energy in this case is −2.06 eV (see Fig. 4c).
Finally, Au5Rh and Au6 have the same triangular structure. In the case of Au5Rh the Rh atom occupies the centre of the triangular edge which is bonded to the MgO surface (see Fig. 3b) and they have low magnetic moments both in the gas phase and when supported; see Fig. 5c. The adsorption energies of Au5Rh and Au6 are −1.97 eV and −1.59 eV, respectively. Again these clusters sit perpendicular to the surface to reduce the total energy due to the metal-on-top effect.
As for the smaller clusters, depositing gas phase minima on the MgO(100) surface leads to only minor structural changes and on average these structures are 0.23 eV higher in energy than the surface-optimised clusters generated by BPGA.
If we compare the clusters optimised on the surface by BPGA with those where the gas-phase minima are deposited onto the surface and locally minimized, we observe that the deposited clusters are less stable (i.e. the combined cluster + substrate energy is higher) than those which are optimised on the surface. The geometries, and also the spin multiplicities of these deposited gas phase clusters, differ from the substrate-optimised clusters discussed here.
We can foresee very important applications for clusters exhibiting the so-called metal-on-top effect. This could be very useful for catalytic purposes as their exposed superficial area is maximized. Another interesting result to emphasize is the fact that some pure or doped rhodium clusters may keep their magnetic moments if they do not suffer a drastic structural transformation upon deposition making them very useful for future magnetic applications.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c6cp03735f |
This journal is © the Owner Societies 2016 |