Hydrogen evolution descriptors of 55-atom PtNi nanoclusters and interaction with graphite.

Density functional simulations have been performed for PtnNi55-nclusters (n=0,12,20,28,42,55) to investigate their catalytic properties for the hydrogen evolution reaction (HER). Starting from the icosahedralPt12Ni43, hydrogen adsorption energetics and electronicd-band descriptors indicate HER activity comparable to that of purePt55(distorted reduced core structure). The PtNi clusters accommodate a large number of adsorbed hydrogen before reaching a saturated coverage, corresponding to 3-4 H atoms per icosahedron facet (in total ∼70-80). The differential adsorption free energies are well within the window of|ΔGH|<0.1 eV which is considered optimal for HER. The electronic descriptors show similarities with the platinumd-band, although the uncovered PtNi clusters are magnetic. Increasing hydrogen coverage suppresses magnetism and depletes electron density, resulting in expansion of the PtNi clusters. For a single H atom, the adsorption free energy varies between -0.32 (Pt12Ni43) and -0.59 eV (Pt55). The most stable adsorption site is Pt-Pt bridge for Pt-rich compositions and a hollow site surrounded by three Ni for Pt-poor compositions. A hydrogen molecule dissociates spontaneously on the Pt-rich clusters. The above HER activity predictions can be extended to PtNi on carbon support as the interaction with a graphite model structure (w/o vacancy defect) results in minor changes in the cluster properties only. The cluster-surface interaction is the strongest forPt55due to its large facing facet and associated van der Waals forces.

Three in One: Three Different Molybdates Trapped in a Thiacalix[4]arene Protected Ag72 Nanocluster for Structural Transformation and Photothermal Conversion.

Polyoxometalates (POMs) represent crucial intermediates in the formation of insoluble metal oxides from soluble metal ions, however, the rapid hydrolysis-condensation kinetics of MoVI or WVI makes the direct characterization of coexisted molecular species in a given medium extremely difficult. Silver nanoclusters have shown versatile capacity to encapsulate diverse POMs, which provides an alternative scene to appreciate landscape of POMs in atomic precision. Here, we report a thiacalix[4]arene protected silver nanocluster (Ag72b) that simultaneously encapsulates three kinds of molybdates (MoO4 2- , Mo6 O22 8- and Mo7 O25 8- ) in situ transformed from classic Lindqvist Mo6 O19 2- , providing more deep understanding on the structural diversity and condensation growth route of POMs in solution. Ag72b is the first silver nanocluster trapping so many kinds of molybdates, which in turn exert collective template effect to aggregate silver atoms into a nanocluster. The post-reaction of Ag72b with AgOAc or PhCOOAg produces a discrete Ag24 nanocluster (Ag24a) or an Ag28 nanocluster based 1D chain structure (Ag28a), respectively. Moreover, the post-synthesized Ag28a can be utilized as potential ignition material for further application. This work not only provides an important model for unlocking dynamic features of POMs at atom-precise level but also pioneers a promising approach to synthesize silver nanoclusters from known to unknown.

Sulfur-deficient edges as active sites for hydrogen evolution on MoS2.

A grand-canonical approach is employed to calculate the voltage-dependent activation energy and estimate the kinetics of the hydrogen evolution reaction (HER) on intrinsic sites of MoS2, including edges of varying S-coverage as well as S-vacancies on the basal plane. Certain edge configurations are found to be vastly more active than others, namely S-deficient edges on the Mo-termination where, in the fully S-depleted case, HER can proceed with activation energy below 0.5 eV at an electrode potential of 0 V vs. SHE. There is a clear distinction between the performance of Mo-rich and S-rich adsorption sites, as HER at the latter sites is characterized by large (generally above 1.5 eV) Heyrovsky and Tafel energy barriers despite near-thermoneutral hydrogen adsorption energy. Thus, exposing Mo-atoms on the edges to which hydrogen can directly bind is crucial for efficient hydrogen evolution. While S-vacancies on the basal plane do expose Mo-rich sites, the energy barriers are still significant due to high coordination of the Mo atoms. Kinetic modelling based on the voltage-dependent reaction energetics gives a theoretical overpotential of 0.25 V and 1.09 V for the Mo-edge with no S atoms and the weakly sulfur-deficient (2% S-vacancies) basal plane, respectively, with Volmer-Heyrovsky being the dominant pathway. These values coincide well with reported experimentally measured values of the overpotential for the edges and basal plane. For the partly Mo-exposed edges, the calculated overpotential is 0.6-0.7 V while edges with only S-sites give overpotential exceeding that of the basal plane. These results show that the overpotential systematically decreases with increased sulfur-deficiency and reduced Mo-coordination. The fundamental difference between Mo- and S-rich sites suggests that catalyst design of transition metal dichalcogenides should be focused on facilitating and modifying the metal sites, rather than activating the chalcogen sites.

