Tuning the electronic structure of gold cluster-assembled materials by altering organophosphine ligands.

Superatomic clusters can be assembled to build bulk matter, where the individual characteristics are preserved. The main benefit of these materials over conventional bulk species is the capability to tailor their features by altering the physicochemical identities of individual clusters. Electronic properties of metal clusters can be modified by a protective shell of ligands that attach to the surface and make the whole nanoparticle soluble in organic or aqueous solvents. In the present work, we demonstrate that properly chosen ligands provide not only steric protection from aggregation but also tune the redox activity of metal clusters. We investigate the role of the ligands in electronic structure tunability and ligand-field splitting. Our first-principles calculations agree with the experiments, showing that phosphine-protected gold materials are small gap semiconductors. The obtained bandgaps strongly depend on the ligand used. Hence, using phosphine and organophosphine ligands should be feasible and promising while designing the novel superatom-based materials since the desired range of the bandgap might be achieved (by the proper choice of the ligand).

An Atomic-scale Explanation For The High Selectivity Towards Carbon Dioxide Reduction Observed On Liquid Metal Catalysts.

The low-temperature liquid metals Ga-In and Ga-Sn have previously showcased >95% selectivity towards the electrochemical reduction of CO2 to formate, occuring only when the alloys are melted, not solid. Here, density functional theory molecular dynamics and metadynamics simulations reveal that CO2 does not directly adsorb to the Ga-alloy surface, but instead is reduced indirectly by reaction with an adsorbed hydrogen. The reaction barrier is vastly more favourable when this process occurs at In or Sn sites (average: 0.26 eV), than when it occurs on Ga (average: 0.47 eV). However, there is no difference in barrier between solid and liquid surfaces. Instead, we find that Hads is mobile only on the liquid surface, travelling due to the motion of the liquid beneath. This process drives Hads to In/Sn sites, allowing low-barrier CO2 reduction to occur only on the liquid. Therefore, the dynamic motion of liquid metal catalysts can underpin their unique reactivity. The result has far reaching implications for any protonation reaction conducted with a liquid metal catalyst.

Structural and electronic changes in Ga-In and Ga-Sn alloys on melting.

The melting behaviour of surface slabs of Ga-In and Ga-Sn is studied using periodic density functional theory and ab initio molecular dynamics. Analysis of the structure and electronics of the solid and liquid phases gives insight into the properties of these alloys, and why they may act as promising CO2 reduction catalysts. We report melting points for slabs of hexa-layer Ga-In (386 K) and Ga-Sn (349 K) that are substantially lower than the pure hexa-layer Ga system (433 K), and attribute the difference to the degree to which the dopant (In or Sn) disrupts the layered Ga network. In molecular dynamics trajectories of the liquid structures, we find that dopant tends to migrate from the centre of the slab towards the surface and accumulate there. Bader charge calculations reveal that the surface dopant atoms have increased positive charge, and density of states analyses suggest the liquid alloys maintain metallic electronic behaviour. Thus, surface In and Sn may provide good binding sites for intermediates in CO2 reduction. This work contributes to our understanding of the properties of liquid metal systems, and provides a foundation for modelling catalysis on these materials.

Dynamic Activation of Ga Sites by Pt Dopant in Low-Temperature Liquid-Metal Catalysts.

Liquid GaPt catalysts with Pt concentrations as low as 1×10-4  atomic % have recently been identified as highly active for the oxidation of methanol and pyrogallol under mild reaction conditions. However, almost nothing is known about how liquid state catalysts support these significant improvements in activity. Here, ab initio molecular dynamics simulations are employed to examine GaPt catalysts in isolation and interacting with adsorbates. We find that persistent geometric features can exist in the liquid state, given the correct environment. We postulate that the Pt dopant may not be limited to direct involvement in catalysis of reactions, but rather that its presence can also enable Ga atoms to become catalytically active.

Catalytic Potential of Post-Transition Metal Doped Graphene-Based Single-Atom Catalysts for the CO2 Electroreduction Reaction.

Catalysts are required to ensure electrochemical reduction of CO2 to fuels proceeds at industrially acceptable rates and yields. As such, highly active and selective catalysts must be developed. Herein, a density functional theory study of p-block element and noble metal doped graphene-based single-atom catalysts in two defect sites for the electrochemical reduction of CO2 to CO and HCOOH is systematically undertaken. It is found that on all of the systems considered, the thermodynamic product is HCOOH. Pb/C3 , Pb/N4 and Sn/C3 are identified as having the lowest overpotential for HCOOH production while Al/C3 , Al/N4 , Au/C3 and Ga/C3 are identified as having the potential to form higher order products due to the strength of binding of adsorbed HCOOH.

