Discharging of Ramsdellite MnO2 Cathode in a Lithium-Ion Battery.

We report an application of our unbiased Monte Carlo approach to investigate thermodynamic and electrochemical properties of lithiated manganese oxide in the ramsdellite phase (R-MnO2) to uncover the mechanism of lithium intercalation and understand charging/discharging of R-MnO2 as a cathode material in lithium-ion batteries. The lithium intercalation reaction was computationally explored by modeling thermodynamically significant distributions of lithium and reduced manganese in the R-MnO2 framework for a realistic range of lithium molar fractions 0 < x < 1 in Li x MnO2. We employed interatomic potentials and analyzed the thermodynamics of the resultant grand canonical ensemble. We found ordered or semiordered phases at x = 0.5 and 1.0 in Li x MnO2, verified by configurational entropy changes and simulated X-ray diffraction patterns of partially lithiated R-MnO2. The radial distribution functions show the preference of lithium for homogeneous distributions across the one-dimensional 2 × 1 ramsdellite channels accompanied by alternating reduced/oxidized manganese ions. The occupation of the interstitial sites in the channels is correlated with the calculated voltage profile, showing a sharp voltage drop at x = 0.5, which is explained by the energy penalty of shifting lithium ions from stable tetrahedral to unstable octahedral sites. To facilitate this work, our in-house software, Knowledge Led Master Code (KLMC) was extended to support massive parallelism on high-performance computers.

Environment-Driven Variability in Absolute Band Edge Positions and Work Functions of Reduced Ceria.

The absolute band edge positions and work function (Φ) are the key electronic properties of metal oxides that determine their performance in electronic devices and photocatalysis. However, experimental measurements of these properties often show notable variations, and the mechanisms underlying these discrepancies remain inadequately understood. In this work, we focus on ceria (CeO2), a material renowned for its outstanding oxygen storage capacity, and combine theoretical and experimental techniques to demonstrate environmental modifications of its ionization potential (IP) and Φ. Under O-deficient conditions, reduced ceria exhibits a decreased IP and Φ with significant sensitivity to defect distributions. In contrast, the IP and Φ are elevated in O-rich conditions due to the formation of surface peroxide species. Surface adsorbates and impurities can further augment these variabilities under realistic conditions. We rationalize the shifts in energy levels by separating the individual contributions from bulk and surface factors, using hybrid quantum mechanical/molecular mechanical (QM/MM) embedded-cluster and periodic density functional theory (DFT) calculations supported by interatomic-potential-based electrostatic analyses. Our results highlight the critical role of on-site electrostatic potentials in determining the absolute energy levels in metal oxides, implying a dynamic evolution of band edges under catalytic conditions.

Multiscale Investigation of the Mechanism and Selectivity of CO2 Hydrogenation over Rh(111).

CO2 hydrogenation over Rh catalysts comprises multiple reaction pathways, presenting a wide range of possible intermediates and end products, with selectivity toward either CO or methane being of particular interest. We investigate in detail the reaction mechanism of CO2 hydrogenation to the single-carbon (C1) products on the Rh(111) facet by performing periodic density functional theory (DFT) calculations and kinetic Monte Carlo (kMC) simulations, which account for the adsorbate interactions through a cluster expansion approach. We observe that Rh readily facilitates the dissociation of hydrogen, thus contributing to the subsequent hydrogenation processes. The reverse water-gas shift (RWGS) reaction occurs via three different reaction pathways, with CO hydrogenation to the COH intermediate being a key step for CO2 methanation. The effects of temperature, pressure, and the composition ratio of the gas reactant feed are considered. Temperature plays a pivotal role in determining the surface coverage and adsorbate composition, with competitive adsorption between CO and H species influencing the product distribution. The observed adlayer configurations indicate that the adsorbed CO species are separated by adsorbed H atoms, with a high ratio of H to CO coverage on the Rh(111) surface being essential to promote CO2 methanation.

Bulk and Surface Contributions to Ionisation Potentials of Metal Oxides.

