Making copper, silver and gold fullerene cages breathe.

We show that optical properties change when the fullerene structures of Au32, Cu32and Ag32inflate and deflate. We first observe significant differences in the extinction spectra employing a classical approach based on the Green's dyadic method. By means of real-time time-dependent density functional theory. We continue to calculate the optical spectrum (OP) via aδ-kick simulation, comparing results with the ground-state energetic property the HOMO-LUMO (HL) gap. Red-shift of the OP is expected as the fullerenes inflate, with only ±10% change in the size. As the fullerene breathes, a 0.8 eV shift in the first peak position could be observed in the gold nanoparticle. Ag has a smoother behaviour than both Au and Cu. We have also found changes in the optical spectra can not be directly interpreted as a result of changes in the HL gap.

A kinetic Monte Carlo-blueprint for oxygen reduction on oxide-supported PtNi nanoalloys.

To elucidate the effect of the architecture of supported bimetallic nanocatalysts, we developed a new lattice kinetic Monte Carlo based on the classifying and counting adsorption sites with respect to their generalized coordination number. We employed this tool to estimate the activity of MgO-supported PtNi nanoalloys for oxygen reduction. We demonstrated that the presence of Ni atoms in contact with the substrate massively enhances their activity with at least a 7-order of magnitude increase in the turnover of water production with respect to the case where only Pt lay at the interface. We further discussed how the nanoalloy shape affects the activity showing that truncated octahedra are 102 more active than cuboctahedra of similar size. We explained our results in terms of their distinct distribution and occurrence of the most active sites for oxygen reduction leading to the stabilization of different chemical species during the reaction dynamics. Our results suggest that engineering multifaceted and long edge PtNi-nanoalloys with a certain chemical ordering at the support interface would improve their catalytic performance for the oxygen reduction reaction.

Author Correction: Thermodynamics of CuPt nanoalloys.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Thermodynamics of CuPt nanoalloys.

The control of structural and chemical transitions in bimetallic nanoalloys at finite temperatures is one of the challenges for their use in advanced applications. Comparing Nested Sampling and Molecular Dynamics simulations, we investigate the phase changes of CuPt nanoalloys with the aim to elucidate the role of kinetic effects during their solidification and melting processes. We find that the quasi-thermodynamic limit for the nucleation of (CuPt)309 is 965 ± 10 K, but its prediction is increasingly underestimated when the system is cooled faster than 109 K/s. The solidified nanoparticles, classified following a novel tool based on Steinhardt parameters and the relative orientation of characteristic atomic environments, are then heated back to their liquid phase. We demonstrate the kinetic origin of the hysteresis in the caloric curve as (i) it closes for rates slower than 108 K/s, with a phase change temperature of 970 K ± 25 K, in very good agreement with its quasi-thermodynamic limit; (ii) the process happens simultaneously in the inner and outer layers; (iii) an onion-shell chemical order - Cu-rich surface, Pt-rich sub-surface, and mixed core - is always preserved.

Structure, thermodynamics, and rearrangement mechanisms in gold clusters-insights from the energy landscapes framework.

We consider finite-size and temperature effects on the structure of model AuN clusters (30 ≤ N ≤ 147) bound by the Gupta potential. Equilibrium behaviour is examined in the harmonic superposition approximation, and the size-dependent melting temperature is also bracketed using molecular dynamics simulations. We identify structural transitions between distinctly different morphologies, characterised by various defect features. Reentrant behaviour and trends with respect to cluster size and temperature are discussed in detail. For N = 55, 85, and 147 we visualise the topography of the underlying potential energy landscape using disconnectivity graphs, colour-coded by the cluster morphology; and we use discrete path sampling to characterise the rearrangement mechanisms between competing structures separated by high energy barriers (up to 1 eV). The fastest transition pathways generally involve metastable states with multiple fivefold disclinations and/or a high degree of amorphisation, indicative of melting. For N = 55 we find that reoptimising low-lying minima using density functional theory (DFT) alters their energetic ordering and produces a new putative global minimum at the DFT level; however, the equilibrium structure predicted by the Gupta potential at room temperature is consistent with previous experiments.

Melting of large Pt@MgO(1 0 0) icosahedra.

On the basis of ab initio calculations, we present a new parametrisation of the Vervisch-Mottet-Goniakowski (VMG) potential (Vervisch et al 2002 Phys. Rev. B 24 245411) for modelling the oxide-metal interaction. Applying this model to mimic the finite temperature behaviour of large platinum icosahedra deposited on the pristine MgO(1 0 0), we find the nanoparticle undergoes two solid-solid transitions. At 650 K the 'squarisation' of the interface layer, while a full reshaping towards a fcc architecture takes place above 950 K. In between, a quite long-lived intermediate state with a (1 0 0) interface but with an icosahedral cap is observed. Our approach reproduces experimental observations, including wetting behaviour and the lack of surface diffusion.

Metallic nanoparticles meet metadynamics.