Exceptional Microscale Plasticity in Amorphous Aluminum Oxide at Room Temperature.

Oxide glasses are an elementary group of materials in modern society, but brittleness limits their wider usability at room temperature. As an exception to the rule, amorphous aluminum oxide (a-Al2 O3 ) is a rare diatomic glassy material exhibiting significant nanoscale plasticity at room temperature. Here, it is shown experimentally that the room temperature plasticity of a-Al2 O3 extends to the microscale and high strain rates using in situ micropillar compression. All tested a-Al2 O3 micropillars deform without fracture at up to 50% strain via a combined mechanism of viscous creep and shear band slip propagation. Large-scale molecular dynamics simulations align with the main experimental observations and verify the plasticity mechanism at the atomic scale. The experimental strain rates reach magnitudes typical for impact loading scenarios, such as hammer forging, with strain rates up to the order of 1 000 s-1 , and the total a-Al2 O3 sample volume exhibiting significant low-temperature plasticity without fracture is expanded by 5 orders of magnitude from previous observations. The discovery is consistent with the theoretical prediction that the plasticity observed in a-Al2 O3 can extend to macroscopic bulk scale and suggests that amorphous oxides show significant potential to be used as light, high-strength, and damage-tolerant engineering materials.

Is the doped MoS2 basal plane an efficient hydrogen evolution catalyst? Calculations of voltage-dependent activation energy.

Transition metal dichalcogenides are cheap and earth-abundant candidates for the replacement of precious metals as catalyst materials. Experimental measurements of the hydrogen evolution reaction (HER), for example, have demonstrated significant electrocatalytic activity of MoS2 but there is large variation depending on the preparation method. In order to gain information about the mechanism and active sites for the HER, we have carried out calculations of the reaction and activation energy for HER at the transition metal doped basal plane of MoS2 under electrochemical conditions, i.e. including applied electrode potential and solvent effects. The calculations are based on identifying the relevant saddle points on the energy surface obtained from density functional theory within the generalized gradient approximation, and the information on energetics is used to construct voltage-dependent volcano plots. Doping with 3d-metal atoms as well as Pt is found to enhance hydrogen adsorption onto the basal plane by introducing electronic states within the band gap, and in some cases (Co, Ni, Cu, Pt) significant local symmetry breaking. The Volmer-Heyrovsky mechanism is found to be most likely and the associated energetics show considerable dopant and voltage-dependence. While the binding free energy of hydrogen can be tuned to be seemingly favorable for HER, the calculated activation energy turns out to be significant, at least 0.7 eV at a voltage of -0.5 V vs. SHE, indicating low catalytic activity of the doped basal plane. This suggests that other sites are responsible for the experimental activity, possibly edges or basal plane defects.

Machine-learned model Hamiltonian and strength of spin-orbit interaction in strained Mg2X (X = Si, Ge, Sn, Pb).

Machine-learned multi-orbital tight-binding (MMTB) Hamiltonian models have been developed to describe the electronic characteristics of intermetallic compounds Mg2Si, Mg2Ge, Mg2Sn, and Mg2Pb subject to strain. The MMTB models incorporate spin-orbital mediated interactions and they are calibrated to the electronic band structures calculated via density functional theory by a massively parallelized multi-dimensional Monte-Carlo search algorithm. The results show that a machine-learned five-band tight-binding (TB) model reproduces the key aspects of the valence band structures in the entire Brillouin zone. The five-band model reveals that compressive strain localizes the contribution of the 3sorbital of Mg to the conduction bands and the outer shellporbitals of X (X = Si, Ge, Sn, Pb) to the valence bands. In contrast, tensile strain has a reversed effect as it weakens the contribution of the 3sorbital of Mg and the outer shellporbitals of X to the conduction bands and valence bands, respectively. Theπbonding in the Mg2X compounds is negligible compared to theσbonding components, which follow the hierarchy|σsp|>|σpp|>|σss|, and the largest variation against strain belongs toσpp. The five-band model allows for estimating the strength of spin-orbit coupling (SOC) in Mg2X and obtaining its dependence on the atomic number of X and strain. Further, the band structure calculations demonstrate a significant band gap tuning and band splitting due to strain. A compressive strain of-10%can open a band gap at the Γ point in metallic Mg2Pb, whereas a tensile strain of+10%closes the semiconducting band gap of Mg2Si. A tensile strain of+5%removes the three-fold degeneracy of valence bands at the Γ point in semiconducting Mg2Ge. The presented MMTB models can be extended for various materials and simulations (band structure, transport, classical molecular dynamics), and the obtained results can help in designing devices made of Mg2X.

