Electro-Design of Bimetallic PdTe Electrocatalyst for Ethanol Oxidation: Combined Experimental Approach and Ab Initio Density Functional Theory (DFT)—Based Study

An alternative electrosynthesis of PdTe, using the electrochemical atomic layer deposition (E-ALD) method, is reported. The cyclic voltammetry technique was used to analyze Au substrate in copper (Cu2+), and a tellurous (Te4+) solution was used to identify UPDs and set the E-ALD cycle program. Results obtained using atomic force microscopy (AFM) and scanning electron microscopy (SEM) techniques reveal the nanometer-sized flat morphology of the systems, indicating the epitaxial characteristics of Pd and PdTe nanofilms. The effect of the Pd:Te ratio on the crystalline structure, electronic properties, and magnetic properties was investigated using a combination of density functional theory (DFT) and X-ray diffraction techniques. Te-containing electrocatalysts showed improved peak current response and negative onset potential toward ethanol oxidation (5 mA; −0.49 V) than Pd (2.0 mA; −0.3 V). Moreover, DFT ab initio calculation results obtained when the effect of Te content on oxygen adsorption was studied revealed that the d-band center shifted relative to the Fermi level: −1.83 eV, −1.98 eV, and −2.14 eV for Pd, Pd3Te, and Pd3Te2, respectively. The results signify the weakening of the CO-like species and the improvement in the PdTe catalytic activity. Thus, the electronic and geometric effects are the descriptors of Pd3Te2 activity. The results suggest that Pd2Te2 is a potential candidate electrocatalyst that can be used for the fabrication of ethanol fuel cells.


Introduction
Electrodeposition of thin films and nanostructured metal chalcogenides plays an important part in technological development. Nanostructured metal chalcogenides have been extensively studied, as they are characterized by fascinating properties that are different from the properties of the bulk counterparts. Transition metal chalcogenide-based tellurides have attracted immense interest in the field of catalysis, energy conversion, storage [1], and photo-electrocatalysis, as they exhibit high electrical conductivity and can be used in the field of molecular structure catalysis [2][3][4]. Among electrocatalysis reactions, ethanol (EtOH) oxidation reaction (EOR) has been a subject of interest, as a high theoretical density (8.01 kWhg −1 ) can be realized [5]. Among the direct alcohol fuel cells (DAFCs)

Reagents and Substrate Pretreatment
The precursor solutions for Te and TePd were prepared using deionized water (purity: >18.2 µScm −1 ). Methods previously reported by us were used for sample preparation [11].

Electrochemical Characterization
The electrochemical characterization method was conducted to analyze the deposition potentials and electrocatalyst performance. AutoLab potentiostat (Metrohm-PGSTAT 302, Metrohm Autolab B.V. Kanaalweg 29-G 3526 KM Utrecht, The Netherlands) was used for analysis. The reference electrode was the saturated Ag|AgCl|3 M KCl, the working electrode was gold (Au) coated on glass (2 × 2 cm 2 ), and the counter electrode was an Au wire embedded in PLEXIGLAS®. The performance of the electrode catalyst toward EtOH solution containing 0.5 M KOH and 0.1 M EtOH was determined using the CV technique (the voltage was cycled between −1.0 and +0.2 V). Prior to experimental analysis, the EtOH solution was de-aerated using N 2 (g). CA and EIS techniques were used to examine the stability of the system and study the electron kinetics.

Physical Characterization
The X-ray diffraction (XRD) technique was used for analysis. The X-ray Rigaku Ultima IV X-Ray DI powder system was used to record the profiles. Cu Kα was used as the source to analyze the crystal structure and phase of the samples. The structural morphology of the samples was studied using a scanning electron microscope (SEM). The morphology and surface topography were examined using the atomic force microscopy (AFM) technique (Agilent Technologies AC Mode III instrument, Santa Clara, CA, USA), ambient tapping mode, using the RTESPW tip (Veeco Manufacturing, Inc., Plainview, NY, USA) model.