Molecular crystals vs. superatomic lattice: a case study with superalkali-superhalogen compounds.

Using a first-principles approach, we study the assembly of atomically-precise cluster solids with atomic precision. The aims are to create binary assemblies of clusters through charge transfer between neutral molecular clusters, and employing intercluster electrostatic attraction as a driving force for co-assembly. We combined pairs of complementary clusters in which one cluster is electron-donating (superalkali) and the other is electron-accepting (superhalogen). From the analysis of the binding energy between superatomic counterparts, charge transfer, and the relative size of the clusters, we analyze the resulting structures as either molecular crystals or superatomic lattices. We demonstrate that the substitution of a single atom can result in minor changes to the crystal structure of the binary solids or entirely new packing structures. The [N4Mg6Li]+[AlCl4]-, [N4Mg6Na]+[AlCl4]-, [N4Mg6K]+[AlCl4]-, [N4Mg6Li]+[AlF4]-, [N4Mg6Na]+[AlF4]-, and [N4Mg6K]+[AlF4]- compounds all form the same close-packed superatomic lattice structure through halogen bonding, with subtle differences in the orientation of the superatoms. These salts may also form molecular crystals where clusters are held to one another by electrostatic interactions. Our results emphasize how the structure of superatomic solids can be tuned upon single atom substitution.

Structural, Thermal, and Electronic Properties of Two-Dimensional Gallium Oxide (ß-Ga2 O3 ) from First-Principles Design.

Two-dimensional (2D) materials with exotic electronic, optical and mechanical properties have attracted tremendous attention in the last two decades, due to their potential applications in electronics, energy storage and conversion technologies. However, only a few dozen 2D materials have been successfully synthesized or exfoliated. Motivated by the recent discovery of 2D gallenene, we have explored new 2D allotropes of ß-Ga2 O3 , an emerging wide-band gap transparent conductive oxide (TCO) with a wide range of semiconducting applications. All the possible 2D allotropes of ß-Ga2 O3 with high energetic stability have been predicted using particle swarm optimization, combined with density functional theory calculations. The structural and dynamical stability of the predicted 2D allotropes has been analyzed. Although ß-Ga2 O3 is not a van der Waals material, results predict that one or two allotropes of ß-Ga2 O3 are stable. In addition, the accurate band structures of these 2D semiconducting oxides have been calculated using both the GGA and LDA-1/2 approach. Remarkably, monolayer Ga2 O3 (100) has a larger indirect band gap of 4 eV, demonstrating a new avenue for the discovery of 2D ß-Ga2 O3 based nano-devices with enhanced electronic properties.

A mechanistic understanding of surface Bi enrichment in dilute GaBi systems.

Experiment has shown that dilute GaBi systems produce a range of self-organised nanostructured patterns at the surface [Tang et al., Nat. Nanotechnol., 2021, 16, 431-439]. Using extensive ab initio molecular dynamics simulations, we elucidate the mechanisms underlying the formation of the Bi surface islands in Bi-doped Ga liquid metals. Here, we show that in order for internal Bi atoms to diffuse to the surface a lateral extension of the Ga surface network is required. Furthermore, the absence of surface Bi patterning perturbs the Ga surface network providing a preferred path for an internal Bi to diffuse. By understanding how and why Bi nucleates at a surface, we increase the ability to control, manipulate and design such systems for use in future electronic devices.

Clustering of metal dopants in defect sites of graphene-based materials.

Single-atom catalysts are promising candidates for many industrial reactions. However, making true single-atom catalysts is an experimental dilemma, due to the difficulty of keeping dopant single atoms stable at temperature and under pressure. This difficulty can lead to clustering of the metal dopant atoms in defect sites. However, the electronic and geometric structure of sub-nanoscale clusters in single-atom defects has not yet been explored. Furthermore, recent studies have proven sub-nanoscale clusters of dopants in single-atom defect sites can be equally good or better catalysts than their single-atom counterparts. Here, a comprehensive DFT study is undertaken to determine the geometric and electronic structure effects that influence clustering of noble and p-block dopants in C3- and N4-defect sites in graphene-based systems. We find that the defect site is the primary driver in determining clustering dynamics in these systems.