Publisher Correction: Conductivity control via minimally invasive anti-Frenkel defects in a functional oxide.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Conductivity control via minimally invasive anti-Frenkel defects in a functional oxide.

Utilizing quantum effects in complex oxides, such as magnetism, multiferroicity and superconductivity, requires atomic-level control of the material's structure and composition. In contrast, the continuous conductivity changes that enable artificial oxide-based synapses and multiconfigurational devices are driven by redox reactions and domain reconfigurations, which entail long-range ionic migration and changes in stoichiometry or structure. Although both concepts hold great technological potential, combined applications seem difficult due to the mutually exclusive requirements. Here we demonstrate a route to overcome this limitation by controlling the conductivity in the functional oxide hexagonal Er(Mn,Ti)O3 by using conductive atomic force microscopy to generate electric-field induced anti-Frenkel defects, that is, charge-neutral interstitial-vacancy pairs. These defects are generated with nanoscale spatial precision to locally enhance the electronic hopping conductivity by orders of magnitude without disturbing the ferroelectric order. We explain the non-volatile effects using density functional theory and discuss its universality, suggesting an alternative dimension to functional oxides and the development of multifunctional devices for next-generation nanotechnology.

Highly ductile amorphous oxide at room temperature and high strain rate.

Oxide glasses are an integral part of the modern world, but their usefulness can be limited by their characteristic brittleness at room temperature. We show that amorphous aluminum oxide can permanently deform without fracture at room temperature and high strain rate by a viscous creep mechanism. These thin-films can reach flow stress at room temperature and can flow plastically up to a total elongation of 100%, provided that the material is dense and free of geometrical flaws. Our study demonstrates a much higher ductility for an amorphous oxide at low temperature than previous observations. This discovery may facilitate the realization of damage-tolerant glass materials that contribute in new ways, with the potential to improve the mechanical resistance and reliability of applications such as electronic devices and batteries.

Varying oxygen coverage on Cu55 and its effect on CO oxidation.

Adsorption of molecular oxygen on a Cu55 cluster and the resulting oxidation effects have been investigated by spin-polarized density functional theory (DFT). The optimal structure for each Cu55O2N (N = 1-20) complex has been obtained via a sequential addition of O2 and systematic screening of the preferable adsorption sites. Upon structural optimization, several O2 molecules dissociate readily on Cu55 at different oxygen coverages, and further DFT molecular dynamics simulations at 300 K confirm the instability (small dissociation barrier) of the remaining O2 and a spontaneous movement of some oxygen atoms from the surface sites towards the cluster interior. The Cu55 cluster and its oxidized derivatives have been placed on a γ-Al2O3(100) surface to study the cluster-support interaction, and furthermore, CO oxidation reactions on both Cu55(O)2N and Cu55(O)2N/γ-Al2O3(100) have been studied as a function of oxygen coverage. The CO oxidation reaction barrier is rather insensitive to the oxygen coverage regardless of the support, indicating a small increase in activity with the number of surface oxygen atoms.

CLEASE: a versatile and user-friendly implementation of cluster expansion method.

