Tuning the work function of graphite nanoparticles via edge termination.

Graphite nanoparticles are important in energy materials applications such as lithium-ion batteries (LIBs), supercapacitors and as catalyst supports. Tuning the work function of the nanoparticles allows local control of lithiation behaviour in LIBs, and the potential of zero charge of electrocatalysts and supercapacitors. Using large scale density functional theory (DFT) calculations, we find that the surface termination of multilayer graphene nanoparticles can substantially modify the work function. Calculations in vacuum and in electrolyte show that manipulating the edge termination substantially modifies the potential not only around the edge, but also on the basal plane. Termination with hydrogen or oxygen completely reverses the potential distribution surrounding the basal plane and edges. The trends can be explained based on the work function differences of the edges dependent on termination, and that of the basal plane. Electronic equilibration between different surfaces at the nanoscale allows manipulation of the work function. We demonstrate a link between the area of the graphite basal plane via changing the nanoparticle size, and the work function. We expect that these insights can be utilised for local control of electrochemical functions of graphite nanoparticles prepared under oxidising or reducing conditions.

An Electronic Structure Investigation of PEDOT with AlCl4- Anions-A Promising Redox Combination for Energy Storage Applications.

In this work, we use density functional theory to investigate the electronic structure of poly(3,4-ethylenedioxythiophene) (PEDOT) oligomers with co-located AlCl4- anions, a promising combination for energy storage. The 1980s bipolaron model remains the dominant interpretation of the electronic structure of PEDOT despite recent theoretical progress that has provided new definitions of bipolarons and polarons. By considering the influence of oligomer length, oxidation or anion concentration and spin state, we find no evidence for many of the assertions of the 1980s bipolaron model and so further contribute to a new understanding. No self-localisation of positive charges in PEDOT is found, as predicted by the bipolaron model at the hybrid functional level. Instead, our results show distortions that exhibit a single or a double peak in bond length alternations and charge density. Either can occur at different oxidation or anion concentrations. Rather than representing bipolarons or polaron pairs in the original model, these are electron distributions driven by a range of factors. Distortions can span an arbitrary number of nearby anions. We also contribute a novel conductivity hypothesis. Conductivity in conducting polymers has been observed to reduce at anion concentrations above 0.5. We show that at high anion concentrations, the energy of the localised, non-bonding anionic orbitals approaches that of the system HOMO due to Coulombic repulsion between anions. We hypothesize that with nucleic motion in the macropolymer, these orbitals will interfere with the hopping of charge carriers between sites of similar energy, lowering conductivity.

Machine learning for data-driven design of high-safety lithium metal anode.

Here, we present a protocol for developing an inorganic-organic hybrid interphase layer using the self-assembled monolayers technique to enhance the surface of the lithium metal anode. We describe steps for extracting organic molecules from open-sourced databases and calculating their microscopic properties. We then detail procedures for developing a machine learning model for predicting the ionic diffusion barrier and preparing the inputs for prediction. This protocol enables a cost-effective workflow to identify promising self-assembled monolayers with exceptional performance. For complete details on the use and execution of this protocol, please refer to Zhang et al. (2023).1.

Recent advances in developing multiscale descriptor approach for the design of oxygen redox electrocatalysts.

Oxygen redox electrocatalysis is the crucial electrode reaction among new-era energy sources. The prerequisite to rationally design an ideal electrocatalyst is accurately identifying the structure-activity relationship based on the so-called descriptors which link the catalytic performance with structural properties. However, the quick discovery of those descriptors remains challenging. In recent, the high-throughput computing and machine learning methods were identified to present great prospects for accelerating the screening of descriptors. That new research paradigm improves cognition in the way of oxygen evolution reaction/oxygen reduction reaction activity descriptor and reinforces the understanding of intrinsic physical and chemical features in the electrocatalytic process from a multiscale perspective. This review summarizes those new research paradigms for screening multiscale descriptors, especially from atomic scale to cluster mesoscale and bulk macroscale. The development of descriptors from traditional intermediate to eigen feature parameters has been addressed which provides guidance for the intelligent design of new energy materials.

Twisting Dynamics of Large Lattice-Mismatch van der Waals Heterostructures.

