A Water Solvation Shell Can Transform Gold Metastable Nanoparticles in the Fluxional Regime

Minimum Free-Energy Shapes of Ag Nanocrystals: Vacuum vs Solution

We use two variants of replica-exchange molecular dynamics (MD) simulations, parallel tempering MD and partial replica exchange MD, to probe the minimum free-energy shapes of Ag nanocrystals containing 100-200 atoms in a vacuum, ethylene glycol (EG) solvent, and EG solvent with a PVP polymer containing 100 repeat units. Our simulations reveal a shape intermediate between a Dh and an Ih, a Dh-Ih, that has distinct structural signatures and magic sizes. We find several prominent features associated with entropy: pure FCC nanocrystals are less common than FCC crystals containing stacking faults, and crystals with the minimum potential energy are not always preferred over the range of relevant temperatures. The shapes of the nanocrystals in solution are influenced by the chemical identities of the solution-phase molecules. Comparing Ag nanocrystal shapes in EG to those in an EG+PVP solution, we find more icosahedra in EG and more decahedra in EG+PVP across all of the nanocrystal sizes probed in this study. At certain critical sizes, nanocrystal shapes can change dramatically with the addition and removal of a single atom or with a change in temperature at a fixed size. The information in our study could be useful in efforts to devise processing routes to achieve selective nanocrystal shapes.

Theory of Anisotropic Metal Nanostructures

A significant challenge in the development of functional materials is understanding the growth and transformations of anisotropic colloidal metal nanocrystals. Theory and simulations can aid in the development and understanding of anisotropic nanocrystal syntheses. The focus of this review is on how results from first-principles calculations and classical techniques, such as Monte Carlo and molecular dynamics simulations, have been integrated into multiscale theoretical predictions useful in understanding shape-selective nanocrystal syntheses. Also, examples are discussed in which machine learning has been useful in this field. There are many areas at the frontier in condensed matter theory and simulation that are or could be beneficial in this area and these prospects for future progress are discussed.

Using micro-solvation and generalized coordination numbers to estimate the solvation energies of adsorbed hydroxyl on metal nanoparticles

Solvent-adsorbate interactions have a great impact on catalytic processes in aqueous systems. Implicit solvent calculations are inexpensive but inaccurate toward hydrogen bonds, while a full incorporation of explicit solvation is computationally demanding. Micro-solvation attempts to break this dilemma by including only those solvent molecules directly interacting with the solute and any nearby interfaces, thereby providing a compromise between accuracy and computational expenses. Here, we show that micro-solvation of *OH and its relation to adsorption sites is largely transferable across late transition metal nanoparticles. Solvation energies for *OH on nanoparticles of Ir, Pd, and Pt range from -0.63 ± 0.04 eV to -0.67 ± 0.12 eV, while those on Au and Ag are -0.75 ± 0.07 eV and -1.01 ± 0.05 eV, respectively. These results enable the use of average solvation corrections for *OH on late transition metal nanostructures.

Beyond single-crystal surfaces: The GAL21 water/metal force field

Solvent effects are notoriously difficult to describe for metallic nanoparticles (NPs). Here, we introduce GAL21 which is the first pairwise additive force field that is specifically designed to modulate the near chemisorption energy of water as a function of the coordination numbers of the metallic atoms. We find a quadratic dependence to be most suitable for capturing the dependence of the adsorption energy of water on the generalized coordination number (GCN) of the metal atoms. GAL21 has been fitted against DFT adsorption energies for Cu, Ag, Au, Ni, Pd, Pt, and Co on 500 configurations and validated on about 3000 configurations for each metal, constructed on five surfaces with GCNs varying from 2.5 to 11.25. Depending on the metals, the root mean square deviation is found between 0.7 kcal mol−1 (Au) to 1.6 kcal mol−1 (Ni). Using GAL21, as implemented in the open-source code CP2K, we then evaluate the solvation energy of Au55 and Pt55 NPs in water using thermodynamic integration. The solvation free energy is found to be larger for Pt than for Au and systematically larger than 200 kcal mol−1, demonstrating the large impact of solvent on the surface energetics of NPs. Still, given that the amorphous NPs are both, the most stable and the most solvated ones, we do not predict a change in the preferred morphology between the gas-phase and in water. Finally, based on a linear regression on three sizes of NPs (from 38 to 147), the solvation energy for Au and Pt surface atoms is found to be −5.2 and −9.9 kcal mol−1, respectively.