Materials exhibiting a substitutional disorder such as multicomponent alloys and mixed metal oxides/oxyfluorides are of great importance in many scientific and technological sectors. Disordered materials constitute an overwhelmingly large configurational space, which makes it practically impossible to be explored manually using first-principles calculations such as density functional theory due to the high computational costs. Consequently, the use of methods such as cluster expansion (CE) is vital in enhancing our understanding of the disordered materials. CE dramatically reduces the computational cost by mapping the first-principles calculation results on to a Hamiltonian which is much faster to evaluate. In this work, we present our implementation of the CE method, which is integrated as a part of the atomic simulation environment (ASE) open-source package. The versatile and user-friendly code automates the complex set up and construction procedure of CE while giving the users the flexibility to tweak the settings and to import their own structures and previous calculation results. Recent advancements such as regularization techniques from machine learning are implemented in the developed code. The code allows the users to construct CE on any bulk lattice structure, which makes it useful for a wide range of applications involving complex materials. We demonstrate the capabilities of our implementation by analyzing the two example materials with varying complexities: a binary metal alloy and a disordered lithium chromium oxyfluoride.

The Role of Temperature and Lipid Charge on Intake/Uptake of Cationic Gold Nanoparticles into Lipid Bilayers.

Understanding the molecular mechanisms governing nanoparticle-membrane interactions is of prime importance for drug delivery and biomedical applications. Neutron reflectometry (NR) experiments are combined with atomistic and coarse-grained molecular dynamics (MD) simulations to study the interaction between cationic gold nanoparticles (AuNPs) and model lipid membranes composed of a mixture of zwitterionic di-stearoyl-phosphatidylcholine (DSPC) and anionic di-stearoyl-phosphatidylglycerol (DSPG). MD simulations show that the interaction between AuNPs and a pure DSPC lipid bilayer is modulated by a free energy barrier. This can be overcome by increasing temperature, which promotes an irreversible AuNP incorporation into the lipid bilayer. NR experiments confirm the encapsulation of the AuNPs within the lipid bilayer at temperatures around 55 °C. In contrast, the AuNP adsorption is weak and impaired by heating for a DSPC-DSPG (3:1) lipid bilayer. These results demonstrate that both the lipid charge and the temperature play pivotal roles in AuNP-membrane interactions. Furthermore, NR experiments indicate that the (negative) DSPG lipids are associated with lipid extraction upon AuNP adsorption, which is confirmed by coarse-grained MD simulations as a lipid-crawling effect driving further AuNP aggregation. Overall, the obtained detailed molecular view of the interaction mechanisms sheds light on AuNP incorporation and membrane destabilization.

Molecular-Scale Ligand Effects in Small Gold-Thiolate Nanoclusters.

Because of the small size and large surface area of thiolate-protected Au nanoclusters (NCs), the protecting ligands are expected to play a substantial role in modulating the structure and properties, particularly in the solution phase. However, little is known on how thiolate ligands explicitly modulate the structural properties of the NCs at atomic level, even though this information is critical for predicting the performance of Au NCs in application settings including as a catalyst interacting with small molecules and as a sensor interacting with biomolecular systems. Here, we report a combined experimental and theoretical study, using synchrotron X-ray spectroscopy and quantum mechanics/molecular mechanics simulations, that investigates how the protecting ligands impact the structure and properties of small Au18(SR)14 NCs. Two representative ligand types, smaller aliphatic cyclohexanethiolate and larger hydrophilic glutathione, are selected, and their structures are followed experimentally in both solid and solution phases. It was found that cyclohexanethiolate ligands are significantly perturbed by toluene solvent molecules, resulting in structural changes that cause disorder on the surface of Au18(SR)14 NCs. In particular, large surface cavities in the ligand shell are created by interactions between toluene and cyclohexanethiolate. The appearance of these small molecule-accessible sites on the  NC surface demonstrates the ability of Au NCs to act as a catalyst for organic phase reactions. In contrast, glutathione ligands encapsulate the Au NC core via intermolecular interactions, minimizing structural changes caused by interactions with water molecules. The much better protection from glutathione ligands imparts a rigidified surface and ligand structure, making the NCs desirable for biomedical applications due to the high stability and also offering a structural-based explanation for the enhanced photoluminescence often reported for glutathione-protected Au NCs.

The molecular mechanism of the ligand exchange reaction of an antibody against a glutathione-coated gold cluster.

The labeling of proteins with heavy atom clusters is of paramount importance in biomedical research, but its detailed molecular mechanism remains unknown. Here we uncover it for the particular case of the anti-influenza N9 neuraminidase NC10 antibody against a glutathione-coated gold cluster by means of ab initio QM/MM calculations. We show that the labeling reaction follows an associative double SN2-like reaction mechanism, involving a proton transfer, with low activation barriers only if one of the two distinct peptide/peptidic ligands (the one that occupies the side position) is substituted. Positively charged residues in the vicinity of the incoming thiol result in strong interactions between the antibody and the AuMPC, favoring the ligand exchange reaction for suitable protein mutants. These results pave the way for future investigations aimed at engineering biomolecules to increase their reactivity towards a desired gold atom cluster.