Bimetallic superalkali substitution in the CsPbBr3 perovskite: Pseudocubic phases and tunable bandgap.

Perovskites attract attention as efficient light absorbers for solar cells due to their high-power conversion efficiency (up to 24%). The high photoelectric conversion efficiency is greatly affected by a suitable band structure. Cation substitution can be an effective approach to tune the electronic band structure of lead halide perovskites. In this work, superalkali cations were introduced to replace the Cs+ cation in the CsPbBr3 material. The bimetallic superalkalis (LiMg, NaMg, LiCa, and NaCa) were inserted since they are structurally simple systems and have a strong tendency to lose one electron to achieve a closed-shell cation. The cation substitution in the lead halide perovskite leads to changes in the shape of both valence and conduction bands compared to CsPbBr3. Introducing superalkali cations produces extra electronic states close to the Fermi level, which arise from the formation of alkali earth metal states at the top of the valence band. Our first-principles computations reveal that bimetallic superalkali substitution decreases the bandgap of the perovskite. The bandgaps of MgLi-PbBr3 (1.35 eV) and MgNa-PbBr3 (1.06 eV) are lower than the bandgap of CsPbBr3 (2.48 eV) and within the optimal bandgap (i.e., 1.1-1.4 eV) for single-junction solar cells. Thus, the MgLi-PbBr3 and MgNa-PbBr3 inorganic perovskites are promising candidates for high-efficiency solar cells.

Emergent electronic properties in Co-deposited superatomic clusters.

We report an intercluster compound based on co-deposition of the Au cluster [Au9(PPh3)8](NO3)3 and the fulleride KC60(THF). Electronic properties characteristic for a charge interaction between superatoms emerge within the solid state material [Au9(PPh3)8](NO3)3-x(C60)x, as confirmed by UV-VIS and Raman spectroscopy and I-V measurements. These emergent properties are related to the superatomic electronic states of the initial clusters. The material is characterized by Fourier-transform infrared spectroscopy, x-ray diffraction, Raman spectroscopy, and electrical measurements. Structural optimization and ab initio band structure calculations are performed with density functional theory to interpret the nature of the electronic states in the material; Bader charge calculations assign effective oxidation states in support of the superatomic model of cluster interactions.

Modified Lennard-Jones potentials for nanoscale atoms.

A classical 6-12 Lennard-Jones (LJ) equation has been widely used to model materials and is the potential of choice in studies when the focus is on fundamental issues. Here we report a systematic study comparing the pair interaction potentials within solid-state materials (i.e., [Co6 Se8 (PEt3 )6 ][C60 ]2 , [Cr6 Te8 (PEt3 )6 ][C60 ]2 , [Ni9 Te6 (PEt3 )8 ][C60 ]) using density functional theory (DFT) calculations and LJ parametrization. Both classical (6-12 LJ) and modified LJ (mLJ) models were developed. In the mLJ approach, the exponents 6 and 12 are replaced by different integer number n and 2n, respectively, and an additional parameter (α) is introduced to describe intermolecular distance shift arising within the geometric centers' approach (instead of the shortest interatomic distance between particles). A general LJ approach reexamination reveals that in the case of nanoatoms, the attractive term decays with distance as the inverse fourth power, and the dominating at short distances repulsive term decays as the inverse eighth power. The modification of the LJ equation is even more prominent for interaction profiles, where intermolecular distance corresponds to separation between geometric centers of particles. In this approach, the attractive term decays with distance as the inverse 12th power, while the repulsive term decays rapidly (as the inverse 24th power). Thus, the mLJ models (e.g., 4-8 LJ) rather than the 6-12 classical ones seem to be a better choice for the description of binary interactions of nanoatoms. The developed mLJ models and electronic structure characteristics give an insight into the explanation of the unique physicochemical properties of superatomic-based solid-state materials.

N4Mg6M (M = Li, Na, K) superalkalis for CO2 activation.