Cyclic Voltammetry of Au Coated Glass Recorded in HClO 4 Solution Te 4+ and Cu 2+ Ions
The CV technique was used to examine the UPD chemistry range of the metal. The E-ALD program used to deposit Te, Pd, and TePd on Au was studied, and the probable deposition mechanism was revealed. The CVs of the HTeO 2 + and Cu 2+ precursor solutions were recorded. Figure 1a shows the CVs of Au (in Te 4+ solution), and the initial potential was 1.0 V. The final potential was varied during cycling (−0.2 V, −0.4 V, and −0.6 V). The reduction potential of the tellurous solution proceeds over multiple steps to produce telluride, a soluble (H 2 Te) species. The Te reduction-oxidation reaction is expressed as follows Equation (1): (1) Figure 1a presents four progressive reductive characteristics observed along the negative scan direction. The peaks at +380 mV and −400 mV can be attributed to Te-UPD , and bulk Te (0) is believed to deposit at approximately −500 mV. The peak at −580 mV is attributed to the onset of hydrogen evolution (H 2 -evol), and intense H 2 -evol is observed at −770 mV. Bulk Te-strip contributes to the oxidative current peaks appearing during cathodic scans at +450 mV and +600 mV. These peaks can be attributed to the oxidation of Te-UPD1 and Te-UPD2 , respectively. The redox characteristic corroborates tellurium deposition in acidic electrolytes, as illustrated by the reduction reaction mechanism reported in the literature [20]. Conclusively, the preference for using Te (OPD = −600 mV) can be attributed to the slow kinetics of Te-UPD . Figure 1b presents data on the Au substrate in solution of 1 mM CuSO 4 + 0.1 M HClO 4 (cycled from +0.6 V to −0.2 V). The presence of two prominent peaks (at 0.180 V and −0.099 V) corresponding to the system under study has been reported in the literature [B]. The peak corresponding to Cu-UPD is more intense than the peak corresponding to Cu( 0 ) counterpart, and this can be attributed to the strong metal substrate (M-S) interactions (between Cu-UPD and the metal substrate) [21]. This indicates that the strength of the bond formed between Cu-UPD and the Au substrate is stronger than the strength of the metalmetal (M-M) bond. The results were recorded around the overpotential deposition (OPD) potential [22]. The standard reduction potential and Cu underpotential deposition are well presented in literature reports [11]. The standard reduction potential corresponding to Cu has been presented in the Supplementary Materials Equations (S1)-(S3). The shapes of the CVs and the peaks observed were almost the same as those reported previously [13,[18][19][20][21][22][23][24][25][26].

Electrodeposition of Monometallic (Te, Pd) and Bimetallic (PdTe, in Various Te Ratios) Nanofilms
Previously reported electrodeposition methods were followed for the electrodeposition of thin films [11,13,27]. The process of depositing Te on Au has been reported in the literature [11], and the schematic representation of the deposition process is presented in Figure 2a,b. The TePd E-ALD formation cycle proceeds over several steps: (i) Cu = open-circuit potential (OCP), 10 s; (ii) Cu-UPD = +180 mV, 10 s; (iii) Pd = OCP, 10 s, (iv) the BE = OCP, (v) Te = OCP, 15 s, (vi) held Step 4 quiescent to permit Te bulk deposition of at −500 mV, and (vii) BE = −600 mV (allows the reduction of bulk Te keeping the atomic Te layer intact). The last step is associated with oxidative Te underpotential deposition (oxTe-UPD ). This process was repeated over Te cycles to generate different Pd:Te ratios (Pd 3 Te 1 , 10 Te cycles; Pd 3 Te 2 , 15 Te cycles).  Figure 2 presents the method of formation of Te and PdTe. It can be observed that when the Te solution is flushed at OCP, no current spike is generated. When Te is deposited at OPD (−500 mV), a prominent reductive current spike is generated, which indicates the occurrence of the process of bulk Te deposition, as denoted in Figure 2a. When the potential for Be flow is set at −600 mV, an oxidative current spike is generated. It was believed that small amounts of bulk Te are stripped off under these conditions, and Te adatoms are left behind. These steps corroborate the results obtained by analyzing Equation (1), and this explains the redox couple observed for oxidative Te-UPD and the process of re-reduction of bulk Te(0). The formation of bulk and soluble tellurite (H 2 Te) species is presented in Figure 2b. Figure 2 also presents the SLRRR process occurring at the beginning of the cycles. Similar results have been previously reported by our group [11,28]. A negative (reductive) current shift was observed during the Cu-UPD step at 180 mV. This reflects the deposition of the Cu-atomic layer. The SLRRR process corresponding to the interchange of Cu-UPD for Pd at the open-circuit potential (OCP) has been reported previously by us [11]. The results are consistent with previously reported results [11,28,29]. The E-ALD model was developed following the methods reported in the literature [11]. A distinct change in current between −1000 and 0.1 µA is observed when BE flow is observed at −500 mV. The spike in the peak is believed to correspond to oxidative current response. This indicates the presence of oxidative Te-UPD adatoms on the surface of Pd. The OCP values increase with an increase in the number of cycles associated with Te deposition on Pd. This is revealed by the Labview sequencer 4 software. This potentially indicates that the process of Te deposition results in varying degrees of surface coverage. Deposition occurs on Pd and Te-modified Pd surfaces.  [3]. It was observed that Pd 3 Te 1 consists of two phases (white and gray particles). The generation of the two phases could be attributed to the characteristics of the Pd and Te nanostructures. The white particles show the characteristics of typical Te-dendritic burls. These particles are cauliflower-shaped, and similar results were reported by Ma et al. [30]. Interestingly, Pd 3 Te 2 deposits present flower-shaped (fern fronds) profiles that are present in abundance throughout the surface. The microstructures observed for these samples were similar to the Te-feather-like microstructures reported previously [31]. The dendrites grow from the core (nuclei) in different directions. These consist of distinct branches of different sizes. The characteristics agree well with the characteristics of the deposited dendritic Te particles reported in the literature [32]. Notably, the branches increase in size as the number of Te cycles increases (15 cycles) (Figure 3c). The increase in the sizes of the branches resulted in the formation of a large number of irregular trunks that overlap with each other at angles in the range of 45-60 • . The ends of the sub-branches appear cauliflower-like. These contain distinct layers, the formation of which can be attributed to the layer-by-layer Te-deposition method. Conclusively, Te plays a crucial role in alternating the structural morphology of the Pd surface. This indicates that the surface characteristics depend on the number of ox-Te-UPD cycles. The characteristics are influenced by the surface-directing diffusion process. Furthermore, the incorporation of Te into Pd was confirmed by the morphology distortion process that resulted in the generation of dendritic structures, the formation of rugged surfaces, and the generation of grain boundaries. Figure 3d-f reveals that the Pd samples present flat and uniform surfaces. The surface characteristics can be potentially attributed to the layer-by-layer epitaxial 2D growth process promoted by Pd SLRR. Pd 3 Te 1 grains sit on top of each other, and this indicates early nucleation and 3D growth. These processes result in the formation of an island. Te-containing Pd presents a rough surface, and coalesced particles attributable to the nanostructures of neighboring Te and Pd units were observed. Grain boundaries were observed for the case of Pd 3 Te 2 , and the formation of grain boundaries can be attributed to the deposition of Te at distinct active sites (Pd remnant site and Te-covered Pd sites). The effect of the Te ratio on the morphology was studied, and it was observed that the roughness increases with an increase in the ratio. This can be attributed to the deposition of the overpotential deposition (OPD) steps associated with Te. The nanofilms differ in roughness. The average surface roughness (Sa) is calculated for Pd (Sa = 21.59 nm), Pd 3 Te 1 (Sa = 32.50 nm), and Pd 3 Te 2 (Sa = 45.5 nm). The results reflect the SEM results presented in Figure 3a-c.

Electrochemical Performance of Pd-Based Electrode Catalysts and Their Kinetics during EOR
The Pd-based electrode catalysts examined were Pd, Pd 3 Te 2 , and Pd 3 Te 1 (Figure 4a). The cycling experiments were conducted between −1.0 and +0.2 V. Analysis of the CVs presents the characteristic features of Pd in a solution of KOH. The results agree well with the results reported in the literature [28]. During the process of cathodic scan, the CVs exhibit intense and broad peaks in the regions spanning from −0.2 to −0.18 V and −0.83 V to −0.99 V. The origin of the peaks was attributed to the reduction and oxidation of Pd-O and hydrogen evolution (adsorption-desorption), respectively. Voltammograms recorded under conditions of positive (anodic) scan depict the onset of oxidation. The potential varied from −0.1 V to +0.05 V, and this was ascribed to the process of Pd-O/Pd-OH monolayer formation [33]. Pd 3 Te 1 generates a high current response, whereas Pd 3 Te 2 exhibits a broad current response under conditions of Pd-O reduction, and oxidation peaks are generated. These can be potentially attributed to the strong electronic interactions between Te and Pd. The increased Pd-O current and the negative shift observed for the unmodified Pd counterpart reveal that significantly high degrees of interactions are present between the Te electrons and Pd. The addition of Te in the crystal structure of Pd can also explain the observations. Moreover, the impact of introducing Te on the surface of Pd increases the amount of the oxide layer, and this can be attributed to the fact that Te is an oxygen species [3].
Two different Pd:Te ratios were investigated towards EOR as illustrated in Figure 4b. It was observed that the introduction of Te on Pd electrode catalysts reduced the onset potential and increased the oxidation current during EOR as observed in the forward peak current (i f ). This indicated a reduction in the activation energy and the rapid shuttling of electrons. The processes were induced by the modifications on the Pd surface with Te adatoms. The results can also be explained by the Te oxyphilic features. These features result in the generation of M(OH) x species at lower potential regions [34]. Te, as an oxophilic metal, may help to further oxidize the CO (EOR intermediates) species. Oxidation of this species can be promoted by the synergistic effect between Te and Pd. The generation of the synergistic effect facilitates the bifunctional mechanism and electronic effect, as described in previously reported papers [6]. The sharp reverse peak (backward peak current (i b )) is believed to reflect the high oxygen binding energy on the surface of Pd. The center of the d-band of Pd shifts under these conditions, and eventually, the species is oxidized further, resulting in the accumulation of EOR byproducts. The results are in strong agreement with DFT calculation results. A summary of the electrochemical activities for EOR is presented in Table 1.  The stability and the electron mobility of the best-performing catalyst (Pd 3 Te 2 in 0.5 M KOH) were investigated, as displayed in (Figure 5a). It is observed that the palladium reduction (PdO r ) peak increases with an increase in the scan rate under conditions of the cathodic scan. The results reflect the stability and the occurrence of the diffusion-controlled reactions on the electrocatalyst. The durability of Pd 3 Te 2 was further tested for EOR (0.5 M KOH and 0.1 M EtOH solution), and the results were compared with the results obtained for the Pd counterparts (Figure 5b). A rapid decay in current at the beginning (a few seconds) was observed for the electrocatalysts. This can be attributed to the accumulation of intermediate species (produced during EOR) on the active sites of the surface of the electrocatalyst. The polarization current decay reached a steady state, and at this point, a 2-fold increase in current was observed for Pd 2 Te 2 . The magnitude varied in the order Pd 2 Te 2 > Pd > Pd 3 Te 1 . This suggests that Pd 3 Te 2 is highly tolerant toward the CO ads intermediate. The high tolerance can be potentially attributed to the incorporation of the Te (oxygenated) species. During the last few minutes, Pd 3 T 1 overlaps with the Pd catalyst, resulting in current decay. The observation may be better described by the geometric effect or the low structural stability and inhomogeneous (agglomeration) surface morphology of the samples (AFM; Figure 3f). The results agree well with the results obtained using the DFT technique. The results explain the effect of the rugged surface and the oxygen binding energy recorded on the Pd surface. The results are in strong agreement with the results obtained by analyzing the CVs (Figure 4a,b).

Electrochemical Impedance Spectroscopy (Pd vs. TePd)
The best-performing catalyst (Pd 2 Te 2 ) was placed under scrutiny to examine the interfacial electron transfer process and compare the results with the results obtained by studying the unmodified Pd sample. Figure 6a shows the EIS plots (illustrated using a Nyquist plot). It can be observed that the electrode catalysts were characterized by a distinct size of the impedance arc (DIA). It has been previously reported that the size of DIA predicts the charge transfer resistance [35,36]. The Randles-Sevcik equivalent circuit (Figure 6a: insert) has been used to simulate the EIS elements (R s (R ct CPE)) from the semicircles to obtain the values, and the associated errors obtained were less than 4% [37]. As displayed in Figure 6, the semicircle corresponding to the Pd electrode catalyst is bigger than that recorded for the PdTe counterparts. The size of the semicircles recorded for PdTe indicates a low charge transfer resistance and fast electron transfer process [38]. These results are affected by the modification of the Pd surface (by Te oxygen species). The liberation of Pd from CO (EOR byproduct) promotes the EtOH adsorption/dehydrogenation process, as explained in previously reported papers [11,25] and Equation (S11). A bifunctional mechanism is thought to be associated with the process. This implies that Te oxygenated (Te promoter) species promote the process and help improve the electric conductivity of the bimetallic PdTe electrode catalyst. A rapid electron kinetic transfer is observed for this system, and the process is faster than the process observed for Pd [11]. These findings corroborate the results obtained using the CV and DFT techniques. The results revealed that the Pd-O bond strength could be improved, indicating the high electrocatalytic activity of the system during EOR. Figure 6b illustrates the Bode plot recorded with Pd-and PdTe-based electrode catalysts. The plots were recorded in EtOH in an alkaline media. The PdTe-based electrocatalysts were characterized by a high maximum phase angle (-θ max ) and metallic character (72 • and 84.7 • , close to 90 • ). The capacitance recorded was comparable to the total capacitances. The impact of Te on the surface of Pd was studied, and it was observed that Te modification helped increase the phase angle from −72 • to 84.7 • . This indicated an increased conductivity, a rapid electron transfer process, and a high electrocatalytic activity [39]. In addition, the electron shuttling process is promoted by the incorporation of Te, as revealed by the CVs presented in Figure 5b. Te insertion can potentially increase the concentration of the hydroxyl ( − OH) anions on the Pd surface. This result may also indicate the strong oxygen binding property (unlike the case of Pd) of PdTe. Similar results were obtained using the DFT technique (Figures 6 and 7). This may promote further attack and result in the stripping off of the intermediate poisonous carbonyl (CO) species [32]. This is also reflected by the bifunctional mechanism and the transformation of PdTe. The resistive electron movement observed at low-frequency values changes to capacitive behavior at high-frequency values. The results reflect the results obtained from the Nyquist plots ( Figure 6a) and CVs (Figure 4b). These plots reveal a decrease in the charge transfer resistance (Rct). A high peak current that confirms electron mobility is observed in this case. The results were compared with the results obtained by studying Pd.

Computational Details
The Vienna Ab initio Simulation Package (VASP) [26], based on the DFT technique [40,41], was used to simulate the ethanol oxidation activities of the PdTe alloy catalyst. The Kohn-Sham wavefunctions were solved to arrive at the results. The projector augmented wave (PAW) pseudopotential was used to replace the interaction between the core electrons [39,42]. The exchange-correlation energies of the electrons were described using Perdew, Burke, and Ernzerhof (PBE) [43]. The energies were revised using the Perdew-Burke-Ernzerhof functionals (RPBE) [44] under conditions of generalized gradient approximation (GGA) [45,46]. Pd was characterized by a traditional fcc structure. The system was designed, and the Te of the second (2nd) metal was replaced based on the ratio of Pd and Te. Five bulk structures were generated by substituting Te atoms at the edge and face sites (based on the Pd-to-Te ratio) to determine the configurations of the systems. Before designing the surface of Pd, Pd x Te y , each sample was fully relaxed in the bulk structure. The gamma point mesh was used (13 × 13 × 13; k-points), and the kinetic cutoff energy was set at 520 eV. The fully optimized bulk structures were cleaved into (111), (110), and (100) surfaces. The z-axis vacuum length was set at 20 Å to avoid interactions between the top and bottom surfaces. The bottom three layers were fixed to the surface, and the designed models were simulated using the gamma point (6 × 6 × 1) k-points. After obtaining the most stable surface structure, the oxygen atoms were placed on the surface of the possible configurations. Binding energies (∆E BE ) of oxygen were calculated on the Pd, Pd x Te y (111) surfaces. The DOS values (Figure 7) were used to execute the tetrahedron method with Blöchl corrections [42] to better illustrate the structural interactions.

The First-Principle Studies on Pd and PdTe Surfaces Used for Ethanol Oxidation
Based on the ab initio study, the electronic structure of Pd and Pd x Te y were calculated to reveal the relationship between the experimentally measured specific activity, the oxygen binding energies, and the positions of the d-band centers (Figure 8).   [47][48][49]. The theoretical d-band center of homogeneous Pd 2 Te 2 (111) and Pd 3 Te 1 (111) were found to be −2.14 and −1.97 eV, respectively (Figure 7a,b). These peaks shifted negatively with respect to the peak positions corresponding to Pd (111). Bader charge analyses [50] were conducted, and DOS corresponding to Pd 3 Te 2 (111) and Pd 3 Te 1 (111) were determined. The Bader charge density is calculated using Equation (2) as follows: In order to verify the experimental results, the five bulk structures of PdTe were designed from the conventional Pd fcc structure. The atoms were substituted at the face or edge sites. Three models having the lowest DFT internal energy were selected (Figure 9a-d). The theoretically calculated lattice parameters for x-, y-, and z-axes in the ground state of each Pd alloy were 3.940 Å for Pd bulk, 3.579, 4.422, and 5.989 Å for Pd 2 Te 2 , and 4.092 Å for Pd 3 Te 1 , respectively. The lattice parameters corresponding to Pd 2 Te 2 exhibit reduction along x, while elongation along the yand c-axes can be attributed to the process of replacement of the face site by the Te atoms. Figure 9a presents the theoretically calculated XRD patterns of Pd, Pd 2 Te 2 , and Pd 3 Te 1 . The results agree well with the experimentally obtained results for the synthesized PdTe alloys, expanding lattice parameters with increasing Te contents in the bulk lattice, as shown in Figure S3. As shown in Figure S3, the experimental XRD patterns is recorded for monometallic (Pd) and bimetallic (PdTe). The samples prepared by varying the Pd:Te ratio were analyzed. Four primary (Bragg's diffraction) peaks were observed in the profiles at 38.3 • (111), 51 • (200), 73 • (220), and 80 • (311) (this corresponds to the face-centered cubic (fcc) characteristic of metallic Pd). The results agree well with previously reported results [28,47]. Peaks corresponding to bimetallic nanofilms shifted to lower 2-theta angles ((111) and (220)) (with respect to the Pd counterpart). This change can be attributed to the insertion of Te into the Pd crystal lattice and the strong interactions between the Pd and Te adatoms. These features indicate the characteristic formation process of an alloy, and the results agree well with the results obtained using the DFT studies. The strong Pd/Au (111) peak suggested that the deposits preferred the (111) orientation. The low-intensity peaks corresponding to the PdTe catalyst appearing at the low (25 • , 20 • , 16 • ) 2θ regions can be attributed to the presence of elemental Te. The results reflect previously reported results [31]. These findings were consistent, as displayed in Figure S4 and as reported by this group in the literature [4]. These peaks are parallel to (111), indicating the existence of an epitaxial relationship between the catalyst and the substrate.
Here, σ is the surface energy; A is the surface area exposed to vacuum on both sides of the slab. E slab and E bulk represent the total energies corresponding to the surface and bulk models of Pd, Pd 2 Te 2 , and Pd 3 Te 1, and n is the number of the unit cells in the perfect bulk crystal that matches the total number of atoms in the slab model. For all Pd x Te y systems, the (111) slab models were thermodynamically stable.
The charge density plots corresponding to PdTe (111) and Pd 3 Te 1 (111) along the z-axes are presented in Figure 10c-f. The accumulation of the negative charge is represented by yellow, and the depletion of the negative charge is denoted in sky blue. The surface charge density corresponding to the two slab models was calculated by analyzing the charge density on the surface. Negative charge accumulated on Pd (−1.14 e), which is the catalytic active site as compared to less amount of negative charge accumulated at the depleted Te site (+0.12 e) on the PdTe surface. EOR in alkaline media has been explained by various research groups [11,[51][52][53]. The process can be explained by the following half-cell reaction C-C bond breaking may not be observed for EOR byproducts present at the Pd interface. This can potentially result in the generation of acetyl (CH 3 CO ads a) products and not CO 2 . The mechanism that can best explain the adsorption and generation of species has been presented. The extent of CH 3 CO ads and OH ads coverage realized [54,55] influences the mechanism presented in the Supplementary Materials and literature reports [7,[54][55][56]. The rate-determining step was identified by analyzing Equations (S8)-(S11) (Supplementary Materials) [7,54]. The OH ads units adsorbed on the electrode play an important role during EOR. The adsorption of the byproducts can be prevented in the presence of the adsorbed hydroxyl (-OH) species, and free active sites are regenerated on the catalyst under these conditions. Furthermore, it is important to study the oxygen adsorption energy to understand the catalytic activity for ethanol. The DOS of electrode catalysts at the Fermi level plays a crucial role in the chemisorption of the adsorbates. It significantly influences the performance of the catalyst. Therefore, the binding energy of oxygen was calculated based on the surface (111) structure of Pd, Pd 2 Te 2 , and Pd 3 Te 1 . The oxygen adsorption sites were investigated for possible locations depending on the top, hcp, and fcc structures, as shown in Figure 10. The three structures are presented in Figure 10b-d, and the energies were obtained using Equation (5) as follows: where E PdxTey(111)-O* denotes the oxygen adsorbed on the catalyst surface, E PdxTey is the internal energy of the (111) slab models of Pd and Pd x Te y , and E O is the DFT-based energy of the isolated single oxygen atom. Analysis of the oxygen binding energies revealed that the higher the d-band center positions, the greater the bond energy between Pd and the oxygen species (e.g., OH − ). Comparisons were made with respect to the pure Pd surface. The strong bond formed by oxygen can potentially result in the induction of a catalytic oxidation reaction. We observed that the d-band shifted downward, indicating the formation of weak bonds between oxygen species and Pd surface. Low activity is observed for the Pd 3 Te 1 catalyst. This result shows the discrepancy with the experimentally obtained results. The overall activity of the catalyst is a cumulative effect of the geometric and electronic effects of the catalyst [57][58][59]. When the Te content with bcc structure increased, the lattice of the Pd 2 Te 2 system was significantly distorted, and elongation along the y-and c-axes was observed ( Figure 10c). Elongation results in an increase in the free surface energy of the rugged Pd 3 Te 2 (111) system (Figure 10c).
The geometric effect results in the generation of high binding energy for oxygen ( Figure 8). This phenomenon represents that OH − ions can be readily adsorbed onto Pd 3 Te 2 , and this was not the case for Pd. The OH ads units adsorbed on the electrode will oxidize the ethanol byproducts such as CH 3 CO ads and CH 3 COOH by interacting with the surface OH − ions [55]. For the case of the Pd 3 Te 1 catalysts, it was observed that the farther the d−band center from the Fermi level, the weaker the interaction of the system with the adsorbents (e.g., oxygen). This indicated that the catalytic activity of the system was lower than the catalytic activity of Pd. High binding energy was recorded for Pd 3 Te 2 , and this could be attributed to the geometric effect. However, a low-lying d-band position was recorded for this system. An increase in the Te content indicates a geometric effect. Under these conditions, the geometric effect is more dominant than the effect attributable to electronic structural changes. This reveals the presence of rugged surfaces that can potentially prevent the formation of EOR byproducts and reveal the active sites of the catalysts. These results are consistent with the experimental results that reveal that the catalytic activity increases during ethanol oxidation reaction ( Figure S3).

Conclusions
Electrodeposition of monometallic Pd and bimetallic PdTe (by varying the Pd-to-Te ratio) was successfully realized following the E-ALD method. The E-ALD cycles proceeded through multiple steps of SLRRR. Cu-UPD was used for the deposition of Te in the presence of oxidative Te-UPD. The introduction of the Te adlayer on Pd improved the electro-oxidation property for EOR. This marks the development of the first PGMchalcogenide (PdTe)-based electrocatalyst following the process of E-ALD. The bimetallic PdTe system exhibited improved catalytic activity and electron transfer kinetics (compared to Pd). The results were arrived at using the CV and EIS techniques. Analysis of the XRD patterns confirmed the formation of the PdTe alloy. The negative 2θ shift indicated the interaction between Te and the Pd crystal lattice. Te plays a crucial role in fine-tuning the structural morphology and the electrochemical performance of the Pd-based electrode catalyst. The insertion of Te was confirmed by studying the morphology distortion of the system. Dendritic structures and surface grain boundaries were formed, and these were obtained using the SEM and AFM techniques, respectively. The EOR mechanisms and electrocatalytic activity of Pd (111) and PdTe (111)-based catalysts were successfully investigated using the ab initio and DFT techniques. Strong oxygen bond energy was recorded following the incorporation of Te into Pd systems. This indicates the presence of weak CO-like species and improved activity. The DFT results corroborate the theoretically calculated results. A combination of geometric and electronic effects could effectively describe the effect of the Te content on the performance of the PdTe alloy. The electronic structural effect was found to be the dominant effect. The E-ALD method can be controlled to deposit PGM-based systems doped with chalcogenides for the fabrication of ethanol fuel cells.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/nano12203607/s1, Figure S1: Defined slab models on oxygen adsorbed Pd x Te y (111), (110), (100) designed from bulk structure to find out which structure is the most stable surface. The three layers below were fixed to simulate surface, and identical 24 atoms were used to describe the same condition; Figure S2: Defined slab models on oxygen adsorbed Pd x Te y (111) designed from fcc bulk structure for calculating oxygen adsorption energies. The three layers below were fixed to simulate surface, and identical 24 atoms were used to describe the same condition; Figure S3: Experimentally XRD pattern of Au, Pd, Pd 3 Te 1 and Pd 2 Te 2 ; Figure S4: Experimental XRD pattern of Pd and PdTe nano film, and insert: zoomed Pd/Au (111). Data Availability Statement: The data presented in this manuscript will be available upon request from the corresponding author.