A computational tool to accurately and quickly predict 19F NMR chemical shifts of molecules with fluorine-carbon and fluorine-boron bonds

We report the evaluation of density-functional-theory (DFT) based procedures for predicting ¹⁹F NMR chemical shifts at modest computational cost for a range of molecules with fluorine bonds, to be used as a tool for assisting the characterisation of reaction intermediates and products and as an aid to identifying mechanistic pathways. The results for a balanced learning set of molecules were then checked using two further testing sets, resulting in the recommendation of the ωB97XD/aug-cc-pvdz DFT method and basis set as having the best combination of accuracy and computational time, with a RMS error of 3.57 ppm. Cationic molecules calculated without counter-anion showed normal errors, whilst anionic molecules showed somewhat larger errors. The method was applied to the prediction of the conformationally averaged ¹⁹F chemical shifts of 2,2,3,3,4,4,5,5-octafluoropentan-1-ol, in which gauche stereoelectronic effects involving fluorine dominate and to determining the position of coordination equilibria of fluorinated boranes as an aid to verifying the relative energies of intermediate species involved in catalytic amidation reactions involving boron catalysts.

Realistic Modelling of Dynamics at Nanostructured Interfaces Relevant to Heterogeneous Catalysis

The focus of this short review is directed towards investigations of the dynamics of nanostructured metallic heterogeneous catalysts and the evolution of interfaces during reaction—namely, the metal–gas, metal–liquid, and metal–support interfaces. Indeed, it is of considerable interest to know how a metal catalyst surface responds to gas or liquid adsorption under reaction conditions, and how its structure and catalytic properties evolve as a function of its interaction with the support. This short review aims to offer the reader a birds-eye view of state-of-the-art methods that enable more realistic simulation of dynamical phenomena at nanostructured interfaces by exploiting resource-efficient methods and/or the development of computational hardware and software.

Revealing Size Dependent Structural Transitions in Supported Gold Nanoparticles in Hydrogen at Atmospheric Pressure

The enhancement of the catalytic activity of gold nanoparticles with their decreasing size is often attributed to the increasing proportion of low-coordinated surface sites. This correlation is based on the paradigmatic picture of working gold nanoparticles as perfect crystal forms having complete and static outer surface layers whatever their size. This picture is incomplete as catalysts can dynamically change their structure according to the reaction conditions and as such changes can be eventually size-dependent. In this work, using aberration-corrected environmental electron microscopy, size-dependent crystal structure and morphological evolution in gold nanoparticles exposed to hydrogen at atmospheric pressure, with loss of the face-centered cubic crystal structure of gold for particle size below 4 nm, are revealed for the first time. Theoretical calculations highlight the role of mobile gold atoms in the observed symmetry changes and particle reshaping in the critical size regime. An unprecedented stable surface molecular structure of hydrogenated gold decorating a highly distorted core is identified. By combining atomic scale in situ observations and modeling of nanoparticle structure under relevant reaction conditions, this work provides a fundamental understanding of the size-dependent reactivity of gold nanoparticles with a precise picture of their surface at working conditions.

Effect of the Ligand Binding Strength on the Morphology of Functionalized Gold Nanoparticles

Functionalized gold nanoparticles are investigated by density functional theory calculations in the context of cancer radiotherapy. Several typical experimental shapes, including nanostars, nanospheres, and nanorods, are modeled by optimizing Au clusters covered by organic monolayers composed of hydrated short-chain polyethylene glycol (PEG) ligands. The PEGylation stabilizes significantly the stellation of decahedral Au₅₄ by deforming significantly its geometry at the spikes. The higher stability of the PEG molecules adsorbed on this stellated nanocluster with respect to the more spherical icosahedral Au₅₅ and truncated octahedral Au₇₉ leads to a larger energy cost to desorb them and thus a weaker propensity for the starred nanoparticle to exchange ligands with the cell membrane, in agreement with experiments. These results open interesting possibilities for advancing our understanding of the cellular uptake of gold nanoparticles.