Collective excitations and viscosity in liquid Bi.

The analysis of extensive density functional/molecular dynamics simulations (over 500 atoms, up to 100 ps) of liquid bismuth at four temperatures between 573 K and 1023 K has provided details of the dynamical structure factors, the dispersion of longitudinal and transverse collective modes, and related properties (power spectrum, viscosity, and sound velocity). Agreement with available inelastic x-ray and neutron scattering data and with previous simulations is generally very good. The results show that density functional/molecular dynamics simulations can give dynamical information of good quality without the use of fitting functions, even at long wavelengths.

Coexisting Honeycomb and Kagome Characteristics in the Electronic Band Structure of Molecular Graphene.

We uncover the electronic structure of molecular graphene produced by adsorbed CO molecules on a copper (111) surface by means of first-principles calculations. Our results show that the band structure is fundamentally different from that of conventional graphene, and the unique features of the electronic states arise from coexisting honeycomb and Kagome symmetries. Furthermore, the Dirac cone does not appear at the K-point but at the Γ-point in the reciprocal space and is accompanied by a third, almost flat band. Calculations of the surface structure with Kekulé distortion show a gap opening at the Dirac point in agreement with experiments. Simple tight-binding models are used to support the first-principles results and to explain the physical characteristics behind the electronic band structures.

DFT simulations and microkinetic modelling of 1-pentyne hydrogenation on Cu20 model catalysts.

Adsorption and dissociation of H2 and hydrogenation of 1-pentyne on neutral and anionic Cu20 clusters have been investigated using the density functional theory and microkinetic modelling. Molecular adsorption of H2 is found to occur strictly at atop sites. The H2 dimer is activated upon adsorption, and the dissociation occurs with moderate energy barriers. The dissociated H atoms reside preferentially on 3-fold face and 2-fold edge sites. Based on these results, the reaction paths leading to the partial and total hydrogenation of 1-pentyne have been studied step-by-step. The results suggest that copper clusters can display selective activity on the hydrogenation of alkyne and alkene molecules. The hydrogenated products are more stable than the corresponding initial reactants following an energetic staircase with the number of added H atoms. Stable semi-hydrogenated intermediates are formed before the partial (1-pentene) and total (pentane) hydrogenation stages of 1-pentyne. The microkinetic model analysis shows that C5H10 is the dominant product. Increasing the reactants (C5H8/H2) ratio enhances the formation of products (C5H10 and C5H12).

The electrooxidation-induced structural changes of gold di-superatomic molecules: Au23vs. Au25.

The gold cluster compounds Au38(SC2H4Ph)24 and [Au25(PPh3)10(SC2H4Ph)5Cl2](2+) are known to possess bi-icosahedral Au23 and Au25 cores, respectively, inside their ligand shells. These Au cores can be viewed as quasi-molecules composed of two Au13 superatoms sharing three and one Au(+) atoms, respectively. In the present work, we studied the structural changes of these gold di-superatomic molecules upon electrooxidation via spectroelectrochemical techniques, X-ray absorption fine structure analysis, and density functional theory calculations. The Au23 core was electrochemically stable, but the Au25 core underwent irreversible structural change. This marked difference in the stability of the oxidized states is ascribed to differences in the bonding scheme of Au13 units and/or the bonding nature of the protecting ligands.

Simple metal under tensile stress: layer-dependent herringbone reconstruction of thin potassium films on graphite.

While understanding the properties of materials under stress is fundamentally important, designing experiments to probe the effects of large tensile stress is difficult. Here tensile stress is created in thin films of potassium (up to 4 atomic layers) by epitaxial growth on a rigid support, graphite. We find that this "simple" metal shows a long-range, periodic "herringbone" reconstruction, observed in 2- and 3- (but not 1- and 4-) layer films by low-temperature scanning tunneling microscopy (STM). Such a pattern has never been observed in a simple metal. Density functional theory (DFT)simulations indicate that the reconstruction consists of self-aligned stripes of enhanced atom density formed to relieve the tensile strain. At the same time marked layer-dependent charging effects lead to substantial variation in the apparent STM layer heights.