##### Authors: Micha Polak* and Leonid Rubinovich
Affiliation(s):
Department of Chemistry, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel
Dates:
Received: 27 December, 2016; Accepted: 27 January, 2017; Published: 31 January, 2017
*Corresponding author:
Micha Polak, Department of Chemistry, Ben-Gurion University of the Negev, Beer-Sheva, 84105, Israel, E-mail: @
Citation:
Polak M, Rubinovich L (2017) Prediction of Enhanced Dimerization inside Dilute Alloy Nanoparticles. Int J Nanomater Nanotechnol Nanomed 3(1): 023-026. DOI: 10.17352/2455-3492.000017
Â© 2017 Polak M, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Keywords:
Nanoconfinement; Alloy nanoparticles; Dimerization; Chemical-equilibrium; Mixing-entropy

According to a unique nano-confinement effect of entropic origin, predicted by us several years ago for the equilibrium state of chemical reactions, the equilibrium-constant and extent are greatly enhanced depending on the small number of molecules involved, and for many reactions also on the nano-space size. This work explored the validity of this effect in the case of elemental dimerization reactions within dilute alloy nanoparticles with separation tendency, PdIr cuboctahedra in particular. Employing a simple model for the system energetics, computations based on the exact canonical partition-function reveal nano-confinement induced Ir2 dimer stabilization within Pd surface segregated nanoparticles, reflected e.g., by up to ~ 60% increased dimerization extent and by doubling of the 1n KD vs. 1/T slope, as compared to the macroscopic thermodynamic limit. The dual role of the configurational entropy, namely mixing of Ir/Ir2 vs. Pd/Ir is elucidated. Study based on more elaborate energetic models is desirable as the next step of this research.

### Introduction

As shown by us some years ago [1], the involvement of small numbers of molecules in a closed reaction space can considerably affect the chemical equilibrium. Thus, significant increase of the reaction extent was predicted for the case of exothermic reactions as compared to that of the corresponding macroscopic system (in the thermodynamic limit, TL). This nanoconfinement entropic effect on chemical equilibrium (NCECE) is universal in a sense that it stems only from the limited amount of molecules in a reaction mixture, resulting in a reduced number of mixed reactant-product microstates in the closed system [1,2]. The effect was predicted for the case of nucleotide dimerization within molecular cages [2] and verified for DNA hybridization [3] that was studied experimentally before [4]. The NCECE is relevant to several advanced routes for the synthesis of encapsulated organic molecules, metallic or inorganic nanoclusters, and other nanoscale structures and applications.

While the NCECE is expected mostly in molecular nanosystems, this effect that arises from fundamental mixing-entropy variations can in principle be manifested also in other reaction classes, e.g., in alloy systems. The present study focuses on the distinct case of constituent atom equilibrated dimerization within phase-separating alloy nanoparticles (NPs). Namely, the main question of this study concerns whether confined dimerization can exhibit significant NCECE effects. Being a first stage in phase separation (or the main process in very dilute alloys) that can affect chemical and physical properties, dimerization inside NPs seems to be a sensible issue for theoretical-computational modeling. For this goal fcc-based PdIr NPs have been chosen, because of the strong Pd tendency to segregate at the alloy surface, making the Pd rich core a natural confined space for Ir solute atoms entirely surrounded by 12 Pd nearest-neighbor (NN) Pd atoms. Moreover, PdIr nanoparticles are efficient catalysts, e.g., promoting preferential oxidation of CO in the presence of excess H2 in proton exchange membrane fuel cells [5-7]. Recently, chemical-order in “magic-number” PdIr NPs was studied by means of density functional theory (DFT) computations [8] and compared to those obtained by the Free-Energy Concentration Expansion Method (FCEM) [9], using derived coordination-dependent bond-energy variations (CBEV) [10] and by the Birmingham Cluster Genetic Algorithm with the Gupta potential [11].