van der Waals (vdW) homo/heterostructures are ideal systems for studying interfacial tribological properties such as structural superlubricity. Previous studies concentrated on the mechanism of translational motion in vdW interfaces. However, detailed mechanisms and general properties of the rotational motion are barely explored. Here, we combine experiments and simulations to reveal the twisting dynamics of the MoS2/graphite heterostructure. Unlike the translational friction falling into the superlubricity regime with no twist angle dependence, the dynamic rotational resistances highly depend on twist angles. Our results show that the periodic rotational resistance force originates from structural potential energy changes during the twisting. The structural potential energy of MoS2/graphite heterostructure increases monotonically from 0° to 30° twist angles, and the estimated relative energy barrier is (1.43 ± 0.36) × 10-3 J/m2. The formation of Moiré superstructures in the graphene layer is the key to controlling the structural potential energy of the MoS2/graphene heterostructure. Our results suggest that in twisting 2D heterostructures, even if the interface sliding friction is negligible, the evolving potential energy change results in a nonvanishing rotational resistance force. The structural change of the heterostructure can be an additional pathway for energy dissipation in the rotational motion, further enhancing the rotational friction force.

Design Guidelines for Two-Dimensional Transition Metal Dichalcogenide Alloys.

Two-dimensional (2D) materials and transition metal dichalcogenides (TMD) in particular are at the forefront of nanotechnology. To tailor their properties for engineering applications, alloying strategies-used successfully for bulk metals in the last century-need to be extended to this novel class of materials. Here we present a systematic analysis of the phase behavior of substitutional 2D alloys in the TMD family on both the metal and the chalcogenide site. The phase behavior is quantified in terms of a metastability metric and benchmarked against systematic computational screening of configurational energy landscapes from First-Principles. The resulting Pettifor maps can be used to identify broad trends across chemical spaces and as starting point for setting up rational search strategies in phase space, thus allowing for targeted computational analysis of properties on likely thermodynamically stable compounds. The results presented here also constitute a useful guideline for synthesis of binary metal 2D TMDs alloys via a range of synthesis techniques.

UItra-low friction and edge-pinning effect in large-lattice-mismatch van der Waals heterostructures.

Two-dimensional heterostructures are excellent platforms to realize twist-angle-independent ultra-low friction due to their weak interlayer van der Waals interactions and natural lattice mismatch. However, for finite-size interfaces, the effect of domain edges on the friction process remains unclear. Here we report the superlubricity phenomenon and the edge-pinning effect at MoS2/graphite and MoS2/hexagonal boron nitride van der Waals heterostructure interfaces. We found that the friction coefficients of these heterostructures are below 10-6. Molecular dynamics simulations corroborate the experiments, which highlights the contribution of edges and interface steps to friction forces. Our experiments and simulations provide more information on the sliding mechanism of finite low-dimensional structures, which is vital to understand the friction process of laminar solid lubricants.

Electrochemistry from first-principles in the grand canonical ensemble.

Progress in electrochemical technologies, such as automotive batteries, supercapacitors, and fuel cells, depends greatly on developing improved charged interfaces between electrodes and electrolytes. The rational development of such interfaces can benefit from the atomistic understanding of the materials involved by first-principles quantum mechanical simulations with Density Functional Theory (DFT). However, such simulations are typically performed on the electrode surface in the absence of its electrolyte environment and at constant charge. We have developed a new hybrid computational method combining DFT and the Poisson-Boltzmann equation (P-BE) capable of simulating experimental electrochemistry under potential control in the presence of a solvent and an electrolyte. The charged electrode is represented quantum-mechanically via linear-scaling DFT, which can model nanoscale systems with thousands of atoms and is neutralized by a counter electrolyte charge via the solution of a modified P-BE. Our approach works with the total free energy of the combined multiscale system in a grand canonical ensemble of electrons subject to a constant electrochemical potential. It is calibrated with respect to the reduction potential of common reference electrodes, such as the standard hydrogen electrode and the Li metal electrode, which is used as a reference electrode in Li-ion batteries. Our new method can be used to predict electrochemical properties under constant potential, and we demonstrate this in exemplar simulations of the differential capacitance of few-layer graphene electrodes and the charging of a graphene electrode coupled to a Li metal electrode at different voltages.

Electronic structure calculations in electrolyte solutions: Methods for neutralization of extended charged interfaces.