Superatoms have exciting properties, including diverse functionalization, redox activity, and magnetic ordering, so the resulting cluster-assembled solids hold the promise of high tunability, atomic precision, and robust architectures. By utilizing adamantane-like clusters as building blocks, a new class of superatoms N4Mg6M (M = Li, Na, K) is proposed here. The studied superalkalis feature low adiabatic ionization energies, an antibonding character in the interactions between magnesium and nitrogen atoms, and highly delocalized highest occupied molecular orbital (HOMO). Consequently, the N4Mg6M superalkalis might easily lose their HOMO electrons when interacting with superhalogen electrophiles to form stable superatom [superalkali]+[superhalogen]- compounds. Moreover, the studied superalkalis interact strongly with carbon dioxide, and the resulting N4Mg6M/CO2 systems represent two strongly interacting ionic fragments (i.e., N4Mg6M+ and CO2 -). In turn, the electron affinity of the N2 molecule (of -1.8 eV) is substantially lower than that observed for carbon dioxide (EA = -0.6 eV) and consequently, the N2 was found to form the weakly bound [N4Mg6M][N2] complex rather than the desired ionic [N4Mg6M]+[N2]- product. Thus, the N4Mg6M superalkalis have high selectivity over N2 when it comes to CO2 reduction and also are themselves stable. We believe that the results described within this paper will be useful for understanding CO2 activation, which is the first step for producing fuels from CO2. Moreover, we demonstrate that designing novel superatomic systems and exploring their physicochemical features might be used to create desirable functional materials.

Aluminum oxo-fluoride clusters: A first principle investigation of stability, synthetic considerations, and the interaction with water.

The introduction of the so called fluorolytic sol-gel synthesis in 2003 gave access to previously inaccessible aluminum oxo-fluorides, thus to nanoscopic materials and, more importantly, novel catalysts. The intermediate cluster structures synthesized and stabilized by Kemnitz and coworkers have mainly been protected by iso-propoxide groups. However, since catalytic reactions take place in a large variety of media, hydrophilic analogs of those clusters would be of interest. In this manuscript, we present a computational analysis for the fluorination reaction, which represents the second part of fluorolytic sol-gel synthesis, and a theoretical study of the synthesized Al4 F4 (µ4 -O)(µ-Oi Pr)5 [H(Oi Pr)2 ] nanostructure's conversion to its hydroxylated analog Al4 F4 (µ4 -O)(µ-OH)5 [H(OH)2 ] utilizing the nudged elastic band method. Furthermore, the role of the fluorine atoms of the cluster in an aqueous medium is evaluated by studying the incremental addition of water molecules to the cluster with and without fluorine atoms. In addition, NMR shifts of clusters exhibiting different substituents are compared. It has been found that the inclusion of an explicit solvent is necessary to capture the magnetic response of the individual cluster atoms in an aqueous solvent correctly. © 2018 Wiley Periodicals, Inc.

Cluster assemblies as superatomic solids: a first principles study of bonding & electronic structure.

The synthesis of cluster based materials poses an exciting challenge for experimental chemistry. The main advantage of these materials compared to conventional bulk compounds is the simple tunability of the chemical and physical characteristics of individual clusters. As a consequence, cluster assemblies can theoretically be used for the creation of designer materials exhibiting specifically desired properties. Since superatoms reveal a large intrinsic thermodynamic stability and often very interesting tunable electronic characteristics, they seem to be an excellent choice as building blocks for the bulk. Here, we present a detailed first principles analysis of carefully chosen superatomic cluster binary and bulk assemblies, in order to determine which forces control the attractive interaction in superatomic solids, and how the individual cluster properties affect these assemblies. This study uses the highly tunable and stable Au13(RS(AuSR)2)6 cluster with a variety of dopants as a model system, while the principles are likely transferable to other ligand protected systems with a straightforward superatomic electron count, such as aluminum or sodium clusters. Three different superatomic materials based on doped gold clusters, boranes and C60s are constructed and evaluated. Beyond the verification that superatoms can be used to create materials that reveal emergent atom-based solid like properties, various factors influencing superatomic materials, such as the EA, IP and relative sizes of the clusters, have been identified and critically evaluated.

How robust is the metallicity of two dimensional gallium?

Atomically thin gallium layers have recently been experimentally produced via solid-melt exfoliation, and show promise as robustly metallic 2D materials for electronic applications. However, the extent to which the experimental technique can be extended to other metals relies on understanding how the 2D structures relate to the bulk form of gallium, which is itself unique as an elemental 'molecular metal'. We relate the experimentally formed 2D materials to the theoretically predicted 'bilayer gallium' which has previously been shown to be stable in vacuum at the nanoscale, via density functional theory calculations. We also study the variation of electronic structure with lattice strain to confirm the extent to which the metallicity will be robust on a wide range of substrate materials.