Determining the absolute band edge positions in solid materials is crucial for optimising their performance in wide-ranging applications including photocatalysis and electronic devices. However, obtaining absolute energies is challenging, as seen in CeO2 , where experimental measurements show substantial discrepancies in the ionisation potential (IP). Here, we have combined several theoretical approaches, from classical electrostatics to quantum mechanics, to elucidate the bulk and surface contributions to the IP of metal oxides. We have determined a theoretical bulk contribution to the IP of stoichiometric CeO2 of only 5.38 eV, while surface orientation results in intrinsic IP variations ranging from 4.2 eV to 8.2 eV. Highly tuneable IPs were also found in TiO2 , ZrO2 , and HfO2 , in which surface polarisation plays a pivotal role in long-range energy level shifting. Our analysis, in addition to rationalising the observed range of experimental results, provides a firm basis for future interpretations of experimental and computational studies of oxide band structures.

Computational infrared and Raman spectra by hybrid QM/MM techniques: a study on molecular and catalytic material systems.

Vibrational spectroscopy is one of the most well-established and important techniques for characterizing chemical systems. To aid the interpretation of experimental infrared and Raman spectra, we report on recent theoretical developments in the ChemShell computational chemistry environment for modelling vibrational signatures. The hybrid quantum mechanical and molecular mechanical approach is employed, using density functional theory for the electronic structure calculations and classical forcefields for the environment. Computational vibrational intensities at chemical active sites are reported using electrostatic and fully polarizable embedding environments to achieve more realistic vibrational signatures for materials and molecular systems, including solvated molecules, proteins, zeolites and metal oxide surfaces, providing useful insight into the effect of the chemical environment on the signatures obtained from experiment. This work has been enabled by the efficient task-farming parallelism implemented in ChemShell for high-performance computing platforms.  This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Multiscale QM/MM modelling of catalytic systems with ChemShell.

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

Toward a Consistent Prediction of Defect Chemistry in CeO2.

Polarizable shell-model potentials are widely used for atomic-scale modeling of charged defects in solids using the Mott-Littleton approach and hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) embedded-cluster techniques. However, at the pure MM level of theory, the calculated defect energetics may not satisfy the requirement of quantitative predictions and are limited to only certain charged states. Here, we proposed a novel interatomic potential development scheme that unifies the predictions of all relevant charged defects in CeO2 based on the Mott-Littleton approach and QM/MM electronic-structure calculations. The predicted formation energies of oxygen vacancies accompanied by different excess electron localization patterns at the MM level of theory reach the accuracy of density functional theory (DFT) calculations using hybrid functionals. The new potential also accurately reproduces a wide range of physical properties of CeO2, showing excellent agreement with experimental and other computational studies. These findings provide opportunities for accurate large-scale modeling of the partial reduction and nonstoichiometry in CeO2, as well as a prototype for developing robust interatomic potentials for other defective crystals.

Parallel-Self-Assembling Stack, Center-Capture Effect, and Reactivity-Enhancing Effect of N-Layer (N = 1, 2, 3) Cyclo[18]carbon.

Cyclo[18]carbon (C18) is an captivating allotrope of carbon synthesized recently, which has drawn the attention among scientists. There are still few studies on the dynamic behaviors of C18. To gain knowledge in this area, we systematically explored the stacking behaviors and the oxidation kinetics of C18, as well the electronic transport behaviors of C18 oxides, by density functional theory and nonequilibrium Green's function calculations combined with reactive force field molecular dynamics simulations. The parallel-self-assembling behaviors were observed in the stack of two- or three-layer C18. During the oxidation process of C18, we found an evident center-capture effect in which the hollow rings would preferentially attract an O2 molecule into their centers. Moreover, the adsorption of O2 on the O2-doped rings was dramatically enhanced by the O2 at the center of the ring, showing the reactivity-enhancing effect. The excellent electron transport property of central-O2-doped C18 among 13 types of C18 oxides demonstrates the potential of C18 oxides as promising molecular devices for various applications. This study reveals the dynamic behaviors of C18 and provides theoretical guidance for use of C18 and C18 oxides in molecular devices.

Oxidation of graphene with variable defects: alternately symmetrical escape and self-restructuring of carbon rings.

Variable defects such as vacancies and grain boundaries are unavoidable in the synthesis of graphene, but play a central role in the activation of oxidation. Here, we apply reactive molecular dynamics simulations to reveal the underpinning mechanisms of oxidation in graphene with or without defects at the atomic scale. There exist four oxidation modes generating CO2 or CO in different stages, beginning from a single-atom vacancy, and proceeding until the ordered structure broken down into carbon oxide chains. The oxidation process of the graphene sheets experiences four typical stages, in which alternately symmetrical escape phenomenon is observed. Importantly, disordered rings can self-restructure during the oxidation of grain boundaries. Of all defects, the oxidation of vacancy has the lowest energy barrier and is therefore the easiest point of nucleation. This study demonstrates the crucial role of defects in determining the oxidation kinetics, and provides theoretical guidance for the oxidation prevention of graphene and the production of functionalized graphene.