Density functional theory (DFT) is often used for simulating extended materials such as infinite crystals or surfaces, under periodic boundary conditions (PBCs). In such calculations, when the simulation cell has non-zero charge, electrical neutrality has to be imposed, and this is often done via a uniform background charge of opposite sign ("jellium"). This artificial neutralization does not occur in reality, where a different mechanism is followed as in the example of a charged electrode in electrolyte solution, where the surrounding electrolyte screens the local charge at the interface. The neutralizing effect of the surrounding electrolyte can be incorporated within a hybrid quantum-continuum model based on a modified Poisson-Boltzmann equation, where the concentrations of electrolyte ions are modified to achieve electroneutrality. Among the infinite possible ways of modifying the electrolyte charge, we propose here a physically optimal solution, which minimizes the deviation of concentrations of electrolyte ions from those in open boundary conditions (OBCs). This principle of correspondence of PBCs with OBCs leads to the correct concentration profiles of electrolyte ions, and electroneutrality within the simulation cell and in the bulk electrolyte is maintained simultaneously, as observed in experiments. This approach, which we call the Neutralization by Electrolyte Concentration Shift (NECS), is implemented in our electrolyte model in the Order-N Electronic Total Energy Package (ONETEP) linear-scaling DFT code, which makes use of a bespoke highly parallel Poisson-Boltzmann solver, DL_MG. We further propose another neutralization scheme ("accessible jellium"), which is a simplification of NECS. We demonstrate and compare the different neutralization schemes on several examples.

Exploring the Stability of Twisted van der Waals Heterostructures.

Recent research showed that the rotational degree of freedom in stacking 2D materials yields great changes in the electronic properties. Here, we focus on an often overlooked question: are twisted geometries stable and what defines their rotational energy landscape? Our simulations show how epitaxy theory breaks down in these systems, and we explain the observed behavior in terms of an interplay between flexural phonons and the interlayer coupling, governed by the moiré superlattice. Our argument, applied to the well-studied MoS2/graphene system, rationalizes experimental results and could serve as guidance to design twistronic devices.

Controllable Heteroatom Doping Effects of CrxCo2-xP Nanoparticles: a Robust Electrocatalyst for Overall Water Splitting in Alkaline Solutions.

The effect of doping Cr on the electrocatalytic activity of Co2P supported on carbon black (CrxCo2-xP/CB) for the hydrogen evolution reaction (HER) and oxygen evolution reaction (OER) in alkaline solution was investigated. A beneficial improvement in the performance of Co2P toward HER and OER was discovered. For the HER at -200 mV overpotential, the turnover frequency (TOF) increases almost 6-fold from 0.26 to 1.52 electron siteCo-1 s-1 when Co2P/CB has a small amount of Cr added to form Cr0.2Co1.8P/CB. Similarly, we estimate an increase from 0.205 to 0.585 electron siteCo-1 s-1 for the OER at 1.6 V for the same change in composition. With 10 atom % Cr doping, the Cr0.2Co1.8P/CB catalyst needed 226 mV overpotential to produce a cathodic current density of -100 A gCo-1 and 380 mV overpotential to produce an anodic current density of 100 A gCo-1. Based on both experimental results and theoretical calculations, the activity improvement results from optimization of the electronic properties of Co2P after Cr doping.

A quick and versatile one step metal-organic chemical deposition method for supported Pt and Pt-alloy catalysts.

A simple, modified Metal-Organic Chemical Deposition (MOCD) method for Pt, PtRu and PtCo nanoparticle deposition onto a variety of support materials, including C, SiC, B4C, LaB6, TiB2, TiN and a ceramic/carbon nanofiber, is described. Pt deposition using Pt(acac)2 as a precursor is shown to occur via a mixed solid/liquid/vapour precursor phase which results in a high Pt yield of 90-92% on the support material. Pt and Pt alloy nanoparticles range 1.5-6.2 nm, and are well dispersed on all support materials, in a one-step method, with a total catalyst preparation time of ∼10 hours (2.4-4× quicker than conventional methods). The MOCD preparation method includes moderate temperatures of 350 °C in a tubular furnace with an inert gas supply at 2 bar, a high pressure (2-4 bar) compared to typical MOCVD methods (∼0.02-10 mbar). Pt/C catalysts with Pt loadings of 20, 40 and 60 wt% were synthesised, physically characterised, electrochemically characterised and compared to commercial Pt/C catalysts. TEM, XRD and ex situ EXAFS show similar Pt particle sizes and Pt particle shape identifiers, namely the ratio of the third to first Pt coordination numbers modelled from ex situ EXAFS, between the MOCD prepared catalysts and commercial catalysts. Moreover, electrochemical characterisation of the Pt/C MOCD catalysts obtained ORR mass activities with a maximum of 428 A gPt -1 at 0.9 V, which has similar mass activities to the commercial catalysts (80-160% compared to the commercial Pt/C catalysts).

Non-equilibrium crystallization pathways of manganese oxides in aqueous solution.