Building machine learning force fields for nanoclusters.

We assess Gaussian process (GP) regression as a technique to model interatomic forces in metal nanoclusters by analyzing the performance of 2-body, 3-body, and many-body kernel functions on a set of 19-atom Ni cluster structures. We find that 2-body GP kernels fail to provide faithful force estimates, despite succeeding in bulk Ni systems. However, both 3- and many-body kernels predict forces within an ∼0.1 eV/Šaverage error even for small training datasets and achieve high accuracy even on out-of-sample, high temperature structures. While training and testing on the same structure always provide satisfactory accuracy, cross-testing on dissimilar structures leads to higher prediction errors, posing an extrapolation problem. This can be cured using heterogeneous training on databases that contain more than one structure, which results in a good trade-off between versatility and overall accuracy. Starting from a 3-body kernel trained this way, we build an efficient non-parametric 3-body force field that allows accurate prediction of structural properties at finite temperatures, following a newly developed scheme [A. Glielmo et al., Phys. Rev. B 95, 214302 (2017)]. We use this to assess the thermal stability of Ni19 nanoclusters at a fractional cost of full ab initio calculations.

A Two-Dimensional Liquid Structure Explains the Elevated Melting Temperatures of Gallium Nanoclusters.

Melting in finite-sized materials differs in two ways from the solid-liquid phase transition in bulk systems. First, there is an inherent scaling of the melting temperature below that of the bulk, known as melting point depression. Second, at small sizes changes in melting temperature become nonmonotonic and show a size-dependence that is sensitive to the structure of the particle. Melting temperatures that exceed those of the bulk material have been shown to occur for a very limited range of nanoclusters, including gallium, but have still never been ascribed a convincing physical explanation. Here, we analyze the structure of the liquid phase in gallium clusters based on molecular dynamics simulations that reproduce the greater-than-bulk melting behavior observed in experiments. We observe persistent nonspherical shape distortion indicating a stabilization of the surface, which invalidates the paradigm of melting point depression. This shape distortion suggests that the surface acts as a constraint on the liquid state that lowers its entropy relative to that of the bulk liquid and thus raises the melting temperature.

From the Superatom Model to a Diverse Array of Super-Elements: A Systematic Study of Dopant Influence on the Electronic Structure of Thiolate-Protected Gold Clusters.

The electronic properties of doped thiolate-protected gold clusters are often referred to as tunable, but their study to date, conducted at different levels of theory, does not allow a systematic evaluation of this claim. Here, using density functional theory, the applicability of the superatomic model to these clusters is critically evaluated, and related to the degree of structural distortion and electronic inhomogeneity in the differently doped clusters, with dopant atoms Pd, Pt, Cu, and Ag. The effect of electron number is systematically evaluated by varying the charge on the overall cluster, and the nominal number of delocalized electrons, employed in the superatomic model, is compared to the numbers obtained from Bader analysis of individual atomic charges. We find that the superatomic model is highly applicable to all of these clusters, and is able to predict and explain the changing electronic structure as a function of charge. However, significant perturbations of the model arise due to doping, due to distortions of the core structure of the Au13 [RS(AuSR)2 ]6- cluster. In addition, analysis of the electronic structure indicates that the superatomic character is distributed further across the ligand shell in the case of the doped clusters, which may have implications for the self-assembly of these clusters into materials. The prediction of appropriate clusters for such superatomic solids relies critically on such quantitative analysis of the tunability of the electronic structure.

Interaction of Boron Nitride Nanosheets with Model Cell Membranes.

Boron nitride nanomaterials have attracted attention for biomedical applications, due to their improved biocompatibility when compared with carbon nanomaterials. Recently, graphene and graphene oxide nanosheets have been shown, both experimentally and computationally, to destructively extract phospholipids from Escherichia coli. Boron nitride nanosheets (BNNSs) have exciting potential biological and environmental applications, for example the ability to remove oil from water. These applications are likely to increase the exposure of prokaryotes and eukaryotes to BNNSs. Yet, despite their promise, the interaction between BNNSs and cell membranes has not yet been investigated. Here, all-atom molecular dynamics simulations were used to demonstrate that BNNSs are spontaneously attracted to the polar headgroups of the lipid bilayer. The BNNSs do not passively cross the lipid bilayer, most likely due to the large forces experienced by the BNNSs. This study provides insight into the interaction of BNNSs with cell membranes and may aid our understanding of their improved biocompatibility.