skip to main content


Title: Identification of optimally stable nanocluster geometries via mathematical optimization and density-functional theory
Small nanoparticles, a.k.a. nanoclusters, of transition metals have been studied extensively for a wide range of applications due to their highly tunable properties dependent on size, structure, and composition. For these small particles, there has been considerable effort towards theoretically predicting what is the most energetically favorable arrangement of atoms when forming a nanocluster. In this work, we develop a computational framework that couples density-functional theory calculations with mathematical optimization modeling to identify highly stable, mono-metallic transition metal nanoclusters of various sizes. This is accomplished by devising and solving a rigorous mathematical optimization model that maximizes a general cohesive energy function to obtain nanocluster structures of provably maximal cohesiveness. We then utilize density-functional theory calculations and error term regression to identify model corrections that are necessary to account with better accuracy for different transition metals. This allows us to encode metal-specific, analytical functions for cohesive energy into a mathematical optimization-based framework that can accurately predict which nanocluster geometries will be most cohesive according to density-functional theory calculations. We employ our framework in the context of Ag, Au, Cu, Pd and Pt, and we present sequences of highly cohesive nanoclusters for sizes up to 100 atoms, yielding insights into structures that might be experimentally accessible and/or structures that could be used as model nanoclusters for further study.  more » « less
Award ID(s):
1634880
NSF-PAR ID:
10201311
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Molecular Systems Design & Engineering
Volume:
5
Issue:
1
ISSN:
2058-9689
Page Range / eLocation ID:
232 to 244
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Determining the energetically most favorable structure of nanoparticles is a fundamentally important task towards understanding their stability. In the case of bimetallic nanoclusters, their vast configurational space makes it especially challenging to find the global energy optimum via experimental or computational screening. To that end, this work proposes a two-step optimization-based design framework to address this hard combinatorial problem. Given a nanocluster of fixed shape, a rigorous mixed-integer linear programming model is formulated based on a bond-centric cohesive energy function to identify the most cohesive bimetallic configuration for a given composition. This capability is coupled with a metaheuristic strategy that searches over the space of nanocluster shapes to obtain optimal structures. We apply our proposed methodology on AgCu, AuAg and CuAu systems, quantifying how the size and composition of a nanocluster influences its overall cohesion. Furthermore, we observe various synergistic effects between Cu and Au in promoting cohesive energy, while multiple segregation patterns are identified in all three studied binary systems. Our methodology serves as an efficient computational tool for investigating bimetallic nanoclusters stability properties as well as provides model nanoclusters for further investigations. 
    more » « less
  2. Atomically precise nanoclusters of coinage metals in the 1–3 nm size regime have been intensively pursued in recent years. Such nanoclusters are attractive as they fill the gap between small molecules (<1 nm) and regular nanoparticles (>3 nm). This intermediate identity endows nanoclusters with unique physicochemical properties and provides nanochemists opportunities to understand the fundamental science of nanomaterials. Metal nanoparticles are well known to exhibit plasmon resonances upon interaction with light; however, when the particle size is downscaled to the nanocluster regime, the plasmons fade out and step-like absorption spectra characteristic of cluster sizes are manifested due to strong quantum confinement effects. Recent research has revealed that nanoclusters are commonly composed of a distinctive kernel and a surface-protecting shell (or staple-like metal–ligand motifs). Understanding the kernel configuration and evolution is one of the central topics in nanoscience research. This Review summarizes the recent progress in identifying the growth patterns of atomically precise coinage nanoclusters. Several basic kernel units have been observed, such as the M 4 , M 13 and M 14 polyhedrons (where, M = metal atom). Among them, the tetrahedral M 4 and icosahedral M 13 units are the most common ones, which are adopted as building blocks to construct larger kernel structures via various fusion or aggregation modes, including the vertex- and face-sharing mode, the double-strand and alternate single-strand growth, and cyclic fusion of units, as well as the fcc-based cubic growth pattern. The identification of the kernel growth pathways has led to deeper understanding of the evolution of electronic structure and optic properties. 
    more » « less
  3. DNA-stabilized silver nanoclusters (AgN-DNAs) are a class of nanomaterials comprised of 10-30 silver atoms held together by short synthetic DNA template strands. AgN-DNAs are promising biosensors and fluorophores due to their small sizes, natural compatibility with DNA, and bright fluorescence---the property of absorbing light and re-emitting light of a different color. The sequence of the DNA template acts as a "genome" for AgN-DNAs, tuning the size of the encapsulated silver nanocluster, and thus its fluorescence color. However, current understanding of the AgN-DNA genome is still limited. Only a minority of DNA sequences produce highly fluorescent AgN-DNAs, and the bulky DNA strands and complex DNA-silver interactions make it challenging to use first principles chemical calculations to understand and design AgN-DNAs. Thus, a major challenge for researchers studying these nanomaterials is to develop methods to employ observational data about studied AgN-DNAs to design new nanoclusters for targeted applications. In this work, we present an approach to design AgN-DNAs by employing variational autoencoders (VAEs) as generative models. Specifically, we employ an LSTM-based β-VAE architecture and regularize its latent space to correlate with AgN-DNA properties such as color and brightness. The regularization is adaptive to skewed sample distributions of available observational data along our design axes of properties. We employ our model for design of AgN-DNAs in the near-infrared (NIR) band, where relatively few AgN-DNAs have been observed to date. Wet lab experiments validate that when employed for designing new AgN-DNAs, our model significantly shifts the distribution of AgN-DNA colors towards the NIR while simultaneously achieving bright fluorescence. This work shows that VAE-based generative models are well-suited for the design of AgN-DNAs with multiple targeted properties, with significant potential to advance the promising applications of these nanomaterials for bioimaging, biosensing, and other critical technologies. 
    more » « less
  4. Electrochemical conversion of carbon dioxide (CO 2 ) to chemicals or fuels can effectively promote carbon capture and utilization, and reduce greenhouse gas emission but a serious impediment to the process is to find highly active electrocatalysts that can selectively produce desired products. Herein, we have established the design principles based on the density functional theory calculations to screen the most promising catalysts from the family of coordinately unsaturated/saturated transition metal (TM) embedded into covalent organic frameworks (TM-COFs). An intrinsic descriptor has been discovered to correlate the molecular structures of the active centers with both the activity and selectivity of the catalysts. Among all the catalysts, the coordinately unsaturated Ni-doped covalent triazine framework (Ni-CTF) is identified as one of the best electrocatalysts with the lowest overpotential (0.34 V) for CO 2 reduction toward CO while inhibiting the formation of the side products, H 2 and formic acid. Compared with coordinately saturated TM-COFs and noble metals ( e.g. Au and Ag), TM-CTFs exhibit higher catalytic activity and stronger inhibition of side products. The predictions are supported by previous experimental results. This study provides an effective strategy and predictive tool for developing desired catalysts with high activity and selectivity. 
    more » « less
  5. null (Ed.)
    Thiolate-protected metal nanoclusters (TPNCs) have attracted great interest in the last few decades due to their high stability, atomically precise structure, and compelling physicochemical properties. Among their various applications, TPNCs exhibit excellent catalytic activity for numerous reactions; however, recent work revealed that these systems must undergo partial ligand removal in order to generate active sites. Despite the importance of ligand removal in both catalysis and stability of TPNCs, the role of ligands and metal type in the process is not well understood. Herein, we utilize Density Functional Theory to understand the energetic interplay between metal–sulfur and sulfur–ligand bond dissociation in metal–thiolate systems. We first probe 66 metal–thiolate molecular complexes across combinations of M = Ag, Au, and Cu with twenty-two different ligands (R). Our results reveal that the energetics to break the metal–sulfur and sulfur–ligand bonds are strongly correlated and can be connected across all complexes through metal atomic ionization potentials. We then extend our work to the experimentally relevant [M 25 (SR) 18 ] − TPNC, revealing the same correlations at the nanocluster level. Importantly, we unify our work by introducing a simple methodology to predict TPNC ligand removal energetics solely from calculations performed on metal–ligand molecular complexes. Finally, a computational mechanistic study was performed to investigate the hydrogenation pathways for SCH 3 -based complexes. The energy barriers for these systems revealed, in addition to thermodynamics, that kinetics favor the break of S–R over the M–S bond in the case of the Au complex. Our computational results rationalize several experimental observations pertinent to ligand effects on TPNCs. Overall, our introduced model provides an accelerated path to predict TPNC ligand removal energies, thus aiding towards targeted design of TPNC catalysts. 
    more » « less