In this work, Ir dimerization in cuboctahedron dilute cores of the magic-number size series [12] is studied (namely, 13, 55, 147, 309, 561 and 923 atom cores inside 55, 147, 309, 561, 923 and 1415 atom alloy NPs, respectively) in the framework of the NN interaction approximation. Energetic equivalence of all intra-core Ir atomic sites and of all Ir2 dimer locations is assumed, since, as noted above, they are all surrounded by Pd only (Figure 1). This includes also the subsurface sites, namely possible CBEV effects are ignored in this study. Correspondingly, the energy of the intra-core dimerization “reaction” equals the doubled bulk Effective Pair Interaction (EPI), ∆E = 2V = (wPdPd + wPdIr – 2 wIrIr) ≈ -74 meV, as derived from DFT-computed low-temperature formation enthalpy [13]. Formulas for the equilibrium dimerization extent, equilibrium constant and the NP entropy are derived below based on canonical statistical mechanics and applied in MATLAB computations. Then, the results are evaluated for possible NCECE manifestations by comparison to the corresponding values of a hypothetical macroscopic alloy with the same composition and dimerization energy.

In the employed lattice-gas model [14] microstates correspond to different arrangements of two Ir solute atoms on n available atomic sites, and to the locations of a single Ir2 dimer at available n' bonds in the NP core (see examples in Figure 1b). Concentrations of Ir atoms and dimers are related to the nanoconfined (NC) reaction extent, ${\xi }_{D}^{NC}$ (${\xi }_{D}^{NC}=0$ ) for pure reactants, for pure products),

1. ##### Figure 1:
(a) Illustration of 2Ir ⇔ Ir2 equilibrated dimerization inside the 309-atom core of Pd surface-segregated 561-atom cuboctahedral nanoparticle. (b) 2D schematics of two representative configurations (“microstates”), corresponding to distributions of Ir atoms on atomic sites (left) and of Ir2 dimers as nearest-neighbor bonds.

Since the lattice-gas chemical-potential [14] depends on a “reduced concentration”,$\frac{c}{1-c}$ , the dimerization equilibrium constant is given by,

In order to calculate the above, the number of ways to arrange just two Ir atoms is given by the binomial coefficient, ( 2 n )= n(n1) 2 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbnvMCYL2DLfgDOvMCaeXatLxBI9gBaerbd9wDYLwzYbItLDharuavP1wzZbItLDhis9wBH5garqqtubsr4rNCHbGeaGakY=3j0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr0=vqpWqaaeaabiGaciaacaqabeaadaabauaaaOqaamaabmaabaWaa0baaSqaaiaaikdaaeaacaWGUbaaaaGccaGLOaGaayzkaaGaeyypa0ZaaSaaaeaacaWGUbGaaiikaiaad6gacqGHsislcaaIXaGaaiykaaqaamaaCaaaleqabaGaaGOmaaaaaaaaaa@4678@ . Since arrangements corresponds to dimer locations, it has to be subtracted in the canonical partition function,

Thus, in addition to the energetics the dimerization is controlled by the overall numbers of intra-core atomic sites and NN bonds (straightforwardly computed by MATLAB). Starting from Q, the equilibrium properties of the system can be obtained, such as the reaction extent, namely, the probability of dimer formation,

Then, the corresponding concentrations are obtained by eq.1 in order to compute the equilibrium constant (via eq.2).

For obtaining the reference TL equilibrium constant, the reaction energy and entropy are needed,