Odd-Even Effects on Transport Properties of Polycyclic Arene Molecular Devices with Decreasing Numbers of Benzene Rings.

The electron transport properties of polycyclic aromatic hydrocarbons (PAHs) with different numbers of benzene rings tethered to narrow zigzag graphene nanoribbon (ZGNR) electrodes have been investigated. Results show that the transport properties of PAHs are dependent on whether the number of benzene rings in the width direction is odd or even. This effect is strong for narrow width PAHs, but its strength decreases as the width of the PAH is increased. PAHs with an odd number of rings exhibit poor transport properties, whereas the ones having an even number of rings show excellent transport properties coupled with a negative differential resistance (NDR) effect. Moreover, the linkage points and the structure of the molecules have a noticeable effect on the transport properties of the PAH, making the odd-even effect weaker or disappear entirely. Although the PAH with three benzene rings displays poor transport capabilities, it shows excellent rectification behavior compared to the other examined molecules. These studies present a feasible avenue for designing molecular devices with enhanced performance by the careful manipulation of the PAH molecular structure.

Crystallization behavior of a confined CuZr metallic liquid film with a sandwich-like structure.

Despite the fact that its crystal state is thermodynamically stable, Cu64Zr36 alloy is prone to form metastable glass at a high cooling rate. However, the confinement can induce nano-crystallization with a novel sandwich-like hierarchical structure consisting of pure Cu layers, pure Zr layers and mixed layers by conducting molecular dynamics simulations. The liquid-to-crystal transition temperature and interatomic repulsion softness display abnormal oscillations, instead of monotonous variation, as the wall-wall separation increases. When the confinement size is 10 Å and 12 Å, the transition temperature reaches a maximum, resulting from the pending new sandwich layer. The atomic movement and dynamical heterogeneity are demonstrated to play a vital role in the abnormal oscillation behavior of physical properties of the nano confined metallic glass. The sandwich-like structure can alter the Cu-Zr bond fraction, which eventually influences the liquid-to-crystal transition temperature and interatomic repulsion softness. Our findings provide a deep insight into the hierarchical nanostructures and its liquid-to-crystal transition characteristics under confinement at the atomic level.

Atomistic Origin of the Complex Morphological Evolution of Aluminum Nanoparticles during Oxidation: A Chain-like Oxide Nucleation and Growth Mechanism.

Metal nanoparticles usually show different oxidation dynamics from bulk metals, which results in various oxide nanostructures because of their size-related surface effects. In this work, we have found and investigated the chain-like nucleation and growth of oxides on the aluminum nanoparticle (ANP) surface, using molecular dynamics simulations with the reactive force-field (ReaxFF). After nucleation, the chain-like oxide nuclei could stay on the ANP surface and continue growing into an oxide shell, extend outward from the surface to form longer oxide chains, or detach from the ANP to generate independent oxide clusters, which is highly dependent on the oxygen content, temperature, and nanoparticle size. Our results emphasize the complicated interplay between the surface structure of nanoparticles and the environmental conditions in determining the formation of oxides, which provides insights into the atomic-scale oxidation mechanism of metal nanoparticles.

Density dependent structural phase transition for confined copper: origin of the layering.

Confinement presents the opportunity for novel structural transition scenarios not observed in three-dimensional systems. Here, we report a comprehensive molecular dynamic (MD) study of the structural phase transition induced by density for an ordinary metal copper (Cu) confined between two parallel panel walls. At 4.19 g cm-3 < ρ < 4.66 g cm-3, a notable structural phase transition occurs between the triangle unit cell structure and quasi-square unit cell structure upon densification. Both the bond order parameter (BOP) and angular distribution function (ADF) can provide evidence for the transition. We highlight the fact that when the sole decrease of the atom distance cannot adapt to the further densification, the system starts to adjust the neighboring bond angle and promote the layering transition, thus inducing the structural phase transition. At the metastable coexistence zone, the viscosity exhibits a remarkable drop and the diffusion coefficient shows a notable increase, both facilitating the accomplishment of the structural transition. Our results will trigger more interest on the phase transition under confinement in a metallic system.