Metadynamics coupled with classical molecular dynamics has been successfully applied to sample the configuration space of metallic and bimetallic nanoclusters. We implement a new set of collective variables related to the pair distance distribution function of the nanoparticle to achieve an exhaustive isomer sampling. As paradigmatic examples, we apply our methodology to Ag147, Pt147, and their alloy Ag(shell)Pt(core) at 2:1 and 1:1 chemical compositions. The proposed scheme is able to reproduce the known solid-solid structural transformation pathways, based on the Lipscomb's diamond-square-diamond mechanisms, both in mono and bimetallic nanoparticles. A discussion of the free energy barriers involved in these processes is provided.

Far-infrared absorption of water clusters by first-principles molecular dynamics.

Based on first-principle molecular dynamic simulations, we calculate the far-infrared spectra of small water clusters (H(2)O)(n) (n = 2, 4, 6) at frequencies below 1000 cm(-1) and at 80 K and at atmospheric temperature (T>200 K). We find that cluster size and temperature affect the spectra significantly. The effect of the cluster size is similar to the one reported for confined water. Temperature changes not only the shape of the spectra but also the total strength of the absorption, a consequence of the complete anharmonic nature of the classical dynamics at high temperature. In particular, we find that in the frequency region up to 320 cm(-1), the absorption strength per molecule of the water dimer at 220 K is significantly larger than that of bulk liquid water, while tetramer and hexamer show bulklike strengths. However, the absorption strength of the dimer throughout the far-infrared region is too small to explain the measured vapor absorption continuum, which must therefore be dominated by other mechanisms.

Finite-temperature effects on the stability and infrared spectra of HCl(H2O)6 clusters.

Theoretical studies of the interaction of HCl with small water clusters have so far neglected the effect of temperature, which ranges from a few tens of kelvin in cluster experiments, up to about 250 K in typical atmospheric conditions. We study the dynamical behavior of a selected set of HCl(H2O)6 clusters, representative of undissociated and dissociated configurations, by means of DFT-based first principles molecular dynamics. We find that the thermodynamcal stability of different configurations can be affected by temperature. We also present the infrared spectra of dissociated and undissociated configurations at 200 K and discuss the origin of the spectral features.

Single impurity effect on the melting of nanoclusters.

We show by molecular dynamics simulations that the melting temperature of clusters can be tuned by selective doping. In fact, a single Ni or Cu impurity in Ag icosahedral clusters considerably increases the melting temperature even for sizes of more than a hundred atoms. The upward shift is correlated to the strain relaxation induced by a small central impurity in icosahedral clusters.

Magic polyicosahedral core-shell clusters.

A new family of magic cluster structures is found by genetic global optimization, whose results are confirmed by density functional calculations. These clusters are Ag-Ni and Ag-Cu nanoparticles with an inner Ni or Cu core and an Ag external shell, as experimentally observed for Ag-Ni, and present a polyicosahedral character. The interplay of the core-shell chemical ordering with the polyicosahedral structural arrangement gives high-symmetry clusters of remarkable structural, thermodynamic, and electronic stability, which can have high melting points (they melt higher than pure clusters of the same size), large energy gaps, and (in the case of Ag-Ni) nonzero magnetic moments.

Amorphization mechanism of icosahedral metal nanoclusters.

The amorphization mechanism of icosahedral Pt nanoclusters is investigated by a combination of molecular dynamics simulations and density functional calculations. A general mechanism for amorphization, involving rosettelike structural transformations at fivefold vertices, is proposed. In the rosette, a fivefold vertex is transformed into a hexagonal ring. We show that, for icosahedral Pt nanoclusters, this transformation is associated with an energy gain, so that their most favorable structures have a low symmetry even at icosahedral magic numbers, and that the same mechanism underlies the formation of amorphous structures in gold.

Growth of three-shell onionlike bimetallic nanoparticles.

We show by molecular dynamics simulations on three systems (B/A=Pd/Ag, Cu/Ag, and Ni/Ag) that three-shell metallic nanoparticles made by a core of a metal A, an intermediate shell of metal B and an external shell of metal A (A-B-A nanoparticles) can be grown by deposition of B atoms onto an A core. The growth of the intermediate B shell is triggered by the fact that the most favorable positions for isolated B impurities inside A clusters are located just one layer below the cluster surface.

Evidence of kinetic trapping in clusters of C60 molecules.

Molecular dynamics simulations of the growth of nanoclusters of C60 provide convincing evidence that experimental magic numbers, which are associated with structures based on Mackay icosahedra, are of kinetic origin. This finding resolves a long-standing contradiction between the experimental observations and the theoretically predicted most stable structures. Our results show that, even if a sticky intermolecular potential energetically disfavors icosahedral structures, the latter are frequently produced because the stickiness of the potential itself enhances kinetic trapping phenomena.

Reentrant morphology transition in the growth of free silver nanoclusters

The growth of free silver nanoclusters is studied by molecular dynamics simulations, from a small seed up to sizes N approximately 150. It is shown that the final outcome of the growth process depends crucially on the growth conditions (deposition flux straight phi and temperature T). A reentrant morphology transition is obtained: at intermediate values of T and straight phi a "decahedral window" is found; the window is surrounded by regimes where icosahedra are preferentially grown.