Thus, in the case of very large NPs and very high temperatures the reaction extent obtained by eqs.3-4, ${\xi }_{D}^{\infty }\approx 2n\text{'}/{n}^{2}$ , is substituted in eq.1 giving the concentrations, ${c}_{a}^{\infty }=2/n$ , ${c}_{d}^{\infty }=2/{n}^{2}$ . Then, according to eq.2 the high temperature limit,${K}_{D}^{\infty }=1/2$ , and eq.5 gives $\Delta S=-k\mathrm{ln}2$ . The ${K}_{D}^{TL}$ of the corresponding hypothetical macroscopic alloy computed for any T (by eq.5) together with eqs.1-2 yields,

This equation has to be solved numerically in order to find the TL reaction extent, $\delta \xi =\left({S}^{NC}/{S}^{TL}-1\right)100%$ .

To conclude the statistical-mechanical derivation, the NP configurational entropy per site is given according to the Shannon formula [15] by,

According to the computations, the dimerization extent under nanoconfinement is significantly enhanced compared to the TL value based on the same ${c}_{a}^{{}^{\mathrm{max}}}$ and ${c}_{d}^{{}^{\mathrm{max}}}$ (Figure 2a), and can be expressed by its relative increase, $\delta \xi =\left({\xi }_{D}^{NC}/{\xi }_{D}^{TL}-1\right)100%$ that reaches about 60% for core size n=13 (Figure 2b). This phenomenon is due to the NCECE that gradually weakens with increased NP size (“nanospace size effect” [2]), until for the 923-atom NP core the two extent curves almost coincide (Figure 2a). As shown in Figure 2c, the equilibrium constant is enhanced too and most of the log KD vs. 1/T slopes are twice larger than the TL slope. Consequently, the logarithm ratios converge to 2 at the low temperature limit (Figure 2d). This slope doubling is characteristic e.g., of exothermic addition and dimerization reactions (A + B = C, 2A=A2) in nanoconfined systems, and relates directly to the sum of the reactant stoichiometric coefficients [1]. Regarding the significant role of configurational entropy, it originates from two mixing contributions: (i) iridium (atomic, dimerized) with palladium atoms, and (ii) iridium atoms with dimers. The dissociation backward endothermic reaction. Ir2→2Ir, is promoted by temperature increase and by the NP-size dependent alloy dilution, (Figure 2a), in accordance with the Le Chatelier’s principle. Thus, as can be further seen in Figure 3a, the entropy per site increases with temperature and the region of its most intensive rise (maximum slope) shifts to lower temperatures. Furthermore, the entire entropy plot shifts downwards with increasing NP size due to the dilution. All the above phenomena are reflected also in the TL curves towards which the NP curves approach with growing size. The NCECE stems primarily from deficiency of type (ii) mixing in NPs reflected in lowering of the entropy under nano-confinement, SNC, compared to the STL plot (Figure 3a). Since it is completely absent in the 2-atoms/1-dimer case under consideration, the corresponding relative decrease in entropy, $\delta \xi =\left({S}^{NC}/{S}^{TL}-1\right)100%$ (Figure 3b), is associated mainly with this lack of type (ii) mixing contribution in the NP. The gradual shrinking of δS with increasing NP size is attributed to the growing relative contribution of the nanospace-related mixing entropy, type (i). The δS minima, associated with the above-mentioned most intensive dissociation regions, follow their shift to lower temperatures with NP size increase. Furthermore, the minimum values obey distinct, nearly linear dependence on their temperatures (Figure 3c).

1. ##### Figure 2:
The temperature dependence of the dimerization extent (a), its relative increase, δξ (b), the equilibrium “constants” (c), and the NC/TL ratios (d), computed for two Ir solute atoms in cores of cuboctahedral Pd- Ir NPs containing 13, 55, 147, 309 and 923 atoms (size magic-numbers). The macroscopic TL values are given for comparison (dashed lines in (a) and (c)).

1. ##### Figure 3:
The dependence on temperature computed for different size NPs: (a) the configurational entropy per site (the numbers of core sites are indicated); (b) the corresponding relative decrease in entropy δS; (c) the magnitude of the δS minimum.

To summarize, a theoretical-computational study is introduced concerning elemental dimerization within cuboctahedral PdIr nanoparticles as a model system. Considering two solute Ir atoms confined under segregated Pd surface layer as lattice gas, statistical-mechanical computations reveal for the first time that “nanoconfinement entropic effect on chemical equilibrium” (NCECE) induces significant enhancement of equilibrium solute dimerization, compared to the bulk alloy with the same composition. Being a kind of first stage of phase separation, the enhanced dimerization indicates strengthening of this tendency, unlike the known phenomena of diminished phase separation temperatures in nondilute alloy NPs. The predicted NCECE-induced dimerization enhancement can be relevant to alloy NP catalytic properties due to possible ligand effects such as bandwidth changes via hybridization of d-states of surface atoms with subsurface atoms [16]. Moreover, coordination-dependent bond-energy variations (CBEV) can make certain sites in the subsurface layer energetically preferable for solute atoms and dimers, thereby further affecting the alloy particle catalytic properties. Hence, a next step of the PdIr dimerization research, which is based on this energetic model, is desirable.

1. Polak M, Rubinovich L (2008) Nanochemical Equilibrium Involving a Small Number of Molecules: A Prediction of a Distinct Confinement Effect. Nano Lett 8: 3543-3547. Link: https://goo.gl/lCl70W
2. Polak M, Rubinovich L (2011) Remarkable nanoconfinement effects on chemical equilibrium manifested in nucleotide dimerization and H-D exchange reactions. Phys Chem Chem Phys 13: 16728-16734. Link: https://goo.gl/U6QSjU
3. Rubinovich L, Polak M (2013) The Intrinsic Role of Nanoconfinernent in Chemical Equilibrium: Evidence from DNA Hybridization. Nano Lett 13: 247-2251. Link: https://goo.gl/20hgdS
4. Shon MJ, Cohen AE (2012) Mass Action at the Single-Molecule Level. J Am Chem Soc 134: 14618-14623. Link: https://goo.gl/8XLWOY
5. Shen SY, Zhao TS, Xu JB (2010) Carbon-supported bimetallic Pdlr catalysts for ethanol oxidation in alkaline media. Electrochimica Acta 55: 9179-9184. Link: https://goo.gl/uprBoF
6. Morfin F, Nassreddine S, Rousset JL, Piccolo L (2012) Nanoalloying Effect in the Preferential Oxidation of CO over Ir-Pd Catalysts. Acs Catalysis 2: 2161-2168. Link: https://goo.gl/YGhUNU
7. Zlotea C, Morfin F, Nguyen TS, Nguyen NT, Nelayah J, et al. (2014) Nanoalloying bulk-immiscible iridium and palladium inhibits hydride formation and promotes catalytic performances. Nanoscale 6: 9955-9959. Link: https://goo.gl/k59M8e
8. Davis JB, Johnston RL, Rubinovich L, Polak M (2014) Comparative modelling of chemical ordering in palladium-iridium nanoalloys. Journal of Chemical Physics 141: 224307. Link: https://goo.gl/pLt412
9. Polak M, Rubinovich L (2008) The interplay of surface segregation and atomic order in alloys. Surface Science Reports 38: 127-194. Link: https://goo.gl/Nj3ZC5
10. Rubinovich L, Polak M (2009) Prediction of distinct surface segregation effects due to coordination-dependent bond-energy variations in alloy nanoclusters. Phys Rev B 80: 045404. Link: https://goo.gl/cIqrGq
11. Davis JBA (2014) private communication.
12. Baletto F, Ferrando R (2005) Structural properties of nanoclusters: Energetic, thermodynamic, and kinetic effects. Rev Mod Phys 77: 371-423. Link: https://goo.gl/50uyPV
13. Kolb B, Mueller S, Botts DB, Hart GLW (2006) Ordering tendencies in the binary alloys of Rh, Pd, Ir, and Pt: Density functional calculations. Phys Rev B 74: 144206. Link: https://goo.gl/oNgSDp
14. Hill TL (1986) An Introduction to Statistical Thermodynamics. Courier Dover Publications. Link: https://goo.gl/0Hl2Vg
15. Ford I (2013) Statistical physics: an entropic approach. A John Wiley & Sons, Ltd. Publication. Chichester, United Kingdom. Link: https://goo.gl/uZKAc8
16. Bligaard T, Norskov JK (2007) Ligand effects in heterogeneous catalysis and electrochemistry. Electrochimica Acta 52: 5512-5516. Link: https://goo.gl/3t7bqw