Single Atom Alloys Segregation in the Presence of Ligands.

Single atom alloys (SAAs) have gained remarkable attention due to their tunable properties leading to enhanced catalytic performance, such as high activity and selectivity. The stability of SAAs is dictated by surface segregation, which can be affected by the presence of surface adsorbates. Research efforts have primarily focused on the effect of commonly found catalytic reaction intermediates, such as CO and H, on the stability of SAAs. However, there is a knowledge gap in understanding the effect of ligands from colloidal nanoparticle (NP) synthesis on surface segregation. Herein, we combine density functional theory (DFT) and machine learning to investigate the effect of thiol and amine ligands on the stability of colloidal SAAs. DFT calculations revealed rich segregation energy (Eseg) data of SAAs with d8 (Pt, Pd, Ni) and d9 (Ag, Au, Cu) metals exposing (111) and (100) surfaces, in the presence and absence of ligands. Using these data, we developed an accurate four-feature neural network using a multilayer perceptron regression (NN MLP) model. The model captures the underlying physics behind surface segregation in the presence of adsorbed ligands by incorporating features representing the thermodynamic stability of metals through the bulk cohesive energy, structural effects using the coordination number of the dopant and the ligands, the binding strength of the adsorbate to the metals, strain effects using the Wigner-Seitz radius, and electronic effects through electron affinities. We found that the presence of ligands makes the thermodynamics of segregation milder compared to the bare (nonligated) SAA surfaces. Importantly, the adsorption configuration (e.g., top vs bridge) and the binding strength of the ligand to the SAA surface (e.g., amines vs thiols) play an important role in altering the Eseg trends compared to the bare surface. We also developed an accurate NN MLP model that predicts Eseg in the presence of ligands to find thermodynamically stable SAAs, leading to the rapid and efficient screening of colloidal SAAs. Our model captures several experimental observations and elucidates complex physics governing segregation at nanoscale interfaces.

Demystifying the Chemical Ordering of Multimetallic Nanoparticles.

ConspectusMultimetallic nanoparticles (NPs) have highly tunable properties due to the synergy between the different metals and the wide variety of NP structural parameters such as size, shape, composition, and chemical ordering. The major problem with studying multimetallic NPs is that as the number of different metals increases, the number of possible chemical orderings (placements of different metals) for a NP of fixed size explodes. Thus, it becomes infeasible to explore NP energetic differences with highly accurate computational methods, such as density functional theory (DFT), which has a high computational cost and is typically applied to up to a couple of hundred metal atoms. Here, we demonstrate a methodology advancing NP simulations by effectively exploring the vast materials space of multimetallic NPs and accurately identifying the ones with the most thermodynamically preferred chemical orderings. With accuracies reaching that of DFT, our methodology is applicable to practically any NP size, shape, and metal composition. We achieve this by significantly advancing the bond-centric (BC) model, a physics-based model that has been previously shown to rapidly predict bimetallic NP cohesive energies (CEs). Specifically, the BC model is trained in a way to understand how the bimetallic bond strength changes under different coordination environments present on a NP and how the metal composition of every site affects the detailed coordination environment using fractional coordination numbers. This newly modified BC model leads to an improvement from 0.331 (original model) to 0.089 eV/atom in CE predictions when compared to DFT values on a robust data set of 90 different NPs consisting of PtPd, AuPt, and AuPd NPs with varying compositions and chemical orderings. By incorporating the modified BC model into an in-house-developed genetic algorithm (GA) we can effectively and accurately predict the most stable chemical orderings of large, realistic bimetallic NPs consisting of thousands of metal atoms. This is demonstrated on AuPd bimetallic NPs, a challenging system due to the similarity in the cohesion of the two metals. By training our BC model using a unique DFT calculation on a bimetallic NP (one calculation for two metals combining together), we expand to explore the chemical ordering of multimetallic NPs. We first demonstrate the application of our methodology on a AuPdPt NP and validate our stability predictions with literature data. Then, we effectively explore the vast materials space of multimetallic NPs consisting of combinations of Au, Pt, and Pd as a function of metal composition. Our thermodynamic stability trends are presented in a ternary diagram revealing detailed, and yet, unexpected chemical ordering trends. Our computational framework can aid both experimental and computational researchers toward effectively screening multimetallic NP stability. Moreover, we provide an outlook of how this framework can be applied to catalyst discovery, high-entropy alloys, and single-atom alloys.

Size-dependent shape distributions of platinum nanoparticles.

While it is well established that nanoparticle shape can depend on equilibrium thermodynamics or growth kinetics, recent computational work has suggested the importance of thermal energy in controlling the distribution of shapes in populations of nanoparticles. Here, we used transmission electron microscopy to characterize the shapes of bare platinum nanoparticles and observed a strong dependence of shape distribution on particle size. Specifically, the smallest nanoparticles (<2.5 nm) had a truncated octahedral shape, bound by 〈111〉 and 〈100〉 facets, as predicted by lowest-energy thermodynamics. However, as particle size increased, the higher-energy 〈110〉 facets became increasingly common, leading to a large population of non-equilibrium truncated cuboctahedra. The observed trends were explained by combining atomistic simulations (both molecular dynamics and an empirical square-root bond-cutting model) with Boltzmann statistics. Overall, this study demonstrates experimentally how thermal energy leads to shape variation in populations of metal nanoparticles, and reveals the dependence of shape distributions on particle size. The prevalence of non-equilibrium facets has implications for metal nanoparticles applications from catalysis to solar energy.