Aqueous precipitation of transition metal oxides often proceeds through non-equilibrium phases, whose appearance cannot be anticipated from traditional phase diagrams. Without a precise understanding of which metastable phases form, or their lifetimes, targeted synthesis of specific metal oxides can become a trial-and-error process. Here, we construct a theoretical framework to reveal the nanoscale and metastable energy landscapes of Pourbaix (E-pH) diagrams, providing quantitative insights into the size-dependent thermodynamics of metastable oxide nucleation and growth in water. By combining this framework with classical nucleation theory, we interrogate how solution conditions influence the multistage oxidation pathways of manganese oxides. We calculate that even within the same stability region of a Pourbaix diagram, subtle variations in pH and redox potential can redirect a non-equilibrium crystallization pathway through different metastable intermediates. Our theoretical framework offers a predictive platform to navigate through the thermodynamic and kinetic energy landscape towards the rational synthesis of target materials.

Electronic metal-support interaction enhanced oxygen reduction activity and stability of boron carbide supported platinum.

Catalysing the reduction of oxygen in acidic media is a standing challenge. Although activity of platinum, the most active metal, can be substantially improved by alloying, alloy stability remains a concern. Here we report that platinum nanoparticles supported on graphite-rich boron carbide show a 50-100% increase in activity in acidic media and improved cycle stability compared to commercial carbon supported platinum nanoparticles. Transmission electron microscopy and x-ray absorption fine structure analysis confirm similar platinum nanoparticle shapes, sizes, lattice parameters, and cluster packing on both supports, while x-ray photoelectron and absorption spectroscopy demonstrate a change in electronic structure. This shows that purely electronic metal-support interactions can significantly improve oxygen reduction activity without inducing shape, alloying or strain effects and without compromising stability. Optimizing the electronic interaction between the catalyst and support is, therefore, a promising approach for advanced electrocatalysts where optimizing the catalytic nanoparticles themselves is constrained by other concerns.

Mechano-electrochemistry effects due to deformation of copper oxide films.

In an attempt to elucidate the relationship and underlying processes of metal oxidation under stress, we combined the electrochemical characterisation with Density-Functional-Theory (DFT) calculations to interrogate the (100) surface of copper. The oxidised (100) surface shows a missing-row reconstruction, which is believed to be driven by surface stress. Hence, additional mechanical stimuli might have a significant impact on the onset of Cu oxidation. We find that different surface sites respond differently to strain. Oads at the thermodynamically favoured high-coordination hollow site (O coordinated to four Cu) is stabilised by up to 130 meV by imposing 2% tensile strain onto the surface, while the low-coordination top site (O coordinated to one Cu) shows a markedly different sensitivity. By cramping into the hollow site, Oads induces compressive stress into the (100) surface, an effect that is largely absent for the adsorption at the top site. We also find that the thermodynamic advantage of reconstructive underpotential surface oxidation is diminished under tensile strain. Hence, imposing tensile stress counter-balances the oxygen induced surface stress, which might have an implication on the onset of bulk copper oxidation. Studying Cu(100) single crystal surfaces in perchloric acid using cyclic voltammetry, we were able to confirm sensitivity of the electrochemical response towards the elastic strain.

Electrocatalytic performance of fuel cell reactions at low catalyst loading and high mass transport.

An alternative approach to the rotating disk electrode (RDE) for characterising fuel cell electrocatalysts is presented. The approach combines high mass transport with a flat, uniform, and homogeneous catalyst deposition process, well suited for studying intrinsic catalyst properties at realistic operating conditions of a polymer electrolyte fuel cell (PEFC). Uniform catalyst layers were produced with loadings as low as 0.16 µgPt cm(-2) and thicknesses as low as 200 nm. Such ultra thin catalyst layers are considered advantageous to minimize internal resistances and mass transport limitations. Geometric current densities as high as 5.7 A cm(-2)Geo were experimentally achieved at a loading of 10.15 µgPt cm(-2) for the hydrogen oxidation reaction (HOR) at room temperature, which is three orders of magnitude higher than current densities achievable with the RDE. Modelling of the associated diffusion field suggests that such high performance is enabled by fast lateral diffusion within the electrode. The electrodes operate over a wide potential range with insignificant mass transport losses, allowing the study of the ORR at high overpotentials. Electrodes produced a specific current density of 31 ± 9 mA cm(-2)Spec at a potential of 0.65 V vs. RHE for the oxygen reduction reaction (ORR) and 600 ± 60 mA cm(-2)Spec for the peak potential of the HOR. The mass activity of a commercial 60 wt% Pt/C catalyst towards the ORR was found to exceed a range of literature PEFC mass activities across the entire potential range. The HOR also revealed fine structure in the limiting current range and an asymptotic current decay for potentials above 0.36 V. These characteristics are not visible with techniques limited by mass transport in aqueous media such as the RDE.